A phylogeographical analysis of the range disjunction for foxtail pine (Pinus balfouriana, Pinaceae): the role of Pleistocene glaciation

Biogeographical patterns within the California Floristic Province have been greatly affected by geological and climatic events. Here, we investigate the phylogeography of foxtail pine (Pinus balfouriana) in an effort to date its range disjunction using molecular data and to further our understanding of phylogeographical patterns for plants within the California Floristic Province. The distribution of foxtail pine is characterized by a 500-km disjunction separating populations located in the Klamath Mountains from those in the southern Sierra Nevada. Previous authors suggested that this disjunction occurred approximately 4000-8000 years ago during the Holocene Xerotherm when western North America became warmer and drier. Those dates, however, are inconsistent with the morphological differences that separate regional populations into formally recognized subspecies. Using the coalescent-based isolation with migration model and DNA sequence data from the chloroplast, mitochondrial, and nuclear genomes, we evaluate several hypotheses addressing the timing of this range disjunction and its effects on subsequent patterns of gene flow. Results from all three genomes are largely consistent with Middle to Early Pleistocene divergence dates. Those dates correspond to the Sherwin glaciation, which was the largest Pleistocene glacial episode in the Sierra Nevada. Gene flow, moreover, was only documented using data from the chloroplast genome, suggesting that low levels of long-distance pollen dispersal (N(e)m < 0.5) have occurred since this divergence event. These results are extended to a discussion of the biogeographical development of subalpine forests in California.