The global burden attributable to low bone mineral density

Introduction The Global Burden of Disease Study 2010 estimated the worldwide health burden of 291 diseases and injuries and 67 risk factors by calculating disability-adjusted life years (DALYs). Osteoporosis was not considered as a disease, and bone mineral density (BMD) was analysed as a risk factor for fractures, which formed part of the health burden due to falls. Objectives To calculate (1) the global distribution of BMD, (2) its population attributable fraction (PAF) for fractures and subsequently for falls, and (3) the number of DALYs due to BMD. Methods A systematic review was performed seeking population-based studies in which BMD was measured by dual-energy X-ray absorptiometry at the femoral neck in people aged 50 years and over. Age- and sex-specific mean ± SD BMD values (g/cm2) were extracted from eligible studies. Comparative risk assessment methodology was used to calculate PAFs of BMD for fractures. The theoretical minimum risk exposure distribution was estimated as the age- and sex-specific 90th centile from the Third National Health and Nutrition Examination Survey (NHANES III). Relative risks of fractures were obtained from a previous meta-analysis. Hospital data were used to calculate the fraction of the health burden of falls that was due to fractures. Results Global deaths and DALYs attributable to low BMD increased from 103 000 and 3 125 000 in 1990 to 188 000 and 5 216 000 in 2010, respectively. The percentage of low BMD in the total global burden almost doubled from 1990 (0.12%) to 2010 (0.21%). Around one-third of falls-related deaths were attributable to low BMD. Conclusions Low BMD is responsible for a growing global health burden, only partially representative of the real burden of osteoporosis.

Matters Arising: Zhu et al, "A novel gene variation of TNFα associated with ankylosing spondylitis: A reconfirmed study"

In a recently published study1 involving 79 ankylosing spondylitis (AS) patients and 132 unrelated healthy blood donors, Zhu et al report association of a tumour necrosis factor α (TNFα) single nucleotide polymorphism (SNP), rs1799724, with AS.1 I have many concerns with this paper and its conclusions...

Matters Arising: Role of HLA genes in familial spondyloarthropathy

Said-Nahal and colleagues report an intriguing finding of an association with HLA-DR4 independent of B27 in families with ankylosing spondylitis (AS),1 a finding highlighted by an accompanying editorial.2 The approach of studying B27 positive and B27 negative haplotypes may prove powerful in identifying further cis or trans encoded genes involved in AS. However, the reported association of DR4 with AS is quite a surprising finding given that no difference was noted in B27-DR4 haplotype frequencies in patients and ethnically matched healthy controls. Many previous studies have not reported any such association,3–9 including a similar preliminary study by the same authors.10 Although these studies were mainly case-control studies, population stratification is highly unlikely to cause a false negative finding if the effect size of the reported association with DR4 is as high as Said-Nahal and colleagues describe...

Ankylosing spondylitis in West Africans - evidence for a non-HLA-B27 protective effect

Objective - To determine the prevalence of ankylosing spondylitis in the Fula ethnic group in The Gambia, and relate the disease prevalence to the B27 frequency and subtype distribution of that population. Methods - 215 first degree relatives of 48 B27 positive Fula twin pairs, and 900 adult Fula males were screened for ankylosing spondylitis by clinical and, where appropriate, radiographic means. The B27 prevalence was determined by PCR/sequence specific oligonucleotides on finger prick samples from 100 unrelated Fula, and B27 subtype distribution by SSCP on unrelated B27 positive individuals. This data were then compared with the prevalence of ankylosing spondylitis among B27 positive Caucasians. Results - No case of ankylosing spondylitis was seen. Six per cent of Fula are B27 positive, of which 32% are B*2703 and 68% B*2705. Assuming the penetrance of ankylosing spondylitis in B27 positive Fula is the same as in B27 positive Caucasians, the probability of not observing any cases of ankylosing spondylitis among the Fula examined is remote (P = 6.7 x 10-6). Similarly, the chance of not seeing any cases among those expected to be either B*2705 or B*2703 was small (P = 3.2 x 10-4 for B*2705, and P = 0.02 for B*2703). Conclusions - The risk of developing ankylosing spondylitis in B27 positive Fula is lower than in B27 positive Caucasians. This is not explained by the B27 subtype distribution among Fula, and suggests the presence of some non-B27 protective factor reducing the prevalence of ankylosing spondylitis in this population.

