Inferring the model and onset of natural selection under varying population size from the site frequency spectrum and haplotype structure.

A fundamental question about adaptation in a population is the time of onset of the selective pressure acting on beneficial alleles. Inferring this time, in turn, depends on the selection model. We develop a framework of approximate Bayesian computation (ABC) that enables the use of the full site frequency spectrum and haplotype structure to test the goodness-of-fit of selection models and estimate the timing of selection under varying population size scenarios.

A population genetic interpretation of GWAS findings for human quantitative traits.

Human genome-wide association studies (GWASs) are revealing the genetic architecture of anthropomorphic and biomedical traits, i.e., the frequencies and effect sizes of variants that contribute to heritable variation in a trait.

Extremely high genetic diversity in a single tumor points to prevalence of non-Darwinian cell evolution.

The prevailing view that the evolution of cells in a tumor is driven by Darwinian selection has never been rigorously tested. Because selection greatly affects the level of intratumor genetic diversity, it is important to assess whether intratumor evolution follows the Darwinian or the non-Darwinian mode of evolution. To provide the statistical power, many regions in a single tumor need to be sampled and analyzed much more extensively than has been attempted in previous intratumor studies.

Estimating the Ages of Selection Signals from Different Epochs in Human History.

Genetic variation harbors signatures of natural selection driven by selective pressures that are often unknown. Estimating the ages of selection signals may allow reconstructing the history of environmental changes that shaped human phenotypes and diseases. We have developed an approximate Bayesian computation (ABC) approach to estimate allele ages under a model of selection on new mutations and under demographic models appropriate for human populations.

A New Proof of the Expected Frequency Spectrum under the Standard Neutral Model.

The sample frequency spectrum is an informative and frequently employed approach for summarizing DNA variation data. Under the standard neutral model the expectation of the sample frequency spectrum has been derived by at least two distinct approaches. One relies on using results from diffusion approximations to the Wright-Fisher Model.

The long-term maintenance of a resistance polymorphism through diffuse interactions.

Plant resistance (R) genes are a crucial component in plant defence against pathogens. Although R genes often fail to provide durable resistance in an agricultural context, they frequently persist as long-lived balanced polymorphisms in nature. Standard theory explains the maintenance of such polymorphisms through a balance of the costs and benefits of resistance and virulence in a tightly coevolving host-pathogen pair.

A genetic polymorphism evolving in parallel in two cell compartments and in two clades.

Background: The enzyme phosphoenolpyruvate carboxykinase, PEPCK, occurs in its guanosine-nucleotide-using form in animals and a few prokaryotes. We study its natural genetic variation in Colias (Lepidoptera, Pieridae). PEPCK offers a route, alternative to pyruvate kinase, for carbon skeletons to move between cytosolic glycolysis and mitochondrial Krebs cycle reactions.

Evaluating the evidence for transmission distortion in human pedigrees.

Children of a heterozygous parent are expected to carry either allele with equal probability. Exceptions can occur, however, due to meiotic drive, competition among gametes, or viability selection, which we collectively term "transmission distortion" (TD). Although there are several well-characterized examples of these phenomena, their existence in humans remains unknown.

Two evolutionary histories in the genome of rice: the roles of domestication genes.

