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

Scion genotypes exert long distance control over rootstock transcriptome responses to low phosphate in grafted grapevine

Abstract

Background

Grafting is widely used in horticulture and rootstocks are known to modify scion growth and adaptation to soil conditions. However, the role of scion genotype in regulating rootstock development and functioning has remained largely unexplored. In this study, reciprocal grafts of two grapevine genotypes were produced as well as the corresponding homo-graft controls. These plants were subjected to a low phosphate (LP) treatment and transcriptome profiling by RNA sequencing was done on root samples collected 27 h after the onset of the LP treatment.

Results

A set of transcripts responsive to the LP treatment in all scion/rootstock combinations was identified. Gene expression patterns associated with genetic variation in response to LP were identified by comparing the response of the two homo-grafts. In addition, the scion was shown to modify root transcriptome responses to LP in a rootstock dependent manner. A weighted gene co-expression network analysis identified modules of correlated genes; the analysis of the association of these modules with the phosphate treatment, and the scion and rootstock genotype identified potential hub genes.

Conclusions

This study provides insights into the response of grafted grapevine to phosphate supply and identifies potential shoot-to-root signals that could vary between different grapevine genotypes.

Background

Phosphorus (P) is an essential macronutrient for plant development involved in numerous metabolic and signalling pathways [1]. Plants are able to take up P from the soil only under its free inorganic form, orthophosphate (Pi), although many forms of P unavailable to plant uptake are present in the soil [2]. Plants have evolved several mechanisms to maximize the acquisition of Pi when resources are limiting; mechanisms ranging from modifications of plant morphology to altering gene expression. For example, many plants are able to increase the Pi foraging capacity of the root system by allocating relatively more carbon to root than shoot growth, and consequently altering shoot/root biomass ratio, and by increasing root branching and production of fine roots [3]. Plants can also increase the bioavailability of P in the soil by exuding protons, organic acids, acid phosphatases and ribonucleases to release Pi from inaccessible sources [2]. These changes are accompanied by metabolic changes in both shoots and roots, in general, the concentration of sugars, organic acids, nitrogenous compounds and secondary metabolites increase, whereas the concentration of phosphorylated metabolites decreases [4,5,6]. In addition, internal Pi use can be economized by the replacement of phospholipids of cellular membranes by sulfolipids and galactolipids [7].

In response to limited Pi supply, a number of genes involved in sensing, signalling and responses to Pi starvation have been identified and are called Pi-starvation-inducible (PSI) genes [8]. Lan et al. (2015) performed a meta-analysis of the response of Arabidopsis roots to P supply and identified a core of one hundred PSI genes induced hours or days after the onset of Pi starvation, independently of growth condition, experimental design and method of transcriptome analysis [8]. This list of PSI genes included genes related to lipid metabolism and galactolipid biosynthesis, transcription factors (TFs) containing SYG/PHO81/XPR1 (SPX) domains, Pi transporters (such as those from the PHOSPHATE TRANSPORTER1/PHT1 family), protein kinases and intracellular secreted purple acid phosphatases [8]. The transcriptome response of perennial plants to Pi availability has received little attention despite both short and long term transcriptional responses being widely studied in annuals such as Arabidopsis, maize, rice, white lupin, tomato, bean and wild mustard [9]. To date, the only studies done on perennial crops (Poncirus trifoliate, Pinus massoniana and poplar) have been restricted to the study of only one genotype and its adaptation to a long-term (weeks to months) limiting Pi treatments [10,11,12]).

Comparing the responses to limited Pi supply of different genotypes can provide evidence of genetic variation of adaptive strategies and identify targets for crop improvement [13,14,15]. In grapevine (Vitis spp.), genotypes differ in their growth responses to limited Pi availability [16, 17], and in Pi acquisition and utilization efficiencies [17, 18]. Approximately one third of total plant P is exported from vineyards annually (in the berries harvested and winter pruning [19]), which suggests that grapevine is an important crop to understand the molecular basis of the genetic variation of adaptive responses to Pi supply. Furthermore, in plants amenable to grafting such as grapevine, the responses of different scion/rootstock combinations to limited Pi supply can elucidate shoot and root-specific responses to Pi supply and provide insights into the shoot and/or root signals involved. Studies using grafting on model species such as Arabidopsis have shown that both root and shoot play roles in regulating plant responses to Pi supply; potential long-distance signals include ions, metabolites, proteins, mRNAs, small RNAs and hormones [9, 20]. Some root-to-shoot signals have been identified such as Pi itself (and potentially inositol polyphosphatases and pyrophosphatases) and hormones such as strigolactones and cytokinins [9, 20]. Shoot-to-root signals include sugars, microRNAs and a calcium related signal [9, 20]. Grafting experiments have shown that under limited Pi microRNA 399 can move from the shoot to the root to suppress PHOSPHATE2/PHO2 expression [21, 22]. PHOSPHATE2 mediates the protein degradation of phosphate transporters and controls shoot Pi homeostasis [23, 24]. Similarly, grafting experiments have shown that two vacuolar calcium transporters, CAX1 and CAX3, trigger a shoot-to-root signal that increases Pi uptake, but the exact mechanism remains unknown [25].

The molecular mechanisms of adaptation to low Pi supply have not yet been studied in grafted plants; here we investigate the roles of the scion and the rootstock in regulating root transcriptomic responses to low Pi supply in grafted grapevine. Firstly, we aimed to determine whether a core set of low Pi responsive transcripts could be identified across all scion/rootstock combinations studied. Secondly, we investigated whether there was genetic variation in the root transcriptomic response to low Pi supply. Thirdly, we asked whether the scion genotype could modify rootstock gene expression and its response to limited Pi supply. A weighted gene co-expression network analysis (WGCNA) allowed the identification of co-expression gene modules strongly correlated to P supply, and the rootstock and scion genotype, and identified potential hub genes.

Results

Four scion/rootstock combinations of grapevine (1103P/1103P, 1103P/PN, PN/1103P and PN/PN) were grown in aerated hydroponic culture on HP, and then half of them were transferred to LP while the remaining half continued to receive HP. After 27 h treatment, RNA-seq was used to compare the short-term transcriptome response to LP in the different scion/rootstock combinations (described in detail below). And after 28 d of treatment, RT-qPCR was used to quantify the expression of four PSI genes (two SPX domain containing genes and two PHT1 transporters) in the roots tips of plants grown on LP and HP; the expression of these genes was strongly up-regulated by LP in all scion/rootstock combinations (Fig. 1).

Fig. 1
figure 1

Expression of selected phosphate starvation induced genes in the root tips of grafted grapevine after 28 d of either high or low phosphate treatment (dark or light bars respectively). Gene expression was quantified by RT-qPCR and expressed as normalized relative quantities (NRQs) in reciprocal scion/rootstock combinations of Vitis vinifera cv. Pinot noir (PN) and V. rupestris x V. berlandieri cv. 1103 Paulsen (1103P). Means, standard deviations and results of a two-way analysis of variance (p < 0.05, with Tukey’s Honest Significant Difference test) are shown

A set of genes was up- or down-regulated in response to LP in all scion/rootstock combinations

To characterize short-term gene expression changes in response to LP, RNA-Seq was used to quantify mRNA abundance in the root tips 27 h after transfer to LP in all scion/rootstock combinations. In comparison to the HP, the LP treatment resulted in the up- and down-regulation of the expression of 1834 and 2501 genes respectively across all scion/rootstock combinations (Fig. 2a).

Fig. 2
figure 2

Venn diagram comparison of the number of shared and unique differentially expressed genes (log fold change > 1, false discovery rate adjusted p-value < 0.01) that a.) responded to 27 h low phosphate treatment in different scion/rootstock combinations of grafted grapevine, b.) were differentially expressed between the different the rootstock genotypes studied, and c.) responded to grafting with a non-self-scion. Root gene expression was studied in reciprocal scion/rootstock combinations of Vitis vinifera cv. Pinot noir (PN) and V. rupestris x V. berlandieri cv. 1103 Paulsen (1103P)

