Pathways from epigenomics and glycobiology towards novel biomarkers of addiction and its radical cure

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.

State Trends of Cannabis Liberalization as a Causal Driver of Increasing Testicular Cancer Rates across the USA

BACKGROUND

The cause of the worldwide doubling-tripling of testicular cancer rates (TCRs) in recent decades is unknown. Previous cohort studies associated cannabis use with TCR including dose-response relationships but the contribution of cannabis to TCRs at the population level is unknown. This relationship was tested by analyzing annual trends across US states and formally assessed causality. Four US datasets were linked at state level: age-adjusted TCRs from Centers for Disease Control Surveillance Epidemiology and End Results database; drug use data from annual National Survey of Drug Use and Health including 74.1% response rate; ethnicity and median household income data from the US Census Bureau; and cannabinoid concentration data from Drug Enforcement Agency reports. Data was processed in R in spatiotemporal and causal inference protocols.

RESULTS

Cannabis-use quintile scatterplot-time and boxplots closely paralleled those for TCRs. The highest cannabis-use quintile had a higher TCR than others (3.44 ± 0.05 vs. 2.91 ± 0.2, mean ± S.E.M., t = 10.68, p = 1.29 × 10^-22). A dose-response relationship was seen between TCR and Δ9-tetrahydrocannabinol (THC), cannabinol, cannabigerol, and cannabichromene (6.75 × 10^-9 < p < 1.83 × 10^-142). In a multivariate inverse probability-weighted interactive regression including race and ethnic cannabis exposure (ECE), ECE was significantly related to TCR (β-estimate = 0.89 (95%C.I. 0.36, 2.67), p < 2.2 × 10^-16). In an additive geospatiotemporal model controlling for other drugs, cannabis alone was significant (β-estimate = 0.19 (0.10, 0.28), p = 3.4 × 10^-5). In a full geospatial model including drugs, income and ethnicity cannabinoid exposure was significant (cannabigerol: β-estimate = 1.39 (0.024, 2.53), p = 0.0017); a pattern repeated at two spatial and two temporal lags (cannabigerol: β-estimate = 0.71 (0.05, 1.37), p = 0.0.0350; THC: β-estimate = 23.60 (11.92, 35.29), p = 7.5 × 10^-5). 40/41 e-Values > 1.25 ranged up to 1.4 × 10⁶³ and 10 > 1000 fitting causal relationship criteria. Cannabis liberalization was associated with higher TCRs (ChiSqu. = 312.2, p = 2.64 × 10^-11). Rates of TC in cannabis-legal states were elevated (3.36 ± 0.09 vs. 3.01 ± 0.03, t = 4.69, p = 4.86 × 10^-5).

CONCLUSIONS

Cannabis use is closely and causally associated with TCRs across both time and space and higher in States with liberal cannabis legislation. Strong dose-response effects were demonstrated for THC, cannabigerol, cannabinol, cannabichromene and cannabidiol. Cannabinoid genotoxicity replicates all major steps to testicular carcinogenesis including whole-genome doubling, chromosomal arm excision, generalized DNA demethylation and chromosomal translocations thereby accelerating the pathway to testicular carcinogenesis by several decades.

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.

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

Background

Cannabinoids including cannabidiol have recognized genotoxic activities but their significance has not been studied broadly epidemiologically across the teratological spectrum. We examined these issues including contextual space-time relationships and formal causal inferential analysis in USA.

Methods

State congenital anomaly (CA) rate (CAR) data was taken from the annual reports of the National Birth Defects Prevention Network 2001–2005 to 2011–2015. Substance abuse rates were from the National Survey of Drug Use and Health a nationally representative longitudinal survey of the non-institutionalized US population with 74.1% response rate. Drugs examined were cigarettes, monthly and binge alcohol, monthly cannabis and analgesic and cocaine abuse. Early termination of pregnancy for abortion (ETOPFA) rates were taken from the published literature. Cannabinoid concentrations were from Drug Enforcement Agency. Ethnicity and income data were from the US Census Bureau. Inverse probability weighted (IPW) regressions and geotemporospatial regressions conducted for selected CAs.

