Endocrine lineage biases arise in temporally distinct endocrine progenitors during pancreatic morphogenesis

Endogenous glucocorticoid receptor activation modulates early-stage cell differentiation in pancreatic progenitors of mice and humans

Understanding pancreatic development is instrumental to diabetes research and β-cell replacement therapies. Here, we investigate glucocorticoid receptor (GR) signaling during early pancreas development in mice and humans. Previous reports suggest that glucocorticoids do not play a significant role in mouse pancreas development before the second transition. In this study, we demonstrate that, under physiological conditions, the GR is selectively active in mouse pro-acinar and early endocrine cells from embryonic day 11.5, silenced in bipotent progenitors, and reactivated during endocrine commitment. In mouse pancreatic explants, ectopic GR activation globally promotes acinar fate. Surprisingly, GR activation in human in vitro-derived multipotent pancreatic progenitors steers lineage commitment toward a bipotent/endocrine trajectory and upregulates genes for which expression profiles resemble those of SOX9 and HES1 during human embryonic pancreatic bipotential and endocrine progenitor fate choice. Our combined epigenomic and single-cell transcriptomic analyses suggest that these newly identified marker genes may play important roles in human pancreas development. Taken together, our findings position the GR pathway as an endogenous developmental modulator of early-stage pancreatic progenitor cell differentiation and provide insights into the underlying transcriptional mechanisms involved.

The expression order determines the pioneer functions of NGN3 and NEUROD1 in pancreatic endocrine differentiation

Pioneer transcription factors (TFs) initiate chromatin remodeling, which is crucial for gene regulation and cell differentiation. In this study, we investigated how the sequential expression of neurogenin 3 (NGN3) and NEUROD1 affects their pioneering functions during pancreatic endocrine differentiation. Using a genetically engineered mouse model, we mapped NGN3-binding sites, confirming the pivotal role of this molecule in regulating chromatin accessibility. The pioneering function of NGN3 involves dose tolerance, and low doses are sufficient. Although NEUROD1 generally acts as a conventional TF, it can assume a pioneering role in the absence of NGN3. The sequential expression of NeuroD1 and Ngn3 predominantly drives α cell generation, which may explain the inefficient β cell induction observed in vitro. Our findings demonstrate that pioneer activity is dynamically shaped by temporal TF expression and inter-TF interactions, providing insights into transcriptional regulation and its implications for disease mechanisms and therapeutic targeting and enhancing in vitro differentiation strategies.

ETVs dictate hPSC differentiation by tuning biophysical properties

Stem cells maintain a dynamic dialog with their niche, integrating biochemical and biophysical cues to modulate cellular behavior. Yet, the transcriptional networks that regulate cellular biophysical properties remain poorly defined. Here, we leverage human pluripotent stem cells (hPSCs) and two morphogenesis models - gastruloids and pancreatic differentiation - to establish ETV transcription factors as critical regulators of biophysical parameters and lineage commitment. Genetic ablation of ETV1 or ETV1/ETV4/ETV5 in hPSCs enhances cell-cell and cell-ECM adhesion, leading to aberrant multilineage differentiation including disrupted germ-layer organization, ectoderm loss, and extraembryonic cell overgrowth in gastruloids. Furthermore, ETV1 loss abolishes pancreatic progenitor formation. Single-cell RNA sequencing and follow-up assays reveal dysregulated mechanotransduction via the PI3K/AKT signaling. Our findings highlight the importance of transcriptional control over cell biophysical properties and suggest that manipulating these properties may improve in vitro cell and tissue engineering strategies.

Identification of mouse and human embryonic pancreatic cells with adult Procr+ progenitor transcriptomic and epigenomic characteristics

Background

The quest to find a progenitor cell in the adult pancreas has driven research in the field for decades. Many potential progenitor cell sources have been reported, but so far this is a matter of debate mainly due to reproducibility issues. The existence of adult Procr+ progenitor cells in mice islets has been recently reported. These were shown to comprise ~1% of islet cells, lack expression of Neurog3 and endocrine hormones, and to be capable of differentiating into all endocrine cell types. However, these findings had limited impact, as further evidence supporting the existence and function of Procr+ progenitors has not emerged.

Methods and findings

We report here an unbiased comparison across mouse and human pancreatic samples, including adult islets and embryonic tissue, to track the existence of Procr+ progenitors originally described based on their global gene expression signature. We could not find Procr+ progenitors on other mouse or human adult pancreatic islet samples. Unexpectedly, our results revealed a transcriptionally close mesothelial cell population in the mouse and human embryonic pancreas. These Procr-like mesothelial cells of the embryonic pancreas share the salient transcriptional and epigenomic features of previously reported Procr+ progenitors found in adult pancreatic islets. Notably, we report here that Procr-like transcriptional signature is gradually established in mesothelial cells during mouse pancreas development from E12.5 to E17.5, which has its largest amount. Further supporting a developmentally relevant role in the human pancreas, we additionally report that a transcriptionally similar population is spontaneously differentiated from human pluripotent stem cells cultured in vitro along the pancreatic lineage.

Conclusions

Our results show that, although the previously reported Procr+ progenitor cell population could not be found in other adult pancreatic islet samples, a mesothelial cell population with a closely related transcriptional signature is present in both the mouse and human embryonic pancreas. Several lines of evidence presented in this work support a developmentally relevant function for these Procr-like mesothelial cells.

SPOCK2 controls the proliferation and function of immature pancreatic β-cells through MMP2

Human pluripotent stem cell-derived β-cells (SC-β-cells) represent an alternative cell source for transplantation in diabetic patients. Although mitogens could in theory be used to expand β-cells, adult β-cells very rarely replicate. In contrast, newly formed β-cells, including SC-β-cells, display higher proliferative capacity and distinct transcriptional and functional profiles. Through bidirectional expression modulation and single-cell RNA-seq, we identified SPOCK2, an ECM protein, as an inhibitor of immature β-cell proliferation. Human β-cells lacking SPOCK2 presented elevated MMP2 expression and activity, leading to β-integrin-FAK-c-JUN pathway activation. Treatment with the MMP2 protein resulted in pronounced short- and long-term SC-β-cell expansion, significantly increasing glucose-stimulated insulin secretion in vitro and in vivo. These findings suggest that SPOCK2 mediates fetal β-cell proliferation and maturation. In summary, we identified a molecular mechanism that specifically regulates SC-β-cell proliferation and function, highlighting a unique signaling milieu of SC-β-cells with promise for the robust derivation of fully functional cells for transplantation.

Mapping cells through time and space with moscot

