NetAct: a computational platform to construct core transcription factor regulatory networks using gene activity

Signaling and transcriptional dynamics underlying early adaptation to oncogenic BRAF inhibition

A major contributor to poor sensitivity to anti-cancer kinase inhibitor therapy is drug-induced cellular adaptation, whereby remodeling of signaling and gene regulatory networks permits a drug-tolerant phenotype. Here, we resolve the scale and kinetics of critical subcellular events following oncogenic kinase inhibition and preceding cell cycle re-entry, using mass spectrometry-based phosphoproteomics and RNA sequencing (RNA-seq) to monitor the dynamics of thousands of growth- and survival-related signals over the first minutes, hours, and days of oncogenic BRAF inhibition in human melanoma cells. We observed sustained inhibition of the BRAF-ERK axis, gradual downregulation of cell cycle signaling, and three distinct, reversible phase transitions toward quiescence. Statistical inference of kinetically defined regulatory modules revealed a dominant compensatory induction of SRC family kinase (SFK) signaling, promoted in part by excess reactive oxygen species, rendering cells sensitive to co-treatment with an SFK inhibitor in vitro and in vivo, underscoring the translational potential for assessing early drug-induced adaptive signaling. A record of this paper's transparent peer review process is included in the supplemental information.

Multi-omic analyses reveal aberrant DNA methylation patterns and the associated biomarkers of nasopharyngeal carcinoma and its cancer stem cells

Aberrant DNA methylation is a hallmark of nasopharyngeal carcinoma (NPC) pathogenesis. The aberrant DNA methylation patterns in NPC, particularly in its cancer stem cells (CSCs), and their underlying significance require further elucidation. We integratively performed DNA methylome and transcriptome combined with single-nucleus RNA sequencing to investigate DNA methylation and gene expression patterns of NPC and CSCs. Unlike Epstein-Barr virus (EBV)-negative cells, NPC and CSCs harboring EBV displayed global DNA hypermethylation and they were more oncogenic and immunosuppressive. By correlating DNA methylation and gene expression profiles, we disclosed potential relationships between aberrant DNA methylation, tumorigenesis, metastasis, immunotherapy response, and radiotherapy resistance of NPC. After validating with datasets from GEO and TCGA, we identified aberrant DNA methylation-associated biomarkers including 9 NPC-specific diagnostic markers that had significantly higher DNA methylation levels in NPC than in normal tissues and 8 types of cancers, and 12 potential prognostic markers that were highly correlated to cell cycle dysregulation. Notably, 2 of these potential biomarkers highly expressed in CSCs were validated at the single-cell level. Our study not only identified new potential diagnostic and prognostic biomarkers but also provided new insight into aberrant DNA methylation-associated pathogenesis of NPC, which is beneficial for the development of precision diagnosis and treatment schemes.

Predicting TF-Target Gene Association Using a Heterogeneous Network and Enhanced Negative Sampling

Identifying interactions between transcription factors (TFs) and target genes is crucial for understanding the molecular mechanisms involved in biological processes and diseases. Traditional biological experiments used to determine these interactions are often time-consuming, costly, and limited in scale. Current computational methods mainly predict binding sites rather than direct interactions. Although recent studies have achieved high performance in predicting TF-target gene associations, they still face a significant challenge related to constructing a robust dataset of positive and negative samples. Currently, methods do not adequately focus on selecting negative samples, resulting in incomplete coverage of potential TF-target gene relationships. This article proposes a method to select enhanced negative samples to improve the prediction performance of TF-target gene interactions. Experimental results show that the proposed method achieves an average area under the curve (AUC) value of 0.9024 ± 0.0008 through 5-fold cross-validation. These results demonstrate the model's high efficiency and accuracy, confirming its potential application in predicting TF-target gene interactions across various datasets and paving the way for large-scale biomedical research.

MorPhiC Consortium: towards functional characterization of all human genes

