Estimating barriers to gene flow from distorted isolation by distance patterns

ABSTRACTIn continuous populations with local migration, nearby pairs of individuals have on average more similar genotypes than geographically well separated pairs. A barrier to gene flow distorts this classical pattern of isolation by distance. Genetic similarity is decreased for sample pairs on different sides of the barrier and increased for pairs on the same side near the barrier. Here, we introduce an inference scheme that utilizes this signal to detect and estimate the strength of a linear barrier to gene flow in two-dimensions. We use a diffusion approximation to model the effects of a barrier on the geographical spread of ancestry backwards in time. This approach allows us to calculate the chance of recent coalescence and probability of identity by descent. We introduce an inference scheme that fits these theoretical results to the geographical covariance structure of bialleleic genetic markers. It can estimate the strength of the barrier as well as several demographic parameters. We investigate the power of our inference scheme to detect barriers by applying it to a wide range of simulated data. We also showcase an example application to a Antirrhinum majus (snapdragon) flower color hybrid zone, where we do not detect any signal of a strong genome wide barrier to gene flow.

Scattering and delay time for 1D asymmetric potentials: The step-linear and the step-exponential case

We analyze the quantum-mechanical behavior of a system described by a one-dimensional asymmetric potential constituted by a step plus (i) a linear barrier or (ii) an exponential barrier. We solve the energy eigenvalue equation by means of the integral representation method, classifying the independent solutions as equivalence classes of homotopic paths in the complex plane. We discuss the structure of the bound states as function of the height U 0 of the step and we study the propagation of a sharp-peaked wave packet reflected by the barrier. For both the linear and the exponential barrier we provide an explicit formula for the delay time τ ( E ) as a function of the peak energy E. We display the resonant behavior of τ ( E ) at energies close to U 0. By analyzing the asymptotic behavior for large energies of the eigenfunctions of the continuous spectrum we also show that, as expected, τ ( E ) approaches the classical value for E → ∞, thus diverging for the step-linear case and vanishing for the step-exponential one.

On the Distribution of First Exit Time for Brownian Motion with Double Linear Time-Dependent Barriers

This paper focuses on the first exit time for a Brownian motion with a double linear time-dependent barrier specified by y=a+bt, y=ct, (a>0, b<0, c>0). We are concerned in this paper with the distribution of the Brownian motion hitting the upper barrier before hitting the lower linear barrier. The main method we applied here is the Girsanov transform formula. As a result, we expressed the density of such exit time in terms of a finite series. This result principally provides us an analytical expression for the distribution of the aforementioned exit time and an easy way to compute the distribution of first exit time numerically.