Gene coexpression networks reveal key drivers of phenotypic divergence in lake whitefish

A functional understanding of processes involved in adaptive divergence is one of the awaiting opportunities afforded by high-throughput transcriptomic technologies. Functional analysis of coexpressed genes has succeeded in the biomedical field in identifying key drivers of disease pathways. However, in ecology and evolutionary biology, functional interpretation of transcriptomic data is still limited. Here, we used Weighted Gene Co-Expression Network Analysis (WGCNA) to identify modules of coexpressed genes in muscle and brain tissue of a lake whitefish backcross progeny. Modules were connected to gradients of known adaptive traits involved in the ecological speciation process between benthic and limnetic ecotypes. Key drivers, that is, hub genes of functional modules related to reproduction, growth, and behavior were identified, and module preservation was assessed in natural populations. Using this approach, we identified modules of coexpressed genes involved in phenotypic divergence and their key drivers, and further identified a module part specifically rewired in the backcross progeny. Functional analysis of transcriptomic data can significantly contribute to the understanding of the mechanisms underlying ecological speciation. Our findings point to bone morphogenetic protein and calcium signaling as common pathways involved in coordinated evolution of trophic behavior, trophic morphology (gill rakers), and reproduction. Results also point to pathways implicating hemoglobins and constitutive stress response (HSP70) governing growth in lake whitefish.