A Chromatin Basis for Cell Lineage and Disease Risk in the Human Pancreas

Cell type- and transcription-independent spatial proximity between enhancers and promoters

Cell type-specific enhancers are critically important for lineage specification. The mechanisms that determine cell-type specificity of enhancer activity, however, are not fully understood. Most current models for how enhancers function invoke physical proximity between enhancer elements and their target genes. Here, we use an imaging-based approach to examine the spatial relationship of cell type-specific enhancers and their target genes with single-cell resolution. Using high-throughput microscopy, we measure the spatial distance from target promoters to their cell type-specific active and inactive enhancers in individual pancreatic cells derived from distinct lineages. We find increased proximity of all promoter-enhancer pairs relative to non-enhancer pairs separated by similar genomic distances. Strikingly, spatial proximity between enhancers and target genes was unrelated to tissue-specific enhancer activity. Furthermore, promoter-enhancer proximity did not correlate with the expression status of target genes. Our results suggest that promoter-enhancer pairs exist in a distinctive chromatin environment but that genome folding is not a universal driver of cell-type specificity in enhancer function.

Genome-Wide Analysis Identifies Nuclear Factor 1C as a Novel Transcription Factor and Potential Therapeutic Target in SCLC

INTRODUCTION

Recent insights regarding mechanisms mediating stemness, heterogeneity, and metastatic potential of lung cancers have yet to be fully translated to effective regimens for the treatment of these malignancies. This study sought to identify novel targets for lung cancer therapy.

METHODS

Transcriptomes and DNA methylomes of 14 SCLC and 10 NSCLC lines were compared with normal human small airway epithelial cells (SAECs) and induced pluripotent stem cell (iPSC) clones derived from SAEC. SCLC lines, lung iPSC (Lu-iPSC), and SAEC were further evaluated by DNase I hypersensitive site sequencing (DHS-seq). Changes in chromatin accessibility and depths of transcription factor (TF) footprints were quantified using Bivariate analysis of Genomic Footprint. Standard techniques were used to evaluate growth, tumorigenicity, and changes in transcriptomes and glucose metabolism of SCLC cells after NFIC knockdown and to evaluate NFIC expression in SCLC cells after exposure to BET inhibitors.

RESULTS

Considerable commonality of transcriptomes and DNA methylomes was observed between Lu-iPSC and SCLC; however, this analysis was uninformative regarding pathways unique to lung cancer. Linking results of DHS-seq to RNA sequencing enabled identification of networks not previously associated with SCLC. When combined with footprint depth, NFIC, a transcription factor not previously associated with SCLC, had the highest score of occupancy at open chromatin sites. Knockdown of NFIC impaired glucose metabolism, decreased stemness, and inhibited growth of SCLC cells in vitro and in vivo. ChIP-seq analysis identified numerous sites occupied by BRD4 in the NFIC promoter region. Knockdown of BRD4 or treatment with Bromodomain and extra-terminal domain (BET) inhibitors (BETis) markedly reduced NFIC expression in SCLC cells and SCLC PDX models. Approximately 8% of genes down-regulated by BETi treatment were repressed by NFIC knockdown in SCLC, whereas 34% of genes repressed after NFIC knockdown were also down-regulated in SCLC cells after BETi treatment.

CONCLUSIONS

NFIC is a key TF and possible mediator of transcriptional regulation by BET family proteins in SCLC. Our findings highlight the potential of genome-wide chromatin accessibility analysis for elucidating mechanisms of pulmonary carcinogenesis and identifying novel targets for lung cancer therapy.

PACS allows comprehensive dissection of multiple factors governing chromatin accessibility from snATAC-seq data

