Genetic evidence of paleolithic colonization and neolithic expansion of modern humans on the tibetan plateau

Why and when was lactase persistence selected for? Insights from Central Asian herders and ancient DNA

The genetic adaptation of humans to the consumption of milk from dairying animals is one of the most emblematic cases of recent human evolution. While the phenotypic change under selection, lactase persistence (LP), is known, the evolutionary advantage conferred to persistent individuals remains obscure. One informative but underappreciated observation is that not all populations whose ancestors had access to milk genetically adapted to become lactase persistent. Indeed, Central Asian herders are mostly lactase nonpersistent, despite their significant dietary reliance on dairy products. Investigating the temporal dynamic of the -13.910:C>T Eurasian mutation associated with LP, we found that, after its emergence in Ukraine 5,960 before present (BP), the T allele spread between 4,000 BP and 3,500 BP throughout Eurasia, from Spain to Kazakhstan. The timing and geographical progression of the mutation coincides well with the migration of steppe populations across and outside of Europe. After 3,000 BP, the mutation strongly increased in frequency in Europe, but not in Asia. We propose that Central Asian herders have adapted to milk consumption culturally, by fermentation, and/or by colonic adaptation, rather than genetically. Given the possibility of a nongenetic adaptation to avoid intestinal symptoms when consuming dairy products, the puzzle then becomes this: why has LP been selected for at all?

Genome-Wide Analysis Reveals Changes in Polled Yak Long Non-coding RNAs in Skeletal Muscle Development

Long non-coding RNAs (lncRNAs) have been extensively studied in recent years. Numerous lncRNAs have been identified in mice, rats, and humans, some of which play important roles in muscle formation and development. However, little is known about lncRNA regulators that affect muscle development in yak (Bos grunniens). LncRNA expression during skeletal muscle development in yak was analyzed by RNA sequencing at three development stages: 3 years (group A), 6 months (group M), and 90-day-old fetuses (group E). A total of 1180 lncRNAs were identified in the three development stages. Compared with group E, 154 were upregulated and 130 were downregulated in group A. Compared with group A, 31 were upregulated and 29 were downregulated in group M. Compared with group E, 147 were upregulated and 149 were downregulated in group M (padj < 0.001, |log2FC| > 1.2). In addition, functional annotation analysis based on gene ontology (GO) and the Kyoto protocol encyclopedia of genes and genomes (KEGG) database showed that differentially expressed lncRNAs (DElncRNAs) were cis–trans target genes. The results showed that DElncRNAs were mainly involved in PI3K-Akt signaling pathway, focal adhesion, MAPK signaling pathway, apoptosis, and p53 signaling pathway. Furthermore, RTL1, IGF2, MEF2C, Pax7, and other well-known muscle development regulators were included in a co-expression network of differentially expressed target genes and lncRNAs. These data will help to further clarify the function of lncRNAs in the different stages of skeletal muscle developmental in yak.

Paternal heritage of the Han Chinese in Henan province (Central China): high diversity and evidence of in situ Neolithic expansions

Background: Due to their long history, complex admixture processes and large population sizes, more research is required to explore the fine genetic structure of Han populations from different geographic locations of China.Aim: To characterise the paternal genetic structure of the Han Chinese in Henan province, which was once the central living region of the ancient Huaxia population, the precursors of the Han Chinese.Subjects and methods: We sequenced Y chromosomes of 60 males from Zhengzhou, Henan Province, and reconstructed a phylogenetic tree for these samples with age estimation.Results: We observed high diversity of paternal lineages in our collection. We found that the in situ Neolithic expansion of the "Major lineages" contributed to a large portion of the paternal gene pool of the Han population in Henan Province. We also detected a large number of "Minor lineages" that diverged in the Palaeolithic Age.Conclusion: We suggest that the high genetic diversity in the paternal gene pool of modern Han populations is mainly attributed to the reservation of a larger number of lineages that diverged in the Palaeolithic Age, while the recent expansion of limited lineages contributed to the majority of the gene pool of modern Han populations. We propose that such a structure is a basal characteristic for the genetic structure of modern Han populations.

Genetic diversity and phylogenetic structure of four Tibeto-Burman-speaking populations in Tibetan-Yi corridor revealed by insertion/deletion polymorphisms

BACKGROUND

Insertion/deletion polymorphisms (InDels), combined with all the desirable features of both short tandem repeat and single nucleotide polymorphism, have been used in archaeological and anthropological research, population genetics and forensic application.

METHODS

Thirty InDels in 530 individuals residing in the Tibetan-Yi corridor (142 Dujiangyan Tibetans, 164 Muli Tibetans, 187 Xichang Yis, and 37 Yanyuan Mosuos) were genotyped using the Investigator DIPplex. Forensic parameters and allele frequency spectrum were calculated. Genetic relationships between the investigated populations and worldwide and nationwide populations were assessed based on both the allele frequency distribution and genotype data.

RESULTS

The combined powers of exclusion were 0.9807 (Dujiangyan Tibetan), 0.9880 (Muli Tibetan), 0.9852 (Xichang Yi) and 0.9892 (Yanyuan Mosuo). The combined powers of discrimination were 0.999999999983 (Dujiangyan Tibetan), 0.999999999942 (Muli Tibetan), 0.999999999982 (Xichang Yi) and 0.999999999962 (Yanyuan Mosuo), respectively. The comprehensive population comparisons among worldwide and nationwide populations uniformly illustrated that the investigated populations have a genetically closer relationship with Tibeto-Burman-speaking populations and geographically adjacent populations.

CONCLUSION

These 30 loci can be regarded as an efficient genetic tool in forensic individual identification and as a supplementary tool in paternity testing in Dujiangyan Tibetan, Muli Tibetan, Xichang Yi, and Yanyuan Mosuo. The genetic proximity between the four populations in the Tibetan-Yi corridor and other populations is strongly correlated with the linguistic origin and geographical distance.

Ancient mitogenomes show plateau populations from last 5200 years partially contributed to present-day Tibetans

The clarification of the genetic origins of present-day Tibetans requires an understanding of their past relationships with the ancient populations of the Tibetan Plateau. Here we successfully sequenced 67 complete mitochondrial DNA genomes of 5200 to 300-year-old humans from the plateau. Apart from identifying two ancient plateau lineages (haplogroups D4j1b and M9a1a1c1b1a) that suggest some ancestors of Tibetans came from low-altitude areas 4750 to 2775 years ago and that some were involved in an expansion of people moving between high-altitude areas 2125 to 1100 years ago, we found limited evidence of recent matrilineal continuity on the plateau. Furthermore, deep learning of the ancient data incorporated into simulation models with an accuracy of 97% supports that present-day Tibetan matrilineal ancestry received partial contribution rather than complete continuity from the plateau populations of the last 5200 years.

