# Multilocus Genetic Risk Scores for Venous Thromboembolism Risk Assessment

José Manuel Soria, Pierre-Emmanuel Morange, Joan Vila, Juan Carlos Souto, Manel Moyano, David-Alexandre Trégouët, José Mateo, Noémie Saut, Eduardo Salas, Roberto Elosua

► To cite this version:

José Manuel Soria, Pierre-Emmanuel Morange, Joan Vila, Juan Carlos Souto, Manel Moyano, et al. Multilocus Genetic Risk Scores for Venous Thromboembolism Risk Assessment. Journal of the American Heart Association, 2014, 3 (5), pp.e001060. 10.1161/JAHA.114.001060. hal-01329408

HAL Id: hal-01329408

https://hal.sorbonne-universite.fr/hal-01329408

Submitted on 9 Jun 2016

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Multilocus Genetic Risk Scores for Venous Thromboembolism Risk Assessment

José Manuel Soria, BSc, PhD; Pierre-Emmanuel Morange, MD, PhD; Joan Vila, PhD; Juan Carlos Souto, MD, PhD; Manel Moyano, BSc; David-Alexandre Trégouët, BSc, PhD; José Mateo, MD, PhD; Noémi Saut, BSc, PhD; Eduardo Salas, MD, PhD; Roberto Elosua, MD, PhD

**Background**—Genetics plays an important role in venous thromboembolism (VTE). Factor V Leiden (FVL or rs6025) and prothrombin gene G20210A (PT or rs1799963) are the genetic variants currently tested for VTE risk assessment. We hypothesized that primary VTE risk assessment can be improved by using genetic risk scores with more genetic markers than just FVL-rs6025 and prothrombin gene PT-rs1799963. To this end, we have designed a new genetic risk score called Thrombo inCode (TiC).

**Methods and Results**—TiC was evaluated in terms of discrimination (Δ of the area under the receiver operating characteristic curve) and reclassification (integrated discrimination improvement and net reclassification improvement). This evaluation was performed using 2 age- and sex-matched case-control populations: SANTPAU (248 cases, 249 controls) and the Marseille Thrombosis Association study (MARTHA; 477 cases, 477 controls). TiC was compared with other literature-based genetic risk scores. TiC including F5 rs6025/rs118203906/rs118203905, F2 rs1799963, F12 rs1801020, F13 rs5985, SERPINC1 rs121909548, and SERPINA10 rs2232698 plus the A1 blood group (rs8176719, rs7853989, rs8176743, rs8176750) improved the area under the curve compared with a model based only on F5-rs6025 and F2-rs1799963 in SANTPAU (0.677 versus 0.575, P<0.001) and MARTHA (0.605 versus 0.576, P=0.008). TiC showed good integrated discrimination improvement of 5.49 (P<0.001) for SANTPAU and 0.96 (P=0.045) for MARTHA. Among the genetic risk scores evaluated, the proportion of VTE risk variance explained by TiC was the highest.

**Conclusions**—We conclude that TiC greatly improves prediction of VTE risk compared with other genetic risk scores. TiC should improve prevention, diagnosis, and treatment of VTE. (J Am Heart Assoc. 2014;3:e001060 doi: 10.1161/JAHA.114.001060)

**Key Words:** genetics • risk factors • tests • thrombosis • veins

Thrombosis is the formation of a blood clot inside a blood vessel that obstructs the normal flow of blood. The clinical manifestations of thrombosis are myocardial infarction, stroke, and venous thromboembolism (VTE). The latter includes pulmonary embolism and deep vein thrombosis. VTE is a common cardiovascular illness associated with high mortality¹ that affects ≈0.2% of the US and European

population annually.¹,² Consequently, it is a considerable public health concern with a high economic burden.³,⁴

VTE is a multifactorial, complex disease that results from a combination of genetic and acquired risk factors. The heritability of VTE has been estimated at about 60%.⁵ The genetic factors underlying the risk of thrombosis include some well-established mutations such as the FVL and PT

