Algebraic Comparison of Partial Lists in Bioinformatics

The outcome of a functional genomics pipeline is usually a partial list of genomic features, ranked by their relevance in modelling biological phenotype in terms of a classification or regression model. Due to resampling protocols or just within a meta-analysis comparison, instead of one list it is often the case that sets of alternative feature lists (possibly of different lengths) are obtained. Here we introduce a method, based on the algebraic theory of symmetric groups, for studying the variability between lists ("list stability") in the case of lists of unequal length. We provide algorithms evaluating stability for lists embedded in the full feature set or just limited to the features occurring in the partial lists. The method is demonstrated first on synthetic data in a gene filtering task and then for finding gene profiles on a recent prostate cancer dataset

A simple approach to ranking differentially expressed gene expression time courses through Gaussian process regression.

BACKGROUND: The analysis of gene expression from time series underpins many biological studies. Two basic forms of analysis recur for data of this type: removing inactive (quiet) genes from the study and determining which genes are differentially expressed. Often these analysis stages are applied disregarding the fact that the data is drawn from a time series. In this paper we propose a simple model for accounting for the underlying temporal nature of the data based on a Gaussian process. RESULTS: We review Gaussian process (GP) regression for estimating the continuous trajectories underlying in gene expression time-series. We present a simple approach which can be used to filter quiet genes, or for the case of time series in the form of expression ratios, quantify differential expression. We assess via ROC curves the rankings produced by our regression framework and compare them to a recently proposed hierarchical Bayesian model for the analysis of gene expression time-series (BATS). We compare on both simulated and experimental data showing that the proposed approach considerably outperforms the current state of the art. CONCLUSIONS: Gaussian processes offer an attractive trade-off between efficiency and usability for the analysis of microarray time series. The Gaussian process framework offers a natural way of handling biological replicates and missing values and provides confidence intervals along the estimated curves of gene expression. Therefore, we believe Gaussian processes should be a standard tool in the analysis of gene expression time series

Learning Temporal Patterns of Risk in a Predator-Diverse Environment

Predation plays a major role in shaping prey behaviour. Temporal patterns of predation risk have been shown to drive daily activity and foraging patterns in prey. Yet the ability to respond to temporal patterns of predation risk in environments inhabited by highly diverse predator communities, such as rainforests and coral reefs, has received surprisingly little attention. In this study, we investigated whether juvenile marine fish, Pomacentrus moluccensis (lemon damselfish), have the ability to learn to adjust the intensity of their antipredator response to match the daily temporal patterns of predation risk they experience. Groups of lemon damselfish were exposed to one of two predictable temporal risk patterns for six days. “Morning risk” treatment prey were exposed to the odour of Cephalopholis cyanostigma (rockcod) paired with conspecific chemical alarm cues (simulating a rockcod present and feeding) during the morning, and rockcod odour only in the evening (simulating a rockcod present but not feeding). “Evening risk” treatment prey had the two stimuli presented to them in the opposite order. When tested individually for their response to rockcod odour alone, lemon damselfish from the morning risk treatment responded with a greater antipredator response intensity in the morning than in the evening. In contrast, those lemon damselfish previously exposed to the evening risk treatment subsequently responded with a greater antipredator response when tested in the evening. The results of this experiment demonstrate that P. moluccensis have the ability to learn temporal patterns of predation risk and can adjust their foraging patterns to match the threat posed by predators at a given time of day. Our results provide the first experimental demonstration of a mechanism by which prey in a complex, multi-predator environment can learn and respond to daily patterns of predation risk

Reassessing the trophic role of reef sharks as apex predators on coral reefs

Apex predators often have strong top-down effects on ecosystem components and are therefore a priority for conservation and management. Due to their large size and conspicuous predatory behaviour, reef sharks are typically assumed to be apex predators, but their functional role is yet to be confirmed. In this study, we used stomach contents and stable isotopes to estimate diet, trophic position and carbon sources for three common species of reef shark (Triaenodon obesus, Carcharhinus melanopterus and C. amblyrhynchos) from the Great Barrier Reef (Australia) and evaluated their assumed functional role as apex predators by qualitative and quantitative comparisons with other sharks and large predatory fishes. We found that reef sharks do not occupy the apex of coral reef food chains, but instead have functional roles similar to those of large predatory fishes such as snappers, emperors and groupers, which are typically regarded as high-level mesopredators. We hypothesise that a degree of functional redundancy exists within this guild of predators, potentially explaining why shark-induced trophic cascades are rare or subtle in coral reef ecosystems. We also found that reef sharks participate in multiple food webs (pelagic and benthic) and are sustained by multiple sources of primary production. We conclude that large conspicuous predators, be they elasmobranchs or any other taxon, should not axiomatically be regarded as apex predators without thorough analysis of their diet. In the case of reef sharks, our dietary analyses suggest they should be reassigned to an alternative trophic group such as high-level mesopredators. This change will facilitate improved understanding of how reef communities function and how removal of predators (e.g., via fishing) might affect ecosystem properties