Association of Factor V Leiden With Subsequent Atherothrombotic Events: A GENIUS-CHD Study of Individual Participant Data

BACKGROUND

Studies examining the role of factor V Leiden among patients at higher risk of atherothrombotic events, such as those with established coronary heart disease (CHD), are lacking. Given that coagulation is involved in the thrombus formation stage on atherosclerotic plaque rupture, we hypothesized that factor V Leiden may be a stronger risk factor for atherothrombotic events in patients with established CHD.

METHODS

We performed an individual-level meta-analysis including 25 prospective studies (18 cohorts, 3 case-cohorts, 4 randomized trials) from the GENIUS-CHD (Genetics of Subsequent Coronary Heart Disease) consortium involving patients with established CHD at baseline. Participating studies genotyped factor V Leiden status and shared risk estimates for the outcomes of interest using a centrally developed statistical code with harmonized definitions across studies. Cox proportional hazards regression models were used to obtain age- and sex-adjusted estimates. The obtained estimates were pooled using fixed-effect meta-analysis. The primary outcome was composite of myocardial infarction and CHD death. Secondary outcomes included any stroke, ischemic stroke, coronary revascularization, cardiovascular mortality, and all-cause mortality.

RESULTS

The studies included 69 681 individuals of whom 3190 (4.6%) were either heterozygous or homozygous (n=47) carriers of factor V Leiden. Median follow-up per study ranged from 1.0 to 10.6 years. A total of 20 studies with 61 147 participants and 6849 events contributed to analyses of the primary outcome. Factor V Leiden was not associated with the combined outcome of myocardial infarction and CHD death (hazard ratio, 1.03 [95% CI, 0.92-1.16]; I2=28%; P-heterogeneity=0.12). Subgroup analysis according to baseline characteristics or strata of traditional cardiovascular risk factors did not show relevant differences. Similarly, risk estimates for the secondary outcomes including stroke, coronary revascularization, cardiovascular mortality, and all-cause mortality were also close to identity.

CONCLUSIONS

Factor V Leiden was not associated with increased risk of subsequent atherothrombotic events and mortality in high-risk participants with established and treated CHD. Routine assessment of factor V Leiden status is unlikely to improve atherothrombotic events risk stratification in this population.

PPARA Polymorphism Influences the Cardiovascular Benefit of Fenofibrate in Type 2 Diabetes: Findings From ACCORD-Lipid

The cardiovascular benefits of fibrates have been shown to be heterogeneous and to depend on the presence of atherogenic dyslipidemia. We investigated whether genetic variability in the PPARA gene, coding for the pharmacological target of fibrates (PPAR-α), could be used to improve the selection of patients with type 2 diabetes who may derive cardiovascular benefit from addition of this treatment to statins. We identified a common variant at the PPARA locus (rs6008845, C/T) displaying a study-wide significant influence on the effect of fenofibrate on major cardiovascular events (MACE) among 3,065 self-reported white subjects treated with simvastatin and randomized to fenofibrate or placebo in the ACCORD-Lipid trial. T/T homozygotes (36% of participants) experienced a 51% MACE reduction in response to fenofibrate (hazard ratio 0.49; 95% CI 0.34-0.72), whereas no benefit was observed for other genotypes (P _interaction = 3.7 × 10^-4). The rs6008845-by-fenofibrate interaction on MACE was replicated in African Americans from ACCORD (N = 585, P = 0.02) and in external cohorts (ACCORD-BP, ORIGIN, and TRIUMPH, total N = 3059, P = 0.005). Remarkably, rs6008845 T/T homozygotes experienced a cardiovascular benefit from fibrate even in the absence of atherogenic dyslipidemia. Among these individuals, but not among carriers of other genotypes, fenofibrate treatment was associated with lower circulating levels of CCL11-a proinflammatory and atherogenic chemokine also known as eotaxin (P for rs6008845-by-fenofibrate interaction = 0.003). The GTEx data set revealed regulatory functions of rs6008845 on PPARA expression in many tissues. In summary, we have found a common PPARA regulatory variant that influences the cardiovascular effects of fenofibrate and that could be used to identify patients with type 2 diabetes who would derive benefit from fenofibrate treatment, in addition to those with atherogenic dyslipidemia.

