A common functional regulatory variant at a type 2 diabetes locus upregulates ARAP1 expression in the pancreatic beta cell

Insulin-related traits and prostate cancer: A Mendelian randomization study

Investigating the causal relationship between insulin secretion and prostate cancer (PCa) development is challenging due to the multifactorial nature of PCa, which complicates the isolation of the specific impact of insulin-related factors. We conducted a Mendelian randomization (MR) study to investigate the associations between insulin secretion-related traits and PCa. We used 36, 60, 56, 23, 48, and 49 single nucleotide polymorphisms (SNPs) as instrumental variables for fasting insulin, insulin sensitivity, proinsulin, and proinsulin in nondiabetic individuals, individuals with diabetes, and individuals receiving exogenous insulin, respectively. These SNPs were selected from various genome-wide association studies. To clarify the causal relationship between insulin-related traits and PCa, we utilized a multivariable MR analysis to adjust for obesity and body fat percentage. Additionally, two-step Mendelian randomization was conducted to assess the role of insulin-like growth factor 1 (IGF-1) in the relationship between proinsulin and PCa. Two-sample MR analysis revealed strong associations between genetically predicted fasting insulin, insulin sensitivity, proinsulin, and proinsulin in nondiabetic individuals and the development of PCa. After adjustment for obesity and body fat percentage using multivariable MR analysis, proinsulin remained significantly associated with PCa, whereas other factors were not. Furthermore, two-step MR analysis demonstrated that proinsulin acts as a negative factor in prostate carcinogenesis, largely independent of IGF-1. This study provides evidence suggesting that proinsulin may act as a negative factor contributing to the development of PCa. Novel therapies targeting proinsulin may have potential benefits for PCa patients, potentially reducing the need for unnecessary surgical treatments.

A comprehensive in silico investigation into the pathogenic SNPs in the RTEL1 gene and their biological consequences

The Regulator of Telomere Helicase 1 (RTEL1) gene encodes a critical DNA helicase intricately involved in the maintenance of telomeric structures and the preservation of genomic stability. Germline mutations in the RTEL1 gene have been clinically associated with Hoyeraal-Hreidarsson syndrome, a more severe version of Dyskeratosis Congenita. Although various research has sought to link RTEL1 mutations to specific disorders, no comprehensive investigation has yet been conducted on missense mutations. In this study, we attempted to investigate the functionally and structurally deleterious coding and non-coding SNPs of the RTEL1 gene using an in silico approach. Initially, out of 1392 nsSNPs, 43 nsSNPs were filtered out through ten web-based bioinformatics tools. With subsequent analysis using nine in silico tools, these 43 nsSNPs were further shortened to 11 most deleterious nsSNPs. Furthermore, analyses of mutated protein structures, evolutionary conservancy, surface accessibility, domains & PTM sites, cancer susceptibility, and interatomic interaction revealed the detrimental effect of these 11 nsSNPs on RTEL1 protein. An in-depth investigation through molecular docking with the DNA binding sequence demonstrated a striking change in the interaction pattern for F15L, M25V, and G706R mutant proteins, suggesting the more severe consequences of these mutations on protein structure and functionality. Among the non-coding variants, two had the highest likelihood of being regulatory variants, whereas one variant was predicted to affect the target region of a miRNA. Thus, this study lays the groundwork for extensive analysis of RTEL1 gene variants in the future, along with the advancement of precision medicine and other treatment modalities.

Characterizing Common and Rare Variations in Nontraditional Glycemic Biomarkers Using Multivariate Approaches on Multiancestry ARIC Study

Genetic studies of nontraditional glycemic biomarkers, glycated albumin and fructosamine, can shed light on unknown aspects of type 2 diabetes genetics and biology. We performed a multiphenotype genome-wide association study of glycated albumin and fructosamine from 7,395 White and 2,016 Black participants in the Atherosclerosis Risk in Communities (ARIC) study on common variants from genotyped/imputed data. We discovered two genome-wide significant loci, one mapping to a known type 2 diabetes gene (ARAP1/STARD10) and another mapping to a novel region (UGT1A complex of genes), using multiomics gene-mapping strategies in diabetes-relevant tissues. We identified additional loci that were ancestry- and sex-specific (e.g., PRKCA in African ancestry, FCGRT in European ancestry, TEX29 in males). Further, we implemented multiphenotype gene-burden tests on whole-exome sequence data from 6,590 White and 2,309 Black ARIC participants. Ten variant sets annotated to genes across different variant aggregation strategies were exome-wide significant only in multiancestry analysis, of which CD1D, EGFL7/AGPAT2, and MIR126 had notable enrichment of rare predicted loss of function variants in African ancestry despite smaller sample sizes. Overall, 8 of 14 discovered loci and genes were implicated to influence these biomarkers via glycemic pathways, and most of them were not previously implicated in studies of type 2 diabetes. This study illustrates improved locus discovery and potential effector gene discovery by leveraging joint patterns of related biomarkers across the entire allele frequency spectrum in multiancestry analysis. Future investigation of the loci and genes potentially acting through glycemic pathways may help us better understand the risk of developing type 2 diabetes.

ARTICLE HIGHLIGHTS

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.

Epromoters are new players in the regulatory landscape with potential pleiotropic roles

Precise spatiotemporal control of gene expression during normal development and cell differentiation is achieved by the combined action of proximal (promoters) and distal (enhancers) cis-regulatory elements. Recent studies have reported that a subset of promoters, termed Epromoters, works also as enhancers to regulate distal genes. This new paradigm opened novel questions regarding the complexity of our genome and raises the possibility that genetic variation within Epromoters has pleiotropic effects on various physiological and pathological traits by differentially impacting multiple proximal and distal genes. Here, we discuss the different observations pointing to an important role of Epromoters in the regulatory landscape and summarize the evidence supporting a pleiotropic impact of these elements in disease. We further hypothesize that Epromoter might represent a major contributor to phenotypic variation and disease.

Pleckstrin Homology [PH] domain, structure, mechanism, and contribution to human disease

The pleckstrin homology [PH] domain is a structural fold found in more than 250 proteins making it the 11th most common domain in the human proteome. 25% of family members have more than one PH domain and some PH domains are split by one, or several other, protein domains although still folding to give functioning PH domains. We review mechanisms of PH domain activity, the role PH domain mutation plays in human disease including cancer, hyperproliferation, neurodegeneration, inflammation, and infection, and discuss pharmacotherapeutic approaches to regulate PH domain activity for the treatment of human disease. Almost half PH domain family members bind phosphatidylinositols [PIs] that attach the host protein to cell membranes where they interact with other membrane proteins to give signaling complexes or cytoskeleton scaffold platforms. A PH domain in its native state may fold over other protein domains thereby preventing substrate access to a catalytic site or binding with other proteins. The resulting autoinhibition can be released by PI binding to the PH domain, or by protein phosphorylation thus providing fine tuning of the cellular control of PH domain protein activity. For many years the PH domain was thought to be undruggable until high-resolution structures of human PH domains allowed structure-based design of novel inhibitors that selectively bind the PH domain. Allosteric inhibitors of the Akt1 PH domain have already been tested in cancer patients and for proteus syndrome, with several other PH domain inhibitors in preclinical development for treatment of other human diseases.

Characterizing common and rare variations in non-traditional glycemic biomarkers using multivariate approaches on multi-ancestry ARIC study

