Impact of converging sociocultural and substance-related trends on US autism rates: combined geospatiotemporal and causal inferential analysis.

Whilst cannabis is known to be toxic to brain development, it is unknown if it is driving rising US autism rates (ASMR). A longitudinal epidemiological study was conducted using national autism census data from the US Department of Education Individuals with Disabilities Act (IDEA) 1991-2011 and nationally representative drug exposure (cigarettes, alcohol, analgesic, and cocaine abuse, and cannabis use monthly, daily, and in pregnancy) datasets from National Survey of Drug Use and Health and US Census (income and ethnicity) and CDC Wonder population and birth data. Analysis was conducted in R. 266,950 were autistic of a population of 40,119,464 8-year-olds in 1994-2011. At national level after adjustment, daily cannabis use was significantly related to ASMR (ß estimate = 4.37 (95%C.I. 4.06, 4.68), P < 2.2 × 10-16) as was first pregnancy trimester cannabis exposure (ß estimate = 0.12 (0.08, 0.16), P = 1.7 × 10-12). At state level following adjustment for cannabis, cannabigerol (from ß estimate = - 13.77 (- 19.41, 8.13), P = 1.8 × 10-6) and Δ9-tetrahydrocannabinol (from ß estimate = 1.96 (0.88-3.04), P = 4 × 10-4) were significant. Geospatial state-level modelling showed exponential relationship between ASMR and Δ9-tetrahydrocannabinol and cannabigerol exposure. Exponential coefficients for the relationship between modelled ASMR and Δ9-tetrahydrocannabinol and cannabigerol exposure were 7.053 (6.39-7.71) and 185.334 (167.88-202.79; both P < 2.0 × 10-7). E-values are an instrument related to the evidence for causality in observational studies. High E-values were noted. Dichotomized legal status was linked with elevated ASMR. Data show cannabis use is associated with ASMR, is powerful enough to affect overall trends, and persists after controlling for other major covariates. Cannabinoids are exponentially associated with ASMR. The cannabis-autism relationship satisfies criteria of causal inference.

Geotemporospatial and causal inference epidemiological analysis of US survey and overview of cannabis, cannabidiol and cannabinoid genotoxicity in relation to congenital anomalies 2001-2015.

BACKGROUND: Cannabinoids including cannabidiol have recognized genotoxic activities but their significance has not been studied broadly epidemiologically across the teratological spectrum. We examined these issues including contextual space-time relationships and formal causal inferential analysis in USA. METHODS: State congenital anomaly (CA) rate (CAR) data was taken from the annual reports of the National Birth Defects Prevention Network 2001-2005 to 2011-2015. Substance abuse rates were from the National Survey of Drug Use and Health a nationally representative longitudinal survey of the non-institutionalized US population with 74.1% response rate. Drugs examined were cigarettes, monthly and binge alcohol, monthly cannabis and analgesic and cocaine abuse. Early termination of pregnancy for abortion (ETOPFA) rates were taken from the published literature. Cannabinoid concentrations were from Drug Enforcement Agency. Ethnicity and income data were from the US Census Bureau. Inverse probability weighted (IPW) regressions and geotemporospatial regressions conducted for selected CAs. RESULTS: Data on 18,328,529 births from an aggregated population of 2,377,483,589 for mid-year analyses 2005-2013 comprehending 12,611 CARs for 62 CAs was assembled and ETOPFA-corrected (ETOPFACAR) where appropriate. E-Values for ETOPFACARs by substance trends were elevated for THC (40 CAs), cannabis (35 CAs), tobacco (11 CAs), cannabidiol (8 CAs), monthly alcohol (5 CAs) and binge alcohol (2 CAs) with minimum E-Values descending from 16.55, 1.55x107, 555.10, 7.53x1019, 9.30 and 32.98. Cardiovascular, gastrointestinal, chromosomal, limb reductions, urinary, face and body wall CAs particularly affected. Highest v. lowest substance use quintile CAR prevalence ratios 2.84 (95%C.I. 2.44, 3.31), 4.85 (4.08, 5.77) and 1.92 (1.63, 2.27) and attributable fraction in exposed 0.28 (0.27, 0.28), 0.57 (0.51, 0.62) and 0.47 (0.38, 0.55) for tobacco, cannabis and cannabidiol. Small intestinal stenosis or atresia and obstructive genitourinary defect were studied in detail in lagged IPW pseudo-randomized causal regressions and spatiotemporal models confirmed the causal role of cannabinoids. Spatiotemporal predictive modelling demonstrated strongly sigmoidal non-linear cannabidiol dose-response power-function relationships (P = 2.83x10-60 and 1.61x10-71 respectively). CONCLUSIONS: Data implicate cannabinoids including cannabidiol in a diverse spectrum of heritable CAs. Sigmoidal non-linear dose-response relationships are of grave concern. These transgenerational genotoxic, epigenotoxic, chromosomal-toxic putatively causal teratogenic effects strongly indicate tight restrictions on community cannabinoid penetration.

Cannabinoid exposure as a major driver of pediatric acute lymphoid Leukaemia rates across the USA: combined geospatial, multiple imputation and causal inference study.