Results

Data on 18,328,529 births from an aggregated population of 2,377,483,589 for mid-year analyses 2005–2013 comprehending 12,611 CARs for 62 CAs was assembled and ETOPFA-corrected (ETOPFACAR) where appropriate. E-Values for ETOPFACARs by substance trends were elevated for THC (40 CAs), cannabis (35 CAs), tobacco (11 CAs), cannabidiol (8 CAs), monthly alcohol (5 CAs) and binge alcohol (2 CAs) with minimum E-Values descending from 16.55, 1.55x107, 555.10, 7.53x1019, 9.30 and 32.98. Cardiovascular, gastrointestinal, chromosomal, limb reductions, urinary, face and body wall CAs particularly affected. Highest v. lowest substance use quintile CAR prevalence ratios 2.84 (95%C.I. 2.44, 3.31), 4.85 (4.08, 5.77) and 1.92 (1.63, 2.27) and attributable fraction in exposed 0.28 (0.27, 0.28), 0.57 (0.51, 0.62) and 0.47 (0.38, 0.55) for tobacco, cannabis and cannabidiol. Small intestinal stenosis or atresia and obstructive genitourinary defect were studied in detail in lagged IPW pseudo-randomized causal regressions and spatiotemporal models confirmed the causal role of cannabinoids. Spatiotemporal predictive modelling demonstrated strongly sigmoidal non-linear cannabidiol dose-response power-function relationships (P = 2.83x10−60 and 1.61x10−71 respectively).

Conclusions

Data implicate cannabinoids including cannabidiol in a diverse spectrum of heritable CAs. Sigmoidal non-linear dose-response relationships are of grave concern.

These transgenerational genotoxic, epigenotoxic, chromosomal-toxic putatively causal teratogenic effects strongly indicate tight restrictions on community cannabinoid penetration.

Cannabinoid 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.

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

BACKGROUND

Acute lymphoid leukaemia (ALL) is the commonest childhood cancer whose incidence is rising in many nations. In the USA, between 1975 and 2016, ALL rates (ALLRs) rose 93.51% from 1.91 to 3.70/100,000 <  20 years. ALL is more common in Caucasian-Americans than amongst minorities. The cause of both the rise and the ethnic differential is unclear, however, prenatal cannabis exposure was previously linked with elevated childhood leukaemia rates. We investigated epidemiologically if cannabis use impacted nationally on ALLRs, its ethnic effects, and if the relationship was causal.

METHODS

State data on overall, and ethnic ALLR from the Surveillance Epidemiology and End Results databank of the Centre for Disease Control (CDC) and National Cancer Institute (NCI) were combined with drug (cigarettes, alcoholism, cannabis, analgesics, cocaine) use data from the National Survey of Drug Use and Health; 74.1% response rate. Income and ethnicity data was from the US Census bureau. Cannabinoid concentration was from the Drug Enforcement Agency Data. Data was analyzed in R by robust and spatiotemporal regression.

RESULTS

In bivariate analyses a dose-response relationship was demonstrated between ALLR and Alcohol Use Disorder (AUD), cocaine and cannabis exposure, with the effect of cannabis being strongest (β-estimate = 3.33(95%C.I. 1.97, 4.68), P = 1.92 × 10^- 6). A strong effect of cannabis use quintile on ALLR was noted (Chi.Sq. = 613.79, P = 3.04 × 10^- 70). In inverse probability weighted robust regression adjusted for other substances, income and ethnicity, cannabis was independently significant (β-estimate = 4.75(0.48, 9.02), P = 0.0389). In a spatiotemporal model adjusted for all drugs, income, and ethnicity, cannabigerol exposure was significant (β-estimate = 0.26(0.01, 0.52), P = 0.0444), an effect increased by spatial lagging (THC: β-estimate = 0.47(0.12, 0.82), P = 0.0083). After missing data imputation ethnic cannabis exposure was significant (β-estimate = 0.64(0.55, 0.72), P = 3.1 × 10^- 40). 33/35 minimum e-Values ranged from 1.25 to 3.94 × 10³⁶ indicative of a causal relationship. Relaxation of cannabis legal paradigms had higher ALLR (Chi.Squ.Trend = 775.12, P = 2.14 × 10^- 112). Cannabis legal states had higher ALLR (2.395 ± 0.039 v. 2.127 ± 0.008 / 100,000, P = 5.05 × 10^- 10).

