| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Thrombosis |
From the Twin Research & Genetic Epidemiology Unit (M.d.L., T.A., H.S., T.D.S.), St Thomas Hospital, London, UK; Department of Biological Psychology (M.d.L.), Vrije Universiteit, Amsterdam, The Netherlands; Georgia Prevention Institute (H.S., D.G.), Department of Pediatrics, Medical College Georgia, Augusta; and Academic Unit of Molecular Vascular Medicine (T.S.F., K.S., P.J.G., R.A.S.A.), The Leeds Institute for Genetics, Health and Therapeutics, Faculty of Medicine and Health, University of Leeds, UK.
Correspondence to Robert Ariëns, PhD, Academic Unit of Molecular Vascular Medicine, LIGHT Laboratories, Clarendon Way, University of Leeds, Leeds LS2 9JT, UK. E-mail r.a.s.ariens{at}leeds.ac.uk
| Abstract |
|---|
|
|
|---|
Methods and Results A total of 201 dizygotic twin pairs were genotyped for 1 promoter and all common nonsynonymous coding polymorphisms in the FXIII A-subunit gene: 246G>A, V34L, Y204F, P564L, V650I, and E651Q. Tests of linkage, association, and combined linkage and association were performed using QTDT software. Significant linkage to the V34L polymorphism (P=5x1012) as well as association (P=3x1049) was observed. Adjusting for association while performing linkage made the linkage signal disappear for the V34L polymorphism (from
2=47.55, P=5x1012 to
2=1.30, P=0.25). Only haplotypes containing the 34L allele showed association with FXIII activation.
Conclusion Testing multiple SNPs in the FXIII A-subunit gene indicates that V34L is the main functional polymorphism influencing FXIII activation.
A joint test of linkage and association of 6 single-nucleotide polymorphisms in the factor XIII-A subunit gene, genotyped in 201 dizygotic female white twin pairs, indicates that a common variant (V34L) in the factor XIII-A subunit gene is the main functional variant influencing factor XIII activation.
Key Words: joint linkage and association analysis functional locus FXIIIV34L FXIII activation twins
| Introduction |
|---|
|
|
|---|
Factor XIII (FXIII), when activated by thrombin, covalently crosslinks fibrin multimer strands, thereby enhancing the mechanical strength of a clot and increasing the resistance to fibrinolysis.5 The prevalence of a common variant in the FXIII-A gene (V34L) was lower in patients with myocardial infarction (MI) than those without, suggesting a protective role for this polymorphism.6,7 Paradoxically, a stepwise increase in FXIII activation with possession of the protective Leu allele has been found, with no differences in FXIII-A and B-subunit levels.8,9 In vitro studies have shown that 34L is cleaved faster by thrombin than the Val variant,10 with an increase in proteolytic efficiency by 2- to 3-fold.10,11 How this relates to the risk of vascular disease remains to be fully elucidated.
We have previously shown moderate to high degree of genetic influence on the variation for FXIII levels, explaining 82% of FXIII activation, 64% and 41% for FXIII-A and B-subunit levels, respectively.12 It is currently unknown to what degree single-nucleotide polymorphisms (SNPs) in the FXIII genes contribute to the genetic variability. This study aimed to investigate which SNP is the main functional polymorphism influencing FXIII activation among 1 common promoter and all common (allele frequency >1%) nonsynonymous polymorphisms in the FXIII A-subunit gene. To this end, we performed combined linkage and association analysis in a sample of 201 dizygotic (DZ) twin pairs.
| Subjects and Methods |
|---|
|
|
|---|
Blood Sampling and Factor XIII Assays
A venous blood sample was taken between 8:00 and 10:00 AM after an overnight fast. The blood was taken within 5 minutes from the cotwin into 0.13 mol/L trisodium citrate vacutainers (Becton Dickinson, France) and kept on ice for tests of fibrinolysis and at room temperature for tests of coagulation as described previously.12 Within 1 hour of collection, the samples were centrifuged at 2560g for 20 minutes to obtain platelet-poor plasma, frozen in aliquots in liquid nitrogen (LN2) and stored at 40°C until analysis. FXIII activation rate (previously called FXIII activity) was analyzed by a kinetic biotinamido pentylamineincorporation assay and FXIII A-subunit and B-subunit by sandwich ELISAs as described.15
DNA Extraction
DNA was extracted at the Twin Research Unit using the BACC2 DNA extraction kit (Nucleon Biosciences, Glasgow, UK) from 10 mL of venous blood collected on 1.6 mg/mL EDTA. Extracted DNA was diluted in Tris-EDTA buffer and stored at 20°C until further analysis. Amplification of DNA was carried out by polymerase chain reaction (PCR) in a final volume of 25 µL, containing 100 ng of genomic DNA. PCR amplification was carried out using a PTC-100 thermal cycler (MJ Research, Watertown, Mass).
Genotyping of Six SNPs in the FXIII A-subunit Gene
The twins were genotyped for a common promoter and all common nonsynonymous coding polymorphisms in the FXIII A-subunit gene (references for genome and messenger sequences respectively: NT 034880 and NM 000129) located on chromosome 6. The following 6 SNPs were genotyped: g.246G>A (rs1024231) in the promoter region (numbered in relation to the transcription start site; 2231 from translation start codon), V34L (c.103G>T; rs5985), Y204F (c.614A>T; rs3024477), P564L (c.1694C>T; rs5982), V650I (c.1951G>A; rs5987), and E651Q (c.1954G>C; rs5988). For the coding SNPs, the amino acid substitutions are numbered for the mature protein and the base substitutions are numbered from the translation start site. The 246G>A promoter polymorphism was determined by PCR amplification and FokI restriction analysis. Restriction products were separated on a 4% Metaphore gel (Flowgen, Nottingham, UK). The G allele resulted in a single band of 97 bp, and the A allele resulted in bands of 85 and 12 bp. V34L is a common G-to-T transition in exon 2 coding for a substitution of valine with leucine at residue 34. The V34L polymorphism genotyping by PCR amplification and DdeI digestion was as described by Kangsadalampai and Board.16 The DdeI digestion products were resolved on a 3% agarose gel. The presence of the G allele produces a 192-bp fragment and the presence of the T allele fragments of 161 and 31 bp. Y204F was analyzed by PCR amplification followed by restriction digestion with RsaI. The phenylalanine allele PCR product is cut into 92- and 20-bp fragments, the tyrosine allele product remains uncut at 112 bp. P564L constitutes a C-to-T transition at +1694 in exon 12 and alters proline 564 to leucine in barrel 1. FXIIIE651Q is a G-to-C transition at +1954 in exon 14 and alters glutamic acid 651 to glutamine. P564L and E651Q were genotyped as described previously.17 V650I was PCR amplified and analyzed by restriction digestion with HinfI. The G allele leads to 98- and 26-bp bands, and the C allele remains uncut at 124 bp. All primer sequences and PCR conditions are available on request.
Preliminary Analysis
FXIII activation was log transformed before linkage and association analysis to obtain a distribution that was approximately normal. This log-transformed variable was adjusted for age, date of blood sampling (ie, seasonal variation), and fasting status as described previously.12 Tests of HardyWeinberg equilibrium were carried out in STATA 8 (StatCorp, College Station, Tex). Pairwise linkage disequilibrium (LD) coefficients were reported as the ratio of the nonstandardized coefficient to the maximal value (D'=D/Dmax).
Linkage and Association
Tests for linkage, allelic association, and joint linkage and association were carried out using QTDT (Quantitative Transmission Disequilibrium Tests) software.18 The simplest test for linkage analysis using sib pairs is the method by Haseman and Elston,19 in which it is postulated that siblings (or twins) who share more alleles identical by descent (IBD) will resemble each other more if that locus (or a locus in close proximity) influences the trait. Using maximum likelihood variance components analysis, the linkage component can be modeled as part of the covariance structure. For this, the estimation of the probability distribution of sharing 0, 1, or 2 alleles IBD for each SNP was accomplished using Genehunter 2 software20 and used in QTDT software to test for linkage, ie, whether sharing alleles IBD affects phenotypic similarity (or covariance) between the DZ pairs.
Tests of overall allelic association, such as implemented in QTDT, compare mean levels of the hemostatic protein between genotype groups while taking into account the nonindependence of the twin pairs.
Fulker et al21 developed a method for simultaneous modeling of association and linkage for quantitative traits using sib pair data that also control for population stratification. Combined linkage and association analysis is a powerful tool for pinpointing functional loci responsible for a linkage signal. Assuming there is significant linkage before modeling association, the extent to which evidence for linkage diminishes in the joint test of linkage and association reflects the proximity of the marker to the functional quantitative trait locus (QTL).21,22 If the linkage signal of a certain locus is entirely explained by modeling association (ie, the linkage signal is no longer significant), that particular locus is either the functional locus or in tight LD with the functional locus.
Apart from close matching of environmental factors, an additional advantage of using family data such as DZ twins is that it is possible to test and control for population stratification as a source of spurious associations by conducting within family association tests such as the quantitative transmission disequilibrium tests implemented in QTDT.22,23
A sample size of 201 DZ twin pairs gave us 80% power (
=0.05) to detect a QTL explaining 19% of the variance with linkage.24 For the within-family association and overall association tests, we had 80% power (
=0.05) to detect loci explaining 5% and 2% of the variance, respectively.24
Details of our approach to test the association of statistically inferred haplotypes with continuous traits have been described previously.25 In short, we adapted haplotype trend regression from Zaykin et al for the analysis of related subjects such as twins26 by replacing the linear regression with the GEE procedure using STATA. The probabilities of haplotype pairs were estimated by PHASE 2.0 software.27
| Results |
|---|
|
|
|---|
|
|
Table 3 presents the differences in levels of FXIII activation for the different SNPs. Largest effect is seen for V34L, where twins homozygote for the Leu allele had 72% higher levels of untransformed FXIII activation compared with the Val homozygotes. In regression analysis, the polymorphism explained 39% of the variation in FXIII activation. Smaller, and opposite, effects were seen for the 246G>A and the P564L genotypes, where the homozygotes of the more frequent allele had 38% and 14.5% higher levels compared with the homozygotes of the less frequent allele, respectively. Only 1 twin had 2 alleles of the less frequent Phe allele of Y204F, but twins homozygote for the more frequent Tyr allele had 21.6% higher levels than the heterozygote twins (Table 3). Significant association to levels of FXIII activation was observed for 246G>A, V34L, Y204F, and P564L but not for V650I or E651Q. Permutation tests in QTDT for association of FXIII activation levels using only the within-family information (to account for possible undetected population stratification) calculated the empirical probability values at 0.0003 for 246G>A, 5x105 for V34L, 0.07 for Y204F, and 0.07 for P564L, with, respectively, 65, 75, 13, and 59 twin pairs informative for analysis.
|
Table 4 shows haplotype frequencies and results of haplotype association tests. Haplotype analyses confirmed individual SNP associations. Only haplotypes 2 and 7, which contained the 34Leu allele, showed significantly higher levels of FXIII activation (P
9x1010) compared with the most common haplotype 1: 1-1-1-1-1-1. The ß coefficient for haplotype 2 on FXIII activation was 66.5, meaning that subjects homozygote for haplotype 2 had FXIII activation levels that were 66.5 higher compared with subjects homozygote for the common haplotype 1. The overall probability value of the haplotype test was highly significant (P=1x1041) with the haplotypes, explaining 31.7% of the variance in FXIII activation levels.
|
Significant linkage of FXIII activation was shown to V34L (P=5x1012) and marginally to E651Q (P=0.03) but not to any of the other SNPs (Table 5). Linkage in the presence of association only showed a drop in fit for the V34L polymorphism (from
2=47.55, P=5x1012, to
2=1.30, P=0.25), indicating that this locus is either the functional polymorphism or in strong LD with the functional locus (Figure and Table 5). No drop in fit was observed for E651Q for which initial linkage was detected (P=0.03), indicating that this SNP is a marker linked to the functional polymorphism.
|
|
| Discussion |
|---|
|
|
|---|
Although there was clear evidence for linkage to V34L, no significant linkage could be seen for either of the flanking markers of V34L. For 246G>A, this was especially unexpected given its close proximity to V34L and its equal informative marker heterozygosity (Table 3). One flanking marker further from V34L, E651Q, showed marginal evidence for linkage. When calculating the empirical probability values for the associations robust to population stratification (ie, using only the within family component), evidence for association of V34L and 246G>A were still significant at P=5x105 and P=0.0003, respectively, but evidence for association to both Y204F and P564L reduced to P=0.07. This is probably attributable to the lower LD observed between V34L and P564L compared with V34L and 246G>A, and the low informative sample size (n=13) for the within family tests for Y204F.
The National Center for Biotechnology Information (NCBI) database (http://www.ncbi.nlm.nih.gov/SNP) lists more SNPs (10 in total) in the FXIII-A gene (F13A1) than those included in our study. The reasons for not including some of these SNPs in our study were that either the allele frequency was smaller than 0.01 (ie, heterozygosity as reported on NCBI is smaller than [0.99x0.01]x2=0.0198) or that there were no population data available for the calculation of allele frequency. It is likely that the SNPs listed on NCBI without population data have an allele frequency of <0.01 in our population from Northwest Europe, as we have extensively analyzed the FXIII-A gene by heteroduplex analysis in this population and have found no other common nonsynonymous SNP than those investigated in our current study (T.S. Futers and P.J. Grant, unpublished data). We made 2 exceptions to these SNP selection criteria: (1) we did not choose to analyze rs3024472, because its NCBI heterozygosity estimate of 0.023 was based on an unreliably small population sample size, and we have not observed this SNP in our population by heteroduplex analysis in the past; and (2) we included promoter polymorphism rs1024231, which is the only promoter variant for which previous studies had indicated a potential effect on FXIII gene expression.30 As changes in gene expression could have influenced our analysis of FXIII activation, this was an important variant to control for. However, our study and methodology was not designed to specifically investigate levels of gene expression, and hence no further noncoding variants were included.
Although we genotyped a promoter polymorphism and all common nonsynonymous coding polymorphisms with an allele frequency >0.01 in the FXIII A-subunit gene in our white population, this does not rule out that other as yet unknown regulatory variants are present elsewhere in the genome. For example, the V34L polymorphism explains 39% of the variation in FXIII activation. If we compare this to the estimated heritability of 82% (the phenotypic variation that can be explained by genetic factors) for FXIII activation, it represents a relatively large contribution of a single SNP to the overall genetic influence. One consequence of this large contribution is that it enabled us to detect linkage to this locus using a relatively modest sample size of 201 DZ twin pairs. However,
50% of the genetic variation remains unexplained and may be accounted for by other (regulatory) genes.
The V34L polymorphism has been associated with decreased cardiovascular risk in a case control study of patients undergoing coronary angiography, stratified by history of MI.6 This finding was confirmed in some, but not all, studies,31 and protection against thrombosis was additionally described in relation to venous thrombosis and ischemic cerebrovascular disease.31 In another study, only the V34L polymorphism (not P564L or E651Q) showed a difference in genotype prevalence in patients with coronary artery disease who had an MI compared with either patients without MI and controls.17 The results from that study supported in vitro studies that indicated that V34L had a direct role in influencing FXIII function, although the conundrum remained that the protective allele (34Leu) was associated with marked increases in FXIII activation,8 a finding apparently at odds with the molecular epidemiology. Subsequently, other studies confirmed that 34Leu enhances FXIII activation.10,11,29 It has been proposed that as FXIII 34Leu has more extensive contact with the binding region for the FXIII activation peptide on thrombin compared with 34Val, which may lead to 34Leu being more susceptible to wasteful conversion to the activated enzyme to provide a protective effect against cardiovascular disease,11 but other mechanisms likely play a role.
Initial studies of fibrin structure and function showed that individuals possessing the protective Leu allele produced fibrin with thinner fibrin fibers and smaller pores10; however, this finding is more in support of a hypothesis of increased CAD risk as observed earlier for these kind of fibrin structures by Fatah et al.32 It was also observed that individuals possessing the protective Leu allele, but who did experience MI, had higher levels of plasminogen activator inhibitor (PAI)-1 and were more insulin resistant, suggesting that additional risk factors for CAD may negate the protective effect of the Leu allele.33 In another study, we observed that effects of 34Leu on fibrin structure were reversed in the presence of high concentrations of fibrinogen, indicating that protection against thrombosis only occurs under conditions of raised fibrinogen.34 These geneenvironment interactions will make univariate (linkage and association) studies difficult to interpret and may represent a clear example of the general reason why there has been relatively little success in finding hemostatic genetic polymorphisms that link to CAD outcome, even if they explain a relatively large proportion of the intermediate (hemostatic) phenotypes (see Lane and Grant4). It is likely that both environmental and genetic risk factors act in concert, and taking these interactions into account may improve gene discovery35,36 and their applied value in future prediction and treatment of common disease. For example, specific genes may have relatively large effects on intermediate phenotypes but, because of interaction with a number of other factors, may have low predictive value on disease outcome. Investigating these genes (or intermediate phenotypes) in isolation may give poor results at both statistical and practical level. Nonetheless, studies into the associations between polymorphisms and phenotype may provide crucial insights into the underlying mechanisms and enhance understanding of disease processes.
Our study indicates that of 6 SNPs in the FXIII A-subunit, V34L is the main functional polymorphism that consistently associates with FXIII activation. The clear effect of this polymorphism on activation and the reports of its cardioprotective effect are an incentive to study the function of FXIIIV34L more closely. These studies could provide new opportunities for the prediction and treatment of atherothrombotic coronary artery disease.
| Acknowledgments |
|---|
Sources of Funding
This study was supported by the British Heart Foundation (FS/99010 and FS/99050) and the Medical Research Council (G/9900904). TwinsUK Registry is supported by the Wellcome Trust, the Arthritis and Rheumatism Campaign, the Framework 5 GenomEUtwin program, and the Framework 6 EuroClot project.
Disclosures
None.
| Footnotes |
|---|
| References |
|---|
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
ATVB Home | Subscriptions | Archives | Feedback | Authors | Help | AHA Journals Home | Search Copyright © 2006 American Heart Association, Inc. All rights reserved. Unauthorized use prohibited. |