Genetics of circulating inflammatory proteins identifies drivers of immune-mediated disease risk and therapeutic targets.

Circulating proteins have important functions in inflammation and a broad range of diseases. To identify genetic influences on inflammation-related proteins, we conducted a genome-wide protein quantitative trait locus (pQTL) study of 91 plasma proteins measured using the Olink Target platform in 14,824 participants. We identified 180 pQTLs (59 cis, 121 trans). Integration of pQTL data with eQTL and disease genome-wide association studies provided insight into pathogenesis, implicating lymphotoxin-α in multiple sclerosis. Using Mendelian randomization (MR) to assess causality in disease etiology, we identified both shared and distinct effects of specific proteins across immune-mediated diseases, including directionally discordant effects of CD40 on risk of rheumatoid arthritis versus multiple sclerosis and inflammatory bowel disease. MR implicated CXCL5 in the etiology of ulcerative colitis (UC) and we show elevated gut CXCL5 transcript expression in patients with UC. These results identify targets of existing drugs and provide a powerful resource to facilitate future drug target prioritization.

Leveraging cross-species transcription factor binding site patterns: from diabetes risk loci to disease mechanisms.

Genome-wide association studies have revealed numerous risk loci associated with diverse diseases. However, identification of disease-causing variants within association loci remains a major challenge. Divergence in gene expression due to cis-regulatory variants in noncoding regions is central to disease susceptibility. We show that integrative computational analysis of phylogenetic conservation with a complexity assessment of co-occurring transcription factor binding sites (TFBS) can identify cis-regulatory variants and elucidate their mechanistic role in disease. Analysis of established type 2 diabetes risk loci revealed a striking clustering of distinct homeobox TFBS. We identified the PRRX1 homeobox factor as a repressor of PPARG2 expression in adipose cells and demonstrate its adverse effect on lipid metabolism and systemic insulin sensitivity, dependent on the rs4684847 risk allele that triggers PRRX1 binding. Thus, cross-species conservation analysis at the level of co-occurring TFBS provides a valuable contribution to the translation of genetic association signals to disease-related molecular 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.

Differences and commonalities in the genetic architecture of protein quantitative trait loci in European and Arab populations.

Polygenic scores (PGS) can identify individuals at risk of adverse health events and guide genetics-based personalized medicine. However, it is not clear how well PGS translate between different populations, limiting their application to well-studied ethnicities. Proteins are intermediate traits linking genetic predisposition and environmental factors to disease, with numerous blood circulating protein levels representing functional readouts of disease-related processes. We hypothesized that studying the genetic architecture of a comprehensive set of blood-circulating proteins between a European and an Arab population could shed fresh light on the translatability of PGS to understudied populations. We therefore conducted a genome-wide association study with whole-genome sequencing data using 1301 proteins measured on the SOMAscan aptamer-based affinity proteomics platform in 2935 samples of Qatar Biobank and evaluated the replication of protein quantitative traits (pQTLs) from European studies in an Arab population. Then, we investigated the colocalization of shared pQTL signals between the two populations. Finally, we compared the performance of protein PGS derived from a Caucasian population in a European and an Arab cohort. We found that the majority of shared pQTL signals (81.8%) colocalized between both populations. About one-third of the genetic protein heritability was explained by protein PGS derived from a European cohort, with protein PGS performing ~20% better in Europeans when compared to Arabs. Our results are relevant for the translation of PGS to non-Caucasian populations, as well as for future efforts to extend genetic research to understudied populations.

Role of human plasma metabolites in prediabetes and type 2 diabetes from the IMI-DIRECT study.

