Novel genomic signals of recent selection in an Ethiopian population

Understanding the contemporary high obesity rate from an evolutionary genetic perspective

The topic of obesity is gaining increasing popularity globally. From an evolutionary genetic perspective, it is believed that the main cause of the high obesity rate is the mismatch between environment and genes after people have shifted toward a modern high-calorie diet. However, it has been debated for over 60 years about how obesity-related genes become prevalent all over the world. Here, we review the three most influential hypotheses or viewpoints, i.e., the thrifty gene hypothesis, the drifty gene hypothesis, and the maladaptation viewpoint. In particular, genome-wide association studies in the recent 10 years have provided rich findings and evidence to be considered for a better understanding of the evolutionary genetic mechanisms of obesity. We anticipate this brief review to direct further studies and inspire the future application of precision medicine in obesity treatment.

Including diverse and admixed populations in genetic epidemiology research

The inclusion of ancestrally diverse participants in genetic studies can lead to new discoveries and is important to ensure equitable health care benefit from research advances. Here, members of the Ethical, Legal, Social, Implications (ELSI) committee of the International Genetic Epidemiology Society (IGES) offer perspectives on methods and analysis tools for the conduct of inclusive genetic epidemiology research, with a focus on admixed and ancestrally diverse populations in support of reproducible research practices. We emphasize the importance of distinguishing socially defined population categorizations from genetic ancestry in the design, analysis, reporting, and interpretation of genetic epidemiology research findings. Finally, we discuss the current state of genomic resources used in genetic association studies, functional interpretation, and clinical and public health translation of genomic findings with respect to diverse populations.

Detecting Genetic Ancestry and Adaptation in the Taiwanese Han People

The Taiwanese people are composed of diverse indigenous populations and the Taiwanese Han. About 95% of the Taiwanese identify themselves as Taiwanese Han, but this may not be a homogeneous population because they migrated to the island from various regions of continental East Asia over a period of 400 years. Little is known about the underlying patterns of genetic ancestry, population admixture, and evolutionary adaptation in the Taiwanese Han people. Here, we analyzed the whole-genome single-nucleotide polymorphism genotyping data from 14,401 individuals of Taiwanese Han collected by the Taiwan Biobank and the whole-genome sequencing data for a subset of 772 people. We detected four major genetic ancestries with distinct geographic distributions (i.e., Northern, Southeastern, Japonic, and Island Southeast Asian ancestries) and signatures of population mixture contributing to the genomes of Taiwanese Han. We further scanned for signatures of positive natural selection that caused unusually long-range haplotypes and elevations of hitchhiked variants. As a result, we identified 16 candidate loci in which selection signals can be unambiguously localized at five single genes: CTNNA2, LRP1B, CSNK1G3, ASTN2, and NEO1. Statistical associations were examined in 16 metabolic-related traits to further elucidate the functional effects of each candidate gene. All five genes appear to have pleiotropic connections to various types of disease susceptibility and significant associations with at least one metabolic-related trait. Together, our results provide critical insights for understanding the evolutionary history and adaption of the Taiwanese Han population.

The evolution of skin pigmentation-associated variation in West Eurasia

Some of the genes responsible for the evolution of light skin pigmentation in Europeans show signals of positive selection in present-day populations. Recently, genome-wide association studies have highlighted the highly polygenic nature of skin pigmentation. It is unclear whether selection has operated on all of these genetic variants or just a subset. By studying variation in over a thousand ancient genomes from West Eurasia covering 40,000 y, we are able to study both the aggregate behavior of pigmentation-associated variants and the evolutionary history of individual variants. We find that the evolution of light skin pigmentation in Europeans was driven by frequency changes in a relatively small fraction of the genetic variants that are associated with variation in the trait today.

Skin pigmentation is a classic example of a polygenic trait that has experienced directional selection in humans. Genome-wide association studies have identified well over a hundred pigmentation-associated loci, and genomic scans in present-day and ancient populations have identified selective sweeps for a small number of light pigmentation-associated alleles in Europeans. It is unclear whether selection has operated on all of the genetic variation associated with skin pigmentation as opposed to just a small number of large-effect variants. Here, we address this question using ancient DNA from 1,158 individuals from West Eurasia covering a period of 40,000 y combined with genome-wide association summary statistics from the UK Biobank. We find a robust signal of directional selection in ancient West Eurasians on 170 skin pigmentation-associated variants ascertained in the UK Biobank. However, we also show that this signal is driven by a limited number of large-effect variants. Consistent with this observation, we find that a polygenic selection test in present-day populations fails to detect selection with the full set of variants. Our data allow us to disentangle the effects of admixture and selection. Most notably, a large-effect variant at SLC24A5 was introduced to Western Europe by migrations of Neolithic farming populations but continued to be under selection post-admixture. This study shows that the response to selection for light skin pigmentation in West Eurasia was driven by a relatively small proportion of the variants that are associated with present-day phenotypic variation.

