Pancreatic islet chromatin accessibility and conformation reveals distal enhancer networks of type 2 diabetes risk

Palmitic acid alters enhancers/super-enhancers near inflammatory and efferocytosis-associated genes in human monocytes

Free fatty acids like palmitic acid (PA) are elevated in obesity and diabetes and dysregulate monocyte and macrophage functions, contributing to enhanced inflammation in these cardiometabolic diseases. Epigenetic mechanisms regulating enhancer functions play key roles in inflammatory gene expression, but their role in PA-induced monocyte/macrophage dysfunction is unknown. We found that PA treatment altered the epigenetic landscape of enhancers and super-enhancers (SEs) in human monocytes. Integration with RNA-seq data revealed that PA-induced enhancers/SEs correlated with PA-increased expression of inflammatory and immune response genes, while PA-inhibited enhancers correlated with downregulation of phagocytosis and efferocytosis genes. These genes were similarly regulated in macrophages from mouse models of diabetes and accelerated atherosclerosis, human atherosclerosis, and infectious agents. PA-regulated enhancers/SEs harbored SNPs associated with diabetes, obesity, and body mass index indicating disease relevance. We verified increased chromatin interactions between PA-regulated enhancers/SEs and inflammatory gene promoters and reduced interactions at efferocytosis genes. PA-induced gene expression was reduced by inhibitors of BRD4, and NF-κB. PA treatment inhibited phagocytosis and efferocytosis in human macrophages. Together, our findings demonstrate that PA-induced enhancer dynamics at key monocyte/macrophage enhancers/SEs regulate inflammatory and immune genes and responses. Targeting these PA-regulated epigenetic changes could provide novel therapeutic opportunities for cardiometabolic disorders.

DconnLoop: a deep learning model for predicting chromatin loops based on multi-source data integration

BACKGROUND

Chromatin loops are critical for the three-dimensional organization of the genome and gene regulation. Accurate identification of chromatin loops is essential for understanding the regulatory mechanisms in disease. However, current mainstream detection methods rely primarily on single-source data, such as Hi-C, which limits these methods' ability to capture the diverse features of chromatin loop structures. In contrast, multi-source data integration and deep learning approaches, though not yet widely applied, hold significant potential.

RESULTS

In this study, we developed a method called DconnLoop to integrate Hi-C, ChIP-seq, and ATAC-seq data to predict chromatin loops. This method achieves feature extraction and fusion of multi-source data by integrating residual mechanisms, directional connectivity excitation modules, and interactive feature space decoders. Finally, we apply density estimation and density clustering to the genome-wide prediction results to identify more representative loops. The code is available from https://github.com/kuikui-C/DconnLoop .

CONCLUSIONS

The results demonstrate that DconnLoop outperforms existing methods in both precision and recall. In various experiments, including Aggregate Peak Analysis and peak enrichment comparisons, DconnLoop consistently shows advantages. Extensive ablation studies and validation across different sequencing depths further confirm DconnLoop's robustness and generalizability.

Exploring and mitigating shortcomings in single-cell differential expression analysis with a new statistical paradigm

BACKGROUND

Differential expression analysis is pivotal in single-cell transcriptomics for unraveling cell-type-specific responses to stimuli. While numerous methods are available to identify differentially expressed genes in single-cell data, recent evaluations of both single-cell-specific methods and methods adapted from bulk studies have revealed significant shortcomings in performance. In this paper, we dissect the four major challenges in single-cell differential expression analysis: excessive zeros, normalization, donor effects, and cumulative biases. These "curses" underscore the limitations and conceptual pitfalls in existing workflows.

RESULTS

To address the limitations of current single-cell differential expression analysis methods, we propose GLIMES, a statistical framework that leverages UMI counts and zero proportions within a generalized Poisson/Binomial mixed-effects model to account for batch effects and within-sample variation. We rigorously benchmarked GLIMES against six existing differential expression methods using three case studies and simulations across different experimental scenarios, including comparisons across cell types, tissue regions, and cell states. Our results demonstrate that GLIMES is more adaptable to diverse experimental designs in single-cell studies and effectively mitigates key shortcomings of current approaches, particularly those related to normalization procedures. By preserving biologically meaningful signals, GLIMES offers improved performance in detecting differentially expressed genes.

CONCLUSIONS

By using absolute RNA expression rather than relative abundance, GLIMES improves sensitivity, reduces false discoveries, and enhances biological interpretability. This paradigm shift challenges existing workflows and highlights the need for careful consideration of normalization strategies, ultimately paving the way for more accurate and robust single-cell transcriptomic analyses.

Construction of cynomolgus monkey type 2 diabetes models by combining genetic prediction model with high-energy diet

BACKGROUND

Type 2 diabetes mellitus (T2D) is a significant health concern. Research using non-human primates, which develop T2D with similar symptoms and pancreatic changes as humans, is crucial but limited by long timelines and low success rates.

RESULTS

We targeted capture sequenced 61 normal and 81 T2D cynomolgus monkeys using a primer panel that captured 269 potential regulatory regions potentially associated with T2D in the cynomolgus monkey genome. 80 variants were identified to be associated with T2D and were used to construct a genetic prediction model. Among 8 machine learning algorithms tested, we found that the best prediction performance was achieve when the model using support vector machine with polynomial kernel as the machine learning algorithm (AUC = 0.933). Including age and sex in this model did not significantly improve the prediction performance. Using the genetic prediction model, we further screened 22 monkeys and found 13 were high risk while 9 were low risk. After feeding the 22 monkeys with high-energy food for 32 weeks, we found all the 9 low risk monkeys did not develop T2D while 4 out of 13 high risk monkeys (31 %) develop T2D.

CONCLUSIONS

This method greatly increased the success rate of establishing T2D monkey models while decreased the time needed compared to traditional methods. Therefore, we developed a new high-efficiency method to establish T2D monkey models by combining the genetic prediction model and high-energy diet, which will greatly contribute to the research on the clinical characteristics, pathogenesis, complications and potential new treatments.

Genetic coupling of enhancer activity and connectivity in gene expression control

Gene enhancers often form long-range contacts with promoters, but it remains unclear if the activity of enhancers and their chromosomal contacts are mediated by the same DNA sequences and recruited factors. Here, we study the effects of expression quantitative trait loci (eQTLs) on enhancer activity and promoter contacts in primary monocytes isolated from 34 male individuals. Using eQTL-Capture Hi-C and a Bayesian approach considering both intra- and inter-individual variation, we initially detect 19 eQTLs associated with enhancer-eGene promoter contacts, most of which also associate with enhancer accessibility and activity. Capitalising on these shared effects, we devise a multi-modality Bayesian strategy, identifying 629 "trimodal QTLs" jointly associated with enhancer accessibility, eGene promoter contact, and gene expression. Causal mediation analysis and CRISPR interference reveal causal relationships between these three modalities. Many detected QTLs overlap disease susceptibility loci and influence the predicted binding of myeloid transcription factors, including SPI1, GABPB and STAT3. Additionally, a variant associated with PCK2 promoter contact directly disrupts a CTCF binding motif and impacts promoter insulation from downstream enhancers. Jointly, our findings suggest an inherent genetic coupling of enhancer activity and connectivity in gene expression control relevant to human disease and highlight the regulatory role of genetically determined chromatin boundaries.

