Who Watches the Watchmen? An Appraisal of Benchmarks for Multiple Sequence Alignment

Multiple sequence alignment (MSA) is a fundamental and ubiquitous technique in bioinformatics used to infer related residues among biological sequences. Thus alignment accuracy is crucial to a vast range of analyses, often in ways difficult to assess in those analyses. To compare the performance of different aligners and help detect systematic errors in alignments, a number of benchmarking strategies have been pursued. Here we present an overview of the main strategies--based on simulation, consistency, protein structure, and phylogeny--and discuss their different advantages and associated risks. We outline a set of desirable characteristics for effective benchmarking, and evaluate each strategy in light of them. We conclude that there is currently no universally applicable means of benchmarking MSA, and that developers and users of alignment tools should base their choice of benchmark depending on the context of application--with a keen awareness of the assumptions underlying each benchmarking strategy.Comment: Revie

Partitioning of on-demand electron pairs

We demonstrate the high fidelity splitting of electron pairs emitted on demand from a dynamic quantum dot by an electronic beam splitter. The fidelity of pair splitting is inferred from the coincidence of arrival in two detector paths probed by a measurement of the partitioning noise. The emission characteristic of the on-demand electron source is tunable from electrons being partitioned equally and independently to electron pairs being split with a fidelity of 90%. For low beam splitter transmittance we further find evidence of pair bunching violating statistical expectations for independent fermions

Evolutionary distances in the twilight zone -- a rational kernel approach

Phylogenetic tree reconstruction is traditionally based on multiple sequence alignments (MSAs) and heavily depends on the validity of this information bottleneck. With increasing sequence divergence, the quality of MSAs decays quickly. Alignment-free methods, on the other hand, are based on abstract string comparisons and avoid potential alignment problems. However, in general they are not biologically motivated and ignore our knowledge about the evolution of sequences. Thus, it is still a major open question how to define an evolutionary distance metric between divergent sequences that makes use of indel information and known substitution models without the need for a multiple alignment. Here we propose a new evolutionary distance metric to close this gap. It uses finite-state transducers to create a biologically motivated similarity score which models substitutions and indels, and does not depend on a multiple sequence alignment. The sequence similarity score is defined in analogy to pairwise alignments and additionally has the positive semi-definite property. We describe its derivation and show in simulation studies and real-world examples that it is more accurate in reconstructing phylogenies than competing methods. The result is a new and accurate way of determining evolutionary distances in and beyond the twilight zone of sequence alignments that is suitable for large datasets.Comment: to appear in PLoS ON

Arterial oxygen content is precisely maintained by graded erythrocytotic responses in settings of high/normal serum iron levels, and predicts exercise capacity: an observational study of hypoxaemic patients with pulmonary arteriovenous malformations.

Oxygen, haemoglobin and cardiac output are integrated components of oxygen transport: each gram of haemoglobin transports 1.34 mls of oxygen in the blood. Low arterial partial pressure of oxygen (PaO2), and haemoglobin saturation (SaO2), are the indices used in clinical assessments, and usually result from low inspired oxygen concentrations, or alveolar/airways disease. Our objective was to examine low blood oxygen/haemoglobin relationships in chronically compensated states without concurrent hypoxic pulmonary vasoreactivity.165 consecutive unselected patients with pulmonary arteriovenous malformations were studied, in 98 cases, pre/post embolisation treatment. 159 (96%) had hereditary haemorrhagic telangiectasia. Arterial oxygen content was calculated by SaO2 x haemoglobin x 1.34/100.There was wide variation in SaO2 on air (78.5-99, median 95)% but due to secondary erythrocytosis and resultant polycythaemia, SaO2 explained only 0.1% of the variance in arterial oxygen content per unit blood volume. Secondary erythrocytosis was achievable with low iron stores, but only if serum iron was high-normal: Low serum iron levels were associated with reduced haemoglobin per erythrocyte, and overall arterial oxygen content was lower in iron deficient patients (median 16.0 [IQR 14.9, 17.4]mls/dL compared to 18.8 [IQR 17.4, 20.1]mls/dL, p<0.0001). Exercise tolerance appeared unrelated to SaO2 but was significantly worse in patients with lower oxygen content (p<0.0001). A pre-defined athletic group had higher Hb:SaO2 and serum iron:ferritin ratios than non-athletes with normal exercise capacity. PAVM embolisation increased SaO2, but arterial oxygen content was precisely restored by a subsequent fall in haemoglobin: 86 (87.8%) patients reported no change in exercise tolerance at post-embolisation follow-up.Haemoglobin and oxygen measurements in isolation do not indicate the more physiologically relevant oxygen content per unit blood volume. This can be maintained for SaO2 ≥78.5%, and resets to the same arterial oxygen content after correction of hypoxaemia. Serum iron concentrations, not ferritin, seem to predict more successful polycythaemic responses

Cryo-EM structure of a helicase loading intermediate containing ORC-Cdc6-Cdt1-MCM2-7 bound to DNA

In eukaryotes, the Cdt1-bound replicative helicase core MCM2-7 is loaded onto DNA by the ORC-Cdc6 ATPase to form a prereplicative complex (pre-RC) with an MCM2-7 double hexamer encircling DNA. Using purified components in the presence of ATP-γS, we have captured in vitro an intermediate in pre-RC assembly that contains a complex between the ORC-Cdc6 and Cdt1-MCM2-7 heteroheptamers called the OCCM. Cryo-EM studies of this 14-subunit complex reveal that the two separate heptameric complexes are engaged extensively, with the ORC-Cdc6 N-terminal AAA+ domains latching onto the C-terminal AAA+ motor domains of the MCM2-7 hexamer. The conformation of ORC-Cdc6 undergoes a concerted change into a right-handed spiral with helical symmetry that is identical to that of the DNA double helix. The resulting ORC-Cdc6 helicase loader shows a notable structural similarity to the replication factor C clamp loader, suggesting a conserved mechanism of action

