Cigarette smoking trajectories from adolescence to young adulthood: first report from the Middle East

PDF version

Hasti Masihay-Akbar1, Soha Razmjouei1, Elaheh Ainy2, Leila Cheraghi1,3, Fereidoun Azizi4, Parisa Amiri

1. Research Center for Social Determinants of Health, Research Institute for Endocrine Sciences, Shahid Beheshti University of Medical Sciences, Tehran, Islamic Republic of Iran

2. Department of Vice Chancellor Research Affairs, Shahid Beheshti University of Medical Sciences, Tehran, Islamic Republic of Iran

3. Department of Biostatistics and Epidemiology, Research Institute for Endocrine Sciences, Shahid Beheshti University of Medical Sciences, Tehran, Islamic Republic of Iran

4. Endocrine Research Center, Research Institute for Endocrine Sciences, Shahid Beheshti University of Medical Sciences, Tehran, Islamic Republic of Iran

Correspondence:

Paris Amiri

Research Center for Social Determinants of Health

Research Institute for Endocrine Sciences

Shahid Beheshti University of Medical Sciences

PO box 19395-4763, Tehran, Islamic Republic of Iran

Tel.: +98 21 240 9309, Fax: +98 21 240 2463

E-mail: عنوان البريد الإلكتروني هذا محمي من روبوتات السبام. يجب عليك تفعيل الجافاسكربت لرؤيته.

Abstract

Background: This study is the first effort in the Middle East to identify cigarette-smoking trajectories and their predictors, from adolescence to young adulthood.

Methods: Using data from the Tehran Lipid and Glucose Study, 1169 adolescents (12–18 years old) were followed into their young adulthood (28–32 years old), from 2002 to 2016. Cigarette smoking (combination of quantity and frequency) was the outcome variable used for group-based trajectory modelling. After detecting the trajectories, the effects of independent variables (individual employment, education, physical activity, and paternal smoking, employment and education) on the trajectories were investigated.

