Phylogeny of Salix subgenus Salix s.l. (Salicaceae): delimitation, biogeography, and reticulate evolution

Phylogenetic insights into the Salicaceae: The evolution of willows and beyond

The Salicaceae includes approximately 54 genera and over 1,400 species with a cosmopolitan distribution. Members of the family are well-known for their diverse secondary plant metabolites, and they play crucial roles in tropical and temperate forest ecosystems. Phylogenetic reconstruction of the Salicaceae has been historically challenging due to the limitations of molecular markers and the extensive history of hybridization and polyploidy within the family. Our study employs whole-genome sequencing of 74 species to generate an extensive phylogeny of the Salicaceae. We generated two RAD-Seq enriched whole-genome sequence datasets and extracted two additional gene sets corresponding to the universal Angiosperms353 and Salicaceae-specific targeted-capture arrays. We reconstructed maximum likelihood-based molecular phylogenies using supermatrix and coalescent-based supertree approaches. Our fossil-calibrated phylogeny estimates that the Salicaceae originated around 128 million years ago and unravels the complex taxonomic relationships within the family. Our findings confirm the non-monophyly of the subgenus Salix s.l. and further support the merging of subgenera Chamaetia and Vetrix, both of which exhibit intricate patterns within and among different sections. Overall, our study not only enhances our understanding of the evolution of the Salicaceae, but also provides valuable insights into the complex relationships within the family.

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.

Challenge accepted: Evolutionary lineages versus taxonomic classification of North American shrub willows (Salix)

PREMISE

The huge diversity of Salix subgenus Chamaetia/Vetrix clade in North America and the lack of phylogenetic resolution within this clade has presented a difficult but fascinating challenge for taxonomists to resolve. Here we tested the existing taxonomic classification with molecular tools.

METHODS

In this study, 132 samples representing 46 species from 22 described sections of shrub willows from the United States and Canada were analyzed and combined with 67 samples from Eurasia. The ploidy levels of the samples were determined using flow cytometry and nQuire. Sequences were produced using a RAD sequencing approach and subsequently analyzed with ipyrad, then used for phylogenetic reconstructions (RAxML, SplitsTree), dating analyses (BEAST, SNAPPER), and character evolution analyses of 14 selected morphological traits (Mesquite).

RESULTS

The RAD sequencing approach allowed the production of a well-resolved phylogeny of shrub willows. The resulting tree showed an exclusively North American (NA) clade in sister position to a Eurasian clade, which included some North American endemics. The NA clade began to diversify in the Miocene. Polyploid species appeared in each observed clade. Character evolution analyses revealed that adaptive traits such as habit and adaxial nectaries evolved multiple times independently.

CONCLUSIONS

The diversity in shrub willows was shaped by an evolutionary radiation in North America. Most species were monophyletic, but the existing sectional classification could not be supported by molecular data. Nevertheless, monophyletic lineages share several morphological characters, which might be useful in the revision of the taxonomic classification of shrub willows.

Phylogeny, Genetic Diversity and Population Structure of Fritillaria cirrhosa and Its Relatives Based on Chloroplast Genome Data

Fritillaria cirrhosa and its relatives have been utilized in traditional Chinese medicine for many years and are under priority protection in China. Despite their medicinal and protective value, research on their phylogeny, genetic diversity, and divergence remains limited. Here, we investigate the chloroplast genome variation architecture of 46 samples of F. cirrhosa and its relatives collected from various regions, encompassing the majority of wild populations across diverse geographical areas. The results indicate abundant variations in 46 accessions including 1659 single-nucleotide polymorphisms and 440 indels. Six variable markers (psbJ, ndhD, ycf1, ndhG, trnT-trnL, and rpl32-trnL) were identified. Phylogenetic and network analysis, population structure analysis, and principal component analysis showed that the 46 accessions formed five clades with significant divergence, which were related to their geographical distribution. The regions spanning from the southern Hengduan Mountains to the Qinghai-Tibet Plateau exhibited the highest levels of genetic diversity. F. cirrhosa and its relatives may have suffered a genetic bottleneck and have a relatively low genetic diversity level. Moreover, geographical barriers and discrete patches may have accelerated population divergence. The study offers novel perspectives on the phylogeny, genetic diversity, and population structure of F. cirrhosa and its relatives, information that can inform conservation and utilization strategies in the future.

Detection of Hybrids in Willows (Salix, Salicaceae) Using Genome-Wide DArTseq Markers

The genus Salix, comprising some 400-500 species, is important in various alluvial or wet habitats of the northern hemisphere. It is a promising crop for applications such as biomass production, biofuels, or environmental projects. Clear species delimitation is crucial in ecology, biotechnology, and horticulture. DArTseq markers, a genome-wide technique, were tested for species and hybrid identification. A total of 179 willow samples were analysed, including six species of Salix subgen. Salix and four species of Salix subgen. Vetrix, including those used in biomass crop production, representing important European taxa. Identification of species-specific markers, clustering analyses (principal coordinate analysis, neighbor-joining) and Bayesian methods (Structure) unambiguously identified putative hybrids. In addition to demonstrating the high efficiency of DArT-seq markers in identifying willow hybrids, we also opened-up new questions about hybridisation processes and systematics. We detected unidirectional hybridisation between S. alba and S. fragilis, forming backcross hybrids, and we rejected the hypothesis that S. fragilis does not occur naturally in Europe. Further, the isolated position of Salix triandra within the genus was confirmed.

Phylogenomics reveals patterns of ancient hybridization and differential diversification that contribute to phylogenetic conflict in willows, poplars, and close relatives

Despite the economic, ecological, and scientific importance of the genera Salix L. (willows) and Populus L. (poplars, cottonwoods, and aspens) Salicaceae, we know little about the sources of differences in species diversity between the genera and of the phylogenetic conflict that often confounds estimating phylogenetic trees. Salix subgenera and sections, in particular, have been difficult to classify, with one recent attempt termed a "spectacular failure" due to a speculated radiation of the subgenera Vetrix and Chamaetia. Here, we use targeted sequence capture to understand the evolutionary history of this portion of the Salicaceae plant family. Our phylogenetic hypothesis was based on 787 gene regions and identified extensive phylogenetic conflict among genes. Our analysis supported some previously described subgeneric relationships and confirmed the polyphyly of others. Using an fbranch analysis, we identified several cases of hybridization in deep branches of the phylogeny, which likely contributed to discordance among gene trees. In addition, we identified a rapid increase in diversification rate near the origination of the Vetrix-Chamaetia clade in Salix. This region of the tree coincided with several nodes that lacked strong statistical support, indicating a possible increase in incomplete lineage sorting due to rapid diversification. The extraordinary level of both recent and ancient hybridization in both Salix and Populus have played important roles in the diversification and diversity in these two genera.

