Genome-wide analysis in over 1 million individuals of European ancestry yields improved polygenic risk scores for blood pressure traits.

Hypertension affects more than one billion people worldwide. Here we identify 113 novel loci, reporting a total of 2,103 independent genetic signals (P < 5 × 10-8) from the largest single-stage blood pressure (BP) genome-wide association study to date (n = 1,028,980 European individuals). These associations explain more than 60% of single nucleotide polymorphism-based BP heritability. Comparing top versus bottom deciles of polygenic risk scores (PRSs) reveals clinically meaningful differences in BP (16.9 mmHg systolic BP, 95% CI, 15.5-18.2 mmHg, P = 2.22 × 10-126) and more than a sevenfold higher odds of hypertension risk (odds ratio, 7.33; 95% CI, 5.54-9.70; P = 4.13 × 10-44) in an independent dataset. Adding PRS into hypertension-prediction models increased the area under the receiver operating characteristic curve (AUROC) from 0.791 (95% CI, 0.781-0.801) to 0.826 (95% CI, 0.817-0.836, ∆AUROC, 0.035, P = 1.98 × 10-34). We compare the 2,103 loci results in non-European ancestries and show significant PRS associations in a large African-American sample. Secondary analyses implicate 500 genes previously unreported for BP. Our study highlights the role of increasingly large genomic studies for precision health research.

Genetic imputation of kidney transcriptome, proteome and multi-omics illuminates new blood pressure and hypertension targets.

Genetic mechanisms of blood pressure (BP) regulation remain poorly defined. Using kidney-specific epigenomic annotations and 3D genome information we generated and validated gene expression prediction models for the purpose of transcriptome-wide association studies in 700 human kidneys. We identified 889 kidney genes associated with BP of which 399 were prioritised as contributors to BP regulation. Imputation of kidney proteome and microRNAome uncovered 97 renal proteins and 11 miRNAs associated with BP. Integration with plasma proteomics and metabolomics illuminated circulating levels of myo-inositol, 4-guanidinobutanoate and angiotensinogen as downstream effectors of several kidney BP genes (SLC5A11, AGMAT, AGT, respectively). We showed that genetically determined reduction in renal expression may mimic the effects of rare loss-of-function variants on kidney mRNA/protein and lead to an increase in BP (e.g., ENPEP). We demonstrated a strong correlation (r = 0.81) in expression of protein-coding genes between cells harvested from urine and the kidney highlighting a diagnostic potential of urinary cell transcriptomics. We uncovered adenylyl cyclase activators as a repurposing opportunity for hypertension and illustrated examples of BP-elevating effects of anticancer drugs (e.g. tubulin polymerisation inhibitors). Collectively, our studies provide new biological insights into genetic regulation of BP with potential to drive clinical translation in hypertension.

Genomic insights into the comorbidity between type 2 diabetes and schizophrenia.

Multimorbidity represents an increasingly important public health challenge with far-reaching implications for health management and policy. Mental health and metabolic diseases have a well-established epidemiological association. In this study, we investigate the genetic intersection between type 2 diabetes and schizophrenia. We use Mendelian randomization to examine potential causal relationships between the two conditions and related endophenotypes. We report no compelling evidence that type 2 diabetes genetic liability potentially causally influences schizophrenia risk and vice versa. Our findings show that increased body mass index (BMI) has a protective effect against schizophrenia, in contrast to the well-known risk-increasing effect of BMI on type 2 diabetes risk. We identify evidence of colocalization of association signals for these two conditions at 11 genomic loci, six of which have opposing directions of effect for type 2 diabetes and schizophrenia. To elucidate these colocalizing signals, we integrate multi-omics data from bulk and single-cell gene expression studies, along with functional information. We identify putative effector genes and find that they are enriched for homeostasis and lipid-related pathways. We also highlight drug repurposing opportunities including N-methyl-D-aspartate (NMDA) receptor antagonists. Our findings provide insights into shared biological mechanisms for type 2 diabetes and schizophrenia, highlighting common factors that influence the risk of the two conditions in opposite directions and shedding light on the complex nature of this comorbidity.

