Reproductive Performance of the Alpine Plant Species Ranunculus kuepferi in a Climatic Elevation Gradient: Apomictic Tetraploids Do Not Show a General Fitness Advantage over Sexual Diploids.

Previous studies on the mountain plant Ranunculus kuepferi concluded that apomictic self-compatible tetraploids have experienced a niche shift toward a colder climate during the Holocene, which suggests a fitness advantage over the sexual, self-sterile diploid parents under cold and stressful high-mountain conditions. However, there is still a lack of information on whether reproductive development would be advantageous for tetraploids. Here, we report on microsporogenesis, megagametogenesis, the dynamics of flower and seed development, and the consequences for reproductive success in a common garden experiment along a 1000 m climatic elevation gradient and in natural populations. Flower buds were initiated in the year preceding anthesis and passed winter in a pre-meiotic stage. Flower morphology differed in the known cytotype-specific way in that tetraploid flowers produced about twice as many carpels and fewer petals, stamens, and pollen grains than diploid flowers. Tetraploids developed precociously aposporous embryo sacs and showed a high rate of developmental disturbances. Sexual seed formation prevailed in diploids and pseudogamous apomixis in tetraploids. Along the elevation gradient, stigma pollen load, pollen performance, and seed output decreased. Combinations of reproductive traits, namely, bypass of meiosis irregularities and uniparental reproduction, might have promoted the vast expansion of apomictic R. kuepferi lines across the European Alps.

Allopolyploidization from two dioecious ancestors leads to recurrent evolution of sex chromosomes.

Polyploidization presents an unusual challenge for species with sex chromosomes, as it can lead to complex combinations of sex chromosomes that disrupt reproductive development. This is particularly true for allopolyploidization between species with different sex chromosome systems. Here, we assemble haplotype-resolved chromosome-level genomes of a female allotetraploid weeping willow (Salix babylonica) and a male diploid S. dunnii. We show that weeping willow arose from crosses between a female ancestor from the Salix-clade, which has XY sex chromosomes on chromosome 7, and a male ancestor from the Vetrix-clade, which has ancestral XY sex chromosomes on chromosome 15. We find that weeping willow has one pair of sex chromosomes, ZW on chromosome 15, that derived from the ancestral XY sex chromosomes in the male ancestor of the Vetrix-clade. Moreover, the ancestral 7X chromosomes from the female ancestor of the Salix-clade have reverted to autosomal inheritance. Duplicated intact ARR17-like genes on the four homologous chromosomes 19 likely have contributed to the maintenance of dioecy during polyploidization and sex chromosome turnover. Taken together, our results suggest the rapid evolution and reversion of sex chromosomes following allopolyploidization in weeping willow.

Sex chromosome turnover plays an important role in the maintenance of barriers to post-speciation introgression in willows.

Almost all species in the genus Salix (willow) are dioecious and willows have variable sex-determining systems, the role of this variation in maintaining species barriers is relatively untested. We first analyzed the sex determination systems (SDS) of two species, Salix cardiophylla and Salix interior, whose positions in the Salix phylogeny make them important for understanding a sex chromosome turnover that has been detected in their relatives, and that changed the system from male (XX/XY) to female (ZW/ZZ) heterogamety. We show that both species have male heterogamety, with sex-linked regions (SLRs) on chromosome 15 (termed a 15XY system). The SLRs occupy 21.3% and 22.8% of the entire reference chromosome, respectively. By constructing phylogenetic trees, we determined the phylogenetic positions of all the species with known SDSs. Reconstruction of ancestral SDS character states revealed that the 15XY system is likely the ancestral state in willows. Turnovers of 15XY to 15ZW and 15XY to 7XY likely contributed to early speciation in Salix and gave rise to major groups of the Vetrix and Salix clades. Finally, we tested introgression among species in the phylogenetic trees based on both autosomes and SLRs separately. Frequent introgression was observed among species with 15XY, 15ZW, and 7XY on autosomes, in contrast to the SLR datasets, which showed less introgression, and in particular no gene flow between 15ZW and 7XY species. We argue that, although SDS turnovers in willow speciation may not create complete reproductive barriers, the evolution of SLRs plays important roles in preventing introgression and maintaining species boundaries.

