HLA class I associations of ankylosing spondylitis in the white population in the United Kingdom

Objective - To investigate the HLA class I associations of ankylosing spondylitis (AS) in the white population, with particular reference to HLA-B27 subtypes. Methods - HLA-B27 and -B60 typing was performed in 284 white patients with AS. Allele frequencies of HLA-B27 and HLA-B60 from 5926 white bone marrow donors were used for comparison. HLA-B27 subtyping was performed by single strand conformation polymorphism (SSCP) in all HLA-B27 positive AS patients, and 154 HLA-B27 positive ethnically matched blood donors. Results - The strong association of HLA-B27 and AS was confirmed (odds ratio (OR) 171, 95% confidence interval (CI) 135 to 218; p < 10-99). The association of HLA-B60 with AS was confirmed in HLA-B27 positive cases (OR 3.6, 95% CI 2.1 to 6.3; p < 5 x 10-5), and a similar association was demonstrated in HLA-B27 negative AS (OR 3.5, 95% CI 1.1 to 11.4; p < 0.05). No significant difference was observed in the frequencies of HLA-B27 allelic subtypes in patients and controls (HLA-B*2702, three of 172 patients v five of 154 controls; HLA-B*2705, 169 of 172 patients v 147 of 154 controls; HkA-B*2708, none of 172 patients v two of 154 controls), and no novel HLA-B27 alleles were detected. Conclusion - HLA-B27 and -B60 are associated with susceptibility to AS, but differences in BLA-B27 subtype do not affect susceptibility to AS in this white population.

Gene expression profiling reveals a downregulation in immune-associated genes in patients with AS

Objective: To identify differentially expressed genes in peripheral blood mononuclear cells (PBMCs) from patients with ankylosing spondylitis (AS) compared with healthy individuals. Methods: RNA was extracted from PBMCs collected from 18 patients with active disease and 18 gender-matched and age-matched controls. Expression profiles of these cells were determined using microarray. Candidate genes with differential expressions were confirmed in the same samples using quantitative reverse transcription-PCR (qRT-PCR). These genes were then validated in a different sample cohort of 35 patients with AS and 18 controls by qRT-PCR. Results: Microarray analysis identified 452 genes detected with 485 probes which were differentially expressed between patients with AS and controls. Underexpression of NR4A2, tumour necrosis factor AIP3 (TNFAIP3) and CD69 was confirmed. These genes were further validated in a different sample group in which the patients with AS had a wider range of disease activity. Predictive algorithms were also developed from the expression data using receiver-operating characteristic curves, which demonstrated that the three candidate genes have ∼80% power to predict AS according to their expression levels. Conclusions: The findings show differences in global gene expression patterns between patients with AS and controls, suggesting an immunosuppressive phenotype in the patients. Furthermore, downregulated expression of three immune-related genes was confirmed. These candidate genes were also shown to be strong predictive markers for AS.

Analysis of killer immunoglobulin-like receptor genes in ankylosing spondylitis

Objectives: To assess the possible association of killer immunoglobulin-like receptor (KIR) genes, specifically KIR3DL1, KIR3DS1 and KIR3DL2, with ankylosing spondylitis (AS). Methods: 14 KIR genes were genotyped in 200 UK patients with AS and 405 healthy controls using multiplex polymerase chain reaction. Sequence-specific oligonucleotide probes were used to subtype 368 cases with AS and 366 controls for 12 KIR3DL2 alleles. Differences in KIR genotypes and KIR3DL2 allele frequencies were assessed using the χp2p test. Results: KIR3DL1 and KIR3DS1 gene frequencies were very similar in cases with AS and controls (odds ratio = 1.5, 95% confidence interval 0.8 to 3.0, and odds ratio = 1.02, 95% confidence interval 0.2 to 5.3, respectively). KIR3DL2 allele frequencies were not significantly different between cases with AS and controls. Conclusions: Neither the KIR gene content of particular KIR haplotypes nor KIR3DL2 polymorphisms contribute to AS.

Genetic testing for haemochromatosis in patients with chondrocalcinosis