Psoriatic arthritis: the role of self-reported non-adherence, non-trough drug levels, immunogenicity and conventional synthetic DMARD co-therapy in adalimumab and etanercept response.

Objective: The aim of this study was to assess the relationship between self-reported non-adherence, non-trough drug levels, immunogenicity and conventional synthetic DMARD (csDMARD) co-therapy in TNF inhibitor (TNF-i) drug response in PsA. Methods: Serum samples and adherence questionnaires were collected at baseline, 3, 6 and 12 months for PsA patients prescribed TNF-i. Non-trough adalimumab (ADL) and etanercept (ETN) drug levels were measured at 3 and 6 months using commercially available ELISAs. Clinical response was assessed using PsA response criteria (PsARC) and change in 28-joint DAS (ΔDAS28) between baseline and 3, 6 and 12 months. Results: In 244 PsA patients (52.5% ADL and 47.5% ETN), self-reported non-adherence was associated with PsARC non-response over 12 months using generalized estimating equation (GEE) modelling (P = 0.037). However, there was no significant difference between non-trough ADL or ETN drug levels based on self-reported non-adherence. Higher ETN levels at 3 months were associated with PsARC response at 3 (P = 0.015), 6 (P = 0.037) and 12 months (P = 0.015) and over 12 months using GEE modelling (P = 0.026). Increased ADL drug levels at 3 months were associated with greater ΔDAS28 at 3 months (P = 0.019). ADL anti-drug antibody-positive status was significantly associated with lower 3- and 6-month ADL levels (P < 0.001) and ΔDAS28 and PsARC response at 3, 6 and 12 months. Meanwhile, MTX co-therapy was associated with a reduction in immunogenicity at 3 and 6 months (P = 0.008 and P = 0.024). Conclusion: Although both were associated with reduced response, the objectively measured non-trough drug levels showed more significant associations with drug response than self-reported non-adherence measures.

Genetic Association and Transferability for Urinary Albumin-Creatinine Ratio as a Marker of Kidney Disease in four Sub-Saharan African Populations and non-continental Individuals of African Ancestry.

BACKGROUND: Genome-wide association studies (GWAS) have predominantly focused on populations of European and Asian ancestry, limiting our understanding of genetic factors influencing kidney disease in Sub-Saharan African (SSA) populations. This study presents the largest GWAS for urinary albumin-to-creatinine ratio (UACR) in SSA individuals, including 8,970 participants living in different African regions and an additional 9,705 non-resident individuals of African ancestry from the UK Biobank and African American cohorts. METHODS: Urine biomarkers and genotype data were obtained from two SSA cohorts (AWI-Gen and ARK), and two non-resident African-ancestry studies (UK Biobank and CKD-Gen Consortium). Association testing and meta-analyses were conducted, with subsequent fine-mapping, conditional analyses, and replication studies. Polygenic scores (PGS) were assessed for transferability across populations. RESULTS: Two genome-wide significant (P<5x10-8) UACR-associated loci were identified, one in the BMP6 region on chromosome 6, in the meta-analysis of resident African individuals, and another in the HBB region on chromosome 11 in the meta-analysis of non-resident SSA individuals, as well as the combined meta-analysis of all studies. Replication of previous significant results confirmed associations in known UACR-associated regions, including THB53, GATM, and ARL15. PGS estimated using previous studies from European ancestry, African ancestry, and multi-ancestry cohorts exhibited limited transferability of PGS across populations, with less than 1% of observed variance explained. CONCLUSION: This study contributes novel insights into the genetic architecture of kidney disease in SSA populations, emphasizing the need for conducting genetic research in diverse cohorts. The identified loci provide a foundation for future investigations into the genetic susceptibility to chronic kidney disease in underrepresented African populations Additionally, there is a need to develop integrated scores using multi-omics data and risk factors specific to the African context to improve the accuracy of predicting disease outcomes. METHODS: Urine biomarkers and genotype data were obtained from two SSA cohorts (AWI-Gen and ARK), and two non-resident African-ancestry studies (UK Biobank and CKD-Gen Consortium). Association testing and meta-analyses were conducted, with subsequent fine-mapping, conditional analyses, and replication studies. Polygenic scores (PGS) were assessed for transferability across populations. RESULTS: Two genome-wide significant (P<5x10-8) UACR-associated loci were identified, one in the BMP6 region on chromosome 6, in the meta-analysis of resident African individuals, and another in the HBB region on chromosome 11 in the meta-analysis of non-resident SSA individuals, as well as the combined meta-analysis of all studies. Replication of previous significant results confirmed associations in known UACR-associated regions, including THB53, GATM, and ARL15. PGS estimated using previous studies from European ancestry, African ancestry, and multi-ancestry cohorts exhibited limited transferability of PGS across populations, with less than 1% of observed variance explained. CONCLUSION: This study contributes novel insights into the genetic architecture of kidney function in SSA populations, emphasizing the need for conducting genetic research in diverse cohorts. The identified loci provide a foundation for future investigations into the genetic susceptibility to chronic kidney disease in underrepresented African populations.

