Genetic signature of natural selection in first Americans

Ancient genomes reveal insights into ritual life at Chichén Itzá

The ancient city of Chichén Itzá in Yucatán, Mexico, was one of the largest and most influential Maya settlements during the Late and Terminal Classic periods (AD 600-1000) and it remains one of the most intensively studied archaeological sites in Mesoamerica1-4. However, many questions about the social and cultural use of its ceremonial spaces, as well as its population's genetic ties to other Mesoamerican groups, remain unanswered2. Here we present genome-wide data obtained from 64 subadult individuals dating to around AD 500-900 that were found in a subterranean mass burial near the Sacred Cenote (sinkhole) in the ceremonial centre of Chichén Itzá. Genetic analyses showed that all analysed individuals were male and several individuals were closely related, including two pairs of monozygotic twins. Twins feature prominently in Mayan and broader Mesoamerican mythology, where they embody qualities of duality among deities and heroes5, but until now they had not been identified in ancient Mayan mortuary contexts. Genetic comparison to present-day people in the region shows genetic continuity with the ancient inhabitants of Chichén Itzá, except at certain genetic loci related to human immunity, including the human leukocyte antigen complex, suggesting signals of adaptation due to infectious diseases introduced to the region during the colonial period.

Sex-biased gene regulation varies across human populations as a result of adaptive evolution

OBJECTIVES

Studies of human sexual dimorphism and gender disparities in health focus on ostensibly universal molecular sex differences, such as sex chromosomes and circulating hormone levels, while ignoring the extraordinary diversity in biology, behavior, and culture acquired by different human populations over their unique evolutionary histories.

MATERIALS AND METHODS

Using RNA-Seq data and whole genome sequences from 1000G and HGDP, we investigate variation in sex-biased gene expression across 11 human populations and test whether population-level variation in sex-biased expression may have resulted from adaptive evolution in regions containing sex-specific regulatory variants.

RESULTS

We find that sex-biased gene expression in humans is highly variable, mostly population-specific, and demonstrates between population reversals. Expression quantitative trait locus mapping reveals sex-specific regulatory regions with evidence of recent positive natural selection, suggesting that variation in sex-biased expression may have evolved as an adaptive response to ancestral environments experienced by human populations.

DISCUSSION

These results indicate that sex-biased gene expression is more flexible than previously thought and is not generally shared among human populations. Instead, molecular phenotypes associated with sex depend on complex interactions between population-specific molecular evolution and physiological responses to contemporary socioecologies.

Human adaptations to diet: Biological and cultural coevolution

Modern humans evolved in Africa some 200,000 years ago, and since then, human populations have expanded and diversified to occupy a broad range of habitats and use different subsistence modes. This has resulted in different adaptations, such as differential responses to diseases and different abilities to digest or tolerate certain foods. The shift from a subsistence strategy based on hunting and gathering during the Palaeolithic to a lifestyle based on the consumption of domesticated animals and plants in the Neolithic can be considered one of the most important dietary transitions of Homo sapiens. In this text, we review four examples of gene-culture coevolution: (i) the persistence of the enzyme lactase after weaning, which allows the digestion of milk in adulthood, related to the emergence of dairy farming during the Neolithic; (ii) the population differences in alcohol susceptibility, in particular the ethanol intolerance of Asian populations due to the increased accumulation of the toxic acetaldehyde, related to the spread of rice domestication; (iii) the maintenance of gluten intolerance (celiac disease) with the subsequent reduced fitness of its sufferers, related to the emergence of agriculture and (iv) the considerable variation in the biosynthetic pathway of long-chain polyunsaturated fatty acids in native populations with extreme diets.

Peopling of the Americas: A new approach to assessing dental morphological variation in Asian and Native American populations

OBJECTIVES

Through biodistance analyses, anthropologists have used dental morphology to elucidate how people moved into and throughout the Americas. Here, we apply a method that focuses on individuals rather than sample frequencies through the application rASUDAS2, based on a naïve Bayes' algorithm.

MATERIALS AND METHODS

Using the database of C.G. Turner II, we calculated the probability that an individual could be assigned to one of seven biogeographic groups (American Arctic, North & South America, East Asia, Southeast Asia & Polynesia, Australo-Melanesia, Western Eurasia, & Sub-Saharan Africa) through rASUDAS2. The frequency of classifications for each biogeographic group was determined for 1418 individuals from six regions across Asia and the Americas.

RESULTS

Southeast Asians show mixed assignments but rarely to American Arctic or "American Indian." East Asians are assigned to East Asia half the time while 30% are assigned as Native American. People from the American Arctic and North & South America are assigned to Arctic America or non-Arctic America 75%-80% of the time, with 10%-15% classified as East Asian.

DISCUSSION

All Native American groups have a similar degree of morphological affinity to East Asia, as 10%-15% are classified as East Asian. East Asians are classified as Native American in 30% of cases. Individuals in the Western Hemisphere are decreasingly classified as Arctic the farther south they are located. Equivalent levels of classification as East Asian across all Native American groups suggests one divergence between East Asians and the population ancestral to all Native Americans. Non-arctic Native American groups are derived from the Arctic population, which represents the Native American founder group.

Ancestry-driven metabolite variation provides insights into disease states in admixed populations

BACKGROUND

Metabolic pathways are related to physiological functions and disease states and are influenced by genetic variation and environmental factors. Hispanics/Latino individuals have ancestry-derived genomic regions (local ancestry) from their recent admixture that have been less characterized for associations with metabolite abundance and disease risk.

METHODS

We performed admixture mapping of 640 circulating metabolites in 3887 Hispanic/Latino individuals from the Hispanic Community Health Study/Study of Latinos (HCHS/SOL). Metabolites were quantified in fasting serum through non-targeted mass spectrometry (MS) analysis using ultra-performance liquid chromatography-MS/MS. Replication was performed in 1856 nonoverlapping HCHS/SOL participants with metabolomic data.

RESULTS

By leveraging local ancestry, this study identified significant ancestry-enriched associations for 78 circulating metabolites at 484 independent regions, including 116 novel metabolite-genomic region associations that replicated in an independent sample. Among the main findings, we identified Native American enriched genomic regions at chromosomes 11 and 15, mapping to FADS1/FADS2 and LIPC, respectively, associated with reduced long-chain polyunsaturated fatty acid metabolites implicated in metabolic and inflammatory pathways. An African-derived genomic region at chromosome 2 was associated with N-acetylated amino acid metabolites. This region, mapped to ALMS1, is associated with chronic kidney disease, a disease that disproportionately burdens individuals of African descent.

CONCLUSIONS

Our findings provide important insights into differences in metabolite quantities related to ancestry in admixed populations including metabolites related to regulation of lipid polyunsaturated fatty acids and N-acetylated amino acids, which may have implications for common diseases in populations.

The role of Beringia in human adaptation to Arctic conditions based on results of genomic studies of modern and ancient populations

