Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Cannabis use is associated with reduced prevalence of non-alcoholic fatty liver disease: A cross-sectional study

Abstract

Cannabis use is associated with reduced prevalence of obesity and diabetes mellitus (DM) in humans and mouse disease models. Obesity and DM are a well-established independent risk factor for non-alcoholic fatty liver disease (NAFLD), the most prevalent liver disease globally. The effects of cannabis use on NAFLD prevalence in humans remains ill-defined. Our objective is to determine the relationship between cannabis use and the prevalence of NAFLD in humans. We conducted a population-based case-control study of 5,950,391 patients using the 2014 Healthcare Cost and Utilization Project (HCUP), Nationwide Inpatient Survey (NIS) discharge records of patients 18 years and older. After identifying patients with NAFLD (1% of all patients), we next identified three exposure groups: non-cannabis users (98.04%), non-dependent cannabis users (1.74%), and dependent cannabis users (0.22%). We adjusted for potential demographics and patient related confounders and used multivariate logistic regression (SAS 9.4) to determine the odds of developing NAFLD with respects to cannabis use. Our findings revealed that cannabis users (dependent and non-dependent) showed significantly lower NAFLD prevalence compared to non-users (AOR: 0.82[0.76–0.88]; p<0.0001). The prevalence of NAFLD was 15% lower in non-dependent users (AOR: 0.85[0.79–0.92]; p<0.0001) and 52% lower in dependent users (AOR: 0.49[0.36–0.65]; p<0.0001). Among cannabis users, dependent patients had 43% significantly lower prevalence of NAFLD compared to non-dependent patients (AOR: 0.57[0.42–0.77]; p<0.0001). Our observations suggest that cannabis use is associated with lower prevalence of NAFLD in patients. These novel findings suggest additional molecular mechanistic studies to explore the potential role of cannabis use in NAFLD development.

Introduction

Cannabis is the most abused substance which is still illicit in most countries globally [1]. Besides alcohol and tobacco, it is the most widely used drug for recreational purposes [2]. It is currently estimated that between 119–224 million people worldwide [1], and over 22 million individuals in the United States of America (USA), aged 12 years and above use cannabis [2]. With the recent rise in cannabis legalization across different states in the USA, this number is expected to grow even further. While cannabis use (CU) is associated with increased prevalence of disorders such as schizophrenia [3], cyclical hyperemesis syndrome [4], and pulmonary diseases [5], it has also been linked to a reduction in obesity [6] and diabetes mellitus (DM) [79]. The prevalence of these disorders is expected to be impacted as CU increases over the coming years.

Non-Alcoholic Fatty Liver Disease (NAFLD) is a metabolic disorder characterized by excess fat accumulation in the liver and is the most common liver disease in the world [10]. About one-third of US adults have been diagnosed with NAFLD [11,12]. While excessive fat accumulation in the liver has previously been thought to be harmless, studies including ours, have demonstrated that NAFLD can progress to steatohepatitis, fibrosis, cirrhosis and even hepatocellular carcinoma [13,14]. By 2020, it is estimated that progressive NAFLD would become the leading cause of advance liver disease requiring liver transplantation in patients [15]. Obesity and DM are among the known risk factors for NAFLD [16]. Since CU has been associated with reduced prevalence of both obesity [6] and DM [79], we surmised that CU might also have modulatory effects on the prevalence of NAFLD.

Our study was aimed at testing the hypothesis that CU is associated with reduced prevalence of NAFLD given its suppressive effect against obesity and diabetes in humans. Our novel studies revealed that CU was associated with reduced prevalence of NAFLD in patients. However, co-commitment NDCU and dependent abusive alcohol consumption, which can independently induce progressive fatty liver disease, negated the reduced prevalence of NAFLD observed in cannabis-only users.

Materials and methods

Data source and study population

We conducted a population-based case-control study to determine the relationship between CU and NAFLD. We used data from the Healthcare Cost and Utilization Project Nationwide Inpatient Sample (HCUP-NIS) database, which is produced by the Agency for Healthcare Research and Quality (AHRQ) [17]. HCUP-NIS is the largest inpatient database, comprising of over 3000 hospitals sampled from more than 40 states representing over 7 million annual hospital discharges in the U.S.A. For this study, we used the 2014 HCUP-NIS data was collected from January 1st to December 31st, 2014. The data contained 7,071,762 patient discharge records. Each discharge record contained up to 30 different clinical diagnoses using the International Classification of Diseases, Ninth Edition, Clinical Modification (ICD-9-CM). The clinical discharge records also included patient demographics, insurance type, clinical procedures, hospital costs, comorbidities, and hospital characteristics. The HCUP-NIS has been used in numerous studies to provide reliable estimates of conditions within the U.S.A [1822]. Given that the HCUP-NIS is a completely de-identified and publicly available data, no Institutional Review Board approval was required. We evaluated clinical records for patients 18–90 years old (n = 5,950,391).

Characterization of outcomes and exposure variables

Our primary outcome was a diagnosis of NAFLD, defined using ICD-9-CM codes (Table 1). Code 571.8 represents, "other chronic nonalcoholic liver disease". This code is for nonalcoholic liver diseases and includes both Non-Alcoholic Steatohepatitis (NASH) and NAFLD. Recent studies have used this code as an indicator for NAFLD diagnosis [23,24]. Our primary exposure variable was based on current CU, which was stratified into two main groups to investigate dose-dependent CU relationships in NAFLD, namely; dependent CU (DCU, ICD-9-CM codes 304.3x) and non-dependent CU (NDCU, ICD-9-CM codes 305.2x). According to ICD-9-CM, substances coded as CU include patient use of Indian hemp, marijuana and other varieties of cannabis and cannabinoids.

