Population Genomics of Human Adaptation

A novel expectation-maximization approach to infer general diploid selection from time-series genetic data

Detecting and quantifying the strength of selection is a main objective in population genetics. Since selection acts over multiple generations, many approaches have been developed to detect and quantify selection using genetic data sampled at multiple points in time. Such time series genetic data is commonly analyzed using Hidden Markov Models, but in most cases, under the assumption of additive selection. However, many examples of genetic variation exhibiting non-additive mechanisms exist, making it critical to develop methods that can characterize selection in more general scenarios. Thus, we extend a previously introduced expectation-maximization algorithm for the inference of additive selection coefficients to the case of general diploid selection, in which heterozygote and homozygote fitnesses are parameterized independently. We furthermore introduce a framework to identify bespoke modes of diploid selection from given data, as well as a procedure for aggregating data across linked loci to increase power and robustness. Using extensive simulation studies, we find that our method accurately and efficiently estimates selection coefficients for different modes of diploid selection across a wide range of scenarios; however, power to classify the mode of selection is low unless selection is very strong. We apply our method to ancient DNA samples from Great Britain in the last 4,450 years, and detect evidence for selection in six genomic regions, including the well-characterized LCT locus. Our work is the first genome-wide scan characterizing signals of general diploid selection.

Pearl Memorial Lecture. Humans at the extremes: Exploring human adaptation to ecological and social stressors

The field of human biology has long explored how human populations have adapted to extreme environmental circumstances. Yet, it has become increasingly clear that conditions of social stress, poverty, and lifestyle change play equally important roles in shaping human biological variation and health. In this paper, I provide a brief background on the foundational human adaptability research of the International Biological Programme (IBP) from the 1960s, highlighting how its successes and critiques have shaped current research directions in the field. I then discuss and reflect on my own field research that has examined the influence of both environmental and social stresses on human populations living in different ecosystems: the Peruvian Andes, the Siberian arctic, and the Bolivian rainforest. Finally, I consider how the papers in this special issue advance our understanding of human adaptability to extreme conditions and offer directions for future research. Drawing on our field's distinctive evolutionary and biocultural perspectives, human biologists are uniquely positioned to examine how the interplay between social and ecological domains influences the human condition.

Evolutionary Insights from a Large-Scale Survey of Population-Genomic Variation

The field of genomics has ushered in new methods for studying molecular-genetic variation in natural populations. However, most population-genomic studies still rely on small sample sizes (typically, <100 individuals) from single time points, leaving considerable uncertainties with respect to the behavior of relatively young (and rare) alleles and, owing to the large sampling variance of measures of variation, to the specific gene targets of unusually strong selection. Genomic sequences of ∼1,700 haplotypes distributed over a 10-year period from a natural population of the microcrustacean Daphnia pulex reveal evolutionary-genomic features at a refined scale, including previously hidden information on the behavior of rare alleles predicted by recent theory. Background selection, resulting from the recurrent introduction of deleterious alleles, appears to strongly influence the dynamics of neutral alleles, inducing indirect negative selection on rare variants and positive selection on common variants. Temporally fluctuating selection increases the persistence of nonsynonymous alleles with intermediate frequencies, while reducing standing levels of variation at linked silent sites. Combined with the results from an equally large metapopulation survey of the study species, classes of genes that are under strong positive selection can now be confidently identified in this key model organism. Most notable among rapidly evolving Daphnia genes are those associated with ribosomes, mitochondrial functions, sensory systems, and lifespan determination.

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.

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.

Evolutionary Insights from a Large-scale Survey of Population-genomic Variation

Results from data on > 1000 haplotypes distributed over a nine-year period from a natural population of the microcrustacean Daphnia pulex reveal evolutionary-genomic features at a refined scale, including key population-genetic properties that are obscured in studies with smaller sample sizes. Background selection, resulting from the recurrent introduction of deleterious alleles, appears to strongly influence the dynamics of neutral alleles, inducing indirect negative selection on rare variants and positive selection on common variants. Fluctuating selection increases the persistence of nonsynonymous alleles with intermediate frequencies, while reducing standing levels of variation at linked silent sites. Combined with the results from an equally large metapopulation survey of the study species, regions of gene structure that are under strong purifying selection and classes of genes that are under strong positive selection in this key species can be confidently identified. Most notable among rapidly evolving Daphnia genes are those associated with ribosomes, mitochondrial functions, sensory systems, and lifespan determination.

Balancing selection on genomic deletion polymorphisms in humans

A key question in biology is why genomic variation persists in a population for extended periods. Recent studies have identified examples of genomic deletions that have remained polymorphic in the human lineage for hundreds of millennia, ostensibly owing to balancing selection. Nevertheless, genome-wide investigation of ancient and possibly adaptive deletions remains an imperative exercise. Here, we demonstrate an excess of polymorphisms in present-day humans that predate the modern human-Neanderthal split (ancient polymorphisms), which cannot be explained solely by selectively neutral scenarios. We analyze the adaptive mechanisms that underlie this excess in deletion polymorphisms. Using a previously published measure of balancing selection, we show that this excess of ancient deletions is largely owing to balancing selection. Based on the absence of signatures of overdominance, we conclude that it is a rare mode of balancing selection among ancient deletions. Instead, more complex scenarios involving spatially and temporally variable selective pressures are likely more common mechanisms. Our results suggest that balancing selection resulted in ancient deletions harboring disproportionately more exonic variants with GWAS (genome-wide association studies) associations. We further found that ancient deletions are significantly enriched for traits related to metabolism and immunity. As a by-product of our analysis, we show that deletions are, on average, more deleterious than single nucleotide variants. We can now argue that not only is a vast majority of common variants shared among human populations, but a considerable portion of biologically relevant variants has been segregating among our ancestors for hundreds of thousands, if not millions, of years.

