Population Genetics: Global clues to the nature of genomic mutations in humans

An analysis of worldwide human genetic variation reveals the footprints of ancient changes in genomic mutation processes.
  1. Aylwyn Scally  Is a corresponding author
  1. University of Cambridge, United Kingdom

All children inherit a mixture of chromosomal sequences from their parents, and although the copying process involved is extremely accurate, some errors occur. We refer to such errors as de novo (new) germline mutations, and when they are passed on to subsequent generations, these mutations are the raw material on which natural selection works and the source of all genetic differences between populations and species. Most have negligible or minor effects, but some on very rare occasions are responsible for serious genetic disease (Veltman and Brunner, 2012). The accumulation of genetic differences through mutation is also a primary source of information about human evolution. Thus it is important to understand the nature of genomic mutation, the rate at which it occurs and the factors causing it.

A key question is to what extent mutation processes differ between individuals, either in the total number of de novo mutations bequeathed to offspring or the places in the genome where they occur. Such differences could be genetic in origin, given that the proteins involved in DNA replication are themselves encoded in the genome, and there could also be environmental effects associated with where and how an individual lives. In both cases these factors may reflect recent evolutionary events, particularly the divergence of human populations and their global dispersal within the last 100,000 years.

One way to address this question would be to collect genomic and other data for thousands of families, identify de novo mutations in each of the offspring, and then analyse the factors contributing to them. However, this approach is extremely demanding in terms of the resources needed. Now, in eLife, Kelley Harris and Jonathan Pritchard of Stanford University report how they have taken an alternative approach (Harris and Pritchard, 2017) that involved using a dataset of whole-genome sequences for 2504 individuals from 26 different populations around the world (Auton et al., 2015). Consider that any genetic variant in this dataset, even if it is found today in many individuals, was once a de novo mutation in a single ancestor. Thus if different genetic or environmental factors have affected mutation processes in different human populations, we might expect to find evidence in the distribution of variants in these populations today.

Harris and Pritchard categorised single-nucleotide variants in this dataset by their ancestral and derived alleles (i.e. the version before and after mutation) and their sequence context as represented by the two flanking nucleotides. For example, the category AGC → ATC represents a mutation from G to T with flanking nucleotides A and C. After counting the number of variants in each category within every individual, the researchers found that the distribution of counts, termed the mutation spectrum, differs between populations to the extent that it is possible to identify an individual’s continent of origin based solely on the spectrum of mutations they carry. In general the differences between spectra comprise a multitude of small discrepancies, rather than large discrepancies in a few categories. However some categories do stand out, most notably an increased abundance of TCC → TTC mutations in European and South Asian populations – a signal also seen in other recent studies based on similar data (Mathieson and Reich, 2017; Narasimhan et al., 2016).

Harris and Pritchard then looked at how this signal changes with time. Mutations themselves have no timestamp, but on average a mutation that is rare in the population is likely to have arisen more recently than a mutation that is common. Using frequency as a proxy for age, Harris and Pritchard found that the enrichment for TCC → TTC was evident mainly in variants of intermediate age, and not in recent or very old variants. The data fit a model in which there was a pulse of mutation between about 2,000 and 15,000 years ago in the ancestors of present-day Europeans and South Asians, during which these mutations were 30–40% more likely.

One possible explanation is that this is the legacy of a mutator allele (that is, a mutation that increases the rate of de novo mutation in some or all categories) that appeared and survived in the population for several thousand years before going extinct (Figure 1). It is not known how frequently such alleles arise, or to what extent they drive differences in the mutation spectrum: however it does seem that under certain conditions they can survive for a long time and have lasting effects (Seoighe and Scally, 2017). In addition to a mutator allele that would have increased the relative rate of the TCC → TTC mutation, it is possible that other mutator alleles that had smaller effects and/or survived for shorter times in different populations may be responsible for less prominent differences in the mutation spectrum.

Changes in the genetic mutation spectrum in humans.

Modern humans originated in Africa (top left) and spread into Eurasia in a series of migrations 50,000–80,000 years ago. Subsequent migrations within Eurasia (such as the migration between West Eurasia and South Asia around 3000 years ago that is shown here; Reich et al., 2009) led to secondary genetic contact and admixture. The mutation spectrum (represented abstractly here) evolved separately in each population after divergence, perhaps due to the effect of mutator alleles: Harris and Pritchard propose that such a mutator allele could have been responsible for an increase in TCC → TTC mutations between 2,000 and 15,000 years ago (red shading) which influenced the mutation spectra of present-day Europeans and South Asians.

