comorbidPGS: an R package assessing shared predisposition between Phenotypes using Polygenic Scores.

Introduction Polygenic Score (PGS) is a valuable method for assessing the estimated genetic liability to a given outcome or genetic variability contributing to a quantitative trait. While PRSs are widely used for complex traits, their application in uncovering shared genetic predisposition between phenotypes, i.e. when genetic variants influence more than one phenotype, remains limited. Methods We developed an R package, comorbidPGS, which facilitates a systematic evaluation of shared genetic effects among (cor)related phenotypes using PGSs. The comorbidPGS package takes as input a set of Single Nucleotide Polymorphisms (SNPs) along with their established effects on the original phenotype (Po), referred to as Po-PGS. It generates a comprehensive summary of effect(s) of Po-PGS on target phenotype(s) (Pt) with customisable graphical features. Results We applied comorbidPGS to investigate the shared genetic predisposition between phenotypes defining elevated blood pressure (Systolic Blood Pressure, SBP; Diastolic Blood Pressure, DBP; Pulse Pressure, PP) and several cancers (Breast Cancer, BrC; Pancreatic Cancer, PanC; Kidney Cancer, KidC; Prostate Cancer, PrC; Colorectal Cancer, CrC) using the European ancestry UK Biobank individuals and GWAS meta-analyses summary statistics from independent set of European ancestry individuals. We report a significant association between elevated DBP and the genetic risk of PrC (ß (SE)=0.066 (0.017), P-value=9.64×10^(-5)), as well as between CrC PGS and both, lower SBP (ß (SE)=-0.10 [0.029], P-value=3.83×10^(-4))) and lower DBP (ß (SE)=-0.055 [0.017], P-value=1.05×10^(-3)). Our analysis highlights two nominally significant relationships for individuals with genetic predisposition to elevated SBP leading to higher risk of KidC (OR [95%CI]=1.04 [1.0039-1.087], P-value=2.82×10^(-2)) and PrC (OR [95%CI]=1.02 [1.003-1.041], P-value=2.22×10^(-2)). Conclusion Using comorbidPGS, we underscore mechanistic relationships between blood pressure regulation and susceptibility to three comorbid malignancies. This package offers valuable means to evaluate shared genetic susceptibility between (cor)related phenotypes through polygenic scores.

A framework for conducting time-varying genome-wide association studies: An application to body mass index across childhood in six multiethnic cohorts.

Genetic effects on changes in human traits over time are understudied and may have important pathophysiological impact. We propose a framework that enables data quality control, implements mixed models to evaluate trajectories of change in traits, and estimates phenotypes to identify age-varying genetic effects in genome-wide association studies (GWASs). Using childhood body mass index (BMI) as an example, we included 71,336 participants from six cohorts and estimated the slope and area under the BMI curve within four time periods (infancy, early childhood, late childhood and adolescence) for each participant, in addition to the age and BMI at the adiposity peak and the adiposity rebound. GWAS on each of the estimated phenotypes identified 28 genome-wide significant variants at 13 loci across the 12 estimated phenotypes, one of which was novel (in DAOA) and had not been previously associated with childhood or adult BMI. Genetic studies of changes in human traits over time could uncover novel biological mechanisms influencing quantitative traits.

Multiple genetic variants at the SLC30A8 locus affect local super-enhancer activity and influence pancreatic ß-cell survival and function.

Variants at the SLC30A8 locus are associated with type 2 diabetes (T2D) risk. The lead variant, rs13266634, encodes an amino acid change, Arg325Trp (R325W), at the C-terminus of the secretory granule-enriched zinc transporter, ZnT8. Although this protein-coding variant was previously thought to be the sole driver of T2D risk at this locus, recent studies have provided evidence for lowered expression of SLC30A8 mRNA in protective allele carriers. In the present study, we examined multiple variants that influence SLC30A8 allele-specific expression. Epigenomic mapping has previously identified an islet-selective enhancer cluster at the SLC30A8 locus, hosting multiple T2D risk and cASE associations, which is spatially associated with the SLC30A8 promoter and additional neighboring genes. Here, we show that deletion of variant-bearing enhancer regions using CRISPR-Cas9 in human-derived EndoC-ßH3 cells lowers the expression of SLC30A8 and several neighboring genes and improves glucose-stimulated insulin secretion. While downregulation of SLC30A8 had no effect on beta cell survival, loss of UTP23, RAD21, or MED30 markedly reduced cell viability. Although eQTL or cASE analyses in human islets did not support the association between these additional genes and diabetes risk, the transcriptional regulator JQ1 lowered the expression of multiple genes at the SLC30A8 locus and enhanced stimulated insulin secretion.