Notes on the Taxonomy of Salix vitellina (Salicaceae)

Salix vitellina L., or golden willow, was described by C. Linnaeus in 1753. It was later considered to be affiliated with S. alba, and its taxonomic rank has been changed to variety, subspecies, and form. A recent proposal designated it as a form of S. alba × S. fragilis. The goal of this study was to verify the taxonomic designation of S. vitellina using morphological characteristics including ovule number. A few specimens of S. vitellina from Europe and North America, including the lectotype LINN1158.13, were analyzed. It was recorded that S. vitellina has an ovule index of 6-10, with most valves with four and five ovules and less than 50% of valves with five ovules. These ovule parameters were similar to those of S. alba. The other floral characteristics also indicated that S. vitellina is associated with S. alba. No signs of androgyny or flower aberrations, commonly occurring in willow hybrids, were found in the specimens of S. vitellina. Thus, the analyses did not corroborate the hybrid origin of S. vitellina. The ovule analysis also confirmed that f. chermesina with orange-red stems is also a taxon of S. alba, which differs from f. vitellina by a greater ovule index of 12-16.

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.

De Novo Assembly and Annotation of 11 Diverse Shrub Willow (Salix) Genomes Reveals Novel Gene Organization in Sex-Linked Regions

Poplar and willow species in the Salicaceae are dioecious, yet have been shown to use different sex determination systems located on different chromosomes. Willows in the subgenus Vetrix are interesting for comparative studies of sex determination systems, yet genomic resources for these species are still quite limited. Only a few annotated reference genome assemblies are available, despite many species in use in breeding programs. Here we present de novo assemblies and annotations of 11 shrub willow genomes from six species. Copy number variation of candidate sex determination genes within each genome was characterized and revealed remarkable differences in putative master regulator gene duplication and deletion. We also analyzed copy number and expression of candidate genes involved in floral secondary metabolism, and identified substantial variation across genotypes, which can be used for parental selection in breeding programs. Lastly, we report on a genotype that produces only female descendants and identified gene presence/absence variation in the mitochondrial genome that may be responsible for this unusual inheritance.

The Ovule Number Variation Provides New Insights into Taxa Delimitation in Willows (Salix subgen. Salix; Salicaceae)

Salix babylonica, S. alba and S. fragilis are closely related species characterized by the lanceolate, acuminate and serrulate leaves. The boundaries between them are defined by relatively few diagnostic characters, and their identification is not fully solved. Recent studies have demonstrated that the number of ovules present in the ovaries of the willow flower can assist in the identification of the species. The detailed ovule data, characteristic for flowers of each species, S. babylonica, S. alba and S. fragilis, and variation in the number of ovules per ovary were documented using many representatives of these species from various geographic regions. The data included the minimum and maximum number of ovules per valve and per ovary and the percentages of valves with a specific number of ovules in a catkin. Some intermediate genotypes and clusters with similar ovule indexes were observed. The important character for the identification of S. babylonica was the presence of valves with 1 or 2 ovules in the ovaries; S. fragilis had valves with 3 ovules while S. alba had the greater number (4-12).

Genomic evidence for contrasting patterns of host-associated genetic differentiation across shared host-plant species in leaf- and bud-galling sawflies

Resource specialization and ecological speciation arising through host-associated genetic differentiation (HAD) are frequently invoked as an explanation for the high diversity of plant-feeding insects and other organisms with a parasitic lifestyle. While genetic studies have demonstrated numerous examples of HAD in insect herbivores, the rarity of comparative studies means that we still lack an understanding of how deterministic HAD is, and whether patterns of host shifts can be predicted over evolutionary timescales. We applied genome-wide single nucleotide polymorphism and mitochondrial DNA sequence data obtained through genome resequencing to define species limits and to compare host-plant use in population samples of leaf- and bud-galling sawflies (Hymenoptera: Tenthredinidae: Nematinae) collected from seven shared willow (Salicaceae: Salix) host species. To infer the repeatability of long-term cophylogenetic patterns, we also contrasted the phylogenies of the two galler groups with each other as well as with the phylogeny of their Salix hosts estimated based on RADseq data. We found clear evidence for host specialization and HAD in both of the focal galler groups, but also that leaf gallers are more specialized to single host species compared with most bud gallers. In contrast to bud gallers, leaf gallers also exhibited statistically significant cophylogenetic signal with their Salix hosts. The observed discordant patterns of resource specialization and host shifts in two related galler groups that have radiated in parallel across a shared resource base indicate a lack of evolutionary repeatability in the focal system, and suggest that short- and long-term host use and ecological diversification in plant-feeding insects are dominated by stochasticity and/or lineage-specific effects.

Repeated turnovers keep sex chromosomes young in willows

Background

Salicaceae species have diverse sex determination systems and frequent sex chromosome turnovers. However, compared with poplars, the diversity of sex determination in willows is poorly understood, and little is known about the evolutionary forces driving their turnover. Here, we characterized the sex determination in two Salix species, S. chaenomeloides and S. arbutifolia, which have an XY system on chromosome 7 and 15, respectively.

Results

Based on the assemblies of their sex determination regions, we found that the sex determination mechanism of willows may have underlying similarities with poplars, both involving intact and/or partial homologs of a type A cytokinin response regulator (RR) gene. Comparative analyses suggested that at least two sex turnover events have occurred in Salix, one preserving the ancestral pattern of male heterogamety, and the other changing heterogametic sex from XY to ZW, which could be partly explained by the “deleterious mutation load” and “sexually antagonistic selection” theoretical models. We hypothesize that these repeated turnovers keep sex chromosomes of willow species in a perpetually young state, leading to limited degeneration.

