Abstract

Background. Vitamin K epoxide reductase complex subunit 1 (VKORC1) gene encodes a key enzyme with multiple cellular activities, namely, the reduction of vitamin K to its active form. VKORC1-1639G>A (rs9923231) is a common single nucleotide polymorphism with a crucial impact on warfarin dosing and possibly other physiological functions. This study aimed at investigating the frequencies of VKORC1-1639G>A alleles and genotypes in Syrian healthy subjects and patients on warfarin for different indications. Methods. A total of 138 individuals were enrolled in this cross-sectional study. Genomic DNA was extracted from both patients on warfarin and healthy subjects, and polymerase chain reaction (PCR) specific amplicons were genotyped via standard sequencing which also allowed the detection of rs397509427. Comparisons of -1639G>A frequency with other populations were drawn. Results. Of 94 patients on warfarin, 53 (56.38%) were with idiopathic venous thromboembolism (VTE). Despite comparable frequencies of the -1639A allele (47% and 50%), the AA and GA genotypes were at disparate frequencies of 93.2% versus 79.8% in the healthy subjects () versus patients on warfarin, respectively. Carriers of the GG genotype were at a four-fold increased risk of VTE in comparison with those of the AA and GA genotypes (odds ratio , 95% , ). All study subjects were wild-type for the rs397509427 variant. Conclusions. Our results prove a high -1639A prevalence in Syrian healthy subjects and patients on warfarin at frequencies comparable to other Mediterranean and Middle Eastern populations. The A allele carriers are at a lower VTE risk, whereas a global prevalence gradient of the G allele is suggested to be associated with VTE risk.

1. Background

Vitamin K-dependent proteins (VKDPs) have crucial roles in several physiological functions, including blood coagulation, arterial calcification prevention, prenatal heart development, calcium regulation, bone homeostasis, and multiple cellular activities such as growth and inflammatory response [1, 2]. A posttranslational modification step, mediated by gamma-glutamyl carboxylase (GGCX) enzyme that utilizes vitamin K hydroquinone (KH2) as a cofactor, must take place in order to activate VKDPs [1]. Concomitantly, the KH2 is converted to its oxidized form vitamin K 2,3-epoxide (KO), which induces its recycling to the reduced form to compensate for the rapid consumption of KH2 [3]. Vitamin K epoxide reductase (VKOR) controls the rate of the carboxylation reaction by exclusively reducing KO to vitamin K and also assisting the next step of vitamin K to KH2 conversion in this vitamin K cycle [4].

Vitamin K epoxide reductase complex subunit 1 (VKORC1) gene that encodes the VKOR enzyme is located on chromosome 16 [5]. Of the several single nucleotide polymorphisms (SNPs) identified in the VKORC1 gene, “-1639G>A” is the most prominent due to its impact on dosing of warfarin, which exerts its effect by inhibiting VKOR, and subsequently impeding the activation of the coagulation factors II, VII, IX, X [46]. The -1639G>A “rs9923231” is located in the promoter region of VKORC1 gene and the presence of the A allele reduces gene expression leading to fewer mRNA transcripts and eventually fewer copies of the enzyme [7]. Thus, the carriers of the AA genotype must have their therapeutic doses reduced by 70%, according to the drug label released in 2007 by the Food and Drug Administration (FDA) [6]. Owing to the VKORC1-1639G>A significance, previous studies attempted to investigate a plausible association between it and other diseases such as stroke, venous thromboembolism (VTE), cardiovascular diseases, and osteoporosis [810].

A plethora of studies on the frequency of VKORC1-1639G>A among Caucasians and other ethnicities (such as Chinese and Africans) has created a wealth of information that aided in the generation of warfarin dosing algorithms [11]; however, studies on the prevalence of the VKORC1-1639G>A allele in Middle Eastern populations, including Syrians, are relatively scarce. Although Syrians’ genetic makeup is predicted to be similar to that of other populations in the Middle East and Mediterranean basin, nonetheless, the population of Syria exhibits an extremely diverse ethnic composition, consisting of a majority of Arab descents, as well as minorities of Kurds, Assyrians, Armenians, Turkmens, and others. Moreover, Syria has always been a mosaic of diverse faiths and “ethnoreligious” groups; although approximately 90% of all Syrians are Muslim and almost all the rest are Christians, Syrians from both religions are divided up into many sects. Therefore, it is vital to determine the genetic composition of such an ethnically, religiously, and culturally diverse population, as a prerequisite for the implementation of pharmacogenetic-guided dosing approach. Hence, this study aimed at investigating the frequency of VKORC1-1639G>A in Syrian healthy volunteers as well as patients on warfarin for different indications, including idiopathic VTE. Comparisons of -1639G>A frequency with other populations were drawn, and a proposed hypothesis of association between the G allele prevalence with increased VTE risk globally was investigated.

