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

How reliably can we predict the reliability of protein structure predictions?

Background: Comparative methods have been the standard techniques for in silico protein structure prediction. The prediction is based on a multiple alignment that contains both reference sequences with known structures and the sequence whose unknown structure is predicted. Intensive research has been made to improve the quality of multiple alignments, since misaligned parts of the multiple alignment yield misleading predictions. However, sometimes all methods fail to predict the correct alignment, because the evolutionary signal is too weak to find the homologous parts due to the large number of mutations that separate the sequences. Results: Stochastic sequence alignment methods define a posterior distribution of possible multiple alignments. They can highlight the most likely alignment, and above that, they can give posterior probabilities for each alignment column. We made a comprehensive study on the HOMSTRAD database of structural alignments, predicting secondary structures in four different ways. We showed that alignment posterior probabilities correlate with the reliability of secondary structure predictions, though the strength of the correlation is different for different protocols. The correspondence between the reliability of secondary structure predictions and alignment posterior probabilities is the closest to the identity function when the secondary structure posterior probabilities are calculated from the posterior distribution of multiple alignments. The largest deviation from the identity function has been obtained in the case of predicting secondary structures from a single optimal pairwise alignment. We also showed that alignment posterior probabilities correlate with the 3D distances between C α amino acids in superimposed tertiary structures. Conclusion: Alignment posterior probabilities can be used to a priori detect errors in comparative models on the sequence alignment level. </p

A Generic System for the Expression and Purification of Soluble and Stable Influenza Neuraminidase

The influenza surface glycoprotein neuraminidase (NA) is essential for the efficient spread of the virus. Antiviral drugs such as Tamiflu (oseltamivir) and Relenza (zanamivir) that inhibit NA enzyme activity have been shown to be effective in the treatment of influenza infections. The recent ‘swine flu’ pandemic and world-wide emergence of Tamiflu-resistant seasonal human influenza A(H1N1) H274Y have highlighted the need for the ongoing development of new anti-virals, efficient production of vaccine proteins and novel diagnostic tools. Each of these goals could benefit from the production of large quantities of highly pure and stable NA. This publication describes a generic expression system for NAs in a baculovirus Expression Vector System (BEVS) that is capable of expressing milligram amounts of recombinant NA. To construct NAs with increased stability, the natural influenza NA stalk was replaced by two different artificial tetramerization domains that drive the formation of catalytically active NA homotetramers: GCN4-pLI from yeast or the Tetrabrachion tetramerization domain from Staphylothermus marinus. Both recombinant NAs are secreted as FLAG-tagged proteins to allow for rapid and simple purification. The Tetrabrachion-based NA showed good solubility, increased stability and biochemical properties closer to the original viral NA than the GCN4-pLI based construct. The expressed quantities and high quality of the purified recombinant NA suggest that this expression system is capable of producing recombinant NA for a broad range of applications including high-throughput drug screening, protein crystallisation, or vaccine development

Aeromonas spp.-mediated cell-contact cytotoxicity is associated with the presence of type III secretion system

In the study we examined the production of cytotonic and cytotoxic toxins and the presence of a type III secretion system (TTSS) in 64 Aeromonas spp. strains isolated from fecal specimens of patients with gastroenteritis. We observed that contact of the bacteria with host epithelial cells is a prerequisite for their cytotoxicity at 3 h incubation. Cell-contact cytotoxic activity of the strains was strongly associated with the presence of the TTSS. Culture supernatants of the strains induced low cytotoxicity effects at the same time of incubation. Cell-free supernatants of 61 (95%) isolates expressed cytotoxic activity which caused the destruction of HEp-2 cells at 24 h. Moreover, 44% strains were cytotonic towards CHO cells and 46% of strains invaded epithelial cells

Knowledge of Patients with Mechanical Valve Prostheses Concerning Chronic Oral Anticoagulant Therapy