Genetic investigation and phylogenetic analysis of three Chinese ethnic groups using 16 X chromosome STR loci

Background: The value of using X-chromosomal short tandem repeats (X-STRs) as genetic markers in human genetics has been widely recognised. However, the 16 X-STRs in the Goldeneye^® DNA ID System 17X kit have not been thoroughly applied.Aim: To investigate the genetic polymorphisms of 16 X-STRs in three main ethnic minorities (Tibetan, Mongolian and Kazakh) in China and to reveal the phylogenetic relationships of different populations.Subjects and methods: A total of 245 Tibetan, 168 Mongolian and 105 Kazakh individuals were genotyped using this 17X kit. The allelic frequencies and other parameters were calculated. An additional eight Chinese populations and nine global populations were included in genetic comparisons based on 16 or 8 overlapped X-STRs.Results: A total of 147 alleles were observed from 16 X-STRs with allelic frequencies ranging from 0.0024 to 0.7952 in the three studied groups. Based on 16 X-STRs, Tibetans, Kazakhs and Mongolians showed more similarity to each other and were genetically distinct from the Shanghai Han group; based on 8 X-STRs, only the genetic relationships between different nations could be clarified.Conclusions: Our study presents an extensive report on a novel X-STR assay in three Chinese ethnic groups and a comprehensive genetic comparison between different populations based on these X-STRs.

Population genetics, diversity, forensic characteristics of four Chinese populations inferred from X-chromosomal short tandem repeats

We genotyped 19 X-STRs included in the AGCU X19 kit in 712 unrelated Chinese individuals from four populations (Li, Hui, Tibetan, and Han) and then merged with 4156 previously genotyped subjects from 19 populations to investigate genetic relationships and population structure, as well as the association between the genetic affinity and language classification. The combined mean paternity exclusion chances and the combined probabilities of discrimination based on single allele frequencies of 19-X-STRs and haplotype frequencies of seven linkage groups were high, indicating that this set of markers was very polymorphic in the four populations studied. Therefore, this panel can complement autosomal or uniparental markers in kinship analysis and complex deficient paternity testing. Subsequently, population differentiation analyses among 23 populations based on 19 STRs and 15 populations based on over 62 million single nucleotide polymorphisms consistently demonstrated that genetic stratifications exist between the different language-speaking populations, especially Tibeto-Burman-speaking, Tungusic and Turkic-speaking populations. Our newly studied populations are genetically close to ethnolinguistically adjacent populations. Our datasets can and should be used as an allele and haplotype frequencies reference database to facilitate the use of 19-X-STRs panel in routine forensic practice.

Paternal Y chromosomal genotyping reveals multiple large-scale admixtures in the formation of Lolo-Burmese-speaking populations in southwest China

Background: Bai and Yi people are two Tibeto-Burman speaking ethnic groups in Yunnan, southwest China. The genetic structure and history of these two groups are largely unknown due to a lack of available genetic data.Aim: To investigate the paternal genetic structure and population relationship of the Yi and Bai people.Subjects and methods: We collected samples from 278 Bai individuals and 283 Yi individuals from Yunnan and subsequently genotyped 43 phylogenetically relevant Y-SNPs in those samples. We estimated haplogroup frequencies and merged our data with a reference database including 46 representative worldwide populations to infer genetic relationships.Results: Y chromosomal haplogroup O-M175 is the dominant lineage in both Bai and Yi people. The Bai and Yi show a close genetic relationship with other Tibeto-Burman-speaking populations with high frequencies of haplogroup O2a2b1a1-Page23, which is also confirmed by PCA. The frequencies of the Tai-Kadai specific lineage O1a-M119, the southern China widespread lineage O1b-P31 and the eastern China enriched lineage O2a1b-002611, are also relatively high in our studied populations.Conclusions: The paternal Y chromosomal affinity of the Bai and Yi with Tibeto-Burman groups is consistent with the language classification. During the formation of the Bai and Yi populations, there were multiple large-scale admixtures, including the expansion of Neolithic farming populations from northern China, the assimilation of Tai-Kadai-speaking populations in southwest China, the demographic expansion driven by Neolithic agricultural revolution from southern China, and the admixture with populations of military immigration from northern and eastern China.

Transcriptome and DNA Methylation Analyses of the Molecular Mechanisms Underlying with Longissimus dorsi Muscles at Different Stages of Development in the Polled Yak

DNA methylation modifications are implicated in many biological processes. As the most common epigenetic mechanism DNA methylation also affects muscle growth and development. The majority of previous studies have focused on different varieties of yak, but little is known about the epigenetic regulation mechanisms in different age groups of animals. The development of muscles in the different stages of yak growth remains unclear. In this study, we selected the longissimus dorsi muscle tissue at three different growth stages of the yak, namely, 90-day-old fetuses (group E), six months old (group M), and three years old (group A). Using RNA-Seq transcriptome sequencing and methyl-RAD whole-genome methylation sequencing technology, changes in gene expression levels and DNA methylation status throughout the genome were investigated during the stages of yak development. Each group was represented by three biological replicates. The intersections of expression patterns of 7694 differentially expressed genes (DEGs) were identified (padj < 0.01, |log2FC| > 1.2) at each of the three developmental periods. Time-series expression profile clustering analysis indicated that the DEGs were significantly arranged into eight clusters which could be divided into two classes (padj < 0.05), class I profiles that were downregulated and class II profiles that were upregulated. Based on this cluster analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that DEGs from class I profiles were significantly (padj < 0.05) enriched in 21 pathways, the most enriched pathway being the Axon guidance signaling pathway. DEGs from the class II profile were significantly enriched in 58 pathways, the pathway most strongly enriched being Metabolic pathway. After establishing the methylation profiles of the whole genomes, and using two groups of comparisons, the three combinations of groups (M-vs.-E, M-vs.-A, A-vs.-E) were found to have 1344, 822, and 420 genes, respectively, that were differentially methylated at CCGG sites and 2282, 3056, and 537 genes, respectively, at CCWGG sites. The two sets of data were integrated and the negative correlations between DEGs and differentially methylated promoters (DMPs) analyzed, which confirmed that TMEM8C, IGF2, CACNA1S and MUSTN1 were methylated in the promoter region and that expression of the modified genes was negatively correlated. Interestingly, these four genes, from what was mentioned above, perform vital roles in yak muscle growth and represent a reference for future genomic and epigenomic studies in muscle development, in addition to enabling marker-assisted selection of growth traits.

Whole mitochondrial genome analysis of highland Tibetan ethnicity using massively parallel sequencing

