A Common Type 2 Diabetes Risk Variant Potentiates Activity of an Evolutionarily Conserved Islet Stretch Enhancer and Increases C2CD4A and C2CD4B Expression

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.

C2CD4B Evokes Oxidative Stress and Vascular Dysfunction via a PI3K/Akt/PKCα-Signaling Pathway

High glucose-induced endothelial dysfunction is an important pathological feature of diabetic vasculopathy. While genome-wide studies have identified an association between type 2 diabetes mellitus (T2DM) and increased expression of a C2 calcium-dependent domain containing 4B (C2CD4B), no study has yet explored the possible direct effect of C2CD4B on vascular function. Vascular reactivity studies were conducted using a pressure myograph, and nitric oxide and oxidative stress were assessed through difluorofluorescein diacetate and dihydroethidium, respectively. We demonstrate that high glucose upregulated both mRNA and protein expression of C2CD4B in mice mesenteric arteries in a time-dependent manner. Notably, the inhibition of C2CD4B expression by genetic knockdown efficiently prevented hyperglycemia-induced oxidative stress, endothelial dysfunction, and loss of nitric oxide (NO) bioavailability. Recombinant C2CD4B evoked endothelial dysfunction of mice mesenteric arteries, an effect associated with increased reactive oxygen species (ROS) and decreased NO production. In isolated human umbilical vein endothelial cells (HUVECs), C2CD4B increased phosphorylation of endothelial nitric oxide synthase (eNOS) at the inhibitory site Thr495 and reduced eNOS dimerization. Pharmacological inhibitors of phosphoinositide 3-kinase (PI3K), Akt, and PKCα effectively attenuated oxidative stress, NO reduction, impairment of endothelial function, and eNOS uncoupling induced by C2CD4B. These data demonstrate, for the first time, that C2CD4B exerts a direct effect on vascular endothelium via a phosphoinositide 3-kinase (PI3K)/Akt/PKCα-signaling pathway, providing a new perspective on C2CD4B as a promising therapeutic target for the prevention of oxidative stress in diabetes-induced endothelial dysfunction.

Population-scale skeletal muscle single-nucleus multi-omic profiling reveals extensive context specific genetic regulation

Skeletal muscle, the largest human organ by weight, is relevant to several polygenic metabolic traits and diseases including type 2 diabetes (T2D). Identifying genetic mechanisms underlying these traits requires pinpointing the relevant cell types, regulatory elements, target genes, and causal variants. Here, we used genetic multiplexing to generate population-scale single nucleus (sn) chromatin accessibility (snATAC-seq) and transcriptome (snRNA-seq) maps across 287 frozen human skeletal muscle biopsies representing 456,880 nuclei. We identified 13 cell types that collectively represented 983,155 ATAC summits. We integrated genetic variation to discover 6,866 expression quantitative trait loci (eQTL) and 100,928 chromatin accessibility QTL (caQTL) (5% FDR) across the five most abundant cell types, cataloging caQTL peaks that atlas-level snATAC maps often miss. We identified 1,973 eGenes colocalized with caQTL and used mediation analyses to construct causal directional maps for chromatin accessibility and gene expression. 3,378 genome-wide association study (GWAS) signals across 43 relevant traits colocalized with sn-e/caQTL, 52% in a cell-specific manner. 77% of GWAS signals colocalized with caQTL and not eQTL, highlighting the critical importance of population-scale chromatin profiling for GWAS functional studies. GWAS-caQTL colocalization showed distinct cell-specific regulatory paradigms. For example, a C2CD4A/B T2D GWAS signal colocalized with caQTL in muscle fibers and multiple chromatin loop models nominated VPS13C, a glucose uptake gene. Sequence of the caQTL peak overlapping caSNP rs7163757 showed allelic regulatory activity differences in a human myocyte cell line massively parallel reporter assay. These results illuminate the genetic regulatory architecture of human skeletal muscle at high-resolution epigenomic, transcriptomic, and cell state scales and serve as a template for population-scale multi-omic mapping in complex tissues and traits.

A critical review on therapeutic approaches of CRISPR-Cas9 in diabetes mellitus

Diabetes mellitus (D.M.) is a common metabolic disorder caused mainly by combining two primary factors, which are (1) defects in insulin production by the pancreatic β-cells and (2) responsiveness of insulin-sensitive tissues towards insulin. Despite the rapid advancement in medicine to suppress elevated blood glucose levels (hyperglycemia) and insulin resistance associated with this hazard, a demand has undoubtedly emerged to find more effective and curative dimensions in therapeutic approaches against D.M. The administration of diabetes treatment that emphasizes insulin production and sensitivity may result in unfavorable side effects, reduced adherence, and potential treatment ineffectiveness. Recent progressions in genome editing technologies, for instance, in zinc-finger nucleases, transcription activator-like effector nucleases, and clustered regularly interspaced short palindromic repeat (CRISPR-Cas)-associated nucleases, have greatly influenced the gene editing technology from concepts to clinical practices. Improvements in genome editing technologies have also opened up the possibility to target and modify specific genome sequences in a cell directly. CRISPR/Cas9 has proven effective in utilizing ex vivo gene editing in embryonic stem cells and stem cells derived from patients. This application has facilitated the exploration of pancreatic beta-cell development and function. Furthermore, CRISPR/Cas9 enables the creation of innovative animal models for diabetes and assesses the effectiveness of different therapeutic strategies in treating the condition. We, therefore, present a critical review of the therapeutic approaches of the genome editing tool CRISPR-Cas9 in treating D.M., discussing the challenges and limitations of implementing this technology.