Single-cell genomic technologies enable the multimodal profiling of millions of cells across temporal and spatial dimensions. However, experimental limitations hinder the comprehensive measurement of cells under native temporal dynamics and in their native spatial tissue niche. Optimal transport has emerged as a powerful tool to address these constraints and has facilitated the recovery of the original cellular context1-4. Yet, most optimal transport applications are unable to incorporate multimodal information or scale to single-cell atlases. Here we introduce multi-omics single-cell optimal transport (moscot), a scalable framework for optimal transport in single-cell genomics that supports multimodality across all applications. We demonstrate the capability of moscot to efficiently reconstruct developmental trajectories of 1.7 million cells from mouse embryos across 20 time points. To illustrate the capability of moscot in space, we enrich spatial transcriptomic datasets by mapping multimodal information from single-cell profiles in a mouse liver sample and align multiple coronal sections of the mouse brain. We present moscot.spatiotemporal, an approach that leverages gene-expression data across both spatial and temporal dimensions to uncover the spatiotemporal dynamics of mouse embryogenesis. We also resolve endocrine-lineage relationships of delta and epsilon cells in a previously unpublished mouse, time-resolved pancreas development dataset using paired measurements of gene expression and chromatin accessibility. Our findings are confirmed through experimental validation of NEUROD2 as a regulator of epsilon progenitor cells in a model of human induced pluripotent stem cell islet cell differentiation. Moscot is available as open-source software, accompanied by extensive documentation.

Pancreatic organogenesis mapped through space and time

The spatial organization of cells within a tissue is dictated throughout dynamic developmental processes. We sought to understand whether cells geometrically coordinate with one another throughout development to achieve their organization. The pancreas is a complex cellular organ with a particular spatial organization. Signals from the mesenchyme, neurons, and endothelial cells instruct epithelial cell differentiation during pancreatic development. To understand the cellular diversity and spatial organization of the developing pancreatic niche, we mapped the spatial relationships between single cells over time. We found that four transcriptionally unique subtypes of mesenchyme in the developing pancreas spatially coordinate throughout development, with each subtype at fixed locations in space and time in relation to other cells, including beta cells, vasculature, and epithelial cells. Our work provides insight into the mechanisms of pancreatic development by showing that cells are organized in a space and time manner.

Extracellular matrix proteins refine microenvironments for pancreatic organogenesis from induced pluripotent stem cell differentiation

Rationale: The current understanding on manipulating signaling pathways to generate mature human islet organoids with all major hormone-secreting endocrine cell types (i.e., α, β, δ, and γ cells) from induced pluripotent stem cells (iPSCs) is insufficient. However, donor islet shortage necessitates that we produce functional islets in vitro. In this study, we aimed to find decellularized pancreatic extracellular matrix (dpECM) proteins that leverage signaling pathways and promote functional iPSC islet organogenesis. Methods: We performed proteomic analysis to identify key islet promoting factors from porcine and rat dpECM. With this, we identified collagen type II (COL2) as a potential biomaterial cue that endorses islet development from iPSCs. Using global transcriptome profiling, gene set enrichment analysis, immunofluorescence microscopy, flow cytometry, Western blot, and glucose-stimulated hormonal secretion analysis, we examined COL2's role in regulating iPSC pancreatic lineage specification and signaling pathways, critical to islet organogenesis and morphogenesis. Results: We discovered COL2 acts as a functional biomaterial that augments islet development from iPSCs, similar to collagen type V (COL5) as reported in our earlier study. COL2 substantially stimulates the formation of endocrine progenitors and subsequent islet organoids with significantly elevated expressions of pancreatic signature genes and proteins. Furthermore, it enhances islets' glucose sensitivity for hormonal secretion. A cluster of gene expressions associated with various signaling pathways, including but not limited to oxidative phosphorylation, insulin secretion, cell cycle, the canonical WNT, hypoxia, and interferon-γ response, were considerably affected by COL2 and COL5 cues. Conclusion: We demonstrated dpECM's crucial role in refining stem cell differentiation microenvironments for organoid development and maturation. Our findings on biomaterial-stimulated signaling for stem cell specification, organogenesis, and maturation open up a new way to increase the differentiation efficacy of endocrine tissues that can contribute to the production of biologically functional islets.

The SWI/SNF chromatin remodelling complex regulates pancreatic endocrine cell expansion and differentiation in mice in vivo

AIMS/HYPOTHESIS

Strategies to augment functional beta cell mass include directed differentiation of stem cells towards a beta cell fate, which requires extensive knowledge of transcriptional programs governing endocrine progenitor cell differentiation in vivo. We aimed to study the contributions of the Brahma-related gene-1 (BRG1) and Brahma (BRM) ATPase subunits of the SWI/SNF chromatin remodelling complex to endocrine cell development.

METHODS

We generated mice with endocrine progenitor-specific Neurog3-Cre BRG1 removal in the presence of heterozygous (Brg1Δendo;Brm+/-) or homozygous (double knockout: DKOΔendo) BRM deficiency. Whole-body metabolic phenotyping, islet function characterisation, islet quantitative PCR and histological characterisation were performed on animals and tissues postnatally. To test the mechanistic actions of SWI/SNF in controlling gene expression during endocrine cell development, single-cell RNA-seq was performed on flow-sorted endocrine-committed cells from embryonic day 15.5 control and mutant embryos.

RESULTS

Brg1Δendo;Brm+/- mice exhibit severe glucose intolerance, hyperglycaemia and hypoinsulinaemia, resulting, in part, from reduced islet number; diminished alpha, beta and delta cell mass; compromised islet insulin secretion; and altered islet gene expression programs, including reductions in MAFA and urocortin 3 (UCN3). DKOΔendo mice were not recovered at weaning; however, postnatal day 6 DKOΔendo mice were severely hyperglycaemic with reduced serum insulin levels and beta cell area. Single-cell RNA-seq of embryonic day 15.5 lineage-labelled cells revealed endocrine progenitor, alpha and beta cell populations from SWI/SNF mutants have reduced expression of Mafa, Gcg, Ins1 and Ins2, suggesting limited differentiation capacity. Reduced Neurog3 transcripts were discovered in DKOΔendo endocrine progenitor clusters, and the proliferative capacity of neurogenin 3 (NEUROG3)+ cells was reduced in Brg1Δendo;Brm+/- and DKOΔendo mutants.

CONCLUSIONS/INTERPRETATION

Loss of BRG1 from developing endocrine progenitor cells has a severe postnatal impact on glucose homeostasis, and loss of both subunits impedes animal survival, with both groups exhibiting alterations in hormone transcripts embryonically. Taken together, these data highlight the critical role SWI/SNF plays in governing gene expression programs essential for endocrine cell development and expansion.

DATA AVAILABILITY

