Using quantitative wastewater analysis to measure daily usage of conventional and emerging illicit drugs at an annual music festival

Introduction and Aims Wastewater analysis provides a non-intrusive way of measuring drug use within a population. We used this approach to determine daily use of conventional illicit drugs [cannabis, cocaine, methamphetamine and 3,4-methylenedioxymethamphetamine (MDMA)] and emerging illicit psychostimulants (benzylpiperazine, mephedrone and methylone) in two consecutive years (2010 and 2011) at an annual music festival. Design and Methods Daily composite wastewater samples, representative of the festival, were collected from the on-site wastewater treatment plant and analysed for drug metabolites. Data over 2 years were compared using Wilcoxon matched-pair test. Data from 2010 festival were compared with data collected at the same time from a nearby urban community using equivalent methods. Results Conventional illicit drugs were detected in all samples whereas emerging illicit psychostimulants were found only on specific days. The estimated per capita consumption of MDMA, cocaine and cannabis was similar between the two festival years. Statistically significant (P < 0.05; Z = −2.0–2.2) decreases were observed in use of methamphetamine and one emerging illicit psychostimulant (benzyl piperazine). Only consumption of MDMA was elevated at the festival compared with the nearby urban community. Discussion and Conclusions Rates of substance use at this festival remained relatively consistent over two monitoring years. Compared with the urban community, drug use among festival goers was only elevated for MDMA, confirming its popularity in music settings. Our study demonstrated that wastewater analysis can objectively capture changes in substance use at a music setting without raising major ethical issues. It would potentially allow effective assessments of drug prevention strategies in such settings in the future.

Challenges and opportunities in using wastewater analysis to measure drug use in a small prison facility

Introduction and Aims Wastewater analysis (WWA) is intended to be a direct and objective method of measuring substance use in large urban populations. It has also been used to measure prison substance use in two previous studies. The application of WWA in this context has raised questions as to how best it might be used to measure illicit drug use in prisons, and whether it can also be used to measure prescription misuse. We applied WWA to a small regional prison to measure the use of 12 licit and illicit substances. We attempted to measure the non-medical use of methadone and buprenorphine and to compare our findings with the results of the prison's mandatory drug testing (MDT). Design and Methods Representative daily composite samples were collected for two periods of 12 consecutive days in May to July 2013 and analysed for 18 drug metabolites. Prescription data and MDT results were obtained from the prison and compared with the substance use estimates calculated from WWA data. Results Daily use of methamphetamine, methadone, buprenorphine and codeine was detected, while sporadic detection of ketamine and methylone was also observed. Overall buprenorphine misuse appeared to be greater than methadone misuse. Discussion and Conclusions Compared with MDT, WWA provides a more comprehensive picture of prison substance use. WWA also has the potential to measure the misuse of medically prescribed substances. However, a great deal of care must be exercised in quantifying the usage of any substance in small populations, such as in prisons.

Systematic and Day-to-Day Effects of Chemical-Derived Population Estimates on Wastewater-Based Drug Epidemiology

Population size is crucial when estimating population-normalized drug consumption (PNDC) from wastewater-based drug epidemiology (WBDE). Three conceptually different population estimates can be used: de jure (common census, residence), de facto (all persons within a sewer catchment), and chemical loads (contributors to the sampled wastewater). De facto and chemical loads will be the same where all households contribute to a central sewer system without wastewater loss. This study explored the feasibility of determining a de facto population and its effect on estimating PNDC in an urban community over an extended period. Drugs and other chemicals were analyzed in 311 daily composite wastewater samples. The daily estimated de facto population (using chemical loads) was on average 32% higher than the de jure population. Consequently, using the latter would systemically overestimate PNDC by 22%. However, the relative day-to-day pattern of drug consumption was similar regardless of the type of normalization as daily illicit drug loads appeared to vary substantially more than the population. Using chemical loads population, we objectively quantified the total methodological uncertainty of PNDC and reduced it by a factor of 2. Our study illustrated the potential benefits of using chemical loads population for obtaining more robust PNDC data in WBDE.

A Model to Estimate the Population Contributing to the Wastewater Using Samples Collected on Census Day