CONCLUSIONS

Data show that ALLR is associated with cannabis consumption across space-time, is associated with the cannabinoids, THC, cannabigerol, cannabinol, cannabichromene, and cannabidiol, contributes to ethnic differentials, demonstrates prominent quintile effects, satisfies criteria for causality and is exacerbated by cannabis legalization.

Epidemiological overview of multidimensional chromosomal and genome toxicity of cannabis exposure in congenital anomalies and cancer development

Cannabis and cannabinoids are implicated in multiple genotoxic, epigenotoxic and chromosomal-toxic mechanisms and interact with several morphogenic pathways, likely underpinning previous reports of links between cannabis and congenital anomalies and heritable tumours. However the effects of cannabinoid genotoxicity have not been assessed on whole populations and formal consideration of effects as a broadly acting genotoxin remain unexplored. Our study addressed these knowledge gaps in USA datasets. Cancer data from CDC, drug exposure data from National Survey of Drug Use and Health 2003-2017 and congenital anomaly data from National Birth Defects Prevention Network were used. We show that cannabis, THC cannabigerol and cannabichromene exposure fulfill causal criteria towards first Principal Components of both: (A) Down syndrome, Trisomies 18 and 13, Turner syndrome, Deletion 22q11.2, and (B) thyroid, liver, breast and pancreatic cancers and acute myeloid leukaemia, have mostly medium to large effect sizes, are robust to adjustment for ethnicity, other drugs and income in inverse probability-weighted models, show prominent non-linear effects, have 55/56 e-Values > 1.25, and are exacerbated by cannabis liberalization (P = 9.67 × 10^-43, 2.66 × 10^-15). The results confirm experimental studies showing that cannabinoids are an important cause of community-wide genotoxicity impacting both birth defect and cancer epidemiology at the chromosomal hundred-megabase level.

Epidemiological overview of multidimensional chromosomal and genome toxicity of cannabis exposure in congenital anomalies and cancer development

Cannabis and cannabinoids are implicated in multiple genotoxic, epigenotoxic and chromosomal-toxic mechanisms and interact with several morphogenic pathways, likely underpinning previous reports of links between cannabis and congenital anomalies and heritable tumours. However the effects of cannabinoid genotoxicity have not been assessed on whole populations and formal consideration of effects as a broadly acting genotoxin remain unexplored. Our study addressed these knowledge gaps in USA datasets. Cancer data from CDC, drug exposure data from National Survey of Drug Use and Health 2003-2017 and congenital anomaly data from National Birth Defects Prevention Network were used. We show that cannabis, THC cannabigerol and cannabichromene exposure fulfill causal criteria towards first Principal Components of both: (A) Down syndrome, Trisomies 18 and 13, Turner syndrome, Deletion 22q11.2, and (B) thyroid, liver, breast and pancreatic cancers and acute myeloid leukaemia, have mostly medium to large effect sizes, are robust to adjustment for ethnicity, other drugs and income in inverse probability-weighted models, show prominent non-linear effects, have 55/56 e-Values > 1.25, and are exacerbated by cannabis liberalization (P = 9.67 × 10^-43, 2.66 × 10^-15). The results confirm experimental studies showing that cannabinoids are an important cause of community-wide genotoxicity impacting both birth defect and cancer epidemiology at the chromosomal hundred-megabase level.