Recurrence risk modelling of the genetic susceptibility to ankylosing spondylitis

Objectives - It has long been suspected that susceptibility to ankylosing spondylitis (AS) is influenced by genes lying distant to the major histocompatibility complex. This study compares genetic models of AS to assess the most likely mode of inheritance, using recurrence risk ratios in relatives of affected subjects. Methods - Recurrence risk ratios in different degrees of relatives were determined using published data from studies specifically designed to address the question. The methods of Risch were used to determine the expected recurrence risk ratios in different degrees of relatives, assuming equal first degree relative recurrence risk between models. Goodness of fit was determined by χ2 comparison of the expected number of affected subjects with the observed number, given equal numbers of each type of relative studied. Results - The recurrence risks in different degrees of relatives were: monozygotic (MZ) twins 63% (17/27), first degree relatives 8.2% (441/5390), second degree relatives 1.0% (8/834), and third degree relatives 0.7% (7/997). Parent-child recurrence risk (7.9%, 37/466) was not significantly different from the sibling recurrence risk (8.2%, 404/4924), excluding a significant dominance genetic component to susceptibility. Poor fitting models included single gene, genetic heterogeneity, additive, two locus multiplicative, and one locus and residual polygenes (χ2 > 32 (two degrees of freedom), p < 10-6 for all models). The best fitting model studied was a five locus model with multiplicative interaction between loci (χ2 = 1.4 (two degrees of freedom), p = 0.5). Oligogenic multiplicative models were the best fitting over a range of population prevalences and first degree recurrence risk rates. Conclusions - This study suggests that of the genetic models tested, the most likely model operating in AS is an oligogenic model with predominantly multiplicative interaction between loci.

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.

Replication of association of IL1 gene complex members with ankylosing spondylitis in Taiwanese Chinese

Objective: To test the association of interleukin 1 (IL1) gene family members with ankylosing spondylitis (AS), previously reported in Europid subjects, in an ethnically remote population. Methods: 200 Taiwanese Chinese AS patients and 200 ethnically matched healthy controls were genotyped for five single nucleotide polymorphisms (SNPs) and the IL1RN.VNTR, markers previously associated with AS. Allele, genotype, and haplotype frequencies were compared between cases and controls. Results: Association of alleles and genotypes of the markers IL1F10.3, IL1RN.4, and IL1RN.VNTR was observed with AS (p<0.05). Haplotypes of pairs of these markers and of the markers IL1RN.6/1 and IL1RN.6/2 were also significantly associated with AS. The strongest associations observed were with the marker IL1RN.4, and with the two-marker haplotype IL1RN.4-IL1RN.VNTR (both p = 0.004). Strong linkage disequilibrium was observed between all marker pairs except those involving IL1B-511 (D′ 0.4 to 0.9, p<0.01). Conclusions: The IL1 gene cluster is associated with AS in Taiwanese Chinese. This finding provides strong statistical support that the previously observed association of this gene cluster with AS is a true positive finding.

Association study of genes related to bone formation and resorption and the extent of radiographic change in ankylosing spondylitis

