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

Genome diversity and signatures of natural selection in mainland Southeast Asia

Mainland Southeast Asia (MSEA) has rich ethnic and cultural diversity with a population of nearly 300 million1,2. However, people from MSEA are underrepresented in the current human genomic databases. Here we present the SEA3K genome dataset (phase I), generated by deep short-read whole-genome sequencing of 3,023 individuals from 30 MSEA populations, and long-read whole-genome sequencing of 37 representative individuals. We identified 79.59 million small variants and 96,384 structural variants, among which 22.83 million small variants and 24,622 structural variants are unique to this dataset. We observed a high genetic heterogeneity across MSEA populations, reflected by the varied combinations of genetic components. We identified 44 genomic regions with strong signatures of Darwinian positive selection, covering 89 genes involved in varied physiological systems such as physical traits and immune response. Furthermore, we observed varied patterns of archaic Denisovan introgression in MSEA populations, supporting the proposal of at least two distinct instances of Denisovan admixture into modern humans in Asia3. We also detected genomic regions that suggest adaptive archaic introgressions in MSEA populations. The large number of novel genomic variants in MSEA populations highlight the necessity of studying regional populations that can help answer key questions related to prehistory, genetic adaptation and complex diseases.

Population genomics of sika deer reveals recent speciation and genetic selective signatures during evolution and domestication

BACKGROUND

Population genomic analysis can reconstruct the phylogenetic relationship and demographic history, and identify genomic selective signatures of a species. To date, fundamental aspects of population genomic analyses, such as intraspecies taxonomy, evolutionary history, and adaptive evolution, of sika deer have not been systematically investigated. Furthermore, accumulating lines of evidences have illustrated that incorrect species delimitation will mislead conservation decisions, and even lead to irreversible mistakes in threatened species.

RESULTS

In this study, we resequenced 81 wild and 71 domesticated sika deer representing 10 main geographic populations and two farms to clarify the species delimitation, demographic and divergence histories, and adaptive evolution of this species. First, our analyses of whole genomes, Y chromosomes and mitochondrial genomes revealed substantial genetic differentiation between the continental and Japanese lineages of sika deer, representing two phylogenetically distinct species. Second, sika deer in Japan were inferred to have experienced a "divergence-mixing-isolation" evolutionary scenario. Third, we identified four candidate genes (XKR4, NPAS3, CTNNA3, and CNTNAP5) possibly involved in body size regulation of sika deer by selective sweep analysis. Furthermore, we also detected two candidate genes (NRP2 and EDIL3) that may be associated with an important economic trait (antler weight) were under selection during the process of domestication.

CONCLUSION

Population genomic analyses revealed that the continental and Japanese lineages represent distinct phylogenetic species. Moreover, our results provide insights into the genetic selection signatures related to body size differences and a valuable genomic resource for future genetic studies and genomics-informed breeding of sika deer.

Selection scan in Native Americans of Mexico identifies FADS2 rs174616: Evidence of gene-diet interactions affecting lipid levels and Delta-6-desaturase activity

Searching for positive selection signals across genomes has identified functional genetic variants responding to environmental change. In Native Americans of Mexico, we used the fixation index (Fst) and population branch statistic (PBS) to identify SNPs suggesting positive selection. The 103 most differentiated SNPs were tested for associations with metabolic traits, the most significant association was FADS2/rs174616 with body mass index (BMI). This variant lies within a linkage disequilibrium (LD) block independent of previously reported FADS selection signals and has not been clearly associated with metabolic phenotypes. We tested this variant in two independent cohorts with cardiometabolic data. In the Genetics of Atherosclerotic Disease (GEA) cohort, the derived allele (T) was associated with increased BMI, lower LDL-C levels and a decreased risk of subclinical atherosclerosis in women. Significant gene-diet interactions affected lipid, apolipoprotein and adiponectin levels with differences according to sex, involving mainly total and complex dietary carbohydrate%. In the Genotype-related Effects of PUFA trial, the derived allele was associated with lower Δ-6 desaturase activity and erythrocyte membrane dihomo-gamma-linolenic acid (DGLA) levels, and with increased Δ-5 desaturase activity and eicosapentaenoic acid levels. This variant interacted with dietary carbohydrate% affecting Δ-6 desaturase activity. Notably, the relationship of DGLA and other erythrocyte membrane LC-PUFA indices with HOMA-IR differed according to rs174616 genotype, which has implications regarding how these indices should be interpreted. In conclusion, this observational study identified rs174616 as a signal suggesting selection in an independent linkage disequilibrium block, was associated with cardiometabolic and erythrocyte measurements of LC-PUFA in two independent Mexican cohorts and showed significant gene-diet interactions.

Testing for differences in polygenic scores in the presence of confounding

