The Number of Transcription Factors at an Enhancer Determines Switch-like Gene Expression

Roles of Post-Translational Modifications of Transcription Factors Involved in Breast Cancer Hypoxia

BC is the most commonly diagnosed cancer and the second leading cause of cancer death among women worldwide. Cellular stress is a condition that leads to disrupted homeostasis by extrinsic and intrinsic factors. Among other stressors, hypoxia is a driving force for breast cancer (BC) progression and a general hallmark of solid tumors. Thus, intratumoral hypoxia is an important determinant of invasion, metastasis, treatment failure, prognosis, and patient mortality. Acquisition of the epithelial-mesenchymal transition (EMT) phenotype is also a consequence of tumor hypoxia. The cellular response to hypoxia is mainly regulated by the hypoxia signaling pathway, governed by hypoxia-inducible factors (HIFs), mainly HIF1α. HIFs are a family of transcription factors (TFs), which induce the expression of target genes involved in cell survival and proliferation, metabolic reprogramming, angiogenesis, resisting apoptosis, invasion, and metastasis. HIF1α cooperates with a large number of other TFs. In this review, we focused on the crosstalk and cooperation between HIF1α and other TFs involved in the cellular response to hypoxia in BC. We identified a cluster of TFs, proposed as the HIF1α-TF interactome, that orchestrates the transcription of target genes involved in hypoxia, due to their post-translational modifications (PTMs), including phosphorylation/dephosphorylation, ubiquitination/deubiquitination, SUMOylation, hydroxylation, acetylation, S-nitrosylation, and palmitoylation. PTMs of these HIF1α-related TFs drive their stability and activity, degradation and turnover, and the bidirectional translocation between the cytoplasm or plasma membrane and nucleus of BC cells, as well as the transcription/activation of proteins encoded by oncogenes or inactivation of tumor suppressor target genes. Consequently, PTMs of TFs in the HIF1α interactome are crucial regulatory mechanisms that drive the cellular response to oxygen deprivation in BC cells.

NFκB dynamics-dependent epigenetic changes modulate inflammatory gene expression and induce cellular senescence

Upregulation of nuclear factor κB (NFκB) signaling is a hallmark of aging and a major cause of age-related chronic inflammation. However, its effect on cellular senescence remains unclear. Here, we show that alteration of NFκB nuclear dynamics from oscillatory to sustained by depleting a negative feedback regulator of NFκB pathway, NFκB inhibitor alpha (IκBα), in the presence of tumor necrosis factor α (TNFα) promotes cellular senescence. Sustained NFκB activity enhanced inflammatory gene expression through increased NFκB-DNA binding and slowed the cell cycle. IκBα protein was decreased under replicative or oxidative stress in vitro. Furthermore, a decrease in IκBα protein and an increase in DNA-NFκB binding at the transcription start sites of age-associated genes in aged mouse hearts suggested that nuclear NFκB dynamics may play a critical role in the progression of aging. Our study suggests that nuclear NFκB dynamics-dependent epigenetic changes regulated over time in a living system, possibly through a decrease in IκBα, enhance the expression of inflammatory genes to advance the cells to a senescent state.

TWEAK/Fn14 signalling driven super-enhancer reprogramming promotes pro-metastatic metabolic rewiring in triple-negative breast cancer

Triple Negative Breast Cancer (TNBC) is the most aggressive breast cancer subtype suffering from limited targeted treatment options. Following recent reports correlating Fibroblast growth factor-inducible 14 (Fn14) receptor overexpression in Estrogen Receptor (ER)-negative breast cancers with metastatic events, we show that Fn14 is specifically overexpressed in TNBC patients and associated with poor survival. We demonstrate that constitutive Fn14 signalling rewires the transcriptomic and epigenomic landscape of TNBC, leading to enhanced tumour growth and metastasis. We further illustrate that such mechanisms activate TNBC-specific super enhancers (SE) to drive the transcriptional activation of cancer dependency genes via chromatin looping. In particular, we uncover the SE-driven upregulation of Nicotinamide phosphoribosyltransferase (NAMPT), which promotes NAD+ and ATP metabolic reprogramming critical for filopodia formation and metastasis. Collectively, our study details the complex mechanistic link between TWEAK/Fn14 signalling and TNBC metastasis, which reveals several vulnerabilities which could be pursued for the targeted treatment of TNBC patients.

