Relationships and genome evolution of polyploid Salix species revealed by RAD sequencing data

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