Association with HLA-DRß1 position 37 distinguishes juvenile dermatomyositis from adult-onset myositis.

Juvenile dermatomyositis (JDM) is a rare, severe autoimmune disease and the most common idiopathic inflammatory myopathy of children. JDM and adult-onset dermatomyositis (DM) have similar clinical, biological and serological features, although these features differ in prevalence between childhood-onset and adult-onset disease, suggesting that age of disease onset may influence pathogenesis. Therefore, a JDM-focused genetic analysis was performed using the largest collection of JDM samples to date. Caucasian JDM samples (n = 952) obtained via international collaboration were genotyped using the Illumina HumanCoreExome chip. Additional non-assayed human leukocyte antigen (HLA) loci and genome-wide single-nucleotide polymorphisms (SNPs) were imputed. HLA-DRB1*03:01 was confirmed as the classical HLA allele most strongly associated with JDM [odds ratio (OR) 1.66; 95% confidence interval (CI) 1.46, 1.89; P = 1.4 × 10-14], with an independent association at HLA-C*02:02 (OR = 1.74; 95% CI 1.42, 2.13, P = 7.13 × 10-8). Analyses of amino acid positions within HLA-DRB1 indicated that the strongest association was at position 37 (omnibus P = 3.3 × 10-19), with suggestive evidence this association was independent of position 74 (omnibus P = 5.1 × 10-5), the position most strongly associated with adult-onset DM. Conditional analyses also suggested that the association at position 37 of HLA-DRB1 was independent of some alleles of the Caucasian HLA 8.1 ancestral haplotype (AH8.1) such as HLA-DQB1*02:01 (OR = 1.62; 95% CI 1.36, 1.93; P = 8.70 × 10-8), but not HLA-DRB1*03:01 (OR = 1.49; 95% CR 1.24, 1.80; P = 2.24 × 10-5). No associations outside the HLA region were identified. Our findings confirm previous associations with AH8.1 and HLA-DRB1*03:01, HLA-C*02:02 and identify a novel association with amino acid position 37 within HLA-DRB1, which may distinguish JDM from adult DM.

Genome-wide association analysis identifies a susceptibility locus for sporadic vestibular schwannoma at 9p21.

Vestibular schwannomas are benign nerve sheath tumours that arise on the vestibulocochlear nerves. Vestibular schwannomas are known to occur in the context of tumour predisposition syndromes NF2-related and LZTR1-related schwannomatosis. However, the majority of vestibular schwannomas present sporadically without identification of germline pathogenic variants. To identify novel genetic associations with risk of vestibular schwannoma development, we conducted a genome-wide association study in a cohort of 911 sporadic vestibular schwannoma cases collated from the neurofibromatosis type 2 genetic testing service in the north-west of England, UK and 5500 control samples from the UK Biobank resource. One risk locus reached genome-wide significance in our association analysis (9p21.3, rs1556516, P = 1.47 × 10-13, odds ratio = 0.67, allele frequency = 0.52). 9p21.3 is a genome-wide association study association hotspot, and a number of genes are localized to this region, notably CDKN2B-AS1 and CDKN2A/B, also referred to as the INK4 locus. Dysregulation of gene products within the INK4 locus have been associated with multiple pathologies and the genes in this region have been observed to directly impact the expression of one another. Recurrent associations of the INK4 locus with components of well-described oncogenic pathways provides compelling evidence that the 9p21.3 region is truly associated with risk of vestibular schwannoma tumorigenesis.

The Mental Health of First Nations Children in Manitoba: A Population-Based Retrospective Cohort Study Using Linked Administrative Data: La santé mentale des enfants des Premières Nations au Manitoba : une étude de cohorte rétrospective dans la population, à l'aide de données administratives liées.