Polygenic scores have become an important tool in human genetics, enabling the prediction of individuals' phenotypes from their genotypes. Understanding how the pattern of differences in polygenic score predictions across individuals intersects with variation in ancestry can provide insights into the evolutionary forces acting on the trait in question, and is important for understanding health disparities. However, because most polygenic scores are computed using effect estimates from population samples, they are susceptible to confounding by both genetic and environmental effects that are correlated with ancestry. The extent to which this confounding drives patterns in the distribution of polygenic scores depends on patterns of population structure in both the original estimation panel and in the prediction/test panel. Here, we use theory from population and statistical genetics, together with simulations, to study the procedure of testing for an association between polygenic scores and axes of ancestry variation in the presence of confounding. We use a general model of genetic relatedness to describe how confounding in the estimation panel biases the distribution of polygenic scores in a way that depends on the degree of overlap in population structure between panels. We then show how this confounding can bias tests for associations between polygenic scores and important axes of ancestry variation in the test panel. Specifically, for any given test, there exists a single axis of population structure in the GWAS panel that needs to be controlled for in order to protect the test. Based on this result, we propose a new approach for directly estimating this axis of population structure in the GWAS panel. We then use simulations to compare the performance of this approach to the standard approach in which the principal components of the GWAS panel genotypes are used to control for stratification.

Differentiated adaptative genetic architecture and language-related demographical history in South China inferred from 619 genomes from 56 populations

BACKGROUND

The underrepresentation of human genomic resources from Southern Chinese populations limited their health equality in the precision medicine era and complete understanding of their genetic formation, admixture, and adaptive features. Besides, linguistical and genetic evidence supported the controversial hypothesis of their origin processes. One hotspot case was from the Chinese Guangxi Pinghua Han people (GPH), whose language was significantly similar to Southern Chinese dialects but whose uniparental gene pool was phylogenetically associated with the indigenous Tai-Kadai (TK) people. Here, we analyzed genome-wide SNP data in 619 people from four language families and 56 geographically different populations, in which 261 people from 21 geographically distinct populations were first reported here.

RESULTS

We identified significant population stratification among ethnolinguistically diverse Guangxi populations, suggesting their differentiated genetic origin and admixture processes. GPH shared more alleles related to Zhuang than Southern Han Chinese but received more northern ancestry relative to Zhuang. Admixture models and estimates of genetic distances showed that GPH had a close genetic relationship with geographically close TK compared to Northern Han Chinese, supporting their admixture origin hypothesis. Further admixture time and demographic history reconstruction supported GPH was formed via admixture between Northern Han Chinese and Southern TK people. We identified robust signatures associated with lipid metabolisms, such as fatty acid desaturases (FADS) and medically relevant loci associated with Mendelian disorder (GJB2) and complex diseases. We also explored the shared and unique selection signatures of ethnically different but linguistically related Guangxi lineages and found some shared signals related to immune and malaria resistance.

CONCLUSIONS

Our genetic analysis illuminated the language-related fine-scale genetic structure and provided robust genetic evidence to support the admixture hypothesis that can explain the pattern of observed genetic diversity and formation of GPH. This work presented one comprehensive analysis focused on the population history and demographical adaptative process, which provided genetic evidence for personal health management and disease risk prediction models from Guangxi people. Further large-scale whole-genome sequencing projects would provide the entire landscape of southern Chinese genomic diversity and their contributions to human health and disease traits.

More than a decade of genetic research on the Denisovans

Denisovans, a group of now extinct humans who lived in Eastern Eurasia in the Middle and Late Pleistocene, were first identified from DNA sequences just over a decade ago. Only ten fragmentary remains from two sites have been attributed to Denisovans based entirely on molecular information. Nevertheless, there has been great interest in using genetic data to understand Denisovans and their place in human history. From the reconstruction of a single high-quality genome, it has been possible to infer their population history, including events of admixture with other human groups. Additionally, the identification of Denisovan DNA in the genomes of present-day individuals has provided insights into the timing and routes of dispersal of ancient modern humans into Asia and Oceania, as well as the contributions of archaic DNA to the physiology of present-day people. In this Review, we synthesize more than a decade of research on Denisovans, reconcile controversies and summarize insights into their population history and phenotype. We also highlight how our growing knowledge about Denisovans has provided insights into our own evolutionary history.

The relationship of fatty acids to ischaemic heart disease and lifespan in men and women using Mendelian randomization

BACKGROUND

Observationally, polyunsaturated fatty acids (PUFAs) have health benefits compared with saturated fatty acids (SFAs); randomized controlled trials suggest fewer benefits. We used uni- and multi-variable Mendelian randomization to assess the association of major fatty acids and their sub-species with ischaemic heart disease (IHD) overall and sex-specifically and with lifespan sex-specifically, given differing lifespan by sex.

