 Methodology
 Open Access
 Open Peer Review
 Published:
Visualising associations between paired ‘omics’ data sets
BioData Mining volume 5, Article number: 19 (2012)
Abstract
Background
Each omics platform is now able to generate a large amount of data. Genomics, proteomics, metabolomics, interactomics are compiled at an ever increasing pace and now form a core part of the fundamental systems biology framework. Recently, several integrative approaches have been proposed to extract meaningful information. However, these approaches lack of visualisation outputs to fully unravel the complex associations between different biological entities.
Results
The multivariate statistical approaches ‘regularized Canonical Correlation Analysis’ and ‘sparse Partial Least Squares regression’ were recently developed to integrate two types of highly dimensional ‘omics’ data and to select relevant information. Using the results of these methods, we propose to revisit few graphical outputs to better understand the relationships between two ‘omics’ data and to better visualise the correlation structure between the different biological entities. These graphical outputs include Correlation Circle plots, Relevance Networks and Clustered Image Maps. We demonstrate the usefulness of such graphical outputs on several biological data sets and further assess their biological relevance using gene ontology analysis.
Conclusions
Such graphical outputs are undoubtedly useful to aid the interpretation of these promising integrative analysis tools and will certainly help in addressing fundamental biological questions and understanding systems as a whole.
Availability
The graphical tools described in this paper are implemented in the freely available R package mixOmics and in its associated web application.
Introduction
‘Omics’ data now form a core part of systems biology by enabling researchers to understand the integrated functions of a living organism. However, the available abundance of such data (genomics, proteomics, metabolomics, interactomics...) is not a guarantee of obtaining useful information in the investigated system if the data are not properly processed and analyzed to highlight this useful information. A major challenge with the integration of omics data is therefore the extraction of discernable biological meaning from multiple omics data.
Recently, several authors have further improved statistical methodologies to integrate two highly dimensional data sets. Such methodologies include regularized and sparse variants of Canonical Correlation Analysis (CCA) [1–5] and Partial Least Squares (PLS) regression [6, 7]  also referred as projectionbased methods. These multivariate approaches aim at unravelling the correlation structure between two sets of data measured on the same samples. In addition, they achieve dimension reduction by summarizing the data into a small number of components or variates, which are linear combinations of the original variables. These exploratory approaches aim at exploiting coexpression between disparate types of biological measures instead of differential expression. The assumption relies on the fact that similar expression patterns across a set of samples are hypothesized to have a functional relationship [8]. In order to better understand the link between the different biological entities from highly dimensional data sets, several clustering techniques have been proposed in the literature. One category of approaches include simple criteria matching, which order the variables according to foldchange or univariate statistical tests for a given threshold. These variables are considered to be ‘clustered’ [9]. Other methods, such as selforganizing maps, use Euclidian distances. However, they are known to encounter difficulties in finding variables “negatively” (oppositely) associated with each other [10]. Another way is to comprehensively compare all variables against each other using a similarity measure, such as Pearson correlation coefficient [10, 11], or mutual information [9]. Once these associations are graphically represented, the aim is to obtain fresh insights into the different biological functional levels, which then act as a foundation for new hypotheses.
So far, most the statistical integrative projectionbased approaches cited above have been limited to numerical results, and little attention has been paid to either the interpretation of the results or the graphical outputs. In this article, we propose to revisit some graphical outputs mostly dedicated to exploratory approaches to highlight associations between two different types of biological entities. We have improved Correlation Circles plots, Relevance Networks and Clustered Image Maps (CIM) to be specifically adapted to the results of our previously published CCA or PLS methods [1, 4, 7]. These graphical outputs are implemented in the R package mixOmics ^{a} that is dedicated to the integrative analysis of ‘omics’ data [12]. For users not familiar with the R programming language, a web application is also available at http://mixomics.qfab.org.
In the following ‘Background’ Section, we first describe the three graphical outputs used in mixOmics to visualise pairwise associations between two types of biological variables. In the ‘Results and discussion’ Section, we assess the relevance of the proposed CIM and Relevance Networks on a simulation study. On two real data sets, we provide a thorough biological interpretation of the results obtained and compare the inferred statistical networks to known biological networks using data knowledge driven analyses. The ‘Methods’ Section describes how to compute the pairwise similarity matrix to construct the graphical representations proposed.
Background
We first briefly introduce PLS and CCA methodologies and their associated variants recently developed for the highly dimensional case. More details about the approaches are given in the ‘Methods’ Section. We review the three main graphical outputs proposed in mixOmics: Correlation Circle plots, Relevance Networks and Clustered Image Maps to visualise pairwise associations between disparate biological entities highlighted by CCA or PLS.
Integrative approaches
The twoblock data matrices to be integrated are denoted X(n×p) and Y(n×q), where p and q are the total number of variables measured on the same n subjects. For example X is a gene expression matrix and Y contains metabolites concentrations, both transcripts and metabolites being measured on the same patients. CCA and PLS search for the largest correlation and covariance respectively between orthogonal components, also called variates, which are linear combinations of the X and Y variables. The number of chosen dimensions or components in CCA or PLS is d, with d≤ min(p,q) for CCA and d≤p for PLS.
In classical CCA and PLS regression, all variables from both data sets are included in the fitted linear combinations or variates. However, in the context of high throughput biological data, the number of variables often exceeds tens of thousands. In this case, linear combinations of the entire set of features make biological interpretability difficult as they contain too many variables to perform further tests or to generate biological hypotheses. Most importantly, the high dimensionality and the insufficient sample size lead to computational problems as CCA requires the computation of the inverse of the covariance matrices of X and Y. To circumvent this problem, regularized CCA (rCCA) has been recently proposed by [1] when dealing with illconditioned covariance matrices by adding a regularization term on their diagonal. Sparse PLS (sPLS) has been recently proposed to perform simultaneous variable selection in the two data sets [4, 7]. sPLS includes Lasso penalization terms on the loading vectors (the vectors which weight are used in the determination of the PLS variates) to shrink some of the coefficients towards zero.
rCCA and sPLS are both implemented in the R package mixOmics[12]. These methodologies require to choose or tune the number of dimensions d, the regularization parameters for rCCA and the number of variables to select in both data sets for sPLS. Guidelines to choose these parameters are discussed in [1, 7]. The biological relevancy of rCCA and sPLS has been recently demonstrated in several biological studies [13–18].
Correlation Circle plots
Correlation Circle plots were primarily used for PCA outputs to visualise the relationship between variates and variables of the same type, where one single omics data set is analysed [19–23]. The use of such a graphical tool was then generalised to represent variables of two different types using statistical integrative approaches such as Canonical Correlation Analysis and Partial Least Squares regression [24].
Although not very well known, this plot is an enlightening tool for data interpretation, as it enables a graphical examination of the relationships between variables and variates. In this plot, the coordinates of the variables are obtained by calculating the correlation between each original variable and their associated component (see Figure 1(a)). Because variables are usually centered and standardized, the correlation between each variable and a component is simply the projection of the variable on the axis defined by the component.
In this plot, variables can be represented as vectors (see Figure 1(b)) and the relationship (correlation) between the two types of variables can be approximated by the inner product between the associated vectors. The inner product is defined as the product of the two vectors lengths and their cosine angle. Thus, the nature of the correlation between two variables can be visualised through the angles between two vectors (Figure 1(b)): if the angle is sharp, the correlation is positive, if the angle is obtuse the correlation is negative and if the angle is right the correlation is null.
The centered and standardized variables are projected onto the space spanned by the two chosen components, inside a circle of radius 1. Thus, from the inner product definition, the longer the distance to the origin, the stronger the relationship between the variables. The variables closely located to the circumference of radius 1 can be directly interpreted, since the closeness on the plane corresponds to the closeness in the ddimensional variables space.
For variables closely located to the origin, it means that some information can be carried on other axes and, it might be necessary to visualise the Correlation Circles plots in the subsequent dimensions (see example in Figure 2).
For the sake of interpretability, variables are not represented as vectors but as the end points of the vectors in mixOmics. Two circles are usually represented, of radii 0.5 and 1, to better visualise the ‘important’ variables. Figure 2 gives an example of the different scenarios that can be encountered when visualising the correlation structure between two data sets. The data come from the simulation study of the following ‘Results and discussion’ Section.
The variables or groups of variables strongly positively correlated are projected closely to each other on the Correlation Circle. This is the case, for instance, with X _{ B } and Y _{ B } on the dimension 2 (Figure 2(ab)) and X _{ C } and ${Y}_{C}^{1}$ on the dimension 3 (Figure 2(b)). When the correlation is strongly negative, the groups of variables are projected at diametrically opposite places on the Correlation Circle. This occurs, for instance, with X _{ A } and Y _{ A } on the dimension 1 (Figure 2(a)) and X _{ C } and ${Y}_{C}^{2}$ on the dimension 3 (Figure 2(b)). The variables or groups of variables that are not correlated are situated 90° one from the other in the circle (for instance, X _{ A } and Y _{ B } (Figure 2(a)) and X _{ C } and Y _{ B } (Figure 2(b))).
Correlation Circle plots were found to supplement pair wise correlation approaches [25]. In the high dimensional case, the interpretation of the correlation structure between variables from two data sets can be difficult, and a threshold can be chosen to remove some weaker associations.
Relevance networks
A conceptually simple approach for modelling netlike correlation structures between two data sets is to use Relevance Networks. This concept was introduced by [10] as a tool to study associations between pair of variables coming from several types of genomic data. This method generates a graph where nodes represent variables, and edges represent variable associations. The Relevance Network is built in a simple manner. First, the correlation matrix is inferred from the data. Second, for every estimated correlation coefficients exceeding (in absolute value) a prespecified threshold between two variables (say 0.6), an edge is drawn between these two variables, otherwise, no edge is drawn and these two variables are considered not associated for this threshold, and the variables/nodes with no link are not represented in the graph.
The construction of biological networks (genegene, proteinprotein, etc.) with direct connections within a variable set is of considerable interest amongst biologists, and has been extensively used in the literature. Therefore, we will not consider this case and focus rather on the representation between variables of two different types. We will thus display Relevance Networks through the use of bipartite graph (or bigraph), where variables/nodes from X can only be connected to variables/nodes from Y.
Instead of computing the Pearson correlation coefficients between each pair of variables as was proposed by [10], bipartite networks are inferred using a pairwise similarity matrix directly obtained from the outputs of the integrative approaches (regularized) CCA and (sparse) PLS. The values in the similarity matrix are computed as the correlation between the two types of projected variables onto the space spanned by the first components retained in the analysis. The values in the similarity matrix can be seen as a robust approximation of the Pearson correlation (see Section ‘Methods’).
The advantage of relevance networks is their ability to simultaneously represent positive and negative correlations, which are missed by methods using Euclidian distances or mutual information. Another advantage is their ability to represent genes in several pathways, and, most importantly for our purpose, to represent correlations across disparate biological measures. One of the main limitation of relevance networks is that it requires extensive computing ressources as mentioned by [26] to compute the comprehensive pairwise associations when the underlying network is fully connected, i.e. when there is an edge between any pair of two types X and Y variables (see also [27] who recently proposed an R package for fast computation of the correlations). In the case of the sparse PLS, extensive computing is not required since the pairwise association are only computed for the variables selected by the approach. A threshold is also proposed to remove some weaker associations.
Since the relevance networks are visual representations of the correlations between variables, one looks for clusters or subnetworks of subsets of variables, where the edge colors indicate the nature of the correlation (positive, negative, strong or weak). Each of these clusters often highlight a specific correlation structure between the features. More details about the relevance networks interpretation can be found in the Section ‘Results and discussion’.
Clustered Image Maps
Clustered Image Maps (CIM), also called ‘clustered correlation’ or ‘heatmaps’ were first introduced by [11, 28, 29] to represent either the expression value of a single data set, or the Pearson correlation between two matched data sets [11, 30]. This type of representation is based on a hierarchical clustering simultaneously operating on the rows and columns of a realvalued similarity matrix. This is graphically represented as a 2dimensional coloured image, where each entry of the matrix is coloured on the basis of its value, and where the rows and columns are reordered according to the hierarchical clustering. Dendrograms (tree diagrams) illustrating the arrangement of the clusters produced by the hiearchical clustering are added to the left (or right) side and to the top (or bottom) of the image. The color in the heatmap indicates the nature of the correlation between subsets of variables (positive, negative, strong or weak), while the dendrogram indicate the proximity between correlated variables. In practice, one looks for well defined large rectangles or squares of the same color corresponding to long branches of the dendrograms. More details about the CIM interpretation can be found in the the Section ‘Results and discussion’.
The similarity matrix represented by the CIM is the same as in the relevance networks described above. CIM is a visualisation tool that complements well the Correlation Circles plots and the Relevance Networks as clusters of subsets of variables of the same type correlated with subsets of variables of the other type can be observed. This complementarity of three graphical outputs is illustrated in the Section ‘Results and discussion’ on the Nutrimouse case study.
Implementation in mixOmics
Correlation Circles plots, Relevance Networks and Clustered Image Maps are implemented in the R package mixOmics[12] to be applied to a variety of integrative approaches implemented in the package, such as rCCA and sPLS methodologies. Full tutorials on how to analyse data sets with different methodologies and how to obtain specific graphical outputs with desired legends and colors are available on the website http://www.math.univtoulouse.fr/~biostat/mixOmics. For users not familiar with the R programming language, an associated web application is available at http://mixomics.qfab.org and provides a Cytoscape plugin to display the Relevance Networks in an attractive manner.
Results and discussion
We investigate the relevance of Correlation Circle plot, Relevance Networks and CIM representations, firstly on a simulated data set to assess if the proposed graphical outputs are able to highlight pairwise association structure between two data sets, and secondly on two biological data sets to assess the biological relevance of such graphical tools.
Simulated data
Data sets
We generated two data sets X and Y with an equal number of 30 observations in each data set. A subset of relevant variables in X were associated with a subset of relevant variables in Y according to the model described below, and the remaining variables were simulated as noise. This simulation study enables to assess if the proposed graphical representations allow differentiate the associated groups of relevant variables from the noisy variables.

