Coupling of 12 chromosomal inversions maintains a strong barrier to gene flow between ecotypes

Chromosomal rearrangements lead to the coupling of reproductive barriers, but whether and how they contribute to completion of speciation remains unclear. Marine snails of the genus Littorina repeatedly form hybrid zones between taxa segregating for inversion arrangements, providing opportunities to study this question. Here, we analysed two adjacent transects across hybrid zones between large and dwarf ecotypes of Littorina fabalis covering wave exposure gradients in a Swedish island. Applying whole-genome sequences we found 12 putative inversions reaching near differential fixation between the opposite ends of each transect, and being in strong linkage disequilibrium. These inversions cover 20% of the genome and carry 93% of divergent SNPs. Bimodal hybrid zones in both transects indicate that the two ecotypes of Littorina fabalis maintain their genetic and phenotypic integrity following contact due to strong coupling between inversion clines that strengthened the reproductive barrier. The bimodality resulting from the linked inversions extends into collinear parts of the genome, suggesting a genome-wide coupling. Demographic inference suggests that the coupling built up during a period of allopatry, and has been maintained for more than 1K generations after secondary contact. Overall, this study shows that the coupling of multiple chromosomal inversions contributes to strong reproductive isolation. Importantly, two of the inversions overlap with inverted genomic regions associated with ecotype differences in a closely-related species (L. saxatilis), suggesting the same regions, with similar structural variants, repeatedly contribute to ecotype evolution in distinct species.