BACKGROUND: Acute lymphoid leukaemia (ALL) is the commonest childhood cancer whose incidence is rising in many nations. In the USA, between 1975 and 2016, ALL rates (ALLRs) rose 93.51% from 1.91 to 3.70/100,000 <  20 years. ALL is more common in Caucasian-Americans than amongst minorities. The cause of both the rise and the ethnic differential is unclear, however, prenatal cannabis exposure was previously linked with elevated childhood leukaemia rates. We investigated epidemiologically if cannabis use impacted nationally on ALLRs, its ethnic effects, and if the relationship was causal. METHODS: State data on overall, and ethnic ALLR from the Surveillance Epidemiology and End Results databank of the Centre for Disease Control (CDC) and National Cancer Institute (NCI) were combined with drug (cigarettes, alcoholism, cannabis, analgesics, cocaine) use data from the National Survey of Drug Use and Health; 74.1% response rate. Income and ethnicity data was from the US Census bureau. Cannabinoid concentration was from the Drug Enforcement Agency Data. Data was analyzed in R by robust and spatiotemporal regression. RESULTS: In bivariate analyses a dose-response relationship was demonstrated between ALLR and Alcohol Use Disorder (AUD), cocaine and cannabis exposure, with the effect of cannabis being strongest (ß-estimate = 3.33(95%C.I. 1.97, 4.68), P = 1.92 × 10- 6). A strong effect of cannabis use quintile on ALLR was noted (Chi.Sq. = 613.79, P = 3.04 × 10- 70). In inverse probability weighted robust regression adjusted for other substances, income and ethnicity, cannabis was independently significant (ß-estimate = 4.75(0.48, 9.02), P = 0.0389). In a spatiotemporal model adjusted for all drugs, income, and ethnicity, cannabigerol exposure was significant (ß-estimate = 0.26(0.01, 0.52), P = 0.0444), an effect increased by spatial lagging (THC: ß-estimate = 0.47(0.12, 0.82), P = 0.0083). After missing data imputation ethnic cannabis exposure was significant (ß-estimate = 0.64(0.55, 0.72), P = 3.1 × 10- 40). 33/35 minimum e-Values ranged from 1.25 to 3.94 × 1036 indicative of a causal relationship. Relaxation of cannabis legal paradigms had higher ALLR (Chi.Squ.Trend = 775.12, P = 2.14 × 10- 112). Cannabis legal states had higher ALLR (2.395 ± 0.039 v. 2.127 ± 0.008 / 100,000, P = 5.05 × 10- 10). CONCLUSIONS: Data show that ALLR is associated with cannabis consumption across space-time, is associated with the cannabinoids, THC, cannabigerol, cannabinol, cannabichromene, and cannabidiol, contributes to ethnic differentials, demonstrates prominent quintile effects, satisfies criteria for causality and is exacerbated by cannabis legalization.

A geospatiotemporal and causal inference epidemiological exploration of substance and cannabinoid exposure as drivers of rising US pediatric cancer rates.

BACKGROUND: Age-adjusted US total pediatric cancer incidence rates (TPCIR) rose 49% 1975-2015 for unknown reasons. Prenatal cannabis exposure has been linked with several pediatric cancers which together comprise the majority of pediatric cancer types. We investigated whether cannabis use was related spatiotemporally and causally to TPCIR. METHODS: State-based age-adjusted TPCIR data was taken from the CDC Surveillance, Epidemiology and End Results cancer database 2003-2017. Drug exposure was taken from the nationally-representative National Survey of Drug Use and Health, response rate 74.1%. Drugs included were: tobacco, alcohol, cannabis, opioid analgesics and cocaine. This was supplemented by cannabinoid concentration data from the Drug Enforcement Agency and ethnicity and median household income data from US Census. RESULTS: TPCIR rose while all drug use nationally fell, except for cannabis which rose. TPCIR in the highest cannabis use quintile was greater than in the lowest (ß-estimate = 1.31 (95%C.I. 0.82, 1.80), P = 1.80 × 10- 7) and the time:highest two quintiles interaction was significant (ß-estimate = 0.1395 (0.82, 1.80), P = 1.00 × 10- 14). In robust inverse probability weighted additive regression models cannabis was independently associated with TPCIR (ß-estimate = 9.55 (3.95, 15.15), P = 0.0016). In interactive geospatiotemporal models including all drug, ethnic and income variables cannabis use was independently significant (ß-estimate = 45.67 (18.77, 72.56), P = 0.0009). In geospatial models temporally lagged to 1,2,4 and 6 years interactive terms including cannabis were significant. Cannabis interactive terms at one and two degrees of spatial lagging were significant (from ß-estimate = 3954.04 (1565.01, 6343.09), P = 0.0012). The interaction between the cannabinoids THC and cannabigerol was significant at zero, 2 and 6 years lag (from ß-estimate = 46.22 (30.06, 62.38), P = 2.10 × 10- 8). Cannabis legalization was associated with higher TPCIR (ß-estimate = 1.51 (0.68, 2.35), P = 0.0004) and cannabis-liberal regimes were associated with higher time:TPCIR interaction (ß-estimate = 1.87 × 10- 4, (2.9 × 10- 5, 2.45 × 10- 4), P = 0.0208). 33/56 minimum e-Values were > 5 and 6 were infinite. CONCLUSION: Data confirm a close relationship across space and lagged time between cannabis and TPCIR which was robust to adjustment, supported by inverse probability weighting procedures and accompanied by high e-Values making confounding unlikely and establishing the causal relationship. Cannabis-liberal jurisdictions were associated with higher rates of TPCIR and a faster rate of TPCIR increase. Data inform the broader general consideration of cannabinoid-induced genotoxicity.

