Abstract
Cascading speciation is predicted to occur when multiple interacting species diverge in parallel as a result of divergence in one species promoting adaptive differentiation in other species. However, there are few examples where ecological interactions among taxa have been shown to result in speciation that cascades across multiple trophic levels. Here, we test for cascading speciation occurring among the western pine beetle (Dendroctonus brevicomis), its primary host tree (Pinus ponderosa), and the beetle's fungal mutualists (Ceratocystiopsis brevicomi and Entomocorticium sp. B). We assembled genomes for the beetle and a fungal symbiont and then generated reduced representation genomic data (RADseq) from range-wide samples of these three interacting species. Combined with published data for the host tree, we present clear evidence that the tree, the beetle, and the fungal symbionts are all genetically structured into at least two distinct groups that have strongly codiverged with geographical isolation. We then combine our genomic results with diverse population and laboratory-based data to show evidence for reproductive isolation at each level of the cascade and for coevolution of both antagonistic and mutualistic species interactions within this complex network.