Genome-wide association studies and musculoskeletal diseases

Bone and joint diseases are major causes of morbidity and mortality worldwide, and their prevalence is increasing as the average population age increases. Most common musculoskeletal diseases show significant heritability, and few have treatments that prevent disease or can induce true treatment-free, disease-free remission. Furthermore, despite valiant efforts of hypothesis-driven research, our understanding of the etiopathogenesis of these conditions is, with few exceptions, at best moderate. Therefore, there has been a long-standing interest in genetics research in musculoskeletal disease as a hypothesis-free method for investigating disease etiopathogenesis. Important contributions have been made through the identification of monogenic causes of disease, but the holy grail of human genetics research has been the identification of the genes responsible for common diseases. The development of genome-wide association (GWA) studies has revolutionized this field, and led to an explosion in the number of genes identified that are definitely involved in musculoskeletal disease pathogenesis. However, this approach will not identify all common disease genes, and although the current progress is exciting and proves the potential of this research discipline, other approaches will be required to identify many of the types of genetic variation likely to be involved.

Breakthroughs in genetic studies of ankylosing spondylitis

Ankylosing spondylitis (AS), the prototypic seronegative arthropathy, is known to be highly heritable, with >90% of the risk of developing the disease determined genetically. As with most common heritable diseases, progress in identifying the genes involved using family-based or candidate gene approaches has been slow. The recent development of the genome-wide association study approach has revolutionized genetic studies of such diseases. Early studies in ankylosing spondylitis have produced two major breakthroughs in the identification of genes contributing roughly one third of the population attributable risk of the disease, and pointing directly to a potential therapy. These exciting findings highlight the potential of future more comprehensive genetic studies of determinants of disease risk and clinical manifestations, and are the biggest advance in our understanding of the causation of the disease since the discovery of the association with HLA-B27.

Investigation of the role of ENPP1 and TNAP genes in chondrocalcinosis

Objectives. Extracellular inorganic pyrophosphate (ePPi) inhibits certain forms of pathological mineralization while promoting others. Three molecules involved in ePPi regulation are important candidates for the development of calcium pyrophosphate dihydrate chondrocalcinosis (CPPD CC). These include ANKH, ectonucleotide pyrophosphatase (ENPP1) and TNAP. We have previously showed that genetic variation in ANKH is a cause of autosomal dominant familial CC and also some sporadic cases of CPPD CC. We now investigate the possible role of ENPP1 and TNAP in CPPD CC. Methods. Exons, untranslated regions (UTR) and exon-intron boundaries of ENPP1 and TNAP were sequenced using ABI Big Dye chemistry on automated sequencers. Sixteen variants were identified (3 in ENPP1 and 13 in TNAP) and were subsequently genotyped in 128 sporadic Caucasian CPPD CC patients and 600 healthy controls using a combination of polymerase chain reaction/restriction fragment-length polymorphism analysis or using Taqman. Allele and genotype frequencies were compared between cases and controls using the χ 2 test. Linkage disequilibrium, haplotype and the single nucleotide polymorphism-specific analyses were also performed. This study had 80% power to detect an odds ratio of 2.2 or more at these loci. Results. No difference was observed in the allele or genotype frequencies between patients and controls at either ENPP1 or TNAP. Conclusions. Polymorphisms of ENPP1 and TNAP are not major determinants of susceptibility to CC in the population studied. Further studies of the aetiology of sporadic CPPD CC are required to determine its causes.

Combined analysis of three whole genome linkage scans for Ankylosing Spondylitis

Objective. Ankylosing spondylitis (AS) is a debilitating chronic inflammatory condition with a high degree of familiality (λs=82) and heritability (>90%) that primarily affects spinal and sacroiliac joints. Whole genome scans for linkage to AS phenotypes have been conducted, although results have been inconsistent between studies and all have had modest sample sizes. One potential solution to these issues is to combine data from multiple studies in a retrospective meta-analysis. Methods: The International Genetics of Ankylosing Spondylitis Consortium combined data from three whole genome linkage scans for AS (n=3744 subjects) to determine chromosomal markers that show evidence of linkage with disease. Linkage markers typed in different centres were integrated into a consensus map to facilitate effective data pooling. We performed a weighted meta-analysis to combine the linkage results, and compared them with the three individual scans and a combined pooled scan. Results: In addition to the expected region surrounding the HLA-B27 gene on chromosome 6, we determined that several marker regions showed significant evidence of linkage with disease status. Regions on chromosome 10q and 16q achieved 'suggestive' evidence of linkage, and regions on chromosomes 1q, 3q, 5q, 6q, 9q, 17q and 19q showed at least nominal linkage in two or more scans and in the weighted meta-analysis. Regions previously associated with AS on chromosome 2q (the IL-1 gene cluster) and 22q (CYP2D6) exhibited nominal linkage in the meta-analysis, providing further statistical support for their involvement in susceptibility to AS. Conclusion: These findings provide a useful guide for future studies aiming to identify the genes involved in this highly heritable condition. . Published by on behalf of the British Society for Rheumatology.

