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.

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.

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

Highly multiplexed 2-dimensional imaging mass cytometry analysis of HBV-infected liver.

Studies of human hepatitis B virus (HBV) immune pathogenesis are hampered by limited access to liver tissues and technologies for detailed analyses. Here, utilizing imaging mass cytometry (IMC) to simultaneously detect 30 immune, viral, and structural markers in liver biopsies from patients with hepatitis B e antigen+ (HBeAg+) chronic hepatitis B, we provide potentially novel comprehensive visualization, quantitation, and phenotypic characterizations of hepatic adaptive and innate immune subsets that correlated with hepatocellular injury, histological fibrosis, and age. We further show marked correlations between adaptive and innate immune cell frequencies and phenotype, highlighting complex immune interactions within the hepatic microenvironment with relevance to HBV pathogenesis.

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.

The Dynamic Chromatin Architecture of the Regenerating Liver.

BACKGROUND & AIMS: The adult liver is the main detoxification organ and routinely is exposed to environmental insults but retains the ability to restore its mass and function upon tissue damage. However, extensive injury can lead to liver failure, and chronic injury causes fibrosis, cirrhosis, and hepatocellular carcinoma. Currently, the transcriptional regulation of organ repair in the adult liver is incompletely understood. METHODS: We isolated nuclei from quiescent as well as repopulating hepatocytes in a mouse model of hereditary tyrosinemia, which recapitulates the injury and repopulation seen in toxic liver injury in human beings. We then performed the assay for transposase accessible chromatin with high-throughput sequencing specifically in repopulating hepatocytes to identify differentially accessible chromatin regions and nucleosome positioning. In addition, we used motif analysis to predict differential transcription factor occupancy and validated the in silico results with chromatin immunoprecipitation followed by sequencing for hepatocyte nuclear factor 4α (HNF4α) and CCCTC-binding factor (CTCF). RESULTS: Chromatin accessibility in repopulating hepatocytes was increased in the regulatory regions of genes promoting proliferation and decreased in the regulatory regions of genes involved in metabolism. The epigenetic changes at promoters and liver enhancers correspond with the regulation of gene expression, with enhancers of many liver function genes showing a less accessible state during the regenerative process. Moreover, increased CTCF occupancy at promoters and decreased HNF4α binding at enhancers implicate these factors as key drivers of the transcriptomic changes in replicating hepatocytes that enable liver repopulation. CONCLUSIONS: Our analysis of hepatocyte-specific epigenomic changes during liver repopulation identified CTCF and HNF4α as key regulators of hepatocyte proliferation and regulation of metabolic programs. Thus, liver repopulation in the setting of toxic injury makes use of both general transcription factors (CTCF) for promoter activation, and reduced binding by a hepatocyte-enriched factor (HNF4α) to temporarily limit enhancer activity. All sequencing data in this study were deposited to the Gene Expression Omnibus database and can be downloaded with accession number GSE109466.

TRAP-seq identifies cystine/glutamate antiporter as a driver of recovery from liver injury.

Understanding the molecular basis of the regenerative response following hepatic injury holds promise for improved treatment of liver diseases. Here, we report an innovative method to profile gene expression specifically in the hepatocytes that regenerate the liver following toxic injury. We used the Fah-/- mouse, a model of hereditary tyrosinemia, which conditionally undergoes severe liver injury unless fumarylacetoacetate hydrolase (FAH) expression is reconstituted ectopically. We used translating ribosome affinity purification followed by high-throughput RNA sequencing (TRAP-seq) to isolate mRNAs specific to repopulating hepatocytes. We uncovered upstream regulators and important signaling pathways that are highly enriched in genes changed in regenerating hepatocytes. Specifically, we found that glutathione metabolism, particularly the gene Slc7a11 encoding the cystine/glutamate antiporter (xCT), is massively upregulated during liver regeneration. Furthermore, we show that Slc7a11 overexpression in hepatocytes enhances, and its suppression inhibits, repopulation following toxic injury. TRAP-seq allows cell type-specific expression profiling in repopulating hepatocytes and identified xCT, a factor that supports antioxidant responses during liver regeneration. xCT has potential as a therapeutic target for enhancing liver regeneration in response to liver injury.

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

p120 Catenin Suppresses Basal Epithelial Cell Extrusion in Invasive Pancreatic Neoplasia.

