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.

Genome-Wide Characterization of the sHsp Gene Family in Salix suchowensis Reveals Its Functions under Different Abiotic Stresses.

Small heat shock proteins (sHsps) function mainly as molecular chaperones that play vital roles in response to diverse stresses, especially high temperature. However, little is known about the molecular characteristics and evolutionary history of the sHsp family in Salix suchowensis, an important bioenergy woody plant. In this study, 35 non-redundant sHsp genes were identified in S. suchowensis, and they were divided into four subfamilies (C, CP, PX, and MT) based on their phylogenetic relationships and predicted subcellular localization. Though the gene structure and conserved motif were relatively conserved, the sequences of the Hsp20 domain were diversified. Eight paralogous pairs were identified in the Ssu-sHsp family, in which five pairs were generated by tandem duplication events. Ka/Ks analysis indicated that Ssu-sHsps had undergone purifying selection. The expression profiles analysis showed Ssu-Hsps tissue-specific expression patterns, and they were induced by at least one abiotic stress. The expression correlation between two paralogous pairs (Ssu-sHsp22.2-CV/23.0-CV and 23.8-MT/25.6-MT) were less than 0.6, indicating that they were divergent during the evolution. Various cis-acting elements related to stress responses, hormone or development, were detected in the promoter of Ssu-sHsps. Furthermore, the co-expression network revealed the potential mechanism of Ssu-sHsps under stress tolerance and development. These results provide a foundation for further functional research on the Ssu-sHsp gene family in S. suchowensis.

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 F1 and later generation hybrids (F2 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.

Genetic and morphological evidence for introgression between three species of willows.

BACKGROUND: Hybridization and introgression are said to occur relatively frequently in plants, and in particular among different species of willows. However, data on the actual frequency of natural hybridization and introgression is rare. Here, we report the first fine-scale genetic analysis of a contact zone shared between the three basket willow species, Salix dasyclados, S. schwerinii and S. viminalis in the vicinity of the Lake Baikal in Southern Siberia. Individuals were sampled in fourteen populations and classified as pure species or hybrids based on a set of morphological characters. They were then genotyped at 384 nuclear SNP and four chloroplast SSR loci. The STRUCTURE and NewHybrids softwares were used to estimate the frequency and direction of hybridization using genotypic data at the nuclear SNP loci. RESULTS: As many as 19 % of the genotyped individuals were classified as introgressed individuals and these were mainly encountered in the centre of the contact zone. All introgressed individuals were backcrosses to S. viminalis or S. schwerinii and no F1 or F2 hybrids were found. The rest of the genotyped individuals were classified as pure species and formed two clusters, one with S. schwerinii individuals and the other with S. viminalis and S. dasyclados individuals. The two clusters were significantly genetically differentiated, with F ST = 0.333 (0.282-0.382, p < 0.001). In contrast, for the chloroplast haplotypes, no genetic differentiation was observed as they were completely shared between the species. Based on morphological classification only 5 % of the individuals were classified as introgressed individuals, which was much less than what was detected using genotypic data. CONCLUSIONS: We have discovered a new willow hybrid zone with relatively high frequency of introgressed individuals. The low frequency of F1 hybrids indicates that ongoing hybridization is limited, which could be because of the presence of reproductive barriers or simply because the conditions are not favorable for hybridization. We further conclude that in order to get a complete picture of the species composition of a hybrid zone it is necessary to use a combination of morphological characters and genetic data from both nuclear and chloroplast markers.

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.

Consistency between molecular phylogeny and morphological classification of the Salix matsudana Koidz. complex (Salicaceae).

The morphological species concept is based on morpho-logical traits, which are often subject to subjectivity or artifact. Molecular evidence is needed to test the reliability of morphological classification of taxa that are controversial and to provide appropriate taxonomic de-limitation. In this study, we used 15 single-copy nuclear loci and 2 chloroplast fragments to verify the morphological classification of the Salix matsudana Koidz. complex using phylogenetic approaches. Complete sequence alignment showed slight diversification in nuclear sequences and no variety in chloroplast DNA fragments. Phylogenetic trees revealed a monophyletic group consisting of all individuals of S. matsudana and 2 clades within this group, with a 100% bootstrap support value and 1.00 posterior probability. The topology of the phylogenetic trees was highly consistent with the morphological classification of the S. matsudana complex. Verifying the genetic background of these classification units based on remarkable morphological differences will provide a foundation for future studies of Salix and the breeding of new horticultural varieties.

A novel set of single-copy nuclear DNA markers for the genetic study of Salicaceae.

Species of Populus are widely distributed worldwide, playing a significant role in both ecology and economy. However, the lack of single-copy nuclear markers limits knowledge about the phylogeny and population genetics of this genus. In the present study, primer pairs of 15 single-copy nuclear markers were developed through bioinformatic methods based on complete genomic sequences of Populus trichocarpa and Salix arbutifolia. Twenty individuals of Populus davidiana Dode and Salix matsudana Koidz were used to evaluate the basic application of these markers with respect to marker length and diversity indices, respectively. The utility of single-copy nuclear markers is anticipated to facilitate further studies about the phylogeny, population genetics, and phylogeography of this genus, in addition to providing information about the evolutionary dynamics of Salicaceae.

Bioaccessibility in vitro of nutraceuticals from bark of selected Salix species.

The aim of this study was to investigate and to compare the extractability, bioaccessibility, and bioavailability in vitro of antioxidative compounds from bark of selected Salix species: S. alba (SA), S. daphnoides (SD), S. purpurea (SP), and S. daphnoides x purpurea (SDP) hybrid willow clones originating from their natural habitats and cultivated on the sandy soil. The highest amount of phenolic glycosides was found in the bark of SDP and SD. The best source of phenolics was bark of SDP. The highest content of flavonoids were found in SD bark samples, whereas the highest concentration of bioaccessible and bioavailable phenolic acids was determined in SDP bark. Bark of all tested Salix species showed significant antiradical activity. This properties are strongly dependent on extraction system and genetic factors. Regardless of Salix genotypes, the lowest chelating power was found for chemically-extractable compounds. Bark of all Salix species contained ethanol-extractable compounds with reducing ability. Besides this, high bioaccessibility and bioavailability in vitro of Salix bark phytochemicals were found. Obtained results indicate that extracts from bark tested Salix genotypes can provide health promoting benefits to the consumers; however, this problem requires further study.

[Phylogenetic relations of Salix L. subg. Salix species (Salicaceae) according to sequencing data of intergenic spacers of chloroplasic genomes and ITS rDNA].

A phylogenetic analysis based on a comparison of nucleotide sequences of six regions from cpDNA and ITS rDNA (petN-psbM, trnD-trnT, trnC-petN, psaA-ycf3, petG-trnP, and rpoB-trnC) allowed for elucidating the relationship among species and sections belonging to the Salix subgenus and, more generally, to the Salix genus, as well as revealing the relations of the Chosenia genus. The definition of the subgenera of Pleuradenia (including the Urbanianea genus and the Chosenia genus), Salix (without the Triandrae section), Triandrae, and Longifoliae is essentially consistent with current classification schemes of the Salix genus. The previously defined genera of Chosenia and Toisusu (Urbanianea) are not only joined with the Salix genus but are also closely related between themselves. The Protitea subgenus only corresponds to the American species of the Humboldtianae section (S. humboldtiana, S. amygdaloides, S. gooddingii). The relationship of S. chaenomeloides, which is a nomenclatural type of this subgenus, as well as the relationship of the Wilsonia section, remains unresolved. The Humboldtianae section should be interpreted more narrowly, apparently, separating the Acmophillae and Tetraspermae sections from it. The monotypic American Floridanae section is related to the Salix, Salicaster, Tetraspermae, and Wilsonia.

Phylogenetic relationships between four Salix L. species based on DArT markers.

The objectives of this study were to evaluate the usefulness of DArT markers in genotypic identification of willow species and describe genetic relationships between four willow species: Salix viminalis, S. purpurea, S. alba and S. triandra. The experimental plant material comprised 53 willow genotypes of these four species, which are popularly grown in Poland. DArT markers seem to identify Salix species with a high degree of accuracy. As a result, the examined species were divided into four distinct groups which corresponded to the four analyzed species. In our study, we observed that S. triandra was very different genetically from the other species, including S. alba which is generally classified into the same subgenus of Salix. The above corroborates the findings of other authors who relied on molecular methods to reveal that the classification of S. triandra to the subgenus Salix was erroneous. The Principal Coordinate Analysis (PCoA) and the neighbor-joining dendrogram also confirmed the clear division of the studied willow genotypes into four clusters corresponding to individual species. This confirmed the usefulness of DArT markers in taxonomic analyses and identification of willow species.

HPLC of flavanones and chalcones in different species and clones of Salix.

The SPE-HPLC method was developed to determine an isosalipurposide (5) and its derivative, 6"-O-p-coumaroyl ester (6) in the bark of eight taxa (I-VIII) belonging to three species of the genus Salix and originating from a natural habitat or cultivated for pharmaceutical purposes. The chalcones were separated by HPLC under gradient elution with the concentration of ACN increasing from 20% to 50% in 0.1% aqueous H3PO, (tG 15 min). The content of both compounds was determined by an external standardization with the use of isoliquiritigenin (7) as a reference substance - a commercially available chalcone, and also isosalipurposide (5) and its derivative, 6"-O-p-coumaroyl ester (6). The latter compound was isolated from the bark of Salix daphnoides (IV) by a CC and semi-preparative HPLC and its structure was elucidated by MS and NMR spectra. It was stated that 6"-O-p-coumaroylisosalipurposide (6), in addition to isosalipurposide (5), is a characteristic flavonoid for the S. daphnoides species. Moreover, the presence of these two chalcones was confirmed in the bark of S. acutifolia (I). Differences were observed in the results obtained from a quantitative analysis due to the type of reference substance used. The content of chalcones was varied and dependent on the species selected for analysis, namely from 22.01/21.08 mg/g in S. daphnoides clone 1095 (III) to 2.47/2.44 mg/g in S. daphnoides (II), collected from a natural habitat. Isosalipurposide (5) was determined in all the investigated species and clones of Salix, besides a number of naringenin derivatives. Separation of all flavonoids: flavanones - naringenin (1), naringenin (+)-5-O-glucoside (2), (-)-5-O-glucoside (3), 7-O-glucoside (4) and chalcones (5 and 6) was performed under gradient elution with the same solvents and changes in ACN concentration from 2% to 37% (tG 60 min). The total amount of flavanones ranged from 4.69 mg/g in S. purpurea clone 1132 (VII) to 41.93 mg/g in S. purpuea (VIII) from Herbapol Wroclaw.

Leaf epidermal features of Salix species (Salicaceae) and their systematic significance.

PREMISE OF THE STUDY: The classification of the genus Salix has historically been intrinsically difficult due to its propensity toward plasticity and high variation in diagnostic morphological characters. We investigated leaf epidermal characteristics, focusing on the stomatal apparatus because it may provide critical insights into the evolution and taxonomy of Salix and its closely related genera. METHODS: Light microscopy was used to examine the leaf epidermal features in 32 taxa of Salix. KEY RESULTS: Characters such as shape, size, and density of stomatal complexes were very useful in differentiating Salix species. Variation in features of stomatal apparatus in Salix is wider than previously known. Moreover, the type of stomatal complex proved to be very helpful in discriminating Chosenia as members of the genus Salix. CONCLUSIONS: The results of the present study support the placement of Chosenia within Salix and the combining of subgenera Chamaetia and Vetrix because of similarities in their unique stomatal apparatus.

Compound-specific responses of phenolic metabolites in the bark of drought-stressed Salix daphnoides and Salix purpurea.

The bark of willows (Salix spp.) is rich in bioactive phenolic compounds from different compound classes and is therefore used as an herbal remedy. The accumulation of these secondary plant metabolites is influenced by environmental factors, including the availability of water. To analyze the influence of drought stress on the profile of phenolic metabolites in willow bark, a pot experiment with Salix daphnoides Vill. and Salix purpurea L. was conducted. Plants were subjected to three irrigation treatments for four and ten weeks: 65-75% field capacity (well-watered), 33-38% field capacity (moderate drought), and 17-22% field capacity (severe drought). Shoot biomass and proline content were assessed as drought-sensitive traits. Contents of phenolic compounds were analyzed by high-performance liquid chromatography. Drought stress reduced shoot biomass and led to an increase of the bark proline content. The particular effects on phenolics depended on the individual compound, Salix species and drought stress duration. Whereas salicylates were not affected, some flavonoids and phenolic acid derivatives, as well as salireposide indicated treatment effects. The effects comprised decreasing as well as increasing contents. However, beyond the impact of drought stress, the observed responses are assumed to be superimposed by seasonal changes in the content of phenolics. Regarding the yield of willow shoots, the impairment of growth under water shortage seems to be more decisive than drought-induced changes of the bark metabolite content.

Morphological and physiological responses of two willow species from different habitats to salt stress.

Plant salt tolerance is a complex mechanism, and different plant species have different strategies for surviving salt stress. In the present study, we analyzed and compared the morphological and physiological responses of two willow species (Salix linearistipularis and Salix matsudana) from different habitats to salt stress. S. linearistipularis exhibited higher seed germination rates and seedling root Na+ efflux than S. matsudana under salt stress. After salt treatment, S. linearistipularis leaves exhibited less Na+ accumulation, loss of water and chlorophyll, reduction in photosynthetic capacity, and damage to leaf cell structure than leaves of S. matsudana. Scanning electron microscopy combined with gas chromatography mass spectrometry showed that S. linearistipularis leaves had higher cuticular wax loads than S. matsudana leaves. Overall, our results showed that S. linearistipularis had higher salt tolerance than S. matsudana, which was associated with different morphological and physiological responses to salt stress. Furthermore, our study suggested that S. linearistipularis could be a promising tree species for saline-alkali land greening and improvement.

Application of high-performance liquid chromatography for research of salicin in bark of different varieties of Salix.

Willow (Salix L.) species are widely spread in Lithuanian natural dendroflora. Willow bark contains active substances known for anti-inflammatory properties and is known as a phytotherapeutic precursor of aspirin. Bark extracts are components of analgesic and antirheumatic preparations. Therapeutic effectiveness is associated with salicin (2-(hydroxymethyl)phenyl-beta-D-glucopyranoside), which turns into salicylic acid. Increasing attention to natural preparations gives primary importance to research of plants. This study focused on 12 willow taxa and employed routine pharmacopoeia methods. High-performance liquid chromatography method was applied for the analysis of bark extractions. The investigation revealed that not all willow species accumulated a therapeutically sufficient amount of salicin. Bark samples were investigated after 1- and 2-year growth in autumn and spring. Salicin content ranged from 0.08 to 12.6%. Higher contents of active materials were determined in autumn and in 2-year-old willows. Certain willow taxa (Salix alba L., Salix mollissima L., Salix triandra L., Salix viminalis "Americana", Salix dasyclados L.) possessed extremely low salicin amounts. In the second year, analysis covered 32 willow species. Results indicated striking differences in salicin amounts (from 0.04% in Salix viminalis "Americana" to 12.06% in Salix acutifolia). Willow species, plant age, and season should be considered when collecting medicinal plant material. The amount of salicylates in 2-year-old willow bark collected in autumn exceeded by 25% that in 1-year-old willow bark collected in spring. Bark of some analyzed willow species contained the amount of salicylates too low for using as anti-inflammatory or antipyretic remedy.

Genome-wide analysis of Cushion willow provides insights into alpine plant divergence in a biodiversity hotspot.

The Hengduan Mountains (HDM) biodiversity hotspot exhibits exceptional alpine plant diversity. Here, we investigate factors driving intraspecific divergence within a HDM alpine species Salix brachista (Cushion willow), a common component of subnival assemblages. We produce a high-quality genome assembly for this species and characterize its genetic diversity, population structure and pattern of evolution by resequencing individuals collected across its distribution. We detect population divergence that has been shaped by a landscape of isolated sky island-like habitats displaying strong environmental heterogeneity across elevational gradients, combined with population size fluctuations that have occurred since approximately the late Miocene. These factors are likely important drivers of intraspecific divergence within Cushion willow and possibly other alpine plants with a similar distribution. Since intraspecific divergence is often the first step toward speciation, the same factors can be important contributors to the high alpine species diversity in the HDM.

Phytoremediation of chromium using Salix species: cloning ESTs and candidate genes involved in the Cr response.

In this research a differential display based on the detection of cDNA-AFLP markers was used to identify candidate genes potentially involved in the regulation of the response to chromium in four different willow species (Salix alba, Salix eleagnos, Salix fragilis and Salix matsudana) chosen on the basis of their suitability in phytoremediation techniques. Our approach enabled the assay of a large set of mRNA-related fragments and increased the reliability of amplification-based transcriptome analysis. The vast majority of transcript-derived fragments were shared among samples within species and thus attributable to constitutively expressed genes. However, a number of differentially expressed mRNAs were scored in each species and a total of 68 transcripts displaying an altered expression in response to Cr were isolated and sequenced. Public database querying revealed that 44.1% and 4.4% of the cloned ESTs score significant similarity with genes encoding proteins having known or putative function, or with genes coding for unknown proteins, respectively, whereas the remaining 51.5% did not retrieve any homology. Semi-quantitative RT-PCR analysis of seven candidate genes fully confirmed the expression patterns obtained by cDNA-AFLP. Our results indicate the existence of common mechanisms of gene regulation in response to Cr, pathogen attack and senescence-mediated programmed cell death, and suggest a role for the genes isolated in the cross-talk of the signaling pathways governing the adaptation to biotic and abiotic stresses.

Integrating coalescent and ecological niche modeling in comparative phylogeography.

Understanding the factors that contribute to the formation of population genetic structure is a central goal of phylogeographic research, but achieving this goal can be complicated by the stochastic variance inherent to genetic processes. Statistical approaches to testing phylogeographic hypotheses accommodate this stochasticity by evaluating competing models of putative historical population structure, often by simulating null distributions of the expected variance. The effectiveness of these tests depends on the biological realism of the models. Information from the fossil record can aid in reconstructing the historical distributions of some taxa. However, for the majority of taxa, which lack sufficient fossils, paleodistributional modeling can provide valuable spatial-geographic data concerning ancestral distributions. Paleodistributional models are generated by projecting ecological niche models, which predict the current distribution of each species, onto a model of past climatic conditions. Here, we generate paleodistributional models describing the suitable habitat during the last glacial maximum for lineages from the mesic forests of the Pacific Northwest of North America, and use these models to generate alternative phylogeographic hypotheses. Coalescent simulations are then used to test these hypotheses to improve our understanding of the historical events that promoted the formation of population genetic structure in this ecosystem. Results from Pacific Northwest mesic forest organisms demonstrate the utility of these combined approaches. Paleodistribution models and population genetic structure are congruent across three amphibian lineages, suggesting that they have responded in a concerted manner to environmental change. Two other species, a willow and a water vole, despite being currently codistributed and having similar population genetic structure, were predicted by the paleodistributional model to have had markedly different distributions during the last glacial maximum. This suggests that congruent phylogeographic patterns can arise from incongruent ancestral distributions. Paleodistributional models introduce a much-needed spatial-geographic perspective to statistical phylogeography. In conjunction with coalescent models of population genetic structure, they have the potential to improve our understanding of the factors that promote population divergence and ultimately produce regional patterns of biodiversity.

Historical factors shaped species diversity and composition of Salix in eastern Asia.

Ambient energy, niche conservatism, historical climate stability and habitat heterogeneity hypothesis have been proposed to explain the broad-scale species diversity patterns and species compositions, while their relative importance have been controversial. Here, we assessed the relative contributions of contemporary climate, historical climate changes and habitat heterogeneity in shaping Salix species diversity and species composition in whole eastern Asia as well as mountains and lowlands using linear regressions and distance-based redundancy analyses, respectively. Salix diversity was negatively related with mean annual temperature. Habitat heterogeneity was more important than contemporary climate in shaping Salix diversity patterns, and their relative contributions were different in mountains and lowlands. In contrast, the species composition was strongly influenced by contemporary climate and historical climate change than habitat heterogeneity, and their relative contributions were nearly the same both in mountains and lowlands. Our findings supported niche conservatism and habitat heterogeneity hypotheses, but did not support ambient energy and historical climate stability hypotheses. The diversity pattern and species composition of Salix could not be well-explained by any single hypothesis tested, suggesting that other factors such as disturbance history and diversification rate may be also important in shaping the diversity pattern and composition of Salix species.