Applying an evolutionary mismatch framework to understand disease susceptibility.

Noncommunicable diseases (NCDs) are on the rise worldwide. Obesity, cardiovascular disease, and type 2 diabetes are among a long list of "lifestyle" diseases that were rare throughout human history but are now common. The evolutionary mismatch hypothesis posits that humans evolved in environments that radically differ from those we currently experience; consequently, traits that were once advantageous may now be "mismatched" and disease causing. At the genetic level, this hypothesis predicts that loci with a history of selection will exhibit "genotype by environment" (GxE) interactions, with different health effects in "ancestral" versus "modern" environments. To identify such loci, we advocate for combining genomic tools in partnership with subsistence-level groups experiencing rapid lifestyle change. In these populations, comparisons of individuals falling on opposite extremes of the "matched" to "mismatched" spectrum are uniquely possible. More broadly, the work we propose will inform our understanding of environmental and genetic risk factors for NCDs across diverse ancestries and cultures.

Human subsistence and signatures of selection on chemosensory genes.

Chemosensation (olfaction, taste) is essential for detecting and assessing foods, such that dietary shifts elicit evolutionary changes in vertebrate chemosensory genes. The transition from hunting and gathering to agriculture dramatically altered how humans acquire food. Recent genetic and linguistic studies suggest agriculture may have precipitated olfactory degeneration. Here, we explore the effects of subsistence behaviors on olfactory (OR) and taste (TASR) receptor genes among rainforest foragers and neighboring agriculturalists in Africa and Southeast Asia. We analyze 378 functional OR and 26 functional TASR genes in 133 individuals across populations in Uganda (Twa, Sua, BaKiga) and the Philippines (Agta, Mamanwa, Manobo) with differing subsistence histories. We find no evidence of relaxed selection on chemosensory genes in agricultural populations. However, we identify subsistence-related signatures of local adaptation on chemosensory genes within each geographic region. Our results highlight the importance of culture, subsistence economy, and drift in human chemosensory perception.

An integrative skeletal and paleogenomic analysis of stature variation suggests relatively reduced health for early European farmers.

Human culture, biology, and health were shaped dramatically by the onset of agriculture ∼12,000 y B.P. This shift is hypothesized to have resulted in increased individual fitness and population growth as evidenced by archaeological and population genomic data alongside a decline in physiological health as inferred from skeletal remains. Here, we consider osteological and ancient DNA data from the same prehistoric individuals to study human stature variation as a proxy for health across a transition to agriculture. Specifically, we compared "predicted" genetic contributions to height from paleogenomic data and "achieved" adult osteological height estimated from long bone measurements for 167 individuals across Europe spanning the Upper Paleolithic to Iron Age (∼38,000 to 2,400 B.P.). We found that individuals from the Neolithic were shorter than expected (given their individual polygenic height scores) by an average of −3.82 cm relative to individuals from the Upper Paleolithic and Mesolithic (P = 0.040) and −2.21 cm shorter relative to post-Neolithic individuals (P = 0.068), with osteological vs. expected stature steadily increasing across the Copper (+1.95 cm relative to the Neolithic), Bronze (+2.70 cm), and Iron (+3.27 cm) Ages. These results were attenuated when we additionally accounted for genome-wide genetic ancestry variation: for example, with Neolithic individuals −2.82 cm shorter than expected on average relative to pre-Neolithic individuals (P = 0.120). We also incorporated observations of paleopathological indicators of nonspecific stress that can persist from childhood to adulthood in skeletal remains into our model. Overall, our work highlights the potential of integrating disparate datasets to explore proxies of health in prehistory.

A chromosome-level genome assembly for the eastern fence lizard (Sceloporus undulatus), a reptile model for physiological and evolutionary ecology.