Converging forms: an examination of sub-Arctic, circumarctic, and Central Asian Ranunculus auricomus agg. populations.

Introduction: Phenotypic complexity in species complexes and recently radiated lineages has resulted in a diversity of forms that have historically been classified into separate taxa. Increasingly, with the proliferation of high-throughput sequencing methods, additional layers of complexity have been recognized, such as frequent hybridization and reticulation, which may call into question the previous morphological groupings of closely related organisms. Methods: We investigated Northern European, Asian, and Beringian populations of Ranunculus auricomus agg. with phylogenomic analysis of 736 genes and 27,586 SNPs in order to deduce the interrelatedness and hybrid origin of this phenotypically and taxonomically complicated group from Europe characterized by a history of hybridization, polyploidy, apomixis, and recent radiation. The ploidy levels and the reproductive mode of the Northern European populations were assessed via flow cytometric seed screening. In addition, in order to examine the phenotypic plasticity of the dwarf forms previously described as species and summarized as the Ranunculus monophyllus group, we conducted climate chamber experiments under cold (northern) and warm (temperate) conditions. Results: The Northern European populations are tetra- to hexaploid and propagate primarily through apomixis. The complex is characterized by highly reticulate relationships. Genetic differentiation of the main clusters has occurred between the above-mentioned geographical regions. We find evidence for the hybrid origin of the taxa in these areas with differing genomic contributions from the geographically nearest European sexual progenitor species. Furthermore, polyphyly in the taxa of the R. monophyllus group is supported. Experiments show low lability in the traits associated with the R. monophyllus group. Discussion: We conclude that multiple adaptations of hybrids to colder climates and shorter vegetation periods have shaped the phenotypes of the R. monophyllus group, and we suggest a formal classification as nothotaxa within the R. auricomus group.

Species delimitation 4.0: integrative taxonomy meets artificial intelligence.

Although species are central units for biological research, recent findings in genomics are raising awareness that what we call species can be ill-founded entities due to solely morphology-based, regional species descriptions. This particularly applies to groups characterized by intricate evolutionary processes such as hybridization, polyploidy, or asexuality. Here, challenges of current integrative taxonomy (genetics/genomics + morphology + ecology, etc.) become apparent: different favored species concepts, lack of universal characters/markers, missing appropriate analytical tools for intricate evolutionary processes, and highly subjective ranking and fusion of datasets. Now, integrative taxonomy combined with artificial intelligence under a unified species concept can enable automated feature learning and data integration, and thus reduce subjectivity in species delimitation. This approach will likely accelerate revising and unraveling eukaryotic biodiversity.

Gap-free X and Y chromosome assemblies of Salix arbutifolia reveal an evolutionary change from male to female heterogamety in willows, without a change in the position of the sex-determining locus.

In the Vetrix clade of Salix, a genus of woody flowering plants, sex determination involves chromosome 15, but an XY system has changed to a ZW system. We studied the detailed genetic changes involved. We used genome sequencing, with chromosome conformation capture (Hi-C) and PacBio HiFi reads to assemble chromosome level gap-free X and Y of Salix arbutifolia, and distinguished the haplotypes in the 15X- and 15Y-linked regions, to study the evolutionary history of the sex-linked regions (SLRs). Our sequencing revealed heteromorphism of the X and Y haplotypes of the SLR, with the X-linked region being considerably larger than the corresponding Y region, mainly due to accumulated repetitive sequences and gene duplications. The phylogenies of single-copy orthogroups within the SLRs indicate that S. arbutifolia and Salix purpurea share an ancestral SLR within a repeat-rich region near the chromosome 15 centromere. During the change in heterogamety, the X-linked region changed to a W-linked one, while the Z was derived from the Y.

