Genome-wide analysis of cold adaptation in indigenous Siberian populations - PubMed (original) (raw)

. 2014 May 21;9(5):e98076.

doi: 10.1371/journal.pone.0098076. eCollection 2014.

Luca Pagani 2, Tiago Antao 3, Daniel J Lawson 4, Christina A Eichstaedt 1, Bryndis Yngvadottir 1, Ma Than Than Shwe 5, Joseph Wee 5, Irene Gallego Romero 6, Srilakshmi Raj 7, Mait Metspalu 8, Richard Villems 8, Eske Willerslev 9, Chris Tyler-Smith 10, Boris A Malyarchuk 11, Miroslava V Derenko 11, Toomas Kivisild 12

Affiliations

Genome-wide analysis of cold adaptation in indigenous Siberian populations

Alexia Cardona et al. PLoS One. 2014.

Abstract

Following the dispersal out of Africa, where hominins evolved in warm environments for millions of years, our species has colonised different climate zones of the world, including high latitudes and cold environments. The extent to which human habitation in (sub-)Arctic regions has been enabled by cultural buffering, short-term acclimatization and genetic adaptations is not clearly understood. Present day indigenous populations of Siberia show a number of phenotypic features, such as increased basal metabolic rate, low serum lipid levels and increased blood pressure that have been attributed to adaptation to the extreme cold climate. In this study we introduce a dataset of 200 individuals from ten indigenous Siberian populations that were genotyped for 730,525 SNPs across the genome to identify genes and non-coding regions that have undergone unusually rapid allele frequency and long-range haplotype homozygosity change in the recent past. At least three distinct population clusters could be identified among the Siberians, each of which showed a number of unique signals of selection. A region on chromosome 11 (chr11:66-69 Mb) contained the largest amount of clustering of significant signals and also the strongest signals in all the different selection tests performed. We present a list of candidate cold adaption genes that showed significant signals of positive selection with our strongest signals associated with genes involved in energy regulation and metabolism (CPT1A, LRP5, THADA) and vascular smooth muscle contraction (PRKG1). By employing a new method that paints phased chromosome chunks by their ancestry we distinguish local Siberian-specific long-range haplotype signals from those introduced by admixture.

PubMed Disclaimer

Conflict of interest statement

Competing Interests: The authors have declared that no competing interests exist.

Figures

Figure 1

Figure 1. Population structure analyses in Siberian populations.

(A) Averaged sampling locations of the Siberian populations genotyped in this study. (B) Principal component analysis of Siberian populations and reference populations from West Asia (European) and Southeast Asia (Vietnamese). Each dot in the plot represents an individual. The PC axes were rotated 180 degrees anti-clockwise to emphasize the similarity to the geographic map of Eurasia. (C) ADMIXTURE analysis at K = 4. (D) Coancestry heatmap for the Siberian individuals and reference populations (Europe, Vietnamese) output by ChromoPainter/fineSTRUCTURE. The heatmap shows the number of shared genetic chunks between the individuals. The raw data is shown on the bottom left and the aggregated data is shown on the upper right of the heatmap. Adjacent to the heatmap is also the ADMIXTURE plot of the respective individuals. To the left is the maximum a posteriori (MAP) tree generated by fineSTRUCTURE which shows the groupings of the different populations. The following abbreviations are used in the Figure: ALT, Altai-Kizhi; BUR, Buryats; CEU, European; CHB, Han Chinese; CHK, Chukchi; E/SE Asia, East/Southeast Asia; ESK, Eskimo; EVN, Evens; EVK, Evenks; KRK, Koryaks; SHR, Shors; TEL, Teleuts; VTN, Vietnamese; X1, Northeastern Siberian admixed individuals; X2, Southern and Central Siberian admixed individuals; YKT, Yakuts. Reference populations are labelled in italics. Outliers removed in downstream analysis are blacked out in the tree.

Figure 2

Figure 2. Zooming on the region on chromosome 11 that showed the strongest signals in the Northeastern Siberian populations.

The |iHS|, XP-EHH and PBS scores that showed strong signals (amongst top 10 ranking windows) in the Northeastern Siberian populations are shown in the upper panels for the different Siberian populations. The pink highlights mark the windows present in the top 10 ranking windows in the Northeastern Siberian populations. The rankings are marked on the respective windows in the respective test panels. Protein-coding genes present in the 4 Mb region are shown under the test plots with the genes present in our predefined cold adaptation list marked in bold font (CPT1A and LRP5). The position on the chromosome is given in Mb. The paintings of the phased chromosomes for the region in the Northeastern and Central Siberian individuals are shown underneath the Position legend. The aggregated ancestral probabilities from ChromoPainter for the Northeastern Siberian individuals are displayed below the paintings. The dotted red line shows the threshold for the upper 5% tail (0.75) and the black dotted line shows the mean of the genome-wide probabilities distribution of the Northeastern Siberian ancestor (red).

Figure 3

Figure 3. Genome-wide distribution of significant selection windows.

The distribution of the significant windows of the three different selection tests results (A) XP-EHH, (B) iHS and (C) PBS for the Northeastern Siberian populations. The x-axis show the count of significant windows (Tables S4–S6 in File S2) per 3 Mb region and the y-axis show the amount of 3 Mb regions that had the respective count displayed on the x-axis. In all the different selection tests, the region on chromosome 11:66–69 Mb is the region with the most significant selection windows (highlighted in red).

Similar articles

Cited by

References

    1. Goebel T (1999) Pleistocene human colonization of Siberia and peopling of the Americas: An ecological approach. Evolutionary Anthropology: Issues, News, and Reviews 8: 208–227.
    1. Goebel T, Waters MR, O'Rourke DH (2008) The Late Pleistocene Dispersal of Modern Humans in the Americas. Science 319: 1497–1502. - PubMed
    1. Pitulko VV, Nikolsky PA, Girya EY, Basilyan AE, Tumskoy VE, et al. (2004) The Yana RHS site: humans in the Arctic before the last glacial maximum. Science 303: 52–56. - PubMed
    1. Richards MP, Pettitt PB, Stiner MC, Trinkaus E (2001) Stable isotope evidence for increasing dietary breadth in the European mid-Upper Paleolithic. Proceedings of the National Academy of Sciences 98: 6528–6532. - PMC - PubMed
    1. Clark PU, Dyke AS, Shakun JD, Carlson AE, Clark J, et al. (2009) The Last Glacial Maximum. Science 325: 710–714. - PubMed

Publication types

MeSH terms

LinkOut - more resources