1. Introduction
Psoriasis and psoriatic arthritis (PsA) are complex entities with a high clinical heterogeneity in which skin-nail manifestations overlap, to varying degrees, with musculoskeletal manifestations. Thus, it is not uncommon for patients to present over time different combinations of peripheral and axial arthritis, dactylitis and/or enthesitis, in addition to cutaneous manifestations. Along with this, other manifestations may occur such as uveitis or inflammatory bowel disease, as well as a wide variety of comorbidities, mainly cardiometabolic [
1]. This clinical complexity has supported the modern concept of psoriatic disease as a systemic entity that goes beyond skin and joint manifestations to encompass all the aforementioned aspects, ultimately conditioning the prognosis, sometimes gloomy, and in many cases the therapeutic approach [
2].
Psoriatic disease has a variegated and complex genetic basis, where the HLA-Cw6 allele stands out, not only as the main genetic determinant of the disease, but also as one of the key responsible for its varied phenotypic expression [
3]. Classically, the presence or absence of this marker has been linked to specific disease endotypes. In addition to being associated with severe, early-onset, familial psoriasis (type I disease), HLA-Cw6 has been associated with longer psoriasis-arthritis latency, a better response to certain therapies (ustekinumab and methotrexate), and very recently, to a better cardiometabolic risk profile [
4,
5,
6]. In this last sense, two independent groups of researchers have related the positivity of this biomarker to a lower risk of hypertension and less visceral adiposity [
4,
5]. Also, very recently, Spanish researchers have found a link between HLA-Cw6 and a lower risk of liver fibrosis in PsA [
7]. In summary, it seems that carriers of this allele would have a better cardiometabolic risk profile, which has sparked renewed interest in the usefulness of routinely requesting this biomarker in the management of psoriatic disease [
8].
However, studies that point to a cardiometabolic benefit of this marker are still few, and the unequal distribution of this allele in different psoriasis populations worldwide demands more studies that expand and support the HLA-Cw6-linked endotype. On the other hand, since HLA-Cw6 is associated with early-onset skin disease, it would be interesting to investigate whether other recently described genetic biomarkers of early-onset disease have the potential to be linked to the cardiometabolic risk so characteristic of the disease [
9]. In the present report, the distribution and possible connection between the HLA-Cw6 allele and IFIH1/MDA5 genetic variants with the different cardiometabolic risk factors is analyzed in a large cohort of patients with psoriatic disease.
2. Materials and Methods
2.1. Study population
Two-center, observational and cross-sectional study of 572 patients with psoriatic disease, of which 30% presented CASPAR criteria for PsA. All of them were over 18 years of age and of Caucasian ancestry, residing in the Principality of Asturias (northwestern Spain, total population 1,000,000 inhabitants). Study cohort was recruited between January 2007 and August 2017 from the dermatology and rheumatology departments of two hospitals. All patients were informed of the purposes of the study and gave their informed consent to participate in this research. The anonymity of the study subjects has been preserved at all times, and the rules of good clinical practice of the Declaration of Helsinki have been respected throughout the study. The study protocol was approved by the research ethics committee of the Principality of Asturias. This patient´s cohort was registered as a Biobank collection by the Spanish Instituto de Salud Carlos III (reference C.0003441).
2.2. Sociodemographic and clinical variables
From the sociodemographic point of view, age, sex, and educational level were collected. Age at disease onset, family history of disease, type of psoriasis and its severity, PASI, nail involvement and treatments received were also recorded. Severe disease was regarded as PASI ≥10. In the case of PsA, and for the purposes of this study, we just recorded its presence according to CASPAR criteria with subsequent diagnostic confirmation by a rheumatologist experienced in this pathology. Smoking habits and alcohol consumption were recorded. Throughout the study period, the presence/absence of high blood pressure, diabetes and its types, dyslipidemia, fatty liver disease, weight, height, BMI, and waist circumference were determined. Basal blood glucose level, total cholesterol, HDL, LDL, and triglycerides were recorded. For the purposes of the study, early disease was regarded as disease onset ≤ 40 years (n: 440), while late onset was reserved for patients with disease onset > 40 years (n: 132).
2.3. Disease risk genotyping
All patients were genotyped for the HLA-Cw6*0602 allele (SNP rs1050414 C/G). The DNA was obtained from whole blood leukocytes, and all the individuals were genotyped for the next IFIH1 SNPs: rs1990760 (g.162267541C>T, c.2719G>A, p.Ala907Thr), rs35337543 (g.162279995C>G, c.1641+1G>C), and rs35744605 (g.162277580C>A, c.1879G>T, p.Glu627Ter). The three were determined by real time PCR Taqman assays (Fisher Scientific, Waltham, MA, USA): rs35337543, assay id C_25985625_10; rs35744605, assay id C_25982959_10; rs1990760, C_2780299_30. The quality of the genotyping method was confirmed by sequencing PCR fragments with different genotypes. We determined the population frequency of the IFIH1 variants in a total of 200 controls aged 21–75 years (50% aged ≤40 years). They were individuals from our general population without other inclusion or exclusion criteria. They were recruited with the main purpose of determining the allele/genotype frequencies of the IFIH1 SNPs in our population. These controls were genotyped for the three SNPs with the above referenced Taqman assays [
9].
2.4. Statistical analysis
For the univariate descriptive analysis, the mean, standard deviation, median, minimum and maximum value of the quantitative variables, and the absolute and relative frequencies (percentages) of the categorical variables are presented. A descriptive and comparative analysis of the variables included based on sex and the HLA-Cw6 marker is also presented. To estimate the crude effect of the genetic markers on the different outcomes, simple linear regression models were constructed when the dependent variable was a quantitative variable, and simple logistic regression models when the dependent variable was categorical. To estimate the adjusted effect, multiple regression models were created (linear or logistic as appropriate) where covariates such as age, sex, time of evolution, and the presence of arthritis were introduced. We performed a logistic regression to compare the difference in sex, severity (PASI ≥ 10), PsA, and the HLA-Cw6 positivity, stratified by age at onset. The allele and genotype frequencies between groups were compared with the chi2 and Fisher's exact tests. For the genetic comparisons, we considered alleles previously related to increased IFIH1 expression or function as putative risk factors (rs1990760 T, rs35337543 G, and rs35744605 C). The threshold for statistical significance was set at p<0.05. Data were analyzed using R software (4.3.1 “Beagle Scouts”).
3. Results
3.1. Study population characteristics
The study cohort was composed of 309 men and 263 women, mean age of 46.7 years (SD 14.5), with mean disease duration of 19.4 years (SD 14.8). The median PASI was 14 (min: 1, max: 49), and 294 (51.4%) patients had a PASI ≥ 10. Nail dystrophy was diagnosed in 328 patients (57.3%), while in 22 (3.8%) this aspect was not recorded in the clinical history. In 222 (38.8%) patients, nail involvement was not detected in any of the visits recorded. A family history of psoriasis and/or PsA was found in 211 patients (36.9%). Regarding the anthropometric variables, the average weight was 77.9 Kg (SD 16.4) with an average BMI of 27.6 (SD 5.02), with the average waist circumference being 98 cm (min: 61, max: 138). Just over a third (34.3%) of the patients were smokers, while the median alcohol consumption according to standard drink units (SDU) was 0 (min: 0, max: 30). Of the study patients, fatty liver disease was ruled out in 405 (70.8%), confirmed in 129 (22.6%) and no information was available in 38 (6.6%). Twenty percent of the patients had hypertension and another 20% had dyslipidemia. Thirty-three patients had adverse coronary events (5.8%). Twenty-two patients were type I diabetics (3.8%), while 45 (7.9%) were type 2. Two hundred and forty-one patients expressed the HLA-Cw6 allele (42.1%). The distribution of
IFIH1/MDA5 genetic variants of the study population has been published elsewhere [
9]. A total of 440 patients had an onset age of ≤40 years compared to 132 cases at age >40 years. Severe disease (PASI ≥ 10), family history of disease, and the presence of HLA-Cw6 were significantly more frequent in patients with onset age ≤40 years (p < 0.001). Upon multiple logistic regression analyses for the presence of PsA, we found significant associations with psoriasis severity (OR: 2.14, 95%CI: 1.46–3.16), female sex (OR: 1.63, 95%CI: 1.12–2.38), and the
IFIH1/MDA5 rs1990760 TT genotype (OR: 1.62, 95%CI: 1.11–2.37), while the HLA- Cw6 allele resulted protective (OR: 0.65, 95%CI: 0.44–0.95).
Table 1 and
Table 2 represent the distribution of the study variables stratified by sex and HLA-Cw6, respectively.
3.2. Relationship between cardiometabolic risk factors and HLA-Cw6
Regarding waist circumference, in the crude model, subjects carrying HLA-Cw6 had an average of 2.37 cm less (95%CI: -4.82, 0.07), p= 0.06. This difference was lost in the model adjusted for age, sex, duration of disease and arthritis [0.05 (-2.19, 2.29), p=0.965]. Regarding BMI, in the crude model, HLA-Cw6 positives had an average of 0.81 units less (95% CI -1.67, 0.05), p=0.06, but in the model adjusted for the previous variables, this difference was no longer apparent [-0.11 (95%CI: - 0.96, 0.74), p=0.80]. HLA-Cw6 subjects reduced the risk of hypertension by 48% (95%CI: 20-66). However, the effect adjusted for age, sex, disease duration and presence of arthritis was a 17% reduction, without being statistically significant (95%CI: 0.49-1.40, p=0.49). Cw6 patients had a lower risk of diabetes, OR: 0.36, 95%CI: 0.19-0.63, p< 0.001, however, after adjustment for age this association was not maintained, OR: 0.56, 95%CI: 0.28-1.05, p=0.08. There was no relationship, in either the crude or adjusted models, between HLA-Cw6, obesity, dyslipidemia, fatty liver disease, or adverse coronary events.
3.3. Relationship between cardiometabolic factors and IFIH1/MDA5 genetic variants
Patients carrying the IFIH1/MDA5 rs1990760 variant had a disease that began on average almost 6 years earlier than subjects without this biomarker (95%CI: -9.93, -1.87, p=0.004). However, there was no association between IFIH1/MDA5 genetic variants and cardiometabolic risk factors, neither in the crude nor in the adjusted regression models.
4. Discussion
In this study we have investigated the connections between two of the genetic biomarkers most consistently linked to the risk of psoriatic disease (HLA-Cw6 and variants of the IFIH1 gene) and the presence of the cardiometabolic comorbidity so typical of this condition. We confirmed the known relationship between the HLA-Cw6 allele and type I disease (early, familial, severe), while PsA was independently and positively associated with female sex, severe psoriasis, and the IFIH1/MDA5 rs1990760 TT variant, but negatively with HLA-Cw6. In the unadjusted regression models, we detected a trend of HLA-Cw6 positive patients towards a better cardiometabolic profile (lower BMI, less visceral adiposity, lower risk of hypertension and lower risk of diabetes). However, in adjusted multiple regression, these differences were no longer apparent. No variant of the IFIH1 gene was found associated with cardiovascular (CV) comorbidity.
It has been known for decades that psoriatic patients carrying the HLA-Cw6 allelic variant present specific features of the disease (HLA-Cw6-linked endotype). This essentially refers to the presence of a more severe psoriasis, with a greater familial inheritance burden and an earlier age of disease onset (type I psoriasis) [
3]. There are also repeated reports regarding the relationship between this allele and a better response to certain therapies [
3,
10], as well as links with other characteristics of psoriatic disease, such as the negative relationship with nail dystrophy [
3,
5], also found in our study (OR 0.40, p<0.001). Although it is known of a certain connection between HLA-Cw6 and a lower presence of diabetes [
11], it has not been until very recent times (last two years) that studies have begun to appear that connect this allele with a better overall cardiometabolic risk profile [
4,
5,
6]. In fact, at least two recent independent studies carried out in populations of European Caucasian ancestry have found associations between HLA-Cw6 and a smaller waist circumference, less visceral adiposity, and in general, a “protective” tendency against the development of certain CV risk factors such as hypertension. Thus, in these studies, the reduction in the risk of hypertension conferred by HLA-Cw6 has been around 30% [
4,
5]. A study carried out in a Spanish population also found a reducing effect of HLA-Cw6 on the risk of liver fibrosis in subjects with PsA, while curiously enough, the presence of hypertension increased this risk [
7]. A fourth study carried out in a Scandinavian population found less visceral adiposity in HLA-Cw6 psoriasis subjects [
6]. However, although our results presented a similar trend to the aforementioned data in the unadjusted models, we could not confirm these findings in the adjustments based on age, sex, disease duration, and the presence of arthritis. Therefore, our findings seem more in line with those works that independently find a greater prevalence of cardiometabolic factors in subjects with PsA in relation to patients with psoriasis who do not develop arthritis [
12,
13,
14]. This is especially interesting if we take into account that, as we were also able to confirm here, the associations between HLA-Cw6 and arthritis (negative association) seem to go in the opposite direction to the same relationship when analyzed in patients with psoriasis without arthritis (positive association). In other words, the possible relationship between the negativity for this allele and a greater presence of cardiometabolic comorbidity seems to be better explained by age and the presence of PsA.
What then could be the reason for the disparity in the results found between our work and others from recent literature? First of all, we must consider that, for example, in the work by Douroudis et al [
5], more than half of the patients were HLA-Cw6 positive, while in our cohort it was slightly higher than 40%, and in the HUNT4 study by Solvin et al. it not exceeded 30% [
6]. Furthermore, the prevalence of psoriasis in the HUNT4 study was greater than 6% of the population [
6], a figure that clearly exceeds the prevalence of the disease in our latitude [
15]. Another possibility is the different frequency that these comorbid factors (hypertension or diabetes, for example) may have in these studies, as well as the different duration of the disease. Regarding this last point, the average duration of the disease in our cohort did not exceed 20 years, while in the studies by Douroudis and that of Solvin it clearly exceeded this threshold [
5,
6]. As our study is a punctual observation in time, we cannot rule out that with a longer follow-up period there may be a higher prevalence of some of the comorbidities analyzed and that this in turn represents a relevant change in our results.
Age and disease duration turned out to be explanatory factors for some of the associations found here between the disease and cardiometabolic factors. Therefore, it seemed a priori interesting to include other genetic markers that have been related to the age of onset of the disease (variants in the IFIH1/MDA5 gene) and to verify whether these variants contributed to a differential distribution of CV comorbidity factors. However, although the disease began an average of 6 years earlier in subjects carrying the IFIH1/MDA5 rs1990760 variant, neither this nor any other allelic variant of the IFIH1/MDA5 gene explained a different distribution of CV comorbid factors in our population.
Although in the work of Macia-Villa et al. HLA-Cw6 was protective against liver fibrosis (an advanced stage of fatty liver disease) in subjects with PsA [
7], we could not find any relationship between fatty liver disease and the different allelic variants analyzed here. In any case, it is important to highlight that the prevalence of fatty liver in our cohort was well below what would be expected according to the figures that have been previously published in this regard [
16], so our findings at this specific point should be taken with caution.
The female population in our study had a higher frequency of arthritis and HLA-Cw6 positivity, similar to what was found by Douroudis et al [
6]. However, none of the differences between genders in the distribution of CV comorbidity factors (
Table 1) could be explained by the genetic markers included in this study. In any case, it is necessary to highlight that the genetic associations reported with respect to CV risk in psoriatic disease are quite heterogeneous, which makes it difficult to have a clear position on the matter [
17,
18,
19,
20].
Therefore, what explained the presence of cardiometabolic comorbidity regardless of HLA-Cw6 in our study? Starting with waist circumference, a surrogate marker of visceral adiposity, for each decade of age this circumference increased by 3 cm (p<0.001); women had an average of almost 10 cm less waist circumference (p<0.001); evolution time, so that individuals with a short evolution time had a lower waist circumference (p=0.021); arthritis, so that patients who developed PsA had 3.68 cm more waist circumference than patients without arthritis (p=0.003). The BMI analysis showed results in line with the above. Regarding hypertension, for each decade of age, the risk increased by 11% (p<0.001), and arthritis increased this risk by 65% (p<0.05). Diabetes was essentially explained by older age, so that for each decade this risk increased by 8% (p<0.001). For each decade of age, the risk of dyslipidemia increased by 7% (p<0.001), while arthritis increased this risk by 86% (p=0.008). The risk of fatty liver disease increased by 2% for each decade of age (p=0.02), and was reduced by 74% in women compared to men (p<0.001). Adverse coronary events increased by 11% per decade (p<0.001). Interestingly, the finding of a relationship between PsA and the risk of certain factors such as hypertension or dyslipidemia had also been collected by other studies [
12,
13,
14,
21,
22].
In a previous work by our group, we had detected a cardiometabolic endotype linked to HLA-Cw6 [
4]. However, that study was carried out for clinical phenotyping purposes (in fact, the characterization of the PsA population was much finer), while the present cohort was constructed for genotyping studies, so the collection of variables was much more selective and exhaustive (for example, in the first cohort no information on fatty liver disease was collected, while in the cohort of the present study it was). Therefore, both studies are not an extension of each other, but rather very different.
The main limiting factor of our findings is the punctual and cross-sectional nature of the observation, which logically prevents causal inferences. However, we included a substantial number of well-characterized patients with a relatively long mean follow-up. Being a cohort designed for disease risk genotyping purposes, all variables have been carefully checked. The loss of data has also been very small. As we have mentioned previously, we do not know if with a longer follow-up observation period a greater number of comorbidities could have been detected, and, therefore, a change in the results of the present study. We have not included important aspects such as tobacco or alcohol consumption among the corrective factors of the regression models; however, we do not believe there is any connection between the genetic markers analyzed and these harmful habits. On the other hand, arthritis has been introduced as a dichotomous variable without us being able to provide a finer characterization of these patients.
Author Contributions
Conceptualization, Rubén Queiro, Pablo Coto-Segura and Eliecer Coto; Data curation, Estefanía Pardo, Stefanie Burger, Marta Loredo; Formal analysis, Ignacio Braña, Estefanía Pardo, Marta Loredo, Stefanie Burger, Pablo González del Pozo, Eva Fernández-Bretón, Pablo Coto-Segura and Paula Alvarez; Funding acquisition, Pablo Coto-Segura; Investigation, Ignacio Braña, Estefanía Pardo, Marta Loredo, Pablo González del Pozo, Paula Alvarez, Eva Fernández-Bretón, Eliecer Coto and Pablo Coto; Methodology, Rubén Queiro, Marta Loredo, Eva Fernández-Bretón and Eliecer Coto; Resources, Rubén Queiro; Supervision, Eliecer Coto; Validation, Estefanía Pardo, Stefanie Burger, Pablo González del Pozo, Paula Alvarez; Visualization, Ignacio Braña, Pablo González del Pozo, Paula Alvarez, Rubén Queiro; Writing – original draft, Rubén Queiro; Writing – review & editing, Rubén Queiro, Ignacio Braña, Marta Loredo, Stefanie Burger and Eva Fernández-Bretón. All authors have read and agreed to the published version of the manuscript.