Article Text

Download PDFPDF

Original research
Use of risk chart algorithms for the identification of psoriatic arthritis patients at high risk for cardiovascular disease: findings derived from the project CARMA cohort after a 7.5-year follow-up period
  1. Jessica Polo y La Borda1,
  2. Santos Castañeda2,
  3. Elena Heras-Recuero3,
  4. Fernando Sánchez-Alonso4,
  5. Zulema Plaza5,
  6. Carmen García Gómez6,
  7. Ivan Ferraz-Amaro7,8,
  8. Jesús Tomás Sanchez-Costa9,
  9. Olga Carmen Sánchez-González10,
  10. Ana Isabel Turrión-Nieves11,
  11. Ana Perez-Alcalá12,
  12. Carolina Pérez-García13,
  13. Carlos González-Juanatey14,15,
  14. Javier Llorca16,17 and
  15. Miguel Angel Gonzalez-Gay3,18
  16. On behalf of the CARMA Project Collaborative Group
    1. 1 Rheumatology, Hospital Universitario Rey Juan Carlos, Madrid, Spain
    2. 2 Rheumatology Department, Hospital Universitario La Princesa, IIS-Princesa, Madrid, Spain
    3. 3 Rheumatology, ISS Fundacion Jimenez Diaz, Madrid, Spain
    4. 4 Research Unit of Spanish Society of Rheumatology, Madrid, Spain
    5. 5 Research Unit, Spanish Society of Rheumatology, Madrid, Spain
    6. 6 Rheumatology, Consorci Sanitari de Terrassa, Barcelona, Spain
    7. 7 Rheumatology Unit, Hospital Universitario de Canarias, Santa Cruz de Tenerife, Spain
    8. 8 Internal Medicine, Universidad de La Laguna (ULL), Tenerife, Spain
    9. 9 Rheumatology, Research Unit, Spanish Society of Rheumatology, Madrid, Spain
    10. 10 Rheumatology, Hospital Universitario del Sureste, Madrid, Spain
    11. 11 Rheumatology Department, Hospital Complejo Asistencial Universitario de Salamanca, Salamanca, Spain
    12. 12 Rheumatology, Hospital Universitario Príncipe Asturias, Alcalá de Henares, Madrid, Spain
    13. 13 Rheumatology Department, Consorci Parc de Salut MAR de Barcelona, Barcelona, Spain
    14. 14 Cardiology, Hospital Universitario Lucus Augusti, Lugo, Spain
    15. 15 Biodiscovery HULA-USC Group, Instituto de Investigación Sanitaria de Santiago de Compostela IDIS, Santiago de Compostela, Spain
    16. 16 Epidemiology, CIBER Epidemiología y Salud Pública (CIBERESP) and Department of Medical and Surgical Sciences, University of Cantabria, Santander, Spain
    17. 17 CARMA Project Collaborative Group, Madrid, Spain
    18. 18 Department of Medicine and Psychiatry, Medicine, University of Cantabria, Santander, Spain
    1. Correspondence to Dr Miguel Angel Gonzalez-Gay; miguelaggay{at}hotmail.com

    Abstract

    Objective To assess the predictive value of four cardiovascular (CV) risk algorithms for identifying high-risk psoriatic arthritis (PsA) patients.

    Methods Evaluation of patients with PsA enrolled in the Spanish prospective project CARdiovascular in RheuMAtology. Baseline data of 669 PsA patients with no history of CV events at the baseline visit, who were followed in rheumatology outpatient clinics at tertiary centres for 7.5 years, were retrospectively analysed to test the performance of the Systematic Coronary Risk Assessment (SCORE), the modified version (mSCORE) European Alliance of Rheumatology Associations (EULAR) 2015/2016, the SCORE2 algorithm (the updated and improved version of SCORE) and the QRESEARCH risk estimator version 3 (QRISK3).

    Results Over 4790 years of follow-up, there were 34 CV events, resulting in a linearised rate of 7.10 per 1000 person-years (95% CI 4.92 to 9.92). The four CV risk scales showed strong correlations and all showed significant associations with CV events (p<0.001). SCORE, mSCORE EULAR 2015/2016 and QRISK3 effectively differentiated between low and high CV risk patients, although the cumulative rate of CV events observed over 7.5 years was lower than expected based on the frequency predicted by these risk scales. Additionally, model improvement was observed when combining QRISK3 with any other scale, particularly the combination of QRISK3 and SCORE2, which yielded the lowest Akaike information criterion (411.15) and Bayesian information criterion (420.10), making it the best predictive model.

    Conclusions Risk chart algorithms are very useful for discriminating PsA at low and high CV risk. An integrated model featuring QRISK3 and SCORE2 yielded the optimal synergy of QRISK3’s discrimination ability and SCORE2’s calibration accuracy.

    • Arthritis, Psoriatic
    • Atherosclerosis
    • Cardiovascular Diseases
    • Mortality

    Data availability statement

    No data are available.

    http://creativecommons.org/licenses/by-nc/4.0/

    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

    Request Permissions

    If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

    WHAT IS ALREADY KNOWN ON THIS TOPIC

    • Cardiovascular (CV) mortality and events are more prevalent in patients with psoriatic arthritis (PsA) compared with the general population, resembling rates seen in rheumatoid arthritis patients.

    • Several studies have highlighted that Systematic Coronary Risk Assessment (SCORE) risk charts consistently underestimate CV risk in patients with PsA, emphasising the imperative need to develop new algorithms that align with modern treatments and effective management of the disease.

    WHAT THIS STUDY ADDS

    • This study compares, for the first time, the strengths and limitations of four distinct scales (SCORE, modified SCORE European Alliance of Rheumatology Associations 2015/2016, SCORE2 and QRISK3) to evaluate CV risk in patients with PsA, followed prospectively in rheumatology units over 7.5 years.

    • Combining QRISK3 with any other risk chart algorithm improves the ability to identify PsA patients at increased CV risk. In particular, a model integrating QRISK3 and SCORE2 emerged as the most effective approach to discriminate between PsA patients with low and high CV risk.

    HOW THIS STUDY MIGHT AFFECT RESEARCH, PRACTICE OR POLICY

    • Understanding the real risk of CV events in patients with PsA treated in referral rheumatology units is crucial to formulating effective health strategies and reducing CV complications.

    • In assessing CV risk in patients with PsA, the combined application of CV risk scales enhances their effectiveness compared with using each scale individually.

    • Early assessment of CV risk enables timely detection of CV complications and the implementation of suitable measures to prevent the occurrence of future CV events in patients with PsA.

    Introduction

    Psoriatic arthritis (PsA) is a chronic inflammatory condition with a broad range of musculoskeletal manifestations, including peripheral arthritis, spondylitis, dactylitis and enthesitis, often linked with cutaneous psoriasis.1 Patients with PsA may also present other extra-articular manifestations such as uveitis or inflammatory bowel disease.2

    PsA correlates with a higher incidence and prevalence of cardiometabolic comorbidities compared with the general population.3 4 The elevated risk of cardiovascular (CV) disease in patients with PsA is due in part to a higher frequency of classic CV risk factors and metabolic syndrome.3 4 Chronic inflammation also plays an important role, leading to a higher frequency of endothelial dysfunction, early-stage atherogenesis and subclinical atherosclerotic disease, even in the absence of traditional CV risk factors.5 6 Additionally, certain medications, such as glucocorticoids and non-steroidal anti-inflammatory drugs, and other treatments such as conventional synthetic disease-modifying antirheumatic drugs (DMARDs), biological DMARDs and targeted synthetic DMARDs, may have an impact on CV risk, whether positive or negative, in these patients.7

    The high morbidity and mortality due to CV events among patients with PsA compared with the general population,8–10 together with that observed in other inflammatory arthritis such as rheumatoid arthritis (RA), led to the European Alliance of Rheumatology Associations (EULAR) an advocate for periodic CV disease risk assessments at least every 5 years for these individuals.7 11

    Taking into account the aforementioned factors, it is imperative to optimise prevention strategies, emphasising the crucial need to identify PsA patients at high risk of CV events. However, there is a notable absence of specific indices designed to establish the CV risk of individuals with rheumatic diseases. In this sense, the Systematic Coronary Risk Assessment (SCORE) has been recommended to evaluate CV risk in the European population.12 13 However, it is important to note that this algorithm, originally designed for the general population, estimates the 10-year risk of death from CV disease based on conventional CV risk factors (age, sex, systolic blood pressure, total cholesterol and smoking).13

    Based on growing concern about the increased CV risk observed in patients with inflammatory arthritis, the EULAR Study Group on CV Diseases in Rheumatic Diseases introduced and endorsed the modified SCORE (mSCORE) EULAR 2015/2016 specifically for patients with RA and other inflammatory arthropathies. This adaptation of the SCORE involved a multiplication factor of 1.5 to compensate for the excess CV mortality observed in these patients compared with the general population.7 11

    In 2021, a notable evolution occurred in the CV risk assessment landscape with the introduction of the SCORE2 predictive model, which marked an update to the existing SCORE CV risk algorithm.14 SCORE2 underwent meticulous calibration and validation processes, specifically designed to predict the 10-year risk of first-time CV disease events in European populations. SCORE2 differs from the traditional SCORE in several aspects. Specifically, it provides risk estimates that encompass fatal and non-fatal CV disease events. Moreover, SCORE2 incorporates consideration of the competing risks posed by non-CV deaths, a feature absent in the original SCORE framework.

    SCORE2’s recalibration strategy further differentiates it as it now serves four distinct European regions, unlike the two levels of regional stratification offered by SCORE.14 Additionally, the SCORE2 recalibration is systematically aligned with contemporary CV disease rates, which differentiated from the original SCORE model that was based on data collected before 1986. This transformative update improves the accuracy and relevance of CV risk assessment, aligning it more closely with the current dynamics of CV disease in the European population.

    An alternative tool in this scenario is the QRESEARCH risk estimator version 3 (QRISK3), designed to measure the probability of suffering fatal and non-fatal CV events over a 10-year period.15 QRISK3, which emerged from collaboration between clinicians and academics associated with the UK National Health Service, has gained recognition as a reliable 10-year predictor of CV disease in British and international cohorts.15 By distinguishing itself, QRISK3 goes beyond traditional CV risk factors, encompassing diabetes mellitus and chronic kidney disease in its considerations. Additionally, it provides the option to take into account the presence of RA or systemic lupus erythematosus in the patient profile.15

    Concerns have been raised regarding the possible limitation of indices developed for the general population in accurately predicting CV risk in people with rheumatic diseases.10 This discrepancy was highlighted in studies in patients with RA16 and ankylosing spondylitis (AS).17 In this context, a study focused on patients with PsA, which used carotid ultrasound to identify subclinical atherosclerosis, highlighted that the SCORE and similar indices may underestimate the actual CV risk of these patients.18 Similarly, a recent study highlighted that the SCORE algorithm tends to underestimate CV risk, and the mSCORE EULAR 2015/2016 does not provide any improvement in the prediction of CV disease in patients with PsA.19

    In line with the above, a study in patients with RA disclosed that employing both the mSCORE EULAR 2015/2016 index and the QRISK3 index in combination enabled the identification of individuals exhibiting an increased likelihood of developing atheromatous plaques in the carotid artery, indicative of elevated CV risk.20 Similarly, in keeping with this perspective, another study within the RA patient population indicated that QRISK3 displayed a heightened sensitivity in recognising individuals with a high or very high CV risk compared with the mSCORE EULAR 2015/2016.21

    Taking all these considerations into account, we aimed to determine the predictive value of SCORE, mSCORE, EULAR 2015/20167 SCORE2 and QRISK3 risk chart algorithms to identify PsA patients at high risk of CV disease. Furthermore, we sought to explore whether combining two risk chart algorithms could improve the identification of PsA patients with high-risk CV disease. To achieve this objective, we evaluated patients enrolled in the Spanish prospective CARdiovascular in RheuMAtology (CARMA) project. This project included a cohort of patients with PsA prospectively, followed in rheumatology outpatient clinics over a period of 7.5 years.

    Patients and methods

    Population

    The CARMA project is a prospective cohort study with the primary goal of delineating the CV disease risk profile in individuals with chronic inflammatory rheumatic diseases over a 10-year observation period. The study enrolled patients diagnosed with AS, PsA and RA, along with a comparative cohort of individuals without inflammatory diseases. The recruitment process occurred across 67 Spanish hospitals from July 2010 to January 2012.22 This report focuses specifically on data from patients with PsA at 7.5 years from the start of the study.

    The initial (baseline) recruitment included 721 patients with PsA.22 As previously reported,22 all met Moll and Wright’s criteria for PsA.23 The reason for using these criteria was that at the time of the first discussion of the project in 2007, the 2006 Classification Criteria for Psoriatic Arthritis (CASPAR)24 were not widespread or used in all centres included in the project.

    Data from a set of 669 Spanish patients diagnosed with PsA in the CARMA project without a history of CV events at the time of recruitment were evaluated.22 In this sense, 52 patients with a history of CV events (ischaemic heart disease, stroke, peripheral arterial disease or heart failure) before the start of the study (before the baseline visit) were excluded. Most patients (n=552; 82.5%) continued to periodically attend outpatient clinics at each centre included in the project. At 7.5 years, information on all patients included in the initial cohort was assessed. It was obtained by consulting their medical records or by calling patients or family members directly. When it was not available, we requested information from the Spanish National Statistics Institute, which is the institution in charge of keeping records of vital events, such as births and deaths in Spain. This allowed us to determine information on 669 patients with PsA who at the time of the initial visit had no history of CV events.

    Information on this cohort is included in online supplemental table 1.

    Supplemental material

    The retrospective examination involved testing of four risk algorithms, namely SCORE, mSCORE EULAR 2015/2016, SCORE2 and QRISK3, using baseline data acquired at the start of the patient’s participation in the CARMA project.

    Variable specifications and operative definitions

    CV events and CV mortality

    The analysis focused on the incidence of CV events and CV mortality. The spectrum of CV events encompassed ischaemic heart disease, stroke, peripheral arterial disease and heart failure. The operational definitions outlining the parameters of the variables under analysis were detailed in a separate report.22

    Calculation of CV risk algorithms

    The determination of SCORE and SCORE2 followed the previously described methodologies.12 14 The SCORE calculation involved factors such as age, sex, smoking, systolic blood pressure and total cholesterol. On the other hand, SCORE2 was obtained using age, sex, smoking, systolic blood pressure and non-HDL cholesterol.

    The SCORE2 calculation was executed using the ‘score2risk’ command for Stata, accessible at http://www.phpc.cam.ac.uk/ceu/erfc/programs//.14 25 The SCORE calculation followed the methodology described by Conroy et al.12 SCORE evaluates the 10-year risk of death from CV disease. However, recognising that combining morbidity with mortality provides a more comprehensive perspective on the total burden of atherosclerotic CV disease (ASCVD), SCORE2 now estimates an individual’s 10-year risk of fatal and non-fatal CV disease events for people from 40 to 69 years old. For people aged ≥70 years, the SCORE2-OP (older people) algorithm estimates the 5-year and 10-year risk of fatal and non-fatal CV disease events.

    The determination of mSCORE EULAR 2015/2016 adhered to the updated recommendations from the EULAR group in 2015/2016.6 7

    As for QRISK3, its calculation was aligned with previously described procedures.15 This encompassed the analysis of several variables, including conventional CV risk factors, diabetes mellitus, chronic kidney disease, emerging CV risk factors and other inflammatory rheumatic diseases. QRISK3 measures an individual’s risk of having a heart attack or stroke over the next 10 years. It reflects the average risk of individuals who share the same risk factors as those entered for that specific person.15 The algorithm establishes a threshold of 10%, classifying CV risk as high-very high at or above this level (≥10%), and below that level, classifying it as low-moderate CV risk.

    Statistical analysis

    Patient demographic and clinical variables were presented as mean and SD, or median and IQR for skewed quantitative variables, and as absolute numbers and relative frequencies (%) for qualitative variables. The correlation between the CV risk scales was quantified using the Pearson correlation coefficient.

    Survival analysis was approached in two ways:

    1. In the first step, we defined failure as the combination of a CV event or death from any cause. Patients without a CV event and still alive at the end of the follow-up period were considered censored. Follow-up time was the duration until the occurrence of the first CV event or death for failures and until the last visit for censored patients. The discriminative effectiveness of any CV scale in identifying failures in this analysis was determined using Cox regression.

    2. Subsequently, we designated any CV event as a failure and death resulting from any non-CV cause was treated as a competing risk. Follow-up time mirrored that of analysis (a). The correlation between the CV scale and CV events was explored using the method proposed by Fine and Gray. No adjustments were made to any model for age or sex, given their inclusion in the CV scales.

    Results from the survival analysis were portrayed as HR, accompanied by their 95% CI and two-tailed p values. Each of the survival analyses (a)–(b) involved the creation of seven models. Models 1–4 incorporated an individual CV scale each. Since SCORE, mSCORE EULAR 2015/2016 and SCORE2 are variants of each other, we abstained from examining their combination. Subsequently, to assess if augmenting QRISK3 with an additional CV risk scale would enhance the model, models 5–7 were constructed. Therefore, models 5–7 were confined to combinations of QRISK3 with any scale from the ‘SCORE family’. Due to the strong correlation between QRISK3 and any other CV scale, the ‘residual method’ was employed. Although we exemplify it here for the QRISK3–SCORE combination, the method is analogous for QRISK3–mSCORE EULAR 2015/2016 and QRISK3–SCORE2 combinations.

    The residual method aimed to dissect SCORE into two components: one linearly dependent on QRISK3 (‘linear prediction’) and another independent of QRISK3 (‘residuals’). To achieve this, we executed a linear regression model with SCORE as the Y variable and QRISK3 as the regressor. The residuals derived from this regression represented the portion of SCORE independent of QRISK3. Subsequently, the survival analysis incorporated both QRISK3 and the residual of SCORE as regressors.

    For a comprehensive comparison of all seven models generated in each survival analysis, we computed the Akaike information criterion (AIC) and the Bayesian information criterion (BIC). Both criteria share a common principle: they penalise the introduction of new variables, thereby favouring parsimonious models unless the inclusion of a new variable significantly enhances the model likelihood. In essence, models with lower AIC or BIC are deemed more favourable.

    All statistical analyses were conducted by using Stata V.18/SE software (StataCorp). Statistical significance was considered with p values<0.05.

    Results

    Correlation among CV risk scales

    The correlation analysis revealed a robust association among all four CV risk scales (table 1). It is noteworthy that SCORE and mSCORE EULAR 2015/2016, differing solely by a multiplying factor of 1.5, essentially function as a simple change of units. Consequently, their correlation with each other attains a perfect value of 1.

    Table 1

    Correlation between CV risk scales in patients with psoriatic arthritis from the CARMA cohort

    Evaluation of CV risk scales: compound events (CV event or death by any cause)

    A total of 41 compound events, encompassing CV events or deaths by any cause, were observed during the 4790 years of follow-up. The linearised rate was calculated at 8.56 per 1000 person-years, with a 95% CI of 6.14 to 11.61.

    The prediction of CV events or deaths by any cause is detailed in table 2 and illustrated in figure 1A–D. In comparing various scales, two essential characteristics must be considered: discrimination and calibration.

    Table 2

    CV risk scales performance with the compound event death by any cause or CV event in the psoriatic arthritis CARMA cohort after 7.5 years of follow-up

    Figure 1

    Cardiovascular event or death from any cause. Prediction using: (A) QRISK3, (B) SCORE, (C) mSCORE EULAR 2015/2016, (D) SCORE2. Notice that numerical data in the vertical axis ranges from 0 to 0.15 for subgraphics A–C, and from 0 to 0.35 for subgraphic D. EULAR, European Alliance of Rheumatology Associations; mSCORE, modified SCORE; QRISK3, QRESEARCH risk estimator version 3; SCORE, Systematic Coronary Risk Assessment.

    Discrimination, which denotes the ability to distinguish between individuals at high and low CV risk, can be measured through the risk indices provided in table 2. A higher risk index means superior discrimination. Additionally, discrimination can be assessed visually by examining the spacing between lines in figure 1A–D. A wider separation indicates greater discrimination.

    Calibration involves verifying the accuracy of the predicted risk. For example, if a scale predicts a 10% CV event risk over 10 years for a group of patients, it would be considered well calibrated if approximately 10% of the patients in that group experienced a CV event over the 10 years period of follow-up. Conversely, it would be miscalibrated if the rate of CV events observed at follow-up differed significantly, such as occurring in 5% or 15% of patients.

    Calibration evaluation can be performed by examining figure 1A–D and comparing each line to the background grid. For example, if the red line (representing individuals with a predicted risk of 10%) aligns with the 10% line on the grid at the end of the follow-up period, the scale is considered well calibrated. Conversely, if the red line deviates significantly from the 10% line on the grid, the scale is considered poorly calibrated.

    All four scales demonstrated a significant association with events (p<0.001) as presented in table 2. The HRs were 1.07 (95% CI 1.05 to 1.08) for QRISK3, 1.06 (95% CI 1.04 to 1.09) for SCORE, 1.04 (95% CI 1.03 to 1.06) for mSCORE EULAR 2015/2016 and 1.17 (95% CI 1.12 to 1.23) for SCORE2.

    QRISK3 effectively differentiated between individuals with lower and higher risk, as illustrated in figure 1A. However, the percentage of events accumulated over 7.5 years was notably lower than expected based on the risk established by QRISK3. For example, individuals with a QRISK3 of 20% in 10 years experienced only about 9% of events after 7.5 years, as depicted by the black line in figure 1A.

    SCORE also effectively differentiated between individuals with lower and higher risk, as illustrated in figure 1B. However, after 7.5 years of follow-up, the cumulative risk for individuals with SCORE=5% exceeded the expected percentage of events. Additionally, for those with SCORE=10%, 15% or 20%, the percentage of events by year 7.5 of follow-up was lower than predicted by this index.

    mSCORE EULAR 2015/2016 demonstrated comparable discrimination ability to SCORE, although the predicted percentage of events was notably higher than the actual percentage. For example, individuals with mSCORE EULAR 2015/2016=20% experienced approximately 8% of events (figure 1C, black line).

    SCORE2 demonstrated robust discrimination ability in terms of CV risk, as evident from the well-separated lines in figure 1D. For individuals with SCORE2=5%, 10% or 15%, the percentage of events by year 7.5 of follow-up aligned closely with predictions (blue, red and green lines in figure 1D, respectively). However, SCORE2 underestimated the risk for individuals with SCORE2=20% black line in figure 1D, where the cumulative hazard reached 33% in 7.5 years).

    Combining QRISK3 with any other scale enhanced the model, as indicated in table 2 (AIC and BIC for any model with two scales were lower than for any model with just one scale). The combination of QRISK3 and SCORE2 achieved the lowest AIC and BIC, making it the most effective predictive model.

    Performance assessment of CV risk scales considering CV events (with death by non-CV causes as a competitive event)

    Table 3 and figure 2A–D illustrate the prediction of CV events when death by non-CV causes was treated as a competitive event. Over 4790 years of follow-up, there were 34 CV events, resulting in a linearised rate of 7.10 per 1000 person-years (95% CI 4.92 to 9.92). All four scales demonstrated discriminative power, as depicted in figure 2A–D, with p values lower than 0.001 in table 3. Notably, SCORE2 exhibited significantly lower AIC and BIC values. Consequently, among the four scales, SCORE2 stands out as the preferred choice. The discriminative ability of QRISK3 improved when combined with any other scale, evident in lower AIC and BIC values for combined models compared with the QRISK3-alone model. A model that integrated QRISK3 and SCORE2 demonstrated the lowest AIC (411.15) and BIC (420.10) values, making this combination the preferable choice.

    Table 3

    CV risk scales performance with the CV event in the psoriatic arthritis CARMA cohort after 7.5 years of follow-up

    Figure 2

    Cardiovascular (CV) event (death from non-CV causes as a competitive event). Prediction using: (A) QRISK3, (B) SCORE, (C) mSCORE EULAR 2015/2016, (D) SCORE2. Notice that numerical data in the vertical axis ranges from 0 to 0.15 for subgraphics A–C, and from 0 to 0.30 for subgraphic D. EULAR, European Alliance of Rheumatology Associations; mSCORE, modified SCORE; QRISK3, QRESEARCH risk estimator version 3; SCORE, Systematic Coronary Risk Assessment.

    Discussion

    In this study, we provide information derived from a large cohort of patients with PsA participating in the Spanish prospective CARMA project, focusing on the CV outcomes of people with inflammatory arthritis. This cohort analysis incorporates data from 669 patients with PsA collected over a period of 7.5 years after enrollment. Our research involved the evaluation of four CV risk algorithms. In addition to several versions of the European SCORE, the algorithm traditionally used to determine CV risk in the Spanish population, we also analysed in our cohort the QRISK3, used in the UK to evaluate CV risk.

    The study findings suggest that SCORE, mSCORE EULAR 2015/2016, SCORE2 and QRISK3 are useful to assess CV risk PsA patients. In each instance, these risk algorithms aptly differentiated between PsA patients with high and low CV risk. The performance of SCORE2 was better than the other equations. It was evident when SCORE2 was compared with SCORE. This observation suggests that non-HDL cholesterol might play a crucial role in influencing CV events in PsA, given the prevalent metabolic risk in these patients. However, we feel that this insight may be extrapolated to the general population.

    Despite strong correlations among the four scales, each failed to encompass all the necessary information about CV risk. Therefore, combining two algorithms could offer a more comprehensive understanding of the CV risk profile in PsA patients. In this regard, a model that integrated QRISK3 and SCORE2 yielded the best approach to identify patients with PsA at high risk for CV events.

    There is limited information on the usefulness of CV risk algorithms in patients with PsA. In this context, Shen et al conducted an assessment of four CV risk scales: Framingham Risk Score (FRS), QRISK2, the 10-year ASCVD risk algorithm proposed by the American College of Cardiology and the American Heart Association, and the mSCORE (SCORE multiplied by a factor of 1.5) as proposed by the EULAR group.26 The study included 146 Chinese PsA patients, with 42.5% displaying carotid atheromatous plaque. These patients with carotid plaques were notably older, with higher systolic blood pressure values and higher LDL-cholesterol levels. FRS, QRISK2, ASCVD and mSCORE algorithms had moderate discrimination in recognising PsA patients with and without atheromatous plaque in the carotid artery. Notably, the predetermined high-risk cut-off values (FRS>10%, QRISK2>20%, SCORE>5% and ASCVD>7.5%) for all indices underestimated the risk by failing to identify atheromatous plaques.26 These findings differ from our results, likely due to differences in the number of patients included in the study and possibly to differences in the genetic background within both studies. More importantly, Shen et al employed carotid atheromatous plaque as a surrogate marker for CV risk, whereas our study focused on evaluating the scales’ performance in identifying CV events.

    Navarini et al 19 conducted another study in patients diagnosed with PsA according to the CASPAR criteria.24 These authors performed a retrospective analysis using prospectively collected data from two Italian rheumatology centres. The data set, comprising 155 PsA patients with no history of CV events at baseline in November 2007, was reviewed and analysis was performed in November 2017. Over the follow-up period, 21 patients with PsA experienced a CV event. The study evaluated the performance of five risk algorithms: FRS, SCORE, QRISK2, Reynolds risk scores (RRS) and the Italian individual index Progetto CUORE. The findings revealed that all five algorithms demonstrated relatively good discrimination between PsA patients with and without CV events, with areas under the ROC curves ranging from 0.7183 for RRS to 0.8660 for QRISK2. However, the calibration of these algorithms ranged from poor to moderate. Notably, CUORE, SCORE and RRS models exhibited a poor fit, as the distribution of observed events significantly differed from the predicted values. In contrast, FRS and QRISK2 showed a more aligned distribution. Moreover, the adaptation according to EULAR recommendations performing multiplication by the factor of 1.5 of CV risk algorithms in patients with inflammatory arthritis did not enhance the algorithms’ discriminative ability or calibration.19 27

    There some possible explanations for the apparent contradiction between the results by Navarini et al and the results of the CARMA cohort after 7.5 years of follow-up. In this respect, Navarini et al confirmed that the persistence of elevated C reactive protein (CRP) and high disease activity may be considered predictive factors of CV disease in patients with axial spondyloarthritis.28 These authors also assessed 295 patients with axial spondyloarthritis without personal history of CV disease and observed that persistency of increased CRP levels at each visit and high values of clinical parameters of disease activity help to identify patients with axial spondyloarthritis at higher risk of CV disease.28 We entirely agree with this fact as in our experience disease activity influenced the reclassification of patients with RA into the category of very high CV risk.29

    Baseline data of the PsA included in the CARMA cohort showed that 42.5% of them were receiving biological DMARDs and 74% conventional DMARDs at the time of recruitment. The baseline median CRP of the PsA from the CARMA cohort was 3 mg/L. Moreover, the median Health Assessment Questionnaire, considered as the functionality score, was<0.5 (median 0.4 with an IQR between 0.0 and 0.9). These results indicate that patients with PsA included in the CARMA registry and followed prospectively at rheumatology units underwent tight control of the disease and consequently had low disease activity. Therefore, most of them achieved some of the objectives proposed by Coates et al of minimal disease activity of the disease.30

    In previous studies on carotid ultrasound in patients with PsA without CV risk factors or clinically evident CV disease, we disclosed a high prevalence of macrovascular disease compared with ethnically matched controls.6 Yet, when examining patients with RA, Karpouzas et al disclosed that the utilisation of biological therapy correlated with a decrease in CV disease among these patients.31 Thus, this intervention played a preventive role in CV events by inhibiting coronary plaque formation and stabilising high-risk coronary lesions.31 Additionally, Karpouzas et al showed that the administration of statins altered the influence of inflammation on the development of new coronary plaque. This intervention also predicted both regression and calcification of non-calcified lesions in patients with RA, leading to a long-term reduction in CV risk for these individuals.32 In this context, as highlighted in our previous discussion on mortality after a 5-year follow-up,33 individuals within the CARMA cohort were closely followed at tertiary referral centres with special attention to CV risk factors. Consequently, the tight control of the disease, incorporating frequent use of biological therapy and providing control of traditional CV risk factors, may have contributed to a diminished atherogenic burden among our patients with PsA. This could potentially render them somewhat comparable in terms of CV risk to individuals in the general population.

    In our study, the follow-up of CARMA cohort was 7.5 years, which is a limitation to evaluate CV score estimations as these scores were based on 10 years follow-up period. In this sense, we do not know if the fact that the follow-up was 7.5 years compared with the 10 years as applies to the CV risk equations could have contributed to finding fewer CV events than predicted. Nevertheless, our study has several strengths, primarily attributed to its prospective design. Notably, we examined a substantial PsA patient cohort undergoing regular follow-up visits, ensuring standardised data collection and minimising the likelihood of data loss. Moreover, we wish to highlight the potential significance of using two scores in clinical practice. From our perspective, this approach may be clinically relevant for the patients and it is also feasible. Combining two risk chart algorithms may assist in more effectively identifying PsA patients at a high risk of CV events. This information could empower us to proactively implement measures for preventing future CV events in these individuals.

    In conclusion, in PsA patients under close observation in rheumatology units included in the prospective CARMA project, risk chart algorithms prove highly valuable for distinguishing individuals at low and high CV risk. The integration of QRISK3 and SCORE2 in a comprehensive model demonstrated an optimal combination, leveraging QRISK3’s discrimination capability and SCORE2’s calibration accuracy.

    Data availability statement

    No data are available.

    Ethics statements

    Patient consent for publication

    Ethics approval

    The study, conducted in accordance with the principles outlined in the Declaration of Helsinki, prioritised ethical considerations. Full written informed consent was obtained from all subjects prior to their integration into the project. Approval of the study was granted by the Clinical Research Ethics Committee of Lugo, Galicia (Spain) according to protocol no. 2009/077, and in parallel, approval was requested and obtained from the ethics committee of each participating hospital.

    Acknowledgments

    This manuscript is dedicated to Dr. José L. Fernández Sueiro, passionate rheumatologist of PsA, who took part in the initial design of this project and passed away in 2012. The authors thank the Research Unit from the Spanish Foundation of Rheumatology and Spanish Society of Rheumatology (Madrid, Spain) for helping with statistics and data management. The authors thank also all the health professionals and patients who generously participated in this study.

    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

    • SC, JL and MAG-G are joint senior authors.

    • X @ivanferrazamaro

    • Collaborators CARMA Project Collaborative Group: The members of the CARMA Project Collaborative Group include: José L Fernández Sueiro, José A Pinto-Tasende and Eugenia González de Rábago (Complejo Hospitalario A Coruña, Xubias de Arriba, A Coruña); María J González Fernández, Ramón Huguet Codina, Beatriz Yoldi and Mercedes Ramentol (Instituto Dexeus, Barcelona); Sara Marsal, Gabriela Ávila and Cayetano Alegre (Hospital Universitari Vall d’Hebron, Barcelona); Martina Steiner, Tatiana Cobo and Santiago Muñoz (Hospital Infanta Sofía, Madrid); Fernando Gamero, José García Torón and María P Moreno Gil (Hospital S Pedro de Alcántara, Cáceres); Antonio J Mas, Pilar Espiño, Inmaculada Ros and Mónica Ibáñez (Hospital Son Llatzer, Palma de Mallorca); Jesús Tornero and José A Piqueras (Hospital Univ. de Guadalajara); Raimon Sanmartí and Horacio Berman (Hospital Clinic i Provincial, Barcelona); Oscar Fontseré Patón, Benjamín Fernández Gutiérrez and Lydia Abasolo (Hospital Clínico Univ. San Carlos, Madrid); José M Pina Salvador and María D Fábregas (Hospital de Barbastro, Huesca); Montserrat Romera and Joan M Nolla (Hospital Univ. de Bellvitge, Barcelona); Miren Uriarte-Ecenarro, Jesús A García Vadillo and Rosario García de Vicuña (Hospital Univ. de La Princesa, Madrid); Eugenio Chamizo Carmona (Hospital de Mérida, Badajoz); Antonio Fernández. Nebro, Inmaculada Ureña, María A Belmonte and María V Irigoyen (Hospital Universitario Regional, Málaga); Olga Martínez González (Hospital General Virgen de la Concha); Rebeca Belmonte Gómez, Pastora Granados Bautista, Azucena Hernández Sanz and José Santos Rey (Hospital Virgen de la Salud, Toledo); Carmen O Sánchez González (Hospital del Sureste, Madrid); Javier Bachiller and Antonio Zea (Hospital Ramón y Cajal, Madrid); Francisco J Manero, Chesús Beltrán Audera, Marta Medrano and Ángela Pecondón (Hospital Univ. Miguel Servet, Zaragoza); Celia Erausquin, Soledad Ojeda and Carlos Rodríguez-Lozano (Hospital Dr. Negrín, Las Palmas de Gran Canaria); Jesús Babío Herráez (Hospital de Cabueñes, Gijón); Francisco J López Longo and Indalecio Monteagudo (Hospital Gregorio Marañón, Madrid); Javier del Pino and Ruth López González (Hospital Univ. de Salamanca); Miguel A González-Gay, Alfonso Corrales and María Enriqueta Peiró (Hospital Univ. Marqués de Valdecilla, Santander); José M Senabre and José C Rosas (Hospital de la Marina Baixa, Alicante); Isabel Rotés, Estefanía Moreno and Alba Erra (Hospital de San Rafael, Barcelona); Javier Calvo and Amalia Rueda (Hospital General. Universitario, Valencia); Ingrid Möller and Isabel Rodríguez (Instituto Poal, Barcelona); Carmen Barbadillo (Hospital Universitario Puerta de Hierro, Madrid); Enrique Raya, Pilar Morales and Ana Nieto (Hospital Clínico Univ. San Cecilio, Granada); Ana Ruibal Escribano (Hospital Santiago Apóstol, Vitoria-Gasteiz); Carmen García Gómez (Consorci Sanitari de Terrassa, Terrassa); Sergio Ros Expósito (Hospital de Viladecans, Barcelona); Ginés Sánchez Nievas, Enrique Júdez Navarro and Manuela Sianes Fernández (Hospital General de Albacete); Silvia Martínez Pardo and Manel Pujol (Hospital Mutua Terrassa, Terrassa); Beatriz Gónzález Alvarez and Alberto Cantabrana (Hospital Ntra. Sra. de Candelaria, Santa Cruz de Tenerife); Sagrario Bustabad and Esmeralda Delgado (Hospital Univ. de Canarias, La Laguna, Tenerife); Alejandro Muñoz and Sergio Rodríguez Montero (Hospital Univ. de Valme, Sevilla); Javier Rivera Redondo and Teresa González Hernández (Instituto Provincial de Rehabilitación, Madrid); Francisco J González. Polo (Hospital de Cantoblanco, Madrid); Raúl Menor Almagro (Hospital de Jerez de la Frontera, Jerez de la Frontera, Cádiz); José M Moreno and Emilio Giner Serret (Hospital Obispo Polanco, Teruel); Laura Cebrián Méndez and María Teresa Navío (Hospital Infanta Leonor, Madrid); Teresa Pedraz Penalva and Cristina Fernández Carballido (Hospital General de Elda, Alicante); Encarnación Pagán and Pablo Mesa del Castillo (Hospital Los Arcos, Murcia); Esperanza Naredo and Ana Cruz (Hospital Severo Ochoa, Madrid); Ana Turrión and Ana Pérez Gómez (Hospital Príncipe de Asturias, Madrid); Julio Sánchez and María Galindo (Hospital Univ. 12 de Octubre, Madrid); Eduardo Collantes, Desireé Ruíz and Pilar Font (Hospital Univ. Reina Sofía, Córdoba); Gema Bonilla (Hospital Univ. La Paz, Madrid); Antonio López Meseguer (Hospital Gutiérrez Ortega, Valdepeñas, Ciudad Real); Manuel J Moreno and Luis F Linares (Hospital Virgen de la Arrixaca, Murcia); Mercedes Morcillo and María L González Gómez (Hospital del Escorial, Madrid); José M Aramburu, Natalia A Rivera and Olaia Fernández. Berrizbeitia (Hospital de Basurto, Bilbao); Hye Sang Park, Manel Riera and Yolanda María León (Hospital Dos de Maig, Barcelona); Carolina Pérez-García, Joan Maymó and Miriam Amirall (Hospital del Mar, Barcelona); Jordi Fiter, Julia Fernández Melón and Luis Espadaler (Hospital Universitario Son Espases, Palma de Mallorca); Joaquín Belzunegui and Inmaculada Bañegil (Hospital de Donostia, Donostia); César Díaz (Hospital de la Santa Creu i Sant Pau, Barcelona); Ramón Valls (Hospital de Palamós, Gerona); Iván Castellví and María Bonet (Hospital Comarcal de L’Alt Penedés, Vilafranca del Penedès, Barcelona); Jaime Calvo Alen (Hospital Sierrallana, Torrelavega); Trinidad Pérez Sandoval (Complejo Asistencial de León); Eva Revuelta Evrard (Hospital General de Ciudad Real); Javier R Godo (Hospital General de Móstoles, Madrid); José Antonio González Fernández (Hospital General Universitario de Elche, Alicante); José A Miranda-Filloy (Hospital Xeral Calde, Lugo).

    • Contributors Study conception and design: JPyLB, JL and MAG-G. Data collection: all authors including members of the CARMA Project Collaborative Group. Statistical analysis: FS-A, ZP and JL. Manuscript guarantor: MAG-G and Spanish Society of Rheumatology. All authors were involved in interpretation of data and drafting the article or revising it for intellectual content. All authors approved the final version to be published.

    • Funding This work is promoted by the Spanish Society of Rheumatology. MAG-G research is supported by the Spanish Ministry of Health, Instituto de Salud Carlos III (ISCIII), Spanish Red de Investigación RICORS (RD21/0002/0025).

    • Competing interests None declared.

    • 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.