A set of 301 differentially expressed (DE) genes down-regulated in response to LP was identified across all scion/rootstock combinations (Additional File 1), which was enriched in the MapMan BINs WRKY and ETHYLENE RESPONSE FACTOR/ERF TFs, members of the TYROSINE KINASE-LIKE (TKL) superfamily, PHENYLALANINE AMMONIA-LYASES/PALs and naringenin-chalcone synthases (particularly STILBENE SYNTHASES/STSs, Table 1). The WRKY TFs in this list included WRKY6 (Vitvi10g00063), WRKY18 (Vitvi04g00510), WRKY33 (Vitvi06g00741 and Vitvi08g00793), WRKY40 (Vitvi04g00511 and Vitvi09g01122), WRKY51 (Vitvi04g00760 and Vitvi04g01985), WRKY53 (Vitvi15g01003 and Vitvi16g01213), and WRKY75 (Vitvi01g01680). There were 16 ERF TFs in this list as well as other genes associated with ethylene signalling such as 1-AMINOCYCLOPROPANE-1-CARBOXYLATE OXIDASE 1/ETHYLENE-FORMING ENZYME (Vitvi12g00445) and the ethylene responsive gene, INFLORESCENCE DEFICIENT IN ABSCISSION (Vitvi01g01924). In addition, this list included genes related to secondary metabolism, such as PREPHENATE DEHYDRATASE (Vitvi06g01946), which is involved in the biosynthesis of aromatic amino acids such as phenylalanine, and two CINNAMATE-4-HYDROXYLASES/C4Hs (Vitvi11g00924 and Vitvi11g01045) and two 4-COUMAROYL COA-LIGASES/4CLs (Vitvi14g01757 and Vitvi02g00717) which are involved in phenylpropanoid metabolism.

Table 1 MapMan BINs enriched in the 301 genes down-regulated in the roots in response to 27 h of a low phosphate treatment in the four scion/rootstock combinations of grafted grapevine studied (PN/PN, PN/1103P, 1103P/PN and 1103P/1103P). PN: Vitis vinifera cv. Pinot noir; 1103P: V. rupestris x V. berlandieri cv. 1103 Paulsen

Only 8 DE genes were up-regulated in response to LP in all the scion/rootstock combinations, these genes were: an ABC transporter (Vitvi14g01865), a receptor kinase (FERONIA, Vitvi14g01365), a disease resistance protein (Vitvi18g01755), two transferases (Vitvi12g02284 and Vitvi07g03041), ISOFLAVONE 2′-HYDROXYLASE/CYP81D8 (Vitvi07g01651) and a putative transcriptional regulator (the transposase-like protein Vitvi13g02203) (Additional File 1). In summary, the DE genes down-regulated in response to LP were well conserved between the scion/rootstock genotypes studied, but the DE genes up-regulated in response to LP were largely scion/rootstock combination specific.

Genes differentially regulated between the two homo-grafts in response to LP

Many genes were DE between the two rootstock genotypes under both HP and LP (Fig. 2b, Additional File 2). Under the control HP conditions, the genes DE between the two rootstock genotypes were enriched in the MapMan BINs related to enzymes and nucleotide-binding domain and leucine-rich repeat-containing (NLR) effector receptors (which perceive pathogen infection to trigger immunity mechanisms). In addition, the MapMan BINs PAL and chalcone synthase were also enriched in the genes more highly expressed in PN (Additional File 3).

The comparison of the transcriptome response to LP between the two homo-grafts (PN/PN and 1103P/1103P) clearly demonstrated that more genes were DE in PN/PN than in 1103P/1103P (Fig. 2a). To determine whether the transcriptome response to LP was different between the homo-grafts, interactions were calculated (i.e. the genes DE between (PN/PN LP versus PN/PN HP) versus (1103P/1103P LP versus 1103P/1103P HP), Additional File 4). The expression of 310 genes differently responded to LP between the two homo-grafts (Additional File 4), these 310 genes were enriched in the MapMan BINs NAC and WRKY TFs, solute transport (VACUOLAR IRON TRANSPORTER/VIT and ALUMINIUM ACTIVATED MALATE TRANSPORTER families), biotic stress, oxidoreductases (mainly orthologues of LACCASE 14, which is involved in lignin degradation) and transferases (that transfer a P containing group) (Additional File 5). The BINs from the biotic stress category were NLR effector receptors and the ERN1 transcription factor (which regulates rhizobial infection in Lotus japonicas [26]). To characterize the expression patterns of these 310 genes showing a significant interaction in response to LP between the homo-grafts, the log fold change of gene expression in response to LP in 1103P/1103P was plotted against PN/PN (Fig. 3a). Globally, the DE genes down-regulated in response to LP were qualitatively well-conserved between the two homo-grafts, with 67 DE genes significantly down-regulated in both PN/PN and 1103/1103P, but the fold change was higher for PN/PN. The DE genes more highly down-regulated in response to LP in PN/PN than 1103P/1103P included many TFs, kinases, biotic stress associated genes and genes related to calcium signalling. The DE genes up-regulated in response to LP were less conserved between the homo-grafts with only one gene (a methyl jasmonate esterase, Vitvi07g02517) significantly up-regulated in response to LP in both homo-grafts, but with a higher fold change in PN/PN. The DE genes more highly up-regulated in response to LP in PN/PN than 1103P/1103P included WRKY19 (Vitvi18g01743), a sulphate transporter (Vitvi09g00484), the Snf1 protein kinase KIN10 (Vitvi08g00935) and a large subunit of ADP-glucose pyrophosphorylase (Vitvi07g00544) which regulates starch synthesis as well as genes associated with jasmonate, cytokinins and abscisic acid (ABA) signalling. Only one gene (Vitvi08g01380) was up-regulated in response to LP in 1103P/1103P and down-regulated in PN/PN, and vice versa only one gene (Vitvi18g02447) was down-regulated in response to LP in 1103P/1103P and up-regulated in PN/PN. In conclusion, the transcriptome response to LP of PN/PN was more extensive than that of 1103P/1103P, but the qualitative change in gene expression is similar between the two homo-grafts, particularly for the down-regulated genes.

Fig. 3
figure 3

Relationship between a.) the transcriptional responses to 27 h low phosphate treatment in PN/PN and 1103P/1103P of the genes showing a significant interaction (log fold change > 1, false discovery rate adjusted p-value < 0.01) between homograft genotype and phosphate treatment, and b.) the transcriptional responses to 27 h low phosphate treatment in PN/1103P and 1103P/1103P of the genes showing a significant interaction between the scion genotype and the phosphate treatment. The scion/rootstock combinations studied were grafts of PN (Vitis vinifera cv. Pinot noir) and 1103P (V. rupestris x V. berlandieri cv. 1103 Paulsen)

Genes differentially regulated in response to grafting with a non-self-scion

Grafting with a non-self-scion triggered the differential expression of many genes; however one gene, a BURP domain-containing protein RD22-like isoform X2, VViBURP17 (Vitvi04g00357) [27] was DE in response to hetero-grafting in both rootstocks and P treatments (Fig. 2c, Additional File 6). BURP domain containing proteins are involved in abiotic stress tolerance in several plant species and are generally responsive to ABA application [27].

When comparing 1103P/PN and PN/PN, only a small number of DE genes responded to grafting with a non-self-scion in both LP and HP in the PN rootstock; 6 and 24 DE genes were up and down-regulated respectively (Fig. 2c). The six transcripts that increased in abundance in the 1103P/PN relative to PN/PN included CAROTENOID CLEAVAGE DIOXYGENASE 8 (Vitvi04g00298), an enzyme involved in the synthesis of strigolactones (Additional File 7). The 24 transcripts that decreased in abundance were enriched in the MapMan BINs meiotic exit regulator (involved in cell cycle), protein modification and iron uptake (Additional File 7).

When comparing PN/1103P and 1103P/1103P, a large transcriptional response to the grafting with a PN scion was detected with 661 and 505 DE genes up- and down-regulated respectively under both LP and HP (Fig. 2c, Additional File 8). The 661 DE genes that were up-regulated in PN/1103P compared to 1103P/1103P were enriched in the MapMan BINs DNA-binding with one finger, NAC, basic helix-loop-helix/bHLH and homeobox TFs, TKL kinases, subtilisin-like proteases (which are associated with biotic stress responses), proteins involved in cellulose synthesis, heat shock proteins and galactinol synthase (Additional File 8). Galactinol synthase mediates raffinose accumulation, an osmoprotectant that promotes stress tolerance [28]. The 505 DE genes that were down-regulated in PN/1103P compared to 1103P/1103P were enriched in the MapMan BINs GATA TFs, bHLH TFs, GROWTH REGULATING FACTORS, serine carboxypeptidases, hydrolases and precursors of cysteine rich peptides (Additional File 8). Cysteine rich peptides are involved in regulating a range of developmental responses in plants and can be involved in cell-to-cell communication, particularly in response to the establishment of plant-bacteria symbiosis and root development [29].

