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.

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.

Development, structure and evolutionary significance of seed appendages in Salix matsudana (Salicaceae).

The seeds of Salix and Populus (Salicaceae) are characterized by having numerous long hairs which loosely accompanying the seeds and a small annular appendage which surrounding the base of the seed along with tufted hairs. In this study, the complete development and detailed structure of the hairs and annular appendage in Salix matsudana were investigated using standard techniques for plant anatomy and histochemistry. The results show that the hairs originate successively from the single epidermal cells of the placenta (in megaspore mother cell phase) and funiculus (in eight-nucleate phase), and that their development consists of a progressive increase in cell size and an absence of cell division. The annular appendage is initiated from four to five rows of cells at the distal end of the funiculus in octant proembryo phase and its development is characterized by reactivated meristematic activity and a size increase of these cells. The initiation and development of the hairs are irrelevant to ovule development but fertilization and a developed embryo is necessary for the annular appendage to occur. Considering the reliable fossils, we inferred that the feature of seeds surrounded by long hairs is an ancestral character, and that the detachment of hairs from the funiculus and the occurrence of an annular appendage with tufts of hairs may be the more derived states for seed dispersal in Salix and Populus.

Slow evolution of sex-biased genes in the reproductive tissue of the dioecious plant Salix viminalis.

The relative rate of evolution for sex-biased genes has often been used as a measure of the strength of sex-specific selection. In contrast to studies in a wide variety of animals, far less is known about the molecular evolution of sex-biased genes in plants, particularly in dioecious angiosperms. Here, we investigate the gene expression patterns and evolution of sex-biased genes in the dioecious plant Salix viminalis. We observe lower rates of sequence evolution for male-biased genes expressed in the reproductive tissue compared to unbiased and female-biased genes. These results could be partially explained by the lower codon usage bias for male-biased genes leading to elevated rates of synonymous substitutions compared to unbiased genes. However, the stronger haploid selection in the reproductive tissue of plants, together with pollen competition, would also lead to higher levels of purifying selection acting to remove deleterious variation. Future work should focus on the differential evolution of haploid- and diploid-specific genes to understand the selective dynamics acting on these loci.

Both Volatiles and Cuticular Plant Compounds Determine Oviposition of the Willow Sawfly Nematus oligospilus on Leaves of Salix spp. (Salicaceae).

Plant volatile organic compounds play a role in selection of host plants by herbivorous insects. Once the insect reaches the plant, contact cues determine host acceptance. Although the willow sawfly Nematus oligospilus (Hymenoptera: Tenthredinidae) can differentiate among willow genotypes, no knowledge is available on the cues used by this insect to seek and accept the host plant. In this study, we recorded behavioral orientation in a Y-tube olfactometer of willow sawfly females to volatiles of the highly preferred genotype Salix nigra and the non-preferred genotype S. viminalis. The volatiles released by undamaged willows of each genotype were analyzed by coupled gas chromatography-mass spectrometry. Contact cues were evaluated first by oviposition preference bioassays after selective leaf wax removal, and then by studying the micromorphology of abaxial and adaxial leaf surfaces and their chemical composition. Willow sawfly females oriented preferentially to S. nigra volatiles, which contained more than 3 times the amount of volatiles than that collected from S. viminalis. Analysis of volatiles showed significant differences in amounts of (Z) and (E)-ß-ocimene, undecane, decanal, and ß-caryophyllene. The adaxial leaf surface of S. nigra was less preferred after wax removal, suggesting a role of cuticular waxes for oviposition acceptance. No differences were found among the micromorphology of leaf surfaces between preferred and non-preferred genotypes. The chemical analysis of cuticular waxes showed that the abaxial leaf surface of S. viminalis, which is completely avoided for oviposition, possessed 97% of alkanes. The accepted leaf surfaces contained a more diverse wax profile including alcohols, acids, and esters. Thus, non-alkane wax compounds might be related to oviposition. In sum, our study suggests that several cues act in concert to provide oviposition cues for the sawfly N. oligospilus: females are attracted to volatiles from a distance, and once alighting on the plant, they seek specific chemical contact cues in order to lay eggs.

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.

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.

