Structural variants (SVs) may play important roles in human adaptation to extreme environments such as high altitude but have been under-investigated. Here, combining long-read sequencing with multiple scaffolding techniques, we assembled a high-quality Tibetan genome (ZF1), with a contig N50 length of 24.57 mega-base pairs (Mb) and a scaffold N50 length of 58.80 Mb. The ZF1 assembly filled 80 remaining N-gaps (0.25 Mb in total length) in the reference human genome (GRCh38). Markedly, we detected 17 900 SVs, among which the ZF1-specific SVs are enriched in GTPase activity that is required for activation of the hypoxic pathway. Further population analysis uncovered a 163-bp intronic deletion in the MKL1 gene showing large divergence between highland Tibetans and lowland Han Chinese. This deletion is significantly associated with lower systolic pulmonary arterial pressure, one of the key adaptive physiological traits in Tibetans. Moreover, with the use of the high-quality de novo assembly, we observed a much higher rate of genome-wide archaic hominid (Altai Neanderthal and Denisovan) shared non-reference sequences in ZF1 (1.32%–1.53%) compared to other East Asian genomes (0.70%–0.98%), reflecting a unique genomic composition of Tibetans. One such archaic hominid shared sequence—a 662-bp intronic insertion in the SCUBE2 gene—is enriched and associated with better lung function (the FEV1/FVC ratio) in Tibetans. Collectively, we generated the first high-resolution Tibetan reference genome, and the identified SVs may serve as valuable resources for future evolutionary and medical studies.

Nitric oxide (NO) is an important molecule for vasomotor tone, and elevated NO signaling was previously hypothesized as a unique and adaptive physiological change in highland Tibetans. However, there has been lack of NO data from Tibetans living at low altitude and lowlander immigrants living at high altitude, which is crucial to test this hypothesis. Here, through cross-altitude (1990–5018 m) and cross-population (Tibetans and Han Chinese) analyses of serum NO metabolites (NOx) of 2086 individuals, we demonstrate that although Tibetans have a higher serum NOx level compared to lowlanders, Han Chinese immigrants living at high altitude show an even higher level than Tibetans. Consequently, our data contradict the previous proposal of increased NO signaling as the unique adaptive strategy in Tibetans. Instead, Tibetans have a relatively lower circulating NOx level at high altitude. This observation is further supported by data from the hypoxic experiments using human umbilical vein endothelial cells and gene knockout mice. No difference is detected between Tibetans and Han Chinese for endothelial nitric oxide synthase (eNOS), the key enzyme for circulating NO synthesis, suggesting that eNOS itself is unlikely to be the cause. We show that other NO synthesis-related genes (e.g. GCH1) carry Tibetan-enriched mutations significantly associated with the level of circulating NOx in Tibetans. Furthermore, gene network analysis revealed that the downregulation and upregulation of NOx is possibly achieved through distinct pathways. Collectively, our findings provide novel insights into the physiological and genetic mechanisms of the evolutionary adaptation of Tibetans to high-altitude hypoxia.