Exploring metabolic pathway disruption in the subchronic phencyclidine model of schizophrenia with the Generalized Singular Value Decomposition
 Xiaolin Xiao^{1},
 Neil Dawson^{2, 3, 4}Email author,
 Lynsey MacIntyre^{3},
 Brian J Morris^{2, 5},
 Judith A Pratt^{2, 3, 4},
 David G Watson^{3, 4} and
 Desmond J Higham^{1, 4}
https://doi.org/10.1186/17520509572
© Xiao et al; licensee BioMed Central Ltd. 2011
Received: 29 October 2010
Accepted: 16 May 2011
Published: 16 May 2011
Abstract
Background
The quantification of experimentallyinduced alterations in biological pathways remains a major challenge in systems biology. One example of this is the quantitative characterization of alterations in defined, established metabolic pathways from complex metabolomic data. At present, the disruption of a given metabolic pathway is inferred from metabolomic data by observing an alteration in the level of one or more individual metabolites present within that pathway. Not only is this approach open to subjectivity, as metabolites participate in multiple pathways, but it also ignores useful information available through the pairwise correlations between metabolites. This extra information may be incorporated using a higherlevel approach that looks for alterations between a pair of correlation networks. In this way experimentallyinduced alterations in metabolic pathways can be quantitatively defined by characterizing group differences in metabolite clustering. Taking this approach increases the objectivity of interpreting alterations in metabolic pathways from metabolomic data.
Results
We present and justify a new technique for comparing pairs of networksin our case these networks are based on the same set of nodes and there are two distinct types of weighted edges. The algorithm is based on the Generalized Singular Value Decomposition (GSVD), which may be regarded as an extension of Principle Components Analysis to the case of two data sets. We show how the GSVD can be interpreted as a technique for reordering the two networks in order to reveal clusters that are exclusive to only one. Here we apply this algorithm to a new set of metabolomic data from the prefrontal cortex (PFC) of a translational model relevant to schizophrenia, rats treated subchronically with the NmethylDAspartic acid (NMDA) receptor antagonist phencyclidine (PCP). This provides us with a means to quantify which predefined metabolic pathways (Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolite pathway database) were altered in the PFC of PCPtreated rats. Several significant changes were discovered, notably: 1) neuroactive ligands active at glutamate and GABA receptors are disrupted in the PFC of PCPtreated animals, 2) glutamate dysfunction in these animals was not limited to compromised glutamatergic neurotransmission but also involves the disruption of metabolic pathways linked to glutamate; and 3) a specific series of purine reactions Xanthine ← Hypoxyanthine ↔ Inosine ← IMP → adenylosuccinate is also disrupted in the PFC of PCPtreated animals.
Conclusions
Network reordering via the GSVD provides a means to discover statistically validated differences in clustering between a pair of networks. In practice this analytical approach, when applied to metabolomic data, allows us to quantify the alterations in metabolic pathways between two experimental groups. With this new computational technique we identified metabolic pathway alterations that are consistent with known results. Furthermore, we discovered disruption in a novel series of purine reactions that may contribute to the PFC dysfunction and cognitive deficits seen in schizophrenia.
Keywords
Background
Background in neuroscience and metabolomics
Schizophrenia is characterized by deficits in cognition known to be dependent upon the functional integrity of the prefrontal cortex (PFC). Furthermore, compromised PFC function in schizophrenia is supported by a multitude of neuroimaging studies reporting hypometabolism ('hypofrontality'), as evidenced by decreased blood flow or glucose utilization [1, 2]. While the pathophysiological basis of PFC dysfunction in schizophrenia is not completely understood, a central role for NMDA receptor hypofunction is widely supported. For example, subchronic exposure to the NMDA receptor antagonist phencyclidine (PCP) induces cognitive deficits and a 'hypofrontality' which directly parallels that seen in schizophrenia [3–5]. Furthermore, subchronic PCP exposure induces alterations in GABAergic cell markers and 5HT receptor expression in the PFC similar to those seen in this disorder [3, 6, 7]. While this evidence places NMDA receptor hypofunction central to the pathophysiology of PFC dysfunction in schizophrenia, the mechanisms through which NMDA hypofunction promotes PFC dysfunction are poorly understood.
Metabolomics is the comprehensive analysis of small molecule metabolites in biological systems [8]. It involves the study of the metabolome which is defined as all of the small molecular weight compounds within a sample that are required for metabolism, whose roles include growth and functionality [9–11]. Sample sources include bacteria, parasites, animals and humans and sample types can include biofluids, cells or tissue extracts. Metabolomics can be utilized as a tool for the characterization and quantification of all of the metabolites in a biological system. Its applications include profiling disease biomarkers [12, 13], monitoring disease progression [14], investigating xenobiotic metabolism [15], investigating druginduced toxicity [16, 17] and investigating metabolism in genetically modified animals [18]. Mass spectrometry (MS) has been employed extensively as an analytical platform for metabolomics studies [19–21]. The popularity of this approach has increased over the last decade, in part due to the advent of high resolution Fourier transform mass spectrometers which offer improved reproducibility, accuracy and sensitivity. This makes mass spectrometry suitable for high throughput metabolomics studies [22]. In addition, the Orbitrap mass spectrometers that are now available offer similar performance to FTMS systems without the need for a high strength magnetic field [23]. HILIC chromatography has been utilized as a separation technique prior to MS detection of polar metabolites in aqueous biofluids such as urine, serum and plasma [24–30].
Additionally, it has also been used for the detection of multiple neurotransmitters in primate cerebral cortex [31]. HILIC chromatography has been chosen for metabolomic studies as it is useful for the analysis of highly polar metabolites which are poorly retained on reverse phase columns [32]. Detailed reviews of the principles and applications of HILIC have been previously outlined [25, 33]. Here, HILICchromatography is utilized in combination with an LTQOrbitrap for metabolic profiling of metabolite extracts from the PFC of control and PCPtreated rats.
Metabolomics represents a robust approach through which alterations in diverse metabolic pathways may be determined at a biological systems level. In this way a metabolomics approach may prove useful in further elucidating the pathophysiological mechanisms contributing to PFC dysfunction in schizophrenia. Furthermore, this approach may also allow for the identification of PFC metabolic biomarkers for the cognitive deficits in this disorder. While the metabolomics approach can provide a rich and comprehensive set of data, the appropriate quantitative analysis of this data has not been adequately developed. In particular, the identification of statistical differences in metabolic pathways between experimental groups rather than the identification of statistical differences in individual metabolites alone represents a major challenge to quantitatively identifying metabolic alterations at a systems level from metabolomic data. One method through which statistical differences in metabolic pathways can be identified from metabolomic data involves the representation of this data as a large, complex network of nodes (single metabolites) connected by realvalue edges (the correlation coefficient between two metabolites). This form of representation has high face validity as the relationship between two metabolites, in a given pathway, is governed by a single or series of enzymatic reactions that can be viewed as being represented by the correlation between the concentrations of the two metabolites. Another advantage is that metabolomic data consist of a range of metabolites detected in both of the experimental groups of interest meaning that these data can be expressed as two complex networks based upon the same set of nodes. This data structure is amenable to analysis through the application of the Generalized Singular Value Decomposition (GSVD) algorithm.
Background in network science and spectral methods
Large, complex interaction networks arise across many applications in science and technology [34–36]. Spectral methods, based on information computed from eigenvectors or singular vectors, have been used successfully to reveal fundamental network properties. For example, we may wish to cluster objects into groups [37], put objects into order [38] or discover specific patterns of connectivity within subgroups [39–42]. In this work, we look at the case where two interaction data sets are available and the aim is discover differences between the two sets in the form of mutually exclusive clusters. For example, a given group of biologically defined entities, such as genes, proteins, metabolites or brain regions, may contain a subgroup that behaves in a coordinated manner under one condition, or in one organism, but not in anotherthe network with respect to one type of interaction contains a cluster that is not present in the other. We will show that the Generalized Singular Value Decomposition, which is becoming more widely used in computational biology [43, 44] can be justified as the basis of a network reordering approach. We also consider how to quantify the statistical significance of network patterns that are uncovered.
Overall, this work develops and applies a novel algorithm in network science and shows that it reveals meaningful insights when applied to cuttingedge metabolomic data.
Results
Derivation of new algorithm
for x ∈ ℝ ^{ N } . Here ·_{2} denotes the Euclidean norm and is one way to generalize the concept of outdegree to the case of a weighted network. Suppose we wish to split the nodes into two groups such that nodes within each group are wellconnected but nodes across different groups are poorly connected. We could use an indicator vector x ∈ ℝ ^{ N } to denote such a partition, with x_{ s } = 1 if node s is placed in group 1 and x_{ s } = 1 if node s is placed in group 2.
 1.
