Extensive NEUROG3 occupancy in the human pancreatic endocrine gene regulatory network

Generation and application of novel hES cell reporter lines for the differentiation and maturation of hPS cell-derived islet-like clusters

The significant advances in the differentiation of human pluripotent stem (hPS) cells into pancreatic endocrine cells, including functional β-cells, have been based on a detailed understanding of the underlying developmental mechanisms. However, the final differentiation steps, leading from endocrine progenitors to mono-hormonal and mature pancreatic endocrine cells, remain to be fully understood and this is reflected in the remaining shortcomings of the hPS cell-derived islet cells (SC-islet cells), which include a lack of β-cell maturation and variability among different cell lines. Additional signals and modifications of the final differentiation steps will have to be assessed in a combinatorial manner to address the remaining issues and appropriate reporter lines would be useful in this undertaking. Here we report the generation and functional validation of hPS cell reporter lines that can monitor the generation of INS+ and GCG+ cells and their resolution into mono-hormonal cells (INSeGFP, INSeGFP/GCGmCHERRY) as well as β-cell maturation (INSeGFP/MAFAmCHERRY) and function (INSGCaMP6). The reporter hPS cell lines maintained strong and widespread expression of pluripotency markers and differentiated efficiently into definitive endoderm and pancreatic progenitor (PP) cells. PP cells from all lines differentiated efficiently into islet cell clusters that robustly expressed the corresponding reporters and contained glucose-responsive, insulin-producing cells. To demonstrate the applicability of these hPS cell reporter lines in a high-content live imaging approach for the identification of optimal differentiation conditions, we adapted our differentiation procedure to generate SC-islet clusters in microwells. This allowed the live confocal imaging of multiple SC-islets for a single condition and, using this approach, we found that the use of the N21 supplement in the last stage of the differentiation increased the number of monohormonal β-cells without affecting the number of α-cells in the SC-islets. The hPS cell reporter lines and the high-content live imaging approach described here will enable the efficient assessment of multiple conditions for the optimal differentiation and maturation of SC-islets.

Directed differentiation of pancreatic δ cells from human pluripotent stem cells

Dysfunction of pancreatic δ cells contributes to the etiology of diabetes. Despite their important role, human δ cells are scarce, limiting physiological studies and drug discovery targeting δ cells. To date, no directed δ-cell differentiation method has been established. Here, we demonstrate that fibroblast growth factor (FGF) 7 promotes pancreatic endoderm/progenitor differentiation, whereas FGF2 biases cells towards the pancreatic δ-cell lineage via FGF receptor 1. We develop a differentiation method to generate δ cells from human stem cells by combining FGF2 with FGF7, which synergistically directs pancreatic lineage differentiation and modulates the expression of transcription factors and SST activators during endoderm/endocrine precursor induction. These δ cells display mature RNA profiles and fine secretory granules, secrete somatostatin in response to various stimuli, and suppress insulin secretion from in vitro co-cultured β cells and mouse β cells upon transplantation. The generation of human pancreatic δ cells from stem cells in vitro would provide an unprecedented cell source for drug discovery and cell transplantation studies in diabetes.

Hippo pathway-mediated YAP1/TAZ inhibition is essential for proper pancreatic endocrine specification and differentiation

