Conditional meta-analysis stratifying on detailed HLA genotypes identifies a novel type 1 diabetes locus around TCF19 in the MHC

The human leukocyte antigen (HLA) class II genes HLA-DRB1, -DQA1 and -DQB1 are the strongest genetic factors for type 1 diabetes (T1D). Additional loci in the major histocompatibility complex (MHC) are difficult to identify due to the region’s high gene density and complex linkage disequilibrium (LD). To facilitate the association analysis, two novel algorithms were implemented in this study: one for phasing the multi-allelic HLA genotypes in trio families, and one for partitioning the HLA strata in conditional testing. Screening and replication were performed on two large and independent datasets: the Wellcome Trust Case–Control Consortium (WTCCC) dataset of 2,000 cases and 1,504 controls, and the T1D Genetics Consortium (T1DGC) dataset of 2,300 nuclear families. After imputation, the two datasets have 1,941 common SNPs in the MHC, of which 22 were successfully tested and replicated based on the statistical testing stratifying on the detailed DRB1 and DQB1 genotypes. Further conditional tests using the combined dataset confirmed eight novel SNP associations around 31.3 Mb on chromosome 6 (rs3094663, p = 1.66 × 10−11 and rs2523619, p = 2.77 × 10−10 conditional on the DR/DQ genotypes). A subsequent LD analysis established TCF19, POU5F1, CCHCR1 and PSORS1C1 as potential causal genes for the observed association

Investigating the Causal Relationship of C-Reactive Protein with 32 Complex Somatic and Psychiatric Outcomes: A Large-Scale Cross-Consortium Mendelian Randomization Study

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

Additive and interaction effects at three amino acid positions in HLA-DQ and HLA-DR molecules drive type 1 diabetes risk.

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.This research utilizes resources provided by the Type 1 Diabetes Genetics Consortium, a collaborative clinical study sponsored by the National Institute of Diabetes and Digestive and Kidney Diseases (NIDDK), National Institute of Allergy and Infectious Diseases (NIAID), National Human Genome Research Institute (NHGRI), National Institute of Child Health and Human Development (NICHD), and Juvenile Diabetes Research Foundation International (JDRF) and supported by U01 DK062418. This work is supported in part by funding from the National Institutes of Health (5R01AR062886-02 (PIdB), 1R01AR063759 (SR), 5U01GM092691-05 (SR), 1UH2AR067677-01 (SR), R01AR065183 (PIWdB)), a Doris Duke Clinical Scientist Development Award (SR), the Wellcome Trust (JAT) and the National Institute for Health Research (JAT and JMMH), and a Vernieuwingsimpuls VIDI Award (016.126.354) from the Netherlands Organization for Scientific Research (PIWdB). TLL was supported by the German Research Foundation (LE 2593/1-1 and LE 2593/2-1).This is the accepted manuscript. The final version is available at http://www.nature.com/ng/journal/v47/n8/full/ng.3353.html

A Genome-Wide Association Study of the Protein C Anticoagulant Pathway

The Protein C anticoagulant pathway regulates blood coagulation by preventing the inadequate formation of thrombi. It has two main plasma components: protein C and protein S. Individuals with protein C or protein S deficiency present a dramatically increased incidence of thromboembolic disorders. Here, we present the results of a genome-wide association study (GWAS) for protein C and protein S plasma levels in a set of extended pedigrees from the Genetic Analysis of Idiopathic Thrombophilia (GAIT) Project. A total number of 397 individuals from 21 families were typed for 307,984 SNPs using the Infinium® 317 k Beadchip (Illumina). Protein C and protein S (free, functional and total) plasma levels were determined with biochemical assays for all participants. Association with phenotypes was investigated through variance component analysis. After correcting for multiple testing, two SNPs for protein C plasma levels (rs867186 and rs8119351) and another two for free protein S plasma levels (rs1413885 and rs1570868) remained significant on a genome-wide level, located in and around the PROCR and the DNAJC6 genomic regions respectively. No SNPs were significantly associated with functional or total protein S plasma levels, although rs1413885 from DNAJC6 showed suggestive association with the functional protein S phenotype, possibly indicating that this locus plays an important role in protein S metabolism. Our results provide evidence that PROCR and DNAJC6 might play a role in protein C and free protein S plasma levels in the population studied, warranting further investigation on the role of these loci in the etiology of venous thromboembolism and other thrombotic diseases