Apomixis and the paradox of sex in plants.

BACKGROUND: The predominance of sex in eukaryotes, despite the high costs of meiosis and mating, remains an evolutionary enigma. Many theories have been proposed, none of them being conclusive on its own, and they are, in part, not well applicable to land plants. Sexual reproduction is obligate in embryophytes for the great majority of species. SCOPE: This review compares the main forms of sexual and asexual reproduction in ferns and angiosperms, based on the generation cycling of sporophyte and gametophyte (leaving vegetative propagation aside). The benefits of sexual reproduction for maintenance of genomic integrity in comparison to asexuality are discussed in the light of developmental, evolutionary, genetic and phylogenetic studies. CONCLUSIONS: Asexual reproduction represents modifications of the sexual pathway, with various forms of facultative sexuality. For sexual land plants, meiosis provides direct DNA repair mechanisms for oxidative damage in reproductive tissues. The ploidy alternations of meiosis-syngamy cycles and prolonged multicellular stages in the haploid phase in the gametophytes provide a high efficiency of purifying selection against recessive deleterious mutations. Asexual lineages might buffer effects of such mutations via polyploidy and can purge the mutational load via facultative sexuality. The role of organelle-nuclear genome compatibility for maintenance of genome integrity is not well understood. In plants in general, the costs of mating are low because of predominant hermaphroditism. Phylogenetic patterns in the archaeplastid clade suggest that high frequencies of sexuality in land plants are concomitant with a stepwise increase of intrinsic and extrinsic stress factors. Furthermore, expansion of genome size in land plants would increase the potential mutational load. Sexual reproduction appears to be essential for keeping long-term genomic integrity, and only rare combinations of extrinsic and intrinsic factors allow for shifts to asexuality.

Phylogenomics of Southern European Taxa in the Ranunculus auricomus Species Complex: The Apple Doesn't Fall Far from the Tree.

The taxonomic status of many Southern European taxa of the Ranunculus auricomus complex remains uncertain despite this region's proximity to the native ranges of the sexual progenitor species of the complex. We investigated whether additional sexual progenitor species are present in the Mediterranean region. Utilizing target enrichment of 736 single-copy nuclear gene regions and flow cytometry, we analyzed phylogenomic relationships, the ploidy level, and the reproductive mode in representatives of 16 populations in Southern Europe, with additional sequence data from herbarium collections. Additionally, phased sequence assemblies from suspected nothotaxa were mapped to previously described sexual progenitor species in order to determine hybrid ancestry. We found the majority of Mediterranean taxa to be tetraploid, with hybrid populations propagating primarily via apomixis. Phylogenomic analysis revealed that except for the progenitor species, the Mediterranean taxa are often polyphyletic. Most apomictic taxa showed evidence of mixed heritage from progenitor species, with certain progenitor genotypes having mapped more to the populations from adjacent geographical regions. Geographical trends were found in phylogenetic distance, roughly following an east-to-west longitudinal demarcation of the complex, with apomicts extending to the southern margins. Additionally, we observed post-hybridization divergence between the western and eastern populations of nothotaxa in Southern Europe. Our results support a classification of apomictic populations as nothotaxa, as previously suggested for Central Europe.

Contrasting levels of ß-diversity and underlying phylogenetic trends indicate different paths to chemical diversity in highland and lowland willow species.

