Article Text
Abstract
Objectives This study aimed to investigate the association between age-specific and sex-specific continuous metabolic syndrome severity score (cMetS-S) and the risk of developing type 2 diabetes mellitus (T2DM). Additionally, the study aimed to assess the added value of cMetS-S in predicting T2DM compared with traditional MetS criteria.
Design The study used a longitudinal cohort design, following participants for 18 years.
Setting The research was conducted within the Tehran Lipid and Glucose Study, a community-based study in Tehran, Iran.
Participants A total of 6957 participants aged 20–60 years were included in the study.
Interventions/exposures The cMetS-S of each participant was determined using age-specific and sex-specific equations and Cox proportional hazard regression models were used to analyse the association between cMetS-S and T2DM using continuous and quantile approaches.
Primary and secondary outcome measures The outcome measure was the association between cMetS-S and the development of T2DM during the 18-year follow-up.
Results A total of 1124 T2DM cases were recorded over 18 years of follow-up. In the fully adjusted model, a 1-SD increase in the cMetS-S was associated with future T2DM (HR 1.72; 95% CI 1.54 to 1.91). Men and women had HRs of 1.65 (95% CI 1.40 to 1.95) and 1.83 (95% CI 1.59 to 2.10) for T2DM per 1-SD increase in cMetS-S, respectively. Higher cMetS-S was associated with increased risk of diabetes in both prediabetic (HR 1.42;95% CI 1.23 to 1.64) and normoglycaemic individuals (HR 2.11;95% CI 1.76 to 2.54); this association was more significant in normoglycaemic individuals. Unlike the traditional-based MetS definitions, the cMetS-S improved diabetes prediction (p<0.001).
Conclusions The cMetS-S is strongly associated with future diabetes in prediabetic and normoglycaemic individuals independent of MetS components during a long term. As the relationship between cMetS-S and T2DM is more pronounced in normoglycaemic individuals than in those with pre-diabetes, implementing the evaluation of cMetS-S can serve as an early identification tool for individuals at risk of T2DM prior to the onset of pre-diabetes.
- Diabetes & endocrinology
- EPIDEMIOLOGIC STUDIES
- EPIDEMIOLOGY
Data availability statement
Data are available on reasonable request. Datasets generated during and/or analysed during the current study are not publicly available due to institutional policies but are available from the corresponding author on reasonable request.
This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.
Statistics from Altmetric.com
STRENGTHS AND LIMITATIONS OF THIS STUDY
We evaluated our newly developed metabolic syndrome severity score (MetS-S) in a large population-based sample of Iranian adults who were followed up for almost two decades.
This is the first study to investigate the application of age-specific and sex-specific continuous MetS-S (cMetS-S) in type 2 diabetes mellitus (T2DM) prediction.
We displayed the association of standardised cMetS-S with T2DM on a continuous scale (for 1-SD increase), allowing for comparison and generalisation of the results from different studies.
We only included non-elderly participants, aged between 20 and 60 years.
This study was conducted solely among the urban residents of Tehran, therefore, its applicability to the rural population remains unknown and can be considered a limitation.
Introduction
Type 2 diabetes mellitus (T2DM) is a major global health issue affecting more than 400 million individuals worldwide. Despite being preventable, T2DM poses a significant burden on the health system and economy, with an annual global health expenditure of US$966 billion.1 2 Metabolic syndrome (MetS), a constellation of five cardiometabolic risk factors, is associated with an approximately fivefold increased risk for T2DM.3 According to the traditional joint interim statement (JIS) MetS criteria,4 an individual with abnormal values for three or more of the following components are categorised as having MetS: abdominal obesity, increased blood pressure, increased triglycerides (TG), low high-density lipoprotein cholesterol (HDL-C) and increased fasting plasma glucose (FPG). This risk factor cluster is seemingly driven by underlying physiological alterations involving insulin resistance, oxidative stress, systemic inflammation, adipocyte dysfunction and cellular function abnormalities.5
The traditional definition, which classifies the individuals as having MetS or not, has limited its clinical utility. This dichotomous definition does not exhibit the extent of abnormality of the MetS components, making it impossible to compare the metabolic abnormality within the populations with or without MetS. In addition, with minimal changes in one or two MetS components’ values, an individual may be identified as having MetS or not.6 The term pre-MetS (preMetS) was developed for these individuals who do not fit the criteria of MetS but have a significant risk for adverse health outcomes.7 8 Last but not least, tracking changes (aggravation or improvement) of MetS is not possible with the traditional criteria and without a definition for severity.
To address these gaps, a few researchers developed continuous MetS severity scores (cMetS-S) using confirmatory factor analysis (CFA) and have shown their association with diabetes.9 10 We developed an age-specific and sex-specific cMetS-S score using CFA to demonstrate MetS severity by considering the weighted contribution of MetS components to MetS in the West Asian population.11 In the current study, we aim to verify the association of our newly developed MetS-S with T2DM in individuals with pre-diabetes and normal glucose metabolism independent of the MetS components and to find whether cMets-S has added value in predicting T2DM over traditional MetS definition.
Methods
Study design and data source
This study was conducted within the framework of the Tehran Lipid and Glucose Study (TLGS), a prospective population-based study performed on a representative sample of residents of district 13 of Tehran with the aim of determining the prevalence of risk factors for non-communicable disease. A total of 15 005 individuals aged 3 years and over who were residence of district No.13 of Tehran and were under the coverage of three medical health centres, were selected using multistage cluster random sampling method (February 1999–August 2001). All members of each family, including those not having risk factors, were invited for baseline measurements and will be followed every 3 years for 18 years. This population was selected through a multistage stratified cluster random sampling technique from the population of district 13 in Tehran. At the time, Tehran was composed of 20 urban districts and made up a population of 6.7 million (Iran National Census 1996). District 13 was chosen mainly because city-wide data showed a high rate of stability in that district. Also, the age distribution in district 13 was representative of the overall population in Tehran. In the current study, out of 15 005 study population, 10 813 participants aged 20–60 years were included.
Participants with diabetes (n=626), cancer (n=39), use of corticosteroid drugs (n=112) and pregnant women (n=78) were excluded at baseline. Moreover, we excluded subjects with missing baseline covariates (n=1704) and missing follow-up (n=1297). Thus, a total of 6957 participants were included and followed by approximately 3-year intervals for incident T2DM for 18 years (online supplemental figure 1).
Supplemental material
Clinical variables and data extraction
Information on age, sex, education level, medication use, medical history and family history of T2DM was collected by interview-based questionnaires. Examinations and anthropometric and biochemical measurements were performed by healthcare professionals.12 Using standard protocols, weight, height and waist circumference (WC) measurements were performed to the nearest 0.1 kg, 0.1 cm and 0.1 cm in the clinic, respectively. Weight was measured, while subjects were minimally clothed without shoes, using digital scales (Seca 707, Seca, Hanover, Maryland, USA; range 0.1–150 kg) and recorded to the nearest 100 g. Height was measured in a standing position without shoes, using tape metre while shoulders were in a normal alignment. Body mass index (BMI) was calculated as weight (kg) divided by square of height (m2). WC was measured at umbilical level, using an unstretched tape metre, without any pressure to body surface and HC at the maximal level over light clothing. Systolic blood pressure and diastolic blood pressure (SBP/DBP) were measured twice using a standardised mercury sphygmomanometer on the right arm, with at least a 30 s interval after 15 min resting in the seated position; the mean of the two readings was reported as the individual’s blood pressure.
Venous blood samples for laboratory tests were drawn from individuals after 12–14 hour of overnight fasting between 07:00 and 09:00 hours FPG was measured using an enzymatic colourimetric method with glucose oxidase. TG and HDL-C were assayed using enzymatic colourimetric methods. A 2‐hour postchallenge plasma glucose (2‐h PCPG) test was performed with 75 g glucose for all participants who did not use any glucose-lowering medication. Total cholesterol (TC) was assayed using the enzymatic colourimetric method with cholesterol esterase and cholesterol oxidase. The analyses were carried out using Pars Azmon kits (Pars Azmon, Tehran, Iran) and a Vitalab Selectra E autoanalyzer (Vital Scientific, Spankeren, The Netherlands). The intra-assay and interassay coefficients of variation were both o2.2% for FPG, and 0.5 and 2% for TC, respectively. A detailed description of protocols and laboratory techniques was previously published
Definitions
Smoking status was classified into smokers (use of any tobacco product daily or occasionally) and non-smokers (never and ex-smokers). Education was categorised based on the highest level of education completed: <6 years, 6–12 years and ≥12 years. The family history of diabetes was considered positive when at least one first-degree relative had T2DM. Individuals with less than 600 min per week of activity in the metabolic equivalents of task scale (METS) were described as low physically active. Hypertension was defined as SBP≥140 mm Hg, DBP≥90 mm Hg or antihypertensive medication use. Pre-diabetes status was defined when 100≤FPG<126 mg/dL and 140≤2 h-PCPG<200 mg/dL. T2DM was defined by FPG≥126 mg/dL, 2 h-PCPG≥200 mg/dL or the use of glucose-lowering drugs.
The cMetS-S was calculated for each individual at baseline using age-based and sex-based formulas to indicate MetS severity (online supplemental table S1). These formulas were derived using the CFA approach for the five MetS components (WC, TG, HDL-C, SBP and FPG) by considering the weighted contribution of each component to the latent MetS factor on an age-specific and sex-specific basis.11 The resulting cMetS-S was then standardised at mean=0 with SD=1, reflecting z-scores for each individual, with higher scores showing higher MetS severity.
Patient and public involvement
Patients or the public were not involved in the research’s design, conduct, reporting or dissemination plans.
Statistical analysis
The characteristics of the individuals in the study are presented as mean values±SD for continuous and frequencies with corresponding percentages for categorical variables. An independent t-test was conducted for continuous, and a χ2 test was performed for categorical variables to compare the groups. The incident rate for T2DM was calculated by dividing the total number of cases by the overall follow-up duration in person-years. The association between cMetS-S and T2DM was analysed using continuous and tertile approaches with Cox proportional hazard regression models in total, prediabetic and normoglycaemic subjects. Model 1 was adjusted for age and sex. Model 2 was adjusted for age, sex, family history of diabetes, physical activity and obesity (BMI ≥30 kg/m2). Model 3 included adjustments in model 2 as well as adjustments for the antihypertension and lipid-lowering drugs and model 4 further adjusted for MetS components, including high WC, SBP, TG, FPG and low HDL-C. Restricted cubic spline regression analysis was employed to demonstrate the relationship between cMetS-S and the risk of T2DM and to check the linearity of this relationship. We also compared the model performance of JIS and international diabetes federation (IDF)-based MetS and cMetS-S for T2DM using the Akaike information criterion (AIC) and log-likelihood ratio. The statistical analyses were conducted by using STATA V.14 (Stata) and R V.3.0.3 (R Foundation for Statistical Computing, Vienna, Austria). A two-sided p<0.05 was considered statistically significant.
Development of cMetS-S: cMetS-S had been developed previously using CFA on the five identified MetS components: WC, FPG, SBP, TG and HDL-C to consider the weighted contribution of these components to the unobserved latent variable of MetS.11 Several CFAs were performed on the eligible TLGS participants aged 20–60 years, and the results were presented in total and on age-specific (20–39 years and 40–60 years) and sex-specific basis. All the five MetS components in the models were standardised at mean=0 and SD=1 over the entire sample. One-factor model CFA was performed, and it was assumed that the measurement errors of the five components were not correlated. The factor loadings were indicative of the magnitude of the association between each component and the unobserved latent variable of MetS. The factor loadings >0.3 were considered to show a moderate correlation. Models were developed in the total population and age and sex subgroups with and without the assumption of the equality of factor loadings across the age and sex subgroups, respectively. Factor scores were produced using proper linear combinations of the variables. We calculated factor scores and cMetS-S using linear regression analysis with unstandardised MetS components to allow for potentially higher standardised scores within sex and age groups. The standardised factor coefficients in the final models were applied to calculate the cMetS-S for each individual. For ease of use in clinical settings, the MetS components and their factor coefficients obtained from the CFA have been back transformed so that actual values of MetS components can be placed in the equations. The resulting cMetS-S score can be standardised in each population and interpreted as z-scores (mean=0; SD=1), with greater values representing higher MetS severity. The performance of the overall and age-specific and sex-specific models was compared with various fit indices.
Results
A total of 6957 participants (42.17% male) with a mean age of 37.03±10.67 were included in the current study. Table 1 displays the baseline characteristics of the study population by cMetS-S quartiles. The prevalence of MetS was 31.02% and 26.41% according to the JIS and IDF criteria, respectively, with individuals in the upper quartiles having a higher prevalence of MetS. The mean values of age, BMI, WC, FPG, blood pressure and lipid indices (except HDL) increased on upper cMetS-S quartiles (p<0.001). The proportion of smokers and individuals who used antihypertension and lowering-lipid drugs was also greater towards the fourth cMetS-S quartile. In contrast, the proportion of individuals with low physical activity and a family history of diabetes did not differ across the cMetS-S quartiles. Baseline characteristics of the individuals are presented according to the presence of pre-diabetes (pre-diabetes vs normoglycaemia) in online supplemental table S2.
Baseline characteristics of the study population according to continuous metabolic syndrome severity score (cMetS-S) quartiles
A total of 1124 T2DM cases (16.15%) were recorded over 18 years of follow-up. T2DM incidence rates were about 12 per 1000 person-years in both men and women, with the highest cMetS-S quartile having an approximate sevenfold incidence rate of T2DM compared with the lowest quartile (table 2). After controlling for potential confounding variables, the risk of future T2DM increased with an HR of 1.87 (95% CI 1.74 to 2.00) per each SD increase in cMetS-S. Even after adjusting for MetS components (high TG, FPG, WC, SBP and low HDL-C), the cMetS-S was still associated with incident T2DM (HR 1.72; 95% CI 1.54 to 1.91 per 1-SD increase). This association was present in both men (HR 1.65; 95% CI 1.40 to 1.95) and women (HR 1.83; 95% CI 1.59 to 2.10). The results were consistent when assessing this association with the quartile approach, with the risk of T2DM increasing incrementally toward the upper cMetS-S quartile in the fully adjusted model (HR 6.45; 95% CI 4.48 to 9.27).
Cox proportional HRs for incidence of type 2 diabetes in total and according to sex
As depicted in table 3, cMetS-S was associated with T2DM in both normoglycaemic and prediabetic participants. A 1-SD increase in cMetS-S increased the risk of T2DM in prediabetic (HR 1.42; 95% CI 1.23 to 1.64) and normoglycaemic (HR 2.11; 95% CI 1.76 to 2.54) participants in the fully adjusted model. This association was consistent when investigating cMetS-S based on quartile groups, with incrementally increased risk for T2DM in higher cMetS-S quartiles; the highest quartile of cMetS-S in normoglycaemic and prediabetic individuals had HRs of 6.12 (95% CI 3.75 to 9.99) and 3.84 (95% CI 2.17 to 6.82), respectively.
Cox proportional HRs for incidence of type 2 diabetes in prediabetic and normoglycaemic individuals
The restricted cubic spline regression graphs were plotted according to the estimated HRs of participants for incident T2DM (figure 1). There was a non-linear positive association between the cMetS-S and future T2DM in all subgroups. The Kaplan-Meier curves showed poorer diabetes-free probability in higher cMetS quartiles over 18 years (online supplemental figure 2).
Restricted cubic spline regression analysis of continuous metabolic syndrome severity score (cMetS-S) with incident type 2 diabetes mellitus. Heavy central lines represent the estimated fully adjusted HRs in the logarithmic scale, with shaded ribbons denoting 95% CIs.
We conducted model fit analyses to determine whether cMetS-S could better predict T2DM than the traditional MetS definitions (JIS and IDF) and MetS components (table 4). Therefore, we created the initial model by including MetS components (high TG, WC, FPG, SBP and low HDL-C) and established T2DM risk factors. We then separately added cMetS-S, JIS and IDF-based MetS to the initial model in models 2, 3 and 4. Unlike the JIS- and IDF-based MetS definitions, adding cMetS-S to the initial model showed significant outperformance by reducing the AIC value by 94.5 on df=1 for diabetes (p<0.001).
The added value of continuous metabolic syndrome (MetS) severity score, MetS defined by JIS and IDF compared with when MetS components are in the model
Discussion
In this study, we, for the first time, evaluated the potential applicability of age-specific and sex-specific MetS-S to predict T2DM incidence. We found that an increase in cMetS-S significantly increases the risk of T2DM even after adjusting for potential confounders independent of the MetS components. The cMetS-S offers additional predictive value for T2DM beyond the MetS components. The association of the cMetS-S and future diabetes was established in men and women and both normoglycaemic and prediabetic groups. The strength of the association was notably higher in normoglycaemic individuals.
The traditional criteria of MetS have limited its clinical use due to the inability to define the extent of MetS severity and the risk of unfavourable health outcomes among individuals who do not have MetS. MetS is an assembly of distinct yet related pathophysiological conditions, and the association of MetS with diabetes can be extended to any individual or specific combinations of the MetS components.7 The term preMetS, defined as the presence of one or two MetS components, is associated with an increased risk of T2DM, which is not foreseen in the traditional criteria. Some researchers argue that the preMetS state is a better stage for initiating medical interventions.13 Multiple phenotypes and indices were created using MetS components, such as the triglyceride-glucose index,14 the hyper triglyceridaemic-waist phenotype15 and a classification based on metabolic health and obesity status.16 The multiplicity of these phenotypes and indices based on MetS components results in confusion since they fail to consider all the components involved in MetS. To overcome the limitations, cMetS-S was formulated as a continuous measure of severity based on CFA to consider the weighted contribution of all the MetS components. The cMetS-S enables clinicians to assess the severity of metabolic abnormality, track its change over time and assess the risk of unfavourable outcomes over time in the population regardless of MetS status.
Our findings demonstrate that the risk of T2DM was increased by 72% per 1-SD increase of the cMetS-S after adjusting for age, sex, education, physical activity, family history of diabetes, antihypertensive and lipid-lowering drug use and even the MetS components. We also conducted the analysis based on quantile classification in which the highest quartile of cMetS-S had a 3.84-fold increased risk for future T2DM. It also appeared that cMetS-S has a greater ability to predict T2DM than the MetS definitions based on JIS and IDF, as it offers additional value beyond the components of MetS. Similarly, Gurka et al found the highest quartile of MetS severity Z-score to be associated with a 2.24–5.30 fold increased risk of T2DM in white and black Americans, respectively, after adjusting for age, sex and MetS components.17 The congruent results from the two studies in the literature indicate that MetS severity defined as z-scores can not only help predict future T2DM independent of the MetS components but also simultaneously show the racial differences in attributed risk of T2DM.
Notably, the current and the American study were different in terms of methodology, adjustments and the studied population. The current study has the advantage of using age and sex-specific MetS severity equations to consider the age and sex variation of MetS component values and their attributed risk for T2DM, while in the study by Gurka et al, MetS severity was solely based on variations by sex. Pathophysiological and epidemiological studies have revealed age-dependent variations in the value and attributed risk of MetS components for diabetes.18–21
According to our results, the relationship between cMetS-S with T2DM was non-significantly higher in women. In the current study, the cMetS-S was linked to a 1.65-fold and 1.83-fold higher risk of T2DM in men and women, respectively. The underlying reason may be the sex variation in the baseline MetS components values, their contribution to MetS22 or their associated risk for diabetes.23 24 The difference in body fat composition, sex hormones and lipid metabolism are some potential explanations.25
Comparison within normoglycaemic and prediabetic subgroups showed the association of cMetS-S with T2DM with higher HRs in normoglycaemic group. Normoglycaemic individuals with a metabolic unhealthy profile may already harbour underlying pathophysiological processes that predispose them to a higher risk of developing diabetes, potentially explaining the stronger association observed in the current study. Studies have shown that insulin resistance can occur before the onset of pre-diabetes and diabetes, and it is often more pronounced in normoglycaemic individuals with metabolic abnormalities. This could lead to a stronger association between metabolic unhealthy profile and diabetes in this population. Research has indicated that individuals with a metabolic unhealthy profile, even in the absence of pre-diabetes, may have higher levels of ectopic fat deposition, such as visceral adipose tissue and hepatic fat content. These ectopic fat depots are known to contribute to insulin resistance and impaired glucose metabolism, potentially increasing the risk of developing diabetes in normoglycaemic individuals. Other contributing factors in normoglycaemic individuals are inflammatory markers, endothelial dysfunction due to dyslipidaemia and hypertension and genetic predispositions.26 27 The recent study and a meta-analysis of 61 studies included 71 196 participants with 11 771 T2D events, 412 metabolites from 37 metabolic pathways associated with T2D were reported among these amino acids, carbohydrates and lipid metabolite are most significant. These metabolites would perform significant mediation effects of BMI, dyslipidaemia, inflammatory markers on T2DM.28
To our knowledge, this is the first study to investigate the association of a cMetS-S with future T2DM in normoglycaemic and prediabetic subgroups. The cMetS-S was predictive of T2DM regardless of glycaemic status. The association between cMetS-S and future T2DM in normoglycaemic individuals was significantly higher (HR 2.11 vs HR 1.42, a 1.5-fold difference) than in individuals with pre-diabetes, suggesting the early-onset increased risk of MetS for T2DM independent of glucose metabolism impairment. Previous reports also suggested an association between MetS (excluding FPG as a component) with T2DM in both subgroups of individuals with normal fasting glucose (NFG) and impaired fasting glucose (IFG), which tended to increase incrementally on an increase in the number of MetS traits. The association was non-significantly higher in the NFG group.24 29 A clinical interpretation of the results from the previous studies is not feasible due to different possible combinations of MetS components. In the current study, we took a step further by investigating the MetS-S instead of the number of MetS components. We also used more specific definitions for glucose impairment by classifying the population into prediabetic and normoglycaemic individuals instead of individuals with NFG and IFG. Our findings suggest that the cMetS-S presents a unified quantitative measure of all MetS components that can identify the risk of future T2DM among prediabetic and normoglycaemic individuals in clinical and epidemiologic settings over long-term follow-up. This association is even significantly stronger in normoglycaemic individuals than in those with pre-diabetes. This score overcomes the limitations of the traditional MetS criteria as it was not defined on the outcome prediction basis (T2DM and cardiovascular disease (CVD)) and instead it was developed based on how MetS components cluster together as MetS by assigning weights to each MetS component using CFA and considering the variations of the MetS components among sex and age groups and the severity of MetS on a continuous scale. The study has some limitations. First, since there are variations in the severity of MetS across different racial groups, it is difficult to generalise the findings to other populations. Also to eliminate potential confounding variables that may be more prevalent among older adults, we confined our study to individuals aged between 20 and 60 years. Additionally, we did not consider the potential impact of lifestyle modifications on the severity of MetS over time, as this would have required additional data and resources.
The newly developed age-specific and sex-specific cMetS-S strongly correlates with future T2DM independent of MetS components over long-term period especially in normoglycaemic individuals, indicating the potential clinical applicability of the cMetS-S in detecting those at risk of diabetes before the onset of pre-diabetes. The cMetS-S also provides additional value over MetS components for diabetes risk prediction contrary to the tradition-based MetS definition. The cMetS-S has significant potentials for widespread use in clinical settings as a tool with significant ability to predict future health outcomes, track the MetS severity changes over time and monitor response to treatment in all individuals regardless of MetS status.
Data availability statement
Data are available on reasonable request. Datasets generated during and/or analysed during the current study are not publicly available due to institutional policies but are available from the corresponding author on reasonable request.
Ethics statements
Patient consent for publication
Ethics approval
This study conformed to the ethical guidelines of the Declaration of Helsinki. It was approved by the National Research Council of the Islamic Republic of Iran, the Human Research Review Committee of the Endocrine Research Center of Shahid Beheshti University of Medical Sciences, Tehran, Iran (IR.SBMU.ENDOCRINE.REC.1401.081). All participants provided informed written subjects. It was not appropriate or possible to involve patients or the public in the design, or conduct, or reporting, or dissemination plans of our research.
Acknowledgments
The authors wish to thank the study participants and the TLGS research team for their passionate support.
References
Supplementary materials
Supplementary Data
This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.
Footnotes
Correction notice This article has been corrected since it was published. Author name ‘Mohammadjavad’ was misspelled.
Contributors LM is the guarantor. All authors contributed to writing the manuscript. AA: project administration and editing the final draft. MH: designing, methodology, data analysis and preparation of the original draft; SM: data analysis; SA: writing the original draft; FA: project supervision and editing the final draft; LM: designing, methodology and project administration. All authors reviewed the manuscript.
Funding This study was supported by Grant No 43009357-4 in Shahid Beheshti University of Medical Sciences.
Competing interests None declared.
Patient and public involvement Patients and/or the public were not involved in the design, or conduct, or reporting, or dissemination plans of this research.
Provenance and peer review Not commissioned; externally peer reviewed.
Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.