thumbnail
Table 1. ICD-9-CM codes for identifying NAFLD, cannabis use (dependent and non-dependent), and other risk factors.

https://doi.org/10.1371/journal.pone.0176416.t001

We focused on the following NAFLD associated risk factors: demographics (age, gender, race, socioeconomic status, insurance type) and patient specific conditions (hypertension, DM, dyslipidemia, metabolic syndrome, obesity, tobacco). All ICD-9-CM codes used in identifying patient conditions are presented in Table 1. Given increased prevalence of NAFLD among the middle age group, we categorized this group into three sub-groups: <40yrs, 40-60yrs, and >60 years. We additionally categorized racial origins into four sub-groups: Non-Hispanic Whites, Non-Hispanic Blacks, Hispanics, Others (Asians and Native-Americans). Based on the average annual household income of where the patient resides (zip code), we categorized socio-economic status into four groups: lowest quartile, second quartile, third quartile and highest quartile. Patient's gender was either male or female. We classified health insurance types into five sub-categories: Medicare, Medicaid, Private, Self-pay and Others (no charge, other government insurance, worker's compensation, miscellaneous).

Overweight (pre-obese) and obesity diagnosis included patient discharge records with a BMI of 25–29.9 and > = 30 kg/m2 respectively. Tobacco use was defined as having a history or current use of tobacco products; diabetes as having a diagnosis of DM (type I or type II), and alcohol use as having a history or current abusive use of alcohol. We defined hyperlipidemia as having a lipid metabolism disorder including hypercholesterolemia, hyperglyceridemia, hyperchylomicronemia, lipodystrophy, lipoprotein deficiencies, lipidoses, mixed hyperlipidemias and other disorders resulting in elevated plasma lipids. We defined hypertension as having a clinical diagnosis of essential hypertension.

Drinking levels of alcohol is defined clinically following the National Institute on Alcohol Abuse and Alcoholism (NIAAA) and the Substance Abuse and Mental Health Services Administration (SAMHSA) standards following guidelines from the U.S. Department of Health and Human Services and U.S. Department of Agriculture.

Data analysis

All statistical analyses were performed using the Statistical Analysis System (SAS) [Release Version 9.4, SAS Institute Inc, Cary, NC, US]. We characterized patient discharge records with descriptive statistics and reported the mean and standard deviation (SD) for continuous variables and the percentages for categorical variables. We used Chi-squared (X2) test to analyze differences between the sub-groups of categorical variables. We additionally used SURVEYFREQ procedure in SAS, with the STRATA and WEIGHT statements to account for the HCUP-NIS sampling method.

We evaluated the crude's odds ratio (COR) by running a univariate binary logistic regression to assess the association between NAFLD (present/absent) and various risk factors. These factors included CU (DCU and NDCU), age, tobacco, DM, hypertension, obesity, race, dyslipidemia, metabolic syndrome, insurance status and household income. Next, we evaluated the adjusted odds ratio (AOR) by running a multivariate logistic regression between NAFLD and cannabis use for the above-listed risk factors as covariates. For regression analysis, we used the SURVEYLOGISTIC procedure in SAS, with the STRATA and WEIGHT statements to account for the HCUP-NIS sampling method.

Using the final adjusted model described above, we further assessed for effect modification by introducing the interaction term between CU and the each of the other risk factors for NAFLD individually. For the identified risk factors with a significant interaction with CU, we then ran a stratified analysis and measured the effect of CU on NAFLD prevalence.

A pvalue of <0.05 was considered statistically significant for all our statistical tests used in this study.

Results

Study population in relation to cannabis use

Of total of 5,950,391 patient discharge records, aged 18 years and above analyzed, we found that 5,833,812 (98.0%) were non-cannabis users, 103,675 (1.74%) were NDCU and 12,904 (0.22%) were DCU (Table 2).

thumbnail
Table 2. NAFLD, socio-demographic characteristics, and other risk factors among patient discharge records (18 years and above) classified by marijuana use (dependent and non-dependent) and non-use.

https://doi.org/10.1371/journal.pone.0176416.t002

The unadjusted prevalence of NAFLD was lower among CU when compared to non-users (DCU:0.38%, NDCU:0.77% vs. non-users:0.91%). NDU and DCU were more likely to be between the ages 18–40 years and 40–60 years compared to >60 years old. Cannabis users (DCU and NDCU) were more likely to be male, white, have Medicare insurance, and be in the lowest quartile of income. When compared to non-cannabis users, DCU and NDCU were less likely to have hyperlipidemia (8.3%, 12.5% vs. 30.4%) and diabetes (7.9%, 13.6% vs. 25.4%). Though they were slightly more likely to be overweight (0.62%, 0.49% vs. 0.32%), they were however less likely to be obese (7.7%, 9.7% vs. 13.6%) or have metabolic syndrome (0.08%, 0.11% vs. 0.17%). DCU and NDCU were more likely to abuse alcohol (40.1%, 28.8% vs. 5.6%) and tobacco (54.9%, 60.1% vs. 28.2%).

Relationship between NAFLD with cannabis use and other risk factors for NAFLD in humans

The development of NAFLD is associated with numerous risk factors. To determine how CU relates to other risk factors of NAFLD development, we estimated the crude odds ratio (COR) for NAFLD with CU, and similarly with other known risk factors of NAFLD (Table 3). We found that CU was associated with 20% reduced prevalence of NAFLD (AOR: 0.80[0.75–0.87]). When compared to non-users, the prevalence of NAFLD was respectively 68% and 15% less among DCU and NDCU (AOR: 0.42[0.32–0.55] and 0.85[0.79–0.92]). We observed a dose-response effect of cannabis use as DCU had 47% less prevalence of NAFLD when compared to NDCU (AOR 0.53[0.39–0.72]). After adjusting for the potential confounding risk factors listed in Table 3 in our statistical modeling, the association between cannabis use and NAFLD did not change significantly (<10%). These observations suggested that the effect of CU on NAFLD might be mediated through alternative factors other than those evaluated by our modeling. It is worth mentioning that all known factors for NAFLD were included in our final model.