Glycated hemoglobin, fasting glucose, glycated albumin, and fructosamine are biomarkers that reflect different aspects of the glycemic process. Genetic studies of these glycemic biomarkers can shed light on unknown aspects of type 2 diabetes genetics and biology. While there exists several GWAS of glycated hemoglobin and fasting glucose, very few GWAS have focused on glycated albumin or fructosamine. We performed a multi-phenotype GWAS of glycated albumin and fructosamine from 7,395 White and 2,016 Black participants in the Atherosclerosis Risk in Communities (ARIC) study on the common variants from genotyped/imputed data. We found 2 genome-wide significant loci, one mapping to known type 2 diabetes gene (ARAP1/STARD10, p = 2.8 × 10-8) and another mapping to a novel gene (UGT1A, p = 1.4 × 10-8) using multi-omics gene mapping strategies in diabetes-relevant tissues. We identified additional loci that were ancestry-specific (e.g., PRKCA from African ancestry individuals, p = 1.7 × 10-8) and sex-specific (TEX29 locus in males only, p = 3.0 × 10-8). Further, we implemented multi-phenotype gene-burden tests on whole-exome sequence data from 6,590 White and 2,309 Black ARIC participants. Eleven genes across different rare variant aggregation strategies were exome-wide significant only in multi-ancestry analysis. Four out of 11 genes had notable enrichment of rare predicted loss of function variants in African ancestry participants despite smaller sample size. Overall, 8 out of 15 loci/genes were implicated to influence these biomarkers via glycemic pathways. This study illustrates improved locus discovery and potential effector gene discovery by leveraging joint patterns of related biomarkers across entire allele frequency spectrum in multi-ancestry analyses. Most of the loci/genes we identified have not been previously implicated in studies of type 2 diabetes, and future investigation of the loci/genes potentially acting through glycemic pathways may help us better understand risk of developing type 2 diabetes.

Bioinformatics and Next-Generation Data Analysis for Identification of Genes and Molecular Pathways Involved in Subjects with Diabetes and Obesity

Background and Objectives: A subject with diabetes and obesity is a class of the metabolic disorder. The current investigation aimed to elucidate the potential biomarker and prognostic targets in subjects with diabetes and obesity. Materials and Methods: The next-generation sequencing (NGS) data of GSE132831 was downloaded from Gene Expression Omnibus (GEO) database. Functional enrichment analysis of DEGs was conducted with ToppGene. The protein-protein interactions network, module analysis, target gene-miRNA regulatory network and target gene-TF regulatory network were constructed and analyzed. Furthermore, hub genes were validated by receiver operating characteristic (ROC) analysis. A total of 872 DEGs, including 439 up-regulated genes and 433 down-regulated genes were observed. Results: Second, functional enrichment analysis showed that these DEGs are mainly involved in the axon guidance, neutrophil degranulation, plasma membrane bounded cell projection organization and cell activation. The top ten hub genes (MYH9, FLNA, DCTN1, CLTC, ERBB2, TCF4, VIM, LRRK2, IFI16 and CAV1) could be utilized as potential diagnostic indicators for subjects with diabetes and obesity. The hub genes were validated in subjects with diabetes and obesity. Conclusion: This investigation found effective and reliable molecular biomarkers for diagnosis and prognosis by integrated bioinformatics analysis, suggesting new and key therapeutic targets for subjects with diabetes and obesity.

Loci Associated with Negative Heterosis for Viability and Meat Productivity in Interspecific Sheep Hybrids

Negative heterosis can occur on different economically important traits, but the exact biological mechanisms of this phenomenon are still unknown. The present study focuses on determining the genetic factors associated with negative heterosis in interspecific hybrids between domestic sheep (Ovis aries) and argali (Ovis ammon). One locus (rs417431015) associated with viability and two loci (rs413302370, rs402808951) associated with meat productivity were identified. One gene (ARAP2) was prioritized for viability and three for meat productivity (PDE2A, ARAP1, and PCDH15). The loci associated with meat productivity were demonstrated to fit the overdominant inheritance model and could potentially be involved int negative heterosis mechanisms.

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.

Modeling common and rare genetic risk factors of neuropsychiatric disorders in human induced pluripotent stem cells

Recent genome-wide association studies (GWAS) and whole-exome sequencing of neuropsychiatric disorders, especially schizophrenia, have identified a plethora of common and rare disease risk variants/genes. Translating the mounting human genetic discoveries into novel disease biology and more tailored clinical treatments is tied to our ability to causally connect genetic risk variants to molecular and cellular phenotypes. When combined with the Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)/CRISPR-associated (Cas) nuclease-mediated genome editing system, human induced pluripotent stem cell (hiPSC)-derived neural cultures (both 2D and 3D organoids) provide a promising tractable cellular model for bridging the gap between genetic findings and disease biology. In this review, we first conceptualize the advances in understanding the disease polygenicity and convergence from the past decade of iPSC modeling of different types of genetic risk factors of neuropsychiatric disorders. We then discuss the major cell types and cellular phenotypes that are most relevant to neuropsychiatric disorders in iPSC modeling. Finally, we critically review the limitations of iPSC modeling of neuropsychiatric disorders and outline the need for implementing and developing novel methods to scale up the number of iPSC lines and disease risk variants in a systematic manner. Sufficiently scaled-up iPSC modeling and a better functional interpretation of genetic risk variants, in combination with cutting-edge CRISPR/Cas9 gene editing and single-cell multi-omics methods, will enable the field to identify the specific and convergent molecular and cellular phenotypes in precision for neuropsychiatric disorders.

The landscape of GWAS validation; systematic review identifying 309 validated non-coding variants across 130 human diseases

Background

The remarkable growth of genome-wide association studies (GWAS) has created a critical need to experimentally validate the disease-associated variants, 90% of which involve non-coding variants.

Methods

To determine how the field is addressing this urgent need, we performed a comprehensive literature review identifying 36,676 articles. These were reduced to 1454 articles through a set of filters using natural language processing and ontology-based text-mining. This was followed by manual curation and cross-referencing against the GWAS catalog, yielding a final set of 286 articles.

Results

We identified 309 experimentally validated non-coding GWAS variants, regulating 252 genes across 130 human disease traits. These variants covered a variety of regulatory mechanisms. Interestingly, 70% (215/309) acted through cis-regulatory elements, with the remaining through promoters (22%, 70/309) or non-coding RNAs (8%, 24/309). Several validation approaches were utilized in these studies, including gene expression (n = 272), transcription factor binding (n = 175), reporter assays (n = 171), in vivo models (n = 104), genome editing (n = 96) and chromatin interaction (n = 33).

Conclusions

This review of the literature is the first to systematically evaluate the status and the landscape of experimentation being used to validate non-coding GWAS-identified variants. Our results clearly underscore the multifaceted approach needed for experimental validation, have practical implications on variant prioritization and considerations of target gene nomination. While the field has a long way to go to validate the thousands of GWAS associations, we show that progress is being made and provide exemplars of validation studies covering a wide variety of mechanisms, target genes, and disease areas.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01216-w.

Shared Genetic Basis and Causal Relationship Between Television Watching, Breakfast Skipping and Type 2 Diabetes: Evidence From a Comprehensive Genetic Analysis

Background

Epidemiological investigations have established unhealthy lifestyles, such as excessive leisurely sedentary behavior (especially TV/television watching) and breakfast skipping, increase the risk of type 2 diabetes (T2D), but the causal relationship is unclear. We aimed to understand how single nucleotide variants contribute to the co-occurrence of unhealthy lifestyles and T2D, thereby providing meaningful insights into disease mechanisms.

Methods

Combining summary statistics from genome-wide association studies (GWAS) on TV watching (N = 422218), breakfast skipping (N = 193860) and T2D (N = 159208) in European pedigrees, we conducted comprehensive pairwise genetic analysis, including high-definition likelihood (HDL-method), cross-phenotype association studies (CPASSOC), GWAS-eQTL colocalization analysis and transcriptome-wide association studies (TWAS), to understand the genetic overlap between them. We also performed bidirectional two-sample Mendelian randomization (MR) analysis for causal inference using genetic instrumental variables, and two-step MR mediation analysis was used to assess any effects explained by body mass index, lipid traits and glycemic traits.