Adherence to chronic Oral Anticoagulant Therapy (OAT) is directly related to patients&#39; understanding of the therapy. This study verified the knowledge of patients with mechanical valve prostheses concerning OAT. This is a contemporary cross-sectional study, the sample is composed of patients with mechanical valve prostheses (MVP) in outpatient follow-up. A 10-question instrument was used; answers were ‘know&#39; (1 point), ‘know partially&#39; (half point), or ‘do not know&#39; (zero). Patients were grouped according to the result obtained: &#8804; 4 points was considered insufficient knowledge; > 4 &#8804; 8 equated to moderate knowledge; and > 8 was considered appropriate knowledge. Of the 110 patients, 61.8% presented moderate knowledge, 40.9% were not able to name at least one factor that may alter the INR (International Normalized Ratio) and 37.3% were not able to report their INR target range. The majority of patients presented moderate knowledge concerning the treatment. Strategies to improve knowledge on the topic should be implemented to minimize risks.El objetivo de este estudio fue verificar el conocimiento de pacientes con prótesis valvular sobre la terapia de anticoagulación oral. La adhesión a la terapia de anticoagulación oral crónica está directamente relacionada al entendimiento del paciente sobre esta terapia. Se trata de un estudio transversal contemporáneo, en el cual se incluyeron pacientes de ambulatorio con prótesis valvular mecánica (PVM). Se aplicó un cuestionario con 10 preguntas; las respuestas fueron "sabe" (1 punto), "sabe parcialmente" (medio punto), "no sabe" (cero puntos). Los pacientes fueron agrupados según la puntuación obtenida: &#8804; 4 puntos, conocimiento insuficiente; > 4 &#8804; 8 puntos, conocimiento regular; > 8 puntos, conocimiento adecuado. De los 110 pacientes incluidos, el 61,8% demostró conocimiento regular, 40,9% no supieron citar un factor que altera el REI (razón estandarizada internacional) y 37,3% no supieron informar su REI objetivo. Se concluyó que la mayoría de los pacientes demostró conocimiento regular sobre el tratamiento. Se recomienda implantar estrategias para mejorar el conocimiento y minimizar los riesgos de la terapia.A adesão à terapia com anticoagulação oral crônica está diretamente relacionada ao entendimento dos pacientes sobre essa terapia. O objetivo deste estudo foi verificar o conhecimento dos pacientes, com prótese valvar mecânica, sobre terapia de anticoagulação oral. Como método, usou-se o estudo transversal contemporâneo. Incluem-se pacientes com prótese valvar mecânica (PVM) em acompanhamento ambulatorial. Aplicou-se instrumento com 10 questões. As respostas eram sabe (um ponto), sabe parcialmente (meio ponto), ou não sabe (zero). Os pacientes foram agrupados conforme a pontuação obtida. Considerou-se &#8804;4 pontos conhecimento insuficiente, >4 &#8804;8 conhecimento regular e >8 conhecimento adequado. Os resultados mostram que, dos 110 pacientes, 61,8% apresentaram conhecimento regular, 40,9% não souberam citar ao menos um fator que alterasse a razão normatizada internacional (RNI) e 37,3% não souberam informar o seu RNI alvo. Conclui-se que a maioria dos pacientes apresentou conhecimento regular sobre o tratamento. Estratégias devem ser implementadas para melhorar o conhecimento e, consequentemente, minimizar os riscos dessa terapêutica

Causes and Correlates of Calf Mortality in Captive Asian Elephants (Elephas maximus)

Juvenile mortality is a key factor influencing population growth rate in density-independent, predation-free, well-managed captive populations. Currently at least a quarter of all Asian elephants live in captivity, but both the wild and captive populations are unsustainable with the present fertility and calf mortality rates. Despite the need for detailed data on calf mortality to manage effectively populations and to minimize the need for capture from the wild, very little is known of the causes and correlates of calf mortality in Asian elephants. Here we use the world's largest multigenerational demographic dataset on a semi-captive population of Asian elephants compiled from timber camps in Myanmar to investigate the survival of calves (n = 1020) to age five born to captive-born mothers (n = 391) between 1960 and 1999. Mortality risk varied significantly across different ages and was higher for males at any age. Maternal reproductive history was associated with large differences in both stillbirth and liveborn mortality risk: first-time mothers had a higher risk of calf loss as did mothers producing another calf soon (<3.7 years) after a previous birth, and when giving birth at older age. Stillbirth (4%) and pre-weaning mortality (25.6%) were considerably lower than those reported for zoo elephants and used in published population viability analyses. A large proportion of deaths were caused by accidents and lack of maternal milk/calf weakness which both might be partly preventable by supplementary feeding of mothers and calves and work reduction of high-risk mothers. Our results on Myanmar timber elephants with an extensive keeping system provide an important comparison to compromised survivorship reported in zoo elephants. They have implications for improving captive working elephant management systems in range countries and for refining population viability analyses with realistic parameter values in order to predict future population size of the Asian elephant

Protocol of a Randomized Controlled Trial of Culturally Sensitive Interventions to Improve African Americans' and Non-African Americans' Early, Shared, and Informed Consideration of Live Kidney Transplantation: The talking about Live Kidney Donation (TALK) study

<p>Abstract</p> <p>Background</p> <p>Live kidney transplantation (LKT) is underutilized, particularly among ethnic/racial minorities. The effectiveness of culturally sensitive educational and behavioral interventions to encourage patients' early, shared (with family and health care providers) and informed consideration of LKT and ameliorate disparities in consideration of LKT is unknown.</p> <p>Methods/Design</p> <p>We report the protocol of the Talking About Live Kidney Donation (TALK) Study, a two-phase study utilizing qualitative and quantitative research methods to design and test culturally sensitive interventions to improve patients' shared and informed consideration of LKT. Study Phase 1 involved the evidence-based development of culturally sensitive written and audiovisual educational materials as well as a social worker intervention to encourage patients' engagement in shared and informed consideration of LKT. In Study Phase 2, we are currently conducting a randomized controlled trial in which participants with progressing chronic kidney disease receive: 1) usual care by their nephrologists, 2) usual care plus the educational materials, or 3) usual care plus the educational materials and the social worker intervention. The primary outcome of the randomized controlled trial will include patients' self-reported rates of consideration of LKT (including family discussions of LKT, patient-physician discussions of LKT, and identification of an LKT donor). We will also assess differences in rates of consideration of LKT among African Americans and non-African Americans.</p> <p>Discussion</p> <p>The TALK Study rigorously developed and is currently testing the effectiveness of culturally sensitive interventions to improve patients' and families' consideration of LKT. Results from TALK will provide needed evidence on ways to enhance consideration of this optimal treatment for patients with end stage renal disease.</p> <p>Trial Registration</p> <p>ClinicalTrials.gov number, <a href="http://www.clinicaltrials.gov/ct2/show/NCT00932334">NCT00932334</a></p

