Novel Approach for Pancreas Transcriptomics Reveals the Cellular Landscape in Homeostasis and Acute Pancreatitis.

BACKGROUND & AIMS: Acinar cells produce digestive enzymes that impede transcriptomic characterization of the exocrine pancreas. Thus, single-cell RNA-sequencing studies of the pancreas underrepresent acinar cells relative to histological expectations, and a robust approach to capture pancreatic cell responses in disease states is needed. We sought to innovate a method that overcomes these challenges to accelerate study of the pancreas in health and disease. METHODS: We leverage FixNCut, a single-cell RNA-sequencing approach in which tissue is reversibly fixed with dithiobis(succinimidyl propionate) before dissociation and single-cell preparation. We apply FixNCut to an established mouse model of acute pancreatitis, validate findings using GeoMx whole transcriptome atlas profiling, and integrate our data with prior studies to compare our method in both mouse and human pancreas datasets. RESULTS: FixNCut achieves unprecedented definition of challenging pancreatic cells, including acinar and immune populations in homeostasis and acute pancreatitis, and identifies changes in all major cell types during injury and recovery. We define the acinar transcriptome during homeostasis and acinar-to-ductal metaplasia and establish a unique gene set to measure deviation from normal acinar identity. We characterize pancreatic immune cells, and analysis of T-cell subsets reveals a polarization of the homeostatic pancreas toward type-2 immunity. We report immune responses during acute pancreatitis and recovery, including early neutrophil infiltration, expansion of dendritic cell subsets, and a substantial shift in the transcriptome of macrophages due to both resident macrophage activation and monocyte infiltration. CONCLUSIONS: FixNCut preserves pancreatic transcriptomes to uncover novel cell states during homeostasis and following pancreatitis, establishing a broadly applicable approach and reference atlas for study of pancreas biology and disease.

Chronic metabolic stress drives developmental programs and loss of tissue functions in non-transformed liver that mirror tumor states and stratify survival.

Under chronic stress, cells must balance competing demands between cellular survival and tissue function. In metabolic dysfunction-associated steatotic liver disease (MASLD, formerly NAFLD/NASH), hepatocytes cooperate with structural and immune cells to perform crucial metabolic, synthetic, and detoxification functions despite nutrient imbalances. While prior work has emphasized stress-induced drivers of cell death, the dynamic adaptations of surviving cells and their functional repercussions remain unclear. Namely, we do not know which pathways and programs define cellular responses, what regulatory factors mediate (mal)adaptations, and how this aberrant activity connects to tissue-scale dysfunction and long-term disease outcomes. Here, by applying longitudinal single-cell multi -omics to a mouse model of chronic metabolic stress and extending to human cohorts, we show that stress drives survival-linked tradeoffs and metabolic rewiring, manifesting as shifts towards development-associated states in non-transformed hepatocytes with accompanying decreases in their professional functionality. Diet-induced adaptations occur significantly prior to tumorigenesis but parallel tumorigenesis-induced phenotypes and predict worsened human cancer survival. Through the development of a multi -omic computational gene regulatory inference framework and human in vitro and mouse in vivo genetic perturbations, we validate transcriptional (RELB, SOX4) and metabolic (HMGCS2) mediators that co-regulate and couple the balance between developmental state and hepatocyte functional identity programming. Our work defines cellular features of liver adaptation to chronic stress as well as their links to long-term disease outcomes and cancer hallmarks, unifying diverse axes of cellular dysfunction around core causal mechanisms.

The complex, dynamic SpliceOme of the small GTPase transcripts altered by technique, sex, genetics, tissue specificity, and RNA base editing.