A Genetic Locus on Chromosome 2q24 Predicting Peripheral Neuropathy Risk in Type 2 Diabetes: Results From the ACCORD and BARI 2D Studies

Genetic factors have been postulated to be involved in the etiology of diabetic peripheral neuropathy (DPN), but their identity remains mostly unknown. The aim of this study was to conduct a systematic search for genetic variants influencing DPN risk using two well-characterized cohorts. A genome-wide association study (GWAS) testing 6.8 million single nucleotide polymorphisms was conducted among participants of the Action to Control Cardiovascular Risk in Diabetes (ACCORD) clinical trial. Included were 4,384 white case patients with type 2 diabetes (T2D) and prevalent or incident DPN (defined as a Michigan Neuropathy Screening Instrument clinical examination score >2.0) and 784 white control subjects with T2D and no evidence of DPN at baseline or during follow-up. Replication of significant loci was sought among white subjects with T2D (791 DPN-positive case subjects and 158 DPN-negative control subjects) from the Bypass Angioplasty Revascularization Investigation in Type 2 Diabetes (BARI 2D) trial. Association between significant variants and gene expression in peripheral nerves was evaluated in the Genotype-Tissue Expression (GTEx) database. A cluster of 28 SNPs on chromosome 2q24 reached GWAS significance (P < 5 × 10^-8) in ACCORD. The minor allele of the lead SNP (rs13417783, minor allele frequency = 0.14) decreased DPN odds by 36% (odds ratio [OR] 0.64, 95% CI 0.55-0.74, P = 1.9 × 10^-9). This effect was not influenced by ACCORD treatment assignments (P for interaction = 0.6) or mediated by an association with known DPN risk factors. This locus was successfully validated in BARI 2D (OR 0.57, 95% CI 0.42-0.80, P = 9 × 10^-4; summary P = 7.9 × 10^-12). In GTEx, the minor, protective allele at this locus was associated with higher tibial nerve expression of an adjacent gene (SCN2A) coding for human voltage-gated sodium channel NaV1.2 (P = 9 × 10^-4). To conclude, we have identified and successfully validated a previously unknown locus with a powerful protective effect on the development of DPN in T2D. These results may provide novel insights into DPN pathogenesis and point to a potential target for novel interventions.

Subsequent Event Risk in Individuals With Established Coronary Heart Disease

Supplemental Digital Content is available in the text.

Background:

The Genetics of Subsequent Coronary Heart Disease (GENIUS-CHD) consortium was established to facilitate discovery and validation of genetic variants and biomarkers for risk of subsequent CHD events, in individuals with established CHD.

Methods:

The consortium currently includes 57 studies from 18 countries, recruiting 185 614 participants with either acute coronary syndrome, stable CHD, or a mixture of both at baseline. All studies collected biological samples and followed-up study participants prospectively for subsequent events.

Results:

Enrollment into the individual studies took place between 1985 to present day with a duration of follow-up ranging from 9 months to 15 years. Within each study, participants with CHD are predominantly of self-reported European descent (38%–100%), mostly male (44%–91%) with mean ages at recruitment ranging from 40 to 75 years. Initial feasibility analyses, using a federated analysis approach, yielded expected associations between age (hazard ratio, 1.15; 95% CI, 1.14–1.16) per 5-year increase, male sex (hazard ratio, 1.17; 95% CI, 1.13–1.21) and smoking (hazard ratio, 1.43; 95% CI, 1.35–1.51) with risk of subsequent CHD death or myocardial infarction and differing associations with other individual and composite cardiovascular endpoints.

Conclusions:

GENIUS-CHD is a global collaboration seeking to elucidate genetic and nongenetic determinants of subsequent event risk in individuals with established CHD, to improve residual risk prediction and identify novel drug targets for secondary prevention. Initial analyses demonstrate the feasibility and reliability of a federated analysis approach. The consortium now plans to initiate and test novel hypotheses as well as supporting replication and validation analyses for other investigators.

Association of Chromosome 9p21 With Subsequent Coronary Heart Disease Events

Supplemental Digital Content is available in the text.

Background:

Genetic variation at chromosome 9p21 is a recognized risk factor for coronary heart disease (CHD). However, its effect on disease progression and subsequent events is unclear, raising questions about its value for stratification of residual risk.

