GWASs and the age of human as the model organism for autoimmune genetic research

Humans should perhaps be the model organism for autoimmune disease studies, considering the success of GWAS studies uncovering risk alleles

Rare protection against type 1 diabetes

Next-generation DNA sequencing reveals rare alleles protective from type 1 diabetes

Genetic polymorphisms in PTPN22, PADI-4, and CTLA-4 and risk for rheumatoid arthritis in two longitudinal cohort studies: evidence of gene-environment interactions with heavy cigarette smoking

Introduction: PTPN22, PADI-4, and CTLA-4 have been associated with risk for rheumatoid arthritis (RA). We investigated whether polymorphisms in these genes were associated with RA in Caucasian women included in two large prospective cohorts, adjusting for confounding factors and testing for interactions with smoking. Methods: We studied RA risk associated with PTPN22 (rs2476601), PADI-4 (rs2240340), and CTLA-4 (rs3087243) in the Nurses' Health Study (NHS) and NHSII. Participants in NHS were aged 30 to 55 years at entry in 1976; those in NHSII were aged 25 to 42 years at entry in 1989. We confirmed incident RA cases through to 2002 in NHS and to 2003 in NHSII by questionnaire and medical record review. We excluded reports not confirmed as RA. In a nested case-control design involving participants for whom there were samples for genetic analyses (45% of NHS and 25% of NHSII), each incident RA case was matched to a participant without RA by year of birth, menopausal status, and postmenopausal hormone use. Genotyping was performed using Taqman single nucleotide polymorphism allelic discrimination on the ABI 7900 HT (Applied Biosystems, 850 Lincoln Centre Drive, Foster City, CA 94404 USA) with published primers. Human leukocyte antigen shared epitope (HLA-SE) genotyping was performed at high resolution. We employed conditional logistic regression analyses, adjusting for smoking and reproductive factors. We tested for additive and multiplicative interactions between each genotype and smoking. Results: A total of 437 incident RA cases were matched to healthy female control individuals. Mean (± standard deviation) age at RA diagnosis was 55 (± 10), 57% of RA cases were rheumatoid factor (RF) positive, and 31% had radiographic erosions at diagnosis. PTPN22 was associated with increased RA risk (pooled odds ratio in multivariable dominant model = 1.46, 95% confidence interval [CI] = 1.02 to 2.08). The risk was stronger for RF-positive than for RF-negative RA. A significant multiplicative interaction between PTPN22 and smoking for more than 10 pack-years was observed (P = 0.04). CTLA-4 and PADI-4 genotypes were not associated with RA risk in the pooled results (pooled odds ratios in multivariable dominant models: 1.27 [95% CI = 0.88 to 1.84] for CTLA-4 and 1.04 [95% CI = 0.77 to 1.40] for PADI-4). No gene-gene interaction was observed between PTPN22 and HLA-SE. Conclusion: After adjusting for smoking and reproductive factors, PTPN22 was associated with RA risk among Caucasian women in these cohorts. We found both additive and multiplicative interactions between PTPN22 and heavy cigarette smoking

Genetic Risk Score Predicting Risk of Rheumatoid Arthritis Phenotypes and Age of Symptom Onset

Cumulative genetic profiles can help identify individuals at high-risk for developing RA. We examined the impact of 39 validated genetic risk alleles on the risk of RA phenotypes characterized by serologic and erosive status.We evaluated single nucleotide polymorphisms at 31 validated RA risk loci and 8 Human Leukocyte Antigen alleles among 542 Caucasian RA cases and 551 Caucasian controls from Nurses' Health Study and Nurses' Health Study II. We created a weighted genetic risk score (GRS) and evaluated it as 7 ordinal groups using logistic regression (adjusting for age and smoking) to assess the relationship between GRS group and odds of developing seronegative (RF- and CCP-), seropositive (RF+ or CCP+), erosive, and seropositive, erosive RA phenotypes. In separate case only analyses, we assessed the relationships between GRS and age of symptom onset. In 542 RA cases, 317 (58%) were seropositive, 163 (30%) had erosions and 105 (19%) were seropositive with erosions. Comparing the highest GRS risk group to the median group, we found an OR of 1.2 (95% CI = 0.8-2.1) for seronegative RA, 3.0 (95% CI = 1.9-4.7) for seropositive RA, 3.2 (95% CI = 1.8-5.6) for erosive RA, and 7.6 (95% CI = 3.6-16.3) for seropositive, erosive RA. No significant relationship was seen between GRS and age of onset.Results suggest that seronegative and seropositive/erosive RA have different genetic architecture and support the importance of considering RA phenotypes in RA genetic studies

Allele-Specific Methylation Occurs at Genetic Variants Associated with Complex Disease

We hypothesize that the phenomenon of allele-specific methylation (ASM) may underlie the phenotypic effects of multiple variants identified by Genome-Wide Association studies (GWAS). We evaluate ASM in a human population and document its genome-wide patterns in an initial screen at up to 380,678 sites within the genome, or up to 5% of the total genomic CpGs. We show that while substantial inter-individual variation exists, 5% of assessed sites show evidence of ASM in at least six samples; the majority of these events (81%) are under genetic influence. Many of these cis-regulated ASM variants are also eQTLs in peripheral blood mononuclear cells and monocytes and/or in high linkage-disequilibrium with variants linked to complex disease. Finally, focusing on autoimmune phenotypes, we extend this initial screen to confirm the association of cis-regulated ASM with multiple complex disease-associated variants in an independent population using next-generation bisulfite sequencing. These four variants are implicated in complex phenotypes such as ulcerative colitis and AIDS progression disease (rs10491434), Celiac disease (rs2762051), Crohn's disease, IgA nephropathy and early-onset inflammatory bowel disease (rs713875) and height (rs6569648). Our results suggest cis-regulated ASM may provide a mechanistic link between the non-coding genetic changes and phenotypic variation observed in these diseases and further suggests a route to integrating DNA methylation status with GWAS results

Identification of the NF-κB activating protein-like locus as a risk locus for rheumatoid arthritis

Objective: To fine-map the NF-κB activating protein-like (NKAPL) locus identified in a prior genome-wide study as a possible rheumatoid arthritis (RA) risk locus and thereby delineate additional variants with stronger and/or independent disease association. Methods: Genotypes for 101 SNPs across the NKAPL locus on chromosome 6p22.1 were obtained on 1368 Canadian RA cases and 1471 controls. Single marker associations were examined using logistic regression and the most strongly associated NKAPL locus SNPs then typed in another Canadian and a US-based RA case/control cohort. Results: Fine-mapping analyses identified six NKAPL locus variants in a single haplotype block showing association with p≤5.6×10−8 in the combined Canadian cohort. Among these SNPs, rs35656932 in the zinc finger 193 gene and rs13208096 in the NKAPL gene remained significant after conditional logistic regression, contributed independently to risk for disease, and were replicated in the US cohort (Pcomb=4.24×10−10 and 2.44×10−9, respectively). These associations remained significant after conditioning on SNPs tagging the HLA-shared epitope (SE) DRB1*0401 allele and were significantly stronger in the HLA-SE negative versus positive subgroup, with a significant negative interaction apparent between HLA-DRB1 SE and NKAPL risk alleles. Conclusions: By illuminating additional NKAPL variants with highly significant effects on risk that are distinct from, but interactive with those arising from the HLA-DRB1 locus, our data conclusively identify NKAPL as an RA susceptibility locus