Genome-Wide Association Study and Gene Expression Analysis Identifies CD84 as a Predictor of Response to Etanercept Therapy in Rheumatoid Arthritis

Anti-tumor necrosis factor alpha (anti-TNF) biologic therapy is a widely used treatment for rheumatoid arthritis (RA). It is unknown why some RA patients fail to respond adequately to anti-TNF therapy, which limits the development of clinical biomarkers to predict response or new drugs to target refractory cases. To understand the biological basis of response to anti-TNF therapy, we conducted a genome-wide association study (GWAS) meta-analysis of more than 2 million common variants in 2,706 RA patients from 13 different collections. Patients were treated with one of three anti-TNF medications: etanercept (n = 733), infliximab (n = 894), or adalimumab (n = 1,071). We identified a SNP (rs6427528) at the 1q23 locus that was associated with change in disease activity score (ΔDAS) in the etanercept subset of patients (P = 8×10-8), but not in the infliximab or adalimumab subsets (P>0.05). The SNP is predicted to disrupt transcription factor binding site motifs in the 3′ UTR of an immune-related gene, CD84, and the allele associated with better response to etanercept was associated with higher CD84 gene expression in peripheral blood mononuclear cells (P = 1×10-11 in 228 non-RA patients and P = 0.004 in 132 RA patients). Consistent with the genetic findings, higher CD84 gene expression correlated with lower cross-sectional DAS (P = 0.02, n = 210) and showed a non-significant trend for better ΔDAS in a subset of RA patients with gene expression data (n = 31, etanercept-treated). A small, multi-ethnic replication showed a non-significant trend towards an association among etanercept-treated RA patients of Portuguese ancestry (n = 139, P = 0.4), but no association among patients of Japanese ancestry (n = 151, P = 0.8). Our study demonstrates that an allele associated with response to etanercept therapy is also associated with CD84 gene expression, and further that CD84 expression correlates with disease activity. These findings support a model in which CD84 genotypes and/or expression may serve as a useful biomarker for response to etanercept treatment in RA patients of European ancestry. © 2013 Cui et al

A Replication Study of the Association between Rheumatoid Arthritis and Deletion of the Late Cornified Envelope Genes LCE3B and LCE3C

OBJECTIVE: Two recent studies, in a Spanish and a Chinese population, point to an association between rheumatoid arthritis (RA) risk and the deletion of the Late Cornified Envelope (LCE) 3B and 3C genes (LCE3C_LCE3B-del), a known risk factor for psoriasis. We aimed to replicate these studies in a large Dutch cohort. METHODS: 1039 RA cases and 759 controls were genotyped for LCE3C_LCE3B-del. Association analysis was performed for the complete cohort and after stratification for the serologic markers anti-cyclic citrullinated peptide and rheumatoid factor. A meta-analysis was performed combining our data with the Spanish and Chinese datasets, resulting in an analysis including 2466 RA cases and 2438 controls. RESULTS: In the Dutch cohort we did not observe a significant association of LCE3C_LCE3B-del (p = 0.093) with RA risk. A stratified analysis for the serologic positive and negative group did not show an association between the genetic variant and disease risk, either. The meta-analysis, however, confirmed a significant association (p<0.0001, OR = 1.31, 95% confidence interval 1.16-1.47). CONCLUSION: Our meta-analysis confirms the association of the LCE3 deletion with RA, suggesting that LCE3C_LCE3B-del is a common risk factor for (auto)immune diseases

Validation Study of Existing Gene Expression Signatures for Anti-TNF Treatment in Patients with Rheumatoid Arthritis

So far, there are no means of identifying rheumatoid arthritis (RA) patients who will fail to respond to tumour necrosis factor blocking agents (anti-TNF), prior to treatment. We set out to validate eight previously reported gene expression signatures predicting therapy outcome. Genome-wide expression profiling using Affymetrix GeneChip Exon 1.0 ST arrays was performed on RNA isolated from whole blood of 42 RA patients starting treatment with infliximab or adalimumab. Clinical response according to EULAR criteria was determined at week 14 of therapy. Genes that have been reported to be associated with anti-TNF treatment were extracted from our dataset. K-means partition clustering was performed to assess the predictive value of the gene-sets. We performed a hypothesis-driven analysis of the dataset using eight existing gene sets predictive of anti-TNF treatment outcome. The set that performed best reached a sensitivity of 71% and a specificity of 61%, for classifying the patients in the current study. We successfully validated one of eight previously reported predictive expression profile. This replicated expression signature is a good starting point for developing a prediction model for anti-TNF treatment outcome that can be used in a daily clinical setting. Our results confirm that gene expression profiling prior to treatment is a useful tool to predict anti-TNF (non) response

T-cell exhaustion, co-stimulation and clinical outcome in autoimmunity and infection.