Raw and processed data for scRNA-seq have been deposited into the NCBI Gene Expression Omnibus (GEO) database under the accession number GSE248369.

Simulation of CRISPR-Cas9 editing on evolving barcode and accuracy of lineage tracing

We designed a simulation program that mimics the CRISPR-Cas9 editing on evolving barcode and double strand break repair procedure along with cell divisions. Emerging barcode mutations tend to build upon previously existing mutations, occurring sequentially with each generation. This process results in a unique mutation profile in each cell. We sample the barcodes in leaf cells and reconstruct the lineage, comparing it to the original lineage tree to test algorithm accuracy under different parameter settings. Our computational simulations validate the reasonable assumptions deduced from experimental observations, emphasizing that factors such as sampling size, barcode length, multiple barcodes, indel probabilities, and Cas9 activity are critical for accurate and successful lineage tracing. Among the many factors we found that sampling size and indel probabilities are two major ones that affect lineage tracing accuracy. Large segment deletions in early generations could greatly impact lineage accuracy. These simulation results offer insightful recommendations for enhancing the design and analysis of Cas9-mediated molecular barcodes in actual experiments.

scParser: sparse representation learning for scalable single-cell RNA sequencing data analysis

The rapid rise in the availability and scale of scRNA-seq data needs scalable methods for integrative analysis. Though many methods for data integration have been developed, few focus on understanding the heterogeneous effects of biological conditions across different cell populations in integrative analysis. Our proposed scalable approach, scParser, models the heterogeneous effects from biological conditions, which unveils the key mechanisms by which gene expression contributes to phenotypes. Notably, the extended scParser pinpoints biological processes in cell subpopulations that contribute to disease pathogenesis. scParser achieves favorable performance in cell clustering compared to state-of-the-art methods and has a broad and diverse applicability.

Stromal Rigidity Stress Accelerates Pancreatic Intraepithelial Neoplasia Progression and Chromosomal Instability via Nuclear Protein Tyrosine Kinase 2 Localization

Because the mechanotransduction by stromal stiffness stimulates the rupture and repair of the nuclear envelope in pancreatic progenitor cells, accumulated genomic aberrations are under selection in the tumor microenvironment. Analysis of cell growth, micronuclei, and phosphorylated Ser-139 residue of the histone variant H2AX (γH2AX) foci linked to mechanotransduction pressure in vivo during serial orthotopic passages of mouse KrasLSL-G12D/+;Trp53flox/flox;Pdx1-Cre (KPC) cancer cells in the tumor and in migrating through the size-restricted 3-μm micropores. To search for pancreatic cancer cell-of-origin, analysis of single-cell data sets revealed that the extracellular matrix shaped an alternate route of acinar-ductal transdifferentiation of acinar cells into topoisomerase II α (TOP2A)-overexpressing cancer cells and derived subclusters with copy number amplifications in MYC-PTK2 (protein tyrosine kinase 2) locus and PIK3CA. High-PTK2 expression is associated with 171 differentially methylated CpG loci, 319 differentially expressed genes, and poor overall survival in The Cancer Genome Atlas-Pancreatic Adenocarcinoma cohort. Abolished RGD-integrin signaling by disintegrin KG blocked the PTK2 phosphorylation, increased cancer apoptosis, decreased vav guanine nucleotide exchange factor 1 (VAV1) expression, and prolonged overall survival in the KPC mice. Reduction of α-smooth muscle actin deposition in the CD248 knockout KPC mice remodeled the tissue stroma and down-regulated TOP2A expression in the epithelium. In summary, stromal stiffness induced the onset of cancer cells-of-origin by ectopic TOP2A expression, and the genomic amplification of MYC-PTK2 locus via alternative transdifferentiation of pancreatic progenitor cells is the vulnerability useful for disintegrin KG treatment.

Spatiotemporal role of SETD2-H3K36me3 in murine pancreatic organogenesis

Pancreas development is tightly controlled by multilayer mechanisms. Despite years of effort, large gaps remain in understanding how histone modifications coordinate pancreas development. SETD2, a predominant histone methyltransferase of H3K36me3, plays a key role in embryonic stem cell differentiation, whose role in organogenesis remains elusive. Here, by combination of cleavage under targets and tagmentation (CUT&Tag), assay for transposase-accessible chromatin using sequencing (ATAC-seq), and bulk RNA sequencing, we show a dramatic increase in the H3K36me3 level from the secondary transition phase and decipher the related transcriptional alteration. Using single-cell RNA sequencing, we define that pancreatic deletion of Setd2 results in abnormalities in both exocrine and endocrine lineages: hyperproliferative tip progenitor cells lead to abnormal differentiation; Ngn3+ endocrine progenitors decline due to the downregulation of Nkx2.2, leading to insufficient endocrine development. Thus, these data identify SETD2 as a crucial player in embryonic pancreas development, providing a clue to understanding the dysregulation of histone modifications in pancreatic disorders.

Regulation of multiple signaling pathways promotes the consistent expansion of human pancreatic progenitors in defined conditions

The unlimited expansion of human progenitor cells in vitro could unlock many prospects for regenerative medicine. However, it remains an important challenge as it requires the decoupling of the mechanisms supporting progenitor self-renewal and expansion from those mechanisms promoting their differentiation. This study focuses on the expansion of human pluripotent stem (hPS) cell-derived pancreatic progenitors (PP) to advance novel therapies for diabetes. We obtained mechanistic insights into PP expansion requirements and identified conditions for the robust and unlimited expansion of hPS cell-derived PP cells under GMP-compliant conditions through a hypothesis-driven iterative approach. We show that the combined stimulation of specific mitogenic pathways, suppression of retinoic acid signaling, and inhibition of selected branches of the TGFβ and Wnt signaling pathways are necessary for the effective decoupling of PP proliferation from differentiation. This enabled the reproducible, 2000-fold, over 10 passages and 40-45 d, expansion of PDX1+/SOX9+/NKX6-1+ PP cells. Transcriptome analyses confirmed the stabilization of PP identity and the effective suppression of differentiation. Using these conditions, PDX1+/SOX9+/NKX6-1+ PP cells, derived from different, both XY and XX, hPS cell lines, were enriched to nearly 90% homogeneity and expanded with very similar kinetics and efficiency. Furthermore, non-expanded and expanded PP cells, from different hPS cell lines, were differentiated in microwells into homogeneous islet-like clusters (SC-islets) with very similar efficiency. These clusters contained abundant β-cells of comparable functionality as assessed by glucose-stimulated insulin secretion assays. These findings established the signaling requirements to decouple PP proliferation from differentiation and allowed the consistent expansion of hPS cell-derived PP cells. They will enable the establishment of large banks of GMP-produced PP cells derived from diverse hPS cell lines. This approach will streamline SC-islet production for further development of the differentiation process, diabetes research, personalized medicine, and cell therapies.

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.

