Phylogenetic analysis of Tomato spotted wilt virus (TSWV) NSs protein demonstrates the isolated emergence of resistance-breaking strains in pepper

Resurgence of Tomato spotted wilt virus (TSWV) worldwide as well as in Hungary causing heavy economic losses directed the attention to the factors con- tributing to the outbreak of this serious epidemics. The introgression of Tsw resistance gene into various pepper cultivars seemed to solve TSWV control, but widely used resistant pepper cultivars bearing the same, unique resis- tance locus evoked the rapid emergence of resistance- breaking (RB) TSWV strains. In Hungary, the sporadic appearance of RB strains in pepper-producing region was first observed in 2010–2011, but in 2012 it was detected frequently. Previously, the non-structural protein (NSs) encoded by small RNA (S RNA) of TSWV was verified as the avirulence factor for Tsw resistance, therefore we analyzed the S RNA of the Hungarian RB and wild type (WT) isolates and compared to previously analyzed TSWV strains with RB properties from different geographical origins. Phylogenetic analysis demonstrated that the dif- ferent RB strains had the closest relationship with the local WT isolates and there is no conserved mutation present in all the NSs genes of RB isolates from different geograph- ical origins. According to these results, we concluded that the RB isolates evolved separately in geographic point of view, and also according to the RB mechanism

Investigating mitochondrial DNA relationships in Neolithic Western Europe through serial coalescent simulations

Recent ancient DNA studies on European Neolithic human populations have provided persuasive evidence of a major migration of farmers originating from the Aegean, accompanied by sporadic hunter-gatherer admixture into early Neolithic populations, but increasing toward the Late Neolithic. In this context, ancient mitochondrial DNA data collected from the Neolithic necropolis of Gurgy (Paris Basin, France), the largest mitochondrial DNA sample obtained from a single archeological site for the Early/Middle Neolithic period, indicate little differentiation from farmers associated to both the Danubian and Mediterranean Neolithic migration routes, as well as from Western European hunter-gatherers. To test whether this pattern of differentiation could arise in a single unstructured population by genetic drift alone, we used serial coalescent simulations. We explore female effective population size parameter combinations at the time of the colonization of Europe 45000 years ago and the most recent of the Neolithic samples analyzed in this study 5900 years ago, and identify conditions under which population panmixia between hunter-gatherers/Early-Middle Neolithic farmers and Gurgy cannot be rejected. In relation to other studies on the current debate of the origins of Europeans, these results suggest increasing hunter-gatherer admixture into farmers' group migrating farther west in Europe.European Journal of Human Genetics advance online publication, 28 December 2016; doi:10.1038/ejhg.2016.180

The genetic legacy of extreme exploitation in a polar vertebrate

Understanding the effects of human exploitation on the genetic composition of wild populations is important for predicting species persistence and adaptive potential. We therefore investigated the genetic legacy of large-scale commercial harvesting by reconstructing, on a global scale, the recent demographic history of the Antarctic fur seal (Arctocephalus gazella), a species that was hunted to the brink of extinction by 18th and 19th century sealers. Molecular genetic data from over 2,000 individuals sampled from all eight major breeding locations across the species’ circumpolar geographic distribution, show that at least four relict populations around Antarctica survived commercial hunting. Coalescent simulations suggest that all of these populations experienced severe bottlenecks down to effective population sizes of around 150–200. Nevertheless, comparably high levels of neutral genetic variability were retained as these declines are unlikely to have been strong enough to deplete allelic richness by more than around 15%. These findings suggest that even dramatic short-term declines need not necessarily result in major losses of diversity, and explain the apparent contradiction between the high genetic diversity of this species and its extreme exploitation history

Conservation of pattern as a tool for inference on spatial snapshots in ecological data

As climate change and other anthropogenic factors increase the uncertainty of vegetation ecosystem persistence, the ability to rapidly assess their dynamics is paramount. Vegetation and sessile communities form a variety of striking regular spatial patterns such as stripes, spots and labyrinths, that have been used as indicators of ecosystem current state, through qualitative analysis of simple models. Here we describe a new method for rigorous quantitative estimation of biological parameters from a single spatial snapshot. We formulate a synthetic likelihood through consideration of the expected change in the correlation structure of the spatial pattern. This then allows Bayesian inference to be performed on the model parameters, which includes providing parameter uncertainty. The method was validated against simulated data and then applied to real data in the form of aerial photographs of seagrass banding. The inferred parameters were found to be able to reproduce similar patterns to those observed and able to detect strength of spatial competition, competition-induced mortality and the local range of reproduction. This technique points to a way of performing rapid inference of spatial competition and ecological stability from a single spatial snapshots of sessile communities

A GWAS in Latin Americans highlights the convergent evolution of lighter skin pigmentation in Eurasia

We report a genome-wide association scan in >6,000 Latin Americans for pigmentation of skin and eyes. We found eighteen signals of association at twelve genomic regions. These include one novel locus for skin pigmentation (in 10q26) and three novel loci for eye pigmentation (in 1q32, 20q13 and 22q12). We demonstrate the presence of multiple independent signals of association in the 11q14 and 15q13 regions (comprising the GRM5/TYR and HERC2/OCA2 genes, respectively) and several epistatic interactions among independently associated alleles. Strongest association with skin pigmentation at 19p13 was observed for an Y182H missense variant (common only in East Asians and Native Americans) in MFSD12, a gene recently associated with skin pigmentation in Africans. We show that the frequency of the derived allele at Y182H is significantly correlated with lower solar radiation intensity in East Asia and infer that MFSD12 was under selection in East Asians, probably after their split from Europeans

