Genomic heterogeneity of historical gene flow between two species of newts inferred from transcriptome data

The role of gene flow in species formation is a major unresolved issue in speciation biology. Progress in this area requires information on the long‐term patterns of gene flow between diverging species. Here, we used thousands of single‐nucleotide polymorphisms derived from transcriptome resequencing and a method modeling the joint frequency spectrum of these polymorphisms to reconstruct patterns of historical gene flow between two Lissotriton newts: L. vulgaris (Lv) and L. montandoni (Lm). We tested several models of divergence including complete isolation and various scenarios of historical gene flow. The model of secondary contact received the highest support. According to this model, the species split from their common ancestor ca. 5.5 million years (MY) ago, evolved in isolation for ca. 2 MY, and have been exchanging genes for the last 3.5 MY Demographic changes have been inferred in both species, with the current effective population size of ca. 0.7 million in Lv and 0.2 million in Lm. The postdivergence gene flow resulted in two‐directional introgression which affected the genomes of both species, but was more pronounced from Lv to Lm. Interestingly, we found evidence for genomic heterogeneity of interspecific gene flow. This study demonstrates the complexity of long‐term gene flow between distinct but incompletely reproductively isolated taxa which divergence was initiated millions of years ago

Evolutionary game dynamics in phenotype space

Evolutionary dynamics can be studied in well-mixed or structured populations. Population structure typically arises from the heterogeneous distribution of individuals in physical space or on social networks. Here we introduce a new type of space to evolutionary game dynamics: phenotype space. The population is well-mixed in the sense that everyone is equally likely to interact with everyone else, but the behavioral strategies depend on distance in phenotype space. Individuals might behave differently towards those who look similar or dissimilar. Individuals mutate to nearby phenotypes. We study the `phenotypic space walk' of populations. We present analytic calculations that bring together ideas from coalescence theory and evolutionary game dynamics. As a particular example, we investigate the evolution of cooperation in phenotype space. We obtain a precise condition for natural selection to favor cooperators over defectors: for a one-dimensional phenotype space and large population size the critical benefit-to-cost ratio is given by b/c=1+2/sqrt{3}. We derive the fundamental condition for any evolutionary game and explore higher dimensional phenotype spaces.Comment: version 2: minor changes; equivalent to final published versio

Vulnerability to depression is associated with a failure to acquire implicit social appraisals

Major depressive disorder is frequently associated with disrupted relationships with spouses, partners, family and peers. These problems can precipitate the onset of clinical illness, influence severity and the prospects for treatment and recovery. Here, we investigated whether individuals who have recovered from depression use interpersonal signals to form favourable appraisals of others as social partners. Twenty recovered-depressed adults (with at least two adult episodes of major depressive disorder but euthymic and medication-free for six months) and twenty three healthy, never-depressed adults completed a reaction time task in which the gaze direction of some faces reliably cued the location a target (valid faces), whereas the gaze direction of other faces cued the opposite spatial location (invalid faces). None of the participants were aware of this contingency. Following this task, participants judged the trustworthiness of the faces. Both the recovered-depressed and healthy never-depressed participants were significantly faster to categorise targets following valid compared with invalid gaze cueing faces. Whereas the healthy never-depressed participants judged the valid faces to be significantly more trustworthy than the invalid faces; this implicit social appraisal was absent in the recovered-depressed participants. Individuals who have recovered from major depressive disorder are able to respond appropriately to joint attention with other people but appear to not use joint attention to form implicit trust appraisals of others as potential social partners

Ordering in voter models on networks: Exact reduction to a single-coordinate diffusion

We study the voter model and related random-copying processes on arbitrarily complex network structures. Through a representation of the dynamics as a particle reaction process, we show that a quantity measuring the degree of order in a finite system is, under certain conditions, exactly governed by a universal diffusion equation. Whenever this reduction occurs, the details of the network structure and random-copying process affect only a single parameter in the diffusion equation. The validity of the reduction can be established with considerably less information than one might expect: it suffices to know just two characteristic timescales within the dynamics of a single pair of reacting particles. We develop methods to identify these timescales, and apply them to deterministic and random network structures. We focus in particular on how the ordering time is affected by degree correlations, since such effects are hard to access by existing theoretical approaches.Comment: 37 pages, 10 figures. Revised version with additional discussion and simulation results to appear in J Phys

Search for the Proton Decay Mode proton to neutrino K+ in Soudan 2

We have searched for the proton decay mode proton to neutrino K+ using the one-kiloton Soudan 2 high resolution calorimeter. Contained events obtained from a 3.56 kiloton-year fiducial exposure through June 1997 are examined for occurrence of a visible K+ track which decays at rest into mu+ nu or pi+ pi0. We found one candidate event consistent with background, yielding a limit, tau/B > 4.3 10^{31} years at 90% CL with no background subtraction.Comment: 13 pages, Latex, 3 tables and 3 figures, Accepted by Physics Letters

Evaluation of a candidate breast cancer associated SNP in ERCC4 as a risk modifier in BRCA1 and BRCA2 mutation carriers. Results from the Consortium of Investigators of Modifiers of BRCA1/BRCA2 (CIMBA)

Background: In this study we aimed to evaluate the role of a SNP in intron 1 of the ERCC4 gene (rs744154), previously reported to be associated with a reduced risk of breast cancer in the general population, as a breast cancer risk modifier in BRCA1 and BRCA2 mutation carriers. Methods: We have genotyped rs744154 in 9408 BRCA1 and 5632 BRCA2 mutation carriers from the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA) and assessed its association with breast cancer risk using a retrospective weighted cohort approach. Results: We found no evidence of association with breast cancer risk for BRCA1 (per-allele HR: 0.98, 95% CI: 0.93–1.04, P=0.5) or BRCA2 (per-allele HR: 0.97, 95% CI: 0.89–1.06, P=0.5) mutation carriers. Conclusion: This SNP is not a significant modifier of breast cancer risk for mutation carriers, though weak associations cannot be ruled out. A Osorio1, R L Milne2, G Pita3, P Peterlongo4,5, T Heikkinen6, J Simard7, G Chenevix-Trench8, A B Spurdle8, J Beesley8, X Chen8, S Healey8, KConFab9, S L Neuhausen10, Y C Ding10, F J Couch11,12, X Wang11, N Lindor13, S Manoukian4, M Barile14, A Viel15, L Tizzoni5,16, C I Szabo17, L Foretova18, M Zikan19, K Claes20, M H Greene21, P Mai21, G Rennert22, F Lejbkowicz22, O Barnett-Griness22, I L Andrulis23,24, H Ozcelik24, N Weerasooriya23, OCGN23, A-M Gerdes25, M Thomassen25, D G Cruger26, M A Caligo27, E Friedman28,29, B Kaufman28,29, Y Laitman28, S Cohen28, T Kontorovich28, R Gershoni-Baruch30, E Dagan31,32, H Jernström33, M S Askmalm34, B Arver35, B Malmer36, SWE-BRCA37, S M Domchek38, K L Nathanson38, J Brunet39, T Ramón y Cajal40, D Yannoukakos41, U Hamann42, HEBON37, F B L Hogervorst43, S Verhoef43, EB Gómez García44,45, J T Wijnen46,47, A van den Ouweland48, EMBRACE37, D F Easton49, S Peock49, M Cook49, C T Oliver49, D Frost49, C Luccarini50, D G Evans51, F Lalloo51, R Eeles52, G Pichert53, J Cook54, S Hodgson55, P J Morrison56, F Douglas57, A K Godwin58, GEMO59,60,61, O M Sinilnikova59,60, L Barjhoux59,60, D Stoppa-Lyonnet61, V Moncoutier61, S Giraud59, C Cassini62,63, L Olivier-Faivre62,63, F Révillion64, J-P Peyrat64, D Muller65, J-P Fricker65, H T Lynch66, E M John67, S Buys68, M Daly69, J L Hopper70, M B Terry71, A Miron72, Y Yassin72, D Goldgar73, Breast Cancer Family Registry37, C F Singer74, D Gschwantler-Kaulich74, G Pfeiler74, A-C Spiess74, Thomas v O Hansen75, O T Johannsson76, T Kirchhoff77, K Offit77, K Kosarin77, M Piedmonte78, G C Rodriguez79, K Wakeley80, J F Boggess81, J Basil82, P E Schwartz83, S V Blank84, A E Toland85, M Montagna86, C Casella87, E N Imyanitov88, A Allavena89, R K Schmutzler90, B Versmold90, C Engel91, A Meindl92, N Ditsch93, N Arnold94, D Niederacher95, H Deißler96, B Fiebig97, R Varon-Mateeva98, D Schaefer99, U G Froster100, T Caldes101, M de la Hoya101, L McGuffog49, A C Antoniou49, H Nevanlinna6, P Radice4,5 and J Benítez1,3 on behalf of CIMB

Estimating time since infection in early homogeneous HIV-1 samples using a poisson model

<p>Abstract</p> <p>Background</p> <p>The occurrence of a genetic bottleneck in HIV sexual or mother-to-infant transmission has been well documented. This results in a majority of new infections being homogeneous, <it>i.e</it>., initiated by a single genetic strain. Early after infection, prior to the onset of the host immune response, the viral population grows exponentially. In this simple setting, an approach for estimating evolutionary and demographic parameters based on comparison of diversity measures is a feasible alternative to the existing Bayesian methods (<it>e.g</it>., BEAST), which are instead based on the simulation of genealogies.</p> <p>Results</p> <p>We have devised a web tool that analyzes genetic diversity in acutely infected HIV-1 patients by comparing it to a model of neutral growth. More specifically, we consider a homogeneous infection (<it>i.e</it>., initiated by a unique genetic strain) prior to the onset of host-induced selection, where we can assume a random accumulation of mutations. Previously, we have shown that such a model successfully describes about 80% of sexual HIV-1 transmissions provided the samples are drawn early enough in the infection. Violation of the model is an indicator of either heterogeneous infections or the initiation of selection.</p> <p>Conclusions</p> <p>When the underlying assumptions of our model (homogeneous infection prior to selection and fast exponential growth) are met, we are under a very particular scenario for which we can use a forward approach (instead of backwards in time as provided by coalescent methods). This allows for more computationally efficient methods to derive the time since the most recent common ancestor. Furthermore, the tool performs statistical tests on the Hamming distance frequency distribution, and outputs summary statistics (mean of the best fitting Poisson distribution, goodness of fit p-value, etc). The tool runs within minutes and can readily accommodate the tens of thousands of sequences generated through new ultradeep pyrosequencing technologies. The tool is available on the LANL website.</p

Genetic Demixing and Evolutionary Forces in the One-Dimensional Stepping Stone Model

We review and extend results for mutation, selection, genetic drift, and migration in a one-dimensional continuous population. The population is described by a continuous limit of the stepping stone model, which leads to the stochastic Fisher-Kolmogorov-Petrovsky-Piscounov equation with additional terms describing mutations. Although the stepping stone model was first proposed for population genetics, it is closely related to "voter models" of interest in nonequilibrium statistical mechanics. The stepping stone model can also be regarded as an approximation to the dynamics of a thin layer of actively growing pioneers at the frontier of a colony of microorganisms undergoing a range expansion on a Petri dish. We find that the population tends to segregate into monoallelic domains. This segregation slows down genetic drift and selection because these two evolutionary forces can only act at the boundaries between the domains; the effects of mutation, however, are not significantly affected by the segregation. Although fixation in the neutral well-mixed (or "zero dimensional") model occurs exponentially in time, it occurs only algebraically fast in the one-dimensional model. If selection is weak, selective sweeps occur exponentially fast in both well-mixed and one-dimensional populations, but the time constants are different. We also find an unusual sublinear increase in the variance of the spatially averaged allele frequency with time. Although we focus on two alleles or variants, q-allele Potts-like models of gene segregation are considered as well. We also investigate the effects of geometry at the frontier by considering growth of circular colonies. Our analytical results are checked with simulations, and could be tested against recent spatial experiments on range expansions off linear inoculations of Escherichia coli and Saccharomyces cerevisiae.Comment: 29 pages, 20 figures; Reviews of Modern Physics, Volume 82, April-June 201

The Genomics of Speciation in Drosophila: Diversity, Divergence, and Introgression Estimated Using Low-Coverage Genome Sequencing

In nature, closely related species may hybridize while still retaining their distinctive identities. Chromosomal regions that experience reduced recombination in hybrids, such as within inversions, have been hypothesized to contribute to the maintenance of species integrity. Here, we examine genomic sequences from closely related fruit fly taxa of the Drosophila pseudoobscura subgroup to reconstruct their evolutionary histories and past patterns of genic exchange. Partial genomic assemblies were generated from two subspecies of Drosophila pseudoobscura (D. ps.) and an outgroup species, D. miranda. These new assemblies were compared to available assemblies of D. ps. pseudoobscura and D. persimilis, two species with overlapping ranges in western North America. Within inverted regions, nucleotide divergence among each pair of the three species is comparable, whereas divergence between D. ps. pseudoobscura and D. persimilis in non-inverted regions is much lower and closer to levels of intraspecific variation. Using molecular markers flanking each of the major chromosomal inversions, we identify strong crossover suppression in F1 hybrids extending over 2 megabase pairs (Mbp) beyond the inversion breakpoints. These regions of crossover suppression also exhibit the high nucleotide divergence associated with inverted regions. Finally, by comparison to a geographically isolated subspecies, D. ps. bogotana, our results suggest that autosomal gene exchange between the North American species, D. ps. pseudoobscura and D. persimilis, occurred since the split of the subspecies, likely within the last 200,000 years. We conclude that chromosomal rearrangements have been vital to the ongoing persistence of these species despite recent hybridization. Our study serves as a proof-of-principle on how whole genome sequencing can be applied to formulate and test hypotheses about species formation in lesser-known non-model systems