Diffusion Approximations for Demographic Inference: DaDi

Models of demographic history (population sizes, migration rates, and divergence times) inferred from genetic data complement archeology and serve as null models in genome scans for selection. Most current inference methods are computationally limited to considering simple models or non-recombining data. We introduce a method based on a diffusion approximation to the joint frequency spectrum of genetic variation between populations. Our implementation, DaDi, can model up to three interacting populations and scales well to genome-wide data. We have applied DaDi to human data from Africa, Europe, and East Asia, building the most complex statistically well-characterized model of human migration out of Africa to date

Inferring the Joint Demographic History of Multiple Populations from Multidimensional SNP Frequency Data

Demographic models built from genetic data play important roles in illuminating prehistorical events and serving as null models in genome scans for selection. We introduce an inference method based on the joint frequency spectrum of genetic variants within and between populations. For candidate models we numerically compute the expected spectrum using a diffusion approximation to the one-locus, two-allele Wright-Fisher process, involving up to three simultaneous populations. Our approach is a composite likelihood scheme, since linkage between neutral loci alters the variance but not the expectation of the frequency spectrum. We thus use bootstraps incorporating linkage to estimate uncertainties for parameters and significance values for hypothesis tests. Our method can also incorporate selection on single sites, predicting the joint distribution of selected alleles among populations experiencing a bevy of evolutionary forces, including expansions, contractions, migrations, and admixture. We model human expansion out of Africa and the settlement of the New World, using 5 Mb of noncoding DNA resequenced in 68 individuals from 4 populations (YRI, CHB, CEU, and MXL) by the Environmental Genome Project. We infer divergence between West African and Eurasian populations 140 thousand years ago (95% confidence interval: 40–270 kya). This is earlier than other genetic studies, in part because we incorporate migration. We estimate the European (CEU) and East Asian (CHB) divergence time to be 23 kya (95% c.i.: 17–43 kya), long after archeological evidence places modern humans in Europe. Finally, we estimate divergence between East Asians (CHB) and Mexican-Americans (MXL) of 22 kya (95% c.i.: 16.3–26.9 kya), and our analysis yields no evidence for subsequent migration. Furthermore, combining our demographic model with a previously estimated distribution of selective effects among newly arising amino acid mutations accurately predicts the frequency spectrum of nonsynonymous variants across three continental populations (YRI, CHB, CEU).</p

High-resolution and high-accuracy topographic and transcriptional maps of the nucleosome barrier.

Nucleosomes represent mechanical and energetic barriers that RNA Polymerase II (Pol II) must overcome during transcription. A high-resolution description of the barrier topography, its modulation by epigenetic modifications, and their effects on Pol II nucleosome crossing dynamics, is still missing. Here, we obtain topographic and transcriptional (Pol II residence time) maps of canonical, H2A.Z, and monoubiquitinated H2B (uH2B) nucleosomes at near base-pair resolution and accuracy. Pol II crossing dynamics are complex, displaying pauses at specific loci, backtracking, and nucleosome hopping between wrapped states. While H2A.Z widens the barrier, uH2B heightens it, and both modifications greatly lengthen Pol II crossing time. Using the dwell times of Pol II at each nucleosomal position we extract the energetics of the barrier. The orthogonal barrier modifications of H2A.Z and uH2B, and their effects on Pol II dynamics rationalize their observed enrichment in +1 nucleosomes and suggest a mechanism for selective control of gene expression

SNP identification, verification, and utility for population genetics in a non-model genus

<p>Abstract</p> <p>Background</p> <p>By targeting SNPs contained in both coding and non-coding areas of the genome, we are able to identify genetic differences and characterize genome-wide patterns of variation among individuals, populations and species. We investigated the utility of 454 sequencing and MassARRAY genotyping for population genetics in natural populations of the teleost, <it>Fundulus heteroclitus </it>as well as closely related <it>Fundulus </it>species (<it>F. grandis</it>, <it>F. majalis </it>and <it>F. similis</it>).</p> <p>Results</p> <p>We used 454 pyrosequencing and MassARRAY genotyping technology to identify and type 458 genome-wide SNPs and determine genetic differentiation within and between populations and species of <it>Fundulus</it>. Specifically, pyrosequencing identified 96 putative SNPs across coding and non-coding regions of the <it>F. heteroclitus </it>genome: 88.8% were verified as true SNPs with MassARRAY. Additionally, putative SNPs identified in <it>F. heteroclitus </it>EST sequences were verified in most (86.5%) <it>F. heteroclitus </it>individuals; fewer were genotyped in <it>F. grandis </it>(74.4%), <it>F. majalis </it>(72.9%), and <it>F. similis </it>(60.7%) individuals. SNPs were polymorphic and showed latitudinal clinal variation separating northern and southern populations and established isolation by distance in <it>F. heteroclitus </it>populations. In <it>F. grandis</it>, SNPs were less polymorphic but still established isolation by distance. Markers differentiated species and populations.</p> <p>Conclusions</p> <p>In total, these approaches were used to quickly determine differences within the <it>Fundulus </it>genome and provide markers for population genetic studies.</p

A Tale of Two Haplotypes: The \u3cem\u3eEDA2R/AR\u3c/em\u3e Intergenic Region is the most Divergent Genomic Segment between Africans and East Asians in the Human Genome

Single nucleotide polymorphisms (SNPs) with large allele frequency differences between human populations are relatively rare. The longest run of SNPs with an allele frequency difference of one between the Yoruba of Nigeria and the Han Chinese is found on the long arm of the X chromosome in the intergenic region separating the EDA2R and AR genes. It has been proposed that the unusual allele frequency distributions of these SNPs are the result of a selective sweep affecting African populations that occurred after the Out-of-Africa migration. To investigate the evolutionary history of the EDA2R/AR intergenic region, we characterized the haplotype structure of 52 of its highly-differentiated SNPs. Using a publicly-available dataset of 3,000 X chromosomes from 65 human populations, we found that nearly all human X chromosomes carry one of two modal haplotypes for these 52 SNPs. The predominance of two highly divergent haplotypes at this locus was confirmed using a subset of individuals sequenced to high coverage. The first of these haplotypes, the α haplotype, is at high frequencies in most of the African populations surveyed and likely arose prior to the separation of African populations into distinct genetic entities. The second, the β haplotype, is frequent or fixed in all non-African populations and likely arose in East Africa prior to the Out-of-Africa migration. We also observed a small group of rare haplotypes with no clear relationship to the α and β haplotypes. These haplotypes occur at relatively high frequencies in African hunter-gatherer populations, like the San and Mbuti Pygmies. Our analysis indicates that these haplotypes are part of a pool of diverse, ancestral haplotypes that have now been almost entirely replaced by the α and β haplotypes. We suggest that the rise of the α and β haplotypes was the result of the demographic forces that human populations experienced during the formation of modern African populations and the Out-of-Africa migration. However, we also present evidence that this region is the target of selection in the form of positive selection on the α and β haplotypes and of purifying selection against α/β recombinants