Pathogenic, Total Loss-of-Function DYRK1B Variants Cause Monogenic Obesity Associated With Type 2 Diabetes.

OBJECTIVE: Rare variants in DYRK1B have been described in some patients with central obesity, type 2 diabetes, and early-onset coronary disease. Owing to the limited number of conducted studies, the broader impact of DYRK1B variants on a larger scale has yet to be investigated. RESEARCH DESIGN AND METHODS: DYRK1B was sequenced in 9,353 participants from a case-control study for obesity and type 2 diabetes. Each DYRK1B variant was functionally assessed in vitro. Variant pathogenicity was determined using criteria from the American College of Medical Genetics and Genomics (ACMG). The effect of pathogenic or likely pathogenic (P/LP) variants on metabolic traits was assessed using adjusted mixed-effects score tests. RESULTS: Sixty-five rare, heterozygous DYRK1B variants were identified and were not associated with obesity or type 2 diabetes. Following functional analyses, 20 P/LP variants were pinpointed, including 6 variants that exhibited a fully inhibitory effect (P/LP-null) on DYRK1B activity. P/LP and P/LP-null DYRK1B variants were associated with increased BMI and obesity risk; however, the impact was notably more pronounced for the P/LP-null variants (effect of 8.0 ± 3.2 and odds ratio of 7.9 [95% CI 1.2-155]). Furthermore, P/LP-null variants were associated with higher fasting glucose and type 2 diabetes risk (effect of 2.9 ± 1.0 and odds ratio of 4.8 [95% CI 0.85-37]), while P/LP variants had no effect on glucose homeostasis. CONCLUSIONS: P/LP, total loss-of-function DYRK1B variants cause monogenic obesity associated with type 2 diabetes. This study underscores the significance of conducting functional assessments in order to accurately ascertain the tangible effects of P/LP DYRK1B variants.

Pathogenic monoallelic variants in GLIS3 increase type 2 diabetes risk and identify a subgroup of patients sensitive to sulfonylureas.

AIMS/HYPOTHESIS: GLIS3 encodes a transcription factor involved in pancreatic beta cell development and function. Rare pathogenic, bi-allelic mutations in GLIS3 cause syndromic neonatal diabetes whereas frequent SNPs at this locus associate with common type 2 diabetes risk. Because rare, functional variants located in other susceptibility genes for type 2 diabetes have already been shown to strongly increase individual risk for common type 2 diabetes, we aimed to investigate the contribution of rare pathogenic GLIS3 variants to type 2 diabetes. METHODS: GLIS3 was sequenced in 5471 individuals from the Rare Variants Involved in Diabetes and Obesity (RaDiO) study. Variant pathogenicity was assessed following the criteria established by the American College of Medical Genetics and Genomics (ACMG). To address the pathogenic strong criterion number 3 (PS3), we conducted functional investigations of these variants using luciferase assays, focusing on capacity of GLIS family zinc finger 3 (GLIS3) to bind to and activate the INS promoter. The association between rare pathogenic or likely pathogenic (P/LP) variants and type 2 diabetes risk (and other metabolic traits) was then evaluated. A meta-analysis combining association results from RaDiO, the 52K study (43,125 individuals) and the TOPMed study (44,083 individuals) was finally performed. RESULTS: Through targeted resequencing of GLIS3, we identified 105 rare variants that were carried by 395 participants from RaDiO. Among them, 49 variants decreased the activation of the INS promoter. Following ACMG criteria, 18 rare variants were classified as P/LP, showing an enrichment in the last two exons compared with the remaining exons (p<5×10-6; OR>3.5). The burden of these P/LP variants was strongly higher in individuals with type 2 diabetes (p=3.0×10-3; OR 3.9 [95% CI 1.4, 12]), whereas adiposity, age at type 2 diabetes diagnosis and cholesterol levels were similar between variant carriers and non-carriers with type 2 diabetes. Interestingly, all carriers with type 2 diabetes were sensitive to oral sulfonylureas. A total of 7 P/LP variants were identified in both 52K and TOPMed studies. The meta-analysis of association studies obtained from RaDiO, 52K and TOPMed showed an enrichment of P/LP GLIS3 variants in individuals with type 2 diabetes (p=5.6×10-5; OR 2.1 [95% CI 1.4, 2.9]). CONCLUSIONS/INTERPRETATION: Rare P/LP GLIS3 variants do contribute to type 2 diabetes risk. The variants located in the distal part of the protein could have a direct effect on its functional activity by impacting its transactivation domain, by homology with the mouse GLIS3 protein. Furthermore, rare P/LP GLIS3 variants seem to have a direct clinical effect on beta cell function, which could be improved by increasing insulin secretion via the use of sulfonylureas.

