Protein Co-Expression Analysis as a Strategy to Complement a Standard Quantitative Proteomics Approach: Case of a Glioblastoma Multiforme Study

Although correlation network studies from co-expression analysis are increasingly popular, they are rarely applied to proteomics datasets. Protein co-expression analysis provides a complementary view of underlying trends, which can be overlooked by conventional data analysis. The core of the present study is based on Weighted Gene Co-expression Network Analysis applied to a glioblastoma multiforme proteomic dataset. Using this method, we have identified three main modules which are associated with three different membrane associated groups; mitochondrial, endoplasmic reticulum, and a vesicle fraction. The three networks based on protein co-expression were assessed against a publicly available database (STRING) and show a statistically significant overlap. Each of the three main modules were de-clustered into smaller networks using different strategies based on the identification of highly connected networks, hierarchical clustering and enrichment of Gene Ontology functional terms. Most of the highly connected proteins found in the endoplasmic reticulum module were associated with redox activity while a core of the unfolded protein response was identified in addition to proteins involved in oxidative stress pathways. The proteins composing the electron transfer chain were found differently affected with proteins from mitochondrial Complex I being more down-regulated than proteins from Complex III. Finally, the two pyruvate kinases isoforms show major differences in their co-expressed protein networks suggesting roles in different cellular locations.

[1]  A. Krainer,et al.  The alternative splicing repressors hnRNP A1/A2 and PTB influence pyruvate kinase isoform expression and cell metabolism , 2010, Proceedings of the National Academy of Sciences.

[2]  S. Groshen,et al.  Critical role of the stress chaperone GRP78/BiP in tumor proliferation, survival, and tumor angiogenesis in transgene-induced mammary tumor development. , 2008, Cancer research.

[3]  Israel Steinfeld,et al.  BMC Bioinformatics BioMed Central , 2008 .

[4]  Michael Unser,et al.  A chemostat array enables the spatio-temporal analysis of the yeast proteome , 2013, Proceedings of the National Academy of Sciences.

[5]  Bruce J. Aronow,et al.  ToppCluster: a multiple gene list feature analyzer for comparative enrichment clustering and network-based dissection of biological systems , 2010, Nucleic Acids Res..

[6]  T. Onda,et al.  Ubiquitous mitochondrial creatine kinase downregulated in oral squamous cell carcinoma , 2006, British Journal of Cancer.

[7]  Y. Levin The role of statistical power analysis in quantitative proteomics , 2011, Proteomics.

[8]  Davide Heller,et al.  STRING v10: protein–protein interaction networks, integrated over the tree of life , 2014, Nucleic Acids Res..

[9]  J. Searcy,et al.  Proteomic analysis of mitochondria in APOE transgenic mice and in response to an ischemic challenge , 2011, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[10]  A. Herrmann,et al.  Adaptive changes in the neuronal proteome: mitochondrial energy production, endoplasmic reticulum stress, and ribosomal dysfunction in the cellular response to metabolic stress , 2013, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[11]  B. Snel,et al.  STRING: a web-server to retrieve and display the repeatedly occurring neighbourhood of a gene. , 2000, Nucleic acids research.

[12]  P. Shannon,et al.  Cytoscape: a software environment for integrated models of biomolecular interaction networks. , 2003, Genome research.

[13]  Jing Li,et al.  dbDEPC 2.0: updated database of differentially expressed proteins in human cancers , 2011, Nucleic Acids Res..

[14]  T. Le Bihan,et al.  Interactions among mitochondrial proteins altered in glioblastoma , 2014, Journal of Neuro-Oncology.

[15]  A. Tischler,et al.  LANDSCAPE OF THE MITOCHONDRIAL Hsp90 METABOLOME IN TUMORS , 2013, Nature Communications.

[16]  Y. Park,et al.  Induction of glucose-regulated protein 78 by chronic hypoxia in human gastric tumor cells through a protein kinase C-epsilon/ERK/AP-1 signaling cascade. , 2001, Cancer research.

[17]  John D. Storey,et al.  Multiple Locus Linkage Analysis of Genomewide Expression in Yeast , 2005, PLoS biology.

[18]  Y. Leea,et al.  Analysis of oncogenic signaling networks in glioblastoma identifies ASPM as a molecular target , 2006 .

[19]  Ming Yi,et al.  bioDBnet: the biological database network , 2009, Bioinform..

[20]  R. Caprioli,et al.  Proteome analysis of human colon cancer by two‐dimensional difference gel electrophoresis and mass spectrometry , 2004, Proteomics.

