Hierarchical clusterization of genes differentially expressed during autogamy. The 2467 genes most differentially expressed during autogamy (treat model fold-change of 2, FDR < 0.05) were hierarchically clustered as described in Methods. The heatmap displays the samples as columns and the genes as rows. The color code goes from dark blue for the lowest expression to dark red for the highest expression. The clusterization of the samples corresponds to the classification based on nuclear morphology, although the DEV2 and DEV3 samples are not resolved. The 6 clusters of co-expressed genes were obtained by cutting the gene dendrogram on the left of the heatmap as indicated. On the right hand side, the average expression profile and standard deviation for each cluster are drawn as they appear in ParameciumDB .