Alexandra Zhernakova

Genetics of rheumatoid arthritis contributes to biology and drug discovery

A major challenge in human genetics is to devise a systematic strategy to integrate disease-associated variants with diverse genomic and biological data sets to provide insight into disease pathogenesis and guide drug discovery for complex traits such as rheumatoid arthritis (RA). Here we performed a genome-wide association study 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. We discovered 42 novel RA risk loci at a genome-wide level of significance, bringing the total to 101 (refs 2, 3, 4). We devised an in silico pipeline using established bioinformatics methods based on functional annotation, cis-acting expression quantitative trait loci and pathway analyses—as well as novel methods based on genetic overlap with human primary immunodeficiency, haematological cancer somatic mutations and knockout 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.

Dense genotyping identifies and localizes multiple common and rare variant association signals in celiac disease.

Using variants from the 1000 Genomes Project pilot European CEU dataset and data from additional resequencing studies, we densely genotyped 183 non-HLA risk loci previously associated with immune-mediated diseases in 12,041 individuals with celiac disease (cases) and 12,228 controls. We identified 13 new celiac disease risk loci reaching genome-wide significance, bringing the number of known loci (including the HLA locus) to 40. We found multiple independent association signals at over one-third of these loci, a finding that is attributable to a combination of common, low-frequency and rare genetic variants. Compared to previously available data such as those from HapMap3, our dense genotyping in a large sample collection provided a higher resolution of the pattern of linkage disequilibrium and suggested localization of many signals to finer scale regions. In particular, 29 of the 54 fine-mapped signals seemed to be localized to single genes and, in some instances, to gene regulatory elements. Altogether, we define the complex genetic architecture of the risk regions of and refine the risk signals for celiac disease, providing the next step toward uncovering the causal mechanisms of the disease.

Meta-Analysis of Genome-Wide Association Studies in Celiac Disease and Rheumatoid Arthritis Identifies Fourteen Non-HLA Shared Loci

Epidemiology and candidate gene studies indicate a shared genetic basis for celiac disease (CD) and rheumatoid arthritis (RA), but the extent of this sharing has not been systematically explored. Previous studies demonstrate that 6 of the established non-HLA CD and RA risk loci (out of 26 loci for each disease) are shared between both diseases. We hypothesized that there are additional shared risk alleles and that combining genome-wide association study (GWAS) data from each disease would increase power to identify these shared risk alleles. We performed a meta-analysis of two published GWAS on CD (4,533 cases and 10,750 controls) and RA (5,539 cases and 17,231 controls). After genotyping the top associated SNPs in 2,169 CD cases and 2,255 controls, and 2,845 RA cases and 4,944 controls, 8 additional SNPs demonstrated P<5×10−8 in a combined analysis of all 50,266 samples, including four SNPs that have not been previously confirmed in either disease: rs10892279 near the DDX6 gene (Pcombinedu200a=u200a1.2×10−12), rs864537 near CD247 (Pcombinedu200a=u200a2.2×10−11), rs2298428 near UBE2L3 (Pcombinedu200a=u200a2.5×10−10), and rs11203203 near UBASH3A (Pcombinedu200a=u200a1.1×10−8). We also confirmed that 4 gene loci previously established in either CD or RA are associated with the other autoimmune disease at combined P<5×10−8 (SH2B3, 8q24, STAT4, and TRAF1-C5). From the 14 shared gene loci, 7 SNPs showed a genome-wide significant effect on expression of one or more transcripts in the linkage disequilibrium (LD) block around the SNP. These associations implicate antigen presentation and T-cell activation as a shared mechanism of disease pathogenesis and underscore the utility of cross-disease meta-analysis for identification of genetic risk factors with pleiotropic effects between two clinically distinct diseases.

Trans-eQTLs Reveal That Independent Genetic Variants Associated with a Complex Phenotype Converge on Intermediate Genes, with a Major Role for the HLA