The small GTPase family is well-studied in cancer and cellular physiology. With 162 annotated human genes, the family has a broad expression throughout cells of the body. Members of the family have multiple exons that require splicing. Yet, the role of splicing within the family has been underexplored. We have studied the splicing dynamics of small GTPases throughout 41,671 samples by integrating Nanopore and Illumina sequencing techniques. Within this work, we have made several discoveries. 1). Using the GTEx long read data of 92 samples, each small GTPase gene averages two transcripts, with 83 genes (51%) expressing two or more isoforms. 2). Cross-tissue analysis of GTEx from 17,382 samples shows 41 genes (25%) expressing two or more protein-coding isoforms. These include protein-changing transcripts in genes such as RHOA, RAB37, RAB40C, RAB4B, RAB5C, RHOC, RAB1A, RAN, RHEB, RAC1, and KRAS. 3). The isolation and library technique of the RNAseq influences the abundance of non-sense-mediated decay and retained intron transcripts of small GTPases, which are observed more often in genes than appreciated. 4). Analysis of 16,243 samples of "Blood PAXgene" identified seven genes (3.7%; RHOA, RAB40C, RAB4B, RAB37, RAB5B, RAB5C, RHOC) with two or more transcripts expressed as the major isoform (75% of the total gene), suggesting a role of genetics in altering splicing. 5). Rare (ARL6, RAB23, ARL13B, HRAS, NRAS) and common variants (GEM, RHOC, MRAS, RAB5B, RERG, ARL16) can influence splicing and have an impact on phenotypes and diseases. 6). Multiple genes (RAB9A, RAP2C, ARL4A, RAB3A, RAB26, RAB3C, RASL10A, RAB40B, and HRAS) have sex differences in transcript expression. 7). Several exons are included or excluded for small GTPase genes (RASEF, KRAS, RAC1, RHEB, ARL4A, RHOA, RAB30, RHOBTB1, ARL16, RAP1A) in one or more forms of cancer. 8). Ten transcripts are altered in hypoxia (SAR1B, IFT27, ARL14, RAB11A, RAB10, RAB38, RAN, RIT1, RAB9A) with RHOA identified to have a transient 3'UTR RNA base editing at a conserved site found in all of its transcripts. Overall, we show a remarkable and dynamic role of splicing within the small GTPase family that requires future explorations.

Germline Testing for Individuals with Pancreatic Adenocarcinoma and Novel Genetic Risk Factors.

Germline genetic variants implicated in increasing lifetime risk of pancreatic cancer (PDAC) have been identified in ∼4% to 10% of cases. Clinical features such as family history have poor sensitivity in identifying carriers of these risk variants. Genetic testing for these germline variants has potential to guide risk assessment and surveillance recommendations in high-risk individuals to promote prevention and early detection measures. Furthermore, identification of novel germline variants can offer important insights into pathogenesis that may inform precision medicine approaches. This article reviews current understanding of germline mutations associated with PDAC risk and implications of genetic testing.

More than acinar identity? A novel cystic phenotype suggests broader roles for NR5A2 in pancreatic cancer†.

The prognosis for pancreatic ductal adenocarcinoma (PDAC) remains dismal. Multiple genome-wide association studies (GWAS) have implicated the nuclear receptor NR5A2 in modulating PDAC risk, but mechanisms for this association are not understood. NR5A2 is a transcription factor that maintains acinar cell identity, and heterozygous loss of Nr5a2 in mice accelerates oncogenic Kras-driven formation of pancreatic intraepithelial neoplasia (PanIN), a PDAC precursor derived from acinar cells. In a recent issue of The Journal of Pathology, Cobo et al characterize a novel mouse model that uses Ptf1a:Cre to drive oncogenic Kras as well as heterozygous Nr5a2 inactivation. In addition to the expected PanIN lesions, these mice exhibited a surprising phenotype: large pancreatic cystic lesions which have not been previously reported. Comparing expression of oncogenic Kras and heterozygous Nr5a2 in various mouse models reveals several possible explanations for these cystic lesions. Importantly, these differences across mouse models suggest that NR5A2 may contribute to PDAC precursors in ways beyond its previously characterized acinar cell-autonomous role. These observations highlight that pathways implicated by GWAS may have roles in unexpected cell types, and an understanding of these roles will be critical to guide new preventive and treatment strategies for PDAC. © 2021 The Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Mutations in RABL3 alter KRAS prenylation and are associated with hereditary pancreatic cancer.

Pancreatic ductal adenocarcinoma is an aggressive cancer with limited treatment options1. Approximately 10% of cases exhibit familial predisposition, but causative genes are not known in most families2. We perform whole-genome sequence analysis in a family with multiple cases of pancreatic ductal adenocarcinoma and identify a germline truncating mutation in the member of the RAS oncogene family-like 3 (RABL3) gene. Heterozygous rabl3 mutant zebrafish show increased susceptibility to cancer formation. Transcriptomic and mass spectrometry approaches implicate RABL3 in RAS pathway regulation and identify an interaction with RAP1GDS1 (SmgGDS), a chaperone regulating prenylation of RAS GTPases3. Indeed, the truncated mutant RABL3 protein accelerates KRAS prenylation and requires RAS proteins to promote cell proliferation. Finally, evidence in patient cohorts with developmental disorders implicates germline RABL3 mutations in RASopathy syndromes. Our studies identify RABL3 mutations as a target for genetic testing in cancer families and uncover a mechanism for dysregulated RAS activity in development and cancer.

Selenoprotein H is an essential regulator of redox homeostasis that cooperates with p53 in development and tumorigenesis.