METHODS

We obtained strong (P <5x10-8), independent (r2<0.001) genetic predictors of fatty acids from genome-wide association studies (GWAS) in a random subset of 114 999 UK Biobank participants. We applied these genetic predictors to the Cardiogram IHD GWAS (cases = 60 801, controls = 123 504) and to the Finngen consortium GWAS (cases = 31 640, controls = 187 152) for replication and to the UK Biobank for sex-specific IHD and for lifespan based on parental attained age (fathers = 415 311, mothers = 412 937). We used sensitivity analysis and assessed sex differences where applicable.

RESULTS

PUFAs were associated with IHD [odds ratio 1.23, 95% confidence interval (CI) 1.05 to 1.44] and lifespan in men (-0.76 years, 95% CI -1.34 to -0.17) but not women (0.20, 95% CI -0.32 to 0.70). Findings were similar for omega-6 fatty acids and linoleic acid. Independent associations of SFAs, mono-unsaturated fatty acids or omega-3 fatty acids with IHD overall or lifespan in men and women were limited.

CONCLUSIONS

PUFAs, via specific subspecies, may contribute to disparities in lifespan by sex. Sex-specific dietary advice might be a start towards personalized public health and addressing inequities.

Impact of FADS gene variation and dietary fatty acid exposure on biochemical and anthropomorphic phenotypes in a Hispanic/Latino cohort

Introduction

Polyunsaturated fatty acids (PUFA) and highly unsaturated fatty acid (HUFA) synthetic products and their signaling metabolites play vital roles in immunity, inflammation, and brain development/function. Frequency differences of variants within the fatty acid desaturase (FADS) gene cluster affect levels of HUFAs, their biologically active products, and numerous physiological phenotypes. Fundamental questions remain regarding the impact of this genetic variation on the health of Hispanic/Latino populations.

Methods

Data and biospecimens (plasma, red blood cells, buffy coat-derived DNA) from 135 participants (83.7% female) were used to assess the relationship(s) between dietary PUFA levels, a FADS haplotype tagging SNP, rs174537, and the capacity of Hispanic/Latino populations to generate HUFAs in plasma and RBC as well as its potential impact on anthropomorphic phenotypes.

Results

The dietary habits of the cohort showed that participant diets contained a high ratio (9.3 ± 0.2, mean ± SEM) of linoleic acid (n-6) to alpha-linolenic acid (n-3) and also contained extremely low levels of n-3 HUFAs (eicosapentaenoic acid, EPA and docosahexaenoic acid, DHA), both features of the Modern Western Diet. Compared to African and European American cohorts, the frequency of the TT rs174537 genotype was highly enriched (53% of subjects) in this Hispanic/Latino cohort and was strongly associated with lower circulating HUFA levels. For example, plasma levels of arachidonic acid (ARA: 20:4, n-6) and EPA (20:5, n-3) were 37% and 23%, respectively, lower in the TT versus the GG genotype. HUFA biosynthetic efficiency, as determined by metabolic product to precursor ratios, was highly dependent (p < 0.0001) on the rs174537 genotype (GG > GT > TT) for both circulating n-6 and n-3 HUFAs. In contrast, the RBC Omega-3 Index (EPA + DHA) was extremely low (2.89 ± 1.65, mean ± sd) in this population and independent of rs174537 genotype. Importantly, the rs174537 genotype was also related to female height with TT genotype participants being 4.5 cm shorter (p = 0.0001) than the GG + GT participants.

Discussion

Taken together, this study illustrates that dietary PUFA + HUFA × FADS gene- interactions place a large proportion (>50%) of Hispanic/Latino populations at high risk of a deficiency in both circulating and cellular levels of n-3 HUFAs.

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.

Spatial and temporal heterogeneity in human mobility patterns in Holocene Southwest Asia and the East Mediterranean

We present a spatiotemporal picture of human genetic diversity in Anatolia, Iran, Levant, South Caucasus, and the Aegean, a broad region that experienced the earliest Neolithic transition and the emergence of complex hierarchical societies. Combining 35 new ancient shotgun genomes with 382 ancient and 23 present-day published genomes, we found that genetic diversity within each region steadily increased through the Holocene. We further observed that the inferred sources of gene flow shifted in time. In the first half of the Holocene, Southwest Asian and the East Mediterranean populations homogenized among themselves. Starting with the Bronze Age, however, regional populations diverged from each other, most likely driven by gene flow from external sources, which we term "the expanding mobility model." Interestingly, this increase in inter-regional divergence can be captured by outgroup-f₃-based genetic distances, but not by the commonly used F_ST statistic, due to the sensitivity of F_ST, but not outgroup-f₃, to within-population diversity. Finally, we report a temporal trend of increasing male bias in admixture events through the Holocene.