Evolutionary forces in diabetes and hypertension pathogenesis in Africans

Rates of type 2 diabetes (T2D) and hypertension are increasing rapidly in urbanizing sub-Saharan Africa (SSA). While lifestyle factors drive the increases in T2D and hypertension prevalence, evidence across populations shows that genetic variation, which is driven by evolutionary forces including a natural selection that shaped the human genome, also plays a role. Here we report the evidence for the effect of selection in African genomes on mechanisms underlying T2D and hypertension, including energy metabolism, adipose tissue biology, insulin action and salt retention. Selection effects found for variants in genes PPARA and TCF7L2 may have enabled Africans to respond to nutritional challenges by altering carbohydrate and lipid metabolism. Likewise, African-ancestry-specific characteristics of adipose tissue biology (low visceral adipose tissue [VAT], high intermuscular adipose tissue and a strong association between VAT and adiponectin) may have been selected for in response to nutritional and infectious disease challenges in the African environment. Evidence for selection effects on insulin action, including insulin resistance and secretion, has been found for several genes including MPHOSPH9, TMEM127, ZRANB3 and MC3R. These effects may have been historically adaptive in critical conditions, such as famine and inflammation. A strong correlation between hypertension susceptibility variants and latitude supports the hypothesis of selection for salt retention mechanisms in warm, humid climates. Nevertheless, adaptive genomics studies in African populations are scarce. More work is needed, particularly genomics studies covering the wide diversity of African populations in SSA and Africans in diaspora, as well as further functional assessment of established risk loci.

Genomic landscape of the signals of positive natural selection in populations of Northern Eurasia: A view from Northern Russia

Natural selection of beneficial genetic variants played a critical role in human adaptation to a wide range of environmental conditions. Northern Eurasia, despite its severe climate, is home to lots of ethnically diverse populations. The genetic variants associated with the survival of these populations have hardly been analyzed. We searched for the genomic signatures of positive selection in (1) the genome-wide microarray data of 432 people from eight different northern Russian populations and (2) the whole-genome sequences of 250 people from Northern Eurasia from a public repository through testing the extended haplotype homozigosity (EHH) and direct comparison of allele frequency, respectively. The 20 loci with the strongest selection signals were characterized in detail. Among the top EHH hits were the NRG3 and NBEA genes, which are involved in the development and functioning of the neural system, the PTPRM gene, which mediates cell-cell interactions and adhesion, and a region on chromosome 4 (chr4:28.7-28.9 Mb) that contained several loci affiliated with different classes of non-coding RNAs (RN7SL101P, MIR4275, MESTP3, and LINC02364). NBEA and the region on chromosome 4 were novel selection targets that were identified for the first time in Western Siberian populations. Cross-population comparisons of EHH profiles suggested a particular role for the chr4:28.7-28.9 Mb region in the local adaptation of Western Siberians. The strongest selection signal identified in Siberian sequenced genomes was formed by six SNPs on chromosome 11 (chr11:124.9-125.2 Mb). This region included well-known genes SLC37A2 and PKNOX2. SLC37A2 is most-highly expressed in the gut. Its expression is regulated by vitamin D, which is often deficient in northern regions. The PKNOX2 gene is a transcription factor of the homeobox family that is expressed in the brain and many other tissues. This gene is associated with alcohol addiction, which is widespread in many Northern Eurasian populations.

Whole-Genome Sequencing Identifies the Egl Nine Homologue 3 (egln3/phd3) and Protein Phosphatase 1 Regulatory Inhibitor Subunit 2 (PPP1R2P1) Associated with High-Altitude Polycythemia in Tibetans at High Altitude

Background

The hypoxic conditions at high altitudes are great threats to survival, causing pressure for adaptation. More and more high-altitude denizens are not adapted with the condition known as high-altitude polycythemia (HAPC) that featured excessive erythrocytosis. As a high-altitude sickness, the etiology of HAPC is still unclear.

Methods

In this study, we reported the whole-genome sequencing-based study of 10 native Tibetans with HAPC and 10 control subjects followed by genotyping of selected 21 variants from discovered single nucleotide variants (SNVs) in an independent cohort (232 cases and 266 controls).

Results