The way you move

Glial cells in the gut are specialized to fine-tune intestinal function.

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.

Deciphering early human pancreas development at the single-cell level

Understanding pancreas development can provide clues for better treatments of pancreatic diseases. However, the molecular heterogeneity and developmental trajectory of the early human pancreas are poorly explored. Here, we performed large-scale single-cell RNA sequencing and single-cell assay for transposase accessible chromatin sequencing of human embryonic pancreas tissue obtained from first-trimester embryos. We unraveled the molecular heterogeneity, developmental trajectories and regulatory networks of the major cell types. The results reveal that dorsal pancreatic multipotent cells in humans exhibit different gene expression patterns than ventral multipotent cells. Pancreato-biliary progenitors that generate ventral multipotent cells in humans were identified. Notch and MAPK signals from mesenchymal cells regulate the differentiation of multipotent cells into trunk and duct cells. Notably, we identified endocrine progenitor subclusters with different differentiation potentials. Although the developmental trajectories are largely conserved between humans and mice, some distinct gene expression patterns have also been identified. Overall, we provide a comprehensive landscape of early human pancreas development to understand its lineage transitions and molecular complexity.

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.

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.

Single-Cell RNA Sequencing of Sox17-Expressing Lineages Reveals Distinct Gene Regulatory Networks and Dynamic Developmental Trajectories

During early embryogenesis, the transcription factor SOX17 contributes to hepato-pancreato-biliary system formation and vascular-hematopoietic emergence. To better understand Sox17 function in the developing endoderm and endothelium, we developed a dual-color temporal lineage-tracing strategy in mice combined with single-cell RNA sequencing to analyze 6934 cells from Sox17-expressing lineages at embryonic days 9.0-9.5. Our analyses showed 19 distinct cellular clusters combined from all 3 germ layers. Differential gene expression, trajectory and RNA-velocity analyses of endothelial cells revealed a heterogenous population of uncommitted and specialized endothelial subtypes, including 2 hemogenic populations that arise from different origins. Similarly, analyses of posterior foregut endoderm revealed subsets of hepatic, pancreatic, and biliary progenitors with overlapping developmental potency. Calculated gene-regulatory networks predict gene regulons that are dominated by cell type-specific transcription factors unique to each lineage. Vastly different Sox17 regulons found in endoderm versus endothelial cells support the differential interactions of SOX17 with other regulatory factors thereby enabling lineage-specific regulatory actions.

Enteric glial hub cells coordinate intestinal motility

Enteric glia are the predominant cell type in the enteric nervous system yet their identities and roles in gastrointestinal function are not well classified. Using our optimized single nucleus RNA-sequencing method, we identified distinct molecular classes of enteric glia and defined their morphological and spatial diversity. Our findings revealed a functionally specialized biosensor subtype of enteric glia that we call "hub cells." Deletion of the mechanosensory ion channel PIEZO2 from adult enteric glial hub cells, but not other subtypes of enteric glia, led to defects in intestinal motility and gastric emptying in mice. These results provide insight into the multifaceted functions of different enteric glial cell subtypes in gut health and emphasize that therapies targeting enteric glia could advance the treatment of gastrointestinal diseases.

USP7 controls NGN3 stability and pancreatic endocrine lineage development

Understanding the factors and mechanisms involved in beta-cell development will guide therapeutic efforts to generate fully functional beta cells for diabetes. Neurogenin 3 (NGN3) is the key transcription factor that marks endocrine progenitors and drives beta-cell differentiation. Here we screen for binding partners of NGN3 and identify the deubiquitylating enzyme USP7 as a key regulator of NGN3 stability. Mechanistically, USP7 interacts with, deubiquitinates and stabilizes NGN3. In vivo, conditional knockout of Usp7 in the mouse embryonic pancreas causes a dramatic reduction in islet formation and hyperglycemia in adult mice, due to impaired NGN3-mediated endocrine specification during pancreatic development. Furthermore, pharmacological inhibition of USP7 during endocrine specification in human iPSC models of beta-cell differentiation decreases NGN3 expressing progenitor cell numbers and impairs beta cell differentiation. Thus, the USP7-NGN3 axis is an essential mechanism for driving endocrine development and beta-cell differentiation, which can be therapeutically exploited.

Pancreatic islet cell type-specific transcriptomic changes during pregnancy and postpartum

Pancreatic β-cell mass expands during pregnancy and regresses in the postpartum period in conjunction with dynamic metabolic demands on maternal glucose homeostasis. To understand transcriptional changes driving these adaptations in β-cells and other islet cell types, we performed single-cell RNA sequencing on islets from virgin, late gestation, and early postpartum mice. We identified transcriptional signatures unique to gestation and the postpartum in β-cells, including induction of the AP-1 transcription factor subunits and other genes involved in the immediate-early response (IEGs). In addition, we found pregnancy and postpartum-induced changes differed within each endocrine cell type, and in endothelial cells and antigen-presenting cells within islets. Together, our data reveal insights into cell type-specific transcriptional changes responsible for adaptations by islet cells to pregnancy and their resolution postpartum.

Deep into the niche: Deciphering local endoderm-microenvironment interactions in development, homeostasis, and disease of pancreas and intestine

Unraveling molecular and functional heterogeneity of niche cells within the developing endoderm could resolve mechanisms of tissue formation and maturation. Here, we discuss current unknowns in molecular mechanisms underlying key developmental events in pancreatic islet and intestinal epithelial formation. Recent breakthroughs in single-cell and spatial transcriptomics, paralleled with functional studies in vitro, reveal that specialized mesenchymal subtypes drive the formation and maturation of pancreatic endocrine cells and islets via local interactions with epithelium, neurons, and microvessels. Analogous to this, distinct intestinal niche cells regulate both epithelial development and homeostasis throughout life. We propose how this knowledge can be used to progress research in the human context using pluripotent stem cell-derived multilineage organoids. Overall, understanding the interactions between the multitude of microenvironmental cells and how they drive tissue development and function could help us make more therapeutically relevant in vitro models.

ISL1 controls pancreatic alpha cell fate and beta cell maturation

BACKGROUND

Glucose homeostasis is dependent on functional pancreatic α and ß cells. The mechanisms underlying the generation and maturation of these endocrine cells remain unclear.

RESULTS

We unravel the molecular mode of action of ISL1 in controlling α cell fate and the formation of functional ß cells in the pancreas. By combining transgenic mouse models, transcriptomic and epigenomic profiling, we uncover that elimination of Isl1 results in a diabetic phenotype with a complete loss of α cells, disrupted pancreatic islet architecture, downregulation of key ß-cell regulators and maturation markers of ß cells, and an enrichment in an intermediate endocrine progenitor transcriptomic profile.

