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...

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 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.

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.

The chromosome 16q region associated with ankylosing spondylitis includes the candidate gene tumour necrosis factor receptor type 1-associated death domain (TRADD)

Objective: To replicate and refine the reported association of ankylosing spondylitis (AS) with two nonsynonymous single nucleotide polymorphisms (nsSNPs) on chromosome 16q22.1. Methods: Firstly, 730 independent UK patients with AS were genotyped for rs9939768 and rs6979 and allele frequencies were compared with 2879 previously typed historic disease controls. Secondly, the two data sets were combined in meta-analyses. Finally, 5 tagging SNPs, located between rs9939768 and rs6979, were analysed in 1604 cases and 1020 controls. Results: The association of rs6979 with AS was replicated, p=0.03, OR=1.14 (95% CI 1.01 to 1.28), and a trend for association with rs9939768 detected, p=0.06, OR=1.25 (95% CI 0.99 to 1.57). Meta-analyses revealed association of both SNPs with AS, p=0.0008, OR=1.31 (95% CI 1.12 to 1.54) and p=0.0009, OR=1.15 (95% CI 1.06 to 1.23) for rs9939768 and rs6979, respectively. New associations with rs9033 and rs868213 (p=0.00002, OR=1.23 (95% CI 1.12 to 1.36) and p=0.00002 OR=1.45 (95% CI 1.22 to 1.72), respectively, were identified. Conclusions: The region on chromosome 16 that has been replicated in the present work is interesting as the highly plausible candidate gene, tumour necrosis factor receptor type 1 (TNFR1)-associated death domain (TRADD), is located between rs9033 and rs868213. It will require additional work to identify the primary genetic association(s) with AS.

The effect of HLA-DR genes on susceptibility to and severity of ankylosing spondylitis

Objective. To analyze the effect of HLA-DR genes on susceptibility to and severity of ankylosing spondylitis (AS). Methods. Three hundred sixty- three white British AS patients were studied; 149 were carefully assessed for a range of clinical manifestations, and disease severity was assessed using a structured questionnaire. Limited HLA class I typing and complete HLA-DR typing were performed using DNA-based methods. HLA data from 13,634 healthy white British bone marrow donors were used for comparison. Results. A significant association between DR1 and AS was found, independent of HLA-B27 (overall odds ratio [OR] 1.4, 95% confidence interval [95% CI] 1.1-1.8, P = 0.02; relative risk [RR] 2.7, 95% CI 1.5-4.8, P = 6 x 10-4 among homozygotes; RR 2.1, 95% CI 1.5-2.8, P = 5 x 10-6 among heterozygotes). A large but weakly significant association between DR8 and AS was noted, particularly among DR8 homozygotes (RR 6.8, 95% CI 1.6-29.2, P = 0.01 among homozygotes; RR 1.6, 95% CI 1.0-2.7, P = 0.07 among heterozygotes). A negative association with DR12 (OR 0.22, 95% CI 0.09-0.5, P = 0.001) was noted. HLA-DR7 was associated with younger age at onset of disease (mean age at onset 18 years for DR7-positive patients and 23 years for DR7-negative patients; Z score 3.21, P = 0.001). No other HLA class I or class H associations with disease severity or with different clinical manifestations of AS were found. Conclusion. The results of this study suggest that HLA-DR genes may have a weak effect on susceptibility to AS independent of HLA-B27, but do not support suggestions that they affect disease severity or different clinical manifestations.

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.

IgE-Associated IGHV Genes from Venom and Peanut Allergic Individuals Lack Mutational Evidence of Antigen Selection

Antigen selection of B cells within the germinal center reaction generally leads to the accumulation of replacement mutations in the complementarity-determining regions (CDRs) of immunoglobulin genes. Studies of mutations in IgE-associated VDJ gene sequences have cast doubt on the role of antigen selection in the evolution of the human IgE response, and it may be that selection for high affinity antibodies is a feature of some but not all allergic diseases. The severity of IgE-mediated anaphylaxis is such that it could result from higher affinity IgE antibodies. We therefore investigated IGHV mutations in IgE-associated sequences derived from ten individuals with a history of anaphylactic reactions to bee or wasp venom or peanut allergens. IgG sequences, which more certainly experience antigen selection, served as a control dataset. A total of 6025 unique IgE and 5396 unique IgG sequences were generated using high throughput 454 pyrosequencing. The proportion of replacement mutations seen in the CDRs of the IgG dataset was significantly higher than that of the IgE dataset, and the IgE sequences showed little evidence of antigen selection. To exclude the possibility that 454 errors had compromised analysis, rigorous filtering of the datasets led to datasets of 90 core IgE sequences and 411 IgG sequences. These sequences were present as both forward and reverse reads, and so were most unlikely to include sequencing errors. The filtered datasets confirmed that antigen selection plays a greater role in the evolution of IgG sequences than of IgE sequences derived from the study participants.

The transcriptome of the bowhead whale Balaena mysticetus reveals adaptations of the longest‐lived mammal

Mammals vary dramatically in lifespan, by at least two-orders of magnitude, but the molecular basis for this difference remains largely unknown. The bowhead whale Balaena mysticetus is the longest-lived mammal known, with an estimated maximal lifespan in excess of two hundred years. It is also one of the two largest animals and the most cold-adapted baleen whale species. Here, we report the first genome-wide gene expression analyses of the bowhead whale, based on the de novo assembly of its transcriptome. Bowhead whale or cetacean-specific changes in gene expression were identified in the liver, kidney and heart, and complemented with analyses of positively selected genes. Changes associated with altered insulin signaling and other gene expression patterns could help explain the remarkable longevity of bowhead whales as well as their adaptation to a lipid-rich diet. The data also reveal parallels in candidate longevity adaptations of the bowhead whale, naked mole rat and Brandt's bat. The bowhead whale transcriptome is a valuable resource for the study of this remarkable animal, including the evolution of longevity and its important correlates such as resistance to cancer and other diseases.

Non-major-histocompatibility-complex genetics of ankylosing spondylitis

There is strong evidence from twin and family studies indicating that a substantial proportion of the heritability of susceptibility to ankylosing spondylitis (AS) and its clinical manifestations is encoded by non-major-histocompatibility-complex genes. Efforts to identify these genes have included genomewide linkage studies and candidate gene association studies. One region, the interleukin (IL)-1 gene complex on chromosome 2, has been repeatedly associated with AS in both Caucasians and Asians. It is likely that more than one gene in this complex is involved in AS, with the strongest evidence to date implicating IL-1A. Identifying the genes underlying other linkage regions has been difficult due to the lack of obvious candidates and the low power of most studies to date to identify genes of the small to moderate magnitude that are likely to be involved. The field is moving towards genomewide association analysis, involving much larger datasets of unrelated cases and controls. Early successes using this approach in other diseases indicates that it is likely to identify genes in common diseases like AS, but there remains the risk that the common-variant, common-disease hypothesis will not hold true in AS. Nonetheless, it is appropriate for the field to be cautiously optimistic that the next few years will bring great advances in our understanding of the genetics of this condition.

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.

Association of IL23R and ERAP1 genes with ankylosing spondylitis in a Portuguese population

Objective: Association between ankylosing spondylitis (AS) and two genes, ERAP1 and IL23R, has recently been reported in North American and British populations. The population attributable risk fraction for ERAP1 in this study was 25%, and for IL23R, 9%. Confirmation of these findings to ERAP1 in other ethnic groups has not yet been demonstrated. We sought to test the association between single nucleotide polymorphisms (SNPs) in these genes and susceptibility to AS among a Portuguese population. We also investigated the role of these genes in clinical manifestations of AS, including age of symptom onset, the Bath Ankylosing Spondylitis Disease Activity, Metrology and Functional Indices, and the modified Stoke Ankylosing Spondylitis Spinal Score. Methods: The study was conducted on 358 AS cases and 285 ethnically matched Portuguese healthy controls. AS was defined according to the modified New York Criteria. Genotyping of IL23R and ERAP1 allelic variants was carried out with TaqMan allelic discrimination assays. Association analysis was performed using the Cochrane-Armitage and linear regression tests of genotypes as implemented in PLINK for dichotomous and quantitative variables respectively. A meta-analysis for Portuguese and previously published Spanish IL23R data was performed using the StatsDirect® Statistical tools, by fixed and random effects models. Results: A total of 14 nsSNPs markers (8 for IL23R, 5 for ERAPl, 1 for LN-PEP) were analysed. Three markers (2 for IL23R and 1 for ERAP1) showed significant single-locus disease associations, confirming that the association of these genes with AS in the Portuguese population. The strongest associated SNP in IL23R was rs1004819 (OR=1.4, p=0.0049), and in ERAP1 was rs30187 (OR=1.26, p=0.035). The population attributable risk fractions in the Portuguese population for these SNPs are 11% and 9.7% respectively. No association was seen with any SNP in LN-PEP, which flanks ERAP1 and was associated with AS in the British population. No association was seen with clinical manifestations of AS. Conclusions: These results show that IL23R and ERAP1 genes are also associated with susceptibility to AS in the Portuguese population, and that they contribute a significant proportion of the population risk for this disease.

Cuff tear arthropathy: Evidence of functional variation in pyrophosphate metabolism genes

We investigated the role of two genes, ANKH and TNAP, in patients with cuff tear arthropathy. These genes encode proteins which regulate the extracellular concentration of inorganic pyrophosphate, fluctuations of which can lead to calcium crystal formation. Variants were detected by direct sequencing of DNA and their frequencies compared with healthy controls. The effect of variants on protein function was further studied by in vitro approaches. Variant genotypes were observed more frequently in the cases when compared with controls in ANKH (45% and 20%) and TNAP (32% and 9%). Variants in ANKH altered inorganic pyrophosphate (PPi) concentrations in transfected human chondrocytes. There was a higher mean serum concentration of TNAP detected in female patients compared with normal ranges. Cuff tear arthropathy is associated with variants in ANKH and TNAP that alter extracellular inorganic pyrophosphate concentrations causing calcium crystal deposition. This supports a theory that genetic variants predispose patients to primary crystal deposition which when combined with a massive rotator cuff tear leads to the development of arthritis.

Genetic association and gene expression studies suggest that genetic variants in the SYNE1 and TNF genes are related to menstrual migraine

BACKGROUND: Menstrual migraine (MM) encompasses pure menstrual migraine (PMM) and menstrually-related migraine (MRM). This study was aimed at investigating genetic variants that are potentially related to MM, specifically undertaking genotyping and mRNA expression analysis of the ESR1, PGR, SYNE1 and TNF genes in MM cases and non-migraine controls. METHODS: A total of 37 variants distributed across 14 genes were genotyped in 437 DNA samples (282 cases and 155 controls). In addition levels of gene expression were determined in 74 cDNA samples (41 cases and 33 controls). Association and correlation analysis were performed using Plink and RStudio. RESULTS: SNPs rs3093664 and rs9371601 in TNF and SYNE1 genes respectively, were significantly associated with migraine in the MM population (p = 0.008; p = 0.009 respectively). Analysis of qPCR results found no significant difference in levels of gene expression between cases and controls. However, we found a significant correlation between the expression of ESR1 and SYNE1, ESR1 and PGR and TNF and SYNE1 in samples taken during the follicular phase of the menstrual cycle. CONCLUSIONS: Our results show that SNPs rs9371601 and rs3093664 in the SYNE1 and TNF genes respectively, are associated with MM. The present study also provides strong evidence to support the correlation of ESR1, PGR, SYNE1 and TNF gene expression in MM.