Conclusions

Our findings further improve the evolutionary trajectory of sex chromosomes in Salicaceae species, explore the evolutionary forces driving the repeated turnovers of their sex chromosomes, and provide a valuable reference for the study of sex chromosomes in other species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02769-w.

The chloroplast genome of Salix floderusii and characterization of chloroplast regulatory elements

Salix floderusii is a rare alpine tree species in the Salix genus. Unfortunately, no extensive germplasm identification, molecular phylogeny, and chloroplast genomics of this plant have been conducted. We sequenced the chloroplast (cp) genome of S. floderusii for the first time using second-generation sequencing technology. The cp genome was 155,540 bp long, including a large single-copy region (LSC, 84,401 bp), a small single-copy region (SSC, 16,221 bp), and inverted repeat regions (IR, 54,918 bp). A total of 131 genes were identified, including 86 protein genes, 37 tRNA genes, and 8 rRNA genes. The S. floderusii cp genome contains 1 complement repeat, 24 forward repeats, 17 palindromic repeats, and 7 reverse repeats. Analysis of the IR borders showed that the IRa and IRb regions of S. floderusii and Salix caprea were shorter than those of Salix cinerea, which may affect plastome evolution. Furthermore, four highly variable regions were found, including the rpl22 coding region, psbM/trnD-GUC non-coding region, petA/psbJ non-coding region, and ycf1 coding region. These high variable regions can be used as candidate molecular markers and as a reference for identifying future Salix species. In addition, phylogenetic analysis indicated that the cp genome of S. floderusii is sister to Salix cupularis and belongs to the Subgenus Vetrix. Genes (Sf-trnI, Sf-PpsbA, aadA, Sf-TpsbA, Sf-trnA) obtained via cloning were inserted into the pBluescript II SK (+) to yield the cp expression vectors, which harbored the selectable marker gene aadA. The results of a spectinomycin resistance test indicated that the cp expression vector had been successfully constructed. Moreover, the aadA gene was efficiently expressed under the regulation of predicted regulatory elements. The present study provides a solid foundation for establishing subsequent S. floderusii cp transformation systems and developing strategies for the genetic improvement of S. floderusii.

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.

Highly Diverse Shrub Willows (Salix L.) Share Highly Similar Plastomes

Plastome phylogenomics is used in a broad range of studies where single markers do not bear enough information. Phylogenetic reconstruction in the genus Salix is difficult due to the lack of informative characters and reticulate evolution. Here, we use a genome skimming approach to reconstruct 41 complete plastomes of 32 Eurasian and North American Salix species representing different lineages, different ploidy levels, and separate geographic regions. We combined our plastomes with published data from Genbank to build a comprehensive phylogeny of 61 samples (50 species) using RAxML (Randomized Axelerated Maximum Likelihood). Additionally, haplotype networks for two observed subclades were calculated, and 72 genes were tested to be under selection. The results revealed a highly conserved structure of the observed plastomes. Within the genus, we observed a variation of 1.68%, most of which separated subg. Salix from the subgeneric Chamaetia/Vetrix clade. Our data generally confirm previous plastid phylogenies, however, within Chamaetia/Vetrix phylogenetic results represented neither taxonomical classifications nor geographical regions. Non-coding DNA regions were responsible for most of the observed variation within subclades and 5.6% of the analyzed genes showed signals of diversifying selection. A comparison of nuclear restriction site associated DNA (RAD) sequencing and plastome data on a subset of 10 species showed discrepancies in topology and resolution. We assume that a combination of (i) a very low mutation rate due to efficient mechanisms preventing mutagenesis, (ii) reticulate evolution, including ancient and ongoing hybridization, and (iii) homoplasy has shaped plastome evolution in willows.

Sex-biased genes and metabolites explain morphologically sexual dimorphism and reproductive costs in Salix paraplesia catkins

Dioecious species evolved from species with monomorphic sex systems in order to achieve overall fitness gains by separating male and female functions. As reproductive organs, unisexual flowers have different reproductive roles and exhibit conspicuous sexual dimorphism. To date, little is known about the temporal variations in and molecular mechanisms underlying the morphology and reproductive costs of dioecious flowers. We investigated male and female flowers of Salix paraplesia in three flowering stages before pollination (the early, blooming and late stages) via transcriptional sequencing as well as metabolite content and phenotypic analysis. We found that a large number of sex-biased genes, rather than sex-limited genes, were responsible for sexual dimorphism in S. paraplesia flowers and that the variation in gene expression in male flowers intensified this situation throughout flower development. The temporal dynamics of sex-biased genes derived from changes in reproductive function during the different flowering stages. Sexually differentiated metabolites related to respiration and flavonoid biosynthesis exhibited the same bias directions as the sex-biased genes. These sex-biased genes were involved mainly in signal transduction, photosynthesis, respiration, cell proliferation, phytochrome biosynthesis, and phenol metabolism; therefore, they resulted in more biomass accumulation and higher energy consumption in male catkins. Our results indicated that sex-biased gene expression in S. paraplesia flowers is associated with different reproductive investments in unisexual flowers; male flowers require a greater reproductive investment to meet their higher biomass accumulation and energy consumption needs.

Phytoremediation of aniline by Salix babylonica cuttings: Removal, accumulation, and photosynthetic response

Aniline, a synthetic compound widely used in industrial and pesticide production, is a potential environmental pollutant. The removal of aniline is extremely important to minimize threats to human health and the surrounding environment. The objectives of this study were to investigate the removal efficiency and physiological response of Salix. babylonica cuttings to aniline pollution. Photosynthesis, chlorophyll fluorescence, spectral reflectance and the concentration of aniline in leaves, stems and roots were analysed. The experiment showed that S. babylonica has a strong removal effect on aniline wastewater. Cuttings from S. babylonica stems and roots played an important role in accumulating aniline. However, this increase in aniline concentration was dose dependent and was not always linear. With increasing aniline concentration in S. babylonica was increasingly stressed, with negative impacts on photosynthesis, chlorophyll fluorescence and spectral reflectance index in S. babylonica leaves. These results indicate that non-stomatal limitations are the main reason for the reduction in Pn in S. babylonica leaves due to chlorophyll structure destruction under aniline stress. In addition, aniline concentrations result in an unbalanced distribution of excitation energy between the two light systems, thereby hindering photosynthetic electron transfer and restricting the efficient operation of photosynthesis. Salix babylonica can endure moderate concentrations of aniline and has potential for the phyto-management of aniline-polluted wastewater, although further studies are needed using polluted wastewater.