thumbnail
Table 3. Crude and adjusted odds ratio for NAFLD among risk factors for NAFLD among patient.

https://doi.org/10.1371/journal.pone.0176416.t003

All reported estimates were after adjustment for potential confounders (Table 3). When compared to the 18-40-year age group, individuals between the 40–60 and >60-year age groups had 56% higher and 9% lower prevalence of NAFLD (AOR: 1.56[1.51–1.60] and 0.91[0.87–0.94]) respectively. The prevalence of NAFLD was 5% and 26% elevated among females and hypertensives (AOR: 1.05[1.03–1.06] & 1.26[1.23–1.28]), 20% higher among those with hyperlipidemia (AOR: 1.19[1.17–1.22]), 85% higher among DM patients (AOR: 1.85[1.81–1.88]), 120% higher among the pre-obese (AOR: 2.12[1.86–2.41]), 330% higher among the obese (AOR: 3.13[3.02–3.24]), and 170% higher among those with metabolic syndrome (AOR: 2.72[2.49–2.97]). The prevalence of NAFLD was also 5% and 90% higher among people with a past or a current history of tobacco and alcohol use respectively (AOR: 1.05[1.03–1.07] and AOR: 1.93[1.87–1.99]). Compared to non-Hispanic Whites, non-Hispanic Blacks had 45% reduced prevalence of NAFLD (AOR: 0.56[0.55–0.58]), while Hispanics have 30% increased prevalence of NAFLD (AOR: 1.29[1.26–1.33]). There was no statistical difference in the prevalence of NAFLD among Asians/Other races and Whites (AOR: 1.02[0.98–1.06]). Compared to those on Private health insurance, patients on Medicaid, Medicare and others (including uninsured) had 33%, 25%, and 18% less prevalence of NAFLD (AOR: 0.67[0.66–0.69], 0.76[0.74–0.78], and 0.78[0.74–0.82]). Compared to the upper-income bracket, every other lower income quartile had a reduced prevalence of NAFLD. This decrease trended down with household wealth, from the richest to the poorest quartile. When compared to the richest (fourth) quartile, the prevalence of NAFLD was 4%, 8% and 11% less in the third, second and first (poorest) quartile (AOR: 0.94[0.92–0.97], 0.92[0.89–0.95], and 0.89[0.87–0.92]).

Sensitivity and post-hoc analysis

Among all the risk factors analyzed by our study, age, obesity, alcohol use, and hyperlipidemia showed significant statistical interactional association with cannabis in its association with NAFLD (Table 4). We found that the effects of CU on NAFLD was heterogeneous across the age groups (p = 0.0011). NDCU was only associated with a reduced prevalence of NAFLD in individuals 40–60 years (Interaction AOR: 0.77 (0.72–0.84), p<0.0012). DCU was associated with a 2 to 3-fold reduced prevalence of NAFLD among individuals 18–40 and 40–60 years but not among those >60 years, where there was no difference compared to non-users (Interaction AOR: 3.1 (1.77–5.43), p = 0.043). Obesity was associated with about 2 to 3-fold elevation in the prevalence of NAFLD so that CU (DCU and NCDU) was no longer associated with a lower prevalence of NAFLD among the obese (Interaction OR: 2.3(1.68–3.10), p = 0.0073).

thumbnail
Table 4. Modification of the effect of cannabis on NAFLD by age, obesity, alcohol use, hyperlipidemia.

https://doi.org/10.1371/journal.pone.0176416.t004

Alcohol use was associated with a 90% increase in NAFLD prevalence (Table 3, AOR: 1.93(1.87–1.99)). Alcohol use, however, showed a differential modulatory effect on NAFLD prevalence amongst CU. Among both alcohol consumers, DCU was associated with a lower prevalence of NAFLD. Among NDCU, there was a lower prevalence of NAFLD in alcoholics compared to non-alcoholics (Interaction AOR: 0.72 (0.66–0.78), p<0.0001). To further explore the effects of CU-alcohol use interactions, we categorized alcohol use into two-subgroups: non-dependent users, and dependent users, to respectively constitute moderate and severe consumption. Statistical analysis with these variable adjustments revealed that NDCU was associated with a 53% lower prevalence of NAFLD among non-dependent alcohol consumers (AOR: 0.57(0.48–0.68)) but not among dependent (AOR: 0.95(0.87–1.05)). DCU, however, was associated with a 70% and 55% reduction in NAFLD prevalence among non-dependent and dependent consumers of alcohol (AOR: 0.30(0.11–0.79) & 0.45(0.28–0.70)).

Given that NAFLD patients who take alcohol could also develop Alcoholic Steatohepatitis (ASH), we completely excluded patients with any past or current history of alcoholic use (n = 358,507) from our study, and re-estimated the crude and adjusted OR for developing NAFLD among the remaining subjects (n = 5,591,884). NDCU and DCU were respectively associated with 21% and 66% reduced prevalence for NAFLD (COR: 0.79(0.73–0.87) & 0.34(0.23–0.52)) [Table 5]. After accounting for other risk factors of NAFLD, NDCU was only associated with 2% reduced prevalence and was no longer statistically significant, unlike DCU which was associated with a significant 48% reduced prevalence of NAFLD (AOR: 0.52(0.34–0.80)).

thumbnail
Table 5. Crude and adjusted odds ratio for NAFLD after excluding patients with any past or present history of alcohol use.

https://doi.org/10.1371/journal.pone.0176416.t005

Discussion

