Biliary Tract Cancer: Molecular Biology of Precursor Lesions.

Biliary tract cancer is a devastating malignancy of the bile ducts and gallbladder with a dismal prognosis. The study of precancerous lesions has received considerable attention and led to a histopathological classification which, in some respects, remains an evolving field. Consequently, increasing efforts have been devoted to characterizing the molecular pathogenesis of the precursor lesions, with the aim of better understanding the mechanisms of tumor progression, and with the ultimate goal of meeting the challenges of early diagnosis and treatment. This review delves into the molecular mechanisms that initiate and promote the development of precursor lesions of intra- and extrahepatic cholangiocarcinoma and of gallbladder carcinoma. It addresses the genomic, epigenomic, and transcriptomic landscape of these precursors and provides an overview of animal and organoid models used to study them. In conclusion, this review summarizes the known molecular features of precancerous lesions in biliary tract cancer and highlights our fragmentary knowledge of the molecular pathogenesis of tumor initiation.

Axon guidance genes control hepatic artery development.

Earlier data on liver development demonstrated that morphogenesis of the bile duct, portal mesenchyme and hepatic artery is interdependent, yet how this interdependency is orchestrated remains unknown. Here, using 2D and 3D imaging, we first describe how portal mesenchymal cells become organised to form hepatic arteries. Next, we examined intercellular signalling active during portal area development and found that axon guidance genes are dynamically expressed in developing bile ducts and portal mesenchyme. Using tissue-specific gene inactivation in mice, we show that the repulsive guidance molecule BMP co-receptor A (RGMA)/neogenin (NEO1) receptor/ligand pair is dispensable for portal area development, but that deficient roundabout 2 (ROBO2)/SLIT2 signalling in the portal mesenchyme causes reduced maturation of the vascular smooth muscle cells that form the tunica media of the hepatic artery. This arterial anomaly does not impact liver function in homeostatic conditions, but is associated with significant tissular damage following partial hepatectomy. In conclusion, our work identifies new players in development of the liver vasculature in health and liver regeneration.

KRAS protein expression becomes progressively restricted during embryogenesis and in adulthood.

KRAS mutants are common in many cancers and wild-type KRAS is essential in development as its absence causes embryonic lethality. Despite this critical role in development and disease, the normal expression pattern of KRAS protein is still largely unknown at the tissue level due to the lack of valid antibodies. To address this issue, we used the citrine-Kras mouse model in which the Citrine-KRAS (Cit-K) fusion protein functions as a validated surrogate of endogenous KRAS protein that can be detected on tissue sections by immunolabeling with a GFP antibody. In the embryo, we found expression of KRAS protein in a wide range of organs and tissues. This expression tends to decrease near birth, mainly in mesenchymal cells. During transition to the adult stage, the dynamics of KRAS protein expression vary among organs and detection of KRAS becomes restricted to specific cell types. Furthermore, we found that steady state KRAS protein expression is detectable at the cell membrane and in the cytoplasm and that this subcellular partitioning differed among cell types. Our results reveal hitherto unanticipated dynamics in developmental, tissular, cell-specific and subcellular expression of KRAS protein. They provide insight into the reason why specific cell-types are sensitive to KRAS mutations during cancer initiation.

Quantitative modeling identifies critical cell mechanics driving bile duct lumen formation.

Biliary ducts collect bile from liver lobules, the smallest functional and anatomical units of liver, and carry it to the gallbladder. Disruptions in this process caused by defective embryonic development, or through ductal reaction in liver disease have a major impact on life quality and survival of patients. A deep understanding of the processes underlying bile duct lumen formation is crucial to identify intervention points to avoid or treat the appearance of defective bile ducts. Several hypotheses have been proposed to characterize the biophysical mechanisms driving initial bile duct lumen formation during embryogenesis. Here, guided by the quantification of morphological features and expression of genes in bile ducts from embryonic mouse liver, we sharpened these hypotheses and collected data to develop a high resolution individual cell-based computational model that enables to test alternative hypotheses in silico. This model permits realistic simulations of tissue and cell mechanics at sub-cellular scale. Our simulations suggest that successful bile duct lumen formation requires a simultaneous contribution of directed cell division of cholangiocytes, local osmotic effects generated by salt excretion in the lumen, and temporally-controlled differentiation of hepatoblasts to cholangiocytes, with apical constriction of cholangiocytes only moderately affecting luminal size.

Notch-IGF1 signaling during liver regeneration drives biliary epithelial cell expansion and inhibits hepatocyte differentiation.