Respiratory allergies: Salicaceae sensitization (Review)

Seasonal allergic rhinitis (SAR) is one of the most frequent chronic conditions of the modern world. Pollen carried by the wind from pollinated trees is a major source of SAR. Betulaceae, Oleaceae and Platanus are the most important sources of airway sensitization with regard to tree pollen and, therefore, they are included in the official recommendations of skin prick testing by different official societies. Salicaceae pollen is a moderate source of pollen sensitization. Conversely, large areas are covered with poplars and willows around the world. A number of studies from many countries showed that in some particular situations (large and compacted areas covered by Salicaceae, weather conditions, air pollution, urban ornamental vegetation), poplar and willow pollens may become of local importance in producing SAR. The aim of this review was to present a synthesis of information regarding Salicaceae pollen allergy showing that, if various unfavorable aspects are brought together, a minor problem (Salicaceae sensitization) can became a public health problem.

The Evolutionary History, Diversity, and Ecology of Willows (Salix L.) in the European Alps

The genus Salix (willows), with 33 species, represents the most diverse genus of woody plants in the European Alps. Many species dominate subalpine and alpine types of vegetation. Despite a long history of research on willows, the evolutionary and ecological factors for this species richness are poorly known. Here we will review recent progress in research on phylogenetic relationships, evolution, ecology, and speciation in alpine willows. Phylogenomic reconstructions suggest multiple colonization of the Alps, probably from the late Miocene onward, and reject hypotheses of a single radiation. Relatives occur in the Arctic and in temperate Eurasia. Most species are widespread in the European mountain systems or in the European lowlands. Within the Alps, species differ ecologically according to different elevational zones and habitat preferences. Homoploid hybridization is a frequent process in willows and happens mostly after climatic fluctuations and secondary contact. Breakdown of the ecological crossing barriers of species is followed by introgressive hybridization. Polyploidy is an important speciation mechanism, as 40% of species are polyploid, including the four endemic species of the Alps. Phylogenomic data suggest an allopolyploid origin for all taxa analyzed so far. Further studies are needed to specifically analyze biogeographical history, character evolution, and genome evolution of polyploids.

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.

Complete plastome sequence of Xylosma longifolia Clos. (Salicaceae)

Xylosma longifolia is a tree species within Salicaceae and is distributed in Guizhou, Yunnan, Fujian, Guangxi, Guangdong, and Hainan provinces of China as well as in Vietnam, Laos, and India. There are no studies utilizing the complete plastome of Xylosma longifolia in the current literature. Therefore, this report provides a reference for the plastid gene sequence of Xylosma longifolia, and it contributes to the phylogenetic placement and species identification. In this report, we described the complete plastome sequence of Xylosma longifolia. The complete plastome length of Xylosma longifolia is 156,938 bp and has the typical quadripartite structure and gene content of angiosperms, including two inverted repeat (IR) regions of 27,514 bp, a large single-copy (LSC) region of 85,221 bp and a small single-copy (SSC) region of 16,689 bp. The plastome contains 130 genes, including 86 protein coding genes, 36 tRNA genes, eight rRNA genes (5S rRNA, 4.5S rRNA, 16S rRNA, and 23S rRNA). The GC content of the plastome is 36.8%. The complete plastome sequence will be a valuable resource for studies involving the phylogenetic inference of Salicaceae.

Pollen morphology and variability of Polish native species from genus Salix L

The pollen morphology was studied of 24 Salix species native to Poland, which represented two subgenera, 17 sections and five subsections occurring in Poland. The aim of this study was to discover the taxonomical usefulness of the pollen features under analysis, and to investigate the ranges of their interspecific variability. In total, 720 pollen grains were studied. They were analysed with respect to seven quantitative features (length of the polar axis - P, equatorial diameter - E, length of the ectoaperture - Le, exine thickness - Ex, and P/E, Ex/P and Le/P ratios) and the following qualitative ones: pollen outline and exine ornamentation. The most important features were exine ornamentation (muri, lumina and margo) characters. The pollen features should be treated as auxiliary because they allowed to distinguish eight individual Salix species, and five groups of species. Statistical analysis of the studied traits indicated a high variability among the tested species. The most variable biometric features were P, E and Le, while lower variability occurred in P/E, Le/P and d/E.

Genome sequencing and phylogenetic analysis of allotetraploid Salix matsudana Koidz

Polyploidy is a common phenomenon among willow species. In this study, genome sequencing was conducted for Salix matsudana Koidz (also named Chinese willow), an important greening and arbor tree species, and the genome of this species was compared with those of four other tree species in Salicaceae. The total genome sequence of S. matsudana was 655.72 Mb in size, with repeated sequences accounting for 45.97% of the total length. In total, 531.43 Mb of the genome sequence could be mapped onto 38 chromosomes using the published genetic map as a reference. The genome of S. matsudana could be divided into two groups, the A and B genomes, through homology analysis with the genome of Populus trichocarpa, and the A and B genomes contained 23,985 and 25,107 genes, respectively. 4DTv combined transposon analysis predicted that allotetraploidy in S. matsudana appeared ~4 million years ago. The results from this study will help reveal the evolutionary history of S. matsudana and lay a genetic basis for its breeding.

A targeted sequence capture array for phylogenetics and population genomics in the Salicaceae

PREMISE

The family Salicaceae has proved taxonomically challenging, especially in the genus Salix, which is speciose and features frequent hybridization and polyploidy. Past efforts to reconstruct the phylogeny with molecular barcodes have failed to resolve the species relationships of many sections of the genus.

METHODS

We used the wealth of sequence data in the family to design sequence capture probes to target regions of 300-1200 bp of exonic regions of 972 genes.

RESULTS

We recovered sequence data for nearly all of the targeted genes in three species of Populus and three species of Salix. We present a species tree, discuss concordance among gene trees, and present population genomic summary statistics for these loci.

CONCLUSIONS

