Construction and validation of a gene co-expression network in grapevine (Vitis vinifera. L.)

Abstract

Gene co-expression analysis has been widely used for predicting gene functions because genes within modules of a co-expression network may be involved in similar biological processes and exhibit similar biological functions. To detect gene relationships in the grapevine genome, we constructed a grapevine gene co-expression network (GGCN) by compiling a total of 374 publically available grapevine microarray datasets. The GGCN consisted of 557 modules containing a total of 3834 nodes with 13 479 edges. The functions of the subnetwork modules were inferred by Gene ontology (GO) enrichment analysis. In 127 of the 557 modules containing two or more GO terms, 38 modules exhibited the most significantly enriched GO terms, including ‘protein catabolism process’, ‘photosynthesis’, ‘cell biosynthesis process’, ‘biosynthesis of plant cell wall’, ‘stress response’ and other important biological processes. The ‘response to heat’ GO term was highly represented in module 17, which is composed of many heat shock proteins. To further determine the potential functions of genes in module 17, we performed a Pearson correlation coefficient test, analyzed orthologous relationships with Arabidopsis genes and established gene expression correlations with real-time quantitative reverse transcriptase PCR (qRT-PCR). Our results indicated that many genes in module 17 were upregulated during the heat shock and recovery processes and downregulated in response to low temperature. Furthermore, two putative genes, Vit_07s0185g00040 and Vit_02s0025g04060, were highly expressed in response to heat shock and recovery. This study provides insight into GGCN gene modules and offers important references for gene functions and the discovery of new genes at the module level.

References

[1] Ruan J, Dean AK, Zhang W . A general co-expression network-based approach to gene expression analysis: comparison and applications. BMC Syst Biol 2010; 4: 8.

[2] Brazhnik P, Fuente A, Mendes P . Gene networks: how to put the function in genomics. Trends Biotechnol 2002; 20: 467-472.

[3] Aoki K, Ogata Y, Shibata D . Approaches for extracting practical information from gene co-expression networks in plant biology. Plant Cell Physiol 2007; 48: 381-390.

[4] Alm E, Arkin AP . Biological networks. Curr Opin Struct Biol 2003; 13: 193-202.

[5] Luo F, Yang Y, Zhong J et al. Constructing gene co-expression networks and predicting functions of unknown genes by random matrix theory. BMC Bioinformatics 2007; 8: 299.

[6] Ruan J, Perez J, Hernandez B, Sunter G, Sponsel VM . Systematic construction and analysis of co-expression networks for identification of functional modules and cis-regulatory elements. In: Proceedings of the 9th International Workshop on Data Mining in Bioinformatics, in Conjuction with Sigkdd’10; 25-28 July 2010; Washington, DC, USA. 2010, pp. 15-24; Arlington, Virginia, USA: ACM press.

[7] Mao L, van Hemert JL, Dash S, Dickerson JA . Arabidopsis gene co-expression network and its functional modules. BMC Bioinformatics 2009; 10: 346.

[8] Han X, Yin L, Xue H . Co-expression analysis identifies CRC and AP1 the regulator of Arabidopsis fatty acid biosynthesis. J Integr Plant Biol 2012; 54: 486-499.

[9] Heyndrickx KS, Vandepoele K . Systematic identification of functional plant modules through the integration of complementary data sources. Plant Physiol 2012; 3: 884-901.

[10] Spangler JB, Ficklin SP, Luo F et al. Conserved non-coding regulatory signatures in Arabidopsis co-expressed gene modules. PloS ONE 2012; 7: e45041.

[11] Wang S, Yin Y, Ma Q et al. Genome-scale identification of cell-wall related genes in Arabidopsis based on co-expression network analysis. BMC Plant Biol 2012; 12: 138.

[12] Zhang L, Yu S, Zuo K et al. Identification of gene modules associated with drought response in rice by network-based Analysis. PloS ONE 2012; 7: e33748.

[13] Ashutosh P, Prashant M, Mohd PK et al. Co-expression of Arabidopsis transcription factor, AtMYB12, and soybean isoflavone synthase, GmIFS1, genes in tobacco leads to enhanced biosynthesis of isoflavones and flavonols resulting in osteoprotective activity. Plant Biotechnol J 2014; 12: 69-80.

[14] Stuart JM, Segal E, Koller D et al. A gene-coexpression network for global discovery of conserved genetic modules. Science 2003; 302: 249-255.

[15] Jaillon O, Aury JM, Noel B et al. The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature 2007; 449: 463-467.

[16] Wang M, Vannozzi A, Wang G et al. Genome and transcriptome analysis of the grapevine (Vitis vinifera L.) WRKY gene family. Hort Res 2014; 1: 1-16.

[17] Ruan J, Zhang W . Identifying network communities with a high resolution. Phys Rev E 2008; 77: 016104.

[18] Maere S, Heymans K, Kuiper M . BiNGO: a Cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics 2005; 21: 3448-3449.

[19] Obayashi T, Kinoshita K, Nakai K et al. ATTED-II: a database of co-expressed genes and cis elements for identifying co-regulated gene groups in Arabidopsis. Nucleic Acids Res 2007; 35: D863-D869.

[20] Shannon P, Markiel A, Ozier O et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res 2003; 13: 2498-2504.

[21] Wang J, Yuan H, Tadesse MG et al. Integration of multiple data sources for identifying functional modules using Bayesian network. In: Proceedings of IEEE International Workshop on Genomic Signal Processing and Statistics; 2-4 December; Washington, DC, USA. 2012, pp. 13-17; Piscataway, New Jersey, USA: IEEE press.

