Known as the ‘third polar region’, the Qinghai-Tibet Plateau represents one of the harshest highland environments in the world and yet a number of organisms thrive there. Previous studies of birds, animals and humans have focused on well-differentiated populations in later stages of phenotypic divergence. The adaptive processes during the initial phase of highland adaptation remain poorly understood. We studied a human commensal, the Eurasian Tree Sparrow, which has followed human agriculture to the Qinghai-Tibet Plateau. Despite strong phenotypic differentiation at multiple levels, in particular in muscle-related phenotypes, highland and lowland populations show shallow genomic divergence and the colonization event occurred within the past few thousand years. In a one-month acclimation experiment investigating phenotypic plasticity, we exposed adult lowland tree sparrows to a hypoxic environment and did not observe muscle changes. Through population genetic analyses, we identified a signature of polygenic adaptation, whereby shifts in allele frequencies are spread across multiple loci, many of which are associated with muscle-related processes. Our results reveal a case of positive selection in which polygenic adaptation appears to drive rapid phenotypic evolution, shedding light on early stages of adaptive evolution to a novel environment.

Phylogenetic trees based on genome-wide sequence data may not always represent the true evolutionary history for a variety of reasons. One process that can lead to incorrect reconstruction of species phylogenies is gene flow, especially if interspecific gene flow has affected large parts of the genome. We investigated phylogenetic relationships within a clade comprising eight species of passerine birds (Phylloscopidae, Phylloscopus, leaf warblers) using one de novo genome assembly and 78 resequenced genomes. On the basis of hypothesis-exclusion trials based on D-statistics, phylogenetic network analysis and demographic inference analysis, we identified ancient gene flow affecting large parts of the genome between one species and the ancestral lineage of a sister species pair. This ancient gene flow consistently caused erroneous reconstruction of the phylogeny when using large amounts of genome-wide sequence data. In contrast, the true relationships were captured when smaller parts of the genome were analyzed, showing that the “winner-takes-all democratic majority tree” is not necessarily the true species tree. Under this condition, smaller amounts of data may sometimes avoid the effects of gene flow due to stochastic sampling, as hidden reticulation histories are more likely to emerge from the use of larger datasets, especially whole-genome datasets. In addition, we also found that genomic regions affected by ancient gene flow generally exhibited higher genomic differentiation but a lower recombination rate and nucleotide diversity. Our study highlights the importance of considering reticulation in phylogenetic reconstructions in the genomic era.

Species in a shared environment tend to evolve similar adaptations under the influence of their phylogenetic context. Using snowfinches, a monophyletic group of passerine birds (Passeridae), we study the relative roles of ancestral and species-specific adaptations to an extreme high-elevation environment, the Qinghai–Tibet Plateau. Our ancestral trait reconstruction shows that the ancestral snowfinch occupied high elevations and had a larger body mass than most nonsnowfinches in Passeridae. Subsequently, this phenotypic adaptation diversified in the descendant species. By comparing high-quality genomes from representatives of the three phylogenetic lineages, we find that about 95% of genes under positive selection in the descendant species are different from those in the ancestor. Consistently, the biological functions enriched for these species differ from those of their ancestor to various degrees (semantic similarity values ranging from 0.27 to 0.5), suggesting that the three descendant species have evolved divergently from the initial adaptation in their common ancestor. Using a functional assay to a highly selective gene, DTL, we demonstrate that the nonsynonymous substitutions in the ancestor and descendant species have improved the repair capacity of ultraviolet-induced DNA damage. The repair kinetics of the DTL gene shows a twofold to fourfold variation across the ancestor and the descendants. Collectively, this study reveals an exceptional case of adaptive evolution to high-elevation environments, an evolutionary process with an initial adaptation in the common ancestor followed by adaptive diversification of the descendant species.