The association between EPCR gene p.Ser219Gly polymorphism and venous thromboembolism risk: a case–control study, meta-analysis, and a reproducibility study
Dóra Pituk
Tünde Miklós
Ágota Schlammadinger3 
Zsuzsanna Bereczky- 1Division of Clinical Laboratory Sciences, Department of Laboratory Medicine, Faculty of Medicine, Kálmán Laki Doctoral School, University of Debrecen, Debrecen, Hungary
- 2Division of Clinical Laboratory Science, Department of Laboratory Medicine, Faculty of Medicine, University of Debrecen, Debrecen, Hungary
- 3Division of Hematology, Department of Internal Medicine, Faculty of Medicine, University of Debrecen, Debrecen, Hungary
Background: The rs867186 single-nucleotide polymorphism in the PROCR gene (g.6936A > G, c.4600A > G) results in a serine-to-glycine substitution at codon 219 of endothelial protein C receptor (EPCR). We performed a case–control study followed by an updated meta-analysis of the association between this polymorphism and the risk of venous thromboembolism (VTE).
Objective and methods: We enrolled 263 VTE patients and 320 unrelated healthy controls for the case–control study. The total number of cases and controls for the meta-analysis were 5,768 and 30,017, respectively. A new online MetaGenyo Statistical Analysis System software was used to perform the current meta-analysis. Furthermore, a reproducibility study was conducted to validate our results.
Results: Among well-defined thrombosis risk factors, Factor V Leiden was more frequent in the VTE group (p < 0.001), while there was no difference in mutation frequency of prothrombin 20210G>A polymorphism between the two groups. There was no difference in the mutation frequency of Factor V Leiden and prothrombin 20210G>A between cases with and without provoking factors and cases with and without VTE recurrence. The rs867186 “G” carriership did not influence the risk of VTE [odds ratio (OR) 1.339; 95% confidence interval (CI): 0.904–1.984] in our study. No significant differences could be demonstrated among the rs867186 genotype frequencies between VTE cases with and without provoking factors (p = 0.430). PROCR rs867186 was associated with an OR of 1.72 (95% CI: 0.95–3.13, p = 0.075) in terms of VTE recurrence. In the meta-analysis, a significant association was found between EPCR Ser219Gly polymorphism and VTE under the dominant model (OR = 1.27, 95% CI: 1.11–1.46, p = 0.0006), the recessive model (OR = 1.60, 95% CI: 1.26–2.04, p = 0.0001), the GG vs. AA contrast model (OR = 1.64, 95% CI: 1.28–2.09, p = 0.0001), and the GA vs. AA contrast model (OR = 1.24, 95% CI: 1.08–1.43, p = 0.002).
Conclusion: The rs867186 was not associated with the first VTE risk in our case–control study; however, a tendency to VTE recurrence was observed. Based on the results of our reproducibility study, MetaGenyo is acceptable for meta-analysis in case of genetic epidemiology studies. Although the risk conferred by the rs867186 is mild in all meta-analyses, including ours, identifying patients carrying the minor allele might have an impact on personalized VTE risk assessment, risk-score calculation, and patient management.
1. Introduction
Two possible manifestations of venous thromboembolism (VTE) are deep venous thrombosis (DVT) and pulmonary embolism (PE) (1). VTE is a common cardiovascular disease associated with high mortality; accurately, this condition is the third most common cause of death worldwide (2). Based on the American Heart Association (AHA) Heart Disease and Stroke Statistics—2023 report, in the United States, there were approximately 1,036,000 total VTE cases in 2019 (the US population was 328 million in 2019) (3). The incidence rate for VTE varies in European countries, mainly 1.4–3.2 per 1,000 person-years (4). Intrinsic and environmental risk factors are recognized in VTE patients, especially since there are several well-established pathogenic mutations in antithrombin (AT), Protein C (PC), and Protein S (PS) genes (SERPINC1, PROC, and PROS1, respectively) that are suggested to influence the risk of thrombosis and the levels of these natural anticoagulants (5–9). While AT, PC, and PS deficiencies are rare but strong risk factors of VTE, they are included in the routine thrombophilia investigation protocol; the investigation of the endothelial protein C receptor (EPCR)—either its plasma level or the genetic variations in the PROCR gene—is currently not part of the risk assessment. VTE is a typical example of a common complex disease in which several genetic determinants may play a role, and gene–gene and gene–environment interactions have been explored in its background. Therefore, individualized risk stratification is a rational approach to establish the patient's personal thrombosis risk and the risk of recurrence. In the case of such a complex disease, the amount of risk conferred by the individual genes may not be strong one by one; however, these variations may play an enhancing effect on the development of VTE. Moreover, studies on the association of these common gene variations (having minor allele frequency, MAF, of at least 0.1) with VTE may help in the estimation of the risk at the population level.
Protein C (PC) is a natural, endogenous anticoagulant protein with anti-inflammatory properties, which interacts with the EPCR (10). The EPCR enhances the activation of protein C on the endothelial cell surface. In contrast to this transmembrane EPCR, the soluble form (sEPCR) lacks the transmembrane domain and has an opposite effect to that of the transmembrane form. The sEPCR inhibits the activation of PC and also the anticoagulant activity of the activated PC (APC) (11). There are single-nucleotide polymorphisms (SNPs), which affect the levels of circulating protein C. The rs867186 diallelic single-nucleotide polymorphism in the endothelial protein C receptor (PROCR) gene (20q11) causes a p.Ser219Gly substitution (NC_000020.11:g.35176751 A>G, c.655 A>G, p.Ser219Gly, MAF 0.10), which increases the solubility relative to membrane-bound EPCR levels and explains 75% of the variability of plasma-soluble EPCR levels and also has an effect on PC levels (12).
Recently, the results of genome-wide association (GWA) studies are also available (13). The MARTHA project described a significant association of PROCR rs867186, rs6060278 (NC_000020.11:g.35165459 T > C, MAF 0.23), and rs8119351 (NC_000020.11:g.35166602 G > A, MAF 0.09) with PC levels in VTE patients explaining 20% of the variance (14). The rs6060278 was in strong negative Linkage disequilibrium (LD) (r2 = 0.03, D′ = −1) with the p.Ser219Gly according to the study by Germain et al. (15). The adjusted odds ratio (OR) for venous thrombosis (VT) was 1.33 [95% confidence interval (CI): 1.11–1.60] in the presence of the “G” allele of rs867186 in their study. They also found that PROCR rs6088735 (NC_000020.11:g.35157873 C>T, MAF 0.22) elevated the risk of VTE. PROCR rs8119351 G>A leads to elevated PC levels, but its mechanism is unknown (16). While the rs867186 “A” allele could indicate haplotypes H1, H2, or H4, the rs867186 “G” allele exclusively associates with the H3 haplotype of EPCR. The p.Ser219Gly leads to an increased level of soluble EPCR, reflecting an increased shedding of EPCR from the endothelial membrane (17). The mechanism underlying the association between the EPCR p.Ser219Gly variant and higher soluble EPCR and variation in plasma protein C level is unclear (13). Several previously published studies (18–27) investigated the connection between EPCR gene p.Ser219Gly gene polymorphism and VTE risk with heterogeneous results, and furthermore, two meta-analyses (13, 28) were conducted.
Recently, evidence for an association between the EPCR gene p.Ser219Gly gene polymorphism and venous thrombosis is still conflicting. Previously published studies showed that rs6088735 and rs867186 could act additively on the risk of VTE (15, 29).
Therefore, this population-based case–control study and an updated, comprehensive meta-analysis aimed to investigate the potential relationship between EPCR gene p.Ser219Gly gene polymorphism and VTE risk.
2. Materials and methods
2.1. Recruitment of the study population
Our study population consisted of non-related individuals with VTE in their case histories (n = 263) and of apparently healthy unrelated adults (control subjects, n = 320).
The Division of Clinical Laboratory Science, Department of Laboratory Medicine, Faculty of Medicine, University of Debrecen, recruited the VT cases consecutively during a 1-year period between February 2014 and January 2015. Diagnosis of VT was established by standard diagnostic modalities, such as color Doppler ultrasound or phlebography at the Department of Internal Medicine. Thrombosis was described according to the guidance of the International Society of Thrombosis and Haemostasis (30). Briefly, acquired risk factors that occurred within 1 month before the diagnosis of VTE as trauma, surgery, hospitalization due to acute illness, central venous catheters, immobilization, pregnancy, oral anticoncipient use, hormonal treatment, prolonged travel, L-asparaginase treatment, the presence of positive family history, and the presence of acquired risk factors [i.e., chronic situations such as malignancy, paroxysmal nocturnal hemoglobinuria, autoimmune diseases, antiphospholipid syndrome, body mass index (BMI), varicose veins, nephrotic syndrome, heart failure, and long-term immobilization] were registered.
Healthy volunteers recruited by our study group in the same period of time were free from any known hemostasis diseases (e.g., hereditary hemorrhagic diathesis and acquired hemorrhagic disease) or drugs influencing hemostasis (vitamin K antagonist, heparin, novel oral anticoagulants, etc.); there was no venous or arterial thrombosis in their case histories (deep vein thrombosis and pulmonary embolism, myocardial infarction, ischemic stroke, and peripheral arterial disease).
Individuals with antithrombin, protein C, or protein S deficiencies with causative mutations in the background were excluded from the whole study population. Individuals with chronic illnesses (autoimmune, malignancies, chronic inflammation, renal, or hepatic diseases) were also excluded. Hypertension, diabetes mellitus, and their therapy were allowed in the study population. We also let in women taking oral anticontraceptive pills (OAC).
The work was carried out according to the principles laid down in the Declaration of Helsinki and amended in 2008 by the World Medical Assembly in Seoul, Korea. Ethical approval for the study was obtained from the National Ethical Council (3166/2012/HER), and informed consent was obtained from all participants.
2.2. Blood sampling and SNP genotyping
Blood samples were collected from the antecubital vein into vacutainer tubes (Beckton Dickinson, Franklin Lakes, NJ, USA) with anticoagulant (1/10 volume of 0.109 M citrate). DNA was isolated from the buffy coat of citrated blood samples by QIAamp DNA Blood Mini Kit (Qiagen, Hilden, Germany). Factor V Leiden mutation (FVL) and prothrombin (FII) 20210 G>A mutation were analyzed by standard molecular genetic methods (31, 32).
Twelve SNPs, located in the PROC, PROCR, PROS1, and SERPINC1 genes were analyzed using the ABI PRISM SNaPshot Multiplex Kit (Applied Biosystems, Foster City, CA, USA), among which only the effect of EPCR p.Ser219Gly (rs867186) was investigated in this study. The effects of two other PROCR polymorphisms, namely, rs8119351 and rs6088735, were also taken into consideration. Using this procedure, we identified the SNPs by utilizing fluorescently labeled dideoxynucleotides (ddNTPs) during the elongation of template-specific primers in thermocycling. PCR primers for multiplex reactions (provided upon request) and the single base extension primers (SNP-specific primers) were designed to have an annealing temperature of 58°C using Primer3Plus free online software. Primers were synthesized by Integrated DNA Technologies (Munich, Germany). After PCR amplification and purification, multiplex minisequencing was performed in a 10 µl reaction volume using 3 µl purified PCR products, 2 µl SNP-specific primers, and 5 µl of SNaPshot ready reaction mix. Sequence cycling consisted of 40 cycles of denaturation at 96°C for 10 s, primer annealing at 50°C for 5 s, and primer extension at 60°C for 30 s. The amplified products were treated with 0.7 µl SAP (shrimp alkaline phosphatase, Sigma-Aldrich, Steinheim, Germany) to prevent further binding of ddNTPs. The mixture was incubated for 60 min at 37°C followed by 15 min incubation at 75°C. The SAP-treated products (0.5 µl) were mixed with 9 µl HiDi formamide (Life Technologies, Carlsbad, CA, USA) and 0.5 µl of GeneScan 120-LIZ size standard (Applied Biosystems, Foster City, CA, USA), denatured at 95°C for 5 min, and kept on ice for at least 1 min. The fluorescently labeled fragments were separated in POP7 polymer on an ABI 3130 Genetic analyzer (Applied Biosystems). Data analyses were performed using GeneMapper Software v4.1. The detailed protocol is available upon request.
Thrombophilia investigation including plasma assays for antithrombin, protein C, and protein S deficiencies were performed using Siemens reagents (Innovance AT, Berichrom PC, Innovance free PS) on a BCS-XP automated coagulometer (Siemens, Marburg, Germany). Plasma was separated by centrifugation at 1,500 g for 20 min, and 500 μl aliquots were stored at −70°C until determination.
2.3. Statistical analysis for case–control study
The differences in the distribution of demographic characteristics and genotypic frequencies of PROCR rs867186 polymorphism between the cases and controls were examined by using the χ2 test. To assess the associations between PROCR rs867186 polymorphism and VTE risk, unconditional logistic regression for crude OR with 95% CI and adjusted OR with 95% CI were calculated. Continuous variables were expressed as median, minimum–maximum values, and interquartile ranges (IQR). IQR is the region between the 25th and 75th percentile. The statistical analysis for the case–control study was performed using Statistical Package for Social Sciences (SPSS version 23.0, Chicago, IL, USA) software and Stata 17 (Statistical Software: Release 17; Stata Corp LLC, College Station, TX, USA). A p-value of 0.05 or less was considered to indicate statistical significance.
2.4. Meta-analysis
Meta-analysis was executed on high methodology quality based on the recommended list of the Preferred Reporting Items for Systematic Review and Meta-analyses (PRISMA) (33).
2.4.1. Meta-analysis—literature search strategy
A primary search was conducted to find all relevant studies through MEDLINE using the following terms of combination, including “6936A/G,” “rs867186,” “4600AG,” “thromboembolism,” “thrombosis,” “EPCR,” and “PROCR.” Only those studies that investigated the association between EPCR gene polymorphism and venous thromboembolism risk were included in our meta-analysis. This current search was done through May 2023, and our search strategy had no language restrictions.
2.4.2. Eligibility criteria—inclusion and exclusion criteria
All the selected studies met the following criteria: (1) available studies investigate the association between EPCR gene polymorphism and venous thromboembolism risk; (2) case–control or cohort studies; (3) the studies have proper allele and/or genotype frequency information for cases and controls. The exclusion criteria were the following: (1) no properly provided information about the predictor and the outcome variable; (2) reviews, abstracts, or conference materials; (3) duplicated data.
2.4.3. Data extraction
Two authors were responsible for the literature screening and the data extraction procedure from the included studies. The following data were extracted from all included studies: first author, publication year, sample size and the number of cases and controls, ethnicities, additional pieces of information on the study population (age, gender, matching strategy), genotype, and allele distribution of cases and controls.
2.4.4. Data synthesis and statistical analysis for meta-analysis and reproducibility study
Based on the different genetic models, the ORs and their corresponding 95% CIs were used to determine the association between PROCR rs867186 polymorphism and venous thromboembolism risk. We analyzed the following genetic models: AG vs. AA, GG vs. AA, GG + AG vs. AA (dominant model), GG vs. AG + AA (recessive model), and G vs. A (allele contrast). The MAF was considered for the G allele. To analyze between-study heterogeneity, the Cochran Q-statistic and I2-statistic were used. I2 value was used to describe the percentage of variation across studies due to heterogeneity. I2 values of 25%, 50%, and 75% show low, medium, and high between-study heterogeneity, respectively (34). Cochran's Q Test p > 0.1 was considered not statistically significant heterogeneity (35, 36). The inverse variance-weighted effects meta-analysis was used if no statistically significant heterogeneity was observed; otherwise, the random-effect model (DerSimonian–Laird method) was applied (37, 38). The sensitivity analysis was performed to determine the individual effect of each selected study on the pooled analysis. Egger's test assessed the publication bias, and it was visualized by a funnel plot (39). The significance level was set at p < 0.05. The MetaGenyo web tool (40) was used to perform the statistical analysis for meta-analysis. In order to evaluate the reliability of MetaGenyo, we also performed a reproducibility study in which we compared the results obtained by the MetaGenyo web tool with those generated by the STATA statistical software performed by Dennis et al. previously (13). In the reproducibility study, we demonstrated the deviation in percentage between the results obtained by the MetaGenyo and STATA platforms. The compared parameters were as follows: the deviation from Hardy–Weinberg equilibrium (HWE), the test of association (OR and 95% CI), and between-study heterogeneity (I2 and Cohran Q p-value).
3. Results
3.1. Case–control study
After informed consent, n = 263 VTE cases and n = 320 unrelated healthy control patients were included in the study (Table 1A). Both VTE cases and controls fell within the age range of 18–85 years, and median age of VTE cases was significantly higher as compared to controls. There was no difference in gender distribution between the two groups. BMI and the frequency of positive family history of thrombosis were significantly higher in the VTE group. The frequency of smokers was significantly lower among cases. The numbers of individuals among females taking oral contraceptives were n = 18 in the cases group and n = 17 in the control group, in which the difference was not statistically significant. Within the VTE group, 23% had any provoking factors related to thrombosis. Among them, 30% were OAC users, 23.3% had trauma and/or plaster cast in their case histories, 16.6% had surgery, 11.7% had thrombosis in the postpartum period, 5% had long travel before thrombosis, immobilization was registered in 5% of cases, 3.3% had pregnancy-associated thrombosis, 1.7% had hormonal replacement, 1.7% had varicose veins, and 1.6% had thrombosis related to intravenous cannula. Cases with provoking factors in their case histories were significantly younger, and the ratio of females was significantly higher compared to cases without provoking factors (Table 1B). There were no differences between the two groups in terms of smoking, BMI, frequency of positive family history, and recurrent VT and PE. Recurrent VT was found in 28% and the frequency of PE was 27% in the VTE group. Frequencies of diabetes mellitus and hypertension were significantly higher in the cases group; however, there were no differences in these two variables between cases with and without provoking factors of VT. FV Leiden mutation was more frequent in the VTE group (p < 0.001), while there was no difference in mutation frequency of prothrombin 20210G>A polymorphism between the two groups. There was no difference in the mutation frequency of FV Leiden and prothrombin 20210G > A polymorphism between cases with and without provoking factors.
Table 1A. Demographic and clinical characteristics of the case–control study population.
Table 1B. Demographic and clinical characteristics of the cases according to the presence or absence of provoking factors.
The p.Ser219Gly polymorphism and the rs6088735 and rs8119351 in the study population met the criteria of Hardy–Weinberg equilibrium (p = 0.7687, p = 0.5521, p = 0.6701, respectively). The three SNPs showed a strong LD in our study population. Upon investigating the combination of the three SNPs, the most frequent haplotype was “ACG” (frequency 0.659), followed by “ATG” (frequency 0.232), while haplotypes “GCA” (frequency 0.094) and “GCG” (frequency 0.014) were rare variants.
In the whole study population, without any adjustments or subgroup analysis, no significant differences could be demonstrated among the rs867186 genotype frequencies between VTE cases and controls and between cases with and without provoking factors (Tables 1, 2) Subgroup analysis according to either gender or smoking did not show any significant difference in rs867186 genotype frequencies. In the subgroups of individuals with BMI less than 30 kg/m2 or age younger than 60 years (below the 75th percentile of VTE cases), the number of wild-type and mutation carriers of rs867186 did not differ significantly. If only rs6025 (FV Leiden) or rs1799963 (FII2021G>A) wild-type or carrier individuals were taken into consideration, no significant differences in genotype frequencies of the rs867186 were demonstrated. We analyzed our study group in terms of VTE recurrence and PROCR rs867186 and we found an OR of 1.72 (95% CI: 0.95–3.13; p = 0.075), which remained similar even in the adjusted model. FV Leiden and FII 20210A did not influence VTE recurrence risk in our analysis.
Table 2. Comparison of the distribution of rs867186 in cases and controls in the whole study population and in different subgroups.
Since the vast majority of our patients were on vitamin K antagonist therapy at the time of recruitment and no therapy modification was allowed during the study; protein C (PC) anticoagulant activity could not be interpreted in them. However, in healthy individuals, who were not on anticoagulant therapy, we demonstrated that the presence of rs867186 increased the PC activity significantly having extremely high plasma PC activity in homozygote mutants. In wild-type individuals, median PC activity was 112% (IQR24); in heterozygotes, median PC activity was 128.5% (IQR32); and in homozygotes (n = 3), median PC activity was 160%.
3.2. Meta-analysis
Overall, we found 53 potentially relevant candidate publications by searching the given literature database. The study selection process has been depicted in PRISMA diagram (Figure 1). For the current meta-analysis, we finally selected 11 studies (18–27). In the Identification phase, n = 19 studies were excluded due to duplication. In the screening phase, n = 20 studies were excluded. The excluded studies were: two in vitro studies involving cell lines, one study in an animal model, three meta-analyses, one study from which primary data could not be extracted, four studies with other gene mutations, and nine studies with different outcomes (acute myeloid leukemia, cerebral thrombosis, CHD patients, multiple myeloma, ovarian cancer patients, and unprovoked VTE). In the eligibility phase, four additional studies were excluded. Among them, two did not meet HWE (41, 42), and one other study had a similar study population and design as previously published by authors (43). One study was excluded because publication bias was detected (44). After the quality assessment and data extraction process, the total number of cases and controls were 5,768 and 30,017, respectively. (Table 3) After pooling all the selected studies, our current meta-analysis found a statistically significant association between EPCR p.Ser219Gly polymorphism and VTE under all applied genetic models (dominant model: OR = 1.27, 95% CI: 1.11–1.46; recessive model: OR = 1.60, 95% CI: 1.26–2.04; GG vs. AA contrast model: OR = 1.64, 95% CI: 1.28–2.09; GA vs. AA contrast model: OR = 1.24, 95% CI: 1.08–1.43) (Table 4) (Figures 2A–E). The between-study heterogeneity was low under all examined genetic models in our study. However, this may be caused by the small number of GG subjects, which can introduce a bias into this calculation and the true between-study heterogeneity is hidden, as previously published (13, 45). There were no significant differences among the included studies in the Egger test under the allele contrast (G vs. A) (p = 0.07), the recessive model (GG vs. GA + AA) (p = 0.58), the dominant model (GG + GA vs. AA) (p = 0.05), and the GG vs. AA contrast model (p = 0.43). Even the funnel plot—the plot of the log-odds ratio against the reciprocal of its standard error—did not demonstrate publication bias in these genetic models. The Egger regression asymmetry test has shown a significant difference under the GA vs. AA contrast model (p = 0.04), and the funnel plot suggested a publication bias (Figure 3). The sensitivity analysis showed no significant change in the pooled OR upon omitting individual studies (Figure 4).
Figure 1. PRISMA 2009 flow diagram showing the selection process for the PROCR rs867186 polymorphism and venous thromboembolism.
Table 3. Characteristics of the eligible studies of the association between the PROCR rs867186 genotypes and VTE.
Table 4. Summary table of the meta-analysis of association between the PROCR rs867186 genotypes and VTE.
Figure 2. Pooled OR and 95% CI for the included studies for the association between the PROCR rs867186 genotypes and venous thromboembolism (VTE). (A) Allele contrast model (G vs. A). (B) Recessive model (GG vs. GA + AA). (C) Dominant model (GG + GA vs. AA). (D) GG vs. AA model. (E) GA vs. AA model. The square size represents the effect size of individual studies. The location of the squares reflects to the OR of the corresponding individual study and the length of their vertical lines describe the appropriate confidence interval. The diamond represents the overall or summary effect.
Figure 3. Funnel plot of the EPCR p.Ser219Gly polymorphism and VTE risk in the overall study population. (A) Allele contrast model (G vs. A). (B) Recessive model (GG vs. GA + AA). (C) Dominant model (GG + GA vs. AA). (D) GG vs. AA model. (E) GA vs. AA model. The points represent the individual studies. Horizontal axis represents the OR and vertical axis represents the standard error of the individual studies. The points, which are in the upper part of the triangle represent the large studies with small standard error.
Figure 4. Summarized sensitivity analysis of association between PROCR rs867186 polymorphism and VTE. (A) Allele contrast model (G vs. A). (B) Recessive model (GG vs. GA + AA). (C) Dominant model (GG + GA vs. AA). (D) GG vs. AA model. (E) GA vs. AA model. The square size represents the effect size of individual studies. The location of the squares reflects to the OR of the corresponding individual study and the length of their vertical lines describe the appropriate confidence interval.
3.3. Reproducibility study
By definition, reproducibility is the precision of the measurement by different research groups and different experimental settings (46). In our reproducibility study, we wanted to analyze whether the results of the meta-analysis remain the same with a different statistical software. We compared the data obtained by our method (a new, online MetaGenyo Statistical Analysis System software) to those obtained and published by Dennis et al. using STATA software version 11 (13). In the comparison study, we included only those studies that were included in the analysis of Dennis et al. We compared the deviation from HWE, the results of the test of association (OR and 95% CI), and between-study heterogeneity (I2 and Cochran Q p-value). We found no discrepancies in the meta-analysis study results obtained by STATA or MetaGenyo (Tables 5, 6), and the detected deviation in percentage was 0%.
Table 5. Details of the reproducibility study between STATA version 11 software and a new online MetaGenyo Statistical Analysis System software.
Table 6. Comparison table of the previously published meta-analysis of association between the PROCR rs867186 genotypes and VTE and simulated results using MetaGenyo Statistical Analysis System software.
4. Discussion
As thrombophilia itself is a complex condition, the risk of first thrombosis or recurrence after a thrombotic episode should be calculated individually based on several pieces of evidence. One should collect data on the status of the players in the coagulation and anticoagulation processes and of environmental and acquired risk factors; after putting the pieces in the puzzle together, a final conclusion on thrombotic risk could be drawn at the personal level. Identifying the risk of recurrent thrombosis is an unmet need for several thrombophilia to decipher which patient should continue anticoagulant therapy. In our study, an elevated, however, statistically not significant VTE recurrence risk was associated with rs867186, which is worthy of further investigation in larger studies.
It is difficult to measure so many coagulation factors’ levels and perform functional assays of many proteins involved in the coagulation/anticoagulation balance since these assays are time-consuming and expensive. Nowadays, by using next-generation sequencing (NGS) technology, genetic testing has become more and more reasonable and widely available. The whole genome or exome sequencing is still not a reality, even not necessary in thrombophilia if a well-defined panel of candidate genes is created. To save time and cost, and to avoid false-positive interpretations, the so-called “thrombophilia gene panel” should include only those genes, whose alterations have proven biological consequences. In the case of PROCR rs867186, the level of sEPCR is elevated, as it was published by several studies and also confirmed by in vitro studies. Qu et al. demonstrated that in the case of stable cell lines expressing the p.219Gly variant of EPCR, the shedding of EPCR from the cell membrane was five- to sevenfold higher as compared to the p.219Ser variant after phorbol-12-myristate 13-acetate (PMA) stimulation. This led to higher sEPCR levels in the conditioned media of p.219Gly cells. They also demonstrated that the A3 haplotype (which includes the p.219Gly variant) is not only associated with an increased sEPCR level but is also associated with decreased PC activation (47). Ireland et al. showed an increased basal release of p.219Gly sEPCR coupled with higher thrombin generation. This suggested that less membrane-bound EPCR would be available for efficient PC activation (48).
In clinical studies, elevated sEPCR levels were detected in individuals carrying the p.219Gly allele, and it was associated with increased risk of thrombosis. In the study by Anastasiou et al., the sEPCR levels were not only elevated in individuals with the p.219Gly variant but they were also significantly higher in carrier patients with venous thrombosis as compared to control carriers (27). On the other hand, they found that the EPCR gene p.Ser219Gly gene polymorphism was not associated with an increased risk of thrombosis in a Greek population, most probably the gene polymorphism associated with an earlier development phase of thrombosis. In eight studies included in the meta-analysis of Dennis et al., a significant increase in sEPCR levels were associated with the p.Ser219Gly gene polymorphism (18, 19, 23, 41–43, 48, 49). The rs867186 and also rs8119351 variants of PROCR have an effect on PC levels. It was demonstrated by Athanasiadis et al. that they were responsible for 10.27% and 9.56% of the variance in PC plasma levels, respectively (16).
The genotype frequencies of EPCR p.Ser219Gly polymorphism were 75% (AA), 24% (AG), and 1% (GG) in the cases, and 81% (AA), 18% (AG), and 1% (GG) in the controls in our case–control study; and carriership of the “G” allele was without a significant effect on the risk of VTE (OR 1.339; 95% CI: 0.904–1.984). The majority of the case–control studies conducted by others and included in our meta-analysis, however, described a significant risk-increasing effect of rs867186 in their whole population (i.e., not in subgroup analysis). In the small study by Chen et al. investigating n = 65 cases and n = 71 controls, a significant risk of rs867186 was calculated in the “GA” vs. “AA” model (OR 2.75; 95%CI: 1.04–7.30) in the Chinese population (23). Significant VTE risk-increasing effect of this polymorphism was also found in another Chinese study (OR 1.912; 95% CI: 1.064–2.818 in the dominant model) (25). In a French population (n = 338 cases and n = 338 controls), the rs867186 increased the risk of VTE even in the multivariate analysis including age, sex, FV Leiden, and FII20210A (OR 1.8; 95% CI: 1.2–2.6) (18). In an Egyptian population, the carriership of the “G” allele was associated with an increased risk of VTE fourfold in the dominant model by investigating n = 90 cases and n = 90 controls (44).
In a small Turkish population (n = 111 cases and n = 73 controls), the rs867186 was without a significant VTE risk-increasing effect, although the sEPCR level was elevated in “GG” homozygous individuals (26). Similarly, in the study by Pecheniuk et al., including white and non-white subjects (n = 114 cases and n = 114 controls), no effect of this polymorphism on the risk of VTE was described (21). Uitte de Willige et al. did not find a significant association of this polymorphism with the risk of VTE either in a Dutch population involving n = 474 cases and n = 474 sex- and age-matched controls (19). They have described, however, that the low level of sEPCR was protective against VTE in their population. In a nested case–control study involving non-white and white subjects, the rs867186 was without effect on the risk of VTE (OR 0.93; 95% CI: 0.70–1.25) (42). The polymorphism, however, increased the sEPCR level, which was not associated with VTE in their study. Subgroup analysis was not performed in these studies.
In a Spanish study by Medina et al., with n = 95 VTE cases and n = 181 controls, all carriers of FV Leiden were investigated, and the rs867186 did not influence the risk of VTE (20). Even the adjustment for sex, age, FII20210A, and the presence of the other EPCR polymorphism (rs9574) did not modify the ORs. On the contrary, although there was no effect of rs867186 on the risk of VTE in the Greek population, the presence of the “G” allele was more prevalent among patients who developed VTE at younger age below 35 years (27). In a recent study, however, the rs867186 increased the risk of VTE in a large Swedish population, including middle-aged and older adults, where the HR was 1.5 (95% CI: 1.1–1.9) in “GG” homozygous individuals (n = 2,584 VTE cases and n = 26,201 without developing VTE) (24). Saposnik et al. demonstrated that the VTE risk-increasing effect of rs867186 was attributed to its effect in the male subgroup (18). In our study, rs867186 increased the risk of VTE in individuals who were wild type for rs8119351 (OR 3.2; 95% CI: 1.372–7.465), which remained significant in the multivariate model.
Based on these studies, the question is still open, whether to quantify sEPCR levels or simply determine the genotype of PROCR rs867186 from the point of view of first VTE and VTE recurrence. The level of sEPCR might be modified not only by PROCR genotype but also by different genetic and environmental factors. The reason for heterogeneity in different studies investigating the risk conferred by rs867186 and sEPCR levels might be the difference among populations in terms of influencing factors. It would be, therefore, beneficial to determine the factors—genetic and environmental—which have an effect on sEPCR levels. However, it should be taken into consideration that changes in sEPCR levels may also be consequences of different medical conditions, like thrombotic diseases, rather than being a risk factor. Large prospective studies could answer for the question, whether PROCR genotype or sEPCR levels are better predictors for VTE risk.
Due to heterogeneity in the results of the individual studies, we performed a meta-analysis in order to clarify the effect of rs867186 on the risk of VTE. Our meta-analysis of n = 5,768 VTE cases and n = 30,017 controls found a significant association between the EPCR gene p.Ser219Gly and VTE, under the dominant model (OR = 1.27, 95% CI: 1.11–1.46, p = 0.0006), the recessive model (OR = 1.60, 95% CI: 1.26–2.04, p = 0.0001), the GG vs. AA contrast model (OR = 1.64, 95% CI: 1.28–2.09, p = 0.0001), and the GA vs. AA contrast model (OR = 1.24, 95% CI: 1.08–1.43, p = 0.002). In the previous meta-analysis conducted by Dennis et al., including a smaller number of cases and controls (n = 4,821 cases and n = 6,070 controls), the main message was similar to our present analysis (OR 1.22; 95% CI: 1.11–1.33) (13). In the second meta-analysis by Li et al., n = 4,440 cases and n = 5,054 controls were analyzed, and they also found a significantly elevated risk of VTE in carriers of the “G” allele of the rs867186 (OR 1.63; 95% CI: 1.30–2.04). Here, we confirmed the results of these two previous meta-analyses. Moreover, we validated the novel, online free-access MetaGenyo Statistical Analysis System software by a reproducibility study. Based on our findings, this online platform is acceptable for meta-analysis in case of genetic epidemiology studies.
The risk conferred by the rs867186 is significantly increased but mild in all meta-analyses, including ours. However, in the concept of personalized medicine including personalized risk assessment, every factor might have an impact. Collection as many pieces of information including environmental and genetics as it is possible may serve as a basis of the personalized risk management of VTE. The availability of high-throughput genetic methods, like next-generation sequencing technology, supports the acquisition of large-scale genetic data from the patients. Upon performing gene panel sequencing, the simultaneous investigation of candidate genes can be executed. The PROCR is currently not included in the panel of the Thrombogenomics platform (50), even missing from the gene panel list of the ISTH SSC on Genomics in Thrombosis and Hemostasis Tier 1 genes (https://www.isth.org/page/GinTh_GeneLists). Including PROCR in the gene panels may help collect more experience of the role of its variants in VTE not only from the point of view of absolute VTE risk but also from the point of view of VTE recurrence and risk modifying effect in individuals with PC deficiency.
There were only two cases detected with PC deficiency in our study group during the recruitment period. The first patient had the c.811C > T (p.Arg271Trp) and the second carried the c.169C > T (p.Arg57Trp) mutation both in heterozygous forms. These mutations are listed in various databases as pathogenic/likely pathogenic variants and associated with PC deficiency and thrombotic phenotype. Both individuals carried the rs867186 in heterozygous form. The first patient had one thrombotic episode at the age of 29 years; no recurrence has been registered until now. The second had venous thrombosis first at the age of 29 years and one recurrent event at the age of 33 years while on oral anticoagulant therapy. Upon including these two PC-deficient patients in the statistical analysis, no modifications in the results were observed. To study the modifying effect of PROCR polymorphisms on the PC levels and the risk of VTE, a study directly recruiting PC-deficient individuals would be beneficial. In a previous study, Fidalgo et al. investigated Portuguese PC-deficient families and suggested that rs867186 explained differences in PC plasma levels in individuals with the same PROC genotype (51), and very high PC levels were detected in a study by Pintao et al. in case of individuals homozygous for the PROCR rs867186 (52). Homozygosity for the rs867186 was hypothesized to be the underlying condition of extremely high PC levels in a large French cohort (53). These studies, however, did not describe the potential VTE risk modifying effect of EPCR polymorphisms in PC deficiency.
The strength of the current meta-analysis is that it updated the previous results with a large-scale study (24) and added the Hungarian case–control study results. In addition to that, we also performed a reproducibility study to ensure the validity of our results.
Our study had limitations that should be acknowledged. Our case–control study might have limited generalizability through the use of strict exclusion criteria regarding the selection of the control group. Therefore, we are not sure to what extent the absence of chronic diseases and arterial thrombosis might influence the results on the risk of EPCR polymorphism on VTE.
Our meta-analysis detected potential publication bias under the GA vs. AA contrast model. Furthermore, we must consider the potential selection bias and possible misclassification of the genotype and phenotype results regarding the individual studies. We were not able to perform any further subgroup analysis because of the small group number of similar ethnicities. On the other hand, due to the non-standardized data availability of the individual studies included in our meta-analysis, we could also not create a dedicated analysis by other factors such as age, gender, smoking, alcohol consumption, and other environmental or lifestyle factors.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
DP: Data curation, Investigation, Methodology, Software, Visualization, Writing – original draft. TM: Data curation, Investigation, Methodology, Writing – review & editing. ÁS: Data curation, Investigation, Writing – review & editing. KR: Data curation, Investigation, Writing – review & editing. ZB: Data curation, Investigation, Conceptualization, Funding acquisition, Methodology, Writing – review & editing.
Funding
The authors declare financial support was received for the research, authorship, and/or publication of this article.
This work was supported by the Hungarian Scientific Research Fund (OTKA K-139293) and by the Ministry of National Economy, Hungary (GINOP-2.3.2-15-2016-00039). APC was funded by University of Debrecen.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Wilbur J, Shian B. Deep venous thrombosis and pulmonary embolism: current therapy. Am Fam Physician. (2017) 95(5):295–302.28290648
2. Klemen ND, Feingold PL, Hashimoto B, Wang M, Kleyman S, Brackett A, et al. Mortality risk associated with venous thromboembolism: a systematic review and Bayesian meta-analysis. Lancet Haematol. (2020) 7(8):e583–93. doi: 10.1016/S2352-3026(20)30211-8
3. Tsao CW, Aday AW, Almarzooq ZI, Anderson CAM, Arora P, Avery CL, et al. Heart disease and stroke statistics—2023 update: a report from the American Heart Association. Circulation. (2023) 147(8):e93–e621.36695182
4. Pastori D, Cormaci VM, Marucci S, Franchino G, Del Sole F, Capozza A, et al. A comprehensive review of risk factors for venous thromboembolism: from epidemiology to pathophysiology. Int J Mol Sci. (2023) 24(4):3169. doi: 10.3390/ijms24043169
5. Spek CA, Koster T, Rosendaal FR, Bertina RM, Reitsma PH. Genotypic variation in the promoter region of the protein C gene is associated with plasma protein C levels and thrombotic risk. Arterioscler Thromb Vasc Biol. (1995) 15(2):214–8. doi: 10.1161/01.ATV.15.2.214
6. Aiach M, Nicaud V, Alhenc-Gelas M, Gandrille S, Arnaud E, Amiral J, et al. Complex association of protein C gene promoter polymorphism with circulating protein C levels and thrombotic risk. Arterioscler Thromb Vasc Biol. (1999) 19(6):1573–6. doi: 10.1161/01.ATV.19.6.1573
7. Le Cam-Duchez V, Bagan-Triquenot A, Ménard JF, Mihout B, Borg JY. Association of the protein C promoter CG haplotype and the factor II G20210A mutation is a risk factor for cerebral venous thrombosis. Blood Coagul Fibrinolysis. (2005) 16(7):495–500. doi: 10.1097/01.mbc.0000184738.27723.b2
8. Scopes D, Berg LP, Krawczak M, Kakkar VV, Cooper DN. Polymorphic variation in the human protein C (PROC) gene promoter can influence transcriptional efficiency in vitro. Blood Coagul Fibrinolysis. (1995) 6(4):317–21. doi: 10.1097/00001721-199506000-00004
9. Giri TK, Yamazaki T, Sala N, Dahlbäck B, de Frutos PG. Deficient APC-cofactor activity of protein S Heerlen in degradation of factor Va Leiden: a possible mechanism of synergism between thrombophilic risk factors. Blood. (2000) 96(2):523–31. doi: 10.1182/blood.V96.2.523.014k04_523_531
10. Leon G, Rehill AM, Preston RJS. The protein C pathways. Curr Opin Hematol. (2022) 29(5):251–8. doi: 10.1097/MOH.0000000000000726
11. Liaw PC, Neuenschwander PF, Smirnov MD, Esmon CT. Mechanisms by which soluble endothelial cell protein C receptor modulates protein C and activated protein C function. J Biol Chem. (2000) 275(8):5447–52. doi: 10.1074/jbc.275.8.5447
12. Munir MS, Weng LC, Tang W, Basu S, Pankow JS, Matijevic N, et al. Genetic markers associated with plasma protein C level in African Americans: the atherosclerosis risk in communities (ARIC) study. Genet Epidemiol. (2014) 38(8):709–13. doi: 10.1002/gepi.21868
13. Dennis J, Johnson CY, Adediran AS, de Andrade M, Heit JA, Morange PE, et al. The endothelial protein C receptor (PROCR) Ser219Gly variant and risk of common thrombotic disorders: a HuGE review and meta-analysis of evidence from observational studies. Blood. (2012) 119(10):2392–400. doi: 10.1182/blood-2011-10-383448
14. Oudot-Mellakh T, Cohen W, Germain M, Saut N, Kallel C, Zelenika D, et al. Genome wide association study for plasma levels of natural anticoagulant inhibitors and protein C anticoagulant pathway: the MARTHA project. Br J Haematol. (2012) 157(2):230–9. doi: 10.1111/j.1365-2141.2011.09025.x
15. Germain M, Saut N, Greliche N, Dina C, Lambert J-C, Perret C, et al. Genetics of venous thrombosis: insights from a new genome wide association study. PLoS One. (2011) 6(9):e25581. doi: 10.1371/journal.pone.0025581
16. Athanasiadis G, Buil A, Souto JC, Borrell M, López S, Martinez-Perez A, et al. A genome-wide association study of the protein C anticoagulant pathway. PLoS One. (2011) 6(12):e29168. doi: 10.1371/journal.pone.0029168
17. Saposnik B, Lesteven E, Lokajczyk A, Esmon CT, Aiach M, Gandrille S. Alternative mRNA is favored by the A3 haplotype of the EPCR gene PROCR and generates a novel soluble form of EPCR in plasma. Blood. (2008) 111(7):3442–51. doi: 10.1182/blood-2007-08-104968
18. Saposnik B, Reny JL, Gaussem P, Emmerich J, Aiach M, Gandrille S. A haplotype of the EPCR gene is associated with increased plasma levels of sEPCR and is a candidate risk factor for thrombosis. Blood. (2004) 103(4):1311–8. doi: 10.1182/blood-2003-07-2520
19. Uitte de Willige S, Van Marion V, Rosendaal FR, Vos HL, de Visser MC, Bertina RM. Haplotypes of the EPCR gene, plasma sEPCR levels and the risk of deep venous thrombosis. J Thromb Haemost. (2004) 2(8):1305–10. doi: 10.1046/j.1538-7836.2004.00855.x
20. Medina P, Navarro S, Estelles A, Vaya A, Bertina RM, Espana F. Influence of the 4600A/G and 4678G/C polymorphisms in the endothelial protein C receptor (EPCR) gene on the risk of venous thromboembolism in carriers of factor V Leiden. Thromb Haemost. (2005) 94(2):389–94. doi: 10.1160/TH05-02-0089
21. Pecheniuk NM, Elias DJ, Xu X, Griffin JH. Failure to validate association of gene polymorphisms in EPCR, PAR-1, FSAP and protein S Tokushima with venous thromboembolism among Californians of European ancestry. Thromb Haemost. (2008) 99(2):453–5. doi: 10.1160/TH07-10-0607
22. Tregouet DA, Heath S, Saut N, Biron-Andreani C, Schved JF, Pernod G, et al. Common susceptibility alleles are unlikely to contribute as strongly as the FV and ABO loci to VTE risk: results from a GWAS approach. Blood. (2009) 113(21):5298–303. doi: 10.1182/blood-2008-11-190389
23. Chen XD, Tian L, Li M, Jin W, Zhang HK, Zheng CF. Relationship between endothelial cell protein C receptor gene 6936A/G polymorphisms and deep venous thrombosis. Chin Med J (Engl). (2011) 124(1):72–5.21362311
24. Manderstedt E, Hallden C, Lind-Hallden C, Elf J, Svensson PJ, Engstrom G, et al. Thrombotic risk determined by protein C receptor (PROCR) variants among middle-aged and older adults: a population-based cohort study. Thromb Haemost. (2022) 122(8):1326–32. doi: 10.1055/a-1738-1564
25. Yin G, Jin X, Ming H, Zheng X, Zhang D. Endothelial cell protein C receptor gene 6936A/G polymorphism is associated with venous thromboembolism. Exp Ther Med. (2012) 3(6):989–92. doi: 10.3892/etm.2012.510
26. Karabiyik A, Yilmaz E, Egin Y, Akar N. The effects of endothelial protein C receptor gene polymorphisms on the plasma sEPCR level in venous thrombosis patients. Turk J Haematol. (2012) 29(1):55–62.24744624
27. Anastasiou G, Politou M, Rallidis L, Grouzi E, Karakitsos P, Merkouri E, et al. Endothelial protein C receptor gene variants and risk of thrombosis. Clin Appl Thromb Hemost. (2016) 22(2):199–204. doi: 10.1177/1076029614547261
28. Li YY, Wu JJ, Yang XX, Geng HY, Gong G, Kim HJ. EPCR gene Ser219Gly polymorphism and venous thromboembolism: a meta-analysis of 9,494 subjects. Front Physiol. (2017) 8:339. doi: 10.3389/fphys.2017.00339
29. Lindström S, Wang L, Smith EN, Gordon W, van Hylckama Vlieg A, de Andrade M, et al. Genomic and transcriptomic association studies identify 16 novel susceptibility loci for venous thromboembolism. Blood. (2019) 134(19):1645–57. doi: 10.1182/blood.2019000435
30. Kearon C, Ageno W, Cannegieter SC, Cosmi B, Geersing GJ, Kyrle PA. Categorization of patients as having provoked or unprovoked venous thromboembolism: guidance from the SSC of ISTH. J Thromb Haemost. (2016) 14(7):1480–3. doi: 10.1111/jth.13336
31. Koeleman BP, Reitsma PH, Allaart CF, Bertina RM. Activated protein C resistance as an additional risk factor for thrombosis in protein C-deficient families. Blood. (1994) 84(4):1031–5. doi: 10.1182/blood.V84.4.1031.1031
32. Poort SR, Rosendaal FR, Reitsma PH, Bertina RM. A common genetic variation in the 3′-untranslated region of the prothrombin gene is associated with elevated plasma prothrombin levels and an increase in venous thrombosis. Blood. (1996) 88(10):3698–703. doi: 10.1182/blood.V88.10.3698.bloodjournal88103698
33. Liberati A, Altman DG, Tetzlaff J, Mulrow C, Gøtzsche PC, Ioannidis JP, et al. The PRISMA statement for reporting systematic reviews and meta-analyses of studies that evaluate health care interventions: explanation and elaboration. J Clin Epidemiol. (2009) 62(10):e1–34. doi: 10.1016/j.jclinepi.2009.06.006
34. Higgins JP, Thompson SG, Deeks JJ, Altman DG. Measuring inconsistency in meta-analyses. Br Med J. (2003) 327(7414):557–60. doi: 10.1136/bmj.327.7414.557
35. Whitehead A, Whitehead J. A general parametric approach to the meta-analysis of randomized clinical trials. Stat Med. (1991) 10(11):1665–77. doi: 10.1002/sim.4780101105
36. Munafò MR, Flint J. Meta-analysis of genetic association studies. Trends Genet. (2004) 20(9):439–44. doi: 10.1016/j.tig.2004.06.014
37. DerSimonian R, Laird N. Meta-analysis in clinical trials. Control Clin Trials. (1986) 7(3):177–88. doi: 10.1016/0197-2456(86)90046-2
38. Mantel N, Haenszel W. Statistical aspects of the analysis of data from retrospective studies of disease. J Natl Cancer Inst. (1959) 22(4):719–48.13655060
39. Egger M, Davey Smith G, Schneider M, Minder C. Bias in meta-analysis detected by a simple, graphical test. Br Med J. (1997) 315(7109):629–34. doi: 10.1136/bmj.315.7109.629
40. Martorell-Marugan J, Toro-Dominguez D, Alarcon-Riquelme ME, Carmona-Saez P. MetaGenyo: a web tool for meta-analysis of genetic association studies. BMC Bioinformatics. (2017) 18(1):563. doi: 10.1186/s12859-017-1990-4
41. Medina P, Navarro S, Estelles A, Vaya A, Woodhams B, Mira Y, et al. Contribution of polymorphisms in the endothelial protein C receptor gene to soluble endothelial protein C receptor and circulating activated protein C levels, and thrombotic risk. Thromb Haemost. (2004) 91(5):905–11. doi: 10.1160/TH03-10-0657
42. Yamagishi K, Cushman M, Heckbert SR, Tsai MY, Folsom AR. Lack of association of soluble endothelial protein C receptor and PROCR 6936A/G polymorphism with the risk of venous thromboembolism in a prospective study. Br J Haematol. (2009) 145(2):221–6. doi: 10.1111/j.1365-2141.2009.07612.x
43. Navarro S, Medina P, Mira Y, Estellés A, Villa P, Ferrando F, et al. Haplotypes of the EPCR gene, prothrombin levels, and the risk of venous thrombosis in carriers of the prothrombin G20210A mutation. Haematologica. (2008) 93(6):885–91. doi: 10.3324/haematol.12448
44. Zoheir N, Eldanasouri N, Abdel-Aal AA, Hosny KA, Abdel-Ghany WM. Endothelial cell protein C receptor gene 6936A/G and 4678G/C polymorphisms as risk factors for deep venous thrombosis. Blood Coagul Fibrinolysis. (2016) 27(3):259–65. doi: 10.1097/MBC.0000000000000402
45. Michael Borenstein LVH, Higgins JPT, Hannah R. Rothstein introduction to meta-analysis. Chichester: John Wiley & Sons Ltd (2009).
46. Plesser HE. Reproducibility vs. replicability: a brief history of a confused terminology. Front Neuroinform. (2017) 11:76.29403370
47. Qu D, Wang Y, Song Y, Esmon NL, Esmon CT. The Ser219-->Gly dimorphism of the endothelial protein C receptor contributes to the higher soluble protein levels observed in individuals with the A3 haplotype. J Thromb Haemost. (2006) 4(1):229–35. doi: 10.1111/j.1538-7836.2005.01676.x
48. Ireland H, Konstantoulas CJ, Cooper JA, Hawe E, Humphries SE, Mather H, et al. EPCR Ser219gly: elevated sEPCR, prothrombin F1+2, risk for coronary heart disease, and increased sEPCR shedding in vitro. Atherosclerosis. (2005) 183(2):283–92. doi: 10.1016/j.atherosclerosis.2005.02.028
49. Medina P, Navarro S, Corral J, Zorio E, Roldán V, Estellés A, et al. Endothelial protein C receptor polymorphisms and risk of myocardial infarction. Haematologica. (2008) 93(9):1358–63. doi: 10.3324/haematol.13066
50. Simeoni I, Stephens JC, Hu F, Deevi SV, Megy K, Bariana TK, et al. A high-throughput sequencing test for diagnosing inherited bleeding, thrombotic, and platelet disorders. Blood. (2016) 127(23):2791–803. doi: 10.1182/blood-2015-12-688267
51. Fidalgo T, Martinho P, Salvado R, Manco L, Oliveira AC, Pinto CS, et al. Familial thrombotic risk based on the genetic background of protein C deficiency in a Portuguese study. Eur J Haematol. (2015) 95(4):294–307. doi: 10.1111/ejh.12488
52. Pintao MC, Roshani S, de Visser MC, Tieken C, Tanck MW, Wichers IM, et al. High levels of protein C are determined by PROCR haplotype 3. J Thromb Haemost. (2011) 9(5):969–76. doi: 10.1111/j.1538-7836.2011.04256.x
Keywords: endothelial protein C receptor, venous thrombosis, PROCR gene, Ser219Gly, rs867186, polymorphism, single nucleotide, meta-analysis
Citation: Pituk D, Miklós T, Schlammadinger Á, Rázsó K and Bereczky Z (2023) The association between EPCR gene p.Ser219Gly polymorphism and venous thromboembolism risk: a case–control study, meta-analysis, and a reproducibility study. Front. Cardiovasc. Med. 10:1270093. doi: 10.3389/fcvm.2023.1270093
Received: 13 August 2023; Accepted: 6 November 2023;
Published: 22 November 2023.
Edited by:
Francis J. Castellino, University of Notre Dame, United StatesReviewed by:
Nicolas Gendron, Hôpital Européen Georges Pompidou, FranceJohn D. Imig, University of Arkansas for Medical Sciences, United States
Gunnar Engström, Lund University, Sweden
© 2023 Pituk, Miklós, Schlammadinger, Rázsó and Bereczky. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zsuzsanna Bereczky zsbereczky@med.unideb.hu