Anorexia nervosa (AN) is a complex and heritable eating disorder characterized by dangerously low body weight. Neither candidate gene studies nor an initial genome wide association study (GWAS) have yielded significant and replicated results. We performed a GWAS in 2,907 cases with AN from 14 countries (15 sites) and 14,860 ancestrally matched controls as part of the Genetic Consortium for AN (GCAN) and the Wellcome Trust Case Control Consortium 3 (WTCCC3). Individual association analyses were conducted in each stratum and meta-analyzed across all 15 discovery datasets. Seventy-six (72 independent) SNPs were taken forward for in silico (two datasets) or de novo (13 datasets) replication genotyping in 2,677 independent AN cases and 8,629 European ancestry controls along with 458 AN cases and 421 controls from Japan. The final global meta-analysis across discovery and replication datasets comprised 5,551 AN cases and 21,080 controls. AN subtype analyses (1,606 AN restricting; 1,445 AN binge-purge) were performed. No findings reached genome-wide significance. Two intronic variants were suggestively associated: rs9839776 (P=3.01×10-7) in SOX2OT and rs17030795 (P=5.84×10-6) in PPP3CA. Two additional signals were specific to Europeans: rs1523921 (P=5.76×10-6) between CUL3 and FAM124B and rs1886797 (P=8.05×10-6) near SPATA13. Comparing discovery to replication results, 76% of the effects were in the same direction, an observation highly unlikely to be due to chance (P=4×10-6), strongly suggesting that true findings exist but that our sample, the largest yet reported, was underpowered for their detection. The accrual of large genotyped AN case-control samples should be an immediate priority for the field.
anorexia nervosa; eating disorders; GWAS; genome-wide association study; body mass index; metabolic
To identify new genetic associations with juvenile and adult dermatomyositis (DM).
We performed a genome-wide association study (GWAS) of adult and juvenile DM patients of European ancestry (n = 1178) and controls (n = 4724). To assess genetic overlap with other autoimmune disorders, we examined whether 141 single nucleotide polymorphisms (SNPs) outside the major histocompatibility complex (MHC) locus, and previously associated with autoimmune diseases, predispose to DM.
Compared to controls, patients with DM had a strong signal in the MHC region consisting of GWAS-level significance (P < 5x10−8) at 80 genotyped SNPs. An analysis of 141 non-MHC SNPs previously associated with autoimmune diseases showed that three SNPs linked with three genes were associated with DM, with a false discovery rate (FDR) < 0.05. These genes were phospholipase C like 1 (PLCL1, rs6738825, FDR=0.00089), B lymphoid tyrosine kinase (BLK, rs2736340, FDR=0.00031), and chemokine (C-C motif) ligand 21 (CCL21, rs951005, FDR=0.0076). None of these genes was previously reported to be associated with DM.
Our findings confirm the MHC as the major genetic region associated with DM and indicate that DM shares non-MHC genetic features with other autoimmune diseases, suggesting the presence of additional novel risk loci. This first identification of autoimmune disease genetic predispositions shared with DM may lead to enhanced understanding of pathogenesis and novel diagnostic and therapeutic approaches.
dermatomyositis; adult; juvenile; shared autoimmunity genes
To identify novel genetic risk factors for rheumatoid arthritis (RA), we conducted a genome-wide association study (GWAS) meta-analysis of 5,539 autoantibody positive RA cases and 20,169 controls of European descent, followed by replication in an independent set of 6,768 RA cases and 8,806 controls. Of 34 SNPs selected for replication, 7 novel RA risk alleles were identified at genome-wide significance (P<5×10−8) in analysis of all 41,282 samples. The associated SNPs are near genes of known immune function, including IL6ST, SPRED2, RBPJ, CCR6, IRF5, and PXK. We also refined the risk alleles at two established RA risk loci (IL2RA and CCL21) and confirmed the association at AFF3. These new associations bring the total number of confirmed RA risk loci to 31 among individuals of European ancestry. An additional 11 SNPs replicated at P<0.05, many of which are validated autoimmune risk alleles, suggesting that most represent bona fide RA risk alleles.
Treatment strategies blocking tumor necrosis factor (anti-TNF) have proven very successful in patients with rheumatoid arthritis (RA). However, a significant subset of patients does not respond for unknown reasons. Currently there are no means of identifying these patients prior to treatment. This study was aimed at identifying genetic factors predicting anti-TNF treatment outcome in patient with RA using a genome-wide association approach.
We conducted a multi-stage, genome-wide association study with a primary analysis of 2,557,253 single nucleotide polymorphisms (SNPs) in 882 RA patients receiving anti-TNF therapy included through the Dutch Rheumatoid Arthritis Monitoring (DREAM) registry and the database of Apotheekzorg. Linear regression analysis of changes in the Disease Activity Score in 28 joints after 14 weeks of treatment was performed using an additive model. Markers with a p<10−3 were selected for replication in 1,821 RA patients from three independent cohorts. Pathway analysis including all SNPs with a p-value < 10−3 was performed using Ingenuity.
Seven hundred seventy two markers demonstrated evidence of association with treatment outcome in the initial stage. Eight genetic loci showed improved p-value in the overall meta-analysis compared to the first stage, three of which (rs1568885, rs1813443 and rs4411591) showed directional consistency over all four studied cohorts. We were unable to replicate markers previously reported to be associated with anti-TNF outcome. Network analysis indicated strong involvement of biological processes underlying inflammatory response and cell morphology.
Using a multi-stage strategy, we have identified 8 genetic loci associated with response to anti-TNF treatment. Further studies are required to validate these findings in additional patient collections.
anti-TNF; gene polymorphism; pharmacogenetics; rheumatoid arthritis; genome-wide association study
A major challenge in human genetics is to devise a systematic strategy to integrate disease-associated variants with diverse genomic and biological datasets to provide insight into disease pathogenesis and guide drug discovery for complex traits such as rheumatoid arthritis (RA)1. Here, we performed a genome-wide association study (GWAS) meta-analysis in a total of >100,000 subjects of European and Asian ancestries (29,880 RA cases and 73,758 controls), by evaluating ~10 million single nucleotide polymorphisms (SNPs). We discovered 42 novel RA risk loci at a genome-wide level of significance, bringing the total to 1012–4. We devised an in-silico pipeline using established bioinformatics methods based on functional annotation5, cis-acting expression quantitative trait loci (cis-eQTL)6, and pathway analyses7–9 – as well as novel methods based on genetic overlap with human primary immunodeficiency (PID), hematological cancer somatic mutations and knock-out mouse phenotypes – to identify 98 biological candidate genes at these 101 risk loci. We demonstrate that these genes are the targets of approved therapies for RA, and further suggest that drugs approved for other indications may be repurposed for the treatment of RA. Together, this comprehensive genetic study sheds light on fundamental genes, pathways and cell types that contribute to RA pathogenesis, and provides empirical evidence that the genetics of RA can provide important information for drug discovery.
The interleukin 2 (IL-2) and interleukin 21 (IL-21) locus at chromosome 4q27 has been associated with several autoimmune diseases, and both genes are related to immune system functions. The aim of this study was to evaluate the role of the IL-2/IL-21 locus in systemic sclerosis (SSc).
Patients and methods
The case control study included 4493 SSc Caucasian patients and 5856 healthy controls from eight Caucasian populations (Spain, Germany, The Netherlands, USA, Italy, Sweden, UK and Norway). Four single nucleotide polymorphisms (rs2069762, rs6822844, rs6835457 and rs907715) were genotyped using TaqMan allelic discrimination assays.
We observed evidence of association of the rs6822844 and rs907715 variants with global SSc (pc=6.6E-4 and pc=7.2E-3, respectively). Similar statistically significant associations were observed for the limited cutaneous form of the disease. The conditional regression analysis suggested that the most likely genetic variation responsible for the association was the rs6822844 polymorphism. Consistently, the rs2069762A-rs6822844T-rs6835457G-rs907715T allelic combination showed evidence of association with SSc and limited cutaneous SSc subtype (pc=1.7E-03 and pc=8E-4, respectively).
These results suggested that the IL-2/IL-21 locus influences the genetic susceptibility to SSc. Moreover, this study provided further support for the IL-2/IL-21 locus as a common genetic factor in autoimmune diseases.
Integrating genetic data from families with highly penetrant forms of disease together with genetic data from outbred populations represents a promising strategy to uncover the complete frequency spectrum of risk alleles for complex traits such as rheumatoid arthritis (RA). Here, we demonstrate that rare, low-frequency and common alleles at one gene locus, phospholipase B1 (PLB1), might contribute to risk of RA in a 4-generation consanguineous pedigree (Middle Eastern ancestry) and also in unrelated individuals from the general population (European ancestry). Through identity-by-descent (IBD) mapping and whole-exome sequencing, we identified a non-synonymous c.2263G>C (p.G755R) mutation at the PLB1 gene on 2q23, which significantly co-segregated with RA in family members with a dominant mode of inheritance (P = 0.009). We further evaluated PLB1 variants and risk of RA using a GWAS meta-analysis of 8,875 RA cases and 29,367 controls of European ancestry. We identified significant contributions of two independent non-coding variants near PLB1 with risk of RA (rs116018341 [MAF = 0.042] and rs116541814 [MAF = 0.021], combined P = 3.2×10−6). Finally, we performed deep exon sequencing of PLB1 in 1,088 RA cases and 1,088 controls (European ancestry), and identified suggestive dispersion of rare protein-coding variant frequencies between cases and controls (P = 0.049 for C-alpha test and P = 0.055 for SKAT). Together, these data suggest that PLB1 is a candidate risk gene for RA. Future studies to characterize the full spectrum of genetic risk in the PLB1 genetic locus are warranted.
Anorexia nervosa (AN) is a complex and heritable eating disorder characterized by dangerously low body weight. Neither candidate gene studies nor an initial genome wide association study (GWAS) have yielded significant and replicated results. We performed a GWAS in 2,907 cases with AN from 14 countries (15 sites) and 14,860 ancestrally matched controls as part of the Genetic Consortium for AN (GCAN) and the Wellcome Trust Case Control Consortium 3 (WTCCC3). Individual association analyses were conducted in each stratum and meta-analyzed across all 15 discovery datasets. Seventy-six (72 independent) SNPs were taken forward for in silico (two datasets) or de novo (13 datasets) replication genotyping in 2,677 independent AN cases and 8,629 European ancestry controls along with 458 AN cases and 421 controls from Japan. The final global meta-analysis across discovery and replication datasets comprised 5,551 AN cases and 21,080 controls. AN subtype analyses (1,606 AN restricting; 1,445 AN binge-purge) were performed. No findings reached genome-wide significance. Two intronic variants were suggestively associated: rs9839776 (P=3.01×10−7) in SOX2OT and rs17030795 (P=5.84×10−6) in PPP3CA. Two additional signals were specific to Europeans: rs1523921 (P=5.76×10−6) between CUL3 and FAM124B and rs1886797 (P=8.05×10−6) near SPATA13. Comparing discovery to replication results, 76% of the effects were in the same direction, an observation highly unlikely to be due to chance (P= 4×10−6), strongly suggesting that true findings exist but that our sample, the largest yet reported, was underpowered for their detection. The accrual of large genotyped AN case-control samples should be an immediate priority for the field.
anorexia nervosa; eating disorders; GWAS; genome-wide association study; body mass index; metabolic
genetic association study; disease genetics; immunogenetics; liver
Naturally acquired microchimerism may arise in the mother and her child during pregnancy when bidirectional trafficking of cells occurs through the placental barrier. The occurrence of maternal microchimerism (maternal cells in the offspring) has been associated with several autoimmune diseases, especially in children. Systemic Lupus erythematosus (SLE) is an autoimmune disorder with a resemblance to graft-versus-host disease. The aim of this study was to investigate the association between maternal microchimerism in the blood and SLE.
Thirty-two patients with SLE, 17 healthy brothers of the patients, and an additional 12 unrelated healthy men were the subjects in this study. A single-nucleotide polymorphism unique to each mother was identified, and maternal microchimerism in the study group and in the control group was detected using a quantitative real-time polymerase chain reaction technique. No differences in the frequency or the concentration of maternal cells were apparent in the blood of patients with SLE or in that of the controls. Two patients and one control tested positive for maternal microchimerism, but the positive subjects were all negative at a follow-up 16 years later. The sensitivity of the method was estimated to 1/10.000.
These results show no association between SLE and maternal microchimerism. The frequency of maternal microchimerism in the blood of adults overall may be lower than earlier reported.
Recent genome-wide association studies (GWASs) conducted in Asian populations have identified novel risk loci for systemic lupus erythematosus (SLE). Here, we genotyped 10 single-nucleotide polymorphisms (SNPs) in eight such loci and investigated their disease associations in three independent Caucasian SLE case–control cohorts recruited from Sweden, Finland and the United States. The disease associations of the SNPs in ETS1, IKZF1, LRRC18-WDFY4, RASGRP3, SLC15A4, TNIP1 and 16p11.2 were replicated, whereas no solid evidence of association was observed for the 7q11.23 locus in the Caucasian cohorts. SLC15A4 was significantly associated with renal involvement in SLE. The association of TNIP1 was more pronounced in SLE patients with renal and immunological disorder, which is corroborated by two previous studies in Asian cohorts. The effects of all the associated SNPs, either conferring risk for or being protective against SLE, were in the same direction in Caucasians and Asians. The magnitudes of the allelic effects for most of the SNPs were also comparable across different ethnic groups. On the contrary, remarkable differences in allele frequencies between Caucasian and Asian populations were observed for all associated SNPs. In conclusion, most of the novel SLE risk loci identified by GWASs in Asian populations were also associated with SLE in Caucasian populations. We observed both similarities and differences with respect to the effect sizes and risk allele frequencies across ethnicities.
systemic lupus erythematosus; genetic-association study; Asian; Caucasian
Epigenetic mechanisms integrate genetic and environmental causes of disease. Comprehensive genome-wide analyses of epigenetic modifications have not demonstrated robust association with common diseases. Using Illumina HumanMethylation450 arrays on 354 ACPA positive rheumatoid arthritis (RA) cases and 337 controls, we identified two clusters within the MHC region whose differential methylation potentially mediates genetic risk for RA. To reduce confounding hampering previous epigenome-wide studies, we corrected for cellular heterogeneity by estimating and adjusting for cell-type proportions and used mediation analysis to filter out associations likely consequential to disease. Four CpGs also showed association between genotype and variance of methylation in addition to mean. The associations for both clusters replicated at least one CpG (p<0.01), with the rest showing suggestive association, in monocytes in an independent 12 cases and 12 controls. Thus, DNA methylation is a potential mediator of genetic risk.
Using the Immunochip custom single nucleotide polymorphism (SNP) array, designed for dense genotyping of 186 genome wide association study (GWAS) confirmed loci we analysed 11,475 rheumatoid arthritis cases of European ancestry and 15,870 controls for 129,464 markers. The data were combined in meta-analysis with GWAS data from additional independent cases (n=2,363) and controls (n=17,872). We identified fourteen novel loci; nine were associated with rheumatoid arthritis overall and 5 specifically in anti-citrillunated peptide antibody positive disease, bringing the number of confirmed European ancestry rheumatoid arthritis loci to 46. We refined the peak of association to a single gene for 19 loci, identified secondary independent effects at six loci and association to low frequency variants (minor allele frequency <0.05) at 4 loci. Bioinformatic analysis of the data generated strong hypotheses for the causal SNP at seven loci. This study illustrates the advantages of dense SNP mapping analysis to inform subsequent functional investigations.
Although genetic and non-genetic studies in mouse and human implicate the CD40 pathway in rheumatoid arthritis (RA), there are no approved drugs that inhibit CD40 signaling for clinical care in RA or any other disease. Here, we sought to understand the biological consequences of a CD40 risk variant in RA discovered by a previous genome-wide association study (GWAS) and to perform a high-throughput drug screen for modulators of CD40 signaling based on human genetic findings. First, we fine-map the CD40 risk locus in 7,222 seropositive RA patients and 15,870 controls, together with deep sequencing of CD40 coding exons in 500 RA cases and 650 controls, to identify a single SNP that explains the entire signal of association (rs4810485, P = 1.4×10−9). Second, we demonstrate that subjects homozygous for the RA risk allele have ∼33% more CD40 on the surface of primary human CD19+ B lymphocytes than subjects homozygous for the non-risk allele (P = 10−9), a finding corroborated by expression quantitative trait loci (eQTL) analysis in peripheral blood mononuclear cells from 1,469 healthy control individuals. Third, we use retroviral shRNA infection to perturb the amount of CD40 on the surface of a human B lymphocyte cell line (BL2) and observe a direct correlation between amount of CD40 protein and phosphorylation of RelA (p65), a subunit of the NF-κB transcription factor. Finally, we develop a high-throughput NF-κB luciferase reporter assay in BL2 cells activated with trimerized CD40 ligand (tCD40L) and conduct an HTS of 1,982 chemical compounds and FDA–approved drugs. After a series of counter-screens and testing in primary human CD19+ B cells, we identify 2 novel chemical inhibitors not previously implicated in inflammation or CD40-mediated NF-κB signaling. Our study demonstrates proof-of-concept that human genetics can be used to guide the development of phenotype-based, high-throughput small-molecule screens to identify potential novel therapies in complex traits such as RA.
A current challenge in human genetics is to follow-up “hits” from genome-wide association studies (GWAS) to guide drug discovery for complex traits. Previously, we identified a common variant in the CD40 locus as associated with risk of rheumatoid arthritis (RA). Here, we fine-map the CD40 signal of association through a combination of dense genotyping and exonic sequencing in large patient collections. Further, we demonstrate that the RA risk allele is a gain-of-function allele that increases the amount of CD40 on the surface of primary human B lymphocyte cells from healthy control individuals. Based on these observations, we develop a high-throughput assay to recapitulate the biology of the RA risk allele in a system suitable for a small molecule drug screen. After a series of primary screens and counter screens, we identify small molecules that inhibit CD40-mediated NF-kB signaling in human B cells. While this is only the first step towards a more comprehensive effort to identify CD40-specific inhibitors that may be used to treat RA, our study demonstrates a successful strategy to progress from a GWAS to a drug screen for complex traits such as RA.
Anti–tumor necrosis factor α (anti-TNF) therapy is a mainstay of treatment in rheumatoid arthritis (RA). The aim of the present study was to test established RA genetic risk factors to determine whether the same alleles also influence the response to anti-TNF therapy.
A total of 1,283 RA patients receiving etanercept, infliximab, or adalimumab therapy were studied from among an international collaborative consortium of 9 different RA cohorts. The primary end point compared RA patients with a good treatment response according to the European League Against Rheumatism (EULAR) response criteria (n = 505) with RA patients considered to be nonresponders (n = 316). The secondary end point was the change from baseline in the level of disease activity according to the Disease Activity Score in 28 joints (ΔDAS28). Clinical factors such as age, sex, and concomitant medications were tested as possible correlates of treatment response. Thirty-one single-nucleotide polymorphisms (SNPs) associated with the risk of RA were genotyped and tested for any association with treatment response, using univariate and multivariate logistic regression models.
Of the 31 RA-associated risk alleles, a SNP at the PTPRC (also known as CD45) gene locus (rs10919563) was associated with the primary end point, a EULAR good response versus no response (odds ratio [OR] 0.55, P = 0.0001 in the multivariate model). Similar results were obtained using the secondary end point, the ΔDAS28 (P = 0.0002). There was suggestive evidence of a stronger association in autoantibody-positive patients with RA (OR 0.55, 95% confidence interval [95% CI] 0.39–0.76) as compared with autoantibody-negative patients (OR 0.90, 95% CI 0.41–1.99).
Statistically significant associations were observed between the response to anti-TNF therapy and an RA risk allele at the PTPRC gene locus. Additional studies will be required to replicate this finding in additional patient collections.
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.
There are no genetic predictors of response to one of the most widely used classes of drugs in the treatment of rheumatoid arthritis—biological modifiers of the inflammatory cytokine tumor necrosis factor-alpha (or anti-TNF therapy). To identify genetic predictors, we performed the largest genome-wide association study (GWAS) to date as part of an international collaboration. In our study, which included 2,706 RA patients treated with one of three anti-TNF drugs, the most significant finding was restricted to RA patients treated with etanercept (P = 8×10−8), a drug that acts as a soluble receptor to bind circulating cytokine and prevents TNF from binding to its cell surface receptor. The associated variant influences expression of a nearby immune-related gene, CD84, whose expression is correlated with disease activity in RA patients. Together, our data support a model in which genomic factors related to CD84 expression serve as a predictor of disease activity and response to etanercept therapy among RA patients of European ancestry, but not anti-TNF therapies that act through different biological mechanisms or potentially in RA patients of other genetic ancestries.
STAT4 is an important transcription factor that contributes to the incidence and severity of different autoimmune diseases and is implicated in the antiviral immune responses in mice. In this study, we evaluated the role of STAT4 in human and murine herpes simplex virus 2 (HSV-2) infections. We show that STAT4 regulates antiviral gamma interferon (IFN-γ) responses and disease severity during chronic HSV-2 infections in humans and vaccine-induced IFN-γ-mediated protection against HSV-2 infection in mice. In a cohort of 228 HSV-2-infected individuals, representing both patients with recurrent disease and asymptomatic HSV-2 carriers, we found that genetic variations in the STAT4 gene were associated with asymptomatic HSV-2 infection, as well as with increased in vitro secretion of IFN-γ in response to the virus. Mice that lacked STAT4 had impaired HSV-2-specific IFN-γ production and delayed-type hypersensitivity responses following vaccination, which led to impaired viral clearance in the genital tract of vaccinated animals after a genital HSV-2 challenge. We conclude that STAT4 plays an important role in IFN-γ-mediated HSV-2-specific immunity, affecting the severity of genital HSV-2 infection.
The majority of our knowledge regarding disease-related mechanisms of uncontrolled citrullination and anti-citrullinated protein antibody development in rheumatoid arthritis (RA) was investigated in Caucasian populations. However, peptidylarginine deiminase (PADI) type 4 gene polymorphisms are associated with RA in East Asian populations and weak or no association was found in Caucasian populations. This study explores the association between the PADI4 polymorphisms and RA risk in a multiethnic population residing in South East Asia with the goal of elucidating generalizability of association in non-Caucasian populations.
A total of 320 SNPs from the PADI locus (including PADI1, PADI2, PADI3, PADI4 and PADI6 genes) were genotyped in 1,238 RA cases and 1,571 control subjects from the Malaysian Epidemiological Investigation of Rheumatoid Arthritis (MyEIRA) case-control study. Additionally, we conducted meta-analysis of our data together with the previously published studies of RA from East Asian populations.
The overall odds ratio (ORoverall) for the PADI4 (rs2240340) allelic model was 1.11 (95% confidence interval (CI) = 1.00 to 1.23, P = 0.04) and for the genotypic model was 1.20 (95% CI = 1.01 to 1.44, P = 0.04). Haplotype analysis for four selected PADI4 SNPs revealed a significant association of one with susceptibility (P = 0.001) and of another with a protective effect (P = 0.02). The RA susceptibility was further confirmed when combined meta-analysis was performed using these data together with data from five previously published studies from Asia comprising 5,192 RA cases and 4,317 control subjects (ORoverall = 1.23 (95% CI = 1.16 to 1.31, Pheterogeneity = 0.08) and 1.31 (95% CI = 1.20 to 1.44, Pheterogeneity = 0.32) in allele and genotype-based models, respectively). In addition, we also detected a novel association of PADI2 genetic variant rs1005753 with RA (ORoverall = 0.87 (95% CI = 0.77 to 0.99)).
Our study demonstrates an association between PADI4 and RA in the multiethnic population from South East Asia and suggests additional association with a PADI2 gene. The study thus provides further support for the notion that polymorphisms in genes for enzymes responsible for citrullination contribute to RA development in multiple populations of Asian descent.
To determine whether a polymorphism in the Fcγ receptor type IIIA (FCGR3A-F158V), influencing immunoglobulin G binding affinity, relates to the therapeutic efficacy of rituximab in rheumatoid arthritis (RA) patients.
Observational cohort study.
Three university hospital rheumatology units in Sweden.
Patients with established RA (n=177; 145 females and 32 males) who started rituximab (Mabthera) as part of routine care.
Primary outcome measures
Response to rituximab therapy in relation to FCGR3A genotype, including stratification for sex.
The frequency of responders differed significantly across FCGR3A genotypes (p=0.017 in a 3×2 contingency table). Heterozygous patients showed the highest response rate at 83%, as compared with patients carrying 158FF (68%) or 158VV (56%) (p=0.028 and 0.016, respectively). Among 158VV patients, response rates differed between male and female patients (p=0.036), but not among 158FF or 158VF patients (p=0.72 and 0.46, respectively).
Therapeutic efficacy of rituximab in RA patients is influenced by FCGR3A genotype, with the highest response rates found among heterozygous patients. This may suggest that different rituximab mechanisms of action in RA are optimally balanced in FCGR3A-158VF patients. Similar to the previously described associations with RA susceptibility and disease course, the impact of 158VV on rituximab response may be influenced by sex.
Rheumatology; Immunology; Therapeutics
The genetic association of the major histocompatibility complex (MHC) to rheumatoid arthritis risk has commonly been attributed to HLA-DRB1 alleles. Yet controversy persists about the causal variants in HLA-DRB1 and the presence of independent effects elsewhere in the MHC. Using existing genome-wide SNP data in 5,018 seropositive cases and 14,974 controls, we imputed and tested classical alleles and amino acid polymorphisms for HLA-A, B, C, DPA1, DPB1, DQA1, DQB1, and DRB1 along with 3,117 SNPs across the MHC. Conditional and haplotype analyses reveal that three amino acid positions (11, 71 and 74) in HLA-DRβ1, and single amino acid polymorphisms in HLA-B (position 9) and HLA-DPβ1 (position 9), all located in the peptide-binding grooves, almost completely explain the MHC association to disease risk. This study illustrates how imputation of functional variation from large reference panels can help fine-map association signals in the MHC.
Microsomal PGE synthase 1 (mPGES-1) is the terminal enzyme in the induced state of prostaglandin E2 (PGE2) synthesis and constitutes a therapeutic target for rheumatoid arthritis (RA) treatment. We examined the role of the prostaglandin E synthase (PTGES) gene polymorphism in susceptibility to and severity of RA and related variations in the gene to its function. The PTGES gene polymorphism was analyzed in 3081 RA patients and 1900 controls from two study populations: Swedish Epidemiological Investigation of Rheumatoid Arthritis (EIRA) and the Leiden Early Arthritis Clinic (Leiden EAC). Baseline disease activity score (DAS28) was employed as a disease severity measure. mPGES-1 expression was analyzed in synovial tissue from RA patients with known genotypes using immunohistochemistry. In the Swedish study population, among women a significant association with risk for RA was observed for PTGES single-nucleotide polymorphisms (SNPs) in univariate analysis and for the distinct haplotype. These results were substantiated by meta-analysis of data from EIRA and Leiden EAC studies with overall OR 1.31 (95% confidence interval 1.11–1.56). Several PTGES SNPs were associated with earlier onset of disease or with higher DAS28 in women with RA. Patients with the genotype associated with higher DAS28 exhibited significantly higher mPGES-1 expression in synovial tissue. Our data reveal a possible influence of PTGES polymorphism on the pathogenesis of RA and on disease severity through upregulation of mPGES-1 at the sites of inflammation. Genetically predisposed individuals may develop earlier and more active disease owing to this mechanism.
mPGES-1; gene polymorphism; rheumatoid arthritis; gender
Dendritic cell immunoreceptor (DCIR) has been implicated in development of autoimmune disorders in rodent and DCIR polymorphisms were associated with anti-citrullinated proteins antibodies (ACPA)-negative rheumatoid arthritis (RA) in Swedish Caucasians. This study was undertaken to further investigate whether DCIR polymorphisms are also risk factors for the development of RA in four Asian populations originated from China and Malaysia.
We genotyped two DCIR SNPs rs2377422 and rs10840759 in Han Chinese population (1,193 cases, 1,278 controls), to assess their association with RA. Subsequently, rs2377422 was further genotyped in three independent cohorts of Malaysian-Chinese subjects (MY_Chinese, 254 cases, 206 controls), Malay subjects (MY_ Malay, 515 cases, 986 controls), and Malaysian-Indian subjects (MY_Indian, 378 cases, 285 controls), to seek confirmation of association in various ethnic groups. Meta-analysis was preformed to evaluate the contribution of rs2377422 polymorphisms to the development of ACPA-negative RA in distinct ethnic groups. Finally, we carried out association analysis of rs2377422 polymorphisms with DCIR mRNA expression levels.
DCIR rs2377422 was found to be significantly associated with ACPA -negative RA in Han Chinese (OR 1.92, 95% CI 1.27–2.90, P = 0.0020). Meta-analysis confirms DCIR rs2377422 as a risk factor for ACPA-negative RA across distinct ethnic groups (ORoverall = 1.17, 95% CI 1.06–1.30, P = 0.003). The SNP rs2377422 polymorphism showed significant association with DCIR mRNA expression level, i.e. RA-risk CC genotype exhibit a significant increase in the expression of DCIR (P = 0.0023, Kruskal–Wallis).
Our data provide evidence for association between DCIR rs2377422 and RA in non-Caucasian populations and confirm the influence of DCIR polymorphisms on RA susceptibility, especially on ACPA-negative RA.
Rheumatoid arthritis (RA) is a multifactorial autoimmune disease in which genetic and environmental factors interact in the etiology. In this study, we investigated whether smoking and HLA-DRB1 shared-epitope (SE) alleles interact differently in the development of the two major subgroups of rheumatoid arthritis (RA), anti-citrullinated proteins antibody (ACPA)-positive and ACPA-negative disease, in a multiethnic population of Asian descent.
A case-control study comprising early diagnosed RA cases was carried out in Malaysia between 2005 and 2009. In total, 1,076 cases and 1,612 matched controls participated in the study. High-resolution HLA-DRB1 genotyping was performed for shared-epitope (SE) alleles. All participants answered a questionnaire on a broad range of issues, including smoking habits. The odds ratio (OR) of developing ACPA-positive and ACPA-negative disease was calculated for smoking and the presence of any SE alleles separately. Potential interaction between smoking history (defined as "ever" and "never" smoking) and HLA-DRB1 SE alleles also was calculated.
In our multiethnic study, both the SE alleles and smoking were associated with an increased risk of developing ACPA-positive RA (OR SE alleles, 4.7; 95% confidence interval (CI), 3.6 to 6.2; OR smoking, 4.1; 95% CI, 1.9 to 9.2). SE-positive smokers had an odds ratio of ACPA-positive RA of 25.6 (95% CI, 10.4 to 63.4), compared with SE-negative never-smokers. The interaction between smoking and SE alleles was significant (attributable proportion due to interaction (AP) was 0.7 (95% CI, 0.5 to 1.0)). The HLA-DRB1*04:05 SE allele, which is common in Asian populations, but not among Caucasians, was associated with an increased risk of ACPA-positive RA, and this allele also showed signs of interaction with smoking (AP, 0.4; 95% CI, -0.1 to 0.9). Neither smoking nor SE alleles nor their combination was associated with an increased risk of ACPA-negative RA.
The risk of developing ACPA-positive RA is associated with a strong gene-environment interaction between smoking and HLA-DRB1 SE alleles in a Malaysian multiethnic population of Asian descent. This interaction seems to apply also between smoking and the specific HLA-DRB1*04:05 SE allele, which is common in Asian populations but not in Caucasians.
Objective. Genetic heterogeneity and risk factor distribution was analyzed in two previously proposed asthma phenotypes. Method. A sample of 412 subjects was investigated at 7-8, 12-13, and 21-22 years of age with questionnaires, skin prick tests, and genetic analysis of IL-4 receptor (IL4R) single-nucleotide polymorphisms. The sample was subdivided in one group with no asthma, and two groups with asthma separated by age of onset of symptoms, namely, early onset asthma (EOA) and late onset asthma (LOA). Risk factors and IL4R markers were analyzed in respect to asthma phenotypes. Results. EOA and LOA groups were both associated with atopy and a maternal history of asthma. Female gender was more common in LOA, whereas childhood eczema, frequent colds in infancy, and a paternal history of asthma were more common in EOA. The AA genotype of rs2057768 and the GG genotype of rs1805010 were more common in LOA, whereas the GG genotype of rs2107356 was less common in EOA. Conclusion. Our data suggest that early and late onset asthma may be of different endotypes and genotypes.