Causal inference multiple imputation investigation of the impact of cannabinoids and other substances on ethnic differentials in US testicular cancer incidence

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.

Cannabis and Paternal Epigenetic Inheritance

Cannabis is the most widely used illicit drug in Western counties and its abuse is particularly high in male adolescents and young adults. Its main psychotropic component, the cannabinoid delta-9-tetrahydrocannabinol (Δ⁹-THC), interferes in the endogenous endocannabinoid system. This signaling system is involved in the control of many biological activities, including the formation of high-quality male gametes. Direct adverse effects of Δ⁹-THC in male reproduction are well known in both animal models and humans. Nevertheless, the possibility of long-term effects due to epigenetic mechanisms has recently been reported. In this review, we summarize the main advances in the field suggesting the need to pay attention to the possible long-term epigenetic risks for the reproductive health of cannabis users and the health of their offspring.

Endocannabinoid system and epigenetics in spermatogenesis and testicular cancer

In mammals, male germ cell development starts during fetal life and is carried out in postnatal life with the formation of sperms. Spermatogenesis is the complex and highly orderly process during which a group of germ stem cells is set at birth, starts to differentiate at puberty. It proceeds through several stages: proliferation, differentiation, and morphogenesis and it is strictly regulated by a complex network of hormonal, autocrine and paracrine factors and it is associated with a unique epigenetic program. Altered epigenetic mechanisms or inability to respond to these factors can impair the correct process of germ development leading to reproductive disorders and/or testicular germ cell cancer. Among factors regulating spermatogenesis an emerging role is played by the endocannabinoid system (ECS). ECS is a complex system comprising endogenous cannabinoids (eCBs), their synthetic and degrading enzymes, and cannabinoid receptors. Mammalian male germ cells have a complete and active ECS which is modulated during spermatogenesis and that crucially regulates processes such as germ cell differentiation and sperm functions. Recently, cannabinoid receptor signaling has been reported to induce epigenetic modifications such as DNA methylation, histone modifications and miRNA expression. Epigenetic modifications may also affect the expression and function of ECS elements, highlighting the establishment of a complex mutual interaction. Here, we describe the developmental origin and differentiation of male germ cells and testicular germ cell tumors (TGCTs) focusing on the interplay between ECS and epigenetic mechanisms involved in these processes.

Do We Have Structure, Process and Outcomes to Support Cannabis as Supportive Therapy in Cancer?

Cannabis is becoming more popular and more available in the United States. It has been approved for use by multiple states for various conditions and several states now allow recreational cannabis. We explore the structure of cannabis distribution, the process of acquisition, outcomes, and the safety of cannabis in the United States.

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.

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.

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.

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 × 10²⁵ to 59.31; and daily cannabis use was the most powerful covariate (median mEV = 1.91 × 10²⁵).

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⁻⁷⁸. 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 1 - continuous bivariate analysis

Background

The genotoxic and cancerogenic impacts of population-wide cannabinoid exposure remains an open but highly salient question. The present report examines these issues from a continuous bivariate perspective with subsequent reports continuing categorical and detailed analyses.

Methods

Age-standardized state census incidence of 28 cancer types (including “All (non-skin) Cancer”) was sourced using SEER*Stat software from Centres for Disease Control and National Cancer Institute across US states 2001–2017. It was joined with drug exposure data from the nationally representative National Survey of Drug Use and Health conducted annually by the Substance Abuse and Mental Health Services Administration 2003–2017, response rate 74.1%. Cannabinoid data was from Federal seizure data. Income and ethnicity data sourced from the US Census Bureau. Data was 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. Regression lines were charted for cancer-substance exposures for cigarettes, alcohol use disorder (AUD), cannabis, THC, cannabidiol, cannabichromene, cannabinol and cannabigerol. In this substance series positive trends were found for 14, 9, 6, 9, 12, 6, 9 and 7 cancers; with largest minimum E-Values (mEV) of 1.76 × 10⁹, 4.67 × 10⁸, 2.74 × 10⁴, 4.72, 2.34 × 10¹⁸, 2.74 × 10¹⁷, 1.90 × 10⁷, 5.05 × 10⁹; and total sum of exponents of mEV of 34, 32, 13, 0, 103, 58, 25, 31 indicating that cannabidiol followed by cannabichromene are the most strongly implicated in environmental carcinogenesis. Breast cancer was associated with tobacco and all cannabinoids (from mEV = 3.53 × 10⁹); “All Cancer” (non-skin) linked with cannabidiol (mEV = 1.43 × 10¹¹); pediatric AML linked with cannabis (mEV = 19.61); testicular cancer linked with THC (mEV = 1.33). Cancers demonstrating elevated mEV in association with THC were: thyroid, liver, pancreas, AML, breast, oropharynx, CML, testis and kidney. Cancers demonstrating elevated mEV in relation to cannabidiol: prostate, bladder, ovary, all cancers, colorectum, Hodgkins, brain, Non-Hodgkins lymphoma, esophagus, breast and stomach.

Conclusion

Data suggest that cannabinoids including THC and cannabidiol are important community carcinogens exceeding the effects of tobacco or alcohol. Testicular, (prostatic) and ovarian tumours indicate mutagenic corruption of the germline in both sexes; pediatric tumourigenesis confirms transgenerational oncogenesis; quantitative criteria implying causality are fulfilled.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13690-022-00811-8.

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 × 10⁵⁹ and 1.92 × 10¹²⁵. 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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13690-022-00813-6.

Geospatiotemporal and causal inference study of cannabis and other drugs as risk factors for female breast cancer USA 2003-2017