Our sequence capture array has extremely high capture efficiency within the genera Populus and Salix, resulting in abundant phylogenetic information. Additionally, these loci show promise for population genomic studies.

Genetic variation in intraspecific populations of Rabdophaga rosaria (Diptera: Cecidomyiidae) indicating possible diversification scenarios into sibling species along with host range expansion on willows (Salicaceae: Salix)

We analysed the COI gene of mitochondrial DNA extracted from larvae of Rabdophaga (Diptera: Cecidomyiidae) that induce rosette galls on Salix in the Holarctic Region. Genetic data show that they belong to species groups of Rabdophaga rosaria and Rabdophaga strobiloides. A maximum likelihood tree indicates that R. rosaria and R. strobiloides populations are divided into clades 1 and 2, respectively. Clade 1 is divided into subclades 1 and 2, and the latter is further divided into subclades 2A and 2B. Subclade 1 consists of populations on several Salix species of section Cinerella in subgenus Vetrix in Georgia and the UK. Subclade 2A contains populations on Salix alba (section Salix, subgenus Salix) in The Netherlands and the UK. Subclade 2B consists of populations on section Helix in Poland, Phylicifoliae in Alaska and Salix species in the Eastern Palaearctic Region. The genetic differences between populations of subclades 1, 2A and 2B range from 1.06 to 3.46%. We propose a possible diversification scenario of R. rosaria into sibling species through the expansion of host plant ranges. Clade 2 consists of R. strobiloides populations on two Alaskan Salix species of the sections Hastatae and Sitchenses in subgenus Vetrix.

Phylogenomic Relationships and Evolution of Polyploid Salix Species Revealed by RAD Sequencing Data

Polyploidy is common in the genus Salix. However, little is known about the origin, parentage and genomic composition of polyploid species because of a lack of suitable molecular markers and analysis tools. We established a phylogenomic framework including species of all described sections of Eurasian shrub willows. We analyzed the genomic composition of seven polyploid willow species in comparison to putative diploid parental species to draw conclusions on their origin and the effects of backcrossing and post-origin evolution. We applied recently developed programs like SNAPP, HyDe, and SNiPloid to establish a bioinformatic pipeline for unravelling the complexity of polyploid genomes. RAD sequencing revealed 23,393 loci and 320,010 high quality SNPs for the analysis of relationships of 35 species of Eurasian shrub willows (Salix subg. Chamaetia/Vetrix). Polyploid willow species appear to be predominantly of allopolyploid origin. More ancient allopolyploidization events were observed for two hexaploid and one octoploid species, while our data suggested a more recent allopolyploid origin for the included tetraploids and identified putative parental taxa. SNiPloid analyses disentangled the different genomic signatures resulting from hybrid origin, backcrossing, and secondary post-origin evolution in the polyploid species. Our RAD sequencing data demonstrate that willow genomes are shaped by ancient and recent reticulate evolution, polyploidization, and post-origin divergence of species.

Biochemical Functions of Glutathione S-Transferase Family of Salix babylonica

Glutathione S-transferases (GSTs) are ubiquitous enzymes that are encoded by a large gene family, and they contribute to the detoxification of endogenous or xenobiotic compounds and oxidative stress metabolism in plants. Although the GSTs gene family has been reported in many land plants, our knowledge of the evolution and function of the willow GSTs is still limited. In this study, 22 full-length GST genes were cloned from Salix babylonica and divided into three classes based on the conserved domain analysis, phylogenetic tree and gene structure: tau, phi and DHAR. The tissue-specific expression patterns were substantially different among the tau and phi GSTs. The Salix GST proteins showed functional divergences in the substrate specificities, substrate activities and kinetic characteristics. The site-directed mutagenesis studies revealed that a single amino acid mutation (Ile/Val53→Thr53) resulted in the lowest activity of SbGSTU7 among the Salix GSTs. These results suggest that non-synonymous substitution of an amino acid at the putative glutathione-binding site may play an important role in the divergence of enzymatic functions of Salix GST family.

Intergeneric Relationships within the Family Salicaceae s.l. based on Plastid Phylogenomics

Many Salicaceae s.l. plants are recognized for their important role in the production of products such as wood, oils, and medicines, and as a model organism in life studies. However, the difference in plastid sequence, phylogenetic relationships, and lineage diversification of the family Salicaceae s.l. remain poorly understood. In this study, we compare 24 species representing 18 genera of the family. Simple sequence repeats (SSRs) are considered effective molecular markers for plant species identification and population genetics. Among them, a total of 1798 SSRs were identified, among which mononucleotide repeat was the most common with 1455 accounts representing 80.92% of the total. Most of the SSRs are located in the non-coding region. We also identified five other types of repeats, including 1750 tandems, 434 forward, 407 palindromic, 86 reverse, and 30 complementary repeats. The species in Salicaceae s.l. have a conserved plastid genome. Each plastome presented a typical quadripartite structure and varied in size due to the expansion and contraction of the inverted repeat (IR) boundary, lacking major structural variations, but we identified six divergence hotspot regions. We obtained phylogenetic relationships of 18 genera in Salicaceae s.l. and the 24 species formed a highly supported lineage. Casearia was identified as the basal clade. The divergence time between Salicaceae s.l. and the outgroup was estimated as ~93 Mya; Salix, and Populus diverged around 34 Mya, consistent with the previously reported time. Our research will contribute to a better understanding of the phylogenetic relationships among the members of the Salicaceae s.l.

Analyzing and Characterizing the Chloroplast Genome of Salix wilsonii

Salix wilsonii is an important ornamental willow tree widely distributed in China. In this study, an integrated circular chloroplast genome was reconstructed for S. wilsonii based on the chloroplast reads screened from the whole-genome sequencing data generated with the PacBio RSII platform. The obtained pseudomolecule was 155,750 bp long and had a typical quadripartite structure, comprising a large single copy region (LSC, 84,638 bp) and a small single copy region (SSC, 16,282 bp) separated by two inverted repeat regions (IR, 27,415 bp). The S. wilsonii chloroplast genome encoded 115 unique genes, including four rRNA genes, 30 tRNA genes, 78 protein-coding genes, and three pseudogenes. Repetitive sequence analysis identified 32 tandem repeats, 22 forward repeats, two reverse repeats, and five palindromic repeats. Additionally, a total of 118 perfect microsatellites were detected, with mononucleotide repeats being the most common (89.83%). By comparing the S. wilsonii chloroplast genome with those of other rosid plant species, significant contractions or expansions were identified at the IR-LSC/SSC borders. Phylogenetic analysis of 17 willow species confirmed that S. wilsonii was most closely related to S. chaenomeloides and revealed the monophyly of the genus Salix. The complete S. wilsonii chloroplast genome provides an additional sequence-based resource for studying the evolution of organelle genomes in woody plants.