Genealogical patterns in different genomic regions may be different due to the joint influence of gene flow and selection. The existence of two subspecies of cultivated rice provides a unique opportunity for analyzing these effects during domestication. We chose 66 accessions from the three rice taxa (about 22 each from Oryza sativa indica, O.

Human endogenous retrovirus K106 (HERV-K106) was infectious after the emergence of anatomically modern humans.

HERV-K113 and HERV-K115 have been considered to be among the youngest HERVs because they are the only known full-length proviruses that are insertionally polymorphic and maintain the open reading frames of their coding genes. However, recent data suggest that HERV-K113 is at least 800,000 years old, and HERV-K115 even older. A systematic study of HERV-K HML2 members to identify HERVs that may have infected the human genome in the more recent evolutionary past is lacking.

A reduced representation approach to population genetic analyses and applications to human evolution.

Second-generation sequencing technologies allow surveys of sequence variation on an unprecedented scale. However, despite the rapid decrease in sequencing costs, collecting whole-genome sequence data on a population scale is still prohibitive for many laboratories. We have implemented an inexpensive, reduced representation protocol for preparing resequencing targets, and we have developed the analytical tools necessary for making population genetic inferences.

The variance of coalescent time estimates from DNA sequences.

A method for estimating the variance of a simple estimate of the time of the most recent common ancestor of a sample of sequences is described.

Adaptive genic evolution in the Drosophila genomes.

Determining the extent of adaptive evolution at the genomic level is central to our understanding of molecular evolution. A suitable observation for this purpose would consist of polymorphic data on a large and unbiased collection of genes from two closely related species, each having a large and stable population. In this study, we sequenced 419 genes from 24 lines of Drosophila melanogaster and its close relatives.

Evidence that the adaptive allele of the brain size gene microcephalin introgressed into Homo sapiens from an archaic Homo lineage.

At the center of the debate on the emergence of modern humans and their spread throughout the globe is the question of whether archaic Homo lineages contributed to the modern human gene pool, and more importantly, whether such contributions impacted the evolutionary adaptation of our species. A major obstacle to answering this question is that low levels of admixture with archaic lineages are not expected to leave extensive traces in the modern human gene pool because of genetic drift. Loci that have undergone strong positive selection, however, offer a unique opportunity to identify low-level admixture with archaic lineages, provided that the introgressed archaic allele has risen to high frequency under positive selection.

Interrogating multiple aspects of variation in a full resequencing data set to infer human population size changes.

We present an expanded data set of 50 unlinked autosomal noncoding regions, resequenced in samples of Hausa from Cameroon, Italians, and Chinese. We use these data to make inferences about human demographic history by using a technique that combines multiple aspects of genetic data, including levels of polymorphism, the allele frequency spectrum, and linkage disequilibrium. We explore an extensive range of demographic parameters and demonstrate that our method of combining multiple aspects of the data results in a significant reduction of the compatible parameter space.

An evolutionary framework for common diseases: the ancestral-susceptibility model.

Unlike rare mendelian diseases, which are due to new mutations (i.e. derived alleles), several alleles that increase the risk to common diseases are ancestral.

Ongoing adaptive evolution of ASPM, a brain size determinant in Homo sapiens.

The gene ASPM (abnormal spindle-like microcephaly associated) is a specific regulator of brain size, and its evolution in the lineage leading to Homo sapiens was driven by strong positive selection. Here, we show that one genetic variant of ASPM in humans arose merely about 5800 years ago and has since swept to high frequency under strong positive selection. These findings, especially the remarkably young age of the positively selected variant, suggest that the human brain is still undergoing rapid adaptive evolution.

Microcephalin, a gene regulating brain size, continues to evolve adaptively in humans.

The gene Microcephalin (MCPH1) regulates brain size and has evolved under strong positive selection in the human evolutionary lineage. We show that one genetic variant of Microcephalin in modern humans, which arose approximately 37,000 years ago, increased in frequency too rapidly to be compatible with neutral drift. This indicates that it has spread under strong positive selection, although the exact nature of the selection is unknown.

Maximum-likelihood estimation of demographic parameters using the frequency spectrum of unlinked single-nucleotide polymorphisms.

A maximum-likelihood method for demographic inference is applied to data sets consisting of the frequency spectrum of unlinked single-nucleotide polymorphisms (SNPs). We use simulation analyses to explore the effect of sample size and number of polymorphic sites on both the power to reject the null hypothesis of constant population size and the properties of two- and three-dimensional maximum-likelihood estimators (MLEs). Large amounts of data are required to produce accurate demographic inferences, particularly for scenarios of recent growth.

Comparative linkage-disequilibrium analysis of the beta-globin hotspot in primates.

Recombination rates vary both across the genome and between different species, but little information is available about the temporal and physical scales over which such rates change. To shed light on these questions, we performed a high-resolution analysis of a genomic region within the beta-globin gene cluster that is known to experience elevated recombination rates in humans. For this purpose, we developed new linkage disequilibrium-based methods that thoroughly search for subsets of the data with unusually high or unusually low estimated values of the population-recombination parameter (4Nr, where N is the effective population size and r is the crossover rate between adjacent base pairs).

Stochasticity overrules the "three-times rule": genetic drift, genetic draft, and coalescence times for nuclear loci versus mitochondrial DNA.

Palumbi et al. (2001) proposed a "three-times rule" that uses mitochondrial DNA (mtDNA) sequences to predict probabilities of monophyly for nuclear loci (i.e.

Mathematical consequences of the genealogical species concept.

A genealogical species is defined as a basal group of organisms whose members are all more closely related to each other than they are to any organisms outside the group ("exclusivity"), and which contains no exclusive group within it. In practice, a pair of species is so defined when phylogenies of alleles from a sample of loci shows them to be reciprocally monophyletic at all or some specified fraction of the loci. We investigate the length of time it takes to attain this status when an ancestral population divides into two descendant populations of equal size with no gene exchange, and when genetic drift and mutation are the only evolutionary forces operating.

Inferences about human demography based on multilocus analyses of noncoding sequences.

Data from 10 unlinked autosomal noncoding regions, resequenced in 15 individuals from each of three populations, were used in a multilocus analysis to test models of human demography. Each of the 10 regions consisted of approximately 2500 bp. The multilocus analysis, based on summary statistics (average and variance of Tajima's D and Fu and Li's D*), was used to test a family of models with recent population expansion.

Geographic and haplotype structure of candidate type 2 diabetes susceptibility variants at the calpain-10 locus.

Recently, a positional cloning study proposed that haplotypes at the calpain-10 locus (CAPN10) are associated with increased risk of type 2 diabetes, or non-insulin-dependent diabetes mellitus, in Mexican Americans, Finns, and Germans. To inform the interpretation of the original mapping results and to look for evidence for the action of natural selection on CAPN10, we undertook a population-based genotyping survey of the candidate susceptibility variants. First, we genotyped sites 43, 19, and 63 (the haplotype-defining variants previously proposed) and four closely linked SNPs, in 561 individuals from 11 populations from five continents, and we examined the linkage disequilibrium among them.