Overall, the transcriptome of the rootstock responds to grafting with a non-self-scion, but the degree of transcriptome response is scion/rootstock dependent and there does not appear to be a general non-self-scion response of the rootstock transcriptome.

The impact of the scion on the root transcriptome response to LP

The comparison of the response to LP of PN/PN and 1103P/PN showed that the scion can have a slight impact on the transcriptome response of the rootstock to LP (Fig. 2a). However, the response to LP was quantitatively different for only five genes between PN/PN and 1103P/PN (i.e. genes DE between (1103P/PN LP versus 1103P/PN HP) versus (PN/PN LP versus PN/PN HP), Additional File 9). The four genes were more down-regulated in 1103P/PN than PN/PN were two ABCG40 transporters (Vitvi09g00473 and Vitvi09g00478), a gene involved in tyrosine synthesis (Vitvi12g00272) and an UDP-GLYCOSYLTRANSFERASE encoding gene (Vitvi12g01724) assigned to the MapMan BIN ABA synthesis/degradation.

In contrast, the comparison of the response to LP of 1103P/1103P and PN/1103P showed that the PN scion dramatically increased the transcriptome response of the 1103P rootstock to LP; only 59 and 42 transcripts specifically decreased and increased in abundance in 1103P/1103P, whereas 1181 and 1084 transcripts decreased and increased in abundance in PN/1103P respectively. To determine whether the transcriptome response to LP was quantitatively different between the 1103P/1103P and PN/1103P interactions were calculated (i.e. the genes DE between (PN/1103P LP versus PN/1103P HP) versus (1103P/1103P LP versus 1103P/1103P HP), Additional File 10). Only 364 DE genes differently responded to LP, these genes were enriched in the MapMan BINs ERF TFs, NLR effector receptors and phosphoserine aminotransferases (which are involved in serine biosynthesis) (Table 2). As was seen in the comparison between the transcriptional responses to LP of the homo-grafts, there was a conservation of response in the LP down-regulated genes (Fig. 3b), with 113 DE genes significantly down-regulated in PN/1103P and 1103P/1103P, but the gene expression response of PN/1103P was of higher magnitude. However, no common DE genes were significantly up-regulated in response to LP in the genes showing an interaction (Fig. 3b); the DE of 12 genes was reduced in response to LP in PN/1103P and increased in response to LP in 1103P/1103P.

Table 2 MapMan BINs enriched in the 364 genes differentially expressed in the roots between the scion/rootstock combinations PN/1103P and 1103P/1103P in response to 27 h of a low phosphate treatment. PN: Vitis vinifera cv. Pinot noir; 1103P: V. rupestris x V. berlandieri cv. 1103 Paulsen

The 104 DE genes more up-regulated in PN/1103P than 1103P/1103P included a cyclin (Vitvi15g00641), the bHLH TF FIZZY-RELATED 3 (Vitvi12g02634), an ALUMINUM ACTIVATED MALATE TRANSPORTER (Vitvi13g02145), four DICER-LIKE3 genes, FAR-RED IMPAIRED RESPONSE1-RELATED SEQUENCE5 (Vitvi08g01733) and many resistance proteins and genes involved in biotic stress responses and (Additional File 10).

The 260 DE genes more down-regulated in PN/1103P than 1103P/1103P included genes related to amino acid metabolism, development (including a VIT family protein), fermentation and glycolysis (cytosolic branch e.g. PYRUVATE KINASE (Vitvi10g01568), PHOSPHOGLYCERATE KINASE (Vitvi19g01690) and GLYCERALDEHYDE-3-PHOSPHATE DEHYDROGENASE (Vitvi01g01538)) (Additional File 10). This list also included five ERF TFs, a zinc finger homeobox TF (Vitvi12g02000), a NADP-DEPENDENT MALIC ENZYME (Vitvi04g00009), two EXORDIUM-like 2 encoding genes, and a TREHALOSE-6-PHOSPHATE PHOSPHATASE (Vitvi18g00384) (Additional File 10).

In summary, the scion can modify rootstock transcriptome responses to LP, but its effect is rootstock dependent.

Gene co-expression network analysis

A WGCNA approach was used to identify 26 modules of highly correlated genes and these modules were analysed for their association with each experimental variable (the P treatment, and scion and rootstock genotype) (Fig. 4). Genes with highly correlated expression patterns were identified based on kME values (Additional File 11) those with the highest correlation coefficients (> 0.9) could be considered as potential hub genes within the module. The module MEgrey contained the genes which could not be assigned to a co-expression module.

Fig. 4
figure 4

Module-trait relationships from the weighted gene co-expression network analysis. Module-phosphate treatment/scion/rootstock genotype correlation coefficients and corresponding p-values (in parenthesis) are given in each cell. The left panel shows the 26 modules and the number of module member genes. The colour scale on right shows module-trait correlation from −1 (blue) to 1 (red). HP: high phosphate, LP: low phosphate; PN: Vitis vinifera cv. Pinot noir; 1103P: V. rupestris x V. berlandieri cv. 1103 Paulsen (1103P)

Gene co-expression modules related to P supply

The gene co-expression module the most positively correlated to LP was MEyellow; which contained 1704 genes (Fig. 4). The highly correlated genes (> 0.8) were enriched in the MapMan BINs transcriptional activation, particularly C3H zing finger and FAR-RED IMPAIRED RESPONSE 1 TFs (Table 3). The most 144 highly correlated genes (> 0.9) included genes involved in RNA (particularly small RNA) processing (such as DICER-LIKE1 and DICER-LIKE3), chromatin remodelling factors, two histone acetyltransferases (Vitvi12g00328 and Vitvi19g00124), DEMETER (Vitvi08g01515), three genes from the MapMan BIN chromatin structure as well as other biotic and abiotic stress responsive transcripts (Additional File 12). The co-expression module most positively correlated with HP was MElightgreen, which contained 438 transcripts. The positively correlated transcripts (> 0.8) were enriched in the MapMan BINs WRKY and ERF TFs, and chalcone synthase (Table 3).

Table 3 MapMan BINs enriched in the most highly positively correlated (> 0.8) genes in selected modules of the weighted gene co-expression network analysis made from the genes expressed in the roots of different scion/rootstock combinations of grapevine grown under two different phosphate supplies

Gene co-expression modules related to the transcriptome differences between the rootstock genotypes

The co-expression module that was most highly positively correlated with the rootstock 1103P is MEturquoise; this module contained 2455 genes (Fig. 4). The highly correlated genes (> 0.8) were enriched in the MapMan BINS associated with NLR effector receptors and various enzymes (Table 3). The co-expression module that was most highly positively correlated with the rootstock PN is MEgreen, containing 1590 genes and the highly correlated genes (> 0.8) were enriched in the MapMan BINs involved with RNA processing, NLR effector receptors and various enzymes (Table 3).

Gene co-expression modules related to the effect of the scion genotype on the root transcriptome

Fewer gene co-expression modules were significantly correlated with the scion genotype and the correlation coefficients were much lower than for the P treatment or rootstock genotype (Fig. 4). The module MElightyellow was the most highly correlated with the 1103P scion (298 transcripts) but was also significantly correlated with the rootstock genotype. Whereas the module MEgreenyellow was only significantly correlated with the scion genotype so the genes are more likely to be scion response specific. MEgreenyellow was enriched in the MapMan BINs signalling peptides (cysteine rich peptides), transcriptional activation, protein degradation (serine carboxypeptidases) and cell wall (Table 3). The 121 most highly positively correlated transcripts (> 0.9) included five GROWTH REGULATING FACTORS, which are potentially targeted by the microRNA 396 [30] (Additional File 12). The module MEmagenta was the co-expression module most highly correlated with the PN scion (736 transcripts), the MapMan BINs strigolactone signalling, C2C2, NAC, bHLH and BEL TFs, and some TKL receptor kinases were enriched in this module (Table 3).

Discussion

The availability of a reference genome for grapevine is a great resource for RNAseq analysis, but our knowledge is limited to the predicted transcriptome of this reference genome (i.e. PN40024). The literature suggests that genomic variations are significant even at the clonal level [31]. Sequencing the clone of 1103P that was used in this study would help correct for transcript quantification of highly duplicated genes/gene families, but was beyond the scope of this project. However, despite differences between genomes of the same genus, pan-genome studies have showed that there is a high proportion of orthologous genes shared between related species [32, 33].