Dominant PDX1 deficiency causes highly penetrant diabetes at different ages, associated with obesity and exocrine pancreatic deficiency: Lessons for precision medicine.

OBJECTIVE: Heterozygous pathogenic or likely pathogenic (P/LP) PDX1 variants cause monogenic diabetes. We comprehensively examined the phenotypes of carriers of P/LP PDX1 variants, and delineated potential treatments that could be efficient in an objective of precision medicine. METHODS: The study primarily involved a family harboring a novel P/LP PDX1 variant. We then conducted an analysis of documented carriers of P/LP PDX1 variants, from the Human Gene Mutation Database (HGMD), RaDiO study, and Type 2 Diabetes Knowledge Portal (T2DKP) including 87 K participants. RESULTS: Within the family, we identified a P/LP PDX1 variant encoding p.G232S in four relatives. All of them exhibited diabetes, albeit with very different ages of onset (10-40 years), along with caudal pancreatic agenesis and childhood-onset obesity. In the HGMD, 79 % of carriers of a P/LP PDX1 variant displayed diabetes (with differing ages of onset from eight days of life to 67 years), 63 % exhibited pancreatic insufficiency and surprisingly 40 % had obesity. The impact of P/LP PDX1 variants on increased risk of type 2 diabetes mellitus was confirmed in the T2DKP. Dipeptidyl peptidase 4 inhibitor (DPP4i) and glucagon-like peptide-1 receptor agonist (GLP1-RA), enabled good glucose control without hypoglycemia and weight management. CONCLUSIONS: This study reveals diverse clinical presentations among the carriers of a P/LP PDX1 variant, highlighting strong variations in diabetes onset, and unexpectedly high prevalence of obesity and pancreatic development abnormalities. Clinical data suggest that DPP4i and GLP1-RA may be the best effective treatments to manage both glucose and weight controls, opening new avenue in precision diabetic medicine.

Exploring the role of purinergic receptor P2RY1 in type 2 diabetes risk and pathophysiology: Insights from human functional genomics.

OBJECTIVE: Human functional genomics has proven powerful in discovering drug targets for common metabolic disorders. Through this approach, we investigated the involvement of the purinergic receptor P2RY1 in type 2 diabetes (T2D). METHODS: P2RY1 was sequenced in 9,266 participants including 4,177 patients with T2D. In vitro analyses were then performed to assess the functional effect of each variant. Expression quantitative trait loci (eQTL) analysis was performed in pancreatic islets from 103 pancreatectomized individuals. The effect of P2RY1 on glucose-stimulated insulin secretion was finally assessed in human pancreatic beta cells (EndoCßH5), and RNA sequencing was performed on these cells. RESULTS: Sequencing P2YR1 in 9,266 participants revealed 22 rare variants, seven of which were loss-of-function according to our in vitro analyses. Carriers, except one, exhibited impaired glucose control. Our eQTL analysis of human islets identified P2RY1 variants, in a beta-cell enhancer, linked to increased P2RY1 expression and reduced T2D risk, contrasting with variants located in a silent region associated with decreased P2RY1 expression and increased T2D risk. Additionally, a P2RY1-specific agonist increased insulin secretion upon glucose stimulation, while the antagonist led to decreased insulin secretion. RNA-seq highlighted TXNIP as one of the main transcriptomic markers of insulin secretion triggered by P2RY1 agonist. CONCLUSION: Our findings suggest that P2RY1 inherited or acquired dysfunction increases T2D risk and that P2RY1 activation stimulates insulin secretion. Selective P2RY1 agonists, impermeable to the blood-brain barrier, could serve as potential insulin secretagogues.