Objective: To identify genetic associations with severity of radiographic damage in ankylosing spondylitis (AS). Method: We studied 1537 AS cases of European descent; all fulfilled the modified New York Criteria. Radiographic severity was assessed from digitised lateral radiographs of the cervical and lumbar spine using the modified Stoke Ankylosing Spondylitis Spinal Score (mSASSS). A two-phase genotyping design was used. In phase 1, 498 single nucleotide polymorphisms (SNPs) were genotyped in 688 cases; these were selected to capture >90% of the common haplotypic variation in the exons, exon-intron boundaries, and 5 kb flanking DNA in the 5' and 3' UTR of 74 genes involved in anabolic or catabolic bone pathways. In phase 2, 15 SNPs exhibiting p<0.05 were genotyped in a further cohort of 830 AS cases; results were analysed both separately and in combination with the discovery phase data. Association was tested by contingency tables after separating the samples into 'mild' and 'severe' groups, defined as the bottom and top 40% by mSASSS, adjusted for gender and disease duration. Results: Experiment-wise association was observed with the SNP rs8092336 (combined OR 0.32, p=1.2×10-5), which lies within RANK (receptor activator of NF?B), a gene involved in osteoclastogenesis, and in the interaction between T cells and dendritic cells. Association was also found with the SNP rs1236913 in PTGS1 (prostaglandin-endoperoxide synthase 1, cyclooxygenase 1), giving an OR of 0.53 (p=2.6×10-3). There was no observed association between radiographic severity and HLA-B*27. Conclusions: These findings support roles for bone resorption and prostaglandins pathways in the osteoproliferative changes in AS.

Association of variants in MMEL1 and CTLA4 with rheumatoid arthritis in the Han Chinese population

Background: The genome-wide association study era has made great progress in identifying susceptibility genes and genetic loci for rheumatoid arthritis (RA) in populations of White European ancestry. However, few studies have tried to dissect disease aetiopathogenesis in other ethnic populations. Objective: To investigate these associations in the Han Chinese population. Methods: Haplotypes from the HapMap database Chinese population were used to select tag-single-nucleotide polymorphisms (SNPs) (r2 =0.8) across 19 distinct RA genomic regions. A two phase case-control association study was performed, with 169 SNPs genotyped in phase I (n=571 cases, n=880 controls), and 64 SNPs achieving p<0.2 in the first phase being genotyped in phase II (n=464 cases, n=822 controls). Association statistics were calculated using permutation tests both unadjusted and adjusted for the number of markers studied. Results: Robust association was detected for MMEL1 and CTLA4 , and modest association was identified for another six loci: PADI4 , STAT4 , PRDM1 , CDK6 , TRAF1-C5 and KIF5A-PIP4K2C. All three markers genotyped in MMEL1 demonstrated association, with peak signal for rs3890745 (p=2.6×10 -5unadjusted, p=0.003 adjusted, OR=0.79). For CTLA4 , significance was detected for three of five variants showing association, with peak association for marker rs12992492 (p=4.3×10-5 unadjusted, p=0.0021 adjusted, OR=0.77). Lack of association of common variants in PTPN22 with RA in Han Chinese was confirmed. Conclusion: This study identifies MMEL1 and CTLA4 as RA susceptibility genes, provides suggestive evidence of association for a further six loci in the Han Chinese population and confirms lack of PTPN22 association in Asian populations. It also confirms the value of multiethnic population studies to help dissect disease aetiopathogenesis.

Association of STAT3 and TNFRSF1A with ankylosing spondylitis in Han Chinese

