Skip to main content
  • Research article
  • Open access
  • Published:

Structure and transcriptional regulation of the major intrinsic protein gene family in grapevine

Abstract

Background

The major intrinsic protein (MIP) family is a family of proteins, including aquaporins, which facilitate water and small molecule transport across plasma membranes. In plants, MIPs function in a huge variety of processes including water transport, growth, stress response, and fruit development. In this study, we characterize the structure and transcriptional regulation of the MIP family in grapevine, describing the putative genome duplication events leading to the family structure and characterizing the family’s tissue and developmental specific expression patterns across numerous preexisting microarray and RNAseq datasets. Gene co-expression network (GCN) analyses were carried out across these datasets and the promoters of each family member were analyzed for cis-regulatory element structure in order to provide insight into their transcriptional regulation.

Results

A total of 29 Vitis vinifera MIP family members (excluding putative pseudogenes) were identified of which all but two were mapped onto Vitis vinifera chromosomes. In this study, segmental duplication events were identified for five plasma membrane intrinsic protein (PIP) and four tonoplast intrinsic protein (TIP) genes, contributing to the expansion of PIPs and TIPs in grapevine. Grapevine MIP family members have distinct tissue and developmental expression patterns and hierarchical clustering revealed two primary groups regardless of the datasets analyzed. Composite microarray and RNA-seq gene co-expression networks (GCNs) highlighted the relationships between MIP genes and functional categories involved in cell wall modification and transport, as well as with other MIPs revealing a strong co-regulation within the family itself. Some duplicated MIP family members have undergone sub-functionalization and exhibit distinct expression patterns and GCNs. Cis-regulatory element (CRE) analyses of the MIP promoters and their associated GCN members revealed enrichment for numerous CREs including AP2/ERFs and NACs.

Conclusions

Combining phylogenetic analyses, gene expression profiling, gene co-expression network analyses, and cis-regulatory element enrichment, this study provides a comprehensive overview of the structure and transcriptional regulation of the grapevine MIP family. The study highlights the duplication and sub-functionalization of the family, its strong coordinated expression with genes involved in growth and transport, and the putative classes of TFs responsible for its regulation.

Background

Aquaporins are channel-forming transmembrane proteins present in plasma and intracellular membranes in all eukaryotes and most prokaryotes [1]. Initially, aquaporins’ water transport capabilities were discovered and functionally characterized in human red blood cells [2,3,4] and later in plants (Arabidopsis thaliana) with the functional characterization of a vacuolar water-transporting protein, γ-TIP [5]. After the discovery of plant aquaporins, many studies have been conducted in order to elucidate their structure, function, and regulation across numerous plant species [6,7,8]. Aquaporins were first characterized as water channels, but they are also recognized to contribute to the transport of other small neutral molecules (e.g., glycerol, urea, boric acid, silicic acid), gases (e.g. CO2, ammonia), and even ions under certain circumstances [7,8,9,10].

Aquaporins fall within an ancient superfamily of membrane proteins called the major intrinsic proteins (MIPs). The MIP family consists of a large number of homologs, and can be subdivided into four major subfamilies based on sequence similarity, which may also indicate their sub-cellular localizations [11, 12]. The plasma membrane intrinsic proteins (PIPs), the tonoplast intrinsic proteins (TIPs), and the nodulin26-like intrinsic proteins (NIPs), comprise the major subfamilies [6, 8, 13]. These three groups of aquaporins have been intensively studied and well-documented. The small basic intrinsic proteins (SIPs) include only a few isoforms localized in the ER (e.g., 3 homologs in Arabidopsis) [9, 14]. In addition to these four well-conserved subfamilies present in all plant species, several additional novel types of aquaporins have been distinguished but with a less ubiquitous presence among plant species. For example, the uncategorized X intrinsic proteins (XIPs) were recently discovered but are absent in some higher plants, including Arabidopsis [15,16,17]. The GlpF-like intrinsic proteins and the hybrid intrinsic proteins were discovered in moss and algae, but are absent in vascular plants [9, 13].

Aquaporins facilitate water transport through plant cells and tissues and play critical rolls in numerous physiological processes. At the cell level, aquaporins act in osmoregulation, reactive oxygen species signaling, and intracellular transport and storage processes [9]. At the tissue and organ level, aquaporins contribute to plant water uptake in roots [18] and facilitate changes in leaf hydraulic conductance [19]. Additionally, aquaporins modulate changes in plant water relations in response to abiotic stress, including drought, salt, and temperature [9]. In fleshy fruit, there is evidence that aquaporins may contribute to ripening processes in tomato [20] and grape [21, 22].

The structure of the MIP gene family, like many plant gene families, has undergone numerous gene duplications resulting in groups of closely related isogenes [11, 23]. These closely related isogenes can have overlapping patterns of expression, or can have undergone sub-functionalization taking on specific developmental and/or tissue related expression patterns [24]. This is certainly the case for MIP family members where many isogenes display tissue and/or developmentally specific expression patterns. Tissue specific expression of MIP isogenes has been observed in numerous species including poplar [25], corn [1, 26, 27], rice [10, 28], Arabidopsis [29], and tomato [20] among other species. On an even finer scale specific isogenes have been associated with specific cell types within organs [19, 30], although most previous studies were not comprehensive across all MIP family members or across organs/tissues.

Grapevine is a plant species of economic and cultural importance and one of the first to have its genome sequenced [31]. This information allowed for the characterization of large gene families such as the MIP family, and indeed this genome information was immediately utilized to integrate cDNA and genome information in characterizing the MIP family members in grapevine [32]. Since then the original Pinot noir genome has been greatly improved and there has been a wealth of microarray and RNAseq studies examining a plethora of conditions (organ specificity, developmental stages, biotic and abiotic stresses, agronomical practices, etc.). Furthermore, new tools and approaches have been developed for analyzing the nature of genome duplications [33], as well as gene expression and cis-regulatory element structure [34]. These improvements allow for a more comprehensive analysis of the grapevine MIP gene family.

In the current study we utilized new tools and approaches to characterize the structure and transcriptional regulation of the MIP gene family in grapevine. We reassessed the MIP family members with the updated genome information describing the putative genome duplication events leading to the current family structure. The expression of family members was then assessed across numerous preexisting microarray and RNAseq datasets in order to determine their tissue and developmental specific expression patterns. Co-expression analyses were carried out across these datasets to determine relevant co-regulation patterns within the MIP gene family and within the transcriptome as a whole. Finally, the promoters of each family member were analyzed for cis-regulatory element structure in order to provide insight into the possible transcriptional regulation of each member.

Methods

Dendrogram construction and gene duplication classification

