Lee MN, Ye C, Villani A-C, Raj T, Li W, Eisenhaure TM, Imboywa SH, Chipendo PI, Ran AF, Slowikowski K, Ward LD, Raddassi K, McCabe C, Lee MH, Frohlich IY, Hafler DA, Kellis M, Raychaudhuri S, Zhang F, Stranger BE, Benoist CO, De Jager PL, Regev A, Hacohen N. Common genetic variants modulate pathogen-sensing responses in human dendritic cells [Internet]. Science 2014;343(6175):1246980. Publisher's VersionAbstract
Little is known about how human genetic variation affects the responses to environmental stimuli in the context of complex diseases. Experimental and computational approaches were applied to determine the effects of genetic variation on the induction of pathogen-responsive genes in human dendritic cells. We identified 121 common genetic variants associated in cis with variation in expression responses to Escherichia coli lipopolysaccharide, influenza, or interferon-β (IFN-β). We localized and validated causal variants to binding sites of pathogen-activated STAT (signal transducer and activator of transcription) and IRF (IFN-regulatory factor) transcription factors. We also identified a common variant in IRF7 that is associated in trans with type I IFN induction in response to influenza infection. Our results reveal common alleles that explain interindividual variation in pathogen sensing and provide functional annotation for genetic variants that alter susceptibility to inflammatory diseases.
, Consortium AGENTD (AGEN-T2D)2, Consortium SATD (SAT2D)2, Consortium MATD (MAT2D)2, by sequencing in muylti-Ethnic Consortium TDGEN-generation S (T2D-GENES)2, Mahajan A, Go MJ, Zhang W, Below JE, Gaulton KJ, Ferreira T, Horikoshi M, Johnson AD, Ng MCY, Prokopenko I, Saleheen D, Wang X, Zeggini E, Abecasis GR, Adair LS, Almgren P, Atalay M, Aung T, Baldassarre D, Balkau B, Bao Y, Barnett AH, Barroso I, Basit A, Been LF, Beilby J, Bell GI, Benediktsson R, Bergman RN, Boehm BO, Boerwinkle E, Bonnycastle LL, Burtt N, Cai Q, Campbell H, Carey J, Cauchi S, Caulfield M, Chan JCN, Chang L-C, Chang T-J, Chang Y-C, Charpentier G, Chen C-H, Chen H, Chen Y-T, Chia K-S, Chidambaram M, Chines PS, Cho NH, Cho YM, Chuang L-M, Collins FS, Cornelis MC, Couper DJ, Crenshaw AT, van Dam RM, Danesh J, Das D, de Faire U, Dedoussis G, Deloukas P, Dimas AS, Dina C, Doney AS, Donnelly PJ, Dorkhan M, van Duijn C, Dupuis J, Edkins S, Elliott P, Emilsson V, Erbel R, Eriksson JG, Escobedo J, Esko T, Eury E, Florez JC, Fontanillas P, Forouhi NG, Forsen T, Fox C, Fraser RM, Frayling TM, Froguel P, Frossard P, Gao Y, Gertow K, Gieger C, Gigante B, Grallert H, Grant GB, Grrop LC, Groves CJ, Grundberg E, Guiducci C, Hamsten A, Han B-G, Hara K, Hassanali N, Hattersley AT, Hayward C, Hedman AK, Herder C, Hofman A, Holmen OL, Hovingh K, Hreidarsson AB, Hu C, Hu FB, Hui J, Humphries SE, Hunt SE, Hunter DJ, Hveem K, Hydrie ZI, Ikegami H, Illig T, Ingelsson E, Islam M, Isomaa B, Jackson AU, Jafar T, James A, Jia W, Jöckel K-H, Jonsson A, Jowett JBM, Kadowaki T, Kang HM, Kanoni S, Kao WHL, Kathiresan S, Kato N, Katulanda P, Keinanen-Kiukaanniemi KM, Kelly AM, Khan H, Khaw K-T, Khor C-C, Kim H-L, Kim S, Kim YJ, Kinnunen L, Klopp N, Kong A, Korpi-Hyövälti E, Kowlessur S, Kraft P, Kravic J, Kristensen MM, Krithika S, Kumar A, Kumate J, Kuusisto J, Kwak SH, Laakso M, Lagou V, Lakka TA, Langenberg C, Langford C, Lawrence R, Leander K, Lee J-M, Lee NR, Li M, Li X, Li Y, Liang J, Liju S, Lim W-Y, Lind L, Lindgren CM, Lindholm E, Liu C-T, Liu J, Lobbens S, Long J, Loos RJF, Lu W, Luan J'an, Lyssenko V, Ma RCW, Maeda S, Mägi R, Männisto S, Matthews DR, Meigs JB, Melander O, Metspalu A, Meyer J, Mirza G, Mihailov E, Moebus S, Mohan V, Mohlke KL, Morris AD, Mühleisen TW, Müller-Nurasyid M, Musk B, Nakamura J, Nakashima E, Navarro P, Ng P-K, Nica AC, Nilsson PM, Njølstad I, Nöthen MM, Ohnaka K, Ong TH, Owen KR, Palmer CNA, Pankow JS, Park KS, Parkin M, Pechlivanis S, Pedersen NL, Peltonen L, Perry JRB, Peters A, Pinidiyapathirage JM, Platou CG, Potter S, Price JF, Qi L, Radha V, Rallidis L, Rasheed A, Rathman W, Rauramaa R, Raychaudhuri S, Rayner WN, Rees SD, Rehnberg E, Ripatti S, Robertson N, Roden M, Rossin EJ, Rudan I, Rybin D, Saaristo TE, Salomaa V, Saltevo J, Samuel M, Sanghera DK, Saramies J, Scott J, Scott LJ, Scott RA, Segrè AV, Sehmi J, Sennblad B, Shah N, Shah S, Shera SA, Ou Shu X, Shuldiner AR, Sigurđsson G, Sijbrands E, Silveira A, Sim X, Sivapalaratnam S, Small KS, So WY, Stančáková A, Stefansson K, Steinbach G, Steinthorsdottir V, Stirrups K, Strawbridge RJ, Stringham HM, Sun Q, Suo C, Syvänen A-C, Takayanagi R, Takeuchi F, Tay WT, Teslovich TM, Thorand B, Thorleifsson G, Thorsteinsdottir U, Tikkanen E, Trakalo J, Tremoli E, Trip MD, Tsai FJ, Tuomi T, Tuomilehto J, Uitterlinden AG, Valladares-Salgado A, Vedantam S, Veglia F, Voight BF, Wang C, Wareham NJ, Wennauer R, Wickremasinghe AR, Wilsgaard T, Wilson JF, Wiltshire S, Winckler W, Wong TY, Wood AR, Wu J-Y, Wu Y, Yamamoto K, Yamauchi T, Yang M, Yengo L, Yokota M, Young R, Zabaneh D, Zhang F, Zhang R, Zheng W, Zimmet PZ, Altshuler D, Bowden DW, Cho YS, Cox NJ, Cruz M, Hanis CL, Kooner J, Lee J-Y, Seielstad M, Teo YY, Boehnke M, Parra EJ, Chambers JC, Tai SE, McCarthy MI, Morris AP. Genome-wide trans-ancestry meta-analysis provides insight into the genetic architecture of type 2 diabetes susceptibility [Internet]. Nat Genet 2014;46(3):234-44. Publisher's VersionAbstract
To further understanding of the genetic basis of type 2 diabetes (T2D) susceptibility, we aggregated published meta-analyses of genome-wide association studies (GWAS), including 26,488 cases and 83,964 controls of European, east Asian, south Asian and Mexican and Mexican American ancestry. We observed a significant excess in the directional consistency of T2D risk alleles across ancestry groups, even at SNPs demonstrating only weak evidence of association. By following up the strongest signals of association from the trans-ethnic meta-analysis in an additional 21,491 cases and 55,647 controls of European ancestry, we identified seven new T2D susceptibility loci. Furthermore, we observed considerable improvements in the fine-mapping resolution of common variant association signals at several T2D susceptibility loci. These observations highlight the benefits of trans-ethnic GWAS for the discovery and characterization of complex trait loci and emphasize an exciting opportunity to extend insight into the genetic architecture and pathogenesis of human diseases across populations of diverse ancestry.
Reynolds RJ, Ahmed AF, Danila MI, Hughes LB, for the of with Investigators CLEAAERA, Gregersen PK, Raychaudhuri S, Plenge RM, Bridges LS. HLA-DRB1-associated rheumatoid arthritis risk at multiple levels in African Americans: hierarchical classification systems, amino acid positions, and residues [Internet]. Arthritis Rheumatol 2014;66(12):3274-82. Publisher's VersionAbstract
OBJECTIVE: To evaluate HLA-DRB1 genetic risk of rheumatoid arthritis (RA) in African Americans by 3 validated allele classification systems and by amino acid position and residue, and to compare genetic risk between African American and European ancestries. METHODS: Four-digit HLA-DRB1 genotyping was performed on 561 autoantibody-positive African American cases and 776 African American controls. Association analysis was performed on Tezenas du Montcel (TdM), de Vries (DV), and Mattey classification system alleles and separately by amino acid position and individual residues. RESULTS: TdM S2 and S3P alleles were associated with RA (odds ratio [95% confidence interval] 2.8 [2.0-3.9] and 2.1 [1.7-2.7], respectively). The DV (P = 3.2 × 10(-12)) and Mattey (P = 6.5 × 10(-13)) system alleles were both protective in African Americans. Amino acid position 11 (permutation P < 0.00001) accounted for nearly all variability explained by HLA-DRB1, although conditional analysis demonstrated that position 57 was also significant (0.01 ≤ permutation P ≤ 0.05). The valine and aspartic acid residues at position 11 conferred the highest risk of RA in African Americans. CONCLUSION: With some exceptions, the genetic risk conferred by HLA-DRB1 in African Americans is similar to that in individuals of European ancestry at multiple levels: classification system (e.g., TdM), amino acid position (e.g., 11), and residue (Val11). Unlike that reported for individuals of European ancestry, amino acid position 57 was associated with RA in African Americans, but positions 71 and 74 were not. Asp11 (odds ratio 1 in European ancestry) corresponds to the 4-digit classical allele *09:01, which is also a risk allele for RA in Koreans.
Cui J, Taylor KE, Lee YC, Källberg H, Weinblatt ME, Coblyn JS, Klareskog L, Criswell LA, Gregersen PK, Shadick NA, Plenge RM, Karlson EW. The influence of polygenic risk scores on heritability of anti-CCP level in RA [Internet]. Genes Immun 2014;15(2):107-14. Publisher's VersionAbstract
The objective of this study was to study genetic factors that influence quantitative anticyclic citrullinated peptide (anti-CCP) antibody levels in RA patients. We carried out a genome-wide association study (GWAS) meta-analysis using 1975 anti-CCP+ RA patients from three large cohorts, the Brigham Rheumatoid Arthritis Sequential Study (BRASS), North American Rheumatoid Arthritis Consortium (NARAC) and the Epidemiological Investigation of RA (EIRA). We also carried out a genome-wide complex trait analysis (GCTA) to estimate the heritability of anti-CCP levels. GWAS-meta-analysis showed that anti-CCP levels were most strongly associated with the human leukocyte antigen (HLA) region with a P-value of 2 × 10(-11) for rs1980493. There were 112 SNPs in this region that exceeded the genome-wide significance threshold of 5 × 10(-8), and all were in linkage disequilibrium (LD) with the HLA- DRB1*03 allele with LD r(2) in the range of 0.25-0.88. Suggestive novel associations outside of the HLA region were also observed for rs8063248 (near the GP2 gene) with a P-value of 3 × 10(-7). None of the known RA risk alleles (∼52 loci) were associated with anti-CCP level. Heritability analysis estimated that 44% of anti-CCP variation was attributable to genetic factors captured by GWAS variants. In summary, anti-CCP level is a heritable trait, and HLA-DR3 and GP2 are associated with lower anti-CCP levels.
Brownstein CA, Beggs AH, Homer N, Merriman B, Yu TW, Flannery KC, DeChene ET, Towne MC, Savage SK, Price EN, Holm IA, Luquette LJ, Lyon E, Majzoub J, Neupert P, McCallie D, Szolovits P, Willard HF, Mendelsohn NJ, Temme R, Finkel RS, Yum SW, Medne L, Sunyaev SR, Adzhubey I, Cassa CA, de Bakker PIW, Duzkale H, Dworzyński P, Fairbrother W, Francioli L, Funke BH, Giovanni MA, Handsaker RE, Lage K, Lebo MS, Lek M, Leshchiner I, MacArthur DG, McLaughlin HM, Murray MF, Pers TH, Polak PP, Raychaudhuri S, Rehm HL, Soemedi R, Stitziel NO, Vestecka S, Supper J, Gugenmus C, Klocke B, Hahn A, Schubach M, Menzel M, Biskup S, Freisinger P, Deng M, Braun M, Perner S, Smith RJH, Andorf JL, Huang J, Ryckman K, Sheffield VC, Stone EM, Bair T, Black-Ziegelbein AE, Braun TA, Darbro B, DeLuca AP, Kolbe DL, Scheetz TE, Shearer AE, Sompallae R, Wang K, Bassuk AG, Edens E, Mathews K, Moore SA, Shchelochkov OA, Trapane P, Bossler A, Campbell CA, Heusel JW, Kwitek A, Maga T, Panzer K, Wassink T, Van Daele D, Azaiez H, Booth K, Meyer N, Segal MM, Williams MS, Tromp G, White P, Corsmeier D, Fitzgerald-Butt S, Herman G, Lamb-Thrush D, McBride KL, Newsom D, Pierson CR, Rakowsky AT, Maver A, Lovrečić L, Palandačić A, Peterlin B, Torkamani A, Wedell A, Huss M, Alexeyenko A, Lindvall JM, Magnusson M, Nilsson D, Stranneheim H, Taylan F, Gilissen C, Hoischen A, van Bon B, Yntema H, Nelen M, Zhang W, Sager J, Zhang L, Blair K, Kural D, Cariaso M, Lennon GG, Javed A, Agrawal S, Ng PC, Sandhu KS, Krishna S, Veeramachaneni V, Isakov O, Halperin E, Friedman E, Shomron N, Glusman G, Roach JC, Caballero J, Cox HC, Mauldin D, Ament SA, Rowen L, Richards DR, San Lucas AF, Gonzalez-Garay ML, Caskey TC, Bai Y, Huang Y, Fang F, Zhang Y, Wang Z, Barrera J, Garcia-Lobo JM, González-Lamuño D, Llorca J, Rodriguez MC, Varela I, Reese MG, De La Vega FM, Kiruluta E, Cargill M, Hart RK, Sorenson JM, Lyon GJ, Stevenson DA, Bray BE, Moore BM, Eilbeck K, Yandell M, Zhao H, Hou L, Chen X, Yan X, Chen M, Li C, Yang C, Gunel M, Li P, Kong Y, Alexander AC, Albertyn ZI, Boycott KM, Bulman DE, Gordon PMK, Innes MA, Knoppers BM, Majewski J, Marshall CR, Parboosingh JS, Sawyer SL, Samuels ME, Schwartzentruber J, Kohane IS, Margulies DM. An international effort towards developing standards for best practices in analysis, interpretation and reporting of clinical genome sequencing results in the CLARITY Challenge [Internet]. Genome Biol 2014;15(3):R53. Publisher's VersionAbstract
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.
Han B, Kang EY, Raychaudhuri S, de Bakker PIW, Eskin E. Fast pairwise IBD association testing in genome-wide association studies [Internet]. Bioinformatics 2014;30(2):206-13. Publisher's VersionAbstract
MOTIVATION: Recently, investigators have proposed state-of-the-art Identity-by-descent (IBD) mapping methods to detect IBD segments between purportedly unrelated individuals. The IBD information can then be used for association testing in genetic association studies. One approach for this IBD association testing strategy is to test for excessive IBD between pairs of cases ('pairwise method'). However, this approach is inefficient because it requires a large number of permutations. Moreover, a limited number of permutations define a lower bound for P-values, which makes fine-mapping of associated regions difficult because, in practice, a much larger genomic region is implicated than the region that is actually associated. RESULTS: In this article, we introduce a new pairwise method 'Fast-Pairwise'. Fast-Pairwise uses importance sampling to improve efficiency and enable approximation of extremely small P-values. Fast-Pairwise method takes only days to complete a genome-wide scan. In the application to the WTCCC type 1 diabetes data, Fast-Pairwise successfully fine-maps a known human leukocyte antigen gene that is known to cause the disease. AVAILABILITY: Fast-Pairwise is publicly available at:
Okada Y, Wu D, Trynka G, Raj T, Terao C, Ikari K, Kochi Y, Ohmura K, Suzuki A, Yoshida S, Graham RR, Manoharan A, Ortmann W, Bhangale T, Denny JC, Carroll RJ, Eyler AE, Greenberg JD, Kremer JM, Pappas DA, Jiang L, Yin J, Ye L, Su D-F, Yang J, Xie G, Keystone E, Westra H-J, Esko T, Metspalu A, Zhou X, Gupta N, Mirel D, Stahl EA, Diogo D, Cui J, Liao K, Guo MH, Myouzen K, Kawaguchi T, Coenen MJ, van Riel PLCM, van de Laar MAFJ, Guchelaar H-J, Huizinga TWJ, Dieudé P, Mariette X, Bridges LS, Zhernakova A, Toes REM, Tak PP, Miceli-Richard C, Bang S-Y, Lee H-S, Martin J, Gonzalez-Gay MA, Rodriguez-Rodriguez L, Rantapää-Dahlqvist S, Arlestig L, Choi HK, Kamatani Y, Galan P, Lathrop M, Lathrop M, Lathrop M, Eyre S, Bowes J, Barton A, de Vries N, Moreland LW, Criswell LA, Karlson EW, Taniguchi A, Yamada R, Kubo M, Liu JS, Bae S-C, Worthington J, Padyukov L, Klareskog L, Gregersen PK, Raychaudhuri S, Stranger BE, De Jager PL, Franke L, Visscher PM, Brown MA, Yamanaka H, Mimori T, Takahashi A, Xu H, Behrens TW, Siminovitch KA, Momohara S, Matsuda F, Yamamoto K, Plenge RM. Genetics of rheumatoid arthritis contributes to biology and drug discovery [Internet]. Nature 2014;506(7488):376-81. Publisher's VersionAbstract
A major challenge in human genetics is to devise a systematic strategy to integrate disease-associated variants with diverse genomic and biological data sets to provide insight into disease pathogenesis and guide drug discovery for complex traits such as rheumatoid arthritis (RA). Here we performed a genome-wide association study 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. We discovered 42 novel RA risk loci at a genome-wide level of significance, bringing the total to 101 (refs 2 - 4). We devised an in silico pipeline using established bioinformatics methods based on functional annotation, cis-acting expression quantitative trait loci and pathway analyses--as well as novel methods based on genetic overlap with human primary immunodeficiency, haematological cancer somatic mutations and knockout 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.
Okada Y, Diogo D, Greenberg JD, Mouassess F, Achkar WAL, Fulton RS, Denny JC, Gupta N, Mirel D, Gabriel S, Li G, Kremer JM, Pappas DA, Carroll RJ, Eyler AE, Trynka G, Stahl EA, Cui J, Saxena R, Coenen MJ, Guchelaar H-J, Huizinga TWJ, Dieudé P, Mariette X, Barton A, Canhão H, Fonseca JE, de Vries N, Tak PP, Moreland LW, Bridges LS, Miceli-Richard C, Choi HK, Kamatani Y, Galan P, Lathrop M, Raj T, De Jager PL, Raychaudhuri S, Worthington J, Padyukov L, Klareskog L, Siminovitch KA, Gregersen PK, Mardis ER, Arayssi T, Kazkaz LA, Plenge RM. Integration of sequence data from a Consanguineous family with genetic data from an outbred population identifies PLB1 as a candidate rheumatoid arthritis risk gene [Internet]. PLoS One 2014;9(2):e87645. Publisher's VersionAbstract
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.
Raj T, Rothamel K, Mostafavi S, Ye C, Lee MN, Replogle JM, Feng T, Lee M, Asinovski N, Frohlich I, Imboywa S, Von Korff A, Okada Y, Patsopoulos NA, Davis S, McCabe C, Paik H-il, Srivastava GP, Raychaudhuri S, Hafler DA, Koller D, Regev A, Hacohen N, Mathis D, Benoist C, Stranger BE, De Jager PL. Polarization of the effects of autoimmune and neurodegenerative risk alleles in leukocytes [Internet]. Science 2014;344(6183):519-23. Publisher's VersionAbstract
To extend our understanding of the genetic basis of human immune function and dysfunction, we performed an expression quantitative trait locus (eQTL) study of purified CD4(+) T cells and monocytes, representing adaptive and innate immunity, in a multi-ethnic cohort of 461 healthy individuals. Context-specific cis- and trans-eQTLs were identified, and cross-population mapping allowed, in some cases, putative functional assignment of candidate causal regulatory variants for disease-associated loci. We note an over-representation of T cell-specific eQTLs among susceptibility alleles for autoimmune diseases and of monocyte-specific eQTLs among Alzheimer's and Parkinson's disease variants. This polarization implicates specific immune cell types in these diseases and points to the need to identify the cell-autonomous effects of disease susceptibility variants.
Liao KP, Diogo D, Cui J, Cai T, Okada Y, Gainer VS, Murphy SN, Gupta N, Mirel D, Ananthakrishnan AN, Szolovits P, Shaw SY, Raychaudhuri S, Churchill S, Kohane I, Karlson EW, Plenge RM. Association between low density lipoprotein and rheumatoid arthritis genetic factors with low density lipoprotein levels in rheumatoid arthritis and non-rheumatoid arthritis controls [Internet]. Ann Rheum Dis 2014;73(6):1170-5. Publisher's VersionAbstract
OBJECTIVES: While genetic determinants of low density lipoprotein (LDL) cholesterol levels are well characterised in the general population, they are understudied in rheumatoid arthritis (RA). Our objective was to determine the association of established LDL and RA genetic alleles with LDL levels in RA cases compared with non-RA controls. METHODS: Using data from electronic medical records, we linked validated RA cases and non-RA controls to discarded blood samples. For each individual, we extracted data on: first LDL measurement, age, gender and year of LDL measurement. We genotyped subjects for 11 LDL and 44 non-HLA RA alleles, and calculated RA and LDL genetic risk scores (GRS). We tested the association between each GRS and LDL level using multivariate linear regression models adjusted for age, gender, year of LDL measurement and RA status. RESULTS: Among 567 RA cases and 979 controls, 80% were female and mean age at the first LDL measurement was 55 years. RA cases had significantly lower mean LDL levels than controls (117.2 vs 125.6 mg/dl, respectively, p<0.0001). Each unit increase in LDL GRS was associated with 0.8 mg/dl higher LDL levels in both RA cases and controls (p=3.0×10(-7)). Each unit increase in RA GRS was associated with 4.3 mg/dl lower LDL levels in both groups (p=0.01). CONCLUSIONS: LDL alleles were associated with higher LDL levels in RA. RA alleles were associated with lower LDL levels in both RA cases and controls. As RA cases carry more RA alleles, these findings suggest a genetic basis for epidemiological observations of lower LDL levels in RA.
Arking DE, Pulit SL, Crotti L, van der Harst P, Munroe PB, Koopmann TT, Sotoodehnia N, Rossin EJ, Morley M, Wang X, Johnson AD, Lundby A, Gudbjartsson DF, Noseworthy PA, Eijgelsheim M, Bradford Y, Tarasov KV, Dörr M, Müller-Nurasyid M, Lahtinen AM, Nolte IM, Smith AV, Bis JC, Isaacs A, Newhouse SJ, Evans DS, Post WS, Waggott D, Lyytikäinen L-P, Hicks AA, Eisele L, Ellinghaus D, Hayward C, Navarro P, Ulivi S, Tanaka T, Tester DJ, Chatel S, Gustafsson S, Kumari M, Morris RW, Naluai ÅT, Padmanabhan S, Kluttig A, Strohmer B, Panayiotou AG, Torres M, Knoflach M, Hubacek JA, Slowikowski K, Raychaudhuri S, Kumar RD, Harris TB, Launer LJ, Shuldiner AR, Alonso A, Bader JS, Ehret G, Huang H, Kao LWH, Strait JB, Macfarlane PW, Brown M, Caulfield MJ, Samani NJ, Kronenberg F, Willeit J, Willeit J, Willeit J, Smith GJ, Greiser KH, Meyer Zu Schwabedissen H, Werdan K, Carella M, Zelante L, Heckbert SR, Psaty BM, Rotter JI, Kolcic I, Polašek O, Wright AF, Griffin M, Daly MJ, Daly MJ, Arnar DO, Hólm H, Thorsteinsdottir U, Thorsteinsdottir U, Denny JC, Roden DM, Zuvich RL, Emilsson V, Plump AS, Larson MG, O'Donnell CJ, Yin X, Bobbo M, D'Adamo AP, Iorio A, Sinagra G, Carracedo A, Cummings SR, Nalls MA, Jula A, Kontula KK, Marjamaa A, Oikarinen L, Perola M, Porthan K, Erbel R, Hoffmann P, Jöckel K-H, Kälsch H, Nöthen MM, Nöthen MM, den Hoed M, Loos RJF, Thelle DS, Gieger C, Meitinger T, Perz S, Peters A, Prucha H, Sinner MF, Waldenberger M, de Boer RA, Franke L, van der Vleuten PA, Beckmann BM, Martens E, Bardai A, Hofman N, Wilde AAM, Behr ER, Dalageorgou C, Giudicessi JR, Medeiros-Domingo A, Barc J, Kyndt F, Probst V, Ghidoni A, Insolia R, Hamilton RM, Scherer SW, Brandimarto J, Margulies K, Moravec CE, del Greco M F, Fuchsberger C, O'Connell JR, Lee WK, Watt GCM, Campbell H, Wild SH, El Mokhtari NE, Frey N, Asselbergs FW, Mateo Leach I, Navis G, van den Berg MP, van Veldhuisen DJ, Kellis M, Krijthe BP, Franco OH, Hofman A, Kors JA, Uitterlinden AG, Witteman JCM, Kedenko L, Lamina C, Oostra BA, Abecasis GR, Lakatta EG, Mulas A, Orrú M, Schlessinger D, Uda M, Markus MRP, Völker U, Snieder H, Spector TD, Ärnlöv J, Lind L, Sundström J, Syvänen A-C, Kivimaki M, Kähönen M, Mononen N, Raitakari OT, Viikari JS, Adamkova V, Kiechl S, Brion M, Nicolaides AN, Paulweber B, Haerting J, Dominiczak AF, Nyberg F, Whincup PH, Hingorani AD, Schott J-J, Bezzina CR, Ingelsson E, Ferrucci L, Gasparini P, Wilson JF, Rudan I, Franke A, Mühleisen TW, Pramstaller PP, Lehtimäki TJ, Paterson AD, Parsa A, Liu Y, van Duijn CM, Siscovick DS, Gudnason V, Jamshidi Y, Salomaa V, Felix SB, Sanna S, Ritchie MD, Stricker BH, Stefansson K, Boyer LA, Cappola TP, Olsen JV, Lage K, Schwartz PJ, Kääb S, Chakravarti A, Ackerman MJ, Pfeufer A, de Bakker PIW, Newton-Cheh C. Genetic association study of QT interval highlights role for calcium signaling pathways in myocardial repolarization [Internet]. Nat Genet 2014;46(8):826-36. Publisher's VersionAbstract
The QT interval, an electrocardiographic measure reflecting myocardial repolarization, is a heritable trait. QT prolongation is a risk factor for ventricular arrhythmias and sudden cardiac death (SCD) and could indicate the presence of the potentially lethal mendelian long-QT syndrome (LQTS). Using a genome-wide association and replication study in up to 100,000 individuals, we identified 35 common variant loci associated with QT interval that collectively explain ∼8-10% of QT-interval variation and highlight the importance of calcium regulation in myocardial repolarization. Rare variant analysis of 6 new QT interval-associated loci in 298 unrelated probands with LQTS identified coding variants not found in controls but of uncertain causality and therefore requiring validation. Several newly identified loci encode proteins that physically interact with other recognized repolarization proteins. Our integration of common variant association, expression and orthogonal protein-protein interaction screens provides new insights into cardiac electrophysiology and identifies new candidate genes for ventricular arrhythmias, LQTS and SCD.
Schramm EC, Clark SJ, Triebwasser MP, Raychaudhuri S, Seddon JM, Atkinson JP. Genetic variants in the complement system predisposing to age-related macular degeneration: a review [Internet]. Mol Immunol 2014;61(2):118-25. Publisher's VersionAbstract
Age-related macular degeneration (AMD) is a major cause of visual impairment in the western world. It is characterized by the presence of lipoproteinaceous deposits (drusen) in the inner layers of the retina. Immunohistochemistry studies identified deposition of complement proteins in the drusen as well as in the choroid. In the last decade, genetic studies have linked both common and rare variants in genes of the complement system to increased risk of development of AMD. Here, we review the variants described to date and discuss the functional implications of dysregulation of the alternative pathway of complement in AMD.
Lim ET, Liu YP, Chan Y, Tiinamaija T, Käräjämäki AM, Madsen E, Madsen E, Altshuler DM, Raychaudhuri S, Groop L, Flannick J, Hirschhorn JN, Katsanis N, Daly MJ. A novel test for recessive contributions to complex diseases implicates Bardet-Biedl syndrome gene BBS10 in idiopathic type 2 diabetes and obesity [Internet]. Am J Hum Genet 2014;95(5):509-20. Publisher's VersionAbstract
Rare-variant association studies in common, complex diseases are customarily conducted under an additive risk model in both single-variant and burden testing. Here, we describe a method to improve detection of rare recessive variants in complex diseases termed RAFT (recessive-allele-frequency-based test). We found that RAFT outperforms existing approaches when the variant influences disease risk in a recessive manner on simulated data. We then applied our method to 1,791 Finnish individuals with type 2 diabetes (T2D) and 2,657 matched control subjects. In BBS10, we discovered a rare variant (c.1189A>G [p.Ile397Val]; rs202042386) that confers risk of T2D in a recessive state (p = 1.38 × 10(-6)) and would be missed by conventional methods. Testing of this variant in an established in vivo zebrafish model confirmed the variant to be pathogenic. Taken together, these data suggest that RAFT can effectively reveal rare recessive contributions to complex diseases overlooked by conventional association tests.
Gusev A, Lee HS, Trynka G, Finucane H, Vilhjálmsson BJ, Xu H, Zang C, Ripke S, Bulik-Sullivan B, Stahl E, of the Consortium SWGPG, of the Consortium SWGPG, Kähler AK, Hultman CM, Purcell SM, McCarroll SA, Daly M, Pasaniuc B, Sullivan PF, Neale BM, Wray NR, Raychaudhuri S, Price AL, of the Consortium SWGPG, of the Consortium SWGPG. Partitioning heritability of regulatory and cell-type-specific variants across 11 common diseases [Internet]. Am J Hum Genet 2014;95(5):535-52. Publisher's VersionAbstract
Regulatory and coding variants are known to be enriched with associations identified by genome-wide association studies (GWASs) of complex disease, but their contributions to trait heritability are currently unknown. We applied variance-component methods to imputed genotype data for 11 common diseases to partition the heritability explained by genotyped SNPs (hg(2)) across functional categories (while accounting for shared variance due to linkage disequilibrium). Extensive simulations showed that in contrast to current estimates from GWAS summary statistics, the variance-component approach partitions heritability accurately under a wide range of complex-disease architectures. Across the 11 diseases DNaseI hypersensitivity sites (DHSs) from 217 cell types spanned 16% of imputed SNPs (and 24% of genotyped SNPs) but explained an average of 79% (SE = 8%) of hg(2) from imputed SNPs (5.1× enrichment; p = 3.7 × 10(-17)) and 38% (SE = 4%) of hg(2) from genotyped SNPs (1.6× enrichment, p = 1.0 × 10(-4)). Further enrichment was observed at enhancer DHSs and cell-type-specific DHSs. In contrast, coding variants, which span 1% of the genome, explained <10% of hg(2) despite having the highest enrichment. We replicated these findings but found no significant contribution from rare coding variants in independent schizophrenia cohorts genotyped on GWAS and exome chips. Our results highlight the value of analyzing components of heritability to unravel the functional architecture of common disease.
Pillai NE, Okada Y, Saw W-Y, Ong RT-H, Wang X, Tantoso E, Xu W, Peterson TA, Bielawny T, Ali M, Tay K-Y, Poh W-T, Tan LW-L, Koo S-H, Lim W-Y, Soong R, Wenk M, Raychaudhuri S, Little P, Plummer FA, Lee EJD, Chia K-S, Luo M, de Bakker PIW, Teo Y-Y. Predicting HLA alleles from high-resolution SNP data in three Southeast Asian populations. Hum Mol Genet 2014;23(16):4443-51.Abstract
The major histocompatibility complex (MHC) containing the classical human leukocyte antigen (HLA) Class I and Class II genes is among the most polymorphic and diverse regions in the human genome. Despite the clinical importance of identifying the HLA types, very few databases jointly characterize densely genotyped single nucleotide polymorphisms (SNPs) and HLA alleles in the same samples. To date, the HapMap presents the only public resource that provides a SNP reference panel for predicting HLA alleles, constructed with four collections of individuals of north-western European, northern Han Chinese, cosmopolitan Japanese and Yoruba Nigerian ancestry. Owing to complex patterns of linkage disequilibrium in this region, it is unclear whether the HapMap reference panels can be appropriately utilized for other populations. Here, we describe a public resource for the Singapore Genome Variation Project with: (i) dense genotyping across ∼ 9000 SNPs in the MHC; (ii) four-digit HLA typing for eight Class I and Class II loci, in 96 southern Han Chinese, 89 Southeast Asian Malays and 83 Tamil Indians. This resource provides population estimates of the frequencies of HLA alleles at these eight loci in the three population groups, particularly for HLA-DPA1 and HLA-DPB1 that were not assayed in HapMap. Comparing between population-specific reference panels and a cosmopolitan panel created from all four HapMap populations, we demonstrate that more accurate imputation is obtained with population-specific panels than with the cosmopolitan panel, especially for the Malays and Indians but even when imputing between northern and southern Han Chinese. As with SNP imputation, common HLA alleles were imputed with greater accuracy than low-frequency variants.
Roussos P, Mitchell AC, Voloudakis G, Fullard JF, Pothula VM, Tsang J, Stahl EA, Georgakopoulos A, Ruderfer DM, Charney A, Okada Y, Siminovitch KA, Worthington J, Padyukov L, Klareskog L, Gregersen PK, Plenge RM, Raychaudhuri S, Fromer M, Purcell SM, Brennand KJ, Robakis NK, Schadt EE, Akbarian S, Sklar P. A role for noncoding variation in schizophrenia [Internet]. Cell Rep 2014;9(4):1417-29. Publisher's VersionAbstract
A large portion of common variant loci associated with genetic risk for schizophrenia reside within noncoding sequence of unknown function. Here, we demonstrate promoter and enhancer enrichment in schizophrenia variants associated with expression quantitative trait loci (eQTL). The enrichment is greater when functional annotations derived from the human brain are used relative to peripheral tissues. Regulatory trait concordance analysis ranked genes within schizophrenia genome-wide significant loci for a potential functional role, based on colocalization of a risk SNP, eQTL, and regulatory element sequence. We identified potential physical interactions of noncontiguous proximal and distal regulatory elements. This was verified in prefrontal cortex and -induced pluripotent stem cell-derived neurons for the L-type calcium channel (CACNA1C) risk locus. Our findings point to a functional link between schizophrenia-associated noncoding SNPs and 3D genome architecture associated with chromosomal loopings and transcriptional regulation in the brain.
Dunstan SJ, Hue NT, Han B, Li Z, Tram TTB, Sim KS, Parry CM, Chinh NT, Vinh H, Lan NPH, Thieu NTV, Vinh PV, Koirala S, Dongol S, Arjyal A, Karkey A, Shilpakar O, Dolecek C, Foo JN, Phuong LT, Lanh MN, Do T, Aung T, Hon DN, Teo YY, Hibberd ML, Anders KL, Okada Y, Raychaudhuri S, Simmons CP, Baker S, de Bakker PIW, Basnyat B, Hien TT, Farrar JJ, Khor CC. Variation at HLA-DRB1 is associated with resistance to enteric fever. Nat Genet 2014;46(12):1333-6.Abstract
Enteric fever affects more than 25 million people annually and results from systemic infection with Salmonella enterica serovar Typhi or Paratyphi pathovars A, B or C(1). We conducted a genome-wide association study of 432 individuals with blood culture-confirmed enteric fever and 2,011 controls from Vietnam. We observed strong association at rs7765379 (odds ratio (OR) for the minor allele = 0.18, P = 4.5 × 10(-10)), a marker mapping to the HLA class II region, in proximity to HLA-DQB1 and HLA-DRB1. We replicated this association in 595 enteric fever cases and 386 controls from Nepal and also in a second independent collection of 151 cases and 668 controls from Vietnam. Imputation-based fine-mapping across the extended MHC region showed that the classical HLA-DRB1*04:05 allele (OR = 0.14, P = 2.60 × 10(-11)) could entirely explain the association at rs7765379, thus implicating HLA-DRB1 as a major contributor to resistance against enteric fever, presumably through antigen presentation.
Yu Y, Triebwasser MP, Wong EKS, Schramm EC, Thomas B, Reynolds R, Mardis ER, Atkinson JP, Daly M, Raychaudhuri S, Kavanagh D, Seddon JM. Whole-exome sequencing identifies rare, functional CFH variants in families with macular degeneration [Internet]. Hum Mol Genet 2014;23(19):5283-93. Publisher's VersionAbstract
We sequenced the whole exome of 35 cases and 7 controls from 9 age-related macular degeneration (AMD) families in whom known common genetic risk alleles could not explain their high disease burden and/or their early-onset advanced disease. Two families harbored novel rare mutations in CFH (R53C and D90G). R53C segregates perfectly with AMD in 11 cases (heterozygous) and 1 elderly control (reference allele) (LOD = 5.07, P = 6.7 × 10(-7)). In an independent cohort, 4 out of 1676 cases but none of the 745 examined controls or 4300 NHBLI Exome Sequencing Project (ESP) samples carried the R53C mutation (P = 0.0039). In another family of six siblings, D90G similarly segregated with AMD in five cases and one control (LOD = 1.22, P = 0.009). No other sample in our large cohort or the ESP had this mutation. Functional studies demonstrated that R53C decreased the ability of FH to perform decay accelerating activity. D90G exhibited a decrease in cofactor-mediated inactivation. Both of these changes would lead to a loss of regulatory activity, resulting in excessive alternative pathway activation. This study represents an initial application of the whole-exome strategy to families with early-onset AMD. It successfully identified high impact alleles leading to clearer functional insight into AMD etiopathogenesis.
Okada Y, Han B, Tsoi LC, Stuart PE, Ellinghaus E, Tejasvi T, Chandran V, Pellett F, Pollock R, Bowcock AM, Krueger GG, Weichenthal M, Voorhees JJ, Rahman P, Gregersen PK, Franke A, Nair RP, Abecasis GR, Gladman DD, Elder JT, de Bakker PIW, Raychaudhuri S. Fine mapping major histocompatibility complex associations in psoriasis and its clinical subtypes [Internet]. Am J Hum Genet 2014;95(2):162-72. Publisher's VersionAbstract
Psoriasis vulgaris (PsV) risk is strongly associated with variation within the major histocompatibility complex (MHC) region, but its genetic architecture has yet to be fully elucidated. Here, we conducted a large-scale fine-mapping study of PsV risk in the MHC region in 9,247 PsV-affected individuals and 13,589 controls of European descent by imputing class I and II human leukocyte antigen (HLA) genes from SNP genotype data. In addition, we imputed sequence variants for MICA, an MHC HLA-like gene that has been associated with PsV, to evaluate association at that locus as well. We observed that HLA-C(∗)06:02 demonstrated the lowest p value for overall PsV risk (p = 1.7 × 10(-364)). Stepwise analysis revealed multiple HLA-C(∗)06:02-independent risk variants in both class I and class II HLA genes for PsV susceptibility (HLA-C(∗)12:03, HLA-B amino acid positions 67 and 9, HLA-A amino acid position 95, and HLA-DQα1 amino acid position 53; p < 5.0 × 10(-8)), but no apparent risk conferred by MICA. We further evaluated risk of two major clinical subtypes of PsV, psoriatic arthritis (PsA; n = 3,038) and cutaneous psoriasis (PsC; n = 3,098). We found that risk heterogeneity between PsA and PsC might be driven by HLA-B amino acid position 45 (Pomnibus = 2.2 × 10(-11)), indicating that different genetic factors underlie the overall risk of PsV and the risk of specific PsV subphenotypes. Our study illustrates the value of high-resolution HLA and MICA imputation for fine mapping causal variants in the MHC.
Han B, Diogo D, Eyre S, Kallberg H, Zhernakova A, Bowes J, Padyukov L, Okada Y, González-Gay MA, Rantapää-Dahlqvist S, Martin J, Huizinga TWJ, Plenge RM, Worthington J, Gregersen PK, Klareskog L, de Bakker PIW, Raychaudhuri S. Fine mapping seronegative and seropositive rheumatoid arthritis to shared and distinct HLA alleles by adjusting for the effects of heterogeneity [Internet]. Am J Hum Genet 2014;94(4):522-32. Publisher's VersionAbstract
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.