Experimental and computational analyses identify multiple cell populations during the epithelial-to-mesenchymal transition.

Epithelial-to-mesenchymal transition (EMT) is a key step in development, wound healing, and cancer development. It involves cooperation of signaling pathways, such as transformation growth factor-β (TGF-β), Sonic Hedgehog (SHH), and WNT pathways. These signaling pathways crosstalk to each other and converge to key transcription factors (e.g., SNAIL1) to initialize and maintain the process of EMT. The functional roles of multi-signaling pathway crosstalks in EMT are sophisticated and, thus, remain to be explored. In this review, we focused on three major signal transduction pathways that promote or regulate EMT in carcinoma. We discussed the network structures, and provided a brief overview of the current therapy strategies and drug development targeted to these three signal transduction pathways. Finally, we highlighted systems biology approaches that can accelerate the process of deconstructing complex networks and drug discovery.

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 chromosome structure of the gene. To examine how local gene environment and transcription factor regulation are coupled, we performed a combined analysis of time-course RNA-seq data of TGF-β treated MCF10A cells and related epigenomic and Hi-C data. Using Dynamic Regulatory Events Miner (DREM), we clustered differentially expressed genes based on gene expression profiles and associated transcription factors. Genes in each class have similar temporal gene expression patterns and share common transcription factors. Next, we defined a set of linear and radial distribution functions, as used in statistical physics, to measure the distributions of genes within a class both spatially and linearly along the genomic sequence. Remarkably, genes within the same class despite sometimes being separated by tens of million bases (Mb) along genomic sequence show a significantly higher tendency to be spatially close despite sometimes being separated by tens of Mb along the genomic sequence than those belonging to different classes do. Analyses extended to the process of mouse nervous system development arrived at similar conclusions. Future studies will be able to test whether this spatial organization of chromosomes contributes to concerted gene expression. 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 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, we performed an integrated quantitative and computational analysis on TGF-β induced activation of SNAIL1, a key transcription factor that regulates several subsequent cell fate decisions such as apoptosis and epithelial-to-mesenchymal transition. We demonstrate that crosstalk among multiple TGF-β activated pathways forms a relay from SMAD to GLI1 that initializes and maintains SNAILl expression, respectively. SNAIL1 functions as a key integrator of information from TGF-β signaling distributed through upstream divergent pathways. The intertwined network serves as a temporal checkpoint, so that cells can generate a transient or sustained expression of SNAIL1 depending on TGF-β duration. Furthermore, we observed that TGF-β treatment leads to an unexpected accumulation of GSK3 molecules in an enzymatically active tyrosine phosphorylation form in Golgi apparatus and ER, followed by accumulation of GSK3 molecules in an enzymatically inhibitive serine phosphorylation in the nucleus. Subsequent model analysis and inhibition experiments revealed that the initial localized increase of GSK3 enzymatic activity couples to the positive feedback loop of the substrate Gli1 to form a network motif with multi-objective functions. That is, the motif is robust against stochastic fluctuations, and has a narrow distribution of response time that is insensitive to initial conditions. Specifically for TGF-β signaling, the motif ensures a smooth relay from SMAD to GLI1 on regulating SNAIL1 expression.

Clustered regularly interspaced short palindromic repeats (CRISPR)-based gene editing techniques find applications in many fields, such as molecular biology, cancer biology, and disease modeling. In contrast to the knock-out procedure, a key step of CRISPR knock-in experiments is the homology-directed repair process that requires donor constructs as repair templates. Therefore, it is desirable to generate a series of donor templates efficiently and cost-effectively. In this study, we developed a new strategy that combines (i) Gibson assembly reaction, (ii) a linker pair composed of eight in silico screened restriction enzyme sites, and (iii) a hierarchical framework, to remarkably improve the efficiency of producing donor constructs for common genes as well as for the genes containing unbalanced guanine-cytosine content and requiring a selectable marker. Furthermore, the approach provides the ability of inserting additional elements into the donor templates, such as single guide RNA recognition sites that have been reported to enhance the efficiency of homology-directed repair. Conclusively, our modularized process is simple, fast, and cost-effective for making donor constructs and benefits the application of CRISPR knock-in methods.