The Hippo pathway plays a central role in tissue development and homeostasis. However, the function of Hippo in pancreatic endocrine development remains obscure. Here, we generated novel conditional genetically engineered mouse models to examine the roles of Hippo pathway-mediated YAP1/TAZ inhibition in the development stages of endocrine specification and differentiation. While YAP1 protein was localized to the nuclei in bipotent progenitor cells, Neurogenin 3 expressing endocrine progenitors completely lost YAP1 expression. Using genetically engineered mouse models, we found that inactivation of YAP1 requires both an intact Hippo pathway and Neurogenin 3 protein. Gene deletion of Lats1 and 2 kinases (Lats1&2) in endocrine progenitor cells of developing mouse pancreas using Neurog3Cre blocked endocrine progenitor cell differentiation and specification, resulting in reduced islets size and a disorganized pancreas at birth. Loss of Lats1&2 in Neurogenin 3 expressing cells activated YAP1/TAZ transcriptional activity and recruited macrophages to the developing pancreas. These defects were rescued by deletion of Yap1/Wwtr1 genes, suggesting that tight regulation of YAP1/TAZ by Hippo signaling is crucial for pancreatic endocrine specification. In contrast, deletion of Lats1&2 using β-cell-specific Ins1CreER resulted in a phenotypically normal pancreas, indicating that Lats1&2 are indispensable for differentiation of endocrine progenitors but not for that of β-cells. Our results demonstrate that loss of YAP1/TAZ expression in the pancreatic endocrine compartment is not a passive consequence of endocrine specification. Rather, Hippo pathway-mediated inhibition of YAP1/TAZ in endocrine progenitors is a prerequisite for endocrine specification and differentiation.

Targeting β-Cell Plasticity: A Promising Approach for Diabetes Treatment

The β-cells within the pancreas play a pivotal role in insulin production and secretion, responding to fluctuations in blood glucose levels. However, factors like obesity, dietary habits, and prolonged insulin resistance can compromise β-cell function, contributing to the development of Type 2 Diabetes (T2D). A critical aspect of this dysfunction involves β-cell dedifferentiation and transdifferentiation, wherein these cells lose their specialized characteristics and adopt different identities, notably transitioning towards progenitor or other pancreatic cell types like α-cells. This process significantly contributes to β-cell malfunction and the progression of T2D, often surpassing the impact of outright β-cell loss. Alterations in the expressions of specific genes and transcription factors unique to β-cells, along with epigenetic modifications and environmental factors such as inflammation, oxidative stress, and mitochondrial dysfunction, underpin the occurrence of β-cell dedifferentiation and the onset of T2D. Recent research underscores the potential therapeutic value for targeting β-cell dedifferentiation to manage T2D effectively. In this review, we aim to dissect the intricate mechanisms governing β-cell dedifferentiation and explore the therapeutic avenues stemming from these insights.

Interpretable trajectory inference with single-cell Linear Adaptive Negative-binomial Expression (scLANE) testing

The rapid proliferation of trajectory inference methods for single-cell RNA-seq data has allowed researchers to investigate complex biological processes by examining underlying gene expression dynamics. After estimating a latent cell ordering, statistical models are used to determine which genes exhibit changes in expression that are significantly associated with progression through the biological trajectory. While a few techniques for performing trajectory differential expression exist, most rely on the flexibility of generalized additive models in order to account for the inherent nonlinearity of changes in gene expression. As such, the results can be difficult to interpret, and biological conclusions often rest on subjective visual inspections of the most dynamic genes. To address this challenge, we propose scLANE testing, which is built around an interpretable generalized linear model and handles nonlinearity with basis splines chosen empirically for each gene. In addition, extensions to estimating equations and mixed models allow for reliable trajectory testing under complex experimental designs. After validating the accuracy of scLANE under several different simulation scenarios, we apply it to a set of diverse biological datasets and display its ability to provide novel biological information when used downstream of both pseudotime and RNA velocity estimation methods.

Integrating single-cell imaging and RNA sequencing datasets links differentiation and morphogenetic dynamics of human pancreatic endocrine progenitors

Basic helix-loop-helix genes, particularly proneural genes, are well-described triggers of cell differentiation, yet information on their dynamics is limited, notably in human development. Here, we focus on Neurogenin 3 (NEUROG3), which is crucial for pancreatic endocrine lineage initiation. By monitoring both NEUROG3 gene expression and protein in single cells using a knockin dual reporter in 2D and 3D models of human pancreas development, we show an approximately 2-fold slower expression of human NEUROG3 than that of the mouse. We observe heterogeneous peak levels of NEUROG3 expression and reveal through long-term live imaging that both low and high NEUROG3 peak levels can trigger differentiation into hormone-expressing cells. Based on fluorescence intensity, we statistically integrate single-cell transcriptome with dynamic behaviors of live cells and propose a data-mapping methodology applicable to other contexts. Using this methodology, we identify a role for KLK12 in motility at the onset of NEUROG3 expression.

