Geotemporospatial and causal inferential epidemiological overview and survey of USA cannabis, cannabidiol and cannabinoid genotoxicity expressed in cancer incidence 2003-2017: part 1 - continuous bivariate analysis

Key insights into cannabis-cancer pathobiology and genotoxicity

Whilst mitochondrial inhibition and micronuclear fragmentation are well established features of the cannabis literature mitochondrial stress and dysfunction has recently been shown to be a powerful and direct driver of micronucleus formation and chromosomal breakage by multiple mechanisms. In turn genotoxic damage can be expected to be expressed as increased rates of cancer, congenital anomalies and aging; pathologies which are increasingly observed in modern continent-wide studies. Whilst cannabinoid genotoxicity has long been essentially overlooked it may in fact be all around us through the rapid induction of aging of eggs, sperm, zygotes, foetus and adult organisms with many lines of evidence demonstrating transgenerational impacts. Indeed this multigenerational dimension of cannabinoid genotoxicity reframes the discussion of cannabis legalization within the absolute imperative to protect the genomic and epigenomic integrity of multiple generations to come.

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.

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.

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.

Congenital Gastrointestinal Anomalies in Europe 2010–2019: A Geo-Spatiotemporal and Causal Inferential Study of Epidemiological Patterns in Relationship to Cannabis- and Substance Exposure

Introduction: Congenital anomalies (CA’s) of most of the gastrointestinal tract have been linked causally with prenatal or community cannabis exposure. Therefore, we studied this relationship in Europe. Methods: CA data were from Eurocat. Drug-use data were sourced from the European Monitoring Centre for Drugs and Drug Addiction. Income data were taken from the World Bank. Results: When countries with increasing rates of daily cannabis use were compared with those which were not, the overall rate of gastrointestinal CA’s (GCA’s) was higher in the former group (p = 0.0032). The five anomalies which were related to the metrics of cannabis exposure on bivariate analysis were bile duct atresia, Hirschsprungs, digestive disorders, annular pancreas and anorectal stenosis or atresia. The following sequence of GCA’s was significantly linked with cannabis metrics at inverse-probability-weighted-panel modelling, as indicated: esophageal stenosis or atresia, bile duct atresia, small intestinal stenosis or atresia, anorectal stenosis or atresia, Hirschsprungs disease: p = 1.83 × 10−5, 0.0046, 3.55 × 10−12, 7.35 × 10−6 and 2.00 × 10−12, respectively. When this GCA series was considered in geospatial modelling, the GCA’s were significantly cannabis-related from p = 0.0003, N.S., 0.0086, 6.652 × 10−5, 0.0002, 71.4% of 35 E-value estimates and 54.3% minimum E-values (mEVv’s) > 9 (high zone) and 100% and 97.1% > 1.25 (causality threshold). The order of cannabis sensitivity by median mEVv was Hirschsprungs > esophageal atresia > small intestinal atresia > anorectal atresia > bile duct atresia. Conclusions: Seven of eight GCA’s were related to cannabis exposure and fulfilled the quantitative criteria for epidemiologically causal relationships. Penetration of cannabinoids into the community should be carefully scrutinized and controlled to protect against exponential and multigenerational genotoxicity ensuing from multiple cannabinoids.

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.

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.

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.

Epigenomic and Other Evidence for Cannabis-Induced Aging Contextualized in a Synthetic Epidemiologic Overview of Cannabinoid-Related Teratogenesis and Cannabinoid-Related Carcinogenesis

BACKGROUND

Twelve separate streams of empirical data make a strong case for cannabis-induced accelerated aging including hormonal, mitochondriopathic, cardiovascular, hepatotoxic, immunological, genotoxic, epigenotoxic, disruption of chromosomal physiology, congenital anomalies, cancers including inheritable tumorigenesis, telomerase inhibition and elevated mortality.

METHODS

Results from a recently published longitudinal epigenomic screen were analyzed with regard to the results of recent large epidemiological studies of the causal impacts of cannabis. We also integrate theoretical syntheses with prior studies into these combined epigenomic and epidemiological results.

