A prime number approach to biological sequencing

Computational sequencing of nucleic acid and amino acid sequences is placing increasing demands on computer resources. The use of prime numbers is explored as a convenient means of improving program speed and reducing storage requirements. It is concluded that the application of the prime number approach leads to significant increases in speed and some reduction in storage requirements.     

A novel conceptual approach to read-filtering in high-throughput amplicon sequencing studies

Adequate read filtering is critical when processing high-throughput data in marker-gene-based studies. Sequencing errors can cause the mis-clustering of otherwise similar reads, artificially increasing the number of retrieved Operational Taxonomic Units (OTUs) and therefore leading to the overestimation of microbial diversity. Sequencing errors will also result in OTUs that are not accurate reconstructions of the original biological sequences. Herein we present the Poisson binomial filtering algorithm (PBF), which minimizes both problems by calculating the error-probability distribution of a sequence from its quality scores. In order to validate our method, we quality-filtered 37 publicly available datasets obtained by sequencing mock and environmental microbial communities with the Roche 454, Illumina MiSeq and IonTorrent PGM platforms, and compared our results to those obtained with previous approaches such as the ones included in mothur, QIIME and USEARCH. Our algorithm retained substantially more reads than its predecessors, while resulting in fewer and more accurate OTUs. This improved sensitiveness produced more faithful representations, both quantitatively and qualitatively, of the true microbial diversity present in the studied samples. Furthermore, the method introduced in this work is computationally inexpensive and can be readily applied in conjunction with any existent analysis pipeline.     

Triadic conceptual structure of the maximum entropy approach to evolution

Many problems in evolutionary theory are cast in dyadic terms, such as the polar oppositions of organism and environment. We argue that a triadic conceptual structure offers an alternative perspective under which the information generating role of evolution as a physical process can be analyzed, and propose a new diagrammatic approach. Peirce's natural philosophy was deeply influenced by his reception of both Darwin's theory and thermodynamics. Thus, we elaborate on a new synthesis which puts together his theory of signs and modern Maximum Entropy approaches to evolution in a process discourse. Following recent contributions to the naturalization of Peircean semiosis, pointing towards ‘physiosemiosis’ or ‘pansemiosis’, we show that triadic structures involve the conjunction of three different kinds of causality, efficient, formal and final. In this, we accommodate the state-centered thermodynamic framework to a process approach. We apply this on Ulanowicz's analysis of autocatalytic cycles as primordial patterns of life. This paves the way for a semiotic view of thermodynamics which is built on the idea that Peircean interpretants are systems of physical inference devices evolving under natural selection. In this view, the principles of Maximum Entropy, Maximum Power, and Maximum Entropy Production work together to drive the emergence of information carrying structures, which at the same time maximize information capacity as well as the gradients of energy flows, such that ultimately, contrary to Schrödinger's seminal contribution, the evolutionary process is seen to be a physical expression of the Second Law.     

Towards a new conceptual approach to "parasitoproteomics"

Many parasitologists are betting heavily on proteomic studies to explain biochemical host-parasite interactions and, thus, to contribute to disease control. However, many "parasitoproteomic" studies are performed with powerful techniques but without a conceptual approach to determine whether the host genomic responses during a parasite infection represent a nonspecific response that might be induced by any parasite or any other stress. In this article, a new conceptual approach, based on evolutionary concepts of immune responses of a host to a parasite, is suggested for parasitologists to study the host proteome reaction after parasite invasion. Also, this new conceptual approach can be used to study other host-parasite interactions such as behavioral manipulation.     

The conceptual approach to quantitative modeling of guard cells

Much of the 70% of global water usage associated with agriculture passes through stomatal pores of plant leaves. The guard cells, which regulate these pores, thus have a profound influence on photosynthetic carbon assimilation and water use efficiency of plants. We recently demonstrated how quantitative mathematical modeling of guard cells with the OnGuard modeling software yields detail sufficient to guide phenotypic and mutational analysis. This advance represents an all-important step toward applications in directing “reverse-engineering” of guard cell function for improved water use efficiency and carbon assimilation. OnGuard is nonetheless challenging for those unfamiliar with a modeler’s way of thinking. In practice, each model construct represents a hypothesis under test, to be discarded, validated or refined by comparisons between model predictions and experimental results. The few guidelines set out here summarize the standard and logical starting points for users of the OnGuard software.     

A segmental maximum a posteriori approach to genome-wide copy number profiling

