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.

Interaction between ERAP1 and HLA-B27 in ankylosing spondylitis implicates peptide handling in the mechanism for HLA-B27 in disease susceptibility

Ankylosing spondylitis is a common form of inflammatory arthritis predominantly affecting the spine and pelvis that occurs in approximately 5 out of 1,000 adults of European descent. Here we report the identification of three variants in the RUNX3, LTBR-TNFRSF1A and IL12B regions convincingly associated with ankylosing spondylitis (P < 5 × 10-8 in the combined discovery and replication datasets) and a further four loci at PTGER4, TBKBP1, ANTXR2 and CARD9 that show strong association across all our datasets (P < 5 × 10-6 overall, with support in each of the three datasets studied). We also show that polymorphisms of ERAP1, which encodes an endoplasmic reticulum aminopeptidase involved in peptide trimming before HLA class I presentation, only affect ankylosing spondylitis risk in HLA-B27-positive individuals. These findings provide strong evidence that HLA-B27 operates in ankylosing spondylitis through a mechanism involving aberrant processing of antigenic peptides.

Genetics and genomics of ankylosing spondylitis

Ankylosing spondylitis (AS) is a common, highly heritable arthropathy, the pathogenesis of which is poorly understood. The mechanism by which the main gene for the disease, HLA-B27, leads to AS is unknown. Genetic and genomic studies have demonstrated involvement of the interleukin-23 (IL-23) signaling pathway in AS, a finding which has stimulated much new research into the disease and has led to therapeutic trials. Several other genes and genetic regions, including further major histocompatibility complex (MHC) and non-MHC loci, have been shown to be involved in the disease, but it is not clear yet how they actually induce the condition. These findings have shown that there is a strong genetic overlap between AS and Crohn's disease in particular, although there are also major differences in the genes involved in the two conditions, presumably explaining their different presentations. Genomic and proteomic studies are in an early phase but have potential both as diagnostic/prognostic tools and as a further hypothesis-free tool to investigate AS pathogenesis. Given the slow progress in studying the mechanism of association of HLA-B27 with AS, these may prove to be more fruitful approaches to investigating the pathogenesis of the disease. © 2009 John Wiley & Sons A/S.

Genomics of ankylosing spondylitis

Ankylosing spondylitis (AS) is the prototypic and most prevalent and debilitating spondyloarthropathy, a group of arthritides where the spine and pelvis are specifically targeted. Unlike many other forms of arthritis in which joint damage is mediated through tissue destruction, in AS uncontrolled bone formation occurs, frequently resulting in joint fusion and consequently significant disability. It is estimated that there are 2.4 million spondyloarthritis sufferers in the U.S., twice as many as rheumatoid arthritis. The pathogenesis of AS is very poorly understood and both genetics and gene expression profiling approaches have been utilized to elucidate the underlying mechanisms and pathways that drive the disease. Using powerful genome-wide association study approaches a number of candidate genes have been found to be associated with AS. However, although such approaches can identify genes that can contribute to the disease process, they do not inform us of the actual changes in gene/cell activity at any point in the disease process. Expression profiling allows us to take a "snapshot" of cellular activity and what gene activity changes are underlying those changes. A number of expression profiling studies have been undertaken in AS, looking at both circulating cells and tissues from affected joints. The results to date have been somewhat disappointing with little consensus on gene activity changes due to the low power of the studies undertaken. Some more recent better powered studies have identified diagnostic expression profiles that do point to a possible role for expression profiling in early AS diagnosis. Future studies will require collaborative approaches to target specific disease stages and sites with larger numbers of samples.

Genome-wide association study of ankylosing spondylitis identifies non-MHC susceptibility loci

