Association of ERAP1, but not IL23R, with ankylosing spondylitis in a Han Chinese population

Objective The results of a recent genome-wide association study have shown that ERAP1 and IL23R are associated with ankylosing spondylitis (AS) in Caucasian populations from North America and the UK. Based on these findings, we undertook the current study to investigate whether single-nucleotide polymorphisms (SNPs) covering the genes ERAP1 and IL23R are associated with AS in a Han Chinese population. Methods A case-control study was performed in Han Chinese patients with AS (n = 527) and controls (n = 945) from Shanghai and Nanjing. All patients met the modified New York criteria for AS. The Sequenom iPlex platform was used to genotype cases and controls for 21 tag SNPs covering IL23R and 38 tag SNPs covering ERAP1. Statistical analysis was performed using the Cochran-Armitage test for trend. Results Multiple SNPs in ERAP1 were significantly associated with AS (for rs27980, P = 0.0048; for rs7711564, P = 0.0081). However, no association was observed between IL23R and AS (for all SNPs, P > 0.1). The nonsynonymous SNP in IL23R, rs11209026, widely thought to be the primary AS-associated SNP in IL23R in Europeans, was found not to be polymorphic in Chinese. Conclusion Our results demonstrate that genetic polymorphisms in ERAP1 are associated with AS in Han Chinese, suggesting a common pathogenic mechanism for the disease in Chinese and Caucasian populations, and that IL23R is not associated with AS in Chinese, indicating a difference in the mechanism of disease pathogenesis between Chinese and Caucasian populations. This may result from the fact that rs11209026, the nonsynonymous SNP in IL23R, is not polymorphic in Chinese patients, providing further evidence that rs11209026 is the key polymorphism associated with AS (and likely inflammatory bowel disease and psoriasis) in this gene.

Enrichment of circulating interleukin-17-secreting interleukin-23 receptor-positive γ/δ T cells in patients with active ankylosing spondylitis

Objective Ankylosing spondylitis (AS) is a common inflammatory arthritis affecting primarily the axial skeleton. IL23R is genetically associated with AS. This study was undertaken to investigate and characterize the role of interleukin-23 (IL-23) signaling in AS pathogenesis. Methods The study population consisted of patients with active AS (n = 17), patients with psoriatic arthritis (n = 8), patients with rheumatoid arthritis, (n = 9), and healthy subjects (n = 20). IL-23 receptor (IL-23R) expression in T cells was determined in each subject group, and expression levels were compared. Results The proportion of IL-23R-expressing T cells in the periphery was 2-fold higher in AS patients than in healthy controls, specifically driven by a 3-fold increase in IL-23R-positive γ/δ T cells in AS patients. The proportions of CD4+ and CD8+ cells that were positive for IL-17 were unchanged. This increased IL-23R expression on γ/δ T cells was also associated with enhanced IL-17 secretion, with no observable IL-17 production from IL-23R-negative γ/δ T cells in AS patients. Furthermore, γ/δ T cells from AS patients were heavily skewed toward IL-17 production in response to stimulation with IL-23 and/or anti-CD3/CD28. Conclusion Recently, mouse models have shown IL-17-secreting γ/δ T cells to be pathogenic in infection and autoimmunity. Our data provide the first description of a potentially pathogenic role of these cells in a human autoimmune disease. Since IL-23 is a maturation and growth factor for IL-17-producing cells, increased IL-23R expression may regulate the function of this putative pathogenic γ/δ T cell population.

Interleukin-23 mediates the intestinal response to microbial β-1,3-glucan and the development of spondyloarthritis pathology in SKG mice