AIMS/HYPOTHESIS: Type 2 diabetes is a chronic condition that is caused by hyperglycaemia. Our aim was to characterise the metabolomics to find their association with the glycaemic spectrum and find a causal relationship between metabolites and type 2 diabetes. METHODS: As part of the Innovative Medicines Initiative - Diabetes Research on Patient Stratification (IMI-DIRECT) consortium, 3000 plasma samples were measured with the Biocrates AbsoluteIDQ p150 Kit and Metabolon analytics. A total of 911 metabolites (132 targeted metabolomics, 779 untargeted metabolomics) passed the quality control. Multivariable linear and logistic regression analysis estimates were calculated from the concentration/peak areas of each metabolite as an explanatory variable and the glycaemic status as a dependent variable. This analysis was adjusted for age, sex, BMI, study centre in the basic model, and additionally for alcohol, smoking, BP, fasting HDL-cholesterol and fasting triacylglycerol in the full model. Statistical significance was Bonferroni corrected throughout. Beyond associations, we investigated the mediation effect and causal effects for which causal mediation test and two-sample Mendelian randomisation (2SMR) methods were used, respectively. RESULTS: In the targeted metabolomics, we observed four (15), 34 (99) and 50 (108) metabolites (number of metabolites observed in untargeted metabolomics appear in parentheses) that were significantly different when comparing normal glucose regulation vs impaired glucose regulation/prediabetes, normal glucose regulation vs type 2 diabetes, and impaired glucose regulation vs type 2 diabetes, respectively. Significant metabolites were mainly branched-chain amino acids (BCAAs), with some derivatised BCAAs, lipids, xenobiotics and a few unknowns. Metabolites such as lysophosphatidylcholine a C17:0, sum of hexoses, amino acids from BCAA metabolism (including leucine, isoleucine, valine, N-lactoylvaline, N-lactoylleucine and formiminoglutamate) and lactate, as well as an unknown metabolite (X-24295), were associated with HbA1c progression rate and were significant mediators of type 2 diabetes from baseline to 18 and 48 months of follow-up. 2SMR was used to estimate the causal effect of an exposure on an outcome using summary statistics from UK Biobank genome-wide association studies. We found that type 2 diabetes had a causal effect on the levels of three metabolites (hexose, glutamate and caproate [fatty acid (FA) 6:0]), whereas lipids such as specific phosphatidylcholines (PCs) (namely PC aa C36:2, PC aa C36:5, PC ae C36:3 and PC ae C34:3) as well as the two n-3 fatty acids stearidonate (18:4n3) and docosapentaenoate (22:5n3) potentially had a causal role in the development of type 2 diabetes. CONCLUSIONS/INTERPRETATION: Our findings identify known BCAAs and lipids, along with novel N-lactoyl-amino acid metabolites, significantly associated with prediabetes and diabetes, that mediate the effect of diabetes from baseline to follow-up (18 and 48 months). Causal inference using genetic variants shows the role of lipid metabolism and n-3 fatty acids as being causal for metabolite-to-type 2 diabetes whereas the sum of hexoses is causal for type 2 diabetes-to-metabolite. Identified metabolite markers are useful for stratifying individuals based on their risk progression and should enable targeted interventions.

Stool and blood metabolomics in the metabolic syndrome: a cross-sectional study.

INTRODUCTION/OBJECTIVES: Changes in the stool metabolome have been poorly studied in the metabolic syndrome (MetS). Moreover, few studies have explored the relationship of stool metabolites with circulating metabolites. Here, we investigated the associations between stool and blood metabolites, the MetS and systemic inflammation. METHODS: We analyzed data from 1,370 participants of the KORA FF4 study (Germany). Metabolites were measured by Metabolon, Inc. (untargeted) in stool, and using the AbsoluteIDQ® p180 kit (targeted) in blood. Multiple linear regression models, adjusted for dietary pattern, age, sex, physical activity, smoking status and alcohol intake, were used to estimate the associations of metabolites with the MetS, its components and high-sensitivity C-reactive protein (hsCRP) levels. Partial correlation and Multi-Omics Factor Analysis (MOFA) were used to investigate the relationship between stool and blood metabolites. RESULTS: The MetS was significantly associated with 170 stool and 82 blood metabolites. The MetS components with the highest number of associations were triglyceride levels (stool) and HDL levels (blood). Additionally, 107 and 27 MetS-associated metabolites (in stool and blood, respectively) showed significant associations with hsCRP levels. We found low partial correlation coefficients between stool and blood metabolites. MOFA did not detect shared variation across the two datasets. CONCLUSIONS: The MetS, particularly dyslipidemia, is associated with multiple stool and blood metabolites that are also associated with systemic inflammation. Further studies are necessary to validate our findings and to characterize metabolic alterations in the MetS. Although our analyses point to weak correlations between stool and blood metabolites, additional studies using integrative approaches are warranted.

Association of transketolase polymorphisms with diabetic polyneuropathy in the general population: The KORA F4 study.

