Transcriptome profiling of immune tissues reveals habitat‐specific gene expression between lake and river sticklebacks

The observation of habitat‐specific phenotypes suggests the action of natural selection. The three‐spined stickleback (Gasterosteus aculeatus) has repeatedly colonized and adapted to diverse freshwater habitats across the northern hemisphere since the last glaciation, while giving rise to recurring phenotypes associated with specific habitats. Parapatric lake and river populations of sticklebacks harbour distinct parasite communities, a factor proposed to contribute to adaptive differentiation between these ecotypes. However, little is known about the transcriptional response to the distinct parasite pressure of those fish in a natural setting. Here, we sampled wild‐caught sticklebacks across four geographical locations from lake and river habitats differing in their parasite load. We compared gene expression profiles between lake and river populations using 77 whole‐transcriptome libraries from two immune‐relevant tissues, the head kidney and the spleen. Differential expression analyses revealed 139 genes with habitat‐specific expression patterns across the sampled population pairs. Among the 139 differentially expressed genes, eight are annotated with an immune function and 42 have been identified as differentially expressed in previous experimental studies in which fish have been immune challenged. Together, these findings reinforce the hypothesis that parasites contribute to adaptation of sticklebacks in lake and river habitats.