Cell-type-resolved chromatin accessibility in the human intestine identifies complex regulatory programs and clarifies genetic associations in Crohn's disease

Crohn's disease (CD) is a complex inflammatory bowel disease resulting from an interplay of genetic, microbial, and environmental factors. Cell-type-specific contributions to CD etiology and genetic risk are incompletely understood. Here we built a comprehensive atlas of cell-type- resolved chromatin accessibility comprising 557,310 candidate cis-regulatory elements (cCREs) in terminal ileum and ascending colon from patients with active and inactive CD and healthy controls. Using this atlas, we identified cell-type-, anatomic location-, and context-specific cCREs and characterized the regulatory programs underlying inflammatory responses in the intestinal mucosa of CD patients. Genetic variants that disrupt binding motifs of cell-type-specific transcription factors significantly affected chromatin accessibility in specific mucosal cell types. We found that CD heritability is primarily enriched in immune cell types. However, using fine- mapped non-coding CD variants we identified 29 variants located within cCREs several of which were accessible in epithelial and stromal cells implicating cell types from additional lineages in mediating CD risk in some loci. Our atlas provides a comprehensive resource to study gene regulatory effects in CD and health, and highlights the cellular complexity underlying CD risk.

Circ_0000284 Is Involved in Arsenite-Induced Hepatic Insulin Resistance Through Blocking the Plasma Membrane Translocation of GLUT4 in Hepatocytes via IGF2BP2/PPAR-γ

Arsenic exposure can induce liver insulin resistance (IR) and diabetes (DM), but the underlying mechanisms are not yet clear. Circular RNAs (circRNAs) are involved in the regulation of the onset of diabetes, especially in the progression of IR. This study aimed to investigate the role of circRNAs in arsenic-induced hepatic IR and its underlying mechanism. Male C57BL/6J mice were given drinking water containing sodium arsenite (0, 0.5, 5, or 50 ppm) for 12 months. The results show that sodium arsenite increased circ_0000284 expression, decreased insulin-like growth factor 2 mRNA binding protein 2 (IGF2BP2) and peroxisome proliferator-activated receptor-γ (PPAR-γ), and inhibited cell membrane protein levels of insulin-responsive glucose transporter protein 4 (GLUT4) in the mouse livers, indicating that arsenic exposure causes liver damage and disruptions to glucose metabolism. Furthermore, sodium arsenite reduced glucose consumption and glycogen levels, increased the expression of circ_0000284, reduced the protein levels of IGF2BP2 and PPAR-γ, and inhibited GLUT4 protein levels in the cell membranes of insulin-treated HepG2 cells. However, a circ_0000284 inhibitor reversed arsenic exposure-induced reductions in IGF2BP2, PPAR-γ, and GLUT4 levels in the plasma membrane. These results indicate that circ_0000284 is involved in arsenite-induced hepatic insulin resistance through blocking the plasma membrane translocation of GLUT4 in hepatocytes via IGF2BP2/PPAR-γ. This study provides a scientific basis for finding early biomarkers for the control of arsenic exposure and type 2 diabetes mellitus (T2DM), and discovering new prevention and control measures.

Chromatin accessibility: biological functions, molecular mechanisms and therapeutic application

The dynamic regulation of chromatin accessibility is one of the prominent characteristics of eukaryotic genome. The inaccessible regions are mainly located in heterochromatin, which is multilevel compressed and access restricted. The remaining accessible loci are generally located in the euchromatin, which have less nucleosome occupancy and higher regulatory activity. The opening of chromatin is the most important prerequisite for DNA transcription, replication, and damage repair, which is regulated by genetic, epigenetic, environmental, and other factors, playing a vital role in multiple biological progresses. Currently, based on the susceptibility difference of occupied or free DNA to enzymatic cleavage, solubility, methylation, and transposition, there are many methods to detect chromatin accessibility both in bulk and single-cell level. Through combining with high-throughput sequencing, the genome-wide chromatin accessibility landscape of many tissues and cells types also have been constructed. The chromatin accessibility feature is distinct in different tissues and biological states. Research on the regulation network of chromatin accessibility is crucial for uncovering the secret of various biological processes. In this review, we comprehensively introduced the major functions and mechanisms of chromatin accessibility variation in different physiological and pathological processes, meanwhile, the targeted therapies based on chromatin dynamics regulation are also summarized.

Beta cells are essential drivers of pancreatic ductal adenocarcinoma development

Pancreatic endocrine-exocrine crosstalk plays a key role in normal physiology and disease. For instance, endocrine islet beta (β) cell secretion of insulin or cholecystokinin (CCK) promotes progression of pancreatic adenocarcinoma (PDAC), an exocrine cell-derived tumor. However, the cellular and molecular mechanisms that govern endocrine-exocrine signaling in tumorigenesis remain incompletely understood. We find that β cell ablation impedes PDAC development in mice, arguing that the endocrine pancreas is critical for exocrine tumorigenesis. Conversely, obesity induces β cell hormone dysregulation, alters CCK-dependent peri-islet exocrine cell transcriptional states, and enhances islet proximal tumor formation. Single-cell RNA-sequencing, in silico latent-space archetypal and trajectory analysis, and genetic lineage tracing in vivo reveal that obesity stimulates postnatal immature β cell expansion and adaptation towards a pro-tumorigenic CCK+ state via JNK/cJun stress-responsive signaling. These results define endocrine-exocrine signaling as a driver of PDAC development and uncover new avenues to target the endocrine pancreas to subvert exocrine tumorigenesis.

A ONECUT1 regulatory, non-coding region in pancreatic development and diabetes

In a patient with permanent neonatal syndromic diabetes clinically similar to cases with ONECUT1 biallelic mutations, we identified a disease-causing deletion located upstream of ONECUT1. Through genetic, genomic, and functional studies, we identified a crucial regulatory region acting as an enhancer of ONECUT1 specifically during pancreatic development. This enhancer region contains a low-frequency variant showing a strong association with type 2 diabetes and other glycemic traits, thus extending the contribution of this region to common forms of diabetes. Clinical relevance is provided by experimentally tailored therapy options for patients carrying ONECUT1 coding or regulatory mutations.

Predictive Prioritization of Enhancers Associated with Pancreas Disease Risk

