Rheumatoid arthritis (RA) is the most common inflammatory arthritis and exhibits genetic overlap with other autoimmune and inflammatory disorders. Although predominant associations with the HLA-DRB1 locus have been known for decades, recent data have revealed additional insight into the likely causative variants within HLA-DRB1 as well as within other HLA loci that contribute to disease risk. In addition, more than 100 common variants in non-HLA loci have been implicated in disease susceptibility. Genetic factors are involved not only in the development of RA, but also with various disease subphenotypes, including production and circulating levels of autoantibodies and joint destruction. The major current challenge is to integrate these new data into a precise understanding of disease pathogenesis, including the critical cell types and molecular networks involved as well as interactions with environmental factors. We predict that delineating the functional effects of genetic variants is likely to drive new diagnostic and therapeutic approaches to the disease.

The genetic architecture of age-related macular degeneration (AMD) involves numerous genetic variants, both common and rare, in the coding region of complement factor H (CFH). While these variants explain high disease burden in some families, they fail to explain the pathology in all. We selected families whose AMD was unexplained by known variants and performed whole exome sequencing to probe for other rare, highly penetrant variants. We identified four rare loss-of-function variants in CFH associated with AMD. Missense variant CFH 1:196646753 (C192F) segregated perfectly within a family characterized by advanced AMD and drusen temporal to the macula. Two families, each comprising a pair of affected siblings with extensive extramacular drusen, carried essential splice site variant CFH 1:196648924 (IVS6+1G>A) or missense variant rs139360826 (R175P). In a fourth family, missense variant rs121913058 (R127H) was associated with AMD. Most carriers had early onset bilateral advanced AMD and extramacular drusen. Carriers tended to have low serum Factor H levels, especially carriers of the splice variant. One missense variant (R127H) has been previously shown not to be secreted. The two other missense variants were produced recombinantly: compared to wild type, one (R175P) had no functional activity and the other (C192F) had decreased secretion.

Despite the progress in human leukocyte antigen (HLA) causal variant mapping, independent localization of major histocompatibility complex (MHC) risk from classical HLA genes is challenging. Here, we conducted a large-scale MHC fine-mapping analysis of rheumatoid arthritis (RA) in a Japanese population (6,244 RA cases and 23,731 controls) population by using HLA imputation, followed by a multi-ethnic validation study including east Asian and European populations (n = 7,097 and 23,149, respectively). Our study identified an independent risk of a synonymous mutation at HLA-DOA, a non-classical HLA gene, on anti-citrullinated protein autoantibody (ACPA)-positive RA risk (p = 1.4 × 10(-9)), which demonstrated a cis-expression quantitative trait loci (cis-eQTL) effect on HLA-DOA expression. Trans-ethnic comparison revealed different linkage disequilibrium (LD) patterns in HLA-DOA and HLA-DRB1, explaining the observed HLA-DOA variant risk heterogeneity among ethnicities, which was most evident in the Japanese population. Although previous HLA fine-mapping studies have identified amino acid polymorphisms of the classical HLA genes as driving genetic susceptibility to disease, our study additionally identifies the dosage contribution of a non-classical HLA gene to disease etiology. Our study contributes to the understanding of HLA immunology in human diseases and suggests the value of incorporating additional ancestry in MHC fine-mapping.

