FIGURE 2: Differential expression between species is more robust than between foraging environments. (A) Volcano plot of the pairwise comparison between MZ and TRC across environments, showing a large number of DEGs (n=5318) with roughly equal numbers of significantly upregulated genes between species (n=2667 MZ; n=2651 TRC). Given the nature of the comparison, genes considered upregulated in MZ are downregulated in TRC, and vice versa. Red indicates upregulated genes for MZ, blue depicts genes upregulated in TRC, while black represents genes that do not meet the significance threshold of <0.05 FDR. (B) Heatmap of the top 500 most variable genes in the RNA-seq dataset is shown. Individuals from the pelagic foraging treatment are labeled with gray triangles, while benthic individuals are labeled with black triangles. Species cluster together, but there is less obvious structuring by foraging environment, although MZ segregate by environment more so than TRC. Photographs courtesy of Ad Koning at Cichlid Press. (C) Enriched GO terms associated with genes upregulated in the MZ:TRC comparison across environments. Colors are representative of the -log10(p-value). Gray bars indicate no significance. Gene counts are given as a total for both species along with the corresponding GO term.