1. a_{ ik }a_{ il } large and positive ⇒ try to choose x_{ k }x_{ l } = +1,
 2.
2. a_{ ik }a_{ il } small or negative ⇒ try to choose x_{ k }x_{ l } = 1.
Returning to the righthand side of (1), we see that is independent of the choice of indicator vector, and gives a measure of how successfully we have incorporated the (possibly conflicting) desiderata in points 1 and 2 over all pairs k, l and third parties i. So we could judge the quality of an indicator vector by its ability to produce a large value of , provided other constraints, such as balanced group sizes, were satisfied.
is a good basis for picking out strong clusters in A that are not present in B.
where U ∈ ℝ^{m×m}and V ∈ ℝ^{p×p}are both orthogonal, X ∈ ℝ^{n×n}is invertible, C ∈ ℝ^{m×n}and S ∈ ℝ^{p×n}are diagonal with nonnegative entries such that C = diag(c_{1}, c_{2},..,c_{ n } ) and S = diag(s_{1}, s_{2},..., s_{ q } ) with q = min(p, n), and 0 ≤ c_{1} ≤ c_{2} ≤ ··· ≤ c_{ n } and s_{1} ≥ s_{2} ≥ ··· ≥ s_{ q } ≥ 0 [49]. The ratios λ_{ i } = c_{ i } /s_{ i } are the generalized singular values of A and B.
A key property of the GSVD is that the columns of X are stationary points of the function f :ℝ ^{ n } ↦ ℝ given by f(x) =  Ax _{2} /Bx _{2}, with the generalized singular values λ_{ i } giving the corresponding stationary values. Hence, we may tackle the problem (3) through the GSVD. Columns 1, 2, 3,... of X are candidates for finding good clusters in B that are poor clusters in A and, analogously, columns N, N  1, N  2,... of X are candidates for finding good clusters in A that are poor clusters in B.
In this way, the existence or lack of clusters in each network becomes apparent from inspection of the heat map of the reordered matrix. This is the approach that we use. We will also show that pvalues can be computed to quantify the statistical significance of the results. The issue of fully automating the choice of cluster size is left as future work.
A variant of the algorithm
which can be solved through the standard SVD.
It is known from spectral graph theory that the dominant singular vectors give good directions in which to look for clusters [37, 50]. Inverting the weight matrix reverses their importance (the singular valueσ becomes σ ^{1}) and hence a spectral clustering approach applied to A^{1} will typically find the opposite of good clusterspoorly connected nodes will be grouped together [51]. So, intuitively, forming AB^{1} in (5) should produce a data matrix for which the SVD approach finds good clusters for A and poor clusters for B. Analogously, the opposite holds for BA^{1}.
We may interpret (6) from the point of view that making B^{1}x large encourages poor clusters for B, while making A^{1}x small encourages good clusters for A. In this case, we would base our algorithm on the GSVD of A^{1} and B^{1}.
Then we may appeal to the arguments given previously and use columns from the inverse of the third factor in the GSVD as the basis for reordering. With this approach we use columns of X ^{T}rather than columns of X. We emphasize that although this heuristic derivation used an assumption that A and B are invertible, the GSVD, and hence the final algorithm, applies in the noninvertible case. Also, the algorithms that we use do not require the computation of matrix inverses.
In tests on both synthetic and real network pairs, we found that this version of the algorithm was more effective, [52]. Hence, in this work we focus on the approach of reordering networks pairs via columns of X ^{T}. In summary, the first few columns of X ^{T}should give orderings that favor clusters in B rather than A and vice versa for the final few columns. In our computational examples, we used the gsvd routine built in to MATLAB http://www.mathworks.com/.
Synthetic test on binary networks
Cluster validation
Suppose we find τ nodes giving a good cluster s for B but a poor cluster for A when the graphs are reordered by column v from X ^{T}. Is this type of substructure likely to arise "by chance"? The following general approach can be used in order to determine a p value, where we will regard a value below 0.05 as indicating statistical significance.
Initialization: Compute a measure of cluster quality, c(A, B), for the promising substructure consisting of those τ nodes in networks A and B reordered by column v.
Step 1: Randomize the networks and obtain new data sets and .
Step 2: Compute the GSVD for the randomized networks and and obtain a matrix .
Step 3: Compute the measure for the τ node 'cluster' in and reordered by column v from .
p value After performing M loops over Steps 1 to 3, compute a p value as the proportion of samples that exceed c(A, B).
Here, E(s) represents the actual number of edges in the object block s, and s is the maximum possible number of edges.
Here, w(s) denotes the average weight in block s. We note the denominator s cancels when ratios are computed in the pvalue algorithm.
In Figure 3, we see that eight nodes 7,9,10,15,14,11,6,13 form a cluster in A, but not in B, when the synthetic data is reordered with the final column of X ^{T}. Applying the procedure above, using permutation to randomize the networks M = 1000 times as described below, we obtained a pvalue of 0.007. Applying the same procedure, we also obtained a pvalue of 0.029 for the first 6 nodes 18, 20, 16, 15, 19, 17 when the synthetic data is reordered with the first column of X ^{T}, which visually form a cluster in B, but not in A. These pvalues (< 0.05) both indicate that the results are statistically significant. As a further test, we arbitrarily selected the subnetworks of A and B composed of nodes 2,4,12,16,1,3,18, which correspond to the 12th to 18th components of the sorted final column from X ^{T}. In this case, we would not expect to find a significant result. This is reflected in the large pvalue of 0.844. In more exhaustive experiments, three randomization methods were tested [52]:
• ErdösRényi: generate a classical random graph with the appropriate number of edges.
• Redistribution: redistribute the entries in each row and each column of A, and perform the same operations on B.
• Permutation: reorder the nodes in A and B and choose the first τ nodes in this new ordering. In this case, recomputation of the GSVD in Step 2 is not necessary, due to the permutation invariance of the factorization.
Of those three approaches, ErdösRényi may be the most commonly used method to randomize a binary network, whereas permutation extends most naturally to the case of weighted edges, so we used permutation in the test shown here. We also tested another simple cluster quality measure which is the ratio of density of edges within the cluster in one graph and that in the other graph.
These variations were studied within this general methodology on both real and synthetic data sets [52]. In all cases, comparable p values were produced.
Synthetic test on correlation networks
Having tested the algorithm on binary networks, we now consider the case where weighted edges arise as correlation coefficients.

D_{ a }∈ ℝ^{20×50}: the first 5 rows are linear combinations of v^{[1]}, v^{[2]}, v^{[3]}, v^{[4]}, v^{[5]}, v^{[6]} and v^{[7]}. Rows 6 to 15 are combinations of v^{[7]} and v^{[8]}. The remaining rows (rows 16 to 20) are Gaussian pseudorandom numbers.

D_{ b }∈ ℝ^{20×50}: the first 5 rows are linear combinations of v^{[1]}, v^{[2]}, v^{[3]}, v^{[4]}, v^{[5]}, v^{[6]} and v^{[7]}. Rows 6 to 14 are Gaussian pseudorandom numbers. The remaining rows (rows 15 to 20) are combinations of v^{[4]} and v^{[9]}.
Building up the rows from the underlying signals in this manner allowed us to construct the correlation patterns seen in Figure 4.
In summary, this additional synthetic test illustrates that our GSVD based algorithm can be extended to reveal the pattern difference between two relative correlation matrices in terms of clustering.
Quantitative determination of metabolic pathways disrupted in the prefrontal cortex of PCPtreated animals
PCPinduced alterations in PFC metabolite levels as determined by SIEVE analysis
Formula  Metabolite  Metaboite KEGG ID  KEGG Pathways  P value  Ratio 

C _{9} H _{11} NO _{3}  LTyrosine  c00082  ko00350, ko00360, ko00400  0.001  0.584 
C _{10} H _{17} N _{3} O _{6}  gamma Glutamylglutamine  NA  NA  0.007  0.673 
C _{6} H _{13} N _{3} O _{3}  LCitrulline  c00327  ko00330  0.007  0.709 
C _{3} H _{7} NO _{2} S  LCysteine  c00097  ko00260, ko00270, ko00430, ko00480, ko00730, ko00770, ko00920  0.012  0.445 
C _{8} H _{9} NO  2Phenylacetamide  c02505  ko00360  0.015  0.561 
C _{9} H _{8} O _{3}  Phenylpyruvate  c00166  ko00360, ko00400  0.016  0.57 
C _{4} H _{6} O _{2}  2,3Butanedione  c00741  map00650  0.017  0.786 
C _{4} H _{5} N _{3} O  Cytosine  c00380  ko00240  0.019  0.665 
C _{04} H _{9} NO _{2}  GABA  c00334  ko00250, ko00330, ko00410, ap00650, ko04080  0.021  0.804 
C _{9} H _{17} NO _{4}  OAcetylcarnitine  c02571  ko00250  0.022  2.649 
C _{14} H _{18} N _{5} O _{11} P  Adenylosuccinate  c03794  ko00230, ko00250  0.029  3.276 
C _{5} H _{5} N _{5} O  Guanine  c00242  ko00230  0.035  0.593 
C _{7} H _{16} NO _{3}  Carnitine  c00487  ko00310  0.037  0.819 
In the context of this study the aim of applying the GSVD algorithm to metabolomic data from control and PCPtreated animals was to quantitatively determine which predefined metabolic pathways were altered in PCPtreated animals. The intermetabolite Pearson's correlation coefficient (partial correlation) was used as the metric of the functional association between each pair of metabolites and was generated from the metabolite peak intensities, as determined by Liquid Chromatography Mass Spectrometry (LCMS), across all animals within the same experimental group (i.e. either control or PCPtreated). These correlations were Fisher transformed to give the correlation data a normal distribution. This resulted in a pair of symmetric, square, realvalued {98 × 98} partial correlation matrices (Control animals: Additional File 2 PCPtreated animals: Additional File 3). Each withingroup matrix represents the specific association strength between each of the 9506 possible pairs of metabolites in that experimental group. In the simplest biological case the correlation coefficient between two metabolites (nodes) in the matrix represents the series of enzymatic reactions responsible for converting one metabolite into another. However, it should be noted that this simple interpretation does not account for the complex relationships that may influence the correlation between two metabolites, such as the involvement of metabolites in alternative, often parallel, metabolic pathways. There are important limitations that must be recognized when modeling metabolomic data as a complex network of interactions between metabolites (as defined by the correlation that exists between them) such as the potential for correlations to exist between metabolites that are not biologically relevant. The impact of such erroneous associations on the interpretation of the data as outlined in this paper will be limited by the approach of characterizing alterations at the level of metabolic pathways, involving multiple metabolites (the approach taken in this study), rather than considering the disruption of single correlation coefficient between two metabolites.
Metabolite identities and their relevant KEGG pathways in the top cluster of Figure 11
Formula  Metabolite  Metaboite KEGG ID  KEGG Pathways  Pvalue  Ratio 

C _{5} H _{10} N _{2} O _{3}  LGlutamine  c00064  Ko00230, ko00240, ko00250, ko00330  0.522  0.959 
H _{3} PO _{4}  Phosphoric acid  c00009  ko00190  0.254  0.915 
C _{5} H _{7} NO _{3}  1Pyrroline4hydroxy2 carboxylate  c04282  ko00330  0.781  0.981 
C _{4} H _{9} N _{3} O _{2}  Creatine  c00300  ko00330, ko00260  0.551  0.953 
C _{ 4 } H _{ 9 } NO _{ 2 }  GABA  c00334  ko00250, ko00330, ko00410, ko04080, map00650  0.021  0.804 
C _{4} H _{7} NO _{4}  LAspartate  c00049  ko00250, ko00260, ko00270, map00300, ko00330, ko00340, ko00410, ko00760, ko00770, ko04080  0.319  0.916 
C _{4} H _{7} NO _{2}  1Aminocyclopropane1carboxylate  c01234  ko00270, ko00640  0.590  0.951 
C _{ 5 } H _{ 5 } N _{ 5 } O  Guanine  c00242  ko00230  0.035  0.593 
C _{5} H _{9} NO _{4}  Glutamate  c00025  ko00250, ko00330, ko00340, ko00471, ko04080, ko00480, map00650  0.845  0.985 
C _{4} H _{7} NO  Hydroxymethylpropanitrile  NA  NA  0.098  0.842 
C_{6}H_{6}N_{2}O  Nicotinamide  c00153  ko00760  0.440  0.917 
C _{ 4 } H _{ 6 } O _{ 2 }  2,3Butanedione  c00741  map00650  0.017  0.786 
C _{6} H _{12} O _{4}  Pantoate  c00552  ko00770  0.722  0.963 
C _{15} H _{23} N _{5} O _{14} P _{2}  ADPribose  c00301  ko00230  0.058  677.029 
C _{3} H _{7} NO _{3}  LSerine  c00065  ko00260, ko00270, ko00600, ko00920, ko00680  0.316  0.856 
C _{ 4 } H _{ 5 } N _{ 3 } O  Cytosine  c00380  ko00240  0.019  0.665 
C _{2} H _{7} NO _{3} S  Taurine  c00245  ko00430, ko04080  0.936  0.995 
C _{4} H _{5} NO _{3}  Maleamate  c01596  ko00760  0.372  0.927 
C _{2} H _{8} NO _{4} P  Ethanolamine phosphate  c00346  ko00260, ko00564, ko00600  0.373  0.889 
Unknown ID  NA  NA  0.271  1.395  
C _{5} H _{11} NO _{3}  Hydroxyvaline  NA  NA  0.585  0.946 
C _{ 6 } H _{ 13 } N _{ 3 } O _{ 3 }  LCitrulline  c00327  ko00330  0.007  0.709 
Metabolite identities and their relevant KEGG pathways in the bottom cluster of Figure 11
Formula  Metabolite  Metaboite KEGG ID  KEGG Pathways  Pvalue  Ratio 

C _{5} H _{4} N _{4} O _{2}  Xanthine  c00385  ko00230  0.339  0.508 
C _{10} H _{16} N _{2} O _{7}  Gamma glutamylglutamic acid  NA  NA  0.143  0.54 
C _{14} H _{26} O _{2}  Myristoleic acid  c08322  NA  0.689  0.623 
C _{5} H _{4} N _{4} O  Hypoxanthine  c00262  ko00230  0.115  0.569 
C _{17} H _{37} NO _{2}  Heptadecasphinganine  NA  NA  0.733  0.769 
C _{10} H _{13} N _{4} O _{8} P  Inosine monophosphate  c00130  ko00230  0.461  0.73 
C _{10} H _{17} N _{3} O _{6}  Peptide fragment (ArgArgGln)  NA  NA  0.775  1.183 
C _{6} H _{15} NO _{3}  Triethanolamine  c06771  ko00564  0.691  1.207 
C _{9} H _{14} N _{4} O _{3}  Carnosine  c00386  ko00340, ko00410  0.872  1.128 
C _{10} H _{12} N _{4} O _{5}  Inosine  c00294  ko00230  0.090  0.6 
C _{15} H _{12} O _{5}  Narigenin  c00509  NA  0.196  0.862 
C _{10} H _{17} N _{3} O _{6}  gamma Glutamylglutamine  NA  NA  0.007  0.673 
C _{26} H _{42} N _{7} O _{20} P _{3} S  2HydroxyglutarylCoA  c03058  map00650  0.179  0.715 
C _{31} H _{54} N _{7} O _{17} P _{3} S  DecanoylCoA  c05274  ko00071  0.410  1.312 
C _{25} H _{44} NO _{7} P  2 Aminoethylphosphocholate  c05683  ko00440  0.243  0.662 
C _{22} H _{26} O _{6}  Eudesmin  NA  NA  0.084  0.493 
C _{3} H _{7} NO _{2} S  LCysteine  c00097  ko00260, ko00270, ko00430, ko00480, ko00730, ko00770, map00920  0.012  0.445 
C_{3}H_{7}O_{6}P )  Glycerone phosphate  c00111  ko00010, ko00051, ko00052, ko00561, ko00562, ko00564, ko00620  0.063  0.381 
Hypergeometric probability of KEGG defined metabolic pathways in the top cluster in Figure 11
KEGG Pathway Identity  KEGG Pathway  Number of metabolites in cluster(A)  Total number of pathway metabolites detected (B)  Hypergeometric Probability (P(X) ≥ k) 

ko00190  Oxidative phosphorylation  1  1  0.224 
ko00230  Purine metabolism  3  13  0.598 
ko00240  Pyrimidine metabolism  2  6  0.406 
ko00250  Alanine, Aspartate and Glutamate metabolism  4  7  0.043 
ko00260  Glycine, Serine and Threonine metabolism  4  7  0.043 
ko00270  Cysteine and Methionine metabolism  3  7  0.186 
map00300  Lysine biosynthesis  1  3  0.538 
ko00330  Arginine and Proline metabolism  7  10  0.001 
ko00340  Histidine metabolism  2  5  0.312 
ko00410  betaAlanine metabolism  2  5  0.312 
ko00430  Taurine and Hypotaurine metabolism  1  3  0.538 
ko00471  Dglutamine and Dglutamate metabolism  1  1  0.224 
ko00480  Glutathione metabolism  1  5  0.728 
ko00564  Glycerophospholipid metabolism  1  11  0.949 
ko00600  Sphingolipid metabolism  2  3  0.126 
ko00640  Propanoate metabolism  1  2  0.400 
map00650  Butanoate metabolism  3  4  0.034 
ko00680  Methane metabolism  1  1  0.224 
ko00760  Nicotinate and Nicotinamide metabolism  3  4  0.034 
ko00770  Pantothenate and CoA biosynthesis  2  5  0.312 
ko00920  Sulphur metabolism  1  3  0.538 
ko04080  Neuroactive ligandreceptor interaction  4  7  0.043 
Hypergeometric probability of KEGG defined metabolic pathways in bottom cluster in Figure 11
KEGG Pathway Identity  KEGG Pathway  Number of metabolites in cluster(A)  Total number of pathway metabolites detected (B)  Hypergeometric Probability (P(X) ≥ k) 

ko00010  Glycolysis/Gluconeogenesis  1  1  0.184 
ko00051  Fructose and Mannose metabolism  1  1  0.184 
ko00052  Galactose metabolism  1  1  0.184 
ko00071  Fatty acid metabolism  1  1  0.184 
ko00230  Purine metabolism  4  13  0.191 
ko00260  Glycine, Serine and Threonine metabolism  1  7  0.770 
ko00270  Cysteine and Methionine metabolism  1  7  0.770 
ko00340  Histidine metabolism  1  5  0.646 
ko00410  betaAlanine metabolism  1  5  0.646 
ko00430  Taurine and Hypotaurine metabolism  1  3  0.460 
ko00440  Phosphonate and Phosphinate metabolism  1  2  0.335 
ko00480  Glutathione metabolism  1  5  0.646 
ko00561  Glycerolipid metabolism  1  2  0.335 
ko00562  Inositol Phosphate metabolism  1  2  0.335 
ko00564  Glycerphopholipid metabolism  2  11  0.642 
ko00620  Pyruvate metabolism  1  2  0.335 
map00650  Butanoate metabolism  1  4  0.562 
ko00730  Thiamine metabolism  1  1  0.184 
ko00770  Pantothenate and CoA biosynthesis  1  5  0.646 
map00920  Sulphur metabolism  1  3  0.460 
Discussion
Through its application to metabolomic data we have clearly demonstrated the added value that can be gained from applying the GSVD algorithm to two sets of complex, network data based upon the same set of nodes. In particular, we have demonstrated that the combined application of the GSVD algorithm with hypergeometric probability analysis provides an analytical framework by which statistical alterations in predefined metabolic pathways between experimental groups can be defined from complex metabolomic data. There is a great unmet need for this type of analytical approach in metabolomics, as well as in the other omics fields (e.g. transcriptomics), which allows the quantification of alterations at the biological systems (pathways) level rather than simply identifying significant alterations of discrete measures (i.e. single metabolites).
Through the application of this analytical approach we identified statistically significant alterations in specific, predefined metabolic pathways (KEGG database pathways) that may contribute to PFC dysfunction in PCPtreated animals, and so in schizophrenia. This included the disruption of the (1) Alanine, Aspartate and Glutamate [ko00250], (2) Arginine and Proline [ko00330], (3) Butanoate [ko00650], (4) Nicotinate and Nicotinamide [ko00760], (5) Glycine, Serine and Threonine metabolic pathways as well as an imbalance in (6) metabolites active as neurotransmitter ligands [ko04080]. The disruption of metabolic pathways involving glutamate in the PFC of PCPtreated rats seems particularly pertinent given the reported alterations in extracellular glutamate availability in the PFC following repeated PCP treatment [53] and the central hypothesis of hypofunctional glutamatergic PFC neurotransmission in schizophrenia [54, 55]. In addition to altered glutamate metabolism there was also evidence to support an imbalance in multiple metabolites known to be active at glutamate receptors. This included an imbalance in the relationship between glutamate, Laspartate and Tauring (Table 2) which are all known to be active at glutamate receptors. Furthermore, evidence for the disruption of glycine, serine and threonine metabolism may suggest that glycine and serine activity as coagonists at the NMDA receptors may be disrupted in the PFC of PCPtreated animals. However, it is important to note that we failed to detect glycine levels in our samples and serine levels appear to be overtly unchanged. The possibility of altered glycine levels in the PFC of PCPtreated rats warrants further investigation given the ability of glycine and NMDA receptor glycine site agonists to reverse subchronic PCPinduced alterations in PFC dopaminergic neurotransmission [56, 57], which may be central to the impact of subchronic PCP treatment on cognition. Altered glycine, serine and threonine metabolism in the PFC of PCPtreated animals is also consistent with the hypothesis that glycine and serine represent potential therapeutic targets for the treatment of schizophrenia [58]. In addition, we found evidence to suggest that GABA neurotransmission was also significantly decreased in the PFC of PCPtreated rats, which may relate to the compromised integrity of GABAergic interneurones in these animals [3, 6], which closely resemble the GABAergic interneuron alterations seen in schizophrenia. The imbalance in glutamate, glutamine and GABA levels identified in the PFC of PCPtreated rats may directly contribute to the hypofrontality (glucose hypometabolism) seen in these animals, as detected using the ^{14}C2deoxyglucose imaging technique [4], as all of these metabolites are intimately linked through metabolic pathways and have a central role in regulating the coupling of neuronal activity to cerebral glucose metabolism [59, 60].
Our results also suggest that glutamatergic dysfunction in the PFC of PCPtreated rats is not limited to the disruption of glutamatergic neurotransmission but also involves the disruption of the metabolic pathways in which glutatmate is involved. For example, altered glutamate metabolism may directly contribute to the disruption of the ArginineProline metabolic pathway in the PFC of PCPtreated animals. The significant disruption of the Arginine pathway in PCPtreated animals suggests that prolonged NMDA receptor hypofunction may result in disrupted nitric oxide (NO) signalling in the PFC. There is increasing evidence that NO signalling is directly linked to NMDA receptor activity through regulation of the enzyme nitric oxide synthase (NOS) [61] and that NO signaling contributes to the deficits in cognition that arise from acute NMDA receptor blockade [62, 63]. The finding that Citrulline levels, a metabolite in the ArginineProline pathway, are significantly decreased in the PFC of the PCPtreated rats in this study further supports the suggestion that NOS activity is altered in the PFC of these animals, as this metabolite is formed by NOS when it releases NO from Larginine. This suggests that NMDA receptor hypofunction may underlie the decreased NOS activity and protein expression levels reported in the PFC of schizophrenia patients [64] and may contribute to the cognitive deficits seen in this disorder.
This result suggests that the activity of adenylosuccinate synthase (ADSS), the enzyme responsible for the conversion of IMP to adenylosuccinate, may be significantly increased in the PFC of PCPtreated animals. An increase in the functional activity of this enzyme could result in both the increased level of adenylosuccinate and the altered balance in the enzyme's downstream metabolites (IMP, Inosine, Hypoxanthine, Xanthine) seen in the PFC of PCPtreated animals. While the influence of prolonged NMDA receptor hypofunction on the functional activity of this specific enzyme remains to be confirmed, and clearly warrants further systematic investigation, the recent finding of altered ADSS gene expression in schizophrenia [65] and the association of ADSS gene polymorphisms with schizophrenia [66] further highlights a potential role for this metabolic pathway in this disorder. In addition, a role for this metabolic pathway in cognition and schizophrenia is supported by the observation that inherited deficiency in the enzyme responsible for the breakdown of adenylosuccinate (ASL) results in mental retardation and autistic features [67, 68]. Furthermore, the ASL gene maps to chromosome 22q13.1q13.2 in humans [69] and these chromosomal loci have been repeatedly linked to schizophrenia [70–72]. The disruption of this metabolic pathway may also contribute to the reduced rate of cerebral glucose metabolism in the PFC of PCPtreated animals [3, 4] as ASL deficiency results in hypometabolism in frontal cortical structures [73]. Overall, these results suggest that the potential role of this specific series of metabolic reactions and its enzymes in cognition and schizophrenia warrants further investigation.
Conclusions
This work addresses the scenario where a pair of networks describes two different patterns of connection between a common set of nodes. We argued from first principles that the Generalized Singular Value Decomposition (equation (4)) can form the basis of a very useful computational tool. In practice, we have shown that this new computational network reordering technique was able to identify alterations in metabolic pathways in the PFC of rats treated subchronically with PCP that may contribute to the PFC dysfunction and cognitive deficits seen in these animals. Furthermore, the metabolic pathways identified as being disrupted in the PFC of PCPtreated rats trough the application of this new computation technique clearly overlap with those metabolic species known to be disrupted in schizophrenia. Applying this new algorithm in this way also identified novel pathways that may also be relevant to schizophrenia. In this way we identified alterations in glutamate metabolism and metabolic pathways central to glutamatergic neurotransmission, alterations in arginine and proline metabolism and the disruption of a novel series of purine reactions that may contribute to the PFC dysfunction and cognitive deficits seen in schizophrenia.
Methods
Chemicals
The solvents used for the study were purchased from the following sources: Acetonitrile, methanol and chloroform (Fisher Scientific, Leicestershire, UK) and formic acid (VWR, Poole, UK). All chemicals used were of analytical reagent grade. A Direct Q3^{®} water purification system (Millipore, Watford, UK) was used to produce HPLC grade water which was used in all analysis. Standards for 90 common biomolecules were also purchased which were used to characterize the ZICHILIC column (Sigma Aldrich, Dorset UK).
Animals
All experiments were completed using male Lister Hooded rats (HarlanOlac, UK) housed under standard conditions (21°C, 4565% humidity, 12h dark/light cycle (lights on 0600h) with food and drinking water available ad libitum). All manipulations were carried out at least 1 week after entry into the facility and all experiments were carried out under the Animals (Scientific Procedures) Act 1986. Animals received either subchronic treatment with vehicle (0.9% saline, i.p., n = 5) or 2.58mg.kg^{1} PCP.HCl (i.p., Sigma Aldrich, UK) once daily for five consecutive days (n = 5). At 72 hours after the final drug treatment dose animals were sacrificed and the brain rapidly dissected out and frozen in isopentane (40°C) and stored at 80°C until sectioning. Frozen brains were sectioned (20 μM) in the coronal plane in a cryostat (20°C). Tissue sections from the prefrontal cortex (PFC, Bregma +4.70mm to Bregma +3.20mm) were collected in 4ml glass vials with reference to a stereotactic rat brain atlas [74] and stored at 80°C until further preparation for LCMS analysis.
Extraction of polar metabolites from brain samples for LCMS analysis
Extraction of polar metabolites from brain tissue was carried out using the twostep extraction method described previously [75], using methanol, water and chloroform for the optimal extraction of polar metabolites. A hand held homogenizer was used to homogenize the samples once in solution. For preparation of samples for LCMS analysis 200 μl of the collected polar extract was added to 600 μl of 1 : 1 acetonitrile:water solution to produce a final solvent:sample ratio of 3 : 1. The samples were then filtered using Acrodisc 13mm syringe filters with 0.2 μm nylon membrane (Sigma Aldrich) before LCMS analysis.
LCMS analysis of polar metabolites
Experiments were carried out using a Finnigan LTQ Orbitrap (Thermo Fisher, Hemel Hempstead, UK) using 30000 resolution. Analysis was carried out in positive mode over a mass range of 601000 m/z. The capillary temperature was set at 250°C and in positive ionization mode the ion spray voltage was 4.5 kV , the capillary voltage 30 V and the tube lens voltage 105 V . The sheath and auxiliary gas flow rates were 45 and 15, respectively (units not specified by manufacturer). A ZICHILIC column (5 μm, 150 × 4.6 mm; HiChrom, Reading, UK) was used in all analysis and a binary gradient method was developed which produced good polar metabolite separation. Solvent A was 0.1% v/v formic acid in HPLC grade water and solvent B was 0.1% v/v formic acid in acetonitrile. A flow rate of 0.3 ml/min. was used and the injection volume was 10 μl. The gradient programme used was 80% B at 0 min. to 50% B at 12 min. to 20% B at 28 min. to 80% B at 37 min., with total run time of 45 minutes. The instrument was externally calibrated before analysis and internally calibrated using lock masses at m/z 83.06037 and m/z 195.08625. Samples were analysed sequentially and the vial tray temperature was set at a constant temperature of 4°C.
Data preparation and analysis
Determination of overt alterations in metabolite levels between experimental groups
The software program Xcalibur (version 2.0) was used to acquire the LCMS data. The raw Xcalibur data files from version 1.2 (Thermo Fisher, Hemel Hempstead, UK). SIEVE software (ThermoFisher Scientific) was used to identify all metabolites affected by drug treatment by calculating a pvalue and ratio based on the difference in average intensities of individual peaks, which correspond to different metabolites, between PCPtreated and control animals. A significant difference in the level of each metabolite between groups was set at p value < 0.05 and/or ratio less than 0.5 for downregulated metabolites and greater than 2 for upregulated metabolites. The ratio is the fold change in average peak intensities from control and treatment groups. For metabolite identification the masses of the polar metabolites were compared to the exact masses of 6000 biomolecules using an inhouse developed macro (Excel, Microsoft 2007).
Hypergeometric probability testing
Significant overrepresentation of a given functional group in any GSVD defined significant cluster was set by a hypergeometric probability threshold of 0.05.
Declarations
Acknowledgements
XX is supported by Engineering and Physical Sciences Research Council grant EP/E049370/1.
ND is supported by Psychiatric Research Institute in Glasgow (PsyRING), a joint initiative between the Universities of Glasgow and Strathclyde and the National Health Service of Greater Glasgow and Clyde. LM and DGW are supported by the Scottish Universities Life Sciences Alliance (SULSA).
DJH is supported by Engineering and Physical Sciences Research Council grant EP/E049370/1 and Medical Research Council grant G0601353.
Authors’ Affiliations
References
 Davidson LL, Heinrichs RW: Quantification of frontal and temporal lobe brainimaging findings in schizophrenia: a metaanalysis. Psychiatry Research: Neuroimaging. 2003, 122: 6987. 10.1016/S09254927(02)00118X.View ArticlePubMedGoogle Scholar
 Hill K, Mann L, Laws KR, Stephenson CME, NimmoSmith I, McKenna PJ: Hypofrontality in schizophrenia: a metaanalysis of functional imaging studies. Acta Psychiatrica Scandinavica. 2004, 110: 243256. 10.1111/j.16000447.2004.00376.x.View ArticlePubMedGoogle Scholar
 Cochran SM, Kennedy M, McKerchar CE, Steward LJ, Pratt JA, Morris BJ: Induction of Metabolic Hypofunction and Neurochemical Deficits after Chronic Intermittent Exposure to Phencyclidine: Differential Modulation by Antipsychotic Drugs. Neuropsychopharmacology. 2003, 28: 265275. 10.1038/sj.npp.1300031.View ArticlePubMedGoogle Scholar
 Dawson N, Thompson RJ, McVie A, Thomson DM, Morris BJ, Pratt JA: Modafinil reverses phencyclidine (PCP)induced deficits in cognitive flexibility, cerebral metabolism and functional brain connectivity. Schizophrenia Bulletin.Google Scholar
 Egerton A, Reid L, McGregor S, Cochran SM, Morris BJ, Pratt JA: Subchronic and chronic PCP treatment produces temporally distinct deficits in attentional set shifting and prepulse inhibition in rats. Psychopharmacology. 2008, 198: 3749. 10.1007/s0021300810715.View ArticlePubMedGoogle Scholar
 Egerton A, Reid L, McKerchar CE, Morris BJ, Pratt JA: Impairment in perceptual attentional setshifting following PCP administration: a rodent model of setshifting deficits in schizophrenia. Psychopharmacology. 2005, 179: 7784. 10.1007/s002130042109y.View ArticlePubMedGoogle Scholar
 Steward LJ, Kennedy MD, Morris BJ, Pratt JA: The atypical antipsychotic drug clozapine enhances chronic PCPinduced regulation of prefrontal cortex 5HT2A receptors. Neuropharmacology. 2004, 47: 527537. 10.1016/j.neuropharm.2004.04.020.View ArticlePubMedGoogle Scholar
 Fiehn O: Combining genomics, metabolome analysis, and biochemical modelling to understand metabolic networks. Comparative and Functional Genomics. 2001, 2: 155168. 10.1002/cfg.82.PubMed CentralView ArticlePubMedGoogle Scholar
 Goodacre R, Vaidyanathan S, Dunn WB, Harrigan GG, Kell DB: Metabolomics by numbers: acquiring and understanding global metabolite data. Trends in Biotechnology. 2004, 22: 245252. 10.1016/j.tibtech.2004.03.007.View ArticlePubMedGoogle Scholar
 Goodacre R, York EV, Heald JK, Scott IM: Chemometric discrimination of unfractionated plant extracts analyzed by electrospray mass spectrometry. Phytochemistry. 2003, 62: 859863. 10.1016/S00319422(02)007185.View ArticlePubMedGoogle Scholar
 Oliver SG, Winson MK, Kell DB, Baganz F: Systematic functional analysis of the yeast genome. Trends in Biotechnology. 1998, 16: 373378. 10.1016/S01677799(98)012141.View ArticlePubMedGoogle Scholar
 KaddurahDaouk R, Krishnan KRR: Metabolomics: A Global Biochemical Approach to the Study of Central Nervous System Diseases. Neuropsychopharmacology. 2009, 34: 173186. 10.1038/npp.2008.174.View ArticlePubMedGoogle Scholar
 Quinones MP, KaddurahDaouk R: Metabolomics tools for identifying biomarkers for neuropsychiatric diseases. Neurobiology of Disease. 2009, 35: 165176. 10.1016/j.nbd.2009.02.019.View ArticlePubMedGoogle Scholar
 Greenberg N, Grassano A, Thambisetty M, Lovestone S, LegidoQuigley C: A proposed metabolic strategy for monitoring disease progression in Alzheimer's disease. Electrophoresis. 2009, 30: 12351239. 10.1002/elps.200800589.View ArticlePubMedGoogle Scholar
 Chen C, Gonzalez FJ, Idle JR: LCMSbased metabolomics in drug metabolism. Drug Metab Rev. 2007, 39: 581597. 10.1080/03602530701497804.PubMed CentralView ArticlePubMedGoogle Scholar
 Lindon JC, Holmes E, Nicholson JK: Metabonomics in pharmaceutical research and development. FEBS Journal. 2007, 274: 11401151. 10.1111/j.17424658.2007.05673.x.View ArticlePubMedGoogle Scholar
 Robertson DG: Metabonomics in Toxicology: A Review. Toxicological Sciences. 2005, 85: 809822. 10.1093/toxsci/kfi102.View ArticlePubMedGoogle Scholar
 Cheng KK, Benson GM, Grimsditch DC, Reid DG, Connor SC, Griffin JL: Metabolomic study of the LDL receptor null mouse fed a highfat diet reveals profound perturbations in choline metabolism that are shared with ApoE null mice. Physiological Genomics. 2010, 41: 224231. 10.1152/physiolgenomics.00188.2009.View ArticlePubMedGoogle Scholar
 Dettmer K, Aronov PA, Hammock BD: Mass spectrometrybased metabolomics. Mass Spectrometry Reviews. 2007, 26: 5178. 10.1002/mas.20108.PubMed CentralView ArticlePubMedGoogle Scholar
 Lu W, Bennett BD, Rabinowitz JD: Analytical strategies for LCMSbased targeted metabolomics. Journal of Chromatography B. 2008, 871: 236242. 10.1016/j.jchromb.2008.04.031.View ArticleGoogle Scholar
 Wilson ID, Plumb R, Granger J, Major H, Williams R, Lenz EM: HPLCMSbased methods for the study of metabonomics. Journal of Chromatography B. 2005, 817: 6776. 10.1016/j.jchromb.2004.07.045.View ArticleGoogle Scholar
 Ohta D, Kanaya S, Suzuki H: Application of Fouriertransform ion cyclotron resonance mass spectrometry to metabolic profiling and metabolite identification. Current Opinion in Biotechnology. 2010, 21: 3544. 10.1016/j.copbio.2010.01.012.View ArticlePubMedGoogle Scholar
 Hu Q, Noll RJ, Li H, Makarov A, Hardman M, Graham Cooks R: The Orbitrap: a new mass spectrometer. Journal of Mass Spectrometry. 2005, 40: 430443. 10.1002/jms.856.View ArticlePubMedGoogle Scholar
 Calvano CD, Zambonin CG, Jensen ON: Assessment of lectin and HILIC based enrichment protocols for characterization of serum glycoproteins by mass spectrometry. Journal of Proteomics. 2008, 71: 304317. 10.1016/j.jprot.2008.06.013.View ArticlePubMedGoogle Scholar
 Cubbon S, Bradbury T, Wilson J, ThomasOates J: Hydrophilic Interaction Chromatography for Mass Spectrometric Metabonomic Studies of Urine. Analytical Chemistry. 2007, 79: 89118918. 10.1021/ac071008v.View ArticlePubMedGoogle Scholar
 Gika HG, Theodoridis GA, Wilson ID: Hydrophilic interaction and reversedphase ultraperformance liquid chromatography TOFMS for metabonomic analysis of Zucker rat urine. Journal of Separation Science. 2008, 31 (9): 15981608. 10.1002/jssc.200700644.View ArticlePubMedGoogle Scholar
 Idborg H, Zamani L, Edlund PO, SchuppeKoistinen I, Jacobsson SP: Metabolic fingerprinting of rat urine by LC/MS: Part 1. Analysis by hydrophilic interaction liquid chromatographyelectrospray ionization mass spectrometry. Journal of Chromatography B. 2005, 828: 913. 10.1016/j.jchromb.2005.07.031.View ArticleGoogle Scholar
 Kind T, Tolstikov V, Fiehn O, Weiss RH: A comprehensive urinary metabolomic approach for identifying kidney cancer. Analytical Biochemistry. 2007, 363: 185195. 10.1016/j.ab.2007.01.028.View ArticlePubMedGoogle Scholar
 Nordstrom A, Want E, Northen T, Lehtio J, Siuzdak G: Multiple Ionization Mass Spectrometry Strategy Used To Reveal the Complexity of Metabolomics. Analytical Chemistry. 2008, 80: 421429. 10.1021/ac701982e.View ArticlePubMedGoogle Scholar
 Paek IB, Moon Y, Ji HY, Kim HH, Lee HW, Lee YB, Lee HS: Hydrophilic interaction liquid chromatographytandem mass spectrometry for the determination of levosulpiride in human plasma. Journal of Chromatography B. 2004, 809: 345350.View ArticleGoogle Scholar
 Zhang X, Rauch A, Lee H, Xiao H, Rainer G, Logothetis NK: Capillary hydrophilic interaction chromatography/mass spectrometry for simultaneous determination of multiple neurotransmitters in primate cerebral cortex. Rapid Communications in Mass Spectrometry. 2007, 21: 36213628. 10.1002/rcm.3251.View ArticlePubMedGoogle Scholar
 Dunn WB, Bailey NJC, Johnson HE: Measuring the metabolome: current analytical technologies. The Analyst. 2005, 130: 606625. 10.1039/b418288j.View ArticlePubMedGoogle Scholar
 Spagou K, Tsoukali H, Raikos N, Gika H, Wilson ID, Theodoridis G: Hydrophilic interaction chromatography coupled to MS for metabonomic/metabolomic studies. Journal of Separation Science. 2010, 33: 716727. 10.1002/jssc.200900803.View ArticlePubMedGoogle Scholar
 Albert R, Barabási AL: Statistical mechanics of complex networks. Reviews of Modern Physics. 2002, 74: 4797. 10.1103/RevModPhys.74.47.View ArticleGoogle Scholar
 Newman MEJ: The structure and function of complex networks. SIAM Review. 2003, 45: 167256. 10.1137/S003614450342480.View ArticleGoogle Scholar
 Strogatz SH: Exploring complex networks. Nature. 2001, 410: 268276. 10.1038/35065725.View ArticlePubMedGoogle Scholar
 Shi J, Malik J: Normalized Cuts and Image Segmentation. IEEE Transactions on Pattern Analysis and Machine Intelligence. 2000, 22: 888905. 10.1109/34.868688.View ArticleGoogle Scholar
 Higham DJ: Spectral reordering of a rangedependent weighted random graph. IMA J. Numer. Anal. 2005, 25: 443457. 10.1093/imanum/dri003.View ArticleGoogle Scholar
 Estrada E: Protein bipartivity and essentiality in the yeast proteinprotein interaction network. J. Proteome Res. 2006, 5: 21772184. 10.1021/pr060106e.View ArticlePubMedGoogle Scholar
 Estrada E, Higham DJ, Hatano N: Communicability and multipartite structures in complex networks at negative absolute temperatures. Phys. Rev. E. 2008, 77: 026102View ArticleGoogle Scholar
 Morrison JL, Breitling R, Higham DJ, Gilbert DR: A LockandKey Model for ProteinProtein Interactions. Bioinformatics. 2006, 2: 20122019.View ArticleGoogle Scholar
 Thomas A, Cannings R, Monk NAM, Cannings C: On the structure of proteinprotein interaction networks. Biochemical Soc. Trans. 2003, 31: 14911496. 10.1042/BST0311491.View ArticleGoogle Scholar
 Alter O, Brown PO, Botstein D: Generalized Singular Value Decomposition for Comparative Analysis of GenomeScale Expression Datasets of Two Different Organisms. Proceedings of the National Academy of Sciences. 2003, 100: 33513356. 10.1073/pnas.0530258100.View ArticleGoogle Scholar
 Schreiber A, Shirley NJ, Burton RA, Fincher GB: Combining transcriptional datasets using the generalized singular value decomposition. BMC Bioinformatics. 2008, 9: 33510.1186/147121059335.PubMed CentralView ArticlePubMedGoogle Scholar
 Higham DJ, Kalna G, Kibble M: Spectral clustering and its use in bioinformatics. J. Computational and Applied Math. 2007, 204: 2537. 10.1016/j.cam.2006.04.026.View ArticleGoogle Scholar
 Fiedler M: A property of eigenvectors of nonnegative symmetric matrices and its application to graph theory. Czechoslovak Mathematical Journal. 1975, 25: 619633.View ArticleGoogle Scholar
 Atkins JE, Boman EG, Hendrickson B: A spectral algorithm for seriation and the consecutive ones problem. SIAM Journal on Computing. 1998, 28: 297310. 10.1137/S0097539795285771.View ArticleGoogle Scholar
 RoblesKelly A, Hancock ER: Graph Edit Distance from Spectral Seriation. IEEE Transactions on Pattern Analysis and Machine Intelligence. 2005, 27: 365378.View ArticlePubMedGoogle Scholar
 Golub GH, Van Loan CF: Matrix Computations. 1996, Baltimore: Johns Hopkins University PressGoogle Scholar
 Strang G: Computational Science and Engineering. 2008, WellesleyCambridge PressGoogle Scholar
 Estrada E: Topological structural classes of complex networks. Physical Review E. 2007, 75: 016103View ArticleGoogle Scholar
 Xiao X: Complex Networks and the Generalized Singular Value Decomposition. PhD thesis. 2010, University of Strathclyde, Department of Mathematics and StatisticsGoogle Scholar
 Murai R, Noda Y, Matsui K, Kamei H, Mouri A, Matsuba K, Nitta A, Furukawa H, Nabeshima T: Hypofunctional glutamatergic neurotransmission in the prefrontal cortex is involved in the emotional deficit induced by repeated treatment with phencyclidine in mice: Implications for abnormalities of glutamate release and NMDACaMKII signaling. Behavioural Brain Research. 2007, 180: 152160. 10.1016/j.bbr.2007.03.003.View ArticlePubMedGoogle Scholar
 Hahn CG, Wang HY, Cho DS, Talbot K, Gur RE, Berrettini WH, Bakshi K, Kamins J, BorgmannWinter KE, Siegel SJ, Gallop RJ, Arnold SE: Altered neuregulin 1erbB4 signaling contributes to NMDA receptor hypofunction in schizophrenia. Nature Medicine. 2006, 12: 824828. 10.1038/nm1418.View ArticlePubMedGoogle Scholar
 Lewis DA, Moghaddam B: Cognitive Dysfunction in Schizophrenia: Convergence of gammaAminobutyric Acid and Glutamate Alterations. Arch Neurol. 2006, 63: 13721376. 10.1001/archneur.63.10.1372.View ArticlePubMedGoogle Scholar
 Javitt DC, Balla A, Burch S, Suckow R, Xie S, Sershen H: Reversal of PhencyclidineInduced Dopaminergic Dysregulation by NMethylDAspartate Receptor//Glycinesite Agonists. Neuropsychopharmacology. 2003, 29: 300307.View ArticleGoogle Scholar
 Sershen H, Balla A, Aspromonte JM, Xie S, Cooper TB, Javitt DC: Characterization of interactions between phencyclidine and amphetamine in rodent prefrontal cortex and striatum: Implications in NMDA/glycinesitemediated dopaminergic dysregulation and dopamine transporter function. Neurochemistry International. 2008, 52: 119129. 10.1016/j.neuint.2007.07.011.View ArticlePubMedGoogle Scholar
 Yang CR, Svensson KA: Allosteric modulation of NMDA receptor via elevation of brain glycine and dserine: The therapeutic potentials for schizophrenia. Pharmacology & Therapeutics. 2008, 120: 317332. 10.1016/j.pharmthera.2008.08.004.View ArticleGoogle Scholar
 Magistretti PJ: Role of glutamate in neuronglia metabolic coupling. Am J Clin Nutr. 2009, 90: 875S880S. 10.3945/ajcn.2009.27462CC.View ArticlePubMedGoogle Scholar
 Patel AB, de Graaf RA, Mason GF, Rothman DL, Shulman RG, Behar KL: The contribution of GABA to glutamate/glutamine cycling and energy metabolism in the rat cortex in vivo. Proceedings of the National Academy of Sciences. 2005, 102: 55885593. 10.1073/pnas.0501703102.View ArticleGoogle Scholar
 Rameau GA, Tukey DS, GarcinHosfield ED, Titcombe RF, Misra C, Khatri L, Getzoff ED, Ziff EB: Biphasic Coupling of Neuronal Nitric Oxide Synthase Phosphorylation to the NMDA Receptor Regulates AMPA Receptor Trafficking and Neuronal Cell Death. The Journal of Neuroscience. 2007, 27: 34453455. 10.1523/JNEUROSCI.479906.2007.View ArticlePubMedGoogle Scholar
 Fejgin K, Palsson E, Wass C, Svensson L, Klamer D: Nitric Oxide Signaling in the Medial Prefrontal Cortex is Involved in the Biochemical and Behavioral Effects of Phencyclidine. Neuropsychopharmacology. 2008, 33: 18741833. 10.1038/sj.npp.1301587.View ArticlePubMedGoogle Scholar
 Wass C, Svensson L, Fejgin K, Pålsson E, Archer T, Engel JA, Klamer D: Nitric oxide synthase inhibition attenuates phencyclidineinduced disruption of cognitive flexibility. Pharmacology Biochemistry and Behavior. 2008, 89: 352359. 10.1016/j.pbb.2008.01.011.View ArticleGoogle Scholar
 Xing G, Chavko M, Zhang LX, Yang S, Post RM: Decreased calciumdependent constitutive nitric oxide synthase (cNOS) activity in prefrontal cortex in schizophrenia and depression. Schizophrenia Research. 2002, 58: 2130. 10.1016/S09209964(01)003887.View ArticlePubMedGoogle Scholar
 Tsuang MT, Nossova N, Yager T, Tsuang MM, Guo SC, Shyu KG, Glatt SJ, Liew C: Assessing the validity of bloodbased gene expression profiles for the classification of schizophrenia and bipolar disorder: A preliminary report. American Journal of Medical Genetics Part B: Neuropsychiatric Genetics. 2005, 133B: 15. 10.1002/ajmg.b.30161.View ArticleGoogle Scholar
 Zhang F, Xu Y, Liu P, Fan H, Huang X, Sun G, Song Y, Sham P: Association analyses of the interaction between the ADSS and ATM genes with schizophrenia in a Chinese population. BMC Medical Genetics. 2008, 9: 119125. 10.1186/147123509119.PubMed CentralView ArticlePubMedGoogle Scholar
 Spiegel EK, Colman RF, Patterson D: Adenylosuccinate lyase deficiency. Molecular Genetics and Metabolism. 2006, 89 (12): 1931. 10.1016/j.ymgme.2006.04.018.View ArticlePubMedGoogle Scholar
 Stone RL, Aimi J, Barshop BA, Jaeken J, Van den Berghe G, Zalkin H, Dixon JE: A mutation in adenylosuccinate lyase associated with mental retardation and autistic features. Nature Genetics. 1992, 1: 5963. 10.1038/ng049259.View ArticlePubMedGoogle Scholar
 Fon E, Demczuk S, Delattre O, Thomas G, Rouleau G: Mapping of the human adenylosuccinate lyase (ADSL) gene to chromosome 22q13.1→q13.2. Cytogenetics and Cell Genetics. 1993, 64: 201203. 10.1159/000133575.View ArticlePubMedGoogle Scholar
 Jorgensen T, Børglum A, Mors O, Wang A, Pinaud M, Flint T, Dahl H, Vang M, Kruse T, Ewald H: Search for common haplotypes on chromosome 22q in patients with schizophrenia or bipolar disorder from the Faroe Islands. American Journal of Medical Genetics. 2002, 114: 245252. 10.1002/ajmg.10191.View ArticlePubMedGoogle Scholar
 Pulver AE, Karayiorgou M, Wolyniec PS, Lasseter VK, Kasch L, Nestadt G, Antonarakis S, Housman D, Kazazian HH, Meyers D, Ott J, Lamacz M, Liang KY, Hanfelt J, Ullrich G, DeMarchi N, Ramu E, McHugh PR, Adler L, Thomas M, Carpenter WT, Manschreck T, Gordon CT, Kimberland M, Babb R, Puck J, Childs B: Sequential strategy to identify a susceptibility gene for schizophrenia: Report of potential linkage on chromosome 22q12q13.1: Part 1. American Journal of Medical Genetics. 1994, 54: 3643. 10.1002/ajmg.1320540108.View ArticlePubMedGoogle Scholar
 Saleem Q, Dash D, Gandhi C, Kishore A, Benegal V, Sherrin T, Mukherjee O, Jain S, Brahmachari SK: Association of CAG repeat loci on chromosome 22 with schizophrenia and bipolar disorder. Molecular Psychiatry. 2001, 6: 694700. 10.1038/sj.mp.4000924.View ArticlePubMedGoogle Scholar
 De Volder AG, Jaeken J, Den Berghe GV, Bol A, Michel C, Cogneau M, Goffinet AM: Regional Brain Glucose Utilization in AdenylosuccinaseDeficient Patients Measured by Positron Emission Tomography. Pediatric Research. 1998, 24: 238242.View ArticleGoogle Scholar
 Paxinos G, Watson C: The Rat Brain in Stereotaxic Coordinates. 1998, Academic PressGoogle Scholar
 Wu H, Southam AD, Hines A, Viant MR: Highthroughput tissue extraction protocol for NMR and MSbased metabolomics. Analytical Biochemistry. 2008, 372: 204212. 10.1016/j.ab.2007.10.002.View ArticlePubMedGoogle Scholar
Copyright
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.