The uplift of the Qinghai–Tibet Plateau and glacial oscillations triggered the diversification of Tetraogallus (Galliformes, Phasianidae)

Abstract The Qinghai–Tibet Plateau (QTP) plays an important role in avian diversification. To reveal the relationship between the QTP uplift and avian diversification since the Late Cenozoic, here, we analyzed the phylogenetic relationship and biogeographical pattern of the genus Tetraogallus (Galliformes, Phasianidae) and the probable factors of speciation in the period of the QTP uplift inferred from concatenated data of four nuclear and five mitochondrial genes using the method of the Bayesian inference. Phylogenetic analysis indicated that T. himalayensis had a close relationship with T. altaicus and conflicted with the previous taxonomy of dark‐bellied and white‐bellied groups. The molecular clock showed that the speciation of Tetraogallus was profoundly affected by the uplift of the QTP and glacial oscillations. Biogeographic analysis suggested that the extant snowcocks originated from the QTP, and the QTP uplift and glacial oscillations triggered the diversification of Tetraogallus ancestor. Specifically, the uplift of the mountain provided a prerequisite for the colonization of snowcocks Tetraogallus as a result of the collision between the Indian and the Arab plates and the Eurasian plate, in which ecological isolation (the glacial and interglacial periods alternate) and geographical barrier had accelerated the Tetraogallus diversification process. Interestingly, we discovered hybrids between T. tibetanus and T. himalayensis for the first time and suggested that T. tibetanus and T. himalayensis hybridized after a second contact during the glacial period. Here, we proposed that the hybrid offspring was the ancestor of the T. altaicus. In conclusion, the uplift of QTP and glacial oscillations triggered the snowcocks colonization, and then, isolation and introgression hybridization promoted diversification.