Importance of incomplete lineage sorting and introgression in the origin of shared genetic variation between two closely related pines with overlapping distributions

Genetic variation shared between closely related species may be due to retention of ancestral polymorphisms because of incomplete lineage sorting (ILS) and/or introgression following secondary contact. It is challenging to distinguish ILS and introgression because they generate similar patterns of shared genetic diversity, but this is nonetheless essential for inferring accurately the history of species with overlapping distributions. To address this issue, we sequenced 33 independent intron loci across the genome of two closely related pine species (Pinus massoniana Lamb. and Pinus hwangshanensis Hisa) from Southeast China. Population structure analyses revealed that the species showed slightly more admixture in parapatric populations than in allopatric populations. Levels of interspecific differentiation were lower in parapatry than in allopatry. Approximate Bayesian computation suggested that the most likely speciation scenario explaining this pattern was a long period of isolation followed by a secondary contact. Ecological niche modeling suggested that a gradual range expansion of P. hwangshanensis during the Pleistocene climatic oscillations could have been the cause of the overlap. Our study therefore suggests that secondary introgression, rather than ILS, explains most of the shared nuclear genomic variation between these two species and demonstrates the complementarity of population genetics and ecological niche modeling in understanding gene flow history. Finally, we discuss the importance of contrasting results from markers with different dynamics of migration, namely nuclear, chloroplast and mitochondrial DNA

Modelling the impact of larviciding on the population dynamics and biting rates of Simulium damnosum (s.l.): implications for vector control as a complementary strategy for onchocerciasis elimination in Africa

Background: In 2012, the World Health Organization set goals for the elimination of onchocerciasis transmission by 2020 in selected African countries. Epidemiological data and mathematical modelling have indicated that elimination may not be achieved with annual ivermectin distribution in all endemic foci. Complementary and alternative treatment strategies (ATS), including vector control, will be necessary. Implementation of vector control will require that the ecology and population dynamics of Simulium damnosum sensu lato be carefully considered. Methods: We adapted our previous SIMuliid POPulation dynamics (SIMPOP) model to explore the impact of larvicidal insecticides on S. damnosum (s.l.) biting rates in different ecological contexts and to identify how frequently and for how long vector control should be continued to sustain substantive reductions in vector biting. SIMPOP was fitted to data from large-scale aerial larviciding trials in savannah sites (Ghana) and small-scale ground larviciding trials in forest areas (Cameroon). The model was validated against independent data from Burkina Faso/Côte d’Ivoire (savannah) and Bioko (forest). Scenario analysis explored the effects of ecological and programmatic factors such as pre-control daily biting rate (DBR) and larviciding scheme design on reductions and resurgences in biting rates. Results: The estimated efficacy of large-scale aerial larviciding in the savannah was greater than that of ground-based larviciding in the forest. Small changes in larvicidal efficacy can have large impacts on intervention success. At 93% larvicidal efficacy (a realistic value based on field trials), 10 consecutive weekly larvicidal treatments would reduce DBRs by 96% (e.g. from 400 to 16 bites/person/day). At 70% efficacy, and for 10 weekly applications, the DBR would decrease by 67% (e.g. from 400 to 132 bites/person/day). Larviciding is more likely to succeed in areas with lower water temperatures and where blackfly species have longer gonotrophic cycles. Conclusions: Focal vector control can reduce vector biting rates in settings where a high larvicidal efficacy can be achieved and an appropriate duration and frequency of larviciding can be ensured. Future work linking SIMPOP with onchocerciasis transmission models will permit evaluation of the impact of combined anti-vectorial and anti-parasitic interventions on accelerating elimination of the disease

A framework for parameter estimation and model selection from experimental data in systems biology using approximate Bayesian computation.

As modeling becomes a more widespread practice in the life sciences and biomedical sciences, researchers need reliable tools to calibrate models against ever more complex and detailed data. Here we present an approximate Bayesian computation (ABC) framework and software environment, ABC-SysBio, which is a Python package that runs on Linux and Mac OS X systems and that enables parameter estimation and model selection in the Bayesian formalism by using sequential Monte Carlo (SMC) approaches. We outline the underlying rationale, discuss the computational and practical issues and provide detailed guidance as to how the important tasks of parameter inference and model selection can be performed in practice. Unlike other available packages, ABC-SysBio is highly suited for investigating, in particular, the challenging problem of fitting stochastic models to data. In order to demonstrate the use of ABC-SysBio, in this protocol we postulate the existence of an imaginary reaction network composed of seven interrelated biological reactions (involving a specific mRNA, the protein it encodes and a post-translationally modified version of the protein), a network that is defined by two files containing 'observed' data that we provide as supplementary information. In the first part of the PROCEDURE, ABC-SysBio is used to infer the parameters of this system, whereas in the second part we use ABC-SysBio's relevant functionality to discriminate between two different reaction network models, one of them being the 'true' one. Although computationally expensive, the additional insights gained in the Bayesian formalism more than make up for this cost, especially in complex problems