A sampling algorithm for segregation analysis

Methods for detecting Quantitative Trait Loci (QTL) without markers have generally used iterative peeling algorithms for determining genotype probabilities. These algorithms have considerable shortcomings in complex pedigrees. A Monte Carlo Markov chain (MCMC) method which samples the pedigree of the whole population jointly is described. Simultaneous sampling of the pedigree was achieved by sampling descent graphs using the Metropolis-Hastings algorithm. A descent graph describes the inheritance state of each allele and provides pedigrees guaranteed to be consistent with Mendelian sampling. Sampling descent graphs overcomes most, if not all, of the limitations incurred by iterative peeling algorithms. The algorithm was able to find the QTL in most of the simulated populations. However, when the QTL was not modeled or found then its effect was ascribed to the polygenic component. No QTL were detected when they were not simulated.     

Estimating Variable Effective Population Sizes from Multiple Genomes: A Sequentially Markov Conditional Sampling Distribution Approach

Throughout history, the population size of modern humans has varied considerably due to changes in environment, culture, and technology. More accurate estimates of population size changes, and when they occurred, should provide a clearer picture of human colonization history and help remove confounding effects from natural selection inference. Demography influences the pattern of genetic variation in a population, and thus genomic data of multiple individuals sampled from one or more present-day populations contain valuable information about the past demographic history. Recently, Li and Durbin developed a coalescent-based hidden Markov model, called the pairwise sequentially Markovian coalescent (PSMC), for a pair of chromosomes (or one diploid individual) to estimate past population sizes. This is an efficient, useful approach, but its accuracy in the very recent past is hampered by the fact that, because of the small sample size, only few coalescence events occur in that period. Multiple genomes from the same population contain more information about the recent past, but are also more computationally challenging to study jointly in a coalescent framework. Here, we present a new coalescent-based method that can efficiently infer population size changes from multiple genomes, providing access to a new store of information about the recent past. Our work generalizes the recently developed sequentially Markov conditional sampling distribution framework, which provides an accurate approximation of the probability of observing a newly sampled haplotype given a set of previously sampled haplotypes. Simulation results demonstrate that we can accurately reconstruct the true population histories, with a significant improvement over the PSMC in the recent past. We apply our method, called diCal, to the genomes of multiple human individuals of European and African ancestry to obtain a detailed population size change history during recent times.     

Analyses of genetic ancestry enable key insights for molecular ecology

Gene flow and recombination in admixed populations produce genomes that are mosaic combinations of chromosome segments inherited from different source populations, that is, chromosome segments with different genetic ancestries. The statistical problem of estimating genetic ancestry from DNA sequence data has been widely studied, and analyses of genetic ancestry have facilitated research in molecular ecology and ecological genetics. In this review, we describe and compare different model‐based statistical methods used to infer genetic ancestry. We describe the conceptual and mathematical structure of these models and highlight some of their key differences and shared features. We then discuss recent empirical studies that use estimates of genetic ancestry to analyse population histories, the nature and genetic basis of species boundaries, and the genetic architecture of traits. These diverse studies demonstrate the breadth of applications that rely on genetic ancestry estimates and typify the genomics‐enabled research that is becoming increasingly common in molecular ecology. We conclude by identifying key research areas where future studies might further advance this field.     

Three-dimensional chromatin ensemble reconstruction via stochastic embedding

1. Download : Download high-res image (359KB)
2. Download : Download full-size image

     

The genealogical process of neutral genes with mutation in geographically structured populations

The genealogical process of neutral genes with mutation in geographically structured populations is investigated. Following Watterson [24], the sampled genes are partitioned into two types, old equivalence classes and new equivalence classes. The model is described by a bivariate continuous time Markov chain as an interactive particle system. Some results are obtained in the two-population model and the stepping stone model with symmetric nearest-neighbour migration.     

The strong-migration limit for the genealogical process in geographically structured populations

In this article we assume that the entire population is subdivided into a finite number of panmictic colonies, each of which consists of a respective number of haploid individuals. We also assume that random genetic drift occurs in each colony and migration among colonies, which is independent of time and ergodic. We study the genealogical process of sampled genes from geographically structured populations. We prove that if the actual total population number is replaced by the effective population number, the mean coalescence time converges to that in a panmictic population in the strong migration limit. We also obtain the geographical distribution of the common ancestor.