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.

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.

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.

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.

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.

Mosaic PPM1D mutations are associated with predisposition to breast and ovarian cancer

Improved sequencing technologies offer unprecedented opportunities for investigating the role of rare genetic variation in common disease. However, there are considerable challenges with respect to study design, data analysis and replication. Using pooled next-generation sequencing of 507 genes implicated in the repair of DNA in 1,150 samples, an analytical strategy focused on protein-truncating variants (PTVs) and a large-scale sequencing case-control replication experiment in 13,642 individuals, here we show that rare PTVs in the p53-inducible protein phosphatase PPM1D are associated with predisposition to breast cancer and ovarian cancer. PPM1D PTV mutations were present in 25 out of 7,781 cases versus 1 out of 5,861 controls (P = 1.12 × 10-5), including 18 mutations in 6,912 individuals with breast cancer (P = 2.42 × 10-4) and 12 mutations in 1,121 individuals with ovarian cancer (P = 3.10 × 10-9). Notably, all of the identified PPM1D PTVs were mosaic in lymphocyte DNA and clustered within a 370-base-pair region in the final exon of the gene, carboxy-terminal to the phosphatase catalytic domain. Functional studies demonstrate that the mutations result in enhanced suppression of p53 in response to ionizing radiation exposure, suggesting that the mutant alleles encode hyperactive PPM1D isoforms. Thus, although the mutations cause premature protein truncation, they do not result in the simple loss-of-function effect typically associated with this class of variant, but instead probably have a gain-of-function effect. Our results have implications for the detection and management of breast and ovarian cancer risk. More generally, these data provide new insights into the role of rare and of mosaic genetic variants in common conditions, and the use of sequencing in their identification.

Mapping the regulatory sequences controlling 93 breast cancer-associated miRNA genes leads to the identification of two functional promoters of the Hsa-mir-200b cluster, methylation of which is associated with metastasis or hormone receptor status in advanced breast cancer

MicroRNAs (miRNAs) are small non-coding RNAs of 20 nt in length that are capable of modulating gene expression post-transcriptionally. Although miRNAs have been implicated in cancer, including breast cancer, the regulation of miRNA transcription and the role of defects in this process in cancer is not well understood. In this study we have mapped the promoters of 93 breast cancer-associated miRNAs, and then looked for associations between DNA methylation of 15 of these promoters and miRNA expression in breast cancer cells. The miRNA promoters with clearest association between DNA methylation and expression included a previously described and a novel promoter of the Hsa-mir-200b cluster. The novel promoter of the Hsa-mir-200b cluster, denoted P2, is located 2 kb upstream of the 5′ stemloop and maps within a CpG island. P2 has comparable promoter activity to the previously reported promoter (P1), and is able to drive the expression of miR-200b in its endogenous genomic context. DNA methylation of both P1 and P2 was inversely associated with miR-200b expression in eight out of nine breast cancer cell lines, and in vitro methylation of both promoters repressed their activity in reporter assays. In clinical samples, P1 and P2 were differentially methylated with methylation inversely associated with miR-200b expression. P1 was hypermethylated in metastatic lymph nodes compared with matched primary breast tumours whereas P2 hypermethylation was associated with loss of either oestrogen receptor or progesterone receptor. Hypomethylation of P2 was associated with gain of HER2 and androgen receptor expression. These data suggest an association between miR-200b regulation and breast cancer subtype and a potential use of DNA methylation of miRNA promoters as a component of a suite of breast cancer biomarkers.