Positive Selection In Humans Research Articles

Rare variation in non-coding regions with evolutionary signatures contributes to autism spectrum disorder risk

Little is known about the role of non-coding regions in the etiology of autism spectrum disorder (ASD). We examined three classes of non-coding regions: human accelerated regions (HARs), which show signatures of positive selection in humans; experimentally validated neural VISTA enhancers (VEs); and conserved regions predicted to act as neural enhancers (CNEs). Targeted and whole-genome analysis of >16,600 samples and >4,900 ASD probands revealed that likely recessive, rare, inherited variants in HARs, VEs, and CNEs substantially contribute to ASD risk in probands whose parents share ancestry, which enriches for recessive contributions, but modestly contribute, if at all, in simplex family structures. We identified multiple patient variants in HARs near IL1RAPL1 and in VEs near OTX1 and SIM1 and showed that they change enhancer activity. Our results implicate both human-evolved and evolutionarily conserved non-coding regions in ASD risk and suggest potential mechanisms of how regulatory changes can modulate social behavior.

Human genetic adaptation related to cellular zinc homeostasis.

SLC30A9 encodes a ubiquitously zinc transporter (ZnT9) and has been consistently suggested as a candidate for positive selection in humans. However, no direct adaptive molecular phenotype has been demonstrated. Our results provide evidence for directional selection operating in two major complementary haplotypes in Africa and East Asia. These haplotypes are associated with differential gene expression but also differ in the Met50Val substitution (rs1047626) in ZnT9, which we show is found in homozygosis in the Denisovan genome and displays accompanying signatures suggestive of archaic introgression. Although we found no significant differences in systemic zinc content between individuals with different rs1047626 genotypes, we demonstrate that the expression of the derived isoform (ZnT9 50Val) in HEK293 cells shows a gain of function when compared with the ancestral (ZnT9 50Met) variant. Notably, the ZnT9 50Val variant was found associated with differences in zinc handling by the mitochondria and endoplasmic reticulum, with an impact on mitochondrial metabolism. Given the essential role of the mitochondria in skeletal muscle and since the derived allele at rs1047626 is known to be associated with greater susceptibility to several neuropsychiatric traits, we propose that adaptation to cold may have driven this selection event, while also impacting predisposition to neuropsychiatric disorders in modern humans.

Evolutionary relevance of single nucleotide variants within the forebrain exclusive human accelerated enhancer regions

BackgroundHuman accelerated regions (HARs) are short conserved genomic sequences that have acquired significantly more nucleotide substitutions than expected in the human lineage after divergence from chimpanzees. The fast evolution of HARs may reflect their roles in the origin of human-specific traits. A recent study has reported positively-selected single nucleotide variants (SNVs) within brain-exclusive human accelerated enhancers (BE-HAEs) hs1210 (forebrain), hs563 (hindbrain) and hs304 (midbrain/forebrain). By including data from archaic hominins, these SNVs were shown to be Homo sapiens-specific, residing within transcriptional factors binding sites (TFBSs) for SOX2 (hs1210), RUNX1/3 (hs563), and FOS/JUND (hs304). Although these findings suggest that the predicted modifications in TFBSs may have some role in present-day brain structure, work is required to verify the extent to which these changes translate into functional variation.ResultsTo start to fill this gap, we investigate the SOX2 SNV, with both forebrain expression and strong signal of positive selection in humans. We demonstrate that the HMG box of SOX2 binds in vitro with Homo sapiens-specific derived A-allele and ancestral T-allele carrying DNA sites in BE-HAE hs1210. Molecular docking and simulation analysis indicated highly favourable binding of HMG box with derived A-allele containing DNA site when compared to site carrying ancestral T-allele.ConclusionThese results suggest that adoptive changes in TF affinity within BE-HAE hs1210 and other HAR enhancers in the evolutionary history of Homo sapiens might have brought about changes in gene expression patterns and have functional consequences on forebrain formation and evolution.MethodsThe present study employ electrophoretic mobility shift assays (EMSA) and molecular docking and molecular dynamics simulations approaches.

Localizing Post-Admixture Adaptive Variants with Object Detection on Ancestry-Painted Chromosomes.