Single nucleus ATAC-seq (snATAC-seq) experimental designs have become increasingly complex with multiple factors that might affect chromatin accessibility, including genotype, cell type, tissue of origin, sample location, batch, etc., whose compound effects are difficult to test by existing methods. In addition, current snATAC-seq data present statistical difficulties due to their sparsity and variations in individual sequence capture. To address these problems, we present a zero-adjusted statistical model, Probability model of Accessible Chromatin of Single cells (PACS), that can allow complex hypothesis testing of factors that affect accessibility while accounting for sparse and incomplete data. For differential accessibility analysis, PACS controls the false positive rate and achieves on average a 17% to 122% higher power than existing tools. We demonstrate the effectiveness of PACS through several analysis tasks including supervised cell type annotation, compound hypothesis testing, batch effect correction, and spatiotemporal modeling. We apply PACS to several datasets from a variety of tissues and show its ability to reveal previously undiscovered insights in snATAC-seq data.

Glucagon-producing α-cell transcriptional identity and reprogramming towards insulin production

β-Cell replacement by in situ reprogramming of non-β-cells is a promising diabetes therapy. Following the observation that near-total β-cell ablation in adult mice triggers the reprogramming of pancreatic α-, δ-, and γ-cells into insulin (INS)-producing cells, recent studies are delving deep into the mechanisms controlling adult α-cell identity. Systematic analyses of the α-cell transcriptome and epigenome have started to pinpoint features that could be crucial for maintaining α-cell identity. Using different transgenic and chemical approaches, significant advances have been made in reprogramming α-cells in vivo into INS-secreting cells in mice. The recent reprogramming of human α-cells in vitro is an important step forward that must now be complemented with a comprehensive molecular dissection of the mechanisms controlling α-cell identity.

PancrESS - a meta-analysis resource for understanding cell-type specific expression in the human pancreas

BACKGROUND

The human pancreas is composed of specialized cell types producing hormones and enzymes critical to human health. These specialized functions are the result of cell type-specific transcriptional programs which manifest in cell-specific gene expression. Understanding these programs is essential to developing therapies for pancreatic disorders. Transcription in the human pancreas has been widely studied by single-cell RNA technologies, however the diversity of protocols and analysis methods hinders their interpretability in the aggregate.

RESULTS

In this work, we perform a meta-analysis of pancreatic single-cell RNA sequencing data. We present a database for reference transcriptome abundances and cell-type specificity metrics. This database facilitates the identification and definition of marker genes within the pancreas. Additionally, we introduce a versatile tool which is freely available as an R package, and should permit integration into existing workflows. Our tool accepts count data files generated by widely-used single-cell gene expression platforms in their original format, eliminating an additional pre-formatting step. Although we designed it to calculate expression specificity of pancreas cell types, our tool is agnostic to the biological source of count data, extending its applicability to other biological systems.

CONCLUSIONS

Our findings enhance the current understanding of expression specificity within the pancreas, surpassing previous work in terms of scope and detail. Furthermore, our database and tool enable researchers to perform similar calculations in diverse biological systems, expanding the applicability of marker gene identification and facilitating comparative analyses.

Inferring regulators of cell identity in the human adult pancreas

Cellular identity during development is under the control of transcription factors that form gene regulatory networks. However, the transcription factors and gene regulatory networks underlying cellular identity in the human adult pancreas remain largely unexplored. Here, we integrate multiple single-cell RNA-sequencing datasets of the human adult pancreas, totaling 7393 cells, and comprehensively reconstruct gene regulatory networks. We show that a network of 142 transcription factors forms distinct regulatory modules that characterize pancreatic cell types. We present evidence that our approach identifies regulators of cell identity and cell states in the human adult pancreas. We predict that HEYL, BHLHE41 and JUND are active in acinar, beta and alpha cells, respectively, and show that these proteins are present in the human adult pancreas as well as in human induced pluripotent stem cell (hiPSC)-derived islet cells. Using single-cell transcriptomics, we found that JUND represses beta cell genes in hiPSC-alpha cells. BHLHE41 depletion induced apoptosis in primary pancreatic islets. The comprehensive gene regulatory network atlas can be explored interactively online. We anticipate our analysis to be the starting point for a more sophisticated dissection of how transcription factors regulate cell identity and cell states in the human adult pancreas.

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.

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.

Mechanism of activation and the rewired network: New drug design concepts