The relevant X and Y variables were generated according to a normal distribution with zero mean and covariance matrix Σ defined by :
$$\Sigma =\left[\phantom{\rule{0.3em}{0ex}}\begin{array}{ll}{\Sigma}_{X\phantom{\rule{0.1em}{0ex}}X}& \phantom{\rule{0.1em}{0ex}}{\Sigma}_{{}_{X\phantom{\rule{0.1em}{0ex}}Y}}\\ {\Sigma}_{X\phantom{\rule{0.1em}{0ex}}Y}^{\prime}& \phantom{\rule{0.1em}{0ex}}{\Sigma}_{{}_{Y\phantom{\rule{0.3em}{0ex}}Y}}\end{array}\phantom{\rule{0.3em}{0ex}}\right],\phantom{\rule{1em}{0ex}}\text{with}\phantom{\rule{1em}{0ex}}{\Sigma}_{\mathit{\text{XY}}}=\left[\begin{array}{ll}{A}_{\mathit{\text{XY}}}\phantom{\rule{1em}{0ex}}0& \phantom{\rule{3.7em}{0ex}}0\\ 0\phantom{\rule{2em}{0ex}}{B}_{\mathit{\text{XY}}}& \phantom{\rule{3.7em}{0ex}}0\\ 0\phantom{\rule{2.5em}{0ex}}0& \phantom{\rule{2em}{0ex}}{C}_{\mathit{\text{XY}}}\end{array}\phantom{\rule{2.3em}{0ex}}\right].$$ 
Details about the covariance matrices can be found in Additional file 1.

