Genetic variant effects on gene expression in human pancreatic islets and their implications for T2D

Untangling the genetics of beta cell dysfunction and death in type 1 diabetes

BACKGROUND

Type 1 diabetes (T1D) is a complex multi-system disease which arises from both environmental and genetic factors, resulting in the destruction of insulin-producing pancreatic beta cells. Over the past two decades, human genetic studies have provided new insight into the etiology of T1D, including an appreciation for the role of beta cells in their own demise.

SCOPE OF REVIEW

Here, we outline models supported by human genetic data for the role of beta cell dysfunction and death in T1D. We highlight the importance of strong evidence linking T1D genetic associations to bona fide candidate genes for mechanistic and therapeutic consideration. To guide rigorous interpretation of genetic associations, we describe molecular profiling approaches, genomic resources, and disease models that may be used to construct variant-to-gene links and to investigate candidate genes and their role in T1D.

MAJOR CONCLUSIONS

We profile advances in understanding the genetic causes of beta cell dysfunction and death at individual T1D risk loci. We discuss how genetic risk prediction models can be used to address disease heterogeneity. Further, we present areas where investment will be critical for the future use of genetics to address open questions in the development of new treatment and prevention strategies for T1D.

Multi-omics characterization of type 2 diabetes associated genetic variation

Discerning the mechanisms driving type 2 diabetes (T2D) pathophysiology from genome-wide association studies (GWAS) remains a challenge. To this end, we integrated omics information from 16 multi-tissue and multi-ancestry expression, protein, and metabolite quantitative trait loci (QTL) studies and 46 multi-ancestry GWAS for T2D-related traits with the largest, most ancestrally diverse T2D GWAS to date. Of the 1,289 T2D GWAS index variants, 716 (56%) demonstrated strong evidence of colocalization with a molecular or T2D-related trait, implicating 657 cis-effector genes, 1,691 distal-effector genes, 731 metabolites, and 43 T2D-related traits. We identified 773 of these cis- and distal-effector genes using either expression QTL data from understudied ancestry groups or inclusion of T2D index variants enriched in underrepresented populations, emphasizing the value of increasing population diversity in functional mapping. Linking these variants, genes, metabolites, and traits into a network, we elucidated mechanisms through which T2D-associated variation may impact disease risk. Finally, we showed that drugs targeting effector proteins were enriched in those approved to treat T2D, highlighting the potential of these results to prioritize drug targets for T2D. These results represent a leap in the molecular characterization of T2D-associated genetic variation and will aid in translating genetic findings into novel therapeutic strategies.

The Apo gene's genetic variants: hidden role in Asian vascular risk

Vascular risk factors, including diabetes, hypertension, hyperlipidemia, and obesity, pose significant health threats with implications extending to neuropsychiatric disorders such as stroke and Alzheimer's disease. The Asian population, in particular, appears to be disproportionately affected due to unique genetic predispositions, as well as epigenetic factors such as dietary patterns and lifestyle habits. Existing management strategies often fall short of addressing these specific needs, leading to greater challenges in prevention and treatment. This review highlights a significant gap in our understanding of the impact of genetic screening in the early detection and tailored treatment of vascular risk factors among the Asian population. Apolipoprotein, a key player in cholesterol metabolism, is primarily associated with dyslipidemia, yet emerging evidence suggests its involvement in conditions such as diabetes, hypertension, and obesity. While genetic variants of vascular risk are ethnic-dependent, current evidence indicates that epigenetics also exhibits ethnic specificity. Understanding the interplay between Apolipoprotein and genetics, particularly within diverse ethnic backgrounds, has the potential to refine risk stratification and enhance precision in management. For Caucasian carrying the APOA5 rs662799 C variant, pharmacological interventions are recommended, as dietary interventions may not be sufficient. In contrast, for Asian populations with the same genetic variant, dietary modifications are initially advised. Should dyslipidemia persist, the consideration of pharmaceutical agents such as statins is recommended.

HumanIslets: An integrated platform for human islet data access and analysis

Comprehensive molecular and cellular phenotyping of human islets can enable deep mechanistic insights for diabetes research. We established the Human Islet Data Analysis and Sharing (HI-DAS) consortium to advance goals in accessibility, usability, and integration of data from human islets isolated from donors with and without diabetes at the Alberta Diabetes Institute (ADI) IsletCore. Here we introduce HumanIslets.com, an open resource for the research community. This platform, which presently includes data on 547 human islet donors, allows users to access linked datasets describing molecular profiles, islet function and donor phenotypes, and to perform various statistical and functional analyses at the donor, islet and single-cell levels. As an example of the analytic capacity of this resource we show a dissociation between cell culture effects on transcript and protein expression, and an approach to correct for exocrine contamination found in hand-picked islets. Finally, we provide an example workflow and visualization that highlights links between type 2 diabetes status, SERCA3b Ca2+-ATPase levels at the transcript and protein level, insulin secretion and islet cell phenotypes. HumanIslets.com provides a growing and adaptable set of resources and tools to support the metabolism and diabetes research community.

From whole bodies to single cells: A guide to transcriptomic approaches for ecology and evolutionary biology

RNA sequencing (RNAseq) methodology has experienced a burst of technological developments in the last decade, which has opened up opportunities for studying the mechanisms of adaptation to environmental factors at both the organismal and cellular level. Selecting the most suitable experimental approach for specific research questions and model systems can, however, be a challenge and researchers in ecology and evolution are commonly faced with the choice of whether to study gene expression variation in whole bodies, specific tissues, and/or single cells. A wide range of sometimes polarised opinions exists over which approach is best. Here, we highlight the advantages and disadvantages of each of these approaches to provide a guide to help researchers make informed decisions and maximise the power of their study. Using illustrative examples of various ecological and evolutionary research questions, we guide the readers through the different RNAseq approaches and help them identify the most suitable design for their own projects.

Unique features of β-cell metabolism are lost in type 2 diabetes

Pancreatic β cells play an essential role in the control of systemic glucose homeostasis as they sense blood glucose levels and respond by secreting insulin. Upon stimulating glucose uptake in insulin-sensitive tissues post-prandially, this anabolic hormone restores blood glucose levels to pre-prandial levels. Maintaining physiological glucose levels thus relies on proper β-cell function. To fulfill this highly specialized nutrient sensor role, β cells have evolved a unique genetic program that shapes its distinct cellular metabolism. In this review, the unique genetic and metabolic features of β cells will be outlined, including their alterations in type 2 diabetes (T2D). β cells selectively express a set of genes in a cell type-specific manner; for instance, the glucose activating hexokinase IV enzyme or Glucokinase (GCK), whereas other genes are selectively "disallowed", including lactate dehydrogenase A (LDHA) and monocarboxylate transporter 1 (MCT1). This selective gene program equips β cells with a unique metabolic apparatus to ensure that nutrient metabolism is coupled to appropriate insulin secretion, thereby avoiding hyperglycemia, as well as life-threatening hypoglycemia. Unlike most cell types, β cells exhibit specialized bioenergetic features, including supply-driven rather than demand-driven metabolism and a high basal mitochondrial proton leak respiration. The understanding of these unique genetically programmed metabolic features and their alterations that lead to β-cell dysfunction is crucial for a comprehensive understanding of T2D pathophysiology and the development of innovative therapeutic approaches for T2D patients.

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.

