A Single-Cell Atlas of the Murine Pancreatic Ductal Tree Identifies Novel Cell Populations With Potential Implications in Pancreas Regeneration and Exocrine Pathogenesis.

BACKGROUND & AIMS: Pancreatic ducts form an intricate network of tubules that secrete bicarbonate and drive acinar secretions into the duodenum. This network is formed by centroacinar cells, terminal, intercalated, intracalated ducts, and the main pancreatic duct. Ductal heterogeneity at the single-cell level has been poorly characterized; therefore, our understanding of the role of ductal cells in pancreas regeneration and exocrine pathogenesis has been hampered by the limited knowledge and unexplained diversity within the ductal network. METHODS: We used single cell RNA sequencing to comprehensively characterize mouse ductal heterogeneity at single-cell resolution of the entire ductal epithelium from centroacinar cells to the main duct. Moreover, we used organoid cultures, injury models, and pancreatic tumor samples to interrogate the role of novel ductal populations in pancreas regeneration and exocrine pathogenesis. RESULTS: We have identified the coexistence of 15 ductal populations within the healthy pancreas and characterized their organoid formation capacity and endocrine differentiation potential. Cluster isolation and subsequent culturing let us identify ductal cell populations with high organoid formation capacity and endocrine and exocrine differentiation potential in vitro, including a Wnt-responsive population, a ciliated population, and Flrt3+ cells. Moreover, we have characterized the location of these novel ductal populations in healthy pancreas, chronic pancreatitis, and tumor samples. The expression of Wnt-responsive, interferon-responsive, and epithelial-to-mesenchymal transition population markers increases in chronic pancreatitis and tumor samples. CONCLUSIONS: In light of our discovery of previously unidentified ductal populations, we unmask potential roles of specific ductal populations in pancreas regeneration and exocrine pathogenesis. Thus, novel lineage-tracing models are needed to investigate ductal-specific populations in vivo.

NKX2-2 based nuclei sorting on frozen human archival pancreas enables the enrichment of islet endocrine populations for single-nucleus RNA sequencing.

BACKGROUND: Current approaches to profile the single-cell transcriptomics of human pancreatic endocrine cells almost exclusively rely on freshly isolated islets. However, human islets are limited in availability. Furthermore, the extensive processing steps during islet isolation and subsequent single cell dissolution might alter gene expressions. In this work, we report the development of a single-nucleus RNA sequencing (snRNA-seq) approach with targeted islet cell enrichment for endocrine-population focused transcriptomic profiling using frozen archival pancreatic tissues without islet isolation. RESULTS: We cross-compared five nuclei isolation protocols and selected the citric acid method as the best strategy to isolate nuclei with high RNA integrity and low cytoplasmic contamination from frozen archival human pancreata. We innovated fluorescence-activated nuclei sorting based on the positive signal of NKX2-2 antibody to enrich nuclei of the endocrine population from the entire nuclei pool of the pancreas. Our sample preparation procedure generated high-quality single-nucleus gene-expression libraries while preserving the endocrine population diversity. In comparison with single-cell RNA sequencing (scRNA-seq) library generated with live cells from freshly isolated human islets, the snRNA-seq library displayed comparable endocrine cellular composition and cell type signature gene expression. However, between these two types of libraries, differential enrichments of transcripts belonging to different functional classes could be observed. CONCLUSIONS: Our work fills a technological gap and helps to unleash frozen archival pancreatic tissues for molecular profiling targeting the endocrine population. This study opens doors to retrospective mappings of endocrine cell dynamics in pancreatic tissues of complex histopathology. We expect that our protocol is applicable to enrich nuclei for transcriptomics studies from various populations in different types of frozen archival tissues.

Author Correction: Subepithelial telocytes are an important source of Wnts that supports intestinal crypts.

Change history: In this Letter, the surname of author Efi E. Massasa was misspelled 'Massassa'. This error has been corrected online.