Gene flow between previously differentiated populations during the founding of an admixed or hybrid population has the potential to introduce adaptive alleles into the new population. If the adaptive allele is common in one source population, but not the other, then as the adaptive allele rises in frequency in the admixed population, genetic ancestry from the source containing the adaptive allele will increase nearby as well. Patterns of genetic ancestry have therefore been used to identify post-admixture positive selection in humans and other animals, including examples in immunity, metabolism, and animal coloration. A common method identifies regions of the genome that have local ancestry "outliers" compared with the distribution across the rest of the genome, considering each locus independently. However, we lack theoretical models for expected distributions of ancestry under various demographic scenarios, resulting in potential false positives and false negatives. Further, ancestry patterns between distant sites are often not independent. As a result, current methods tend to infer wide genomic regions containing many genes as under selection, limiting biological interpretation. Instead, we develop a deep learning object detection method applied to images generated from local ancestry-painted genomes. This approach preserves information from the surrounding genomic context and avoids potential pitfalls of user-defined summary statistics. We find the method is robust to a variety of demographic misspecifications using simulated data. Applied to human genotype data from Cabo Verde, we localize a known adaptive locus to a single narrow region compared with multiple or long windows obtained using two other ancestry-based methods.

Differentially Active and Conserved Neural Enhancers Define Two Forms of Adaptive Noncoding Evolution in Humans.

The human and chimpanzee genomes are strikingly similar, but our neural phenotypes are very different. Many of these differences are likely driven by changes in gene expression, and some of those changes may have been adaptive during human evolution. Yet, the relative contributions of positive selection on regulatory regions or other functional regulatory changes are unclear. Where are these changes located throughout the human genome? Are functional regulatory changes near genes or are they in distal enhancer regions? In this study, we experimentally combined both human and chimpanzee cis-regulatory elements (CREs) that showed either (1) signs of accelerated evolution in humans or (2) that have been shown to be active in the human brain. Using a massively parallel reporter assay, we tested the ability of orthologous human and chimpanzee CREs to activate transcription in induced pluripotent stem-cell-derived neural progenitor cells and neurons. With this assay, we identified 179 CREs with differential activity between human and chimpanzee; in contrast, we found 722 CREs with signs of positive selection in humans. Selection and differentially expressed CREs strikingly differ in level of expression, size, and genomic location. We found a subset of 69 CREs in loci with genetic variants associated with neuropsychiatric diseases, which underscores the consequence of regulatory activity in these loci for proper neural development and function. By combining CREs that either experienced recent selection in humans or CREs that are functional brain enhancers, presents a novel way of studying the evolution of noncoding elements that contribute to human neural phenotypes.

A Positively Selected MAGEE2 LoF Allele Is Associated with Sexual Dimorphism in Human Brain Size and Shows Similar Phenotypes in Magee2 Null Mice.

A nonsense allele at rs1343879 in human MAGEE2 on chromosome X has previously been reported as a strong candidate for positive selection in East Asia. This premature stop codon causing ∼80% protein truncation is characterized by a striking geographical pattern of high population differentiation: common in Asia and the Americas (up to 84% in the 1000 Genomes Project East Asians) but rare elsewhere. Here, we generated a Magee2 mouse knockout mimicking the human loss-of-function mutation to study its functional consequences. The Magee2 null mice did not exhibit gross abnormalities apart from enlarged brain structures (13% increased total brain area, P = 0.0022) in hemizygous males. The area of the granular retrosplenial cortex responsible for memory, navigation, and spatial information processing was the most severely affected, exhibiting an enlargement of 34% (P = 3.4×10−6). The brain size in homozygous females showed the opposite trend of reduced brain size, although this did not reach statistical significance. With these insights, we performed human association analyses between brain size measurements and rs1343879 genotypes in 141 Chinese volunteers with brain MRI scans, replicating the sexual dimorphism seen in the knockout mouse model. The derived stop gain allele was significantly associated with a larger volume of gray matter in males (P = 0.00094), and smaller volumes of gray (P = 0.00021) and white (P = 0.0015) matter in females. It is unclear whether or not the observed neuroanatomical phenotypes affect behavior or cognition, but it might have been the driving force underlying the positive selection in humans.

Transmission of SARS-CoV-2 in domestic cats imposes a narrow bottleneck.