Tracing Bai-Yue Ancestry in Aboriginal Li People on Hainan Island

As the most prevalent aboriginal group on Hainan Island located between South China and the mainland of Southeast Asia, the Li people are believed to preserve some unique genetic information due to their isolated circumstances, although this has been largely uninvestigated. We performed the first whole-genome sequencing of 55 Hainan Li (HNL) individuals with high coverage (∼30-50×) to gain insight into their genetic history and potential adaptations. We identified the ancestry enriched in HNL (∼85%) is well preserved in present-day Tai-Kadai speakers residing in South China and North Vietnam, that is, Bai-Yue populations. A lack of admixture signature due to the geographical restriction exacerbated the bottleneck in the present-day HNL. The genetic divergence among Bai-Yue populations began ∼4,000-3,000 years ago when the proto-HNL underwent migration and the settling of Hainan Island. Finally, we identified signatures of positive selection in the HNL, some outstanding examples included FADS1 and FADS2 related to a diet rich in polyunsaturated fatty acids. In addition, we observed that malaria-driven selection had occurred in the HNL, with population-specific variants of malaria-related genes (e.g., CR1) present. Interestingly, HNL harbors a high prevalence of malaria leveraged gene variants related to hematopoietic function (e.g., CD3G) that may explain the high incidence of blood disorders such as B-cell lymphomas in the present-day HNL. The results have advanced our understanding of the genetic history of the Bai-Yue populations and have provided new insights into the adaptive scenarios of the Li people.

The contribution of Neanderthal introgression to modern human traits

Neanderthals, our closest extinct relatives, lived in western Eurasia from 400,000 years ago until they went extinct around 40,000 years ago. DNA retrieved from ancient specimens revealed that Neanderthals mated with modern human contemporaries. As a consequence, introgressed Neanderthal DNA survives scattered across the human genome such that 1-4% of the genome of present-day people outside Africa are inherited from Neanderthal ancestors. Patterns of Neanderthal introgressed genomic sequences suggest that Neanderthal alleles had distinct fates in the modern human genetic background. Some Neanderthal alleles facilitated human adaptation to new environments such as novel climate conditions, UV exposure levels and pathogens, while others had deleterious consequences. Here, we review the body of work on Neanderthal introgression over the past decade. We describe how evolutionary forces shaped the genomic landscape of Neanderthal introgression and highlight the impact of introgressed alleles on human biology and phenotypic variation.

Population structure, genomic diversity and demographic history of Komodo dragons inferred from whole-genome sequencing

Population and conservation genetics studies have greatly benefited from the development of new techniques and bioinformatic tools associated with next-generation sequencing. Analysis of extensive data sets from whole-genome sequencing of even a few individuals allows the detection of patterns of fine-scale population structure and detailed reconstruction of demographic dynamics through time. In this study, we investigated the population structure, genomic diversity and demographic history of the Komodo dragon (Varanus komodoensis), the world's largest lizard, by sequencing the whole genomes of 24 individuals from the five main Indonesian islands comprising the entire range of the species. Three main genomic groups were observed. The populations of the Island of Komodo and the northern coast of Flores, in particular, were identified as two distinct conservation units. Degrees of genomic divergence among island populations were interpreted as a result of changes in sea level affecting connectivity across islands. Demographic inference suggested that Komodo dragons probably experienced a relatively steep population decline over the last million years, reaching a relatively stable N_e during the Saalian glacial cycle (400-150 thousand years ago) followed by a rapid N_e decrease. Genomic diversity of Komodo dragons was similar to that found in endangered or already extinct reptile species. Overall, this study provides an example of how whole-genome analysis of a few individuals per population can help define population structure and intraspecific demographic dynamics. This is particularly important when applying population genomics data to conservation of rare or elusive endangered species.

Allele frequency differentiation at height-associated SNPs among continental human populations

Methods to detect polygenic adaptation have recently been shown to be sensitive to uncorrected stratification in GWAS, thereby casting doubts on whether polygenic adaptation is prevalent among humans. Consistent with a signal of adaptation at human height loci, the mean F_ST among African, East Asian, and European populations was shown to be significantly higher at height-associated SNPs than that at non-associated SNPs. This conclusion was reached, however, using height-associated SNPs ascertained from a GWAS design impacted by residual confounding due to uncorrected stratification. Specifically, we show here that the estimated effect sizes are significantly correlated with population structure across continents, potentially explaining the elevated differentiation previously reported. We alleviated these concerns of confounding by ascertaining height-associated SNPs from two biobank GWAS (UK Biobank, UKB, and Biobank Japan, BBJ), where measures to control for confounding in GWAS are more effective. Consistent with a global signature of polygenic adaptation, we found that compared to non-associated SNPs, frequencies of height-associated SNPs are indeed significantly more differentiated among continental populations from both the 1000 Genomes Project (p = 0.0012 for UKB and p = 0.0265 for BBJ), and the Human Genome Diversity Project (p = 0.0225 for UKB and p = 0.0032 for BBJ). However, we found no significant difference among continental populations in polygenic height scores. Through simulations, we found that polygenic score-based statistics could lose power in detecting polygenic adaptation in presence of independent converging selections, thereby potentially explaining the inconsistent results based on F_ST and polygenic scores.

