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

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.

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.

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

Introduction. Of the many congenital anomalies (CAs) recently linked with community cannabis exposure, arguably the most concerning are neurological CAs (NCAs). We therefore conducted a detailed study of this in fourteen European nations. Methods. Congenital anomaly data were from Eurocat. Drug exposure data were from European Monitoring Centre for Drugs and Drug Addiction. Income from World bank. Results. The Netherlands, Spain, France and Bulgaria reported increasing rates of many NCAs. The NCA rate (NCAR) was higher in nations with increasing daily cannabis use when compared to those without (p = 0.0204, minimum E-value (mEV) = 1.35). At bivariate analysis, the mEVs of the following NCAs were significantly cannabis related: severe microcephaly 2.14 × 1013, craniosynostosis 5.27 × 1011, nervous system 4.87 × 1011, eye 2.73 × 107, microphthalmos 4.07 × 106, anencephalus 710.37, hydrocephalus 245.64, spina bifida 14.86 and neural tube defects 13.15. At inverse probability, weighted panel regression terms including cannabis were significantly related to the following series of anomalies: nervous system, anencephalus, severe microcephalus, microphthalmos, neural tube defect and spina bifida from p = 5.09 × 10−8, <2.2 × 10−16, <2.2 × 10−16, 4.84 × 10−11, <2.2 × 10−16 and 9.69 × 10−7. At geospatial regression, this same series of anomalies had terms including cannabis significant from p = 0.0027, 1.53 × 10−7, 3.65 × 10−6, 2.13 × 10−8, 0.0002 and 9.76 × 10−12. 88.0% of 50 E-value estimates and 72.0% of mEVs > 9. This analysis therefore demonstrates both close association of cannabis exposure with multiple NCAs across space-time and also fulfills the quantitative criteria of causal inferential analysis. Conclusions. Nine NCARs on bivariate and six NCARs on multivariable regression were cannabis related and fulfilled quantitative epidemiological criteria for causality and are consistent with other series. Particular concerns relate to exponential dose−response effects demonstrated in the laboratory and epidemiological studies. Great caution with community cannabinoid penetration is warranted. Data indicate that cannabis is a significant environmental teratogen and thus imply that cannabinoids should be regulated similarly to the manner in which all other important genotoxins are carefully controlled by communities for their self-sustaining longevity and the protection of generations yet to come.

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.

Socioeconomic, Ethnocultural, Substance- and Cannabinoid-Related Epidemiology of Down Syndrome USA 1986-2016: Combined Geotemporospatial and Causal Inference Investigation

BACKGROUND

Down syndrome (DS) is the commonest of the congenital genetic defects whose incidence has been rising in recent years for unknown reasons. This study aims to assess the impact of substance and cannabinoid use on the DS Rate (DSR) and assess their possible causal involvement.

METHODS

An observational population-based epidemiological study 1986-2016 was performed utilizing geotemporospatial and causal inferential analysis. Participants included all patients diagnosed with DS and reported to state based registries with data obtained from National Birth Defects Prevention Network of Centers for Disease Control. Drug exposure data was from the National Survey of Drug Use and Health (NSDUH) a nationally representative sample interviewing 67,000 participants annually. Drug exposures assessed were: cigarette consumption, alcohol abuse, analgesic/opioid abuse, cocaine use and last month cannabis use. Covariates included ethnicity and median household income from US Census Bureau; maternal age of childbearing from CDC births registries; and cannabinoid concentrations from Drug Enforcement Agency.

RESULTS

NSDUH reports 74.1% response rate. Other data was population-wide. DSR was noted to rise over time and with cannabis use and cannabis-use quintile. In the optimal geospatial model lagged to four years terms including Δ9-tetrahydrocannabinol and cannabigerol were significant (from β-est. = 4189.96 (95%C.I. 1924.74, 6455.17), p = 2.9 × 10-4). Ethnicity, income, and maternal age covariates were not significant. DSR in states where cannabis was not illegal was higher than elsewhere (β-est. = 2.160 (1.5, 2.82), R.R. = 1.81 (1.51, 2.16), p = 4.7 × 10-10). In inverse probability-weighted mixed models terms including cannabinoids were significant (from β-estimate = 18.82 (16.82, 20.82), p < 0.0001). 62 E-value estimates ranged to infinity with median values of 303.98 (IQR 2.50, 2.75 × 107) and 95% lower bounds ranged to 1.1 × 1071 with median values of 10.92 (IQR 1.82, 7990).

CONCLUSIONS