Methods:

A variant at chromosome 9p21 (rs1333049) was tested for association with subsequent events during follow-up in 103 357 Europeans with established CHD at baseline from the GENIUS-CHD (Genetics of Subsequent Coronary Heart Disease) Consortium (73.1% male, mean age 62.9 years). The primary outcome, subsequent CHD death or myocardial infarction (CHD death/myocardial infarction), occurred in 13 040 of the 93 115 participants with available outcome data. Effect estimates were compared with case/control risk obtained from the CARDIoGRAMplusC4D consortium (Coronary Artery Disease Genome-wide Replication and Meta-analysis [CARDIoGRAM] plus The Coronary Artery Disease [C4D] Genetics) including 47 222 CHD cases and 122 264 controls free of CHD.

Results:

Meta-analyses revealed no significant association between chromosome 9p21 and the primary outcome of CHD death/myocardial infarction among those with established CHD at baseline (GENIUS-CHD odds ratio, 1.02; 95% CI, 0.99–1.05). This contrasted with a strong association in CARDIoGRAMPlusC4D odds ratio 1.20; 95% CI, 1.18–1.22; P for interaction <0.001 compared with the GENIUS-CHD estimate. Similarly, no clear associations were identified for additional subsequent outcomes, including all-cause death, although we found a modest positive association between chromosome 9p21 and subsequent revascularization (odds ratio, 1.07; 95% CI, 1.04–1.09).

Conclusions:

In contrast to studies comparing individuals with CHD to disease-free controls, we found no clear association between genetic variation at chromosome 9p21 and risk of subsequent acute CHD events when all individuals had CHD at baseline. However, the association with subsequent revascularization may support the postulated mechanism of chromosome 9p21 for promoting atheroma development.

A novel genetic marker of decreased inflammation and improved survival after acute myocardial infarction