Our study explored the 2014 NIS data set to evaluate the association between CU and NAFLD among discharged patients aged 18 years and above in the USA. We observed a strong dose-dependent reduction in the prevalence of NAFLD with cannabis use (Fig 1) suggesting that cannabis use might suppress or reverse NAFLD development. Additionally, increased prevalence of NAFLD was associated with known factors including age 40–60 years, female gender, hyperlipidemia, hypertension, alcohol, diabetes, and metabolic syndrome. We observed racial differences in NAFLD, whereby non-Hispanic Whites and Hispanics had higher odds for disease compared to non-Hispanic Blacks. Patients with public health insurance (Medicare and Medicaid) had a lower prevalence of NAFLD compared to privately insured and uninsured individuals (Health insurance and type of insurance relate broadly to access and quality of health care provided to patients respectively). This novel revelations contrast commonly held notions and suggest that individuals on government health insurance might have a reduced prevalence of a chronic disorder such as NAFLD due to improved emphasis on preventive services since the implementation of the Affordable Care Act [25,26].

thumbnail
Fig 1. Cannabis use is associated with reduced non-alcoholic fatty liver disease.

Dependent and Non-Dependent Cannabis Use (DCU & NDCU) are associated with reduced NAFLD when compared to non-cannabis users (NCU). These observations suggest that dependent cannabis use suppresses or reverses the development and progression of NAFLD to advance liver disease (non-alcoholic steatohepatitis [NASH]). Illustrated schematics made use of some motifolio templates (www.motifolio.com).

https://doi.org/10.1371/journal.pone.0176416.g001

We found that alcohol use was associated with a higher prevalence of NAFLD. Though CU reduced the prevalence of NAFLD in patients, only non-dependent alcohol use maintained some of this beneficial health effect (Fig 2). NDCU only had lower odds for NAFLD among non-dependent alcohol individuals, unlike DCU which maintained lower odds for NAFLD irrespective of the volume of alcohol use (non-dependent & dependent) (Fig 2). Studies have reported that DCU more likely consume cannabis at a larger quantity and frequency than NDCU [27,28]. Therefore, with higher consumption of cannabis, DCU might have a higher dose and longer duration of cannabis acting on the hepatic cannabinoid receptors, and offering a protective effect when compared to NDCU. A recent study has indeed shown that cannabis act directly on hepatic cells and protects against liver disease [29]. Taken together, these observations suggest that chronic dependent alcohol negated the beneficial effects of NDCU in modestly reducing the prevalence of NAFLD in patients. Our CU-non-dependent alcohol user interaction observations suggests that moderate alcohol intake, irrespective of the volume of cannabis used, might have beneficial effects in reducing the prevalence of NAFLD. Our novel observations are somehow in concert with published findings which demonstrated that moderate alcohol consumption provided some beneficial effects in NAFLD patients [28].

thumbnail
Fig 2. Dependent alcohol use abolishes reduced NAFLD prevalence observed in non-dependent cannabis users.

Bar graph described from the leftmost to the rightmost: Taking non-cannabis users/non-alcohol users (NCU+NAU) as the reference group, non-dependent (moderate) cannabis use/non-alcohol use (NDCU+NAU) caused a slight reduction in NAFLD prevalence, though not statistically significant. Non-dependent (moderate) alcohol use among NDCU (NDCU+NDA) resulted in a reduced NAFLD prevalence, though this protection was lost with dependent (abusive) alcohol use (NDCU+DA). However dependent (high volume) cannabis use (DCU) was always associated with a reduction in NAFLD prevalence among all categories of alcohol users: non-alcohol (DCU+NAU), moderate/non-dependent alcohol users (DCU+NDA) and copious/dependent alcohol users (DCU+DA). Illustrated schematics made use of some motifolio templates (www.motifolio.com).

https://doi.org/10.1371/journal.pone.0176416.g002

CU is negatively associated with many risk factors of NAFLD characterized by reduced prevalence of DM [8,9,29,30], with overweight/obesity [6], and hyperlipidemia [7]. Cannabis might act on fatty deposits in obese tissues through omega-3 (n-3) fatty acids [31]. For example, preclinical trials have revealed that endogenously produced cannabinoids mediate the ability of (n-3) fatty acids in reducing ectopic fats deposits and fatty inflammation [30]. A meta-analysis in humans further confirms that dietary (n-3) supplementation resulted in a reduction in steatosis and improved liver enzyme profile [31]. Another study demonstrated that CU resulted in reduced fasting insulin, insulin resistance, smaller waist circumference and higher levels of HDL [32]. Taken together, these observations suggest that cannabis use has modulatory effects on NAFLD risk factors.

Given that the effects of CU on NAFLD persisted even after accounting for the other risk factors for NAFLD and the proposed mediators (obesity, hyperlipidemia, DM) in our statistical modeling, we surmise that cannabis probably acts on the liver via alternate mediators and mechanisms. In addition to other ingredients, cannabis contains tetrahydrocannabinol (THC) and cannabidiol (CBD) [33,34]. THC and CBD act on CB-1 and CB-2 receptors [35]. While THC has more effect on the CB-1, CBD act preferentially on CB-2. CB-1 has been shown to be more pro-inflammatory and psychotropic, unlike CB-2, which is more anti-inflammatory and non-psychotropic [35,36]. Pre-clinical studies have demonstrated that hepatic CB-1 receptors up-regulate hepatic fatty acid production and participate in diet-induced obesity [37]. Hepatic CB-1 receptors also promote liver fibrosis in both murine and human subjects [38]. In contrast, CB-2 agonists can suppress obesity, steatohepatitis and can additionally protect the liver from ischemic reperfusion injury. Further, CB-2 agonist decrease hepatic immune cells infiltration, kupffer cell activation, oxidative stress, hepatic injury and fibrogenesis [3942]. Taken together, our novel observations and previous findings suggest that CB-2 agonist and hemp which has a higher CBD might offer more protection from NAFLD [43,44].