Co-occurrence across time and space of drug- and cannabinoid- exposure and adverse mental health outcomes in the National Survey of Drug Use and Health: combined geotemporospatial and causal inference analysis.

BACKGROUND: Whilst many studies have linked increased drug and cannabis exposure to adverse mental health (MH) outcomes their effects on whole populations and geotemporospatial relationships are not well understood. METHODS: Ecological cohort study of National Survey of Drug Use and Health (NSDUH) geographically-linked substate-shapefiles 2010-2012 and 2014-2016 supplemented by five-year US American Community Survey. Drugs: cigarettes, alcohol abuse, last-month cannabis use and last-year cocaine use. MH: any mental illness, major depressive illness, serious mental illness and suicidal thinking. DATA ANALYSIS: two-stage, geotemporospatial, robust generalized linear regression and causal inference methods in R. RESULTS: 410,138 NSDUH respondents. Average response rate 76.7%. When drug and sociodemographic variables were combined in geospatial models significant terms including tobacco, alcohol, cannabis exposure and various ethnicities remained in final models for all four major mental health outcomes. Interactive terms including cannabis were related to any mental illness (ß-estimate = 1.97 (95%C.I. 1.56-2.37), P <  2.2 × 10- 16), major depressive episode (ß-estimate = 2.03 (1.54-2.52), P = 3.6 × 10- 16), serious mental illness (SMI, ß-estimate = 2.04 (1.48-2.60), P = 1.0 × 10- 12), suicidal ideation (ß-estimate = 1.99 (1.52-2.47), P <  2.2 × 10- 16) and in each case cannabis alone was significantly associated (from ß-estimate = - 3.43 (- 4.46 - -2.42), P = 3.4 × 10- 11) with adverse MH outcomes on complex interactive regression surfaces. Geospatial modelling showed a monotonic upward trajectory of SMI which doubled (3.62 to 7.06%) as cannabis use increased. Extrapolated to whole populations cannabis decriminalization (4.26%, (4.18, 4.34%)), Prevalence Ratio (PR) = 1.035(1.034-1.036), attributable fraction in the exposed (AFE) = 3.28%(3.18-3.37%), P < 10- 300) and legalization (4.75% (4.65, 4.84%), PR = 1.155 (1.153-1.158), AFE = 12.91% (12.72-13.10%), P < 10- 300) were associated with increased SMI vs. illegal status (4.26, (4.18-4.33%)). CONCLUSIONS: Data show all four indices of mental ill-health track cannabis exposure across space and time and are robust to multivariable adjustment for ethnicity, socioeconomics and other drug use. MH deteriorated with cannabis legalization. Cannabis use-MH data are consistent with causal relationships in the forward direction and include dose-response and temporal-sequential relationships. Together with similar international reports and numerous mechanistic studies preventative action to reduce cannabis use is indicated.

Contemporary epidemiology of rising atrial septal defect trends across USA 1991-2016: a combined ecological geospatiotemporal and causal inferential study.

BACKGROUND: Cardiovascular anomalies are the largest group of congenital anomalies and the major cause of death in young children, with various data linking rising atrial septal defect incidence (ASDI) with prenatal cannabis exposure. Objectives / Hypotheses. Is cannabis associated with ASDI in USA? Is this relationship causal? METHODS: Geospatiotemporal cohort study, 1991-2016. Census populations of adults, babies, congenital anomalies, income and ethnicity. Drug exposure data on cigarettes, alcohol abuse, past month cannabis use, analgesia abuse and cocaine taken from National Survey of Drug Use and Health (78.9% response rate). Cannabinoid concentrations from Drug Enforcement Agency. Inverse probability weighted (ipw) regressions. Analysis conducted in R. RESULTS: ASDI rose nationally three-fold from 27.4 to 82.8 / 10,000 births 1991-2014 during a period when tobacco and alcohol abuse were falling but cannabis was rising. States including Nevada, Kentucky, Mississippi and Tennessee had steeply rising epidemics (Time: Status ß-estimate = 10.72 (95%C.I. 8.39-13.05), P < 2.0 × 10 - 16). ASDI was positively related to exposure to cannabis and most cannabinoids. Drug exposure data was near-complete from 2006 thus restricting spatial modelling from 2006 to 2014, N = 282. In geospatial regression models cannabis: alcohol abuse term was significant (ß-estimate = 19.44 (9.11, 29.77), P = 2.2 × 10 - 4); no ethnic or income factors survived model reduction. Cannabis legalization was associated with a higher ASDI (Time: Status ß-estimate = 0.03 (0.01, 0.05), P = 1.1 × 10 -3). Weighted panel regression interactive terms including cannabis significant (from ß-estimate = 1418, (1080.6, 1755.4), P = 7.3 × 10 -15). Robust generalized linear models utilizing inverse probability weighting interactive terms including cannabis appear (from ß-estimate = 78.88, (64.38, 93.38), P = 1.1 × 10 -8). Marginal structural models with machine-aided SuperLearning association of ASDI with high v. low cannabis exposure R.R. = 1.32 (1.28, 1.36). Model e-values mostly > 1.5. CONCLUSIONS: ASDI is associated with cannabis use, frequency, intensity and legalization in a spatiotemporally significant manner, robust to socioeconomicodemographic adjustment and fulfilled causal criteria, consistent with multiple biological mechanisms and similar reports from Hawaii, Colorado, Canada and Australia. Not only are these results of concern in themselves, but they further imply that our list of the congenital teratology of cannabis is as yet incomplete, and highlight in particular cardiovascular toxicology of prenatal cannabinoid and drug exposure.

