Genetic structure in the Sherpa and neighboring Nepalese populations

Ethno-geographic distribution and histopathological classification of nasopharyngeal carcinoma in a single center in Nepal

INTRODUCTION

Nasopharyngeal carcinoma (NPC) shows geographic and ethnic variation with specific etiopathogenesis. This study characterized the distribution of NPC patients stratified by ethnicity, geography, and histology in a tertiary-level cancer center in Nepal.

METHODS

A descriptive retrospective study was designed to analyze NPC cases from different regions among patients visiting the hospital from 2016 to 2021. Demographic and clinical information was obtained from medical records. Ethical approval was granted by the Nepal Health Research Council (NHRC). Data analyses and visualization were carried out with R software.

RESULTS

During the six-year period, a total of 65 individuals were diagnosed with NPC, comprising 42 males and 23 females. Patient median age was 43 years (range 11-85 years). A bimodal age distribution of cases was observed with peaks in patients aged 30-39 years and 50-59 years. Of the NPC patients studied, 29 were from Koshi Province, with 7 cases from Ilam district and 6 cases from Morang district. There were 18 patients in Bagmati Province, and Kathmandu district had the highest number of cases within this region, with 8 patients. The highest proportion of cases were observed among patients of Janajati ethnicity (60%), including Rai, Limbu, and Sherpa people. Histologically, undifferentiated non-keratinizing NPC was the most commonly observed subtype, accounting for 43.1% of cases, followed by 20% differentiated non-keratinizing NPC and 4.6% keratinizing NPC across the entire sample population. The majority of patients (75.3%) were diagnosed at an advanced stage (stage III or IV) with none diagnosed at stage I.

CONCLUSIONS

In our study, most cases of NPC occurred in patients from provinces in eastern Nepal (Koshi province), and of the Janajati ethnic community. The most common histological subtype was undifferentiated non-keratinizing carcinoma. Further epidemiological studies could address differences in prevalence and the challenge of late presentation of NPC patients in Nepal.

Maternal genetic history of ancient Tibetans over the past 4000 years

The settlement of the Tibetan Plateau epitomizes human adaptation to a high-altitude environment that poses great challenges to human activity. Here, we reconstruct a 4000-year maternal genetic history of Tibetans using 128 ancient mitochondrial genome data from 37 sites in Tibet. The phylogeny of haplotypes M9a1a, M9a1b, D4g2, G2a'c, and D4i show that ancient Tibetans share the most recent common ancestor with ancient Middle and Upper Yellow River populations around the Early and Middle Holocene. In addition, the connections between Tibetans and Northeastern Asians vary over the past 4000 years, with a stronger matrilineal connection between the two during 4000 BP-3000 BP, and a weakened connection after 3000 BP, that are coincident with climate change, followed by a reinforced connection after the Tubo period (1400 BP-1100 BP). Besides, an over 4000-year matrilineal continuity is observed in some of the maternal lineages. We also find the maternal genetic structure of ancient Tibetans is correlated to the geography and interactions between ancient Tibetans and ancient Nepal and Pakistan populations. Overall, the maternal genetic history of Tibetans can be characterized as a long-term matrilineal continuity with frequent internal and external population interactions that are dynamically shaped by geography, climate changes, as well as historical events.

The matrilineal ancestry of Nepali populations

The Tibetan plateau and high mountain ranges of Nepal are one of the challenging geographical regions inhabited by modern humans. While much of the ethnographic and population-based genetic studies were carried out to investigate the Tibetan and Sherpa highlanders, little is known about the demographic processes that enabled the colonization of the hilly areas of Nepal. Thus, the present study aimed to investigate the past demographic events that shaped the extant Nepalese genetic diversity using mitochondrial DNA (mtDNA) variations from ethnic Nepalese groups. We have analyzed mtDNA sequences of 999 Nepalese and compared data with 38,622 published mtDNA sequences from rest of the world. Our analysis revealed that the genomic landscapes of prehistoric Himalayan settlers of Nepal were similar to that of the low-altitude extant Nepalese (LAN), especially Newar and Magar population groups, but differ from contemporary high-altitude Sherpas. LAN might have derived their East Eurasian ancestry mainly from low-altitude Tibeto-Burmans, who likely have migrated from East Asia and assimilated across the Eastern Himalayas extended from the Eastern Nepal to the North-East of India, Bhutan, Tibet and Northern Myanmar. We also identified a clear genetic sub-structure across different ethnic groups of Nepal based on mtDNA haplogroups and ectodysplasin-A receptor (EDAR) gene polymorphism. Our comprehensive high-resolution mtDNA-based genetic study of Tibeto-Burman communities reconstructs the maternal origins of prehistoric Himalayan populations and sheds light on migration events that have brought most of the East Eurasian ancestry to the present-day Nepalese population.

