Quantitative lineage analysis identifies a hepato-pancreato-biliary progenitor niche.

Studies based on single cells have revealed vast cellular heterogeneity in stem cell and progenitor compartments, suggesting continuous differentiation trajectories with intermixing of cells at various states of lineage commitment and notable degrees of plasticity during organogenesis1-5. The hepato-pancreato-biliary organ system relies on a small endoderm progenitor compartment that gives rise to a variety of different adult tissues, including the liver, pancreas, gall bladder and extra-hepatic bile ducts6,7. Experimental manipulation of various developmental signals in the mouse embryo has underscored important cellular plasticity in this embryonic territory6. This is reflected in the existence of human genetic syndromes as well as congenital malformations featuring multi-organ phenotypes in liver, pancreas and gall bladder6. Nevertheless, the precise lineage hierarchy and succession of events leading to the segregation of an endoderm progenitor compartment into hepatic, biliary and pancreatic structures have not yet been established. Here we combine computational modelling approaches with genetic lineage tracing to accurately reconstruct the hepato-pancreato-biliary lineage tree. We show that a multipotent progenitor subpopulation persists in the pancreato-biliary organ rudiment, contributing cells not only to the pancreas and gall bladder but also to the liver. Moreover, using single-cell RNA sequencing and functional experiments we define a specialized niche that supports this subpopulation in a multipotent state for an extended time during development. Together these findings indicate sustained plasticity underlying hepato-pancreato-biliary development that might also explain the rapid expansion of the liver while attenuating pancreato-biliary growth.

A Quantitative Modular Modeling Approach Reveals the Effects of Different A20 Feedback Implementations for the NF-kB Signaling Dynamics.

Signaling pathways involve complex molecular interactions and are controled by non-linear regulatory mechanisms. If details of regulatory mechanisms are not fully elucidated, they can be implemented by different, equally reasonable mathematical representations in computational models. The study presented here focusses on NF-κB signaling, which is regulated by negative feedbacks via IκBα and A20. A20 inhibits NF-κB activation indirectly through interference with proteins that transduce the signal from the TNF receptor complex to activate the IκB kinase (IKK) complex. A number of pathway models has been developed implementing the A20 effect in different ways. We here focus on the question how different A20 feedback implementations impact the dynamics of NF-κB. To this end, we develop a modular modeling approach that allows combining previously published A20 modules with a common pathway core module. The resulting models are fitted to a published comprehensive experimental data set and therefore show quantitatively comparable NF-κB dynamics. Based on defined measures for the initial and long-term behavior we analyze the effects of a wide range of changes in the A20 feedback strength, the IκBα feedback strength and the TNFα stimulation strength on NF-κB dynamics. This shows similarities between the models but also model-specific differences. In particular, the A20 feedback strength and the TNFα stimulation strength affect initial and long-term NF-κB concentrations differently in the analyzed models. We validated our model predictions experimentally by varying TNFα concentrations applied to HeLa cells. These time course data indicate that only one of the A20 feedback models appropriately describes the impact of A20 on the NF-κB dynamics in this cell type.

N-Myc-induced metabolic rewiring creates novel therapeutic vulnerabilities in neuroblastoma.

N-Myc is a transcription factor that is aberrantly expressed in many tumor types and is often correlated with poor patient prognosis. Recently, several lines of evidence pointed to the fact that oncogenic activation of Myc family proteins is concomitant with reprogramming of tumor cells to cope with an enhanced need for metabolites during cell growth. These adaptions are driven by the ability of Myc proteins to act as transcriptional amplifiers in a tissue-of-origin specific manner. Here, we describe the effects of N-Myc overexpression on metabolic reprogramming in neuroblastoma cells. Ectopic expression of N-Myc induced a glycolytic switch that was concomitant with enhanced sensitivity towards 2-deoxyglucose, an inhibitor of glycolysis. Moreover, global metabolic profiling revealed extensive alterations in the cellular metabolome resulting from overexpression of N-Myc. Limited supply with either of the two main carbon sources, glucose or glutamine, resulted in distinct shifts in steady-state metabolite levels and significant changes in glutathione metabolism. Interestingly, interference with glutamine-glutamate conversion preferentially blocked proliferation of N-Myc overexpressing cells, when glutamine levels were reduced. Thus, our study uncovered N-Myc induction and nutrient levels as important metabolic master switches in neuroblastoma cells and identified critical nodes that restrict tumor cell proliferation.

