categoryCompare, an analytical tool based on feature annotations

Assessment of high-throughput—omics data initially focuses on relative or raw levels of a particular feature, such as an expression value for a transcript, protein, or metabolite. At a second level, analyses of annotations including known or predicted functions and associations of each individual feature, attempt to distill biological context. Most currently available comparative- and meta-analyses methods are dependent on the availability of identical features across data sets, and concentrate on determining features that are differentially expressed across experiments, some of which may be considered “biomarkers.” The heterogeneity of measurement platforms and inherent variability of biological systems confounds the search for robust biomarkers indicative of a particular condition. In many instances, however, multiple data sets show involvement of common biological processes or signaling pathways, even though individual features are not commonly measured or differentially expressed between them. We developed a methodology, categoryCompare, for cross-platform and cross-sample comparison of high-throughput data at the annotation level. We assessed the utility of the approach using hypothetical data, as well as determining similarities and differences in the set of processes in two instances: (1) denervated skin vs. denervated muscle, and (2) colon from Crohn's disease vs. colon from ulcerative colitis (UC). The hypothetical data showed that in many cases comparing annotations gave superior results to comparing only at the gene level. Improved analytical results depended as well on the number of genes included in the annotation term, the amount of noise in relation to the number of genes expressing in unenriched annotation categories, and the specific method in which samples are combined. In the skin vs. muscle denervation comparison, the tissues demonstrated markedly different responses. The Crohn's vs. UC comparison showed gross similarities in inflammatory response in the two diseases, with particular processes specific to each disease

A Model for the Evolution of Nucleotide Polymerase Directionality

Background: In all known living organisms, every enzyme that synthesizes nucleic acid polymers does so by adding nucleotide 59-triphosphates to the 39-hydroxyl group of the growing chain. This results in the well known 5’?3’ directionality of all DNA and RNA Polymerases. The lack of any alternative mechanism, e.g. addition in a 3’?5 ’ direction, may indicate a very early founder effect in the evolution of life, or it may be the result of a selective pressure against such an alternative. Methodology/Principal Findings: In an attempt to determine whether the lack of an alternative polymerase directionality is the result of a founder effect or evolutionary selection, we have constructed a basic model of early polymerase evolution. This model is informed by the essential chemical properties of the nucleotide polymerization reaction. With this model, we are able to simulate the growth of organisms with polymerases that synthesize either 5’?3 ’ or 3’?5 ’ in isolation or in competition with each other. Conclusions/Significance: We have found that a competition between organisms with 5’?3 ’ polymerases and 3’?5’ polymerases only results in a evolutionarily stable strategy under certain conditions. Furthermore, we have found that mutations lead to a much clearer delineation between conditions that lead to a stable coexistence of these populations and conditions which ultimately lead to success for the 5’?3 ’ form. In addition to presenting a plausible explanation for th

Free energy estimation of short DNA duplex hybridizations

<p>Abstract</p> <p>Background</p> <p>Estimation of DNA duplex hybridization free energy is widely used for predicting cross-hybridizations in DNA computing and microarray experiments. A number of software programs based on different methods and parametrizations are available for the theoretical estimation of duplex free energies. However, significant differences in free energy values are sometimes observed among estimations obtained with various methods, thus being difficult to decide what value is the accurate one.</p> <p>Results</p> <p>We present in this study a quantitative comparison of the similarities and differences among four published DNA/DNA duplex free energy calculation methods and an extended Nearest-Neighbour Model for perfect matches based on triplet interactions. The comparison was performed on a benchmark data set with 695 pairs of short oligos that we collected and manually curated from 29 publications. Sequence lengths range from 4 to 30 nucleotides and span a large GC-content percentage range. For perfect matches, we propose an extension of the Nearest-Neighbour Model that matches or exceeds the performance of the existing ones, both in terms of correlations and root mean squared errors. The proposed model was trained on experimental data with temperature, sodium and sequence concentration characteristics that span a wide range of values, thus conferring the model a higher power of generalization when used for free energy estimations of DNA duplexes under non-standard experimental conditions.</p> <p>Conclusions</p> <p>Based on our preliminary results, we conclude that no statistically significant differences exist among free energy approximations obtained with 4 publicly available and widely used programs, when benchmarked against a collection of 695 pairs of short oligos collected and curated by the authors of this work based on 29 publications. The extended Nearest-Neighbour Model based on triplet interactions presented in this work is capable of performing accurate estimations of free energies for perfect match duplexes under both standard and non-standard experimental conditions and may serve as a baseline for further developments in this area of research.</p

Electrochemical activation and inhibition of neuromuscular systems through modulation of ion concentrations with ion-selective membranes