AIMS: We recently reported that genetic variability in the TKT gene encoding transketolase, a key enzyme in the pentose phosphate pathway, is associated with measures of diabetic sensorimotor polyneuropathy (DSPN) in recent-onset diabetes. Here, we aimed to substantiate these findings in a population-based KORA F4 study. MATERIALS AND METHODS: In this cross-sectional study, we assessed seven single nucleotide polymorphisms (SNPs) in the transketolase gene in 952 participants from the KORA F4 study with normal glucose tolerance (NGT; n = 394), prediabetes (n = 411), and type 2 diabetes (n = 147). DSPN was defined by the examination part of the Michigan Neuropathy Screening Instrument (MNSI) using the original MNSI > 2 cut-off and two alternative versions extended by touch/pressure perception (TPP) (MNSI > 3) and by TPP plus cold perception (MNSI > 4). RESULTS: After adjustment for sex, age, BMI, and HbA1c, in type 2 diabetes participants, four out of seven transketolase SNPs were associated with DSPN for all three MNSI versions (all p ≤ 0.004). The odds ratios of these associations increased with extending the MNSI score, for example, OR (95% CI) for SNP rs62255988 with MNSI > 2: 1.99 (1.16-3.41), MNSI > 3: 2.27 (1.26-4.09), and MNSI > 4: 4.78 (2.22-10.26); SNP rs9284890 with MNSI > 2: 2.43 (1.42-4.16), MNSI > 3: 3.46 (1.82-6.59), and MNSI > 4: 4.75 (2.15-10.51). In contrast, no associations were found between transketolase SNPs and the three MNSI versions in the NGT and prediabetes groups. CONCLUSIONS: The link of genetic variation in transketolase enzyme to diabetic polyneuropathy corroborated at the population level strengthens the concept suggesting an important role of pathways metabolising glycolytic intermediates in the evolution of diabetic polyneuropathy.

Circulatory levels of lysophosphatidylcholine species in obese adolescents: Findings from cross-sectional and prospective lipidomics analyses.

BACKGROUND AND AIMS: Obesity has reached epidemic proportions, emphasizing the importance of reliable biomarkers for detecting early metabolic alterations and enabling early preventative interventions. However, our understanding of the molecular mechanisms and specific lipid species associated with childhood obesity remains limited. Therefore, the aim of this study was to investigate plasma lipidomic signatures as potential biomarkers for adolescent obesity. METHODS AND RESULTS: A total of 103 individuals comprising overweight/obese (n = 46) and normal weight (n = 57) were randomly chosen from the baseline ORANGE (Obesity Reduction and Noncommunicable Disease Awareness through Group Education) cohort, having been followed up for a median of 7.1 years. Plasma lipidomic profiling was performed using the UHPLC-HRMS method. We used three different models adjusted for clinical covariates to analyze the data. Clustering methods were used to define metabotypes, which allowed for the stratification of subjects into subgroups with similar clinical and metabolic profiles. We observed that lysophosphatidylcholine (LPC) species like LPC.16.0, LPC.18.3, LPC.18.1, and LPC.20.3 were significantly (p < 0.05) associated with baseline and follow-up BMI in adolescent obesity. The association of LPC species with BMI remained consistently significant even after adjusting for potential confounders. Moreover, applying metabotyping using hierarchical clustering provided insights into the metabolic heterogeneity within the normal and obese groups, distinguishing metabolically healthy individuals from those with unhealthy metabolic profiles. CONCLUSION: The specific LPC levels were found to be altered and increased in childhood obesity, particularly during the follow-up. These findings suggest that LPC species hold promise as potential biomarkers of obesity in adolescents, including healthy and unhealthy metabolic profiles.

Human genome diversity data reveal that L564P is the predominant TPC2 variant and a prerequisite for the blond hair associated M484L gain-of-function effect.