CONCLUSIONS

Mechanistically, apart from the altered transcriptome of pancreatic endocrine cells, Isl1 elimination results in altered silencing H3K27me3 histone modifications in the promoter regions of genes that are essential for endocrine cell differentiation. Our results thus show that ISL1 transcriptionally and epigenetically controls α cell fate competence, and ß cell maturation, suggesting that ISL1 is a critical component for generating functional α and ß cells.

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.

TrackSOM: Mapping immune response dynamics through clustering of time-course cytometry data

Mapping the dynamics of immune cell populations over time or disease-course is key to understanding immunopathogenesis and devising putative interventions. We present TrackSOM, a novel method for delineating cellular populations and tracking their development over a time- or disease-course cytometry datasets. We demonstrate TrackSOM-enabled elucidation of the immune response to West Nile Virus infection in mice, uncovering heterogeneous subpopulations of immune cells and relating their functional evolution to disease severity. TrackSOM is easy to use, encompasses few parameters, is quick to execute, and enables an integrative and dynamic overview of the immune system kinetics that underlie disease progression and/or resolution.

Lessons from single-cell RNA sequencing of human islets

Islet dysfunction is central in type 2 diabetes and full-blown type 2 diabetes develops first when the beta cells lose their ability to secrete adequate amounts of insulin in response to raised plasma glucose. Several mechanisms behind beta cell dysfunction have been put forward but many important questions still remain. Furthermore, our understanding of the contribution of each islet cell type in type 2 diabetes pathophysiology has been limited by technical boundaries. Closing this knowledge gap will lead to a leap forward in our understanding of the islet as an organ and potentially lead to improved treatments. The development of single-cell RNA sequencing (scRNAseq) has led to a breakthrough for characterising the transcriptome of each islet cell type and several important observations on the regulation of cell-type-specific gene expression have been made. When it comes to identifying type 2 diabetes disease mechanisms, the outcome is still limited. Several studies have identified differentially expressed genes, although there is very limited consensus between the studies. As with all new techniques, scRNAseq has limitations; in addition to being extremely expensive, genes expressed at low levels may not be detected, noise may not be appropriately filtered and selection biases for certain cell types are at hand. Furthermore, recent advances suggest that commonly used computational tools may be suboptimal for analysis of scRNAseq data in small-scale studies. Fortunately, development of new computational tools holds promise for harnessing the full potential of scRNAseq data. Here we summarise how scRNAseq has contributed to increasing the understanding of various aspects of islet biology as well as type 2 diabetes disease mechanisms. We also focus on challenges that remain and propose steps to promote the utilisation of the full potential of scRNAseq in this area.

Identification and implication of tissue-enriched ligands in epithelial-endothelial crosstalk during pancreas development

Development of the pancreas is driven by an intrinsic program coordinated with signals from other cell types in the epithelial environment. These intercellular communications have been so far challenging to study because of the low concentration, localized production and diversity of the signals released. Here, we combined scRNAseq data with a computational interactomic approach to identify signals involved in the reciprocal interactions between the various cell types of the developing pancreas. This in silico approach yielded 40,607 potential ligand-target interactions between the different main pancreatic cell types. Among this vast network of interactions, we focused on three ligands potentially involved in communications between epithelial and endothelial cells. BMP7 and WNT7B, expressed by pancreatic epithelial cells and predicted to target endothelial cells, and SEMA6D, involved in the reverse interaction. In situ hybridization confirmed the localized expression of Bmp7 in the pancreatic epithelial tip cells and of Wnt7b in the trunk cells. On the contrary, Sema6d was enriched in endothelial cells. Functional experiments on ex vivo cultured pancreatic explants indicated that tip cell-produced BMP7 limited development of endothelial cells. This work identified ligands with a restricted tissular and cellular distribution and highlighted the role of BMP7 in the intercellular communications contributing to vessel development and organization during pancreas organogenesis.

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.

Spatial and transcriptional heterogeneity of pancreatic beta cell neogenesis revealed by a time-resolved reporter system

AIMS/HYPOTHESIS

While pancreatic beta cells have been shown to originate from endocrine progenitors in ductal regions, it remains unclear precisely where beta cells emerge from and which transcripts define newborn beta cells. We therefore investigated characteristics of newborn beta cells extracted by a time-resolved reporter system.

METHODS

We established a mouse model, 'Ins1-GFP; Timer', which provides spatial information during beta cell neogenesis with high temporal resolution. Single-cell RNA-sequencing (scRNA-seq) was performed on mouse beta cells sorted by fluorescent reporter to uncover transcriptomic profiles of newborn beta cells. scRNA-seq of human embryonic stem cell (hESC)-derived beta-like cells was also performed to compare newborn beta cell features between mouse and human.

RESULTS

Fluorescence imaging of Ins1-GFP; Timer mouse pancreas successfully dissected newly generated beta cells as green fluorescence-dominant cells. This reporter system revealed that, as expected, some newborn beta cells arise close to the ducts (β^duct); unexpectedly, the others arise away from the ducts and adjacent to blood vessels (β^vessel). Single-cell transcriptomic analyses demonstrated five distinct populations among newborn beta cells, confirming spatial heterogeneity of beta cell neogenesis such as high probability of glucagon-positive β^duct, musculoaponeurotic fibrosarcoma oncogene family B (MafB)-positive β^duct and musculoaponeurotic fibrosarcoma oncogene family A (MafA)-positive β^vessel cells. Comparative analysis with scRNA-seq data of mouse newborn beta cells and hESC-derived beta-like cells uncovered transcriptional similarity between mouse and human beta cell neogenesis including microsomal glutathione S-transferase 1 (MGST1)- and synaptotagmin 13 (SYT13)-highly-expressing state.

CONCLUSIONS/INTERPRETATION

The combination of time-resolved histological imaging with single-cell transcriptional mapping demonstrated novel features of spatial and transcriptional heterogeneity in beta cell neogenesis, which will lead to a better understanding of beta cell differentiation for future cell therapy.

DATA AVAILABILITY

Raw and processed single-cell RNA-sequencing data for this study has been deposited in the Gene Expression Omnibus under accession number GSE155742.

Human pancreatic microenvironment promotes β-cell differentiation via non-canonical WNT5A/JNK and BMP signaling