HLA-DRB1-DQB1 Haplotypes Confer Susceptibility and Resistance to Multiple Sclerosis in Sardinia

Introduction: Genetic predisposition to multiple sclerosis (MS) in Sardinia (Italy) has been associated with five DRB1*-DQB1* haplotypes of the human leukocyte antigen (HLA). Given the complexity of these associations, an in-depth re-analysis was performed with the specific aims of confirming the haplotype associations; establishing the independence of the associated haplotypes; and assessing patients ’ genotypic risk of developing MS. Methods and Results: A transmission disequilibrium test (TDT) of the DRB1*-DQB1 * haplotypes in 943 trio families

Epilepsy, hippocampal sclerosis and febrile seizures linked by common genetic variation around SCN1A

Epilepsy comprises several syndromes, amongst the most common being mesial temporal lobe epilepsy with hippocampal sclerosis. Seizures in mesial temporal lobe epilepsy with hippocampal sclerosis are typically drug-resistant, and mesial temporal lobe epilepsy with hippocampal sclerosis is frequently associated with important co-morbidities, mandating the search for better understanding and treatment. The cause of mesial temporal lobe epilepsy with hippocampal sclerosis is unknown, but there is an association with childhood febrile seizures. Several rarer epilepsies featuring febrile seizures are caused by mutations in SCN1A, which encodes a brain-expressed sodium channel subunit targeted by many anti-epileptic drugs. We undertook a genome-wide association study in 1018 people with mesial temporal lobe epilepsy with hippocampal sclerosis and 7552 control subjects, with validation in an independent sample set comprising 959 people with mesial temporal lobe epilepsy with hippocampal sclerosis and 3591 control subjects. To dissect out variants related to a history of febrile seizures, we tested cases with mesial temporal lobe epilepsy with hippocampal sclerosis with (overall n = 757) and without (overall n = 803) a history of febrile seizures. Meta-analysis revealed a genome-wide significant association for mesial temporal lobe epilepsy with hippocampal sclerosis with febrile seizures at the sodium channel gene cluster on chromosome 2q24.3 [rs7587026, within an intron of the SCN1A gene, P = 3.36 × 10(-9), odds ratio (A) = 1.42, 95% confidence interval: 1.26-1.59]. In a cohort of 172 individuals with febrile seizures, who did not develop epilepsy during prospective follow-up to age 13 years, and 6456 controls, no association was found for rs7587026 and febrile seizures. These findings suggest SCN1A involvement in a common epilepsy syndrome, give new direction to biological understanding of mesial temporal lobe epilepsy with hippocampal sclerosis with febrile seizures, and open avenues for investigation of prognostic factors and possible prevention of epilepsy in some children with febrile seizures

The Oxytocin Receptor (OXTR) Contributes to Prosocial Fund Allocations in the Dictator Game and the Social Value Orientations Task

Background: Economic games observe social decision making in the laboratory that involves real money payoffs. Previously we have shown that allocation of funds in the Dictator Game (DG), a paradigm that illustrates costly altruistic behavior, is partially determined by promoter-region repeat region variants in the arginine vasopressin 1a receptor gene (AVPR1a). In the current investigation, the gene encoding the related oxytocin receptor (OXTR) was tested for association with the DG and a related paradigm, the Social Values Orientation (SVO) task. Methodology/Principal Findings: Association (101 male and 102 female students) using a robust-family based test between 15 single tagging SNPs (htSNPs) across the OXTR was demonstrated with both the DG and SVO. Three htSNPs across the gene region showed significant association with both of the two games. The most significant association was observed with rs1042778 (p = 0.001). Haplotype analysis also showed significant associations for both DG and SVO. Following permutation test adjustment, significance was observed for 2–5 locus haplotypes (p,0.05). A second sample of 98 female subjects was subsequently and independently recruited to play the dictator game and was genotyped for the three significant SNPs found in the first sample. The rs1042778 SNP was shown to be significant for the second sample as well (p = 0.004, Fisher’s exact test). Conclusions: The demonstration that genetic polymorphisms for the OXTR are associated with human prosocial decisio