Conventional functional electrical stimulation aims to restore functional motor activity of patients with disabilities resulting from spinal cord injury or neurological disorders. However, intervention with functional electrical stimulation in neurological diseases lacks an effective implantable method that suppresses unwanted nerve signals. We have developed an electrochemical method to activate and inhibit a nerve by electrically modulating ion concentrations in situ along the nerve. Using ion-selective membranes to achieve different excitability states of the nerve, we observe either a reduction of the electrical threshold for stimulation by up to approximately 40%, or voluntary, reversible inhibition of nerve signal propagation. This low-threshold electrochemical stimulation method is applicable in current implantable neuroprosthetic devices, whereas the on-demand nerve-blocking mechanism could offer effective clinical intervention in disease states caused by uncontrolled nerve activation, such as epilepsy and chronic pain syndromes.Massachusetts Institute of Technology. Faculty Discretionary Research FundNational Institutes of Health (U.S.) (Award UL1 RR 025758)Harvard Catalyst (Grant

P2X receptor-mediated purinergic sensory pathways to the spinal cord dorsal horn

P2X receptors are expressed on different functional groups of primary afferent fibers. P2X receptor-mediated sensory inputs can be either innocuous or nociceptive, depending on which dorsal horn regions receive these inputs. We provide a brief review of P2X receptor-mediated purinergic sensory pathways to different regions in the dorsal horn. These P2X purinergic pathways are identified in normal animals, which provides insights into their physiological functions. Future studies on P2X purinergic pathways in animal models of pathological conditions may provide insights on how P2X receptors play a role in pathological pain states

Thermodynamic Basis for the Emergence of Genomes during Prebiotic Evolution

The RNA world hypothesis views modern organisms as descendants of RNA molecules. The earliest RNA molecules must have been random sequences, from which the first genomes that coded for polymerase ribozymes emerged. The quasispecies theory by Eigen predicts the existence of an error threshold limiting genomic stability during such transitions, but does not address the spontaneity of changes. Following a recent theoretical approach, we applied the quasispecies theory combined with kinetic/thermodynamic descriptions of RNA replication to analyze the collective behavior of RNA replicators based on known experimental kinetics data. We find that, with increasing fidelity (relative rate of base-extension for Watson-Crick versus mismatched base pairs), replications without enzymes, with ribozymes, and with protein-based polymerases are above, near, and below a critical point, respectively. The prebiotic evolution therefore must have crossed this critical region. Over large regions of the phase diagram, fitness increases with increasing fidelity, biasing random drifts in sequence space toward ‘crystallization.’ This region encloses the experimental nonenzymatic fidelity value, favoring evolutions toward polymerase sequences with ever higher fidelity, despite error rates above the error catastrophe threshold. Our work shows that experimentally characterized kinetics and thermodynamics of RNA replication allow us to determine the physicochemical conditions required for the spontaneous crystallization of biological information. Our findings also suggest that among many potential oligomers capable of templated replication, RNAs may have evolved to form prebiotic genomes due to the value of their nonenzymatic fidelity

Spinal afferent neurons projecting to the rat lung and pleura express acid sensitive channels

BACKGROUND: The acid sensitive ion channels TRPV1 (transient receptor potential vanilloid receptor-1) and ASIC3 (acid sensing ion channel-3) respond to tissue acidification in the range that occurs during painful conditions such as inflammation and ischemia. Here, we investigated to which extent they are expressed by rat dorsal root ganglion neurons projecting to lung and pleura, respectively. METHODS: The tracer DiI was either injected into the left lung or applied to the costal pleura. Retrogradely labelled dorsal root ganglion neurons were subjected to triple-labelling immunohistochemistry using antisera against TRPV1, ASIC3 and neurofilament 68 (marker for myelinated neurons), and their soma diameter was measured. RESULTS: Whereas 22% of pulmonary spinal afferents contained neither channel-immunoreactivity, at least one is expressed by 97% of pleural afferents. TRPV1(+)/ASIC3(- )neurons with probably slow conduction velocity (small soma, neurofilament 68-negative) were significantly more frequent among pleural (35%) than pulmonary afferents (20%). TRPV1(+)/ASIC3(+ )neurons amounted to 14 and 10% respectively. TRPV1(-)/ASIC3(+ )neurons made up between 44% (lung) and 48% (pleura) of neurons, and half of them presumably conducted in the A-fibre range (larger soma, neurofilament 68-positive). CONCLUSION: Rat pleural and pulmonary spinal afferents express at least two different acid-sensitive channels that make them suitable to monitor tissue acidification. Patterns of co-expression and structural markers define neuronal subgroups that can be inferred to subserve different functions and may initiate specific reflex responses. The higher prevalence of TRPV1(+)/ASIC3(- )neurons among pleural afferents probably reflects the high sensitivity of the parietal pleura to painful stimuli