Diverse specialised metabolites contributed to the success of vascular plants in colonising most terrestrial habitats. Understanding how distinct aspects of chemical diversity arise through heterogeneous environmental pressures can help us understand the effects of abiotic and biotic stress on plant evolution and community assembly. We examined highland and lowland willow species within a phylogenetic framework to test for trends in their chemical α-diversity (richness) and ß-diversity (variation among species sympatric in elevation). We show that differences in chemistry among willows growing at different elevations occur mainly through shifts in chemical ß-diversity and due to convergence or divergence among species sharing their elevation level. We also detect contrasting phylogenetic trends in concentration and α-diversity of metabolites in highland and lowland willow species. The resulting elevational patterns contribute to the chemical diversity of willows and suggest that variable selective pressure across ecological gradients may, more generally, underpin complex changes in plant chemistry.

Insights into the Taxonomically Challenging Hexaploid Alpine Shrub Willows of Salix Sections Phylicifoliae and Nigricantes (Salicaceae).

The complex genomic composition of allopolyploid plants leads to morphologically diverse species. The traditional taxonomical treatment of the medium-sized, hexaploid shrub willows distributed in the Alps is difficult based on their variable morphological characters. In this study, RAD sequencing data, infrared-spectroscopy, and morphometric data are used to analyze the phylogenetic relationships of the hexaploid species of the sections Nigricantes and Phylicifoliae in a phylogenetic framework of 45 Eurasian Salix species. Both sections comprise local endemics as well as widespread species. Based on the molecular data, the described morphological species appeared as monophyletic lineages (except for S. phylicifolia s.str. and S. bicolor, which are intermingled). Both sections Phylicifoliae and Nigricantes are polyphyletic. Infrared-spectroscopy mostly confirmed the differentiation of hexaploid alpine species. The morphometric data confirmed the molecular results and supported the inclusion of S. bicolor into S. phylicifolia s.l., whereas the alpine endemic S. hegetschweileri is distinct and closely related to species of the section Nigricantes. The genomic structure and co-ancestry analyses of the hexaploid species revealed a geographical pattern for widespread S. myrsinifolia, separating the Scandinavian from the alpine populations. The newly described S. kaptarae is tetraploid and is grouped within S. cinerea. Our data reveal that both sections Phylicifoliae and Nigricantes need to be redefined.

Geometric Morphometric Versus Genomic Patterns in a Large Polyploid Plant Species Complex.

Plant species complexes represent a particularly interesting example of taxonomically complex groups (TCGs), linking hybridization, apomixis, and polyploidy with complex morphological patterns. In such TCGs, mosaic-like character combinations and conflicts of morphological data with molecular phylogenies present a major problem for species classification. Here, we used the large polyploid apomictic European Ranunculus auricomus complex to study relationships among five diploid sexual progenitor species and 75 polyploid apomictic derivate taxa, based on geometric morphometrics using 11,690 landmarked objects (basal and stem leaves, receptacles), genomic data (97,312 RAD-Seq loci, 48 phased target enrichment genes, 71 plastid regions) from 220 populations. We showed that (1) observed genomic clusters correspond to morphological groupings based on basal leaves and concatenated traits, and morphological groups were best resolved with RAD-Seq data; (2) described apomictic taxa usually overlap within trait morphospace except for those taxa at the space edges; (3) apomictic phenotypes are highly influenced by parental subgenome composition and to a lesser extent by climatic factors; and (4) allopolyploid apomictic taxa, compared to their sexual progenitor, resemble a mosaic of ecological and morphological intermediate to transgressive biotypes. The joint evaluation of phylogenomic, phenotypic, reproductive, and ecological data supports a revision of purely descriptive, subjective traditional morphological classifications.

Evolutionary origin and establishment of a dioecious diploid-tetraploid complex.