In the adult liver, a population of facultative progenitor cells called biliary epithelial cells (BECs) proliferate and differentiate into cholangiocytes and hepatocytes after injury, thereby restoring liver function. In mammalian models of chronic liver injury, Notch signaling is essential for bile duct formation from these cells. However, the continual proliferation of BECs and differentiation of hepatocytes in these models have limited their use for determining whether Notch signaling is required for BECs to replenish hepatocytes after injury in the mammalian liver. Here, we used a temporally restricted model of hepatic repair in which large-scale hepatocyte injury and regeneration are initiated through the acute loss of Mdm2 in hepatocytes, resulting in the rapid, coordinated proliferation of BECs. We found that transient, early activation of Notch1- and Notch3-mediated signaling and entrance into the cell cycle preceded the phenotypic expansion of BECs into hepatocytes. Notch inhibition reduced BEC proliferation, which resulted in failure of BECs to differentiate into hepatocytes, indicating that Notch-dependent expansion of BECs is essential for hepatocyte regeneration. Notch signaling increased the abundance of the insulin-like growth factor 1 receptor (IGF1R) in BECs, and activating IGFR signaling increased BEC numbers but suppressed BEC differentiation into hepatocytes. These results suggest that different signaling mechanisms control BEC expansion and hepatocyte differentiation.

Building consensus on definition and nomenclature of hepatic, pancreatic, and biliary organoids.

Hepatic, pancreatic, and biliary (HPB) organoids are powerful tools for studying development, disease, and regeneration. As organoid research expands, the need for clear definitions and nomenclature describing these systems also grows. To facilitate scientific communication and consistent interpretation, we revisit the concept of an organoid and introduce an intuitive classification system and nomenclature for describing these 3D structures through the consensus of experts in the field. To promote the standardization and validation of HPB organoids, we propose guidelines for establishing, characterizing, and benchmarking future systems. Finally, we address some of the major challenges to the clinical application of organoids.

Temporal dynamics of a CSF1R signaling gene regulatory network involved in epilepsy.

Colony Stimulating Factor 1 Receptor (CSF1R) is a potential target for anti-epileptic drugs. However, inhibition of CSF1R is not well tolerated by patients, thereby prompting the need for alternative targets. To develop a framework for identification of such alternatives, we here develop a mathematical model of a pro-inflammatory gene regulatory network (GRN) involved in epilepsy and centered around CSF1R. This GRN comprises validated transcriptional and post-transcriptional regulations involving STAT1, STAT3, NFκB, IL6R, CSF3R, IRF8, PU1, C/EBPα, TNFR1, CSF1 and CSF1R. The model was calibrated on mRNA levels of all GRN components in lipopolysaccharide (LPS)-treated mouse microglial BV-2 cells, and allowed to predict that STAT1 and STAT3 have the strongest impact on the expression of the other GRN components. Microglial BV-2 cells were selected because, the modules from which the GRN was deduced are enriched for microglial marker genes. The function of STAT1 and STAT3 in the GRN was experimentally validated in BV-2 cells. Further, in silico analysis of the GRN dynamics predicted that a pro-inflammatory stimulus can induce irreversible bistability whereby the expression level of GRN components occurs as two distinct states. The irreversibility of the switch may enforce the need for chronic inhibition of the CSF1R GRN in order to achieve therapeutic benefit. The cell-to-cell heterogeneity driven by the bistability may cause variable therapeutic response. In conclusion, our modeling approach uncovered a GRN controlling CSF1R that is predominantly regulated by STAT1 and STAT3. Irreversible inflammation-induced bistability and cell-to-cell heterogeneity of the GRN provide a theoretical foundation to the need for chronic GRN control and the limited potential for disease modification via inhibition of CSF1R.

A Mouse Model of Cholangiocarcinoma Uncovers a Role for Tensin-4 in Tumor Progression.

