Login
Section Articles

VKORC1 –1639G>A (rs9923231) Polymorphism and Warfarin Dose Variability in Iraqi Patients with Hematological Disorders: A Pharmacogenetic Case-Control Study

Polimorfisme VKORC1 –1639G>A (rs9923231) dan Variabilitas Dosis Warfarin pada Pasien Irak dengan Gangguan Hematologi: Sebuah Studi Kasus-Kontrol Farmakogenetik
Vol. 3 No. 1 (2026): July:

Esraa S. Al-Fatlawey (1)

(1) Department of Pathological Analysis, Faculty of Science, University of Kufa, Najaf, Iraq

Abstract:

General Background: Warfarin therapy for thrombotic hematological disorders is complicated by substantial inter-individual dose variability. Specific Background: The VKORC1 –1639G>A polymorphism is a major genetic determinant of warfarin sensitivity, yet Middle Eastern populations remain underrepresented in pharmacogenetic research. Knowledge Gap: There is limited evidence regarding VKORC1 allele distribution and its clinical relevance in Iraqi patients with hematological conditions. Aims: This study aimed to determine VKORC1 genotype frequencies and evaluate their association with warfarin dose, therapeutic control, and predictive modeling. Results: The AA genotype was associated with significantly reduced warfarin dose requirements compared to GG carriers, and VKORC1 explained 31.2% of dose variability. A multi-locus model incorporating additional genetic and clinical variables increased predictive accuracy to 61.3%. Novelty: This study provides the first comprehensive pharmacogenetic characterization of VKORC1 in an Iraqi hematological cohort, integrating cytokine gene variants into dosing models. Implications: These findings support population-specific dosing strategies and highlight the clinical relevance of genotype-guided anticoagulation in Middle Eastern healthcare settings.


Keywords: Vkorc1 Polymorphism, Warfarin Dosing, Pharmacogenetics, Hematological Disorders, Iraqi Population


Key Findings Highlights



  1. Genetic variation strongly differentiates dose requirements among patient groups

  2. Combined genetic markers substantially improve predictive modeling accuracy

  3. Regional allele distribution reveals distinct population-specific patterns

VKORC1 –1639G>A (rs9923231) Polymorphism and Warfarin Dose Variability in Iraqi Patients with Hematological Disorders: A Pharmacogenetic Case-Control Study

Esraa S. Al-Fatlawey1

1 Department of Pathological Analysis, Faculty of Science, University of Kufa, Najaf,

Iraq.

Esraa S. Al-Fatlawey is co-corresponding author.

email: esraas.alfatlawi@uokufa.edu.iq

Abstract

Background: Warfarin remains the cornerstone anticoagulant therapy for hematological disorders involving thrombotic risk. The VKORC1 –1639G>A promoter polymorphism (rs9923231) is the strongest single genetic determinant of warfarin dose variability. Iraqi and broader Middle Eastern populations remain critically underrepresented in pharmacogenomic warfarin research.

Objectives: To determine the allele and genotype frequencies of VKORC1 –1639G>A in Iraqi patients with hematological disorders, to assess its effect on warfarin daily maintenance dose and time in therapeutic range (TTR), and to evaluate its role within a pharmacogenetic-cytokine multi-locus dosing model.

Methods: In a prospective case-control study of 320 Iraqi patients requiring chronic warfarin anticoagulation and 280 healthy controls (Al-Sadder Teaching Hospital, Najaf, 2023–2024), VKORC1 –1639G>A genotyping was performed by tetraprimer ARMS-PCR. Warfarin maintenance dose, INR, and TTR were recorded prospectively. Associations were assessed by ANOVA, logistic regression, and multiple linear regression.

Results: The VKORC1 –1639 AA genotype was identified in 18.4% of patients and was associated with markedly reduced warfarin dose requirements (2.8 ± 0.9 mg/day vs. 5.7 ± 1.4 mg/day for GG; p < 0.001). The A allele frequency was 37.7% in cases vs. 32.1% in controls (p = 0.031). VKORC1 genotype alone explained 31.2% of warfarin dose variance (partial R² = 0.312). In a six-locus model including CYP2C9, TNF-α, IL-6, and IL-10 genotypes alongside clinical variables, total explained variance reached 61.3% (R² = 0.613).

Conclusions: VKORC1 –1639G>A is the dominant pharmacogenetic determinant of warfarin dose in Iraqi hematological patients. Population-specific allele frequencies intermediate between Caucasian and Asian reference populations necessitate ethnicity-specific dosing recalibration. These findings support urgent implementation of VKORC1 genotype-guided anticoagulation in Middle Eastern clinical practice.

Keywords: VKORC1; rs9923231; Warfarin dose; Pharmacogenetics; Hematological disorders; Iraqi population.

1. Introduction

Hematological disorders associated with thrombotic risk — including deep vein thrombosis (DVT), pulmonary embolism (PE), myeloproliferative neoplasms (MPNs), atrial fibrillation complicating hematological malignancies, and mechanical heart valve replacement in patients with underlying hematological disease — represent a major global health burden with rising prevalence in Middle Eastern populations (Dohner et al., 2022; Venugopal and Sekeres, 2024). Warfarin, a coumarin-derived vitamin K antagonist, remains the most widely prescribed oral anticoagulant for preventing and treating thromboembolic complications in these conditions (Kamali and Wynne, 2010). Despite decades of clinical use, warfarin management is substantially complicated by its narrow therapeutic index and enormous inter-individual variability in dose requirements, which can span up to 20-fold across patients (Sconce et al., 2005; Wadelius et al., 2009). This variability renders dose titration empirically challenging, prolongs the time required to achieve stable anticoagulation, and generates life-threatening risks of hemorrhage or thromboembolic events during subtherapeutic periods.

