Inference of cellular level signaling networks using single-cell gene expression data in Caenorhabditis elegans reveals mechanisms of cell fate specification

Cell fate specification plays a key role to generate distinct cell types during metazoan development. However, most of the underlying signaling networks at cellular level are not well understood. Availability of time lapse single-cell gene expression data collected throughout Caenorhabditis elegans embryogenesis provides an excellent opportunity for investigating signaling networks underlying cell fate specification at systems, cellular and molecular levels. We propose a framework to infer signaling networks at cellular level by exploring the single-cell gene expression data. Through analyzing the expression data of nhr-25 , a hypodermis-specific transcription factor, in every cells of both wild-type and mutant C.elegans embryos through RNAi against 55 genes, we have inferred a total of 23 genes that regulate (activate or inhibit) nhr-25 expression in cell-specific fashion. We also infer the signaling pathways consisting of each of these genes and nhr-25 based on a probabilistic graphical model for the selected five founder cells, 'ABarp', 'ABpla', 'ABpra', 'Caa' and 'Cpa', which express nhr-25 and mostly develop into hypodermis. By integrating the inferred pathways, we reconstruct five signaling networks with one each for the five founder cells. Using RNAi gene knockdown as a validation method, the inferred networks are able to predict the effects of the knockdown genes. These signaling networks in the five founder cells are likely to ensure faithful hypodermis cell fate specification in C.elegans at cellular level. All source codes and data are available at the github repository https://github.com/xthuang226/Worm_Single_Cell_Data_and_Codes.git . zhuyuan@cug.edu.cn. Supplementary data are available at Bioinformatics online.