Many PSI genes from Arabidopsis were not induced in response to 27 h of LP in grapevine

In this study on grapevine, we identified a set of 301 and eight DE genes that were down and up-regulated respectively in all scion/rootstock combinations. Surprisingly, there were few of the PSI genes of Arabidopsis [8] in this dataset, this is partly because the grapevine orthologues have not been identified (e.g. INDUCED BY PHOSPHATE STARVATION1IPS1 and IPS2) or they were expressed at very low levels (e.g. MITOGEN-ACTIVATED PROTEIN KINASE KINASE KINASE 19). However, some orthologues of Arabidopsis PSI genes were not DE, e.g. PHOSPHATE STARVATION RESPONSE1-LIKE1 (Vitvi07g00666 and Vitvi14g00736) and the two orthologues of PHO2 (Vitvi00g02237 and Vitvi10g01764, which were not significantly DE, but were slightly up-regulated under LP particularly in PN/PN and PN/1103P). Orthologues of WRKY6, WRKY18, WRKY40 and WRKY75 TFs were present in the set of 301 genes that were down-regulated in response to LP in all scion/rootstock combinations, although these genes are PSI genes in Arabidopsis and have roles in regulating Pi starvation responses [34,35,36]. However, after 28 d of LP, the expression of orthologues of four Arabidopsis PSI genes was up-regulated relative to HP suggesting that the response to LP at 27 h represents the early stress-related changes in gene expression.

A set of 309 genes that were DE in response to LP in all scion/rootstock genotypes was identified