BACKGROUND: High-quality genomic resources facilitate investigations into behavioral ecology, morphological and physiological adaptations, and the evolution of genomic architecture. Lizards in the genus Sceloporus have a long history as important ecological, evolutionary, and physiological models, making them a valuable target for the development of genomic resources. FINDINGS: We present a high-quality chromosome-level reference genome assembly, SceUnd1.0 (using 10X Genomics Chromium, HiC, and Pacific Biosciences data), and tissue/developmental stage transcriptomes for the eastern fence lizard, Sceloporus undulatus. We performed synteny analysis with other snake and lizard assemblies to identify broad patterns of chromosome evolution including the fusion of micro- and macrochromosomes. We also used this new assembly to provide improved reference-based genome assemblies for 34 additional Sceloporus species. Finally, we used RNAseq and whole-genome resequencing data to compare 3 assemblies, each representing an increased level of cost and effort: Supernova Assembly with data from 10X Genomics Chromium, HiRise Assembly that added data from HiC, and PBJelly Assembly that added data from Pacific Biosciences sequencing. We found that the Supernova Assembly contained the full genome and was a suitable reference for RNAseq and single-nucleotide polymorphism calling, but the chromosome-level scaffolds provided by the addition of HiC data allowed synteny and whole-genome association mapping analyses. The subsequent addition of PacBio data doubled the contig N50 but provided negligible gains in scaffold length. CONCLUSIONS: These new genomic resources provide valuable tools for advanced molecular analysis of an organism that has become a model in physiology and evolutionary ecology.

Evolutionary medicine.

In recognition that evolutionary theory is critical for understanding modern human health, eLife is publishing a special issue on evolutionary medicine to showcase recent research in this growing and increasingly interdisciplinary field.

Evolutionary and phylogenetic insights from a nuclear genome sequence of the extinct, giant, "subfossil" koala lemur Megaladapis edwardsi.

No endemic Madagascar animal with body mass >10 kg survived a relatively recent wave of extinction on the island. From morphological and isotopic analyses of skeletal "subfossil" remains we can reconstruct some of the biology and behavioral ecology of giant lemurs (primates; up to ∼160 kg) and other extraordinary Malagasy megafauna that survived into the past millennium. Yet, much about the evolutionary biology of these now-extinct species remains unknown, along with persistent phylogenetic uncertainty in some cases. Thankfully, despite the challenges of DNA preservation in tropical and subtropical environments, technical advances have enabled the recovery of ancient DNA from some Malagasy subfossil specimens. Here, we present a nuclear genome sequence (∼2× coverage) for one of the largest extinct lemurs, the koala lemur Megaladapis edwardsi (∼85 kg). To support the testing of key phylogenetic and evolutionary hypotheses, we also generated high-coverage nuclear genomes for two extant lemurs, Eulemur rufifrons and Lepilemur mustelinus, and we aligned these sequences with previously published genomes for three other extant lemurs and 47 nonlemur vertebrates. Our phylogenetic results confirm that Megaladapis is most closely related to the extant Lemuridae (typified in our analysis by E. rufifrons) to the exclusion of L. mustelinus, which contradicts morphology-based phylogenies. Our evolutionary analyses identified significant convergent evolution between M. edwardsi and an extant folivore (a colobine monkey) and an herbivore (horse) in genes encoding proteins that function in plant toxin biodegradation and nutrient absorption. These results suggest that koala lemurs were highly adapted to a leaf-based diet, which may also explain their convergent craniodental morphology with the small-bodied folivore Lepilemur.

Patterns of recent natural selection on genetic loci associated with sexually differentiated human body size and shape phenotypes.

Levels of sex differences for human body size and shape phenotypes are hypothesized to have adaptively reduced following the agricultural transition as part of an evolutionary response to relatively more equal divisions of labor and new technology adoption. In this study, we tested this hypothesis by studying genetic variants associated with five sexually differentiated human phenotypes: height, body mass, hip circumference, body fat percentage, and waist circumference. We first analyzed genome-wide association (GWAS) results for UK Biobank individuals (~194,000 females and ~167,000 males) to identify a total of 114,199 single nucleotide polymorphisms (SNPs) significantly associated with at least one of the studied phenotypes in females, males, or both sexes (P<5x10-8). From these loci we then identified 3,016 SNPs (2.6%) with significant differences in the strength of association between the female- and male-specific GWAS results at a low false-discovery rate (FDR<0.001). Genes with known roles in sexual differentiation are significantly enriched for co-localization with one or more of these SNPs versus SNPs associated with the phenotypes generally but not with sex differences (2.73-fold enrichment; permutation test; P = 0.0041). We also confirmed that the identified variants are disproportionately associated with greater phenotype effect sizes in the sex with the stronger association value. We then used the singleton density score statistic, which quantifies recent (within the last ~3,000 years; post-agriculture adoption in Britain) changes in the frequencies of alleles underlying polygenic traits, to identify a signature of recent positive selection on alleles associated with greater body fat percentage in females (permutation test; P = 0.0038; FDR = 0.0380), directionally opposite to that predicted by the sex differences reduction hypothesis. Otherwise, we found no evidence of positive selection for sex difference-associated alleles for any other trait. Overall, our results challenge the longstanding hypothesis that sex differences adaptively decreased following subsistence transitions from hunting and gathering to agriculture.

