Bayesian modeling of recombination events in bacterial populations

Background: We consider the discovery of recombinant segments jointly with their origins within multilocus DNA sequences from bacteria representing heterogeneous populations of fairly closely related species. The currently available methods for recombination detection capable of probabilistic characterization of uncertainty have a limited applicability in practice as the number of strains in a data set increases. Results: We introduce a Bayesian spatial structural model representing the continuum of origins over sites within the observed sequences, including a probabilistic characterization of uncertainty related to the origin of any particular site. To enable a statistically accurate and practically feasible approach to the analysis of large-scale data sets representing a single genus, we have developed a novel software tool (BRAT, Bayesian Recombination Tracker) implementing the model and the corresponding learning algorithm, which is capable of identifying the posterior optimal structure and to estimate the marginal posterior probabilities of putative origins over the sites. Conclusion: A multitude of challenging simulation scenarios and an analysis of real data from seven housekeeping genes of 120 strains of genus Burkholderia are used to illustrate the possibilities offered by our approach. The software is freely available for download at URL http://web.abo.fi/fak/ mnf//mate/jc/software/brat.html

Quantifying Adaptive Evolution in the Drosophila Immune System

It is estimated that a large proportion of amino acid substitutions in Drosophila have been fixed by natural selection, and as organisms are faced with an ever-changing array of pathogens and parasites to which they must adapt, we have investigated the role of parasite-mediated selection as a likely cause. To quantify the effect, and to identify which genes and pathways are most likely to be involved in the host–parasite arms race, we have re-sequenced population samples of 136 immunity and 287 position-matched non-immunity genes in two species of Drosophila. Using these data, and a new extension of the McDonald-Kreitman approach, we estimate that natural selection fixes advantageous amino acid changes in immunity genes at nearly double the rate of other genes. We find the rate of adaptive evolution in immunity genes is also more variable than other genes, with a small subset of immune genes evolving under intense selection. These genes, which are likely to represent hotspots of host–parasite coevolution, tend to share similar functions or belong to the same pathways, such as the antiviral RNAi pathway and the IMD signalling pathway. These patterns appear to be general features of immune system evolution in both species, as rates of adaptive evolution are correlated between the D. melanogaster and D. simulans lineages. In summary, our data provide quantitative estimates of the elevated rate of adaptive evolution in immune system genes relative to the rest of the genome, and they suggest that adaptation to parasites is an important force driving molecular evolution

Ultra-violet radiation is responsible for the differences in global epidemiology of chickenpox and the evolution of varicella-zoster virus as man migrated out of Africa

<p>Abstract</p> <p>Background</p> <p>Of the eight human herpes viruses, varicella-zoster virus, which causes chickenpox and zoster, has a unique epidemiology. Primary infection is much less common in children in the tropics compared with temperate areas. This results in increased adult susceptibility causing outbreaks, for example in health-care workers migrating from tropical to temperate countries. The recent demonstration that there are different genotypes of varicella-zoster virus and their geographic segregation into tropical and temperate areas suggests a distinct, yet previously unconsidered climatic factor may be responsible for both the clinical and molecular epidemiological features of this virus infection.</p> <p>Presentation of the hypothesis</p> <p>Unlike other human herpes viruses, varicella-zoster virus does not require intimate contact for infection to occur indicating that transmission may be interrupted by a geographically restricted climatic factor. The factor with the largest difference between tropical and temperate zones is ultra-violet radiation. This could reduce the infectiousness of chickenpox cases by inactivating virus in vesicles, before or after rupture. This would explain decreased transmissibility in the tropics and why the peak chickenpox incidence in temperate zones occurs during winter and spring, when ultra-violet radiation is at its lowest. The evolution of geographically restricted genotypes is also explained by ultra-violet radiation driving natural selection of different virus genotypes with varying degrees of resistance to inactivation, tropical genotypes being the most resistant. Consequently, temperate viruses should be more sensitive to its effects. This is supported by the observation that temperate genotypes are found in the tropics only in specific circumstances, namely where ultra-violet radiation has either been excluded or significantly reduced in intensity.</p> <p>Testing the Hypothesis</p> <p>The hypothesis is testable by exposing different virus genotypes to ultra-violet radiation and quantifying virus survival by plaque forming units or quantitative mRNA RT-PCR.</p> <p>Implications of the hypothesis</p> <p>The ancestral varicella-zoster virus, most probably a tropical genotype, co-migrated with man as he left Africa approximately 200,000 years ago. For this virus to have lost the selective advantage of resistance to ultra-violet radiation, the hypothesis would predict that the temperate, ultra-violet sensitive virus should have acquired another selective advantage as an evolutionary trade-off. One obvious advantage could be an increased reactivation rate as zoster to set up more rounds of chickenpox transmission. If this were so, the mechanism responsible for resistance to ultra-violet radiation might also be involved in reactivation and latency. This could then provide the first insight into a genetic correlate of the survival strategy of this virus.</p

