% Use only LaTeX2e, calling the article.cls class and 12-point type.

\documentclass[12pt]{article}
\usepackage{amssymb}
\usepackage{amsmath}
\usepackage{graphicx}
\usepackage[margin=1in]{geometry}
\usepackage{setspace}

\usepackage{morefloats}
\usepackage{placeins}

\usepackage{lineno}
\linenumbers

%% make sure you have the nature.cls and naturemag.bst files where
%% LaTeX can find them

\bibliographystyle{naturemag}



\title{Rapid evolution of the human mutation spectrum.}
%\title{Temporary pulses of accelerated mutagenesis in human and great ape population history}
%\title{Rapid evolution of the mutation spectrum in humans and other apes}

%Supplementary Figure references
\newcommand{\allHumanHeatmap}{1}
\newcommand{\admixedPCA}{2}
\newcommand{\SGDPheatmaps}{3}
\newcommand{\SGDPregressions}{4}
\newcommand{\noBGC}{5}
\newcommand{\TCCpulseAllPops}{6}
\newcommand{\pulseSeriesStart}{7}
\newcommand{\pulseSeriesEnd}{8}
\newcommand{\fineHeatmap}{9}
\newcommand{\AmrPCA}{10}
\newcommand{\AfrPCA}{11}
\newcommand{\EasPCA}{12}
\newcommand{\SasPCA}{13}
\newcommand{\EurPCA}{14}
\newcommand{\allApeHeatmaps}{15}

%Main figure references
\newcommand{\overall}{1}

\author
{Kelley Harris$^{1}$ and Jonathan K. Pritchard$^{1,2,3}$\\
\\
\normalsize{$^{1}$Department of Genetics, Stanford University}\\
\normalsize{$^{2}$Department of Biology, Stanford University}\\
\normalsize{$^{3}$Howard Hughes Medical Institute, Stanford University}\\
\\
\normalsize{Correspondence:  kelleyh@stanford.edu, pritch@stanford.edu.}
}

% Include the date command, but leave its argument blank.

\date{\today}



%%%%%%%%%%%%%%%%% END OF PREAMBLE %%%%%%%%%%%%%%%%



\begin{document} 
% Double-space the manuscript.
\baselineskip 20pt %24pt

\maketitle 
\clearpage

\begin{abstract}

\noindent DNA is a remarkably precise medium for copying and storing biological information. This high fidelity results from 
the action of hundreds of genes involved in replication, proofreading, and damage repair. Evolutionary theory suggests that in such a system, selection has limited ability to remove genetic variants that change mutation rates by small amounts or in specific sequence contexts.
Consistent with this, using SNV variation as a proxy for mutational input, we report here that mutational spectra differ substantially among species, human continental groups and even some closely-related 
populations. Close examination of one signal, an increased TCC$\to$TTC mutation rate in Europeans, indicates a burst of mutations from about 15,000 to 2,000 years ago, perhaps due to the appearance, drift, and ultimate elimination of a genetic modifier of mutation rate.  Our results suggest that mutation rates can evolve markedly over short evolutionary timescales and suggest the possibility of mapping mutational modifiers.


\end{abstract}

\clearpage

\section*{Main Text}

Germline mutations provide the raw material for evolution, but also generate genetic load and inherited disease.  Indeed, the vast majority of mutations that affect fitness are deleterious, and hence biological systems have evolved elaborate mechanisms for accurate DNA replication and repair of diverse types of spontaneous damage.  Due to the combined action of hundreds of genes, mutation rates are extremely low--in humans, about 1 point mutation per 100MB or about 60 genome-wide per generation \cite{kong2012, segurel2014}.

While the precise roles of most of the relevant genes have not been fully elucidated, research on somatic mutations in cancer has shown that defects in particular genes can lead to increased mutation rates within very specific sequence contexts \cite{alexandrov2013, helleday2014}.  For example, mutations in the proofreading exonuclease domain of DNA polymerase $\epsilon$ cause TCT$\to$TAT and TCG$\to$TTG mutations on the leading DNA strand \cite{Shinbrot2014}.  Mutational shifts of this kind have been referred to as ``mutational signatures''.  Specific signatures may also be caused by nongenetic factors such as chemical mutagens, UV damage, or guanine oxidation \cite{ohno2014}.  

Together, these observations imply a high degree of specialization of individual genes involved in DNA proofreading and repair.  While the repair system has evolved to be extremely accurate overall, theory suggests that in such a system, natural selection may have limited ability to fine-tune the efficacy of individual genes \cite{lynch2011, sung2012}.  If a variant in a repair gene increases or decreases the overall mutation rate by a small amount--for example, only in a very specific sequence context--then the net effect on fitness may fall below the threshold at which natural selection is effective.  (Drift tends to dominate selection when the change in fitness is less than the inverse of effective population size).  The limits of selection on mutation rate modifiers are especially acute in recombining organisms such as humans because a variant that increases the mutation rate can recombine away from deleterious mutations it generates elsewhere in the genome.

Given these theoretical predictions, we hypothesized that there may be substantial scope for modifiers of mutation rates to segregate within human populations, or between closely related species. Most triplet sequence contexts have mutation rates that vary across the evolutionary tree of mammals \cite{hwang2004}, but evolution of the mutation spectrum over short time scales has been less well described.  Weak natural mutators have recently been observed in yeast \cite{bui2017} and inferred from human haplotype data \cite{seoighe2017}; if such mutators affect specific pathways of proofreading or repair, then we may expect shifts in the abundance of mutations within particular sequence contexts.  Indeed, one of us has recently identified a candidate signal of this type, namely an increase in TCC$\to$TTC transitions in Europeans, relative to other populations \cite{harris2015}; this was recently replicated
%  A recent preprint has replicated that result%, and reported an additional mutational signature enriched in South American populations 
\cite{mathieson2016}. Here we show that mutation spectrum change is much more widespread than these initial studies suggested: although the TCC$\to$TTC rate increase in Europeans was unusually dramatic, smaller-scale changes are so commonplace that almost every great ape species and human continental group has its own distinctive mutational spectrum. 

\section*{Results}

To investigate the mutational processes in different human populations, we classified each single nucleotide variants (SNV) in the 1000 Genomes Phase 3 data \cite{1000Genomes2015} in terms of its ancestral allele, derived allele, and 5' and 3' flanking nucleotides. We collapsed strand complements together to obtain 96 SNV categories. Since the detection of singletons may vary across samples, and because some singletons may result from cell-line or somatic mutations, we only considered variants seen in more than one copy. We further excluded
variants in annotated repeats (since read mapping error rates may be higher in such regions) and in PhyloP conserved regions (to avoid selectively constrained regions) \cite{pollard2010}.  From the remaining sites, we calculated the distribution of derived SNVs carried by each Phase 3 individual. We used this as a proxy for the mutational input spectrum in the ancestors of each individual. %Although it is conceivable that different populations could experience selective pressures that differ systematically by SNV triplet context, this seems \emph{a priori} less plausible than mutation spectrum evolution. 

To explore global patterns of the mutation spectrum, we performed principal component analysis (PCA) in which each individual was characterized simply by the fraction of their derived alleles in each of the 96 SNV categories (Fig.~\overall A). PCA is commonly applied to individual-level genotypes, in which case the PCs are usually highly correlated with geography \cite{novembre2008}.  Although the triplet mutation spectrum is an extremely compressed summary statistic compared to typical genotype arrays, we found that it contains sufficient information to reliably classify individuals by continent of origin. The first principal component separated Africans from non-Africans, and the second separated Europeans from East Asians, with South Asians and admixed native Americans (Figure 1--Figure Supplement 2) appearing intermediate between the two. 

\begin{figure}
\centering
\includegraphics[width=\textwidth]{../Figure2/Figure2.pdf}

\caption{ \it{\bf Global patterns of variation in SNV spectra.} {\bf A.} Principal Component Analysis of individuals according to the fraction of derived alleles that each individual carries in each of 96 mutational types.  {\bf B.} Heatmaps showing, for pairs of continental groups, the ratio of the proportions of SNVs in each of the 96 mutational types. Each block corresponds to one mutation type; within blocks, rows indicate the 5' nucleotide, and columns indicate the 3' nucleotide. Red colors indicate a greater fraction of a given mutation type in the first-listed group relative to the second. Points indicate significant contrasts at p $<10^{-5}$. See Figure Supplements 1, 2, and 3 for heatmap comparisons between additional population pairs as well as a description of PCA loadings and the $p$-values of all mutation class enrichments. Figure Supplement 4 demonstrates that these patterns are unlikely to be driven by biased gene conversion. In Figure Supplement 5, we see that this mutation spectrum structure replicates on both strands of the transcribed genome as well as the non-transcribed portion of the genome. Figure Supplements 6, 7, and 8 show that most of this structure replicates across multiple chromatin states and varies little with replication timing. }\label{overall}
\end{figure}

Remarkably, we found that the mutation spectrum differences among continental groups are composed of small shifts in the abundance of many different mutation types (Fig.~\overall B). For example, comparing Africans and Europeans, 43 of the 96 mutation types are significant at a $p<10^{-5}$ threshold using a forward variable selection procedure.   The previously described TCC$\to$TTC signature partially drives the difference between Europeans and the other groups, but most other shifts are smaller in magnitude and appear to be spread over more diffuse sets of related mutation types. East Asians have excess A$\to$T transversions in most sequence contexts, as well as about 10\% more *AC$\to$*CC mutations than any other group. Compared to Africans, all Eurasians have proportionally fewer C$\to$* mutations relative to A$\to$* mutations. 

\paragraph{Replication of mutation spectrum shifts.}
One possible concern is that batch effects or other sequencing artifacts might contribute to differences in mutational spectra. Therefore we replicated our analysis using 201 genomes from the Simons Genome Diversity Project \cite{mallick2016}.  The SGDP genomes were sequenced at high coverage, independently from 1000 Genomes, using an almost non-overlapp\-ing panel of samples. We found extremely strong agreement between the mutational shifts in the two data sets (Fig. \ref{simons_comparison}).  For example, all of the 43 mutation types with a significant difference between Africa and Europe (at $p<10^{-5}$) in 1000 Genomes also show a frequency difference in the same direction in SGDP (comparing Africa and West Eurasia). 
In both 1000 Genomes and SGDP, the enrichment of *AC$\to$*CC mutations in East Asia is larger in magnitude than any other signal aside from the previously described TCC$\to$TTC imbalance. 

