Polynomial algorithms for the Maximal Pairing Problem: efficient phylogenetic targeting on arbitrary trees

Background: The Maximal Pairing Problem (MPP) is the prototype of a class of combinatorial optimization problems that are of considerable interest in bioinformatics: Given an arbitrary phylogenetic tree T and weights ωxy for the paths between any two pairs of leaves (x, y), what is the collection of edge-disjoint paths between pairs of leaves that maximizes the total weight? Special cases of the MPP for binary trees and equal weights have been described previously; algorithms to solve the general MPP are still missing, however. Results: We describe a relatively simple dynamic programming algorithm for the special case of binary trees. We then show that the general case of multifurcating trees can be treated by interleaving solutions to certain auxiliary Maximum Weighted Matching problems with an extension of this dynamic programming approach, resulting in an overall polynomial-time solution of complexity (n^4 log n) w.r.t. the number n of leaves. The source code of a C implementation can be obtained under the GNU Public License from http://www.bioinf.uni-leipzig.de/Software/Targeting. For binary trees, we furthermore discuss several constrained variants of the MPP as well as a partition function approach to the probabilistic version of the MPP. Conclusions: The algorithms introduced here make it possible to solve the MPP also for large trees with high-degree vertices. This has practical relevance in the field of comparative phylogenetics and, for example, in the context of phylogenetic targeting, i.e., data collection with resource limitations.Human Evolutionary Biolog

Single amino acid change in gp41 region of HIV-1 alters bystander apoptosis and CD4 decline in humanized mice

<p>Abstract</p> <p>Background</p> <p>The mechanism by which HIV infection leads to a selective depletion of CD4 cells leading to immunodeficiency remains highly debated. Whether the loss of CD4 cells is a direct consequence of virus infection or bystander apoptosis of uninfected cells is also uncertain.</p> <p>Results</p> <p>We have addressed this issue in the humanized mouse model of HIV infection using a HIV variant with a point mutation in the gp41 region of the Env glycoprotein that alters its fusogenic activity. We demonstrate here that a single amino acid change (V38E) altering the cell-to-cell fusion activity of the Env minimizes CD4 loss in humanized mice without altering viral replication. This differential pathogenesis was associated with a lack of bystander apoptosis induction by V38E virus even in the presence of similar levels of infected cells. Interestingly, immune activation was observed with both WT and V38E infection suggesting that the two phenomena are likely not interdependent in the mouse model.</p> <p>Conclusions</p> <p>We conclude that Env fusion activity is one of the determinants of HIV pathogenesis and it may be possible to attenuate HIV by targeting gp41.</p

Modeling and Analysis of the Molecular Basis of Pain in Sensory Neurons

Intracellular calcium dynamics are critical to cellular functions like pain transmission. Extracellular ATP plays an important role in modulating intracellular calcium levels by interacting with the P2 family of surface receptors. In this study, we developed a mechanistic mathematical model of ATP-induced P2 mediated calcium signaling in archetype sensory neurons. The model architecture, which described 90 species connected by 162 interactions, was formulated by aggregating disparate molecular modules from literature. Unlike previous models, only mass action kinetics were used to describe the rate of molecular interactions. Thus, the majority of the 252 unknown model parameters were either association, dissociation or catalytic rate constants. Model parameters were estimated from nine independent data sets taken from multiple laboratories. The training data consisted of both dynamic and steady-state measurements. However, because of the complexity of the calcium network, we were unable to estimate unique model parameters. Instead, we estimated a family or ensemble of probable parameter sets using a multi-objective thermal ensemble method. Each member of the ensemble met an error criterion and was located along or near the optimal trade-off surface between the individual training data sets. The model quantitatively reproduced experimental measurements from dorsal root ganglion neurons as a function of extracellular ATP forcing. Hypothesized architecture linking phosphoinositide regulation with P2X receptor activity explained the inhibition of P2X-mediated current flow by activated metabotropic P2Y receptors. Sensitivity analysis using individual and the whole system outputs suggested which molecular subsystems were most important following P2 activation. Taken together, modeling and analysis of ATP-induced P2 mediated calcium signaling generated qualitative insight into the critical interactions controlling ATP induced calcium dynamics. Understanding these critical interactions may prove useful for the design of the next generation of molecular pain management strategies

Patterns of co-speciation and host switching in primate malaria parasites

