Natural genetic variants at the phosphoglucose isomerase, PGI, gene differ in spatial patterning of their polymorphism among species complexes of Colias butterflies in North America. In both lowland and alpine complexes, molecular-functional properties of the polymorphic genotypes can be used to predict genotype-specific adult flight performances and resulting large genotypic differences in adult fitness components. In the lowland species complex, there is striking uniformity of PGI polymorph frequencies at a number of sites across the American West; this fits with earlier findings of strong, similar differences in fitness components over this range. In an alpine complex, Colias meadii shows similar uniformity of PGI frequencies within habitat types, either montane steppe or alpine tundra, over several hundred kilometres in the absence of dispersal. At the same time, large shifts (10-20%) in frequency of the most common alleles occur between steppe and tundra populations, whether these are isolated or, as in some cases, are in contact and exchange many dispersing adults each generation. Data on male mating success of common C. meadii PGI genotypes in steppe and tundra show heterozygote advantage in both habitat types, with shifts in relative homozygote disadvantage between habitats which are consistent with observed frequency differences. Nonadaptive explanations for this situation are rejected, and alternative, thermal-ecology-based adaptive hypotheses are proposed for later experimental test. These findings show that strong local selection may dominate dispersal as an evolutionary agent, whether or not dispersal is present, and that selection may often be the major force promoting 'cohesion' of species over long distances. This case offers new opportunities for integrating studies of molecular structure and function with ecological aspects of natural selection in the wild, both within and among species.