Introduction

One of the world’s top three causes of mortality are cardiovascular diseases (CVDs). In 2019, of the 55 million deaths worldwide, non-communicable diseases (NCDs) accounted for 71%. Cardiovascular diseases (CVDs) contributed the largest number of deaths worldwide (17.9 million) among NCDs1. According to the Global Burden of CVD, years lived with disability increased from 17.7 million to 34.4 million between 1990 and 20192. The burden of CVDs is particularly high in low- and middle-income countries (LMICs), and because of limited human and economic resources deaths occurs at younger ages than in high-income countries (HIC)3, 4. HICs have reported a decrease in CVD mortality rates over the last decades but 80% of the morbidity burden and 50% of the mortality burden due to CVD is reported to occur in LMICs including the Eastern Mediterranean Region (EMR)5.

Recently, many tools and clinical practice guidelines have been developed for CVD risk prediction, including the Framingham risk score, European Society of Cardiology/European Society of Hypertension guidelines, American College of Cardiology and American Heart Association (ACC/AHA guidelines, Systematic COronary Risk Evaluation (SCORE), and SCORE26,7,8,9,10. These tools have shown benefit for risk prediction, however, they were developed and have mainly been tested in Europe and America11. Because the levels of CVD predictors vary across populations, CVD risk prediction models that have been developed for specific countries cannot automatically be used in other countries12, 13. In the recent decade, two CVD risk equations; WHO and Globorisk were developed that are more specific to each population.

In 2015, Hajifathalian et al. developed the Globorisk CVD risk prediction equations. Globorisk is a new CVD risk calculator that is country-specific. Globorisk has been validated and calibrated using data from 182 countries and its external validity has been examined in different regions14, 15. In 2019, WHO revised previous CVD risk prediction models that were created in 2007. The new WHO CVD risk equations estimated risk for 21 regions of GBD. Also, revised WHO models estimate risk for LMICs16.

WHO and Globorisk CVD risk prediction equations have a laboratory-based and non-laboratory-based version. Laboratory-based models estimate CVD risk using age, sex, SBP, smoking, diabetes, and cholesterol. Non-laboratory-based estimated CVD risk using age, sex, smoking, SBP, and BMI15, 16.

It is important to evaluate the agreement between two 10-year risk prediction models of CVDs, including the WHO model (a region-specific model) and Globorisk (a country-specific model). Some previous studies in Iran determined the 10-year risk of CVDs with Framingham and WHO models17,18,19. However, no study has been conducted to evaluate the agreement and correlation between Globorisk and WHO risk prediction models. Therefore, the present study was conducted with the aim of examining the agreement between the laboratory-based versions of the WHO and Globorisk models, as well as the agreement between the non-laboratory-based versions of these two models in a large Iranian population.

Methods

This cross-sectional study is a part of Fasa cohort study. The details of the Fasa cohort study have been published previously20, 21. Briefly, the study population in the Fasa cohort study was 10,138 persons aged ≥ 35 years old, resident in the Sheshdeh region in Fasa county. Data were collected from 2015 to 2016. Before beginning the study, interviewers were trained to collect data. Trained interviewers collected and registered demographic characteristics, medical history (diabetes, cardiovascular disease, cancer, and etc.), physical activity, smoking status, opium use and alcohol drinking.

To record the anthropometric characteristics, the trained staff measured the height, weight and waist circumference of the participants. Laboratory staff took samples for biochemical tests, including blood sugar and blood lipids including cholesterol, triglycerides (TG), high density lipoprotein (HDL), and low density lipoprotein (LDL). In this study, people who had a history of cardiovascular diseases and stroke were excluded. Also, because the 10-year CVD risk of Globorisk equations was created for people aged ≥ 40 years, participants  74 years old were additionally excluded. Finally, 6796 persons aged 40–74 who had no history of CVD and stroke were included in this study.

CVD risk assessment

In this study, fatal plus non-fatal WHO and Globorisk CVDs risk prediction models were used to determine 10-year CVD risk. For WHO and Globorisk models, cardiovascular outcomes are a 10-year risk of fatal and nonfatal CVD, CHD, or stroke22. WHO and Globorisk equations have two versions; laboratory-based and non-laboratory-based. In the WHO laboratory-based model, CVD risk determined by factors including age, sex, smoking status, diabetes status, SBP and cholesterol (mmol/L). In the WHO non-laboratory-based model, cholesterol and diabetes is replaced by BMI and risk is calculated based on age, sex, smoking status, SBP, and BMI. CVD risk in WHO models is determined according 21 regions of GBD. According to GBD, Iran is in North Africa and Middle East region. Afghanistan, Algeria, Bahrain, Egypt, Iraq, Jordan, Kuwait, Lebanon, Libya, Morocco, Occupied Palestinian territory, Oman, Qatar, Saudi Arabia, Sudan, Syrian Arab Republic, Tunisia, Turkey, United Arab Emirates, and Yemen are included in this region22.

Globorisk is country-specific model. Globorisk also has laboratory-based and non-laboratory-based versions. In the laboratory-based version, risk factors include age, sex, SBP, smoking status, diabetes status, and cholesterol (mmol/L). The office-based version includes age, sex, SBP, smoking status and BMI15.

In the Fasa cohort study, all measurements were done according PERSIAN cohort guidelines23. Smoking status was determined by questions about smoking. For biochemical index measures, each person had to have fasted for 12–16 h. Diabetes was determined by self-reported history of diabetes, treatment or fasting blood sugar ≥ 126 mg/dl. Blood pressure was measured two times with an interval of 5 min after a five-minute rest in sitting position using a standard calibrated sphygmomanometer and the mean of the two measurements was recorded. Weight (kg) and height (m) were measured and BMI calculated (weight in kilograms divided by height in meters squared).

Ethics approval and consent to participate

This study was approved by the Ethics Committee of Jahrom University of Medical Sciences (IR.JUMS.REC.1400.071). The informed consent of participants was obtained if they intended to participate in the study and all information of the participants was collected anonymously. Also, all methods were carried out in accordance with relevant guidelines and regulations.

Statistical analysis

Percentages were used for categorical variables and means ± standard deviations (SDs) were used for continuous variables were used. Chi-square and t-test were used to test for differences across categories. The analysis was done using different methods according to whether the risk is grouped or quantitative. Globorisk and WHO risk scores were divided into three categorized (low, moderate, and high) and percent agreement and kappa statistics were estimated. Kappa statics range is − 1 to 1. However, most is between 0 to 1. According to the suggestion of Landis and Koch, the strength of agreement was interpreted as follows: kappa statistic ≤ 0 be considered ‘poor’ agreement, kappa statistic = 0.00–0.20 be considered ‘slight’ agreement, kappa statistic = 0.21–0.40 be considered ‘fair’ agreement, kappa statistic = 0.41–0.60 be considered ‘moderate’ agreement, kappa statistic = 0.61–0.80 be considered ‘ substantial’, and kappa statistic = 0.81–1.00 be considered ‘almost perfect’24.

In WHO risk models, risk categorized as < 5%, 5% to < 10%, 10% to < 20%, 20% to < 30% and ≥ 30% was defined as very low-, low-, moderate-, high-, and very high-risk16. For comparison WHO and Globorisk models, WHO risk scores was categorized in three groups < 10%, 10–< 20% and ≥ 20% which were considered low-, moderate-, and high-risk. In Globorisk models, risk grouping has not been specified the same way as other CVDs risk prediction models. In the study by Ueda et al. from an Iranian population a threshold ≥ 20% risk was considered a high-risk group15. So, In Globorisk models, risk categorized as < 10%, 10% to < 20%, and ≥ 20% were considered low-, moderate-, and high-risk.

The correlation between WHO and Globorisk models was calculated using correlation coefficients and scatter plots. A scatter plot shows the relationship between two quantitative variables measured for the same people. The general pattern of a scatter plot can be described by the direction, form and strength of the association. The correlation coefficient measures the strength of that association. The value of the correlation coefficient ranges from − 1 to 125. According to the suggestion of Evans, a correlation coefficient between two variables is considered to have a very weak correlation (0.00–< 0.20), weak correlation (0.20–< 0.40), moderate correlation (0.40–< 0.60), strong correlation (0.60–< 0.80), and very strong correlation (0.80–1.0)26.

Bland–Altman plots were presented for determination agreement between risk scores at the individual level. This method was used to assess agreement between two quantitative measurements. A Bland–Altman plot is a scatter plot that plots the difference between two paired of measurements on the Y axis and the mean of the two measurements on the X axis27.

All statistical analyses were conducted using Statistical Package for Social Science (IBM SPSS Statistics for Windows, Version 23.0) and Stata Statistical Software (Stata for windows, Version 14). P-values < 0.05 were considered as significant.

Results

In this study, 6796 participants were included. Table 1 shows a summary of participants’ characteristics. In brief, 53.5% were female, and mean (SD) age of the participants was 51.0 (7.8) years. Also, 52.3% were illiterate. The prevalence of smoking was higher in males than in females, while the prevalence of hypertension and diabetes was higher in females than in males. SBP, DPB, HDL, cholesterol, and BMI was higher in females than males.

Table 1 Reporting of the participants’ characteristics. DBP, diastolic blood pressure; SBP: systolic blood pressure; HDL: high density lipoprotein; Chol: cholesterol; BMI: body mass index, * chi-square test, ** t-test.
Full size table

The mean of laboratory-based WHO risk score was 7.4 ± 5.4 and the mean of laboratory-based Globorisk score was 6.1 ± 6.9. Also, the mean of non-laboratory-based WHO and Globorisk risk scores was 7.2 ± 4.9 and 6.2 ± 6.4, respectively. In the WHO and Globorisk CVDs risk scores, the mean laboratory-based and non-laboratory-based models was higher in males than females (Table 2).

Table 2 The distribution of the 10-year WHO and Globorisk CVD risk scores.
Full size table

Figure 1 showed that in the WHO and Globorisk laboratory-based models 3.3% and 4.8% were in high-risk group, respectively. Also, in the WHO and Globorisk non-laboratory-based models 2.4% and 4.7% were in high-risk group, respectively.

Figure 1

The percentage of the 10-years WHO and Globorisk CVD risks classified (a) laboratory-based models; (b) non-laboratory-based models.

Full size image

Tables 3 and 4 shows the agreement of grouped risk between the WHO and Globorisk risk scores according to the laboratory-based and non-laboratory-based models, respectively. In the total population, agreement between the laboratory-based WHO and Globorisk risk groups was 91.4% (kappa = 0.75). The agreement was almost the same in females (kappa = 0.74) and males (kappa = 0.76).

Table 3 Agreement between the laboratory-based WHO and Globorisk cardiovascular disease risk scores models according to the grouped risk.
Full size table
Table 4 Agreement between the non-laboratory-based WHO and Globorisk cardiovascular disease risk scores models according to the grouped risk.
Full size table

In the total population, agreement between the non-laboratory-based WHO and Globorisk risk groups was 92.6% (kappa = 0.78). The agreement was better in males than females (kappa: 0.82 vs. 0.73).

Correlation coefficients between the WHO and Globorisk risk scores are shown in Table 5. In the overall population and in males and females separately, the correlations between laboratory-based WHO and Globorisk CVD risk scores were very strong. There were significant positive correlations between the laboratory-based WHO and Globorisk risk scores in males and females. Also, in the overall population and in males and females separately, the correlations between non-laboratory-based WHO and Globorisk CVD risk scores were very strong. There were significant positive correlations between the non-laboratory-based WHO and Globorisk risk scores in males and females (Figs. 2 and 3).

Table 5 Correlation coefficients between WHO and Globorisk cardiovascular disease risk scores models according to laboratory-based and non-laboratory-based model.
Full size table
Figure 2
figure 2

Scatter plots and Bland–Altman plots of individual-level risk of laboratory-based WHO and Globorisk cardiovascular disease risk scores models. (a) scatter plot for total population, (b) scatter plot for males, (c) scatter plot for females, (d) Bland–Altman plot for total population, (e) Bland–Altman plot for males, (f) Bland–Altman plot for females.

Full size image
Figure 3
figure 3

Scatter plots and Bland–Altman plots of individual-level risk of non-laboratory-based WHO and Globorisk cardiovascular disease risk scores models. (a) scatter plot for total population, (b) scatter plot for males, (c) scatter plot for females, (d) Bland–Altman plot for total population, (e) Bland–Altman plot for males, (f) Bland–Altman plot for females.

Full size image

Bland–Altman plots which show the agreement between the Globorisk and WHO risk scores on the individual level according to laboratory-based and non-laboratory-based models are shown in Figs. 2 and 3. In laboratory-based models, the limit of agreement was better for males than females (2.1 to − 4.2% vs. 4.3 to − 7.3%). Also, in non-laboratory-based models, the limit of agreement was better for males than females (2.9 to − 4.0% vs. 3.2 to − 6.1%) (Table 5).

Discussion

In this study, two CVD risk prediction models were assessed including WHO risk model which is region-specific and Globorisk that is country-specific. To our knowledge, this study is the first study to evaluate and compare laboratory-based WHO and Globorisk and also non-laboratory-based versions of WHO and Globorisk separately in a large population. The agreement and correlation beetween these modeles were determined according to laboratory-based and non-laboratory-based models, separately.

Recently for this population, the agreement between laboratory-based and office-based Globorisk and also the agreement between laboratory-based and non-laboratory-based WHO CVD risk were estimated separately28, 29. The results showed that there was substantial agreement between WHO and Globorisk laboratory-based risk grouping and also between WHO and Globorisk non-laboratory-based risk grouping. Very strong positive correlations were observed between the WHO and Globorisk laboratory-based risk scores and also between the WHO and Globorisk non-laboratory-based risk scores.

The results showed substantial agreement between Globorisk and WHO risk grouping. Other studies determined agreement between CVDs risk models. Rezaei et al. assessed the agreement between laboratory-based and non-laboratory-based versions of the Framingham risk score and WHO risk equations in the Pars cohort and found lower agreement between the two than in the current analysis18, 19. In Bangladesh the agreement between non-laboratory-based Framingham risk score, Globorisk and WHO/ISH models was estimated, and was also lower than in the current study30. However, we are not aware of any studies that have assessed the agreement between the laboratory-based versions of the WHO and Globorisk models and the non-laboratory-based versions of the two models separately.

The results from this study showed that 3.3% when using the WHO laboratory-based model and 2.4% when using the WHO non-laboratory-based model were in the high risk group (≥ 20%). In another study, 5.0% and 4.7% of participants were in the the high-risk group based on the WHO laboratory-based model and the WHO non-laboratory-based model, respectively19. Although both studies were conducted in Iran using the revised WHO version, due to the difference in the population structure and risk factors, the risk classification differed slightly in the two studies.

Another CVD risk assessment tool that was used was the Globorisk equations. The results showed that 4.8% and 4.5% in the Globorisk laboratory- based model and non-laboratory-based model were in the high- risk group (≥ 20%), respectively. In Malaysia, fewer people were in the low-risk group, likely reflecting differences in the distribution of risk factors in that population compared to the current Iranian study. In the Globorisk laboratory-based model, 45% had a risk of < 10%, 14% ≥ 20%, and 11% ≥ 30%, and in the non-laboratory-based model 51.1% had a risk of < 10%, 11% ≥ 20%, and 4.9% ≥ 30%31.

In this study, the agreement between the laboratory-based WHO and Globorisk groups risk was substantial both overall and in females and males separately. The non-laboratory-based models is a highly effective screening tool at lower cost in settings with limited resources where laboratory tests are unreachable or expensive32. In Bangladesh, the concordance between non-laboratory-based World Health Organization/International Society of Hypertension (WHO/ISH) and Globorisk CVD risk prediction models was 0.3730, which was very different when compared to our study, where agreement was substantial. One explanation could be that the study from Bangladesh used the 2007 version of WHO/ISH risk model, while we used the revised WHO risk model from 201916. In addition, differences in risk factor levels between the populations may also have contributed to differences in agreement.

In this study, there was a very strong positive correlation between Globorisk and WHO risk scores. On the other hand, in both laboratoty-based and non-laboratory-based models, there were more high-risk individuals in the Globorisk model than WHO model. Thus, for the important function of identifying people at elevated risk for interventions the two scores perform differently. So, Bland–Altman plots were used for determination agreement between risk scores. The results showed that in laboratory-based and also in non-laboratory-based models the agreement between the two risk scores was better among males than females. In a study carried out in Iran, the agreement between laboratory-based and BMI-based Framingham models was examined with Bland–Altman plot and showed that the agreement between the two risk scores was better among the yonger females and males and was wider among the older males18. A study in Peru showed that there were no substantial differences between the mean CVD risk computed with the laboratory-based model and the non-laboratory-based model33. It is important to mention that acceptable limits should be defined based on various factors, including clinical and biological factors, as well as other considerations that have already been defined27, 34, but we had no information about these factors.

Totally, this study showed that for men the limit of agreement was appropriate especially in low-risk scores. For females, although the limit of agreement was wider statistically in the laboratory and non-laboratory-based models, it seems that these differences are not very important clinically. The scatter plots and Bland–Altman plots demonstrate systematic differences between the two scores that vary according to the level of risk and are more problematic in females than in males. Almost all females at the lower end of moderate risk by WHO are classified as low risk by Globorisk. The majority of females at the lower end of high risk on Globorisk models are only moderate risk according to WHO models. This is completely consistent with the limited overlap of these categories reported in Tables 3 and 4. So, it may be necessary to modify the cut points of risk groups for the two models. However, for this population, modification of the risk grouping cut-point was not the objective of this study.

The Fasa PERSIAN cohort study has various ethnicities including Fars, Arab, and Turk that have different behaviors and socioeconomic status. Differences in ethnicity, socioeconomic status, and genetics of people in various geographic regions may also contribute to differences in CVD risk estimates. It is important to mention that most of the studies determined the agreement of the Framingham model with other models, and until the writing of this manuscript, the authors have not seen other papers that evaluated the agreement and correlation of the Globorisk model with the WHO updated version from 2019.

Strengths and limitations

The main strength of this study was the large sample and the use of carefully collected data from a population-based study. The major limitation of the current study is the cross-sectional design using baseline data of a cohort study, to validate both laboratory-based and non-laboratory-based versions of the Globorisk and WHO risk models. Further studies using a longitudinal design and with longer follow-up is needed. Another limitation is that we only determined the agreement between the models. Further analyses using CVD incidence will allow us to evaluate how well these scores actually predict CVD outcomes.

Conclusion

In this study, an acceptable agreement was observed between the classified risk of laboratory-based WHO and Globorisk models and the non-laboratory-based WHO and Globorisk models. Also, strong positive correlations were observed between the risk scores of these models. However, it should be noted that although a good agreement was observed between the models, a large number of people who were classified as moderate risk by the WHO models were low-risk and high-risk according to Globorisk models. Therefore, it may be necessary to modify the cut points of risk groups for the two models. For WHO and Globorisk models, the ability to determine risk groups must be confirmed separately. Further longitudinal studies with 10-year follow-up are needed to assess actual fatal plus non-fatal CVD events.