In vitro derivation of pancreatic β-cells from human pluripotent stem cells holds promise as diabetes treatment. Despite recent progress, efforts to generate physiologically competent β-cells are still hindered by incomplete understanding of the microenvironment's role in β-cell development and maturation. Here, we analyze the human mesenchymal and endothelial primary cells from weeks 9-20 fetal pancreas and identify a time point-specific microenvironment that permits β-cell differentiation. Further, we uncover unique factors that guide in vitro development of endocrine progenitors, with WNT5A markedly improving human β-cell differentiation. WNT5A initially acts through the non-canonical (JNK/c-JUN) WNT signaling and cooperates with Gremlin1 to inhibit the BMP pathway during β-cell maturation. Interestingly, we also identify the endothelial-derived Endocan as a SST+ cell promoting factor. Overall, our study shows that the pancreatic microenvironment-derived factors can mimic in vivo conditions in an in vitro system to generate bona fide β-cells for translational applications.

The progress of pluripotent stem cell-derived pancreatic β-cells regeneration for diabetic therapy

Diabetes is a complex metabolic disorder of carbohydrate metabolism, characterized by high blood glucose levels either due to an absolute deficiency of insulin secretion or an ineffective response of cells to insulin, a hormone synthetized by β-cells in the pancreas. Despite the current substantial progress of new drugs and strategies to prevent and treat diabetes, we do not understand precisely the exact cause of the failure and impairment of β-cells. Therefore, there is an urgent need to find new methods to restore β-cells. In recent years, pluripotent stem cells (PSCs) such as embryonic stem cells (ESCs) and induced pluripotent stem cells (iPSC) can serve as an ideal alternative source for the pancreatic β-cells. In this review, we systematically summarize the current progress and protocols of generating pancreatic β-cells from human PSCs. Meanwhile, we also discuss some challenges and future perspectives of human PSCs treatments for diabetes.

Recent advances in tissue stem cells

Stem cells are undifferentiated cells capable of self-renewal and differentiation, giving rise to specialized functional cells. Stem cells are of pivotal importance for organ and tissue development, homeostasis, and injury and disease repair. Tissue-specific stem cells are a rare population residing in specific tissues and present powerful potential for regeneration when required. They are usually named based on the resident tissue, such as hematopoietic stem cells and germline stem cells. This review discusses the recent advances in stem cells of various tissues, including neural stem cells, muscle stem cells, liver progenitors, pancreatic islet stem/progenitor cells, intestinal stem cells, and prostate stem cells, and the future perspectives for tissue stem cell research.

Restoring normal islet mass and function in type 1 diabetes through regenerative medicine and tissue engineering

Type 1 diabetes is characterised by autoimmune-mediated destruction of pancreatic β-cell mass. With the advent of insulin therapy a century ago, type 1 diabetes changed from a progressive, fatal disease to one that requires lifelong complex self-management. Replacing the lost β-cell mass through transplantation has proven successful, but limited donor supply and need for lifelong immunosuppression restricts widespread use. In this Review, we highlight incremental advances over the past 20 years and remaining challenges in regenerative medicine approaches to restoring β-cell mass and function in type 1 diabetes. We begin by summarising the role of endocrine islets in glucose homoeostasis and how this is altered in disease. We then discuss the potential regenerative capacity of the remaining islet cells and the utility of stem cell-derived β-like cells to restore β-cell function. We conclude with tissue engineering approaches that might improve the engraftment, function, and survival of β-cell replacement therapies.

Monogenic Diabetes Modeling: In Vitro Pancreatic Differentiation From Human Pluripotent Stem Cells Gains Momentum

The occurrence of diabetes mellitus is characterized by pancreatic β cell loss and chronic hyperglycemia. While Type 1 and Type 2 diabetes are the most common types, rarer forms involve mutations affecting a single gene. This characteristic has made monogenic diabetes an interesting disease group to model in vitro using human pluripotent stem cells (hPSCs). By altering the genotype of the original hPSCs or by deriving human induced pluripotent stem cells (hiPSCs) from patients with monogenic diabetes, changes in the outcome of the in vitro differentiation protocol can be analyzed in detail to infer the regulatory mechanisms affected by the disease-associated genes. This approach has been so far applied to a diversity of genes/diseases and uncovered new mechanisms. The focus of the present review is to discuss the latest findings obtained by modeling monogenic diabetes using hPSC-derived pancreatic cells generated in vitro. We will specifically focus on the interpretation of these studies, the advantages and limitations of the models used, and the future perspectives for improvement.

Pancreatic cell fate specification: insights into developmental mechanisms and their application for lineage reprogramming

Diabetes is a group of metabolic disorders, which results from insufficient functional pancreatic β-cell mass either due to the autoimmune destruction of insulin producing β-cells, or their death or de-differentiation as compensation for insulin resistance. The ability to reprogram cell types within close developmental proximity to β-cells offers a strategy to replenish β-cell mass and a future possible treatment of diabetes. Here, we review recent advances in the fields of pancreas development and lineage reprogramming. We also probe the possibility of using reprogrammed cells as an approach by which to further understand developmental mechanisms, in particular roadblocks to changing cell identity. Finally, we highlight fundamental challenges that need to be overcome to advance lineage reprogramming for generating pancreatic cells.

Endocrine Pancreas Development and Dysfunction Through the Lens of Single-Cell RNA-Sequencing

A chronic inability to maintain blood glucose homeostasis leads to diabetes, which can damage multiple organs. The pancreatic islets regulate blood glucose levels through the coordinated action of islet cell-secreted hormones, with the insulin released by β-cells playing a crucial role in this process. Diabetes is caused by insufficient insulin secretion due to β-cell loss, or a pancreatic dysfunction. The restoration of a functional β-cell mass might, therefore, offer a cure. To this end, major efforts are underway to generate human β-cells de novo, in vitro, or in vivo. The efficient generation of functional β-cells requires a comprehensive knowledge of pancreas development, including the mechanisms driving cell fate decisions or endocrine cell maturation. Rapid progress in single-cell RNA sequencing (scRNA-Seq) technologies has brought a new dimension to pancreas development research. These methods can capture the transcriptomes of thousands of individual cells, including rare cell types, subtypes, and transient states. With such massive datasets, it is possible to infer the developmental trajectories of cell transitions and gene regulatory pathways. Here, we summarize recent advances in our understanding of endocrine pancreas development and function from scRNA-Seq studies on developing and adult pancreas and human endocrine differentiation models. We also discuss recent scRNA-Seq findings for the pathological pancreas in diabetes, and their implications for better treatment.

Maternal obesity as a risk factor for developing diabetes in offspring: An epigenetic point of view