Subepithelial telocytes are an important source of Wnts that supports intestinal crypts.

Tissues that undergo rapid cellular turnover, such as the mammalian haematopoietic system or the intestinal epithelium, are dependent on stem and progenitor cells that proliferate to provide differentiated cells to maintain organismal health. Stem and progenitor cells, in turn, are thought to rely on signals and growth factors provided by local niche cells to support their function and self-renewal. Several cell types have been hypothesized to provide the signals required for the proliferation and differentiation of the intestinal stem cells in intestinal crypts1-6. Here we identify subepithelial telocytes as an important source of Wnt proteins, without which intestinal stem cells cannot proliferate and support epithelial renewal. Telocytes are large but rare mesenchymal cells that are marked by expression of FOXL1 and form a subepithelial plexus that extends from the stomach to the colon. While supporting the entire epithelium, FOXL1+ telocytes compartmentalize the production of Wnt ligands and inhibitors to enable localized pathway activation. Conditional genetic ablation of porcupine (Porcn), which is required for functional maturation of all Wnt proteins, in mouse FOXL1+ telocytes causes rapid cessation of Wnt signalling to intestinal crypts, followed by loss of proliferation of stem and transit amplifying cells and impaired epithelial renewal. Thus, FOXL1+ telocytes are an important source of niche signals to intestinal stem cells.

A negative reciprocal regulatory axis between cyclin D1 and HNF4α modulates cell cycle progression and metabolism in the liver.

Hepatocyte nuclear factor 4α (HNF4α) is a master regulator of liver function and a tumor suppressor in hepatocellular carcinoma (HCC). In this study, we explore the reciprocal negative regulation of HNF4α and cyclin D1, a key cell cycle protein in the liver. Transcriptomic analysis of cultured hepatocyte and HCC cells found that cyclin D1 knockdown induced the expression of a large network of HNF4α-regulated genes. Chromatin immunoprecipitation-sequencing (ChIP-seq) demonstrated that cyclin D1 inhibits the binding of HNF4α to thousands of targets in the liver, thereby diminishing the expression of associated genes that regulate diverse metabolic activities. Conversely, acute HNF4α deletion in the liver induces cyclin D1 and hepatocyte cell cycle progression; concurrent cyclin D1 ablation blocked this proliferation, suggesting that HNF4α maintains proliferative quiescence in the liver, at least, in part, via repression of cyclin D1. Acute cyclin D1 deletion in the regenerating liver markedly inhibited hepatocyte proliferation after partial hepatectomy, confirming its pivotal role in cell cycle progression in this in vivo model, and enhanced the expression of HNF4α target proteins. Hepatocyte cyclin D1 gene ablation caused markedly increased postprandial liver glycogen levels (in a HNF4α-dependent fashion), indicating that the cyclin D1-HNF4α axis regulates glucose metabolism in response to feeding. In AML12 hepatocytes, cyclin D1 depletion led to increased glucose uptake, which was negated if HNF4α was depleted simultaneously, and markedly elevated glycogen synthesis. To summarize, mutual repression by cyclin D1 and HNF4α coordinately controls the cell cycle machinery and metabolism in the liver.

A High-Content Screen Identifies MicroRNAs That Regulate Liver Repopulation After Injury in Mice.