Advancing equity in human genomics through tissue-specific multi-ancestry molecular data.

There is a pressing need to generate molecular data from diverse tissues across global populations. These currently missing data are necessary to resolve genome-wide association study loci, identify effector genes, and move the translational genomics needle beyond European-ancestry individuals and the minority of diseases for which blood is the relevant tissue.

Multi-ancestry meta-analysis of host genetic susceptibility to tuberculosis identifies shared genetic architecture.

The heritability of susceptibility to tuberculosis (TB) disease has been well recognized. Over 100 genes have been studied as candidates for TB susceptibility, and several variants were identified by genome-wide association studies (GWAS), but few replicate. We established the International Tuberculosis Host Genetics Consortium to perform a multi-ancestry meta-analysis of GWAS, including 14,153 cases and 19,536 controls of African, Asian, and European ancestry. Our analyses demonstrate a substantial degree of heritability (pooled polygenic h2 = 26.3%, 95% CI 23.7-29.0%) for susceptibility to TB that is shared across ancestries, highlighting an important host genetic influence on disease. We identified one global host genetic correlate for TB at genome-wide significance (p<5 × 10-8) in the human leukocyte antigen (HLA)-II region (rs28383206, p-value=5.2 × 10-9) but failed to replicate variants previously associated with TB susceptibility. These data demonstrate the complex shared genetic architecture of susceptibility to TB and the importance of large-scale GWAS analysis across multiple ancestries experiencing different levels of infection pressure.

Non-trough serum drug levels of adalimumab and etanercept are associated with response in patients with psoriatic arthritis.

OBJECTIVES: Up to 40% of psoriatic arthritis (PsA) patients experience first-line Tumour Necrosis Factor inhibitors (TNF-i) failure. Lower serum drug levels (SDL) have been associated with lower response in autoimmune conditions. This study aimed to: (i) establish the relationship between adalimumab (ADL) and etanercept (ETN) SDL and 3-month response; and (ii) identify optimal non-trough SDL thresholds in PsA. METHODS: PsA patients commencing ADL or ETN were recruited to the UK observational study OUTPASS. Patients were seen pre-TNF-i and at 3 months when response was measured, and non-trough serum samples collected. Response was defined according to the PsARC or EULAR criteria. Descriptive statistics and concentration-effect curves established differences in SDL based on response. Receiver operating characteristics and regression identified optimal SDL thresholds. RESULTS: PsA ETN (n = 97) PsARC and EULAR good responders had significantly higher 3-month SDL compared with non-responders (p= 0.006 and p= 0.020 respectively). Non-trough 3-month ETN SDL discriminated PsARC responders from non-responders (AUC = 0.70), with a threshold of 1.8 µg/ml being 63% specific and 69% sensitive. EULAR good and non-/moderate responders were discriminated with an AUC of 0.65 with a threshold of 2.0 µg/ml being 57% specific and 69% sensitive. ADL prescribed (n = 104) EULAR good responders had significantly higher 3-month SDL (p= 0.049). Non-trough 3-month ADL SDL discriminated EULAR good and non-/moderate responders (AUC = 0.63) with a threshold of 3.6 µg/ml being 48% specific and 81% sensitive. CONCLUSION: Higher 3-month SDL were detected in responders. Interventions to optimise SDL may improve treatment response earlier. This study suggests 3-month SDL thresholds which may be useful in clinical practice to optimise treatment response.