Population history and genome wide association studies of birth weight in a native high altitude Ladakhi population

Pathological low birth weight due to fetal growth restriction (FGR) is an important predictor of adverse obstetric and neonatal outcomes. It is more common amongst native lowlanders when gestating in the hypoxic environment of high altitude, whilst populations who have resided at high altitude for many generations are relatively protected. Genetic study of pregnant populations at high altitude permits exploration of the role of hypoxia in FGR pathogenesis, and perhaps of FGR pathogenesis more broadly. We studied the umbilical cord blood DNA of 316 neonates born to pregnant women managed at the Sonam Norboo Memorial Hospital, Ladakh (altitude 3540m) between February 2017 and January 2019. Principal component, admixture and genome wide association studies (GWAS) were applied to dense single nucleotide polymorphism (SNP) genetic data, to explore ancestry and genetic predictors of low birth weight. Our findings support Tibetan ancestry in the Ladakhi population, with subsequent admixture with neighboring Indo-Aryan populations. Fetal growth protection was evident in Ladakhi neonates. Although no variants achieved genome wide significance, we observed nominal association of seven variants across genes (ZBTB38, ZFP36L2, HMGA2, CDKAL1, PLCG1) previously associated with birthweight.

The blunted vascular endothelial growth factor-A (VEGF-A) response to high-altitude hypoxia and genetic variants in the promoter region of the VEGFA gene in Sherpa highlanders

Background

Sherpa highlanders demonstrate extraordinary tolerance to hypoxia at high altitudes, which may be achieved by mechanisms promoting microcirculatory blood flow and capillary density at high altitudes for restoring oxygen supply to tissues. Vascular endothelial growth factors (VEGFs) are important signaling proteins involved in vasculogenesis and angiogenesis which are stimulated by hypoxia. We hypothesize that the VEGF-A, the major member of the VEGF family, and the gene encoding VEGF-A (VEGFA) play a part in the adaptation to high-altitude hypoxia in Sherpa highlanders.

Methods

Fifty-one Sherpa highlanders in Namche Bazaar village at a high altitude of 3,440 meters (m) above sea level and 76 non-Sherpa lowlanders in Kathmandu city at 1,300 m in Nepal were recruited for the study. Venous blood was sampled to obtain plasma and extract DNA from each subject. The plasma VEGF-A concentrations were measured and five single-nucleotide polymorphisms (SNPs, rs699947, rs833061, rs1570360, rs2010963, and rs3025039) in the VEGFA were genotyped. The VEGF-A levels and allelic frequencies of the SNPs were compared between the two populations.

Results

A significant difference in oxygen saturation (SpO₂) was observed between the two ethnic groups locating at different elevations (93.7 ± 0.2% in Sherpas at 3,440 m vs. 96.7 ± 0.2% in non-Sherpas at 1,300 m, P < 0.05). The plasma VEGF-A concentration in the Sherpas at high altitude was on the same level as that in the non-Sherpas at low altitude (262.8 ± 17.9 pg/ml vs. 266.8 ± 21.8 pg/ml, P = 0.88). This result suggested that the plasma VEGF-A concentration in Sherpa highlanders was stable despite a high-altitude hypoxic stimulus and that therefore the Sherpas exhibited a phenotype of blunted response to hypoxic stress. Moreover, the allele frequencies of the SNPs rs699947, rs833061, and rs2010963 in the promoter region of the VEGFA were different between the Sherpa highlanders and non-Sherpa lowlanders (corrected P values = 3.30 ×10^-5, 4.95 ×10^-4, and 1.19 ×10^-7, respectively).

Conclusions

Sherpa highlanders exhibited a blunted VEGF-A response to hypoxia at high altitudes, which was speculated to be associated with the distinctive genetic variations of the SNPs and haplotype in the promoter region of VEGFA in Sherpa highlanders.