Hereditary haemochromatosis (HH) is the most common lethal monogenic human disease, affecting roughly 1 in 300 white northern Europeans. Homozygosity for the C282Y polymorphism within the HFE gene causes more than 80% of cases, with compound heterozygosity of the C282Y and H63D polymorphism also increasing susceptibility to disease. The aim of this study was to determine the frequency of the C282Y and H63D polymorphisms in the disease, and to assess the risk of HH in heterozygotes for the C282Y polymorphism. 128 patients were recruited because of either radiographic chondrocalcinosis (at least bicompartmental knee disease or joints other than the knee involved) or CPPD pseudogout. Genotyping of the HFE C282Y and H63D mutations was performed using PCR/SSP and genotypes for the C282Y polymorphism confirmed by PCR/RFLP. Historical white European control data were used for comparison. Two previously undiagnosed C282Y homozygotes (1.6%), and 16 C282Y heterozygotes (12.5%), including four (3.1%) C282Y/ H63D compound heterozygotes were identified. This represents a significant overrepresentation of C282Y homozygotes (relative risk 3.4, p-0.037), but the number of heterozygotes was not significantly increased. At a cost per test of £1 for each subject, screening all patients with chondrocalcinosis using the above ascertainment criteria costs only £64 for each case of haemochromatosis identified, clearly a highly cost effective test given the early mortality associated with untreated haemochromatosis. Routine screening for haemochromatosis in patients with appreciable chondrocatcinosis is recommended.

Susceptibility to ankylosing spondylitis in twins: The role of genes, HLA, and the environment

Objective To determine the relative effects of genetic and environmental factors in susceptibility to ankylosing spondylitis (AS). Methods Twins with AS were identified from the Royal National Hospital for Rheumatic Diseases database. Clinical and radiographic examinations were performed to establish diagnoses, and disease severity was assessed using a combination of validated scoring systems. HLA typing for HLA-B27, HLA-B60, and HLA-DR1 was performed by polymerase chain reaction with sequence- specific primers, and zygosity was assessed using microsatellite markers. Genetic and environmental variance components were assessed with the program Mx, using data from this and previous studies of twins with AS. Results Six of 8 monozygotic (MZ) twin pairs were disease concordant, compared with 4 of 15 B27-positive dizygotic (DZ) twin pairs (27%) and 4 of 32 DZ twin pairs overall (12.5%). Nonsignificant increases in similarity with regard to age at disease onset and all of the disease severity scores assessed were noted in disease-concordant MZ twins compared with concordant DZ twins. HLA-B27 and B60 were associated with the disease in probands, and the rate of disease concordance was significantly increased among DZ twin pairs in which the co- twin was positive for both B27 and DR1. Additive genetic effects were estimated to contribute 97% of the population variance. Conclusion Susceptibility to AS is largely genetically determined, and the environmental trigger for the disease is probably ubiquitous. HLA-B27 accounts for a minority of the overall genetic susceptibility to AS.

A linkage study across the T cell receptor A and T cell receptor B loci in families with rheumatoid arthritis

Objective. To examine whether the T cell receptor (TCR) A or TCRB loci exhibit linkage with disease in multiplex rheumatoid arthritis (RA) families. Methods. A linkage study was performed in 184 RA families from the UK Arthritis and Rheumatism Council Repository, each containing at least 1 affected sibpair. The microsatellites D14S50, TCRA, and D14S64 spanning the TCRA locus and D7S509, Vβ6.7, and D7S688 spanning the TCRB locus were used as DNA markers. The subjects were genotyped using a semiautomated polymerase chain reaction-based method. Two-point and multipoint linkage analyses were performed. Results. Nonparametric single-marker likelihood odds (LOD) scores were 0.49 (P = 0.07) for D14S50, 0.65 (P = 0.04) for TCRA, 0.07 (P = 0.29) for D14S64, 0.01 (P = 0.43) for D7S509, 0.0 (P = 0.50) for Vβ6.7, and 0.0 (P = 0.50) for D7S688. By multipoint analysis, there was no evidence of linkage at TCRB (LOD score 0), and the maximum LOD score at the TCRA locus was 0.37 (at D14S50). The presence of a susceptibility locus (LOD score < -2.0) was excluded, with lambda ≤ 1.8 at TCRA and ≤1.4 at TCRB. Conclusion. These linkage studies provide no significant evidence of a major germline-encoded TCRA or TCRB component of susceptibility to RA.

Is disease severity in ankylosing spondylitis genetically determined?