For many complex traits, genetic variants have been found associated. However, it is still mostly unclear through which downstream mechanism these variants cause these phenotypes. Knowledge of these intermediate steps is crucial to understand pathogenesis, while also providing leads for potential pharmacological intervention. Here we relied upon natural human genetic variation to identify effects of these variants on trans-gene expression (expression quantitative trait locus mapping, eQTL) in whole peripheral blood from 1,469 unrelated individuals. We looked at 1,167 published trait- or disease-associated SNPs and observed trans-eQTL effects on 113 different genes, of which we replicated 46 in monocytes of 1,490 different individuals and 18 in a smaller dataset that comprised subcutaneous adipose, visceral adipose, liver tissue, and muscle tissue. HLA single-nucleotide polymorphisms (SNPs) were 10-fold enriched for trans-eQTLs: 48% of the trans-acting SNPs map within the HLA, including ulcerative colitis susceptibility variants that affect plausible candidate genes AOAH and TRBV18 in trans. We identified 18 pairs of unlinked SNPs associated with the same phenotype and affecting expression of the same trans-gene (21 times more than expected, P<10−16). This was particularly pronounced for mean platelet volume (MPV): Two independent SNPs significantly affect the well-known blood coagulation genes GP9 and F13A1 but also C19orf33, SAMD14, VCL, and GNG11. Several of these SNPs have a substantially higher effect on the downstream trans-genes than on the eventual phenotypes, supporting the concept that the effects of these SNPs on expression seems to be much less multifactorial. Therefore, these trans-eQTLs could well represent some of the intermediate genes that connect genetic variants with their eventual complex phenotypic outcomes.

Differential association of the PTPN22 coding variant with autoimmune diseases in a Dutch population.

Protein tyrosine phosphatase PTPN22 is involved in the negative regulation of T-cell responsiveness. Recently, the association of a coding variant of the PTPN22 gene-R620W(1858C>T) with a number of autoimmune diseases has been described. Therefore, we tested the association of PTPN22 1858*T allele in Dutch early onset type 1 diabetes (T1D) and rheumatoid arthritis (RA) patients, as well as celiac disease (CD) patients, for which no previous study of PTPN22 has been reported. The PTPN22 variant was strongly associated with T1D in cases vs controls (P=2 × 10−7, OR=2.3, 95% CI=1.7–3.1) as well as in a transmission disequilibrium test in nuclear trios (P=9 × 10−9, OR=3.3, CI=2.1–5.0), RA (case/control: P=0.003, OR=1.8 CI =1.2–2.6), but not CD, in spite of a trend of increased homozygosity (P=0.05) and early age at onset (P=0.01). PTPN22 is not generally associated with T-cell mediated autoimmune diseases, although it might play a role in the CD patients with early clinical manifestation.

Evolutionary and Functional Analysis of Celiac Risk Loci Reveals SH2B3 as a Protective Factor against Bacterial Infection

Celiac disease (CD) is an intolerance to dietary proteins of wheat, barley, and rye. CD may have substantial morbidity, yet it is quite common with a prevalence of 1%-2% in Western populations. It is not clear why the CD phenotype is so prevalent despite its negative effects on human health, especially because appropriate treatment in the form of a gluten-free diet has only been available since the 1950s, when dietary gluten was discovered to be the triggering factor. The high prevalence of CD might suggest that genes underlying this disease may have been favored by the process of natural selection. We assessed signatures of selection for ten confirmed CD-associated loci in several genome-wide data sets, comprising 8154 controls from four European populations and 195 individuals from a North African population, by studying haplotype lengths via the integrated haplotype score (iHS) method. Consistent signs of positive selection for CD-associated derived alleles were observed in three loci: IL12A, IL18RAP, and SH2B3. For the SH2B3 risk allele, we also show a difference in allele frequency distribution (Fst) between HapMap phase II populations. Functional investigation of the effect of the SH2B3 genotype in response to lipopolysaccharide and muramyl dipeptide revealed that carriers of the SH2B3 rs3184504*A risk allele showed stronger activation of the NOD2 recognition pathway. This suggests that SH2B3 plays a role in protection against bacteria infection, and it provides a possible explanation for the selective sweep on SH2B3, which occurred sometime between 1200 and 1700 years ago.