Multi-omic QTL mapping in early developmental tissues reveals phenotypic and temporal complexity of regulatory variants underlying GWAS loci

Most GWAS loci are presumed to affect gene regulation, however, only ∼43% colocalize with expression quantitative trait loci (eQTLs). To address this colocalization gap, we identify eQTLs, chromatin accessibility QTLs (caQTLs), and histone acetylation QTLs (haQTLs) using molecular samples from three early developmental (EDev) tissues. Through colocalization, we annotate 586 GWAS loci for 17 traits by QTL complexity, QTL phenotype, and QTL temporal specificity. We show that GWAS loci are highly enriched for colocalization with complex QTL modules that affect multiple elements (genes and/or peaks). We also demonstrate that caQTLs and haQTLs capture regulatory variations not associated with eQTLs and explain ∼49% of the functionally annotated GWAS loci. Additionally, we show that EDev-unique QTLs are strongly depleted for colocalizing with GWAS loci. By conducting one of the largest multi-omic QTL studies to date, we demonstrate that many GWAS loci exhibit phenotypic complexity and therefore, are missed by traditional eQTL analyses.

Proinsulin folding and trafficking defects trigger a common pathological disturbance of endoplasmic reticulum homeostasis

Primary defects in folding of mutant proinsulin can cause dominant-negative proinsulin accumulation in the endoplasmic reticulum (ER), impaired anterograde proinsulin trafficking, perturbed ER homeostasis, diminished insulin production, and β-cell dysfunction. Conversely, if primary impairment of ER-to-Golgi trafficking (which also perturbs ER homeostasis) drives misfolding of nonmutant proinsulin-this might suggest bi-directional entry into a common pathological phenotype (proinsulin misfolding, perturbed ER homeostasis, and deficient ER export of proinsulin) that can culminate in diminished insulin storage and diabetes. Here, we've challenged β-cells with conditions that impair ER-to-Golgi trafficking, and devised an accurate means to assess the relative abundance of distinct folded/misfolded forms of proinsulin using a novel nonreducing SDS-PAGE/immunoblotting protocol. We confirm abundant proinsulin misfolding upon introduction of a diabetogenic INS mutation, or in the islets of db/db mice. Whereas blockade of proinsulin trafficking in Golgi/post-Golgi compartments results in intracellular accumulation of properly-folded proinsulin (bearing native disulfide bonds), impairment of ER-to-Golgi trafficking (regardless whether such impairment is achieved by genetic or pharmacologic means) results in decreased native proinsulin with more misfolded proinsulin. Remarkably, reversible ER-to-Golgi transport defects (such as treatment with brefeldin A or cellular energy depletion) upon reversal quickly restore the ER folding environment, resulting in the disappearance of pre-existing misfolded proinsulin while preserving proinsulin bearing native disulfide bonds. Thus, proper homeostatic balance of ER-to-Golgi trafficking is linked to a more favorable proinsulin folding (as well as trafficking) outcome.

Transcriptome-wide association studies associated with Crohn's disease: challenges and perspectives

Crohn's disease (CD) is regarded as a lifelong progressive disease affecting all segments of the intestinal tract and multiple organs. Based on genome-wide association studies (GWAS) and gene expression data, transcriptome-wide association studies (TWAS) can help identify susceptibility genes associated with pathogenesis and disease behavior. In this review, we overview seven reported TWASs of CD, summarize their study designs, and discuss the key methods and steps used in TWAS, which affect the prioritization of susceptibility genes. This article summarized the screening of tissue-specific susceptibility genes for CD, and discussed the reported potential pathological mechanisms of overlapping susceptibility genes related to CD in a certain tissue type. We observed that ileal lipid-related metabolism and colonic extracellular vesicles may be involved in the pathogenesis of CD by performing GO pathway enrichment analysis for susceptibility genes. We further pointed the low reproducibility of TWAS associated with CD and discussed the reasons for these issues, strategies for solving them. In the future, more TWAS are needed to be designed into large-scale, unified cohorts, unified analysis pipelines, and fully classified databases of expression trait loci.

Genomic insights into the comorbidity between type 2 diabetes and schizophrenia

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

Multiplexed CRISPR gene editing in primary human islet cells with Cas9 ribonucleoprotein

Successful genome editing in primary human islets could reveal features of the genetic regulatory landscape underlying β cell function and diabetes risk. Here, we describe a CRISPR-based strategy to interrogate functions of predicted regulatory DNA elements using electroporation of a complex of Cas9 ribonucleoprotein (Cas9 RNP) and guide RNAs into primary human islet cells. We successfully targeted coding regions including the PDX1 exon 1, and non-coding DNA linked to diabetes susceptibility. CRISPR-Cas9 RNP approaches revealed genetic targets of regulation by DNA elements containing candidate diabetes risk SNPs, including an in vivo enhancer of the MPHOSPH9 gene. CRISPR-Cas9 RNP multiplexed targeting of two cis-regulatory elements linked to diabetes risk in PCSK1, which encodes an endoprotease crucial for Insulin processing, also demonstrated efficient simultaneous editing of PCSK1 regulatory elements, resulting in impaired β cell PCSK1 regulation and Insulin secretion. Multiplex CRISPR-Cas9 RNP provides powerful approaches to investigate and elucidate human islet cell gene regulation in health and diabetes.

COLOCdb: a comprehensive resource for multi-model colocalization of complex traits