OBJECTIVE: First Nations children face a greater risk of experiencing mental disorders than other children from the general population because of family and societal factors, yet there is little research examining their mental health. This study compares diagnosed mental disorders and suicidal behaviours of First Nations children living on-reserve and off-reserve to all other children living in Manitoba. METHOD: The research team, which included First Nations and non-First Nations researchers, utilized population-based administrative data that linked de-identified individual-level records from the 2016 First Nations Research File to health and social information for children living in Manitoba. Adjusted rates and rate ratios of mental disorders and suicide behaviours were calculated using a generalized linear modelling approach to compare First Nations children (n = 40,574) and all other children (n = 197,109) and comparing First Nations children living on- and off-reserve. RESULTS: Compared with all other children, First Nations children had a higher prevalence of schizophrenia (adjusted rate ratio (aRR): 4.42, 95% confidence interval (CI), 3.36 to 5.82), attention-deficit hyperactivity disorder (ADHD; aRR: 1.21, 95% CI, 1.09 to 1.33), substance use disorders (aRR: 5.19; 95% CI, 4.25 to 6.33), hospitalizations for suicide attempts (aRR: 6.96; 95% CI, 4.36 to 11.13) and suicide deaths (aRR: 10.63; 95% CI, 7.08 to 15.95). The prevalence of ADHD and mood/anxiety disorders was significantly higher for First Nations children living off-reserve compared with on-reserve; in contrast, hospitalization rates for suicide attempts were twice as high on-reserve than off-reserve. When the comparison cohort was restricted to only other children in low-income areas, a higher prevalence of almost all disorders remained for First Nations children. CONCLUSION: Large disparities were found in mental health indicators between First Nations children and other children in Manitoba, demonstrating that considerable work is required to improve the mental well-being of First Nations children. Equitable access to culturally safe services is urgently needed and these services should be self-determined, planned, and implemented by First Nations people.

EULAR study group on 'MHC-I-opathy': identifying disease-overarching mechanisms across disciplines and borders.

The 'MHC-I (major histocompatibility complex class I)-opathy' concept describes a family of inflammatory conditions with overlapping clinical manifestations and a strong genetic link to the MHC-I antigen presentation pathway. Classical MHC-I-opathies such as spondyloarthritis, Behçet's disease, psoriasis and birdshot uveitis are widely recognised for their strong association with certain MHC-I alleles and gene variants of the antigen processing aminopeptidases ERAP1 and ERAP2 that implicates altered MHC-I peptide presentation to CD8+T cells in the pathogenesis. Progress in understanding the cause and treatment of these disorders is hampered by patient phenotypic heterogeneity and lack of systematic investigation of the MHC-I pathway.Here, we discuss new insights into the biology of MHC-I-opathies that strongly advocate for disease-overarching and integrated molecular and clinical investigation to decipher underlying disease mechanisms. Because this requires transformative multidisciplinary collaboration, we introduce the EULAR study group on MHC-I-opathies to unite clinical expertise in rheumatology, dermatology and ophthalmology, with fundamental and translational researchers from multiple disciplines such as immunology, genomics and proteomics, alongside patient partners. We prioritise standardisation of disease phenotypes and scientific nomenclature and propose interdisciplinary genetic and translational studies to exploit emerging therapeutic strategies to understand MHC-I-mediated disease mechanisms. These collaborative efforts are required to address outstanding questions in the etiopathogenesis of MHC-I-opathies towards improving patient treatment and prognostication.

Optimization of polygenic risk scores in BRCA1/2 pathogenic variant heterozygotes in epithelial ovarian cancer.