To identify susceptibility loci for ankylosing spondylitis, we undertook a genome-wide association study in 2,053 unrelated ankylosing spondylitis cases among people of European descent and 5,140 ethnically matched controls, with replication in an independent cohort of 898 ankylosing spondylitis cases and 1,518 controls. Cases were genotyped with Illumina HumHap370 genotyping chips. In addition to strong association with the major histocompatibility complex (MHC; P 10 800), we found association with SNPs in two gene deserts at 2p15 (rs10865331; combined P = 1.9 × 10 19) and 21q22 (rs2242944; P = 8.3 × 10 20), as well as in the genes ANTXR2 (rs4333130; P = 9.3 × 10 8) and IL1R2 (rs2310173; P = 4.8 × 10 7). We also replicated previously reported associations at IL23R (rs11209026; P = 9.1 × 10 14) and ERAP1 (rs27434; P = 5.3 × 10 12). This study reports four genetic loci associated with ankylosing spondylitis risk and identifies a major role for the interleukin (IL)-23 and IL-1 cytokine pathways in disease susceptibility. © 2010 Nature America, Inc. All rights reserved.

Inflammation-driven bone formation in a mouse model of ankylosing spondylitis: Sequential not parallel processes

Background Ankylosing spondylitis (AS) is an immune-mediated arthritis particularly targeting the spine and pelvis and is characterised by inflammation, osteoproliferation and frequently ankylosis. Current treatments that predominately target inflammatory pathways have disappointing efficacy in slowing disease progression. Thus, a better understanding of the causal association and pathological progression from inflammation to bone formation, particularly whether inflammation directly initiates osteoproliferation, is required. Methods The proteoglycan-induced spondylitis (PGISp) mouse model of AS was used to histopathologically map the progressive axial disease events, assess molecular changes during disease progression and define disease progression using unbiased clustering of semi-quantitative histology. PGISp mice were followed over a 24-week time course. Spinal disease was assessed using a novel semi-quantitative histological scoring system that independently evaluated the breadth of pathological features associated with PGISp axial disease, including inflammation, joint destruction and excessive tissue formation (osteoproliferation). Matrix components were identified using immunohistochemistry. Results Disease initiated with inflammation at the periphery of the intervertebral disc (IVD) adjacent to the longitudinal ligament, reminiscent of enthesitis, and was associated with upregulated tumor necrosis factor and metalloproteinases. After a lag phase, established inflammation was temporospatially associated with destruction of IVDs, cartilage and bone. At later time points, advanced disease was characterised by substantially reduced inflammation, excessive tissue formation and ectopic chondrocyte expansion. These distinct features differentiated affected mice into early, intermediate and advanced disease stages. Excessive tissue formation was observed in vertebral joints only if the IVD was destroyed as a consequence of the early inflammation. Ectopic excessive tissue was predominantly chondroidal with chondrocyte-like cells embedded within collagen type II- and X-rich matrix. This corresponded with upregulation of mRNA for cartilage markers Col2a1, sox9 and Comp. Osteophytes, though infrequent, were more prevalent in later disease. Conclusions The inflammation-driven IVD destruction was shown to be a prerequisite for axial disease progression to osteoproliferation in the PGISp mouse. Osteoproliferation led to vertebral body deformity and fusion but was never seen concurrent with persistent inflammation, suggesting a sequential process. The findings support that early intervention with anti-inflammatory therapies will be needed to limit destructive processes and consequently prevent progression of AS.

Disease-associated polymorphisms in ERAP1 do not alter endoplasmic reticulum stress in patients with ankylosing spondylitis