Time-of-day-dependent variation of the human liver transcriptome and metabolome is disrupted in MASLD.

Background & Aims: Liver homeostasis is ensured in part by time-of-day-dependent processes, many of them being paced by the molecular circadian clock. Liver functions are compromised in metabolic dysfunction-associated steatotic liver disease (MASLD) and metabolic dysfunction-associated steatohepatitis (MASH), and clock disruption increases susceptibility to MASLD progression in rodent models. We therefore investigated whether the time-of-day-dependent transcriptome and metabolome are significantly altered in human steatotic and MASH livers. Methods: Liver biopsies, collected within an 8 h-window from a carefully phenotyped cohort of 290 patients and histologically diagnosed to be either normal, steatotic or MASH hepatic tissues, were analyzed by RNA sequencing and unbiased metabolomic approaches. Time-of-day-dependent gene expression patterns and metabolomes were identified and compared between histologically normal, steatotic and MASH livers. Results: Herein, we provide a first-of-its-kind report of a daytime-resolved human liver transcriptome-metabolome and associated alterations in MASLD. Transcriptomic analysis showed a robustness of core molecular clock components in steatotic and MASH livers. It also revealed stage-specific, time-of-day-dependent alterations of hundreds of transcripts involved in cell-to-cell communication, intracellular signaling and metabolism. Similarly, rhythmic amino acid and lipid metabolomes were affected in pathological livers. Both TNFα and PPARγ signaling were predicted as important contributors to altered rhythmicity. Conclusion: MASLD progression to MASH perturbs time-of-day-dependent processes in human livers, while the differential expression of core molecular clock components is maintained. Impact and implications: This work characterizes the rhythmic patterns of the transcriptome and metabolome in the human liver. Using a cohort of well-phenotyped patients (n = 290) for whom the time-of-day at biopsy collection was known, we show that time-of-day variations observed in histologically normal livers are gradually perturbed in liver steatosis and metabolic dysfunction-associated steatohepatitis. Importantly, these observations, albeit obtained across a restricted time window, provide further support for preclinical studies demonstrating alterations of rhythmic patterns in diseased livers. On a practical note, this study indicates the importance of considering time-of-day as a critical biological variable which may significantly affect data interpretation in animal and human studies of liver diseases.

Genetic variants of interferon-response factor 5 are associated with the incidence of chronic kidney disease: the D.E.S.I.R. study.

Inflammation has been associated with renal diseases. The Interferon Regulatory Factor (IRF)-5 is a key transcription factor in the pro-inflammatory polarization of M1-like macrophages. GWAS have reported that the IRF5 locus is associated with autoimmune diseases and with the estimated glomerular filtration rate (eGFR). We study whether allelic variations in IRF5 are associated with the incidence of chronic kidney disease (CKD) in a general population. We genotyped eleven IRF5 SNPs in the French D.E.S.I.R. cohort from the general population (n = 4820). Associations of SNPs with baseline renal parameters were assessed. Data were analyzed for three endpoints during a 9-year follow-up, incidence of:at least stage 3 CKD, the KDIGO criterion "certain drop in eGFR", and incidence of micro/macro albuminuria. In the cross-sectional analysis, rs10954213 and rs10954214 were associated with eGFR and rs1874328 with urinary albumin/creatinine ratio (ACR). Rs3807306, rs11761199, rs78658945, rs1874328, rs10954213 and rs11770589 were associated with the incidence of stage 3 CKD in multi-adjusted models. Rs4731532, rs3807306, and rs11761199 were associated with the incidence of CKD defined by the KDIGO. Rs4731532, rs3807306, rs11761199 and rs79288514 were associated with the incidence of micro/macro albuminuria. Our results support the hypothesis of the importance of IRF5 mediated macrophage polarization in the etiology of CKD.