X-ray micro-computed tomography in willow reveals tissue patterning of reaction wood and delay in programmed cell death.

BACKGROUND: Variation in the reaction wood (RW) response has been shown to be a principle component driving differences in lignocellulosic sugar yield from the bioenergy crop willow. The phenotypic cause(s) behind these differences in sugar yield, beyond their common elicitor, however, remain unclear. Here we use X-ray micro-computed tomography (µCT) to investigate RW-associated alterations in secondary xylem tissue patterning in three dimensions (3D). RESULTS: Major architectural alterations were successfully quantified in 3D and attributed to RW induction. Whilst the frequency of vessels was reduced in tension wood tissue (TW), the total vessel volume was significantly increased. Interestingly, a delay in programmed-cell-death (PCD) associated with TW was also clearly observed and readily quantified by µCT. CONCLUSIONS: The surprising degree to which the volume of vessels was increased illustrates the substantial xylem tissue remodelling involved in reaction wood formation. The remodelling suggests an important physiological compromise between structural and hydraulic architecture necessary for extensive alteration of biomass and helps to demonstrate the power of improving our perspective of cell and tissue architecture. The precise observation of xylem tissue development and quantification of the extent of delay in PCD provides a valuable and exciting insight into this bioenergy crop trait.

Climate change alters leaf anatomy, but has no effects on volatile emissions from Arctic plants.

Biogenic volatile organic compound (BVOC) emissions are expected to change substantially because of the rapid advancement of climate change in the Arctic. BVOC emission changes can feed back both positively and negatively on climate warming. We investigated the effects of elevated temperature and shading on BVOC emissions from arctic plant species Empetrum hermaphroditum, Cassiope tetragona, Betula nana and Salix arctica. Measurements were performed in situ in long-term field experiments in subarctic and high Arctic using a dynamic enclosure system and collection of BVOCs into adsorbent cartridges analysed by gas chromatography-mass spectrometry. In order to assess whether the treatments had resulted in anatomical adaptations, we additionally examined leaf anatomy using light microscopy and scanning electron microscopy. Against expectations based on the known temperature and light-dependency of BVOC emissions, the emissions were barely affected by the treatments. In contrast, leaf anatomy of the studied plants was significantly altered in response to the treatments, and these responses appear to differ from species found at lower latitudes. We suggest that leaf anatomical acclimation may partially explain the lacking treatment effects on BVOC emissions at plant shoot-level. However, more studies are needed to unravel why BVOC emission responses in arctic plants differ from temperate species.

To each its own: differential response of specialist and generalist herbivores to plant defence in willows.

Plant-insect food webs tend to be dominated by interactions resulting from diffuse co-evolution between plants and multiple lineages of herbivores rather than by reciprocal co-evolution and co-cladogenesis. Plants therefore require defence strategies effective against a broad range of herbivore species. In one extreme, plants could develop a single universal defence effective against all herbivorous insects, or tailor-made strategies for each herbivore species. The evolution and ecology of plant defence has to be studied with entire insect assemblages, rather than small subsets of pairwise interactions. The present study examines whether specialists and generalists in three coexisting insect lineages, forming the leaf-chewing guild, respond uniformly to plant phylogeny, secondary metabolites, nutrient content and mechanical antiherbivore defences of their hosts, thus permitting universal plant defence strategies against specialized and generalist folivorous insects from various taxa. The extensive data on folivorous assemblages comprising three insect orders and 193 species are linked with plant phylogeny, secondary chemistry (salicylates, flavonoids and tannins), leaf morphological traits [specific leaf area (SLA) and trichome coverage], nutrient (C : N) content and growth form of eight willow (Salix) and one aspen (Populus) species growing in sympatry. Generalists responded to overall host plant chemistry and trichomes, whilst specialists responded to host plant phylogeny and secondary metabolites that are unique to willows and that are capable of being utilized as an antipredator protection. We did not find any significant impact of other plant traits, that is SLA, C : N ratio, flavonoids, tannins and growth form, on the composition of leaf-chewing communities. Our results show that the response to plant traits is differential among specialists and generalists. This finding constrains the ability of plants to develop defensive traits universally effective against herbivores and may lead to diversification of plant defensive mechanisms into several complementary syndromes, required for effective protection against generalists and specialists from multiple insect taxa comprising most leaf-chewing assemblages. These results point to the necessity of broad studies of plant-herbivore interactions, across multiple insect taxa and guilds.