Elevation of the ACTH/cortisol ratio in female opioid dependent patients: A biomarker of aging and correlate of metabolic and immune activation.

BACKGROUND: Whilst the hypothalamic-pituitary-adrenal (HPA) Axis is a major stress axis, and is necessarily perturbed in opioid dependency, and stress is a major contributor to aging mechanisms, the HPA axis has not been studied in opioid dependency in an age-dependent manner. OBJECTIVE: Hypothesis - Differences in age dependent levels of HPA components. DESIGN: Cross-sectional comparison of general medical and opioid dependent patients (ODP, GMP). Setting - Primary Care. Patients - 51 GMC, 233 ODP. Ages 37.92+1.95 v. 37.12+0.62 years (P - N.S.) and 33.33% v. 71.67% male (p<0.0001). Intervention(s) - Measurement ACTH, cortisol and their ratio (ACR). Main Outcome Measure(s) - Pre-planned analysis ACR. SECONDARY OUTCOMES: Impact of immune and metabolic markers. RESULTS: ACTH/cortisol was a negative biomarker for age in female patients. Whilst the mean ACR were not different, the (log) ACTH/cortisol showed a positive relationship with age:sex:status (p=0.0396) and age:status (p=0.0437). The effect of addictive status was confined to hepatitis C (HCV) positive female ODP (p=0.0355), and the age:status interaction was also stronger in female HCV+ (p=0.0075) compared to HCV - (p=0.0667) patients. Multiple regression of ACR against age, status, ALT, CRP, and Globulins confirmed many significant interactions. ACTH/cortisol ratio interacted significantly from p=0.0008 in males and p=0.0079 in females, and in both sexes four terms included addictive status. CONCLUSIONS: These data establish the ACTH/cortisol ratio as a negative biomarker of aging in females, and show that this decline is more pronounced in ODP an effect which is partly related to HCV seropositivity, immune and metabolic factors. Dementias are one of the most serious health and socioeconomic issues. Multi-infarct dementia (MID) and Alzheimer´s type dementia (AD) exhibit differences in cerebrovascular blood flow velocity profiles and in presence of microemboli, detected by transcranial Doppler sonography.

Elevated IGF1 in clinical opiate dependence.

OBJECTIVE: To compare IGF1 levels in opiate dependent and general patients both absolutely and by age. DESIGN: A naturalistic observational study was undertaken of opiate dependent and general medical patients. SETTING: Primary care. PATIENTS: 74 opiate substance use dependent (SUD) patients were compared with 262 non-SUD (NSUD) patients. RESULTS: (1) Comparative IGF1 levels; (2) age and sex corrected IGF1 levels; (3) IGF1 levels corrected for age, sex, and hepatic and immune biomarkers. MAIN FINDINGS: The SUD patients were younger (32.60+0.89 vs. 42.49+0.96 years, mean+S.E.M., p<0.0001) and had more males (72.9% and 39.3%, p<0.0001) than the NSUD patients. Restriction of the age range to 15-45 years (70 vs. 153 patients) made the difference in ages non-significant (31.27+0.71 vs. 32.32+0.61 years, p=0.47) but IGF1 remained elevated in SUD (26.56+1.21 vs. 22.65+0.57nmol/L, p=0.0039). When multiple regression was used to correct for the age and sex disparities, the age: addiction interaction remained significantly elevated (p=0.0003). In an additive model opiate dependence showed a 23.8% elevation in IGF1. When the interactive model was further adjusted by the inclusion of ALT and CRP as indices of hepatic inflammation and immune activation respectively, addictive status remained significant both alone (p=0.0134) and in 2-, 3- and 4-way interactions with age, male sex, and ALT (all p<0.0255). CONCLUSION: These data demonstrate that serum IGF1 is elevated in opiate dependence both absolutely and after adjustment for age, sex, and markers of immune and hepatic activation.