Genetic and epigenetic variations in regulatory enhancer elements increase susceptibility to a range of pathologies. Despite recent advances, linking enhancer elements to target genes and predicting transcriptional outcomes of enhancer dysfunction remain significant challenges. Using 3D chromatin conformation assays, we generated an extensive enhancer interaction dataset for the human pancreas, encompassing more than 20 donors and five major cell types, including both exocrine and endocrine compartments. We employed a network approach to parse chromatin interactions into enhancer-promoter tree models, facilitating a quantitative, genome-wide analysis of enhancer connectivity. With these tree models, we developed a machine learning algorithm to estimate the impact of enhancer perturbations on cell type-specific gene expression in the human pancreas. Orthogonal to our computational approach, we perturbed enhancer function in primary human pancreas cells using CRISPR interference and quantified the effects at the single-cell level through RNA FISH coupled with high-throughput imaging. Our enhancer tree models enabled the annotation of common germline risk variants associated with pancreas diseases, linking them to putative target genes in specific cell types. For pancreatic ductal adenocarcinoma, we found a stronger enrichment of disease susceptibility variants within acinar cell regulatory elements, despite ductal cells historically being assumed as the primary cell-of-origin. Our integrative approach-combining cell type-specific enhancer-promoter interaction mapping, computational models, and single-cell enhancer perturbation assays-produced a robust resource for studying the genetic basis of pancreas disorders.

Identifying patterns differing between high-dimensional datasets with generalized contrastive PCA

High-dimensional data have become ubiquitous in the biological sciences, and it is often desirable to compare two datasets collected under different experimental conditions to extract low-dimensional patterns enriched in one condition. However, traditional dimensionality reduction techniques cannot accomplish this because they operate on only one dataset. Contrastive principal component analysis (cPCA) has been proposed to address this problem, but it has seen little adoption because it requires tuning a hyperparameter resulting in multiple solutions, with no way of knowing which is correct. Moreover, cPCA uses foreground and background conditions that are treated differently, making it ill-suited to compare two experimental conditions symmetrically. Here we describe the development of generalized contrastive PCA (gcPCA), a flexible hyperparameter-free approach that solves these problems. We first provide analyses explaining why cPCA requires a hyperparameter and how gcPCA avoids this requirement. We then describe an open-source gcPCA toolbox containing Python and MATLAB implementations of several variants of gcPCA tailored for different scenarios. Finally, we demonstrate the utility of gcPCA in analyzing diverse high-dimensional biological data, revealing unsupervised detection of hippocampal replay in neurophysiological recordings and heterogeneity of type II diabetes in single-cell RNA sequencing data. As a fast, robust, and easy-to-use comparison method, gcPCA provides a valuable resource facilitating the analysis of diverse high-dimensional datasets to gain new insights into complex biological phenomena.

Single cell multiome profiling of pancreatic islets reveals physiological changes in cell type-specific regulation associated with diabetes risk

Physiological variability in pancreatic cell type gene regulation and the impact on diabetes risk is poorly understood. In this study we mapped gene regulation in pancreatic cell types using single cell multiomic (joint RNA-seq and ATAC-seq) profiling in 28 non-diabetic donors in combination with single cell data from 35 non-diabetic donors in the Human Pancreas Analysis Program. We identified widespread associations with age, sex, BMI, and HbA1c, where gene regulatory responses were highly cell type- and phenotype-specific. In beta cells, donor age associated with hypoxia, apoptosis, unfolded protein response, and external signal-dependent transcriptional regulators, while HbA1c associated with inflammatory responses and gender with chromatin organization. We identified 10.8K loci where genetic variants were QTLs for cis regulatory element (cRE) accessibility, including 20% with lineage- or cell type-specific effects which disrupted distinct transcription factor motifs. Type 2 diabetes and glycemic trait associated variants were enriched in both phenotype- and QTL-associated beta cell cREs, whereas type 1 diabetes showed limited enrichment. Variants at 226 diabetes and glycemic trait loci were QTLs in beta and other cell types, including 40 that were statistically colocalized, and annotating target genes of colocalized QTLs revealed genes with putatively novel roles in disease. Our findings reveal diverse responses of pancreatic cell types to phenotype and genotype in physiology, and identify pathways, networks, and genes through which physiology impacts diabetes risk.

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

BACKGROUND

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

SCOPE OF REVIEW

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

MAJOR CONCLUSIONS

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

Chromosome 20p11.2 deletions cause congenital hyperinsulinism via the loss of FOXA2 or its regulatory elements

Persistent congenital hyperinsulinism (HI) is a rare genetically heterogeneous condition characterised by dysregulated insulin secretion leading to life-threatening hypoglycaemia. For up to 50% of affected individuals screening of the known HI genes does not identify a disease-causing variant. Large deletions have previously been used to identify novel regulatory regions causing HI. Here, we used genome sequencing to search for novel large (>1 Mb) deletions in 180 probands with HI of unknown cause and replicated our findings in a large cohort of 883 genetically unsolved individuals with HI using off-target copy number variant calling from targeted gene panels. We identified overlapping heterozygous deletions in five individuals (range 3-8 Mb) spanning chromosome 20p11.2. The pancreatic beta-cell transcription factor gene, FOXA2, a known cause of HI was deleted in two of the five individuals. In the remaining three, we found a minimal deleted region of 2.4 Mb adjacent to FOXA2 that encompasses multiple non-coding regulatory elements that are in conformational contact with FOXA2. Our data suggests that the deletions in these three children may cause disease through the dysregulation of FOXA2 expression. These findings provide new insights into the regulation of FOXA2 in the beta-cell and confirm an aetiological role for chromosome 20p11.2 deletions in syndromic HI.

Senescence of human pancreatic beta cells enhances functional maturation through chromatin reorganization and promotes interferon responsiveness

Senescent cells can influence the function of tissues in which they reside, and their propensity for disease. A portion of adult human pancreatic beta cells express the senescence marker p16, yet it is unclear whether they are in a senescent state, and how this affects insulin secretion. We analyzed single-cell transcriptome datasets of adult human beta cells, and found that p16-positive cells express senescence gene signatures, as well as elevated levels of beta-cell maturation genes, consistent with enhanced functionality. Senescent human beta-like cells in culture undergo chromatin reorganization that leads to activation of enhancers regulating functional maturation genes and acquisition of glucose-stimulated insulin secretion capacity. Strikingly, Interferon-stimulated genes are elevated in senescent human beta cells, but genes encoding senescence-associated secretory phenotype (SASP) cytokines are not. Senescent beta cells in culture and in human tissue show elevated levels of cytoplasmic DNA, contributing to their increased interferon responsiveness. Human beta-cell senescence thus involves chromatin-driven upregulation of a functional-maturation program, and increased responsiveness of interferon-stimulated genes, changes that could increase both insulin secretion and immune reactivity.

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

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