Altitude-Induced Sleep Apnea Is Highly Dependent on Ethnic Background (Sherpa Vs. Tamang)

Heiniger, Grégory, Simon Walbaum, Claudio Sartori, Alban Lovis, Marco Sazzini, Andrew Wellman, and Raphael Heinzer. Altitude-Induced Sleep Apnea Is Highly Dependent on Ethnic Background (Sherpa Vs. Tamang). High Alt Med Biol. 23:165-172, 2022. Rationale: High altitude-induced hypocapnic alkalosis generates central sleep apnea (CSA). In Nepal, two ethnic groups live at medium-to-high altitude: Tamangs originate from low-altitude Tibeto-Burman populations, whereas Sherpas descend from high-altitude Tibetans. Objective: To compare apnea severity at low and high altitude between Sherpas and Tamangs. Methods: Polygraphy recordings, including airflow and oxygen saturation, were performed in Nepal at "low" (2,030 m) and "high" (4,380 m) altitudes. Resting ventilation (V̇E) and mixed-exhaled CO₂ (F_ECO₂) were also measured at the same altitudes. Differences in apnea-hypopnea index (AHI), oxygen desaturation index (ODI), and % of nocturnal periodic breathing (NPB) at the two altitudes were compared between ethnicities. Measurements and Main Results: Twenty Sherpas and 20 Tamangs were included (males, median [interquartile range] age: 24.5 [21.5-27.8] years vs. 26.0 [21.5-39.8] years, body mass index: 23.9 [22.1-26.1] kg/m² vs. 25.21 [20.6-27.6] kg/m²). Compared with Tamangs, Sherpas showed a lower increase in AHI (+7.5 [2.6-17.2]/h vs. +31.5 [18.2-57.3]/h, p < 0.001), ODI (+13.8 [5.5-28.2]/h vs. +42.0 [22.6-77.6]/h, p < 0.001), and NPB proportion (+0.9 [0-3.5]% vs. +12.8 [3.1-27.4]%, p < 0.001) from low to high altitude. Resting V̇E was higher in Sherpas versus Tamangs at both low (8.45 [6.89-10.70] l/min vs. 6.3 [4.9-8.3] l/min, p = 0.005) and high (9.7 [8.5-11] l/min vs. 8.74 [7.39-9.73] l/min, p = 0.020) altitudes, whereas the mean ± standard deviation F_ECO₂ decrease between low and high altitude was greater in Tamangs versus Sherpas (-0.50% ± 0.44% vs. -0.80% ± 0.33%, p < 0.023). Conclusion: Overall, altitude-adapted Sherpas showed a 3.2-times smaller increase in sleep-disordered breathing between low and high altitude compared with Tamangs, and higher ventilation and a smaller drop in F_ECO₂ at high altitude. These data suggest that genetic differences in breathing control can be protective against CSA.

Ancient genomes from the Himalayas illuminate the genetic history of Tibetans and their Tibeto-Burman speaking neighbors

Present-day Tibetans have adapted both genetically and culturally to the high altitude environment of the Tibetan Plateau, but fundamental questions about their origins remain unanswered. Recent archaeological and genetic research suggests the presence of an early population on the Plateau within the past 40 thousand years, followed by the arrival of subsequent groups within the past 10 thousand years. Here, we obtain new genome-wide data for 33 ancient individuals from high elevation sites on the southern fringe of the Tibetan Plateau in Nepal, who we show are most closely related to present-day Tibetans. They derive most of their ancestry from groups related to Late Neolithic populations at the northeastern edge of the Tibetan Plateau but also harbor a minor genetic component from a distinct and deep Paleolithic Eurasian ancestry. In contrast to their Tibetan neighbors, present-day non-Tibetan Tibeto-Burman speakers living at mid-elevations along the southern and eastern margins of the Plateau form a genetic cline that reflects a distinct genetic history. Finally, a comparison between ancient and present-day highlanders confirms ongoing positive selection of high altitude adaptive alleles.

Genomic history and forensic characteristics of Sherpa highlanders on the Tibetan Plateau inferred from high-resolution InDel panel and genome-wide SNPs