Selenium, an essential micronutrient known for its cancer prevention properties, is incorporated into a class of selenocysteine-containing proteins (selenoproteins). Selenoprotein H (SepH) is a recently identified nucleolar oxidoreductase whose function is not well understood. Here we report that seph is an essential gene regulating organ development in zebrafish. Metabolite profiling by targeted LC-MS/MS demonstrated that SepH deficiency impairs redox balance by reducing the levels of ascorbate and methionine, while increasing methionine sulfoxide. Transcriptome analysis revealed that SepH deficiency induces an inflammatory response and activates the p53 pathway. Consequently, loss of seph renders larvae susceptible to oxidative stress and DNA damage. Finally, we demonstrate that seph interacts with p53 deficiency in adulthood to accelerate gastrointestinal tumor development. Overall, our findings establish that seph regulates redox homeostasis and suppresses DNA damage. We hypothesize that SepH deficiency may contribute to the increased cancer risk observed in cohorts with low selenium levels.

Yap reprograms glutamine metabolism to increase nucleotide biosynthesis and enable liver growth.

The Hippo pathway is an important regulator of organ size and tumorigenesis. It is unclear, however, how Hippo signalling provides the cellular building blocks required for rapid growth. Here, we demonstrate that transgenic zebrafish expressing an activated form of the Hippo pathway effector Yap1 (also known as YAP) develop enlarged livers and are prone to liver tumour formation. Transcriptomic and metabolomic profiling identify that Yap1 reprograms glutamine metabolism. Yap1 directly enhances glutamine synthetase (glul) expression and activity, elevating steady-state levels of glutamine and enhancing the relative isotopic enrichment of nitrogen during de novo purine and pyrimidine biosynthesis. Genetic or pharmacological inhibition of GLUL diminishes the isotopic enrichment of nitrogen into nucleotides, suppressing hepatomegaly and the growth of liver cancer cells. Consequently, Yap-driven liver growth is susceptible to nucleotide inhibition. Together, our findings demonstrate that Yap1 integrates the anabolic demands of tissue growth during development and tumorigenesis by reprogramming nitrogen metabolism to stimulate nucleotide biosynthesis.

Iterative use of nuclear receptor Nr5a2 regulates multiple stages of liver and pancreas development.

The stepwise progression of common endoderm progenitors into differentiated liver and pancreas organs is regulated by a dynamic array of signals that are not well understood. The nuclear receptor subfamily 5, group A, member 2 gene nr5a2, also known as Liver receptor homolog-1 (Lrh-1) is expressed in several tissues including the developing liver and pancreas. Here, we interrogate the role of Nr5a2 at multiple developmental stages using genetic and chemical approaches and uncover novel pleiotropic requirements during zebrafish liver and pancreas development. Zygotic loss of nr5a2 in a targeted genetic null mutant disrupted the development of the exocrine pancreas and liver, while leaving the endocrine pancreas intact. Loss of nr5a2 abrogated exocrine pancreas markers such as trypsin, while pancreas progenitors marked by ptf1a or pdx1 remained unaffected, suggesting a role for Nr5a2 in regulating pancreatic acinar cell differentiation. In the developing liver, Nr5a2 regulates hepatic progenitor outgrowth and differentiation, as nr5a2 mutants exhibited reduced hepatoblast markers hnf4α and prox1 as well as differentiated hepatocyte marker fabp10a. Through the first in vivo use of Nr5a2 chemical antagonist Cpd3, the iterative requirement for Nr5a2 for exocrine pancreas and liver differentiation was temporally elucidated: chemical inhibition of Nr5a2 function during hepatopancreas progenitor specification was sufficient to disrupt exocrine pancreas formation and enhance the size of the embryonic liver, suggesting that Nr5a2 regulates hepatic vs. pancreatic progenitor fate choice. Chemical inhibition of Nr5a2 at a later time during pancreas and liver differentiation was sufficient to block the formation of mature acinar cells and hepatocytes. These findings define critical iterative and pleiotropic roles for Nr5a2 at distinct stages of pancreas and liver organogenesis, and provide novel perspectives for interpreting the role of Nr5a2 in disease.

Ferritinophagy via NCOA4 is required for erythropoiesis and is regulated by iron dependent HERC2-mediated proteolysis.