Recent advances in functional genomics and human cellular models have substantially enhanced our understanding of the structure and regulation of the human genome. However, our grasp of the molecular functions of human genes remains incomplete and biased towards specific gene classes. The Molecular Phenotypes of Null Alleles in Cells (MorPhiC) Consortium aims to address this gap by creating a comprehensive catalogue of the molecular and cellular phenotypes associated with null alleles of all human genes using in vitro multicellular systems. In this Perspective, we present the strategic vision of the MorPhiC Consortium and discuss various strategies for generating null alleles, as well as the challenges involved. We describe the cellular models and scalable phenotypic readouts that will be used in the consortium's initial phase, focusing on 1,000 protein-coding genes. The resulting molecular and cellular data will be compiled into a catalogue of null-allele phenotypes. The methodologies developed in this phase will establish best practices for extending these approaches to all human protein-coding genes. The resources generated-including engineered cell lines, plasmids, phenotypic data, genomic information and computational tools-will be made available to the broader research community to facilitate deeper insights into human gene functions.

Global transcription machinery engineering in Yarrowia lipolytica

Global transcription machinery engineering (gTME) is a strategy for optimizing complex phenotypes in microbes by manipulating transcription factors (TFs) and their downstream transcriptional regulatory networks (TRN). In principle, gTME leads to a focused but comprehensive optimization of a microbe, also enabling the engineering of nonpathway functionalities, like stress resistance, protein expression, or growth rate. A link between a TF and a desired phenotype is to be established for a rationally designed gTME. For use in a high-throughput format with extensive libraries of TRN-engineered clones tested under multiple conditions, well-developed culturing and analytical protocols are needed, to reveal the pleiotropic effects of the TFs. This mini-review summarizes the gTME strategies and TFs described under different contexts in Yarrowia lipolytica. The outcomes of the gTME strategy application are also addressed, demonstrating its effectiveness in engineering complex, industrially relevant traits in Y. lipolytica.

Endothelial Dysfunction in the Tubule Area Accelerates the Progression of Early Diabetic Kidney Disease

Diabetic kidney disease (DKD) is the leading cause of end-stage renal disease. Therefore, understanding the molecular regulatory mechanisms underlying the pathogenesis of DKD is imperative. In this study, we aimed to explore the molecular mechanisms of tubule region endothelial dysfunction in early DKD. Early-stage DKD model was established in 16-week-old female db/db mice for 16 weeks. Body weight, glucose level, and urine albumin-to-creatinine ratio (UACR) were measured. Hematoxylin and eosin (H&E) and periodic acid-Schiff (PAS) staining were performed to evaluate pathological lesions. RNA sequencing data of the kidneys and integrated publicly available single-cell and spatial transcriptome datasets were used to investigate the mechanism of endothelial dysfunction. There was a significant increase in body weight (p = 0.001), glucose levels (p=0.0008), and UACR (p=0.006) in db/db mice compared with db/m mice. H&E and PAS staining showed that vacuolar lesions and protein casts of tubules were the major histopathological changes observed in early-stage DKD mice. The apoptotic pathway in endothelial cells was notably activated in DKD, and Thbs1 was identified as the central gene involved in this apoptotic process. Deconvolution of the cell composition in the RNA sequencing data showed a decrease in the proportion of endothelial cells in the DKD mice. Further analysis of the activity and regulatory network of transcription factors showed that Creb1 was activated in both mouse and human early-stage DKD, suggesting that Creb1 activation may be involved in early kidney injury. The endothelial cell apoptotic pathway is activated in DKD, and the proportion of endothelial cells was reduced in the DKD mice, which is significantly associated with Thbs1. Keywords: Diabetic kidney disease, Endothelial dysfunction, RNA sequencing,Thbs1, Creb1.

OneSC: a computational platform for recapitulating cell state transitions

MOTIVATION

Computational modeling of cell state transitions has been a great interest of many in the field of developmental biology, cancer biology, and cell fate engineering because it enables performing perturbation experiments in silico more rapidly and cheaply than could be achieved in a lab. Recent advancements in single-cell RNA-sequencing (scRNA-seq) allow the capture of high-resolution snapshots of cell states as they transition along temporal trajectories. Using these high-throughput datasets, we can train computational models to generate in silico "synthetic" cells that faithfully mimic the temporal trajectories.

RESULTS