RESULTS

Cannabis dependence not only recapitulates many of the key features of aging, but is characterized by both age-defining and age-generating illnesses including immunomodulation, hepatic inflammation, many psychiatric syndromes with a neuroinflammatory basis, genotoxicity and epigenotoxicity. DNA breaks, chromosomal breakage-fusion-bridge morphologies and likely cycles, and altered intergenerational DNA methylation and disruption of both the histone and tubulin codes in the context of increased clinical congenital anomalies, cancers and heritable tumors imply widespread disruption of the genome and epigenome. Modern epigenomic clocks indicate that, in cannabis-dependent patients, cannabis advances cellular DNA methylation age by 25-30% at age 30 years. Data have implications not only for somatic but also stem cell and germ line tissues including post-fertilization zygotes. This effect is likely increases with the square of chronological age.

CONCLUSION

Recent epigenomic studies of cannabis exposure provide many explanations for the broad spectrum of cannabis-related teratogenicity and carcinogenicity and appear to account for many epidemiologically observed findings. Further research is indicated on the role of cannabinoids in the aging process both developmentally and longitudinally, from stem cell to germ cell to blastocystoids to embryoid bodies and beyond.

Letter to the editor on “Geotemporospatial and causal inferential epidemiological overview and survey of USA cannabis, cannabidiol and cannabinoid genotoxicity expressed in cancer incidence 2003–2017, parts 1–3”

We would like to thank authors Reece and Hulse (2022) for their three-part article titled “Geotemporospatial and causal inferential epidemiological overview and survey of USA cannabis, cannabidiol and cannabinoid genotoxicity expressed in cancer incidence 2003-2017”, in which the authors infer that cannabis use has a causal role in the development of various cancer types. While the authors use reputable datasets and a well-established epidemiological methodology, the authors’ conclusion of a causal association is limited due to biases inherent in ecological epidemiological studies. Though the researchers attempt to overcome these biases through validation and statistical manipulations, their approaches are insufficient to create conditions suitable for causal inferencing upon examination. There are also concerns in the practical and conceptual application of the studies’ dataset that further question the validity of the authors’ inferences. Further research exploring the potential benefits and harm of cannabinoids in the context of cancer must be performed before a distinct relationship can be defined.

Epidemiological Patterns of Cannabis- and Substance- Related Congenital Uronephrological Anomalies in Europe: Geospatiotemporal and Causal Inferential Study

INTRODUCTION