Data show that the association between DSR and substance- and cannabinoid- exposure is robust to multivariable geotemporospatial adjustment, implicate particularly cannabigerol and Δ9-tetrahydrocannabinol, and fulfil quantitative epidemiological criteria for causality. Nevertheless, detailed experimental studies would be required to formally demonstrate causality. Cannabis legalization was associated with elevated DSR's at both bivariate and multivariable analysis. Findings are consistent with those from Hawaii, Colorado, Canada, Australia and Europe and concordant with several cellular mechanisms. Given that the cannabis industry is presently in a rapid growth-commercialization phase the present findings linking cannabis use with megabase scale genotoxicity suggest unrecognized DS risk factors, are of public health importance and suggest that re-focussing the cannabis debate on multigenerational health concerns is prudent.

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 × 1063 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.

Patent Ductus Arteriosus: A Contemporary Perspective for the Pediatric and Adult Cardiac Care Provider

The burden of patent ductus arteriosus (PDA) continues to be significant. In view of marked differences in preterm infants versus more mature, term counterparts (viewed on a continuum with adolescent and adult patients), mechanisms regulating ductal patency, genetic contributions, clinical consequences, and diagnostic and treatment thresholds are discussed separately, when appropriate. Among both preterm infants and older children and adults, a range of hemodynamic profiles highlighting the markedly variable consequences of the PDA are provided. In most contemporary settings, transcatheter closure is preferable over surgical ligation, but data on longer-term outcomes, particularly among preterm infants, are lacking. The present review provides recommendations to identify gaps in PDA diagnosis, management, and treatment on which subsequent research can be developed. Ultimately, the combination of refined diagnostic thresholds and expanded treatment options provides the best opportunities to address the burden of PDA. Although fundamental gaps remain unanswered, the present review provides pediatric and adult cardiac care providers with a contemporary framework in PDA care to support the practice of evidence-based medicine.

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

INTRODUCTION

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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

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

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

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

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSION

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

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSION

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

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.

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.

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

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

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.

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.

Cannabis-related diagnosis in pregnancy and adverse maternal and infant outcomes

BACKGROUND

Cannabis use and cannabis use disorders are increasing in prevalence, including among pregnant women. The objective was to evaluate the association of a cannabis-related diagnosis (CRD) in pregnancy and adverse maternal and infant outcomes.

METHODS

We queried an administrative birth cohort of singleton deliveries in California between 2011-2017 linked to maternal and infant hospital discharge records. We classified pregnancies with CRD from International Classification of Disease codes. We identified nicotine and other substance-related diagnoses (SRD) in the same manner. Outcomes of interest included maternal (hypertensive disorders) and infant (prematurity, small for gestational age, NICU admission, major structural malformations) adverse outcomes.

RESULTS

From 3,067,069 pregnancies resulting in live births, 29,112 (1.0 %) had a CRD. CRD was associated with an increased risk of all outcomes studied; the strongest risks observed were for very preterm birth (aRR 1.4, 95 % CI 1.3, 1.6) and small for gestational age (aRR 1.4, 95 % CI 1.3, 1.4). When analyzed with or without co-exposure diagnoses, CRD alone conferred increased risk for all outcomes compared to no use. The strongest effects were seen for CRD with other SRD (preterm birth aRR 2.3, 95 % CI 2.2, 2.5; very preterm birth aRR 2.6, 95 % CI 2.3, 3.0; gastrointestinal malformations aRR 2.0, 95 % CI 1.6, 2.6). The findings were generally robust to unmeasured confounding and misclassification analyses.

CONCLUSIONS

CRD in pregnancy was associated with increased risk of adverse maternal and infant outcomes. Providing education and effective treatment for women with a CRD during prenatal care may improve maternal and infant health.

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

Background

Ethnic differences in testicular cancer rates (TCRs) are recognized internationally. Cannabis is a known risk factor for testicular cancer (TC) in multiple studies with dose-response effects demonstrated, however the interaction between ancestral and environmental mutagenic effects has not been characterized. We examined the effects of this presumed gene-environment interaction across US states.

Methods

State based TCR was downloaded from the Surveillance Epidemiology and End Results (SEER) website via SEERStat. Drug use data for cigarettes, alcohol use disorder, analgesics, cannabis and cocaine was taken from the National Survey of Drug Use and Health a nationally representative study conducted annually by the Substance Abuse and Mental Health Services Administration (SAMHSA) with a 74.1% response rate. Cannabinoid concentrations derived from Drug Enforcement Agency publications. Median household income and ethnicity data (Caucasian-American, African-American, Hispanic-American, Asian-American, American-Indian-Alaska-Native-American, Native-Hawaiian-Pacific-Islander-American) was from the US Census Bureau. Data were processed in R using instrumental regression, causal inference and multiple imputation.

Results