BACKGROUND: C-reactive protein (CRP) is associated with immune, cardiometabolic, and psychiatric traits and diseases. Yet it is inconclusive whether these associations are causal. METHODS AND FINDINGS: We performed Mendelian randomization (MR) analyses using two genetic risk scores (GRSs) as instrumental variables (IVs). The first GRS consisted of four single nucleotide polymorphisms (SNPs) in the CRP gene (GRSCRP), and the second consisted of 18 SNPs that were significantly associated with CRP levels in the largest genome-wide association study (GWAS) to date (GRSGWAS). To optimize power, we used summary statistics from GWAS consortia and tested the association of these two GRSs with 32 complex somatic and psychiatric outcomes, with up to 123,865 participants per outcome from populations of European ancestry. We performed heterogeneity tests to disentangle the pleiotropic effect of IVs. A Bonferroni-corrected significance level of less than 0.0016 was considered statistically significant. An observed p-value equal to or less than 0.05 was considered nominally significant evidence for a potential causal association, yet to be confirmed. The strengths (F-statistics) of the IVs were 31.92-3,761.29 and 82.32-9,403.21 for GRSCRP and GRSGWAS, respectively. CRP GRSGWAS showed a statistically significant protective relationship of a 10% genetically elevated CRP level with the risk of schizophrenia (odds ratio [OR] 0.86 [95% CI 0.79-0.94]; p < 0.001). We validated this finding with individual-level genotype data from the schizophrenia GWAS (OR 0.96 [95% CI 0.94-0.98]; p < 1.72 × 10-6). Further, we found that a standardized CRP polygenic risk score (CRPPRS) at p-value thresholds of 1 × 10-4, 0.001, 0.01, 0.05, and 0.1 using individual-level data also showed a protective effect (OR < 1.00) against schizophrenia; the first CRPPRS (built of SNPs with p < 1 × 10-4) showed a statistically significant (p < 2.45 × 10-4) protective effect with an OR of 0.97 (95% CI 0.95-0.99). The CRP GRSGWAS showed that a 10% increase in genetically determined CRP level was significantly associated with coronary artery disease (OR 0.88 [95% CI 0.84-0.94]; p < 2.4 × 10-5) and was nominally associated with the risk of inflammatory bowel disease (OR 0.85 [95% CI 0.74-0.98]; p < 0.03), Crohn disease (OR 0.81 [95% CI 0.70-0.94]; p < 0.005), psoriatic arthritis (OR 1.36 [95% CI 1.00-1.84]; p < 0.049), knee osteoarthritis (OR 1.17 [95% CI 1.01-1.36]; p < 0.04), and bipolar disorder (OR 1.21 [95% CI 1.05-1.40]; p < 0.007) and with an increase of 0.72 (95% CI 0.11-1.34; p < 0.02) mm Hg in systolic blood pressure, 0.45 (95% CI 0.06-0.84; p < 0.02) mm Hg in diastolic blood pressure, 0.01 ml/min/1.73 m2 (95% CI 0.003-0.02; p < 0.005) in estimated glomerular filtration rate from serum creatinine, 0.01 g/dl (95% CI 0.0004-0.02; p < 0.04) in serum albumin level, and 0.03 g/dl (95% CI 0.008-0.05; p < 0.009) in serum protein level. However, after adjustment for heterogeneity, neither GRS showed a significant effect of CRP level (at p < 0.0016) on any of these outcomes, including coronary artery disease, nor on the other 20 complex outcomes studied. Our study has two potential limitations: the limited variance explained by our genetic instruments modeling CRP levels in blood and the unobserved bias introduced by the use of summary statistics in our MR analyses. CONCLUSIONS: Genetically elevated CRP levels showed a significant potentially protective causal relationship with risk of schizophrenia. We observed nominal evidence at an observed p < 0.05 using either GRSCRP or GRSGWAS-with persistence after correction for heterogeneity-for a causal relationship of elevated CRP levels with psoriatic osteoarthritis, rheumatoid arthritis, knee osteoarthritis, systolic blood pressure, diastolic blood pressure, serum albumin, and bipolar disorder. These associations remain yet to be confirmed. We cannot verify any causal effect of CRP level on any of the other common somatic and neuropsychiatric outcomes investigated in the present study. This implies that interventions that lower CRP level are unlikely to result in decreased risk for the majority of common complex outcomes.

We simultaneously investigated the genetic landscape of ankylosing spondylitis, Crohn's disease, psoriasis, primary sclerosing cholangitis and ulcerative colitis to investigate pleiotropy and the relationship between these clinically related diseases. Using high-density genotype data from more than 86,000 individuals of European ancestry, we identified 244 independent multidisease signals, including 27 new genome-wide significant susceptibility loci and 3 unreported shared risk loci. Complex pleiotropy was supported when contrasting multidisease signals with expression data sets from human, rat and mouse together with epigenetic and expressed enhancer profiles. The comorbidities among the five immune diseases were best explained by biological pleiotropy rather than heterogeneity (a subgroup of cases genetically identical to those with another disease, possibly owing to diagnostic misclassification, molecular subtypes or excessive comorbidity). In particular, the strong comorbidity between primary sclerosing cholangitis and inflammatory bowel disease is likely the result of a unique disease, which is genetically distinct from classical inflammatory bowel disease phenotypes.