Results

HDL-method showed that T2D shared a strong genetic correlation with TV watching (rg = 0.26; P = 1.63×10-29) and skipping breakfast (rg = 0.15; P =2.02×10-6). CPASSOC identifies eight independent SNPs shared between T2D and TV watching, including one novel shared locus. TWAS and CPASSOC showed that shared genes were enriched in lung, esophageal, adipose, and thyroid tissues and highlighted potential shared regulatory pathways for lipoprotein metabolism, pancreatic β-cell function, cellular senescence and multi-mediator factors. MR showed TV watching had a causal effect on T2D (βIVW = 0.629, PIVW = 1.80×10-10), but no significant results were observed between breakfast skipping and T2D. Mediation analysis provided evidence that body mass index, fasting glucose, hemoglobin A1c and high-density lipoprotein are potential factors that mediate the causal relationship between TV and T2D.

Conclusions

Our findings provide strong evidence of shared genetics and causation between TV watching and T2D and facilitate our identification of common genetic architectures shared between them.

Association of prenatal acetaminophen use and acetaminophen metabolites with DNA methylation of newborns: analysis of two consecutive generations of the Isle of Wight birth cohort

Acetaminophen is used by nearly two-thirds of pregnant women. Although considered safe, studies have demonstrated associations between prenatal acetaminophen use and adverse health outcomes in offspring. Since DNA methylation (DNAm) at birth may act as an early indicator of later health, assessments on whether DNAm of newborns is associated with gestational acetaminophen use or its metabolites are needed. Using data from three consecutive generations of the Isle of Wight cohort (F0-grandmothers, F1-mothers, and F2-offspring) we investigated associations between acetaminophen metabolites in F0 serum at delivery with epigenome-wide DNAm in F1 (Guthrie cards) and between acetaminophen use of F1 and F2-cord-serum levels with F2 cord blood DNAm. In epigenome-wide screening, we eliminated non-informative DNAm sites followed by linear regression of informative sites. Based on repeated pregnancies, indication bias analyses tested whether acetaminophen indicated maternal diseases or has a risk in its own right. Considering that individuals with similar intake process acetaminophen differently, metabolites were clustered to distinguish metabolic exposures. Finally, metabolite clusters from F1-maternal and F2-cord sera were tested for their associations with newborn DNAm (F1 and F2). Twenty-one differential DNAm sites in cord blood were associated with reported maternal acetaminophen intake in the F2 generation. For 11 of these cytosine-phosphate-guanine (CpG) sites, an indication bias was excluded and five were replicated in F2 with metabolite clusters. In addition, metabolite clusters showed associations with 25 CpGs in the F0-F1 discovery analysis, of which five CpGs were replicated in the F2-generation. Our results suggest that prenatal acetaminophen use, measured as metabolites, may influence DNAm in newborns.

Deep Resequencing of 9 Candidate Genes Identifies a Role for ARAP1 and IGF2BP2 in Modulating Insulin Secretion Adjusted for Insulin Resistance in Obese Southern Europeans

Type 2 diabetes is characterized by impairment in insulin secretion, with an established genetic contribution. We aimed to evaluate common and low-frequency (1–5%) variants in nine genes strongly associated with insulin secretion by targeted sequencing in subjects selected from the extremes of insulin release measured by the disposition index. Collapsing data by gene and/or function, the association between disposition index and nonsense variants were significant, also after adjustment for confounding factors (OR = 0.25, 95% CI = 0.11–0.59, p = 0.001). Evaluating variants individually, three novel variants in ARAP1, IGF2BP2 and GCK, out of eight reaching significance singularly, remained associated after adjustment. Constructing a genetic risk model combining the effects of the three variants, only carriers of the ARAP1 and IGF2BP2 variants were significantly associated with a reduced probability to be in the lower, worst, extreme of insulin secretion (OR = 0.223, 95% CI = 0.105–0.473, p < 0.001). Observing a high number of normal glucose tolerance between carriers, a regression posthoc analysis was performed. Carriers of genetic risk model variants had higher probability to be normoglycemic, also after adjustment (OR = 2.411, 95% CI = 1.136–5.116, p = 0.022). Thus, in our southern European cohort, nonsense variants in all nine candidate genes showed association with better insulin secretion adjusted for insulin resistance, and we established the role of ARAP1 and IGF2BP2 in modulating insulin secretion.

Stem Cell-Derived β Cells: A Versatile Research Platform to Interrogate the Genetic Basis of β Cell Dysfunction

Pancreatic β cell dysfunction is a central component of diabetes progression. During the last decades, the genetic basis of several monogenic forms of diabetes has been recognized. Genome-wide association studies (GWAS) have also facilitated the identification of common genetic variants associated with an increased risk of diabetes. These studies highlight the importance of impaired β cell function in all forms of diabetes. However, how most of these risk variants confer disease risk, remains unanswered. Understanding the specific contribution of genetic variants and the precise role of their molecular effectors is the next step toward developing treatments that target β cell dysfunction in the era of personalized medicine. Protocols that allow derivation of β cells from pluripotent stem cells, represent a powerful research tool that allows modeling of human development and versatile experimental designs that can be used to shed some light on diabetes pathophysiology. This article reviews different models to study the genetic basis of β cell dysfunction, focusing on the recent advances made possible by stem cell applications in the field of diabetes research.

Identification of direct transcriptional targets of NFATC2 that promote β cell proliferation

The transcription factor NFATC2 induces β cell proliferation in mouse and human islets. However, the genomic targets that mediate these effects have not been identified. We expressed active forms of Nfatc2 and Nfatc1 in human islets. By integrating changes in gene expression with genomic binding sites for NFATC2, we identified approximately 2200 transcriptional targets of NFATC2. Genes induced by NFATC2 were enriched for transcripts that regulate the cell cycle and for DNA motifs associated with the transcription factor FOXP. Islets from an endocrine-specific Foxp1, Foxp2, and Foxp4 triple-knockout mouse were less responsive to NFATC2-induced β cell proliferation, suggesting the FOXP family works to regulate β cell proliferation in concert with NFATC2. NFATC2 induced β cell proliferation in both mouse and human islets, whereas NFATC1 did so only in human islets. Exploiting this species difference, we identified approximately 250 direct transcriptional targets of NFAT in human islets. This gene set enriches for cell cycle-associated transcripts and includes Nr4a1. Deletion of Nr4a1 reduced the capacity of NFATC2 to induce β cell proliferation, suggesting that much of the effect of NFATC2 occurs through its induction of Nr4a1. Integration of noncoding RNA expression, chromatin accessibility, and NFATC2 binding sites enabled us to identify NFATC2-dependent enhancer loci that mediate β cell proliferation.

Functional characterization of T2D-associated SNP effects on baseline and ER stress-responsive β cell transcriptional activation

Genome-wide association studies (GWAS) have linked single nucleotide polymorphisms (SNPs) at >250 loci in the human genome to type 2 diabetes (T2D) risk. For each locus, identifying the functional variant(s) among multiple SNPs in high linkage disequilibrium is critical to understand molecular mechanisms underlying T2D genetic risk. Using massively parallel reporter assays (MPRA), we test the cis-regulatory effects of SNPs associated with T2D and altered in vivo islet chromatin accessibility in MIN6 β cells under steady state and pathophysiologic endoplasmic reticulum (ER) stress conditions. We identify 1,982/6,621 (29.9%) SNP-containing elements that activate transcription in MIN6 and 879 SNP alleles that modulate MPRA activity. Multiple T2D-associated SNPs alter the activity of short interspersed nuclear element (SINE)-containing elements that are strongly induced by ER stress. We identify 220 functional variants at 104 T2D association signals, narrowing 54 signals to a single candidate SNP. Together, this study identifies elements driving β cell steady state and ER stress-responsive transcriptional activation, nominates causal T2D SNPs, and uncovers potential roles for repetitive elements in β cell transcriptional stress response and T2D genetics.