Herbivores influence the growth, reproduction, and morphology of a widespread Arctic willow.

Shrubs have expanded in Arctic ecosystems over the past century, resulting in significant changes to albedo, ecosystem function, and plant community composition. Willow and rock ptarmigan (Lagopus lagopus, L. muta) and moose (Alces alces) extensively browse Arctic shrubs, and may influence their architecture, growth, and reproduction. Furthermore, these herbivores may alter forage plants in such a way as to increase the quantity and accessibility of their own food source. We estimated the effect of winter browsing by ptarmigan and moose on an abundant, early-successional willow (Salix alaxensis) in northern Alaska by comparing browsed to unbrowsed branches. Ptarmigan browsed 82-89% of willows and removed 30-39% of buds, depending on study area and year. Moose browsed 17-44% of willows and browsed 39-55% of shoots. Browsing inhibited apical dominance and activated axillary and adventitious buds to produce new vegetative shoots. Ptarmigan- and moose-browsed willow branches produced twice the volume of shoot growth but significantly fewer catkins the following summer compared with unbrowsed willow branches. Shoots on browsed willows were larger and produced 40-60% more buds compared to unbrowsed shoots. This process of shoot production at basal parts of the branch is the mechanism by which willows develop a highly complex "broomed" architecture after several years of browsing. Broomed willows were shorter and more likely to be re-browsed by ptarmigan, but not moose. Ptarmigan likely benefit from the greater quantity and accessibility of buds on previously browsed willows and may increase the carrying capacity of their own habitat. Despite the observed tolerance of willows to browsing, their vertical growth and reproduction were strongly inhibited by moose and ptarmigan. Browsing by these herbivores therefore needs to be considered in future models of shrub expansion in the Arctic.

Early selection of novel triploid hybrids of shrub willow with improved biomass yield relative to diploids.

BACKGROUND: Genetic improvement of shrub willow (Salix), a perennial energy crop common to temperate climates, has led to the development of new cultivars with improved biomass yield, pest and disease resistance, and biomass composition suitable for bioenergy applications. These improvements have largely been associated with species hybridization, yet little is known about the genetic mechanisms responsible for improved yield and performance of certain willow species hybrids. RESULTS: The top performing genotypes in this study, representing advanced pedigrees compared with those in previous studies, were mostly triploid in nature and outperformed current commercial cultivars. Of the genotypes studied, the diploids had the lowest mean yield of 8.29 oven dry Mg ha-1 yr-1, while triploids yielded 12.65 Mg ha-1 yr-1, with the top five producing over 16 Mg ha-1 yr-1. Triploids had high stem area and height across all three years of growth in addition to greatest specific gravity. The lowest specific gravity was observed among the tetraploid genotypes. Height was the early trait most correlated with and the best predictor of third-year yield. CONCLUSIONS: These results establish a paradigm for future breeding and improvement of Salix bioenergy crops based on the development of triploid species hybrids. Stem height and total stem area are effective traits for early prediction of relative yield performance.

Floral reward, advertisement and attractiveness to honey bees in dioecious Salix caprea.