High morbidity and mortality in children with untreated congenital deficiency of leptin or its receptor.

The long-term clinical outcomes of severe obesity due to leptin signaling deficiency are unknown. We carry out a retrospective cross-sectional investigation of a large cohort of children with leptin (LEP), LEP receptor (LEPR), or melanocortin 4 receptor (MC4R) deficiency (n = 145) to evaluate the progression of the disease. The affected individuals undergo physical, clinical, and metabolic evaluations. We report a very high mortality in children with LEP (26%) or LEPR deficiency (9%), mainly due to severe pulmonary and gastrointestinal infections. In addition, 40% of surviving children with LEP or LEPR deficiency experience life-threatening episodes of lung or gastrointestinal infections. Although precision drugs are currently available for LEP and LEPR deficiencies, as yet, they are not accessible in Pakistan. An appreciation of the severe impact of LEP or LEPR deficiency on morbidity and early mortality, educational attainment, and the attendant stigmatization should spur efforts to deliver the available life-saving drugs to these children as a matter of urgency.

Evidence of a causal and modifiable relationship between kidney function and circulating trimethylamine N-oxide.

The host-microbiota co-metabolite trimethylamine N-oxide (TMAO) is linked to increased cardiovascular risk but how its circulating levels are regulated remains unclear. We applied "explainable" machine learning, univariate, multivariate and mediation analyses of fasting plasma TMAO concentration and a multitude of phenotypes in 1,741 adult Europeans of the MetaCardis study. Here we show that next to age, kidney function is the primary variable predicting circulating TMAO, with microbiota composition and diet playing minor, albeit significant, roles. Mediation analysis suggests a causal relationship between TMAO and kidney function that we corroborate in preclinical models where TMAO exposure increases kidney scarring. Consistent with our findings, patients receiving glucose-lowering drugs with reno-protective properties have significantly lower circulating TMAO when compared to propensity-score matched control individuals. Our analyses uncover a bidirectional relationship between kidney function and TMAO that can potentially be modified by reno-protective anti-diabetic drugs and suggest a clinically actionable intervention for decreasing TMAO-associated excess cardiovascular risk.

Human GLP1R variants affecting GLP1R cell surface expression are associated with impaired glucose control and increased adiposity.

The glucagon-like peptide 1 receptor (GLP1R) is a major drug target with several agonists being prescribed in individuals with type 2 diabetes and obesity1,2. The impact of genetic variability of GLP1R on receptor function and its association with metabolic traits are unclear with conflicting reports. Here, we show an unexpected diversity of phenotypes ranging from defective cell surface expression to complete or pathway-specific gain of function (GoF) and loss of function (LoF), after performing a functional profiling of 60 GLP1R variants across four signalling pathways. The defective insulin secretion of GLP1R LoF variants is rescued by allosteric GLP1R ligands or high concentrations of exendin-4/semaglutide in INS-1 823/3 cells. Genetic association studies in 200,000 participants from the UK Biobank show that impaired GLP1R cell surface expression contributes to poor glucose control and increased adiposity with increased glycated haemoglobin A1c and body mass index. This study defines impaired GLP1R cell surface expression as a risk factor for traits associated with type 2 diabetes and obesity and provides potential treatment options for GLP1R LoF variant carriers.

Identification of myokines susceptible to improve glucose homeostasis after bariatric surgery.