2. Methods

2.1. Study Subjects

A total of 138 unrelated Syrian individuals were enrolled in the study from three major governorates (Damascus, Damascus Countryside, and Homs) from October 2018 to September 2020. Ninety-four outpatients treated with warfarin for different indications, including 53 (56.38%) with idiopathic VTE, were recruited from several hospitals and medical laboratories. The exact number of patients for every indication is shown in Supplementary File 1. Of the 94 patients, 54 were males and 40 females, with a deviation (SD) of . Forty-four healthy unrelated subjects matched with the patients in terms of age and gender (20 males and 24 females, with a of ) were enrolled for the purpose of allele and genotype frequency comparison. Informed consents were obtained from all participants and the legal guardian of one 17-year-old male patient. This cross-sectional study was approved by the Scientific Research Bioethics Committee at the Faculty of Pharmacy, Damascus University; and all methods were performed in accordance with its regulation.

2.2. Genotyping

Genomic DNA was extracted from peripheral blood samples using a solution-based DNA isolation kit (Wizard® Genomic DNA Purification Kit, Promega®, USA) according to the manufacturer’s protocol. -1639G>A genotyping was performed via standard sequencing of specific polymerase chain reaction (PCR) products using specific primers. The amplicon’s length allowed for the detection of another SNP (rs397509427), which is recognized by the insertion of cytosine nucleotide [12]. The reactions contained 30-100 ng of gDNA, 20 pmol of each primer, and 20 μL of master mix (OnePCR™, GeneDirex®, Taiwan); thereafter, the volume was completed with nuclease-free water to 40 μL.

A pair of primers previously published in relevant literature [13, 14] for the amplification of VKORC1-1639G>A containing region was tried at first. Then, we designed another pair using Geneious Pro 4.8.4 (Biomatters®, Auckland, New Zealand). The specificity of our de novo designed forward primer (5-GCCAGCAGGAGAGGGAAATATCA-3) and reverse primer (5-TGGTGTCACCAAGACGCTAGAC-3) was confirmed using BLAST (NCBI, USA) and MFEPrimer 3.1(iGeneTech®, China) and synthesized at Macrogen® Inc. (Seoul, South Korea). The cycling profile is summarized in Table 1.

All reactions yielded the expected product that was confirmed by the appearance of a single band of 457 base pair (pb) length via ethidium bromide-stained 1.5% agarose gel electrophoresis. Amplicons sequencing was conducted at Macrogen®, using the ABI PRISM® BigDye™ Terminator Cycle Sequencing Kit with an ABI PRISM® 3730XI DNA Analyzer, after an enzymatic clean-up method using EnzSAP™ PCR Clean-up Reagent. A representative gel image with representative sequencing chromatograms are shown in Supplementary File 2.

2.3. Statistical Analysis

Statistical analysis was performed using SPSS® (version 25) and GraphPad Prism® (Version 8) software. The observed frequencies of VKORC1-1639G>A genotypes were compared to the expected frequencies from the Hardy–Weinberg equilibrium (HWE) using the chi-square goodness of fit test. The difference between patients’ and healthy subjects’ genotypes was assessed using the chi-square test, whereas Fisher’s exact test was applied to compare the frequencies of the A carriers between VTE patients and healthy subjects and estimate the odds ratio with its corresponding 95% confidence intervals (CI). A value less than 0.05 was considered to be statistically significant in all performed tests.

3. Results

Forty-four healthy subjects and 97 patients were enrolled; however, three patients were nonassessable due to inadequate blood samples () and unsuccessful genotyping (). Thus, the final number of assessable patients was 94. We screened for both rs9923231 (-1639G>A) and rs397509427; however, only the -1639G>A was identified in the study subjects. The heterozygous genotype (GA) constituted the highest frequency genotype in both study groups, and the A allele carriers (GA+AA) represented 79.8% and 93.2% of the patients and healthy subjects, respectively (). Genotypes and allele frequencies for study groups are presented in Table 2. Both groups exhibited a significant deviation from HWE with values of 0.002 and less than 0.001 for the patients and healthy subjects, respectively. Expected values from HWE are shown in Supplementary File 3.

