Open Access

The lp13.3 genomic region -rs599839- is associated with endothelial dysfunction in patients with rheumatoid arthritis

  • Raquel López-Mejías1,
  • Carlos González-Juanatey2,
  • Mercedes García-Bermúdez3,
  • Santos Castañeda4,
  • José A Miranda-Filloy5,
  • Ricardo Blanco1,
  • Javier Llorca6,
  • Javier Martín3 and
  • Miguel A González-Gay1Email author
Contributed equally
Arthritis Research & Therapy201214:R42

DOI: 10.1186/ar3755

Received: 11 October 2011

Accepted: 1 March 2012

Published: 1 March 2012



Rheumatoid arthritis (RA) is an inflammatory disease associated with accelerated atherosclerosis and high risk of cardiovascular (CV) disease. Since genome-wide association studies demonstrated association between rs599839 polymorphism and coronary artery disease, in the present study we assessed the potential association of this polymorphism with endothelial dysfunction, an early step in atherogenesis.


A total of 128 RA patients without history of CV events were genotyped for rs599839 A/G polymorphism. The presence of endothelial dysfunction was assessed by brachial ultrasonography (brachial flow-mediated endothelium-dependent (FMD)).


Patients carrying the allele G exhibited more severe endothelial dysfunction (FMD%: 4.61 ± 3.94%) than those carrying the wild allele A (FMD%: 6.01 ± 5.15%) (P = 0.08). Adjustment for gender, age at the time of study, follow-up time and classic CV risk factors disclosed a significant association between the rs599839 polymorphism and FMD (G vs. A: P = 0.0062).


Our results confirm an association of the rs599839 polymorphism with endothelial dysfunction in RA.


atherosclerosis cardiovascular disease endothelial dysfunction rs599839 rheumatoid arthritis


Rheumatoid arthritis (RA) is a complex polygenic autoimmune inflammatory disease with high risk of cardiovascular (CV) complications [1]. This is a consequence of accelerated atherosclerosis [1]. Besides classic CV risk factors and chronic systemic inflammation, recent studies have emphasized the relevance of several genetic polymorphisms, such as HLA-DRB1 and TNF, in the susceptibility to CV disease in RA [2, 3].

A major issue in the process of accelerated atherosclerosis in RA is the development of endothelial dysfunction, an early step in the development of atherosclerosis. An important step forward might be to identify high-risk RA patients who would benefit from active therapy to prevent clinical disease. Several noninvasive imaging techniques provide the opportunity to study the relationship of surrogate markers to the development of atherosclerosis. Among them, ultrasound techniques based on flow velocity are considered efficient ways to measure subclinical atherosclerosis. Using brachial artery ultrasonography assessment, we and others have disclosed the presence of endothelial dysfunction expressed by abnormal levels of flow-mediated endothelium-dependent vasodilatation (FMD) in patients without clinically evident CV disease who had either long-standing RA [4] or early-onset RA [5]. We also disclosed that the presence of endothelial dysfunction in RA was at least in part genetically determined [4].

Genome-wide association studies (GWAS) aimed to predict coronary artery disease (CAD) revealed several novel putative loci that may increase the risk of CAD in the general population. In this regard, the polymorphism rs599839 (A > G) seems to be associated with CAD [6, 7] and with higher plasma total and low-density lipoprotein (LDL) cholesterol levels [810]. The non-coding rs599839 variant is located near three coding genes on chromosome lp13.3 genomic region [11]. Physically, the closest genes are PSRC1, encoding for "proline/serine-rich coiled coil protein 1", and CELSR2, encoding for the "cadherin EGF LAG seven-pass G-type receptor 2". Also, located at chromosome 1p13.3 is the SORT1 gene, encoding for a cell surface receptor (sortilin) with multiligand capabilities that has been implicated in insulin mediated glucose uptake. Although rs599839 is located in a non-codificant region, there is strong linkage disequilibrium between rs599839 and variants located within the 3' end of the adjacent CELSR2 gene, including variants with potential functional effects [7]. Although it is postulated that this protein is a receptor involved in contact-mediated communication, its specific function has not been fully clarified.

Results obtained by GWAS, need to be validated with replication studies in different cohorts to confirm these findings [12].

Taking all these considerations together we aimed to determine, for the first time, the potential role of rs599839 polymorphism in the development of endothelial dysfunction in a cohort of RA patients without clinically evident CV disease.