Large-scale genome-wide association studies (GWAS) have provided profound insights into complex traits and diseases. Yet, deciphering the fine-scale molecular mechanisms of how genetic variants manifest to cause the phenotypes remains a daunting task. Here, we present COLOCdb (https://ngdc.cncb.ac.cn/colocdb), a comprehensive genetic colocalization database by integrating more than 3000 GWAS summary statistics and 13 types of xQTL to date. By employing two representative approaches for the colocalization analysis, COLOCdb deposits results from three key components: (i) GWAS-xQTL, pair-wise colocalization between GWAS loci and different types of xQTL, (ii) GWAS-GWAS, pair-wise colocalization between the trait-associated genetic loci from GWASs and (iii) xQTL-xQTL, pair-wise colocalization between the genetic loci associated with molecular phenotypes in xQTLs. These results together represent the most comprehensive colocalization analysis, which also greatly expands the list of shared variants with genetic pleiotropy. We expect that COLOCdb can serve as a unique and useful resource in advancing the discovery of new biological mechanisms and benefit future functional studies.

Genetic risk converges on regulatory networks mediating early type 2 diabetes

Type 2 diabetes mellitus (T2D), a major cause of worldwide morbidity and mortality, is characterized by dysfunction of insulin-producing pancreatic islet β cells1,2. T2D genome-wide association studies (GWAS) have identified hundreds of signals in non-coding and β cell regulatory genomic regions, but deciphering their biological mechanisms remains challenging3-5. Here, to identify early disease-driving events, we performed traditional and multiplexed pancreatic tissue imaging, sorted-islet cell transcriptomics and islet functional analysis of early-stage T2D and control donors. By integrating diverse modalities, we show that early-stage T2D is characterized by β cell-intrinsic defects that can be proportioned into gene regulatory modules with enrichment in signals of genetic risk. After identifying the β cell hub gene and transcription factor RFX6 within one such module, we demonstrated multiple layers of genetic risk that converge on an RFX6-mediated network to reduce insulin secretion by β cells. RFX6 perturbation in primary human islet cells alters β cell chromatin architecture at regions enriched for T2D GWAS signals, and population-scale genetic analyses causally link genetically predicted reduced RFX6 expression with increased T2D risk. Understanding the molecular mechanisms of complex, systemic diseases necessitates integration of signals from multiple molecules, cells, organs and individuals, and thus we anticipate that this approach will be a useful template to identify and validate key regulatory networks and master hub genes for other diseases or traits using GWAS data.

Functional interrogation of twenty type 2 diabetes-associated genes using isogenic human embryonic stem cell-derived β-like cells

Genetic studies have identified numerous loci associated with type 2 diabetes (T2D), but the functional roles of many loci remain unexplored. Here, we engineered isogenic knockout human embryonic stem cell lines for 20 genes associated with T2D risk. We examined the impacts of each knockout on β cell differentiation, functions, and survival. We generated gene expression and chromatin accessibility profiles on β cells derived from each knockout line. Analyses of T2D-association signals overlapping HNF4A-dependent ATAC peaks identified a likely causal variant at the FAIM2 T2D-association signal. Additionally, the integrative association analyses identified four genes (CP, RNASE1, PCSK1N, and GSTA2) associated with insulin production, and two genes (TAGLN3 and DHRS2) associated with β cell sensitivity to lipotoxicity. Finally, we leveraged deep ATAC-seq read coverage to assess allele-specific imbalance at variants heterozygous in the parental line and identified a single likely functional variant at each of 23 T2D-association signals.

eQTL mapping in fetal-like pancreatic progenitor cells reveals early developmental insights into diabetes risk

The impact of genetic regulatory variation active in early pancreatic development on adult pancreatic disease and traits is not well understood. Here, we generate a panel of 107 fetal-like iPSC-derived pancreatic progenitor cells (iPSC-PPCs) from whole genome-sequenced individuals and identify 4065 genes and 4016 isoforms whose expression and/or alternative splicing are affected by regulatory variation. We integrate eQTLs identified in adult islets and whole pancreas samples, which reveal 1805 eQTL associations that are unique to the fetal-like iPSC-PPCs and 1043 eQTLs that exhibit regulatory plasticity across the fetal-like and adult pancreas tissues. Colocalization with GWAS risk loci for pancreatic diseases and traits show that some putative causal regulatory variants are active only in the fetal-like iPSC-PPCs and likely influence disease by modulating expression of disease-associated genes in early development, while others with regulatory plasticity likely exert their effects in both the fetal and adult pancreas by modulating expression of different disease genes in the two developmental stages.

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.

Guidelines and Recommendations for Laboratory Analysis in the Diagnosis and Management of Diabetes Mellitus

BACKGROUND

Numerous laboratory tests are used in the diagnosis and management of diabetes mellitus. The quality of the scientific evidence supporting the use of these assays varies substantially.

APPROACH

An expert committee compiled evidence-based recommendations for laboratory analysis in screening, diagnosis, or monitoring of diabetes. The overall quality of the evidence and the strength of the recommendations were evaluated. The draft consensus recommendations were evaluated by invited reviewers and presented for public comment. Suggestions were incorporated as deemed appropriate by the authors (see Acknowledgments). The guidelines were reviewed by the Evidence Based Laboratory Medicine Committee and the Board of Directors of the American Association for Clinical Chemistry and by the Professional Practice Committee of the American Diabetes Association.

CONTENT

Diabetes can be diagnosed by demonstrating increased concentrations of glucose in venous plasma or increased hemoglobin A1c (HbA1c) in the blood. Glycemic control is monitored by the people with diabetes measuring their own blood glucose with meters and/or with continuous interstitial glucose monitoring (CGM) devices and also by laboratory analysis of HbA1c. The potential roles of noninvasive glucose monitoring, genetic testing, and measurement of ketones, autoantibodies, urine albumin, insulin, proinsulin, and C-peptide are addressed.

SUMMARY

The guidelines provide specific recommendations based on published data or derived from expert consensus. Several analytes are found to have minimal clinical value at the present time, and measurement of them is not recommended.

Multiplexed CRISPR gene editing in primary human islet cells with Cas9 ribonucleoprotein

Successful genome editing in primary human islets could reveal features of the genetic regulatory landscape underlying β cell function and diabetes risk. Here, we describe a CRISPR-based strategy to interrogate functions of predicted regulatory DNA elements using electroporation of a complex of Cas9 ribonucleoprotein (Cas9 RNP) and guide RNAs into primary human islet cells. We successfully targeted coding regions including the PDX1 exon 1, and non-coding DNA linked to diabetes susceptibility. CRISPR/Cas9 RNP approaches revealed genetic targets of regulation by DNA elements containing candidate diabetes risk SNPs, including an in vivo enhancer of the MPHOSPH9 gene. CRISPR/Cas9 RNP multiplexed targeting of two cis-regulatory elements linked to diabetes risk in PCSK1, which encodes an endoprotease crucial for insulin processing, also demonstrated efficient simultaneous editing of PCSK1 regulatory elements, resulting in impaired β cell PCSK1 regulation and insulin secretion. Multiplex CRISPR/Cas9 RNP provides powerful approaches to investigate and elucidate human islet cell gene regulation in health and diabetes.

Delineating mouse β-cell identity during lifetime and in diabetes with a single cell atlas

Although multiple pancreatic islet single-cell RNA-sequencing (scRNA-seq) datasets have been generated, a consensus on pancreatic cell states in development, homeostasis and diabetes as well as the value of preclinical animal models is missing. Here, we present an scRNA-seq cross-condition mouse islet atlas (MIA), a curated resource for interactive exploration and computational querying. We integrate over 300,000 cells from nine scRNA-seq datasets consisting of 56 samples, varying in age, sex and diabetes models, including an autoimmune type 1 diabetes model (NOD), a glucotoxicity/lipotoxicity type 2 diabetes model (db/db) and a chemical streptozotocin β-cell ablation model. The β-cell landscape of MIA reveals new cell states during disease progression and cross-publication differences between previously suggested marker genes. We show that β-cells in the streptozotocin model transcriptionally correlate with those in human type 2 diabetes and mouse db/db models, but are less similar to human type 1 diabetes and mouse NOD β-cells. We also report pathways that are shared between β-cells in immature, aged and diabetes models. MIA enables a comprehensive analysis of β-cell responses to different stressors, providing a roadmap for the understanding of β-cell plasticity, compensation and demise.

Modeling islet enhancers using deep learning identifies candidate causal variants at loci associated with T2D and glycemic traits

Genetic association studies have identified hundreds of independent signals associated with type 2 diabetes (T2D) and related traits. Despite these successes, the identification of specific causal variants underlying a genetic association signal remains challenging. In this study, we describe a deep learning (DL) method to analyze the impact of sequence variants on enhancers. Focusing on pancreatic islets, a T2D relevant tissue, we show that our model learns islet-specific transcription factor (TF) regulatory patterns and can be used to prioritize candidate causal variants. At 101 genetic signals associated with T2D and related glycemic traits where multiple variants occur in linkage disequilibrium, our method nominates a single causal variant for each association signal, including three variants previously shown to alter reporter activity in islet-relevant cell types. For another signal associated with blood glucose levels, we biochemically test all candidate causal variants from statistical fine-mapping using a pancreatic islet beta cell line and show biochemical evidence of allelic effects on TF binding for the model-prioritized variant. To aid in future research, we publicly distribute our model and islet enhancer perturbation scores across ~67 million genetic variants. We anticipate that DL methods like the one presented in this study will enhance the prioritization of candidate causal variants for functional studies.

Genetic analysis of blood molecular phenotypes reveals common properties in the regulatory networks affecting complex traits

We evaluate the shared genetic regulation of mRNA molecules, proteins and metabolites derived from whole blood from 3029 human donors. We find abundant allelic heterogeneity, where multiple variants regulate a particular molecular phenotype, and pleiotropy, where a single variant associates with multiple molecular phenotypes over multiple genomic regions. The highest proportion of share genetic regulation is detected between gene expression and proteins (66.6%), with a further median shared genetic associations across 49 different tissues of 78.3% and 62.4% between plasma proteins and gene expression. We represent the genetic and molecular associations in networks including 2828 known GWAS variants, showing that GWAS variants are more often connected to gene expression in trans than other molecular phenotypes in the network. Our work provides a roadmap to understanding molecular networks and deriving the underlying mechanism of action of GWAS variants using different molecular phenotypes in an accessible tissue.

Inferring causal genes at type 2 diabetes GWAS loci through chromosome interactions in islet cells

Background: Resolving causal genes for type 2 diabetes at loci implicated by genome-wide association studies (GWAS) requires integrating functional genomic data from relevant cell types. Chromatin features in endocrine cells of the pancreatic islet are particularly informative and recent studies leveraging chromosome conformation capture (3C) with Hi-C based methods have elucidated regulatory mechanisms in human islets. However, these genome-wide approaches are less sensitive and afford lower resolution than methods that target specific loci. Methods: To gauge the extent to which targeted 3C further resolves chromatin-mediated regulatory mechanisms at GWAS loci, we generated interaction profiles at 23 loci using next-generation (NG) capture-C in a human beta cell model (EndoC-βH1) and contrasted these maps with Hi-C maps in EndoC-βH1 cells and human islets and a promoter capture Hi-C map in human islets. Results: We found improvements in assay sensitivity of up to 33-fold and resolved ~3.6X more chromatin interactions. At a subset of 18 loci with 25 co-localised GWAS and eQTL signals, NG Capture-C interactions implicated effector transcripts at five additional genetic signals relative to promoter capture Hi-C through physical contact with gene promoters. Conclusions: High resolution chromatin interaction profiles at selectively targeted loci can complement genome- and promoter-wide maps.

Identifying novel regulatory effects for clinically relevant genes through the study of the Greek population

BACKGROUND

Expression quantitative trait loci (eQTL) studies provide insights into regulatory mechanisms underlying disease risk. Expanding studies of gene regulation to underexplored populations and to medically relevant tissues offers potential to reveal yet unknown regulatory variants and to better understand disease mechanisms. Here, we performed eQTL mapping in subcutaneous (S) and visceral (V) adipose tissue from 106 Greek individuals (Greek Metabolic study, GM) and compared our findings to those from the Genotype-Tissue Expression (GTEx) resource.

RESULTS

We identified 1,930 and 1,515 eGenes in S and V respectively, over 13% of which are not observed in GTEx adipose tissue, and that do not arise due to different ancestry. We report additional context-specific regulatory effects in genes of clinical interest (e.g. oncogene ST7) and in genes regulating responses to environmental stimuli (e.g. MIR21, SNX33). We suggest that a fraction of the reported differences across populations is due to environmental effects on gene expression, driving context-specific eQTLs, and suggest that environmental effects can determine the penetrance of disease variants thus shaping disease risk. We report that over half of GM eQTLs colocalize with GWAS SNPs and of these colocalizations 41% are not detected in GTEx. We also highlight the clinical relevance of S adipose tissue by revealing that inflammatory processes are upregulated in individuals with obesity, not only in V, but also in S tissue.

CONCLUSIONS

By focusing on an understudied population, our results provide further candidate genes for investigation regarding their role in adipose tissue biology and their contribution to disease risk and pathogenesis.

Genetic underpinning of the comorbidity between type 2 diabetes and osteoarthritis

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

Guidelines and Recommendations for Laboratory Analysis in the Diagnosis and Management of Diabetes Mellitus

BACKGROUND

Numerous laboratory tests are used in the diagnosis and management of diabetes mellitus. The quality of the scientific evidence supporting the use of these assays varies substantially.

APPROACH

An expert committee compiled evidence-based recommendations for laboratory analysis in screening, diagnosis, or monitoring of diabetes. The overall quality of the evidence and the strength of the recommendations were evaluated. The draft consensus recommendations were evaluated by invited reviewers and presented for public comment. Suggestions were incorporated as deemed appropriate by the authors (see Acknowledgments). The guidelines were reviewed by the Evidence Based Laboratory Medicine Committee and the Board of Directors of the American Association of Clinical Chemistry and by the Professional Practice Committee of the American Diabetes Association.

CONTENT

Diabetes can be diagnosed by demonstrating increased concentrations of glucose in venous plasma or increased hemoglobin A1c (Hb A1c) in the blood. Glycemic control is monitored by the people with diabetes measuring their own blood glucose with meters and/or with continuous interstitial glucose monitoring (CGM) devices and also by laboratory analysis of Hb A1c. The potential roles of noninvasive glucose monitoring, genetic testing, and measurement of ketones, autoantibodies, urine albumin, insulin, proinsulin, and C-peptide are addressed.

SUMMARY

The guidelines provide specific recommendations based on published data or derived from expert consensus. Several analytes are found to have minimal clinical value at the present time, and measurement of them is not recommended.

Transcriptome- and proteome-wide association studies nominate determinants of kidney function and damage

BACKGROUND

The pathophysiological causes of kidney disease are not fully understood. Here we show that the integration of genome-wide genetic, transcriptomic, and proteomic association studies can nominate causal determinants of kidney function and damage.

RESULTS

Through transcriptome-wide association studies (TWAS) in kidney cortex, kidney tubule, liver, and whole blood and proteome-wide association studies (PWAS) in plasma, we assess for effects of 12,893 genes and 1342 proteins on kidney filtration (glomerular filtration rate (GFR) estimated by creatinine; GFR estimated by cystatin C; and blood urea nitrogen) and kidney damage (albuminuria). We find 1561 associations distributed among 260 genomic regions that are supported as putatively causal. We then prioritize 153 of these genomic regions using additional colocalization analyses. Our genome-wide findings are supported by existing knowledge (animal models for MANBA, DACH1, SH3YL1, INHBB), exceed the underlying GWAS signals (28 region-trait combinations without significant GWAS hit), identify independent gene/protein-trait associations within the same genomic region (INHBC, SPRYD4), nominate tissues underlying the associations (tubule expression of NRBP1), and distinguish markers of kidney filtration from those with a role in creatinine and cystatin C metabolism. Furthermore, we follow up on members of the TGF-beta superfamily of proteins and find a prognostic value of INHBC for kidney disease progression even after adjustment for measured glomerular filtration rate (GFR).

CONCLUSION

In summary, this study combines multimodal, genome-wide association studies to generate a catalog of putatively causal target genes and proteins relevant to kidney function and damage which can guide follow-up studies in physiology, basic science, and clinical medicine.

Body Fat Distribution, Fasting Insulin Levels, and Insulin Secretion: A Bidirectional Mendelian Randomization Study

CONTEXT

Hyperinsulinemia and adiposity are associated with one another, but the directionality of this relation is debated.

OBJECTIVE

Here, we tested the direction of the causal effects of fasting insulin (FI) levels and body fat accumulation/distribution using 2-sample bidirectional Mendelian randomization (MR).

METHODS

We included summary statistics from large-scale genome-wide association studies for body mass index (BMI, n = 806 834), waist to hip ratio adjusted for BMI (WHRadjBMI, n = 694 649), abdominal subcutaneous, visceral and gluteofemoral adipose tissue (n = 38 965), FI levels (n = 98 210), pancreatic islets gene expression (n = 420), and hypothalamus gene expression (n = 155). We used inverse variance-weighted and robust MR methods that relied on statistically and biologically driven genetic instruments.

RESULTS

Both BMI and WHRadjBMI were positively associated with FI. Results were consistent across all robust MR methods and when variants mapped to the hypothalamus (presumably associated with food behavior) were included. In multivariable MR analyses, when waist circumference and BMI were mutually adjusted, the direct effect of waist circumference on FI was 2.43 times larger than the effect of BMI on FI. FI was not associated with adiposity. By contrast, using genetic instruments mapped to gene expression in pancreatic islets (presumably more specific to insulin secretion), insulin was positively associated with BMI and abdominal subcutaneous and gluteofemoral adipose tissue, but not with visceral adipose tissue.

CONCLUSION

Although these results will need to be supported by experimental investigations, results of this MR study suggest that abdominal adiposity may be a key determinant of circulating insulin levels. Alternatively, insulin secretion may promote peripheral adipose tissue accumulation.

Functional interrogation of twenty type 2 diabetes-associated genes using isogenic hESC-derived β-like cells

Genetic studies have identified numerous loci associated with type 2 diabetes (T2D), but the functional role of many loci has remained unexplored. In this study, we engineered isogenic knockout human embryonic stem cell (hESC) lines for 20 genes associated with T2D risk. We systematically examined β-cell differentiation, insulin production and secretion, and survival. We performed RNA-seq and ATAC-seq on hESC-β cells from each knockout line. Analyses of T2D GWAS signals overlapping with HNF4A-dependent ATAC peaks identified a specific SNP as a likely causal variant. In addition, we performed integrative association analyses and identified four genes ( CP, RNASE1, PCSK1N and GSTA2 ) associated with insulin production, and two genes ( TAGLN3 and DHRS2 ) associated with sensitivity to lipotoxicity. Finally, we leveraged deep ATAC-seq read coverage to assess allele-specific imbalance at variants heterozygous in the parental hESC line, to identify a single likely functional variant at each of 23 T2D GWAS signals.

Physiological and pathological roles of lipogenesis

Lipids are essential metabolites, which function as energy sources, structural components and signalling mediators. Most cells are able to convert carbohydrates into fatty acids, which are often converted into neutral lipids for storage in the form of lipid droplets. Accumulating evidence suggests that lipogenesis plays a crucial role not only in metabolic tissues for systemic energy homoeostasis but also in immune and nervous systems for their proliferation, differentiation and even pathophysiological roles. Thus, excessive or insufficient lipogenesis is closely associated with aberrations in lipid homoeostasis, potentially leading to pathological consequences, such as dyslipidaemia, diabetes, fatty liver, autoimmune diseases, neurodegenerative diseases and cancers. For systemic energy homoeostasis, multiple enzymes involved in lipogenesis are tightly controlled by transcriptional and post-translational modifications. In this Review, we discuss recent findings regarding the regulatory mechanisms, physiological roles and pathological importance of lipogenesis in multiple tissues such as adipose tissue and the liver, as well as the immune and nervous systems. Furthermore, we briefly introduce the therapeutic implications of lipogenesis modulation.

Effect of tissue-grouped regulatory variants associated to type 2 diabetes in related secondary outcomes

Genome-wide association studies have identified over five hundred loci that contribute to variation in type 2 diabetes (T2D), an established risk factor for many diseases. However, the mechanisms and extent through which these loci contribute to subsequent outcomes remain elusive. We hypothesized that combinations of T2D-associated variants acting on tissue-specific regulatory elements might account for greater risk for tissue-specific outcomes, leading to diversity in T2D disease progression. We searched for T2D-associated variants acting on regulatory elements and expression quantitative trait loci (eQTLs) in nine tissues. We used T2D tissue-grouped variant sets as genetic instruments to conduct 2-Sample Mendelian Randomization (MR) in ten related outcomes whose risk is increased by T2D using the FinnGen cohort. We performed PheWAS analysis to investigate whether the T2D tissue-grouped variant sets had specific predicted disease signatures. We identified an average of 176 variants acting in nine tissues implicated in T2D, and an average of 30 variants acting on regulatory elements that are unique to the nine tissues of interest. In 2-Sample MR analyses, all subsets of regulatory variants acting in different tissues were associated with increased risk of the ten secondary outcomes studied on similar levels. No tissue-grouped variant set was associated with an outcome significantly more than other tissue-grouped variant sets. We did not identify different disease progression profiles based on tissue-specific regulatory and transcriptome information. Bigger sample sizes and other layers of regulatory information in critical tissues may help identify subsets of T2D variants that are implicated in certain secondary outcomes, uncovering system-specific disease progression.

Fine mapping spatiotemporal mechanisms of genetic variants underlying cardiac traits and disease

The causal variants and genes underlying thousands of cardiac GWAS signals have yet to be identified. Here, we leverage spatiotemporal information on 966 RNA-seq cardiac samples and perform an expression quantitative trait locus (eQTL) analysis detecting eQTLs considering both eGenes and eIsoforms. We identify 2,578 eQTLs associated with a specific developmental stage-, tissue- and/or cell type. Colocalization between eQTL and GWAS signals of five cardiac traits identified variants with high posterior probabilities for being causal in 210 GWAS loci. Pulse pressure GWAS loci are enriched for colocalization with fetal- and smooth muscle- eQTLs; pulse rate with adult- and cardiac muscle- eQTLs; and atrial fibrillation with cardiac muscle- eQTLs. Fine mapping identifies 79 credible sets with five or fewer SNPs, of which 15 were associated with spatiotemporal eQTLs. Our study shows that many cardiac GWAS variants impact traits and disease in a developmental stage-, tissue- and/or cell type-specific fashion.

Type 2 diabetes candidate genes, including PAX5, cause impaired insulin secretion in human pancreatic islets

Type 2 diabetes (T2D) is caused by insufficient insulin secretion from pancreatic β cells. To identify candidate genes contributing to T2D pathophysiology, we studied human pancreatic islets from approximately 300 individuals. We found 395 differentially expressed genes (DEGs) in islets from individuals with T2D, including, to our knowledge, novel (OPRD1, PAX5, TET1) and previously identified (CHL1, GLRA1, IAPP) candidates. A third of the identified expression changes in islets may predispose to diabetes, as expression of these genes associated with HbA1c in individuals not previously diagnosed with T2D. Most DEGs were expressed in human β cells, based on single-cell RNA-Seq data. Additionally, DEGs displayed alterations in open chromatin and associated with T2D SNPs. Mouse KO strains demonstrated that the identified T2D-associated candidate genes regulate glucose homeostasis and body composition in vivo. Functional validation showed that mimicking T2D-associated changes for OPRD1, PAX5, and SLC2A2 impaired insulin secretion. Impairments in Pax5-overexpressing β cells were due to severe mitochondrial dysfunction. Finally, we discovered PAX5 as a potential transcriptional regulator of many T2D-associated DEGs in human islets. Overall, we have identified molecular alterations in human pancreatic islets that contribute to β cell dysfunction in T2D pathophysiology.

Human pancreatic islet microRNAs implicated in diabetes and related traits by large-scale genetic analysis

Genetic studies have identified ≥240 loci associated with the risk of type 2 diabetes (T2D), yet most of these loci lie in non-coding regions, masking the underlying molecular mechanisms. Recent studies investigating mRNA expression in human pancreatic islets have yielded important insights into the molecular drivers of normal islet function and T2D pathophysiology. However, similar studies investigating microRNA (miRNA) expression remain limited. Here, we present data from 63 individuals, the largest sequencing-based analysis of miRNA expression in human islets to date. We characterized the genetic regulation of miRNA expression by decomposing the expression of highly heritable miRNAs into cis- and trans-acting genetic components and mapping cis-acting loci associated with miRNA expression [miRNA-expression quantitative trait loci (eQTLs)]. We found i) 84 heritable miRNAs, primarily regulated by trans-acting genetic effects, and ii) 5 miRNA-eQTLs. We also used several different strategies to identify T2D-associated miRNAs. First, we colocalized miRNA-eQTLs with genetic loci associated with T2D and multiple glycemic traits, identifying one miRNA, miR-1908, that shares genetic signals for blood glucose and glycated hemoglobin (HbA1c). Next, we intersected miRNA seed regions and predicted target sites with credible set SNPs associated with T2D and glycemic traits and found 32 miRNAs that may have altered binding and function due to disrupted seed regions. Finally, we performed differential expression analysis and identified 14 miRNAs associated with T2D status-including miR-187-3p, miR-21-5p, miR-668, and miR-199b-5p-and 4 miRNAs associated with a polygenic score for HbA1c levels-miR-216a, miR-25, miR-30a-3p, and miR-30a-5p.

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.