Mitochondrial DNA (mtDNA) is a key player in numerous multifaceted and intricate biological processes and plays a pivotal role in dissecting the peopling of different populations, due to its maternally inherited property and comparatively high mutation rate. In this study, 119 Tibetan individuals from the Muli Tibetan Autonomous County of China (average altitude above 3,000 m) were employed in mitochondrial genome (mitogenome) sequencing by massively parallel sequencing (MPS) techniques using the Precision ID mtDNA Whole Genome Panel on an Ion S5XL system. The dataset presented 88 distinct haplotypes, resulting in the haplotype diversity of 0.9909. The majority of haplotypes were assigned to East Asian lineages and the distribution of haplogroups of Muli Tibetan significantly differed from reference Tibetan populations. The maximum parsimony phylogeny reconstructed by 119 newly generated mitogenomes revealed 12 major Muli Tibetan lineages. Intriguingly, a Sherpa-specific sub-haplogroup A15c1 with the lack of mutations at 4216 and 15,924 was discerned in our dataset, which suggested that the maternal gene pool of Sherpas may derive from Tibetan populations. The shared haplogroups between Muli Tibetan and lowland Han Chinese hinted that these lineages may derive from non-Tibetans and have already differentiated before their arrival on the Tibetan Plateau. Furthermore, extensive pairwise population comparisons displayed that Muli Tibetan had a closer genetic relationship with ethnically or linguistically close Nyingtri Tibetan, Nyingtri Lhoba and Chamdo Tibetan populations. Genetic affinity was also observed between the Muli Tibetan and North Han Chinese. Collectively, the results generated in this study enriched the existing forensic mtDNA database and raised additional interest in the application of whole mitogenome sequencing in forensic investigations.

De novo assembly of a Tibetan genome and identification of novel structural variants associated with high-altitude adaptation

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.

Population History and Altitude-Related Adaptation in the Sherpa

The first ascent of Mount Everest by Tenzing Norgay and Sir Edmund Hillary in 1953 brought global attention to the Sherpa people and human performance at altitude. The Sherpa inhabit the Khumbu Valley of Nepal, and are descendants of a population that has resided continuously on the Tibetan plateau for the past ∼25,000 to 40,000 years. The long exposure of the Sherpa to an inhospitable environment has driven genetic selection and produced distinct adaptive phenotypes. This review summarizes the population history of the Sherpa and their physiological and genetic adaptation to hypoxia. Genomic studies have identified robust signals of positive selection across EPAS1, EGLN1, and PPARA, that are associated with hemoglobin levels, which likely protect the Sherpa from altitude sickness. However, the biological underpinnings of other adaptive phenotypes such as birth weight and the increased reproductive success of Sherpa women are unknown. Further studies are required to identify additional signatures of selection and refine existing Sherpa-specific adaptive phenotypes to understand how genetic factors have underpinned adaptation in this population. By correlating known and emerging signals of genetic selection with adaptive phenotypes, we can further reveal hypoxia-related biological mechanisms of adaptation. Ultimately this work could provide valuable information regarding treatments of hypoxia-related illnesses including stroke, heart failure, lung disease and cancer.

The paternal genetic structure of Jingpo and Dai in southwest China

Yunnan province harbours substantial genetic, cultural and linguistic diversity, with the largest number of Aborigines in China, but the relationship among these Aborigines remains enigmatic. This study genotyped 45 Y chromosomal single nucleotide polymorphisms (SNPs) of 500 males from two aboriginal cross-border populations, Jingpo and Dai, from Dehong, Yunnan. It is reported that Haplogroup O2a2b1a1-M117 is the dominant lineage in both Jingpo and Dai. The Jingpo people show affinity with Tibeto-Burman speaking populations with a relatively high frequency of Haplogroup D-M174, and the Dai people are generally genetically similar with Tai-Kadai speaking populations with high frequencies of Haplogroup O1a-M119 and O1b1a1a-M95, which is consistent with their language classification.

Neolithic millet farmers contributed to the permanent settlement of the Tibetan Plateau by adopting barley agriculture

The permanent human settlement of the Tibetan Plateau (TP) has been suggested to have been facilitated by the introduction of barley agriculture ∼3.6 kilo-years ago (ka). However, how barley agriculture spread onto the TP remains unknown. Given that the lower altitudes in the northeastern TP were occupied by millet cultivators from 5.2 ka, who also adopted barley farming ∼4 ka, it is highly possible that it was millet farmers who brought barley agriculture onto the TP ∼3.6 ka. To test this hypothesis, we analyzed mitochondrial DNA (mtDNA) from 8277 Tibetans and 58 514 individuals from surrounding populations, including 682 newly sequenced whole mitogenomes. Multiple lines of evidence, together with radiocarbon dating of cereal remains at different elevations, supports the scenario that two haplogroups (M9a1a1c1b1a and A11a1a), which are common in contemporary Tibetans (20.9%) and were probably even more common (40–50%) in early Tibetans prior to historical immigrations to the TP, represent the genetic legacy of the Neolithic millet farmers. Both haplogroups originated in northern China between 10.0–6.0 ka and differentiated in the ancestors of modern Tibetans ∼5.2–4.0 ka, matching the dispersal history of millet farming. By showing that substantial genetic components in contemporary Tibetans can trace their ancestry back to the Neolithic millet farmers, our study reveals that millet farmers adopted and brought barley agriculture to the TP ∼3.6–3.3 ka, and made an important contribution to the Tibetan gene pool.

Genetic structure and forensic characteristics of Tibeto-Burman-speaking Ü-Tsang and Kham Tibetan Highlanders revealed by 27 Y-chromosomal STRs

Culturally diverse Tibetans (Ü-Tsang, Kham and Ando) harboring a unique molecular mechanism that allows them to successfully adapt to hypoxic environments in the Qinghai-Tibet Plateau have been a subject of great interest in medical genetics, linguistics, archeology and forensic science. However, forensic characteristics and genetic variations of the Y-chromosomal 27-marker haplotype included in the Yfiler Plus system in the Ü-Tsang and Kham Tibeto-Burman-speaking Tibetans remain unexplored. Thus, we genotyped 27 Y-STRs in 230 Shigatse Ü-Tsang Tibetans (SUT) and 172 Chamdo Kham Tibetans (CKT) to investigate the forensic characterization and genetic affinity of Chinese Tibetan Highlanders. The haplotype diversities were 0.999962028 in SUT and 0.999796002 in CKT. Forensic diversity measures indicated that this 27-Y-STR amplification system is appropriate for routine forensic applications, such as identifying and separating unrelated males in deficiency paternity cases, male disaster victims and missing person identification and determining male components in sexual assault cases. Moreover, the genetic relationships among 63 worldwide populations (16,282 individuals), 16 Asian populations, and 21 Chinese populations were analyzed and reconstructed using principal component analysis, multidimensional scaling plots and a phylogenetic tree. Considerable genetic differences were observed between Tibetan populations and other geographically/ethnically diverse populations (Han Chinese). Our studied SUT and CKT have a genetically closer relationship with Gansu Ando Tibetans than with other Asians. In total, our analyses indicated that subpopulation structures exist among Asian and Chinese populations, and population-specific reference databases should be established for forensic applications.

