Peter R. Wilton

Leveraging Distant Relatedness to Quantify Human Mutation and Gene-Conversion Rates

The rate at which human genomes mutate is a central biological parameter that has many implications for our ability to understand demographic and evolutionary phenomena. We present a method for inferring mutation and gene-conversion rates by using the number of sequence differences observed in identical-by-descent (IBD) segments together with a reconstructed model of recent population-size history. This approach is robust to, and can quantify, the presence of substantial genotyping error, as validated in coalescent simulations. We applied the method to 498 trio-phased sequenced Dutch individuals and inferred a point mutation rate of 1.66 × 10(-8) per base per generation and a rate of 1.26 × 10(-9) for <20 bp indels. By quantifying how estimates varied as a function of allele frequency, we inferred the probability that a site is involved in non-crossover gene conversion as 5.99 × 10(-6). We found that recombination does not have observable mutagenic effects after gene conversion is accounted for and that local gene-conversion rates reflect recombination rates. We detected a strong enrichment of recent deleterious variation among mismatching variants found within IBD regions and observed summary statistics of local sharing of IBD segments to closely match previously proposed metrics of background selection; however, we found no significant effects of selection on our mutation-rate estimates. We detected no evidence of strong variation of mutation rates in a number of genomic annotations obtained from several recent studies. Our analysis suggests that a mutation-rate estimate higher than that reported by recent pedigree-based studies should be adopted in the context of DNA-based demographic reconstruction.

Effects of the population pedigree on genetic signatures of historical demographic events

Genetic variation among loci in the genomes of diploid biparental organisms is the result of mutation and genetic transmission through the genealogy, or population pedigree, of the species. We explore the consequences of this for patterns of variation at unlinked loci for two kinds of demographic events: the occurrence of a very large family or a strong selective sweep that occurred in the recent past. The results indicate that only rather extreme versions of such events can be expected to structure population pedigrees in such a way that unlinked loci will show deviations from the standard predictions of population genetics, which average over population pedigrees. The results also suggest that large samples of individuals and loci increase the chance of picking up signatures of these events, and that very large families may have a unique signature in terms of sample distributions of mutant alleles.

The SMC′ Is a Highly Accurate Approximation to the Ancestral Recombination Graph

Two sequentially Markov coalescent models (SMC and SMC′) are available as tractable approximations to the ancestral recombination graph (ARG). We present a Markov process describing coalescence at two fixed points along a pair of sequences evolving under the SMC′. Using our Markov process, we derive a number of new quantities related to the pairwise SMC′, thereby analytically quantifying for the first time the similarity between the SMC′ and the ARG. We use our process to show that the joint distribution of pairwise coalescence times at recombination sites under the SMC′ is the same as it is marginally under the ARG, which demonstrates that the SMC′ is, in a particular well-defined, intuitive sense, the most appropriate first-order sequentially Markov approximation to the ARG. Finally, we use these results to show that population size estimates under the pairwise SMC are asymptotically biased, while under the pairwise SMC′ they are approximately asymptotically unbiased.

A Population Phylogenetic View of Mitochondrial Heteroplasmy

The mitochondria contained within the human body are genetically diverse. This type of variation, called heteroplasmy, is emerging as an important factor... The mitochondrion has recently emerged as an active player in myriad cellular processes. Additionally, it was recently shown that >200 diseases are known to be linked to variants in mitochondrial DNA or in nuclear genes interacting with mitochondria. This has reinvigorated interest in its biology and population genetics. Mitochondrial heteroplasmy, or genotypic variation of mitochondria within an individual, is now understood to be common in humans and important in human health. However, it is still not possible to make quantitative predictions about the inheritance of heteroplasmy and its proliferation within the body, partly due to the lack of an appropriate model. Here, we present a population-genetic framework for modeling mitochondrial heteroplasmy as a process that occurs on an ontogenetic phylogeny, with genetic drift and mutation changing heteroplasmy frequencies during the various developmental processes represented in the phylogeny. Using this framework, we develop a Bayesian inference method for inferring rates of mitochondrial genetic drift and mutation at different stages of human life. Applying the method to previously published heteroplasmy frequency data, we demonstrate a severe effective germline bottleneck comprised of the cumulative genetic drift occurring between the divergence of germline and somatic cells in the mother, and the separation of germ layers in the offspring. Additionally, we find that the two somatic tissues we analyze here undergo tissue-specific bottlenecks during embryogenesis, less severe than the effective germline bottleneck, and that these somatic tissues experience little additional genetic drift during adulthood. We conclude with a discussion of possible extensions of the ontogenetic phylogeny framework and its possible applications to other ontogenetic processes in addition to mitochondrial heteroplasmy.