We discovered the egl nine homologue 3 (egln3/phd3) (14q13.1, rs1346902, P = 1.91 × 10⁻⁵) and PPP1R2P1 (Protein Phosphatase 1 Regulatory Inhibitor Subunit 2) gene (6p21.32, rs521539, P = 0.012). Our results indicated an unbiased framework to identify etiological mechanisms of HAPC and showed that egln3/phd3 and PPP1R2P1 may be associated with the susceptibility to HAPC. Egln3/phd3b is associated with hypoxia-inducible factor subunit α (HIFα). Protein Phosphatase 1 Regulatory Inhibitor is associated with reactive oxygen species (ROS) and oxidative stress.

Conclusions

Our genome sequencing conducted in Tibetan HAPC patients identified egln3/phd3 and PPP1R2P1 associated with HAPC.

Genome-Wide Characterization of Arabian Peninsula Populations: Shedding Light on the History of a Fundamental Bridge between Continents

The Arabian Peninsula (AP) was an important crossroad between Africa, Asia, and Europe, being the cradle of the structure defining these main human population groups, and a continuing path for their admixture. The screening of 741,000 variants in 420 Arabians and 80 Iranians allowed us to quantify the dominant sub-Saharan African admixture in the west of the peninsula, whereas South Asian and Levantine/European influence was stronger in the east, leading to a rift between western and eastern sides of the Peninsula. Dating of the admixture events indicated that Indian Ocean slave trade and Islamization periods were important moments in the genetic makeup of the region. The western-eastern axis was also observable in terms of positive selection of diversity conferring lactose tolerance, with the West AP developing local adaptation and the East AP acquiring the derived allele selected in European populations and existing in South Asia. African selected malaria resistance through the DARC gene was enriched in all Arabian genomes, especially in the western part. Clear European influences associated with skin and eye color were equally frequent across the Peninsula.

Combination of 247 Genome-Wide Association Studies Reveals High Cancer Risk as a Result of Evolutionary Adaptation

Analysis of GLOBOCAN-2012 data shows clearly here that cancer incidence worldwide is highly related with low average annual temperatures and extreme low temperatures. This applies for all cancers together or separately for many frequent or rare cancer types (all cancers P = 9.49×10-18). Supporting fact is that Inuit people, living at extreme low temperatures, have the highest cancer rates today. Hypothesizing an evolutionary explanation, 240 cancer genome-wide association studies, and seven genome-wide association studies for cold and high-altitude adaptation were combined. A list of 1,377 cancer-associated genes was created to initially investigate whether cold selected genes are enriched with cancer-associated genes. Among Native Americans, Inuit and Eskimos, the highest association was observed for Native Americans (P = 6.7×10-5). An overall or a meta-analysis approach confirmed further this result. Similar approach for three populations living at extreme high altitude, revealed high association for Andeans-Tibetans (P = 1.3×10-11). Overall analysis or a meta-analysis was also significant. A separate analysis showed special selection for tumor suppressor genes. These results can be viewed along with those of previous functional studies that showed that reduced apoptosis potential due to specific p53 variants (the most important tumor suppressor gene) is beneficial in high-altitude and cold environments. In conclusion, this study shows that genetic variants selected for adaptation at extreme environmental conditions can increase cancer risk later on age. This is in accordance with antagonistic pleiotropy hypothesis.

The roles of vitamin D and cutaneous vitamin D production in human evolution and health

Most of the vitamin D necessary for the maintenance of human health and successful reproduction is made in the skin under the influence of a narrow portion of the electromagnetic spectrum emitted from the sun, namely ultraviolet B radiation (UVB). During the course of human evolution, skin pigmentation has evolved to afford protection against high levels of UVR while still permitting cutaneous production of vitamin D. Similar pigmentation phenotypes evolved repeatedly as the result of independent genetic events when isolated human populations dispersed into habitats of extremely low or high UVB. The gradient of skin color seen in modern human populations is evidence of the operation of two clines, one favoring photoprotection near the equator, the other favoring vitamin D production nearer the poles. Through time, human adaptations to different solar regimes have become more cultural than biological. Rapid human migrations, increasing urbanization, and changes in lifestyle have created mismatches between skin pigmentation and environmental conditions leading to vitamin D deficiency. The prevalence and significance for health of vitamin D deficiencies, and the definition of optimal levels of vitamin D in the bloodstream are subjects of intense research and debate, but two of the causes of vitamin D deficiency - lack of sun exposure and abandonment of vitamin D rich foods in the diet - are traceable to changes in human lifestyles accompanying urbanization in prehistory.

Homo sapiens-Specific Binding Site Variants within Brain Exclusive Enhancers Are Subject to Accelerated Divergence across Human Population