Fine Mapping Seronegative and Seropositive Rheumatoid Arthritis to Shared and Distinct HLA Alleles by Adjusting for the Effects of Heterogeneity

Despite progress in defining human leukocyte antigen (HLA) alleles for anti-citrullinated-protein-autoantibody-positive (ACPA(+)) rheumatoid arthritis (RA), identifying HLA alleles for ACPA-negative (ACPA(-)) RA has been challenging because of clinical heterogeneity within clinical cohorts. We imputed 8,961 classical HLA alleles, amino acids, and SNPs from Immunochip data in a discovery set of 2,406 ACPA(-) RA case and 13,930 control individuals. We developed a statistical approach to identify and adjust for clinical heterogeneity within ACPA(-) RA and observed independent associations for serine and leucine at position 11 in HLA-DRβ1 (p = 1.4 × 10(-13), odds ratio [OR] = 1.30) and for aspartate at position 9 in HLA-B (p = 2.7 × 10(-12), OR = 1.39) within the peptide binding grooves. These amino acid positions induced associations at HLA-DRB1(∗)03 (encoding serine at 11) and HLA-B(∗)08 (encoding aspartate at 9). We validated these findings in an independent set of 427 ACPA(-) case subjects, carefully phenotyped with a highly sensitive ACPA assay, and 1,691 control subjects (HLA-DRβ1 Ser11+Leu11: p = 5.8 × 10(-4), OR = 1.28; HLA-B Asp9: p = 2.6 × 10(-3), OR = 1.34). Although both amino acid sites drove risk of ACPA(+) and ACPA(-) disease, the effects of individual residues at HLA-DRβ1 position 11 were distinct (p < 2.9 × 10(-107)). We also identified an association with ACPA(+) RA at HLA-A position 77 (p = 2.7 × 10(-8), OR = 0.85) in 7,279 ACPA(+) RA case and 15,870 control subjects. These results contribute to mounting evidence that ACPA(+) and ACPA(-) RA are genetically distinct and potentially have separate autoantigens contributing to pathogenesis. We expect that our approach might have broad applications in analyzing clinical conditions with heterogeneity at both major histocompatibility complex (MHC) and non-MHC regions.

CTLA4 is differentially associated with autoimmune diseases in the Dutch population.

Cytotoxic T lymphocyte-associated antigen 4 (CTLA4) is an important negative regulator of T-cell response and its genetic association with type 1 diabetes (T1D) has recently been demonstrated. The frequent co-association of autoimmune diseases (AID) and the implication from multiple genome scans that the CTLA4 gene region is a general autoimmune region, led us to study the role of CTLA4 in independent cohorts of T1D, coeliac disease (CD) and rheumatoid arthritis (RA) patients. We present independent data that confirm the association of CTLA4 in Dutch patients with juvenile onset T1D and show differential association of CTLA4 with CD and RA. The CTLA4 gene polymorphisms were tested for association in 350 T1D, 310 CD, 520 RA patients and 900 controls. In addition, 218 families were tested by the transmission disequilibrium test (TDT). T1D patients showed the highest association with the MH30*G: −1147*C: +49*G: CT60*G: JO37_3*G (haplotype 2) in both a case/control cohort (P=0.002, OR=1.42) and by TDT (P=0.02, OR=1.43). In contrast, this haplotype showed no association in the RA and CD cohorts. However, we observed an increased frequency of the MH30*G: −1147*T: +49*A: CT60*G: JO37_3*A (haplotype 3) in the CD patients diagnosed at a young age (OR=1.6, P=0.026, Pc=0.052). Furthermore, when T1D and CD patients were stratified based on the HLA risk, the T1D susceptible CTLA4 haplotype 2 was over-represented in the high HLA-risk T1D and CD groups. In conclusion, we confirmed association between CTLA4 haplotype 2 and T1D in the Dutch population. Association with another CTLA4 haplotype (haplotype 3) was confirmed for CD, but only in those patients who had an early age of expression. No effect was found between RA and CTLA4. The association of the CTLA4 haplotype 2 with the high-risk HLA genotype in T1D and CD, which share DQ2 as the one of high-risk alleles, might provide a clue to understanding the common genetic background of AID.