Leveraging information between multiple population groups and traits improves fine-mapping resolution.

Statistical fine-mapping helps to pinpoint likely causal variants underlying genetic association signals. Its resolution can be improved by (i) leveraging information between traits; and (ii) exploiting differences in linkage disequilibrium structure between diverse population groups. Using association summary statistics, MGflashfm jointly fine-maps signals from multiple traits and population groups; MGfm uses an analogous framework to analyse each trait separately. We also provide a practical approach to fine-mapping with out-of-sample reference panels. In simulation studies we show that MGflashfm and MGfm are well-calibrated and that the mean proportion of causal variants with PP > 0.80 is above 0.75 (MGflashfm) and 0.70 (MGfm). In our analysis of four lipids traits across five population groups, MGflashfm gives a median 99% credible set reduction of 10.5% over MGfm. MGflashfm and MGfm only require summary level data, making them very useful fine-mapping tools in consortia efforts where individual-level data cannot be shared.

Genomic insights into the comorbidity between type 2 diabetes and schizophrenia.

Multimorbidity represents an increasingly important public health challenge with far-reaching implications for health management and policy. Mental health and metabolic diseases have a well-established epidemiological association. In this study, we investigate the genetic intersection between type 2 diabetes and schizophrenia. We use Mendelian randomization to examine potential causal relationships between the two conditions and related endophenotypes. We report no compelling evidence that type 2 diabetes genetic liability potentially causally influences schizophrenia risk and vice versa. Our findings show that increased body mass index (BMI) has a protective effect against schizophrenia, in contrast to the well-known risk-increasing effect of BMI on type 2 diabetes risk. We identify evidence of colocalization of association signals for these two conditions at 11 genomic loci, six of which have opposing directions of effect for type 2 diabetes and schizophrenia. To elucidate these colocalizing signals, we integrate multi-omics data from bulk and single-cell gene expression studies, along with functional information. We identify high-confidence effector genes and find that they are enriched for homeostasis and lipid-related pathways. We also highlight drug repurposing opportunities including N-methyl-D-aspartate (NMDA) receptor antagonists. Our findings provide insights into shared biological mechanisms for type 2 diabetes and schizophrenia, highlighting common factors that influence the risk of the two conditions in opposite directions and shedding light on the complex nature of this comorbidity.

Multi-trait discovery and fine-mapping of lipid loci in 125,000 individuals of African ancestry.

Most genome-wide association studies (GWAS) for lipid traits focus on the separate analysis of lipid traits. Moreover, there are limited GWASs evaluating the genetic variants associated with multiple lipid traits in African ancestry. To further identify and localize loci with pleiotropic effects on lipid traits, we conducted a genome-wide meta-analysis, multi-trait analysis of GWAS (MTAG), and multi-trait fine-mapping (flashfm) in 125,000 individuals of African ancestry. Our meta-analysis and MTAG identified four and 14 novel loci associated with lipid traits, respectively. flashfm yielded an 18% mean reduction in the 99% credible set size compared to single-trait fine-mapping with JAM. Moreover, we identified more genetic variants with a posterior probability of causality >0.9 with flashfm than with JAM. In conclusion, we identified additional novel loci associated with lipid traits, and flashfm reduced the 99% credible set size to identify causal genetic variants associated with multiple lipid traits in African ancestry.