Transcriptional signature of islet neogenesis-associated protein peptide-treated rat pancreatic islets reveals induction of novel long non-coding RNAs

Background

Diabetes mellitus is characterized by chronic hyperglycemia with loss of β-cell function and mass. An attractive therapeutic approach to treat patients with diabetes in a non-invasive way is to harness the innate regenerative potential of the pancreas. The Islet Neogenesis-Associated Protein pentadecapeptide (INGAP-PP) has been shown to induce β-cell regeneration and improve their function in rodents. To investigate its possible mechanism of action, we report here the global transcriptional effects induced by the short-term INGAP-PP in vitro treatment of adult rat pancreatic islets.

Methods and findings

Rat pancreatic islets were cultured in vitro in the presence of INGAP-PP for 4 days, and RNA-seq was generated from triplicate treated and control islet samples. We performed a de novo rat gene annotation based on the alignment of RNA-seq reads. The list of INGAP-PP-regulated genes was integrated with epigenomic data. Using the new gene annotation generated in this work, we quantified RNA-seq data profiled in INS-1 cells treated with IL1β, IL1β+Calcipotriol (a vitamin D agonist) or vehicle, and single-cell RNA-seq data profiled in rat pancreatic islets. We found 1,669 differentially expressed genes by INGAP-PP treatment, including dozens of previously unannotated rat transcripts. Genes differentially expressed by the INGAP-PP treatment included a subset of upregulated transcripts that are associated with vitamin D receptor activation. Supported by epigenomic and single-cell RNA-seq data, we identified 9 previously unannotated long noncoding RNAs (lncRNAs) upregulated by INGAP-PP, some of which are also differentially regulated by IL1β and vitamin D in β-cells. These include Ri-lnc1, which is enriched in mature β-cells.

Conclusions

Our results reveal the transcriptional program that could explain the enhancement of INGAP-PP-mediated physiological effects on β-cell mass and function. We identified novel lncRNAs that are induced by INGAP-PP in rat islets, some of which are selectively expressed in pancreatic β-cells and downregulated by IL1β treatment of INS-1 cells. Our results suggest a relevant function for Ri-lnc1 in β-cells. These findings are expected to provide the basis for a deeper understanding of islet translational results from rodents to humans, with the ultimate goal of designing new therapies for people with diabetes.

Selective Transcription Factor Blockade Reduces Human Retinal Endothelial Cell Expression of Intercellular Adhesion Molecule-1 and Leukocyte Binding

The interaction between leukocytes and cytokine-activated retinal endothelium is an initiating step in non-infectious uveitis involving the posterior eye, mediated by cell adhesion molecules. However, because cell adhesion molecules are required for immune surveillance, therapeutic interventions would ideally be employed indirectly. Using 28 primary human retinal endothelial cell isolates, this study sought to identify transcription factor targets for reducing levels of the key retinal endothelial cell adhesion molecule, intercellular adhesion molecule (ICAM)-1, and limiting leukocyte binding to the retinal endothelium. Five candidate transcription factors-C2CD4B, EGR3, FOSB, IRF1, and JUNB-were identified by differential expression analysis of a transcriptome generated from IL-1β- or TNF-α-stimulated human retinal endothelial cells, interpreted in the context of the published literature. Further filtering involved molecular studies: of the five candidates, C2CD4B and IRF1 consistently demonstrated extended induction in IL-1β- or TNF-α-activated retinal endothelial cells and demonstrated a significant decrease in both ICAM-1 transcript and ICAM-1 membrane-bound protein expression by cytokine-activated retinal endothelial cells following treatment with small interfering RNA. RNA interference of C2CD4B or IRF1 significantly reduced leukocyte binding in a majority of human retinal endothelial cell isolates stimulated by IL-1β or TNF-α. Our observations suggest that the transcription factors C2CD4B and IRF1 may be potential drug targets for limiting leukocyte-retinal endothelial cell interactions in non-infectious uveitis involving the posterior eye.

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.

Multidimensional chromatin profiling of zebrafish pancreas to uncover and investigate disease-relevant enhancers

The pancreas is a central organ for human diseases. Most alleles uncovered by genome-wide association studies of pancreatic dysfunction traits overlap with non-coding sequences of DNA. Many contain epigenetic marks of cis-regulatory elements active in pancreatic cells, suggesting that alterations in these sequences contribute to pancreatic diseases. Animal models greatly help to understand the role of non-coding alterations in disease. However, interspecies identification of equivalent cis-regulatory elements faces fundamental challenges, including lack of sequence conservation. Here we combine epigenetic assays with reporter assays in zebrafish and human pancreatic cells to identify interspecies functionally equivalent cis-regulatory elements, regardless of sequence conservation. Among other potential disease-relevant enhancers, we identify a zebrafish ptf1a distal-enhancer whose deletion causes pancreatic agenesis, a phenotype previously found to be induced by mutations in a distal-enhancer of PTF1A in humans, further supporting the causality of this condition in vivo. This approach helps to uncover interspecies functionally equivalent cis-regulatory elements and their potential role in human disease.