The endo-lysosomal two-pore channel (TPC2) has been established as an intracellular cation channel of significant physiological and pathophysiological relevance in recent years. For example, TPC2-/- mice show defects in cholesterol degradation, leading to hypercholesterinemia; TPC2 absence also results in mature-onset obesity, and a role in glucagon secretion and diabetes has been proposed. Infections with bacterial toxins or viruses e.g., cholera toxin or Ebola virus result in reduced infectivity rates in the absence of TPC2 or after pharmacological blockage, and TPC2-/- cancer cells lose their ability to migrate and metastasize efficiently. Finally, melanin production is affected by changes in hTPC2 activity, resulting in pigmentation defects and hair color variation. Here, we analyzed several publicly available genome variation data sets and identified multiple variations in the TPC2 protein in distinct human populations. Surprisingly, one variation, L564P, was found to be the predominant TPC2 isoform on a global scale. By applying endo-lysosomal patch-clamp electrophysiology, we found that L564P is a prerequisite for the previously described M484L gain-of-function effect that is associated with blond hair. Additionally, other gain-of-function variants with distinct geographical and ethnic distribution were discovered and functionally characterized. A meta-analysis of genome-wide association studies was performed, finding the polymorphisms to be associated with both distinct and overlapping traits. In sum, we present the first systematic analysis of variations in TPC2. We functionally characterized the most common variations and assessed their association with various disease traits. With TPC2 emerging as a novel drug target for the treatment of various diseases, this study provides valuable insights into ethnic and geographical distribution of TPC2 polymorphisms and their effects on channel activity.

Genome-wide association study of circulating interleukin 6 levels identifies novel loci.

Interleukin 6 (IL-6) is a multifunctional cytokine with both pro- and anti-inflammatory properties with a heritability estimate of up to 61%. The circulating levels of IL-6 in blood have been associated with an increased risk of complex disease pathogenesis. We conducted a two-staged, discovery and replication meta genome-wide association study (GWAS) of circulating serum IL-6 levels comprising up to 67 428 (ndiscovery = 52 654 and nreplication = 14 774) individuals of European ancestry. The inverse variance fixed effects based discovery meta-analysis, followed by replication led to the identification of two independent loci, IL1F10/IL1RN rs6734238 on chromosome (Chr) 2q14, (Pcombined = 1.8 × 10-11), HLA-DRB1/DRB5 rs660895 on Chr6p21 (Pcombined = 1.5 × 10-10) in the combined meta-analyses of all samples. We also replicated the IL6R rs4537545 locus on Chr1q21 (Pcombined = 1.2 × 10-122). Our study identifies novel loci for circulating IL-6 levels uncovering new immunological and inflammatory pathways that may influence IL-6 pathobiology.

A randomization-based causal inference framework for uncovering environmental exposure effects on human gut microbiota.

Statistical analysis of microbial genomic data within epidemiological cohort studies holds the promise to assess the influence of environmental exposures on both the host and the host-associated microbiome. However, the observational character of prospective cohort data and the intricate characteristics of microbiome data make it challenging to discover causal associations between environment and microbiome. Here, we introduce a causal inference framework based on the Rubin Causal Model that can help scientists to investigate such environment-host microbiome relationships, to capitalize on existing, possibly powerful, test statistics, and test plausible sharp null hypotheses. Using data from the German KORA cohort study, we illustrate our framework by designing two hypothetical randomized experiments with interventions of (i) air pollution reduction and (ii) smoking prevention. We study the effects of these interventions on the human gut microbiome by testing shifts in microbial diversity, changes in individual microbial abundances, and microbial network wiring between groups of matched subjects via randomization-based inference. In the smoking prevention scenario, we identify a small interconnected group of taxa worth further scrutiny, including Christensenellaceae and Ruminococcaceae genera, that have been previously associated with blood metabolite changes. These findings demonstrate that our framework may uncover potentially causal links between environmental exposure and the gut microbiome from observational data. We anticipate the present statistical framework to be a good starting point for further discoveries on the role of the gut microbiome in environmental health.

Mouse mutant phenotyping at scale reveals novel genes controlling bone mineral density.

The genetic landscape of diseases associated with changes in bone mineral density (BMD), such as osteoporosis, is only partially understood. Here, we explored data from 3,823 mutant mouse strains for BMD, a measure that is frequently altered in a range of bone pathologies, including osteoporosis. A total of 200 genes were found to significantly affect BMD. This pool of BMD genes comprised 141 genes with previously unknown functions in bone biology and was complementary to pools derived from recent human studies. Nineteen of the 141 genes also caused skeletal abnormalities. Examination of the BMD genes in osteoclasts and osteoblasts underscored BMD pathways, including vesicle transport, in these cells and together with in silico bone turnover studies resulted in the prioritization of candidate genes for further investigation. Overall, the results add novel pathophysiological and molecular insight into bone health and disease.

Epigenome-wide association study of incident type 2 diabetes: a meta-analysis of five prospective European cohorts.