Human genetic admixture through the lens of population genomics

Over the past 50 years, geneticists have made great strides in understanding how our species' evolutionary history gave rise to current patterns of human genetic diversity classically summarized by Lewontin in his 1972 paper, 'The Apportionment of Human Diversity'. One evolutionary process that requires special attention in both population genetics and statistical genetics is admixture: gene flow between two or more previously separated source populations to form a new admixed population. The admixture process introduces ancestry-based structure into patterns of genetic variation within and between populations, which in turn influences the inference of demographic histories, identification of genetic targets of selection and prediction of complex traits. In this review, we outline some challenges for admixture population genetics, including limitations of applying methods designed for populations without recent admixture to the study of admixed populations. We highlight recent studies and methodological advances that aim to overcome such challenges, leveraging genomic signatures of admixture that occurred in the past tens of generations to gain insights into human history, natural selection and complex trait architecture. This article is part of the theme issue 'Celebrating 50 years since Lewontin's apportionment of human diversity'.

Adaptive eQTLs reveal the evolutionary impacts of pleiotropy and tissue-specificity while contributing to health and disease

Large numbers of expression quantitative trait loci (eQTLs) have recently been identified in humans, and many of these regulatory variants have large allele frequency differences between populations. Here, we conducted genome-wide scans of selection to identify adaptive eQTLs (i.e., eQTLs with large population branch statistics). We then tested if tissue pleiotropy affects whether eQTLs are more or less likely to be adaptive and identified tissues that have been key targets of positive selection during the last 100,000 years. Top adaptive eQTL outliers include rs1043809, rs66899053, and rs2814778 (a SNP that is associated with malaria resistance). We found that effect sizes of eQTLs were negatively correlated with population branch statistics and that adaptive eQTLs affect two-thirds as many tissues as do non-adaptive eQTLs. Because the tissue breadth of an eQTL can be viewed as a measure of pleiotropy, these results imply that pleiotropy inhibits adaptation. The proportion of eQTLs that are adaptive varies by tissue, and we found that eQTLs that regulate expression in testis, thyroid, blood, or sun-exposed skin are enriched for signatures of positive selection. By contrast, eQTLs that regulate expression in the cerebrum or female-specific tissues have a relative lack of adaptive outliers. Scans of selections also reveal that many adaptive eQTLs are closely linked to disease-associated loci. Taken together, our results indicate that eQTLs have played an important role in recent human evolution.

Controlled Ovarian Hyperstimulation Protocol in Infertile Patients During the COVID-19 Pandemic

Objectives: To explore the appropriate controlled ovarian hyperstimulation (COH) protocols in infertility patients who received the in vitro fertilization (IVF)/intracytoplasmic sperm injection (ICSI) treatments during the COVID-19 pandemic.

Materials and Methods: This retrospective cohort study evaluated the efficiency of the early follicular-phase long-acting GnRH-agonist long (EFLL) protocol (a new protocol developed by Chinese clinicians), prolonged pituitary down-regulation of EFLL protocol (Pro-EFLL), and the GnRH-ant protocol for couples meeting the study criteria between February 2020 and June 2020 who were treated by the First Affiliated Hospital of Zhengzhou University during the COVID-19 pandemic, and compared the pregnancy rates and miscarriage rates per fresh transfer cycle, number of retrieved oocytes, endometrial thickness on the day of hCG injection and the number of fertilized oocytes, mature oocytes, fertilized oocytes, and transferable embryos among the three protocols.

Results: We found that the prolonged pituitary down-regulation during the COVID-19 pandemic by utilizing a full-dose of GnRH-a administrated in infertility patients were no differences in clinical outcomes than other protocols, The prolonged pituitary down-regulation protocol and EFLL protocol were associated with a higher Endometrial thickness on the day of hCG injection (12.67 ± 2.21 vs. 12.09 ± 2.35 vs. 10.79 ± 2.38, P < 0.001), retrieved oocytes (14.49 ± 6.30 vs. 15.02 ± 7.93 vs. 10.06 ± 7.63, P < 0.001), mature oocytes (11.60 ± 5.71 vs. 11.96 ± 6.00 vs. 7.63 ± 6.50, P < 0.001), fertilized oocytes (9.14 ± 5.43 vs. 8.44 ± 5.34 vs. 5.42 ± 5.20, P < 0.001), and transferable embryos (4.87 ± 2.96 vs. 6.47 ± 5.12 vs. 3.00 ± 3.28 vs. P < 0.001) in the GnRH-antagonist protocol.

Conclusion: We recommend that patients start Gn injections 33–42 days after a pituitary downregulated full dose (3.75 mg) of gonadotropin-releasing hormone agonist during the COVID-19 pandemic, even a delay of 2–4 weeks does not affect the implantation rate. The study can provide a more detailed estimate and clinical management strategies for infertile couples during the COVID-19 pandemic.

Detecting Genetic Ancestry and Adaptation in the Taiwanese Han People

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

Human biological variation and the "normal"

Anatomically modern human being is a relatively young species (~300 000 years old) with small amounts of genetic variation contained within them. The vast majority of its existence was spent in Eastern Africa, migration out of the region began around 100 000 YBP. Sub-Saharan African populations have the greatest amount of human genetic variation. However, migration allowed populations to accumulate genomic variation associated with living in the arctic, higher altitudes, disease resistance, living on high fat or starchy foods, surviving toxic arsenic-rich environments, lactase persistence, changing skin pigmentation, gaining thicker hair, and changing height and body mass index. Understanding these aspects of human evolution forces us to reconsider our notion of the "normal." Thus, normal for our species includes having dark melanic skin, brown eyes, and brown tightly curled hair. Derived features include lighter skin (~10 000 YBP), blue eyes (~6000 YBP), and blond straight hair (~6000 YBP). Yet in reality, "normal" has no meaning for a species that inhabits such a broad geographic range. Natural selection and genetic drift have genetically differentiated human populations in ways that impact our morphological and physiological traits. The genomic differentiation is small and does not allow any unambiguous classification of human populations into biological races. Despite these now well-established facts of human variation, significant confusion associated with Eurocentric notions of the normal still persist in both the lay public and various professions such as biomedical research and clinical practice.