Ethnogenetic analysis reveals that Kohistanis of Pakistan were genetically linked to west Eurasians by a probable ancestral genepool from Eurasian steppe in the bronze age

Despite the unique geographic, ethnic, social and cultural features of Kohistan in Pakistan, the origin and descent of Kohistanis remain still obscure. In an effort to address questions concerning the genetic structure, origin and genetic affinities of Kohistanis, we herein applied an ethnogenetic approach consisting on mitochondrial DNA (mtDNA) analysis and dental morphology analysis. We sequenced HVS1 of mtDNA, observed 14 haplotypes and assigned a total of 9 haplogroups belonging to macrolineages M (17%) and N (83%). Genetic diversity estimates in Kohistanis (Hd = 0.910 ± 0.014; Pi = 0.019 ± 0.001; θw = 0.019 ± 0.006) were similar to that of previous studies in other Pakistani populations. Overall, the analyses of dental morphology and mtDNA profile of Kohistanis resulted in similar findings. All the analyses indicate that Kohistanis share affinities to populations from Europe, Near East, Central Asia and South Asia. The Kohistani HVS1 haplotype 2 shares 100% identity to HVS1 haplotypes across the Europe. These results in light of recent insights into ancient genomics lead us to conclude that ancestry from Eurasian Steppe genetically linked Kohistanis to all these populations in the Bronze Age. This is consistent with linguistic evidence and also with the Indo-Aryan migration model for the peopling of South Asia.

The Kobresia pygmaea ecosystem of the Tibetan highlands - Origin, functioning and degradation of the world's largest pastoral alpine ecosystem: Kobresia pastures of Tibet

With 450,000 km²Kobresia (syn. Carex) pygmaea dominated pastures in the eastern Tibetan highlands are the world's largest pastoral alpine ecosystem forming a durable turf cover at 3000-6000 m a.s.l. Kobresia's resilience and competitiveness is based on dwarf habit, predominantly below-ground allocation of photo assimilates, mixture of seed production and clonal growth, and high genetic diversity. Kobresia growth is co-limited by livestock-mediated nutrient withdrawal and, in the drier parts of the plateau, low rainfall during the short and cold growing season. Overstocking has caused pasture degradation and soil deterioration over most parts of the Tibetan highlands and is the basis for this man-made ecosystem. Natural autocyclic processes of turf destruction and soil erosion are initiated through polygonal turf cover cracking, and accelerated by soil-dwelling endemic small mammals in the absence of predators. The major consequences of vegetation cover deterioration include the release of large amounts of C, earlier diurnal formation of clouds, and decreased surface temperatures. These effects decrease the recovery potential of Kobresia pastures and make them more vulnerable to anthropogenic pressure and climate change. Traditional migratory rangeland management was sustainable over millennia, and possibly still offers the best strategy to conserve and possibly increase C stocks in the Kobresia turf.

The Transcriptomic Landscape of Yaks Reveals Molecular Pathways for High Altitude Adaptation

Yak is one of the largest native mammalian species at the Himalayas, the highest plateau area in the world with an average elevation of >4,000 m above the sea level. Yak is well adapted to high altitude environment with a set of physiological features for a more efficient blood flow for oxygen delivery under hypobaric hypoxia. Yet, the genetic mechanism underlying its adaptation remains elusive. We conducted a cross-tissue, cross-altitude, and cross-species study to characterize the transcriptomic landscape of domestic yaks. The generated multi-tissue transcriptomic data greatly improved the current yak genome annotation by identifying tens of thousands novel transcripts. We found that among the eight tested tissues (lung, heart, kidney, liver, spleen, muscle, testis, and brain), lung and heart are two key organs showing adaptive transcriptional changes and >90% of the cross-altitude differentially expressed genes in lung display a nonlinear regulation. Pathways related to cell survival and proliferation are enriched, including PI3K-Akt, HIF-1, focal adhesion, and ECM–receptor interaction. These findings, in combination with the comprehensive transcriptome data set, are valuable to understanding the genetic mechanism of hypoxic adaptation in yak.

Population Genetic Diversity and Phylogenetic Characteristics for High-Altitude Adaptive Kham Tibetan Revealed by DNATyperTM 19 Amplification System

Tibetans residing in the high-altitude inhospitable environment have undergone significant natural selection of their genetic architecture. Recently, highly mutational autosomal short tandem repeats were widely used not only in the anthropology and population genetics to investigate the genetic structure and relationships, but also in the medical genetics to explore the pathogenesis of multiple genetic diseases and in the forensic science to identify individual and parentage relatedness. However, genetic variants and forensic efficiency of DNATyper^TM 19 amplification system and genetic background of Kham Tibetan remain uncharacterized. Thus, we genotyped 19 forensic genetic markers in 11,402 Kham Tibetans to gain insight into the genetic diversity of Chinese high-altitude adaptive population. Highly discriminating and polymorphic forensic measures were observed, which indicated that this new-developed DNATyper 19 PCR amplification is suitable for routine forensic identification purposes and Chinese national DNA database establishment. Pairwise genetic distances among the comprehensive population comparisons suggested that this high-altitude adaptive Kham Tibetan has genetically closer relationships with lowlanders of Tibeto-Burman-speaking populations (Chengdu Tibetan, Liangshan Tibetan, and Liangshan Yi). Genetic substructure analyses via phylogenetic reconstruction, principal component analysis, and multidimensional scaling analysis in both nationwide and worldwide contexts suggested that the genetic proximity exists along the linguistic, ethnic, and continental geographical boundary. Further studies with whole-genome sequencing of modern or archaic Kham Tibetans would be useful in reconstructing the Tibetan population history.

Investigating the genetic diversity and affinities of historical populations of Tibet