Obscured phylogeny and possible recombinational dormancy in Escherichia coli

<p>Abstract</p> <p>Background</p> <p><it>Escherichia coli </it>is one of the best studied organisms in all of biology, but its phylogenetic structure has been difficult to resolve with current data and analytical techniques. We analyzed single nucleotide polymorphisms in chromosomes of representative strains to reconstruct the topology of its emergence.</p> <p>Results</p> <p>The phylogeny of <it>E. coli </it>varies according to the segment of chromosome analyzed. Recombination between extant <it>E. coli </it>groups is largely limited to only three intergroup pairings.</p> <p>Conclusions</p> <p>Segment-dependent phylogenies most likely are legacies of a complex recombination history. However, <it>E. coli </it>are now in an epoch in which they no longer broadly share DNA. Using the definition of species as organisms that freely exchange genetic material, this recombinational dormancy could reflect either the end of <it>E. coli </it>as a species, or herald the coalescence of <it>E. coli </it>groups into new species.</p

Tyrannobdella rex N. Gen. N. Sp. and the Evolutionary Origins of Mucosal Leech Infestations

BACKGROUND: Leeches have gained a fearsome reputation by feeding externally on blood, often from human hosts. Orificial hirudiniasis is a condition in which a leech enters a body orifice, most often the nasopharyngeal region, but there are many cases of leeches infesting the eyes, urethra, vagina, or rectum. Several leech species particularly in Africa and Asia are well-known for their propensity to afflict humans. Because there has not previously been any data suggesting a close relationship for such geographically disparate species, this unnerving tendency to be invasive has been regarded only as a loathsome oddity and not a unifying character for a group of related organisms. PRINCIPAL FINDINGS: A new genus and species of leech from Perú was found feeding from the nasopharynx of humans. Unlike any other leech previously described, this new taxon has but a single jaw with very large teeth. Phylogenetic analyses of nuclear and mitochondrial genes using parsimony and Bayesian inference demonstrate that the new species belongs among a larger, global clade of leeches, all of which feed from the mucosal surfaces of mammals. CONCLUSIONS: This new species, found feeding from the upper respiratory tract of humans in Perú, clarifies an expansion of the family Praobdellidae to include the new species Tyrannobdella rex n. gen. n. sp., along with others in the genera Dinobdella, Myxobdella, Praobdella and Pintobdella. Moreover, the results clarify a single evolutionary origin of a group of leeches that specializes on mucous membranes, thus, posing a distinct threat to human health

Observation of associated near-side and away-side long-range correlations in √sNN=5.02 TeV proton-lead collisions with the ATLAS detector

Two-particle correlations in relative azimuthal angle (Δϕ) and pseudorapidity (Δη) are measured in √sNN=5.02  TeV p+Pb collisions using the ATLAS detector at the LHC. The measurements are performed using approximately 1  μb-1 of data as a function of transverse momentum (pT) and the transverse energy (ΣETPb) summed over 3.1<η<4.9 in the direction of the Pb beam. The correlation function, constructed from charged particles, exhibits a long-range (2<|Δη|<5) “near-side” (Δϕ∼0) correlation that grows rapidly with increasing ΣETPb. A long-range “away-side” (Δϕ∼π) correlation, obtained by subtracting the expected contributions from recoiling dijets and other sources estimated using events with small ΣETPb, is found to match the near-side correlation in magnitude, shape (in Δη and Δϕ) and ΣETPb dependence. The resultant Δϕ correlation is approximately symmetric about π/2, and is consistent with a dominant cos⁡2Δϕ modulation for all ΣETPb ranges and particle pT

Search for direct pair production of the top squark in all-hadronic final states in proton-proton collisions at s√=8 TeV with the ATLAS detector

The results of a search for direct pair production of the scalar partner to the top quark using an integrated luminosity of 20.1fb−1 of proton–proton collision data at √s = 8 TeV recorded with the ATLAS detector at the LHC are reported. The top squark is assumed to decay via t˜→tχ˜01 or t˜→ bχ˜±1 →bW(∗)χ˜01 , where χ˜01 (χ˜±1 ) denotes the lightest neutralino (chargino) in supersymmetric models. The search targets a fully-hadronic final state in events with four or more jets and large missing transverse momentum. No significant excess over the Standard Model background prediction is observed, and exclusion limits are reported in terms of the top squark and neutralino masses and as a function of the branching fraction of t˜ → tχ˜01 . For a branching fraction of 100%, top squark masses in the range 270–645 GeV are excluded for χ˜01 masses below 30 GeV. For a branching fraction of 50% to either t˜ → tχ˜01 or t˜ → bχ˜±1 , and assuming the χ˜±1 mass to be twice the χ˜01 mass, top squark masses in the range 250–550 GeV are excluded for χ˜01 masses below 60 GeV