Positive Selection in Gene Regulatory Factors Suggests Adaptive Pleiotropic Changes During Human Evolution

Gene regulatory factors (GRFs), such as transcription factors, co-factors and histone-modifying enzymes, play many important roles in modifying gene expression in biological processes. They have also been proposed to underlie speciation and adaptation. To investigate potential contributions of GRFs to primate evolution, we analyzed GRF genes in 27 publicly available primate genomes. Genes coding for zinc finger (ZNF) proteins, especially ZNFs with a Krüppel-associated box (KRAB) domain were the most abundant TFs in all genomes. Gene numbers per TF family differed between all species. To detect signs of positive selection in GRF genes we investigated more than 3,000 human GRFs with their more than 70,000 orthologs in 26 non-human primates. We implemented two independent tests for positive selection, the branch-site-model of the PAML suite and aBSREL of the HyPhy suite, focusing on the human and great ape branch. Our workflow included rigorous procedures to reduce the number of false positives: excluding distantly similar orthologs, manual corrections of alignments, and considering only genes and sites detected by both tests for positive selection. Furthermore, we verified the candidate sites for selection by investigating their variation within human and non-human great ape population data. In order to approximately assign a date to positively selected sites in the human lineage, we analyzed archaic human genomes. Our work revealed with high confidence five GRFs that have been positively selected on the human lineage and one GRF that has been positively selected on the great ape lineage. These GRFs are scattered on different chromosomes and have been previously linked to diverse functions. For some of them a role in speciation and/or adaptation can be proposed based on the expression pattern or association with human diseases, but it seems that they all contributed independently to human evolution. Four of the positively selected GRFs are KRAB-ZNF proteins, that induce changes in target genes co-expression and/or through arms race with transposable elements. Since each positively selected GRF contains several sites with evidence for positive selection, we suggest that these GRFs participated pleiotropically to phenotypic adaptations in humans.

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.

Rapid adaptation to malaria facilitated by admixture in the human population of Cabo Verde

Humans have undergone large migrations over the past hundreds to thousands of years, exposing ourselves to new environments and selective pressures. Yet, evidence of ongoing or recent selection in humans is difficult to detect. Many of these migrations also resulted in gene flow between previously separated populations. These recently admixed populations provide unique opportunities to study rapid evolution in humans. Developing methods based on distributions of local ancestry, we demonstrate that this sort of genetic exchange has facilitated detectable adaptation to a malaria parasite in the admixed population of Cabo Verde within the last ~20 generations. We estimate that the selection coefficient is approximately 0.08, one of the highest inferred in humans. Notably, we show that this strong selection at a single locus has likely affected patterns of ancestry genome-wide, potentially biasing demographic inference. Our study provides evidence of adaptation in a human population on historical timescales.

Insertion variants missing in the human reference genome are widespread among human populations

Background

Structural variants comprise diverse genomic arrangements including deletions, insertions, inversions, and translocations, which can generally be detected in humans through sequence comparison to the reference genome. Among structural variants, insertions are the least frequently identified variants, mainly due to ascertainment bias in the reference genome, lack of previous sequence knowledge, and low complexity of typical insertion sequences. Though recent developments in long-read sequencing deliver promise in annotating individual non-reference insertions, population-level catalogues on non-reference insertion variants have not been identified and the possible functional roles of these hidden variants remain elusive.

Results

To detect non-reference insertion variants, we developed a pipeline, InserTag, which generates non-reference contigs by local de novo assembly and then infers the full-sequence of insertion variants by tracing contigs from non-human primates and other human genome assemblies. Application of the pipeline to data from 2535 individuals of the 1000 Genomes Project helped identify 1696 non-reference insertion variants and re-classify the variants as retention of ancestral sequences or novel sequence insertions based on the ancestral state. Genotyping of the variants showed that individuals had, on average, 0.92-Mbp sequences missing from the reference genome, 92% of the variants were common (allele frequency > 5%) among human populations, and more than half of the variants were major alleles. Among human populations, African populations were the most divergent and had the most non-reference sequences, which was attributed to the greater prevalence of high-frequency insertion variants. The subsets of insertion variants were in high linkage disequilibrium with phenotype-associated SNPs and showed signals of recent continent-specific selection.

Conclusions

Non-reference insertion variants represent an important type of genetic variation in the human population, and our developed pipeline, InserTag, provides the frameworks for the detection and genotyping of non-reference sequences missing from human populations.

Supplementary information

Supplementary information accompanies this paper at 10.1186/s12915-020-00894-1.

Elevation as a selective force on mitochondrial respiratory chain complexes of the Phrynocephalus lizards in the Tibetan plateau

Animals living in extremely high elevations have to adapt to low temperatures and low oxygen availability (hypoxia), but the underlying genetic mechanisms associated with these adaptations are still unclear. The mitochondrial respiratory chain can provide >95% of the ATP in animal cells, and its efficiency is influenced by temperature and oxygen availability. Therefore, the respiratory chain complexes (RCCs) could be important molecular targets for positive selection associated with respiratory adaptation in high-altitude environments. Here, we investigated positive selection in 5 RCCs and their assembly factors by analyzing sequences of 106 genes obtained through RNA-seq of all 15 Chinese Phrynocephalus lizard species, which are distributed from lowlands to the Tibetan plateau (average elevation >4,500 m). Our results indicate that evidence of positive selection on RCC genes is not significantly different from assembly factors, and we found no difference in selective pressures among the 5 complexes. We specifically looked for positive selection in lineages where changes in habitat elevation happened. The group of lineages evolving from low to high altitude show stronger signals of positive selection than lineages evolving from high to low elevations. Lineages evolving from low to high elevation also have more shared codons under positive selection, though the changes are not equivalent at the amino acid level. This study advances our understanding of the genetic basis of animal respiratory metabolism evolution in extreme high environments and provides candidate genes for further confirmation with functional analyses.