The mechanism by which human leukocyte antigen B27 (HLA-B27) contributes to ankylosing spondylitis (AS) remains unclear. Genetic studies demonstrate that association with and interaction between polymorphisms of endoplasmic reticulum aminopeptidase 1 (ERAP1) and HLA-B27 influence the risk of AS. It has been hypothesised that ERAP1-mediated HLA-B27 misfolding increases endoplasmic reticulum (ER) stress, driving an interleukin (IL) 23-dependent, pro-inflammatory immune response. We tested the hypothesis that AS-risk ERAP1 variants increase ER-stress and concomitant pro-inflammatory cytokine production in HLA-B27 + but not HLA-B27-AS patients or controls. Forty-nine AS cases and 22 healthy controls were grouped according to HLA-B27 status and AS-associated ERAP1 rs30187 genotypes: HLA-B27 + ERAP1 risk, HLA-B27 + ERAP1 protective, HLA-B27-ERAP1 risk and HLA-B27-ERAP1 protective. Expression levels of ER-stress markers GRP78 (8 kDa glucose-regulated protein), CHOP (C/EBP-homologous protein) and inflammatory cytokines were determined in peripheral blood mononuclear cell and ileal biopsies. We found no differences in ER-stress gene expression between HLA-B27 + and HLA-B27-cases or healthy controls, or between cases or controls stratified by carriage of ERAP1 risk or protective alleles in the presence or absence of HLA-B27. No differences were observed between expression of IL17A or TNF (tumour necrosis factor) in HLA-B27 + ERAP1 risk, HLA-B27 + ERAP1 protective and HLA-B27-ERAP1 protective cases. These data demonstrate that aberrant ERAP1 activity and HLA-B27 carriage does not alter ER-stress levels in AS, suggesting that ERAP1 and HLA-B27 may influence disease susceptibility through other mechanisms. © 2015 Macmillan Publishers Limited.

ERAP2 functional knockout in humans does not alter surface heavy chains or HLA-B27, inflammatory cytokines or endoplasmic reticulum stress markers

Introduction Single nucleotide polymorphisms in ERAP2 are strongly associated with ankylosing spondylitis (AS). One AS-associated single nucleotide polymorphism, rs2248374, causes a truncated ERAP2 protein that is degraded by nonsense-mediated decay. Approximately 25% of the populations of European ancestry are therefore natural ERAP2 knockouts. We investigated the effect of this associated variant on HLA class I allele presentation, surface heavy chains, endoplasmic reticulum (ER) stress markers and cytokine gene transcription in AS. Methods Patients with AS and healthy controls with either AA or GG homozygous status for rs2248374 were studied. Antibodies to CD14, CD19-ECD, HLA-A-B-C, Valpha7.2, CD161, anti-HC10 and anti-HLA-B27 were used to analyse peripheral blood mononuclear cells. Expression levels of ER stress markers (GRP78 and CHOP) and proinflammatory genes (tumour necrosis factor (TNF), IL6, IL17 and IL22) were assessed by qPCR. Results There was no significant difference in HLAclass I allele presentation or major histocompatibility class I heavy chains or ER stress markers GRP78 and CHOP or proinflammatory gene expression between genotypes for rs2248374 either between cases, between cases and controls, and between controls. Discussion Large differences were not seen in HLAB27 expression or cytokine levels between subjects with and without ERAP2 in AS cases and controls. This suggests that ERAP2 is more likely to influence AS risk through other mechanisms.

A brilliant breakthrough in OI type V

Interferon-induced transmembrane protein 5 or bone-restricted i ifitm-like gene (Bril) was first identified as a bone gene in 2008, although no in vivo role was identified at that time. A role in human bone has now been demonstrated with a number of recent studies identifying a single point mutation in Bril as the causative mutation in osteogenesis imperfecta type V (OI type V). Such a discovery suggests a key role for Bril in skeletal regulation, and the completely novel nature of the gene raises the possibility of a new regulatory pathway in bone. Furthermore, the phenotype of OI type V has unique and quite divergent features compared with other forms of OI involving defects in collagen biology. Currently it appears that the underlying genetic defect in OI type V may be unrelated to collagen regulation, which also raises interesting questions about the classification of this form of OI. This review will discuss current knowledge of OI type V, the function of Bril, and the implications of this recent discovery.

Whole-exome re-sequencing in a family quartet identifies POP1 mutations as the cause of a novel skeletal dysplasia