A multiscale spatial modeling framework for the germinal center response

The germinal center response or reaction (GCR) is a hallmark event of adaptive humoral immunity. Unfolding in the B cell follicles of the secondary lymphoid organs, a GC culminates in the production of high-affinity antibody-secreting plasma cells along with memory B cells. By interacting with follicular dendritic cells (FDC) and T follicular helper (Tfh) cells, GC B cells exhibit complex spatiotemporal dynamics. Driving the B cell dynamics are the intracellular signal transduction and gene regulatory network that responds to cell surface signaling molecules, cytokines, and chemokines. As our knowledge of the GC continues to expand in depth and in scope, mathematical modeling has become an important tool to help disentangle the intricacy of the GCR and inform novel mechanistic and clinical insights. While the GC has been modeled at different granularities, a multiscale spatial simulation framework - integrating molecular, cellular, and tissue-level responses - is still rare. Here, we report our recent progress toward this end with a hybrid stochastic GC framework developed on the Cellular Potts Model-based CompuCell3D platform. Tellurium is used to simulate the B cell intracellular molecular network comprising NF-κB, FOXO1, MYC, AP4, CXCR4, and BLIMP1 that responds to B cell receptor (BCR) and CD40-mediated signaling. The molecular outputs of the network drive the spatiotemporal behaviors of B cells, including cyclic migration between the dark zone (DZ) and light zone (LZ) via chemotaxis; clonal proliferative bursts, somatic hypermutation, and DNA damage-induced apoptosis in the DZ; and positive selection, apoptosis via a death timer, and emergence of plasma cells in the LZ. Our simulations are able to recapitulate key molecular, cellular, and morphological GC events, including B cell population growth, affinity maturation, and clonal dominance. This novel modeling framework provides an open-source, customizable, and multiscale virtual GC simulation platform that enables qualitative and quantitative in silico investigations of a range of mechanistic and applied research questions on the adaptive humoral immune response in the future.

The Blueprint of Logical Decisions in a NF-κB Signaling System

Nearly identical cells can exhibit substantially different responses to the same stimulus that causes phenotype diversity. Such interplay between phenotype diversity and the architecture of regulatory circuits is crucial since it determines the state of a biological cell. Here, we theoretically analyze how the circuit blueprints of NF-κB in cellular environments are formed and their role in determining the cells' metabolic state. The NF-κB is a collective name for a developmental conserved family of five different transcription factors that can form homodimers or heterodimers and often promote DNA looping to reprogram the inflammatory gene response. The NF-κB controls many biological functions, including cellular differentiation, proliferation, migration, and survival. Our model shows that nuclear localization of NF-κB differentially promotes logic operations such as AND, NAND, NOR, and OR in its regulatory network. Through the quantitative thermodynamic model of transcriptional regulation and systematic variation of promoter-enhancer interaction modes, we can account for the origin of various logic gates as formed in the NF-κB system. We further show that the interconversion or switching of logic gates yielded under systematic variations of the stimuli activity and DNA looping parameters. Such computation occurs in regulatory and signaling pathways in individual cells at a molecular scale, which one can exploit to design a biomolecular computer.

SMARCA4 is a haploinsufficient B cell lymphoma tumor suppressor that fine-tunes centrocyte cell fate decisions

