Genetic, morphological and ecological variation across a sharp hybrid zone between two alpine butterfly species

Identifying the mechanisms involved in the formation and maintenance of species is a central question in evolutionary biology, and distinguishing the selective drivers of populations' divergence from demographic processes is of particular interest to better understand the speciation process. Hybrid zones are recognized to provide ideal places to investigate the genetic architecture of speciation and to identify the mechanisms allowing diverging species to maintain their integrity in the face of gene flow. Here, we studied two alpine butterfly species, Coenonympha macromma and C. gardetta, which can be found flying together and hybridizing in narrow contact zones in the southern French Alps. We characterized the genomic composition of individuals, their morphology and their local habitat requirements, within and around a hybrid zone. Genetic diversity analysis at 794 SNPs revealed that all individuals within the hybrid zone were highly admixed, which was not the case outside the hybrid zone. Cline analysis showed that, despite ongoing hybridization, 56 out of 122 loci differentially fixed or nearly so between the two species were impermeable to introgression across the sharp hybrid zone (9 km wide). We also found concordance in cline position and width among genetic, morphological and environmental variation, suggesting a coupling of different reproductive barriers. Habitat characteristics such as the presence of trees and shrubs and the start of the growing season were strongly associated with the genetic variation, and we found evidence of divergence at genetic markers associated with morphology and physiology, putatively involved in visual or environmental reproductive isolation. We discuss the various behavioural and ecological factors that might interplay to maintain current levels of divergence and gene flow between this species pair.