The complete chloroplast genome sequence of rose-gold pussy willow, Salix gracilistyla Miq. (Salicaceae)

To understand genetic background of Salix gracilistyla Miq., we presented its complete chloroplast genome which is 155,557 bp and has four sub regions: 84,530 bp of large single copy (LSC) and 16,218 bp of small single copy (SSC) regions are separated by 27,405 bp of inverted repeat (IR) regions including 130 genes (84 protein-coding gene, eight rRNAs, and 38 tRNAs). The overall GC content of the chloroplast genome is 36.7% and those in the LSC, SSC, and IR regions are 34.5%, 31.0%, and 41.9%, respectively. Phylogenetic trees show phylogenetic position of S. gracilistyla with low level of inter-species sequence variations.

Species of the Genus Salix L.: Biochemical Screening and Molecular Docking Approach to Potential Acetylcholinesterase Inhibitors

The genus Salix includes about 500 different, mainly woody species with potentially significant medicinal values. The aim of this study was to evaluate the chemical composition and antioxidant activity of little-studied bark and leaves extracts of seven different species of the genus Salix, and to examine the acetylcholinesterase (AChE) inhibitory potential of selected compounds. The extracts were characterized by High Pressure Liquid Chromatography (HPLC). Total phenolics and flavonoids content was determined spectrophotometrically and the antioxidant activity by 2,2-diphenyl-1-picrylhydrazyl (DPPH•) and hydroxyl radical (•OH) scavenging assays. Molecular docking studies were conducted in order to elucidate the interaction and binding affinity between selected compounds of willow bark and leaves against AChE. The major components in bark and leaves of most of the species were rutin (1.26–22.09 mg/g), salicin (1.62–17.33 mg/g), chlorogenic acid (0.74–7.53 mg/g) and epicatechin (0.71–4.83 mg/g). The latter three compounds demonstrated significant inhibitory potential against AChE in docking studies. All extracts exhibited notable antioxidant activity as scavengers of both DPPH• and •OH. The obtained results indicate that willow species other than those in commercial use, and not only bark, but willow leaves as well, could be utilized as sources of valuable phytocompounds with antioxidant and neuroprotective properties.

Comparative genomics and transcriptomics analysis reveals evolution patterns of selection in the Salix phylogeny

BACKGROUND

Willows are widely distributed in the northern hemisphere and have good adaptability to different living environment. The increasing of genome and transcriptome data provides a chance for comparative analysis to study the evolution patterns with the different origin and geographical distributions in the Salix phylogeny.

RESULTS

Transcript sequences of 10 Salicaceae species were downloaded from public databases. All pairwise of orthologues were identified by comparative analysis in these species, from which we constructed a phylogenetic tree and estimated the rate of diverse. Divergence times were estimated in the 10 Salicaceae using comparative transcriptomic analysis. All of the fast-evolving positive selection sequences were identified, and some cold-, drought-, light-, universal-, and heat- resistance genes were discovered.

CONCLUSIONS

The divergence time of subgenus Vetrix and Salix was about 17.6-16.0 Mya during the period of Middle Miocene Climate Transition (21-14 Mya). Subgenus Vetrix diverged to migratory and resident groups when the climate changed to the cool and dry trend by 14 Mya. Cold- and light- stress genes were involved in positive selection among the resident Vetrix, and which would help them to adapt the cooling stage. Universal- stress genes exhibited positive selection among the migratory group and subgenus Salix. These data are useful for comprehending the adaptive evolution and speciation in the Salix lineage.

Climatic Change Can Influence Species Diversity Patterns and Potential Habitats of Salicaceae Plants in China

Salicaceae is a family of temperate woody plants in the Northern Hemisphere that are highly valued, both ecologically and economically. China contains the highest species diversity of these plants. Despite their widespread human use, how the species diversity patterns of Salicaceae plants formed remains mostly unknown, and these may be significantly affected by global climate warming. Using past, present, and future environmental data and 2673 georeferenced specimen records, we first simulated the dynamic changes in suitable habitats and population structures of Salicaceae. Based on this, we next identified those areas at high risk of habitat loss and population declines under different climate change scenarios/years. We also mapped the patterns of species diversity by constructing niche models for 215 Salicaceae species, and assessed the driving factors affecting their current diversity patterns. The niche models showed Salicaceae family underwent extensive population expansion during the Last Inter Glacial period but retreated to lower latitudes during and since the period of the Last Glacial Maximum. Looking ahead, as climate warming intensifies, suitable habitats will shift to higher latitudes and those at lower latitudes will become less abundant. Finally, the western regions of China harbor the greatest endemism and species diversity of Salicaceae, which are significantly influenced by annual precipitation and mean temperature, ultraviolet-B (UV-B) radiation, and the anomaly of precipitation seasonality. From these results, we infer water–energy dynamic equilibrium and historical climate change are both the main factors likely regulating contemporary species diversity and distribution patterns. Nevertheless, this work also suggests that other, possibly interacting, factors (ambient energy, disturbance history, soil condition) influence the large-scale pattern of Salicaceae species diversity in China, making a simple explanation for it unlikely. Because Southwest China likely served as a refuge for Salicaceae species during the Last Glacial Maximum, it is a current hotspot for endemisms. Under predicted climate change, Salicaceae plants may well face higher risks to their persistence in southwest China, so efforts to support their in-situ conservation there are urgently needed.

Sex-biased gene expression in flowers, but not leaves, reveals secondary sexual dimorphism in Populus balsamifera