Polyploids recurrently emerge in angiosperms, but most polyploids are likely to go extinct before establishment due to minority cytotype exclusion, which may be specifically a constraint for dioecious plants. Here we test the hypothesis that a stable sex-determination system and spatial/ecological isolation facilitate the establishment of dioecious polyploids. We determined the ploidy levels of 351 individuals from 28 populations of the dioecious species Salix polyclona, and resequenced 190 individuals of S. polyclona and related taxa for genomic diversity analyses. The ploidy survey revealed a frequency 52% of tetraploids in S. polyclona, and genomic k-mer spectra analyses suggested an autopolyploid origin for them. Comparisons of diploid male and female genomes identified a female heterogametic sex-determining factor on chromosome 15, which probably also acts in the dioecious tetraploids. Phylogenetic analyses revealed two diploid clades and a separate clade/grade of tetraploids with a distinct geographic distribution confined to western and central China, where complex mountain systems create higher levels of environmental heterogeneity. Fossil-calibrated phylogenies showed that the polyploids emerged during 7.6-2.3 million years ago, and population demographic histories largely matched the geological and climatic history of the region. Our results suggest that inheritance of the sex-determining system from the diploid progenitor as intrinsic factor and spatial isolation as extrinsic factor may have facilitated the preservation and establishment of polyploid dioecious populations.

Geographical Parthenogenesis in Alpine and Arctic Plants.

The term "Geographical parthenogenesis" describes the phenomenon that asexual organisms usually occupy larger and more northern distribution areas than their sexual relatives, and tend to colonize previously glaciated areas. Several case studies on alpine and arctic plants confirm the geographical pattern, but the causal factors behind the phenomenon are still unclear. Research of the last decade in several plant families has shed light on the question and evaluated some of the classical evolutionary theories. Results confirmed, in general, that the advantages of uniparental reproduction enable apomictic plants to re-colonize faster in larger and more northern distribution areas. Associated factors like polyploidy seem to contribute mainly to the spatial separation of sexual and asexual cytotypes. Ecological studies suggest a better tolerance of apomicts to colder climates and temperate extremes, whereby epigenetic flexibility and phenotypic plasticity play an important role in occupying ecological niches under harsh conditions. Genotypic diversity appears to be of lesser importance for the distributional success of asexual plants. Classical evolutionary theories like a reduced pressure of biotic interactions in colder climates and hence an advantage to asexuals (Red Queen hypothesis) did not gain support from studies on plants. However, it is also still enigmatic why sexual outcrossing remains the predominant mode of reproduction also in alpine floras. Constraints for the origin of apomixis might play a role. Interestingly, some studies suggest an association of sexuality with abiotic stresses. Light stress in high elevations might explain why most alpine plants retain sexual reproduction despite other environmental factors that would favor apomixis. Directions for future research will be given.

Evolution of a hybrid zone of two willow species (Salix L.) in the European Alps analyzed by RAD-seq and morphometrics.

Natural hybridization of plants can result in many outcomes with several evolutionary consequences, such as hybrid speciation and introgression. Natural hybrid zones can arise in mountain systems as a result of fluctuating climate during the exchange of glacial and interglacial periods, where species retract and expand their territories, resulting in secondary contacts. Willows are a large genus of woody plants with an immense capability of interspecific crossing. In this study, the sympatric area of two diploid sister species, S. foetida and S. waldsteiniana in the eastern European Alps, was investigated to study the genomic structure of populations within and outside their contact zone and to analyze congruence of morphological phenotypes with genetic data. Eleven populations of the two species were sampled across the Alps and examined using phylogenetic network and population genetic structure analyses of RAD Seq data and morphometric analyses of leaves. The results showed that a homoploid hybrid zone between the two species was established within their sympatric area. Patterns of genetic admixture in homoploid hybrids indicated introgression with asymmetric backcrossing to not only one of the parental species but also one hybrid population forming a separate lineage. The lack of F1 hybrids indicated a long-term persistence of the hybrid populations. Insignificant isolation by distance suggests that gene flow can act over large geographical scales. Morphometric characteristics of hybrids supported the molecular data and clearly separated populations of the parental species, but showed intermediacy in the hybrid zone populations with a bias toward S. waldsteiniana. The homoploid hybrid zone might have been established via secondary contact hybridization, and its establishment was fostered by the low genetic divergence of parental species and a lack of strong intrinsic crossing barriers. Incomplete ecological separation and the ability of long-distance dispersal of willows could have contributed to the spatial expansion of the hybrid zone.

