Phylogeny, diversification rates and species boundaries of Mesoamerican firs (Abies, Pinaceae) in a genus-wide context

[Display omitted] ► Repetitive and non-repetitive cpDNA regions were used to infer the phylogeny of firs. ► The genus Abies is monophyletic and composed by five main phylogeographic groups. ► Only two significantly divergent taxa were observed within the Mesoamerican lineage. ► Divergence within this lineage was prompted by stochastic forces on isolated stands. ► This subtropical clade did not show faster diversification than the temperate ones. The genus Abies is distributed discontinuously in the temperate and subtropical montane forests of the northern hemisphere. In Mesoamerica (Mexico and northern Central America), modern firs originated from the divergence of isolated mountain populations of migrating North American taxa. However, the number of ancestral species, migratory waves and diversification speed of these taxa is unknown. Here, variation in repetitive (Pt30204, Pt63718, and Pt71936) and non-repetitive (rbcL, rps18-rpl20 and trnL-trnF) regions of the chloroplast genome was used to reconstruct the phylogenetic relationships of the Mesoamerican Abies in a genus-wide context. These phylogenies and two fossil-calibrated scenarios were further employed to estimate divergence dates and diversification rates within the genus, and to test the hypothesis that, as in many angiosperms, conifers may exhibit accelerated speciation rates in the subtropics. All phylogenies showed five main clusters that mostly agreed with the currently recognized sections of Abies and with the geographic distribution of species. The Mesoamerican taxa formed a single group with species from southwestern North America of sections Oiamel and Grandis. However, populations of the same species were not monophyletic within this group. Divergence of this whole group dated back to the late Paleocene and the early Miocene depending on the calibration used, which translated in very low diversification rates (r0.0=0.026–0.054, r0.9=0.009–0.019sp/Ma). Such low rates were a constant along the entire genus, including both the subtropical and temperate taxa. An extended phylogeographic analysis on the Mesoamerican clade indicated that Abies flinckii and A. concolor were the most divergent taxa, while the remaining species (A. durangensis, A. guatemalensis, A. hickelii, A. religiosa and A. vejari) formed a single group. Altogether, these results show that divergence of Mesoamerican firs coincides with a model of environmental stasis and decreased extinction rate, being probably prompted by a