Single cell RNA-Seq analysis holds great promise for elucidating the networks and pathways controlling cellular differentiation and disease. However, the analysis of time series single cell RNA-Seq data raises several new computational challenges. Cells at each time point are often sampled from a mixture of cell types, each of which may be a progenitor of one, or several, specific fates making it hard to determine which cells should be used to reconstruct temporal trajectories. In addition, cells, even from the same time point, may be unsynchronized making it hard to rely on the measured time for determining these trajectories.
Carnegie Mellon University researchers present TASIC a new method for determining temporal trajectories, branching and cell assignments in single cell time series experiments. Unlike prior approaches TASIC uses on a probabilistic graphical model to integrate expression and time information making it more robust to noise and stochastic variations. Applying TASIC to in vitro myoblast differentiation and in-vivo lung development data the researchers show that it accurately reconstructs developmental trajectories from single cell experiments. The reconstructed models enabled them to identify key genes involved in cell fate determination and to obtain new insights about a specific type of lung cells and its role in development.
Graphical representation of an example HMM structure with 5 states
gs and ts correspond to the gene expression profiles and average time associated with state s. yi is the expression profile of cell Ci and vi is the temporal value of cell Ci. P(S5/S3) and P(S4/S3) are the transition probabilities. See text for a detailed description of the parameters of the model.
Availability – The TASIC software package is posted in the supporting website: https://www.andrew.cmu.edu/user/sabrinar/TASIC
Contact – [email protected]