Estimation and implications of the genetic architecture of fasting and non-fasting blood glucose

The genetic regulation of post-prandial glucose levels is poorly understood. Here, we characterise the genetic architecture of blood glucose variably measured within 0 and 24 h of fasting in 368,000 European ancestry participants of the UK Biobank. We found a near-linear increase in the heritability of non-fasting glucose levels over time, which plateaus to its fasting state value after 5 h post meal (h² = 11%; standard error: 1%). The genetic correlation between different fasting times is > 0.77, suggesting that the genetic control of glucose is largely constant across fasting durations. Accounting for heritability differences between fasting times leads to a ~16% improvement in the discovery of genetic variants associated with glucose. Newly detected variants improve the prediction of fasting glucose and type 2 diabetes in independent samples. Finally, we meta-analysed summary statistics from genome-wide association studies of random and fasting glucose (N = 518,615) and identified 156 independent SNPs explaining 3% of fasting glucose variance. Altogether, our study demonstrates the utility of random glucose measures to improve the discovery of genetic variants associated with glucose homeostasis, even in fasting conditions.

Insights from multi-omics integration in complex disease primary tissues

Genome-wide association studies (GWAS) have provided insights into the genetic basis of complex diseases. In the next step, integrative multi-omics approaches can characterize molecular profiles in relevant primary tissues to reveal the mechanisms that underlie disease development. Here, we highlight recent progress in four examples of complex diseases generated by integrative studies: type 2 diabetes (T2D), osteoarthritis, Alzheimer's disease (AD), and systemic lupus erythematosus (SLE). High-resolution methodologies such as single-cell and spatial omics techniques will become even more important in the future. Furthermore, we emphasize the urgent need to include as yet understudied cell types and increase the diversity of studied populations.