Integration of genetic fine-mapping and multi-omics data reveals candidate effector genes for hypertension.

Genome-wide association studies of blood pressure (BP) have identified >1,000 loci, but the effector genes and biological pathways at these loci are mostly unknown. Using published association summary statistics, we conducted annotation-informed fine-mapping incorporating tissue-specific chromatin segmentation and colocalization to identify causal variants and candidate effector genes for systolic BP, diastolic BP, and pulse pressure. We observed 532 distinct signals associated with ≥2 BP traits and 84 with all three. For >20% of signals, a single variant accounted for >75% posterior probability, 65 were missense variants in known (SLC39A8, ADRB2, and DBH) and previously unreported BP candidate genes (NRIP1 and MMP14). In disease-relevant tissues, we colocalized >80 and >400 distinct signals for each BP trait with cis-eQTLs and regulatory regions from promoter capture Hi-C, respectively. Integrating mouse, human disorder, gene expression and tissue abundance data, and literature review, we provide consolidated evidence for 436 BP candidate genes for future functional validation and discover several potential drug targets.

Integrated proteomics and genomics analysis of paradoxical eczema in psoriasis patients treated with biologics.

BACKGROUND: Few studies have explored the immunology and genetic risk of paradoxical eczema occurring as an adverse event of biologic therapy in patients with psoriasis. OBJECTIVES: We sought to describe the systemic inflammatory signature of paradoxical eczema using proteomics and explore whether this is genetically mediated. METHODS: This study used the Olink Target 96 Inflammation panel on 256 serum samples from 71 patients with psoriasis and paradoxical eczema, and 75 controls with psoriasis to identify differentially expressed proteins and enriched gene sets. Case samples from 1 or more time points (T1 prebiologic, T2 postbiologic, and T3 postparadoxical eczema) were matched 1:1 with control samples. Genes contributing to enriched gene sets were selected, and functional single nucleotide polymorphisms used to create polygenic risk scores in a genotyped cohort of 88 paradoxical eczema cases and 3124 psoriasis controls. RESULTS: STAMBP expression was lower in cases at T1 than in controls (log-fold change: -0.44; adjusted P = .022); no other proteins reached statistical significance at equivalent time points. Eleven gene sets including cytokine and chemokine pathways were enriched in cases at T2 and 10 at T3. Of the 39 proteins contributing to enriched gene sets, the majority are associated with the atopic dermatitis serum proteome. A polygenic risk score including 38 functional single nucleotide polymorphisms linked to enriched gene sets was associated with paradoxical eczema (adjusted P = .046). CONCLUSIONS: The paradoxical eczema systemic inflammatory proteome trends toward atopic dermatitis at a gene-set level and is detectable before onset of the phenotype. This signature could be genetically determined.

Genome-wide association analysis of cystatin-C kidney function in continental Africa.

BACKGROUND: Chronic kidney disease is becoming more prevalent in Africa, and its genetic determinants are poorly understood. Creatinine-based estimated glomerular filtration rate (eGFR) is commonly used to estimate kidney function, modelling the excretion of the endogenous biomarker (creatinine). However, eGFR based on creatinine has been shown to inadequately detect individuals with low kidney function in Sub-Saharan Africa, with eGFR based on cystatin-C (eGFRcys) exhibiting significantly superior performance. Therefore, we opted to conduct a GWAS for eGFRcys. METHODS: Using the Uganda Genomic Resource, we performed a genome-wide association study (GWAS) of eGFRcys in 5877 Ugandans and evaluated replication in independent studies. Subsequently, putative causal variants were screened through Bayesian fine-mapping. Functional annotation of the GWAS loci was performed using Functional Mapping and Annotation (FUMA). FINDINGS: Three independent lead single nucleotide polymorphisms (SNPs) (P-value <5 × 10-8 (based on likelihood ratio test (LRT))) were identified; rs59288815 (ANK3), rs4277141 (OR51B5) and rs911119 (CST3). From fine-mapping, rs59288815 and rs911119 each had a posterior probability of causality of >99%. The rs911119 SNP maps to the cystatin C gene and has been previously associated with eGFRcys among Europeans. With gene-set enrichment analyses of the olfactory receptor family 51 overlapping genes, we identified an association with the G-alpha-S signalling events. INTERPRETATION: Our study found two previously unreported associated SNPs for eGFRcys in continental Africans (rs59288815 and rs4277141) and validated a previously well-established SNP (rs911119) for eGFRcys. The identified gene-set enrichment for the G-protein signalling pathways relates to the capacity of the kidney to readily adapt to an ever-changing environment. Additional GWASs are required to represent the diverse regions in Africa. FUNDING: Wellcome (220740/Z/20/Z).