Objective. To assess the role of genes and the environment in determining the severity of ankylosing spondylitis. Methods: One hundred seventy-three families with >1 case of ankylosing spondylitis were recruited (120 affected sibling pairs, 26 affected parent-child pairs, 20 families with both first- and second-degree relatives affected, and 7 families with only second-degree relatives affected), comprising a total of 384 affected individuals. Disease severity was assessed by the Bath Ankylosing Spondylitis Disease Activity Index (BASDAI) and functional impairment was determined using the Bath Ankylosing Spondylitis Functional Index (BASFI). Disease duration and age at onset were also studied. Variance-components modeling was used to determine the genetic and environmental components Contributing to familiality of the traits examined, and complex segregation analysis was performed to assess different disease models. Results. Both the disease activity and functional capacity as assessed by the BASDAI and the BASFI, respectively, were found to be highly familial (BASDAI familiality 0.51 [P = 10-4], BASFI familiality 0,68 [P = 3 × 10-7]). No significant shared environmental component was demonstrated to be associated with either the BASDAI or the BASFI. Including age at disease onset and duration of disease as covariates made no difference in the heritability assessments. A strong correlation was noted between the BASDAI and the BASFI (genetic correlation 0.9), suggesting the presence of shared determinants of these 2 measures. However, there was significant residual heritability for each measure independent of the other (BASFI residual heritability 0.48, BASDAI 0,36), perhaps indicating that not all genes influencing disease activity influence chronicity. No significant heritability of age at disease onset was found (heritability 0.18; P = 0.2). Segregation studies suggested the presence of a single major gene influencing the BASDAI and the BASFI. Conclusion. This study demonstrates a major genetic contribution to disease severity in ankylosing spondylitis. As with susceptibility to ankylosing spondylitis, shared environmental factors play little role in determining the disease severity.

Whole-genome linkage analysis of rheumatoid arthritis susceptibility loci in 252 affected sibling pairs in the United Kingdom

Objective. To undertake a systematic wholegenome screen to identify regions exhibiting genetic linkage to rheumatoid arthritis (RA). Methods. Two hundred fifty-two RA-affected sibling pairs from 182 UK families were genotyped using 365 highly informative microsatellite markers. Microsatellite genotyping was performed using fluorescent polymerase chain reaction primers and semiautomated DNA sequencing technology. Linkage analysis was undertaken using MAPMAKER/SIBS for single-point and multipoint analysis. Results. Significant linkage (maximum logarithm of odds score 4.7 [P = 0.000003] at marker D6S276, 1 cM from HLA-DRB1) was identified around the major histocompatibility complex (MHC) region on chromosome 6. Suggestive linkage (P < 7.4 × 10-4) was identified on chromosome 6q by single- and multipoint analysis. Ten other sites of nominal linkage (P < 0.05) were identified on chromosomes 3p, 4q, 7p, 2 regions of 10q, 2 regions of 14q, 16p, 21q, and Xq by single-point analysis and on 3 sites (1q, 14q, and 14q) by multipoint analysis. Conclusion. Linkage to the MHC region was confirmed. Eleven non-HLA regions demonstrated evidence of suggestive or nominal linkage, but none reached the genome-wide threshold for significant linkage (P = 2.2 × 10-5). Results of previous genome screens have suggested that 6 of these regions may be involved in RA susceptibility.

Dissection of class III major histocompatibility complex haplotypes associated with rheumatoid arthritis

Objective. We have previously identified a single-nucleotide polymorphism (SNP) haplotype involving the lymphotoxin α (LTA) and tumor necrosis factor (TNF) loci (termed haplotype LTA-TNF2) on chromosome 6 that shows differential association with rheumatoid arthritis (RA) on HLA-DRB1*0404 and *0401 haplotypes, suggesting the presence of additional non-HLA-DRB1 RA susceptibility genes on these haplotypes. To refine this association, we performed a case-control association study using both SNPs and microsatellite markers in haplotypes matched either for HLA-DRB1*0404 or for HLA-DRB1*0401. Methods. Fourteen SNPs lying between HLA-DRB1 and LTA were genotyped in 87 DRB1*04-positive families. High-density microsatellite typing was performed using 24 markers spanning 2,500 kb centered around the TNF gene in 305 DRB1*0401 or *0404 cases and 400 DRB1*0401 or *0404 controls. Single-marker, 2-marker, and 3-marker minihaplotypes were constructed and their frequencies compared between the DRB1*0401 and DRB1*0404 matched case and control haplotypes. Results. Marked preservation of major histocompatibility complex haplotypes was seen, with chromosomes carrying LTA-TNF2 and either DRB1*0401 or DRB1*0404 both carrying an identical SNP haplotype across the 1-Mb region between TNF and HLA-DRB1. Using microsatellite markers, we observed two 3-marker minihaplotypes that were significantly overrepresented in the DRB1*0404 case haplotypes (P = 0.00024 and P = 0.00097). Conclusion. The presence of a single extended SNP haplotype between LTA-TNF2 and both DRB1*0401 and DRB1*0404 is evidence against this region harboring the genetic effects in linkage disequillbrium with LTA-TNF2. Two RA-associated haplotypes on the background of DRB1*0404 were identified in a 126-kb region surrounding and centromeric to the TNF locus.

