Parametric construction of episode networks from pseudoperiodic time series based on mutual information

Recently, the construction of networks from time series data has gained widespread interest. In this paper, we develop this area further by introducing a network construction procedure for pseudoperiodic time series. We call such networks episode networks, in which an episode corresponds to a temporal interval of a time series, and which defines a node in the network. Our model includes a number of features which distinguish it from current methods. First, the proposed construction procedure is a parametric model which allows it to adapt to the characteristics of the data; the length of an episode being the parameter. As a direct consequence, networks of minimal size containing the maximal information about the time series can be obtained. In this paper, we provide an algorithm to determine the optimal value of this parameter. Second, we employ estimates of mutual information values to define the connectivity structure among the nodes in the network to exploit efficiently the nonlinearities in the time series. Finally, we apply our method to data from electroencephalogram (EEG) experiments and demonstrate that the constructed episode networks capture discriminative information from the underlying time series that may be useful for diagnostic purposes.     

Fusing time series expression data through hybrid aggregation and hierarchical merge

SUMMARY: A novel integration approach targeting the combination of multi-experiment time series expression data is proposed. A recursive hybrid aggregation algorithm is initially employed to extract a set of genes, which are eventually of interest for the biological phenomenon under study. Next, a hierarchical merge procedure is specifically developed for the purpose of fusing together the multiple-experiment expression pro.les of the selected genes. This employs dynamic time warping alignment techniques in order to account adequately for the potential phase shift between the different experiments. We subsequently demonstrate that the resulting gene expression pro.les consistently re.ect the behavior of the original expression pro.les in the different experiments. SUPPLEMENTARY INFORMATION: Supplementary data are available athttp://www.tu-plovdiv.bg/Container/bi/DataIntegration/     

Interocular yoking in human saccades examined by mutual information analysis

Background

Saccadic eye movements align the two eyes precisely to foveate a target. Trial-by-trial variance of eye movement is always observed within an identical experimental condition. This has often been treated as experimental error without addressing its significance. The present study examined statistical linkages between the two eyes’ movements, namely interocular yoking, for the variance of eye position and velocity.

Methods

Horizontal saccadic movements were recorded from twelve right-eye-dominant subjects while they decided on saccade direction in Go-Only sessions and on both saccade execution and direction in Go/NoGo sessions. We used infrared corneal reflection to record simultaneously and independently the movement of each eye. Quantitative measures of yoking were provided by mutual information analysis of eye position or velocity, which is sensitive to both linear and non-linear relationships between the eyes’ movements. Our mutual information analysis relied on the variance of the eyes movements in each experimental condition. The range of movements for each eye varies for different conditions so yoking was further studied by comparing GO-Only vs. Go/NoGo sessions, leftward vs. rightward saccades.

Results

Mutual information analysis showed that velocity yoking preceded positional yoking. Cognitive load increased trial variances of velocity with no increase in velocity yoking, suggesting that cognitive load may alter neural processes in areas to which oculomotor control is not tightly linked. The comparison between experimental conditions showed that interocular linkage in velocity variance of the right eye lagged that of the left eye during saccades.

Conclusions

We conclude quantitative measure of interocular yoking based on trial-to-trial variance within a condition, as well as variance between conditions, provides a powerful tool for studying the binocular movement mechanism.

     

Influence of statistical estimators of mutual information and data heterogeneity on the inference of gene regulatory networks

The inference of gene regulatory networks from gene expression data is a difficult problem because the performance of the inference algorithms depends on a multitude of different factors. In this paper we study two of these. First, we investigate the influence of discrete mutual information (MI) estimators on the global and local network inference performance of the C3NET algorithm. More precisely, we study 4 different MI estimators (Empirical, Miller-Madow, Shrink and Schürmann-Grassberger) in combination with 3 discretization methods (equal frequency, equal width and global equal width discretization). We observe the best global and local inference performance of C3NET for the Miller-Madow estimator with an equal width discretization. Second, our numerical analysis can be considered as a systems approach because we simulate gene expression data from an underlying gene regulatory network, instead of making a distributional assumption to sample thereof. We demonstrate that despite the popularity of the latter approach, which is the traditional way of studying MI estimators, this is in fact not supported by simulated and biological expression data because of their heterogeneity. Hence, our study provides guidance for an efficient design of a simulation study in the context of network inference, supporting a systems approach.     

Integrated analysis of microarray data and gene function information

Microarray data should be interpreted in the context of existing biological knowledge. Here we present integrated analysis of microarray data and gene function classification data using homogeneity analysis. Homogeneity analysis is a graphical multivariate statistical method for analyzing categorical data. It converts categorical data into graphical display. By simultaneously quantifying the microarray-derived gene groups and gene function categories, it captures the complex relations between biological information derived from microarray data and the existing knowledge about the gene function. Thus, homogeneity analysis provides a mathematical framework for integrating the analysis of microarray data and the existing biological knowledge.     

