University of Oxford, Institute of Molecular Medicine, Oxford,
1 Nuffield Orthopaedic Centre, Oxford, UK,
2 Departments of Biomathematics and Biostatistics, University of California, Los Angeles, USA and
3 Musgrave Park Hospital, Belfast, UK
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Methods. A total of 481 families have been ascertained in which at least two siblings have had joint replacement surgery of the hip, or knee, or hip and knee for idiopathic OA. Each candidate gene was targeted using one or more intragenic or closely linked microsatellite marker. The linkage data were analysed unstratified and following stratification by sex and by joint replaced (hip or knee).
Results. The analyses revealed suggestive linkage of the type IX collagen gene COL9A1 (6q12-q13) to a subset of 132 families that contained affected females who were concordant for hip OA (female-hip) with a P-value of 0.00053 and logarithm of the odds (LOD) score of 2.33 [corrected P-value of 0.0016, corrected LOD score of 1.85].
Conclusions.COL9A1 may therefore be a susceptibility locus for female hip OA. In addition, there was weak evidence of linkage to HLA/COL11A2 (6p21.3) in female hip OA with a corrected P-value of 0.016.
KEY WORDS: Osteoarthritis, Linkage, Affected sibling pairs, COL9A1.
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
It has often been noted in epidemiological studies that there is a female preponderance for OA. This may be accounted for by differential effects on the two sexes of environmental factors. However, a Finnish twin study has suggested that genetic susceptibility may be greater in women than men [4] and this result has been supported by a segregation analysis [5]. These studies have tended to concentrate on patients with hand, knee, or hand and knee OA. Some studies have attempted to quantify the risks of developing OA for specific joints and differences in the heritability values between joint groups have been reported [3, 6, 7]. These results suggest that there may be genetic heterogeneity of OA between the sexes and between different joint groups.
A number of candidate genes have been proposed as potential susceptibility loci for idiopathic OA. These genes predominantly encode structural proteins of the extracellular matrix. Particular emphasis has been placed on the gene for type II collagen (COL2A1) which is the major cartilage collagen. There are two principal reasons for this: (i) COL2A1 encodes the most abundant protein of the cartilage extracellular matrix and (ii) mutations in COL2A1 have been found in severa osteo- chondrodysplasias. Osteo-chondrodysplasias are rare monogenic diseases of cartilage and bone which often have OA as one of their principal phenotypic components. The involvement of COL2A1 in the osteo-chondrodysplasia ocular-Stickler syndrome is particularly interesting as this disease involves precocious, early-onset OA [8]. Less severe mutations of the COL2A1 gene could account for some cases of idiopathic OA. Indeed, linkage has been established to COL2A1 in a couple of rare families in which early-onset OA segregates as a Mendelian trait [9]. Furthermore, null alleles of COL2A1 have been reported in some individuals with idiopathic OA [10], whilst two studies have suggested association of polymorphisms in the vitamin D receptor gene (which is syntenic with COL2A1) with OA [11, 12]. However, other studies have excluded COL2A1 as a high impact susceptibility locus [13].
Other candidates proposed as potential susceptibility loci include genes that encode additional structural proteins of the cartilage extracellular matrix, particularly those that, like COL2A1, have been implicated in osteo-chondrodysplasias. These include the type XI collagen genes COL11A1 and COL11A2 (non-ocular Stickler disease) [14, 15], the type IX collagen gene COL9A2 (multiple epiphyseal dysplasia) [16] and the cartilage oligomeric matrix protein gene, COMP (multiple epiphyseal dysplasia and pseudoachondroplasia) [17].
Two loci which do not encode cartilage structural proteins but which have been considered as potential candidates are the HLA gene cluster on chromosome 6p and the 1-antitrypsin gene on chromosome 14q. These loci have been proposed on the basis that a subset of OA which is characterized by the presence of multiple Heberden's nodes on the distal interphalangeal joints (nodal generalized OA) may have an inflammatory component. Tentative association was established in a small group of unrelated affected individuals [18].
Since idiopathic OA rarely segregates as a Mendelian trait, and as both genetic and environmental factors play a role in the aetiology of the disease, its genetic dissection is not readily amenable to conventional family-based parametric linkage analysis. Instead, it is more suitable to the affected sibling pair linkage approach which does not specify the mode of inheritance.
There are a number of candidates which merit analysis as susceptibility loci for OA and in this paper we describe the results of an analysis of 11 such candidate genes. To test for linkage we used an affected sibling pair approach and recruited affected pairs using joint replacement surgery for idiopathic OA as our ascertainment criterion. Our aim was to use families whose idiopathic OA was severe and is therefore more likely to have involved a genetic component. We used a two-stage approach. In the first stage all candidate microsatellites were genotyped in 297 of our 481 families. Any candidate that demonstrated evidence for linkage at a P-value 0.05 in the 297 families or in a stratified subset of the families (stratified by sex and by hip or knee joint replacement) was then analysed in the remaining 184 families.
![]() |
Patients and methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Candidate genes
The candidate genes examined are shown in Table 1, which also lists the polymorphic markers used to target each candidate. Where available, intragenic markers were used. Other markers were chosen that were 5 cM or less from the candidate gene. For the COL2A1/VDR locus, the variable number of tandem repeat (VNTR) marker located immediately downstream of the COL2A1 gene [28] and two additional markers were used to cover chromosome 12q13-q14 (D12S368 and D12S87).
|
Polymerase chain reaction (PCR), gel electrophoresis and genotyping
Either the forward or the reverse primer in a PCR pair was fluorescently labelled and amplification products were electrophoresed through Applied Biosystems 373 or 377 Automated DNA Sequencers as described by Reed et al. [29]. Alleles were sized using Applied Biosystems Genescan version 2.0.2 and Genotyper version 1.1 software.
Linkage analysis strategy
The linkage analysis strategy was to genotype all markers in a first stage containing 297 of our 481 families. Any marker that had a P-value 0.05 in stage 1 would then be examined in the remaining 184 families (stage 2). The aim of this strategy was to take through to stage 2 only those markers that demonstrated reasonable evidence of linkage. In stage 2 we were not necessarily expecting to repeat any positive linkage results of stage 1 but instead were looking for further evidence of linkage, even if only moderate: if a marker's P-value for stages 1 and 2 combined was less than the P-value for stage 1 only then this would support linkage at that marker. There was no difference in the ascertainment criteria between the stage 1 and stage 2 familiesonce a reasonably large number of families had been collected we began stage 1 of our linkage strategy.
As well as analysing all the stage 1 families together, these families were also stratified into those containing affected females-only (132 families), affected males-only (60 families), hips-only (male and/or female) (194 families), knees-only (male and/or female) (34 families), affected females-only who had undergone hip replacement but not knee replacement (female-hip) (85 families) and affected males-only who had undergone hip replacement but not knee replacement (male-hip) (44 families). We did not stratify by female-knee or male-knee as the number of families was too low (16 and four, respectively) to allow reliable inference of linkage. Again, any marker that had a P-value 0.05 following stratification in stage 1 would then be examined in the stage 2 families.
Stratification criteria
We stratified by sex, by joint replaced (hip or knee) and by sex combined with joint replaced. The stratification criteria have been published elsewhere [19].
Statistical analysis
Linkage analysis was performed using the SIBPAIR module of the ANALYZE package (ftp://linkage.cpmc.columbia.edu). This module is able to use data from siblings to determine identity by descent (IBD) allele transmittance. In the linkage analysis siblings who had not undergone joint replacement were given a clinical status of unknown. The SIBPAIR module produces a single-point logarithm of the odds (LOD) score and its asymptotic P-value. Allele frequencies were calculated from the input data using GAS (http: //users.ox.ac.uk/ayoung/gas.html). As well as the unstratified data we tested six strata in stage 1, three strata in stage 2 and three strata when stages 1 and 2 were combined. We therefore adjusted the stratified LOD scores and P-values to correct for the number of strata tested [30]. For stage 1 we deducted log 6 = 0.78 from original LOD scores and multiplied original P-values by 6, whilst for stage 2, and for stage 1 and 2 combined, we deducted log 3 = 0.48 from original LOD scores and multiplied original P-values by 3.
![]() |
Results |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
Stratified linkage analysis.
Table 3 lists the uncorrected P-values and LOD scores for the stratified data obtained for each marker in the 297 families of stage 1. Twenty of the 114 P-values were
0.05. Following correction for multiple testing [30], five of these 20 P-values were still
0.05, representing three of the 11 candidates: COL9A1 marker 5098B2 for hips-only (P = 0.042) and for female-hips (P = 0.028), HLA/COL11A2 marker D6S265 for hips-only (P = 0.024) and HLA/COL11A2 marker D6S273 for females-only (P = 0.018) and for female-hips (P = 0.014) (Table 4
).
|
|
Stage 2
Unstratified linkage analysis.
None of the six markers had a P-value 0.05 in the unstratified families of stage 2 (Table 5
). After combining data for screens 1 and 2, only DTDST marker (CA)n had a P-value
0.05 (P = 0.022). This compares with a P-value of 0.0093 for stage 1 onlyincreasing the number of families studied had reduced the evidence for linkage at this marker.
|
Stratified linkage analysis.
In the corrected results from the stage 1 stratification COL9A1 and HLA/COL11A2 markers were significant at P 0.05 in one or more of the following strata: females-only, hips-only or female-hips. We therefore restricted our stratification analysis of COL9A1 and HLA/COL11A2 in stage 2, and stages 1 and 2 combined, to these three strata.
In stage 2, the COL9A1 marker 5098B2 had uncorrected P-values 0.05 for females-only (P = 0.034) and for female-hips (P = 0.024) (Table 6
). However, following correction for the three tests performed, none of the four markers targeting COL9A1 and HLA/COL11A2 had a P-value
0.05 in any of the three strata, although 5098B2 did approach significance in females-only (P = 0.10) and female-hips (P = 0.072). When stages 1 and 2 were combined (Table 7
) 5098B2 had a corrected P-value
0.05 in females-only (P = 0.01) and female-hips (P = 0.0016) and approached significance in hips-only (P = 0.051). 50912B1 was also significant for female-hips (P = 0.042). The HLA/COL11A2 marker D6S273 had a corrected P-value
0.05 in females-only (P = 0.01) and female-hips (P = 0.016). D6S265 was also significant in female-hips (P = 0.025) and approached significance in females-only (P = 0.054).
|
|
|
Nodal OA
As mentioned in the Introduction, the association that has previously been reported for both the HLA gene cluster and the 1-antitrypsin gene is to nodal generalized OA [18]. We therefore analysed the D6S265, D6S273 and the intragenic
1-antitrypsin markers in those families from stage 1 which contained at least two affected siblings who were concordant for the presence of three or more Heberden's nodes. Of the 297 families in stage 1, 29 were concordant for this phenotype and the P-values obtained for each marker were >0.1 (uncorrected).
![]() |
Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Of the 11 genes targeted, two demonstrated evidence of suggestive linkage at a P-value of 0.05 when the 297 families of stage 1 were analysed without stratification: the HLA/COL11A2 locus and DTDST. Analysis of these two genes in an additional 184 families demonstrated no evidence of linkage and when the data for stages 1 and 2 were combined only the DTDST gene still supported linkage at P
0.05. However, the P-value for the DTDST (CA)n marker was greater in the combined stages than in stage 1 only, suggesting that the positive result in stage 1 was probably a false positive.
When we stratified our stage 1 results and corrected for multiple testing, the COL9A1 marker 5098B2 was significant in hips-only and female-hips, the HLA/COL11A2 marker D6S265 was significant in hips-only and female-hips and the HLA/COL11A2 marker D6S273 was significant in females-only. No markers targeting the remaining eight candidate genes were significant in the six strata tested. Despite correction, these positive stratification results could still be false. The two-stage nature of our linkage analysis strategy enabled us to test this possibility in that we were able to genotype the COL9A1 and HLA/COL11A2 markers in additional families. After correction, none of the four markers genotyped in stage 2 were significant in any of the three strata tested. However, this was not unexpected: to have repeated the significant results in an independent data set (which had lower numbers in each of the strata tested) would have been surprising [31]. Nevertheless, 5098B2 did approach significance in females-only and female-hips. A more informative analysis was the comparison between the P-values obtained in stage 1 with those of stages 1 and 2 combined: a reduction in P-values when more families were genotyped (i.e. when the power of the study was increased) would support the positive results of stage 1. These comparisons revealed a reduced P-value in the combined vs stage 1 for 5098B2 in females-only and female-hips (with a very moderate increase in P-value for hips-only), a reduced P-value for the second COL9A1 marker 50912B1 in female-hips, a reduced P-value for the two HLA/COL11A2 markers (D6S265 and D6S273) in females-only with a reduced P-value for D6S265 in female-hips and only a very moderate increase in the P-value for D6S273 in this strata. Overall, therefore, increasing the number of families increased the evidence for linkage of the COL9A1 and the HLA/COL11A2 markers in females-only and female-hips, suggesting that these positive results are probably real.
Why the linkage should be restricted to females-only deserves comment, particularly so for COL9A1 where the difference between females-only and males-only is particularly striking. The obvious possibility is that the susceptibility at COL9A1 is mediated through a pathway that involves a female-specific hormonal element. Females are at a greater risk than males of developing OA, and this risk increases after the menopause [1]. Indeed, hormone replacement therapy can affect the development and/or progression of OA in females [32]. Detailed analysis of the regulation of the COL9A1 gene, with particular emphasis on controlling elements that have a hormonal component, may be helpful in dissecting the apparent female bias in susceptibility at this locus.
Type IX collagen is a heterotrimeric protein composed of 1(IX),
2(IX) and
3(IX) polypeptide chains that are encoded by the COL9A1, COL9A2 and COL9A3 genes, respectively. Type IX collagen is a quantitatively minor cartilage collagen that decorates the type II collagen fibrils and this interaction serves to regulate the growth of the type II fibril. Type IX is composed of three collagenous (COL13) and four non-collagenous domains (NC14) with the NC4 domain encoded by COL9A1 being significantly larger than the NC4 domains encoded by COL9A2 or COL9A3. This large COL9A1 NC4 domain may interact with cartilage proteoglycans and could help stabilize the cartilage extracellular matrix. Two mouse models have shed light on to the function of type IX collagen, and more specifically the COL9A1-encoded
1(IX) chain. In the first model, a truncated form of col9a1 resulted in a mild osteo-chondrodysplasia with OA [33] whilst in the second model, a col9a1 knock-out mouse had no congenital abnormality but developed severe OA that was comparable in timing and pathology to human idiopathic OA [34]. These models highlighted the importance of the
1(IX) chain in stabilizing cartilage and demonstrated that mutations in Col9a1 can result in an OA phenotype.
Overall, we have identified two genetic loci that demonstrate evidence of suggestive linkage to severe hip OA in females. The COL9A1 gene sequence is available in public databases and this can be used for the identification of coding and regulatory single nucleotide DNA variants that would be a starting point for detailed association analysis of this gene with OA.
Stratification increases the level of genetic homogeneity and can therefore assist in the mapping of loci for complex traits. Our analysis highlights the potential utility of this approach for OA.
Our analysis cannot exclude as susceptibility loci all the genes that were negative in this study. It may be that one or more of these genes is of low impact in the severe end-stage form of OA examined and that we lacked the power to detect this effect. Alternatively, they may encode susceptibility for different forms of OA (i.e. hand or disc) that we did not examine.
![]() |
Acknowledgments |
---|
![]() |
Notes |
---|
![]() |
References |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|