The male-heterogametic sex determination system on chromosome 15 of Salix triandra and Salix arbutifolia reveals ancestral male heterogamety and subsequent turnover events in the genus Salix.

Dioecious Salix evolved more than 45 million years ago, but have homomorphic sex chromosomes, suggesting that turnover event(s) prevented major differentiation. Sex chromosome turnover events have been inferred in the sister genus Populus. The genus Salix includes two main clades, Salix and Vetrix, with several previously studied Vetrix clade species having female-heterogametic (ZW) or male-heterogametic (XY) sex-determining systems (SDSs) on chromosome 15, while three Salix clade species have XY SDSs on chromosome 7. We here studied two basal taxa of the Vetrix clade, S. arbutifolia and S. triandra using S. purpurea as the reference genome. Analyses of whole genome resequencing data for genome-wide associations (GWAS) with the sexes and genetic differentiation between the sexes (FST values) showed that both species have male heterogamety with a sex-determining locus on chromosome 15, suggesting an early turnover event within the Vetrix clade, perhaps promoted by sexually antagonistic or (and) sex-ratio selection. Changepoint analysis based on FST values identified small sex-linked regions of ~3.33 Mb and ~2.80 Mb in S. arbutifolia and S. triandra, respectively. The SDS of S. arbutifolia was consistent with recent results that used its own genome as reference. Ancestral state reconstruction of SDS suggests that at least two turnover events occurred in Salix.

Evolution of Transcriptomes in Early-Generation Hybrids of the Apomictic Ranunculus auricomus Complex (Ranunculaceae).

Hybridisation in plants may cause a shift from sexual to asexual seed formation (apomixis). Indeed, natural apomictic plants are usually hybrids, but it is still unclear how hybridisation could trigger the shift to apomixis. The genome evolution of older apomictic lineages is influenced by diverse processes such as polyploidy, mutation accumulation, and allelic sequence divergence. To disentangle the effects of hybridisation from these other factors, we analysed the transcriptomes of flowering buds from artificially produced, diploid F2 hybrids of the Ranunculus auricomus complex. The hybrids exhibited unreduced embryo sac formation (apospory) as one important component of apomixis, whereas their parental species were sexual. We revealed 2915 annotated single-copy genes that were mostly under purifying selection according to dN/dS ratios. However, pairwise comparisons revealed, after rigorous filtering, 79 genes under diversifying selection between hybrids and parents, whereby gene annotation assigned ten of them to reproductive processes. Four genes belong to the meiosis-sporogenesis phase (ASY1, APC1, MSP1, and XRI1) and represent, according to literature records, candidate genes for apospory. We conclude that hybridisation could combine novel (or existing) mutations in key developmental genes in certain hybrid lineages, and establish (together with altered gene expression profiles, as observed in other studies) a heritable regulatory mechanism for aposporous development.

Does polyploidy inhibit sex chromosome evolution in angiosperms?

Dioecy is rare in flowering plants (5-6% of species), but is often controlled genetically by sex-linked regions (SLRs). It has so far been unclear whether, polyploidy affects sex chromosome evolution, as it does in animals, though polyploidy is quite common in angiosperms, including in dioecious species. Plants could be different, as, unlike many animal systems, degenerated sex chromosomes, are uncommon in plants. Here we consider sex determination in plants and plant-specific factors, and propose that constraints created at the origin of polyploids limit successful polyploidization of species with SLRs. We consider the most likely case of a polyploid of a dioecious diploid with an established SLR, and discuss the outcome in autopolyploids and allopolyploids. The most stable system possibly has an SLR on just one chromosome, with a strongly dominant genetic factor in the heterogametic sex (e.g., xxxY male in a tetraploid). If recombination occurs with its homolog, this will prevent Y chromosome degeneration. Polyploidy may also allow for reversibility of multiplied Z or X chromosomes into autosomes. Otherwise, low dosage of Y-linked SLRs compared to their multiple homologous x copies may cause loss of reliable sex-determination at higher ploidy levels. We discuss some questions that can be studied using genome sequencing, chromosome level-assemblies, gene expression studies and analysis of loci under selection.