The CHRNA5 gene encodes a neurotransmitter receptor subunit involved in multiple processes, including cholinergic autonomic nerve activity and inflammation. Common variants in CHRNA5 have been linked with atherosclerotic cardiovascular disease. Association of variation in CHRNA5 and specific haplotypes with cardiovascular outcomes has not been described. The aim of this study was to examine the association of CHRNA5 haplotypes with gene expression and mortality among patients with acute myocardial infarction (AMI) and explore potential mechanisms of this association. Patients (N = 2054) hospitalized with AMI were genotyped for two common variants in CHRNA5. Proportional hazard models were used to estimate independent association of CHRNA5 haplotype with 1-year mortality. Both individual variants were associated with mortality (p = 0.0096 and 0.0004, respectively) and were in tight LD (D' = 0.99). One haplotype, HAP3, was associated with decreased mortality one year after AMI (adjusted HR = 0.42, 95% CI 0.26, 0.68; p = 0.0004). This association was validated in an independent cohort (N = 637) of post-MI patients (adjusted HR = 0.23, 95% CI 0.07, 0.79; p = 0.019). Differences in CHRNA5 expression by haplotype were investigated in human heart samples (n = 28). Compared with non-carriers, HAP3 carriers had threefold lower cardiac CHRNA5 mRNA expression (p = 0.023). Circulating levels of the inflammatory marker hsCRP were significantly lower in HAP3 carriers versus non-carriers (3.43 ± 4.2 versus 3.91 ± 5.1; p = 0.0379). Activation of the inflammasome, an important inflammatory complex involved in cardiovascular disease that is necessary for release of the pro-inflammatory cytokine IL-1 β, was assessed in bone marrow-derived macrophages (BMDM) from CHRNA5 knockout mice and wild-type controls. In BMDM from CHRNA5 knockout mice, IL-1β secretion was reduced by 50% compared to wild-type controls (p = 0.004). Therefore, a common haplotype of CHRNA5 that results in reduced cardiac expression of CHRNA5 and attenuated macrophage inflammasome activation is associated with lower mortality after AMI. These results implicate CHRNA5 and the cholinergic anti-inflammatory pathway in survival following AMI.

A high throughput, functional screen of human Body Mass Index GWAS loci using tissue-specific RNAi Drosophila melanogaster crosses

Human GWAS of obesity have been successful in identifying loci associated with adiposity, but for the most part, these are non-coding SNPs whose function, or even whose gene of action, is unknown. To help identify the genes on which these human BMI loci may be operating, we conducted a high throughput screen in Drosophila melanogaster. Starting with 78 BMI loci from two recently published GWAS meta-analyses, we identified fly orthologs of all nearby genes (± 250KB). We crossed RNAi knockdown lines of each gene with flies containing tissue-specific drivers to knock down (KD) the expression of the genes only in the brain and the fat body. We then raised the flies on a control diet and compared the amount of fat/triglyceride in the tissue-specific KD group compared to the driver-only control flies. 16 of the 78 BMI GWAS loci could not be screened with this approach, as no gene in the 500-kb region had a fly ortholog. Of the remaining 62 GWAS loci testable in the fly, we found a significant fat phenotype in the KD flies for at least one gene for 26 loci (42%) even after correcting for multiple comparisons. By contrast, the rate of significant fat phenotypes in RNAi KD found in a recent genome-wide Drosophila screen (Pospisilik et al. (2010) is ~5%. More interestingly, for 10 of the 26 positive regions, we found that the nearest gene was not the one that showed a significant phenotype in the fly. Specifically, our screen suggests that for the 10 human BMI SNPs rs11057405, rs205262, rs9925964, rs9914578, rs2287019, rs11688816, rs13107325, rs7164727, rs17724992, and rs299412, the functional genes may NOT be the nearest ones (CLIP1, C6orf106, KAT8, SMG6, QPCTL, EHBP1, SLC39A8, ADPGK /ADPGK-AS1, PGPEP1, KCTD15, respectively), but instead, the specific nearby cis genes are the functional target (namely: ZCCHC8, VPS33A, RSRC2; SPDEF, NUDT3; PAGR1; SETD1, VKORC1; SGSM2, SRR; VASP, SIX5; OTX1; BANK1; ARIH1; ELL; CHST8, respectively). The study also suggests further functional experiments to elucidate mechanism of action for genes evolutionarily conserved for fat storage.

Human Genome Wide Association Studies have successfully found thousands of novel genetic variants associated with many diseases. While these undoubtedly point to new biology, the field has been slowed in exploiting these new findings to reach a better understanding of exactly how they confer increased risk. Many, if not most, appear to be regulatory not coding variants, so their immediate consequence is not obvious. A real rate limiting step is even identifying which gene these variants might be regulating, and in what tissues they are operating to increase disease risk. In the absence of any other information, a first order assumption is that they may be more likely to be regulating a nearby gene, and such variants are often initially annotated by the “nearest” gene until their function is more definitively validated. Exploiting the idea that many genes may have conserved function across species, we conducted a high-throughput screen of fruit-fly orthologs of human genes nearby 78 well validated GWAS variants for human obesity, in order to more precisely identify the gene(s) of action. We systematically knocked down the function of each of these nearby genes in the brain and fat-body of the flies, raised them on a standard diet, and compared their percent body fat with control flies, in order to validate which genes showed a fat response. 43% of the time when fly orthologs existed in the region, we were able to identify the causal gene. Interestingly, nearly half the time (46%), it was not the nearest gene but another nearby one that regulated fat.

Relations between lipoprotein(a) concentrations, LPA genetic variants, and the risk of mortality in patients with established coronary heart disease: a molecular and genetic association study

BACKGROUND

Lipoprotein(a) concentrations in plasma are associated with cardiovascular risk in the general population. Whether lipoprotein(a) concentrations or LPA genetic variants predict long-term mortality in patients with established coronary heart disease remains less clear.

METHODS

We obtained data from 3313 patients with established coronary heart disease in the Ludwigshafen Risk and Cardiovascular Health (LURIC) study. We tested associations of tertiles of lipoprotein(a) concentration in plasma and two LPA single-nucleotide polymorphisms ([SNPs] rs10455872 and rs3798220) with all-cause mortality and cardiovascular mortality by Cox regression analysis and with severity of disease by generalised linear modelling, with and without adjustment for age, sex, diabetes diagnosis, systolic blood pressure, BMI, smoking status, estimated glomerular filtration rate, LDL-cholesterol concentration, and use of lipid-lowering therapy. Results for plasma lipoprotein(a) concentrations were validated in five independent studies involving 10 195 patients with established coronary heart disease. Results for genetic associations were replicated through large-scale collaborative analysis in the GENIUS-CHD consortium, comprising 106 353 patients with established coronary heart disease and 19 332 deaths in 22 studies or cohorts.

FINDINGS

The median follow-up was 9·9 years. Increased severity of coronary heart disease was associated with lipoprotein(a) concentrations in plasma in the highest tertile (adjusted hazard radio [HR] 1·44, 95% CI 1·14-1·83) and the presence of either LPA SNP (1·88, 1·40-2·53). No associations were found in LURIC with all-cause mortality (highest tertile of lipoprotein(a) concentration in plasma 0·95, 0·81-1·11 and either LPA SNP 1·10, 0·92-1·31) or cardiovascular mortality (0·99, 0·81-1·2 and 1·13, 0·90-1·40, respectively) or in the validation studies.

INTERPRETATION

In patients with prevalent coronary heart disease, lipoprotein(a) concentrations and genetic variants showed no associations with mortality. We conclude that these variables are not useful risk factors to measure to predict progression to death after coronary heart disease is established.

FUNDING

Seventh Framework Programme for Research and Technical Development (AtheroRemo and RiskyCAD), INTERREG IV Oberrhein Programme, Deutsche Nierenstiftung, Else-Kroener Fresenius Foundation, Deutsche Stiftung für Herzforschung, Deutsche Forschungsgemeinschaft, Saarland University, German Federal Ministry of Education and Research, Willy Robert Pitzer Foundation, and Waldburg-Zeil Clinics Isny.

Cytochrome p450 gene variants, race, and mortality among clopidogrel-treated patients after acute myocardial infarction

BACKGROUND

Clopidogrel is recommended after acute myocardial infarction but has variable efficacy and safety, in part related to the effect of cytochrome P450 (CYP) polymorphisms on its metabolism. The effect of CYP polymorphisms on cardiovascular events among clopidogrel-treated patients after acute myocardial infarction remains controversial, and no studies to date have investigated the association of CYP variants with outcomes in black patients.

METHODS AND RESULTS

Subjects (2732: 2062 whites; 670 blacks) hospitalized with acute myocardial infarction enrolled in the prospective, multicenter TRIUMPH study were genotyped for CYP polymorphisms. The majority of whites (79%) and blacks (64.4%) were discharged on clopidogrel. Among whites, carriers of the loss-of-function CYP2C19*2 allele had significantly increased 1-year mortality (adjusted hazards ratio [HR]: 1.70; confidence interval [CI]: 1.01-2.86; P=0.046) and a trend toward increased rate of recurrent MI (adjusted HR: 2.10; CI: 0.95-4.63; P=0.066). Among blacks, increased 1-year mortality was associated with the gain-of-function CYP2C19*17 allele (adjusted HR for *1/*17 versus *1/*1: 2.02; CI: 0.92-4.44; *17/*17 versus *1/*1: 8.97; CI: 3.34-24.10; P<0.0001) and the CYP1A2*1C allele (adjusted HR for *1/*1C versus *1/*1: 1.89; CI: 0.85-4.22; *1C/*1C versus *1/*1: 4.96; CI: 1.69-14.56; P=0.014). Bleeding events were significantly more common among black carriers of CYP2C19*17 or CYP1A2*1C.

CONCLUSIONS

Both loss-of-function and gain-of-function CYP polymorphisms affecting clopidogrel metabolism are associated with increased mortality among clopidogrel-treated patients after acute myocardial infarction; the specific polymorphism and the putative mechanism vary according to race.

CHRNA5 variant predicts smoking cessation in patients with acute myocardial infarction

INTRODUCTION

While smoking is a major modifiable risk factor for secondary prevention of myocardial infarction (MI), active smoking is common among patients hospitalized with acute MI. Recent studies suggest that nicotinic receptor variants, and specifically the high-risk CHRNA5 rs16969968 A allele, are associated with cessation failure among noncardiac patients. This study investigates the association between CHRNA5 rs16969968 and smoking cessation in patients hospitalized with acute MI.

METHODS

Using data from the TRIUMPH study, we ascertained smoking status at the time of index hospitalization for acute MI and 1 year after hospitalization. After adjusting for age and sex, we used logistic regression to model the association between smoking cessation and CHRNA5 rs16969968.

RESULTS

At index admission, 752 Caucasian subjects were active smokers and 699 were former smokers. Among these ever-smokers, the A allele was associated with significantly decreased abstinence (45.0% abstinence for A allele carriers vs. 51.7% for GG homozygotes; odds ratio [OR] = 0.70, 95% confidence interval [CI] = 0.56-0.88, p = .0027). The A allele was also significantly associated with decreased abstinence at 1 year (69.1% abstinence for A allele carriers vs. 76.0% for GG homozygotes; OR = 0.70, 95% CI = 0.53-0.94, p = .0185).

CONCLUSIONS

Among patients who have smoked and who are hospitalized with acute MI, the high-risk CHRNA5 allele was associated with lower likelihood of quitting before hospitalization and significantly less abstinence 1 year after hospitalization with MI. The CHRNA5 rs16969968 genotype may therefore identify patients who would benefit from aggressive, personalized smoking cessation intervention.

Pharmacogenetic warfarin dose refinements remain significantly influenced by genetic factors after one week of therapy

By guiding initial warfarin dose, pharmacogenetic (PGx) algorithms may improve the safety of warfarin initiation. However, once international normalised ratio (INR) response is known, the contribution of PGx to dose refinements is uncertain. This study sought to develop and validate clinical and PGx dosing algorithms for warfarin dose refinement on days 6-11 after therapy initiation. An international sample of 2,022 patients at 13 medical centres on three continents provided clinical, INR, and genetic data at treatment days 6-11 to predict therapeutic warfarin dose. Independent derivation and retrospective validation samples were composed by randomly dividing the population (80%/20%). Prior warfarin doses were weighted by their expected effect on S-warfarin concentrations using an exponential-decay pharmacokinetic model. The INR divided by that "effective" dose constituted a treatment response index . Treatment response index, age, amiodarone, body surface area, warfarin indication, and target INR were associated with dose in the derivation sample. A clinical algorithm based on these factors was remarkably accurate: in the retrospective validation cohort its R(2) was 61.2% and median absolute error (MAE) was 5.0 mg/week. Accuracy and safety was confirmed in a prospective cohort (N=43). CYP2C9 variants and VKORC1-1639 G→A were significant dose predictors in both the derivation and validation samples. In the retrospective validation cohort, the PGx algorithm had: R(2)= 69.1% (p<0.05 vs. clinical algorithm), MAE= 4.7 mg/week. In conclusion, a pharmacogenetic warfarin dose-refinement algorithm based on clinical, INR, and genetic factors can explain at least 69.1% of therapeutic warfarin dose variability after about one week of therapy.

A polymorphism in the VKORC1 regulator calumenin predicts higher warfarin dose requirements in African Americans

Warfarin demonstrates a wide interindividual variability in response that is mediated partly by variants in cytochrome P450 2C9 (CYP2C9) and vitamin K 2,3-epoxide reductase complex subunit 1 (VKORC1). It is not known whether variants in calumenin (CALU) (vitamin K reductase regulator) have an influence on warfarin dose requirements. We resequenced CALU regions in a discovery cohort of dose outliers: patients with high (>90th percentile, n = 55) or low (<10th percentile, n = 53) warfarin dose requirements (after accounting for known genetic and nongenetic variables). One CALU variant, rs339097, was associated with high doses (P = 0.01). We validated this variant as a predictor of higher warfarin doses in two replication cohorts: (i) 496 patients of mixed ethnicity and (ii) 194 African-American patients. The G allele of rs339097 (the allele frequency was 0.14 in African Americans and 0.002 in Caucasians) was associated with the requirement for a 14.5% (SD +/- 7%) higher therapeutic dose (P = 0.03) in the first replication cohort and a higher-than-predicted dose in the second replication cohort (allele frequency 0.14, one-sided P = 0.03). CALU rs339097 A>G is associated with higher warfarin dose requirements, independent of known genetic and nongenetic predictors of warfarin dose in African Americans.

Laboratory and clinical outcomes of pharmacogenetic vs. clinical protocols for warfarin initiation in orthopedic patients

BACKGROUND

Warfarin is commonly prescribed for prophylaxis and treatment of thromboembolism after orthopedic surgery. During warfarin initiation, out-of-range International Normalized Ratio (INR) values and adverse events are common.

METHODS

In orthopedic patients beginning warfarin therapy, we developed and prospectively validated pharmacogenetic and clinical dose refinement algorithms to revise the estimated therapeutic dose after 4 days of therapy.

RESULTS

The pharmacogenetic algorithm used the cytochrome P450 (CYP) 2C9 genotype, smoking status, peri-operative blood loss, liver disease, INR values and dose history to predict the therapeutic dose. The R(2) was 82% in a derivation cohort (n = 86) and 70% when used prospectively (n = 146). The R(2) of the clinical algorithm that used INR values and dose history to predict the therapeutic dose was 57% in a derivation cohort (n = 178) and 48% in a prospective validation cohort (n = 146). In 1 month of prospective follow-up, the percent time spent in the therapeutic range was 7% higher (95% CI: 2.7-11.7) in the pharmacogenetic cohort. The risk of a laboratory or clinical adverse event was also significantly reduced in the pharmacogenetic cohort (Hazard Ratio 0.54; 95% CI: 0.30-0.97).

CONCLUSIONS

Warfarin dose adjustments that incorporate genotype and clinical variables available after four warfarin doses are accurate. In this non-randomized, prospective study, pharmacogenetic dose refinements were associated with more time spent in the therapeutic range and fewer laboratory or clinical adverse events. To facilitate gene-guided warfarin dosing we created a non-profit website, http://www.WarfarinDosing.org.

Optimal initial dose adjustment of warfarin in orthopedic patients

BACKGROUND

Warfarin sodium is commonly prescribed for the prophylaxis and treatment of venous thromboembolism. Dosing algorithms have not been widely adopted because they require a fixed initial warfarin dose (eg, 5 mg) and are not tailored to other factors that may affect the international normalized ratio (INR).

OBJECTIVE

To develop an algorithm that could predict a therapeutic warfarin dose based on drug interactions, INR response after the initial warfarin doses, and other clinical factors.

METHODS

We used stepwise regression to quantify the relationship between these factors in patients beginning prophylactic warfarin therapy immediately prior to joint replacement. In the derivation cohort (n = 271), we separately modeled the therapeutic dose after 2 and 3 initial doses. We prospectively validated these 2 models in an independent cohort (n = 105).

RESULTS

About half of the therapeutic dose variability was predictable after 3 days of therapy: R2 was 53% in the derivation cohort and 42% in the validation cohort. INR response after 3 warfarin doses (INR3) inversely correlated with therapeutic dose (p < 0.001). Intraoperative blood loss transiently, but significantly, elevated the postoperative INR values. Other significant (p < 0.03) predictors were the first and second warfarin doses (+7% and +6%, respectively, per 1 mg), and statin use (-15.0%). The model derived after 2 warfarin doses explained 32% of the variability in therapeutic dose.

CONCLUSIONS

We developed and validated algorithms that estimate therapeutic warfarin doses based on clinical factors and INR response available after 2-3 days of warfarin therapy. The algorithms are implemented online at www.WarfarinDosing.org.

Genetic-based dosing in orthopedic patients beginning warfarin therapy

High variability in drug response and a narrow therapeutic index complicate warfarin therapy initiation. No existing algorithm provides recommendations on refining the initial warfarin dose based on genetic variables, clinical data, and international normalized ratio (INR) values. Our goal was to develop such an algorithm. We studied 92 patients undergoing primary or revision total hip or knee replacement. From each patient we collected a blood sample, clinical variables, current medications, and preoperative and postoperative laboratory values. We genotyped for polymorphisms in the cytochrome P450 (CYP) 2C9 and vitamin K epoxide reductase (VKORC1) genes. Using stepwise regression, we developed a model for refining the warfarin dose after the third warfarin dose. The algorithm explained four fifths of the variability in therapeutic dose (R(2)(adj) of 79%). Significant (P > .05) predictors were INR value after 3 doses (47% reduction per 0.25-unit rise), first warfarin dose (+7% per 1 mg), CYP2C9*3 and CYP2C9*2 genotype (-38% and -17% per allele), estimated blood loss (interacting with INR(3)), smoking status (+20% in current smokers), and VKORC1 (-11% per copy of haplotype A). If validated, this model should provide a safer, more effective process for initiating warfarin therapy.