AIMS/HYPOTHESIS: Type 2 diabetes is a complex metabolic disease with increasing prevalence worldwide. Improving the prediction of incident type 2 diabetes using epigenetic markers could help tailor prevention efforts to those at the highest risk. The aim of this study was to identify predictive methylation markers for incident type 2 diabetes by combining epigenome-wide association study (EWAS) results from five prospective European cohorts. METHODS: We conducted a meta-analysis of EWASs in blood collected 7-10 years prior to type 2 diabetes diagnosis. DNA methylation was measured with Illumina Infinium Methylation arrays. A total of 1250 cases and 1950 controls from five longitudinal cohorts were included: Doetinchem, ESTHER, KORA1, KORA2 and EPIC-Norfolk. Associations between DNA methylation and incident type 2 diabetes were examined using robust linear regression with adjustment for potential confounders. Inverse-variance fixed-effects meta-analysis of cohort-level individual CpG EWAS estimates was performed using METAL. The methylGSA R package was used for gene set enrichment analysis. Confirmation of genome-wide significant CpG sites was performed in a cohort of Indian Asians (LOLIPOP, UK). RESULTS: The meta-analysis identified 76 CpG sites that were differentially methylated in individuals with incident type 2 diabetes compared with control individuals (p values <1.1 × 10-7). Sixty-four out of 76 (84.2%) CpG sites were confirmed by directionally consistent effects and p values <0.05 in an independent cohort of Indian Asians. However, on adjustment for baseline BMI only four CpG sites remained genome-wide significant, and addition of the 76 CpG methylation risk score to a prediction model including established predictors of type 2 diabetes (age, sex, BMI and HbA1c) showed no improvement (AUC 0.757 vs 0.753). Gene set enrichment analysis of the full epigenome-wide results clearly showed enrichment of processes linked to insulin signalling, lipid homeostasis and inflammation. CONCLUSIONS/INTERPRETATION: By combining results from five European cohorts, and thus significantly increasing study sample size, we identified 76 CpG sites associated with incident type 2 diabetes. Replication of 64 CpGs in an independent cohort of Indian Asians suggests that the association between DNA methylation levels and incident type 2 diabetes is robust and independent of ethnicity. Our data also indicate that BMI partly explains the association between DNA methylation and incident type 2 diabetes. Further studies are required to elucidate the underlying biological mechanisms and to determine potential causal roles of the differentially methylated CpG sites in type 2 diabetes development.

Comparison of genetic risk prediction models to improve prediction of coronary heart disease in two large cohorts of the MONICA/KORA study.

It is still unclear how genetic information, provided as single-nucleotide polymorphisms (SNPs), can be most effectively integrated into risk prediction models for coronary heart disease (CHD) to add significant predictive value beyond clinical risk models. For the present study, a population-based case-cohort was used as a trainingset (451 incident cases, 1488 noncases) and an independent cohort as testset (160 incident cases, 2749 noncases). The following strategies to quantify genetic information were compared: A weighted genetic risk score including Metabochip SNPs associated with CHD in the literature (GRSMetabo ); selection of the most predictive SNPs among these literature-confirmed variants using priority-Lasso (PLMetabo ); validation of two comprehensive polygenic risk scores: GRSGola based on Metabochip data, and GRSKhera (available in the testset only) based on cross-validated genome-wide genotyping data. We used Cox regression to assess associations with incident CHD. C-index, category-free net reclassification index (cfNRI) and relative integrated discrimination improvement (IDIrel ) were used to quantify the predictive performance of genetic information beyond Framingham risk score variables. In contrast to GRSMetabo and PLMetabo , GRSGola significantly improved the prediction (delta C-index [95% confidence interval]: 0.0087 [0.0044, 0.0130]; IDIrel : 0.0509 [0.0131, 0.0894]; cfNRI improved only in cases: 0.1761 [0.0253, 0.3219]). GRSKhera yielded slightly worse prediction results than GRSGola .

A Family and a Genome-Wide Polygenic Risk Score Are Independently Associated With Stroke in a Population-Based Study.