Scanning the human genome for "signatures" of positive selection: Transformative opportunities and ethical obligations.

The relationship history of evolutionary anthropology and genetics is complex. At best, genetics is a beautifully integrative part of the discipline. Yet this integration has also been fraught, with punctuated, disruptive challenges to dogma, periodic reluctance by some members of the field to embrace results from analyses of genetic data, and occasional over-assertions of genetic definitiveness by geneticists. At worst, evolutionary genetics has been a tool for reinforcing racism and colonialism. While a number of genetics/genomics papers have disproportionately impacted evolutionary anthropology, here we highlight the 2002 presentation of an elegantly powerful approach for identifying "signatures" of past positive selection from haplotype-based patterns of genetic variation. Together with technological advances in genotyping methods, this article transformed our field by facilitating genome-wide "scans" for signatures of past positive selection in human populations. This approach helped researchers test longstanding evolutionary anthropology hypotheses while simultaneously providing opportunities to develop entirely new ones. Genome-wide scans for signatures of positive selection have since been conducted in diverse worldwide populations, with striking findings of local adaptation and convergent evolution. Yet there are ethical considerations with respect to the ubiquity of these studies and the cross-application of the genome-wide scan approach to existing datasets, which we also discuss.

Modern, archaeological, and paleontological DNA analysis of a human-harvested marine gastropod (Strombus pugilis) from Caribbean Panama.

Although protocols exist for the recovery of ancient DNA from land snail and marine bivalve shells, marine conch shells have yet to be studied from a paleogenomic perspective. We first present reference assemblies for both a 623.7 Mbp nuclear genome and a 15.4 kbp mitochondrial genome for Strombus pugilis, the West Indian fighting conch. We next detail a method to extract and sequence DNA from conch shells and apply it to conch from Bocas del Toro, Panama across three time periods: recently-eaten and discarded (n = 3), Late Holocene (984-1258 before present [BP]) archaeological midden (n = 5), and mid-Holocene (5711-7187 BP) paleontological fossil coral reef (n = 5). These results are compared to control DNA extracted from live-caught tissue and fresh shells (n = 5). Using high-throughput sequencing, we were able to obtain S. pugilis nuclear sequence reads from shells across all age periods: up to 92.5 thousand filtered reads per sample in live-caught shell material, 4.57 thousand for modern discarded shells, 12.1 thousand reads for archaeological shells, and 114 reads in paleontological shells. We confirmed authenticity of the ancient DNA recovered from the archaeological and paleontological shells based on 5.7× higher average frequency of deamination-driven misincorporations and 15% shorter average read lengths compared to the modern shells. Reads also mapped to the S. pugilis mitochondrial genome for all but the paleontological shells, with consistent ratios of mitochondrial to nuclear mapped reads across sample types. Our methods can be applied to diverse archaeological sites to facilitate reconstructions of the long-term impacts of human behaviour on mollusc evolutionary biology.

The genomics of ecological flexibility, large brains, and long lives in capuchin monkeys revealed with fecalFACS.

Ecological flexibility, extended lifespans, and large brains have long intrigued evolutionary biologists, and comparative genomics offers an efficient and effective tool for generating new insights into the evolution of such traits. Studies of capuchin monkeys are particularly well situated to shed light on the selective pressures and genetic underpinnings of local adaptation to diverse habitats, longevity, and brain development. Distributed widely across Central and South America, they are inventive and extractive foragers, known for their sensorimotor intelligence. Capuchins have among the largest relative brain size of any monkey and a lifespan that exceeds 50 y, despite their small (3 to 5 kg) body size. We assemble and annotate a de novo reference genome for Cebus imitator Through high-depth sequencing of DNA derived from blood, various tissues, and feces via fluorescence-activated cell sorting (fecalFACS) to isolate monkey epithelial cells, we compared genomes of capuchin populations from tropical dry forests and lowland rainforests and identified population divergence in genes involved in water balance, kidney function, and metabolism. Through a comparative genomics approach spanning a wide diversity of mammals, we identified genes under positive selection associated with longevity and brain development. Additionally, we provide a technological advancement in the use of noninvasive genomics for studies of free-ranging mammals. Our intra- and interspecific comparative study of capuchin genomics provides insights into processes underlying local adaptation to diverse and physiologically challenging environments, as well as the molecular basis of brain evolution and longevity.

Factors influencing taxonomic unevenness in scientific research: a mixed-methods case study of non-human primate genomic sequence data generation.

Scholars have noted major disparities in the extent of scientific research conducted among taxonomic groups. Such trends may cascade if future scientists gravitate towards study species with more data and resources already available. As new technologies emerge, do research studies employing these technologies continue these disparities? Here, using non-human primates as a case study, we identified disparities in massively parallel genomic sequencing data and conducted interviews with scientists who produced these data to learn their motivations when selecting study species. We tested whether variables including publication history and conservation status were significantly correlated with publicly available sequence data in the NCBI Sequence Read Archive (SRA). Of the 179.6 terabases (Tb) of sequence data in SRA for 519 non-human primate species, 135 Tb (approx. 75%) were from only five species: rhesus macaques, olive baboons, green monkeys, chimpanzees and crab-eating macaques. The strongest predictors of the amount of genomic data were the total number of non-medical publications (linear regression; r 2 = 0.37; p = 6.15 × 10-12) and number of medical publications (r 2 = 0.27; p = 9.27 × 10-9). In a generalized linear model, the number of non-medical publications (p = 0.00064) and closer phylogenetic distance to humans (p = 0.024) were the most predictive of the amount of genomic sequence data. We interviewed 33 authors of genomic data-producing publications and analysed their responses using grounded theory. Consistent with our quantitative results, authors mentioned their choice of species was motivated by sample accessibility, prior published work and relevance to human medicine. Our mixed-methods approach helped identify and contextualize some of the driving factors behind species-uneven patterns of scientific research, which can now be considered by funding agencies, scientific societies and research teams aiming to align their broader goals with future data generation efforts.

How human behavior can impact the evolution of genetically-mediated behavior in wild non-human species.

Humans intensely modify the ecosystems we inhabit. Many of the impacts that this behavior can have on other species also sharing these spaces are obvious. A prime example is the devastating current extinction crisis. Yet some populations of non-human, non-domesticated species survive or even appear to thrive in heavily disturbed or human-built habitats. Theoretically, this apparent paradox could be facilitated partly by the evolution of genetically-mediated trait adaptations to the impacts of human behavior. At the least, persistence in strongly modified habitats would provide requisite selection pressures for this process to potentially occur in the future. In fact, we have a growing number of well-characterized examples of morphological trait adaptations to human behavior. However, our knowledge of genetically-mediated behavioral adaptations in similar contexts is less well developed. In this review I set up and discuss several evolutionary scenarios by which human behavior might have impacted the evolution of genetically mediated behavior in non-human, non-domestic species and highlight several approaches that could be used in future studies of this process.

Studying human and nonhuman primate evolutionary biology with powerful in vitro and in vivo functional genomics tools.

In recent years, tools for functional genomic studies have become increasingly feasible for use by evolutionary anthropologists. In this review, we provide brief overviews of several exciting in vitro techniques that can be paired with "-omics" approaches (e.g., genomics, epigenomics, transcriptomics, proteomics, and metabolomics) for potentially powerful evolutionary insights. These in vitro techniques include ancestral protein resurrection, cell line experiments using primary, immortalized, and induced pluripotent stem cells, and CRISPR-Cas9 genetic manipulation. We also discuss how several of these methods can be used in vivo, for transgenic organism studies of human and nonhuman primate evolution. Throughout this review, we highlight example studies in which these approaches have already been used to inform our understanding of the evolutionary biology of modern and archaic humans and other primates while simultaneously identifying future opportunities for anthropologists to use this toolkit to help answer additional outstanding questions in evolutionary anthropology.