In the presence of population structure: From genomics to candidate genes underlying local adaptation

Understanding the genomic signatures, genes, and traits underlying local adaptation of organisms to heterogeneous environments is of central importance to the field evolutionary biology. To identify loci underlying local adaptation, models that combine allelic and environmental variation while controlling for the effects of population structure have emerged as the method of choice. Despite being evaluated in simulation studies, there has not been a thorough investigation of empirical evidence supporting local adaptation across these alleles. To evaluate these methods, we use 875 Arabidopsis thaliana Eurasian accessions and two mixed models (GEMMA and LFMM) to identify candidate SNPs underlying local adaptation to climate. Subsequently, to assess evidence of local adaptation and function among significant SNPs, we examine allele frequency differentiation and recent selection across Eurasian populations, in addition to their distribution along quantitative trait loci (QTL) explaining fitness variation between Italy and Sweden populations and cis-regulatory/nonsynonymous sites showing significant selective constraint. Our results indicate that significant LFMM/GEMMA SNPs show low allele frequency differentiation and linkage disequilibrium across locally adapted Italy and Sweden populations, in addition to a poor association with fitness QTL peaks (highest logarithm of odds score). Furthermore, when examining derived allele frequencies across the Eurasian range, we find that these SNPs are enriched in low-frequency variants that show very large climatic differentiation but low levels of linkage disequilibrium. These results suggest that their enrichment along putative functional sites most likely represents deleterious variation that is independent of local adaptation. Among all the genomic signatures examined, only SNPs showing high absolute allele frequency differentiation (AFD) and linkage disequilibrium (LD) between Italy and Sweden populations showed a strong association with fitness QTL peaks and were enriched along selectively constrained cis-regulatory/nonsynonymous sites. Using these SNPs, we find strong evidence linking flowering time, freezing tolerance, and the abscisic-acid pathway to local adaptation.

In the presence of population structure: From genomics to candidate genes underlying local adaptation

Understanding the genomic signatures, genes, and traits underlying local adaptation of organisms to heterogeneous environments is of central importance to the field evolutionary biology. To identify loci underlying local adaptation, models that combine allelic and environmental variation while controlling for the effects of population structure have emerged as the method of choice. Despite being evaluated in simulation studies, there has not been a thorough investigation of empirical evidence supporting local adaptation across these alleles. To evaluate these methods, we use 875 Arabidopsis thaliana Eurasian accessions and two mixed models (GEMMA and LFMM) to identify candidate SNPs underlying local adaptation to climate. Subsequently, to assess evidence of local adaptation and function among significant SNPs, we examine allele frequency differentiation and recent selection across Eurasian populations, in addition to their distribution along quantitative trait loci (QTL) explaining fitness variation between Italy and Sweden populations and cis-regulatory/nonsynonymous sites showing significant selective constraint. Our results indicate that significant LFMM/GEMMA SNPs show low allele frequency differentiation and linkage disequilibrium across locally adapted Italy and Sweden populations, in addition to a poor association with fitness QTL peaks (highest logarithm of odds score). Furthermore, when examining derived allele frequencies across the Eurasian range, we find that these SNPs are enriched in low-frequency variants that show very large climatic differentiation but low levels of linkage disequilibrium. These results suggest that their enrichment along putative functional sites most likely represents deleterious variation that is independent of local adaptation. Among all the genomic signatures examined, only SNPs showing high absolute allele frequency differentiation (AFD) and linkage disequilibrium (LD) between Italy and Sweden populations showed a strong association with fitness QTL peaks and were enriched along selectively constrained cis-regulatory/nonsynonymous sites. Using these SNPs, we find strong evidence linking flowering time, freezing tolerance, and the abscisic-acid pathway to local adaptation.

A genomics approach reveals the global genetic polymorphism, structure, and functional diversity of ten accessions of the marine model diatom Phaeodactylum tricornutum

Diatoms emerged in the Mesozoic period and presently constitute one of the main primary producers in the world's ocean and are of a major economic importance. In the current study, using whole genome sequencing of ten accessions of the model diatom Phaeodactylum tricornutum, sampled at broad geospatial and temporal scales, we draw a comprehensive landscape of the genomic diversity within the species. We describe strong genetic subdivisions of the accessions into four genetic clades (A-D) with constituent populations of each clade possessing a conserved genetic and functional makeup, likely a consequence of the limited dispersal of P. tricornutum in the open ocean. We further suggest dominance of asexual reproduction across all the populations, as implied by high linkage disequilibrium. Finally, we show limited yet compelling signatures of genetic and functional convergence inducing changes in the selection pressure on many genes and metabolic pathways. We propose these findings to have significant implications for understanding the genetic structure of diatom populations in nature and provide a framework to assess the genomic underpinnings of their ecological success and impact on aquatic ecosystems where they play a major role. Our work provides valuable resources for functional genomics and for exploiting the biotechnological potential of this model diatom species.

Nutrition, the visceral immune system, and the evolutionary origins of pathogenic obesity