Investigation of the role of ANKH in ankylosing spondylitis

Objective The ank/ank mouse develops a phenotype similar to ankylosing spondylitis (AS) in humans. ANKH, the human homolog of the mutated gene in the ank/ank mouse, has been implicated in familial autosomal-dominant chondrocalcinosis and autosomal-dominant craniometaphyseal dysplasia. This study was undertaken to investigate the role of ANKH in susceptibility to and clinical manifestations of AS. Methods Sequence variants were identified by genomic sequencing of the 12 ANKH exons and their flanking splice sites in 48 AS patients; variants were then screened in 233 patients and 478 controls. Linkage to the ANKH locus was assessed in 185 affected-sibling-pair families. Results Five single-nucleotide polymorphisms were identified within the coding region and flanking splice sites. No association between either susceptibility to AS or its clinical manifestations and these novel polymorphisms, or between disease susceptibility and 3 known promoter variants, was seen. No linkage between the ANKH locus and AS was observed. Multipoint exclusion mapping rejected the hypothesis of a locus of a magnitude λ≥1.4 (logarithm of odds score <-2) (equivalent to a genetic contribution of >10% to the AS sibling recurrence risk ratio) within this area contributing to AS. Conclusion These findings indicate that ANKH is not significantly involved in susceptibility to or clinical manifestations of AS.

Excessive bone formation in a mouse model of ankylosing spondylitis is associated with decreases in Wnt pathway inhibitors

Introduction: Ankylosing spondylitis (AS) is unique in its pathology where inflammation commences at the entheses before progressing to an osteoproliferative phenotype generating excessive bone formation that can result in joint fusion. The underlying mechanisms of this progression are poorly understood. Recent work has suggested that changes in Wnt signalling, a key bone regulatory pathway, may contribute to joint ankylosis in AS. Using the proteoglycan-induced spondylitis (PGISp) mouse model which displays spondylitis and eventual joint fusion following an initial inflammatory stimulus, we have characterised the structural and molecular changes that underlie disease progression. Methods: PGISp mice were characterised 12 weeks after initiation of inflammation using histology, immunohistochemistry (IHC) and expression profiling. Results: Inflammation initiated at the periphery of the intervertebral discs progressing to disc destruction followed by massively excessive cartilage and bone matrix formation, as demonstrated by toluidine blue staining and IHC for collagen type I and osteocalcin, leading to syndesmophyte formation. Expression levels of DKK1 and SOST, Wnt signalling inhibitors highly expressed in joints, were reduced by 49% and 63% respectively in the spine PGISp compared with control mice (P < 0.05) with SOST inhibition confirmed by IHC. Microarray profiling showed genes involved in inflammation and immune-regulation were altered. Further, a number of genes specifically involved in bone regulation including other members of the Wnt pathway were also dysregulated. Conclusions: This study implicates the Wnt pathway as a likely mediator of the mechanism by which inflammation induces bony ankylosis in spondyloarthritis, raising the potential that therapies targeting this pathway may be effective in preventing this process.

Whole blood transcriptional profiling in ankylosing spondylitis identifies novel candidate genes that might contribute to the inflammatory and tissue-destructive disease aspects

