Next Article in Journal
Clinical and Laboratory Manifestation of Gastrointestinal Involvement in MIS-C: A Single-Center Observational Study
Previous Article in Journal
Glycemic Abnormalities in Pancreatic Cystic Lesions—A Single-Center Retrospective Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sociodemographically Stratified Exploration of Pancreatic Cancer Incidence in Younger US Patients: Implication of Cannabis Exposure as a Risk Factor

by
Albert Stuart Reece
1,2,* and
Gary Kenneth Hulse
1,2
1
Division of Psychiatry, University of Western Australia, Crawley, WA 6009, Australia
2
School of Medical and Health Sciences, Edith Cowan University, Joondalup, WA 6027, Australia
*
Author to whom correspondence should be addressed.
Gastroenterol. Insights 2023, 14(2), 204-235; https://doi.org/10.3390/gastroent14020016
Submission received: 1 March 2023 / Revised: 1 May 2023 / Accepted: 7 May 2023 / Published: 11 May 2023
(This article belongs to the Section Pancreas)

Abstract

:
Introduction. The aetiology for the recent increase in pancreatic cancer incidence (PCI) in the US is unknown. This paper provides an epidemiological investigation of the exponential increase in PCI in young people aged 15–34 years, particularly amongst females, with a focus on the exponential rise amongst African American females, and its relationship to substance use. Methods. National pancreatic cancer data from recent reports. Tobacco, alcohol and daily cannabis use data taken from the annual nationally representative National Survey of Drug Use and Health, response rate = 74%. Results. Amongst the 15–34-year-aged cohort, PCI was found to be significantly more common in females (females: β-est. = 0.1749 p = 0.0005). African American females are noted to have the highest rates of daily cannabis use amongst females in the 26–34 and 35–49-year groups. The relationship between PCI and daily cannabis use was strongly positive across all ethnicities and in both sexes. In African American females, the Pearson correlation between daily cannabis use and PCI was R = 0.8539, p = 0.0051. In an additive multivariable model for each sex and race, cannabis was the only significant term remaining in the final model in the 15–34-year-aged cohort and thus out-performed alcohol as a risk factor. The most significant term in multivariate models was the alcohol:cannabis interaction which was highly significant in all ethnicities from p = 2.50 × 10−7 for Caucasian American females and the highest E-value pair was for Hispanic American females (E-value estimate = 1.26 × 10102 and E-value lower bound 2.20 × 1074). Conclusion. These data show that cannabis fulfills quantitative criteria of causality in all age, sex and ethnicity cohorts, and thus explains both the recent surge in PCI and its ethnocentric predominance. Cannabis interacts powerfully genotoxically and cancerogenically with alcohol, with increases in cannabis use driving the current PCI surge. These results raise the important question as to how much cannabis might be responsible for the modern renaissance in cancer rates amongst younger people.

1. Introduction

Pancreatic cancer is a common and serious condition with a five-year survival rate of only 12.5% in 2013–2019 [1]. Pancreatic cancer is becoming more common in many places [1] for reasons which are presently obscure. Indeed, the rising incidence of many different cancers occurring in younger people has been noted across many nations [2]. This phenomenon begins to challenge the medical dogma that cancer is primarily a disease of older age. A detailed review of the whole USA caseload of pancreatic cancer 2001–2018 was recently published numbering in total around 734,761 cases [3]. Very concerningly, it showed that the disease is growing quickly across the USA and is more common in females in many age cohorts. Indeed, in some groups, such as African American females aged 15–34 years and for localized disease, exponential effects were noted. Some of the curves drawn in that report suggest that a break point occurred in pancreatic incidence around 2006–2009. The authors suggest that an environmental intoxicant may be acting whose effects were most marked amongst females [3] and it appears to have become more widespread about the time of the decadal change.
In this context, recent reports from the USA and Europe that cannabis is associated with pancreatic cancer [4,5,6,7] become both interesting and highly relevant. Moreover, as it is well known that cannabis use is being popularized in the USA and its daily use has doubled in that country in recent years [8,9], it is theoretically possible that the use of a known carcinogen is becoming widespread and commonplace and indeed normalized. Interestingly, the role of cannabis in hepatocarcinogenic changes is increasingly being realized and the American Association of Liver Diseases has recently released a statement to this effect [10]. This conclusion solidly implicates cannabinoids in gastrointestinal system oncogenic pathways.
Several powerful epigenomic studies of cannabis exposure in humans and rodents have recently appeared [11,12,13,14,15,16] which have greatly increased our understanding of cannabinoid pathophysiology. Global hypomethylation of DNA has been reported by many workers [15,17,18] and this is a change which is characteristic of aging [19,20,21] and also has the effect of removing the control of mobile elements from the genome which is an oncogenic and pro-aging change [22,23]. Many cannabinoids (including cannabidiol) cause single- and double-stranded DNA breaks [24,25,26,27,28,29,30] and this perturbation has been shown to cause aging [21] by re-arranging the epigenomic machinery and causing DNA hypomethylation [21]. Since epigenomic markers are concentrated at tissue defining superenhancers, and since superenhancers play a critical and irreplaceable role in cell lineage definition and control of differentiation state [31,32,33,34,35], this major change weakens cell identity and primes them towards de-differentiation and malignant transformation in the pancreas and also in many other tissues [36,37,38,39,40,41,42]. Moreover, cannabis exposure also weakens the CTCF boundary markers [17] including superanchors which control active gene transcription generally and superenhancer activation in particular, which further disrupts the normal control of gene expression. Many cancers are caused by superenhancers gaining aberrant access to the promoter regions of oncogenes such as KRAS, a process which is common in pancreatic neocarcinogenesis [43,44,45].
Modern multi-channel sequencing methods are increasingly being applied in order to advance our understanding of the pathogenesis of pancreatic cancer [45,46,47,48]. It has recently been shown that pancreatic acinar cells remember past episodes of inflammation epigenomically which is a change which both allows them to mobilize their genes more rapidly the next time they are required for a similar episode and moves them along the cancerogenic spectrum to the point where a widespread premalignant field change can be induced [45,48]. This disruption will be exacerbated by the known frequently pro-inflammatory [49,50,51,52,53,54,55] and generally pro-aging [56,57,58] effect of cannabinoids.
This is believed to be the underlying mechanism by which known predisposing factors such as gall stones and alcohol act to raise the risk of pancreatic tumours which is through chronic or recurrent inflammation in that organ [45,48].
It should be noted in passing that cancerogenesis is just one part of the broader subject of cannabinoid genotoxicity more generally. The use of significantly genotoxic compounds may be expected to be reflected in higher rates of cancer, congenital anomalies (birth defects) and aging. Increased evidence of such outcomes has appeared in the recent literature [7,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82] which becomes important background information to the present study and is highly relevant to the consideration of causal issues. The rising rates of cannabis use have recently been shown to be driving breast cancer, acute lymphoid leukaemia, hepatocellular carcinoma, testicular cancer and total paediatric cancer in the USA [10,83,84,85,86,87,88,89,90].
Accordingly, this study tested the hypothesis that the rising rate of high intensity daily cannabis use in the USA may be an important driver of higher rates of pancreatic cancer. We particularly wanted to test this hypothesis in population subgroups defined by age, sex and ethnicity. The main focus of this study was on the pancreatic cancer incidence in younger patients of 15–34 years of age and in localized disease, which are the two groups shown in the Abboud study where the disease appeared to be growing fastest and exponentially [3]. These dual hypotheses were formulated prior to the commencement of this study. This study also represents a quantitative follow-up of an epidemiological investigation which has recently been suggested [91].

2. Methods

2.1. Pancreatic Cancer Rates

Pancreatic cancer rates were taken from two of the graphs in the report of Abboud et. al. [6] which were Figure 4C for males and females aged 15–34 and Figure 6C for localized tumours presenting in males and females. The senior study author disclosed that data had been provided to them on a confidential basis only and they were not permitted to release it for secondary analysis (Gaddam S, personal communication). The data were taken from the graphs very carefully using highly exploded images of the graphs using WebDigitizer [92] which is a highly accurate platform for performing such tasks. These two sex-specific rates were log transformed in the interests of improving compliance with normality assumptions as indicated by the Shapiro–Wilks test.

2.2. Drug Use Rates

The data on national levels of drug use by age group, sex and ethnicity were derived from the Substance Abuse and Mental Health Services Administration (SAMHSA) 2020 National Survey of Drug Use and Health (NSDUH) [11]. This is a nationally representative study of the non-institutionalized US population conducted each year by SAMHSA. Consideration was restricted to the major ethnicities Caucasian American (“White”), African American (“Black”) and Hispanic American (“Hispanic). The annual results of NSDUH appear as a Public Data Archive System (PDAS) on the Substance Abuse and Mental Health Data Archive (SAMHDA) and can be searched manually. The major terms of interest were for age (CATAG3), ethnicity and sex combined (SEXRACE) and then substance. The three substances of interest and their metrics were for tobacco, the percentage which reported tobacco consumption in the previous month (TOBMON); for alcohol, the percentage which reported alcohol consumption in the previous month (ALCMON); and for cannabis, the percentage reporting daily or near daily use for 20–30 days per month (MRJMDAYS). In each case, cross-tabulations were performed with the sex–race metric as the columns, substance exposure as the rows and the age category as the control variable. In the analysis, these three substance exposure rates were not log transformed as indicated by the Shapiro–Wilks test.

2.3. Data Analysis

The data analysis was performed in February 2023 in R Studio v 2002.12.0 Build 353 based on R version 4.2.2 (2002-10-31 ucrt)). Multiple regression and correlation matrices were constructed using the R Base and stats modules. Graphs were drawn in ggplot2 from the tidyverse [93]. Multiple graphs were arranged using the R packages ggpubr and patchwork [94,95]. Multiple regression models were reduced by the classical technique of removal of the least significant term. p-values were adjusted for multiple testing using the false discovery rate or adjustment method described by Holm. Multiple linear or exponential models were analysed simultaneously using the purrr-broom workflow described in the tidyverse [93,96,97]. E-Values were calculated using the package EValue in R [98,99,100]. Correlograms including correlograms of significance levels were drawn with the R Package corrplot [101]. p-values less than 0.05 were considered significant.

2.4. Data Extrapolation

Data on the rate of pancreas cancer incidence were available from the source paper by Abboud from 2001 to 2018 [3]. Drug use data from SAMHDA were available from 2002 to 2019. To match up this slight discrepancy in the time periods of the two datasets, an extra year was added in 2001 to the SAMSHA dataset and in 2019 for the pancreatic cancer dataset. Extrapolation was performed with the predict function on the log model as detailed in the accompanying R Code. Each dataset was processed separately (pancreatic cancer 15–34 years, localized pancreatic disease, tobacco, alcohol and daily cannabis use rates).

3. Results

3.1. Input Data

Pancreatic cancer (PC) rates were taken from the figures in the paper by Abboud and colleagues as described [3]. Our main focus of interest was on the exponentially rising cancer rates in younger patients 15–34 years of age. The sex differential of this rise was first published in 2021 from the SEER dataset [102]. Abboud et.al. extended this finding by confirming the sex differential in the data of the US Cancer Statistics (USCS) dataset without the data from states contributing to the SEER data [3]. Whilst both the USCS and SEER data are available publicly, the data used by the Abboud group which omits data from the states contributing the SEER data are not publicly available. This is why the data employed had to be derived from the figures of the later publication [3]. The sex differential of the three datasets are compared directly in Supplementary Figure S1 (present publication) and are found to be qualitatively similar across all three datasets. Data on disease incidence by stage are not publicly available so this also had to be similarly derived from the published figures.
Drug use data were taken from online databases at NSDUH as discussed in the Methods section. This survey has a reported completion rate of 74% [103].
The rates of PC by sex in younger patients of 15–34-years-of-age are shown graphically in Figure 1 on linear and log scales with both quadratic and loess curves fitted. As suggested in these figures, the rates were significantly different between the two sexes when sex was considered in an additive exponential model with time (females: β-est. = 0.1749, t = 3.838, p = 0.0005; model Adj.R.Squ. = 0.90, F = 1598.1, df = 2.33, model p < 2.2 × 10−16). When sex was considered as an interactive term with time, both sex as a factor (females: β-est. = −35.99, t = −2.148, p = 0.0394; model Adj.R.Squ. = 0.9103, F = 119.4, df = 3.32, model p < 2.2 × 10−16) and in a sex:time interaction (Females: β-est. = 0.0180, t = 2.159, p = 0.0385) were significant.
Figure 2 shows the rates of tobacco consumption by the ethnicity, gender and age cohorts over time and a general decline is observed consistent with the overall trend to reduced tobacco consumption in the community.
Figure 3 is a similar representation of alcohol consumption. This pattern also shows an earlier rise followed by a later decline in all ethnicities and age cohorts.
Daily or near daily cannabis use trends are shown in Figure 4 by age, sex and ethnicity. Strong rises are noted in all sub-groups. African American females are noted to have the highest rates of daily cannabis use amongst females in the 26–34 year group and in the 35–49 year group. Indeed, amongst African American females in this age cohort, daily or near-daily use rose from 3.11% in 2002 to 9.59% in 2019 which represents more than a tripling in rate (308% rise). This compares to rises amongst Hispanic American females from 1.05% to 5.02% (334%) and in Caucasian American females from 3.19% to 7.72% (242%).
The rises in near daily cannabis use in males in this age cohort across this period for these three ethnicities were for African Americans 7.08% to 15.9% (224%), Caucasian Americans 7.51% to 12.8% (170%) and for Hispanic Americans 3.45% to 10.1% (292%).
In Figure 5, NSDUH respondents have been dichotomized into patients older than 50 years and those who are younger. There is a clear trend to increased substance use for all types in younger patients. This difference is particularly marked in the case of cannabis where monotonic rises are also observed (Figure 5E,F).
Figure 6 sets out the relationship of pancreatic cancer incidence (PCI) to the various substance exposures. The relationship of PCI with tobacco appears to be inverse (lower panels). The relationship with alcohol is apparently non-descript (upper panels). However, for all ethnicities and in both sexes the relationship with daily cannabis use is clearly strongly positively correlated (middle panels).

3.2. Correlations