Recent advances in DNA sequencing have enabled mapping of genes for monogenic traits in families with small pedigrees and even in unrelated cases. We report the identification of disease-causing mutations in a rare, severe, skeletal dysplasia, studying a family of two healthy unrelated parents and two affected children using whole-exome sequencing. The two affected daughters have clinical and radiographic features suggestive of anauxetic dysplasia (OMIM 607095), a rare form of dwarfism caused by mutations of RMRP. However, mutations of RMRP were excluded in this family by direct sequencing. Our studies identified two novel compound heterozygous loss-of-function mutations in POP1, which encodes a core component of the RNase mitochondrial RNA processing (RNase MRP) complex that directly interacts with the RMRP RNA domains that are affected in anauxetic dysplasia. We demonstrate that these mutations impair the integrity and activity of this complex and that they impair cell proliferation, providing likely molecular and cellular mechanisms by which POP1 mutations cause this severe skeletal dysplasia. © 2011 Glazov et al.

Computer planning of stereotactic iodine-125 seed brachytherapy for recurrent malignant gliomas

At St Thomas' Hospital, we have developed a computer program on a Titan graphics supercomputer to plan the stereotactic implantation of iodine-125 seeds for the palliative treatment of recurrent malignant gliomas. Use of the Gill-Thomas-Cosman relocatable frame allows planning and surgery to be carried out at different hospitals on different days. Stereotactic computed tomography (CT) and positron emission tomography (PET) scans are performed and the images transferred to the planning computer. The head, tumour and frame fiducials are outlined on the relevant images, and a three-dimensional model generated. Structures which could interfere with the surgery or radiotherapy, such as major vessels, shunt tubing etc., can also be outlined and included in the display. Catheter target and entry points are set using a three-dimensional cursor controlled by a set of dials attached to the computer. The program calculates and displays the radiation dose distribution within the target volume for various catheter and seed arrangements. The CT co-ordinates of the fiducial rods are used to convert catheter co-ordinates from CT space to frame space and to calculate the catheter insertion angles and depths. The surgically implanted catheters are after-loaded the next day and the seeds left in place for between 4 and 6 days, giving a nominal dose of 50 Gy to the edge of the target volume. 25 patients have been treated so far.

Years lived with disability (YLDs) for 1160 sequelae of 289 diseases and injuries 1990–2010 : a systematic analysis for the Global Burden of Disease Study 2010

Background Non-fatal health outcomes from diseases and injuries are a crucial consideration in the promotion and monitoring of individual and population health. The Global Burden of Disease (GBD) studies done in 1990 and 2000 have been the only studies to quantify non-fatal health outcomes across an exhaustive set of disorders at the global and regional level. Neither effort quantified uncertainty in prevalence or years lived with disability (YLDs). Methods Of the 291 diseases and injuries in the GBD cause list, 289 cause disability. For 1160 sequelae of the 289 diseases and injuries, we undertook a systematic analysis of prevalence, incidence, remission, duration, and excess mortality. Sources included published studies, case notification, population-based cancer registries, other disease registries, antenatal clinic serosurveillance, hospital discharge data, ambulatory care data, household surveys, other surveys, and cohort studies. For most sequelae, we used a Bayesian meta-regression method, DisMod-MR, designed to address key limitations in descriptive epidemiological data, including missing data, inconsistency, and large methodological variation between data sources. For some disorders, we used natural history models, geospatial models, back-calculation models (models calculating incidence from population mortality rates and case fatality), or registration completeness models (models adjusting for incomplete registration with health-system access and other covariates). Disability weights for 220 unique health states were used to capture the severity of health loss. YLDs by cause at age, sex, country, and year levels were adjusted for comorbidity with simulation methods. We included uncertainty estimates at all stages of the analysis. Findings Global prevalence for all ages combined in 2010 across the 1160 sequelae ranged from fewer than one case per 1 million people to 350 000 cases per 1 million people. Prevalence and severity of health loss were weakly correlated (correlation coefficient −0·37). In 2010, there were 777 million YLDs from all causes, up from 583 million in 1990. The main contributors to global YLDs were mental and behavioural disorders, musculoskeletal disorders, and diabetes or endocrine diseases. The leading specific causes of YLDs were much the same in 2010 as they were in 1990: low back pain, major depressive disorder, iron-deficiency anaemia, neck pain, chronic obstructive pulmonary disease, anxiety disorders, migraine, diabetes, and falls. Age-specific prevalence of YLDs increased with age in all regions and has decreased slightly from 1990 to 2010. Regional patterns of the leading causes of YLDs were more similar compared with years of life lost due to premature mortality. Neglected tropical diseases, HIV/AIDS, tuberculosis, malaria, and anaemia were important causes of YLDs in sub-Saharan Africa. Interpretation Rates of YLDs per 100 000 people have remained largely constant over time but rise steadily with age. Population growth and ageing have increased YLD numbers and crude rates over the past two decades. Prevalences of the most common causes of YLDs, such as mental and behavioural disorders and musculoskeletal disorders, have not decreased. Health systems will need to address the needs of the rising numbers of individuals with a range of disorders that largely cause disability but not mortality. Quantification of the burden of non-fatal health outcomes will be crucial to understand how well health systems are responding to these challenges. Effective and affordable strategies to deal with this rising burden are an urgent priority for health systems in most parts of the world. Funding Bill & Melinda Gates Foundation.