Recent reports linking prenatal and community cannabis exposure to elevated uronephrological congenital anomaly (UCA) rates (UCAR's) raise the question of its European epidemiology given recent increases in community cannabinoid penetration there.

METHODS

UCAR data from Eurocat. Drug use data from European Monitoring Centre for Drugs and Drug Addiction. Income from World bank.

RESULTS

UCAR increased across Spain, Netherlands, Poland and France. UCAR's and cannabis resin THC increased simultaneously in France, Spain, Netherlands and Bulgaria. At bivariate analysis all UCA's were related to cannabis herb and resin THC concentrations. All UCAR's were bivariately related to cannabis metrics ordered by median minimum E-value (mEV) as hypospadias > multicystic renal disease > bilateral renal agenesis > UCA's > hydronephrosis > posterior urethral valve > bladder exstrophy/epispadias. At inverse probability weighted multivariable analysis terms including cannabis were significant for the following series of anomalies: UCA's, multicystic renal disease, bilateral renal agenesis, hydronephrosis, congenital posterior urethral valves from P = 1.91 × 10-5, 2.61 × 10-8, 4.60 × 10-15, 4.60 × 10-15 and 2.66 × 10-10. At geospatial analysis the same series of UCA's were significantly related to cannabis from P = 7.84 × 10-15, 7.72 × 10-5, 0.0023, 6.95 × 10-5, and 8.82 × 10-5. 45/51 (88.2%) of E-value estimates and 31/51 (60.8%) of mEV's >9.

CONCLUSION

Analysis confirms a close relationship between cannabis metrics and all seven UCA's and fulfill formal criteria for quantitative causal inference. Given the exponential cannabinoid genotoxicity dose-response relationship results provide a powerful stimulus to constrain community cannabinoid exposure including protection of the food chain to preserve the genome and epigenome of coming generations.

Cannabis- and Substance-Related Epidemiological Patterns of Chromosomal Congenital Anomalies in Europe: Geospatiotemporal and Causal Inferential Study

INTRODUCTION

Laboratory data link cannabinoid exposure to chromosomal mis-segregation errors. Recent epidemiological reports confirm this link and raise concern that elevated chromosomal congenital anomaly rates (CCAR) may be occurring in Europe which is experiencing increased cannabis use, daily intensity of use and cannabinoid potency.

METHODS

CCAR data from Eurocat. Drug use data from the European Monitoring Centre for Drugs and Drug Addiction. Income from World Bank. Bivariate, multivariate, panel and geotemporospatial regressions analyzed. Inverse probability weighting of panel models and E-values used as major quantitative causal inferential methodologies.

RESULTS

In countries where daily cannabis use was rising the trend for CCA's was upwards whereas in those where daily use was declining it was usually downwards (p = 0.0002). In inverse probability weighted panel models terms for cannabis metrics were significant for chromosomal disorders, trisomies 21 and 13 and Klinefelters syndrome from p < 2.2 × 10-16. In spatiotemporal models cannabis terms were positive and significant for chromosomal disorders, genetic disorders, trisomies 21, 18 and 13, Turners and Klinefelters syndromes from 4.28 × 10-6, 5.79 × 10-12, 1.26 × 10-11, 1.12 × 10-7, 7.52 × 10-9, 7.19 × 10-7 and 7.27 × 10-7. 83.7% of E-value estimates and 74.4% of minimum E-values (mEV) > 9 including four values each at infinity. Considering E-values: the sensitivity of the individual disorders was trisomy 13 > trisomy 21 > Klinefelters > chromosomal disorders > Turners > genetic syndromes > trisomy 18 with mEV's 1.91 × 1025 to 59.31; and daily cannabis use was the most powerful covariate (median mEV = 1.91 × 1025).

CONCLUSIONS

Data indicate that, consistent with reports from Hawaii, Canada, Colorado, Australia and USA, CCARs are causally and spatiotemporally related to metrics and intensity of cannabis exposure, directly impact 645 MB (21.5%) of the human genome and may implicate epigenomic-centrosomal mechanisms.

European Epidemiological Patterns of Cannabis- and Substance-Related Body Wall Congenital Anomalies: Geospatiotemporal and Causal Inferential Study

As body wall congenital anomalies (BWCAs) have a long history of being associated with prenatal or community cannabis exposure (CCE), it was of interest to investigate these epidemiological relationships in Europe given the recent increases in cannabis use prevalence, daily intensity, and Δ9-tetrahydrocannabinol (THC) potency. Methods: This study makes use of BWCA data from Eurocat, drug exposure data from the European Monitoring Centre for Drugs and Drug Addiction, and income from the World Bank. Results: The mapping analysis showed that BWCARs increased in France, Spain, and the Netherlands. The bivariate mapping analysis showed that the BWCA rates (BWCAR) and the cannabis resin THC concentration rose simultaneously in France, the Netherlands, Bulgaria, Sweden, and Norway. The bivariate ranking of the BWCARs by median minimum E-value (mEV) was omphalocele > diaphragmatic hernia > abdominal wall defects > gastroschisis. With inverse probability weighted multivariable panel regression, the series of BWCAs, including gastroschisis, omphalocele, and diaphragmatic hernia, was positively related to various metrics of cannabis use from p = 2.45 × 10−14, 4.77 × 10−7 and <2.2 × 10−16. With geospatial regression, the same series of BWCAs was related to cannabis metrics from p = 0.0016, 5.28 × 10−6 and 4.88 × 10−9. Seventeen out of twenty-eight (60.7%) of the E-value estimates were >9 (high range), as were 14/28 (50.0%) of the mEVs. Conclusion: The data confirm the close relationship of the BWCARs with the metrics of CCE, fulfill the quantitative criteria of causal inference, and underscore the salience of the public health impacts of cannabinoid teratogenicity. Of major concern is the rising CCE impacting exponential cannabinoid genotoxic dose-response relationships. CCE should be carefully restricted to protect the food chain, the genome, and the epigenome of coming generations.

Epidemiology of Δ8THC-Related Carcinogenesis in USA: A Panel Regression and Causal Inferential Study

The use of Δ8THC is increasing at present across the USA in association with widespread cannabis legalization and the common notion that it is "legal weed". As genotoxic actions have been described for many cannabinoids, we studied the cancer epidemiology of Δ8THC. Data on 34 cancer types was from the Centers for Disease Control Atlanta Georgia, substance abuse data from the Substance Abuse and Mental Health Services Administration, ethnicity and income data from the U.S. Census Bureau, and cannabinoid concentration data from the Drug Enforcement Agency, were combined and processed in R. Eight cancers (corpus uteri, liver, gastric cardia, breast and post-menopausal breast, anorectum, pancreas, and thyroid) were related to Δ8THC exposure on bivariate testing, and 18 (additionally, stomach, Hodgkins, and Non-Hodgkins lymphomas, ovary, cervix uteri, gall bladder, oropharynx, bladder, lung, esophagus, colorectal cancer, and all cancers (excluding non-melanoma skin cancer)) demonstrated positive average marginal effects on fully adjusted inverse probability weighted interactive panel regression. Many minimum E-Values (mEVs) were infinite. p-values rose from 8.04 × 10-78. Marginal effect calculations revealed that 18 Δ8THC-related cancers are predicted to lead to a further 8.58 cases/100,000 compared to 7.93 for alcoholism and -8.48 for tobacco. Results indicate that between 8 and 20/34 cancer types were associated with Δ8THC exposure, with very high effect sizes (mEVs) and marginal effects after adjustment exceeding tobacco and alcohol, fulfilling the epidemiological criteria of causality and suggesting a cannabinoid class effect. The inclusion of pediatric leukemias and testicular cancer herein demonstrates heritable malignant teratogenesis.

Geotemporospatial and causal inferential epidemiological overview and survey of USA cannabis, cannabidiol and cannabinoid genotoxicity expressed in cancer incidence 2003-2017: part 2 - categorical bivariate analysis and attributable fractions

BACKGROUND

As the cannabis-cancer relationship remains an important open question epidemiological investigation is warranted to calculate key metrics including Rate Ratios (RR), Attributable Fractions in the Exposed (AFE) and Population Attributable Risks (PAR) to directly compare the implicated case burden between emerging cannabinoids and the established carcinogen tobacco.

METHODS

SEER*Stat software from Centres for Disease Control was used to access age-standardized state census incidence of 28 cancer types (including "All (non-skin) Cancer") from National Cancer Institute in US states 2001-2017. Drug exposures taken from the National Survey of Drug Use and Health 2003-2017, response rate 74.1%. Federal seizure data provided cannabinoid exposure. US Census Bureau furnished income and ethnicity. Exposure dichotomized as highest v. lowest exposure quintiles. Data processed in R.

RESULTS

Nineteen thousand eight hundred seventy-seven age-standardized cancer rates were returned. Based on these rates and state populations this equated to 51,623,922 cancer cases over an aggregated population 2003-2017 of 124,896,418,350. Fifteen cancers displayed elevated E-Values in the highest compared to the lowest quintiles of cannabidiol exposure, namely (in order): prostate, melanoma, Kaposi sarcoma, ovarian, bladder, colorectal, stomach, Hodgkins, esophagus, Non-Hodgkins lymphoma, All cancer, brain, lung, CLL and breast. Eleven cancers were elevated in the highest THC exposure quintile: melanoma, thyroid, liver, AML, ALL, pancreas, myeloma, CML, breast, oropharynx and stomach. Twelve cancers were elevated in the highest tobacco quintile confirming extant knowledge and study methodology. For cannabidiol RR declined from 1.397 (95%C.I. 1.392, 1.402), AFE declined from 28.40% (28.14, 28.66%), PAR declined from 15.3% (15.1, 15.5%) and minimum E-Values declined from 2.13. For THC RR declined from 2.166 (95%C.I. 2.153, 2.180), AFE declined from 53.8% (53.5, 54.1%); PAR declined from 36.1% (35.9, 36.4%) and minimum E-Values declined from 3.72. For tobacco, THC and cannabidiol based on AFE this implies an excess of 93,860, 91,677 and 48,510 cases; based on PAR data imply an excess of 36,450, 55,780 and 14,819 cases.

CONCLUSION

Data implicate 23/28 cancers as being linked with THC or cannabidiol exposure with epidemiologically-causal relationships comparable to those for tobacco. AFE-attributable cases for cannabinoids (91,677 and 48,510) compare with PAR-attributable cases for tobacco (36,450). Cannabinoids constitute an important multivalent community carcinogen.

Geotemporospatial and causal inferential epidemiological overview and survey of USA cannabis, cannabidiol and cannabinoid genotoxicity expressed in cancer incidence 2003-2017: part 3 - spatiotemporal, multivariable and causal inferential pathfinding and exploratory analyses of prostate and ovarian cancers

BACKGROUND

The epidemiology of cannabinoid-related cancerogenesis has not been studied with cutting edge epidemiological techniques. Building on earlier bivariate papers in this series we aimed to conduct pathfinding studies to address this gap in two tumours of the reproductive tract, prostate and ovarian cancer.

METHODS

Age-standardized cancer incidence data for 28 tumour types (including "All (non-skin) Cancer") was sourced from Centres for Disease Control and National Cancer Institute using SEER*Stat software across US states 2001-2017. Drug exposure was sourced from the nationally representative household survey National Survey of Drug Use and Health conducted annually by the Substance Abuse and Mental Health Services Administration 2003-2017 with response rate 74.1%. Federal seizure data provided cannabinoid concentration data. US Census Bureau provided income and ethnicity data. Inverse probability weighted mixed effects, robust and panel regression together with geospatiotemporal regression analyses were conducted in R. E-Values were also calculated.

RESULTS

19,877 age-standardized cancer rates were returned. Based on these rates and state populations this equated to 51,623,922 cancer cases over an aggregated population 2003-2017 of 124,896,418,350. Inverse probability weighted regressions for prostate and ovarian cancers confirmed causal associations robust to adjustment. Cannabidiol alone was significantly associated with prostate cancer (β-estimate = 1.61, (95%C.I. 0.99, 2.23), P = 3.75 × 10- 7). In a fully adjusted geospatiotemporal model at one spatial and two temporal years lags cannabidiol was significantly independently associated with prostate cancer (β-estimate = 2.08, (1.19, 2.98), P = 5.20 × 10- 6). Cannabidiol alone was positively associated with ovarian cancer incidence in a geospatiotemporal model (β-estimate = 0.36, (0.30, 0.42), P <  2.20 × 10- 16). The cigarette: THC: cannabidiol interaction was significant in a fully adjusted geospatiotemporal model at six years of temporal lag (β-estimate = 1.93, (1.07, 2.78), P = 9.96 × 10- 6). Minimal modelled polynomial E-Values for prostate and ovarian cancer ranged up to 5.59 × 1059 and 1.92 × 10125. Geotemporospatial modelling of these tumours showed that the cannabidiol-carcinogenesis relationship was supra-linear and highly sigmoidal (P = 1.25 × 10- 45 and 12.82 × 10- 52 for linear v. polynomial models).

CONCLUSION

Cannabinoids including THC and cannabidiol are therefore important community carcinogens additive to the effects of tobacco and greatly exceeding those of alcohol. Reproductive tract carcinogenesis necessarily implies genotoxicity and epigenotoxicity of the germ line with transgenerational potential. Pseudoexponential and causal dose-response power functions are demonstrated.

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 × 10⁶⁵ to 7.06 and median mEVs from 6.14 × 10³³ 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.