X contains three independent sets of respectively 10, 10 and 3 crosscorrelated variables: ${X}_{A}=\left[{X}_{A}^{1},\dots ,{X}_{A}^{10}\right]$, ${X}_{B}=\left[{X}_{B}^{1},\dots ,{X}_{B}^{10}\right]$ and ${X}_{C}=\left[{X}_{C}^{1},{X}_{C}^{2},{X}_{C}^{3}\right]$; and Y contains three independent sets of respectively 10, 5 and 2 crosscorrelated variables: ${Y}_{A}=\left[{Y}_{A}^{1},\dots ,{Y}_{A}^{10}\right]$, ${Y}_{B}=\left[{Y}_{B}^{1},\dots ,{Y}_{B}^{5}\right]$ and ${Y}_{C}=\left[{Y}_{C}^{1},{Y}_{C}^{2}\right]$. These groups of variables are associated with each other according to the crosscorrelation matrix Σ _{ X Y }.

The relevant variables in X _{ A } and Y _{ A } were generated with a negatively crosscorrelation varying between −0.93 and −0.51. The variables in X _{ B } and Y _{ B } were generated with a positive crosscorrelation varying between 0.5 and 0.85; and the variables in X _{ C } and Y _{ C } were generated with an absolute crosscorrelation varying between 0.81 and 0.93, X _{ C } is positively correlated with ${Y}_{C}^{1}$ and is negatively correlated with ${Y}_{C}^{2}$.