SMARCA4 encodes one of two mutually exclusive ATPase subunits in the BRG/BRM associated factor (BAF) complex that is recruited by transcription factors (TFs) to drive chromatin accessibility and transcriptional activation. SMARCA4 is among the most recurrently mutated genes in human cancer, including ∼30% of germinal center (GC)-derived Burkitt lymphomas. In mice, GC-specific Smarca4 haploinsufficiency cooperated with MYC over-expression to drive lymphomagenesis. Furthermore, monoallelic Smarca4 deletion drove GC hyperplasia with centroblast polarization via significantly increased rates of centrocyte recycling to the dark zone. Mechanistically, Smarca4 loss reduced the activity of TFs that are activated in centrocytes to drive GC-exit, including SPI1 (PU.1), IRF family, and NF-κB. Loss of activity for these factors phenocopied aberrant BCL6 activity within murine centrocytes and human Burkitt lymphoma cells. SMARCA4 therefore facilitates chromatin accessibility for TFs that shape centrocyte trajectories, and loss of fine-control of these programs biases toward centroblast cell-fate, GC hyperplasia and lymphoma.

A robust ultrasensitive transcriptional switch in noisy cellular environments

Ultrasensitive transcriptional switches enable sharp transitions between transcriptional on and off states and are essential for cells to respond to environmental cues with high fidelity. However, conventional switches, which rely on direct repressor-DNA binding, are extremely noise-sensitive, leading to unintended changes in gene expression. Here, through model simulations and analysis, we discovered that an alternative design combining three indirect transcriptional repression mechanisms, sequestration, blocking, and displacement, can generate a noise-resilient ultrasensitive switch. Although sequestration alone can generate an ultrasensitive switch, it remains sensitive to noise because the unintended transcriptional state induced by noise persists for long periods. However, by jointly utilizing blocking and displacement, these noise-induced transitions can be rapidly restored to the original transcriptional state. Because this transcriptional switch is effective in noisy cellular contexts, it goes beyond previous synthetic transcriptional switches, making it particularly valuable for robust synthetic system design. Our findings also provide insights into the evolution of robust ultrasensitive switches in cells. Specifically, the concurrent use of seemingly redundant indirect repression mechanisms in diverse biological systems appears to be a strategy to achieve noise-resilience of ultrasensitive switches.

Mechanosensitive super-enhancers regulate genes linked to atherosclerosis in endothelial cells

Vascular homeostasis and pathophysiology are tightly regulated by mechanical forces generated by hemodynamics. Vascular disorders such as atherosclerotic diseases largely occur at curvatures and bifurcations where disturbed blood flow activates endothelial cells while unidirectional flow at the straight part of vessels promotes endothelial health. Integrated analysis of the endothelial transcriptome, the 3D epigenome, and human genetics systematically identified the SNP-enriched cistrome in vascular endothelium subjected to well-defined atherosclerosis-prone disturbed flow or atherosclerosis-protective unidirectional flow. Our results characterized the endothelial typical- and super-enhancers and underscored the critical regulatory role of flow-sensitive endothelial super-enhancers. CRISPR interference and activation validated the function of a previously unrecognized unidirectional flow-induced super-enhancer that upregulates antioxidant genes NQO1, CYB5B, and WWP2, and a disturbed flow-induced super-enhancer in endothelium which drives prothrombotic genes EDN1 and HIVEP in vascular endothelium. Our results employing multiomics identify the cis-regulatory architecture of the flow-sensitive endothelial epigenome related to atherosclerosis and highlight the regulatory role of super-enhancers in mechanotransduction mechanisms.

A Multiscale Spatial Modeling Framework for the Germinal Center Response