Genetic underpinning of the comorbidity between type 2 diabetes and osteoarthritis.

Multimorbidity is a rising public health challenge with important implications for health management and policy. The most common multimorbidity pattern is the combination of cardiometabolic and osteoarticular diseases. Here, we study the genetic underpinning of the comorbidity between type 2 diabetes and osteoarthritis. We find genome-wide genetic correlation between the two diseases and robust evidence for association-signal colocalization at 18 genomic regions. We integrate multi-omics and functional information to resolve the colocalizing signals and identify high-confidence effector genes, including FTO and IRX3, which provide proof-of-concept insights into the epidemiologic link between obesity and both diseases. We find enrichment for lipid metabolism and skeletal formation pathways for signals underpinning the knee and hip osteoarthritis comorbidities with type 2 diabetes, respectively. Causal inference analysis identifies complex effects of tissue-specific gene expression on comorbidity outcomes. Our findings provide insights into the biological basis for the type 2 diabetes-osteoarthritis disease co-occurrence.

Bridging the diversity gap: Analytical and study design considerations for improving the accuracy of trans-ancestry genetic prediction.

Genetic prediction of common complex disease risk is an essential component of precision medicine. Currently, genome-wide association studies (GWASs) are mostly composed of European-ancestry samples and resulting polygenic scores (PGSs) have been shown to poorly transfer to other ancestries partly due to heterogeneity of allelic effects between populations. Fixed-effects (FETA) and random-effects (RETA) trans-ancestry meta-analyses do not model such ancestry-related heterogeneity, while ancestry-specific (AS) scores may suffer from low power due to low sample sizes. In contrast, trans-ancestry meta-regression (TAMR) builds ancestry-aware PGS that account for more complex trans-ancestry architectures. Here, we examine the predictive performance of these four PGSs under multiple genetic architectures and ancestry configurations. We show that the predictive performance of FETA and RETA is strongly affected by cross-ancestry genetic heterogeneity, while AS PGS performance decreases in under-represented target populations. TAMR PGS is also impacted by heterogeneity but maintains good prediction performance in most situations, especially in ancestry-diverse scenarios. In simulations of human complex traits, TAMR scores currently explain 25% more phenotypic variance than AS in triglyceride levels and 33% more phenotypic variance than FETA in type 2 diabetes in most non-European populations. Importantly, a high proportion of non-European-ancestry individuals is needed to reach prediction levels that are comparable in those populations to the one observed in European-ancestry studies. Our results highlight the need to rebalance the ancestral composition of GWAS to enable accurate prediction in non-European-ancestry groups, and demonstrate the relevance of meta-regression approaches for compensating some of the current population biases in GWAS.

Atopic Polygenic Risk Score Is Associated with Paradoxical Eczema Developing in Patients with Psoriasis Treated with Biologics.