Manifold implications of forgotten hyperglycemia in clinical opiate dependence.

Both chronic opiate dependence and diabetes are widespread public health problems. The historically causal relationship between opiate pharmacology and hyperglycemia would appear to have fallen from current medical knowledge. The aims of this study were to assess the severity of hyperglycemia in our opiate substance use disorder (SUD) patient population and review its pathophysiological basis. A review was undertaken of our clinical pathology database, comparing SUD patients and general medical controls (NSUD). A total of 1,602 SUD patients were compared with 2,858 NSUD patients 15-50 years of age. Mean ages were 31.29 ± 7.40 and 31.55 ± 9.17 years, respectively (P = 0.42) and were 68.60 and 52.76% male (P < 0.0001). The glycosylated hemoglobin (HbA1c) level in diabetics was 8.14 ± 2.76 and 6.29 ± 1.88 (P = 0.0032). In a nested case-control subset, random glucose, HbA1c, and fructosamine levels were all significantly elevated (P < 0.05). Adjusted for age, the HbA1c level was elevated among the SUD group (P < 0.0001) the age at addiction interaction was also significant (P < 0.0001). Serum glucose, HbA1c, fructosamine, and microalbumin levels were shown to be biomarkers of age. At multiple regression analysis, interactions between age, addictive status, glycemia, and the inflammatory markers, ethrocyte sedimentation rate and C-reactive protein, were significant (P < 0.05). At age 50, diabetic SUD patients had a HbA1c of 8.31%, which was equivalent to that of NSUD patients at age 140.54 years, a 181.07% elevation. These results confirm that hyperglycemia is higher in both an absolute and an age-dependent sense among diabetic SUD patients. These results are consistent with an acceleration of age-related and degenerative pathologies in many tissues and therefore carry far-reaching implications for the safety of indefinite opiate maintenance.

Perturbation of 3D nuclear architecture, epigenomic dysregulation and aging, and cannabinoid synaptopathy reconfigures conceptualization of cannabinoid pathophysiology: part 1-aging and epigenomics.

Much recent attention has been directed toward the spatial organization of the cell nucleus and the manner in which three-dimensional topologically associated domains and transcription factories are epigenetically coordinated to precisely bring enhancers into close proximity with promoters to control gene expression. Twenty lines of evidence robustly implicate cannabinoid exposure with accelerated organismal and cellular aging. Aging has recently been shown to be caused by increased DNA breaks. These breaks rearrange and maldistribute the epigenomic machinery to weaken and reverse cellular differentiation, cause genome-wide DNA demethylation, reduce gene transcription, and lead to the inhibition of developmental pathways, which contribute to the progressive loss of function and chronic immune stimulation that characterize cellular aging. Both cell lineage-defining superenhancers and the superanchors that control them are weakened. Cannabis exposure phenocopies the elements of this process and reproduces DNA and chromatin breakages, reduces the DNA, RNA protein and histone synthesis, interferes with the epigenomic machinery controlling both DNA and histone modifications, induces general DNA hypomethylation, and epigenomically disrupts both the critical boundary elements and the cohesin motors that create chromatin loops. This pattern of widespread interference with developmental programs and relative cellular dedifferentiation (which is pro-oncogenic) is reinforced by cannabinoid impairment of intermediate metabolism (which locks in the stem cell-like hyper-replicative state) and cannabinoid immune stimulation (which perpetuates and increases aging and senescence programs, DNA damage, DNA hypomethylation, genomic instability, and oncogenesis), which together account for the diverse pattern of teratologic and carcinogenic outcomes reported in recent large epidemiologic studies in Europe, the USA, and elsewhere. It also accounts for the prominent aging phenotype observed clinically in long-term cannabis use disorder and the 20 characteristics of aging that it manifests. Increasing daily cannabis use, increasing use in pregnancy, and exponential dose-response effects heighten the epidemiologic and clinical urgency of these findings. Together, these findings indicate that cannabinoid genotoxicity and epigenotoxicity are prominent features of cannabis dependence and strongly indicate coordinated multiomics investigations of cannabinoid genome-epigenome-transcriptome-metabolome, chromatin conformation, and 3D nuclear architecture. Considering the well-established exponential dose-response relationships, the diversity of cannabinoids, and the multigenerational nature of the implications, great caution is warranted in community cannabinoid penetration.

Perturbation of 3D nuclear architecture, epigenomic aging and dysregulation, and cannabinoid synaptopathy reconfigures conceptualization of cannabinoid pathophysiology: part 2-Metabolome, immunome, synaptome.