An important uncertainty when estimating per capita consumption of, for example, illicit drugs by means of wastewater analysis (sometimes referred to as “sewage epidemiology”) relates to the size and variability of the de facto population in the catchment of interest. In the absence of a day-specific direct population count any indirect surrogate model to estimate population size lacks a standard to assess associated uncertainties. Therefore, the objective of this study was to collect wastewater samples at a unique opportunity, that is, on a census day, as a basis for a model to estimate the number of people contributing to a given wastewater sample. Mass loads for a wide range of pharmaceuticals and personal care products were quantified in influents of ten sewage treatment plants (STP) serving populations ranging from approximately 3500 to 500 000 people. Separate linear models for population size were estimated with the mass loads of the different chemical as the explanatory variable: 14 chemicals showed good, linear relationships, with highest correlations for acesulfame and gabapentin. De facto population was then estimated through Bayesian inference, by updating the population size provided by STP staff (prior knowledge) with measured chemical mass loads. Cross validation showed that large populations can be estimated fairly accurately with a few chemical mass loads quantified from 24-h composite samples. In contrast, the prior knowledge for small population sizes cannot be improved substantially despite the information of multiple chemical mass loads. In the future, observations other than chemical mass loads may improve this deficit, since Bayesian inference allows including any kind of information relating to population size.

Estimating daily and diurnal variations of illicit drug use in Hong Kong: A pilot study of using wastewater analysis in an Asian metropolitan city

The measurement of illicit drug metabolites in raw wastewater is increasingly being adopted as an approach to objectively monitor population-level drug use, and is an effective complement to traditional epidemiological methods. As such, it has been widely applied in western countries. In this study, we utilised this approach to assess drug use patterns over nine days during April 2011 in Hong Kong. Raw wastewater samples were collected from the largest wastewater treatment plant serving a community of approximately 3.5 million people and analysed for excreted drug residues including cocaine, ketamine, methamphetamine, 3,4-methylenedioxymethamphetamine (MDMA) and key metabolites using liquid chromatography coupled with tandem mass spectrometry. The overall drug use pattern determined by wastewater analysis was consistent with that have seen amongst people coming into contact with services in relation to substance use; among our target drugs, ketamine (estimated consumption: 1400–1600 mg/day/1000 people) was the predominant drug followed by methamphetamine (180–200 mg/day/1000 people), cocaine (160–180 mg/day/1000 people) and MDMA (not detected). The levels of these drugs were relatively steady throughout the monitoring period. Analysing samples at higher temporal resolution provided data on diurnal variations of drug residue loads. Elevated ratios of cocaine to benzoylecgonine were identified unexpectedly in three samples during the evening and night, providing evidence for potential dumping events of cocaine. This study provides the first application of wastewater analysis to quantitatively evaluate daily drug use in an Asian metropolitan community. Our data reinforces the benefit of wastewater monitoring to health and law enforcement authorities for strategic planning and evaluation of drug intervention strategies.

Could wastewater analysis be a useful tool for China? — A review

Analysing wastewater samples is an innovative approach that overcomes many limitations of traditional surveys to identify and measure a range of chemicals that were consumed by or exposed to people living in a sewer catchment area. First conceptualised in 2001, much progress has been made to make wastewater analysis (WWA) a reliable and robust tool for measuring chemical consumption and/or exposure. At the moment, the most popular application of WWA, sometimes referred as sewage epidemiology, is to monitor the consumption of illicit drugs in communities around the globe, including China. The approach has been largely adopted by law enforcement agencies as a device to monitor the temporal and geographical patterns of drug consumption. In the future, the methodology can be extended to other chemicals including biomarkers of population health (e.g. environmental or oxidative stress biomarkers, lifestyle indicators or medications that are taken by different demographic groups) and pollutants that people are exposed to (e.g. polycyclic aromatic hydrocarbons, perfluorinated chemicals, and toxic pesticides). The extension of WWA to a huge range of chemicals may give rise to a field called sewage chemical-information mining (SCIM) with unexplored potentials. China has many densely populated cities with thousands of sewage treatment plants which are favourable for applying WWA/SCIM in order to help relevant authorities gather information about illicit drug consumption and population health status. However, there are some prerequisites and uncertainties of the methodology that should be addressed for SCIM to reach its full potential in China.

