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.
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.
HLA-DRB1 shared epitope (SE) alleles are the strongest genetic determinants for autoantibody positive rheumatoid arthritis (RA). One of the key regulators in expression of HLA class II receptors is MHC class II transactivator (CIITA). A variant of the CIITA gene has been found to associate with inflammatory diseases.
We wanted to explore whether the risk variant rs3087456 in the CIITA gene interacts with the HLA-DRB1 SE alleles regarding the risk of developing RA. We tested this hypothesis in a case-control study with 11767 individuals from four European Caucasian populations (6649 RA cases and 5118 controls).
We found no significant additive interaction for risk alleles among Swedish Caucasians with RA (n = 3869, attributable proportion due to interaction (AP) = 0.2, 95%CI: −0.2–0.5) or when stratifying for anti-citrullinated protein antibodies (ACPA) presence (ACPA positive disease: n = 2945, AP = 0.3, 95%CI: −0.05–0.6, ACPA negative: n = 2268, AP = −0.2, 95%CI: −1.0–0.6). We further found no significant interaction between the main subgroups of SE alleles (DRB1*01, DRB1*04 or DRB1*10) and CIITA. Similar analysis of three independent RA cohorts from British, Dutch and Norwegian populations also indicated an absence of significant interaction between genetic variants in CIITA and SE alleles with regard to RA risk.
Our data suggest that risk from the CIITA locus is independent of the major risk for RA from HLA-DRB1 SE alleles, given that no significant interaction between rs3087456 and SE alleles was observed. Since a biological link between products of these genes is evident, the genetic contribution from CIITA and class II antigens in the autoimmune process may involve additional unidentified factors.
Earlier studies have demonstrated that smoking and genetic risk factors interact in providing an increased risk for Rheumatoid Arthritis (RA). Less is known on how smoking contributes to RA in the context of genetic variability, and what proportion of RA that may be caused by smoking.
To determine the association between amount of smoking and risk of RA in the context of different HLA-DRB1 shared epitope (SE) alleles, and to estimate proportions of RA cases attributed to smoking.
Design, Setting, and Participants
Data from the Swedish Epidemiological Investigation of Rheumatoid Arthritis (EIRA) case-control study encompassing 1204 cases and 871 controls were analysed.
Main Outcome Measure
Estimated odds ratio to develop RA and excess fraction of cases attributable to smoking according to amount of smoking and genotype.
Smoking was estimated to be responsible for 35 % of the ACPA+ cases. For each HLA-DRB1 SE genotype, smoking was dose-dependently associated with increased risk of ACPA+ RA (p-trend<0.001). In individuals carrying two copies of the HLA-DRB1 shared epitope, 55 % of ACPA-positive RA were attributable to smoking.
Smoking is a preventable risk factor for RA. The increased risk due to smoking is dependent on amount of smoking and genotype.
Rheumatoid Arthritis; Smoking; HLA-DRB1; Case-control studies; Epidemiology
Selective IgA deficiency (IgAD; serum IgA<0.07 g/l) is the most common form of human primary immune deficiency, affecting approximately 1∶600 individuals in populations of Northern European ancestry. The polygenic nature of IgAD is underscored by the recent identification of several new risk genes in a genome-wide association study. Among the characterized susceptibility loci, the association with specific HLA haplotypes represents the major genetic risk factor for IgAD. Despite the robust association, the nature and location of the causal variants in the HLA region remains unknown. To better characterize the association signal in this region, we performed a high-density SNP mapping of the HLA locus and imputed the genotypes of common HLA-B, -DRB1, and -DQB1 alleles in a combined sample of 772 IgAD patients and 1,976 matched controls from 3 independent European populations. We confirmed the complex nature of the association with the HLA locus, which is the result of multiple effects spanning the entire HLA region. The primary association signal mapped to the HLA-DQB1*02 allele in the HLA Class II region (combined P = 7.69×10−57; OR = 2.80) resulting from the combined independent effects of the HLA-B*0801-DRB1*0301-DQB1*02 and -DRB1*0701-DQB1*02 haplotypes, while additional secondary signals were associated with the DRB1*0102 (combined P = 5.86×10−17; OR = 4.28) and the DRB1*1501 (combined P = 2.24×10−35; OR = 0.13) alleles. Despite the strong population-specific frequencies of HLA alleles, we found a remarkable conservation of these effects regardless of the ethnic background, which supports the use of large multi-ethnic populations to characterize shared genetic association signals in the HLA region. We also provide evidence for the location of association signals within the specific extended haplotypes, which will guide future sequencing studies aimed at characterizing the precise functional variants contributing to disease pathogenesis.
The human leukocyte antigen (HLA) locus is robustly associated with many immune-mediated conditions. However, identification of the genetic variants contributing to the disease pathophysiology has been greatly hampered by the extensive chromosomal conservation within this genomic region. To better understand the association of the HLA locus in selective IgA deficiency (IgAD), we used an extensive genotyping database from a recent genome-wide association study (GWAS) to generate a high-density SNP map of this region in a combined sample of >2,700 individuals from 3 independent European populations. In addition, we took advantage of recent methodological advances to impute the more common HLA-B, -DRB1, and -DQB1 alleles in all subjects. We confirmed the strong disease-association of the HLA locus and identified several different signals located in specific conserved HLA haplotypes contributing independent risk or protection for IgAD. Further analysis of the chromosomal sequences associated with the associated HLA alleles allowed us to refine the mapping of the susceptibility variants. These findings represent the most comprehensive high-density SNP mapping of the HLA locus in IgAD to date and provide important new information as to the location of the genetic variants contributing to this common immune deficiency.
The R620W variant in protein tyrosine phosphatase non-receptor 22 (PTPN22) is associated with rheumatoid arthritis (RA). The PTPN22 gene has alternatively spliced transcripts and at least two of the splice forms have been confirmed to encode different PTPN22 (LYP) proteins, but detailed information regarding expression of these is lacking, especially with regard to autoimmune diseases.
We have investigated the mRNA expression of known PTPN22 splice forms with TaqMan real-time PCR in relation to ZNF592 as an endogenous reference in peripheral blood cells from three independent cohorts with RA patients (n = 139) and controls (n = 111) of Caucasian origin. Polymorphisms in the PTPN22 locus (25 SNPs) and phenotypic data (gender, disease activity, ACPA and RF status) were used for analysis. Additionally, we addressed possible effects of methotrexate treatment on PTPN22 expression.
We found consistent differences in the expression of the PTPN22 splice forms in unstimulated peripheral blood mononuclear cells between RA patients and normal controls. This difference was more pronounced when comparing the ratio of splice forms and was not affected by methotrexate treatment.
Our data show that RA patients and healthy controls have a shift in balance of expression of splice forms derived from the PTPN22 gene. This balance seems not to be caused by treatment and may be of importance during immune response due to great structural differences in the encoded PTPN22 proteins.
A report of the 60th annual meeting of the American Society of Human Genetics in Washington, DC, USA, 2-6 November 2010.
Glutathione S-transferase (GST) and heme oxygenase-1 (HMOX1) genes encode enzymes that detoxify carcinogens and protect against oxidative stress. We studied gene-smoking interactions on rheumatoid arthritis (RA) susceptibility.
549 matched Caucasian RA cases and controls were selected from the Nurses' Health Study. Genotyping by TaqMan and BioTrove identified GSTM1, GSTT1 homozygous deletions (null) and GSTP1 (rs1695), HMOX1 (rs2071746) alleles, respectively. We studied gene-smoking interactions on the risk of all RA and serologic RA phenotypes in separate logistic models adjusted for covariates. We assessed multiplicative interactions with product terms in the logistic models and additive interactions with the attributable proportion due to interaction (AP). For replication, we repeated significant analyses in an independent case-control sample from the Epidemiologic Investigation of RA.
For all RA risk, we observed: multiplicative (p=0.05) and additive (AP=0.53, p=0.0005) interactions between GSTT1-null and smoking and multiplicative interaction (p=0.05) between HMOX1 and smoking. For seropositive RA risk, we found multiplicative (p=0.01) and additive (AP=0.63, p<0.0001) interactions between GSTT1-null and smoking and additive interaction (AP=0.41, p=0.03) between HMOX1 and smoking. After correction for multiple comparisons, additive interactions for GSTT1-null and smoking remained significant. GSTM1-null and GSTP1 did not show significant interactions. No associations were seen with seronegative RA. In replication analyses, we observed significant multiplicative (p=0.04) and additive (AP=0.32, p=0.02) interactions between GSTT1-null and smoking for ACPA-positive RA risk.
We observed significant gene-environment interactions between GSTT1-null and heavy smoking on RA risk. Future studies are needed to assess the impact of these interactions on RA prediction.
To discover novel RA risk loci, we systematically examined 370 SNPs from 179 independent loci with p<0.001 in a published meta-analysis of RA GWAS of 3,393 cases and 12,462 controls1. We used GRAIL2, a computational method that applies statistical text mining to PubMed abstracts, to score these 179 loci for functional relationships to genes in 16 established RA disease loci1,3-11. We identified 22 loci with a significant degree of functional connectivity. We genotyped 22 representative SNPs in an independent set of 7,957 cases and 11,958 matched controls. Three validate convincingly: CD2/CD58 (rs11586238, p=1×10−6 replication, p=1×10−9 overall), and CD28 (rs1980422, p=5×10−6 replication, p=1×10−9 overall), PRDM1 (rs548234, p=1×10−5 replication, p=2×10−8 overall). An additional four replicate (p<0.0023): TAGAP (rs394581, p=0.0002 replication, p=4×10−7 overall), PTPRC (rs10919563, p=0.0003 replication, p=7×10−7 overall), TRAF6/RAG1 (rs540386, p=0.0008 replication, p=4×10−6 overall), and FCGR2A (rs12746613, p=0.0022 replication, p=2×10−5 overall). Many of these loci are also associated to other immunologic diseases.
It has repeatedly been suggested that the development of complex diseases can be elucidated by gene–gene interactions. Recently, we found that HTR2A, a member of the serotonin receptor family, is associated with rheumatoid arthritis (RA). This study was aimed to investigate the possibility of a gene–gene interaction between HTR2A and the major genetic risk factor for RA, HLA-DRB1 shared epitope (SE) alleles. We studied 4095 RA cases and 3223 controls from three different populations – from Sweden, the United States and the Netherlands – to test for interaction between the protective HTR2A haplotype and HLA-DRB1 SE alleles. Further, we analyzed mRNA and/or protein expression of HTR2A and HLA-DR in biopsy samples and in synovial fibroblasts from RA patients. The interaction was defined as departure from additivity of effects using attributable proportion due to interaction. First, we could demonstrate and further replicate an interaction between a protective haplotype in HTR2A and HLA-DRB1 SE alleles regarding risk of developing autoantibody-positive RA. Second, we could show that both genes are constitutively expressed in fibroblasts from synovial tissue of RA patients, and, by double immunofluorescence staining, we demonstrated that these two proteins are colocalized in these cells. In conclusion, our data demonstrate a statistical interaction between HTR2A and HLA-DRB1 SE alleles and colocalization of the product of these two genes in inflamed synovial tissue, which suggest a possible biological relationship between these two proteins. This finding may lead to the development of treatment based on enhancing the protective features of 5-HT2A in individuals with a certain HLA genotype.
serotonin receptor; MHC class II; RA; gene–gene interaction
To investigate the associations between HLA-DRB1 shared epitope (SE) alleles and rheumatoid arthritis in subsets of rheumatoid arthritis defined by autoantibodies in three Asian populations from Malaysia.
1,079 rheumatoid arthritis patients and 1,470 healthy controls were included in the study. Levels of antibodies to citrullinated proteins (ACPA) and rheumatoid factors were assessed and the PCR-SSO method was used for HLA-DRB1 genotyping.
The proportion of ACPA positivity among Malay, Chinese and Indian rheumatoid arthritis patients were 62.9%, 65.2% and 68.6%, respectively. An increased frequency of SE alleles was observed in ACPA-positive rheumatoid arthritis among the three Asian ethnic groups. HLA-DRB1*10 was highly associated with rheumatoid arthritis susceptibility in these Asian populations. HLA-DRB1*0405 was significantly associated with susceptibility to rheumatoid arthritis in Malays and Chinese, but not in Indians. HLA-DRB1*01 did not show any independent effect as a risk factor for rheumatoid arthritis in this study and HLA-DRB1*1202 was protective in Malays and Chinese. There was no association between SE alleles and ACPA- negative rheumatoid arthritis in any of the three Asian ethnic groups.
The HLA-DRB1 SE alleles increase the risk of ACPA-positive rheumatoid arthritis in all three Asian populations from Malaysia.
GWAS owe their popularity to the expectation that they will make a major impact on diagnosis, prognosis and management of disease by uncovering genetics underlying clinical phenotypes. The dominant paradigm in GWAS data analysis so far consists of extensive reliance on methods that emphasize contribution of individual SNPs to statistical association with phenotypes. Multivariate methods, however, can extract more information by considering associations of multiple SNPs simultaneously. Recent advances in other genomics domains pinpoint multivariate causal graph-based inference as a promising principled analysis framework for high-throughput data. Designed to discover biomarkers in the local causal pathway of the phenotype, these methods lead to accurate and highly parsimonious multivariate predictive models. In this paper, we investigate the applicability of causal graph-based method TIE* to analysis of GWAS data. To test the utility of TIE*, we focus on anti-CCP positive rheumatoid arthritis (RA) GWAS datasets, where there is a general consensus in the community about the major genetic determinants of the disease.
Application of TIE* to the North American Rheumatoid Arthritis Cohort (NARAC) GWAS data results in six SNPs, mostly from the MHC locus. Using these SNPs we develop two predictive models that can classify cases and disease-free controls with an accuracy of 0.81 area under the ROC curve, as verified in independent testing data from the same cohort. The predictive performance of these models generalizes reasonably well to Swedish subjects from the closely related but not identical Epidemiological Investigation of Rheumatoid Arthritis (EIRA) cohort with 0.71-0.78 area under the ROC curve. Moreover, the SNPs identified by the TIE* method render many other previously known SNP associations conditionally independent of the phenotype.
Our experiments demonstrate that application of TIE* captures maximum amount of genetic information about RA in the data and recapitulates the major consensus findings about the genetic factors of this disease. In addition, TIE* yields reproducible markers and signatures of RA. This suggests that principled multivariate causal and predictive framework for GWAS analysis empowers the community with a new tool for high-quality and more efficient discovery.
This article was reviewed by Prof. Anthony Almudevar, Dr. Eugene V. Koonin, and Prof. Marianthi Markatou.
The GEIRA (Gene-Environment and Gene–Gene Interaction Research Application) algorithm and subsequent program is dedicated to genome-wide gene-environment and gene–gene interaction analysis. It implements concepts of both additive and multiplicative interaction as well as calculations based on dominant, recessive and co-dominant genetic models, respectively. Estimates of interactions are incorporated in a single table to make the output easily read. The algorithm is coded in both SAS and R. GEIRA is freely available to non-commercial users at http://www.epinet.se. Additional information, including user’s manual and example datasets is available online at http://www.epinet.se.
GWAS; Interaction; Program