Pooled ecotype sequencing reveals candidate genetic mechanisms for adaptive differentiation and reproductive isolation.

The early stages of speciation are often characterized by the formation of partially reproductively isolated ecotypes, which evolve as a by-product of divergent selective forces that are endemic to different habitats. Identifying the genomic regions, genes and ultimately functional polymorphisms that are involved in the processes of ecotype formation is inherently challenging, as there are likely to be many different loci involved in the process. To localize candidate regions of the genome contributing to ecotype formation, we conducted whole-genome pooled sequencing (pool-seq) with 47 coastal perennial and 50 inland annual populations of the yellow monkeyflower, Mimulus guttatus. Coastal perennial and inland annual ecotypes of M. guttatus have previously been shown to be ecologically reproductively isolated and highly locally adapted to their respective habitats. Our pool-seq results found allelic differentiation between the ecotypes for two chromosomal inversions, suggesting that frequencies of inversion heterokaryotypes are strongly differentiated between the ecotypes. Further, there were elevated levels of nonsynonymous change across chromosomal inversions. Across the genome, we identified multiple strong candidate genes potentially driving the morphological, life history and salt tolerance differences between the ecotypes. Several candidate genes coincide with previously identified quantitative trait locus regions and also show a signature of recent natural selection. Overall, the results of our study add to growing support for a major role of chromosomal inversions in adaptation and speciation and provide new insights into the genetic mechanisms underlying classic plant ecotype adaptations to wet and dry habitats.