MOTIVATION: Copy number profiling methods aim at assigning DNA copy numbers to chromosomal regions using measurements from microarray-based comparative genomic hybridizations. Among the proposed methods to this end, Hidden Markov Model (HMM)-based approaches seem promising since DNA copy number transitions are naturally captured in the model. Current discrete-index HMM-based approaches do not, however, take into account heterogeneous information regarding the genomic overlap between clones. Moreover, the majority of existing methods are restricted to chromosome-wise analysis. RESULTS: We introduce a novel Segmental Maximum A Posteriori approach, SMAP, for DNA copy number profiling. Our method is based on discrete-index Hidden Markov Modeling and incorporates genomic distance and overlap between clones. We exploit a priori information through user-controllable parameterization that enables the identification of copy number deviations of various lengths and amplitudes. The model parameters may be inferred at a genome-wide scale to avoid overfitting of model parameters often resulting from chromosome-wise model inference. We report superior performances of SMAP on synthetic data when compared with two recent methods. When applied on our new experimental data, SMAP readily recognizes already known genetic aberrations including both large-scale regions with aberrant DNA copy number and changes affecting only single features on the array. We highlight the differences between the prediction of SMAP and the compared methods and show that SMAP accurately determines copy number changes and benefits from overlap consideration.     

A stereologic approach to estimate the number of immunogold-labeled molecules in cells of tubules

OBJECTIVE: The quantity of molecules can be measured very precisely by molecular biological methods, but the capabilities of these are limited to measure only the total mass of tissue. For estimating the number of molecules at the cell level, it is necessary to combine an immunohistochemical protocol with designed-based principles of stereology at the level of electron microscopy (EM). This article focuses on the problems and practical solutions of fitting together immunohistochemistry, stereology, and electron microscopy for the estimation of the number of angiotensin II AT1 receptors in rat kidney arterioles. STUDY DESIGN: We performed the preembedding immunostaining of angiotensin II AT1 receptors using the silver-enhanced immunogold labeling system at EM level on serial sections of renal arterioles from 5 rats. RESULTS: Using this method the number of molecules can be estimated along the renal arterioles separately on the cell's surface, in cytoplasm, in nucleus, or in any subcellular location. CONCLUSION: For estimating the number of AT1 receptors, we designed a protocol that took into account the requirements for both immuno-EM and stereology. This method can be applied for estimating any molecule number in different types of cells in tubules.     

A conceptual contribution to battles in the brain

Badcock and Crespi have advanced the hypothesis that autism and schizophrenia are caused by imbalanced imprinting in the brain. They argue that an imbalance between the effects of paternally and maternally expressed genes on brain development results in either an extreme paternal (autism) or maternal brain (schizophrenia). In this paper their conceptual model is discussed and criticized since it presupposes an incoherent distinction between observable physical and hidden mental phenomena. An alternative model is discussed that may be more fruitful for investigating the possible role of imprinted genes in the development of social behaviour. The development of crying and reactive crying and behaviours necessary for collaborative action are discussed as a promising research area for understanding the effects of imprinted genes.     

A distance-type measure approach to the analysis of copy number variation in DNA sequencing data

Background

The next generation sequencing technology allows us to obtain a large amount of short DNA sequence (DNA-seq) reads at a genome-wide level. DNA-seq data have been increasingly collected during the recent years. Count-type data analysis is a widely used approach for DNA-seq data. However, the related data pre-processing is based on the moving window method, in which a window size need to be defined in order to obtain count-type data. Furthermore, useful information can be reduced after data pre-processing for count-type data.

Results

In this study, we propose to analyze DNA-seq data based on the related distance-type measure. Distances are measured in base pairs (bps) between two adjacent alignments of short reads mapped to a reference genome. Our experimental data based simulation study confirms the advantages of distance-type measure approach in both detection power and detection accuracy. Furthermore, we propose artificial censoring for the distance data so that distances larger than a given value are considered potential outliers. Our purpose is to simplify the pre-processing of DNA-seq data. Statistically, we consider a mixture of right censored geometric distributions to model the distance data. Additionally, to reduce the GC-content bias, we extend the mixture model to a mixture of generalized linear models (GLMs). The estimation of model can be achieved by the Newton-Raphson algorithm as well as the Expectation-Maximization (E-M) algorithm. We have conducted simulations to evaluate the performance of our approach. Based on the rank based inverse normal transformation of distance data, we can obtain the related z-values for a follow-up analysis. For an illustration, an application to the DNA-seq data from a pair of normal and tumor cell lines is presented with a change-point analysis of z-values to detect DNA copy number alterations.