Effects of PTPN22 C1858T polymorphism on susceptibility and clinical characteristics of British Caucasian rheumatoid arthritis patients

Objectives. To confirm the association of a functional single-nucleotide polymorphism (SNP), C1858T (rs2476601), in the PTPN22 gene of British Caucasian rheumatoid arthritis (RA) patients and to evaluate its influence on the RA phenotype. Methods. A total of 686 RA patients and 566 healthy volunteers, all of British Caucasian origin, were genotyped for C1858T polymorphism by PCR-restriction fragment length polymorphism assay. Data were analysed using SPSS software and the χ 2 test as applicable. Results. The PTPN22 1858T risk allele was more prevalent in the RA patients (13.9%) compared with the healthy controls (10.3%) (P = 0.008, odds ratio 1.4, 95% confidence interval 1.09-1.79). The association of the T allele was restricted to those with rheumatoid factor (RF)-positive disease (n = 524, 76.4%) (P = 0.004, odds ratio 1.5, 95% confidence interval 1.1-1.9). We found no association between PTPN22 and the presence of the HLA-DRB1 shared epitope or clinical characteristics. Conclusions. We confirmed the previously reported association of PTPN22 with RF-positive RA, which was independent from the HLA-DRB1 genotype.

Genetic and genomic studies of PADI4 in rheumatoid arthritis