Because sexual dimorphism in plants is often less morphologically conspicuous than in animals, studies of sex-biased gene expression may provide a quantitative metric to better address their commonality, molecular pathways, consistency across tissues and taxa, and evolution. The presence of sex-biased gene expression in tissues other than the androecium or gynoecium, termed secondary sexual characters, suggests that these traits arose after the initial evolution of dioecy. Patterns of sequence evolution may provide evidence of positive selection that drove sexual specialization. We compared gene expression in male and female flowers and leaves of Populus balsamifera to assess the extent of sex-biased expression, and tested whether sex-biased genes exhibit elevated rates of protein evolution. Sex-biased expression was pervasive in floral tissue, but nearly absent in leaf tissue. Female-biased genes in flowers were associated with photosynthesis, whereas male-biased genes were associated with mitochondrial function. Sex-biased genes did not exhibit elevated rates of protein evolution, contrary to results from other studies in animals and plants. Our results suggest that the ecological and physiological constraints associated with the energetics of flowering, rather than sexual conflict, have probably shaped the differences in male and female gene expression in P. balsamifera.

Molecular phylogenetics and historical biogeography of the tribe Lilieae (Liliaceae): bi-directional dispersal between biodiversity hotspots in Eurasia

Background and Aims

The role played by the Qinghai-Tibet Plateau (QTP) in the organismal diversification and biogeography of plants in the Northern Hemisphere has attracted much attention from evolutionary biologists. Here we use tribe Lilieae (Liliaceae), including primarily temperate and alpine lineages with disjunct distributions in the North Temperate Zone, as a case study to shed light upon these processes.

Methods

Using 191 taxa (five outgroup taxa) comprising more than 60 % of extant Lilieae species across the entire geographical range, we analyse phylogenetic relationships based on three plastid markers (matK, rbcL, rpl16) and nuclear ITS. Divergence time estimation and ancestral range reconstruction were further inferred.

Key Results

The results support a monophyletic Lilieae divided into four clades. Lilium is nested within Fritillaria, which is paraphyletic and partitioned into two clades, New World and Old World, in the chloroplast DNA (cpDNA) analysis. Incongruences between the ITS and cpDNA trees may be explained by divergent ITS paralogues and hybridization. Lilieae originated around 40-49 (28-67) Mya and probably diversified in the QTP region with four major clades that were established during the Oligocene and the Early Miocene. Uplift of the QTP and climatic changes probably drove early diversification of Lilieae in the QTP region. A rapid radiation occurred during the Late Miocene and the Pleistocene, coinciding temporally with recent orogenic process in the QTP region and climatic oscillations. Several lineages dispersed out of the QTP.

Conclusions

Lineage persistence and explosive radiation were important processes for establishing high species diversity of Lilieae in the QTP region. Both long-distance dispersal and migration across Beringia probably contributed to the modern distribution range of Lilieae. Our study shows that biotic interchanges between the QTP region and Irano-Turanian region and the Mediterranean Basin were bi-directional, suggesting the latter was a secondary centre of diversity.

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.

Variability in the composition of phenolic compounds in winter-dormant Salix pyrolifolia in relation to plant part and age

The phenolic phytochemicals of winter-dormant Salix pyrolifolia were determined from the vegetative buds, and the bark and wood of different-aged twigs by HPLC-DAD and UHPLC-QTOF-MS analyses. All the plant parts were composed of salicylate glucosides and the other Salix-specific, simple phenolic glucosides as well as of phenolic acids, flavonoids and the high molecular-weight condensed tannins. The flavonoid composition was most diverse in buds and they also contained a large amount of chlorogenic acid (5-caffeoylquinic acid IUPAC), while salicylate glucosides and simple phenolic glucosides predominated in bark. The wooden interior part of the twigs contained fewer components and the lowest concentrations of compounds. Salicortin was the main compound in winter-dormant S. pyrolifolia (over 10% of bark biomass), but the concentrations of picein, salireposide, isosalipurposide, catechin and condensed tannins were also high. The flavonoid composition was highly naringenin- and quercetin-biassed. The composition of phytochemicals was organ-specific and remained relatively similar between different-aged trees. However, there were compound-specific fluctuations in the concentrations of phytochemicals with the age of the trees and within plant parts. Generally, the one-year-old plants differed from the older trees in their high concentration of condensed tannins in all the plant parts studied and in the highest concentration of isosalipurposide in bark, while the total amounts of salicylate glucosides in plant parts, and of naringenin glucosides in buds, tended to be highest in 20 year-old-trees.

Plastome phylogeny and lineage diversification of Salicaceae with focus on poplars and willows

Phylogenetic relationships and lineage diversification of the family Salicaceae sensu lato (s.l.) remain poorly understood. In this study, we examined phylogenetic relationships between 42 species from six genera based on the complete plastomes. Phylogenetic analyses of 77 protein coding genes of the plastomes produced good resolution of the interrelationships among most sampled species and the recovered clades. Of the sampled genera from the family, Flacourtia was identified as the most basal and the successive clades comprised both Itoa and Poliothyrsis, Idesia, two genera of the Salicaceae sensu stricto (s.s.) (Populus and Salix). Five major subclades were recovered within the Populus clade. These subclades and their interrelationships are largely inconsistent with morphological classifications and molecular phylogeny based on nuclear internal transcribed spacer sequence variations. Two major subclades were identified for the Salix clade. Molecular dating suggested that species diversification of the major subclades in the Populus and Salix clades occurred mainly within the recent Pliocene. In addition, we found that the rpl32 gene was lost and the rps7 gene evolved into a pseudogene multiple times in the sampled genera of the Salicaceae s.l. Compared with previous studies, our results provide a well-resolved phylogeny from the perspective of the plastomes.

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.

Phylogenomic Analysis and Dynamic Evolution of Chloroplast Genomes in Salicaceae