The results of studies in Quaternary geology, archeology, paleoanthropology and human genetics demonstrate that the ancestors of Native Americans arrived in mid-latitude North America mainly along the Pacific Northwest Coast, but had previously inhabited the Arctic and during the last glacial maximum were in a refugium in Beringia, a land bridge connecting Eurasia and North America. The gene pool of Native Americans is represented by unique haplogroups of mitochondrial DNA and the Y chromosome, the evolutionary age of which ranges from 13 to 22 thousand years. The results of a paleogenomic analysis also show that during the last glacial maximum Beringia was populated by human groups that had arisen as a result of interaction between the most ancient Upper Paleolithic populations of Northern Eurasia and newcomer groups from East Asia. Approximately 20 thousand years ago the Beringian populations began to form, and the duration of their existence in relative isolation is estimated at about 5 thousand years. Thus, the adaptation of the Beringians to the Arctic conditions could have taken several millennia. The adaptation of Amerindian ancestors to high latitudes and cold climates is supported by genomic data showing that adaptive genetic variants in Native Americans are associated with various metabolic pathways: melanin production processes in the skin, hair and eyes, the functioning of the cardiovascular system, energy metabolism and immune response characteristics. Meanwhile, the analysis of the existing hypotheses about the selection of some genetic variants in the Beringian ancestors of the Amerindians in connection with adaptation to the Arctic conditions (for example, in the FADS, ACTN3, EDAR genes) shows the ambiguity of the testing results, which may be due to the loss of some traces of the "Beringian" adaptation in the gene pools of modern Native Americans. The most optimal strategy for further research seems to be the search for adaptive variant.

Dietary oily fish intake reduces the risk of all-cause mortality in frequent fish consumers of Amerindian ancestry living in coastal Ecuador: the Atahualpa project

PURPOSE

To assess the relationship between dietary oily fish intake and all-cause mortality in a population of frequent fish consumers of Amerindian ancestry living in rural Ecuador.

METHODS

Individuals aged ≥ 40 years enrolled in the prospective population-based Atahualpa Project cohort received annual questionnaires to estimate their dietary oily fish intake. Only fish served broiled or cooked in the soup were included for analysis. Poisson regression and Cox-proportional hazards models adjusted for demographics, education level and cardiovascular risk factors were obtained to estimate mortality risk according to the amount of oily fish intake stratified in tertiles.

RESULTS

Analysis included 909 individuals (mean age: 55.1 ± 12.8 years) followed by a median of 7.5 ± 3 years. Mean oily fish intake was 9.4 ± 5.7 servings per week. A total of 142 (16%) individuals died during the follow-up. The mortality rate for individuals in the first tertile de oily fish intake (0.0-6.29 servings) was 2.87 per 100 person-years, which decreased to 1.78 for those in the third tertile (10.59-35.0 servings). An adjusted Cox-proportional hazards model showed that individuals allocated to the second (HR 0.61; 95% CI 0.41-0.92) and third (HR 0.60; 95% CI 0.40-0.91) tertiles of dietary oily fish intake had significantly lower mortality risk than those in the first tertile.

CONCLUSION

Sustained oily fish intake of more than six servings per week reduces mortality risk in middle-aged and older adults of Amerindian ancestry.

Indigenous people from Amazon show genetic signatures of pathogen-driven selection

Ecological conditions in the Amazon rainforests are historically favorable for the transmission of numerous tropical diseases, especially vector-borne diseases. The high diversity of pathogens likely contributes to the strong selective pressures for human survival and reproduction in this region. However, the genetic basis of human adaptation to this complex ecosystem remains unclear. This study investigates the possible footprints of genetic adaptation to the Amazon rainforest environment by analyzing the genomic data of 19 native populations. The results based on genomic and functional analysis showed an intense signal of natural selection in a set of genes related to Trypanosoma cruzi infection, which is the pathogen responsible for Chagas disease, a neglected tropical parasitic disease native to the Americas that is currently spreading worldwide.

Blocks identical by descent in the genomes of the indigenous population of Siberia demonstrate genetic links between populations

The gene pool of the indigenous population of Siberia is a unique system for studying population and evolutionary genetic processes, analyzing genetic diversity, and reconstructing the genetic history of populations. High ethnic diversity is a feature of Siberia, as one of the regions of the peripheral settlement of modern human. The vast expanses of this region and the small number of aboriginal populations contributed to the formation of significant territorial and genetic subdivision. About 40 indigenous peoples are settled on the territory of the Siberian historical and ethnographic province. Within the framework of this work, a large-scale population study of the gene pool of the indigenous peoples of Siberia was carried out for the first time at the level of high-density biochips. This makes it possible to fill in a significant gap in the genogeographic picture of the Eurasian population. For this, DNA fragments were analyzed, which had been inherited without recombination by each pair of individuals from their recent common ancestor, that is, segments (blocks) identical by descent (IBD). The distribution of IBD blocks in the populations of Siberia is in good agreement with the geographical proximity of the populations and their linguistic affiliation. Among the Siberian populations, the Chukchi, Koryaks, and Nivkhs form a separate cluster from the main Siberian group, with the Chukchi and Koryaks being more closely related. Separate subclusters of Evenks and Yakuts, Kets and Chulyms are formed within the Siberian cluster. Analysis of SNPs that fell into more IBD segments of the analyzed populations made it possible to compile a list of 5358 genes. According to the calculation results, biological processes enriched with these genes are associated with the detection of a chemical stimulus involved in the sensory perception of smell. Enriched for the genes found, molecular pathways are associated with the metabolism of linoleic, arachidonic, tyrosic acids and by olfactory transduction. At the same time, an analysis of the literature data showed that some of the selected genes, which were found in a larger number of IBD blocks in several populations at once, can play a role in genetic adaptation to environmental factors.

Beringia and the peopling of the Western Hemisphere

Did Beringian environments represent an ecological barrier to humans until less than 15 000 years ago or was access to the Americas controlled by the spatial-temporal distribution of North American ice sheets? Beringian environments varied with respect to climate and biota, especially in the two major areas of exposed continental shelf. The East Siberian Arctic Shelf ('Great Arctic Plain' (GAP)) supported a dry steppe-tundra biome inhabited by a diverse large-mammal community, while the southern Bering-Chukchi Platform ('Bering Land Bridge' (BLB)) supported mesic tundra and probably a lower large-mammal biomass. A human population with west Eurasian roots occupied the GAP before the Last Glacial Maximum (LGM) and may have accessed mid-latitude North America via an interior ice-free corridor. Re-opening of the corridor less than 14 000 years ago indicates that the primary ancestors of living First Peoples, who already had spread widely in the Americas at this time, probably dispersed from the NW Pacific coast. A genetic 'arctic signal' in non-arctic First Peoples suggests that their parent population inhabited the GAP during the LGM, before their split from the former. We infer a shift from GAP terrestrial to a subarctic maritime economy on the southern BLB coast before dispersal in the Americas from the NW Pacific coast.

Unraveling Signatures of Local Adaptation among Indigenous Groups from Mexico

Few studies have addressed how selective pressures have shaped the genetic structure of the current Native American populations, and they have mostly limited their inferences to admixed Latin American populations. Here, we searched for local adaptation signals, based on integrated haplotype scores and population branch statistics, in 325 Mexican Indigenous individuals with at least 99% Native American ancestry from five previously defined geographical regions. Although each region exhibited its own local adaptation profile, only PPARG and AJAP1, both negative regulators of the Wnt/β catenin signaling pathway, showed significant adaptation signals in all the tested regions. Several signals were found, mainly in the genes related to the metabolic processes and immune response. A pathway enrichment analysis revealed the overrepresentation of selected genes related to several biological phenotypes/conditions, such as the immune response and metabolic pathways, in agreement with previous studies, suggesting that immunological and metabolic pressures are major drivers of human adaptation. Genes related to the gut microbiome measurements were overrepresented in all the regions, highlighting the importance of studying how humans have coevolved with the microbial communities that colonize them. Our results provide a further explanation of the human evolutionary history in response to environmental pressures in this region.

Repeated genetic adaptation to altitude in two tropical butterflies