Skeletal remains of a Pleistocene modern human (Homo sapiens) from Sulawesi

Major gaps remain in our knowledge of the early history of Homo sapiens in Wallacea. By 70-60 thousand years ago (ka), modern humans appear to have entered this distinct biogeographical zone between continental Asia and Australia. Despite this, there are relatively few Late Pleistocene sites attributed to our species in Wallacea. H. sapiens fossil remains are also rare. Previously, only one island in Wallacea (Alor in the southeastern part of the archipelago) had yielded skeletal evidence for pre-Holocene modern humans. Here we report on the first Pleistocene human skeletal remains from the largest Wallacean island, Sulawesi. The recovered elements consist of a nearly complete palate and frontal process of a modern human right maxilla excavated from Leang Bulu Bettue in the southwestern peninsula of the island. Dated by several different methods to between 25 and 16 ka, the maxilla belongs to an elderly individual of unknown age and sex, with small teeth (only M1 to M3 are extant) that exhibit severe occlusal wear and related dental pathologies. The dental wear pattern is unusual. This fragmentary specimen, though largely undiagnostic with regards to morphological affinity, provides the only direct insight we currently have from the fossil record into the identity of the Late Pleistocene people of Sulawesi.

Philippine Ayta possess the highest level of Denisovan ancestry in the world

Multiple lines of evidence show that modern humans interbred with archaic Denisovans. Here, we report an account of shared demographic history between Australasians and Denisovans distinctively in Island Southeast Asia. Our analyses are based on ∼2.3 million genotypes from 118 ethnic groups of the Philippines, including 25 diverse self-identified Negrito populations, along with high-coverage genomes of Australopapuans and Ayta Magbukon Negritos. We show that Ayta Magbukon possess the highest level of Denisovan ancestry in the world-∼30%-40% greater than that of Australians and Papuans-consistent with an independent admixture event into Negritos from Denisovans. Together with the recently described Homo luzonensis, we suggest that there were multiple archaic species that inhabited the Philippines prior to the arrival of modern humans and that these archaic groups may have been genetically related. Altogether, our findings unveil a complex intertwined history of modern and archaic humans in the Asia-Pacific region, where distinct Islander Denisovan populations differentially admixed with incoming Australasians across multiple locations and at various points in time.

Genome structural variation in human evolution

Structural variation (SV) is a large difference (typically >100 bp) in the genomic structure of two genomes and includes both copy number variation and variation that does not change copy number of a genomic region, such as an inversion. Improved reference genomes, combined with widespread genome sequencing using short-read sequencing technology, and increasingly using long-read sequencing, have reignited interest in SV. Recent large-scale studies and functional focused analyses have highlighted the role of SV in human evolution. In this review, we highlight human-specific SVs involved in changes in the brain, population-specific SVs that affect response to the environment, including adaptation to diet and infectious diseases, and summarise the contribution of archaic hominin admixture to present-day human SV.

Our Tangled Family Tree: New Genomic Methods Offer Insight into the Legacy of Archaic Admixture

The archaic ancestry present in the human genome has captured the imagination of both scientists and the wider public in recent years. This excitement is the result of new studies pushing the envelope of what we can learn from the archaic genetic information that has survived for over 50,000 years in the human genome. Here, we review the most recent ten years of literature on the topic of archaic introgression, including the current state of knowledge on Neanderthal and Denisovan introgression, as well as introgression from other as-yet unidentified archaic populations. We focus this review on four topics: 1) a reimagining of human demographic history, including evidence for multiple admixture events between modern humans, Neanderthals, Denisovans, and other archaic populations; 2) state-of-the-art methods for detecting archaic ancestry in population-level genomic data; 3) how these novel methods can detect archaic introgression in modern African populations; and 4) the functional consequences of archaic gene variants, including how those variants were co-opted into novel function in modern human populations. The goal of this review is to provide a simple-to-access reference for the relevant methods and novel data, which has changed our understanding of the relationship between our species and its siblings. This body of literature reveals the large degree to which the genetic legacy of these extinct hominins has been integrated into the human populations of today.

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.

Anthropometric, biochemical and hormonal profiles of the partially admixed pygmoid group in Rampasasa (Flores, Indonesia)