The grapevine MIP gene family sequences were retrieved from the ORCAE 12× grapevine annotation V2 (http://bioinformatics.psb.ugent.be/orcae/) through a combination of keyword and BLAST searches (using default parameters). For the truncated sequences the surrounding regions were visually inspected for sequence homology to ensure the predicted open reading frames were correct. Gene nomenclature was created following the guidelines established in Grimplet et al. 2014 [35]. Orthology assignment between predicted grapevine MIPs with Arabidopsis proteins was performed using the Conditional Reciprocal Best (CRB)-BLAST method using default settings [36].

Multiple sequence alignments and dendrogram constructions were carried out with Phylogeny.fr [37]. The family was split into sub-families for alignments in order to avoid artifacts caused by aligning large groups [38]. Sequences were aligned with MUSCLE (v3.8.31) using the highest accuracy default settings. After alignment gaps and/or poorly aligned regions were removed using Gblocks (v0.91b) using the following parameters: minimum length of a block after gap cleaning = 5, no gap positions were allowed in the final alignment, all segments with contiguous nonconserved positions bigger than 8 were rejected, minimum number of sequences for a flank position = 55%. Dendrograms were reconstructed using the maximum likelihood method implemented in the PhyML program (v3.1/3.0 aLRT) using default settings. Reliability for internal branch was assessed using the bootstrapping method (100 bootstrap replicates). Dendrograms were drawn with TreeDyn (v198.3).

Analysis of genome structure and duplication analysis was performed using MCScanX [33] using previously established parameters [39]. Information pertaining to the gene duplication type (i.e. singleton, dispersed, proximal, tandem, and segmental; for definition see http://chibba.pgml.uga.edu/mcscan2/), detected collinear pairs, and tandem/proximal gene duplicate groups were further analyzed. Briefly, all genes are initially assigned as ‘singletons’ and ranked (in ascending order) following their positions along chromosomes. Next, all-vs-all BLASTP is performed and results evaluated. The genes with BLASTP hits to other genes are assigned with ‘dispersed’ duplicates. Any two genes are assigned ‘proximal’ duplicates if the difference between gene ranks are < 20 while a rank = 1 between two genes are assigned as ‘tandem’ duplicates. Anchor genes within collinear blocks are assigned as ‘WGD/segmental’ duplicates. In the event where a gene have multiple BLASTP hits, assignment of duplication mode will be in the order of priority beginning with WGD/segmental followed by tandem, proximal, and finally dispersed duplication.

RNA-seq data analysis

Publicly available grapevine next generation sequencing datasets were downloaded from NCBI Sequence Read Archive (http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/sra). Raw fastq reads (single- and paired-end) were extracted using SRA toolkit fastq-dump. Read trimming and quality filtering of reads (single- and paired-end) were performed with Trimmomatic v0.36 [40], with the following parameters; LEADING:20, TRAILING:20 SLIDINGWINDOW:4:20, MINLEN:40, AVGQUAL:20. Alignment of filtered reads towards the 12× grapevine reference genome [31] was performed using HISAT2 v2.0.5 [41] with default parameters. Gene-level count summarization was performed using featureCounts [42] using the grapevine 12× v1 (http://genomes.cribi.unipd.it/) reference annotation and subsequent transcript abundance, expressed as fragments per kilobase of transcript per million mapped reads (FPKM), estimated with edgeR [43].

Gene co-expression network analysis

Two mutual rank (MR) [44] gene co-expression networks (GCN) were constructed, one based on RNA-seq data analyzed in this study and another based on the 29 K NimbleGen whole-genome microarray data. RNA-seq GCN was constructed using log-transformed FPKM values of 29, 970 genes × 237 conditions obtained in this study. Experiment accessions and publication references of analyzed data can be found in Additional file 1: Table S11. Microarray GCN was constructed from an updated input matrix of Wong et al. 2016 [39] containing 29, 000 genes × 358 conditions, an additional of 139 conditions compared to the previous study. Gene-centric co-expression clusters were created for each MIP gene from both RNA-seq and microarray GCNs by considering their top 100 co-expressed genes (ranked by MR value). Visualization of the various MIP networks was carried out in Cytoscape v3.0 [45]. Enrichment of MapMan BIN categories within co-expression clusters were evaluated for enrichment using Fisher’s exact test adjusted with false discovery rate (FDR) for multiple hypothesis correction according to Wong et al. 2016 and 2017 [34, 39]. MapMan BIN categories were considered significantly enriched within co-expression clusters with a FDR < 0.05.

Cis-regulatory element analysis in promoter region

The frequencies and position information of selected cis-regulatory elements (CREs) within 1 kb promoter region from the transcription start site of MIP genes were obtained from Wong et al. 2017 [34] and further analyzed for position bias Z-score considering MIP gene family as a whole/only [46]. The Z-score for each CRE was determined using the equation: Z-score = (L/2 + p)/ √[((L- l + 1)^2–1)/n]. This strategy takes into account the length of the promoter, L; length of the CRE, l; total number of CRE hits present in all promoters, n; and mean position from all identified CRE hits, p. Consideration of these well-established criteria as a whole improves the likelihood of identifying bona fide CREs in selected promoter groups [46].

Results

Family structure

A total of 33 Vitis vinifera MIP family members were identified (Fig. 1; Additional file 1: Table S1). Of these 33 family members we designated 4 of them (VviPIP1–2b, VviPIP2–9, VviNIP9-1a and b) as putative pseudogenes (shown in red in Figs. 1 and 2) because they were both truncated and not expressed in any of the RNAseq datasets we analyzed. These 4 genes were excluded from subsequent analyses in this work. Direct orthologous relationships between Vitis vinifera, poplar, and Arabidopsis are extremely difficult to establish as evidenced by the numerous collapsed dendrogram branches (Fig. 2). We performed additional reciprocal BLAST analyses between the Arabidopsis and Vitis vinifera genes to aid in orthology identification, but again in many cases the orthology could not be resolved (Additional file 1: Table S1; column J “ambiguous”).

Fig. 1
figure 1

Protein sequence relationships within the Vitis vinifera MIP family. The six major MIP sub-families are shown: PIP1s, PIP2s, TIPs, NIPs, SIPs, and XIPs. Red numbers represent bootstrap values (100 bootstrap replicates). Putative pseudogenes are shown in red. Detailed accession, homology, and duplication information is presented in Additional file 1: Table S1

Fig. 2
figure 2

Protein sequence relationships between the Vitis vinifera, Arabidopsis, and poplar MIP families. Six major MIP sub-families: PIP1s (a), PIP2s (b), TIPs (c), NIPs (d), SIPs (e), and XIPs (f). Red numbers represent bootstrap values and the tree was collapsed for all bootstrap values under 50 (100 bootstrap replicates). Linked proteins represent gene duplications for Vitis vinifera (green links) and poplar (blue links as detailed in [25]). Putative pseudogenes are shown in red. Detailed accession, homology, and duplication information is presented in Additional file 1: Table S1

We examined the nature of duplication events contributing to the size of the grapevine MIP gene family (Fig. 2; Additional file 1: Table S1). A total of 27 of 29 grapevine MIP genes were successfully mapped on all 19 grapevine chromosomes. Location of the remaining two, VviTIP2–2 and VviTIP2–3, remains unresolved based on the current 12× genome assembly. In this study, segmental (9 of 29) duplication events were identified for five PIP (VviPIP1–2a, VviPIP1–3, VviPIP2–3, VviPIP2–4, VviPIP2–5) and four TIP (VviTIP1–1, VviTIP1–2, VviTIP1–3, VviTIP1–4) genes, contributing to the expansion of PIPs and TIPs in grapevine. For example, VviPIP2–4 is collinear to both VviPIP2–3 and VviPIP2–5 and VviTIP1–1 is collinear with VviTIP1–2. The PIP duplicates are located on collinear blocks on chromosomes 2/15 and 6/8/13 while TIP duplicates are on chromosome 6/8/13. Meanwhile, tandem duplication was observed for VviPIP1–4 and VviPIP1–2a, where the latter is also a segmental duplicate with VviPIP1–3. Proximal duplication was observed for VviXIP2–1 and VviXIP2–2 where the two are separated by a disease resistance protein. The remaining were classified as dispersed (16 of 29) duplicates whereby the specific mode of duplication is unclear (i.e. other than segmental, tandem, and proximal duplication) and no MIPs were identified as singletons.

Tissue and developmental specific expression and sub-functionalization

Tissue and developmental specific expression profiles of the MIP family members were assessed by examining their expression profiles across the nimblegen grapevine expression atlas [47] (Fig. 3a; Additional file 1: Table S2) and a wide range of existing RNA-seq datasets (Fig. 3b and c; Additional file 1: Table S3). Grapevine MIP family members have distinct tissue and developmental expression patterns. Hierarchical clustering revealed two primary groups (Fig. 3 groups 1 and 2) that were similar regardless of the datasets analyzed. Comparing just the expression atlas (Fig. 3a) with grape berry RNAseq datasets (Fig. 3b), the composition of several subgroups are nearly identical (Fig. 3 sub-groups 3–6).

Fig. 3
figure 3

Expression of the grapevine MIP gene family across the NimbleGen grapevine expression atlas (a) and various other RNAseq datasets in berries (b) and other organs (c; note only “control” states are compared). Colored bars group like tissues or genotypes. Heatmap represents the Z-score according to the scale depicted. Like groupings are numbered for clarity (Groups 1& 2, numbered blue circles, Sub-groups 3–6, numbered white circles). RNAseq experiment accessions and publication references of analyzed data can be found in Additional file 1: Table S11. Raw Z-score values can be found in Additional file 1: File S2 and S3

Generally speaking MIP family members are ubiquitously expressed across tissues, although their expression differs across developmental stages (Fig. 3a). This is true within subfamilies as well with particular isogenes being expressed in almost all tissues, again at specific developmental stages. The inflorescence and flower parts tend to have high levels of MIP expression across the whole family. The primary groups described above (Fig. 3a groups 1 and 2) generally differ in that group 1 is more highly expressed.

Expression across berry development was examined more closely because of the lack of information on aquaporins’ role in fruit development as well as the wealth of datasets available. Of the two primary groups (Fig. 3b groups 1 and 2), group 1 has a much more dynamic expression pattern across berry development regardless of tissue or genotype. In most cases these family members are highly expressed early in berry development and down-regulated as development progresses. However, several members of sub-groups 3 and 4 are up-regulated at the onset of ripening and later during maturation of the berry (Fig. 3b; e.g. VviPIP2–3, VviPIP2–5, VviTIP1–2, VviTIP1–3). In contrast to group 1, group 2 is less dynamic across berry development with a few exceptions, most notably a cluster of family members that exhibit pericarp-specific expression (Fig. 3a sub-group 5).

Duplicated MIP family members exhibit sub-functionalization, at least at the level of their transcriptional regulation. For example, VviXIP2–1 and VviXIP2–2 have distinct expression patterns across a variety of datasets (Fig. 3). This is also true for other examples such as for VviPIP1–4 and VviPIP1–2a. However some duplicated family members exhibit less distinct expression patterns such as VviTIP1–1 and VviTIP1–2.

Enriched functional categories in grapevine MIP gene co-expression networks

To infer the most representative biological functions of this mid-sized gene family, we queried two condition-independent gene co-expression network (GCNs) using individual MIP genes as ‘guides’ separately (Additional file 1: Table S4 and Additional file 1: Table S5) and analyzed their top 100 correlators in detail for biological pathway enrichment (Additional file 1: Table 6 and Additional file 1: Table S7).

In a composite microarray and RNA-seq GCN highlighting MIP genes and their high-level (BIN depth ≤ 1) enriched functional categories (Fig. 4), BIN categories such as transport (BIN34), cell wall (BIN10), miscellaneous enzyme reactions (BIN26) were commonly enriched in MIP co-expression networks. Conversely, categories such as stress (BIN20), cell organization (BIN31), protein metabolism (BIN29), and development (BIN33) were only enriched in specific MIP networks. Specific categories within transport, especially Major Intrinsic Proteins (BIN34.19), were enriched in 17 MIP co-expression networks. Other categories within transport such as ABC/multi-drug transporter (BIN34.16), phosphate (BIN34.7), nitrate (BIN34.4), metal (BIN34.12) were enriched in one (VviNIP7–1), one (VviXIP2–1), two (VviTIP2–3, VviTIP2–2), and three (VviPIP2–4, VviTIP5–1, VviNIP4–1) MIP GCNs, respectively (Additional file 1: Table S6 and Additional file 1: Table S7).

Fig. 4
figure 4

Enriched functional categories identified from composite microarray and RNA-seq MIP gene co-expression networks. Node color represents the parent functional BIN category at BIN depth ≤ 1. Circle size represents the frequency of microarray and/or RNA-seq MIP gene co-expression networks enriched with the corresponding functional BIN category. Solid and dashed edges represent enriched functional BIN category in relevant microarray and RNA-seq MIP gene co-expression networks, respectively

Enrichment of cell wall (BIN10) categories observed in many MIP co-expression networks was unexpected (Fig. 4, Additional file 1: Table S6 and Additional file 1: Table S7). In particular, genes encoding cell wall proteins (BIN10.5), degrading enzymes (BIN10.6), pectin esterases (BIN10.8), cellulose synthesis (BIN10.2), and modification (BIN10.7) belong to categories that were enriched in three (VviTIP5–2, VviPIP2–7, VviTIP1–1), four (VviNIP7–1, VviNIP8–1, VviTIP5–1, VviTIP1–3), three (VviTIP5–1, VviTIP1–1, VviTIP1–3), two (VviPIP2–7, VviPIP2–5), and one (VviTIP1–3) MIP GCNs, respectively. The two categories that deserve attention within the miscellaneous (BIN26) category relate to the enrichment of glutathione-S-transferase and peroxidase co-expressed genes in four (VviTIP2–3, VviTIP2–2, VviPIP2–4, VviPIP1–1) and five (VviTIP1–4, VviNIP1–2, VviTIP2–3, VviTIP2–2, VviTIP5–1) MIP GCNs, respectively. Meanwhile, seven MIP GCNs (VviNIP8–1, VviXIP2–1, VviTIP1–4, VviTIP2–1, VviPIP1–3, VviTIP2–2, VviTIP2–3) enriched with abiotic stress related genes (BIN20.2) are also of interest.

Divergence of enriched functional categories in gene co-expression networks of grapevine MIP duplicates

As a significant proportion of grapevine MIP duplicates showed sub-functionalization of gene expression in a tissue-specific manner (Fig. 3), we compared the enriched functional categories in the GCNs (Additional file 1: Table S6 and Additional file 1: Table S7) of MIP duplicate pairs (Additional file 1: Table S1). GCNs of duplicates such as VviXIP2–1 and VviXIP2–2 have totally distinct enriched categories. While major intrinsic proteins (BIN34.19) and abiotic stress (BIN20.2) genes were enriched in VviXIP2–1 GCNs, the latter two categories were absent in VviXIP2–2 GCN. Instead, genes involved in the light reaction of photosynthesis (BIN1.1), isoprenoid metabolism (BIN16.1), auxin metabolism (BIN17.2), and CYP450-coding genes (BIN26.1) were enriched in VviXIP2–2 GCNs. In another example, GCNs of VviTIP1–3 and VviTIP1–4 duplicate pairs have in common enrichment for major intrinsic proteins (BIN34.19), however, enrichment of cell wall pectin esterases (BIN10.8) and modification (BIN10.7) was observed for VviTIP1–3 while abiotic stress (BIN20.2), cell cycle (BIN31.3), and hormone (i.e. JA and ABA) metabolism (BIN17) functional categories were among the many categories enriched in the GCN of its duplicate VviTIP1–4.

Conversely, duplicated family members that exhibit less divergent expression profiles such as VviPIP2–4 and VviPIP2–5 showed more commonalities. Both VviPIP2–4 and VviPIP2–5 have a common enrichment for cell wall (BIN10) and major intrinsic proteins (BIN34.19), albeit some differences were apparent such as enrichment for glutathione S transferases (BIN26.9) and cell organization (BIN31.1) in VviPIP2–4 and VviPIP2–5, respectively (Additional file 1: Table S6 and Additional file 1: Table S7). Similarly, VviTIP1–1 and VviTIP1–2 share enrichment for cell wall (BIN10) related genes, but categories related to light reaction (BIN1.1) and cell organization (BIN31.1) were enriched in VviTIP1–1 and VviTIP1–2, respectively.

Cis-regulatory element structure of grapevine MIP promoters

Genome-wide analysis in grapevine promoters have highlighted many CREs possessing strong position bias towards the transcription start site (TSS) which were implicated in a variety of grapevine development and stress responses [34]. To determine which CREs are biologically relevant for the regulation of grapevine MIPs we extracted the distribution patterns of 222 CREs (6- to 8-mer in length) in the promoter region for grapevine MIPs (Additional file 1: Table S8) selected from Wong et al. 2017 [34]. The frequency of occurrence, the median position of occurrence, and position bias Z-score were evaluated. On these subset of MIP genes, 6-mer and 7-mer CREs namely RYCGAC, YAACKG, TTRCGT, and ACGTGKC were amongst top 10 most highly ranked CREs based Z-score (Fig. 5; Additional file 1: Table S9). The most highly ranked CRE, the RYCGAC – part/variant of the dehydration-responsive element (DRE)/C-repeat elements/low-temperature-responsive element [48] – were present in 14 MIP promoters (∑ hits: 23, M position: 262) followed by YAACKG CRE – part/variant of the type I R2R3-MYB recognition sites [49] – that were present in 19 MIP promoters (∑ hits: 35, M position: 315). The TTRCGT CRE – the major NAC TF recognition sites [50] – was also ranked highly and was present in 9 MIP promoters (∑ hits: 13, M position: 286). Longer CREs such as ACGTGKC [51] – a well-known ABA-responsive element (ABRE) – were present in 8 MIP promoter (∑ hits: 12, M position: 226). For most of these CREs, a position bias towards the TSS (M position < 300) was also observed considering MIP genes.

Fig. 5
figure 5

Number and location of highlighted cis-regulatory elements in the promoters of grapevine MIP family members. The AP2/ERF (orange, blue, pink, green), bZIP (turquoise), NAC (yellow), and R2R3–MYB (beige) were amongst top 10 (of 222) most highly ranked CREs (6- to 8-mer in length) based Z-score. Each occurrence of the CRE is noted at its position with the appropriate colored line. Complete promoter CRE data can be found in Additional file 1: Table 8 and Additional file 1: Table S9

Enrichment of known cis-regulatory elements in grapevine MIP gene co-expression networks

In this study, promoters of genes within MIP GCNs were also tested for enrichment for known CREs in order to identify putative shared TF families within the MIP GCNs that may be responsible for their transcriptional regulation. Nineteen MIP GCNs displayed significant enrichment (FDR < 0.01) for at least one CRE tested (Fig. 6; Additional file 1: Table S10). Of these, six MIP GCNs (i.e. VviNIP8–1, VviPIP1–1, VviXIP2–1, VviTIP1–4, VviTIP2–2, and VviTIP2–3) were commonly enriched for the PHR1-binding sequence (P1BS, GNATATNC). Several of these MIP GCNs were also co-enriched with other CREs. For example, AP2/ERF (GCCGGC) and R2R3-MYB (GKTKGTTR) related CREs were observed in the VviNIP8–1 GCN along with related genes such as AP2/ERF TFs, VviTOE3 (VIT_01s0026g01690), VviERF1L4 (VIT_07s0005g03270), and VviMYB82C (VIT_11s0016g05690). Promoters of the VviXIP2–1 GCN are also enriched for the GCCGGC CRE correlating with the presence of two AP2/ERF TFs, VviTOE2 (VIT_14s0108g00050) and VviTOE3 (VIT_01s0026g01690). Other co-enriched CREs of interest include HB (CAATWATT) and extended DRE elements (DEAR4, CRCCGACA) in promoters of the VviPIP1–4 GCN, coincident with three HB TFs (VIT_08s0007g01290, VIT_16s0100g00670, VIT_18s0001g08410) and two AP2/ERF TFs VviERF061 (VIT_02s0025g01360) and VviERF022 (VIT_18s0001g05850).

Fig. 6
figure 6

Enriched cis-regulatory elements within shared MIP co-expression networks. Enriched CREs between 6- and 8-mer are depicted as separate panels. Circle opacity represents the enrichment score (−log10 FDR values) of the corresponding enriched CRE and circle size of represents the total number of genes containing the enriched CRE

Interestingly, VviTIP3–1 and VviPIP1 − 2c were the two GCNs enriched with many of the CREs tested. For example, many ACGT-related (e.g. ACGTABREMOTIFA2OSEM, GADOWNAT, CACGTGMOTIF, ABREATCONSENSUS) and RY-related (i.e. RYREPEATBNNAPA, RYREPEATLEGUMINBOX) CREs, but also others involved in Ca2+/calmodulin signaling (CGCGBOXAT) and lateral organ boundary TF binding (TCCGGA) were enriched in the VviTIP3–1 GCN specifically. A suite of TFs whose homologs target the latter CREs were also present reaffirming the biological relevance of this broad enrichment pattern (Fig. 6; Additional file 1: Table S10). This includes two B3 TFs, VviABI3 (VIT_07s0005g05400) and VviFUS3 (VIT_14s0068g01290) whose homologs/orthologs in Arabidopsis target the RY motif [52] and VviABI5/VvbZIP25 (VIT_08s0007g03420) that targets various ACGT-related CREs.

Discussion

Expansion and sub-functionalization of the grapevine MIP family

The number of grapevine MIPs identified (33) is similar to the number identified in earlier versions of the grape genome assembly (29 MIPs, [32]), Arabidopsis (35 MIPs, [11]), and rice (33 MIPs, [10]). The MIP family is highly conserved and although many orthologous grapevine-Arabidopsis pairs were identified more than half of the orthologous relationships were impossible to resolve. The annotation and gene names presented here differ at times (8 of 33) from those established in Shelden et al. (2009) [32]. This is primarily due to a much improved genome assembly which allowed for the identification of previously unidentified family members, and improved computational methods for identifying the most likely Arabidopsis orthologs [36].

The grapevine MIP gene family has undergone a number of duplication events consistent with the highly duplicated nature of plant genomes and grapevine specifically [31]. The duplication events concerning segmental and tandem duplications identified in this study have also been reported for Arabidopsis [23] and rice [53, 54]. Nonetheless, novel duplication events involving VviXIP2–1 and VviXIP2–2 may be grapevine-specific.

It is commonplace that duplicated genes often take on different expression patterns, with respect to specific portions of development and/or location [24]. In the current study some duplicated MIP gene family members have distinct patterns of expression. It is likely that these duplicates have a similar protein function yet function in different contexts, for example VviXIP2–2 in leaves and VviXIP2–1 in roots (see Fig. 3c). In grapevine, several other gene families have a similar history of duplication and sub-functionalization [39, 55, 56]. Concerning fruit specifically, the expression of most MIP family members decreases as berry development progresses consistent with earlier studies [21, 22]. Grape berries become increasingly hydraulically buffered from the parent plant during ripening and this buffering is thought to result in part from decreases in hydraulic conductivity [21, 57]. This general downregulation of MIP family members during ripening may contribute to these decreases in berry hydraulic conductivity. In contrast, some specific isogenes (e.g. VviPIP2–3 and VviPIP2–5; note PIP2–5 was previously referred to as PIP2–1) show significant expression and even up-regulation throughout the later stages of berry development [21, 22]. Their role in fruit ripening remains unknown, but some have speculated that they may facilitate small ion transport and/or osmoregulation [58]. Fleshy fruits like grape berries undergo rapid growth and sugar accumulation during ripening and the role of aquaporins in mediating grape berry water relations is certainly worthy of further study [59].

The grapevine expression atlas [47] is a powerful dataset for examining tissue and developmental specific expression patterns however caution is warranted especially when examining highly conserved gene families. Microarray based expression analyses can be biased via cross-hybridization [60], and this is why it is important to include RNAseq based analyses as well. The results of the MIP family members presented here show strong parallels between both approaches suggesting that any potential cross-hybridization did not lead to erroneous results in the case of the expression atlas.

Grapevine MIP co-regulation networks

Based on the ‘guilt-by-association’ principle, genes involved in related processes often share parallel expression dynamics across a wide range conditions including different organ/cell types, developmental stages, stress, and hormonal perturbations [61]. Gene co-expression networks (GCNs) analyses, which are built upon the ‘guilt-by-association’ principle, have been particularly useful for ascribing the most representative biological functions to both individual gene(s) [62,63,64,65] and large gene families [39, 66] in grapevine. This study highlights the strong co-expression relationships within the MIP family itself, and between MIP family members and genes involved various processes such as growth, cell-division, and cell redox homeostasis.

One of the strongest GCN relationships revealed in this study was that between the MIP family and genes involved in growth and transport processes, namely cell wall modification and cell expansion. Aquaporins have been implicated in the growth of rose flower petals and are part of a GCN associated with petal cell expansion [67]. In grape berries, targeted analyses of a limited number of aquaporins and cell wall metabolic genes were shown to have similar patterns of expression that correlated with growth [68]. The treatment of grape berries with exogenous ethylene stimulated growth and associated micro-array analyses revealed coordinated changes in the expression of both aquaporin and cell wall metabolic genes [69]. Among the cell wall metabolic genes identified by Schlosser et al. (2008) [68] and Chervin et al. (2008) [69] were the pectin esterases (BIN10.8) and cellulose synthesis (BIN10.2) identified in this study. The congruence between these previous studies and the more global a priori approach utilized here provides robust evidence for a functional link between these groups of genes.

Our GCN analyses also revealed a strong link between the MIP family and cell division (cell cycle, BIN31.3, and cell organization, BIN31.1). This is a relationship that has not been studied in plants apart from a few studies. Over-expression of tobacco NtTIP1;1 in cell culture enhanced cellular expansion and cell-division [70] and specific aquaporin isoforms have been associated with rapidly proliferating tissues in roots [71, 72]. Cell proliferation and growth involves the regulation of source-sink relationships intersecting with turgor driven growth, and one could speculate an important role for MIP family members in both of these processes. Outside of plants there is a growing body of work linking aquaporin function with the regulation of cell proliferation [73].

Another interesting GCN highlighted in this study was between the MIP family and cell redox homeostasis. The most obvious link between MIPs and redox homeostasis is the fact that many MIP isoforms transport hydrogen peroxide [74]. Therefore perhaps it should not come as a surprise that MIP family members would be among coordinated redox homeostasis genes. Links between aquaporin function and redox homeostasis are involved in the regulation of root water uptake under stress [75,76,77] but not necessarily through a transcriptional mechanism [78] and the same is true for pathogen responses [79]. Perhaps one of the most interesting observations is the nexus between cell expansion, cell division, and redox homeostasis [80], where aquaporins may play a cornerstone role in coordinating water fluxes and redox homeostasis in the control of growth.

The diversity of bona fide cis-regulatory elements in grapevine MIP promoters

Regulation of plant MIP genes is still poorly understood. This study represents a first attempt of characterizing the CRE structure of the grapevine MIP family and identifying putative TFs responsible for its regulation. As limitations exists even for well-established statistical measures used for prioritizing CREs, combining several metrics may overcome potential caveats of each approach [34].

Recent studies have shown that the DRE and GCC-box (GCCGCC) core sequences are critical for the regulation of MIP genes by members of AP2/ERF subgroups I, IV, and V in several plants [81,82,83,84]. This is consistent with the highly prioritized DRE in MIP promoters among all other CREs (Fig. 5) and the co-regulation with AP2/ERF TFs including several predicted grapevine subgroup I, IV, and V members. Some of these regulatory relationships are conserved while many others are novel. Known relationships include co-regulation of a closely related grapevine homolog of Arabidopsis RAP2.11 (VIT_02s0025g03170) with VviTIP2–1, and co-regulation of VviTIP3–1 with grapevine homolog of Arabidopsis DREB2D. These examples of co-regulation are consistent with known and predicted targets of Arabidopsis RAP2.11 [85] and DREB2D [86]. The DRE sites within VviTIP2–1 and VviTIP3–1 promoter may be important for its regulation in grapevine. Unexpectedly, GCC-box elements and other GCC-related CREs (GCCGGC, GCCGTC) were not found within most MIP promoters within 1000 bases from the TSS (Additional file 1: Table S8 and Additional file 1: Table S9). This observation might indicate potential divergence in AP2/ERF transcriptional regulatory networks involving MIP genes between plant species and the DRE may be preferred in grapevines, and/or that GCC-box and related CREs are located beyond the promoter regions analyzed in this study.

Several bZIP, NAC, and R2R3-MYB transcription factors have been shown to regulate specific MIP genes [67, 87, 88] consistent with many highly co-regulated TFs of these families in MIP subnetworks. Differences in the distribution of CREs present in MIP promoters were also observed (Fig. 5). PIP and TIP promoters contain mostly AP2/ERF and bZIP-related CREs while NIP promoters contain mostly NAC and R2R3-MYB-related CREs, suggesting some degree of transcriptional regulation specificity in grapevine aquaporin regulation.

Promoter analysis suggests that hormone metabolic pathways such as ABA and ethylene play an important role in the regulation of MIP genes. There is evidence for ethylene-regulated aquaporin expression in rose petals [67, 89] and aquaporin genes are among those regulated by exogenous ethylene treatment in grape berries [69]. Several studies demonstrate that ABA regulates the expression of numerous MIP family members [90,91,92]. However, it is important to point out that short-term modulation of aquaporin activity via ABA, and possibly other hormones such as ethylene, likely occurs at the post-translational level [93, 94]. In grapevine, ABA has been shown to differentially regulate the same aquaporin isogene (VviTIP1–1, VIT_06s0061g00730) depending on the organ [95]. These complex relationships involved in the hormonal regulation of aquaporin gene expression require further study.

The promoters of genes in six MIP GCNs were also commonly enriched for the PHR1-binding sequence (P1BS, GNATATNC). The cognate sequence, such as GAATATTC, is known to be bound by members of the MYB (GARP, G2) TF and is related to the regulation of transcriptional repressors [96] and accordingly no homologs of MYB (GARP, G2) TFs were represented in respective MIP GCNs. The enrichment of this CRE may suggest a potential role of large-scale transcriptional repression in the regulation of MIPs. Conversely, for many other CREs enrichment profiles were often accompanied by the presence of TF families known to target them (Fig. 6, Additional file 1: Table S10) suggesting a role of transcriptional activation of MIP and co-regulated genes. The diversity of enriched CREs also highlights that in addition to the those shown to be directly implicated in MIP regulation such as AP2/ERF, bZIP, NAC, and R2R3-MYB TFs, regulation of MIPs may involve more TF families than previously described. Several genes that belong to HB, LBD, and B3 TF families may also represent novel candidate regulators of grapevine MIP and co-regulated genes.

Conclusions

The current work utilized the most high quality and up-to-date genome information in characterizing the grapevine MIP gene family, its structure, and the putative duplication events involved in its evolution. When paired with the GCN analyses conducted here we identified those MIP family members that have undergone duplication and sub-functionalization through characterizing the tissue and developmental specific expression patterns across the family. GCN analyses revealed several interesting relationships between MIP family members and genes involved in cell expansion, cell division, and transport processes. Characterizing the cis-regulatory elements in grapevine MIP promoters along with associated GCN members identified the putative classes of TFs responsible for the regulation of the family and their associated GCNs. Combining phylogenetic analyses, gene expression profiling, GCN analyses, and CRE enrichment, this study provides a comprehensive overview of the structure and transcriptional regulation of the grapevine MIP family. These results can help guide future studies aimed at understanding the role of specific transcription factors in controlling the diverse expression patterns within the MIP family.

Abbreviations

CRE:

Cis-Regulatory Element

DRE:

Dehydration-Responsive Element

GCN:

Gene Co-expression Network

MIP:

Major Intrinsic Protein

NIP:

Nodulin26-like Intrinsic Proteins

PIP:

Plasma membrane Intrinsic Protein

SIP:

Small basic Intrinsic Proteins

TIP:

Tonoplast Intrinsic Protein

TSS:

Transcription Start Site

XIP:

X Intrinsic Proteins

References

  1. Chaumont F, Barrieu F, Wojcik E, Chrispeels MJ, Jung R. Aquaporins constitute a large and highly divergent protein family in maize. Plant Physiol. 2001;125:1206–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Benga G, Popescu O, Borza V, Pop VI, Muresan A, Mocsy I, et al. Water permeability in human erythrocytes: identification of membrane proteins involved in water transport. Eur J Cell Biol. 1986;41:252–62.

    CAS  PubMed  Google Scholar 

  3. Denker BM, Smith BL, Kuhajda FP, Agre P. Identification, purification, and partial characterization of a novel Mr 28,000 integral membrane protein from erythrocytes and renal tubules. J Biol Chem. 1988;263:15634–42.

    CAS  PubMed  Google Scholar 

  4. Preston GM, Carroll TP, Guggino WB, Agre P. Appearance of water channels in Xenopus oocytes expressing red cell CHIP28 protein. Science. 1992;256:385–7.

    Article  CAS  PubMed  Google Scholar 

  5. Maurel C, Reizer J, Schroeder JI, Chrispeels MJ. The vacuolar membrane protein gamma-TIP creates water specific channels in Xenopus oocytes. EMBO J. 1993;12:2241–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  6. Chaumont F, Tyerman SD. Aquaporins: highly regulated channels controlling plant water relations. Plant Physiol. 2014;164:1600–18.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Tyerman SD, Niemietz CM, Bramley H. Plant aquaporins: multifunctional water and solute channels with expanding roles. Plant Cell Environ. 2002;25:173–94.

    Article  CAS  PubMed  Google Scholar 

  8. Maurel C, Verdoucq L, Luu D-T, Santoni V. Plant aquaporins: membrane channels with multiple integrated functions. Annu Rev Plant Biol. 2008;59:595–624.

    Article  CAS  PubMed  Google Scholar 

  9. Maurel C, Boursiac Y, Luu D-T, Santoni V, Shahzad Z, Verdoucq L. Aquaporins in plants. Physiol Rev. 2015;95:1321–58.

    Article  CAS  PubMed  Google Scholar 

  10. Sakurai J, Ishikawa F, Yamaguchi T, Uemura M, Maeshima M. Identification of 33 rice aquaporin genes and analysis of their expression and function. Plant Cell Physiol. 2005;46:1568–77.

    Article  CAS  PubMed  Google Scholar 

  11. Johanson U, Karlsson M, Johansson I, Gustavsson S, Sjövall S, Fraysse L, et al. The complete set of genes encoding major intrinsic proteins in Arabidopsis provides a framework for a new nomenclature for major intrinsic proteins in plants. Plant Physiol. 2001;126:1358–69.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Alexandersson E, Fraysse L, Sjövall-Larsen S, Gustavsson S, Fellert M, Karlsson M, et al. Whole gene family expression and drought stress regulation of Aquaporins. Plant Mol Biol. 2005;59:469–84.

    Article  CAS  PubMed  Google Scholar 

  13. Li G, Santoni V, Maurel C. Plant aquaporins: roles in plant physiology. Biochim. Biophys. Acta - Gen. Subj. 1840;2014:1574–82.

    Google Scholar 

  14. Ishikawa F, Suga S, Uemura T, Sato MH, Maeshima M. Novel type aquaporin SIPs are mainly localized to the ER membrane and show cell-specific expression in Arabidopsis thaliana. FEBS Lett. 2005;579:5814–20.

    Article  CAS  PubMed  Google Scholar 

  15. Gupta A, Sankararamakrishnan R. Genome-wide analysis of major intrinsic proteins in the tree plant Populus trichocarpa: characterization of XIP subfamily of aquaporins from evolutionary perspective. BMC Plant Biol. 2009;9:134.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Danielson JÅ, Johanson U. Unexpected complexity of the aquaporin gene family in the moss Physcomitrella patens. BMC Plant Biol. 2008;8:45.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Lopez D, Bronner G, Brunel N, Auguin D, Bourgerie S, Brignolas F, et al. Insights into Populus XIP aquaporins: evolutionary expansion, protein functionality, and environmental regulation. J Exp Bot. 2012;63:2217–30. Oxford University Press

    Article  CAS  PubMed  Google Scholar 

  18. Gambetta GA, Knipfer T, Fricke W, McElrone AJ. Aquaporins and root water uptake. Cham: Springer; 2017. p. 133–53.

    Google Scholar 

  19. Heinen RB, Ye Q, Chaumont FX. Role of aquaporins in leaf physiology. J Exp Bot. 2009;60:2971–85.

    Article  CAS  PubMed  Google Scholar 

  20. Reuscher S, Akiyama M, Mori C, Aoki K, Shibata D, Shiratake K. Genome-Wide Identification and Expression Analysis of Aquaporins in Tomato. Boudko D, editor. PLoS One. Public Libr Sci; 2013;8:e79052.

  21. Choat B, Gambetta GA, Shackel KA, Matthews MA. Vascular function in grape berries across development and its relevance to apparent hydraulic isolation. Plant Physiol. 2009;151:1677–87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Fouquet R, Léon C, Ollat N, Barrieu F. Identification of grapevine aquaporins and expression analysis in developing berries. Plant Cell Rep. 2008;27:1541–50.

    Article  CAS  PubMed  Google Scholar 

  23. Cannon SB, Mitra A, Baumgarten A, Young ND, May G. The roles of segmental and tandem gene duplication in the evolution of large gene families in Arabidopsis thaliana. BMC Plant Biol. 2004;4:10. BioMed Central

    Article  PubMed  PubMed Central  Google Scholar 

  24. Adams KL, Wendel JF. Novel patterns of gene expression in polyploid plants. Sci Trends Genet. 1999;298:2157–67.

    Google Scholar 

  25. Cohen D, Bogeat-Triboulot M-B, Vialet-Chabrand S, Merret R, Courty P-E, Moretti S, et al. Developmental and environmental regulation of aquaporin gene expression across Populus species: divergence or redundancy? Blazquez MA, editor. PLoS One. 2013;8:e55506. Public Library of Science

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Gaspar M, Bousser A, Sissoeff I, Roche O, Hoarau J, Mahe A. Cloning and characterization of ZmPIP1-5b, an aquaporin transporting water and urea. Plant Sci. 2003;165:21–31.

    Article  CAS  Google Scholar 

  27. Opitz N, Marcon C, Paschold A, Malik WA, Lithio A, Brandt R, et al. Extensive tissue-specific transcriptomic plasticity in maize primary roots upon water deficit. J Exp Bot. 2016;67:1095–107.

    Article  CAS  PubMed  Google Scholar 

  28. Sakurai J, Ahamed A, Murai M, Maeshima M, Uemura M. Tissue and cell-specific localization of rice aquaporins and their water transport activities. Plant Cell Physiol. 2008;49:30–9.

    Article  CAS  PubMed  Google Scholar 

  29. Weig A, Deswarte C, Chrispeels MJ. The major intrinsic protein family of Arabidopsis has 23 members that form three distinct groups with functional aquaporins in each group. Plant Physiol. 1997;114:1347–57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Kirch H-H, Vera-Estrella R, Golldack D, Quigley F, Michalowski CB, Barkla BJ, et al. Expression of Water Channel proteins in Mesembryanthemum crystallinum. Plant Physiol. 2000;123:111–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Jaillon O, Aury J-M, Noel B, Policriti A, Clepet C, Casagrande A, et al. The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007;449:463–7.

    Article  CAS  PubMed  Google Scholar 

  32. Shelden MC, Howitt SM, Kaiser BN, Tyerman SD. Identification and functional characterisation of aquaporins in the grapevine, Vitis vinifera. Funct Plant Biol. 2009;36:1065–78.

    Article  CAS  Google Scholar 

  33. Wang Y, Tang H, Debarry JD, Tan X, Li J, Wang X, et al. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012;40:e49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Wong DCJ, Lopez Gutierrez R, Gambetta GA, Castellarin SD. Genome-wide analysis of cis-regulatory element structure and discovery of motif-driven gene co-expression networks in grapevine. DNA Res. 2017;24:311–26.

    CAS  PubMed  PubMed Central  Google Scholar 

  35. Grimplet J, Adam-Blondon A-F, Bert P-F, Bitz O, Cantu D, Davies C, et al. The grapevine gene nomenclature system. BMC Genomics. 2014;15:1077.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Aubry S, Kelly S, Kümpers BMC, Smith-Unna RD, Hibberd JM. Deep evolutionary comparison of gene expression identifies parallel recruitment of trans-factors in two independent origins of C4 photosynthesis. Bomblies K, editor. PLoS Genet. 2014;10:e1004365. Public Library of Science

    Article  PubMed  PubMed Central  Google Scholar 

  37. Dereeper A, Guignon V, Blanc G, Audic S, Buffet S, Chevenet F, et al. Phylogeny.Fr: robust phylogenetic analysis for the non-specialist. Nucleic Acids Res. 2008;36:W465–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Boyce K, Sievers F, Higgins DG. Instability in progressive multiple sequence alignment algorithms. Algorithms Mol Biol. 2015;10:26.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Wong DCJ, Schlechter R, Vannozzi A, Höll J, Hmmam I, Bogs J, et al. A systems-oriented analysis of the grapevine R2R3-MYB transcription factor family uncovers new insights into the regulation of stilbene accumulation. DNA Res. 2016;23:451–66.

    Article  CAS  PubMed Central  Google Scholar 

  40. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12:357–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30:923–30.

    Article  CAS  PubMed  Google Scholar 

  43. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.

    Article  CAS  PubMed  Google Scholar 

  44. Obayashi T, Kinoshita K. Rank of correlation coefficient as a comparable measure for biological significance of gene coexpression. DNA Res. 2009;16:249–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Ma S, Shah S, Bohnert HJ, Snyder M, Dinesh-Kumar SP. Incorporating motif analysis into gene co-expression networks reveals novel modular expression pattern and new signaling pathways. Copenhaver GP, editor. PLoS Genet. 2013;9:e1003840.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Fasoli M, Dal Santo S, Zenoni S, Tornielli GB, Farina L, Zamboni A, et al. The grapevine expression atlas reveals a deep transcriptome shift driving the entire plant into a maturation program. Plant Cell. 2012;24:3489–505.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Yamaguchi-Shinozaki K, Shinozaki K. Transcriptional regulatory networks in cellular responses and tolerance to dehydration and cold stresses. Annu Rev Plant Biol. 2006;57:781–803.

    Article  CAS  PubMed  Google Scholar 

  49. Prouse MB, Campbell MM. The interaction between MYB proteins and their target DNA binding sites. Biochim Biophys Acta - Gene Regul Mech. 1819;2012:67–77.

    Google Scholar 

  50. Lindemose S, Jensen MK, Van de Velde J, O’Shea C, Heyndrickx KS, Workman CT, et al. A DNA-binding-site landscape and regulatory network analysis for NAC transcription factors in Arabidopsis thaliana. Nucleic Acids Res. 2014;42:7681–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Narusaka Y, Nakashima K, Shinwari ZK, Sakuma Y, Furihata T, Abe H, et al. Interaction between two cis-acting elements, ABRE and DRE, in ABA-dependent expression of Arabidopsis rd29A gene in response to dehydration and high-salinity stresses. Plant J. 2003;34:137–48.

    Article  CAS  PubMed  Google Scholar 

  52. Reidt W, Wohlfarth T, Ellerström M, Czihal A, Tewes A, Ezcurra I, et al. Gene regulation during late embryogenesis: the RY motif of maturation-specific gene promoters is a direct target of the FUS3 gene product. Plant J. 2000;21:401–8.

    Article  CAS  PubMed  Google Scholar 

  53. Guyot R, Keller B. Ancestral genome duplication in rice. Genome. 2004;47:610–4.

    Article  CAS  PubMed  Google Scholar 

  54. Yu J, Wang J, Lin W, Li S, Li H, Zhou J, et al. The genomes of Oryza sativa: a history of duplications. Bennetzen J, editor. PLoS Biol Sinauer. 2005;3:e38.

    Article  Google Scholar 

  55. Falginella L, Castellarin SD, Testolin R, Gambetta GA, Morgante M, Di Gaspero G. Expansion and subfunctionalisation of flavonoid 3′,5′-hydroxylases in the grapevine lineage. BMC Genomics. 2010;11:562.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Vannozzi A, Dry IB, Fasoli M, Zenoni S, Lucchin M, Henikoff S, et al. Genome-wide analysis of the grapevine stilbene synthase multigenic family: genomic organization and expression profiles upon biotic and abiotic stresses. BMC Plant Biol. 2012;12:130. BioMed Central

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Knipfer T, Fei J, Gambetta GA, McElrone AJ, Shackel KA, Matthews MA. Water transport properties of the grape pedicel during fruit development: insights into xylem anatomy and function using microtomography. Plant Physiol. 2015;168:1590–602. American Society of Plant Biologists

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Rogiers SY, Coetzee ZA, Walker RR, Deloire A, Tyerman SD. Potassium in the grape (Vitis vinifera L.) berry: transport and function. Front. Plant Sci. 2017;8:1629.

    Google Scholar 

  59. Tyerman SD, Chaves MM, Barrieu F. Water Relations of the Grape Berry and Aquaporins. In: Gerós H, Chaves MM, Delrot S, editors. Biochem. grape berry. Potomac: Bentham Science; 2012. p. 3–22.

  60. Cramer GR, Ghan R, Schlauch KA, Tillett RL, Heymann H, Ferrarini A, et al. Transcriptomic analysis of the late stages of grapevine (Vitis vinifera cv. Cabernet sauvignon) berry ripening reveals significant induction of ethylene signaling and flavor pathways in the skin. BMC Plant Biol. 2014;14:370. BioMed Central

    Article  PubMed  PubMed Central  Google Scholar 

  61. Gillis J, Pavlidis P. “Guilt by association” is the exception rather than the rule in gene networks. Rzhetsky a, editor. PLoS Comput Biol. 2012;8:e1002444.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Loyola R, Herrera D, Mas A, Wong DCJ, Höll J, Cavallini E, et al. 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. 2016;67:5429–45.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Amato A, Cavallini E, Zenoni S, Finezzo L, Begheldo M, Ruperti B, et al. A grapevine TTG2-like WRKY transcription factor is involved in regulating vacuolar transport and flavonoid biosynthesis. Front Plant Sci. 2017;7:1979. Frontiers.

    Article  PubMed  PubMed Central  Google Scholar 

  64. Malacarne G, Coller E, Czemmel S, Vrhovsek U, Engelen K, Goremykin V, et al. The grapevine VvibZIPC22 transcription factor is involved in the regulation of flavonoid biosynthesis. J Exp Bot. 2016;67:3509–22. Oxford University Press.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Vannozzi A, Wong DCJ, Höll J, Hmmam I, Matus JT, Bogs J, et al. Combinatorial Regulation of Stilbene Synthase Genes by WRKY and MYB Transcription Factors in Grapevine (Vitis vinifera L.). Plant Cell Physiol. 2018.

  66. Wong DCJ, Ariani P, Castellarin S, Polverari A, Vandelle E. Co-expression network analysis and cis-regulatory element enrichment determine putative functions and regulatory mechanisms of grapevine ATL E3 ubiquitin ligases. Sci Rep. 2018;8:3151.

    Article  PubMed  PubMed Central  Google Scholar 

  67. Pei H, Ma N, Tian J, Luo J, Chen J, Li J, et al. An NAC transcription factor controls ethylene-regulated cell expansion in flower petals. Plant Physiol. 2013;163:775–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Schlosser J, Olsson N, Weis M, Reid K, Peng F, Lund S, et al. Cellular expansion and gene expression in the developing grape (Vitis vinifera L.). Protoplasma. 2008;232:255–65.

    Article  CAS  PubMed  Google Scholar 

  69. Chervin C, Tira-umphon A, Terrier N, Zouine M, Severac D, Roustan J-P. Stimulation of the grape berry expansion by ethylene and effects on related gene transcripts, over the ripening phase. Physiol Plant. 2008;134:534–46.

    Article  CAS  PubMed  Google Scholar 

  70. Okubo-Kurihara E, Sano T, Higaki T, Kutsuna N, Hasezawa S. Acceleration of vacuolar regeneration and cell growth by overexpression of an aquaporin NtTIP1;1 in tobacco BY-2 cells. Plant Cell Physiol. 2009;50:151–60. Oxford University Press

    Article  CAS  PubMed  Google Scholar 

  71. Chaumont F, Barrieu F, Herman EM, Chrispeels MJ. Characterization of a maize tonoplast aquaporin expressed in zones of cell division and elongation. Plant Physiol. 1998;117:1143–52. American Society of Plant Biologists

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Gambetta GA, Fei J, Rost TL, Knipfer T, Matthews MA, Shackel KA, et al. Water uptake along the length of grapevine fine roots: developmental anatomy, tissue-specific aquaporin expression, and pathways of water transport. Plant Physiol. 2013;163:1254–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Galán-Cobo A, Ramírez-Lorca R, Echevarría M. Role of aquaporins in cell proliferation: what else beyond water permeability? Channels. 2016;10:185–201.

    Article  PubMed  PubMed Central  Google Scholar 

  74. Bienert GP, Chaumont F. Aquaporin-facilitated transmembrane diffusion of hydrogen peroxide. Biochim Biophys Acta - Gen Subj. 1840;2014:1596–604.

    Google Scholar 

  75. Benabdellah K, Ruiz-Lozano JM, Aroca R. Hydrogen peroxide effects on root hydraulic properties and plasma membrane aquaporin regulation in Phaseolus vulgaris. Plant Mol Biol. 2009;70:647–61. Springer Netherlands

    Article  CAS  PubMed  Google Scholar 

  76. Aroca R, Amodeo G, Fernández-Illescas S, Herman EM, Chaumont F, Chrispeels MJ. The role of aquaporins and membrane damage in chilling and hydrogen peroxide induced changes in the hydraulic conductance of maize roots. Plant Physiol. 2005;137:341–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  77. Lee SH, Chung GC, Steudle E. Gating of aquaporins by low temperature in roots of chilling-sensitive cucumber and chilling-tolerant figleaf gourd. J Exp Bot. 2005;56:985–95.

    Article  CAS  PubMed  Google Scholar 

  78. Boursiac Y, Boudet J, Postaire O, Luu D-T, Tournaire-Roux C, Maurel C. Stimulus-induced downregulation of root water transport involves reactive oxygen species-activated cell signalling and plasma membrane intrinsic protein internalization. Plant J. 2008;56:207–18.

    Article  CAS  PubMed  Google Scholar 

  79. Tian S, Wang X, Li P, Wang H, Ji H, Xie J, et al. Plant aquaporin AtPIP1;4 links Apoplastic H 2 O 2 induction to disease immunity pathways. Plant Physiol. 2016;171:1635–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  80. Schmidt R, Kunkowska AB, Schippers JHM. Role of reactive oxygen species during cell expansion in leaves. Plant Physiol. 2016;172:2098–106. American Society of Plant Biologists

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  81. Rae L, Lao NT, Kavanagh TA. Regulation of multiple aquaporin genes in Arabidopsis by a pair of recently duplicated DREB transcription factors. Planta. 2011;234:429–44.

    Article  CAS  PubMed  Google Scholar 

  82. Zhu D, Wu Z, Cao G, Li J, Wei J, Tsuge T, et al. TRANSLUCENT GREEN, an ERF family transcription factor, controls water balance in Arabidopsis by activating the expression of aquaporin genes. Mol Plant. 2014;7:601–15.

    Article  CAS  PubMed  Google Scholar 

  83. Hichri I, Muhovski Y, Clippe A, Žižková E, Dobrev PI, Motyka V, et al. SlDREB2, a tomato dehydration-responsive element-binding 2 transcription factor, mediates salt stress tolerance in tomato and Arabidopsis. Plant Cell Environ. 2016;39:62–79.

    Article  CAS  PubMed  Google Scholar 

  84. Liao X, Guo X, Wang Q, Wang Y, Zhao D, Yao L, et al. Overexpression of MsDREB6.2 results in cytokinin-deficient developmental phenotypes and enhances drought tolerance in transgenic apple plants. Plant J. 2017;89:510–26.

    Article  CAS  PubMed  Google Scholar 

  85. Kim MJ, Ruzicka D, Shin R, Schachtman DP. The Arabidopsis AP2/ERF transcription factor RAP2.11 modulates plant response to low-potassium conditions. Mol Plant. 2012;5:1042–57.

    Article  CAS  PubMed  Google Scholar 

  86. González-Morales SI, Chávez-Montes RA, Hayano-Kanashiro C, Alejo-Jacuinde G, Rico-Cambron TY, de Folter S, et al. Regulatory network analysis reveals novel regulators of seed desiccation tolerance in Arabidopsis thaliana. Proc Natl Acad Sci U S A. 2016;113:E5232–41.

    Article  PubMed  PubMed Central  Google Scholar 

  87. Borg M, Brownfield L, Khatab H, Sidorova A, Lingaya M, Twell D. The R2R3 MYB transcription factor DUO1 activates a male germline-specific regulon essential for sperm cell differentiation in Arabidopsis. Plant Cell. 2011;23:534–49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  88. Nicolas P, Lecourieux D, Kappel C, Cluzet S, Cramer G, Delrot S, et al. The basic leucine zipper transcription factor ABSCISIC ACID RESPONSE ELEMENT-BINDING FACTOR2 is an important transcriptional regulator of abscisic acid-dependent grape berry ripening processes. Plant Physiol. 2014;164:365–83.

    Article  CAS  PubMed  Google Scholar 

  89. Ma N, Xue J, Li Y, Liu X, Dai F, Jia W, et al. Rh-PIP2;1, a rose aquaporin gene, is involved in ethylene-regulated petal expansion. Plant Physiol. 2008;148:894–907. American Society of Plant Biologists

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  90. Olaetxea M, Mora V, Bacaicoa E, Baigorri R, Garnica M, Fuentes M, et al. ABA-regulation of root hydraulic conductivity and aquaporin gene- expression is crucial to the plant shoot rise caused by rhizosphere humic acids. Plant Physiol. 2015.

  91. Jang JY, Kim DG, Kim YO, Kim JS, Kang H. An expression analysis of a gene family encoding plasma membrane aquaporins in response to abiotic stresses in Arabidopsis thaliana. Plant Mol Biol. 2004;54:713–25.

    Article  CAS  PubMed  Google Scholar 

  92. Aroca R, Ferrante A, Vernieri P, Chrispeels MJ. Drought, abscisic acid and transpiration rate effects on the regulation of PIP aquaporin gene expression and abundance in Phaseolus vulgaris plants. Ann Bot. 2006;98:1301–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  93. Morillon R, Chrispeels MJ. The role of ABA and the transpiration stream in the regulation of the osmotic water permeability of leaf cells. Proc Natl Acad Sci U S A. National Academy of Sciences; 2001;98:14138–14143.

  94. Grondin A, Rodrigues O, Verdoucq L, Merlot S, Leonhardt N, Maurel C. Aquaporins contribute to ABA-triggered stomatal closure through OST1-mediated phosphorylation. Plant Cell Online. 2015;

  95. Rattanakon S, Ghan R, Gambetta GA, Deluc LG, Schlauch KA, Cramer GR. Abscisic acid transcriptomic signaling varies with grapevine organ. BMC Plant Biol. 2016;16:72.

    Article  PubMed  PubMed Central  Google Scholar 

  96. Wu G, Lin W-C, Huang T, Poethig RS, Springer PS, Kerstetter RA. KANADI1 regulates adaxial-abaxial polarity in Arabidopsis by directly repressing the transcription of ASYMMETRIC LEAVES2. Proc Natl Acad Sci U S A. 2008;105:16392–7. National Academy of Sciences

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors would like to thank Grant Cramer for his help and technical advice and the grapevine research community for making various RNA-seq and microarray datasets publicly available.

Funding

This study has been carried out with financial support from the University of Bordeaux’s Initiative of Excellence (IdEx) program, doctoral school of life and health sciences, and Cluster of Excellence COTE (ANR-10-LABX-45, within the Water Stress project), as well as the Canada Research Chairs Program, Genome British Columbia (10R21188), and the Natural Sciences and Engineering Research Council of Canada (10R23082). The funding bodies themselves were not involved in the design of the study or the collection, analysis, and interpretation of data, or in writing the manuscript.

Availability of data and materials

All data analyzed in this study are available from NCBI, https://www.ncbi.nlm.nih.gov/. The grapevine MIP gene family sequences were retrieved from the ORCAE 12× grapevine annotation V2, http://bioinformatics.psb.ugent.be/orcae/. All specific gene/protein accessions and transcriptomic dataset accessions can be found in Additional file 1: Table S1 and Additional file 1: Table S11, respectively.

Author information

Authors and Affiliations

Authors

Contributions

DCJW, SC, GG designed the experiments; DCJW, IM, SC, and GG carried out the experiments and analyses; DCJW, IM, and GG analyzed the results; and DCJW, LZ, SC, and GG wrote the manuscript. All authors participated in editing and approving the final work.

Corresponding author

Correspondence to Gregory A. Gambetta.

Ethics declarations

Ethics approval and consent to participate

Not applicable; this study has not directly involved humans, animals or plants.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional file

Additional file 1:

Supplementary Tables. (XLSX 748 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wong, D.C.J., Zhang, L., Merlin, I. et al. Structure and transcriptional regulation of the major intrinsic protein gene family in grapevine. BMC Genomics 19, 248 (2018). https://0-doi-org.brum.beds.ac.uk/10.1186/s12864-018-4638-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12864-018-4638-5

Keywords