Evolution of Gene Expression and Splicing in Parallel Cold-Adapted Fly Populations

Changes in gene regulation at multiple levels may comprise an important share of the molecular changes underlying adaptive evolution in nature. However, few studies have assayed within- and between-population variation in gene regulatory traits at a transcriptomic scale, and therefore inferences about the characteristics of adaptive regulatory changes have been elusive. Here, we assess quantitative trait differentiation in gene expression levels and alternative splicing (intron usage) between three closely-related pairs of natural populations of Drosophila melanogaster from contrasting thermal environments that reflect three separate instances of cold tolerance evolution. The cold-adapted populations were known to show population genetic evidence for parallel evolution at the SNP level, and here we find significant although somewhat limited evidence for parallel expression evolution between them, and less evidence for parallel splicing evolution. We find that genes with mitochondrial functions are particularly enriched among candidates for adaptive expression evolution. We also develop a method to estimate cis- versus trans-encoded contributions to expression or splicing differences that does not rely on the presence of fixed differences between parental strains. Applying this method, we infer important roles of both cis- and trans-regulation among our putatively adaptive expression and splicing differences. The apparent contributions of cis- versus trans-regulation to adaptive evolution vary substantially among population pairs, with an Ethiopian pair showing pervasive trans-effects, suggesting that basic characteristics of regulatory evolution may depend on biological context. These findings expand our knowledge of adaptive gene regulatory evolution and our ability to make inferences about this important and widespread process.