Breast cancer (BC) is the commonest human cancer and its incidence (BC incidence, BCI) is rising worldwide. Whilst both tobacco and alcohol have been linked to BCI genotoxic cannabinoids have not been investigated. Age-adjusted state-based BCI 2003-2017 was taken from the Surveillance Epidemiology and End Results database of the Centers for Disease Control. Drug use from the National Survey of Drug Use and Health, response rate 74.1%. Median age, median household income and ethnicity were from US census. Inverse probability weighted (ipw) multivariable regression conducted in R. In bivariate analysis BCI was shown to be significantly linked with rising cannabis exposure {β-est. = 3.93 [95% confidence interval 2.99, 4.87], P = 1.10 × 10^-15}. At 8 years lag cigarettes:cannabis [β-est. = 2660 (2150.4, 3169.3), P = 4.60 × 10^-22] and cannabis:alcoholism [β-est. = 7010 (5461.6, 8558.4), P = 1.80 × 10^-17] were significant in ipw-panel regression. Terms including cannabidiol [CBD; β-est. = 16.16 (0.39, 31.93), P = 0.446] and cannabigerol [CBG; β-est. = 6.23 (2.06, 10.39), P = 0.0034] were significant in spatiotemporal models lagged 1:2 years, respectively. Cannabis-liberal paradigms had higher BCI [67.50 ± 0.26 v. 65.19 ± 0.21/100 000 (mean ± SEM), P = 1.87 × 10^-11; β-est. = 2.31 (1.65, 2.96), P = 9.09 × 10^-12]. 55/58 expected values >1.25 and 13/58 >100. Abortion was independently and causally significant in space-time models. Data show that exposure to cannabis and the cannabinoids Δ9-tetrahydrocannabinol, CBD, CBG and alcoholism fulfil quantitative causal criteria for BCI across space and time. Findings are robust to adjustment for age and several known sociodemographic, socio-economic and hormonal risk factors and establish cannabinoids as an additional risk factor class for breast carcinogenesis. BCI is higher under cannabis-liberal legal paradigms.

Cannabinoid and substance relationships of European congenital anomaly patterns: a space-time panel regression and causal inferential study

With reports from Australia, Canada, USA, Hawaii and Colorado documenting a link between cannabis and congenital anomalies (CAs), this relationship was investigated in Europe. Data on 90 CAs were accessed from Eurocat. Tobacco and alcohol consumption and median household income data were from the World Bank. Amphetamine, cocaine and last month and daily use of cannabis from the European Monitoring Centre for Drugs and Drug Addiction. Cannabis herb and resin Δ9-tetrahydrocannabinol concentrations were from published reports. Data were processed in R. Twelve thousand three hundred sixty CA rates were sourced across 16 nations of Europe. Nations with a higher or increasing rate of daily cannabis use had a 71.77% higher median CA rates than others [median ± interquartile range 2.13 (0.59, 6.30) v. 1.24 (0.15, 5.14)/10 000 live births (P = 4.74 × 10^-17; minimum E-value (mEV) = 1.52]. Eighty-nine out of 90 CAs in bivariate association and 74/90 CAs in additive panel inverse probability weighted space-time regression were cannabis related. In inverse probability weighted interactive panel models lagged to zero, two, four and six years, 76, 31, 50 and 29 CAs had elevated mEVs (< 2.46 × 10³⁹) for cannabis metrics. Cardiovascular, central nervous, gastrointestinal, genital, uronephrology, limb, face and chromosomalgenetic systems along with the multisystem VACTERL syndrome were particularly vulnerable targets. Data reveal that cannabis is related to many CAs and fulfil epidemiological criteria of causality. The triple convergence of rising cannabis use prevalence, intensity of daily use and Δ9-tetrahydrocannabinol concentration in herb and resin is powerfully implicated as a primary driver of European teratogenicity, confirming results from elsewhere.

Congenital anomaly epidemiological correlates of Δ8THC across USA 2003-16: panel regression and causal inferential study

Δ8-Tetrahydrocannabinol (Δ8THC) is marketed in many US states as 'legal weed'. Concerns exist relating to class-wide genotoxic cannabinoid effects. We conducted an epidemiological investigation of Δ8THC-related genotoxicity expressed as 57 congenital anomaly (CA) rates (CARs) in the USA. CARs were taken from the Centers for Disease Control, Atlanta, Georgia. Drug exposure data were taken from the National Survey of Drug Use and Health, with a response rate of 74.1%. Ethnicity and income data were taken from the US Census Bureau. National cannabinoid exposure was taken from Drug Enforcement Agency publications and multiplied by state cannabis use data to derive state-based estimates of Δ8THC exposure. At bivariate continuous analysis, Δ8THC was associated with 23 CAs on raw CA rates, 33 CARs after correction for early termination for anomaly estimates and 41 on a categorical analysis comparing the highest and lowest exposure quintiles. At inverse probability weighted multivariable additive and interactive models lagged to 0, 2 and 4 years, Δ8THC was linked with 39, 8, 4 and 9 CAs. Chromosomal, cardiovascular, gastrointestinal, genitourinary, limb, central nervous system (CNS) and face systems were particularly affected. The minimum E-values ranged to infinity. Both the number of anomalies implicated and the effect sizes demonstrated were much greater for Δ8THC than for tobacco and alcohol combined. Δ8THC appears epidemiologically to be more strongly associated with many CAs than for tobacco and alcohol and is consistent with a cannabinoid class genotoxic/epigenotoxic effect. Quantitative causality criteria were fulfilled, and causal relationships either for Δ8THC or for cannabinoid/s, for which it is a surrogate marker, may be in operation.