Monitoring substance use in prisons: Assessing the potential value of wastewater analysis

Prison substance use is a major concern for prison authorities and the wider community. Australia has responded to this problem by implementing the National Corrections Drug Strategy. Across Australia, the true extent of prison substance use cannot be determined. As a result, the effectiveness of the interventions employed as part of this strategy cannot be properly assessed. This has important implications for the allocation of corrective services resources and future policy development. This article explores the benefits and limitations, as well as the ethical and practical issues in using wastewater analysis (WWA) to measure levels of substance use in prisons. It reports results from the first application of WWA to an Australian prison, which supports the use of WWA in this context. Given the increasing concern for prescription misuse in prisons, we also highlight the novel use of WWA to measure the extent of prescription misuse by prisoners. The article concludes that as a result of its objectivity, sensitivity and cost-effectiveness, the use of WWA in prisons warrants further consideration in Australia.

The first application of wastewater-based drug epidemiology in five South Korean cities

Illicit drug consumption in five cities in South Korea was estimated by analyzing 17 drug residues in untreated wastewater samples collected during the Christmas and New Year period of 2012-13. Only methamphetamine, amphetamine, and codeine were detected at concentrations of tens of nanograms per liter or even lower concentrations in more than 90% of the samples. Other illicit drug residues (including cocaine, methadone, and benzoylecgonine) that have been detected frequently in wastewater from other countries were not found in this study. Methamphetamine was found to be the most widely used illicit drug in South Korea, and the estimated average consumption rate was 22 mg d−1 (1000 people)−1. This rate is, for example, 2-5 times lower than the estimated average consumption rates in Hong Kong and other parts of China and 4-80 times lower than the estimated average consumption rates in cities in Western countries. It should be noted that the wastewater samples analyzed in this study were collected during a holiday season, when daily consumption of illicit drugs is often higher than on an average day. The methamphetamine usage rates were calculated for different cities in South Korea, and the usage rates in smaller cities was higher (2-4 times) than the average.

Measuring drug use patterns in Queensland through wastewater analysis

Estimating the use of illicit drugs in the general community is an important task with ramifications for law enforcement agencies, as well as health portfolios. Australia has four ongoing drug monitoring systems, including the AIC’s DUMA program, the National Drug Strategy Household Survey, the Illicit Drug Reporting System and the Ecstasy and Related Drug Reporting System. The systems vary in methods, but broadly they are reliant upon self-report data and may be subject to selection biases. The present study employed a completely different method. By chemically analysing sewerage water, the study produced daily estimates of consumption of methamphetamine, MDMA and cocaine. Samples were collected in November 2009 and November 2010 from a municipality in Queensland, with an population of over 150,000 people. Estimates were made of the average daily dose and average daily street value per 1,000 people. On the basis of estimated dose and price, the methamphetamine market appeared considerably stronger than either MDMA or cocaine. This paper explains the strengths and weaknesses of wastewater analysis. It considers the potential value of wastewater analysis in measuring net consumption of illicit drugs and the effectiveness of law enforcement agency strategies.

Effects of sewer conditions on the degradation of selected illicit drug residues in wastewater

The stability of five illicit drug markers in wastewater was tested under different sewer conditions using laboratory-scale sewer reactors. Wastewater was spiked with deuterium labelled isotopes of cocaine, benzoyl ecgonine, methamphetamine, MDMA and 6-acetyl morphine to avoid interference from the native isotopes already present in the wastewater matrix. The sewer reactors were operated at 20 °C and pH 7.5, and wastewater was sampled at 0, 0.25, 0.5, 1, 2, 3, 6, 9 and 12 h to measure the transformation/degradation of these marker compounds. The results showed that while methamphetamine, MDMA and benzoyl ecgonine were stable in the sewer reactors, cocaine and 6-acetyl morphine degraded quickly. Their degradation rates are significantly higher than the values reportedly measured in wastewater alone (without biofilms). All the degradation processes followed first order kinetics. Benzoyl ecgonine and morphine were also formed from the degradation of cocaine and 6-acetyl morphine, respectively, with stable formation rates throughout the test. These findings suggest that, in sewage epidemiology, it is essential to have relevant information of the sewer system (i.e. type of sewer, hydraulic retention time) in order to accurately back-estimate the consumption of illicit drugs. More research is required to look into detailed sewer conditions (e.g. temperature, pH and ratio of biofilm area to wastewater volume among others) to identify their effects on the fate of illicit drug markers in sewer systems.

