“…Microsatellites are efficient, high-resolution genetic markers for both population-and species-level studies (Hodel et al, 2016;Li et al, 2017), and have also been used with success when addressing hybridization events and polyploid evolution (e.g., Sampson and Byrne, 2012;Zozomová-Lihová et al, 2014aFeulner et al, 2019;Hu et al, 2019). The Hyb-Seq approach, which captures target exons with flanking intronic and intergenic regions (Weitemier et al, 2014), allows to resolve relationships at various evolutionary levels from the population level up (Villaverde et al, 2018;Tomasello et al, 2020), including polyploids (e.g., Carter et al, 2019). Our specific aims were to (1) resolve the cytotypic and genetic structure of the species under study, (2) identify intraspecific genetic lineages, potentially reflecting descendants from distinct glacial refugia and postglacial colonization routes, (3) explore if and to what extent the lineages correlate with the cytotypic variation and display different environmental niches, (4) resolve single or polytopic polyploid origins, (5) explore the patterns and extent of hybridization between species and lineages, and, finally, (6) to determine how the structure resolved corresponds with previous taxonomic treatments.…”