How Does Spatial Study Design Influence Density Estimates from Spatial Capture-Recapture Models?

When estimating population density from data collected on non-invasive detector arrays, recently developed spatial capture-recapture (SCR) models present an advance over non-spatial models by accounting for individual movement. While these models should be more robust to changes in trapping designs, they have not been well tested. Here we investigate how the spatial arrangement and size of the trapping array influence parameter estimates for SCR models. We analysed black bear data collected with 123 hair snares with an SCR model accounting for differences in detection and movement between sexes and across the trapping occasions. To see how the size of the trap array and trap dispersion influence parameter estimates, we repeated analysis for data from subsets of traps: 50% chosen at random, 50% in the centre of the array and 20% in the South of the array. Additionally, we simulated and analysed data under a suite of trap designs and home range sizes. In the black bear study, we found that results were similar across trap arrays, except when only 20% of the array was used. Black bear density was approximately 10 individuals per 100 km2. Our simulation study showed that SCR models performed well as long as the extent of the trap array was similar to or larger than the extent of individual movement during the study period, and movement was at least half the distance between traps. SCR models performed well across a range of spatial trap setups and animal movements. Contrary to non-spatial capture-recapture models, they do not require the trapping grid to cover an area several times the average home range of the studied species. This renders SCR models more appropriate for the study of wide-ranging mammals and more flexible to design studies targeting multiple species

Split-based computation of majority-rule supertrees

<p>Abstract</p> <p>Background</p> <p>Supertree methods combine overlapping input trees into a larger supertree. Here, I consider split-based supertree methods that first extract the split information of the input trees and subsequently combine this split information into a phylogeny. Well known split-based supertree methods are matrix representation with parsimony and matrix representation with compatibility. Combining input trees on the same taxon set, as in the consensus setting, is a well-studied task and it is thus desirable to generalize consensus methods to supertree methods.</p> <p>Results</p> <p>Here, three variants of majority-rule (MR) supertrees that generalize majority-rule consensus trees are investigated. I provide simple formulas for computing the respective score for bifurcating input- and supertrees. These score computations, together with a heuristic tree search minmizing the scores, were implemented in the python program PluMiST (Plus- and Minus SuperTrees) available from <url>http://www.cibiv.at/software/plumist</url>. The different MR methods were tested by simulation and on real data sets. The search heuristic was successful in combining compatible input trees. When combining incompatible input trees, especially one variant, MR(-) supertrees, performed well.</p> <p>Conclusions</p> <p>The presented framework allows for an efficient score computation of three majority-rule supertree variants and input trees. I combined the score computation with a heuristic search over the supertree space. The implementation was tested by simulation and on real data sets and showed promising results. Especially the MR(-) variant seems to be a reasonable score for supertree reconstruction. Generalizing these computations to multifurcating trees is an open problem, which may be tackled using this framework.</p

Efficient representation of uncertainty in multiple sequence alignments using directed acyclic graphs

Background A standard procedure in many areas of bioinformatics is to use a single multiple sequence alignment (MSA) as the basis for various types of analysis. However, downstream results may be highly sensitive to the alignment used, and neglecting the uncertainty in the alignment can lead to significant bias in the resulting inference. In recent years, a number of approaches have been developed for probabilistic sampling of alignments, rather than simply generating a single optimum. However, this type of probabilistic information is currently not widely used in the context of downstream inference, since most existing algorithms are set up to make use of a single alignment. Results In this work we present a framework for representing a set of sampled alignments as a directed acyclic graph (DAG) whose nodes are alignment columns; each path through this DAG then represents a valid alignment. Since the probabilities of individual columns can be estimated from empirical frequencies, this approach enables sample-based estimation of posterior alignment probabilities. Moreover, due to conditional independencies between columns, the graph structure encodes a much larger set of alignments than the original set of sampled MSAs, such that the effective sample size is greatly increased. Conclusions The alignment DAG provides a natural way to represent a distribution in the space of MSAs, and allows for existing algorithms to be efficiently scaled up to operate on large sets of alignments. As an example, we show how this can be used to compute marginal probabilities for tree topologies, averaging over a very large number of MSAs. This framework can also be used to generate a statistically meaningful summary alignment; example applications show that this summary alignment is consistently more accurate than the majority of the alignment samples, leading to improvements in downstream tree inference. Implementations of the methods described in this article are available at http://statalign.github.io/WeaveAlign webcite