The second part of this paper builds upon and expands the epigenomic-aging perspective presented in Part 1 to describe the metabolomic and immunomic bases of the epigenomic-aging changes and then considers in some detail the application of these insights to neurotoxicity, neuronal epigenotoxicity, and synaptopathy. Cannabinoids are well-known to have bidirectional immunomodulatory activities on numerous parts of the immune system. Immune perturbations are well-known to impact the aging process, the epigenome, and intermediate metabolism. Cannabinoids also impact metabolism via many pathways. Metabolism directly impacts immune, genetic, and epigenetic processes. Synaptic activity, synaptic pruning, and, thus, the sculpting of neural circuits are based upon metabolic, immune, and epigenomic networks at the synapse, around the synapse, and in the cell body. Many neuropsychiatric disorders including depression, anxiety, schizophrenia, bipolar affective disorder, and autistic spectrum disorder have been linked with cannabis. Therefore, it is important to consider these features and their complex interrelationships in reaching a comprehensive understanding of cannabinoid dependence. Together these findings indicate that cannabinoid perturbations of the immunome and metabolome are important to consider alongside the well-recognized genomic and epigenomic perturbations and it is important to understand their interdependence and interconnectedness in reaching a comprehensive appreciation of the true nature of cannabinoid pathophysiology. For these reasons, a comprehensive appreciation of cannabinoid pathophysiology necessitates a coordinated multiomics investigation of cannabinoid genome-epigenome-transcriptome-metabolome-immunome, chromatin conformation, and 3D nuclear architecture which therefore form the proper mechanistic underpinning for major new and concerning epidemiological findings relating to cannabis exposure.

Patterns of Cannabis- and Substance-Related Congenital General Anomalies in Europe: A Geospatiotemporal and Causal Inferential Study.

INTRODUCTION: Recent series of congenital anomaly (CA) rates (CARs) have showed the close and epidemiologically causal relationship of cannabis exposure to many CARs. We investigated these trends in Europe where similar trends have occurred. METHODS: CARs from EUROCAT. Drug use from European Monitoring Centre for Drugs and Drug Addiction. Income data from World Bank. RESULTS: CARs were higher in countries with increasing daily use overall (p = 9.99 × 10-14, minimum E-value (mEV) = 2.09) and especially for maternal infections, situs inversus, teratogenic syndromes and VACTERL syndrome (p = 1.49 × 10-15, mEV = 3.04). In inverse probability weighted panel regression models the series of anomalies: all anomalies, VACTERL, foetal alcohol syndrome, situs inversus (SI), lateralization (L), and teratogenic syndromes (TS; AAVFASSILTS) had cannabis metric p-values from: p < 2.2 × 10-16, 1.52 × 10-12, 1.44 × 10-13, 1.88 × 10-7, 7.39 × 10-6 and <2.2 × 10-16. In a series of spatiotemporal models this anomaly series had cannabis metric p-values from: 8.96 × 10-6, 6.56 × 10-6, 0.0004, 0.0019, 0.0006, 5.65 × 10-5. Considering E-values, the cannabis effect size order was VACTERL > situs inversus > teratogenic syndromes > FAS > lateralization syndromes > all anomalies. 50/64 (78.1%) E-value estimates and 42/64 (65.6%) mEVs > 9. Daily cannabis use was the strongest predictor for all anomalies. CONCLUSION: Data confirmed laboratory, preclinical and recent epidemiological studies from Canada, Australia, Hawaii, Colorado and USA for teratological links between cannabis exposure and AAVFASSILTS anomalies, fulfilled epidemiological criteria for causality and underscored importance of cannabis teratogenicity. VACTERL data are consistent with causation via cannabis-induced Sonic Hedgehog inhibition. TS data suggest cannabinoid contribution. SI&L data are consistent with results for cardiovascular CAs. Overall, these data show that cannabis is linked across space and time and in a manner which fulfills epidemiological criteria for causality not only with many CAs, but with several multiorgan teratologic syndromes. The major clinical implication of these results is that access to cannabinoids should be tightly restricted in the interests of safeguarding the community's genetic heritage to protect and preserve coming generations, as is done for all other major genotoxins.

Geospatiotemporal and Causal Inferential Study of European Epidemiological Patterns of Cannabis- and Substance-Related Congenital Orofacial Anomalies.

INTRODUCTION: Since high rates of congenital anomalies (CAs), including facial CAs (FCAs), causally attributed to antenatal and community cannabis use have been reported in several recent series, it was of interest to examine this subject in detail in Europe. METHODS: CA data were taken from the EUROCAT database. Drug exposure data were downloaded from the European Monitoring Centre for Drugs and Drug Addiction (EMCDDA). Income was taken from the World Bank's online sources. RESULTS: On the bivariate maps of both orofacial clefts and holoprosencephaly against resin, the Δ9-tetrahydrocannabinol concentration rates of both covariates increased together in France, Bulgaria, and the Netherlands. In the bivariate analysis, the anomalies could be ranked by the minimum E-value (mEV) as congenital glaucoma > congenital cataract > choanal atresia > cleft lip ± cleft palate > holoprosencephaly > orofacial clefts > ear, face, and neck anomalies. When nations with increasing daily use were compared to those without, the former had generally higher rates of FCAs (p = 0.0281). In the inverse probability weighted panel regression, the sequence of anomalies-orofacial clefts, anotia, congenital cataract, and holoprosencephaly-had positive and significant cannabis coefficients of p = 2.65 × 10-5, 1.04 × 10-8, 5.88 × 10-16, and 3.21 × 10-13, respectively. In the geospatial regression, the same series of FCAs had positive and significant regression terms for cannabis of p = 8.86 × 10-9, 0.0011, 3.36 × 10-8, and 0.0015, respectively. Some 25/28 (89.3%) E-value estimates and 14/28 (50%) mEVs were >9 (considered to be in the high range), and 100% of both were >1.25 (understood to be in the causal range). CONCLUSION: Rising cannabis use is associated with all the FCAs and fulfils the epidemiological criteria for causality. The data indicate particular concerns relating to brain development and exponential genotoxic dose-responses, urging caution with regard to community cannabinoid penetration.