Precision oncology benefits from effective early phase drug discovery decisions. Recently, drugging inactive protein conformations has shown impressive successes, raising the cardinal questions of which targets can profit and what are the principles of the active/inactive protein pharmacology. Cancer driver mutations have been established to mimic the protein activation mechanism. We suggest that the decision whether to target an inactive (or active) conformation should largely rest on the protein mechanism of activation. We next discuss the recent identification of double (multiple) same-allele driver mutations and their impact on cell proliferation and suggest that like single driver mutations, double drivers also mimic the mechanism of activation. We further suggest that the structural perturbations of double (multiple) in cis mutations may reveal new surfaces/pockets for drug design. Finally, we underscore the preeminent role of the cellular network which is deregulated in cancer. Our structure-based review and outlook updates the traditional Mechanism of Action, informs decisions, and calls attention to the intrinsic activation mechanism of the target protein and the rewired tumor-specific network, ushering innovative considerations in precision medicine.

Role of Enhancers in Development and Diseases

Enhancers are cis-regulatory elements containing short DNA sequences that serve as binding sites for pioneer/regulatory transcription factors, thus orchestrating the regulation of genes critical for lineage determination. The activity of enhancer elements is believed to be determined by transcription factor binding, thus determining the cell state identity during development. Precise spatio-temporal control of the transcriptome during lineage specification requires the coordinated binding of lineage-specific transcription factors to enhancers. Thus, enhancers are the primary determinants of cell identity. Numerous studies have explored the role and mechanism of enhancers during development and disease, and various basic questions related to the functions and mechanisms of enhancers have not yet been fully answered. In this review, we discuss the recently published literature regarding the roles of enhancers, which are critical for various biological processes governing development. Furthermore, we also highlight that altered enhancer landscapes provide an essential context to understand the etiologies and mechanisms behind numerous complex human diseases, providing new avenues for effective enhancer-based therapeutic interventions.

The Human Islet: Mini-Organ With Mega-Impact

This review focuses on the human pancreatic islet-including its structure, cell composition, development, function, and dysfunction. After providing a historical timeline of key discoveries about human islets over the past century, we describe new research approaches and technologies that are being used to study human islets and how these are providing insight into human islet physiology and pathophysiology. We also describe changes or adaptations in human islets in response to physiologic challenges such as pregnancy, aging, and insulin resistance and discuss islet changes in human diabetes of many forms. We outline current and future interventions being developed to protect, restore, or replace human islets. The review also highlights unresolved questions about human islets and proposes areas where additional research on human islets is needed.

Epigenetic Regulation of β Cell Identity and Dysfunction

β cell dysfunction and failure are driving forces of type 2 diabetes mellitus (T2DM) pathogenesis. Investigating the underlying mechanisms of β cell dysfunction may provide novel targets for the development of next generation therapy for T2DM. Epigenetics is the study of gene expression changes that do not involve DNA sequence changes, including DNA methylation, histone modification, and non-coding RNAs. Specific epigenetic signatures at all levels, including DNA methylation, chromatin accessibility, histone modification, and non-coding RNA, define β cell identity during embryonic development, postnatal maturation, and maintain β cell function at homeostatic states. During progression of T2DM, overnutrition, inflammation, and other types of stress collaboratively disrupt the homeostatic epigenetic signatures in β cells. Dysregulated epigenetic signatures, and the associating transcriptional outputs, lead to the dysfunction and eventual loss of β cells. In this review, we will summarize recent discoveries of the establishment and disruption of β cell-specific epigenetic signatures, and discuss the potential implication in therapeutic development.

A Transcriptionally Distinct Subpopulation of Healthy Acinar Cells Exhibit Features of Pancreatic Progenitors and PDAC