Materials and methods

Patients and study protocol

A series of 128 Spanish RA patients recruited from Lugo (NW Spain) with no previous history of CV disease were included in the present study. The study was approved by the ethics committee of the Hospital Xeral-Calde (Lugo) and a subject's written consent was obtained in all the cases. All patients fulfilled the 1987 American College of Rheumatology criteria for the classification of RA [13]. Information on the main characteristics and CV risk factors of the patients enrolled in the study is shown in Table 1.
Table 1

Main characteristics of the RA patients included in the study

Clinical Feature

% (n/N)



Main characteristics


   Age at the time of disease onset (years, mean ± SD)

50 ± 13.5

   Follow-up (years, mean ± SD)

12.9 ± 7.8

   Percentage of women


   Rheumatoid factor positive

76.5 (98/128)

   Anti-CCP antibodies positive

67.2 (86/128)

   Shared epitope positive

70.3 (90/128)

Cardiovascular risk factors



26.5 (34/128)

   Diabetes mellitus

7.8 (10/128)


19.5 (25/128)


4.6 (6/128)

   Smoking habit

10.9 (14/128)

Anti-CCP antibodies, anti-cyclic citrullinated peptide antibodies; RA, rheumatoid arthritis; SD, Standard deviation.


DNA from patients was obtained from peripheral blood using standard methods.

The rs599839 A/G polymorphism was genotyped with TaqMan SNP genotyping assays (C____972962_10) in a 7900 HT real-time polymerase chain reaction (PCR) system, according to the conditions recommended by the manufacturer (Applied Biosystem, Foster City, CA, USA). Negative controls and duplicate samples were included to check the accuracy of genotyping.

Brachial artery reactivity

Endothelial function was determined using high-sensitivity brachial ultrasonography according to the guidelines for the ultrasound assessment of endothelial-dependent FMD% [4, 14]. B-mode scan of the right brachial artery, in a longitudinal section 2 to 12 cm proximal to the antecubital fossa, was performed in supine participants using a vascular software for two-dimensional imaging, color and spectral Doppler, an internal electrocardiogram (EKG) monitor, and a 7.5-MHz phased-array transducer Hewlett-Packard SONOS 5500 system (Hewlett-Packard, Palo Alto, CA, USA). The anterior and posterior intima-media interfaces were used to define the baseline artery diameter, calculated as the average of measurements made during four cardiac cycles at end diastole. Timing of each image frame with respect to the cardiac cycle was determined with simultaneous EKG recordings on the ultrasound system digital monitor. During image acquisition, anatomic landmarks were noted to maintain the same image of the artery throughout the study using a specific stereotactic clamp. The forearm blood pressure cuff was inflated on the ipsilateral wrist to at least 50 mm Hg above resting systolic blood pressure for five minutes, and then was released. FMD% (an increase in brachial artery diameter) was measured 30 to 60 seconds after cuff release. To assess endothelium-independent vasodilatation (NTG%), we used 400 micrograms of sublingual nitroglycerin, which acts directly on vessel smooth muscle to cause vasodilatation. NTG% was measured four minutes after nitroglycerin intake. In all cases a cardiologist (CG-J) analyzed all of the ultrasound data offline. A FMD value < 7% was considered pathologic, indicating the presence of endothelial dysfunction [14]. Intraobserver variability for FMD and NTG was 1.3% and 1.9%, respectively, based on repeat brachial ultrasonography in 32 individuals. Assessment of the endothelial function of RA patients undergoing anti-TNF-α therapy was performed 24 to 48 hours before drug administration.

Statistical analysis

The association between the genotypes of the rs599839 polymorphism and surrogate markers of subclinical atherosclerosis was tested using unpaired t test to compare between two groups, and one-way analysis of variance (ANOVA) to compare among more than two groups. We also tested the association between these parameters and alleles using analysis of covariance (ANCOVA) adjusting for gender, age, duration of the disease at the time of the ultrasonographic study and traditional CV risk factors (hypertension, diabetes mellitus, dyslipidemia, obesity and a smoking habit).

Statistical significance was defined as P < 0.05. All analyses were performed with STATA statistical software 9.1 (Stata Corp., College Station, TX, USA).


