High-throughput sequencing of IL23R reveals a low-frequency, nonsynonymous single-nucleotide polymorphism that is associated with ankylosing spondylitis in a Han Chinese population

Objective Ankylosing spondylitis (AS) is a highly heritable common inflammatory arthritis that targets the spine and sacroiliac joints of the pelvis, causing pain and stiffness and leading eventually to joint fusion. Although previous studies have shown a strong association of IL23R with AS in white Europeans, similar studies in East Asian populations have shown no association with common variants of IL23R, suggesting either that IL23R variants have no role or that rare genetic variants contribute. The present study was undertaken to screen IL23R to identify rare variants associated with AS in Han Chinese. Methods A 170-kb region containing IL23R and its flanking regions was sequenced in 50 patients with AS and 50 ethnically matched healthy control subjects from a Han Chinese population. In addition, the 30-kb region of peak association in white Europeans was sequenced in 650 patients with AS and 1,300 healthy controls. Validation genotyping was undertaken in 846 patients with AS and 1,308 healthy controls. Results We identified 1,047 variants, of which 729 were not found in the dbSNP genomic build 130. Several potentially functional rare variants in IL23R were identified, including one nonsynonomous single-nucleotide polymorphism (nsSNP), Gly149Arg (position 67421184 GA on chromosome 1). Validation genotyping showed that the Gly149Arg variant was associated with AS (odds ratio 0.61, P = 0.0054). Conclusion This is the first study to implicate rare IL23R variants in the pathogenesis of AS. The results identified a low-frequency nsSNP with predicted loss-of-function effects that was protectively associated with AS in Han Chinese, suggesting that decreased function of the interleukin-23 (IL-23) receptor protects against AS. These findings further support the notion that IL-23 signaling has an important role in the pathogenesis of AS.

Omission of author name in the article by Davidson et al (Arthritis Rheum, July 2013)

This article corrects: Brief Report: High-Throughput Sequencing of IL23R Reveals a Low-Frequency, Nonsynonymous Single-Nucleotide Polymorphism That Is Associated With Ankylosing Spondylitis in a Han Chinese Population Vol. 65, Issue 7, 1747–1752, Article first published online: 2 JUL 2013

Association of sporadic chondrocalcinosis with a -4-basepair G-to-A transition in the 5′-untranslated region of ANKH that promotes enhanced expression of ANKH protein and excess generation of extracellular inorganic pyrophosphate

Objective Certain mutations in ANKH, which encodes a multiple-pass transmembrane protein that regulates inorganic pyrophosphate (PPi) transport, are linked to autosomal-dominant familial chondrocalcinosis. This study investigated the potential for ANKH sequence variants to promote sporadic chondrocalcinosis. Methods ANKH variants identified by genomic sequencing were screened for association with chondrocalcinosis in 128 patients with severe sporadic chondrocalcinosis or pseudogout and in ethnically matched healthy controls. The effects of specific variants on expression of common markers were evaluated by in vitro transcription/translation. The function of these variants was studied in transfected human immortalized CH-8 articular chondrocytes. Results Sporadic chondrocalcinosis was associated with a G-to-A transition in the ANKH 5′-untranslated region (5′-UTR) at 4 bp upstream of the start codon (in homozygotes of the minor allele, genotype relative risk 6.0, P = 0.0006; overall genotype association P = 0.02). This -4-bp transition, as well as 2 mutations previously linked with familial and sporadic chondrocalcinosis (+14 bp C-to-T and C-terminal GAG deletion, respectively), but not the French familial chondrocalcinosis kindred 143-bp T-to-C mutation, increased reticulocyte ANKH transcription/ANKH translation in vitro. Transfection of complementary DNA for both the wild-type ANKH and the -4-bp ANKH protein variant promoted increased extracellular PPi in CH-8 cells, but unexpectedly, these ANKH mutants had divergent effects on the expression of extracellular PPi and the chondrocyte hypertrophy marker, type X collagen. Conclusion A subset of sporadic chondrocalcinosis appears to be heritable via a -4-bp G-to-A ANKH 5′-UTR transition that up-regulates expression of ANKH and extracellular PPi in chondrocyte cells. Distinct ANKH mutations associated with heritable chondrocalcinosis may promote disease by divergent effects on extracellular PPi and chondrocyte hypertrophy, which is likely to mediate differences in the clinical phenotypes and severity of the disease.

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.

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.