Pancreatic ductal adenocarcinoma (PDAC) tumors can originate either from acinar or ductal cells in the adult pancreas. We re-analyze multiple pancreas and PDAC single-cell RNA-seq datasets and find a subset of nonmalignant acinar cells, which we refer to as acinar edge (AE) cells, whose transcriptomes highly diverge from a typical acinar cell in each dataset. Genes upregulated among AE cells are enriched for transcriptomic signatures of pancreatic progenitors, acinar dedifferentiation, and several oncogenic programs. AE-upregulated genes are upregulated in human PDAC tumors, and consistently, their promoters are hypomethylated. High expression of these genes is associated with poor patient survival. The fraction of AE-like cells increases with age in healthy pancreatic tissue, which is not explained by clonal mutations, thus pointing to a nongenetic source of variation. The fraction of AE-like cells is also significantly higher in human pancreatitis samples. Finally, we find edge-like states in lung, liver, prostate, and colon tissues, suggesting that subpopulations of healthy cells across tissues can exist in pre-neoplastic states. SIGNIFICANCE: These findings show "edge" epithelial cell states with oncogenic transcriptional activity in human organs without oncogenic mutations. In the pancreas, the fraction of acinar cells increases with age.

Interpreting type 1 diabetes risk with genetics and single-cell epigenomics

Genetic risk variants that have been identified in genome-wide association studies of complex diseases are primarily non-coding1. Translating these risk variants into mechanistic insights requires detailed maps of gene regulation in disease-relevant cell types2. Here we combined two approaches: a genome-wide association study of type 1 diabetes (T1D) using 520,580 samples, and the identification of candidate cis-regulatory elements (cCREs) in pancreas and peripheral blood mononuclear cells using single-nucleus assay for transposase-accessible chromatin with sequencing (snATAC-seq) of 131,554 nuclei. Risk variants for T1D were enriched in cCREs that were active in T cells and other cell types, including acinar and ductal cells of the exocrine pancreas. Risk variants at multiple T1D signals overlapped with exocrine-specific cCREs that were linked to genes with exocrine-specific expression. At the CFTR locus, the T1D risk variant rs7795896 mapped to a ductal-specific cCRE that regulated CFTR; the risk allele reduced transcription factor binding, enhancer activity and CFTR expression in ductal cells. These findings support a role for the exocrine pancreas in the pathogenesis of T1D and highlight the power of large-scale genome-wide association studies and single-cell epigenomics for understanding the cellular origins of complex disease.

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.

A multilayered post-GWAS assessment on genetic susceptibility to pancreatic cancer

BACKGROUND

Pancreatic cancer (PC) is a complex disease in which both non-genetic and genetic factors interplay. To date, 40 GWAS hits have been associated with PC risk in individuals of European descent, explaining 4.1% of the phenotypic variance.

METHODS

We complemented a new conventional PC GWAS (1D) with genome spatial autocorrelation analysis (2D) permitting to prioritize low frequency variants not detected by GWAS. These were further expanded via Hi-C map (3D) interactions to gain additional insight into the inherited basis of PC. In silico functional analysis of public genomic information allowed prioritization of potentially relevant candidate variants.

RESULTS

We identified several new variants located in genes for which there is experimental evidence of their implication in the biology and function of pancreatic acinar cells. Among them is a novel independent variant in NR5A2 (rs3790840) with a meta-analysis p value = 5.91E-06 in 1D approach and a Local Moran's Index (LMI) = 7.76 in 2D approach. We also identified a multi-hit region in CASC8-a lncRNA associated with pancreatic carcinogenesis-with a lowest p value = 6.91E-05. Importantly, two new PC loci were identified both by 2D and 3D approaches: SIAH3 (LMI = 18.24), CTRB2/BCAR1 (LMI = 6.03), in addition to a chromatin interacting region in XBP1-a major regulator of the ER stress and unfolded protein responses in acinar cells-identified by 3D; all of them with a strong in silico functional support.

CONCLUSIONS

This multi-step strategy, combined with an in-depth in silico functional analysis, offers a comprehensive approach to advance the study of PC genetic susceptibility and could be applied to other diseases.

SIX2 and SIX3 coordinately regulate functional maturity and fate of human pancreatic β cells