Introduction: A number of genetic-association studies have identified genes contributing to ankylosing spondylitis (AS) susceptibility but such approaches provide little information as to the gene activity changes occurring during the disease process. Transcriptional profiling generates a 'snapshot' of the sampled cells' activity and thus can provide insights into the molecular processes driving the disease process. We undertook a whole-genome microarray approach to identify candidate genes associated with AS and validated these gene-expression changes in a larger sample cohort. Methods: A total of 18 active AS patients, classified according to the New York criteria, and 18 gender- and age-matched controls were profiled using Illumina HT-12 whole-genome expression BeadChips which carry cDNAs for 48,000 genes and transcripts. Class comparison analysis identified a number of differentially expressed candidate genes. These candidate genes were then validated in a larger cohort using qPCR-based TaqMan low density arrays (TLDAs). Results: A total of 239 probes corresponding to 221 genes were identified as being significantly different between patients and controls with a P-value <0.0005 (80% confidence level of false discovery rate). Forty-seven genes were then selected for validation studies, using the TLDAs. Thirteen of these genes were validated in the second patient cohort with 12 downregulated 1.3- to 2-fold and only 1 upregulated (1.6-fold). Among a number of identified genes with well-documented inflammatory roles we also validated genes that might be of great interest to the understanding of AS progression such as SPOCK2 (osteonectin) and EP300, which modulate cartilage and bone metabolism. Conclusions: We have validated a gene expression signature for AS from whole blood and identified strong candidate genes that may play roles in both the inflammatory and joint destruction aspects of the disease.

Mapping genes for osteoporosis-Old dogs and new tricks

In stark contrast to its horticultural origins, modern genetics is an extremely technology-driven field. Almost all the major advances in the field over the past 20 years have followed technological developments that have permitted change in study designs. The development of PCR in the 1980s led to RFLP mapping of monogenic diseases. The development of fluorescent-tagged genotyping methods led to linkage mapping approaches for common diseases that dominated the 1990s. The development of microarray SNP genotyping has led to the genome-wide association study era of the new millennium. And now the development of next-generation sequencing technologies is about to open up a new era of gene-mapping, enabling many potential new study designs. This review aims to present the strengths and weaknesses of the current approaches, and present some new ideas about gene-mapping approaches that are likely to advance our knowledge of the genes involved in heritable bone traits such as bone mineral density (BMD) and fracture.

Genetic studies of common rheumatological diseases

Editorial

The dominant 55kDa allergen of the subtropical Bahia grass (Paspalum notatum) pollen is a group 13 pollen allergen, Pas n 13

Bahia grass, Paspalum notatum, is an important pollen allergen source with a long season of pollination and wide distribution in subtropical and temperate regions. We aimed to characterize the 55. kDa allergen of Bahia grass pollen (BaGP) and ascertain its clinical importance. BaGP extract was separated by 2D-PAGE and immunoblotted with serum IgE of a grass pollen-allergic patient. The amino-terminal protein sequence of the predominant allergen isoform at 55. kDa had similarity with the group 13 allergens of Timothy grass and maize pollen, Phl p 13 and Zea m 13. Four sequences obtained by rapid amplification of the allergen cDNA ends represented multiple isoforms of Pas n 13. The predicted full length cDNA for Pas n 13 encoded a 423 amino acid glycoprotein including a signal peptide of 28 residues and with a predicted pI of 7.0. Tandem mass spectrometry of tryptic peptides of 2D gel spots identified peptides specific to the deduced amino acid sequence for each of the four Pas n 13 cDNA, representing 47% of the predicted mature protein sequence of Pas n 13. There was 80.6% and 72.6% amino acid identity with Zea m 13 and Phl p 13, respectively. Reactivity with a Phl p 13-specific monoclonal antibody AF6 supported designation of this allergen as Pas n 13. The allergen was purified from BaGP extract by ammonium sulphate precipitation, hydrophobic interaction and size exclusion chromatography. Purified Pas n 13 reacted with serum IgE of 34 of 71 (48%) grass pollen-allergic patients and specifically inhibited IgE reactivity with the 55. kDa band of BaGP for two grass pollen-allergic donors. Four isoforms of Pas n 13 from pI 6.3-7.8 had IgE-reactivity with grass pollen allergic sera. The allergenic activity of purified Pas n 13 was demonstrated by activation of basophils from whole blood of three grass pollen-allergic donors tested but not control donors. Pas n 13 is thus a clinically relevant pollen allergen of the subtropical Bahia grass likely to be important in eliciting seasonal allergic rhinitis and asthma in grass pollen-allergic patients.