Inferring the past and present connectivity across the range of a North American leaf beetle: combining ecological niche modeling and a geographically explicit model of coalescence.

The leaf beetle Chrysomela aeneicollis occurs across Western North America, either at high elevation or in small, isolated populations along the coast, and thus has a highly fragmented distribution. DNA sequence data (three loci) were collected from five regions across the species range. Population connectivity was examined using traditional ecological niche modeling, which suggested that gene flow could occur among regions now and in the past. We developed geographically explicit coalescence models of sequence evolution that incorporated a two-dimensional representation of the hypothesized ranges suggested by the niche-modeling estimates. We simulated sequence data according to these models and compared them to observed sequences to identify most probable scenarios regarding the migration history of C. aeneicollis. Our results disagreed with initial niche-modeling estimates by clearly rejecting recent connectivity among regions, and were instead most consistent with a long period of range fragmentation, extending well beyond the last glacial maximum. This application of geographically explicit models of coalescence has highlighted some limitations of the use of climatic variables for predicting the present and past range of a species and has explained aspects of the Pleistocene evolutionary history of a cold-adapted organism in Western North America.