The clinical course of autoimmune and infectious disease varies greatly, even between individuals with the same condition. An understanding of the molecular basis for this heterogeneity could lead to significant improvements in both monitoring and treatment. During chronic infection the process of T-cell exhaustion inhibits the immune response, facilitating viral persistence. Here we show that a transcriptional signature reflecting CD8 T-cell exhaustion is associated with poor clearance of chronic viral infection, but conversely predicts better prognosis in multiple autoimmune diseases. The development of CD8 T-cell exhaustion during chronic infection is driven both by persistence of antigen and by a lack of accessory 'help' signals. In autoimmunity, we find that where evidence of CD4 T-cell co-stimulation is pronounced, that of CD8 T-cell exhaustion is reduced. We can reproduce the exhaustion signature by modifying the balance of persistent stimulation of T-cell antigen receptors and specific CD2-induced co-stimulation provided to human CD8 T cells in vitro, suggesting that each process plays a role in dictating outcome in autoimmune disease. The 'non-exhausted' T-cell state driven by CD2-induced co-stimulation is reduced by signals through the exhaustion-associated inhibitory receptor PD-1, suggesting that induction of exhaustion may be a therapeutic strategy in autoimmune and inflammatory disease. Using expression of optimal surrogate markers of co-stimulation/exhaustion signatures in independent data sets, we confirm an association with good clinical outcome or response to therapy in infection (hepatitis C virus) and vaccination (yellow fever, malaria, influenza), but poor outcome in autoimmune and inflammatory disease (type 1 diabetes, anti-neutrophil cytoplasmic antibody-associated vasculitis, systemic lupus erythematosus, idiopathic pulmonary fibrosis and dengue haemorrhagic fever). Thus, T-cell exhaustion plays a central role in determining outcome in autoimmune disease and targeted manipulation of this process could lead to new therapeutic opportunities

The Evolution of Compact Binary Star Systems

We review the formation and evolution of compact binary stars consisting of white dwarfs (WDs), neutron stars (NSs), and black holes (BHs). Binary NSs and BHs are thought to be the primary astrophysical sources of gravitational waves (GWs) within the frequency band of ground-based detectors, while compact binaries of WDs are important sources of GWs at lower frequencies to be covered by space interferometers (LISA). Major uncertainties in the current understanding of properties of NSs and BHs most relevant to the GW studies are discussed, including the treatment of the natal kicks which compact stellar remnants acquire during the core collapse of massive stars and the common envelope phase of binary evolution. We discuss the coalescence rates of binary NSs and BHs and prospects for their detections, the formation and evolution of binary WDs and their observational manifestations. Special attention is given to AM CVn-stars -- compact binaries in which the Roche lobe is filled by another WD or a low-mass partially degenerate helium-star, as these stars are thought to be the best LISA verification binary GW sources.Comment: 105 pages, 18 figure

Survival dimensionality reduction (SDR): development and clinical application of an innovative approach to detect epistasis in presence of right-censored data