BACKGROUND AND AIMS: Earlier diagnosis and treatment of intrahepatic cholangiocarcinoma (iCCA) are necessary to improve therapy, yet limited information is available about initiation and evolution of iCCA precursor lesions. Therefore, there is a need to identify mechanisms driving formation of precancerous lesions and their progression toward invasive tumors using experimental models that faithfully recapitulate human tumorigenesis. APPROACH AND RESULTS: To this end, we generated a mouse model which combines cholangiocyte-specific expression of KrasG12D with 3,5-diethoxycarbonyl-1,4-dihydrocollidine (DDC) diet-induced inflammation to mimic iCCA development in patients with cholangitis. Histological and transcriptomic analyses of the mouse precursor lesions and iCCA were performed and compared with human analyses. The function of genes overexpressed during tumorigenesis was investigated in human cell lines. We found that mice expressing KrasG12D in cholangiocytes and fed a DDC diet developed cholangitis, ductular proliferations, intraductal papillary neoplasms of bile ducts (IPNBs), and, eventually, iCCAs. The histology of mouse and human IPNBs was similar, and mouse iCCAs displayed histological characteristics of human mucin-producing, large-duct-type iCCA. Signaling pathways activated in human iCCA were also activated in mice. The identification of transition zones between IPNB and iCCA on tissue sections, combined with RNA-sequencing analyses of the lesions supported that iCCAs derive from IPNBs. We further provide evidence that tensin-4 (TNS4), which is stimulated by KRASG12D and SRY-related HMG box transcription factor 17, promotes tumor progression. CONCLUSIONS: We developed a mouse model that faithfully recapitulates human iCCA tumorigenesis and identified a gene cascade which involves TNS4 and promotes tumor progression.

Dynamic Regulation of Expression of KRAS and Its Effectors Determines the Ability to Initiate Tumorigenesis in Pancreatic Acinar Cells.

Pancreatic acinar cells are a cell type of origin for pancreatic cancer that become progressively less sensitive to tumorigenesis induced by oncogenic Kras mutations after birth. This sensitivity is increased when Kras mutations are combined with pancreatitis. Molecular mechanisms underlying these observations are still largely unknown. To identify these mechanisms, we generated the first CRISPR-edited mouse models that enable detection of wild-type and mutant KRAS proteins in vivo. Analysis of these mouse models revealed that more than 75% of adult acinar cells are devoid of detectable KRAS protein. In the 25% of acinar cells expressing KRAS protein, transcriptomic analysis highlighted a slight upregulation of the RAS and MAPK pathways. However, at the protein level, only marginal pancreatic expression of essential KRAS effectors, including C-RAF, was observed. The expression of KRAS and its effectors gradually decreased after birth. The low sensitivity of adult acinar cells to Kras mutations resulted from low expression of KRAS and its effectors and the subsequent lack of activation of RAS/MAPK pathways. Pancreatitis triggered expression of KRAS and its effectors as well as subsequent activation of downstream signaling; this induction required the activity of EGFR. Finally, expression of C-RAF in adult pancreas was required for pancreatic tumorigenesis. In conclusion, our study reveals that control of the expression of KRAS and its effectors regulates the sensitivity of acinar cells to transformation by oncogenic Kras mutations. SIGNIFICANCE: This study generates new mouse models to study regulation of KRAS during pancreatic tumorigenesis and highlights a novel mechanism through which pancreatitis sensitizes acinar cells to Kras mutations.

Regeneration Defects in Yap and Taz Mutant Mouse Livers Are Caused by Bile Duct Disruption and Cholestasis.

BACKGROUND AND AIMS: The Hippo pathway and its downstream effectors YAP and TAZ (YAP/TAZ) are heralded as important regulators of organ growth and regeneration. However, different studies provided contradictory conclusions about their role during regeneration of different organs, ranging from promoting proliferation to inhibiting it. Here we resolve the function of YAP/TAZ during regeneration of the liver, where Hippo's role in growth control has been studied most intensely. METHODS: We evaluated liver regeneration after carbon tetrachloride toxic liver injury in mice with conditional deletion of Yap/Taz in hepatocytes and/or biliary epithelial cells, and measured the behavior of different cell types during regeneration by histology, RNA sequencing, and flow cytometry. RESULTS: We found that YAP/TAZ were activated in hepatocytes in response to carbon tetrachloride toxic injury. However, their targeted deletion in adult hepatocytes did not noticeably impair liver regeneration. In contrast, Yap/Taz deletion in adult bile ducts caused severe defects and delay in liver regeneration. Mechanistically, we showed that Yap/Taz mutant bile ducts degenerated, causing cholestasis, which stalled the recruitment of phagocytic macrophages and the removal of cellular corpses from injury sites. Elevated bile acids activated pregnane X receptor, which was sufficient to recapitulate the phenotype observed in mutant mice. CONCLUSIONS: Our data show that YAP/TAZ are practically dispensable in hepatocytes for liver development and regeneration. Rather, YAP/TAZ play an indirect role in liver regeneration by preserving bile duct integrity and securing immune cell recruitment and function.

