Fig. 3

Gene expression modules constructed by the WGCNA. a the clustering dendrogram of genes expressed in the XX and XY gonads is shown. WGCNA identifies gene modules using the cutreeDynamic function, which detects clusters in a dendrogram depending on their shape [32]. Original modules of very similar coexpression are merged. Color name-based module labels are generated automatically by WGCNA and are shown along with the number of genes belonging to each module. b Association between modules and traits (such as Genotype, Phenotype, Developmental stages, and Key genes) is depicted. Rows correspond to module eigengenes (shown on the left with the colored boxes), columns to a trait (indicated below each column). Each cell contains the P value of the trait-module association. Large positive values indicate a strong correlation, while large negative values indicate a strong negative association