A recent study led by Dr. Xu Shuhua from the CAS-MPG Partner Institute for Computational Biology, Shanghai Institute of Nutrition and Health of the Chinese Academy of Sciences (CAS), was published online in National Science Review on August 7, 2019. It was entitled “Prioritizing natural selection signals from the deep-sequencing genomic data suggests multi-variant adaptation in Tibetan highlanders”.

This study constructed a map of adaptive genetic variants of Tibetan highlanders on the whole-genome scale. It is the very first effort for such a purpose based on the whole-genome deep sequencing data, and provides the first comprehensive panel of functional variants closely related to the adaptive evolution of the Tibetan people.

Signature of local adaptation at TMEM247-rs116983452-T and association analysis of its functional and phenotypic consequences [Image: Dr. Xu Shuhua’s team]

The last decade witnessed intensive studies of Tibetan high-altitude adaptation. EPAS1 is so far the most well-recognized candidate gene for high-altitude adaptation of Tibetans, as well as other human and non-human species on the highlands. However, no adaptive functional variant has been identified in EPAS1, leaving an unsolved problem in this field, and raising some more fundamental questions: (1) How many genes or genetic variants have driven the adaptive evolution of humans in high altitudes? (2) Does a causal variant located within the EPAS1 gene or anywhere else underlie the association between EPAS1 and the high-altitude adaptation? (3) Are there any other genes with greater contribution than EPAS1 to human adaptation to high altitudes?

To address these questions, researchers conducted a systematic analysis of the Tibetan genomes. Taking advantage of the whole-genome deep sequencing, they constructed a map of adaptive genetic variants of high-altitude adaptation, including 1,877 key variants with known functions: 63 missense, 7 loss-of-function, 1,298 evolutionarily conserved variants and 509 expression quantitative traits loci. These functional variants are very likely to have contributed to the adaptive evolution of the Tibetans, but not necessarily in a direct way. High-altitude adaptation involves a wide range of phenotypic variations driven by enormously large numbers of variants and genes; it could be even more complicated than some complex diseases. Moreover, the research team developed a statistic (FIS, functional importance score) to prioritize these identified adaptive genetic variants, and found that the top signal is not the well-known EPAS1, but TMEM247, a transmembrane protein coding gene. Especially, a missense variant (rs116983452) located in TMEM247 showed rather high frequency specifically in Tibetans. This key variant may largely increase the genetic differentiation between the lowlanders and Tibetan highlanders, by transforming Alanine, the common wild type in the lowlanders, to Valine, the Tibetan-specific mutant type. This mutant is carried by around 94 percent of Tibetans, but is in low frequency or even missing in other world-wide modern human populations. It is so far the most highly differentiated missense variant between Tibetans and the lowlanders. Interestingly, this mutant also presents as a homozygote in a 50,000-year-old human genome that was found in the Denisova Cave in Siberia. Researchers estimated the age of the adaptive Tibetan sequences carrying TMEM247-rs116983452-T to be about 60,000 years, implying that this Tibetan-specific variant could be inherited from early inhabitants of the Tibet Plateau with archaic ancestry.

In fact, it was a long journey for human to conquer the Tibet Plateau. A previous study by Xu’s team estimated that the genetic origin of the Tibetan highlanders could be traced back to around 40,000 to 60,000 years ago, in the middle-late Paleolithic. The early migrants to the plateau had extensive genetic admixture with each other, and had further gene flows with the latecomers, leading to admixed descendants with very complex genetic makeup -- inherited from ancestral lineages of modern human and archaic hominins (e.g. Altai Neanderthal, Denisovan and other unknown archaic species). During this process, some archaic genomic segments that were advantageous to the high-altitude adaptation were retained and accumulated to high frequencies through natural selection. TMEM247-rs116983452-T is a typical case of such a scenario.

The team found that the frequency of TMEM247-rs116983452-T is strongly and positively correlated with altitude. Further analyses revealed that this variant is significantly associated with the expression levels of EPAS1 and TMEM247, and could possibly regulate the concentrations of hemoglobin and red blood cells of Tibetans under hypoxia. For most lowland people, a long period of exposure to hypoxia may induce the production of red blood cells to increase blood oxygen, but may finally lead to polycythemia. In contrast, the hemoglobin and red blood cells are kept in relatively lower levels in the Tibetan highlanders, and TMEM247-rs116983452-T could be one of the key genetic factors of this protective mechanism. By statistical modeling, the team found that TMEM247-rs116983452 shows greater effect size and better predicts the phenotypic outcome than any EPAS1 variants in the association with adaptive traits in Tibetans, but interactions were also observed between TMEM247-rs116983452 and EPAS1 variants, indicating that multiple variants may jointly be responsible for the fitness of the Tibetans on the Plateau, where a complex model is needed to elucidate the adaptive evolution mechanism. The adaptive genetic variants map provided by this study helps to narrow down the targets for further investigations on the genetic basis and molecular mechanism of high-altitude adaptation of Tibetans, and provides a new perspective for unrevealing the mystery of human conquest of extreme environments at high altitudes. Notably, this study proposed that multiple variants may jointly be the cause of the fitness of the Tibetans on the Plateau, where a complex model is needed to elucidate the adaptive evolution mechanism.

This work was conducted by Dr. Deng Lian, Dr. Zhang Chao, Dr. Yuan Kai, and PhD students Gao Yang (from ShanghaiTech University) and Pan Yuwen from Dr. Xu Shuhua’s team, in collaboration with researchers from Kunming Institute of Zoology of CAS, Wenzhou Medical University, Fudan University, Xizang Minzu University, and other institutions.

For more information, please contact:

Ms Wang Jin of the Shanghai Institute of Nutrition and Health,

Chinese Academy of Sciences

Email: sibssc@sibs.ac.cn

Source: Shanghai Institute of Nutrition and Health, CAS