Differential impact of the ERBB receptors EGFR and ERBB2 on the initiation of precursor lesions of pancreatic ductal adenocarcinoma.

Earlier diagnosis of pancreatic ductal adenocarcinoma (PDAC) requires better understanding of the mechanisms driving tumorigenesis. In this context, depletion of Epidermal Growth Factor Receptor (EGFR) is known to impair development of PDAC-initiating lesions called acinar-to-ductal metaplasia (ADM) and Pancreatic Intraepithelial Neoplasia (PanIN). In contrast, the role of v-erb-b2 erythroblastic leukemia viral oncogene homolog 2 (ERBB2), the preferred dimerization partner of EGFR, remains poorly understood. Here, using a mouse model with inactivation of Erbb2 in pancreatic acinar cells, we found that Erbb2 is dispensable for inflammation- and KRasG12D-induced development of ADM and PanIN. A mathematical model of EGFR/ERBB2-KRAS signaling, which was calibrated on mouse and human data, supported the observed roles of EGFR and ERBB2. However, this model also predicted that overexpression of ERBB2 stimulates ERBB/KRAS signaling; this prediction was validated experimentally. We conclude that EGFR and ERBB2 differentially impact ERBB signaling during PDAC tumorigenesis, and that the oncogenic potential of ERBB2 is only manifested when it is overexpressed. Therefore, the level of ERBB2, not only its mere presence, needs to be considered when designing therapies targeting ERBB signaling.

Development of the Intrahepatic and Extrahepatic Biliary Tract: A Framework for Understanding Congenital Diseases.

The involvement of the biliary tract in the pathophysiology of liver diseases and the increased attention paid to bile ducts in the bioconstruction of liver tissue for regenerative therapy have fueled intense research into the fundamental mechanisms of biliary development. Here, I review the molecular, cellular and tissular mechanisms driving differentiation and morphogenesis of the intrahepatic and extrahepatic bile ducts. This review focuses on the dynamics of the transcriptional and signaling modules that promote biliary development in human and mouse liver and discusses studies in which the use of zebrafish uncovered unexplored processes in mammalian biliary development. The review concludes by providing a framework for interpreting the mechanisms that may help us understand the origin of congenital biliary diseases.

Kras and Lkb1 mutations synergistically induce intraductal papillary mucinous neoplasm derived from pancreatic duct cells.

OBJECTIVE: Pancreatic cancer can arise from precursor lesions called intraductal papillary mucinous neoplasms (IPMN), which are characterised by cysts containing papillae and mucus-producing cells. The high frequency of KRAS mutations in IPMN and histological analyses suggest that oncogenic KRAS drives IPMN development from pancreatic duct cells. However, induction of Kras mutation in ductal cells is not sufficient to generate IPMN, and formal proof of a ductal origin of IPMN is still missing. Here we explore whether combining oncogenic KrasG12D mutation with an additional gene mutation known to occur in human IPMN can induce IPMN from pancreatic duct cells. DESIGN: We created and phenotyped mouse models in which mutations in Kras and in the tumour suppressor gene liver kinase B1 (Lkb1/Stk11) are conditionally induced in pancreatic ducts using Cre-mediated gene recombination. We also tested the effect of ß-catenin inhibition during formation of the lesions. RESULTS: Activating KrasG12D mutation and Lkb1 inactivation synergised to induce IPMN, mainly of gastric type and with malignant potential. The mouse lesions shared several features with human IPMN. Time course analysis suggested that IPMN developed from intraductal papillae and glandular neoplasms, which both derived from the epithelium lining large pancreatic ducts. ß-catenin was required for the development of glandular neoplasms and subsequent development of the mucinous cells in IPMN. Instead, the lack of ß-catenin did not impede formation of intraductal papillae and their progression to papillary lesions in IPMN. CONCLUSION: Our work demonstrates that IPMN can result from synergy between KrasG12D mutation and inactivation of a tumour suppressor gene. The ductal epithelium can give rise to glandular neoplasms and papillary lesions, which probably both contribute to IPMN formation.

Genetically engineered animal models of biliary tract cancers.