Repeated evolution can provide insight into the mechanisms that facilitate adaptation to novel or changing environments. Here we study adaptation to altitude in two tropical butterflies, Heliconius erato and H. melpomene, which have repeatedly and independently adapted to montane habitats on either side of the Andes. We sequenced 518 whole genomes from altitudinal transects and found many regions differentiated between highland (~ 1200 m) and lowland (~ 200 m) populations. We show repeated genetic differentiation across replicate populations within species, including allopatric comparisons. In contrast, there is little molecular parallelism between the two species. By sampling five close relatives, we find that a large proportion of divergent regions identified within species have arisen from standing variation and putative adaptive introgression from high-altitude specialist species. Taken together our study supports a role for both standing genetic variation and gene flow from independently adapted species in promoting parallel local adaptation to the environment.

Here, the authors study adaptation to altitude in 518 whole genomes from two species of tropical butterflies. They find repeated genetic differentiation within species, little molecular parallelism between these species, and introgression from closely related species, concluding that standing genetic variation promotes parallel local adaptation.

Uncovering signals of positive selection in Peruvian populations from three ecological regions

Perú hosts extremely diverse ecosystems which can be broadly classified into three major ecoregions: the Pacific desert coast, the Andean highlands, and the Amazon rainforest. Since its initial peopling approximately 12,000 years ago, the populations inhabiting such ecoregions might have differentially adapted to their contrasting environmental pressures. Previous studies have described several candidate genes underlying adaptation to hypobaric hypoxia among Andean highlanders. However, the adaptive genetic diversity of coastal and rainforest populations has been less studied. Here, we gathered genome-wide SNP-array data from 286 Peruvians living across the three ecoregions and analysed signals of recent positive selection through population differentiation and haplotype-based selection scans. Among highland populations, we identify candidate genes related to cardiovascular function (TLL1, DUSP27, TBX5, PLXNA4, SGCD), to the Hypoxia-Inducible Factor pathway (TGFA, APIP), to skin pigmentation (MITF), as well as to glucose (GLIS3) and glycogen metabolism (PPP1R3C, GANC). In contrast, most signatures of adaptation in coastal and rainforest populations comprise candidate genes related to the immune system (including SIGLEC8, TRIM21, CD44 and ICAM1 in the coast; CBLB and PRDM1 in rainforest and the BRD2- HLA-DOA- HLA-DPA1 region in both), possibly as a result of strong pathogen-driven selection. This study identifies candidate genes related to human adaptation to the diverse environments of South America.

Evolutionary adaptation highlights the interconnection of fatty acids, sunlight, inflammation and epithelial adhesion

Gene variants that influence human biology today reflect thousands of years of evolution. Genetic effects on infant health are a major point of selective pressure, given that childhood survival is essential to evolutionary success. Knowledge of this evolutionary history can have implications for paediatric research. CONCLUSION: An episode of human adaptation to the extremely low ultraviolet radiation environment of the Arctic 20,000 years ago implicates the Ectodysplasin A Receptor (EDAR) and the Fatty Acid Desaturases (FADS) in human lactation and epithelial inflammation.

Genetic differentiation in East African ethnicities and its relationship with endurance running success

Since the 1960s, East African athletes, mainly from Kenya and Ethiopia, have dominated long-distance running events in both the male and female categories. Further demographic studies have shown that two ethnic groups are overrepresented among elite endurance runners in each of these countries: the Kalenjin, from Kenya, and the Oromo, from Ethiopia, raising the possibility that this dominance results from genetic or/and cultural factors. However, looking at the life history of these athletes or at loci previously associated with endurance athletic performance, no compelling explanation has emerged. Here, we used a population approach to identify peaks of genetic differentiation for these two ethnicities and compared the list of genes close to these regions with a list, manually curated by us, of genes that have been associated with traits possibly relevant to endurance running in GWAS studies, and found a significant enrichment in both populations (Kalenjin, P = 0.048, and Oromo, P = 1.6x10^-5). Those traits are mainly related to anthropometry, circulatory and respiratory systems, energy metabolism, and calcium homeostasis. Our results reinforce the notion that endurance running is a systemic activity with a complex genetic architecture, and indicate new candidate genes for future studies. Finally, we argue that a deterministic relationship between genetics and sports must be avoided, as it is both scientifically incorrect and prone to reinforcing population (racial) stereotyping.

Disentangling Signatures of Selection Before and After European Colonization in Latin Americans

Throughout human evolutionary history, large-scale migrations have led to intermixing (i.e., admixture) between previously separated human groups. Although classical and recent work have shown that studying admixture can yield novel historical insights, the extent to which this process contributed to adaptation remains underexplored. Here, we introduce a novel statistical model, specific to admixed populations, that identifies loci under selection while determining whether the selection likely occurred post-admixture or prior to admixture in one of the ancestral source populations. Through extensive simulations, we show that this method is able to detect selection, even in recently formed admixed populations, and to accurately differentiate between selection occurring in the ancestral or admixed population. We apply this method to genome-wide SNP data of ∼4,000 individuals in five admixed Latin American cohorts from Brazil, Chile, Colombia, Mexico, and Peru. Our approach replicates previous reports of selection in the human leukocyte antigen region that are consistent with selection post-admixture. We also report novel signals of selection in genomic regions spanning 47 genes, reinforcing many of these signals with an alternative, commonly used local-ancestry-inference approach. These signals include several genes involved in immunity, which may reflect responses to endemic pathogens of the Americas and to the challenge of infectious disease brought by European contact. In addition, some of the strongest signals inferred to be under selection in the Native American ancestral groups of modern Latin Americans overlap with genes implicated in energy metabolism phenotypes, plausibly reflecting adaptations to novel dietary sources available in the Americas.

Novel alleles gained during the Beringian isolation period

During the Last Glacial Maximum, a small band of Siberians entered the Beringian corridor, where they persisted, isolated from gene flow, for several thousand years before expansion into the Americas. The ecological features of the Beringian environment, coupled with an extended period of isolation at small population size, would have provided evolutionary opportunity for novel genetic variation to arise as both rare standing variants and new mutations were driven to high frequency through both neutral and directed processes. Here we perform a full genome investigation of Native American populations in the Thousand Genomes Project Phase 3 to identify unique high frequency alleles that can be dated to an origin in Beringia. Our analyses demonstrate that descendant populations of Native Americans harbor 20,424 such variants, which is on a scale comparable only to Africa and the Out of Africa bottleneck. This is consistent with simulations of a serial founder effects model. Tests for selection reveal that some of these Beringian variants were likely driven to high frequency by adaptive processes, and bioinformatic analyses suggest possible phenotypic pathways that were under selection during the Beringian Isolation period. Specifically, pathways related to cardiac processes and melanocyte function appear to be enriched for selected Beringian variants.

Haplotypes of single cancer driver genes and their local ancestry in a highly admixed long-lived population of Northeast Brazil

Admixed populations have not been examined in detail in cancer genetic studies. Here, we inferred the local ancestry of cancer-associated single nucleotide polymorphisms (SNPs) and haplotypes of a highly admixed Brazilian population. SNP array was used to genotype 73 unrelated individuals aged 80-102 years. Local ancestry inference was performed by merging genotyped regions with phase three data from the 1000 Genomes Project Consortium using RFmix. The average ancestry tract length was 9.12-81.71 megabases. Strong linkage disequilibrium was detected in 48 haplotypes containing 35 SNPs in 10 cancer driver genes. All together, 19 risk and eight protective alleles were identified in 23 out of 48 haplotypes. Homozygous individuals were mainly of European ancestry, whereas heterozygotes had at least one Native American and one African ancestry tract. Native-American ancestry for homozygous individuals with risk alleles for HNF1B, CDH1, and BRCA1 was inferred for the first time. Results indicated that analysis of SNP polymorphism in the present admixed population has a high potential to identify new ancestry-associated alleles and haplotypes that modify cancer susceptibility differentially in distinct human populations. Future case-control studies with populations with a complex history of admixture could help elucidate ancestry-associated biological differences in cancer incidence and therapeutic outcomes.