Here we present OneSC, a platform that can simulate cell state transitions using systems of stochastic differential equations govern by a regulatory network of core transcription factors (TFs). Different from many current network inference methods, OneSC prioritizes on generating Boolean network that produces faithful cell state transitions and terminal cell states that mimic real biological systems. Applying OneSC to real data, we inferred a core TF network using a mouse myeloid progenitor scRNA-seq dataset and showed that the dynamical simulations of that network generate synthetic single-cell expression profiles that faithfully recapitulate the four myeloid differentiation trajectories going into differentiated cell states (erythrocytes, megakaryocytes, granulocytes, and monocytes). Finally, through the in silico perturbations of the mouse myeloid progenitor core network, we showed that OneSC can accurately predict cell fate decision biases of TF perturbations that closely match with previous experimental observations.

AVAILABILITY AND IMPLEMENTATION

OneSC is implemented as a Python package on GitHub (https://github.com/CahanLab/oneSC) and on Zenodo (https://zenodo.org/records/14052421).

Network modeling links kidney developmental programs and the cancer type-specificity of VHL mutations

Elucidating the molecular dependencies behind the cancer-type specificity of driver mutations may reveal new therapeutic opportunities. We hypothesized that developmental programs would impact the transduction of oncogenic signaling activated by a driver mutation and shape its cancer-type specificity. Therefore, we designed a computational analysis framework by combining single-cell gene expression profiles during fetal organ development, latent factor discovery, and information theory-based differential network analysis to systematically identify transcription factors that selectively respond to driver mutations under the influence of organ-specific developmental programs. After applying this approach to VHL mutations, which are highly specific to clear cell renal cell carcinoma (ccRCC), we revealed important regulators downstream of VHL mutations in ccRCC and used their activities to cluster patients with ccRCC into three subtypes. This classification revealed a more significant difference in prognosis than the previous mRNA profile-based method and was validated in an independent cohort. Moreover, we found that EP300, a key epigenetic factor maintaining the regulatory network of the subtype with the worst prognosis, can be targeted by a small inhibitor, suggesting a potential treatment option for a subset of patients with ccRCC. This work demonstrated an intimate relationship between organ development and oncogenesis from the perspective of systems biology, and the method can be generalized to study the influence of other biological processes on cancer driver mutations.

Dissecting reversible and irreversible single cell state transitions from gene regulatory networks

Understanding cell state transitions and their governing regulatory mechanisms remains one of the fundamental questions in biology. We develop a computational method, state transition inference using cross-cell correlations (STICCC), for predicting reversible and irreversible cell state transitions at single-cell resolution by using gene expression data and a set of gene regulatory interactions. The method is inspired by the fact that the gene expression time delays between regulators and targets can be exploited to infer past and future gene expression states. From applications to both simulated and experimental single-cell gene expression data, we show that STICCC-inferred vector fields capture basins of attraction and irreversible fluxes. By connecting regulatory information with systems' dynamical behaviors, STICCC reveals how network interactions influence reversible and irreversible state transitions. Compared to existing methods that infer pseudotime and RNA velocity, STICCC provides complementary insights into the gene regulation of cell state transitions.

An Interplay between Transcription Factors and Recombinant Protein Synthesis in Yarrowia lipolytica at Transcriptional and Functional Levels-The Global View

Transcriptional regulatory networks (TRNs) associated with recombinant protein (rProt) synthesis in Yarrowia lipolytica are still under-described. Yet, it is foreseen that skillful manipulation with TRNs would enable global fine-tuning of the host strain's metabolism towards a high-level-producing phenotype. Our previous studies investigated the transcriptomes of Y. lipolytica strains overproducing biochemically different rProts and the functional impact of transcription factors (TFs) overexpression (OE) on rProt synthesis capacity in this species. Hence, much knowledge has been accumulated and deposited in public repositories. In this study, we combined both biological datasets and enriched them with further experimental data to investigate an interplay between TFs and rProts synthesis in Y. lipolytica at transcriptional and functional levels. Technically, the RNAseq datasets were extracted and re-analyzed for the TFs' expression profiles. Of the 140 TFs in Y. lipolytica, 87 TF-encoding genes were significantly deregulated in at least one of the strains. The expression profiles were juxtaposed against the rProt amounts from 125 strains co-overexpressing TF and rProt. In addition, several strains bearing knock-outs (KOs) in the TF loci were analyzed to get more insight into their actual involvement in rProt synthesis. Different profiles of the TFs' transcriptional deregulation and the impact of their OE or KO on rProts synthesis were observed, and new engineering targets were pointed.

Assessing next-generation sequencing-based computational methods for predicting transcriptional regulators with query gene sets

This article provides an in-depth review of computational methods for predicting transcriptional regulators (TRs) with query gene sets. Identification of TRs is of utmost importance in many biological applications, including but not limited to elucidating biological development mechanisms, identifying key disease genes, and predicting therapeutic targets. Various computational methods based on next-generation sequencing (NGS) data have been developed in the past decade, yet no systematic evaluation of NGS-based methods has been offered. We classified these methods into two categories based on shared characteristics, namely library-based and region-based methods. We further conducted benchmark studies to evaluate the accuracy, sensitivity, coverage, and usability of NGS-based methods with molecular experimental datasets. Results show that BART, ChIP-Atlas, and Lisa have relatively better performance. Besides, we point out the limitations of NGS-based methods and explore potential directions for further improvement.

Computational Reconstruction of the Transcription Factor Regulatory Network Induced by Auxin in Arabidopsis thaliana L

In plant hormone signaling, transcription factor regulatory networks (TFRNs), which link the master transcription factors to the biological processes under their control, remain insufficiently characterized despite their crucial function. Here, we identify a TFRN involved in the response to the key plant hormone auxin and define its impact on auxin-driven biological processes. To reconstruct the TFRN, we developed a three-step procedure, which is based on the integrated analysis of differentially expressed gene lists and a representative collection of transcription factor binding profiles. Its implementation is available as a part of the CisCross web server. With the new method, we distinguished two transcription factor subnetworks. The first operates before auxin treatment and is switched off upon hormone application, the second is switched on by the hormone. Moreover, we characterized the functioning of the auxin-regulated TFRN in control of chlorophyll and lignin biosynthesis, abscisic acid signaling, and ribosome biogenesis.

OneSC: A computational platform for recapitulating cell state transitions

Computational modelling of cell state transitions has been a great interest of many in the field of developmental biology, cancer biology and cell fate engineering because it enables performing perturbation experiments in silico more rapidly and cheaply than could be achieved in a wet lab. Recent advancements in single-cell RNA sequencing (scRNA-seq) allow the capture of high-resolution snapshots of cell states as they transition along temporal trajectories. Using these high-throughput datasets, we can train computational models to generate in silico 'synthetic' cells that faithfully mimic the temporal trajectories. Here we present OneSC, a platform that can simulate synthetic cells across developmental trajectories using systems of stochastic differential equations govern by a core transcription factors (TFs) regulatory network. Different from the current network inference methods, OneSC prioritizes on generating Boolean network that produces faithful cell state transitions and steady cell states that mimic real biological systems. Applying OneSC to real data, we inferred a core TF network using a mouse myeloid progenitor scRNA-seq dataset and showed that the dynamical simulations of that network generate synthetic single-cell expression profiles that faithfully recapitulate the four myeloid differentiation trajectories going into differentiated cell states (erythrocytes, megakaryocytes, granulocytes and monocytes). Finally, through the in-silico perturbations of the mouse myeloid progenitor core network, we showed that OneSC can accurately predict cell fate decision biases of TF perturbations that closely match with previous experimental observations.

Data-driven modeling of core gene regulatory network underlying leukemogenesis in IDH mutant AML

Acute myeloid leukemia (AML) is characterized by uncontrolled proliferation of poorly differentiated myeloid cells, with a heterogenous mutational landscape. Mutations in IDH1 and IDH2 are found in 20% of the AML cases. Although much effort has been made to identify genes associated with leukemogenesis, the regulatory mechanism of AML state transition is still not fully understood. To alleviate this issue, here we develop a new computational approach that integrates genomic data from diverse sources, including gene expression and ATAC-seq datasets, curated gene regulatory interaction databases, and mathematical modeling to establish models of context-specific core gene regulatory networks (GRNs) for a mechanistic understanding of tumorigenesis of AML with IDH mutations. The approach adopts a new optimization procedure to identify the top network according to its accuracy in capturing gene expression states and its flexibility to allow sufficient control of state transitions. From GRN modeling, we identify key regulators associated with the function of IDH mutations, such as DNA methyltransferase DNMT1, and network destabilizers, such as E2F1. The constructed core regulatory network and outcomes of in-silico network perturbations are supported by survival data from AML patients. We expect that the combined bioinformatics and systems-biology modeling approach will be generally applicable to elucidate the gene regulation of disease progression.

Assessing NGS-based computational methods for predicting transcriptional regulators with query gene sets

This article provides an in-depth review of computational methods for predicting transcriptional regulators with query gene sets. Identification of transcriptional regulators is of utmost importance in many biological applications, including but not limited to elucidating biological development mechanisms, identifying key disease genes, and predicting therapeutic targets. Various computational methods based on next-generation sequencing (NGS) data have been developed in the past decade, yet no systematic evaluation of NGS-based methods has been offered. We classified these methods into two categories based on shared characteristics, namely library-based and region-based methods. We further conducted benchmark studies to evaluate the accuracy, sensitivity, coverage, and usability of NGS-based methods with molecular experimental datasets. Results show that BART, ChIP-Atlas, and Lisa have relatively better performance. Besides, we point out the limitations of NGS-based methods and explore potential directions for further improvement.

Key points

An introduction to available computational methods for predicting functional TRs from a query gene set.A detailed walk-through along with practical concerns and limitations.A systematic benchmark of NGS-based methods in terms of accuracy, sensitivity, coverage, and usability, using 570 TR perturbation-derived gene sets.NGS-based methods outperform motif-based methods. Among NGS methods, those utilizing larger databases and adopting region-centric approaches demonstrate favorable performance. BART, ChIP-Atlas, and Lisa are recommended as these methods have overall better performance in evaluated scenarios.

Characterizing neuroinflammation and identifying prenatal diagnostic markers for neural tube defects through integrated multi-omics analysis

BACKGROUND

Neural Tube Defects (NTDs) are congenital malformations of the central nervous system resulting from the incomplete closure of the neural tube during early embryonic development. Neuroinflammation refers to the inflammatory response in the nervous system, typically resulting from damage to neural tissue. Immune-related processes have been identified in NTDs, however, the detailed relationship and underlying mechanisms between neuroinflammation and NTDs remain largely unclear. In this study, we utilized integrated multi-omics analysis to explore the role of neuroinflammation in NTDs and identify potential prenatal diagnostic markers using a murine model.

METHODS

Nine public datasets from Gene Expression Omnibus (GEO) and ArrayExpress were mined using integrated multi-omics analysis to characterize the molecular landscape associated with neuroinflammation in NTDs. Special attention was given to the involvement of macrophages in neuroinflammation within amniotic fluid, as well as the dynamics of macrophage polarization and their interactions with neural cells at single-cell resolution. We also used qPCR assay to validate the key TFs and candidate prenatal diagnostic genes identified through the integrated analysis in a retinoic acid-induced NTDs mouse model.

RESULTS

Our analysis indicated that neuroinflammation is a critical pathological feature of NTDs, regulated both transcriptionally and epigenetically within central nervous system tissues. Key alterations in gene expression and pathways highlighted the crucial role of STATs molecules in the JAK-STAT signaling pathway in regulating NTDs-associated neuroinflammation. Furthermore, single-cell resolution analysis revealed significant polarization of macrophages and their interaction with neural cells in amniotic fluid, underscoring their central role in mediating neuroinflammation associated with NTDs. Finally, we identified a set of six potential prenatal diagnostic genes, including FABP7, CRMP1, SCG3, SLC16A10, RNASE6 and RNASE1, which were subsequently validated in a murine NTDs model, indicating their promise as prospective markers for prenatal diagnosis of NTDs.

CONCLUSIONS

Our study emphasizes the pivotal role of neuroinflammation in the progression of NTDs and underlines the potential of specific inflammatory and neural markers as novel prenatal diagnostic tools. These findings provide important clues for further understanding the underlying mechanisms between neuroinflammation and NTDs, and offer valuable insights for the future development of prenatal diagnostics.

Single-cell transcriptome analysis reveals the association between histone lactylation and cisplatin resistance in bladder cancer

Patients with bladder cancer (BCa) frequently acquires resistance to platinum-based chemotherapy, particularly cisplatin. This study centered on the mechanism of cisplatin resistance in BCa and highlighted the pivotal role of lactylation in driving this phenomenon. Utilizing single-cell RNA sequencing, we delineated the single-cell landscape of Bca, pinpointing a distinctive subset of BCa cells that exhibit marked resistance to cisplatin with association with glycolysis metabolism. Notably, we observed that H3 lysine 18 lactylation (H3K18la) plays a crucial role in activating the transcription of target genes by enriching in their promoter regions. Targeted inhibition of H3K18la effectively restored cisplatin sensitivity in these cisplatin-resistant epithelial cells. Furthermore, H3K18la-driven key transcription factors YBX1 and YY1 promote cisplatin resistance in BCa. These findings enhance our understanding of the mechanisms underlying cisplatin resistance, offering valuable insights for identifying novel intervention targets to overcome drug resistance in Bca.

Signaling and transcriptional dynamics underlying early adaptation to oncogenic BRAF inhibition

A major contributor to poor sensitivity to anti-cancer kinase inhibitor therapy is drug-induced cellular adaptation, whereby remodeling of signaling and gene regulatory networks permits a drug-tolerant phenotype. Here, we resolve the scale and kinetics of critical subcellular events following oncogenic kinase inhibition and preceding cell cycle re-entry, using mass spectrometry-based phosphoproteomics and RNA sequencing to capture molecular snapshots within the first minutes, hours, and days of BRAF kinase inhibitor exposure in a human BRAF -mutant melanoma model of adaptive therapy resistance. By enriching specific phospho-motifs associated with mitogenic kinase activity, we monitored the dynamics of thousands of growth- and survival-related protein phosphorylation events under oncogenic BRAF inhibition and drug removal. We observed early and sustained inhibition of the BRAF-ERK axis, gradual downregulation of canonical cell cycle-dependent signals, and three distinct and reversible phase transitions toward quiescence. Statistical inference of kinetically-defined signaling and transcriptional modules revealed a concerted response to oncogenic BRAF inhibition and a dominant compensatory induction of SRC family kinase (SFK) signaling, which we found to be at least partially driven by accumulation of reactive oxygen species via impaired redox homeostasis. This induction sensitized cells to co-treatment with an SFK inhibitor across a panel of patient-derived melanoma cell lines and in an orthotopic mouse xenograft model, underscoring the translational potential for measuring the early temporal dynamics of signaling and transcriptional networks under therapeutic challenge.

Transcriptional state dynamics lead to heterogeneity and adaptive tumor evolution in urothelial bladder carcinoma

Intra-tumor heterogeneity contributes to treatment failure and poor survival in urothelial bladder carcinoma (UBC). Analyzing transcriptome from a UBC cohort, we report that intra-tumor transcriptomic heterogeneity indicates co-existence of tumor cells in epithelial and mesenchymal-like transcriptional states and bi-directional transition between them occurs within and between tumor subclones. We model spontaneous and reversible transition between these partially heritable states in cell lines and characterize their population dynamics. SMAD3, KLF4 and PPARG emerge as key regulatory markers of the transcriptional dynamics. Nutrient limitation, as in the core of large tumors, and radiation treatment perturb the dynamics, initially selecting for a transiently resistant phenotype and then reconstituting heterogeneity and growth potential, driving adaptive evolution. Dominance of transcriptional states with low PPARG expression indicates an aggressive phenotype in UBC patients. We propose that phenotypic plasticity and dynamic, non-genetic intra-tumor heterogeneity modulate both the trajectory of disease progression and adaptive treatment response in UBC.

Identifying key factors in cell fate decisions by machine learning interpretable strategies

Cell fate decisions and transitions are common in almost all developmental processes. Therefore, it is important to identify the decision-making mechanisms and important individual molecules behind the fate decision processes. In this paper, we propose an interpretable strategy based on systematic perturbation, unsupervised hierarchical cluster analysis (HCA), machine learning (ML), and Shapley additive explanation (SHAP) analysis for inferring the contribution and importance of individual molecules in cell fate decision and transition processes. In order to verify feasibility of the approach, we apply it to the core epithelial to mesenchymal transition (EMT)-metastasis network. The key factors identified in EMT-metastasis are consistent with relevant experimental observations. The approach presented here can be applied to other biological networks to identify important factors related to cell fate decisions and transitions.

Quantifying cancer cell plasticity with gene regulatory networks and single-cell dynamics

Phenotypic plasticity of cancer cells can lead to complex cell state dynamics during tumor progression and acquired resistance. Highly plastic stem-like states may be inherently drug-resistant. Moreover, cell state dynamics in response to therapy allow a tumor to evade treatment. In both scenarios, quantifying plasticity is essential for identifying high-plasticity states or elucidating transition paths between states. Currently, methods to quantify plasticity tend to focus on 1) quantification of quasi-potential based on the underlying gene regulatory network dynamics of the system; or 2) inference of cell potency based on trajectory inference or lineage tracing in single-cell dynamics. Here, we explore both of these approaches and associated computational tools. We then discuss implications of each approach to plasticity metrics, and relevance to cancer treatment strategies.