Canadian Cannabis Consumption and Patterns of Congenital Anomalies: An Ecological Geospatial Analysis

OBJECTIVES

Cannabis is a known teratogen. Data availability addressing both major congenital anomalies and cannabis use allowed us to explore their geospatial relationships.

METHODS

Data for the years 1998 to 2009 from Canada Health and Statistics Canada was analyzed in R. Maps have been drawn and odds ratios, principal component analysis, correlation matrices, least squares regression and geospatial regression analyses have been conducted using the R packages base, dplyr, epiR, psych, ggplot2, colorplaner and the spml and spreml functions from package splm.

RESULTS

Mapping showed cannabis use was more common in the northern Territories of Canada in the Second National Survey of Cannabis Use 2018. Total congenital anomalies, all cardiovascular defects, orofacial clefts, Downs syndrome and gastroschisis were all found to be more common in these same regions and rose as a function of cannabis exposure. When Canada was dichotomized into high and low cannabis use zones by Provinces v Territories the Territories had a higher rate of total congenital anomalies 450.026 v 390.413 (O.R. = 1.16 95%C.I. 1.08-1.25, P = 0.000058; attributable fraction in exposed 13.25%, 95%C.I. 7.04-19.04%). In geospatial analysis in a spreml spatial error model cannabis was significant both alone as a main effect (P < 2.0 × 10) and in all its first and second order interactions with both tobacco and opioids from P < 2.0 × 10.

CONCLUSION

These results show that the northern Territories of Canada share a higher rate of cannabis use together with elevated rates of total congenital anomalies, all cardiovascular defects, Down's syndrome and gastroschisis. This is the second report of a significant association between cannabis use and both total defects and all cardiovascular anomalies and the fourth published report of a link with Downs syndrome and thereby direct major genotoxicity. The correlative relationships described in this paper are confounded by many features of social disadvantage in Canada's northern territories. However, in the context of a similar broad spectrum of defects described both in animals and in epidemiological reports from Hawaii, Colorado, USA and Australia they are cause for particular concern and indicate further research.

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

BACKGROUND

Cardiovascular anomalies are the largest group of congenital anomalies and the major cause of death in young children, with various data linking rising atrial septal defect incidence (ASDI) with prenatal cannabis exposure. Objectives / Hypotheses. Is cannabis associated with ASDI in USA? Is this relationship causal?

METHODS

Geospatiotemporal cohort study, 1991-2016. Census populations of adults, babies, congenital anomalies, income and ethnicity. Drug exposure data on cigarettes, alcohol abuse, past month cannabis use, analgesia abuse and cocaine taken from National Survey of Drug Use and Health (78.9% response rate). Cannabinoid concentrations from Drug Enforcement Agency. Inverse probability weighted (ipw) regressions. Analysis conducted in R.

RESULTS

ASDI rose nationally three-fold from 27.4 to 82.8 / 10,000 births 1991-2014 during a period when tobacco and alcohol abuse were falling but cannabis was rising. States including Nevada, Kentucky, Mississippi and Tennessee had steeply rising epidemics (Time: Status β-estimate = 10.72 (95%C.I. 8.39-13.05), P < 2.0 × 10 ^- 16). ASDI was positively related to exposure to cannabis and most cannabinoids. Drug exposure data was near-complete from 2006 thus restricting spatial modelling from 2006 to 2014, N = 282. In geospatial regression models cannabis: alcohol abuse term was significant (β-estimate = 19.44 (9.11, 29.77), P = 2.2 × 10 ^- 4); no ethnic or income factors survived model reduction. Cannabis legalization was associated with a higher ASDI (Time: Status β-estimate = 0.03 (0.01, 0.05), P = 1.1 × 10 ^-3). Weighted panel regression interactive terms including cannabis significant (from β-estimate = 1418, (1080.6, 1755.4), P = 7.3 × 10 ^-15). Robust generalized linear models utilizing inverse probability weighting interactive terms including cannabis appear (from β-estimate = 78.88, (64.38, 93.38), P = 1.1 × 10 ^-8). Marginal structural models with machine-aided SuperLearning association of ASDI with high v. low cannabis exposure R.R. = 1.32 (1.28, 1.36). Model e-values mostly > 1.5.