Infant Metabolome in Relation to Prenatal DHA Supplementation and Maternal Single-Nucleotide Polymorphism rs174602: Secondary Analysis of a Randomized Controlled Trial in Mexico

BACKGROUND

Although DHA (22:6n-3) is critical for fetal development, results from randomized controlled trials (RCTs) of prenatal DHA supplementation report inconsistent effects on offspring health. Variants in fatty acid desaturase (FADS) genes that regulate the conversion of n-3 and n-6 essential fatty acids into their biologically active derivatives may explain this heterogeneity.

OBJECTIVES

We investigated the effect of prenatal DHA supplementation on the offspring metabolome at age 3 mo and explored differences by maternal FADS single-nucleotide polymorphism (SNP) rs174602.

METHODS

Data were obtained from a double-blind RCT in Mexico [POSGRAD (Prenatal Omega-3 Fatty Acid Supplementation and Child Growth and Development)] in which women (18-35 y old) received DHA (400 mg/d) or placebo from mid-gestation until delivery. Using high-resolution MS with LC, untargeted metabolomics was performed on 112 offspring plasma samples. Discriminatory metabolic features were selected via linear regression (P < 0.05) with false discovery rate (FDR) correction (q = 0.2). Interaction by SNP rs174602 was assessed using 2-factor ANOVA. Stratified analyses were performed, where the study population was grouped into carriers (TT, TC; n = 70) and noncarriers (CC; n = 42) of the minor allele. Pathway enrichment analysis was performed with Mummichog (P < 0.05).

RESULTS

After FDR correction, there were no differences in metabolic features between infants whose mothers received prenatal DHA (n = 58) and those whose mothers received placebo (n = 54). However, we identified 343 differentially expressed features in the interaction analysis after FDR correction. DHA supplementation positively enriched amino acid and aminosugars metabolism pathways and decreased fatty acid metabolism pathways among offspring of minor allele carriers and decreased metabolites within the tricarboxylic acid cycle and galactose metabolism pathways among offspring of noncarriers.

CONCLUSIONS

Our findings demonstrate differences in infant metabolism in response to prenatal DHA supplementation by maternal SNP rs174602 and further support the need to incorporate genetic analysis of FADS polymorphisms into DHA supplementation trials.This trial was registered at clinicaltrials.gov as NCT00646360.

Sex differences in adipose development in a hunter-gatherer population

OBJECTIVE

Humans are unusually sexually dimorphic in body composition, with adult females having on average nearly twice the fat mass as males. The development of adipose sex differences has been well characterized for children growing up in food-abundant environments, with less known about cross-cultural variation, particularly in populations without exposure to market foods, mechanized technologies, schooling, vaccination, or other medical interventions.

METHODS

To add to the existing cross-cultural data, we fit multiple growth curves to body composition and anthropometric data to describe adipose development for the Savanna Pumé, South American hunter-gatherers.

RESULTS

(1) Little evidence is found for an adiposity 'rebound' at the end of early childhood among either Savanna Pumé girls or boys. (2) Rather, fat deposition fluctuates during childhood, from age ~4 to ~9 years, with no appreciable accumulation until the onset of puberty, a pattern also observed among Congo Baka hunter-gatherers. (3) Body fat fluctuations are more pronounced for girls than boys. (4) The age of peak skeletal, weight, and adipose gains are staggered to a much greater extent among the Savanna Pumé compared to the National Health and Nutrition Examination Survey (NHANES III) reference, suggesting this is an important developmental strategy in lean populations.

CONCLUSION

Documenting growth patterns under diverse preindustrial energetic conditions provides an important baseline for understanding sex differences in body fat emerging today under food abundance.

How natural selection shapes genetic differentiation in the MHC region: A case study with Native Americans

The Human Leukocyte Antigen (HLA) loci are extremely well documented targets of balancing selection, yet few studies have explored how selection affects population differentiation at these loci. In the present study we investigate genetic differentiation at HLA genes by comparing differentiation at microsatellites distributed genomewide to those in the MHC region. Our study uses a sample of 494 individuals from 30 human populations, 28 of which are Native Americans, all of whom were typed for genomewide and MHC region microsatellites. We find greater differentiation in the MHC than in the remainder of the genome (FST-MHC = 0.130 and FST-Genomic = 0.087), and use a permutation approach to show that this difference is statistically significant, and not accounted for by confounding factors. This finding lies in the opposite direction to the expectation that balancing selection reduces population differentiation. We interpret our findings as evidence that selection favors different sets of alleles in distinct localities, leading to increased differentiation. Thus, balancing selection at HLA genes simultaneously increases intra-population polymorphism and inter-population differentiation in Native Americans.

The evolutionary ecology of fatty-acid variation: Implications for consumer adaptation and diversification

The nutritional diversity of resources can affect the adaptive evolution of consumer metabolism and consumer diversification. The omega-3 long-chain polyunsaturated fatty acids eicosapentaenoic acid (EPA; 20:5n-3) and docosahexaenoic acid (DHA; 22:6n-3) have a high potential to affect consumer fitness, through their widespread effects on reproduction, growth and survival. However, few studies consider the evolution of fatty acid metabolism within an ecological context. In this review, we first document the extensive diversity in both primary producer and consumer fatty acid distributions amongst major ecosystems, between habitats and amongst species within habitats. We highlight some of the key nutritional contrasts that can shape behavioural and/or metabolic adaptation in consumers, discussing how consumers can evolve in response to the spatial, seasonal and community-level variation of resource quality. We propose a hierarchical trait-based approach for studying the evolution of consumers' metabolic networks and review the evolutionary genetic mechanisms underpinning consumer adaptation to EPA and DHA distributions. In doing so, we consider how the metabolic traits of consumers are hierarchically structured, from cell membrane function to maternal investment, and have strongly environment-dependent expression. Finally, we conclude with an outlook on how studying the metabolic adaptation of consumers within the context of nutritional landscapes can open up new opportunities for understanding evolutionary diversification.

How ancestral subsistence strategies solve salmon starvation and the "protein problem" of Pacific Rim resources

This article provides a theoretical treatment of hunter-gatherer diet and physiology. Through a synthesis of nutritional studies, informed by ethno-archaeological data, we examine the risk of protein-rich diets for human survival, and how societies circumvent "salmon starvation" in the northeastern Pacific Rim. Fundamental nutritional constraints associated with salmon storage and consumption counter long-standing assumptions about the engine of cultural evolution in the region. Excess consumption of lean meat can lead to protein poisoning, termed by early explorers "rabbit starvation." While consumption of fats and carbohydrates is widely portrayed as a pathway to "offsetting" protein thresholds, there are true limits to the amount of protein individuals can consume, and constraints are most extreme for smaller individuals, children, and pregnant/nursing mothers. While this problem is not usually perceived as associated with fish, the risk of protein poisoning limits the amount of low-fat fish that people can eat safely. Compared with smaller, mass-harvested species (e.g., eulachon), dried salmon are exceedingly lean. Under certain circumstances fattier foods (small forage fish, marine mammals, whales, and even bears) or carbohydrate-rich plants may have been preferred not just for taste but to circumvent this "dietary protein ceiling." Simply put, "salmon specialization" cannot evolve without access to complimentary caloric energy through fat-rich or carbohydrate-rich resources. By extension, the evolution of storage-based societies requires this problem be solved prior to or in tandem with-salmon intensification. Without such solutions, increased mortality and reproductive rates would have made salmon reliance unsustainable. This insight is in line with genomic research suggesting protein toxicity avoidance was a powerful evolutionary force, possibly linked to genetic adaptations among First Americans. It is also relevant to evaluating the plausibility of other purportedly "focal" economies and informs understanding of the many solutions varied global societies have engineered to overcome physiological protein limits.