Biologic therapies for psoriasis can cause paradoxical eczema. The role of genetic factors in its pathogenesis is unknown. To identify risk variants, we conducted a GWAS of 3,212 patients with psoriasis, of whom 88 developed paradoxical eczema. Two lead SNPs reached genome-wide significance (P ≤ 5 × 10-8) for association with paradoxical eczema: rs192705221 (near UNC5B, P = 9.52 × 10-10) and rs72925168 (within SLC1A2, P = 1.66 × 10-9). Genome-wide significant SNPs from published GWAS were used to generate polygenic risk scores (PRSs) for atopic eczema, general atopic disease, or a combination, which were tested for association with paradoxical eczema. Improvement over a clinical risk model was assessed by the area under the curve. All three atopy polygenic risk scores were associated with paradoxical eczema (P < 0.05); polygenic risk score for a combination of atopic eczema and general atopic disease had the strongest association (OR = 1.83, 95% CI = 1.17-2.84, P = 0.0078). Including atopic polygenic risk scores in the multivariable model, which included age, sex, atopic background, and psoriatic arthritis history, increased the area under the curve from 0.671 to 0.681-0.686. Atopic genetic burden is associated with paradoxical eczema occurring in biologic-treated patients with psoriasis, indicating shared underlying mechanisms. Incorporating genetic risk may improve treatment outcome prediction models for psoriasis.

Genome-Wide Gene-Set Analysis Identifies Molecular Mechanisms Associated with ALS.

Amyotrophic lateral sclerosis (ALS) is a fatal late-onset motor neuron disease characterized by the loss of the upper and lower motor neurons. Our understanding of the molecular basis of ALS pathology remains elusive, complicating the development of efficient treatment. Gene-set analyses of genome-wide data have offered insight into the biological processes and pathways of complex diseases and can suggest new hypotheses regarding causal mechanisms. Our aim in this study was to identify and explore biological pathways and other gene sets having genomic association to ALS. Two cohorts of genomic data from the dbGaP repository were combined: (a) the largest available ALS individual-level genotype dataset (N = 12,319), and (b) a similarly sized control cohort (N = 13,210). Following comprehensive quality control pipelines, imputation and meta-analysis, we assembled a large European descent ALS-control cohort of 9244 ALS cases and 12,795 healthy controls represented by genetic variants of 19,242 genes. Multi-marker analysis of genomic annotation (MAGMA) gene-set analysis was applied to an extensive collection of 31,454 gene sets from the molecular signatures database (MSigDB). Statistically significant associations were observed for gene sets related to immune response, apoptosis, lipid metabolism, neuron differentiation, muscle cell function, synaptic plasticity and development. We also report novel interactions between gene sets, suggestive of mechanistic overlaps. A manual meta-categorization and enrichment mapping approach is used to explore the overlap of gene membership between significant gene sets, revealing a number of shared mechanisms.

Loci for insulin processing and secretion provide insight into type 2 diabetes risk.

Insulin secretion is critical for glucose homeostasis, and increased levels of the precursor proinsulin relative to insulin indicate pancreatic islet beta-cell stress and insufficient insulin secretory capacity in the setting of insulin resistance. We conducted meta-analyses of genome-wide association results for fasting proinsulin from 16 European-ancestry studies in 45,861 individuals. We found 36 independent signals at 30 loci (p value < 5 × 10-8), which validated 12 previously reported loci for proinsulin and ten additional loci previously identified for another glycemic trait. Half of the alleles associated with higher proinsulin showed higher rather than lower effects on glucose levels, corresponding to different mechanisms. Proinsulin loci included genes that affect prohormone convertases, beta-cell dysfunction, vesicle trafficking, beta-cell transcriptional regulation, and lysosomes/autophagy processes. We colocalized 11 proinsulin signals with islet expression quantitative trait locus (eQTL) data, suggesting candidate genes, including ARSG, WIPI1, SLC7A14, and SIX3. The NKX6-3/ANK1 proinsulin signal colocalized with a T2D signal and an adipose ANK1 eQTL signal but not the islet NKX6-3 eQTL. Signals were enriched for islet enhancers, and we showed a plausible islet regulatory mechanism for the lead signal in the MADD locus. These results show how detailed genetic studies of an intermediate phenotype can elucidate mechanisms that may predispose one to disease.