In dioecious, zoophilous plants potential pollinators have to be attracted to both sexes and switch between individuals of both sexes for pollination to occur. It often has been suggested that males and females require different numbers of visits for maximum reproductive success because male fertility is more likely limited by access to mates, whereas female fertility is rather limited by resource availability. According to sexual selection theory, males therefore should invest more in pollinator attraction (advertisement, reward) than females. However, our knowledge on the sex specific investment in floral rewards and advertisement, and its effects on pollinator behaviour is limited. Here, we use an approach that includes chemical, spectrophotometric, and behavioural studies i) to elucidate differences in floral nectar reward and advertisement (visual, olfactory cues) in dioecious sallow, Salix caprea, ii) to determine the relative importance of visual and olfactory floral cues in attracting honey bee pollinators, and iii) to test for differential attractiveness of female and male inflorescence cues to honey bees. Nectar amount and sugar concentration are comparable, but sugar composition varies between the sexes. Olfactory sallow cues are more attractive to honey bees than visual cues; however, a combination of both cues elicits the strongest behavioural responses in bees. Male flowers are due to the yellow pollen more colourful and emit a higher amount of scent than females. Honey bees prefer the visual but not the olfactory display of males over those of females. In all, the data of our multifaceted study are consistent with the sexual selection theory and provide novel insights on how the model organism honey bee uses visual and olfactory floral cues for locating host plants.

The effect of limited availability of N or water on C allocation to fine roots and annual fine root turnover in Alnus incana and Salix viminalis.

The effect of limited nitrogen (N) or water availability on fine root growth and turnover was examined in two deciduous species, Alnus incana L. and Salix viminalis L., grown under three different regimes: (i) supply of N and water in amounts which would not hamper growth, (ii) limited N supply and (iii) limited water supply. Plants were grown outdoors during three seasons in covered and buried lysimeters placed in a stand structure and filled with quartz sand. Computer-controlled irrigation and fertilization were supplied through drip tubes. Production and turnover of fine roots were estimated by combining minirhizotron observations and core sampling, or by sequential core sampling. Annual turnover rates of fine roots <1 mm (5-6 year(-1)) and 1-2 mm (0.9-2.8 year(-1)) were not affected by changes in N or water availability. Fine root production (<1 mm) differed between Alnus and Salix, and between treatments in Salix; i.e., absolute length and biomass production increased in the order: water limited < unlimited < N limited. Few treatment effects were detected for fine roots 1-2 mm. Proportionally more C was allocated to fine roots (≤2 mm) in N or water-limited Salix; 2.7 and 2.3 times the allocation to fine roots in the unlimited regime, respectively. Estimated input to soil organic carbon increased by ca. 20% at N limitation in Salix. However, future studies on fine root decomposition under various environmental conditions are required. Fine root growth responses to N or water limitation were less pronounced in Alnus, thus indicating species differences caused by N-fixing capacity and slower initial growth in Alnus, or higher fine root plasticity in Salix. A similar seasonal growth pattern across species and treatments suggested the influence of outer stimuli, such as temperature and light.

Is active management the key to the conservation of saproxylic biodiversity? Pollarding promotes the formation of tree hollows.

Trees with hollows are key features sustaining biodiversity in wooded landscapes. They host rich assemblages of often highly specialised organisms. Hollow trees, however, have become rare and localised in Europe. Many of the associated biota is thus declining or endangered. The challenge of its conservation, therefore, is to safeguard the presence of hollow trees in sufficient numbers. Populations of numerous species associated with tree hollows and dead wood are often found in habitats that were formed by formerly common traditional silvicultural practices such as coppicing, pollarding or pasture. Although it has been occasionally mentioned that such practices increase the formation of hollows and the availability of often sun-exposed dead wood, their effect has never been quantified. Our study examined the hollow incidence in pollard and non-pollard (unmanaged) willows and the effect of pollarding on incremental growth rate by tree ring analysis. The probability of hollow occurrence was substantially higher in pollard than in non-pollard trees. Young pollards, especially, form hollows much more often than non-pollards; for instance, in trees of 50 cm DBH, the probability of hollow ocurrence was ∼0.75 in pollards, but only ∼0.3 in non-pollards. No difference in growth rate was found. Pollarding thus leads to the rapid formation of tree hollows, a habitat usually associated with old trees. It is therefore potentially a very important tool in the restoration of saproxylic habitats and conservation of hollow-dependent fauna. If applied along e.g. roads and watercourses, pollarding could also be used to increase landscape connectivity for saproxylic organisms. In reserves where pollarding was formerly practiced, its restoration would be necessary to prevent loss of saproxylic biodiversity. Our results point to the importance of active management measures for maintaining availability, and spatial and temporal continuity of deadwood microhabitats.

