Fig. 7

Correlated gene networks form modules related to tissue type and mating status. a-b Correlations between WGCNA modules and sample traits including tissue (AMY, HT, or VP/NAc), mating status (Mated), collection time point as a continuous variable (Time) or individual post-mating time points (0.5 h, 2 h, or 12 h) for a prairie voles and b meadow voles. Color scale bar indicates Pearson r value for correlation. Red indicates positive r value and blue indicates negative r value. Each cell in matrix shows correlation p-value. See Supplementary Tables 8 & 9 for Pearson r values for significant correlations. c Matrix of distances between sample module membership scores (kME). Color bar indicates sample distances with darker shades indicating more similar samples. Dendrograms show hierarchical clustering of samples based on distance. d Matrix of correlations between sample module membership scores (kME). Color bar indicates sample correlations with reds indicating positively correlated samples and blues indicating negatively correlated samples. Dendrograms show hierarchical clustering of samples based on distance. Color legends in c and d indicate module species (black = prairie vole, white = meadow vole) and the brain region most positively correlated with the module (blue = AMY, green = HT, purple = VP/NAc, grey = no significant positive correlation with brain region)