PURPOSE: A third of familial epithelial ovarian cancer (EOC) is explained by BRCA1/2 pathogenic variants. Polygenic risk scores (PRSs) for BRCA1/2 heterozygotes associated with EOC have been created, but impact of combination with clinical and hormonal risk factors is unclear. METHODS: We genotyped 300 cases and 355 controls and constructed modified PRSs based on those validated by Barnes et al. Model discrimination and EOC risk was assessed by area under the curve (AUC) values and difference between lowest and highest quintile odds ratios (ORs). We investigated model optimization using logistic regression to combine models with clinical and hormonal data. RESULTS: Unadjusted AUC values ranged from 0.526 to 0.551 with 2.2- to 2.3-fold increase in OR between lowest and highest quintiles (BRCA1 heterozygotes) and 0.574 to 0.585 AUC values with a 6.3- to 7.7-fold increase (BRCA2 heterozygotes). The optimized model (parity, age at menarche, menopause, and first full-term pregnancy) estimated AUC values of 0.872 to 0.876 and 21- to 23-fold increase in OR (BRCA1 heterozygotes) and AUC values of 0.857 to 0.867 and 40- to 41-fold increase (BRCA2 heterozygotes). CONCLUSION: The combination of PRS with age, family history, and hormonal factors significantly improved the EOC risk discrimination ability. However, the contribution of the PRS was small. Larger prospective studies are needed to assess if combined-PRS models could provide information to inform risk-reducing decisions.

Validation of lung cancer polygenic risk scores in a high-risk case-control cohort.

PURPOSE: Screening with low-dose computed tomography reduces lung cancer (LC) mortality. Risk prediction models used for screening selection do not include genetic variables. Here, we investigated the performance of previously published polygenic risk scores (PRSs) for LC, considering their potential to improve screening selection. METHODS: We validated 9 PRSs in a high-risk case-control cohort, comprising genotype data from 652 surgical patients with LC and 550 cancer-free, high-risk (PLCOM2012 score ≥ 1.51%) participants of the Manchester Lung Health Check, a community-based LC screening program (n = 550). Discrimination (area under the curve [AUC]) between cases and controls was assessed for each PRS independently and alongside clinical risk factors. RESULTS: Median age was 67 years, 53% were female, 46% were current smokers, and 76% were National Lung Screening Trial eligible. Median PLCOM2012 score among controls was 3.4%, 80% of cases were early stage. All PRSs significantly improved discrimination, AUC increased between +0.002 (P = .02) and +0.015 (P < .0001), compared with clinical risk factors alone. The best-performing PRS had an independent AUC of 0.59. Two novel loci, in the DAPK1 and MAGI2 genes, were significantly associated with LC risk. CONCLUSION: PRSs may improve LC risk prediction and screening selection. Further research, particularly examining clinical utility and cost-effectiveness, is required.

Association between psoriatic disease and lifestyle factors and comorbidities: cross-sectional analysis and Mendelian randomization.

OBJECTIVES: To examine associations between PsA and psoriasis vs lifestyle factors and comorbidities by triangulating observational and genetic evidence. METHODS: We analysed cross-sectional data from the UK Biobank (1836 PsA, 8995 psoriasis, 36 000 controls) to describe the association between psoriatic disease and lifestyle factors (including BMI and smoking) and 15 comorbidities [including diabetes and coronary artery disease (CAD)] using logistic models adjusted for age, sex and lifestyle factors. We applied bidirectional Mendelian randomization (MR) to genome-wide association data (3609 PsA and 7804 psoriasis cases, up to 1.2 million individuals for lifestyle factors and 757 601 for comorbidities) to examine causal direction, using the inverse-variance weighted method. RESULTS: BMI was cross-sectionally associated with risk of PsA (OR 1.31 per 5 kg/m2 increase; 95% CI 1.26, 1.37) and psoriasis (OR 1.23; 1.20, 1.26), with consistent MR estimates (PsA OR 1.38; 1.14, 1.67; psoriasis OR 1.36; 1.18, 1.58). In both designs, smoking was more strongly associated with psoriasis than PsA. PsA and psoriasis were cross-sectionally associated with diabetes (OR 1.35 and 1.39, respectively) and CAD (OR 1.56 and 1.38, respective). Genetically predicted glycated haemoglobin (surrogate for diabetes) increased PsA risk (OR 1.18 per 6.7 mmol/mol increase; 1.02, 1.36) but not psoriasis. Genetic liability to PsA (OR 1.05; 1.003, 1.09) and psoriasis (OR 1.03; 1.001, 1.06) were associated with increased risk of CAD. CONCLUSION: Observational and genetic evidence converge to suggest that BMI and glycaemic control are associated with increased psoriatic disease risk, while psoriatic disease is associated with increased CAD risk. Further research is needed to understand the mechanism of these associations.