The germinal center response or reaction (GCR) is a hallmark event of adaptive humoral immunity. Unfolding in the B cell follicles of the secondary lymph organs, a GC culminates in the production of high-affinity antibody-secreting plasma cells along with memory B cells. By interacting with follicular dendritic cells (FDC) and T follicular helper (Tfh) cells, GC B cells exhibit complex spatiotemporal dynamics. Driving the B cell dynamics are the intracellular signal transduction and gene regulatory network that responds to cell surface signaling molecules, cytokines, and chemokines. As our knowledge of the GC continues to expand in depth and in scope, mathematical modeling has become an important tool to help disentangle the intricacy of the GCR and inform novel mechanistic and clinical insights. While the GC has been modeled at different granularities, a multiscale spatial simulation framework - integrating molecular, cellular, and tissue-level responses - is still rare. Here, we report our recent progress toward this end with a hybrid stochastic GC framework developed on the Cellular Potts Model-based CompuCell3D platform. Tellurium is used to simulate the B cell intracellular molecular network comprising NF-κB, FOXO1, MYC, AP4, CXCR4, and BLIMP1 that responds to B cell receptor (BCR) and CD40-mediated signaling. The molecular outputs of the network drive the spatiotemporal behaviors of B cells, including cyclic migration between the dark zone (DZ) and light zone (LZ) via chemotaxis; clonal proliferative bursts, somatic hypermutation, and DNA damage-induced apoptosis in the DZ; and positive selection, apoptosis via a death timer, and emergence of plasma cells in the LZ. Our simulations are able to recapitulate key molecular, cellular, and morphological GC events including B cell population growth, affinity maturation, and clonal dominance. This novel modeling framework provides an open-source, customizable, and multiscale virtual GC simulation platform that enables qualitative and quantitative in silico investigations of a range of mechanic and applied research questions in future.

BRD9 determines the cell fate of hematopoietic stem cells by regulating chromatin state

ATP-dependent chromatin remodeling SWI/SNF complexes exist in three subcomplexes: canonical BAF (cBAF), polybromo BAF (PBAF), and a newly described non-canonical BAF (ncBAF). While cBAF and PBAF regulate fates of multiple cell types, roles for ncBAF in hematopoietic stem cells (HSCs) have not been investigated. Motivated by recent discovery of disrupted expression of BRD9, an essential component of ncBAF, in multiple cancers, including clonal hematopoietic disorders, we evaluate here the role of BRD9 in normal and malignant HSCs. BRD9 loss enhances chromatin accessibility, promoting myeloid lineage skewing while impairing B cell development. BRD9 significantly colocalizes with CTCF, whose chromatin recruitment is augmented by BRD9 loss, leading to altered chromatin state and expression of myeloid-related genes within intact topologically associating domains. These data uncover ncBAF as critical for cell fate specification in HSCs via three-dimensional regulation of gene expression and illuminate roles for ncBAF in normal and malignant hematopoiesis.

Transcriptional kinetic synergy: A complex landscape revealed by integrating modeling and synthetic biology

Transcription factors (TFs) control gene expression, often acting synergistically. Classical thermodynamic models offer a biophysical explanation for synergy based on binding cooperativity and regulated recruitment of RNA polymerase. Because transcription requires polymerase to transition through multiple states, recent work suggests that "kinetic synergy" can arise through TFs acting on distinct steps of the transcription cycle. These types of synergy are not mutually exclusive and are difficult to disentangle conceptually and experimentally. Here, we model and build a synthetic circuit in which TFs bind to a single shared site on DNA, such that TFs cannot synergize by simultaneous binding. We model mRNA production as a function of both TF binding and regulation of the transcription cycle, revealing a complex landscape dependent on TF concentration, DNA binding affinity, and regulatory activity. We use synthetic TFs to confirm that the transcription cycle must be integrated with recruitment for a quantitative understanding of gene regulation.

Quantitative Imaging Analysis of NF-κB for Mathematical Modeling Applications

Mathematical models can integrate different types of experimental datasets, reconstitute biological systems in silico, and identify previously unknown molecular mechanisms. Over the past decade, mathematical models have been developed based on quantitative observations, such as live-cell imaging and biochemical assays. However, it is difficult to directly integrate next-generation sequencing (NGS) data. Although highly dimensional, NGS data mostly only provides a "snapshot" of cellular states. Nevertheless, the development of various methods for NGS analysis has led to much more accurate predictions of transcription factor activity and has revealed various concepts regarding transcriptional regulation. Therefore, fluorescence live-cell imaging of transcription factors can help alleviate the limitations in NGS data by supplementing temporal information, linking NGS to mathematical modeling. This chapter introduces an analytical method for quantifying dynamics of nuclear factor kappaB (NF-κB) which forms aggregates in the nucleus. The method may also be applicable to other transcription factors regulated in a similar fashion.