BACKGROUND: Positive family history and genetic risk scores have been shown to independently capture those individuals with high risk for stroke. The aim of our study was to evaluate the amount of shared information between family history and genetic risk and to investigate their combined effect on the association with prevalent and incident stroke cases. METHODS: We obtained a family risk score (FamRS), weighted for disease onset and family size as well as genome-wide polygenic risk score (PGS) including over 3.2 million single-nucleotide polymorphisms in the population-based prospective KORA F3 (Cooperative Health Research in the Region of Augsburg) study (n=3071) from Southern Germany. FamRS and PGS were evaluated separately and combined. The measures were once treated as continuous variables but also divided in the highest 20%, 10%, 5%, and 1% percentiles. Odds ratios via logistic regression and hazard ratios via Cox regression were estimated. A stroke event was defined as a hospitalization for stroke that was self-reported in a standardized interview by certified and supervised personnel. RESULTS: The FamRS outperformed other simplified family measures such as affected parents or number of affected family members. FamRS and PGS were not correlated, and no individuals were observed with both very high FamRS and very high PGS (top 1% percentile). In a combined model, both FamRS and PGS were independently from each other associated with risk of stroke, also independent of other traditional risk factors (p [FamRS]=0.02, p [PGS]=0.005). Individuals in the top 1% of either FamRS or PGS were found to have >5-fold risk for stroke (odds ratios, 5.82 [95% CI, 2.08-14]; P=0.0002). The results for incident stroke events showed the same trend but were not significant. CONCLUSIONS: Our study shows that a family risk score and PGS capture different information concerning individual stroke risk. Combining the risk measures FamRS and PGS increases predictive power, as demonstrated in a population-based study.

Exome-Derived Adiponectin-Associated Variants Implicate Obesity and Lipid Biology.

Circulating levels of adiponectin, an adipocyte-secreted protein associated with cardiovascular and metabolic risk, are highly heritable. To gain insights into the biology that regulates adiponectin levels, we performed an exome array meta-analysis of 265,780 genetic variants in 67,739 individuals of European, Hispanic, African American, and East Asian ancestry. We identified 20 loci associated with adiponectin, including 11 that had been reported previously (p < 2 × 10-7). Comparison of exome array variants to regional linkage disequilibrium (LD) patterns and prior genome-wide association study (GWAS) results detected candidate variants (r2 > .60) spanning as much as 900 kb. To identify potential genes and mechanisms through which the previously unreported association signals act to affect adiponectin levels, we assessed cross-trait associations, expression quantitative trait loci in subcutaneous adipose, and biological pathways of nearby genes. Eight of the nine loci were also associated (p < 1 × 10-4) with at least one obesity or lipid trait. Candidate genes include PRKAR2A, PTH1R, and HDAC9, which have been suggested to play roles in adipocyte differentiation or bone marrow adipose tissue. Taken together, these findings provide further insights into the processes that influence circulating adiponectin levels.

Colocalization analysis of pancreas eQTLs with risk loci from alcoholic and novel non-alcoholic chronic pancreatitis GWAS suggests potential disease causing mechanisms.

BACKGROUND: Previous genome-wide association studies (GWAS) identified genome-wide significant risk loci in chronic pancreatitis and investigated underlying disease causing mechanisms by simple overlaps with expression quantitative trait loci (eQTLs), a procedure which may often result in false positive conclusions. METHODS: We conducted a GWAS in 584 non-alcoholic chronic pancreatitis (NACP) patients and 6040 healthy controls. Next, we applied Bayesian colocalization analysis of identified genome-wide significant risk loci from both, our recently published alcoholic chronic pancreatitis (ACP) and the novel NACP dataset, with pancreas eQTLs from the GTEx V8 European cohort to prioritize candidate causal genes and extracted credible sets of shared causal variants. RESULTS: Variants at the CTRC (p = 1.22 × 10-21) and SPINK1 (p = 6.59 × 10-47) risk loci reached genome-wide significance in NACP. CTRC risk variants colocalized with CTRC eQTLs in ACP (PP4 = 0.99, PP4/PP3 = 95.51) and NACP (PP4 = 0.99, PP4/PP3 = 95.46). For both diseases, the 95% credible set of shared causal variants consisted of rs497078 and rs545634. CLDN2-MORC4 risk variants colocalized with CLDN2 eQTLs in ACP (PP4 = 0.98, PP4/PP3 = 42.20) and NACP (PP4 = 0.67, PP4/PP3 = 7.18), probably driven by the shared causal variant rs12688220. CONCLUSIONS: A shared causal CTRC risk variant might unfold its pathogenic effect in ACP and NACP by reducing CTRC expression, while the CLDN2-MORC4 shared causal variant rs12688220 may modify ACP and NACP risk by increasing CLDN2 expression.