This study elucidates Y chromosome distribution patterns in the three general provincial populations of historical Tibet, Amdo (n = 88), Dotoe (n = 109) and U-Tsang (n = 153) against the backdrop of 37 Asian reference populations. The central aim of this study is to investigate the genetic affinities of the three historical Tibetan populations among themselves and to neighboring populations. Y-SNP and Y-STR profiles were assessed in these historical populations. Correspondence analyses (CA) were generated with Y-SNP haplogroup data. Y-STR haplotypes were determined and employed to generate multidimensional scaling (MDS) plots based on Rst distances. Frequency contour maps of informative Y haplogroups were constructed to visualize the distributions of specific chromosome types. Network analyses based on Y-STR profiles of individuals under specific Y haplogroups were generated to examine the genetic heterogeneity among populations. Average gene diversity values and other parameters of population genetics interest were estimated to characterize the populations. The Y chromosomal results generated in this study indicate that using two sets of markers (Y-SNP, and Y-STR) the three Tibetan populations are genetically distinct. In addition, U-Tsang displays the highest gene diversity, followed by Amdo and Dotoe. The results of this transcontinental biogeographical investigation also indicate various degrees of paternal genetic affinities among these three Tibetan populations depending on the type of loci (Y-SNP or Y-STR) analyzed. The CA generated with Y-SNP haplogroup data demonstrates that Amdo and U-Tsang are closer to each other than to any neighboring non-Tibetan group. In contrast, the MDS plot based on Y-STR haplotypes displays Rst distances that are much shorter between U-Tsang and its geographic nearby populations of Ladakh, Punjab, Kathmandu and Newar than between it and Amdo. Moreover, although Dotoe is isolated from all other groups using both types of marker systems, it lies nearer to the other Tibetan collections in the Y-SNP CA than in the Y-STR MDS plot. High resolution and shallow evolutionary time frames engendered by Y-STR based analyses may reflect a more recent demographic history than that delineated by the more conserved Y-SNP markers.

Demographic History and Genetic Adaptation in the Himalayan Region Inferred from Genome-Wide SNP Genotypes of 49 Populations

We genotyped 738 individuals belonging to 49 populations from Nepal, Bhutan, North India, or Tibet at over 500,000 SNPs, and analyzed the genotypes in the context of available worldwide population data in order to investigate the demographic history of the region and the genetic adaptations to the harsh environment. The Himalayan populations resembled other South and East Asians, but in addition displayed their own specific ancestral component and showed strong population structure and genetic drift. We also found evidence for multiple admixture events involving Himalayan populations and South/East Asians between 200 and 2,000 years ago. In comparisons with available ancient genomes, the Himalayans, like other East and South Asian populations, showed similar genetic affinity to Eurasian hunter-gatherers (a 24,000-year-old Upper Palaeolithic Siberian), and the related Bronze Age Yamnaya. The high-altitude Himalayan populations all shared a specific ancestral component, suggesting that genetic adaptation to life at high altitude originated only once in this region and subsequently spread. Combining four approaches to identifying specific positively selected loci, we confirmed that the strongest signals of high-altitude adaptation were located near the Endothelial PAS domain-containing protein 1 and Egl-9 Family Hypoxia Inducible Factor 1 loci, and discovered eight additional robust signals of high-altitude adaptation, five of which have strong biological functional links to such adaptation. In conclusion, the demographic history of Himalayan populations is complex, with strong local differentiation, reflecting both genetic and cultural factors; these populations also display evidence of multiple genetic adaptations to high-altitude environments.

Haplotype diversity in mitochondrial DNA reveals the multiple origins of Tibetan horse

The Tibetan horse is a species endemic to the Tibetan plateau, with considerable economic value in the region. However, we currently have little genetic evidence to verify whether the breed originated in Tibet or if it entered the area via an ancient migratory route. In the present study, we analyzed the hypervariable segment I sequences of mitochondrial DNA (mtDNA) in 2,050 horses, including 290 individuals from five Tibetan populations and 1,760 from other areas across Asia. Network analysis revealed multiple maternal lineages in the Tibetan horse. Component analysis of sub-lineage F3 indicated that it decreased in frequency from east to west, a trend reflected both southward and northward from Inner Mongolia. Analysis of population genetics showed that the Deqen horse of eastern Tibet was more closely related to the Ningqiang horse of northern China than to other Tibetan horses or the Yunnan horse. These results indicated that the Tibetan horse migrated first from Central Asia to Mongolia, moved south to eastern Tibet (near Deqen), then finally westward to other regions of Tibet. We also identified a novel lineage K that mainly comprises Tibetan and Yunnan horses, suggesting autochthonous domesticated origin for some Tibetan horse breeds from local wild horses. In conclusion, our study demonstrated that modern Tibetan horse breeds originated from the introgression of local wild horses with exotic domesticated populations outside China.

GCH1 plays a role in the high-altitude adaptation of Tibetans

Tibetans are well adapted to high-altitude hypoxia. Previous genome-wide scans have reported many candidate genes for this adaptation, but only a few have been studied. Here we report on a hypoxia gene ( GCH1, GTP-cyclohydrolase I), involved in maintaining nitric oxide synthetase (NOS) function and normal blood pressure, that harbors many potentially adaptive variants in Tibetans. We resequenced an 80.8 kb fragment covering the entire gene region of GCH1 in 50 unrelated Tibetans. Combined with previously published data, we demonstrated many GCH1 variants showing deep divergence between highlander Tibetans and lowlander Han Chinese. Neutrality tests confirmed a signal of positive Darwinian selection on GCH1 in Tibetans. Moreover, association analysis indicated that the Tibetan version of GCH1 was significantly associated with multiple physiological traits in Tibetans, including blood nitric oxide concentration, blood oxygen saturation, and hemoglobin concentration. Taken together, we propose that GCH1 plays a role in the genetic adaptation of Tibetans to high altitude hypoxia.

Human phylogeography and diversity

Homo sapiens phylogeography begins with the species' origin nearly 200 kya in Africa. First signs of the species outside Africa (in Arabia) are from 125 kya. Earliest dates elsewhere are now 100 kya in China, 45 kya in Australia and southern Europe (maybe even 60 kya in Australia), 32 kya in northeast Siberia, and maybe 20 kya in the Americas. Humans reached arctic regions and oceanic islands last-arctic North America about 5 kya, mid- and eastern Pacific islands about 2-1 kya, and New Zealand about 700 y ago. Initial routes along coasts seem the most likely given abundant and easily harvested shellfish there as indicated by huge ancient oyster shell middens on all continents. Nevertheless, the effect of geographic barriers-mountains and oceans-is clear. The phylogeographic pattern of diasporas from several single origins-northeast Africa to Eurasia, southeast Eurasia to Australia, and northeast Siberia to the Americas-allows the equivalent of a repeat experiment on the relation between geography and phylogenetic and cultural diversity. On all continents, cultural diversity is high in productive low latitudes, presumably because such regions can support populations of sustainable size in a small area, therefore allowing a high density of cultures. Of course, other factors operate. South America has an unusually low density of cultures in its tropical latitudes. A likely factor is the phylogeographic movement of peoples from the Old World bringing novel and hence, lethal diseases to the New World, a foretaste, perhaps, of present day global transport of tropical diseases.

The genomic landscape of Nepalese Tibeto-Burmans reveals new insights into the recent peopling of Southern Himalayas