3.1. Association between the -1639G>A SNP and VTE Risk

There was a significant difference between the healthy and patients’ groups in terms of the genotype distribution (). Therefore, a comparison was drawn between the healthy subjects and the patients diagnosed with idiopathic (unprovoked) VTE (). Carriers of GG genotype had shown a 4-fold increased risk of VTE compared to the A allele carriers (GA+AA) (odds ratio , 95% , ) as summarized in Table 3. Nevertheless, no statistically significant difference () was noticed in the GG genotype frequencies (6.8% and 16.2%, respectively) in healthy subjects versus patients () on warfarin for other than VTE indications (i.e., heart valve replacement, ischemic stroke, and atrial fibrillation).

4. Discussion

-1639G>A SNP in VKORC1 gene has been linked to a wide variety of diseases in addition to their crucial role in warfarin dosing, in different populations and ethnicities [2]. This study is the first, to the best of our knowledge, that investigated the frequency of VKORC1-1639G>A in Syrians, both patients on warfarin and healthy volunteers. Our findings demonstrate a high prevalence of the -1639A variant allele in the study populations (47% and 50% in patients and healthy Syrians, respectively). These frequencies were comparable to the frequencies in neighboring populations and the Mediterranean basin populations, whilst greater differences were evident with those of South-African and East-Asian populations [15]. The frequencies of -1639A allele in various populations and countries are summarized in Table 4. The A allele seems to follow a frequency gradient, in which the highest intensity is seen in East Asian populations e.g. Japanese and Chinese, where it is thought to be originated. It starts to decrease towards the southwest of the old-world map, where it becomes relatively low in the south of the African continent (Figure 1).

On the other hand, North and South Americas have proved a discrepancy in the A allele frequency according to the ethnicities of the individuals: Caucasians, Latin Americans, Asians, or Africans [15]. Syria and the Levantine countries in general have both alleles nearly equally prevalent, which could be explained by its diverse genetic makeup of the populations. Its pivotal location at the heart of the old world, between the eastern Mediterranean coast and the Euphrates Valley, has placed Syria at the crossroads of the Eurasian silk and spice trade routes, active from the second century BCE until the mid-15th century. Additionally, wars and multihuman migration waves throughout the history of the Levant have resulted in remarkable ethnic, cultural, and genetic diversity [31].

The heterozygous genotype GA reached higher percentages than expected, based on HWE, in both study groups. This deviation from HWE is less likely to be attributed to genotyping error as standard sequencing was used for the genotype assignment. However, a possible explanation of this deviance might stem from the recent Syrian war that caused dislocations and displacement of millions of Syrians who have fled the country throughout the ten-year War. It is speculated that these tragic circumstances have caused a demographic skewness in the Syrian population causing a subtle kind of stratification [3234].

Another probable reason for such a deviation from HWE stems from the cultural and religious norms, including the high rate of consanguineous marriages, in this Middle Eastern population. Nevertheless, all enrolled subjects in this study were unrelated. A third likable justification is the small sample size since one of the basic assumptions of HWE is an infinitely large population where genetic drift can be ignored. Therefore, a large sample size is preferred in order to exclude selection bias and make the sample well-representative to the community, although this is less likely since the size of our study population is comparable to those in several other frequency studies [17, 23, 33].

Our patients and healthy subjects had similar allele frequencies. To the contrary, a significantly different genotype distribution between healthy subjects and patients with VTE was evident, suggesting an association between the GG genotype with thrombophilia. VTE is a multifactorial disease with a complicated interaction between genetic and environmental factors, and it is estimated that half of the cases occurred without identifiable risk factors [35]. In our study, only unprovoked VTE patients were compared with the healthy subject to negate any other evident causes of the disease such as surgery or pregnancy. The carriers of the A allele comprising genotypes (GA+AA) had a lower risk (75%) for developing VTE than the GG genotype carriers which could be ascribed to the decrease in gene expression of the VKORC1, leading to a reduced protein synthesis (by 44% in homozygotes). The decline in VKORC1 enzyme might affect both the activation status and plasma levels of vitamin K-dependent coagulation factors [9]. Carrying the A allele can be considered as a lifelong administration of anticoagulant at a low dose to prevent or alleviate thrombotic accidents.