PURPOSE OF REVIEW: Biliary tract cancers which include intrahepatic and extrahepatic cholangiocarcinomas and gallbladder cancer, are characterized by poor outcome. Therefore, identifying the molecular mechanisms of the disease has become a priority. However, such identification has to cope with extreme heterogeneity of the disease, which results from the variable anatomical location, the numerous cell types of origin and the high number of known genetic alterations. RECENT FINDINGS: Animal models can develop invasive and metastatic tumours that recapitulate as faithfully as possible the molecular features of the human tumours. To generate animal models of cholangiocarcinoma, investigators resorted to the administration of carcinogens, induction of cholestasis, grafting of tumour cells and induction of genetic modifications. SUMMARY: Here, we summarize the currently available genetically engineered animal models, and focus on mice and zebrafish. The experimental strategies that were selected to induce cholangiocarcinoma in a time-controlled and cell-type-specific manner are critically examined. We discuss their strengths and limitations while considering their relevance to human pathophysiology.

Modeling the Dynamics of Let-7-Coupled Gene Regulatory Networks Linking Cell Proliferation to Malignant Transformation.

Let-7 microRNA controls the expression of proteins that belong to two distinct gene regulatory networks, namely, a cyclin-dependent kinase (Cdk) network driving the cell cycle and a cell transformation network that can undergo an epigenetic switch between a non-transformed and a malignant transformed cell state. Using mathematical modeling and transcriptomic data analysis, we here investigate how Let-7 controls the Cdk-dependent cell cycle network, and how it couples the latter with the transformation network. We also assess the consequence of this coupling on cancer progression. Our analysis shows that the switch from a quiescent to a proliferative state depends on the relative levels of Let-7 and several cell cycle activators. Numerical simulations further indicate that the Let-7-coupled cell cycle and transformation networks mutually control each other, and our model identifies key players for this mutual control. Transcriptomic data analysis from The Cancer Genome Atlas (TCGA) suggests that the two networks are activated in cancer, in particular in gastrointestinal cancers, and that the activation levels vary significantly among patients affected by a same cancer type. Our mathematical model, when applied to a heterogeneous cell population, suggests that heterogeneity among tumors may in part result from stochastic switches between a non-transformed cell state with low proliferative capability and a transformed cell state with high proliferative property. The model further predicts that Let-7 may reduce tumor heterogeneity by decreasing the occurrence of stochastic switches toward a transformed, proliferative cell state. In conclusion, we identified the key components responsible for the qualitative dynamics of two networks interconnected by Let-7. The two networks are heterogeneously activated in several cancers, thereby stressing the need to consider patient's specific characteristics to optimize therapeutic strategies.

Dynamics and predicted drug response of a gene network linking dedifferentiation with beta-catenin dysfunction in hepatocellular carcinoma.

BACKGROUND & AIMS: Alterations of individual genes variably affect the development of hepatocellular carcinoma (HCC). Thus, we aimed to characterize the function of tumor-promoting genes in the context of gene regulatory networks (GRNs). METHODS: Using data from The Cancer Genome Atlas, from the LIRI-JP (Liver Cancer - RIKEN, JP project), and from our transcriptomic, transfection and mouse transgenic experiments, we identify a GRN which functionally links LIN28B-dependent dedifferentiation with dysfunction of ß-catenin (CTNNB1). We further generated and validated a quantitative mathematical model of the GRN using human cell lines and in vivo expression data. RESULTS: We found that LIN28B and CTNNB1 form a GRN with SMARCA4, Let-7b (MIRLET7B), SOX9, TP53 and MYC. GRN functionality is detected in HCC and gastrointestinal cancers, but not in other cancer types. GRN status negatively correlates with HCC prognosis, and positively correlates with hyperproliferation, dedifferentiation and HGF/MET pathway activation, suggesting that it contributes to a transcriptomic profile typical of the proliferative class of HCC. The mathematical model predicts how the expression of GRN components changes when the expression of another GRN member varies or is inhibited by a pharmacological drug. The dynamics of GRN component expression reveal distinct cell states that can switch reversibly in normal conditions, and irreversibly in HCC. The mathematical model is available via a web-based tool which can evaluate the GRN status of HCC samples and predict the impact of therapeutic agents on the GRN. CONCLUSIONS: We conclude that identification and modelling of the GRN provide insights into the prognosis of HCC and the mechanisms by which tumor-promoting genes impact on HCC development. LAY SUMMARY: Hepatocellular carcinoma (HCC) is a heterogeneous disease driven by the concomitant deregulation of several genes functionally organized as networks. Here, we identified a gene regulatory network involved in a subset of HCCs. This subset is characterized by increased proliferation and poor prognosis. We developed a mathematical model which uncovers the dynamics of the network and allows us to predict the impact of a therapeutic agent, not only on its specific target but on all the genes belonging to the network.