Contains fulltext : 89126.pdf (publisher's version ) (Open Access)BACKGROUND: Epistasis is recognized as a fundamental part of the genetic architecture of individuals. Several computational approaches have been developed to model gene-gene interactions in case-control studies, however, none of them is suitable for time-dependent analysis. Herein we introduce the Survival Dimensionality Reduction (SDR) algorithm, a non-parametric method specifically designed to detect epistasis in lifetime datasets. RESULTS: The algorithm requires neither specification about the underlying survival distribution nor about the underlying interaction model and proved satisfactorily powerful to detect a set of causative genes in synthetic epistatic lifetime datasets with a limited number of samples and high degree of right-censorship (up to 70%). The SDR method was then applied to a series of 386 Dutch patients with active rheumatoid arthritis that were treated with anti-TNF biological agents. Among a set of 39 candidate genes, none of which showed a detectable marginal effect on anti-TNF responses, the SDR algorithm did find that the rs1801274 SNP in the Fc gamma RIIa gene and the rs10954213 SNP in the IRF5 gene non-linearly interact to predict clinical remission after anti-TNF biologicals. CONCLUSIONS: Simulation studies and application in a real-world setting support the capability of the SDR algorithm to model epistatic interactions in candidate-genes studies in presence of right-censored data. Availability: http://sourceforge.net/projects/sdrproject/

CD6 and Syntaxin Binding Protein 6 Variants and Response to Tumor Necrosis Factor Alpha Inhibitors in Danish Patients with Rheumatoid Arthritis

<div><h3>Background</h3><p>TNFα inhibitor therapy has greatly improved the treatment of patients with rheumatoid arthritis, however at least 30% do not respond. We aimed to investigate insertions and deletions (INDELS) associated with response to TNFα inhibitors in patients with rheumatoid arthritis (RA).</p> <h3>Methodology and Principal Findings</h3><p>In the DANBIO Registry we identified 237 TNFα inhibitor naïve patients with RA (81% women; median age 56 years; disease duration 6 years) who initiated treatment with infliximab (n = 160), adalimumab (n = 56) or etanercept (n = 21) between 1999 and 2008 according to national treatment guidelines. Clinical response was assessed at week 26 using EULAR response criteria. Based on literature, we selected 213 INDELS potentially related to RA and treatment response using the GeneVa® (Compugen) <em>in silico</em> database of 350,000 genetic variations in the human genome. Genomic segments were amplified by polymerase chain reaction (PCR), and genotyped by Sanger sequencing or fragment analysis. We tested the association between genotypes and EULAR good response versus no response, and EULAR good response versus moderate/no response using Fisher’s exact test. At baseline the median DAS28 was 5.1. At week 26, 68 (29%) patients were EULAR good responders, while 81 (34%) and 88 (37%) patients were moderate and non-responders, respectively. A 19 base pair insertion within the CD6 gene was associated with EULAR good response vs. no response (OR = 4.43, 95% CI: 1.99–10.09, p = 7.211×10⁻⁵) and with EULAR good response vs. moderate/no response (OR = 4.54, 95% CI: 2.29–8.99, p = 3.336×10⁻⁶). A microsatellite within the syntaxin binding protein 6 (STXBP6) was associated with EULAR good response vs. no response (OR = 4.01, 95% CI: 1.92–8.49, p = 5.067×10⁻⁵).</p> <h3>Conclusion</h3><p>Genetic variations within CD6 and STXBP6 may influence response to TNFα inhibitors in patients with RA.</p> </div

Multi-model seascape genomics identifies distinct environmental drivers of selection among sympatric marine species

Background As global change and anthropogenic pressures continue to increase, conservation and management increasingly needs to consider species’ potential to adapt to novel environmental conditions. Therefore, it is imperative to characterise the main selective forces acting on ecosystems, and how these may influence the evolutionary potential of populations and species. Using a multi-model seascape genomics approach, we compare putative environmental drivers of selection in three sympatric southern African marine invertebrates with contrasting ecology and life histories: Cape urchin (Parechinus angulosus), Common shore crab (Cyclograpsus punctatus), and Granular limpet (Scutellastra granularis). Results Using pooled (Pool-seq), restriction-site associated DNA sequencing (RAD-seq), and seven outlier detection methods, we characterise genomic variation between populations along a strong biogeographical gradient. Of the three species, only S. granularis showed significant isolation-by-distance, and isolation-by-environment driven by sea surface temperatures (SST). In contrast, sea surface salinity (SSS) and range in air temperature correlated more strongly with genomic variation in C. punctatus and P. angulosus. Differences were also found in genomic structuring between the three species, with outlier loci contributing to two clusters in the East and West Coasts for S. granularis and P. angulosus, but not for C. punctatus. Conclusion The findings illustrate distinct evolutionary potential across species, suggesting that species-specific habitat requirements and responses to environmental stresses may be better predictors of evolutionary patterns than the strong environmental gradients within the region. We also found large discrepancies between outlier detection methodologies, and thus offer a novel multi-model approach to identifying the principal environmental selection forces acting on species. Overall, this work highlights how adding a comparative approach to seascape genomics (both with multiple models and species) can elucidate the intricate evolutionary responses of ecosystems to global change

Transcriptome Analysis Describing New Immunity and Defense Genes in Peripheral Blood Mononuclear Cells of Rheumatoid Arthritis Patients

Background: Large-scale gene expression profiling of peripheral blood mononuclear cells from Rheumatoid Arthritis (RA) patients could provide a molecular description that reflects the contribution of diverse cellular responses associated with this disease. The aim of our study was to identify peripheral blood gene expression profiles for RA patients, using Illumina technology, to gain insights into RA molecular mechanisms. Methodology/Principal Findings: The Illumina Human-6v2 Expression BeadChips were used for a complete genome-wide transcript profiling of peripheral blood mononuclear cells (PBMCs) from 18 RA patients and 15 controls. Differential analysis per gene was performed with one-way analysis of variance (ANOVA) and P values were adjusted to control the False Discovery Rate (FDR < 5%). Genes differentially expressed at significant level between patients and controls were analyzed using Gene Ontology (GO) in the PANTHER database to identify biological processes. A differentially expression of 339 Reference Sequence genes (238 down-regulated and 101 up-regulated) between the two groups was observed. We identified a remarkably elevated expression of a spectrum of genes involved in Immunity and Defense in PBMCs of RA patients compared to controls. This result is confirmed by GO analysis, suggesting that these genes could be activated systemically in RA. No significant down-regulated ontology groups were found. Microarray data were validated by real time PCR in a set of nine genes showing a high degree of correlation. Conclusions/Significance: Our study highlighted several new genes that could contribute in the identification of innovative clinical biomarkers for diagnostic procedures and therapeutic interventions