14Despite the general progress in using next generation sequencing techniques for 15 evolutionary research questions, the analysis of polyploid species is still hampered by the lack 16 of suitable analytical tools and the statistical difficulties of dealing with more than two alleles 17 per locus. Polyploidization and especially allopolyploidy leads to new combinations of traits 18 by combining genomes of two or more parental species. This enhances the adaptive potential 19 and often results in speciation. However, multiple origins of polyploids, backcrossing to the 20 parental species and post-origin evolution can strongly influence the genome composition of 21 polyploid species. Here, we used RAD sequencing, which revealed 23,393 loci and 320,010 22 high quality SNPs, to analyze the relationships and origin of seven polyploid species of the 23 diverse genus Salix by utilizing a phylogenomic and a network approach, as well as analyzing 24 the genetic structure and composition of the polyploid genome in comparison to putative 25 parental species. We adapted the SNiPloid pipeline that was originally developed to analyse 26 SNP composition of recently established allotetraploid crop lineages to RAD sequencing data 27 by using concatenated RAD loci as reference. Our results revealed a well-resolved phylogeny 28 of 35 species of Eurasian shrub willows (Salix subg. Chamaetia/Vetrix), including 28 diploid 29 and 7 polyploid species. Polyploidization in willows appears to be predominantly connected 30 to hybridization, i.e. to an allopolyploid origin of species. More ancient allopolyploidization 31 events involving hybridization of more distantly related, ancestral lineages were observed for 32 two hexaploid and one octoploid species. Our data suggested a more recent allopolyploid 33 origin for the included tetraploids within the major subclades and identified putative parental 34 taxa that appear to be plausible in the context of geographical, morphological and ecological 35 patterns. SNiPloid and HyDe analyses disentangled the different genomic signatures resulting 36 from hybrid origin, backcrossing, and secondary post-origin evolution in the polyploid 37 species. All tetraploids showed a considerable post-origin, species-specific proportion of 38 Origin of polyploid Salix species 3 SNPs. The amount of extant hybridization appears to be related to the degree of geographical 39 and ecological isolation of species. Our data demonstrate that high-quality RAD sequencing 40 data are suitable and highly informative for the analysis of the origin and relationships of 41 polyploid species. The combination of the traditional tools RAxML, STRUCTURE, 42 SplitsTree and recently developed programs like SNAPP, HyDe and SNiPloid established a 43 bioinformatic pipeline for unraveling the complexity of polyploid genomes.44