According to the developmental origin of health and disease concept, the risk of many age-related diseases is not only determined by genetic and adult lifestyle factors but also by factors acting during early development. In particular, maternal obesity and neonatal accelerated growth predispose offspring to overweight and type 2 diabetes (T2D) in adulthood. This concept mainly relies on the developmental plasticity of adipose tissue and pancreatic β-cell programming in response to suboptimal milieu during the perinatal period. These changes result in unhealthy hypertrophic adipocytes with decreased capacity to store fat, low-grade inflammation and loss of insulin-producing pancreatic β-cells. Over the past years, many efforts have been made to understand how maternal obesity induces long-lasting adipose tissue and pancreatic β-cell dysfunction in offspring and what are the molecular basis of the transgenerational inheritance of T2D. In particular, rodent studies have shed light on the role of epigenetic mechanisms in linking maternal nutritional manipulations to the risk for T2D in adulthood. In this review, we discuss epigenetic adipocyte and β-cell remodeling during development in the progeny of obese mothers and the persistence of these marks as a basis of obesity and T2D predisposition.

Functional maturation of immature β cells: A roadblock for stem cell therapy for type 1 diabetes

Type 1 diabetes mellitus (T1DM) is a chronic autoimmune disease caused by the specific destruction of pancreatic islet β cells and is characterized as the absolute insufficiency of insulin secretion. Current insulin replacement therapy supplies insulin in a non-physiological way and is associated with devastating complications. Experimental islet transplantation therapy has been proven to restore glucose homeostasis in people with severe T1DM. However, it is restricted by many factors such as severe shortage of donor sources, progressive loss of donor cells, high cost, etc. As pluripotent stem cells have the potential to give rise to all cells including islet β cells in the body, stem cell therapy for diabetes has attracted great attention in the academic community and the general public. Transplantation of islet β-like cells differentiated from human pluripotent stem cells (hPSCs) has the potential to be an excellent alternative to islet transplantation. In stem cell therapy, obtaining β cells with complete insulin secretion in vitro is crucial. However, after much research, it has been found that the β-like cells obtained by in vitro differentiation still have many defects, including lack of adult-type glucose stimulated insulin secretion, and multi-hormonal secretion, suggesting that in vitro culture does not allows for obtaining fully mature β-like cells for transplantation. A large number of studies have found that many transcription factors play important roles in the process of transforming immature to mature human islet β cells. Furthermore, PDX1, NKX6.1, SOX9, NGN3, PAX4, etc., are important in inducing hPSC differentiation in vitro. The absent or deficient expression of any of these key factors may lead to the islet development defect in vivo and the failure of stem cells to differentiate into genuine functional β-like cells in vitro. This article reviews β cell maturation in vivo and in vitro and the vital roles of key molecules in this process, in order to explore the current problems in stem cell therapy for diabetes.

Functional maturation of immature β cells: A roadblock for stem cell therapy for type 1 diabetes

Type 1 diabetes mellitus (T1DM) is a chronic autoimmune disease caused by the specific destruction of pancreatic islet β cells and is characterized as the absolute insufficiency of insulin secretion. Current insulin replacement therapy supplies insulin in a non-physiological way and is associated with devastating complications. Experimental islet transplantation therapy has been proven to restore glucose homeostasis in people with severe T1DM. However, it is restricted by many factors such as severe shortage of donor sources, progressive loss of donor cells, high cost, etc. As pluripotent stem cells have the potential to give rise to all cells including islet β cells in the body, stem cell therapy for diabetes has attracted great attention in the academic community and the general public. Transplantation of islet β-like cells differentiated from human pluripotent stem cells (hPSCs) has the potential to be an excellent alternative to islet transplantation. In stem cell therapy, obtaining β cells with complete insulin secretion in vitro is crucial. However, after much research, it has been found that the β-like cells obtained by in vitro differentiation still have many defects, including lack of adult-type glucose stimulated insulin secretion, and multi-hormonal secretion, suggesting that in vitro culture does not allows for obtaining fully mature β-like cells for transplantation. A large number of studies have found that many transcription factors play important roles in the process of transforming immature to mature human islet β cells. Furthermore, PDX1, NKX6.1, SOX9, NGN3, PAX4, etc., are important in inducing hPSC differentiation in vitro. The absent or deficient expression of any of these key factors may lead to the islet development defect in vivo and the failure of stem cells to differentiate into genuine functional β-like cells in vitro. This article reviews β cell maturation in vivo and in vitro and the vital roles of key molecules in this process, in order to explore the current problems in stem cell therapy for diabetes.

Groucho co-repressor proteins regulate β cell development and proliferation by repressing Foxa1 in the developing mouse pancreas

Groucho-related genes (GRGs) are transcriptional co-repressors that are crucial for many developmental processes. Several essential pancreatic transcription factors are capable of interacting with GRGs; however, the in vivo role of GRG-mediated transcriptional repression in pancreas development is still not well understood. In this study, we used complex mouse genetics and transcriptomic analyses to determine that GRG3 is essential for β cell development, and in the absence of Grg3 there is compensatory upregulation of Grg4 Grg3/4 double mutant mice have severe dysregulation of the pancreas gene program with ectopic expression of canonical liver genes and Foxa1, a master regulator of the liver program. Neurod1, an essential β cell transcription factor and predicted target of Foxa1, becomes downregulated in Grg3/4 mutants, resulting in reduced β cell proliferation, hyperglycemia, and early lethality. These findings uncover novel functions of GRG-mediated repression during pancreas development.

Sequential progenitor states mark the generation of pancreatic endocrine lineages in mice and humans

The pancreatic islet contains multiple hormone⁺ endocrine lineages (α, β, δ, PP and ε cells), but the developmental processes that underlie endocrinogenesis are poorly understood. Here, we generated novel mouse lines and combined them with various genetic tools to enrich all types of hormone⁺ cells for well-based deep single-cell RNA sequencing (scRNA-seq), and gene coexpression networks were extracted from the generated data for the optimization of high-throughput droplet-based scRNA-seq analyses. These analyses defined an entire endocrinogenesis pathway in which different states of endocrine progenitor (EP) cells sequentially differentiate into specific endocrine lineages in mice. Subpopulations of the EP cells at the final stage (EP4^early and EP4^late) show different potentials for distinct endocrine lineages. ε cells and an intermediate cell population were identified as distinct progenitors that independently generate both α and PP cells. Single-cell analyses were also performed to delineate the human pancreatic endocrinogenesis process. Although the developmental trajectory of pancreatic lineages is generally conserved between humans and mice, clear interspecies differences, including differences in the proportions of cell types and the regulatory networks associated with the differentiation of specific lineages, have been detected. Our findings support a model in which sequential transient progenitor cell states determine the differentiation of multiple cell lineages and provide a blueprint for directing the generation of pancreatic islets in vitro.