BACKGROUND & AIMS: Liver regeneration is impaired in mice with hepatocyte-specific deficiencies in microRNA (miRNA) processing, but it is not clear which miRNAs regulate this process. We developed a high-throughput screen to identify miRNAs that regulate hepatocyte repopulation after toxic liver injury using fumarylacetoacetate hydrolase-deficient mice. METHODS: We constructed plasmid pools encoding more than 30,000 tough decoy miRNA inhibitors (hairpin nucleic acids designed to specifically inhibit interactions between miRNAs and their targets) to target hepatocyte miRNAs in a pairwise manner. The plasmid libraries were delivered to hepatocytes in fumarylacetoacetate hydrolase-deficient mice at the time of liver injury via hydrodynamic tail-vein injection. Integrated transgene-containing transposons were quantified after liver repopulation via high-throughput sequencing. Changes in polysome-bound transcripts after miRNA inhibition were determined using translating ribosome affinity purification followed by high-throughput sequencing. RESULTS: Analyses of tough decoy abundance in hepatocyte genomic DNA and input plasmid pools identified several thousand miRNA inhibitors that were significantly depleted or increased after repopulation. We classified a subset of miRNA binding sites as those that have strong effects on liver repopulation, implicating the targeted hepatocyte miRNAs as regulators of this process. We then generated a high-content map of pairwise interactions between 171 miRNA-binding sites and identified synergistic and redundant effects. CONCLUSIONS: We developed a screen to identify miRNAs that regulate liver repopulation after injury in live mice.

Combinatorial genetics in liver repopulation and carcinogenesis with a in vivo CRISPR activation platform.

Clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated 9 activation (CRISPRa) systems have enabled genetic screens in cultured cell lines to discover and characterize drivers and inhibitors of cancer cell growth. We adapted this system for use in vivo to assess whether modulating endogenous gene expression levels can result in functional outcomes in the native environment of the liver. We engineered the catalytically dead CRISPR-associated 9 (dCas9)-positive mouse, cyclization recombination-inducible (Cre) CRISPRa system for cell type-specific gene activation in vivo. We tested the capacity for genetic screening in live animals by applying CRISPRa in a clinically relevant model of liver injury and repopulation. We targeted promoters of interest in regenerating hepatocytes using multiple single guide RNAs (gRNAs), and employed high-throughput sequencing to assess enrichment of gRNA sequences during liver repopulation and to link specific gRNAs to the initiation of carcinogenesis. All components of the CRISPRa system were expressed in a cell type-specific manner and activated endogenous gene expression in vivo. Multiple gRNA cassettes targeting a proto-oncogene were significantly enriched following liver repopulation, indicative of enhanced division of cells expressing the proto-oncogene. Furthermore, hepatocellular carcinomas developed containing gRNAs that activated this oncogene, indicative of cancer initiation events. Also, we employed our system for combinatorial cancer genetics in vivo as we found that while clonal hepatocellular carcinomas were dependent on the presence of the oncogene-inducing gRNAs, they were depleted for multiple gRNAs activating tumor suppressors. CONCLUSION: The in vivo CRISPRa platform developed here allows for parallel and combinatorial genetic screens in live animals; this approach enables screening for drivers and suppressors of cell replication and tumor initiation. (Hepatology 2017).

Pericentral hepatocytes produce insulin-like growth factor-2 to promote liver regeneration during selected injuries in mice.

Liver regeneration (LR) happens after various types of injuries. Unlike the well-studied LR caused by partial hepatectomy (PHx), there is accumulating evidence suggesting that LR during other injuries may result from unknown mechanisms. In this study, we found that insulin-like growth factor 2 (IGF-2) was drastically induced following the liver injuries caused by tyrosinemia or long-term treatments of CCl4 . However, this was not observed during the early phase of acute liver injuries after PHx or single treatment of CCl4 . Remarkably, most IGF-2-expressing hepatocytes were located at the histological area around the central vein of the liver lobule after the liver injuries caused either in fumarylacetoacetate hydrolase-deficient mice or in CCl4 chronically treated mice. Hepatocyte proliferation in vivo was significantly promoted by induced IGF-2 overexpression, which could be inhibited by adeno-associated virus-delivered IGF-2 short hairpin RNAs or linsitinib, an inhibitor of IGF-2 signaling. Proliferating hepatocytes in vivo responded to IGF-2 through both insulin receptor and IGF-1 receptor. IGF-2 also significantly promoted DNA synthesis of primary hepatocytes in vitro. More interestingly, the significantly induced IGF-2 was also found to colocalize with glutamine synthetase in the region enriched with proliferating hepatocytes for the liver samples from patients with liver fibrosis. CONCLUSION: IGF-2 is produced by pericentral hepatocytes to promote hepatocyte proliferation and repair tissue damage in the setting of chronic liver injury, which is distinct from the signaling that occurs post-PHx. (Hepatology 2017;66:2002-2015).