scTIGER: A Deep-Learning Method for Inferring Gene Regulatory Networks from Case versus Control scRNA-seq Datasets

Inferring gene regulatory networks (GRNs) from single-cell RNA-seq (scRNA-seq) data is an important computational question to find regulatory mechanisms involved in fundamental cellular processes. Although many computational methods have been designed to predict GRNs from scRNA-seq data, they usually have high false positive rates and none infer GRNs by directly using the paired datasets of case-versus-control experiments. Here we present a novel deep-learning-based method, named scTIGER, for GRN detection by using the co-differential relationships of gene expression profiles in paired scRNA-seq datasets. scTIGER employs cell-type-based pseudotiming, an attention-based convolutional neural network method and permutation-based significance testing for inferring GRNs among gene modules. As state-of-the-art applications, we first applied scTIGER to scRNA-seq datasets of prostate cancer cells, and successfully identified the dynamic regulatory networks of AR, ERG, PTEN and ATF3 for same-cell type between prostatic cancerous and normal conditions, and two-cell types within the prostatic cancerous environment. We then applied scTIGER to scRNA-seq data from neurons with and without fear memory and detected specific regulatory networks for BDNF, CREB1 and MAPK4. Additionally, scTIGER demonstrates robustness against high levels of dropout noise in scRNA-seq data.

Data-driven modeling of core gene regulatory network underlying leukemogenesis in IDH mutant AML

Acute myeloid leukemia (AML) is characterized by uncontrolled proliferation of poorly differentiated myeloid cells, with a heterogenous mutational landscape. Mutations in IDH1 and IDH2 are found in 20% of the AML cases. Although much effort has been made to identify genes associated with leukemogenesis, the regulatory mechanism of AML state transition is still not fully understood. To alleviate this issue, here we develop a new computational approach that integrates genomic data from diverse sources, including gene expression and ATAC-seq datasets, curated gene regulatory interaction databases, and mathematical modeling to establish models of context-specific core gene regulatory networks (GRNs) for a mechanistic understanding of tumorigenesis of AML with IDH mutations. The approach adopts a novel optimization procedure to identify the optimal network according to its accuracy in capturing gene expression states and its flexibility to allow sufficient control of state transitions. From GRN modeling, we identify key regulators associated with the function of IDH mutations, such as DNA methyltransferase DNMT1, and network destabilizers, such as E2F1. The constructed core regulatory network and outcomes of in-silico network perturbations are supported by survival data from AML patients. We expect that the combined bioinformatics and systems-biology modeling approach will be generally applicable to elucidate the gene regulation of disease progression.