Objective Spondyloarthritides (SpA) occur in 1% of the population and include ankylosing spondylitis (AS) and arthropathy of inflammatory bowel disease (IBD), with characteristic spondylitis, arthritis, enthesitis, and IBD. Genetic studies implicate interleukin-23 (IL-23) receptor signaling in the development of SpA and IBD, and IL-23 overexpression in mice is sufficient for enthesitis, driven by entheseal-resident T cells. However, in genetically prone individuals, it is not clear where IL-23 is produced and how it drives the SpA syndrome, including IBD or subclinical gut inflammation of AS. Moreover, it is unclear why specific tissue involvement varies between patients with SpA. We undertook this study to determine the location of IL-23 production and its role in SpA pathogenesis in BALB/c ZAP-70W163C-mutant (SKG) mice injected intraperitoneally with β-1,3-glucan (curdlan). Methods Eight weeks after curdlan injection in wild-type or IL-17A-/- SKG or BALB/c mice, pathology was scored in tissue sections. Mice were treated with anti-IL-23 or anti-IL-22. Cytokine production and endoplasmic reticulum (ER) stress were determined in affected organs. Results In curdlan-treated SKG mice, arthritis, enthesitis, and ileitis were IL-23 dependent. Enthesitis was specifically dependent on IL-17A and IL-22. IL-23 was induced in the ileum, where it amplified ER stress, goblet cell dysfunction, and proinflammatory cytokine production. IL-17A was pathogenic, while IL-22 was protective against ileitis. IL-22+CD3- innate-like cells were increased in lamina propria mononuclear cells of ileitis-resistant BALB/c mice, which developed ileitis after curdlan injection and anti-IL-22. Conclusion In response to systemic β-1,3-glucan, intestinal IL-23 provokes local mucosal dysregulation and cytokines driving the SpA syndrome, including IL-17/IL-22-dependent enthesitis. Innate IL-22 production promotes ileal tolerance.

Serum levels of soluble receptor for advanced glycation end products and of S100 proteins are associated with inflammatory, autoantibody, and classical risk markers of joint and vascular damage in rheumatoid arthritis

Introduction: The receptor for advanced glycation end products (RAGE) is a member of the immunoglobulin superfamily of cell surface receptor molecules. High concentrations of three of its putative proinflammatory ligands, S100A8/A9 complex (calprotectin), S100A8, and S100A12, are found in rheumatoid arthritis (RA) serum and synovial fluid. In contrast, soluble RAGE (sRAGE) may prevent proinflammatory effects by acting as a decoy. This study evaluated the serum levels of S100A9, S100A8, S100A12 and sRAGE in RA patients, to determine their relationship to inflammation and joint and vascular damage. Methods: Serum sRAGE, S100A9, S100A8 and S100A12 levels from 138 patients with established RA and 44 healthy controls were measured by ELISA and compared by unpaired t test. In RA patients, associations with disease activity and severity variables were analyzed by simple and multiple linear regressions. Results: Serum S100A9, S100A8 and S100A12 levels were correlated in RA patients. S100A9 levels were associated with body mass index (BMI), and with serum levels of S100A8 and S100A12. S100A8 levels were associated with serum levels of S100A9, presence of anti-citrullinated peptide antibodies (ACPA), and rheumatoid factor (RF). S100A12 levels were associated with presence of ACPA, history of diabetes, and serum S100A9 levels. sRAGE levels were negatively associated with serum levels of C-reactive protein (CRP) and high-density lipoprotein (HDL), history of vasculitis, and the presence of the RAGE 82Ser polymorphism. Conclusions: sRAGE and S100 proteins were associated not just with RA inflammation and autoantibody production, but also with classical vascular risk factors for end-organ damage. Consistent with its role as a RAGE decoy molecule, sRAGE had the opposite effects to S100 proteins in that S100 proteins were associated with autoantibodies and vascular risk, whereas sRAGE was associated with protection against joint and vascular damage. These data suggest that RAGE activity influences co-development of joint and vascular disease in rheumatoid arthritis patients.

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.

The enigma of IgE+ B-cell memory in human subjects