Estimates on the coverage of parameter space using populations of models

In this paper we provide estimates for the coverage of parameter space when using Latin Hypercube Sampling, which forms the basis of building so-called populations of models. The estimates are obtained using combinatorial counting arguments to determine how many trials, k, are needed in order to obtain specified parameter space coverage for a given value of the discretisation size n. In the case of two dimensions, we show that if the ratio (Ø) of trials to discretisation size is greater than 1, then as n becomes moderately large the fractional coverage behaves as 1-exp-ø. We compare these estimates with simulation results obtained from an implementation of Latin Hypercube Sampling using MATLAB.

Catalyst engineering for lithium ion batteries: The catalytic role of Ge in enhancing the electrochemical performance of SnO2(GeO2)0.13/G anodes

The catalytic role of germanium (Ge) was investigated to improve the electrochemical performance of tin dioxide grown on graphene (SnO(2)/G) nanocomposites as an anode material of lithium ion batteries (LIBs). Germanium dioxide (GeO(20) and SnO(2) nanoparticles (<10 nm) were uniformly anchored on the graphene sheets via a simple single-step hydrothermal method. The synthesized SnO(2)(GeO(2))0.13/G nanocomposites can deliver a capacity of 1200 mA h g(-1) at a current density of 100 mA g(-1), which is much higher than the traditional theoretical specific capacity of such nanocomposites (∼ 702 mA h g(-1)). More importantly, the SnO(2)(GeO(2))0.13/G nanocomposites exhibited an improved rate, large current capability (885 mA h g(-1) at a discharge current of 2000 mA g(-1)) and excellent long cycling stability (almost 100% retention after 600 cycles). The enhanced electrochemical performance was attributed to the catalytic effect of Ge, which enabled the reversible reaction of metals (Sn and Ge) to metals oxide (SnO(2) and GeO(2)) during the charge/discharge processes. Our demonstrated approach towards nanocomposite catalyst engineering opens new avenues for next-generation high-performance rechargeable Li-ion batteries anode materials.

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.

The correlation between reading and mathematics ability at age twelve has a substantial genetic component

Dissecting how genetic and environmental influences impact on learning is helpful for maximizing numeracy and literacy. Here we show, using twin and genome-wide analysis, that there is a substantial genetic component to children’s ability in reading and mathematics, and estimate that around one half of the observed correlation in these traits is due to shared genetic effects (so-called Generalist Genes). Thus, our results highlight the potential role of the learning environment in contributing to differences in a child’s cognitive abilities at age twelve.

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.

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

Genome-wide study of familial juvenile hyperuricaemic (gouty) nephropathy (FJHN) indicates a new locus, FJHN3, linked to chromosome 2p22.1-p21

Familial juvenile hyperuricaemic (gouty) nephropathy (FJHN), is an autosomal dominant disease associated with a reduced fractional excretion of urate, and progressive renal failure. FJHN is genetically heterogeneous and due to mutations of three genes: uromodulin (UMOD), renin (REN) and hepatocyte nuclear factor-1beta (HNF-1β) on chromosomes 16p12, 1q32.1, and 17q12, respectively. However, UMOD, REN or HNF-1β mutations are found in only ~45% of FJHN probands, indicating the involvement of other genetic loci in ~55% of probands. To identify other FJHN loci, we performed a single nucleotide polymorphism (SNP)-based genome-wide linkage analysis, in six FJHN families in whom UMOD, HNF-1β and REN mutations had been excluded. Parametric linkage analysis using a 'rare dominant' model established linkage in five of the six FJHN families, with a LOD score >+3, at 0% recombination, between FJHN and SNPs at chromosome 2p22.1-p21. Analysis of individual recombinants in two unrelated affected individuals defined a ~5.5 Mbp interval, flanked telomerically by SNP RS372139 and centromerically by RS896986 that contained the locus, designated FJHN3. The interval contains 28 genes, and DNA sequence analysis of the most likely candidate, solute carrier family 8 member 1 (SLC8A1), did not identify any abnormalities in the FJHN3 probands. FJHN3 is likely located within a ~5.5 Mbp interval on chromosome 2p22.1-p21, and identifying the genetic abnormality will help to further elucidate mechanisms predisposing to gout and renal failure.