scholarly journals Estimating barriers to gene flow from distorted isolation by distance patterns

2017 ◽  
Author(s):  
Harald Ringbauer ◽  
Alexander Kolesnikov ◽  
David Field ◽  
Nicholas H. Barton

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.

PeerJ ◽  
2021 ◽  
Vol 9 ◽  
pp. e11425
Author(s):  
Sofia G. Seabra ◽  
Ana S.B. Rodrigues ◽  
Sara E. Silva ◽  
Ana Carina Neto ◽  
Francisco Pina-Martins ◽  
...  

Understanding patterns of population differentiation and gene flow in insect vectors of plant diseases is crucial for the implementation of management programs of disease. We investigated morphological and genome-wide variation across the distribution range of the spittlebug Philaenus spumarius (Linnaeus, 1758) (Hemiptera, Auchenorrhyncha, Aphrophoridae), presently the most important vector of the plant pathogenic bacterium Xylella fastidiosa Wells et al., 1987 in Europe. We found genome-wide divergence between P. spumarius and a very closely related species, P. tesselatus Melichar, 1899, at RAD sequencing markers. The two species may be identified by the morphology of male genitalia but are not differentiated at mitochondrial COI, making DNA barcoding with this gene ineffective. This highlights the importance of using integrative approaches in taxonomy. We detected admixture between P. tesselatus from Morocco and P. spumarius from the Iberian Peninsula, suggesting gene-flow between them. Within P. spumarius, we found a pattern of isolation-by-distance in European populations, likely acting alongside other factors restricting gene flow. Varying levels of co-occurrence of different lineages, showing heterogeneous levels of admixture, suggest other isolation mechanisms. The transatlantic populations of North America and Azores were genetically closer to the British population analyzed here, suggesting an origin from North-Western Europe, as already detected with mitochondrial DNA. Nevertheless, these may have been produced through different colonization events. We detected SNPs with signatures of positive selection associated with environmental variables, especially related to extremes and range variation in temperature and precipitation. The population genomics approach provided new insights into the patterns of divergence, gene flow and adaptation in these spittlebugs and led to several hypotheses that require further local investigation.


2019 ◽  
Author(s):  
Constantin Ammar ◽  
Markus Gruber ◽  
Gergely Csaba ◽  
Ralf Zimmer

1AbstractMass spectrometry based proteomics is the method of choice for quantifying genome-wide differential changes of proteins in a wide range of biological and biomedical applications. Protein changes need to be reliably derived from a large number of measured peptide intensities and their corresponding fold changes. These fold changes vary considerably for a given protein.Numerous instrumental setups aim to reduce this variability, while current computational methods only implicitly account for this problem. We introduce a new method, MS-EmpiRe (github.com/zimmerlab/MS-EmpiRe), which explicitly accounts for the noise underlying peptide fold changes. We derive dataset-specific, intensity-dependent empirical error distributions, which are used for individual weighing of peptide fold changes to detect differentially abundant proteins. The method requires only peptide intensities mapped to proteins and, thus, can be applied to any common quantitative proteomics setup. In a recently published proteome-wide benchmarking dataset, MS-EmpiRe doubles the number of correctly identified changing proteins at a correctly estimated FDR cutoff in comparison to state-of-the-art tools. We confirm the superior performance of MS-EmpiRe on simulated data. MS-EmpiRe provides rapid processing (< 2min) and is an easy to use, general-purpose tool.


2019 ◽  
Vol 36 (9) ◽  
pp. 1849-1861 ◽  
Author(s):  
Mário Vicente ◽  
Mattias Jakobsson ◽  
Peter Ebbesen ◽  
Carina M Schlebusch

Abstract Southern African indigenous groups, traditionally hunter-gatherers (San) and herders (Khoekhoe), are commonly referred to as “Khoe-San” populations and have a long history in southern Africa. Their ancestors were largely isolated up until ∼2,000 years ago before the arrival of pastoralists and farmers in southern Africa. Assessing relationships among regional Khoe-San groups has been challenging due to admixture with immigrant populations that obscure past population affinities and gene flow among these autochthonous communities. We re-evaluate a combined genome-wide data set of previously published southern Africa Khoe-San populations in conjunction with novel data from Khoe-San individuals collected in Xade (Central Kalahari Game Reserve, Botswana) prior to their resettlement outside the reserve. After excluding regions in the genome that trace their ancestry to recent migrant groups, the genetic diversity of 20 Khoe-San groups fitted an isolation-by-distance model. Even though isolation-by-distance explained most genetic affinities between the different autochthonous groups, additional signals of contact between Khoe-San groups could be detected. For instance, we found stronger genetic affinities, than what would be explained by isolation-by-distance gene flow, between the two geographically separated Khoe-San groups, who speak branches of the Kx’a-language family (ǂHoan and Ju). We also scanned the genome-wide data for signals of adaptive gene flow from farmers/herders into Khoe-San groups and identified a number of genomic regions potentially introduced by the arrival of the new groups. This study provides a comprehensive picture of affinities among Khoe-San groups, prior to the arrival of recent migrants, and found that these affinities are primarily determined by the geographic landscape.