Results: Three trajectories were detected: non-smokers (79%), experimenters (12%) and escalators (9%). Compared to girls, boys were approximately three times (OR=2.94, 95%CI: 2.32-3.24, P

Conclusion: Iranian adolescents follow three trajectories – non-smokers, experimenters and escalators – into young adulthood. Male sex, employment after high school, and living with a smoker father are associated with unfavourable smoking patterns. These findings provide valuable insight that can be translated into practical interventions tailored to the local context.

Keywords: Smoking trajectories, cigarette smoking, adolescence, young adulthood, longitudinal

Citation: Masihay-Akbar H, Razmjouei S, Ainy E, Cheraghi L, Azizi F and Amiri P. Cigarette smoking trajectories from adolescence to young adulthood: first report from the Middle East. East Mediterr Health J. 2023;29(11):xxx–xxx. https://doi.org/10.26719/emhj.23.115 Received: 25/01/23; Accepted: 24/07/23

Copyright: © Authors; Licensee: World Health Organization. EMHJ is an open access journal. All papers published in EMHJ are available under the Creative Commons Attribution Non-Commercial ShareAlike 3.0 IGO licence (CC BY-NC-SA 3.0 IGO; https://creativecommons.org/licenses/by-nc-sa/3.0/igo).


Declarations

Ethics approval and consent to participate

The study was approved by the Ethical Committee of Research Institute for Endocrine Sciences and the National Research Council of the Islamic Republic of Iran (EC 121). Informed consent was obtained from all individual participants included in the study. All procedures were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards.

Consent for publication

Not applicable.

Availability of data and materials

The data sets used and/or analysed during the study are available from the corresponding author on reasonable request.

Conflict of interest

The authors declare that they have no competing interests.

Funding

The authors received no financial support for the research, authorship and/or publication of this article.

Authors’ contributions

HM-A and PA designed the study. FA and PA contributed to the design of the TLGS. LC conducted the statistical analysis. HM-A, SR, EA and LC contributed to interpretation of data. HM-A and SR contributed to writing the original draft. PA, FA and EA reviewed and revised the draft manuscript. PA and FA supervised drafting. All authors read and approved the final manuscript.

Acknowledgements

The authors would like to thank the TLGS participants and personnel for their collaboration.

Background

The negative health consequences of cigarette smoking and its attributed morbidity and mortality are well known (1). As the most common type of tobacco use, cigarette smoking continues to be a public health problem. Globally, 22.3% of the population uses tobacco, 36.7% of all men and 7.8% of the world’s women (2). In the Islamic Republic of Iran, the prevalence of smoking and daily cigarette smoking among adults in 2016 was 14.1 and 9.7%, respectively. The prevalence of smokers in adolescents was 3.4% (3,4).

Cigarette smoking is a complex and dynamic behaviour that follows a progressive pattern throughout the life course (5). Most adult smokers start smoking in adolescence (6), making it a critical period in developing smoking behaviour. The importance of this period has made it the target of most smoking-prevention interventions for many years. Although effective in preventing and delaying smoking initiation, these interventions often fail to stop smoking progression (7).

In adulthood, when smoking behaviour has already been established, the focus shifts towards cessation interventions. Therefore, the path between experimentation with tobacco products and established smoking behaviour that occurs during young adulthood is often neglected (8,9).

The dynamicity of smoking behaviour makes it logical to examine its longitudinal pattern to understand how and when this behaviour forms in different subpopulations. Trajectory analysis can help describe the longitudinal trend of cigarette use, identify subgroups at risk of sustained and heavy smoking, and pinpoint the optimum time to intervene (10,11).

Many North American and European cohort studies have used this method to identify smoking trajectories (11). While most of them traced smoking into adolescence, few have examined smoking behaviours from adolescence through young adulthood (12–17). Previous studies also varied in duration of follow-up (between 1.5 and 23 years) and measurement intervals (between three months and 4.5 years) (11). Researchers used different indicators for smoking in previous articles, including smoking frequency (number of days of smoking in a period of time), smoking intensity (number of cigarettes smoked in a defined time) and a combination of intensity and frequency.

Given the methodological and analytical disparities, the observed trajectories in prior research diverged in terms of quantity, ranging from two to six. Studies with fewer data points detected a lower number of trajectories. In terms of typology (the visual inspection of the curves), trajectories with lowest level of cigarette consumption were the most prevalent, followed by increasing pattern.

In the existing body of literature, various factors have been identified as predictors of cigarette smoking-trajectory membership. These include age, sex/gender, race, parental education and smoking, peers’ smoking, other substance use, and depression. While some of these predictors remain consistent over time, such as race and sex, others, such as socioeconomic and behavioural predictors, are subject to change throughout the lifespan. However, these results are not generalizable to the countries with socio-cultural differences.

Results from longitudinal studies in Middle Eastern countries, including Iran, indicate high progression rates from the primary stages of the smoking behaviour spectrum to established consumption in school years (18–20). Yet, with the maximum follow-up time of three years, few of them have delved into changes in smoking behaviour in the transition from adolescence to young adulthood.

The Tehran Lipid and Glucose Study (TLGS), one of the unique longitudinal cohorts in the Middle East, provides an opportunity to conduct advanced analysis on tobacco-related data and prepares preliminary data for designing age- and stage-sensitive interventions in the region. Leveraging the TLGS data, this study is designed to fill the gaps mentioned above and enrich the smoking-related literature on the topic.

The first objective of the study is to identify the developmental trajectories of cigarette smoking from adolescence to young adulthood. The second objective is to identify the individual and parental predictors of the membership of each trajectory.

Methods

Sampling and participants

Data is drawn from TLGS, an ongoing family-based cohort study to capture non-communicable disease prevalence and risk factors in an urban population of Iran. The multistage cluster random sampling method was used to select the target sample. In stage one, three of the 20 health care centres in district 13 of Tehran were selected, based on geographical location and data availability.

In the second stage, the data of 15 005 residents (≥3 years) were randomly collected from those health care centres. TLGS included a cross-sectional phase (1999–2001) with five follow-up examinations at three-year intervals (2002–2004; 2005–2007; 2008–2010; 2011–2013; and 2014 –2016). Details of TLGS design and sampling were previously published (21).

This study was restricted to participants aged 12–18 years (n=1567), who were followed from 2002 to 2016. After excluding 396 individuals with at least three missing in the smoking variable and two individuals who died over the follow-ups, analysis was conducted on the data of the remaining 1169 participants. The median follow-up for the whole population was 12.5 years.

Measurements

Smoking: As per previous use in the literature (11,22), in order to concisely capture diverse patterns of smoking behaviour (especially occasional smoking, whose pattern may not be succinctly captured solely with frequency and intensity variables), we created a metric for smoking by combining intensity and frequency. At each follow-up measurement, smoking data was acquired using standard questionnaires. Adolescents (≤18 years) reported their quantity (number of cigarettes in one day) and frequency of smoking (number of days in the past 30 days).

We categorized each variable as follows: quantity:

Individual covariates: Employment status was defined as follows: before age 18, when attending school is compulsory, all participants were considered unemployed; and after 18 years of age (adulthood), anyone who reported a change in employment status to having a job was categorized as employed; otherwise, they were considered unemployed. A highly educated person was defined as someone who has studied for more than 12 years or reports university education after finishing high school. Binary sex was also added as covariate.

The Persian translation of the Modifiable Activity Questionnaire was used to assess physical activity (PA) in adolescents (23). PA type (leisure and occupational activities), frequency and duration over the past 12 months were reported. The number of minutes/weeks for each activity was multiplied to its metabolic equivalent. Subsequently, an individual’s total PA was calculated by adding up domains of PA, and levels were defined as low (Parental covariates: For each parent, employment was defined based on having or not having a job. Regarding education, parents were categorized as illiterate/primary, secondary/diploma, and higher. For their smoking status, parents were classified as current smokers (smoking always or occasionally) and non-smokers.

Statistical analysis

We used group-based trajectory models (GBTM) to identify cigarette-smoking trajectories from age 12 to 32. The analytic sample was 12–18-year-olds of baseline who became 28–32-year-olds in the last follow-up. We estimated zero-inflated Poisson (ZIP) trajectory models with a user-written programme called TRAJ (24). The zero-inflated model is used when analysing count data that has excess zero counts (zero-inflated). The excess zeros are suggested to be generated by a separate process from the count values and can be modelled independently. Therefore, the ZIP model combines two distributional forms: a Poisson count model and a Logit model for predicting excess zeros (25).

To assess the appropriateness and suitability of the model, we created and examined visual representations depicting the distribution of the smoking index at baseline and subsequent follow-up examinations. Through the visual representations, we identified a significant number of zero values in the data, which aligned with our expectations given the age range of the participants, especially at baseline. These zero values correspond to individuals who can be categorized as non-smokers, indicating the absence of reported smoking.

At the beginning of the study, 98% of individuals had a smoking index value of zero. As the follow-up examinations progressed, the percentage of individuals with a zero smoking index decreased gradually, measuring 95%, 90%, 86% and 81% in the first through the last follow-up examinations, respectively. Subsequently, we investigated the possible number of latent trajectories by a series of models considering several linear and polynomial (cubic, quadratic) specifications of the smoking index as a function of age.

For better model detection, age was centralized at 20 years. GBTM (10) was used to identify subgroups of participants who shared similar underlying trajectories of the outcome variables. In a GBTM framework, the latent growth factors, that is, intercepts and slopes, determine each group’s trajectory. The intercept refers to the initial smoking index for 20-year-old participants in phase two, and the slopes correspond to the rate of linear or non-linear change in the smoking trajectory across assessments. Adding quadratic or cubic growth parameters to the model helped to capture nonlinear trajectories.

We began with a single model consisting of one group and then increased the number of groups until the number of trajectories that best fit the data was identified based on the Bayesian Information Criterion (BIC), the average posterior probability of group membership (APP) and the odds of correct classification (OCC).

We tested models (see Table 1) with one to four groups and different specifications of smoking index as a function of age for each group and for zero-inflation probability. The three-group model with linear function of age for the first group and cubic function of age for the second and third groups was selected (model 12). The optimal model was selected as the lowest BIC (BIC= -2509.77), highest APP greater than 0.70, OCC more than five, and the minimum of 5% of the total sample for the size of each class. The same model with three groups and a quadratic function of age for zero-inflation probability (model 26) had a slight decrease in BIC value; however, as this model was more complex and model 12 was better interpretable, the former was chosen.

[ADD Table 1: Model fit indices for the trajectory model analysis NEAR HERE]

In this analysis, the BIC value increased for models with more than three trajectory groups. The models with four or more groups did not fit well and, as a result, the three-group model was considered optimal. For each subject, the model provides the probability of belonging to each of the identified trajectories and then assigns the subject to the group where it has the highest chance. Based on the existing literature, the selected model’s substantive and theoretical interpretability were also considered.

While classifying the GBTM model, there was a possibility to add predictors to increase the accuracy of fit indices and shapes of trajectories. We included the effect of some previously studied covariates associated with smoking behaviours as time-stable predictors (individual sex, employment, education, physical activity, and paternal smoking, employment and education) on the defined trajectory groups. These covariates were included as time-invariant.

Missing data was handled using the STATA programme (24) in this study. The TRAJ programme uses the maximum likelihood method to estimate parameters, including group sizes and shapes of trajectories. When the data is assumed missing at random, this method generates an asymptotically unbiased parameter. Subjects with missing data are included in the analysis but only available data for each subject are used (10).

To determine whether “missingness” affected the results (whether missing data is random or not), we explored missing data patterns and their effects on investigated smoking trajectories. A variable was created to represent the number of missing data among five data points for each participant and compared among three smoking-trajectory groups. In this study, there were potentially five data points for each participant.

In running the trajectory analysis, only participants with three or more available measurements of smoking are included; therefore, included participants may have zero, one or two missing smoking measurements. Having no or one missing data point in smoking was compared with having two missing data points, in prediction of trajectory membership. In addition, four dummy variables were made, representing missing data in the middle of follow-up examinations, and investigating the associations between any of these dummy variables and smoking-trajectory groups.

We conducted the trajectory analysis, including identification of the optimal trajectories and determining associated variables with the explored trajectories in STATA software version 16 and the remained analysis in IBM SPSS Statistics version 26; two-sided P values of

Results

The tested models in the GBTM are shown in Table 1. The APP of each group in the selected model ranged from 0.93 to 0.95, the estimated group sizes were close to the actual ones, and the OCC exceeded five, suggesting a good fit and accurate group assignment. The three groups were named as non-smokers (n=918, 79%), experimenters (n=146, 12%) and escalators (n=105, 9%).

As presented quantitatively in Table 2 and visually in Figure 1, the non-smokers represented the largest group in the sample, which followed intercept and a very slight linear trajectory.

[ADD Table 2: Parameter estimates of the trajectory model with a zero-inflated Poisson distributional form NEAR HERE]

[ADD Figure 1 NEAR HERE]

Both the experimenter and escalator groups had positive slopes, with negative quadratic function and positive cubic function of age. However, the magnitude of estimated parameters for slope, quadratic and cubic functions was smaller in the escalator group, causing a more-steeply increasing graph overall than the experimenter group.

The experimenters showed a smoking index similar to the non-smokers until 16 years of age and gradually increased their cigarette use into a young adulthood steady phase, which started to rise again at the age of 30. The escalators diverged from non-smokers and experimenters at age 14 and increased their cigarette smoking frequency/quantity with a negative quadratic trajectory followed by a late positive cubic trajectory.

Socio-demographic and behavioural characteristics of adolescents and their parents are represented in Table 3. The mean baseline age of adolescents was 15.14±1.97, 15.23±1.81 and 15.21±2.05 years, in non-smokers, experimenters and escalators, respectively. While boys comprised 35% of non-smokers, they formed the majority of the experimenters (70%) and escalators (89%). Experimenters were more physically active. No significant difference was observed in parental characteristics between trajectories, except for paternal smoking (P

[ADD Table 3: Baseline characteristics of adolescents and their parents according to smoking-trajectory group NEAR HERE]

Table 4 presents the effects of individual covariates in predicting trajectory membership, including participants’ sex, education, employment, PA and missingness as well as paternal smoking status and missingness. Boys were three times and 25 times more likely than girls, respectively, to be experimenters and escalator smokers rather than non-smokers.

[ADD Table 4: Associated factors of the identified smoking trajectories using zero-inflated Poisson model NEAR HERE]

Receiving university education decreased the odds of placing in the escalator trajectory by 18% (P=0.001). Getting employed showed undesirable effects, as it doubled the odds of experimenter and escalator group membership (respectively, OR=2.00 and OR=2.33, P

Our model shows that, compared to those with no or one missing data point, adolescents with two missing data points are more likely to be in the escalator group (OR=2.03, 95% CI: 1.41-2.65, P=0.02). Paternal smoking (at any point during the study) is associated with 1.88- and 2.23-times higher odds of following an experimenter and escalator trajectory, respectively.

Discussion

Using data from TLGS, a population-based cohort study in the Islamic Republic of Iran, we followed 1169 Iranian adolescents for 15 years and identified three developmental trajectories of cigarette smoking: non-smokers (79%), experimenters (12%) and escalators (9%). Male sex, employment after age 18 and having a smoker father increased the likelihood of experimentation with smoking and continuing at a higher intensity. Participants with a university education were less likely to be escalator smokers. Physical activity level had no predictive value in trajectory-group membership.

In this study, using GBTM, we detected three distinct trajectories of cigarette smoking from adolescence to young adulthood. The majority of our sample was non-smokers. The second prevalent group was experimenters, with approximately one-sixth of the prevalence of non-smokers.

Experimenters started smoking at about 16 years of age at a low rate, increased their consumption from 18, and maintained a low-steady state between 24 and 30 years of age. This group seemed to increase cigarette smoking from the early fourth decade of their life.

Escalators, with less than 10% prevalence, started smoking increasingly from early adolescence (before 14 years of age). This increasing trend became steeper in young adulthood, highlighting the importance of this period.

Longitudinal studies from the Middle East region indicate high progression rates from the primary stages of the smoking behaviour spectrum to established school-age consumption. A large-scale study in Jordan showed a 38% increase in the frequency and intensity of smoking during three school years (18).

In Iran, the transition rate from never smoking to experimentation has been reported to be between 10% and 14% in different adolescent samples. Moreover, 16-17% of novice smokers reported continued regular smoking after one to three years (19,20). However, investigating the developmental pathways of smoking behaviour using different trajectory analysis methods has only been conducted in North American and European populations. Variations in the density of measurements, smoking indicators and time axis resulted in various numbers of trajectories between two and six (11).

In line with our results, the most common finding in these studies was a large non-smoker group, which made up the majority of the sample. Occasional or light smoker, early-onset, and late-onset stable or increasing groups were other identified trajectories. Previous studies with relatively larger sample sizes, more data points or shorter time intervals detected additional smoking patterns, such as quitters or decliners (15,17,26–28). These minority groups include adolescents who, despite smoking initiation, did not become established smokers.

In contrast, we did not detect such decreasing trends in our sample. The absence of a quitter group may be influenced by the characteristics of our study sample. It is possible that the prevalence of individuals who successfully quit smoking during the observed time period was relatively low or that the quitting patterns did not meet the criteria for forming a separate trajectory group. In this case, targeted interventions are needed in Iran to prevent smoking acceleration as mid- and late-adolescents move towards young adulthood.

On the other hand, the identification of specific trajectory groups depends on the underlying distribution of smoking behaviours and the fit of the statistical model to the observed data. Therefore, it is also plausible that the patterns of smoking cessation in our study did not exhibit clear trajectory profiles that could be distinguished from other trajectory groups. In this regard, future studies are needed to examine quitting trajectories and provide a comprehensive understanding of smoking cessation patterns.

Generally, participants’ sex, education, employment and paternal smoking were predictive of trajectory-group membership in this study. In line with previous studies, we found that boys are more likely than girls to follow the experimenter and escalator trajectory. However, the significant difference between girls and boys is consistent with previous studies in Iran, which might be related to the taboo of female smoking, and underreporting.

In this study, continuing education beyond high school is associated with reduced odds of smoking escalation. Education has been previously examined in relation to smoking trajectories with different indicators and definitions, such as academic performance through school years and years/level of education (14,16). In line with our results, many studies found that education is a predicting factor in trajectory-group membership. In contrast, others have investigated education as an outcome of long-term smoking, showing that those who start smoking sooner and with increasing intensity have lower education and academic performance (17).

Our data also showed that being employed after 18 is associated with membership of both riskier trajectory groups (experimenter and escalator). This finding contradicts previous studies that found unemployment increases the chance of being in high-risk groups. It might be interpreted that, in Iranian youth, employment does not act as a direct socioeconomic indicator but as a factor that could interfere with continued education after high school or cause low academic performance. As stated in previous studies, young Iranians in wealthier neighbourhoods smoke more than those in poorer areas (29,30).

Our results show that paternal smoking was associated with 1.88 times higher odds of being in the experimenter group and doubled the odds of following an escalating trajectory. This is consistent with many studies confirming the association between parental smoking and adolescent smoking initiation and intensity (11,31–33).

The intergenerational transmission could be due to more and easier access to cigarettes and weaker smoking policies at home. Although in previous studies, mothers’ smoking had a stronger effect on offspring tobacco use (32,33), studies of the TLGS sample did not find maternal smoking to be important in their children’s smoking behaviours. This could be due to Iranian women underreporting their smoking and acting more conservatively in the family in this regard.

To the best of our knowledge, this study is the first effort to identify smoking developmental pathways in the Middle East region. With a large sample size and long follow-up time, the TLGS cohort study allowed us to examine trajectories of smoking behaviours from adolescence to young adulthood in association with some influential factors.

In this study, we created a metric for smoking by combining intensity and frequency to better capture diverse patterns of smoking behaviour (especially among occasional and non-daily smokers). This method has been used in previous studies investigating trajectories of cigarette smoking (11,22). While our study’s predictors may align with previous findings, its significance lies in its focus on the Middle East region, and its potential to inform targeted interventions and policies that are culturally appropriate and address the challenges and dynamics of the Middle Eastern context.

Our findings have the potential to guide targeted and tailored interventions that address the needs and challenges faced by individuals within each identified trajectory. Given that the majority of the sample consisted of non-smokers, it is crucial to maintain and reinforce this non-smoking behaviour. Thus, prevention programmes should focus on promoting and strengthening the non-smoking trajectory.

Recognizing their trajectory, early intervention can be designed to target experimenters during late adolescence and young adulthood. These interventions should emphasize the risks associated with continued smoking and provide support for smoking cessation and behaviour change. Future studies are needed to investigate the triggers or reasons behind the increase in smoking during the early fourth decade of life to prevent escalation.

Targeted intervention for escalators who are at high risk for long-term smoking and associated health consequences should be developed and implemented to address this group’s needs, focusing on preventing the acceleration of smoking behaviour in young adulthood.

Limitations

In the analysis of this study, we were unable to capture the effects of variables that could vary over time due to data availability. Although these covariates were included as time-invariant, we defined their categories to detect one-point changes throughout the study. For example, a change in employment status was defined as a change from unemployment to having a job in any follow-up after age 18. More prospective research is needed to examine the effect of time-varying data on smoking trajectories.

Due to low prevalence of smoking in girls, we could not examine parental smoking influences on boys and girls separately in this study. Future studies with adequate sample size are needed to assess the sex-specific smoking trajectories in Iran. This study is also limited by self-reported questionnaires, which may have been prone to underreporting and recall bias, especially in women.

As the data were unavailable, other tobacco-related factors, such as mental health, family structure, peer smoking, attitudes and beliefs towards smoking, other risky behaviours and use of other substances like alcohol were not considered in this analysis.

Another limitation of our study is the potential presence of sampling bias. Our sample was drawn from urban areas of the capital city of Tehran, which may not be fully representative of the broader population in Iran or the Middle East. It is important to consider this potential bias when interpreting the results and applying them to wider populations. To mitigate this bias, future studies should aim for more diverse and representative samples that encompass different regions, socioeconomic backgrounds and cultural contexts within the Middle East. This would enhance the external validity and generalizability of the findings and provide a more comprehensive understanding of cigarette-smoking trajectories in the region.

Additionally, our study focused on adolescents and young adults, and the findings may not fully capture the smoking trajectories and predictors beyond this age range. Future research should consider extending the study duration and including participants across a broader age range to examine the long-term trajectories and predictors of smoking behaviour.

Conclusion

Iranian adolescents follow three trajectories – non-smokers, experimenters and escalators – to young adulthood. Male sex, employment after high school, and living with a smoker father are associated with unfavourable smoking behaviour patterns. On the other hand, continued education has protective effects against increasingly high levels of cigarette use.

References

1. National Center for Chronic Disease Prevention and Health Promotion (US) Office on Smoking and Health. The Health Consequences of Smoking—50 Years of Progress: A Report of the Surgeon General. Atlanta (GA): Centers for Disease Control and Prevention (US). 2014. PMID: 24455788.

2. World Health Organization‎. WHO report on the global tobacco epidemic, 2021: addressing new and emerging products. Geneva: World Health Organization. 2021. (https://apps.who.int/iris/handle/10665/343287).

3. WHO EMRO; 2018. Global Youth Tobacco Survey (GYTS), Fact Sheet Islamic Republic of Iran 2016. (https://portal-uat.who.int/fctcapps/sites/default/files/2023-04/iran_2018_annex-1_GYTS_factsheet_2016.pdf).

4. Sohrabi M-R, Abbasi-Kangevari M, Kolahi A-A. Current Tobacco Smoking Prevalence Among Iranian Population: A Closer Look at the STEPS Surveys. Front Public Health. 2020;8:571062. DOI: 10.3389/fpubh.2020.571062.

5. Mayhew KP, Flay BR, Mott JA. Stages in the development of adolescent smoking. Drug Alcohol Depend. 2000;59 Suppl 1:S61-81. DOI: 10.1016/s0376-8716(99)00165-9.

6. National Center for Chronic Disease Prevention and Health Promotion (US) Office on Smoking and Health. Preventing Tobacco Use Among Youth and Young Adults: A Report of the Surgeon General. Atlanta (GA): Centers for Disease Control and Prevention (US); 2012. PMID: 22876391.

7. Masihay-Akbar H, Amiri P, Cheraghi L, Momenan AA, Azizi F. The Cigarette Smoking Initiation and Continuation in Adolescents Undergoing a Long-Term Behavioral Intervention. Nicotine Tob Res. 2021;23(4):702–710. DOI: 10.1093/ntr/ntaa208.

8. Rath JM, Villanti AC, Abrams DB, Vallone DM. Patterns of tobacco use and dual use in US young adults: the missing link between youth prevention and adult cessation. J Environ Public Health. 2012;2012:679134. DOI: 10.1155/2012/679134.

9. Villanti AC, Niaura RS, Abrams DB, Mermelstein R. Preventing Smoking Progression in Young Adults: the Concept of Prevescalation. Prev Sci. 2019;20(3):377–384. DOI: 10.1007/s11121-018-0880-y.

10. Nagin D. Group-Based Modeling of Development. Boston (MA): Harvard University Press. 2005.

11. Ahun MN, Lauzon B, Sylvestre M-P, Bergeron-Caron C, Eltonsy S, O’Loughlin J. A systematic review of cigarette smoking trajectories in adolescents. Int J Drug Policy. 2020;83:102838. DOI: 10.1016/j.drugpo.2020.102838.

12. White HR, Pandina RJ, Chen P-H. Developmental trajectories of cigarette use from early adolescence into young adulthood. Drug and alcohol dependence. 2002;65(2):167–178.

13. Orlando M, Tucker JS, Ellickson PL, Klein DJ. Developmental trajectories of cigarette smoking and their correlates from early adolescence to young adulthood. J Consult Clin Psychol. 2004;72(3):400–410. DOI: 10.1037/0022-006X.72.3.400.

14. Brook DW, Brook JS, Zhang C, Whiteman M, Cohen P, Finch SJ. Developmental trajectories of cigarette smoking from adolescence to the early thirties: Personality and behavioral risk factors. Nicotine Tob Res. 2008;10(8):1283–1291. DOI: 10.1080/14622200802238993.

15. Dutra LM, Glantz SA, Lisha NE, Song AV. Beyond experimentation: Five trajectories of cigarette smoking in a longitudinal sample of youth. PLoS One. 2017;12(2):e0171808. DOI: 10.1371/journal.pone.0171808.

16. Clergue-Duval V, Mary-Krause M, Bolze C, Fombonne E, Melchior M. Early Predictors of Trajectories of Tobacco Use Level from Adolescence to Young Adulthood: A 16-Year Follow-Up of the TEMPO Cohort Study (1999-2015). Eur Addict Res. 2019;25(1):2-9. DOI: 10.1159/000495759.

17. Lenk KM, Erickson DJ, Forster JL. Trajectories of Cigarette Smoking From Teens to Young Adulthood: 2000 to 2013. Am J Health Promot. 2018;32(5):1214–1220. DOI: 10.1177/0890117117696358.

18. Jaber R, Mzayek F, Madhivanan P, Khader Y, Maziak W. Predictors of Cigarette Smoking Progression Among a School-Based Sample of Adolescents in Irbid, Jordan: A Longitudinal Study (2008–2011). Nicotine Tob Res. 2016;18(4):403–409. DOI: 10.1093/ntr/ntv100.

19. Mohammadpoorasl A, Fakhari A, Shamsipour M, Rostami F, Rashidian H. Transitions between the stages of smoking in Iranian adolescents. Prev Med. 2011;52(2):136–138. DOI: 10.1016/j.ypmed.2010.11.024.

20. Mohammadpoorasl A, Nedjat S, Fakhari A, Yazdani K, Fotouhi A. Predictors of transition in smoking stages in Iranian adolescents: latent transition analysis. East Mediterr Health J. 2014;20(5):330–339. PMID: 24952291.

21. Azizi F, Ghanbarian A, Momenan AA, Hadaegh F, Mirmiran P, Hedayati M, et al. Prevention of non-communicable disease in a population in nutrition transition: Tehran Lipid and Glucose Study phase II. Trials. 2009;10(1):5. DOI: 10.1186/1745-6215-10-5.

22. Costello DM, Dierker LC, Jones BL, Rose JS. Trajectories of smoking from adolescence to early adulthood and their psychosocial risk factors. Health Psychol. 2008;27(6):811–818. DOI: 10.1037/0278-6133.27.6.811.

23. Delshad M, Ghanbarian A, Ghaleh NR, Amirshekari G, Askari S, Azizi F. Reliability and validity of the modifiable activity questionnaire for an Iranian urban adolescent population. Int J Prev Med. 2015;6:3. DOI: 10.4103/2008-7802.151433.

24. Jones BL, Nagin DS. A Note on a Stata Plugin for Estimating Group-based Trajectory Models. Sociological Methods & Research. 2013;42(4):608–613. DOI: 10.1177/0049124113503141.

25. Grimm KJ, Stegmann G. Modeling change trajectories with count and zero-inflated outcomes: Challenges and recommendations. Addict Behav. 2019;94:4–15. DOI: 10.1016/j.addbeh.2018.09.016.

26. Bradshaw M, Kent BV, Davidson JC, De Leon S. Parents, Peers, and Trajectories of Cigarette Smoking: A Group-Based Approach. Youth Soc. 2021;53(4):676–694. DOI: 10.1177/0044118x19862450.

27. Karp I, O’Loughlin J, Paradis G, Hanley J, Difranza J. Smoking trajectories of adolescent novice smokers in a longitudinal study of tobacco use. Ann Epidemiol. 2005;15(6):445–452. DOI: 10.1016/j.annepidem.2004.10.002.

28. Oura P, Rissanen I, Junno J-A, Harju T, Paananen M. Lifelong smoking trajectories of Northern Finns are characterized by sociodemographic and lifestyle differences in a 46-year follow-up. Sci Rep. 2020;10(1):16365. DOI: 10.1038/s41598-020-73334-3.

29. Rezaei F, Nedjat S, Golestan B, Majdzadeh R. Comparison of onset age and pattern of male adolescent smoking in two different socioeconomic districts of Tehran, Iran. Int J Prev Med. 2011;2(4):224–228.

30. Kelishadi R, Shahsanai A, Qorbani M, Esmaeil Motlagh M, Jari M, Ardalan G et al. Tobacco Use and Influencing Factors Among Iranian Children and Adolescents at National and Subnational Levels, According to Socioeconomic Status: The Caspian-IV Study. Iran Red Crescent Med J. 2016;18(5):e21858. DOI: 10.5812/ircmj.21858.

31. Mays D, Gilman SE, Rende R, Luta G, Tercyak KP, Niaura RS. Parental smoking exposure and adolescent smoking trajectories. Pediatrics. 2014;133(6):983–991. DOI: 10.1542/peds.2013-3003.

32. Gilman SE, Rende R, Boergers J, Abrams DB, Buka SL, Clark MA et al. Parental smoking and adolescent smoking initiation: an intergenerational perspective on tobacco control. Pediatrics. 2009;123(2):e274–281. DOI: 10.1542/peds.2008-2251.

33. Leonardi-Bee J, Jere ML, Britton J. Exposure to parental and sibling smoking and the risk of smoking uptake in childhood and adolescence: a systematic review and meta-analysis. Thorax. 2011;66(10):847–855. DOI: 10.1136/thx.2010.153379.