|
|
||||||||
Original Research |
Michael D. Caldwell, MD, PhD, Department of Surgery, Marshfield Clinic, 1000 North Oak Avenue, Marshfield, Wisconsin 54449
Richard L. Berg, MS, Biomedical Informatics Research Center, Marshfield Clinic Research Foundation, 1000 North Oak Avenue, Marshfield, Wisconsin 54449
Kai Qi Zhang, BS, Center for Human Genetics, Marshfield Clinic Research Foundation 1000 North Oak Avenue Marshfield, Wisconsin 54449
Ingrid Glurich, PhD, Office of Research Facilitation, Marshfield Clinic Research Foundation, 1000 North Oak Avenue, Marshfield, Wisconsin 54449
John R. Schmelzer, PhD, Health Services Research Center, Marshfield Clinic Research Foundation, 1000 North Oak Avenue, Marshfield, Wisconsin 54449
Steven H. Yale, MD, Department of Internal Medicine, Marshfield Clinic and Clinical Research Center, Marshfield Clinic Research Foundation, 1000 North Oak Avenue, Marshfield, Wisconsin 54449
Humberto J. Vidaillet, MD, Department of Cardiology, Marshfield Clinic and Administration, Marshfield Clinic Research Foundation, 1000 North Oak Avenue, Marshfield, Wisconsin 54449
James K. Burmester, PhD, Center for Human Genetics, Marshfield Clinic Research Foundation, 1000 North Oak Avenue, Marshfield, Wisconsin 54449
Reprint Requests: James K. Burmester, PhD, Center for Human Genetics, Marshfield Clinic Research Foundation, 1000 North Oak Avenue, Marshfield, WI 54449, Telephone: 715-389-4368, Fax: 715-389-3808, Email: burmester.jim{at}mcrf.mfldclin.edu
Received: September 26, 2006.
Revised: November 8, 2006.
Accepted: January 8, 2006.
| Abstract |
|---|
|
|
|---|
Method: A model was constructed by evaluating the contribution to dosing variability of the following clinical factors: age, gender, body surface area, and presence or absence of prosthetic heart valves or diabetes. The model was then sequentially expanded by incorporating polymorphisms of cytochrome P450 (CYP) 2C9; vitamin K 2,3 epoxide reductase complex, subunit 1 (VKORC1); gamma carboxylase; factor VII; and apolipoprotein (Apo) E genes.
Results: Of genetic factors evaluated in the model, CYP2C9 and VKORC1 each contributed substantially to dose variability, and together with clinical factors explained 56% of the individual variability in stable warfarin dose. In contrast, gamma carboxylase, factor VII and Apo E polymorphisms contributed little to dose variability.
Conclusion: The importance of CYP2C9 and VKORC1 to patient-specific dose of warfarin has been confirmed, while polymorphisms of gamma carboxylase, factor VII and Apo E genes did not substantially contribute to predictive models for stable warfarin dose.
Key Words: Bleeding Coagulation Vitamin K
Warfarin is an inexpensive, commonly prescribed anticoagulant that inhibits activity of vitamin K dependent blood coagulation factors. Warfarin is prescribed for the treatment of deep venous thrombosis and pulmonary embolism and also prophylactically for prevention of thromboembolic events.13
Despite its proven efficacy and low cost, warfarin is underutilized among individuals who could benefit from this therapy due to the high rate of adverse events, the lack of precise warfarin dosing parameters and the need to monitor patients on a regular basis.46 Strategies to more reliably achieve optimal anticoagulation and maintain patients in therapeutic range would increase the use of warfarin by practitioners and offer the potential to improve patient safety and reduce adverse events.
Warfarin is a racemic mixture of S-warfarin and R-warfarin. S-warfarin is approximately 5 times more potent than R-warfarin. Metabolism of S-warfarin occurs through the cytochrome P450 (CYP) 2C9 enzyme, while metabolism of R-warfarin occurs through CYP1A2 and CYP3A4. Three prominent genetic variants have been identified in CYP2C9, resulting in reduced enzymatic activity for warfarin metabolism. The *2/*2 homozygous mutant leads to a reduction to approximately 12% of CYP2C9 activity and the *3/*3 homozygous mutant has <5% of wild type CYP2C9 activity.7 Additionally, dosing variability observed among subjects with at least one CYP2C9 *1 allele may be partly attributable to genetic polymorphisms in the C1 subunit of the vitamin K 2,3 epoxide reductase complex (VKORC1).813 This enzyme has recently been identified as the rate-limiting step in the warfarin sensitive, vitamin K-dependent gamma carboxylation system.14 Because the initiation of warfarin treatment is problematic and is dependent on complex interactions between physiological, environmental and genetic factors, we previously developed a warfarin-dosing algorithm designed to more accurately initiate warfarin therapy.15 That algorithm predicts an initiation dose that more closely approximates stable dosing requirements by formally incorporating the contributions of genetic variants in CYP2C9. In this study, we extend our earlier model to incorporate variability in stable dose contributed by VKORC1, and we investigate the putative contribution of known polymorphisms in gamma carboxylase, factor VII and apolipoprotein (Apo) E.
| Materials and Methods |
|---|
|
|
|---|
Patients were excluded from this study if they were known to have underlying conditions that influence drug dosage, such as cancer, renal or hepatic insufficiency, or congestive heart failure. Seven patients taking CYP2C9 inducer medications were excluded.
The majority of patients (533/570) were recruited for a prior study of the impact of CYP2C9 genotype on warfarin dosing, and data from 453 of these participants were reported previously.15 Additional patients (37/570) were recruited for a prospective study,17 which evaluated the predictive model developed in the first study in a randomized controlled pilot trial. A subset of the study subjects (501/570) subsequently consented to additional genetic testing to evaluate other promising genetic markers. All study subjects were consented and enrolled according to institutional protocols approved by the Institutional Review Board of Marshfield Clinic. Age, gender, body surface area, clinical attributes such as presenting condition (e.g., valve replacement), presence or absence of diabetes, anticoagulation data and prescription medications were abstracted from the clinical records of patients.
Genetic Testing
The CYP2C9 assay was performed as previously described.15 For VKORC1, we tested single nucleotide polymorphism (SNP) rs8050894 that is described in dbSNP (http://www.ncbi.nlm.nih.gov/SNP/index.html). We performed this testing prior to publications that tested other SNPs in VKORC1. This SNP corresponds to position 6853 of GenBank sequence AY587020 and is a variant that tags the major haplotypes present in VKORC1.12 A validated TaqMan assay for detection of this polymorphism was purchased from Applied Biosystems (Foster City, CA), catalog number C_2847860_10.
Testing of Apo E was performed using TaqMan genotyping assays purchased from Applied Biosystems. These assays tested for the SNPs that encode Arg to Cys amino acid substitutions at position 112 in rs429358 and position 158 in rs7412 (catalog numbers C_3084793_20 and C_904973_10, respectively). Possible haplotypes include E2=Cys(112), Cys(158); E3=Cys(112), Arg(158); E4= Arg(112), Arg(158). These data were summarized for statistical analysis as the number of E4 alleles.
Testing of factor VII was done to determine if variants at positions 402 and 401 of its promoter region affect warfarin dosing.18,19 Since these positions are adjacent to each other, it is difficult to design direct assays for these positions where the testing of one polymorphism would not be hindered by variation at the adjacent polymorphism. To overcome this testing problem, we assayed the SNP at position 1704 rs762635, which is in complete linkage disequilibrium with position 1972 (aka 402) rs570317, and the SNP at position 2456 rs6039 that tags the common haplotypes for position 1973 (aka 401) rs7981123 through linkage disequilibrium with a discordance of 2%.20 Custom designed TaqMan assays were used for these SNP analyses (a TaqMan assay for rs6039 is now commercially available from Applied Biosystems under catalog number C_783182_10). The two SNPs were considered separately as potential predictors of warfarin metabolism for statistical analyses.
To test the microsatellite polymorphism in gamma carboxylase (rs10654848), we performed polymerase chain reaction (PCR) amplification followed by electrophoresis through a 6% acrylamide gel and analysis on a PhosphorImager (Molecular Dynamics, Sunnyvale, CA). Primers used for amplification were 5'-agtcagctcctcagacaaca and 5'-cagctaggcgtcggtggca. Thermal cycling conditions were 95°C for 5 min, followed by 35 cycles of 94°C for 1 min, 50°C for 1 min, 72°C for 1 min with a final extension at 72°C for 6 min. Allele sizes of 10, 11, 13 and 14 were confirmed by sequencing PCR products and were used as standards for gel electrophoresis. For statistical analysis, repeat size was summarized by the maximum allele size up to 14 with alleles of 10,13 or 11,13 classified as 13 (e.g., 10,10 was classified as 10; and 10,14 was classified as 14).
Statistical Analysis
The clinical parameters considered here were those reported previously:15 age, gender, body surface area, treatment indication (heart valve or not), and presence or absence of diabetes. Most of our subjects had an INR target range of 2 to 3, with the exception of heart valve patients for whom the range is generally 2.5 to 3.5. This association implies that the heart valve parameter in our models is to some degree a surrogate for INR target range.
Because not all study subjects consented to additional genetic testing, and because genetic test results are sometimes inconclusive for any given DNA sample, fewer than 570 observations were available for analyses of genetic data other than CYP2C9. In all analyses, all informative cases were included, and a consistent model fitting strategy was applied. Whenever models involving different predictor variables were compared, the same study subjects were used in each analysis. For example, when comparing a model that included clinical parameters and CYP2C9 to a model with only clinical parameters, all 570 subjects were included, but only 431 subjects with valid VKORC1 results were used to evaluate the effects of adding VKORC1 to the set of clinical and CYP2C9 predictors.
All models were fit with the log of stable warfarin dose as the response. Clinically, it was felt that the effect of warfarin was best reflected by a proportional (multiplicative) scale.21 This implies, for example, that the adverse effects of 15 mg/week in a patient needing 10 mg/week would be equivalent to the adverse effects of 60 mg/week in another patient needing 40 mg/week, since in each case the dose is 50% higher than needed. Statistically, model assumptions of constant variance and normality were observed to be better met under models fit to the log of dose.
Although there is ample evidence in the literature of an important association between CYP2C9 and stable warfarin dose, even this large study did not include a sufficient number of subjects to reliably estimate the independent impact of clinical and genetic predictors of stable dose in those with no wild type CYP2C9 alleles (e.g., no *1 alleles). Therefore, we conducted multiple regression analyses only on the *1 genotypes and used the observed mean dose (log scale) as the predicted value for those with only non-wild CYP2C9 alleles (*2/*2, *2/*3, and *3/*3). Accordingly, our models that include genetic variables are composites of multiple regression and simple mean-based estimation approaches.
We employed a modified step-down approach in developing our regression models by first estimating a "full" model. Full models included all potential clinical and genetic predictors pertinent to that model, plus all two-factor interactions and quadratic terms for age and body surface area. Terms were then cautiously excluded from this full model only if they contributed very little. The regression models were finalized after evaluating residuals and excluding cases with large residuals as outliers (|Studentized Residual| >3). At most, 5 cases (<1% of the cohort) were excluded as outliers in the genetic models.
In order to evaluate potential improvements to dosing based solely on clinical factors like age and body size, a series of genetic models was evaluated. Since subject numbers were too limited to consider all genetic factors simultaneously, models were built sequentially by modeling dose with and without each of the individual genetic terms under consideration. The evidence in the literature is strongest for CYP2C9 and VKORC1; therefore, these were the first factors considered. The explanatory power for competing models was assessed graphically and by using the adjusted R-squared (R2adj) statistic which is based on the squared difference between the observed and predicted dose (log scale). The R2adj statistic measures the proportion of the total variability explained by the model but with adjustment for the number of parameters in the model. The statistic was directly calculated to use all cases under analysis by pooling across both parts of the composite models and by applying the final regression model to any outliers dropped when estimating the parameters in the regression models.
| Results |
|---|
|
|
|---|
|
|
|
We depict the impact of CYP2C9, VKORC1 and the additional genetic variants in figure 2
. The paired plots in this figure show changes in the relationship between observed stable dose and dose predicted by the model as we serially add each genetic factor. CYP2C9 contributed substantially and significantly (P<0.001) to a predictive model for stable warfarin dose which included only clinical factors (figure 2
, panels AB), increasing the R2adj from 15% to 28%. VKORC1 contributed substantially and significantly (P<0.001) to a predictive model for stable warfarin dose which included both CYP2C9 and clinical factors (figure 2
, panels CD), increasing the R2adj from 34% (in the subset with valid VKORC1 results) to 56%. Our results confirm that both VKORC1 and CYP2C9 are important factors associated with stable warfarin dose.
|
The results from our model could be used to construct visual aids, such as nomograms, for personalized warfarin dosing. Figure 3
illustrates a series of nomograms for predicting the warfarin dose based on a simplified version of our model which includes the patients CYP2C9 and VKORC1 genotype, age and gender. Each nomogram represents a different CYP2C9 genotype with separate results presented for each VKORC1 genotype, specified as homozygous GG, heterozygous CG or homozygous CC. Patient age is referenced on the horizontal scale in the nomograms.
|
| Discussion |
|---|
|
|
|---|
Our model explains about 56% of the variation in stable dose. We acknowledge that this model requires validation in other populations. We expect some reduction in predictive ability, even in our own population, for cases not used in the model development. We also note a relatively small number of our current patients for whom dose is not well predicted by the model. It is possible that these patients carry additional genetic mutations in genes involved in the coagulation pathway. We fully expect that the model will continue to evolve as we discover new factors that contribute to variability in dosing and acquire additional data from cases with rare genotypes.
We have illustrated dosing nomograms that are simple to use but allow for a more accurate and personalized approach to warfarin dosing. These nomograms graphically display and highlight the wide variation in warfarin stable dose, and allow physicians the ability to visually determine a more precise dose when genetic, gender and age parameters are known. Our ultimate goal is to develop a computer-based application that will be both easy to use and easy to update as the model is improved.
As methods for personalized dosing of warfarin are improved, it is likely that physician comfort with this drug will increase. Although the use of genetic testing and dosing algorithms is currently feasible in the United States, these tools are not currently available to clinical practice. Two studies have demonstrated the feasibility of genetic testing prior to warfarin dosing.17,37 Both studies concluded that genetic testing is technically feasible, acceptable to the patient and physician, and likely to improve the clinical management of anticoagulation.
| Acknowledgments |
|---|
|
|
|---|
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. W. Linder, M. Bon Homme, K. K. Reynolds, B. F. Gage, C. Eby, N. Silvestrov, and R. Valdes Jr. Interactive Modeling for Ongoing Utility of Pharmacogenetic Diagnostic Testing: Application for Warfarin Therapy Clin. Chem., October 1, 2009; 55(10): 1861 - 1868. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. G. McDonald, M. J. Rieder, M. Nakano, C. K. Hsia, and A. E. Rettie CYP4F2 Is a Vitamin K1 Oxidase: An Explanation for Altered Warfarin Dose in Carriers of the V433M Variant Mol. Pharmacol., June 1, 2009; 75(6): 1337 - 1346. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Li, U. I. Schwarz, M. D. Ritchie, D. M. Roden, C. M. Stein, and D. Kurnik Relative contribution of CYP2C9 and VKORC1 genotypes and early INR response to the prediction of warfarin sensitivity during initiation of therapy Blood, April 23, 2009; 113(17): 3925 - 3930. [Abstract] [Full Text] [PDF] |
||||
![]() |
The International Warfarin Pharmacogenetics Consor Estimation of the Warfarin Dose with Clinical and Pharmacogenetic Data N. Engl. J. Med., February 19, 2009; 360(8): 753 - 764. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. P. Gulseth, G. R. Grice, and W. E. Dager Pharmacogenomics of warfarin: Uncovering a piece of the warfarin mystery Am. J. Health Syst. Pharm., January 15, 2009; 66(2): 123 - 133. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. A. Johnson Ethnic Differences in Cardiovascular Drug Response: Potential Contribution of Pharmacogenetics Circulation, September 23, 2008; 118(13): 1383 - 1393. [Full Text] [PDF] |
||||
![]() |
D. Wang, H. Chen, K. M. Momary, L. H. Cavallari, J. A. Johnson, and W. Sadee Regulatory polymorphism in vitamin K epoxide reductase complex subunit 1 (VKORC1) affects gene expression and warfarin dose requirement Blood, August 15, 2008; 112(4): 1013 - 1021. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. M. Cooper, J. A. Johnson, T. Y. Langaee, H. Feng, I. B. Stanaway, U. I. Schwarz, M. D. Ritchie, C. M. Stein, D. M. Roden, J. D. Smith, et al. A genome-wide scan for common genetic variants with a large influence on warfarin maintenance dose Blood, August 15, 2008; 112(4): 1022 - 1027. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. D. Caldwell, T. Awad, J. A. Johnson, B. F. Gage, M. Falkowski, P. Gardina, J. Hubbard, Y. Turpaz, T. Y. Langaee, C. Eby, et al. CYP4F2 genetic variant alters required warfarin dose Blood, April 15, 2008; 111(8): 4106 - 4112. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. A.R. Doi Initiation of Warfarin in Hospitals Clin. Med. Res., December 1, 2007; 5(4): 207 - 208. [Full Text] [PDF] |
||||
![]() |
H. Skirton Genetics Issues in Nursing Practice Journal of Research in Nursing, November 1, 2007; 12(6): 595 - 596. [PDF] |
||||
![]() |
P. Herijgers and P. Verhamme Improving the quality of anticoagulant therapy in patients with mechanical heart valves: what are we waiting for? Eur. Heart J., October 2, 2007; 28(20): 2424 - 2426. [Full Text] [PDF] |
||||
![]() |
G. Lippi, G. L. Salvagno, and G. C. Guidi Genetic Factors for Warfarin Dose Prediction Clin. Chem., September 1, 2007; 53(9): 1721 - 1722. [Full Text] [PDF] |
||||
![]() |
G. Lippi MD, G. L. Salvagno MD, and G. C. Guidi MD Genetic analysis to prevent warfarin complications Can. Med. Assoc. J., August 14, 2007; 177(4): 377 - 377. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |