- Research article
- Open Access
The application of selective reaction monitoring confirms dysregulation of glycolysis in a preclinical model of schizophrenia
- Daniel Martins-de-Souza†1Email author,
- Murtada Alsaif†1,
- Agnes Ernst†1,
- Laura W Harris1,
- Nancy Aerts2,
- Ilse Lenaerts2,
- Pieter J Peeters2,
- Bob Amess1,
- Hassan Rahmoune1,
- Sabine Bahn1, 3 and
- Paul C Guest1
© Martins-de-Souza et al; licensee BioMed Central Ltd. 2012
Received: 17 December 2011
Accepted: 15 March 2012
Published: 15 March 2012
Establishing preclinical models is essential for novel drug discovery in schizophrenia. Most existing models are characterized by abnormalities in behavioral readouts, which are informative, but do not necessarily translate to the symptoms of the human disease. Therefore, there is a necessity of characterizing the preclinical models from a molecular point of view. Selective reaction monitoring (SRM) has already shown promise in preclinical and clinical studies for multiplex measurement of diagnostic, prognostic and treatment-related biomarkers.
We have established an SRM assay for multiplex analysis of 7 enzymes of the glycolysis pathway which is already known to be affected in human schizophrenia and in the widely-used acute PCP rat model of schizophrenia. The selected enzymes were hexokinase 1 (Hk1), aldolase C (Aldoc), triosephosphate isomerase (Tpi1), glyceraldehyde-3-phosphate dehydrogenase (Gapdh), phosphoglycerate mutase 1 (Pgam1), phosphoglycerate kinase 1 (Pgk1) and enolase 2 (Eno2). The levels of these enzymes were analyzed using SRM in frontal cortex from brain tissue of PCP treated rats.
Univariate analyses showed statistically significant altered levels of Tpi1 and alteration of Hk1, Aldoc, Pgam1 and Gapdh with borderline significance in PCP rats compared to controls. Most interestingly, multivariate analysis which considered the levels of all 7 enzymes simultaneously resulted in generation of a bi-dimensional chart that can distinguish the PCP rats from the controls.
This study not only supports PCP treated rats as a useful preclinical model of schizophrenia, but it also establishes that SRM mass spectrometry could be used in the development of multiplex classification tools for complex psychiatric disorders such as schizophrenia.
Schizophrenia is a psychiatric disorder of uncertain etiology that affects about 1% of the world's population. It is likely to be caused by a complex crosstalk between genetic, neurodevelopmental and environmental factors  resulting in the differential expression of several genes and proteins. The multi-factorial characteristic of schizophrenia not only hinders its biochemical understanding and characterization, but makes it difficult to establish preclinical models, which are needed urgently for testing of existing and novel therapies . The current evaluation of preclinical models for psychiatric disorders is mostly based on readouts following behavioral tests which simulate particular features or symptoms of the disease. However, characteristic schizophrenia symptoms, such as hallucinations, delusions and disorganized thoughts cannot be assessed using these tests. In contrast, characterization of these models using molecular fingerprinting approaches can be translated directly between preclinical and clinical studies.
The blockade of NMDA receptors using drugs such as phencyclidine (PCP) and ketamine is widely used to mimic schizophrenia in rodents. The psychosis-inducing effects of these drugs in humans and the widespread nature of the NMDA receptor in the brain result in multiple effects which may resemble the multifactorial features of schizophrenia . Previous studies have shown that PCP-treated rats present altered glucose utilization in the limbic regions of their brains [4, 5], consistent with findings from post mortem brain studies of schizophrenia patients .
Most targeted protein expression analyses are performed using either Western blotting (WB) or enzyme-linked immunosorbent assay (ELISA) approaches. Although these methods are effective, both are dependent on antibody availability and neither are suitable for large-scale validation in their standard formats. Moreover, depending on the quality of the antibody employed, WB analysis is only considered to be a semi-quantitative method . An alternative for targeted protein expression analysis is selective reaction monitoring (SRM) mass spectrometry [or multiple reaction monitoring (MRM)]. SRM is a mass spectrometry-based technique which can accurately measure the relative or absolute concentrations of selected peptides/proteins down to the attomolar range. In addition, SRM can be multiplexed, allowing the measurement of peptides from several proteins in a single experiment. Therefore, this method has the capability of measuring all components in a metabolic pathway in a single sample.
In this study, we have performed SRM mass spectrometry analysis of the seven glycolysis enzymes listed above in brain tissues of PCP-treated (1.0 and 2.5 mg/kg) and control rats. The resulting protein levels were then assessed using univariate and multivariate analyses, which showed that the acute PCP model reproduces some of the glycolytic deficits seen in schizophrenia at the molecular level.
Results and discussion
As a multifactorial disorder, schizophrenia is likely to be a consequence of several molecular differences of small effect  triggered by environmental factors . These subtle and sometimes even non-significant molecular differences may be one of the reasons why it has been difficult to replicate and even determine particular genes or biomarkers responsible for schizophrenia . However, the combination of these small molecular differences may lead to major biological disturbances, especially if the differentially expressed molecules belong to the same biochemical pathway. We were therefore interested in determining whether the combination of the differential expression of glycolytic enzymes can be used to differentiate between controls and PCP treated groups despite the lack of statistical significance at the univariate level. In addition, we aim to check whether SRM is a suitable proteomic tool for this task.
The translation of academic findings to the clinic is now a major objective of biomarker validation studies, especially in support of drug discovery. The Food and Drug Administration (FDA) has now called for efforts to modernize methods, tools and techniques for the purpose of delivering more efficacious and safer drugs [25, 26]. SRM mass spectrometry can be optimized for performing high throughput multiplex analyses, which is robust and user friendly for use in the clinical environment, as stated above. Moreover, we have used the multiplex SRM mass spectrometry approach here to determine molecular phenotypes of a pre-clinical model of schizophrenia. This demonstrated the potential of translating this model for use in drug discovery with a set of companion biomarkers representing the glycolytic pathway to facilitate this process.
The validation of large scale proteomic profiling results has been a bottle neck in biomarker discovery. This is due in part to the lack of a robust technological pipeline connecting discovery with validation and translation, which is applicable to a clinical setting . Here, we have shown that SRM mass spectrometry is a useful platform for achieving this goal. This approach has already been used in clinical studies such as those screening for the levels of dihydroartemisinin (DHA) in plasma of malaria patients , for determining the levels of apolipoproteins in human plasma  and for cancer biomarker screening . Finally, although in a small number of samples, the results of the analysis support a role for the PCP rat as a good model for some aspects of schizophrenia, including metabolic dysfunction, as this has shown dysregulation of the glycolysis-related enzymes as previously shown in brain tissue from schizophrenia patients.
Limitations of our study include the number of analyzed samples. Replication for validating the biological interpretation here is suggested, which can also confirm the validity of our assessment using PLS-DA. Methodologically, Tpi1 was only analyzed by one peptide. All other possible Tpi1 peptides were neither proteotypic  nor quantotypic  as well as not prone to post-translational modifications as verified by PHOSIDA http://www.phosida.com.
Glycolytic enzymes candidates and peptides analyzed
Our findings warrant further biological research on the analysis of glycolytic enzymes in other schizophrenia models, including investigating the effects of drug treatments. Further technical developments should be directed towards the use of SRM as a multiplex tool for diagnostic purposes. In addition, pathways representing other aspects of schizophrenia could be developed using the SRM format, especially in cases of models where other pathways are known to be modified.
Male, adult Lister Hooded rats (Harlan, Horst, The Netherlands), weighing 180-220 g were housed in groups of three to four under standard laboratory conditions at a temperature of 21°C (± 1°C). Rats were maintained on a 12-h/12-h light/dark cycle (lights on at 08:00). Experiments were conducted during the light cycle. Food and water were available ad libitum. All experiments were performed in accordance with the European Communities Council Directives (86/609/EEC) and were approved by the policies of the ethical committee of Johnson and Johnson Pharma (Beerse, Belgium). Rats received an acute subcutaneous injection of either 1 mg/kg (n = 9) or 2.5 mg/kg (n = 7) PCP as these dosages have been shown previously to significantly enhance locomotion in open field testing . PCP (Johnson and Johnson Pharma, Beerse, Belgium) was dissolved in 0.9% saline. The control group (n = 5) received an equivalent volume of saline (Baxter, Belgium). Sixty minutes after PCP or saline administration, rats were decapitated and brains were dissected immediately on ice to obtain frontal cortex (FC) samples, which were stored at -80°C until further use.
Brain tissue samples were sonicated for 2 minutes in a buffer containing 7 M urea, 2 M thiourea, 4% CHAPS, 2% ASB14, 70 mM dithiothreitol (DTT) (all from Sigma) at a 5:1 (v:w) ratio  for proteome extraction. Samples were vortexed for 30 minutes and centrifuged for 3 minutes at 17000 g at 4°C. Protein concentrations in the supernatants were determined using a Bradford assay (Biorad, UK). Extracted proteins (100 mg) were precipitated overnight using cold acetone (-20°C) and centrifuged for 10 min at 17000 g at 4°C. The supernatants were discarded and the pellets dissolved in 50 mM ammonium bicarbonate. Next, sulfhydryl groups were reduced with 5 mM DTT at 60°C for 30 minutes followed by alkylation with 10 mM iodacetamide in the dark at 37°C for 30 minutes. Each resultant proteome was digested using trypsin at a 1:50 (w:w) ratio for 17 hours at 37°C. Reactions were stopped by addition of 8.8 M HCl at a 1:60 (w:w) ratio. The concentrations of the digested proteins (also referred to as the peptidome) were adjusted to 0.1 μg/μL by addition of the appropriate amounts of 0.1% formic acid. Samples were spiked with 25 fmol of yeast enolase for data normalization.
Mass spectrometry experiments
Peptidomes were injected in duplicate into a nanoUltra performance liquid chromatography instrument containing a BEH-130 C18 column (75 um × 200 mm) at a flow rate of 0.3 μL/min, coupled online to a Xevo-triple-quadrupole mass spectrometer (Waters; Milford, MA, USA). The separation buffers were A) 0.1% formic acid and B) acetonitrile in 0.1% formic acid. For separation of peptides, the following 48 min gradient was applied: 97/3% (A/B) to 60/40%B in 30 min; 60/40% to 15/85% in 2 min; 5 min at 15/85%; returning to the initial condition in 1 min. Selected peptides (Table 1) eluting from the nanoLC were measured in SRM mode using an electrospray voltage of 22 kV and a cone voltage of 30 V. All SRM functions had a 2 min window of the predicted retention time and scan times were 20 milliseconds. At least 2 SRM transitions per peptide were measured. The collision energy for each transition was optimized using Skyline software (MacCoss Lab Software; Seattle, WA, USA)  based on the equation: CE = 0.034*m/z +3.314. Acquired data were processed and quantified using TargetLynx (Waters). Mass spectral intensities of peptides from the target proteins were averaged and normalized against yeast enolase intensities.
To account for potential non-Gaussian distributions, data were assessed using the Kruskal-Wallis test to find if any differences were present among the three groups. Significant differences between the two dosages of PCP and control groups were calculated using corrected Mann Whitney U p-values, applying the Benjamini-Hochberg method for multiple testing. All tests were two-tailed and conducted in Prism 5 (Graphpad Software Inc; La Jolla, CA, USA) or R 2.13 (The R Foundation for Statistical Computing, Vienna, Austria).
For the multivariate analyses, partial least squares discriminant analysis (PLS-DA) was used (Simca-P + 12; Umetrics; Umea, Sweden). PLS-DA is a supervised method for fine tuning a principle component analysis (PCA) model to maximize separation between two groups of observations. The protein intensities were uploaded to the software which assembled these linearly to form PLS components that were plotted in multidimensional space. These plots were then reduced to the principle components which are presented as a two-dimensional PLS-DA scores plot . In addition to this 2D plot, the software calculates the variable influence on projection (VIP) which assesses the contribution of each variable on the separation observed on the PLS-DA scores plot . Lastly, the loading plot (S-plot) is presented. The x-axis represents the contribution of each variable to the separation). Variables positioned furthest from the centre have the strongest contribution. The data plotted along the y-axis show the correlation of each variable from sample to sample (the further the variable is from the centre, the stronger the correlation ).
We thank the Stanley Medical Research Institute for all support. DMS specially thanks Dr. Simon Cubbon and Dr. Jonathan Fox for the continuous support with TQ-Xevo. MA was supported by a scholarship from the Ministry of Higher Education, Saudi Arabia. BA was funded by Autism Speaks. This research was supported by the IMI NEWMEDS European Framework Programme (FP7).
- Freedman R: Schizophrenia. N Engl J Med. 2003, 349 (18): 1738-1749. 10.1056/NEJMra035458.PubMedView ArticleGoogle Scholar
- Young JW, Zhou X, Geyer MA: Animal models of schizophrenia. Curr Top Behav Neurosci. 2010, 4: 391-433. 10.1007/7854_2010_62.PubMedView ArticleGoogle Scholar
- Kantrowitz JT, Javitt DC: Thinking glutamatergically: changing concepts of schizophrenia based upon changing neurochemical models. Clin Schizophr Relat Psychoses. 2010, 4 (3): 189-200. 10.3371/CSRP.4.3.6.PubMedView ArticleGoogle Scholar
- Gao XM, Shirakawa O, Du F, Tamminga CA: Delayed regional metabolic actions of phencyclidine. Eur J Pharmacol. 1993, 241 (1): 7-15. 10.1016/0014-2999(93)90926-9.PubMedView ArticleGoogle Scholar
- Iltis I, Koski DM, Eberly LE, Nelson CD, Deelchand DK, Valette J, Ugurbil K, Lim KO, Henry PG: Neurochemical changes in the rat prefrontal cortex following acute phencyclidine treatment: an in vivo localized (1)H MRS study. NMR Biomed. 2009, 22 (7): 737-744. 10.1002/nbm.1385.PubMedPubMed CentralView ArticleGoogle Scholar
- Vyas NS, Patel NH, Nijran KS, Al-Nahhas A, Puri BK: The use of PET imaging in studying cognition, genetics and pharmacotherapeutic interventions in schizophrenia. Expert Rev Neurother. 2011, 11 (1): 37-51. 10.1586/ern.10.160.PubMedView ArticleGoogle Scholar
- Horvath S, Janka Z, Mirnics K: Analyzing schizophrenia by DNA microarrays. Biol Psychiatry. 2011, 69 (2): 157-162. 10.1016/j.biopsych.2010.07.017.PubMedPubMed CentralView ArticleGoogle Scholar
- Martins De-Souza D, Dias-Neto E: RNA Biomarkers in Schizophrenia. Biomarkers for Psychiatric Disorders. 2009, Turck CW: Springer, 97-127. 1Google Scholar
- Martins-De-Souza D, Dias-Neto E, Schmitt A, Falkai P, Gormanns P, Maccarrone G, Turck CW, Gattaz WF: Proteome analysis of schizophrenia brain tissue. World J Biol Psychiatry. 2010, 11 (2): 110-120. 10.3109/15622970903490626.PubMedView ArticleGoogle Scholar
- English JA, Pennington K, Dunn MJ, Cotter DR: The neuroproteomics of schizophrenia. Biol Psychiatry. 2011, 69 (2): 163-172. 10.1016/j.biopsych.2010.06.031.PubMedView ArticleGoogle Scholar
- Martins-de-Souza D, Harris LW, Guest PC, Bahn S: The role of energy metabolism dysfunction and oxidative stress in schizophrenia revealed by proteomics. Antioxid Redox Signal. 2011, 15 (7): 2067-2079. 10.1089/ars.2010.3459.PubMedView ArticleGoogle Scholar
- Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28 (1): 27-30. 10.1093/nar/28.1.27.PubMedPubMed CentralView ArticleGoogle Scholar
- Andreasen NC, Arndt S, Swayze V, Cizadlo T, Flaum M, O'Leary D, Ehrhardt JC, Yuh WT: Thalamic abnormalities in schizophrenia visualized through magnetic resonance image averaging. Science. 1994, 266 (5183): 294-298. 10.1126/science.7939669.PubMedView ArticleGoogle Scholar
- Fucetola R, Newcomer JW, Craft S, Melson AK: Age- and dose-dependent glucose-induced increases in memory and attention in schizophrenia. Psychiatry Res. 1999, 88 (1): 1-13. 10.1016/S0165-1781(99)00063-3.PubMedView ArticleGoogle Scholar
- Tamminga CA, Thaker GK, Buchanan R, Kirkpatrick B, Alphs LD, Chase TN, Carpenter WT: Limbic system abnormalities identified in schizophrenia using positron emission tomography with fluorodeoxyglucose and neocortical alterations with deficit syndrome. Arch Gen Psychiatry. 1992, 49 (7): 522-530. 10.1001/archpsyc.1992.01820070016003.PubMedView ArticleGoogle Scholar
- Rezin GT, Amboni G, Zugno AI, Quevedo J, Streck EL: Mitochondrial dysfunction and psychiatric disorders. Neurochem Res. 2009, 34 (6): 1021-1029. 10.1007/s11064-008-9865-8.PubMedView ArticleGoogle Scholar
- Martins-de-Souza D, Guest PC, Harris LW, Vanattou-Saifoudine N, Webster MJ, Rahmoune H, Bahn S: Identification of proteomic signatures associated with depression and psychotic depression in post-mortem brains from major depression patients. Transl Psychiatry. 2012, 2: e87-10.1038/tp.2012.13.PubMedPubMed CentralView ArticleGoogle Scholar
- Mann M: Can proteomics retire the western blot?. J Proteome Res. 2008, 7 (8): 3065-10.1021/pr800463v.PubMedView ArticleGoogle Scholar
- Prabakaran S, Swatton JE, Ryan MM, Huffaker SJ, Huang JT, Griffin JL, Wayland M, Freeman T, Dudbridge F, Lilley KS: Mitochondrial dysfunction in schizophrenia: evidence for compromised brain metabolism and oxidative stress. Mol Psychiatry. 2004, 9 (7): 684-697. 643.PubMedView ArticleGoogle Scholar
- Martins-de-Souza D, Gattaz WF, Schmitt A, Maccarrone G, Hunyadi-Gulyas E, Eberlin MN, Souza GH, Marangoni S, Novello JC, Turck CW: Proteomic analysis of dorsolateral prefrontal cortex indicates the involvement of cytoskeleton, oligodendrocyte, energy metabolism and new potential markers in schizophrenia. J Psychiatr Res. 2009, 43 (11): 978-986. 10.1016/j.jpsychires.2008.11.006.PubMedView ArticleGoogle Scholar
- Johnston-Wilson NL, Sims CD, Hofmann JP, Anderson L, Shore AD, Torrey EF, Yolken RH: Disease-specific alterations in frontal cortex brain proteins in schizophrenia, bipolar disorder, and major depressive disorder. The Stanley Neuropathology Consortium. Mol Psychiatry. 2000, 5 (2): 142-149. 10.1038/sj.mp.4000696.PubMedView ArticleGoogle Scholar
- Cloninger CR: The discovery of susceptibility genes for mental disorders. Proc Natl Acad Sci USA. 2002, 99 (21): 13365-13367. 10.1073/pnas.222532599.PubMedPubMed CentralView ArticleGoogle Scholar
- van Os J, Rutten BP, Poulton R: Gene-environment interactions in schizophrenia: review of epidemiological findings and future directions. Schizophr Bull. 2008, 34 (6): 1066-1082. 10.1093/schbul/sbn117.PubMedPubMed CentralView ArticleGoogle Scholar
- Collier DA: Schizophrenia: the polygene princess and the pea. Psychol Med. 2008, 38 (12): 1687-1691. 10.1017/S0033291708003668. discussion 1818-1620PubMedView ArticleGoogle Scholar
- Marson B: ''Critical Path'' is on the road forward; FDA reports industry activity is high. The Pink Sheet. 2007, 69: 29-Google Scholar
- Owens J: Funding for accelerating drug development initiative critical. Nat Rev Drug Discov. 2006, 5 (4): 271-10.1038/nrd2015.PubMedView ArticleGoogle Scholar
- Plymoth A, Hainaut P: Proteomics beyond proteomics: toward clinical applications. Curr Opin Oncol. 2011, 23 (1): 77-82. 10.1097/CCO.0b013e32834179c1.PubMedView ArticleGoogle Scholar
- Wiesner L, Govender K, Meredith SA, Norman J, Smith PJ: A liquid-liquid LC/MS/MS assay for the determination of artemether and DHA in malaria patient samples. J Pharm Biomed Anal. 2011, 55 (2): 373-378. 10.1016/j.jpba.2011.01.036.PubMedView ArticleGoogle Scholar
- Agger SA, Marney LC, Hoofnagle AN: Simultaneous quantification of apolipoprotein A-I and apolipoprotein B by liquid-chromatography-multiple- reaction-monitoring mass spectrometry. Clin Chem. 2010, 56 (12): 1804-1813. 10.1373/clinchem.2010.152264.PubMedPubMed CentralView ArticleGoogle Scholar
- Ang CS, Nice EC: Targeted in-gel MRM: a hypothesis driven approach for colorectal cancer biomarker discovery in human feces. J Proteome Res. 2010, 9 (9): 4346-4355. 10.1021/pr100509e.PubMedView ArticleGoogle Scholar
- Kuster B, Schirle M, Mallick P, Aebersold R: Scoring proteomes with proteotypic peptide probes. Nat Rev Mol Cell Biol. 2005, 6 (7): 577-583. 10.1038/nrm1683.PubMedView ArticleGoogle Scholar
- Brownridge P, Holman SW, Gaskell SJ, Grant CM, Harman VM, Hubbard SJ, Lanthaler K, Lawless C, O'Cualain R, Sims P: Global absolute quantification of a proteome: Challenges in the deployment of a QconCAT strategy. Proteomics. 2011, 11 (15): 2957-2970. 10.1002/pmic.201100039.PubMedView ArticleGoogle Scholar
- Gnad F, Gunawardena J, Mann M: PHOSIDA 2011: the posttranslational modification database. Nucleic Acids Res. 2011, 39: 253-260. 10.1093/nar/gkq1159.View ArticleGoogle Scholar
- Amann LC, Gandal MJ, Halene TB, Ehrlichman RS, White SL, McCarren HS, Siegel SJ: Mouse behavioral endophenotypes for schizophrenia. Brain Res Bull. 2010, 83 (3-4): 147-161. 10.1016/j.brainresbull.2010.04.008.PubMedView ArticleGoogle Scholar
- Neill JC, Barnes S, Cook S, Grayson B, Idris NF, McLean SL, Snigdha S, Rajagopal L, Harte MK: Animal models of cognitive dysfunction and negative symptoms of schizophrenia: focus on NMDA receptor antagonism. Pharmacol Ther. 2010, 128 (3): 419-432. 10.1016/j.pharmthera.2010.07.004.PubMedView ArticleGoogle Scholar
- Kalinichev M, Robbins MJ, Hartfield EM, Maycox PR, Moore SH, Savage KM, Austin NE, Jones DN: Comparison between intraperitoneal and subcutaneous phencyclidine administration in Sprague-Dawley rats: a locomotor activity and gene induction study. Prog Neuropsychopharmacol Biol Psychiatry. 2008, 32 (2): 414-422. 10.1016/j.pnpbp.2007.09.008.PubMedView ArticleGoogle Scholar
- Martins-de-Souza DE, Menezes-de-Oliveira B, Dos-Santos Farias A, Horiuchi RS, Crepaldi Domingues C, De-Paula E, Marangoni S, Gattaz WF, Dias-Neto E, Camillo Novello J: The use of ASB-14 in combination with CHAPS is the best for solubilization of human brain proteins for two-dimensional gel electrophoresis. Brief Funct Genomic Proteomic. 2007, 6 (1): 70-75. 10.1093/bfgp/elm009.PubMedView ArticleGoogle Scholar
- MacLean B, Tomazela DM, Shulman N, Chambers M, Finney GL, Frewen B, Kern R, Tabb DL, Liebler DC, MacCoss MJ: Skyline: an open source document editor for creating and analyzing targeted proteomics experiments. Bioinformatics. 2010, 26 (7): 966-968. 10.1093/bioinformatics/btq054.PubMedPubMed CentralView ArticleGoogle Scholar
- Schwarz E, Guest PC, Rahmoune H, Harris LW, Wang L, Leweke FM, Rothermundt M, Bogerts B, Koethe D, Kranaster L: Identification of a biological signature for schizophrenia in serum. Mol Psychiatry. 2011, doi:10.1038/mp.2011.42,Google Scholar
- Rossinia KVS, Carioua V, Qannaria EM, Fogliatto FS: PLS discriminant analysis applied to conventional sensory profiling data. Food Quality and Preference. 2012, 23: 18-24. 10.1016/j.foodqual.2011.01.005.View ArticleGoogle Scholar
- Using LC-oa-TOF MSE with a Multivariate Statistical Sample Statistical Sample Profiling Strategy to Distinguish Chinese Red Ginseng from Korean Red Ginseng. [http://chromatographyonline.findanalytichem.com/lcgc/article/articleDetail.jsp?id=632624&sk=&date=&pageID=5]
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.