Figure 7 is a correlogram showing the correlation between the two different pancreatic cancer rates (15–34 years and localized disease (HOP) shown in third and fourth rows) and tobacco, alcohol and cannabis in males of a Caucasian background 18–34 years. The left-hand panel depicts the correlation coefficients and the panel on the right shows the statistical significance. The left-hand panel indicates the strong correlation of both pancreatic cancer and localized disease with cannabis exposure (Pearson R = 0.83, 0.70) which are both significant (p = 0.006 and 0.023), respectively. Negative correlations of pancreatic cancer with both alcohol and tobacco are observed. A weak and non-significant positive correlation between daily cannabis and alcohol use is noted here (R = 0.1657, p = 0.1369).
In African American males of 18–34 years, a similar pattern is observed (Figure 8). In this case, the Pearson correlation coefficients for pancreatic cancer and for localized disease are 0.83 and 0.64 which are again significant (0.007 and 0.038).
In Hispanic American males, a very similar pattern is again observed (Supplementary Figure S2).
Figure 9 shows the situation in white females 15–34 years of age. There, the correlation between PCI and daily cannabis exposure is 0.9309 and between localized disease and cannabis is 0.9040. These correlations are significant at the p = 0.0008 and 0.0005 levels, respectively.
When African American females are considered, there is again a strong relationship (Figure 10). The Pearson correlation coefficients for PCI and localized disease are 0.8539 and 0.8421 which are significant at the p = 0.0051 and 0.0056 levels, respectively. In this figure, a weak correlation and non-significant correlation between daily cannabis use and alcohol use is noted (R = 0.59, p = 0.68).
Similar patterns are observed when Hispanic American females are considered (Supplementary Figure S3). The correlation in this group with daily cannabis exposure is 0.8982 and 0.9133 which are significant at the 0.0008 and 0.0005 levels, respectively.
Thus, similar changes are observed in all age and ethnicity sub-groups.
The exact values of these various correlation coefficients are shown in tabular form in Table 1 and their corresponding significance levels are given in Table 2. If one compares these tables closely, it emerges that the only consistent positive correlations of interest are between cannabis and PCI and localized disease. Most of these correlations are 0.8–0.9 which indicate statistical significance mostly in the range 0.001–0.01.

3.3. Bivariate Regressions

Table 3 sets out the results of the regression studies for PCI against the three substances of interest considered separately. The rate of pancreatic cancer which was regressed as the dependent variable was the sex-specific rate taken from the Abboud paper [3]. The main positive and significant terms in this table relate to cannabis in both sexes and alcohol in females. If one compares the size of the β-estimates, they are clearly much greater for cannabis (19–32 in females) than alcohol (5–6 in females). Interestingly, the relationship between sex and alcohol seems to be somewhat different.
These various regression coefficients are associated with applicable E-Values as indicated in the Table. The E-Value estimates are very high for cannabis in females (2 × 1028 to 3.46 × 1053) and also in males (3.74 × 1021 to 2.97 × 1038). Whilst the E-Values for alcohol are also elevated, they are much less elevated than those for cannabis.
Indeed, considering the lower bounds of these E-Values, those for cannabis are all strongly positive and highly significant (6.88 × 1019 to 1.20 × 1042 in females and 2.53 × 1015 to 4.73 × 1030 in males). Two of the four lower E-Values for alcohol are significant in this table, and they are both for females.
Thus, the uniform finding from this table is that cannabis is a much more powerful covariate for pancreatic cancer incidence than alcohol.
Table 4 sets out the findings for localized disease with again very similar findings.

3.4. Multivariable Regressions

A multivariable model with additive terms for the three substances of concern was tested during exploratory investigations. In each subgroup, only cannabis was left as the significant term in the final model. This is illustrated in Table 5 where the only significant terms in the table relate to cannabis exposure in all subgroups. As shown in this table, the levels of significance are high and the E-Values are also very elevated. This again shows that, statistically, cannabis is a more powerful covariate for pancreatic cancer incidence than alcohol and indeed cannabis terms would knock alcohol terms out of final regression models were model reduction to proceed by the classical method.
When a model which was interactive between the three substance exposure terms was tested during exploratory data mining, in each case the alcohol:cannabis interaction was the most significant and indeed only remaining term in the final model and excluded all other terms including daily cannabis use. For these reasons, this interaction was tested against PCI specifically across all combinations of race and sex in this 15–34-year cohort. As shown in Table 6, in every subgroup the alcohol:cannabis interaction was highly significant which in general terms was more marked amongst females. The lowest p-value shown in this table was for Caucasian American females (p = 2.50 × 10−7) and the highest E-Value pair was for Hispanic American females (E-Value estimate = 1.26 × 10102 and E-Value lower bound of 2.20 × 1074).
When localized disease was considered, similar results were obtained (Table 7).

4. Discussion

4.1. Main Results

As the most detailed data to date on the relationship between cannabis and pancreatic cancer, this study explores the sociodemographic profiling of younger pancreatic cancer patients in greater detail than has appeared previously. Of the many striking results described herein, arguably the most impressive is the clear direct linear relationship between daily cannabis exposure and pancreatic cancer incidence across all sex and ethnic cohorts in the 15–34-year age group as shown unequivocally in Figure 6 and Table 1, Table 2, Table 3, Table 4 and Table 5. This high performance is reflected in the unexpected supersession of cannabis over alcohol as a risk factor for pancreatic disease of any type and pancreatic malignant disease in particular. The predominant picture which emerges from this study is that the heightened and exponentially rising rate of PCI in young African American females appears to relate largely to their tripled environmental exposure to cannabinoids. The determination of any genetic contribution to this predisposition conferred by ethnicity must await further research.
The effect of increased cannabis use was so strong that in an additive multivariable model, cannabis knocked out the other terms for tobacco and alcohol in all sociodemographic strata in this age cohort. Given that alcohol is a long-established risk factor for pancreatic disease of any type, this is a most remarkable finding; albeit, it could be predicted from the shape of the curves in Figure 6. At multivariable interactive modelling, the cannabis:alcohol interaction was shown to be the only thing more powerful than the effect of the cannabis variable itself. This was also predominant across all sociodemographic strata clearly demonstrating the co-carcinogenicity of these two agents and the importance that consideration of this co-genotoxicity should take in public health planning in this field.
Importantly, Figure 5 shows clearly that the level of cannabis use in populations aged older than 50 years is much less than younger groups. This has two implications. Firstly, it implies that since cannabis is a significant pancreatic carcinogen this effect will be less marked in older populations, which is indeed what is observed because the rise in pancreatic cancer seems to be disproportionate amongst younger patients. Secondly, many studies show that both cannabinoid genotoxicity [25,104,105,106,107,108,109] and cannabinoid metabolic toxicity [110,111,112,113,114,115] follow exponential dose-response curves. We note again that metabolomotoxicity and epigenotoxicity are closely and inherently interrelated in numerous ways. It is also clear that the rise in the prevalence of cannabis use, the intensity of people smoking cannabis daily and the THC concentration of cannabis herb and resin and the many other cannabinoid products widely available on the market are all occurring at the same time. From a public health perspective, this can be expected to launch the community relatively abruptly into the high cannabis use zones where genotoxic outcomes become much more commonplace. It seems clear from the evidence that this is exactly what is occurring in younger populations with cannabis use prevalence and the intensity of use such as in younger African American females in the USA over the past decade.

4.2. Interpretation

Our interpretation of these data is that in this younger patient cohort, cannabis is clearly linked with pancreatic carcinogenesis across all sex and ethnic cohorts considered in this younger patients age cohort. Given that younger patients are normally considered to be at a low risk for pancreatic carcinogenesis, this is a signal finding indeed. The findings from this and similar studies are consistent with those from Europe [4,5,6,7].

4.3. Causal Inference

4.3.1. Qualitative Causal Inference

In 1965, after his experience with the tobacco-induced lung cancer epidemic of the 1950s, renowned epidemiologist A.B. Hill laid out nine criteria which were required to be fulfilled to confidently denote any particular epidemiological relationship as being causal in nature. These criteria were: strength of association, consistency amongst studies, specificity, temporality, coherence with known data, biological plausibility, a dose-response curve, similar analogies elsewhere and experimental confirmation. The current data clearly indicate that the cannabis–pancreatic cancer link fulfils the first eight of these criteria. Experimental confirmation is still lacking but only because laboratory studies have not yet been performed to the best of our knowledge at the time of writing. However, given the centrality of the biological pathways and mechanisms to supporting the causality of the cannabis–pancreatic cancerogenesis causal argument, discussion on this topic is expanded below.

4.3.2. Quantitative Causal Inference

One of the major limitations of observational studies is that they are subject to a failure to include the full universe of covariates in their analytical models. This error of uncontrolled confounding is known by several names. For an unknown covariate to disrupt a relationship which at first appears to be causal in nature, it must be correlated with both the exposure of interest and the outcome of concern. This co-correlation is quantified as the E-Value which is a mathematical re-formulation of the relative risk. As noted above, an E-Value greater than 1.25 generally denotes a causal relationship [100]. E-Values in excess of nine, such as those between lung cancer and tobacco smoking, are said to be high [116]. Studies of E-Values of 1.32 × 1055, 3.46 × 1053 and 2.04 × 1028 for Hispanic American, Caucasian American and African American females, respectively (Table 3), are far in excess of these levels, and are observed in causal relationships [100]. Moreover, the E-Value has a 95% lower bound which signifies the bottom end of the applicable confidence interval. In Table 3, the lower bounds for these E-Values are given as 1.20 × 1042, 8.41 × 1041 and 6.88 × 1019. These are also far in excess of the typical cut-offs required to show causality and provide robust arithmetical support to the main conclusions of the study. For these reasons, the study data provide epidemiological evidence of causality; albeit, this is not the same thing as experimental evidence of causality.

4.4. Mechanisms

4.4.1. Epigenomics

Cannabis has been found to have a large epigenomic footprint and changes the methylation status of up to 9% of the whole genome [15]. This has far-reaching implications on virtually every cellular system including the epigenomic machinery for DNA methylation and demethylation, for histone methylation and acetylation and their reversal and for the active repositioning of the nucleosomes to allow new gene transcription to occur [17,77,78,79,117]. Cannabis causes DNA breaks [24,25,26,27,28,29,30] which ages the genome [21]. Cannabis interferes with actin and tubulin and thus the cytoskeleton and machinery of the mitotic spindle [21] and the kinetochore–centrosomal machinery of chromosomal segregation at mitotic and meiotic division [17,118]. This interference with mitotic and meiotic division sets up micronucleus formation and chromosomal shattering (chromothripsis) [73] which is a powerful engine for major and widespread genomic rearrangement and instability [73,119,120,121,122,123,124,125,126,127,128,129,130] and also powerfully stimulates innate immunity through cyclic AMP-GMP (cGAS)–STimulator of Interferon Gamma (STING) signalling [43,130,131,132,133,134,135,136]. cGAS-STING has been shown to be the prime driver of the senescence-associated secretory phenotype (SASP) [135,137] which both drives innate immunity with the release of TNFα (tumour necrosis factor α), and interferon IL1β and interleukin 6 (IL6) [132]. Together, this drives senescence induction in nearby cells systemically [135,138] and also induces the inflammaging which is known to drive DNA damage and further exacerbate cellular aging [128,130,131,132,133,136,137,138,139,140,141,142,143,144,145,146,147]. The growth factors and cytokines and chemokines of the SASP are also known to drive malignant transformation [137,138,140,141]. cGAS-STING activation both increases the degree of malignancy and drives metastasis [148] which also provides an explanation for the many reports of increased incidence of cancer in young patients and their rapid demise due to a widespread and aggressive disease at presentation [149,150,151,152].
The action of cannabinoids is often directly immunostimulatory especially when acting via the type 1 cannabinoid receptor (CB1R) [90,153,154,155,156,157,158]. Cannabis disrupts both the tubulin code and the histone code [17,26,118,159,160,161,162,163]. By inhibiting mitochondrial metabolism [111,112,164,165,166,167,168,169,170], cannabinoids interfere with the supply of substrates to the epigenomic machinery and the supply of energy for both genomic and epigenomic maintenance. Due to mitochondrial inhibition, the cellular levels of lactate rises [171,172,173,174,175] which induces both passive and enzymically mediated lactylation of key metabolic enzymes [176,177]. This further reinforces the increase in cellular dependence on glycolysis and a stem-like cellular de-differentiated state and locks cells into a premalignant metabolomic–epigenomic state [176,177,178,179].
Low doses of cannabidiol have also been implicated directly in genotoxicity and have been shown to induce both DNA breaks and DNA base oxidation [25].
Since all of these changes are pro-carcinogenic, the link between cannabis and pancreatic cancer is not surprising from a mechanistic perspective.

4.4.2. Alcohol: Cannabis Interaction

