IQCELL: A platform for predicting the effect of gene perturbations on developmental trajectories using single-cell RNA-seq data

Inferring Gene Regulatory Networks and Predicting the Effect of Gene Perturbations via IQCELL

IQCELL is a platform that infers Boolean gene regulatory networks from single-cell RNA sequencing data. Boolean networks can be simulated under normal and perturbed conditions. In this chapter, we provide a detailed guideline for implementing IQCELL from a raw dataset. The steps include processing data, inferring informative genes, inferring gene regulatory network, and simulating the resulted network under normal and perturbed conditions.

Normalizing Input-Output Relationships of Cancer Networks for Reversion Therapy

Accumulated genetic alterations in cancer cells distort cellular stimulus-response (or input-output) relationships, resulting in uncontrolled proliferation. However, the complex molecular interaction network within a cell implicates a possibility of restoring such distorted input-output relationships by rewiring the signal flow through controlling hidden molecular switches. Here, a system framework of analyzing cellular input-output relationships in consideration of various genetic alterations and identifying possible molecular switches that can normalize the distorted relationships based on Boolean network modeling and dynamics analysis is presented. Such reversion is demonstrated by the analysis of a number of cancer molecular networks together with a focused case study on bladder cancer with in vitro experiments and patient survival data analysis. The origin of reversibility from an evolutionary point of view based on the redundancy and robustness intrinsically embedded in complex molecular regulatory networks is further discussed.

Identification of transcriptional programs using dense vector representations defined by mutual information with GeneVector

Deciphering individual cell phenotypes from cell-specific transcriptional processes requires high dimensional single cell RNA sequencing. However, current dimensionality reduction methods aggregate sparse gene information across cells, without directly measuring the relationships that exist between genes. By performing dimensionality reduction with respect to gene co-expression, low-dimensional features can model these gene-specific relationships and leverage shared signal to overcome sparsity. We describe GeneVector, a scalable framework for dimensionality reduction implemented as a vector space model using mutual information between gene expression. Unlike other methods, including principal component analysis and variational autoencoders, GeneVector uses latent space arithmetic in a lower dimensional gene embedding to identify transcriptional programs and classify cell types. In this work, we show in four single cell RNA-seq datasets that GeneVector was able to capture phenotype-specific pathways, perform batch effect correction, interactively annotate cell types, and identify pathway variation with treatment over time.

Systems bioengineering approaches for developmental toxicology

Developmental toxicology is the field of study that examines the effects of chemical and physical agents on developing organisms. By using principles of systems biology and bioengineering, a systems bioengineering approach could be applied to study the complex interactions between developing organisms, the environment, and toxic agents. This approach would result in a holistic understanding of the effects of toxic agents on organisms, by considering the interactions between different biological systems and the impacts of toxicants on those interactions. It would be useful in identifying key biological pathways and mechanisms affected by toxic agents, as well as in the development of predictive models to assess potential risks of exposure to toxicants during development. In this review, we discuss the relevance of systems bioengineering to the field of developmental toxicity and provide up-to-date examples that illustrate the use of engineering principles for this application.

Transcriptional networks of transient cell states during human prefrontal cortex development

The human brain is divided into various anatomical regions that control and coordinate unique functions. The prefrontal cortex (PFC) is a large brain region that comprises a range of neuronal and non-neuronal cell types, sharing extensive interconnections with subcortical areas, and plays a critical role in cognition and memory. A timely appearance of distinct cell types through embryonic development is crucial for an anatomically perfect and functional brain. Direct tracing of cell fate development in the human brain is not possible, but single-cell transcriptome sequencing (scRNA-seq) datasets provide the opportunity to dissect cellular heterogeneity and its molecular regulators. Here, using scRNA-seq data of human PFC from fetal stages, we elucidate distinct transient cell states during PFC development and their underlying gene regulatory circuitry. We further identified that distinct intermediate cell states consist of specific gene regulatory modules essential to reach terminal fate using discrete developmental paths. Moreover, using in silico gene knock-out and over-expression analysis, we validated crucial gene regulatory components during the lineage specification of oligodendrocyte progenitor cells. Our study illustrates unique intermediate states and specific gene interaction networks that warrant further investigation for their functional contribution to typical brain development and discusses how this knowledge can be harvested for therapeutic intervention in challenging neurodevelopmental disorders.

Virtual cells in a virtual microenvironment recapitulate early development-like patterns in human pluripotent stem cell colonies

The mechanism by which morphogenetic signals engage the regulatory networks responsible for early embryonic tissue patterning is incompletely understood. Here, we developed a minimal gene regulatory network (GRN) model of human pluripotent stem cell (hPSC) lineage commitment and embedded it into "cellular" agents that respond to a dynamic morphogenetic signaling microenvironment. Simulations demonstrated that GRN wiring had significant non-intuitive effects on tissue pattern order, composition, and dynamics. Experimental perturbation of GRN connectivities supported model predictions and demonstrated the role of OCT4 as a master regulator of peri-gastrulation fates. Our so-called GARMEN strategy provides a multiscale computational platform to understand how single-cell-based regulatory interactions scale to tissue domains. This foundation provides new opportunities to simulate the impact of network motifs on normal and aberrant tissue development.

spliceJAC: transition genes and state-specific gene regulation from single-cell transcriptome data

Extracting dynamical information from single-cell transcriptomics is a novel task with the promise to advance our understanding of cell state transition and interactions between genes. Yet, theory-oriented, bottom-up approaches that consider differences among cell states are largely lacking. Here, we present spliceJAC, a method to quantify the multivariate mRNA splicing from single-cell RNA sequencing (scRNA-seq). spliceJAC utilizes the unspliced and spliced mRNA count matrices to constructs cell state-specific gene-gene regulatory interactions and applies stability analysis to predict putative driver genes critical to the transitions between cell states. By applying spliceJAC to biological systems including pancreas endothelium development and epithelial-mesenchymal transition (EMT) in A549 lung cancer cells, we predict genes that serve specific signaling roles in different cell states, recover important differentially expressed genes in agreement with pre-existing analysis, and predict new transition genes that are either exclusive or shared between different cell state transitions.

OUP accepted manuscript

Recent advances in direct cell reprogramming have made possible the conversion of one cell type to another cell type, offering a potential cell-based treatment to many major diseases. Despite much attention, substantial roadblocks remain including the inefficiency in the proportion of reprogrammed cells of current experiments, and the requirement of a significant amount of time and resources. To this end, several computational algorithms have been developed with the goal of guiding the hypotheses to be experimentally validated. These approaches can be broadly categorized into two main types: transcription factor identification methods which aim to identify candidate transcription factors for a desired cell conversion, and transcription factor perturbation methods which aim to simulate the effect of a transcription factor perturbation on a cell state. The transcription factor perturbation methods can be broken down into Boolean networks, dynamical systems and regression models. We summarize the contributions and limitations of each method and discuss the innovation that single cell technologies are bringing to these approaches and we provide a perspective on the future direction of this field.