Pan-Transcriptome Analysis of Willow Species from Diverse Geographic Distributions

Willows, in the genus Salix, are widespread on the earth with significant ecological and economic values for humans. Although about 500 Salix species have been estimated, the genomic foundation of their adaptations to environments with diverse stresses has been underexplored. Here, we applied a pan-transcriptome approach to investigate the phylogenetic relationships and genetic variations among 16 willow species. A pan-transcriptome of 29,668 gene families was assembled, 69% of which exhibited presence/absence variation across the analyzed species. In comparison to core genes present in all species, shell gene families absent in at least one species were enriched with genes in pathways of signaling transduction and response to stimuli, suggesting their functions in the interaction with diverse environmental factors. A phylogenetic tree of 16 willow species was constructed with high confidence based on 870 single-copy orthologous genes, providing detailed evolutionary relationships of willow sections. The willow species were further assigned into four species clusters using the gene numbers in each family. The diversity of gene family size and gene expression levels among the willow species are closely associated with their geographical distributions. The gene family members involved in DNA repair and cellular response to DNA damage stimuli were expanded in willow species from high-altitude regions in southwestern China, which may contribute to their tolerance to ultraviolet radiation stress. Our study generates a comprehensive pan-transcriptome resource for a large set of Salix species and provides insights into the adaptations of willows to diverse environments, which will be valuable for comparative analysis with other related woody and herbaceous plants.