The evolutionary mechanisms by which SARS-CoV-2 viruses adapt to mammalian hosts and, potentially, undergo antigenic evolution depend on the ways genetic variation is generated and selected within and between individual hosts. Using domestic cats as a model, we show that SARS-CoV-2 consensus sequences remain largely unchanged over time within hosts, while dynamic sub-consensus diversity reveals processes of genetic drift and weak purifying selection. We further identify a notable variant at amino acid position 655 in Spike (H655Y), which was previously shown to confer escape from human monoclonal antibodies. This variant arises rapidly and persists at intermediate frequencies in index cats. It also becomes fixed following transmission in two of three pairs. These dynamics suggest this site may be under positive selection in this system and illustrate how a variant can quickly arise and become fixed in parallel across multiple transmission pairs. Transmission of SARS-CoV-2 in cats involved a narrow bottleneck, with new infections founded by fewer than ten viruses. In RNA virus evolution, stochastic processes like narrow transmission bottlenecks and genetic drift typically act to constrain the overall pace of adaptive evolution. Our data suggest that here, positive selection in index cats followed by a narrow transmission bottleneck may have instead accelerated the fixation of S H655Y, a potentially beneficial SARS-CoV-2 variant. Overall, our study suggests species- and context-specific adaptations are likely to continue to emerge. This underscores the importance of continued genomic surveillance for new SARS-CoV-2 variants as well as heightened scrutiny for signatures of SARS-CoV-2 positive selection in humans and mammalian model systems.

The Role of Human Evolution In The Development of Schizophrenia And Cognitive Traits Applying Polygenic Analytical Approaches

Psychiatric disorders have like other complex human phenotypes a polygenic architecture and are unlikely to undergo selection in the form of major selective sweeps. It is more likely to develop due to subtle shifts in allele frequencies which are hard to detect using traditional statistical methods. Our recently developed false discovery rate (FDR) framework exploits ancillary information to improve the power of genome-wide association studies (GWAS). This method can assess enrichment of polygenic effects on various phenotypes in specific genic annotation categories. In a series of studies, we applied this method to investigate the effect of evolutionary metrics using data from GWASs of psychiatric and cognitive phenotypes, and evolutionary proxies representing various time periods of human evolution.The Neanderthal selective sweep score, which measures positive selection in humans after the human-Neanderthal split approximately 50,000 years ago, was used to test the hypothesis that the schizophrenia is related to the development of the complex human brain and concurrently, language skills and creative reasoning. We show that regions that underwent positive selection in humans after divergence from Neanderthals are more likely to be associated with schizophrenia.We also investigated if Neanderthal selective sweep score are enriched in cognitive phenotypes including education attainment, college completion, and more specific cognitive functions including reaction time and verbal numerical reasoning. The findings indicate that SNPs in the swept regions are more likely to be associated with education attainment, college completion, and to some extent to general cognitive function.Less recently evolved regions like the human accelerated regions (HAR), segmental duplication regions (SD), and ohnologous regions (Ohno) were analyzed for association with schizophrenia. These represent time periods before 200,000 years ago when hominid and chimpanzee branches split (HAR), 35-40 million years ago during the great ape-hominid split (SD), and the period of separation of vertebrate and invertebrate lineages 500 million years ago (Ohno). We could not detect significant enrichment for any of these evolutionary markers.These findings help us understand the evolutionary history of schizophrenia and cognitive phenotypes. As humans evolved, our ability to think in abstract terms, use language to communicate with one another and organize complex rule-based societies improved, but such development probably made us more vulnerable to mental pathologies like schizophrenia.

Selection Pressures on RNA Sequences and Structures.

With the discovery of increasingly more functional noncoding RNAs (ncRNAs), it becomes eminent to more strongly consider them as important players during species evolution. Although tests for negative selection of ncRNAs already exist since the beginning of this century, the SSS-test is the first one for also investigating positive selection. When analyzing selection in ncRNAs, it should be taken into account that selection pressures can independently act on sequence and structure. We applied the SSS-test to explore the evolution of ncRNAs in primates and identified more than 100 long noncoding RNAs (lncRNAs) that might evolve under positive selection in humans. With this test, it is now possible to more thoroughly include ncRNAs into evolutionary studies.

Enrichment of genetic markers of recent human evolution in educational and cognitive traits

