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.

The Alzheimer's disease risk gene BIN1 regulates activity-dependent gene expression in human-induced glutamatergic neurons.

Bridging Integrator 1 (BIN1) is the second most important Alzheimer's disease (AD) risk gene, but its physiological roles in neurons and its contribution to brain pathology remain largely elusive. In this work, we show that BIN1 plays a critical role in the regulation of calcium homeostasis, electrical activity, and gene expression of glutamatergic neurons. Using single-cell RNA-sequencing on cerebral organoids generated from isogenic BIN1 wild type (WT), heterozygous (HET) and homozygous knockout (KO) human-induced pluripotent stem cells (hiPSCs), we show that BIN1 is mainly expressed by oligodendrocytes and glutamatergic neurons, like in the human brain. Both BIN1 HET and KO cerebral organoids show specific transcriptional alterations, mainly associated with ion transport and synapses in glutamatergic neurons. We then demonstrate that BIN1 cell-autonomously regulates gene expression in glutamatergic neurons by using a novel protocol to generate pure culture of hiPSC-derived induced neurons (hiNs). Using this system, we also show that BIN1 plays a key role in the regulation of neuronal calcium transients and electrical activity via its interaction with the L-type voltage-gated calcium channel Cav1.2. BIN1 KO hiNs show reduced activity-dependent internalization and higher Cav1.2 expression compared to WT hiNs. Pharmacological blocking of this channel with clinically relevant doses of nifedipine, a calcium channel blocker, partly rescues electrical and gene expression alterations in BIN1 KO glutamatergic neurons. Further, we show that transcriptional alterations in BIN1 KO hiNs that affect biological processes related to calcium homeostasis are also present in glutamatergic neurons of the human brain at late stages of AD pathology. Together, these findings suggest that BIN1-dependent alterations in neuronal properties could contribute to AD pathophysiology and that treatment with low doses of clinically approved calcium blockers should be considered as an option to slow disease-onset and progression.

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 polygenic risk scores 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 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) and several cancers (breast cancer; 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.

Towards the recognition of oligogenic forms of type 2 diabetes.

The demarcation between monogenic and polygenic type 2 diabetes (T2D) is less distinct than previously believed. Notably, recent research has highlighted a new entity, that we suggest calling oligogenic forms of T2D, serving as a genetic link between these two forms. In this opinion article, we have reviewed scientific advances that suggest categorizing genes involved in oligogenic T2D. Research focused on polygenic T2D has faced challenges in deepening our comprehension of the pathophysiology of T2D due to the inability to directly establish causal links between a signal and the molecular mechanisms underlying the disease. However, the study of oligogenic forms of T2D has illuminated distinct causal connections between genes and disease risk, thereby indicating potential new drug targets.

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.

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.

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.

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.

PNLIPRP1 hypermethylation in exocrine pancreas links type 2 diabetes and cholesterol metabolism.

We postulated that T2D predisposes to exocrine pancreatic diseases through (epi)genetic mechanisms. We explored the methylome (methylationEPIC arrays) of the exocrine pancreas of 141 donors, assessing the impact of T2D. Epigenome-wide association study (EWAS) for T2D identified a hypermethylation in an enhancer of the Pancreatic-Lipase-Related-Protein 1 (PNLIPRP1) gene, associated with decreased PNLIPRP1 expression. PNLIPRP1 null variants (in 191K participants of the UKbiobank) associated with elevated glycemia and LDL-cholesterol. Mendelian Randomisation using 2.5M SNP OmniArrays in 111 donors evidenced that T2D was causal of PNLIPRP1 hypermethylation, which was causal for LDL-cholesterol. Further AR42J rat exocrine cell studies demonstrated that Pnliprp1 knockdown induced acinar-to-ductal metaplasia, a known pre-pancreatic cancer state, and increased cholesterol levels, reversible with statin. This (epi)genetic study suggests a role for PNLIPRP1 in human metabolism and on exocrine pancreas function with potential implications for pancreatic diseases.

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.

Functional genetics reveals the contribution of delta opioid receptor to type 2 diabetes and beta-cell function.

Functional genetics has identified drug targets for metabolic disorders. Opioid use impacts metabolic homeostasis, although mechanisms remain elusive. Here, we explore the OPRD1 gene (encoding delta opioid receptor, DOP) to understand its impact on type 2 diabetes. Large-scale sequencing of OPRD1 and in vitro analysis reveal that loss-of-function variants are associated with higher adiposity and lower hyperglycemia risk, whereas gain-of-function variants are associated with lower adiposity and higher type 2 diabetes risk. These findings align with studies of opium addicts. OPRD1 is expressed in human islets and beta cells, with decreased expression under type 2 diabetes conditions. DOP inhibition by an antagonist enhances insulin secretion from human beta cells and islets. RNA-sequencing identifies pathways regulated by DOP antagonism, including nerve growth factor, circadian clock, and nuclear receptor pathways. Our study highlights DOP as a key player between opioids and metabolic homeostasis, suggesting its potential as a therapeutic target for type 2 diabetes.