NCOA4 is a selective cargo receptor for the autophagic turnover of ferritin, a process critical for regulation of intracellular iron bioavailability. However, how ferritinophagy flux is controlled and the roles of NCOA4 in iron-dependent processes are poorly understood. Through analysis of the NCOA4-FTH1 interaction, we demonstrate that direct association via a key surface arginine in FTH1 and a C-terminal element in NCOA4 is required for delivery of ferritin to the lysosome via autophagosomes. Moreover, NCOA4 abundance is under dual control via autophagy and the ubiquitin proteasome system. Ubiquitin-dependent NCOA4 turnover is promoted by excess iron and involves an iron-dependent interaction between NCOA4 and the HERC2 ubiquitin ligase. In zebrafish and cultured cells, NCOA4 plays an essential role in erythroid differentiation. This work reveals the molecular nature of the NCOA4-ferritin complex and explains how intracellular iron levels modulate NCOA4-mediated ferritinophagy in cells and in an iron-dependent physiological setting.

Cannabinoid Receptor-2 Regulates Embryonic Hematopoietic Stem Cell Development via Prostaglandin E2 and P-Selectin Activity.

Cannabinoids (CB) modulate adult hematopoietic stem and progenitor cell (HSPCs) function, however, impact on the production, expansion, or migration of embryonic HSCs is currently uncharacterized. Here, using chemical and genetic approaches targeting CB-signaling in zebrafish, we show that CB receptor (CNR) 2, but not CNR1, regulates embryonic HSC development. During HSC specification in the aorta-gonad-mesonephros (AGM) region, CNR2 stimulation by AM1241 increased runx1;cmyb(+) HSPCs, through heightened proliferation, whereas CNR2 antagonism decreased HSPC number; FACS analysis and absolute HSC counts confirmed and quantified these effects. Epistatic investigations showed AM1241 significantly upregulated PGE2 synthesis in a Ptgs2-dependent manner to increase AGM HSCs. During the phases of HSC production and colonization of secondary niches, AM1241 accelerated migration to the caudal hematopoietic tissue (CHT), the site of embryonic HSC expansion, and the thymus; however these effects occurred independently of PGE2. Using a candidate approach for HSC migration and retention factors, P-selectin was identified as the functional target of CNR2 regulation. Epistatic analyses confirmed migration of HSCs into the CHT and thymus was dependent on CNR2-regulated P-selectin activity. Together, these data suggest CNR2-signaling optimizes the production, expansion, and migration of embryonic HSCs by modulating multiple downstream signaling pathways.

Estrogen defines the dorsal-ventral limit of VEGF regulation to specify the location of the hemogenic endothelial niche.

Genetic control of hematopoietic stem and progenitor cell (HSPC) function is increasingly understood; however, less is known about the interactions specifying the embryonic hematopoietic niche. Here, we report that 17ß-estradiol (E2) influences production of runx1+ HSPCs in the AGM region by antagonizing VEGF signaling and subsequent assignment of hemogenic endothelial (HE) identity. Exposure to exogenous E2 during vascular niche development significantly disrupted flk1+ vessel maturation, ephrinB2+ arterial identity, and specification of scl+ HE by decreasing expression of VEGFAa and downstream arterial Notch-pathway components; heat shock induction of VEGFAa/Notch rescued E2-mediated hematovascular defects. Conversely, repression of endogenous E2 activity increased somitic VEGF expression and vascular target regulation, shifting assignment of arterial/venous fate and HE localization; blocking E2 signaling allowed venous production of scl+/runx1+ cells, independent of arterial identity acquisition. Together, these data suggest that yolk-derived E2 sets the ventral boundary of hemogenic vascular niche specification by antagonizing the dorsal-ventral regulatory limits of VEGF.

Genetic markers of malignant transformation in intraductal papillary mucinous neoplasm of the pancreas: a meta-analysis.

OBJECTIVES: The objective of this study was to determine the relationship between specific genetic alterations and malignant transformation in intraductal papillary mucinous neoplasm (IPMN) of the pancreas. METHODS: Quantitative meta-analysis was conducted of studies through October 2010 that adhered to the 1996 World Health Organization guidelines for distinguishing adenoma and borderline IPMN versus carcinoma in surgically resected specimens using a random-effects model. We developed a 6-point scoring system to assess study quality. RESULTS: Thirty-nine studies (1235 IPMN samples) satisfied the inclusion criteria, and we conducted pooled analysis of 8 genetic markers: MUC1, MUC2, MUC5AC, kRas, p53, hTERT (human telomerase reverse transcriptase), cyclooxygenase 2, and Shh (Sonic hedgehog). Markers having the strongest association with malignant IPMN were hTERT (odds ratio [OR], 11.4; 95% confidence interval [CI], 3.5-36.7) and Shh (OR, 6.9; 95% CI, 2.4-20.2), whereas MUC5AC (OR, 1.0; 95% CI, 0.1-13.9) and kRas (OR, 2.0; 95% CI, 1.0-4.3) showed weak association with IPMN histologic progression. CONCLUSIONS: Expression of hTERT is strongly associated with malignant transformation in IPMN, consistent with up-regulation of hTERT as a key step in progression of IPMN to cancer. Expression of kRas and MUC5AC is common but not strongly associated with IPMN histologic progression. The quality criteria used here may guide future reporting of genetic markers related to malignant transformation of IPMN.