IMPORTANCE AND OBJECTIVE: The identification of myokines susceptible to improve glucose homeostasis following bariatric surgery could lead to new therapeutic approaches for type 2 diabetes. METHODS: Changes in the homeostasis model assessment (HOMA) test were assessed in patients before and 3 months after bariatric surgery. Changes in myokines expression and circulating levels were assessed using real-time quantitative polymerase chain reaction (RT-qPCR) and enzyme-linked immunosorbent assay (ELISA). Myokines known to regulate glucose homeostasis were identified using literature (targeted study) and putative myokines using RNA-sequencing (untargeted study). A linear regression analysis adjusted for age and sex was used to search for associations between changes in the HOMA test and changes in myokines. RESULTS: In the targeted study, brain-derived neurotrophic factor (BDNF) expression was upregulated (+30%, P = .006) while BDNF circulating levels were decreased (-12%, P = .001). Upregulated BDNF expression was associated with decreased HOMA of insulin resistance (HOMA-IR) (adjusted estimate [95% confidence interval {CI}]: -0.51 [-0.88 to -0.13], P = .010). Decreased BDNF serum levels were associated with decreased HOMA of beta-cell function (HOMA-B) (adjusted estimate [95% CI] = 0.002 [0.00002-0.0031], P = .046). In the untargeted study, upregulated putative myokines included XYLT1 (+64%, P < .001), LGR5 (+57, P< .001), and SPINK5 (+46%, P < .001). Upregulated LGR5 was associated with decreased HOMA-IR (adjusted estimate [95% CI] = -0.50 [-0.86 to -0.13], P = .009). Upregulated XYLT1 and SPINK5 were associated with increased HOMA of insulin sensitivity (HOMA-S) (respectively, adjusted estimate [95% CI] = 109.1 [28.5-189.8], P = .009 and 16.5 [0.87-32.19], P = .039). CONCLUSIONS: Improved glucose homeostasis following bariatric surgery is associated with changes in myokines expression and circulating levels. In particular, upregulation of BDNF, XYLT1, SPINK5, and LGR5 is associated with improved insulin sensitivity. These results suggest that these myokines could contribute to improved glucose homeostasis following bariatric surgery. STUDY REGISTRATION: NCT03341793 on ClinicalTrials.gov (https://clinicaltrials.gov/).

Multiple genetic variants at the SLC30A8 locus affect local super-enhancer activity and influence pancreatic ß-cell survival and function.

Variants at the SLC30A8 locus are associated with type 2 diabetes (T2D) risk. The lead variant, rs13266634, encodes an amino acid change, Arg325Trp (R325W), at the C-terminus of the secretory granule-enriched zinc transporter, ZnT8. Although this protein-coding variant was previously thought to be the sole driver of T2D risk at this locus, recent studies have provided evidence for lowered expression of SLC30A8 mRNA in protective allele carriers. In the present study, combined allele-specific expression (cASE) analysis in human islets revealed multiple variants that influence SLC30A8 expression. Epigenomic mapping identified an islet-selective enhancer cluster at the SLC30A8 locus, hosting multiple T2D risk and cASE associations, which is spatially associated with the SLC30A8 promoter and additional neighbouring genes. Deletions of variant-bearing enhancer regions using CRISPR-Cas9 in human-derived EndoC-ßH3 cells lowered the expression of SLC30A8 and several neighbouring genes, and improved insulin secretion. Whilst down-regulation of SLC30A8 had no effect on beta cell survival, loss of UTP23, RAD21 or MED30 markedly reduced cell viability. Although eQTL or cASE analyses in human islets did not support the association between these additional genes and diabetes risk, the transcriptional regulator JQ1 lowered the expression of multiple genes at the SLC30A8 locus and enhanced stimulated insulin secretion.

Pharmacological HDAC inhibition impairs pancreatic ß-cell function through an epigenome-wide reprogramming.

Histone deacetylases enzymes (HDACs) are chromatin modifiers that regulate gene expression through deacetylation of lysine residues within specific histone and non-histone proteins. A cell-specific gene expression pattern defines the identity of insulin-producing pancreatic ß cells, yet molecular networks driving this transcriptional specificity are not fully understood. Here, we investigated the HDAC-dependent molecular mechanisms controlling pancreatic ß-cell identity and function using the pan-HDAC inhibitor trichostatin A through chromatin immunoprecipitation assays and RNA sequencing experiments. We observed that TSA alters insulin secretion associated with ß-cell specific transcriptome programming in both mouse and human ß-cell lines, as well as on human pancreatic islets. We also demonstrated that this alternative ß-cell transcriptional program in response to HDAC inhibition is related to an epigenome-wide remodeling at both promoters and enhancers. Our data indicate that HDAC activity could be required to protect against loss of ß-cell identity with unsuitable expression of genes associated with alternative cell fates.

Bidirectional Mendelian Randomization and Multiphenotype GWAS Show Causality and Shared Pathophysiology Between Depression and Type 2 Diabetes.

OBJECTIVE: Depression is a common comorbidity of type 2 diabetes. We assessed the causal relationships and shared genetics between them. RESEARCH DESIGN AND METHODS: We applied two-sample, bidirectional Mendelian randomization (MR) to assess causality between type 2 diabetes and depression. We investigated potential mediation using two-step MR. To identify shared genetics, we performed 1) genome-wide association studies (GWAS) separately and 2) multiphenotype GWAS (MP-GWAS) of type 2 diabetes (19,344 case subjects, 463,641 control subjects) and depression using major depressive disorder (MDD) (5,262 case subjects, 86,275 control subjects) and self-reported depressive symptoms (n = 153,079) in the UK Biobank. We analyzed expression quantitative trait locus (eQTL) data from public databases to identify target genes in relevant tissues. RESULTS: MR demonstrated a significant causal effect of depression on type 2 diabetes (odds ratio 1.26 [95% CI 1.11-1.44], P = 5.46 × 10-4) but not in the reverse direction. Mediation analysis indicated that 36.5% (12.4-57.6%, P = 0.0499) of the effect from depression on type 2 diabetes was mediated by BMI. GWAS of type 2 diabetes and depressive symptoms did not identify shared loci. MP-GWAS identified seven shared loci mapped to TCF7L2, CDKAL1, IGF2BP2, SPRY2, CCND2-AS1, IRS1, CDKN2B-AS1. MDD has not brought any significant association in either GWAS or MP-GWAS. Most MP-GWAS loci had an eQTL, including single nucleotide polymorphisms implicating the cell cycle gene CCND2 in pancreatic islets and brain and the insulin signaling gene IRS1 in adipose tissue, suggesting a multitissue and pleiotropic underlying mechanism. CONCLUSIONS: Our results highlight the importance to prevent type 2 diabetes at the onset of depressive symptoms and the need to maintain a healthy weight in the context of its effect on depression and type 2 diabetes comorbidity.