2019 ◽  
Vol 12 (1) ◽  
Author(s):  
Maysa Tiemi Motoki ◽  
Dina Madera Fonseca ◽  
Elliott Frederic Miot ◽  
Bruna Demari-Silva ◽  
Phoutmany Thammavong ◽  
...  

Abstract Background The Asian tiger mosquito, Aedes (Stegomyia) albopictus (Skuse) is an important worldwide invasive species and can be a locally important vector of chikungunya, dengue and, potentially, Zika. This species is native to Southeast Asia where populations thrive in both temperate and tropical climates. A better understanding of the population structure of Ae. albopictus in Lao PDR is very important in order to support the implementation of strategies for diseases prevention and vector control. In the present study, we investigated the genetic variability of Ae. albopictus across a north-south transect in Lao PDR. Methods We used variability in a 1337-bp fragment of the mitochondrial cytochrome c oxidase subunit 1 gene (cox1), to assess the population structure of Ae. albopictus in Lao PDR. For context, we also examined variability at the same genetic locus in samples of Ae. albopictus from Thailand, China, Taiwan, Japan, Singapore, Italy and the USA. Results We observed very high levels of genetic polymorphism with 46 novel haplotypes in Ae. albopictus from 9 localities in Lao PDR and Thailand populations. Significant differences were observed between the Luangnamtha population and other locations in Lao PDR. However, we found no evidence of isolation by distance. There was overall little genetic structure indicating ongoing and frequent gene flow among populations or a recent population expansion. Indeed, the neutrality test supported population expansion in Laotian Ae. albopictus and mismatch distribution analyses showed a lack of low frequency alleles, a pattern often seen in bottlenecked populations. When samples from Lao PDR were analyzed together with samples from Thailand, China, Taiwan, Japan, Singapore, Italy and the USA, phylogenetic network and Bayesian cluster analysis showed that most populations from tropical/subtropical regions are more genetically related to each other, than populations from temperate regions. Similarly, most populations from temperate regions are more genetically related to each other, than those from tropical/subtropical regions. Conclusions Aedes albopictus in Lao PDR are genetically related to populations from tropical/subtropical regions (i.e. Thailand, Singapore, and California and Texas in the USA). The extensive gene flow among locations in Lao PDR indicates that local control is undermined by repeated introductions from untreated sites.


2021 ◽  
Vol 0 (0) ◽  
Author(s):  
Boris Kargoll ◽  
Alexander Dorndorf ◽  
Mohammad Omidalizarandi ◽  
Jens-André Paffenholz ◽  
Hamza Alkhatib

Abstract In this contribution, a vector-autoregressive (VAR) process with multivariate t-distributed random deviations is incorporated into the Gauss-Helmert model (GHM), resulting in an innovative adjustment model. This model is versatile since it allows for a wide range of functional models, unknown forms of auto- and cross-correlations, and outlier patterns. Subsequently, a computationally convenient iteratively reweighted least squares method based on an expectation maximization algorithm is derived in order to estimate the parameters of the functional model, the unknown coefficients of the VAR process, the cofactor matrix, and the degree of freedom of the t-distribution. The proposed method is validated in terms of its estimation bias and convergence behavior by means of a Monte Carlo simulation based on a GHM of a circle in two dimensions. The methodology is applied in two different fields of application within engineering geodesy: In the first scenario, the offset and linear drift of a noisy accelerometer are estimated based on a Gauss-Markov model with VAR and multivariate t-distributed errors, as a special case of the proposed GHM. In the second scenario real laser tracker measurements with outliers are adjusted to estimate the parameters of a sphere employing the proposed GHM with VAR and multivariate t-distributed errors. For both scenarios the estimated parameters of the fitted VAR model and multivariate t-distribution are analyzed for evidence of auto- or cross-correlations and deviation from a normal distribution regarding the measurement noise.


2012 ◽  
Vol 696 ◽  
pp. 228-262 ◽  
Author(s):  
A. Kourmatzis ◽  
J. S. Shrimpton