Limited Evidence for Selection at the FADS Locus in Native American Populations

The FADS locus contains the genes FADS1 and FADS2 that encode enzymes involved in the synthesis of long-chain polyunsaturated fatty acids. This locus appears to have been a repeated target of selection in human evolution, likely because dietary input of long-chain polyunsaturated fatty acids varied over time depending on environment and subsistence strategy. Several recent studies have identified selection at the FADS locus in Native American populations, interpreted as evidence for adaptation during or subsequent to the passage through Beringia. Here, we show that these signals are confounded by independent selection—postdating the split from Native Americans—in the European and, possibly, the East Asian populations used in the population branch statistic test. This is supported by direct evidence from ancient DNA that one of the putatively selected haplotypes was already common in Northern Eurasia at the time of the separation of Native American ancestors. An explanation for the present-day distribution of the haplotype that is more consistent with the data is that Native Americans retain the ancestral state of Paleolithic Eurasians. Another haplotype at the locus may reflect a secondary selection signal, although its functional impact is unknown.

Adaptation and co-adaptation of skin pigmentation and vitamin D genes in native Americans

We carried out an exhaustive review regarding human skin color variation and how much it may be related to vitamin D metabolism and other photosensitive molecules. We discuss evolutionary contexts that modulate this variability and hypotheses postulated to explain them; for example, a small amount of melanin in the skin facilitates vitamin D production, making it advantageous to have fair skin in an environment with little radiation incidence. In contrast, more melanin protects folate from degradation in an environment with a high incidence of radiation. Some Native American populations have a skin color at odds with what would be expected for the amount of radiation in the environment in which they live, a finding challenging the so-called "vitamin D-folate hypothesis." Since food is also a source of vitamin D, dietary habits should also be considered. Here we argue that a gene network approach provides tools to explain this phenomenon since it indicates potential alleles co-evolving in a compensatory way. We identified alleles of the vitamin D metabolism and pigmentation pathways segregated together, but in different proportions, in agriculturalists and hunter-gatherers. Finally, we highlight how an evolutionary approach can be useful to understand current topics of medical interest.

The role of the microbiota in human genetic adaptation

As human populations spread across the world, they adapted genetically to local conditions. So too did the resident microorganism communities that everyone carries with them. However, the collective influence of the diverse and dynamic community of resident microbes on host evolution is poorly understood. The taxonomic composition of the microbiota varies among individuals and displays a range of sometimes redundant functions that modify the physicochemical environment of the host and may alter selection pressures. Here we review known human traits and genes for which the microbiota may have contributed or responded to changes in host diet, climate, or pathogen exposure. Integrating host–microbiota interactions in human adaptation could offer new approaches to improve our understanding of human health and evolution.

Polyunsaturated fatty acid biosynthesis pathway and genetics. implications for interindividual variability in prothrombotic, inflammatory conditions such as COVID-19✰,✰✰,★,★★

COVID-19 symptoms vary from silence to rapid death, the latter mediated by both a cytokine storm and a thrombotic storm. SARS-CoV (2003) induces Cox-2, catalyzing the synthesis, from highly unsaturated fatty acids (HUFA), of eicosanoids and docosanoids that mediate both inflammation and thrombosis. HUFA balance between arachidonic acid (AA) and other HUFA is a likely determinant of net signaling to induce a healthy or runaway physiological response. AA levels are determined by a non-protein coding regulatory polymorphisms that mostly affect the expression of FADS1, located in the FADS gene cluster on chromosome 11. Major and minor haplotypes in Europeans, and a specific functional insertion-deletion (Indel), rs66698963, consistently show major differences in circulating AA (>50%) and in the balance between AA and other HUFA (47-84%) in free living humans; the indel is evolutionarily selective, probably based on diet. The pattern of fatty acid responses is fully consistent with specific genetic modulation of desaturation at the FADS1-mediated 20:3→20:4 step. Well established principles of net tissue HUFA levels indicate that the high linoleic acid and low alpha-linoleic acid in populations drive the net balance of HUFA for any individual. We predict that fast desaturators (insertion allele at rs66698963; major haplotype in Europeans) are predisposed to higher risk and pathological responses to SARS-CoV-2 could be reduced with high dose omega-3 HUFA.

Postadmixture Selection on Chileans Targets Haplotype Involved in Pigmentation, Thermogenesis and Immune Defense against Pathogens

Detection of positive selection signatures in populations around the world is helping to uncover recent human evolutionary history as well as the genetic basis of diseases. Most human evolutionary genomic studies have been performed in European, African, and Asian populations. However, populations with Native American ancestry have been largely underrepresented. Here, we used a genome-wide local ancestry enrichment approach complemented with neutral simulations to identify postadmixture adaptations underwent by admixed Chileans through gene flow from Europeans into local Native Americans. The top significant hits (P = 2.4×10-7) are variants in a region on chromosome 12 comprising multiple regulatory elements. This region includes rs12821256, which regulates the expression of KITLG, a well-known gene involved in lighter hair and skin pigmentation in Europeans as well as in thermogenesis. Another variant from that region is associated with the long noncoding RNA RP11-13A1.1, which has been specifically involved in the innate immune response against infectious pathogens. Our results suggest that these genes were relevant for adaptation in Chileans following the Columbian exchange.

Omega-3 PUFA Responders and Non-Responders and the Prevention of Lipid Dysmetabolism and Related Diseases

The long-chain omega-3 polyunsaturated fatty acids (LC-omega-3 PUFAs) eicosapentaenoic acid and docosahexaenoic acid are the most popular dietary supplements recommended for the prevention/management of lipid dysmetabolisms and related diseases. However, remarkable inconsistencies exist among the outcomes of the human intervention studies in this field, which contrast with the impressive homogeneity of positive results of most of the preclinical studies. In the present review, we will firstly examine a series of factors-such as background diet composition, gut microbiota and genetic/epigenetic variants, which may lie beneath these inconsistencies. Moreover, we will discuss the recent advance in the knowledge of possible specific biomarkers (genetic-, epigenetic- and microbiota-related) that are being investigated with the goal to apply them in a personalized supplementation with omega-3 PUFAs. We will also consider the possibility of using already available parameters (Omega-3 index, Omega-6 PUFA/Omega-3 PUFA ratio) able to predict the individual responsiveness to these fatty acids and will discuss the optimal timing for their use. Finally, we will critically examine the results of those human studies that have already adopted the distinction of the subjects into omega-3 PUFA responders and non-responders and will discuss the advantage of using such an approach.

Prospective clinical trial examining the impact of genetic variation in FADS1 on the metabolism of linoleic acid- and ɣ-linolenic acid-containing botanical oils

BACKGROUND

Unexplained heterogeneity in clinical trials has resulted in questions regarding the effectiveness of ɣ-linolenic acid (GLA)-containing botanical oil supplements. This heterogeneity may be explained by genetic variation within the fatty acid desaturase (FADS) gene cluster that is associated with circulating and tissue concentrations of arachidonic acid (ARA) and dihomo-ɣ-linolenic acid (DGLA), both of which may be synthesized from GLA and result in proinflammatory and anti-inflammatory metabolites, respectively.

OBJECTIVES

The objective of this study was to prospectively compare the capacity of a non-Hispanic white cohort, stratified by FADS genotype at the key single-nucleotide polymorphism (SNP) rs174537, to metabolize 18-carbon omega-6 (n-6) PUFAs in borage oil (BO) and soybean oil (SO) to GLA, DGLA, and ARA.

METHODS