Inflammation drives alternative first exon usage to regulate immune genes including a novel iron-regulated isoform of Aim2

Determining the layers of gene regulation within the innate immune response is critical to our understanding of the cellular responses to infection and dysregulation in disease. We identified a conserved mechanism of gene regulation in human and mouse via changes in alternative first exon (AFE) usage following inflammation, resulting in changes to the isoforms produced. Of these AFE events, we identified 95 unannotated transcription start sites in mice using a de novo transcriptome generated by long-read native RNA-sequencing, one of which is in the cytosolic receptor for dsDNA and known inflammatory inducible gene, Aim2. We show that this unannotated AFE isoform of Aim2 is the predominant isoform expressed during inflammation and contains an iron-responsive element in its 5′UTR enabling mRNA translation to be regulated by iron levels. This work highlights the importance of examining alternative isoform changes and translational regulation in the innate immune response and uncovers novel regulatory mechanisms of Aim2.

Genome-wide discovery of genetic loci that uncouple excess adiposity from its comorbidities

Obesity is a major risk factor for cardiometabolic diseases. Nevertheless, a substantial proportion of individuals with obesity do not suffer cardiometabolic comorbidities. The mechanisms that uncouple adiposity from its cardiometabolic complications are not fully understood. Here, we identify 62 loci of which the same allele is significantly associated with both higher adiposity and lower cardiometabolic risk. Functional analyses show that the 62 loci are enriched for genes expressed in adipose tissue, and for regulatory variants that influence nearby genes that affect adipocyte differentiation. Genes prioritized in each locus support a key role of fat distribution (FAM13A, IRS1 and PPARG) and adipocyte function (ALDH2, CCDC92, DNAH10, ESR1, FAM13A, MTOR, PIK3R1 and VEGFB). Several additional mechanisms are involved as well, such as insulin-glucose signalling (ADCY5, ARAP1, CREBBP, FAM13A, MTOR, PEPD, RAC1 and SH2B3), energy expenditure and fatty acid oxidation (IGF2BP2), browning of white adipose tissue (CSK, VEGFA, VEGFB and SLC22A3) and inflammation (SH2B3, DAGLB and ADCY9). Some of these genes may represent therapeutic targets to reduce cardiometabolic risk linked to excess adiposity.

LncRNA ARAP1-AS2 promotes high glucose-induced human proximal tubular cell injury via persistent transactivation of the EGFR by interacting with ARAP1

The persistent transactivation of epidermal growth factor receptor (EGFR) causes subsequent activation of the TGF-β/Smad3 pathway, which is closely associated with fibrosis and cell proliferation in diabetic nephropathy (DN), but the exact mechanism of persistent EGFR transactivation in DN remains unclear. ARAP1, a susceptibility gene for type 2 diabetes, can regulate the endocytosis and ubiquitination of membrane receptors, but the effect of ARAP1 and its natural antisense long non-coding RNA (lncRNA), ARAP1-AS2, on the ubiquitination of EGFR in DN is not clear. In this study, we verified that the expression of ARAP1 and ARAP1-AS2 was significantly up-regulated in high glucose-induced human proximal tubular epithelial cells (HK-2 cells). Moreover, we found that overexpression or knockdown of ARAP1-AS2 could regulate fibrosis and HK-2 cell proliferation through EGFR/TGF-β/Smad3 signalling. RNA pulldown assays revealed that ARAP1-AS2 directly interacts with ARAP1. Coimmunoprecipitation, dual-immunofluorescence and ubiquitination assays showed that ARAP1 may maintain persistent EGFR activation by reducing EGFR ubiquitination through competing with Cbl for CIN85 binding. Taken together, our results suggest that the lncRNA ARAP1-AS2 may promote high glucose-induced proximal tubular cell injury via persistent EGFR/TGF-β/Smad3 pathway activation by interacting with ARAP1.

Pancreatic β-cells in type 1 and type 2 diabetes mellitus: different pathways to failure

Loss of functional β-cell mass is the key mechanism leading to the two main forms of diabetes mellitus - type 1 diabetes mellitus (T1DM) and type 2 diabetes mellitus (T2DM). Understanding the mechanisms behind β-cell failure is critical to prevent or revert disease. Basic pathogenic differences exist in the two forms of diabetes mellitus; T1DM is immune mediated and T2DM is mediated by metabolic mechanisms. These mechanisms differentially affect early β-cell dysfunction and eventual fate. Over the past decade, major advances have been made in the field, mostly delivered by studies on β-cells in human disease. These advances include studies of islet morphology and human β-cell gene expression in T1DM and T2DM, the identification and characterization of the role of T1DM and T2DM candidate genes at the β-cell level and the endoplasmic reticulum stress signalling that contributes to β-cell failure in T1DM (mostly IRE1 driven) and T2DM (mostly PERK-eIF2α dependent). Here, we review these new findings, focusing on studies performed on human β-cells or on samples obtained from patients with diabetes mellitus.

Maya gene variants related to the risk of type 2 diabetes in a family-based association study

Mexican Maya populations have a notably high prevalence of type 2 diabetes (T2D) as a consequence of the interaction between environmental factors and a genetic component. To assess the impact of 24 single nucleotide variants (SNVs) located in 18 T2D risk genes, we conducted a family-based association evaluation in samples from Maya communities with a high incidence of the disease. A total of four hundred individuals were recruited from three Maya communities with a high T2D incidence. Family pedigrees (100) and 49 nuclear families were included. Genotyping was performed by allelic discrimination with TaqMan probes. This study also included the family-based association test (FBAT) statistic U to assess the genetic associations with T2D, and the multivariate statistical and haplotype analyses. A positive association with TD2 risk was found for WFS1 rs6446482 (p = 0.046, Z = 1.994) under an additive model, and SIRT1 rs7896005 (p = 0.038, Z = 2.073) under the dominant model. Multivariate model analysis, including T2D status, age, and body mass index (BMI), displayed significant covariance in PPARGC-1α rs8192678; SIRT1 rs7896005; TCF7L2 rs7903146 and rs122243326; UCP3 rs3781907; and HHEX rs1111875 with a P < 0.05. This study revealed an association of SIRT1 and WFS1 with T2D risk.

From Genetic Association to Molecular Mechanisms for Islet-cell Dysfunction in Type 2 Diabetes

Genome-wide association studies (GWAS) have identified over 400 signals robustly associated with risk for type 2 diabetes (T2D). At the vast majority of these loci, the lead single nucleotide polymorphisms (SNPs) reside in noncoding regions of the genome, which hampers biological inference and translation of genetic discoveries into disease mechanisms. The study of these T2D risk variants in normoglycemic individuals has revealed that a significant proportion are exerting their disease risk through islet-cell dysfunction. The central role of the islet is also demonstrated by numerous studies, which have shown an enrichment of these signals in islet-specific epigenomic annotations. In recent years the emergence of authentic human beta-cell lines, and advances in genome-editing technologies coupled with improved protocols differentiating human pluripotent stem cells into beta-like cells has opened up new opportunities for T2D disease modeling. Here we review the current understanding on the genetic basis of T2D focusing on approaches, which have facilitated the identification of causal variants and their effector transcripts in human islets. We will present examples of functional studies based on animal and conventional cellular systems and highlight the potential of novel stem cell-based T2D disease models.