AbstractThe fundamental mechanisms responsible for the creation of electrohydrodynamically driven roll structures in free electroconvection between two plates are analysed with reference to traditional Rayleigh–Bénard convection (RBC). Previously available knowledge limited to two dimensions is extended to three-dimensions, and a wide range of electric Reynolds numbers is analysed, extending into a fully inherently three-dimensional turbulent regime. Results reveal that structures appearing in three-dimensional electrohydrodynamics (EHD) are similar to those observed for RBC, and while two-dimensional EHD results bear some similarities with the three-dimensional results there are distinct differences. Analysis of two-point correlations and integral length scales show that full three-dimensional electroconvection is more chaotic than in two dimensions and this is also noted by qualitatively observing the roll structures that arise for both low (${\mathit{Re}}_{E} = 1$) and high electric Reynolds numbers (up to ${\mathit{Re}}_{E} = 120$). Furthermore, calculations of mean profiles and second-order moments along with energy budgets and spectra have examined the validity of neglecting the fluctuating electric field ${ E}_{i}^{\ensuremath{\prime} } $ in the Reynolds-averaged EHD equations and provide insight into the generation and transport mechanisms of turbulent EHD. Spectral and spatial data clearly indicate how fluctuating energy is transferred from electrical to hydrodynamic forms, on moving through the domain away from the charging electrode. It is shown that ${ E}_{i}^{\ensuremath{\prime} } $ is not negligible close to the walls and terms acting as sources and sinks in the turbulent kinetic energy, turbulent scalar flux and turbulent scalar variance equations are examined. Profiles of hydrodynamic terms in the budgets resemble those in the literature for RBC; however there are terms specific to EHD that are significant, indicating that the transfer of energy in EHD is also attributed to further electrodynamic terms and a strong coupling exists between the charge flux and variance, due to the ionic drift term.


2020 ◽  
Vol 146 ◽  
pp. 03004
Author(s):  
Douglas Ruth

The most influential parameter on the behavior of two-component flow in porous media is “wettability”. When wettability is being characterized, the most frequently used parameter is the “contact angle”. When a fluid-drop is placed on a solid surface, in the presence of a second, surrounding fluid, the fluid-fluid surface contacts the solid-surface at an angle that is typically measured through the fluid-drop. If this angle is less than 90°, the fluid in the drop is said to “wet” the surface. If this angle is greater than 90°, the surrounding fluid is said to “wet” the surface. This definition is universally accepted and appears to be scientifically justifiable, at least for a static situation where the solid surface is horizontal. Recently, this concept has been extended to characterize wettability in non-static situations using high-resolution, two-dimensional digital images of multi-component systems. Using simple thought experiments and published experimental results, many of them decades old, it will be demonstrated that contact angles are not primary parameters – their values depend on many other parameters. Using these arguments, it will be demonstrated that contact angles are not the cause of wettability behavior but the effect of wettability behavior and other parameters. The result of this is that the contact angle cannot be used as a primary indicator of wettability except in very restricted situations. Furthermore, it will be demonstrated that even for the simple case of a capillary interface in a vertical tube, attempting to use simply a two-dimensional image to determine the contact angle can result in a wide range of measured values. This observation is consistent with some published experimental results. It follows that contact angles measured in two-dimensions cannot be trusted to provide accurate values and these values should not be used to characterize the wettability of the system.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Raju Bheemanahalli ◽  
Montana Knight ◽  
Cherryl Quinones ◽  
Colleen J. Doherty ◽  
S. V. Krishna Jagadish

AbstractHigh night temperatures (HNT) are shown to significantly reduce rice (Oryza sativa L.) yield and quality. A better understanding of the genetic architecture of HNT tolerance will help rice breeders to develop varieties adapted to future warmer climates. In this study, a diverse indica rice panel displayed a wide range of phenotypic variability in yield and quality traits under control night (24 °C) and higher night (29 °C) temperatures. Genome-wide association analysis revealed 38 genetic loci associated across treatments (18 for control and 20 for HNT). Nineteen loci were detected with the relative changes in the traits between control and HNT. Positive phenotypic correlations and co-located genetic loci with previously cloned grain size genes revealed common genetic regulation between control and HNT, particularly grain size. Network-based predictive models prioritized 20 causal genes at the genetic loci based on known gene/s expression under HNT in rice. Our study provides important insights for future candidate gene validation and molecular marker development to enhance HNT tolerance in rice. Integrated physiological, genomic, and gene network-informed approaches indicate that the candidate genes for stay-green trait may be relevant to minimizing HNT-induced yield and quality losses during grain filling in rice by optimizing source-sink relationships.


Sign in / Sign up

Export Citation Format

Share Document