The irrelevant (noisy) variables were simulated with a normal distribution with zero mean and covariance identity matrices and were added to the sets such that final data set contained 100 variables for X and 50 variables for Y. These variables are independent within the sets X and Y and with each other.
Analysis process
PLS canonical mode (PLScan) was applied to these data sets and the graphical representations Correlation Circle plots, CIM and Relevance Networks resulting from the statistical approach were plotted. The first three dimensions were chosen for these graphical displays (correlation values between latent variables equal to 0.97, 0.94 and 0.95 respectively on each dimension, before decreasing for the following dimensions, Section ‘Methods’ gives the original references to the parameters tuning for the different approaches in mixOmics).
Graphical outputs
Figure 2 displays the corresponding correlation circle plots. Figure 2(a) highlights the strong negative correlation between clusters X _{ A } and Y _{ A } on dimension 1 and the strong positive correlation between clusters X _{ B } and Y _{ B } on dimension 2. Figure 2(b) underlines a strong positive correlation between X _{ C } and ${Y}_{C}^{1}$ and a strong negative correlation between the latter and ${Y}_{C}^{2}$ on dimension 3. For that last case, Figure 2(c) represents the expression profiles of the variables across the 30 individuals and illustrates in a more intuitive manner the nature of the correlation between the variables.
The pairwise similarity matrix was then computed using our proposed method (see Section ‘Methods’) for the first three PLS dimensions in order to display the CIM (Figure 3). The Euclidian distance and the Ward method were used for the hierarchical clustering. In the CIM display, each coloured block represents an association between subsets of the Xvariables and the Yvariables. The green colour indicates that the X and Y clusters are positively correlated (cluster X _{ B } and Y _{ B }, and cluster X _{ C } and ${Y}_{C}^{1}$), and the red colour indicates a negative correlation in the XY cluster (cluster X _{ A } and Y _{ A }, and cluster X _{ C } and ${Y}_{C}^{2}$), whereas yellow indicate weaker correlation values. The dendrograms on the top and the left hand side of the map indicate how the clusters join, the longer the distance, the sharper the boundary between the coloured blocks.
The variables with blank names indicate variables with weak correlations (irrelevant variables). The CIM details in a more comprehensive manner than the correlation circle plots the correlations between all variables.
The Relevance Networks obtained with PLScan are displayed in Figure 4. Similarly to CIM representation, the pairwise similarity matrix was computed for the first three dimensions (see Section ‘Methods’). Three relevant components were obtained setting a threshold to 0.5, linking the corresponding correlated subsets: X _{ A } with Y _{ A }, X _{ B } with Y _{ B } and X _{ C } with Y _{ C }. Note that none of the irrelevant variables were displayed in the network, demonstrating the good ability of the PLS approach to estimate the real simulated correlations.
Comparison with rCCA
The rCCA approach was also applied to these data sets with regularization parameters λ _{1}=0.889 and λ _{2}=0.889 for the first three dimensions (canonical values obtained were of 0.959, 0.925, and 0.881 on each dimension respectively, followed by much lower values). As expected, the graphical outputs were very similar to those with PLScan.
In this simulation setting, accuracy of inferred networks from rCCA and PLScan were similar, while networks inferred with Pearson correlation gave a much higher error rate, especially for a small number of samples (see Additional file 2).
This simulation study shows that Correlation Circle plots, Relevance Networks and CIM are able to highlight the relevant variables amongst the noisy ones and pinpoint the pairwise association structure between the two data sets. In the following, we illustrate the use of such graphical outputs on real data sets and discuss the biological relevancy of the obtained results.
Biological data
Data sets
These data sets are publicly available in the mixOmics package [12] and provide good illustrative examples for this Section. However, much larger biological data sets could be analysed through mixOmics as the integrative approaches rCCA and sPLS have been specifically developed to handle large data sets (several thousands of variables in both data sets).
Nutrimouse data. The data come from a nutrigenomic study [31] in which 40 mice from two genotypes (wildtype and PPAR α / deficient) were fed with five diets with contrasted fatty acid compositions. Oils used for experimental diets preparation were corn and colza oils (50/50) for a reference diet (REF), hydrogenated coconut oil for a saturated fatty acid diet (COC), sunflower oil for an Omega6 fatty acid rich diet (SUN), linseed oil for an Omega3 rich diet (LIN) and corn/colza/enriched fish oils (43/43/14) for the FISH diet. Expression of 120 genes in liver cells were acquired through microarray experiment and concentrations of 21 hepatic fatty acids were measured by gas chromatography. The study therefore includes two data matrices of size (40 × 120) for the gene expression and (40 × 21) for the fatty acids measurements.
Liver toxicity data. The data come from a liver toxicity study [32] in which 64 male rats of the inbred strain Fisher F344/N were exposed to low (50 mg/kg or 150 mg/kg) or to high (1500 mg/kg or 2000 mg/kg) doses of acetaminophen (paracetamol) in a controlled experiment. Necropsies were performed at 6, 18, 24 and 48 hours after exposure and the mRNA from the liver was extracted. Ten clinical chemistry measurements of variables containing markers for liver injury are available for each subject and the serum enzymes levels are numerically measured. The study therefore includes two data matrices of size (64 × 3116) for the gene expression and (64 × 10) for the clinical measurements.
Analysis process
Several methodologies are implemented in mixOmics to integrate data from two different types, to model the relationships between the two types of features in an appropriate manner (see Section ‘Methods’ for a brief description of the models). In the Nutrimouse data, we applied the methodology sPLScan as the aim is to highlight highly correlated subsets of genes and hepatic fatty acids in the two data sets. This study was previously analysed with another approach (rCCA [1]). In this paper, the aim is to illustrate the usefulness of combining the three graphical outputs to interpret such results, and assessing the biological relevance of the Relevance Networks obtained.
In the Liver Toxicity data, we applied the methodology sPLSreg as the aim is to highlight a subset of correlated genes which expression can predict the clinical chemistry measurements [33]. This analysis was performed in a previous paper to demonstrate the numerical good results of the sPLSreg approach but no focus was made on the biological relevance of the results or on the use of variable graphical outputs. In this paper, we focus instead on the biological relevancy of the resulting Relevance Networks. In both studies, using these integrative methodologies and associated graphical outputs, the biological questions we ask are: which subsets of variables from both types are strongly positively or negatively correlated with each other? Do these selected features bring any relevant insight in relation to system under study?
Two parameters need to be tune in sPLS: the number of dimensions and the number of variables to select on each dimension. For both data sets, three dimensions were chosen (see numerical results presented in [1, 7]). To illustrate the use of the proposed graphical outputs, we arbitrarily chose to select 50 transcripts or genes on each dimension. This rather large selection size (150 transcript or genes) is justified by the Gene Ontology (GO) analysis which require a sufficient number of variables to assess their biological relevance. The similarity matrices were computed from the sPLS method on the basis of the selected variables.
To highlight the strongest variable associations only, variables with an association score greater than 0.6 in absolute value were chosen to infer the Relevance Networks. This threshold was arbitrarily chosen in order to obtain biologically interpretable networks that were neither too sparse nor too dense. The obtained networks were then used as an input to Cytoscape [34] for visualization and GeneGo [35] and topGO [36, 37] were used to assess the biological relevancy of the inferred associations between the different types of variables (see Additional file 3 for the R script used and how to export the network to a Cytoscape file format). This analysis is similar to the one performed by [38] who assessed the results of rCCA in a metabolic syndrome study. We then compared the obtained inferred networks to known biological networks through data driven and knowledge driven biological analyses.
Application to Nutrimouse data
Preliminary analysis comparing the different graphical outputs. In order to illustrate the usefulness of the variable graphical outputs in a real case study, we first discuss the outputs obtained on the first two components, where 50 genes were selected on each dimension. The Correlation Circle plot (Figure 5) displays all fatty acids and the genes selected on each component (a 100 in total in this plot). It highlights subsets of variables that are important to define each component. For example C18:2 ω 6, C20:2 ω 6 and C16:0 are the fatty acids which variation mainly participate to the definition of the sPLS component 2 (top and bottom of the yaxis). Similarly, genes such as CAR1, ACOTH, SIAT4C, SR.BI, Ntop are positively correlated to each other, and to the fatty acid C16:1 ω 9 and their variation participate to the definition of the sPLS component 1 (lefthand side of the xaxis).
While the CIM better highlights different clusters of variables and their degree of correlation (indicated by the colour code) than the Correlation Circle plots (Figure 6), the visualisation of the correlation within variables sets is more difficult to observe. For example, the Correlation Circle plot highlights a negative correlation between [C18:2 ω 6, C20:2 ω 6] and C16:0, which is less obvious in the CIM.
Finally, the relevance network representation (Figure 7) adds another layer of information as it allows the visualisation of variable groups in the network. In this case, the network highlights two main subsets of genes and fatty acids (top and bottom) which seem to contain very specific information in each of these groups. This information is slightly suggested on the CIM after a careful interpretation of the dendrograms, but is barely observable in the Correlation Circle plot. This comparison demonstrates the usefulness of such graphical outputs, as well as their complementarity to unravel the complex relationship structure between these different biological features. In the following, we discuss the biological relevance of the subsets of genes and fatty acids highlighted by the full sPLScan analysis (with 3 sPLS components) using GeneGo and topGO. For an easier visualisation, the cytoscape software was used to represent the networks, but the similarity matrix was estimated with mixOmics.
GeneGO analysis. The Relevance Network generated for the Nutrimouse data at a threshold 0.6 highlighted two subsets of fatty acids, and three subsets of genes (Figure 8). Considering first the fatty acids, the yellow group on the lefthand side contained all the ω 6 fatty acids from the data set (C18:2 ω 6, C20:2 ω 6, C20:4 ω 6, C20:3 ω 6, C22:5 ω 6, and C22:4 ω 6). The second group of fatty acids consisted of those in the ω 9, ω 7, and saturated fatty acid groups, along with the two ω 3 fatty acids included in the data set. These groups made sense in the context of lipid biosynthetic pathways – one biosynthetic pathway leads to the production of ω 6 lipids, while the ω 9, ω 7 and saturated lipids are the product of an alternative lipid biosynthetic pathway (orange nodes). The ω 3 group was the exception in our analysis – it was generated by a pathway related to the ω 6 pathway (yellow nodes), but based on the connectivity in our network, these fatty acids partitionned with the ω 7, ω 9 and saturated fatty acid group [39].
The three gene sets defined by network topology were: (1) a set of genes that were negatively correlated with only the ω 6 lipid group; (2) a set of genes that were negatively correlated with the ω 6 group, but largely positively correlated with the other lipid group; and (3) a gene set that was only associated with the second lipid group, with positive correlations to the ω 3, ω 7, ω 9, and saturated fatty acids C14:0 and C16:0, but negatively correlated with the C18:0.
The ω 6 group showed only negative correlations with genes selected by sPLScan. This was consistent with the observations made by [31] that feeding mice a diet rich in ω 6 fatty acids lead to the down regulation of several genes on the array.
The second group of genes contained many targets of PPAR α, a nuclear receptor transcription factor associated with the highlevel regulation lipid metabolism (dark blue nodes). PPAR α targets are expected to be associated with longchain polyunsaturated fatty acids from the ω 3 family, while the final subset of genes involved in lipid biosynthesis is expected to be closely associated with the saturated and monosaturated fatty acids of the ω 7 and ω 9 families. Both of these associations were apparent in the network. An indepth analysis of the Nutrimouse data is behind the scope of this article. The reader can refer to [31, 39] for more details about the underlying biological interpretation.
Application to Liver Toxicity data
Visualization of the association between variables Relevance Networks for the Liver Toxicity data were generated from the results obtained with the sPLSreg method. The selected variables with a pairwise association score greater than 0.6 in absolute value were used as an input to Cytoscape (Figure 9), green (red) edge color represent a positive (negative) correlation. This network contained three groups of clinical chemistry measurements (white nodes) denoted A, B and C and four groups of genes (colored nodes) denoted 1, 2, 3 and 4. Considering first the chemistry measurements, groups A and B only consisted of albumin [ALB] and cholesterol [CHOLE] levels respectively. Group C contained indicators of liver injury (Alanine Aminotransferase [ALT] and Aspartate aminotransferase [AST]), indication of renal injury (urea nitrogen [BUN]), and assessment of cholestasis – bile flow interruption (total bile acids [TBA]).
The four gene subsets defined by network topology (Figure 9) were: group 1: a set of genes that were positively correlated with the cholesterol levels B but negatively correlated with C; group 2, a set of genes that were negatively correlated with ALB levels only (group A); group 3, a set of genes positively correlated with group C and negatively correlated with group B; and group 4 a gene set with only positive correlations with group C.
Biological relevance of the extracted genes. Hierarchical clustering (heatmap) of the biological samples on the extracted genes is displayed in Figure 10. This clustering highlights the groups of rats which were treated with different doses of acetaminophen (also found in [32]). Clusters labelled (coloured at the top of the heatmap) with either no (violet), moderate (cyan) or severe (magenta) necrosis of the centrilobular region of the rat liver were obtained by using the expression values of the genes extracted from the network. Levels of the clinical chemistry measurements on each group of samples are given in Additional file 4. Figure 10 also highlights the differences in gene expression profiles between each gene cluster (coloured in dark brown, brown, orange and yellow at the left side of the heatmap). Gene expression differences are clearly observed between the clusters.
Note that this heatmap represents the relationship between the samples (the rats) and the variables (gene expression) and the color inside the heatmap indicates the expression values of the genes, whereas the CIM represents the relationship between two groups of variables (clinical variables and gene expression) and the color inside the heatmap indicates the correlation between the variables.
The extracted genes were uploaded into topGO [36, 37]. A Gene Ontology (GO) enrichment analysis from the gene list was then performed. GO terms significantly enriched include biological processes related to nitric oxide metabolism and cellular stress responses, including responses to unfolded proteins. The top GO molecular functions enriched in the gene set relate to protein binding, nucleotide binding, and enzyme activity (eg. hydrolase, phosphatase, decarboxylase). Cellular component GO terms enriched in the set mostly relate to very general locations, however both an endopeptidase complex and the peroxisome are also present in the list, reinforcing the association of the selected gene products with proteolysis and the response to stress and unfolded proteins.
The individual gene clusters in the sPLSreg network (Figure 9) may also be examined for GO enrichment, as we have done for the larger cluster 4. For example, while examining the biological process terms associated this cluster, we saw an enrichment for processes involving xenobiotic transport, and interesting functional enrichments such as positive regulation of mesenchymal cell proliferation, a process that was previously observed to occur in other tissues in response to epithelial damage signalling to the underlying mesenchyme to initiate proliferation and tissue remodelling [40], and negative regulation of CREB transcription factor activity, interesting due to the previous association of CREB transcription factor with responses to cytotoxic stress [41, 42], particularly in renal tubular cells [43].
Analysis of the gene list using the GeneGo [35] network analysis algorithm identified a total of 14 networks with a significant enrichment of genes in the Relevance Network. The top five networks were (i) regulation of programmed cell death in response to stress; (ii) cell cycle and regulation of metabolism; (iii) cholesterol and sterol metabolism; (iv) regulation of programmed cell death in response to organic substances; (v) response to stress and presentation of endogenous antigens. A summary of these networks can be found in Additional file 5.
Conclusions
Since several methodologies have been recently proposed to jointly analyse two data sets, the development or the improvement of graphical tools is now crucial to better visualise and understand complex associations between biological entities. In the omics era in particular, the deluge of data can make the interpretation of the results extremely difficult. In our R package mixOmics, we have proposed such graphical tools to ease the interpretation of the implemented integrative methodologies dedicated to the analysis of large biological data sets.
In this paper, we revisited and further developed three types of graphical displays to better understand and interpret the results obtained with CCA and PLS related methods. We thoroughly described how to interpret Correlation Circles plots, which are very insightful graphical outputs to represent the associations between two types of variables. To complement the Correlation Circle plots, we proposed two types of graphical displays: both CIM and Relevance Networks are graphical representations of a pairwise similarity matrix directly estimated from the results of the integrative methodologies implemented in mixOmics. The results obtained on simulated and real data sets illustrated very well the usefulness of these graphical tools to further explore the relationships between two omics data sets. The thorough biological interpretation of the obtained inferred networks using geneGO analysis demonstrated the relevancy of the approach.
Full tutorials are available on http://www.math.univtoulouse.fr/~biostat/mixOmics to use all the methodologies and graphical outputs implemented in mixOmics. An associated web application is also available at http://mixomics.qfab.org and provides a Cytoscape plugin to display the Relevance Networks in an attractive manner.
Methods
We revisit and further develop graphical outputs to visualise correlation structures between two data sets. Correlation Circle plots, CIM and Relevance Networks all use as input byproducts of the integrative approaches implemented in the mixOmics package. Both CIM and the Relevance Networks require the estimation of large scale association or pairwise similarity matrix M as an input. Previously, several similarity measures have been proposed, including Pearson correlation coefficient [10, 29, 44, 45], entropy and mutual information [9]. We propose instead a novel approach to estimate a pairwise similarity matrix using the results of either PLS or CCA approaches.
We briefly describe the PLS and CCA methodologies and associated variants recently developed for the highly dimensional case, more details about these approaches can be found in [1, 7]. We then describe how to estimate the pairwise similarity matrix to construct Relevance Networks and CIM.
CCA and PLS based methods
Notations
We focus on twoblock data matrices denoted X(n×p) and Y(n×q) where the p variables X ^{j} and q variables Y ^{k} are of two types and are measured on the same samples or observations n, for $j=1,\dots ,p$ and $k=1,\dots ,q$. We adopt the following notation: ${M}_{k}^{j}$ represents the element of the k th row and j th column of the matrix M.
CCA
CCA [46] looks for the largest correlation between a linear combination of the variables in the first set X and a linear combination of the variables in the second set Y. The first pair maximizes the correlation ρ _{1}=cor(X a ^{1},Y b ^{1}) subject to var(X a ^{1})=var(Y b ^{1})=1. The subsequent pairs (X a ^{l},Y b ^{l}), (l=2,…, min(p,q)) maximize the residual correlation with the additional requirements that each pair is to be uncorrelated with the previous pairs. In the following, we will refer to a ^{l} and b ^{l} as the canonical loadings (or weights). The resulting variables U ^{l}=X a ^{l} and V ^{l}=Y b ^{l} are called the canonical variates and ρ _{ l } are known as the canonical correlations.
PLS
PLS [47] searches for the largest covariance between linear combinations of the X and the Y variable. The first pair maximises cov(X a ^{1},Y b ^{1}) subject to a ^{1}=b ^{1}=1. Similar to CCA, the subsequent pairs (X a ^{l},Y b ^{l}), (l=2,…,q) maximize the residual covariance with the additional requirements that each pair is to be uncorrelated with the previous pairs. The algorithm of PLS differ from CCA as it is solved in an iterative manner by sequentially decomposing the data matrices. The a ^{l} and b ^{l} are referred to loading vectors (or weights) and the linear combinations U ^{l}=X a ^{l} and V ^{l}=Y b ^{l} as the latent variables (or variates). Several PLS algorithms have been proposed in the literature, for different shapes of data (SIMPLS [48], PLS1 and PLS2 [47], PLSSVD [49]), as well as for different modelling aims (predictive like PLS2, or modelling like PLSmode A, see [2, 23, 50]). In the present paper, we will refer to a PLS approach with two different aims. PLSreg (for PLSregression mode) is used to model an ‘asymmetric’ or unidirectional relationship between the two data sets, i.e. we want to predict the matrix Y with the data X. In that case, the model is Y=A X where A is the matrix of the regression coefficients. PLScan (for PLScanonical mode) is used to model an ‘symmetric’ way and therefore models a bidirectional relationship. In that case, we would like to model A Y=B X where A and B are the matrices of the regression coefficients. PLScan and CCA have very similar purposes.
Regularized and sparse based methods
rCCA. The high dimensionality and the insufficient sample size lead to computational problems as CCA requires the computation of the inverse of matrices X ^{′} X and Y ^{′} Y. To circumvent this problem, [1] developed a regularized (or ridge) extension of CCA (rCCA). rCCA solves the instability of the loadings due to multicollinearity by adding a regularization term on the diagonal of the illconditionned matrices, i.e. the covariance matrices. Thus, highly correlated variables get similar loadings, resulting in a grouping effect. The regularization terms λ _{1} and λ _{2} associated to each data set are chosen by crossvalidation in order to maximize the first canonical correlation.
sPLS. Several sparse PLS have been proposed in the literature to select variables [6, 7]. These approaches introduce l _{1} (Lasso) penalization terms on the loading vectors to shrink some of the coefficients towards zero, thus allowing for simultaneous variables selection in the two data sets. The sparse PLS therefore solves the problem of interpretability by selecting variables from both sets and therefore providing sparse sets of associated variables. In the article, we consider the sparse PLS proposed by [7] since both regression (sPLSreg) and canonical mode (sPLScan, [4]) are available. For practical purposes, the two penalization parameters associated to each data set were replaced by the number of variable to select on each data set and on each sPLS dimension.
Parameters tuning. Both rCCA and sPLS are implemented in mixOmics. These approaches require to choose the number of dimensions d and the regularization/penalization parameters associated to X and Y. For rCCA, the choice of these parameters is based on crossvalidation (see [1] for more details). For sPLS, depending on the modes, several criteria are available to choose these parameters. They are based on the Q ^{2} criterion for the regression mode, or on the maximisation of the correlation for the canonical mode (as was also proposed by [2, 5], see the original article [7] for more details).
Pairwise variable associations for CCA
The similarity measure that we propose to use is analogous to a correlation coefficient. Firstly, similar to a Correlation Circle output, the X ^{j} and Y ^{k} variables are projected onto a low dimensional space. Let d≤ min(p,q) the chosen dimensions to adequately account for the data association, and let Z ^{l}=U ^{l}+V ^{l} the equiangular vector between the canonical variates U ^{l} and V ^{l} (l=1,…,d). The coordinates of the variable X ^{j} and Y ^{k} are obtained by projecting them on the axes defined by Z ^{l}. The projection on the Z axes seems the most natural as X and Y are symmetrically analysed in CCA. Furthermore, [22] showed that the Z variables have the property to be the closest to X and Y, i.e. the sum of their squared multiple correlation coefficients with X and with Y is maximal.
Let ${\mathbf{x}}^{j}={({x}_{1}^{j},\dots ,{x}_{d}^{j})}^{\prime}$ and ${\mathbf{y}}^{k}={({y}_{1}^{k},\dots ,{y}_{d}^{k})}^{\prime}$ the coordinates of the variable X ^{j} and Y ^{k} respectively on the axes defined by Z ^{1},…,Z ^{d}. These coordinates are obtained by computing the scalar inner product ${x}_{l}^{j}=\u3008{X}^{j},{Z}^{l}\u3009$ and ${y}_{l}^{k}=\u3008{Y}^{k},{Z}^{l}\u3009$ (l=1,…,d). As the variables X ^{j} and Y ^{k} are assumed to be of unit variance, the inner product is equal to the correlation between the variables X (or Y) and Z: ${x}_{l}^{j}=\text{cor}({X}^{j},{Z}^{l})$ and ${y}_{l}^{k}=\text{cor}({Y}^{k},{Z}^{l})$.
Then, for any two variables X ^{j} and Y ^{k}, a similarity score can be computed as follows:
where $0\le \left{M}_{j}^{k}\right\le 1$. The matrix M can be factorized as M=x y ^{′} with x and y matrices of order (p×d) and (q×d) respectively. When d=2, M is represented in the Correlation Circle by plotting the rows of x and the rows of y as vectors in a 2dimensional Cartesian coordinate system. Therefore, the inner product of the X ^{j} and Y ^{k} coordinates is an approximation of their association score.
Pairwise variable associations for PLS
For PLSreg, the association score ${M}_{k}^{j}$ between the variables X ^{j} and Y ^{k} can be obtained from an approximation of their correlation coefficient. Let r the rank of the matrix X, according to [51], PLSreg allows for the decomposition of X and Y by:
where ϕ ^{l} and φ ^{l}, are the regression coefficients on the variates U ^{1},…,U ^{r}, and E ^{(r)} is the residual matrix (l=1,…,r). By denoting u _{ l } the standard deviation of U ^{l}, using the orthogonal properties of the variates and the decompositions in (2) and (3), we obtain ${x}_{l}^{j}=\text{cor}({X}^{j},{U}^{l})={u}_{l}{\varphi}_{j}^{l}$ and ${y}_{l}^{k}=\text{cor}({Y}^{k},{U}^{l})={u}_{l}{\phi}_{k}^{l}$. Let d<r the number of components selected to adequately account for the variable association, then for any two variables X ^{j} and Y ^{k}, the similarity score is defined by:
where ${\mathbf{x}}^{j}={({x}_{1}^{j},\dots ,{x}_{d}^{j})}^{\prime}$ and ${\mathbf{y}}^{k}={({y}_{1}^{k},\dots ,{y}_{d}^{k})}^{\prime}$ are the coordinates of the variable X ^{j} and Y ^{k} respectively on the axes defined by U ^{1},…,U ^{d}. When d=2, a Correlation Circle representation is obtained by plotting x ^{j} and x ^{k} as points in a 2dimensional Cartesian coordinate system.
For PLScan, the association score ${M}_{k}^{j}$ is calculated by substituting ${y}_{l}^{k}=\text{cor}({Y}^{k},{V}^{l})$ in (4) for l=1,…,d, as in this case the decomposition of Y is given by:
where φ ^{l} (l=1,…,r), are the regression coefficients on the variates V ^{1},…,V ^{r}. Then,
where ${\sigma}_{l}^{2}$ is the variance of V ^{l}.
Constructing Relevance Networks
Bipartite networks are inferred using the pairwise similarity matrix M defined in (1) and (4) for (r)CCA and (s)PLS results respectively. Entry ${M}_{k}^{j}$ in the matrix M represents the association score between X ^{j} and Y ^{k} variables. Then, by setting a userdefined score threshold, the pairs of variables X ^{j} and Y ^{k} with a $\left{M}_{k}^{j}\right$ value greater than the threshold will be aggregated in the Relevance Network. By changing this threshold, the user can choose to include or exclude relationships in the Relevance Network. This option is proposed in an interactive manner in the mixOmics package [12].
Relevance Networks for (r)CCA assume that the underlying network is fully connected, i.e. that there is an edge between any pair of X and Y variables. For sPLSreg and sPLScan, Relevance Networks are solely represented for the variables selected in the model. In this case, ${M}_{k}^{j}$ pairwise associations are calculated based on the selected variables.
Displaying CIM
CIM or heatmaps were introduced in [11, 29] to represent data resulting from gene expression profiles. This type of representation is based on a hierarchical clustering simultaneously operating on the rows and columns of a realvalued similarity matrix M. The initial matrix is graphically represented as a 2dimensional coloured image, where each entry of the matrix is coloured on the basis of its value, and where the rows and columns are reordered according to a hierarchical clustering. Dendrograms resulting of the clustering are added to the left (or right) side and to the top (or bottom) of the image. With (r)CCA, (s)PLScan and (s)PLSreg, we chose to display CIM based on the pairwise similarity matrix M defined in (1) and in (4).
Endnotes
Authors contributions
IG performed the statistical analysis, the network analysis, wrote the R functions and drafted the manuscript. KALC performed the statistical analysis and drafted the manuscript. MD performed the some of the network analysis. SD participated in the design of the manuscript and helped to draft the manuscript. All authors read and approved the final manuscript.
References
 1.
