Summer Rainfall Drives Adaptation with Gene Flow in a Widespread Butterfly

Understanding how environmental variation interacts with gene flow to shape population genomic patterns is a central goal in evolutionary biology. We investigated how geographic and environmental differences impact genetic variation in the clouded sulfur butterfly (Colias philodice eriphyle) by conducting whole-genome resequencing across replicated transects consisting of paired high- and low-elevation sites on both sides of a major mountain range. Despite sampling across steep environmental gradients, we found no evidence of discrete population structure, providing support for high connectivity across the region or recent divergence with limited drift. Nonetheless, significant isolation by distance-strongest in eastern populations-revealed that geographic distance may impose limits on gene flow, and genetic diversity was also elevated in the east. Genotype-environment association analyses identified more than 16,000 loci associated with elevation, precipitation, and solar radiation. Our redundancy analysis identified precipitation as the strongest predictor of adaptive genomic differentiation, and candidate genes included those linked to melanization and thermoregulation (e.g., TH and yellow). These results demonstrate that even in a largely panmictic population, environmental variation can maintain regional-scale signals of local adaptation. Because insects are declining globally and remain underrepresented in genomic monitoring, conducting whole-genome analyses in a widespread species provides valuable context for assessing how insects today persist across such diverse landscapes and their potential for withstanding future environmental change.