Genetic lineage tracing analysis of the cell of origin of hepatotoxin-induced liver tumors in mice.

UNLABELLED: The expression of biliary/progenitor markers by hepatocellular carcinoma (HCC) is often associated with poor prognosis and stem cell-like behaviors of tumor cells. Hepatocellular adenomas (HCAs) also often express biliary/progenitor markers and frequently act as precursor lesions for HCC. However, the cell of origin of HCA and HCC that expresses these markers remains unclear. Therefore, to evaluate if mature hepatocytes give rise to HCA and HCC tumors and to understand the molecular pathways involved in tumorigenesis, we lineage-labeled hepatocytes by injecting adeno-associated virus containing thyroxine-binding globulin promoter-driven causes recombination (AAV-TBG-Cre) into Rosa(YFP) mice. Yellow fluorescent protein (YFP) was present in >96% of hepatocytes before exposure to carcinogens. We treated AAV-TBG-Cre; Rosa(YFP) mice with diethylnitrosamine (DEN), followed by multiple injections of carbon tetrachloride to induce carcinogenesis and fibrosis and found that HCA and HCC nodules were YFP(+) lineage-labeled; positive for osteopontin, SRY (sex determining region Y)-box 9, and epithelial cell adhesion molecule; and enriched for transcripts of biliary/progenitor markers such as prominin 1, Cd44, and delta-like 1 homolog. Next, we performed the converse experiment and lineage-labeled forkhead box protein L1(Foxl1)-positive hepatic progenitor cells simultaneously with exposure to carcinogens. None of the tumor nodules expressed YFP, indicating that Foxl1-expressing cells are not the origin for hepatotoxin-induced liver tumors. We confirmed that HCA and HCC cells are derived from mature hepatocytes and not from Foxl1-Cre-marked cells in a second model of toxin-induced hepatic neoplasia, using DEN and 3,3',5,5'-tetrachloro-1,4-bis(pyridyloxy)benzene (TCPOBOP). CONCLUSION: Hepatocytes are the cell of origin of HCA and HCC in DEN/carbon tetrachloride and DEN/TCPOBOP induced liver tumors. (Hepatology 2016;64:1163-1177).

ß-Cells are not uniform after all-Novel insights into molecular heterogeneity of insulin-secreting cells.

While the ß-cells of the endocrine pancreas are defined as cells with high levels of insulin production and tight stimulus-secretion coupling, the existence of functional heterogeneity among them has been known for decades. Recent advances in molecular technologies, in particular single-cell profiling on both the protein and messenger RNA level, have uncovered that ß-cells exist in several antigenically and molecularly definable states. Using antibodies to cell surface markers or multidimensional clustering of ß-cells using more than 20 protein markers by mass cytometry, 4 distinct groups of ß-cells could be differentiated. However, whether these states represent permanent cell lineages or are readily interconvertible from one group to another remains to be determined. Nevertheless, future analysis of the pathogenesis of type 1 and type 2 diabetes will certainly benefit from a growing appreciation of ß-cell heterogeneity. Here, we aim to summarize concisely the recent advances in the field and their possible impact on our understanding of ß-cell physiology and pathophysiology.

p120 Catenin is required for normal tubulogenesis but not epithelial integrity in developing mouse pancreas.