A plethora of genetic mutations and polymorphisms are known to be risk factors (i.e., deficiencies of proteins C, S, and antithrombin III; factor V Leiden (FVL); prothrombin (PT) 20210A mutation; and non-O blood groups) for developing VTE [36]. Nevertheless, interindividual variability could be partially explained by these factors; as only 30% of unprovoked VTE cases have been linked to a well-characterized genetic risk factor. Furthermore, the interethnic variations in VTE incidence manifest as a greater dilemma. Although the highest incidence rates of thrombophilia or hypercoagulability (i.e., VTE) have been reported in Africans and African descendants followed by Caucasians and Asians, respectively, the frequencies of the aforementioned familiar genetic risk factors do not correlate with such between-ethnicity disparity in VTE prevalence. For instance, FVL and PT mutations are well-known prothrombotic factors, with FVL being the most common hereditary thrombophilia; however, FVL and PT mutations have the least frequencies in Africa. Furthermore, the frequencies of protein C & S deficiencies in Africa and Europe are comparable. Remarkably, the frequencies of the VKORC1-1639G risk allele are reasonably correlated with VTE incidence around the world, where its highest frequency recorded in Africa (88%) followed by Europe (58%) and the least in Asia (10%) [3739]. The main limitation of this study is the sample size. Therefore, further prospective studies on larger number of subjects are called for to prove or refute the causality association between VKORC1-1639GG genotype and VTE incidence.

5. Conclusions

This study proves a high prevalence of the VKORC1-1639A allele in Syrians, with almost equal frequencies in both patients on warfarin (50%) and healthy subjects (47%). Our findings shed light on the frequency gradient of this allele across the globe. It also suggests a possible protective effect for the -1639A allele against VTE. Moreover, our study may help reveal another piece of the complicated VTE genetic predisposition puzzle.

Abbreviations

bp:Base pairs
CI:Confidence intervals
FVL:Factor V Leiden
FDA:Food and Drug Administration
GGCX:Gamma-glutamyl carboxylase
HWE:Hardy–Weinberg equilibrium
KO:Vitamin K 2,3-epoxide
KH2:Vitamin K hydroquinone
OR:Odds ratio
PCR:Polymerase chain reaction
PT:Prothrombin
SD:Standard deviation
SNP:Single nucleotide polymorphism
VKDPs:Vitamin K-dependent proteins
VKOR:Vitamin K epoxide reductase
VKORC1:Vitamin K epoxide reductase complex subunit 1
VTE:Venous thromboembolism.

Data Availability

The datasets generated and analysed during the current study are available in the ClinVar repository, with the accession number SCV003836537.

Ethical Approval

The study protocol was approved by the Scientific Research Ethics Committee at the Faculty of Pharmacy, Damascus University (Damascus, Number: 6; Date: July 26, 2018) in accordance with the Declaration of Helsinki (1964). All the study subjects gave signed informed consent.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this article.

Authors’ Contributions

LAY conceived and designed the study, and was involved in data interpretation, cowriting, and revisions of the manuscript. YA was responsible for sampling, data acquisition, laboratory work, statistical analyses, data interpretation, and drafting the manuscript. All authors had read and approved the final version of the manuscript.

Acknowledgments

This work was, in part, supported by a special funding provided by Damascus University. We are grateful to Dr. Weam Aldiban and Ms. Marwa Hahi for their help in patients’ eligibility assessment and sampling. Our gratitude goes to Professor Majd N. Aljamali for his invaluable technical advice. We also acknowledge the Department of Pharmaceutical Biotechnology at the National Commission for Biotechnology for hosting us in their laboratories.

Supplementary Materials

Supplementary 1. Includes Table S1. that summarizes classification of patients () according to warfarin indications.

Supplementary 2. Includes Figure S2.1. that provides a representative gel electrophoresis of PCR products encompassing the VKORC1-1639G>A polymorphism, Figure S2.2. that presents sequencing chromatograms showing examples of wild-type GG (A); homozygote AA (B), and heterozygote GA (C) of the VKORC1-1639G>A polymorphism, and Figure S2.3. that displays a sequencing chromatogram for VKORC1 rs397509427.

Supplementary 3. Includes Table S2. entitled “Frequencies of the observed VKORC1-1639G>A genotypes and those expected from Hardy-Weinberg Equilibrium (HWE) in patients on warfarin and healthy subjects.”