No evidence that genetic predictors of susceptibility predict changes in core outcomes in JIA.

OBJECTIVES: The clinical progression of JIA is unpredictable. Knowing who will develop severe disease could facilitate rapid intensification of therapies. We use genetic variants conferring susceptibility to JIA to predict disease outcome measures. METHODS: A total of 713 JIA patients with genotype data and core outcome variables (COVs) at diagnosis (baseline) and 1 year follow-up were identified from the Childhood Arthritis Prospective Study (CAPS). A weighted genetic risk score (GRS) was generated, including all single nucleotide polymorphisms (SNPs) previously associated with JIA susceptibility (P-value < 5×10-08). We used multivariable linear regression to test the GRS for association with COVS (limited joint count, active joint count, physician global assessment, parent/patient general evaluation, childhood HAQ and ESR) at baseline and change in COVS from baseline to 1 year, adjusting for baseline COV and International League of Associations of Rheumatology (ILAR) category. The GRS was split into quintiles to identify high (quintile 5) and low (quintile 1) risk groups. RESULTS: Patients in the high-risk group for the GRS had a younger age at presentation (median low risk 7.79, median high risk 3.51). No association was observed between the GRS and any outcome measures at 1 year follow-up or baseline. CONCLUSION: For the first time we have used all known JIA genetic susceptibility loci (P=<5×10-08) in a GRS to predict changes in disease outcome measured over time. Genetic susceptibility variants are poor predictors of changes in core outcome measures, it is likely that genetic factors predicting disease outcome are independent to those predicting susceptibility. The next step will be to conduct a genome-wide association analysis of JIA outcome.

Extending the allelic spectrum at noncoding risk loci of orofacial clefting.

Genome-wide association studies (GWAS) have generated unprecedented insights into the genetic etiology of orofacial clefting (OFC). The moderate effect sizes of associated noncoding risk variants and limited access to disease-relevant tissue represent considerable challenges for biological interpretation of genetic findings. As rare variants with stronger effect sizes are likely to also contribute to OFC, an alternative approach to delineate pathogenic mechanisms is to identify private mutations and/or an increased burden of rare variants in associated regions. This report describes a framework for targeted resequencing at selected noncoding risk loci contributing to nonsyndromic cleft lip with/without cleft palate (nsCL/P), the most frequent OFC subtype. Based on GWAS data, we selected three risk loci and identified candidate regulatory regions (CRRs) through the integration of credible SNP information, epigenetic data from relevant cells/tissues, and conservation scores. The CRRs (total 57 kb) were resequenced in a multiethnic study population (1061 patients; 1591 controls), using single-molecule molecular inversion probe technology. Combining evidence from in silico variant annotation, pedigree- and burden analyses, we identified 16 likely deleterious rare variants that represent new candidates for functional studies in nsCL/P. Our framework is scalable and represents a promising approach to the investigation of additional congenital malformations with multifactorial etiology.

Combined genetic analysis of juvenile idiopathic arthritis clinical subtypes identifies novel risk loci, target genes and key regulatory mechanisms.

OBJECTIVES: Juvenile idiopathic arthritis (JIA) is the most prevalent form of juvenile rheumatic disease. Our understanding of the genetic risk factors for this disease is limited due to low disease prevalence and extensive clinical heterogeneity. The objective of this research is to identify novel JIA susceptibility variants and link these variants to target genes, which is essential to facilitate the translation of genetic discoveries to clinical benefit. METHODS: We performed a genome-wide association study (GWAS) in 3305 patients and 9196 healthy controls, and used a Bayesian model selection approach to systematically investigate specificity and sharing of associated loci across JIA clinical subtypes. Suggestive signals were followed-up for meta-analysis with a previous GWAS (2751 cases/15 886 controls). We tested for enrichment of association signals in a broad range of functional annotations, and integrated statistical fine-mapping and experimental data to identify target genes. RESULTS: Our analysis provides evidence to support joint analysis of all JIA subtypes with the identification of five novel significant loci. Fine-mapping nominated causal single nucleotide polymorphisms with posterior inclusion probabilities ≥50% in five JIA loci. Enrichment analysis identified RELA and EBF1 as key transcription factors contributing to disease risk. Our integrative approach provided compelling evidence to prioritise target genes at six loci, highlighting mechanistic insights for the disease biology and IL6ST as a potential drug target. CONCLUSIONS: In a large JIA GWAS, we identify five novel risk loci and describe potential function of JIA association signals that will be informative for future experimental works and therapeutic strategies.