Performance of deep-rooted phreatophytic trees at a site containing total petroleum hydrocarbons.

Poplar and willow tree stands were installed in 2003 at a site in Raleigh, North Carolina containing total petroleum hydrocarbon - contaminated groundwater. The objective was groundwater uptake and plume control. The water table was 5 to 6 m below ground surface (bgs) and therefore methods were used to encourage deep root development. Growth rates, rooting depth and sap flow were measured for trees in Plot A located in the center of the plume and in Plot B peripheral to the plume. The trees were initially sub-irrigated with vertically installed drip-lines and by 2005 had roots 4 to 5 m bgs. Water balance calculations suggested groundwater uptake. In 2007, the average sap flow was higher for Plot B (approximately 59 L per day per tree) than for Plot A (approximately 23 L per day per tree), probably as a result of TPH-induced stress in Plot A. Nevertheless, the estimated rate of groundwater uptake for Plot A was sufficient, relative to the calculated rate of groundwater flux beneath the stand, that a high level of plume control was achieved based on MODFLOW modeling results. Down-gradient groundwater monitoring wells installed in late 2011 should provide quantitative data for plume control.

Recovery performance in xylem hydraulic conductivity is correlated with cavitation resistance for temperate deciduous tree species.

Woody species hydraulically vulnerable to xylem cavitation may experience daily xylem embolism. How such species cope with the possibility of accumulated embolism is unclear. In this study, we examined seven temperate woody species to assess the hypothesis that low cavitation resistance (high vulnerability to cavitation) is compensated by high recovery performance via vessel refilling. We also evaluated leaf functional and xylem structural traits. The xylem recovery index (XRI), defined as the ratio of xylem hydraulic conductivity in plants rewatered after soil drought to that in plants under moist conditions, varied among species. The xylem water potential causing 50% loss of hydraulic conductivity (Ψ50) varied among the species studied, whereas only a slight difference was detected with respect to midday xylem water potential (Ψmin), indicating smaller hydraulic safety margins (Ψmin - Ψ50) for species more vulnerable to cavitation. Cavitation resistance (|Ψ50|) was negatively correlated with XRI across species, with cavitation-vulnerable species showing a higher performance in xylem recovery. Wood density was positively correlated with cavitation resistance and was negatively correlated with XRI. These novel results reveal that coordination exists between cavitation resistance and xylem recovery performance, in association with wood functional traits such as denser wood for cavitation-resistant xylem and less-dense but water-storable wood for refillable xylem. These findings provide insights into long-term maintenance of water transport in tree species growing under variable environmental conditions.

Accumulation and effects of sulfadimethoxine in Salix fragilis L. plants: a preliminary study to phytoremediation purposes.

The application of manure to fertilize arable lands is one of the major means through which veterinary sulfonamides (SAs) enter the environment. Little is known about the capacity of woody plants to phytoremediate this class of antibiotics. To this purpose we performed preliminary studies to evaluate Salix fragilis L. response to sulfadimethoxine (SDM) by investigating both its ability to absorb and tolerate doses of SDM found in fresh faeces of treated calves. Forty cuttings were exposed to either 0, 0.5, 1, or 2 mM of SDM for one month. Decreases in photosynthetic electron transport rate and net CO2 assimilation after 25 days for the higher SDM concentrations were noticed. Moreover, alterations in root morphology of treated plants were observed and further investigated through electron microscopy. However, collected data revealed high root accumulation potential. These preliminary results are promising as they demonstrate that Salix fragilis L. can both absorb and tolerate high concentrations of SAs.

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.