“…Few other studies have addressed the issue of ecotype formation using model-based inference methods. When such methods are used, they are commonly focused on testing whether divergence occurred with or without gene flow (Klütsch, Manseau, Trim, Polfus, & Wilson, 2016), whether there is ongoing gene flow (Bakovic et al, 2021), as well as on finding the most likely population tree for a given set of sampled populations (Louis et al, 2014) or estimating relevant demographic parameters (Andrew, Kane, Baute, Grassa, & Rieseberg, 2013). Thus, model-based inference methods were rarely used to explicitly contrast different demographic scenarios of ecotype formation, but some examples exist using coalescent-based approaches (Hume, Recknagel, Bean, Adams, & Mable, 2018) coupled with maximum composite-likelihoods (Le Moan, Gagnaire, & Bonhomme, 2016) or ABC (Butlin et al, 2014).…”