The greatest discrepancies between 1000 Genomes and SGDP involve transversions at CpG sites, which are among the rarest mutational classes. These discrepancies might result from data processing differences or random sampling variation, but might also reflect differences in the fine-scale ethnic composition of the two panels. 

\begin{figure}
\includegraphics[width=\textwidth]{../Figure_Simons/updated_simons_figure.pdf}
\caption{\it{\bf Concordance of mutational shifts in 1000 Genomes vs.~SGDP.} Each panel shows natural-log mutation spectrum ratios between a pair of continental groups, based on 1000 Genomes (x-axis) and SGDP (y-axis) data. Data points encoded by (+) symbols denote mutation types that are not significantly enriched in either population in the Figure \overall~ 1000 Genomes analysis ($p<10^{-5}$). These heatmaps use the same labeling and color scale as in Figure \overall.  All 1000 Genomes ratios in this figure were estimated after projecting the 1000 Genomes site frequency spectrum down to the sample size of SGDP.  See Figure Supplements 1 and 2 for a complete set of SGDP heatmaps and regressions versus 1000 Genomes.}\label{simons_comparison}
\end{figure}

\paragraph{Evidence for a pulse of TCC$\to$TTC mutations in Europe and South Asia.}

To investigate the timescale over which the mutation spectrum change occurred, we analyzed the allele frequency distribution of TCC$\to$TTC mutations, which are highly enriched in Europeans (Fig.~3A; $p<1\times10^{-300}$ for Europe vs. Africa) and to a lesser extent in South Asians. We calculated allele frequencies both in 1000 Genomes and in the larger UK10K genome panel \cite{UK10K2015}.  As expected for a signal that is primarily European, we found particular enrichment of these mutations at low frequencies.  But surprisingly, the enrichment peaks around 0.6\% frequency in UK10K, and there is practically no enrichment among the very lowest frequency variants (Figure~3B and Figure 3--Figure Supplement 1). C$\to$T mutations on other backgrounds, namely within TCT, CCC and ACC contexts, 
are also enriched in Europe and South Asia and show a similar enrichment around 0.6\% frequency that declines among rarer variants (Fig.~3C). This suggests that these four mutation types comprise the signature of a single mutational pulse that is no longer active. No other mutation types show such a pulse-like distribution in UK10K, though several types show evidence of monotonic rate change over time (Figure 3--Figure Supplements 3,4 and 5).

\begin{figure}
\centering
\includegraphics[width=1.1\textwidth]{../Figure1/Figure1.pdf}

\caption{\it{\bf Geographic distribution and age of the TCC mutation pulse.} {\bf (A)} Observed frequencies of TCC$\to$TTC variants in 1000 Genomes populations. {\bf (B)} Fraction of TCC$\to$TTC variants as a function of allele frequency in different samples indicates that these peak around 1\%. See Figure Supplement 1 for distributions of TCC$\to$TTC allele frequency within all 1000 Genomes populations, and see Figure Supplement 2 for the replication of this result in the Exome Aggregation Consortium Data. In the UK10K data, which has the largest sample size, the peak occurs at 0.6\% allele frequency. {\bf (C)} Other enriched C$\to$T mutations with similar context also peak at 0.6\% frequency in UK10K. See Figure Supplements 3, 4 and 5 for labeled allele frequency distributions of all 96 mutation types (most represented here as unlabeled grey lines). See Figure Supplement 6 for heatmap comparisons of the 1000 Genomes populations partitioned by allele frequency, which provide a different view of these patterns. {\bf (D)} A population genetic model supports a pulse of TCC$\to$TTC mutations from 15,000--2,000 years ago. Inset shows the observed and predicted frequency distributions of this mutation under the inferred model.}
\label{TCC}
\end{figure}

%\paragraph{Inferring the age of the European mutation pulse.}
We used the enrichment of TCC$\to$TTC mutations as a function of allele frequency to estimate when this mutation pulse was active.  Assuming a simple piecewise-constant model, we infer that the rate of TCC$\to$TTC mutations increased dramatically $\sim$15,000 years ago and decreased again $\sim$2,000 years ago.  This time-range is consistent with results showing this signal in a pair of prehistoric European samples from 7,000 and 8,000 years ago, respectively \cite{mathieson2016}. We hypothesize that this mutation pulse may have been caused by a mutator allele that drifted up in frequency starting 15,000 years ago, but that is now rare or absent from present day populations.


Although low frequency allele calls often contain a higher proportion of base calling errors than higher frequency allele calls do, it is not plausible that base-calling errors could be responsible for the pulse we have described. In the UK10K data, a minor allele present at 0.6\% frequency corresponds to a derived allele that is present in 23 out of 3854 sampled haplotypes and supported by 80 short reads on average (assuming 7x coverage per individual). When independently generated datasets of different sizes are projected down to the same sample size, the TCC$\to$TTC pulse spans the same range of allele frequencies in both datasets (Figure 3--Figure Supplements 1 and 2), which would not be the case if the shape of the curve were a function of low frequency errors.

%Although low frequency allele calls often contain a higher proportion of base calling errors than higher frequency allele calls do, it is not plausible that base-calling errors could be responsible for the pulse we have described. In the UK10K data, a minor allele present at 0.6\% frequency corresponds to a derived allele that is present in 23 out of 3854 sampled haplotypes and supported by 80 short reads on average (assuming 7x coverage per individual). Alleles that are only present at 0.3\% frequency, which have a significantly lower fraction of TCC$\to$TTC variants, are still supported by 40 short reads on average and thus unlikely to contain a high percentage of sequencing errors. When UK10K and 1000 Genomes are projected down to the same sample size, the TCC$\to$TTC pulse spans the same range of allele frequencies in both datasets (Figure 3--Figure Supplement 1), which would not be the case if the shape of the curve were a function of low-frequency errors. Furthermore, the same pulse pattern can be observed in the Exome Aggregation Consortium dataset of 60,706 exomes (Figure 3--Figure Supplement 2) \cite{lek2016}. When we partition the 1000 Genomes variants by allele frequency and generate separate heatmaps for rare and more common alleles, the lowest frequency variants ($<$0.5\%) are not less differentiated between populations than alleles of frequency 0.5\%--1\%, as they would be if most of these variants were artifactual. For example, the rarest alleles are more enriched in East Asia for CC$\to$AC, CG$\to$TG, and AC$\to$TC than is observed at higher frequencies (Figure 3--Figure Supplement 6). 

\paragraph{Fine-scale mutation spectrum variation in other populations.}
Encouraged by these results, we sought to find other signatures of recent mutation pulses.  We generated heatmaps and PCA plots of mutation spectrum variation within each continental group, looking for fine-scale differences between closely related populations (Figure 4 and Figure 4--Figure Supplements 1 through 6). In some cases mutational spectra differ even between very closely related populations. For example, the *AC$\to$*CC mutations with elevated rates in East Asia appear to be distributed heterogeneously within that group, with most of the load carried by a subset of the Japanese individuals. These individuals also have elevated rates of ACA$\to$AAA and TAT$\to$TTT mutations (Figure 4A and Figure Supplement 4). This signature appears to be present in only a handful of Chinese individuals, and no Kinh or Dai individuals.  As seen for the European TCC mutation, the enrichment of these mutation types peaks at low frequencies, i.e., $\sim$1\%. Given the availability of only 200 Japanese individuals in 1000 Genomes, it is hard to say whether the true peak is at a frequency much lower than 1\%.

PCA reveals relatively little fine-scale structure within the mutational spectra of Europeans or South Asians (Figure 4--Figure Supplements 5 and 6). However, Africans exhibit some substructure (Figure 4--Figure Supplement 3), with the Luhya exhibiting the most distinctive mutational spectrum. Unexpectedly, a closer examination of PC loadings reveals that the Luhya outliers are enriched for the same mutational signature identified in the Japanese. Even in Europeans and South Asians, the first PC is heavily weighted toward *AC$\to$*CC, ACA$\to$AAA, and TAT$\to$TTT, although this signature explains less of the mutation spectrum variance within these more homogeneous populations. The sharing of this signature may suggest either parallel increases of a shared mutation modifier, or a shared aspect of environment or life history that affects the mutation spectrum.

\begin{figure}
%\vskip -2cm
%\hskip -4cm
\includegraphics[width=0.8\textwidth]{../Figure3/Figure3.pdf}
%\vskip -1cm
\caption{ \it{\bf Mutational variation among east Asian populations.} {\bf (A)} PCA of east Asian samples from 1000 Genomes, based on the relative proportions of each of the 96 mutational types. See Figure Supplements 2 through 6 for other finescale population PCAs. {\bf (B)} Heatmaps showing, for pairs of east Asian samples, the ratio of the proportions of SNVs in each of the 96 mutational types. Points indicate significant contrasts at p $<10^{-5}$. See Figure Supplement 1 for additional finescale heatmaps. {\bf (C)} and {\bf (D)} Relative enrichment of each mutational type in Japanese and Dai, respectively as a function of allele frequency.  Six mutation types that are enriched in JPT are indicated. Populations: CDX=Dai, CHB=Han (Beijing); CHS=Han (south China); KHV=Kinh; JPT=Japanese.}
\label{Japanese}
\end{figure}


\paragraph{Mutation spectrum variation among the great apes.}

Finally, given our finding of extensive fine-scale variation in mutational spectra between human populations, we hypothesized that mutational variation between species is likely to be even greater. %Mutation spectrum differences between diverged orders of mammals are so pronounced that they were detectable in a 1.7 megabase genomic region before the invention of next-generation sequencing \cite{hwang2004}. More recent analyses of primate reference genomes have revealed pervasive heterogeneity of the mutation rate and spectrum \cite{moorjani2016}, confirming previous indications that species closely related to humans have lower mutation rates than species more distantly related to humans \cite{goodman1961,li1987,scally2012}. 
To compare the mutation spectra of the great apes in more detail, we obtained SNV data from the Great Ape Diversity Panel, which includes 78 whole genome sequences from six great ape species including human  \cite{prado2013}. Overall, we find dramatic variation in mutational spectra among the great ape species (Figure 5 and Figure 5--Figure Supplement 1).  


%\begin{figure}
%%\hskip -1.5cm
%\includegraphics[width=1.0\textwidth]{../Figure4/Figure4.pdf}
%%\vspace{-5cm}
%%\vskip 1cm
%\caption{\it {\bf Mutational differences among the great apes.} {\bf (A)} Relative abundance of SNV types in 5 ape species compared to Bornean Orangutan; data from \cite{prado2013}. Boxes indicate labels in (B).  For additional comparisons see Figure Supplement 1. {\bf (B)} Schematic phylogeny of the great apes highlighting notable changes in SNV abundances.}
%\label{apes}
%\end{figure}

