VTCdb: a gene co-expression database for the crop species Vitis vinifera (grapevine)

CoGTEx: Unscaled system-level coexpression estimation from GTEx data forecast novel functional gene partners

MOTIVATION

Coexpression estimations are helpful for analysis of pathways, cofactors, regulators, targets, and human health and disease. Ideally, coexpression estimations should consider as many diverse cell types as possible and consider that available data is not uniform across tissues. Importantly, the coexpression estimations accessible today are performed on a "tissue level", which is based on cell type standardized formulations. Little or no attention is paid to overall gene expression levels. The tissue-level estimation assumes that variance expression levels are more important than mean expression levels. Here, we challenge this assumption by estimating a coexpression calculation at the "system level", which is estimated without standardization by tissue, and show that it provides valuable information. We made available a resource to view, download, and analyze both, tissue- and system-level coexpression estimations from GTEx human data.

METHODS

GTEx v8 expression data was globally normalized, batch-processed, and filtered. Then, PCA, clustering, and tSNE stringent procedures were applied to generate 42 distinct and curated tissue clusters. Coexpression was estimated from these 42 tissue clusters computing the correlation of 33,445 genes by sampling 70 samples per tissue cluster to avoid tissue overrepresentation. This process was repeated 20 times, extracting the minimum value provided as a robust estimation. Three metrics were calculated (Pearson, Spearman, and G-statistic) in two data processing modes, at the system-level (TPM scale) and tissue levels (z-score scale).

RESULTS

We first validate our tissue-level estimations compared with other databases. Then, by specific analyses in several examples and literature validations of predictions, we show that system-level coexpression estimation differs from tissue-level estimations and that both contain valuable information reflected in biological pathways. We also show that coexpression estimations are associated to transcriptional regulation. Finally, we present CoGTEx, a valuable resource for viewing and analyzing coexpressed genes in human adult tissues from GTEx v8 data. We introduce our web resource to list, view and explore the coexpressed genes from GTEx data.

CONCLUSION

We conclude that system-level coexpression is a novel and interesting coexpression metric capable of generating plausible predictions and biological hypotheses; and that CoGTEx is a valuable resource to view, compare, and download system- and tissue- level coexpression estimations from GTEx data.

AVAILABILITY

The web resource is available at http://bioinformatics.mx/cogtex.

Genome-wide characterisation of HD-Zip transcription factors and functional analysis of PbHB24 during stone cell formation in Chinese white pear (Pyrus bretschneideri)

BACKGROUND

The homodomain-leucine zipper (HD-Zip) is a conserved transcription factor family unique to plants that regulate multiple developmental processes including lignificaion. Stone cell content is a key determinant negatively affecting pear fruit quality, which causes a grainy texture of fruit flesh, because of the lignified cell walls.

RESULTS

In this study, a comprehensive bioinformatics analysis of HD-Zip genes in Chinese white pear (Pyrus bretschneideri) (PbHBs) was performed. Genome-wide identification of the PbHB gene family revealed 67 genes encoding PbHB proteins, which could be divided into four subgroups (I, II, III, and IV). For some members, similar intron/exon structural patterns support close evolutionary relationships within the same subgroup. The functions of each subgroup of the PbHB family were predicted through comparative analysis with the HB genes in Arabidopsis and other plants. Cis-element analysis indicated that PbHB genes might be involved in plant hormone signalling and external environmental responses, such as light, stress, and temperature. Furthermore, RNA-sequencing data and quantitative real-time PCR (RT-qPCR) verification revealed the regulatory roles of PbHB genes in pear stone cell formation. Further, co-expression network analysis revealed that the eight PbHB genes could be classified into different clusters of co-expression with lignin-related genes. Besides, the biological function of PbHB24 in promoting stone cell formation has been demonstrated by overexpression in fruitlets.

CONCLUSIONS

This study provided the comprehensive analysis of PbHBs and highlighted the importance of PbHB24 during stone cell development in pear fruits.

Aggregated gene co-expression networks predict transcription factor regulatory landscapes in grapevine

Gene co-expression networks (GCNs) have not been extensively studied in non-model plants. However, the rapid accumulation of transcriptome datasets in certain species represents an opportunity to explore underutilized network aggregation approaches. In fact, aggregated GCNs (aggGCNs) highlight robust co-expression interactions and improve functional connectivity. We applied and evaluated two different aggregation methods on public grapevine RNA-Seq datasets from three different tissues (leaf, berry, and 'all organs'). Our results show that co-occurrence-based aggregation generally yielded the best-performing networks. We applied aggGCNs to study several transcription factor gene families, showing their capacity for detecting both already-described and novel regulatory relationships between R2R3-MYBs, bHLH/MYC, and multiple specialized metabolic pathways. Specifically, transcription factor gene- and pathway-centered network analyses successfully ascertained the previously established role of VviMYBPA1 in controlling the accumulation of proanthocyanidins while providing insights into its novel role as a regulator of p-coumaroyl-CoA biosynthesis as well as the shikimate and aromatic amino acid pathways. This network was validated using DNA affinity purification sequencing data, demonstrating that co-expression networks of transcriptional activators can serve as a proxy of gene regulatory networks. This study presents an open repository to reproduce networks in other crops and a GCN application within the Vitviz platform, a user-friendly tool for exploring co-expression relationships.

RoseAP: an analytical platform for gene function of Rosa rugosa