González I, Déjean S, Martin P, Gonçalves O, Besse P, Baccini A: Highlighting relationships between heteregeneous biological data through graphical displays based on regularized Canonical Correlation Analysis. J Biol Syst. 2009, 17 (2): 173199. 10.1142/S0218339009002831.
 2.
Waaijenborg S, Verselewel de Witt Hamer PC, Zwinderman A: Quantifying the association between gene expressions and dnamarkers by penalized canonical correlation analysis. Stat App Genet Mol Biol. 2008, 7: Article 3
 3.
Parkhomenko E, Tritchler D, Beyene J: Sparse canonical correlation analysis with application to genomic data integration. Stat App Genet Mol Biol. 2009, 8: 134.
 4.
Lê Cao KA, Martin P, RobertGranié C, Besse P: Sparse canonical methods for biological data integration: application to a crossplatform study. BMC Bioinformatics. 2009, 10: 3410.1186/147121051034.
 5.
Witten DM, Tibshirani R, Hastie T: A penalized matrix decomposition, with applications to sparse principal components and canonical correlation analysis. Biostatistics. 2009, 10 (3): 515534. 10.1093/biostatistics/kxp008.
 6.
Chun H, Keles S: Sparse Partial Least Squares Regression with an Application to Genome Scale Transcription Factor Analysis. Technical report, Department of Statistics, University of Wisconsin, Madison, USA 2007
 7.