BACKGROUND: Myocardial mass is a key determinant of cardiac muscle function and hypertrophy. Myocardial depolarization leading to cardiac muscle contraction is reflected by the amplitude and duration of the QRS complex on the electrocardiogram (ECG). Abnormal QRS amplitude or duration reflect changes in myocardial mass and conduction, and are associated with increased risk of heart failure and death. OBJECTIVES: This meta-analysis sought to gain insights into the genetic determinants of myocardial mass. METHODS: We carried out a genome-wide association meta-analysis of 4 QRS traits in up to 73,518 individuals of European ancestry, followed by extensive biological and functional assessment. RESULTS: We identified 52 genomic loci, of which 32 are novel, that are reliably associated with 1 or more QRS phenotypes at p < 1 × 10(-8). These loci are enriched in regions of open chromatin, histone modifications, and transcription factor binding, suggesting that they represent regions of the genome that are actively transcribed in the human heart. Pathway analyses provided evidence that these loci play a role in cardiac hypertrophy. We further highlighted 67 candidate genes at the identified loci that are preferentially expressed in cardiac tissue and associated with cardiac abnormalities in Drosophila melanogaster and Mus musculus. We validated the regulatory function of a novel variant in the SCN5A/SCN10A locus in vitro and in vivo. CONCLUSIONS: Taken together, our findings provide new insights into genes and biological pathways controlling myocardial mass and may help identify novel therapeutic targets.

OBJECTIVES: Juvenile idiopathic arthritis (JIA) is a heterogeneous group of diseases, comprising seven categories. Genetic data could potentially be used to help redefine JIA categories and improve the current classification system. The human leucocyte antigen (HLA) region is strongly associated with JIA. Fine-mapping of the region was performed to look for similarities and differences in HLA associations between the JIA categories and define correspondences with adult inflammatory arthritides. METHODS: Dense genotype data from the HLA region, from the Immunochip array for 5043 JIA cases and 14 390 controls, were used to impute single-nucleotide polymorphisms, HLA classical alleles and amino acids. Bivariate analysis was performed to investigate genetic correlation between the JIA categories. Conditional analysis was used to identify additional effects within the region. Comparison of the findings with those in adult inflammatory arthritic diseases was performed. RESULTS: We identified category-specific associations and have demonstrated for the first time that rheumatoid factor (RF)-negative polyarticular JIA and oligoarticular JIA are genetically similar in their HLA associations. We also observe that each JIA category potentially has an adult counterpart. The RF-positive polyarthritis association at HLA-DRB1 amino acid at position 13 mirrors the association in adult seropositive rheumatoid arthritis (RA). Interestingly, the combined oligoarthritis and RF-negative polyarthritis dataset shares the same association with adult seronegative RA. CONCLUSIONS: The findings suggest the value of using genetic data in helping to classify the categories of this heterogeneous disease. Mapping JIA categories to adult counterparts could enable shared knowledge of disease pathogenesis and aetiology and facilitate transition from paediatric to adult services.

