UNLABELLED: : Commonly used multiplicity adjustments fail to control the error rate for reported findings in many expression quantitative trait loci (eQTL) studies. TreeQTL implements a hierarchical multiple testing procedure which allows control of appropriate error rates defined relative to a grouping of the eQTL hypotheses. AVAILABILITY AND IMPLEMENTATION: The R package TreeQTL is available for download at http://bioinformatics.org/treeqtl CONTACT: [email protected] SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Spatial memory is often studied in the Morris Water Maze, where the animal's spatial orientation has been shown to be mainly shaped by distal visual cues. Cognition-related behavior has also been described along "well-trodden paths"--spatial habits established by animals in the wild and in captivity reflecting a form of spatial memory. In the present study we combine the study of Open Field behavior with the study of behavior on well-trodden paths, revealing a form of locational memory that appears to correlate with spatial memory. The tracked path of the mouse is used to examine the dynamics of visiting behavior to locations. A visit is defined as either progressing through a location or stopping there, where progressing and stopping are computationally defined. We then estimate the probability of stopping at a location as a function of the number of previous visits to that location, i.e., we measure the effect of visiting history to a location on stopping in it. This can be regarded as an estimate of the familiarity of the mouse with locations. The recently wild-derived inbred strain CZECHII shows the highest effect of visiting history on stopping, C57 inbred mice show a lower effect, and DBA mice show no effect. We employ a rarely used, bottom-to-top computational approach, starting from simple kinematics of movement and gradually building our way up until we end with (emergent) locational memory. The effect of visiting history to a location on stopping in it can be regarded as an estimate of the familiarity of the mouse with locations, implying memory of these locations. We show that the magnitude of this estimate is strain-specific, implying a genetic influence. The dynamics of this process reveal that locations along the mouse's trodden path gradually become places of attraction, where the mouse stops habitually.

UNLABELLED: Gene expression and phenotypic functionality can best be associated when they are measured quantitatively within the same experiment. The analysis of such a complex experiment is presented, searching for associations between measures of exploratory behavior in mice and gene expression in brain regions. The analysis of such experiments raises several methodological problems. First and foremost, the size of the pool of potential discoveries being screened is enormous yet only few biologically relevant findings are expected, making the problem of multiple testing especially severe. We present solutions based on screening by testing related hypotheses, then testing the hypotheses of interest. In one variant the subset is selected directly, in the other one a tree of hypotheses is tested hierarchical; both variants control the False Discovery Rate (FDR). Other problems in such experiments are in the fact that the level of data aggregation may be different for the quantitative traits (one per animal) and gene expression measurements (pooled across animals); in that the association may not be linear; and in the resolution of interest only few replications exist. We offer solutions to these problems as well. The hierarchical FDR testing strategies presented here can serve beyond the structure of our motivating example study to any complex microarray study. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

MOTIVATION: DNA microarrays have recently been used for the purpose of monitoring expression levels of thousands of genes simultaneously and identifying those genes that are differentially expressed. The probability that a false identification (type I error) is committed can increase sharply when the number of tested genes gets large. Correlation between the test statistics attributed to gene co-regulation and dependency in the measurement errors of the gene expression levels further complicates the problem. In this paper we address this very large multiplicity problem by adopting the false discovery rate (FDR) controlling approach. In order to address the dependency problem, we present three resampling-based FDR controlling procedures, that account for the test statistics distribution, and compare their performance to that of the naïve application of the linear step-up procedure in Benjamini and Hochberg (1995). The procedures are studied using simulated microarray data, and their performance is examined relative to their ease of implementation. RESULTS: Comparative simulation analysis shows that all four FDR controlling procedures control the FDR at the desired level, and retain substantially more power then the family-wise error rate controlling procedures. In terms of power, using resampling of the marginal distribution of each test statistics substantially improves the performance over the naïve one. The highest power is achieved, at the expense of a more sophisticated algorithm, by the resampling-based procedures that resample the joint distribution of the test statistics and estimate the level of FDR control. AVAILABILITY: An R program that adjusts p-values using FDR controlling procedures is freely available over the Internet at www.math.tau.ac.il/~ybenja.