Genetics of Low Polygenic Risk Score Type 1 Diabetes Patients: rare variants in 22 novel loci

Jingchun Qu, Hui-Qi Qu, Jonathan Bradfield, Joseph Glessner, Xiao Chang, Lifeng Tian, Michael March, Jeffrey D Roizen, Patrick Sleiman, Hakon Hakonarson
doi: https://doi.org/10.1101/2020.10.13.20211987

Abstract

With polygenic risk score (PRS) for autoimmune type 1 diabetes (T1D), this study identified T1D cases with low T1D PRS and searched for susceptibility loci in these cases. Our hypothesis is that genetic effects (likely mediated by relatively rare genetic variants) of non-mainstream (or non-autoimmune) T1D might have been diluted in the previous studies on T1D cases in general. Two cohorts for the PRS modeling and testing respectively were included. The first cohort consisted of 3,356 T1D cases and 6,203 controls, and the independent second cohort consisted of 3,355 T1D cases and 6,203 controls. Cases with low T1D PRS were identified using PRSice-2 and compared to controls with low T1D PRS by genome-wide association (GWA) test. Twenty-six genetic loci with SNPs/SNVs associated with low PRS T1D at genome-wide significance (P≤5.0xE-08) were identified, including 4 established T1D loci, as well as 22 novel loci represented by rare SNVs. For the 22 novel loci, 12 regions have been reported of association with obesity related traits by previous GWA studies. Five loci encoding long intergenic non-protein coding RNAs (lncRNA), two loci involved in N-linked glycosylation, two loci encoding GTPase activators, and two ciliopathy genes, are also highlighted in this study.

  • Genome-wide Association Study
  • N-linked glycosylation
  • Non-autoimmune
  • Long intergenic non-protein coding RNA
  • Polygenic Risk Score
  • Type 1 Diabetes

Introduction