Untying Gordian knots: unraveling reticulate polyploid plant evolution by genomic data using the large Ranunculus auricomus species complex.

Speciation via hybridization and polyploidization is a major evolutionary force in plant evolution but is still poorly understood for neopolyploid groups. Challenges are attributed to high heterozygosity, low genetic divergence, and missing information on progenitors, ploidy, and reproduction. We study the large Eurasian Ranunculus auricomus species complex and use a comprehensive workflow integrating reduced-representation sequencing (RRS) genomic data to unravel reticulate evolution, genome diversity and composition of polyploids. We rely on 97 312 restriction site-associated DNA sequencing (RAD-Seq) loci, 576 targeted nuclear genes (48 phased), and 71 plastid regions derived from 78 polyploid apomictic taxa and four diploid and one tetraploid putative sexual progenitor species. We applied (phylo)genomic structure, network, and single nucleotide polymorphism (SNP)-origin analyses. Results consistently showed only 3-5 supported and geographically structured polyploid genetic groups, each containing extant sexual and one unknown progenitor species. Combined analyses demonstrated predominantly allopolyploid origins, each involving 2-3 different diploid sexual progenitor species. Young allotetraploids were characterized by subgenome dominance and nonhybrid SNPs, suggesting substantial post-origin but little lineage-specific evolution. The biodiversity of neopolyploid complexes can result from multiple hybrid origins involving different progenitors and substantial post-origin evolution (e.g. homoeologous exchanges, hybrid segregation, gene flow). Reduced-representation sequencing genomic data including multi-approach information is efficient to delimit shallow reticulate relationships.

The phylogeny of Salix revealed by whole genome re-sequencing suggests different sex-determination systems in major groups of the genus.

BACKGROUND AND AIMS: The largest genus of Salicaceae sensu lato, Salix, has been shown to consist of two main clades: clade Salix, in which species have XY sex-determination systems (SDSs) on chromosome 7, and clade Vetrix including species with ZW SDSs on chromosome 15. Here, we test the utility of whole genome re-sequencing (WGR) for phylogenomic reconstructions of willows to infer changes between different SDSs. METHODS: We used more than 1 TB of WGR data from 70 Salix taxa to ascertain single nucleotide polymorphisms on the autosomes, the sex-linked regions (SLRs) and the chloroplast genomes, for phylogenetic and species tree analyses. To avoid bias, we chose reference genomes from both groups, Salix dunnii from clade Salix and S. purpurea from clade Vetrix. KEY RESULTS: Two main largely congruent groups were recovered: the paraphyletic Salix grade and the Vetrix clade. The autosome dataset trees resolved four subclades (C1-C4) in Vetrix. C1 and C2 comprise species from the Hengduan Mountains and adjacent areas and from Eurasia, respectively. Section Longifoliae (C3) grouped within the Vetrix clade but fell into the Salix clade in trees based on the chloroplast dataset analysis. Salix triandra from Eurasia (C4) was revealed as sister to the remaining species of clade Vetrix. In Salix, the polyploid group C5 is paraphyletic to clade Vetrix and subclade C6 is consistent with Argus's subgenus Protitea. Chloroplast datasets separated both Vetrix and Salix as monophyletic, and yielded C5 embedded within Salix. Using only diploid species, both the SLR and autosomal datasets yielded trees with Vetrix and Salix as well-supported clades. CONCLUSION: WGR data are useful for phylogenomic analyses of willows. The different SDSs may contribute to the isolation of the two major groups, but the reproductive barrier between them needs to be studied.