The curses of performing differential expression analysis using single-cell data

Differential expression analysis is pivotal in single-cell transcriptomics for unraveling cell-type- specific responses to stimuli. While numerous methods are available to identify differentially expressed genes in single-cell data, recent evaluations of both single-cell-specific methods and methods adapted from bulk studies have revealed significant shortcomings in performance. In this paper, we dissect the four major challenges in single-cell DE analysis: normalization, excessive zeros, donor effects, and cumulative biases. These "curses" underscore the limitations and conceptual pitfalls in existing workflows. In response, we introduce a novel paradigm addressing several of these issues.

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.

An Enhancer Within Abcb11 Regulates G6pc2 in C57BL/6 Mouse Pancreatic Islets

G6PC2 is predominantly expressed in pancreatic islet β-cells where it encodes a glucose-6-phosphatase catalytic subunit that modulates the sensitivity of insulin secretion to glucose by opposing the action of glucokinase, thereby regulating fasting blood glucose (FBG). Prior studies have shown that the G6pc2 promoter alone is unable to confer sustained islet-specific gene expression in mice, suggesting the existence of distal enhancers that regulate G6pc2 expression. Using information from both mice and humans and knowledge that single nucleotide polymorphisms (SNPs) both within and near G6PC2 are associated with variations in FBG in humans, we identified several putative enhancers 3' of G6pc2. One region, herein referred to as enhancer I, resides in the 25th intron of Abcb11 and binds multiple islet-enriched transcription factors. CRISPR-mediated deletion of enhancer I in C57BL/6 mice had selective effects on the expression of genes near the G6pc2 locus. In isolated islets, G6pc2 and Spc25 expression were reduced ∼50%, and Gm13613 expression was abolished, whereas Cers6 and nostrin expression were unaffected. This partial reduction in G6pc2 expression enhanced islet insulin secretion at basal glucose concentrations but did not affect FBG or glucose tolerance in vivo, consistent with the absence of a phenotype in G6pc2 heterozygous C57BL/6 mice.

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, combined allele-specific expression (cASE) analysis in human islets revealed multiple variants that influence SLC30A8 expression. Epigenomic mapping identified an islet-selective enhancer cluster at the SLC30A8 locus, hosting multiple T2D risk and cASE associations, which is spatially associated with the SLC30A8 promoter and additional neighbouring genes. Deletions of variant-bearing enhancer regions using CRISPR-Cas9 in human-derived EndoC-βH3 cells lowered the expression of SLC30A8 and several neighbouring genes, and improved insulin secretion. Whilst down-regulation of SLC30A8 had no effect on beta cell survival, loss of UTP23, RAD21 or MED30 markedly reduced cell viability. Although eQTL or cASE analyses in human islets did not support the association between these additional genes and diabetes risk, the transcriptional regulator JQ1 lowered the expression of multiple genes at the SLC30A8 locus and enhanced stimulated insulin secretion.

Regulation of CTCF loop formation during pancreatic cell differentiation

Transcription reprogramming during cell differentiation involves targeting enhancers to genes responsible for establishment of cell fates. To understand the contribution of CTCF-mediated chromatin organization to cell lineage commitment, we analyzed 3D chromatin architecture during the differentiation of human embryonic stem cells into pancreatic islet organoids. We find that CTCF loops are formed and disassembled at different stages of the differentiation process by either recruitment of CTCF to new anchor sites or use of pre-existing sites not previously involved in loop formation. Recruitment of CTCF to new sites in the genome involves demethylation of H3K9me3 to H3K9me2, demethylation of DNA, recruitment of pioneer factors, and positioning of nucleosomes flanking the new CTCF sites. Existing CTCF sites not involved in loop formation become functional loop anchors via the establishment of new cohesin loading sites containing NIPBL and YY1 at sites between the new anchors. In both cases, formation of new CTCF loops leads to strengthening of enhancer promoter interactions and increased transcription of genes adjacent to loop anchors. These results suggest an important role for CTCF and cohesin in controlling gene expression during cell differentiation.

CHARGE syndrome-associated CHD7 acts at ISL1-regulated enhancers to modulate second heart field gene expression

AIMS

Haploinsufficiency of the chromo-domain protein CHD7 underlies most cases of CHARGE syndrome, a multisystem birth defect including congenital heart malformation. Context specific roles for CHD7 in various stem, progenitor, and differentiated cell lineages have been reported. Previously, we showed severe defects when Chd7 is absent from cardiopharyngeal mesoderm (CPM). Here, we investigate altered gene expression in the CPM and identify specific CHD7-bound target genes with known roles in the morphogenesis of affected structures.

METHODS AND RESULTS

We generated conditional KO of Chd7 in CPM and analysed cardiac progenitor cells using transcriptomic and epigenomic analyses, in vivo expression analysis, and bioinformatic comparisons with existing datasets. We show CHD7 is required for correct expression of several genes established as major players in cardiac development, especially within the second heart field (SHF). We identified CHD7 binding sites in cardiac progenitor cells and found strong association with histone marks suggestive of dynamically regulated enhancers during the mesodermal to cardiac progenitor transition of mESC differentiation. Moreover, CHD7 shares a subset of its target sites with ISL1, a pioneer transcription factor in the cardiogenic gene regulatory network, including one enhancer modulating Fgf10 expression in SHF progenitor cells vs. differentiating cardiomyocytes.

CONCLUSION

We show that CHD7 interacts with ISL1, binds ISL1-regulated cardiac enhancers, and modulates gene expression across the mesodermal heart fields during cardiac morphogenesis.

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

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

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

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

Prioritization of genes associated with type 2 diabetes mellitus for functional studies

Existing therapies for type 2 diabetes mellitus (T2DM) show limited efficacy or have adverse effects. Numerous genetic variants associated with T2DM have been identified, but progress in translating these findings into potential drug targets has been limited. Here, we describe the tools and platforms available to identify effector genes from T2DM-associated coding and non-coding variants and prioritize them for functional studies. We discuss QSER1 and SLC12A8 as examples of genes that have been identified as possible T2DM candidate genes using these tools and platforms. We suggest further approaches, including the use of sequencing data with increased sample size and ethnic diversity, single-cell omics data for analyses, glycaemic trait associations to predict gene function and, potentially, human induced pluripotent stem cell 'village' cultures, to strengthen current gene functionalization workflows. Effective prioritization of T2DM-associated genes for experimental validation could expedite our understanding of the genetic mechanisms responsible for T2DM to facilitate the use of precision medicine in its treatment.

Chromatin accessibility differences between alpha, beta, and delta cells identifies common and cell type-specific enhancers

BACKGROUND