Enhanced transcriptional heterogeneity mediated by NF-κB super-enhancers

The transcription factor NF-κB, which plays an important role in cell fate determination, is involved in the activation of super-enhancers (SEs). However, the biological functions of the NF-κB SEs in gene control are not fully elucidated. We investigated the characteristics of NF-κB-mediated SE activity using fluorescence imaging of RelA, single-cell transcriptome and chromatin accessibility analyses in anti-IgM-stimulated B cells. The formation of cell stimulation-induced nuclear RelA foci was abolished in the presence of hexanediol, suggesting an underlying process of liquid-liquid phase separation. The gained SEs induced a switch-like expression and enhanced cell-to-cell variability in transcriptional response. These properties were correlated with the number of gained cis-regulatory interactions, while switch-like gene induction was associated with the number of NF-κB binding sites in SE. Our study suggests that NF-κB SEs have an important role in the transcriptional regulation of B cells possibly through liquid condensate formation consisting of macromolecular interactions.

NF-κB produces an all-or-none activation response upon the activation of B cell receptors. These dynamics modulate the amplitude and frequency of target mRNA induction in cell populations. In this research, we performed single-cell assessment of chromatin accessibility and RNA expression, coupled with fluorescence imaging to characterize the nuclear dynamics of NF-κB proteins in B cell upon receptor stimulation. We found that upon cellular activation, NF-κB-mediated long-range activation of enhancers cooperatively evoked RNA production. In addition, predicted DNA contacts brought by open chromatin led to the high heterogeneity of RNA levels in cell populations. Stimuli-dependent NF-κB foci formation was further inhibited by 1,6-hexanediol (liquid-liquid phase separation inhibitor) and JQ1 (coactivator protein BRD4 inhibitor). We thus propose that nuclear NF-κB plays an important role in the transcriptional regulation of B cell development possibly through the formation of liquid condensates.

Systems approaches to investigate the role of NF-κB signaling in aging

The nuclear factor-κB (NF-κB) signaling pathway is one of the most well-studied pathways related to inflammation, and its involvement in aging has attracted considerable attention. As aging is a complex phenomenon and is the result of a multi-step process, the involvement of the NF-κB pathway in aging remains unclear. To elucidate the role of NF-κB in the regulation of aging, different systems biology approaches have been employed. A multi-omics data-driven approach can be used to interpret and clarify unknown mechanisms but cannot generate mechanistic regulatory structures alone. In contrast, combining this approach with a mathematical modeling approach can identify the mechanistics of the phenomena of interest. The development of single-cell technologies has also helped clarify the heterogeneity of the NF-κB response and underlying mechanisms. Here, we review advances in the understanding of the regulation of aging by NF-κB by focusing on omics approaches, single-cell analysis, and mathematical modeling of the NF-κB network.

OCT4-mediated inflammation induces cell reprogramming at the origin of cardiac valve development and calcification

Cell plasticity plays a key role in embryos by maintaining the differentiation potential of progenitors. Whether postnatal somatic cells revert to an embryonic-like naïve state regaining plasticity and redifferentiate into a cell type leading to a disease remains intriguing. Using genetic lineage tracing and single-cell RNA sequencing, we reveal that Oct4 is induced by nuclear factor κB (NFκB) at embyronic day 9.5 in a subset of mouse endocardial cells originating from the anterior heart forming field at the onset of endocardial-to-mesenchymal transition. These cells acquired a chondro-osteogenic fate. OCT4 in adult valvular aortic cells leads to calcification of mouse and human valves. These calcifying cells originate from the Oct4 embryonic lineage. Genetic deletion of Pou5f1 (Pit-Oct-Unc, OCT4) in the endocardial cell lineage prevents aortic stenosis and calcification of ApoE^−/− mouse valve. We established previously unidentified self-cell reprogramming NFκB- and OCT4-mediated inflammatory pathway triggering a dose-dependent mechanism of valve calcification.