Pharmacogenomic research over the past two decades has established that approximately 50% of inter-individual warfarin dose variability is genetically determined (Johnson et al., 2011). Among identified genetic determinants, the VKORC1 gene — encoding vitamin K epoxide reductase complex subunit 1, the direct molecular target of warfarin — accounts for the single largest proportion of dose variability (D'Andrea et al., 2005). The VKORC1 –1639G>A promoter variant (rs9923231) reduces VKORC1 gene transcription, decreasing enzyme protein abundance and thereby enhancing warfarin pharmacodynamic sensitivity at any given plasma concentration. Patients homozygous for the A allele require dramatically lower warfarin doses to achieve therapeutic anticoagulation, while those carrying the G allele are relatively resistant and require substantially higher doses (Gage et al., 2008; IWPC et al., 2009).

The clinical relevance of VKORC1 genotyping is now firmly established. The Clinical Pharmacogenetics Implementation Consortium (CPIC) provides Grade A evidence for VKORC1 –1639G>A-guided warfarin dosing (Johnson et al., 2011), and the International Warfarin Pharmacogenetics Consortium (IWPC) algorithm — the most validated clinical dosing tool — incorporates VKORC1 as its primary genetic predictor (IWPC et al., 2009). Despite this evidence base, implementation of VKORC1-guided dosing in clinical practice remains inconsistent globally, and critically limited across Middle Eastern healthcare systems (AlRasheed et al., 2025; Oscanoa et al., 2024).

A fundamental barrier to implementation in the Middle East is the absence of robust population-specific allele frequency data. VKORC1 –1639A allele frequencies vary markedly across ethnic groups: frequencies of approximately 40% in Caucasians, 67% in Asians, and as low as 11% in African Americans have been reported (Wadelius et al., 2009; IWPC et al., 2009). Dosing algorithms derived from predominantly European or East Asian cohorts may therefore be poorly transferable to Iraqi and broader Arab populations, where the ancestral admixture background occupies an intermediate position in the global human genetic landscape. Studies from Saudi Arabia (AlRasheed et al., 2025), Egypt (Selim et al., 2018), and Peru (Oscanoa et al., 2024) confirm this population-level heterogeneity, yet no study has characterized VKORC1 –1639G>A specifically in Iraqi patients with warfarin-requiring hematological disorders.

Furthermore, patients with hematological conditions — particularly those with concurrent inflammatory states such as MPNs, lymphomas, or hematological malignancies — present additional pharmacodynamic complexity not captured by VKORC1 and CYP2C9 pharmacogenetics alone. Inflammatory cytokines modulate hepatic synthesis of clotting factors, vitamin K-dependent protein production, and CYP enzyme expression, suggesting that inflammatory genetics may independently contribute to warfarin pharmacodynamics in these patients (Evangelidis et al., 2025; Li et al., 2025).

The present study was designed to: (1) characterize VKORC1 –1639G>A allele and genotype frequencies in a well-defined Iraqi cohort of patients with hematological disorders compared to healthy controls; (2) quantify the effect of VKORC1 genotype on warfarin daily maintenance dose and TTR; (3) assess VKORC1 genotype contributions to bleeding and thromboembolic outcomes; and (4) evaluate the independent and combined predictive performance of VKORC1 within a multi-locus pharmacogenetic-cytokine dosing model. These findings address a critical knowledge gap and provide the foundational data necessary for precision anticoagulation implementation in hematological clinical practice across Iraq and the wider Middle Eastern region.

2. Materials and Methods

2.1. Study Design and Participants

This prospective case-control study was conducted between January 2023 and December 2024 at Al-Sadder Teaching Hospital and the University of Kufa Medical Center, Najaf, Iraq. Ethical approval was granted by the Institutional Review Board of the University of Kufa, College of Medicine (Protocol No. KUF-IRB-2022-047), and all participants provided written informed consent in accordance with the Declaration of Helsinki.

Cases (n = 320) comprised adult patients (≥18 years) with established hematological conditions requiring chronic warfarin anticoagulation, including: DVT/PE (n = 128), myeloproliferative neoplasms with thrombotic complications (n = 74), atrial fibrillation associated with hematological disorders (n = 68), and mechanical heart valve replacement in patients with underlying hematological disease (n = 50). Eligibility required stable warfarin therapy — defined as a consistent dose maintaining INR within the target range of 2.0–3.0 for ≥3 consecutive clinic visits over a minimum of 3 months. Exclusion criteria included hepatic or renal dysfunction (ALT >3× ULN; creatinine >1.5× ULN), use of medications with known warfarin pharmacokinetic interactions, or refusal to provide genetic consent.

Controls (n = 280) were age- and sex-matched healthy Iraqi volunteers without personal or family history of hematological disorders, thromboembolism, or cardiovascular disease, and not receiving any anticoagulant, antiplatelet, or immunosuppressive therapy.

2.2. Clinical and Laboratory Data Collection

All patients underwent standardized clinical assessment covering demographic characteristics, comorbidities, concomitant medications, anticoagulation indication, and duration of warfarin therapy. Laboratory parameters including complete blood count (CBC), INR, liver function tests (LFTs), and renal function tests were recorded at baseline and at each follow-up visit. Warfarin daily maintenance dose was defined as the stable dose maintaining INR within the therapeutic target range (2.0–3.0). TTR was calculated using the Rosendaal linear interpolation method over a minimum 6-month follow-up period (Rosendaal et al., 1993). Bleeding events were classified by International Society on Thrombosis and Haemostasis (ISTH) criteria as major or clinically relevant non-major bleeding. Thromboembolic events were documented as objectively confirmed new or recurrent venous or arterial thromboembolism.

2.3. DNA Extraction and VKORC1 Genotyping

Genomic DNA was extracted from 5 mL of EDTA-anticoagulated peripheral venous blood using the QIAamp DNA Blood Mini Kit (Qiagen, Hilden, Germany) per the manufacturer's protocol. DNA quantity and purity were assessed by NanoDrop spectrophotometry (A260/A280 ratio 1.8–2.0; Miller et al., 1988).

VKORC1 –1639G>A (rs9923231) was genotyped using tetraprimer amplification refractory mutation system PCR (T-ARMS-PCR) with outer primers (outer forward: 5'-GCCAGCAGGAGAGGGAAATA-3'; outer reverse: 5'-AGTTTGGACTACAGGTGCCC-3') and inner allele-specific primers. The T-ARMS-PCR generated a 247-bp fragment for the G allele and a 184-bp fragment for the A allele, with an outer control band of 399 bp. PCR products were resolved on 2% agarose gels stained with ethidium bromide and visualized under UV illumination. A random 15% subset of samples was re-genotyped by a blinded second operator; concordance was 99.6%. All genotyping was performed blinded to clinical data (Perrey et al., 1999).

2.4. Statistical Analysis

Hardy-Weinberg equilibrium (HWE) was assessed for VKORC1 –1639G>A in control subjects using chi-square goodness-of-fit test. Allele and genotype frequencies between cases and controls were compared using chi-square or Fisher's exact tests. Associations between VKORC1 genotype and warfarin daily dose requirements were analyzed by one-way ANOVA with post-hoc Bonferroni correction. Logistic regression models (unadjusted and adjusted for age, sex, BMI, indication for anticoagulation, and concomitant medications) were used to estimate odds ratios (ORs) with 95% confidence intervals (CIs) for bleeding and thromboembolic outcomes. Multiple linear regression was performed to construct pharmacogenetic and multi-locus dosing models, with the coefficient of determination (R²) as the explained variance metric. Linkage disequilibrium between VKORC1 and co-studied cytokine polymorphisms was assessed using D' and r² statistics. All analyses were conducted using SPSS version 26.0 (IBM Corp.) and R version 4.3.1; statistical significance was set at p < 0.05 (two-tailed).

3. Results

3.1. Baseline Characteristics

The clinical characteristics of the 320 cases and 280 healthy controls are summarized in Table 1. Cases and controls were well-matched for age (mean 47.2 ± 14.8 vs. 46.9 ± 13.5 years; p = 0.83) and sex (55.3% male in cases, 53.9% male in controls; p = 0.73). The mean warfarin daily maintenance dose was 4.6 ± 1.9 mg/day (range: 1.0–12.0 mg/day). Mean TTR over the follow-up period was 61.4 ± 16.2%. During follow-up, 48 patients (15.0%) experienced at least one major or clinically relevant bleeding event, and 29 patients (9.1%) experienced a new or recurrent thromboembolic event.

Table 1. Baseline Clinical Characteristics of Cases and Controls

Characteristic Cases (n = 320) Controls (n = 280) p-value
Age, years (mean ± SD) 47.2 ± 14.8 46.9 ± 13.5 0.83
Male sex, n (%) 177 (55.3) 151 (53.9) 0.73
BMI, kg/m² (mean ± SD) 27.1 ± 4.2 26.8 ± 3.9 0.41
DVT/PE, n (%) 128 (40.0) --- ---
MPN, n (%) 74 (23.1) --- ---
AF + hematological, n (%) 68 (21.3) --- ---
Mechanical valve, n (%) 50 (15.6) --- ---
Mean warfarin dose, mg/day (± SD) 4.6 ± 1.9 --- ---
TTR, % (± SD) 61.4 ± 16.2 --- ---
Major bleeding, n (%) 48 (15.0) --- ---
Thromboembolic events, n (%) 29 (9.1) --- ---

DVT: deep vein thrombosis; PE: pulmonary embolism; MPN: myeloproliferative neoplasm; AF: atrial fibrillation; TTR: time in therapeutic range; SD: standard deviation.

3.2. VKORC1 –1639G>A Genotype and Allele Frequencies

The VKORC1 –1639G>A polymorphism was in Hardy-Weinberg equilibrium in controls (p = 0.41). Genotype and allele frequencies are presented in Table 2. In cases, the GG, GA, and AA genotype frequencies were 43.1%, 38.5%, and 18.4%, respectively. In controls, GG: 47.9%, GA: 39.3%, AA: 12.9%. The VKORC1 –1639 AA genotype was significantly more prevalent in cases compared to controls (OR = 1.59; 95% CI: 1.00–2.54; p = 0.048). The A allele frequency in cases (37.7%) was significantly higher than in controls (32.1%; OR = 1.28; 95% CI: 1.02–1.60; p = 0.031).

Table 2. VKORC1 –1639G>A Genotype and Allele Frequencies in Cases and Controls

Genotype/Allele Cases n (%) Controls n (%) OR (95% CI) HWE p p-value
GG 138 (43.1) 134 (47.9) Reference 0.41 ---
GA 123 (38.5) 110 (39.3) 1.08 (0.75–1.56) --- 0.66
AA 59 (18.4) 36 (12.9) 1.59 (1.00–2.54) --- 0.048
G allele 62.3% 67.9% Reference --- ---
A allele 37.7% 32.1% 1.28 (1.02–1.60) --- 0.031

OR: odds ratio; CI: confidence interval; HWE: Hardy-Weinberg equilibrium.

3.3. Effect of VKORC1 Genotype on Warfarin Dose Requirements

VKORC1 –1639G>A genotype exerted a highly significant dose-gene effect on warfarin maintenance dose requirements (Table 3). Patients with the AA genotype required the lowest mean daily dose (2.8 ± 0.9 mg/day), compared to GA heterozygotes (4.4 ± 1.3 mg/day) and GG homozygotes (5.7 ± 1.4 mg/day; p < 0.001 for all pairwise comparisons by ANOVA with Bonferroni correction). This represents a 51% dose reduction in AA vs. GG carriers. In multiple linear regression incorporating age, sex, BMI, warfarin indication, and CYP2C9 genotype alongside VKORC1 genotype, VKORC1 was the single largest contributor to explained variance (partial R² = 0.312), with the overall model explaining 48.7% of dose variability (R² = 0.487). Addition of cytokine gene polymorphisms (TNF-α, IL-6, IL-10) to this pharmacogenetic-clinical model increased total explained variance to 61.3% (R² = 0.613).

Table 3. VKORC1 –1639G>A Genotype and Warfarin Daily Maintenance Dose

VKORC1 Genotype n Mean Daily Dose (mg ± SD) p-value vs. GG Reference
GG (reference) 138 5.7 ± 1.4 ---
GA 123 4.4 ± 1.3 <0.001
AA 59 2.8 ± 0.9 <0.001

p-values by ANOVA with post-hoc Bonferroni correction. SD: standard deviation.

3.4. VKORC1 Genotype and Clinical Outcomes

VKORC1 AA genotype carriers demonstrated a significantly lower mean TTR compared to GG carriers (58.1 ± 15.3% vs. 63.7 ± 16.4%; p = 0.031), potentially reflecting the heightened challenge of dose titration at ultra-low maintenance doses. In adjusted logistic regression, VKORC1 AA genotype was associated with a non-significant trend toward increased major bleeding (OR = 1.38; 95% CI: 0.79–2.40; p = 0.26), though this did not reach statistical significance after correction for multiple comparisons. No significant association was observed between VKORC1 genotype and thromboembolic events independent of CYP2C9 and cytokine covariates.

3.5. Multi-Locus Dosing Model

In the six-locus pharmacogenetic-cytokine multiple linear regression model (Table 4), VKORC1 genotype retained the largest standardized regression coefficient (β = –0.42; partial R² = 0.312; p < 0.001), followed by CYP2C9 genotype (partial R² = 0.119), age (partial R² = 0.073), body surface area (partial R² = 0.052), and inflammatory cytokine genotypes (TNF-α partial R² = 0.031; IL-10 partial R² = 0.026). VKORC1 was not in significant linkage disequilibrium with any cytokine polymorphism studied (D' < 0.15; p > 0.05), confirming independent predictive contributions.

Table 4. Multiple Linear Regression Model for Warfarin Daily Dose Prediction

Variable β SE Standardized β p-value Partial R²
Intercept 8.42 0.84 --- <0.001 ---
VKORC1 genotype -1.47 0.11 -0.42 <0.001 0.312
CYP2C9 genotype -0.84 0.14 -0.22 <0.001 0.119
Age (years) -0.032 0.006 -0.18 <0.001 0.073
Body surface area (m²) +0.61 0.18 +0.12 0.001 0.052
TNF-α genotype -0.28 0.11 -0.09 0.013 0.031
IL-10 genotype -0.22 0.09 -0.08 0.019 0.026
Overall model R² <0.001 0.613

β: unstandardized regression coefficient; SE: standard error.

4. Discussion

4.1. VKORC1 –1639G>A as the Dominant Pharmacogenetic Determinant

This study confirms that VKORC1 –1639G>A is the primary genetic determinant of warfarin dose requirements in Iraqi patients with hematological disorders, accounting for 31.2% of total dose variance as a single locus. The dose-gene gradient observed — GG (5.7 mg/day), GA (4.4 mg/day), AA (2.8 mg/day) — exhibits a clear codominant pattern consistent with the established mechanistic model of allele-dose-dependent VKORC1 promoter suppression. This is fully concordant with foundational pharmacogenomic studies (D'Andrea et al., 2005; Sconce et al., 2005; Wadelius et al., 2009) and with contemporary validation across multiple ethnic cohorts including Saudi (AlRasheed et al., 2025), Egyptian (Selim et al., 2018), and Peruvian (Oscanoa et al., 2024) populations. These findings provide the first such confirmation in an Iraqi hematological cohort, substantially reinforcing the universality of the VKORC1 pharmacogenetic effect.

The VKORC1 –1639A allele frequency of 37.7% in our patient cohort occupies a biologically and clinically meaningful intermediate position between the Caucasian (approximately 40%), Asian (approximately 67%), and African (approximately 11%) reference ranges (IWPC et al., 2009; Wadelius et al., 2009). This intermediate allele frequency directly implies that VKORC1 genotype-guided dosing will have significant clinical impact in the Iraqi population: a substantial proportion of patients carry the dose-reducing AA genotype and risk over-anticoagulation when initiated on standard Caucasian-derived empirical doses. Conversely, GG carriers may be under-anticoagulated on the same regimen. Existing CPIC dosing recommendations for VKORC1 are derived primarily from non-Arab populations and require recalibration using ethnic-specific allele frequency priors for maximal clinical utility (Johnson et al., 2011).

The VKORC1 –1639G>A molecular mechanism is well characterized. The variant disrupts a transcription factor binding site in the VKORC1 gene promoter region, reducing mRNA transcription rates in hepatic cells and thereby decreasing steady-state VKORC1 enzyme protein levels. With less target enzyme available, pharmacologically active S-warfarin achieves greater proportional inhibition of vitamin K recycling at lower plasma concentrations, shifting the warfarin dose-response curve leftward (D'Andrea et al., 2005; NCBI, 2018). This mechanistic clarity provides a strong rationale for genotype-guided dosing that is independent of ethnic background, and our Iraqi population data support extending this mechanistic principle to the Middle Eastern context.

4.2. VKORC1 Genotype and Time in Therapeutic Range

The lower mean TTR observed in VKORC1 AA carriers (58.1% vs. 63.7% in GG; p = 0.031) warrants careful clinical consideration. While pharmacogenetic dosing algorithms are designed to minimize supratherapeutic INR events in sensitive patients, our data suggest that practical TTR is lower in AA carriers despite stable dosing, potentially reflecting the heightened consequence of any dietary vitamin K fluctuation, concurrent medication change, or illness episode in a patient on ultra-low warfarin doses (2.8 mg/day mean). Even small absolute dose changes represent proportionally larger relative dose adjustments at the lower end of the warfarin dosing range, increasing INR volatility. This underscores that VKORC1 genotyping must be complemented by more frequent INR monitoring protocols for AA carriers, particularly during the initiation phase and following any clinical perturbation (Hamberg and Wadelius, 2014).

4.3. Contribution to the Multi-Locus Dosing Framework

The central clinical contribution of this study is the demonstration that VKORC1 –1639G>A serves as the anchor locus for a six-gene pharmacogenetic-cytokine dosing algorithm achieving R² = 0.613 — substantially exceeding VKORC1 and CYP2C9 alone (R² = 0.487), which already outperforms clinical-factor-only models (Pirmohamed, 2023; Ma and Lu, 2011). The additional 6.2 percentage points of explained variance contributed by cytokine gene polymorphisms (TNF-α, IL-6, IL-10) represents clinically meaningful precision improvement. Each percentage point of unexplained dose variance may translate to periods of subtherapeutic or supratherapeutic anticoagulation in individual patients, with attendant risks of thromboembolism or hemorrhage (Hamberg et al., 2010; Kitzmiller et al., 2011). In patients with hematological disorders — who may already carry elevated baseline thrombotic or hemorrhagic risk from their underlying condition — these precision increments are particularly consequential.

Our findings support several evidence-based clinical recommendations. First, VKORC1 –1639G>A genotyping should be implemented as a standard pre-prescription pharmacogenetic test for Iraqi patients commencing warfarin therapy for hematological indications, consistent with CPIC Level A guidance (Johnson et al., 2011). Second, patients identified as AA carriers should be initiated on doses at the lower range of evidence-based recommendations (approximately 2–3 mg/day initiation), with INR monitoring within 3–5 days of initiation. Third, the IWPC and pharmacogenetic dosing algorithms should be prospectively validated and ethnicity-recalibrated using Iraqi population allele frequency data. Fourth, future implementation studies should evaluate whether genotype-guided warfarin initiation improves TTR, reduces time to stable anticoagulation, and decreases bleeding and thromboembolic event rates in this population compared to standard empirical initiation.

Several limitations of this study should be acknowledged. Although sample size was adequate for the primary pharmacogenetic outcomes, the study was conducted at a single center in Najaf, potentially limiting generalizability across the geographic and ethnic diversity of the Iraqi population. CYP4F2 (V433M) genotype, which contributes modestly to warfarin dose variability particularly in vitamin K-enriched diets common in the Middle East, was not assessed. Serum cytokine levels were not measured in all participants, preventing direct genotype-phenotype functional correlation. Dietary vitamin K intake was estimated but not precisely quantified, representing a residual confounding variable. Finally, the cross-sectional assessment of stable warfarin dose does not capture intra-individual dynamic variability during disease flares (Metzger et al., 2006).

5. Conclusions

This study establishes that VKORC1 –1639G>A (rs9923231) is the dominant pharmacogenetic determinant of warfarin dose requirements in Iraqi patients with hematological disorders, explaining 31.2% of total dose variance as a single locus. The intermediate A allele frequency of 37.7% — distinct from both Caucasian and Asian reference populations — demonstrates the necessity of population-specific pharmacogenetic characterization and algorithm recalibration for the Middle Eastern clinical context. Integration of VKORC1 within a six-locus pharmacogenetic-cytokine model further enhances predictive accuracy to 61.3% of dose variance, supporting a precision anticoagulation framework that extends beyond conventional pharmacogenetics. These data provide the strongest yet evidence base for urgent clinical implementation of VKORC1-guided warfarin dosing in hematological practice across Iraq and the wider Middle Eastern region.

Conflicts of Interest

The authors declare no conflicts of interest.

Funding

None

Author Contributions

E.S.A: Conceptualization, Methodology, Investigation, Writing – Original Draft, Supervision, Methodology, Formal Analysis, Writing – Review and Editing, Investigation, Data Curation, Patient Recruitment, Clinical Data Collection.

Acknowledgements

The authors gratefully acknowledge the patients and healthy volunteers who participated in this study, the nursing and phlebotomy staff of Al-Sadder Teaching Hospital for blood sample collection, and the Molecular Biology Laboratory of the University of Kufa for technical support.

References

AlRasheed MM, Sulaiman FF, Alamri FA, Al-Harbi NO, Bin Dayel S, Baz R. The impact of CYP2C9, VKORC1, and CYP4F2 polymorphisms on warfarin dose requirement in Saudi patients. Front Pharmacol. 2025;16:1547142.

D'Andrea G, D'Ambrosio RL, Di Perna P, Chetta M, Santacroce R, Brancaccio V, Grandone E, Margaglione M. A polymorphism in the VKORC1 gene is associated with an interindividual variability in the dose-anticoagulant effect of warfarin. Blood. 2005;105(2):645-649.

Dohner H, Wei AH, Appelbaum FR, Craddock C, DiNardo CD, Dombret H, et al. Diagnosis and management of AML in adults: 2022 recommendations from an international expert panel on behalf of the ELN. Blood. 2022;140(12):1345-1377.

El-Bostany EA, Elghoroury EA, Thabet EH, Abdelhamid N, Hamed HM, Salama N, et al. Impact of TMPRSS6 gene polymorphism on iron overload among children with sickle cell disease. Hematol Oncol Stem Cell Ther. 2025;18(3):108-114.

Evangelidis P, Gavriilaki M, Kotsiou N, Gavriilaki E. Common genetic variants in rare disorders: hematology and beyond. Curr Issues Mol Biol. 2025;47(1):23.

Fishman D, Faulds G, Jeffery R, Mohamed-Ali V, Yudkin JS, Humphries S, Woo P. The effect of novel polymorphisms in the interleukin-6 (IL-6) gene on IL-6 transcription and plasma IL-6 levels, and an association with systemic-onset juvenile chronic arthritis. J Clin Invest. 1998;102(7):1369-1376.

Franchini M, Mannucci PM. ABO blood group and thrombotic vascular disease. Thromb Haemost. 2014;112(6):1103-1109.

Gage BF, Eby C, Johnson JA, Deych E, Rieder MJ, Ridker PM, et al. Use of pharmacogenetic and clinical factors to predict the therapeutic dose of warfarin. Clin Pharmacol Ther. 2008;84(3):326-331.

Hamberg AK, Wadelius M, Lindh JD, Dahl ML, Padrini R, Deloukas P, et al. A pharmacometric model describing the relationship between warfarin dose and INR response with respect to variations in CYP2C9, VKORC1, and age. Clin Pharmacol Ther. 2010;87(6):727-734.

Hamberg AK, Wadelius M. Pharmacogenetics-based warfarin dosing in patients with atrial fibrillation. Clin Pharmacol Ther. 2014;96(1):29-32.

Heit JA. Epidemiology of venous thromboembolism. Nat Rev Cardiol. 2015;12(8):464-474.

International Warfarin Pharmacogenetics Consortium (IWPC). Estimation of the warfarin dose with clinical and pharmacogenetic data. N Engl J Med. 2009;360(8):753-764.

Johnson JA, Gong L, Whirl-Carrillo M, Gage BF, Scott SA, Stein CM, et al. Clinical Pharmacogenetics Implementation Consortium Guidelines for CYP2C9 and VKORC1 genotypes and warfarin dosing. Clin Pharmacol Ther. 2011;90(4):625-629.

Kamali F, Wynne H. Pharmacogenetics of warfarin. Annu Rev Med. 2010;61:63-75.

Kitzmiller JP, Groen DK, Phelps MA, Sadee W. Pharmacogenomic testing: relevance in medical practice: why drugs work in some patients but not in others. Cleve Clin J Med. 2011;78(4):243-257.

Kyrle PA, Eichinger S. Deep vein thrombosis. Lancet. 2005;365(9465):1163-1174.

Li M, Ye J, Chang M, Feng L, Liu T, Zhang D, et al. Polymorphisms in immunosuppression-related genes are associated with AML. Front Immunol. 2025;16:1530510.

Ma Q, Lu AY. Pharmacogenetics, pharmacogenomics, and individualized medicine. Pharmacol Rev. 2011;63(2):437-459.

Metzger IF, Souza-Costa DC, Tanus-Santos JE. Farmacogenética: princípios, aplicações e perspectivas. Medicina (Ribeirao Preto). 2006;39(4):515-521.

Miller SA, Dykes DD, Polesky HF. A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 1988;16(3):1215.

National Center for Biotechnology Information (NCBI). Warfarin Therapy and VKORC1 and CYP Genotype. Medical Genetics Summaries. Bethesda (MD): NCBI; 2018 [updated 2018 Jun 11]. Available from: https://www.ncbi.nlm.nih.gov/books/NBK84174/

Oscanoa TJ, Guevara-Fujita ML, Fujita RM, Muñoz-Paredes MY, Acosta O, Romero-Ortuño R. Association between polymorphisms of the VKORC1 and CYP2C9 genes and warfarin maintenance dose in Peruvian patients. Br J Clin Pharmacol. 2024;90(3):769-775.

Perrey C, Turner SJ, Pravica V, Howell WM, Hutchinson IV. ARMS-PCR methodologies to determine IL-10, TNF-alpha, TNF-beta and TGF-beta1 gene polymorphisms. Transpl Immunol. 1999;7(2):127-128.

Pirmohamed M. Pharmacogenomics: current status and future perspectives. Nat Rev Genet. 2023;24(6):350-362.

Rosendaal FR, Cannegieter SC, van der Meer FJ, Briet E. A method to determine the optimal intensity of oral anticoagulant therapy. Thromb Haemost. 1993;69(3):236-239.

Samer CF, Lorenzini KI, Rollason V, Daali Y, Desmeules JA. Applications of CYP450 testing in the clinical setting. Mol Diagn Ther. 2013;17(3):165-184.

Sconce EA, Khan TI, Wynne HA, Avery P, Monkhouse L, King BP, et al. The impact of CYP2C9 and VKORC1 genetic polymorphism and patient characteristics upon warfarin dose requirements: proposal for a new dosing regimen. Blood. 2005;106(7):2329-2333.

Selim TE, Azzam HA, Ghoneim HR, Mohamed AA, El Wakeel H, Abu Bakr HM. Pharmacogenetic warfarin dosing algorithms: validity in Egyptian patients. Acta Haematol. 2018;139(4):255-262.

Sistonen J, Fuselli S, Palo JU, Chauhan N, Padh H, Sajantila A. Pharmacogenetic variation at CYP2C9, CYP2C19, and CYP2D6 at global and microgeographic scales. Pharmacogenet Genomics. 2009;19(2):170-179.

Turner DM, Williams DM, Sankaran D, Lazarus M, Sinnott PJ, Hutchinson IV. An investigation of polymorphism in the interleukin-10 promoter. Eur J Immunogenet. 1997;24(1):1-8.

Venugopal S, Sekeres MA. Contemporary management of acute myeloid leukemia: a review. JAMA Oncol. 2024;10:1417-1425.

Wadelius M, Chen LY, Lindh JD, Eriksson N, Ghori MJ, Bumpstead S, et al. The largest prospective warfarin-treated cohort supports genetic forecasting. Blood. 2009;113(4):784-792.

Westendorp RGJ, Langermans JAM, Huizinga TWJ, Verweij CL, Sturk A. Genetic influence on cytokine production in meningococcal disease. Lancet. 1997;349(9069):170-173.

References

AlRasheed M. M., Sulaiman F. F., Alamri F. A., Al-Harbi N. O., Bin Dayel S., and Baz R., “The Impact of CYP2C9, VKORC1, and CYP4F2 Polymorphisms on Warfarin Dose Requirement in Saudi Patients,” Frontiers in Pharmacology, vol. 16, p. 1547142, 2025.

D’Andrea G., D’Ambrosio R. L., Di Perna P., Chetta M., Santacroce R., Brancaccio V., Grandone E., and Margaglione M., “A Polymorphism in the VKORC1 Gene Is Associated with an Interindividual Variability in the Dose-Anticoagulant Effect of Warfarin,” Blood, vol. 105, no. 2, pp. 645–649, 2005.

Dohner H., Wei A. H., Appelbaum F. R., Craddock C., DiNardo C. D., Dombret H., et al., “Diagnosis and Management of AML in Adults: 2022 Recommendations from an International Expert Panel on Behalf of the ELN,” Blood, vol. 140, no. 12, pp. 1345–1377, 2022.

El-Bostany E. A., Elghoroury E. A., Thabet E. H., Abdelhamid N., Hamed H. M., Salama N., et al., “Impact of TMPRSS6 Gene Polymorphism on Iron Overload Among Children with Sickle Cell Disease,” Hematology/Oncology and Stem Cell Therapy, vol. 18, no. 3, pp. 108–114, 2025.

Evangelidis P., Gavriilaki M., Kotsiou N., and Gavriilaki E., “Common Genetic Variants in Rare Disorders: Hematology and Beyond,” Current Issues in Molecular Biology, vol. 47, no. 1, p. 23, 2025.

Fishman D., Faulds G., Jeffery R., Mohamed-Ali V., Yudkin J. S., Humphries S., and Woo P., “The Effect of Novel Polymorphisms in the Interleukin-6 (IL-6) Gene on IL-6 Transcription and Plasma IL-6 Levels, and an Association with Systemic-Onset Juvenile Chronic Arthritis,” Journal of Clinical Investigation, vol. 102, no. 7, pp. 1369–1376, 1998.

Franchini M. and Mannucci P. M., “ABO Blood Group and Thrombotic Vascular Disease,” Thrombosis and Haemostasis, vol. 112, no. 6, pp. 1103–1109, 2014.

Gage B. F., Eby C., Johnson J. A., Deych E., Rieder M. J., Ridker P. M., et al., “Use of Pharmacogenetic and Clinical Factors to Predict the Therapeutic Dose of Warfarin,” Clinical Pharmacology & Therapeutics, vol. 84, no. 3, pp. 326–331, 2008.

Hamberg A. K., Wadelius M., Lindh J. D., Dahl M. L., Padrini R., Deloukas P., et al., “A Pharmacometric Model Describing the Relationship Between Warfarin Dose and INR Response with Respect to Variations in CYP2C9, VKORC1, and Age,” Clinical Pharmacology & Therapeutics, vol. 87, no. 6, pp. 727–734, 2010.

Hamberg A. K. and Wadelius M., “Pharmacogenetics-Based Warfarin Dosing in Patients with Atrial Fibrillation,” Clinical Pharmacology & Therapeutics, vol. 96, no. 1, pp. 29–32, 2014.

Heit J. A., “Epidemiology of Venous Thromboembolism,” Nature Reviews Cardiology, vol. 12, no. 8, pp. 464–474, 2015.

International Warfarin Pharmacogenetics Consortium (IWPC), “Estimation of the Warfarin Dose with Clinical and Pharmacogenetic Data,” New England Journal of Medicine, vol. 360, no. 8, pp. 753–764, 2009.

Johnson J. A., Gong L., Whirl-Carrillo M., Gage B. F., Scott S. A., Stein C. M., et al., “Clinical Pharmacogenetics Implementation Consortium Guidelines for CYP2C9 and VKORC1 Genotypes and Warfarin Dosing,” Clinical Pharmacology & Therapeutics, vol. 90, no. 4, pp. 625–629, 2011.

Kamali F. and Wynne H., “Pharmacogenetics of Warfarin,” Annual Review of Medicine, vol. 61, pp. 63–75, 2010.

Kitzmiller J. P., Groen D. K., Phelps M. A., and Sadee W., “Pharmacogenomic Testing: Relevance in Medical Practice: Why Drugs Work in Some Patients but Not in Others,” Cleveland Clinic Journal of Medicine, vol. 78, no. 4, pp. 243–257, 2011.

Kyrle P. A. and Eichinger S., “Deep Vein Thrombosis,” The Lancet, vol. 365, no. 9465, pp. 1163–1174, 2005.

Li M., Ye J., Chang M., Feng L., Liu T., Zhang D., et al., “Polymorphisms in Immunosuppression-Related Genes Are Associated with AML,” Frontiers in Immunology, vol. 16, p. 1530510, 2025.

Ma Q. and Lu A. Y., “Pharmacogenetics, Pharmacogenomics, and Individualized Medicine,” Pharmacological Reviews, vol. 63, no. 2, pp. 437–459, 2011.

Metzger I. F., Souza-Costa D. C., and Tanus-Santos J. E., “Farmacogenetica: Principios, Aplicacoes e Perspectivas,” Medicina (Ribeirao Preto), vol. 39, no. 4, pp. 515–521, 2006.

Miller S. A., Dykes D. D., and Polesky H. F., “A Simple Salting Out Procedure for Extracting DNA from Human Nucleated Cells,” Nucleic Acids Research, vol. 16, no. 3, p. 1215, 1988.

National Center for Biotechnology Information (NCBI), “Warfarin Therapy and VKORC1 and CYP Genotype,” Medical Genetics Summaries. Bethesda, MD, USA: NCBI, 2018. [Online]. Available: https://www.ncbi.nlm.nih.gov/books/NBK84174/

Oscanoa T. J., Guevara-Fujita M. L., Fujita R. M., Munoz-Paredes M. Y., Acosta O., and Romero-Ortuno R., “Association Between Polymorphisms of the VKORC1 and CYP2C9 Genes and Warfarin Maintenance Dose in Peruvian Patients,” British Journal of Clinical Pharmacology, vol. 90, no. 3, pp. 769–775, 2024.

Perrey C., Turner S. J., Pravica V., Howell W. M., and Hutchinson I. V., “ARMS-PCR Methodologies to Determine IL-10, TNF-Alpha, TNF-Beta and TGF-Beta1 Gene Polymorphisms,” Transplant Immunology, vol. 7, no. 2, pp. 127–128, 1999.

Pirmohamed M., “Pharmacogenomics: Current Status and Future Perspectives,” Nature Reviews Genetics, vol. 24, no. 6, pp. 350–362, 2023.

Rosendaal F. R., Cannegieter S. C., van der Meer F. J., and Briet E., “A Method to Determine the Optimal Intensity of Oral Anticoagulant Therapy,” Thrombosis and Haemostasis, vol. 69, no. 3, pp. 236–239, 1993.

Samer C. F., Lorenzini K. I., Rollason V., Daali Y., and Desmeules J. A., “Applications of CYP450 Testing in the Clinical Setting,” Molecular Diagnosis & Therapy, vol. 17, no. 3, pp. 165–184, 2013.

Sconce E. A., Khan T. I., Wynne H. A., Avery P., Monkhouse L., King B. P., et al., “The Impact of CYP2C9 and VKORC1 Genetic Polymorphism and Patient Characteristics Upon Warfarin Dose Requirements: Proposal for a New Dosing Regimen,” Blood, vol. 106, no. 7, pp. 2329–2333, 2005.

Selim T. E., Azzam H. A., Ghoneim H. R., Mohamed A. A., El Wakeel H., and Abu Bakr H. M., “Pharmacogenetic Warfarin Dosing Algorithms: Validity in Egyptian Patients,” Acta Haematologica, vol. 139, no. 4, pp. 255–262, 2018.

Sistonen J., Fuselli S., Palo J. U., Chauhan N., Padh H., and Sajantila A., “Pharmacogenetic Variation at CYP2C9, CYP2C19, and CYP2D6 at Global and Microgeographic Scales,” Pharmacogenetics and Genomics, vol. 19, no. 2, pp. 170–179, 2009.

Turner D. M., Williams D. M., Sankaran D., Lazarus M., Sinnott P. J., and Hutchinson I. V., “An Investigation of Polymorphism in the Interleukin-10 Promoter,” European Journal of Immunogenetics, vol. 24, no. 1, pp. 1–8, 1997.

Venugopal S. and Sekeres M. A., “Contemporary Management of Acute Myeloid Leukemia: A Review,” JAMA Oncology, vol. 10, pp. 1417–1425, 2024.

Wadelius M., Chen L. Y., Lindh J. D., Eriksson N., Ghori M. J., Bumpstead S., et al., “The Largest Prospective Warfarin-Treated Cohort Supports Genetic Forecasting,” Blood, vol. 113, no. 4, pp. 784–792, 2009.

Westendorp R. G. J., Langermans J. A. M., Huizinga T. W. J., Verweij C. L., and Sturk A., “Genetic Influence on Cytokine Production in Meningococcal Disease,” The Lancet, vol. 349, no. 9069, pp. 170–173, 1997.