Alterations in cis-regulatory elements (CREs) can contribute to pancreatic diseases. Here the authors combine chromatin profiling and interaction points with in vivo reporter assays in zebrafish to uncover functionally equivalent human CREs, helping to predict disease-relevant enhancers.

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.

The novel circSLC6A6/miR-1265/C2CD4A axis promotes colorectal cancer growth by suppressing p53 signaling pathway

Background

Colorectal cancer (CRC) is one of the most frequent malignancy and a leading cause of cancer-related deaths. Therefore, further researches are required to identify novel and more effective diagnoses and to identify molecular targets in treatment of CRC.

Methods

C2CD4A expression in CRC tissues and cell lines was detected by qRT-PCR and western blot. The biological functions of C2CD4A were performed both in vitro and in vivo. Western blot, cDNA array, IP-MS, Co-immunoprecipitation assay, and Ubiquitination assay were used to analyze the interaction between C2CD4A and p53. Bioinformatics analysis, FISH, RNA sequencing, luciferase reporter assay, RNA immunoprecipitation, RNA pull-down and rescue experiments, were deployed to detect upstream regulation mechanism of C2CD4A.

Results

C2CD4A was elevated in CRC tissues compared with adjacent normal colorectal tissues. C2CD4A knockdown significantly promoted cell apoptosis and with inhibited proliferation in vitro, and tumorigenicity in vivo, whereas C2CD4A overexpression led to opposite effects. Moreover, circSLC6A6 was upregulated and shown to positively regulate C2CD4A expression via sponging miR-1265. Fundamentally, C2CD4A inhibited p53 signaling pathway through interacting with p53 and increasing its ubiquitination and degradation.

Conclusion

Our results identified that circSLC6A6/miR-1265/C2CD4A axis, which was involved in CRC via the p53 signaling pathway, may serve as a therapeutic target for CRC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13046-021-02126-y.

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.

Mouse aging cell atlas analysis reveals global and cell type-specific aging signatures

Aging is associated with complex molecular and cellular processes that are poorly understood. Here we leveraged the Tabula Muris Senis single-cell RNA-seq data set to systematically characterize gene expression changes during aging across diverse cell types in the mouse. We identified aging-dependent genes in 76 tissue-cell types from 23 tissues and characterized both shared and tissue-cell-specific aging behaviors. We found that the aging-related genes shared by multiple tissue-cell types also change their expression congruently in the same direction during aging in most tissue-cell types, suggesting a coordinated global aging behavior at the organismal level. Scoring cells based on these shared aging genes allowed us to contrast the aging status of different tissues and cell types from a transcriptomic perspective. In addition, we identified genes that exhibit age-related expression changes specific to each functional category of tissue-cell types. Altogether, our analyses provide one of the most comprehensive and systematic characterizations of the molecular signatures of aging across diverse tissue-cell types in a mammalian system.

Sexually dimorphic roles for the type 2 diabetes-associated C2cd4b gene in murine glucose homeostasis

AIMS/HYPOTHESIS

Variants close to the VPS13C/C2CD4A/C2CD4B locus are associated with altered risk of type 2 diabetes in genome-wide association studies. While previous functional work has suggested roles for VPS13C and C2CD4A in disease development, none has explored the role of C2CD4B.

METHODS

CRISPR/Cas9-induced global C2cd4b-knockout mice and zebrafish larvae with c2cd4a deletion were used to study the role of this gene in glucose homeostasis. C2 calcium dependent domain containing protein (C2CD)4A and C2CD4B constructs tagged with FLAG or green fluorescent protein were generated to investigate subcellular dynamics using confocal or near-field microscopy and to identify interacting partners by mass spectrometry.

RESULTS

Systemic inactivation of C2cd4b in mice led to marked, but highly sexually dimorphic changes in body weight and glucose homeostasis. Female C2cd4b mice displayed unchanged body weight compared with control littermates, but abnormal glucose tolerance (AUC, p = 0.01) and defective in vivo, but not in vitro, insulin secretion (p = 0.02). This was associated with a marked decrease in follicle-stimulating hormone levels as compared with wild-type (WT) littermates (p = 0.003). In sharp contrast, male C2cd4b null mice displayed essentially normal glucose tolerance but an increase in body weight (p < 0.001) and fasting blood glucose (p = 0.003) after maintenance on a high-fat and -sucrose diet vs WT littermates. No metabolic disturbances were observed after global inactivation of C2cd4a in mice, or in pancreatic beta cell function at larval stages in C2cd4a null zebrafish. Fasting blood glucose levels were also unaltered in adult C2cd4a-null fish. C2CD4B and C2CD4A were partially localised to the plasma membrane, with the latter under the control of intracellular Ca2+. Binding partners for both included secretory-granule-localised PTPRN2/phogrin.

CONCLUSIONS/INTERPRETATION