Aberrant regulation of cellular extrusion can promote invasion and metastasis. Here, we identify molecular requirements for early cellular invasion using a premalignant mouse model of pancreatic cancer with conditional knockout of p120 catenin (Ctnnd1). Mice with biallelic loss of p120 catenin progressively develop high-grade pancreatic intraepithelial neoplasia (PanIN) lesions and neoplasia accompanied by prominent acute and chronic inflammatory processes, which is mediated, in part, through NF-κB signaling. Loss of p120 catenin in the context of oncogenic Kras also promotes remarkable apical and basal epithelial cell extrusion. Abundant single epithelial cells exit PanIN epithelium basally, retain epithelial morphology, survive, and display features of malignancy. Similar extrusion defects are observed following p120 catenin knockdown in vitro, and these effects are completely abrogated by the activation of S1P/S1pr2 signaling. In the context of oncogenic Kras, p120 catenin loss significantly reduces expression of genes mediating S1P/S1pr2 signaling in vivo and in vitro, and this effect is mediated at least, in part, through activation of NF-κB. These results provide insight into mechanisms controlling early events in the metastatic process and suggest that p120 catenin and S1P/S1pr2 signaling enhance cancer progression by regulating epithelial cell invasion. Cancer Res; 76(11); 3351-63. ©2016 AACR.

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.

Dicer is required for maintenance of adult pancreatic acinar cell identity and plays a role in Kras-driven pancreatic neoplasia.

The role of miRNA processing in the maintenance of adult pancreatic acinar cell identity and during the initiation and progression of pancreatic neoplasia has not been studied in detail. In this work, we deleted Dicer specifically in adult pancreatic acinar cells, with or without simultaneous activation of oncogenic Kras. We found that Dicer is essential for the maintenance of acinar cell identity. Acinar cells lacking Dicer showed increased plasticity, as evidenced by loss of polarity, initiation of epithelial-to-mesenchymal transition (EMT) and acinar-to-ductal metaplasia (ADM). In the context of oncogenic Kras activation, the initiation of ADM and pancreatic intraepithelial neoplasia (PanIN) were both highly sensitive to Dicer gene dosage. Homozygous Dicer deletion accelerated the formation of ADM but not PanIN. In contrast, heterozygous Dicer deletion accelerated PanIN initiation, revealing complex roles for Dicer in the regulation of both normal and neoplastic pancreatic epithelial identity.

In vitro screening for compounds that enhance human L1 mobilization.

The Long interspersed element 1 (LINE1 or L1) retrotransposon constitutes 17% of the human genome. There are currently 80-100 human L1 elements that are thought to be active in any diploid human genome. These elements can mobilize into new locations of the genome, resulting in changes in genomic information. Active L1s are thus considered to be a type of endogenous mutagen, and L1 insertions can cause disease. Certain stresses, such as gamma radiation, oxidative stress, and treatment with some agents, can induce transcription and/or mobilization of retrotransposons. In this study, we used a reporter gene assay in HepG2 cells to screen compounds for the potential to enhance the transcription of human L1. We assessed 95 compounds including genotoxic agents, substances that induce cellular stress, and commercially available drugs. Treatment with 15 compounds increased the L1 promoter activity by >1.5-fold (p<0.05) after 6 or 24 hours of treatment. In particular, genotoxic agents (benzo[a]pyrene, camptothecin, cytochalasin D, merbarone, and vinblastine), PPARα agonists (bezafibrate and fenofibrate), and non-steroidal anti-inflammatory drugs (diflunisal, flufenamic acid, salicylamide, and sulindac) induced L1 promoter activity. To examine their effects on L1 retrotransposition, we developed a high-throughput real-time retrotransposition assay using a novel secreted Gaussia luciferase reporter cassette. Three compounds (etomoxir, WY-14643, and salicylamide) produced a significant enhancement in L1 retrotransposition. This is the first study to report the effects of a wide variety of compounds on L1 transcription and retrotransposition. These results suggest that certain chemical- and drug-induced stresses might have the potential to cause genomic mutations by inducing L1 mobilization. Thus, the risk of induced L1 transcription and retrotransposition should be considered during drug safety evaluation and environmental risk assessments of chemicals.

Sphere-forming assays for assessment of benign and malignant pancreatic stem cells.

Sphere-forming assays are an in vitro technique to assay both normal and neoplastic cells for clonogenic growth potential. Currently, the identification of adult progenitors in the pancreas remains an area of intense investigation. The use of sphere-forming assays provides a critical step to identify new cell types in the pancreas that are capable of clonogenic growth and differentiation. In the field of cancer biology, cancer stem cells have been defined functionally by two major criteria: their ability to undergo self-renewal and their ability to produce differentiated progeny, two conditions which satisfy the criteria of stem cells. Here we briefly review both the capabilities of pancreatosphere and pancreatic tumorsphere assays, discuss important caveats regarding their use, and provide detailed protocols for the assay of both normal and neoplastic cells.