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

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.

Causal relationship between cannabis use and cancer: a genetically informed perspective

PURPOSE

Cannabis use is increasing legally worldwide, while its impact on cancer risk is unclear. This study was performed to investigate the relationship between cannabis use and the risk of several types of cancer.

METHODS

We conducted a two-sample Mendelian randomization (MR) study to explore the causality of cannabis use on 9 site-specific types of cancer including breast cancer, cervical cancer, melanoma, colorectal cancer, laryngeal cancer, oral cancer, oropharyngeal cancer, esophageal cancer, and glioma. Genome-wide significant genetic instruments (P < 5E-06) for cannabis use were extracted from a large-scale genome-wide association meta-analysis of European ancestry, whereas cancer genetic instruments were extracted from the UK Biobank (UKB) cohort and GliomaScan consortium in the OpenGWAS database. The inverse-variance weighted (IVW) was considered the main method for MR analysis, and sensitivity analyses including MR-Egger, weighted median, MR pleiotropy residual sum, and outlier test (MR-PRESSO) were conducted to evaluate the robustness of the results.

RESULTS

Cannabis use was a significant promoting factor for cervical cancer (OR = 1.001265, 95% CI 1.000375-1.002155, P = 0.0053). And we also detected suggestive evidence of the causality of cannabis use on laryngeal cancer (OR = 1.000350, 95% CI 1.000027-1.000672, P = 0.0336) and breast cancer (OR = 1.003741, 95% CI 1.000052-1.007442, P = 0.0467). No evidence of a causal association of cannabis use with other site-specific types of cancer was detected. Additionally, no pleiotropy or heterogeneity was found in the sensitivity analysis.

CONCLUSION

This study indicates a causative association of cannabis use on cervical cancer, while cannabis use may increase the odds of breast cancer and laryngeal cancer, which require further evaluation in large-scale population-based studies.

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.

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.

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.