Premium
Putative chromosomal rearrangements are associated primarily with ecotype divergence rather than geographic separation in an intertidal, poorly dispersing snail
Author(s) -
Kess Tony,
Brachmann Matthew,
Boulding Elizabeth G.
Publication year - 2021
Publication title -
journal of evolutionary biology
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.289
H-Index - 128
eISSN - 1420-9101
pISSN - 1010-061X
DOI - 10.1111/jeb.13724
Subject(s) - ecotype , biology , gene flow , reproductive isolation , local adaptation , evolutionary biology , genetic architecture , genetics , genetic divergence , ecological speciation , cline (biology) , quantitative trait locus , gene , genetic variation , population , genetic diversity , demography , sociology
Littorina saxatilis is becoming a model system for understanding the genomic basis of ecological speciation. The parallel formation of crab‐adapted ecotypes that exhibit partial reproductive isolation from wave‐adapted ecotypes has enabled genomic investigation of conspicuous shell traits. Recent genomic studies suggest that chromosomal rearrangements may enable ecotype divergence by reducing gene flow. However, the genomic architecture of traits that are divergent between ecotypes remains poorly understood. Here, we use 11,504 single nucleotide polymorphism (SNP) markers called using the recently released L. saxatilis genome to genotype 462 crab ecotype, wave ecotype and phenotypically intermediate Spanish L. saxatilis individuals with scored phenotypes. We used redundancy analysis to study the genetic architecture of loci associated with shell shape, shape corrected for size, shell size and shell ornamentation, and to compare levels of co‐association among different traits. We discovered 341 SNPs associated with shell traits. Loci associated with trait divergence between ecotypes were often located inside putative chromosomal rearrangements recently characterized in Swedish L. saxatilis . In contrast, we found that shell shape corrected for size varied primarily by geographic site rather than by ecotype and showed little association with these putative rearrangements. We conclude that genomic regions of elevated divergence inside putative rearrangements were associated with divergence of L. saxatilis ecotypes along steep environmental axes—consistent with models of adaptation with gene flow—but were not associated with divergence among the three geographical sites. Our findings support predictions from models indicating the importance of genomic regions of reduced recombination allowing co‐association of loci during ecological speciation with ongoing gene flow.