Repair of ischemic tissue injury by Wnt-mediated self-renewal and myofibroblasts expansion repres... more Repair of ischemic tissue injury by Wnt-mediated self-renewal and myofibroblasts expansion represents a trade-off because myofibroblasts can promote chronic fibrosis. We used a cell interaction network model to analyze the double-edged effect of Wnt-signaling in mouse ischemic kidney disease. The model predicts multistable dynamics that mapped to three distinct outcomes observed in surviving mice as a function of insult severity: full recovery, recovery with residual and chronic fibrosis. Furthermore we predicted and experimentally verified protective effect of preconditioning by mild ischemia at the cost of fibrosis risk. The model allowed for multi-objective optimization to overcome the trade-off between rapid repair and chronification and guided a successful biphasic treatment with Wnt-agonist followed by Wnt-antagonist. The results underscore the importance of non-linear disease dynamics and propose a non-monotonic targeted therapy by phasic blockade and activation of the same t...
Biological systems need robust mechanism to recover from injury. Increasing evidences suggest tha... more Biological systems need robust mechanism to recover from injury. Increasing evidences suggest that incomplete recovery of acute kidney injury (AKI) increases the risk of chronic kidney diseases (CKD), but it remains ambiguous how AKI and CKD, traditionally treated as two distinct pathological entities, are mechanistically related. With different durations of ischemia-reperfusion injury in mice, we identified a unified repair process that leads to perfect or imperfect adaptive response, or maladaptive response and even death depending on injury severity. Our mathematical modeling revealed that these multi-state outcomes result from a trade-off between two functional requirements imposed on the repair system, fast recovery to reduce the risk of immediate organ failure and organism death, and complete resolution to reduce the risk of CKD. Subsequent experiments confirmed this trade-off mechanism, and demonstrated that targeting the temporal dynamics, instead of traditional monotonic me...
Phenotype transition takes place in many biological processes such as differentiation, and unders... more Phenotype transition takes place in many biological processes such as differentiation, and understanding how a cell reprograms its global gene expression profile is a problem of rate theories. A cell phenotype transition accompanies with switching of expression rates of clusters of genes, analogous to domain flipping in an Ising system. Here through analyzing single cell RNA sequencing data in the framework of transition path theory, we set to study how such a genome-wide expression program switching proceeds in three different cell transition processes. For each process after reconstructing a Markov transition model in the cell state space, we formed an ensemble of shortest paths connecting the initial and final cell states, reconstructed a reaction coordinate describing the transition progression, and inferred the gene regulation network (GRN) along the reaction coordinate. In all three processes we observed common pattern that the frustration of gene regulatory network (GRN), def...
How do mammalian cells that share the same genome exist in notably distinct phenotypes, exhibitin... more How do mammalian cells that share the same genome exist in notably distinct phenotypes, exhibiting differences in morphology, gene expression patterns, and epigenetic chromatin statuses? Furthermore how do cells of different phenotypes differentiate reproducibly from a single fertilized egg? These fundamental questions are closely related to a deeply rooted paradigm in developmental biology that cell differentiation is irreversible. Yet, recently a growing body of research suggests the possibility of cell reprogramming, which offers the potential for us to convert one type of cell into another. Despite the significance of quantitative understandings of cell reprogramming, theoretical efforts often suffer from the complexity of large circuits maintaining cell phenotypes coupled at many different epigenetic and gene regulation levels. To capture the global architecture of cell phenotypes, we propose an "epigenetic state network" approach that translates the classical concept...
Live-cell imaging in multidimensional feature space reveals parallel transition paths of epitheli... more Live-cell imaging in multidimensional feature space reveals parallel transition paths of epithelial-to-mesenchymal transition.
Tissue fibrosis compromises organ function and occurs as a potential long-term outcome in respons... more Tissue fibrosis compromises organ function and occurs as a potential long-term outcome in response to acute tissue injuries. Currently, lack of mechanistic understanding prevents effective prevention and treatment of the progression from acute injury to fibrosis. Here, we combined quantitative experimental studies with a mouse kidney injury model and a computational approach to determine how the physiological consequences are determined by the severity of ischemia injury, and to identify how to manipulate Wnt signaling to accelerate repair of ischemic tissue damage while minimizing fibrosis. The study reveals that Wnt-mediated memory of prior injury contributes to fibrosis progression, and ischemic preconditioning reduces the risk of death but increases the risk of fibrosis. Furthermore, we validated the prediction that sequential combination therapy of initial treatment with a Wnt agonist followed by treatment with a Wnt antagonist can reduce both the risk of death and fibrosis in response to acute injuries.
ABSTRACTRecent advances in single cell techniques catalyze an emerging field of studying how cell... more ABSTRACTRecent advances in single cell techniques catalyze an emerging field of studying how cells convert from one phenotype to another, in a step-by-step process. Two grand technical challenges, however, impede further development of the field. Fixed cell-based approaches can provide genome-wide snapshots of cell status but have fundamental limits on revealing temporal information, and fluorescence-based live cell imaging approaches provide temporal information but are technically challenging for multiplex long- term imaging. We present a live-cell imaging platform, Multiplex Trajectory Recording and Analysis of Cellular Kinetics, or M-TRACK, that tracks cellular status change through combining endogenous fluorescent labeling that minimizes perturbation to cell physiology and live cell imaging of high-dimensional cell morphological and texture features. To test the functionality of our platform, we used the A549 VIM-RFP EMT reporter line (ATCC® CCL-185EMT) to explore cell phenotyp...
ABSTRACTCRISPR-based gene knock-in at endogenous sites is desirable in multiple fields such as qu... more ABSTRACTCRISPR-based gene knock-in at endogenous sites is desirable in multiple fields such as quantitative studies of signal transduction pathways and gene regulation, synthetic biology, and disease modeling. Contrasting the knock-out procedure, a key step of CRISPR knock-in procedure relies on the homology-directed repairing (HDR) process that requires a donor construct as repair template. Therefore, it is desirable to generate a series of donor DNA constructs efficiently and cost-effectively. In this study, we developed a general Gibson assembly procedure that combines strengths of a Modular Overlap-Directed Assembly with Linkers (MODAL) strategy and a restriction enzyme based hierarchical framework. This procedure also allows fusing sgRNAs to the constructs for enhanced homology-directed repairing efficiency. Experimental tests on multiple constructs achieved from 3-8 folds of increase in assembly efficiency to high yield of constructs that failed to make with conventional Gibso...
Many cellular responses to surrounding cues require temporally concerted transcriptional regulati... more Many cellular responses to surrounding cues require temporally concerted transcriptional regulation of multiple genes. In prokaryotic cells, a single-input-module motif with one transcription factor regulating multiple target genes can generate coordinated gene expression. In eukaryotic cells, transcriptional activity of a gene is affected by not only transcription factors but also the epigenetic modifications and three-dimensional Author Summary Cellular responses to environmental stimulation are often accompanied by changes in gene expression patterns. Genes are linearly arranged along chromosomal DNA, which folds into a three-dimensional structure. The chromosome structure affects gene expression activities and is regulated by multiple events such as histone modifications and DNA binding of transcription factors. A basic question is how these mechanisms work together to regulate gene expression. In this study, we analyzed temporal gene expression patterns in the context of chromosome structure both in a human cell line under TGF-β treatment and during mouse nervous system development. In both cases, we observed that genes regulated by common transcription factors have an enhanced tendency to be spatially close. Our analysis suggests that spatial co-localization of genes may facilitate the concerted gene expression.
The transition between epithelial and mesenchymal (EMT) is a fundamental cellular process that pl... more The transition between epithelial and mesenchymal (EMT) is a fundamental cellular process that plays critical roles in development, cancer metastasis, and tissue wound healing. EMT is not a binary process but involves multiple partial EMT states that give rise to a high degree of cell state plasticity. Here, we first reviewed several studies on theoretical predictions and experimental verification of these intermediate states, the role of partial EMT on kidney fibrosis development, and how quantitative signaling information controls cell commitment to partial or full EMT upon transient signals. Next, we summarized existing knowledge and open questions on the coupling between EMT and other biological processes, such as the cell cycle, epigenetic regulation, stemness, and apoptosis. Taken together, EMT is a model system that has attracted increasing interests for quantitative experimental and theoretical studies.
Journal of the American Society of Nephrology : JASN, Jan 17, 2018
AKI is a devastating condition with high morbidity and mortality. The pathologic features of AKI ... more AKI is a devastating condition with high morbidity and mortality. The pathologic features of AKI are characterized by tubular injury, inflammation, and vascular impairment. Whether fibroblasts in the renal interstitium have a role in the pathogenesis of AKI is unknown. In this study, we investigated the role of fibroblast-specific β-catenin signaling in dictating the outcome of AKI, using conditional knockout mice in which β-catenin was specifically ablated in fibroblasts (Gli1-β-cat-/-). After ischemia-reperfusion injury (IRI), Gli1-β-cat-/- mice had lower serum creatinine levels and less morphologic injury than Gli1-β-cat+/+ littermate controls. Moreover, we detected fewer apoptotic cells, as well as decreased cytochrome C release; reduced expression of Bax, FasL, and p53; and increased phosphorylation of Akt, in the Gli1-β-cat-/- kidneys. Gli1-β-cat-/- kidneys also exhibited upregulated expression of proliferating cell nuclear antigen and Ki-67, which are markers of cell prolifer...
Targeting microenvironmental factors that foster migratory cell phenotypes is a promising strateg... more Targeting microenvironmental factors that foster migratory cell phenotypes is a promising strategy for halting tumor migration. However, lack of mechanistic understanding of the process impedes pharmaceutical drug development. Using a novel 3D microtumor model with tight control over tumor size, we recapitulated tumor size-induced hypoxic microenvironment and emergence of migratory phenotypes in epithelial T47D breast microtumors as well as those of patient-derived primary metastatic breast cancer cells, mesothelioma cells and lung cancer xenograft cells (PDX). The microtumor models from various patient-derived tumor cells and PDX cells revealed upregulation of tumor secretome, matrix metalloproteinase-9 (MMP9), fibronectin (FN), and soluble E-cadherin (sE-CAD) consistent with the clinically reported elevated levels of FN and MMP9 in the patient breast tumors compared to healthy mammary gland. We further showed that the tumor secretome induces migratory phenotype in non-hypoxic, non...
The detection and transmission of the temporal quality of intracellular and extracellular signals... more The detection and transmission of the temporal quality of intracellular and extracellular signals is an essential cellular mechanism. It remains largely unexplored how cells interpret the duration information of a stimulus. In this paper, through an integrated quantitative and computational approach we demonstrate that crosstalk among multiple TGF-β activated pathways forms a relay from SMAD to GLI1 that initializes and maintains SNAILl expression, respectively. This transaction is smoothed and accelerated by another temporal switch from elevated cytosolic GSK3 enzymatic activity to reduced nuclear GSK3 enzymatic activity. The intertwined network places SNAIL1 as a key integrator of information from TGF-β signaling subsequently distributed through upstream divergent pathways; essentially cells generate a transient or sustained expression of SNAIL1 depending on TGF-β duration. Other signaling pathways may use similar network structure to encode duration information.
miRNAs serve as crucial post-transcriptional regulators in a variety of essential cell fate decis... more miRNAs serve as crucial post-transcriptional regulators in a variety of essential cell fate decisions. However, the contribution of the mRNA-miRNA mutual regulation to bistability is not fully understood. Here, we built a set of mathematical models of mRNA-miRNA interactions and systematically analyzed the sensitivity of the response curves under various conditions. Our findings indicate that mRNA-miRNA reciprocal regulation could manifest ultrasensitivity to subserve the generation of bistability when equipped with a positive feedback loop. We also find that the region of bistability is expanded by a stronger competing endogenous mRNA (ceRNA). Interestingly, bistability can be generated without a feedback loop if multiple miRNA binding sites exist on a target mRNA. Thus, we demonstrate the importance of simple mRNA-miRNA reciprocal regulation in cell fate decisions. This article is protected by copyright. All rights reserved.
Multiple-objective optimization is common in biological systems. In the mammalian olfactory syste... more Multiple-objective optimization is common in biological systems. In the mammalian olfactory system, each sensory neuron stochastically expresses only one out of up to thousands of olfactory receptor (OR) gene alleles; at organism level the types of expressed ORs need to be maximized. Existing models focus only on monoallele activation, and cannot explain recent observations in mutants, especially the reduced global diversity of expressed ORs in G9a/GLP knockouts. In this work we integrated existing information on OR expression, and constructed a comprehensive model that has all its components based on physical interactions. Analyzing the model reveals an evolutionarily optimized three-layer regulation mechanism, which includes zonal segregation, epigenetic barrier crossing coupled to a negative feedback loop that mechanistically differs from previous theoretical proposals, and a previously unidentified enhancer competition step. This model not only recapitulates monoallelic OR expression, but also elucidates how the olfactory system maximizes and maintains the diversity of OR expression, and has multiple predictions validated by existing experimental results. Through making analogy to a physical system with thermally activated barrier crossing and comparative reverse engineering analyses, the study reveals that the olfactory receptor selection system is optimally designed, and particularly underscores cooperativity and synergy as a general design principle for multi-objective optimization in biology.
Proceedings of the National Academy of Sciences, 2016
Multiple-objective optimization is common in biological systems. In the mammalian olfactory syste... more Multiple-objective optimization is common in biological systems. In the mammalian olfactory system, each sensory neuron stochastically expresses only one out of up to thousands of olfactory receptor (OR) gene alleles; at the organism level, the types of expressed ORs need to be maximized. Existing models focus only on monoallele activation, and cannot explain recent observations in mutants, especially the reduced global diversity of expressed ORs in G9a/GLP knockouts. In this work we integrated existing information on OR expression, and constructed a comprehensive model that has all its components based on physical interactions. Analyzing the model reveals an evolutionarily optimized three-layer regulation mechanism, which includes zonal segregation, epigenetic barrier crossing coupled to a negative feedback loop that mechanistically differs from previous theoretical proposals, and a previously unidentified enhancer competition step. This model not only recapitulates monoallelic OR ...
How do mammalian cells that share the same genome exist in notably distinct phenotypes, exhibitin... more How do mammalian cells that share the same genome exist in notably distinct phenotypes, exhibiting differences in morphology, gene expression patterns, and epigenetic chromatin statuses? Furthermore how do cells of different phenotypes differentiate reproducibly from a single fertilized egg? These are fundamental problems in developmental biology. Epigenetic histone modifications play an important role in the maintenance of different cell phenotypes. The exact molecular mechanism for inheritance of the modification patterns over cell generations remains elusive. The complexity comes partly from the number of molecular species and the broad time scales involved. In recent years mathematical modeling has made significant contributions on elucidating the molecular mechanisms of DNA methylation and histone Book chapter in Epigenetic Technological Applications (Elsevier), in press 2 covalent modification inheritance. We will pedagogically introduce the typical procedure and some technical details of performing a mathematical modeling study, and discuss future developments.
Experimental and computational analyses identify multiple cell populations during the epithelial-... more Experimental and computational analyses identify multiple cell populations during the epithelial-to-mesenchymal transition.
A focused theme in systems biology is to uncover design principles of biological networks, that i... more A focused theme in systems biology is to uncover design principles of biological networks, that is, how specific network structures yield specific systems properties. For this purpose, we have previously developed a reverse engineering procedure to identify network topologies with high likelihood in generating desired systems properties. Our method searches the continuous parameter space of an assembly of network topologies, without enumerating individual network topologies separately as traditionally done in other reverse engineering procedures. Here we tested this CPSS (continuous parameter space search) method on a previously studied problem: the resettable bistability of an Rb-E2F gene network in regulating the quiescence-to-proliferation transition of mammalian cells. From a simplified Rb-E2F gene network, we identified network topologies responsible for generating resettable bistability. The CPSS-identified topologies are consistent with those reported in the previous study based on individual topology search (ITS), demonstrating the effectiveness of the CPSS approach. Since the CPSS and ITS searches are based on different mathematical formulations and different algorithms, the consistency of the results also helps cross-validate both approaches. A unique advantage of the CPSS approach lies in its applicability to biological networks with large numbers of nodes. To aid the application of the CPSS approach to the study of other biological systems, we have developed a computer package that is available in Information S1.
Repair of ischemic tissue injury by Wnt-mediated self-renewal and myofibroblasts expansion repres... more Repair of ischemic tissue injury by Wnt-mediated self-renewal and myofibroblasts expansion represents a trade-off because myofibroblasts can promote chronic fibrosis. We used a cell interaction network model to analyze the double-edged effect of Wnt-signaling in mouse ischemic kidney disease. The model predicts multistable dynamics that mapped to three distinct outcomes observed in surviving mice as a function of insult severity: full recovery, recovery with residual and chronic fibrosis. Furthermore we predicted and experimentally verified protective effect of preconditioning by mild ischemia at the cost of fibrosis risk. The model allowed for multi-objective optimization to overcome the trade-off between rapid repair and chronification and guided a successful biphasic treatment with Wnt-agonist followed by Wnt-antagonist. The results underscore the importance of non-linear disease dynamics and propose a non-monotonic targeted therapy by phasic blockade and activation of the same t...
Biological systems need robust mechanism to recover from injury. Increasing evidences suggest tha... more Biological systems need robust mechanism to recover from injury. Increasing evidences suggest that incomplete recovery of acute kidney injury (AKI) increases the risk of chronic kidney diseases (CKD), but it remains ambiguous how AKI and CKD, traditionally treated as two distinct pathological entities, are mechanistically related. With different durations of ischemia-reperfusion injury in mice, we identified a unified repair process that leads to perfect or imperfect adaptive response, or maladaptive response and even death depending on injury severity. Our mathematical modeling revealed that these multi-state outcomes result from a trade-off between two functional requirements imposed on the repair system, fast recovery to reduce the risk of immediate organ failure and organism death, and complete resolution to reduce the risk of CKD. Subsequent experiments confirmed this trade-off mechanism, and demonstrated that targeting the temporal dynamics, instead of traditional monotonic me...
Phenotype transition takes place in many biological processes such as differentiation, and unders... more Phenotype transition takes place in many biological processes such as differentiation, and understanding how a cell reprograms its global gene expression profile is a problem of rate theories. A cell phenotype transition accompanies with switching of expression rates of clusters of genes, analogous to domain flipping in an Ising system. Here through analyzing single cell RNA sequencing data in the framework of transition path theory, we set to study how such a genome-wide expression program switching proceeds in three different cell transition processes. For each process after reconstructing a Markov transition model in the cell state space, we formed an ensemble of shortest paths connecting the initial and final cell states, reconstructed a reaction coordinate describing the transition progression, and inferred the gene regulation network (GRN) along the reaction coordinate. In all three processes we observed common pattern that the frustration of gene regulatory network (GRN), def...
How do mammalian cells that share the same genome exist in notably distinct phenotypes, exhibitin... more How do mammalian cells that share the same genome exist in notably distinct phenotypes, exhibiting differences in morphology, gene expression patterns, and epigenetic chromatin statuses? Furthermore how do cells of different phenotypes differentiate reproducibly from a single fertilized egg? These fundamental questions are closely related to a deeply rooted paradigm in developmental biology that cell differentiation is irreversible. Yet, recently a growing body of research suggests the possibility of cell reprogramming, which offers the potential for us to convert one type of cell into another. Despite the significance of quantitative understandings of cell reprogramming, theoretical efforts often suffer from the complexity of large circuits maintaining cell phenotypes coupled at many different epigenetic and gene regulation levels. To capture the global architecture of cell phenotypes, we propose an "epigenetic state network" approach that translates the classical concept...
Live-cell imaging in multidimensional feature space reveals parallel transition paths of epitheli... more Live-cell imaging in multidimensional feature space reveals parallel transition paths of epithelial-to-mesenchymal transition.
Tissue fibrosis compromises organ function and occurs as a potential long-term outcome in respons... more Tissue fibrosis compromises organ function and occurs as a potential long-term outcome in response to acute tissue injuries. Currently, lack of mechanistic understanding prevents effective prevention and treatment of the progression from acute injury to fibrosis. Here, we combined quantitative experimental studies with a mouse kidney injury model and a computational approach to determine how the physiological consequences are determined by the severity of ischemia injury, and to identify how to manipulate Wnt signaling to accelerate repair of ischemic tissue damage while minimizing fibrosis. The study reveals that Wnt-mediated memory of prior injury contributes to fibrosis progression, and ischemic preconditioning reduces the risk of death but increases the risk of fibrosis. Furthermore, we validated the prediction that sequential combination therapy of initial treatment with a Wnt agonist followed by treatment with a Wnt antagonist can reduce both the risk of death and fibrosis in response to acute injuries.
ABSTRACTRecent advances in single cell techniques catalyze an emerging field of studying how cell... more ABSTRACTRecent advances in single cell techniques catalyze an emerging field of studying how cells convert from one phenotype to another, in a step-by-step process. Two grand technical challenges, however, impede further development of the field. Fixed cell-based approaches can provide genome-wide snapshots of cell status but have fundamental limits on revealing temporal information, and fluorescence-based live cell imaging approaches provide temporal information but are technically challenging for multiplex long- term imaging. We present a live-cell imaging platform, Multiplex Trajectory Recording and Analysis of Cellular Kinetics, or M-TRACK, that tracks cellular status change through combining endogenous fluorescent labeling that minimizes perturbation to cell physiology and live cell imaging of high-dimensional cell morphological and texture features. To test the functionality of our platform, we used the A549 VIM-RFP EMT reporter line (ATCC® CCL-185EMT) to explore cell phenotyp...
ABSTRACTCRISPR-based gene knock-in at endogenous sites is desirable in multiple fields such as qu... more ABSTRACTCRISPR-based gene knock-in at endogenous sites is desirable in multiple fields such as quantitative studies of signal transduction pathways and gene regulation, synthetic biology, and disease modeling. Contrasting the knock-out procedure, a key step of CRISPR knock-in procedure relies on the homology-directed repairing (HDR) process that requires a donor construct as repair template. Therefore, it is desirable to generate a series of donor DNA constructs efficiently and cost-effectively. In this study, we developed a general Gibson assembly procedure that combines strengths of a Modular Overlap-Directed Assembly with Linkers (MODAL) strategy and a restriction enzyme based hierarchical framework. This procedure also allows fusing sgRNAs to the constructs for enhanced homology-directed repairing efficiency. Experimental tests on multiple constructs achieved from 3-8 folds of increase in assembly efficiency to high yield of constructs that failed to make with conventional Gibso...
Many cellular responses to surrounding cues require temporally concerted transcriptional regulati... more Many cellular responses to surrounding cues require temporally concerted transcriptional regulation of multiple genes. In prokaryotic cells, a single-input-module motif with one transcription factor regulating multiple target genes can generate coordinated gene expression. In eukaryotic cells, transcriptional activity of a gene is affected by not only transcription factors but also the epigenetic modifications and three-dimensional Author Summary Cellular responses to environmental stimulation are often accompanied by changes in gene expression patterns. Genes are linearly arranged along chromosomal DNA, which folds into a three-dimensional structure. The chromosome structure affects gene expression activities and is regulated by multiple events such as histone modifications and DNA binding of transcription factors. A basic question is how these mechanisms work together to regulate gene expression. In this study, we analyzed temporal gene expression patterns in the context of chromosome structure both in a human cell line under TGF-β treatment and during mouse nervous system development. In both cases, we observed that genes regulated by common transcription factors have an enhanced tendency to be spatially close. Our analysis suggests that spatial co-localization of genes may facilitate the concerted gene expression.
The transition between epithelial and mesenchymal (EMT) is a fundamental cellular process that pl... more The transition between epithelial and mesenchymal (EMT) is a fundamental cellular process that plays critical roles in development, cancer metastasis, and tissue wound healing. EMT is not a binary process but involves multiple partial EMT states that give rise to a high degree of cell state plasticity. Here, we first reviewed several studies on theoretical predictions and experimental verification of these intermediate states, the role of partial EMT on kidney fibrosis development, and how quantitative signaling information controls cell commitment to partial or full EMT upon transient signals. Next, we summarized existing knowledge and open questions on the coupling between EMT and other biological processes, such as the cell cycle, epigenetic regulation, stemness, and apoptosis. Taken together, EMT is a model system that has attracted increasing interests for quantitative experimental and theoretical studies.
Journal of the American Society of Nephrology : JASN, Jan 17, 2018
AKI is a devastating condition with high morbidity and mortality. The pathologic features of AKI ... more AKI is a devastating condition with high morbidity and mortality. The pathologic features of AKI are characterized by tubular injury, inflammation, and vascular impairment. Whether fibroblasts in the renal interstitium have a role in the pathogenesis of AKI is unknown. In this study, we investigated the role of fibroblast-specific β-catenin signaling in dictating the outcome of AKI, using conditional knockout mice in which β-catenin was specifically ablated in fibroblasts (Gli1-β-cat-/-). After ischemia-reperfusion injury (IRI), Gli1-β-cat-/- mice had lower serum creatinine levels and less morphologic injury than Gli1-β-cat+/+ littermate controls. Moreover, we detected fewer apoptotic cells, as well as decreased cytochrome C release; reduced expression of Bax, FasL, and p53; and increased phosphorylation of Akt, in the Gli1-β-cat-/- kidneys. Gli1-β-cat-/- kidneys also exhibited upregulated expression of proliferating cell nuclear antigen and Ki-67, which are markers of cell prolifer...
Targeting microenvironmental factors that foster migratory cell phenotypes is a promising strateg... more Targeting microenvironmental factors that foster migratory cell phenotypes is a promising strategy for halting tumor migration. However, lack of mechanistic understanding of the process impedes pharmaceutical drug development. Using a novel 3D microtumor model with tight control over tumor size, we recapitulated tumor size-induced hypoxic microenvironment and emergence of migratory phenotypes in epithelial T47D breast microtumors as well as those of patient-derived primary metastatic breast cancer cells, mesothelioma cells and lung cancer xenograft cells (PDX). The microtumor models from various patient-derived tumor cells and PDX cells revealed upregulation of tumor secretome, matrix metalloproteinase-9 (MMP9), fibronectin (FN), and soluble E-cadherin (sE-CAD) consistent with the clinically reported elevated levels of FN and MMP9 in the patient breast tumors compared to healthy mammary gland. We further showed that the tumor secretome induces migratory phenotype in non-hypoxic, non...
The detection and transmission of the temporal quality of intracellular and extracellular signals... more The detection and transmission of the temporal quality of intracellular and extracellular signals is an essential cellular mechanism. It remains largely unexplored how cells interpret the duration information of a stimulus. In this paper, through an integrated quantitative and computational approach we demonstrate that crosstalk among multiple TGF-β activated pathways forms a relay from SMAD to GLI1 that initializes and maintains SNAILl expression, respectively. This transaction is smoothed and accelerated by another temporal switch from elevated cytosolic GSK3 enzymatic activity to reduced nuclear GSK3 enzymatic activity. The intertwined network places SNAIL1 as a key integrator of information from TGF-β signaling subsequently distributed through upstream divergent pathways; essentially cells generate a transient or sustained expression of SNAIL1 depending on TGF-β duration. Other signaling pathways may use similar network structure to encode duration information.
miRNAs serve as crucial post-transcriptional regulators in a variety of essential cell fate decis... more miRNAs serve as crucial post-transcriptional regulators in a variety of essential cell fate decisions. However, the contribution of the mRNA-miRNA mutual regulation to bistability is not fully understood. Here, we built a set of mathematical models of mRNA-miRNA interactions and systematically analyzed the sensitivity of the response curves under various conditions. Our findings indicate that mRNA-miRNA reciprocal regulation could manifest ultrasensitivity to subserve the generation of bistability when equipped with a positive feedback loop. We also find that the region of bistability is expanded by a stronger competing endogenous mRNA (ceRNA). Interestingly, bistability can be generated without a feedback loop if multiple miRNA binding sites exist on a target mRNA. Thus, we demonstrate the importance of simple mRNA-miRNA reciprocal regulation in cell fate decisions. This article is protected by copyright. All rights reserved.
Multiple-objective optimization is common in biological systems. In the mammalian olfactory syste... more Multiple-objective optimization is common in biological systems. In the mammalian olfactory system, each sensory neuron stochastically expresses only one out of up to thousands of olfactory receptor (OR) gene alleles; at organism level the types of expressed ORs need to be maximized. Existing models focus only on monoallele activation, and cannot explain recent observations in mutants, especially the reduced global diversity of expressed ORs in G9a/GLP knockouts. In this work we integrated existing information on OR expression, and constructed a comprehensive model that has all its components based on physical interactions. Analyzing the model reveals an evolutionarily optimized three-layer regulation mechanism, which includes zonal segregation, epigenetic barrier crossing coupled to a negative feedback loop that mechanistically differs from previous theoretical proposals, and a previously unidentified enhancer competition step. This model not only recapitulates monoallelic OR expression, but also elucidates how the olfactory system maximizes and maintains the diversity of OR expression, and has multiple predictions validated by existing experimental results. Through making analogy to a physical system with thermally activated barrier crossing and comparative reverse engineering analyses, the study reveals that the olfactory receptor selection system is optimally designed, and particularly underscores cooperativity and synergy as a general design principle for multi-objective optimization in biology.
Proceedings of the National Academy of Sciences, 2016
Multiple-objective optimization is common in biological systems. In the mammalian olfactory syste... more Multiple-objective optimization is common in biological systems. In the mammalian olfactory system, each sensory neuron stochastically expresses only one out of up to thousands of olfactory receptor (OR) gene alleles; at the organism level, the types of expressed ORs need to be maximized. Existing models focus only on monoallele activation, and cannot explain recent observations in mutants, especially the reduced global diversity of expressed ORs in G9a/GLP knockouts. In this work we integrated existing information on OR expression, and constructed a comprehensive model that has all its components based on physical interactions. Analyzing the model reveals an evolutionarily optimized three-layer regulation mechanism, which includes zonal segregation, epigenetic barrier crossing coupled to a negative feedback loop that mechanistically differs from previous theoretical proposals, and a previously unidentified enhancer competition step. This model not only recapitulates monoallelic OR ...
How do mammalian cells that share the same genome exist in notably distinct phenotypes, exhibitin... more How do mammalian cells that share the same genome exist in notably distinct phenotypes, exhibiting differences in morphology, gene expression patterns, and epigenetic chromatin statuses? Furthermore how do cells of different phenotypes differentiate reproducibly from a single fertilized egg? These are fundamental problems in developmental biology. Epigenetic histone modifications play an important role in the maintenance of different cell phenotypes. The exact molecular mechanism for inheritance of the modification patterns over cell generations remains elusive. The complexity comes partly from the number of molecular species and the broad time scales involved. In recent years mathematical modeling has made significant contributions on elucidating the molecular mechanisms of DNA methylation and histone Book chapter in Epigenetic Technological Applications (Elsevier), in press 2 covalent modification inheritance. We will pedagogically introduce the typical procedure and some technical details of performing a mathematical modeling study, and discuss future developments.
Experimental and computational analyses identify multiple cell populations during the epithelial-... more Experimental and computational analyses identify multiple cell populations during the epithelial-to-mesenchymal transition.
A focused theme in systems biology is to uncover design principles of biological networks, that i... more A focused theme in systems biology is to uncover design principles of biological networks, that is, how specific network structures yield specific systems properties. For this purpose, we have previously developed a reverse engineering procedure to identify network topologies with high likelihood in generating desired systems properties. Our method searches the continuous parameter space of an assembly of network topologies, without enumerating individual network topologies separately as traditionally done in other reverse engineering procedures. Here we tested this CPSS (continuous parameter space search) method on a previously studied problem: the resettable bistability of an Rb-E2F gene network in regulating the quiescence-to-proliferation transition of mammalian cells. From a simplified Rb-E2F gene network, we identified network topologies responsible for generating resettable bistability. The CPSS-identified topologies are consistent with those reported in the previous study based on individual topology search (ITS), demonstrating the effectiveness of the CPSS approach. Since the CPSS and ITS searches are based on different mathematical formulations and different algorithms, the consistency of the results also helps cross-validate both approaches. A unique advantage of the CPSS approach lies in its applicability to biological networks with large numbers of nodes. To aid the application of the CPSS approach to the study of other biological systems, we have developed a computer package that is available in Information S1.
CRISPR-based gene knock-in at endogenous sites is desirable in multiple fields such as quantitati... more CRISPR-based gene knock-in at endogenous sites is desirable in multiple fields such as quantitative studies of signal transduction pathways and gene regulation, synthetic biology, and disease modeling. Contrasting the knock-out procedure, a key step of CRISPR knock-in procedure relies on the homology-directed repairing (HDR) process that requires a donor construct as repair template. Therefore, it is desirable to generate a series of donor DNA constructs efficiently and cost-effectively. In this study, we developed a general Gibson assembly procedure that combines strengths of a Modular Overlap-Directed Assembly with Linkers (MODAL) strategy and a restriction enzyme based hierarchical framework. This procedure also allows fusing sgRNAs to the constructs for enhanced homology-directed repairing efficiency. Experimental tests on multiple constructs achieved from 3-8 folds of increase in assembly efficiency to high yield of constructs that failed to make with conventional Gibson assembly. The modularized procedure is simple, fast and cost-effective while making multiple constructs, and a computer package is provided for customized design.
Targeting microenvironmental factors that foster migratory cell phenotypes is a promising strateg... more Targeting microenvironmental factors that foster migratory cell phenotypes is a promising strategy for halting tumor migration. However, lack of mechanistic understanding of the process impedes pharmaceutical drug development. Using a novel 3D microtumor model with tight control over tumor size, we recapitulated tumor size-induced hypoxic microenvironment and emergence of migratory phenotypes in epithelial T47D breast microtumors as well as those of patient-derived primary metastatic breast cancer cells, mesothelioma cells and lung cancer xenograft cells (PDX). The microtumor models from various patient-derived tumor cells and PDX cells revealed upregulation of tumor secretome, matrix metalloproteinase-9 (MMP9), fibronectin (FN), and soluble E-cadherin (sE-CAD) consistent with the clinically reported elevated levels of FN and MMP9 in the patient breast tumors compared to healthy mammary gland. We further showed that the tumor secretome induces migratory phenotype in non-hypoxic, non-migratory small microtumors. Subsequent mathematical model analysis identified a two-stage microtumor progression and migration mechanism, i.e., hypoxia induces migratory phenotype in the early initialization stage, which then becomes self-sustained through positive feedback loop established among the secretome. Both computational and experimental studies showed that inhibition of tumor secretome effectively halts microtumor migration despite tumor heterogeneity, while inhibition of the hypoxia is effective only within a time window and is compromised by tumor-to-tumor variation of the growth dynamics, supporting our notion that hypoxia initiates migratory phenotypes but does not sustain it. In summary, we show that targeting temporal dynamics of evolving microenvironments during tumor progression can halt and bypass major hurdle of tumor heterogeneity.
11 The detection and transmission of the strength and temporal quality of intracellular and 12 ex... more 11 The detection and transmission of the strength and temporal quality of intracellular and 12 extracellullar signals is an essential cellular mechanism. While TGF-β signaling is one of the 13 most thoroughly studied signaling pathways, the mechanisms by which cells translate TGF-β 14 signals remain unclear. In this paper, through an integrated quantitative and computational 15 approach we demonstrate that crosstalk among multiple TGF-β activated pathways forms a relay 16 from SMAD to GLI1 that initializes and maintains SNAILl expression, respectively. This 17 transaction is smoothed and accelerated by another temporal switch from elevated cytosolic 18 GSK3 enzymatic activity to reduced nuclear GSK3 enzymatic activity. This nested relay 19 mechanism places SNAIL1 as a key integrator of information from TGF-β signaling 20 subsequently distributed through divergent pathways; essentially cells generate a transient or 21 peer-reviewed)
Uploads
Papers by Jianhua Xing