Objectives. Strong genetic association of rheumatoid arthritis (RA) with PADI4 (peptidyl arginine deiminase) has previously been described in Japanese, although this was not confirmed in a subsequent study in the UK. We therefore undertook a further study of genetic association between PADI4 and RA in UK Caucasians and also studied expression of PADI4 in the peripheral blood of patients with RA. Methods. Seven single-nucleotide polymorphisms (SNP) were genotyped using polymerase chain reaction (PCR)-restriction fragment length polymorphism in 111 RA cases and controls. A marker significantly associated with RA (PADI4_100, rs#2240339) in this first data set (P = 0.03) was then tested for association in a larger group of 439 RA patients and 428 controls. PADI4 transcription was also assessed by real-time quantitative PCR using RNA extracted from peripheral blood mononuclear cells from 13 RA patients and 11 healthy controls. Results. A single SNP was weakly associated with RA (P = 0.03) in the initial case-control study, a single SNP (PADI4_100) and a two marker haplotype of that SNP and the neighbouring SNP (PADI4_04) were significantly associated with RA (P = 0.02 and P = 0.03 respectively). PADI4_100 was not associated with RA in a second sample set. PADI4 expression was four times greater in cases than controls (P = 0.004), but expression levels did not correlate with the levels of markers of inflammation. Conclusion. PADI4 is significantly overexpressed in the blood of RA patients but genetic variation within PADI4 is not a major risk factor for RA in Caucasians.

The effect of transforming growth factor β1 gene polymorphisms in ankylosing spondylitis

Objectives. To determine whether genetic polymorphisms in or near the transforming growth factor β1 (TGFB1) locus were associated d with susceptibility to or severity of ankylosing spondylitis (AS). Methods. Five intragenic single-nucleotide polymorphisms (SNP) and three microsatellite markers flanking the TGFB1 locus were genotyped. Seven hundred and sixty-two individuals from 184 multiplex families were genotyped for the microsatellite markers and two of the promoter SNPs. One thousand and two individuals from 212 English and 170 Finnish families with AS were genotyped for all five intragenic SNPs. A structured questionnaire was used to assess the age of symptom onset, disease duration and disease severity scores, including the BASDAI (Bath Ankylosing Spondylitis Disease Activity Index) and BASFI (Bath Ankylosing Spondylitis Functional Index). Results. A weak association was noted between the rare TGFB1 + 1632 T allele and AS in the Finnish population (P = 0.04) and in the combined data set (P = 0.03). No association was noted between any other SNPs or SNP haplotype and AS, even among those families with positive non-parametric linkage scores. The TGFB1 +1632 polymorphism was also associated with a younger age of symptom onset (English population, allele 2 associated with age of onset greater by 4.2 yr, P = 0.05; combined data set, allele 2 associated with age of onset greater by 3.2 yr, P = 0.02). A haplotype of coding region SNPs (TGFB1 +869/ +915+1632 alleles 2/1/2) was associated with age of symptom onset in both the English parent-case trios and the combined data set (English data set, haplotype 2/1/2 associated with age of onset greater by 4.9 yr, P = 0.03; combined data set, haplotype 2/1/2 associated with greater age of onset by 4.2 yr, P = 0.006). Weak linkage with AS susceptibility was noted and the peak LOD score was 1.3 at distance 2 cM centromeric to the TGFB1 gene. No other linkage or association was found between quantitative traits and the markers. Conclusion. This study suggests that the polymorphisms within the TGFB1 gene play at most a small role in AS and that other genes encoded on chromosome 19 are involved in susceptibility to the disease.

Concordance of disease severity among family members with ankylosing spondylitis?

Objective. The heritability of disease activity and function in ankylosing spondylitis (AS) have been estimated at 0.51 and 0.63 (i.e., 51% and 63%), respectively. We examined the concordance of disease severity among family members in terms of disease activity, function, radiological change, prevalence of iritis, and juvenile onset. Methods. Disease activity and functional impairment due to AS were studied using the Bath AS Disease Activity Index (BASDAI) and Functional Index (BASFI) self-administered questionnaires; radiographic involvement was measured using the Bath AS Radiology Index (BASRI) scale. Familial correlation of BASDAI and BASFI was assessed in 406 families with 2 or more cases, using the program PAP. Parent-child and sibling-sibling concordance for iritis and juvenile AS were also studied in these families. Heritability of radiological disease severity based on the BASRI was assessed in 29 families containing 60 affected individuals using the program SOLAR. Results. Correlations between parent-child pairs for disease activity and function were 0.07 for both. Correlations between sibling pairs for disease activity and function were 0.27 and 0.36, respectively. The children of AS parents with iritis were more likely to develop iritis [27/71 (38%)] than children of non-iritis AS parents [13/70 (19%)] (p = 0.01). Parents with JAS were more likely to have children with JAS [17/30 (57%) compared to non-JAS parents 34/111 (30%)] (p = 0.002). The heritability of radiological disease severity based on the BASRI was 0.62. Conclusion. While correlation in severity between parent and child is poor, siblings do resemble each other in terms of severity, supporting the findings of segregation studies indicating significant genetic dominance in the heritable component of disease activity. Significant parent-child concordance for iritis and juvenile disease onset suggest that there are genetic risk factors for these traits independent of those determining the risk of AS itself. The finding of significant heritability of radiological change (BASRI) provides support using an objective measure for the observed heritability of the questionnaire-assessed disease severity scores, ASDAI and BASFI.

Osteoclast formation from circulating precursors in osteoporosis

Objective: An imbalance between bone formation and bone resorption is thought to underlie the pathogenesis of reduced bone mass in osteoporosis. Bone resorption is carried out by osteoclasts, which are formed from marrow-derived cells that circulate in the monocyte fraction. Ihe aim of this study was to determine the role of osteoclast formation in the pathogenesis of bone loss in osteoporosis. Methods: The proportion of circulating osteoclast precursors and their relative sensitivity to the osteoclastogenic effects of M-CSF, 1,25(OH)2D3 and RANKL were assessed in primary osteoporosis patients and normal controls. Results: Although there was no difference in the number of circulating osteoclast precursors in osteoporosis patients and normal controls, osteoclasts formed from osteoporosis patients exhibited substantially increased resorptive activity relative to normal controls. Although no increased sensitivity to the osteoclastogenic effects of 1,25(OH)2D3 or M-CSF was noted, increased bone resorption was found in osteoporosis peripheral blood mononuclear cell (PBMC) cultures to which these factors were added. Conclusion: Our findings suggest that osteoclast functional activity rather than formation is increased in primary involutional osteoporosis and that dexamethasone acts to increase osteoclast formation.

Non-inherited maternal HLA alleles are associated with rheumatoid arthritis

Background. Rheumatoid arthritis (RA) is strongly associated with a series of HLA-DRB1 alleles that encode a conserved sequence of amino acids (70Q/R K/R R A A74) in the DRβ1 chain, known as the shared epitope (SE). However 30% of patients are negative for DRB1*04 and 15% are SE-negative. Exposure to these alleles as non-inherited maternal antigens (NIMA) might explain this discrepancy. We undertook a family study to investigate the role of NIMA in RA. Methods. One hundred families, including the RA proband and both parents, were recruited. HLA-DRB1 genotyping was performed using an allele-specific polymerase chain reaction by standard methods. The frequencies of NIMA and non-inherited paternal antigens (NIPA) were compared using contingency tables and a two-tailed P test. We then reviewed four previously published studies of NIMA in RA and conducted an analysis of the combined data Results. We identified 36 families in which the proband was DRB1*04-negative and 13 in which the proband lacked the SE. There was an excess of DRB1*04 and SE NIMA (P=0.05) compared with NIPA. Combined analysis with previous studies showed that 53/231 mothers (23%) versus 25/205 fathers (12%) had a non-inherited DRB1*04 (P=0.003) and 30/99 mothers versus 18/101 fathers had a non-inherited SE allele (P=0.03). Conclusion. A role for HLA NIMA in RA is suggested by these results.

Genetic studies of disorders of calcium crystal deposition

In summary, although many factors are likely to be involved in regulating calcification and ossification processes, studies of the causation of articular chondrocalcinosis and disorders of spinal ossification, such as DISH and OPLL, implicate control over inorganic pyrophosphate levels as being one of the most important factors in their aetiopathogenesis. The findings of these studies may prove relevant to other rheumatic diseases in which ectopic ossification occurs, such as AS.

HLA-DR/DQ haplotype in rheumatoid arthritis: Novel allelic associations in UK Caucasians

Objective. To elucidate the relative importance of the HLA-DR and HLA-DQ loci in conferring genetic predisposition to rheumatoid arthritis (RA). Methods. HLA-DRB1 and HLA-DQB1 alleles were typed in a set of 685 patients with RA using sequence-specific polymerase chain reaction. Allele and phenotype frequencies were compared with those in 2 large sets of historical, ethnically matched healthy controls, using the relative predispositional effect method. Results. Positive association was confirmed with the shared epitope positive HLA-DRB1 alleles associated with RA in Caucasians. A significant susceptibility effect was observed with HLA-DRB1*09, described in other ethnically diverse populations but not in Caucasians. A significant underrepresentation of the HLA-DRB1*0103 variant was noted among the RA cases, supporting the proposed protective role of the DERAA motif at residues 70-74 of the DRβ molecule. No HLA-DRB1 independent association of the HLA-DQB1 alleles, implicated in predisposing to RA, was evident. Conclusion. These data corroborate the shared epitope hypothesis of susceptibility to RA and provide strong evidence for the DRB1 locus as the primary RA susceptibility factor in the HLA region.

Susceptibility to ankylosing spondylitis

Sir The association between HLA‐B27 (B27) and ankylosing spondylitis (AS) has been known for 25 yr. Familial aggregation in AS is well established, and first‐degree relatives of AS patients have been shown to be at increased risk of developing the disease. The recurrence risk in siblings of AS patients is quite uncertain, previous studies have variously reported recurrence risks between 6.9 and 27% [1, 2]. Accurate knowledge of the sibling recurrence risk is important both to advise families of the likelihood of disease recurrence, and in genetic statistical analyses utilizing Risch's recurrence risk ratio [3]. This study was designed to determine the risk of developing AS in siblings and to determine the role of the major histocompatibility complex in familial recurrence of AS....

Secretor status and ankylosing spondylitis

Objective. To assess the genetic association between secretor status and ankylosing spondylitis (AS). Methods. A restriction digest method for determining secretor status was developed; 166 patients with AS and 216 healthy British white controls were typed for secretor status using this method. Results. The frequency of nonsecretors among patients with AS (47/166, 28%) was not significantly different from controls (72/216, 33%). Conclusion. Secretor status does not influence susceptibility to AS.