Clinical Epigenomic Explanation of the Epidemiology of Cannabinoid Genotoxicity Manifesting as Transgenerational Teratogenesis, Cancerogenesis and Aging Acceleration.

As global interest in the therapeutic potential of cannabis and its' derivatives for the management of selected diseases increases, it is increasingly imperative that the toxic profile of cannabinoids be thoroughly understood in order to correctly assess the balance between the therapeutic risks and benefits. Modern studies across a number of jurisdictions, including Canada, Australia, the US and Europe have confirmed that some of the most worrying and severe historical reports of both congenital anomalies and cancer induction following cannabis exposure actually underestimate the multisystem thousand megabase-scale transgenerational genetic damage. These findings from teratogenic and carcinogenic literature are supported by recent data showing the accelerated patterns of chronic disease and the advanced DNA methylation epigenomic clock age in cannabis exposed patients. Together, the increased multisystem carcinogenesis, teratogenesis and accelerated aging point strongly to cannabinoid-related genotoxicity being much more clinically significant than it is widely supposed and, thus, of very considerable public health and multigenerational impact. Recently reported longitudinal epigenome-wide association studies elegantly explain many of these observed effects with considerable methodological sophistication, including multiple pathways for the inhibition of the normal chromosomal segregation and DNA repair, the inhibition of the basic epigenetic machinery for DNA methylation and the demethylation and telomerase acceleration of the epigenomic promoter hypermethylation characterizing aging. For cancer, 810 hits were also noted. The types of malignancy which were observed have all been documented epidemiologically. Detailed epigenomic explications of the brain, heart, face, uronephrological, gastrointestinal and limb development were provided, which amply explained the observed teratological patterns, including the inhibition of the key morphogenic gradients. Hence, these major epigenomic insights constituted a powerful new series of arguments which advanced both our understanding of the downstream sequalae of multisystem multigenerational cannabinoid genotoxicity and also, since mechanisms are key to the causal argument, inveighed strongly in favor of the causal nature of the relationship. In this introductory conceptual overview, we present the various aspects of this novel synthetic paradigmatic framework. Such concepts suggest and, indeed, indicate numerous fields for further investigation and basic science research to advance the exploration of many important issues in biology, clinical medicine and population health. Given this, it is imperative we correctly appraise the risk-benefit ratio for each potential cannabis application, considering the potency, severity of disease, stage of human development and duration of use.

Cannabis- and Substance-Related Carcinogenesis in Europe: A Lagged Causal Inferential Panel Regression Study.

Recent European data facilitate an epidemiological investigation of the controversial cannabis-cancer relationship. Of particular concern were prior findings associating high-dose cannabis use with reproductive problems and potential genetic impacts. Cancer incidence data age-standardised to the world population was obtained from the European Cancer Information System 2000-2020 and many European national cancer registries. Drug use data were obtained from the European Monitoring Centre for Drugs and Drug Addiction. Alcohol and tobacco consumption was sourced from the WHO. Median household income was taken from the World bank. Cancer rates in high-cannabis-use countries were significantly higher than elsewhere (ß-estimate = 0.4165, p = 3.54 × 10-115). Eighteen of forty-one cancers (42,675 individual rates) were significantly associated with cannabis exposure at bivariate analysis. Twenty-five cancers were linked in inverse-probability-weighted multivariate models. Temporal lagging in panel models intensified these effects. In multivariable models, cannabis was a more powerful correlate of cancer incidence than tobacco or alcohol. Reproductive toxicity was evidenced by the involvement of testis, ovary, prostate and breast cancers and because some of the myeloid and lymphoid leukaemias implicated occur in childhood, indicating inherited intergenerational genotoxicity. Cannabis is a more important carcinogen than tobacco and alcohol and fulfills epidemiological qualitative and quantitative criteria for causality for 25/41 cancers. Reproductive and transgenerational effects are prominent. These findings confirm the clinical and epidemiological salience of cannabis as a major multigenerational community carcinogen.

Novel Insights into Potential Cannabis-Related Cancerogenesis from Recent Key Whole Epigenome Screen of Cannabis Dependence and Withdrawal: Epidemiological Commentary and Explication of Schrott et al.