NEUROD2 function is dispensable for human pancreatic β cell specification

Introduction

The molecular programs regulating human pancreatic endocrine cell induction and fate allocation are not well deciphered. Here, we investigated the spatiotemporal expression pattern and the function of the neurogenic differentiation factor 2 (NEUROD2) during human endocrinogenesis.

Methods

Using Crispr-Cas9 gene editing, we generated a reporter knock-in transcription factor (TF) knock-out human inducible pluripotent stem cell (iPSC) line in which the open reading frame of both NEUROD2 alleles are replaced by a nuclear histone 2B-Venus reporter (NEUROD2nVenus/nVenus).

Results

We identified a transient expression of NEUROD2 mRNA and its nuclear Venus reporter activity at the stage of human endocrine progenitor formation in an iPSC differentiation model. This expression profile is similar to what was previously reported in mice, uncovering an evolutionarily conserved gene expression pattern of NEUROD2 during endocrinogenesis. In vitro differentiation of the generated homozygous NEUROD2nVenus/nVenus iPSC line towards human endocrine lineages uncovered no significant impact upon the loss of NEUROD2 on endocrine cell induction. Moreover, analysis of endocrine cell specification revealed no striking changes in the generation of insulin-producing b cells and glucagon-secreting a cells upon lack of NEUROD2.

Discussion

Overall, our results suggest that NEUROD2 is expendable for human b cell formation in vitro.

NEUROD1 reinforces endocrine cell fate acquisition in pancreatic development

NEUROD1 is a transcription factor that helps maintain a mature phenotype of pancreatic β cells. Disruption of Neurod1 during pancreatic development causes severe neonatal diabetes; however, the exact role of NEUROD1 in the differentiation programs of endocrine cells is unknown. Here, we report a crucial role of the NEUROD1 regulatory network in endocrine lineage commitment and differentiation. Mechanistically, transcriptome and chromatin landscape analyses demonstrate that Neurod1 inactivation triggers a downregulation of endocrine differentiation transcription factors and upregulation of non-endocrine genes within the Neurod1-deficient endocrine cell population, disturbing endocrine identity acquisition. Neurod1 deficiency altered the H3K27me3 histone modification pattern in promoter regions of differentially expressed genes, which resulted in gene regulatory network changes in the differentiation pathway of endocrine cells, compromising endocrine cell potential, differentiation, and functional properties.

Characterization of cell-fate decision landscapes by estimating transcription factor dynamics

Time-specific modulation of gene expression during differentiation by transcription factors promotes cell diversity. However, estimating their dynamic regulatory activity at the single-cell level and in a high-throughput manner remains challenging. We present FateCompass, an integrative approach that utilizes single-cell transcriptomics data to identify lineage-specific transcription factors throughout differentiation. By combining a probabilistic framework with RNA velocities or differentiation potential, we estimate transition probabilities, while a linear model of gene regulation is employed to compute transcription factor activities. Considering dynamic changes and correlations of expression and activities, FateCompass identifies lineage-specific regulators. Our validation using in silico data and application to pancreatic endocrine cell differentiation datasets highlight both known and potentially novel lineage-specific regulators. Notably, we uncovered undescribed transcription factors of an enterochromaffin-like population during in vitro differentiation toward ß-like cells. FateCompass provides a valuable framework for hypothesis generation, advancing our understanding of the gene regulatory networks driving cell-fate decisions.

Single-cell chromatin accessibility of developing murine pancreas identifies cell state-specific gene regulatory programs