These findings are also relevant to the question of whether or not the overall rate of genomic mutation varies between populations and over time. There is evidence that the mean rate has probably not changed in at least the last 50,000 years (Fu et al., 2014), but we also know that the genome-wide mutation rate has slowed down when measured over timescales of millions of years (Moorjani et al., 2016). The results of Harris and Pritchard involve the relative rates of mutation in different sequence contexts, and so strictly speaking they do not tell us about variation in the overall rate, but they are suggestive of differences perhaps on the order of a few percent.

The next few years will see a substantial increase in the amount of available de novo mutation data, particularly from large-scale sequencing projects aimed at understanding the causes of rare genetic disease. These data will enable direct exploration of the factors determining mutation and its evolution within human populations, and shed further light on the questions addressed and raised by Harris and Pritchard.


Article and author information

Author details

  1. Aylwyn Scally

    Department of Genetics, University of Cambridge, Cambridge, United Kingdom
    For correspondence
    Competing interests
    The author declares that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-0807-1167

Publication history

  1. Version of Record published: May 17, 2017 (version 1)


© 2017, Scally

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.


  • 2,152
  • 257
  • 2

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Aylwyn Scally
Population Genetics: Global clues to the nature of genomic mutations in humans
eLife 6:e27605.
  1. Further reading

Further reading

    1. Biochemistry and Chemical Biology
    2. Evolutionary Biology
    Eva Pyrihová, Martin S King ... Edmund RS Kunji
    Research Article

    Stramenopiles form a clade of diverse eukaryotic organisms, including multicellular algae, the fish and plant pathogenic oomycetes, such as the potato blight Phytophthora, and the human intestinal protozoan Blastocystis. In most eukaryotes, glycolysis is a strictly cytosolic metabolic pathway that converts glucose to pyruvate, resulting in the production of NADH and ATP (Adenosine triphosphate). In contrast, stramenopiles have a branched glycolysis in which the enzymes of the pay-off phase are located in both the cytosol and the mitochondrial matrix. Here, we identify a mitochondrial carrier in Blastocystis that can transport glycolytic intermediates, such as dihydroxyacetone phosphate and glyceraldehyde-3-phosphate, across the mitochondrial inner membrane, linking the cytosolic and mitochondrial branches of glycolysis. Comparative analyses with the phylogenetically related human mitochondrial oxoglutarate carrier (SLC25A11) and dicarboxylate carrier (SLC25A10) show that the glycolytic intermediate carrier has lost its ability to transport the canonical substrates malate and oxoglutarate. Blastocystis lacks several key components of oxidative phosphorylation required for the generation of mitochondrial ATP, such as complexes III and IV, ATP synthase, and ADP/ATP carriers. The presence of the glycolytic pay-off phase in the mitochondrial matrix generates ATP, which powers energy-requiring processes, such as macromolecular synthesis, as well as NADH, used by mitochondrial complex I to generate a proton motive force to drive the import of proteins and molecules. Given its unique substrate specificity and central role in carbon and energy metabolism, the carrier for glycolytic intermediates identified here represents a specific drug and pesticide target against stramenopile pathogens, which are of great economic importance.

    1. Evolutionary Biology
    2. Genetics and Genomics
    Brian PH Metzger, Yeonwoo Park ... Joseph W Thornton
    Research Article

    A protein’s genetic architecture – the set of causal rules by which its sequence produces its functions – also determines its possible evolutionary trajectories. Prior research has proposed that the genetic architecture of proteins is very complex, with pervasive epistatic interactions that constrain evolution and make function difficult to predict from sequence. Most of this work has analyzed only the direct paths between two proteins of interest – excluding the vast majority of possible genotypes and evolutionary trajectories – and has considered only a single protein function, leaving unaddressed the genetic architecture of functional specificity and its impact on the evolution of new functions. Here, we develop a new method based on ordinal logistic regression to directly characterize the global genetic determinants of multiple protein functions from 20-state combinatorial deep mutational scanning (DMS) experiments. We use it to dissect the genetic architecture and evolution of a transcription factor’s specificity for DNA, using data from a combinatorial DMS of an ancient steroid hormone receptor’s capacity to activate transcription from two biologically relevant DNA elements. We show that the genetic architecture of DNA recognition consists of a dense set of main and pairwise effects that involve virtually every possible amino acid state in the protein-DNA interface, but higher-order epistasis plays only a tiny role. Pairwise interactions enlarge the set of functional sequences and are the primary determinants of specificity for different DNA elements. They also massively expand the number of opportunities for single-residue mutations to switch specificity from one DNA target to another. By bringing variants with different functions close together in sequence space, pairwise epistasis therefore facilitates rather than constrains the evolution of new functions.