Genomic Evidence for Local Adaptation of Hunter-Gatherers to the African Rainforest.

African rainforests support exceptionally high biodiversity and host the world's largest number of active hunter-gatherers [1-3]. The genetic history of African rainforest hunter-gatherers and neighboring farmers is characterized by an ancient divergence more than 100,000 years ago, together with recent population collapses and expansions, respectively [4-12]. While the demographic past of rainforest hunter-gatherers has been deeply characterized, important aspects of their history of genetic adaptation remain unclear. Here, we investigated how these groups have adapted-through classic selective sweeps, polygenic adaptation, and selection since admixture-to the challenging rainforest environments. To do so, we analyzed a combined dataset of 566 high-coverage exomes, including 266 newly generated exomes, from 14 populations of rainforest hunter-gatherers and farmers, together with 40 newly generated, low-coverage genomes. We find evidence for a strong, shared selective sweep among all hunter-gatherer groups in the regulatory region of TRPS1-primarily involved in morphological traits. We detect strong signals of polygenic adaptation for height and life history traits such as reproductive age; however, the latter appear to result from pervasive pleiotropy of height-associated genes. Furthermore, polygenic adaptation signals for functions related to responses of mast cells to allergens and microbes, the IL-2 signaling pathway, and host interactions with viruses support a history of pathogen-driven selection in the rainforest. Finally, we find that genes involved in heart and bone development and immune responses are enriched in both selection signals and local hunter-gatherer ancestry in admixed populations, suggesting that selection has maintained adaptive variation in the face of recent gene flow from farmers.

Natural selection contributed to immunological differences between hunter-gatherers and agriculturalists.

The shift from a hunter-gatherer to an agricultural mode of subsistence is believed to have been associated with profound changes in the burden and diversity of pathogens across human populations. Yet, the extent to which the advent of agriculture affected the evolution of the human immune system remains unknown. Here we present a comparative study of variation in the transcriptional responses of peripheral blood mononuclear cells to bacterial and viral stimuli between Batwa rainforest hunter-gatherers and Bakiga agriculturalists from Uganda. We observed increased divergence between hunter-gatherers and agriculturalists in the early transcriptional response to viruses compared with that for bacterial stimuli. We demonstrate that a significant fraction of these transcriptional differences are under genetic control and we show that positive natural selection has helped to shape population differences in immune regulation. Across the set of genetic variants underlying inter-population immune-response differences, however, the signatures of positive selection were disproportionately observed in the rainforest hunter-gatherers. This result is counter to expectations on the basis of the popularized notion that shifts in pathogen exposure due to the advent of agriculture imposed radically heightened selective pressures in agriculturalist populations.

Impact of DNA degradation on massively parallel sequencing-based autosomal STR, iiSNP, and mitochondrial DNA typing systems.

Biological samples, including skeletal remains exposed to environmental insults for extended periods of time, exhibit increasing levels of DNA damage and fragmentation. Human forensic identification methods typically use a combination of mitochondrial (mt) DNA sequencing and short tandem repeat (STR) analysis, which target segments of DNA ranging from 80 to 500 base pairs (bps). Larger templates are often unavailable as skeletal samples age and the associated DNA degrades. Single-nucleotide polymorphism (SNP) loci target shorter templates and may serve as a solution to the problem. Recently developed assays for STR and SNP analysis using a massively parallel sequencing approach, such as the ForenSeq kit (Verogen, San Diego, CA), offer a means for generating results from degraded samples as they target templates down to 60 to 170 bps. We performed a modeling study that demonstrates that SNPs can increase the significance of an identification when analyzing DNA down to an average size of 100 bps for input amounts between 0.375 and 1 ng of nuclear DNA. Observations from this study were then compared with human skeletal material results (n = 14, ninth to eighteenth centuries), which further demonstrated the utility of the ForenSeq kit for degraded samples. The robustness of the Promega PowerSeq™ Mito System was also tested with human skeletal remains (n = 70, ninth to eighteenth centuries), resulting in successful coverage of 99.29% of the mtDNA control region at 50× coverage or more. This was accompanied by modifications to a mainstream DNA extraction technique for skeletal remains that improved recovery of shorter templates.

Horse Paleogenomes and Human-Animal Interactions in Prehistory.