| INTRODUC TI ON
The uplift of Qinghai-Tibet Plateau (QTP) profoundly impacts speciation and diversification in the plateau birds (Lei, Qu, & Song, 2014). The uplift of the QTP results from the collision of the Indian Plate with Eurasia during the Eocene (55-50 million years ago [Mya]) and then has experienced different stages of growth to reach the current altitude, which has created the unique topography, complex climate, and diversified habitats of the QTP, and makes the QTP an area of worldwide importance for biodiversity Favre et al., 2015;Li & Fang, 1999). Moreover, the unique geomorphological configuration provides multitudinous refugia for animals to survive in harsh environments during the series of major ice ages of the Quaternary (2.4 Mya to the present) (Hewitt, 2000;Lei et al., 2014). Most of the currently endemic birds of the QTP have went through speciation since the Late Cenozoic in response to the QTP uplift, such as partridge Perdix (Bao et al., 2010), monal-partridge Tetraophasis , and ground tit Pseudopodoces humilis (James et al., 2003).
Therefore, the study on Tetraogallus diversification in the QTP can help reveal the geographic history of the QTP and the speciation mechanism. Snowcocks inhabit and evolve in the high-altitude mountain land (Liu, 1998;Shen & Wang, 1963), and make it an ideal animal to investigate the relationship between geological history and speciation. However, the diversification of Tetraogallus has long been ignored and its biogeographical history is poorly understood so far.
T. himalayensis inhabit the Tien Shan, Pamir Plateau, and parts of the QTP, and the distribution area overlaps with T. tibetanus on the QTP (Liu, 1998;Zheng et al., 1978). T. altaicus are distributed in the Altai-Sayan-Hangay Mountains (Liu, 1998;Potapov, 1987;Zabelin, 2007), including the Altai Mountains in China (Huang, Mi, & Shao, 1992). T. caspius are distributed in the southern part of the Caspian Sea from the eastern Anatolian Plateau to the Iranian Plateau, and T. caucasicus only survive in the Caucasus Mountains (Dementiev & Gladkov, 1967;Liu, 1998). Bianki (1898) once thought that extant five snowcocks can classify into dark-bellied (T. himalayensis, T. caspius, and T. caucasicus) and white-bellied (T. tibetanus and T. altaicus) groups based mainly on adult abdomen plumage color. Liu (1998) suggested that T. tibetanus and T. himalayensis first split from the Tetraogallus ancestor based on their morphological characteristics and represent white-bellied and white-bellied groups. T. altaicus split from T. tibetanus; both T. caspius and T. caucasicus split from T. himalayensis and speculated that T. altaicus is likely to hybridize with T. himalayensis (Liu, 1998). Liu (1998) systematically elaborated the phylogenetic relationship of Tetraogallus, but the aforementioned conclusion has not been validated by molecular phylogenetic approach up to now. Stein, Brown, and Mooers (2015) reconstructed the phylogenetic relationship of Tetraogallus except for T. caucasicus using molecular data and the result indicated that T. altaicus has a close relationship with T. caspius; however, this result did not explain the morphological characteristics and the distribution pattern of extant snowcocks. Thus, a further phylogenetic study on the genus Tetraogallus remains critical to understand the speciation in relation to the geographic history.
The origin of Tetraogallus has been controversial over the past few decades (Liu, 1998). In the aspect of origin time of the Tetraogallus ancestor, Koslova (1952) and Baziev (1978) suggested that snowcocks originate in the Early Pleistocene (about 2.5 Mya) based on the orogeneses of the QTP. Potapov (1992) disagreed with the aforementioned point and proposed that the ancestor of Tetraogallus occurs in the first glacier of the Pleistocene (Minde glacial period) because the climate of the Pliocene/Pleistocene boundary is warm and does not have the conditions for the origin of the Tetraogallus, but that is too late for snowcocks to have evolved obviously. As such, Liu (1998) pointed out that the snowcock first occurs in the Quaternary Early glaciation (Hongya glaciation, about 3.5-2.6 Mya). However, recent molecular phylogenetic studies have greatly advanced the origin of snowcocks. Stein et al. (2015) suggested that the ancestor of Tetraogallus can be dated back to the Middle Oligocene (about 29.32 Mya) and diversification begins in the Late Miocene (about 7.31 Mya) inferred from the concatenation dataset of nuclear and mitochondrial genes. In the aspect of origin area of the Tetraogallus ancestor, Koslova (1952) pointed out that the origin of the snowcock-like ancestors is related to the uplift of the QTP, and they may be originated from the jungles of the Eastern Kunlun Mountains and the Hengduan Mountains (Western Sichuan Province, China). Baziev (1978) believed that the snowcock-like ancestor is not alpine birds, and they inhabit the hilly areas from the Caucasus to Central China; the snowcock-like ancestor will rise with the process of mountain uplift, so each mountain system has a kind of snowcock.
Consequently, the distribution area of extant snowcocks is the origin of ancestor. However, Potapov (1992) found a paradox that snowcocks are unlikely to inhabit the Hengduan Mountains during the Pliocene/Pleistocene boundary because of lower altitude and suggest that snowcocks originated in the Pamir Plateau, Tien Shan, and Kunlun mountains. Fossil data can be used to estimate the origin of snowcocks, in addition, however snowcocks have few fossils available, and precious fossil data are only found in the Altai and Caucasus Mountains, which can date back to the Middle and Late Pleistocene, and these are not sufficient to infer the origin of Tetraogallus because time is too short (Panteleev, 2002;Potapov, 1992). In summary, the origin of the Tetraogallus ancestor has not yet been determined in terms of spatial and temporal.
To gain knowledge on the speciation mechanisms in snowcocks, here, we investigated the phylogeny, biogeographical history, and diversification rate using a concatenated DNA dataset.
As a consequence, the present investigation aimed at solving three questions: 1. When and where did snowcocks originate? 2. Was the diversification of snowcocks affected by the QTP uplift and glacial oscillations? 3. Was hybridization promoting the speciation of T. altaicus?
Here, T. caucasicus was not taken into account in the ingroup taxon set because its DNA data were not available up to now. All DNA sequence data were downloaded from the GenBank database (downloaded on or before December 30, 2018; Table A1) in the present study. Two outgroup taxa were selected from the genus Alectoris based on the previous study (Stein et al., 2015), including A. chukar and A. rufa. We assembled a DNA data matrix which was composed of 3 protein-coding genes (COX1, CYTB, and ND2) and 2 nonprotein-coding genes (12S and D-loop) mitochondrial loci combined with 4 nonprotein-coding genes (CLTC, CLTCL1, RHO, and EEF2) nuclear loci (Table A1). The DNA matrix was 6,975 base pairs (bp) long at its maximum extent (including gaps), and the alignment data were sparse and average coverage was approximately 85% across DNA markers (mitochondrial locus: 67%-100%; nuclear locus: 83%).
In order to improve the resolution ability of the phylogenetic inference, here, we chose to mark as high coverage as possible rather than adopted all publicly available nucleotide sequences. In addition, D-loop sequences were downloaded from the GenBank database and were used to construct phylogenetic analysis of populations between T. tibetanus and T. himalayensis. Sequence information please see Table A2.