As mentioned above, alcohol is a classical cause of acute pancreatitis and the mechanisms by which it predisposes to cancer has been shown several times to arise from the incomplete resolution of the inflammatory changes in the epigenome of the acinar stem cells [45,48]. Such epigenomic changes are an example of antagonistic pleiotropy in that they help the cells reactivate subsequent to another episode, but also reset the epigenome towards cancer and aging based on this alteration [45,48,180,181,182]. The fact that cannabis way outperformed the recognized pancreatic disease risk factor of alcohol in the study results presented above is one of the most remarkable findings to emerge from the present study.
As noted above, the effect of cannabinoids in tissues is often pro-inflammatory especially when acting through CB1R [90,153,154,155,156,157,158]. Cannabinoids interfere with the machinery of DNA methylation and the net effect of this is to cause DNA hypomethylation which is an aging change and also tends to weaken epigenomic cell linage specification and therefore the epigenomic hills between the “Waddington valleys” (to borrow Waddington’s now famous analogy [183]) and thus drive the cells both towards aging and de-differentiation [21]. Cannabinoids (including cannabidiol) also cause DNA breakages [24,25,26,27,28,29,30] which attract the large multiprotein epigenomic complexes to the sites of DNA breakage and away from their usual stations on the genome which directly lead to DNA demethylation and cellular aging [21]. By inducing chromosomal mis-segregation, cannabis drives chromosomal shattering releasing fragments of double-stranded DNA into the cytoplasm where it is a powerful stimulus for cytosolic sensors and a driver of both intracellular and extracellular innate immune effectors including cGAS-STING, some of the toll-like receptors (TLRs 3, 7 and 8), the RIG1 sensors, the AIM1 (Absent in Melanoma 1) and ASC1 (Apoptosis-associated Speck-like protein containing a CARD (caspase activation and recruitment domain) inflammasomes, the SASP [128,131,132,136,137,138,139,184,185,186] and a long list of cytosolic sensors of DNA including RNA polymerase III, DAI (DNA-dependent activator of interferons or Z-DNA binding protein), IFI16 (interferon gamma inducible protein 16), DDX4 (DEAD-Box helicase 4), LSm14A (LSM14A MRNA Processing Body Assembly Factor), LRRFIP1 (LRR Binding FLII Interacting Protein 1), Sox2 (SRY-Box Transcription Factor 2), DHX9/36 (DExH-Box Helicase 9/36) and Ku70 (X-Ray Repair Cross Complementing 6) [187].
Therefore, it is not at all surprising that there should be a powerful carcinogenic synergism between the two genotoxins, cannabis and alcohol, as documented above in our multivariate regression models both for PCI and for localized disease. However, it was surprising to us that this interactive effect was so powerful statistically as to lead to the complete deletion of cannabis from the final interactive models across all sociodemographic categories.

4.5. Future Research Directions

As this study is the most detailed investigation of the cannabis–pancreatic cancer link to have been published, it raises many further research questions both epidemiological and mechanistic in nature. Epidemiologically, with a combined sample of 734,761 cancers this dataset almost certainly lends itself to small area–space–time studies. Epidemiological studies would be facilitated by the availability from SAMHSA of NSDUH results in 405 small substate areas across the USA. Assuming that pancreatic cancer incidence could be similarly mapped at such fine spatial resolution across a series of aggregated years, this would lend itself to further space–time epidemiological models which is conceptually very powerful. Ideally such space–time modelling could also incorporate inverse probability weighting to allow causal geospatial models to be developed, although to our knowledge this methodology has not yet been invented.
From a mechanistic perspective, many studies in the basic sciences should be undertaken to explore cannabinoid-induced pancreatic tumourigenesis in molecular detail. The above-described convergence of pro-inflammatory, pro-aging, de-differentiating, direct genotoxic and various epigenotoxic pathways should all be dissected to define the route to carcinogenesis induced by cannabis. It needs to be resolved how many cannabinoids are implicated in pancreaticocarcinogenesis especially in the light of earlier demonstrations that the genotoxic moiety resides with the olivetol nucleus on the C-ring which is common to all cannabinoids [26]. It is common for the profusion of names for modern cannabinoids such as Δ8-, Δ0-, Δ10, Δ11-tetrahydrocannabinol etc. to confuse lay people. However, as these designations only refer to the carbon atom on the A ring which is hydroxylated and the genotoxic moiety is the olevitol nucleus of the C-ring, none of these modifications have been found to materially impact the genotoxic actions after extensive in vitro testing [26] and nor would they really be expected to do so (because it is the wrong benzene ring). Dose-response and timing studies need to be performed for cannabinoids in clinically relevant concentrations. It may be that as such processes are understood, agents acting either in the endocannabinoid system or perhaps in other pathways may be developed which will become lead candidates for further therapeutic development for this very serious disease.
Cannabinoids are known to negatively impact the mitochondria in many ways [111,112,164,165,166,167,168,169,170]. This will induce genotoxic stress through mitonuclear signalling and also impede the genomic and epigenomic machinery which are both dependent on mitochondrially derived substrates and energy and homeostatic mechanisms to ensure genomic stability. Therefore, the cannabinoid-induced mitochondriopathy can reasonably be expected to interact powerfully and negatively with the cannabinoid-induced epigenotoxicity. Since as described the lactylation state is modified by cannabinoids and is involved in numerous cancerogenic pathways, it is likely that the impact of cannabinoids on the lactylation machinery and metabolome presents a key investigative window which will likely find downstream utility in areas such as stem cell and regenerative medicine and cancer biology. These mechanisms need to be investigated both in the mitochondria and in the epigenome and genome in a coordinated and comprehensive manner.
As described above, cannabinoids generally have a heavy epigenetic footprint. The general footprint of cannabinoids is DNA hypomethylation which is a classical age-related change. Thus, phytocannabinoids can be expected to relatively de-differentiate cells. This putative cannabinoid-induced de-differentiation needs to be compared with premalignant de-differentiation and also that occurring normally in inflammatory and post-inflammatory states to investigate how this might compare to these better characterized states. In addition, the epigenome is known to be further impacted by immune and metabolic processes, all of which are disrupted by cannabinoid exposure.
The deep learning algorithms of artificial intelligence may find a place in the epidemiological and other studies of pancreatic cancer in years to come. Some of these directions have been explored in a recent paper on the application of integrative medicine to the pathobiology of cancer of the pancreas [188].
Thus, the present study strongly suggests many lines for further research investigations both in epidemiology and in the basic sciences.

4.6. Generalizability

Given that the sample on which this study is based is the whole US population of pancreatic cancer from 2001 to 2018 combined with SAMHSA’s highly reputable nationally representative annual NSDUH study including its high response rate, and given also that the present results are closely concordant with other published reports and with recent analyses from Europe, we are confident that these results are indeed generally applicable [4,5,6,7]. This view is reinforced by the very high E-Values reported throughout which exclude significant effects of unmeasured confounder covariates. The E-Values reported in the present study are far in excess of the threshold for causality which is usually said to be 1.25 [100] and also exceed the E-Value of the relationship of tobacco use with lung cancer incidence which is said to be nine and is typically regarded as being in the high range [116]. For these reasons, it is clear that the relationship between cannabis use and pancreatic cancer is much stronger than the tobacco–lung cancer relationship which is widely acknowledged to be causal. The present findings are also supported by the robust observation that the effects of cannabis interact with, and are exacerbated by, those of alcohol which is a known carcinogen in the pancreas [45,46,48].

4.7. Strengths and Limitations

This study has a number of strengths. Cancer incidence data were based on a national sample across the USA from 2001 to 2018. Drug use data were taken from the nationally representative survey with a high response rate which was SAMHSA’s NSDUH. The analysis was stratified separately for each sex and major ethnic group with the single 15–34-year age group. Highly concordant and consistent results were obtained from bivariate, additive multivariable and interactive multivariable models. Furthermore, the findings for cannabis were noted to positively interact with alcohol which is a known pancreatic carcinogen. For these reasons, we regard the present analysis as being powerful in its elegant analytical simplicity. Importantly, the high incidence rate of pancreatic cancer in young African American females was found to be likely related mainly to their lifestyle exposures. Study shortcomings are that, in common with many other epidemiological studies, individual participant data were not available to the present investigators. Stratified data on income were not available to the present research team. Controlling for these and other covariates remains for future research projects. However, in view of the very high E-Values reported herein, its consistency with other studies both in the USA and Europe, its concordance with a variety of basic aetiological mechanisms from laboratory and cellular science, and of similar studies which do take income into account [4,5,6], and its relationship with the known carcinogen ethanol, we are confident that controlling for these external covariates will not materially impact the major study conclusions.

5. Conclusions

The main results from this study are to confirm the opening hypotheses that cannabis is involved causally in pancreatic carcinogenesis in all sex and ethnic cohorts in the 15–35-year age bracket and thus it appears to be at least one of the previously unknown principal driver/s of the present renaissance in PCI across the USA. The present results are supported by other studies from both the USA and Europe and by much experimental research [4,5,6,7]. Cannabis was observed to powerfully interact at the epidemiological level with alcohol exposure and clearly this cooperative carcinogenicity requires further basic science and epidemiological research. Importantly, this report adds cancer of the pancreas to the growing list of cancers which cannabis has been shown to be driving including testicular, liver and breast cancer, childhood leukaemia and total paediatric cancer. Clearly, further research is indicated of the pathways to cannabinoid-induced tumourigenesis particularly in the metabolomic, immunomic and epigenomic fields. However, in the light of now-robust epidemiological evidence implicating cannabis in yet another significant cancer renaissance nationwide across the USA [10,83,84,85,86,87,88,89,90], it becomes essential and imperative that communities take cannabinoid carcinogenesis seriously and restrict access of their populations to cannabinoid genotoxicity just as it is for many other significantly genotoxic compounds. This is not only for cancer prevention and public health concerns in regard to extant populations, but also to protect the generations to come from the implicit and now well-defined multigenerational cannabinoid epigenotoxicity with which carcinogenicity is so often intricately and intimately involved [77,78,79,117,189,190,191,192,193].

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/gastroent14020016/s1. Supplementary Figures S1–S3 have been made available online together with this paper on the journal website.

Author Contributions

A.S.R. assembled the data, designed and conducted the analyses, and wrote the first manuscript draft. G.K.H. provided technical and logistic support, co-wrote the paper, assisted with gaining ethical approval, provided advice on manuscript preparation and general guidance to study conduct. A.S.R. had the idea for the article, performed the literature search, wrote the first draft and is the guarantor for the article. All authors have read and agreed to the published version of the manuscript.

Funding

No funding was provided for this study. No funding organization played any role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.

Institutional Review Board Statement

The Human Research Ethics Committee of the University of Western Australia provided ethical approval for the study to be undertaken 24 September 2021 (No. RA/4/20/4724).

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated or analysed during this study are included in this published article and its supplementary information files. Data along with the relevant R code have been made publicly available on the Mendeley Database Repository and can be accessed from these URLs: Raw Data: https:/data.mendeley.com/datasets/wjb5f622n4 (accessed on 19 February 2023) https://doi.org/10.17632/wjb5f622n4.1. Pre-processed Data: https:/data.mendeley.com/datasets/4s87k4v5n7 (accessed on 19 February 2023) https://doi.org/10.17632/4s87k4v5n7.1. Major datasets: https:/data.mendeley.com/datasets/y32txvjhpf (accessed on 19 February 2023) https://doi.org/10.17632/y32txvjhpf.1. Computational Code: https:/data.mendeley.com/datasets/shcnmv5dhx (accessed on 19 February 2023) https://doi.org/10.17632/shcnmv5dhx.1.

Acknowledgments

All authors had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis.

Conflicts of Interest

The authors declare that they have no competing interest.

References

  1. SEER Explorer. Available online: https://seer.cancer.gov/explorer/application.html (accessed on 19 February 2023).
  2. Ugai, T.; Sasamoto, N.; Lee, H.Y.; Ando, M.; Song, M.; Tamimi, R.M.; Kawachi, I.; Campbell, P.T.; Giovannucci, E.L.; Weiderpass, E.; et al. Is early-onset cancer an emerging global epidemic? Current evidence and future implications. Nat. Rev. Clin. Oncol. 2022, 19, 656–673. [Google Scholar] [CrossRef] [PubMed]
  3. Abboud, Y.; Samaan, J.S.; Oh, J.; Jiang, Y.; Randhawa, N.; Lew, D.; Ghaith, J.; Pala, P.; Leyson, C.; Watson, R.; et al. Increasing Pancreatic Cancer Incidence in Young Women in the US: A Population-Based Time-Trend Analysis, 2001–2018. Gastroenterology, 2023; in press. [Google Scholar]
  4. Reece, A.S.; Hulse, G.K. 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. Arch. Public. Health 2022, 80, 99–133. [Google Scholar] [CrossRef] [PubMed]
  5. Reece, A.S.; Hulse, G.K. 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. Arch. Public. Health 2022, 80, 100–135. [Google Scholar] [CrossRef] [PubMed]
  6. Reece, A.S.; Hulse, G.K. 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. Arch. Public Health 2022, 80, 100–136. [Google Scholar] [CrossRef] [PubMed]
  7. Reece, A.S.; Hulse, G.K. Epidemiological Overview of Multidimensional Chromosomal and Genome Toxicity of Cannabis Exposure in Congenital Anomalies and Cancer Development. Sci. Rep. 2021, 11, 13892–13912. [Google Scholar] [CrossRef]
  8. Key Substance Use and Mental Health Indicators in the United States: Results from the 2020 National Survey on Drug Use and Health (NSDUH). Available online: https://www.samhsa.gov/data/sites/default/files/reports/rpt35325/NSDUHFFRPDFWHTMLFiles2020/2020NSDUHFFR1PDFW102121.pdf (accessed on 20 January 2023).
  9. United National Office of Drugs and Crime. World Drug Report 2019; World Health Organization Office of Drugs and Crime, Ed.; United National World Health Organization: Geneva, Switzerland, 2019; Volume 1–5, Available online: https://wdr.unodc.org/wdr2019/index.html (accessed on 20 January 2023).
  10. Patsenker, E.; Stickel, F. Cannabinoids in liver diseases. Clin. Liver Dis. 2016, 7, 21–25. [Google Scholar] [CrossRef]
  11. DiNieri, J.A.; Wang, X.; Szutorisz, H.; Spano, S.M.; Kaur, J.; Casaccia, P.; Dow-Edwards, D.; Hurd, Y.L. Maternal cannabis use alters ventral striatal dopamine D2 gene regulation in the offspring. Biol. Psychiatry 2011, 70, 763–769. [Google Scholar] [CrossRef]
  12. Ellis, R.J.; Bara, A.; Vargas, C.A.; Frick, A.L.; Loh, E.; Landry, J.; Uzamere, T.O.; Callens, J.E.; Martin, Q.; Rajarajan, P.; et al. Prenatal Δ(9)-Tetrahydrocannabinol Exposure in Males Leads to Motivational Disturbances Related to Striatal Epigenetic Dysregulation. Biol. Psychiatry 2021, 92, 127–138. [Google Scholar] [CrossRef]
  13. Szutorisz, H.; Hurd, Y.L. Epigenetic Effects of Cannabis Exposure. Biol. Psychiatry 2016, 79, 586–594. [Google Scholar] [CrossRef]
  14. Szutorisz, H.; DiNieri, J.A.; Sweet, E.; Egervari, G.; Michaelides, M.; Carter, J.M.; Ren, Y.; Miller, M.L.; Blitzer, R.D.; Hurd, Y.L. Parental THC exposure leads to compulsive heroin-seeking and altered striatal synaptic plasticity in the subsequent generation. Neuropsychopharmacology 2014, 39, 1315–1323. [Google Scholar] [CrossRef] [PubMed]
  15. Watson, C.T.; Szutorisz, H.; Garg, P.; Martin, Q.; Landry, J.A.; Sharp, A.J.; Hurd, Y.L. Genome-Wide DNA Methylation Profiling Reveals Epigenetic Changes in the Rat Nucleus Accumbens Associated with Cross-Generational Effects of Adolescent THC Exposure. Neuropsychopharmacology 2015, 40, 2993–3005. [Google Scholar] [CrossRef] [PubMed]
  16. Szutorisz, H.; Hurd, Y.L. High times for cannabis: Epigenetic imprint and its legacy on brain and behavior. Neurosci. Biobehav. Rev. 2018, 85, 93–101. [Google Scholar] [CrossRef] [PubMed]
  17. Schrott, R.; Murphy, S.K.; Modliszewski, J.L.; King, D.E.; Hill, B.; Itchon-Ramos, N.; Raburn, D.; Price, T.; Levin, E.D.; Vandrey, R.; et al. Refraining from use diminishes cannabis-associated epigenetic changes in human sperm. Environ. Epigenetics 2021, 7, dvab009. [Google Scholar] [CrossRef] [PubMed]
  18. Schrott, R.; Greeson, K.W.; King, D.; Symosko Crow, K.M.; Easley, C.A.T.; Murphy, S.K. Cannabis alters DNA methylation at maternally imprinted and autism candidate genes in spermatogenic cells. Syst. Biol. Reprod. Med. 2022, 68, 357–369. [Google Scholar] [CrossRef]
  19. Zykovich, A.; Hubbard, A.; Flynn, J.M.; Tarnopolsky, M.; Fraga, M.F.; Kerksick, C.; Ogborn, D.; MacNeil, L.; Mooney, S.D.; Melov, S. Genome-wide DNA methylation changes with age in disease-free human skeletal muscle. Aging Cell 2014, 13, 360–366. [Google Scholar] [CrossRef]
  20. Greenblatt, M.B.; Debnath, S. A stem-cell basis for skeletal ageing. Nature 2021, 597, 182–183. [Google Scholar] [CrossRef]
  21. Yang, J.H.; Hayano, M.; Griffin, P.T.; Amorim, J.A.; Bonkowski, M.S.; Apostolides, J.K.; Salfati, E.L.; Blanchette, M.; Munding, E.M.; Bhakta, M.; et al. Loss of epigenetic information as a cause of mammalian aging. Cell 2023, 186, 305–326.E27. [Google Scholar] [CrossRef]
  22. Childs, B.; van Deursen, J. Inhibition of ‘jumping genes’ promotes healthy ageing. Nature 2019, 566, 46–48. [Google Scholar] [CrossRef]
  23. Hoyt, S.J.; Storer, J.M.; Hartley, G.A.; Grady, P.G.S.; Gershman, A.; de Lima, L.G.; Limouse, C.; Halabian, R.; Wojenski, L.; Rodriguez, M.; et al. From telomere to telomere: The transcriptional and epigenetic state of human repeat elements. Science 2022, 376, eabk3112. [Google Scholar] [CrossRef]
  24. Stenchever, M.A.; Kunysz, T.J.; Allen, M.A. Chromosome breakage in users of marihuana. Am. J. Obs. Gynecol. 1974, 118, 106–113. [Google Scholar] [CrossRef] [PubMed]
  25. Russo, C.; Ferk, F.; Mišík, M.; Ropek, N.; Nersesyan, A.; Mejri, D.; Holzmann, K.; Lavorgna, M.; Isidori, M.; Knasmüller, S. Low doses of widely consumed cannabinoids (cannabidiol and cannabidivarin) cause DNA damage and chromosomal aberrations in human-derived cells. Arch. Toxicol. 2019, 93, 179–188. [Google Scholar] [CrossRef] [PubMed]
  26. Nahas, G.G.; Morishima, A.; Desoize, B. Effects of cannabinoids on macromolecular synthesis and replication of cultured lymphocytes. Fed. Proc. 1977, 36, 1748–1752. [Google Scholar] [PubMed]
  27. Dalterio, S.L.; de Rooij, D.G. Maternal cannabinoid exposure. Effects on spermatogenesis in male offspring. Int. J. Androl. 1986, 9, 250–258. [Google Scholar] [CrossRef] [PubMed]
  28. Leuchtenberger, C.; Leuchtenberger, R. Morphological and cytochemical effects of marijuana cigarette smoke on epithelioid cells of lung explants from mice. Nature 1971, 234, 227–229. [Google Scholar] [CrossRef]
  29. Leuchtenberger, C.; Leuchtenberger, R.; Schneider, A. Effects of marijuana and tobacco smoke on human lung physiology. Nature 1973, 241, 137–139. [Google Scholar] [CrossRef]
  30. Dalterio, S.; Badr, F.; Bartke, A.; Mayfield, D. Cannabinoids in male mice: Effects on fertility and spermatogenesis. Science 1982, 216, 315–316. [Google Scholar] [CrossRef]
  31. Hnisz, D.; Abraham, B.J.; Lee, T.I.; Lau, A.; Saint-Andre, V.; Sigova, A.A.; Hoke, H.A.; Young, R.A. Super-enhancers in the control of cell identity and disease. Cell 2013, 155, 934–947. [Google Scholar] [CrossRef]
  32. Hnisz, D.; Schuijers, J.; Lin, C.Y.; Weintraub, A.S.; Abraham, B.J.; Lee, T.I.; Bradner, J.E.; Young, R.A. Convergence of developmental and oncogenic signaling pathways at transcriptional super-enhancers. Mol. Cell 2015, 58, 362–370. [Google Scholar] [CrossRef]
  33. Kloetgen, A.; Thandapani, P.; Ntziachristos, P.; Ghebrechristos, Y.; Nomikou, S.; Lazaris, C.; Chen, X.; Hu, H.; Bakogianni, S.; Wang, J.; et al. Three-dimensional chromatin landscapes in T cell acute lymphoblastic leukemia. Nat. Genet. 2020, 52, 388–400. [Google Scholar] [CrossRef]
  34. Yang, H.; Zhang, H.; Luan, Y.; Liu, T.; Yang, W.; Roberts, K.G.; Qian, M.X.; Zhang, B.; Yang, W.; Perez-Andreu, V.; et al. Noncoding genetic variation in GATA3 increases acute lymphoblastic leukemia risk through local and global changes in chromatin conformation. Nat. Genet. 2022, 54, 170–179. [Google Scholar] [CrossRef] [PubMed]
  35. Xu, J.; Song, F.; Lyu, H.; Kobayashi, M.; Zhang, B.; Zhao, Z.; Hou, Y.; Wang, X.; Luan, Y.; Jia, B.; et al. Subtype-specific 3D genome alteration in acute myeloid leukaemia. Nature 2022, 611, 387–398. [Google Scholar] [CrossRef] [PubMed]
  36. Ishoey, M.; Chorn, S.; Singh, N.; Jaeger, M.G.; Brand, M.; Paulk, J.; Bauer, S.; Erb, M.A.; Parapatics, K.; Muller, A.C.; et al. Translation Termination Factor GSPT1 Is a Phenotypically Relevant Off-Target of Heterobifunctional Phthalimide Degraders. ACS Chem. Biol. 2018, 13, 553–560. [Google Scholar] [CrossRef] [PubMed]
  37. Duan, Q.; Mao, X.; Xiao, Y.; Liu, Z.; Wang, Y.; Zhou, H.; Zhou, Z.; Cai, J.; Xia, K.; Zhu, Q.; et al. Super enhancers at the miR-146a and miR-155 genes contribute to self-regulation of inflammation. Biochim. Biophys. Acta 2016, 1859, 564–571. [Google Scholar] [CrossRef]
  38. Ounzain, S.; Pedrazzini, T. Super-enhancer lncs to cardiovascular development and disease. Biochim. Biophys. Acta 2016, 1863, 1953–1960. [Google Scholar] [CrossRef] [PubMed]
  39. Kohnken, R.; Wen, J.; Mundy-Bosse, B.; McConnell, K.; Keiter, A.; Grinshpun, L.; Hartlage, A.; Yano, M.; McNeil, B.; Chakravarti, N.; et al. Diminished microRNA-29b level is associated with BRD4-mediated activation of oncogenes in cutaneous T-cell lymphoma. Blood 2018, 131, 771–781. [Google Scholar] [CrossRef] [PubMed]
  40. Ott, C.J.; Federation, A.J.; Schwartz, L.S.; Kasar, S.; Klitgaard, J.L.; Lenci, R.; Li, Q.; Lawlor, M.; Fernandes, S.M.; Souza, A.; et al. Enhancer Architecture and Essential Core Regulatory Circuitry of Chronic Lymphocytic Leukemia. Cancer Cell 2018, 34, 982–995.E7. [Google Scholar] [CrossRef]
  41. Yokoyama, Y.; Zhu, H.; Lee, J.H.; Kossenkov, A.V.; Wu, S.Y.; Wickramasinghe, J.M.; Yin, X.; Palozola, K.C.; Gardini, A.; Showe, L.C.; et al. BET Inhibitors Suppress ALDH Activity by Targeting ALDH1A1 Super-Enhancer in Ovarian Cancer. Cancer Res. 2016, 76, 6320–6330. [Google Scholar] [CrossRef]
  42. Whyte, W.A.; Orlando, D.A.; Hnisz, D.; Abraham, B.J.; Lin, C.Y.; Kagey, M.H.; Rahl, P.B.; Lee, T.I.; Young, R.A. Master transcription factors and mediator establish super-enhancers at key cell identity genes. Cell 2013, 153, 307–319. [Google Scholar] [CrossRef]
  43. Skoulidis, F.; Heymach, J.V.; Cascone, T. A STING operation to expose KRAS and STK11 co-mutated lung cancers. Cancer Cell 2022, 40, 1073–1076. [Google Scholar] [CrossRef]
  44. Kortlever, R.M.; Sodir, N.M.; Wilson, C.H.; Burkhart, D.L.; Pellegrinet, L.; Brown Swigart, L.; Littlewood, T.D.; Evan, G.I. Myc Cooperates with Ras by Programming Inflammation and Immune Suppression. Cell 2017, 171, 1301–1315.E14. [Google Scholar] [CrossRef] [PubMed]
  45. Neuhöfer, P.; Roake, C.M.; Kim, S.J.; Lu, R.J.; West, R.B.; Charville, G.W.; Artandi, S.E. Acinar cell clonal expansion in pancreas homeostasis and carcinogenesis. Nature 2021, 597, 715–719. [Google Scholar] [CrossRef] [PubMed]
  46. Park, W.; Chawla, A.; O’Reilly, E.M. Pancreatic Cancer: A Review. Jama 2021, 326, 851–862. [Google Scholar] [CrossRef] [PubMed]
  47. Lan, L.; Evan, T.; Li, H.; Hussain, A.; Ruiz, E.J.; Zaw Thin, M.; Ferreira, R.M.M.; Ps, H.; Riising, E.M.; Zen, Y.; et al. GREM1 is required to maintain cellular heterogeneity in pancreatic cancer. Nature 2022, 607, 163–168. [Google Scholar] [CrossRef]
  48. Del Poggetto, E.; Ho, I.L.; Balestrieri, C.; Yen, E.Y.; Zhang, S.; Citron, F.; Shah, R.; Corti, D.; Diaferia, G.R.; Li, C.Y.; et al. Epithelial memory of inflammation limits tissue damage while promoting pancreatic tumorigenesis. Science 2021, 373, eabj0486. [Google Scholar] [CrossRef]
  49. Klein, T.W.; Newton, C.; Larsen, K.; Lu, L.; Perkins, I.; Nong, L.; Friedman, H. The cannabinoid system and immune modulation. J. Leukoc. Biol. 2003, 74, 486–496. [Google Scholar] [CrossRef]
  50. Bindukumar, B.; Mahajan, S.D.; Reynolds, J.L.; Hu, Z.; Sykes, D.E.; Aalinkeel, R.; Schwartz, S.A. Genomic and proteomic analysis of the effects of cannabinoids on normal human astrocytes. Brain Res. 2008, 1191, 1–11. [Google Scholar] [CrossRef]
  51. Eisenstein, T.K.; Meissler, J.J. Effects of Cannabinoids on T-cell Function and Resistance to Infection. J. Neuroimmune Pharm. 2015, 10, 204–216. [Google Scholar] [CrossRef]
  52. Zumbrun, E.E.; Sido, J.M.; Nagarkatti, P.S.; Nagarkatti, M. Epigenetic Regulation of Immunological Alterations Following Prenatal Exposure to Marijuana Cannabinoids and its Long Term Consequences in Offspring. J. Neuroimmune Pharm. 2015, 10, 245–254. [Google Scholar] [CrossRef]
  53. Chiurchiu, V. Endocannabinoids and Immunity. Cannabis Cannabinoid Res. 2016, 1, 59–66. [Google Scholar] [CrossRef]
  54. Gallily, R.; Yekhtin, Z. Avidekel Cannabis extracts and cannabidiol are as efficient as Copaxone in suppressing EAE in SJL/J mice. Inflammopharmacology 2018, 27, 167–173. [Google Scholar] [CrossRef] [PubMed]
  55. Kaplan, B.L.F. Evaluation of Marijuana Compounds on Neuroimmune Endpoints in Experimental Autoimmune Encephalomyelitis. Curr. Protoc. Toxicol. 2018, 75, 11 25 11–11 25 22. [Google Scholar] [CrossRef] [PubMed]
  56. Allen, J.P.; Danoff, J.S.; Costello, M.A.; Hunt, G.L.; Hellwig, A.F.; Krol, K.M.; Gregory, S.G.; Giamberardino, S.N.; Sugden, K.; Connelly, J.J. Lifetime marijuana use and epigenetic age acceleration: A 17-year prospective examination. Drug Alcohol Depend. 2022, 233, 109363. [Google Scholar] [CrossRef] [PubMed]
  57. Phillips, K.T.; Pedula, K.L.; Choi, N.G.; Tawara, K.K.; Simiola, V.; Satre, D.D.; Owen-Smith, A.; Lynch, F.F.; Dickerson, J. Chronic health conditions, acute health events, and healthcare utilization among adults over age 50 in Hawai’i who use cannabis: A matched cohort study. Drug Alcohol Depend. 2022, 234, 109387. [Google Scholar] [CrossRef] [PubMed]
  58. Reece, A.S.; Norman, A.; Hulse, G.K. Cannabis exposure as an interactive cardiovascular risk factor and accelerant of organismal ageing: A longitudinal study. BMJ Open 2016, 6, e011891–e011901. [Google Scholar] [CrossRef] [PubMed]
  59. Forrester, M.B.; Merz, R.D. Risk of selected birth defects with prenatal illicit drug use, Hawaii, 1986–2002. J. Toxicol. Environ. Health 2007, 70, 7–18. [Google Scholar] [CrossRef]
  60. Reece, A.S.; Hulse, G.K. Cannabinoid- and Substance- Relationships of European Congenital Anomaly Patterns: A Space-Time Panel Regression and Causal Inferential Study. Environ. Epigenetics 2022, 8, dvab015. [Google Scholar] [CrossRef]
  61. Reece, A.S.; Hulse, G.K. Cannabinoid Genotoxicity and Congenital Anomalies: A Convergent Synthesis of European and USA Datasets. In Cannabis, Cannabinoids and Endocannabinoids; Preedy, V., Patel, V., Eds.; Elsevier: London, UK, 2022; Volume 1, in press. [Google Scholar]
  62. Reece, A.S.; Hulse, G.K. Congenital Anomaly Epidemiological Correlates of Δ8THC Across USA 2003–2016: Panel Regression and Causal Inferential Study. Environ. Epigenetics 2022, 8, dvac012. [Google Scholar] [CrossRef]
  63. Reece, A.S.; Hulse, G.K. European Epidemiological Patterns of Cannabis- and Substance- Related Congenital Body Wall Anomalies: Geospatiotemporal and Causal Inferential Study. Int. J. Environ. Res. Public Health 2022, 19, 9027. [Google Scholar] [CrossRef]
  64. Reece, A.S.; Hulse, G.K. Cannabis- and Substance- Related Epidemiological Patterns of Chromosomal Congenital Anomalies in Europe: Geospatiotemporal and Causal Inferential Study. Int. J. Environ. Res. Public Health 2022, 19, 11208. [Google Scholar] [CrossRef]
  65. Reece, A.S.; Hulse, G.K. European Epidemiological Patterns of Cannabis- and Substance- Related Congenital Cardiovascular Anomalies: Geospatiotemporal and Causal Inferential Study. Environ. Epigenetics 2022, 8, dvac015. [Google Scholar] [CrossRef]
  66. Reece, A.S.; Hulse, G.K. European Epidemiological Patterns of Cannabis- and Substance- Related Congenital Neurological Anomalies: Geospatiotemporal and Causal Inferential Study. Int. J. Environ. Res. Public Health 2022, 20, 441. [Google Scholar] [CrossRef]
  67. Reece, A.S.; Hulse, G.K. Effects of Cannabis on Congenital Limb Anomalies in 14 European Nations: A Geospatiotemporal and Causal Inferential Study. Environ. Epigenetics 2022, 8, dvac016. [Google Scholar] [CrossRef]
  68. Reece, A.S.; Hulse, G.K. European Epidemiological Patterns of Cannabis- and Substance- Related Congenital Uronephrological Anomalies: Geospatiotemporal and Causal Inferential Study. Int. J. Environ. Res. Public Health 2022, 19, 13769. [Google Scholar] [CrossRef]
  69. Reece, A.S.; Hulse, G.K. Geospatiotemporal and Causal Inferential Study of European Epidemiological Patterns of Cannabis- and Substance- Related Congenital Orofacial Anomalies. J. Xenobiotics 2023, 13, 42–74. [Google Scholar] [CrossRef]
  70. Reece, A.S.; Hulse, G.K. Congenital Gastrointestinal Anomalies in Europe 2010–2019: A Geospatiotemporal and Causal Inferential Study of Epidemiological Patterns in Relationship to Cannabis- and Substance- Exposure. Gastrointest. Insights 2023, 14, 64–109. [Google Scholar] [CrossRef]
  71. Reece, A.S.; Hulse, G.K. Patterns of Cannabis- and Substance- Related Congenital General Anomalies in Europe: A Geospatiotemporal and Causal Inferential Study. Pediatr. Rep. 2023, 15, 69–121. [Google Scholar] [CrossRef]
  72. Reece, A.S.; Hulse, G.K. Chapter 3: Geospatiotemporal and Causal Inferential Analysis of United States Congenital Anomalies as a Function of Multiple Cannabinoid- and Substance- Exposures: Phenocopying Thalidomide and Hundred Megabase-Scale Genotoxicity. In Epidemiology of Cannabis: Genotoxicity and Neurotoxicity, Epigenomics and Aging; Elsevier: New York, NY, USA, 2023; Volume 1, in press. [Google Scholar]
  73. Reece, A.S.; Hulse, G.K. Chromothripsis and epigenomics complete causality criteria for cannabis- and addiction-connected carcinogenicity, congenital toxicity and heritable genotoxicity. Mutat. Res. 2016, 789, 15–25. [Google Scholar] [CrossRef]
  74. Reece, A.S.; Hulse, G.K. Cannabis Teratology Explains Current Patterns of Coloradan Congenital Defects: The Contribution of Increased Cannabinoid Exposure to Rising Teratological Trends. Clin. Pediatr. 2019, 58, 1085–1123. [Google Scholar] [CrossRef]
  75. Reece, A.S.; Hulse, G.K. Canadian Cannabis Consumption and Patterns of Congenital Anomalies: An Ecological Geospatial Analysis. J. Addict. Med. 2020, 14, e195–e210. [Google Scholar] [CrossRef] [PubMed]
  76. Reece, A.S.; Hulse, G.K. Geotemporospatial and causal inference epidemiological analysis of US survey and overview of cannabis, cannabidiol and cannabinoid genotoxicity in relation to congenital anomalies 2001–2015. BMC Pediatr. 2022, 22, 47–124. [Google Scholar] [CrossRef] [PubMed]
  77. Reece, A.S.; Hulse, G.K. Epigenomic and Other Evidence for Cannabis-Induced Aging Contextualized in a Synthetic Epidemiologic Overview of Cannabinoid-Related Teratogenesis and Cannabinoid-Related Carcinogenesis. Int. J. Environ. Res. Public Health 2022, 19, 16721. [Google Scholar] [CrossRef]
  78. Reece, A.S.; Hulse, G.K. Clinical Epigenomic Explanation of the Epidemiology of Cannabinoid Genotoxicity Manifesting as Transgenerational Teratogenesis, Cancerogenesis and Aging Acceleration. Int. J. Environ. Res. Public Health 2023, 20, 3360. [Google Scholar] [CrossRef]
  79. Reece, A.S.; Hulse, G.K. Chapter 5: Multivalent Cannabinoid Epigenotoxicities and Multigenerational Aging. In Epidemiology of Cannabis: Genotoxicity, Neurotoxicity, Epigenomics and Aging; Elsevier: New York, NY, USA, 2023; Volume 1, in press. [Google Scholar]
  80. Reece, A.S. Disruption of Interlocking Synchrony between Metabolome and Epigenome Key to Understanding Widespread Embyrotoxicity and Carcinogenicity of Diverse Cannabinoids. Available online: https://www.bmj.com/content/377/bmj.o1567/rr-0 (accessed on 20 January 2023).
  81. Reece, A.S.; Hulse, G.K. Chapter 4: Geospatiotemporal and Causal Inferential Epidemiological Survey and Exploration of Cannabinoid- and Substance- Related Carcinogenesis in USA 2003–2017. In Epidemiology of Cannabis: Genotoxicity and Neurotoxicity, Epigenomics and Aging; Elsevier: New York, NY, USA, 2023; Volume 1, in press: 2500. [Google Scholar]
  82. Reece, A.S.; Hulse, G.K. Epidemiology of Δ8THC-Related Carcinogenesis in USA: A Panel Regression and Causal Inferential Study. Int. J. Environ. Res. Public Health 2022, 19, 7726. [Google Scholar] [CrossRef]
  83. Reece, A.S.; Hulse, G.K. State Trends of Cannabis Liberalization as a Causal Driver of Increasing Testicular Cancer Rates across the USA. Int. J. Environ. Res. Public Health 2022, 19, 12759. [Google Scholar] [CrossRef]
  84. Reece, A.S.; Hulse, G.K. Causal inference multiple imputation investigation of the impact of cannabinoids and other substances on ethnic differentials in US testicular cancer incidence. BMC Pharm. Toxicol. 2021, 22, 40–71. [Google Scholar] [CrossRef]
  85. Reece, A.S.; Hulse, G.K. Geospatiotemporal and Causal Inference Study of Cannabis and Other Drugs as Risk Factors for Female Breast Cancer USA 2003–2017. Environ. Epigenetics 2022, 8, dvac006. [Google Scholar] [CrossRef]
  86. Reece, A.S.; Hulse, G.K. Cannabinoid exposure as a major driver of pediatric acute lymphoid Leukaemia rates across the USA: Combined geospatial, multiple imputation and causal inference study. BMC Cancer 2021, 21, 984–1017. [Google Scholar] [CrossRef]
  87. Patsenker, E.; Stoll, M.; Millonig, G.; Agaimy, A.; Wissniowski, T.; Schneider, V.; Mueller, S.; Brenneisen, R.; Seitz, H.K.; Ocker, M.; et al. Cannabinoid receptor type I modulates alcohol-induced liver fibrosis. Mol. Med. 2011, 17, 1285–1294. [Google Scholar] [CrossRef]
  88. Mukhopadhyay, B.; Cinar, R.; Yin, S.; Liu, J.; Tam, J.; Godlewski, G.; Harvey-White, J.; Mordi, I.; Cravatt, B.F.; Lotersztajn, S.; et al. Hyperactivation of anandamide synthesis and regulation of cell-cycle progression via cannabinoid type 1 (CB1) receptors in the regenerating liver. Proc. Natl. Acad. Sci. USA 2011, 108, 6323–6328. [Google Scholar] [CrossRef] [PubMed]
  89. Mukhopadhyay, B.; Liu, J.; Osei-Hyiaman, D.; Godlewski, G.; Mukhopadhyay, P.; Wang, L.; Jeong, W.I.; Gao, B.; Duester, G.; Mackie, K.; et al. Transcriptional regulation of cannabinoid receptor-1 expression in the liver by retinoic acid acting via retinoic acid receptor-gamma. J. Biol. Chem. 2010, 285, 19002–19011. [Google Scholar] [CrossRef] [PubMed]
  90. Mukhopadhyay, B.; Schuebel, K.; Mukhopadhyay, P.; Cinar, R.; Godlewski, G.; Xiong, K.; Mackie, K.; Lizak, M.; Yuan, Q.; Goldman, D.; et al. Cannabinoid receptor 1 promotes hepatocellular carcinoma initiation and progression through multiple mechanisms. Hepatology 2015, 61, 1615–1626. [Google Scholar] [CrossRef]
  91. Reece, A.S.; Hulse, G.K. Cannabis Could be the Missing Environmental Pancreatic Carcinogen Hiding in Plain View. Gastroenterology, 2023; in press. [Google Scholar]
  92. WebDigitizer. Available online: https://automeris.io/WebPlotDigitizer/ (accessed on 20 January 2023).
  93. Wickham, H.; Averick, M.; Bryan, J.; Chang, W.; McGowan, L.D.; Francios, R.; Groelmund, G.; Hayes, A.; Henry, L.; Hester, J.; et al. Welcome to the Tidyverse. J. Open. Source Softw. 2019, 4, 1686–1691. [Google Scholar] [CrossRef]
  94. ggpubr: ‘ggplot2’ Based Publication Ready Plots. Available online: https://CRAN.R-project.org/package=ggpubr (accessed on 20 January 2023).
  95. Patchwork: The Composer of Plots_, R Package Version 1.1.2 CRAN. Available online: https://CRAN.R-project.org/package=patchwork (accessed on 20 January 2023).
  96. Broom.Mixed: Tidying Methods for Mixed Models. Available online: http://github.com/bbolker/broom.mixed (accessed on 20 January 2023).
  97. Broom: Convert Statistical Objects into Tidy Tibbles. Available online: https://CRAN.R-project.org/package=broom (accessed on 20 January 2023).
  98. VanderWeele, T.J.; Ding, P. Sensitivity Analysis in Observational Research: Introducing the E-Value. Ann. Intern. Med. 2017, 167, 268–274. [Google Scholar] [CrossRef]
  99. VanderWeele, T.J.; Martin, J.N.; Mathur, M.B. E-values and incidence density sampling. Epidemiology 2020, 31, e51–e52. [Google Scholar] [CrossRef]
  100. VanderWeele, T.J.; Ding, P.; Mathur, M. Technical Considerations in the Use of the E-Value. J. Causal Inference 2019, 7, 1–11. [Google Scholar] [CrossRef]
  101. R Package “Corrplot”: Visualization of a Correlation Matrix. Available online: https://github.com/taiyun/corrplot (accessed on 20 January 2023).
  102. Gaddam, S.; Abboud, Y.; Oh, J.; Samaan, J.S.; Nissen, N.N.; Lu, S.C.; Lo, S.K. Incidence of Pancreatic Cancer by Age and Sex in the US, 2000–2018. Jama 2021, 326, 2075–2077. [Google Scholar] [CrossRef]
  103. National Survey of Drug Use and Health (NSDUH 2018). Available online: https://www.datafiles.samhsa.gov/study/national-survey-drug-use-and-health-nsduh-2018-nid18757# (accessed on 20 January 2023).
  104. Vela, G.; Martin, S.; Garcia-Gil, L.; Crespo, J.A.; Ruiz-Gayo, M.; Fernandez-Ruiz, J.J.; Garcia-Lecumberri, C.; Pelaprat, D.; Fuentes, J.A.; Ramos, J.A.; et al. Maternal exposure to delta9-tetrahydrocannabinol facilitates morphine self-administration behavior and changes regional binding to central mu opioid receptors in adult offspring female rats. Brain Res. 1998, 807, 101–109. [Google Scholar] [CrossRef]
  105. Shoyama, Y.; Sugawa, C.; Tanaka, H.; Morimoto, S. Cannabinoids act as necrosis-inducing factors in Cannabis sativa. Plant. Signal. Behav. 2008, 3, 1111–1112. [Google Scholar] [CrossRef]
  106. Koller, V.J.; Auwarter, V.; Grummt, T.; Moosmann, B.; Misik, M.; Knasmuller, S. Investigation of the In Vitro toxicological properties of the synthetic cannabimimetic drug CP-47,497-C8. Toxicol. Appl. Pharmacol. 2014, 277, 164–171. [Google Scholar] [CrossRef] [PubMed]
  107. Koller, V.J.; Ferk, F.; Al-Serori, H.; Misik, M.; Nersesyan, A.; Auwarter, V.; Grummt, T.; Knasmuller, S. Genotoxic properties of representatives of alkylindazoles and aminoalkyl-indoles which are consumed as synthetic cannabinoids. Food Chem. Toxicol. 2015, 80, 130–136. [Google Scholar] [CrossRef] [PubMed]
  108. Fish, E.W.; Murdaugh, L.B.; Zhang, C.; Boschen, K.E.; Boa-Amponsem, O.; Mendoza-Romero, H.N.; Tarpley, M.; Chdid, L.; Mukhopadhyay, S.; Cole, G.J.; et al. Cannabinoids Exacerbate Alcohol Teratogenesis by a CB1-Hedgehog Interaction. Sci. Rep. 2019, 9, 16057–16075. [Google Scholar] [CrossRef]
  109. Hölzel, B.N.; Pfannkuche, K.; Allner, B.; Allner, H.T.; Hescheler, J.; Derichsweiler, D.; Hollert, H.; Schiwy, A.; Brendt, J.; Schaffeld, M.; et al. Following the adverse outcome pathway from micronucleus to cancer using H2B-eGFP transgenic healthy stem cells. Arch. Toxicol. 2020, 94, 3265–3280. [Google Scholar] [CrossRef]
  110. Tahir, S.K.; Zimmerman, A.M. Influence of marihuana on cellular structures and biochemical activities. Pharmacol. Biochem. Behav. 1991, 40, 617–623. [Google Scholar] [CrossRef]
  111. Sarafian, T.A.; Kouyoumjian, S.; Khoshaghideh, F.; Tashkin, D.P.; Roth, M.D. Delta 9-tetrahydrocannabinol disrupts mitochondrial function and cell energetics. Am. J. Physiol. 2003, 284, L298–L306. [Google Scholar] [CrossRef]
  112. Sarafian, T.A.; Habib, N.; Oldham, M.; Seeram, N.; Lee, R.P.; Lin, L.; Tashkin, D.P.; Roth, M.D. Inhaled marijuana smoke disrupts mitochondrial energetics in pulmonary epithelial cells In Vivo. Am. J. Physiol. 2006, 290, L1202–L1209. [Google Scholar] [CrossRef]
  113. Morimoto, S.; Tanaka, Y.; Sasaki, K.; Tanaka, H.; Fukamizu, T.; Shoyama, Y.; Shoyama, Y.; Taura, F. Identification and characterization of cannabinoids that induce cell death through mitochondrial permeability transition in Cannabis leaf cells. J. Biol. Chem. 2007, 282, 20739–20751. [Google Scholar] [CrossRef]
  114. Fisar, Z.; Singh, N.; Hroudova, J. Cannabinoid-induced changes in respiration of brain mitochondria. Toxicol. Lett. 2014, 231, 62–71. [Google Scholar] [CrossRef]
  115. Singh, N.; Hroudova, J.; Fisar, Z. Cannabinoid-Induced Changes in the Activity of Electron Transport Chain Complexes of Brain Mitochondria. J. Mol. Neurosci. 2015, 56, 926–931. [Google Scholar] [CrossRef]
  116. Pearl, J.; Mackenzie, D. The Book of Why: The New Science of Cause and Effect; Basic Books: New York, NY, USA, 2018; Volume 1. [Google Scholar]
  117. Reece, A.S.; Hulse, G.K. Novel Insights into Potential Cannabis-Related Cancerogenesis from Recent Key Whole Epigenome Screen of Cannabis Dependence and Withdrawal: Epidemiological Comment and Explication of Schrott et.al. Genes 2022, 14, 32. [Google Scholar] [CrossRef]
  118. Wang, J.; Yuan, W.; Li, M.D. Genes and pathways co-associated with the exposure to multiple drugs of abuse, including alcohol, amphetamine/methamphetamine, cocaine, marijuana, morphine, and/or nicotine: A review of proteomics analyses. Mol. Neurobiol. 2011, 44, 269–286. [Google Scholar] [CrossRef]
  119. de Pagter, M.S.; van Roosmalen, M.J.; Baas, A.F.; Renkens, I.; Duran, K.J.; van Binsbergen, E.; Tavakoli-Yaraki, M.; Hochstenbach, R.; van der Veken, L.T.; Cuppen, E.; et al. Chromothripsis in healthy individuals affects multiple protein-coding genes and can result in severe congenital abnormalities in offspring. Am. J. Hum. Genet. 2015, 96, 651–656. [Google Scholar] [CrossRef]
  120. Forment, J.V.; Kaidi, A.; Jackson, S.P. Chromothripsis and cancer: Causes and consequences of chromosome shattering. Nat. Rev. Cancer 2012, 12, 663–670. [Google Scholar] [CrossRef]
  121. Kloosterman, W.P.; Guryev, V.; van Roosmalen, M.; Duran, K.J.; de Bruijn, E.; Bakker, S.C.; Letteboer, T.; van Nesselrooij, B.; Hochstenbach, R.; Poot, M.; et al. Chromothripsis as a mechanism driving complex de novo structural rearrangements in the germline. Hum. Mol. Genet. 2011, 20, 1916–1924. [Google Scholar] [CrossRef]
  122. Kloosterman, W.P.; Hoogstraat, M.; Paling, O.; Tavakoli-Yaraki, M.; Renkens, I.; Vermaat, J.S.; van Roosmalen, M.J.; van Lieshout, S.; Nijman, I.J.; Roessingh, W.; et al. Chromothripsis is a common mechanism driving genomic rearrangements in primary and metastatic colorectal cancer. Genome Biol. 2011, 12, R103. [Google Scholar] [CrossRef]
  123. Liu, S.; Kwon, M.; Mannino, M.; Yang, N.; Renda, F.; Khodjakov, A.; Pellman, D. Nuclear envelope assembly defects link mitotic errors to chromothripsis. Nature 2018, 561, 551–555. [Google Scholar] [CrossRef]
  124. Shoshani, O.; Brunner, S.F.; Yaeger, R.; Ly, P.; Nechemia-Arbely, Y.; Kim, D.H.; Fang, R.; Castillon, G.A.; Yu, M.; Li, J.S.Z.; et al. Chromothripsis drives the evolution of gene amplification in cancer. Nature 2021, 591, 137–141. [Google Scholar] [CrossRef]
  125. Terradas, M.; Martín, M.; Genescà, A. Detection of Impaired DNA Replication and Repair in Micronuclei as Indicators of Genomic Instability and Chromothripsis. Methods Mol. Biol. 2018, 1769, 197–208. [Google Scholar] [CrossRef]
  126. Waldron, D. Genome stability: Chromothripsis and micronucleus formation. Nat. Rev. Genet. 2015, 16, 376–377. [Google Scholar] [CrossRef]
  127. Zhang, C.Z.; Spektor, A.; Cornils, H.; Francis, J.M.; Jackson, E.K.; Liu, S.; Meyerson, M.; Pellman, D. Chromothripsis from DNA damage in micronuclei. Nature 2015, 522, 179–184. [Google Scholar] [CrossRef] [PubMed]
  128. Mackenzie, K.J.; Carroll, P.; Martin, C.A.; Murina, O.; Fluteau, A.; Simpson, D.J.; Olova, N.; Sutcliffe, H.; Rainger, J.K.; Leitch, A.; et al. cGAS surveillance of micronuclei links genome instability to innate immunity. Nature 2017, 548, 461–465. [Google Scholar] [CrossRef] [PubMed]
  129. Zhang, C.; Shang, G.; Gui, X.; Zhang, X.; Bai, X.C.; Chen, Z.J. Structural basis of STING binding with and phosphorylation by TBK1. Nature 2019, 567, 394–398. [Google Scholar] [CrossRef] [PubMed]
  130. Zhang, X.; Bai, X.C.; Chen, Z.J. Structures and Mechanisms in the cGAS-STING Innate Immunity Pathway. Immunity 2020, 53, 43–53. [Google Scholar] [CrossRef] [PubMed]
  131. Domizio, J.D.; Gulen, M.F.; Saidoune, F.; Thacker, V.V.; Yatim, A.; Sharma, K.; Nass, T.; Guenova, E.; Schaller, M.; Conrad, C.; et al. The cGAS-STING pathway drives type I IFN immunopathology in COVID-19. Nature 2022, 603, 145–151. [Google Scholar] [CrossRef]
  132. Hong, C.; Schubert, M.; Tijhuis, A.E.; Requesens, M.; Roorda, M.; van den Brink, A.; Ruiz, L.A.; Bakker, P.L.; van der Sluis, T.; Pieters, W.; et al. cGAS-STING drives the IL-6-dependent survival of chromosomally instable cancers. Nature 2022, 607, 366–373. [Google Scholar] [CrossRef]
  133. Li, S.; Mirlekar, B.; Johnson, B.M.; Brickey, W.J.; Wrobel, J.A.; Yang, N.; Song, D.; Entwistle, S.; Tan, X.; Deng, M.; et al. STING-induced regulatory B cells compromise NK function in cancer immunity. Nature 2022, 610, 373–380. [Google Scholar] [CrossRef]
  134. Tang, F.; Xu, D.; Wang, S.; Wong, C.K.; Martinez-Fundichely, A.; Lee, C.J.; Cohen, S.; Park, J.; Hill, C.E.; Eng, K.; et al. Chromatin profiles classify castration-resistant prostate cancers suggesting therapeutic targets. Science 2022, 376, eabe1505. [Google Scholar] [CrossRef]
  135. Wan, D.; Jiang, W.; Hao, J. Research Advances in How the cGAS-STING Pathway Controls the Cellular Inflammatory Response. Front. Immunol. 2020, 11, 615. [Google Scholar] [CrossRef]
  136. Xiao, T.S.; Fitzgerald, K.A. The cGAS-STING pathway for DNA sensing. Mol. Cell 2013, 51, 135–139. [Google Scholar] [CrossRef]
  137. Yang, H.; Wang, H.; Ren, J.; Chen, Q.; Chen, Z.J. cGAS is essential for cellular senescence. Proc. Natl. Acad. Sci. USA 2017, 114, E4612–E4620. [Google Scholar] [CrossRef] [PubMed]
  138. Gluck, S.; Guey, B.; Gulen, M.F.; Wolter, K.; Kang, T.W.; Schmacke, N.A.; Bridgeman, A.; Rehwinkel, J.; Zender, L.; Ablasser, A. Innate immune sensing of cytosolic chromatin fragments through cGAS promotes senescence. Nat. Cell Biol. 2017, 19, 1061–1070. [Google Scholar] [CrossRef] [PubMed]
  139. Sladitschek-Martens, H.L.; Guarnieri, A.; Brumana, G.; Zanconato, F.; Battilana, G.; Xiccato, R.L.; Panciera, T.; Forcato, M.; Bicciato, S.; Guzzardo, V.; et al. YAP/TAZ activity in stromal cells prevents ageing by controlling cGAS-STING. Nature 2022, 607, 790–798. [Google Scholar] [CrossRef]
  140. Dou, Z.; Ghosh, K.; Vizioli, M.G.; Zhu, J.; Sen, P.; Wangensteen, K.J.; Simithy, J.; Lan, Y.; Lin, Y.; Zhou, Z.; et al. Cytoplasmic chromatin triggers inflammation in senescence and cancer. Nature 2017, 550, 402–406. [Google Scholar] [CrossRef] [PubMed]
  141. Wang, T.W.; Johmura, Y.; Suzuki, N.; Omori, S.; Migita, T.; Yamaguchi, K.; Hatakeyama, S.; Yamazaki, S.; Shimizu, E.; Imoto, S.; et al. Blocking PD-L1-PD-1 improves senescence surveillance and ageing phenotypes. Nature 2022, 611, 358–364. [Google Scholar] [CrossRef] [PubMed]
  142. Olivieri, F.; Bonafe, M.; Giovagnetti, S.; Stecconi, R.; Cardelli, M.; Cavallone, L.; Spazzafumo, L.; Marchegiani, F.; Carrieri, G.; Mugianesi, E.; et al. In Vitro IL-6 production by EBV-immortalized B lymphocytes from young and elderly people genotyped for -174 C/G polymorphism in IL-6 gene: A model to study the genetic basis of inflamm-aging. Mech. Ageing Dev. 2003, 124, 549–553. [Google Scholar] [CrossRef]
  143. Salvioli, S.; Capri, M.; Valensin, S.; Tieri, P.; Monti, D.; Ottaviani, E.; Franceschi, C. Inflamm-aging, cytokines and aging: State of the art, new hypotheses on the role of mitochondria and new perspectives from systems biology. Curr. Pharm. Des. 2006, 12, 3161–3171. [Google Scholar] [CrossRef]
  144. Chazaud, B.; Mouchiroud, G. Inflamm-aging: STAT3 signaling pushes muscle stem cells off balance. Cell Stem Cell 2014, 15, 401–402. [Google Scholar] [CrossRef]
  145. Morrisette-Thomas, V.; Cohen, A.A.; Fulop, T.; Riesco, E.; Legault, V.; Li, Q.; Milot, E.; Dusseault-Belanger, F.; Ferrucci, L. Inflamm-aging does not simply reflect increases in pro-inflammatory markers. Mech. Ageing Dev. 2014, 139, 49–57. [Google Scholar] [CrossRef]
  146. Kovtonyuk, L.V.; Fritsch, K.; Feng, X.; Manz, M.G.; Takizawa, H. Inflamm-Aging of Hematopoiesis, Hematopoietic Stem Cells, and the Bone Marrow Microenvironment. Front. Immunol. 2016, 7, 502. [Google Scholar] [CrossRef]
  147. Pangrazzi, L.; Meryk, A.; Naismith, E.; Koziel, R.; Lair, J.; Krismer, M.; Trieb, K.; Grubeck-Loebenstein, B. “Inflamm-aging” influences immune cell survival factors in human bone marrow. Eur. J. Immunol. 2016, 47, 481–492. [Google Scholar] [CrossRef] [PubMed]
  148. Bakhoum, S.F.; Ngo, B.; Laughney, A.M.; Cavallo, J.A.; Murphy, C.J.; Ly, P.; Shah, P.; Sriram, R.K.; Watkins, T.B.K.; Taunk, N.K.; et al. Chromosomal instability drives metastasis through a cytosolic DNA response. Nature 2018, 553, 467–472. [Google Scholar] [CrossRef] [PubMed]
  149. Donald, P.J. Marijuana smoking--possible cause of head and neck carcinoma in young patients. Otolaryngol. Head. Neck Surg. 1986, 94, 517–521. [Google Scholar] [CrossRef] [PubMed]
  150. Graef, S.; Choo, C.G.; Warfield, A.; Cullen, M.; Woolhouse, I. Small cell lung cancer in a 26-year-old man with significant Cannabis exposure. J. Thorac. Oncol. 2011, 6, 218–219. [Google Scholar] [CrossRef]
  151. Kothadia, J.P.; Chhabra, S.; Marcus, A.; May, M.; Saraiya, B.; Jabbour, S.K. Anterior mediastinal mass in a young marijuana smoker: A rare case of small-cell lung cancer. Case Rep. Med. 2012, 2012, 754231. [Google Scholar] [CrossRef]
  152. Van Hoozen, B.E.; Cross, C.E. Marijuana. Respiratory tract effects. Clin. Rev. Allergy Immunol. 1997, 15, 243–269. [Google Scholar] [CrossRef]
  153. Cabral, G.A.; Raborn, E.S.; Griffin, L.; Dennis, J.; Marciano-Cabral, F. CB(2) receptors in the brain: Role in central immune function. Br. J. Pharm. 2008, 153, 240–251. [Google Scholar] [CrossRef]
  154. De Meij, J.; Alfanek, Z.; Morel, L.; Decoeur, F.; Leyrolle, Q.; Picard, K.; Carrier, M.; Aubert, A.; Séré, A.; Lucas, C.; et al. Microglial Cannabinoid Type 1 Receptor Regulates Brain Inflammation in a Sex-Specific Manner. Cannabis Cannabinoid Res. 2021, 6, 488–507. [Google Scholar] [CrossRef]
  155. Ge, Q.; Maury, E.; Rycken, L.; Gerard, J.; Noel, L.; Detry, R.; Navez, B.; Brichard, S.M. Endocannabinoids regulate adipokine production and the immune balance of omental adipose tissue in human obesity. Int. J. Obes. 2013, 37, 874–880. [Google Scholar] [CrossRef]
  156. Mandolesi, G.; Bullitta, S.; Fresegna, D.; Gentile, A.; De Vito, F.; Dolcetti, E.; Rizzo, F.R.; Strimpakos, G.; Centonze, D.; Musella, A. Interferon-gamma causes mood abnormalities by altering cannabinoid CB1 receptor function in the mouse striatum. Neurobiol. Dis. 2017, 108, 45–53. [Google Scholar] [CrossRef]
  157. Pacher, P.; Steffens, S.; Hasko, G.; Schindler, T.H.; Kunos, G. Cardiovascular effects of marijuana and synthetic cannabinoids: The good, the bad, and the ugly. Nat. Rev. Cardiol. 2018, 15, 151–166. [Google Scholar] [CrossRef] [PubMed]
  158. Vujic, N.; Schlager, S.; Eichmann, T.O.; Madreiter-Sokolowski, C.T.; Goeritzer, M.; Rainer, S.; Schauer, S.; Rosenberger, A.; Woelfler, A.; Doddapattar, P.; et al. Monoglyceride lipase deficiency modulates endocannabinoid signaling and improves plaque stability in ApoE-knockout mice. Atherosclerosis 2016, 244, 9–21. [Google Scholar] [CrossRef] [PubMed]
  159. Mon, M.J.; Jansing, R.L.; Doggett, S.; Stein, J.L.; Stein, G.S. Influence of delta9-tetrahydrocannabinol on cell proliferation and macromolecular biosynthesis in human cells. Biochem. Pharmacol. 1978, 27, 1759–1765. [Google Scholar] [CrossRef] [PubMed]
  160. Mon, M.J.; Haas, A.E.; Stein, J.L.; Stein, G.S. Influence of psychoactive and nonpsychoactive cannabinoids on cell proliferation and macromolecular biosynthesis in human cells. Biochem. Pharmacol. 1981, 30, 31–43. [Google Scholar] [CrossRef] [PubMed]
  161. Mon, M.J.; Haas, A.E.; Stein, J.L.; Stein, G.S. Influence of psychoactive and nonpsychoactive cannabinoids on chromatin structure and function in human cells. Biochem. Pharmacol. 1981, 30, 45–58. [Google Scholar] [CrossRef]
  162. Blevins, R.D.; Regan, J.D. delta-9-Tetrahydrocannabinol: Effect on macromolecular synthesis in human and other mammalian cells. Arch. Toxicol. 1976, 35, 127–135. [Google Scholar] [CrossRef]
  163. McClean, D.K.; Zimmerman, A.M. Action of delta 9-tetrahydrocannabinol on cell division and macromolecular synthesis in division-synchronized protozoa. Pharmacology 1976, 14, 307–321. [Google Scholar] [CrossRef]
  164. Costa, B.; Colleoni, M. Changes in rat brain energetic metabolism after exposure to anandamide or Delta(9)-tetrahydrocannabinol. Eur. J. Pharmacol. 2000, 395, 1–7. [Google Scholar] [CrossRef]
  165. Rossato, M.; Ion Popa, F.; Ferigo, M.; Clari, G.; Foresta, C. Human sperm express cannabinoid receptor Cb1, the activation of which inhibits motility, acrosome reaction, and mitochondrial function. J. Clin. Endocrinol. Metab. 2005, 90, 984–991. [Google Scholar] [CrossRef]
  166. Badawy, Z.S.; Chohan, K.R.; Whyte, D.A.; Penefsky, H.S.; Brown, O.M.; Souid, A.K. Cannabinoids inhibit the respiration of human sperm. Fertil. Steril. 2009, 91, 2471–2476. [Google Scholar] [CrossRef]
  167. Volkow, N.D.; Baler, R.D.; Compton, W.M.; Weiss, S.R. Adverse Health Effects of Marijuana Use. N. Engl. J. Med. 2014, 371, 878–879. [Google Scholar] [CrossRef] [PubMed]
  168. Wolff, V.; Rouyer, O.; Schlagowski, A.; Zoll, J.; Raul, J.S.; Marescaux, C. Study of the effect of THC on mitochondrial respiration of the rat brain. Rev. Neurol. Neurol. Rev. 2014, 170, A19–A20. [Google Scholar] [CrossRef]
  169. Wolff, V.; Schlagowski, A.I.; Rouyer, O.; Charles, A.L.; Singh, F.; Auger, C.; Schini-Kerth, V.; Marescaux, C.; Raul, J.S.; Zoll, J.; et al. Tetrahydrocannabinol induces brain mitochondrial respiratory chain dysfunction and increases oxidative stress: A potential mechanism involved in cannabis-related stroke. Biomed. Res. Int. 2015, 2015, 323706. [Google Scholar] [CrossRef] [PubMed]
  170. Rupprecht, A.; Theisen, U.; Wendt, F.; Frank, M.; Hinz, B. The Combination of Δ(9)-Tetrahydrocannabinol and Cannabidiol Suppresses Mitochondrial Respiration of Human Glioblastoma Cells via Downregulation of Specific Respiratory Chain Proteins. Cancers 2022, 14, 3129. [Google Scholar] [CrossRef] [PubMed]
  171. Alasmari, F.; Assiri, M.A.; Ahamad, S.R.; Aljumayi, S.R.; Alotaibi, W.H.; Alhamdan, M.M.; Alhazzani, K.; Alharbi, M.; Alqahtani, F.; Alasmari, A.F. Serum Metabolomic Analysis of Male Patients with Cannabis or Amphetamine Use Disorder. Metabolites 2022, 12, 179. [Google Scholar] [CrossRef]
  172. Jimenez-Blasco, D.; Busquets-Garcia, A.; Hebert-Chatelain, E.; Serrat, R.; Vicente-Gutierrez, C.; Ioannidou, C.; Gómez-Sotres, P.; Lopez-Fabuel, I.; Resch-Beusher, M.; Resel, E.; et al. Glucose metabolism links astroglial mitochondria to cannabinoid effects. Nature 2020, 583, 603–608. [Google Scholar] [CrossRef]
  173. Magistretti, P.J. How lactate links cannabis to social behaviour. Nature 2020, 583, 526–527. [Google Scholar] [CrossRef]
  174. Zhou, X.; Zhang, N.; Zhang, J.; Gu, Q.; Dong, C.; Lin, B.; Zou, C. Microbiome and fermentation parameters in the rumen of dairy buffalo in response to ingestion associated with a diet supplemented with cysteamine and hemp seed oil. J. Anim. Physiol. Anim. Nutr. 2022, 106, 471–484. [Google Scholar] [CrossRef]
  175. Papadakis, D.P.; Michael, C.M.; Kephalas, T.A.; Miras, C.J. Effects of cannabis smoking in blood lactic acid and glucose in humans. Experientia 1974, 30, 1183–1184. [Google Scholar] [CrossRef]
  176. Li, X.; Yang, Y.; Zhang, B.; Lin, X.; Fu, X.; An, Y.; Zou, Y.; Wang, J.-X.; Wang, Z.; Yu, T. Lactate metabolism in human health and disease. Signal. Transduct. Target. Ther. 2022, 7, 305. [Google Scholar] [CrossRef]
  177. Yang, Z.; Yan, C.; Ma, J.; Peng, P.; Ren, X.; Cai, S.; Shen, X.; Wu, Y.; Zhang, S.; Wang, X.; et al. Lactylome analysis suggests lactylation-dependent mechanisms of metabolic adaptation in hepatocellular carcinoma. Nat. Metab. 2023, 5, 61–79. [Google Scholar] [CrossRef] [PubMed]
  178. Flores, A.; Schell, J.; Krall, A.S.; Jelinek, D.; Miranda, M.; Grigorian, M.; Braas, D.; White, A.C.; Zhou, J.L.; Graham, N.A.; et al. Lactate dehydrogenase activity drives hair follicle stem cell activation. Nat. Cell Biol. 2017, 19, 1017–1026. [Google Scholar] [CrossRef] [PubMed]
  179. Martin, L.J.; Cairns, E.A.; Heblinski, M.; Fletcher, C.; Krycer, J.R.; Arnold, J.C.; McGregor, I.S.; Bowen, M.T.; Anderson, L.L. Cannabichromene and Δ(9)-Tetrahydrocannabinolic Acid Identified as Lactate Dehydrogenase-A Inhibitors by In Silico and In Vitro Screening. J. Nat. Prod. 2021, 84, 1469–1477. [Google Scholar] [CrossRef] [PubMed]
  180. Stern, M. Evidence that a mitochondrial death spiral underlies antagonistic pleiotropy. Aging Cell 2017, 16, 435–443. [Google Scholar] [CrossRef] [PubMed]
  181. Goto, M. Inflammaging (inflammation + aging): A driving force for human aging based on an evolutionarily antagonistic pleiotropy theory? Biosci. Trends 2008, 2, 218–230. [Google Scholar]
  182. Lauc, G.; Huffman, J.E.; Pucic, M.; Zgaga, L.; Adamczyk, B.; Muzinic, A.; Novokmet, M.; Polasek, O.; Gornik, O.; Kristic, J.; et al. Loci associated with N-glycosylation of human immunoglobulin G show pleiotropy with autoimmune diseases and haematological cancers. PLoS Genet. 2013, 9, e1003225. [Google Scholar] [CrossRef]
  183. Waddington, C.H. Organizers and Genes; Cambridge University Press: Cambridge, UK, 1940; Volume 1. [Google Scholar]
  184. Andreeva, L.; Hiller, B.; Kostrewa, D.; Lassig, C.; de Oliveira Mann, C.C.; Jan Drexler, D.; Maiser, A.; Gaidt, M.; Leonhardt, H.; Hornung, V.; et al. cGAS senses long and HMGB/TFAM-bound U-turn DNA by forming protein-DNA ladders. Nature 2017, 549, 394–398. [Google Scholar] [CrossRef]
  185. Luecke, S.; Holleufer, A.; Christensen, M.H.; Jønsson, K.L.; Boni, G.A.; Sørensen, L.K.; Johannsen, M.; Jakobsen, M.R.; Hartmann, R.; Paludan, S.R. cGAS is activated by DNA in a length-dependent manner. EMBO Rep. 2017, 18, 1707–1715. [Google Scholar] [CrossRef]
  186. Ablasser, A.; Chen, Z.J. cGAS in action: Expanding roles in immunity and inflammation. Science 2019, 363, eaat8657. [Google Scholar] [CrossRef]
  187. Xia, P.; Wang, S.; Gao, P.; Gao, G.; Fan, Z. DNA sensor cGAS-mediated immune recognition. Protein Cell 2016, 7, 777–791. [Google Scholar] [CrossRef]
  188. Gumbs, A.A.; Gogol, M.; Spolverato, G.; Taher, H.; Chouillard, E.K. Systematic Review of the Integrative Medicine Recommendations for Patients with Pancreatic Cancer. Surgeries 2021, 2, 216–230. [Google Scholar] [CrossRef]
  189. Reece, A.S. Cannabinoid Epigenotoxicity and Genotoxicity. JAMA 2022, E3. [Google Scholar] [CrossRef]
  190. Reece, A.S.; Hulse, G.K. Extending the “Paracentral Dogma” of biology with the metabolome: Implications for understanding genomic-glycomic-metabolic-epigenomic synchronization. Engineering, 2022; in press. [Google Scholar] [CrossRef]
  191. Reece, A.S.; Hulse, G.K. Impacts of cannabinoid epigenetics on human development: Reflections on Murphy et al. ‘cannabinoid exposure and altered DNA methylation in rat and human sperm’ epigenetics 2018; 13: 1208–1221. Epigenetics 2019, 14, 1041–1056. [Google Scholar] [CrossRef]
  192. Reece, A.S.; Hulse, G.K. Epidemiology of Cannabis: Genotoxicity, Neurotoxicity, Epigenomics and Aging; Elsevier: New York, NY, USA, 2023; Volume 1. [Google Scholar]
  193. Reece, A.S.; Hulse, G.K. Cannabis, Cannabidiol, Cannabinoids and Multigenerational Policy. Engineering, 2022; in press. [Google Scholar] [CrossRef]
Figure 1. Annual incidence of pancreatic cancer with (A) linear scale and loess fitted curves, (B) linear scale and quadratic fitted curve, (C) log scale and loess fitted curve and (D) log scale and quadratic fitted curves.
Figure 1. Annual incidence of pancreatic cancer with (A) linear scale and loess fitted curves, (B) linear scale and quadratic fitted curve, (C) log scale and loess fitted curve and (D) log scale and quadratic fitted curves.
Gastroent 14 00016 g001
Figure 2. Tobacco Consumption by Sociodemographic Strata.
Figure 2. Tobacco Consumption by Sociodemographic Strata.
Gastroent 14 00016 g002
Figure 3. Alcohol Consumption by Sociodemographic Strata.
Figure 3. Alcohol Consumption by Sociodemographic Strata.
Gastroent 14 00016 g003
Figure 4. Near Daily Cannabis Use by Sociodemographic Strata.
Figure 4. Near Daily Cannabis Use by Sociodemographic Strata.
Gastroent 14 00016 g004
Figure 5. Substance Use by Age Group as younger or older than 50 years.
Figure 5. Substance Use by Age Group as younger or older than 50 years.
Gastroent 14 00016 g005
Figure 6. Substance Exposure by Sociodemographic Strata.
Figure 6. Substance Exposure by Sociodemographic Strata.
Gastroent 14 00016 g006
Figure 7. Correlogram for Caucasian American males 15–34 years (left) Pearson correlation coefficient and (right) significance levels.
Figure 7. Correlogram for Caucasian American males 15–34 years (left) Pearson correlation coefficient and (right) significance levels.
Gastroent 14 00016 g007
Figure 8. Correlogram for African American males 15–34 years (left) Pearson correlation coefficient and (right) significance levels.
Figure 8. Correlogram for African American males 15–34 years (left) Pearson correlation coefficient and (right) significance levels.
Gastroent 14 00016 g008
Figure 9. Correlogram for Caucasian American females 15–34 years (left) Pearson correlation coefficient and (right) significance levels.
Figure 9. Correlogram for Caucasian American females 15–34 years (left) Pearson correlation coefficient and (right) significance levels.
Gastroent 14 00016 g009
Figure 10. Correlogram for African American females 15–34 years (left) Pearson correlation coefficient and (right) significance levels.
Figure 10. Correlogram for African American females 15–34 years (left) Pearson correlation coefficient and (right) significance levels.
Gastroent 14 00016 g010
Table 1. Stratified Correlation Matrix.
Table 1. Stratified Correlation Matrix.
CovariatePDAC.RateHOPTobaccoAlcoholCannabis
Males, White
PDAC.Rate10.8225−0.7732−0.27120.8308
HOP0.82251−0.6628−0.32570.7015
Tobacco−0.7732−0.662810.7601−0.3816
Alcohol−0.2712−0.32570.760110.1657
Cannabis0.83080.7015−0.38160.16571
Males, Black
PDAC.Rate10.8225−0.6820−0.25240.8312
HOP0.82251−0.4851−0.25690.6359
Tobacco−0.6820−0.485110.5705−0.3691
Alcohol−0.2524−0.25690.570510.0738
Cannabis0.83120.6359−0.36910.07381
Males, Hispanic
PDAC.Rate10.8225−0.81910.14640.8980
HOP0.82251−0.71250.10520.7840
Tobacco−0.8191−0.712510.2834−0.6757
Alcohol0.14640.10520.283410.3221
Cannabis0.89800.7840−0.67570.32211
Females, White
PDAC.Rate10.9350−0.88650.34750.9309
HOP0.93501−0.80520.44270.9040
Tobacco−0.8865−0.805210.0211−0.7555
Alcohol0.34750.44270.021110.5538
Cannabis0.93090.9040−0.75550.55381
Females, Black
PDAC.Rate10.9350−0.38780.44600.8539
HOP0.93501−0.28960.52530.8421
Tobacco−0.3878−0.289610.3687−0.1939
Alcohol0.44600.52530.368710.5877
Cannabis0.85390.8421−0.19390.58771
Females, Hispanic
PDAC.Rate10.9350−0.81510.63110.8982
HOP0.93501−0.71380.69120.9133
Tobacco−0.8151−0.71381−0.4312−0.6870
Alcohol0.63110.6912−0.431210.7699
Cannabis0.89820.9133−0.68700.76991
Table 2. Significance Levels (p-Values) of Correlation Matrix.
Table 2. Significance Levels (p-Values) of Correlation Matrix.
CovariatePDAC.RateHOPTobaccoAlcoholCannabis
Males, White
PDAC.Rate00.00290.00230.04390.0063
HOP0.002900.00150.02100.0229
Tobacco0.00230.001500.01150.0303
Alcohol0.04390.02100.011500.1369
Cannabis0.00630.02290.03030.13690
Males, Black
PDAC.Rate00.00520.00020.06070.0068
HOP0.005200.00410.03040.0376
Tobacco0.00020.004100.05530.0139
Alcohol0.06070.03040.055300.1487
Cannabis0.00680.03760.01390.14870
Males, Hispanic
PDAC.Rate00.00280.00140.63340.0009
HOP0.002800.00210.56260.0068
Tobacco0.00140.002100.47030.0073
Alcohol0.63340.56260.470300.7493
Cannabis0.00090.00680.00730.74930
Females, White
PDAC.Rate00.00020.00100.45400.0008
HOP0.000200.00280.40890.0005
Tobacco0.00100.002800.60000.0068
Alcohol0.45400.40890.600000.3435
Cannabis0.00080.00050.00680.34350
Females, Black
PDAC.Rate00.00030.00130.83240.0051
HOP0.000300.00380.77780.0056
Tobacco0.00130.003800.98410.0147
Alcohol0.83240.77780.984100.6839
Cannabis0.00510.00560.01470.68390
Females, Hispanic
PDAC.Rate00.00010.00000.03440.0008
HOP0.000100.00020.02800.0005
Tobacco0.00000.000200.03510.0010
Alcohol0.03440.02800.035100.0166
Cannabis0.00080.00050.00100.01660
Table 3. Bivariate Regression Results for 15–34-Year Cohort.
Table 3. Bivariate Regression Results for 15–34-Year Cohort.
SubsRaceSexβ-Estimate (S.E.)p-ValueP.Adj.HolmP.Adj.FDRE-Value EstimateE-Value Lower Bound
CannabisWhiteFemale29.05 (22.7, 35.4)7.00 × 10−81.19 × 10−66.30 × 10−73.46 × 10538.41 × 1041
CannabisHispanicFemale32.1 (24.45, 39.75)2.50 × 10−74.00 × 10−61.50 × 10−61.32 × 10551.20 × 1042
CannabisBlackFemale19.67 (13.71, 25.63)5.80 × 10−66.96 × 10−51.49 × 10−52.04 × 10286.88 × 1019
AlcoholHispanicFemale6.28 (2.62, 9.94)0.00370.02960.00603.19 × 106785.71
AlcoholBlackFemale4.76 (0.84, 8.68)0.02920.17500.04043.26 × 10410.81
AlcoholWhiteFemale6.01 (−0.36, 12.37)0.08170.40870.10512.20 × 1051.00
TobaccoBlackFemale−6.67 (−16.72, 3.38)0.21070.84280.25294.53 × 105-
TobaccoHispanicFemale−12.6 (−17.52, −7.69)1.04 × 10−41.04 × 10−32.08 × 10−43.27 × 1015-
TobaccoWhiteFemale−8.37 (−10.83, −5.91)3.97 × 10−65.16 × 10−51.19 × 10−53.02 × 1012-
CannabisHispanicMale15.44 (12.28, 18.6)3.00 × 10−85.40 × 10−75.40 × 10−72.97 × 10384.73 × 1030
CannabisWhiteMale14.66 (10.68, 18.64)1.43 × 10−62.15 × 10−56.44 × 10−61.64 × 10292.42 × 1021
CannabisBlackMale11.35 (8.05, 14.64)3.43 × 10−64.80 × 10−51.19 × 10−53.74 × 10212.53 × 1015
AlcoholHispanicMale2.57 (−2.55, 7.7)0.33891.00000.3813770.98971.00
AlcoholBlackMale−1.09 (−4.91, 2.73)0.58381.00000.583823.24-
AlcoholWhiteMale−1.82 (−6.48, 2.83)0.45301.00000.4797129.86-
TobaccoBlackMale−6.21 (−9.93, −2.49)0.00450.03150.00671.14 × 108-
TobaccoWhiteMale−4.83 (−7.1, −2.56)6.47 × 10−40.00580.00121.02 × 107-
TobaccoHispanicMale−6.51 (−8.83, −4.2)3.72 × 10−50.00048.36 × 10−58.94 × 1010-
Table 4. Bivariate Regression Results for Localized Disease.
Table 4. Bivariate Regression Results for Localized Disease.
SubstanceRaceSexβ-Estimate (S.E.)p-ValueP.Adj.HolmP.Adj.FDRE-Value EstimateE-Value Lower Bound
CannabisHispanicFemale24.47 (18.25, 30.69)6.00 × 10−71.08 × 10−51.08 × 10−55.10 × 10514.63 × 1038
CannabisWhiteFemale21.48 (15.68, 27.29)1.34 × 10−62.28 × 10−51.21 × 10−52.41 × 10435.86 × 1031
CannabisBlackFemale14.69 (9.72, 19.66)2.15 × 10−50.00030.00012.56 × 10258.62 × 1016
AlcoholHispanicFemale5.16 (2.46, 7.86)0.00160.01790.00371.59 × 1073.91 × 103
AlcoholBlackFemale4.23 (1.36, 7.09)0.01010.07070.01522.61 × 10590.28
AlcoholWhiteFemale5.73 (1.07, 10.38)0.02750.16480.03807.45 × 10634.35
TobaccoBlackFemale−3.32 (−11.31, 4.67)0.42670.90560.45184.51 × 103-
TobaccoHispanicFemale−8.13 (−12.71, −3.56)0.00280.02540.00517.14 × 1010-
TobaccoWhiteFemale−5.64 (−8.06, −3.22)2.76 × 10−40.00390.00104.34 × 108-
CannabisHispanicMale7.06 (4.35, 9.77)8.74 × 10−50.00130.00044.62 × 10207.36 × 1012
CannabisWhiteMale6.37 (3.28, 9.47)8.60 × 10−40.01030.00222.84 × 10164.20 × 108
CannabisBlackMale4.49 (1.84, 7.15)0.00410.03300.00675.60 × 10103.79 × 104
AlcoholHispanicMale0.56 (−2.27, 3.39)0.70320.90560.703220.35451.00
AlcoholBlackMale−1.09 (−3.11, 0.92)0.30190.90560.3396223.3997-
AlcoholWhiteMale−1.79 (−4.2, 0.62)0.16450.65800.19745.40 × 103-
TobaccoBlackMale−2.61 (−4.85, −0.36)0.03610.18050.04644.95 × 105-
TobaccoWhiteMale−2.42 (−3.73, −1.11)0.00210.02110.00421.35 × 106-
TobaccoHispanicMale−3.23 (−4.64, −1.82)3.28 × 10−40.00430.00108.74 × 108-
Table 5. Additive Multivariable Regression for Substances for 15–34-Year Age Group.
Table 5. Additive Multivariable Regression for Substances for 15–34-Year Age Group.
RaceSexSubstanceβ-Estimate (S.E.)p-ValueP.Adj.HolmP.Adj.FDRE-Value EstimateE-Value Lower Bound
BlackFemaleCannabis18.17 (7.87, 28.46)0.00380.04610.00995.59 × 10252.27 × 1011
BlackFemaleAlcohol1.36 (−2.59, 5.31)0.510510.57431.60 × 1021.00
BlackFemaleTobacco−5.37 (−12.79, 2.04)0.17750.88760.22826.77 × 107-
HispanicFemaleCannabis27.84 (13.14, 42.54)0.00230.03250.00844.81 × 10501.34 × 1024
HispanicFemaleAlcohol−0.77 (−3.99, 2.44)0.64431.00000.68224.99 × 101-
HispanicFemaleTobacco−4.95 (−9.64, −0.27)0.05720.51470.10291.85 × 109-
WhiteFemaleCannabis14.28 (−1.19, 29.76)0.09190.73500.13939.44 × 10371.00
WhiteFemaleAlcohol2.71 (−1.74, 7.17)0.252310.30282.88 × 1071.00
WhiteFemaleTobacco−5.99 (−9.58, −2.39)0.00570.06240.01281.23 × 1016-
BlackMaleCannabis9.63 (6.34, 12.91)5.14 × 10−59.24 × 10−44.63 × 10−44.20 × 10212.26 × 1014
BlackMaleAlcohol0.11 (−2.17, 2.39)0.928410.92842.83 × 1001.00
BlackMaleTobacco−3.96 (−7.09, −0.84)0.02620.26170.05231.19 × 109-
HispanicMaleCannabis8.97 (4.05, 13.89)0.00310.03990.00921.33 × 10311.75 × 1014
HispanicMaleAlcohol2.16 (−0.28, 4.6)0.10470.73500.14495.25 × 1071.00
HispanicMaleTobacco−4.38 (−6.53, −2.23)0.00130.02020.00602.24 × 1015-
WhiteMaleCannabis9.9 (6.52, 13.29)0.00010.00090.00052.84 × 10382.80 × 1025
WhiteMaleAlcohol2.65 (−0.23, 5.53)0.09290.734975760.13933.23 × 10101.00
WhiteMaleTobacco−5.18 (−7.37, −2.99)3.81 × 10−40.00609760.00231.83 × 1020-
Table 6. Multivariable Interactive Regression Results for 15–34-Year Age Group.
Table 6. Multivariable Interactive Regression Results for 15–34-Year Age Group.
RaceSexß-Estimate (S.E.)p-ValueP.Adj.HolmP.Adj.FDRE-Value EstimateE-Value Lower Bound
BlackFemale37.76 (25.24, 50.28)2.19 × 10−56.57 × 10−53.29 × 10−51.60 × 10542.33 × 1036
HispanicFemale61.97 (45.04, 78.9)2.21 × 10−68.84 × 10−64.42 × 10−61.26 × 101022.20 × 1074
WhiteFemale48.15 (37.05, 59.26)2.50 × 10−71.50 × 10−61.41 × 10−65.70 × 10908.62 × 1069
BlackMale17.46 (10.63, 24.29)1.28 × 10−42.56 × 10−41.38 × 10−48.62 × 10285.71178 × 1017
HispanicMale25.9 (19.63, 32.17)4.70 × 10−72.35 × 10−61.41 × 10−61.38 × 10604.84 × 1045
WhiteMale18.26 (11.07, 25.46)1.38 × 10−42.56 × 10−41.38 × 10−41.32 × 10302.48 × 1018
Table 7. Multivariable Interactive Regression Results for Localized Disease.
Table 7. Multivariable Interactive Regression Results for Localized Disease.
RaceSexß-Estimate (S.E.)p-ValueP.Adj.HolmP.Adj.FDRE-Value EstimateE-Value Lower Bound
BlackFemale27.2 (18.43, 35.98)1.24 × 10−54.94 × 10−52.47 × 10−51.98 × 10486.99 × 1032
HispanicFemale46.53 (34.52, 58.54)7.40 × 10−74.44 × 10−63.09 × 10−67.76 × 10961.02 × 1072
WhiteFemale32.85 (24.15, 41.54)1.03 × 10−65.15 × 10−63.09 × 10−61.59 × 10673.41 × 1049
BlackMale6.94 (1.77, 12.11)0.01740.01740.01741.67 × 10151.37 × 104
HispanicMale11.95 (7.33, 16.57)9.42 × 10−52.83 × 10−41.41 × 10−44.26 × 10342.46 × 1021
WhiteMale8.14 (3.28, 12.99)0.00440.00870.00521.44 × 10198.63 × 107
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Reece, A.S.; Hulse, G.K. Sociodemographically Stratified Exploration of Pancreatic Cancer Incidence in Younger US Patients: Implication of Cannabis Exposure as a Risk Factor. Gastroenterol. Insights 2023, 14, 204-235. https://doi.org/10.3390/gastroent14020016

AMA Style

Reece AS, Hulse GK. Sociodemographically Stratified Exploration of Pancreatic Cancer Incidence in Younger US Patients: Implication of Cannabis Exposure as a Risk Factor. Gastroenterology Insights. 2023; 14(2):204-235. https://doi.org/10.3390/gastroent14020016

Chicago/Turabian Style

Reece, Albert Stuart, and Gary Kenneth Hulse. 2023. "Sociodemographically Stratified Exploration of Pancreatic Cancer Incidence in Younger US Patients: Implication of Cannabis Exposure as a Risk Factor" Gastroenterology Insights 14, no. 2: 204-235. https://doi.org/10.3390/gastroent14020016

APA Style

Reece, A. S., & Hulse, G. K. (2023). Sociodemographically Stratified Exploration of Pancreatic Cancer Incidence in Younger US Patients: Implication of Cannabis Exposure as a Risk Factor. Gastroenterology Insights, 14(2), 204-235. https://doi.org/10.3390/gastroent14020016

Article Metrics

Back to TopTop