# Transcriptional Characterization of a Widely-Used Grapevine Rootstock Genotype under Different Iron-Limited Conditions

Alessandro Vannozzi, Silvia Donnini, Gianpiero Vigani, Massimiliano Corso, Giorgio Valle, Nicola Vitulo, Claudio Bonghi, Graziano Zocchi, Margherita Lucchin

<figure><img src="image_1.png" /></figure>

To cite this version:

Alessandro Vannozzi, Silvia Donnini, Gianpiero Vigani, Massimiliano Corso, Giorgio Valle, et al. Transcriptional Characterization of a Widely-Used Grapevine Rootstock Genotype under Different Iron-Limited Conditions. Frontiers in Plant Science, 2017, 7, 10.3389/fpls.2016.01994. [https://doi.org/10.3389/fpls.2016.01994](https://doi.org/10.3389/fpls.2016.01994)

HAL Id: hal-02861101

https://hal.inrae.fr/hal-02861101v1

Submitted on 23 Oct 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Transcriptional Characterization of a Widely-Used Grapevine Rootstock Genotype under Different Iron-Limited Conditions

Alessandro Vannozzi¹,²*, Silvia Donnini³, Gianpiero Vigani³, Massimiliano Corso¹,²†,
Giorgio Valle⁴, Nicola Vito¹⁰†, Claudio Bonghi¹,², Graziano Zocchi³ and
Margherita Lucchin¹,²

¹ Dipartimento di Agronomia Animali Alimenti Risorse Naturali e Ambiente, Università di Padova, Legnaro, Italy, ² Centro Interdipartimentale per la Ricerca in Viticoltura ed Enologia, Conegliano, Italy, ³ Dipartimento di Scienze Agrarie e Ambientali, Università di Milano, Milano, Italy, ⁴ Centro di Ricerca Interdipartimentale per le Biotecnologie Innovative, Padova, Italy

OPEN ACCESS

**Edited by:**
Patrick Brown,
University of California, Davis, USA

**Reviewed by:**
Ümit Banş Kutman,
Konya Food and Agriculture
University, Turkey
Stefania Astolfi,
Tuscia University, Italy

*Correspondence:
Vannozzi Alessandro
alessandro.vannozzi@unipd.it

†Present Address:
Massimiliano Corso,
Laboratory of Plant Physiology and
Molecular Genetics, Université Libre
de Bruxelles, Brussels, Belgium
Nicola Vito¹,
Department of Biotechnology,
University of Verona, Verona, Italy

**Specialty section:**
This article was submitted to
Plant Nutrition,
a section of the journal
Frontiers in Plant Science

**Received:** 03 October 2016
**Accepted:** 16 December 2016

**Published:** 05 January 2017

**Citation:**
Vannozzi A, Donnini S, Vigani G, Corso M, Valle G, Vito N, Bonghi C, Zocchi G and Lucchin M (2017) Transcriptional Characterization of a Widely-Used Grapevine Rootstock Genotype under Different Iron-Limited Conditions. Front. Plant Sci. 7:1994. doi: 10.3389/fpls.2016.01994

Iron chlorosis is a serious deficiency that affects orchards and vineyards reducing quality and yield production. Chlorotic plants show abnormal photosynthesis and yellowing shoots. In grapevine iron uptake and homeostasis are most likely controlled by a mechanism known as “Strategy I,” characteristic of non-graminaceous plants and based on a system of soil acidification, iron reduction and transporter-mediated uptake. Nowadays, grafting of varieties of economic interest on tolerant rootstocks is widely used practice against many biotic and abiotic stresses. Nevertheless, many interspecific rootstocks, and in particular those obtained by crossing exclusively non-vinifera genotypes, can show limited nutrient uptake and transport, in particular for what concerns iron. In the present study, 101.14, a commonly used rootstock characterized by susceptibility to iron chlorosis was subjected to both Fe-absence and Fe-limiting conditions. Grapevine plantlets were grown in control, Fe-deprived, and bicarbonate-supplemented hydroponic solutions. Whole transcriptome analyses, via mRNA-Seq, were performed on root apices of stressed and unstressed plants. Analysis of differentially expressed genes (DEGs) confirmed that Strategy I is the mechanism responsible for iron uptake in grapevine, since many orthologs genes to the *Arabidopsis* "ferrome" were differentially regulated in stressed plant. Molecular differences in the plant responses to Fe absence and presence of bicarbonate were also identified indicating the two treatments are able to induce response-mechanisms only partially overlapping. Finally, we measured the expression of a subset of genes differentially expressed in 101.14 (such as *IRT1*, *FERRITIN1*, *bHLH38/39*) or known to be fundamental in the "strategy I" mechanism (*AHA2* and *FRO2*) also in a tolerant rootstock (M1) finding important differences which could be responsible for the different degrees of tolerance observed.

**Keywords:** Vitis, ferrome, mRNA-Seq, strategy I, micronutrients
INTRODUCTION

Nowadays, much of the European viticulture is based primarily
on grafting, where the scion is represented by a *Vitis vinifera*
cultivar and the rootstock is either a North American *Vitis*
species or an interspecific *Vitis* hybrid carrying the resistance to
phyloxera (Weaver, 1976). Although, at first, the replacement of
the entire root system of “vinifera” varieties with “non-vinifera”
or “American” species was seen as a sort of contamination of the
purity and quality of wine, subsequently, the use of grapevine
rootstocks spread also over those countries where viticulture
constituted a young and new challenge and where phylloxera
did not represent a menace yet (Arrigo and Arnold, 2007;
Gregory et al., 2013). The explanation of this trend relies on
the fact that rootstocks were found not only to confer resistance
to diseases, but also to imply a larger range of advantages,
controlling numerous physiological processes at the level of scion
such as biomass accumulation (Gregory et al., 2013), fruit quality
(Walker et al., 2002, 2004), grape berry development and ripening
(Corso et al., 2016), and the ability to respond to many biotic
and abiotic stresses (Marguerit et al., 2012; Meggio et al., 2014)
including salinity (Fisarakis et al., 2001) and drought (Gambetta
et al., 2012; Marguerit et al., 2012; Corso et al., 2015). Despite
this, since their first introduction in the agronomical practice,
one of the problems winegrowers had to cope with was the
adaptation of these American species to European soils. In fact,
sometimes interspecific rootstocks, and in particular those ones
obtained by crossing exclusively *non-vinifera* species such as the
widely-used 101.14 genotype (*V. riparia* × *V. rupestris*), can
show several problems related to the nutrient uptake (Bavaresco
et al., 1993) in particular for what concerns iron (Fe). Although
this element is quite abundant in the earth's crust as well as in
cultivated soils, its acquisition by plants is often impaired since it
is insoluble in well-oxygenated and alkaline environments (such
as calcareous soils) being mainly accumulated as oxi-hydroxide
compounds (Tagliavini and Rombolà, 2001). Among the essential
micronutrients in plants, Fe is required in the greatest abundance.
Because its capability to change redox state and being the
key component of both hemoproteins and iron sulfur proteins
(including cytochromes, catalases, peroxidases, and ferredoxin) it
plays a pivotal role in many important cellular processes such as
photosynthesis, respiration, chlorophyll biosynthesis, nitrogen,
and sulfur assimilation, hormone biosynthesis, etc. (Hänsch
and Mendel, 2009; Kobayashi and Nishizawa, 2012; Balk and
Schaedler, 2014; Briat et al., 2015). To give an idea of the global
importance of Fe fine tune regulation in all these fundamental
biologic processes, on the one hand it was demonstrated to be
a limiting factor for biomass and seed yield in different crops
including tomato (*Solanum lycopersicum*; Jin et al., 2009), spinach
(*Spinacea oleracea*; Jin et al., 2013), and rice (*Oryza sativa*;
Takahashi et al., 2001). On the other hand, its excess is toxic to
plant cells due to Fe ability to change its oxidation status that
leads to ROS production (Thomine and Lanquar, 2011).

In order to counteract the low availability of Fe in the soil,
during their evolution plants have adopted two main uptake
mechanisms, proposed for the first time by Römheld and
Marschner (1986), and known as Strategy I in non-graminaceous
plants and Strategy II in graminaceous plants. Since the end

of the 1990s, many different genes involved in these strategies
have been identified and almost perfectly fit the model scheme
proposed by Römheld and Marschner. The Strategy I consists
of three main steps: (a) the reduction of Fe(III)-Fe(II) by a
ferric reductase oxidase enzyme (FRO2) located in the apoplast
(Robinson et al., 1999); (b) the import of the resulting Fe(II)
ions through the plasmalemma by the activity of a divalent iron-
regulated transporter (IRT1; Eide et al., 1996; Vert et al., 2002)
belonging to the ZIP multigenic family (Guerinot, 2000) and
(c) the acidification of local rhizosphere mediated by plasma
membrane H⁺-ATPase (AHA2/7) proton extrusion (Dell'Orto
et al., 2002; Santi and Schmidt, 2009). Orthologs of Arabidopsis
AHA2, FRO2, and IRT1 genes have been identified and studied
in various plant species including bean (*Pisum sativum*; Waters
et al., 2002), tomato (*Solanum lycopersicon*; Zamboni et al.,
2012), peanut (*Arachis hypogea*; Xiong et al., 2014), cucumber
(*Cucumis sativus*; Santi et al., 2005). Although the key steps
of Strategy I iron uptake have been exhaustively described,
many other processes seem to be involved in the mechanisms
of Fe uptake such as the production and subsequent excretion
of iron-chelators. The type of chelators released by the roots
appears to be species specific, for example flavins for tobacco
*(Nicotiana tabacum)* and *Medicago trucantula* (Vorwieger et al.,
2007; Rodríguez-Celma et al., 2013), and phenolic compounds
for Arabidopsis and red clover (*Trifolium pratense*; Jin et al., 2007;
Fourcroy et al., 2014; Schmidt et al., 2014). Plant responses to Fe
deficiency have been recently analyzed on the basis of large-scale
changes not only in transcriptome (Thimm et al., 2001; O'Rourke
et al., 2007a,b, 2009; Buckhout et al., 2009; Forner-Giner et al.,
2010; Yang et al., 2010), but also in proteome (Brumbarova et al.,
2008; Li et al., 2008; Donnini et al., 2010; Lan et al., 2011) and
metabolome (Rellán-Álvarez et al., 2009). Recently, a set of 92
transcripts that robustly reflect the transcriptional response of
*Arabidopsis* to Fe deficiency has been described as the “ferrome”
by Schmidt and Buckhout (2011) and consists of a list of genes
considered to be involved in the basic response to iron deficiency.

Focusing on grapevine, the only information available on the
plant response to Fe deficiency derives from several physiological
and biochemical studies (Brancadoro et al., 1995; Nikolic et al.,
2000; Bertamini et al., 2002; Piagnani et al., 2003; Bavaresco
et al., 2006; Russo et al., 2010), most of which performed on
*Vinifera* varieties, not considering that in the majority of cases,
roots, which are the first organ perceiving the stress and actively
responding to it, belong to American or interspecific hybrids.

In the present study we performed a physiological and
genome-wide transcriptional characterization of 101.14
(*V. riparia* × *V. rupestris*), a commonly used grapevine
rootstock susceptible to iron chlorosis, in response to Fe
deficiency. Not only we analyzed the transcriptome effect of
iron deprivation itself, but we also compared it to a condition
of low iron availability determined by high concentration of
bicarbonate, a phenomenon quite common in nature, which
better mimics what happens in the vineyard. Our approach led
to the identification of a set of genes differentially expressed
(DEG) in response to these nutritional stresses. Many DEGs
were orthologs of the so-called *Arabidopsis* ferrome, confirming
the grapevine roots behave similarly to *Arabidopsis* under
iron shortage, undertaking a Strategy I mechanism of iron
homeostasis. Finally, we compared the expression of a subset of genes found to be differentially regulated in 101.14 (such as *IRT1*, *FERRITIN1*, *bHLH38/39*) or known to be fundamental in the “strategy I” mechanism (*AHA2* and *FRO2*) both in 101.14 and in a tolerant experimental rootstock (M1), finding important differences which could be responsible for the different degrees of resistance observed.

# MATERIALS AND METHODS

## Plant Material and Stress Treatments

Plants of 101.14 (*Vitis riparia* × *Vitis rupestris*) and M1 (106/8 [*V. riparia* × (*V. cordifolia* × *V. rupestris*)] × *V. berlandieri* cv. Ressequier n.1) rootstock genotypes were obtained from rooted canes and grown in hydroponic culture in controlled conditions, with a relative humidity of ~65–75%, photoperiod of 16 h light/8 h dark, PAR of 200 mmol m⁻¹s⁻¹ and temperature ranging from 25°C day/21°C night. Hydroponic cultures were set up into plastic light-reflecting containers, to avoid formation of algae within the nutritive solution. Standard solution, obtained from a modification of the Hoagland solution, was composed of 2 mM Ca(NO₃)₂, 0.75 mM K₂SO₄, 0.5 mM KH₂PO₄, 0.65 mM MgSO₄, 0.1 mM Fe-ethylenediaminetetraacetic acid (EDTA), 10 µM H₃BO₃, 0.5 µM CuSO₄, 0.5 µM ZnSO₄, 1 µM MnSO₄, and 0.05 µM (NH₄)₆Mo₇O₂₄. The experiments were carried out using six rooted rootstocks per tank (20 L). The solution was aerated bubbling filtered air and was changed every week. Plants were left into a half-nutritive solution for 10 days, in order to allow them to adapt to the growth conditions. Subsequently, stress treatments were imposed according to Donnini et al. (2012). Acclimated plants were transferred into containers containing a modified nutritive standard solution according to the different conditions applied: (i) control (+Fe): 100 µM Fe-EDTA (pH 6 ~ 6.2); (ii) absence of iron (−Fe): 0 µM Fe-EDTA (pH 6 ~ 6.2); (iii) presence of Bicarbonate (FeBic): 100 µM Fe-EDTA, 10 mM NaHCO₃, and 0.5 g/L CaCO₃ (pH 8.2 ~ 8.3). The choice to culture plants within a liquid instead of solid-agar culture was taken in order to avoid perturbations during the experimental procedure, such as the physical injury of plants during transferring from normal to stress conditions and the roots carry-over of residual amounts of Fe-containing agarose. The system allowed the bulk transfer of plants with intact roots into Fe-sufficient (+Fe), Fe-deficient (−Fe), and Fe-limiting (FeBic) media within seconds and with minimal mechanical damage. After 10 days of hydroponic culture in the modified solutions, the root apices of FIVE plants per condition (+Fe, −Fe, FeBic) were collected and pooled together for biochemical and molecular analyses. Only actively growing root apices were sampled, trying to collect only those ones that were exposed to the entire kinetic of stress and avoiding the new ones generated during the experimental trial. Sampling was performed in three biological replicates consisting of three independent experiments.

## Physiological and Biochemical Analyses of 101.14 under Different Growth Conditions

Leaves and roots of 101.14 plants grown in the above mentioned conditions were analyzed by mean of Inductively Coupled Plasma-Mass Spectrometry (ICP-MS) to determine their content

in macro- and microelements. Sampled tissues were dried and then mineralized in HNO₃ by using a Microwave Digestion System (Multiwave ECO). In order to reduce Fe apoplast contamination, roots were collected and rinsed with a solution containing 0,5 mM CaSO₄, 10 mM of Na EDTA for 30 min and washed with distilled water several times. Ionome content was determined by ICP-MS (Aurora M90 BRUKER). At the shoot level, gas exchanges were measured with CIRAS portable instrument according to Rolli et al. (2015) and chlorophyll content were determined by using DUALEX system according to Cerovic et al. (2012) on three expanded leaves per plants from three different plants for each treatment. Data collected for both ionome profile as well as *in vivo* leaf physiological analyses are from three independent experiments.

## RNA Extraction and Libraries Construction

Total RNA for both mRNA-Seq and RT-PCR was extracted from ~80 mg of frozen root apices sampled 10 days after imposition of the stress using the “SpectrumTM Plant total RNA Kit” according to manufacturer’s instructions (Sigma®). RNA quality and concentration were determined by both UV spectrophotometry (NanodropTM, Thermo Scientific™) and gel electrophoresis analyses. NGS analysis was performed in 101.14 plants on two biological replicates per each condition. mRNA was purified from the total RNA using the Dynabeads mRNA Direct kit (Invitrogen pn 610.12). A variable quantity of mRNA ranging from 0.4 to 1.6% with respect to the amount total of RNA was obtained. Samples for ligation sequencing were prepared according to the “SOLiD Whole transcriptome library preparation” protocol (pn 4452437 Rev. B). Before RNase III digestion samples were purified with Purelink RNA micro kit columns (Invitrogen™), digested from 3' to 10' depending on the starting amount of mRNA, retro-transcribed, size-selected using Agencourt AMPure XP beads (Beckman Coulter pn A63881) and barcoded during the final amplification. Obtained libraries (2.5 samples per lane) were sequenced using Applied Biosystems, SOLiD 5500 XL, producing paired-end reads of 75 and 35 nucleotides for the forward and reverse sequences, respectively. Reads were aligned to the v1 prediction of grapevine PN40024 reference genome (http://genomes.cribi.unipd.it/grape) using PASS aligner (Campagna et al., 2009). The percentage identity was set to 90% and one gap was allowed whereas the quality filtering parameters were set automatically by PASS. A minimum reads length cut-off of 50 and 30 nt was set for the forward sequences and reverse reads, respectively. The PASS-pair tool of the PASS package was used to perform the pairing between the forward and the reverse reads and to selected only those sequences that were uniquely aligned. Finally htseq-counts program (http://www.huber.embl.de/users/anders/HTSeq/doc/count.html) was used to quantify gene abundance.

## mRNA-Seq Analyses and Bioinformatics

Statistical analyses for discovering differentially expressed genes (DEGs) were performed with DEseq2 R package (http://www.r-project.org/; Love et al., 2014) according to an FDR-corrected *p* < 0.05. Functional and Gene Ontology associations were inferred from an improved version of the VitisNet Gene Annotation Reference file (Grimplet et al., 2012). Graphical
representation of genes expression, hierarchical, and k-Means (KPM) clustering were performed using the Multi Experiment Viewer software (MeV; Saeed et al., 2006) using the Pearson correlation coefficient. Venn diagrams and identification of common and specific DEGs were performed using Venn diagram plotter (https://omics.pnl.gov/software/venn-diagram-plotter).

## Real Time Reverse Transcriptase-PCR

For quantitative real-time PCR analyses (qPCRs) cDNA was synthesized starting from 1 µg of total RNA using the “Transcriptor First Strand cDNA Synthesis Kit” (Roche) with the oligo (dT)₂₀ primer according to manufacturer’s instructions. qPCRs were carried out in three technical replicates for each sample. Gene specific primers were designed using Primer3 Plus software (Untergasser et al., 2007) and are the following ones: *LDOX* (VIT_06s0004g00770), 5'- AAGCGA GGTCATCGGTCCTT-3' and 5'- TCTTCCCATCGTGTGC CTTC-3'; *FERRITIN1* (VIT_13s0067g01840), 5'- GCAGGT GGAAGCCATTAAAA-3' and 5'- GCAGCAACAACTCCT CCATT-3'; *NODULIN* (VIT_19s0093g00320), 5'- GGTAA GAGGTGCATGGAGGA-3' and 5'- TTCCTGTAATCCAC CACCAAG-3'; *AHA2* (VIT_11s0052g00620), 5'-GACTCTC CATGGCCTTCAAC-3' and 5'- ACCCTTCAGCTTGACC ACTG-3'; *IRT1* (VIT_10s0042g01100), 5'- CGGAAATCGA AGTTGCAGAT-3' and 5'- GCTGCCACAAGAGGCTTTA T-3'; *FRO2* (VIT_15s0046g01900), 5'- TCCGCTTTCGGGT GCCCCTT-3' and 5'- TGGCGCCAAGGCATAGCCAG-3'; *bHLH38/39* (VIT_08s0007g05300), 5'- TTCAGCAAACAG

GCTCAGTG-3' and 5'- TCACTCAACACCTCGCACTC-3'; *POPEYE* (VIT_17s0000g06000), 5'- TGCCTCAAAAAGGA GAATGC-3' and 5'- GGAGCAGAGTTGAGGTCAGG-3'; *BRUTUS* (VIT_17s0000g03440), 5'- TTGCCTGGGGATGAA GTTAG-3' and 5'- GACCAACAAAGCATCAAGCA-3'. Melt curve analysis, agarose gel electrophoresis, and DNA sequencing validated the absence of illegitimate cross-amplification of other genes. Real-time PCR were performed using a Sybr green method on a StepOne Plus Real-Time PCR System (Applied Biosystems™) thermal cycler. Each 10 µl PCR reaction contained 0.6 µl of each primer (10 µM), 1 µl of diluted cDNA, 2X Fast Start Sybr green (Applied Biosystems™) and sterile water. The thermal cycling conditions used were 94°C for 10 min followed by 40 cycles of: 95°C for 30 s, 60°C for 30 s, and 72°C for 30 s, followed by a melt cycle with 1°C increments from 55 to 96°C. Real time PCR data processing was performed using the Δ(ΔCT) method. After testing the suitability of ubiquitin, actin and elongation factor for use of reference genes, elongation factor (VIT_12s0035g01120) was selected for normalization of all samples analyzed.

## RESULTS AND DISCUSSION

### Characterization of 101.14 Morphological and Physiological Parameters under +Fe, -Fe, and FeBic Conditions

After 10 days of culture in both Fe impaired conditions (—Fe and FeBic), both leaves and roots of 101.14 rootstocks showed the

<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Shoot and root apparatus of 101.14 plants grown under different Fe-supply conditions. Leaf detail of 101.14 plants grown upon —Fe (A), +Fe (D), and FeBic conditions (G). Root apparatus and relative detail of —Fe (B,C), +Fe (E,F), and FeBic (H,I).</figcaption></figure>
typical symptoms of chlorosis compared to control (Figure 1). At the shoot level, the widespread yellowing of leaf internerval tissue was coupled to the decrease of net photosynthesis (Pn), stomatal conductance (Gs), and transpiration rate (E) (Figure 2A), together with a reduction in chlorophyll content determined as chlorophyll meter readings (SPAD index; Figure 2B). Intracellular carbon dioxide concentration (Ci) was significantly increased only in plants grown in -Fe conditions whereas it maintained levels similar to control (+Fe) in plants grown in presence of bicarbonate (Figure 2A). According to previous observations on other plant species, such as tomato (Zamboni et al., 2012), the root apparatus of Fe-deprived 101.14 plants developed more lateral roots (Figure 1).

## Uptake and Translocation of Iron and Major Macro- and Micronutrients in 101.14

The scarcity or excess of a nutrient such as Fe in the substrate not only determines an alteration of the concentration of that particular element both in root and leaves tissues, but can significantly alter the profile of many other micro- and macronutrients within the plant resulting in important changes in primary and secondary metabolism. In fact, the minerals or ions taken up by plants can be used by the cells as building blocks for the synthesis of organic molecules, or as catalyzers in enzymatic reactions. Moreover, they play important roles as

counter-ion in the transport of ionized molecules by the plant, in the osmotic regulation, and in maintaining the electrochemical potential gradient across the cell. Based on this awareness, leaves and roots of 101.14 plants grown in the above mentioned conditions were analyzed by mean of ICP-MS to determine their content in macro- and microelements (Table 1). A first obvious observation, which leapt out analyzing the root micronutrient profile, was the drop of Fe concentration in both -Fe (848.5 µM g⁻¹ DW) and FeBic (989.1 µg g⁻¹ DW) conditions with respect to control (+Fe, 2104.2 µg g⁻¹ DW). In FeBic treatment, roots still maintained a basal capacity to take up Fe²⁺ with respect to the condition of total Fe-absence showing values slightly higher than in -Fe, but this was not reflected in leaf tissues, which in FeBic conditions showed values even lower than in -Fe (71.2 µg g⁻¹ DW; Table 1), suggesting that bicarbonate could impair not only those processes involved in Fe-uptake at the root level but also in its translocation to the shoot. In most of cases all other micro-nutrients, including manganese (Mn), copper (Cu), zinc (Zn), and molybdenum (Mo) showed an opposite behavior with respect to Fe, being accumulated at higher concentration in roots of plants grown in -Fe conditions. Generally, this trend was not observed in roots of those plants cultured in presence of bicarbonate (FeBic), apart from Cu that was the only microelement showing an increased concentration with respect to control. In all other cases the ionome homeostasis was not significantly affected by bicarbonate compared to control, except for Mo, whose concentration was markedly diminished with respect to control (Table 1). We recently demonstrated that Fe and Mo mutually impact their homeostasis in cucumber plant (Vigani et al., 2016). We indeed observed a strong accumulation of Mo in root of Fe-deficient plants. Particularly, we observed Fe accumulation in Mo-deficient mitochondria as well as Mo accumulation in Fe-deficient mitochondria, suggesting that such

A

<table><thead><tr><th rowspan="2">Leaf gas exchanges</th><th colspan="3">101.14 genotype</th></tr><tr><th>-Fe</th><th>+Fe</th><th>FeBic</th></tr></thead><tbody><tr><td>Pn (µmol CO<sub>2</sub> s<sup>-1</sup> m<sup>-2</sup>)</td><td>0.63<sup>c</sup></td><td>2.53<sup>a</sup></td><td>1.31<sup>b</sup></td></tr><tr><td>E (mmol H<sub>2</sub>O s<sup>-1</sup> m<sup>-2</sup>)</td><td>0.83<sup>b</sup></td><td>1.42<sup>a</sup></td><td>1.01<sup>b</sup></td></tr><tr><td>G<sub>s</sub> (mmol CO<sub>2</sub> s<sup>-1</sup> m<sup>-2</sup>)</td><td>61.60<sup>b</sup></td><td>121.03<sup>a</sup></td><td>86.32<sup>b</sup></td></tr><tr><td>Ci (ppm)</td><td>301.40<sup>a</sup></td><td>275.33<sup>b</sup></td><td>276.02<sup>b</sup></td></tr><tr><td>VPD (mb)</td><td>13.06</td><td>12.46</td><td>14.73</td></tr></tbody></table>

B

Chlorophyll content

<figure><img src="image_4.png" /><figcaption>FIGURE 2 | 101.14 physiological parameters upon +Fe, -ig, and FeBic.</figcaption></figure>

(A) Leaf gas exchange parameters were measured by means of portable CIRAS instrument. Values followed by different letters (a, b, and c) are statistically different (test ANOVA, *p* < 0.05). Pn, net photosynthesis; E, evapotranspiration; Gs, Stomatal conductance; Ci, CO2 internal concentration; VPD, Vapor Pressure Deficit. (B) Chlorophyll content (SPAD index) upon -Fe, Fe, and FeBic conditions.

<table><caption>TABLE 1 | Macro- and micronutrient concentrations (µg/g DW) in leaves and roots of 101.14 plants upon -Fe, +Fe, and FeBic.</caption><thead><tr><th rowspan="2">Macro nutrients</th><th colspan="3">Leaves</th><th colspan="3">Roots</th></tr><tr><th>-Fe</th><th>+Fe</th><th>FeBic</th><th>-Fe</th><th>+Fe</th><th>FeBic</th></tr></thead><tbody><tr><td>P</td><td>62637b</td><td>73119a</td><td>22883c</td><td>6470b</td><td>10744a</td><td>5421c</td></tr><tr><td>Mg</td><td>4647b</td><td>4736c</td><td>3818a</td><td>6617b</td><td>9741a</td><td>6410b</td></tr><tr><td>K</td><td>34311b</td><td>51468c</td><td>51544a</td><td>15072c</td><td>39565a</td><td>34457b</td></tr><tr><td>Ca</td><td>12550a</td><td>10583c</td><td>11906b</td><td>15072c</td><td>16356b</td><td>33464a</td></tr></tbody></table>

<table><thead><tr><th rowspan="2">Micro nutrients</th><th colspan="3">Leaves</th><th colspan="3">Roots</th></tr><tr><th>-Fe</th><th>+Fe</th><th>FeBic</th><th>-Fe</th><th>+Fe</th><th>FeBic</th></tr></thead><tbody><tr><td>Mn</td><td>63a</td><td>76a</td><td>39b</td><td>71a</td><td>39b</td><td>34b</td></tr><tr><td>Fe</td><td>95b</td><td>145a</td><td>71c</td><td>849c</td><td>2104a</td><td>989b</td></tr><tr><td>Cu</td><td>13a</td><td>8b</td><td>14a</td><td>231a</td><td>39c</td><td>137b</td></tr><tr><td>Zn</td><td>30</td><td>25</td><td>31</td><td>111a</td><td>38b</td><td>40b</td></tr><tr><td>Mo</td><td>3c</td><td>4b</td><td>41a</td><td>60a</td><td>26b</td><td>18c</td></tr></tbody></table>

Content in macro and microelements determined by mean of ICP-MS in leaves and roots of 101.14 plants grown upon -Fe, +Fe, and FeBic conditions. Values followed by different letters are statistically different (ANOVA, *p* < 0.05).
cellular compartment plays a central role in the Fe and Mo interaction. We demonstrated that the accumulation of Mo in Fe-deficient plant might be related to the activity induction or protein accumulation of cytosolic-localized (such as nitrate reductase) and mitochondrial-localized (mARC) Mo-containing enzymes.

Moving up the leaf compartment in both the Fe impaired conditions (–Fe and FeBic) Fe showed a decrease in concentration coupled to an increase in Cu and Zn similarly to what observed in roots. The most interesting observations regarded leaves of plants growth under FeBic conditions, with Mn concentration significantly dropping compared to +Fe and –Fe, and, conversely Mo showing a remarkable increase in its concentration with respect to the other two conditions considered.

For what concerns the major macro-nutrients roots showed a reduced accumulation of phosphorus (P) and magnesium (Mg) compared to control, and this behavior was partially maintained also at the leaf level, although with less significant differences for what concerns Mg. Interestingly in leaves there was a huge drop in P concentration under FeBic growing conditions. A similar pattern of nutrients accumulation between the apical and basal organs of 101.14 plants was also observed for potassium (K), whose concentration was reduced only in –Fe conditions in both organs. Finally, of particular relevance was the increase in calcium (Ca) uptake at the root level. This was particularly evident in plants grown in presence of bicarbonate and was only partially correlated with measures obtained in leaves.

## mRNA-Seq Analyses of 101.14 Rootstock Genotype Subjected to Iron-Deficient and Iron-Limited Conditions

The whole transcriptomes of root apices collected from 101.14 plants hydroponically grown for 10 days in +Fe, –Fe, and FeBic were sequenced by mean of NGS technology using the SOLiD 5500xl platform. Recently, the genome resequencing of 101.14 indicated that, on average, there is a frequency of approximately one nucleotide variant (SNP) every 200 bases with respect to the PN40024 reference genome (Jaillon et al., 2007; Vitulo et al., 2014; Corso et al., 2015), suggesting the latter represents a suitable reference for read mapping. Thus, all mRNA-Seq reads that passed the quality control test (filtered based on read length after low-quality bases trimming) were mapped on the PN40024 12X v1 grape reference genome (http://genomes.cribi.unipd.it/grape/), producing, on average, ~21 million unique mapping reads per sample (ranging from 18 to 24 million reads depending on the sample). The Principal Component Analysis (PCA) plot (Supplementary Image 1) illustrating the correlation between biological replicates based on mRNA-Seq row filtered reads, showed a good correlation, although those replicates obtained from plants grown under –Fe appeared to be less closely related with respect to those collected from plants grown in +Fe and FeBic. Differentially expressed genes (DEGs) between stressed (–Fe and FeBic) and unstressed (+Fe) conditions are listed in Supplementary Data Sheets 1, 2. As a first general observation, the –Fe and FeBic exposure affected the expression

of a similar number of genes (365 DEGs against 331 DEGs, respectively). Eighty-three DE genes (13.5%) were common to both stress treatments, 282 DEGs were specifically regulated in root apices grown in –Fe and 248 in plants grown in presence of bicarbonate (Figure 3; Supplementary Data Sheet 3). In –Fe conditions there was a substantial balance between up and down regulated genes (180 up-regulated and 185 down regulated genes) whereas in FeBic conditions the majority of DE genes were down regulated (Figure 3). A similar trend was observed also for those DE genes that were common between –Fe and FeBic, with the majority of them being down-regulated with respect to the control (Figure 3; Supplementary Data Sheet 3). These results are in agreement with other studies (Thimm et al., 2001; Buckhout et al., 2009; Schmidt and Buckhout, 2011; Zamboni et al., 2012) showing that plant transcriptional response to Fe shortage is based on the modulation of a relative small set of genes compared to other stresses such as water deficit or high salinity (Corso et al., 2015).

## A Core Set of Genes Is Differentially Expressed in Both –Fe and FeBic Conditions and Can Be Related to the Arabidopsis Ferrome

The majority of DEGs in common between –Fe and FeBic showed the same behavior in both conditions compared to +Fe condition (Figure 4A). Nevertheless, the response observed in complete absence of Fe appeared to be more drastic, especially for stress-induced DE genes, which reached fold-changes (FC) often higher than FeBic (Supplementary Data Sheet 3). Figure 4B shows a manually curated distribution of DE genes in Gene Ontology (GO) categories. Several categories, such as those related to “Response to stress,” “Secondary metabolic process,” and “Transport” were clearly over-represented, with 15 (18%), 13 (16%), and 10 (12%) related-genes, respectively. Whereas the majority of genes belonging to the “Response to

<figure><img src="image_5.png" /><figcaption>**FIGURE 3** | Venn diagram illustrating Differentially Expressed Genes in grapevine roots in –Fe and FeBic conditions. Numbers of differentially regulated genes upon –Fe, FeBic or upon both –Fe and FeBic conditions are indicated, together with the relative contribute of up- (red) and down- (green) regulated genes to each subgroup (Fe-specific, common and FeBic-specific).</figcaption></figure>
<figure><img src="image_6.png" /></figure>

stress” and “Secondary metabolic process” were down-regulated (13 genes out of 15 in the “Response to stress” category and 8 genes out of 13 in the “Secondary metabolic process” one) in stressed roots compared to the control ones, the majority of genes belonging to the “Transport” category were induced (Figure 4B). Among them was VIT_10s0042g01100, encoding for IRT1 (IRON REGULATED TRANSPORTER 1), the main transporter involved in Fe²⁺ uptake in Strategy I plants (Lin et al., 2009). Together with IRT1-like transcript we identified other genes encoding for transporters such as VIT_00s0259g00140, VIT_00s0259g00120, and VIT_01s0026g01490, all annotated as oligopeptide transporters (OPTs), plasma membrane-embedded proteins, which are involved in the uptake of substrates from the apoplasm and the external environment and in their loading in the vascular system. In particular, VIT_00s0259g00120 is the orthologous of AtOPT3, which was demonstrated to transport transition ions *in vitro* and to be involved in phloem in Fe-loading, Fe recirculation from the xylem to the phloem, shoot-to-root Fe signaling, and Fe redistribution from mature to developing tissues (Stacey et al., 2006; Lubkowitz, 2011). VIT_05s0049g00930 was another transcript encoding for a transporter induced in both -Fe and FeBic conditions. In this case the deduced protein is an inorganic phosphate transporter whose role in iron-deficiency has not been described before to our knowledge. Not surprisingly all up-regulated genes

encoding for transporters can be associated to related genes in the so called *Arabidopsis* ferrome,” a core list of 95 transcripts which was identified by Schmidt and Buckhout (2011) and that robustly reflects the transcriptional response of *Arabidopsis* under Fe-deficient conditions. Together with them, many other genes belonging to other ontology categories were related to the ferrome. Among them are genes encoding for two Germin-like proteins (VIT_14s0060g00120 and VIT_14s0128g00650), a purple acid phosphatase (VIT_08s0007g05970), a zinc-finger protein (VIT_09s0054g01740) a cytochrome p450 gene (VIT_13s0067g00110), and most interesting, a bHLH transcription factor (TF), which is orthologous to the *Arabidopsis bHLH38*. Phylogenetic analyses (data not shown) on the amino acidic sequences of predicted bHLH members of the PN40024 genome (Jaillon et al., 2007) and *Arabidopsis* showed that there are three grapevine prediction which cluster with AtbHLH38 (AT3G56970) and AtbHLH39 (At3G56980): namely VIT_08s0007g05300, VIT_13s0073g00400, and VIT_13s0047g00450. Only VIT_08s0007g05300 was induced in both -Fe and FeBic. The other two were listed exclusively in -Fe specifically induced DE genes. In *Arabidopsis* bHLH38/39, together with other members of the bHLH family such as bHLH100 and bHLH101, act together with the binding partner FIT (FER-like Iron Deficiency Induced; Bauer et al., 2007), making heterodimers and activating transcription of
downstream iron-related genes such as *IRT1* and *FRO* and thus
inducing expression of the Strategy I iron-uptake machinery.
Interestingly *FIT* orthologous was not listed on DE genes
neither in –Fe and FeBic, being scarcely induced under stress
(FC = 1.3 in –Fe and 1.2 in FeBic). A possible explanation
of this observation is the fact that *FIT* expression is mainly
controlled post-transcriptionally as already demonstrated for its
orthologous gene in tomato (Brumbarova and Bauer, 2005).

Amongst those predictions associated to the *Arabidopsis*
ferrome were also transcripts undergoing down-regulation,
such as those encoding for ferritin 1 (VIT_13s0067g01840)
and vacuolar iron transporters-like 5 and 2 (VIT)
(VIT_12s0055g00860, VIT_18s0076g00320) proteins. Ferritins
are a class of ubiquitous proteins, generally located in chloroplast
and mitochondria, which are involved in iron storage and
chelation (Petit et al., 2001) and prevent the generation of ROS
due to Fe oxidation (Ravet et al., 2009). In the light of this,
their down-regulation is most likely due to the lower content
of Fe in the –Fe and FeBic samples. VIT (or nodulins-like)
are a group of proteins whose function is still object of debate,
although recently they have been described to be putatively
involved in Fe transport and storage under metal cation
sufficiency (Gollhofer et al., 2011). The down-regulation of
ferritin transcripts is in line with previous observations in roots
of Fe-depleted *Arabidopsis* plants (Thimm et al., 2001; Buckhout
et al., 2009; Yang et al., 2010; Stein and Waters, 2011) but also in
other species such as tomato (Zamboni et al., 2012). Similarly,
Nodulins were described to be down regulated under iron
shortage in *Arabidopsis* (Gollhofer et al., 2011), NODULIN-like
1 and 21 in particular (Stein and Waters, 2011). An interesting
aspect was the up-regulation of a nodulin-like gene in both –Fe
and FeBic conditions. Although this behavior was unexpected,
something similar was already observed in tomato roots, where
a positive modulation of two transcripts encoding for proteins
with a nodulin-like domain was described (Zamboni et al., 2012).

## Genes Specifically Expressed in Absence of Iron

DESeq2 R 0 package led to the identification of 282 genes specifically regulated under –Fe. As previously mentioned the behavior of most of –Fe specific DEGs was conserved also in FeBic, although with a weaker response compared to –Fe. This observation is clear looking at the heat map represented in **Figure 5A**, showing the behavior of –Fe specific DEGs in both conditions analyzed. A k-mean cluster analysis (**Figure 5B**) showed that among the six clusters that were identified, four were composed of genes showing the same behavior in both conditions (clusters Cl2/Cl3 and Cl5/Cl6 group genes respectively up and down regulated in response to both stresses). More intriguing are clusters Cl1 and Cl4, which contain several genes showing different patterns of expression between –Fe and FeBic. In most cases these genes were induced under –Fe and scarcely modulated or even repressed in FeBic conditions.

**Figure 5C** illustrates the relative contribute of up- and
down-regulated genes to the most-represented gene ontology
categories. Of particular interest was the category related to

“Secondary metabolic process,” which is the largest in size (up to 38 genes) and encompasses a large proportion of repressed genes in both –Fe and FeBic (Clusters Cl5 and Cl6). Among these, the majority belongs to the grapevine stilbene synthase (*VvSTS*) multigenic family, a large group of structural genes, which have been detected in a limited number of plant species and which are involved in the biosynthesis of resveratrol (three-hydroxy-*trans*-stilbene; Vannozzi et al., 2012; Parage et al., 2012), one of the major grapevine phytoalexins. The extended down-regulation of *VvSTSs* was also coupled to the repression of *VvMYB14* (VIT_07s0005g03340), an R2R3-MYB TF that was demonstrated to actively regulate the expression of *VvSTSs* (Höll et al., 2013) and to the up regulation of several genes belonging to the flavonoid branch, including two flavonoid 3′5′-hydroxylase (*F3′H*), a leucoanthocyanidin dioxygenase (*LDOX*) and two UDP-glucose: anthocyanidin 3-O-D-glucosyltransferase (*UGGT*; **Supplementary Data Sheet 3**). The opposite behavior between genes belonging to the flavonoid and stilbene competing pathways in response to stress was already reported by Vannozzi et al. (2012). In that study, it was observed that both downy mildew (*Plasmopara viticola*) infection and UV-C exposure led to an induction of *VvSTS* genes associated with the repression of three grapevine *VvCHSs* (chalcone synthase), the key genes leading to the flavonoid scaffold.

The fact that 101.14 rootstock appears to drive the flow of
carbon trough the synthesis of flavonoid rather than stilbenes
under iron deprivation stresses the role for these compounds
in the root response to iron deficiency. As a matter of fact, the
production of phenolics as root exudates is a well-documented
response of strategy I plants to Fe deficiency and include a wide
range of molecules, including organic acids, phenolics, flavins,
and flavonoids, among others (Schmidt, 1999; Dakora and
Phillips, 2002; Cesco et al., 2010; Abadía et al., 2011; Rodríguez-
Celma et al., 2013). The composition of these root exudates
appears to be genetically determined (Dakora and Phillips, 2002)
and it was speculated that excreted phenolic compounds could
contribute to the obligatory reduction of iron. Unfortunately, the
present study lacks a characterization and quantification of main
phenolic compounds secreted by 101.14 roots in the growing
media, but this certainly represents an interesting issue to focus
on in the future.

As observed for common DEGs, transport represents a
GO category particularly enriched in –Fe DEGs, with at
least 31 members. As observed for those DE genes encoding
for transporters in common with FeBic, also these –Fe
specific ones were induced in the majority of cases. Together
with VIT_10s0042g01100, we identified VIT_10s0042g01120
and VIT_03s0017g02170, orthologous to *AtIRT1* and *AtZIP12*.
Interestingly the latter was amongst those genes showing
opposite behavior between –Fe and FeBic, being repressed in –Fe
and induced in FeBic. Interestingly, *AtZIP12* was never described
as an induced gene in Fe-deprived roots but was identified as a
constitutively-expressed gene in the Cd/Zn hyper accumulator
*A. halleri* with respect to the non-tolerant *A. thaliana* (Roosens
et al., 2008). Many studies indicate that plant metal transporters
have wide substrate specificity with limited ion-selectivity and
probably the effective combinatorial control of metal transport
<figure><img src="image_7.png" /><figcaption>**FIGURE 5** | Heat-map of DE genes specifically regulated in roots grown in absence of iron. (A) The mean unique-read values of DE genes identified specifically in absence of iron were graphically represented by mean of MeV software. Hierarchical clustering was performed using MeV HCL tool with the Pearson correlation coefficient. Up-regulated genes are indicated in red whereas down-regulated genes are colored in green. (B) Expression clusters identified using the k-mean clustering using the Pearson correlation coefficient; (C) Main gene ontology categories related to DEGs according to VitisNet gene prediction update (Grimplet et al., 2012). For each category the contribute of up- (red) and -down (green) regulated genes is reported.</figcaption></figure>

in a particular deficiency situation must be seen as the result of the expression of a specific mix of transporters with limited selectivity together with high selective transporters such as *IRT1*. Moreover, post-transcriptional regulation should be taken into consideration when analyzing regulation of metal homeostasis, since this regulative process was previously described for *IRT1* in *Arabidopsis* (Connolly et al., 2002). VIT_18s0001g02130 and VIT_07s0129g00620 were other two transcripts showing significant induction under -Fe. They encode for proteins belonging to a different family of membrane-bound metal transporters known as NRAMPs (Natural Resistance-Associated Macrophage Protein; Forbes and Gros, 2001). In *Arabidopsis* AtNRAMP1, AtNRAMP3, and AtNRAMP4 encode functional plant metal transporters and are induced under low iron supply (Thomine et al., 2000). Although the specific role of these transporters in Fe homeostasis remains unclear, it appears likely that they localize to an intracellular compartment, probably to the plastids, as already confirmed for NRAMP1 by Curie et al. (2000). The induction of many aspecific ion transporters under iron shortage could explain the increase in concentration of other micronutrients that we observed thorough ICP-MS (Table 1).

Finally, among the other key genes composing the *Arabidopsis* ferrome we also identified VIT_17s0000g03440, the orthologous of *BRUTUS* (BTS; AT3G18290), which encodes a protein with a putative E3 ubiquitin-protein ligase domain and 6 haemerythrin cation binding motifs, which could potentially bind Fe. This protein, BTS, was shown to be co-regulated with PYE (POPEYE; Long et al., 2010), a newly identified TF which plays a critical role in Fe homeostasis and, although its mode of action has not been uncovered yet, it seems to act as a repressor of the iron deficiency stress response since roots of *bts* mutant plants were found to have a higher acidification response when compared to the wild-type and were less susceptible to Fe starvation. Interestingly we were not able to identify *POPEYE* among the DE genes.

## Genes Specifically Expressed in Plants Grown in Presence of Bicarbonate

Differing from what observed for Fe-specific DEGs, most of those genes differentially regulated in presence of bicarbonate did not show the same behavior under -Fe, confirming the hypothesis that the presence of bicarbonate leads to a specific response in root apices and should not be simply considered a "less severe" form of iron deprivation. **Figures 6A,B** clearly show that there
<figure><img src="image_8.png" /></figure>

are only a few genes, corresponding to those belonging to cluster Cl 2, which are repressed in both stress conditions (anyway the repression observed in -Fe is relatively low compared to FeBic). The remaining DE genes are divided in cluster Cl 1, composed of transcript induced only in FeBic conditions and cluster Cl 3, encompassing genes which are down-regulated only in FeBic. This observation deserves a thought, since many studies on plant response to iron deficiency are based on experiments aimed at mimicking the field conditions by subtracting iron from the growing media. Based on our observations the roots response to a total absence of Fe from the nutritive solution is only partially overlapping with the response observed in presence of Fe plus bicarbonate, and considering the latter condition is more similar to what happens in the field, we believe future research should be oriented at unrevealing the biologic processes involved in the plant response to this particular condition rather than to a complete Fe deprivation. As observed for common and -Fe specific DEGs, also in this case the most represented GO categories were related to "Secondary metabolic Processes," "Transport" and "transcription factor activity (Figure 6C)."

Among the secondary metabolic processes, several genes encoding for enzyme related to the coumarins synthesis

were down regulated (i.e., ferulate 5-hydroxylase; Cinnamoyl alcohol dehydrogenase). As recently demonstrated by several authors, coumarins represent important compound that plants root released into rhizosphere to implement Fe solubility and availability. However, such effect of coumarins has been demonstrated only in Arabidopsis plants, while other plants, such as *Medicago trucantula* are able to release out of the root mainly flavins instead of coumarins (Rodríguez-Celma et al., 2013). Such evidence suggested that the root exudates profiles of Fe-deficient plants might be species-specific. Our result might suggest that response to low Fe availability in *Vitis* spp. might not involve coumarins. Indeed, some other genes involved in the flavonoid biosynthetic pathways (VIT_00s0555g00020, F3'5'H; VIT_07s0005g03150, F3H-like) were up regulated under +FeBic suggesting that other compounds such as myricetin might be likely synthetized in +FeBic-treated plants. Other than secondary metabolites, root exudates released by Fe-deficient plants might involve also carboxylate compounds, such as citrate. It is known that carboxylates are actively released by plants through secondary active transport mediated by protein carriers at the plasma membrane, and the exudation of citrate most probably relies on the activity of MATE transporters (Magalhaes et al.,
2007; Valentinuzzi et al., 2015). Interestingly, in the Transport GO category we observed a down regulation of three genes encoding for MATE transporters (Supplementary Data Sheet 3) in +FeBic treated plants, suggesting that probably in 101.14, also the exudation of citrate should not be induced under low Fe availability caused by the presence of bicarbonate. Nevertheless, it's worth to note that MATE transporters are not only involved in the citrate transport but they have a vast range of roles in the mobilization of many other phenolics.

The decrease of both coumarins-synthesis related genes and MATE transporter are more pronounced in +FeBic than in -Fe, suggesting that the exudation of root compounds likely changes when bicarbonate is present in the external medium. Interestingly, with respect to the -Fe, +FeBic treatment strongly induce the expression of some genes related to phosphate homeostasis: VTI_11s0016g05330 (SPX2), VIT_07s0005g03290 (Inorganic phosphate transporter 1-4), and VIT_01s0011g06290 (Purple acid phosphatase 3). It is known that calcareous soil limit P availability for plant, therefore 101.14 plant growing under +FeBic by facing also to P deficiency might induce some P-responsive genes.

## Comparison between 101.14 and M1 a New Experimental Tolerant Genotype

The expression of set of selected genes, which were differentially expressed in our mRNA-Seq analyses, was checked by quantitative RT-PCR (qPCR) in a third biological replicate constituted by and experiment performed independently from those used for NGS whole transcriptome analyses. We also included the expression of genes which surprisingly were not listed within the DE genes but that are known to be part of the Arabidopsis ferrome and to be involved in the Strategy I mechanism of iron uptake such as *PYE* (*POPEYE*), *AHA2* and *FRO2*. In this new experimental trial, stresses were imposed not only to 101.14 plants, but also to another rootstock genotype, named M1. This genotype is characterized by a high degree of tolerance to iron chlorosis, high performances to grafting, reduced vigor and medium resistance to salinity and was recently registered to the National Register of grapevine varieties (http://catalogoviti.politicheagricole.it/crediti.php). Differently from 101.14, the physiological data obtained in M1 under Fe-limited availability pointed out that *Pn*, *GS* and *E* are similar to control (Supplementary Image 2) without evident symptoms of chlorosis. **Figure 7** shows the normalized transcript level of nine iron-related genes including *IRT1*, a *FERRITIN1*-like gene, *bHLH38/39*, a *NODULIN*-like gene, and *VvSTS41* in both the genotypes. For what concerns 101.14, results confirmed what observed in mRNA-Seq analyses, but also highlighted several differences between the two genotypes under study. For example, in both genotypes the *IRT1* transcript was strongly induced in both -Fe and FeBic conditions, but the expression in +Fe condition in M1 genotype was almost six times higher than that observed in 101.14. Also the induction registered under stress was stronger than 101.14. This observation, together with the stronger expression *FERRITIN*-like gene in M1 genotype with respect to 101.14, could be related to a higher capacity of M1

to uptake and to store iron already in unstressed condition. All these results are also supported by the biochemical analysis that pointed out a higher accumulation of Fe under FeBic conditions in M1 than in 101.14. Looking at the expression of genes encoding for transcriptional regulators such as *BRUTUS* and *POPEYE*, both the genotypes showed similar behaviors, although M1 showed a stronger induction of *PYE* in -Fe. Also stilbene synthase and leucoanthocyanidin oxidase genes (*VvSTS41* & *LDOX*) confirmed what observed in the mRNA-Seq data for what concerns 101.14 and highlighted some differences between the two rootstock genotypes under study, at least in terms of entity of the response, although the patterns of expression were comparable. In fact, in both genotypes *VvSTS41* was down regulated under -Fe and was coupled to a remarkable up-regulation of *LDOX* confirming that there must be some involvement of the flavonoid pathway in the plant adaptation to iron-impaired conditions.

An aspect of particular interest regards the transcript accumulation of genes encoding for *AHA2* and *FRO2*, two key genes involved in the Strategy I iron uptake mechanism. These genes were considered separately from the other ones analyzed by qPCR since we associated to their expression (Figure 8A) also other biochemical analyses such as the measure of the enzymatic activity (Figure 8B) and the evaluation of the root ability to acidify the substrate and to reduce Fe in agar inclusions (Figure 8C). The first interesting observation is the fact that the expression of *AHA2* in M1 genotype is much higher compared to 101.14 in all conditions considered, suggesting a better ability of this genotype to extrude protons in the apoplast. This observation was confirmed at the level of enzymatic activity, which ranged between 0.6 and 0.9 µmol NADH min⁻¹ mg⁻¹ prot in M1 (in +Fe and -Fe, respectively) and between 0.1 and 0.15 µmol NADH min⁻¹ mg⁻¹ prot in 101.14 and by the higher capability of M1 to acidify the rhizosphere. All together these observations suggest that the typical strategy I responses are strongly active in M1 while are weak, more delocalized, and diffused in 101.14.

Something similar happened for *FRO2*, the gene encoding for an iron deficiency-inducible ferric chelate reductase responsible for reduction of iron at the root surface. Also in this case the measure of the enzyme activity and the *in vivo* evaluation measured using BPDS confirmed a higher activity of M1 genotype compared to 101.14. Nevertheless, the correlation between biochemical and molecular data was less clear since the expression of *FRO2* gene showed a different pattern respect to the enzyme activity, with a total absence of transcript in M1 in FeBic. A possible explanation relies on the fact that, although focusing on VIT_15s0046g01900, which was the closest *AtFRO2* orthologous in grapevine, we identified at least other three predictions, which cluster with *AtFRO2* (Atg01580). Thus, the effect observed in the biochemical analyses could be due to the synergic action of all or part of these enzymes although at the transcriptional level they could show different patterns of expression. In other words, the activity observed in FeBic conditions could be due to the action of other grapevine *AtFRO2* orthologs with a higher level of expression than VIT_15s0046g01900.
<figure><img src="image_9.png" /><figcaption>FIGURE 7 | Real-time RT-PCR validation of a set of genes differentially expressed in mRNA-Seq analyses in 101.14 and in an experimental tolerant rootstock (M1). A subset of DE genes was screened by qPCR upon -Fe, +Fe, and FeBic conditions both in 101.14 M1 genotypes. Transcript were normalized to the expression of elongation factor (EF1) and plotted as normalized transcript level. Bars indicate standard error (SE) in three technical replicates. Empty columns represent 101.14 genotype whereas dark gray bars represent M1. The gene considered in qPCR analyses were LDOX (VIT_06s0004g00770), FERRITIN (VIT_13s0067g01840), NODULIN (VIT_19s0093g00320), AHA2 (VIT_11s0052g00620, IRT1 (VIT_10s0042g01100), FRO2 (VIT_15s0046g01900), bHLH38/39 (VIT_08s0007g05300), POPEYE (VIT_17s0000g06000), BRUTUS (VIT_17s0000g03440).</figcaption></figure>

CONCLUSIONS

Low iron bioavailability is a primary constraint to plant growth in many ecosystems, particularly in calcareous soils commonly found in the Mediterranean basin. Thus, iron chlorosis is one of the most frequent nutritional problems in grapevine production. The present study was performed in order to provide a first general description of the transcriptome profile in root apices of a commonly used rootstock genotype under iron deficiency. Not only we analyzed the transcriptome effect of iron deprivation itself, but we also compared it to a condition of low iron availability determined by high concentration of bicarbonate, a phenomenon quite common in nature, which better mimics what happens in the vineyard. To our knowledge this is the first

study based on a whole transcriptome analysis describing the
response of a grapevine rootstock to iron deficiency. Thanks to
this approach we identified many genes potentially involved in
the Fe homeostasis in grapevine, some of which differentially
modulated in response to both -Fe and FeBic conditions, some
others specific to a particular one, suggesting that the root
responses to Fe absence and presence of bicarbonate are only
partially overlapping and could induce different mechanisms
related to different gene regulatory networks. It is important to
consider that FeBic treatment, by mimicking a calcareous-like
condition, would affect the availability of other nutrients, such as
P. Therefore, the partial overlap of results between -Fe and FeBic
treatments provides the specific Fe deficiency-induced responses
of grape rootstock.
A

<figure><img src="image_10.png" /></figure>

B

<figure><img src="image_11.png" /></figure>

C

<figure><img src="image_12.png" /></figure>

Acididification

<figure><img src="image_13.png" /></figure>

FRO2

<figure><img src="image_14.png" /></figure>

FRO2 activity

<figure><img src="image_15.png" /></figure>

101.14

<figure><img src="image_16.png" /></figure>

Reduction

M1

-Fe +Fe FeBic

<caption>FIGURE 8 | Gene expression and enzymatic activity of AHA2 and FRO2 in 101.14 and M1 genotypes upon -Fe, +Fe and FeBic conditions. (A)</caption>

Expression analyses of AHA2 and FRO2 transcript in both 101.14 and M1 rootstocks; (B) measure of the H+/ATPase and FeIII-chelate reductase activity; (C) evaluation of the root ability to acidify the substrate and reduce iron in agar inclusions.
Many DEGs were orthologs of the so-called Arabidopsis ferrome, confirming the grapevine roots behave similarly to Arabidopsis under iron shortage undertaking a Strategy I mechanism of iron homeostasis and representing candidate targets for more focused functional characterization studies. Our results indicated that, together with genes involved in transport, transcriptional regulation, and response to stress, those one related to secondary metabolism are strongly affected by the stresses, as observed for the biosynthetic pathways related to stilbene biosynthesis, which was strongly impaired by stress, and the flavonoid one, some of whose genes were highly induced. Although the 101.14 rootstock genotype is characterized by a high susceptibility to iron chlorosis, the induction of several genes related to the biosynthesis of phenols under stress highlights, once again, the transversal role of these class of compounds in the plant response and adaptation to a vast range of conditions (Ramakrishna and Ravishankar, 2011). These compounds, whose chemical nature has not been described in this study, could be accumulated in the roots with the scope to increase the solubility of iron and its mobilization from insoluble reservoirs and could represent a determinant of tolerance in genotypes characterized by different level of susceptibility to Fe chlorosis. From an agronomical point of view, one of the most interesting outcomes of our study is the different behavior observed for genes involved in the acidification and Fe-reduction processes between 101.14 and M1, another rootstock genotype showing high tolerance to iron chlorosis. In these two genotypes the capability of modulating the strategy I response at the root level, well-correlate with the manifestation of symptoms at the leaf level, suggesting that the activity of AHA2 and FRO2, possibly in association with the capability to accumulate phenolics and transcripts of their biosynthetic genes, could be candidate parameters to be used as functional markers for the selection of new tolerant genotypes.

## AUTHOR CONTRIBUTIONS

AV and ML conceived the design of this study and participated in its coordination; AV performed the bioinformatic data analysis

and wrote the paper; SD planned and conducted most of the lab experiments; MC strongly contributed to the statistical analyses of NGS data and in writing the manuscript; GVi and GZ contributed to the interpretation of results and in writing of manuscript; GVa and NV provided NGS platforms and contributed to preliminary bioinformatic analyses; CB and ML contributed to the writing of the paper and critically revised the manuscript. All authors have read and approved the final manuscript.

## FUNDING

This study was supported by the AGER "SERRES" project, grant n° 2010-2015.

## ACKNOWLEDGMENTS

Authors would like to thank Dr. Patrick Guerretta and Dr. Igor Zen for their technical contribution in the experimental procedures and in data analyses.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2016.01994/full#supplementary-material

<caption>Supplementary Image 1 | PCA plot of mRNA-Seq unique reads.</caption>

<caption>Supplementary Image 2 | M1 genotype physiological data upon –Fe, +Fe, and FeBic conditions.</caption>

<caption>Supplementary Data Sheet 1 | List of Differentially Expressed Genes identified in root apices of 101.14 rootstock genotype upon –Fe growing conditions.</caption>

<caption>Supplementary Data Sheet 2 | List of Differentially Expressed Genes identified in root apices of 101.14 rootstock genotype upon FeBic growing conditions.</caption>

<caption>Supplementary Data Sheet 3 | List of common (–Fe + FeBic), Fe-specific and FeBic-specific Differentially Expressed Genes identified in root apices of 101.14 rootstock genotype upon –Fe and FeBic growing conditions.</caption>

## REFERENCES

Abadía, J., Vázquez, S., Rellán-Álvarez, R., El-Jendoubi, H., Abadía, A., Álvarez-Fernández, A., et al. (2011). Towards a knowledge-based correction of iron chlorosis. *Plant Physiol. Biochem.* 49, 471–482. doi: 10.1016/j.plaphy.2011.01.026

Arrigo, N., and Arnold, C. (2007). Naturalised vitis rootstocks in Europe and consequences to native wild grapevine. *PLoS ONE* 2:e521. doi: 10.1371/journal.pone.0000521

Balk, J., and Schaedler, T. A. (2014). Iron cofactor assembly in plants. *Annu. Rev. Plant Biol.* 65, 125–153. doi: 10.1146/annurev-arplant-050213-035759

Bauer, P., Ling, H. Q., and Guerinot, M. L. (2007). FIT, the FER-LIKE IRON DEFICIENCY INDUCED TRANSCRIPTION FACTOR in Arabidopsis. *Plant Physiol. Biochem.* 45, 260–261. doi: 10.1016/j.plaphy.2007.03.006

Bavaresco, L., Bertamini, M., and Iacono, F. (2006). Lime-induced chlorosis and physiological responses in grapevine (*Vitis vinifera* L. cv. Pinot blanc) leaves. 45, 45–46.

Bavaresco, L., Fraschini, P., and Perino, A. (1993). Effect of the rootstock on the occurrence of lime-induced chlorosis of potted *Vitis vinifera* L. cv. 'Pinot blanc'. *Plant Soil* 157, 305–311. doi: 10.1007/BF00011058

Bertamini, M., Muthuchelian, K., and Nedunchezhian, N. (2002). Iron deficiency induced changes on the donor side of PS II in field grown grapevine (*Vitis vinifera* L. cv. Pinot noir) leaves. *Plant Sci.* 162, 599–605. doi: 10.1016/S0168-9452(01)00604-5

Brancadoro, L., Rabotti, G., Scienza, A., and Zocchi, G. (1995). Mechanisms of Fe-efficiency in roots of *Vitis* spp. in response to iron deficiency stress. *Plant Soil* 171, 229–234. doi: 10.1007/BF00010276

Briat, J.-F., Dubos, C., and Gaymard, F. (2015). Iron nutrition, biomass production, and plant product quality. *Trends Plant Sci.* 20, 33–40. doi: 10.1016/j.tplants.2014.07.005

Brumbarova, T., and Bauer, P. (2005). Iron-mediated control of the basic helix-loop-helix protein FER, a regulator of iron uptake in tomato. *Plant Physiol.* 137, 1018–1026. doi: 10.1104/pp.104.054270

Brumbarova, T., Matros, A., Mock, H. P., and Bauer, P. (2008). A proteomic study showing differential regulation of stress, redox regulation and peroxidase
proteins by iron supply and the transcription factor FER. *Plant J.* 54, 321–334.
doi: 10.1111/j.1365-313X.2008.03421.x

Buckhout, T. J., Yang, T. J., and Schmidt, W. (2009). Early iron-deficiency-induced transcriptional changes in Arabidopsis roots as revealed by microarray analyses. *BMC Genomics* 10:147. doi: 10.1186/1471-2164-10-147

Campagna, D., Albiero, A., Bilardi, A., Caniato, E., Forcato, C., Manavski, S., et al. (2009). PASS: a program to align short sequences. *Bioinformatics* 25, 967–968. doi: 10.1093/bioinformatics/btp087

Cerovic, Z. G., Masdoumier, G., Ghozlen, N. B., and Latouche, G. (2012). A new optical leaf-clip meter for simultaneous non-destructive assessment of leaf chlorophyll and epidermal flavonoids. *Physiol. Plant.* 146, 251–260. doi: 10.1111/j.1399-3054.2012.01639.x

Cesco, S., Neumann, G., Tomasi, N., Pinton, R., and Weisskopf, L. (2010). Release of plant-borne flavonoids into the rhizosphere and their role in plant nutrition. *Plant Soil* 329, 1–25. doi: 10.1007/s11104-009-0266-9

Connolly, E. L., Fett, J. P., and Guerinot, M. L. (2002). Expression of the IRT1 metal transporter is controlled by metals at the levels of transcript and protein accumulation. *Plant Cell* 14, 1347–1357. doi: 10.1105/tpc.001263

Corso, M., Vannozzi, A., Maza, E., Vitulo, N., Meggio, F., Pitacco, A., et al. (2015). Comprehensive transcript profiling of two grapevine rootstock genotypes contrasting in drought susceptibility links the phenylpropanoid pathway to enhanced tolerance. *J. Exp. Bot.* 66, 5739–5752. doi: 10.1093/jxb/erv274

Corso, M., Vannozzi, A., Ziliotto, F., Zouine, M., Maza, E., Nicolato, T., et al. (2016). Grapevine rootstocks differentially affect the rate of ripening and modulate auxin-related genes in cabernet sauvignon berries. *Front. Plant Sci.* 7:69. doi: 10.3389/fpls.2016.00069

Curie, C., Alonso, J. M., Le Jean, M., Ecker, J. R., and Briat, J. F. (2000). Involvement of NRAMP1 from *Arabidopsis thaliana* in iron transport. *Biochem. J.* 347 (Pt 3),
749–55. doi: 10.1042/bj3470749

Dakora, F. D., and Phillips, D. A. (2002). Root exudates as mediators of mineral acquisition in low-nutrient environments. *Plant Soil* 245, 35–47. doi: 10.1023/A:1020809400075

Dell'Orto, M., Pirovano, L., Villalba, J. M., González Reyes, J. A., and Zocchi, G. (2002). Localization of the plasma membrane H+-ATPase in Fe-deficient cucumber roots by immunodetection. *Plant Soil* 241, 11–17. doi: 10.1023/A:1016030514200

Donnini, S., de Nisi, P., Gabotti, D., Tato, L., and Zocchi, G. (2012). Adaptive strategies of *Parietaria diffusa* (M.&K.) to calcareous habitat with limited iron availability. *Plant Cell Environ.* 35, 1171–1184. doi: 10.1111/j.1365-3040.2012.02481.x

Donnini, S., Prinsi, B., Negri, A. S., Vigani, G., Espen, L., and Zocchi, G. (2010). Proteomic characterization of iron deficiency responses in *Cucumis sativus* L. roots. *BMC Plant Biol.* 10:268. doi: 10.1186/1471-2229-10-268

Eide, D., Broderius, M., Fett, J., and Guerinot, M. L. (1996). A novel iron-regulated metal transporter from plants identified by functional expression in yeast. *Proc. Natl. Acad. Sci. U.S.A.* 93, 5624–5628. doi: 10.1073/pnas.93.11.5624

Fisarakis, I., Chartzoulakis, K., and Stavrakas, D. (2001). Response of Sultana vines (*V. vinifera* L.) on six rootstocks to NaCl salinity exposure and recovery. *Agric. Water Manage.* 51, 13–27. doi: 10.1016/S0378-3774(01)00115-9

Forbes, J. R., and Gros, P. (2001). Divalent-metal transport by NRAMP proteins at the interface of host-pathogen interactions. *Trends Microbiol.* 9, 397–403. doi: 10.1016/S0966-842X(01)02098-4

Forner-Giner, M. A., Llosá, M. J., Carrasco, J. L., Perez-Amador, M. A., Navarro, L., and Ancillo, G. (2010). Differential gene expression analysis provides new insights into the molecular basis of iron deficiency stress response in the citrus rootstock *Poncirus trifoliata* (L.) Raf. *J. Exp. Bot.* 61, 483–490. doi: 10.1093/jxb/erp328

Fourcroy, P., Sisó-Terraza, P., Sudre, D., Savirón, M., Reyt, G., Gaymard, F., et al. (2014). Involvement of the ABCG37 transporter in secretion of scopoletin and derivatives by Arabidopsis roots in response to iron deficiency. *New Phytol.* 201, 155–167. doi: 10.1111/nph.12471

Gambetta, G. A., Manuck, C. M., Drucker, S. T., Shaghasi, T., Fort, K., Matthews, M. A., et al. (2012). The relationship between root hydraulics and scion vigour across Vitis rootstocks: what role do root aquaporins play? *J. Exp. Bot.* 63, 6445–6455. doi: 10.1093/jxb/ers312

Gollhofer, J., Schlwicke, C., Jungnick, N., Schmidt, W., and Buckhout, T. J. (2011). Members of a small family of nodulin-like genes are regulated under iron

deficiency in roots of *Arabidopsis thaliana*. *Plant Physiol. Biochem.* 49, 557–564.
doi: 10.1016/j.plaphy.2011.02.011

Gregory, P. J., Atkinson, C. J., Bengough, A. G., Else, M. A., Fernández-Fernández, F., Harrison, R. J., et al. (2013). Contributions of roots and rootstocks to sustainable, intensified crop production. *J. Exp. Bot.* 64, 1209–1222. doi: 10.1093/jxb/ers385

Grimplet, J., Van Hemert, J., Carbonell-Bejerano, P., Díaz-Riquelme, J., Dickerson, J., Fennell, A., et al. (2012). Comparative analysis of grapevine whole-genome gene predictions, functional annotation, categorization and integration of the predicted gene sequences. *BMC Res.* 5:213. doi: 10.1186/1756-0500-5-213

Guerinot, M. L. (2000). The ZIP family of metal transporters. *Biochim. Biophys. Acta Biomembr.* 1465, 190–198. doi: 10.1016/S0005-2736(00)00138-3

Hänsch, R., and Mendel, R. R. (2009). Physiological functions of mineral micronutrients (Cu, Zn, Mn, Fe, Ni, Mo, B, Cl). *Curr. Opin. Plant Biol.* 12, 259–266. doi: 10.1016/j.pbi.2009.05.006

Höll, J., Vannozzi, A., Czemmel, S., D'Onofrio, C., Walker, A. R., Rausch, T., et al. (2013). The R2R3-MYB transcription factors MYB14 and MYB15 regulate stilbene biosynthesis in *Vitis vinifera*. *Plant Cell* 25, 4135–4149. doi: 10.1105/tpc.113.117127

Jaillon, O., Aury, J.-M., Noel, B., Policriti, A., Clepet, C., Casagrande, A., et al. (2007). The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. *Nature* 449, 463–467. doi: 10.1038/nature06148

Jin, C. W., Du, S. T., Chen, W. W., Li, G. X., Zhang, Y. S., and Zheng, S. J. (2009). Elevated carbon dioxide improves plant iron nutrition through enhancing the iron-deficiency-induced responses under iron-limited conditions in tomato. *Plant Physiol.* 150, 272–280. doi: 10.1104/pp.109.136721

Jin, C. W., Liu, Y., Mao, Q. Q., Wang, Q., and Du, S. T. (2013). Mild Fe-deficiency improves biomass production and quality of hydroponic-cultivated spinach plants (*Spinacia oleracea* L.). *Food Chem.* 138, 2188–2194. doi: 10.1016/j.foodchem.2012.12.025

Jin, C. W., You, G. Y., He, Y. F., Tang, C., Wu, P., and Zheng, S. J. (2007). Iron deficiency-induced secretion of phenolics facilitates the reutilization of root apoplastic iron in red clover. *Plant Physiol.* 144, 278–285. doi: 10.1104/pp.107.095794

Kobayashi, T., and Nishizawa, N. K. (2012). Iron uptake, translocation, and regulation in higher plants. *Annu. Rev. Plant Biol.* 63, 131–152. doi: 10.1146/annurev-arplant-042811-105522

Lan, P., Li, W., Wen, T.-N., Shiau, J.-Y., Wu, Y.-C., Lin, W., et al. (2011). iTRAQ protein profile analysis of Arabidopsis roots reveals new aspects critical for iron homeostasis. *Plant Physiol.* 155, 821–834. doi: 10.1104/pp.110.169508

Li, J., Wu, X. D., Hao, S. T., Wang, X. J., and Ling, H. Q. (2008). Proteomic response to iron deficiency in tomato root. *Proteomics* 8, 2299–2311. doi: 10.1002/pmic.200700942

Lin, Y. F., Liang, H. M., Yang, S. Y., Boch, A., Clemens, S., Chen, C. C., et al. (2009). Arabidopsis IRT3 is a zinc-regulated and plasma membrane localized zinc/iron transporter. *New Phytol.* 182, 392–404. doi: 10.1111/j.1469-8137.2009.02766.x

Long, T. A., Tsukagoshi, H., Busch, W., Lahner, B., Salt, D. E., and Benfey, P. N. (2010). The bHLH transcription factor POPEYE regulates response to iron deficiency in Arabidopsis roots. *Plant Cell* 22, 2219–2236. doi: 10.1105/tpc.110.074096

Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. *Genome Biol.* 15, 1–34. doi: 10.1186/s13059-014-0550-8

Lukkowitz, M. (2011). The oligopeptide transporters: a small gene family with a diverse group of substrates and functions? *Mol. Plant* 4, 407–415. doi: 10.1093/mp/ssr004

Magalhaes, J. V, Liu, J., Guimarães, C. T., Lana, U. G. P., Alves, V. M. C., Wang, Y.-H., et al. (2007). A gene in the multidrug and toxic compound extrusion (MATE) family confers aluminum tolerance in sorghum. *Nat. Genet.* 39, 1156–1161. doi: 10.1038/ng2074

Marguerit, E., Brendel, O., Lebon, E., Van Leeuwen, C., and Ollat, N. (2012). Rootstock control of scion transpiration and its acclimation to water deficit are controlled by different genes. *New Phytol.* 194, 416–429. doi: 10.1111/j.1469-8137.2012.04059.x

Meggio, F., Prinsi, B., Negri, A. S., Simone Di Lorenzo, G., Lucchini, G., Pitacco, A., Failla, O., et al. (2014). Biochemical and physiological responses of two grapevine rootstock genotypes to drought and salt treatments. *Aust. J. Grape Wine Res.* 20, 310–323. doi: 10.1111/ajgw.12071
Nikolic, M., Römheld, V., and Merkt, N. (2000). Effect of bicarbonate on uptake and translocation of 59Fe in two grapevine rootstocks differing in their resistance to Fe deficiency chlorosis. *Vitis* 39, 145-149.

O'Rourke, J. A., Charlson, D. V., Gonzalez, D. O., Vodkin, L. O., Graham, M. A., Cianzio, S. R., et al. (2007a). Microarray analysis of iron deficiency chlorosis in near-isogenic soybean lines. *BMC Genomics* 8:476. doi: 10.1186/1471-2164-8-476

O'Rourke, J. A., Graham, M. A., Vodkin, L., Gonzalez, D. O., Cianzio, S. R., and Shoemaker, R. C. (2007b). Recovering from iron deficiency chlorosis in near-isogenic soybeans: a microarray study. *Plant Physiol. Biochem.* 45, 287-292. doi: 10.1016/j.plaphy.2007.03.008

O'Rourke, J. A., Nelson, R. T., Grant, D., Schmutz, J., Grimwood, J., Cannon, S., et al. (2009). Integrating microarray analysis and the soybean genome to understand the soybeans iron deficiency response. *BMC Genomics* 10:376. doi: 10.1186/1471-2164-10-376

Parage, C., Tavares, R., Réty, S., Baltenweck-Guyot, R., Poutaraud, A., Renault, L., et al. (2012). Structural, functional, and evolutionary analysis of the unusually large stilbene synthase gene family in grapevine. *Plant Physiol.* 160, 1407-1419. doi: 10.1104/pp.112.202705

Petit, J.-M., Briat, J.-F., and Lobréaux, S. (2001). Structure and differential expression of the four members of the *Arabidopsis thaliana* ferritin gene family. *Biochem. J.* 359:575. doi: 10.1042/bj3590575

Piagnani, C., De Nisi, P., Espen, L., and Zocchi, G. (2003). Adaptive responses to iron-deficiency in callus cultures of two cultivars of *Vitis* spp. *J. Plant Physiol.* 160, 865-870. doi: 10.1078/0176-1617-00993

Ramakrishna, A., and Ravishankar, G. A. (2011). Influence of abiotic stress signals on secondary metabolites in plants. *Plant Signal. Behav.* 6, 1720-1731. doi: 10.4161/psb.6.11.17613

Ravet, K., Touraine, B., Boucherez, J., Briat, J. F., Gaymard, F., and Cellier, F. (2009). Ferritins control interaction between iron homeostasis and oxidative stress in Arabidopsis. *Plant J.* 57, 400-412. doi: 10.1111/j.1365-313X.2008.03698.x

Rellán-Álvarez, R., Andaluz, S., López-Millán, A., Fiehn, O., Álvarez-Fernández, A., and Abadia, J. (2009). "Changes in the proteomic and metabolomic profiles of beta vulgaris root tips in response to iron deficiency and resupply," in *Proc Int Plant Nutr Colloq XVI*, University California, Davis.

Robinson, N. J., Procter, C. M., Connolly, E. L., and Guerinot, M. L. (1999). A ferric-chelate reductase for iron uptake from soils. *Nature* 397, 694-697. doi: 10.1038/17800

Rodríguez-Celma, J. Lin, W. D., Fu, G. M., Abadia, J., López-Millán, A. F., and Schmidt W (2013). Mutually exclusive alterations in secondary metabolism are critical for the uptake of insoluble iron compounds by Arabidopsis and *Medicago truncatula*. *Plant Physiol.* 162, 1473-1485. doi: 10.1104/pp.113.220426

Rolli, E., Marasco, R., Vigani, G., Ettoumi, B., Mapelli, F., Deangelis, M. L., et al. (2015). Improved plant resistance to drought is promoted by the root-associated microbiome as a water stress-dependent trait. *Environ. Microbiol.* 17, 316-331. doi: 10.1111/1462-2920.12439

Römheld, V., and Marschner, H. (1986). Evidence for a specific uptake system for iron phytosiderophores in roots of grasses. *Plant Physiol.* 80, 175-180. doi: 10.1104/pp.80.1.175

Roosens, N. H. Willems, G., and Saumitou-Laprade, P. (2008). Using Arabidopsis to explore zinc tolerance and hyperaccumulation. *Trend Plant Sci.* 13, 208-215. doi: 10.1016/j.tplants.2008.02.006

Russo, M. A., Sambuco, F., and Belligno, A. (2010). The response to iron deficiency of two sensitive grapevine cultivars grafted on a tolerant rootstock. *Afr. J. Biochem. Res.* 4, 33-42.

Saeed, A. I., Bhagabati, N. K., Braisted, J. C., Liang, W., Sharov, V., Howe, E. A., et al. (2006). TM4 Microarray Software Suite. *Methods Enzymol.* 411, 134-193. doi: 10.1016/S0076-6879(06)11009-5

Santi, S., and Schmidt, W. (2009). Dissecting iron deficiency-induced proton extrusion in Arabidopsis roots. *New Phytol.* 183, 1072-1084. doi: 10.1111/j.1469-8137.2009.02908.x

Santi, S., Cesco, S., Varanini, Z., and Pinton, R. (2005). Two plasma membrane H+-ATPase genes are differentially expressed in iron-deficient cucumber plants. *Plant Physiol. Biochem.* 43, 287-292. doi: 10.1016/j.plaphy.2005.02.007

Schmidt, H., Günther, C., Weber, M., Spörlein, C., Loscher, S., Böttcher, C., et al. (2014). Metabolome analysis of *Arabidopsis thaliana* roots identifies

a key metabolic pathway for iron acquisition. *PLoS ONE* 9:e0102444. doi: 10.1371/journal.pone.0102444

Schmidt, W. (1999). Mechanisms and regulation of reduction-based iron uptake in plants. *New Phytol.* 141, 1-26. doi: 10.1046/j.1469-8137.1999.00331.x

Schmidt, W., and Buckhout, T. J. (2011). A hitchhiker's guide to the Arabidopsis ferrome. *Plant Physiol. Biochem.* 49, 462-470. doi: 10.1016/j.plaphy.2010.12.001

Stacey, M. G., Osawa, H., Patel, A., Gassmann, W., and Stacey, G. (2006). Expression analyses of Arabidopsis oligopeptide transporters during seed germination, vegetative growth and reproduction. *Planta* 223, 291-305. doi: 10.1007/s00425-005-0087-x

Stein, R. J., and Waters, B. (2011). Use of natural variation reveals core genes in the transcriptome of iron-deficient *Arabidopsis thaliana* roots. *J. Exp. Bot.* 63, 1039-1055. doi: 10.1093/jxb/err343

Tagliavini, M., and Rombolà, A. D. (2001). Iron deficiency and chlorosis in orchard and vineyard ecosystems. *Eur. J. Agron.* 15, 71-92. doi: 10.1016/S1161-0301(01)00125-3

Takahashi, M., Nakanishi, H., Kawasaki, S., Nishizawa, N. K., and Mori, S. (2001). Enhanced tolerance of rice to low iron availability in alkaline soils using barley nicotianamine aminotransferase genes. *Nat Biotechnol.* 19, 466-469. doi: 10.1038/88143

Thimm, O., Essigmann, B., Kloska, S., Altmann, T., and Buckhout, T. J. (2001). Response of Arabidopsis to iron deficiency stress as revealed by microarray analysis. *Plant Physiol.* 127, 1030-1043. doi: 10.1104/pp.010191

Thomine, S., and Lanquar, V. (2011). "Iron transport and signaling in plants" in *Transporters and Pumps in Plant Signaling*, eds M. Geisler, and K. Venema, (Berlin; Heidelberg: Springer-Verlag), 99-131. doi: 10.1007/978-3-642-14369-4_4

Thomine, S., Wang, R., Ward, J. M., Crawford, N. M., and Schroeder, J. I. (2000). Cadmium and iron transport by members of a plant metal transporter family in Arabidopsis with homology to NRAMP genes. *Proc. Natl. Acad. Sci. U.S.A.* 97, 4991-4996. doi: 10.1073/pnas.97.9.4991

Untergasser, A., Nijveen, H., Rao, X., Bisseling, T., Geurts, R., and Leunissen, J. A. (2007). Primer3Plus, an enhanced web interface to Primer3. *Nucleic Acids Res.* 35(suppl 2), W71-W74. doi: 10.1093/nar/gkm306

Valentinuzzi, F., Pii, Y., Vigani, G., Lehmann, M., Cesco, S., and Mimmo, T. (2015). Phosphorous and iron deficiencies induce a metabolic reprogramming and affect the exudation traits of the woody plant *Fragaria x ananassa*. *J. Exp. Bot.* 66, 6483-6495. doi: 10.1093/jxb/erv364

Vannozzi, A., Dry, I. B., Fasoli, M., Zenoni, S., and Lucchin, M. (2012). Genome-wide analysis of the grapevine stilbene synthase multigenic family: genomic organization and expression profiles upon biotic and abiotic stresses. *BMC Plant Biol.* 12:130. doi: 10.1186/1471-2229-12-130

Vert, G., Grotz, N., Dédalchéchamp, F., Gaymard, F., Guerinot, M. L., Briat, J. F., et al. (2002). IRT1, an Arabidopsis transporter essential for iron uptake from the soil and for plant growth. *Plant Cell* 14, 1223-1233. doi: 10.1105/tpc.001388

Vigani, G., De Silvestre, D., Agresta, A. M., Donnini, S., Mauri, P., Gehl, C., et al. (2016). Molybdenum and iron mutually impact homeostasis in cucumber (*Cucumis sativus*) plants. *New Phytol.* doi: 10.1111/nph.14214. [Epub ahead of print].

Vitulo, N., Forcato, C., Carpinelli, E. C., Telatin, A., Campagna, D., D'Angelo, M., et al. (2014). A deep survey of alternative splicing in grape reveals changes in the splicing machinery related to tissue, stress condition and genotype. *BMC Plant Biol.* 14:99. doi: 10.1186/1471-2229-14-99

Vorwieger, A., Gryczka, C., Czihal, A., Douchkov, D., Tiedemann, J., Mock, H. P., et al. (2007). Iron assimilation and transcription factor controlled synthesis of riboflavin in plants. *Planta* 226, 147-158. doi: 10.1007/s00425-006-0476-9

Walker, R. R., Blackmore, D. H., Clingeleffer, P. R., and Correll, R. L. (2002). Rootstock effects on salt tolerance of irrigated field-grown grapevines (*Vitis vinifera* L. cv. Sultana) 1. Yield and vigour inter- relationships. *Aust. J. Grape Wine Res.* 8, 3-14. doi: 10.1111/j.1755-0238.2002.tb00206.x

Walker, R. R., Blackmore, D. H., Clingeleffer, P. R., and Correll, R. L. (2004). Rootstock effects on salt tolerance of irrigated field-grown grapevines (*Vitis vinifera* L. cv. Sultana) 2. Ion concentrations in leaves and juice. *Aust. J. Grape Wine Res.* 10, 90-99. doi: 10.1111/j.1755-0238.2004.tb00011.x

Waters, B. M., Blevins, D. G., and Eide, D. J. (2002). Characterization of FRO1, a pea ferric-chelate reductase involved in root
iron acquisition. *Plant Physiol.* 129, 85-94. doi: 10.1104/pp.010829

Weaver, R. J. (1976). *Grape Growing*. New York, NY: A Wiley Interscience Publication.

Xiong, H., Guo, X., Kobayashi, T., Kakei, Y., Nakanishi, H., Nozoye, T., et al. (2014). Expression of peanut Iron Regulated Transporter 1 in tobacco and rice plants confers improved iron nutrition. *Plant Physiol. Biochem.* 80, 83-89. doi: 10.1016/j.plaphy.2014.03.021

Yang, T. J., Lin, W. D., and Schmidt, W. (2010). Transcriptional profiling of the Arabidopsis iron deficiency response reveals conserved transition metal homeostasis networks. *Plant Physiol.* 152, 2130-2141. doi: 10.1104/pp.109.152728

Zamboni, A., Zanin, L., Tomasi, N., Pezzotti, M., Pinton, R., Varanini, Z., et al. (2012). Genome-wide microarray analysis of tomato roots

showed defined responses to iron deficiency. BMC Genomics 13:101. doi: 10.1186/1471-2164-13-101

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2017 Vannozzi, Donnini, Vigani, Corso, Valle, Vito, Bonghi, Zocchi and Lucchin. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.