While much research attention has focused on demographic processes that enabled human diffusion on the Tibetan plateau, little is known about more recent colonization of Southern Himalayas. In particular, the history of migrations, admixture and/or isolation of populations speaking Tibeto-Burman languages, which is supposed to be quite complex and to have reshaped patterns of genetic variation on both sides of the Himalayan arc, remains only partially elucidated. We thus described the genomic landscape of previously unsurveyed Tibeto-Burman (i.e. Sherpa and Tamang) and Indo-Aryan communities from remote Nepalese valleys. Exploration of their genomic relationships with South/East Asian populations provided evidence for Tibetan admixture with low-altitude East Asians and for Sherpa isolation. We also showed that the other Southern Himalayan Tibeto-Burmans derived East Asian ancestry not from the Tibetan/Sherpa lineage, but from low-altitude ancestors who migrated from China plausibly across Northern India/Myanmar, having experienced extensive admixture that reshuffled the ancestral Tibeto-Burman gene pool. These findings improved the understanding of the impact of gene flow/drift on the evolution of high-altitude Himalayan peoples and shed light on migration events that drove colonization of the southern Himalayan slopes, as well as on the role played by different Tibeto-Burman groups in such a complex demographic scenario.

Evolutionary significance of selected EDAR variants in Tibetan high-altitude adaptations

Humans have been exposed to many environmental challenges since their evolutionary origins in Africa and subsequent migrations to the rest of the world. A severe environmental challenge to human migrants was hypoxia caused by low barometric oxygen pressure at high altitudes. Several genome-wide scans have elucidated the genetic basis of human high-altitude adaptations. However, the dearth of functional variant information has led to the successful association of only a few candidate genes. In the present study, we employed a candidate gene approach and re-sequenced the EDAR locus in 45 Tibetan individuals to identify mutations involved in hypoxia adaptation. We identified 10 and five quantitative trait-associated mutations for oxygen saturation (SaO₂) and blood platelet count, respectively, at the EDAR locus. Among these, rs10865026 and rs3749110 (associated with SaO₂ and platelet count, respectively) were identified as functional candidate targets. These data demonstrate that EDAR has undergone natural selection in recent human history and indicate an important role of EDAR variants in Tibetan high-altitude adaptations.

Permanent human occupation of the central Tibetan Plateau in the early Holocene

Current models of the peopling of the higher-elevation zones of the Tibetan Plateau postulate that permanent occupation could only have been facilitated by an agricultural lifeway at ~3.6 thousand calibrated carbon-14 years before present. Here we report a reanalysis of the chronology of the Chusang site, located on the central Tibetan Plateau at an elevation of ~4270 meters above sea level. The minimum age of the site is fixed at ~7.4 thousand years (thorium-230/uranium dating), with a maximum age between ~8.20 and 12.67 thousand calibrated carbon-14 years before present (carbon-14 assays). Travel cost modeling and archaeological data suggest that the site was part of an annual, permanent, preagricultural occupation of the central plateau. These findings challenge current models of the occupation of the Tibetan Plateau.

Down-Regulation of EPAS1 Transcription and Genetic Adaptation of Tibetans to High-Altitude Hypoxia

Tibetans are well adapted to the hypoxic environments at high altitude, yet the molecular mechanism of this adaptation remains elusive. We reported comprehensive genetic and functional analyses of EPAS1, a gene encoding hypoxia inducible factor 2α (HIF-2α) with the strongest signal of selection in previous genome-wide scans of Tibetans. We showed that the Tibetan-enriched EPAS1 variants down-regulate expression in human umbilical endothelial cells and placentas. Heterozygous EPAS1 knockout mice display blunted physiological responses to chronic hypoxia, mirroring the situation in Tibetans. Furthermore, we found that the Tibetan version of EPAS1 is not only associated with the relatively low hemoglobin level as a polycythemia protectant, but also is associated with a low pulmonary vasoconstriction response in Tibetans. We propose that the down-regulation of EPAS1 contributes to the molecular basis of Tibetans’ adaption to high-altitude hypoxia.

EPAS1 variants in high altitude Tibetan wolves were selectively introgressed into highland dogs

BACKGROUND

Admixture can facilitate adaptation. For example, black wolves have obtained the variant causing black coat color through past hybridization with domestic dogs and have higher fitness than gray colored wolves. Another recent example of the transfer of adaptive variation between the two species has been suggested by the similarity between high altitude Tibetan mastiffs and wolves at the EPAS1 gene, a transcription factor induced in low oxygen environments.

METHODS

Here, we investigate the directionality of admixture in EPAS1 between 28 reference highland gray wolves, 15 reference domestic dogs, and 21 putatively admixed highland wolves. This experimental design represents an expanded sample of Asian dogs and wolves from previous studies. Admixture was inferred using 17,709 publicly available SNP genotypes on canine chromosome 10. We additionally conducted a scan for positive selection in the highland dog genome.

RESULTS

We find an excess of highland gray wolf ancestry at the EPAS1 locus in highland domestic dogs, suggesting adaptive introgression from wolves to dogs. The signal of admixture is limited in genomic extent to a small region on chromosome 10, indicating that it is the focus of selection in an oxygen-limited environment.

DISCUSSION

Our results suggest that an adaptive variant of EPAS1 in highland wolves was transferred to highland dogs, carrying linked variants that potentially function in hypoxia response at high elevation. The intertwined history of dogs and wolves ensures a unique evolutionary dynamic where variants that have appeared in the history of either species can be tested for their effects on fitness under natural and artificial selection. Such coupled evolutionary histories may be key to the persistence of wild canines and their domesticated kin given the increasing anthropogenic modifications that characterize the future of both species.

Differentiated demographic histories and local adaptations between Sherpas and Tibetans

BACKGROUND

The genetic relationships reported by recent studies between Sherpas and Tibetans are controversial. To gain insights into the population history and the genetic basis of high-altitude adaptation of the two groups, we analyzed genome-wide data in 111 Sherpas (Tibet and Nepal) and 177 Tibetans (Tibet and Qinghai), together with available data from present-day human populations.

RESULTS

Sherpas and Tibetans show considerable genetic differences and can be distinguished as two distinct groups, even though the divergence between them (~3200-11,300 years ago) is much later than that between Han Chinese and either of the two groups (~6200-16,000 years ago). Sub-population structures exist in both Sherpas and Tibetans, corresponding to geographical or linguistic groups. Differentiation of genetic variants between Sherpas and Tibetans associated with adaptation to either high-altitude or ultraviolet radiation were identified and validated by genotyping additional Sherpa and Tibetan samples.

CONCLUSIONS

Our analyses indicate that both Sherpas and Tibetans are admixed populations, but the findings do not support the previous hypothesis that Tibetans derive their ancestry from Sherpas and Han Chinese. Compared to Tibetans, Sherpas show higher levels of South Asian ancestry, while Tibetans show higher levels of East Asian and Central Asian/Siberian ancestry. We propose a new model to elucidate the differentiated demographic histories and local adaptations of Sherpas and Tibetans.