IMPORTANCE: Advances have been made in identifying genetic susceptibility loci for autoimmune diseases, but evidence is needed regarding their association with prognosis and treatment response. OBJECTIVE: To assess whether specific HLA-DRB1 haplotypes associated with rheumatoid arthritis (RA) susceptibility are also associated with radiological severity, mortality, and response to tumor necrosis factor (TNF) inhibitor drugs. DESIGN, SETTING, AND PARTICIPANTS: The Norfolk Arthritis Register (NOAR; 1691 patients and 2811 radiographs; recruitment: 1989-2008; 2008 as final follow-up) was used as a discovery cohort and the Early Rheumatoid Arthritis Study (421 patients and 3758 radiographs; recruitment: 1986-1999; 2005 as final follow-up) as an independent replication cohort for studies of radiographic outcome. Mortality studies were performed in the NOAR cohort (2432 patients; recruitment: 1990-2007; 2011 as final follow-up) and studies of treatment response in the Biologics in Rheumatoid Arthritis Genetics and Genomics Study Syndicate cohort (1846 patients enrolled at initiation of TNF inhibitor; recruitment: 2006-2010; 2011 as final follow-up). Longitudinal statistical modeling was performed to integrate multiple radiograph records per patient over time. All patients were from the United Kingdom and had self-reported white ancestry. EXPOSURES: Sixteen HLA-DRB1 haplotypes defined by amino acids at positions 11, 71, and 74. MAIN OUTCOMES AND MEASURES: Radiological outcome using the Larsen score (range: 0 [none] to 200 [severe joint damage]) and erosions of the hands and feet on radiographs, all-cause mortality, and treatment response measured by change in Disease Activity Score based on 28 joint counts and European League Against Rheumatism (EULAR) response. RESULTS: Patients with RA and valine at position 11 of HLA-DRB1 had the strongest association with radiological damage (OR, 1.75 [95% CI, 1.51-2.05], P = 4.6E-13). By year 5, the percentages of patients with erosions of the hands and feet were 48% of noncarriers (150/314) of valine at position 11, 61% of heterozygote carriers (130/213), and 74% of homozygote carriers (43/58). Valine at position 11 also was associated with higher all-cause mortality in patients with inflammatory polyarthritis (hazard ratio, 1.16 [95% CI, 1.03-1.31], P = .01) (noncarriers: 319 deaths in 1398 patients over 17,196 person-years, mortality rate of 1.9% per year; carriers: 324 deaths in 1116 patients in 13,208 person-years, mortality rate of 2.5% per year) and with better EULAR response to TNF inhibitor therapy (OR, 1.14 [95% CI, 1.01-1.30], P = .04) (noncarriers: 78% [439/561 patients] with moderate or good EULAR response; heterozygote carriers: 81% [698/866]; and homozygote carriers: 86% [277/322]). The risk hierarchy defined by HLA-DRB1 haplotypes was correlated between disease susceptibility, severity, and mortality, but inversely correlated with TNF inhibitor treatment response. CONCLUSIONS AND RELEVANCE: Among patients with RA, the HLA-DRB1 locus, which is associated with disease susceptibility, was also associated with radiological severity, mortality, and treatment response. Replication of these findings in other cohorts is needed as a next step in evaluating the role of HLA-DRB1 haplotype analysis for management of RA.

Variation in the human leukocyte antigen (HLA) genes accounts for one-half of the genetic risk in type 1 diabetes (T1D). Amino acid changes in the HLA-DR and HLA-DQ molecules mediate most of the risk, but extensive linkage disequilibrium complicates the localization of independent effects. Using 18,832 case-control samples, we localized the signal to 3 amino acid positions in HLA-DQ and HLA-DR. HLA-DQβ1 position 57 (previously known; P = 1 × 10(-1,355)) by itself explained 15.2% of the total phenotypic variance. Independent effects at HLA-DRβ1 positions 13 (P = 1 × 10(-721)) and 71 (P = 1 × 10(-95)) increased the proportion of variance explained to 26.9%. The three positions together explained 90% of the phenotypic variance in the HLA-DRB1-HLA-DQA1-HLA-DQB1 locus. Additionally, we observed significant interactions for 11 of 21 pairs of common HLA-DRB1-HLA-DQA1-HLA-DQB1 haplotypes (P = 1.6 × 10(-64)). HLA-DRβ1 positions 13 and 71 implicate the P4 pocket in the antigen-binding groove, thus pointing to another critical protein structure for T1D risk, in addition to the HLA-DQ P9 pocket.

Identifying genomic annotations that differentiate causal from trait-associated variants is essential to fine mapping disease loci. Although many studies have identified non-coding functional annotations that overlap disease-associated variants, these annotations often colocalize, complicating the ability to use these annotations for fine mapping causal variation. We developed a statistical approach (Genomic Annotation Shifter [GoShifter]) to assess whether enriched annotations are able to prioritize causal variation. GoShifter defines the null distribution of an annotation overlapping an allele by locally shifting annotations; this approach is less sensitive to biases arising from local genomic structure than commonly used enrichment methods that depend on SNP matching. Local shifting also allows GoShifter to identify independent causal effects from colocalizing annotations. Using GoShifter, we confirmed that variants in expression quantitative trail loci drive gene-expression changes though DNase-I hypersensitive sites (DHSs) near transcription start sites and independently through 3' UTR regulation. We also showed that (1) 15%-36% of trait-associated loci map to DHSs independently of other annotations; (2) loci associated with breast cancer and rheumatoid arthritis harbor potentially causal variants near the summits of histone marks rather than full peak bodies; (3) variants associated with height are highly enriched in embryonic stem cell DHSs; and (4) we can effectively prioritize causal variation at specific loci.