The effect of lncRNA-ARAP1-AS2/ARAP1 on high glucose-induced cytoskeleton rearrangement and epithelial-mesenchymal transition in human renal tubular epithelial cells

The epithelial-mesenchymal transition (EMT) plays an important role in diabetic renal fibrosis. The ARAP1 gene is located near risk alleles for Type 2 diabetes, and its function has been linked to cytoskeleton rearrangement, Golgi apparatus remodeling, and endocytic trafficking of membrane receptors. The role of ARAP1 and its antisense RNA, ARAP1-AS2, in the pathogenesis of diabetes is unclear. To clarify the roles of ARAP1 and its antisense RNA in diabetes and related complications, we examined if the expression of these transcripts changed under high glucose (HG) conditions. To do this, we examined transcript levels in HK-2 cells, and explored the roles of ARAP1 and ARAP1-AS2 in the EMT process in HK-2 cells. We found increased expression of ARAP1-AS2 and ARAP1 in HK-2 cells under HG condition, and observed that the overexpression of ARAP1-AS2 significantly increased the EMT process. In addition, HG upregulated Cdc42-GTP levels in HK-2 cells, and increased cytoskeleton rearrangement, cell viability, and migration. After knockdown of ARAP1, the level of Cdc42-GTP was decreased; cytoskeleton reorganization, cell viability, and migration processes were decreased; and EMT and expression of fibrosis marker protein. Overall, our results indicated that ARAP1-AS2/ARAP1 may participate in cytoskeleton rearrangement and EMT processes in HK-2 cells through increased Cdc42-GTP levels.

Identifying the association between single nucleotide polymorphisms in KCNQ1, ARAP1, and KCNJ11 and type 2 diabetes mellitus in a Chinese population

Background: Type 2 diabetes mellitus (T2DM) has a high global prevalence, and insufficient insulin secretion is one of the major reasons for its development. Therefore, investigating the association between T2DM and the single nucleotide polymorphisms (SNPs) in genes associated with insulin secretion is necessary. Methods: T2DM (1,194) and nondiabetic (NDM) (1,292) subjects were enrolled and the ten single nucleotide polymorphisms (SNPs) in KCNQ1, ARAP1, and KCNJ11 associated with insulin secretion were genotyped in a Chinese population. Results: Our data revealed that the rs2237897T allele in KCNQ1 is the protective allele for T2DM (P<0.001, OR=0.793; 95%CI: 0.705-0.893). However, the A allele of rs1552224 in ARAP1 may be a risk factor for T2DM (P=0.002, OR=12.070; 95% CI: 1.578-92.337). The haplotype analysis revealed that rs151290-rs2237892CC and rs2237895-rs2237897CC in KCNQ1 constitute the risk haplotype in T2DM development (P=0.010, OR=1.160; 95% CI: 1.037-1.299 and P=0.004, OR=1.192; 95% CI: 1.057-1.344). Moreover, rs2237895-rs2237897AT in KCNQ1 constitutes the protective haplotype in T2DM (P=0.001, OR=0.819; 95% CI: 0.727-0.923). In the inheritance models analysis, the rs2283228 (C/A-C/C) genotype is the protective factor compared to the A/A genotype (P=0.005, OR=0.79; 95% CI: 0.68-0.93). For rs2237897, the C/T-T/T genotype is the protective factor compared to the C/C genotype (P<0.001, OR=0.74; 95% CI: 0.63-0.87). Furthermore, when compared with the rs2237897 (C/T-T/T) genotype, rs2237897C/C genotype showed higher HbA1C levels (8.731±2.697 vs 9.282±2.921, P=0.001). Conclusion: Our results revealed that genetic variations in KCNQ1 and ARAP1 were associated with T2DM susceptibility in a Chinese population.

PAX proteins and their role in pancreas

Gene regulatory factors that govern the expression of heritable information come in an array of flavors, chiefly with transcription factors, the proteins which bind to regions of specific genes and modulate gene transcription, subsequently altering cellular function. PAX transcription factors are sequence-specific DNA-binding proteins exerting its regulatory activity in many tissues. Notably, three members of the PAX family namely PAX2, PAX4 and PAX6 have emerged as crucial players at multiple steps of pancreatic development and differentiation and also play a pivotal role in the regulation of pancreatic islet hormones synthesis and secretion. Providing a comprehensive outline of these transcription factors and their primordial and divergent roles in the pancreas is far-reaching in contemporary diabetes research. Accordingly, this review furnishes an outline of the role of pancreatic specific PAX regulators in the development of the pancreas and its associated disorders.

Open Chromatin Profiling in Adipose Tissue Marks Genomic Regions with Functional Roles in Cardiometabolic Traits

Identifying the regulatory mechanisms of genome-wide association study (GWAS) loci affecting adipose tissue has been restricted due to limited characterization of adipose transcriptional regulatory elements. We profiled chromatin accessibility in three frozen human subcutaneous adipose tissue needle biopsies and preadipocytes and adipocytes from the Simpson Golabi-Behmel Syndrome (SGBS) cell strain using an assay for transposase-accessible chromatin (ATAC-seq). We identified 68,571 representative accessible chromatin regions (peaks) across adipose tissue samples (FDR < 5%). GWAS loci for eight cardiometabolic traits were enriched in these peaks (P < 0.005), with the strongest enrichment for waist-hip ratio. Of 110 recently described cardiometabolic GWAS loci colocalized with adipose tissue eQTLs, 59 loci had one or more variants overlapping an adipose tissue peak. Annotated variants at the SNX10 waist-hip ratio locus and the ATP2A1-SH2B1 body mass index locus showed allelic differences in regulatory assays. These adipose tissue accessible chromatin regions elucidate genetic variants that may alter adipose tissue function to impact cardiometabolic traits.

Functional characteristics of novel pancreatic Pax6 regulatory elements

Complex diseases, such as diabetes, are influenced by comprehensive transcriptional networks. Genome-wide association studies have revealed that variants located in regulatory elements for pancreatic transcription factors are linked to diabetes, including those functionally linked to the paired box transcription factor Pax6. Pax6 deletions in adult mice cause rapid onset of classic diabetes, but the full spectrum of pancreatic Pax6 regulators is unknown. Using a regulatory element discovery approach, we identified two novel Pax6 pancreatic cis-regulatory elements in a poorly characterized regulatory desert. Both new elements, Pax6 pancreas cis-regulatory element 3 (PE3) and PE4, are located 50 and 100 kb upstream and interact with different parts of the Pax6 promoter and nearby non-coding RNAs. They drive expression in the developing pancreas and brain and code for multiple pancreas-related transcription factor-binding sites. PE3 binds CCCTC-binding factor (CTCF) and is marked by stem cell identity markers in embryonic stem cells, whilst a common variant located in the PE4 element affects binding of Pax4, a known pancreatic regulator, altering Pax6 gene expression. To determine the ability of these elements to regulate gene expression, synthetic transcriptional activators and repressors were targeted to PE3 and PE4, modulating Pax6 gene expression, as well as influencing neighbouring genes and long non-coding RNAs, implicating the Pax6 locus in pancreas function and diabetes.

Widespread Enhancer Activity from Core Promoters

Gene expression in higher eukaryotes is precisely regulated in time and space through the interplay between promoters and gene-distal regulatory regions, known as enhancers. The original definition of enhancers implies the ability to activate gene expression remotely, while promoters entail the capability to locally induce gene expression. Despite the conventional distinction between them, promoters and enhancers share many genomic and epigenomic features. One intriguing finding in the gene regulation field comes from the observation that many core promoter regions display enhancer activity. Recent high-throughput reporter assays along with clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9-related approaches have indicated that this phenomenon is common and might have a strong impact on our global understanding of genome organisation and gene expression regulation.