From the Unitat de Genòmica de Malalties Complexes, IIB-Sant Pau, Barcelona, Spain (J.M.S.); Inserm UMR_S 1062, F-13385, Aix-Marseille Université, Marseille, France (P.-E.M.); Grupo de Epidemiología y, Genética Cardiovascular, IMIM (Institut Hospital del Mar d'Investigacions Mèdiques), Barcelona, Spain (J.V., R.E.); CIBER Epidemiología y Salud Pública (CIBERESP), Barcelona, Spain (J.V.); Unitat d'Hemostasia i Trombosis IIB-Sant Pau, Barcelona, Spain (J.C.S., J.M.); Gendiag.exe, Barcelona, Spain (M.M., E.S.); Inserm Unité Mixte de Recherche en Santé (UMR_S) 937, ICAN Institute for Cardiometabolism and Nutrition, Université Pierre et Marie Curie Paris 6, Paris, France (D.-A.T., N.S.).

Presented in part in poster form at the XXIV ISTH Congress & 59th Annual SSC Meeting, Amsterdam, July 4, 2013, at the 46th Nordic Coagulation Meeting, Tromsø, Norway, June 7, 2013, and at the EuroThrombosis Summit 2011, October 7, 2011.

Accompanying Tables S1 through S3 are available at http://jaha.ahajournals.org/content/3/5/e001060/suppl/DC1

**Correspondence to:** José Manuel Soria, BSc, PhD, Unit of Complex Disease Genomics, IIB-Hospital de la Santa Creu i Sant Pau (IIB Sant Pau), C/ Sant Antoni Maria Claret, 167, 08025 Barcelona, Spain. E-mail: jsoria@santpau.cat

Received May 28, 2014; accepted August 29, 2014.

© 2014 The Authors. Published on behalf of the American Heart Association, Inc., by Wiley Blackwell. This is an open access article under the terms of the Creative Commons Attribution-NonCommercial License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited and is not used for commercial purposes.
mutations, which give rise to deficient anticoagulant or gain-of-function proteins.⁶ In addition, recent evidence shows several genetic variants that predispose someone to VTE by modifying components in the coagulation pathway.⁷⁻¹⁵ Moreover, several common low-penetrance gene variants and unsuspected genes have been identified by genomewide association studies (GWASs) as contributing to a risk of thromboembolic disease.¹⁶⁻¹⁸ These new variants still require proper clinical validation.

This new knowledge of the genetic profile of VTE could be used to increase the ability to more accurately predict the risk of a thrombotic event. In current clinical practice, only *FVL* and *PT* mutations are used as markers to assess a patient's risk of VTE. In a genetic risk score (GRS) described by de Haan et al,¹⁹ 5 of 31 single nucleotide polymorphisms (SNPs) linked to VTE were found to improve the predictive capacity of clinical factors including family history assessment. Moreover, the similar discriminative capacity observed by these authors with the use of 5 of 31 SNPs associated with VTE clearly indicates a need to identify an appropriate panel of genetic variants and to determine the predictive capacity of such a panel as a risk assessment score for this complex disease.

Our study was designed to compare the predictive capacity of a new GRS, Thrombo inCode (TiC), with a risk score based on family history alone and another based on *FVL* and *PT* mutations.

## Methods

### Study Populations

The SANTPAU case-control study, conducted on a Spanish population, has been described extensively.²⁰ Briefly, 248 consecutive unrelated patients who had been referred to or had visited the Hospital de la Santa Creu i Sant Pau (Barcelona, Spain) for thrombophilia screening were recruited over the period from November 1997 to April 2002. The inclusion criterion was having suffered a first thrombotic event at an age younger than 68 years. Patients were excluded if they had cancer or a history of chronic or acute liver disease or nephrotic syndrome. A medical history was obtained for each patient, including the site of thrombosis and acquired predisposing factors. The diagnosis of deep vein thrombosis of the lower limbs was established objectively by ultrasonography or ascending venography. Pulmonary embolism was diagnosed by ventilation-perfusion lung scanning, pulmonary angiography, or spiral computed tomography. Intracranial venous thrombosis was diagnosed by magnetic resonance imaging. A patient's family history was scored as positive if at least 1 other first- or second-degree family member had a venous thrombosis. As

controls, 249 unrelated, asymptomatic, and apparently healthy persons were recruited with no personal history of VTE or use of oral anticoagulants. The control group was matched to the patient group for age and sex. To avoid genetic stratification, both the case and control groups were recruited from the same geographical region; all participants were white, and all their family names were Spanish. The study protocol was approved by the hospital's institutional review board, and signed informed consent was obtained from each patient.

The Marseille Thrombosis Association study (MARTHA)²¹ is a case-control study including 1150 patients and 801 controls. The patients were unrelated and white; were recruited consecutively from the Thrombophilia Center, Hôpital de la Timone (Marseille, France) over the period from January 1994 to October 2005; had VTE; and were without known risk factors including antithrombin, protein C or protein S deficiency, homozygosity for *FVL* or for *PT*, or the presence of lupus anticoagulant. Thrombotic events including deep vein thrombosis and pulmonary embolism were documented by venography, Doppler ultrasound, spiral computed tomographic scanning angiography, and/or ventilation-perfusion lung scan. The control group was comprised of 2 subgroups: one included 475 healthy French white subjects with no personal history of cardiovascular disease (VTE was also considered) who were from the Marseille area and another group made up of 326 healthy French white heterozygotes for the FV Leiden or FII 20210A variants. In the MARTHA case-control study, controls were older than patients (47.4 versus 38.0 years, respectively), and the proportion of men was also higher among controls than among patients (47.8% versus 30.1%, respectively). To avoid the paradoxical (and confounded, by design) protective association between age and VTE and between sex and control, we matched patients and controls 1:1 by age and sex. Patients and controls were randomly matched 1:1 for the same sex and similar age (±5 years). In the original MARTHA study, 1148 patients and 801 controls were included, but for our study, we considered 477 cases and 477 controls.

### SNP Selection and Genotyping

We performed a systematic review and meta-analysis to select genetic variants that contribute to VTE risk (Table 1). Based on this information, we defined a panel, TiC (Table 1), with the variants rs6025 (*F5*, Factor V Leiden), rs118203906 (*F5*, Factor V Hong Kong), rs118203905 (*F5*, Factor V Cambridge), rs1799963 (*F2*, G20210A), rs5985 (*F13*, V34L), rs121909548 (*SERPINC1*, 384 Ala>Ser), rs2232698 (*SERPINA10*, 67 ARG>Stop), and rs1801020 (*F12*, 46 C>T) and the A1 carriers rs8176719, rs7853989, rs8176743, and rs8176750.²² It is important to note that all of these
<table><caption>**Table 1.** Genetic Variants Included in the Different Genetic Risk Scores Assessed and Coefficients (Weights) Assigned to Each Risk Factor</caption><thead><tr><th rowspan="2">SNP</th><th rowspan="2">Gene</th><th rowspan="2">Mutation</th><th rowspan="2">Risk Coefficient Assigned (B)</th><th colspan="4">GRS*</th></tr><tr><th>1</th><th>2</th><th>3</th><th>4</th></tr></thead><tbody><tr><td>rs6025, FV Leiden</td><td>F5</td><td>R506Q</td><td>1.589</td><td></td><td></td><td></td><td></td></tr><tr><td>rs118203905, FV Hong Kong</td><td>F5</td><td>R306G</td><td>1.589</td><td></td><td></td><td></td><td></td></tr><tr><td>rs118203906, FV Cambridge</td><td>F5</td><td>R306T</td><td>1.589</td><td></td><td></td><td></td><td></td></tr><tr><td>rs1799963</td><td>F2</td><td>G20210A</td><td>0.293</td><td></td><td></td><td></td><td></td></tr><tr><td>ABO</td><td>ABO</td><td>A1 carriers</td><td>0.956</td><td></td><td></td><td></td><td></td></tr><tr><td>rs8176719</td><td>ABO</td><td></td><td>—</td><td></td><td></td><td></td><td></td></tr><tr><td>rs1801020</td><td>F12</td><td>C46T</td><td>1.633</td><td></td><td></td><td></td><td></td></tr><tr><td>rs5985</td><td>F13</td><td>V34L</td><td>0.198</td><td></td><td></td><td></td><td></td></tr><tr><td>rs2232698</td><td>SERPINE10</td><td>R67X</td><td>1.358</td><td></td><td></td><td></td><td></td></tr><tr><td>rs121909548</td><td>SERPINC1</td><td>A384S</td><td>2.277</td><td></td><td></td><td></td><td></td></tr><tr><td>rs2036914</td><td>F11</td><td></td><td>0.293<br>0.519</td><td></td><td></td><td></td><td></td></tr><tr><td>rs2066865</td><td>FGG</td><td></td><td>0.344</td><td></td><td></td><td></td><td></td></tr><tr><td>rs710446</td><td>KNG1</td><td></td><td>0.182</td><td></td><td></td><td></td><td></td></tr><tr><td>rs2289252</td><td>F11</td><td></td><td>0.315<br>0.577</td><td></td><td></td><td></td><td></td></tr></tbody></table>

GRS indicates genetic risk score; SNP, single nucleotide polymorphism.

*GRS 1, FVL+PT based on F5 and rs1799963 (F2, prothrombin); GRS 2: Thrombo inCode based on F5, rs1799963 (F2, prothrombin), ABO-A1 carriers (rs8176719, rs7853989, rs8176743, rs8176750), rs1801020, rs5985, rs2232698, and rs121909548; GRS 3, de Haan et al based on F5, rs1799963 (F2, prothrombin), ABO (rs8176719), rs2036914, and rs2066865; GRS 4, expanded based on F5, rs1799963 (F2, prothrombin), ABO-A1 carriers (rs8176719, rs7853989, rs8176743, rs8176750), rs1801020, rs5985, rs2232698, rs121909548, rs2036914, rs2066865, and rs2289252.

genetic variants have functional effects on the coagulation cascade.⁷⁻¹⁵ All except the A1 carriers are gain- or loss-of-function variants.

In addition, 3 different panels of genetic variants were defined (Table 1):

1. FVL+PT: rs6025 (F5, Factor V Leiden) and rs1799963 (F2, 20210 G>A). This panel represents the genetic variants most commonly tested in current clinical practice.

2. de Haan et al panel: rs6025 (F5, Factor V Leiden), rs1799963 (F2, 20210 G>A), rs8176719 (ABO), rs2066865 (FGG, 10034 C>T), and rs2036914 (F11, 7872 C>T).

3. Extended panel: All the TiC variants plus rs2289252 (F11, 22771 T>C), rs2036914 (F11, 7872 C>T), rs710446 (KNG1, Ile581Thr), and rs2066865 (FGG, 10034 C>T). These SNPs were added because of their relationships to VTE that were detected recently by GWAS.

DNA samples from the 2 populations were genotyped. With the SANTPAU samples, the Thrombo inCode kit (Ferer inCode) was used to identify the variants included in this panel, and the remaining variants were detected by Taqman assays run in an ABI 7500 instrument. The MARTHA samples were genotyped by allele-specific polymerase chain reaction.

## Genetic Risk Score

To take into account the association strengths between the selected SNPs and VTE, we created a weighted GRS for each of the panels described (Table 1). The weights assigned to each SNP were defined a priori and based on the results of published meta-analyses in the case of FVL²³ and PT²⁴ or data from individual reports for rs2232698,¹³,²⁵ rs1801020,²⁰,²⁶,²⁷ rs5985²,⁸,²⁸,²⁹ rs2289252,³⁰,³¹ and rs2036914³¹,³² or meta-GWAS for the variants rs2066865,¹⁹,³³,³⁴ ABO,³⁵ rs121909548,¹⁴ and rs710446³⁶ (Tables S1 and S2). For FV Cambridge and Hong Kong, the same weights as FV Leiden were assigned. For the panel described by de Haan et al,¹⁹ we used the weights cited by the authors.

The genetic variants were introduced considering the genetic risk, and all were weighted in the same direction. The only genetic variant with a minor allele associated with lower thrombosis risk and odds was rs5985; in this case, we considered the common homozygote group as the risk.

## Family History

For SANTPAU populations, data were compiled on the family history of VTE.
## Statistical Analysis

Continuous variables are designated as means and standard deviations, and categorical variables are designated as proportions. Odds ratios (ORs) for the different variables linked to VTE and their 95% confidence intervals were calculated by conditional logistic regression.

We constructed different predictive models based on the *FVL+PT* GRS, the TiC GRS, the de Haan et al GRS, the extended GRS, or family history and combinations of these.

As previously mentioned, the weights assigned to all variables included in the models were defined a priori based on prior evidence.

The different scores were assessed according to the scientific statement of the American Heart Association, which describes the steps to be taken to evaluate novel risk markers in the cardiovascular field.³⁷

To assess whether the strengths of association between clinical or genetic factors and thrombosis were different for those observed in the literature (expected) and those we observed in our study (observed), we compared coefficients (ie, logarithm of the ORs) by the z test statistic; in the equation z=(b[E]-b[0])/SE, b[E] and b[0] are, respectively, the coefficients expected and observed, whereas SE is the standard error of the difference in the coefficient.

We used different measures of performance to test the quality of fit in the GRS models. *Discrimination* measures the ability of the model to discriminate between participants who will and will not have a VTE. We quantified this by calculating the area under the receiver operating characteristic curve.³⁸ This value represents an estimate of the probability that a model assigns a higher risk to those participants who will have a VTE than to those who will not have a VTE.

*Reclassification* measures how the inclusion of a new marker classifies as highest risk those participants with a VTE and as lowest risk those without a VTE. We used the methods described by Pencina et al.³⁹,⁴⁰ Integrated discrimination improvement (IDI) considers changes in the estimated VTE prediction probabilities as a continuous variable. The IDI increases when a new marker is added and thus enhances the estimate of the risk in those with VTE and decreases in those without VTE. Similarly, net reclassification improvement (NRI) requires the classification of the participants in risk categories and considers changes in the predicted probabilities of estimated VTE that imply a change from one category to another. Risk categories required to estimate the net reclassification improvement were established according to risk tertiles.

The sensitivity and specificity of the different GRSs were calculated⁴¹ using the cut points giving the highest sensitivity.

All tests were performed using R statistical software (version 3.0.1).⁴²

## Results

The sociodemographic, clinical, and genetic characteristics of the participants are listed in Table 2. In both populations, patients showed higher scores than controls in all of the GRSs examined.

Regression coefficients (ie, logarithm of the ORs) and their standard errors for associations between the different variables and VTE are shown in Table S3. In our case, not all SNPs were associated with VTE, although we included all SNPs in the GRS calculations because of consistent reports in the literature of their correlation with VTE. In Table S3, we also provide the expected regression coefficients based on a literature review and our own meta-analysis and *P* values for differences between observed and expected coefficients. We did not observe a significant difference between expected and observed coefficients except for mutations in the gene for *FVL* in the MARTHA population (expected coefficient 1.589 versus observed 0.805, *P*=0.028).

As shown in Table 3, in the SANTPAU population, the predictive model based on FVL and PT mutations showed an area under the curve of 0.575 (95% CI 0.547 to 0.604), which increased significantly when the model was based on the *TiC* GRS (0.575 versus 0.677, *P*<0.001). The de Haan et al GRS also was significantly better than the predicted model based on FVL and PT mutations (0.575 versus 0.645, *P*=0.015); however, this GRS did not improve the discriminative capacity of TiC (0.677 versus 0.645; *P*=0.346). Moreover, when we extended the TiC score (extended GRS) by adding 4 common VTE-associated SNPs (*F11* rs2289252 and rs2036914, *KNG1* rs710446, and *FGG* rs2066865), no improvement was observed over the TiC GRS area under the curve (0.677 versus 0.671, *P*=0.848).

The same approach was used to assess the validity and predictive improvement capacity of the different models when *FVL* and *PT* mutations (often used in clinical practice) were considered in addition to family history of VTE. The results in Table 4 indicate the capacity of each of the GRS models in addition to family history of VTE to improve the discrimination capacity when compared with family history of VTE alone. The discriminative capacity of TiC plus family history of VTE was not improved by the de Haan et al or extended GRSs plus family history of VTE.

Reclassification was improved by TiC, extended, or de Haan et al GRS when compared with the *FVL+PT* model, as measured by IDI (5.49, *P*<0.001; 2.56, *P*=0.009; and 2.43, *P*=0.015, respectively) (Table 3). The only GRS showing an improvement in reclassification (net reclassification improvement) over the simple *FVL+PT* model was the TiC GRS (19.17, *P*=0.002). Similar results were obtained when all GRSs plus family history of VTE were compared with family history of VTE alone (Table 4).
<table><caption>**Table 2.** Main Sociodemographic, Clinical, and Genetic Characteristics of the Study Participants</caption><thead><tr><th></th><th>Controls</th><th>Cases</th><th>P Value</th></tr></thead><tbody><tr><td>SANTPAU</td><td>n=249</td><td>n=248</td><td></td></tr><tr><td>Sex (male), n (%)</td><td>109 (44.0)</td><td>111 (44.6)</td><td>0.960</td></tr><tr><td>Age (y), mean (SD)</td><td>49.0 (14.9)</td><td>47.1 (14.0)</td><td>0.145</td></tr><tr><td>Smoker, n (%)</td><td>101 (40.7)</td><td>108 (43.7)</td><td>0.559</td></tr><tr><td>Diabetes, n (%)</td><td>9 (3.7)</td><td>14 (5.7)</td><td>0.404</td></tr><tr><td>Oral contraceptives, n (%)</td><td>74 (29.8)</td><td>83 (33.5)</td><td>0.440</td></tr><tr><td>Family history, n (%)</td><td>45 (23.2)</td><td>97 (40.9)</td><td>&lt;0.001</td></tr><tr><td>F5*†, n (%)</td><td>5 (2.02)</td><td>32 (12.9)</td><td>&lt;0.001</td></tr><tr><td>F2: rs1799963*, n (%)</td><td>7 (2.82)</td><td>19 (7.63)</td><td>0.027</td></tr><tr><td>ABO-A1 carriers/ABO*, n (%)</td><td>87 (35.7)</td><td>147 (59.0)</td><td>&lt;0.001</td></tr><tr><td>F12: rs1801020*, n (%)</td><td>5 (2.02)</td><td>15 (6.02)</td><td>0.041</td></tr><tr><td>F13: rs5985*, n (%)</td><td>139 (56.5)</td><td>146 (58.6)</td><td>0.698</td></tr><tr><td>SERPINE10: rs2232698*, n (%)</td><td>4 (1.61)</td><td>10 (4.02)</td><td>0.178</td></tr><tr><td>SERPINC1: rs121909548*, n (%)</td><td>1 (0.40)</td><td>4 (1.61)</td><td>0.372</td></tr><tr><td>F11: rs2036914</td><td></td><td></td><td></td></tr><tr><td>Hetero, n (%)</td><td>111 (46.2)</td><td>119 (48.0)</td><td>0.77</td></tr><tr><td>Homo, n (%)</td><td>54 (22.5)</td><td>43 (17.3)</td><td>0.189</td></tr><tr><td>FGG: rs2066865*, n (%)</td><td>92 (37.9)</td><td>98 (39.4)</td><td>0.804</td></tr><tr><td>KNG1: rs710446b‡, n (%)</td><td>47 (19.5)</td><td>49 (19.8)</td><td>0.966</td></tr><tr><td>F11: rs2289252</td><td></td><td></td><td></td></tr><tr><td>Hetero, n (%)</td><td>113 (47.1)</td><td>122 (49.4)</td><td>0.675</td></tr><tr><td>Homo, n (%)</td><td>39 (16.2)</td><td>46 (18.6)</td><td>0.568</td></tr><tr><td>GRS 1§, mean (SD)</td><td>0.04 (0.23)</td><td>0.23 (0.54)</td><td>&lt;0.001</td></tr><tr><td>GRS 2§, mean (SD)</td><td>0.56 (0.61)</td><td>1.10 (0.89)</td><td>&lt;0.001</td></tr><tr><td>GRS 3§, mean (SD)</td><td>0.78 (0.56)</td><td>1.16 (0.76)</td><td>&lt;0.001</td></tr><tr><td>GRS 4§, mean (SD)</td><td>1.23 (0.66)</td><td>1.76 (0.95)</td><td>&lt;0.001</td></tr><tr><td>MARTHA</td><td>n=477</td><td>N=477</td><td></td></tr><tr><td>Sex (male), n (%)</td><td>198 (41.5)</td><td>198 (41.5)</td><td>1.000</td></tr><tr><td>Age (years), mean (SD)</td><td>44.2 (13.6)</td><td>43.9 (14.0)</td><td>0.681</td></tr><tr><td>Smoker, n (%)</td><td>143 (30.2)</td><td>124 (27.7)</td><td>0.447</td></tr><tr><td>BMI (kg/m²), mean (SD)</td><td>23.8 (3.8)</td><td>25.0 (4.2)</td><td>&lt;0.001</td></tr><tr><td>Oral contraceptives, n (%)</td><td>105 (22.1)</td><td>187 (39.4)</td><td>&lt;0.001</td></tr><tr><td>F5*†, n (%)</td><td>103 (21.6)</td><td>168 (35.2)</td><td>&lt;0.001</td></tr><tr><td>F2: rs1799963*, n (%)</td><td>92 (19.3)</td><td>86 (18.0)</td><td>0.678</td></tr><tr><td>ABO-A1 carriers/ABO*, n (%)</td><td>28 (5.87)</td><td>47 (9.85)</td><td>0.030</td></tr><tr><td>F12: rs1801020*, n (%)</td><td>20 (4.19)</td><td>29 (6.08)</td><td>0.241</td></tr><tr><td>F13: rs5985*, n (%)</td><td>255 (53.5)</td><td>283 (59.3)</td><td>0.078</td></tr></tbody></table>

Continued

<table><caption>**Table 2.** Continued</caption><thead><tr><th></th><th>Controls</th><th>Cases</th><th>P Value</th></tr></thead><tbody><tr><td>SERPINE10: rs2232698*, n (%)</td><td>8 (1.68)</td><td>15 (3.14)</td><td>0.205</td></tr><tr><td>SERPINC1: rs121909548*, n (%)</td><td>3 (0.63)</td><td>1 (0.21)</td><td>0.324</td></tr><tr><td>F11: rs2036914</td><td></td><td></td><td></td></tr><tr><td>Hetero, n (%)</td><td>231 (49.7)</td><td>236 (50.1)</td><td>0.948</td></tr><tr><td>Homo, n (%)</td><td>118 (25.4)</td><td>100 (21.2)</td><td>0.155</td></tr><tr><td>FGG: rs2066865*, n (%)</td><td>178 (38.3)</td><td>220 (49.5)</td><td>0.001</td></tr><tr><td>KNG1: rs710446b‡, n (%)</td><td>92 (19.6)</td><td>82 (18.2)</td><td>0.638</td></tr><tr><td>F11: rs2289252</td><td></td><td></td><td></td></tr><tr><td>Hetero, n (%)</td><td>231 (48.4)</td><td>225 (47.2)</td><td>0.746</td></tr><tr><td>Homo, n (%)</td><td>72 (15.1)</td><td>126 (26.4)</td><td>&lt;0.001</td></tr><tr><td>GRS 1§, mean (SD)</td><td>0.40 (0.63)</td><td>0.61 (0.73)</td><td>&lt;0.001</td></tr><tr><td>GRS 2§, mean (SD)</td><td>0.67 (0.83)</td><td>0.97 (0.92)</td><td>&lt;0.001</td></tr><tr><td>GRS 3§, mean (SD)</td><td>0.85 (0.73)</td><td>1.12 (0.83)</td><td>&lt;0.001</td></tr><tr><td>GRS 4§, mean (SD)</td><td>1.34 (0.86)</td><td>1.70 (0.94)</td><td>&lt;0.001</td></tr></tbody></table>

BMI indicates body mass index; GRS, genetic risk score; hetero, heterozygosis; homo, homozygosis.

*Carriers of the risk allele.

†Carrier of any risk allele (Leiden, Hong Kong, or Cambridge).

‡Homozygotes for the risk allele.

§GRS 1, FVL+PT based on F5 and rs1799963 (F2, prothrombin); GRS 2, Thrombo inCode based on F5, rs1799963 (F2, prothrombin), ABO-A1 carriers, rs1801020, rs5985, rs2232698, and rs121909548; GRS 3, de Haan et al based on F5, rs1799963 (F2, prothrombin), ABO, rs2036914, rs2066865; GRS 4, expanded based on F5, rs1799963 (F2, prothrombin), ABO, rs1801020, rs5985, rs2232698, rs121909548, rs2036914, rs2066865, and rs2289252; ABO-A1 carriers: rs8176719, rs7853989, rs8176743, rs8176750; ABO: rs8176719.

In addition, using the GRS cut points to obtain maximal sensitivity, the sensitivity of TiC in the SANTPAU population was significantly higher than that of FVL+PT (0.85 versus 0.20%, respectively) (Table 5). The specificity of the FVL+PT GRS was higher than that of TiC (0.95 versus 0.25, respectively).

More important, the variance in VTE risk explained by the different GRSs in the SANTPAU population were 7.1%, 15.1%, 9.9%, and 13.0% for FVL+PT, TiC, de Haan et al, and extended GRSs, respectively. It is noteworthy that the TiC score explained the greatest amount of variance in thrombotic risk; in fact, it was >2-fold the variance explained by the conventional FVL+PT model.

Similar results were observed for the MARTHA population (Table 3). The AUC increased significantly with respect to the model based on FVL and PT alone when we used the TiC genetic variants (0.576 versus 0.605, P=0.008) and the extended genetic model (0.576 versus 0.629, P=0.037); however, the extended score did not improve the discriminative capacity of TiC (0.605 versus 0.629, P=0.361). Moreover, the de Haan et al GRS offered no improvement over the discriminative capacity of FVL+PT (0.576 versus 0.594,
<table><caption>**Table 3.** Predictive Capacities of the Different Models and Improvements Observed Including Different Genetic Variants Compared With the Simplest Model (FVL+PT)</caption><thead><tr><th></th><th>GRS 1<br>FVL+PT (95% CI)</th><th>GRS 2<br>TiC (95% CI)</th><th>GRS 3<br>de Haan et al (95% CI)</th><th>GRS 4<br>Extended (95% CI)</th></tr></thead><tbody><tr><td><strong>SANTPAU</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Discrimination</td><td></td><td></td><td></td><td></td></tr><tr><td>AUC</td><td>0.575 (0.547; 0.604)</td><td>0.677 (0.631; 0.724)</td><td>0.645 (0.596; 0.694)</td><td>0.671 (0.623; 0.719)</td></tr><tr><td><em>P</em> value ΔAUC</td><td>NA</td><td>&lt;0.001</td><td>0.015</td><td>&lt;0.001</td></tr><tr><td>Reclassification</td><td></td><td></td><td></td><td></td></tr><tr><td>IDI</td><td>NA</td><td>5.49 (3.35; 7.63)</td><td>2.43 (0.47; 4.39)</td><td>2.57 (0.65; 4.49)</td></tr><tr><td><em>P</em> value</td><td>NA</td><td>&lt;0.001</td><td>0.015</td><td>0.009</td></tr><tr><td>NRI</td><td>NA</td><td>19.17 (7.01; 31.33)</td><td>-5.76 (-21.84; 10.32)</td><td>-8.66 (-25.61; 8.29)</td></tr><tr><td><em>P</em> value</td><td>NA</td><td>0.002</td><td>0.483</td><td>0.317</td></tr><tr><td><strong>MARTHA</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Discrimination</td><td></td><td></td><td></td><td></td></tr><tr><td>AUC</td><td>0.576 (0.544; 0.609)</td><td>0.605 (0.570; 0.640)</td><td>0.594 (0.557; 0.631)</td><td>0.629 (0.592; 0.665)</td></tr><tr><td><em>P</em> value ΔAUC</td><td>NA</td><td>0.008</td><td>0.478</td><td>0.037</td></tr><tr><td>Reclassification</td><td></td><td></td><td></td><td></td></tr><tr><td>IDI</td><td>NA</td><td>0.96 (0.02; 1.90)</td><td>-0.06 (-0.86;0.75)</td><td>0.19 (-0.88; 1.25)</td></tr><tr><td><em>P</em> value</td><td>NA</td><td>0.045</td><td>0.889</td><td>0.730</td></tr><tr><td>NRI</td><td>NA</td><td>4.94 (-1.46; 11.33)</td><td>-5.98 (-13.98; 2.02)</td><td>-7.11 (-16.35; 2.13)</td></tr><tr><td><em>P</em> value</td><td>NA</td><td>0.130</td><td>0.143</td><td>0.131</td></tr></tbody></table>

AUC indicates area under the receiver operating characteristic curve; GRS, genetic risk score; IDI, integrated discrimination improvement; NA, not applicable; NRI, net reclassification improvement; TiC, Thrombo inCode.

P=0.47). In addition, when analyzing the reclassification, only the TiC score improved the reclassification capacity of FVL+PT, as assessed by the IDI (0.96, P=0.045).

The sensitivity and specificity of TiC scores were similar in the MARTHA population (0.850 and 0.264, respectively); however, the sensitivity and specificity of the FVL+PT GRS in

the MARTHA population differed from that found in the SANTPAU population (sensitivity 0.532 and specificity 0.591 in MARTHA).

It should be emphasized that in the MARTHA study, the extended GRS explained a high proportion of the variance in VTE risk (5.3%), whereas the FVL+PT, TiC, and de Haan et al

<table><caption>**Table 4.** Predictive Capacity of the Different Models and Improvements Observed When Including Different Genetic Variants With Respect to Family History in the SANTPAU Population</caption><thead><tr><th></th><th>Family History (95% CI)</th><th>GRS 1<br>FVL+PT (95% CI)</th><th>GRS 2<br>TiC (95% CI)</th><th>GRS 3<br>de Haan (95% CI)</th><th>GRS 4<br>Extended (95% CI)</th></tr></thead><tbody><tr><td><strong>Discrimination</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>AUC</td><td>0.589 (0.545; 0.632)</td><td>0.647 (0.602; 0.691)</td><td>0.701 (0.652; 0.749)</td><td>0.684 (0.633; 0.734)</td><td>0.700 (0.649; 0.750)</td></tr><tr><td><em>P</em> value ΔAUC</td><td>NA</td><td>&lt;0.001</td><td>&lt;0.001</td><td>0.005</td><td>0.001</td></tr><tr><td><strong>Reclassification</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>IDI</td><td>NA</td><td>3.43 (2.10; 4.76)</td><td>6.63 (4.45; 8.82)</td><td>4.28 (2.18; 6.38)</td><td>2.57 (0.40; 4.74)</td></tr><tr><td><em>P</em> value</td><td>NA</td><td>&lt;0.001</td><td>&lt;0.001</td><td>&lt;0.001</td><td>0.020</td></tr><tr><td>NRI</td><td>NA</td><td>16.04 (9.50; 22.57)</td><td>29.42 (14.33; 44.53)</td><td>6.92 (-11.40; 25.24)</td><td>0.65 (-18.57; 19.88)</td></tr><tr><td><em>P</em> value</td><td>NA</td><td>&lt;0.001</td><td>&lt;0.001</td><td>0.459</td><td>0.947</td></tr></tbody></table>

AUC indicates area under the receiver operating characteristic curve; GRS, genetic risk score; IDI, integrated discrimination improvement; NA, not applicable; NRI, net reclassification improvement; TiC, Thrombo inCode.
<table><caption>**Table 5.** Clinical Utility (Measured as the Sensitivity) and Specificity of TiC Compared With *FVL+PT* in the SANTPAU and MARTHA Populations</caption><thead><tr><th rowspan="2"></th><th rowspan="2">Selected GRS Cut Points</th><th colspan="2">SANTPAU</th><th colspan="2">MARTHA</th></tr><tr><th>Sensitivity</th><th>Specificity</th><th>Sensitivity</th><th>Specificity</th></tr></thead><tbody><tr><th scope="row">FVL+PT</th><td>0.147</td><td>0.20</td><td>0.95</td><td>0.53</td><td>0.59</td></tr><tr><th scope="row">TiC</th><td>0.099</td><td>0.85</td><td>0.25</td><td>0.85</td><td>0.26</td></tr></tbody></table>

GRS indicates genetic risk scores; TiC, Thrombo inCode.

GRSs explained 3.1%, 3.9%, and 3.7%, respectively. This observation could be attributable to the fact that the extended panel included the genetic variants included in TiC plus 4 common SNPs (*F11* rs2289252 and *rs2036914*, *KNG1* rs710446, and *FGG* rs2066865) reported to be associated with VTE in the MARTHA study.

## Discussion

As knowledge of disease improves, new biomarkers and new tests are changing the traditional concept of risk assessment. Given that thrombosis is the final outcome of many systemic disorders, it is not surprising that interest in this disease is increasing rapidly; however, accurately predicting a person's risk of developing a complex disease is very difficult. This difficulty results, in large measure, from the many risk factors that exist for a given disease. Most of these factors and their interactions are unknown. VTE is a case in point because it has a large number of risk factors related to genetic variability.⁵ Despite this substantial genetic component of VTE and the new knowledge generated by GWASs, only 2 of these variants—*FVL* and *PT*—are used conventionally in clinical settings worldwide.

In our study, we assessed the predictive validity of 3 GRSs in 2 independent populations (SANTPAU and MARTHA) according to the American Heart Association's guidelines for the evaluation of novel cardiovascular risk markers.³⁷ The proof of concept derived from our study is that the *FVL+PT* model can be greatly improved by using new genomic information. More important, by comparing the use of 3 GRSs on the same populations, we were able to show conclusively that their predictive capacities are greatly augmented using the TiC GRS.

In an initial step, we performed a systematic review of the literature and meta-analysis to select genetic variants that contribute to VTE risk, and we assigned these variants a corresponding VTE risk coefficient (Tables 1, S1, and S2). The genetic variants selected have a direct functional effect on the blood-clotting proteins,⁷⁻¹⁵ highlighting the role of coagulation in thrombosis risk. Based on these selected

genetic variants, we defined 3 weighted GRSs and determined whether they could assess the risk of VTE better than a model based on family history or *FVL+PT* alone. The weights assigned to each genetic variant were defined a priori and based on the literature.

In the second step of our study, we selected 2 independent populations in which to compare the performance of these genetic scores: the Spanish population examined in the SANTPAU case-control study and the population of the French MARTHA study, with cases and controls that were rich in FV Leiden or FII 20210A mutations. Accordingly, the SANTPAU study participants better represented the general white population and formed the basis of our study, whereas the MARTHA study participants were considered a stricter sample in which to assess the role of new genetic markers because the effects of FV Leiden and FII mutations were overrepresented.

In this study, we were able to confirm the prior finding that *FVL+PT* improves the VTE predictive capacity of family history.⁴³ We were also able to demonstrate that TiC, a GRS including 12 low-frequency, high-impact genetic risk factors, was the only score of the 3 examined that significantly improved the VTE predictive capacity of family history of VTE (Table 4) and *FVL+PT* (Table 3). This improvement consisted of better discrimination and reclassification in both the SANTPAU and MARTHA study populations.

One of the GRSs that we examined was reported by de Haan et al.¹⁹ Remarkably, these authors observed that a risk score including 31 SNPs independently associated with a VTE risk through GWAS showed a similar predictive capacity to a score including only 5 of these SNPs (*FVL*, *PT*, ABO blood group, *FGG* gene rs2066865, and *F11* gene rs2289252). Three of these 5 genetic variants (*FVL*, *PT*, ABO blood group) reported by de Haan et al were included in the TiC score. Importantly, these matching markers are low-frequency variants that show the highest individual ORs for VTE risk. This fact clearly indicates the limited value of risk scores composed of common SNPs that show low individual ORs, such as those identified by GWAS, when low-frequency variants with high ORs are also considered. In effect, the addition of 4 common SNPs (rs710446, rs2289252, rs2066865 and rs2036914) identified through GWAS to TiC (extended GRS) failed to improve the area under the curve.

This point is important regarding GWASs because the findings of both the study by de Haan et al and our study indicate the better predictive capacity of a GRS that considers low-frequency variants that returns high ORs. It should be underscored that the TiC includes rare variants featuring high individual ORs (eg, *SERPINC1* gene rs121909548 and *SERPIN10* gene rs2232698 polymorphisms). Based on these results and in agreement with a previous GWAS,²¹ it is
unlikely that common risk alleles identified through GWAS (showing a Minon Allele Frequency [MAF] >0.05 and a modestly increased VTE risk with an OR in the range 1.10 to 1.35) alone account for a large proportion of the familial risk of VTE and its clinical variability, as observed for most other human diseases investigated through GWASs.⁴⁴

These findings emphasize the need to pay special attention to low-frequency or rare variants rendering high ORs. In a recent study, the benefits of an in-depth sequencing strategy were emphasized because it permitted the identification of a rare mutation responsible for familial cases of early-onset VTE.⁴⁵,⁴⁶ The important predisposing role played by low-frequency or rare variants in VTE was highlighted in a study based on multigenerational data by Zöller et al,⁴⁷ who observed a high familial risk of VTE in a small number of siblings, suggesting segregation of rare but strong genetic risk factors.

In our study, the clinical utility (measured as the sensitivity of the GRSs) of the TiC GRS compared with *FVL+PT* was also examined. Using the cut point obtaining the highest sensitivity for *FVL+PT*, the sensitivity of *FVL+PT* was 0.20 and 0.53 of patients with VTE in SANTPAU and MARTHA, respectively, whereas using TiC at the cut point for maximal sensitivity, the sensitivity of TiC was 0.85 and 0.85 in SANTPAU and MARTHA, respectively. High sensitivity is important when the test is used to identify a serious but preventable or treatable disease such as VTE.⁴¹ This is especially true when low specificity may represent a genetic predisposition to VTE that will lead to VTE only when pro-VTE clinical conditions are also present. The identification of genetic thrombophilia has several clinical applications because most guidelines⁴⁸ advocate assigning a moderate risk of VTE in patients with thrombophilia.

Although our results point to a clear association between GRS and the risk of VTE, they should be interpreted with caution because of the limitations of our study. First, we did not include clinical data (apart from family history) in the scores because, as indicated by the guidelines,⁴⁶ several clinical scenarios exist for VTE but the genetic basis of the disease is common to all of them. Consequently, we focused our study on genetic factors and will be examining the use of the TiC GRS in combination with several clinical variables in future studies. This strategy of combining genetic and clinical data, as observed by de Haan et al,¹⁹ is likely to further improve the performance of TiC, especially its sensitivity, specificity, and the area under the receiver operating characteristic curve.

Given the lack of universally acceptable risk categories for VTE, in our net reclassification improvement analysis, we established categories based on risk tertiles, which could be an overly demanding process. Nevertheless, this limitation was resolved using the IDI index, which measures risk as a

continuum and thus should be considered a more powerful indicator for reclassification in our study.

Although we have described the differences in the design of the 2 populations, the differences in the predictive capacities of the TiC GRS in those 2 populations may raise concerns about the replicability of TiC in other populations. As mentioned, MARTHA has a special design to put the new markers under additional stress to evaluate whether they add value to *FVL+PT*. It is a case-control study enriched in carriers of *FVL* and *PT* mutations, for which >50% of participants are carriers of 1 of these genetic variants. It is expected that the area under the receiver operating characteristic curve and the IDI and net reclassification improvement results could be modest when comparing TiC with *FVL+PT* because these 2 mutations account for an important part of the risk, especially if half of the participants are carriers. Despite this big effect of these 2 mutations (due to the study design), it is important to emphasize that TiC also significantly improved the predictive capacity of *FVL+PT* in MARTHA population. Moreover, we must consider the sensitivity and specificity of TiC because those classification functions, contrary to other statistical functions, are more strongly associated with the capacity of the test than with the specific characteristic of the population in which it is tested.⁴¹ The sensitivity and sensibility of TiC were the same in the SANTPAU and MARTHA populations (sensitivity 0.85 for both; specificity 0.25 versus 0.26 in SANTPAU and MARTHA, respectively). This supports the idea that TiC will have similar performance in other populations. Considering the number, weight, and pathological relevance of the genetic variants included in TiC, the significant improvement in area under the received operating character curve and IDI in MARTHA, and the similar sensitivity and specificity in both populations, it is reasonable to believe that TiC will be replicable.

Finally, replication of the present results in other populations, evaluation of the clinical utility (measured as reduction in VTE events by the use of TiC GRS), and analysis of the cost-effectiveness of TiC remain to be accomplished.

In summary, to the best of our knowledge, this study is the first to assess the efficiency of new genetic markers of VTE, such as TiC, according to the recommendations of experts in the field of cardiovascular disease.³⁷ By determining the discrimination and reclassification capacity of these markers and their clinical utility (measured as sensitivity) in comparison to conventional recognized risk models (*FVL+PT* or family history), we were able to conclude that by using new genetic markers, especially TiC, conventional VTE risk assessment algorithms are substantially improved. New risk scores such as the TiC GRS, proposed in this paper, should allow for more tailored thromboprophylaxis strategies and improve estimates of a patient's risk of VTE.
## Sources of Funding

Supported by Spanish grants RD12/0042/0032, RD12/0042/0061, PI 11/0184, PI 12/0612, and FI12/00322.

## Disclosures

Moyano and Salas are employees of the company Gendiag.exe, S.L. assigned of the patent covering TiC GRS. Soria is a member of the Thrombo inCode Scientific Committee, advising Ferrer in Code.

## References

1. Cohen AT, Agnelli G, Anderson FA, Arcelus JI, Berqvist D, Brech JG, Greer IA, Heit JA, Hutchinson JL, Kakkar AK, Mottier D, Oger E, Samana MM, Spannagl M; for the VTE Impact Assessment Group in Europe (VITAE). Venous thromboembolism (VTE) in Europe. *Thromb Haemost.* 2007;98:756–764.

2. Anderson FA Jr, Zayarany M, Heit JA, Fidan D, Cohen AT. Estimated annual numbers of US acute-care hospital patients at risk for venous thromboembolism. *Am J Hematol.* 2007;82:777–782.

3. Stein PD, Matta F. Epidemiology and incidence: the scope of the problem and risk factors for development of venous thromboembolism. *Crit Care Clin.* 2011;27:907–932.

4. Dobesh PP. Economic burden of venous thromboembolism in hospitalized patients. *Pharmacotherapy.* 2009;29:943–953.

5. Souto JC, Almasy L, Borrell M, Blanco-Vaca F, Mateo J, Soria JM, Coll I, Felices R, Stone W, Fontcuberta J, Blangero J. Genetic susceptibility to thrombosis and its relationship to physiological risk factors: the GAIT study. Genetic analysis of idiopathic thrombophilia. *Am J Human Genet.* 2000;67:1452–1459.

6. Reitsma PH. Genetic heterogeneity in hereditary thrombophilia. *Haemostasis.* 2000;30(suppl 2): 1–10.

7. Kanaji T, Okamura T, Osaki K, Kuroiwa M, Shimoda K, Hamasaki N, Niho Y. A common genetic polymorphism (46 C to T substitution) in the 5'-untranslated region of the coagulation factor XII gene is associated with low translation efficiency and decrease in plasma factor XII level. *Blood.* 1998;91:2010–2014.

8. Van HylckamaVlieg A, Komansin N, Ariëns RAS, Poort SR, Grant PJ, Bertina RM, Rosendaal FR. Factor XIII Val34Leu polymorphism, factor XIII antigen levels and activity and the risk of deep venous thrombosis. *BJH.* 2002; 119:169–175.

9. Williamson D., Brown K, Luddington R, Baglin C, Baglin T. Factor V Cambridge: a new mutation (Arg306Thr) associated with resistance to activated protein C. *Blood* 1998;91:1140–1144.

10. Norstrom E, Thorelli E, Dahlbäck B. Functional characterization of recombinant FV Hong Kong and FV Cambridge. *Blood.* 2002; 100:524–530.

11. Bertina R, Koeleman BPC, Koster T, Rosendaal FR, Driven RJ, De Ronde H, Van der Velden PA, Reitsma PH. Mutation in blood coagulation factor V associated with resistance to activated protein C. *Nature.* 1994;369:64–67.

12. 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:3698–3703.

13. Corral J, Gonzalez-Conejero R, Soria JM, Gonzalez-Porras JR, Pérez-Ceballos E, Lecumberri R, Roldán V, Souto JC, Minano A, Hernández-Espinosa D, Alberca I, Fontcuberta J, Vicente V. A nonsense polymorphism in the protein Z-dependent protease inhibitor increases the risk for venous thrombosis. *Blood.* 2006;108:177–183.

14. Corral J, Hernandez-Espinosa D, Soria JM, González-Conejero R, Ordonez A, González-Porras JR, Pérez-Ceballos E, Lecumberri R, Sánchez I, Roldan V, Mateo J, Minano A, González M, Alberca I, Fontcuberta J, Vicente V. Antithrombin Cambridge II (A384S): an nderestimated genetic risk factor for venous thrombosis. *Blood.* 2007;109:4258–4263.

15. O'Donnell J, Boulton FE, Manning RA, Laffan MA. Amount of H antigen expressed on circulating von Willebrand factor is modified by ABO blood group genotype and is a major determinant of plasma von Willebrand factor antigen levels. *Arterioscler Thromb Vasc Biol.* 2002;22:335–341.

16. Morange PE, Tregouet DA. Lessons from genome-wide association studies in venous thrombosis. *J Thromb Haemost.* 2011;9 (suppl 1): 258–264.

17. Tang W, Teichert M, Chasman DI, Heit JA, Morange PE, Li G, Pankratz N, Leebeek FW, Paré G, de Andrade M, Tzourio C, Psaty BM, Basu S, Ruiter R,

Rose L, Armasu SM, Lmley T, Heckbert SR, Uitterlinde AG, Lathrop M, Rice KM, Cushman M, Hofman A, Lambert JC, Glazer NL, Pankow JS, Witterman JC, Amouyel P, Bis JC, Bovill EG, Kong X, Tracy RP, Boerwinkle E, Rotter JI, Trégouët DA, Loth DW, Stricker BHC, Ridker PM, Folsom AR, Smith NL. A genome-wide association study for venous thromboembolism: the extended cohorts for heart ad aging research in genomic epidemiology (CHARGE) consortium. *Genet Epidemiol.* 2013;37:512–521.

18. Sabater-Lleal M, Martinez-Perez A, Buil A, Folkersen L, Souto JC, Bruzelius M, Borrell M, Odeberg J, Silveira A, Eriksson P, Almasy L, Hamsten A, Soria JM. A genome-wide association study identifies KNG1 as a genetic determinant of plasma factor XI level and activated partial thromboplastin time. *Arterioscler Thromb Vasc Biol.* 2012;32:2008–2016.

19. De Haan HG, Bezemer ID, Doggen CJM, Le Cessie S, Reitsma PH, Arellano AR, Tong CH, Devlin JJ, Bare LA, Rosendaal FR, Vossen CY. Multiple SNP testing improves risk prediction of first venous thrombosis. *Blood.* 2012;120:656–663.

20. Tirado I, Soria JM, Mateo J, Oliver A, Souto JC, Santamaria A, Felices R, Borrell M, Fontcuberta J. Association after linkage analysis indicates that homozygosity for the 46C>T polymorphism in the F12 gene is a genetic risk factor for venous thrombosis. *Thromb Haemost.* 2004;91:899–904.

21. Trégouët DA, Heath S, Saut N, Biron-Andreani C, Schved JF, Pernod G, Galan P, Drouet L, Zelenika D, Juan-Vague I, Alessi MC, Tiret L, Lathrop M, Emmerich J, Morange PE. Common susceptibility alleles are unlikely to contribute as strongly as the FV and AB0 loci to VTE risk: results from a GWAS approach. *Blood.* 2009;113:5298–5303.

22. Yip SP. Single-tube multiplex PCR-SSCP analysis distinguishes 7 common ABO alleles and readily identifies new alleles. *Blood.* 2000;95:1487–1492.

23. Juul K, Tybjaerg-Hansen A, Schnohr P, Nordestgaard BG. Factor V Leiden and the risk for venous thromboembolism in the adult Danish population. *Ann Intern Med.* 2004;140:330–337.

24. Wu O, Robertson L, Twaddle S, Lowe G, Clark P, Walker I, Brenkel I, Greaves M, Langhorne P, Regan L, Greer I; for The Thrombosis: Risk and Economic Assessment of Thrombophilia Screening—TREATS-Study. Screening for thrombophilia in high-risk situations: a meta-analysis and cost-effectiveness analysis. *Br J Haematol.* 2005;131:80–90.

25. de Van Water N, Tan T, Ashton F, O'Grady A, Day T, Browett P, Ockelford P, Harper P. Mutations within the protein Z-dependent protease inhibitor gene are associated with venous thromboembolic disease: a new form of thrombophilia. *Br J Haematol.* 2004;127:190–194.

26. Reuner KH, Jenetzky E, Aleu A, Litfin F, Mellado P, Kloss M, Jüttler E, Grau AJ, Rickmann H, Patscheke H, Lichy C. Factor XII c46T gene polymorphism and the risk of cerebral venous thrombosis. *Neurology.* 2008;70:129–132.

27. Cochery-Nouvellon E, Mercier E, Lissalde-Lavigne G, Daurès J-P, Quéré I, Dauzat M, Marès P, Gris C. Homozygosity for the C46T polymorphism of the F12 gene is a risk factor for venous thrombosis during the first pregnancy. *J Thromb Haemost.* 2007;5:700–707.

28. Catto AJ, Kohler HP, Coore J, Mansfield MW, Stickland MH, Grant PJ. Association of a common polymorphism in the factor XIII gene with venous thrombosis. *Blood.* 1999;93:906–908.

29. Wells PS, Anderson JL, Scarvelis DK, Doucette SP, Gagnon F. Factor XIII Val34Leu variant is protective against venous thromboembolism: a HuGE review and meta-analysis. *Am J Epidemiol.* 2006;164:101–109.

30. Arellano AR, Bezemer ID, Tong CH, Catanese JJ, Devlin JJ, Reitsma PH, Bare LA, Rosendaal FR. Gene variants associated with venous thrombosis: confirmation in the MEGA study. *J Thromb Haemost.* 2010;8:1132–1134.

31. Austin H, De Staercke C, Lally C, Bezemer ID, Rosendaal FR, Hooper WC. New gene variants associated with venous thrombosis: a replication study in White and black Americans. *J Thromb Haemost.* 2011;9:489–495.

32. Bezemer ID, Bare LA, Doggen CJ, Arellano AR, Tong C, Rowland CM, Catanese J, Young BA, Reitsma PH, Devlin JJ, Rosendaal FR. Gene variants associated with deep vein thrombosis. *JAMA.* 2008;299:1306–1314.

33. Uitte de Willige S, de Visser MC, Houwing-Duistermaat JJ, Rosendaal FR, Vos HL, Bertina RM. Genetic variation in the fibrinogen gamma gene increases the risk of deep venous thrombosis by reducing plasma fibrinogen gamma levels. *Blood.* 2005;106:4176–7183.

34. Grünbacher G, Weger W, Max-Neuhold E, Pilger E, Köppel H, Wascher T, März W, Renner W. The fibrinogen gamma (FGG) 10034C>T polymorphism is associated with venous thrombosis. *Thromb Res.* 2007;121:33–36.

35. Tirado I, Mateo J, Soria JM, Oliver A, Martínez-Sánchez E, Vallvé C, Borrell M, Urruti T, Fontcuberta J. The ABO blood group genotype and factor VIII levels as independent risk factors for venous thromboembolism. *Thromb Haemost.* 2005;93:468–474.

36. Morange PE, Oudot-Mellakh T, Cohen W, Germain M, Saut N, Antoni G, Alessi MC, Bertrand M, Dupuy AM, Letenneur L, Lathrop M, Lopez LM, Lambert JC,
Emmerich J, Amouyel P, trégoüet DA. KNG1 ile581Thr and susceptibility to venous thrombosis. *Blood*. 2011;117:3692–3694.

37. Hlatky MA, Greenland P, Arnett DK, Ballantyne CM, Criqui MH, Elkind MSV, Go AS, Harrell FE Jr, Hong Y, Howard BV, Howard VJ, Hsue PY, Kramer CM, McConnell JP, Normand SLT, O'Donnell CJ, Smith SC Jr, Wilson PWF; on behalf of the American Heart Association Expert Panel on Subclinical Atherosclerotic Disease and Emerging risk factors and the Stroke Council. Criteria for evaluation of novel cardiovascular risk: a scientific statement from the American Heart Association. *Circulation*. 2009;119:2408–2416.

38. Hanley JA, Haijian-Tilaki KO. Sampling variability of nonparametric estimates of the areas under receiver operating characteristic curves: an update. *Acad Radiol*. 1997;4:49–58.

39. Pencina MJ, D'Agostino RB Sr, D'Agostino RB Jr, Vasan RS. Evaluating the added predictive ability of a new marker: from area under the ROC curve to reclassification and beyond. *Stat Med*. 2008;27:157–172.

40. Pencina MJ, D'Agostino RB, Pencina KM, Janssens AC, Greeland P. Interpreting incremental value of markers added to risk prediction models. *Am J Epidemiol*. 2012;176:473–481.

41. Lalkhen AG, McCluskey A. Clinical tests: sensitivity and specificity. *Cont Educ Anaesth Crit Care Pain*. 2008;8:221–223.

42. R Development Core Team. *R: a Language and Environment for Statistical Computing*. Vienna, Austria: R Foundation for Statistical Computing; 2010. Available at: http://www.R-project.org/. Accessed October 17, 2014.

43. Varga EA, Kujovich JL. Management of inherited thrombophilia: guide for genetics professionals. *Clin Genet*. 2012;81:7–17.

44. Manolio TA, Collins FS, Cox NJ, Goldstein DB, Hindorff LA, Hunter DJ, McCarthy MI, Ramos EM, Cardon LR, Chakravarti A, Cho JH, Guttmacher AE, Kong A, Kruglyak L, Mardis E, Rotimi CN, Slatkin M, Valle D, Whittemore AS, Boehnke M, Clark AG, Eichler EE, Gibson G, Haines JL, Mackay TF, McCarroll SA, Visscher PM. Finding the missing heritability of complex diseases. *Nature*. 2009;461:747–753.

45. Antoni G, Morange PE, Luo Y, Saut N, Burgos G, Heath S, Germain M, Biron-Andreani C, Schved JF, Pernod G, Galan P, Zelanika D, Alessi MA, Drouet L, Visvikis-Siest S, Wells PS, Lathrop M, Emmerich J, Tregouet DA, Gagnon F. A multi-stage multi-design strategy provides strong evidence that the BAI3 locus is associated with early-onset venous thromboembolism. *J Thromb Hemost*. 2010;8:2671–2679.

46. Wu C, Dwivedi DJ, Pepler L, Lysov Z, Waye J, Julian J, Desch K, Ginsburg D, Weitz JI, Kearon C, Liaw P. Targeted gene sequencing identifies variants in the protein C and endothelial protein C receptor genes in patients with unprovoked venous thromboembolism. *Arterioscler Thromb Vasc Biol*. 2013;33:2674–2681.

47. Zöller B, Li X, Sundquist J, Sundquist K. Age- and gender-specific familial risks for venous thromboembolism: a nationwide epidemiological study base don hospitalizations in Sweden. *Circulation*. 2011;124:1012–1020.

48. Antithrombotic therapy and prevention of thrombosis, 9th ed: American College of Chest Physicians Evidence-Based Clinical Practice Guidelines. *Chest*. 2012;14(2 suppl):e1S–e801S.
# SUPPLEMENTAL MATERIAL

## MULTILOCUS GENETIC RISK SCORES FOR VENOUS THROMBOEMBOLISM RISK ASSESSMENT

Soria et al

<table><caption>**Table 1.** Genetic variants included in the different genetic risk scores (GRS) assessed: coefficients (weights) assigned to each genetic risk factor, and allele frequency and Hardy-Weinberg equilibrium (H-W) data obtained in the two studies, Sant Pau and MARTHA.</caption><thead><tr><th rowspan="2">SNP</th><th rowspan="2">Gene</th><th rowspan="2">Mutation</th><th rowspan="2">Genetic risk group*</th><th rowspan="2">Risk coefficient assigned</th><th colspan="4">GRS</th><th colspan="2">SANT PAU</th><th colspan="2">MARTHA</th></tr><tr><th>1</th><th>2</th><th>3</th><th>4</th><th>Rare allele Frequency</th><th>H-W</th><th>Rare allele frequency</th><th>H-W</th></tr></thead><tbody><tr><td>rs6025</td><td>F5 Leiden</td><td>R506</td><td>A carriers</td><td>1.589</td><td></td><td></td><td></td><td></td><td>0.01</td><td>0.873</td><td>0.11</td><td>0.008</td></tr><tr><td>rs118203905</td><td>F5 Hong Kong</td><td>F5 Hong Kong</td><td>A carriers</td><td>1.589</td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.00</td><td>0.982</td></tr><tr><td>rs118203906</td><td>F5 Cambridge</td><td>F5 Cambridge</td><td>C carriers</td><td>1.589</td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.00</td><td>---</td></tr><tr><td>rs1799963</td><td>F2</td><td>G20210A</td><td>A carriers</td><td>0.293</td><td></td><td></td><td></td><td></td><td>0.01</td><td>0.822</td><td>0.10</td><td>0.020</td></tr><tr><td>ABO**</td><td>ABO</td><td></td><td>A1 carriers</td><td>0.956</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>rs8176719</td><td>ABO</td><td></td><td></td><td>--</td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.36</td><td>0.262</td></tr><tr><td>rs1801020</td><td>F12</td><td>C46T</td><td>TT homoz.</td><td>1.633</td><td></td><td></td><td></td><td></td><td>0.21</td><td>0.033</td><td>0.23</td><td>0.223</td></tr><tr><td>rs5985</td><td>F13</td><td>V34L</td><td>GG homoz.</td><td>0.198</td><td></td><td></td><td></td><td></td><td>0.25</td><td>0.832</td><td>0.27</td><td>0.929</td></tr><tr><td>rs2232698</td><td>SERPINE10</td><td>R67X</td><td>T carriers</td><td>1.358</td><td></td><td></td><td></td><td></td><td>0.01</td><td>0.898</td><td>0.01</td><td>0.854</td></tr><tr><td>rs121909548</td><td>SERPINC1</td><td>A384S</td><td>T carriers</td><td>2.277</td><td></td><td></td><td></td><td></td><td>0.00</td><td>0.975</td><td>0.00</td><td>0.945</td></tr><tr><td rowspan="2">rs2036914</td><td rowspan="2">FXI</td><td rowspan="2"></td><td>CT heter.</td><td>0.293</td><td></td><td></td><td></td><td></td><td rowspan="2">0.46</td><td rowspan="2">0.293</td><td rowspan="2">0.50</td><td rowspan="2">0.890</td></tr><tr><td>TT homoz.</td><td>0.519</td><td></td><td></td><td></td><td></td></tr><tr><td>rs2066865</td><td>FGG</td><td></td><td>TT homoz.</td><td>0.344</td><td></td><td></td><td></td><td></td><td>0.21</td><td>0.739</td><td>0.23</td><td>0.028</td></tr><tr><td>rs710446</td><td>KNG1</td><td></td><td>CC homoz.</td><td>0.182</td><td></td><td></td><td></td><td></td><td>0.44</td><td>0.831</td><td>0.42</td><td>0.131</td></tr><tr><td rowspan="2">rs2289252</td><td rowspan="2">FXI</td><td rowspan="2"></td><td>CT heter.</td><td>0.315</td><td></td><td></td><td></td><td></td><td rowspan="2">0.40</td><td rowspan="2">0.788</td><td rowspan="2">0.40</td><td rowspan="2">0.655</td></tr><tr><td>TT homoz.</td><td>0.577</td><td></td><td></td><td></td><td></td></tr></tbody></table>

<footnote>*Homoz.=Homozygous; Heter.=Heterozygous. ** ABO: A1 carriers (rs8176719, rs7853989, rs8176743, rs8176750)</footnote>
<table><caption>**Table 2:** Data sources and results of the meta-analyses undertaken by the authors to assign a risk coefficient to each variant included in the genetic risk score.</caption><thead><tr><th></th><th>Gen</th><th>Reference</th><th>Number of cases</th><th>Number of controls</th><th>OR</th><th>Lower confidence interval</th><th>Upper confidence interval</th><th>Risk Coefficient assigned</th></tr></thead><tbody><tr><td>rs6025-FVL</td><td>F5</td><td>1</td><td></td><td></td><td></td><td></td><td></td><td>1.589</td></tr><tr><td>rs118203905-F5 Hong-Kong</td><td>F5</td><td></td><td></td><td></td><td></td><td></td><td></td><td>1.589</td></tr><tr><td>rs118203906-F5 Cambridge</td><td>F5</td><td></td><td></td><td></td><td></td><td></td><td></td><td>1.589</td></tr><tr><td>rs1799963</td><td>F2</td><td>2</td><td></td><td></td><td></td><td></td><td></td><td>0.293</td></tr><tr><td>AB0</td><td>AB0</td><td>3</td><td></td><td></td><td></td><td></td><td></td><td>0.956</td></tr><tr><td>rs8176719</td><td>AB0</td><td>3,4</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>rs1801020</td><td>F12</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>5</td><td>250</td><td>250</td><td>4.82</td><td>1.50</td><td>15.60</td><td></td></tr><tr><td></td><td></td><td>6</td><td>--</td><td>--</td><td>4.57</td><td>1.55</td><td>13.40</td><td></td></tr><tr><td></td><td></td><td>7</td><td>32463</td><td>32463</td><td>5.99</td><td>2.10</td><td>17.30</td><td></td></tr><tr><td>Meta-analysis</td><td></td><td></td><td></td><td></td><td>5.12</td><td></td><td></td><td>1.633</td></tr><tr><td>rs5985</td><td>F13</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>8</td><td>221</td><td>254</td><td>0.63</td><td>0.38</td><td>0.82</td><td></td></tr><tr><td></td><td></td><td>9</td><td>3165</td><td>4909</td><td>0.85</td><td>0.77</td><td>0.95</td><td></td></tr><tr><td></td><td></td><td>10</td><td>475</td><td>475</td><td>0.80</td><td>0.60</td><td>1.00</td><td></td></tr><tr><td>Meta-analysis</td><td></td><td></td><td></td><td></td><td>0.82</td><td></td><td></td><td>-0.198 (0.198)*</td></tr><tr><td>rs2232698</td><td>SERPINE10</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>11</td><td>1018</td><td>1018</td><td>3.35</td><td>1.30</td><td>8.60</td><td></td></tr><tr><td></td><td></td><td>12</td><td>250</td><td>250</td><td>5.70</td><td>1.25</td><td>26.0</td><td></td></tr><tr><td>Meta-analysis</td><td></td><td></td><td></td><td></td><td>3.89</td><td></td><td></td><td>1.358</td></tr><tr><td>rs121909548</td><td>SERPINC1</td><td>13</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>rs2036914</td><td>F11</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Bezemer (2008) : Heteroz. Homoz.</td><td></td><td>14</td><td>1314</td><td>2877</td><td>1.38<br>1.71</td><td>1.17<br>1.43</td><td>1.64<br>2.05</td><td></td></tr><tr><td>Austin (2011) : Heteroz. Homoz.</td><td></td><td>15</td><td>1076</td><td>1239</td><td>1.20<br>1.60</td><td>0.92<br>1.20</td><td>1.70<br>2.30</td><td></td></tr><tr><td>Meta-analysis: Heteroz. Homoz.</td><td></td><td></td><td></td><td></td><td>1.34<br>1.68</td><td></td><td></td><td>0.293<br>0.519</td></tr><tr><td>rs2066865</td><td>FGG</td><td>16, 17, 18</td><td></td><td></td><td></td><td></td><td></td><td>0.344</td></tr><tr><td>rs710446</td><td>KNG1</td><td>19</td><td></td><td></td><td></td><td></td><td></td><td>0.182</td></tr><tr><td>rs2289252</td><td>F11</td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.315</td></tr></tbody></table>
<table><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.577</td></tr><tr><td>Arellano (2010): Heteroz.</td><td></td><td>20</td><td>3921</td><td>4634</td><td>1.37</td><td>1.24</td><td>1.51</td><td></td></tr><tr><td>Homoz.</td><td></td><td></td><td></td><td></td><td>1.78</td><td>1.57</td><td>2.01</td><td></td></tr><tr><td>Austin (2011): Heteroz.</td><td></td><td>21</td><td>1076</td><td>1239</td><td>1.40</td><td>1.10</td><td>1.80</td><td></td></tr><tr><td>Homoz.</td><td></td><td></td><td></td><td></td><td>1.80</td><td>1.30</td><td>2.50</td><td></td></tr><tr><td>Meta-analysis: Heteroz.</td><td></td><td></td><td></td><td></td><td>1.37</td><td></td><td></td><td>0.315</td></tr><tr><td>Homoz.</td><td></td><td></td><td></td><td></td><td>1.78</td><td></td><td></td><td>0.577</td></tr></table>

<footnote>*In the case of the variant rs5985, the rare variant was associated with lower probability of presenting a VTE. To build the GRS the common variant was considered to be associated with a higher probability of presenting of VTE.</footnote>
<table><caption>**Table 3.** Observed and expected strengths of associations between selected clinical or genetic factors and venous thromboembolic events and p-values for the differences between observed and expected values.</caption><thead><tr><th rowspan="3"></th><th rowspan="3">Expected<br>coefficients</th><th colspan="3">SANT PAU</th><th colspan="3">MARTHA</th></tr><tr><th colspan="3">Observed</th><th colspan="3">Observed</th></tr><tr><th>Coefficients</th><th>Standard<br>Error</th><th>p-value*</th><th>Coefficients</th><th>Standard<br>Error</th><th>p-value*</th></tr></thead><tbody><tr><td>Family History</td><td>1.185</td><td>0.711</td><td>0.240</td><td>0.323</td><td>NA</td><td>NA</td><td>NA</td></tr><tr><td>F5**†</td><td>1.589</td><td>2.803</td><td>0.755</td><td>0.421</td><td>0.805</td><td>0.178</td><td>0.028</td></tr><tr><td>F2- rs1799963†</td><td>0.293</td><td>0.773</td><td>0.489</td><td>0.624</td><td>0.378</td><td>0.199</td><td>0.831</td></tr><tr><td>ABO (A1)†</td><td>0.956</td><td>0.997</td><td>0.224</td><td>0.927</td><td>0.428</td><td>0.280</td><td>0.346</td></tr><tr><td>F12- rs1801020†</td><td>1.633</td><td>0.859</td><td>0.621</td><td>0.533</td><td>0.305</td><td>0.345</td><td>0.054</td></tr><tr><td>F13- rs5985†</td><td>0.198</td><td>-0.078</td><td>0.224</td><td>0.538</td><td>0.171</td><td>0.152</td><td>0.929</td></tr><tr><td>SERPINE10- rs2232698†</td><td>1.358</td><td>0.937</td><td>0.753</td><td>0.780</td><td>0.572</td><td>0.502</td><td>0.434</td></tr><tr><td>SERPINC1- rs121909548†</td><td>2.277</td><td>1.476</td><td>1.174</td><td>0.733</td><td>-0.983</td><td>1.264</td><td>0.197</td></tr><tr><td>F11- rs2036914 Hetero</td><td>0.293</td><td>-0.039</td><td>0.273</td><td>0.543</td><td>0.228</td><td>0.207</td><td>0.875</td></tr><tr><td>Homo</td><td>0.519</td><td>0.028</td><td>0.395</td><td>0.534</td><td>0.206</td><td>0.258</td><td>0.544</td></tr></tbody></table>
<table><tr><td>FGG- rs2066865 †</td><td>0.344</td><td>-0.095</td><td>0.227</td><td>0.334</td><td>0.328</td><td>0.151</td><td>0.958</td></tr><tr><td>KNG1- rs710446b ‡</td><td>0.182</td><td>-0.206</td><td>0.278</td><td>0.485</td><td>-0.086</td><td>0.190</td><td>0.481</td></tr><tr><td>F11- rs2289252 Hetero</td><td>0.315</td><td>0.186</td><td>0.299</td><td>0.829</td><td>0.191</td><td>0.190</td><td>0.744</td></tr><tr><td>Homo</td><td>0.577</td><td>-0.017</td><td>0.403</td><td>0.461</td><td>0.965</td><td>0.262</td><td>0.459</td></tr></table>

*for the difference between observed and expected coefficients.

** Carrier of any risk allele (Leiden, Hong-Kong or Cambridge).

NA indicates not applicable.
References:

1. Juul K, Tybjaerg-Hansen A, Schnohr P, Nordestgaard BG. Factor V Leiden and the risk for venous thromboembolism in the adult Danish population. Ann Intern Med. 2004;140:330-337.

2. Wu O, Robertson L, Twaddle S, Lowe G, Clark P, Walker I, Brenkel I, Greaves M, Langhorne P, Regan L, Greer I; Thrombosis: Risk and Economic Assessment of Thrombophilia Screening (TREATS) Study. Screening for thrombophilia in high-risk situations: a meta-analysis and cost-effectiveness analysis. Br J Haematol. 2005;131:80-90.

3. Tirado I, Mateo J, Soria JM, Oliver A, Martínez-Sánchez E, Vallvé C, Borrell M, Urrutia T, Fontcuberta J. The ABO blood group genotype and factor VIII levels as independent risk factors for venous thromboembolism. Thromb Haemost. 2005;93:468-474.

4. De Haan HG, Bezemer ID, Doggen CJM, Le Cessie S, Reitsma PH, Arellano AR, Tong CH, Devlin JJ, Bare LA, Rosendaal FR, Vossen CY. Multiple SNP testing improves risk prediction of first venous thrombosis. Blood. 2012;120:656-663.

5. Tirado I, Soria JM, Mateo J, Oliver A, Souto JC, Santamaría A, Felices R, Borrell M, Fontcuberta J. Association after linkage analysis indicates that homozygosity for the 46C>T polymorphism in the F12 gene is a genetic risk factor for venous thrombosis. Thromb Haemost. 2004;91:899-904

6. Reuner KH, Jenetzky E, Aleu A, Litfin F, Mellado P, Kloss M, Jüttler E, Grau AJ, Rickmann H, Patscheke H, Lichy C. Factor XII c46T genepolymorphism and the risk of cerebral venous thrombosis. Neurology. 2008;70:129-132.

7. Cochery-Nouvellon E, Mercier E, Lissalde-Lavigne G, Daurès JP, Quéré I, Dauzat M, Marès P, Gris JC.
Homozygosity for the C46T polymorphism of the F12 gene is a risk factor for venous thrombosis during the first
pregnancy. J Thromb Haemost. 2007;5:700-707.

8. Catto AJ, Kohler HP, Coore J, Mansfield MW, Stickland MH, Grant PJ. Association of a common polymorphism in the factor XIII gene with venous thrombosis. Blood. 1999;93:906-908.

9. Wells PS, Anderson JL, Scarvelis DK, Doucette SP, Gagnon F. Factor XIII Val34Leu variant is protective against venous thromboembolism: a HuGE review and meta-analysis. Am J Epidemiol. 2006;164:101-109.

10. Van HylckamaVlieg A, Komansin N, Ariëns RAS, Poort SR, Grant PJ, Bertina RM, Rossendaal FR. Factor XIII Val34Leu polymorphism, factor XIII antigen levels and activity and the risk of deep venous thrombosis. BJH. 2002;119:169-175.

11. Corral J, Gonzalez-Conejero R, Soria JM, Gonzalez-Porras JR, Pérez-Ceballos E, Lecumberri R, Roldán V, Souto JC, Miñano A, Hernández-Espinosa D, Alberca I, Fontcuberta J, Vicente V. A nonsense polymorphism in the protein Z-dependent protease inhibitor increases the risk for venous thrombosis. Blood. 2006;108:177-183.

12. Van de Water N, Tan T, Ashton F, O'Grady A, Day T, Browett P, Ockelford P, Harper P. Mutations within the protein Z-dependent protease inhibitor gene are associated with venous thromboembolic disease: a new form of thrombophilia. Br J Haematol. 2004;127:190-194.

13. Corral J, Hernandez-Espinosa D, Soria JM. Antithrombin Cambridge II (A384S): an underestimated genetic risk factor for venous thrombosis. Blood. 2007;109:4258-4263.

14. Austin H, De Staercke C, Lally C, Bezemer ID, Rosendaal FR, Hooper WC. New gene variants associated with venous thrombosis: a replication study in White and black Americans. J Thromb Haemost. 2011;9:489-495.

15. Bezemer ID, Bare LA, Doggen CJ, Arellano AR, Tong C, Rowland CM, Catanese J, Young BA, Reitsma PH,
Devlin JJ, Rosendaal FR. Gene variants associated with deep vein thrombosis. JAMA. 2008;299:1306-1314.

16. De Haan HG, Bezemer ID, Doggen CJM, Le Cessie S, Reitsma PH, Arellano AR, Tong CH, Devlin JJ, Bare LA, Rosendaal FR, Vossen CY. Multiple SNP testing improves risk prediction of first venous thrombosis. Blood. 2012;120:656-663.

17. Uitte de Willige S, de Visser MC, Houwing-Duistermaat JJ, Rosendaal, FR, Vos HL, Bertina RM. Genetic variation in the fibrinogen gamma gene increases the risk of deep venous thrombosis by reducing plasma fibrinogen gamma levels. Blood. 2005;106:4176-7183.

18. Grünbacher G, Weger W, Max-Neuhold E, Pilger E, Köppel H, Wascher T, März W, Renner W. The fibrinogen gamma (FGG) 10034C>T polymorphism is associated with venous thrombosis. Thromb Res. 2007;121:33-36.

19. Morange PE, Oudot-Mellakh T, Cohen W, Germain M, Saut N, Antoni G, Alessi MC, Bertrand M, Dupuy AM, Letenneur L, Lathrop M, Lopez LM, Lambert JC, Emmerich J, Amouyel P, trégoüet DA. KNG1 ile581Thr and susceptibility to venous thrombosis. Blood. 2011;117:3692-3694.

20. Arellano AR, Bezemer ID, Tong CH, Catanese JJ, Devlin JJ, Reitsma PH, Bare LA, Rosendaal FR. Gene variants associated with venous thrombosis: confirmation in the MEGA study. J Thromb Haemost. 2010;8:1132-1134.

21. Austin H, De Staercke C, Lally C, Bezemer ID, Rosendaal FR, Hooper WC. New gene variants associated with venous thrombosis: a replication study in White and black Americans. J Thromb Haemost. 2011;9:489-495.
<figure><img src="image_2.png" /></figure>

# Multilocus Genetic Risk Scores for Venous Thromboembolism Risk Assessment

José Manuel Soria, Pierre-Emmanuel Morange, Joan Vila, Juan Carlos Souto, Manel Moyano, David-Alexandre Trégouët, José Mateo, Noémi Saut, Eduardo Salas and Roberto Elosua

J Am Heart Assoc. 2014;3:e001060; originally published October 23, 2014; doi: 10.1161/JAHA.114.001060

The Journal of the American Heart Association is published by the American Heart Association, 7272 Greenville Avenue, Dallas, TX 75231

Online ISSN: 2047-9980

The online version of this article, along with updated information and services, is located on the World Wide Web at:

http://jaha.ahajournals.org/content/3/5/e001060

Data Supplement (unedited) at:

http://jaha.ahajournals.org/content/suppl/2014/10/23/jah3707.DC1.html
# SUPPLEMENTAL MATERIAL

## MULTILOCUS GENETIC RISK SCORES FOR VENOUS THROMBOEMBOLISM RISK ASSESSMENT

Soria et al

<table><caption>**Table 1.** Genetic variants included in the different genetic risk scores (GRS) assessed: coefficients (weights) assigned to each genetic risk factor, and allele frequency and Hardy-Weinberg equilibrium (H-W) data obtained in the two studies, Sant Pau and MARTHA.</caption><thead><tr><th rowspan="2">SNP</th><th rowspan="2">Gene</th><th rowspan="2">Mutation</th><th rowspan="2">Genetic risk group*</th><th rowspan="2">Risk coefficient assigned</th><th colspan="4">GRS</th><th colspan="2">SANT PAU</th><th colspan="2">MARTHA</th></tr><tr><th>1</th><th>2</th><th>3</th><th>4</th><th>Rare allele Frequency</th><th>H-W</th><th>Rare allele frequency</th><th>H-W</th></tr></thead><tbody><tr><td>rs6025</td><td>F5 Leiden</td><td>R506</td><td>A carriers</td><td>1.589</td><td></td><td></td><td></td><td></td><td>0.01</td><td>0.873</td><td>0.11</td><td>0.008</td></tr><tr><td>rs118203905</td><td>F5 Hong Kong</td><td>F5 Hong Kong</td><td>A carriers</td><td>1.589</td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.00</td><td>0.982</td></tr><tr><td>rs118203906</td><td>F5 Cambridge</td><td>F5 Cambridge</td><td>C carriers</td><td>1.589</td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.00</td><td>---</td></tr><tr><td>rs1799963</td><td>F2</td><td>G20210A</td><td>A carriers</td><td>0.293</td><td></td><td></td><td></td><td></td><td>0.01</td><td>0.822</td><td>0.10</td><td>0.020</td></tr><tr><td>ABO**</td><td>ABO</td><td></td><td>A1 carriers</td><td>0.956</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>rs8176719</td><td>ABO</td><td></td><td></td><td>--</td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.36</td><td>0.262</td></tr><tr><td>rs1801020</td><td>F12</td><td>C46T</td><td>TT homoz.</td><td>1.633</td><td></td><td></td><td></td><td></td><td>0.21</td><td>0.033</td><td>0.23</td><td>0.223</td></tr><tr><td>rs5985</td><td>F13</td><td>V34L</td><td>GG homoz.</td><td>0.198</td><td></td><td></td><td></td><td></td><td>0.25</td><td>0.832</td><td>0.27</td><td>0.929</td></tr><tr><td>rs2232698</td><td>SERPINE10</td><td>R67X</td><td>T carriers</td><td>1.358</td><td></td><td></td><td></td><td></td><td>0.01</td><td>0.898</td><td>0.01</td><td>0.854</td></tr><tr><td>rs121909548</td><td>SERPINC1</td><td>A384S</td><td>T carriers</td><td>2.277</td><td></td><td></td><td></td><td></td><td>0.00</td><td>0.975</td><td>0.00</td><td>0.945</td></tr><tr><td rowspan="2">rs2036914</td><td rowspan="2">FXI</td><td rowspan="2"></td><td>CT heter.</td><td>0.293</td><td></td><td></td><td></td><td></td><td rowspan="2">0.46</td><td rowspan="2">0.293</td><td rowspan="2">0.50</td><td rowspan="2">0.890</td></tr><tr><td>TT homoz.</td><td>0.519</td><td></td><td></td><td></td><td></td></tr><tr><td>rs2066865</td><td>FGG</td><td></td><td>TT homoz.</td><td>0.344</td><td></td><td></td><td></td><td></td><td>0.21</td><td>0.739</td><td>0.23</td><td>0.028</td></tr><tr><td>rs710446</td><td>KNG1</td><td></td><td>CC homoz.</td><td>0.182</td><td></td><td></td><td></td><td></td><td>0.44</td><td>0.831</td><td>0.42</td><td>0.131</td></tr><tr><td rowspan="2">rs2289252</td><td rowspan="2">FXI</td><td rowspan="2"></td><td>CT heter.</td><td>0.315</td><td></td><td></td><td></td><td></td><td rowspan="2">0.40</td><td rowspan="2">0.788</td><td rowspan="2">0.40</td><td rowspan="2">0.655</td></tr><tr><td>TT homoz.</td><td>0.577</td><td></td><td></td><td></td><td></td></tr></tbody></table>

<footnote>*Homoz.=Homozygous; Heter.=Heterozygous. ** ABO: A1 carriers (rs8176719, rs7853989, rs8176743, rs8176750)</footnote>
<table><caption>**Table 2:** Data sources and results of the meta-analyses undertaken by the authors to assign a risk coefficient to each variant included in the genetic risk score.</caption><thead><tr><th></th><th>Gen</th><th>Reference</th><th>Number of cases</th><th>Number of controls</th><th>OR</th><th>Lower confidence interval</th><th>Upper confidence interval</th><th>Risk Coefficient assigned</th></tr></thead><tbody><tr><td>rs6025-FVL</td><td>F5</td><td>1</td><td></td><td></td><td></td><td></td><td></td><td>1.589</td></tr><tr><td>rs118203905-F5 Hong-Kong</td><td>F5</td><td></td><td></td><td></td><td></td><td></td><td></td><td>1.589</td></tr><tr><td>rs118203906-F5 Cambridge</td><td>F5</td><td></td><td></td><td></td><td></td><td></td><td></td><td>1.589</td></tr><tr><td>rs1799963</td><td>F2</td><td>2</td><td></td><td></td><td></td><td></td><td></td><td>0.293</td></tr><tr><td>AB0</td><td>AB0</td><td>3</td><td></td><td></td><td></td><td></td><td></td><td>0.956</td></tr><tr><td>rs8176719</td><td>AB0</td><td>3,4</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>rs1801020</td><td>F12</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>5</td><td>250</td><td>250</td><td>4.82</td><td>1.50</td><td>15.60</td><td></td></tr><tr><td></td><td></td><td>6</td><td>--</td><td>--</td><td>4.57</td><td>1.55</td><td>13.40</td><td></td></tr><tr><td></td><td></td><td>7</td><td>32463</td><td>32463</td><td>5.99</td><td>2.10</td><td>17.30</td><td></td></tr><tr><td>Meta-analysis</td><td></td><td></td><td></td><td></td><td>5.12</td><td></td><td></td><td>1.633</td></tr><tr><td>rs5985</td><td>F13</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>8</td><td>221</td><td>254</td><td>0.63</td><td>0.38</td><td>0.82</td><td></td></tr><tr><td></td><td></td><td>9</td><td>3165</td><td>4909</td><td>0.85</td><td>0.77</td><td>0.95</td><td></td></tr><tr><td></td><td></td><td>10</td><td>475</td><td>475</td><td>0.80</td><td>0.60</td><td>1.00</td><td></td></tr><tr><td>Meta-analysis</td><td></td><td></td><td></td><td></td><td>0.82</td><td></td><td></td><td>-0.198 (0.198)*</td></tr><tr><td>rs2232698</td><td>SERPINE10</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>11</td><td>1018</td><td>1018</td><td>3.35</td><td>1.30</td><td>8.60</td><td></td></tr><tr><td></td><td></td><td>12</td><td>250</td><td>250</td><td>5.70</td><td>1.25</td><td>26.0</td><td></td></tr><tr><td>Meta-analysis</td><td></td><td></td><td></td><td></td><td>3.89</td><td></td><td></td><td>1.358</td></tr><tr><td>rs121909548</td><td>SERPINC1</td><td>13</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>rs2036914</td><td>F11</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Bezemer (2008) : Heteroz. Homoz.</td><td></td><td>14</td><td>1314</td><td>2877</td><td>1.38<br>1.71</td><td>1.17<br>1.43</td><td>1.64<br>2.05</td><td></td></tr><tr><td>Austin (2011) : Heteroz. Homoz.</td><td></td><td>15</td><td>1076</td><td>1239</td><td>1.20<br>1.60</td><td>0.92<br>1.20</td><td>1.70<br>2.30</td><td></td></tr><tr><td>Meta-analysis: Heteroz. Homoz.</td><td></td><td></td><td></td><td></td><td>1.34<br>1.68</td><td></td><td></td><td>0.293<br>0.519</td></tr><tr><td>rs2066865</td><td>FGG</td><td>16, 17, 18</td><td></td><td></td><td></td><td></td><td></td><td>0.344</td></tr><tr><td>rs710446</td><td>KNG1</td><td>19</td><td></td><td></td><td></td><td></td><td></td><td>0.182</td></tr><tr><td>rs2289252</td><td>F11</td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.315</td></tr></tbody></table>
<table><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.577</td></tr><tr><td>Arellano (2010): Heteroz.<br/>Homoz.</td><td></td><td>20</td><td>3921</td><td>4634</td><td>1.37<br/>1.78</td><td>1.24<br/>1.57</td><td>1.51<br/>2.01</td><td></td></tr><tr><td>Austin (2011): Heteroz.<br/>Homoz.</td><td></td><td>21</td><td>1076</td><td>1239</td><td>1.40<br/>1.80</td><td>1.10<br/>1.30</td><td>1.80<br/>2.50</td><td></td></tr><tr><td>Meta-analysis: Heteroz.<br/>Homoz.</td><td></td><td></td><td></td><td></td><td>1.37<br/>1.78</td><td></td><td></td><td>0.315<br/>0.577</td></tr></table>

<footnote>*In the case of the variant rs5985, the rare variant was associated with lower probability of presenting a VTE. To build the GRS the common variant was considered to be associated with a higher probability of presenting of VTE.</footnote>
<table><caption>**Table 3.** Observed and expected strengths of associations between selected clinical or genetic factors and venous thromboembolic events and p-values for the differences between observed and expected values.</caption><thead><tr><th rowspan="3"></th><th rowspan="3">Expected<br>coefficients</th><th colspan="3">SANT PAU</th><th colspan="3">MARTHA</th></tr><tr><th colspan="3">Observed</th><th colspan="3">Observed</th></tr><tr><th>Coefficients</th><th>Standard<br>Error</th><th>p-value*</th><th>Coefficients</th><th>Standard<br>Error</th><th>p-value*</th></tr></thead><tbody><tr><td>Family History</td><td>1.185</td><td>0.711</td><td>0.240</td><td>0.323</td><td>NA</td><td>NA</td><td>NA</td></tr><tr><td>F5**†</td><td>1.589</td><td>2.803</td><td>0.755</td><td>0.421</td><td>0.805</td><td>0.178</td><td>0.028</td></tr><tr><td>F2- rs1799963†</td><td>0.293</td><td>0.773</td><td>0.489</td><td>0.624</td><td>0.378</td><td>0.199</td><td>0.831</td></tr><tr><td>ABO (A1)†</td><td>0.956</td><td>0.997</td><td>0.224</td><td>0.927</td><td>0.428</td><td>0.280</td><td>0.346</td></tr><tr><td>F12- rs1801020†</td><td>1.633</td><td>0.859</td><td>0.621</td><td>0.533</td><td>0.305</td><td>0.345</td><td>0.054</td></tr><tr><td>F13- rs5985†</td><td>0.198</td><td>-0.078</td><td>0.224</td><td>0.538</td><td>0.171</td><td>0.152</td><td>0.929</td></tr><tr><td>SERPINE10- rs2232698†</td><td>1.358</td><td>0.937</td><td>0.753</td><td>0.780</td><td>0.572</td><td>0.502</td><td>0.434</td></tr><tr><td>SERPINC1- rs121909548†</td><td>2.277</td><td>1.476</td><td>1.174</td><td>0.733</td><td>-0.983</td><td>1.264</td><td>0.197</td></tr><tr><td>F11- rs2036914 Hetero</td><td>0.293</td><td>-0.039</td><td>0.273</td><td>0.543</td><td>0.228</td><td>0.207</td><td>0.875</td></tr><tr><td>Homo</td><td>0.519</td><td>0.028</td><td>0.395</td><td>0.534</td><td>0.206</td><td>0.258</td><td>0.544</td></tr></tbody></table>
<table><tr><td>FGG- rs2066865 †</td><td>0.344</td><td>-0.095</td><td>0.227</td><td>0.334</td><td>0.328</td><td>0.151</td><td>0.958</td></tr><tr><td>KNG1- rs710446b ‡</td><td>0.182</td><td>-0.206</td><td>0.278</td><td>0.485</td><td>-0.086</td><td>0.190</td><td>0.481</td></tr><tr><td>F11- rs2289252 Hetero</td><td>0.315</td><td>0.186</td><td>0.299</td><td>0.829</td><td>0.191</td><td>0.190</td><td>0.744</td></tr><tr><td>Homo</td><td>0.577</td><td>-0.017</td><td>0.403</td><td>0.461</td><td>0.965</td><td>0.262</td><td>0.459</td></tr></table>

*for the difference between observed and expected coefficients.

** Carrier of any risk allele (Leiden, Hong-Kong or Cambridge).

NA indicates not applicable.
References:

1. Juul K, Tybjaerg-Hansen A, Schnohr P, Nordestgaard BG. Factor V Leiden and the risk for venous thromboembolism in the adult Danish population. Ann Intern Med. 2004;140:330-337.

2. Wu O, Robertson L, Twaddle S, Lowe G, Clark P, Walker I, Brenkel I, Greaves M, Langhorne P, Regan L, Greer I; Thrombosis: Risk and Economic Assessment of Thrombophilia Screening (TREATS) Study. Screening for thrombophilia in high-risk situations: a meta-analysis and cost-effectiveness analysis. Br J Haematol. 2005;131:80-90.

3. Tirado I, Mateo J, Soria JM, Oliver A, Martínez-Sánchez E, Vallvé C, Borrell M, Urrutia T, Fontcuberta J. The ABO blood group genotype and factor VIII levels as independent risk factors for venous thromboembolism. Thromb Haemost. 2005;93:468-474.

4. De Haan HG, Bezemer ID, Doggen CJM, Le Cessie S, Reitsma PH, Arellano AR, Tong CH, Devlin JJ, Bare LA, Rosendaal FR, Vossen CY. Multiple SNP testing improves risk prediction of first venous thrombosis. Blood. 2012;120:656-663.

5. Tirado I, Soria JM, Mateo J, Oliver A, Souto JC, Santamaría A, Felices R, Borrell M, Fontcuberta J. Association after linkage analysis indicates that homozygosity for the 46C>T polymorphism in the F12 gene is a genetic risk factor for venous thrombosis. Thromb Haemost. 2004;91:899-904

6. Reuner KH, Jenetzky E, Aleu A, Litfin F, Mellado P, Kloss M, Jüttler E, Grau AJ, Rickmann H, Patscheke H, Lichy C. Factor XII c46T genepolymorphism and the risk of cerebral venous thrombosis. Neurology. 2008;70:129-132.

7. Cochery-Nouvellon E, Mercier E, Lissalde-Lavigne G, Daurès JP, Quéré I, Dauzat M, Marès P, Gris JC. Homozygosity for the C46T polymorphism of the F12 gene is a risk factor for venous thrombosis during the first
pregnancy. J Thromb Haemost. 2007;5:700-707.

8. Catto AJ, Kohler HP, Coore J, Mansfield MW, Stickland MH, Grant PJ. Association of a common polymorphism in the factor XIII gene with venous thrombosis. Blood. 1999;93:906-908.

9. Wells PS, Anderson JL, Scarvelis DK, Doucette SP, Gagnon F. Factor XIII Val34Leu variant is protective against venous thromboembolism: a HuGE review and meta-analysis. Am J Epidemiol. 2006;164:101-109.

10. Van HylckamaVlieg A, Komansin N, Ariëns RAS, Poort SR, Grant PJ, Bertina RM, Rossendaal FR. Factor XIII Val34Leu polymorphism, factor XIII antigen levels and activity and the risk of deep venous thrombosis. BJH. 2002;119:169-175.

11. Corral J, Gonzalez-Conejero R, Soria JM, Gonzalez-Porras JR, Pérez-Ceballos E, Lecumberri R, Roldán V, Souto JC, Miñano A, Hernández-Espinosa D, Alberca I, Fontcuberta J, Vicente V. A nonsense polymorphism in the protein Z-dependent protease inhibitor increases the risk for venous thrombosis. Blood. 2006;108:177-183.

12. Van de Water N, Tan T, Ashton F, O'Grady A, Day T, Browett P, Ockelford P, Harper P. Mutations within the protein Z-dependent protease inhibitor gene are associated with venous thromboembolic disease: a new form of thrombophilia. Br J Haematol. 2004;127:190-194.

13. Corral J, Hernandez-Espinosa D, Soria JM. Antithrombin Cambridge II (A384S): an underestimated genetic risk factor for venous thrombosis. Blood. 2007;109:4258-4263.

14. Austin H, De Staercke C, Lally C, Bezemer ID, Rosendaal FR, Hooper WC. New gene variants associated with venous thrombosis: a replication study in White and black Americans. J Thromb Haemost. 2011;9:489-495.

15. Bezemer ID, Bare LA, Doggen CJ, Arellano AR, Tong C, Rowland CM, Catanese J, Young BA, Reitsma PH,
Devlin JJ, Rosendaal FR. Gene variants associated with deep vein thrombosis. JAMA. 2008;299:1306-1314.

16. De Haan HG, Bezemer ID, Doggen CJM, Le Cessie S, Reitsma PH, Arellano AR, Tong CH, Devlin JJ, Bare LA, Rosendaal FR, Vossen CY. Multiple SNP testing improves risk prediction of first venous thrombosis. Blood. 2012;120:656-663.

17. Uitte de Willige S, de Visser MC, Houwing-Duistermaat JJ, Rosendaal, FR, Vos HL, Bertina RM. Genetic variation in the fibrinogen gamma gene increases the risk of deep venous thrombosis by reducing plasma fibrinogen gamma levels. Blood. 2005;106:4176-7183.

18. Grünbacher G, Weger W, Max-Neuhold E, Pilger E, Köppel H, Wascher T, März W, Renner W. The fibrinogen gamma (FGG) 10034C>T polymorphism is associated with venous thrombosis. Thromb Res. 2007;121:33-36.

19. Morange PE, Oudot-Mellakh T, Cohen W, Germain M, Saut N, Antoni G, Alessi MC, Bertrand M, Dupuy AM, Letenneur L, Lathrop M, Lopez LM, Lambert JC, Emmerich J, Amouyel P, trégoüet DA. KNG1 ile581Thr and susceptibility to venous thrombosis. Blood. 2011;117:3692-3694.

20. Arellano AR, Bezemer ID, Tong CH, Catanese JJ, Devlin JJ, Reitsma PH, Bare LA, Rosendaal FR. Gene variants associated with venous thrombosis: confirmation in the MEGA study. J Thromb Haemost. 2010;8:1132-1134.

21. Austin H, De Staercke C, Lally C, Bezemer ID, Rosendaal FR, Hooper WC. New gene variants associated with venous thrombosis: a replication study in White and black Americans. J Thromb Haemost. 2011;9:489-495.