[21]  H. Tsuda,et al.  Overexpression of YWHAZ relates to tumor cell proliferation and malignant outcome of gastric carcinoma , 2013, British Journal of Cancer.

[22]  Eugenia G. Giannopoulou,et al.  Search for potential markers for prostate cancer diagnosis, prognosis and treatment in clinical tissue specimens using amine-specific isobaric tagging (iTRAQ) with two-dimensional liquid chromatography and tandem mass spectrometry. , 2008, Journal of proteome research.

[23]  Ru Wei,et al.  The M2 splice isoform of pyruvate kinase is important for cancer metabolism and tumour growth , 2008, Nature.

[24]  Xiaohang Zhao,et al.  Quantitative proteomic signature of liver cancer cells: tissue transglutaminase 2 could be a novel protein candidate of human hepatocellular carcinoma. , 2008, Journal of proteome research.

[25]  W. Lukiw,et al.  An Exploratory Analysis of Conservation of Co-Expressed Genes across Alzheimer's disease Progression , 2013 .

[26]  L. Alldridge,et al.  Proteome profiling of breast tumors by gel electrophoresis and nanoscale electrospray ionization mass spectrometry. , 2008, Journal of proteome research.

[27]  Ziyan Wang,et al.  Repairing DNA damage by XRCC6/KU70 reverses TLR4-deficiency-worsened HCC development via restoring senescence and autophagic flux , 2013, Autophagy.

[28]  T. Simmen,et al.  Where the endoplasmic reticulum and the mitochondrion tie the knot: the mitochondria-associated membrane (MAM). , 2013, Biochimica et biophysica acta.

[29]  Yingyi Wang,et al.  NF-κB/RelA-PKM2 mediates inhibition of glycolysis by fenofibrate in glioblastoma cells , 2015, Oncotarget.

[30]  Richard D. Smith,et al.  Quantitative proteome analysis of breast cancer cell lines using 18O‐labeling and an accurate mass and time tag strategy , 2006, Proteomics.

[31]  J. Bähler,et al.  Meta-analysis of genome regulation and expression variability across hundreds of environmental and genetic perturbations in fission yeast. , 2010, Molecular bioSystems.

[32]  Hongyu Zhao,et al.  Statistical Analysis Reveals Co-Expression Patterns of Many Pairs of Genes in Yeast Are Jointly Regulated by Interacting Loci , 2013, PLoS genetics.

[33]  Yongqiang Chen,et al.  Mitochondrial electron-transport-chain inhibitors of complexes I and II induce autophagic cell death mediated by reactive oxygen species , 2007, Journal of Cell Science.

[34]  Albert-László Barabási,et al.  Scale-free networks , 2008, Scholarpedia.

[35]  Sylvain Brohée Using the NeAT toolbox to compare networks to networks, clusters to clusters, and network to clusters. , 2012, Methods in molecular biology.

[36]  T. Ueda,et al.  Synaptic Vesicle-bound Pyruvate Kinase can Support Vesicular Glutamate Uptake , 2009, Neurochemical Research.

[37]  K. Krause,et al.  Proteomic profiling of cold thyroid nodules. , 2007, Endocrinology.

[38]  Steve Horvath,et al.  WGCNA: an R package for weighted correlation network analysis , 2008, BMC Bioinformatics.

[39]  Gary D. Bader,et al.  An automated method for finding molecular complexes in large protein interaction networks , 2003, BMC Bioinformatics.

[40]  H. Lehrach,et al.  No evidence for a shift in pyruvate kinase PKM1 to PKM2 expression during tumorigenesis , 2011, Oncotarget.

[41]  W. Hancock,et al.  Proteomic analysis of high-grade dysplastic cervical cells obtained from ThinPrep slides using laser capture microdissection and mass spectrometry. , 2007, Journal of proteome research.

[42]  Richard D. Smith,et al.  Protein co-expression network analysis (ProCoNA) , 2013, Journal of Clinical Bioinformatics.

[43]  Hsuan-Cheng Huang,et al.  Quantitative proteomic and genomic profiling reveals metastasis-related protein expression patterns in gastric cancer cells. , 2006, Journal of proteome research.

[44]  Thanura R. Elvitigala,et al.  Global Proteomics Reveal an Atypical Strategy for Carbon/Nitrogen Assimilation by a Cyanobacterium Under Diverse Environmental Perturbations* , 2010, Molecular & Cellular Proteomics.

[45]  G. Kroemer,et al.  The protein disulfide isomerases PDIA4 and PDIA6 mediate resistance to cisplatin-induced cell death in lung adenocarcinoma , 2014, Cell Death and Differentiation.