Abdominal obesity is a more important causal risk factor for pancreatic cancer than overall obesity.

Obesity and type 2 diabetes (T2D) are associated with increased risk of pancreatic cancer. Here we assessed the relationship between pancreatic cancer and two distinct measures of obesity, namely total adiposity, using BMI, versus abdominal adiposity, using BMI adjusted waist-to-hip ratio (WHRadjBMI) by utilising polygenic scores (PGS) and Mendelian randomisation (MR) analyses. We constructed z-score weighted PGS for BMI and WHRadjBMI using publicly available data and tested for their association with pancreatic cancer defined in UK biobank (UKBB). Using publicly available summary statistics, we then performed bi-directional MR analyses between the two obesity traits and pancreatic cancer. PGSBMI was significantly (multiple testing-corrected) associated with pancreatic cancer (OR[95%CI] = 1.0804[1.025-1.14], P = 0.0037). The significance of association declined after T2D adjustment (OR[95%CI] = 1.073[1.018-1.13], P = 0.00904). PGSWHRadjBMI association with pancreatic cancer was at the margin of statistical significance (OR[95%CI] = 1.047[0.99-1.104], P = 0.086). T2D adjustment effectively lost any suggestive association of PGSWHRadjBMI with pancreatic cancer (OR[95%CI] = 1.039[0.99-1.097], P = 0.14). MR analyses showed a nominally significant causal effect of WHRadjBMI on pancreatic cancer (OR[95%CI] = 1.00095[1.00011-1.0018], P = 0.027) but not for BMI on pancreatic cancer. Overall, we show that abdominal adiposity measured using WHRadjBMI, may be a more important causal risk factor for pancreatic cancer compared to total adiposity, with T2D being a potential driver of this relationship.

ß-Cell-Specific E2f1 Deficiency Impairs Glucose Homeostasis, ß-Cell Identity, and Insulin Secretion.