Metabolic syndrome and the plasma proteome: from association to causation.

BACKGROUND: The metabolic syndrome (MetS), defined by the simultaneous clustering of cardio-metabolic risk factors, is a significant worldwide public health burden with an estimated 25% prevalence worldwide. The pathogenesis of MetS is not entirely clear and the use of molecular level data could help uncover common pathogenic pathways behind the observed clustering. METHODS: Using a highly multiplexed aptamer-based affinity proteomics platform, we examined associations between plasma proteins and prevalent and incident MetS in the KORA cohort (n = 998) and replicated our results for prevalent MetS in the HUNT3 study (n = 923). We applied logistic regression models adjusted for age, sex, smoking status, and physical activity. We used the bootstrap ranking algorithm of least absolute shrinkage and selection operator (LASSO) to select a predictive model from the incident MetS associated proteins and used area under the curve (AUC) to assess its performance. Finally, we investigated the causal effect of the replicated proteins on MetS using two-sample Mendelian randomization. RESULTS: Prevalent MetS was associated with 116 proteins, of which 53 replicated in HUNT. These included previously reported proteins like leptin, and new proteins like NTR domain-containing protein 2 and endoplasmic reticulum protein 29. Incident MetS was associated with 14 proteins in KORA, of which 13 overlap the prevalent MetS associated proteins with soluble advanced glycosylation end product-specific receptor (sRAGE) being unique to incident MetS. The LASSO selected an eight-protein predictive model with an (AUC = 0.75; 95% CI = 0.71-0.79) in KORA. Mendelian randomization suggested causal effects of three proteins on MetS, namely apolipoprotein E2 (APOE2) (Wald-Ratio = - 0.12, Wald-p = 3.63e-13), apolipoprotein B (APOB) (Wald-Ratio = - 0.09, Wald-p = 2.54e-04) and proto-oncogene tyrosine-protein kinase receptor (RET) (Wald-Ratio = 0.10, Wald-p = 5.40e-04). CONCLUSIONS: Our findings offer new insights into the plasma proteome underlying MetS and identify new protein associations. We reveal possible casual effects of APOE2, APOB and RET on MetS. Our results highlight protein candidates that could potentially serve as targets for prevention and therapy.

Genome-wide meta-analysis identifies novel determinants of circulating serum progranulin.

Progranulin is a secreted protein with important functions in processes including immune and inflammatory response, metabolism and embryonic development. The present study aimed at identification of genetic factors determining progranulin concentrations. We conducted a genome-wide association meta-analysis for serum progranulin in three independent cohorts from Europe: Sorbs (N = 848) and KORA (N = 1628) from Germany and PPP-Botnia (N = 335) from Finland (total N = 2811). Single nucleotide polymorphisms (SNPs) associated with progranulin levels were replicated in two additional German cohorts: LIFE-Heart Study (Leipzig; N = 967) and Metabolic Syndrome Berlin Potsdam (Berlin cohort; N = 833). We measured mRNA expression of genes in peripheral blood mononuclear cells (PBMC) by micro-arrays and performed mRNA expression quantitative trait and expression-progranulin association studies to functionally substantiate identified loci. Finally, we conducted siRNA silencing experiments in vitro to validate potential candidate genes within the associated loci. Heritability of circulating progranulin levels was estimated at 31.8% and 26.1% in the Sorbs and LIFE-Heart cohort, respectively. SNPs at three loci reached study-wide significance (rs660240 in CELSR2-PSRC1-MYBPHL-SORT1, rs4747197 in CDH23-PSAP and rs5848 in GRN) explaining 19.4%/15.0% of the variance and 61%/57% of total heritability in the Sorbs/LIFE-Heart Study. The strongest evidence for association was at rs660240 (P = 5.75 × 10-50), which was also associated with mRNA expression of PSRC1 in PBMC (P = 1.51 × 10-21). Psrc1 knockdown in murine preadipocytes led to a consecutive 30% reduction in progranulin secretion. In conclusion, the present meta-GWAS combined with mRNA expression identified three loci associated with progranulin and supports the role of PSRC1 in the regulation of progranulin secretion.