CONCLUSIONS

ASDI is associated with cannabis use, frequency, intensity and legalization in a spatiotemporally significant manner, robust to socioeconomicodemographic adjustment and fulfilled causal criteria, consistent with multiple biological mechanisms and similar reports from Hawaii, Colorado, Canada and Australia. Not only are these results of concern in themselves, but they further imply that our list of the congenital teratology of cannabis is as yet incomplete, and highlight in particular cardiovascular toxicology of prenatal cannabinoid and drug exposure.

Broad Spectrum epidemiological contribution of cannabis and other substances to the teratological profile of northern New South Wales: geospatial and causal inference analysis

Background

Whilst cannabis commercialization is occurring rapidly guided by highly individualistic public narratives, evidence that all congenital anomalies (CA) increase alongside cannabis use in Canada, a link with 21 CA’s in Hawaii, and rising CA’s in Colorado indicate that transgenerational effects can be significant and impact public health. It was therefore important to study Northern New South Wales (NNSW) where cannabis use is high.

Methods

Design: Cohort. 2008–2015. Setting: NNSW and Queensland (QLD), Australia. Participants. Whole populations. Exposures. Tobacco, alcohol, cannabis. Source: National Drug Strategy Household Surveys 2010, 2013. Main Outcomes. CA Rates. NNSW-QLD comparisons. Geospatial and causal regression.

Results

Cardiovascular, respiratory and gastrointestinal anomalies rose with falling tobacco and alcohol but rising cannabis use rates across Queensland. Maternal age NNSW-QLD was not different (2008–2015: 4265/22084 v. 96,473/490514 > 35 years/total, Chi.Sq. = 1.687, P = 0.194). A higher rate of NNSW cannabis-related than cannabis-unrelated defects occurred (prevalence ratio (PR) = 2.13, 95%C.I. 1.80–2.52, P = 3.24 × 10^− 19). CA’s rose more potently with rising cannabis than with rising tobacco or alcohol use. Exomphalos and gastroschisis had the highest NNSW:QLD PR (6.29(2.94–13.48) and 5.85(3.54–9.67)) and attributable fraction in the exposed (84.11%(65.95–92.58%) and 82.91%(71.75–89.66%), P = 2.83 × 10^− 8 and P = 5.62 × 10^− 15). In multivariable geospatial models cannabis was significantly linked with cardiovascular (atrial septal defect, ventricular septal defect, tetralogy of Fallot, patent ductus arteriosus), genetic (chromosomal defects, Downs syndrome), gastrointestinal (small intestinal atresia), body wall (gastroschisis, diaphragmatic hernia) and other (hypospadias) (AVTPCDSGDH) CA’s. In linear modelling cannabis use was significantly linked with anal stenosis, congenital hydrocephalus and Turner syndrome (ACT) and was significantly linked in borderline significant models (model P < 0.1) with microtia, microphthalmia, and transposition of the great vessels. At robust and mixed effects inverse probability weighted multivariable regression cannabis was related to 18 defects. 16/17 E-Values in spatial models were > 1.25 ranging up to 5.2 × 10¹³ making uncontrolled confounding unlikely.

Conclusions

These results suggest that population level CA’s react more strongly to small rises in cannabis use than tobacco or alcohol; cardiovascular, chromosomal, body wall and gastrointestinal CA’s rise significantly with small increases in cannabis use; that cannabis is a bivariate correlate of AVTPCDSGDH and ACT anomalies, is robust to adjustment for other substances; and is causal.

Supplementary information

Supplementary information accompanies this paper at 10.1186/s40360-020-00450-1.