Abstract
Cannabis is the most widely used substance in the world. This study aimed to estimate the prevalence of cannabis lifetime use (CLU) in high school and college students of Iran and also to determine factors related to changes in prevalence. A systematic review of literature on cannabis use in Iran was conducted according to MOOSE guideline. Domestic scientific databases, PubMed/Medline, ISI Web of Knowledge, and Google Scholar, relevant reference lists, and relevant journals were searched up to April, 2014. Prevalences were calculated using the variance stabilizing double arcsine transformation and confidence intervals (CIs) estimated using the Wilson method. Heterogeneity was assessed by Cochran’s Q statistic and I2 index and causes of heterogeneity were evaluated using meta-regression model. In electronic database search, 4,000 citations were retrieved, producing a total of 33 studies. CLU was reported with a random effects pooled prevalence of 4.0% (95% CI = 3.0% to 5.0%). In subgroups of high school and college students, prevalences were 5.0% (95% CI = 3.0% to −7.0%) and 2.0% (95% CI = 2.0% to −3.0%), respectively. Meta-regression model indicated that prevalence is higher in college students (β = 0.089, p < .001), male gender (β = 0.017, p < .001), and is lower in studies with sampling versus census studies (β = −0.096, p < .001). This study reported that prevalence of CLU in Iranian students are lower than industrialized countries. In addition, gender, level of education, and methods of sampling are highly associated with changes in the prevalence of CLU across provinces.
Introduction
Cannabis (marijuana) is the most widely used substance in the world (Asbridge, Hayden, & Cartwright, 2012). Cannabis has been known to Homo sapiens in history and concern has been raised during the past centuries regarding its side effects (Anthony, 2012). Several serious side effects are enumerated for cannabis use such as suppression of rapid eye movement and diffuse slowing of background electroencephalographic activity (Greydanus, Hawver, Greydanus, & Merrick, 2013), chronic cough, sinusitis, conjunctivitis, pharyngitis, and bronchitis (Mallaret, Dal’Bo-Rohrer, & Dematteis, 2005). Increased mortality related to its side effects is reported in countries where the prevalence of cannabis use is high (Greydanus et al., 2013). Also, cannabis use has been associated with an increased risk of deadly motor vehicle crashes (Asbridge et al., 2012).
According to United Nations Office on Drugs and Crime, the majority of consumers of cannabis were young and the prevalence has been increasing in many countries (Meimandi, Nakhaee, Divsalar, & Dabiri, 2005). Meanwhile, the geographic distribution is very diverse. In Armenia, it is estimated that the prevalence is about 3.0% and in Czech Republic, 45.0% (Greydanus et al., 2013). The distribution in countries of Middle East and western part of Asia is not well studied.
It is established that two of the major cannabis production and distribution routes are in Pakistan and Afghanistan, implicating major problems for the countries around, especially for Iran (United Nations Office on Drugs and Crime, 2003). Studies in Iran have reported different prevalences across the country. In the north and northwest of the country, a lifetime prevalence of 22.2% (Mohammad Poorasl, Vahidi, Fakhari, Rostami, & Dastghiri, 2007) and a point prevalence of 3.6% (Najafi et al., 2007) have been reported.
No systematic review or meta-analyses have been conducted in Iran to estimate a combined prevalence of cannabis use according to different administrative divisions. In addition, the causes of changes in the prevalence across provinces in a country have not been studied. The aim of this study is to systematically review all cross-sectional studies conducted during the past 30 years in order to estimate cannabis lifetime use (CLU) pooled prevalence in high school and college students of Iran as well as determining factors related to changes in prevalence across provinces using meta-regression model.
Materials and Method
Search Strategy
The search was conducted for cross-sectional studies or surveys that estimated the prevalence of CLU in each province of Iran that was published since 1979. All domestic scientific databases including Iranmedex, Scientific Information Database, Magiran, Irandoc, Medlib, and IranPsych, as well as international databases, including PubMed/Medline, ISI Web of Knowledge, and Google Scholar, relevant reference lists, and relevant journals were searched up to April 2014. Communication with authors and Iranian experts was done for any additional information or unpublished studies.
Two researchers independently scanned the titles of all retrieved citations, removed duplicates, and identified potentially relevant articles for inclusion. Abstracts from selected articles were then independently reviewed by two researchers for further relevance, and full text manuscripts retrieved when appropriate. In the case of disagreement, a third assessor acted as a mediator.
The search strategy was limited to the Persian and/or English. For the databases search, the combinations of keywords relating to CLU were used. These include the following: “Marijuana,” “Marijuana Smoking,” “Marijuana abuse,” “Hashish,” “Cannabis,” “Substance,” “Drug,” “Illicit drug,” “Adolescent,” “Youth,” “School students,” “high school students,” “student,” and “Iran.” The search strategy in PubMed was (Marijuana OR “Marijuana Smoking” OR “Marijuana abuse” OR Hashish OR Cannabis OR Substance OR Drug OR “Illicit drug”) AND (Adolescent OR Youth OR “School students”) AND (Students) AND (Iran). Also, this strategy was modified for each database. EndNote X4 software was employed to manage and scan citations.
Inclusion and Exclusion Criteria and Quality Assessment
The following exclusion criteria were applied: (a) irrelevant study design such as case-control or clinical trials; (b) inadequate reporting of results, that is, studies not reporting prevalence for CLU; (c) poor quality score (scores below five stars) based on the Newcastle−Ottawa scale (Wells et al., 2014; This scale assigns a maximum of ten stars to each study—five stars for representativeness of the sample, justification of sample size, nonrespondents rate, and ascertainment of the exposure; two stars for comparability [control for confounder and additional factors]; and three stars for ascertainment of the outcomes); (d) articles that were conducted before 1979, and (e) review articles or meta-analyses.
For inclusion into the systematic review, the studies had to meet the following criteria: (a) having clear definition of CLU; (b) having standardized diagnostic criteria for drug abuse/use using the International Classification of Diseases, Classification of Mental and Behavioral Disorders, or Diagnostic and Statistical Manual of Mental Disorders, or direct question in questionnaire and interview; (c) having prevalence proportion for CLU; and (d) place of study should be in Iran.
Data Extraction and Measures
Information on first author, year of study, method of sampling, sample size, average age, diagnostic instrument, diagnostic criteria, participation rate, geographical location, level of education, and prevalence during lifetime, 1 month later, 1 week later, and current use were extracted independently from every eligible study.
Average age categorized as ≥18 versus ≤18 years; level of education: college students and high school students; gender: male, female, and mixed; sampling methods: census, multistage, stratified, random, and haphazard sampling; geographical location: name of provinces. Other variables were considered as continuous in meta-regression analysis.
Statistical Analysis
Prevalence estimates were calculated using the variance stabilizing double arcsine transformation (Freeman & Tukey, 1950) because the use of the inverse variance weight is inappropriate when dealing with binary data with low prevalence. Also, with this method, studies with zero prevalence can be included in the meta-analyses and their weight considered in combined estimation. In addition, due to zero to very low prevalence of CLU in some studies, confidence intervals (CIs) were calculated using the Wilson method (Newcombe, 1998; Wilson, 1927). The presence of heterogeneity was assessed by Cochran’s Q statistic (p < .1), combined with I2 statistic for estimates of inconsistency within the studies. A value of 0% indicates no observed heterogeneity while 100% indicates significant heterogeneity. The I2 values above 75% were determined as significant heterogeneity warranting analysis with a random effect model to adjust for the observed variability (Higgins, Thompson, Deeks, & Altman, 2003). Potential source of heterogeneity was further explored through subgroup analyses and meta-regression model using likelihood ratio method. In the meta-regression model, variables were examined both in univariable and in a multiple variable models at a significance level of <.05. All variables were entered in the model (Model 2), because all factors in univariable model (Model 1) were significant. To avoid model instability, only those variables that were significant in Model 2 were entered in speared multiple regression model (Model 3). Egger’s test was employed to examine potential publication bias. All analyses were conducted using Stata version 11.2 (Stata Corp LP, College Station, TX) with “metan,” “metareg,” “metabias,” and “propcii” commands.
Results
Descriptive Statistics
In an electronic database search, 4,000 citations were retrieved, producing a total of 33 studies from 12 provinces of Iran (Ahmadi, Alishahi, & Alavi, 2004; Ahmadi & Hasani, 2003; Ahmadi & Javadpour, 2001; Ahmadi, Javadpour, & Pridmore, 2009; Ahmadi, Maharlooy, & Alishahi, 2004; Ahmadi & Yazdanfar, 2002; Attari, Asgary, Shahrokhi, Naderi, & Shariatirad, 2012; Dehghani, Zare, Dehghani, Sedghi, & Pormovahed, 2009; Farhadinasab, Allahverdipour, Bashirian, & Mahjoub, 2008; Ghavidel et al., 2012; Jodati, Shakurie, Nazari, & Raufie, 2007; Meimandi et al., 2005; Mohamadi, Shobeiri, & Mahgoob, 2012; Mohammad Poorasl et al., 2007; Mohammadkhani, 2012; Mohammadpoorasl et al., 2012; Mohtasham-Amiri, Jafari-Shakib, & Khalili-Moosavi, 2011; Momtazi & Rawson, 2010; Mortazavi Moghadam, Madarshahian, Tabiei, Pejmankhah, & Sadeghi, 2009; Mustafa et al., 2010; Najafi et al., 2007; Najafi, Zarrabi, Shirazi, Fekri, & Mohseni, 2009; Nakhaee, Ziaaddini, & Karimzadeh, 2009; Nazarzadeh, Bidel, & Carson, 2014; Sajjadi & Sajjadi, 2012; Sarajzadeh & Feizi, 2003; Shams Alizadeh, Moghadam, Mohsenpour, & Rostami Gooran, 2008; Talaei, Mokhber, Fayyazi Bordbar, Javanbakht, & Samari, 2008; Taremian, Bolhari, Peirovi, & Tabatabaei, 2005; Zardkhaneh et al., 2011; Zarrabi et al., 2009; Ziaaddini, Zare-zadeh, & Heshmati, 2006; Ziaaddini, Ziaaddini, & Nakhaee, 2013). The characteristics of each included study are reported in Tables 1 to 3 and Figure 1 outlining the details related to the selection process of studies. The studies included a total of 30,972 adolescents and youth (12,931 high school students and 18,041 college students). Fourteen articles reported prevalence in men (n = 12,788), nine included mixed gender samples (n = 9,177), and nine investigated women (n = 9,007). The average age of men (reported in 13 reports) and women (from 7 reports) were 19.6 and 20.8 years, respectively. Average age of the mixed samples was 20.3 years (from 27 studies). Overall, quality assessment of all studies produces medium quality.
Description and Characteristics of Included Studies (Only Male).
Note. NM = not mentioned; DSM-IV = Diagnostic and Statistical Manual of Mental Disorders, fourth edition; C-D = census of dormitory; C-A = census of all student in the city; Q = questionnaire; dash (—) in prevalence column indicates not mentioned.
Survey includes five universities from whole of Iran: Tehran University, Razi University of Kermanshah, Ferdowsi University of Mashhad, Isfahan University of Technology, and Shahid Bahonar University of Kerman.
Description and Characteristics of Included Studies (Only Female).
Note. NM = not mentioned; DSM = Diagnostic and Statistical Manual of Mental Disorders, fourth edition; C-A = census of all student in the city; Q = questionnaire; dash (—) in prevalence column indicates not mentioned.
Survey include five universities from whole of Iran: Tehran University, Razi University of Kermanshah, Ferdosi University of Mashhad, Isfahan University of Technology, and Shahid Bahonar University of Kerman.
Description and Characteristics of Included Studies (Studies in Which Men and Women Were Not Separated).
Note. NM = not mentioned; Q = questionnaire; DSM = Diagnostic and Statistical Manual of Mental Disorders; ELISA = Enzyme–Linked Immunosorbent Assay; em dash (—) in prevalence column indicates not mentioned.