A genome-wide CRISPR screen identifies CALCOCO2 as a regulator of beta cell function influencing type 2 diabetes risk

Identification of the genes and processes mediating genetic association signals for complex diseases represents a major challenge. As many of the genetic signals for type 2 diabetes (T2D) exert their effects through pancreatic islet-cell dysfunction, we performed a genome-wide pooled CRISPR loss-of-function screen in a human pancreatic beta cell line. We assessed the regulation of insulin content as a disease-relevant readout of beta cell function and identified 580 genes influencing this phenotype. Integration with genetic and genomic data provided experimental support for 20 candidate T2D effector transcripts including the autophagy receptor CALCOCO2. Loss of CALCOCO2 was associated with distorted mitochondria, less proinsulin-containing immature granules and accumulation of autophagosomes upon inhibition of late-stage autophagy. Carriers of T2D-associated variants at the CALCOCO2 locus further displayed altered insulin secretion. Our study highlights how cellular screens can augment existing multi-omic efforts to support mechanistic understanding and provide evidence for causal effects at genome-wide association studies loci.

The landscape of expression and alternative splicing variation across human traits

Understanding the consequences of individual transcriptome variation is fundamental to deciphering human biology and disease. We implement a statistical framework to quantify the contributions of 21 individual traits as drivers of gene expression and alternative splicing variation across 46 human tissues and 781 individuals from the Genotype-Tissue Expression project. We demonstrate that ancestry, sex, age, and BMI make additive and tissue-specific contributions to expression variability, whereas interactions are rare. Variation in splicing is dominated by ancestry and is under genetic control in most tissues, with ribosomal proteins showing a strong enrichment of tissue-shared splicing events. Our analyses reveal a systemic contribution of types 1 and 2 diabetes to tissue transcriptome variation with the strongest signal in the nerve, where histopathology image analysis identifies novel genes related to diabetic neuropathy. Our multi-tissue and multi-trait approach provides an extensive characterization of the main drivers of human transcriptome variation in health and disease.