Refining the accuracy of validated target identification through coding variant fine-mapping in type 2 diabetes

We aggregated coding variant data for 81,412 type 2 diabetes cases and 370,832 controls of diverse ancestry, identifying 40 coding variant association signals (P < 2.2 × 10-7); of these, 16 map outside known risk-associated loci. We make two important observations. First, only five of these signals are driven by low-frequency variants: even for these, effect sizes are modest (odds ratio ≤1.29). Second, when we used large-scale genome-wide association data to fine-map the associated variants in their regional context, accounting for the global enrichment of complex trait associations in coding sequence, compelling evidence for coding variant causality was obtained for only 16 signals. At 13 others, the associated coding variants clearly represent 'false leads' with potential to generate erroneous mechanistic inference. Coding variant associations offer a direct route to biological insight for complex diseases and identification of validated therapeutic targets; however, appropriate mechanistic inference requires careful specification of their causal contribution to disease predisposition.

Integrating ChIP-seq with other functional genomics data

Transcription is regulated by transcription factor (TF) binding at promoters and distal regulatory elements and histone modifications that control the accessibility of these elements. Chromatin immunoprecipitation followed by sequencing (ChIP-seq) has become the standard assay for identifying genome-wide protein-DNA interactions in vitro and in vivo. As large-scale ChIP-seq data sets have been collected for different TFs and histone modifications, their potential to predict gene expression can be used to test hypotheses about the mechanisms of gene regulation. In addition, complementary functional genomics assays provide a global view of chromatin accessibility and long-range cis-regulatory interactions that are being combined with TF binding and histone remodeling to study the regulation of gene expression. Thus, ChIP-seq analysis is now widely integrated with other functional genomics assays to better understand gene regulatory mechanisms. In this review, we discuss advances and challenges in integrating ChIP-seq data to identify context-specific chromatin states associated with gene activity. We describe the overall computational design of integrating ChIP-seq data with other functional genomics assays. We also discuss the challenges of extending these methods to low-input ChIP-seq assays and related single-cell assays.

Applications of the 1000 Genomes Project resources

The 1000 Genomes Project created a valuable, worldwide reference for human genetic variation. Common uses of the 1000 Genomes dataset include genotype imputation supporting Genome-wide Association Studies, mapping expression Quantitative Trait Loci, filtering non-pathogenic variants from exome, whole genome and cancer genome sequencing projects, and genetic analysis of population structure and molecular evolution. In this article, we will highlight some of the multiple ways that the 1000 Genomes data can be and has been utilized for genetic studies.

Impact of Genetic Variants on the Individual Potential for Body Fat Loss

The past decade has witnessed the discovery of obesity-related genetic variants and their functions through genome-wide association studies. Combinations of risk alleles can influence obesity phenotypes with different degrees of effectiveness across various individuals by interacting with environmental factors. We examined the interaction between genetic variation and changes in dietary habits or exercise that influences body fat loss from a large Korean cohort (n = 8840). Out of 673 obesity-related SNPs, a total of 100 SNPs (37 for carbohydrate intake; 19 for fat intake; 44 for total calories intake; 25 for exercise onset) identified to have gene-environment interaction effect in generalized linear model were used to calculate genetic risk scores (GRS). Based on the GRS distribution, we divided the population into four levels, namely, “very insensitive”, “insensitive”, “sensitive”, and “very sensitive” for each of the four categories, “carbohydrate intake”, “fat intake”, “total calories intake”, and “exercise”. Overall, the mean body fat loss became larger when the sensitivity level was increased. In conclusion, genetic variants influence the effectiveness of dietary regimes for body fat loss. Based on our findings, we suggest a platform for personalized body fat management by providing the most suitable and effective nutrition or activity plan specific to an individual.

Functional Variants Identified Efficiently through an Integrated Transcriptome and Epigenome Analysis

Although genome-wide association studies (GWAS) have identified numerous genetic loci associated with complex diseases, the underlying molecular mechanisms of how these loci contribute to disease pathogenesis remain largely unknown, due to the lack of an efficient strategy to identify these risk variants. Here, we proposed a new strategy termed integrated transcriptome and epigenome analysis (iTEA) to identify functional genetic variants in non-coding elements. We considered type 2 diabetes mellitus as a model and identified a well-known diabetic risk variant rs35767 using iTEA. Furthermore, we discovered a new functional SNP, rs815815, involved in glucose metabolism. Our study provides an approach to directly and quickly identify functional genetic variants in type 2 diabetes mellitus, and this approach can be extended to study other complex diseases.

Genetic aetiology of glycaemic traits: approaches and insights

Glycaemic traits such as fasting and post-challenge glucose and insulin measures, as well as glycated haemoglobin (HbA1c), are used to diagnose and monitor diabetes. These traits are risk factors for cardiovascular disease even below the diabetic threshold, and their study can additionally yield insights into the pathophysiology of type 2 diabetes. To date, a diverse set of genetic approaches have led to the discovery of over 97 loci influencing glycaemic traits. In this review, we will focus on recent advances in the genetic aetiology of glycaemic traits, and the resulting biological insights. We will provide a brief overview of results ranging from common, to low- and rare-frequency variant-trait association studies, studies leveraging the diversity across populations, and studies harnessing the power of genetic and genomic approaches to gain insights into the biological underpinnings of these traits.

Trans-ancestry Fine Mapping and Molecular Assays Identify Regulatory Variants at the ANGPTL8 HDL-C GWAS Locus

Recent genome-wide association studies (GWAS) have identified variants associated with high-density lipoprotein cholesterol (HDL-C) located in or near the ANGPTL8 gene. Given the extensive sharing of GWAS loci across populations, we hypothesized that at least one shared variant at this locus affects HDL-C. The HDL-C–associated variants are coincident with expression quantitative trait loci for ANGPTL8 and DOCK6 in subcutaneous adipose tissue; however, only ANGPTL8 expression levels are associated with HDL-C levels. We identified a 400-bp promoter region of ANGPTL8 and enhancer regions within 5 kb that contribute to regulating expression in liver and adipose. To identify variants functionally responsible for the HDL-C association, we performed fine-mapping analyses and selected 13 candidate variants that overlap putative regulatory regions to test for allelic differences in regulatory function. Of these variants, rs12463177-G increased transcriptional activity (1.5-fold, P = 0.004) and showed differential protein binding. Six additional variants (rs17699089, rs200788077, rs56322906, rs3760782, rs737337, and rs3745683) showed evidence of allelic differences in transcriptional activity and/or protein binding. Taken together, these data suggest a regulatory mechanism at the ANGPTL8 HDL-C GWAS locus involving tissue-selective expression and at least one functional variant.

A Type 2 Diabetes-Associated Functional Regulatory Variant in a Pancreatic Islet Enhancer at the ADCY5 Locus

Molecular mechanisms remain unknown for most type 2 diabetes genome-wide association study identified loci. Variants associated with type 2 diabetes and fasting glucose levels reside in introns of ADCY5, a gene that encodes adenylate cyclase 5. Adenylate cyclase 5 catalyzes the production of cyclic AMP, which is a second messenger molecule involved in cell signaling and pancreatic β-cell insulin secretion. We demonstrated that type 2 diabetes risk alleles are associated with decreased ADCY5 expression in human islets and examined candidate variants for regulatory function. rs11708067 overlaps a predicted enhancer region in pancreatic islets. The type 2 diabetes risk rs11708067-A allele showed fewer H3K27ac ChIP-seq reads in human islets, lower transcriptional activity in reporter assays in rodent β-cells (rat 832/13 and mouse MIN6), and increased nuclear protein binding compared with the rs11708067-G allele. Homozygous deletion of the orthologous enhancer region in 832/13 cells resulted in a 64% reduction in expression level of Adcy5, but not adjacent gene Sec22a, and a 39% reduction in insulin secretion. Together, these data suggest that rs11708067-A risk allele contributes to type 2 diabetes by disrupting an islet enhancer, which results in reduced ADCY5 expression and impaired insulin secretion.