OBJECTIVES

We performed a cross-sectional study on anthropometric and laboratory characteristics of inhabitants of Rampasasa (Flores, Indonesia). Adults were categorised according to ancestry into three groups: pygmoid (P/P, offspring of pygmoid parents, n=8), mixed pygmoid (P/N, offspring of pygmoid and non-pygmoid parents, n=12) and non-pygmoid (N/N, n=10). Children (n=28) were P/N.

METHODS

Measurements included height, weight, sitting height, arm span, head circumference, haematological analysis and serum albumin, calcium, vitamin D, insulin-like growth factor-I (IGF-I) and IGF binding protein 3 (IGFBP-3). Pubertal stage and bone age was assessed in children. Anthropometric data were expressed as standard deviation score (SDS) for age. IGF-I, IGFBP-3 and IGF-I/IGFBP-3 ratio were expressed as SDS for age, bone age and pubertal stage.

RESULTS

Mean height SDS showed a gradient from P/P (-4.0) via P/N (-3.2) to N/N (-2.3) (-3.4, -3.1 and -2.2 adjusted for age-associated shrinking). Sitting height and head circumference showed similar gradients. Serum IGF-I SDS was similar among groups (approximately -1 SDS). IGFBP-3 SDS tended toward a gradient from P/P (-1.9) via P/N (-1.5) to N/N (-1.1), but IGF-I/IGFBP-3 ratio was normal in all groups. In P/P and P/N, mean head circumference SDS was >2 SD greater than mean height SDS. Children showed a progressive growth failure and bone age delay, delayed female pubertal onset and an initial low serum IGF-I, normal IGFBP-3 and low IGF-I/IGFBP-3 ratio.

CONCLUSIONS

P/P showed proportionate short stature with relative macrocephaly and relatively low IGFBP-3; P/N presented an intermediate pattern. P/N children were progressively short, showed delayed skeletal maturation, delayed puberty in girls and low IGF-I and IGF-I/IGFBP-3.

Widespread Denisovan ancestry in Island Southeast Asia but no evidence of substantial super-archaic hominin admixture

The hominin fossil record of Island Southeast Asia (ISEA) indicates that at least two endemic 'super-archaic' species-Homo luzonensis and H. floresiensis-were present around the time anatomically modern humans arrived in the region >50,000 years ago. Intriguingly, contemporary human populations across ISEA carry distinct genomic traces of ancient interbreeding events with Denisovans-a separate hominin lineage that currently lacks a fossil record in ISEA. To query this apparent disparity between fossil and genetic evidence, we performed a comprehensive search for super-archaic introgression in >400 modern human genomes, including >200 from ISEA. Our results corroborate widespread Denisovan ancestry in ISEA populations, but fail to detect any substantial super-archaic admixture signals compatible with the endemic fossil record of ISEA. We discuss the implications of our findings for the understanding of hominin history in ISEA, including future research directions that might help to unlock more details about the prehistory of the enigmatic Denisovans.

Associations of Arachidonic Acid Synthesis with Cardiovascular Risk Factors and Relation to Ischemic Heart Disease and Stroke: A Univariable and Multivariable Mendelian Randomization Study

Arachidonic acid (AA), a major long-chain omega-6 polyunsaturated fatty acid, is associated with ischemic heart disease (IHD) and stroke. We assessed bi-directional associations of AA synthesis reflected by plasma phospholipid AA with CVD risk factors, and identified mediators of associations of AA with IHD and stroke using Mendelian randomization (MR). We used two-sample MR to assess bi-directional associations of AA synthesis with lipids, blood pressure, adiposity, and markers of inflammation and coagulation. We used multivariable MR to assess mediators of associations of AA with IHD and stroke. Genetically predicted AA (% of total fatty acids increase) was positively associated with apolipoprotein B (ApoB, 0.022 standard deviations (SD), 95% confidence interval (CI) 0.010, 0.034), high-density (0.030 SD, 95% CI 0.012, 0.049) and low-density lipoprotein cholesterol (LDL-C, 0.016 SD, 95% CI 0.004, 0.027) and lower triglycerides (-0.031 SD, 95% CI -0.049, -0.012) but not with other traits. Genetically predicted these traits gave no association with AA. The association of AA with IHD was attenuated adjusting for ApoB or LDL-C. Genetically predicted AA was associated with lipids but not other traits. Given ApoB is thought to be the key lipid in IHD, the association of AA with IHD is likely mediated by ApoB.

The distinct morphological phenotypes of Southeast Asian aborigines are shaped by novel mechanisms for adaptation to tropical rainforests