Type 1 diabetes risk genes mediate pancreatic beta cell survival in response to proinflammatory cytokines

We combined functional genomics and human genetics to investigate processes that affect type 1 diabetes (T1D) risk by mediating beta cell survival in response to proinflammatory cytokines. We mapped 38,931 cytokine-responsive candidate cis-regulatory elements (cCREs) in beta cells using ATAC-seq and snATAC-seq and linked them to target genes using co-accessibility and HiChIP. Using a genome-wide CRISPR screen in EndoC-βH1 cells, we identified 867 genes affecting cytokine-induced survival, and genes promoting survival and up-regulated in cytokines were enriched at T1D risk loci. Using SNP-SELEX, we identified 2,229 variants in cytokine-responsive cCREs altering transcription factor (TF) binding, and variants altering binding of TFs regulating stress, inflammation, and apoptosis were enriched for T1D risk. At the 16p13 locus, a fine-mapped T1D variant altering TF binding in a cytokine-induced cCRE interacted with SOCS1, which promoted survival in cytokine exposure. Our findings reveal processes and genes acting in beta cells during inflammation that modulate T1D risk.

The human batokine EPDR1 regulates β-cell metabolism and function

OBJECTIVE

Ependymin-Related Protein 1 (EPDR1) was recently identified as a secreted human batokine regulating mitochondrial respiration linked to thermogenesis in brown fat. Despite that EPDR1 is expressed in human pancreatic β-cells and that glucose-stimulated mitochondrial metabolism is critical for stimulus-secretion coupling in β-cells, the role of EPDR1 in β-cell metabolism and function has not been investigated.