Chloroplast genomes of plants are highly conserved in both gene order and gene content. Analysis of the whole chloroplast genome is known to provide much more informative DNA sites and thus generates high resolution for plant phylogenies. Here, we report the complete chloroplast genomes of three Salix species in family Salicaceae. Phylogeny of Salicaceae inferred from complete chloroplast genomes is generally consistent with previous studies but resolved with higher statistical support. Incongruences of phylogeny, however, are observed in genus Populus, which most likely results from homoplasy. By comparing three Salix chloroplast genomes with the published chloroplast genomes of other Salicaceae species, we demonstrate that the synteny and length of chloroplast genomes in Salicaceae are highly conserved but experienced dynamic evolution among species. We identify seven positively selected chloroplast genes in Salicaceae, which might be related to the adaptive evolution of Salicaceae species. Comparative chloroplast genome analysis within the family also indicates that some chloroplast genes are lost or became pseudogenes, infer that the chloroplast genes horizontally transferred to the nucleus genome. Based on the complete nucleus genome sequences from two Salicaceae species, we remarkably identify that the entire chloroplast genome is indeed transferred and integrated to the nucleus genome in the individual of the reference genome of P. trichocarpa at least once. This observation, along with presence of the large nuclear plastid DNA (NUPTs) and NUPTs-containing multiple chloroplast genes in their original order in the chloroplast genome, favors the DNA-mediated hypothesis of organelle to nucleus DNA transfer. Overall, the phylogenomic analysis using chloroplast complete genomes clearly elucidates the phylogeny of Salicaceae. The identification of positively selected chloroplast genes and dynamic chloroplast-to-nucleus gene transfers in Salicaceae provide resources to better understand the successful adaptation of Salicaceae species.

Conflicting phylogenomic signals reveal a pattern of reticulate evolution in a recent high-Andean diversification (Asteraceae: Astereae: Diplostephium)

High-throughput sequencing is helping biologists to overcome the difficulties of inferring the phylogenies of recently diverged taxa. The present study analyzes the phylogenetic signal of genomic regions with different inheritance patterns using genome skimming and ddRAD-seq in a species-rich Andean genus (Diplostephium) and its allies. We analyzed the complete nuclear ribosomal cistron, the complete chloroplast genome, a partial mitochondrial genome, and a nuclear-ddRAD matrix separately with phylogenetic methods. We applied several approaches to understand the causes of incongruence among datasets, including simulations and the detection of introgression using the D-statistic (ABBA-BABA test). We found significant incongruence among the nuclear, chloroplast, and mitochondrial phylogenies. The strong signal of hybridization found by simulations and the D-statistic among genera and inside the main clades of Diplostephium indicate reticulate evolution as a main cause of phylogenetic incongruence. Our results add evidence for a major role of reticulate evolution in events of rapid diversification. Hybridization and introgression confound chloroplast and mitochondrial phylogenies in relation to the species tree as a result of the uniparental inheritance of these genomic regions. Practical implications regarding the prevalence of hybridization are discussed in relation to the phylogenetic method.

Evolution in situ: hybrid origin and establishment of willows (Salix L.) on alpine glacier forefields

Little attention has been paid to the evolutionary consequences of the colonizing dynamics and succession processes following glacier retreat. Here we studied hybrid populations that have recently formed and established on glacier forefields of the European Alps owing to secondary contact of a lowland colonizer with a subalpine species. We analyzed the composition of two hybrid populations between Salix purpurea and Salix helvetica with nine microsatellite markers by using Bayesian methods (structure and NewHybrids), and simulations. We also studied niche differentiation between the hybrids and the parental species based on indicator values, soil pH and water retention potential measurements. Allelic structure of hybrids confirms the assumed parentage and in situ origin of the crosses on two independent sites within the last decades. Both hybrid populations comprised F₁ and later generation hybrids (F₂ and backcrosses), confirming hybrid fertility. The parental species showed significant differences in niche characteristics for temperature, soil pH, nutrients and moisture. Remarkably, the hybrids exhibited a higher tolerance to cold temperatures, nutrient-poor and acidic soils than either parent. Our results show that willow hybrids originated after glacier retreat and have established persistent populations within a few decades. One factor contributing to hybrid establishment in sympatry with their parents is their ability to occupy more extreme niches than either parental species within a mosaic-like pattern of microhabitats on the forefield. Introgression and/or transgressive segregation may have resulted in novel genotypes that are able to expand the ecological spectrum of either parent.

Phylogeny of Salix subgenus Salix s.l. (Salicaceae): delimitation, biogeography, and reticulate evolution

Background

The taxonomy and systematics of Salix subgenus Salix s.l. is difficult. The reliability and evolutionary implications of two important morphological characters (number of stamens, and morphology of bud scales) used in subgeneric classification within Salix remain untested, and a disjunct Old–New World distribution pattern of a main clade of subgenus Salix s.l., revealed by a previous study, lacks a reasonable explanation. To study these questions, we conducted phylogenetic analyses based on 4,688 bp of sequence data from four plastid (rbcL, trnD–T, matK, and atpB–rbcL) and two nuclear markers (ETS and ITS) covering all subgenera of Salix, and all sections of subgenus Salix s.l.

Results

Subgenus Salix came out as para- or polyphyletic in both nrDNA and plastid trees. The plastid phylogeny successfully resolved relationships among the major clades of Salix, but resolution within subgenus Salix s.l. remained low. Nevertheless, three monophyletic groups were identifiable in subgenus Salix s.l.: the ‘main clade’ of subgenus Salix s.l., with New and Old World species being reciprocally monophyletic; the section Triandroides clade; and the subgenus Pleuradenia clade. While nrDNA regions showed higher resolution within subgenus Salix s.l., they failed to resolve subgeneric relationships. Extensive, statistically significant gene-tree incongruence was detected across nrDNA–plastid as well as nrDNA ETS–ITS phylogenies, suggesting reticulate evolution or hybridization within the group. The results were supported by network analyses. Ancestral-state reconstructions indicated that multiple stamens and free bud scales represent the plesiomorphic states within Salix, and that several significant shifts in stamen number and bud scale morphology have occurred.

Conclusions

Subgenus Salix s.l. is not monophyletic, and the evolutionary history of the subgenus has involved multiple reticulation events that may mainly be due to hybridization. The delimitation of subgenus Salix s.l. should be redefined by excluding section Triandrae and subgenus Pleuradenia from it. The evolutionary lability of bud-scale morphology and stamen number means that these characters are unreliable bases for classification. The disjunct Old–New World distribution of subgenus Salix s.l. appears to be linked to the profound climatic cooling during the Tertiary, which cut off gene exchange between New and Old World lineages.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0311-7) contains supplementary material, which is available to authorized users.