Colias eurytheme and C. philodice are sister species with broad sympatry in North America. They hybridize frequently and likely share a significant portion of their genomes through introgression. Both taxa have been ecologically well characterized and exploited to address a broad spectrum of evolutionary issues. Using AFLP markers, we constructed the first linkage map of Colias butterflies. The map is composed of 452 markers spanning 2541.7 cM distributed over 51 linkage groups (40 major groups and 11 small groups with 2-4 markers). Statistical tests indicate that these AFLP markers tend to cluster over the map, with the coefficient of variation of interval sizes being 1.236 (95% C.I. is 1.234-1.240). This nonrandom marker distribution can account for the nonequivalence between the number of linkage groups and the actual haploid chromosome number (N = 31). This study presents the initial step for further marker-assisted research on Colias butterflies, including QTL and introgression analyses. Further investigation of the genomes will help us understand better the roles of introgression and natural selection in the evolution of hybridizing species and devise more appropriate strategies to control these pests.