Global, regional, and national disability-adjusted life years (DALYs) for 306 diseases and injuries and healthy life expectancy (HALE) for 188 countries, 1990–2013: Quantifying the epidemiological transition

Background The Global Burden of Disease Study 2013 (GBD 2013) aims to bring together all available epidemiological data using a coherent measurement framework, standardised estimation methods, and transparent data sources to enable comparisons of health loss over time and across causes, age–sex groups, and countries. The GBD can be used to generate summary measures such as disability-adjusted life-years (DALYs) and healthy life expectancy (HALE) that make possible comparative assessments of broad epidemiological patterns across countries and time. These summary measures can also be used to quantify the component of variation in epidemiology that is related to sociodemographic development. Methods We used the published GBD 2013 data for age-specific mortality, years of life lost due to premature mortality (YLLs), and years lived with disability (YLDs) to calculate DALYs and HALE for 1990, 1995, 2000, 2005, 2010, and 2013 for 188 countries. We calculated HALE using the Sullivan method; 95% uncertainty intervals (UIs) represent uncertainty in age-specific death rates and YLDs per person for each country, age, sex, and year. We estimated DALYs for 306 causes for each country as the sum of YLLs and YLDs; 95% UIs represent uncertainty in YLL and YLD rates. We quantified patterns of the epidemiological transition with a composite indicator of sociodemographic status, which we constructed from income per person, average years of schooling after age 15 years, and the total fertility rate and mean age of the population. We applied hierarchical regression to DALY rates by cause across countries to decompose variance related to the sociodemographic status variable, country, and time. Findings Worldwide, from 1990 to 2013, life expectancy at birth rose by 6·2 years (95% UI 5·6–6·6), from 65·3 years (65·0–65·6) in 1990 to 71·5 years (71·0–71·9) in 2013, HALE at birth rose by 5·4 years (4·9–5·8), from 56·9 years (54·5–59·1) to 62·3 years (59·7–64·8), total DALYs fell by 3·6% (0·3–7·4), and age-standardised DALY rates per 100 000 people fell by 26·7% (24·6–29·1). For communicable, maternal, neonatal, and nutritional disorders, global DALY numbers, crude rates, and age-standardised rates have all declined between 1990 and 2013, whereas for non–communicable diseases, global DALYs have been increasing, DALY rates have remained nearly constant, and age-standardised DALY rates declined during the same period. From 2005 to 2013, the number of DALYs increased for most specific non-communicable diseases, including cardiovascular diseases and neoplasms, in addition to dengue, food-borne trematodes, and leishmaniasis; DALYs decreased for nearly all other causes. By 2013, the five leading causes of DALYs were ischaemic heart disease, lower respiratory infections, cerebrovascular disease, low back and neck pain, and road injuries. Sociodemographic status explained more than 50% of the variance between countries and over time for diarrhoea, lower respiratory infections, and other common infectious diseases; maternal disorders; neonatal disorders; nutritional deficiencies; other communicable, maternal, neonatal, and nutritional diseases; musculoskeletal disorders; and other non-communicable diseases. However, sociodemographic status explained less than 10% of the variance in DALY rates for cardiovascular diseases; chronic respiratory diseases; cirrhosis; diabetes, urogenital, blood, and endocrine diseases; unintentional injuries; and self-harm and interpersonal violence. Predictably, increased sociodemographic status was associated with a shift in burden from YLLs to YLDs, driven by declines in YLLs and increases in YLDs from musculoskeletal disorders, neurological disorders, and mental and substance use disorders. In most country-specific estimates, the increase in life expectancy was greater than that in HALE. Leading causes of DALYs are highly variable across countries. Interpretation Global health is improving. Population growth and ageing have driven up numbers of DALYs, but crude rates have remained relatively constant, showing that progress in health does not mean fewer demands on health systems. The notion of an epidemiological transition—in which increasing sociodemographic status brings structured change in disease burden—is useful, but there is tremendous variation in burden of disease that is not associated with sociodemographic status. This further underscores the need for country-specific assessments of DALYs and HALE to appropriately inform health policy decisions and attendant actions.