1975–2017 TCR rose 41% in African-Americans and 78.1% in Caucasian-Americans; 2003–2017 TCR rose 36.1% in Hispanic-Americans and 102.9% in Asian-Pacific-Islander-Americans. Ethnicity-based scatterplot-time and boxplots for cannabis use and TCR closely mirrored each other. At inverse probability-weighted interactive robust regression including drugs, income and ethnicity, ethnic THC exposure was the most significant factor and was independently significant (β-estimate = 4.72 (2.04, 7.41), P = 0.0018). In a similar model THC, and cannabigerol were also significant (both β-estimate = 13.87 (6.33, 21.41), P = 0.0017). In additive instrumental models the interaction of ethnic THC exposure with Asian-American, Hispanic-American, and Native-Hawaiian-Pacific-Islander-American ethnicities was significant (β-estimate = − 0.63 (− 0.74, − 0.52), P = 3.6 × 10^− 29, β-estimate = − 0.25 (− 0.32, − 0.18), P = 4.2 × 10^− 13, β-estimate = − 0.19 (− 0.25, − 0.13), P = 3.4 × 10^− 9). After multiple imputation, ethnic THC exposure became more significant (β-estimate = 0.68 (0.62, 0.74), P = 1.80 × 10^− 92). 25/33 e-Values > 1.25 ranging up to 1.07 × 10⁵. Liberalization of cannabis laws was linked with higher TCR’s in Caucasian-Americans (β-estimate = 0.09 (0.06, 0.12), P = 6.5 × 10^− 10) and African-Americans (β-estimate = 0.22 (0.12, 0.32), P = 4.4 × 10^− 5) and when dichotomized to illegal v. others (t = 6.195, P = 1.18 × 10^− 9 and t = 4.50, P = 3.33 × 10^− 5).

Conclusion

Cannabis is shown to be a TC risk factor for all ethnicities including Caucasian-American and African-American ancestries, albeit at different rates. For both ancestries cannabis legalization elevated TCR. Dose-response and causal relationships are demonstrated.

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.

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

BACKGROUND

Age-adjusted US total pediatric cancer incidence rates (TPCIR) rose 49% 1975-2015 for unknown reasons. Prenatal cannabis exposure has been linked with several pediatric cancers which together comprise the majority of pediatric cancer types. We investigated whether cannabis use was related spatiotemporally and causally to TPCIR.

METHODS

State-based age-adjusted TPCIR data was taken from the CDC Surveillance, Epidemiology and End Results cancer database 2003-2017. Drug exposure was taken from the nationally-representative National Survey of Drug Use and Health, response rate 74.1%. Drugs included were: tobacco, alcohol, cannabis, opioid analgesics and cocaine. This was supplemented by cannabinoid concentration data from the Drug Enforcement Agency and ethnicity and median household income data from US Census.

RESULTS

TPCIR rose while all drug use nationally fell, except for cannabis which rose. TPCIR in the highest cannabis use quintile was greater than in the lowest (β-estimate = 1.31 (95%C.I. 0.82, 1.80), P = 1.80 × 10^- 7) and the time:highest two quintiles interaction was significant (β-estimate = 0.1395 (0.82, 1.80), P = 1.00 × 10^- 14). In robust inverse probability weighted additive regression models cannabis was independently associated with TPCIR (β-estimate = 9.55 (3.95, 15.15), P = 0.0016). In interactive geospatiotemporal models including all drug, ethnic and income variables cannabis use was independently significant (β-estimate = 45.67 (18.77, 72.56), P = 0.0009). In geospatial models temporally lagged to 1,2,4 and 6 years interactive terms including cannabis were significant. Cannabis interactive terms at one and two degrees of spatial lagging were significant (from β-estimate = 3954.04 (1565.01, 6343.09), P = 0.0012). The interaction between the cannabinoids THC and cannabigerol was significant at zero, 2 and 6 years lag (from β-estimate = 46.22 (30.06, 62.38), P = 2.10 × 10^- 8). Cannabis legalization was associated with higher TPCIR (β-estimate = 1.51 (0.68, 2.35), P = 0.0004) and cannabis-liberal regimes were associated with higher time:TPCIR interaction (β-estimate = 1.87 × 10^- 4, (2.9 × 10^- 5, 2.45 × 10^- 4), P = 0.0208). 33/56 minimum e-Values were > 5 and 6 were infinite.

CONCLUSION

Data confirm a close relationship across space and lagged time between cannabis and TPCIR which was robust to adjustment, supported by inverse probability weighting procedures and accompanied by high e-Values making confounding unlikely and establishing the causal relationship. Cannabis-liberal jurisdictions were associated with higher rates of TPCIR and a faster rate of TPCIR increase. Data inform the broader general consideration of cannabinoid-induced genotoxicity.