The association between age and NAFLD remains controversial based on published literature. While some studies showed that the prevalence of NAFLD rises with age [4548], many have reported otherwise [4952]. It is viewed that since the risk factors for NAFLD including DM, obesity and hyperlipidemia all increase with age, NAFLD prevalence should equally follow this trend. Given that numerous studies have however demonstrated that NAFLD prevalence and age follows an "inverted U shaped curve" [53], we reason that this could be due to two factors: One could be that with increasing age, there is relatively more visceral fat deposition [54] and reduced subcutaneous and hepatic fat deposits [55]. This argument might however be faulted because visceral fat correlates with hepatic fat and NAFLD [56]. Another reason could be because individuals with NAFLD also have a higher cardio-metabolic risks, they may die early from other associated diseases, resulting in fewer elderly patients with NAFLD in the study population[57,58].

Strikingly, low-income patients had the lowest prevalence of NAFLD compared to patients in the higher income bracket. Being in a lower median household income corresponds to having a lower socioeconomic status (SES) or poverty. It is well established that low SES is associated with obesity, hypertension, diabetes, hyperlipidemia [5961], which are all predisposing factors for NAFLD. It would be expected that the association of these conditions with low SES would result in a higher NAFLD among the low income population. We did additional posthoc analysis on our data, where we built statistical models for each of the aforementioned four conditions. Besides hyperlipidemia were the odds for NAFLD was reduced and further decreased with reducing income, the odds for obesity, hypertension and diabetes was higher and trended upwards as the median income level reduced (S1 Table). Despite three of the four risk factors for NAFLD associated with decreasing income, the odds for NAFLD was lower with reducing income both in our crude and adjusted models. The reasons for this new observation remain unclear, but suggest that hyperlipidemia might be a more important predisposing risk factor for NAFLD development. A pediatric study revealed no association between maternal poverty levels and developing NAFLD in children [62]. We suspect that the slight reduction in odds for NAFLD among individuals of low socioeconomic status might be due to a higher non-exercise related physical activity from engaging in more manual occupations. Given that the liver is responsible for energy (gluconeogenisis and lipogenesis) production and storage in high fed states, and the body first depletes hepatic energy stores before resulting to extra-hepatic reserves during manual activities. Individuals with more blue-collar jobs might be more likely to expend the hepatic stores and prevent/reduce long term hepatic fat accumulation. A recent study demonstrated that NAFLD subjects undergoing any form of physical activity experience significant reduction in their hepatic fat deposits even without any change in total body weight [6365]. A more comprehensive meta-analysis revealed that different forms of exercise, ranging from low-impact to vigorous exertion, might all be beneficial for NAFLD [6668]. More prospective studies are however needed to determine the reasons for this observation.

Strengths and limitations

Our study is limited by its cross-sectional study design [69]. These include recall bias in reporting exposures, not well-established sensitivity and specificity of ICD-9-CM codes for cannabis use. Indeed, the ICD-9-CM code does not differentiate the particular type of cannabis plant used. We, however, expect the recall bias, and less specificity and sensitivity of ICD-9-CM cannabis codes to have a similar chance of occurring among those with NAFLD as well as those without NAFLD. With this imprecision in the coding, our results would likely be biased towards the null, that is no association, and the true effect of CU on NAFLD might be underestimated.

The lack of temporal observations in cross-sectional studies does inhibit our ability to draw direct causal relationships between CU and NAFLD. Failure to measure and or record other possible confounding factors might limit our inferences on the effects of CU and NAFLD. To strengthen our observations, we reviewed recent literature to identify additional risk factor for NAFLD which were not retrievable from the NIS dataset [16].

Different strains of the cannabis plant also contain different ratios of anti-inflammatory CBD and pro-psychotic/pro-inflammatory THC agents [34,36,43,44]. For example, hemp, cannabis oil, hash oil, sativa, and indica all have different contents of CBD and THC and the various contents of these active ingredients could not be assessed from the data set analyzed.

Despite the shortcomings of cross-sectional studies, our novel observations are necessary given the increased legalization of cannabis in the USA and increased prevalence of NAFLD globally. The fact that we used the NIS data archives, which is a very reliable large population-based database, allows for under-reported outcomes, such as cannabis use and NAFLD with good statistical power.

Conclusions

To the best of our knowledge, this is the first population-based cross-sectional study of hospitalized patients to explore the associations between CU and NAFLD. Our analyses revealed a strong relationship between cannabis use and reduced prevalence of NAFLD in patients. Due to our inability to draw direct causation effects from our cross-sectional studies, we suggest prospective basic and human studies to decipher the mechanistic details of how the various active ingredients in cannabis modulate NAFLD development.

Supporting information

S1 Table. Adjusted odds ratio for having a diagnosis of NAFLD among different income levels with obesity, hypertension, diabetes and hyperdilipemia.

*p-value<0.0001; ** Annual income stratified by residence zip-code, 1st quartile:$1-$39,999, 2nd quartile:$40,000-$50,999, 3rd quartile:$51,000-$65,999, 4th quartile:$66,000+.

https://doi.org/10.1371/journal.pone.0176416.s001

(DOCX)

Acknowledgments

This work was funded by a start-up grant to TNB by INRS-Institut Armand-Frappier, Institut National de la Reserche Scientifique, 531 Boulevard des Prairies, Laval (Quebec) H7V 1B7 Canada. The funders had no role in the project design, execution, data interpretation or decision to publish.