The physiological functions of many vital tissues and organs continue to mature after birth, but the genetic mechanisms governing this postnatal maturation remain an unsolved mystery. Human pancreatic β cells produce and secrete insulin in response to physiological cues like glucose, and these hallmark functions improve in the years after birth. This coincides with expression of the transcription factors SIX2 and SIX3, whose functions in native human β cells remain unknown. Here, we show that shRNA-mediated SIX2 or SIX3 suppression in human pancreatic adult islets impairs insulin secretion. However, transcriptome studies revealed that SIX2 and SIX3 regulate distinct targets. Loss of SIX2 markedly impaired expression of genes governing β-cell insulin processing and output, glucose sensing, and electrophysiology, while SIX3 loss led to inappropriate expression of genes normally expressed in fetal β cells, adult α cells, and other non-β cells. Chromatin accessibility studies identified genes directly regulated by SIX2. Moreover, β cells from diabetic humans with impaired insulin secretion also had reduced SIX2 transcript levels. Revealing how SIX2 and SIX3 govern functional maturation and maintain developmental fate in native human β cells should advance β-cell replacement and other therapeutic strategies for diabetes.

A gene-environment-induced epigenetic program initiates tumorigenesis

Tissue damage increases the risk of cancer through poorly understood mechanisms1. In mouse models of pancreatic cancer, pancreatitis associated with tissue injury collaborates with activating mutations in the Kras oncogene to markedly accelerate the formation of early neoplastic lesions and, ultimately, adenocarcinoma2,3. Here, by integrating genomics, single-cell chromatin assays and spatiotemporally controlled functional perturbations in autochthonous mouse models, we show that the combination of Kras mutation and tissue damage promotes a unique chromatin state in the pancreatic epithelium that distinguishes neoplastic transformation from normal regeneration and is selected for throughout malignant evolution. This cancer-associated epigenetic state emerges within 48 hours of pancreatic injury, and involves an 'acinar-to-neoplasia' chromatin switch that contributes to the early dysregulation of genes that define human pancreatic cancer. Among the factors that are most rapidly activated after tissue damage in the pre-malignant pancreatic epithelium is the alarmin cytokine interleukin 33, which recapitulates the effects of injury in cooperating with mutant Kras to unleash the epigenetic remodelling program of early neoplasia and neoplastic transformation. Collectively, our study demonstrates how gene-environment interactions can rapidly produce gene-regulatory programs that dictate early neoplastic commitment, and provides a molecular framework for understanding the interplay between genetic and environmental cues in the initiation of cancer.

Aggressive PDACs Show Hypomethylation of Repetitive Elements and the Execution of an Intrinsic IFN Program Linked to a Ductal Cell of Origin

Pancreatic ductal adenocarcinoma (PDAC) is characterized by extensive desmoplasia, which challenges the molecular analyses of bulk tumor samples. Here we FACS-purified epithelial cells from human PDAC and normal pancreas and derived their genome-wide transcriptome and DNA methylome landscapes. Clustering based on DNA methylation revealed two distinct PDAC groups displaying different methylation patterns at regions encoding repeat elements. Methylation^low tumors are characterized by higher expression of endogenous retroviral transcripts and double-stranded RNA sensors, which lead to a cell-intrinsic activation of an interferon signature (IFNsign). This results in a protumorigenic microenvironment and poor patient outcome. Methylation^low/IFNsign^high and Methylation^high/IFNsign^low PDAC cells preserve lineage traits, respective of normal ductal or acinar pancreatic cells. Moreover, ductal-derived Kras ^G12D/Trp53 ^-/- mouse PDACs show higher expression of IFNsign compared with acinar-derived counterparts. Collectively, our data point to two different origins and etiologies of human PDACs, with the aggressive Methylation^low/IFNsign^high subtype potentially targetable by agents blocking intrinsic IFN signaling. SIGNIFICANCE: The mutational landscapes of PDAC alone cannot explain the observed interpatient heterogeneity. We identified two PDAC subtypes characterized by differential DNA methylation, preserving traits from normal ductal/acinar cells associated with IFN signaling. Our work suggests that epigenetic traits and the cell of origin contribute to PDAC heterogeneity.This article is highlighted in the In This Issue feature, p. 521.