Super Enhancer Profiles Identify Key Cell Identity Genes During Differentiation From Embryonic Stem Cells to Trophoblast Stem Cells Super Enhencers in Trophoblast Differentiation

Trophoblast stem cells (TSCs) are derived from blastocysts and the extra-embryonic ectoderm (ExE) of post-implantation embryos and play a significant role in fetal development, but the roles that TSCs play in the earlier status of fetal diseases need further exploration. Super enhancers (SEs) are dense clusters of stitched enhancers that control cell identity determination and disease development and may participate in TSC differentiation. We identified key cell identity genes regulated by TSC-SEs via integrated analysis of H3K27ac and H3K4me1 chromatin immunoprecipitation sequencing (ChIP-seq), RNA-sequencing (RNA-seq) and ATAC-sequencing (ATAC-seq) data. The identified key TSC identity genes regulated by SEs, such as epidermal growth factor receptor (EGFR), integrin β5 (ITGB5) and Paxillin (Pxn), were significantly upregulated during TSC differentiation, and the transcription network mediated by TSC-SEs enriched in terms like focal adhesion and actin cytoskeleton regulation related to differentiation of TSCs. Additionally, the increased chromatin accessibility of the key cell identity genes verified by ATAC-seq further demonstrated the regulatory effect of TSC-SEs on TSC lineage commitment. Our results illustrated the significant roles of the TSC-SE-regulated network in TSC differentiation, and identified key TSC identity genes EGFR, ITGB5 and Pxn, providing novel insight into TSC differentiation and lays the foundation for future studies on embryo implantation and related diseases.

The transcription factor reservoir and chromatin landscape in activated plasmacytoid dendritic cells

Background

Transcription factors (TFs) control gene expression by direct binding to regulatory regions of target genes but also by impacting chromatin landscapes and modulating DNA accessibility for other TFs. In recent years several TFs have been defined that control cell fate decisions and effector functions in the immune system. Plasmacytoid dendritic cells (pDCs) are an immune cell type with the unique capacity to produce high amounts of type I interferons quickly in response to contact with viral components. Hereby, this cell type is involved in anti-infectious immune responses but also in the development of inflammatory and autoimmune diseases. To date, the global TF reservoir in pDCs early after activation remains to be fully characterized.

Results

To fill this gap, we have performed a comprehensive analysis in naïve versus TLR9-activated murine pDCs in a time course study covering early timepoints after stimulation (2 h, 6 h, 12 h) integrating gene expression (RNA-Seq) and chromatin landscape (ATAC-Seq) studies. To unravel the biological processes underlying the changes in TF expression on a global scale gene ontology (GO) analyses were performed. We found that 70% of all genes annotated as TFs in the mouse genome (1014 out of 1636) are expressed in pDCs for at least one stimulation time point and are covering a wide range of TF classes defined by their specific DNA binding mechanisms. GO analysis revealed involvement of TLR9-induced TFs in epigenetic modulation, NFκB and JAK-STAT signaling, and protein production in the endoplasmic reticulum. pDC activation predominantly “turned on” the chromatin regions associated with TF genes. Our in silico analyses pointed at the AP-1 family of TFs as less noticed but possibly important players in these cells after activation. AP-1 family members exhibit (1) increased gene expression, (2) enhanced chromatin accessibility in their promoter region, and (3) a TF DNA binding motif that is globally enriched in genomic regions that were found more accessible in pDCs after TLR9 activation.

Conclusions

In this study we define the complete set of TLR9-regulated TFs in pDCs. Further, this study identifies the AP-1 family of TFs as potentially important but so far less well characterized regulators of pDC function.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-021-00991-2.

Potential roles of super enhancers in inflammatory gene transcription