Higher cognitive functions are regarded as one of the main distinctive traits of humans. Evidence for the cognitive evolution of human beings is mainly based on fossil records of an expanding cranium and an increasing complexity of material culture artefacts. However, the molecular genetic factors involved in the evolution are still relatively unexplored. Here, we investigated whether genomic regions that underwent positive selection in humans after divergence from Neanderthals are enriched for genetic association with phenotypes related to cognitive functions. We used genome wide association data from a study of college completion (N = 111,114), one of educational attainment (N = 293,623) and two different studies of general cognitive ability (N = 269,867 and 53,949). We found nominally significant polygenic enrichment of associations with college completion (p = 0.025), educational attainment (p = 0.043) and general cognitive ability (p = 0.015 and 0.025, respectively), suggesting that variants influencing these phenotypes are more prevalent in evolutionarily salient regions. The enrichment remained significant after controlling for other known genetic enrichment factors, and for affiliation to genes highly expressed in the brain. These findings support the notion that phenotypes related to higher order cognitive skills typical of humans have a recent genetic component that originated after the separation of the human and Neanderthal lineages.

Abstract 1028: Patient-derived xenografts undergo mouse-specific tumor evolution

Abstract Patient-derived xenografts (PDXs) have become a prominent cancer model system, as they are presumed to faithfully represent the genomic features of primary tumors. Here we monitored the dynamics of copy number alterations (CNAs) in 1,110 PDX samples across 24 cancer types. We observed rapid accumulation of CNAs during PDX passaging, often due to selection of pre-existing minor clones. CNA acquisition in PDXs was correlated with the tissue-specific levels of aneuploidy and genetic heterogeneity observed in primary tumors. However, the particular CNAs acquired during PDX passaging differed from those acquired during tumor evolution in patients. Several CNAs recurrently observed in primary tumors gradually disappeared in PDXs, indicating that events undergoing positive selection in humans can become dispensable during propagation in mice. Importantly, the genomic stability of PDXs was associated with their response to chemotherapy and targeted drugs. These findings have important implications for PDX-based modeling of human cancer. Citation Format: Uri Ben-David, Gavin Ha, Yuen-Yi Tseng, Noah F. Greenwald, Coyin Oh, Juliann Shih, James M. McFarland, Bang Wong, Jesse S. Boehm, Rameen Beroukhim, Todd R. Golub. Patient-derived xenografts undergo mouse-specific tumor evolution [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 1028.

Human-Specific Mutations and Positively Selected Sites in MARCO Confer Functional Changes.

Macrophage Receptor with COllagenous structure (MARCO) is a class A scavenger receptor that binds, phagocytoses, and modifies inflammatory responses to bacterial pathogens. Multiple candidate gene approach studies have shown that polymorphisms in MARCO are associated with susceptibility or resistance to Mycobacterium tuberculosis infection, but how these variants alter function is not known. To complement candidate gene approach studies, we previously used phylogenetic analyses to identify a residue, glutamine 452 (Q452), within the ligand-binding Scavenger Receptor Cysteine Rich domain as undergoing positive selection in humans. Herein, we show that Q452 is found in Denisovans, Neanderthals, and extant humans, but all other nonprimate, terrestrial, and aquatic mammals possess an aspartic acid (D452) residue. Further analysis of hominoid sequences of MARCO identified an additional human-specific mutation, phenylalanine 282 (F282), within the collagenous domain. We show that residue 282 is polymorphic in humans, but only 17% of individuals (rs6761637) possess the ancestral serine residue at position 282. We show that rs6761637 is in linkage disequilibrium with MARCO polymorphisms that have been previously linked to susceptibility to pulmonary tuberculosis. To assess the functional importance of sites Q452 and F282 in humans, we cloned the ancestral residues and loss-of-function mutations and investigated the role of these residues in binding and internalizing polystyrene microspheres and Escherichia coli. Herein, we show that the residues at sites 452 and 282 enhance receptor function.

Patient-derived xenografts undergo mouse-specific tumor evolution.

Patient-derived xenografts (PDXs) have become a prominent cancer model system, as they are presumed to faithfully represent the genomic features of primary tumors. Here we monitored the dynamics of copy number alterations (CNAs) in 1,110 PDX samples across 24 cancer types. We observed rapid accumulation of CNAs during PDX passaging, often due to selection of pre-existing minor clones. CNA acquisition in PDXs was correlated with the tissue-specific levels of aneuploidy and genetic heterogeneity observed in primary tumors. However, the particular CNAs acquired during PDX passaging differed from those acquired during tumor evolution in patients. Several CNAs recurrently observed in primary tumors gradually disappeared in PDXs, indicating that events undergoing positive selection in humans can become dispensable during propagation in mice. Importantly, the genomic stability of PDXs was associated with their response to chemotherapy and targeted drugs. These findings have important implications for PDX-based modeling of human cancer.

REST, a master regulator of neurogenesis, evolved under strong positive selection in humans and in non human primates

The transcriptional repressor REST regulates many neuronal genes by binding RE1 motifs. About one third of human RE1s are recently evolved and specific to primates. As changes in the activity of a transcription factor reverberate on its downstream targets, we assessed whether REST displays fast evolutionary rates in primates. We show that REST was targeted by very strong positive selection during primate evolution. Positive selection was also evident in the human lineage, with six selected sites located in a region that surrounds a VNTR in exon 4. Analysis of expression data indicated that REST brain expression peaks during aging in humans but not in other primates. Because a REST coding variant (rs3796529) was previously associated with protection from hippocampal atrophy in elderly subjects with mild cognitive impairment (MCI), we analyzed a cohort of Alzheimer disease (AD) continuum patients. Genotyping of two coding variants (rs3796529 and rs2227902) located in the region surrounding the VNTR indicated a role for rs2227902 in modulation of hippocampal volume loss, indirectly confirming a role for REST in neuroprotection. Experimental studies will be instrumental to determine the functional effect of positively selected sites in REST and the role of REST variants in neuropreservation/neurodegeneration.

Detecting ancient positive selection in humans using extended lineage sorting.

Natural selection that affected modern humans early in their evolution has likely shaped some of the traits that set present-day humans apart from their closest extinct and living relatives. The ability to detect ancient natural selection in the human genome could provide insights into the molecular basis for these human-specific traits. Here, we introduce a method for detecting ancient selective sweeps by scanning for extended genomic regions where our closest extinct relatives, Neandertals and Denisovans, fall outside of the present-day human variation. Regions that are unusually long indicate the presence of lineages that reached fixation in the human population faster than expected under neutral evolution. Using simulations, we show that the method is able to detect ancient events of positive selection and that it can differentiate those from background selection. Applying our method to the 1000 Genomes data set, we find evidence for ancient selective sweeps favoring regulatory changes and present a list of genomic regions that are predicted to underlie positively selected human specific traits.

Ancient selection for derived alleles at a GDF5 enhancer influencing human growth and osteoarthritis risk.

Variants in GDF5 are associated with human arthritis and height reduction, but the causal mutations are still unknown. We surveyed GDF5 for regulatory regions in transgenic mice, and fine-mapped separate enhancers controlling expression in joints versus growing ends of long bones. A large downstream regulatory region harbors a novel growth enhancer (GROW1), which is required for normal GDF5 expression at ends of developing bones, and for normal bone lengths in vivo. Human GROW1 contains a common base pair change that reduces enhancer activity and colocalizes with peaks of positive selection in humans. The derived allele is rare in Africa, but common in Eurasia and found in Neandertals and Denisovans. Our studies suggest that an ancient regulatory variant in GROW1 has been repeatedly selected in northern environments, and that past selection on growth phenotypes explains the high frequency of a GDF5 haplotype that also increases arthritis susceptibility in many human populations.

Population genetic analysis of the DARC locus (Duffy) reveals adaptation from standing variation associated with malaria resistance in humans.

The human DARC (Duffy antigen receptor for chemokines) gene encodes a membrane-bound chemokine receptor crucial for the infection of red blood cells by Plasmodium vivax, a major causative agent of malaria. Of the three major allelic classes segregating in human populations, the FY*O allele has been shown to protect against P. vivax infection and is at near fixation in sub-Saharan Africa, while FY*B and FY*A are common in Europe and Asia, respectively. Due to the combination of strong geographic differentiation and association with malaria resistance, DARC is considered a canonical example of positive selection in humans. Despite this, details of the timing and mode of selection at DARC remain poorly understood. Here, we use sequencing data from over 1,000 individuals in twenty-one human populations, as well as ancient human genomes, to perform a fine-scale investigation of the evolutionary history of DARC. We estimate the time to most recent common ancestor (TMRCA) of the most common FY*O haplotype to be 42 kya (95% CI: 34–49 kya). We infer the FY*O null mutation swept to fixation in Africa from standing variation with very low initial frequency (0.1%) and a selection coefficient of 0.043 (95% CI:0.011–0.18), which is among the strongest estimated in the human genome. We estimate the TMRCA of the FY*A mutation in non-Africans to be 57 kya (95% CI: 48–65 kya) and infer that, prior to the sweep of FY*O, all three alleles were segregating in Africa, as highly diverged populations from Asia and ≠Khomani San hunter-gatherers share the same FY*A haplotypes. We test multiple models of admixture that may account for this observation and reject recent Asian or European admixture as the cause.

Ancient DNA Investigation of a Medieval German Cemetery Confirms Long-Term Stability of CCR5-Δ32 Allele Frequencies in Central Europe.

The CCR5-Δ32 mutation present in European populations is among the most prominently debated cases of recent positive selection in humans. This allele, a 32-bp deletion that renders the T-cell CCR5 receptor nonfunctional, has important epidemiological and public health significance, as homozygous carriers are resistant to several HIV strains. However, although the function of this allele in preventing HIV infection is now well described, its human evolutionary origin is poorly understood. Initial attempts to determine the emergence of the CCR5-Δ32 allele pointed to selection during the 14th-century Black Death pandemic; however, subsequent analyses suggest that the allele rose in frequency more than 5,000 years ago, possibly through drift. Recently, three studies have identified populations predating the 14th century CE that are positive for the CCR5-Δ32 allele, supporting the claim for a more ancient origin. However, these studies also suggest poorly understood regional differences in the recent evolutionary history of the CCR5-Δ32 allele. Here a new hydrolysis-probe-based real-time PCR assay was designed to ascertain CCR5 allele frequency in 53 individuals from a 10th- to 12th-century CE church and convent complex in central Germany that predates outbreaks of the Black Death pandemic. High-confidence genotypes were obtained for 32 individuals, and results show that CCR5-Δ32 allele frequency has remained unchanged in this region of Central Europe over the last millennium, suggesting that there has been no strong positive selective pressure over this time period and confirming a more ancient origin for the allele.

Changes in selective pressures associated with human population expansion may explain metabolic and immune related pathways enriched for signatures of positive selection.

BackgroundThe study of local adaptation processes is a very important research topic in the field of population genomics. There is a particular interest in the study of human populations because they underwent a process of rapid spatial expansion and faced important environmental changes that translated into changes in selective pressures. New mutations may have been selected for in the new environment and previously existing genetic variants may have become detrimental. Immune related genes may have been released from the selective pressure exerted by pathogens in the ancestral environment and new variants may have been positively selected due to pathogens present in the newly colonized habitat. Also, variants that had a selective advantage in past environments may have become deleterious in the modern world due to external stimuli including climatic, dietary and behavioral changes, which could explain the high prevalence of some polygenic diseases such as diabetes and obesity.ResultsWe performed an enrichment analysis to identify gene sets enriched for signals of positive selection in humans. We used two genome scan methods, XPCLR and iHS to detect selection using a dense coverage of SNP markers combined with two gene set enrichment approaches. We identified immune related gene sets that could be involved in the protection against pathogens especially in the African population. We also identified the glycolysis & gluconeogenesis gene set, related to metabolism, which supports the thrifty genotype hypothesis invoked to explain the current high prevalence of diseases such as diabetes and obesity. Extending our analysis to the gene level, we found signals for 23 candidate genes linked to metabolic syndrome, 13 of which are new candidates for positive selection.ConclusionsOur study provides a list of genes and gene sets associated with immunity and metabolic syndrome that are enriched for signals of positive selection in three human populations (Europeans, Africans and Asians). Our results highlight differences in the relative importance of pathogens as drivers of local adaptation in different continents and provide new insights into the evolution and high incidence of metabolic syndrome in modern human populations.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-2783-2) contains supplementary material, which is available to authorized users.

Genetic Evidence of Human Adaptation to a Cooked Diet

Humans have been argued to be biologically adapted to a cooked diet, but this hypothesis has not been tested at the molecular level. Here, we combine controlled feeding experiments in mice with comparative primate genomics to show that consumption of a cooked diet influences gene expression and that affected genes bear signals of positive selection in the human lineage. Liver gene expression profiles in mice fed standardized diets of meat or tuber were affected by food type and cooking, but not by caloric intake or consumer energy balance. Genes affected by cooking were highly correlated with genes known to be differentially expressed in liver between humans and other primates, and more genes in this overlap set show signals of positive selection in humans than would be expected by chance. Sequence changes in the genes under selection appear before the split between modern humans and two archaic human groups, Neandertals and Denisovans, supporting the idea that human adaptation to a cooked diet had begun by at least 275,000 years ago.