Fig. 3
Differential gene expression analysis of 7dpf Danio rerio enteric neurons. a and b. Plot showing the correlation between the two biological replicates of the GFP-positive neuronal population (a) and of the GFP-negative non-neuronal population (b). The x-axis represents the 1st biological replicate while the y-axis represents the 2nd biological replicate. The correlation coefficient was calculated using FPKM. The correlation coefficient for the neuronal and non-neuronal cell population is 0.95 and 0.96, respectively. c Dispersion plot showing the log2foldchange of individual reads as compared to the mean normalized counts using DESeq with the adjusted p-value ≤ 0.01. The red dots represent the upregulated/enriched and the blue dots represent downregulated/depleted genes in the neuronal cells of the zebrafish intestine. d. Heat map of the normalized counts for the 4418 DEGs between neuronal and non-neuronal cell populations. The thresholds selected are FDR ≤ 0.01 and log2foldchange ≥ 1.5. The first two columns represent biological replicates (a and b) of the non-neuronal cell population and the second two columns represent biological replicates (a and b) of the neuronal population. Each row denotes a single gene and their expression pattern across the two different samples. The blue regions in the neuronal population represent the 57.99% of the upregulated genes and the red region represents the downregulated genes. The dendogram on the x-axis shows the grouping of same samples (biological replicates) while the y-axis shows clustering of (genes with) similar expression patterns