Results of the comparison between the different genotypes and alleles of rs599839 polymorphism according to surrogate markers of subclinical atherosclerosis are shown in Table 2. Patients carrying the allele G exhibited more severe endothelial dysfunction (FMD%: 4.61 ± 3.94%) than those carrying the wild allele A (FMD%: 6.01 ± 5.15%). However, the difference was slightly out of the range of significance (P = 0.08). Likewise, values of NTG% showed a similar trend with a marginal decrease of NTG% in patients carrying the allelic variant G compared with those carrying allele A (P = 0.06). Also, RA patients carrying the GG and AG genotypes had lower FMD% values (1.94 ± 3.98% and 5.35 ± 3.70%, respectively) than those homozygous for the AA genotype (6.15 ± 5.42%). However, the difference did not achieve statistical significance. It was also the case when genotypes were assessed according to NTG% results (Table 2).
Table 2

Comparison of FMD and NTG values in RA patients according to rs599839 polymorphism



mean ± SD (n)



mean ± SD (n)







6.01 ± 5.15 (210)


16.95 ± 8.14 (210)



4.61 ± 3.94 (46)


14.52 ± 6.92 (46)





6.15 ± 5.42 (87)


17.34 ± 8.34 (87)



5.35 ± 3.70 (36)


15.04 ± 7.01 (36)



1.94 ± 3.98 (5)


12.66 ± 6.99 (5)



4.93 ± 3.85 (41)


14.74 ± 6.96 (41)


FMD, flow-mediated endothelium-dependent (post-ischemia) vasodilatation; NTG, flow-mediated endothelial independent (post-nitroglycerin) vasodilatation; RA, rheumatoid arthritis; SD, standard deviation

Since gender, age at the time of ultrasonography study, follow-up time and classic CV risk factors may act as potential confounders of the results derived from the ultrasonography assessment; adjustment for these potential confounders was performed. Following this procedure, a comparison of FMD and NTG values in RA patients according to rs599839 alleles in an adjusted ANCOVA model yielded a significant association between the rs599839 A/G polymorphism and FMD and NTG (G versus A: FMD% P = 0.0062 and NTG% P = 0.041, respectively).


GWA studies have become a powerful approach to rapidly identifying genetic variants that influence susceptibility to common complex diseases. Novel putative loci, such as the rs599839 A/G polymorphism (chromosome 1p13.3), that seem to increase the risk to CAD have been described by this method [6, 7]. This rs599839 polymorphism has also been implicated in the presence of higher plasma total and LDL cholesterol levels [810].

Since GWAS have to be validated by replication studies in different cohorts and endothelial dysfunction, an early step in the atherogenesis, has been described in patients with RA, we assessed for the first time the association between the rs599839 A/G polymorphism and the presence of endothelial dysfunction in a series of RA patients without clinically evident CV disease. Interestingly, an adjusted analysis disclosed that the presence of the mutant allele G was associated with the presence of endothelial dysfunction.

Impaired FMD of the brachial artery due to endothelial dysfunction has been associated with both CV risk factors and future CV morbidity and mortality in the general population [15]. In addition, endothelial dysfunction manifested by impaired FMD was observed in both long-standing RA patients [4] and early-onset RA patients [5] without clinically evident CV disease. These observations support a potential role of FMD in establishing the presence of endothelial dysfunction as a subclinical marker of atherosclerotic disease in RA. The results derived from this study suggest a potential implication of the rs599839 polymorphism in the development of endothelial dysfunction in RA. They also support the claim of a genetic influence in the development of the atherosclerotic disease in RA.


Our results showed that rs599839 polymorphism is associated with the presence of endothelial dysfunction in patients with RA.

Conflicting interests

The authors declare that they have no competing interests.


We thank Rodrigo Ochoa, Sofía Vargas, M. Luisa López, M. Jesús Ibañez and Sara Olavarria for their technical assistance. This study was supported by two grants from "Fondo de Investigaciones Sanitarias" PI06-0024 and PI09/007/48 (Spain). This work was partially supported by RETICS Program, RD08/0075 (RIER) from "Instituto de Salud Carlos III" (ISCIII). MGB is a beneficiary of a grant from Fundación Española de Reumatología (FER).




analysis of covariance


analysis of variance


coronary artery disease


cadherin EGF LAG seven-pass G-type receptor 2






flow-mediated endothelium-dependent (post-ischemia) vasodilatation.


genome-wide association studies


human leukocyte antigen


Low-density lipoprotein


flow-mediated endothelial independent (post-nitroglycerin) vasodilatation