Time series analysis of aerobic bacterial flora during Miso fermentation

AIMS: This article reports a microbiological study of aerobic mesophilic bacteria that are present during the fermentation process of Miso. METHODS AND RESULTS: Aerobic bacteria were enumerated and isolated from Miso during fermentation and divided into nine groups using traditional phenotypic tests. The strains were identified by biochemical analysis and 16S rRNA sequence analysis. They were identified as Bacillus subtilis, B. amyloliquefaciens, Kocuria kristinae, Staphylococcus gallinarum and S. kloosii. All strains were sensitive to the bacteriocins produced by the lactic acid bacteria isolated from Miso. CONCLUSIONS: The dominant species among the undesirable species throughout the fermentation process were B. subtilis and B. amyloliquefaciens. It is suggested that bacteriocin-producing lactic acid bacteria are effective in the growth prevention of aerobic bacteria in Miso. SIGNIFICANCE AND IMPACT OF THE STUDY: This study has provided useful information for controlling of bacterial flora during Miso fermentation.     

Genome information management and integrated data analysis with HaloLex

HaloLex is a software system for the central management, integration, curation, and web-based visualization of genomic and other -omics data for any given microorganism. The system has been employed for the manual curation of three haloarchaeal genomes, namely Halobacterium salinarum (strain R1), Natronomonas pharaonis, and Haloquadratum walsbyi. HaloLex, in particular, enables the integrated analysis of genome-wide proteomic results with the underlying genomic data. This has proven indispensable to generate reliable gene predictions for GC-rich genomes, which, due to their characteristically low abundance of stop codons, are known to be hard targets for standard gene finders, especially concerning start codon assignment. The proteomic identification of more than 600 N-terminal peptides has greatly increased the reliability of the start codon assignment for Halobacterium salinarum. Application of homology-based methods to the published genome of Haloarcula marismortui allowed to detect 47 previously unidentified genes (a problem that is particularly serious for short protein sequences) and to correct more than 300 start codon misassignments.     

Judgment algorithm for periodicity of time series data based on Bayesian information criterion

Judgment periodicity of biological time series data is important and done widely to find the circadian expression of genes, monthly change of hormones, etc. To keep complete reproducibility of judgment is a problem because popular judgment methods such as curve fitting, Fourier analysis, etc. need judgment criteria determined by analysts considering experimental conditions and results (level, S/N, distribution, etc.) based on their experience. Judgment results are often affected by analysts' subjects. Reproducible criterion determination is therefore strongly needed. We propose introducing the information criterion to replace analysts' criteria. A judgment algorithm by combining Bayesian information criterion (BIC) and discrete Fourier transform (DFT) has been developed and has proved its ability through application to mice microarray data and finding of circadian genes. Our method, named "Piccolo", shows higher sensitivity than the simple DFT (without BIC) method with reproducibility, and can be fully automated.     

Time series in analysis of yerba-mate biennial growth modified by environment

To assess differences in the lag-effect pattern in the relationship between yerba-mate biennial growth and environmental factors, a time-series analysis was performed. A generalized Poisson regression model was used to control time trends, temperature, growing degree days (GDD), rainfalls and night length (NL). It was hypothesized that the active growth and growth pauses in yerba-mate are controlled endogenously and modified by environment, and that genders would respond differently to environmental modifications. The patterns in the lag effect from the distributed-lag models were similar to those of time-series models with meteorological data means with lag = 0. GDD and NL were principal factors affecting biennial yerba-mate shoot elongation and the number of green leaves of females grown in monoculture, besides their significant effects on metamer emission and leaf area in males grown in monoculture. NL also had a significant influence on shoot elongation and leaf area of both genders grown in forest understorey (FUS), indicating that yerba-mate growth is synchronized by an internal clock sensitive to temperature adjustments. The morphological plasticity and the adaptation efforts of yerba-mate were more pronounced in monoculture than in FUS. Sexual dimorphism was expressed—males were more sensitive to environmental changes than females, especially in monoculture. Growth modifications were much more intense when plants were grown in a cultivation system that is less like yerba-mate natural habitat (monoculture) than in one resembling its natural habitat (FUS). Our data support the ecological specialization theory.     

Computing and visually analyzing mutual information in molecular co-evolution

Background  

Selective pressure in molecular evolution leads to uneven distributions of amino acids and nucleotides. In fact one observes correlations among such constituents due to a large number of biophysical mechanisms (folding properties, electrostatics, ...). To quantify these correlations the mutual information -after proper normalization - has proven most effective. The challenge is to navigate the large amount of data, which in a study for a typical protein cannot simply be plotted.