The loss of pancreatic ß-cell identity has emerged as an important feature of type 2 diabetes development, but the molecular mechanisms are still elusive. Here, we explore the cell-autonomous role of the cell-cycle regulator and transcription factor E2F1 in the maintenance of ß-cell identity, insulin secretion, and glucose homeostasis. We show that the ß-cell-specific loss of E2f1 function in mice triggers glucose intolerance associated with defective insulin secretion, altered endocrine cell mass, downregulation of many ß-cell genes, and concomitant increase of non-ß-cell markers. Mechanistically, epigenomic profiling of the promoters of these non-ß-cell upregulated genes identified an enrichment of bivalent H3K4me3/H3K27me3 or H3K27me3 marks. Conversely, promoters of downregulated genes were enriched in active chromatin H3K4me3 and H3K27ac histone marks. We find that specific E2f1 transcriptional, cistromic, and epigenomic signatures are associated with these ß-cell dysfunctions, with E2F1 directly regulating several ß-cell genes at the chromatin level. Finally, the pharmacological inhibition of E2F transcriptional activity in human islets also impairs insulin secretion and the expression of ß-cell identity genes. Our data suggest that E2F1 is critical for maintaining ß-cell identity and function through sustained control of ß-cell and non-ß-cell transcriptional programs. ARTICLE HIGHLIGHTS: ß-Cell-specific E2f1 deficiency in mice impairs glucose tolerance. Loss of E2f1 function alters the ratio of α- to ß-cells but does not trigger ß-cell conversion into α-cells. Pharmacological inhibition of E2F activity inhibits glucose-stimulated insulin secretion and alters ß- and α-cell gene expression in human islets. E2F1 maintains ß-cell function and identity through control of transcriptomic and epigenetic programs.

Heterozygous pathogenic variants in POMC are not responsible for monogenic obesity: Implication for MC4R agonist use.

PURPOSE: Recessive deficiency of proopiomelanocortin (POMC) causes childhood-onset severe obesity. Cases can now benefit from the melanocortin 4 receptor agonist setmelanotide. Furthermore, a phase 3 clinical trial is evaluating setmelanotide in heterozygotes for POMC. We performed a large-scale genetic analysis to assess the effect of heterozygous, pathogenic POMC variants on obesity. METHODS: A genetic analysis was performed in a family including 2 cousins with childhood-onset obesity. We analyzed the obesity status of heterozygotes for pathogenic POMC variants in the Human Gene Mutation Database. The association between heterozygous pathogenic POMC variants and obesity risk was assessed using 190,000 exome samples from UK Biobank. RESULTS: The 2 cousins carried a compound heterozygous pathogenic variant in POMC. Six siblings were heterozygotes; only 1 of them had obesity. In Human Gene Mutation Database, we identified 60 heterozygotes for pathogenic POMC variants, of whom 14 had obesity. In UK Biobank, heterozygous pathogenic POMC variants were not associated with obesity risk, but they modestly increased body mass index levels. CONCLUSION: Heterozygous pathogenic POMC variants do not contribute to monogenic obesity, but they slightly increase body mass index. Setmelanotide use in patients with obesity, which would only be based on the presence of a heterozygous POMC variant, can be questioned.

Biallelic Mutations in P4HTM Cause Syndromic Obesity.

We previously demonstrated that 50% of children with obesity from consanguineous families from Pakistan carry pathogenic variants in known monogenic obesity genes. Here, we have discovered a novel monogenetic recessive form of severe childhood obesity using an in-house computational staged approach. The analysis included whole-exome sequencing data of 366 children with severe obesity, 1,000 individuals of the Pakistan Risk of Myocardial Infarction Study (PROMIS) study, and 200,000 participants of the UK Biobank to prioritize genes harboring rare homozygous variants with putative effect on human obesity. We identified five rare or novel homozygous missense mutations predicted deleterious in five consanguineous families in P4HTM encoding prolyl 4-hydroxylase transmembrane (P4H-TM). We further found two additional homozygous missense mutations in children with severe obesity of Indian and Moroccan origin. Molecular dynamics simulation suggested that these mutations destabilized the active conformation of the substrate binding domain. Most carriers also presented with hypotonia, cognitive impairment, and/or developmental delay. Three of the five probands died of pneumonia during the first 2 years of the follow-up. P4HTM deficiency is a novel form of syndromic obesity, affecting 1.5% of our children with obesity associated with high mortality. P4H-TM is a hypoxia-inducible factor that is necessary for survival and adaptation under oxygen deprivation, but the role of this pathway in energy homeostasis and obesity pathophysiology remains to be elucidated.