Healthy adults (n = 64) participated in a randomized, double-blind, crossover intervention. Individuals received encapsulated BO (Borago officinalis L.; 37% LA and 23% GLA) or SO [Glycine max (L.) Merr.; 50% LA and 0% GLA] for 4 wk, followed by an 8-wk washout period, before consuming the opposite oil for 4 wk. Serum lipids and markers of inflammation (C-reactive protein) were assessed for both oil types at baseline and during weeks 2 and 4 of the intervention.

RESULTS

SO supplementation failed to alter circulating concentrations of any n-6 long-chain PUFAs. In contrast, a modest daily dose of BO elevated serum concentrations of GLA and DGLA in an rs174537 genotype-dependent manner. In particular, DGLA increased by 57% (95% CI: 0.38, 0.79) in GG genotype individuals, but by 141% (95% CI: 1.03, 2.85) in TT individuals. For ARA, baseline concentrations varied substantially by genotype and increased modestly with BO supplementation, suggesting a key role for FADS variation in the balance of DGLA and ARA.

CONCLUSIONS

The results of this study clearly suggest that personalized and population-based approaches considering FADS genetic variation may be necessary to optimize the design of future clinical studies with GLA-containing oils. This trial was registered at clinicaltrials.gov as NCT02337231.

Inuit metabolism revisited: what drove the selective sweep of CPT1a L479?

This article reassesses historical studies of Inuit metabolism in light of recent developments in evolutionary genetics. It discusses the possible selective advantage of a variant of CPT1a, which encodes the rate limiting enzyme in hepatic fatty acid oxidation. The L479 variant of CPT1a underwent one of the strongest known selective sweeps in human history and is specific to Inuit and Yu'pik populations. Recent hypotheses predict that this variant may have been selected in response to possible detrimental effects of chronic ketosis in communities with very low carbohydrate consumption. Assessing these hypotheses alongside several alternative explanations of the selective sweep, this article challenges the notion that the selection of L479 is linked to predicted detrimental effects of ketosis. Bringing together for the first time data from biochemical, metabolic, and physiological studies inside and outside the Inuit sphere, it aims to provide a broader interpretative framework and a more comprehensive way to understand the selective sweep. It suggests that L479 may have provided a selective advantage in glucose conservation as part of a metabolic adaptation to very low carbohydrate and high protein consumption, but not necessarily a ketogenic state, in an extremely cold environment. A high intake of n-3 fatty acids may be linked to selection through the mitigation of a detrimental effect of the mutation that arises in the fasted state. The implications of these conclusions for our broader understanding of very low carbohydrate metabolism, and for dietary recommendations for Inuit and non-Inuit populations, are discussed.

Adaptation to Extreme Environments in an Admixed Human Population from the Atacama Desert

Inorganic arsenic (As) is a toxic xenobiotic and carcinogen associated with severe health conditions. The urban population from the Atacama Desert in northern Chile was exposed to extremely high As levels (up to 600 µg/l) in drinking water between 1958 and 1971, leading to increased incidence of urinary bladder cancer (BC), skin cancer, kidney cancer, and coronary thrombosis decades later. Besides, the Andean Native-American ancestors of the Atacama population were previously exposed for millennia to elevated As levels in water (∼120 µg/l) for at least 5,000 years, suggesting adaptation to this selective pressure. Here, we performed two genome-wide selection tests—PBS_n1 and an ancestry-enrichment test—in an admixed population from Atacama, to identify adaptation signatures to As exposure acquired before and after admixture with Europeans, respectively. The top second variant selected by PBS_n1 was associated with LCE4A-C1orf68, a gene that may be involved in the immune barrier of the epithelium during BC. We performed association tests between the top PBS_n1 hits and BC occurrence in our population. The strongest association (P = 0.012) was achieved by the LCE4A-C1orf68 variant. The ancestry-enrichment test detected highly significant signals (P = 1.3 × 10⁻⁹) mapping MAK16, a gene with important roles in ribosome biogenesis during the G1 phase of the cell cycle. Our results contribute to a better understanding of the genetic factors involved in adaptation to the pathophysiological consequences of As exposure.

An ancestry informative marker panel design for individual ancestry estimation of Hispanic population using whole exome sequencing data

Background

Europeans and American Indians were major genetic ancestry of Hispanics in the U.S. These ancestral groups have markedly different incidence rates and outcomes in many types of cancers. Therefore, the genetic admixture may cause biased genetic association study with cancer susceptibility variants specifically in Hispanics. For example, the incidence rate of liver cancer has been shown with substantial disparity between Hispanic, Asian and non-Hispanic white populations. Currently, ancestry informative marker (AIM) panels have been widely utilized with up to a few hundred ancestry-informative single nucleotide polymorphisms (SNPs) to infer ancestry admixture. Notably, current available AIMs are predominantly located in intron and intergenic regions, while the whole exome sequencing (WES) protocols commonly used in translational research and clinical practice do not cover these markers. Thus, it remains challenging to accurately determine a patient’s admixture proportion without additional DNA testing.

Results

In this study we designed an unique AIM panel that infers 3-way genetic admixture from three distinct and selective continental populations (African (AFR), European (EUR), and East Asian (EAS)) within evolutionarily conserved exonic regions. Initially, about 1 million exonic SNPs from selective three populations in the 1000 Genomes Project were trimmed by their linkage disequilibrium (LD), restricted to biallelic variants, and finally we optimized to an AIM panel with 250 SNP markers, or the UT-AIM250 panel, using their ancestral informativeness statistics. Comparing to published AIM panels, UT-AIM250 performed better accuracy when we tested with three ancestral populations (accuracy: 0.995 ± 0.012 for AFR, 0.997 ± 0.007 for EUR, and 0.994 ± 0.012 for EAS). We further demonstrated the performance of the UT-AIM250 panel to admixed American (AMR) samples of the 1000 Genomes Project and obtained similar results (AFR, 0.085 ± 0.098; EUR, 0.665 ± 0.182; and EAS, 0.250 ± 0.205) to previously published AIM panels (Phillips-AIM34: AFR, 0.096 ± 0.127, EUR, 0.575 ± 0.290, and EAS, 0.330 ± 0.315; Wei-AIM278: AFR, 0.070 ± 0.096, EUR, 0.537 ± 0.267, and EAS, 0.393 ± 0.300). Subsequently, we applied the UT-AIM250 panel to a clinical dataset of 26 self-reported Hispanic patients in South Texas with hepatocellular carcinoma (HCC). We estimated the admixture proportions using WES data of adjacent non-cancer liver tissues (AFR, 0.065 ± 0.043; EUR, 0.594 ± 0.150; and EAS, 0.341 ± 0.160). Similar admixture proportions were identified from corresponding tumor tissues. In addition, we estimated admixture proportions of The Cancer Genome Atlas (TCGA) collection of hepatocellular carcinoma (TCGA-LIHC) samples (376 patients) using the UT-AIM250 panel. The panel obtained consistent admixture proportions from tumor and matched normal tissues, identified 3 possible incorrectly reported race/ethnicity, and/or provided race/ethnicity determination if necessary.

Conclusions

Here we demonstrated the feasibility of using evolutionarily conserved exonic regions to infer admixture proportions and provided a robust and reliable control for sample collection or patient stratification for genetic analysis. R implementation of UT-AIM250 is available at https://github.com/chenlabgccri/UT-AIM250.

Evolution of Hominin Polyunsaturated Fatty Acid Metabolism: From Africa to the New World