Our understanding of the origin and fate of the IgE-switched B cell has been markedly improved by studies in mouse models. The immediate precursor of the IgE-switched B cell is either a relatively naive nonswitched B cell or a mature IgG-switched B cell. These 2 routes are referred to as the direct and indirect pathways, respectively. IgE responses derived from each pathway differ significantly, largely reflecting the difference in time spent in a germinal center and thus time for clonal expansion, somatic hypermutation, affinity maturation, and acquisition of a memory phenotype. The clinical and therapeutic implications for IgE responses in human subjects are still a matter of debate, largely because the immunization procedures used in the animal models are significantly different from classical atopic sensitization to allergens from pollen and mites. On the basis of the limited information available, it seems likely that these atopic IgE responses are characterized by a relatively low IgG/IgE ratio, low B-cell memory, and modest affinity maturation, which fits well with the direct switching pathway. It is still unresolved how the IgE response evolves to cover a wide epitope repertoire involving many epitopes per allergen, as well as many different allergens from a single allergen source. © 2013 American Academy of Allergy, Asthma & Immunology.

Mixture models of nucleotide sequence evolution that account for heterogeneity in the substitution process across sites and across lineages

Molecular phylogenetic studies of homologous sequences of nucleotides often assume that the underlying evolutionary process was globally stationary, reversible, and homogeneous (SRH), and that a model of evolution with one or more site-specific and time-reversible rate matrices (e.g., the GTR rate matrix) is enough to accurately model the evolution of data over the whole tree. However, an increasing body of data suggests that evolution under these conditions is an exception, rather than the norm. To address this issue, several non-SRH models of molecular evolution have been proposed, but they either ignore heterogeneity in the substitution process across sites (HAS) or assume it can be modeled accurately using the distribution. As an alternative to these models of evolution, we introduce a family of mixture models that approximate HAS without the assumption of an underlying predefined statistical distribution. This family of mixture models is combined with non-SRH models of evolution that account for heterogeneity in the substitution process across lineages (HAL). We also present two algorithms for searching model space and identifying an optimal model of evolution that is less likely to over- or underparameterize the data. The performance of the two new algorithms was evaluated using alignments of nucleotides with 10 000 sites simulated under complex non-SRH conditions on a 25-tipped tree. The algorithms were found to be very successful, identifying the correct HAL model with a 75% success rate (the average success rate for assigning rate matrices to the tree's 48 edges was 99.25%) and, for the correct HAL model, identifying the correct HAS model with a 98% success rate. Finally, parameter estimates obtained under the correct HAL-HAS model were found to be accurate and precise. The merits of our new algorithms were illustrated with an analysis of 42 337 second codon sites extracted from a concatenation of 106 alignments of orthologous genes encoded by the nuclear genomes of Saccharomyces cerevisiae, S. paradoxus, S. mikatae, S. kudriavzevii, S. castellii, S. kluyveri, S. bayanus, and Candida albicans. Our results show that second codon sites in the ancestral genome of these species contained 49.1% invariable sites, 39.6% variable sites belonging to one rate category (V1), and 11.3% variable sites belonging to a second rate category (V2). The ancestral nucleotide content was found to differ markedly across these three sets of sites, and the evolutionary processes operating at the variable sites were found to be non-SRH and best modeled by a combination of eight edge-specific rate matrices (four for V1 and four for V2). The number of substitutions per site at the variable sites also differed markedly, with sites belonging to V1 evolving slower than those belonging to V2 along the lineages separating the seven species of Saccharomyces. Finally, sites belonging to V1 appeared to have ceased evolving along the lineages separating S. cerevisiae, S. paradoxus, S. mikatae, S. kudriavzevii, and S. bayanus, implying that they might have become so selectively constrained that they could be considered invariable sites in these species.

Estimating the potential impact of entry fees for marine parks on dive tourism in South East Asia

Marine reserves are increasingly being established as a mechanism to protect marine biodiversity and sensitive habitats. As well as providing conservation benefits, marine reserves provide benefits to recreational scuba divers who dive within the reserve, as well as to recreational and commercial fishers outside the reserve through spill-over effects. To ensure benefits are being realised, management of marine reserves requires ongoing monitoring and surveillance. These are not costless, and many marine reserve managers impose an entry fee. In some countries, dive tourism is major income source to coastal industries, and a concern is that high entry fees may dissuade divers. In this paper, the price elasticity of demand for dive tourism in three countries in South East Asia – Indonesia, Thailand and Malaysia – is estimated using a travel-cost model. From the model, the total non-market use value associated with diving in the area is estimated to be in the order of US$4.5 billion a year. The price elasticity of demand in the region is highly inelastic, such that increasing the cost of diving through a management levy would have little impact on total diver numbers.

Modeling approaches for pesticide exposure assessment in rice paddies

Pesticide use in paddy rice production may contribute to adverse ecological effects in surface waters. Risk assessments conducted for regulatory purposes depend on the use of simulation models to determine predicted environment concentrations (PEC) of pesticides. Often tiered approaches are used, in which assessments at lower tiers are based on relatively simple models with conservative scenarios, while those at higher tiers have more realistic representations of physical and biochemical processes. This chapter reviews models commonly used for predicting the environmental fate of pesticides in rice paddies. Theoretical considerations, unique features, and applications are discussed. This review is expected to provide information to guide model selection for pesticide registration, regulation, and mitigation in rice production areas.

Environmental health risk assessment of dioxin in foods at the two most severe dioxin hot spots in Vietnam

Background Bien Hoa and Da Nang airbases were bulk storages for Agent Orange during the Vietnam War and currently are the two most severe dioxin hot spots. Objectives This study assesses the health risk of exposure to dioxin through foods for local residents living in seven wards surrounding these airbases. Methods This study follows the Australian Environmental Health Risk Assessment Framework to assess the health risk of exposure to dioxin in foods. Forty-six pooled samples of commonly consumed local foods were collected and analyzed for dioxin/furans. A food frequency and Knowledge–Attitude–Practice survey was also undertaken at 1000 local households, various stakeholders were involved and related publications were reviewed. Results Total dioxin/furan concentrations in samples of local “high-risk” foods (e.g. free range chicken meat and eggs, ducks, freshwater fish, snail and beef) ranged from 3.8 pg TEQ/g to 95 pg TEQ/g, while in “low-risk” foods (e.g. caged chicken meat and eggs, seafoods, pork, leafy vegetables, fruits, and rice) concentrations ranged from 0.03 pg TEQ/g to 6.1 pg TEQ/g. Estimated daily intake of dioxin if people who did not consume local high risk foods ranged from 3.2 pg TEQ/kg bw/day to 6.2 pg TEQ/kg bw/day (Bien Hoa) and from 1.2 pg TEQ/kg bw/day to 4.3 pg TEQ/kg bw/day (Da Nang). Consumption of local high risk foods resulted in extremely high dioxin daily intakes (60.4–102.8 pg TEQ/kg bw/day in Bien Hoa; 27.0–148.0 pg TEQ/kg bw/day in Da Nang). Conclusions Consumption of local “high-risk” foods increases dioxin daily intakes far above the WHO recommended TDI (1–4 pg TEQ/kg bw/day). Practicing appropriate preventive measures is necessary to significantly reduce exposure and health risk.

Dioxin risk reduction programme at the most severe dioxin hot spot in Vietnam: Residual knowledge, attitude and practices five years post intervention

Bien Hoa Airbase was one of the bulk storage and supply facilities for defoliants during the Vietnam War. Environmental and biological samples taken around the airbase have elevated levels of dioxin. In 2007, a pre-intervention knowledge, attitude and practice (KAP) survey of local residents living in Trung Dung and Tan Phong wards was undertaken regarding appropriate strategies to reduce dioxin exposure. A risk reduction programme was implemented in 2008 and post-intervention KAP surveys were undertaken in 2009 and 2013 to evaluate the longer term impacts. Quantitative assessment was undertaken via a KAP survey in 2013 among 600 local residents randomly selected from the two intervention wards and one control ward (Buu Long). Eight in-depth interviews and two focus group discussions were also undertaken for qualitative assessment. Most programme activities had ceased and dioxin risk communication activities had not been integrated into local routine health education programmes; however, main results generally remained and were better than that in Buu Long. In total, 48.2% of households undertook measures to prevent exposure, higher than those in pre- and post-intervention surveys (25.8% and 39.7%) and the control ward (7.7%). Migration and the sensitive nature of dioxin issues were the main challenges for the programme's sustainability

Discovery of candidate serum proteomic and metabolomic biomarkers in ankylosing spondylitis

Ankylosing Spondylitis (AS) is a common inflammatory rheumatic disease with a predilection for the axial skeleton, affecting 0.2% of the population. Current diagnostic criteria rely on a composite of clinical and radiological changes, with a mean time to diagnosis of 5 to 10 years. In this study we employed nano liquid-chromatography mass spectrometry analysis to detect and quantify proteins and small compounds including endogenous peptides and metabolites in serum from 18 AS patients and nine healthy individuals. We identified a total of 316 proteins in serum, of which 22 showed significant up- or down-regulation (p < 0.05) in AS patients. Receiver operating characteristic analysis of combined levels of serum amyloid P component and inter-α-trypsin inhibitor heavy chain 1 revealed high diagnostic value for Ankylosing Spondylitis (area under the curve = 0.98). We also depleted individual sera of proteins to analyze endogenous peptides and metabolic compounds. We detected more than 7000 molecular features in patients and healthy individuals. Quantitative MS analysis revealed compound profiles that correlate with the clinical assessment of disease activity. One molecular feature identified as a Vitamin D3 metabolite-(23S,25R)-25-hydroxyvitamin D3 26,23-peroxylactone-was down-regulated in AS. The ratio of this vitamin D metabolite versus vitamin D binding protein serum levels was also altered in AS as compared with controls. These changes may contribute to pathological skeletal changes in AS. Our study is the first example of an integration of proteomic and metabolomic techniques to find new biomarker candidates for the diagnosis of Ankylosing Spondylitis

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.

Genetic analyses in a sample of individuals with high or low BMD shows association with multiple Wnt pathway genes

INTRODUCTION Although the high heritability of BMD variation has long been established, few genes have been conclusively shown to affect the variation of BMD in the general population. Extreme truncate selection has been proposed as a more powerful alternative to unselected cohort designs in quantitative trait association studies. We sought to test these theoretical predictions in studies of the bone densitometry measures BMD, BMC, and femoral neck area, by investigating their association with members of the Wnt pathway, some of which have previously been shown to be associated with BMD in much larger cohorts, in a moderate-sized extreme truncate selected cohort (absolute value BMD Z-scores = 1.5-4.0; n = 344). MATERIALS AND METHODS Ninety-six tag-single nucleotide polymorphism (SNPs) lying in 13 Wnt signaling pathway genes were selected to tag common genetic variation (minor allele frequency [MAF] > 5% with an r(2) > 0.8) within 5 kb of all exons of 13 Wnt signaling pathway genes. The genes studied included LRP1, LRP5, LRP6, Wnt3a, Wnt7b, Wnt10b, SFRP1, SFRP2, DKK1, DKK2, FZD7, WISP3, and SOST. Three hundred forty-four cases with either high or low BMD were genotyped by Illumina Goldengate microarray SNP genotyping methods. Association was tested either by Cochrane-Armitage test for dichotomous variables or by linear regression for quantitative traits. RESULTS Strong association was shown with LRP5, polymorphisms of which have previously been shown to influence total hip BMD (minimum p = 0.0006). In addition, polymorphisms of the Wnt antagonist, SFRP1, were significantly associated with BMD and BMC (minimum p = 0.00042). Previously reported associations of LRP1, LRP6, and SOST with BMD were confirmed. Two other Wnt pathway genes, Wnt3a and DKK2, also showed nominal association with BMD. CONCLUSIONS This study shows that polymorphisms of multiple members of the Wnt pathway are associated with BMD variation. Furthermore, this study shows in a practical trial that study designs involving extreme truncate selection and moderate sample sizes can robustly identify genes of relevant effect sizes involved in BMD variation in the general population. This has implications for the design of future genome-wide studies of quantitative bone phenotypes relevant to osteoporosis.