Genomic risk scores for juvenile idiopathic arthritis and its subtypes.

OBJECTIVES: Juvenile idiopathic arthritis (JIA) is an autoimmune disease and a common cause of chronic disability in children. Diagnosis of JIA is based purely on clinical symptoms, which can be variable, leading to diagnosis and treatment delays. Despite JIA having substantial heritability, the construction of genomic risk scores (GRSs) to aid or expedite diagnosis has not been assessed. Here, we generate GRSs for JIA and its subtypes and evaluate their performance. METHODS: We examined three case/control cohorts (UK, US-based and Australia) with genome-wide single nucleotide polymorphism (SNP) genotypes. We trained GRSs for JIA and its subtypes using lasso-penalised linear models in cross-validation on the UK cohort, and externally tested it in the other cohorts. RESULTS: The JIA GRS alone achieved cross-validated area under the receiver operating characteristic curve (AUC)=0.670 in the UK cohort and externally-validated AUCs of 0.657 and 0.671 in the US-based and Australian cohorts, respectively. In logistic regression of case/control status, the corresponding odds ratios (ORs) per standard deviation (SD) of GRS were 1.831 (1.685 to 1.991) and 2.008 (1.731 to 2.345), and were unattenuated by adjustment for sex or the top 10 genetic principal components. Extending our analysis to JIA subtypes revealed that the enthesitis-related JIA had both the longest time-to-referral and the subtype GRS with the strongest predictive capacity overall across data sets: AUCs 0.82 in UK; 0.84 in Australian; and 0.70 in US-based. The particularly common oligoarthritis JIA also had a GRS that outperformed those for JIA overall, with AUCs of 0.72, 0.74 and 0.77, respectively. CONCLUSIONS: A GRS for JIA has potential to augment clinical JIA diagnosis protocols, prioritising higher-risk individuals for follow-up and treatment. Consistent with JIA heterogeneity, subtype-specific GRSs showed particularly high performance for enthesitis-related and oligoarthritis JIA.

Using functional genomics to advance the understanding of psoriatic arthritis.

Psoriatic arthritis (PsA) is a complex disease where susceptibility is determined by genetic and environmental risk factors. Clinically, PsA involves inflammation of the joints and the skin, and, if left untreated, results in irreversible joint damage. There is currently no cure and the few treatments available to alleviate symptoms do not work in all patients. Over the past decade, genome-wide association studies (GWAS) have uncovered a large number of disease-associated loci but translating these findings into functional mechanisms and novel targets for therapeutic use is not straightforward. Most variants have been predicted to affect primarily long-range regulatory regions such as enhancers. There is now compelling evidence to support the use of chromatin conformation analysis methods to discover novel genes that can be affected by disease-associated variants. Here, we will review the studies published in the field that have given us a novel understanding of gene regulation in the context of functional genomics and how this relates to the study of PsA and its underlying disease mechanism.

Imputation of orofacial clefting data identifies novel risk loci and sheds light on the genetic background of cleft lip ± cleft palate and cleft palate only.

