Theory predicts that local adaptation should favor the evolution of a concentrated genetic architecture, where the alleles driving adaptive divergence are tightly clustered on chromosomes. Adaptation to marine versus freshwater environments in threespine stickleback has resulted in an architecture that seems consistent with this prediction: divergence among populations is mainly driven by a few genomic regions harboring multiple quantitative trait loci for environmentally adapted traits, as well as candidate genes with well-established phenotypic effects. One theory for the evolution of these “genomic islands” is that rearrangements remodel the genome to bring causal loci into tight proximity, but this has not been studied explicitly. We tested this theory using synteny analysis to identify micro- and macro-rearrangements in the stickleback genome and assess their potential involvement in the evolution of genomic islands. To identify rearrangements, we conducted a de novo assembly of the closely related tubesnout (Aulorhyncus flavidus) genome and compared this to the genomes of threespine stickleback and two other closely related species. We found that small rearrangements, within-chromosome duplications, and lineage-specific genes (LSGs) were enriched around genomic islands, and that all three chromosomes harboring large genomic islands have experienced macro-rearrangements. We also found that duplicates and micro-rearrangements are 9.9× and 2.9× more likely to involve genes differentially expressed between marine and freshwater genotypes. While not conclusive, these results are consistent with the explanation that strong divergent selection on candidate genes drove the recruitment of rearrangements to yield clusters of locally adaptive loci.

A main goal of evolutionary biology is to understand the genetic basis of adaptive evolution. Although the genes that underlie some adaptive phenotypes are now known, the molecular pathways and regulatory mechanisms mediating the phenotypic effects of those genes often remain a black box. Unveiling this black box is necessary to fully understand the genetic basis of adaptive phenotypes, and to understand why particular genes might be used during phenotypic evolution. Here, we investigated which genes and regulatory mechanisms are mediating the phenotypic effects of the Eda haplotype, a locus responsible for the loss of lateral plates and changes in the sensory lateral line of freshwater threespine stickleback (Gasterosteus aculeatus) populations. Using a combination of RNAseq and a cross design that isolated the Eda haplotype on a fixed genomic background, we found that the Eda haplotype affects both gene expression and alternative splicing of genes related to bone development, neuronal development and immunity. These include genes in conserved pathways, like the BMP, netrin and bradykinin signalling pathways, known to play a role in these biological processes. Furthermore, we found that differentially expressed and differentially spliced genes had different levels of connectivity and expression, suggesting that these factors might influence which regulatory mechanisms are used during phenotypic evolution. Taken together, these results provide a better understanding of the mechanisms mediating the effects of an important adaptive locus in stickleback and suggest that alternative splicing could be an important regulatory mechanism mediating adaptive phenotypes.