| Phylogenetic analysis, genetic distance, nucleotide mutation rate, and estimation of divergence time
Multiple-sequence alignments of mitochondrial and nuclear genes were performed using BioEdit v. 7.1.3.0 program (Hall, 1999) with default parameters. Each gene was aligned separately and manually concatenated these sequences using Sequence Matrix v. 1.7.8 (Vaidya, Lohman, & Meier, 2011). All concatenated sequences were format converted for further analyses using Geneious v. 9.1.4 program. The genetic distance was calculated by MEGA 6.0 (Tamura, Stecher, Peterson, Filipski, & Kumar, 2013) with the Kimura 2-parameter (K2P) model of nucleotide substitution (Kimura, 1980) using concatenated data. Phylogenetic analysis and divergence date on the genus Tetraogallus species were estimated by BEAST version 1.7.4 program (Drummond & Rambaut, 2007) using the Markov Chain Monte Carlo (MCMC) to estimate the divergence times of the genus Tetraogallus. The best-fit model of substitution and best partition schemes for the dataset were identified with the corrected Akaike information criterion (AICc) (Akaike, 1974), implemented in Modeltest 3.7 program (Posada & Crandall, 1998).
To acquire as accurate divergence time as possible, we estimated the number of nucleotide substitutions per site (d) from comparisons of the focal species and an outgroup species using the formula: d = (tv + tvR)/m, where tv is the number of transversions between the genus Tetraogallus and outgroup taxa, R is the transition/transversion F I G U R E 1 Geographical distribution of Tetraogallus species. Species distribution information is cited from Liu (1998) ratio within the genus Tetraogallus, and m is the sequence length (Nei, 1992;Rooney, Honeycutt, & Derr, 2001). Transition and transversion values were calculated in the program MEGA 6.0. The rate of nucleotide substitutions per site per lineage per year is λ = d/2T when an estimate of d was obtained, where T is the divergence time between the ingroup and outgroup species (Rooney et al., 2001).
The mutation rate per nucleotide site per generation is μ = λg, where g is the generation time (g = 3 years of snowcock; Huang, Ma, Shao, & Jiang, 1990). Here, d was 0.084 (tv = 84, R = 5.96, m = 6,975) for combined genes. The rate of nucleotide substitution per site per lineage per year (λ) was about 0.14 × 10 -8 (T = 29.32 Mya, quoted from Stein et al. (2015)), and the mutation rate of per generation (μ) was about 0.40% Mya for concatenation DNA sequence of snowcocks.
Similarly, we used the same method to calculate the rate of nucle-