Lê Cao KA, Rossouw D, RobertGranié C, Besse P: A sparse PLS for variable selection when integrating omics data. Stat App Genet Mol Biol. 2008, 7 (1): 35
 8.
Lee H, Hsu A, Sajdak J, Qin J, Pavlidis P: Coexpression analysis of human genes across many microarray data sets. Genome res. 2004, 14 (6): 10851094. 10.1101/gr.1910904.
 9.
Butte AJ, Kohane IS: Mutual Information Relevance Networks: Functional Genomic Clustering Using Pairwise Entropy Measurement. Genome research, Volume 5. 2000, 415426.
 10.
Butte AJ, Tamayo P, Slonim D, Golub TR, Kohane IS: Discovering functional relationships between RNA expression and chemotherapeutic susceptibility using relevance networks. Proc Nat Acad Sci USA. 2000, 97: 1218212186. 10.1073/pnas.220392197.
 11.
Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome–wide expression patterns. Proc Nat Acad Sci USA. 1998, 95: 1486314868. 10.1073/pnas.95.25.14863.
 12.
Lê Cao KA, González ISD: integrOmics: an R package to unravel relationships between two omics data sets. Bioinformatics. 2009, 25 (21): 28552856. 10.1093/bioinformatics/btp515.
 13.
Bertini I, Calabro A, De Carli V, Luchinat C, Nepi S, Porfirio B, Renzi D, Saccenti E, Tenori L: The metabonomic signature of celiac disease. J proteome res. 2009, 8: 170177. 10.1021/pr800548z.
 14.
