Rearranging and completing the puzzle: Phylogenomic analysis of bark beetles Dendroctonus reveals new hypotheses about genus diversification

[Display omitted] •A reference genome improves the precision and accuracy of phylogenomic inference.•Dendroctonus is composed of two major clades, each containing five and four subclades.•Dendroctonus is North American, from this region, some lineages colonized southern regions, Europe and Asia.•Two important colonization events via Beringia occurred during the Oligocene and Miocene•Dendroctonus and pine diversification seem to be associated. Studies carried out on bark beetles within Dendroctonus have been extensive and revealed diverse information in different areas of their natural history, taxonomy, evolution, and interactions, among others. Despite these efforts, phylogenetic hypotheses have remained obscured mainly due to limited information analyzed (taxonomic, gene sampling, or both) in studies focused on obtaining evolutionary hypotheses for this genus. With the aim of filling these gaps in the evolutionary history for Dendroctonus, we analyzed ∼1800 loci mapped to a reference genome obtained for 20 of the 21 species recognized to date, minimizing the impact of missing information and improving the assumption of orthology in a phylogenomic framework. We obtained congruent phylogenetic topologies from two phylogenomic inference strategies: loci concatenation (ML framework) and a multispecies coalescent model (MSC) through the analysis of site pattern frequencies (SNPs). Dendroctonus is composed of two major clades (A and B), each containing five and four subclades, respectively. According to our divergence dating analysis, the MRCA for Dendroctonus dates back to the early Eocene, while the MRCA for each major clade diverged in the mid-Eocene. Interestingly, most of the speciation events of extant species occurred during the Miocene, which could be correlated with the diversification of pine trees (Pinus). The MRCA for Dendroctonus inhabited large regions of North America, with all ancestors and descendants of clade A having diversified within this region. The Mexican Transition Zone is important in the diversification processes for the majority of clade A species. For clade B, we identified two important colonization events to the Old World from America: the first in the early Oligocene from the Arctic to Asia (via Beringia), and the second during the Miocene from the Arctic-Western-Alleghany region to Europe and Siberia (also via Beringia). Our genomic analyses also supported the existence of hidden structured lineages within the frontalis complex