Controlling Nuclear NF-κB Dynamics by ß-TrCP-Insights from a Computational Model.

The canonical nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) signaling pathway regulates central processes in mammalian cells and plays a fundamental role in the regulation of inflammation and immunity. Aberrant regulation of the activation of the transcription factor NF-κB is associated with severe diseases such as inflammatory bowel disease and arthritis. In the canonical pathway, the inhibitor IκB suppresses NF-κB's transcriptional activity. NF-κB becomes active upon the degradation of IκB, a process that is, in turn, regulated by the ß-transducin repeat-containing protein (ß-TrCP). ß-TrCP has therefore been proposed as a promising pharmacological target in the development of novel therapeutic approaches to control NF-κB's activity in diseases. This study explores the extent to which ß-TrCP affects the dynamics of nuclear NF-κB using a computational model of canonical NF-κB signaling. The analysis predicts that ß-TrCP influences the steady-state concentration of nuclear NF-κB, as well as changes characteristic dynamic properties of nuclear NF-κB, such as fold-change and the duration of its response to pathway stimulation. The results suggest that the modulation of ß-TrCP has a high potential to regulate the transcriptional activity of NF-κB.

Paracrine and autocrine regulation of gene expression by Wnt-inhibitor Dickkopf in wild-type and mutant hepatocytes.

BACKGROUND: Cells are able to communicate and coordinate their function within tissues via secreted factors. Aberrant secretion by cancer cells can modulate this intercellular communication, in particular in highly organised tissues such as the liver. Hepatocytes, the major cell type of the liver, secrete Dickkopf (Dkk), which inhibits Wnt/ ß-catenin signalling in an autocrine and paracrine manner. Consequently, Dkk modulates the expression of Wnt/ ß-catenin target genes. We present a mathematical model that describes the autocrine and paracrine regulation of hepatic gene expression by Dkk under wild-type conditions as well as in the presence of mutant cells. RESULTS: Our spatial model describes the competition of Dkk and Wnt at receptor level, intra-cellular Wnt/ ß-catenin signalling, and the regulation of target gene expression for 21 individual hepatocytes. Autocrine and paracrine regulation is mediated through a feedback mechanism via Dkk and Dkk diffusion along the porto-central axis. Along this axis an APC concentration gradient is modelled as experimentally detected in liver. Simulations of mutant cells demonstrate that already a single mutant cell increases overall Dkk concentration. The influence of the mutant cell on gene expression of surrounding wild-type hepatocytes is limited in magnitude and restricted to hepatocytes in close proximity. To explore the underlying molecular mechanisms, we perform a comprehensive analysis of the model parameters such as diffusion coefficient, mutation strength and feedback strength. CONCLUSIONS: Our simulations show that Dkk concentration is elevated in the presence of a mutant cell. However, the impact of these elevated Dkk levels on wild-type hepatocytes is confined in space and magnitude. The combination of inter- and intracellular processes, such as Dkk feedback, diffusion and Wnt/ ß-catenin signal transduction, allow wild-type hepatocytes to largely maintain their gene expression.

Different promoter affinities account for specificity in MYC-dependent gene regulation.

Enhanced expression of the MYC transcription factor is observed in the majority of tumors. Two seemingly conflicting models have been proposed for its function: one proposes that MYC enhances expression of all genes, while the other model suggests gene-specific regulation. Here, we have explored the hypothesis that specific gene expression profiles arise since promoters differ in affinity for MYC and high-affinity promoters are fully occupied by physiological levels of MYC. We determined cellular MYC levels and used RNA- and ChIP-sequencing to correlate promoter occupancy with gene expression at different concentrations of MYC. Mathematical modeling showed that binding affinities for interactions of MYC with DNA and with core promoter-bound factors, such as WDR5, are sufficient to explain promoter occupancies observed in vivo. Importantly, promoter affinity stratifies different biological processes that are regulated by MYC, explaining why tumor-specific MYC levels induce specific gene expression programs and alter defined biological properties of cells.

A Multiplexed NMR-Reporter Approach to Measure Cellular Kinase and Phosphatase Activities in Real-Time.

Cell signaling is governed by dynamic changes in kinase and phosphatase activities, which are difficult to assess with discontinuous readout methods. Here, we introduce an NMR-based reporter approach to directly identify active kinases and phosphatases in complex physiological environments such as cell lysates and to measure their individual activities in a semicontinuous fashion. Multiplexed NMR profiling of reporter phosphorylation states provides unique advantages for kinase inhibitor studies and reveals reversible modulations of cellular enzyme activities under different metabolic conditions.

