Abstract:Recent technological advancements have made time-resolved, quantitative, multi-omics data available for many model systems, which could be integrated for systems pharmacokinetic use. Here, we present large-scale simulation modeling (LASSIM), which is a novel mathematical tool for performing large-scale inference using mechanistically defined ordinary differential equations (ODE) for gene regulatory networks (GRNs). LASSIM integrates structural knowledge about regulatory interactions and non-linear equations with multiple steady state and dynamic response expression datasets. The rationale behind LASSIM is that biological GRNs can be simplified using a limited subset of core genes that are assumed to regulate all other gene transcription events in the network. The LASSIM method is implemented as a general-purpose toolbox using the PyGMO Python package to make the most of multicore computers and high performance clusters, and is available at https://gitlab.com/Gustafsson-lab/lassim. As a method, LASSIM works in two steps, where it first infers a non-linear ODE system of the pre-specified core gene expression. Second, LASSIM in parallel optimizes the parameters that model the regulation of peripheral genes by core system genes. We showed the usefulness of this method by applying LASSIM to infer a large-scale non-linear model of naïve Th2 cell differentiation, made possible by integrating Th2 specific bindings, time-series together with six public and six novel siRNA-mediated knock-down experiments. ChIP-seq showed significant overlap for all tested transcription factors. Next, we performed novel time-series measurements of total T-cells during differentiation towards Th2 and verified that our LASSIM model could monitor those data significantly better than comparable models that used the same Th2 bindings. In summary, the LASSIM toolbox opens the door to a new type of model-based data analysis that combines the strengths of reliable mechanistic models with truly systems-level data. We demonstrate the power of this approach by inferring a mechanistically motivated, genome-wide model of the Th2 transcription regulatory system, which plays an important role in several immune related diseases.

GLMsim: a GLM-based single cell RNA-seq simulator incorporating batch and biological effects

SimCH: simulation of single-cell RNA sequencing data by modeling cellular heterogeneity at gene expression level

scReadSim: a single-cell RNA-seq and ATAC-seq read simulator

simCAS: an embedding-based method for simulating single-cell chromatin accessibility sequencing data

SIMLR: a Tool for Large-Scale Single-Cell Analysis by Multi-Kernel Learning.

The Severn Estuary and Bristol Channel: a 25 year critical review.

A sandbox for prediction and integration of DNA, RNA, and proteins in single cells

Assessment of batch-correction methods for scRNA-seq data with a new test metric

MOSim: bulk and single-cell multi-layer regulatory network simulator

Multi-Experiment Nonlinear Mixed Effect Modeling of Single-Cell Translation Kinetics after Transfection

NGSNGS: next-generation simulator for next-generation sequencing data

AsaruSim: a single-cell and spatial RNA-Seq Nanopore long-reads simulation workflow

SimCD: Simultaneous Clustering and Differential expression analysis for single-cell transcriptomic data

Multi-task benchmarking of spatially resolved gene expression simulation models

A benchmark of batch-effect correction methods for single-cell RNA sequencing data

Sensitive spectrofluorimetric determination of tetracycline residues in bovine milk.

SIMLR: A Tool for Large-Scale Genomic Analyses by Multi-Kernel Learning

An Automated Microwell Platform for Large-Scale Single Cell RNA-Seq

LTMG: a novel statistical modeling of transcriptional expression states in single-cell RNA-Seq data

The shaky foundations of simulating single-cell RNA sequencing data

LASSIM-A network inference toolbox for genome-wide mechanistic modeling