Flowchart of the reviewing process for prevalence of cannabis lifetime use in Iran.
Pooled Prevalence and Subgroup Analysis
A potential publication bias was detected for CLU (Egger’s test β0 = −1.4; p < .001). CLU was reported in 33 surveys with a random effects pooled prevalence of 4.0% (95% CI = 3.0% to 5.0%; Figure 2). Estimates ranged from 0.0% to 23.0% with substantial heterogeneity (χ2 = 1,424.5, p < .001, I2 = 97.8%). Higher prevalences were identified in smaller studies with smaller weight (Figure 3). Moreover, subgroup analyses based on gender produce evidence of significant difference between three subgroups (Figure 4). In mixed, male, and female subgroups, prevalence were 2.0% (95% CI = 1.0% to 3.0%), 6.0% (95% CI = 4.0% to 9.0%), and 3.0% (95% CI = 2.0% to 6.0%), respectively. Further exploration according to level of education and method of sampling were conducted. For high school and college students subgroups prevalences were 5.0% (95% CI = 3.0% to 7.0%), and 2.0% (95% CI = 2.0% to 3.0%), respectively (Figure 2) and in multistage, random, stratified sampling and census, pooled prevalences were 4.0% (95% CI = 3.0% to 5.0%), 1.0% (95% CI = 0.0% to 3.0%), 2.0% (95% CI = 0.0% to 4.0%), and 11.0% (95% CI = 3.0% to 19.0%), respectively (Figure 5).