Sherpa people, one of the high-altitude hypoxic adaptive populations, mainly reside in Nepal and the southern Tibet Autonomous Region. The genetic origin and detailed evolutionary profiles of Sherpas remain to be further explored and comprehensively characterized. Here we analyzed the newly-generated InDel genotype data from 628 Dingjie Sherpas by merging with 4222 worldwide InDel profiles and collected genome-wide SNP data (approximately 600K SNPs) from 1612 individuals in 191 modern and ancient populations to explore and reconstruct the fine-scale genetic structure of Sherpas and their relationships with nearby modern and ancient East Asians based on the shared alleles and haplotypes. The forensic parameters of 57 autosomal InDels (A-InDels) included in our used new-generation InDel amplification system showed that this focused InDel panel is informative and polymorphic in Dingjie Sherpas, suggesting that it can be used as the supplementary tool for forensic personal identification and parentage testing in Dingjie Sherpas. Descriptive findings from the PCA, ADMIXTURE, and TreeMix-based phylogenies suggested that studied Nepal Sherpas showed excess allele sharing with neighboring Tibeto-Burman Tibetans. Furthermore, patterns of allele sharing in f-statistics demonstrated that Nepal Sherpas had a different evolutionary history compared with their neighbors from Nepal (Newar and Gurung) but showed genetic similarity with 2700-year-old Chokhopani and modern Tibet Tibetans. QpAdm/qpGraph-based admixture sources and models further showed that Sherpas, core Tibetans, and Chokhopani formed one clade, which could be fitted as having the main ancestry from late Neolithic Qijia millet farmers and other deep ancestries from early Asians. Chromosome painting profiles and shared IBD fragments inferred from fineSTRUCTURE and ChromoPainter not only confirmed the abovementioned genomic affinity patterns but also revealed the fine-scale genetic microstructures among Sino-Tibetan speakers. Finally, natural-selection signals revealed via iHS, nSL and iHH12 showed natural selection signatures associated with disease susceptibility in Sherpas. Generally, we provided the comprehensive landscape of admixture and evolutionary history of Sherpa people based on the shared alleles and haplotypes from the InDel-based genotype data and high-density genome-wide SNP data. The more detailed genetic landscape of Sherpa people should be further confirmed and characterized via ancient genomes or single-molecule real-time sequencing technology.

Tracing the Genetic Legacy of the Tibetan Empire in the Balti

The rise and expansion of Tibetan Empire in the 7th to 9th centuries AD affected the course of history across East Eurasia, but the genetic impact of Tibetans on surrounding populations remains undefined. We sequenced 60 genomes for four populations from Pakistan and Tajikistan to explore their demographic history. We showed that the genomes of Balti people from Baltistan comprised 22.6–26% Tibetan ancestry. We inferred a single admixture event and dated it to about 39–21 generations ago, a period that postdated the conquest of Baltistan by the ancient Tibetan Empire. The analyses of mitochondrial DNA, Y, and X chromosome data indicated that both ancient Tibetan males and females were involved in the male-biased dispersal. Given the fact that the Balti people adopted Tibetan language and culture in history, our study suggested the impact of Tibetan Empire on Baltistan involved dominant cultural and minor demic diffusion.

Death Zone Weather Extremes Mountaineers Have Experienced in Successful Ascents

Background

Few data are available on mountaineers’ survival prospects in extreme weather above 8000 m (the Death Zone). We aimed to assess Death Zone weather extremes experienced in climbing-season ascents of Everest and K2, all winter ascents of 8000 m peaks (8K) in the Himalayas and Karakoram, environmental records of human survival, and weather extremes experienced with and without oxygen support.

Materials and Methods

We analyzed 528 ascents of 8K peaks: 423 non-winter ascents without supplemental oxygen (Everest–210, K2–213), 76 ascents in winter without oxygen, and 29 in winter with oxygen. We assessed environmental conditions using the ERA5 dataset (1978–2021): barometric pressure (BP), temperature (Temp), wind speed (Wind), wind chill equivalent temperature (WCT), and facial frostbite time (FFT).

Results

