Fig. 2

Weighted correlation network analysis of the scRNA-seq data. A Plot of scale-free topology and mean connectivity with regards to soft-thresholding power for samples from scRNA-seq data. Red line indicates an R2 cut-off of 0.87. B Visualization of the eigengene network representing the relationships among the modules and the topology of genes. C Clustering dendrogram showing correlations between modules from scRNA-seq data. D A gene co-expression network of single cells from the limb bud was constructed by weighted correlation network analysis. The hierarchical clustering and module assignment of genes are shown along the left side and top. E Gene expression modules showing correlations in relation to cell type. Positive correlations are red while negative correlations are green. The module and the sample traits are related to the heat map, the row represents the module, and the column represents the trait. The values in the box indicate correlation and FDR. F Relative numbers of DEGs identified in each module. Black bars represent |correlation coefficients|> 0.6 and FDR < 0.05 modules, and gray bars represent all other modules. G The box plot shows the relative number of DEGs in significant (sig.) modules and not significant (not sig.) modules compared with the blue module (including Hoxd13)