Rosa rugosa, a perennial shrub belonging to family Rosaceae, is a well-known ornamental plant. Its petals contain an abundance of essential oils and anthocyanins with enormous economic and health benefits when used as edible or cosmetic ingredients. The whole genome of R. rugosa was sequenced in 2021, which provided opportunities and challenges for gene regulation. However, many gene functions remain unknown. Therefore, an analytical platform named RoseAP (http://www.gzybioinformatics.cn/RoseAP/index.php) for the functional analysis of R. rugosa genes was constructed. It improved the gene annotation rate by integrating and analyzing genomic and transcriptomic datasets. First, 38,815 genes, covering 97.76% of the coding genes, were annotated functionally and structurally using a variety of algorithms and rules. Second, a total of 33 transcriptome samples were integrated, including 23 samples from our lab and 10 samples from the SRA database. A co-expression network containing approximately 29,657 positive or negative gene pairs, covering 74.7% of the coding genes, was constructed based on PCC and MR algorithms. Network analysis revealed that the DFR function was closely related to anthocyanin metabolism. It demonstrated the reliability of the network. Several SAUR genes of R. rugosa shared similar expression patterns. RoseAP was used to determine the sequence, structure, functional annotation, expression profile, regulatory network, and functional modules at the transcriptional and protein levels by inputting gene IDs. In addition, auxiliary analytical tools, including BLAST, gene set enrichment, orthologue conversion, gene sequence extraction, gene expression value extraction, and JBrowse, were utilized. Regular updates to RoseAP are expected to facilitate mining of gene function and promote genetic improvement in R. rugosa.

Transcriptomic and metabolomic integration as a resource in grapevine to study fruit metabolite quality traits

Transcriptomics and metabolomics are methodologies being increasingly chosen to perform molecular studies in grapevine (Vitis vinifera L.), focusing either on plant and fruit development or on interaction with abiotic or biotic factors. Currently, the integration of these approaches has become of utmost relevance when studying key plant physiological and metabolic processes. The results from these analyses can undoubtedly be incorporated in breeding programs whereby genes associated with better fruit quality (e.g., those enhancing the accumulation of health-promoting compounds) or with stress resistance (e.g., those regulating beneficial responses to environmental transition) can be used as selection markers in crop improvement programs. Despite the vast amount of data being generated, integrative transcriptome/metabolome meta-analyses (i.e., the joint analysis of several studies) have not yet been fully accomplished in this species, mainly due to particular specificities of metabolomic studies, such as differences in data acquisition (i.e., different compounds being investigated), unappropriated and unstandardized metadata, or simply no deposition of data in public repositories. These meta-analyses require a high computational capacity for data mining a priori, but they also need appropriate tools to explore and visualize the integrated results. This perspective article explores the universe of omics studies conducted in V. vinifera, focusing on fruit-transcriptome and metabolome analyses as leading approaches to understand berry physiology, secondary metabolism, and quality. Moreover, we show how omics data can be integrated in a simple format and offered to the research community as a web resource, giving the chance to inspect potential gene-to-gene and gene-to-metabolite relationships that can later be tested in hypothesis-driven research. In the frame of the activities promoted by the COST Action CA17111 INTEGRAPE, we present the first grapevine transcriptomic and metabolomic integrated database (TransMetaDb) developed within the Vitis Visualization (VitViz) platform (https://tomsbiolab.com/vitviz). This tool also enables the user to conduct and explore meta-analyses utilizing different experiments, therefore hopefully motivating the community to generate Findable, Accessible, Interoperable and Reusable (F.A.I.R.) data to be included in the future.

Multi-Omics Approaches and Resources for Systems-Level Gene Function Prediction in the Plant Kingdom

In higher plants, the complexity of a system and the components within and among species are rapidly dissected by omics technologies. Multi-omics datasets are integrated to infer and enable a comprehensive understanding of the life processes of organisms of interest. Further, growing open-source datasets coupled with the emergence of high-performance computing and development of computational tools for biological sciences have assisted in silico functional prediction of unknown genes, proteins and metabolites, otherwise known as uncharacterized. The systems biology approach includes data collection and filtration, system modelling, experimentation and the establishment of new hypotheses for experimental validation. Informatics technologies add meaningful sense to the output generated by complex bioinformatics algorithms, which are now freely available in a user-friendly graphical user interface. These resources accentuate gene function prediction at a relatively minimal cost and effort. Herein, we present a comprehensive view of relevant approaches available for system-level gene function prediction in the plant kingdom. Together, the most recent applications and sought-after principles for gene mining are discussed to benefit the plant research community. A realistic tabulation of plant genomic resources is included for a less laborious and accurate candidate gene discovery in basic plant research and improvement strategies.

NetREx: Network-based Rice Expression Analysis Server for abiotic stress conditions

Recent focus on transcriptomic studies in food crops like rice, wheat and maize provide new opportunities to address issues related to agriculture and climate change. Re-analysis of such data available in public domain supplemented with annotations across molecular hierarchy can be of immense help to the plant research community, particularly co-expression networks representing transcriptionally coordinated genes that are often part of the same biological process. With this objective, we have developed NetREx, a Network-based Rice Expression Analysis Server, that hosts ranked co-expression networks of Oryza sativa using publicly available messenger RNA sequencing data across uniform experimental conditions. It provides a range of interactable data viewers and modules for analysing user-queried genes across different stress conditions (drought, flood, cold and osmosis) and hormonal treatments (abscisic and jasmonic acid) and tissues (root and shoot). Subnetworks of user-defined genes can be queried in pre-constructed tissue-specific networks, allowing users to view the fold change, module memberships, gene annotations and analysis of their neighbourhood genes and associated pathways. The web server also allows querying of orthologous genes from Arabidopsis, wheat, maize, barley and sorghum. Here, we demonstrate that NetREx can be used to identify novel candidate genes and tissue-specific interactions under stress conditions and can aid in the analysis and understanding of complex phenotypes linked to stress response in rice. Database URL: https://bioinf.iiit.ac.in/netrex/index.html.

Gene Co-Expression Network Tools and Databases for Crop Improvement

Transcriptomics has significantly grown as a functional genomics tool for understanding the expression of biological systems. The generated transcriptomics data can be utilised to produce a gene co-expression network that is one of the essential downstream omics data analyses. To date, several gene co-expression network databases that store correlation values, expression profiles, gene names and gene descriptions have been developed. Although these resources remain scattered across the Internet, such databases complement each other and support efficient growth in the functional genomics area. This review presents the features and the most recent gene co-expression network databases in crops and summarises the present status of the tools that are widely used for constructing the gene co-expression network. The highlights of gene co-expression network databases and the tools presented here will pave the way for a robust interpretation of biologically relevant information. With this effort, the researcher would be able to explore and utilise gene co-expression network databases for crops improvement.

VviAGL11 self-regulates and targets hormone- and secondary metabolism-related genes during seed development

VviAGL11, the Arabidopsis SEEDSTICK homolog, has been proposed to have a causative role in grapevine stenospermocarpy. An association between a mutation in the coding sequence (CDS) and the seedless phenotype was reported, however, no working mechanisms have been demonstrated yet. We performed a deep investigation of the full VviAGL11 gene sequence in a collection of grapevine varieties belonging to several seedlessness classes that revealed three different promoter-CDS combinations. By investigating the expression of the three VviAGL11 alleles, and by evaluating their ability to activate the promoter region, we observed that VviAGL11 self-activates in a specific promoter-CDS combination manner. Furthermore, by transcriptomic analyses on ovule and developing seeds in seeded and seedless varieties and co-expression approaches, candidate VviAGL11 targets were identified and further validated through luciferase assay and in situ hybridization. We demonstrated that VviAGL11 Wild Type CDS activates Methyl jasmonate esterase and Indole-3-acetate beta-glucosyltransferase, both involved in hormone signaling and Isoflavone reductase, involved in secondary metabolism. The dominant-negative effect of the mutated CDS was also functionally ectopically validated in target induction. VviAGL11 was shown to co-localize with its targets in the outer seed coat integument, supporting its direct involvement in seed development, possibly by orchestrating the crosstalk among MeJA, auxin, and isoflavonoids synthesis. In conclusion, the VviAGL11 expression level depends on the promoter-CDS allelic combination, and this will likely affect its ability to activate important triggers of the seed coat development. The dominant-negative effect of the mutated VviAGL11 CDS on the target genes activation was molecularly validated. A new regulatory mechanism correlating VviAGL11 haplotype assortment and seedlessness class in grapevine is proposed.

A COMPASS for VESPUCCI: A FAIR Way to Explore the Grapevine Transcriptomic Landscape

Successfully integrating transcriptomic experiments is a challenging task with the ultimate goal of analyzing gene expression data in the broader context of all available measurements, all from a single point of access. In its second major release VESPUCCI, the integrated database of gene expression data for grapevine, has been updated to be FAIR-compliant, employing standards and created with open-source technologies. It includes all public grapevine gene expression experiments from both microarray and RNA-seq platforms. Transcriptomic data can be accessed in multiple ways through the newly developed COMPASS GraphQL interface, while the expression values are normalized using different methodologies to flexibly satisfy different analysis requirements. Sample annotations are manually curated and use standard formats and ontologies. The updated version of VESPUCCI provides easy querying and analyzing of integrated grapevine gene expression (meta)data and can be seamlessly embedded in any analysis workflow or tools. VESPUCCI is freely accessible and offers several ways of interaction, depending on the specific goals and purposes and/or user expertise; an overview can be found at https://vespucci.readthedocs.io/.

Gene Co-expression Network Analysis and Linking Modules to Phenotyping Response in Plants

Environmental factors, including different stresses, can have an impact on the expression of genes and subsequently the phenotype and development of plants. Since a large number of genes are involved in response to the perturbation of the environment, identifying groups of co-expressed genes is meaningful. The gene co-expression network models can be used for the exploration, interpretation, and identification of genes responding to environmental changes. Once a gene co-expression network is constructed, one can determine gene modules and the association of gene modules to the phenotypic response. To link modules to phenotype, one approach is to find the correlated eigengenes of given modules or to integrate all eigengenes in regularized linear model. This manuscript describes the method from construction of co-expression network, module discovery, association between modules and phenotypic data, and finally to annotation/visualization.

Identification of GH17 gene family in Vitis vinifera and expression analysis of GH17 under various adversities

Glycoside hydrolase (GH, EC 3.2.1) is a group of enzymes that hydrolyzes glycosidic bonds and play a role in the hydrolysis and synthesis of sugars in living organisms. Vitis vinifera is an important fruit crop and it harbors GH17 gene family however, their function in grapes has not been systematically investigated. In this study, a total of 870 GH17 genes were identified from 14 plant species and their structural domain, sequence alignment, phylogenetic tree, collinear analysis, with the expression profiles of VviGH17 gene family was performed. The promoter analysis of VviGH17 gene showed the presence of cis-acting elements, which are responsive to plant growth and development. In addition, elements for plant hormones were found that are triggered in response to abiotic/biological stress. Transcriptomic data led to the identification of several VviGH17 genes, which are associated with bud dormancy and in response to abiotic stress. Transcript analysis was carried out for some of the selected VviGH17 genes RT-qPCR. VviGH17-16 and VviGH17-30 genes were differentially expressed during bud dormancy, fruit development and different abiotic stresses. Moreover, VviGH17-37 and VviGH17-44 were differentially expressed at fruit development, in response to abiotic stress. In addition, subcellular localization predicts that the VviGH17-16, VviGH17-30, and VviGH17-37 genes were located in the cell membrane, while VviGH17-44 gene was located in the vacuole. In conclusion, our study led to the identification of several GH17s and their probable role in development and stress.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01014-1.

Network aggregation improves gene function prediction of grapevine gene co-expression networks

Aggregation across multiple networks highlights robust co-expression interactions and improves the functional connectivity of grapevine gene co-expression networks. In recent years, the rapid accumulation of transcriptome datasets from diverse experimental conditions has enabled the widespread use of gene co-expression network (GCN) analysis in plants. In grapevine, GCN analysis has shown great promise for gene function prediction, however, measurable progress is currently lacking. Using accumulated microarray datasets from the grapevine whole-genome array (33 experiments, 1359 samples), we explored how meta-analysis through aggregation influences the functional connectivity (performance) of derived networks using guilt-by-association neighbor voting. Two annotation schemes, i.e. MapMan BIN and Pfam, at two sparsity thresholds, i.e. top 100 (stringent) and 300 (relaxed) ranked genes were evaluated. We observed that aggregating across multiple networks improves performance dramatically, with the aggregate outperforming the majority of functional terms across individual networks. Network sparsity and size (i.e. the number of samples and aggregates) were key factors influencing performance while the choice of annotation scheme had little. Systematic comparison with various state-of-the-art microarray and RNA-seq networks was also performed, however, none outperformed the aggregate microarray network despite having good predictive performance. Repeating these series of tests using a functional enrichment-based performance metric also showed remarkably consistent findings with guilt-by-association neighbor voting. To demonstrate its functionality, we explore the function and transcriptional regulation of grapevine EXPANSIN genes. We envisage that network aggregation will offer new and unique opportunities for gene function prediction in future grapevine functional genomics studies. To this end, we make the aggregate networks and associated metadata publicly available at VTC-Agg (https://sites.google.com/view/vtc-agg).

Recent advances in biotechnological studies on wild grapevines as valuable resistance sources for smart viticulture

Cultivated grapevines, Vitis vinifera subsp. sativa, are thought to have been domesticated from wild populations of Vitis vinifera subsp. sylvestris in Central Asia. V. vinifera subsp. sativa is one of the most economically important fruit crops worldwide. Since cultivated grapevines are susceptible to multiple biotic and abiotic soil factors, they also need to be grafted on resistant rootstocks that are mostly developed though hybridization between American wild grapevine species (V. berlandieri, V. riparia, and V. rupestris). Therefore, wild grapevine species are essential genetic materials for viticulture to face biotic and abiotic stresses in both cultivar and rootstock parts. Actually, viticulture faces several environmental constraints that are further intensified by climate change. Recently, several reports on biotic and abiotic stresses-response in wild grapevines revealed accessions tolerant to different constraints. The emergence of advanced techniques such as omics technologies, marker-assisted selection (MAS), and functional analysis tools allowed a more detailed characterization of resistance mechanisms in these wild grapevines and suggest a number of species (V. rotundifolia, V. rupestris, V. riparia, V. berlandieri and V. amurensis) have untapped potential for new resistance traits including disease resistance loci and key tolerance genes. The present review reports on the importance of different biotechnological tools in exploring and examining wild grapevines tolerance mechanisms that can be employed to promote elite cultivated grapevines under climate change conditions.

Gene network modules associated with abiotic stress response in tolerant rice genotypes identified by transcriptome meta-analysis

Abiotic stress tolerance is a complex trait regulated by multiple genes and gene networks in plants. A range of abiotic stresses are known to limit rice productivity. Meta-transcriptomics has emerged as a powerful approach to decipher stress-associated molecular network in model crops. However, retaining specificity of gene expression in tolerant and susceptible genotypes during meta-transcriptome analysis is important for understanding genotype-dependent stress tolerance mechanisms. Addressing this aspect, we describe here "abiotic stress tolerant" (ASTR) genes and networks specifically and differentially expressing in tolerant rice genotypes in response to different abiotic stress conditions. We identified 6,956 ASTR genes, key hub regulatory genes, transcription factors, and functional modules having significant association with abiotic stress-related ontologies and cis-motifs. Out of the 6956 ASTR genes, 73 were co-located within the boundary of previously identified abiotic stress trait-related quantitative trait loci. Functional annotation of 14 uncharacterized ASTR genes is proposed using multiple computational methods. Around 65% of the top ASTR genes were found to be differentially expressed in at least one of the tolerant genotypes under different stress conditions (cold, salt, drought, or heat) from publicly available RNAseq data comparison. The candidate ASTR genes specifically associated with tolerance could be utilized for engineering rice and possibly other crops for broad-spectrum tolerance to abiotic stresses.

Facing Climate Change: Biotechnology of Iconic Mediterranean Woody Crops

The Mediterranean basin is especially sensitive to the adverse outcomes of climate change and especially to variations in rainfall patterns and the incidence of extremely high temperatures. These two concurring adverse environmental conditions will surely have a detrimental effect on crop performance and productivity that will be particularly severe on woody crops such as citrus, olive and grapevine that define the backbone of traditional Mediterranean agriculture. These woody species have been traditionally selected for traits such as improved fruit yield and quality or alteration in harvesting periods, leaving out traits related to plant field performance. This is currently a crucial aspect due to the progressive and imminent effects of global climate change. Although complete genome sequence exists for sweet orange (Citrus sinensis) and clementine (Citrus clementina), olive tree (Olea europaea) and grapevine (Vitis vinifera), the development of biotechnological tools to improve stress tolerance still relies on the study of the available genetic resources including interspecific hybrids, naturally occurring (or induced) polyploids and wild relatives under field conditions. To this respect, post-genomic era studies including transcriptomics, metabolomics and proteomics provide a wide and unbiased view of plant physiology and biochemistry under adverse environmental conditions that, along with high-throughput phenotyping, could contribute to the characterization of plant genotypes exhibiting physiological and/or genetic traits that are correlated to abiotic stress tolerance. The ultimate goal of precision agriculture is to improve crop productivity, in terms of yield and quality, making a sustainable use of land and water resources under adverse environmental conditions using all available biotechnological tools and high-throughput phenotyping. This review focuses on the current state-of-the-art of biotechnological tools such as high throughput -omics and phenotyping on grapevine, citrus and olive and their contribution to plant breeding programs.

PbrMYB169 positively regulates lignification of stone cells in pear fruit

Stone cells negatively affect fruit quality because of their firm and lignified cell walls, so are targets for reduction in pear breeding programmes. However, there is only limited knowledge of the molecular mechanisms underlying the formation of stone cells. Here, we show that PbrMYB169, an R2R3 MYB transcription factor, of Pyrus bretschneideri positively regulates lignification of stone cells in pear fruit. PbrMYB169 was shown to be co-expressed with lignin biosynthesis genes during pear fruit development, and this co-expression pattern was coincident with stone cell formation in the fruit of Pyrus bretschneideri 'Dangshansuli'. The PbrMYB169 expression level was also positively correlated with stone cell content in 36 pear cultivars tested. PbrMYB169 protein significantly activated the promoter of lignin genes C3H1, CCR1, CCOMT2, CAD, 4CL1, 4CL2, HCT2, and LAC18 via binding to AC elements [ACC(T/A)ACC] in these promoters. Furthermore, overexpression of PbrMYB169 in transgenic Arabidopsis plants enhanced the expression of lignin genes, and increased lignin deposition and cell wall thickness of vessel elements, but did not change the ratio of syringyl and guaiacyl lignin monomers. In conclusion, PbrMYB169 appears to be a transcriptional activator of lignin biosynthesis and regulates secondary wall formation in fruit stone cells. This study advances the understanding of the regulation of lignin biosynthesis and provides valuable molecular genetic information for reducing stone cell content in pear fruit.

AppleMDO: A Multi-Dimensional Omics Database for Apple Co-Expression Networks and Chromatin States

As an economically important crop, apple is one of the most cultivated fruit trees in temperate regions worldwide. Recently, a large number of high-quality transcriptomic and epigenomic datasets for apple were made available to the public, which could be helpful in inferring gene regulatory relationships and thus predicting gene function at the genome level. Through integration of the available apple genomic, transcriptomic, and epigenomic datasets, we constructed co-expression networks, identified functional modules, and predicted chromatin states. A total of 112 RNA-seq datasets were integrated to construct a global network and a conditional network (tissue-preferential network). Furthermore, a total of 1,076 functional modules with closely related gene sets were identified to assess the modularity of biological networks and further subjected to functional enrichment analysis. The results showed that the function of many modules was related to development, secondary metabolism, hormone response, and transcriptional regulation. Transcriptional regulation is closely related to epigenetic marks on chromatin. A total of 20 epigenomic datasets, which included ChIP-seq, DNase-seq, and DNA methylation analysis datasets, were integrated and used to classify chromatin states. Based on the ChromHMM algorithm, the genome was divided into 620,122 fragments, which were classified into 24 states according to the combination of epigenetic marks and enriched-feature regions. Finally, through the collaborative analysis of different omics datasets, the online database AppleMDO (http://bioinformatics.cau.edu.cn/AppleMDO/) was established for cross-referencing and the exploration of possible novel functions of apple genes. In addition, gene annotation information and functional support toolkits were also provided. Our database might be convenient for researchers to develop insights into the function of genes related to important agronomic traits and might serve as a reference for other fruit trees.

miRVIT: A Novel miRNA Database and Its Application to Uncover Vitis Responses to Flavescence dorée Infection

Micro(mi)RNAs play crucial roles in plant developmental processes and in defense responses to biotic and abiotic stresses. In the last years, many works on small RNAs in grapevine (Vitis spp.) were published, and several conserved and putative novel grapevine-specific miRNAs were identified. In order to reorganize the high quantity of available data, we produced "miRVIT," the first database of all novel grapevine miRNA candidates characterized so far, and still not deposited in miRBase. To this aim, each miRNA accession was renamed, repositioned in the last version of the grapevine genome, and compared with all the novel and conserved miRNAs detected in grapevine. Conserved and novel miRNAs cataloged in miRVIT were then used for analyzing Vitis vinifera plants infected by Flavescence dorée (FD), one of the most severe phytoplasma diseases affecting grapevine. The analysis of small RNAs from healthy, recovered (plants showing spontaneous and stable remission of symptoms), and FD-infected "Barbera" grapevines showed that FD altered the expression profiles of several miRNAs, including those involved in cell development and photosynthesis, jasmonate signaling, and disease resistance response. The application of miRVIT in a biological context confirmed the effectiveness of the followed approach, especially for the identification of novel miRNA candidates in grapevine. miRVIT database is available at http://mirvit.ipsp.cnr.it. Highlights: The application of the newly produced database of grapevine novel miRNAs to the analysis of plants infected by Flavescence dorée reveals key roles of miRNAs in photosynthesis and jasmonate signaling.

Combinatorial Regulation of Stilbene Synthase Genes by WRKY and MYB Transcription Factors in Grapevine (Vitis vinifera L.)

Stilbene synthase (STS) is the key enzyme leading to the biosynthesis of resveratrol. Recently we reported two R2R3-MYB transcription factor (TF) genes that regulate the stilbene biosynthetic pathway in grapevine: VviMYB14 and VviMYB15. These genes are strongly co-expressed with STS genes under a range of stress and developmental conditions, in agreement with the specific activation of STS promoters by these TFs. Genome-wide gene co-expression analysis using two separate transcriptome compendia based on microarray and RNA sequencing data revealed that WRKY TFs were the top TF family correlated with STS genes. On the basis of correlation frequency, four WRKY genes, namely VviWRKY03, VviWRKY24, VviWRKY43 and VviWRKY53, were further shortlisted and functionally validated. Expression analyses under both unstressed and stressed conditions, together with promoter-luciferase reporter assays, suggested different hierarchies for these TFs in the regulation of the stilbene biosynthetic pathway. In particular, VviWRKY24 seems to act as a singular effector in the activation of the VviSTS29 promoter, while VviWRKY03 acts through a combinatorial effect with VviMYB14, suggesting that these two regulators may interact at the protein level as previously reported in other species.

Co-expression network analysis and cis-regulatory element enrichment determine putative functions and regulatory mechanisms of grapevine ATL E3 ubiquitin ligases

Arabidopsis thaliana Toxicos en Levadura (ATL) proteins are a subclass of the RING-H2 zinc finger binding E3 ubiquitin ligases. The grapevine (Vitis vinifera) ATL family was recently characterized, revealing 96 members that are likely to be involved in several physiological processes through protein ubiquitination. However, the final targets and biological functions of most ATL E3 ligases are still unknown. We analyzed the co-expression networks among grapevine ATL genes across a set of transcriptomic data related to defense and abiotic stress, combined with a condition-independent dataset. This revealed strong correlations between ATL proteins and diverse signal transduction components and transcriptional regulators, in particular those involved in immunity. An enrichment analysis of cis-regulatory elements in ATL gene promoters and related co-expressed genes highlighted the importance of hormones in the regulation of ATL gene expression. Our work identified several ATL proteins as candidates for further studies aiming to decipher specific grapevine resistance mechanisms activated in response to pathogens.

Discovering Causal Relationships in Grapevine Expression Data to Expand Gene Networks. A Case Study: Four Networks Related to Climate Change

In recent years the scientific community has been heavily engaged in studying the grapevine response to climate change. Final goal is the identification of key genetic traits to be used in grapevine breeding and the setting of agronomic practices to improve climatic resilience. The increasing availability of transcriptomic studies, describing gene expression in many tissues and developmental, or treatment conditions, have allowed the implementation of gene expression compendia, which enclose a huge amount of information. The mining of transcriptomic data represents an effective approach to expand a known local gene network (LGN) by finding new related genes. We recently published a pipeline based on the iterative application of the PC-algorithm, named NES²RA, to expand gene networks in Escherichia coli and Arabidopsis thaliana. Here, we propose the application of this method to the grapevine transcriptomic compendium Vespucci, in order to expand four LGNs related to the grapevine response to climate change. Two networks are related to the secondary metabolic pathways for anthocyanin and stilbenoid synthesis, involved in the response to solar radiation, whereas the other two are signaling networks, related to the hormones abscisic acid and ethylene, possibly involved in the regulation of cell water balance and cuticle transpiration. The expansion networks produced by NES²RA algorithm have been evaluated by comparison with experimental data and biological knowledge on the identified genes showing fairly good consistency of the results. In addition, the algorithm was effective in retaining only the most significant interactions among the genes providing a useful framework for experimental validation. The application of the NES²RA to Vitis vinifera expression data by means of the BOINC-based implementation is available upon request (valter.cavecchia@cnr.it).

Potassium in the Grape (Vitis vinifera L.) Berry: Transport and Function

K+ is the most abundant cation in the grape berry. Here we focus on the most recent information in the long distance transport and partitioning of K+ within the grapevine and postulate on the potential role of K+ in berry sugar accumulation, berry water relations, cellular growth, disease resistance, abiotic stress tolerance and mitigating senescence. By integrating information from several different plant systems we have been able to generate new hypotheses on the integral functions of this predominant cation and to improve our understanding of how these functions contribute to grape berry growth and ripening. Valuable contributions to the study of K+ in membrane stabilization, turgor maintenance and phloem transport have allowed us to propose a mechanistic model for the role of this cation in grape berry development.

Genome-wide analysis of cis-regulatory element structure and discovery of motif-driven gene co-expression networks in grapevine

Coordinated transcriptional and metabolic reprogramming ensures a plant’s continued growth and survival under adverse environmental conditions. Transcription factors (TFs) act to modulate gene expression through complex cis-regulatory element (CRE) interactions. Genome-wide analysis of known plant CREs was performed for all currently predicted protein-coding gene promoters in grapevine (Vitis vinifera L.). Many CREs such as abscisic acid (ABA)-responsive, drought-responsive, auxin-responsive, and evening elements, exhibit bona fide CRE properties such as strong position bias towards the transcription start site (TSS) and over-representation when compared with random promoters. Genes containing these CREs are enriched in a large repertoire of plant biological pathways. Large-scale transcriptome analyses also show that these CREs are highly implicated in grapevine development and stress response. Numerous CRE-driven modules in condition-specific gene co-expression networks (GCNs) were identified and many of these modules were highly enriched for plant biological functions. Several modules corroborate known roles of CREs in drought response, pathogen defense, cell wall metabolism, and fruit ripening, whereas others reveal novel functions in plants. Comparisons with Arabidopsis suggest a general conservation in promoter architecture, gene expression dynamics, and GCN structure across species. Systems analyses of CREs provide insights into the grapevine cis-regulatory code and establish a foundation for future genomic studies in grapevine.

Connectivity in gene coexpression networks negatively correlates with rates of molecular evolution in flowering plants

Gene coexpression networks are a useful tool for summarizing transcriptomic data and providing insight into patterns of gene regulation in a variety of species. Though there has been considerable interest in studying the evolution of network topology across species, less attention has been paid to the relationship between network position and patterns of molecular evolution. Here, we generated coexpression networks from publicly available expression data for seven flowering plant taxa (Arabidopsis thaliana, Glycine max, Oryza sativa, Populus spp., Solanum lycopersicum, Vitis spp., and Zea mays) to investigate the relationship between network position and rates of molecular evolution. We found a significant negative correlation between network connectivity and rates of molecular evolution, with more highly connected (i.e., “hub”) genes having significantly lower nonsynonymous substitution rates and dN/dS ratios compared to less highly connected (i.e., “peripheral”) genes across the taxa surveyed. These findings suggest that more centrally located hub genes are, on average, subject to higher levels of evolutionary constraint than are genes located on the periphery of gene coexpression networks. The consistency of this result across disparate taxa suggests that it holds for flowering plants in general, as opposed to being a species-specific phenomenon.

Co-expressed Pathways DataBase for Tomato: a database to predict pathways relevant to a query gene

BACKGROUND

Gene co-expression, the similarity of gene expression profiles under various experimental conditions, has been used as an indicator of functional relationships between genes, and many co-expression databases have been developed for predicting gene functions. These databases usually provide users with a co-expression network and a list of strongly co-expressed genes for a query gene. Several of these databases also provide functional information on a set of strongly co-expressed genes (i.e., provide biological processes and pathways that are enriched in these strongly co-expressed genes), which is generally analyzed via over-representation analysis (ORA). A limitation of this approach may be that users can predict gene functions only based on the strongly co-expressed genes.

RESULTS

In this study, we developed a new co-expression database that enables users to predict the function of tomato genes from the results of functional enrichment analyses of co-expressed genes while considering the genes that are not strongly co-expressed. To achieve this, we used the ORA approach with several thresholds to select co-expressed genes, and performed gene set enrichment analysis (GSEA) applied to a ranked list of genes ordered by the co-expression degree. We found that internal correlation in pathways affected the significance levels of the enrichment analyses. Therefore, we introduced a new measure for evaluating the relationship between the gene and pathway, termed the percentile (p)-score, which enables users to predict functionally relevant pathways without being affected by the internal correlation in pathways. In addition, we evaluated our approaches using receiver operating characteristic curves, which concluded that the p-score could improve the performance of the ORA.

CONCLUSIONS

We developed a new database, named Co-expressed Pathways DataBase for Tomato, which is available at http://cox-path-db.kazusa.or.jp/tomato . The database allows users to predict pathways that are relevant to a query gene, which would help to infer gene functions.

Constructing Integrated Networks for Identifying New Secondary Metabolic Pathway Regulators in Grapevine: Recent Applications and Future Opportunities

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.

A Concise Review on Multi-Omics Data Integration for Terroir Analysis in Vitis vinifera

Vitis vinifera (grapevine) is one of the most important fruit crops, both for fresh consumption and wine and spirit production. The term terroir is frequently used in viticulture and the wine industry to relate wine sensory attributes to its geographic origin. Although, it can be cultivated in a wide range of environments, differences in growing conditions have a significant impact on fruit traits that ultimately affect wine quality. Understanding how fruit quality and yield are controlled at a molecular level in grapevine in response to environmental cues has been a major driver of research. Advances in the area of genomics, epigenomics, transcriptomics, proteomics and metabolomics, have significantly increased our knowledge on the abiotic regulation of yield and quality in many crop species, including V. vinifera. The integrated analysis of multiple 'omics' can give us the opportunity to better understand how plants modulate their response to different environments. However, 'omics' technologies provide a large amount of biological data and its interpretation is not always straightforward, especially when different 'omic' results are combined. Here we examine the current strategies used to integrate multi-omics, and how these have been used in V. vinifera. In addition, we also discuss the importance of including epigenomics data when integrating omics data as epigenetic mechanisms could play a major role as an intermediary between the environment and the genome.

Omics Approaches for Understanding Grapevine Berry Development: Regulatory Networks Associated with Endogenous Processes and Environmental Responses

Grapevine fruit development is a dynamic process that can be divided into three stages: formation (I), lag (II), and ripening (III), in which physiological and biochemical changes occur, leading to cell differentiation and accumulation of different solutes. These stages can be positively or negatively affected by multiple environmental factors. During the last decade, efforts have been made to understand berry development from a global perspective. Special attention has been paid to transcriptional and metabolic networks associated with the control of grape berry development, and how external factors affect the ripening process. In this review, we focus on the integration of global approaches, including proteomics, metabolomics, and especially transcriptomics, to understand grape berry development. Several aspects will be considered, including seed development and the production of seedless fruits; veraison, at which anthocyanin accumulation begins in the berry skin of colored varieties; and hormonal regulation of berry development and signaling throughout ripening, focusing on the transcriptional regulation of hormone receptors, protein kinases, and genes related to secondary messenger sensing. Finally, berry responses to different environmental factors, including abiotic (temperature, water-related stress and UV-B radiation) and biotic (fungi and viruses) stresses, and how they can significantly modify both, development and composition of vine fruit, will be discussed. Until now, advances have been made due to the application of Omics tools at different molecular levels. However, the potential of these technologies should not be limited to the study of single-level questions; instead, data obtained by these platforms should be integrated to unravel the molecular aspects of grapevine development. Therefore, the current challenge is the generation of new tools that integrate large-scale data to assess new questions in this field, and to support agronomical practices.

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

Grapevine (Vitis vinifera L.) is a species well known for its adaptation to radiation. However, photomorphogenic factors related to UV-B responses have not been molecularly characterized. We cloned and studied the role of UV-B RECEPTOR (UVR1), ELONGATED HYPOCOTYL 5 (HY5), and HY5 HOMOLOGUE (HYH) from V. vinifera We performed gene functional characterizations, generated co-expression networks, and tested them in different environmental conditions. These genes complemented the Arabidopsis uvr8 and hy5 mutants in morphological and secondary metabolic responses to radiation. We combined microarray and RNA sequencing (RNA-seq) data with promoter inspections to identify HY5 and HYH putative target genes and their DNA binding preferences. Despite sharing a large set of common co-expressed genes, we found different hierarchies for HY5 and HYH depending on the organ and stress condition, reflecting both co-operative and partially redundant roles. New candidate UV-B gene markers were supported by the presence of HY5-binding sites. These included a set of flavonol-related genes that were up-regulated in a HY5 transient expression assay. We irradiated in vitro plantlets and fruits from old potted vines with high and low UV-B exposures and followed the accumulation of flavonols and changes in gene expression in comparison with non-irradiated conditions. UVR1, HY5, and HYH expression varied with organ, developmental stage, and type of radiation. Surprisingly, UVR1 expression was modulated by shading and temperature in berries, but not by UV-B radiation. We propose that the UV-B response machinery favours berry flavonol accumulation through the activation of HY5 and HYH at different developmental stages at both high and low UV-B exposures.

Genome-wide identification and characterization of the NF-Y gene family in grape (vitis vinifera L.)

Background

Nuclear factor Y (NF-Y) transcription factor is composed of three distinct subunits: NF-YA, NF-YB and NF-YC. Many members of NF-Y family have been reported to be key regulators in plant development, phytohormone signaling and drought tolerance. However, the function of the NF-Y family is less known in grape (Vitis vinifera L.).

Results

A total of 34 grape NF-Y genes that distributed unevenly on grape (V. vinifera) chromosomes were identified in this study. Phylogenetic analysis was performed to predict functional similarities between Arabidopsis thaliana and grape NF-Y genes. Comparison of the structures of grape NF-Y genes (VvNF-Ys) revealed their functional conservation and alteration. Furthermore, we investigated the expression profiles of VvNF-Ys in response to various stresses, phytohormone treatments, and in leaves and grape berries with various sugar contents at different developmental stages. The relationship between VvNF-Y transcript levels and sugar content was examined to select candidates for exogenous sugar treatments. Quantitative real-time PCR (qPCR) indicated that many VvNF-Ys responded to different sugar stimuli with variations in transcript abundance. qPCR and publicly available microarray data suggest that VvNF-Ys exhibit distinct expression patterns in different grape organs and developmental stages, and a number of VvNF-Ys may participate in responses to multiple abiotic and biotic stresses, phytohormone treatments and sugar accumulation or metabolism.

Conclusions

In this study, we characterized 34 VvNF-Ys based on their distributions on chromosomes, gene structures, phylogenetic relationship with Arabidopsis NF-Y genes, and their expression patterns. The potential roles of VvNF-Ys in sugar accumulation or metabolism were also investigated. Altogether, the data provide significant insights on VvNF-Ys, and lay foundations for further functional studies of NF-Y genes in grape.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2989-3) contains supplementary material, which is available to authorized users.

A systems-oriented analysis of the grapevine R2R3-MYB transcription factor family uncovers new insights into the regulation of stilbene accumulation

R2R3-MYB transcription factors (TFs) belong to a large and functionally diverse protein superfamily in plants. In this study, we explore the evolution and function of this family in grapevine (Vitis vinifera L.), a high-value fruit crop. We identified and manually curated 134 genes using RNA-Seq data, and named them systematically according to the Super-Nomenclature Committee. We identified novel genes, splicing variants and grapevine/woody-specific duplicated subgroups, suggesting possible neo- and sub-functionalization events. Regulatory network analysis ascribed biological functions to uncharacterized genes and validated those of known genes (e.g. secondary cell wall biogenesis and flavonoid biosynthesis). A comprehensive analysis of different MYB binding motifs in the promoters of co-expressed genes predicted grape R2R3-MYB binding preferences and supported evidence for putative downstream targets. Enrichment of cis-regulatory motifs for diverse TFs reinforced the notion of transcriptional coordination and interaction between MYBs and other regulators. Analysis of the network of Subgroup 2 showed that the resveratrol-related VviMYB14 and VviMYB15 share common co-expressed STILBENE SYNTHASE genes with the uncharacterized VviMYB13. These regulators have distinct expression patterns within organs and in response to biotic and abiotic stresses, suggesting a pivotal role of VviMYB13 in regulating stilbene accumulation in vegetative tissues and under biotic stress conditions.

Combined physiological, transcriptome, and cis-regulatory element analyses indicate that key aspects of ripening, metabolism, and transcriptional program in grapes (Vitis vinifera L.) are differentially modulated accordingly to fruit size

BACKGROUND

In wine grape production, management practices have been adopted to optimize grape and wine quality attributes by producing, or screening for, berries of smaller size. Fruit size and composition are influenced by numerous factors that include both internal (e.g. berry hormone metabolism) and external (e.g. environment and cultural practices) factors. Combined physiological, biochemical, and transcriptome analyses were performed to improve our current understanding of metabolic and transcriptional pathways related to berry ripening and composition in berries of different sizes.

RESULTS

The comparison of berry physiology between small and large berries throughout development (from 31 to 121 days after anthesis, DAA) revealed significant differences in firmness, the rate of softening, and sugar accumulation at specific developmental stages. Small berries had significantly higher skin to berry weight ratio, lower number of seeds per berry, and higher anthocyanin concentration compared to large berries. RNA-sequencing analyses of berry skins at 47, 74, 103, and 121 DAA revealed a total of 3482 differentially expressed genes between small and large berries. Abscisic acid, auxin, and ethylene hormone pathway genes were differentially modulated between berry sizes. Fatty acid degradation and stilbenoid pathway genes were upregulated at 47 DAA while cell wall degrading and modification genes were downregulated at 74 DAA in small compared to large berries. In the late ripening stage, concerted upregulation of the general phenylpropanoid and stilbenoid pathway genes and downregulation of flavonoid pathway genes were observed in skins of small compared to large berries. Cis-regulatory element analysis of differentially expressed hormone, fruit texture, flavor, and aroma genes revealed an enrichment of specific regulatory motifs related to bZIP, bHLH, AP2/ERF, NAC, MYB, and MADS-box transcription factors.

CONCLUSIONS

The study demonstrates that physiological and compositional differences between berries of different sizes parallel transcriptome changes that involve fruit texture, flavor, and aroma pathways. These results suggest that, in addition to direct effects brought about by differences in size, key aspects involved in the regulation of ripening likely contribute to different quality profiles between small and large berries.

Transcriptome and metabolite profiling reveals that prolonged drought modulates the phenylpropanoid and terpenoid pathway in white grapes (Vitis vinifera L.)

BACKGROUND

Secondary metabolism contributes to the adaptation of a plant to its environment. In wine grapes, fruit secondary metabolism largely determines wine quality. Climate change is predicted to exacerbate drought events in several viticultural areas, potentially affecting the wine quality. In red grapes, water deficit modulates flavonoid accumulation, leading to major quantitative and compositional changes in the profile of the anthocyanin pigments; in white grapes, the effect of water deficit on secondary metabolism is still largely unknown.

RESULTS

In this study we investigated the impact of water deficit on the secondary metabolism of white grapes using a large scale metabolite and transcript profiling approach in a season characterized by prolonged drought. Irrigated grapevines were compared to non-irrigated grapevines that suffered from water deficit from early stages of berry development to harvest. A large effect of water deficit on fruit secondary metabolism was observed. Increased concentrations of phenylpropanoids, monoterpenes, and tocopherols were detected, while carotenoid and flavonoid accumulations were differentially modulated by water deficit according to the berry developmental stage. The RNA-sequencing analysis carried out on berries collected at three developmental stages-before, at the onset, and at late ripening-indicated that water deficit affected the expression of 4,889 genes. The Gene Ontology category secondary metabolic process was overrepresented within up-regulated genes at all the stages of fruit development considered, and within down-regulated genes before ripening. Eighteen phenylpropanoid, 16 flavonoid, 9 carotenoid, and 16 terpenoid structural genes were modulated by water deficit, indicating the transcriptional regulation of these metabolic pathways in fruit exposed to water deficit. An integrated network and promoter analyses identified a transcriptional regulatory module that encompasses terpenoid genes, transcription factors, and enriched drought-responsive elements in the promoter regions of those genes as part of the grapes response to drought.

CONCLUSION

Our study reveals that grapevine berries respond to drought by modulating several secondary metabolic pathways, and particularly, by stimulating the production of phenylpropanoids, the carotenoid zeaxanthin, and of volatile organic compounds such as monoterpenes, with potential effects on grape and wine antioxidant potential, composition, and sensory features.

Identification of Gene Modules Associated with Low Temperatures Response in Bambara Groundnut by Network-Based Analysis

Bambara groundnut (Vigna subterranea (L.) Verdc.) is an African legume and is a promising underutilized crop with good seed nutritional values. Low temperature stress in a number of African countries at night, such as Botswana, can effect the growth and development of bambara groundnut, leading to losses in potential crop yield. Therefore, in this study we developed a computational pipeline to identify and analyze the genes and gene modules associated with low temperature stress responses in bambara groundnut using the cross-species microarray technique (as bambara groundnut has no microarray chip) coupled with network-based analysis. Analyses of the bambara groundnut transcriptome using cross-species gene expression data resulted in the identification of 375 and 659 differentially expressed genes (p<0.01) under the sub-optimal (23°C) and very sub-optimal (18°C) temperatures, respectively, of which 110 genes are commonly shared between the two stress conditions. The construction of a Highest Reciprocal Rank-based gene co-expression network, followed by its partition using a Heuristic Cluster Chiseling Algorithm resulted in 6 and 7 gene modules in sub-optimal and very sub-optimal temperature stresses being identified, respectively. Modules of sub-optimal temperature stress are principally enriched with carbohydrate and lipid metabolic processes, while most of the modules of very sub-optimal temperature stress are significantly enriched with responses to stimuli and various metabolic processes. Several transcription factors (from MYB, NAC, WRKY, WHIRLY & GATA classes) that may regulate the downstream genes involved in response to stimulus in order for the plant to withstand very sub-optimal temperature stress were highlighted. The identified gene modules could be useful in breeding for low-temperature stress tolerant bambara groundnut varieties.

A Grapevine TTG2-Like WRKY Transcription Factor Is Involved in Regulating Vacuolar Transport and Flavonoid Biosynthesis

A small set of TTG2-like homolog proteins from different species belonging to the WRKY family of transcription factors were shown to share a similar mechanism of action and to control partially conserved biochemical/developmental processes in their native species. In particular, by activating P-ATPases residing on the tonoplast, PH3 from Petunia hybrida promotes vacuolar acidification in petal epidermal cells whereas TTG2 from Arabidopsis thaliana enables the accumulation of proanthocyanidins in the seed coat. In this work we functionally characterized VvWRKY26 identified as the closest grapevine homolog of PhPH3 and AtTTG2. When constitutively expressed in petunia ph3 mutant, VvWRKY26 can fulfill the PH3 function in the regulation of vacuolar pH and restores the wild type pigmentation phenotype. By a global correlation analysis of gene expression and by transient over-expression in Vitis vinifera, we showed transcriptomic relationships of VvWRKY26 with many genes related to vacuolar acidification and transport in grapevine. Moreover, our results indicate an involvement in flavonoid pathway possibly restricted to the control of proanthocyanidin biosynthesis that is consistent with its expression pattern in grape berry tissues. Overall, the results show that, in addition to regulative mechanisms and biological roles shared with TTG2-like orthologs, VvWRKY26 can play roles in fleshy fruit development that have not been previously reported in studies from dry fruit species. This study paves the way toward the comprehension of the regulatory network controlling vacuolar acidification and flavonoid accumulation mechanisms that contribute to the final berry quality traits in grapevine.

VESPUCCI: Exploring Patterns of Gene Expression in Grapevine

Large-scale transcriptional studies aim to decipher the dynamic cellular responses to a stimulus, like different environmental conditions. In the era of high-throughput omics biology, the most used technologies for these purposes are microarray and RNA-Seq, whose data are usually required to be deposited in public repositories upon publication. Such repositories have the enormous potential to provide a comprehensive view of how different experimental conditions lead to expression changes, by comparing gene expression across all possible measured conditions. Unfortunately, this task is greatly impaired by differences among experimental platforms that make direct comparisons difficult. In this paper, we present the Vitis Expression Studies Platform Using COLOMBOS Compendia Instances (VESPUCCI), a gene expression compendium for grapevine which was built by adapting an approach originally developed for bacteria, and show how it can be used to investigate complex gene expression patterns. We integrated nearly all publicly available microarray and RNA-Seq expression data: 1608 gene expression samples from 10 different technological platforms. Each sample has been manually annotated using a controlled vocabulary developed ad hoc to ensure both human readability and computational tractability. Expression data in the compendium can be visually explored using several tools provided by the web interface or can be programmatically accessed using the REST interface. VESPUCCI is freely accessible at http://vespucci.colombos.fmach.it.

Regulation of flavonol content and composition in (Syrah×Pinot Noir) mature grapes: integration of transcriptional profiling and metabolic quantitative trait locus analyses

Flavonols are a ubiquitous class of flavonoids that accumulate preferentially in flowers and mature berries. Besides their photo-protective function, they play a fundamental role during winemaking, stabilizing the colour by co-pigmentation with anthocyanins and contributing to organoleptic characteristics. Although the general flavonol pathway has been genetically and biochemically elucidated, the genetic control of flavonol content and composition at harvest is still not clear. To this purpose, the grapes of 170 segregating F1 individuals from a 'Syrah'×'Pinot Noir' population were evaluated at the mature stage for the content of six flavonol aglycons in four seasons. Metabolic data in combination with genetic data enabled the identification of 16 mQTLs (metabolic quantitative trait loci). For the first time, major genetic control by the linkage group 2 (LG 2)/MYBA region on flavonol variation, in particular of tri-hydroxylated flavonols, is demonstrated. Moreover, seven regions specifically associated with the fine control of flavonol biosynthesis are identified. Gene expression profiling of two groups of individuals significantly divergent for their skin flavonol content identified a large set of differentially modulated transcripts. Among these, the transcripts coding for MYB and bZIP transcription factors, methyltranferases, and glucosyltranferases specific for flavonols, proteins, and factors belonging to the UV-B signalling pathway and co-localizing with the QTL regions are proposed as candidate genes for the fine regulation of flavonol content and composition in mature grapes.

New candidate genes for the fine regulation of the colour of grapes

In the last decade, great progress has been made in clarifying the main determinants of anthocyanin accumulation in grape berry skin. However, the molecular details of the fine variation among cultivars, which ultimately contributes to wine typicity, are still not completely understood. To shed light on this issue, the grapes of 170 F1 progeny from the cross 'Syrah'×'Pinot Noir' were characterized at the mature stage for the content of 15 anthocyanins during four growing seasons. This huge data set was used in combination with a dense genetic map to detect genomic regions controlling the anthocyanin pathway both at key enzymatic points and at particular branches. Genes putatively involved in fine tuning the global regulation of anthocyanin biosynthesis were identified by exploring the gene predictions in the QTL (quantitative trait locus) confidence intervals and their expression profile during berry development in offspring with contrasting anthocyanin accumulation. New information on some aspects which had scarcely been investigated so far, such as anthocyanin transport into the vacuole, or completely neglected, such as acylation, is provided. These genes represent a valuable resource in grapevine molecular-based breeding programmes to improve both fruit and wine quality and to tailor wine sensory properties according to consumer demand.

A knowledge base for Vitis vinifera functional analysis

Background

Vitis vinifera (Grapevine) is the most important fruit species in the modern world. Wine and table grapes sales contribute significantly to the economy of major wine producing countries. The most relevant goals in wine production concern quality and safety. In order to significantly improve the achievement of these objectives and to gain biological knowledge about cultivars, a genomic approach is the most reliable strategy. The recent grapevine genome sequencing offers the opportunity to study the potential roles of genes and microRNAs in fruit maturation and other physiological and pathological processes. Although several systems allowing the analysis of plant genomes have been reported, none of them has been designed specifically for the functional analysis of grapevine genomes of cultivars under environmental stress in connection with microRNA data.

Description

Here we introduce a novel knowledge base, called BIOWINE, designed for the functional analysis of Vitis vinifera genomes of cultivars present in Sicily. The system allows the analysis of RNA-seq experiments of two different cultivars, namely Nero d'Avola and Nerello Mascalese. Samples were taken under different climatic conditions of phenological phases, diseases, and geographic locations. The BIOWINE web interface is equipped with data analysis modules for grapevine genomes. In particular users may analyze the current genome assembly together with the RNA-seq data through a customized version of GBrowse. The web interface allows users to perform gene set enrichment by exploiting third-party databases.

Conclusions

BIOWINE is a knowledge base implementing a set of bioinformatics tools for the analysis of grapevine genomes. The system aims to increase our understanding of the grapevine varieties and species of Sicilian products focusing on adaptability to different climatic conditions, phenological phases, diseases, and geographic locations.

New insights into the evolutionary history of plant sorbitol dehydrogenase

BACKGROUND

Sorbitol dehydrogenase (SDH, EC 1.1.1.14) is the key enzyme involved in sorbitol metabolism in higher plants. SDH genes in some Rosaceae species could be divided into two groups. L-idonate-5-dehydrogenase (LIDH, EC 1.1.1.264) is involved in tartaric acid (TA) synthesis in Vitis vinifera and is highly homologous to plant SDHs. Despite efforts to understand the biological functions of plant SDH, the evolutionary history of plant SDH genes and their phylogenetic relationship with the V. vinifera LIDH gene have not been characterized.

RESULTS

A total of 92 SDH genes were identified from 42 angiosperm species. SDH genes have been highly duplicated within the Rosaceae family while monocot, Brassicaceae and most Asterid species exhibit singleton SDH genes. Core Eudicot SDHs have diverged into two phylogenetic lineages, now classified as SDH Class I and SDH Class II. V. vinifera LIDH was identified as a Class II SDH. Tandem duplication played a dominant role in the expansion of plant SDH family and Class II SDH genes were positioned in tandem with Class I SDH genes in several plant genomes. Protein modelling analyses of V. vinifera SDHs revealed 19 putative active site residues, three of which exhibited amino acid substitutions between Class I and Class II SDHs and were influenced by positive natural selection in the SDH Class II lineage. Gene expression analyses also demonstrated a clear transcriptional divergence between Class I and Class II SDH genes in V. vinifera and Citrus sinensis (orange).

CONCLUSIONS

Phylogenetic, natural selection and synteny analyses provided strong support for the emergence of SDH Class II by positive natural selection after tandem duplication in the common ancestor of core Eudicot plants. The substitutions of three putative active site residues might be responsible for the unique enzyme activity of V. vinifera LIDH, which belongs to SDH Class II and represents a novel function of SDH in V. vinifera that may be true also of other Class II SDHs. Gene expression analyses also supported the divergence of SDH Class II at the expression level. This study will facilitate future research into understanding the biological functions of plant SDHs.

Network-assisted crop systems genetics: network inference and integrative analysis

Although next-generation sequencing (NGS) technology has enabled the decoding of many crop species genomes, most of the underlying genetic components for economically important crop traits remain to be determined. Network approaches have proven useful for the study of the reference plant, Arabidopsis thaliana, and the success of network-based crop genetics will also require the availability of a genome-scale functional networks for crop species. In this review, we discuss how to construct functional networks and elucidate the holistic view of a crop system. The crop gene network then can be used for gene prioritization and the analysis of resequencing-based genome-wide association study (GWAS) data, the amount of which will rapidly grow in the field of crop science in the coming years.

A co-expression gene network associated with developmental regulation of apple fruit acidity

Apple fruit acidity, which affects the fruit's overall taste and flavor to a large extent, is primarily determined by the concentration of malic acid. Previous studies demonstrated that the major QTL malic acid (Ma) on chromosome 16 is largely responsible for fruit acidity variations in apple. Recent advances suggested that a natural mutation that gives rise to a premature stop codon in one of the two aluminum-activated malate transporter (ALMT)-like genes (called Ma1) is the genetic causal element underlying Ma. However, the natural mutation does not explain the developmental changes of fruit malate levels in a given genotype. Using RNA-seq data from the fruit of 'Golden Delicious' taken at 14 developmental stages from 1 week after full-bloom (WAF01) to harvest (WAF20), we characterized their transcriptomes in groups of high (12.2 ± 1.6 mg/g fw, WAF03-WAF08), mid (7.4 ± 0.5 mg/g fw, WAF01-WAF02 and WAF10-WAF14) and low (5.4 ± 0.4 mg/g fw, WAF16-WAF20) malate concentrations. Detailed analyses showed that a set of 3,066 genes (including Ma1) were expressed not only differentially (P FDR < 0.05) between the high and low malate groups (or between the early and late developmental stages) but also in significant (P < 0.05) correlation with malate concentrations. The 3,066 genes fell in 648 MapMan (sub-) bins or functional classes, and 19 of them were significantly (P FDR < 0.05) co-enriched or co-suppressed in a malate dependent manner. Network inferring using the 363 genes encompassed in the 19 (sub-) bins, identified a major co-expression network of 239 genes. Since the 239 genes were also differentially expressed between the early (WAF03-WAF08) and late (WAF16-WAF20) developmental stages, the major network was considered to be associated with developmental regulation of apple fruit acidity in 'Golden Delicious'.

Transcriptional Regulatory Network Analysis of MYB Transcription Factor Family Genes in Rice

MYB transcription factor (TF) is one of the largest TF families and regulates defense responses to various stresses, hormone signaling as well as many metabolic and developmental processes in plants. Understanding these regulatory hierarchies of gene expression networks in response to developmental and environmental cues is a major challenge due to the complex interactions between the genetic elements. Correlation analyses are useful to unravel co-regulated gene pairs governing biological process as well as identification of new candidate hub genes in response to these complex processes. High throughput expression profiling data are highly useful for construction of co-expression networks. In the present study, we utilized transcriptome data for comprehensive regulatory network studies of MYB TFs by "top-down" and "guide-gene" approaches. More than 50% of OsMYBs were strongly correlated under 50 experimental conditions with 51 hub genes via "top-down" approach. Further, clusters were identified using Markov Clustering (MCL). To maximize the clustering performance, parameter evaluation of the MCL inflation score (I) was performed in terms of enriched GO categories by measuring F-score. Comparison of co-expressed cluster and clads analyzed from phylogenetic analysis signifies their evolutionarily conserved co-regulatory role. We utilized compendium of known interaction and biological role with Gene Ontology enrichment analysis to hypothesize function of coexpressed OsMYBs. In the other part, the transcriptional regulatory network analysis by "guide-gene" approach revealed 40 putative targets of 26 OsMYB TF hubs with high correlation value utilizing 815 microarray data. The putative targets with MYB-binding cis-elements enrichment in their promoter region, functional co-occurrence as well as nuclear localization supports our finding. Specially, enrichment of MYB binding regions involved in drought-inducibility implying their regulatory role in drought response in rice. Thus, the co-regulatory network analysis facilitated the identification of complex OsMYB regulatory networks, and candidate target regulon genes of selected guide MYB genes. The results contribute to the candidate gene screening, and experimentally testable hypotheses for potential regulatory MYB TFs, and their targets under stress conditions.

Metabolic effects of elevated temperature on organic acid degradation in ripening Vitis vinifera fruit

Berries of the cultivated grapevine Vitis vinifera are notably responsive to temperature, which can influence fruit quality and hence the future compatibility of varieties with their current growing regions. Organic acids represent a key component of fruit organoleptic quality and their content is significantly influenced by temperature. The objectives of this study were to (i) manipulate thermal regimes to realistically capture warming-driven reduction of malate content in Shiraz berries, and (ii) investigate the mechanisms behind temperature-sensitive malate loss and the potential downstream effects on berry metabolism. In the field we compared untreated controls at ambient temperature with longer and milder warming (2-4 °C differential for three weeks; Experiment 1) or shorter and more severe warming (4-6 °C differential for 11 days; Experiment 2). We complemented field trials with control (25/15 °C) and elevated (35/20 °C) day/night temperature controlled-environment trials using potted vines (Experiment 3). Elevating maximum temperatures (4-10 °C above controls) during pre-véraison stages led to higher malate content, particularly with warmer nights. Heating at véraison and ripening stages reduced malate content, consistent with effects typically seen in warm vintages. However, when minimum temperatures were also raised by 4-6 °C, malate content was not reduced, suggesting that the regulation of malate metabolism differs during the day and night. Increased NAD-dependent malic enzyme activity and decreased phosphoenolpyruvate carboxylase and pyruvate kinase activities, as well as the accumulation of various amino acids and γ-aminobutyric acid, suggest enhanced anaplerotic capacity of the TCA cycle and a need for coping with decreased cytosolic pH in heated fruit.

ImmuCo: a database of gene co-expression in immune cells

Current gene co-expression databases and correlation networks do not support cell-specific analysis. Gene co-expression and expression correlation are subtly different phenomena, although both are likely to be functionally significant. Here, we report a new database, ImmuCo (http://immuco.bjmu.edu.cn), which is a cell-specific database that contains information about gene co-expression in immune cells, identifying co-expression and correlation between any two genes. The strength of co-expression of queried genes is indicated by signal values and detection calls, whereas expression correlation and strength are reflected by Pearson correlation coefficients. A scatter plot of the signal values is provided to directly illustrate the extent of co-expression and correlation. In addition, the database allows the analysis of cell-specific gene expression profile across multiple experimental conditions and can generate a list of genes that are highly correlated with the queried genes. Currently, the database covers 18 human cell groups and 10 mouse cell groups, including 20 283 human genes and 20 963 mouse genes. More than 8.6 × 10⁸ and 7.4 × 10⁸ probe set combinations are provided for querying each human and mouse cell group, respectively. Sample applications support the distinctive advantages of the database.

Annotation of gene function in citrus using gene expression information and co-expression networks

BACKGROUND

The genus Citrus encompasses major cultivated plants such as sweet orange, mandarin, lemon and grapefruit, among the world's most economically important fruit crops. With increasing volumes of transcriptomics data available for these species, Gene Co-expression Network (GCN) analysis is a viable option for predicting gene function at a genome-wide scale. GCN analysis is based on a "guilt-by-association" principle whereby genes encoding proteins involved in similar and/or related biological processes may exhibit similar expression patterns across diverse sets of experimental conditions. While bioinformatics resources such as GCN analysis are widely available for efficient gene function prediction in model plant species including Arabidopsis, soybean and rice, in citrus these tools are not yet developed.

RESULTS

We have constructed a comprehensive GCN for citrus inferred from 297 publicly available Affymetrix Genechip Citrus Genome microarray datasets, providing gene co-expression relationships at a genome-wide scale (33,000 transcripts). The comprehensive citrus GCN consists of a global GCN (condition-independent) and four condition-dependent GCNs that survey the sweet orange species only, all citrus fruit tissues, all citrus leaf tissues, or stress-exposed plants. All of these GCNs are clustered using genome-wide, gene-centric (guide) and graph clustering algorithms for flexibility of gene function prediction. For each putative cluster, gene ontology (GO) enrichment and gene expression specificity analyses were performed to enhance gene function, expression and regulation pattern prediction. The guide-gene approach was used to infer novel roles of genes involved in disease susceptibility and vitamin C metabolism, and graph-clustering approaches were used to investigate isoprenoid/phenylpropanoid metabolism in citrus peel, and citric acid catabolism via the GABA shunt in citrus fruit.

CONCLUSIONS

Integration of citrus gene co-expression networks, functional enrichment analysis and gene expression information provide opportunities to infer gene function in citrus. We present a publicly accessible tool, Network Inference for Citrus Co-Expression (NICCE, http://citrus.adelaide.edu.au/nicce/home.aspx), for the gene co-expression analysis in citrus.