Southeast Asian aborigines, the hunter-gatherer populations living in tropical rainforests, exhibit distinct morphological phenotypes, including short stature, dark skin, curly hair and a wide and snub nose. The underlying genetic architecture and evolutionary mechanism of these phenotypes remain a long-term mystery. We conducted whole genome deep sequencing of 81 Cambodian aborigines from eight ethnic groups. Through a genome-wide scan of selective sweeps, we discovered key genes harboring Cambodian-enriched mutations that may contribute to their phenotypes, including two hair morphogenesis genes (TCHH and TCHHL1), one nasal morphology gene (PAX3) and a set of genes (such as ENTPD1-AS1) associated with short stature. The identified new genes and novel mutations suggest an independent origin of the distinct phenotypes in Cambodian aborigines through parallel evolution, refuting the long-standing argument on the common ancestry of these phenotypes among the worldwide rainforest hunter-gatherers. Notably, our discovery reveals that various types of molecular mechanisms, including antisense transcription and epigenetic regulation, contribute to human morphogenesis, providing novel insights into the genetics of human environmental adaptation.

Origins of modern human ancestry

New finds in the palaeoanthropological and genomic records have changed our view of the origins of modern human ancestry. Here we review our current understanding of how the ancestry of modern humans around the globe can be traced into the deep past, and which ancestors it passes through during our journey back in time. We identify three key phases that are surrounded by major questions, and which will be at the frontiers of future research. The most recent phase comprises the worldwide expansion of modern humans between 40 and 60 thousand years ago (ka) and their last known contacts with archaic groups such as Neanderthals and Denisovans. The second phase is associated with a broadly construed African origin of modern human diversity between 60 and 300 ka. The oldest phase comprises the complex separation of modern human ancestors from archaic human groups from 0.3 to 1 million years ago. We argue that no specific point in time can currently be identified at which modern human ancestry was confined to a limited birthplace, and that patterns of the first appearance of anatomical or behavioural traits that are used to define Homo sapiens are consistent with a range of evolutionary histories.

Evidence of Polygenic Adaptation in Sardinia at Height-Associated Loci Ascertained from the Biobank Japan

Adult height is one of the earliest putative examples of polygenic adaptation in humans. However, this conclusion was recently challenged because residual uncorrected stratification from large-scale consortium studies was considered responsible for the previously noted genetic difference. It thus remains an open question whether height loci exhibit signals of polygenic adaptation in any human population. We re-examined this question, focusing on one of the shortest European populations, the Sardinians, in addition to mainland European populations. We utilized height-associated loci from the Biobank Japan (BBJ) dataset to further alleviate concerns of biased ascertainment of GWAS loci and showed that the Sardinians remain significantly shorter than expected under neutrality (∼0.22 standard deviation shorter than Utah residents with ancestry from northern and western Europe [CEU] on the basis of polygenic height scores, p = 3.89 × 10-4). We also found the trajectory of polygenic height scores between the Sardinian and the British populations diverged over at least the last 10,000 years (p = 0.0082), consistent with a signature of polygenic adaptation driven primarily by the Sardinian population. Although the polygenic score-based analysis showed a much subtler signature in mainland European populations, we found a clear and robust adaptive signature in the UK population by using a haplotype-based statistic, the trait singleton density score (tSDS), driven by the height-increasing alleles (p = 9.1 × 10-4). In summary, by ascertaining height loci in a distant East Asian population, we further supported the evidence of polygenic adaptation at height-associated loci among the Sardinians. In mainland Europeans, the adaptive signature was detected in haplotype-based analysis but not in polygenic score-based analysis.

Methods for detecting introgressed archaic sequences

Analysis of genome sequences from archaic and modern humans have revealed multiple episodes of admixture between highly-diverged population groups. Statistical methods that attempt to localize DNA segments introduced by these events offer a powerful tool to investigate recent human evolution. We review recent advances in methods for detecting introgressed sequences.

Human adaptation over the past 40,000 years

Over the past few years several methodological and data-driven advances have greatly improved our ability to robustly detect genomic signatures of selection in humans. New methods applied to large samples of present-day genomes provide increased power, while ancient DNA allows precise estimation of timing and tempo. However, despite these advances, we are still limited in our ability to translate these signatures into understanding about which traits were actually under selection, and why. Combining information from different populations and timescales may allow interpretation of selective sweeps. Other modes of selection have proved more difficult to detect. In particular, despite strong evidence of the polygenicity of most human traits, evidence for polygenic selection is weak, and its importance in recent human evolution remains unclear. Balancing selection and archaic introgression seem important for the maintenance of potentially adaptive immune diversity, but perhaps less so for other traits.

Evolutionary perspectives on polygenic selection, missing heritability, and GWAS