Type 1 diabetes (T1D) is caused by T-cell mediated autoimmune destruction of pancreatic β-cells(1). There is no cure for T1D to date. The molecular mechanisms underlying T1D are complex and not completely understood. Human genetic studies have uncovered multiple T1D genes that contribute to our understanding of the pathogenesis ofT1D(2-7). With the rapid advances in human genomics technology in recent years, over 70 T1D loci have been identified(8) (https://www.ebi.ac.uk/gwas/). While these discoveries of T1D-associated genes have greatly increased our knowledge of T1D, our current genetic knowledge on T1D is far from complete, and a large number of T1D genes remain uncovered(9). A key bottleneck for the GWAS approach is limitation of sample size even with the presense of collaborative international consortia(10). The phenotype of type 1 diabetes has been regarded as heterogeneous. While the majority of T1D patients have autoimmune disease, 5–10% of Caucasian diabetic subjects with recent-onset T1D do not have islet cell antibodies, often referred to as T1bD(11). Due to different pathogenesis, T1bD cases may be associated with different genetic loci from autoimmune T1D, or T1aD. However, the smaller proportion of T1bD cases suggests that T1bD-related genetic effects have been diluted in the previous studies with T1D cases studied in general. Besides T1bD, the non-autoimmune and monogenic form of pediatric diabetes, maturity-onset diabetes of the young (MODY) cases, may be misdiagnosed as T1D(12), which further contributes to the heterogeneity of the T1D phenotype.

With numerous genetic loci for many human complex diseases identified to date, polygenic risk scores (PRS) aggregate the effects of many genetic variants across the human genome into a single score, an approach that has been shown of improve disease prediction and differential diagnosis(13). The T1D loci identified by the GWAS studies to date are mainly associated with the genetic susceptibility of the major component of the heterogeneous T1D phenotype, i.e. T1aD, while the genetic susceptibility of the minor non-autoimmune components (e.g. T1bD and misdiagnosed MODY) are undere-represented in those results likely as a result of being diluted In this study, we propose that a high T1D PRS score predicts or suggests a T1aD case, whereas a low T1D PRS score in a T1D case suggests the opposite and represents our major interest in this study. Our aim in this study is to identify low PRS T1D cases and to run a seprate GWAS in an attempt to uncover genetic loci associated with T1bD patients.

Methods

Subjects: 6,711 European T1D cases and 12,406 European controls were included in this study. The T1D cases were from the Children’s Hospital of Philadelphia (CHOP)(14), The Diabetes Control and Complications Trial – Epidemiology of Diabetes Interventions and Complications (DCCT-EDIC) cohort (http://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000086.v2.p1), the Type 1 Diabetes Genetics Consortium (T1DGC, http://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000180.v1.p1), and later recruited subjects at CHOP, respectively. The genotyping was done by the Illumina Human Hap550 Genotyping BeadChip or a newer version of Illumina Genotyping BeadChip. Other demographic, phenotypic and genotypic details about these individuals were described in our previous publication(15). Imputation of 39,131,579 single nucleotide polymorphisms (SNP) on auto-chromosomes was done using the Sanger Imputation Service (https://www.sanger.ac.uk/tool/sanger-imputation-service/) based on the Haplotype Reference Consortium (HRC) r1.1 reference panel (HRC.r1-1.GRCh37.wgs.mac5.sites.tab), with the quality filters of R2 ≥ 0.4. Altogether, 32,251,301 autosomal single nucleotide variants (SNV) with quality R2 ≥ 0.4 were included in this study. Population stratification was assessed by principal component analysis (PCA), and genetic association tests were corrected by the first 10 principal components (PC). The association test was done using PLINK1.9 software(16).

Polygenic risk scores (PRS): To avoid the issue of overfitting for PRS scoring, the subjects were randomly splitted into two independent cohorts without duplication, i.e. the PRS training cohort including 3,356 T1D cases and 6,203 controls, and the PRS testing cohort including 3,355 T1D cases and 6,203 controls. PRSs of the test cohort were calculated using the Polygenic Risk Score software (PRSice-2)(17), based on the statistics of the training group. The performance of a series of cutoff of T1D association P-values (including 10−10, 10−9, 10−8, 10−7, 10−6, 10−5, 10−4, 0.001, 0.01, 0.05, 0.1, 0.2, and 1) for selection of SNP markers was assessed by the Area Under the ROC Curve (AUC). The P-value cutoff with the largest AUC was adopted.

GWAS of T1D patients with low PRS: According to the PRS values, the T1D patients were separated into two groups, i.e. a low PRS group and a high PRS group. The PRS cutoff was determined by the maximum Matthews correlation coefficient (MCC). Using the same PRS cutoff, health controls with low T1D PRS were identified. The GWAS of T1D patients with low PRS was performed by comparing to health controls with low T1D PRS. The Manhattan plots were done using the web-based FUMA platform(18). Genetic association signals within each locus were plotted by LocusZoom(19).

Data and Resource Availability

The datasets generated during and/or analyzed during the current study are available from the corresponding author upon reasonable request.

Results

AUC of different cutoffs of T1D association P-values for SNP selection and PRS

The AUCs of different cutoffs of T1D association P-values for selection of SNP sets are shown in Table 1a. The best AUC (0.8607) is seen at the cutoff of P-value≤1E-05, which suggests that stricter cutoff may cause the missing of informative SNPs, while looser may introduce noise by including SNPs with spurious T1D association. Based on the SNP markers with T1D association P-value≤1E-05, a PRS score was acquired for each individual in the independent test cohort. By the maximum MCC (Supplementary Table 1), a PRS cutoff of 1.11E-03 has the maximum MCC (0.6294). A PRS≤1.11E-03 was defined as low risk, and a PRS>1.11E-03 was defined as high risk. With this threshold, the sensitivity (True positive rate, TPR) for T1D prediction is 75.9%, and the specificity (True negative rate, TFR) for T1D prediction is 86.4%. By PRS≤1.11E-03, 810 (24.1%, including 408 males, 400 females, and 2 cases with undetermined sex) out of 3,355 T1D cases had low PRS; and 5,358 (86.4%, including 2,893 males, 2,453 females, and 12 cases with undetermined sex) out of 6,203 controls had low PRS.

View this table:
Table 1

The AUCs of different cutoffs of T1D association P-values

GWAS of T1D patients with low PRS

The GWAS of T1D patients with low T1D PRS compared to controls with low T1D PRS identified a large number of SNPs associated with T1D with genome-wide significance (P≤5.0xE-08), from 7 genetic loci (Supplementary Table 2, Figure 1). Among these 7 genetic loci, 3 loci have been established of T1D association by previous studies, including HLA, INS, and PTPN22 (Table 2a). By looking at the established leading T1D signal of each locus, the frequencies of the predisposing alleles of HLA and PTPN22 were lower in the low T1D PRS cohort, while the protective allele of INS were higher in the low T1D PRS cohort. The effect sizes of HLA (P=2.72E-06) and PTPN22 (P=0.047) were significantly smaller in the low PRS cases. Besides these 3 established T1D loci, 4 novel loci associated with low PRS T1D were identified (Table 3a). LocusZoom plots for genetic association signals within each locus are shown in Supplementary Figure 1-4. The association signals of these loci are only seen in low PRS T1D cases, but not in the T1D cases overall, and were missed previously due to diluted genetic effects.

View this table:
Table 2

Leading SNPs at three loci have been established of T1D association

View this table:
Table 3

Novel loci associated with low PRS T1D

Figure 1.
Figure 1.

The Manhattan plots of the first cohort. (a) The plot of the GWAS of T1D patients with low T1D PRS compared to controls with low T1D PRS (810 cases vs. 5358 controls); (b) The plot of the GWAS of all T1D patients compared to all controls (3355 cases vs. 6203 controls).

Replication of the PRS model and additional novel loci

Consequently, we switched the two cohorts, i.e. using the second cohort for the statistics of PRS modelling, then we tested the PRS models in the first cohort. The AUCs of different cutoffs of T1D association P-values for selection of SNP sets are shown in Table 1b. The best AUC (0.8654) is seen at the cutoff of P-value≤1E-05, which repeated the PRS model in the above step. Based on the SNP markers with T1D association P-value≤1E-05, a PRS score was acquired for each individual in the independent test cohort. By the maximum MCC (Supplementary Table 3), a PRS cutoff of 1.24E-03 has the maximum MCC (0.6294). A PRS≤7.18E-04 was defined as low risk, and a PRS>7.18E-04 was defined as high risk. With this threshold, the sensitivity (True positive rate, TPR) for T1D prediction is 66.0%, and the specificity (True negative rate, TFR) for T1D prediction is 93.6%. By PRS≤7.18E-04, 918 (27.4%, including 437 males, 479 females, and 2 cases with undetermined sex) out of 3,356 T1D cases had low PRS; and 5,585 (90.0%, including 3,008 males, 2,565 females, and 12 cases with undetermined sex) out of 6,203 controls had low PRS.

As expected from the above results, in the switched cohort, the GWAS of T1D patients with low T1D PRS compared to controls with low T1D PRS identified a large number of SNPs associated with T1D with genome-wide significance (P≤5.0xE-08) as well (Supplementary Table 4, Figure 2). Among these loci, 4 loci have been established of T1D association by previous studies, including HLA, INS, PTPN22, and the IKZF4/RPS26/ERBB3 locus (Table 2b). Consistent to the first GWAS results listed above, by looking at the established leading T1D signal of each locus, the frequencies of the predisposing alleles of HLA, PTPN22 and IKZF4 were lower in the low T1D PRS cohort, while the protective allele of INS were higher in the low T1D PRS cohort. The effect size of the leading HLA SNP was significantly smaller in the low PRS cases (P=1.05E-11). Besides these established T1D loci, 18 novel loci associated with low PRS T1D were identified in this cohort (Table 3b). LocusZoom plots for genetic association signals within each locus are shown in Supplementary Figure 5-22.

Figure 2.
Figure 2.

The Manhattan plots of the second cohort. (a) The plot of the GWAS of T1D patients with low T1D PRS compared to controls with low T1D PRS (918 cases vs. 5585 controls); (b) The plot of the GWAS of all T1D patients compared to all controls (3356 cases vs. 6203 controls).

Discussion

Altogether, rare variants (MAF<5%) from 22 novel loci were identified in the low PRS T1D cases with genome-wide significance (P<5.00E-08), in addition to the 4 established T1D loci with smaller genetic effects in these cases. The association signals of these loci are only seen in low PRS T1D cases, but not in the T1D cases overall, and were missed previously due to rare allele frequencies and diluted genetic effects in the general T1D cohort. Among the 22 loci, two genetic regions have been reported of association with diabetes, i.e. the region containing the DLL1/ FAM120B locus associated with type 1 diabetes in Caucasian by our previous study(20), and the region containing the TICRR locus associated Type 2 diabetes in African population(21). In addition, a number of genetic associations with body mass index (BMI), obesity, and autoimmunity, have been reported in the flanking regions of 300kb on each side of the new loci according to the GWAS Catalog (https://www.ebi.ac.uk/gwas/, Supplementary materials). Further details on these 22 loci are described below.

LINC01865/LINC01874 tagged by rs186500234

The long intergenic non-protein coding RNA 1865 gene (LINC01865) has low expression observed in testis, brain, and duodenum. The long intergenic non-protein coding RNA 1874 gene (LINC01874) has restricted expression toward kidney(22). This genetic region has been reported of association with body mass index (BMI) by previous study(23).

LOC730100 tagged by rs28957087

LOC730100 encodes a long non-coding RNA (ncRNA), a competing endogenous RNA for human microRNA 760 (miR-760)(24). The latter inhibits the expression of the Forkhead Box A1 gene (FOXA1). As a hepatocyte nuclear factor, FOXA1, also known as HNF3A or TCF3A, regulates tissue-specific gene expression in liver and many other tissues(25). FoxA1 is essential for normal pancreatic and ß-cell function and a negative regulator of the hepatocyte nuclear factor-1 (HNF1) homeobox A gene (HNF1A) and the hepatocyte nuclear factor 4, alpha gene (HNF4A)(26) (27). HNF1A and HNF4A are established genes causing maturity-onset diabetes of the young (MODY). The FOXA1 mutation Ser448Asn has been suggested of association with impaired glucose homeostasis(27).

B3GNT2/TMEM17 tagged by rs75634056

The UDP-GlcNAc:betaGal beta-1,3-N-acetylglucosaminyltransferase 2 gene (B3GNT2) encodes an enzyme involved in the biosynthesis of poly-N-acetyllactosamine chains. The gene plays an important role in immunological biofunctions, and its deficiency causes hyperactivation of lymphocytes in mice(28). The transmembrane protein 17 gene (TMEM17) encodes a critical component of a protein complex at the base of cilia. Previous GWAS studies have reported association with Crohn’s disease, ankylosing spondylitis, and hypothyroidism in this genetic region.

FAM136A/TGFA tagged by rs77418738

The family with sequence similarity 136 member A gene (FAM136A) encodes a mitochondrially localized protein. The transforming growth factor alpha gene (TGFA) mediates cell-cell adhesion and activates cell proliferation, differentiation and development. This region has been reported of association with obesity-related traits(29).

GCC2/EDAR tagged by rs922452

The GRIP and coiled-coil domain containing 2 (GCC2) encodes a long coiled-coil protein, also known as GCC185, which is localized to the trans-Golgi network with critical function in maintaining Golgi structure and tethering transport vesicle(30). The ectodysplasin A receptor gene (EDAR) encodes a member of the tumor necrosis factor receptor family with a key role in ectodermal differentiation. Association with low birth weight at this region has been reported(31).

SEL1L3 tagged by rs6842426

The locus SEL1L family member 3 gene (SEL1L3) is a paralog of the SEL1L adaptor subunit of ERAD E3 ubiquitin ligase gene (SEL1L). The latter is highly expressed in pancreas and thyroid, and is crucial for misfolded proteins in the endoplasmic reticulum being discharged into the cytosol and degraded by the proteasome(32). This gene region has been reported association with obesity-related traits(29) and non-alcoholic fatty liver disease(33).

TBC1D1/LINC01258 tagged by rs4833044

This genetic locus contains two genes, the TBC1 domain family member 1 gene (TBC1D1) and the long intergenic non-protein coding RNA 1258 gene (LINC01258). A common variant it this locus has been reported to be associated with childhood obesity(29; 34), triacylglycerol 54:5 levels(35), lymphocyte percentage of leukocytes(36) by previous studies. Acting as a GTPase activator, the TBC1D1 protein plays a role in regulating cell growth and differentiation. Rare mutations in TBC1D1 have been reported to be associated with congenital anomalies of the kidney and urinary tract.

LINC02432/IL15 tagged by rs9790756

The long intergenic non-protein coding RNA 2432 gene (LINC02432) has higher expression in kidney and pancreas. Interleukin 15 (IL-15) encoded by the gene IL15 is essential for regulating activation and proliferation of T and natural killer cells, and supporting lymphoid homeostasis. IL-15 and interleukine 2 (IL-2) share many biological activities and receptor components with IL-2. IL-2 is a powerful growth factor for both T and B lymphocytes. Both IL2 and the α chain of the IL2 receptor complex gene (IL2RA) has been established of genetic association with T1D by previous studies(37-39).

DEK/RNF144B tagged by rs16880565

The DEK proto-oncogene gene (DEK) encodes a site-specific DNA binding protein and a component of the pre-mRNA splicing complex, and is involved in transcriptional regulation and pre-mRNA splicing. DEK encoded protein is also an autoantigen in patients with pauciarticular onset juvenile rheumatoid arthritis. The ring finger protein 144B gene (RNF144B) encoded protein inhibits LPS-induced inflammatory responses by binding with TANK binding kinase 1 (TBK1) and causing interferon regulatory factor 3 (IRF3) dephosphorylation and interferon β (IFN-β) reduction. This region has been reported of association with BMI by previous studies(23; 40).

RGS17 tagged by rs80292134

The regulator of G protein signaling 17 gene (RGS17) encodes a member of the regulator of G-protein signaling family. This genetic region has been established association with BMI by previous studies(23; 40).

DLL1/FAM120B tagged by rs3800237

The delta like canonical Notch ligand 1 (DLL1) encodes a Notch ligand with a role in cell-fate decision processes in lymphopoiesis. This Notch ligand can completely inhibit the differentiation of human hematopoietic progenitors into the B cell lineage while promoting the generation of T cell/natural killer (NK) precursors(41). The family with sequence similarity 120B gene (FAM120B) encodes a constitutive coactivator of peroxisome proliferator-activated receptor γ (PPARγ, a major therapeutic target for insulin sensitivity) and promotes adipogenesis(42). The region containing the DLL1/ FAM120B genes has been reported of association with T1D in Caucasian by our previous study(20).

NME8/GPR141 tagged by rs12532321

The NME/NM23 family member 8 gene (NME8) encodes an axoneme protein, and its mutation may cause primary ciliary dyskinesia. The G protein-coupled receptor 141 gene (GPR141) at the upstream of NME8 is highly expressed in bone marrow. This genetic region has been reported of association with obesity-related traits in Hispanic children(29).

CALN1 tagged by rs118182411

The calneuron 1 gene (CALN1), encoding a protein with high similarity to the calcium-binding proteins of calmodulin, is highly expressed in brain and adrenal. This genetic region has established association with BMI by previous studies(23; 40).

ZNF804B tagged by rs77205087

The zinc finger protein 804B gene (ZNF804B) has been reported of association with N-linked glycosylation of human immunoglobulin G (IgG), which modulates its binding to Fc receptors(43). N-glycosylation of cytokines and proteases is also a regulatory mechanism in inflammation and autoimmunity(44). Changes in N-glycosylation have been associated with different autoimmune diseases, including rheumatoid arthritis(45), type 1 diabetes(46), Crohn’s disease(47).

NFIB tagged by rs10961435

The nuclear factor I B gene (NFIB) encodes a transcription factor in the FOXA1 transcription factor network. NFIB has been shown to play critical roles in lung and brain development. A previous study has shown that NFIB can bind with FoxA1 and modulate the transcriptional activity of FoxA1(48), while the later has been suggested to play a role in pancreatic and ß-cell function and non-autoimmune diabetes as discussed above.

TBC1D2/GABBR2 tagged by rs11559334

This genetic locus contains two protein-coding genes, the TBC1 domain family member 2 gene (TBC1D2) and the gamma-aminobutyric acid type B receptor subunit 2 gene (GABBR2, encoding a member of the G-protein coupled receptor 3 family). As discussed above, this study identified an association signal in the TBC1D1 region, and the TBC1D1 locus has been reported of association with childhood obesity(29; 34).

LINC00841/C10orf142 tagged by rs746298

The two genes at this locus, LINC00841/C10orf142, encode two long intergenic non-protein coding RNAs (lincRNA). While the function of these two genes remain unknown, this locus has been reported of association with obesity-related traits(29).

SYT10/ALG10 tagged by rs10506114

The synaptotagmin 10 gene (SYT10) encodes a membrane protein of secretory vesicles expressed in pancreas, lung and kidney(49). The ALG10 alpha-1,2-glucosyltransferase gene (ALG10) encodes a membrane-associated protein that adds the third glucose residue to the lipid-linked oligosaccharide precursor for N-glycosylation in endoplasmic reticulum (ER)(50). As discussed above in the ZNF804B locus, N-glycosylation of IgG, cytokines and proteases is also a regulatory mechanism in inflammation and autoimmunity(43; 44) associated with different autoimmune diseases. This region has established association with waist-hip ratio by previous study(40).

CHST11 tagged by rs75438334

The carbohydrate sulfotransferase 11 gene (CHST11) encodes a member of the sulfotransferase 2 family catalyzing chondroitin sulfate synthesis. This genetic region has been reported of association with waist circumference adjusted for body mass index by previous study(51).

CHFR/LOC101928530/ZNF605 tagged by rs12230138

The checkpoint with forkhead and ring finger domains gene (CHFR) encodes an E3 ubiquitin-protein ligase and is involved in the DNA damage response and checkpoint regulation. The structure and function of the gene LOC101928530 is still uncharacterized. The function of the zinc finger protein 605 gene (ZNF605) may be related to Herpes Simplex Virus 1 infection (https://pathcards.genecards.org/card/herpes_simplex_virus_1_infection). This region has been reported of association with BMI by previous study(52).

TICRR/ KIF7 tagged by rs2197053

The TOPBP1 interacting checkpoint and replication regulator gene (TICRR) encodes Treslin, which is involved in triggering the initiation of DNA replication. The kinesin family member 7 gene (KIF7) in this region encodes a cilia-associated protein of the kinesin family, with its mutations causing ciliopathies. The region containing the TICRR gene has been reported of association with T2D in African population(21), BMI(23; 52), and obesity-related traits(29) by previous studies.

LINC01695/LINC00161 tagged by rs7278151

Function of the long intergenic non-protein coding RNA 1695 gene (LINC01695) is still uncharacterized. The long intergenic non-protein coding RNA 161 gene (LINC00161) encodes a functional RNA that regulates Mitogen-activated protein kinase 1 (MAPK1) expression. The MAPK1/STAT3 pathway has been proposed as a novel diabetes target for its critical role in glucose homeostasis(53).

In summary, in the genetic regions containing the 22 novel loci disclosed by this study, more than half of these regions have been reported of association with obesity-related traits, BMI, or waist circumference. The correlation with obesity related traits or impaired glucose homeostasis is in keeping with non-autoimmune roles in the diabetes patients with low T1D PRS. Interestingly, genes related N-linked glycosylation, e.g. ZNF804B and ALG10, are highlighted in this study, which may suggest the role of N-glycosylation bridging impaired glucose homeostasis and autoimmune diabetes. N-glycosylation is commonly altered in diabetes(54). This particular locus supports an interesting hypothesis of T1D pathogenesis, i.e. the accelerator hypothesis, which implies that increasing obesity-associated insulin resistance accelerates the disease process of type 1 diabetes(55; 56). Insulin resistance-related mechanisms might thus be able to serve as potential novel therapeutic targets for these patients with low T1D PRS.

In addition, 5 loci encoding long intergenic non-protein coding RNAs (lncRNA) identified in this study emphasize the importance of lncRNAs in these diabetes patients. This study identified 2 loci containing TBC1D1 and TBC1D2 respectively, encoding two GTPase activators. TBC1D1 has been suggested as a novel obesity gene by previous study(34). Two loci containing the TMEM17 and KIF7 genes corrected with ciliopathies suggest a role of primary cilia in diabetes(57). However, we admit that this study has limitations related to the bottleneck of sample size and data resources. The novel loci reported in this study still need replication in independent samples. In addition, the functional mechanisms of these genetic loci in diabetes warrant experimental investigation.

Data Availability

The datasets generated during and/or analyzed during the current study are available from the corresponding author upon reasonable request.

Competing interests

none to declare.

Acknowledgement

The authors apologize that many important references listed in the supplementary materials cannot be cited in the main text because of page limitation. The study was supported by Institutional Development Funds from the Children’s Hospital of Philadelphia to the Center for Applied Genomics and The Children’s Hospital of Philadelphia Endowed Chair in Genomic Research to HH. Dr. Hakon Hakonarson is the guarantor of this work and, as such, had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.

Reference

  1. 1.↵
    Atkinson MA, Eisenbarth GS, Michels AW: Type 1 diabetes. The Lancet 2014;383:69–82
    OpenUrl
  2. 2.↵
    Todd JA, Bell JI, McDevitt HO: HLA-DQ[beta] gene contributes to susceptibility and resistance to insulin-dependent diabetes mellitus. Nature 1987;329:599–604
    OpenUrlCrossRefPubMedWeb of Science
  3. 3.
    Baisch JM, Weeks T, Giles R, Hoover M, Stastny P, Capra JD: Analysis of HLA-DQ genotypes and susceptibility in insulin-dependent diabetes mellitus. N Engl J Med 1990;322:1836–1841
    OpenUrlCrossRefPubMedWeb of Science
  4. 4.
    Todd JA: Genetic Analysis of Type 1 Diabetes Using Whole Genome Approaches. PNAS 1995;92:8560–8565
    OpenUrlAbstract/FREE Full Text
  5. 5.
    Noble JA, Valdes AM, Cook M, Klitz W, Thomson G, Erlich HA: The role of HLA class II genes in insulin-dependent diabetes mellitus: molecular analysis of 180 Caucasian, multiplex families. Am J Hum Genet 1996;59:1134–1148
    OpenUrlPubMedWeb of Science
  6. 6.
    she J-X: Susceptibility to type I diabetes: HLA-DQ and DR revisited. Immunology Today 1996;17:323
    OpenUrlCrossRefPubMedWeb of Science
  7. 7.↵
    Bell GI, Horita S, Karam JH: A polymorphic locus near the human insulin gene is associated with insulin-dependent diabetes mellitus. Diabetes 1984;33:176–183
    OpenUrlAbstract/FREE Full Text
  8. 8.↵
    Onengut-Gumuscu S, Chen W-M, Burren O, Cooper NJ, Quinlan AR, Mychaleckyj JC, Farber E, Bonnie JK, Szpak M, Schofield E: Fine mapping of type 1 diabetes susceptibility loci and evidence for colocalization of causal variants with lymphoid gene enhancers. Nature genetics 2015;47:381–386
    OpenUrlCrossRefPubMed
  9. 9.↵
    Polychronakos C, Li Q: Understanding type 1 diabetes through genetics: advances and prospects. Nature Reviews Genetics 2011;12:781–792
    OpenUrlCrossRefPubMed
  10. 10.↵
    Rich SS, Concannon P, Erlich H, Julier C, Morahan G, Nerup J, Pociot F, Todd JA: The type 1 diabetes genetics consortium. Annals of the New York Academy of Sciences 2006;1079:1–8
    OpenUrlCrossRefPubMedWeb of Science
  11. 11.↵
    Leslie RD, Atkinson MA, Notkins AL: Autoantigens IA-2 and GAD in Type I (insulin-dependent) diabetes. Diabetologia 1999;42:3–14
    OpenUrlCrossRefPubMedWeb of Science
  12. 12.↵
    Ehtisham S, Hattersley A, Dunger D, Barrett T: First UK survey of paediatric type 2 diabetes and MODY. Archives of disease in childhood 2004;89:526–529
    OpenUrlAbstract/FREE Full Text
  13. 13.↵
    Lambert SA, Abraham G, Inouye M: Towards clinical utility of polygenic risk scores. Human Molecular Genetics 2019;28:R133–R142
    OpenUrlCrossRefPubMed
  14. 14.↵
    Hakonarson H, Grant SF, Bradfield JP, Marchand L, Kim CE, Glessner JT, Grabs R, Casalunovo T, Taback SP, Frackelton EC, Lawson ML, Robinson LJ, Skraban R, Lu Y, Chiavacci RM, Stanley CA, Kirsch SE, Rappaport EF, Orange JS, Monos DS, Devoto M, Qu HQ, Polychronakos C: A genome-wide association study identifies KIAA0350 as a type 1 diabetes gene. Nature 2007;448:591–594
    OpenUrlCrossRefPubMedWeb of Science
  15. 15.↵
    Bradfield JP, Qu H-Q, Wang K, Zhang H, Sleiman PM, Kim CE, Mentch FD, Qiu H, Glessner JT, Thomas KA: A genome-wide meta-analysis of six type 1 diabetes cohorts identifies multiple associated loci. PLoS genetics 2011;7:e1002293
    OpenUrl
  16. 16.↵
    Chang CC, Chow CC, Tellier LC, Vattikuti S, Purcell SM, Lee JJ: Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience 2015;4:s13742-13015-10047-13748
  17. 17.↵
    Choi SW, O’Reilly PF: PRSice-2: Polygenic Risk Score software for biobank-scale data. GigaScience 2019;8
  18. 18.↵
    Watanabe K, Taskesen E, Van Bochoven A, Posthuma D: Functional mapping and annotation of genetic associations with FUMA. Nature communications 2017;8:1–11
    OpenUrlCrossRef
  19. 19.↵
    Pruim RJ, Welch RP, Sanna S, Teslovich TM, Chines PS, Gliedt TP, Boehnke M, Abecasis GR, Willer CJ: LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics (Oxford, England) 2010;26:2336–2337
    OpenUrlCrossRefPubMedWeb of Science
  20. 20.↵
    Bradfield JP, Qu HQ, Wang K, Zhang H, Sleiman PM, Kim CE, Mentch FD, Qiu H, Glessner JT, Thomas KA, Frackelton EC, Chiavacci RM, Imielinski M, Monos DS, Pandey R, Bakay M, Grant SF, Polychronakos C, Hakonarson H: A genome-wide meta-analysis of six type 1 diabetes cohorts identifies multiple associated loci. PLoS Genet 2011;7:e1002293
    OpenUrlCrossRefPubMed
  21. 21.↵
    Chen J, Sun M, Adeyemo A, Pirie F, Carstensen T, Pomilla C, Doumatey AP, Chen G, Young EH, Sandhu M, Morris AP, Barroso I, McCarthy MI, Mahajan A, Wheeler E, Rotimi CN, Motala AA: Genome-wide association study of type 2 diabetes in Africa. Diabetologia 2019;62:1204–1211
    OpenUrl
  22. 22.↵
    Fagerberg L, Hallström BM, Oksvold P, Kampf C, Djureinovic D, Odeberg J, Habuka M, Tahmasebpoor S, Danielsson A, Edlund K, Asplund A, Sjöstedt E, Lundberg E, Szigyarto CA, Skogs M, Takanen JO, Berling H, Tegel H, Mulder J, Nilsson P, Schwenk JM, Lindskog C, Danielsson F, Mardinoglu A, Sivertsson A, von Feilitzen K, Forsberg M, Zwahlen M, Olsson I, Navani S, Huss M, Nielsen J, Ponten F, Uhlén M: Analysis of the human tissue-specific expression by genome-wide integration of transcriptomics and antibody-based proteomics. Molecular & cellular proteomics : MCP 2014;13:397–406
    OpenUrl
  23. 23.↵
    Zhu Z, Guo Y, Shi H, Liu CL, Panganiban RA, Chung W, O’Connor LJ, Himes BE, Gazal S, Hasegawa K, Camargo CA, Jr., Qi L, Moffatt MF, Hu FB, Lu Q, Cookson WOC, Liang L: Shared genetic and experimental links between obesity-related traits and asthma subtypes in UK Biobank. The Journal of allergy and clinical immunology 2020;145:537–549
    OpenUrl
  24. 24.↵
    Li Q, Lu J, Xia J, Wen M, Wang C: Long non-coding RNA LOC730100 enhances proliferation and invasion of glioma cells through competitively sponging miR-760 from FOXA1 mRNA. Biochemical and biophysical research communications 2019;512:558–563
    OpenUrl
  25. 25.↵
    Lee CS, Friedman JR, Fulmer JT, Kaestner KH: The initiation of liver development is dependent on Foxa transcription factors. Nature 2005;435:944–947
    OpenUrlCrossRefPubMedWeb of Science
  26. 26.↵
    Duncan SA, Navas MA, Dufort D, Rossant J, Stoffel M: Regulation of a transcription factor network required for differentiation and metabolism. Science (New York, NY) 1998;281:692–695
    OpenUrlAbstract/FREE Full Text
  27. 27.↵
    Navas MA, Vaisse C, Boger S, Heimesaat M, Kollee LA, Stoffel M: The human HNF-3 genes: cloning, partial sequence and mutation screening in patients with impaired glucose homeostasis. Human heredity 2000;50:370–381
    OpenUrlCrossRefPubMed
  28. 28.↵
    Togayachi A, Kozono Y, Kuno A, Ohkura T, Sato T, Hirabayashi J, Ikehara Y, Narimatsu H: Beta3GnT2 (B3GNT2), a major polylactosamine synthase: analysis of B3GNT2-deficient mice. Methods in enzymology 2010;479:185–204
    OpenUrlCrossRefPubMedWeb of Science
  29. 29.↵
    Comuzzie AG, Cole SA, Laston SL, Voruganti VS, Haack K, Gibbs RA, Butte NF: Novel genetic loci identified for the pathophysiology of childhood obesity in the Hispanic population. PLoS One 2012;7:e51954
    OpenUrlCrossRefPubMed
  30. 30.↵
    Brown FC, Schindelhaim CH, Pfeffer SR: GCC185 plays independent roles in Golgi structure maintenance and AP-1-mediated vesicle tethering. The Journal of cell biology 2011;194:779–787
    OpenUrlAbstract/FREE Full Text
  31. 31.↵
    Plotnikov D, Williams C, Guggenheim JA: Association between birth weight and refractive error in adulthood: a Mendelian randomisation study. The British journal of ophthalmology 2020;104:214–219
    OpenUrlAbstract/FREE Full Text
  32. 32.↵
    Mueller B, Klemm EJ, Spooner E, Claessen JH, Ploegh HL: SEL1L nucleates a protein complex required for dislocation of misfolded glycoproteins. Proc Natl Acad Sci U S A 2008;105:12325–12330
    OpenUrlAbstract/FREE Full Text
  33. 33.↵
    Chalasani N, Guo X, Loomba R, Goodarzi MO, Haritunians T, Kwon S, Cui J, Taylor KD, Wilson L, Cummings OW, Chen YD, Rotter JI: Genome-wide association study identifies variants associated with histologic features of nonalcoholic Fatty liver disease. Gastroenterology 2010;139:1567–1576, 1576.e1561-1566
    OpenUrlCrossRefPubMed
  34. 34.↵
    Stone S, Abkevich V, Russell DL, Riley R, Timms K, Tran T, Trem D, Frank D, Jammulapati S, Neff CD, Iliev D, Gress R, He G, Frech GC, Adams TD, Skolnick MH, Lanchbury JS, Gutin A, Hunt SC, Shattuck D: TBC1D1 is a candidate for a severe obesity gene and evidence for a gene/gene interaction in obesity predisposition. Hum Mol Genet 2006;15:2709–2720
    OpenUrlCrossRefPubMedWeb of Science
  35. 35.↵
    Rhee EP, Ho JE, Chen MH, Shen D, Cheng S, Larson MG, Ghorbani A, Shi X, Helenius IT, O’Donnell CJ, Souza AL, Deik A, Pierce KA, Bullock K, Walford GA, Vasan RS, Florez JC, Clish C, Yeh JR, Wang TJ, Gerszten RE: A genome-wide association study of the human metabolome in a community-based cohort. Cell metabolism 2013;18:130–143
    OpenUrl
  36. 36.↵
    Astle WJ, Elding H, Jiang T, Allen D, Ruklisa D, Mann AL, Mead D, Bouman H, Riveros-Mckay F, Kostadima MA, Lambourne JJ, Sivapalaratnam S, Downes K, Kundu K, Bomba L, Berentsen K, Bradley JR, Daugherty LC, Delaneau O, Freson K, Garner SF, Grassi L, Guerrero J, Haimel M, Janssen-Megens EM, Kaan A, Kamat M, Kim B, Mandoli A, Marchini J, Martens JHA, Meacham S, Megy K, O’Connell J, Petersen R, Sharifi N, Sheard SM, Staley JR, Tuna S, van der Ent M, Walter K, Wang SY, Wheeler E, Wilder SP, Iotchkova V, Moore C, Sambrook J, Stunnenberg HG, Di Angelantonio E, Kaptoge S, Kuijpers TW, Carrillo-de-Santa-Pau E, Juan D, Rico D, Valencia A, Chen L, Ge B, Vasquez L, Kwan T, Garrido-Martín D, Watt S, Yang Y, Guigo R, Beck S, Paul DS, Pastinen T, Bujold D, Bourque G, Frontini M, Danesh J, Roberts DJ, Ouwehand WH, Butterworth AS, Soranzo N: The Allelic Landscape of Human Blood Cell Trait Variation and Links to Common Complex Disease. Cell 2016;167:1415–1429.e1419
    OpenUrlCrossRefPubMed
  37. 37.↵
    Vella A, Cooper JD, Lowe CE, Walker N, Nutland S, Widmer B, Jones R, Ring SM, McArdle W, Pembrey ME, Strachan DP, Dunger DB, Twells RC, Clayton DG, Todd JA: Localization of a type 1 diabetes locus in the IL2RA/CD25 region by use of tag single-nucleotide polymorphisms. Am J Hum Genet 2005;76:773–779
    OpenUrlCrossRefPubMedWeb of Science
  38. 38.
    Qu H-Q, Montpetit A, Ge B, Hudson TJ, Polychronakos C: Toward Further Mapping of the Association Between the IL2RA Locus and Type 1 Diabetes. Diabetes 2007;56:1174–1176
    OpenUrlAbstract/FREE Full Text
  39. 39.↵
    Plagnol V, Howson JM, Smyth DJ, Walker N, Hafler JP, Wallace C, Stevens H, Jackson L, Simmonds MJ, Bingley PJ, Gough SC, Todd JA: Genome-wide association analysis of autoantibody positivity in type 1 diabetes cases. PLoS Genet 2011;7:e1002216
    OpenUrlCrossRefPubMedWeb of Science
  40. 40.↵
    Pulit SL, Stoneman C, Morris AP, Wood AR, Glastonbury CA, Tyrrell J, Yengo L, Ferreira T, Marouli E, Ji Y, Yang J, Jones S, Beaumont R, Croteau-Chonka DC, Winkler TW, Hattersley AT, Loos RJF, Hirschhorn JN, Visscher PM, Frayling TM, Yaghootkar H, Lindgren CM: Meta-analysis of genome-wide association studies for body fat distribution in 694 649 individuals of European ancestry. Hum Mol Genet 2019;28:166–174
    OpenUrl
  41. 41.↵
    Jaleco AC, Neves H, Hooijberg E, Gameiro P, Clode N, Haury M, Henrique D, Parreira L: Differential effects of Notch ligands Delta-1 and Jagged-1 in human lymphoid differentiation. J Exp Med 2001;194:991-1002
  42. 42.↵
    Li D, Kang Q, Wang DM: Constitutive coactivator of peroxisome proliferator-activated receptor (PPARgamma), a novel coactivator of PPARgamma that promotes adipogenesis. Molecular endocrinology (Baltimore, Md) 2007;21:2320–2333
    OpenUrlCrossRefPubMed
  43. 43.↵
    Lauc G, Huffman JE, Pučić M, Zgaga L, Adamczyk B, Mužinić A, Novokmet M, Polašek O, Gornik O, Krištić J, Keser T, Vitart V, Scheijen B, Uh HW, Molokhia M, Patrick AL, McKeigue P, Kolčić I, Lukić IK, Swann O, van Leeuwen FN, Ruhaak LR, Houwing-Duistermaat JJ, Slagboom PE, Beekman M, de Craen AJ, Deelder AM, Zeng Q, Wang W, Hastie ND, Gyllensten U, Wilson JF, Wuhrer M, Wright AF, Rudd PM, Hayward C, Aulchenko Y, Campbell H, Rudan I: Loci associated with N-glycosylation of human immunoglobulin G show pleiotropy with autoimmune diseases and haematological cancers. PLoS Genet 2013;9:e1003225
    OpenUrlCrossRefPubMed
  44. 44.↵
    Van den Steen P, Rudd PM, Dwek RA, Van Damme J, Opdenakker G: Cytokine and protease glycosylation as a regulatory mechanism in inflammation and autoimmunity. In Glycoimmunology 2, Springer, 1998, p. 133–143
    OpenUrl
  45. 45.↵
    Nakagawa H, Hato M, Takegawa Y, Deguchi K, Ito H, Takahata M, Iwasaki N, Minami A, Nishimura S: Detection of altered N-glycan profiles in whole serum from rheumatoid arthritis patients. Journal of chromatography B, Analytical technologies in the biomedical and life sciences 2007;853:133–137
    OpenUrlCrossRefPubMedWeb of Science
  46. 46.↵
    Bermingham ML, Colombo M, McGurnaghan SJ, Blackbourn LAK, Vučković F, Pučić Baković M, Trbojević- Akmačić I, Lauc G, Agakov F, Agakova AS, Hayward C, Klarić L, Palmer CNA, Petrie JR, Chalmers J, Collier A, Green F, Lindsay RS, Macrury S, McKnight JA, Patrick AW, Thekkepat S, Gornik O, McKeigue PM, Colhoun HM: N-Glycan Profile and Kidney Disease in Type 1 Diabetes. Diabetes Care 2018;41:79–87
    OpenUrlAbstract/FREE Full Text
  47. 47.↵
    Trbojević Akmačić I, Ventham NT, Theodoratou E, Vučković F, Kennedy NA, Krištić J, Nimmo ER, Kalla R, Drummond H, Štambuk J, Dunlop MG, Novokmet M, Aulchenko Y, Gornik O, Campbell H, Pučić Baković M, Satsangi J, Lauc G: Inflammatory bowel disease associates with proinflammatory potential of the immunoglobulin G glycome. Inflammatory bowel diseases 2015;21:1237–1247
    OpenUrlCrossRefPubMed
  48. 48.↵
    Boachie AM, Degraff D, Yu X, Sun Q, Friedman D, Gronostajski R, Matusik R: Abstract 1231: Nuclear Factor I family members interact with FoxA1 to regulate androgen responsive promoters. Cancer Research 2010;70:1231–1231
    OpenUrlCrossRef
  49. 49.↵
    Zhao E, Li Y, Fu X, Zeng L, Zeng H, Jin W, Chen J, Yin G, Qian J, Ying K, Xie Y, Zhao RC, Mao Y: Cloning and characterization of human synaptotagmin 10 gene. DNA sequence : the journal of DNA sequencing and mapping 2003;14:393–398
    OpenUrl
  50. 50.↵
    Burda P, Aebi M: The ALG10 locus of Saccharomyces cerevisiae encodes the alpha-1,2 glucosyltransferase of the endoplasmic reticulum: the terminal glucose of the lipid-linked oligosaccharide is required for efficient N-linked glycosylation. Glycobiology 1998;8:455–462
    OpenUrlCrossRefPubMedWeb of Science
  51. 51.↵
    Tachmazidou I, Süveges D, Min JL, Ritchie GRS, Steinberg J, Walter K, Iotchkova V, Schwartzentruber J, Huang J, Memari Y, McCarthy S, Crawford AA, Bombieri C, Cocca M, Farmaki AE, Gaunt TR, Jousilahti P, Kooijman MN, Lehne B, Malerba G, Männistö S, Matchan A, Medina-Gomez C, Metrustry SJ, Nag A, Ntalla I, Paternoster L, Rayner NW, Sala C, Scott WR, Shihab HA, Southam L, St Pourcain B, Traglia M, Trajanoska K, Zaza G, Zhang W, Artigas MS, Bansal N, Benn M, Chen Z, Danecek P, Lin WY, Locke A, Luan J, Manning AK, Mulas A, Sidore C, Tybjaerg-Hansen A, Varbo A, Zoledziewska M, Finan C, Hatzikotoulas K, Hendricks AE, Kemp JP, Moayyeri A, Panoutsopoulou K, Szpak M, Wilson SG, Boehnke M, Cucca F, Di Angelantonio E, Langenberg C, Lindgren C, McCarthy MI, Morris AP, Nordestgaard BG, Scott RA, Tobin MD, Wareham NJ, Burton P, Chambers JC, Smith GD, Dedoussis G, Felix JF, Franco OH, Gambaro G, Gasparini P, Hammond CJ, Hofman A, Jaddoe VWV, Kleber M, Kooner JS, Perola M, Relton C, Ring SM, Rivadeneira F, Salomaa V, Spector TD, Stegle O, Toniolo D, Uitterlinden AG, Barroso I, Greenwood CMT, Perry JRB, Walker BR, Butterworth AS, Xue Y, Durbin R, Small KS, Soranzo N, Timpson NJ, Zeggini E: Whole-Genome Sequencing Coupled to Imputation Discovers Genetic Signals for Anthropometric Traits. Am J Hum Genet 2017;100:865–884
    OpenUrlCrossRef
  52. 52.↵
    Kichaev G, Bhatia G, Loh PR, Gazal S, Burch K, Freund MK, Schoech A, Pasaniuc B, Price AL: Leveraging Polygenic Functional Enrichment to Improve GWAS Power. Am J Hum Genet 2019;104:65–75
    OpenUrlCrossRefPubMed
  53. 53.↵
    Kinoshita T, Doi K, Sugiyama H, Kinoshita S, Wada M, Naruto S, Tomonaga A: Knowledge-Based Identification of the ERK2/STAT3 Signal Pathway as a Therapeutic Target for Type 2 Diabetes and Drug Discovery. Chemical Biology & Drug Design 2011;78:471–476
    OpenUrl
  54. 54.↵
    Rudman N, Gornik O, Lauc G: Altered N-glycosylation profiles as potential biomarkers and drug targets in diabetes. FEBS Letters 2019;593:1598–1615
    OpenUrlCrossRefPubMed
  55. 55.
    Wilkin TJ: The accelerator hypothesis: weight gain as the missing link between Type I and Type II diabetes. Diabetologia 2001;44:914–922
    OpenUrlCrossRefPubMedWeb of Science
  56. 56.
    Kibirige M, Metcalf B, Renuka R, Wilkin T: Testing the accelerator hypothesis: the relationship between body mass and age at diagnosis of type 1 diabetes. Diabetes care 2003;26:2865–2870
    OpenUrlAbstract/FREE Full Text
  57. 57.↵
    Volta F, Gerdes JM: The role of primary cilia in obesity and diabetes. Annals of the New York Academy of Sciences 2017;1391:71–84
    OpenUrlCrossRef
Back to top
Posted October 14, 2020.
Download PDF

Supplementary Material

Data/Code
Email
Genetics of Low Polygenic Risk Score Type 1 Diabetes Patients: rare variants in 22 novel loci
Jingchun Qu, Hui-Qi Qu, Jonathan Bradfield, Joseph Glessner, Xiao Chang, Lifeng Tian, Michael March, Jeffrey D Roizen, Patrick Sleiman, Hakon Hakonarson
medRxiv 2020.10.13.20211987; doi: https://doi.org/10.1101/2020.10.13.20211987
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area