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.

RAD-seq reveals genetic structure of the F2-generation of natural willow hybrids (Salix L.) and a great potential for interspecific introgression.

BACKGROUND: Hybridization of species with porous genomes can eventually lead to introgression via repeated backcrossing. The potential for introgression between species is reflected by the extent of segregation distortion in later generation hybrids. Here we studied a population of hybrids between Salix purpurea and S. helvetica that has emerged within the last 30 years on a glacier forefield in the European Alps due to secondary contact of the parental species. We used 5758 biallelic SNPs produced by RAD sequencing with the aim to ascertain the predominance of backcrosses (F1 hybrid x parent) or F2 hybrids (F1 hybrid x F1 hybrid) among hybrid offspring. Further, the SNPs were used to study segregation distortion in the second hybrid generation. RESULTS: The analyses in STRUCTURE and NewHybrids revealed that the population consisted of parents and F1 hybrids, whereas hybrid offspring consisted mainly of backcrosses to either parental species, but also some F2 hybrids. Although there was a clear genetic differentiation between S. purpurea and S. helvetica (FST = 0.24), there was no significant segregation distortion in the backcrosses or the F2 hybrids. Plant height of the backcrosses resembled the respective parental species, whereas F2 hybrids were more similar to the subalpine S. helvetica. CONCLUSIONS: The co-occurrence of the parental species and the hybrids on the glacier forefield, the high frequency of backcrossing, and the low resistance to gene flow via backcrossing make a scenario of introgression in this young hybrid population highly likely, potentially leading to the transfer of adaptive traits. We further suggest that this willow hybrid population may serve as a model for the evolutionary processes initiated by recent global warming.

Chromosome-scale assembly of the genome of Salix dunnii reveals a male-heterogametic sex determination system on chromosome 7.

Sex determination systems in plants can involve either female or male heterogamety (ZW or XY, respectively). Here we used Illumina short reads, Oxford Nanopore Technologies (ONT) long reads and Hi-C reads to assemble the first chromosome-scale genome of a female willow tree (Salix dunnii), and to predict genes using transcriptome sequences and available databases. The final genome sequence of 328 Mb in total was assembled in 29 scaffolds, and includes 31,501 predicted genes. Analyses of short-read sequence data that included female and male plants suggested a male heterogametic sex-determining factor on chromosome 7, implying that, unlike the female heterogamety of most species in the genus Salix, male heterogamety evolved in the subgenus Salix. The S. dunnii sex-linked region occupies about 3.21 Mb of chromosome 7 in females (representing its position in the X chromosome), probably within a pericentromeric region. Our data suggest that this region is enriched for transposable element insertions, and about one-third of its 124 protein-coding genes were gained via duplications from other genome regions. We detect purifying selection on the genes that were ancestrally present in the region, though some have been lost. Transcriptome data from female and male individuals show more male- than female-biased genes in catkin and leaf tissues, and indicate enrichment for male-biased genes in the pseudo-autosomal regions. Our study provides valuable genomic resources for further studies of sex-determining regions in the family Salicaceae, and sex chromosome evolution.

RAD sequencing resolved phylogenetic relationships in European shrub willows (Salix L. subg. Chamaetia and subg. Vetrix) and revealed multiple evolution of dwarf shrubs.

The large and diverse genus Salix L. is of particular interest for decades of biological research. However, despite the morphological plasticity, the reconstruction of phylogenetic relationships was so far hampered by the lack of informative molecular markers. Infrageneric classification based on morphology separates dwarf shrubs (subg. Chamaetia) and taller shrubs (subg. Vetrix), while previous phylogenetic studies placed species of these two subgenera just in one largely unresolved clade. Here we want to test the utility of genomic RAD sequencing markers for resolving relationships at different levels of divergence in Salix. Based on a sampling of 15 European species representing 13 sections of the two subgenera, we used five different RAD sequencing datasets generated by ipyrad to conduct phylogenetic analyses. Additionally we reconstructed the evolution of growth form and analyzed the genetic composition of the whole clade. The results showed fully resolved trees in both ML and BI analysis with high statistical support. The two subgenera Chamaetia and Vetrix were recognized as nonmonophyletic, which suggests that they should be merged. Within the Vetrix/Chamaetia clade, a division into three major subclades could be observed. All species were confirmed to be monophyletic. Based on our data, arctic-alpine dwarf shrubs evolved four times independently. The structure analysis showed five mainly uniform genetic clusters which are congruent in sister relationships observed in the phylogenies. Our study confirmed RAD sequencing as a useful genomic tool for the reconstruction of relationships on different taxonomic levels in the genus Salix.