Silent mutation rate estimates for Pinus vary 50-fold, ranging from angiosperm-like to among the slowest reported for plants. These differences either reflect extraordinary genomic processes or inconsistent fossil calibration, and they have important consequences for population and biogeographical inferences. Here we estimate mutation rates from 4 Pinus species that represent the major lineages using 11 nuclear and 4 chloroplast loci. Calibration was tested at the divergence of Pinus subgenera with the oldest leaf fossil from subg. Strobus (Eocene; 45 MYA) or a recently published subg. Strobus wood fossil (Cretaceous; 85 MYA). These calibrations place the origin of Pinus 190-102 MYA and give absolute silent rate estimates of 0.70-1.31x10(-9) and 0.22-0.42x10(-9).site-1.year-1 for the nuclear and chloroplast genomes, respectively. These rates are approximately 4- to 20-fold slower than angiosperms, but unlike many previous estimates, they are more consistent with the high per-generation deleterious mutation rates observed in pines. Chronograms from nuclear and chloroplast genomes show that the divergence of subgenera accounts for about half of the time since Pinus diverged from Picea, with subsequent radiations occurring more recently. By extending the sampling to encompass the phylogenetic diversity of Pinus, we predict that most extant subsections diverged during the Miocene. Moreover, subsect. Australes, Ponderosae, and Contortae, containing over 50 extant species, radiated within a 5 Myr time span starting as recently as 18 MYA. An Eocene divergence of pine subgenera (using leaf fossils) does not conflict with fossil-based estimates of the Pinus-Picea split, but a Cretaceous divergence using wood fossils accommodates Oligocene fossils that may represent modern subsections. Because homoplasy and polarity of character states have not been tested for fossil pine assignments, the choice of fossil and calibration node represents a significant source of uncertainty. Based on several lines of evidence (including agreement with ages inferred using calibrations outside of Pinus), we conclude that the 85 MYA calibration at the divergence of pine subgenera provides a reasonable lower bound and that further refinements in age and mutation rate estimates will require a synthetic examination of pine fossil history.

Uncertainties in the age and phylogenetic position of Pinaceae fossils present significant obstacles to our understanding of the timing of diversification in the family. We demonstrate that simultaneous phylogenetic analyses of chloroplast DNA (matK and rbcL) and nonmolecular characters that include both extant genera and a limited number of fossil taxa provide useful hypotheses for calibrating molecular trees. Root placements varied for Pinaceae, with Bayesian analyses recovering mutually monophyletic subfamilies Pinoideae and Abietoideae and parsimony analyses recovering Abietoideae as paraphyletic by placing the root between Cedrus and the remaining genera. The inferred phylogenetic positions of fossil taxa Pityostrobus bernissartensis as the sister group to Pinus and Pseudolarix erensis as the sister group to extant Pseudolarix were used to guide divergencetime calibrations; these calibrations yielded an Early Cretaceous and an Early Jurassic age for crown-group Pinaceae, respectively. The older age estimates based on Pseudolarix erensis are supported by weaker evidence from the fossil record but are consistent with recent reports of Early Cretaceous leaf fossils that appear to coincide with extant genera. There remains a great need to characterize the anatomy of extant and fossil species and to code additional nonmolecular characters.