Person:

Raychaudhuri, Soumya

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 (P({combined}) = 1.2×10(^{−12})), rs864537 near CD247 (P({combined}) = 2.2×10(^{−11})), rs2298428 near UBE2L3 (P({combined}) = 2.5×10(^{−10})), and rs11203203 near UBASH3A (P({combined}) = 1.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.

Bone mineral density (BMD) is the most important predictor of fracture risk. We performed the largest meta-analysis to date on lumbar spine and femoral neck BMD, including 17 genome-wide association studies and 32,961 individuals of European and East Asian ancestry. We tested the top-associated BMD markers for replication in 50,933 independent subjects and for risk of low-trauma fracture in 31,016 cases and 102,444 controls. We identified 56 loci (32 novel)associated with BMD atgenome-wide significant level (P<5×10−8). Several of these factors cluster within the RANK-RANKL-OPG, mesenchymal-stem-cell differentiation, endochondral ossification and the Wnt signalling pathways. However, we also discovered loci containing genes not known to play a role in bone biology. Fourteen BMD loci were also associated with fracture risk (P<5×10−4, Bonferroni corrected), of which six reached P<5×10−8 including: 18p11.21 (C18orf19), 7q21.3 (SLC25A13), 11q13.2 (LRP5), 4q22.1 (MEPE), 2p16.2 (SPTBN1) and 10q21.1 (DKK1). These findings shed light on the genetic architecture and pathophysiological mechanisms underlying BMD variation and fracture susceptibility.

Anti-tumor necrosis factor alpha (anti-TNF) biologic therapy is a widely used treatment for rheumatoid arthritis (RA). It is unknown why some RA patients fail to respond adequately to anti-TNF therapy, which limits the development of clinical biomarkers to predict response or new drugs to target refractory cases. To understand the biological basis of response to anti-TNF therapy, we conducted a genome-wide association study (GWAS) meta-analysis of more than 2 million common variants in 2,706 RA patients from 13 different collections. Patients were treated with one of three anti-TNF medications: etanercept (n = 733), infliximab (n = 894), or adalimumab (n = 1,071). We identified a SNP (rs6427528) at the 1q23 locus that was associated with change in disease activity score (ΔDAS) in the etanercept subset of patients (P = 8×10−8), but not in the infliximab or adalimumab subsets (P>0.05). The SNP is predicted to disrupt transcription factor binding site motifs in the 3′ UTR of an immune-related gene, CD84, and the allele associated with better response to etanercept was associated with higher CD84 gene expression in peripheral blood mononuclear cells (P = 1×10−11 in 228 non-RA patients and P = 0.004 in 132 RA patients). Consistent with the genetic findings, higher CD84 gene expression correlated with lower cross-sectional DAS (P = 0.02, n = 210) and showed a non-significant trend for better ΔDAS in a subset of RA patients with gene expression data (n = 31, etanercept-treated). A small, multi-ethnic replication showed a non-significant trend towards an association among etanercept-treated RA patients of Portuguese ancestry (n = 139, P = 0.4), but no association among patients of Japanese ancestry (n = 151, P = 0.8). Our study demonstrates that an allele associated with response to etanercept therapy is also associated with CD84 gene expression, and further that CD84 expression correlates with disease activity. These findings support a model in which CD84 genotypes and/or expression may serve as a useful biomarker for response to etanercept treatment in RA patients of European ancestry.

Through genome-wide association meta-analyses of up to 133,010 individuals of European ancestry without diabetes, including individuals newly genotyped using the Metabochip, we have raised the number of confirmed loci influencing glycemic traits to 53, of which 33 also increase type 2 diabetes risk (q < 0.05). Loci influencing fasting insulin showed association with lipid levels and fat distribution, suggesting impact on insulin resistance. Gene-based analyses identified further biologically plausible loci, suggesting that additional loci beyond those reaching genome-wide significance are likely to represent real associations. This conclusion is supported by an excess of directionally consistent and nominally significant signals between discovery and follow-up studies. Functional follow-up of these newly discovered loci will further improve our understanding of glycemic control.

A major challenge in human genetics is to devise a systematic strategy to integrate disease-associated variants with diverse genomic and biological datasets to provide insight into disease pathogenesis and guide drug discovery for complex traits such as rheumatoid arthritis (RA)1. Here, we performed a genome-wide association study (GWAS) 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 (SNPs). We discovered 42 novel RA risk loci at a genome-wide level of significance, bringing the total to 1012–4. We devised an in-silico pipeline using established bioinformatics methods based on functional annotation5, cis-acting expression quantitative trait loci (cis-eQTL)6, and pathway analyses7–9 – as well as novel methods based on genetic overlap with human primary immunodeficiency (PID), hematological cancer somatic mutations and knock-out 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.

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.

Summary: We created a fast, robust and general C++ implementation of a single-nucleotide polymorphism (SNP) set enrichment algorithm to identify cell types, tissues and pathways affected by risk loci. It tests trait-associated genomic loci for enrichment of specificity to conditions (cell types, tissues and pathways). We use a non-parametric statistical approach to compute empirical P-values by comparison with null SNP sets. As a proof of concept, we present novel applications of our method to four sets of genome-wide significant SNPs associated with red blood cell count, multiple sclerosis, celiac disease and HDL cholesterol. Availability and implementation: http://broadinstitute.org/mpg/snpsea Contact: soumya@broadinstitute.org Supplementary information: Supplementary data are available at Bioinformatics online.

Integrating genetic data from families with highly penetrant forms of disease together with genetic data from outbred populations represents a promising strategy to uncover the complete frequency spectrum of risk alleles for complex traits such as rheumatoid arthritis (RA). Here, we demonstrate that rare, low-frequency and common alleles at one gene locus, phospholipase B1 (PLB1), might contribute to risk of RA in a 4-generation consanguineous pedigree (Middle Eastern ancestry) and also in unrelated individuals from the general population (European ancestry). Through identity-by-descent (IBD) mapping and whole-exome sequencing, we identified a non-synonymous c.2263G>C (p.G755R) mutation at the PLB1 gene on 2q23, which significantly co-segregated with RA in family members with a dominant mode of inheritance (P = 0.009). We further evaluated PLB1 variants and risk of RA using a GWAS meta-analysis of 8,875 RA cases and 29,367 controls of European ancestry. We identified significant contributions of two independent non-coding variants near PLB1 with risk of RA (rs116018341 [MAF = 0.042] and rs116541814 [MAF = 0.021], combined P = 3.2×10−6). Finally, we performed deep exon sequencing of PLB1 in 1,088 RA cases and 1,088 controls (European ancestry), and identified suggestive dispersion of rare protein-coding variant frequencies between cases and controls (P = 0.049 for C-alpha test and P = 0.055 for SKAT). Together, these data suggest that PLB1 is a candidate risk gene for RA. Future studies to characterize the full spectrum of genetic risk in the PLB1 genetic locus are warranted.

To extend understanding of the genetic architecture and molecular basis of type 2 diabetes (T2D), we conducted a meta-analysis of genetic variants on the Metabochip involving 34,840 cases and 114,981 controls, overwhelmingly of European descent. We identified ten previously unreported T2D susceptibility loci, including two demonstrating sex-differentiated association. Genome-wide analyses of these data are consistent with a long tail of further common variant loci explaining much of the variation in susceptibility to T2D. Exploration of the enlarged set of susceptibility loci implicates several processes, including CREBBP-related transcription, adipocytokine signalling and cell cycle regulation, in diabetes pathogenesis.

Background: There is tremendous potential for genome sequencing to improve clinical diagnosis and care once it becomes routinely accessible, but this will require formalizing research methods into clinical best practices in the areas of sequence data generation, analysis, interpretation and reporting. The CLARITY Challenge was designed to spur convergence in methods for diagnosing genetic disease starting from clinical case history and genome sequencing data. DNA samples were obtained from three families with heritable genetic disorders and genomic sequence data were donated by sequencing platform vendors. The challenge was to analyze and interpret these data with the goals of identifying disease-causing variants and reporting the findings in a clinically useful format. Participating contestant groups were solicited broadly, and an independent panel of judges evaluated their performance. Results: A total of 30 international groups were engaged. The entries reveal a general convergence of practices on most elements of the analysis and interpretation process. However, even given this commonality of approach, only two groups identified the consensus candidate variants in all disease cases, demonstrating a need for consistent fine-tuning of the generally accepted methods. There was greater diversity of the final clinical report content and in the patient consenting process, demonstrating that these areas require additional exploration and standardization. Conclusions: The CLARITY Challenge provides a comprehensive assessment of current practices for using genome sequencing to diagnose and report genetic diseases. There is remarkable convergence in bioinformatic techniques, but medical interpretation and reporting are areas that require further development by many groups.