Empirical assessments of human accelerated noncoding DNA frgaments have delineated presence of many cis-regulatory elements. Enhancers make up an important category of such accelerated cis-regulatory elements that efficiently control the spatiotemporal expression of many developmental genes. Establishing plausible reasons for accelerated enhancer sequence divergence in Homo sapiens has been termed significant in various previously published studies. This acceleration by including closely related primates and archaic human data has the potential to open up evolutionary avenues for deducing present-day brain structure. This study relied on empirically confirmed brain exclusive enhancers to avoid any misjudgments about their regulatory status and categorized among them a subset of enhancers with an exceptionally accelerated rate of lineage specific divergence in humans. In this assorted set, 13 distinct transcription factor binding sites were located that possessed unique existence in humans. Three of 13 such sites belonging to transcription factors SOX2, RUNX1/3, and FOS/JUND possessed single nucleotide variants that made them unique to H. sapiens upon comparisons with Neandertal and Denisovan orthologous sequences. These variants modifying the binding sites in modern human lineage were further substantiated as single nucleotide polymorphisms via exploiting 1000 Genomes Project Phase3 data. Long range haplotype based tests laid out evidence of positive selection to be governing in African population on two of the modern human motif modifying alleles with strongest results for SOX2 binding site. In sum, our study acknowledges acceleration in noncoding regulatory landscape of the genome and highlights functional parts within it to have undergone accelerated divergence in present-day human population.

Loci associated with skin pigmentation identified in African populations

Despite the wide range of skin pigmentation in humans, little is known about its genetic basis in global populations. Examining ethnically diverse African genomes, we identify variants in or near SLC24A5, MFSD12, DDB1, TMEM138, OCA2, and HERC2 that are significantly associated with skin pigmentation. Genetic evidence indicates that the light pigmentation variant at SLC24A5 was introduced into East Africa by gene flow from non-Africans. At all other loci, variants associated with dark pigmentation in Africans are identical by descent in South Asian and Australo-Melanesian populations. Functional analyses indicate that MFSD12 encodes a lysosomal protein that affects melanogenesis in zebrafish and mice, and that mutations in melanocyte-specific regulatory regions near DDB1/TMEM138 correlate with expression of ultraviolet response genes under selection in Eurasians.

High-altitude adaptation in humans: from genomics to integrative physiology

About 1.2 to 33% of high-altitude populations suffer from Monge's disease or chronic mountain sickness (CMS). Number of factors such as age, sex, and population of origin (older, male, Andean) contribute to the percentage reported from a variety of samples. It is estimated that there are around 83 million people who live at altitudes > 2500 m worldwide and are at risk for CMS. In this review, we focus on a human "experiment in nature" in various high-altitude locations in the world-namely, Andean, Tibetan, and Ethiopian populations that have lived under chronic hypoxia conditions for thousands of years. We discuss the adaptive as well as mal-adaptive changes at the genomic and physiological levels. Although different genes seem to be involved in adaptation in the three populations, we can observe convergence at genetic and signaling, as well as physiological levels. What is important is that we and others have shown that lessons learned from the genes mined at high altitude can be helpful in better understanding and treating diseases that occur at sea level. We discuss two such examples: EDNRB and SENP1 and their role in cardiac tolerance and in the polycythemic response, respectively.

Diversity and inclusion in genomic research: why the uneven progress?

Conducting genomic research in diverse populations has led to numerous advances in our understanding of human history, biology, and health disparities, in addition to discoveries of vital clinical significance. Conducting genomic research in diverse populations is also important in ensuring that the genomic revolution does not exacerbate health disparities by facilitating discoveries that will disproportionately benefit well-represented populations. Despite the general agreement on the need for genomic research in diverse populations in terms of equity and scientific progress, genomic research remains largely focused on populations of European descent. In this article, we describe the rationale for conducting genomic research in diverse populations by reviewing examples of advances facilitated by their inclusion. We also explore some of the factors that perpetuate the disproportionate attention on well-represented populations. Finally, we discuss ongoing efforts to ameliorate this continuing bias. Collaborative and intensive efforts at all levels of research, from the funding of studies to the publication of their findings, will be necessary to ensure that genomic research does not conserve historical inequalities or curtail the contribution that genomics could make to the health of all humanity.

Examining How Our Shared Evolutionary History Shapes Future Disease Outcomes