Regulation of Gremlin expression in the posterior limb bud.

Proper outgrowth of the limb bud requires a positive feedback loop between Sonic hedgehog (Shh) in the zone of polarizing activity (ZPA) and Fgfs in the overlying apical ectodermal ridge. The Bmp antagonist Gremlin is expressed in a domain anterior to the ZPA and is thought to act as a signaling intermediate between Shh and Fgf. It is currently unclear whether Shh acts directly or indirectly to initiate and maintain Gremlin. In this study, we confirm that Bmp activity is necessary and sufficient for induction of Gremlin. Beads soaked in the Bmp antagonist Noggin downregulate Gremlin, while beads soaked in Bmp2 cause its upregulation. Furthermore, Bmp2 is also capable of upregulating Gremlin in oligozeugodactyly mutant limbs that lack Shh activity, demonstrating that Gremlin expression does not depend on the combined exposure to both these factors. In spite of the ability of Bmp2 to induce Gremlin, beads soaked in high concentrations of Bmp2 downregulate Gremlin around the bead without apparent induction of cell death, whereas another target gene Msx2 is upregulated around the bead. Consistent with this concentration-dependent effect, we find that low concentrations of Bmp2 upregulate Gremlin while high concentrations of Bmp2 downregulate Gremlin in limb mesenchyme cultures. These data implicate Bmp activity as a required intermediate in the Shh-Fgf4 signaling loop. Though we show that Bmp activity is sufficient to upregulate Gremlin, Gremlin expression is excluded from a posterior domain of the limb, and expansion of this domain as limb outgrowth proceeds is important in terminating the Shh-Fgf4 signaling loop. We find that the posterior limb is refractory to Gremlin induction in response to Bmp2, suggesting that termination of the Shh-Fgf4 signaling loop results from inability of Bmp activity to induce Gremlin in the posterior. In contrast, in the oligozeugodactyly limb, we find that beads soaked in Bmp2 can induce Gremlin in the posterior, demonstrating that Shh activity is required for exclusion of Gremlin in the posterior. Finally, by blocking Shh activity with cyclopamine, we find evidence that continued Shh activity is also required to maintain refractoriness to Gremlin expression in response to Bmp activity.

The Hedgehog-inducible ubiquitin ligase subunit WSB-1 modulates thyroid hormone activation and PTHrP secretion in the developing growth plate.

WSB-1 is a SOCS-box-containing WD-40 protein of unknown function that is induced by Hedgehog signalling in embryonic structures during chicken development. Here we show that WSB-1 is part of an E3 ubiquitin ligase for the thyroid-hormone-activating type 2 iodothyronine deiodinase (D2). The WD-40 propeller of WSB-1 recognizes an 18-amino-acid loop in D2 that confers metabolic instability, whereas the SOCS-box domain mediates its interaction with a ubiquitinating catalytic core complex, modelled as Elongin BC-Cul5-Rbx1 (ECS(WSB-1)). In the developing tibial growth plate, Hedgehog-stimulated D2 ubiquitination via ECS(WSB-1) induces parathyroid hormone-related peptide (PTHrP), thereby regulating chondrocyte differentiation. Thus, ECS(WSB-1) mediates a mechanism by which 'systemic' thyroid hormone can effect local control of the Hedgehog-PTHrP negative feedback loop and thus skeletogenesis.

Evidence for an expansion-based temporal Shh gradient in specifying vertebrate digit identities.

The zone of polarizing activity (ZPA) in the posterior limb bud produces Sonic Hedgehog (Shh) protein, which plays a critical role in establishing distinct fates along the anterior-posterior axis. This activity has been modeled as a concentration-dependent response to a diffusible morphogen. Using recombinase base mapping in the mouse, we determine the ultimate fate of the Shh-producing cells. Strikingly, the descendants of the Shh-producing cells encompass all cells in the two most posterior digits and also contribute to the middle digit. Our analysis suggests that, while specification of the anterior digits depends upon differential concentrations of Shh, the length of time of exposure to Shh is critical in the specification of the differences between the most posterior digits. Genetic studies of the effects of limiting accessibility of Shh within the limb support this model, in which the effect of the Shh morphogen is dictated by a temporal as well as a spatial gradient.