(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.

In vivo studies of glucagon secretion by human islets transplanted in mice

Little is known about regulated glucagon secretion by human islet α-cells compared to insulin secretion from β-cells, despite conclusive evidence of dysfunction in both cell types in diabetes mellitus. Distinct insulins in humans and mice permit in vivo studies of human β-cell regulation after human islet transplantation in immunocompromised mice, whereas identical glucagon sequences prevent analogous in vivo measures of glucagon output from human α-cells. Here, we use CRISPR-Cas9 editing to remove glucagon codons 2-29 in immunocompromised NSG mice, preserving the production of other proglucagon-derived hormones. Glucagon knockout NSG (GKO-NSG) mice have metabolic, liver and pancreatic phenotypes associated with glucagon-signalling deficits that revert after transplantation of human islets from non-diabetic donors. Glucagon hypersecretion by transplanted islets from donors with type 2 diabetes revealed islet-intrinsic defects. We suggest that GKO-NSG mice provide an unprecedented resource to investigate human α-cell regulation in vivo.

Are Parallel Proliferation Pathways Redundant?

Are the receptor tyrosine kinase (RTK) and JAK-STAT-driven proliferation pathways 'parallel' or 'redundant'? And what about those of K-Ras4B versus N-Ras? 'Parallel' proliferation pathways accomplish a similar drug resistance outcome. Thus, are they 'redundant'? In this paper, it is argued that there is a fundamental distinction between 'parallel' and 'redundant'. Cellular proliferation pathways are influenced by the genome sequence, 3D organization and chromatin accessibility, and determined by protein availability prior to cancer emergence. In the opinion presented, if they operate the same downstream protein families, they are redundant; if evolutionary-independent, they are parallel. Thus, RTK and JAK-STAT-driven proliferation pathways are parallel; those of Ras isoforms are redundant. Our Precision Medicine Call to map cancer proliferation pathways is vastly important since it can expedite effective therapeutics.

Molecular and genetic regulation of pig pancreatic islet cell development

Reliance on rodents for understanding pancreatic genetics, development and islet function could limit progress in developing interventions for human diseases such as diabetes mellitus. Similarities of pancreas morphology and function suggest that porcine and human pancreas developmental biology may have useful homologies. However, little is known about pig pancreas development. To fill this knowledge gap, we investigated fetal and neonatal pig pancreas at multiple, crucial developmental stages using modern experimental approaches. Purification of islet β-, α- and δ-cells followed by transcriptome analysis (RNA-seq) and immunohistology identified cell- and stage-specific regulation, and revealed that pig and human islet cells share characteristic features that are not observed in mice. Morphometric analysis also revealed endocrine cell allocation and architectural similarities between pig and human islets. Our analysis unveiled scores of signaling pathways linked to native islet β-cell functional maturation, including evidence of fetal α-cell GLP-1 production and signaling to β-cells. Thus, the findings and resources detailed here show how pig pancreatic islet studies complement other systems for understanding the developmental programs that generate functional islet cells, and that are relevant to human pancreatic diseases.

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.

Enhancer signatures stratify and predict outcomes of non-functional pancreatic neuroendocrine tumors

Most pancreatic neuroendocrine tumors (PNETs) do not produce excess hormones and are therefore considered 'non-functional'^1-3. As clinical behaviors vary widely and distant metastases are eventually lethal^2,4, biological classifications might guide treatment. Using enhancer maps to infer gene regulatory programs, we find that non-functional PNETs fall into two major subtypes, with epigenomes and transcriptomes that partially resemble islet α- and β-cells. Transcription factors ARX and PDX1 specify these normal cells, respectively^5,6, and 84% of 142 non-functional PNETs expressed one or the other factor, occasionally both. Among 103 cases, distant relapses occurred almost exclusively in patients with ARX⁺PDX1^- tumors and, within this subtype, in cases with alternative lengthening of telomeres. These markedly different outcomes belied similar clinical presentations and histology and, in one cohort, occurred irrespective of MEN1 mutation. This robust molecular stratification provides insight into cell lineage correlates of non-functional PNETs, accurately predicts disease course and can inform postoperative clinical decisions.