Cardiometabolic diseases are major contributors to mortality and morbidity, and their burden displays global and regional disparities. Gene-environment interactions contribute to the pathogenesis of cardiometabolic diseases. Population differences in genetic structure, ancient environmental pressures that shape the human genome, and early life environmental adversities (e.g., in utero conditions) all contribute to observed disparities in global cardiometabolic diseases. The genetic and sociocultural diversity of global populations presents opportunities for discovering genomic loci that influence cardiometabolic diseases as illustrated by a few genetic, epigenetic, and population-genetic discoveries leading to notable understanding of disease mechanisms. However, African, Latin American and Hispanic, and indigenous peoples represent <4% of all genome-wide association study samples analyzed to date. Using examples of recent studies in African populations, we discuss the crucial importance of conducting genomic studies in ancestrally diverse populations to understand disease mechanisms and to prepare fertile ground for future delivery of precise health care to all individuals.

The African diaspora: history, adaptation and health

The trans-Atlantic slave trade brought millions of Africans to the New World. Advances in genomics are providing novel insights into the history and health of Africans and the diasporan populations. Recent examples reviewed here include the unraveling of substantial hunter-gatherer and 'Eurasian' admixtures across sub-Saharan Africa, expanding our understanding of ancestral African genetics; the global ubiquity of mixed ancestry; the revealing of African ancestry in Latin Americans that likely derived from the slave trade; and understanding of the ancestral backgrounds of APOL1 and LPL found to influence kidney disease and lipid levels, respectively, providing specific insights into disease etiology and health disparities.

A probabilistic method for testing and estimating selection differences between populations

Human populations around the world encounter various environmental challenges and, consequently, develop genetic adaptations to different selection forces. Identifying the differences in natural selection between populations is critical for understanding the roles of specific genetic variants in evolutionary adaptation. Although numerous methods have been developed to detect genetic loci under recent directional selection, a probabilistic solution for testing and quantifying selection differences between populations is lacking. Here we report the development of a probabilistic method for testing and estimating selection differences between populations. By use of a probabilistic model of genetic drift and selection, we showed that logarithm odds ratios of allele frequencies provide estimates of the differences in selection coefficients between populations. The estimates approximate a normal distribution, and variance can be estimated using genome-wide variants. This allows us to quantify differences in selection coefficients and to determine the confidence intervals of the estimate. Our work also revealed the link between genetic association testing and hypothesis testing of selection differences. It therefore supplies a solution for hypothesis testing of selection differences. This method was applied to a genome-wide data analysis of Han and Tibetan populations. The results confirmed that both the EPAS1 and EGLN1 genes are under statistically different selection in Han and Tibetan populations. We further estimated differences in the selection coefficients for genetic variants involved in melanin formation and determined their confidence intervals between continental population groups. Application of the method to empirical data demonstrated the outstanding capability of this novel approach for testing and quantifying differences in natural selection.

Predicting Carriers of Ongoing Selective Sweeps without Knowledge of the Favored Allele

Methods for detecting the genomic signatures of natural selection have been heavily studied, and they have been successful in identifying many selective sweeps. For most of these sweeps, the favored allele remains unknown, making it difficult to distinguish carriers of the sweep from non-carriers. In an ongoing selective sweep, carriers of the favored allele are likely to contain a future most recent common ancestor. Therefore, identifying them may prove useful in predicting the evolutionary trajectory—for example, in contexts involving drug-resistant pathogen strains or cancer subclones. The main contribution of this paper is the development and analysis of a new statistic, the Haplotype Allele Frequency (HAF) score. The HAF score, assigned to individual haplotypes in a sample, naturally captures many of the properties shared by haplotypes carrying a favored allele. We provide a theoretical framework for computing expected HAF scores under different evolutionary scenarios, and we validate the theoretical predictions with simulations. As an application of HAF score computations, we develop an algorithm (PreCIOSS: Predicting Carriers of Ongoing Selective Sweeps) to identify carriers of the favored allele in selective sweeps, and we demonstrate its power on simulations of both hard and soft sweeps, as well as on data from well-known sweeps in human populations.

Methods for detecting the genomic signatures of natural selection have been heavily studied, and they have been successful in identifying genomic regions under positive selection. However, methods that detect positive selective sweeps do not typically identify the favored allele, or even the haplotypes carrying the favored allele. The main contribution of this paper is the development and analysis of a new statistic (the HAF score), assigned to individual haplotypes. Using both theoretical analyses and simulations, we describe how the HAF scores differ for carriers and non-carriers of the favored allele, and how they change dynamically during a selective sweep. We also develop an algorithm, PreCIOSS, for separating carriers and non-carriers. Our tool has broad applicability as carriers of the favored allele are likely to contain a future most recent common ancestor. Therefore, identifying them may prove useful in predicting the evolutionary trajectory—for example, in contexts involving drug-resistant pathogen strains or cancer subclones.