[22] Ulitsky I, Shamir R . Identification of functional modules using network topology and high-throughput data. BMC Syst Biol 2007; 1: 8.

[23] Murashige T, Skoog F . A revised medium for rapid growth and bio assays with tobacco tissue cultures. Physiol Plantarum 1962; 15: 473-497.

[24] Zhang J, Du X, Wang Q et al. Expression of pathogenesis related genes in response to salicylic acid, methyl jasmonate and 1-aminocyclopropane-1-carboxylic acid in Malus hupehensis (Pamp.) Rehd. BMC Res Notes 2010; 3: 208.

[25] Ye X, Yuan S, Guo H et al. Evolution and divergence in the coding and promoter regions of the Populus gene family encoding xyloglucan endotransglycosylase/hydrolases. Tree Genet Genomes 2012; 8: 177-194.

[26] Jacoby WG . LOESS: a nonparametric, graphical tool for depicting relationships between variables. Elect Stud 2000; 19: 577-613.

[27] Hagquist C, Stenbeck M . Goodness of fit in regression analysis-R² and G² reconsidered . Qual Quant 1998; 32: 229-245.

[28] Norusis MJ . IBM SPSS Statistics 19 Advanced Statistical Procedures Companion. Upper Saddle River, NJ: Pearson Education, Prentice Hall Press, 2012.

[29] Timperio AM, Egidi MG, Zolla L . Proteomics applied on plant abiotic stresses: role of heat shock proteins (HSP). J Proteomics 2008; 71: 391-411.

[30] Kotak S, Larkindale J, Lee U et al. Complexity of the heat stress response in plants. Curr Opin Plant Biol 2007; 10: 310-316.

[31] Suzuki N, Bajad S, Shuman J et al. The transcriptional co-activator MBF1c is a key regulator of thermotolerance in Arabidopsis thaliana. J Biol Chem 2008; 283: 9269-9275.

[32] Taji T, Ohsumi C, Iuchi S et al. Important roles of drought-and cold-inducible genes for galactinol synthase in stress tolerance in Arabidopsis thaliana. Plant J 2002; 29: 417-426.

[33] Liu GT, Wang JF, Cramer G et al. Transcriptomic analysis of grape (Vitis vinifera L.) leaves during and after recovery from heat stress. BMC Plant Biol 2012; 12: 174.

[34] Craig E, Huang P, Aron R et al. The diverse roles of J-proteins, the obligate Hsp70 co-chaperone. Rev Physiol Biochem Pharmacol 2006; 36: 1-21.

[35] Panikulangara TJ, Eggers-Schumacher G, Wunderlich M et al. Galactinol synthase1. A novel heat shock factor target gene responsible for heat-induced synthesis of raffinose family oligosaccharides in Arabidopsis. Plant Physiol 2004; 136: 3148-3158.

[36] Wei H, Persson S, Mehta T et al. Transcriptional coordination of the metabolic network in Arabidopsis. Plant Physiol 2006; 142: 762-774.

[37] Lajoie M, Gascuel O, Lefort V et al. Computational discovery of regulatory elements in a continuous expression space. Genome Biol 2012; 13: 1-17.

[38] Allocco DJ, Kohane IS, Butte AJ . Quantifying the relationship between co-expression, co-regulation and gene function. BMC Bioinformatics 2004; 5: 18.

[39] Hu H, Yan X, Huang Y, Han J, Zhou XJ . Mining coherent dense subgraphs across massive biological networks for functional discovery. Bioinformatics 2005; 21: i213-i221.

[40] Mittler R, Finka A, Goloubinoff P . How do plants feel the heat? Trends Biochem Sci 2012; 37: 118-125.

[41] Wang W, Vinocur B, Shoseyov O, Altman A . Role of plant heat-shock proteins and molecular chaperones in the abiotic stress response. Trends Plant Sci 2004; 9: 244-252.

[42] Georgopoulos C, Welch W . Role of the major heat shock proteins as molecular chaperones. Annul Rev Cell Biol 1993; 9: 601-634.

[43] Janska A, Marsik P, Zelenkova S, Ovesna J . Cold stress and acclimation-what is important for metabolic adjustment? Plant Biol 2010; 12: 395-405.

[44] Kosova K, Vitamvas P, Urban MO, Prasil IT . Plant proteome responses to salinity stress-comparison of glycophytes and halophytes. Funct Plant Biol 2013; 40: 775-786.

[45] Grigorova B, Vaseva I, Demirevska K, Feller U . Combined drought and heat stress in wheat: changes in some heat shock proteins. Biol Plantarum 2011; 55: 105-111.

[46] Xoconostle-Cazares B, Ramirez-Ortega FA, Flores-Elenes L, Ruiz-Medrano R . Drought tolerance in crop plants. Am J Plant Physiol 2010; 5: 241-256.

[47] Carvalho LC, Vilela BJ, Mullineaux PM, Amancio S . Comparative transcriptomic profiling of Vitis vinifera under high light using a custom-made array and the Affymetrix GeneChip. Mol Plant 2011; 4: 1038-1051.

[48] Sokolowska I, Woods AG, Wagner J et al. Mass spectrometry for proteomics-based investigation of oxidative stress and heat shock proteins. In: Oxidative Stress: Diagnostics, Prevention, Therapy. Washington, DC: ACS, 2011: 369-411.

[49] Hossain Z, Komatsu S . Contribution of proteomic studies towards understanding plant heavy metal stress response. Front Plant Sci 2012; 3: 1-12.