As noted previously \cite{moorjani2016b}, one major trend is a higher proportion of CpG mutations among the species closest to human, possibly reflecting lengthening generation time along the human lineage, consistent with 
previous indications that species closely related to humans have lower mutation rates than more distant species \cite{goodman1961,li1987,scally2012}. 
 However, most other differences are not obviously related to known processes such as biased gene conversion and generation time change.  The A$\to$T mutation rate appears to have sped up in the common ancestor of humans, chimpanzees, and bonobos, a change that appears consistent with a mutator variant that was fixed instead of lost. It is unclear whether this ancient A$\to$T speedup is related to the A$\to$T speedup in East Asians. Other mutational signatures appear on only a single branch of the great ape tree, such as a slowdown of A$\to$C mutations in gorillas. 

\section*{Discussion}

The widespread differences captured in Figures \overall ~and \ref{simons_comparison} may be footprints of allele frequency shifts affecting different mutator alleles.  But in principle, other genetic and non-genetic processes may also impact the observed mutational spectrum. First, biased gene conversion (BGC) tends to favor C/G alleles over A/T, and BGC is potentially more efficient in populations of large effective size compared to populations of smaller effective size \cite{galtier2001}. %Demographic forces like population bottlenecks and changes in generation time might be responsible for some of these mutation spectrum shifts. However, it is hard to imagine how these forces could be solely responsible for rich landscape of differences visualized in Figure 1B. Population bottlenecks increase the proportion of C/G$\to$A/T alleles relative to A/T$\to$C/G alleles by reducing the efficiency of biased gene conversion (BGC). Like all types of natural selection, BGC is potentially more efficient in populations of large effective size compared to populations of smaller effective size \cite{galtier2001}. 
However, despite the bottlenecks that are known to have affected Eurasian diversity, there is no clear trend of an increased fraction of C/G$\to$A/T relative to A/T$\to$C/G in non-Africans vs. Africans, or with distance from Africa (Figure 1--Figure Supplement 7), and previous studies have also found little evidence for a strong genome-wide effect of BGC on the mutational spectrum in humans and great apes \cite{do2015, moorjani2016b}. For these reasons, we think that evolution of the mutational process is a better explanation than BGC or selection for differences that have been observed between the spectra of ultra-rare singleton variants and older human genetic variation \cite{carlson2017}; 

It is also known that shifts in generation time or other life-history traits may affect mutational spectra, particularly for CpG transitions \cite{martin1993, amster2016}. Most CpG transitions result from spontaneous methyl-cytosine deamination as opposed to errors in DNA replication. Hence %, if deaminated cytosines are not efficiently repaired, 
the rate of CpG transitions is less affected by generation time than other mutations \cite{hwang2004,moorjani2016, gao2016}. We observe that Europeans have a lower fraction of CpG variants compared to Africans, East Asians and South Asians (Fig.~1B), consistent with a recent report of a lower rate of \emph{de novo} CCG$\to$CTG mutations in European individuals compared to Pakistanis \cite{narasimhan2016}. Such a pattern may be consistent with a shorter average generation time in Europeans \cite{moorjani2016}, though it is unclear that a plausible shift in generation-time could produce such a large effect. Apart from this, the other patterns evident in Figure 1 do not seem explainable by known processes. 

In summary, we report here that, mutational spectra differ significantly among closely related human populations, and that they differ greatly among the great ape species.  Our work shows that subtle, concerted shifts in the frequencies of many different mutation types are more widespread than dramatic jumps in the rate of single mutation types, although the existence of the European TCC$\to$TTC pulse shows that both modes of evolution do occur \cite{harris2015, moorjani2016, mathieson2016}.  

At this time, we cannot exclude a role for nongenetic factors such as changes in life history or mutagen exposure in driving these signals.  However given the sheer diversity of the effects reported here, it seems parsimonious to us to propose that most of this variation is driven by the appearance and drift of genetic modifiers of mutation rate.  This situation is perhaps reminiscent of the earlier observation that genome-wide recombination patterns are variable among individuals \cite{coop2008}, and ultimate discovery of PRDM9 \cite{baudat2010}; although in this case it is unlikely that a single gene is responsible for all signals seen here.  As large datasets of {\it de novo} mutations become available, it should be possible to map mutator loci genome-wide. 
%are sequenced in the near future, we will have greater power to test for mutator activity in cancer-associated candidates linked to APOBEC and other genes.  
In summary, our results suggest the likelihood that mutational modifiers are an important part of the landscape of human genetic variation. 

\section*{Acknowledgements}
This work was funded by NIH grants GM116381 and HG008140,  and by the Howard Hughes Medical Institute. We thank Jeffrey Spence and Yun S. Song for technical assistance. 
%sharing code for simulating genetic data given the assumption of a variable mutation rate.  
We also thank 
Ziyue Gao, Arbel Harpak, Molly Przeworski, Joshua Schraiber, and Aylwyn Scally for comments and discussion, as well as two anonymous reviewers. 

\newpage
\section*{Methods}

%\tableofcontents
%\newpage

\subsection*{Data Availability.}
All datasets analyzed here are publicly available at the following websites:

\begin{center}
\begin{tabular}{|c|c|}
\hline
1000 Genomes Phase 3 & http://www.1000genomes.org/category/phase-3/ \\
UK10K & http://www.uk10k.org/data-access.html \\
Simons Genome Diversity Panel & https://www.simonsfoundation.org/life-sciences/ \\
& simons-genome-diversity-project-dataset/ \\
\hline
\end{tabular}
\end{center}

\subsection*{Human Mutation Spectrum Processing.}

Mutation spectra were computed using 1000 Genomes Phase 3 SNPs \cite{1000Genomes2015} that are biallelic, pass all 1000 Genomes quality filters, and are not adjacent to any N's in the hg19 reference sequence. Ancestral states were assigned using the UCSC Genome Browser alignment of hg19 to the PanTro2 chimpanzee reference genome; SNPs were discarded if neither the reference nor alternate allele matched the chimpanzee reference. To minimize the potential impact of ancestral misidentification errors, SNPs with derived allele frequency higher than 0.98 were discarded. We also filtered out regions annotated as ``conserved" based on the 100-way PhyloP conservation score \cite{pollard2010}, download from http://hgdownload.cse. ucsc.edu/goldenPath/ hg19/phastCons100way/, as well as regions annotated as repeats by RepeatMasker \cite{smit2013}, downloaded from\\ 
http://hgdownload.cse.ucsc.edu/ goldenpath/hg19/database/nestedRepeats.txt.gz . To be counted as part of the mutation spectrum of population $P$ (which can be either a continental group or a finer-scale population from one city), a SNP should not be a singleton within population $P$--at least two copies of the ancestral and derived alleles must be present within that population.

An identical approach was used to extract the mutation spectrum of the UK10K ALSPAC panel \cite{UK10K2015}, which is not subdivided into smaller populations. The data were filtered as described in \cite{field2016}. The filtering procedure performed by Field, et al. reduces the ALSPAC sample size to 1927 individuals. 

We also computed mutation spectra of the Simons Genome Diversity Panel (SGDP) populations \cite{mallick2016}. Four of the SGDP populations, West Eurasia, East Asia, South Asia, and Africa, were compared to their direct counterparts in the 1000 Genomes data. Three additional SGDP populations, Central Asia and Siberia, Oceania, and America, had no close 1000 Genomes counterparts and were not analyzed here (although each project contained a panel of people from the Americans, the composition of the American panels was extremely different, with the 1000 Genomes populations being much more admixed with Europeans and Africans). SGDP sites with more than 20\% missing data were not utilized. All other data processing was done the same way described for the 1000 Genomes data. 

The following table gives the same size of each population panel, as well as the total number of SNPs segregating in the panel that are used to compute mutation type ratios:

\begin{center}
\begin{tabular}{|c|c|c|c|}
\hline
Dataset & Population & Number of individuals & Number of SNPs \\
\hline
1kg & Africa & 504 & 16,870,400 \\
1kg & Europe& 503 & 8,508,040 \\
1kg & East Asia & 504 & 7,895,925 \\
1kg & South Asia & 489 & 9,552,781 \\
SGDP & Africa & 45 & 6,569,658 \\
SGDP & West Eurasia & 69 & 4,201,571 \\
SGDP & East Asia & 49 & 3,312,645 \\
SGDP & South Asia & 38 & 3,449,624 \\
\hline
\end{tabular}
\end{center}

\subsection*{Great Ape Diversity Panel Data Processing.}

%\begin{minipage}{0.5\textwidth}
Biallelic great ape SNPs were extracted from the Great Ape Diversity Panel VCF \cite{prado2013}, which is aligned to the hg18 human reference sequence. Ancestral states were assigned using the Great Ape Genetic Diversity project annotation, which used the Felsenstein pruning algorithm to assign allelic states to internal nodes in the great ape tree. In the Great Ape Diversity Panel, the most recent common ancestor (MRCA) of the human species is labeled as node 18; the MRCAs of chimpanzees, bonobos, gorillas, and orangutans, respectively, are labeled as node 16, node 17, node 19, and node 15. We extracted the state of each MRCA at each SNP in the alignment and used it to polarize the ancestral and derived allele at that site; a SNP was discarded whenever the ancestral node was assigned an uncertain or polymorphic ancestral state. As with the human data, SNPs with derived allele frequency higher than 0.98 were not used, and both repeats and PhyloP-annotated conserved regions were filtered away. 
%\end{minipage}
%\begin{minipage}{0.49\textwidth}

%\end{minipage}