Objectives: Recent association studies by the Australo-Anglo-American Spondyloarthritis Consortium (TASC) in Caucasian European populations from Australia, North America and the UK have identified a number of genes as being associated with ankylosing spondylitis (AS). A candidate gene study in a Han Chinese population was performed based on these findings to identify associated genes in this population. Methods: A case-control study was performed in a Han Chinese population of patients with AS (n=775) and controls (n=1587) from Shanghai and Nanjing. All patients met the modified New York criteria for AS. The cases and controls were genotyped for 115 single nucleotide polymorphisms (SNPs) tagging IL23R, ERAP1, STAT3, JAK2, TNFRSF1A and TRADD, as well as other confirmation SNPs from the TASC study, using the Sequenom iPlex and the ABI OpenArray platforms. Statistical analysis of genotyped SNPs was performed using the Cochran - Armitage test for trend and meta-analysis was performed using METAL. SNPs in AS-associated genes in this study were then imputed using MaCH, and association with AS tested by logistic regression. Results: SNPs in TNFRSF1A (rs4149577, p=8.2×10-4), STAT3 (rs2293152, p=0.0015; rs1053005, p=0.017) and ERAP1 (rs27038, p=0.0091; rs27037, p=0.0092) were significantly associated with AS in Han Chinese. Association was also observed between AS and the intergenic region 2p15 (rs10865331, p=0.023). The lack of association between AS and IL23R in Han Chinese was confirmed (all SNPs p>0.1). Conclusions: The study results demonstrate for the first time that genetic polymorphisms in STAT3, TNFRSF1A and 2p15 are associated with AS in Han Chinese, suggesting common pathogenic mechanisms for the disease in Chinese and Caucasian European populations. Furthermore, previous findings demonstrating that ERAP1, but not IL23R, is associated with AS in Chinese patients were confirmed.

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.

Matters Arising: HLA-DR-DQ haplotypes and genotypes in Finnish patients with rheumatoid arthritis

We read with great interest the paper by Laivoranta-Nyman et al.1 They studied Finnish patients with rheumatoid arthritis (RA) and controls, and suggested that there exist susceptibility, neutral and protective HLA-DR-DQ haplotypes that do not match the predictions of current hypotheses for the mechanism of association of the HLA class II region and RA...

Letter: Brain natriuretic peptide is a potentially useful screening tool for the detection of cardiovascular disease in patients with rheumatoid arthritis

Patients with rheumatoid arthritis (RA) have a significantly higher risk of coronary heart disease, despite being less likely to report symptoms of angina, and are more likely to experience unrecognised myocardial infarction and sudden cardiac death than non-RA controls.1 Furthermore, left ventricular diastolic dysfunction has been described in up to 40% of patients with RA.2...

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.

Finnish HLA studies confirm the increased risk conferred by HLA-B27 homozygosity in ankylosing spondylitis

Objective: To determine the influence of HLA-B27 homozygosity and HLA-DRB1 alleles in the susceptibility to, and severity of, ankylosing spondylitis in a Finnish population. Methods: 673 individuals from 261 families with ankylosing spondylitis were genotyped for HLA-DRB1 alleles and HLA-B27 heterozygosity/ homozygosity. The frequencies of HLA-B27 homozygotes in probands from these families were compared with the expected number of HLA-B27 homozygotes in controls under Hardy-Weinberg equilibrium (HWE). The effect of HLA-DRB1 alleles was assessed using a logistic regression procedure conditioned on HLA-B27 and case-control analysis. Results: HLA-B27 was detected in 93% of cases of ankylosing spondylitis. An overrepresentation of HLA-B27 homozygotes was noted in ankylosing spondylitis (11%) compared with the expected number of HLA-B27 homozygotes under HWE (4%) (odds ratio (OR) = 3.3 (95% confidence interval, 1.6 to 6.8), p = 0.002). HLA-B27 homozygosity was marginally associated with reduced BASDAI (HLA-B27 homozygotes, 4.5 (1.6); HLA-B27 heterozygotes, 5.4 (1.8) (mean (SD)), p = 0.05). Acute anterior uveitis (AAU) was present in significantly more HLA-B27 positive cases (50%) than HLA-B27 negative cases (16%) (OR = 5.4 (1.7 to 17), p<0.004). HLA-B27 positive cases had a lower average age of symptom onset (26.7 (8.0) years) compared with HLA-B27 negative cases (35.7 (11.2) years) (p<0.0001). Conclusions: HLA-627 homozygosity is associated with a moderately increased risk of ankylosing spondylitis compared with HLA-β27 heterozygosity. HLA-B27 positive cases had an earlier age of onset of ankylosing spondylitis than HLA-B27 negative cases and were more likely to develop AAU. HLA-DRB1 alleles may influence the age of symptom onset of ankylosing spondylitis.