Nonsyndromic cleft lip with or without cleft palate (nsCL/P) is among the most common human birth defects with multifactorial etiology. Here, we present results from a genome-wide imputation study of nsCL/P in which, after adding replication cohort data, four novel risk loci for nsCL/P are identified (at chromosomal regions 2p21, 14q22, 15q24 and 19p13). On a systematic level, we show that the association signals within this high-density dataset are enriched in functionally-relevant genomic regions that are active in both human neural crest cells (hNCC) and mouse embryonic craniofacial tissue. This enrichment is also detectable in hNCC regions primed for later activity. Using GCTA analyses, we suggest that 30% of the estimated variance in risk for nsCL/P in the European population can be attributed to common variants, with 25.5% contributed to by the 24 risk loci known to date. For each of these, we identify credible SNPs using a Bayesian refinement approach, with two loci harbouring only one probable causal variant. Finally, we demonstrate that there is no polygenic component of nsCL/P detectable that is shared with nonsyndromic cleft palate only (nsCPO). Our data suggest that, while common variants are strongly contributing to risk for nsCL/P, they do not seem to be involved in nsCPO which might be more often caused by rare deleterious variants. Our study generates novel insights into both nsCL/P and nsCPO etiology and provides a systematic framework for research into craniofacial development and malformation.

Focused HLA analysis in Caucasians with myositis identifies significant associations with autoantibody subgroups.

OBJECTIVES: Idiopathic inflammatory myopathies (IIM) are a spectrum of rare autoimmune diseases characterised clinically by muscle weakness and heterogeneous systemic organ involvement. The strongest genetic risk is within the major histocompatibility complex (MHC). Since autoantibody presence defines specific clinical subgroups of IIM, we aimed to correlate serotype and genotype, to identify novel risk variants in the MHC region that co-occur with IIM autoantibodies. METHODS: We collected available autoantibody data in our cohort of 2582 Caucasian patients with IIM. High resolution human leucocyte antigen (HLA) alleles and corresponding amino acid sequences were imputed using SNP2HLA from existing genotyping data and tested for association with 12 autoantibody subgroups. RESULTS: We report associations with eight autoantibodies reaching our study-wide significance level of p<2.9×10-5. Associations with the 8.1 ancestral haplotype were found with anti-Jo-1 (HLA-B*08:01, p=2.28×10-53 and HLA-DRB1*03:01, p=3.25×10-9), anti-PM/Scl (HLA-DQB1*02:01, p=1.47×10-26) and anti-cN1A autoantibodies (HLA-DRB1*03:01, p=1.40×10-11). Associations independent of this haplotype were found with anti-Mi-2 (HLA-DRB1*07:01, p=4.92×10-13) and anti-HMGCR autoantibodies (HLA-DRB1*11, p=5.09×10-6). Amino acid positions may be more strongly associated than classical HLA associations; for example with anti-Jo-1 autoantibodies and position 74 of HLA-DRB1 (p=3.47×10-64) and position 9 of HLA-B (p=7.03×10-11). We report novel genetic associations with HLA-DQB1 anti-TIF1 autoantibodies and identify haplotypes that may differ between adult-onset and juvenile-onset patients with these autoantibodies. CONCLUSIONS: These findings provide new insights regarding the functional consequences of genetic polymorphisms within the MHC. As autoantibodies in IIM correlate with specific clinical features of disease, understanding genetic risk underlying development of autoantibody profiles has implications for future research.

The prevalence of co-morbidities and their impact on physical activity in people with inflammatory rheumatic diseases compared with the general population: results from the UK Biobank.