Rare, Low-Frequency, and Common Variants in the Protein-Coding Sequence of Biological Candidate Genes from GWASs Contribute to Risk of Rheumatoid Arthritis.

The extent to which variants in the protein-coding sequence of genes contribute to risk of rheumatoid arthritis (RA) is unknown. In this study, we addressed this issue by deep exon sequencing and large-scale genotyping of 25 biological candidate genes located within RA risk loci discovered by genome-wide association studies (GWASs). First, we assessed the contribution of rare coding variants in the 25 genes to the risk of RA in a pooled sequencing study of 500 RA cases and 650 controls of European ancestry. We observed an accumulation of rare nonsynonymous variants exclusive to RA cases in IL2RA and IL2RB (burden test: p = 0.007 and p = 0.018, respectively). Next, we assessed the aggregate contribution of low-frequency and common coding variants to the risk of RA by dense genotyping of the 25 gene loci in 10,609 RA cases and 35,605 controls. We observed a strong enrichment of coding variants with a nominal signal of association with RA (p < 0.05) after adjusting for the best signal of association at the loci (p(enrichment) = 6.4 × 10(-4)). For one locus containing CD2, we found that a missense variant, rs699738 (c.798C>A [p.His266Gln]), and a noncoding variant, rs624988, reside on distinct haplotypes and independently contribute to the risk of RA (p = 4.6 × 10(-6)). Overall, our results indicate that variants (distributed across the allele-frequency spectrum) within the protein-coding portion of a subset of biological candidate genes identified by GWASs contribute to the risk of RA. Further, we have demonstrated that very large sample sizes will be required for comprehensively identifying the independent alleles contributing to the missing heritability of RA.

Genetic variants of RANTES are associated with serum RANTES level and protection for type 1 diabetes

RANTES (regulated on activation, normal T-cell expressed and secreted) is a T-helper type 1 (Th1) chemokine that promotes T-cell activation and proliferation. RANTES is genetically associated with asthma, sarcoidosis and multiple sclerosis. The concentration of RANTES is increased at inflammation sites in different autoimmune diseases. Type 1 diabetes (T1D) is a Th1-mediated disease with complex genetic predisposition. We tested RANTES as a candidate gene for association with T1D using three single-nucleotide polymorphism (SNP) variants (rs4251719, rs2306630 and rs2107538) to capture haplotype information. The minor alleles of all SNPs were transmitted less frequently to T1D offspring (transmission rates 37.3% (P=0.002), 38.7% (P=0.007) and 41.0% (P=0.01)) and were less frequently present in patients compared to controls (P=0.009, 0.03 and 0.04, respectively). A similar protective effect was observed for the haplotype carrying three minor alleles (transmission disequilibrium test (TDT): P=0.003; odds ratio (OR)=0.55; confidence interval (CI): 0.37–0.83; case/control: P=0.03; OR=0.74; CI: 0.55–0.98). Both patients and controls carrying the protective haplotype express significantly lower serum levels of RANTES compared to non-carriers. Subsequently, we tested a cohort of 310 celiac disease patients, but failed to detect association. RANTES SNPs are significantly associated with RANTES serum concentration and development of T1D. The rs4251719*A–rs2306630*A–rs2107538*A haplotype associated with low RANTES production confers protection from T1D. Our data imply that RANTES is associated with T1D both genetically and functionally, and contributes to diabetes-prone Th1 cytokine profile.