Forest plot of studies related to prevalence of cannabis lifetime use in young Iranian population. Data are presented separately for college students and high school students.

Forest plot of studies related to prevalence of cannabis lifetime use in young Iranian population. Studies are sorted by increasing weight in a random effects model.

Forest plot of studies related to prevalence of cannabis lifetime use in young Iranian population. Data are presented separately for male, female, and mixed.

Forest plot of studies related to prevalence of cannabis lifetime use in young Iranian population. Data are presented separately for sampling methods.
Current, last week, and last month prevalences reported only in few surveys, consequently pooled prevalences were not calculated for them. (Tables 1-3).
Meta-Regression Model
In univariable meta-regression analyses, newer studies had higher prevalences of CLU (β for study date [year] = 0.0008, p < .001), and where the average age was ≥18 years (vs. ≤18), the prevalence was higher (β = 0.007, p < .001). The college versus high school students had higher prevalence (β = 0.008, p < .001). Also in category of ≥1,000 sample size versus <1,000, prevalence was higher (β = 0.009, p < .001), but where gender was male versus female and mixed, and in multistage, stratified, and random sampling versus census, prevalences were lower (β = −0.003, p < .001; β = −0.111, p < .001), respectively (Table 4, Model 1). In a multiple variable meta-regression model, all of the variables were significant except average age (Table 4, Model 2). Consequently, average age was excluded from the Model 3. Model 3 has identified that the prevalence is significantly higher in college students (β = 0.089, p < .001), and in male gender (β = 0.017, p < .001), and is lower in studies with sampling versus census studies (β = −0.096, p < .001).
Meta-Regression Analysis a of Cannabis Prevalence Heterogeneity Suspected Variables.
Note. Significant associations (p ≤ 0.05) are indicated in bold typeface; Ref = reference group.SE = standard error
Likelihood ratio method.
Univariable model.
Multiple variable model.
Variables entered into the model if p ≤ 0.05.
Discussion
This systematic review of CLU prevalence in Iranian youth identified 33 surveys including 30,972 subjects. There are four main findings: (a) prevalence of CLU in young Iranian people was 4.0%, (b) higher prevalence in male compared with female and in college students versus high school students, (c) higher prevalence in census based studies versus sampling based studies, (d) little changes to prevalence during past 30 years.
Much research has been conducted to confirm high use of cannabis among adolescents and young adults across communities. The 2007 European School Survey Project on Alcohol and other Drugs had reported that CLU among students aged 15 to 16 years in Europe is 19.0% in average (among 35 countries) and ranged from 3.0% in Armenia to 45.0% in the Czech Republic (Greydanus et al., 2013). In addition, in Canadian college students, the reported prevalence of cannabis abuse was approximately 30.0% (Fischer et al., 2013). According to the United States Centers for Disease and Prevention’s Youth Risk Behavioral Surveillance reports, CLU prevalence among high school students was 31.3% in 1991, 47.2% in 1999, and 36.8% in 2009 (Eaton et al., 2010). Compared with these international reports, the prevalence of CLU was considerably less in Iran. A primary explanation for this small prevalence might be legal limitations of cannabis and other substances in Iran. In addition, because of the historical background of opium and derivatives use, the prevalence of opium use is higher than other substances such as cannabis in Iran (Siassi & Fozouni, 1980). On the other hand, most of the included studies in these meta-analyses assessed CLU using questionnaires or direct questions. Consequently, false negative answers can be considered as potential low value for pooled prevalence.
This study’s finding that the sampling method of studies was significantly associated with prevalences suggests that researchers need to interpret their cross-sectional findings with regard to their methodology specifically with regard to sampling method for which we reported some underestimation. One of the interesting results of this study were the differences in subgroup analyses and meta-regression results on prevalence between high school and college students. Subgroup analysis indicated that prevalence in high school students is higher than for college students but the meta-regression model had indicated contrast result. It seems possible that these contrasts are due to the confounding role of other variables that were not considered in subgroup analysis. After adjustment for all possible confounders in meta-regression, it was revealed that the prevalence is significantly higher in college students versus high school students. Consequently, the results for subgroup analysis in meta-analyses studies should be interpreted with consideration of other possible confounders. One finding of this review is that future research in prevalence of any drug should provide clear information about study questionnaires and their validity and reliability.
This review is important for enhancing our methodological knowledge of the epidemiology of CLU in adolescents and young adults and provides a more definitive statement on the direction of prevalence changes. Also, it provides better estimation for prevalences in national studies such as the burden of diseases.
This study does have several limitations, as all pooled analyses contained significant heterogeneity and subsequently should be interpreted with caution. Most studies failed to report questionnaire validity, reliability, and response rate. These characteristics may be a significant source of heterogeneity (Fazel, Khosla, Doll, & Geddes, 2008). Moreover, a lack of comprehensive coverage for gray literature containing university databases and research projects, further limit the search.
In conclusion, the results of these meta-analyses indicate that prevalence of CLU in Iranian high school and college students are lower than for industrialized countries. In addition, gender, level of education, and methods of sampling are highly associated with variation in estimation of CLU prevalence.
Footnotes
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was financially supported by the Ilam University of Medical Sciences, and a research grant (Grant No. 909230) from the Prevention of Psychosocial Injuries Research Center.