Whilst the cannabis-cancer link has been traditionally described as controversial recent whole nation and whole continent studies have demonstrated that well documented laboratory-based multimodal cannabinoid genotoxicity is indeed reflected in numerous cancer types in larger epidemiological series. A recent longitudinal human sperm epigenome-wide DNA methylation screen in both cannabis dependence and cannabis withdrawal has revealed remarkable insights into the manner in which widespread perturbations of DNA methylation may lead to cancerogenic changes in both the exposed and subsequent generations as a result of both cannabis exposure and withdrawal. These results therefore powerfully strengthen and further robustify the causal nature of the relationship between cannabinoid exposure and cancerous outcomes well beyond the previously published extensive mechanistic literature on cannabinoid genotoxicity. The reported epigenomic results are strongly hypothesis generating and call powerfully for further work to investigate oncogenic mechanisms in many tissues, organs and preclinical models. These epigenomic results provide an extraordinarily close predictive account for the epidemiologically observed pattern of cannabis-related malignant disease and indicate that malignant and multigenerational cannabinoid epigenotoxicity is potentially a significant and major public health concern.

European epidemiological patterns of cannabis- and substance-related congenital cardiovascular anomalies: geospatiotemporal and causal inferential study.

As prenatal and community cannabis exposures have recently been linked with congenital heart disease (CHD), it was of interest to explore these associations in Europe in a causal framework and space-time context. Congenital anomaly data from Eurocat, drug-use data from the European Monitoring Centre for Drugs and Drug Addiction, and income from the World Bank. Countries with rising daily cannabis use had in general higher congenital anomaly rates over time than those without (time: status interaction: ß-Est. = 0.0267, P = 0.0059). At inverse probability-weighted panel regression, cannabis terms were positive and significant for CHD, severe CHD, atrial septal defect, ventricular septal defect, atrioventricular septal defect, patent ductus arteriosus, tetralogy of Fallot, vascular disruptions, double outlet right ventricle, transposition of the great vessels, hypoplastic right heart, and mitral valve anomalies from 1.75 × 10-19, 4.20 × 10-11, <2.2 × 10-16, <2.2 × 10-16, 1.58 × 10-12, 4.30 × 10-9, 4.36 × 10-16, 3.50 × 10-8, 5.35 × 10-12, <2.2 × 10-16, 5.65 × 10-5 and 6.06 × 10-10. At spatial regression, terms including cannabis were positive and significant for this same list of anomalies from 0.0038, 1.05 × 10-10, 0.0215, 8.94 × 10-6, 1.23 × 10-5, 2.05 × 10-5, 1.07 × 10-6, 8.77 × 10-5, 9.11 × 10-6, 0.0001, 3.10 × 10-7 and 2.17 × 10-7. 92.6% and 75.2% of 149 E-value estimates and minimum E-values were in high zone >9; 100.0% and 98.7% >1.25. Data show many congenital cardiac anomalies exhibit strong bivariate relationships with metrics of cannabis exposure. Causal inferential modelling for the twelve anomalies selected demonstrated convincing evidence of robust relationships to cannabis which survived adjustment and fulfilled epidemiological criteria for causal relationships. Space-time regression was similarly confirmatory. Epigenomic pathways constitute viable potential mechanisms. Given exponential genotoxic dose-response effects, careful and astute control of cannabinoid penetration is indicated.

Effects of cannabis on congenital limb anomalies in 14 European nations: A geospatiotemporal and causal inferential study.

Cannabinoid exposure is increasing in some European nations. Europe therefore provides an interesting test environment for the recently reported link between cannabis exposure and congenital limb anomaly (CLA) rates (CLARs). Exponential genotoxic dose-response relationships make this investigation both intriguing and imperative. Annual CLAR in 14 nations were from Epidemiological Surveillance of Congenital Anomalies. Drug use rates were from European Monitoring Centre for Drugs and Drug Dependency. Median household income was from the World Bank. E-values provide a quantitative measure of robustness of results to confounding by extraneous covariates. Inverse probability weighting is an important technique for equalizing exposures across countries and removing sources of bias. Rates of CLA, hip dysplasia and the whole group of limb anomalies were higher in countries with increasing daily cannabis use (P = 1.81 × 10-16, 0.0005 and 2.53 × 10-6, respectively). In additive inverse-probability-weighted panel models, the limb reduction-resin Δ9-tetrahydrocannabinol (THC) concentration E-value estimate was 519.93 [95% lower bound (mEV) 49.56], order Resin > Herb ≫ Tobacco > Alcohol. Elevations were noted in 86% E-value estimates and 70.2% of mEVs from 57 E-value pairs from inverse-probability-weighted panel models and from spatial models. As judged by the mEV the degree of association with metrics of cannabis exposure was hip dysplasia > polydactyly > syndactyly > limb anomalies > limb reductions with median E-value estimates from 3.40 × 1065 to 7.06 and median mEVs from 6.14 × 1033 to 3.41. Daily cannabis use interpolated was a more powerful metric of cannabis exposure than herb or resin THC exposure. Data indicate that metrics of cannabis exposure are closely linked with CLAR and satisfy epidemiological criteria for causality. Along with Hawaii and the USA, Europe now forms the third international population in which this causal link has been demonstrated. Cannabis as a predictor of limb anomalies was more potent than tobacco or alcohol. Cannabinoid access should be restricted to protect public health and the community genome/epigenome transgenerationally.