Abstract
Understanding how species adapt to extreme environments is an extension of the main goals of evolutionary biology. While alpine plants are an ideal system for investigating the genetic basis of high-altitude adaptation, genomic resources in these species are still limited. In the present study, we generated reference-level transcriptomic data of five Saussurea species through high-throughput sequencing and de novo assembly. Three of them are located in the highland of the Qinghai-Tibet Plateau (QTP), and the other two are close relatives distributed in the lowland. A series of comparative and evolutionary genomics analyses were conducted to explore the genetic signatures of adaptive evolution to high-altitude environments. Estimation of divergence time using single-copy orthologs revealed that Saussurea species diversified during the Miocene, a period with extensive tectonic movement and climatic fluctuation on the QTP. We characterized gene families specific to the alpine species, including genes involved in oxidoreductase activity, pectin catabolic process, lipid transport, and polysaccharide metabolic process, which may play important roles in defense of hypoxia and freezing temperatures of the QTP. Furthermore, in a phylogenetic context with the branch model, we identified hundreds of genes with signatures of positive selection. These genes are involved in DNA repair, membrane transport, response to UV-B and hypoxia, and reproductive processes, as well as some metabolic processes associated with nutrient intake, potentially responsible for Saussurea adaptation to the harsh environments of high altitude. Overall, our study provides valuable genomic resources for alpine species and gained helpful insights into the genomic basis of plants adapting to extreme environments.