Conclusion

Our distance-type measure approach is novel. It does not require either a fixed or a sliding window procedure for generating count-type data. Its advantages have been demonstrated by our simulation studies and its practical usefulness has been illustrated by an experimental data application.

     

A straightforward conceptual approach for evaluating spatial conservation priorities under climate change

Despite wide evidence of a quickly changing world, systematic conservation planning analyses are usually static assuming that the biodiversity being preserved in sites do not change through time. Here we generated a comprehensive ensemble forecasting experiment for 444 amphibian species inhabiting the Atlantic Forest Biodiversity Hotspot. Models were based on four methods for modeling ecological niches, and three future climate simulations. Combinations of these models were used to estimate species occurrences. We used species occurrences to optimize the current and future representation of amphibians with different conservation targets based on their geographic range size. We compared spatial priority outcomes (variance of site selection frequency scores) under dynamic conditions, using a bi-dimensional plot in which the relative importance of each site in achieving conservation targets was assessed both for current time and to 2050. Projections for 2050 show that species richness pattern will remain approximately constant, whereas high turnover rates are forecasted. Selection frequency of several locations varied widely, with recurrent sites located at the north and southeast of the biome. As for 2050, spatial priorities concentrate in the northern part of the biome. Thirty-three sites have high priority for conservation as they play an important role now and will still stand as priority locations in 2050. We present a conceptual model for dynamic spatial conservation prioritization that helps to identify priority sites under climate change. We also call attention to sites in which risk of investment is high, and to those that may become interesting options in the future.     

A new approach to estimation of the number of central synapse(s) included in the H-reflex

Background  

Among the main clinical applications of the H-reflex are the evaluation of the S1 nerve root conductivity such as radiculopathy and measurement of the excitability of the spinal motoneurons in neurological conditions. An attempt has been made to reduce the pathway over which H-reflex can be obtained in a hope to localize a lesion to the S1 nerve root, so the S1 central loop has been suggested. The main goal of this study is the estimation of the H-reflex number of synapse(s) for better understanding of the physiology of this practical reflex.     

A topological approach to nucleosome structure and dynamics: the linking number paradox and other issues.

The linking number paradox of DNA in chromatin (two negative crossings around the octamer, associated with a unit linking number reduction), which is 21 years old this year, has come of age. After stirring much debate in the past, the initially hypothetical explanation of the paradox by DNA overtwisting on the nucleosome surface is now presented as a hard fact in recent textbooks. The first part of this article presents a historical perspective of the problem and details the numerous attempts to measure DNA local periodicity, which in one remarkable example sowed the seeds for the discovery of DNA bending. The second part is devoted to the DNA minicircle system, which has been developed in the author's laboratory as an alternative to the local-periodicity-measurement approach. It offers a simple proposal: a unit linking number reduction associated with a single crossing. This conclusion is contrasted with the latest high-resolution crystallographic data of the nucleosome in the third part of the article, and the fourth part examines the available evidence supporting an extension of these results to nucleosomes in chromatin. The last part addresses another basic question pertaining to nucleosome dynamics, the conformational flexibility of the histone tetramer.     

Determinants of the eradicability of filarial infections: a conceptual approach

Lymphatic filariasis and onchocerciasis are subject to major intervention programs by the WHO. The Onchocerciasis Control Programme in West Africa was launched 30 years ago and has led to considerable insights into the control of this infection. The Global Alliance to Eliminate Lymphatic Filariasis is a relatively recent control program with ambitious targets concerning its efficacy and its schedule. These expectations, however, are based on certain assumptions about the density-dependent processes of limitation and facilitation which determine eradicability: the levels of transmission thresholds and breakpoints. Here, we review these processes operating in filarial infections and show their impact on the persistence of the parasite, as well as pointing out those issues where more information is required to develop sound predictions about the eradicability of these infections.     

Positron emission tomography: the conceptual idea using a multidisciplinary approach

Positron emission tomography (PET) is a method for quantitatively measuring biochemical and physiological processes in vivo by using radiopharmaceuticals labeled with positron-emitting radionuclides such as 11C, 13N, 15O, and 18F and by measuring the annihilation radiation using a coincidence technique. This technique is also used for measurement of the pharmacokinetics of labeled drugs and measurement of the effects of drugs on metabolism. Deviations from normal metabolism can be measured and insight into biological processes responsible for diseases can be obtained.