Fig. 5From: Comparative neurotranscriptomics reveal widespread species differences associated with bondingVole brain gene expression varies over time by species. a-c Volcano plots showing significance (−log10(p-value)) and magnitude of difference (log2 fold change, LFC) in expression for each gene in species:time interaction in a AMY, b HT, and c VP/NAc. Each point represents a single gene. Orange colored points pass significance threshold of FDR < 0.1. Darker shaded points pass FDR cutoff and have LFC > 2. Genes with positive LFC values are enriched in prairie voles, genes with negative LFC values are enriched in meadow voles. d-f Pattern of expression over time for genes with lowest adjusted p-value (FDR) in each regional model. Plots show normalized counts for samples collected at each pre-mating (0 h) and post-mating (0.5, 2, or 12 h) collection point. Each point represents one sample. Darker shades and triangles represent prairie vole samples. Lighter shades and circles represent meadow vole samples. d Npas4 in AMY samples. e Aqp3 in HT samples. f Spata16 in VP/NAc samples. g-i Enriched GO terms for species:time interaction in g AMY, h HT, and i VP/NAc. Hierarchical clustering tree shows relationship between GO categories based on shared genes. Branches with length of zero are subsets of one another. Fractions preceding GO terms indicate proportion of “good” genes that have raw p-value< 0.05 compared to total number of genes in the category. P-values are corrected using Benjamini-Hochberg false discovery rate procedure [35]. Red terms enriched in prairie voles, blue terms enriched in meadow voles. Bold text indicates adjusted p < 0.01, plain text indicates adjusted p < 0.05, and italicized text indicates adjusted p < 0.1 for term. n.s. indicates no significant terms in the categoryBack to article page