The most extreme conditions that climbers have experienced with and without supplemental oxygen were: BP 320 hPa (winter Everest) vs. 329 hPa (non-winter Everest); Temp –41°C (winter Everest) vs. –45°C (winter Nanga Parbat); Wind 46 m⋅s^–1 (winter Everest) vs. 48 m⋅s^–1 (winter Kangchenjunga). The most extreme combined conditions of BP ≤ 333 hPa, Temp ≤ −30°C, Wind ≥ 25 m⋅s^–1, WCT ≤ −54°C and FFT ≤ 3 min were encountered in 14 ascents of Everest, two without oxygen (late autumn and winter) and 12 oxygen-supported in winter. The average extreme conditions experienced in ascents with and without oxygen were: BP 326 ± 3 hPa (winter Everest) vs. 335 ± 2 hPa (non-winter Everest); Temp −40 ± 0°C (winter K2) vs. −38 ± 5°C (winter low Karakoram 8K peaks); Wind 36 ± 7 m⋅s^–1 (winter Everest) vs. 41 ± 9 m⋅s^–1 (winter high Himalayan 8K peaks).

Conclusions

Cross-Species Insights Into Genomic Adaptations to Hypoxia

Over millions of years, vertebrate species populated vast environments spanning the globe. Among the most challenging habitats encountered were those with limited availability of oxygen, yet many animal and human populations inhabit and perform life cycle functions and/or daily activities in varying degrees of hypoxia today. Of particular interest are species that inhabit high-altitude niches, which experience chronic hypobaric hypoxia throughout their lives. Physiological and molecular aspects of adaptation to hypoxia have long been the focus of high-altitude populations and, within the past decade, genomic information has become increasingly accessible. These data provide an opportunity to search for common genetic signatures of selection across uniquely informative populations and thereby augment our understanding of the mechanisms underlying adaptations to hypoxia. In this review, we synthesize the available genomic findings across hypoxia-tolerant species to provide a comprehensive view of putatively hypoxia-adaptive genes and pathways. In many cases, adaptive signatures across species converge on the same genetic pathways or on genes themselves [i.e., the hypoxia inducible factor (HIF) pathway). However, specific variants thought to underlie function are distinct between species and populations, and, in most cases, the precise functional role of these genomic differences remains unknown. Efforts to standardize these findings and explore relationships between genotype and phenotype will provide important clues into the evolutionary and mechanistic bases of physiological adaptations to environmental hypoxia.

Trends and clinicopathological characteristics of oral squamous cell carcinomas reported at a tertiary cancer hospital in Nepal during 1999 to 2009

OBJECTIVE

Reliable data describing the trends and clinicopathological characteristics of oral squamous cell carcinoma (OSCC) in the Nepalese population are very limited. The current study aimed to examine the demographics, trend, and clinicopathological characteristics of OSCC reported to the main referral/tertiary cancer hospital, the B.P. Koirala Memorial Cancer Hospital (BPKMCH) in Nepal for a period of 11 years (1999-2009).

MATERIAL AND METHODS

This is a cross-sectional study. Data were retrieved retrospectively from hospital register maintained in the Department of Ear, Nose, Throat at BPKMCH, categorized into demographic and clinicopathological variables and SPSS (V25) was used for the analysis.

RESULTS

In a period of 11 years, 3,452 cases of head and neck cancer were registered at the Department of Ear, Nose, Throat, BPKMCH. Out of 1,111 oral cancer cases, 1,081 (97.3%) were OSCC. A trend for increasing number of OSCCs presenting to BPKMCH was observed during that period. OSCC was found to be more common among males (73.0%), Brahmin/Chhetri ethnic groups (33.0%), in age group of 51-60 years (31.9%), and in Terai region (62.0%). Tongue (42.8%) was the most common site, followed by buccal mucosa (27.2%). Nevertheless, when stratified with respect to the geographical location and ethnicity, buccal mucosa was the most common site for OSCC in Terai region (63.9%, p = .002) and in Madhesi ethnic group (34.2%, p < .001). Majority of OSCC cases were diagnosed at advanced stage (49.7%, Stage IV) and received a combination therapy (42.0%).

CONCLUSIONS

Hospital-based records can provide valuable information on disease characteristics in countries like Nepal. This study revealed that the clinicopathological characteristics of OSCC in Nepal follow the global trend. Nevertheless, relationship between specific intraoral sites for OSCC with geographic location and ethnic groups is an interesting observation and requires further population-based studies to clarify these findings.

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.

Convergent evolution in human and domesticate adaptation to high-altitude environments