| Biogeographical analysis
To infer the possible origin area of the genus Tetraogallus, we used the program with the method of BBM (Bayesian binary MCMC) in RASP (Reconstruct Ancestral State in Phylogenies) (Yu, Harris, Blair, & He, 2015) to reconstruct the possible ancestral region of snowcocks on the phylogenetic tree. Here, the biogeographic distribution range of extant species of Tetraogallus was divided into four sections according to the available literature (Liu, 1998;Potapov, 1987;Zheng et al., 1978)  was used for BBM analysis with null root distribution. The possible ancestral ranges were determined at each node on a selected tree.

| Diversification analysis
To reconstruct the macroevolutionary dynamics of Tetraogallus over time, we implemented a diversification rate analysis. BAMM software used a Bayesian framework to account for rate variation through time and among lineages in phylogenetic trees, and we employed it to simulate posterior distributions of rate-shift configurations (Rabosky, Donnellan, Grundler, & Lovette, 2014a). We then used BAMMtools package in R to conduct rate-through-time analysis and to identify and visualize diversification rate shift along branches (the BEAST output result) (R Core Team, 2012; Rabosky, Grundler, et al., 2014b). BAMM was run by setting four independent MCMC running for 10 million generations and sampled every 1,000 generations. After removing 25% of trees as burn-in, the BAMM output was analyzed with employing BAMMtools and the 95% credible rate-shift configurations were estimated using Bayes factors. The best shift configuration with the highest maximum a posteriori probability was estimated in this analysis. In addition, we demonstrated a relationship between diversification rate shifts and ecological opportunity processes, which if ecological opportunity occurred in new regions could trigger rapid radiation and cause speciation rate increase and colonization of a new region (Schenk, Rowe, & Steppan, 2013;Yoder et al., 2010). The phytools package in R was used to estimate a gamma statistic (γ) with a phylogenetic tree to test region-specific patterns in the diversification rates and to visualize lineage-through-time (LTT) plot (Revell, 2012). Herein, if γ value was positive (γ > 0), it suggested a late outbreak of speciation rate, and conversely, a negative value for γ (γ < 0) indicated that this branch has a process of decreasing the speciation rate (Pybus & Harvey, 2000).

| Phylogenetic analysis, divergence time, and genetic distance
Phylogenetic relationships of the genus Tetraogallus were reconstructed using the concatenated DNA data with the method of Bayesian inference. The results indicated that T. tibetanus and T. himalayensis were first split from the Tetraogallus ancestor (PP = 1) and dated back to the Late Miocene (about 5.91 Mya) with 95% confidence intervals (95% CI) of 3.09-8.87 Mya (Figure 2 Figure 2).

| Biogeographical analysis
BBM analysis indicated that the ancestor of the genus Tetraogallus originated from the area A at node 7 with high marginal probability (94.4%) (Figure 3a). Similarly, the ancestral reconstruction of BBM suggested that A was the ancestral area of nodes 6 and 5, with 50.1% and 55.1% marginal probability, respectively. BBM results suggested that the speciation between T. tibetanus and T. himalayensis occurred in the QTP during the Late Miocene (node 7), and then they underwent dispersal and isolation that may be driven by the geographic and climatic events, as evident from ancestral ranges at nodes 6 and 5. BBM detected two vicariance events and three dispersal events at nodes 6 and 5, while node 7 without any biogeographic event occurred. Moreover, BBM estimated the possible colonization route of snowcocks (nodes 7, 6, and 5) as the result that the origin and spread of snowcocks were related to the QTP uplift. In addition, the time curve of dispersal and vicariance events illustrated that a fastigium of biogeographical events occurred in the Early Pleistocene (about 2.28 Mya) (Figure 3b). Herein, it was clear that this discontinuous distribution pattern was related to geological and climatic changes.