Mathematical modelling suggests a differential impact of ß-transducin repeat-containing protein paralogues on Wnt/ß-catenin signalling dynamics.

The Wnt/ß-catenin signalling pathway is involved in the regulation of a multitude of cellular processes by controlling the concentration of the transcriptional regulator ß-catenin. Proteasomal degradation of ß-catenin is mediated by two ß-transducin repeat-containing protein paralogues, homologous to Slimb protein (HOS) and F-box/WD repeat-containing protein 1A (FWD1), which are functionally interchangeable and thereby considered to function redundantly in the pathway. HOS and FWD1 are both regulated by Wnt/ß-catenin signalling, albeit in opposite directions, thus establishing interlocked negative and positive feedback loops. The functional relevance of the opposite regulation of HOS and FWD1 by Wnt/ß-catenin signalling in conjunction with their redundant activities in proteasomal degradation of ß-catenin remains unresolved. Using a detailed ordinary differential equation model, we investigated the specific influence of each individual feedback mechanism and their combination on Wnt/ß-catenin signal transduction under wild-type and cancerous conditions. We found that, under wild-type conditions, the signalling dynamics are predominantly affected by the HOS feedback as a result of a higher concentration of HOS than FWD1. Transcriptional up-regulation of FWD1 by other signalling pathways reduced the impact of the HOS feedback. The opposite regulation of HOS and FWD1 expression by Wnt/ß-catenin signalling allows the FWD1 feedback to be employed as a compensation mechanism against aberrant pathway activation as a result of a reduced HOS concentration. By contrast, the FWD1 feedback provides no protection against aberrant activation in adenomatous polyposis coli protein mutant cancer cells.

Modeling Wnt/ß-Catenin Target Gene Expression in APC and Wnt Gradients Under Wild Type and Mutant Conditions.

The Wnt/ß-catenin pathway is involved in the regulation of a multitude of physiological processes by controlling the differential expression of target genes. In certain tissues such as the adult liver, the Wnt/ß-catenin pathway can attain different levels of activity due to gradients of Wnt ligands and/or intracellular pathway components like APC. How graded pathway activity is converted into regionally distinct patterns of Wnt/ß-catenin target gene expression is largely unknown. Here, we apply a mathematical modeling approach to investigate the impact of different regulatory mechanisms on target gene expression within Wnt or APC concentration gradients. We develop a minimal model of Wnt/ß-catenin signal transduction and combine it with various mechanisms of target gene regulation. In particular, the effects of activation, inhibition, and an incoherent feedforward loop (iFFL) are compared. To specify activation kinetics, we analyze experimental data that quantify the response of ß-catenin/TCF reporter constructs to different Wnt concentrations, and demonstrate that the induction of these constructs occurs in a cooperative manner with Hill coefficients between 2 and 5. In summary, our study shows that the combination of specific gene regulatory mechanisms with a time-independent gradient of Wnt or APC is sufficient to generate distinct target gene expression patterns as have been experimentally observed in liver. We find that cooperative gene activation in combination with a TCF feedback can establish sharp borders of target gene expression in Wnt or APC gradients. In contrast, the iFFL renders gene expression independent of gradients of the upstream signaling components. Our subsequent analysis of carcinogenic pathway mutations reveals that their impact on gene expression is determined by the gene regulatory mechanism and the APC concentration of the cell in which the mutation occurs.

Digital NFATc2 activation per cell transforms graded T cell receptor activation into an all-or-none IL-2 expression.

The expression of interleukin-2 (IL-2) is a key event in T helper (Th) lymphocyte activation, controlling both, the expansion and differentiation of effector Th cells as well as the activation of regulatory T cells. We demonstrate that the strength of TCR stimulation is translated into the frequency of memory Th cells expressing IL-2 but not into the amount of IL-2 per cell. This molecular switch decision for IL-2 expression per cell is located downstream of the cytosolic Ca2+ level. Here we show that in a single activated Th cell, NFATc2 activation is digital but NF-kappaB activation is graded after graded T cell receptor (TCR) signaling. Subsequently, NFATc2 translocates into the nucleus in an all-or-none fashion per cell, transforming the strength of TCR-stimulation into the number of nuclei positive for NFATc2 and IL-2 transcription. Thus, the described NFATc2 switch regulates the number of Th cells actively participating in an immune response.