METHODS

EPDR1 mRNA levels in human pancreatic islets from non-diabetic (ND) and type 2 diabetes (T2D) subjects were assessed. Human islets, EndoC-βH1 and INS1 832/13 cells were transfected with scramble (control) and EPDR1 siRNAs (EPDR1-KD) or treated with human EPDR1 protein, and glucose-stimulated insulin secretion (GSIS) assessed by ELISA. Mitochondrial metabolism was investigated by extracellular flux analyzer, confocal microscopy and mass spectrometry-based metabolomics analysis.

RESULTS

EPDR1 mRNA expression was upregulated in human islets from T2D and obese donors and positively correlated to BMI of donors. In T2D donors, EPDR1 mRNA levels negatively correlated with HbA1c and positively correlated with GSIS. EPDR1 silencing in human islets and β-cell lines reduced GSIS whereas treatment with human EPDR1 protein increased GSIS. Epdr1 silencing in INS1 832/13 cells reduced glucose- and pyruvate- but not K⁺-stimulated insulin secretion. Metabolomics analysis in Epdr1-KD INS1 832/13 cells suggests diversion of glucose-derived pyruvate to lactate production and decreased malate-aspartate shuttle and the tricarboxylic acid (TCA) cycle activity. The glucose-stimulated rise in mitochondrial respiration and ATP/ADP-ratio was impaired in Epdr1-deficient cells.

CONCLUSION

These results suggests that to maintain glucose homeostasis in obese people, upregulation of EPDR1 may improve β-cell function via channelling glycolysis-derived pyruvate to the mitochondrial TCA cycle.

Zmiz1 is required for mature β-cell function and mass expansion upon high fat feeding

OBJECTIVE

Identifying the transcripts which mediate genetic association signals for type 2 diabetes (T2D) is critical to understand disease mechanisms. Studies in pancreatic islets support the transcription factor ZMIZ1 as a transcript underlying a T2D GWAS signal, but how it influences T2D risk is unknown.

METHODS

β-Cell-specific Zmiz1 knockout (Zmiz1^βKO) mice were generated and phenotypically characterised. Glucose homeostasis was assessed in Zmiz1^βKO mice and their control littermates on chow diet (CD) and high fat diet (HFD). Islet morphology and function were examined by immunohistochemistry and in vitro islet function was assessed by dynamic insulin secretion assay. Transcript and protein expression were assessed by RNA sequencing and Western blotting. In islets isolated from genotyped human donors, we assessed glucose-dependent insulin secretion and islet insulin content by static incubation assay.

RESULTS

Male and female Zmiz1^βKO mice were glucose intolerant with impaired insulin secretion, compared with control littermates. Transcriptomic profiling of Zmiz1^βKO islets identified over 500 differentially expressed genes including those involved in β-cell function and maturity, which we confirmed at the protein level. Upon HFD, Zmiz1^βKO mice fail to expand β-cell mass and become severely diabetic. Human islets from carriers of the ZMIZ1-linked T2D-risk alleles have reduced islet insulin content and glucose-stimulated insulin secretion.

CONCLUSIONS

β-Cell Zmiz1 is required for normal glucose homeostasis. Genetic variation at the ZMIZ1 locus may influence T2D-risk by reducing islet mass expansion upon metabolic stress and the ability to maintain a mature β-cell state.

HiChIPdb: a comprehensive database of HiChIP regulatory interactions