The metabolic conversion of dietary omega-3 and omega-6 18 carbon (18C) to long chain (>20 carbon) polyunsaturated fatty acids (LC-PUFAs) is vital for human life. The rate-limiting steps of this process are catalyzed by fatty acid desaturase (FADS) 1 and 2. Therefore, understanding the evolutionary history of the FADS genes is essential to our understanding of hominin evolution. The FADS genes have two haplogroups, ancestral and derived, with the derived haplogroup being associated with more efficient LC-PUFA biosynthesis than the ancestral haplogroup. In addition, there is a complex global distribution of these haplogroups that is suggestive of Neanderthal introgression. We confirm that Native American ancestry is nearly fixed for the ancestral haplogroup, and replicate a positive selection signal in Native Americans. This positive selection potentially continued after the founding of the Americas, although simulations suggest that the timing is dependent on the allele frequency of the ancestral Beringian population. We also find that the Neanderthal FADS haplotype is more closely related to the derived haplogroup and the Denisovan clusters closer to the ancestral haplogroup. Furthermore, the derived haplogroup has a time to the most recent common ancestor of 688,474 years before present. These results support an ancient polymorphism, as opposed to Neanderthal introgression, forming in the FADS region during the Pleistocene with possibly differential selection pressures on both haplogroups. The near fixation of the ancestral haplogroup in Native American ancestry calls for future studies to explore the potential health risk of associated low LC-PUFA levels in these populations.

FADS1 and FADS2 Polymorphisms Modulate Fatty Acid Metabolism and Dietary Impact on Health

Variants in the FADS gene cluster modify the activity of polyunsaturated fatty acid (PUFA) desaturation and the lipid composition in human blood and tissue. FADS variants have been associated with plasma lipid concentrations, risk of cardiovascular diseases, overweight, eczema, pregnancy outcomes, and cognitive function. Studies on variations in the FADS genecluster provided some of the first examples for marked gene-diet interactions in modulating complex phenotypes, such as eczema, asthma, and cognition. Genotype distribution differs markedly among ethnicities, apparently reflecting an evolutionary advantage of genotypes enabling active long-chain PUFA synthesis when the introduction of agriculture provided diets rich in linoleic acid but with little arachidonic and eicosapentaenoic acids. Discovering differential effects of PUFA supply that depend on variation of FADS genotypes could open new opportunities for developing precision nutrition strategies based either on an individual's genotype or on genotype distributions in specific populations.

Footprints of adaptive evolution revealed by whole Z chromosomes haplotypes in flycatchers

Detecting positive selection using genomic data is critical to understanding the role of adaptive evolution. Of particular interest in this context is sex chromosomes since they are thought to play a special role in local adaptation and speciation. We sought to circumvent the challenges associated with statistical phasing when using haplotype-based statistics in sweep scans by benefitting from that whole chromosome haplotypes of the sex chromosomes can be obtained by resequencing of individuals of the hemizygous sex. We analyzed whole Z chromosome haplotypes from 100 females from several populations of four black and white flycatcher species (in birds, females are ZW and males ZZ). Based on integrated haplotype score (iHS) and number of segregating sites by length (nSL) statistics, we found strong and frequent haplotype structure in several regions of the Z chromosome in each species. Most of these sweep signals were population-specific, with essentially no evidence for regions under selection shared among species. Some completed sweeps were revealed by the cross-population extended haplotype homozygosity (XP-EHH) statistic. Importantly, by using statistically phased Z chromosome data from resequencing of males, we failed to recover the signals of selection detected in analyses based on whole chromosome haplotypes from females; instead, what likely represent false signals of selection were frequently seen. This highlights the power issues in statistical phasing and cautions against conclusions from selection scans using such data. The detection of frequent selective sweeps on the avian Z chromosome supports a large role of sex chromosomes in adaptive evolution.

Nutrition and its role in human evolution

Our understanding of human evolution has improved rapidly over recent decades, facilitated by large-scale cataloguing of genomic variability amongst both modern and archaic humans. It seems clear that the evolution of the ancestors of chimpanzees and hominins separated 7-9 million years ago with some migration out of Africa by the earlier hominins; Homo sapiens slowly emerged as climate change resulted in drier, less forested African conditions. The African populations expanded and evolved in many different conditions with slow mutation and selection rates in the human genome, but with much more rapid mutation occurring in mitochondrial DNA. We now have evidence stretching back 300 000 years of humans in their current form, but there are clearly four very different large African language groups that correlate with population DNA differences. Then, about 50 000-100 000 years ago a small subset of modern humans also migrated out of Africa resulting in a persistent signature of more limited genetic diversity amongst non-African populations. Hybridization with archaic hominins occurred around this time such that all non-African modern humans possess some Neanderthal ancestry and Melanesian populations additionally possess some Denisovan ancestry. Human populations both within and outside Africa also adapted to diverse aspects of their local environment including altitude, climate, UV exposure, diet and pathogens, in some cases leaving clear signatures of patterns of genetic variation. Notable examples include haemoglobin changes conferring resistance to malaria, other immune changes and the skin adaptations favouring the synthesis of vitamin D. As humans migrated across Eurasia, further major mitochondrial changes occurred with some interbreeding with ancient hominins and the development of alcohol intolerance. More recently, an ability to retain lactase persistence into adulthood has evolved rapidly under the environmental stimulus of pastoralism with the ability to husband lactating ruminants. Increased amylase copy numbers seem to relate to the availability of starchy foods, whereas the capacity to desaturase and elongate monounsaturated fatty acids in different societies seems to be influenced by whether there is a lack of supply of readily available dietary sources of long-chain polyunsaturated fatty acids. The process of human evolution includes genetic drift and adaptation to local environments, in part through changes in mitochondrial and nuclear DNA. These genetic changes may underlie susceptibilities to some modern human pathologies including folate-responsive neural tube defects, diabetes, other age-related pathologies and mental health disorders.

The Genetic Basis of Metabolic Disease

Recent developments in genetics and genomics are providing a detailed and systematic characterization of the genetic underpinnings of common metabolic diseases and traits, highlighting the inherent complexity within systems for homeostatic control and the many ways in which that control can fail. The genetic architecture underlying these common metabolic phenotypes is complex, with each trait influenced by hundreds of loci spanning a range of allele frequencies and effect sizes. Here, we review the growing appreciation of this complexity and how this has fostered the implementation of genome-scale approaches that deliver robust mechanistic inference and unveil new strategies for translational exploitation.

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.

FADS1 and the Timing of Human Adaptation to Agriculture

Variation at the FADS1/FADS2 gene cluster is functionally associated with differences in lipid metabolism and is often hypothesized to reflect adaptation to an agricultural diet. Here, we test the evidence for this relationship using both modern and ancient DNA data. We show that almost all the inhabitants of Europe carried the ancestral allele until the derived allele was introduced ∼8,500 years ago by Early Neolithic farming populations. However, we also show that it was not under strong selection in these populations. We find that this allele, and other proposed agricultural adaptations at LCT/MCM6 and SLC22A4, were not strongly selected until much later, perhaps as late as the Bronze Age. Similarly, increased copy number variation at the salivary amylase gene AMY1 is not linked to the development of agriculture although, in this case, the putative adaptation precedes the agricultural transition. Our analysis shows that selection at the FADS locus was not tightly linked to the initial introduction of agriculture and the Neolithic transition. Further, it suggests that the strongest signals of recent human adaptation in Europe did not coincide with the Neolithic transition but with more recent changes in environment, diet, or efficiency of selection due to increases in effective population size.

Reconstructing the Deep Population History of Central and South America

We report genome-wide ancient DNA from 49 individuals forming four parallel time transects in Belize, Brazil, the Central Andes, and the Southern Cone, each dating to at least ∼9,000 years ago. The common ancestral population radiated rapidly from just one of the two early branches that contributed to Native Americans today. We document two previously unappreciated streams of gene flow between North and South America. One affected the Central Andes by ∼4,200 years ago, while the other explains an affinity between the oldest North American genome associated with the Clovis culture and the oldest Central and South Americans from Chile, Brazil, and Belize. However, this was not the primary source for later South Americans, as the other ancient individuals derive from lineages without specific affinity to the Clovis-associated genome, suggesting a population replacement that began at least 9,000 years ago and was followed by substantial population continuity in multiple regions.