Genome-wide association studies (GWAS) have successfully identified many trait-associated variants, but there is still much we do not know about the genetic basis of complex traits. Here, we review recent theoretical and empirical literature regarding selection on complex traits to argue that "missing heritability" is as much an evolutionary problem as it is a statistical problem. We discuss empirical findings that suggest a role for selection in shaping the effect sizes and allele frequencies of causal variation underlying complex traits, and the limitations of these studies. We then use simulations of selection, realistic genome structure, and complex human demography to illustrate the results of recent theoretical work on polygenic selection, and show that statistical inference of causal loci is sharply affected by evolutionary processes. In particular, when selection acts on causal alleles, it hampers the ability to detect causal loci and constrains the transferability of GWAS results across populations. Last, we discuss the implications of these findings for future association studies, and suggest that future statistical methods to infer causal loci for genetic traits will benefit from explicit modeling of the joint distribution of effect sizes and allele frequencies under plausible evolutionary models.

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.

Somewhere beyond the sea: Human cranial remains from the Lesser Sunda Islands (Alor Island, Indonesia) provide insights on Late Pleistocene peopling of Island Southeast Asia

The migration of anatomically modern humans (AMH) from Africa to every inhabitable continent included their dispersal through Island Southeast Asia (ISEA) to Australia. Significantly, this involved overwater dispersal through the Lesser Sunda Islands between Sunda (continental Southeast Asia) and Sahul (Australia and New Guinea). However, the timing and direction of this movement is still debated. Here, we report on human skeletal material recovered from excavations at two rockshelters, known locally as Tron Bon Lei, on Alor Island, Indonesia. The remains, dated to the Late Pleistocene, are the first anatomically modern human remains recovered in Wallacea dated to this period and are associated with cultural material demonstrating intentional burial. The human remains from Tron Bon Lei represent a population osteometrically distinct from Late Pleistocene Sunda and Sahul AMH. Instead, morphometrically, they appear more similar to Holocene populations in the Lesser Sundas. Thus, they may represent the remains of a population originally from Sunda whose Lesser Sunda Island descendants survived into the Holocene.

Population genomics perspectives on convergent adaptation

Convergent adaptation is the independent evolution of similar traits conferring a fitness advantage in two or more lineages. Cases of convergent adaptation inform our ideas about the ecological and molecular basis of adaptation. In judging the degree to which putative cases of convergent adaptation provide an independent replication of the process of adaptation, it is necessary to establish the degree to which the evolutionary change is unexpected under null models and to show that selection has repeatedly, independently driven these changes. Here, we discuss the issues that arise from these questions particularly for closely related populations, where gene flow and standing variation add additional layers of complexity. We outline a conceptual framework to guide intuition as to the extent to which evolutionary change represents the independent gain of information owing to selection and show that this is a measure of how surprised we should be by convergence. Additionally, we summarize the ways population and quantitative genetics and genomics may help us address questions related to convergent adaptation, as well as open new questions and avenues of research. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

Genetic variants influencing obesity-related traits in Japanese population

Context: Adipose tissue is the main organ that stores energy and participates in adaptive thermogenesis of the human body. The adipose tissue content in an individual is determined by a combination of genetic factors and lifestyle related factors. While Japanese people, along with the closely related East Asians, are generally thinner than individuals of European ancestry, they are prone to accumulating visceral adipose tissues. Genome-wide discovery of loci influencing obesity-related traits, and application of the genome sequence data to assess natural selection, provides evidence that the obesity-related traits in East Asians might be shaped by natural selection. Objective: This review aims to summarise health and evolutionary implications of genetic variants influencing obesity-related traits in Japanese. Methods: This study gathered recently published papers of medical, genetic and evolutionary studies regarding obesity-related traits in the Japanese and closely related East Asians. Results and conclusion: A high susceptibility to central obesity of Japanese and closely related East Asians might have been shaped by natural selection favouring thrifty genotypes. Moreover, natural selection favouring higher thermogenic activity of brown adipose tissues would contribute to increased non-thrifty alleles in ancestors of East Asians.

Using hominin introgression to trace modern human dispersals

The dispersal of anatomically modern human populations out of Africa and across much of the rest of the world around 55 to 50 thousand years before present (ka) is recorded genetically by the multiple hominin groups they met and interbred with along the way, including the Neandertals and Denisovans. The signatures of these introgression events remain preserved in the genomes of modern-day populations, and provide a powerful record of the sequence and timing of these early migrations, with Asia proving a particularly complex area. At least 3 different hominin groups appear to have been involved in Asia, of which only the Denisovans are currently known. Several interbreeding events are inferred to have taken place east of Wallace's Line, consistent with archaeological evidence of widespread and early hominin presence in the area. However, archaeological and fossil evidence indicates archaic hominins had not spread as far as the Sahul continent (New Guinea, Australia, and Tasmania), where recent genetic evidence remains enigmatic.