GSVA: gene set variation analysis for microarray and RNA-seq data

Gene set enrichment (GSE) analysis is a popular framework for condensing information from gene expression profiles into a pathway or signature summary. The strengths of this approach over single gene analysis include noise and dimension reduction, as well as greater biological interpretability. As molecular profiling experiments move beyond simple case-control studies, robust and flexible GSE methodologies are needed that can model pathway activity within highly heterogeneous data sets. To address this challenge, we introduce Gene Set Variation Analysis (GSVA), a GSE method that estimates variation of pathway activity over a sample population in an unsupervised manner. We demonstrate the robustness of GSVA in a comparison with current state of the art sample-wise enrichment methods. Further, we provide examples of its utility in differential pathway activity and survival analysis. Lastly, we show how GSVA works analogously with data from both microarray and RNA-seq experiments. GSVA provides increased power to detect subtle pathway activity changes over a sample population in comparison to corresponding methods. While GSE methods are generally regarded as end points of a bioinformatic analysis, GSVA constitutes a starting point to build pathway-centric models of biology. Moreover, GSVA contributes to the current need of GSE methods for RNA-seq data. GSVA is an open source software package for R which forms part of the Bioconductor project and can be downloaded at http://www.bioconductor.org.

Frontshear and backshear instabilities of the mean longshore current

An analytical model based on Bowen and Holman [1989] is used to prove the existence of instabilities due to the presence of a second extremum of the background vorticity at the front side of the longshore current. The growth rate of the so-called frontshear waves depends primarily upon the frontshear but also upon the backshear and the maximum and the width of the current. Depending on the values of these parameters, either the frontshear or the backshear instabilities may dominate. Both types of waves have a cross-shore extension of the order of the width of the current, but the frontshear modes are localized closer to the coast than are the backshear modes. Moreover, under certain conditions both unstable waves have similar growth rates with close wave numbers and angular frequencies, leading to the possibility of having modulated shear waves in the alongshore direction. Numerical analysis performed on realistic current profiles confirm the behavior anticipated by the analytical model. The theory has been applied to a current profile fitted to data measured during the 1980 Nearshore Sediment Transport Studies experiment at Leadbetter Beach that has an extremum of background vorticity at the front side of the current. In this case and in agreement with field observations, the model predicts instability, whereas the theory based only on backshear instability fai led to do so.