Our studies suggest that C2cd4b may act centrally in the pituitary to influence sex-dependent circuits that control pancreatic beta cell function and glucose tolerance in rodents. However, the absence of sexual dimorphism in the impact of diabetes risk variants argues for additional roles for C2CD4A or VPS13C in the control of glucose homeostasis in humans.

DATA AVAILABILITY

The datasets generated and/or analysed during the current study are available in the Biorxiv repository ( www.biorxiv.org/content/10.1101/2020.05.18.099200v1 ). RNA-Seq (GSE152576) and proteomics (PXD021597) data have been deposited to GEO ( www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE152576 ) and ProteomeXchange ( www.ebi.ac.uk/pride/archive/projects/PXD021597 ) repositories, respectively.

Single-cell chromatin accessibility identifies pancreatic islet cell type- and state-specific regulatory programs of diabetes risk

Single-nucleus assay for transposase-accessible chromatin using sequencing (snATAC-seq) creates new opportunities to dissect cell type-specific mechanisms of complex diseases. Since pancreatic islets are central to type 2 diabetes (T2D), we profiled 15,298 islet cells by using combinatorial barcoding snATAC-seq and identified 12 clusters, including multiple alpha, beta and delta cell states. We cataloged 228,873 accessible chromatin sites and identified transcription factors underlying lineage- and state-specific regulation. We observed state-specific enrichment of fasting glucose and T2D genome-wide association studies for beta cells and enrichment for other endocrine cell types. At T2D signals localized to islet-accessible chromatin, we prioritized variants with predicted regulatory function and co-accessibility with target genes. A causal T2D variant rs231361 at the KCNQ1 locus had predicted effects on a beta cell enhancer co-accessible with INS and genome editing in embryonic stem cell-derived beta cells affected INS levels. Together our findings demonstrate the power of single-cell epigenomics for interpreting complex disease genetics.

Personizing the prediction of future susceptibility to a specific disease