Jozefczuk S, Klie S, Catchpole G, Szymanski J, CuadrosInostroza A, Steinhauser D, Selbig J, Willmitzer L: Metabolomic and transcriptomic stress response of Escherichia coli. Mol Syst Biol. 2010,, 6: 364
 15.
Koenig J, Spor A, Scalfone N, Fricker A, Stombaugh J, Knight R, Angenent L, Ley R, Succession of microbial consortia in the developing infant gut microbiome: Proc Nat Acad Sci. 2011, 108 (Supplement 1): 4578
 16.
Morine M, McMonagle J, Toomey S, Reynolds C, Moloney A, Gormley I, Gaora P, Roche H: Bidirectional gene set enrichment and canonical correlation analysis identify key dietsensitive pathways and biomarkers of metabolic syndrome. BMC bioinformatics. 2010, 11: 49910.1186/1471210511499.
 17.
Rose M, Rose T, PariascaTanaka J, Wissuwa M: Revisiting the role of organic acids in the bicarbonate tolerance of zincefficient rice genotypes. Funct Plant Biol. 2011, 38 (6): 493504.
 18.
Hamid J, Meaney C, Crowcroft N, Granerod J, Beyene J: Potential risk factors associated with human encephalitis: application of canonical correlation analysis. BMC Med Res Methodology. 2011, 11: 12010.1186/1471228811120.
 19.
Caillez F, Pages JP: Introduction à l’analyse des données. 1976, Paris, SMASH: Mathématiques et sciences humaines
 20.
van der Burg E, de Leeuw J: Non–linear canonical correlation. Br J Math Stat Psychology. 1983, 36: 5480. 10.1111/j.20448317.1983.tb00765.x.
 21.
van der Geer JP: Introduction to linear multivariate data analysis, Vol. 1. 1986, Leiden, The Netherlands: DSWO Press
 22.
Saporta G: Probabilités analyse des données et statistique. 2006, Paris: Ed. Technip
 23.
Tenenhaus M: La régression PLS: théorie et pratique. 1998, Paris: Ed. Technip
 24.
