The explosion of next-generation sequencing technologies has allowed researchers to move from studying single genes, to thousands of genes, and thereby to also consider the relationships within gene networks. Like others, we are interested in understanding how developmental and evolutionary forces shape the expression of individual genes, as well as the interactions among genes. To this end, we characterized the effects of genetic background and developmental environment on brain gene coexpression in two parallel, independent evolutionary lineages of Trinidadian guppies ( Poecilia reticulata ). We asked whether connectivity patterns among genes differed based on genetic background and rearing environment, and whether a gene's connectivity predicted its propensity for expression divergence. In pursuing these questions, we confronted the central challenge that standard approaches fail to control the Type I error and/or have low power in the presence of high dimensionality (i.e., large number of genes) and small sample size, as in many gene expression studies. Using our data as a case study, we detail central challenges, discuss sample size guidelines, and provide rigorous statistical approaches for exploring coexpression differences with small sample sizes. Using these approaches, we find evidence that coexpression relationships differ based on both genetic background and rearing environment. We report greater expression divergence in less connected genes and suggest this pattern may arise and be reinforced by selection.