Time spent with cats is never wasted: Lessons learned from feline acromegalic cardiomyopathy, a naturally occurring animal model of the human disease

<div><p>Background</p><p>In humans, acromegaly due to a pituitary somatotrophic adenoma is a recognized cause of increased left ventricular (LV) mass. Acromegalic cardiomyopathy is incompletely understood, and represents a major cause of morbidity and mortality. We describe the clinical, echocardiographic and histopathologic features of naturally occurring feline acromegalic cardiomyopathy, an emerging disease among domestic cats.</p><p>Methods</p><p>Cats with confirmed hypersomatotropism (IGF-1>1000ng/ml and pituitary mass; n = 67) were prospectively recruited, as were two control groups: diabetics (IGF-1<800ng/ml; n = 24) and healthy cats without known endocrinopathy or cardiovascular disease (n = 16). Echocardiography was performed in all cases, including after hypersomatotropism treatment where applicable. Additionally, tissue samples from deceased cats with hypersomatotropism, hypertrophic cardiomyopathy and age-matched controls (n = 21 each) were collected and systematically histopathologically reviewed and compared.</p><p>Results</p><p>By echocardiography, cats with hypersomatotropism had a greater maximum LV wall thickness (6.5mm, 4.1–10.1mm) than diabetic (5.9mm, 4.2–9.1mm; Mann Whitney, p<0.001) or control cats (5.2mm, 4.1–6.5mm; Mann Whitney, p<0.001). Left atrial diameter was also greater in cats with hypersomatotropism (16.6mm, 13.0–29.5mm) than in diabetic (15.4mm, 11.2–20.3mm; Mann Whitney, p<0.001) and control cats (14.0mm, 12.6–17.4mm; Mann Whitney, p<0.001). After hypophysectomy and normalization of IGF-1 concentration (n = 20), echocardiographic changes proved mostly reversible. As in humans, histopathology of the feline acromegalic heart was dominated by myocyte hypertrophy with interstitial fibrosis and minimal myofiber disarray.</p><p>Conclusions</p><p>These results demonstrate cats could be considered a naturally occurring model of acromegalic cardiomyopathy, and as such help elucidate mechanisms driving cardiovascular remodeling in this disease.</p></div

A primary care, multi-disciplinary disease management program for opioid-treated patients with chronic non-cancer pain and a high burden of psychiatric comorbidity

BACKGROUND: Chronic non-cancer pain is a common problem that is often accompanied by psychiatric comorbidity and disability. The effectiveness of a multi-disciplinary pain management program was tested in a 3 month before and after trial. METHODS: Providers in an academic general medicine clinic referred patients with chronic non-cancer pain for participation in a program that combined the skills of internists, clinical pharmacists, and a psychiatrist. Patients were either receiving opioids or being considered for opioid therapy. The intervention consisted of structured clinical assessments, monthly follow-up, pain contracts, medication titration, and psychiatric consultation. Pain, mood, and function were assessed at baseline and 3 months using the Brief Pain Inventory (BPI), the Center for Epidemiological Studies-Depression Scale scale (CESD) and the Pain Disability Index (PDI). Patients were monitored for substance misuse. RESULTS: Eighty-five patients were enrolled. Mean age was 51 years, 60% were male, 78% were Caucasian, and 93% were receiving opioids. Baseline average pain was 6.5 on an 11 point scale. The average CESD score was 24.0, and the mean PDI score was 47.0. Sixty-three patients (73%) completed 3 month follow-up. Fifteen withdrew from the program after identification of substance misuse. Among those completing 3 month follow-up, the average pain score improved to 5.5 (p = 0.003). The mean PDI score improved to 39.3 (p < 0.001). Mean CESD score was reduced to 18.0 (p < 0.001), and the proportion of depressed patients fell from 79% to 54% (p = 0.003). Substance misuse was identified in 27 patients (32%). CONCLUSIONS: A primary care disease management program improved pain, depression, and disability scores over three months in a cohort of opioid-treated patients with chronic non-cancer pain. Substance misuse and depression were common, and many patients who had substance misuse identified left the program when they were no longer prescribed opioids. Effective care of patients with chronic pain should include rigorous assessment and treatment of these comorbid disorders and intensive efforts to insure follow up

Cooperation, Norms, and Revolutions: A Unified Game-Theoretical Approach

Cooperation is of utmost importance to society as a whole, but is often challenged by individual self-interests. While game theory has studied this problem extensively, there is little work on interactions within and across groups with different preferences or beliefs. Yet, people from different social or cultural backgrounds often meet and interact. This can yield conflict, since behavior that is considered cooperative by one population might be perceived as non-cooperative from the viewpoint of another. To understand the dynamics and outcome of the competitive interactions within and between groups, we study game-dynamical replicator equations for multiple populations with incompatible interests and different power (be this due to different population sizes, material resources, social capital, or other factors). These equations allow us to address various important questions: For example, can cooperation in the prisoner's dilemma be promoted, when two interacting groups have different preferences? Under what conditions can costly punishment, or other mechanisms, foster the evolution of norms? When does cooperation fail, leading to antagonistic behavior, conflict, or even revolutions? And what incentives are needed to reach peaceful agreements between groups with conflicting interests? Our detailed quantitative analysis reveals a large variety of interesting results, which are relevant for society, law and economics, and have implications for the evolution of language and culture as well

Structure and evolutionary origin of Ca2+-dependent herring type II antifreeze protein

10.1371/journal.pone.0000548PLoS ONE26

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