Numerous studies have characterized the existence of cell subtypes, along with their corresponding transcriptional profiles, within the developing mouse pancreas. The upstream mechanisms that initiate and maintain gene expression programs across cell states, however, remain largely unknown. Here, we generate single-nucleus ATAC-Sequencing data of developing murine pancreas and perform an integrated, multi-omic analysis of both chromatin accessibility and RNA expression to describe the chromatin landscape of the developing pancreas at both E14.5 and E17.5 at single-cell resolution. We identify candidate transcription factors regulating cell fate and construct gene regulatory networks of active transcription factor binding to regulatory regions of downstream target genes. This work serves as a valuable resource for the field of pancreatic biology in general and contributes to our understanding of lineage plasticity among endocrine cell types. In addition, these data identify which epigenetic states should be represented in the differentiation of stem cells to the pancreatic beta cell fate to best recapitulate in vitro the gene regulatory networks that are critical for progression along the beta cell lineage in vivo.

Heterogeneity of Islet Cells during Embryogenesis and Differentiation

Diabetes is caused by insufficient insulin secretion due to β-cell dysfunction and/or β-cell loss. Therefore, the restoration of functional β-cells by the induction of β-cell differentiation from embryonic stem (ES) and induced-pluripotent stem (iPS) cells, or from somatic non-β-cells, may be a promising curative therapy. To establish an efficient and feasible method for generating functional insulin-producing cells, comprehensive knowledge of pancreas development and β-cell differentiation, including the mechanisms driving cell fate decisions and endocrine cell maturation is crucial. Recent advances in single-cell RNA sequencing (scRNA-seq) technologies have opened a new era in pancreas development and diabetes research, leading to clarification of the detailed transcriptomes of individual insulin-producing cells. Such extensive high-resolution data enables the inference of developmental trajectories during cell transitions and gene regulatory networks. Additionally, advancements in stem cell research have not only enabled their immediate clinical application, but also has made it possible to observe the genetic dynamics of human cell development and maturation in a dish. In this review, we provide an overview of the heterogeneity of islet cells during embryogenesis and differentiation as demonstrated by scRNA-seq studies on the developing and adult pancreata, with implications for the future application of regenerative medicine for diabetes.

Jag1-Notch cis-interaction determines cell fate segregation in pancreatic development

The Notch ligands Jag1 and Dll1 guide differentiation of multipotent pancreatic progenitor cells (MPCs) into unipotent pro-acinar cells (PACs) and bipotent duct/endocrine progenitors (BPs). Ligand-mediated trans-activation of Notch receptors induces oscillating expression of the transcription factor Hes1, while ligand-receptor cis-interaction indirectly represses Hes1 activation. Despite Dll1 and Jag1 both displaying cis- and trans-interactions, the two mutants have different phenotypes for reasons not fully understood. Here, we present a mathematical model that recapitulates the spatiotemporal differentiation of MPCs into PACs and BPs. The model correctly captures cell fate changes in Notch pathway knockout mice and small molecule inhibitor studies, and a requirement for oscillatory Hes1 expression to maintain the multipotent state. Crucially, the model entails cell-autonomous attenuation of Notch signaling by Jag1-mediated cis-inhibition in MPC differentiation. The model sheds light on the underlying mechanisms, suggesting that cis-interaction is crucial for exiting the multipotent state, while trans-interaction is required for adopting the bipotent fate.

Single-cell transcriptome and accessible chromatin dynamics during endocrine pancreas development

Delineating gene regulatory networks that orchestrate cell-type specification is a continuing challenge for developmental biologists. Single-cell analyses offer opportunities to address these challenges and accelerate discovery of rare cell lineage relationships and mechanisms underlying hierarchical lineage decisions. Here, we describe the molecular analysis of mouse pancreatic endocrine cell differentiation using single-cell transcriptomics, chromatin accessibility assays coupled to genetic labeling, and cytometry-based cell purification. We uncover transcription factor networks that delineate β-, α-, and δ-cell lineages. Through genomic footprint analysis, we identify transcription factor-regulatory DNA interactions governing pancreatic cell development at unprecedented resolution. Our analysis suggests that the transcription factor Neurog3 may act as a pioneer transcription factor to specify the pancreatic endocrine lineage. These findings could improve protocols to generate replacement endocrine cells from renewable sources, like stem cells, for diabetes therapy.