Phylogeographic Analysis of Lodgepole Pine (Pinus contorta) Reveals Limited Subspecies Differentiation and Evidence for Glacial Refugia

Lodgepole pine (Pinus contorta Dougl.) exhibits pronounced morphological variation across its range, historically attributed to allopatric differentiation during the Wisconsin glaciation. However, whether genetic divergence aligns with morphological differentiation—a fundamental prediction of allopatric speciation theory—remains untested. We conducted a comprehensive phylogeographic analysis of chloroplast DNA (trnL intron and trnL/trnF spacer) and mitochondrial DNA (nad1 b/c intron) across 31 populations representing all four recognized subspecies to test hypotheses of refugial isolation and to evaluate the genetic basis of current taxonomic classification. Contrary to predictions of allopatric divergence, both organellar genomes showed striking genetic uniformity (π = 0.000178–0.000186; intersubspecific genetic distances: 1.06 × 10−4 to 3.96 × 10−4) with no phylogenetic structure corresponding to morphological boundaries. Significant negative neutrality test values (Tajima’s D = −2.26, p < 0.02; Fu and Li’s D* = −4.52, p < 0.02) suggest recent demographic expansion rather than equilibrium divergence. A distinctive 5 bp indel in coastal populations provides molecular evidence for a northern Pacific refugium, and its occurrence in interior populations is consistent with post-glacial pollen-mediated gene flow, though this directionality remains inferential pending nuclear genomic confirmation. These findings suggest that morphological divergence reflects rapid adaptive evolution in heterogeneous environments rather than deep phylogenetic divisions. This pattern exemplifies gene flow-selection balance, in which divergent selection maintains local adaptation despite extensive gene flow—supporting an ecotypic rather than a phylogenetic interpretation of intraspecific diversity. The persistence of morphological variation despite genetic homogeneity indicates strong selection on ecologically important traits, likely driven by variation in fire regimes, differential seed predation, and climate gradients. These results have critical implications for understanding adaptive evolution rates in widespread conifers and for developing conservation strategies that emphasize adaptive processes over taxonomic categories.