The global obesity epidemic is the subject of an immense, diversely specialized research effort. An evolutionary analysis reveals connections among disparate findings, starting with two well-documented facts: Obesity-associated illnesses (e.g., type-2 diabetes and cardiovascular disease), are especially common in: (i) adults with abdominal obesity, especially enlargement of visceral adipose tissue (VAT), a tissue with important immune functions; and (ii) individuals with poor fetal nutrition whose nutritional input increases later in life. I hypothesize that selection favored the evolution of increased lifelong investment in VAT in individuals likely to suffer lifelong malnutrition because of its importance in fighting intraabdominal infections. Then, when increased nutrition violates the adaptive fetal prediction of lifelong nutritional deficit, preferential VAT investment could contribute to abdominal obesity and chronic inflammatory disease. VAT prioritization may help explain several patterns of nutrition-related disease: the paradoxical increase of chronic disease with increased food availability in recently urbanized and migrant populations; correlations between poor fetal nutrition, improved childhood (catch-up) growth, and adult metabolic syndrome; and survival differences between children with marasmus and kwashiorkor malnutrition. Fats and sugars can aggravate chronic inflammation via effects on intestinal bacteria regulating gut permeability to visceral pathogens. The extremes in a nutrition-sensitive trade-off between visceral (immune-function) vs. subcutaneous (body shape) adiposity may have been favored by selection in highly stratified premedicine societies. Altered adipose allocation in populations with long histories of social stratification and malnutrition may be the result of genetic accommodation of developmental responses to poor maternal/fetal conditions, increasing their vulnerability to inflammatory disease.

Interpreting polygenic scores, polygenic adaptation, and human phenotypic differences

Recent analyses of polygenic scores have opened new discussions concerning the genetic basis and evolutionary significance of differences among populations in distributions of phenotypes. Here, we highlight limitations in research on polygenic scores, polygenic adaptation and population differences. We show how genetic contributions to traits, as estimated by polygenic scores, combine with environmental contributions so that differences among populations in trait distributions need not reflect corresponding differences in genetic propensity. Under a null model in which phenotypes are selectively neutral, genetic propensity differences contributing to phenotypic differences among populations are predicted to be small. We illustrate this null hypothesis in relation to health disparities between African Americans and European Americans, discussing alternative hypotheses with selective and environmental effects. Close attention to the limitations of research on polygenic phenomena is important for the interpretation of their relationship to human population differences.

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

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

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

Searching for ancient balanced polymorphisms shared between Neanderthals and Modern Humans

Hominin evolution is characterized by adaptive solutions often rooted in behavioral and cognitive changes. If balancing selection had an important and long-lasting impact on the evolution of these traits, it can be hypothesized that genes associated with them should carry an excess of shared polymorphisms (trans- SNPs) across recent Homo species. In this study, we investigate the role of balancing selection in human evolution using available exomes from modern (Homo sapiens) and archaic humans (H. neanderthalensis and Denisovan) for an excess of trans-SNP in two gene sets: one associated with the immune system (IMMS) and another one with behavioral system (BEHS). We identified a significant excess of trans-SNPs in IMMS (N=547), of which six of these located within genes previously associated with schizophrenia. No excess of trans-SNPs was found in BEHS, but five genes in this system harbor potential signals for balancing selection and are associated with psychiatric or neurodevelopmental disorders. Our approach evidenced recent Homo trans-SNPs that have been previously implicated in psychiatric diseases such as schizophrenia, suggesting that a genetic repertoire common to the immune and behavioral systems could have been maintained by balancing selection starting before the split between archaic and modern humans.

Adaptation of Arabidopsis thaliana to the Yangtze River basin

Background

Organisms need to adapt to keep pace with a changing environment. Examining recent range expansion aids our understanding of how organisms evolve to overcome environmental constraints. However, how organisms adapt to climate changes is a crucial biological question that is still largely unanswered. The plant Arabidopsis thaliana is an excellent system to study this fundamental question. Its origin is in the Iberian Peninsula and North Africa, but it has spread to the Far East, including the most south-eastern edge of its native habitats, the Yangtze River basin, where the climate is very different.

Results

We sequenced 118 A. thaliana strains from the region surrounding the Yangtze River basin. We found that the Yangtze River basin population is a unique population and diverged about 61,409 years ago, with gene flows occurring at two different time points, followed by a population dispersion into the Yangtze River basin in the last few thousands of years. Positive selection analyses revealed that biological regulation processes, such as flowering time, immune and defense response processes could be correlated with the adaptation event. In particular, we found that the flowering time gene SVP has contributed to A. thaliana adaptation to the Yangtze River basin based on genetic mapping.

Conclusions

A. thaliana adapted to the Yangtze River basin habitat by promoting the onset of flowering, a finding that sheds light on how a species can adapt to locales with very different climates.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1378-9) contains supplementary material, which is available to authorized users.

Detection of Pathways Affected by Positive Selection in Primate Lineages Ancestral to Humans

Gene set enrichment approaches have been increasingly successful in finding signals of recent polygenic selection in the human genome. In this study, we aim at detecting biological pathways affected by positive selection in more ancient human evolutionary history. Focusing on four branches of the primate tree that lead to modern humans, we tested all available protein coding gene trees of the Primates clade for signals of adaptation in these branches, using the likelihood-based branch site test of positive selection. The results of these locus-specific tests were then used as input for a gene set enrichment test, where whole pathways are globally scored for a signal of positive selection, instead of focusing only on outlier "significant" genes. We identified signals of positive selection in several pathways that are mainly involved in immune response, sensory perception, metabolism, and energy production. These pathway-level results are highly significant, even though there is no functional enrichment when only focusing on top scoring genes. Interestingly, several gene sets are found significant at multiple levels in the phylogeny, but different genes are responsible for the selection signal in the different branches. This suggests that the same function has been optimized in different ways at different times in primate evolution.

Frequency of adult type-associated lactase persistence LCT-13910C/T genotypes in the Czech/Slav and Czech Roma/Gypsy populations