Revisiting the male genetic landscape of China: a multi-center study of almost 38,000 Y-STR haplotypes

China has repeatedly been the subject of genetic studies to elucidate its prehistoric and historic demography. While some studies reported a genetic distinction between Northern and Southern Han Chinese, others showed a more clinal picture of small differences within China. Here, we investigated the distribution of Y chromosome variation along administrative as well as ethnic divisions in the mainland territory of the People's Republic of China, including 28 administrative regions and 19 recognized Chinese nationalities, to assess the impact of recent demographic processes. To this end, we analyzed 37,994 Y chromosomal 17-marker haplotype profiles from the YHRD database with respect to forensic diversity measures and genetic distance between groups defined by administrative boundaries and ethnic origin. We observed high diversity throughout all Chinese provinces and ethnicities. Some ethnicities, including most prominently Kazakhs and Tibetans, showed significant genetic differentiation from the Han and other groups. However, differences between provinces were, except for those located on the Tibetan plateau, less pronounced. This discrepancy is explicable by the sizeable presence of Han speakers, who showed high genetic homogeneity all across China, in nearly all studied provinces. Furthermore, we observed a continuous genetic North-South gradient in the Han, confirming previous reports of a clinal distribution of Y chromosome variation and being in notable concordance with the previously observed spatial distribution of autosomal variation. Our findings shed light on the demographic changes in China accrued by a fast-growing and increasingly mobile population.

A longitudinal cline characterizes the genetic structure of human populations in the Tibetan plateau

Indigenous populations of the Tibetan plateau have attracted much attention for their good performance at extreme high altitude. Most genetic studies of Tibetan adaptations have used genetic variation data at the genome scale, while genetic inferences about their demography and population structure are largely based on uniparental markers. To provide genome-wide information on population structure, we analyzed new and published data of 338 individuals from indigenous populations across the plateau in conjunction with worldwide genetic variation data. We found a clear signal of genetic stratification across the east-west axis within Tibetan samples. Samples from more eastern locations tend to have higher genetic affinity with lowland East Asians, which can be explained by more gene flow from lowland East Asia onto the plateau. Our findings corroborate a previous report of admixture signals in Tibetans, which were based on a subset of the samples analyzed here, but add evidence for isolation by distance in a broader geospatial context.

Improved phylogenetic resolution for Y-chromosome Haplogroup O2a1c-002611

Y-chromosome Haplogroup O2a1c-002611 is one of the dominant lineages of East Asians and Southeast Asians. However, its internal phylogeny remains insufficiently investigated. In this study, we genotyped 89 new highly informative single nucleotide polymorphisms (SNPs) in 305 individuals with Haplogroup O2a1c-002611 identified from 2139 Han Chinese males. Two major branches were identified, O2a1c1-F18 and O2a1c2-L133.2 and the first was further divided into two main subclades, O2a1c1a-F11 and O2a1c1b-F449, accounting for 11.13% and 2.20% of Han Chinese, respectively. In Haplogroup O2a1c1a-F11, we also determined seven sublineages with quite different frequency distributions in Han Chinese ranging from 0.187% to 3.553%, implying they might have different demographic history. The reconstructed haplogroup tree for all the major clades within Haplogroup O2a1c-002611 permits better resolution of male lineages in population studies of East Asia and Southeast Asia. The dataset generated in the present study are also valuable for forensic identification and paternity tests in China.

Habitat Variability and Ethnic Diversity in Northern Tibetan Plateau

There are 56 officially-recognized ethnic groups in China. However, the distinct geographic patterns of various ethnic groups in relation to the physical environment in China have rarely been investigated. Based on the geo-referenced physical environmental parameters of 455 Han, Tu, Hui, Salar, Mongolian, and Tibetan communities in Qinghai, we found that the communities could be statistically demarcated by temperature and aridity threshold according to their ethnicity, implying that the geographic distribution of each ethnic group is mediated by the physical environment. We also observed that the habitat of each ethnic group is ecologically compatible with current subsistence strategies. Tibetans settle in cold and humid high-altitude regions owing to the cultivation of highland barley and the breeding of yak, dzo, Tibetan sheep and Tibetan goat. Mongolians survive by animal husbandry in cold and dry grassland areas. Han and Tu people settle in the Huangshui River Valley, which offers relatively humid climate and flat land for agriculture. Hui and Salar people occupy the Yellow River Valley with its relatively arid environment and grassland vegetation suitable for animal breeding. Our findings offer a new perspective in explaining the geographic patterns and the varieties of ethnic groups in China and elsewhere.

Genetic signatures of high-altitude adaptation in Tibetans

Indigenous Tibetan people have lived on the Tibetan Plateau for millennia. There is a long-standing question about the genetic basis of high-altitude adaptation in Tibetans. We conduct a genome-wide study of 7.3 million genotyped and imputed SNPs of 3,008 Tibetans and 7,287 non-Tibetan individuals of Eastern Asian ancestry. Using this large dataset, we detect signals of high-altitude adaptation at nine genomic loci, of which seven are unique. The alleles under natural selection at two of these loci [methylenetetrahydrofolate reductase (MTHFR) and EPAS1] are strongly associated with blood-related phenotypes, such as hemoglobin, homocysteine, and folate in Tibetans. The folate-increasing allele of rs1801133 at the MTHFR locus has an increased frequency in Tibetans more than expected under a drift model, which is probably a consequence of adaptation to high UV radiation. These findings provide important insights into understanding the genomic consequences of high-altitude adaptation in Tibetans.

Genetic structure of Tibetan populations in Gansu revealed by forensic STR loci

The origin and diversification of Sino-Tibetan speaking populations have been long-standing hot debates. However, the limited genetic information of Tibetan populations keeps this topic far from clear. In the present study, we genotyped 15 forensic autosomal short tandem repeats (STRs) from 803 unrelated Tibetan individuals from Gansu Province (635 from Gannan and 168 from Tianzhu) in northwest China. We combined these data with published dataset to infer a detailed population affinities and genetic substructure of Sino-Tibetan populations. Our results revealed Tibetan populations in Gannan and Tianzhu are genetically very similar with Tibetans from other regions. The Tibetans in Tianzhu have received more genetic influence from surrounding lowland populations. The genetic structure of Sino-Tibetan populations was strongly correlated with linguistic affiliations. Although the among-population variances are relatively small, the genetic components for Tibetan, Lolo-Burmese, and Han Chinese were quite distinctive, especially for the Deng, Nu, and Derung of Lolo-Burmese. Han Chinese but not Tibetans are suggested to share substantial genetic component with southern natives, such as Tai-Kadai and Hmong-Mien speaking populations, and with other lowland East Asian populations, which implies there might be extensive gene flow between those lowland groups and Han Chinese after Han Chinese were separated from Tibetans. The dataset generated in present study is also valuable for forensic identification and paternity tests in China.