The intracellular protein p120 catenin aids in maintenance of cell-cell adhesion by regulating E-cadherin stability in epithelial cells. In an effort to understand the biology of p120 catenin in pancreas development, we ablated p120 catenin in mouse pancreatic progenitor cells, which resulted in deletion of p120 catenin in all epithelial lineages of the developing mouse pancreas: islet, acinar, centroacinar, and ductal. Loss of p120 catenin resulted in formation of dilated epithelial tubules, expansion of ductal epithelia, loss of acinar cells, and the induction of pancreatic inflammation. Aberrant branching morphogenesis and tubulogenesis were also observed. Throughout development, the phenotype became more severe, ultimately resulting in an abnormal pancreas comprised primarily of duct-like epithelium expressing early progenitor markers. In pancreatic tissue lacking p120 catenin, overall epithelial architecture remained intact; however, actin cytoskeleton organization was disrupted, an observation associated with increased cytoplasmic PKCζ. Although we observed reduced expression of adherens junction proteins E-cadherin, ß-catenin, and α-catenin, p120 catenin family members p0071, ARVCF, and δ-catenin remained present at cell membranes in homozygous p120(f/f) pancreases, potentially providing stability for maintenance of epithelial integrity during development. Adult mice homozygous for deletion of p120 catenin displayed dilated main pancreatic ducts, chronic pancreatitis, acinar to ductal metaplasia (ADM), and mucinous metaplasia that resembles PanIN1a. Taken together, our data demonstrate an essential role for p120 catenin in pancreas development.

Fate mapping of ptf1a-expressing cells during pancreatic organogenesis and regeneration in zebrafish.

BACKGROUND: Pancreas development in zebrafish shares many features with mammals, including the participation of epithelial progenitor cells expressing pancreas transcription factor 1a (ptf1a). However, to date it has remained unclear whether, as in mammals, ptf1a-expressing zebrafish pancreatic progenitors are able to contribute to multiple exocrine and endocrine lineages. To delineate the lineage potential of ptf1a-expressing cells, we generated ptf1a:creER(T2) transgenic fish and performed genetic-inducible lineage tracing in developmental, regenerating, and ptf1a-deficient zebrafish pancreas. RESULTS: In addition to their contribution to the acinar cell lineage, ptf1a-expressing cells give rise to both pancreatic Notch-responsive-cells (PNCs) as well as small numbers of endocrine cells during pancreatic development. In fish with ptf1a haploinsufficiency, a higher proportion of ptf1a lineage-labeled cells are traced into the PNC and endocrine compartments. Further reduction of ptf1a gene dosage converts pancreatic progenitor cells to gall bladder and other non-pancreatic cell fates. CONCLUSIONS: Our results confirm the presence of multipotent ptf1a-expressing progenitor cells in developing zebrafish pancreas, with reduced ptf1a dosage promoting greater contributions towards non-acinar lineages. As in mammals, loss of ptf1a results in conversion of nascent pancreatic progenitor cells to non-pancreatic cell fates, underscoring the central role of ptf1a in foregut tissue specification.

Comprehensive Characterization of Islet Remodeling in Development and in Diabetes Using Mass Cytometry.

The pancreatic islet is the functional and structural unit of the pancreatic endocrine portion. Islet remodeling occurs in both normal development and pathogenesis of type 1 (T1D) and type 2 diabetes (T2D). However, accurately quantifying changes in islet cellular makeup and hormone expressions poses significant challenges due to large intra- and inter-donor heterogeneity and the limited scalability of traditional methods such as immunostaining. The cytometry by time-of-flight (CyTOF) technology enables simultaneous quantification of more than 30 protein markers at single-cell resolution in a high-throughput fashion. Moreover, with distinct DNA and viability markers, single live cells can be explicitly selected in CyTOF. Here, leveraging the CyTOF data generated by the Human Pancreas Analysis Program, we characterized more than 12 million islet cells from 71 donors. Our data revealed continued age-related changes in islet endocrine cell compositions, but the maturity of endocrine cells is reached by 3 years of age. We also observed significant changes in beta cell numbers and key protein expressions, along with a significant increase in bihormonal cells in T1D donors. In contrast, T2D donors exhibited minimal islet remodeling events. Our data shine a light on the islet dynamics during development and diabetes pathogenesis and suggest divergent pathogenesis processes of T1D and T2D. Our comprehensive approach not only elucidates islet plasticity but also establishes a foundation for integrated CyTOF analysis in islet biology and beyond.