High throughput sequencing has enabled the interrogation of the transcriptomic landscape of glucagon-secreting alpha cells, insulin-secreting beta cells, and somatostatin-secreting delta cells. These approaches have furthered our understanding of expression patterns that define healthy or diseased islet cell types and helped explicate some of the intricacies between major islet cell crosstalk and glucose regulation. All three endocrine cell types derive from a common pancreatic progenitor, yet alpha and beta cells have partially opposing functions, and delta cells modulate and control insulin and glucagon release. While gene expression signatures that define and maintain cellular identity have been widely explored, the underlying epigenetic components are incompletely characterized and understood. However, chromatin accessibility and remodeling is a dynamic attribute that plays a critical role to determine and maintain cellular identity.

RESULTS

Here, we compare and contrast the chromatin landscape between mouse alpha, beta, and delta cells using ATAC-Seq to evaluate the significant differences in chromatin accessibility. The similarities and differences in chromatin accessibility between these related islet endocrine cells help define their fate in support of their distinct functional roles. We identify patterns that suggest that both alpha and delta cells are poised, but repressed, from becoming beta-like. We also identify patterns in differentially enriched chromatin that have transcription factor motifs preferentially associated with different regions of the genome. Finally, we not only confirm and visualize previously discovered common endocrine- and cell specific- enhancer regions across differentially enriched chromatin, but identify novel regions as well. We compiled our chromatin accessibility data in a freely accessible database of common endocrine- and cell specific-enhancer regions that can be navigated with minimal bioinformatics expertise.

CONCLUSIONS

Both alpha and delta cells appear poised, but repressed, from becoming beta cells in murine pancreatic islets. These data broadly support earlier findings on the plasticity in identity of non-beta cells under certain circumstances. Furthermore, differential chromatin accessibility shows preferentially enriched distal-intergenic regions in beta cells, when compared to either alpha or delta cells.

Dynamic chromatin architectures provide insights into the genetics of cattle myogenesis

BACKGROUND

Sharply increased beef consumption is propelling the genetic improvement projects of beef cattle in China. Three-dimensional genome structure is confirmed to be an important layer of transcription regulation. Although genome-wide interaction data of several livestock species have already been produced, the genome structure states and its regulatory rules in cattle muscle are still limited.

RESULTS

Here we present the first 3D genome data in Longissimus dorsi muscle of fetal and adult cattle (Bos taurus). We showed that compartments, topologically associating domains (TADs), and loop undergo re-organization and the structure dynamics were consistent with transcriptomic divergence during muscle development. Furthermore, we annotated cis-regulatory elements in cattle genome during myogenesis and demonstrated the enrichments of promoter and enhancer in selection sweeps. We further validated the regulatory function of one HMGA2 intronic enhancer near a strong sweep region on primary bovine myoblast proliferation.

CONCLUSIONS

Our data provide key insights of the regulatory function of high order chromatin structure and cattle myogenic biology, which will benefit the progress of genetic improvement of beef cattle.

AdaLiftOver: high-resolution identification of orthologous regulatory elements with Adaptive liftOver

MOTIVATION

Elucidating functionally similar orthologous regulatory regions for human and model organism genomes is critical for exploiting model organism research and advancing our understanding of results from genome-wide association studies (GWAS). Sequence conservation is the de facto approach for finding orthologous non-coding regions between human and model organism genomes. However, existing methods for mapping non-coding genomic regions across species are challenged by the multi-mapping, low precision, and low mapping rate issues.

RESULTS

We develop Adaptive liftOver (AdaLiftOver), a large-scale computational tool for identifying functionally similar orthologous non-coding regions across species. AdaLiftOver builds on the UCSC liftOver framework to extend the query regions and prioritizes the resulting candidate target regions based on the conservation of the epigenomic and the sequence grammar features. Evaluations of AdaLiftOver with multiple case studies, spanning both genomic intervals from epigenome datasets across a wide range of model organisms and GWAS SNPs, yield AdaLiftOver as a versatile method for deriving hard-to-obtain human epigenome datasets as well as reliably identifying orthologous loci for GWAS SNPs.

AVAILABILITY AND IMPLEMENTATION

The R package and the data for AdaLiftOver is available from https://github.com/keleslab/AdaLiftOver.

Insights from multi-omics integration in complex disease primary tissues

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

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

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

Human T2D-Associated Gene IMP2/IGF2BP2 Promotes the Commitment of Mesenchymal Stem Cells Into Adipogenic Lineage

Excessive adiposity is the main cause of obesity and type two diabetes (T2D). Variants in human IMP2/IGF2BP2 gene are associated with increased risk of T2D. However, little is known about its role in adipogenesis and in insulin resistance. Here, we investigate the function of IMP2 during adipocyte development. Mice with Imp2 deletion in mesenchymal stem cells (MSC) are resistant to diet-induced obesity without glucose and insulin tolerance affected. Imp2 is essential for the early commitment of adipocyte-derived stem cells (ADSC) into preadipocytes, but the deletion of Imp2 in MSC is not required for the proliferation and terminal differentiation of committed preadipocytes. Mechanistically, Imp2 binds Wnt receptor Fzd8 mRNA and promotes its degradation by recruiting CCR4-NOT deadenylase complex in an mTOR-dependent manner. Our data demonstrate that Imp2 is required for maintaining white adipose tissue homeostasis through controlling mRNA stability in ADSC. However, the contribution of IMP2 to insulin resistance, a main risk of T2D, is not evident.

Lineage-specific rearrangement of chromatin loops and epigenomic features during adipocytes and osteoblasts commitment

Human mesenchymal stem cells (hMSCs) can be differentiated into adipocytes and osteoblasts. The processes are driven by the rewiring of chromatin architectures and transcriptomic/epigenomic changes. Here, we induced hMSCs to adipogenic and osteogenic differentiation, and performed 2 kb resolution Hi-C experiments for chromatin loops detection. We also generated matched RNA-seq, ChIP-seq and ATAC-seq data for integrative analysis. After comprehensively comparing adipogenesis and osteogenesis, we quantitatively identified lineage-specific loops and screened out lineage-specific enhancers and open chromatin. We reveal that lineage-specific loops can activate gene expression and facilitate cell commitment through combining enhancers and accessible chromatin in a lineage-specific manner. We finally proposed loop-mediated regulatory networks and identified the controlling factors for adipocytes and osteoblasts determination. Functional experiments validated the lineage-specific regulation networks towards IRS2 and RUNX2 that are associated with adipogenesis and osteogenesis, respectively. These results are expected to help better understand the chromatin conformation determinants of hMSCs fate commitment.

Non-coding variants disrupting a tissue-specific regulatory element in HK1 cause congenital hyperinsulinism

