Abstract
Objective. C677T and A1298C polymorphisms in the enzyme methylenetetrahydrofolate reductase (MTHFR) have been associated with increased cardiovascular (CV) events in non-rheumatoid arthritis (RA) populations. We investigated potential associations of MTHFR polymorphisms and use of methotrexate (MTX) with time-to-CV event in data from the Veterans Affairs Rheumatoid Arthritis (VARA) registry.
Methods. VARA participants were genotyped for MTHFR polymorphisms. Variables included demographic information, baseline comorbidities, RA duration, autoantibody status, and disease activity. Patients’ comorbidities and outcome variables were defined using International Classification of Diseases-9 and Current Procedural Terminology codes. The combined CV event outcome included myocardial infarction (MI), percutaneous coronary intervention, coronary artery bypass graft surgery, and stroke. Cox proportional hazards regression was used to model the time-to-CV event.
Results. Data were available for 1047 subjects. Post-enrollment CV events occurred in 97 patients (9.26%). Although there was a trend toward reduced risk of CV events, MTHFR polymorphisms were not significantly associated with time-to-CV event. Time-to-CV event was associated with prior stroke (HR 2.01, 95% CI 1.03–3.90), prior MI (HR 1.70, 95% CI 1.06–2.71), hyperlipidemia (HR 1.57, 95% CI 1.01–2.43), and increased modified Charlson-Deyo index (HR 1.23, 95% CI 1.13–1.34). MTX use (HR 0.66, 95% CI 0.44–0.99) and increasing education (HR 0.87, 95% CI 0.80–0.95) were associated with a lower risk for CV events.
Conclusion. Although MTHFR polymorphisms were previously associated with CV events in non-RA populations, we found only a trend toward decreased association with CV events in RA. Traditional risk factors conferred substantial CV risk, while MTX use and increasing years of education were protective.
- RHEUMATOID ARTHRITIS
- METHOTREXATE
- SINGLE-NUCLEOTIDE POLYMORPHISM
- METHYLENETETRAHYDROFOLATE REDUCTASE
- CARDIOVASCULAR DISEASES
In rheumatoid arthritis (RA), cardiovascular (CV) disease represents the most important cause of mortality1. Further, the rate of CV disease in RA exceeds that of the general population2. Traditional risk factors such as age, sex, smoking, hypertension, and diabetes have not fully explained this increased risk of CV disease, leading investigators to search for novel biologic or genetic explanations3. Some of these theories include an increased genetic susceptibility inherent in RA, enhanced systemic inflammation, or a modification of traditional CV risk factors by RA.
Methotrexate (MTX) is the “cornerstone” drug in treatment of RA4. MTX affects the folic acid pathway, and has been associated with hyperhomocysteinemia5. Methylenetetrahydrofolate reductase (MTHFR) is an enzyme that also participates in the folic acid pathway6 and catalyzes the irreversible conversion of 5,10-methylenetetrahydrofolate to 5-methyltetrahydrofolate. Two MTHFR gene polymorphisms, C677T and A1298C [single-nucleotide polymorphism (SNP) identification numbers rs1801133 and rs1801131, respectively], result in an enzyme with reduced activity6. These MTHFR polymorphisms have been linked to the occurrence of MTX adverse drug reactions in a number of reports7,8,9,10.
MTHFR polymorphisms have also been implicated in the occurrence of CV disease in the general population. Specifically, patients with the MTHFR C677T polymorphism may have an increased risk of CV disease11, premature myocardial infarction (MI) in the setting of angiographically normal coronary arteries12, and ischemic stroke11,13. Further, individuals who are homozygous for MTHFR C677T polymorphism demonstrate increased levels of homocysteine6, and elevated concentrations of this amino acid may impart additional risk for CV and cerebrovascular disease11.
Recognizing the biologic plausibility of these associations, reports linking MTHFR polymorphisms with CV risk have been inconsistent in their findings. Of 2 metaanalyses examining MTHFR C677T in non-RA populations, one supported an increased risk of CV disease14, while the other did not15. In RA — a condition known to accelerate CV disease — a single modest-size case-control study (n = 612 RA cases) investigating MTHFR polymorphisms and CV disease supported an association of MTHFR A1298C with increased CV risk16. However, that study did not control for important covariates such as MTX use.
Because CV disease is the most common cause of morbidity and mortality in patients with RA, MTHFR polymorphisms may carry an increased risk of CV events, and MTX may further complicate the effects of the MTHFR polymorphisms, we hypothesized that RA patients with MTHFR C677T and A1298C polymorphisms would be at greater risk for CV events (decreased time-to-CV event) compared to RA patients without these polymorphisms, after controlling for important risk factors.
MATERIALS AND METHODS
Study design
We evaluated a prospective observational registry supplementing with data recorded in patients’ electronic medical records.
Setting
This substudy of the Veterans Affairs Rheumatoid Arthritis (VARA) registry included data from 6 US Department of Veterans Affairs medical centers: Denver, CO; Omaha, NE; Salt Lake City, UT; Washington, DC; Dallas, TX; and Jackson, MS. VARA is a prospective longitudinal clinical registry and biologic repository. Details of the registry are described elsewhere17,18.
Participants
Inclusion criteria consisted of patient participation in the VARA registry, patients meeting the 1987 American College of Rheumatology criteria for RA19, and the availability of genotyping data for MTHFR C677T and/or A1298C polymorphisms (n = 1047).
Each participating VARA site obtained Institutional Review Board approval and all subjects underwent an informed consent procedure. The Colorado Multiple Institutional Review Board granted approval for this VARA substudy, as did the VARA Scientific Ethics Advisory Committee.
Genotyping
DNA samples were derived from whole blood for all study subjects collected at enrollment into VARA. Patients were genotyped for the MTHFR C677T (rs1801133) and A1298C (rs1801131) polymorphisms. For amplification of A1298C, a TaqMan assay was performed on a GeneAmp 9700 machine (Applied Biosystems) with endpoint analysis on a Prism 7900HT Sequence Detection System (Applied Biosystems). For C677T, the genotype was generated by a BeadExpress platform (Illumina). These 2 markers were part of a larger set of 113 markers run on these 2 platforms. Ninety-one samples (8.7%) were run as duplicates to assess genotype quality. Call rates among duplicate samples were 99.98% concordant. In addition, 4 markers were run on both platforms, yielding an inter-platform concordance of 99.95%.
Administrative and clinical data
All administrative data were derived from the Veterans Affairs inpatient and outpatient patient treatment files. Patient characteristics and RA severity measures were derived from the VARA clinical database.
Predictor variables
We assembled covariates from 3 domains: patient characteristics, RA severity/disease activity measures, and genetic characteristics, as follows.
-
Patient and exposure characteristics included age; sex; race (white vs non-white); education level (years); tobacco use (never, former, current) upon enrollment; MTX exposure at time of enrollment (yes/no); glucocorticoid ever-user (yes/no) during the time of enrollment; body mass index (BMI) upon enrollment (kg/m2); a prevalent diagnosis of alcohol abuse, diabetes mellitus, hypertension, or hyperlipidemia [based on International Classification of Diseases, 9th Revision, Clinical Modification (ICD-9-CM) in any coding position in inpatient and outpatient records] in the 7 years prior to enrollment in VARA (Appendix 1); urban vs rural residence (defined using the patient’s zip code and metropolitan statistical area); and modified Charlson-Deyo comorbidity index (excluding diabetes mellitus and CV codes, which were modeled separately)20. We included covariates to indicate a personal history of CV-related conditions/events in the 7 years prior to registry enrollment. These conditions/events consisted of MI, heart failure (HF), peripheral vascular disease (PVD), percutaneous coronary intervention (PCI), coronary artery bypass graft (CABG), or stroke. We defined these baseline comorbidities using ICD-9-CM codes and Current Procedural Terminology (CPT) codes in any position in the inpatient and outpatient VA patient treatment files occurring in the 7 years prior to enrollment in the VARA registry. The coding for each condition/event was based on previously validated algorithms, when available21,22, and if unavailable, algorithms appearing in peer-reviewed publications23,24,25,26. ICD-9-CM and CPT codes used to define the baseline comorbidities are shown in Appendix 1.
-
RA severity/disease activity characteristics consisted of RA duration prior to enrollment (years); radiological changes present in hand, wrist, or foot radiographs (ever); mean 28-joint 4-variable Disease Activity Score (DAS-28) using erythrocyte sedimentation rate (ESR) since the time of enrollment27; rheumatoid factor (RF) positivity (ever); anticitrullinated protein antibody (ACPA) positivity (ever); C-reactive protein (CRP) level at enrollment; and presence of rheumatoid nodules (ever). For this VARA substudy, antibody results were based upon the clinical laboratory results reported at each study site for participating subjects.
-
Genetic characteristics were evaluated by a genotypic analysis with variables for MTHFR 677 TT and 677 CT and MTHFR 1298 CC and 1298 AC. Additionally, we used HLA-DRB-1 shared-epitope (SE) status as a control variable28, and performed DRB-1 genotyping as described29.
Outcome variable
We defined time-to-CV event as the number of days from enrollment into VARA until the first occurrence of any one of the following conditions/events (both fatal and nonfatal): MI, stroke, PCI, or CABG. These events were defined by ICD-9-CM or CPT coding in the primary coding position of either inpatient or outpatient treatment files for MI (ICD-9-CM 410.x) and stroke (ICD-9-CM 433.11, 434.91, 435.x, 438.x), or in any coding position for PCI (ICD-9-procedure: 36.06, 36.07, 0.66; CPT: 92973, 92980, 92995) or CABG (ICD-9-procedure: 36.1x)30. If no event was identified, subjects were censored 1 day after their last inpatient or outpatient observation.
Statistical methods
For missing variables, education (163 subjects) and BMI (103 subjects), we used the SAS Multiple Imputation procedure to impute missing values based on race, sex, and age. Twenty imputation datasets were created and the average value for each subject was used. Patient characteristics for subjects sustaining a CV event and those without a CV event were compared using the Student t test and chi-square test where appropriate. Data were tested such that they met the assumptions of the Cox proportional hazards model by graphing the hazards. Cox proportional hazards regression was used to model time-to-CV event, with the time beginning with enrollment in VARA. We first identified significant predictor variables using backward stepwise selection and an alpha level of 0.2. All predictors identified in the initial model were included in the final backward stepwise selection using a 2-tailed alpha level of 0.05. As the variables of interest, A1298C and C677T polymorphisms were forced into the final model. We also tested for interaction terms between MTHFR polymorphisms and MTX use.
Secondary analyses
Two secondary analyses were performed. The first secondary analysis used Cox proportional hazards regression to predict time-to-CV event, but this analysis examined only those participants without a prior CV event (MI, stroke, CABG, PCI). Because of lower numbers of observations, the MTHFR polymorphisms were analyzed in a dominant model (0 or ≥ 1 copy of the minor allele) for this particular analysis. The second analysis used Cox proportional hazards regression to identify significant predictors of time-to-event all-cause mortality (as specific causes of death were not known). Both models were performed using the methods described above. All analyses were performed using SAS version 9.2 or Stata version 11.
RESULTS
Descriptive data and univariate analysis
Patient characteristics for those included in the analysis are summarized in Table 1. Data were available for 1047 patients covering a mean observation period of 1550 days, or an average of 4.25 years, with a total of 3743 person-years of followup. From this cohort, 97 patients (9.26%) sustained CV events post-enrollment, producing a crude event rate of 2.59 events/100 person-years (95% CI 2.11–3.15). Of the 1047 patients observed, 252 (24.07%) died during the followup period. The cause of death was not known (i.e., all-cause mortality). However, of these 252 patients, 24 (2.29%) died within 30 days of the CV event.
The preponderance of the cohort was male, white, former or current smokers, having a diagnosis of hypertension and hyperlipidemia, with well-established RA. Most were seropositive for RF (89.1%) or ACPA (81.4%). Compared to those without a CV event occurring during observation, the CV event group included a higher percentage of males, lower education status, higher rates of hypertension and hyperlipidemia diagnoses, higher modified Charlson-Deyo index, and higher rate of prior CV events (combined, MI, HF, PVD, and stroke). The CV event group also contributed a longer period of observation (1738 vs 1530 days; p < 0.0001) and demonstrated a higher rate of CC or CT polymorphisms in MTHFR C677T (homozygous major allele and heterozygous allele, respectively).
The homozygous minor allele genotypes, MTHFR 1298 CC and 677 TT, were present in 10.8% and 10.1% of the cohort, respectively (Table 1). While the minor allele frequencies did not differ between the CV event and no-CV event groups for MTHFR A1298C (rs1801131), there was a statistically significant difference between the 2 groups for MTHFR C677T, with a lower frequency of the homozygous minor (TT) genotype in those experiencing a CV event compared to those in the no-CV event group (p = 0.023). The cohort, stratified by ethnicity31, were in Hardy-Weinberg equilibrium for both MTHFR A1298C (white: Pearson chi-square = 3.70, p = 0.054; African American: Pearson chi-square = 0.053, p = 0.818) and MTHFR C677T (white: Pearson chi-square = 0.107, p = 0.744; African American: Pearson chi-square = 0.329, p = 0.566).
Multivariate analysis
Results of the multivariate Cox proportional hazards regression are shown in Table 2. The MTHFR A1298C and C677T variables were not significant in the genetic domain model, but were forced into the final model. In the final model, although MTHFR 677 minor allele homozygotes trended toward significance (HR 0.24, 95% CI 0.06–1.02, p = 0.053), the variable did not achieve significance of p < 0.05. The other MTHFR polymorphisms were nonsignificant. Interaction terms between the MTHFR polymorphisms and MTX were also not statistically significant (data not shown). Prior stroke (HR 2.01, 95% CI 1.03–3.90), prior MI (HR 1.70, 95% CI 1.06–2.71), hyperlipidemia (HR 1.57, 95% CI 1.01–2.43), and increasing modified Charlson-Deyo index (HR 1.23, 95% CI 1.13–1.34) were associated with a higher hazard ratio of a CV event. MTX use (HR 0.66, 95% CI 0.44–1.00) and higher education status (HR 0.87, 95% CI 0.80–0.95) were associated with a lower hazard ratio for CV events.
There was no significant difference between the survival curves of the major homozygous forms (MTHFR 1298 AA and 677 CC) in comparison with MTHFR 1298 CC and 677 TT minor homozygotes.
Secondary analyses
In secondary analyses, we first excluded those who had a history of CV events and performed backward stepwise analysis, as described above. Data were available for 844 patients with a total of 3029 person-years. From this cohort, 55 patients (6.52%) sustained CV events post-enrollment, producing a crude event rate of 1.82 events/100 person-years (95% CI 1.39–2.36). In the multivariate Cox analysis (Table 3), increasing modified Charlson-Deyo index (HR 1.26, 95% CI 1.12–1.42, p < 0.0001) and ACPA positivity (HR 4.59, 95% CI 1.10–19.08, p = 0.036) were associated with a higher hazard ratio of a CV event, while increasing years of education (HR 0.84, 95% CI 0.76–0.94, p = 0.003) and prevalent diagnosis of diabetes (HR 0.43, 95% CI 0.19–0.96, p = 0.039) were associated with a lower hazard ratio for CV events. MTHFR 1298 and 677 minor alleles were not significant in this model (HR 1.28, 95% CI 0.71–2.32, p = 0.414; HR 1.07, 95% CI 0.59–1.93, p = 0.832, respectively).
In a second secondary analysis, we examined predictors of time-to-all-cause mortality. Data were available for 1047 patients, with a total of 3800 person-years. From this cohort, 252 (24.07%) patients died during the followup period, producing a crude rate of 6.63 deaths/100 person-years (95% CI 5.86–7.50). In the multivariate Cox analysis (Table 4), variables associated with increased hazard ratio of death were increasing age (HR 1.08, 95% CI 1.06–1.09, p < 0.0001), increasing modified Charlson-Deyo index (HR 1.13, 95% CI 1.07–1.20, p < 0.0001), increasing DAS-28 (HR 1.38, 95% CI 1.23–1.55, p < 0.0001), increasing CRP (HR 1.01, 95% CI 1.00–1.01, p = 0.011), current smoking (HR 1.93, 95% CI 1.26–2.96, p = 0.003), and prevalent diagnosis of heart failure (HR 1.52, 95% CI 1.13–2.06, p = 0.006). MTX use (HR 0.75, 95% CI 0.58–0.97, p = 0.027) and increasing BMI (HR 0.96, 95% CI 0.93–0.98, p = 0.002) were associated with decreased hazard ratio of all-cause mortality.
DISCUSSION
In our study, we used a prospective registry of US veterans with RA to study potential association between MTHFR polymorphisms and time-to-CV event. After adjusting for multiple patient demographic and RA severity measures, we found that traditional risk factors for CV disease, including hyperlipidemia, prior MI, and prior stroke, were associated with an increased hazard ratio for time-to-combined CV event (MI, stroke, PCI, or CABG), while MTX use was associated with a lower risk for time-to-CV event. We did not find a statistically significant association of the MTHFR polymorphisms A1298C or C677T with time-to-CV events, but rather a trend toward the minor alleles being protective. In 2 secondary analyses, the first excluding those with prior CV events and the second examining all-cause mortality, we again did not find association of the MTHFR polymorphisms with the outcomes of interest (time-to-CV event and time-to-all-cause mortality).
CV disease is of high concern to patients and rheumatologists because it results in greater mortality among patients with RA than any other cause of death1. Several studies have examined the associations between CV risk factors, markers of RA severity, and atherosclerosis. However, most of these studies have not used CV events as an outcome, but rather have used proxies such as carotid intimal thickness and carotid artery plaque3,32. In this study, we chose to use time-to-CV event as the outcome variable, as we believe that this represents a more clinically durable endpoint.
MTHFR polymorphisms are an example of a potentially additive nontraditional CV risk in patients with RA. In the general population, MTHFR C677T has been shown to increase homocysteine6 and has been associated with increased risks of CV disease11, cerebrovascular disease11, premature MI12, and ischemic stroke11,13. It has been theorized that these multiple CV events may result from increased circulating homocysteine levels12. MTX, which interacts with the folic acid pathway, has been shown to improve CV outcomes in patients with RA33,34; however, the interaction between MTHFR polymorphisms, RA, and MTX use has not previously been delineated.
We did not find convincing support for our hypothesis — that MTHFR homozygous polymorphisms are associated with time-to-CV event — nor did we find a decreased time-to-CV event in patients with the MTHFR polymorphisms who were taking MTX. In a posthoc power analysis based on a Cox regression model, a sample of 1047 observations achieved 62% power at an alpha level of 0.05 to detect a hazard ratio of 1.5. Under similar assumptions, 1047 observations achieved 98% power to detect a hazard ratio of 2.0. Thus, this study was well powered to detect moderate-size effects. As our study did not detect an association between MTHFR status and CV events, if there were a true effect, it would likely be quite modest.
To date, only 1 publication has focused specifically on MTHFR polymorphisms in RA patients and CV events16. In that prospective observational study, there was an association between CV events and endothelial dysfunction related to the MTHFR A1298C genotype, but not with C677T. The study cohort differed from our population. The mean age was younger (mean age 53 yrs), fewer patients with RA were followed (n = 612), and the period of followup (14 yrs) was longer. In addition, because the cohort was from Spain, the genetic background of the participants may be different from our US population. Although our 2 studies have arrived at divergent conclusions, with certain considerations, our findings may be similar. The association between the MTHFR A1298C minor allele and CV disease may be associated with early-onset disease, which was not observed in our older cohort. An additional potential difference is that of folic acid supplementation. While we could not accurately track folic acid supplementation in the VA system (as it is typically given over the counter), US cereal products are widely supplemented with folic acid, while European cereal products are not fortified35. This background folic acid supplementation could potentially blunt the possible effects of the MTHFR polymorphisms. The differences in our studies also could be explained by methodological differences, CV outcome variable, sample sizes, or variance in the underlying populations. Replication of these studies in other populations may resolve these questions.
Not surprisingly, CV events were associated with previous MI, prior stroke, hyperlipidemia, and increasing modified Charlson-Deyo comorbidity indices. Previous MI or stroke and hyperlipidemia are traditional CV risk factors3 and are consistent with findings in the general population. These associations are consistent with the existing literature and lend face validity to our findings. As the Charlson-Deyo comorbidity index is a measure of multiple comorbidities, it is unsurprising that an increasing measure of comorbidities would be associated with increasing risk of CV disease. The decreased risk of CV events associated with MTX use might reflect the antiinflammatory effects of MTX, as documented33,34. These reports suggest that MTX not only reduces the inflammation in RA, but improves associated diseases, such as atherosclerosis.
The decreased risk of CV events associated with increasing years of education may be explained in one of many ways. One possibility is that as education status increases, patients become more “health literate” and are thus more likely to adhere to medication regimens, diets, and exercise as prescribed by their doctors36. Alternatively, education status may be a proxy marker for income status. In a recent report by Maynard, et al37, low socioeconomic status was associated with increased CV risk factors and related outcomes in patients with systemic lupus erythematosus. While this may not apply to our inflammatory disease of interest, RA, it is a point to be considered and further investigated.
Of interest, we did not find an association of CV events with measures of disease activity. This is similar to other studies38,39 that found no association of measures of disease activity with carotid intimal thickness, and found no association between disease activity and CV events in a case-control study, respectively.
In our secondary analyses, we first examined time-to-CV event in those without prior CV events. The findings in this secondary analysis were similar to the primary analysis in that increasing modified Charlson-Deyo index was associated with an increased hazard ratio and increasing education was associated with a decreased hazard ratio of CV event. After excluding those with prior CV events, those with positive ACPA had an increased risk of a CV event, which is similar to findings by Hjeltnes, et al40, who associated ACPA with impaired endothelial function. Interestingly, prevalent diagnosis of diabetes was associated with a decreased hazard ratio of CV event. This may be due to preferential removal of “high-risk” (i.e., poorly controlled) diabetic patients from the cohort through removal of those with prior CV events, leaving only those patients with well-controlled disease, who are appropriately managed for CV risk factors.
In the second secondary analysis, we examined time to all-cause mortality. As expected, increasing age, increasing modified Charlson-Deyo index, prevalent diagnosis of heart failure, and current smoking were associated with increased hazard ratio of all-cause mortality. Elevated CRP has been associated with increased risk of all-cause mortality in men with cancer41 and with early mortality in hospitalized nursing home residents42, and thus is not unexpected. Improvement in DAS-28 has previously been associated with reduction in risk of CV event43, and it may be presumed that the reverse may be true, similar to what we found in our study. MTX use was associated with decreased mortality, similar to the findings of Ajeganova, et al in an early disease cohort43. Of note, increasing BMI was associated with decreased mortality. This may be because our outcome was all-cause mortality, which presumably would include cachexia-associated diseases such as cancers, etc.
Our study has several potential limitations. Because our cohort consisted primarily of elderly men with longstanding RA, it is possible that MTHFR-associated CV events preferentially occurred early in life12, thus eliminating from our study those who died from these events (survival bias). However, this population with high risk for CV events represents an excellent cohort in which to study these genetic factors if they are not related to early death. Another potential limitation is that we cannot currently comment on the amount of time required for MTX to affect the CV system. It is possible that there is an interaction between MTHFR polymorphisms and MTX after a certain cumulative dose or after a certain period of taking the medication. Additionally, as our study was primarily in a male cohort, there is the potential for a female-MTHFR interaction that we could not identify with our population. While MTHFR polymorphisms have been associated with such disparate outcomes as recurrent fetal loss44 and breast cancer45 in women, the studies that do exist that associated MTHFR polymorphisms with increased risk for CV disease in women11 were focused on women, and did not compare female to male populations. In a study that contrasts the polymorphisms in male and female populations46, the investigators, to the contrary, found an association between MTHFR C677T polymorphisms and carotid artery compliance in men, but not women. Thus, this would potentially argue that men’s CV system may be more influenced by MTHFR polymorphisms than women’s, and our population was actually a well-selected group to study for these effects. Another potential limitation is that patients may have received care outside the VA system. Patients who may not have been identified in our database are those who received care outside the VA system and died from a CV event. To contain costs, VA policy dictates patients hospitalized at non-VA facilities be transferred to VA facilities. We attempted to control for this with a variable that categorized patients as urban versus rural, under the assumption that the likelihood of being hospitalized at a VA facility may vary according to whether patients reside in rural or urban areas47. A further potential limitation was our inability to track certain medications such as over-the-counter preparations (folic acid and nonsteroidal antiinflammatory drugs; data not collected in the VA system) and certain prescription medications such as tumor necrosis factor inhibitors (data are collected in a different pharmacy system). While it is possible that these medications had a significant effect on the cardiovascular system and potentially interact with the MTHFR polymorphisms, they are not the primary emphasis of this report.
Our study took place within a comprehensive healthcare system, the Veterans Health Affairs, using a prospective registry, thereby increasing the likelihood that the relevant clinical data were captured during the period of observation. Additionally, the sample size for our study was the largest to date to evaluate CV events and MTHFR polymorphisms in patients with RA16. Finally, the outcome we examined represents a more robust clinical outcome (CV events) compared with other studies, which used proxies for CV disease.
We were able to confirm the association of CV disease with multiple established traditional CV risk factors; however, we found no evidence that MTHFR status influences CV events in patients with RA in a time-to-event analysis.
Appendix 1. ICD-9-CM, ICD-9 procedure, and CPT codes used to define baseline comorbidities and cardiovascular (CV) events
Footnotes
-
The views expressed in this article are those of the authors and do not necessarily reflect the position or policy of the Department of Veterans Affairs.
-
Supported by Veterans Affairs Health Services Research & Development (HSR&D) grant SHP 08-172 and the Division of Rheumatology, University of Colorado Denver. Dr. Davis is supported by NIH/NIAMS T32 AR007534-24, National Institutes of Health Loan Repayment Award, the American College of Rheumatology Research and Education Fund Rheumatology Scientist Development Award, and the Colorado Clinical and Translational Science Institute. Dr. Caplan is supported by VA HSR&D Career Development Award CDA 07-221.
- Accepted for publication January 25, 2013.
REFERENCES
- 1.
- 2.
- 3.
- 4.
- 5.
- 6.
- 7.
- 8.
- 9.
- 10.
- 11.
- 12.
- 13.
- 14.
- 15.
- 16.
- 17.
- 18.
- 19.
- 20.
- 21.
- 22.
- 23.
- 24.
- 25.
- 26.
- 27.
- 28.
- 29.
- 30.
- 31.
- 32.
- 33.
- 34.
- 35.
- 36.
- 37.
- 38.
- 39.
- 40.
- 41.
- 42.
- 43.
- 44.
- 45.
- 46.
- 47.