Objective. Rheumatoid arthritis (RA) is the most common autoimmune arthropathy worldwide. The increased prevalence of cardiovascular disease (CVD) in RA is not fully explained by classic risk factors. The aim of this study was to determine the influence of rs1058587 SNP within GDF15(MIC1) gene on the risk of CVD in a Colombian RA population. Methods. This was a cross-sectional analytical study in which 310 consecutive Colombian patients with RA and 228 age- and sex-matched controls were included and assessed for variables associated with CVD. The mixed cluster methodology based on multivariate descriptive methods such as principal components analysis and multiple correspondence analyses and regression tree (CART) predictive model were performed. Results. Of the 310 patients, 87.4% were women and CVD was reported in 69.5%. Significant differences concerning GDF15 polymorphism were not observed between patients and controls. Mean arterial pressure, current smoking, and some clusters were significantly associated with CVD. Conclusion. GDF15 (rs1058587) does not influence the development of CVD in the population studied.
1. Introduction
Rheumatoid arthritis (RA) is the most common autoimmune arthropathy worldwide. It is a chronic, multifactorial, and systemic disease characterized by the presence of long-standing inflammation [1, 2]. As with most autoimmune diseases, RA predominantly affects women [2, 3], and its prevalence had been reported to be between 0.3 and 1.6% [4]. In addition to diarthrodial joints, RA can damage virtually any organ thus leading to potential extra-articular manifestations (EAMs), including cardiovascular disease (CVD) [5, 6]. CVD is the major predictor of poor prognosis and represents the main cause of death in this population [7, 8]. It accounts for 30–50% of all deaths in RA patients [9] and the prevalence of CVD in RA Latin American population has been estimated to be around 35% [8].
In the etiological and pathogenic mechanisms of CVD in RA the inflammatory activity plays an important role. Several studies assessing traditional risk factors for CVD in RA have been published [10]. Nevertheless, the increased prevalence of cardiovascular events in RA is not fully explained by these classic risk factors. Nontraditional risk factors have also been identified and categorized into three groups: genetic, RA-related, and others [7]. Recently, an association between a nonsynonymous single nucleotide polymorphism (SNP) at growth differentiation factor 15 (GDF15) gene (i.e., rs1058587) and CVD in RA patients was reported [11].
GDF15, also known as macrophage inhibitory cytokine-1 (MIC1), depending on the tissue that was described, is a growth and differentiation factor which belongs to the superfamily of transforming growth factor-beta (TGF-β). GDF15 plays multiple roles in various pathologies such as CVD, obesity, inflammation, and cancer due to its anti-inflammatory, antiproliferative, and antitumoral properties [11–13].
Its mechanism of action is not fully elucidated but it is believed that MIC1 activates the TGF-β receptors. Specifically in inflammation, it can limit the later phases of macrophage activation, inhibit the production of tumoral necrosis factor-α (TNF-α) induced by lipopolysaccharides, and also regulate the role of the proinflammatory cytokine interleukin-6 (IL-6) [10]. Additionally, MIC1 influences metabolism of carbohydrates and lipids. For instance, patients with obesity and type 2 diabetes mellitus (T2DM) have the highest serum concentrations of this factor [14].
Previous studies have found that some GDF15(MIC1) polymorphisms are associated with different atherothrombotic manifestations such as stroke, transient ischemic attack, deep vein thrombosis, and pulmonary embolism. GDF15(MIC1) had also been associated with RA per se, proposing that it acts as a proinflammatory cytokine as well as a common etiologic factor for RA and CVD [11, 12]. The aim of this study was to determine the influence of rs1058587 SNP within GDF15(MIC1) gene on the risk of CVD in a Colombian RA population.
2. Methods2.1. Study Population
This was a cross-sectional analytical study in which 310 consecutive Colombian patients with RA and 228 age- and sex-matched controls were included. The sample size was not calculated. This was a nonprobability sample. The subjects were seen at the Center for Autoimmune Diseases Research (CREA) in Medellín and Bogotá, Colombia. Patients fulfilled the 1987 American College of Rheumatology (ACR) classification criteria for RA [15]. This study was undertaken between February 2008 and April 2010 and done in compliance with Act 008430/1993 by Ministry of Health of the Republic of Colombia. The institutional review board of the Universidad del Rosario approved the study design.
Each patient was evaluated by a rheumatologist. The information on patient sociodemographic and cumulative clinical and laboratory data was obtained by interview, physical examination, and chart review. Household description by questionnaire and clinical evaluation of the affected family members was done with the same methodology as above. All data were collected in an electronic and secure database.
Sociodemographic variables included age at RA onset, disease duration, educational status, socioeconomic status (SES), current occupation, smoking habits, and physical activity. Definitions of these variables are as follows. Age at onset (AOD) is the age at which patients began to suffer from pain, typical morning stiffness (more than 1 hour), and symmetrical inflammation of hand and/or foot joints. Disease duration is the difference between AOD and the date of first participation in the study. Educational level was recorded as the number of years of education and was also divided into two groups—more or less than 9 years of education based on the General Law of Education in Colombia laws [16, 17]. SES was categorized on the basis of national legislation and was divided into high status, medium status, and low status. Information on current occupational status was also requested. Familial autoimmunity (FA) was defined as the presence of any AD in first degree relatives (FDRs) of the proband [18].
Erosions were defined as having at least one unequivocal cortical bone defect evaluated by a rheumatologist [19]. EAMs were defined as the presence of at least one of the following: skin and lung nodules, skin ulcerations, episcleritis, vasculitis, neuropathy, pleural effusion, pulmonary hypertension, or embolism. CVD was categorized as positive if any of the following variables were present: hypertension, coronary artery disease (CAD), occlusive arterial disease, carotid disease, or thrombosis [20]. Regarding medical treatment, current or past use of methotrexate and other conventional disease modifying antirheumatic drugs (DMARDs) such as sulfasalazine, D-penicillamine, gold salts, and leflunomide; biological therapy (etanercept, infliximab, adalimumab, abatacept, tocilizumab, and rituximab); and azathioprine and cyclosporine were also assessed. In addition, steroid therapy (prednisolone, methylprednisolone, and deflazacort) and antimalarials (chloroquine and hydroxychloroquine) were taken into account.
Autoantibodies including rheumatoid factor (RF), anti-citrullinated protein (ACPA), and anti-nuclear antibodies (ANA) were extracted from the patient’s clinical record.
2.2. Assessment of Traditional Risk Factors for CVD
Patients were assessed for traditional CVD risk factors including current age (≥45 and ≥55 years for men and women, resp.) [21] and history of premature CAD in FDR [22], and all individuals were asked about smoking and physical activity [23]. T2DM was defined as having a fasting plasma glucose level ≥7 mmol/L (126 mg/dL) or taking any antidiabetic agents at the time of assessment [24]. A diagnosis of dyslipidemia was given if the patient had (a) hypercholesterolemia, defined as taking lipid-lowering medication as a surrogate or having a fasting plasma total cholesterol ≥200 mg/dL, (b) HDL ≤40 mg/dL, (c) hypertriglyceridemia (triglycerides ≥150 mg/dL), or (d) elevated LDL (≥100 mg/dL) [25] and hypertension (defined as having a blood pressure ≥140/90 mmHg or using any antihypertensive medication). Systolic and diastolic blood pressures were measured twice with at least a 15-minute interval in between and the averages were recorded [26].
2.3. Anthropometric Measurements
A body mass index (BMI) ≥25 kg/m2 (overweight and obesity) was considered abnormal [27]. Values of waist circumference (≥102 cm for men and ≥88 cm for women) and waist-to-hip ratio (WHR) ≥0.9 for men and ≥0.85 for women were considered indicators of abdominal obesity. Waist circumference was measured around the midpoint between the lowest rib and the iliac crest after exhaling and viewed from the front. Hip circumference was measured at the point of maximum extension of the buttocks when viewed from the side [28]. Abnormal WHR values are consistent with National Cholesterol Education Program Adult Treatment Panel III and World Health Organization definitions [29, 30]. The above measurements were used to find the most reliable predictors of vascular involvement.
2.4. CVD and Assessment of Carotid Intima-Media Thickness (IMT)
CVD was categorized as positive if any of the following subphenotypes were present: hypertension, coronary artery disease (CAD), occlusive arterial disease, carotid disease, or thrombosis [20]. In addition, an expert blind to the patients’ clinical records performed ultrasound measurements of the IMT from the common carotid artery. High-resolution, 2-dimensional images were obtained using an ultrasound machine (Agilent-Hewlett-Packard, Santa Clara, CA). A 7.5 MHz linear-array transducer with an axial resolution of 0.15 mm and a penetration depth from 1.0 to 5.0 cm was used [31]. The patient rested on the examination table for 15 minutes before the initial carotid ultrasound scan. The measurement was applied to the far wall of the right and left carotid arteries. Following R. Salonen and J. T. Salonen [32], a B-mode screening method was used with electronic calipers within 10 mm proximal to the common carotid bifurcation in a temperature-controlled room (22°C to 24°C) [33]. Ultrasound images were recorded on videotape (Sony MD385). IMT was measured at the site of the greatest thickness and at 2 additional points: 1 cm upstream and 1 cm downstream from this site. The average of these 3 values was computed. Since carotid IMT greater than 0.90 mm is included among the definitions of subclinical organ damage [34], we established this point as having severe subclinical atherosclerosis (AT). The reproducibility of the IMT was evaluated in 12 volunteers, by taking 2 measurements one month apart; obtaining an intraclass correlation of 0.98.
2.5. Genotyping
Genomic DNA was extracted from buffy-coat cells using standard methods. All the identified genomic DNA samples were genotyped for GDF15(MIC1) rs1058587 at the Oklahoma Medical Research Foundation using the TaqMan allelic discrimination assay (Applied Biosystems). The GDF15(MIC1) gene, located at band p13.11 on chromosome 19, has two exons that encode the 308-amino acid GDF15(MIC1) polypeptide, consisting of a 29-amino acid signal peptide, a 167-amino acid propeptide, and a 112-amino acid mature protein. Cleavage of the propeptide allows the mature protein to be secreted as a disulfide-linked homodimer. rs1058587 SNP is at codon 202 (CAC to GAC) and results in a histidine to aspartic acid substitution at position 6 of the mature protein (H6D) (http://www.ncbi.nlm.nih.gov/gene/9518).
2.6. Statistical Analysis
Hardy-Weinberg equilibrium was checked in each genotypic marker. Univariate analyses were performed as follows: the categorical variables were analyzed by the frequencies. The quantitative continuous variables were expressed as the mean and standard deviation, as well as the median and range. Bivariate analysis of CVD outcome (subphenotypes and IMT) versus genetic, clinical, and demographic variables was assessed by means of chi square and Kruskal-Wallis test. The mixed cluster methodology proposed by Morineau et al. [35] based on multivariate descriptive methods such as principal components analysis (PCA) and multiple correspondence analysis (MCA) was performed to resume some sets of variables that have strong associations. For example, AOD and duration of disease were analyzed in this setting to derive three groups that resemble the associations of these variables (i.e., time cluster). The same was done for SES and educational level deriving four groups (i.e., sod-cluster) and, for Sjögren’s syndrome (SS), EAM, and comorbidity three groups were determined (i.e., clinical cluster). This allows the integration of these sets of variables that are confounded, via the new clusters, in association models like logistic regression and classification and regression trees (CART).
CART were used to find predictive factors for CVD in RA. The CART model was adjusted in the Salford Predictive Modeler software v7 using different splitting rules and cross-validation samples to accurately estimate the relative errors of the classification tree. As the independent factors, the model included the variables that were statistically significant in bivariate analyses and those variables that were biologically plausible, as well as the new clusters derived. Bivariate and multivariate analyses were performed in R 3.0.2 [36].
3. Results
Of the 310 patients, 87.4% were women and CVD was reported in 69.5%. The SNP rs1058587 was in Hardy-Weinberg equilibrium (P>0.05). Allelic frequencies corresponding to rs1058587 SNP C and G were 82.7% and 17.5%, respectively. The most frequent genotype was CC (68.4%). The frequencies of RF and ACPA were 68.3.% and 80.4%, respectively. Characteristics of the cohort are illustrated in Table 1. Significant differences concerning GDF15 polymorphism (rs1058587) were not observed between patients and controls (Table 2).
Demographic and clinical characteristics of 310 patients with RA.
Variable
Median (range)
Age (y)
37 (20–54)
Age at onset (y)
50 (26–67)
Duration of the disease (y)
12 (1–20)
BMI
29 (19–37)
EAM duration (m; media ± SD)
11.7 ± 3.3
Genetics
n/N (%)
rs1058587 CC
212/310 (68.4)
rs1058587 CG
88/310 (28.4)
rs1058587 GG
10/310 (3.22)
Allele C
513/620 (82.7)
Allele G
109/620 (17.5)
Sociodemographic
n/N (%)
Female
271/310 (87.4)
Civil status: married
174/282 (61.7)
Occupation: household duties
109/287 (38)
Socioeconomic status: medium
195/278 (70.1)
High educational level
169/212 (79.7)
Current smoking
49/272 (18)
RA-related
n/N (%)
Polyautoimmunity: SS
16/310 (5.1)
Familial autoimmunity
47/272 (17.3)
EAM&
117/303 (38.6)
Comorbidity#
225/305 (73.7)
Rheumatoid factor (+)
162/237 (68.3)
ACPA (+)
107/133 (80.4)
ANA (+)
31/75 (41.3)
Erosive (+)
99/110 (90)
Methotrexate
248/273 (90.8)
DMARDs, ever¥
83/117 (71)
Antimalarials, ever£
158/256 (61.7)
Steroids, everæ
250/269 (92.9)
Biological therapy, ever€
37/256 (14.4)
Cardiovascular
n/N (%)
Cardiovascular disease₳
130/187 (69.5)
cIMT abnormal
66/118 (55.9)
Familial cardiovascular disease
19/245 (7.7)
BMI: obesity
18/118 (15.2)
BMI: overweight
42/118 (35.6)
Physical activity
20/121 (16.5)
Abnormal waist-hip ratio
90/117 (76.9)
ACPA: anti-citrullinated protein antibodies; ANA: anti-nuclear antibodies; BMI: body mass index; cIMT: carotid intima-media thickness; DMARDs: disease modifying antirheumatic drugs; EAMs: extra-articular manifestations; HLA: human leukocyte antigen; m: months; MIC1: macrophage inhibitory cytokine-1; RA: rheumatoid arthritis; SD: standard deviation; SS: Sjögren’s syndrome; y: years.
&EAM is defined as the presence of at least one of the following: nodules, skin ulcerations, episcleritis, scleritis, vasculitis, neuropathy, pleural effusion, pulmonary nodules, or pulmonary hypertension.
#It is defined as the presence of at least one of the following: type 2 diabetes mellitus, dyslipidemia, kidney disease, anemia, osteoporosis, depression, fibromyalgia, acid peptic disease, epilepsy, or infections (e.g., hepatitises A, B, and C, tuberculosis, or malaria).
¥If patient had used at least one of the following: sulfasalazine, D-penicillamine, gold salts, and leflunomide.
£If patient had used at least one of the following: chloroquine and hydroxychloroquine.
æIf patient had used at least one of the following: prednisolone, methylprednisolone, and deflazacort.
€If patient had used at least one of the following: etanercept, infliximab, adalimumab, abatacept, tocilizumab, and rituximab.
₳It was categorized as positive if any of the following variables were present: hypertension, coronary artery disease, occlusive arterial disease, carotid disease, or thrombosis.
GDF15 polymorphism (rs1058587) in Colombian patients with RA.
Genotype
RAN=310 (%)
CTR*N=228 (%)
CC
212 (68.4)
151 (66.2)
CG
88 (28.4)
70 (30.7)
GG
10 (3.2)
7 (3.1)
Allele
2N=620
2N=456
C
512 (82.6)
372 (81.6)
G
108 (17.4)
84 (18.4)
∗Significant differences were not observed.
Data analysed by PCA and MCA were resumed in clusters under clinical, time, and sociodemographic categories. The clinical clusters was gathered into 3 groups: (1) the highest frequency comorbidities (e.g., T2DM, dyslipidemia, osteoporosis, acid peptic disease, kidney disease, depression, periodontal disease, and fibromyalgia), not having polyautoimmunity with SS and moderate frequency of EAM; (2) presence of comorbidities and polyautoimmunity with SS and the highest frequency of EAM; (3) absence of comorbidities, not having polyautoimmunity with SS and lower frequency of EAM. The time clusters was categorized into 3 groups: (1) later AOD and shorter duration of the disease; (2) earlier AOD and longer duration of the disease; (3) earlier AOD and shorter duration of the disease. Finally, the sod-clusters was classified into 4 groups: (1) medium SES and the highest educational level; (2) medium SES and lower educational level; (3) the highest SES and higher educational level; and (4) the lowest SES and higher educational level.
The patients with CVD were significantly older and had a later AOD of RA. Duration of EAM, mean arterial pressure, current smoking, clinical cluster number 1, and time cluster number 1 were all risk factors significantly associated with CVD. Associations between GDF15 rs1058587 SNP and CVD (subphenotypes) in RA patients (Table 3) or with IMT were not found.
Characteristics associated with CVD in RA patients.
&Clinical cluster 1 corresponds to patients with the highest frequency of comorbidity, with moderate frequency of EAM, and without Sjögren’s syndrome.
#Time cluster 1 corresponds to patients with later age at onset and shorter duration of the disease.
In Tables 4 and 5 are mentioned the variables associated with the presence of allele C and allele G of the MIC1 rs1058587 SNP, respectively. Once again, we did not find associations between this SNP and CVD in those patients.
Characteristics associated with allele C of the GDF15 rs1058587 SNP in RA patients.
&Clinical cluster 1 corresponds to patients with the highest frequency of comorbidity, with moderate frequency of EAM, and without Sjögren’s syndrome.
#Time cluster 1 corresponds to patients with later age at onset and shorter duration of the disease.
In the genotype analysis, the homozygotes G/G were significantly older (median 54, IQR 11.7 years), had later AOD (median 43.7, IQR 9.92 years), and had the highest diastolic blood pressure (median 88.33, IQR 8.66 mmHg). On the other hand, the genotype CC was found as protective factor to have positive ANA (OR 0.27; 95% C.I. 0.11–0.85, P=0.019).
3.1. CART Predictive Model
The CART model displayed a high proportion of patients with RA through clinical clusters 1 and 2 who also develop CVD (nodes 1 and 2). However, there was association between CVD and sod-cluster, mean arterial pressure (MAP), and current smoking; in fact the major association with CVD is in sod-clusters 2 and 3. However, sod-clusters 1 and 4 show a better discrimination for CVD (node 3). Moreover, in the same group of patients, CVD was higher in those with high MAP (i.e., more than 100 mmHg) (node 4). Finally, a greater proportion of CVD was observed in patients in node 4 with low MAP but presenting with current tobacco exposure (node 5) (Figure 1). This model had an AUC of 0.73 (cross-validated AUC), which represents an adequate predictive performance.
CART predictive model: the main predictive factors significantly associated with the development of CVD in RA are clinical and sociodemographic cluster, mean arterial pressure (MAP), and current smoking, with variable importance scores higher than 0.73.
4. Discussion
Several previous studies of gene polymorphisms in patients with RA who had CVD had demonstrated genetic factors implicated in the development of this comorbidity, such as HLA-DRB1 shared epitope alleles [37, 38]. Moreover, the choice of genes for analysis in those studies was based on a relationship with RA per se or with inflammation. Herein, we were unable to replicate the influence of GDF15(MIC1) rs1058587 SNP in Colombian patients with RA and CVD as was previously observed in Swedish and Australian patients [11, 14]. Differences in admixture patterns and allele frequencies among populations may account for the observed differences in the influence of genetic factors on disease. Our population is highly admixed with Amerindian, African, and Caucasian ancestries as compared with Australians and Swedish who are mainly Caucasians. The rs1058587 G allele frequency was significantly lower in Colombians than in Australians and Swedish [11, 14].
The general prevalence of CVD in our cohort was 69.5%. Some associations between RA patients with CVD were found. Variables that were significantly associated include smoking, higher MAP, longer duration of disease with EAM, clinical cluster number 1, and time cluster number 1. It is important to highlight that these associations encompass only clinical variables that include modifiable traditional risk factors for CVD like smoking and nontraditional risk factors including duration of disease and EAM [5].
GDF15(MIC1) rs1058587 C allele was nevertheless associated with minor frequency of smoking and obesity, lower BMI, and lower measurements of diastolic blood pressure (Table 4), whereas G allele was associated with earlier AOD, clinical cluster number 1, and positive ANA and it was inversely associated with SS and time cluster number 1. Although modest, these findings contrast with reports showing a specific interaction between the G allele and smoking, with an increase of risk for stroke [11].
At the genotype level a statistically significant association was found between G/G genotype and high diastolic blood pressure, while C/C genotype was inversely associated with positive ANA, suggesting a protective role for polyautoimmunity. We were not able to find other studies reporting associations between the rs1058587 homozygote forms and CVD or RA per se. As it is shown, these associations do not confirm a direct relationship between MIC1 and CVD.
Many reasons exist why no significant association was elicited within the cohort and the gen evaluated with the presence of CVD. One of them is the sample size calculation since it was a nonprobability sample and all RA patients from this region were not included in the study. Besides, the heterogeneity of the cohort may prevent us from finding a direct association. A group of patients studied constitute a closed subpopulation into the admixed Colombian community (i.e., Medellín) and another group represents a major variation in ancestry. Otherwise, previous publications addressing the influence of GDF15(MIC1) in CVD were made in general population, not in RA patients. Therefore, the interaction among autoimmunity, admixture population, and genetic, epigenetic, and environmental factors could influence the behavior of GDF15(MIC1) in our study.
5. Limitations of the Study
The aim of this study was to determine the influence of rs1058587 SNP within GDF15(MIC1) gene on the risk of CVD in a Colombian RA population and we are aware of our study limitations. First of all, selection bias could be present in our analysis as not all patients from this region with RA were systematically included. Secondly, although the study sample size is not negligible, it would have been more valuable to have had an appropriate follow-up to establish valid associations between CVD, GDF15(MIC1), and RA. This, in turn, could have improved both internal and external validities. Finally, the cross-sectional nature of the study does not allow us to infer causality. Although we were able to evaluate some associations between GDF15(MIC1) and risk factors for CVD, it would have been more valuable to have a follow-up, as it is possible in cohort multicenter studies, to validate and establish more associations between this gene, CVD, and RA.
6. Conclusions
Contrary to prior reports, we could not find an association between the SNP rs1058587 of the gene GDF15 and the development of CVD (measured through subphenotypes or IMT) in the population studied. There was an association between this SNP and traditional and environmental factors for CVD (i.e., age, obesity, and smoking) and some nontraditional risk factors (i.e., ANA and SS). However, we had not found direct association between the polymorphism and CVD. In fact, MAP, EAM, smoking, duration of the disease, and AOD are highlighting the pivotal role of these factors in our population. CVD in RA patients should be approached following the 5P rule: (1) predicting patients at risk of developing CVD; (2) preventing CVD in susceptible patients; (3) personalizing the approach of CVD based on the premise that susceptibility and severity of the disease are unique in each patient; (4) allowing the patient to participate in the decisions made in order to prevent and treat CVD; and (5) incorporating all data and policies according to the population, since the characteristics and natural history of diseases are population-specific [39].
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors are grateful to the members of the Center for Autoimmune Diseases Research (CREA) for fruitful discussions and to Parvathi Viswanathan for technical support. This work was supported by the Universidad del Rosario, Bogotá, Colombia.
SanmartíR.Ruiz-EsquideV.HernándezM. V.Rheumatoid arthritis: a clinical overview of new diagnostic and treatment approaches201313669870410.2174/156802661131399900922-s2.0-84878773082ScottD. L.WolfeF.HuizingaT. W. J.Rheumatoid arthritis201037697461094110810.1016/S0140-6736(10)60826-42-s2.0-77957196913Barragán-MartínezC.Amaya-AmayaJ.Pineda-TamayoR.MantillaR. D.Castellanos-De la HozJ.Bernal-MacíasS.Rojas-VillarragaA.AnayaJ.-M.Gender differences in latin-american patients with rheumatoid arthritis201296490.e5510.e510.1016/j.genm.2012.10.0052-s2.0-84870488556ChopraA.The WHO ILAR COPCORD Latin America: consistent with the world and setting a new perspective201218416716910.1097/rhu.0b013e31825d929b2-s2.0-84862206059Ortega-HernandezO.-D.Pineda-TamayoR.PardoA. L.Rojas-VillarragaA.AnayaJ.-M.Cardiovascular disease is associated with extra-articular manifestations in patients with rheumatoid arthritis200928776777510.1007/s10067-009-1145-82-s2.0-67349175002Rojas-VillarragaA.Ortega-HernandezO.-D.GomezL. F.PardoA. L.López-GuzmánS.Arango-FerreiraC.HincapieM.-E.BetancurJ. F.Pineda-TamayoR.DiazF. J.AnayaJ.-M.Risk factors associated with different stages of atherosclerosis in Colombian patients with rheumatoid arthritis2008382718210.1016/j.semarthrit.2008.01.0192-s2.0-52149121712Amaya-AmayaJ.Sarmiento-MonroyJ. C.MantillaR.-D.Pineda-TamayoR.Rojas-VillarragaA.AnayaJ.-M.Novel risk factors for cardiovascular disease in rheumatoid arthritis2013562-326728610.1007/s12026-013-8398-72-s2.0-84879549613Sarmiento-MonroyJ. C.Amaya-AmayaJ.Espinosa-SernaJ. S.Herrera-DíazC.AnayaJ.-M.Rojas-VillarragaA.Cardiovascular disease in rheumatoid arthritis: a systematic literature review in Latin America201220121737190910.1155/2012/371909SandooA.van ZantenJ. J. C. S. V.MetsiosG. S.CarrollD.KitasG. D.Vascular function and morphology in rheumatoid arthritis: a systematic review20115011212521392-s2.0-8455517798710.1093/rheumatology/ker275GkaliagkousiE.GavriilakiE.DoumasM.PetidisK.AslanidisS.StellaD.Cardiovascular risk in rheumatoid arthritis: pathogenesis, diagnosis, and management201218842243010.1097/rhu.0b013e31827846b12-s2.0-84871767508ÄrlestigL.Rantapää-DahlqvistS.Polymorphisms of the genes encoding CD40 and growth differentiation factor 15 and in the 9p21.3 region in patients with rheumatoid arthritis and cardiovascular disease201239593994510.3899/jrheum.1113362-s2.0-84860711206UnsickerK.SpittauB.KrieglsteinK.The multiple facets of the TGF-β family cytokine growth/differentiation factor-15/macrophage inhibitory cytokine-1201324437338410.1016/j.cytogfr.2013.05.0032-s2.0-84883809336MehtaR. S.SongM.BezawadaN.WuK.Garcia-AlbenizX.MorikawaT.FuchsC. S.OginoS.GiovannucciE. L.ChanA. T.A Prospective study of macrophage inhibitory cytokine-1 (MIC-1/GDF15) and risk of colorectal cancer20141064dju01610.1093/jnci/dju0162-s2.0-84898920630BrownD. A.MooreJ.JohnenH.SmeetsT. J.BauskinA. R.KuffnerT.WeedonH.MillikenS. T.TakP. P.SmithM. D.BreitS. N.Serum macrophage inhibitory cytokine 1 in rheumatoid arthritis: a potential marker of erosive joint destruction20075637537642-s2.0-3394715389110.1002/art.22410ArnettF. C.EdworthyS. M.BlochD. A.McShaneD. J.FriesJ. F.CooperN. S.HealeyL. A.KaplanS. R.LiangM. H.LuthraH. S.MedsgerT. A.Jr.MitchellD. M.NeustadtD. H.PinalsR. S.SchallerJ. G.SharpJ. T.WilderR. L.HunderG. G.The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis198831331532410.1002/art.17803103022-s2.0-0023945481Law 115, General on Education, 8th February 1994Law 30General on Higher Education, p. 4, 1992Cárdenas-RoldánJ.Rojas-VillarragaA.AnayaJ.-M.How do autoimmune diseases cluster in families? A systematic review and meta-analysis2013111, article 7310.1186/1741-7015-11-732-s2.0-84875011344van der HeijdeD.How to read radiographs according to the Sharp/van der Heijde method19992637437452-s2.0-0032981611British Cardiac SocietyBritish Hyperlipidaemia AssociationBritish Hypertension SocietyBritish Diabetic AssociationJoint British recommendations on prevention of coronary heart disease in clinical practice: summary200032070570810.1136/bmj.320.7236.705TalwalkarP. G.SreenivasC. G.GulatiA.BaxiH.Journey in guidelines for lipid management: from adult treatment panel (ATP)-I to ATP-III and what to expect in ATP-IV201317628635VeronesiG.GianfagnaF.GiampaoliS.ChamblessL. E.ManciaG.CesanaG.FerrarioM. M.Improving long-term prediction of first cardiovascular event: the contribution of family history of coronary heart disease and social status201464758010.1016/j.ypmed.2014.04.0072-s2.0-84899704612FranklinB. A.DurstineJ. L.RobertsC. K.BarnardR. J.Impact of diet and exercise on lipid management in the modern era201428340542110.1016/j.beem.2014.01.0052-s2.0-84901270124American Diabetes AssociationStandards of medical care in diabetes—2012201235supplement 1S11S6310.2337/dc12-s011ReinerZ.CatapanoA. L.de BackerG.GrahamI.TaskinenM.-R.WiklundO.AgewallS.AlegriaE.ChapmanM. J.DurringtonP.ErdineS.HalcoxJ.HobbsR.KjekshusJ.FilardiP. P.RiccardiG.StoreyR. F.WoodD.ESC/EAS guidelines for the management of dyslipidaemias: the task force for the management of dyslipidaemias of the European Society of Cardiology (ESC) and the European Atherosclerosis Society (EAS)201132141769181810.1093/eurheartj/ehr158JonesD. W.HallJ. E.Seventh report of the Joint National Committee on Prevention, Detection, Evaluation, and Treatment of High Blood Pressure and evidence from new hypertension trials20044311310.1161/01.hyp.0000110061.06674.ca2-s2.0-0347418917LiaoY.KwonS.ShaughnessyS.WallaceP.HuttoA.JenkinsA. J.KleinR. L.GarveyW. T.Critical evaluation of adult treatment panel III criteria in identifying insulin resistance with dyslipidemia200427497898310.2337/diacare.27.4.9782-s2.0-1842485893KleinS.AllisonD. B.HeymsfieldS. B.KelleyD. E.LeibelR. L.NonasC.KahnR.Waist circumference and cardiometabolic risk: a consensus statement from shaping America's health: Association for weight management and obesity prevention; NAASO, the obesity society; the American society for nutrition; and the American diabetes association20071551061106710.1038/oby.2007.6322-s2.0-34250682510DayC.Metabolic syndrome, or what you will: definitions and epidemiology200741323810.3132/dvdr.2007.003DeedwaniaP. C.MakiK. C.DicklinM. R.StoneN. J.BallantyneC. M.DavidsonM. H.Application of recent definitions of the metabolic syndrome to survey data from the National Cholesterol Education Program Evaluation Project Utilizing Novel E-Technology (NEPTUNE II)20061529530010.1111/j.1559-4564.2006.05902.x2-s2.0-34548065212WendelhagI.WiklundO.WikstrandJ.Arterial wall thickness in familial hypercholesterolemia: ultrasound measurement of intima-media thickness in the common carotid artery1992121707710.1161/01.atv.12.1.702-s2.0-0026499899SalonenR.SalonenJ. T.Determinants of carotid intima-media thickness: a population-based ultrasonography study in Eastern Finnish men1991229322523110.1111/j.1365-2796.1991.tb00336.x2-s2.0-0025980304PignoliP.TremoliE.PoliA.OresteP.PaolettiR.Intimal plus medial thickness of the arterial wall: a direct measurement with ultrasound imaging19867461399140610.1161/01.cir.74.6.13992-s2.0-0023002525ManciaG.de BackerG.DominiczakA.CifkovaR.FagardR.GermanoG.GrassiG.HeagertyA. M.KjeldsenS. E.2007 guidelines for the management of arterial hypertension: the Task Force for the Management of Arterial Hypertension of the European Society of Hypertension (ESH) and of the European Society of Cardiology (ESC)20072814621536MorineauL.LebartA.PironM.1995Paris, FranceDunodR Development Core Team2013Vienna, AustriaR Foundation for Statistical ComputingGonzalez-GayM. A.Garcia-PorruaC.HajeerA. H.Influence of human leukocyte antigen-DRB1 on the susceptibility and severity of rheumatoid arthritis20023163553602-s2.0-003630505110.1053/sarh.2002.32552KarlsonE. W.ChibnikL. B.CuiJ.PlengeR. M.GlassR. J.MaherN. E.ParkerA.RoubenoffR.IzmailovaE.CoblynJ. S.WeinblattM. E.ShadickN. A.Associations between Human leukocyte antigen, PTPN22, CTLA4 genotypes and rheumatoid arthritis phenotypes of autoantibody status, age at diagnosis and erosions in a large cohort study200867335836310.1136/ard.2007.0716622-s2.0-39549114063AuffrayC.CharronD.HoodL.Predictive, preventive, personalized and participatory medicine: back to the future201028, article 572-s2.0-7865032367110.1186/gm178