Genetic structure in the Sherpa and neighboring Nepalese populations

Background

We set out to describe the fine-scale population structure across the Eastern region of Nepal. To date there is relatively little known about the genetic structure of the Sherpa residing in Nepal and their genetic relationship with the Nepalese. We assembled dense genotype data from a total of 1245 individuals representing Nepal and a variety of different populations resident across the greater Himalayan region including Tibet, China, India, Pakistan, Kazakhstan, Uzbekistan, Tajikistan and Kirghizstan. We performed analysis of principal components, admixture and homozygosity.

Results

We identified clear substructure across populations resident in the Himalayan arc, with genetic structure broadly mirroring geographical features of the region. Ethnic subgroups within Nepal show distinct genetic structure, on both admixture and principal component analysis. We detected differential proportions of ancestry from northern Himalayan populations across Nepalese subgroups, with the Nepalese Rai, Magar and Tamang carrying the greatest proportions of Tibetan ancestry.

Conclusions

We show that populations dwelling on the Himalayan plateau have had a clear impact on the Northern Indian gene pool. We illustrate how the Sherpa are a remarkably isolated population, with little gene flow from surrounding Nepalese populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3469-5) contains supplementary material, which is available to authorized users.

Characterization of mitochondrial DNA polymorphisms in the Han population in Liaoning Province, Northeast China

This study characterized the genetic variations of mitochondrial DNA (mtDNA) to elucidate the maternal genetic structure of Liaoning Han Chinese. A total of 317 blood samples of unrelated individuals were collected for analysis in Liaoning Province. The mtDNA samples were analyzed using two distinct methods: sequencing of the hypervariable sequences I and II (HVSI and HVSII), and polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) analysis of the coding region. The results indicated a high gene diversity value (0.9997 ± 0.0003), a high polymorphism information content (0.99668) and a random match probability (0.00332). These samples were classified into 305 haplotypes, with 9 shared haplotypes. The most common haplogroup was D4 (12.93%). The principal component analysis map, the phylogenetic tree map, and the genetic distance matrix all indicated that the genetic distance of the Liaoning Han population from the Tibetan group was distant, whereas that from the Miao group was relatively close.

Characterising private and shared signatures of positive selection in 37 Asian populations

The Asian Diversity Project (ADP) assembled 37 cosmopolitan and ethnic minority populations in Asia that have been densely genotyped across over half a million markers to study patterns of genetic diversity and positive natural selection. We performed population structure analyses of the ADP populations and divided these populations into four major groups based on their genographic information. By applying a highly sensitive algorithm haploPS to locate genomic signatures of positive selection, 140 distinct genomic regions exhibiting evidence of positive selection in at least one population were identified. We examined the extent of signal sharing for regions that were selected in multiple populations and observed that populations clustered in a similar fashion to that of how the ancestry clades were phylogenetically defined. In particular, populations predominantly located in South Asia underwent considerably different adaptation as compared with populations from the other geographical regions. Signatures of positive selection present in multiple geographical regions were predicted to be older and have emerged prior to the separation of the populations in the different regions. In contrast, selection signals present in a single population group tended to be of lower frequencies and thus can be attributed to recent evolutionary events.

Carriers of human mitochondrial DNA macrohaplogroup M colonized India from southeastern Asia

Background

From a mtDNA dominant perspective, the exit from Africa of modern humans to colonize Eurasia occurred once, around 60 kya, following a southern coastal route across Arabia and India to reach Australia short after. These pioneers carried with them the currently dominant Eurasian lineages M and N. Based also on mtDNA phylogenetic and phylogeographic grounds, some authors have proposed the coeval existence of a northern route across the Levant that brought mtDNA macrohaplogroup N to Australia. To contrast both hypothesis, here we reanalyzed the phylogeography and respective ages of mtDNA haplogroups belonging to macrohaplogroup M in different regions of Eurasia and Australasia.

Results

The macrohaplogroup M has a historical implantation in West Eurasia, including the Arabian Peninsula. Founder ages of M lineages in India are significantly younger than those in East Asia, Southeast Asia and Near Oceania. Moreover, there is a significant positive correlation between the age of the M haplogroups and its longitudinal geographical distribution. These results point to a colonization of the Indian subcontinent by modern humans carrying M lineages from the east instead the west side.

Conclusions

The existence of a northern route, previously proposed for the mtDNA macrohaplogroup N, is confirmed here for the macrohaplogroup M. Both mtDNA macrolineages seem to have differentiated in South East Asia from ancestral L3 lineages. Taking this genetic evidence and those reported by other disciplines we have constructed a new and more conciliatory model to explain the history of modern humans out of Africa.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0816-8) contains supplementary material, which is available to authorized users.

Molecular Genealogy of a Mongol Queen's Family and Her Possible Kinship with Genghis Khan

Members of the Mongol imperial family (designated the Golden family) are buried in a secret necropolis; therefore, none of their burial grounds have been found. In 2004, we first discovered 5 graves belonging to the Golden family in Tavan Tolgoi, Eastern Mongolia. To define the genealogy of the 5 bodies and the kinship among them, SNP and/or STR profiles of mitochondria, autosomes, and Y chromosomes were analyzed. Four of the 5 bodies were determined to carry the mitochondrial DNA haplogroup D4, while the fifth carried haplogroup CZ, indicating that this individual had no kinship with the others. Meanwhile, Y-SNP and Y-STR profiles indicate that the males examined belonged to the R1b-M343 haplogroup. Thus, their East Asian D4 or CZ matrilineal and West Eurasian R1b-M343 patrilineal origins reveal genealogical admixture between Caucasoid and Mongoloid ethnic groups, despite a Mongoloid physical appearance. In addition, Y chromosomal and autosomal STR profiles revealed that the four D4-carrying bodies bore the relationship of either mother and three sons or four full siblings with almost the same probability. Moreover, the geographical distribution of R1b-M343-carrying modern-day individuals demonstrates that descendants of Tavan Tolgoi bodies today live mainly in Western Eurasia, with a high frequency in the territories of the past Mongol khanates. Here, we propose that Genghis Khan and his family carried Y-haplogroup R1b-M343, which is prevalent in West Eurasia, rather than the Y-haplogroup C3c-M48, which is prevalent in Asia and which is widely accepted to be present in the family members of Genghis Khan. Additionally, Tavan Tolgoi bodies may have been the product of marriages between the lineage of Genghis Khan’s Borjigin clan and the lineage of either the Ongud or Hongirad clans, indicating that these individuals were members of Genghis Khan’s immediate family or his close relatives.