Lactase non-persistence (leading to primary lactose intolerance) is a genetically dependent inability to digest lactose in adulthood. As part of the human adaptation to dairying, the human lactase LCT-13910C/T mutation (which propagates adult expression of lactase) developed, spread and participated in the adaptation to dairying. This variant is associated with lactase activity persistence, and its carriers are able to digest lactose. We compared the frequencies of lactase 13910C/T (rs4988235) genotypes in Czechs/Slavs (N = 288) and Czech Gypsies/Roma (N = 300), two ethnically different groups where this polymorphism has not yet been analysed. Allelic frequencies significantly differed between the populations (p < 0.0001). In Czechs/Slavs, the lactase persistence T allele was present in 76% of the individuals, which is in agreement with frequencies among geographically neighbouring populations. In the Czech Gypsy/Roma population, only 27% of the adults were carriers of at least one lactase persistence allele, similar to the Indian population. In agreement with this result, dairy product consumption was reported by 70.5% of Czechs/Slavs and 39.0% of the Czech Gypsy/Roma population. Both in the Czech Gypsy/Roma and in the Czech/Slavs populations, the presence of carriers of the lactase persistence allele was similar in subjects self-reporting the consumption of unfermented/fresh milk, in comparison to the others.

Autism spectrum disorder in sub-saharan africa: A comprehensive scoping review

Autism spectrum disorder (ASD) is recognized as a global public health concern, yet almost everything we know about ASD comes from high-income countries. Here we performed a scoping review of all research on ASD ever published in sub-Saharan Africa (SSA) in order to identify ASD knowledge gaps in this part of the world. Fifty-three publications met inclusion criteria. Themes included the phenotype, genetics and risk factors for ASD in SSA, screening and diagnosis, professional knowledge, interventions for ASD, parental perceptions, and social-cognitive neuroscience. No epidemiological, early intervention, school-based or adult studies were identified. For each identified theme, we aimed to summarize results and make recommendations to fill the knowledge gaps. The quality of study methodologies was generally not high. Few studies used standardized diagnostic instruments, and intervention studies were typically small-scale. Overall, findings suggest a substantial need for large-scale clinical, training, and research programmes to improve the lives of people who live with ASD in SSA. However, SSA also has the potential to make unique and globally-significant contributions to the etiology and treatments of ASD through implementation, interventional, and comparative genomic science. Autism Res 2017, 10: 723-749. © 2017 International Society for Autism Research, Wiley Periodicals, Inc.

Going global by adapting local: A review of recent human adaptation

The spread of modern humans across the globe has led to genetic adaptations to diverse local environments. Recent developments in genomic technologies, statistical analyses, and expanded sampled populations have led to improved identification and fine-mapping of genetic variants associated with adaptations to regional living conditions and dietary practices. Ongoing efforts in sequencing genomes of indigenous populations, accompanied by the growing availability of "-omics" and ancient DNA data, promises a new era in our understanding of recent human evolution and the origins of variable traits and disease risks.

Comparison of Human Neonatal and Adult Blood Leukocyte Subset Composition Phenotypes

The human peripheral leukocyte subset composition depends on genotype variation and pre-natal and post-natal environmental influence diversity. We quantified this composition in adults and neonates, and compared the median values and dispersal ranges of various subsets in them. We confirmed higher frequencies of monocytes and regulatory T cells (Tregs), similar frequencies of neutrophils, and lower frequencies of CD8 T cells, NKT cells, B1 B cells and gamma-delta T cells in neonatal umbilical cord blood. Unlike previous reports, we found higher frequencies of eosinophils and B cells, higher CD4:CD8 ratios, lower frequencies of T cells and iNKT cells, and similar frequencies of CD4 T cells and NK cells in neonates. We characterized monocyte subsets and dendritic cell (DC) subsets in far greater detail than previously reported, using recently described surface markers and gating strategies and observed that neonates had lower frequencies of patrolling monocytes and lower myeloid dendritic cell (mDC):plasmacytoid DC (pDC) ratios. Our data contribute to South Asian reference values for these parameters. We found that dispersal ranges differ between different leukocyte subsets, suggesting differential determination of variation. Further, some subsets were more dispersed in adults than in neonates suggesting influences of postnatal sources of variation, while some show the opposite pattern suggesting influences of developmental process variation. Together, these data and analyses provide interesting biological possibilities for future exploration.

Development of a Comprehensive Genotype-to-Fitness Map of Adaptation-Driving Mutations in Yeast

Adaptive evolution plays a large role in generating the phenotypic diversity observed in nature, yet current methods are impractical for characterizing the molecular basis and fitness effects of large numbers of individual adaptive mutations. Here, we used a DNA barcoding approach to generate the genotype-to-fitness map for adaptation-driving mutations from a Saccharomyces cerevisiae population experimentally evolved by serial transfer under limiting glucose. We isolated and measured the fitness of thousands of independent adaptive clones and sequenced the genomes of hundreds of clones. We found only two major classes of adaptive mutations: self-diploidization and mutations in the nutrient-responsive Ras/PKA and TOR/Sch9 pathways. Our large sample size and precision of measurement allowed us to determine that there are significant differences in fitness between mutations in different genes, between different paralogs, and even between different classes of mutations within the same gene.

Comparative Methylome Analyses Identify Epigenetic Regulatory Loci of Human Brain Evolution

How do epigenetic modifications change across species and how do these modifications affect evolution? These are fundamental questions at the forefront of our evolutionary epigenomic understanding. Our previous work investigated human and chimpanzee brain methylomes, but it was limited by the lack of outgroup data which is critical for comparative (epi)genomic studies. Here, we compared whole genome DNA methylation maps from brains of humans, chimpanzees and also rhesus macaques (outgroup) to elucidate DNA methylation changes during human brain evolution. Moreover, we validated that our approach is highly robust by further examining 38 human-specific DMRs using targeted deep genomic and bisulfite sequencing in an independent panel of 37 individuals from five primate species. Our unbiased genome-scan identified human brain differentially methylated regions (DMRs), irrespective of their associations with annotated genes. Remarkably, over half of the newly identified DMRs locate in intergenic regions or gene bodies. Nevertheless, their regulatory potential is on par with those of promoter DMRs. An intriguing observation is that DMRs are enriched in active chromatin loops, suggesting human-specific evolutionary remodeling at a higher-order chromatin structure. These findings indicate that there is substantial reprogramming of epigenomic landscapes during human brain evolution involving noncoding regions.