Author Contributions

  1. Conceptualization: ACA TNB.
  2. Data curation: ACA SA TOA TNB.
  3. Formal analysis: ACA SA TOA KA TNB.
  4. Funding acquisition: TNB.
  5. Investigation: ACA SA TOA KA NO OMA MD AA.
  6. Methodology: ACA SA TNB.
  7. Project administration: TNB.
  8. Resources: ACA TNB.
  9. Software: ACA TNB.
  10. Supervision: TNB.
  11. Validation: ACA SA TOA KA TNB.
  12. Writing – original draft: ACA TNB.
  13. Writing – review & editing: ACA SA TOA OMA KA MD TNB.

References

  1. 1. United Nations Office on Drugs and Crime. World drug report 2012. [Internet]. New York: United Nations; 2012. http://search.ebscohost.com/login.aspx?direct=true&scope=site&db=nlebk&db=nlabk&AN=464981
  2. 2. Substance Abuse and Mental Health Services Administration. Results from the 2014 National Survey on Drug Use and Health: National Findings. HHS Publication No. (SMA) 14–4863. Center for Behavioral Health Statistics and Quality, Rockville, MD:.;
  3. 3. Müller-Vahl KR, Emrich HM. Cannabis and schizophrenia: towards a cannabinoid hypothesis of schizophrenia. Expert Rev Neurother. 2008;8: 1037–1048. pmid:18590475
  4. 4. Allen JH. Cannabinoid hyperemesis: cyclical hyperemesis in association with chronic cannabis abuse. Gut. 2004;53: 1566–1570. pmid:15479672
  5. 5. Martinasek MP, McGrogan JB, Maysonet A. A Systematic Review of the Respiratory Effects of Inhalational Marijuana. Respir Care. 2016;61: 1543–1551. pmid:27507173
  6. 6. Le Strat Y, Le Foll B. Obesity and Cannabis Use: Results From 2 Representative National Surveys. Am J Epidemiol. 2011;174: 929–933. pmid:21868374
  7. 7. Rajavashisth TB, Shaheen M, Norris KC, Pan D, Sinha SK, Ortega J, et al. Decreased prevalence of diabetes in marijuana users: cross-sectional data from the National Health and Nutrition Examination Survey (NHANES) III. BMJ Open. 2012;2: e000494. pmid:22368296
  8. 8. Sidney S. Marijuana Use and Type 2 Diabetes Mellitus: a Review. Curr Diab Rep. 2016;16.
  9. 9. Dagogo-Jack S. Diabetes risks from prescription and non-prescription drugs: mechanisms and approaches to risk reduction. Alexandria: American Diabetes Association; 2016.
  10. 10. Vernon G, Baranova A, Younossi ZM. Systematic review: the epidemiology and natural history of non-alcoholic fatty liver disease and non-alcoholic steatohepatitis in adults. Aliment Pharmacol Ther. 2011;34: 274–285. pmid:21623852
  11. 11. Loomba R, Sanyal AJ. The global NAFLD epidemic. Nat Rev Gastroenterol Hepatol. 2013;10: 686–690. pmid:24042449
  12. 12. Bellentani S, Scaglioni F, Marino M, Bedogni G. Epidemiology of Non-Alcoholic Fatty Liver Disease. Dig Dis. 2010;28: 155–161. pmid:20460905
  13. 13. Ganz M, Bukong TN, Csak T, Saha B, Park J-K, Ambade A, et al. Progression of non-alcoholic steatosis to steatohepatitis and fibrosis parallels cumulative accumulation of danger signals that promote inflammation and liver tumors in a high fat–cholesterol–sugar diet model in mice. J Transl Med. 2015;13: 193. pmid:26077675
  14. 14. Nonalcoholic Fatty Liver Disease: Pathology and Pathogenesis. Annu Rev Pathol Mech Dis. 2010;5: 145–171.
  15. 15. Charlton M. Cirrhosis and liver failure in nonalcoholic fatty liver disease: Molehill or mountain? Hepatology. 2008;47: 1431–1433. pmid:18393323
  16. 16. Bellentani S. The epidemiology of non-alcoholic fatty liver disease. Liver Int. 2017;37: 81–84. pmid:28052624
  17. 17. Citations for HCUP Databases and Tools [Internet]. [cited 6 Feb 2017]. https://www.hcup-us.ahrq.gov/tech_assist/citations.jsp
  18. 18. HCUP-US NIS Overview [Internet]. [cited 6 Feb 2017]. https://www.hcup-us.ahrq.gov/nisoverview.jsp
  19. 19. Houchen R, Ross D, Elixhauser A, Jiang J. Nationwide Inpatient Sample (NIS) Redesign Final Report. 2014. HCUP Methods Series Report # 2014–04 ONLINE [Internet]. U.S. Agency for Healthcare Research and Quality; https://www.hcup-us.ahrq.gov/reports/methods/methods.jsp
  20. 20. Olsen MA, Young-Xu Y, Stwalley D, Kelly CP, Gerding DN, Saeed MJ, et al. The burden of clostridium difficile infection: estimates of the incidence of CDI from U.S. Administrative databases. BMC Infect Dis. 2016;16.
  21. 21. Schoenman JA, Sutton JP, Elixhauser A, Love D. Understanding and Enhancing the Value of Hospital Discharge Data. Med Care Res Rev. 2007;64: 449–468. pmid:17684112
  22. 22. Okafor PN, Chiejina M, de Pretis N, Talwalkar JA. Secondary analysis of large databases for hepatology research. J Hepatol. 2016;64: 946–956. pmid:26739689
  23. 23. Loomis AK, Kabadi S, Preiss D, Hyde C, Bonato V, St. Louis M, et al. Body Mass Index and Risk of Nonalcoholic Fatty Liver Disease: Two Electronic Health Record Prospective Studies. J Clin Endocrinol Metab. 2016;101: 945–952. pmid:26672639
  24. 24. Cohen D, Gonzales-Pacheco D, Myers O. Relationships Between Alanine Aminotransferase, Serum Triglycerides, Body Mass Index and Nonalcoholic Fatty Liver Disease in an Outpatient Pediatric Clinic Population. J Pediatr Nurs. 2016;31: 152–158. pmid:26690717
  25. 25. 04 A, 2015. Preventive Services Covered by Private Health Plans under the Affordable Care Act [Internet]. [cited 9 Mar 2017]. http://kff.org/health-reform/fact-sheet/preventive-services-covered-by-private-health-plans/
  26. 26. Shaw FE, Asomugha CN, Conway PH, Rein AS. The Patient Protection and Affordable Care Act: opportunities for prevention and public health. The Lancet. 2014;384: 75–82.
  27. 27. Chen K, Kandel DB, Davies M. Relationships between frequency and quantity of marijuana use and last year proxy dependence among adolescents and adults in the United States1. Drug Alcohol Depend. 1997;46: 53–67. pmid:9246553
  28. 28. Swift W, Coffey C, Carlin JB, Degenhardt L, Calabria B, Patton GC. Are adolescents who moderate their cannabis use at lower risk of later regular and dependent cannabis use? Addiction. 2009;104: 806–814. pmid:19344439
  29. 29. Trebicka J, Racz I, Siegmund SV, Cara E, Granzow M, Schierwagen R, et al. Role of cannabinoid receptors in alcoholic hepatic injury: steatosis and fibrogenesis are increased in CB2 receptor-deficient mice and decreased in CB1 receptor knockouts. Liver Int Off J Int Assoc Study Liver. 2011;31: 860–870.
  30. 30. Batetta B, Griinari M, Carta G, Murru E, Ligresti A, Cordeddu L, et al. Endocannabinoids May Mediate the Ability of (n-3) Fatty Acids to Reduce Ectopic Fat and Inflammatory Mediators in Obese Zucker Rats. J Nutr. 2009;139: 1495–1501. pmid:19549757
  31. 31. Parker HM, Johnson NA, Burdon CA, Cohn JS, O’Connor HT, George J. Omega-3 supplementation and non-alcoholic fatty liver disease: A systematic review and meta-analysis. J Hepatol. 2012;56: 944–951. pmid:22023985
  32. 32. Penner EA, Buettner H, Mittleman MA. The Impact of Marijuana Use on Glucose, Insulin, and Insulin Resistance among US Adults. Am J Med. 2013;126: 583–589. pmid:23684393
  33. 33. McPartland JM, Russo EB. Cannabis and Cannabis Extracts: Greater Than the Sum of Their Parts? J Cannabis Ther. 2001;1: 103–132.
  34. 34. Andre CM, Hausman J-F, Guerriero G. Cannabis sativa: The Plant of the Thousand and One Molecules. Front Plant Sci. 2016;7.
  35. 35. Burstein S. Cannabidiol (CBD) and its analogs: a review of their effects on inflammation. Bioorg Med Chem. 2015;23: 1377–1385. pmid:25703248
  36. 36. Piomelli D, Russo EB. The Cannabis sativa Versus Cannabis indica Debate: An Interview with Ethan Russo, MD. Cannabis Cannabinoid Res. 2016;1: 44–46.
  37. 37. Osei-Hyiaman D, DePetrillo M, Pacher P, Liu J, Radaeva S, Bátkai S, et al. Endocannabinoid activation at hepatic CB1 receptors stimulates fatty acid synthesis and contributes to diet-induced obesity. J Clin Invest. 2005;115: 1298–1305. pmid:15864349
  38. 38. Teixeira-Clerc F, Julien B, Grenard P, Van Nhieu JT, Deveaux V, Li L, et al. CB1 cannabinoid receptor antagonism: a new strategy for the treatment of liver fibrosis. Nat Med. 2006;12: 671–676. pmid:16715087
  39. 39. Rossi F, Bellini G, Luongo L, Manzo I, Tolone S, Tortora C, et al. Cannabinoid Receptor 2 as Antiobesity Target: Inflammation, Fat Storage, and Browning Modulation. J Clin Endocrinol Metab. 2016;101: 3469–3478. pmid:27294325
  40. 40. Batkai S, Osei-Hyiaman D, Pan H, El-Assal O, Rajesh M, Mukhopadhyay P, et al. Cannabinoid-2 receptor mediates protection against hepatic ischemia/reperfusion injury. FASEB J. 2007;21: 1788–1800. pmid:17327359
  41. 41. Patsenker E, Stickel F. Cannabinoids in liver diseases: Endocannabinoids and Liver. Clin Liver Dis. 2016;7: 21–25.
  42. 42. Julien B, Grenard P, Teixeira-Clerc F, Van Nhieu JT, Li L, Karsak M, et al. Antifibrogenic role of the cannabinoid receptor CB2 in the liver. Gastroenterology. 2005;128: 742–755. pmid:15765409
  43. 43. Giroud C. Analysis of Cannabinoids in Hemp Plants. Chim Int J Chem. 2002;56: 80–83.
  44. 44. Brown DT. Cannabis: the genus Cannabis [Internet]. 1998. http://www.crcnetbase.com/isbn/9789057022913
  45. 45. Kojima S-I, Watanabe N, Numata M, Ogawa T, Matsuzaki S. Increase in the prevalence of fatty liver in Japan over the past 12 years: analysis of clinical background. J Gastroenterol. 2003;38: 954–961. pmid:14614602
  46. 46. Eguchi Y, Hyogo H, Ono M, Mizuta T, Ono N, Fujimoto K, et al. Prevalence and associated metabolic factors of nonalcoholic fatty liver disease in the general population from 2009 to 2010 in Japan: a multicenter large retrospective study. J Gastroenterol. 2012;47: 586–595. pmid:22328022
  47. 47. Shen L, Fan J-G, Shao Y, Zeng M-D, Wang J-R, Luo G-H, et al. Prevalence of nonalcoholic fatty liver among administrative officers in Shanghai: an epidemiological survey. World J Gastroenterol. 2003;9: 1106–1110. pmid:12717867
  48. 48. Foster T, Anania FA, Li D, Katz R, Budoff M. The prevalence and clinical correlates of nonalcoholic fatty liver disease (NAFLD) in African Americans: the multiethnic study of atherosclerosis (MESA). Dig Dis Sci. 2013;58: 2392–2398. pmid:23546700
  49. 49. Lonardo A, Lombardini S, Scaglioni F, Ballestri S, Verrone A-M, Bertolotti M, et al. Fatty liver, carotid disease and gallstones: a study of age-related associations. World J Gastroenterol. 2006;12: 5826–5833. pmid:17007049
  50. 50. Chen C-H, Huang M-H, Yang J-C, Nien C-K, Yang C-C, Yeh Y-H, et al. Prevalence and risk factors of nonalcoholic fatty liver disease in an adult population of taiwan: metabolic significance of nonalcoholic fatty liver disease in nonobese adults. J Clin Gastroenterol. 2006;40: 745–752. pmid:16940890
  51. 51. Younossi ZM, Stepanova M, Negro F, Hallaji S, Younossi Y, Lam B, et al. Nonalcoholic fatty liver disease in lean individuals in the United States. Medicine (Baltimore). 2012;91: 319–327.
  52. 52. Nomura H, Kashiwagi S, Hayashi J, Kajiyama W, Tani S, Goto M. Prevalence of fatty liver in a general population of Okinawa, Japan. Jpn J Med. 1988;27: 142–149. pmid:3047469
  53. 53. Bertolotti M, Lonardo A, Mussi C, Baldelli E, Pellegrini E, Ballestri S, et al. Nonalcoholic fatty liver disease and aging: Epidemiology to management. World J Gastroenterol WJG. 2014;20: 14185–14204. pmid:25339806
  54. 54. Hunter GR, Gower BA, Kane BL. Age Related Shift in Visceral Fat. Int J Body Compos Res. 2010;8: 103–108. pmid:24834015
  55. 55. N T, N T, Y J, S A, I Y, Y T. Age is a Negative, and Visceral Fat Accumulation is a Positive, Contributor to Hepatic Steatosis, Regardless of the Fibrosis Progression in Non-alcoholic Fatty Liver Disease. J Gastroenterol Hepatol Res. 2012;1: 315–319.
  56. 56. Eguchi Y, Eguchi T, Mizuta T, Ide Y, Yasutake T, Iwakiri R, et al. Visceral fat accumulation and insulin resistance are important factors in nonalcoholic fatty liver disease. J Gastroenterol. 2006;41: 462–469. pmid:16799888
  57. 57. Azzam H, Malnick S. Non-alcoholic fatty liver disease—the heart of the matter. World J Hepatol. 2015;7: 1369–1376. pmid:26052382
  58. 58. Lonardo A, Sookoian S, Chonchol M, Loria P, Targher G. Cardiovascular and systemic risk in nonalcoholic fatty liver disease—atherosclerosis as a major player in the natural course of NAFLD. Curr Pharm Des. 2013;19: 5177–5192. pmid:23432668
  59. 59. Hruby A, Hu FB. The Epidemiology of Obesity: A Big Picture. PharmacoEconomics. 2015;33: 673–689. pmid:25471927
  60. 60. Levine JA. Poverty and Obesity in the U.S. Diabetes. 2011;60: 2667–2668. pmid:22025771
  61. 61. Bird Y, Lemstra M, Rogers M, Moraros J. The relationship between socioeconomic status/income and prevalence of diabetes and associated conditions: A cross-sectional population-based study in Saskatchewan, Canada. Int J Equity Health. 2015;14: 93. pmid:26458543
  62. 62. Katz K, Brar PC, Parekh N, Liu Y-H, Weitzman M. Suspected Nonalcoholic Fatty Liver Disease Is Not Associated with Vitamin D Status in Adolescents after Adjustment for Obesity. J Obes. 2010;2010.
  63. 63. Keating SE, Hackett DA, Parker HM, O’Connor HT, Gerofi JA, Sainsbury A, et al. Effect of aerobic exercise training dose on liver fat and visceral adiposity. J Hepatol. 2015;63: 174–182. pmid:25863524
  64. 64. Loomba R, Cortez-Pinto H. Exercise and improvement of NAFLD: Practical recommendations. J Hepatol. 2015;63: 10–12. pmid:25863525
  65. 65. Hallsworth K, Fattakhova G, Hollingsworth KG, Thoma C, Moore S, Taylor R, et al. Resistance exercise reduces liver fat and its mediators in non-alcoholic fatty liver disease independent of weight loss. Gut. 2011;60: 1278–1283. pmid:21708823
  66. 66. Oh S, Shida T, Sawai A, Maruyama T, Eguchi K, Isobe T, et al. Acceleration training for managing nonalcoholic fatty liver disease: a pilot study. Ther Clin Risk Manag. 2014;10: 925–936. pmid:25404857
  67. 67. Whitsett M, VanWagner LB. Physical activity as a treatment of non-alcoholic fatty liver disease: A systematic review. World J Hepatol. 2015;7: 2041–2052. pmid:26261693
  68. 68. Kawaguchi T, Shiba N, Maeda T, Matsugaki T, Takano Y, Itou M, et al. Hybrid training of voluntary and electrical muscle contractions reduces steatosis, insulin resistance, and IL-6 levels in patients with NAFLD: a pilot study. J Gastroenterol. 2011;46: 746–757. pmid:21340530
  69. 69. Levin KA. Study design III: Cross-sectional studies. Evid Based Dent. 2006;7: 24–25. pmid:16557257