A single-cell atlas of the murine pancreatic ductal tree identifies novel cell populations with potential implications in pancreas regeneration and exocrine pathogenesis.

Background and aims: Pancreatic ducts form an intricate network of tubules that secrete bicarbonate and drive acinar secretions into the duodenum. This network is formed by centroacinar cells, terminal, intercalated, intracalated ducts, and the main pancreatic duct. Ductal heterogeneity at the single-cell level has been poorly characterized; therefore, our understanding of the role of ductal cells in pancreas regeneration and exocrine pathogenesis has been hampered by the limited knowledge and unexplained diversity within the ductal network. Methods: We used scRNA-seq to comprehensively characterize mouse ductal heterogeneity at single-cell resolution of the entire ductal epithelium from centroacinar cells to the main duct. Moreover, we used organoid cultures, injury models and pancreatic tumor samples to interrogate the role of novel ductal populations in pancreas regeneration and exocrine pathogenesis. Results: We have identified the coexistence of 15 ductal populations within the healthy pancreas and characterized their organoid formation capacity and endocrine differentiation potential. Cluster isolation and subsequent culturing let us identify ductal cell populations with high organoid formation capacity and endocrine and exocrine differentiation potential in vitro , including Wnt-responsive-population, ciliated-population and FLRT3 + cells. Moreover, we have characterized the location of these novel ductal populations in healthy pancreas, chronic pancreatitis, and tumor samples, highlighting a putative role of WNT-responsive, IFN-responsive and EMT-populations in pancreatic exocrine pathogenesis as their expression increases in chronic pancreatitis and PanIN lesions. Conclusions: In light of our discovery of previously unidentified ductal populations, we unmask the potential roles of specific ductal populations in pancreas regeneration and exocrine pathogenesis.

AnnoSpat annotates cell types and quantifies cellular arrangements from spatial proteomics.

Cellular composition and anatomical organization influence normal and aberrant organ functions. Emerging spatial single-cell proteomic assays such as Image Mass Cytometry (IMC) and Co-Detection by Indexing (CODEX) have facilitated the study of cellular composition and organization by enabling high-throughput measurement of cells and their localization directly in intact tissues. However, annotation of cell types and quantification of their relative localization in tissues remain challenging. To address these unmet needs for atlas-scale datasets like Human Pancreas Analysis Program (HPAP), we develop AnnoSpat (Annotator and Spatial Pattern Finder) that uses neural network and point process algorithms to automatically identify cell types and quantify cell-cell proximity relationships. Our study of data from IMC and CODEX shows the higher performance of AnnoSpat in rapid and accurate annotation of cell types compared to alternative approaches. Moreover, the application of AnnoSpat to type 1 diabetic, non-diabetic autoantibody-positive, and non-diabetic organ donor cohorts recapitulates known islet pathobiology and shows differential dynamics of pancreatic polypeptide (PP) cell abundance and CD8+ T cells infiltration in islets during type 1 diabetes progression.

AnnoSpat annotates cell types and quantifies cellular arrangements from spatial proteomics.

Cellular composition and anatomical organization influence normal and aberrant organ functions. Emerging spatial single-cell proteomic assays such as Image Mass Cytometry (IMC) and Co-Detection by Indexing (CODEX) have facilitated the study of cellular composition and organization by enabling high-throughput measurement of cells and their localization directly in intact tissues. However, annotation of cell types and quantification of their relative localization in tissues remain challenging. To address these unmet needs, we developed AnnoSpat (Annotator and Spatial Pattern Finder) that uses neural network and point process algorithms to automatically identify cell types and quantify cell-cell proximity relationships. Our study of data from IMC and CODEX show the superior performance of AnnoSpat in rapid and accurate annotation of cell types compared to alternative approaches. Moreover, the application of AnnoSpat to type 1 diabetic, non-diabetic autoantibody-positive, and non-diabetic organ donor cohorts recapitulated known islet pathobiology and showed differential dynamics of pancreatic polypeptide (PP) cell abundance and CD8+ T cells infiltration in islets during type 1 diabetes progression.