Evidence of Polygenic Adaptation in the Systems Genetics of Anthropometric Traits

Many signals of natural selection have been identified in the human genome. However, except for some single-locus mechanisms, most molecular processes generating these adaptation signals are still unknown. We developed an approach that integrates datasets related to genome-wide association studies (GWAS) with information about systems biology and genetic signatures of natural selection to identify evidence of polygenic adaptation. Specifically, we focused on five anthropometric measurements: body mass index (BMI), height, waist-to-hip ratio adjusted for BMI (WHR), and waist circumference adjusted for BMI (WC), and sex differences for WHR and WC. We performed an enrichment analysis for signals of natural selection in protein interaction networks associated with anthropometric traits in European populations. The adaptation signals-enriched gene networks associated highlighted epistatic interactions in the context of polygenic selection for the investigated traits. These polygenic mechanisms indicated intriguing selective mechanisms related to the anthropometric traits: adult locomotory behavior for BMI, infection resistance for height, interplay between lipid transport and immune systems for WHR, and female-specific polygenic adaptation for WHR and WC. In conclusion, we observed evidence of polygenic adaptation in the context of systems genetics of anthropometric traits that indicates polygenic mechanisms related to the natural selection in European populations.

Population-genetic properties of differentiated copy number variations in cattle

While single nucleotide polymorphism (SNP) is typically the variant of choice for population genetics, copy number variation (CNV) which comprises insertion, deletion and duplication of genomic sequence, is an informative type of genetic variation. CNVs have been shown to be both common in mammals and important for understanding the relationship between genotype and phenotype. However, CNV differentiation, selection and its population genetic properties are not well understood across diverse populations. We performed a population genetics survey based on CNVs derived from the BovineHD SNP array data of eight distinct cattle breeds. We generated high resolution results that show geographical patterns of variations and genome-wide admixture proportions within and among breeds. Similar to the previous SNP-based studies, our CNV-based results displayed a strong correlation of population structure and geographical location. By conducting three pairwise comparisons among European taurine, African taurine, and indicine groups, we further identified 78 unique CNV regions that were highly differentiated, some of which might be due to selection. These CNV regions overlapped with genes involved in traits related to parasite resistance, immunity response, body size, fertility, and milk production. Our results characterize CNV diversity among cattle populations and provide a list of lineage-differentiated CNVs.

The Relevance of Genomic Signatures at Adhesion GPCR Loci in Humans

Adhesion G protein-coupled receptors (aGPCRs) have a long evolutionary history dating back to very basal unicellular eukaryotes. Almost every vertebrate is equipped with a set of different aGPCRs. Genomic sequence data of several hundred extinct and extant species allows for reconstruction of aGPCR phylogeny in vertebrates and non-vertebrates in general but also provides a detailed view into the recent evolutionary history of human aGPCRs. Mining these sequence sources with bioinformatic tools can unveil many facets of formerly unappreciated aGPCR functions. In this review, we extracted such information from the literature and open public sources and provide insights into the history of aGPCR in humans. This includes comprehensive analyses of signatures of selection, variability of human aGPCR genes, and quantitative traits at human aGPCR loci. As indicated by a large number of genome-wide genotype-phenotype association studies, variations in aGPCR contribute to specific human phenotypes. Our survey demonstrates that aGPCRs are significantly involved in adaptation processes, phenotype variations, and diseases in humans.

Variation in NAT2 acetylation phenotypes is associated with differences in food-producing subsistence modes and ecoregions in Africa

Background

Dietary changes associated to shifts in subsistence strategies during human evolution may have induced new selective pressures on phenotypes, as currently held for lactase persistence. Similar hypotheses exist for arylamine N-acetyltransferase 2 (NAT2) mediated acetylation capacity, a well-known pharmacogenetic trait with wide inter-individual variation explained by polymorphisms in the NAT2 gene. The environmental causative factor (if any) driving its evolution is as yet unknown, but significant differences in prevalence of acetylation phenotypes are found between hunter-gatherer and food-producing populations, both in sub-Saharan Africa and worldwide, and between agriculturalists and pastoralists in Central Asia. These two subsistence strategies also prevail among sympatric populations of the African Sahel, but knowledge on NAT2 variation among African pastoral nomads was up to now very scarce. Here we addressed the hypothesis of different selective pressures associated to the agriculturalist or pastoralist lifestyles having acted on the evolution of NAT2 by sequencing the gene in 287 individuals from five pastoralist and one agriculturalist Sahelian populations.

Results

We show that the significant NAT2 genetic structure of African populations is mainly due to frequency differences of three major haplotypes, two of which are categorized as decreased function alleles (NAT2*5B and NAT2*6A), particularly common in populations living in arid environments, and one fast allele (NAT2*12A), more frequently detected in populations living in tropical humid environments. This genetic structure does associate more strongly with a classification of populations according to ecoregions than to subsistence strategies, mainly because most Sahelian and East African populations display little to no genetic differentiation between them, although both regions hold nomadic or semi-nomadic pastoralist and sedentary agriculturalist communities. Furthermore, we found significantly higher predicted proportions of slow acetylators in pastoralists than in agriculturalists, but also among food-producing populations living in the Sahelian and dry savanna zones than in those living in humid environments, irrespective of their mode of subsistence.

Conclusion