<p>Abstract</p> <p>Background</p> <p>The evolutionary history of many parasites is dependent on the evolution of their hosts, leading to an association between host and parasite phylogenies. However, frequent host switches across broad phylogenetic distances may weaken this close evolutionary link, especially when vectors are involved in parasites transmission, as is the case for malaria pathogens. Several studies suggested that the evolution of the primate-infective malaria lineages may be constrained by the phylogenetic relationships of their hosts, and that lateral switches between distantly related hosts may have been occurred. However, no systematic analysis has been quantified the degree of phylogenetic association between primates and their malaria parasites.</p> <p>Methods</p> <p>Here phylogenetic approaches have been used to discriminate statistically between events due to co-divergence, duplication, extinction and host switches that can potentially cause historical association between <it>Plasmodium </it>parasites and their primate hosts. A Bayesian reconstruction of parasite phylogeny based on genetic information for six genes served as basis for the analyses, which could account for uncertainties about the evolutionary hypotheses of malaria parasites.</p> <p>Results</p> <p>Related lineages of primate-infective <it>Plasmodium </it>tend to infect hosts within the same taxonomic family. Different analyses testing for congruence between host and parasite phylogenies unanimously revealed a significant association between the corresponding evolutionary trees. The most important factor that resulted in this association was host switching, but depending on the parasite phylogeny considered, co-speciation and duplication may have also played some additional role. Sorting seemed to be a relatively infrequent event, and can occur only under extreme co-evolutionary scenarios. The concordance between host and parasite phylogenies is heterogeneous: while the evolution of some malaria pathogens is strongly dependent on the phylogenetic history of their primate hosts, the congruent evolution is less emphasized for other parasite lineages (e.g. for human malaria parasites). Estimation of ancestral states of host use along the phylogenetic tree of parasites revealed that lateral transfers across distantly related hosts were likely to occur in several cases. Parasites cannot infect all available hosts, and they should preferentially infect hosts that provide a similar environment for reproduction. Marginally significant evidence suggested that there might be a consistent variation within host ranges in terms of physiology.</p> <p>Conclusion</p> <p>The evolution of primate malarias is constrained by the phylogenetic associations of their hosts. Some parasites can preserve a great flexibility to infect hosts across a large phylogenetic distance, thus host switching can be an important factor in mediating host ranges observed in nature. Due to this inherent flexibility and the potential exposure to various vectors, the emergence of new malaria disease in primates including humans cannot be predicted from the phylogeny of parasites.</p

Network analysis of the transcriptional pattern of young and old cells of Escherichia coli during lag phase

Background: The aging process of bacteria in stationary phase is halted if cells are subcultured and enter lag phase and it is then followed by cellular division. Network science has been applied to analyse the transcriptional response, during lag phase, of bacterial cells starved previously in stationary phase for 1 day (young cells) and 16 days (old cells). Results: A genome scale network was constructed for E. coli K-12 by connecting genes with operons, transcription and sigma factors, metabolic pathways and cell functional categories. Most of the transcriptional changes were detected immediately upon entering lag phase and were maintained throughout this period. The lag period was longer for older cells and the analysis of the transcriptome revealed different intracellular activity in young and old cells. The number of genes differentially expressed was smaller in old cells (186) than in young cells (467). Relatively, few genes (62) were up- or down-regulated in both cultures. Transcription of genes related to osmotolerance, acid resistance, oxidative stress and adaptation to other stresses was down-regulated in both young and old cells. Regarding carbohydrate metabolism, genes related to the citrate cycle were up-regulated in young cells while old cells up-regulated the Entner Doudoroff and gluconate pathways and down-regulated the pentose phosphate pathway. In both old and young cells, anaerobic respiration and fermentation pathways were down-regulated, but only young cells up-regulated aerobic respiration while there was no evidence of aerobic respiration in old cells.Numerous genes related to DNA maintenance and replication, translation, ribosomal biosynthesis and RNA processing as well as biosynthesis of the cell envelope and flagellum and several components of the chemotaxis signal transduction complex were up-regulated only in young cells. The genes for several transport proteins for iron compounds were up-regulated in both young and old cells. Numerous genes encoding transporters for carbohydrates and organic alcohols and acids were down-regulated in old cells only. Conclusion: Network analysis revealed very different transcriptional activities during the lag period in old and young cells. Rejuvenation seems to take place during exponential growth by replicative dilution of old cellular components

Population ecology of the sea lamprey (Petromyzon marinus) as an invasive species in the Laurentian Great Lakes and an imperiled species in Europe

The sea lamprey Petromyzon marinus (Linnaeus) is both an invasive non-native species in the Laurentian Great Lakes of North America and an imperiled species in much of its native range in North America and Europe. To compare and contrast how understanding of population ecology is useful for control programs in the Great Lakes and restoration programs in Europe, we review current understanding of the population ecology of the sea lamprey in its native and introduced range. Some attributes of sea lamprey population ecology are particularly useful for both control programs in the Great Lakes and restoration programs in the native range. First, traps within fish ladders are beneficial for removing sea lampreys in Great Lakes streams and passing sea lampreys in the native range. Second, attractants and repellants are suitable for luring sea lampreys into traps for control in the Great Lakes and guiding sea lamprey passage for conservation in the native range. Third, assessment methods used for targeting sea lamprey control in the Great Lakes are useful for targeting habitat protection in the native range. Last, assessment methods used to quantify numbers of all life stages of sea lampreys would be appropriate for measuring success of control in the Great Lakes and success of conservation in the native range