Humans and their domestic animals have lived and thrived in high-altitude environments worldwide for thousands of years. These populations have developed a number of adaptations to survive in a hypoxic environment, and several genomic studies have been conducted to identify the genes that drive these adaptations. Here, we discuss the various adaptations and genetic variants that have been identified as adaptive in human and domestic animal populations and the ways in which convergent evolution has occurred as these populations have adapted to high-altitude environments. We found that human and domesticate populations have adapted to hypoxic environments in similar ways. Specific genes and biological pathways have been involved in high-altitude adaptation for multiple populations, although the specific variants differ between populations. Additionally, we found that the gene EPAS1 is often a target of selection in hypoxic environments and has been involved in multiple adaptive introgression events. High-altitude environments exert strong selective pressures, and human and animal populations have evolved in convergent ways to cope with a chronic lack of oxygen.

This article is part of the theme issue ‘Convergent evolution in the genomics era: new insights and directions'.

Genetic variants at the EGLN1 locus associated with high-altitude adaptation in Tibetans are absent or found at low frequency in highland Andeans

EGLN1 encodes the hypoxia-inducible factor (HIF) pathway prolyl hydroxylase 2 (PHD2) that serves as an oxygen-sensitive regulator of HIF activity. The EGLN1 locus exhibits a signature of positive selection in Tibetan and Andean populations and is associated with hemoglobin concentration in Tibetans. Recent reports provide evidence for functional roles of protein-coding variants within the first exon of EGLN1 (rs186996510, rs12097901) that are linked to an adaptive signal in Tibetans, yet whether these same variants are present and contribute to adaptation in Andean highlanders is unknown. We determined the frequencies of these adaptive Tibetan alleles in Quechua Andeans resident at high altitude (4,350 m) in addition to individuals of Nepali ancestry resident at sea level. The rs186996510 C (minor) allele previously found at high frequency in Tibetans is absent in Andean (G: 100%) and rare among Nepali (C: 11.8%, G: 88.2%) cohorts. The minor G allele of rs12097901 is found at similarly low frequencies in Andeans (G: 12.7%, C: 87.3%) and Nepalis (G: 23.5%, C: 76.5%) compared to Tibetans. These results suggest that adaptation involving EGLN1 in Andeans involves different mechanisms than those described in Tibetans. The precise Andean adaptive variants remain to be determined.

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.

Weighted likelihood inference of genomic autozygosity patterns in dense genotype data

Background

Genomic regions of autozygosity (ROA) arise when an individual is homozygous for haplotypes inherited identical-by-descent from ancestors shared by both parents. Over the past decade, they have gained importance for understanding evolutionary history and the genetic basis of complex diseases and traits. However, methods to infer ROA in dense genotype data have not evolved in step with advances in genome technology that now enable us to rapidly create large high-resolution genotype datasets, limiting our ability to investigate their constituent ROA patterns.

Methods

We report a weighted likelihood approach for inferring ROA in dense genotype data that accounts for autocorrelation among genotyped positions and the possibilities of unobserved mutation and recombination events, and variability in the confidence of individual genotype calls in whole genome sequence (WGS) data.

Results

Forward-time genetic simulations under two demographic scenarios that reflect situations where inbreeding and its effect on fitness are of interest suggest this approach is better powered than existing state-of-the-art methods to infer ROA at marker densities consistent with WGS and popular microarray genotyping platforms used in human and non-human studies. Moreover, we present evidence that suggests this approach is able to distinguish ROA arising via consanguinity from ROA arising via endogamy. Using subsets of The 1000 Genomes Project Phase 3 data we show that, relative to WGS, intermediate and long ROA are captured robustly with popular microarray platforms, while detection of short ROA is more variable and improves with marker density. Worldwide ROA patterns inferred from WGS data are found to accord well with those previously reported on the basis of microarray genotype data. Finally, we highlight the potential of this approach to detect genomic regions enriched for autozygosity signals in one group relative to another based upon comparisons of per-individual autozygosity likelihoods instead of inferred ROA frequencies.

Conclusions

This weighted likelihood ROA inference approach can assist population- and disease-geneticists working with a wide variety of data types and species to explore ROA patterns and to identify genomic regions with differential ROA signals among groups, thereby advancing our understanding of evolutionary history and the role of recessive variation in phenotypic variation and disease.

Electronic supplementary material

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