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.

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.

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.

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.

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.

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.

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.

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.

A novel pathogenic variant in MRAP2 in an obese patient with successful outcome of bariatric surgery.

Mutations in genes encoding proteins located in the leptin/melanocortin pathway have been identified in the rare cases of genetic obesities. Heterozygous variants of MRAP2, encoding a G coupled-protein receptor accessory protein implicated in energy control notably via the melanocortin-4 receptor, have been recently identified. A 24-year-old patient with early-onset severe obesity (body mass index [BMI]: 64 kg/m2) associated with hypertension, respiratory complications, nonalcoholic fatty liver disease, and type 2 diabetes was referred to our department. Sleeve gastrectomy was successful. A new heterozygous variant in MRAP2 (NM_138409.4: c.154G>C/p.G52R) variant was identified in the patient DNA. Functional assessment confirmed that this new variant was pathogenic. We report a new pathogenic loss-of-function mutation in MRAP2 in a patient suffering from a severe multicomplicated obesity. This confirms the metabolic phenotype in patients with this monogenic form of obesity. Longer follow-up will be necessary. Our finding will allow a personalized medicine.

Amyloid-Beta Peptides Trigger Premature Functional and Gene Expression Alterations in Human-Induced Neurons.

Alzheimer's disease (AD) is the most prevalent cause of dementia in the elderly, characterized by the presence of amyloid-beta (Aß) plaques, neurofibrillary tangles, neuroinflammation, synapse loss and neurodegeneration in the brain. The amyloid cascade hypothesis postulates that deposition of Aß peptides is the causative agent of AD pathology, but we still lack comprehensive understanding of the molecular mechanisms connecting Aß peptides to neuronal dysfunctions in AD. In this work, we investigate the early effects of Aß peptide accumulation on the functional properties and gene expression profiles of human-induced neurons (hiNs). We show that hiNs acutely exposed to low concentrations of both cell-secreted Aß peptides or synthetic Aß1-42 exhibit alterations in the frequency of calcium transients suggestive of increased neuronal excitability. Using single-cell RNA sequencing, we also show that cell-secreted Aß up-regulates the expression of several synapse-related genes and down-regulates the expression of genes associated with metabolic stress mainly in glutamatergic neurons and, to a lesser degree, in GABAergic neurons and astrocytes. These neuronal alterations correlate with activation of the SEMA5, EPHA and NECTIN signaling pathways, which are important regulators of synaptic plasticity. Altogether, our findings indicate that slight elevations in Aß concentrations are sufficient to elicit transcriptional changes in human neurons, which can contribute to early alterations in neural network activity.

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.

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.

Heimdall, an alternative protein issued from a ncRNA related to kappa light chain variable region of immunoglobulins from astrocytes: a new player in neural proteome.

The dogma "One gene, one protein" is clearly obsolete since cells use alternative splicing and generate multiple transcripts which are translated into protein isoforms, but also use alternative translation initiation sites (TISs) and termination sites on a given transcript. Alternative open reading frames for individual transcripts give proteins originate from the 5'- and 3'-UTR mRNA regions, frameshifts of mRNA ORFs or from non-coding RNAs. Longtime considered as non-coding, recent in-silico translation prediction methods enriched the protein databases allowing the identification of new target structures that have not been identified previously. To gain insight into the role of these newly identified alternative proteins in the regulation of cellular functions, it is crucial to assess their dynamic modulation within a framework of altered physiological modifications such as experimental spinal cord injury (SCI). Here, we carried out a longitudinal proteomic study on rat SCI from 12 h to 10 days. Based on the alternative protein predictions, it was possible to identify a plethora of newly predicted protein hits. Among these proteins, some presented a special interest due to high homology with variable chain regions of immunoglobulins. We focus our interest on the one related to Kappa variable light chains which is similarly highly produced by B cells in the Bence jones disease, but here expressed in astrocytes. This protein, name Heimdall is an Intrinsically disordered protein which is secreted under inflammatory conditions. Immunoprecipitation experiments showed that the Heimdall interactome contained proteins related to astrocyte fate keepers such as "NOTCH1, EPHA3, IPO13" as well as membrane receptor protein including "CHRNA9; TGFBR, EPHB6, and TRAM". However, when Heimdall protein was neutralized utilizing a specific antibody or its gene knocked out by CRISPR-Cas9, sprouting elongations were observed in the corresponding astrocytes. Interestingly, depolarization assays and intracellular calcium measurements in Heimdall KO, established a depolarization effect on astrocyte membranes KO cells were more likely that the one found in neuroprogenitors. Proteomic analyses performed under injury conditions or under lipopolysaccharides (LPS) stimulation, revealed the expression of neuronal factors, stem cell proteins, proliferation, and neurogenesis of astrocyte convertor factors such as EPHA4, NOTCH2, SLIT3, SEMA3F, suggesting a role of Heimdall could regulate astrocytic fate. Taken together, Heimdall could be a novel member of the gatekeeping astrocyte-to-neuroprogenitor conversion factors.

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.