Insights from single cell studies of human pancreatic islets and stem cell-derived islet cells to guide functional beta cell maturation in vitro

There is now a sizeable number of single cell transcriptomics studies performed on human and rodent pancreatic islets that have shed light on the unique gene signatures and level of heterogeneity within each individual islet cell type. Following closely from these studies, there is also rapidly-growing activity on characterizing islet-like cells derived from in vitro differentiation of human pluripotent stem cells (hPSCs) at the single cell level. The overall consensus across the studies so far suggests that the first few stages of differentiation are largely uniform, whereas during pancreatic endocrine commitment, cell trajectories start to diverge, resulting in multiple end-stage pancreatic cells that include progenitor-like, endocrine and non-endocrine cells. Comprehensive transcriptional profiling is important for understanding how and why islet cells, especially the insulin-secreting beta cells, exist in subpopulations that differ in maturity, proliferation rate, sensitivity to stress, and insulin secretion function. For hPSC-derived beta cells to be used confidently for cell therapy, optimal differentiation and thorough characterization is required. The key questions to address are-What is the trajectory of differentiation? Is heterogeneity a natural occurrence or is it a consequence of imperfect differentiation protocols? Can lessons be drawn from the extensive single cell transcriptomic data to help guide maturation of beta cells in vitro? This book chapter seeks to address some of these questions, and facilitate ongoing efforts in improving the beta cell differentiation pipeline or enriching for desired beta cell populations following differentiation, to make way for better mechanistic studies and future clinical translation.

Molecular mechanisms of transcription factor mediated cell reprogramming: conversion of liver to pancreas

Transdifferentiation is a type of cellular reprogramming involving the conversion of one differentiated cell type to another. This remarkable phenomenon holds enormous promise for the field of regenerative medicine. Over the last 20 years techniques used to reprogram cells to alternative identities have advanced dramatically. Cellular identity is determined by the transcriptional profile which comprises the subset of mRNAs, and therefore proteins, being expressed by a cell at a given point in time. A better understanding of the levers governing transcription factor activity benefits our ability to generate therapeutic cell types at will. One well-established example of transdifferentiation is the conversion of hepatocytes to pancreatic β-cells. This cell type conversion potentially represents a novel therapy in T1D treatment. The identification of key master regulator transcription factors (which distinguish one body part from another) during embryonic development has been central in developing transdifferentiation protocols. Pdx1 is one such example of a master regulator. Ectopic expression of vector-delivered transcription factors (particularly the triumvirate of Pdx1, Ngn3 and MafA) induces reprogramming through broad transcriptional remodelling. Increasingly, complimentary cell culture techniques, which recapitulate the developmental microenvironment, are employed to coax cells to adopt new identities by indirectly regulating transcription factor activity via intracellular signalling pathways. Both transcription factor-based reprogramming and directed differentiation approaches ultimately exploit transcription factors to influence cellular identity. Here, we explore the evolution of reprogramming and directed differentiation approaches within the context of hepatocyte to β-cell transdifferentiation focussing on how the introduction of new techniques has improved our ability to generate β-cells.

Gene Signatures of NEUROGENIN3+ Endocrine Progenitor Cells in the Human Pancreas

NEUROGENIN3+ (NEUROG3+) cells are considered to be pancreatic endocrine progenitors. Our current knowledge on the molecular program of NEUROG3+ cells in humans is largely extrapolated from studies in mice. We hypothesized that single-cell RNA-seq enables in-depth exploration of the rare NEUROG3+ cells directly in humans. We aligned four large single-cell RNA-seq datasets from postnatal human pancreas. Our integrated analysis revealed 10 NEUROG3+ epithelial cells from a total of 11,174 pancreatic cells. Noticeably, human NEUROG3+ cells clustered with mature pancreatic cells and epsilon cells displayed the highest frequency of NEUROG3 positivity. We confirmed the co-expression of NEUROG3 with endocrine markers and the high percentage of NEUROG3+ cells among epsilon cells at the protein level based on immunostaining on pancreatic tissue sections. We further identified unique genetic signatures of the NEUROG3+ cells. Regulatory network inference revealed novel transcription factors including Prospero homeobox protein 1 (PROX1) may act jointly with NEUROG3. As NEUROG3 plays a central role in endocrine differentiation, knowledge gained from our study will accelerate the development of beta cell regeneration therapies to treat diabetes.

Single-cell lineage analysis reveals extensive multimodal transcriptional control during directed beta-cell differentiation

The in vitro differentiation of insulin-producing beta-like cells can model aspects of human pancreatic development. Here, we generate 95,308 single-cell transcriptomes and reconstruct a lineage tree of the entire differentiation process from human embryonic stem cells to beta-like cells to study temporally regulated genes during differentiation. We identify so-called 'switch genes' at the branch point of endocrine/non-endocrine cell fate choice, revealing insights into the mechanisms of differentiation-promoting reagents, such as NOTCH and ROCKII inhibitors, and providing improved differentiation protocols. Over 20% of all detectable genes are activated multiple times during differentiation, even though their enhancer activation is usually unimodal, indicating extensive gene reuse driven by different enhancers. We also identify a stage-specific enhancer at the TCF7L2 locus for diabetes, uncovered by genome-wide association studies, that drives a transient wave of gene expression in pancreatic progenitors. Finally, we develop a web app to visualize gene expression on the lineage tree, providing a comprehensive single-cell data resource for researchers studying islet biology and diabetes.

Integration of single-cell datasets reveals novel transcriptomic signatures of β-cells in human type 2 diabetes

Pancreatic islet β-cell failure is key to the onset and progression of type 2 diabetes (T2D). The advent of single-cell RNA sequencing (scRNA-seq) has opened the possibility to determine transcriptional signatures specifically relevant for T2D at the β-cell level. Yet, applications of this technique have been underwhelming, as three independent studies failed to show shared differentially expressed genes in T2D β-cells. We performed an integrative analysis of the available datasets from these studies to overcome confounding sources of variability and better highlight common T2D β-cell transcriptomic signatures. After removing low-quality transcriptomes, we retained 3046 single cells expressing 27 931 genes. Cells were integrated to attenuate dataset-specific biases, and clustered into cell type groups. In T2D β-cells (n = 801), we found 210 upregulated and 16 downregulated genes, identifying key pathways for T2D pathogenesis, including defective insulin secretion, SREBP signaling and oxidative stress. We also compared these results with previous data of human T2D β-cells from laser capture microdissection and diabetic rat islets, revealing shared β-cell genes. Overall, the present study encourages the pursuit of single β-cell RNA-seq analysis, preventing presently identified sources of variability, to identify transcriptomic changes associated with human T2D and underscores specific traits of dysfunctional β-cells across different models and techniques.