Abstract
Although chronic thromboembolic pulmonary hypertension (CTEPH) is characterised by the persistence of organised thrombus, few pro-thrombotic risk factors have been identified in subjects with the disease. The aim of the present study was to compare the prevalence of eight functionally relevant haemostatic polymorphisms between CTEPH subjects and healthy controls.
Genomic DNA was isolated from 214 CTEPH subjects and 200 healthy controls, and analysed for Factor V Leiden, prothrombin guanine (G) to adenine (A) substitution at nucleotide 20210 (20210G>A), plasminogen activator inhibitor-1 4G/5G, tissue plasminogen activator 7351 cytosine (C)>thymidine (T), Factor XIII 100G>T, fibrinogen Aα substitution of threonine with alanine at position 312 (Thr312Ala), fibrinogen Bβ substitution of arginine with lysine at position 448 (Arg448Lys) and fibrinogen Bβ 455G>A polymorphisms.
A significant difference was demonstrated in fibrinogen Aα Thr312Ala genotype and allele frequencies between CTEPH subjects and controls. The presence of the alanine allele significantly increased the risk of CTEPH.
The fibrinogen Aα alanine 312 allele alters fibrinogen α–α chain cross-linkage and has previously been associated with both increased risk of embolisation and increased resistance to thrombolysis. An association between this polymorphism and chronic thromboembolic pulmonary hypertension, therefore, supports an embolic aetiology for this disease, and may provide a mechanism by which thrombus persists following an acute event.
Chronic thromboembolic pulmonary hypertension (CTEPH) is increasingly being recognised as an important cause of pulmonary hypertension, and is generally considered to represent an uncommon consequence of acute pulmonary embolic (PE) disease 1. CTEPH has been reported to occur after ∼4% of cases of acute PE 2, whereas retrospective studies have described a preceding deep venous thrombosis (DVT) in approximately a third of cases 3, 4. Therefore, it is presumed that, following an acute embolic event, thrombus persists and organises, resulting in chronic obstruction to pulmonary arterial blood flow.
Given the key role of thrombus formation in CTEPH, it would intuitively be expected that there be some aberration of haemostasis, at either a local or systemic level. However, although hereditary thrombophilic conditions are commonly associated with acute PE and DVT 5, 6, no such conditions have been strongly linked with CTEPH 7. Elevated serum levels of antiphospholipids and Factor VIII have been described previously 8, 9, but it is unclear whether these represent primary or secondary phenomena. As a result, no pro-thrombotic pathogenic mechanisms have been clearly shown to predispose to CTEPH. It is possible, however, that other steps in the haemostatic process, such as fibrin stabilisation or fibrinolysis, may be involved in the pathophysiology of the disease 10. A number of single nucleotide polymorphisms (SNPs) have been described that predispose to arterial and venous thrombotic disease, by affecting these more global aspects of haemostasis 11–20. This case–control study, therefore, uses a candidate gene approach to determine the potential relevance of these SNPs in CTEPH.
METHODS
Subjects
All consecutive Caucasian patients attending Papworth Hospital (Papworth Everard, UK) with a diagnosis of CTEPH during 1999–2006 were offered entry to the present study. CTEPH was diagnosed through clinical evaluation, right heart catheterisation and appropriate imaging (chest radiography, ventilation/perfusion scintigraphy, computed tomographic pulmonary angiography, magnetic resonance pulmonary angiography and catheter-directed angiography), following a standardised protocol described elsewhere 4, 21. CTEPH subjects (n = 214) included those with proximal (n = 169) and distal disease (n = 45). Healthy controls (100 males and 100 females) were recruited from a National Blood Service site in West London (UK) and comprised Caucasians aged 18–65 yrs.
The study protocol was approved by the Huntingdon Research Ethics Committee (UK). Written informed consent was obtained from each subject before entry into the study.
Genotyping
Genomic DNA was extracted from whole blood (Nucleon extraction kit; Amersham Biosciences, Little Chalfont, UK) and analysed for SNPs using commercially available primers (Taqman® SNP genotyping assay; Applied Biosystems, Foster City, CA, USA) using an ABI Prism® 7900HT Sequence Detection system (Applied Biosystems). Contract research services for genotyping were provided by Geneservice (Cambridge, UK). The eight SNPs studied were common genomic variations of the following haemostatic factors: Factor V (Factor V Leiden mutation or guanine (G) to adenine (A) substitution at nucleotide 1691 (1691G>A); National Center for Biotechnology Information (Bethesda, MD, USA) SNP database accession code rs6025) 11; prothrombin (20210G>A; rs 1799963) 11; plasminogen activator inhibitor-1 (4G/5G; rs1799768) 11; tissue plasminogen activator (cytosine (C) to thymidine (T) substitution at nucleotide 7351 (7351C>T; rs2020918)) 17; Factor XIII (100G>T; rs5985) 15; fibrinogen Aα Thr312Ala (substitution of threonine (Thr) at position 312 with alanine (Ala); rs6050); and fibrinogen Bβ Arg448Lys (substitution of arginine (Arg) with lysine (Lys) at position 448; rs4220); and fibrinogen Bβ G455A (rs1800790) 22.
Statistics
Genotype and allele frequencies were calculated for each locus. Observed frequencies in controls were compared with those predicted by the Hardy–Weinberg equilibrium equation using the Chi-squared test. Linkage disequilibrium (LD) between SNPs was calculated from maximum likelihood haplotype frequencies estimated using the expectation maximisation algorithm 23.
Genotype and allele frequencies were compared between CTEPH and control subjects using the Chi-squared test, or Fisher’s exact test in the case of low numbers. The risk of disease associated with each genotype (genotype model) was estimated by logistic regression analysis, where SNP genotype was coded in a three-level fashion, i.e. rarest homozygous genotype, heterozygous genotype and commonest homozygous genotype, with the latter being the reference category. For SNPs for which there were no rare homozygous subjects, the genotype was coded as heterozygous and compared with the most common homozygous genotype. The nature of the relationship between each SNP and risk of disease was also examined using models of dominant and recessive inheritance (allele model). In the dominant case, having one or more copies of the rare allele was compared with having no copies. In the recessive case, having two copies of the rare allele was compared with having one or no copies. Multiple logistic regression models of the effect of pairs of SNPs were used to assess whether or not there was any synergistic effect of having two SNPs. In view of sex differences between the two groups, all analyses were also repeated using multiple logistic regression models that adjusted for sex.
Adjustment for multiple tests was not performed for single SNP analyses, since each has previously been shown to be independently associated with an increased risk of thrombosis, and the present study aimed to re-address these hypotheses in the context of CTEPH. As such, calculations for each SNP constituted a separate and distinct analysis.
RESULTS
Subjects and samples
The baseline demographic characteristics of the CTEPH subjects are shown in table 1⇓. Only sex information was available for the 200 healthy controls, although it is known that they were Caucasian, aged 18–65 yrs and fulfilled the criteria required for blood donation.
Baseline demographics of chronic thromboembolic pulmonary hypertension subjects
The mean yield from SNP analyses was 96.3%. All SNPs were in Hardy–Weinberg equilibrium, except for fibrinogen Bβ Arg448Lys (p = 0.045). Only Fibrinogen Bβ Arg448Lys and 455G>A were demonstrated to be in LD (r2 = 0.91).
The genotype and allele frequencies of each of the SNPs in the two groups are shown in table 2⇓.
Genotype and allele frequencies of the eight candidate genes
Polymorphism frequencies
A significant difference was demonstrated in fibrinogen Aα Thr312Ala genotype (p = 0.03) and allele (p = 0.01) frequencies between the CTEPH and control groups (table 2⇑). The three-level genotype model showed an association between the fibrinogen Aα Thr312Ala SNP and CTEPH, although this did not reach significance for homozygotes (odds ratio (OR) 1.65, 95% confidence interval (CI) 1.13–2.49 for heterozygotes, p = 0.02; OR 1.82, 95% CI 0.92–3.61 for homozygotes, p = 0.09). A dominant model of inheritance of the Ala allele (i.e. one or two copies of the Ala allele versus none) showed increased odds of CTEPH (OR 1.68, 95%CI 1.13–2.49, p = 0.01). There was no evidence of confounding of the relationship between the fibrinogen Aα Thr312Ala SNP and CTEPH by sex (dominant model (females versus males): OR 1.06, 95% CI 0.72–1.57, p = 0.77; three-level genotype model (females versus males); OR 1.06, 95% CI 0.72–1.57, p = 0.78)
There was a difference in Factor V Leiden genotype frequencies between the two groups, although this did not reach significance (p = 0.051). The Factor V Leiden mutation was associated with an increased risk of CTEPH (genotype model (no rare homozygotes): OR 2.54, 95% CI 0.97–6.68, p = 0.06), although, again, this was not significant.
There were no other significant differences between the two groups and no evidence of synergism between SNP pairs.
The six SNPs that failed to show an association with CTEPH were also used as a form of genomic control, to test for population stratification. The Kolmogorov–Smirnov goodness-of-fit test of the Chi-squared statistics for the association between each SNP and disease demonstrated that these followed a distribution that was not significantly different from a Chi-squared distribution (p = 0.35). These results suggest a lack of stratification.
DISCUSSION
The present study demonstrates that the presence of the fibrinogen Aα chain Ala312 allele is associated with increased odds of being diagnosed with CTEPH. In addition, a significant association with the Factor V Leiden mutation was also noted. These are the first haemostatic polymorphisms to be reported in a large series of CTEPH subjects, and thus may add to understanding of the pathophysiological mechanisms behind this enigmatic disease.
Although CTEPH is widely considered to be a consequence of acute PE disease 1, 24, the evidence supporting this belief is largely circumstantial. One concern that has fuelled the debate over the aetiology of this disease is the discrepancy that exists between the thrombophilic profiles of DVT patients compared with CTEPH patients 25. However, rather than viewing this as a barrier to the embolic theory of CTEPH, it is worth considering what might distinguish the minority of patients who develop CTEPH following a simple DVT from the vast majority who do not. In this context, it is important to consider not only the haemostatic factors that predispose to thrombus formation but also those that determine how thrombus is subsequently handled. A key factor in this process is clearly the fibrinogen molecule, which provides the essential building blocks for the formation of the fibrin clot. Abnormalities within this molecule may not only promote thrombus formation but also delay its degradation, thus encouraging the switch from thrombus resolution to thrombus organisation.
The fibrinogen molecule comprises three pairs of polypeptide chains (Aα, Bβ and γ), which are linked by disulphide bonds and arranged into three regions: a central (E) region, containing fibrinopeptides A and B and the amino acid termini of all six chains; and two distal (D) regions, each formed by the unpaired carboxyl termini of the Aα, Bβ and γ chains 26. D regions also include a globular Cα domain, formed by the terminal two-thirds of the α-chain. Within a given fibrinogen molecule, the two Cα domains are tethered together in close proximity to the E region 27. Following activation of the coagulation cascade, thrombin cleaves fibrinopeptides A and B from fibrinogen, producing fibrin monomers and initiating clot formation. Thrombin simultaneously activates Factor XIII, promoting cross-linkage between newly formed fibrin oligomers. As part of this process, Cα domains, now untethered from the E region following cleavage of fibrinopeptide B, are encouraged to bind with Cα domains from neighbouring fibrin molecules 28. The resulting α–α linkage strengthens the association between adjacent protofibrils and stabilises the developing fibrin structure. This permits the formation of a fibrin structure that is stronger, more rigid and more resistant to thrombolysis 29.
The pivotal role of the Cα domain in fibrin polymerisation has previously been highlighted by an in vitro study that demonstrated that monoclonal antibodies directed against this domain inhibit the formation of fibrin polymers from monomers 30. Recombinant fibrinogen lacking the Cα domain has also been shown to form clot in vitro that is less stiff and more susceptible to thrombolysis than clot formed in the presence of normal fibrinogen 31. Further information about the role of the Cα domain has also been gained from studying hereditary dysfibrinogenaemias 32. Dusart syndrome, caused by the substitution of Arg fibrinogen Aα for cysteine at position 554 within, is associated with increased α–α chain linkage, and hence enhanced lateral aggregation of fibrin protofibrils. Fibrin clots from patients with this disease tend to be stiff, impermeable and resistant to thrombolysis in vitro 33, 34. Clinically, patients suffer recurrent thrombotic complications and are particularly prone to embolic phenomena 35. Conversely, the dysfibrinogen Caracas II, caused by the substitution of serine for asparagines at position 434 within fibrinogen Aα, is associated with reduced α–α chain linkage compared with normal fibrinogen. Clots are characteristically loosely packed in structure, and show normal stiffness but increased permeability 36. Patients with the disorder are typically asymptomatic and do not suffer any haematological complications 36.
Thr312Ala is a well-recognised polymorphism of the Aα gene, and encodes a Thr-to-Ala substitution at residue 312 of the Cα domain 37. Although the functional relevance of this polymorphism remains unclear, the area immediately surrounding this site is known to be central to both Factor XIII activation and to Factor XIII-dependent α-chain cross-linking and α2-antiplasmin binding 38, 39. Clinically, homozygosity for the Ala allele has been associated with acute PE disease but not DVT 14. The presence of the Ala allele, in both heterozygotes and homozygotes, has also been associated with higher post-stroke mortality in subjects with atrial fibrillation but not those in sinus rhythm 13. It has, therefore, been suggested that the Ala substitution may predispose to embolic phenomena, through the formation of more brittle clot 40. In the case of DVT, such clot would be more likely to fragment and dislodge as it propagates proximally and comes under the influence of oblique flow from conjoining veins. The results of the present study support this hypothesis, and contribute further to the evidence that CTEPH is of embolic aetiology.
Evidence from in vitro studies suggests that, other than simply predisposing patients to embolic phenomena, the Aα Thr312Ala polymorphism may play an additional pathophysiological role in CTEPH. Clot formed from Ala/Ala subjects demonstrates more extensive α-chain cross-linking, increased fibre thickness and greater clot stiffness than is seen in Thr/Thr controls 40. Such clot also shows less permeability to flow, which, in turn, permits less access to fibrinolytic factors in vivo 41. These properties suggest that the Ala substitution enhances lateral aggregation of fibrin protofibrils and encourages formation of a more tightly packed fibrin structure. Although not specifically measured in these studies, such parameters are typically associated with clot that is more resistant to thrombolysis. Moreover, recent data suggest that fibrin derived from CTEPH patients is slower to lyse than that derived from controls 10. This combined evidence implies that delayed fibrin degradation, possibly mediated in some by the Aα Thr312Ala polymorphism, is a pivotal step in the development of organised scar tissue from acute PE material in a subset of patients with CTEPH.
The present study also suggested a difference in the frequency of the Factor V Leiden mutation between CTEPH subjects and controls, although this difference did not reach significance (p = 0.051). The Factor V Leiden mutation, characterised by substitution of arginine for glutamic acid at position 506 is resistant to cleavage by activated protein C (APC) and is, therefore, inactivated more slowly 42. Additionally, the presence of this mutation interferes with the cofactor activity of Factor V on APC-catalysed inactivation of Factor VIIIa 43. Both mechanisms promote a hypercoagulable state that has been significantly associated with the risk of developing DVT 44. Previous small studies have shown that the prevalence of this CTEPH mutation is low and not dissimilar to that of controls 9, 45. The present study examines a much larger cohort of patients and suggests a possible association, although the prevalence remains much lower than that found in patients with DVT. Given these data, a role for Factor V Leiden in CTEPH cannot be fully discounted, particularly since Factor VIII, a substrate of Factor V, has already been implicated in the disease 8.
Although the present study represents the largest of its kind in CTEPH, it does have several limitations. The use of “convenience” controls that were not necessarily drawn from the same population as the CTEPH cases may have confounded the study, although the results of the Kolmogorov–Smirnov test suggest that population stratification was not an issue. Furthermore, although there is evidence to support a pathophysiological role of both the fibrinogen Aα Thr312Ala and Factor V Leiden mutations in thrombotic disease, it is possible that these mutations themselves are not causal but, instead, in LD with other causal loci. Finally, the relatively small numbers of subjects involved meant that the study may not have been powered sufficiently to detect differences in SNP frequencies between the disease and control groups. The differences that were found could have been chance findings in the group of individuals studied, and thus further studies are needed in order to confirm these results. If these future studies were to include a second control group of subjects with a history of completely resolved acute PE, a better understanding of the mechanisms underlying CTEPH could potentially be gained.
In conclusion, the present study demonstrates an association between both the fibrinogen Aα substitution of threonine with alanine at position 312 and Factor V Leiden mutations and chronic thromboembolic pulmonary hypertension that adds to understanding of the disease. The presence of these polymorphisms may not only predispose patients to embolic phenomena but also confer resistance to fibrinolysis that subsequently sets the stage for thrombus organisation.
Statement of interest
None declared.
- Received May 7, 2007.
- Accepted November 21, 2007.
- © ERS Journals Ltd