Prioritising Causal Genes at Type 2 Diabetes Risk Loci

PURPOSE OF REVIEW

Genome-wide association studies (GWAS) for type 2 diabetes (T2D) risk have identified a large number of genetic loci associated with disease susceptibility. However, progress moving from association signals through causal genes to functional understanding has so far been slow, hindering clinical translation. This review discusses the benefits and limitations of emerging, unbiased approaches for prioritising causal genes at T2D risk loci.

RECENT FINDINGS

Candidate causal genes can be identified by a number of different strategies that rely on genetic data, genomic annotations, and functional screening of selected genes. To overcome the limitations of each particular method, integration of multiple data sets is proving essential for establishing confidence in the prioritised genes. Previous studies have also highlighted the need to support these efforts through identification of causal variants and disease-relevant tissues. Prioritisation of causal genes at T2D risk loci by integrating complementary lines of evidence promises to accelerate our understanding of disease pathology and promote translation into new therapeutics.

Mechanisms of Type 2 Diabetes Risk Loci

PURPOSE OF REVIEW

Deciphering the mechanisms of type 2 diabetes (T2DM) risk loci can greatly inform on disease pathology. This review discusses current knowledge of mechanisms through which genetic variants influence T2DM risk and considerations for future studies.

RECENT FINDINGS

Over 100 T2DM risk loci to date have been identified. Candidate causal variants at risk loci map predominantly to non-coding sequence. Physiological, epigenomic and gene expression data suggest that variants at many known T2DM risk loci affect pancreatic islet regulation, although variants at other loci also affect protein function and regulatory processes in adipose, pre-adipose, liver, skeletal muscle and brain. The effects of T2DM variants on regulatory activity in these tissues appear largely, but not exclusively, due to altered transcription factor binding. Putative target genes of T2DM variants have been defined at an increasing number of loci and some, such as FTO, may entail several genes and multiple tissues. Gene networks in islets and adipocytes have been implicated in T2DM risk, although the molecular pathways of risk genes remain largely undefined. Efforts to fully define the mechanisms of T2DM risk loci are just beginning. Continued identification of risk mechanisms will benefit from combining genetic fine-mapping with detailed phenotypic association data, high-throughput epigenomics data from diabetes-relevant tissue, functional screening of candidate genes and genome editing of cellular and animal models.

Genomics of Islet (Dys)function and Type 2 Diabetes

Pancreatic islet dysfunction and beta cell failure are hallmarks of type 2 diabetes mellitus (T2DM) pathogenesis. In this review, we discuss how genome-wide association studies (GWASs) and recent developments in islet (epi)genome and transcriptome profiling (particularly single cell analyses) are providing novel insights into the genetic, environmental, and cellular contributions to islet (dys)function and T2DM pathogenesis. Moving forward, study designs that interrogate and model genetic variation [e.g., allelic profiling and (epi)genome editing] will be critical to dissect the molecular genetics of T2DM pathogenesis, to build next-generation cellular and animal models, and to develop precision medicine approaches to detect, treat, and prevent islet (dys)function and T2DM.

Genetic regulatory signatures underlying islet gene expression and type 2 diabetes

Genome-wide association studies (GWAS) have identified >100 independent SNPs that modulate the risk of type 2 diabetes (T2D) and related traits. However, the pathogenic mechanisms of most of these SNPs remain elusive. Here, we examined genomic, epigenomic, and transcriptomic profiles in human pancreatic islets to understand the links between genetic variation, chromatin landscape, and gene expression in the context of T2D. We first integrated genome and transcriptome variation across 112 islet samples to produce dense cis-expression quantitative trait loci (cis-eQTL) maps. Additional integration with chromatin-state maps for islets and other diverse tissue types revealed that cis-eQTLs for islet-specific genes are specifically and significantly enriched in islet stretch enhancers. High-resolution chromatin accessibility profiling using assay for transposase-accessible chromatin sequencing (ATAC-seq) in two islet samples enabled us to identify specific transcription factor (TF) footprints embedded in active regulatory elements, which are highly enriched for islet cis-eQTL. Aggregate allelic bias signatures in TF footprints enabled us de novo to reconstruct TF binding affinities genetically, which support the high-quality nature of the TF footprint predictions. Interestingly, we found that T2D GWAS loci were strikingly and specifically enriched in islet Regulatory Factor X (RFX) footprints. Remarkably, within and across independent loci, T2D risk alleles that overlap with RFX footprints uniformly disrupt the RFX motifs at high-information content positions. Together, these results suggest that common regulatory variations have shaped islet TF footprints and the transcriptome and that a confluent RFX regulatory grammar plays a significant role in the genetic component of T2D predisposition.

Decreased STARD10 Expression Is Associated with Defective Insulin Secretion in Humans and Mice

Genetic variants near ARAP1 (CENTD2) and STARD10 influence type 2 diabetes (T2D) risk. The risk alleles impair glucose-induced insulin secretion and, paradoxically but characteristically, are associated with decreased proinsulin:insulin ratios, indicating improved proinsulin conversion. Neither the identity of the causal variants nor the gene(s) through which risk is conferred have been firmly established. Whereas ARAP1 encodes a GTPase activating protein, STARD10 is a member of the steroidogenic acute regulatory protein (StAR)-related lipid transfer protein family. By integrating genetic fine-mapping and epigenomic annotation data and performing promoter-reporter and chromatin conformational capture (3C) studies in β cell lines, we localize the causal variant(s) at this locus to a 5 kb region that overlaps a stretch-enhancer active in islets. This region contains several highly correlated T2D-risk variants, including the rs140130268 indel. Expression QTL analysis of islet transcriptomes from three independent subject groups demonstrated that T2D-risk allele carriers displayed reduced levels of STARD10 mRNA, with no concomitant change in ARAP1 mRNA levels. Correspondingly, β-cell-selective deletion of StarD10 in mice led to impaired glucose-stimulated Ca²⁺ dynamics and insulin secretion and recapitulated the pattern of improved proinsulin processing observed at the human GWAS signal. Conversely, overexpression of StarD10 in the adult β cell improved glucose tolerance in high fat-fed animals. In contrast, manipulation of Arap1 in β cells had no impact on insulin secretion or proinsulin conversion in mice. This convergence of human and murine data provides compelling evidence that the T2D risk associated with variation at this locus is mediated through reduction in STARD10 expression in the β cell.

The Metabolic Syndrome in Men study: a resource for studies of metabolic and cardiovascular diseases

The Metabolic Syndrome in Men (METSIM) study is a population-based study including 10,197 Finnish men examined in 2005-2010. The aim of the study is to investigate nongenetic and genetic factors associated with the risk of T2D and CVD, and with cardiovascular risk factors. The protocol includes a detailed phenotyping of the participants, an oral glucose tolerance test, fasting laboratory measurements including proton NMR measurements, mass spectometry metabolomics, adipose tissue biopsies from 1,400 participants, and a stool sample. In our ongoing follow-up study, we have, to date, reexamined 6,496 participants. Extensive genotyping and exome sequencing have been performed for essentially all METSIM participants, and >2,000 METSIM participants have been whole-genome sequenced. We have identified several nongenetic markers associated with the development of diabetes and cardiovascular events, and participated in several genetic association studies to identify gene variants associated with diabetes, hyperglycemia, and cardiovascular risk factors. The generation of a phenotype and genotype resource in the METSIM study allows us to proceed toward a "systems genetics" approach, which includes statistical methods to quantitate and integrate intermediate phenotypes, such as transcript, protein, or metabolite levels, to provide a global view of the molecular architecture of complex traits.