%\begin{figure}
%\centering
%\includegraphics[width=0.5\textwidth]{Ancestral_nodes.jpg}
%\caption[Phylogenetic tree of great ape relationships]{
%Phylogenetic tree of great ape relationships downloaded from the Great Ape Genome Project data server (\texttt{https://eichlerlab.gs.washington.edu /greatape/data/ }). Species name abbreviations: hom = human; ppa = bonobo; ptv, pte, ptt, pts = chimpanzee; ggd, ggg, gbg = gorilla; ponp = Bornean orangutan; pona = Sumatran orangutan.
%}\label{apetree}
%\end{figure}


\subsection*{Visual representation of mutation spectra.}

The mutation type of a SNP is defined in terms of its ancestral allele, its derived allele, and its two immediate 5' and 3' neighbors. Two mutation types are considered equivalent if they are strand-complementary to each other (e.g. ACG$\to$ATG is equivalent to CGT$\to$CAT). This scheme classifies SNPs into 96 different mutation types, each that can be represented with an A or C ancestral allele.  

To compute the frequency $f_P(m)$ of SNP $m$ in population $P$, we count up all SNPs of type $m$ where the derived allele is present in at least one representative of population $P$ (which can be either a specific population such as YRI or a broader continental group such as AFR). After obtaining this count $C_P(m)$, we define $f_P(m)$ to be the ratio $C_P(m)/\sum_{m'}C_P(m')$, where the sum in the denominator ranges over all 96 mutation types $m'$.  The enrichment of mutation type $m$ in population $P_1$ relative to population $P_2$ is defined to be $f_{P_1}(m)/f_{P_2}(m)$; these enrichments are visualized as heat maps in Figures 1B, 3B, and 4A. 

To track changes in the mutational spectrum over time, we compute $f_P(m)$ in bins of restricted allele frequency. This involves counting the number of SNPs of type $m$ that are present at frequency $\phi$ in population $P$ to obtain counts $C_P(m,\phi)$ and frequencies $f_P(m,\phi)=C_P(m,\phi)/\sum_{m'} C_P(m'\phi)$. Deviation of the ratio $f_P(m,\phi)/f_P(m)$ from 1 indicates that the rate of $m$ has fluctuated recently in the history of population $P$. To make the sampling noise approximately uniform across alleles of different frequencies, alleles of derived count greater than 5 were grouped into approximately log-spaced bins that each contained similar numbers of UK10K SNPs. More precisely, we defined a set of bin endpoints $b_1,b_2,...$ such that the total number of SNPs ranging in derived allele count between $b_i$ and $b_{i+1}-1$ is greater than or equal to the number of 5-ton SNPs, while the total number of SNPs ranging in derived allele count from $b_i$ to $b_{i+1}-2$ is less than the number of 5-ton SNPs.  

In some cases, e.g. Figures \ref{simons_comparison}, Figure 2--Figure Supplement 1B, and Figure 3--Figure Supplement 1, site frequency spectra were projected down to a smaller sample size before counting SNPs in order to more accurately compare datasets of different sample sizes. A binomial sampling approach was used to project a sample of $N$ haplotypes does to a smaller sample size $n$. Letting $C_P^{(N)}(m,\phi)$ denote the SNP counts in the large sample of $N$ haplotypes, effective SNP counts $C_P^{(n)}(m,\phi)$ in a sample of $n$ haplotypes are computed as follows:

$$C_P^{(n)}(m,k/n)=\binom{n}{k} \sum_{\ell=1}^{N-1} (\ell/N)^k(1-\ell/N)^{n-k} C_P^{(N)}(m,\ell/N)$$

\subsection*{Significance testing.}

One central goal of this paper is to test whether many mutation types differ in rate between human populations or whether mutation spectrum shifts have been rare events affecting only a small proportion of mutation types. A simple statistical method for answering this question would be to perform 96 separate chi-square tests, one for each triplet-context-dependent mutation type, as follows:

Let $S_i$ denote the total number of SNPs segregating in population $P_i$, and let $S_i^{(m)}$ denote the number of SNPs of mutation type $m$. If mutation type $m$ is more prevalent in population $P_1$ than in population $P_2$, a chi-square test provides a natural way of assessing the significance of this difference. As described in \cite{harris2015}, this test is performed on the following 2-by-2 contingency table:

\begin{center}
\begin{tabular}{|c|c|}
\hline
$S_1^{(m)}$ & $P_1-S_1^{(m)}$\\ 
\hline
$S_2^{(m)}$ & $P_2-S_2^{(m)}$\\ 
\hline
\end{tabular}
\end{center}

It would be appealing to conclude that every mutation type ``passing" this chi-square test is a mutation type that has changed in rate during recent human history. However, if we were to perform the full set of 96 tests, they would not be independent. A sufficiently large increase in the rate of one mutation type $m_1$ in population $P_1$ after divergence from $P_2$ could cause another mutation type $m_2$, whose rate has remained constant, to comprise significantly different fractions of the SNPs from $P_1$ and $P_2$. To minimize this effect, we formulate the following iterative procedure of conditionally independent tests: first, compute a chi-square significance value $p_{\mathrm{unordered}}(m)$ for each mutation type $m$ using the 2-by-2 chi-square table above. We then use these values to order the SNPs from lowest $p$ value to highest and compute a set of ordered $p$ values $p_{\mathrm{ordered}}(m)$. For the mutation type $m_0$ with the lowest unordered $p$ value, $p_{\mathrm{unordered}}(m_0)=p_{\mathrm{ordered}}(m_0)$. For mutation type $m_i$, which has the $i$th lowest unordered $p$ value and $i<96$, $p_{\mathrm{ordered}}(m_i)$ is computed from the following contingency table:

\begin{center}
\begin{tabular}{|c|c|}
\hline
$S_1^{(m_i)}$ & $\sum_{j=i+1}^{96} S_1^{(m_j)}$\\ 
\hline
$S_2^{(m_i)}$ & $\sum_{j=i+1}^{96} S_2^{(m_j)}$ \\ 
\hline
\end{tabular}
\end{center}

For mutation type $m_{96}$, which has the highest unordered $p$ value, the ordered $p$ value is computed from the contingency table 

\begin{center}
\begin{tabular}{|c|c|}
\hline
$S_1^{(m_{96})}$ & $S_1^{(m_{95})}$ \\ 
\hline
$S_2^{(m_{96})}$ &  $S_2^{(m_{95})}$ \\ 
\hline
\end{tabular}
\end{center}

This procedure is guaranteed to find fewer mutation types to differ significantly in rate between populations compared to separate chi-square tests.

\subsection*{Principal component analysis (PCA).}

The python package matplotlib.mlab.PCA was used to perform PCA on the complete set of 1000 Genomes haplotypes, each haplotype $h$ represented by a 96-element vector encoding the mutation frequencies $(f_h(m))_m$ of the non-singleton derived alleles present on that haplotype. In the same way, a separate PCA was performed on each of the 5 continental groups to reveal finescale components of mutation spectrum variation. 

\subsection*{Dating of the TCC$\to$T mutation pulse.}

We estimated the duration and intensity of TCC$\to$T rate acceleration in Europe by fitting a simple piecewise-constant rate model to the UK10K frequency data. To specify the parameters of the model, we divide time into discrete log-spaced intervals bounded by time points $t_1,...,t_d$, assigning each interval a TCC$\to$T mutation rate $r_0,...r_d$. In units of generations before the present, the time discretization points were chosen to be: 20, 40, 200, 400, 800, 1200, 1600, 2000, 2400, 2800, 3200, 3600, 4000, 8000, 12000, 16000, 20000, 24000, 28000, 32000, 36000, 40000. We assume that the total rate $r$ of mutations other than TCC$\to$T stays constant over time (a first-order approximation). 

In terms of these rate variables, we can calculate the expected shape of the TCC$\to$T pulse shown in Figure 2B of the main text. The shape of this curve depends on both the mutation rate parameters $r_i$ and the demographic history of the European population, which determines the joint distribution of allele frequency and allele age. To account for the effects of demography, we use Hudson's ms program to simulate 10,000 random coalescent trees under a realistic European demographic history inferred from allele frequency data \cite{tennessen2012} and condition our inference upon this collection of trees as follows:

Let $A(m,t)$ be the function for which $\int_{t_i}^{t_{i+1}} A(m,t)dt$ equals the coalescent tree branch length, averaged over the sample of simulated trees, that is ancestral to exactly $m$ lineages and falls between time $t_i$ and $t_{i+1}$. Given this function, which can be empirically estimated from a sample of simulated trees, the expected frequency spectrum entry $k/n$ is
$$E(k/n)=\frac{\sum_{i=1}^d\int_{t_{i-1}}^{t_i} A(k,t)dt}{\sum_{j=1}^n\sum_{i=1}^d\int_{t_{i-1}}^{t_i} A(j,t)dt}$$
and the expected fraction of TCC$\to$T mutations in allele frequency bin $k/n$ is
$$E(f_{\mathrm{TCC}\to\mathrm{T}}(k/n))=\frac{\sum_{i=1}^dr_i\int_{t_{i-1}}^{t_i} A(k,t)dt}{r\sum_{i=1}^d\int_{t_{i-1}}^{t_i} A(k,t)dt}.$$
The expected value of the TCC$\to$T enrichment ratio being plotted in Figure 2B is
$$E(r_{\mathrm{TCC}\to\mathrm{T}}(k/n))=\frac{\sum_{i=1}^dr_i\int_{t_{i-1}}^{t_i} A(k,t)dt\cdot \sum_{j=1}^n\sum_{i=1}^d\int_{t_{i-1}}^{t_i} A(j,t)dt}{\sum_{i=1}^d\int_{t_{i-1}}^{t_i} A(k,t)dt\cdot \sum_{j=1}^n\sum_{i=1}^dr_i\int_{t_{i-1}}^{t_i} A(j,t)dt}$$

In Figure 2B, enrichment ratios are not computed for every allele frequency in isolation, but for allele frequency bins that each contain similar numbers of SNPs. Given integers $1\leq k_m<k_{m+1}\leq n$, the expected TCC$\to$T enrichment ratio averaged over all SNPs with allele frequency between $k_m/n$ and $k_{m+1}/n$ is:
$$E(r_{\mathrm{TCC}\to\mathrm{T}}(k_m/n))=\frac{\sum_{i=1}^dr_i\int_{t_{i-1}}^{t_i} \sum_{k=k_m}^{k_{m+1}}A(k,t)dt\cdot \sum_{j=1}^n\sum_{i=1}^d\int_{t_{i-1}}^{t_i} A(j,t)dt}{\sum_{i=1}^d\int_{t_{i-1}}^{t_i} \sum_{k=k_m}^{k_{m+1}}A(k,t)dt\cdot \sum_{j=1}^n\sum_{i=1}^dr_i\int_{t_{i-1}}^{t_i} A(j,t)dt}$$

We optimize the mutation rates $r_1,\ldots,r_d$ using a log-spaced quantization of allele frequencies $k_1/n,\ldots,k_m/n$ defined such that all bins contain similar numbers of SNPs. The chosen allele count endpoints $k_1,\ldots,k_m$ are: 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 200, 300, 400, 500, 600, 700, 800, 900, 1000, 2000, 3000, 4000. Given this quantization of allele frequencies, we optimize $r_1,\ldots,r_d$ by using the BFGS algorithm to minimize the least squares distance $D(r_0,\ldots,r_d)$ between $E(r_{\mathrm{TCC}\to\mathrm{T}}(k_m/n))$ and the empirical ratio $r_{\mathrm{TCC}\to\mathrm{T}}(k_m/n)$ computed from the UK10K data. This optimization is subject to a regularization penalty that minimizes the jumps between adjacent mutation rates $r_i$ and $r_{i+1}$:
$$D(r_0,\ldots,r_d)= \sum_{m=1}^d \left( E(r_{\mathrm{TCC}\to\mathrm{T}}(k_m/n))-r_{\mathrm{TCC}\to\mathrm{T}}(k_m/n)\right)^2+0.25\sqrt{\sum_{i=1}^d (r_{i-1}-r_i)^2}$$

Although the underlying model of mutation rate change assumed here is very simple, it still represents an advance over the method used in \cite{harris2015} to estimate of the timing of the TCC$\to$TTC mutation rate increase. That method relied upon explicit estimates of allele age from a dataset of less than 100 individuals, which are much noisier than integration of a joint distribution of allele age and frequency across a sample of thousands of haplotypes. 

%\subsection*{Testing candidate loci for APOBEC-associated mutagenicity.} 

%A recent study found that two common germline variants affect both cancer risk and the somatic mutational spectrum by altering the regulation of APOBEC DNA-editing enzymes \cite{middlebrooks2016}. We tested each variant for germline mutator activity in the ALSPAC cohort of the UK10K data by looking at mutations that fall within 20 kB of the candidate mutator and have likely not had time to recombine away. We posited that a variant arose on the mutator background if at least 80\% of its derived alleles occur on a haplotype containing the candidate mutator; conversely, we posited that it arose on the ancestral background if less than 20\% of its derived alleles are currently linked to the candidate mutator. In genomes that are heterozygous for the putative mutator allele, singletons were assigned uniformly at random to the ancestral and derived allelic backgrounds. Based on previous descriptions of the APOBEC signature \cite{alexandrov2013,nikzainal2014}, we classified TC$\to$TT and TC$\to$TG mutations as potentially APOBEC-associated, all other mutations as non-APOBEC-associated, and tested for an enrichment of APOBEC-associated variants on the candidate mutator background. 

%The derived allele of SNP rs1014971, in the APOBEC3B regulatory region, is associated with an increased burden of somatic APOBEC mutations as well as increased risk of bladder cancer. We counted 34 APOBEC-type mutations linked to the derived allele in UK10K, as well as 206 non-APOBEC-type mutations. On the ancestral allele background, we counted 15 APOBEC-type mutations and 189 non-APOBEC-type mutations. By a $\chi^2$ test, the derived allele is significantly associated with an increased proportion of APOBEC-type mutations ($p=0.03$).

%A second APOBEC variant is associated with an increased risk of breast cancer: a 30 kB deletion that elides APOBEC3A and APOBEC3B into a single upregulated chimera. We attempted to test this deletion for germline mutagenicity by applying our testing procedure to the tag SNP rs12628403. One complicating factor is that the APOBEC deletion is rare in Europeans compared to East Asians, Amerindians, and Oceanians \cite{kidd2007}. In addition, low coverage sequencing datasets like the UK10K are vulnerable to assembly difficulties and a dearth of SNP calls near sites of structural variation.  As a result, only 19 SNPs, 1 APOBEC-type and 18 non-APOBEC-type, were found linked to the deletion within 20 kB, too few to reliably test for germline mutator activity.  

%\subsection{Code availability.} Contact K.H. (kelleyh@stanford.edu) to request copies of the python scripts that were used to perform the analyses in this paper.

\newpage
\section*{Supplementary Figures}

%\listoffigures
%\newpage

\renewcommand{\figurename}{Supplementary Figure S\hspace{-0.125 cm}}
\setcounter{figure}{0}

\renewcommand{\tablename}{Supplementary Table S\hspace{-0.125 cm}}
\setcounter{table}{0}

%\section{Comparisons between 1kG continental groups}

%Figure 1B of the main text shows heatmap comparisons among selected pairs of 1kG continental groups. For completeness, Figure S\ref{heatmap_all1kg} visualizes all possible comparisons between 1000 Genomes continental groups. Figure S\ref{PCA_wAMR} summarizes the first two principal components of variation within these groups, adding admixed American populations to Figure 1A, accompanied by heatmaps that visualize the loadings of the first two PCs.

%\begin{figure}
\begin{center}
%\includegraphics[width=0.75\textwidth]{../supplementary_material/heatmap_allgroups_pt1.pdf}
%\includegraphics[width=0.75\textwidth]{../supplementary_material/heatmap_allgroups_pt2.pdf}
\includegraphics[width=0.75\textwidth]{../final-eLife-fig-supplements/Figure-1-Supplement-1.pdf}
\end{center}

\vspace{-1.5cm}

\begin{singlespace}
\begin{flushleft}
{\bf Figure 1-- Figure Supplement 1: Pairwise mutation spectrum comparisons among continental groups.} Each of these plots compares the mutation spectra of two populations $P_1$ and $P_2$. Letting $f_i$ denote the fraction of SNVs in population $P_i$ that have a given triplet context, ancestral allele, and derived allele, the corresponding heat map square visualizes the enrichment ratio $f_1/f_2$. Black dots mark mutation types for which the difference between populations has a $\chi^2$ $p$-value less than $10^{-5}$.
\end{flushleft}
\end{singlespace}

%\caption[Pairwise mutation spectrum comparisons among continental groups]{\textbf{Pairwise mutation spectrum comparisons among continental groups.} Each of these plots compares the mutation spectra of two populations $P_1$ and $P_2$. Letting $f_i$ denote the fraction of SNVs in population $P_i$ that have a given triplet context, ancestral allele, and derived allele, the corresponding heat map square visualizes the enrichment ratio $f_1/f_2$. Black dots mark mutation types for which the difference between populations has a $\chi^2$ $p$-value less than $10^{-5}$. }\label{heatmap_all1kg}
%\end{figure}

%\begin{figure}

\begin{center}
%\includegraphics[width=0.49\textwidth]{../supplementary_material/diploid_mut_PCA_nosingle_wAMR.pdf}
%\includegraphics[width=0.24\textwidth]{../supplementary_material/mutation_PC0_nosingle_wAMR.pdf}
%\includegraphics[width=0.24\textwidth]{../supplementary_material/mutation_PC1_nosingle_wAMR.pdf}
\includegraphics[width=\textwidth]{../final-eLife-fig-supplements/Figure-1-supplement-2.pdf}
\end{center}

\vspace{-1.5cm}

\begin{singlespace}
\begin{flushleft}
{\bf Figure 1--Figure Supplement 2: PCA of all 1000 Genomes continental groups.} All admixed North and South American individuals were omitted from Figure 1 in the main text to clarify the separation of other populations along an African vs non-African axis and an East vs West Eurasian axis. Here, admixed Americans are added in black. As expected, some African-Americans group with the Africans, while other admixed Americans fall within the variation of other East and West Eurasians. The accompanying heat maps show the mutation type loadings of the first two principal components, the second of which is heavily weighted toward the European TCC$\to$TTC signature.
\end{flushleft}
\end{singlespace}


%\caption[PCA of all 1000 Genomes continental groups]{\textbf{PCA of all 1000 Genomes continental groups.} All admixed North and South American individuals were omitted from Figure 1 in the main text to clarify the separation of other populations along an African vs non-African axis and an East vs West Eurasian axis. Here, admixed Americans are added in black. As expected, some African-Americans group with the Africans, while other admixed Americans fall within the variation of other East and West Eurasians. The accompanying heat maps show the mutation type loadings of the first two principal components, the second of which is heavily weighted toward the European TCC$\to$TTC signature.}\label{PCA_wAMR}
%\end{figure}

\begin{center}
%\includegraphics[width=0.49\textwidth]{../supplementary_material/ranked_Pvals_EUR_AFR.pdf}
%\includegraphics[width=0.49\textwidth]{../supplementary_material/ranked_Pvals_EUR_AFR_zoomin.pdf}
\includegraphics[width=0.8\textwidth]{../final-eLife-fig-supplements/Figure-1-Supplement-3.pdf}
\end{center}

\vspace{-1cm}

\begin{flushleft}
\begin{singlespace}
{\bf Figure 1--Figure Supplement 3: Mutation spectrum comparison $p$-values.} Each left-hand plot shows all chi-squared $p$-values corresponding to the ratios from Figure 1A. In the absence of recent mutation spectrum evolution, only one out of 96 SNP categories is expected to have a $p$-value below 0.01 (lower dotted line). In contrast, the majority of $p$ values meet the more stringent threshold $p<1e-5$. The corresponding right hand panel shows a closeup of the distribution of $p$-values greater than $1e-5$.
\end{singlespace}
\end{flushleft}

%\caption[Europe vs Africa mutation spectrum $p$-values]{\textbf{Europe vs Africa mutation spectrum $p$-values. } The left-hand plot shows all chi-squared $p$-values corresponding to the ratios from Figure 1A. In the absence of recent mutation spectrum evolution, only one out of 96 SNP categories is expected to have a $p$-value below 0.01 (lower dotted line). In contrast, the majority of $p$ values meet the more stringent threshold $p<1e-5$. The right hand panel shows a closeup of the distribution of $p$-values greater than $1e-5$.}\label{pval_list}
%\end{figure}

%\begin{center}
%\includegraphics[width=0.49\textwidth]{../supplementary_material/ranked_Pvals_EUR_EAS.pdf}
%\includegraphics[width=0.49\textwidth]{../supplementary_material/ranked_Pvals_EUR_EAS_zoomin.pdf}
%\end{center}

%\vspace{-1cm}

%\begin{flushleft}
%\begin{singlespace}
%{\bf Figure 1--Figure Supplement 4: Europe vs East Asia mutation spectrum $p$-values.} 
%\end{singlespace}
%\end{flushleft}

%\begin{center}
%\includegraphics[width=0.49\textwidth]{../supplementary_material/ranked_Pvals_EAS_AFR.pdf}
%\includegraphics[width=0.49\textwidth]{../supplementary_material/ranked_Pvals_EAS_AFR_zoomin.pdf}
%\end{center}

%\begin{flushleft}
%\begin{singlespace}
%{\bf Figure 1--Figure Supplement 5: East Asia vs Africa mutation spectrum $p$-values.} 
%\end{singlespace}
%\end{flushleft}

%\begin{center}
%\includegraphics[width=0.49\textwidth]{../supplementary_material/ranked_Pvals_SAS_EAS.pdf}
%\includegraphics[width=0.49\textwidth]{../supplementary_material/ranked_Pvals_SAS_EAS_zoomin.pdf}
%\end{center}

%\vspace{-1cm}

%\begin{flushleft}
%\begin{singlespace}
%{\bf Figure 1--Figure Supplement 6: South Asia vs East Asia mutation spectrum $p$-values.} 
%\end{singlespace}
%\end{flushleft}

\flushleft

\includegraphics[width=\textwidth]{../supplementary_material/BGC_effects_violin.pdf}

\begin{singlespace}
\noindent {\bf Figure 1--Figure Supplement 4: The effects of biased gene conversion on mutation spectra.} When using segregating variation to study the mutation spectrum, one potential source of bias is that strong-to-weak mutations, where the ancestral allele is G or C and the derived allele is A or T, have a lower fixation probability than weak-to-strong mutations due to biased gene conversion (BGC). If this effect were sufficiently strong, it would inflate the apparent mutation fractions of weak-to-strong mutations, especially in populations with large effective sizes where natural selection is particularly efficient. Within humans, Africans have the largest long-term effective population size, while East Asians and Native Americans have the lowest. Therefore, if BGC has created differences in mutation spectra between populations, the fraction of weak-to-strong SNVs should be highest in Africans, intermediate in Europeans and South Asians, and lowest in East Asians and Native Americans. This violin plot reveals no such pattern, suggesting that BGC is not a strong driver of mutation spectrum differences between human populations. We do not observe either a direct correlation between in strong-to-weak mutation fraction and distance from Africa or an inverse correlation between weak-to-strong mutation fraction and distance from Africa.
\end{singlespace}

%\begin{centering}
%\begin{tabular}{cc}
%A. (+)-Transcribed DNA & B. (-)-Transcribed DNA \\
%\includegraphics[width=0.49\textwidth]{../supplementary_material/heatmap_strand+.pdf} & \includegraphics[width=0.49\textwidth]{../supplementary_material/heatmap_strand-.pdf} \\
%C. Nontranscribed DNA & Enrichment on (+)-strand relative to (-)-strand \\
%\includegraphics[width=0.49\textwidth]{../supplementary_material/heatmap_nontranscribed.pdf} & \includegraphics[width=0.49\textwidth]{../supplementary_material/strandbias_heatmap.pdf}
%\end{tabular}
%\end{centering}

\begin{center}
\includegraphics[width=\textwidth]{../final-eLife-fig-supplements/Figure-1-Supplement-5.pdf}
\end{center}

\begin{flushleft}
\begin{singlespace}
{\bf Figure 1--Figure Supplement 5: Mutation Spectra of Transcribed vs Non-Transcribed DNA.} Using the UCSC Genome Browser annotations of the human reference hg19, we determined whether each SNP occurs in a transcribed or non transcribed region. We further divided SNPs occurring in transcribed regions according to whether the ancestral A or C allele occurs on the (+)-strand or the (-)-strand. Panels A, B, and C all show the same population-specific mutation type enrichments that are observed in Figure 1B. Panel D plots the residuals between panels A and B, highlighting mutation types that show a modest difference in strand bias between populations. 
\end{singlespace}
\end{flushleft}

%\begin{centering}
%\begin{tabular}{cc}
%1. Active Promoter & 2. Weak Promoter\\
%\includegraphics[width=0.49\textwidth]{../supplementary_material/heatmap_1_Active_Promoter_p1e2_nosingletons.pdf} & \includegraphics[width=0.49\textwidth]{../supplementary_material/heatmap_2_Weak_Promoter_p1e2_nosingletons.pdf} \\
%3. Poised Promoter & 4. Strong Enhancer\\
%\includegraphics[width=0.49\textwidth]{../supplementary_material/heatmap_3_Poised_Promoter_p1e2_nosingletons.pdf} & \includegraphics[width=0.49\textwidth]{../supplementary_material/heatmap_4_Strong_Enhancer_p1e2_nosingletons.pdf} \\
%5. Strong Enhancer & 6. Weak Enhancer \\
%\includegraphics[width=0.49\textwidth]{../supplementary_material/heatmap_5_Strong_Enhancer_p1e2_nosingletons.pdf} & \includegraphics[width=0.49\textwidth]{../supplementary_material/heatmap_6_Weak_Enhancer_p1e2_nosingletons.pdf} \\
%\end{tabular}
%\end{centering}

\begin{center}
\includegraphics[width=\textwidth]{../final-eLife-fig-supplements/Figure-1-Supplement-6.pdf}
\end{center}

\begin{flushleft}
\begin{singlespace}
{\bf Figure 1--Figure Supplement 6: Mutation Spectra of ChromHMM chromatin states (Part I of II).} To investigate whether any mutation spectrum shifts might be confined to particular chromatin states, we used chromHMM annotations of the human embryonic stem cell line HESC-H1 \cite{hoffman2013}. Each heat map plots mutation spectrum comparisons for SNPs that are annotated as being part of the same chromatin state, and dots mark mutation types that show a significant enrichment in one population at the level $p<0.01$. Every chromatin state shows enrichment of the TCC$\to$TTC signature in Europe and South Asia. Some heat maps are noisy due to the small sample size of SNPs contained within these regions, but all showcase the same general patterns as Figure 1B.
\end{singlespace}
\end{flushleft}


%\begin{centering}
%\begin{tabular}{cc}
%7. Weak Enhancer & 8. Insulator \\
%\includegraphics[width=0.49\textwidth]{../supplementary_material/heatmap_7_Weak_Enhancer_p1e2_nosingletons.pdf} & \includegraphics[width=0.49\textwidth]{../supplementary_material/heatmap_8_Insulator_p1e2_nosingletons.pdf} \\
%9. Transcription Transition & 10. Transcription Elongation \\
%\includegraphics[width=0.49\textwidth]{../supplementary_material/heatmap_9_Txn_Transition_p1e2_nosingletons.pdf} & \includegraphics[width=0.49\textwidth]{../supplementary_material/heatmap_10_Txn_Elongation_p1e2_nosingletons.pdf} \\
%11. Weak Transcription & 12. Repressed \\
%\includegraphics[width=0.49\textwidth]{../supplementary_material/heatmap_11_Weak_Txn_p1e2_nosingletons.pdf} & \includegraphics[width=0.49\textwidth]{../supplementary_material/heatmap_12_Repressed_p1e2_nosingletons.pdf} \\
%\end{tabular}
%\end{centering}

\begin{center}
\includegraphics[width=\textwidth]{../final-eLife-fig-supplements/Figure-1-Supplement-7.pdf}
\end{center}

\begin{flushleft}
\begin{singlespace}
{\bf Figure 1--Figure Supplement 7: Mutation Spectra of ChromHMM chromatin states (Part II of II).} 
\end{singlespace}
\end{flushleft}

%\centering
%\includegraphics[width=\textwidth]{../supplementary_material/heatmap_EUR_AFR_by_reptiming.pdf}
%\includegraphics[width=\textwidth]{../supplementary_material/heatmap_EAS_AFR_by_reptiming.pdf}

\begin{center}
\includegraphics[width=\textwidth]{../final-eLife-fig-supplements/Figure-1-Supplement-8.pdf}
\end{center}

\begin{flushleft}
\begin{singlespace}
{\bf Figure 1--Figure Supplement 8: Variation of the mutation spectrum with DNA replication timing.} We partitioned the genome into 10 equal replication timing quantiles using data obtained from \cite{woodfine2004}, then computed mutation spectrum differences within each quantile. Although most patterns from Figure 1B replicate within each replication timing bin, there are a few exceptions. CpG transitions, which occur most often in early-replicating regions, vary in population bias depending on replication timing. In addition, the deficit of ACA$\to$AAA and AAA$\to$ATA mutations in Africa compared to Europe and Asia is observed mainly in early-replicating regions. 
\end{singlespace}
\end{flushleft}

\begin{singlespace}
\noindent {\bf Figure 1--Source Data 1.} This text file shows the number of SNPs in each of the 96 mutational categories that passed all filters in each 1000 Genomes continental group. 
\end{singlespace}


%\FloatBarrier

%\section{Replication of evidence for mutation spectrum evolution in the SGDP}
%To verify that our putative signatures of mutation spectrum change are not artifacts that were someone introduced through the 1000 Genomes pipeline, we replicated our human heat map results using the Simons Genomic Diversity Panel (SGDP), a dataset produced independently of the 1000 Genomes panel using different cell lines representing different fine-scale populations. Because the SGDP samples fine-scale populations broadly instead of deeply, with only 2-3 individuals drawn from each population, this dataset was not suited for replicating our results on fine-scale population structure. At the continental scale, however, we observed divergence between the mutation spectra of the SGDP populations that closely matched our 1000 Genomes results. Qualitative agreement between the heat maps is seen in Figure S\ref{simons_heatmap}, and a quantitative regression of SGDP enrichments versus 1000 Genomes enrichments is shown in Figure S\ref{regression}.

%\begin{centering}
%\begin{tabular}{cc}
%A. Simons Genomic Diversity Panel (SGDP) & B. 1000 Genomes Phase 3 Panel (1000G) \\
%\includegraphics[width=0.49\textwidth]{../supplementary_material/heatmap_simons_filtered_nosingletons.pdf} & \includegraphics[width=0.49\textwidth]{../supplementary_material/heatmap_1kg_projected_nosingletons.pdf} \\
%\end{tabular}
%\end{centering}

\begin{center}
\includegraphics[width=\textwidth]{../final-eLife-fig-supplements/Figure-2-Supplement-1.pdf}
\end{center}

\begin{flushleft}
\begin{singlespace}
{\bf Figure 2--Figure Supplement 1: Heatmap comparisons between continental groups in 1000 Genomes and the SGDP.} Here, each 1000 Genomes population is projected down to the sample size of the corresponding SGDP population in order to sample alleles with a similar distribution of ages and frequencies.
\end{singlespace}
\end{flushleft}

%\begin{figure}
%\includegraphics[width=0.49\textwidth]{../supplementary_material/simons_1kg_ratios_WestEurasia_v_Africa_nosingle.pdf}
%\includegraphics[width=0.49\textwidth]{../supplementary_material/simons_1kg_ratios_EastAsia_v_Africa_nosingle.pdf} 
%\includegraphics[width=0.49\textwidth]{../supplementary_material/simons_1kg_ratios_EastAsia_v_WestEurasia_nosingle.pdf}  
%\includegraphics[width=0.49\textwidth]{../supplementary_material/simons_1kg_ratios_EastAsia_v_SouthAsia_nosingle.pdf}  

%\includegraphics[width=0.49\textwidth]{../supplementary_material/projected_simons_1kg_ratios_WestEurasia_v_SouthAsia_nosingle.pdf}  
%\includegraphics[width=0.49\textwidth]{../supplementary_material/projected_simons_1kg_ratios_SouthAsia_v_Africa_nosingle.pdf}  

\begin{center}
\includegraphics[width=\textwidth]{../final-eLife-fig-supplements/Figure-2-Supplement-2.pdf}
\end{center}

\begin{flushleft}
\begin{singlespace}
{\bf Figure 2--Figure Supplement 2: Regression of the SGDP heatmap coefficients versus the corresponding 1000 Genomes heatmap coefficients.} 
\end{singlespace}
\end{flushleft}


%\FloatBarrier

%\section{The mutation spectrum as a function of allele frequency}

%Figures 2B and 2C of the main text show how four mutation types: TCC$\to$TTC, ACC$\to$ATC, TCT$\to$TTT, and CCC$\to$CTC, have pulse-like distributions in Europe and South Asia. Their abundance is not monotonic as a function of allele frequency, but peaks around 1\% in the 1000 Genomes data and around 0.6\% in the larger UK10K dataset. This section puts these results into context by providing a more complete description of how the mutation spectrum varies with allele frequency in these datasets. 

%Figure S\ref{TCC_1kg} shows how TCC$\to$TTC abundance has a pulse-like distribution within all European, South Asian, and American populations,. Among American populations, we can see that the most European-admixed groups (Mexicans, Colombians, and Puerto Ricans) have higher-frequency TCC$\to$TTC mutations than the less admixed Peruvians, African-Americans, and Afro-Caribbeans. TCC$\to$TTC mutations in South Asians are also lower, frequency, on average than those in  European or American populations. Within Europe, the Finns  have fewer high-frequency TCC$\to$TTC mutations than other Europeans do, but a similar quantity of low-frequency TCC$\to$TTC mutations.

%Figures S\ref{UK10K_type_v_freq1} and S\ref{UK10K_type_v_freq2} show a comprehensive picture of how the mutation spectrum varies with allele frequency in the UK10K. Every mutation type $m$ is represented by a line that is expected to fluctuate around $y=1$ unless the rate of $m$ has recently increased or decreased. All 96 mutation types are visualized in every panel, but most corresponding lines are greyed out to enhance readability. Only a few related mutation types are highlighted in color. 



\begin{center}
\includegraphics[width=0.9\textwidth]{../supplementary_material/finescale_fTCC_v_allele_freq_all1kg_UK10K.pdf}
\end{center}

\begin{flushleft}
\begin{singlespace}
{\bf Figure 3--Figure Supplement 1: TCC$\to$TTC mutation fraction as a function of allele frequency in all 1000 Genomes populations.} To enable better comparison with the 1000 Genomes data, the UK10K SNPs have been downsampled to 200 individuals. The age distribution of alleles of a given frequency varies as a function of the number of lineages being sampled--this is why the UK10K pulse peaks around 0.6\% frequency when measured in a dataset of thousands of lineages, but peaks around 2\% in a subsample of only 400 lineages. Some African and East Asian population names have been omitted for clarity since the TCC$\to$TTC mutation fraction is so uniform within these continental groups. Red = European populations; Blue = South Asian; Orange = Americas; Purple = Africa; Green = East Asia.
\end{singlespace}
\end{flushleft}

\centering
\includegraphics[width=0.7\textwidth]{../supplementary_material/fTCC_v_allele_freq_1kg_ExAC_evenbins.pdf}

\begin{flushleft}
\begin{singlespace}
{\bf Figure 3--Figure Supplement 2: Fraction of TCC$\to$TTC mutations as a function of allele frequency in ExAC.} Lek, et al. compiled data from 60,706 exomes to create the Exome Aggregation Consortium dataset, which enables the analysis of ultra-rare human variation \cite{lek2016}. The overall fraction of TCC$\to$TTC mutations is slightly higher in exome data than in whole genome data because exons contain a skewed distribution of triplet contexts, but the pulse pattern from Figure 3B reproduces unmistakably. 
\end{singlespace}
\end{flushleft}

\newpage 

%\includegraphics[width=0.49\textwidth]{../supplementary_material/all_mut_fracs_EURpulse_CT.pdf}
%\includegraphics[width=0.49\textwidth]{../supplementary_material/all_mut_fracs_EURpulse_CpGts.pdf}

\begin{center}
\includegraphics[width=\textwidth]{../final-eLife-fig-supplements/Figure-3-Supplement-3.pdf}
\end{center}

\begin{flushleft}
\begin{singlespace}
{\bf Figure 3--Figure Supplement 3: Mutation type enrichment as a function of allele frequency in UK10K (Part I of III).} The eleven panels in Figure Supplements 2, 3, and 4 show the full dependence of mutation spectrum on allele frequency in the UK10K data. If we let $F(f,m)$ denote the fraction of SNVs of frequency $f$ that are of type $m$ and let $F(m)$ denote the fraction of all mutations that are of type $m$, the enrichment of mutation type $m$ as a function of frequency is $F(f,m)/F(m)$. This function is expected to fluctuate around $y=1$ unless the rate of $m$ has recently increased or decreased. All 96 mutation types are visualized in every panel, but most corresponding lines are greyed out to enhance readability. Some lines deviate from $y=1$ due to the effects of biased gene conversion (BGC)--this occurs when one of the ancestral or derived alleles is a weak base (A or T, abbreviated W) and the other allele is a strong base (G or C, abbreviated S). W$\to$S mutations are more abundant at high allele frequencies, while S$\to$W mutations are more abundant at low frequencies. These effects are visible but modest in panels D, G, H, and I, but much more pronounced in panels B, C, and F, which focus on mutations in the CpG context. Transitions of the type CpA$\to$CpG, which create CpG motifs, are extremely enriched at high frequencies, and this pattern may be an artifact of ancestral misidentification \cite{hernandez2007}. CpG motifs have such high mutation rates that CpG$\to$CpT transitions often happen at the same site in humans and chimps, and these low-frequency double mutations are misclassified as high-frequency CpT$\to$CpG mutations. Although it is not surprising to see a peak of CpT$\to$CpG transitions at high frequencies in panel F, it is somewhat surprising to see CpG$\to$GpG transversions peak in abundance at high frequencies in panel C. This might be a signature of recent declines in the rates of these mutations, since neither ancestral misidentification nor biased gene conversion is thought to produce such a pattern. In addition, neither of these processes can explain the strong enrichment of certain A$\to$T mutations at high frequencies that is observed in panel K.
\end{singlespace}
\end{flushleft}


%\includegraphics[width=0.49\textwidth]{../supplementary_material/all_mut_fracs_EURpulse_CpGtv.pdf}
%\includegraphics[width=0.49\textwidth]{../supplementary_material/all_mut_fracs_EURpulse_CA.pdf}

%
%\includegraphics[width=0.49\textwidth]{../supplementary_material/all_mut_fracs_EURpulse_CG.pdf}
%\includegraphics[width=0.49\textwidth]{../supplementary_material/all_mut_fracs_EURpulse_invCpG.pdf}

%\includegraphics[width=0.49\textwidth]{../supplementary_material/all_mut_fracs_EURpulse_AG.pdf}

\begin{center}
\includegraphics[width=\textwidth]{../final-eLife-fig-supplements/Figure-3-Supplement-4.pdf}
\end{center}

\begin{flushleft}
\begin{singlespace}
{\bf Figure 3--Figure Supplement 4: Mutation type enrichment as a function of allele frequency in UK10K (Part II of III).} The eleven panels in this 3-part figure show the full dependence of mutation spectrum on allele frequency in the UK10K data.
\end{singlespace}
\end{flushleft}

\begin{center}
\includegraphics[width=\textwidth]{../final-eLife-fig-supplements/Figure-3-Supplement-5.pdf}
\end{center}

%\includegraphics[width=0.49\textwidth]{../supplementary_material/all_mut_fracs_EURpulse_AC1.pdf}
%\includegraphics[width=0.49\textwidth]{../supplementary_material/all_mut_fracs_EURpulse_AC2.pdf}

%\includegraphics[width=0.49\textwidth]{../supplementary_material/all_mut_fracs_EURpulse_AT1.pdf}
%\includegraphics[width=0.49\textwidth]{../supplementary_material/all_mut_fracs_EURpulse_AT2.pdf}

\begin{flushleft}
\begin{singlespace}
{\bf Figure 3--Figure Supplement 5: Mutation type enrichment as a function of allele frequency in UK10K (Part III of III).} The eleven panels in this 3-part figure show the full dependence of mutation spectrum on allele frequency in the UK10K data.
\end{singlespace}
\end{flushleft}

\newpage


%\includegraphics[width=0.49\textwidth]{../supplementary_material/mut_spectrum_by_allele_freq_EUR_v_AFR.pdf}
%\includegraphics[width=0.49\textwidth]{../supplementary_material/mut_spectrum_by_allele_freq_EAS_v_AFR.pdf}

%\includegraphics[width=0.49\textwidth]{../supplementary_material/mut_spectrum_by_allele_freq_EUR_v_EAS.pdf}
%\includegraphics[width=0.49\textwidth]{../supplementary_material/mut_spectrum_by_allele_freq_SAS_v_EAS.pdf}

\begin{center}
\includegraphics[width=\textwidth]{../final-eLife-fig-supplements/Figure-3-Supplement-6.pdf}
\end{center}

\begin{flushleft}
\begin{singlespace}
{\bf Figure 3--Figure Supplement 6: Mutation spectrum comparisons partitioned by allele frequency.} Each of these heatmaps shows a subset of the data used to construct Figure 1B, partitioned by allele frequency to show how rare variants are the most highly differentiated between populations. Black dots highlight mutation types that are significantly different in abundance between two populations in a particular frequency class at the $p<10^{-5}$ level according to a chi-square test. 
\end{singlespace}
\end{flushleft}

%\FloatBarrier

%\section{Finescale mutation spectra of 1kG continental groups}

%Figure 3B of the main text shows heat map comparisons between East Asian populations, which display fine-scale differences that are exceptionally well defined. For completeness, this section shows finescale heatmap comparisons within all 1kG continental groups. We can see that CAC$\to$CCC and TAT$\to$TTT are heterogeneously distributed within multiple continents, but to the greatest extent in East Asia. In addition the TCC$\to$TTC signature is somewhat heterogeneously distributed within Europe and South Asia, being depleted in Finns  and enriched in the Punjabi and Gujarati.


%\begin{tabular}{cc}
%A. Finescale African Mutation Spectra &  B. Finescale European Mutation Spectra \\
%\centering
%\includegraphics[width=0.45\textwidth]{../supplementary_material/AFR_heatmap_v_MSL_nosingle.pdf} &
%\includegraphics[width=0.45\textwidth]{../supplementary_material/EUR_heatmap_v_IBS_nosingle.pdf} \\
%C. Finescale East Asian Mutation Spectra &  D. Finescale South Asian Mutation Spectra \\
%\includegraphics[width=0.45\textwidth]{../supplementary_material/ASN_heatmap_v_CHS_nosingle.pdf} &
%\includegraphics[width=0.45\textwidth]{../supplementary_material/SAS_heatmap_v_ITU_nosingle.pdf} 
%\end{tabular}

\begin{center}
\includegraphics[width=\textwidth]{../final-eLife-fig-supplements/Figure-4-Supplement-1.pdf}
\end{center}

\vspace{-1cm}

\begin{flushleft}
\begin{singlespace}
{\bf Figure 4--Figure Supplement 1: Mutation spectrum differences within Africa, Europe, East Asia, and South Asia.} Figure 4B of the main text shows heat map comparisons between East Asian populations, which display fine-scale differences that are exceptionally well defined. For completeness, this figure shows finescale heatmap comparisons within all 1kG continental groups.  We can see that CAC$\to$CCC and TAT$\to$TTT are heterogeneously distributed within multiple continents, but to the greatest extent in East Asia. In addition the TCC$\to$TTC signature is somewhat heterogeneously distributed within Europe and South Asia, being depleted in Finns  and enriched in the Punjabi and Gujarati.
 Each continental group in the 1000 Genomes data is divided into 5 sub-populations. These heat maps compare the mutation spectra of these fine-scale populations to each other.  African populations are: MSL = Mende in Sierra Leone; LWK = Luhya in Webuye, Kenya; YRI = Yoruba in Ibadan, Nigeria; GWD = Gambian in Western Divisions; ESN = Esan in Nigeria. European populations are: IBS = Iberian Population in Spain; TSI = Toscani in Italia; GBR = British in England and Scotland; CEU = Utah Residents (CEPH) with Northern and Western Ancestry; FIN = Finnish in Finland. East Asian populations are: CDX = Chinese Dai in Xishuangbanna, China; JPT = Japanese in Tokyo, Japan; CHB = Han Chinese in Bejing, China; CHS = Southern Han Chinese; KHV = Kinh in Ho Chi Minh City, Vietnam. South Asian populations are: ITU = Indian Telugu from the UK; GIH = Gujarati Indian from Houston, Texas; PJL = Punjabi from Lahore, Pakistan; BEB = Bengali from Bangladesh; STU = Sri Lankan Tamil from the UK.
 \end{singlespace}
 \end{flushleft}
%\FloatBarrier

%\section{Finescale principal component analysis of mutation spectra (1kG)}

%In this section, we show a principal component analysis of individual mutation spectra within each 1000 Genomes continental group, complementing the plot shown in Figure 3A of the main text. Some groups such as South Asians and Europeans show very little substructure, while East Asians and Africans show significant substructure related to the AC$\to$CC signature.  Each PCA plot is accompanied by two heatmaps that visualize the loadings of the first two PCs.



%\includegraphics[width=0.49\textwidth]{../supplementary_material/AMR_mut_PCA_1kg_nosingle_altlegend.pdf}
%\includegraphics[width=0.24\textwidth]{../supplementary_material/AMR_mut_PC0_nosingle.pdf}
%\includegraphics[width=0.24\textwidth]{../supplementary_material/AMR_mut_PC1_nosingle.pdf}

\begin{center}
\includegraphics[width=\textwidth]{../final-eLife-fig-supplements/Figure-4-Supplement-2.pdf}
\end{center}

\begin{singlespace}
\noindent {\bf Figure 4--Figure Supplement 2: PCA of American populations.} Population abbreviations are: CLM = Colombians from Medellin, Colombia; MXL = Mexican Ancestry from Los Angeles, USA; PUR = Puerto Ricans from Puerto Rico; PEL = Peruvians from Lima, Peru; ACB = African Caribbeans in Barbados; ASW = Americans of African Ancestry in SW USA. Admixed populations from the Americans show structure that mirrors the continental groups, with PC1 essentially measuring the ratio between African and non-African ancestry and PC2 measuring the ratio between European and Native American ancestry. The accompanying heat maps show the loadings of the first two principal components.
\end{singlespace}

\begin{center}
\includegraphics[width=\textwidth]{../final-eLife-fig-supplements/Figure-4-Supplement-3.pdf}
\end{center}

%\includegraphics[width=0.49\textwidth]{../supplementary_material/AFR_mut_PCA_1kg_nosingle_altlegend.pdf}
%\includegraphics[width=0.24\textwidth]{../supplementary_material/AFR_mut_PC0_nosingle.pdf}
%\includegraphics[width=0.24\textwidth]{../supplementary_material/AFR_mut_PC1_nosingle.pdf}

\begin{singlespace}
\noindent {\bf Figure 4--Figure Supplement 3: PCA of African populations.} Population abbreviations are: MSL = Mende in Sierra Leone; LWK = Luhya in Webuye, Kenya; YRI = Yoruba in Ibadan, Nigeria; GWD = Gambian in Western Divisions; ESN = Esan in Nigeria. 
\end{singlespace}

\begin{center}
\includegraphics[width=\textwidth]{../final-eLife-fig-supplements/Figure-4-Supplement-4.pdf}
\end{center}

%\includegraphics[width=0.49\textwidth]{../supplementary_material/EAS_mut_PCA_1kg_nosingle_altlegend.pdf}
%\includegraphics[width=0.24\textwidth]{../supplementary_material/EAS_mut_PC0_nosingle.pdf}
%\includegraphics[width=0.24\textwidth]{../supplementary_material/EAS_mut_PC1_nosingle.pdf}

\begin{singlespace}
\noindent {\bf Figure 4--Figure Supplement 4: PCA of East Asian populations.} Population abbreviations are: CDX = Chinese Dai in Xishuangbanna, China; JPT = Japanese in Tokyo, Japan; CHB = Han Chinese in Bejing, China; CHS = Southern Han Chinese; KHV = Kinh in Ho Chi Minh City, Vietnam. 
\end{singlespace}

\begin{center}
\includegraphics[width=\textwidth]{../final-eLife-fig-supplements/Figure-4-Supplement-5.pdf}
\end{center}

%\includegraphics[width=0.49\textwidth]{../supplementary_material/SAS_mut_PCA_1kg_nosingle_altlegend.pdf}
%\includegraphics[width=0.24\textwidth]{../supplementary_material/SAS_mut_PC0_nosingle.pdf}
%\includegraphics[width=0.24\textwidth]{../supplementary_material/SAS_mut_PC1_nosingle.pdf}

\begin{singlespace}
\noindent {\bf Figure 4--Figure Supplement 5: PCA of South Asian populations.} Population abbreviations are:  ITU = Indian Telugu from the UK; GIH = Gujarati Indian from Houston, Texas; PJL = Punjabi from Lahore, Pakistan; BEB = Bengali from Bangladesh; STU = Sri Lankan Tamil from the UK.
\end{singlespace}

\begin{center}
\includegraphics[width=\textwidth]{../final-eLife-fig-supplements/Figure-4-Supplement-6.pdf}
\end{center}

%\includegraphics[width=0.49\textwidth]{../supplementary_material/EUR_mut_PCA_1kg_nosingle_altlegend.pdf}
%\includegraphics[width=0.24\textwidth]{../supplementary_material/EUR_mut_PC0_nosingle.pdf}
%\includegraphics[width=0.24\textwidth]{../supplementary_material/EUR_mut_PC1_nosingle.pdf}

\begin{singlespace}
\noindent {\bf Figure 4--Figure Supplement 6: PCA of European populations.} Population abbreviations are:  IBS = Iberian Population in Spain; TSI = Toscani in Italia; GBR = British in England and Scotland; CEU = Utah Residents (CEPH) with Northern and Western Ancestry; FIN = Finnish in Finland.
\end{singlespace}


\begin{singlespace}
\noindent {\bf Figure 4--Source Data 1.} This text file shows the number of SNPs in each of the 96 mutational categories that passed all filters in each finescale 1000 Genomes population. 
\end{singlespace}

%\FloatBarrier

%\section{Mutation spectra from the Great Ape Genome Diversity Project }

%Figure S\ref{ape_heatmaps} fleshes out the patterns that are visible in Figure 4, which compares all other great apes to the Bornean Orangutan. As expected, we can see that Bonobos and Chimps have more correlated mutation spectra than more distantly related pairs of ape species.
\newpage

\begin{center}
\includegraphics[width=\textwidth]{../final-eLife-fig-supplements/Figure-5-Supplement-1.pdf}
\end{center}

%\includegraphics[width=0.49\textwidth]{../supplementary_material/ape_heatmap_v_Homo_nosingle.pdf}
%\includegraphics[width=0.49\textwidth]{../supplementary_material/ape_heatmap_v_Gorilla_nosingle.pdf}

%\includegraphics[width=0.49\textwidth]{../supplementary_material/ape_heatmap_v_Pan_troglodytes_nosingle.pdf}
%\includegraphics[width=0.49\textwidth]{../supplementary_material/ape_heatmap_v_Pan_paniscus_nosingle.pdf}

%\includegraphics[width=0.49\textwidth]{../supplementary_material/ape_heatmap_v_Pongo_abelii_nosingle.pdf}
%\includegraphics[width=0.49\textwidth]{../supplementary_material/ape_heatmap_v_Pongo_pygmaeus_nosingle.pdf}

\begin{singlespace}
\noindent {\bf Figure 5--Figure Supplement 1: Mutation spectra of great apes.} These heatmap comparisons demonstrate that closely related great apes such as Chimpanzees and Bonobos have more similar mutation spectra than more distantly related apes do.
\end{singlespace}


\FloatBarrier

\newpage

\bibliography{scibib}

%  {scilastnote} is a Science formatted environment to put Acknowledgements as last item in the bibliography
%  to implement this remove the comments and \paragraph{Ack}  -JP



% \end{scilastnote}

% \item[Correspondence] Correspondence and requests for materials
%should be addressed to K.H.~(email: kelleyh@stanford.edu) and J.K.P.~(email:pritch@stanford.edu).


\end{document}