polymerase chain reaction


proline/serine-rich coiled coil protein 1


rheumatoid arthritis


standard deviation


single-nucleotide polymorphism




tumor necrosis factor.


Authors’ Affiliations

Department of Rheumatology, Hospital Universitario Marqués de Valdecilla, IFIMAV
Cardiology Division, Hospital Xeral-Calde
Instituto de Parasitología y Biomedicina López-Neyra, IPBLN-CSIC
Rheumatology Department, Hospital Universitario la Princesa, IIS-Princesa
Division of Rheumatology, Hospital Xeral-Calde
Department of Epidemiology and Computational Biology, School of Medicine, University of Cantabria, and CIBER Epidemiología y Salud Pública (CIBERESP), IFIMAV


  1. Chung CP, Oeser A, Raggi P, Gebretsadik T, Shintani AK, Sokka T, Pincus T, Avalos I, Stein CM: Increased coronary-artery atherosclerosis in rheumatoid arthritis: relationship to disease duration and cardiovascular risk factors. Arthritis Rheum. 2005, 52: 3045-3053. 10.1002/art.21288.View ArticleGoogle Scholar
  2. Gonzalez-Gay MA, Gonzalez-Juanatey C, Lopez-Diaz MJ, Pineiro A, Garcia-Porrua C, Miranda-Filloy JA, Ollier WE, Martin J, Llorca J: HLA-DRB1 and persistent chronic inflammation contribute to cardiovascular events and cardiovascular mortality in patients with rheumatoid arthritis. Arthritis Rheum. 2007, 57: 125-132. 10.1002/art.22482.View ArticleGoogle Scholar
  3. Rodriguez-Rodriguez L, Gonzalez-Juanatey C, Palomino-Morales R, Vazquez-Rodriguez TR, Miranda-Filloy JA, Fernandez-Gutierrez B, Llorca J, Martin J, Gonzalez-Gay MA: TNFA -308 (rs1800629) polymorphism is associated with a higher risk of cardiovascular disease in patients with rheumatoid arthritis. Atherosclerosis. 2011, 216: 125-30. 10.1016/j.atherosclerosis.2010.10.052.View ArticleGoogle Scholar
  4. Gonzalez-Juanatey C, Testa A, Garcia-Castelo A, Garcia-Porrua C, Llorca J, Vidan J, Hajeer AH, Ollier WE, Mattey DL, Gonzalez-Gay MA: HLA-DRB1 status affects endothelial function in treated patients with rheumatoid arthritis. Am J Med. 2003, 114: 647-652. 10.1016/S0002-9343(03)00133-5.View ArticleGoogle Scholar
  5. Vaudo G, Marchesi S, Gerli R, Allegrucci R, Giordano A, Siepi D, Pirro M, Shoenfeld Y, Schillaci G, Mannarino E: Endothelial dysfunction in young patients with rheumatoid arthritis and low disease activity. Ann Rheum Dis. 2004, 63: 31-35. 10.1136/ard.2003.007740.PubMed CentralView ArticleGoogle Scholar
  6. Samani NJ, Braund PS, Erdmann J, Gotz A, Tomaszewski M, Linsel-Nitschke P, Hajat C, Mangino M, Hengstenberg C, Stark K, Ziegler A, Caulfield M, Burton PR, Schunkert H, Tobin MD: The novel genetic variant predisposing to coronary artery disease in the region of the PSRC1 and CELSR2 genes on chromosome 1 associates with serum cholesterol. J Mol Med (Berl). 2008, 86: 1233-1241. 10.1007/s00109-008-0387-2.View ArticleGoogle Scholar
  7. Samani NJ, Erdmann J, Hall AS, Hengstenberg C, Mangino M, Mayer B, Dixon RJ, Meitinger T, Braund P, Wichmann HE, Barrett JH, König IR, Stevens SE, Szymczak S, Tregouet DA, Iles MM, Pahlke F, Pollard H, Lieb W, Cambien F, Fischer M, Ouwehand W, Blankenberg S, Balmforth AJ, Baessler A, Ball SG, Strom TM, Braenne I, Gieger C, Deloukas P, et al: Genomewide association analysis of coronary artery disease. N Engl J Med. 2007, 357: 443-453. 10.1056/NEJMoa072366.PubMed CentralView ArticleGoogle Scholar
  8. Kathiresan S, Melander O, Guiducci C, Surti A, Burtt NP, Rieder MJ, Cooper GM, Roos C, Voight BF, Havulinna AS, Wahlstrand B, Hedner T, Corella D, Tai ES, Ordovas JM, Berglund G, Vartiainen E, Jousilahti P, Hedblad B, Taskinen MR, Newton-Cheh C, Salomaa V, Peltonen L, Groop L, Altshuler DM, Orho-Melander M: Six new loci associated with blood low-density lipoprotein cholesterol, high-density lipoprotein cholesterol or triglycerides in humans. Nat Genet. 2008, 40: 189-197. 10.1038/ng.75.PubMed CentralView ArticleGoogle Scholar
  9. Sandhu MS, Waterworth DM, Debenham SL, Wheeler E, Papadakis K, Zhao JH, Song K, Yuan X, Johnson T, Ashford S, Inouye M, Luben R, Sims M, Hadley D, McArdle W, Barter P, Kesäniemi YA, Mahley RW, McPherson R, Grundy SM, Wellcome Trust Case Control Consortium, Bingham SA, Khaw KT, Loos RJ, Waeber G, Barroso I, Strachan DP, Deloukas P, Vollenweider P, Wareham NJ, Mooser V: LDL-cholesterol concentrations: a genome-wide association study. Lancet. 2008, 371: 483-491. 10.1016/S0140-6736(08)60208-1.PubMed CentralView ArticleGoogle Scholar
  10. Willer CJ, Sanna S, Jackson AU, Scuteri A, Bonnycastle LL, Clarke R, Heath SC, Timpson NJ, Najjar SS, Stringham HM, Strait J, Duren WL, Maschio A, Busonero F, Mulas A, Albai G, Swift AJ, Morken MA, Narisu N, Bennett D, Parish S, Shen H, Galan P, Meneton P, Hercberg S, Zelenika D, Chen WM, Li Y, Scott LJ, Scheet PA, et al: Newly identified loci that influence lipid concentrations and risk of coronary artery disease. Nat Genet. 2008, 40: 161-169. 10.1038/ng.76.View ArticleGoogle Scholar
  11. Nakayama M, Nakajima D, Nagase T, Nomura N, Seki N, Ohara O: Identification of high-molecular-weight proteins with multiple EGF-like motifs by motif-trap screening. Genomics. 1998, 51: 27-34. 10.1006/geno.1998.5341.View ArticleGoogle Scholar
  12. NCI-NHGRI Working Group on Replication in Association Studies, Chanock SJ, Manolio T, Boehnke M, Boerwinkle E, Hunter DJ, Thomas G, Hirschhorn JN, Abecasis G, Altshuler D, Bailey-Wilson JE, Brooks LD, Cardon LR, Daly M, Donnelly P, Fraumeni JF, Freimer NB, Gerhard DS, Gunter C, Guttmacher AE, Guyer MS, Harris EL, Hoh J, Hoover R, Kong CA, Merikangas KR, Morton CC, Palmer LJ, Phimister EG, Rice JP, et al: Replicating genotype-phenotype associations. Nature. 2007, 447: 655-660. 10.1038/447655a.View ArticleGoogle Scholar
  13. Arnett FC, Edworthy SM, Bloch DA, McShane DJ, Fries JF, Cooper NS, Healey LA, Kaplan SR, Liang MH, Luthra HS, Medsger TA, Mitchell DM, Neustadt DH, Pinals RS, Schaller JG, Sharp JT, Wilder RL, Hunder GG: The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum. 1988, 31: 315-324. 10.1002/art.1780310302.View ArticleGoogle Scholar
  14. Gonzalez-Gay MA, Gonzalez-Juanatey C, Vazquez-Rodriguez TR, Martin J, Llorca J: Endothelial dysfunction, carotid intima-media thickness, and accelerated atherosclerosis in rheumatoid arthritis. Semin Arthritis Rheum. 2008, 38: 67-70. 10.1016/j.semarthrit.2008.02.001.View ArticleGoogle Scholar
  15. Ter Avest E, Stalenhoef AF, de Graaf J: What is the role of non-invasive measurements of atherosclerosis in individual cardiovascular risk prediction?. Clin Sci (Lond). 2007, 112: 507-516. 10.1042/CS20060228.View ArticleGoogle Scholar


© López Mejías et al.; licensee BioMed Central Ltd. 2012

This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.