Our results suggest a possible independent influence of both the dietary habits associated with subsistence modes and the chemical environment associated with climatic zones and biomes on the evolution of NAT2 diversity in sub-Saharan African populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0543-6) contains supplementary material, which is available to authorized users.

The epigenomic landscape of African rainforest hunter-gatherers and farmers

The genetic history of African populations is increasingly well documented, yet their patterns of epigenomic variation remain uncharacterized. Moreover, the relative impacts of DNA sequence variation and temporal changes in lifestyle and habitat on the human epigenome remain unknown. Here we generate genome-wide genotype and DNA methylation profiles for 362 rainforest hunter-gatherers and sedentary farmers. We find that the current habitat and historical lifestyle of a population have similarly critical impacts on the methylome, but the biological functions affected strongly differ. Specifically, methylation variation associated with recent changes in habitat mostly concerns immune and cellular functions, whereas that associated with historical lifestyle affects developmental processes. Furthermore, methylation variation—particularly that correlated with historical lifestyle—shows strong associations with nearby genetic variants that, moreover, are enriched in signals of natural selection. Our work provides new insight into the genetic and environmental factors affecting the epigenomic landscape of human populations over time.

Genetic and environmental factors affect genome-wide patterns of epigenetic variation. Here, the authors show that while current habitat and historical lifestyle impact the methylome of rainforest hunter-gatherers and sedentary farmers, the biological functions affected and the degree of genetic control differ.

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

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

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

The dJ/dS Ratio Test Reveals Hundreds of Novel Putative Cancer Drivers

Computational tools with a balanced sensitivity and specificity in identification of candidate cancer drivers are highly desired. In this study, we propose a new statistical test, namely the dJ/dS ratio test, to compute the relative mutation rate of exon/intron junction sites (dJ) to synonymous sites (dS); observation of dJ/dS ratio larger than 1 in cancer indicates positive selection for splicing deregulation, a signature of cancer driver genes. Using this method, we analyzed the data from The Cancer Genome Atlas and identified hundreds of novel putative cancer drivers. Interestingly, these genes are highly enriched in biological processes related to the development and maintenance of multicellularity, paralleling a previous finding that cancer evolves back to be unicellular by knocking down the multicellularity-associated genetic network.

Towards an ecological understanding of morphological evolution

The roots of modern evo-devo can be traced back to the comparative anatomy of the 19th century. Inheriting from this tradition, the field has maintained a mechanistic approach to understanding the origins of distinct animal morphologies. While this focus has produced a valuable body of work, we argue here that a fuller understanding of why species diverge morphologically must be centered on the selective forces driving divergence, and these forces ultimately reside in the ecological context in which organisms live and reproduce. We discuss reasons why we expect many morphological novelties to evolve largely secondarily to, and often as a by-product of, primary selection on life-history traits. By shifting the focus to proximate evolutionary causes, our perspective necessarily prioritises selection experiments as a means of empirical testing. We outline experimental approaches designed to dissect the role of ecological variables in the evolution of animal development and morphology, and we show how methods and advances in fields as diverse as population genomics and ecological stoichiometry can contribute to progress in this direction.

Evidence for negative selection of gene variants that increase dependence on dietary choline in a Gambian cohort

Choline is an essential nutrient, and the amount needed in the diet is modulated by several factors. Given geographical differences in dietary choline intake and disparate frequencies of single-nucleotide polymorphisms (SNPs) in choline metabolism genes between ethnic groups, we tested the hypothesis that 3 SNPs that increase dependence on dietary choline would be under negative selection pressure in settings where choline intake is low: choline dehydrogenase (CHDH) rs12676, methylenetetrahydrofolate reductase 1 (MTHFD1) rs2236225, and phosphatidylethanolamine-N-methyltransferase (PEMT) rs12325817. Evidence of negative selection was assessed in 2 populations: one in The Gambia, West Africa, where there is historic evidence of a choline-poor diet, and the other in the United States, with a comparatively choline-rich diet. We used 2 independent methods, and confirmation of our hypothesis was sought via a comparison with SNP data from the Maasai, an East African population with a genetic background similar to that of Gambians but with a traditional diet that is higher in choline. Our results show that frequencies of SNPs known to increase dependence on dietary choline are significantly reduced in the low-choline setting of The Gambia. Our findings suggest that adequate intake levels of choline may have to be reevaluated in different ethnic groups and highlight a possible approach for identifying novel functional SNPs under the influence of dietary selective pressure.—Silver, M. J., Corbin, K. D., Hellenthal, G., da Costa, K.-A., Dominguez-Salas, P., Moore, S. E., Owen, J., Prentice, A. M., Hennig, B. J., Zeisel, S. H. Evidence for negative selection of gene variants that increase dependence on dietary choline in a Gambian cohort.

Evidence for recent, population-specific evolution of the human mutation rate

As humans dispersed out of Africa they adapted to new environmental challenges, including changes in exposure to mutagenic solar radiation. Humans in temperate latitudes have acquired light skin that is relatively transparent to UV light, and some evidence suggests that their DNA damage response pathways have also experienced local adaptation. This raises the possibility that different populations have experienced different selective pressures affecting genome integrity. Here, I present evidence that the rate of a particular mutation type has recently increased in the European population, rising in frequency by 50% during the 40,000-80,000 y since Europeans began diverging from Asians. A comparison of SNPs private to Africa, Asia, and Europe in the 1000 Genomes data reveals that private European variation is enriched for the transition 5'-TCC-3' → 5'-TTC-3'. Although it is not clear whether UV played a causal role in changing the European mutational spectrum, 5'-TCC-3' → 5'-TTC-3' is known to be the most common somatic mutation present in melanoma skin cancers, as well as the mutation most frequently induced in vitro by UV. Regardless of its causality, this change indicates that DNA replication fidelity has not remained stable even since the origin of modern humans and might have changed numerous times during our recent evolutionary history.