Objectives: To compare the prevalence and incidence of chronic co-morbidities in people with inflammatory rheumatic and musculoskeletal diseases (iRMDs), and to determine whether the prevalent co-morbidities are associated with physical activity levels in people with iRMDs and in those without iRMDs. Methods: Participants were recruited to the UK Biobank; a population-based cohort. Data were collected about demographics, physical activity, iRMDs (RA, PsA, AS, SLE) and other chronic conditions, including angina, myocardial infarction, stroke, hypertension, pulmonary disease, diabetes and depression. The standardized prevalence of co-morbidities in people with iRMDs was calculated. Cox regression was used to determine the relationship between the presence of an iRMD and an incident co-morbidity. The relationship between the presence (versus absence) of a (co-)morbidity and physical activity level (low, moderate, high) in people with iRMDs and in those without was assessed using multinomial logistic regression. Results: A total of 488 991 participants were included. The estimated prevalence of each co-morbidity was increased in participants with an iRMD, compared with in those without, particularly for stroke in participants with SLE (standardized morbidity ratio (95% CI), 4.9 (3.6, 6.6). Compared with people with no iRMD and no morbidity, the odds ratios (95% CI) for moderate physical activity were decreased for: no iRMD and morbidity, 0.87 (0.85, 0.89); iRMD and no co-morbidity, 0.71 (0.64, 0.80); and iRMD and co-morbidity, 0.58 (0.54, 0.63). Conclusion: Having a (co-)morbidity is associated with reduced physical activity in the general population, and to a greater extent in participants with an iRMD. Optimal management of both iRMDs and co-morbidities may help to reduce their impact on physical activity.

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.

A rare coding allele in IFIH1 is protective for psoriatic arthritis.

OBJECTIVES: Psoriatic arthritis (PsA) is an inflammatory arthritis associated with psoriasis. While many common risk alleles have been reported for association with PsA as well as psoriasis, few rare coding alleles have yet been identified. METHODS: To identify rare coding variation associated with PsA risk or protection, we genotyped 41 267 variants with the exome chip and investigated association within an initial cohort of 1980 PsA cases and 5913 controls. Genotype data for an independent cohort of 2234 PsA cases and 5708 controls was also made available, allowing for a meta-analysis to be performed with the discovery dataset. RESULTS: We identified an association with the rare variant rs35667974 (p=2.39x10-6, OR=0.47), encoding an Ile923Val amino acid change in the IFIH1 gene protein product. The association was reproduced in our independent cohort, which reached a high level of significance on meta-analysis with the discovery and replication datasets (p=4.67x10-10). We identified a strong association with IFIH1 when performing multiple-variant analysis (p=6.77x10-6), and found evidence of independent effects between the rare allele and the common PsA variant at the same locus. CONCLUSION: For the first time, we report a rare coding allele in IFIH1 to be protective for PsA. This rare allele has also been identified to have the same direction of effect on type I diabetes and psoriasis. While this association further supports existing evidence for IFIH1 as a causal gene for PsA, mechanistic studies will need to be pursued to confirm that IFIH1 is indeed causal.

Cross-phenotype association mapping of the MHC identifies genetic variants that differentiate psoriatic arthritis from psoriasis.

OBJECTIVES: Psoriatic arthritis (PsA) is a chronic inflammatory arthritis, with a strong heritable component, affecting patients with psoriasis. Here we attempt to identify genetic variants within the major histocompatibility complex (MHC) that differentiate patients with PsA from patients with cutaneous psoriasis alone (PsC). METHODS: 2808 patients with PsC, 1945 patients with PsA and 8920 population controls were genotyped. We imputed SNPs, amino acids and classical HLA alleles across the MHC and tested for association with PsA compared to population controls and the PsC patient group. In addition we investigated the impact of the age of disease onset on associations. RESULTS: HLA-C*06:02 was protective of PsA compared to PsC (p=9.57×10-66, OR 0.37). The HLA-C*06:02 risk allele was associated with a younger age of psoriasis onset in all patients (p=1.01×10-59). After controlling for the age of psoriasis onset no association of PsA to HLA-C*06:02 (p=0.07) was observed; instead, the most significant association was to amino acid at position 97 of HLA-B (p=1.54×10-9) where the presence of asparagine or serine residue increased PsA risk. Asparagine at position 97 of HLA-B defines the HLA-B*27 alleles. CONCLUSIONS: By controlling for the age of psoriasis onset, we show, for the first time, that HLA-C*06:02 is not associated with PsA and that amino acid position 97 of HLA-B differentiates PsA from PsC. This amino acid also represents the largest genetic effect for ankylosing spondylitis, thereby refining the genetic overlap of these two spondyloarthropathies. Correcting for bias has important implications for cross-phenotype genetic studies.

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.