The Genetics of Transcription Factor DNA Binding Variation

Most complex trait-associated variants are located in non-coding regulatory regions of the genome, where they have been shown to disrupt transcription factor (TF)-DNA binding motifs. Variable TF-DNA interactions are therefore increasingly considered as key drivers of phenotypic variation. However, recent genome-wide studies revealed that the majority of variable TF-DNA binding events are not driven by sequence alterations in the motif of the studied TF. This observation implies that the molecular mechanisms underlying TF-DNA binding variation and, by extrapolation, inter-individual phenotypic variation are more complex than originally anticipated. Here, we summarize the findings that led to this important paradigm shift and review proposed mechanisms for local, proximal, or distal genetic variation-driven variable TF-DNA binding. In addition, we discuss the biomedical implications of these findings for our ability to dissect the molecular role(s) of non-coding genetic variants in complex traits, including disease susceptibility.

Single-cell transcriptomes identify human islet cell signatures and reveal cell-type-specific expression changes in type 2 diabetes

Blood glucose levels are tightly controlled by the coordinated action of at least four cell types constituting pancreatic islets. Changes in the proportion and/or function of these cells are associated with genetic and molecular pathophysiology of monogenic, type 1, and type 2 (T2D) diabetes. Cellular heterogeneity impedes precise understanding of the molecular components of each islet cell type that govern islet (dys)function, particularly the less abundant delta and gamma/pancreatic polypeptide (PP) cells. Here, we report single-cell transcriptomes for 638 cells from nondiabetic (ND) and T2D human islet samples. Analyses of ND single-cell transcriptomes identified distinct alpha, beta, delta, and PP/gamma cell-type signatures. Genes linked to rare and common forms of islet dysfunction and diabetes were expressed in the delta and PP/gamma cell types. Moreover, this study revealed that delta cells specifically express receptors that receive and coordinate systemic cues from the leptin, ghrelin, and dopamine signaling pathways implicating them as integrators of central and peripheral metabolic signals into the pancreatic islet. Finally, single-cell transcriptome profiling revealed genes differentially regulated between T2D and ND alpha, beta, and delta cells that were undetectable in paired whole islet analyses. This study thus identifies fundamental cell-type–specific features of pancreatic islet (dys)function and provides a critical resource for comprehensive understanding of islet biology and diabetes pathogenesis.

Targeting pancreatic expressed PAX genes for the treatment of diabetes mellitus and pancreatic neuroendocrine tumors

INTRODUCTION

Four members of the PAX family, PAX2, PAX4, PAX6 and PAX8 are known to be expressed in the pancreas. Accumulated evidences indicate that several pancreatic expressed PAX genes play a significant role in pancreatic development/functionality and alterations in these genes are involved in the pathogenesis of pancreatic diseases. Areas covered: In this review, we summarize the ongoing research related to pancreatic PAX genes in diabetes mellitus and pancreatic neuroendocrine tumors. We dissect the current knowledge at different levels; from mechanistic studies in cell lines performed to understand the molecular processes controlled by pancreatic PAX genes, to in vivo studies using rodent models that over-express or lack specific PAX genes. Finally, we describe human studies associating variants on pancreatic-expressed PAX genes with pancreatic diseases. Expert opinion: Based on the current literature, we propose that future interventions to treat pancreatic neuroendocrine tumors and diabetes mellitus could be developed via the modulation of PAX4 and/or PAX6 regulated pathways.

Heritability and Genome-Wide Association Analyses of Sleep Duration in Children: The EAGLE Consortium

STUDY OBJECTIVES

Low or excessive sleep duration has been associated with multiple outcomes, but the biology behind these associations remains elusive. Specifically, genetic studies in children are scarce. In this study, we aimed to: (1) estimate the proportion of genetic variance of sleep duration in children attributed to common single nucleotide polymorphisms (SNPs), (2) identify novel SNPs associated with sleep duration in children, and (3) investigate the genetic overlap of sleep duration in children and related metabolic and psychiatric traits.

METHODS

We performed a population-based molecular genetic study, using data form the EArly Genetics and Life course Epidemiology (EAGLE) Consortium. 10,554 children of European ancestry were included in the discovery, and 1,250 children in the replication phase.

RESULTS

We found evidence of significant but modest SNP heritability of sleep duration in children (SNP h² 0.14, 95% CI [0.05, 0.23]) using the LD score regression method. A novel region at chromosome 11q13.4 (top SNP: rs74506765, P = 2.27e-08) was associated with sleep duration in children, but this was not replicated in independent studies. Nominally significant genetic overlap was only found (rG = 0.23, P = 0.05) between sleep duration in children and type 2 diabetes in adults, supporting the hypothesis of a common pathogenic mechanism.

CONCLUSIONS

The significant SNP heritability of sleep duration in children and the suggestive genetic overlap with type 2 diabetes support the search for genetic mechanisms linking sleep duration in children to multiple outcomes in health and disease.

Type 2 diabetes: genetic data sharing to advance complex disease research

As with other complex diseases, unbiased association studies followed by physiological and experimental characterization have for years formed a paradigm for identifying genes or processes of relevance to type 2 diabetes mellitus (T2D). Recent large-scale common and rare variant genome-wide association studies (GWAS) suggest that substantially larger association studies are needed to identify most T2D loci in the population. To hasten clinical translation of genetic discoveries, new paradigms are also required to aid specialized investigation of nascent hypotheses. We argue for an integrated T2D knowledgebase, designed for a worldwide community to access aggregated large-scale genetic data sets, as one paradigm to catalyse convergence of these efforts.

Multiple Hepatic Regulatory Variants at the GALNT2 GWAS Locus Associated with High-Density Lipoprotein Cholesterol

Genome-wide association studies (GWASs) have identified more than 150 loci associated with blood lipid and cholesterol levels; however, the functional and molecular mechanisms for many associations are unknown. We examined the functional regulatory effects of candidate variants at the GALNT2 locus associated with high-density lipoprotein cholesterol (HDL-C). Fine-mapping and conditional analyses in the METSIM study identified a single locus harboring 25 noncoding variants (r(2) > 0.7 with the lead GWAS variants) strongly associated with total cholesterol in medium-sized HDL (e.g., rs17315646, p = 3.5 × 10(-12)). We used luciferase reporter assays in HepG2 cells to test all 25 variants for allelic differences in regulatory enhancer activity. rs2281721 showed allelic differences in transcriptional activity (75-fold [T] versus 27-fold [C] more than the empty-vector control), as did a separate 780-bp segment containing rs4846913, rs2144300, and rs6143660 (49-fold [AT(-) haplotype] versus 16-fold [CC(+) haplotype] more). Using electrophoretic mobility shift assays, we observed differential CEBPB binding to rs4846913, and we confirmed this binding in a native chromatin context by performing chromatin-immunoprecipitation (ChIP) assays in HepG2 and Huh-7 cell lines of differing genotypes. Additionally, sequence reads in HepG2 DNase-I-hypersensitivity and CEBPB ChIP-seq signals spanning rs4846913 showed significant allelic imbalance. Allelic-expression-imbalance assays performed with RNA from primary human hepatocyte samples and expression-quantitative-trait-locus (eQTL) data in human subcutaneous adipose tissue samples confirmed that alleles associated with increased HDL-C are associated with a modest increase in GALNT2 expression. Together, these data suggest that at least rs4846913 and rs2281721 play key roles in influencing GALNT2 expression at this HDL-C locus.