Publications by Year: 2016

2016

Li Z, Xia Y, Feng LN, Chen JR, Li HM, Cui J, Cai QQ, Sim KS, Nairismägi ML, Laurensia Y, Meah WY, Liu WS, Guo YM, Chen LZ, Feng QS, Pang CP, Chen LJ, Chew SH, Ebstein R, Foo JN, Liu J, Ha J, Khoo LP, Chin ST, Zeng YX, Aung T, Chowbay B, Diong CP, Zhang F, Liu YH, Tang T, Tao M, Quek R, Mohamad F, Tan SY, Teh BT, Ng SB, Chng WJ, Ong CK, Okada Y, Raychaudhuri S, Lim ST, Tan W, Peng RJ, Khor CC, Bei JX. Genetic risk of extranodal natural killer T-cell lymphoma: a genome-wide association study. Lancet Oncol. 2016;17(9):1240–7.
BACKGROUND: Extranodal natural killer T-cell lymphoma (NKTCL), nasal type, is a rare and aggressive malignancy that occurs predominantly in Asian and Latin American populations. Although Epstein-Barr virus infection is a known risk factor, other risk factors and the pathogenesis of NKTCL are not well understood. We aimed to identify common genetic variants affecting individual risk of NKTCL. METHODS: We did a genome-wide association study of 189 patients with extranodal NKTCL, nasal type (WHO classification criteria; cases) and 957 controls from Guangdong province, southern China. We validated our findings in four independent case-control series, including 75 cases from Guangdong province and 296 controls from Hong Kong, 65 cases and 983 controls from Guangdong province, 125 cases and 1110 controls from Beijing (northern China), and 60 cases and 2476 controls from Singapore. We used imputation and conditional logistic regression analyses to fine-map the associations. We also did a meta-analysis of the replication series and of the entire dataset. FINDINGS: Associations exceeding the genome-wide significance threshold (p<5 × 10(-8)) were seen at 51 single-nucleotide polymorphisms (SNPs) mapping to the class II MHC region on chromosome 6, with rs9277378 (located in HLA-DPB1) having the strongest association with NKTCL susceptibility (p=4·21 × 10(-19), odds ratio [OR] 1·84 [95% CI 1·61-2·11] in meta-analysis of entire dataset). Imputation-based fine-mapping across the class II MHC region suggests that four aminoacid residues (Gly84-Gly85-Pro86-Met87) in near-complete linkage disequilibrium at the edge of the peptide-binding groove of HLA-DPB1 could account for most of the association between the rs9277378*A risk allele and NKTCL susceptibility (OR 2·38, p value for haplotype 2·32 × 10(-14)). This association is distinct from MHC associations with Epstein-Barr virus infection. INTERPRETATION: To our knowledge, this is the first time that a genetic variant conferring an NKTCL risk is noted at genome-wide significance. This finding underlines the importance of HLA-DP antigen presentation in the pathogenesis of NKTCL. FUNDING: Top-Notch Young Talents Program of China, Special Support Program of Guangdong, Specialized Research Fund for the Doctoral Program of Higher Education (20110171120099), Program for New Century Excellent Talents in University (NCET-11-0529), National Medical Research Council of Singapore (TCR12DEC005), Tanoto Foundation Professorship in Medical Oncology, New Century Foundation Limited, Ling Foundation, Singapore National Cancer Centre Research Fund, and the US National Institutes of Health (1R01AR062886, 5U01GM092691-04, and 1R01AR063759-01A1).
Terao C, Raychaudhuri S, Gregersen P. Recent Advances in Defining the Genetic Basis of Rheumatoid Arthritis. Annu Rev Genomics Hum Genet. 2016;17:273–301.
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.
Wagner E, Raychaudhuri S, Villalonga M, Java A, Triebwasser M, Daly M, Atkinson J, Seddon J. Mapping rare, deleterious mutations in Factor H: Association with early onset, drusen burden, and lower antigenic levels in familial AMD. Sci Rep. 2016;6:31531.
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.
Okada Y, Suzuki A, Ikari K, Terao C, Kochi Y, Ohmura K, Higasa K, Akiyama M, Ashikawa K, Kanai M, Hirata J, Suita N, Teo YY, Xu H, Bae SC, Takahashi A, Momozawa Y, Matsuda K, Momohara S, Taniguchi A, Yamada R, Mimori T, Kubo M, Brown M, Raychaudhuri S, Matsuda F, Yamanaka H, Kamatani Y, Yamamoto K. Contribution of a Non-classical HLA Gene, HLA-DOA, to the Risk of Rheumatoid Arthritis. Am J Hum Genet. 2016;99(2):366–74.
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.
Prins B, Abbasi A, Wong A, Vaez A, Nolte I, Franceschini N, Stuart P, Guterriez Achury J, Mistry V, Bradfield J, Valdes A, Bras J, Shatunov A, PAGE Consortium, International Stroke Genetics Consortium, Systemic Sclerosis consortium, Treat OA consortium, DIAGRAM Consortium, CARDIoGRAMplusC4D Consortium, ALS consortium, International Parkinson’s Disease Genomics Consortium, Consortium ASDWGPG, CKDGen consortium, GERAD1 Consortium, Pressure ICB, Consortium SWGPG, Consortium IWGC, Lu C, Han B, Raychaudhuri S, Bevan S, Mayes M, Tsoi L, Evangelou E, Nair R, Grant S, Polychronakos C, Radstake T, Heel D, Dunstan M, Wood N, Al-Chalabi A, Dehghan A, Hakonarson H, Markus H, Elder J, Knight J, Arking D, Spector T, Koeleman B, Duijn C, Martin J, Morris A, Weersma R, Wijmenga C, Munroe P, Perry J, Pouget J, Jamshidi Y, Snieder H, Alizadeh B. Investigating the Causal Relationship of C-Reactive Protein with 32 Complex Somatic and Psychiatric Outcomes: A Large-Scale Cross-Consortium Mendelian Randomization Study. PLoS Med. 2016;13(6):e1001976.
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.
Ellinghaus D, Jostins L, Spain S, Cortes A, Bethune J, Han B, Park YR, Raychaudhuri S, Pouget J, Hübenthal M, Folseraas T, Wang Y, Esko T, Metspalu A, Westra HJ, Franke L, Pers T, Weersma R, Collij V, D’Amato M, Halfvarson J, Jensen AB, Lieb W, Degenhardt F, Forstner A, Hofmann A, International IBD Genetics Consortium (IIBDGC), (IGAS) IGASC, International PSC Study Group (IPSCSG), (GAPC) GAPC, Psoriasis Association Genetics Extension (PAGE), Schreiber S, Mrowietz U, Juran B, Lazaridis K, Brunak S, Dale A, Trembath R, Weidinger S, Weichenthal M, Ellinghaus E, Elder J, Barker J, Andreassen O, McGovern D, Karlsen T, Barrett J, Parkes M, Brown M, Franke A. Analysis of five chronic inflammatory diseases identifies 27 new associations and highlights disease-specific patterns at shared loci. Nat Genet. 2016;48(5):510–8.
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.
Hinks, Bowes, Cobb, Ainsworth, Marion, Comeau, Sudman, Han, Juvenile Arthritis Consortium for Immunochip, Becker, Bohnsack, Bakker, Haas, Hazen, Lovell, Nigrovic, Nordal, Punnaro, Rosenberg, Rygg, Smith, Wise, Videm, Wedderburn, Yarwood, Yeung, Prahalad, Langefeld, Raychaudhuri, Thompson, Thomson. Fine-mapping the MHC locus in juvenile idiopathic arthritis (JIA) reveals genetic heterogeneity corresponding to distinct adult inflammatory arthritic diseases. Ann Rheum Dis. 2016;76(4):765–772.
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.