Mevik BH, Wehrens R: The pls Package: Principal Component and Partial Least Squares Regression in R. J Stat Software. 2007, 18 (2): 124. [http://www.jstatsoft.org/v18/i02]
 25.
Buehler D, Versteegh M, Matson K, Tieleman B: One Problem, Many Solutions: Simple Statistical Approaches Help Unravel the Complexity of the Immune System in an Ecological Context. PloS one. 2011, 6 (4): e1859210.1371/journal.pone.0018592.
 26.
Moriyama M, Hoshida Y, Otsuka M, Nishimura S, Kato N, Goto T, Taniguchi H, Shiratori Y, Seki N, Omata M: Relevance Network between Chemosensitivity and Transcriptome in Human Hepatoma Cells1. Mol Cancer Ther. 2003, 2 (2): 199
 27.
Langfelder P, Horvath S: Fast R Functions for Robust Correlations and Hierarchical Clustering. J Stat Software. 2012, 46 (11): i11
 28.
Weinstein J, Myers T, Buolamwini J, Raghavan K, Van Osdol W, Licht J, Viswanadhan V, Kohn K, Rubinstein L, Koutsoukos A: Predictive statistics and artificial intelligence in the US National Cancer Institute’s drug discovery program for cancer and AIDS. Stem Cells. 1994, 12: 1322. 10.1002/stem.5530120106.
 29.
Weinstein JN, Myers TG, O’Connor PM, Friend SH, Fornace Jr AJ, Kohn KW, Fojo T, Bates SE, Rubinstein LV, Anderson NL, Buolamwini JK, van Osdol WW, Monks AP, Scudiero DA, Sausville EA, Zaharevitz DW, Bunow B, Viswanadhan VN, Johnson GS, Wittes RE, Paull KD: An information–intensive approach to the molecular pharmacology of cancer. Science. 1997, 275: 343349. 10.1126/science.275.5298.343.
 30.
Scherf U, Ross D, Waltham M, Smith L, Lee J, Tanabe L, Kohn K, Reinhold W, Myers T, Andrews D: A gene expression database for the molecular pharmacology of cancer. Nat Genet. 2000, 24 (3): 236244. 10.1038/73439.
 31.
Martin P, Guillou H, Lasserre F, Déjean S, Lan A, Pascussi JM, San Cristobal M, Legrand P, Besse P, Pineau T: Novel aspects of PPARalphamediated regulation of lipid and xenobiotic metabolism revealed through a nutrigenomic study. Hepatology. 2007, 54: 767777.
 32.
Bushel P, Wolfinger RD, Gibson G: Simultaneous clustering of gene expression data with clinical chemistry and pathological evaluations reveals phenotypic prototypes. BMC Syst Biol. 2007, 1: 1510.1186/17520509115.
 33.
Gidskehaug: A framework for significance analysis of gene expression data using dimension reduction methods. BMC Bioinformatics. 2007, 8: 34610.1186/147121058346.
 34.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: A Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res. 2003, 13: 24982504. 10.1101/gr.1239303.
 35.
Ashburner M, Ball C, Blake J, Botstein D, Butler H, Cherry J, Davis A, Dolinski K, Dwight S, Eppig J: Gene Ontology: tool for the unification of biology. Nat Genet. 2000, 25: 2529. 10.1038/75556.
 36.
Alexa A, Rahnenfuhrer J, Lengauer T: Improved scoring of functional groups from gene expression data by decorrelating GO graph structure. Bioinformatics. 2006, 22: 16001607. 10.1093/bioinformatics/btl140.
 37.
Alexa A, Rahnenfuhrer J: topGO: Enrichment analysis for Gene Ontology. 2010, [R package version 2.2.0]
 38.
Morine M, McMonagle J, Toomey S, Reynolds C, Moloney A, Gormley I, Gaora P, Roche H: Bidirectional gene set enrichment and canonical correlation analysis identify key dietsensitive pathways and biomarkers of metabolic syndrome. BMC Bioinformatics. 2010, 11: 49910.1186/1471210511499.
 39.
Guillou H, Zadravec D, Martin PGP, Jacobsson A: The key roles of elongases and desaturases in mammalian fatty acid metabolism: Insights from transgenic mice. Prog Lipid Res. 2010, 49: 186199. 10.1016/j.plipres.2009.12.002.
 40.
Holgate S, Holloway J, Wilson S, Bucchieri F, Puddicombe S, Davies D: Epithelialmesenchymal communication in the pathogenesis of chronic asthma. Proc Am Thoraic Soc. 2004, 1 (2): 9310.1513/pats.2306034.
 41.
Holownia A, Mroz R, Wielgat P, Skiepko A, Sitko E, Jakubow P, Kolodziejczyk A, Braszko J: Propofol protects rat astroglial cells against tertbutyl hydroperoxideinduced cytotoxicity; the effect on histone and cAMPresponseelementbinding protein (CREB) signalling. J Physiol Pharmacol. 2009, 60 (4): 6369.
 42.
Lee B, Cao R, Choi Y, Cho H, Rhee A, Hah C, Hoyt K, Obrietan K: The CREB/CRE transcriptional pathway: protection against oxidative stressmediated neuronal cell death. J Neurochem. 2009, 108 (5): 12511265. 10.1111/j.14714159.2008.05864.x.
 43.
Arany I, Herbert J, Herbert Z, Safirstein R: Restoration of CREB function ameliorates cisplatin cytotoxicity in renal tubular cells. Am J Physiol Renal Physiol. 2008, 294 (3): F57710.1152/ajprenal.00487.2007.
 44.
Scherf U, Ross DT, Waltham M, Smith LH, Lee JK, Tanabe L, Kohn KW, Reinhold WC, Myers TG, Andrews DT, Scudiero DA, Eisen MB, Sausville EA, Pommier Y, Botstein D, Brown PO, Weinstein JN: A Gene Expression Database for the Molecular Pharmacology of Cancer. Nat Genet. 2000, 24: 236244. 10.1038/73439.
 45.
Moriyama M, Hoshida Y, Otsuka M, Nishimura S, Kato N, Goto T, Taniguchi H, Shiratori Y, Seki N, Omata M: Relevance Network between Chemosensitivity and Transcriptome in Human Hepatoma Cells. Mol Cancer Ther. 2003, 2: 199205.
 46.
Hotelling H: Relations between two sets of variates. Biometrika. 1936, 28: 321377.
 47.
Wold H: Estimation of principal components and related models by iterative least squares. Multivariate Anal, Volume 2. Edited by: Krishnaiah P. 1966, New York: Wiley, 391420.
 48.
de Jong: Simpls: An alternative approach to partial least squares regression. Chemom Intell Lab Syst. 1993, 18: 251263. 10.1016/01697439(93)85002X.
 49.
Lorber A, Wangen L, Kowalski B: A theoretical foundation for the PLS algorithm. Chemometrics. 1987, 1: 1931. 10.1002/cem.1180010105.
 50.
Wegelin J: A survey of Partial Least Squares (PLS) methods, with emphasis on the twoblock case. Technical Report 371, Department of Statistics, University of Washington, Seattle. 2000
 51.
Tenenhaus M, Gauchi JP, Ménardo C: Régression PLS et applications. Revue de Statistique Appliquée. 1995, 43: 763.
Acknowledgements
We would like to thank Dr. Pierre Bushel (National Institute of Environmental Health Sciences) for his assistance on the Liver Toxicity study and Dr. Pascal Martin (Institut National de la Recherche Agronomique) for his feedback on the Nutrimouse study. This work was supported, in part, by the Wound Management Innovation CRC to KALC (established and supported under the Australian Government’s Cooperative Research Centres Program).
Author information
Additional information
Competing interests
The authors declare that they have no competing interests.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Received
Accepted
Published
DOI
Keywords
 Partial Little Square
 Canonical Correlation Analysis
 Graphical Output
 Association Score
 Relevance Network