A new analysis of paleogenomic data from 278 ancient horses (Fages et al. Cellhttp://doi.org/10.1016/j.cell.2019.03.049) finds that this animal - crucially important to many ancient and contemporary human societies for subsistence, transportation, conflict, and more - was domesticated in at least two different regions, but with the geographic and cultural origins of the modern domestic horse lineage remaining unknown. By tracing ancient horse population movements and inferring the spatiotemporal trajectories of phenotypic adaptations, this study provides fresh perspectives on past human group interactions and activities.

An environmental DNA sampling method for aye-ayes from their feeding traces.

Noninvasive sampling is an important development in population genetic monitoring of wild animals. Particularly, the collection of environmental DNA (eDNA) which can be collected without needing to encounter the target animal facilitates the genetic analysis of endangered species. One method that has been applied to these sample types is target capture and enrichment which overcomes the issue of high proportions of exogenous (nonhost) DNA from these lower quality samples. We tested whether target capture of mitochondrial DNA from sampled feeding traces of the aye-aye, an endangered lemur species would yield mitochondrial DNA sequences for population genetic monitoring. We sampled gnawed wood where aye-ayes excavate wood-boring insect larvae from trees. We designed RNA probes complementary to the aye-aye's mitochondrial genome and used these to isolate aye-aye DNA from other nontarget DNA in these samples. We successfully retrieved six near-complete mitochondrial genomes from two sites within the aye-aye's geographic range that had not been sampled previously. Our method demonstrates the application of next-generation molecular techniques to species of conservation concern. This method can likely be applied to alternative foraged remains to sample endangered species other than aye-ayes.

Polygenic adaptation and convergent evolution on growth and cardiac genetic pathways in African and Asian rainforest hunter-gatherers.

Different human populations facing similar environmental challenges have sometimes evolved convergent biological adaptations, for example, hypoxia resistance at high altitudes and depigmented skin in northern latitudes on separate continents. The "pygmy" phenotype (small adult body size), characteristic of hunter-gatherer populations inhabiting both African and Asian tropical rainforests, is often highlighted as another case of convergent adaptation in humans. However, the degree to which phenotypic convergence in this polygenic trait is due to convergent versus population-specific genetic changes is unknown. To address this question, we analyzed high-coverage sequence data from the protein-coding portion of the genomes of two pairs of populations: Batwa rainforest hunter-gatherers and neighboring Bakiga agriculturalists from Uganda and Andamanese rainforest hunter-gatherers and Brahmin agriculturalists from India. We observed signatures of convergent positive selection between the rainforest hunter-gatherers across the set of genes with "growth factor binding" functions ([Formula: see text]). Unexpectedly, for the rainforest groups, we also observed convergent and population-specific signatures of positive selection in pathways related to cardiac development (e.g., "cardiac muscle tissue development"; [Formula: see text]). We hypothesize that the growth hormone subresponsiveness likely underlying the adult small body-size phenotype may have led to compensatory changes in cardiac pathways, in which this hormone also plays an essential role. Importantly, in the agriculturalist populations, we did not observe similar patterns of positive selection on sets of genes associated with growth or cardiac development, indicating our results most likely reflect a history of convergent adaptation to the similar ecology of rainforests rather than a more general evolutionary pattern.

The demographic history and mutational load of African hunter-gatherers and farmers.

Understanding how deleterious genetic variation is distributed across human populations is of key importance in evolutionary biology and medical genetics. However, the impact of population size changes and gene flow on the corresponding mutational load remains a controversial topic. Here, we report high-coverage exomes from 300 rainforest hunter-gatherers and farmers of central Africa, whose distinct subsistence strategies are expected to have impacted their demographic pasts. Detailed demographic inference indicates that hunter-gatherers and farmers recently experienced population collapses and expansions, respectively, accompanied by increased gene flow. We show that the distribution of deleterious alleles across these populations is compatible with a similar efficacy of selection to remove deleterious variants with additive effects, and predict with simulations that their present-day additive mutation load is almost identical. For recessive mutations, although an increased load is predicted for hunter-gatherers, this increase has probably been partially counteracted by strong gene flow from expanding farmers. Collectively, our predicted and empirical observations suggest that the impact of the recent population decline of African hunter-gatherers on their mutation load has been modest and more restrained than would be expected under a fully recessive model of dominance.