Global, regional, and national comparative risk assessment of 79 behavioural, environmental and occupational, and metabolic risks or clusters of risks in 188 countries, 1990–2013: a systematic analysis for the Global Burden of Disease Study 2013

Background The Global Burden of Disease, Injuries, and Risk Factor study 2013 (GBD 2013) is the first of a series of annual updates of the GBD. Risk factor quantification, particularly of modifiable risk factors, can help to identify emerging threats to population health and opportunities for prevention. The GBD 2013 provides a timely opportunity to update the comparative risk assessment with new data for exposure, relative risks, and evidence on the appropriate counterfactual risk distribution. Methods Attributable deaths, years of life lost, years lived with disability, and disability-adjusted life-years (DALYs) have been estimated for 79 risks or clusters of risks using the GBD 2010 methods. Risk–outcome pairs meeting explicit evidence criteria were assessed for 188 countries for the period 1990–2013 by age and sex using three inputs: risk exposure, relative risks, and the theoretical minimum risk exposure level (TMREL). Risks are organised into a hierarchy with blocks of behavioural, environmental and occupational, and metabolic risks at the first level of the hierarchy. The next level in the hierarchy includes nine clusters of related risks and two individual risks, with more detail provided at levels 3 and 4 of the hierarchy. Compared with GBD 2010, six new risk factors have been added: handwashing practices, occupational exposure to trichloroethylene, childhood wasting, childhood stunting, unsafe sex, and low glomerular filtration rate. For most risks, data for exposure were synthesised with a Bayesian meta-regression method, DisMod-MR 2.0, or spatial-temporal Gaussian process regression. Relative risks were based on meta-regressions of published cohort and intervention studies. Attributable burden for clusters of risks and all risks combined took into account evidence on the mediation of some risks such as high body-mass index (BMI) through other risks such as high systolic blood pressure and high cholesterol. Findings All risks combined account for 57·2% (95% uncertainty interval [UI] 55·8–58·5) of deaths and 41·6% (40·1–43·0) of DALYs. Risks quantified account for 87·9% (86·5–89·3) of cardiovascular disease DALYs, ranging to a low of 0% for neonatal disorders and neglected tropical diseases and malaria. In terms of global DALYs in 2013, six risks or clusters of risks each caused more than 5% of DALYs: dietary risks accounting for 11·3 million deaths and 241·4 million DALYs, high systolic blood pressure for 10·4 million deaths and 208·1 million DALYs, child and maternal malnutrition for 1·7 million deaths and 176·9 million DALYs, tobacco smoke for 6·1 million deaths and 143·5 million DALYs, air pollution for 5·5 million deaths and 141·5 million DALYs, and high BMI for 4·4 million deaths and 134·0 million DALYs. Risk factor patterns vary across regions and countries and with time. In sub-Saharan Africa, the leading risk factors are child and maternal malnutrition, unsafe sex, and unsafe water, sanitation, and handwashing. In women, in nearly all countries in the Americas, north Africa, and the Middle East, and in many other high-income countries, high BMI is the leading risk factor, with high systolic blood pressure as the leading risk in most of Central and Eastern Europe and south and east Asia. For men, high systolic blood pressure or tobacco use are the leading risks in nearly all high-income countries, in north Africa and the Middle East, Europe, and Asia. For men and women, unsafe sex is the leading risk in a corridor from Kenya to South Africa. Interpretation Behavioural, environmental and occupational, and metabolic risks can explain half of global mortality and more than one-third of global DALYs providing many opportunities for prevention. Of the larger risks, the attributable burden of high BMI has increased in the past 23 years. In view of the prominence of behavioural risk factors, behavioural and social science research on interventions for these risks should be strengthened. Many prevention and primary care policy options are available now to act on key risks.

Identification of a novel fusion transcript between human relaxin-1 (RLN1) and human relaxin-2 (RLN2) in prostate cancer