Increasing insulin measurement throughput by fluorescence anisotropy imaging immunoassays.

Insulin secreted from islets of Langerhans is the main hormone to reduce blood glucose. Examination of insulin secretion patterns at the single islet level reveals functional differences in the timings and patterns of release. This heterogeneous response highlights the importance of developing systems to measure dynamic release from small numbers of islets in parallel. Toward this, we describe fluorescence anisotropy imaging immunoassays as a relatively simple method for increased throughput of islet secretion measurements. In this system, vacuum pressure from a syringe pump pulled perfusate from 12 islet chambers and reagents into 12 parallel mixing channels for a competitive immunoassay. Light from a Xe arc lamp was filtered and polarized prior to focusing on the microfluidic device at the region where the 12 mixing channels converged. Emission was collected and passed through vertical and horizontal emission polarizers housed in an automated filter wheel before being imaged with a sCMOS camera for the determination of anisotropy. This microfluidic system was tested by monitoring insulin release from groups of murine and human islets. Heterogeneity was observed in the islet traces; however, the presence of islets affected the resistance of the islet chambers, hampering insulin quantification. Nonetheless, this microfluidic system is a step towards increasing the throughput of hormone release measurements from islets of Langerhans.

Sex Differences in the Molecular Programs of Pancreatic Cells Contribute to the Differential Risks of Type 2 Diabetes.

Epidemiology studies demonstrate that women are at a significantly lower risk of developing type 2 diabetes (T2D) compared to men. However, the molecular basis of this risk difference is not well understood. In this study, we examined the sex differences in the genetic programs of pancreatic endocrine cells. We combined pancreas perifusion data and single-cell genomic data from our laboratory and from publicly available data sets to investigate multiple axes of the sex differences in the human pancreas at the single-cell type and single-cell level. We systematically compared female and male islet secretion function, gene expression program, and regulatory principles of pancreatic endocrine cells. The perifusion data indicate that female endocrine cells have a higher secretion capacity than male endocrine cells. Single-cell RNA-sequencing analysis suggests that endocrine cells in male controls have molecular signatures that resemble T2D. In addition, we identified genomic elements associated with genome-wide association study T2D loci to have differential accessibility between female and male delta cells. These genomic elements may play a sex-specific causal role in the pathogenesis of T2D. We provide molecular mechanisms that explain the differential risk of T2D between women and men. Knowledge gained from our study will accelerate the development of diagnostics and therapeutics in sex-aware precision medicine for diabetes.

Single-cell multi-omics analysis of human pancreatic islets reveals novel cellular states in type 1 diabetes.

Type 1 diabetes (T1D) is an autoimmune disease in which immune cells destroy insulin-producing beta cells. The aetiology of this complex disease is dependent on the interplay of multiple heterogeneous cell types in the pancreatic environment. Here, we provide a single-cell atlas of pancreatic islets of 24 T1D, autoantibody-positive and nondiabetic organ donors across multiple quantitative modalities including ~80,000 cells using single-cell transcriptomics, ~7,000,000 cells using cytometry by time of flight and ~1,000,000 cells using in situ imaging mass cytometry. We develop an advanced integrative analytical strategy to assess pancreatic islets and identify canonical cell types. We show that a subset of exocrine ductal cells acquires a signature of tolerogenic dendritic cells in an apparent attempt at immune suppression in T1D donors. Our multimodal analyses delineate cell types and processes that may contribute to T1D immunopathogenesis and provide an integrative procedure for exploration and discovery of human pancreatic function.

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.