A traceable biomarker is a member of a disease's molecular pathway. A disease may be associated with several molecular pathways. Each different combination of these molecular pathways, to which detected traceable biomarkers belong, may serve as an indicative of the elicitation of the disease at a different time frame in the future. Based on this notion, we introduce a novel methodology for personalizing an individual's degree of future susceptibility to a specific disease. We implemented the methodology in a working system called Susceptibility Degree to a Disease Predictor (SDDP). For a specific disease d, let S be the set of molecular pathways, to which traceable biomarkers detected from most patients of d belong. For the same disease d, let S' be the set of molecular pathways, to which traceable biomarkers detected from a certain individual belong. SDDP is able to infer the subset S'' ⊆{S-S'} of undetected molecular pathways for the individual. Thus, SDDP can infer undetected molecular pathways of a disease for an individual based on few molecular pathways detected from the individual. SDDP can also help in inferring the combination of molecular pathways in the set {S'+S''}, whose traceable biomarkers collectively is an indicative of the disease. SDDP is composed of the following four components: information extractor, interrelationship between molecular pathways modeler, logic inferencer, and risk indicator. The information extractor takes advantage of the exponential increase of biomedical literature to automatically extract the common traceable biomarkers for a specific disease. The interrelationship between molecular pathways modeler models the hierarchical interrelationships between the molecular pathways of the traceable biomarkers. The logic inferencer transforms the hierarchical interrelationships between the molecular pathways into rule-based specifications. It employs the specification rules and the inference rules for predicate logic to infer as many as possible undetected molecular pathways of a disease for an individual. The risk indicator outputs a risk indicator value that reflects the individual's degree of future susceptibility to the disease. We evaluated SDDP by comparing it experimentally with other methods. Results revealed marked improvement.

SNP-adjacent super enhancer network mediates enhanced osteogenic differentiation of MSCs in ankylosing spondylitis

Ankylosing spondylitis (AS) is a rheumatic disease with pathological osteogenesis that causes bony ankylosis and even deformity over time. Mesenchymal stem cells (MSCs) are multipotent stem cells that are the main source of osteoblasts. We previously demonstrated that enhanced osteogenic differentiation of MSCs from AS patients (ASMSCs) is related to pathological osteogenesis in AS. However, the more concrete mechanism needs further exploration. Super enhancers (SEs) are dense clusters of stitched enhancers that control cell identity determination and disease development. Single-nucleotide polymorphisms (SNPs) regulate the formation and interaction of SEs and denote genes accounting for AS susceptibility. Via integrative analysis of multiomic data, including histone 3 lysine 27 acetylation (H3K27ac), chromatin immunoprecipitation sequencing (ChIP-seq), SNPs and RNA sequencing (RNA-seq) data, we discovered a transcription network mediated by AS SNP-adjacent SEs (SASEs) in ASMSCs and identified key genes, such as Toll-like receptor 4 (TLR4), interleukin 18 receptor 1 (IL18R1), insulin-like growth factor binding protein 4 (IGFBP4), transportin 1 (TNPO1) and proprotein convertase subtilisin/kexin type 5 (PCSK5), which are pivotal in osteogenesis and AS pathogenesis. The SASE-regulated network modulates the enhanced osteogenic differentiation of ASMSCs by synergistically activating the PI3K-Akt, NF-kappaB and Hippo signaling pathways. Our results emphasize the crucial role of the SASE-regulated network in pathological osteogenesis in AS, and the preferential inhibition of ASMSC osteogenic differentiation by JQ1 indicates that SEs may be attractive targets in future treatment for new bone formation in AS.

Persistent or Transient Human β Cell Dysfunction Induced by Metabolic Stress: Specific Signatures and Shared Gene Expression with Type 2 Diabetes

Pancreatic β cell failure is key to type 2 diabetes (T2D) onset and progression. Here, we assess whether human β cell dysfunction induced by metabolic stress is reversible, evaluate the molecular pathways underlying persistent or transient damage, and explore the relationships with T2D islet traits. Twenty-six islet preparations are exposed to several lipotoxic/glucotoxic conditions, some of which impair insulin release, depending on stressor type, concentration, and combination. The reversal of dysfunction occurs after washout for some, although not all, of the lipoglucotoxic insults. Islet transcriptomes assessed by RNA sequencing and expression quantitative trait loci (eQTL) analysis identify specific pathways underlying β cell failure and recovery. Comparison of a large number of human T2D islet transcriptomes with those of persistent or reversible β cell lipoglucotoxicity show shared gene expression signatures. The identification of mechanisms associated with human β cell dysfunction and recovery and their overlap with T2D islet traits provide insights into T2D pathogenesis, fostering the development of improved β cell-targeted therapeutic strategies.

Functional Genomics in Pancreatic β Cells: Recent Advances in Gene Deletion and Genome Editing Technologies for Diabetes Research

The inheritance of variants that lead to coding changes in, or the mis-expression of, genes critical to pancreatic beta cell function can lead to alterations in insulin secretion and increase the risk of both type 1 and type 2 diabetes. Recently developed clustered regularly interspaced short palindromic repeats (CRISPR/Cas9) gene editing tools provide a powerful means of understanding the impact of identified variants on cell function, growth, and survival and might ultimately provide a means, most likely after the transplantation of genetically "corrected" cells, of treating the disease. Here, we review some of the disease-associated genes and variants whose roles have been probed up to now. Next, we survey recent exciting developments in CRISPR/Cas9 technology and their possible exploitation for β cell functional genomics. Finally, we will provide a perspective as to how CRISPR/Cas9 technology may find clinical application in patients with diabetes.

Cortical Transcriptomic Alterations in Association With Appetitive Neuropeptides and Body Mass Index in Posttraumatic Stress Disorder

BACKGROUND

The molecular pathology underlying posttraumatic stress disorder (PTSD) remains unclear mainly due to a lack of human PTSD postmortem brain tissue. The orexigenic neuropeptides ghrelin, neuropeptide Y, and hypocretin were recently implicated in modulating negative affect. Drawing from the largest functional genomics study of human PTSD postmortem tissue, we investigated whether there were molecular changes of these and other appetitive molecules. Further, we explored the interaction between PTSD and body mass index (BMI) on gene expression.

METHODS

We analyzed previously reported transcriptomic data from 4 prefrontal cortex regions from 52 individuals with PTSD and 46 matched neurotypical controls. We employed gene co-expression network analysis across the transcriptomes of these regions to uncover PTSD-specific networks containing orexigenic genes. We utilized Ingenuity Pathway Analysis software for pathway annotation. We identified differentially expressed genes (DEGs) among individuals with and without PTSD, stratified by sex and BMI.

RESULTS

Three PTSD-associated networks (P < .01) contained genes in signaling families of appetitive molecules: 2 in females and 1 in all subjects. We uncovered DEGs (P < .05) between PTSD and control subjects stratified by sex and BMI with especially robust changes in males with PTSD with elevated vs normal BMI. Further, we identified putative upstream regulators (P < .05) driving these changes, many of which were enriched for involvement in inflammation.

CONCLUSIONS

PTSD-associated cortical transcriptomic modules contain transcripts of appetitive genes, and BMI further interacts with PTSD to impact expression. DEGs and inferred upstream regulators of these modules could represent targets for future pharmacotherapies for obesity in PTSD.

A prognostic eight-gene expression signature for patients with breast cancer receiving adjuvant chemotherapy

Breast cancer is a popularly diagnosed malignant tumor. Genomic profiling studies suggest that breast cancer is a disease with heterogeneity. Chemotherapy is one of the chief means to treat breast cancer, while its responses and clinical outcomes vary largely due to the conventional clinicopathological factors and inherent chemosensitivity of breast cancer. Using the least absolute shrinkage and selection operator (LASSO) Cox regression model, our study established a multi-mRNA-based signature model and constructed a relative nomogram in predicting distant-recurrence-free survival for patients receiving surgery and following chemotherapy. We constructed a signature of eight mRNAs (IPCEF1, SYNDIG1, TIGIT, SPESP1, C2CD4A, CLCA2, RLN2, and CCL19) with the LASSO model, which was employed to separate subjects into groups with high- and low-risk scores. Obvious differences of distant-recurrence-free survival were found between these two groups. This eight-mRNA-based signature was independently associated with the prognosis and had better prognostic value than classical clinicopathologic factors according to multivariate Cox regression results. Receiver operating characteristic results demonstrated excellent performance in diagnosing 3-year distant-recurrence by the eight-mRNA signature. A nomogram that combined both the eight-mRNA-based signature and clinicopathological risk factors was constructed. Comparing with an ideal model, the nomograms worked well both in the training and validation sets. Through the results that the eight-mRNA signature effectively classified patients into low- and high-risk of distant recurrence, we concluded that this eight-mRNA-based signature played a promising predictive role in prognosis and could be clinically applied in breast cancer patients receiving adjuvant chemotherapy.

The Role of Noncoding Variants in Heritable Disease

The genetic basis of disease has largely focused on coding regions. However, it has become clear that a large proportion of the noncoding genome is functional and harbors genetic variants that contribute to disease etiology. Here, we review recent examples of inherited noncoding alterations that are responsible for Mendelian disorders or act to influence complex traits. We explore both rare and common genetic variants and discuss the wide range of mechanisms by which they affect gene regulation to promote disease. We also debate the challenges and progress associated with identifying and interpreting the functional and clinical significance of genetic variation in the context of the noncoding regulatory landscape.

A Long Non-coding RNA, LOC157273, Is an Effector Transcript at the Chromosome 8p23.1-PPP1R3B Metabolic Traits and Type 2 Diabetes Risk Locus

AIMS

Causal transcripts at genomic loci associated with type 2 diabetes (T2D) are mostly unknown. The chr8p23.1 variant rs4841132, associated with an insulin-resistant diabetes risk phenotype, lies in the second exon of a long non-coding RNA (lncRNA) gene, LOC157273, located 175 kilobases from PPP1R3B, which encodes a key protein regulating insulin-mediated hepatic glycogen storage in humans. We hypothesized that LOC157273 regulates expression of PPP1R3B in human hepatocytes.

METHODS

We tested our hypothesis using Stellaris fluorescent in situ hybridization to assess subcellular localization of LOC157273; small interfering RNA (siRNA) knockdown of LOC157273, followed by RT-PCR to quantify LOC157273 and PPP1R3B expression; RNA-seq to quantify the whole-transcriptome gene expression response to LOC157273 knockdown; and an insulin-stimulated assay to measure hepatocyte glycogen deposition before and after knockdown.

RESULTS

We found that siRNA knockdown decreased LOC157273 transcript levels by approximately 80%, increased PPP1R3B mRNA levels by 1.7-fold, and increased glycogen deposition by >50% in primary human hepatocytes. An A/G heterozygous carrier (vs. three G/G carriers) had reduced LOC157273 abundance due to reduced transcription of the A allele and increased PPP1R3B expression and glycogen deposition.

CONCLUSION

We show that the lncRNA LOC157273 is a negative regulator of PPP1R3B expression and glycogen deposition in human hepatocytes and a causal transcript at an insulin-resistant T2D risk locus.

(Epi)genomic heterogeneity of pancreatic islet function and failure in type 2 diabetes

Background

Pancreatic Islets of Langerhans are heterogeneous tissues consisting of multiple endocrine cell types that carry out distinct yet coordinated roles to regulate blood glucose homeostasis. Islet dysfunction and specifically failure of the beta cells to secrete adequate insulin are known precursors to type 2 diabetes (T2D) onset. However, the exact genetic, (epi)genomic, and environmental mechanisms that contribute to islet failure, and ultimately to T2D pathogenesis, require further elucidation.

Scope of review

This review summarizes efforts and advances in dissection of the complex genetic underpinnings of islet function and resilience in T2D pathogenesis. In this review, we will highlight results of the latest T2D genome-wide association study (GWAS) and discuss how these data are being combined with clinical measures in patients to uncover putative T2D subtypes and with functional (epi)genomic studies in islets to understand the genetic programming of islet cell identity, function, and adaptation. Finally, we discuss new and important opportunities to address major knowledge gaps in our understanding of islet (dys)function in T2D risk and progression.

Major conclusions

Genetic variation exerts clear effects on the islet epigenome, regulatory element usage, and gene expression. Future (epi)genomic comparative analyses between T2D and normal islets should incorporate genetics to distinguish patient-specific from disease-specific differences. Incorporating genotype information into future analyses and studies will also enable more precise insights into the molecular genetics of islet deficiency and failure in T2D risk, and should ultimately contribute to a stratified view of T2D and more precise treatment strategies. Islet cellular heterogeneity continues to remain a challenge for understanding the associations between islet failure and T2D development. Further efforts to obtain purified islet cell type populations and determine the specific genetic and environmental effects on each will help address this. Beyond observation of islets at steady state conditions, more research of islet stress and stimulation responses are needed to understand the transition of these tissues from a healthy to diseased state. Together, focusing on these objectives will provide more opportunities to prevent, treat, and manage T2D.

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.

Epigenetic Modifications Associated with the Pathogenesis of Type 2 Diabetes Mellitus

BACKGROUND AND OBJECTIVE

Type 2 diabetes mellitus (T2DM) is a multifactorial metabolic disorder. Pancreatic β-cell dysfunction and insulin resistance are the most common and crucial events of T2DM. Increasing evidence suggests the association of epigenetic modifications with the pathogenesis of T2DM through the changes in important biological processes including pancreatic β- cell differentiation, development and maintenance of normal β-cell function. Insulin sensitivity by the peripheral glucose uptake tissues is also changed by the altered epigenetic mechanisms. In this review, we discussed the major epigenetic alterations and their effects on β-cell function, insulin secretion and insulin resistance in context of T2DM.

METHODS

We investigated the presently available epigenetic modifications including DNA methylation, posttranslational histone modifications, ATP-dependent chromatin remodeling and non-coding RNAs related to the pathogenesis of T2DM. Published literatures on this topic were searched both on Google Scholar and Pubmed with related keywords and investigated for relevant information.

RESULTS

The epigenetic modifications introduce changes in gene expression which are essential for appropriate β-cell development and functions, insulin secretion and sensitivity resulting in the pathogenesis of T2DM. Interestingly, T2DM could also be a prominent reason for the mentioned epigenetic alterations.

CONCLUSION

This review article emphasized on the epigenetic modifications associated with T2DM and discussed the consequences in deterioration of the disease condition.

Single-cell ATAC-Seq in human pancreatic islets and deep learning upscaling of rare cells reveals cell-specific type 2 diabetes regulatory signatures

OBJECTIVE

Type 2 diabetes (T2D) is a complex disease characterized by pancreatic islet dysfunction, insulin resistance, and disruption of blood glucose levels. Genome-wide association studies (GWAS) have identified > 400 independent signals that encode genetic predisposition. More than 90% of associated single-nucleotide polymorphisms (SNPs) localize to non-coding regions and are enriched in chromatin-defined islet enhancer elements, indicating a strong transcriptional regulatory component to disease susceptibility. Pancreatic islets are a mixture of cell types that express distinct hormonal programs, so each cell type may contribute differentially to the underlying regulatory processes that modulate T2D-associated transcriptional circuits. Existing chromatin profiling methods such as ATAC-seq and DNase-seq, applied to islets in bulk, produce aggregate profiles that mask important cellular and regulatory heterogeneity.

METHODS

We present genome-wide single-cell chromatin accessibility profiles in >1,600 cells derived from a human pancreatic islet sample using single-cell combinatorial indexing ATAC-seq (sci-ATAC-seq). We also developed a deep learning model based on U-Net architecture to accurately predict open chromatin peak calls in rare cell populations.

RESULTS

We show that sci-ATAC-seq profiles allow us to deconvolve alpha, beta, and delta cell populations and identify cell-type-specific regulatory signatures underlying T2D. Particularly, T2D GWAS SNPs are significantly enriched in beta cell-specific and across cell-type shared islet open chromatin, but not in alpha or delta cell-specific open chromatin. We also demonstrate, using less abundant delta cells, that deep learning models can improve signal recovery and feature reconstruction of rarer cell populations. Finally, we use co-accessibility measures to nominate the cell-specific target genes at 104 non-coding T2D GWAS signals.

CONCLUSIONS

Collectively, we identify the islet cell type of action across genetic signals of T2D predisposition and provide higher-resolution mechanistic insights into genetically encoded risk pathways.

Pancreatic Islet Transcriptional Enhancers and Diabetes

PURPOSE OF REVIEW

Common genetic variants that associate with type 2 diabetes risk are markedly enriched in pancreatic islet transcriptional enhancers. This review discusses current advances in the annotation of islet enhancer variants and their target genes.

RECENT FINDINGS

Recent methodological advances now allow genetic and functional mapping of diabetes causal variants at unprecedented resolution. Mapping of enhancer-promoter interactions in human islets has provided a unique appreciation of the complexity of islet gene regulatory processes and enabled direct association of noncoding diabetes risk variants to their target genes. The recently improved human islet enhancer annotations constitute a framework for the interpretation of diabetes genetic signals in the context of pancreatic islet gene regulation. In the future, integration of existing and yet to come regulatory maps with genetic fine-mapping efforts and in-depth functional characterization will foster the discovery of novel diabetes molecular risk mechanisms.

Molecular mechanism for the multiple sclerosis risk variant rs17594362

Multiple sclerosis (MS) is known as an autoimmune demyelinating disease of the central nervous system. However, its cause remains elusive. Given previous studies suggesting that dysfunctional oligodendrocytes (OLs) may trigger MS, we tested whether single nucleotide polymorphisms (SNPs) associated with MS affect OL enhancers, potentially increasing MS risk by dysregulating gene expression of OL lineage cells. We found that two closely spaced OL enhancers, which are 3 Kb apart on chromosome 13, overlap two MS SNPs in linkage disequilibrium-rs17594362 and rs12429256. Our data revealed that the two MS SNPs significantly up-regulate the associated OL enhancers, which we have named as Rgcc-E1 and Rgcc-E2. Analysis of Hi-C data and epigenome editing experiments shows that Rgcc is the primary target of Rgcc-E1 and Rgcc-E2. Collectively, these data indicate that the molecular mechanism of rs17594362 and rs12429256 is to induce Rgcc overexpression by potentiating the enhancer activity of Rgcc-E1 and Rgcc-E2. Importantly, the dosage of the rs17594362/rs12429256 risk allele is positively correlated with the expression level of Rgcc in the human population, confirming our molecular mechanism. Our study also suggests that Rgcc overexpression in OL lineage cells may be a key cellular mechanism of rs17594362 and rs12429256 for MS.

Identification of C2CD4A as a human diabetes susceptibility gene with a role in β cell insulin secretion

Fine mapping and validation of genes causing β cell failure from susceptibility loci identified in type 2 diabetes genome-wide association studies (GWAS) poses a significant challenge. The VPS13C-C2CD4A-C2CD4B locus on chromosome 15 confers diabetes susceptibility in every ethnic group studied to date. However, the causative gene is unknown. FoxO1 is involved in the pathogenesis of β cell dysfunction, but its link to human diabetes GWAS has not been explored. Here we generated a genome-wide map of FoxO1 superenhancers in chemically identified β cells using 2-photon live-cell imaging to monitor FoxO1 localization. When parsed against human superenhancers and GWAS-derived diabetes susceptibility alleles, this map revealed a conserved superenhancer in C2CD4A, a gene encoding a β cell/stomach-enriched nuclear protein of unknown function. Genetic ablation of C2cd4a in β cells of mice phenocopied the metabolic abnormalities of human carriers of C2CD4A-linked polymorphisms, resulting in impaired insulin secretion during glucose tolerance tests as well as hyperglycemic clamps. C2CD4A regulates glycolytic genes, and notably represses key β cell "disallowed" genes, such as lactate dehydrogenase A We propose that C2CD4A is a transcriptional coregulator of the glycolytic pathway whose dysfunction accounts for the diabetes susceptibility associated with the chromosome 15 GWAS locus.

Integration of multiple epigenomic marks improves prediction of variant impact in saturation mutagenesis reporter assay

The integrative analysis of high-throughput reporter assays, machine learning, and profiles of epigenomic chromatin state in a broad array of cells and tissues has the potential to significantly improve our understanding of noncoding regulatory element function and its contribution to human disease. Here, we report results from the CAGI 5 regulation saturation challenge where participants were asked to predict the impact of nucleotide substitution at every base pair within five disease-associated human enhancers and nine disease-associated promoters. A library of mutations covering all bases was generated by saturation mutagenesis and altered activity was assessed in a massively parallel reporter assay (MPRA) in relevant cell lines. Reporter expression was measured relative to plasmid DNA to determine the impact of variants. The challenge was to predict the functional effects of variants on reporter expression. Comparative analysis of the full range of submitted prediction results identifies the most successful models of transcription factor binding sites, machine learning algorithms, and ways to choose among or incorporate diverse datatypes and cell-types for training computational models. These results have the potential to improve the design of future studies on more diverse sets of regulatory elements and aid the interpretation of disease-associated genetic variation.

CREB Promotes Beta Cell Gene Expression by Targeting Its Coactivators to Tissue-Specific Enhancers

CREB mediates effects of cyclic AMP on cellular gene expression. Ubiquitous CREB target genes are induced following recruitment of CREB and its coactivators to promoter proximal binding sites. We found that CREB stimulates the expression of pancreatic beta cell-specific genes by targeting CBP/p300 to promoter-distal enhancer regions. Subsequent increases in histone acetylation facilitate recruitment of the coactivators CRTC2 and BRD4, leading to release of RNA polymerase II over the target gene body. Indeed, CREB-induced hyperacetylation of chromatin over superenhancers promoted beta cell-restricted gene expression, which is sensitive to inhibitors of CBP/p300 and BRD4 activity. Neurod1 appears critical in establishing nucleosome-free regions for recruitment of CREB to beta cell-specific enhancers. Deletion of a CREB-Neurod1-bound enhancer within the Lrrc10b-Syt7 superenhancer disrupted the expression of both genes and decreased beta cell function. Our results demonstrate how cross talk between signal-dependent and lineage-determining factors promotes the expression of cell-type-specific gene programs in response to extracellular cues.

The Genetic Epidemiology of Type 2 Diabetes: Opportunities for Health Translation

PURPOSE OF REVIEW

Genome-wide association studies have delineated the genetic architecture of type 2 diabetes. While functional studies to identify target transcripts are ongoing, new genetic knowledge can be translated directly to health applications. The review covers several translation directions but focuses on genomic polygenic scores for screening and prevention.

RECENT FINDINGS

Over 400 genomic variants associated with T2D and its related quantitative traits are now known. Genetic scores comprising dozens to millions of associated variants can predict incident T2D. However, measurement of body mass index is more efficient than genetic scores to detect T2D risk groups, and knowledge of T2D genetic risk alone seems insufficient to improve health. Genetically determined metabolic sub-phenotypes can be identified by clustering variants associated with physiological axes like insulin resistance. Genetic sub-phenotyping may be a way forward to identify specific individual phenotypes for prevention and treatment. Genomic polygenic scores for T2D can predict incident diabetes but may not be useful to improve health overall. Genetic detection of T2D sub-phenotypes could be useful to personalize screening and care.

Enhancers: bridging the gap between gene control and human disease

Enhancers are a class of regulatory elements essential for precise spatio-temporal control of gene expression during development and in terminally differentiated cells. This review highlights signature features of enhancer elements as well as new advances that provide mechanistic insights into enhancer-mediated gene control in the context of three-dimensional chromatin. We detail the various ways in which non-coding mutations can instigate aberrant gene control and cause a variety of Mendelian disorders, common diseases and cancer.