Integrated genomic characterization of endometrial carcinoma

We performed an integrated genomic, transcriptomic and proteomic characterization of 373 endometrial carcinomas using array- and sequencing-based technologies. Uterine serous tumours and ∼25% of high-grade endometrioid tumours had extensive copy number alterations, few DNA methylation changes, low oestrogen receptor/progesterone receptor levels, and frequent TP53 mutations. Most endometrioid tumours had few copy number alterations or TP53 mutations, but frequent mutations in PTEN, CTNNB1, PIK3CA, ARID1A and KRAS and novel mutations in the SWI/SNF chromatin remodelling complex gene ARID5B. A subset of endometrioid tumours that we identified had a markedly increased transversion mutation frequency and newly identified hotspot mutations in POLE. Our results classified endometrial cancers into four categories: POLE ultramutated, microsatellite instability hypermutated, copy-number low, and copy-number high. Uterine serous carcinomas share genomic features with ovarian serous and basal-like breast carcinomas. We demonstrated that the genomic features of endometrial carcinomas permit a reclassification that may affect post-surgical adjuvant treatment for women with aggressive tumours.

Analysis of albumin-associated peptides and proteins from ovarian cancer patients

Albumin binds low–molecular-weight molecules, including proteins and peptides, which then acquire its longer half-life, thereby protecting the bound species from kidney clearance. We developed an experimental method to isolate albumin in its native state and to then identify [mass spectrometry (MS) sequencing] the corresponding bound low–molecular-weight molecules. We used this method to analyze pooled sera from a human disease study set (high-risk persons without cancer, n= 40; stage I ovarian cancer, n = 30; stage III ovarian cancer, n = 40) to demonstrate the feasibility of this approach as a discovery method. Methods Albumin was isolated by solid-phase affinity capture under native binding and washing conditions. Captured albumin-associated proteins and peptides were separated by gel electrophoresis and subjected to iterative MS sequencing by microcapillary reversed-phase tandem MS. Selected albumin-bound protein fragments were confirmed in human sera by Western blotting and immunocompetition. Results In total, 1208 individual protein sequences were predicted from all 3 pools. The predicted sequences were largely fragments derived from proteins with diverse biological functions. More than one third of these fragments were identified by multiple peptide sequences, and more than one half of the identified species were in vivo cleavage products of parent proteins. An estimated 700 serum peptides or proteins were predicted that had not been reported in previous serum databases. Several proteolytic fragments of larger molecules that may be cancer-related were confirmed immunologically in blood by Western blotting and peptide immunocompetition. BRCA2, a 390-kDa low-abundance nuclear protein linked to cancer susceptibility, was represented in sera as a series of specific fragments bound to albumin. Conclusion Carrier-protein harvesting provides a rich source of candidate peptides and proteins with potential diverse tissue and cellular origins that may reflect important disease-related information.