Elucidating the role of 3D architecture of DNA in gene regulation is crucial for understanding cell differentiation, tissue homeostasis and disease development. Among various chromatin conformation capture methods, HiChIP has received increasing attention for its significant improvement over other methods in profiling of regulatory (e.g. H3K27ac) and structural (e.g. cohesin) interactions. To facilitate the studies of 3D regulatory interactions, we developed a HiChIP interactions database, HiChIPdb (http://health.tsinghua.edu.cn/hichipdb/). The current version of HiChIPdb contains ∼262M annotated HiChIP interactions from 200 high-throughput HiChIP samples across 108 cell types. The functionalities of HiChIPdb include: (i) standardized categorization of HiChIP interactions in a hierarchical structure based on organ, tissue and cell line and (ii) comprehensive annotations of HiChIP interactions with regulatory genes and GWAS Catalog SNPs. To the best of our knowledge, HiChIPdb is the first comprehensive database that utilizes a unified pipeline to map the functional interactions across diverse cell types and tissues in different resolutions. We believe this database has the potential to advance cutting-edge research in regulatory mechanisms in development and disease by removing the barrier in data aggregation, preprocessing, and analysis.

Genetic regulation of RNA splicing in human pancreatic islets

BACKGROUND

Non-coding genetic variants that influence gene transcription in pancreatic islets play a major role in the susceptibility to type 2 diabetes (T2D), and likely also contribute to type 1 diabetes (T1D) risk. For many loci, however, the mechanisms through which non-coding variants influence diabetes susceptibility are unknown.

RESULTS

We examine splicing QTLs (sQTLs) in pancreatic islets from 399 human donors and observe that common genetic variation has a widespread influence on the splicing of genes with established roles in islet biology and diabetes. In parallel, we profile expression QTLs (eQTLs) and use transcriptome-wide association as well as genetic co-localization studies to assign islet sQTLs or eQTLs to T2D and T1D susceptibility signals, many of which lack candidate effector genes. This analysis reveals biologically plausible mechanisms, including the association of T2D with an sQTL that creates a nonsense isoform in ERO1B, a regulator of ER-stress and proinsulin biosynthesis. The expanded list of T2D risk effector genes reveals overrepresented pathways, including regulators of G-protein-mediated cAMP production. The analysis of sQTLs also reveals candidate effector genes for T1D susceptibility such as DCLRE1B, a senescence regulator, and lncRNA MEG3.

CONCLUSIONS

These data expose widespread effects of common genetic variants on RNA splicing in pancreatic islets. The results support a role for splicing variation in diabetes susceptibility, and offer a new set of genetic targets with potential therapeutic benefit.

Islet Gene View-a tool to facilitate islet research

Characterization of gene expression in pancreatic islets and its alteration in type 2 diabetes (T2D) are vital in understanding islet function and T2D pathogenesis. We leveraged RNA sequencing and genome-wide genotyping in islets from 188 donors to create the Islet Gene View (IGW) platform to make this information easily accessible to the scientific community. Expression data were related to islet phenotypes, diabetes status, other islet-expressed genes, islet hormone-encoding genes and for expression in insulin target tissues. The IGW web application produces output graphs for a particular gene of interest. In IGW, 284 differentially expressed genes (DEGs) were identified in T2D donor islets compared with controls. Forty percent of DEGs showed cell-type enrichment and a large proportion significantly co-expressed with islet hormone-encoding genes; glucagon (<i>GCG</i>, 56%), amylin (<i>IAPP</i>, 52%), insulin (<i>INS</i>, 44%), and somatostatin (<i>SST</i>, 24%). Inhibition of two DEGs, <i>UNC5D</i> and <i>SERPINE2</i>, impaired glucose-stimulated insulin secretion and impacted cell survival in a human β-cell model. The exploratory use of IGW could help designing more comprehensive functional follow-up studies and serve to identify therapeutic targets in T2D.

TCF7L2 plays a complex role in human adipose progenitor biology, which might contribute to genetic susceptibility to type 2 diabetes

INTRODUCTION

Non-coding genetic variation at TCF7L2 is the strongest genetic determinant of type 2 diabetes (T2D) risk in humans. TCF7L2 encodes a transcription factor mediating the nuclear effects of WNT signaling in adipose tissue (AT). In vivo studies in transgenic mice have highlighted important roles for TCF7L2 in adipose tissue biology and systemic metabolism.

OBJECTIVE

To map the expression of TCF7L2 in human AT, examine its role in human adipose cell biology in vitro, and investigate the effects of the fine-mapped T2D-risk allele at rs7903146 on AT morphology and TCF7L2 expression.

METHODS

Ex vivo gene expression studies of TCF7L2 in whole and fractionated human AT. In vitro TCF7L2 gain- and/or loss-of-function studies in primary and immortalized human adipose progenitor cells (APCs) and mature adipocytes (mADs). AT phenotyping of rs7903146 T2D-risk variant carriers and matched controls.

RESULTS

Adipose progenitors (APs) exhibited the highest TCF7L2 mRNA abundance compared to mature adipocytes and adipose-derived endothelial cells. Obesity was associated with reduced TCF7L2 transcript levels in whole subcutaneous abdominal AT but paradoxically increased expression in APs. In functional studies, TCF7L2 knockdown (KD) in abdominal APs led to dose-dependent activation of WNT/β-catenin signaling, impaired proliferation and dose-dependent effects on adipogenesis. Whilst partial KD enhanced adipocyte differentiation, near-total KD impaired lipid accumulation and adipogenic gene expression. Over-expression of TCF7L2 accelerated adipogenesis. In contrast, TCF7L2-KD in gluteal APs dose-dependently enhanced lipid accumulation. Transcriptome-wide profiling revealed that TCF7L2 might modulate multiple aspects of AP biology including extracellular matrix secretion, immune signaling and apoptosis. The T2D-risk allele at rs7903146 was associated with reduced AP TCF7L2 expression and enhanced AT insulin sensitivity.

CONCLUSIONS

TCF7L2 plays a complex role in AP biology and has both dose- and depot-dependent effects on adipogenesis. In addition to regulating pancreatic insulin secretion, genetic variation at TCF7L2 might also influence T2D risk by modulating AP function.

Every islet matters: improving the impact of human islet research

Detailed characterization of human pancreatic islets is key to elucidating the pathophysiology of all forms of diabetes, especially type 2 diabetes. However, access to human pancreatic islets is limited. Pancreatic tissue for islet retrieval can be obtained from brain-dead organ donors or from individuals undergoing pancreatectomy, often referred to as 'living donors'. Different protocols for human islet procurement can substantially impact islet function. This variability, coupled with heterogeneity between individuals and islets, results in analytical challenges to separate genuine disease pathology or differences between human donors from experimental noise. There are currently no international guidelines for human donor phenotyping, islet procurement and functional characterization. This lack of standardization means that substantial investments from multiple international efforts towards improved understanding of diabetes pathology cannot be fully leveraged. In this Perspective, we overview the status of the field of human islet research, highlight the challenges and propose actions that could accelerate research progress and increase understanding of type 2 diabetes to slow its pandemic spreading.

Whole genome sequence association analysis of fasting glucose and fasting insulin levels in diverse cohorts from the NHLBI TOPMed program

The genetic determinants of fasting glucose (FG) and fasting insulin (FI) have been studied mostly through genome arrays, resulting in over 100 associated variants. We extended this work with high-coverage whole genome sequencing analyses from fifteen cohorts in NHLBI's Trans-Omics for Precision Medicine (TOPMed) program. Over 23,000 non-diabetic individuals from five race-ethnicities/populations (African, Asian, European, Hispanic and Samoan) were included. Eight variants were significantly associated with FG or FI across previously identified regions MTNR1B, G6PC2, GCK, GCKR and FOXA2. We additionally characterize suggestive associations with FG or FI near previously identified SLC30A8, TCF7L2, and ADCY5 regions as well as APOB, PTPRT, and ROBO1. Functional annotation resources including the Diabetes Epigenome Atlas were compiled for each signal (chromatin states, annotation principal components, and others) to elucidate variant-to-function hypotheses. We provide a catalog of nucleotide-resolution genomic variation spanning intergenic and intronic regions creating a foundation for future sequencing-based investigations of glycemic traits.

Trends in insulin resistance: insights into mechanisms and therapeutic strategy

The centenary of insulin discovery represents an important opportunity to transform diabetes from a fatal diagnosis into a medically manageable chronic condition. Insulin is a key peptide hormone and mediates the systemic glucose metabolism in different tissues. Insulin resistance (IR) is a disordered biological response for insulin stimulation through the disruption of different molecular pathways in target tissues. Acquired conditions and genetic factors have been implicated in IR. Recent genetic and biochemical studies suggest that the dysregulated metabolic mediators released by adipose tissue including adipokines, cytokines, chemokines, excess lipids and toxic lipid metabolites promote IR in other tissues. IR is associated with several groups of abnormal syndromes that include obesity, diabetes, metabolic dysfunction-associated fatty liver disease (MAFLD), cardiovascular disease, polycystic ovary syndrome (PCOS), and other abnormalities. Although no medication is specifically approved to treat IR, we summarized the lifestyle changes and pharmacological medications that have been used as efficient intervention to improve insulin sensitivity. Ultimately, the systematic discussion of complex mechanism will help to identify potential new targets and treat the closely associated metabolic syndrome of IR.

Powerful and robust inference of complex phenotypes' causal genes with dependent expression quantitative loci by a median-based Mendelian randomization

Isolating the causal genes from numerous genetic association signals in genome-wide association studies (GWASs) of complex phenotypes remains an open and challenging question. In the present study, we proposed a statistical approach, the effective-median-based Mendelian randomization (MR) framework, for inferring the causal genes of complex phenotypes with the GWAS summary statistics (named EMIC). The effective-median method solved the high false-positive issue in the existing MR methods due to either correlation among instrumental variables or noises in approximated linkage disequilibrium (LD). EMIC can further perform a pleiotropy fine-mapping analysis to remove possible false-positive estimates. With the usage of multiple cis-expression quantitative trait loci (eQTLs), EMIC was also more powerful than the alternative methods for the causal gene inference in the simulated datasets. Furthermore, EMIC rediscovered many known causal genes of complex phenotypes (schizophrenia, bipolar disorder, and total cholesterol) and reported many new and promising candidate causal genes. In sum, this study provided an efficient solution to discriminate the candidate causal genes from vast amounts of GWAS signals with eQTLs. EMIC has been implemented in our integrative software platform KGGSEE.