| Diversification rate shift
The distinct shift configuration for diversification was shown in Figure 4a, along with the corresponding consensus phylorate plot and rate acceleration events within the genus Tetraogallus taxa which fell into the Early Pleistocene (about 2.28 Mya). Here, the speciation rate of snowcocks has undergone a rapid increase and then subsequent decrease, and we detected the best rate-shift point which was noted in the phylorate plot (Figure 4a). The LTT plot showed an increase of diversification rate after colonized the F I G U R E 2 A phylochronogy of the genus Tetraogallus species based on concatenated nucleotide data. The divergence time is estimated using the BEAST with the calibration method under the relaxed molecular clock model (using the estimated mutation rate). Branch lengths represent the mean values of the posterior distribution. The posterior probability and divergence time are indicated at each inner node. The node bars indicate the posterior probability distribution of the node age under the 95% CI. Snowcock portraits are quoted from MacKinnon, Phillipps, He, and Lu (2000) and Svensson, Mullarney, Zetterstrom, and Grant (2009) (Figure 4b) and suggested a potential early initial burst at approximately 6 Mya, and after that the speciation rate has undergone a subsequent reduction (γ = −0.082, p > .05) in line with phylorate ( Figure 4a).

| Phylogeographic analysis between T. tibetanus and T. himalayensis
Phylochronology showed two distinct lineages (PP = 1) and the popula-  1.04-2.81 Mya). Here, we speculated that the hybridization area was likely to be located on the sympatric zone of T. tibetanus and T. himalayensis ( Figure 5; Top-left: dim gray shade).  which predated previous studies (Baziev, 1978;Koslova, 1952;Liu, 1998;Potapov, 1992) and can adequately reveal the speciation mechanism of snowcocks because for this timeframe takes about 2 million years for it a vertebrate to evolve into a new species (Avise, Walker, & Johns, 1998 BBM analysis indicated that the Tetraogallus ancestor originated from the QTP (including Pamir Plateau) (Figure 3). Previous studies on the origin area of snowcocks supported our result (Koslova, 1952;Liu, 1998;Potapov, 1992) and suggested that the uplift of the QTP triggered the Tetraogallus diversification. The uplift of the QTP is a complicated geological development process and has underwent several uplifts which is the progressive and heterogeneous uplift of the QTP from south to north (Favre et al., 2015;Mulch & Chamberlain, 2006;Wang et al., 2008).

| D ISCUSS I ON
Subsequently, a series of rapid uplift events give rise to approach its present elevation during the Late Cenozoic (Fang, 2017;Li & Fang, 1999;Shi, Li, & Li, 1998) (Li & Fang, 1999). The progressive extension of the uplift of the QTP is associated with the major mountains uplift and has blocked the northward of the India Ocean warm air, such as the higher Himalayas and the rise of the Tien Shan, which causes the progressive aridification of Central Asia during the Miocene (Miao, Herrmann, Wu, Yan, & Yang, 2012;Miao et al., 2011;Sun, Gong, Tian, Jia, & Windley, 2015). Moreover, the development of the Pleistocene glaciers of the QTP is closely related to the progressive uplift of the plateau and the surrounding mountains (Zheng & Rutter, 1998). Consequently, the QTP after the uplift results in drastic shifts in the distribution of plant communities and major faunal turnover (Deng & Ding, 2015;Li, Fang, Pan, Zhao, & Song, 2001;Sun & Wang, 2005). All these factors can be beneficial snowcocks to acquire suitable ecological resource and promoted adaptive radiation, of which the geographic colonization is an important way to survive (Stroud & Losos, 2016). The colonization of novel territory (e.g., mountains uplift or glacier recession) can provide a release from competition and predation pressures, abundant food resources, or suitable climatic conditions, allowing them to differentiate a variety of species to colonize multiple unexploited ecological resources, such as ground tit (James et al., 2003) and voles (Lv, Xia, Ge, Wu, & Yang, 2016). In summary, this evolutionary model can be characterized by ecological opportunity. Ecological opportunity can affect speciation rate, and animals will soon fill unoccupied niche if rate increase and then decrease tend to be stable finally (Rabosky & Lovette, 2008;Schenk et al., 2013). These theories were in favor of the diversification mechanism of snowcocks. Here, the diversification analysis of snowcocks indicated that the increase of speciation rate can be dated back to the Late Miocene and peaked in the Early Pleistocene (about 2.28 Mya, Figure 4), which coincided with the active period of the Late Cenozoic of the QTP. Obviously, the phases of A and B profoundly influenced on the snowcocks speciation and supported the snowcocks speciation was subjected to ecological opportunity. As such, the LTT analysis showed that the snowcocks speciation was subjected to ecological opportunity even though the gamma value was not significant because of the small number of species. In conclusion, the uplift of the QTP provided ecological opportunity for Tetraogallus diversification.
The T. tibetanus, an endemic to the QTP, has the long evolutionary history among the extant snowcocks. T. tibetanus split from T.
himalayensis that can date back to the Late Miocene to coincide with the geological active period of QTP and Asian interior aridification , namely that the uplift of the QTP promoted speciation of T. tibetanus. Here, similar geological event occurred on the Phrynocephalus theobaldi differentiation (about 5.65 Mya), and the plateau uplift promoted the differentiation of toad-headed lizards that inhabited the QTP during the Late Miocene (Jin, Liu, & Brown, 2017), and maybe, similar geographic events induced the divergent possibility between T. tibetanus and T. himalayensis at that time.
Theoretically, the Gloger's rule believes that melanin is abundant in hot and humid areas, while maroon pigment or tawny pigment is abundant in dry areas (Edward, Burtt, & Jann, 2004;Zheng, 1952).
We speculated that the origin area of T. tibetanus was likely to occur in the Gangdise Mountains and the Himalayas with wetter air and warm climate before the intense uplift of the QTP during the Late Cenozoic (since 3.6 Mya, Li and Fang (1999)) based on black stripes of belly feathers. Moreover, the origin area of T. himalayensis was likely to take place in the Pamir Plateau or the area east of the Pamir Plateau (Western Kunlun Mountains) with drier air and cold climate.
A rational explain was that the Gangdise and Himalayas are the first to uplift as the result that has caused ecological shift in different regions and can provide ecological opportunity for T. tibetanus and T. himalayensis to colonize because the geographical development of the QTP is progressive uplift from south to north (Li & Fang, 1999;Mulch & Chamberlain, 2006). Competitive exclusion and ecological F I G U R E 6 Proposed speciation model of T. altaicus. The hybridization between T. tibetanus and T. himalayensis occurs in the sympatric zone during the Quaternary Early glacial period, and then reproduces F1 generation. The F1 hybridizes with T. himalayensis and reproduces F2 during the interglacial period, and then several generations of introgression allowed hybrids to inherit numerous loci from T. himalayensis, which can be combined via gene flow through hybridization. Competitive exclusion between them gives rise to the introgression zone shift along with the Tien Shan and makes Fn (hybrids) to be an ancestor for T. altaicus. Hybrids (Fn generations) can potentially be adaptive and favored in a new habitat via adaptive introgression and can lead to a separate hybrid taxon (hybrid speciation). Finally, competitive exclusion and ecological opportunity can explain that hybrid species (the ancestor of T. altaicus) colonize the Altai-Sayan-Hangay Mountains during the glacial period, and then, geographical isolation and adaptive radiation have driven hybrids to be a novel species-T. altaicus F I G U R E 5 A chronogram of population divergence between T. tibetanus and T. himalayensis based on D-loop haplotypes. The divergence time is estimated using the BEAST with the calibration method under the strict clock model (using the estimated mutation rate). Branch lengths represent the mean values of the posterior distribution. The posterior probability and divergence time are indicated at the major inner nodes, and population names are provided in Table A2. The node bars indicate the posterior probability distribution of the node age under the 95% CI. Proposed hybridization zone of T. tibetanus and T. himalayensis is marked in dim gray shade (Top-left) isolation caused the differentiation between T. tibetanus and T. himalayensis (Liu, 1998(Liu, , 1999. Ecological opportunity had created a novel habitat for T. tibetanus and T. himalayensis, and then adaptive radiation accelerated them to divergence. Consequently, ecological isolation drove two snowcocks (T. tibetanus and T. himalayensis) to occupy different niches to evolve in their respective directions (Liu, 1999), including altitude, diet, breeding strategy, and breeding time, and finally evolved into two separate species. Geographical or ecological isolation leads to a reduction or disruption of genetic exchange between populations, and alleles change under the influence of genetic drift, which ultimately give rise to reproductive isolation (Hoskin, Higgie, McDonald, & Moritz, 2005;Via, 2012). For example, the uplift of the QTP gives rise to the speciation of the genus Perdix (Bao et al., 2010), and ecological divergent is an important factor for eared pheasants Crossoptilon speciation (Wang et al., 2017). In summary, the uplift of the QTP and niche divergence promoted the speciation of T. tibetanus and T. himalayensis. Liu (1998) (Berberian, 2014;Mouthereau, 2011).
After that, the tectonic movements of Iranian Plateau give rise to the climatic change and the formation extensions of steppe-like landscape during the Middle and Late Miocene (Storch, 2004), in which provides several suitable habitats for snowcock survival and also provides a prerequisite for T. himalayensis colonization. Furthermore, the Northern hemisphere great glacial period has led to the widespread development of ice sheets or glaciers in low-latitude areas Curry, 1966;Wu & Li, 1990) and promoted the T. CI: 1.04-2.81 Mya). When allopatric populations meet in the contact zone, theoretically, hybridization may occur if reproductive isolation is incomplete (Dowling & Secor, 1997). When allopatric taxa became sympatric again, moreover, speciation can accelerate because of secondary contact such as the genus Drosophila (Coyne & Orr, 1997). For instance, the Darwin's Heath (Coenonympha darwiniana) originate through hybridization between the Pearly Heath (C. arcania) and the Alpine Heath (C. gardetta) with different parental contributions (Capblancq, Després, Rioux, & Mavárez, 2015), and hybridization promotes speciation in Coenonympha butterflies.
A possible explanation is that the Quaternary Early ice age (e.g., Danube-Gonzi ice age (about 2.6-1.5 Mya), Hongya ice age (about 2.5 Mya) (Liu, 1998), and Xixiabangma ice age (about 2.5-0.78 Mya) (Zheng, Xu, & Shen, 2002)) caused T. tibetanus and T. himalayensis to secondary contact and hybridize in sympatric zone, and then, the introgressive hybridization between F1 and T. himalayensis occurred in the interglacial period ( Figure 6) (Chen, An, & Liu, 2016;Liu, Wen, Huang, & Hou, 2006), and supported our inference. The niche overlap forced in- Australian parrots. Grant, Grant, and Deutsch (1996) supported that hybridization can contribute to the speciation process by enhancing genetic variation and relaxing genetic constrains on particular directions of evolutionary change, such as some island birds (Darwin's Finches). In conclusion, the genetic exchange between the backcross offspring and T. himalayensis was weakened as a result that the mountains uplift promoted the ancestor of T. altaicus to colonize and Middle-Asia aridification caused geographical barriers, and the interaction of these factors promoted a new species differentiate-T.

| CON CLUS IONS
The extant snowcocks originated from the QTP, and the uplift of the QTP and the Quaternary glacial oscillations accelerated the genus himalayensis during the interglacial period as a result of inheriting many characteristics from T. himalayensis, and glacial dispersal and isolation finally promoted the speciation of T. altaicus.

ACK N OWLED G M ENTS
The work was supported by the National Natural Science Foundation of China (Nos. 30530130). Li Ding and Jicheng Liao dedicate this paper to the late Mr. Naifa Liu.

CO N FLI C T O F I NTE R E S T
The authors declare that they have no conflict of interests.

AUTH O R CO NTR I B UTI O N S
NF L conceived the study; L D was responsible for obtaining and analyzing data and wrote the manuscript with the help of JC L.

DATA AVA I L A B I L I T Y S TAT E M E N T
All data used in this review paper have been published elsewhere.
Genbank accession number please see the Tables A1 and A2.