Related papers: Detecting range expansions from genetic data
We describe a statistical hypothesis test for the presence of a signal. The test allows the researcher to fix the signal location and/or width a priori, or perform a search to find the signal region that maximizes the signal. The background…
We consider in this paper detection of signal regions associated with disease outcomes in whole genome association studies. Gene- or region-based methods have become increasingly popular in whole genome association analysis as a…
Spatial models where growth is limited to the edge of the expansions have been instrumental to understand the population dynamics and the clone size distribution in growing cellular populations, such as microbial colonies and avascular…
Second-generation sequencing technologies have replaced array-based technologies and become the default method for genomics and epigenomics analysis. Second-generation sequencing technologies sequence tens of millions of DNA/cDNA fragments…
As whole genomes become widely available, maximum likelihood and Bayesian phylogenetic methods are demonstrating their limits in meeting the escalating computational demands. Conversely, distance-based phylogenetic methods are efficient,…
The site frequency spectrum (SFS) is a widely used summary statistic of genomic data. Motivated by recent evidence for the role of neutral evolution in cancer, we investigate the SFS of neutral mutations in an exponentially growing…
Understanding the dynamics of genome rearrangements is a major issue of phylogenetics. Phylogenetics is the study of species evolution. A major goal of the field is to establish evolutionary relationships within groups of species, in order…
Gene regulatory network inference (GRNI) aims to discover how genes causally regulate each other from gene expression data. It is well-known that statistical dependencies in observed data do not necessarily imply causation, as spurious…
Pathogen genome data offers valuable structure for spatial models, but its utility is limited by incomplete sequencing coverage. We propose a probabilistic framework for inferring genetic distances between unsequenced cases and known…
Many theoretical and experimental studies suggest that range expansions can have severe consequences for the gene pool of the expanding population. Due to strongly enhanced genetic drift at the advancing frontier, neutral and weakly…
The modeling of time series is becoming increasingly critical in a wide variety of applications. Overall, data evolves by following different patterns, which are generally caused by different user behaviors. Given a time series, we define…
We study a universal object for the genealogy of a sample in populations with mutations: the critical birth-death process with Poissonian mutations, conditioned on its population size at a fixed time horizon. We show how this process arises…
We consider an expanding population on the plane. The genealogy of a sample from the population is modelled by coalescing Brownian motion on the circle. We establish a weak law of large numbers for the site frequency spectrum in this model.…
Phylogenies are commonly used to represent the evolutionary relationships between species, and often these phylogenies are equipped with edge lengths that indicate degrees of evolutionary difference. Given such a phylogeny, a popular…
The paper reviews the results obtained for spatial population models and the evolution of the genealogies of these populations during the last decade by the author and his coworkers. The focus is on their large scale behaviour and on the…
The shape of allele-frequency clines maintained by migration-selection balance depends not only on the properties of migration and selection, but also on the dominance relations among alleles and on linkage to other loci under selection. We…
We introduce a low dimensional function of the site frequency spectrum that is tailor-made for distinguishing coalescent models with multiple mergers from Kingman coalescent models with population growth, and use this function to construct…
A key step in any resonant anomaly detection search is accurate modeling of the background distribution in each signal region. Data-driven methods like CATHODE accomplish this by training separate generative models on the complement of each…
Results on the spectral behavior of random matrices as the dimension increases are applied to the problem of detecting the number of sources impinging on an array of sensors. A common strategy to solve this problem is to estimate the…
Recent microbial experiments suggest that enhanced genetic drift at the frontier of a two-dimensional range expansion can cause genetic sectoring patterns with fractal domain boundaries. Here, we propose and analyze a simple model of…