Selection scan reveals three new loci related to high altitude adaptation in Native Andeans

The Andean Altiplano has been occupied continuously since the late Pleistocene, ~12,000 years ago, which places the Andean natives as one of the most ancient populations living at high altitudes. In the present study, we analyzed genomic data from Native Americans living a long-time at Andean high altitude and at Amazonia and Mesoamerica lowland areas. We have identified three new candidate genes - SP100, DUOX2 and CLC - with evidence of positive selection for altitude adaptation in Andeans. These genes are involved in the TP53 pathway and are related to physiological routes important for high-altitude hypoxia response, such as those linked to increased angiogenesis, skeletal muscle adaptations, and immune functions at the fetus-maternal interface. Our results, combined with other studies, showed that Andeans have adapted to the Altiplano in different ways and using distinct molecular strategies as compared to those of other natives living at high altitudes.

Evolutionary history and adaptation of a human pygmy population of Flores Island, Indonesia

Flores Island, Indonesia, was inhabited by the small-bodied hominin species Homo floresiensis, which has an unknown evolutionary relationship to modern humans. This island is also home to an extant human pygmy population. Here we describe genome-scale single-nucleotide polymorphism data and whole-genome sequences from a contemporary human pygmy population living on Flores near the cave where H. floresiensis was found. The genomes of Flores pygmies reveal a complex history of admixture with Denisovans and Neanderthals but no evidence for gene flow with other archaic hominins. Modern individuals bear the signatures of recent positive selection encompassing the FADS (fatty acid desaturase) gene cluster, likely related to diet, and polygenic selection acting on standing variation that contributed to their short-stature phenotype. Thus, multiple independent instances of hominin insular dwarfism occurred on Flores.

Long-chain polyunsaturated fatty acids, gestation duration, and birth size: a Mendelian randomization study using fatty acid desaturase variants

Background

In randomized trials, supplementation of n-3 (ω-3) long-chain polyunsaturated fatty acids (LC-PUFAs) during pregnancy has resulted in increased size at birth, which is attributable to longer gestation.

Objective

We examined this finding by using a Mendelian randomization approach utilizing fatty acid desaturase (FADS) gene variants affecting LC-PUFA metabolism.

Design

As part of a tri-ethnic mother-offspring cohort in Singapore, 35 genetic variants in FADS1, FADS2, and FADS3 were genotyped in 898 mothers and 1103 offspring. Maternal plasma n-3 and n-6 PUFA concentrations at 26-28 wk of gestation were measured. Gestation duration was derived from an ultrasound dating scan in early pregnancy and from birth date. Birth length and weight were measured. Eight FADS variants were selected through a tagging-SNP approach and examined in association with PUFA concentrations, gestation duration among spontaneous labors, and birth size with the use of ethnicity-adjusted linear regressions and survival models that accounted for the competing risks of induced labor and prelabor cesarean delivery.

Results

Maternal FADS1 variant rs174546, tagging for 8 other variants located on FADS1 and FADS2, was strongly related to plasma n-6 but not n-3 LC-PUFA concentrations. Offspring and maternal FADS3 variants were associated with gestation duration among women who had spontaneous labor: each copy of rs174450 minor allele C was associated with a shorter gestation by 2.2 d (95% CI: 0.9, 3.4 d) and 1.9 d (0.7, 3.0 d) for maternal and offspring variants, respectively. In survival models, rs174450 minor allele homozygotes had reduced time to delivery after spontaneous labor compared with major allele homozygotes [HR (95% CI): 1.51 (1.18, 1.95) and 1.51 (1.20, 1.89) for mothers and offspring, respectively].

Conclusions

With the use of a Mendelian randomization approach, we observed associations between FADS variants and gestation duration. This suggests a potential role of LC-PUFAs in gestation duration. This trial was registered at http://www.clinicaltrials.gov as NCT01174875.

Environmental selection during the last ice age on the mother-to-infant transmission of vitamin D and fatty acids through breast milk

The frequency of the human-specific EDAR V370A isoform is highly elevated in North and East Asian populations. The gene is known to have several pleiotropic effects, among which are sweat gland density and ductal branching in the mammary gland. The former has led some geneticists to argue that the near-fixation of this allele was caused by selection for modulation of thermoregulatory sweating. We provide an alternative hypothesis, that selection instead acted on the allele’s effect of increasing ductal branching in the mammary gland, thereby amplifying the transfer of critical nutrients to infants via mother’s milk. This is likely to have occurred during the Last Glacial Maximum when a human population was genetically isolated in the high-latitude environment of the Beringia.

Because of the ubiquitous adaptability of our material culture, some human populations have occupied extreme environments that intensified selection on existing genomic variation. By 32,000 years ago, people were living in Arctic Beringia, and during the Last Glacial Maximum (LGM; 28,000–18,000 y ago), they likely persisted in the Beringian refugium. Such high latitudes provide only very low levels of UV radiation, and can thereby lead to dangerously low levels of biosynthesized vitamin D. The physiological effects of vitamin D deficiency range from reduced dietary absorption of calcium to a compromised immune system and modified adipose tissue function. The ectodysplasin A receptor (EDAR) gene has a range of pleiotropic effects, including sweat gland density, incisor shoveling, and mammary gland ductal branching. The frequency of the human-specific EDAR V370A allele appears to be uniquely elevated in North and East Asian and New World populations due to a bout of positive selection likely to have occurred circa 20,000 y ago. The dental pleiotropic effects of this allele suggest an even higher occurrence among indigenous people in the Western Hemisphere before European colonization. We hypothesize that selection on EDAR V370A occurred in the Beringian refugium because it increases mammary ductal branching, and thereby may amplify the transfer of critical nutrients in vitamin D-deficient conditions to infants via mothers’ milk. This hypothesized selective context for EDAR V370A was likely intertwined with selection on the fatty acid desaturase (FADS) gene cluster because it is known to modulate lipid profiles transmitted to milk from a vitamin D-rich diet high in omega-3 fatty acids.

Exome Sequencing Provides Evidence of Polygenic Adaptation to a Fat-Rich Animal Diet in Indigenous Siberian Populations

Siberia is one of the coldest environments on Earth and has great seasonal temperature variation. Long-term settlement in northern Siberia undoubtedly required biological adaptation to severe cold stress, dramatic variation in photoperiod, and limited food resources. In addition, recent archeological studies show that humans first occupied Siberia at least 45,000 years ago; yet our understanding of the demographic history of modern indigenous Siberians remains incomplete. In this study, we use whole-exome sequencing data from the Nganasans and Yakuts to infer the evolutionary history of these two indigenous Siberian populations. Recognizing the complexity of the adaptive process, we designed a model-based test to systematically search for signatures of polygenic selection. Our approach accounts for stochasticity in the demographic process and the hitchhiking effect of classic selective sweeps, as well as potential biases resulting from recombination rate and mutation rate heterogeneity. Our demographic inference shows that the Nganasans and Yakuts diverged ∼12,000-13,000 years ago from East-Asian ancestors in a process involving continuous gene flow. Our polygenic selection scan identifies seven candidate gene sets with Siberian-specific signals. Three of these gene sets are related to diet, especially to fat metabolism, consistent with the hypothesis of adaptation to a fat-rich animal diet. Additional testing rejects the effect of hitchhiking and favors a model in which selection yields small allele frequency changes at multiple unlinked genes.