The down-regulation of gene expression in response to LP was well-conserved between the four scion/rootstock combinations studied (301 genes), whereas the expression of only eight genes was up-regulated in response to LP across all scion/rootstock combinations. Messenger RNA stability plays an important role in the regulation of gene expression. Short-lived transcripts (such as TFs and kinases) are often associated with regulatory processes (whereas longer-lived transcripts tend to be associated with protein synthesis and energy balance) [37, 38]. It has previously been shown that transcript degradation makes a large contribution to the rapid response to sucrose in Arabidopsis seedlings [39] and is involved in a wide range of abiotic stress responses (as reviewed by [40]. The set of 301 genes down-regulated in response to LP in all scion/rootstock combinations was enriched in TFs and kinases suggesting that their down-regulation could be related to their rapid turnover. It seems that conservation of genes down-regulated in response to LP between the different scion/rootstock combinations is important and this could be because rapid transcript turnover has the potential to ensure quick responses to changing soil nutrient conditions.

Many genes involved in the pathway of stilbene synthesis were present in the set of 301 DE genes down-regulated in response to LP in all scion/rootstock combinations: for example one PREPHENATE DEHYDRATASE, nine PALs, two C4Hs, two 4CLs and 13 STSs. Many metabolites are known to be differentially accumulated in plants in response to LP, a general response to increase the concentration of secondary metabolites (particularly the concentration of anthocyanins in leaves in response to LP supply [41]), but little is known about whether stilbene synthesis is altered in response to LP because only few plants produce stilbenes. To date, there has only been one study of a plant that produces stilbenes, this study was done on pine trees, after a longer term low Pi treatment and did not highlight the differential expression of genes related to secondary metabolism [11]. Stilbenes are known to have roles in plant defence responses [42], so this down-regulation of expression of genes related to stilbene synthesis 27 h after transfer to a LP treatment may be an early response to redirect metabolism from defence responses to adapting to limited P. Many ethylene related genes were also present in the set of 301 DE genes that were down-regulated in response to LP in all scion/rootstock combinations; these genes are known to have a wide range of roles in both biotic and abiotic stress responses, including regulating root growth responses to P supply (as reviewed by [9]). In the WGCNA analysis, the module most strongly positively correlated with HP was MElightgreen; the highly correlated genes in this module were enriched in the MapMan BINs WRKY and ERF TFs and chalcone (stilbene) synthases. Both WRKY TFs and ethylene have been shown to regulate the expression of STSs in grapevine [43] suggesting that these TFs could be hub genes regulating metabolomic responses to LP. In addition, VviMYB15 (Vitvi05g01733) and VViMYB14 (Vitvi07g00598), positive regulators of STSs [44], were down-regulated in most scion/rootstock combinations. The down-regulation of expression of VviMYB14 and many STSs has also been observed in grapevine in response to iron deficiency, these changes were accompanied by an increase in the expression of genes involved in flavonoid synthesis [45]. Similarly, one of the eight DE genes that were up-regulated in all scion/rootstock combinations in response to LP was ISOFLAVONE 2′-HYDROXYLASE/CYP81D8 (Vitvi07g01651). This suggests that at the transcriptional level the down-regulation of stilbene synthesis and the up-regulation of flavonoid synthesis are common responses to Pi and iron deficiency in grapevine roots.

The expression of only eight genes was up-regulated in response to LP in all scion/rootstock combinations, this suggest that there was more variation in the up- than the down-regulation of gene expression in response to LP of individual scion/rootstock combinations. Despite this variation of transcriptome response, the WGCNA analysis allowed the identification of modules positively correlated with the LP treatment. The most strongly positively correlated module with LP was MEyellow; the genes in this module were enriched in the MapMan BINs transcriptional activation including FAR-RED IMPAIRED RESPONSE1. FAR-RED IMPAIRED RESPONSE1 is known to be involved in phytochrome signalling, but it has recently been shown to regulate the PHR1 expression (along with ethylene) in Arabidopsis [46]. This could suggest that FAR-RED IMPAIRED RESPONSE1 TFs regulate early LP transcriptome responses in grapevine. Many genes involved in chromatin modification and RNA processing such as DICER-LIKE1 and DICER-LIKE3 genes were present in the list of genes most positively correlated with MEyellow. In Arabidopsis, DICER-LIKE1 is involved in the formation of microRNAs, which are known to be differentially produced in response to Pi starvation [9]. DICER-LIKE3 generates siRNAs that direct DNA methylation and has been shown to be inhibited by Pi starvation [47]; the presence of DICER-LIKE3 genes in MEyellow suggests that siRNA-directed DNA methylation may also be important in the response of grapevine to P availability. DNA methylation is known to be important in the regulation of Pi starvation responses in both annual [48] and perennial species [49]. Chromatin modifications are implicated in the regulation of PSI gene expression in Arabidopsis [50] and one transposase-like protein was present in the set LP up-regulated genes in all scion/rootstock combinations suggesting that this may also be important in grapevine.

Do grapevine genotypes differ in their response to LP?

Many genes were DE in the roots between the two different genotypes studied under both HP and LP conditions, and the modules MEturquoise and MEgreen were associated with the rootstocks 1103P and PN respectively. The genes DE between the two rootstocks were associated with various enzymes and NLR effector receptors, similar results have been acquired in other genotype-specific transcriptome comparisons in grapevine [51, 52]. The two homo-grafts also differed in their transcriptome response to P supply; the transcriptome response of PN/PN was more extensive than that of 1103P/1103P after 27 h of LP treatment although the qualitative changes in gene expression were similar between the two genotypes. Genes related to iron, organic acid and sulphate transport, and carbon metabolism were DE between the two homo-grafts along with genes associated with hormone signalling. The large number of TFs was DE between the two homo-grafts in response to LP suggesting that the genotypes differ in the TF networks employed to regulate responses to LP.

Does the scion genotype modify the rootstock transcriptome and its response to P supply?

There have been a number of studies in the literature of how rootstocks alter gene expression in the scion of perennial [53,54,55,56,57] and annual crops [58], and the responses of the scion to different environmental conditions [59, 60]. However, few studies have been done the other way round, to investigate how the scion modifies the transcriptome of the rootstock and its response to the environment. In this study we identified only one transcript, VViBURP17, which increased in abundance in all hetero-grafts in comparison to the corresponding homo-graft control. This suggests that there no conserved root transcriptome response to grafting with a non-self-scion and that each rootstock responds differently.

In grapevine, grafting with non-self-rootstocks triggers the differential expression of genes involved in defence and stress responses (such as genes encoding the biosynthesis of secondary metabolites and receptor kinases), chromatin modification, transcriptional regulation and hormone signalling [53, 54] in the scion. Similar genes were DE in the root in response to grafting with a non-self-scion. In addition, some of the scion-responsive transcripts and potential hub genes identified in the WGCNA suggest that the scions used in this study alter cysteine rich peptide signalling and microRNA production (particularly microRNA 396). It is also possible that some of the scion-responsive transcripts are in fact mRNAs that have moved from the scion to the rootstock, such movement is complex to understand and in tobacco/tomato grafts seems to involve both a regulated and un-regulated components [61]. We assume that these defence and stress-related scion-responsive transcripts are not related to problems of incompatibility because grapevine grafts relatively easily in in vitro culture and the treatments were applied 10 weeks after grafting.

The impact of the scion on the rootstock response to LP was dependent on the rootstock studied; only six genes were DE between 1103P/PN and PN/PN in response to LP, whereas 364 genes were DE between PN/1103P and 1103P/1103P. This is largely because few genes were DE in response to LP in the 1103P/1103P in comparison to the other scion/rootstock combinations. In general, the qualitative change in gene expression of these 364 was similar between the two scion genotypes, particularly for the genes down-regulated in response to LP. The enrichment of MapMan BINs in the 364 genes DE in response to LP between PN/1103P and 1103P/1103P suggests that the PN scion particularly alters ethylene signalling and amino acid biosynthesis.

Conclusion

Here we showed that short-term (27 h) of LP treatment in grafted grapevine resulted in the differential regulation of many genes in the roots; a core set of genes DE in all scion/rootstock in response to LP was identified, which surprisingly contained only few of the PSI genes identified in Arabidopsis (although four of these genes were induced after 28 d of treatment). This core set of DE genes contained many more down- than up-regulated genes and many genes related to the synthesis of stilbenes were present in the list of down-regulated genes. The comparison of the response to LP between the two homo-grafts highlighted the genotype specific variation in transcriptional response. Superimposed upon the genotype specific variation in transcriptional response to LP, the scion also modified the rootstock gene expression response to LP in a rootstock genotype specific fashion. Weighted gene co-expression network analysis identified genes clusters correlated with P supply, and the rootstock and scion genotype. This data set suggests that some potential phloem mobile signals associated with P supply were differentially regulated in the different scion/rootstocks genotypes studied; these signals include metabolites, cysteine rich peptides, hormones such as strigolactones and ethylene, and microRNAs. Although these signals are known to convey messages from the shoot to the root in model species, here we show that genetic variation in a crops species can potentially alter these signalling events in grafted plants. This study highlights the importance of studying grafted plants during rootstock selection and that the scion genotype can modify root responses to the environment.

Methods

Plant material and growing conditions

An American rootstock genotype, the V. berlandieri x V. rupestris hybrid cv. ‘1103 Paulsen’ (1103P, clone 198; Vitis international variety catalogue number: 9023), was collected from a vineyard in Bordeaux, France (according to institutional guidelines). The identification of 1103P was done by the Centre de Ressources Biologiques de la Vigne, collection ampélographique de Vassal, Montpellier, France, by simple sequence repeats markers. The V. vinifera genotype used was obtained from the sexual reproduction from self-pollination of Pinot Noir, named 40,024 (PN), and was provided by INRAE Colmar, France, and its genome was sequenced by a French-Italian consortium [62] and it also grows and grafts well in in vitro culture. No permissions were required to obtain this plant material. Both genotypes were introduced into in vitro culture after surface sterilisation. All four possible scion/rootstock combinations were micro-grafted in vitro using the cleft grafting system, i.e. 1103P/1103P, 1103P/PN, PN/1103P and PN/PN. Plants were cultivated in vitro on McCown Woody Plant Medium (Duchefa) supplemented with 30 g L− 1 sucrose, 0.27 μM 1-naphthalene acetic acid and 0.4% agar in a growth chamber at 22 °C and with a photoperiod of 16 h light/8 h dark at a photon flux density of 55 μmol m− 2 s− 1. Six-week-old plantlets were then acclimated for 4 weeks to perlite-filled pots, irrigated with water, in a growth chamber at 26 °C and with a photoperiod of 16 h light/8 h dark at a photon flux density of 145 μmol m− 2 s− 1 at plant level. Plants were then transferred into hydroponic culture with an aerated solution; each pot contained 2 plants of the same scion/rootstock combination with 700 mL of high P solution (HP, 0.6 mM P). Four days later, half the pots continued receiving the HP solution and the other half were subjected to a low P (LP, 0.001 mM P) treatment. The macronutrient composition was 2.45 mM KNO3, 0.69 mM MgSO4 and 1.27 mM CaCl2 for both the HP and LP solutions; HP solution also contained 0.6 mM KH2PO4 and 0.6 mM CaSO4, whereas the LP solution contained 0.3 mM K2SO4 and 0.3 mM CaSO4. Micronutrients were supplied as 46 μmol H3BO3, 9.1 μmol MnCl2, 2.4 μmol ZnSO4, 0.5 μmol CuSO4 and 14 nmol (NH4)6Mo7O24, and iron was supplied as 8.5 mg L− 1 Sequestrene 138 (i.e. 31.3 μmol ethylenediamine-N,N′-bis (2-hydroxyphenylacetic acid) NaFe, Syngenta Agro S.A.S.). The solutions were changed once a week for 28 d.

RNA extraction

After 27 h and 28 d of LP or HP treatment, three pools of three root tips per plant (~ 15 mm in length) were harvested and immediately snap-frozen in liquid N for RNA sequencing (RNA-Seq) and quantitative reverse transcription polymerase chain reaction (RT-qPCR) analysis respectively. Total RNA of samples was extracted using the Spectrum Plant Total RNA kit (Sigma-Aldrich) with some modifications as previously described [54].

RT-qPCR analysis

Total RNA (1.5 μg) was reverse transcribed into cDNA using the SuperScript III First-Strand Synthesis System (Invitrogen). Quantitative reverse transcription polymerase chain reactions were performed using SYBR Green on an iCycler iQH (Bio-Rad) according to the procedure described by the supplier, with 0.2 μM of primers for each gene. Gene expression was calculated as normalized relative quantities [63] with three reference genes. Primer sequences are listed in Additional File 12. Statistical analysis was done using the software R [64]. Data were analysed using a two-way analysis of variance (p < 0.05, with Tukey’s Honest Significant Difference test).

RNA sequencing

RNA sequencing libraries were generated from 500 ng of total RNA using TruSeq Stranded mRNA LT Sample Preparation Kit (Illumina), according to manufacturer’s instructions. These libraries were then sequenced using Illumina Hiseq 4000 to produce 50 bp paired-end reads following Illumina’s instructions. The RNA-seq data has been deposited in the ArrayExpress database (http://www.ebi.ac.uk/arrayexpress [65]) under the accession number E-MTAB-8678.

Pre-processing of RNA-Seq data

Adapter dimer reads were removed using DimerRemover (https://sourceforge.net/projects/dimerremover/) and quality of each sample was assessed with FastQC v0.11.2 [66]. Transcript abundance was quantified using salmon v0.13.1 [67] with the flags --seqBias --gcBias --fldMean 50 --fldSD 1 --validateMappings -rangeFactorizationBins 4. Quantification results were concatenated to a count matrix using tximport R package [68]. After quality verification, one sample was identified as outlier and discarded from subsequent analysis.

Genome functional annotation

Gene models of the V3 annotation [69] were searched against the Araport11 (release 06.17.16, https://www.araport.org/) protein database with the blastx function of the DIAMOND version 0.9.19 software set to default parameters and reporting alignments in the 1% range of the top alignment score [70]. For each V3 gene model, the best blast hit was kept (1-to-1) and for multiple hits with the same score (1-to-many), the first hit was kept as representative keeping the other hits accessible. The corresponding gene annotations were obtained from the Araport11 gff file (release 06.22.16) and correspondences for V1 IDs were downloaded (https://urgi.versailles.inra.fr/Species/Vitis/Annotations).

Weighted gene co-expression network analysis (WGCNA)

A co-expression gene network was constructed using the WGCNA software package (v1.63) in R [71, 72]. After having filtered the low-expressed genes, 24,190 genes were used. The power β was set at 11. The module eigengenes were used to evaluate the association between the 26 identified gene modules and the experimental variables. Then, a kME value (module eigengene-based connectivity) was calculated for each gene with every module. For each module eigengene, highly correlated genes were filtered with a correlation coefficient > 0.80 and a p-value < 0.01 and used to perform the enrichment analysis described below.

Differential expression analysis and MapMan BIN enrichment.

The R package DESeq2 [73] was used to detect differentially expressed (DE) genes using the following thresholds: False Discovery Rate (FDR) adjusted p-value < 0.01 and |log2 fold change (LFC)| > 1. MapMan BIN codes [74] were attributed to each gene using Mercator 4 (http://plabipd.de/portal/mercator-ii-alpha-version-). Enrichment analysis was performed using Fisher’s test, results with a Bonferroni’s adjusted p-value <= 0.05 and an enrichment > 1 were selected.

Availability of data and materials

The datasets generated during the current study are available in the ArrayExpress database (http://www.ebi.ac.uk/arrayexpress [65]) under the accession number E-MTAB-8678.

Abbreviations

1103P:

Vitis berlandieri x V. rupestris cultivar 1103 Paulsen

4CL:

4-COUMAROYL COA-LIGASE;

ABA:

Abscisic acid

bHLH:

basic helix-loop-helix

C4H:

CINNAMATE-4-HYDROXYLASE

DE:

Differentially expressed

ERF:

ETHYLENE RESPONSE FACTOR

HP:

High phosphorus

IPS:

INDUCED BY PHOSPHATE STARVATION1

LP:

Low phosphorus

NLR:

NUCLEOTIDE-BINDING DOMAIN AND LEUCINE-RICH REPEAT-CONTAINING

P:

Phosphorus

PAL:

PHENYLALANINE AMMONIA-LYASE

PHO2:

PHOSPHATE2

PHR1:

PHOSPHATE STARVATION RESPONSE1

PHT1:

PHOSPHATE TRANSPORTER1

Pi:

Orthophosphate

PN:

V. vinifera cultivar Pinot Noir

PSI:

Phosphate-starvation-inducible

RNA-Seq:

RNA-sequencing

RT-qPCR:

quantitative reverse transcription polymerase chain reaction

SPX:

SYG/PHO81/XPR1

STS:

STILBENE SYNTHASE

TF:

Transcription factor

TKL:

TYROSINE KINASE-LIKE

VIT:

VACUOLAR IRON TRANSPORTER

WGCNA:

weighted gene co-expression network analysis.

References

  1. Chien PS, Chiang CP, Leong SJ, Chiou TJ. Sensing and signaling of phosphate starvation: from local to long distance. Plant Cell Physiol. 2018;59(9):1714–22.

    Article  CAS  PubMed  Google Scholar 

  2. Hinsinger P. Bioavailability of soil inorganic P in the rhizosphere as affected by root-induced chemical changes: a review. Plant Soil. 2001;237(2):173–95.

    Article  CAS  Google Scholar 

  3. Crombez H, Motte H, Beeckman T. Tackling plant phosphate starvation by the roots. Dev Cell. 2019;48(5):599–615.

    Article  CAS  PubMed  Google Scholar 

  4. Pant BD, Pant P, Erban A, Huhman D, Kopka J, Scheible WR. Identification of primary and secondary metabolites with phosphorus status-dependent abundance in Arabidopsis, and of the transcription factor PHR1 as a major regulator of metabolic changes during phosphorus limitation. Plant Cell Environ. 2015;38(1):172–87.

    Article  CAS  PubMed  Google Scholar 

  5. Mo XH, Zhang MK, Liang CY, Cai LY, Tian J. Integration of metabolome and transcriptome analyses highlights soybean roots responding to phosphorus deficiency by modulating phosphorylated metabolite processes. Plant Physiol Biochem. 2019;139:697–706.

    Article  CAS  PubMed  Google Scholar 

  6. Morcuende R, Bari R, Gibon Y, Zheng WM, Pant BD, Blasing O, Usadel B, Czechowski T, Udvardi MK, Stitt M, et al. Genome-wide reprogramming of metabolism and regulatory networks of Arabidopsis in response to phosphorus. Plant Cell Environ. 2007;30(1):85–112.

    Article  CAS  PubMed  Google Scholar 

  7. Nakamura Y. Phosphate starvation and membrane lipid remodeling in seed plants. Prog Lipid Res. 2013;52(1):43–50.

    Article  CAS  PubMed  Google Scholar 

  8. Lan P, Li WF, Schmidt W. 'Omics' approaches towards understanding plant phosphorus acquisition and use. In: Plaxton WC, Lambers H, editors. Phosphorus Metabolism in Plants, vol. 48; 2015. p. 65–97.

    Google Scholar 

  9. Zhang ZL, Liao H, Lucas WJ. Molecular mechanisms underlying phosphate sensing, signaling, and adaptation in plants. J Integr Plant Biol. 2014;56(3):192–220.

    Article  CAS  PubMed  Google Scholar 

  10. Bai FX, Chen CL, An JY, Xiao SY, Deng XX, Pan ZY. Transcriptome responses to phosphate deficiency in Poncirus trifoliata (L.) Raf. Acta Physiol Plant. 2014;36(12):3207–15.

    Article  CAS  Google Scholar 

  11. Fan FH, Cui BW, Zhang T, Qiao G, Ding GJ, Wen XP. The temporal Transcriptomic response of Pinus massoniana seedlings to phosphorus deficiency. PLoS One. 2014;9(8):11.

    Google Scholar 

  12. Kavka M, Polle A. Dissecting nutrient-related co-expression networks in phosphate starved poplars. PLoS One. 2017;12(2):e0171958.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  13. Tyagi W, Rai M. Root transcriptomes of two acidic soil adapted Indica rice genotypes suggest diverse and complex mechanism of low phosphorus tolerance. Protoplasma. 2017;254(2):725–36.

    Article  CAS  PubMed  Google Scholar 

  14. Zhang D, Zhang HY, Chu SS, Li HY, Chi YJ, Triebwasser-Freese D, Lv HY, Yu DY. Integrating QTL mapping and transcriptomics identifies candidate genes underlying QTLs associated with soybean tolerance to low-phosphorus stress. Plant Mol Biol. 2017;93(1–2):137–50.

    Article  CAS  PubMed  Google Scholar 

  15. Du QG, Wang K, Xu C, Zou C, Xie CX, Xu YB, Li WX. Strand-specific RNA-Seq transcriptome analysis of genotypes with and without low-phosphorus tolerance provides novel insights into phosphorus-use efficiency in maize. BMC Plant Biol. 2016;16:222.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  16. Grant RS, Matthews MA. The influence of phosphorus availability, scion, and rootstock on grapevine shoot growth, leaf area, and petiole phosphorus concentration. Am J Enol Vitic. 1996;47(2):217–24.

    Google Scholar 

  17. Grant RS, Matthews MA. The influence of phosphorus availability and rootstock on root system characteristics, phosphorus uptake, phosphorus partitioning, and growth efficiency. Am J Enol Vitic. 1996;47(4):403–9.

    Google Scholar 

  18. Gautier A, Cookson SJ, Hevin C, Vivin P, Lauvergeat V, Mollier A. Phosphorus acquisition efficiency and phosphorus remobilization mediate genotype-specific differences in shoot phosphorus content in grapevine. Tree Physiol. 2018;38(11):1742–51.

    CAS  PubMed  Google Scholar 

  19. Schreiner RP, Scagel CF. Nutrient uptake and distribution in a mature 'Pinot Noir' vineyard. HortScience. 2013;41(2):336–45.

    Article  Google Scholar 

  20. Lin W-Y, Huang T-K, Leong SJ, Chiou T-J. Long-distance call from phosphate: systemic regulation of phosphate starvation responses. J Exp Bot. 2013;65(7):1817–27.

    Article  PubMed  CAS  Google Scholar 

  21. Lin SI, Chiang SF, Lin WY, Chen JW, Tseng CY, Wu PC, Chiou TJ. Regulatory network of microRNA399 and PHO2 by systemic signaling. Plant Physiol. 2008;147(2):732–46.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Pant BD, Buhtz A, Kehr J, Scheible WR. MicroRNA399 is a long-distance signal for the regulation of plant phosphate homeostasis. Plant J. 2008;53(5):731–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Liu TY, Huang TK, Tseng CY, Lai YS, Lin SI, Lin WY, Chen JW, Chiou TJ. PHO2-dependent degradation of PHO1 modulates phosphate homeostasis in Arabidopsis. Plant Cell. 2012;24(5):2168–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Huang TK, Han CL, Lin SI, Chen YJ, Tsai YC, Chen YR, Chen JW, Lin WY, Chen PM, Liu TY, et al. Identification of downstream components of ubiquitin-conjugating enzyme PHOSPHATE2 by quantitative membrane proteomics in Arabidopsis roots. Plant Cell. 2013;25(10):4044–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Liu TY, Aung K, Tseng CY, Chang TY, Chen YS, Chiou TJ. Vacuolar Ca2+/H+ transport activity is required for systemic phosphate homeostasis involving shoot-to-root signaling in Arabidopsis. Plant Physiol. 2011;156(3):1176–89.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Cerri MR, Wang Q, Stolz P, Folgmann J, Frances L, Katzer K, Li X, Heckmann AB, Wang TL, Downie JA, et al. The ERN1 transcription factor gene is a target of the CCaMK/CYCLOPS complex and controls rhizobial infection in Lotus japonicus. New Phytol. 2017;215(1):323–37.

    Article  CAS  PubMed  Google Scholar 

  27. Matus JT, Aquea F, Espinoza C, Vega A, Cavallini E, Santo SD, Cañón P, de la Guardia AR-H, Serrano J, Tornielli GB, et al. Inspection of the grapevine BURP superfamily highlights an expansion of RD22 genes with distinctive expression features in berry development and ABA-mediated stress responses. PLoS One. 2014;9(10):e110372.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  28. Sengupta S, Mukherjee S, Basak P, Majumder AL. Significance of galactinol and raffinose family oligosaccharide synthesis in plants. Front Plant Sci. 2015;6:656.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Marshall E, Costa LM, Gutierrez-Marcos J. Cysteine-rich peptides (CRPs) mediate diverse aspects of cell–cell communication in plant reproduction and development. J Exp Bot. 2011;62(5):1677–86.

    Article  CAS  PubMed  Google Scholar 

  30. Omidbakhshfard Mohammad A, Proost S, Fujikura U, Mueller-Roeber B. Growth-regulating factors (GRFs): a Small transcription factor family with important functions in plant biology. Mol Plant. 2015;8(7):998–1010.

    Article  PubMed  CAS  Google Scholar 

  31. Vondras AM, Minio A, Blanco-Ulate B, Figueroa-Balderas R, Penn MA, Zhou Y, Seymour D, Ye Z, Liang D, Espinoza LK, et al. The genomic diversification of grapevine clones. BMC Genomics. 2019;20(1):972.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Song J-M, Guan Z, Hu J, Guo C, Yang Z, Wang S, Liu D, Wang B, Lu S, Zhou R, et al. Eight high-quality genomes reveal pan-genome architecture and ecotype differentiation of Brassica napus. Nature Plants. 2020;6(1):34–45.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Gao L, Gonda I, Sun H, Ma Q, Bao K, Tieman DM, Burzynski-Chang EA, Fish TL, Stromberg KA, Sacks GL, et al. The tomato pan-genome uncovers new genes and a rare allele regulating fruit flavor. Nat Genet. 2019;51(6):1044–51.

    Article  CAS  PubMed  Google Scholar 

  34. Lin WD, Liao YY, Yang TJW, Pan CY, Buckhout TJ, Schmidt W. Coexpression-based clustering of Arabidopsis root genes predicts functional modules in early phosphate deficiency signaling. Plant Physiol. 2011;155(3):1383–402.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Devaiah BN, Karthikeyan AS, Raghothama KG. WRKY75 Transcription Factor Is a Modulator of Phosphate Acquisition and Root Development in Arabidopsis. Plant Physiol. 2007;143(4):1789–801.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Chen Y-F, Li L-Q, Xu Q, Kong Y-H, Wang H, Wu W-H. The WRKY6 Transcription Factor Modulates PHOSPHATE1 Expression in Response to Low Pi Stress in Arabidopsis. Plant Cell. 2009;21(11):3554–66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Gutierrez RA, Ewing RM, Cherry JM, Green PJ. Identification of unstable transcripts in Arabidopsis by cDNA microarray analysis: rapid decay is associated with a group of touch- and specific clock-controlled genes. Proc Natl Acad Sci U S A. 2002;99(17):11513–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Narsai R, Howell KA, Millar AH, O'Toole N, Small I, Whelan J. Genome-wide analysis of mRNA decay rates and their determinants in Arabidopsis thaliana. Plant Cell. 2007;19(11):3418–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Cookson SJ, Yadav UP, Klie S, Morcuende R, Usadel B, Lunn JE, Stitt M. Temporal kinetics of the transcriptional response to carbon depletion and sucrose readdition in Arabidopsis seedlings. Plant Cell Environ. 2016;39(4):768–86.

    Article  CAS  PubMed  Google Scholar 

  40. Matsui A, Nakaminami K, Seki M. Biological function of changes in RNA metabolism in plant adaptation to abiotic stress. Plant Cell Physiol. 2019;60(9):1897–905.

    Article  CAS  PubMed  Google Scholar 

  41. Marschner H. Mineral nutrition of higher plants. 2nd ed. London: Acedemic Press; 1995.

  42. Chong JL, Poutaraud A, Hugueney P. Metabolism and roles of stilbenes in plants. Plant Sci. 2009;177(3):143–55.

    Article  CAS  Google Scholar 

  43. Vannozzi A, Wong DCJ, Höll J, Hmmam I, Matus JT, Bogs J, Ziegler T, Dry I, Barcaccia G, Lucchin M. Combinatorial regulation of Stilbene synthase genes by WRKY and MYB transcription factors in grapevine (Vitis vinifera L.). Plant Cell Physiol. 2018;59(5):1043–59.

    Article  CAS  PubMed  Google Scholar 

  44. Höll J, Vannozzi A, Czemmel S, D'Onofrio C, Walker AR, Rausch T, Lucchin M, Boss PK, Dry IB, Bogs J. The R2R3-MYB Transcription Factors MYB14 and MYB15 Regulate Stilbene Biosynthesis in Vitis vinifera. Plant Cell. 2013;25(10):4135–49.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  45. Vannozzi A, Donnini S, Vigani G, Corso M, Valle G, Vitulo N, Bonghi C, Zocchi G, Lucchin M. Transcriptional Characterization of a Widely-Used Grapevine Rootstock Genotype under Different Iron-Limited Conditions. Front Plant Sci. 2017;7(1994):1994.

    PubMed  PubMed Central  Google Scholar 

  46. Liu Y, Xie Y, Wang H, Ma X, Yao W, Wang H. Light and Ethylene Coordinately Regulate the Phosphate Starvation Response through Transcriptional Regulation of PHOSPHATE STARVATION RESPONSE1. Plant Cell. 2017;29(9):2269–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Seta A, Tabara M, Nishibori Y, Hiraguri A, Ohkama-Ohtsu N, Yokoyama T, Hara S, Yoshida K, Hisabori T, Fukudome A, et al. Post-translational regulation of the dicing activities of Arabidopsis DICER-LIKE 3 and 4 by inorganic phosphate and the redox state. Plant Cell Physiol. 2017;58(3):485–95.

    CAS  PubMed  Google Scholar 

  48. Yong-Villalobos L, Gonzalez-Morales SI, Wrobel K, Gutierrez-Alanis D, Cervantes-Perez SA, Hayano-Kanashiro C, Oropeza-Aburto A, Cruz-Ramirez A, Martinez O, Herrera-Estrella L. Methylome analysis reveals an important role for epigenetic changes in the regulation of the Arabidopsis response to phosphate starvation. Proc Natl Acad Sci U S A. 2015;112(52):E7293–302.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Schönberger B, Chen X, Mager S, Ludewig U. Site-dependent differences in DNA methylation and their impact on plant establishment and phosphorus nutrition in Populus trichocarpa. PLoS One. 2016;11(12):e0168623.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  50. Secco D, Whelan J, Rouached H, Lister R. Nutrient stress-induced chromatin changes in plants. Curr Opin Plant Biol. 2017;39:1–7.

    Article  CAS  PubMed  Google Scholar 

  51. Cochetel N, Escudie F, Cookson SJ, Dai ZW, Vivin P, Bert PF, Munoz MS, Delrot S, Klopp C, Ollat N, et al. Root transcriptomic responses of grafted grapevines to heterogeneous nitrogen availability depend on rootstock genotype. J Exp Bot. 2017;68(15):4339–55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Livigni S, Lucini L, Sega D, Navacchi O, Pandolfini T, Zamboni A, Varanini Z. The different tolerance to magnesium deficiency of two grapevine rootstocks relies on the ability to cope with oxidative stress. BMC Plant Biol. 2019;19:17.

    Article  Google Scholar 

  53. Chitarra W, Perrone I, Avanzato CG, Minio A, Boccacci P, Santini D, Gilardi G, Siciliano I, Gullino ML, Delledonne M, et al. Grapevine grafting: Scion transcript profiling and defense-related metabolites induced by rootstocks. Front Plant Sci. 2017;8:654.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Cookson SJ, Ollat N. Grafting with rootstocks induces extensive transcriptional re-programming in the shoot apical meristem of grapevine. BMC Plant Biol. 2013;13:147.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  55. Jensen PJ, Halbrendt N, Fazio G, Makalowska I, Altman N, Praul C, Maximova SN, Ngugi HK, Crassweller RM, Travis JW, et al. Rootstock-regulated gene expression patterns associated with fire blight resistance in apple. BMC Genomics. 2012;13:9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Prassinos C, Ko JH, Lang G, Iezzoni AF, Han KH. Rootstock-induced dwarfing in cherries is caused by differential cessation of terminal meristem growth and is triggered by rootstock-specific gene regulation. Tree Physiol. 2009;29(7):927–36.

    Article  CAS  PubMed  Google Scholar 

  57. Liu XY, Li J, Liu MM, Yao Q, Chen JZ. Transcriptome profiling to understand the effect of Citrus rootstocks on the growth of 'Shatangju' mandarin. PLoS One. 2017;12(1):22.

    Google Scholar 

  58. Liu N, Yang JH, Fu XX, Zhang L, Tang K, Guy KM, Hu ZY, Guo SG, Xu Y, Zhang MF. Genome-wide identification and comparative analysis of grafting-responsive mRNA in watermelon grafted onto bottle gourd and squash rootstocks by high-throughput sequencing. Mol Gen Genomics. 2016;291(2):621–33.

    Article  CAS  Google Scholar 

  59. Berdeja M, Nicolas P, Kappel C, Dai ZW, Hilbert G, Peccoux A, Lafontaine M, Ollat N, Gomes E, Delrot S. Water limitation and rootstock genotype interact to alter grape berry metabolism through transcriptome reprogramming. Hortic Res-England. 2015;2:13.

    Google Scholar 

  60. Yildirim K, Yagci A, Sucu S, Tunc S. Responses of grapevine rootstocks to drought through altered root system architecture and root transcriptomic regulations. Plant Physiol Biochem. 2018;127:256–68.

    Article  CAS  PubMed  Google Scholar 

  61. Xia C, Zheng Y, Huang J, Zhou X, Li R, Zha M, Wang S, Huang Z, Lan H, Turgeon R, et al. Elucidation of the mechanisms of long-distance mRNA movement in a Nicotiana benthamiana/tomato heterograft system. Plant Physiol. 2018;177(2):745–58.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Jaillon O, Aury J-M, Noel B, Policriti A, Clepet C, Casagrande A, Choisne N, Aubourg S, Vitulo N, Jubin C, et al. The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007;449(7161):463–7.

    Article  CAS  PubMed  Google Scholar 

  63. Hellemans J, Mortier G, De Paepe A, Speleman F, Vandesompele J. qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biol. 2007;8(2):R19.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  64. R: A language and environment for statistical computing. http://www.R-project.org. Accessed Aug 2018.

  65. Athar A, Fullgrabe A, George N, Iqbal H, Huerta L, Ali A, Snow C, Fonseca NA, Petryszak R, Papatheodorou I, et al. ArrayExpress update - from bulk to single-cell expression data. Nucleic Acids Res. 2019;47(D1):D711–5.

    Article  CAS  PubMed  Google Scholar 

  66. Andrews S: FastQC: a quality control tool for high throughput sequence data. Available online at: http://www.bioinformatics.babraham.ac.uk/projects/fastqc. 2010.

    Google Scholar 

  67. Patro R, Duggal G, Love MI, Irizarry RA, Kingsford C. Salmon provides fast and bias-aware quantification of transcript expression. Nat Methods. 2017;14(4):417–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Soneson C, Love MI, Robinson MD. Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences. F1000Research. 2015;4:1521.

    Article  PubMed  CAS  Google Scholar 

  69. Canaguier A, Grimplet J, Di Gaspero G, Scalabrin S, Duchene E, Choisne N, Mohellibi N, Guichard C, Rombauts S, Le Clainche I, et al. A new version of the grapevine reference genome assembly (12X.v2) and of its annotation (VCost.v3). Genomics Data. 2017;14:56–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015;12(1):59–60.

    Article  CAS  PubMed  Google Scholar 

  71. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  72. Langfelder P, Horvath S. Fast R functions for robust correlations and hierarchical clustering. J Stat Softw. 2012;46(11):1–17.

    Article  Google Scholar 

  73. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  74. Thimm O, Blasing O, Gibon Y, Nagel A, Meyer S, Kruger P, Selbig J, Muller LA, Rhee SY, Stitt M. MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 2004;37(6):914–39.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

RNA-sequencing was performed by the GenomEast platform, a member of the France Génomique consortium (ANR-10-INBS-0009).

Funding

This study has been carried out with financial support from the French National Research Agency (ANR) in the frame of the Investments for the future Program, within the Cluster of Excellence COTE (ANR-10-LABX-45); the Ph.D. studentship of AG and consumable costs were paid by this grant.

Author information

Authors and Affiliations

Authors

Contributions

AG, CH and IM executed the experiments and collected data, NC analyzed the data, AG and SC prepared figures, SC drafted the manuscript, and all authors revised, read and approved the manuscript.

Corresponding author

Correspondence to Sarah J. Cookson.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests except Sarah Jane Cookson who is an Associate Editor of BMC Plant Biology.

Additional information

Publisher’s Note

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

Supplementary information

Additional file 1.

Genes differentially expressed between high or low phosphate treatments (HP and LP respectively) in four scion/rootstock combinations of grapevine.

Additional file 2.

Genes differentially expressed between the roots of two grapevine genotypes grown grafted under high or low phosphate supply (HP and LP respectively).

Additional file 3.

MapMan BINs enriched in the genes differentially expressed in the roots between the two grapevine genotypes grown grafted under high phosphate conditions.

Additional file 4.

Genes with a significant interaction in their gene expression response to low versus high phosphate supply between two homo-grafts of grapevine

Additional file 5.

MapMan BINs enriched in the genes differentially expressed in the roots in response to low phosphate between two homo-grafts of grapevine: Vitis vinifera cv. Pinot noir and V. rupestris x V. berlandieri cv. 1103 Paulsen.

Additional file 6.

Genes differentially expressed between hetero- and homo-grafts of grapevine grown under either high or low phosphate supply.

Additional file 7.

MapMan BINs enriched in the 24 genes down-regulated in the roots between the scion/rootstock combinations 1103P/PN versus PN/PN grown under both high and low phosphate supply. PN: Vitis vinifera cv. Pinot noir; 1103P: V. rupestris x V. berlandieri cv. 1103 Paulsen.

Additional file 8.

MapMan BINs enriched in the 661 and 505 genes up- and down-regulated respectively in the roots between the scion/rootstock combinations PN/1103P versus 1103P/1103P grown under both high and low phosphate supply. PN: Vitis vinifera cv. Pinot noir; 1103P: V. rupestris x V. berlandieri cv. 1103 Paulsen.

Additional file 9.

Genes with a significant interaction in their gene expression response to low versus high phosphate supply between the scion/rootstock combinations 1103P/PN and PN/PN. PN: Vitis vinifera cv. Pinot noir; 1103P: V. rupestris x V. berlandieri cv. 1103 Paulsen.

Additional file 10.

Genes with a significant interaction in their gene expression response to low versus high phosphate supply between the scion/rootstock combinations PN/1103P and 1103P/1103P. PN: Vitis vinifera cv. Pinot noir; 1103P: V. rupestris x V. berlandieri cv. 1103 Paulsen.

Additional file 11.

Gene module membership. Module eigengene-based connectivity values for each gene in each module identified in a WGCNA analysis of genes expressed in the roots of grafted grapevine grown under two phosphate supplies.

Additional file 12.

Primers used for RT-qPCR experiments.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Gautier, A.T., Cochetel, N., Merlin, I. et al. Scion genotypes exert long distance control over rootstock transcriptome responses to low phosphate in grafted grapevine. BMC Plant Biol 20, 367 (2020). https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-020-02578-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-020-02578-y

Keywords