Simultaneous expression of highly homologous RLN1 and RLN2 genes in prostate impairs their accurate delineation. We used PacBio SMRT sequencing and RNA-Seq in LNCaP cells in order to dissect the expression of RLN1 and RLN2 variants. We identified a novel fusion transcript comprising the RLN1 and RLN2 genes and found evidence of its expression in the normal and prostate cancer tissues. The RLN1-RLN2 fusion putatively encodes RLN2 isoform with the deleted secretory signal peptide. The identification of the fusion transcript provided information to determine unique RLN1-RLN2 fusion and RLN1 regions. The RLN1-RLN2 fusion was co-expressed with RLN1 in LNCaP cells, but the two gene products were inversely regulated by androgens. We showed that RLN1 is underrepresented in common PCa cell lines in comparison to normal and PCa tissue. The current study brings a highly relevant update to the relaxin field, and will encourage further studies of RLN1 and RLN2 in PCa and broader.

Analysis of differential gene expression based on Bayesian estimation of variance

Gene expression is arguably the most important indicator of biological function. Thus identifying differentially expressed genes is one of the main aims of high throughout studies that use microarray and RNAseq platforms to study deregulated cellular pathways. There are many tools for analysing differentia gene expression from transciptomic datasets. The major challenge of this topic is to estimate gene expression variance due to the high amount of ‘background noise’ that is generated from biological equipment and the lack of biological replicates. Bayesian inference has been widely used in the bioinformatics field. In this work, we reveal that the prior knowledge employed in the Bayesian framework also helps to improve the accuracy of differential gene expression analysis when using a small number of replicates. We have developed a differential analysis tool that uses Bayesian estimation of the variance of gene expression for use with small numbers of biological replicates. Our method is more consistent when compared to the widely used cyber-t tool that successfully introduced the Bayesian framework to differential analysis. We also provide a user-friendly web based Graphic User Interface for biologists to use with microarray and RNAseq data. Bayesian inference can compensate for the instability of variance caused when using a small number of biological replicates by using pseudo replicates as prior knowledge. We also show that our new strategy to select pseudo replicates will improve the performance of the analysis. - See more at: http://www.eurekaselect.com/node/138761/article#sthash.VeK9xl5k.dpuf

Fusion transcript loci share many genomic features with non-fusion loci

Background Fusion transcripts are found in many tissues and have the potential to create novel functional products. Here, we investigate the genomic sequences around fusion junctions to better understand the transcriptional mechanisms mediating fusion transcription/splicing. We analyzed data from prostate (cancer) cells as previous studies have shown extensively that these cells readily undergo fusion transcription. Results We used the FusionMap program to identify high-confidence fusion transcripts from RNAseq data. The RNAseq datasets were from our (N = 8) and other (N = 14) clinical prostate tumors with adjacent non-cancer cells, and from the LNCaP prostate cancer cell line that were mock-, androgen- (DHT), and anti-androgen- (bicalutamide, enzalutamide) treated. In total, 185 fusion transcripts were identified from all RNAseq datasets. The majority (76 %) of these fusion transcripts were ‘read-through chimeras’ derived from adjacent genes in the genome. Characterization of sequences at fusion loci were carried out using a combination of the FusionMap program, custom Perl scripts, and the RNAfold program. Our computational analysis indicated that most fusion junctions (76 %) use the consensus GT-AG intron donor-acceptor splice site, and most fusion transcripts (85 %) maintained the open reading frame. We assessed whether parental genes of fusion transcripts have the potential to form complementary base pairing between parental genes which might bring them into physical proximity. Our computational analysis of sequences flanking fusion junctions at parental loci indicate that these loci have a similar propensity as non-fusion loci to hybridize. The abundance of repetitive sequences at fusion and non-fusion loci was also investigated given that SINE repeats are involved in aberrant gene transcription. We found few instances of repetitive sequences at both fusion and non-fusion junctions. Finally, RT-qPCR was performed on RNA from both clinical prostate tumors and adjacent non-cancer cells (N = 7), and LNCaP cells treated as above to validate the expression of seven fusion transcripts and their respective parental genes. We reveal that fusion transcript expression is similar to the expression of parental genes. Conclusions Fusion transcripts maintain the open reading frame, and likely use the same transcriptional machinery as non-fusion transcripts as they share many genomic features at splice/fusion junctions.