Gene expression is tightly regulated, with many genes exhibiting cell-specific silencing when their protein product would disrupt normal cellular function1. This silencing is largely controlled by non-coding elements, and their disruption might cause human disease2. We performed gene-agnostic screening of the non-coding regions to discover new molecular causes of congenital hyperinsulinism. This identified 14 non-coding de novo variants affecting a 42-bp conserved region encompassed by a regulatory element in intron 2 of the hexokinase 1 gene (HK1). HK1 is widely expressed across all tissues except in the liver and pancreatic beta cells and is thus termed a 'disallowed gene' in these specific tissues. We demonstrated that the variants result in a loss of repression of HK1 in pancreatic beta cells, thereby causing insulin secretion and congenital hyperinsulinism. Using epigenomic data accessed from public repositories, we demonstrated that these variants reside within a regulatory region that we determine to be critical for cell-specific silencing. Importantly, this has revealed a disease mechanism for non-coding variants that cause inappropriate expression of a disallowed gene.

Systematic analysis of the effects of genetic variants on chromatin accessibility to decipher functional variants in non-coding regions

Genome-wide association study (GWAS) has identified thousands of single nucleotide polymorphisms (SNPs) associated with complex diseases and traits. However, deciphering the functions of these SNPs still faces challenges. Recent studies have shown that SNPs could alter chromatin accessibility and result in differences in tumor susceptibility between individuals. Therefore, systematically analyzing the effects of SNPs on chromatin accessibility could help decipher the functions of SNPs, especially those in non-coding regions. Using data from The Cancer Genome Atlas (TCGA), chromatin accessibility quantitative trait locus (caQTL) analysis was conducted to estimate the associations between genetic variants and chromatin accessibility. We analyzed caQTLs in 23 human cancer types and identified 9,478 caQTLs in breast carcinoma (BRCA). In BRCA, these caQTLs tend to alter the binding affinity of transcription factors, and open chromatin regions regulated by these caQTLs are enriched in regulatory elements. By integrating with eQTL data, we identified 141 caQTLs showing a strong signal for colocalization with eQTLs. We also identified 173 caQTLs in genome-wide association studies (GWAS) loci and inferred several possible target genes of these caQTLs. By performing survival analysis, we found that ~10% caQTLs potentially influence the prognosis of patients. To facilitate access to relevant data, we developed a user-friendly data portal, BCaQTL (http://gong_lab.hzau.edu.cn/caqtl_database), for data searching and downloading. Our work may facilitate fine-map regulatory mechanisms underlying risk loci of cancer and discover the biomarkers or therapeutic targets for cancer prognosis. The BCaQTL database will be an important resource for genetic and epigenetic studies.

Prediction of gene co-expression from chromatin contacts with graph attention network

MOTIVATION

The technology of high-throughput chromatin conformation capture (Hi-C) allows genome-wide measurement of chromatin interactions. Several studies have shown statistically significant relationships between gene-gene spatial contacts and their co-expression. It is desirable to uncover epigenetic mechanisms of transcriptional regulation behind such relationships using computational modeling. Existing methods for predicting gene co-expression from Hi-C data use manual feature engineering or unsupervised learning, which either limits the prediction accuracy or lacks interpretability.

RESULTS

To address these issues, we propose HiCoEx (Hi-C predicts gene co-expression), a novel end-to-end framework for explainable prediction of gene co-expression from Hi-C data based on graph neural network. We apply graph attention mechanism to a gene contact network inferred from Hi-C data to distinguish the importance among different neighboring genes of each gene, and learn the gene representation to predict co-expression in a supervised and task-specific manner. Then, from the trained model, we extract the learned gene embeddings as a model interpretation to distill biological insights. Experimental results show that HiCoEx can learn gene representation from 3D genomics signals automatically to improve prediction accuracy, and make the black box model explainable by capturing some biologically meaningful patterns, e.g., in a gene contact network, the common neighbors of two central genes might contribute to the co-expression of the two central genes through sharing enhancers.

AVAILABILITY AND IMPLEMENTATION

The source code is freely available at https://github.com/JieZheng-ShanghaiTech/HiCoEx.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

3D chromatin maps of the human pancreas reveal lineage-specific regulatory architecture of T2D risk

Three-dimensional (3D) chromatin organization maps help dissect cell-type-specific gene regulatory programs. Furthermore, 3D chromatin maps contribute to elucidating the pathogenesis of complex genetic diseases by connecting distal regulatory regions and genetic risk variants to their respective target genes. To understand the cell-type-specific regulatory architecture of diabetes risk, we generated transcriptomic and 3D epigenomic profiles of human pancreatic acinar, alpha, and beta cells using single-cell RNA-seq, single-cell ATAC-seq, and high-resolution Hi-C of sorted cells. Comparisons of these profiles revealed differential A/B (open/closed) chromatin compartmentalization, chromatin looping, and transcriptional factor-mediated control of cell-type-specific gene regulatory programs. We identified a total of 4,750 putative causal-variant-to-target-gene pairs at 194 type 2 diabetes GWAS signals using pancreatic 3D chromatin maps. We found that the connections between candidate causal variants and their putative target effector genes are cell-type stratified and emphasize previously underappreciated roles for alpha and acinar cells in diabetes pathogenesis.

Recent Progress in the Diagnosis and Management of Type 2 Diabetes Mellitus in the Era of COVID-19 and Single Cell Multi-Omics Technologies

Type 2 diabetes mellitus (T2DM) is one of the world's leading causes of death and life-threatening conditions. Therefore, we review the complex vicious circle of causes responsible for T2DM and risk factors such as the western diet, obesity, genetic predisposition, environmental factors, and SARS-CoV-2 infection. The prevalence and economic burden of T2DM on societal and healthcare systems are dissected. Recent progress on the diagnosis and clinical management of T2DM, including both non-pharmacological and latest pharmacological treatment regimens, are summarized. The treatment of T2DM is becoming more complex as new medications are approved. This review is focused on the non-insulin treatments of T2DM to reach optimal therapy beyond glycemic management. We review experimental and clinical findings of SARS-CoV-2 risks that are attributable to T2DM patients. Finally, we shed light on the recent single-cell-based technologies and multi-omics approaches that have reached breakthroughs in the understanding of the pathomechanism of T2DM.

RNA-Binding Proteins in the Regulation of Adipogenesis and Adipose Function

The obesity epidemic represents a critical public health issue worldwide, as it is a vital risk factor for many diseases, including type 2 diabetes (T2D) and cardiovascular disease. Obesity is a complex disease involving excessive fat accumulation. Proper adipose tissue accumulation and function are highly transcriptional and regulated by many genes. Recent studies have discovered that post-transcriptional regulation, mainly mediated by RNA-binding proteins (RBPs), also plays a crucial role. In the lifetime of RNA, it is bound by various RBPs that determine every step of RNA metabolism, from RNA processing to alternative splicing, nucleus export, rate of translation, and finally decay. In humans, it is predicted that RBPs account for more than 10% of proteins based on the presence of RNA-binding domains. However, only very few RBPs have been studied in adipose tissue. The primary aim of this paper is to provide an overview of RBPs in adipogenesis and adipose function. Specifically, the following best-characterized RBPs will be discussed, including HuR, PSPC1, Sam68, RBM4, Ybx1, Ybx2, IGF2BP2, and KSRP. Characterization of these proteins will increase our understanding of the regulatory mechanisms of RBPs in adipogenesis and provide clues for the etiology and pathology of adipose-tissue-related diseases.

The role of insulin-like growth factor 2 mRNA binding proteins in female reproductive pathophysiology

Insulin-like growth factor 2 (IGF2) mRNA binding proteins (IMPs) family belongs to a highly conserved family of RNA-binding proteins (RBPs) and is responsible for regulating RNA processing including localization, translation and stability. Mammalian IMPs (IMP1-3) take part in development, metabolism and tumorigenesis, where they are believed to play a major role in cell growth, metabolism, migration and invasion. IMPs have been identified that are expressed in ovary, placenta and embryo. The up-to-date evidence suggest that IMPs are involved in folliculogenesis, oocyte maturation, embryogenesis, implantation, and placentation. The dysregulation of IMPs not only contributes to carcinogenesis but also disturbs the female reproduction, and may participate in the pathogenesis of reproductive diseases and obstetric syndromes, such as polycystic ovary syndrome (PCOS), pre-eclampsia (PE), gestational diabetes mellitus (GDM) and gynecological tumors. In this review, we summarize the role of IMPs in female reproductive pathophysiology, and hope to provide new insights into the identification of potential therapeutic targets.

Chromatin accessibility profiling by ATAC-seq

The assay for transposase-accessible chromatin using sequencing (ATAC-seq) provides a simple and scalable way to detect the unique chromatin landscape associated with a cell type and how it may be altered by perturbation or disease. ATAC-seq requires a relatively small number of input cells and does not require a priori knowledge of the epigenetic marks or transcription factors governing the dynamics of the system. Here we describe an updated and optimized protocol for ATAC-seq, called Omni-ATAC, that is applicable across a broad range of cell and tissue types. The ATAC-seq workflow has five main steps: sample preparation, transposition, library preparation, sequencing and data analysis. This protocol details the steps to generate and sequence ATAC-seq libraries, with recommendations for sample preparation and downstream bioinformatic analysis. ATAC-seq libraries for roughly 12 samples can be generated in 10 h by someone familiar with basic molecular biology, and downstream sequencing analysis can be implemented using benchmarked pipelines by someone with basic bioinformatics skills and with access to a high-performance computing environment.

Insulin-like Growth Factor-2 Binding Protein-2 Gene Polymorphisms in Iraqi Patients with Type 2 Diabetes Mellitus

Background: Diabetes mellitus type2 (T2DM) represent a hyperglycemia causing metabolic disease which exists in the peripheral tissues due to incomplete pancreatic insulin secretion or insulin resistance. IGF2BP2 is a protein that is involved in embryogenesis and pancreatic development. Genetic association researches had suggested that the single nucleotide polymorphisms (SNP) spanning IGF2BP2 gene are associated with the progression as well as development of the T2DM. Aim: This study aims to evaluate the association of IGF2BP2 gene polymorphisms (rs4402960 & rs1470579) with T2DM in a sample of Iraqi individuals. Methods: A case-control study has been conducted on 800 participants, they were divided to two equal groups, which are a healthy control group (400) and type 2 diabetic patients (400). Fast blood sugar (FBS), total cholesterol (TC), triglyceride (TG), low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), and HbA1c] measured suitable for both participant groups. IGF2BP2 gene has been genotyped for polymorphisms; rs4402960 and rs1470579 by using the PCR-RFLP technique. Results: There is significant changes in the biochemical parameters in patients group when compared to the control group.The SNP rs4402960 show minor allele frequency of T allele considerably different between the two participating groups (p 0.0013) with 33.6 % in T2DM group. Homo-variant TT shows a significant p <0.0001) odd ratio (4.5) as codominant type. Similarly, dominant and recessive models exert significant (0.02 & <0.0001 respectively) adjusted odd ratio (1.45 & 4.14 respectively). The rs1470579 SNP show a significant (0.024) risk (1.28) of C allele in the patients group than in A allele. The CC genotype in codominant and recessive models show significant (0.03) odd ratio differences (2.03 & 1.96 respectively. The rs1470579 SNP exerts significant differences as codominant model in biochemical features of BMI, FBG, Tgs, VLDL-C, insulin and HOMA-IR. The study power of rs4402960 is 69.5% and rs1470579 is 34.1%. Conclusion: This study confirmed the association of rs4402960 as codominant, dominant and recessive with T2DM significantly. However, rs1470579 is associate as recessive model with T2DM in Iraqi population. 

An effector index to predict target genes at GWAS loci

Drug development and biological discovery require effective strategies to map existing genetic associations to causal genes. To approach this problem, we selected 12 common diseases and quantitative traits for which highly powered genome-wide association studies (GWAS) were available. For each disease or trait, we systematically curated positive control gene sets from Mendelian forms of the disease and from targets of medicines used for disease treatment. We found that these positive control genes were highly enriched in proximity of GWAS-associated single-nucleotide variants (SNVs). We then performed quantitative assessment of the contribution of commonly used genomic features, including open chromatin maps, expression quantitative trait loci (eQTL), and chromatin conformation data. Using these features, we trained and validated an Effector Index (Ei), to map target genes for these 12 common diseases and traits. Ei demonstrated high predictive performance, both with cross-validation on the training set, and an independently derived set for type 2 diabetes. Key predictive features included coding or transcript-altering SNVs, distance to gene, and open chromatin-based metrics. This work outlines a simple, understandable approach to prioritize genes at GWAS loci for functional follow-up and drug development, and provides a systematic strategy for prioritization of GWAS target genes.

Identification of Novel Regulatory Regions Induced by Intrauterine Growth Restriction in Rat Islets

Intrauterine growth restriction (IUGR) leads to the development of type 2 diabetes in adulthood, and the permanent alterations in gene expression implicate an epigenetic mechanism. Using a rat model of IUGR, we performed TrueSeq-HELP Tagging to assess the association of DNA methylation changes and gene dysregulation in islets. We identified 511 differentially methylated regions (DMRs) and 4377 significantly altered single CpG sites. Integrating the methylome and our published transcriptome data sets resulted in the identification of pathways critical for islet function. The identified DMRs were enriched with transcription factor binding motifs, such as Elk1, Etv1, Foxa1, Foxa2, Pax7, Stat3, Hnf1, and AR. In silico analysis of 3-dimensional chromosomal interactions using human pancreas and islet Hi-C data sets identified interactions between 14 highly conserved DMRs and 35 genes with significant expression changes at an early age, many of which persisted in adult islets. In adult islets, there were far more interactions between DMRs and genes with significant expression changes identified with Hi-C, and most of them were critical to islet metabolism and insulin secretion. The methylome was integrated with our published genome-wide histone modification data sets from IUGR islets, resulting in further characterization of important regulatory regions of the genome altered by IUGR containing both significant changes in DNA methylation and specific histone marks. We identified novel regulatory regions in islets after exposure to IUGR, suggesting that epigenetic changes at key transcription factor binding motifs and other gene regulatory regions may contribute to gene dysregulation and an abnormal islet phenotype in IUGR rats.

Epigenomic and Transcriptomic Prioritization of Candidate Obesity-Risk Regulatory GWAS SNPs

Concern about rising rates of obesity has prompted searches for obesity-related single nucleotide polymorphisms (SNPs) in genome-wide association studies (GWAS). Identifying plausible regulatory SNPs is very difficult partially because of linkage disequilibrium. We used an unusual epigenomic and transcriptomic analysis of obesity GWAS-derived SNPs in adipose versus heterologous tissues. From 50 GWAS and 121,064 expanded SNPs, we prioritized 47 potential causal regulatory SNPs (Tier-1 SNPs) for 14 gene loci. A detailed examination of seven loci revealed that four (CABLES1, PC, PEMT, and FAM13A) had Tier-1 SNPs positioned so that they could regulate use of alternative transcription start sites, resulting in different polypeptides being generated or different amounts of an intronic microRNA gene being expressed. HOXA11 and long noncoding RNA gene RP11-392O17.1 had Tier-1 SNPs in their 3' or promoter region, respectively, and strong preferences for expression in subcutaneous versus visceral adipose tissue. ZBED3-AS1 had two intragenic Tier-1 SNPs, each of which could contribute to mediating obesity risk through modulating long-distance chromatin interactions. Our approach not only revealed especially credible novel regulatory SNPs, but also helped evaluate previously highlighted obesity GWAS SNPs that were candidates for transcription regulation.

An adipose lncRAP2-Igf2bp2 complex enhances adipogenesis and energy expenditure by stabilizing target mRNAs

lncRAP2 is a conserved cytoplasmic lncRNA enriched in adipose tissue and required for adipogenesis. Using purification and in vivo interactome analyses, we show that lncRAP2 forms complexes with proteins that stabilize mRNAs and modulate translation, among them Igf2bp2. Surveying transcriptome-wide Igf2bp2 client mRNAs in white adipocytes reveals selective binding to mRNAs encoding adipogenic regulators and energy expenditure effectors, including adiponectin. These same target proteins are downregulated when either Igf2bp2 or lncRAP2 is downregulated, hindering adipocyte lipolysis. Proteomics and ribosome profiling show this occurs predominantly through mRNA accumulation, as lncRAP2-Igf2bp2 complex binding does not impact translation efficiency. Phenome-wide association studies reveal specific associations of genetic variants within both lncRAP2 and Igf2bp2 with body mass and type 2 diabetes, and both lncRAP2 and Igf2bp2 are suppressed in adipose depots of obese and diabetic individuals. Thus, the lncRAP2-Igf2bp2 complex potentiates adipose development and energy expenditure and is associated with susceptibility to obesity-linked diabetes.

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.

Multi-Omic Approaches to Identify Genetic Factors in Metabolic Syndrome

Metabolic syndrome (MetS) is a highly heritable disease and a major public health burden worldwide. MetS diagnosis criteria are met by the simultaneous presence of any three of the following: high triglycerides, low HDL/high LDL cholesterol, insulin resistance, hypertension, and central obesity. These diseases act synergistically in people suffering from MetS and dramatically increase risk of morbidity and mortality due to stroke and cardiovascular disease, as well as certain cancers. Each of these component features is itself a complex disease, as is MetS. As a genetically complex disease, genetic risk factors for MetS are numerous, but not very powerful individually, often requiring specific environmental stressors for the disease to manifest. When taken together, all sequence variants that contribute to MetS disease risk explain only a fraction of the heritable variance, suggesting additional, novel loci have yet to be discovered. In this article, we will give a brief overview on the genetic concepts needed to interpret genome-wide association studies (GWAS) and quantitative trait locus (QTL) data, summarize the state of the field of MetS physiological genomics, and to introduce tools and resources that can be used by the physiologist to integrate genomics into their own research on MetS and any of its component features. There is a wealth of phenotypic and molecular data in animal models and humans that can be leveraged as outlined in this article. Integrating these multi-omic QTL data for complex diseases such as MetS provides a means to unravel the pathways and mechanisms leading to complex disease and promise for novel treatments. © 2022 American Physiological Society. Compr Physiol 12:1-40, 2022.

Enhancing Acsl4 in absence of mTORC2/Rictor drove β-cell dedifferentiation via inhibiting FoxO1 and promoting ROS production

Rapamycin insensitive companion of mechanistic target of Rapamycin (Rictor), the key component of mTOR complex 2 (mTORC2), controls both β-cell proliferation and function. We sought to study whether long chain acyl-CoA synthetase 4 (Acsl4) worked downstream of Rictor/mTORC2 to maintain β-cell functional mass. We found Acsl4 was positively regulated by Rictor at transcriptional and posttranslational levels in mouse β-cell. Infecting adenovirus expressing Acsl4 in β-cell-specific-Rictor-knockout (βRicKO) islets and Min6 cells knocking down Rictor with lentivirus-expressing siRNA-oligos targeting Rictor(siRic), recovered the β-cell dysplasia but not dysfunction. Cell bioenergetic experiment performed with Seahorse XF showed that Acsl4 could not rescue the dampened glucose oxidation in Rictor-lacking β-cell, but further promoted lipid oxidation. Transposase-Accessible Chromatin (ATAC) and H3K27Ac chromatin immunoprecipitation (ChIP) sequencing studies reflected the epigenetic elevated molecular signature for β-cell dedifferentiation and mitigated oxidative defense/response. These results were confirmed by the observations of elevated acetylation and ubiquitination of FoxO1, increased protein levels of Gpx1 and Hif1an, excessive reactive oxygen species (ROS) production and diminished MafA in Acsl4 overexpressed Rictor-lacking β-cells. In these cells, antioxidant treatment significantly recovered MafA level and insulin content. Inducing lipid oxidation alone could not mimic the effect of Acsl4 in Rictor lacking β-cell. Our study suggested that Acsl4 function in β-cell was context dependent and might facilitate β-cell dedifferentiation with attenuated Rictor/mTORC2 activity or insulin signaling via posttranslational inhibiting FoxO1 and epigenetically enhancing ROS induced MafA degradation.

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.