Acute and chronic inflammation is a basic pathological event that contributes to atherosclerosis, cancer, infectious diseases, and immune disorders. Inflammation is an adaptive process to both external and internal stimuli experienced by the human body. Although the mechanism of gene transcription is highly complicated and orchestrated in a timely and spatial manner, recent developments in next-generation sequencing, genome-editing, cryo-electron microscopy, and single cell-based technologies could provide us with insights into the roles of super enhancers (SEs). Initially, SEs were implicated in determining cell fate; subsequent studies have clarified that SEs are associated with various pathological conditions, including cancer and inflammatory diseases. Recent technological advances have unveiled the molecular mechanisms of SEs, which involve epigenetic histone modifications, chromatin three-dimensional structures, and phase-separated condensates. In this review, we discuss the relationship between inflammation and SEs and the therapeutic potential of SEs for inflammatory diseases.

Regulation of Transcription Factor NF-κB in Its Natural Habitat: The Nucleus

Activation of the transcription factor NF-κB elicits an individually tailored transcriptional response in order to meet the particular requirements of specific cell types, tissues, or organs. Control of the induction kinetics, amplitude, and termination of gene expression involves multiple layers of NF-κB regulation in the nucleus. Here we discuss some recent advances in our understanding of the mutual relations between NF-κB and chromatin regulators also in the context of different levels of genome organization. Changes in the 3D folding of the genome, as they occur during senescence or in cancer cells, can causally contribute to sustained increases in NF-κB activity. We also highlight the participation of NF-κB in the formation of hierarchically organized super enhancers, which enable the coordinated expression of co-regulated sets of NF-κB target genes. The identification of mechanisms allowing the specific regulation of NF-κB target gene clusters could potentially enable targeted therapeutic interventions, allowing selective interference with subsets of the NF-κB response without a complete inactivation of this key signaling system.

Cell shape-based chemical screening reveals an epigenetic network mediated by focal adhesions

Adapter proteins CRK and CRKL participate in a variety of signaling pathways, including cell adhesion, and fate regulation of mammalian cells. However, the molecular functions of CRK/CRKL in epigenetic regulation remain largely unknown. Here, we developed a pipeline to evaluate cell morphology using high-content image analysis combined with chemical screening of kinase and epigenetic modulators. We found that CRK/CRKL modulates gene regulatory networks associated with cell morphology through epigenetic alteration in mouse embryonic fibroblasts. Integrated epigenome and transcriptome analyses revealed that CRK/CRKL is involved in super-enhancer activity and upregulation of Cdt1, Rin1, and Spp1 expression for the regulation of cell morphology. Screening of a library of 80 epigenetic inhibitors showed that histone H3 modifiers, euchromatic histone methyltransferase 2 and mitogen- and stress-activated kinase 1, may be important for CRK/CRKL-mediated morphological changes. Taken together, our results indicate that CRK/CRKL plays a critical role in gene regulatory networks through epigenetic modification. DATABASES: Chromatin immunoprecipitation sequencing and RNA sequencing data were deposited in the DNA Data Bank of Japan under DRA011080 and DRA011081 accession numbers, respectively.

Building patient-specific models for receptor tyrosine kinase signaling networks

Cancer progresses due to changes in the dynamic interactions of multidimensional factors associated with gene mutations. Cancer research has actively adopted computational methods, including data-driven and mathematical model-driven approaches, to identify causative factors and regulatory rules that can explain the complexity and diversity of cancers. A data-driven, statistics-based approach revealed correlations between gene alterations and clinical outcomes in many types of cancers. A model-driven mathematical approach has elucidated the dynamic features of cancer networks and identified the mechanisms of drug efficacy and resistance. More recently, machine learning methods have emerged that can be used for mining omics data and classifying patient. However, as the strengths and weaknesses of each method becoming apparent, new analytical tools are emerging to combine and improve the methodologies and maximize their predictive power for classifying cancer subtypes and prognosis. Here, we introduce recent advances in cancer systems biology aimed at personalized medicine, with focus on the receptor tyrosine kinase signaling network.