Constructing Integrated Networks for Identifying New Secondary Metabolic Pathway Regulators in Grapevine: Recent Applications and Future Opportunities
- 1Ecology and Evolution, Research School of Biology, Australian National University, Acton, ACT, Australia
- 2Centre for Research in Agricultural Genomics, CSIC-IRTA-UAB-UB, Barcelona, Spain
Representing large biological data as networks is becoming increasingly adopted for predicting gene function while elucidating the multifaceted organization of life processes. In grapevine (Vitis vinifera L.), network analyses have been mostly adopted to contribute to the understanding of the regulatory mechanisms that control berry composition. Whereas, some studies have used gene co-expression networks to find common pathways and putative targets for transcription factors related to development and metabolism, others have defined networks of primary and secondary metabolites for characterizing the main metabolic differences between cultivars throughout fruit ripening. Lately, proteomic-related networks and those integrating genome-wide analyses of promoter regulatory elements have also been generated. The integration of all these data in multilayered networks allows building complex maps of molecular regulation and interaction. This perspective article describes the currently available network data and related resources for grapevine. With the aim of illustrating data integration approaches into network construction and analysis in grapevine, we searched for berry-specific regulators of the phenylpropanoid pathway. We generated a composite network consisting of overlaying maps of co-expression between structural and transcription factor genes, integrated with the presence of promoter cis-binding elements, microRNAs, and long non-coding RNAs (lncRNA). This approach revealed new uncharacterized transcription factors together with several microRNAs potentially regulating different steps of the phenylpropanoid pathway, and one particular lncRNA compromising the expression of nine stilbene synthase (STS) genes located in chromosome 10. Application of network-based approaches into multi-omics data will continue providing supplementary resources to address important questions regarding grapevine fruit quality and composition.
Network-based Approaches Into Omics Data
Complex biological processes can be studied from a “multi-omics” perspective thanks to the recent improvements in genome-wide techniques and systems biology approaches. Each omics data type is particularly useful in elucidating the constituents and function of a particular cellular domain. Together, they constitute layers of biological complexity. Genomic data generated from genome sequencing projects are commonly used to ascribe molecular function and biological processes based on sequence similarity, while transcriptomics and metabolomics data typically provide a global “snapshot” of gene expression and metabolite dynamics in various biological contexts.
For many omics data, interactions/associations between molecules can be represented as networks, where nodes (genes, proteins, metabolites) are connected by edges. These edges denote an association often inferred from correlational and informational theoretic measures such as Pearson correlation coefficient (PCC) and mutual information (MI), respectively. In the case of gene co-expression networks (GCNs), edges represent similar gene expression behaviors. Based on the “guilt by association” principle, genes involved in related processes share similar gene expression dynamics across a wide range of experiments (Wolfe et al., 2005). However, as functional information could be delimited to a reduced number of interactions within a gene network (Gillis and Pavlidis, 2012), subsequent targeted gene characterizations are needed to prove these relationships. Whether the function of a network is dependent or not on specific interactions, GCN analysis have proven to be a powerful tool for inferring gene function and coordinated biological processes related to plant metabolism (Persson et al., 2005; Itkin et al., 2013).
Other forms of networks constructed from omics datasets do not necessarily rely on abundance or expression levels to establish node relationships. For example, protein-protein interaction networks describe physically interacting protein pairs identified from high-throughput yeast two-hybrid screens (e.g., Arabidopsis Interactome Mapping Consortium, 2011). Also, genome-wide location studies (i.e., by using ChIP-Seq) allow determining regulatory networks for transcription factors (TF) and other DNA-binding proteins. These TF-binding networks have led to the identification of novel components and of new connections that alter the network diagrams originally drawn by genetic and molecular analyses (reviewed by Ferrier et al., 2011).
Recent Application of Metabolite Network-Based Approaches in Grapes
Studies utilizing networks constructed from omics data profiled in the berry are continuously increasing (Table 1). Network analyses involving metabolite datasets (primary and/or secondary metabolites) are by far the most reported. These studies have included networks inferred from single contexts such as berry development and ripening (Zamboni et al., 2010; Dai et al., 2013; Wang et al., in press), or in combination with other factors including environmental influence (Guan et al., 2016; Savoi et al., 2016; Reshef et al., 2017), and/or cultivar differences (Degu et al., 2014; Cuadros-Inostroza et al., 2016). Network topology has also been investigated in detail to reveal critical metabolites and their regulation. For instance, Cuadros-Inostroza et al. (2016) showed that an increase in network connectedness and density (especially regarding primary metabolites) became prevalent at specific berry developmental stages such as fruit set and veraison (i.e., the onset of ripening). The same study, in concordance with Degu et al. (2014), highlighted that berry-metabolite networks from different cultivars could possess contrasting network topologies, albeit with overall network connections generally maintained. Metabolite networks from the cultivars cv. “Merlot” (Cuadros-Inostroza et al., 2016) and cv. “Shiraz” (Degu et al., 2014) were consistently denser to that off cv. “Cabernet Sauvignon.”
Table 1. Studies of grape berry development and composition involving molecular networks approaches.
Rewiring of berry metabolite networks under different environmental conditions or perturbations such as drought (Savoi et al., 2016) and sunlight exposure (Reshef et al., 2017) have also been reported. These studies have shown that higher network connectivity is commonly observed in perturbed networks. Such property could be associated to a tighter metabolic control of the metabolic pathways under investigation. Such is seen for phenylpropanoid and volatile organic compounds (VOC) in berries under prolonged drought compared to non-stress berries (Savoi et al., 2016). Similarly, primary metabolite networks encompassing compounds related to glycolysis, the TCA cycle, and amino acid metabolism showed higher network connectivity in shaded berries compared to fully exposed berries (Reshef et al., 2017).
Some metabolic-network studies have shown that certain metabolites (or classes) could act as important switches in the developmental regulation of metabolism during berry growth and ripening, given their high centrality (number of connections) or degree scores in their network. Dai et al. (2013) showed that trehalose-6-phosphate appeared to be the most connected compound in the primary metabolite network of cv. “Cabernet Sauvignon” grapes, with significant partial correlations to sugar metabolism, glycolysis, and TCA cycle intermediates. Altogether, these compounds may be implicated in coordinating metabolite dynamics during berry development. One recent study highlighted fucose as critical for coordinating metabolic regulation in a stage-specific manner, thus deprioritizing the importance of sugars such as glucose, fructose, and sucrose as a function of network centrality measure (Cuadros-Inostroza et al., 2016). These findings further demonstrate the complexity of berry metabolic regulation during development and ripening.
Gene Co-expression Networks to Study Grape Berry Ripening
The increased ease of transcriptome profiling, combined with availability of datasets shared by the grapevine research community in public repositories, has led to increased attention in the use of gene co-expression networks (GCNs) in the study of berry development and metabolism. GCNs can be classified into “condition-dependent” and “condition-independent” categories (Usadel et al., 2009). In grapes, several studies have focused on condition-independent GCNs (encompassing different cultivars, tissues, developmental stages, stress and vineyard management treatments) as it provides a more convenient and representative (albeit “static”) relationship overview (Table 1). This approach has been useful for ascribing the most representative biological functions of the 134 grapevine R2R3-MYB transcription factors based on their top 100 co-expressed genes (Wong et al., 2016), where VviMYB13 (close homolog of VviMYB14 and VviMYB15) was identified as an additional STILBENE SYNTHASE regulator acting in a tissue- and/or stress-dynamic manner.
Platforms such as the ViTis Co-expression DataBase (VTCdb; Wong et al., 2013) and VESPUCCI (Moretto et al., 2016) have been successfully exploited to study the extent of transcription factor regulatory networks, providing support for targeted functional studies. Such is the case for the bZIP TF VvibZIPC22, which is involved in the regulation of flavonoid biosynthesis in grapes and may be also implicated in carbohydrate and amino acid metabolism, as inferred from VESPUCCI (Malacarne et al., 2016). Two other bZIP TFs (VviHY5 and VviHYH) were shown to co-operatively mediate flavonol accumulation in grapes in response to sunlight and ultraviolet radiation exposure (Loyola et al., 2016). As inferred from VTCdb and GCN analysis, these regulators were potentially implicated in carbohydrate and isoprenoid metabolism in addition to the control of the flavonoid pathway. Similarly, the involvement of the grapevine VviWRKY26 in the regulation of vacuolar transport and flavonoid biosynthesis was demonstrated using a combination of transcriptomic approaches including GCNs (Amato et al., 2017).
Condition-dependent GCNs have been constructed from tissue- or stress-specific datasets, including berry (Zamboni et al., 2010; Palumbo et al., 2014) or abiotic and biotic stresses (Wong et al., 2017). These GCNs provide several advantages over condition-independent networks as inferring gene function is largely simplified, providing a more “dynamic” overview of gene relationships that otherwise could be enhanced or lost in certain conditions (Obayashi et al., 2011). One example of a condition-specific GCN involves the study of the transcriptomes of five black-skinned cultivars across four berry phenological stages (Palumbo et al., 2014). The authors identified “fight-club” nodes and “switch” genes, having the latter unique expression profiles and network topological properties, such as a marked negative correlation connectivity to both neighboring genes and genes grouped to other modules in the network. Genes associated with transcription factor activity; cell wall modification and carbohydrate and secondary metabolism were found as candidate master regulators, potentially inducing large-scale transcriptome reprogramming during berry development (Palumbo et al., 2014).
Finally, miRNA and siRNA-mediated gene regulatory networks have also been constructed from high-throughput small RNA and degradome sequencing and computational target prediction methods (Zhang et al., 2012; Belli Kullan et al., 2015). These networks (not relying in abundance or expression levels) revealed novel modules such as miR156/miR172 regulatory circuits and VviTAS3/4 regulatory cascades, which are implicated in regulating plant growth and development and in the control of flavonoid biosynthesis, respectively.
Toward the Integration of Multi-Omics Data in Grapes
Although individual omic network methods have been widely used, a shift toward multi-omics data and integration is increasingly being adopted in plant biology (Proost and Mutwil, 2016), including grapevine (Table 1). Integration approaches allow building complex maps of molecular regulation and interaction. By these means, complex traits from these networks can be assessed (e.g., plasticity and evolution).
The first systems level study in grapes leveraged transcriptomic, metabolomic, and proteomic technologies to understand berry development and the postharvest withering process (i.e., controlled dehydration) in cv. “Corvina” grapes (Zamboni et al., 2010). Using a combination of hypothesis-free and -driven integration approaches, the authors were able to tease out putative berry stage-specific functional networks. As an outcome, a fully integrated network related to the withering process revealed key phenylpropanoid and stress-responsive genes (i.e., biotic, osmotic, and oxidative), together with proteins involved in oxidative- and osmotic-stress, and secondary metabolites such as acylated anthocyanins and stilbenes. Recently, integration of berry metabolome (primary and secondary) and proteome networks encompassing 12 developmental stages revealed a greater propensity of an energy-linked metabolism in berries prior to veraison (Wang et al., in press). These observations corroborated earlier studies (Dai et al., 2013; Cuadros-Inostroza et al., 2016), demonstrating that pronounced changes in the berry occurs before veraison, characterized by a reduction of many early accumulating primary metabolites. Interestingly, the integrated network also revealed several modules with high node degree for many metabolites (amino acids and organic acids) and corresponding enzymes catalyzing their synthesis (Wang et al., in press).
Characterizing genes that regulate the accumulation of secondary metabolites throughout fruit ripening is key for improving quality traits and for predicting plant behavior in response to the environment. In this sense, transcript-metabolite associations have been used to prioritize candidate genes important for determining berry quality parameters under adverse environmental conditions (Savoi et al., 2016). Integrated transcript-metabolite networks encompassing monoterpenes that are both ripening-related and drought-modulated (e.g., linalool, nerol, α-terpineol) revealed many highly co-regulated transcripts to be involved in terpene and lipid metabolism. The authors further highlighted VviMYB24 as a promising regulatory candidate for monoterpene biosynthesis given consistent correlations with all three monoterpenes in their study.
Cis-regulatory element-driven networks have been recently constructed using integrated information of promoter CRE structure and network connectivity (Wong et al., 2017). Numerous CRE-driven modules inferred using condition-dependent GCNs (development-dynamic and stress-specific) highlighted roles in stress response (e.g., to drought and pathogens) and developmental processes (e.g., fruit ripening). For example, GCC-core sub-modules contained many genes that were highly induced in berries and leaves infected with fungi (Wong et al., 2017).
Cis-regulatory element enrichment maps or transcript information for miRNA target enrichment analysis can be easily integrated into plant GCNs. This approach has been used to prioritize target genes of the entire grape R2R3-MYB family (Wong et al., 2016) and also to explain the expression responses of module genes under prolonged drought stress in berries (Savoi et al., 2016). Enrichment for miRNA targets within GCNs has suggested a pivotal role of these molecules in regulating the expression of “switch” genes in a stage-specific manner (Palumbo et al., 2014). Finally, aggregating several networks into a community network can also be advantageous to effectively reveal discrepancies between individual networks while highlighting associations common across individual networks (Proost and Mutwil, 2016). This approach has been used by Loyola et al. (2016) to identify a set of high confidence targets of HY5 and HYH given by the combination of microarray and RNA-Seq data with genome-wide promoter inspections. It is noteworthy that “condition-independent” and “condition-dependent” approaches are still useful for providing a preliminary insight into co-expression relationships in grapes.
An Illustration for the Integration of Multilayered Networks for Dissecting the Complexity of the Berry's Phenylpropanoid Composition
In grapes, phenylpropanoids influence their organoleptic properties and beneficial attributes to human health, highlighting the importance of their study. Several reports have demonstrated the complex nature of secondary metabolism in grapevine, both at the level of chemical composition and genetic regulation (Dal Santo et al., 2013; Costantini et al., 2015; Malacarne et al., 2015). Among the many phenylpropanoid compounds that influence the quality of grapes and wines, some of the most important are flavonoids (anthocyanins, flavonols and tannins) and stilbenes. These compounds accumulate in a temporal and compartmentalized manner and numerous regulators of their accumulation have been characterized to date (Reviewed by Kuhn et al., 2014; Matus, 2016). One strikingly relevant feature of the grapevine genome is that wine quality-related gene families are expanded in gene number (Martin et al., 2010; Vannozzi et al., 2012), including those related with transcription factor activity (Matus et al., 2008; Wong et al., 2016). Genomics and transcriptomics data originated from these and others studies suggest that the regulation of secondary metabolism in grape is a much more complex trait compared to plant model species. As large-scale omics data are periodically accumulating; there is an enormous potential for gene discovery in relation to grape secondary metabolic pathways.
To demonstrate how various biological networks can be integrated to study berry's phenylpropanoid composition, we gathered networks generated from gene co-expression analyses, predicted miRNA-gene and long non-coding RNA (lncRNA) -gene interactions. First, we re-analyzed a comprehensive berry ripening RNA-Seq transcriptome dataset (five black-skinned cultivars sampled at four developmental stages; Palumbo et al., 2014) and constructed a ripening-specific gene co-expression network (PCC > |0.8|). This ripening-specific GCN was then used as a basis for lncRNA-gene network, which consisted of predicted lncRNAs (Vitulo et al., 2014) that showed strong correlation with a putative “interacting” gene (PCC > | 0.8 |) that was co-located within 100 kb flanking the lncRNA position. Using a comprehensive catalog of grapevine miRNAs (Belli Kullan et al., 2015; Pulvirenti et al., 2015), we also reanalyzed potential miRNA-mRNA interactions using psRNATarget with default parameters (Dai and Zhao, 2011). As the interpretation of each network at a global scale is out of the scope of this perspective, we focused our attention on the early phenylpropanoid and flavonoid (ePP and Fla) pathways and on the potential regulatory genes and their interactions (among genes, miRNAs, and lncRNAs). The resulting network is composed of 112 ePP/Fla pathway genes (differentially expressed during berry development and ripening) together with five miRNA and 14 predicted grapevine lncRNAs (Figure 1). GCN analysis revealed a strong co-regulation within early phenylpropanoid and flavonoid pathway genes maintaining few connections between both sub-pathway genes during the course of berry development and ripening.
Figure 1. An integrated network for phenylpropanoid regulation in the grape berry. Circle, diamond, square, and triangle nodes represent enzyme-coding genes, transcription factors (TFs), long non-coding RNAs (lncRNAs), and micro RNAs (miRNAs), respectively. Gray and light blue solid edges connecting circle and square nodes depict positive and negative RNA-Seq-derived co-expression correlations, respectively. Dashed and dotted lines depict predicted miRNA-gene interaction and lncRNA-gene co-location (within 100 kb). Clusters I, II and III connect genes belonging to the early phenylpropanoid and flavonoid (ePP and Fla) sub-pathways. Cluster IV shows a dense group containing predominantly PHENYLALANINE AMMONIA-LYASE (PAL) and STILBENE SYNTHASE (STS) genes that are largely co-expressed with positive co-expression correlations. Purple edges represent positive co-expression correlations between TFs and early phenylpropanoid pathway genes. Pie chart colors represent the presence of selected TF-binding sites (based on cis-regulatory element enrichment analysis) in promoter regions of the corresponding enzyme-coding genes. Light blue border edges depict STS genes located in chromosome 10.
Three clusters (I, II and III) were observed for Fla pathway genes sharing many positive correlations within each group (Figure 1). Cluster I contained genes mainly involved in the regulation of anthocyanin accumulation such as five flavonoid-3′,5′-hydroxylases (F3′5′H), two anthocyanin-o-methyltransferases (AOMT1-2), the UDP-GLUCOSE:FLAVONOID 3-O-GLUCOSYLTRANSFERASE (UFGT) and ANTHOCYANIN-3-O-GLUCOSIDE-6′′-O-ACYLTRANSFERASE (3AT). Cluster II consisted of genes encoding proanthocyanidin biosynthesis genes including three predicted galloyl glucosyltransferases, ANTHOCYANIDIN REDUCTASE (ANR), and LEUCOANTHOCYANIDIN REDUCTASE (LAR), as well as upstream flavonoid pathway genes such as CHALCONE SYNTHASE (CHS) and CHALCONE ISOMERASE (CHI). One predicted antisense lncRNA (VIT_203s0180n00020) collocated (within 50 kB) and positively correlated with one galloyl glucosyltransferase gene (VIT_03s0180g00200). This cluster also contained genes encoding one 4-coumarate-Co-A ligase (4CL), two flavonol synthases (FLS4-5), one flavanone-3-hydroxylase (F3H), and one caffeic acid 3-o-methyltransferase (COMT), all of which were negatively correlated with genes from cluster I. Furthermore, grapevine miRNAs miR169r/t and grape-m0534 were predicted to target 4CL and FLS4, respectively. Several genes belonging to cluster II and I shared negative correlations (light blue solid edges, Figure 1). This separation is evident whereby the majority of genes from cluster I are ripening-specific (i.e., upregulated from veraison onwards), while many genes from cluster II are mostly expressed during the early-to-mid stages of berry development (and subsequently downregulated as ripening progresses).
As there is much less evidence in the regulation of the early phenylpropanoid and stilbene sub-pathways compared to the regulation of flavonoid biosynthesis, we focused our attention on a fourth, highly connected cluster (IV) holding strong positive correlations within and between the two large PHENYLALANINE AMMONIA-LYASE (PAL) and STILBENE SYNTHASE (STS) gene families (Figure 1). Two cinnamate-4-hydroxylases (C4H) also shared many strong positive correlations with PAL genes and one 4CL was positively correlated with many STS encoding genes. Gene expressions within this cluster were mainly late-ripening specific, with many of them peaking at harvest. Promoters from cluster IV were highly enriched for cis-regulatory elements including those for R2R3-MYB, AP2/ERF, WRKY, bHLH, and bZIP TF binding. In particular, the MYB binding site CCWACC was present in one CCoAMT, two C4H, 10 PAL, and 27 STS genes. The potential regulation of these genes by MYB transcription factors is supported by recent studies showing that several grapevine MYBs may have regulatory roles controlling the levels of small weight phenylpropanoids and stilbenes (Höll et al., 2013; Cavallini et al., 2015). Our approach is novel in suggesting the regulatory roles by other TF families such as WRKY and AP2/ERF. For example, strong co-regulation of nine WRKY TF to 11 PAL and 44 STS genes coincided with the presence of WRKY cis-regulatory elements in many PAL and STS genes. Interestingly, one of the four predicted intergenic lncRNAs (VIT_210s0042n00100) was co-located and strongly co-regulated with all nine STS positioned on chromosome 10. Recent evidence from several functionally characterized lncRNAs in animals and plants suggest that lncRNAs could operate as decoys, guides, signals, and scaffolds, acting as single molecules or complexes regulating pre- and post-transcriptional processes (Wang and Chang, 2011). As such, our observation raises the plausibility of a large-scale regulatory function between this lncRNA and STS genes. This STS-associated lncRNA may fulfill combinatorial roles for the fine-regulation of multiple STS, as signals for transcription activity in a stage-specific way or as guides for chromatin modifiers to the cluster of tandem-positioned STS of chromosome 10, potentially modulating DNA accessibility.
Multi-omics studies incorporating systems biology approaches in grapevine have facilitated the identification of new grape secondary metabolism regulators and have helped in the characterization of genome-wide responses to environmental factors. These studies have brought knowledge and new tools to understand how to modify and improve grape's quality. Additional efforts will still be needed to map protein-DNA and protein-protein landscapes at a large scale. Also, DNAse I hypersensitivity mapping could be useful to identify pioneering transcription factors controlling grape and wine quality traits.
JTM conceived the article and planned its structure. DW and JTM searched and discussed the literature and wrote the manuscript. DW generated new network data. All authors have read and approved the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors wish to acknowledge Dr. Jason Argyris (Centre for Research in Agricultural Genomics, CRAG) for critically reviewing this work.
Amato, A., Cavallini, E., Zenoni, S., Finezzo, L., Begheldo, M., Ruperti, B., et al. (2017). A grapevine TTG2-like WRKY transcription factor is involved in regulating vacuolar transport and flavonoid biosynthesis. Front. Plant Sci. 7:1979. doi: 10.3389/fpls.2016.01979
Belli Kullan, J., Lopes Paim Pinto, D., Bertolini, E., Fasoli, M., Zenoni, S., Tornielli, G. B., et al. (2015). miRVine: a microRNA expression atlas of grapevine based on small RNA sequencing. BMC Genomics 16:393. doi: 10.1186/s12864-015-1610-5
Cavallini, E., Matus, J. T., Finezzo, L., Zenoni, S., Loyola, R., Guzzo, F., et al. (2015). The phenylpropanoid pathway is controlled at different branches by a set of R2R3-MYB C2 repressors in grapevine. Plant Physiol. 167, 1448–1470. doi: 10.1104/pp.114.256172
Costantini, L., Malacarne, G., Lorenzi, S., Troggio, M., Mattivi, F., Moser, C., et al. (2015). New candidate genes for the fine regulation of the colour of grapes. J. Exp. Bot. 66, 4427–4440. doi: 10.1093/jxb/erv159
Cuadros-Inostroza, A., Ruíz-Lara, S., González, E., Eckardt, A., Willmitzer, L., and Peña-Cortés, H. (2016). GC-MS metabolic profiling of Cabernet Sauvignon and Merlot cultivars during grapevine berry development and network analysis reveals a stage- and cultivar-dependent connectivity of primary metabolites. Metabolomics 12, 1–17. doi: 10.1007/s11306-015-0927-z
Dai, Z. W., León, C., Feil, R., Lunn, J. E., Delrot, S., and Gomès, E. (2013). Metabolic profiling reveals coordinated switches in primary carbohydrate metabolism in grape berry (Vitis vinifera L.), a non-climacteric fleshy fruit. J. Exp. Bot. 64, 1345–1355. doi: 10.1093/jxb/ers396
Degu, A., Hochberg, U., Sikron, N., Venturini, L., Buson, G., Ghan, R., et al. (2014). Metabolite and transcript profiling of berry skin during fruit development elucidates differential regulation between Cabernet Sauvignon and Shiraz cultivars at branching points in the polyphenol pathway. BMC Plant Biol. 14:188. doi: 10.1186/s12870-014-0188-4
Ferrier, T., Matus, J. T., Jin, J., and Riechmann, J. L. (2011). Arabidopsis paves the way: genomic and network analyses in crops. Curr. Opin. Biotechnol. 22, 260–270. doi: 10.1016/j.copbio.2010.11.010
Guan, L., Dai, Z., Wu, B. H., Wu, J., Merlin, I., Hilbert, G., et al. (2016). Anthocyanin biosynthesis is differentially regulated by light in the skin and flesh of white-fleshed and teinturier grape berries. Planta 243, 23–41. doi: 10.1007/s00425-015-2391-4
Höll, J., Vannozzi, A., Czemmel, S., D'Onofrio, C., Walker, A. R., Rausch, T., et al. (2013). The R2R3-MYB transcription factors MYB14 and MYB15 regulate stilbene biosynthesis in Vitis vinifera. Plant Cell 25, 4135–4149. doi: 10.1105/tpc.113.117127
Itkin, M., Heinig, U., Tzfadia, O., Bhide, A. J., Shinde, B., Cardenas, P. D., et al. (2013). Biosynthesis of antinutritional alkaloids in solanaceous crops is mediated by clustered genes. Science 341, 175–179. doi: 10.1126/science.1240230
Loyola, R., Herrera, D., Mas, A., Wong, D. C. J., Höll, J., Cavallini, E., et al. (2016). The photomorphogenic factors UV-B RECEPTOR 1, ELONGATED HYPOCOTYL 5, and HY5 HOMOLOGUE are part of the UV-B signalling pathway in grapevine and mediate flavonol accumulation in response to the environment. J. Exp. Bot. 67, 5429–5445. doi: 10.1093/jxb/erw307
Malacarne, G., Costantini, L., Coller, E., Battilana, J., Velasco, R., Vrhovsek, U., et al. (2015). Regulation of flavonol content and composition in (Syrah×Pinot Noir) mature grapes: integration of transcriptional profiling and metabolic quantitative trait locus analyses. J. Exp. Bot. 66, 4441–4453. doi: 10.1093/jxb/erv243
Malacarne, G., Coller, E., Czemmel, S., Vrhovsek, U., Engelen, K., Goremykin, V., et al. (2016). The grapevine VvibZIPC22 transcription factor is involved in the regulation of flavonoid biosynthesis. J. Exp. Bot. 67, 3509–3522. doi: 10.1093/jxb/erw181
Martin, D. M., Aubourg, S., Schouwey, M. B., Daviet, L., Schalk, M., Toub, O., et al. (2010). Functional annotation, genome organization and phylogeny of the grapevine (Vitis vinifera) terpene synthase gene family based on genome assembly, FLcDNA cloning, and enzyme assays. BMC Plant Biol. 10:226. doi: 10.1186/1471-2229-10-226
Matus, J. T. (2016). Transcriptomic and metabolomic setworks in the grape berry illustrate that it takes more than flavonoids to fight against ultraviolet radiation. Front. Plant Sci. 7:1337. doi: 10.3389/fpls.2016.01337
Matus, J. T., Aquea, F., and Arce-Johnson, P. (2008). Analysis of the grape MYB R2R3 subfamily reveals expanded wine quality-related clades and conserved gene structure organization across Vitis and Arabidopsis genomes. BMC Plant Biol. 8:83. doi: 10.1186/1471-2229-8-83
Moretto, M., Sonego, P., Pilati, S., Malacarne, G., Costantini, L., Grzeskowiak, L., et al. (2016). VESPUCCI: exploring patterns of gene expression in grapevine. Front. Plant Sci. 7:633. doi: 10.3389/fpls.2016.00633
Obayashi, T., Nishida, K., Kasahara, K., and Kinoshita, K. (2011). ATTED-II updates: condition-specific gene coexpression to extend coexpression analyses and applications to a broad range of flowering plants. Plant Cell Physiol. 52, 213–219. doi: 10.1093/pcp/pcq203
Palumbo, M. C., Zenoni, S., Fasoli, M., Massonnet, M., Farina, L., Castiglione, F., et al. (2014). Integrated network analysis identifies fight-club nodes as a class of hubs encompassing key putative switch genes that induce major transcriptome reprogramming during grapevine development. Plant Cell Online 26, 4617–4635. doi: 10.1105/tpc.114.133710
Persson, S., Wei, H., Milne, J., Page, G. P., and Somerville, C. R. (2005). Identification of genes required for cellulose synthesis by regression analysis of public microarray data sets. Proc. Natl. Acad. Sci. U.S.A. 102, 8633–8638. doi: 10.1073/pnas.0503392102
Pulvirenti, A., Giugno, R., Distefano, R., Pigola, G., Mongiovi, M., Giudice, G., et al. (2015). A knowledge base for Vitis vinifera functional analysis. BMC Syst. Biol. 9(Suppl. 3):S5. doi: 10.1186/1752-0509-9-S3-S5
Reshef, N., Walbaum, N., Agam, N., and Fait, A. (2017). Sunlight modulates fruit metabolic profile and shapes the spatial pattern of compound accumulation within the grape cluster. Front. Plant Sci. 8, 1–20. doi: 10.3389/fpls.2017.00070
Savoi, S., Wong, D. C., Arapitsas, P., Miculan, M., Bucchetti, B., Peterlunger, E., et al. (2016). Transcriptome and metabolite profiling reveals that prolonged drought modulates the phenylpropanoid and terpenoid pathway in white grapes (Vitis vinifera L.). BMC Plant Biol. 16:67. doi: 10.1186/s12870-016-0760-1
Usadel, B., Obayashi, T., Mutwil, M., Giorgi, F. M., Bassel, G. W., Tanimoto, M., et al. (2009). Co-expression tools for plant biology: opportunities for hypothesis generation and caveats. Plant Cell Environ. 32, 1633–1651. doi: 10.1111/j.1365-3040.2009.02040.x
Vannozzi, A., Dry, I. B., Fasoli, M., Zenoni, S., and Lucchin, M. (2012). Genome-wide analysis of the grapevine stilbene synthase multigenic family: genomic organization and expression profiles upon biotic and abiotic stresses. BMC Plant Biol. 12:130. doi: 10.1186/1471-2229-12-130
Vitulo, N., Forcato, C., Carpinelli, E. C., Telatin, A., Campagna, D., D'Angelo, M., et al. (2014). A deep survey of alternative splicing in grape reveals changes in the splicing machinery related to tissue, stress condition and genotype. BMC Plant Biol. 14:99. doi: 10.1186/1471-2229-14-99
Wang, M., Vannozzi, A., Wang, G., Zhong, Y., Corso, M., Cavallini, E., et al. (2015). A comprehensive survey of the grapevine VQ gene family and its transcriptional correlation with WRKY proteins. Front. Plant Sci. 6:417. doi: 10.3389/fpls.2015.00417
Wang, L., Sun, X., Weiszmann, J., and Weckwerth, W. (in press). System-level granger network analysis of integrated proteomic metabolomic dynamics identifies key points of grape berry development at the interface of primary secondary metabolism. Front. Plant Sci. Rev.
Wen, Y. Q., Zhong, G. Y., Gao, Y., Lan, Y. B., Duan, C. Q., and Pan, Q. H. (2015). Using the combined analysis of transcripts and metabolites to propose key genes for differential terpene accumulation across two regions. BMC Plant Biol. 15:240. doi: 10.1186/s12870-015-0631-1
Wolfe, C. J., Kohane, I. S., and Butte, A. J. (2005). Systematic survey reveals general applicability of guilt-by-association within gene coexpression networks. BMC Bioinformatics 6:227. doi: 10.1186/1471-2105-6-227
Wong, D. C., Sweetman, C., Drew, D. P., and Ford, C. M. (2013). VTCdb: a gene co-expression database for the crop species Vitis vinifera (grapevine). BMC Genomics 14:882. doi: 10.1186/1471-2164-14-882
Wong, D. C. J., Schlechter, R., Vannozzi, A., Höll, J., Hmmam, I., Bogs, J., et al. (2016). A systems-oriented analysis of the grapevine R2R3-MYB transcription factor family uncovers new insights into the regulation of stilbene accumulation. DNA Res. 23, 451–466. doi: 10.1093/dnares/dsw028
Wong, D. C. J., Lopez Gutierrez, R., Gambetta, G. A., and Castellarin, S. D. (2017). Genome-wide analysis of cis-regulatory element structure and discovery of motif-driven gene co-expression networks in grapevine. DNA Res. doi: 10.1093/dnares/dsw061. [Epub ahead of print].
Zamboni, A., Di Carli, M., Guzzo, F., Stocchero, M., Zenoni, S., Ferrarini, A., et al. (2010). Identification of putative stage-specific grapevine berry biomarkers and omics data integration into networks. Plant Physiol. 154, 1439–1459. doi: 10.1104/pp.110.160275
Keywords: Vitis vinifera, stilbenes, flavonoids, ripening, systems biology, genome-wide, multi-omics, data integration
Citation: Wong DCJ and Matus JT (2017) Constructing Integrated Networks for Identifying New Secondary Metabolic Pathway Regulators in Grapevine: Recent Applications and Future Opportunities. Front. Plant Sci. 8:505. doi: 10.3389/fpls.2017.00505
Received: 16 February 2017; Accepted: 22 March 2017;
Published: 12 April 2017.
Edited by:Ashraf El-kereamy, University of California, USA
Reviewed by:Sara Zenoni, University of Verona, Italy
Jose A. Casaretto, University of Guelph, Canada
Copyright © 2017 Wong and Matus. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: José Tomás Matus, email@example.com