Genomic analysis of demographic history and ecological niche modeling in the endangered Chinese Grouse Tetrastes sewerzowi
BMC Genomics volume 21, Article number: 581 (2020)
The Quaternary had worldwide consequences in forming the contemporary diversity of many populations, species and communities, which is characterized by marked climatic oscillations between glacial and interglacial periods. The origin and evolution of biodiversity in mountainous areas are highly dependent on historical orogenesis and associated climatic changes. The Chinese grouse Tetrastes sewerzowi is a forest-dwelling species endemic to the mountains to the east of the Qinghai–Tibet Plateau, which has been listed as Near Threatened with a decreasing trend by the IUCN because of ongoing deforestation and fragmentation of coniferous forests. It is important to place current population status into a broader ecological and evolutionary context to understand their demographic history.
Analyses of the Chinese Grouse genome revealed fluctuations throughout the Pleistocene in effective population size. Populations decreased during early to middle Pleistocene but showed an expansion during late Pleistocene which was then followed by a sharp decline during the last glacial maximum (LGM). Ecological niche modeling indicated that a suitable habitat shift between high altitude regions to low altitude regions was due to a changing climate. This result parallels patterns of population size change in Chinese Grouse estimated from PSMC modelling, which suggested an expansion in population size from the last interglacial period (LIG) and then a peak and a bottleneck occurring at the last glacial maximum (LGM). Furthermore, the present-day distribution of Chinese Grouse is greatly reduced and fragmented. It will likely become even more fragmented in the future since coniferous forest cover is threatened in the region of their distribution and the availability of such habitat restricts their ecological niche.
The Chinese Grouse have experienced substantial population size changes from the beginning to the LIG and reached a peak before the LGM. A sharp decrease and bottleneck occurred during the LGM, when the coniferous forests were subjected to extensive loss. The results inferred from the whole genome sequencing and species distribution models both support historical population fluctuations. The distribution of the Chinese Grouse is strongly dependent on the coniferous forest cover. To protect the fragmented coniferous forests is an essential action to protect the Chinese Grouse.
The Chinese grouse is a forest-dwelling bird species endemic to the Qinghai–Tibetan Plateau (QTP) [1, 2], which has been listed as Near Threatened with a decreasing trend by the International Union for the Conservation of Nature (IUCN) because of ongoing deforestation and fragmentation of coniferous forests . Throughout the range of the species, only the wetter northern slopes of mountain coniferous forests have vegetation that support grouse populations . This results in ongoing significant reductions of population size and area of occupancy when such areas are cut . In addition, there are large-scale deforestation and intensive livestock grazing in QTP area, which exacerbates habitat loss and fragmentation . Furthermore, world climate change possibly threatens the Chinese Grouse not only by increasing habitat loss and fragmentation, but also by influencing range shifts [5, 6].
In mountainous areas, the origin and evolution of biodiversity are highly dependent on the historical orogenesis and past associated climatic changes [7,8,9]. Complex geological environments and Pleistocene glacial events contributes to complicated refugial patterns, species distributions and the demographic history of populations [9,10,11,12,13]. Four significant glaciations have occurred in the QTP during the Pleistocene, the Xixiabangma (Early Pleistocene), Nyanyaxungla (Middle Pleistocene; 0.78–0.50 million years ago), Guxiang (late Middle Pleistocene) and the last glaciation, including two glacial stages with the last corresponding to the LGM [10, 14, 15]. The coniferous forests declines and expansions coincide with the glaciations. Cold and dry weather conditions happened in the glaciation period, which resulted in decline of coniferous forests. On the contrary, warm and wet weather conditions followed the retreat of the glaciers which were ideal for the spread of coniferous forests. Pollen data and analyses of past animal faunas show that the forest and forest steppe including spruce and fir connected the QTP and the Siberian taiga by the loess plateau in northern China during the Mid Pleistocene . The distribution of existing Palearctic realm taiga birds and mammals, such as Black Woodpecker (Dryocopus martius), Spotted Nutcracker (Nucifraga caryocatactes), Eurasian Jay (Garrulus glandarius), Snowy-browed Nuthatch (Sitta villosa), Red Deer (Cervus elaphus) support this past connection. A large altitudinal gradient across the region spanning from 500 to 8848 m is created by the uplift of the QTP . The eastern edge of the QTP are now characterized mainly by a warm and wet climate associated with deep valleys and .
The Quaternary had worldwide consequences in forming the contemporary diversity of many populations, species and communities, which is characterized by marked climatic oscillations between glacial and interglacial periods that [19,20,21,22]. The effect of the Pleistocene climate fluctuations on species distributions and diversity is well known in Europe  and North America . The periodic uplift of the QTP before the Pleistocene has contributed to geologically induced complex environments where Pleistocene glaciations formed many relatively small alpine glaciers and refugia [10, 25, 26]. As a consequence, more complex glacial effects may have influenced biotic speciation, distribution and diversity in the QTP compared to Europe and North America [10, 11]. Before the LGM, the warm and wet weather during the “Greatest lake period” (25,000–40,000 years ago) in this region could have contributed to the expansion of the coniferous forests [27,28,29]. Compared to today, the temperature was 2–4 °C higher and precipitation was 40–100% higher . Furthermore, pollen-based vegetation reconstruction shows that at this time [28, 30] the alpine coniferous forest extended beyond its present western limit about 400–800 km. This expansion of the alpine coniferous forests, the primary habitat for the Giant Panda (Ailuropoda melanoleuca) as well as the Chinese Grouse, have resulted in a second population expansion of Giant Pandas [1, 12]. The Chinese Grouse as well as the other typical forests-dwelling birds and mammals should have had similar demographic histories during this period.
The aim of the present study is to infer the past population dynamics of Chinese Grouse by applying Partial Sequential Markov Chain (PSMC) modelling on whole genome sequence data. The PSMC model infer the time of the most recent common ancestor for each independent DNA segment and infer ancestral effective population size (Ne) at a given epoch through a hidden Markov framework and identifies historical recombination events across a single diploid genome . Combined with glacial event data, we aim at understanding the effects of past climatic changes on past distribution changes of Chinese Grouse via Species Distribution Modelling (SDM). The combination of SDM and population genomics also provides new insights to understand the impact of past climatic changes on population dynamics. At present, human interference activities had significantly changed climate and habitats in the QTP of southwest China, which have already had, and will continue to cause, untold biodiversity loss . We also performed a GAP analysis, a Geographic Approach to Protection of Biological diversity , between present forests and climate driven habitat to inspire future work in the conservation of Chinese Grouse and other forest species in the QTP.
De novo assembly of the Chinese grouse genome
The Illumina sequencing had total 184.34 G PE 150 bp reads and the average depth was 184 X. The Pacbio long-read Sequencing had total length of 36.7 GB and the depth was 36 X. The software DBG2OLC produced a reference de novo genome assembly which comprises 1272 scaffolds that had a total length of 1,003,232,344 bp, corresponding to 91% of the total genome length of the expected 1.1 G. Scaffold N50 was 2.16 Mb, and the longest scaffold measured 5.27 Mb. The BUSCO complete assessment score is 93%. For comparison, the well assembled Chicken (G. gallus) genome is 1054.6 Mb , which is similar with the genome size of Chinese Grouse in the current project.
The simulated demographic history of Chinese Grouse was inferred by plotting the population size fluctuation curve against time with superimposed glacial periods . Chinese Grouse effective population size increased from the time of the last interglacial (around 115 kya ago). Before this time, the effective population size of Chinese Grouse was kept at a high number (20 × 104) from the beginning of the simulations. The results indicated that Chinese Grouse population size reached a peak before the last glaciation Maximum (40–100 kya ago), experienced a significant decrease during the last glaciation Maximum (LGM; ~ 20 kya ago) (Fig. 1), when substantial alpine glaciations (such as Gongga glacial II) would likely have resulted in extensive loss of coniferous forests.
Ecological niche modeling
To assess the role of climate change, we used SDMs and palaeoclimatic data on temperature and precipitation to estimate the suitable habitat of Chinese Grouse during the time periods of 120–140 kya (LIG), 21 kya (LGM), 6 kya (Mid-Holocene), and the present (Fig. 2, Additional file 2). All four SDMs suggest that past climate change may have contributed significantly to the distribution change and population dynamics of Chinese Grouse. The model revealed that suitable distribution during the colder LGM was more widely available than during the warmer LIG and warmer mid Holocene (Fig. 2). Specifically, compared with the LIG period, the area of suitable distribution during the colder LGM period had increased about 11. 9375 thousand km2 and decreased 2.3875 thousand km2 (Fig. 2a, Fig. 3). The expansion of suitable distribution was moved eastwards which was a low elevation area. On the contrary, the area of suitable distribution during LGM period was less compared with the Mid Holocene period (Fig. 2b). The suitable distribution during Mid Holocene moved to westwards which was at high elevation in the Tibetan plateau and the former eastern area was lost (Fig. 2b).
The population size change in Chinese Grouse estimated from the PSMC model suggested an expansion in population size during the LIG and a bottleneck occurring at the LGM. Strikingly, the suitable distribution of Chinese Grouse expanded toward warmer eastern areas during the LGM and moved to higher and western areas during Mid Holocene. Compared with the Mid Holocene period, the predicted present-day distributions is increased and concentrated predominantly in the model only conducted in climate change scenario (Fig. 2c, Fig. 3). But the present-day distribution of Chinese Grouse was greatly reduced and become highly fragmented after the GAP analysis was performed between the predicted present-day distributions and coniferous forest cover (Fig. 4).
The demographic analyses showed the grouse population reached a peak ~ 40,000 years ago followed by sharp decrease with a bottleneck ~ 20,000 years ago. The population expansion can be explained by the warmer weather during the Greatest Lake Period (30,000–40,000 years ago) as the conifer forests, the primary habitat for Chinese Grouse , reached their greatest extent at this time [10, 12]. The Palaeo-distribution of coniferous forests changes in the QTP and adjacent mountains inferred by pollen cores  support the demographic changes inferred of Chinese Grouse [35, 36]. The demographic results and palaeo coniferous forests data show that coniferous forests played a pivotal role in Chinese Grouse demographic history by providing permanent refugias in the eastern QTP during the Quaternary. The QTP with its special intricate geographical environment had a significant influence on biodiversity and changes in the geographic distribution of animals in east of Asia during the Pleistocene and may play a key role on potential responses and feedbacks to global change in the future [11, 37].
Our demographic analyses highlighted that a high Ne of Chinese grouse persisted from the Early Pleistocene (2.43 million years ago-0.73 million years ago) to the Mid Pleistocene (0.73 million years ago - 0.10 million years ago) and that a population expansion occurred after the retreat of the Penultimate Glaciation (0.30–0.13 million years ago). Three significant glaciations have occurred in QTP during this time, the Xixiabangma glaciation in Early Pleistocene, Nyanyaxungla glaciation in Middle Pleistocene and Guxiang glaciation in late Middle Pleistocene [14, 15]. Like previous work, our results suggest that the QTP had several Pleistocene refugia which have affected Ne in this grouse inhabiting coniferous forests . During the Pleistocene, many avian species living in coniferous forest in the boreal taiga and the QTP have diverged and are classified as superspecies or different subspecies, as a result of allopatric divergence and local adaptation [8, 23, 38, 39]. A narrow high altitude boreal forests were preserved in the southeast edge of the QTP by the uplift which comprise one of the key high-altitude biodiversity hotspots in the world . This region harbors one of the world’s richest fauna and floras, because it harbored large and complicated refugial areas in different mountain regions, lowland, rivers and basins [10, 41, 42].
When substantial alpine glaciations (for example, Gongga glacial II) likely resulted in extensive loss of coniferous forests, a population bottleneck occurred during the last glacial maximum (~ 20,000 years ago), . Population bottlenecks in the same region coinciding with the LGM have also been identified in other animals such as the giant panda  and snub-nosed monkeys (Rhinopithecus) , suggesting that the Last Glaciation and palaeo coniferous forests changes had strong impacts on the population sizes of arboreal species. The sharp decline in Ne for Chinese Grouse, as in the Giant panda, snub-nosed monkeys, is consistent with the extreme cooling of the climate during the LGM.
Many previous studies have suggested that the LIG was a favorable period for other avian species, such as black grouse (Tetrao tetrix) because of increased seasonality in temperatures [43, 44]. During this period, black grouse and Greenland rock ptarmigan (Lagopus muta) populations expanded to a high Ne and reached their peaks [43, 45]. It is suggested that the severe reduction in avian population sizes approximately coincide with the beginning of the last glacial period (LGP, 110–12 kya) or occurring during this period . In contrast to previous work, the demographical results of Chinese Grouse show no fluctuations associated with historical climatic changes during Early Pleistocene and Middle Pleistocene . Instead a peak followed by a sharp decrease and a bottleneck occurred during this period. This pattern observed in Chinese Grouse may be related to the special geographical position, climate change and the mass accumulation rate (MAR) of Chinese loess, an index indicating cold and dry or warm and wet climatic periods in China [46, 47].
The SDM results showed that climate change had an influence on the amount and distribution shifts of suitable habitat of Chinese grouse. The model also revealed that suitable habitat during the LGM was more widely available than both during the warmer LIG and mid Holocene. Climate events during the LGM period resulted in an expansion of suitable habitat to low altitude regions in eastward areas compared with the LIG period and a shrink from this low altitude region to a more westward distribution during the mid-Holocene period. The result parallels patterns of population size change in Chinese Grouse estimated from the PSMC model, which suggested an expansion in population size from LIG period and then a peak and a bottleneck occurring at the LGM. Significantly, the suitable habitat of Chinese Grouse expanded toward lower altitudes in eastern areas during the LGM. In contrast, the population moved into high altitudes in the western areas during the mid-Holocene. It is of significance that the geographic distribution of Chinese Grouse populations shifted in an east-west direction from LIG to the LGM and the LGM to mid-Holocene. The reason is that the shift of coniferous forests was significantly affected by seasonal extremes in temperature. In addition, the formation of the loess plateau stopped the forest to spread to the north.
Alarmingly, our results suggest that although suitable areas for the Chinese Grouse in the present-day are larger than during the mid-Holocene, the LGM and the LIG, the Ne has been kept at a low size after the bottleneck during the LGM. Dense coniferous forests are of critical importance for the Chinese Grouse, which are affected by climate change and recent human actions . Our GAP analysis between the present-day distribution of Chinese Grouse and extent of present day coniferous forest shows that the Chinese Grouse populations were severely affected by loss and fragmentation of suitable forests. Forest fragmentation has profound and lasting influences not only on the Earth’s biological diversity, but on ecosystem function and numerous ecosystem services [48, 49]. The forest in the QTP started to become protected from 1998 when a disastrous flooding happened. Because of the increase in density of local people with the high demands for farmland, timber, firewood and roads, deforestation of Chinese grouse’s mountain habitats continued [1, 6]. Previous work on Chinese Grouse, blood pheasant (Ithaginis cruentus) and Sichuan jay (Perisoreus internigrans) have found that their distributions would be decreased and fragmented severely in the future [5, 6, 50]. The Chinese Grouse distributions under climate change scenarios is similar to other montane species , which would shift northward and upward under realistic climate change scenarios . This long term fragmentation of coniferous forests may have effect on gene flow, speciation and divergence of these endemic isolated populations.
The Chinese Grouse have experienced substantial population size changes from the beginning to the LIG and reached a peak before the LGM. A sharp decrease and bottleneck happened during the LGM, when the conifer forests were subjected to extensive loss. The results inferred from the whole genome sequencing and species distribution models both support a history of population size fluctuations. The distribution of the Chinese Grouse were strongly dependent on the coniferous forest cover. To protect the fragmented such forests is an essential action to protect the Chinese Grouse.
A blood sample was collected from one individual from the Lianhuashan Nature Reserve, Gansu province, China. The sample collection was conform to the National Wildlife Conservation Law in China and with permission from the Forestry Department of China. Our sample collection procedures also followed the regulations of the animal experimental and medical ethics committee of the Institute of Zoology, Chinese Academy of Sciences. The bird was released after the blood sample collection.
De novo genome sequencing and data production
We used a Gentra Puregene Blood Kit (Qiagen) to extract the genomic DNA from a whole blood sample according to the manufacturer’s instructions. Then we assessed the quality of DNA by electrophoresis on 1% agarose gel and the quantity of DNA by a BioDrop mLITE pectrophotometer (a total of 15 mg of DNA was quantified using the spectrophotometer). The sequencing included two platforms: Illumina (San Diego, CA) Hiseq4000 and Pacific Biosciences (Menlo Park, CA) PacBio RSII. We used Illumina DNA Sample Prep Kit to construct three paired-end libraries with insert size of 350 bp and SMRT bell library preparation protocol (10 SMRT cells were sequenced) to generated by Pacific Biosciences long reads (> 10 kb).
We removed contaminated reads with Illumina paired-end adaptors using TRIMMOMATIC 0.36  with the options ILLUMINACLIP: TRUSEQ3-PE.FA: 2: 30: 10 MINLEN: 150. Then bad reads were eliminated based on thresh hold of: base quality lower than Phred score of 20; N content above 10%; base quality below 20 proportion above 50%. This allowed us to only keep high-quality paired end reads for the assembly steps.
We estimated the best k-mer length for genome de novo assembly by KmerGenie. We used a hybrid assembly strategy following the protocol of DBG2OLC toolset . The parameter: LD 0 k 57 g 25 NodeCovTh 2 EdgeCovTh 1 GS 4000000000 were employed during SparseAssembler step to process the illumina paired-end reads. K 57 was chosed according to k-mer estimation of KmerGenie. GS was set to about 4 folds of conventional avian genome size 1G for memory usage management. The preliminary contigs file was passed to DBG2OLC to make layout with Pacbio long SMRT reads. LD 1 k 49 AdaptiveTh 0.005 KmerCovTh 5 MinOverlap 30 RemoveChimera 1 were used to generate raw backbone. Finally, blasr and consensus modules were implemented to produce the genome assembly. The final assembled genome was polished by pilon-1.22  to detect possible misassemblies. The final assembly was used for downstream analysis.
Demographic history inference using PSMC
To reconstruct the demographic history of Chinese Grouse, we used the Pairwise Sequentially Markovian Coalescent (PSMC 0.6.5)  model to infer the effective population sizes (Ne) of Chinese Grouse across genome sequences with SNP sites. We used the consensus autosomal sequence (in fastq format) for Chinese Grouse as the input for the PSMC modelling. To run the PSMC method, we used two parameters: the generation time (2.5 years) and the mutation rate per generation (0.3 × 10− 8 per nucleotide per year) of Chinese Grouse. The mutation rate (per nucleotide per year, μ) for Chinese Grouse was selected from studies of willow grouse (0.299 × 10− 8) and rock ptarmigan (0.310 × 10− 8) . PSMC was run for 100 iterations using an initial h = q ratio of 5 and the default time patterning. Bootstrapping was performed as in Li and Durbin (2011)  by resampling 500,000 bp chunks of the genome with replacement to perform 100 bootstrap replicates.
We used the maximum-entropy approach (MaxEnt, ver. 3.3.3 k)  to model Chinese Grouse habitat. A set of layers, or environmental variables, as well as a set of georeferenced occurrence data were used to produce a model of the ancient and present range of Chinese Grouse. We projected the present-day and three other periods Species Distribution Models (SDMs) based on the climate data during the current, mid-Holocene (6000 years BP), last glacial maximum (LGM, 21,000 years BP), and last interglacial (LIG, 120,000–140,000 years BP). Our models performed well (AUC value is 0.991) for the Chinese Grouse (Additional file 1). We obtained all the 19 bio-climatic variables from the WorldClim 1.4 database at a spatial resolution of 2.5 arc-minutes . The current climatic data layers were based on spatially interpolated values of temperature and precipitation from 1960 to 1990, whereas the variables were generated from the climate reconstructions based on CCSM4 for both the period mid-Holocene and LGM , and the variables for LIG were downscaled from simulations with CCSM4 . We used the same data set of grouse location (observation) records as from a reference by Lu et al. [5, 6]. The 41 location records we used to construct the SDMs were set in a spatial distance threshold of 0.083 decimal degrees (i.e. 5′, about 8 km) which can reduced spatial autocorrelation . Further parameter setting refer to Lu et al. . Furthermore, a 10% training presence logistic threshold was used to define the minimum probability of suitable habitat. A GAP  analysis between current distribution and coniferous forests were performed to assess conservation status. All data preparation and raster calculations were performed in ArcGIS 10.5.
Availability of data and materials
Sequencing data for the Chinese Grouse have been deposited in Short Read Archive under project number PRJNA588719. WorldClim – Global Climate Data is free climate data for ecological modeling and GIS at https://worldclim.org/data/v1.4/worldclim14.html and https://worldclim.org/data/v1.4/paleo1.4.html. These data can be used for mapping and spatial modeling. The data set of grouse location (observation) records is from a paper by Lu et al. 2012. The link is https://www.sciencedirect.com/science/article/pii/S0006320712001863.
Species Distribution Modelling
Pairwise Sequentially Markovian Coalescent
Last glacial maximum
Sun YH, Swenson JE, Fang Y, Klaus S, Scherzinger W. Population ecology of the Chinese grouse, Bonasa sewerzowi, in a fragmented landscape. Biol Conserv. 2003;110(2):0–184.
Storch I. Conservation status of grouse worldwide: An update. Wildl Biol. 2007;13(1):5–12.
The IUCN Red List of Threatened Species 2017:e.T22679497A119421341. https://dx.doi.org/10.2305/IUCN.UK.2017-3.RLTS.T22679497A119421341.en. Downloaded on 30 July 2020.
BirdLife International. Species factsheet: Bonasa sewerzowi. 2020. Downloaded from http://www.birdlife.org on 30/07/2020.
Lu N, Jia C-X, Lloyd H, Sun Y-H. Species-specific habitat fragmentation assessment, considering the ecological niche requirements and dispersal capability. Biol Conserv. 2012;152:102–9.
Lyu N, Sun Y-H. Predicting threat of climate change to the Chinese grouse on the Qinghai—Tibet plateau. Wildl Biol. 2014;20(2):73–82.
Hoorn C, Wesselingh F, Ter Steege H, Bermudez M, Mora A, Sevink J, Sanmartín I, Sanchez-Meseguer A, Anderson C, Figueiredo J. Amazonia through time: Andean uplift, climate change, landscape evolution, and biodiversity. Science. 2010;330(6006):927–31.
Ren G, Mateo RG, Liu J, Suchan T, Alvarez N, Guisan A, Conti E, Salamin N. Genetic consequences of quaternary climatic oscillations in the Himalayas: Primula tibetica as a case study based on restriction site-associated DNA sequencing. New Phytol. 2017;213(3):1500–12.
Qu Y, Tian S, Han N, Zhao H, Gao B, Fu J, Cheng Y, Song G, Ericson PG, Zhang YE. Genetic responses to seasonal variation in altitudinal stress: whole-genome resequencing of great tit in eastern Himalayas. Sci Rep. 2015;5(1):14256.
Zhan X, Zheng Y, Wei F, Bruford MW, Jia C. Molecular evidence for Pleistocene refugia at the eastern edge of the Tibetan plateau. Mol Ecol. 2011;20(14):3014–26.
He K, Jiang X. Sky islands of Southwest China. I: an overview of phylogeographic patterns. Chin Sci Bull. 2014;59(7):585–97.
Zhao S, Zheng P, Dong S, Zhan X, Wu Q, Guo X, Hu Y, He W, Zhang S, Fan W, et al. Whole-genome sequencing of giant pandas provides insights into demographic history and local adaptation. Nat Genet. 2013;45(1):67–71.
Zhou X, Wang B, Pan Q, Zhang J, Kumar S, Sun X, Liu Z, Pan H, Lin Y, Liu G, et al. Whole-genome sequencing of the snub-nosed monkey provides insights into folivory and evolutionary history. Nat Genet. 2014;46(12):1303–10.
Shi YF. A suggestion to improve the chronology of quaternary glaciations in China. J Glaciolgy Geocryol. 2002;24(6):687–92.
Zheng B, Xu Q, Shen Y. The relationship between climate change and quaternary glacial cycles on the Qinghai–Tibetan plateau: review and speculation. Quat Int. 2002;97(1):93–101.
Liu N. Studies on the southward migration of Severtzov's hazel grouse and it's Incontinuous distribution with the other species of hazel grouse. J Lanzhou Univ. 1994.
Wu C. Flora of Tibet. Beijing: Science Press; 1987.
Liu J, Michael M, Jim P, Lian-Ming G, Ram Chandra P, De-Zhu L. Geological and ecological factors drive cryptic speciation of yews in a biodiversity hotspot. New Phytol. 2013;199(4):1093–108.
Dynesius M, Jansson R. Evolutionary consequences of changes in species' geographical distributions driven by Milankovitch climate oscillations. Proc Natl Acad Sci. 2000;97(16):9115–20.
Hewitt G. The genetic legacy of the quaternary ice ages. Nature. 2000;405(6789):907.
Sandel B, Arge L, Dalsgaard B, Davies R, Gaston K, Sutherland W, Svenning J-C. The influence of Late Quaternary climate-change velocity on species endemism. Science. 2011;334(6056):660–4.
Weigelt P, Steinbauer MJ, Cabral JS, Kreft H. Late Quaternary climate change shapes island biodiversity. Nature. 2016;532(7597):99.
Hewitt G. Genetic consequences of climatic oscillations in the quaternary. Philosophical Transactions of the Royal Society B: Biological Sciences. 2004;359(1442):183–95.
Soltis DE, Morris AB, McLachlan JS, Manos PS, Soltis PS. Comparative phylogeography of unglaciated eastern North America. Mol Ecol. 2006;15(14):4261–93.
Wang C, Zhao X, Liu Z, Lippert PC, Graham SA, Coe RS, Yi H, Zhu L, Liu S, Li Y. Constraints on the early uplift history of the Tibetan plateau. Proc Natl Acad Sci. 2008;105(13):4987–92.
Chen S, Wu G, Zhang D, Gao Q, Duan Y, Zhang F, Chen S. Potential refugium on the Qinghai–Tibet plateau revealed by the chloroplast DNA phylogeography of the alpine species Metagentiana striata (Gentianaceae). Bot J Linn Soc. 2008;157(1):125–40.
Shi Y, Ge Y, Liu X, Li B, Yao T. Reconstruction of the 30–40 ka bp enhanced Indian monsoon climate based on geological records from the Tibetan plateau. Palaeogeogr Palaeoclimatol Palaeoecol. 2001;169(1):69–83.
Yu G, Gui F, Shi Y, Zheng Y. Late marine isotope stage 3 palaeoclimate for East Asia: A data–model comparison. Palaeogeogr Palaeoclimatol Palaeoecol. 2007;250(1):167–83.
Li B. The last greatest lakes on the Xizang (Tibetan) plateau. Acta Geograph Sin. 2000;55:174–85.
Tang LY. Pollen evidence of climate during the last glacial maximum in eastern Tibetan plateau. J Glaciolgy Geocryol. 1998.
Li H, Durbin R. Inference of human population history from individual whole-genome sequences. Nature. 2011;475(7357):493.
Dai Z, Du J, Li J, Li W, Chen J. Runoff characteristics of the Changjiang River during 2006: effect of extreme drought and the impounding of the three gorges dam. Geophys Res Lett. 2008;35(7):521–39.
Scott JM, Davis F, Csuti B, Noss R, Butterfield B, Groves C, Anderson H, Caicco S, D'Erchia F, Edwards TC Jr. Gap analysis: a geographic approach to protection of biological diversity. Wildl Monogr. 1993;123:3–41.
Hillier LW, Miller W, Birney E, Warren W, Hardison RC, Ponting CP, Bork P, Burt DW, Groenen MAM, Delany ME, et al. Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution. Nature. 2004;432(7018):695–716.
Tong G, Li Y, Yang Z. Records of quaternary palynoflora events in China. Acta Micropalaeontologica Sinica. 2000;17(2):186-97.
Tong G, Zhang J, Fan S. Distribution of quaternary palynoflora in China. Marine Geol Quaternary Geol. 1992;12:45–56.
Huang K, Zhang Y, Zhu J, Liu Y, Zu J, Zhang J. The influences of climate change and human activities on vegetation dynamics in the Qinghai-Tibet plateau. Remote Sens. 2016;8(10):876.
Davis MB, Shaw RG. Range shifts and adaptive responses to quaternary climate change. Science. 2001;292(5517):673–9.
Li L, Abbott RJ, Liu B, Sun Y, Li L, Zou J, Wang X, Miehe G, Liu J. Pliocene intraspecific divergence and Plio-Pleistocene range expansions within Picea likiangensis (Lijiang spruce), a dominant forest tree of the Qinghai-Tibet plateau. Mol Ecol. 2013;22(20):5237–55.
Myers N, Mittermeier RA, Mittermeier CG, Da Fonseca GA, Kent J. Biodiversity hotspots for conservation priorities. Nature. 2000;403(6772):853.
Ursenbacher S, Carlsson M, ., Helfer V, Tegelström H, Fumagalli L. Phylogeography and Pleistocene refugia of the adder (Vipera berus) as inferred from mitochondrial DNA sequence data. Mol Ecol 2006, 15(11):3425–3437.
Xing Y, Ree RH. Uplift-driven diversification in the Hengduan Mountains, a temperate biodiversity hotspot. Proc Natl Acad Sci U S A. 2017;114(17):E3444.
Kozma R. Inferring demographic history and speciation of grouse using whole genome sequences[D]. Acta Universitatis Upsaliensis, 2016.
Nadachowskabrzyska K, Li C, Smeds L, Zhang G, Ellegren H. Temporal dynamics of avian populations during Pleistocene revealed by whole-genome sequences. Curr Biol. 2015;25(10):1375–80.
Kozma R, Lillie M, Benito BM, Svenning JC, Höglund J. Past and potential future population dynamics of three grouse species using ecological and whole genome coalescent modeling. Ecology Evolution. 2018;8(13):6671–81.
Sun Y, An Z. Late Pliocene-Pleistocene changes in mass accumulation rates of eolian deposits on the central Chinese Loess Plateau. J Geophysical Res: Atmospheres. 2005;110(D23101).
Sun Y, An Z. Chinese loess plateau mass accumulation rate data. IGBP PAGES/World Data Center for Paleoclimatology Data Contribution Series; 2006. p. 34.
Kettle CJ, Koh LP. Global forest fragmentation: CABI; 2014.
Pfeifer M, Lefebvre V, Peres C, Banks-Leite C, Wearn O, Marsh C, Butchart S, Arroyo-Rodríguez V, Barlow J, Cerezo A. Creation of forest edges has a global impact on forest vertebrates. Nature. 2017;551(7679):187.
Lu N, Jing Y, Lloyd H, Sun Y-H. Assessing the distributions and potential risks from climate change for the Sichuan jay (Perisoreus internigrans). Condor. 2012;114(2):365–76.
Pounds JA, Fogden MP, Campbell JH. Biological response to climate change on a tropical mountain. Nature. 1999;398(6728):611.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.
Ye C, Hill CM, Wu S, Ruan J, Ma Z. DBG2OLC: efficient assembly of large genomes using long erroneous reads of the third generation sequencing technologies. Sci Rep. 2016;6(X):31900.
Walker BJ, Abeel T, Shea T, Priest M, Abouelliel A, Sakthikumar S, Cuomo CA, Zeng Q, Wortman J, Young SK. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS One. 2014;9(11):e112963.
Kozma R, Melsted P, Magnússon KP, Höglund J. Looking into the past - the reaction of three grouse species to climate change over the last million years using whole genome sequences. Mol Ecol. 2016;25(2):570–80.
Phillips SJ, Anderson RP, Schapire RE. Maximum entropy modeling of species geographic distributions. Ecol Model. 2006;190(3):231–59.
Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A. Very high resolution interpolated climate surfaces for global land areas. Int J Climatol. 2005;25(15):1965–78.
Otto-Bliesner BL, Marshall SJ, Overpeck JT, Miller GH, Aixue H. Simulating Arctic climate warmth and icefield retreat in the last interglaciation. Science. 2006;311(5768):1751–3.
Dormann CM, McPherson JB, Araújo M, Bivand R, Bolliger J, Carl G, Davies RG, Hirzel A, Jetz W, Daniel Kissling W. Methods to account for spatial autocorrelation in the analysis of species distributional data: a review. Ecography. 2007;30(5):609–28.
We thank Annoroad Gene Technology in Beijing for performing the whole genome sequencing.
Funding for this project was provided by National Natural Science Foundation of China (NSFC: 31520103903 to SYH, JH). Field work and de novo sequences were supported by this funding. Open access funding provided by Uppsala University.
Ethics approval and consent to participate
This study was approved by the Animal Ethics Committee of the Institute of Zoology, Chinese Academy of Sciences. The blood sample was collected from a wild Chinese Grouse at Lianhuashan National Nature Reserve which was released to the wild. The procedure of blood collection was in strict accordance with the Animal Ethics Procedures and Guidelines of the People’s Republic of China.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Song, K., Gao, B., Halvarsson, P. et al. Genomic analysis of demographic history and ecological niche modeling in the endangered Chinese Grouse Tetrastes sewerzowi. BMC Genomics 21, 581 (2020). https://doi.org/10.1186/s12864-020-06957-5
- Demographic history
- Qinghai-Tibet plateau
- Chinese grouse
- Tetrastes sewerzowi
- Whole genome sequence