Population structure and coalescence in pedigrees: Comparisons to the structured coalescent and a framework for inference

Contrary to what is often assumed in population genetics, independently segregating loci do not have completely independent ancestries, since all loci are inherited through a single, shared population pedigree. Previous work has shown that the non-independence between gene genealogies of independently segregating loci created by the population pedigree is weak in panmictic populations, and predictions made from standard coalescent theory are accurate for populations that are at least moderately sized. Here, we investigate patterns of coalescence in pedigrees of structured populations. We find that the pedigree creates deviations away from the predictions of the structured coalescent that persist on a longer timescale than in the case of panmictic populations. Nevertheless, we find that the structured coalescent provides a reasonable approximation for the coalescent process in structured population pedigrees so long as migration events are moderately frequent and there are no migration events in the recent pedigree of the sample. When there are migration events in the recent sample pedigree, we find that distributions of coalescence in the sample can be modeled as a mixture of distributions from different initial sample configurations. We use this observation to motivate a maximum-likelihood approach for inferring migration rates and mutation rates jointly with features of the pedigree such as recent migrant ancestry and recent relatedness. Using simulation, we show that our inference framework accurately recovers long-term migration rates in the presence of recent migration events in the sample pedigree.

The effects of Pleistocene climate change on biotic differentiation in a montane songbird clade from Wallacea

The role of Pleistocene Ice Age in tropical diversification is poorly understood, especially in archipelagos, in which glaciation-induced sea level fluctuations may lead to complicated changes in land distribution. To assess how Pleistocene land bridges may have facilitated gene flow in tropical archipelagos, we investigated patterns of diversification in the rarely-collected rusty-bellied fantail Rhipidura teysmanni (Passeriformes: Rhipiduridae) complex from Wallacea using a combination of bioacoustic traits and whole-genome sequencing methods (dd-RADSeq). We report a biogeographic leapfrog pattern in the vocalizations of these birds, and uncover deep genomic divergence among island populations despite the presence of intermittent land connections between some. We demonstrate how rare instances of genetic introgression have affected the evolution of this species complex, and document the presence of double introgressive mitochondrial sweeps, highlighting the dangers of using only mitochondrial DNA in evolutionary research. By applying different tree inference approaches, we demonstrate how concatenation methods can give inaccurate results when investigating divergence in closely-related taxa. Our study highlights high levels of cryptic avian diversity in poorly-explored Wallacea, elucidates complex patterns of Pleistocene climate-mediated diversification in an elusive montane songbird, and suggests that Pleistocene land bridges may have accounted for limited connectivity among montane Wallacean biota.

Pleistocene land bridges act as semipermeable agents of avian gene flow in Wallacea

Cyclical periods of global cooling have been important drivers of biotic differentiation throughout the Quaternary. Ice age-induced sea level fluctuations can lead to changing patterns of land connections, both facilitating and disrupting gene flow. In this study, we test if species with differing life histories are differentially affected by Quaternary land connections. We used genome-wide SNPs in combination with mitochondrial gene sequences to analyse levels of divergence and gene flow between two songbird complexes across two Wallacean islands that have been repeatedly connected during glaciations. Although the two bird complexes are similar in ecological attributes, the forest and edge-inhabiting golden whistler Pachycephala pectoralis is comparatively flexible in its diet and niche requirements as compared to the henna-tailed jungle-flycatcher Cyornis colonus, which is largely restricted to the forest interior. Using population-genomic and coalescent approaches, we estimated levels of gene flow, population differentiation and divergence time between the two island populations. We observed higher levels of differentiation, an approximately two to four times deeper divergence time and near-zero levels of gene flow between the two island populations of the more forest-dependent henna-tailed jungle-flycatcher as compared to the more generalist golden whistler. Our results suggest that Quaternary land bridges act as semipermeable agents of gene flow in Wallacea, allowing only certain taxa to connect between islands while others remain isolated. Quaternary land bridges do not accommodate all terrestrial species equally, differing in suitability according to life history and species biology. More generalist species are likely to use Quaternary land connections as a conduit for gene flow between islands whereas island populations of more specialist species may continue to be reproductively isolated even during periods of Quaternary land bridges.

Coalescent and Models of Identity by Descent

The concept of identity by descent, or IBD, traces back about a century to the early days of population and evolutionary genetics. Since that time, IBD has been employed in a number of ways, but always as a measure of close genetic relationship. Compared to other ways of quantifying genetic variation, IBD focuses its power on the recent ancestry of populations and species. Observations of IBD can be used to make inferences about events of interest to evolutionary biologists studying a wide range of taxa, especially now that the sequencing of genomes has become practical.