Impact Factor 5.753 | CiteScore 8.2
More on impact ›

ORIGINAL RESEARCH article

Front. Plant Sci., 17 February 2021 | https://doi.org/10.3389/fpls.2021.616868

Evidences for a Nutritional Role of Iodine in Plants

  • 1Plant Lab, Institute of Life Sciences, Scuola Superiore Sant’Anna, Pisa, Italy
  • 2Proteomics and Mass Spectrometry Laboratory, Institute for the Animal Production System in the Mediterranean Environment (ISPAAM), National Research Council, Napoli, Italy
  • 3Institute of Clinical Physiology, National Research Council, Pisa, Italy
  • 4SQM International N.V., Antwerpen, Belgium

Little is known about the role of iodine in plant physiology. We evaluated the impact of low concentrations of iodine on the phenotype, transcriptome and proteome of Arabidopsis thaliana. Our experiments showed that removal of iodine from the nutrition solution compromises plant growth, and restoring it in micromolar concentrations is beneficial for biomass accumulation and leads to early flowering. In addition, iodine treatments specifically regulate the expression of several genes, mostly involved in the plant defence response, suggesting that iodine may protect against both biotic and abiotic stress. Finally, we demonstrated iodine organification in proteins. Our bioinformatic analysis of proteomic data revealed that iodinated proteins identified in the shoots are mainly associated with the chloroplast and are functionally involved in photosynthetic processes, whereas those in the roots mostly belong and/or are related to the action of various peroxidases. These results suggest the functional involvement of iodine in plant nutrition.

Introduction

Plants need macro- and micro-nutrients for their growth and development. Nutrients are chemical elements that are components of biological molecules and/or influence essential metabolic functions. The elements that to date are considered as plant nutrients are C, H, O, N, P, K (primary nutrients), Ca, Mg, S (secondary nutrients), and Fe, Zn, Cu, Mn, B, Cl, Mo, Co, and Ni (micro-nutrients) (Mengel and Kirkby, 2001).

Halogens are the least represented chemical group of plant micro-nutrients, chloride being the only micro-nutrient currently recognised in plant physiology, due to its regulatory action in proton-transfer reactions in the photosystem II (Brahmachari et al., 2018). Studying the effect of different concentrations and forms of iodine on the growth of several crops of agricultural importance, Borst Pauwels (1961) referred to iodine as a micro-nutrient for plant, and a similar conclusion was derived by Lehr et al. (1958) working on tomato.

A growing number of recent studies reporting the effect of iodine on plant growth have focused on the benefit of increasing iodine content in plants as a biofortificant in human and animal health (Medrano-Macías et al., 2016; Gonzali et al., 2017). Plant tissues generally increase their iodine content following its exogenous administration. However, the presence of iodine as a trace element/contaminant in the soil/nutrient solution/atmosphere cannot be avoided, thus preventing the effects related to the presence/absence of this element from being easily observed (Fuge and Johnson, 1986; Ashworth, 2009). The functional role of iodine as a plant nutrient might therefore have been masked.

Plants can absorb iodine from roots or above ground structures (stomata and cuticular waxes) (Medrano-Macías et al., 2016; Gonzali et al., 2017), translocate it mainly through the xylematic route and volatilise it as methyl iodide (CH3I) through the action of halide ion-methyltransferase (HMT) and halide/thiol methyltransferase (HTMT) enzymes (Medrano-Macías et al., 2016; Gonzali et al., 2017).

Little is known about the chemical forms of iodine inside plant tissues. Inorganic iodine, in particular iodide (I), however, seems to be predominant (Weng et al., 2008). Plants can also incorporate iodine into organic molecules, such as iodosalicylates, iodobenzoates (Smoleń et al., 2020), monoiodotyrosine (MIT), di-iodotyrosine (DIT) and triiodothyronine (T3) (Eales, 1997; Smoleń et al., 2020). Interestingly, MIT and DIT have a key role in the physiology of vertebrates, as they are precursors of the two thyroid hormones (THs) triiodothyronine (T3) and L-thyroxine (T4) as part of the thyroglobulin protein (Zimmermann et al., 2008).

In plants, the presence of a thyroglobulin-like protein has never been reported, and the metabolic role of MIT, DIT and T3 molecules, if any, and their biosynthetic mechanism are still unknown. Nevertheless, protein iodination has been verified in several seaweed species (Hou et al., 2000; Romarís-Hortas et al., 2014), even if it has not yet been demonstrated in plants.

Iodine is likely involved in several physiological and biochemical processes (Medrano-Macías et al., 2016; Gonzali et al., 2017). The presence of low concentrations of iodine is often associated with beneficial effects on plant growth, production and stress resistance, whereas toxic effects are observed when applying iodine at high concentrations, especially in the I form, which is more phytotoxic than iodate (IO3) (Voogt et al., 2010; Medrano-Macías et al., 2016; Gonzali et al., 2017; Incrocci et al., 2019). Thresholds for beneficial or toxic concentrations have been reported for all micro-nutrients (Welch and Shuman, 1995). Interestingly, the concentrations of iodine added to nutrient solutions that have been associated with positive effects for plants (ranging from approximately 102–104 nM) (Medrano-Macías et al., 2016; Gonzali et al., 2017) are comparable with those generally effective for other elements described as plant micro-nutrients (Sonneveld, 2002), suggesting that iodine may play a similar role in plant nutrition.

We explored the role of iodine as a nutrient for plants using various experimental approaches. Our results showed that iodine, when supplied at a well-defined concentration range, positively affected the phenotype of Arabidopsis thaliana plants, and altered the organism’s transcriptome. Most importantly, protein iodination was observed for the first time. These results are strongly suggestive of the role of iodine as a plant nutrient.

Materials and Methods

Plant Material and Cultivation System

Plants of Arabidopsis thaliana, ecotype Columbia 0, Solanum lycopersicum L., cv. Micro-Tom, Lactuca sativa L., var. crispa, Triticum turgidum L., var. durum, and Zea mays L., var. saccharata, were used in the experiments, as summarised in Supplementary Figure S1.

The cultivation protocol commonly applied in all the experiments is described as follows: seeds of the different species were sown on rockwool plugs and vernalised for 3 days. After this period, plants were hydroponically cultivated in a growth chamber (22°C day/18°C night with a 12 h photoperiod, a quantum irradiance of 100 μmol photons m–2 s–1 and a relative humidity close to 35%), in a floating system. A base nutrient solution, renewed once a week, was prepared minimising iodine contamination by dissolving in MilliQ water the following amounts of ultrapure salts: 1.25 mM KNO3, 1.50 mM Ca(NO3)2, 0.75 mM MgSO4, 0.50 mM KH2PO4, 50 μM KCl, 50 μM H3BO3, 10 μM MnSO4, 2.0 μM ZnSO4, 1.5 μM CuSO4, 0.075 μM (NH4)6Mo7O24, and 72 μM Fe-EDTA. At preparation, the pH and the electrical conductivity (EC) values were 6.0 and 0.6 dS m–1, respectively, whereas the iodine concentration in the nutrient solution was below the detection limit of 8 nM, as determined by ICP-MS analysis. The technical peculiarities of each experiment are described in the devoted sections.

Phenotypical Determinations

Two separate experiments were performed. In both experiments, Arabidopsis plants were initially fed with the base nutrient solution. After 15 days of growth, plants homogeneous in size and leaf number were selected, grouped, and fed with different concentration and/or type of halogen-containing salts added to the nutrient solution. Plants were distributed in nine separate hydroponic trays (three different trays/treatment), and a total number of 90 plants were cultivated for each experimental condition (30 plants/tray).

In the first experiment (exp. 1—phenotype), 0.20 or 10 μM KIO3 was added to the nutrient solution. One month later, during the formation/elongation of the main inflorescence, half of the plants (15 plants/tray) was harvested and characterised according to the main morphological traits, such as rosette and inflorescence fresh weight (FW), dry weight (DW), dry matter content, rosette diameter, and inflorescence length. The remaining half was allowed to complete the growing cycle and was characterised in terms of flowering and seed production. Flowering, defined as the presence of the first open flower on the stem, was recorded at intervals of 3 days and expressed on a percentage basis. The percentage of bloomed plants/tray was calculated at each assessment date. Toward the final part of the plants’ life cycle, a periodical harvesting of the produced/matured seeds was carried out until the complete plant desiccation. Seed production was determined in terms of total seed weight/tray (15 plants/tray), number of seeds/silique and number of siliques/plant.

In the second experiment (exp. 2—phenotype), plants were treated by adding either KI, NaI or KBr (0, 10, or 30 μM) to the nutrient solution. Fifteen days after the salt treatment, half of the plants was characterised in terms of plant FW, DW and dry matter content, while the other half was subsequently characterised in terms of flowering (recorded with intervals of 2 days), as described for experiment 1-phenothype.

Total RNA Extraction and Processing

Gene expression analysis was performed on 3-week-old Arabidopsis plants hydroponically grown on the base nutrient solution (control plants) or in the same medium to which 10 μM of KBr, NaI, or KI was added. Plant material was collected 48 h after the beginning of the treatment. Each sample consisted of a pool of rosettes or roots sampled from three different plants, which were immediately frozen in liquid nitrogen and stored at −80°C until further analysis.

Total RNA from rosettes was extracted as described by Perata et al. (1997), avoiding the use of aurintricarboxylic acid. RNA from roots was extracted using the SpectrumTM Plant Total RNA Kit (Sigma-Aldrich). RNA was subsequently processed for microarray and qPCR analysis. The TURBO DNA-free kit (Thermo Fisher Scientific) was used to remove DNA contaminations and the iScript TM cDNA synthesis kit (Bio-Rad Laboratories) was used for RNA reverse-transcription.

Microarray Analysis

RNA from rosettes and roots was processed and hybridised to Affymetrix GeneChip Arabidopsis ATH1 Genome Arrays as described by Loreti et al. (2005). Normalisation was performed using Microarray Suite 5.0 (MAS5.0). Differentially expressed genes (DEGs) were selected based on the two following criteria: log2FC treated/control ≥ 2 and mas5-Detection p ≤ 0.05. In addition, the absolute expression level ≥ 100 mas5-Signal was chosen to select only well-expressed genes. Rosette and root DEGs resulting from KI, NaI, and KBr treatments were processed and visualised in a Venn diagram. Only DEGs commonly regulated by KI- and NaI-treated plants and not responding to KBr treatments were considered specifically linked with the iodine treatment. This group of DEGs was then subjected to gene set enrichment using Gorilla1 and analysed with Mapman2, whereas the co-expression analysis was performed using Genevestigator3.

Gene Expression Analysis (RT-qPCR)

Quantitative PCR (ABI Prism 7300 Sequence Detection System, Applied Biosystems) was performed using 30 ng cDNA and the iQ SYBR Green Supermix (Bio-Rad Laboratories). UBIQUITIN10 (At4g05320) and TIP4 (At2g25810.1) were used as reference genes. Relative expression levels were calculated using GeNorm4. The list of the primers and their sequences are reported in Supplementary Table S1. Four biological replicates were analysed, each consisting of a pool of rosettes or roots sampled from three different plants.

Feeding With Radioactive Iodine

Two separate experiments were performed by feeding radioactive iodine (125I—NaI, PerkinElmer) to hydroponically grown Arabidopsis thaliana (exp. 1—radioactive) or tomato, lettuce, wheat and maize (exp. 2—radioactive) plants. Treatments were performed on 1-month-old plants. The solution of 125I was prepared by dissolving 60 μl of the commercial radioactive 125I product (2.4 mCi/100 μl—9.41 μM) in 250 ml of base nutrient solution. Plants were individually transferred into plastic tubes and treated with the hydroponic solution (with or without Na125I). Sampling was performed after 48 h of incubation by collecting leaf and root material, which was immediately frozen in liquid nitrogen, and stored at −80°C until the analysis. Control, non-treated plants (no 125I added during their growth) were used in both experiments.

Protein Extraction, Electrophoresis, and Gel Autoradiography

Leaf and root samples from 125I-fed and control plants were ground to fine powder in liquid nitrogen. The protein extraction buffer (50 mM TrisHCl, pH 7.0, 1% w/v SDS, P9599 protease inhibitor cocktail, Sigma-Aldrich) was added to the powder. The resulting solution was vortexed vigorously, and then centrifuged (18,407 g, 30 min, 4°C). Radioactive iodine solution (10 μl; prepared as described above) was added to the control samples during the extraction process to check for the occurrence of false positive signals (technical artifacts), possibly due to unspecific binding of iodine with the protein extract.

Protein extracts were dissolved in a 5 × Laemmli buffer, treated at 95°C for 10 min, and a volume of 20 μl (corresponding, approximately, to 65 or 20 μg of proteins, in shoot and root samples, respectively) was loaded to Invitrogen NuPAGE gels (10% Bis-Tris Midi Gels, Thermo Fisher Scientific), together with a protein marker (Precision Plus ProteinTM Dual Color Standards, Bio-Rad). After electrophoresis, the gel was rinsed in MilliQ water, and the proteins were fixed (40:7:53 ethanol/glacial acetic acid/H2O – 30 min) and then stained (EZ Blue Gel Staining reagent, Sigma-Aldrich – 30 min) on an orbital shaker. After rinsing, gels were exposed to a multipurpose phosphor storage screen (Cyclone Storage Phosphor System, PerkinElmer) in order to obtain a digital image of the radioactivity distribution. Radioactive signals were quantified after 72 h of gel exposure using a Cyclone Phosphor Imaging System (PerkinElmer). In order to prevent the occurrence of any radioactive emissions from the control samples, after each image acquisition, gels were re-exposed for 15 days, and the absence of 125I labelled bands was verified in the newly acquired images.

Database Search for Iodinated Peptides in Protein Datasets From Proteomic Data Repositories

Mass spectrometry data were downloaded from the PRIDE (PRoteomics IDEntification database) archive5 (Perez-Riverol et al., 2019). The PRIDE archive was searched to select A. thaliana datasets based on the analysis of specific plant organs, such as cauline, rosette and roots, and/or subcellular districts, such as chloroplasts and mitochondria. Datasets were excluded if enrichment/immunopurification strategies were used during protein purification. Finally, 21 experimental sets of nano-LC-ESI-MS/MS raw data included in 14 PRIDE repositories (March 2020) were obtained and re-analyzed by database searching. Only raw files corresponding to the analysis of control/non-treated plants were downloaded and the experimental protocols and the search parameters for each different dataset were annotated. Supplementary Table S2 lists the experimental sets, with details on the MS instrument, plant organ and/or subcellular compartment, sample preparation, and proteomic strategies adopted.

Raw files were searched separately using Proteome Discoverer 2.4 (Thermo Fisher Scientific, United States) with the Mascot v. 2.6 search engine (Matrix Science Ltd., United Kingdom) against the TAIR10 database6 (71,567 sequences, accessed May 2017) and a database containing common laboratory contaminants on the MaxQuant website7. Workflows were built for each experimental dataset, considering the specific mass tolerance values used for the original search and reported on the PRIDE repository, or in the publication associated with the dataset.

For all the workflows, Cys carbamidomethylation was set as a fixed modification, while iodination at Tyr and His (Δm = +125.8966 Da), oxidation at Met, protein N-terminal acetylation, deamidation at Asp, and pyroglutamate formation at N-terminal Gln were selected as variable modifications. Isotopic labelling was also considered in the modification parameters when performed for protein quantification.

Trypsin was selected as the proteolytic enzyme and peptides were allowed to have a maximum of two missed cleavages. The minimum peptide length was set at six amino acids. The site probability threshold for peptide modification was set at 75. Only high confidence peptide identifications were retained by setting the target false discovery rate (FDR) for PSM at 0.01 and further filtered to keep only peptides (P < 0.05) with a Mascot Ion score > 30. In addition, the results of the identification analysis were processed by putting together the output of iodinated peptides from all the datasets and further applying a filter to keep only those identified with a Mascot Ion score > 50, in at least one dataset, to limit the identification to peptides with the best scoring matches and corresponding to high certainty. The presence of the MS/MS spectrum of the unmodified counterpart was verified for each iodinated peptide to further validate the identification.

Protein Bioinformatics

Proteins containing iodinated peptides were functionally annotated according to MapMan categories by using the Mercator pipeline8. Final outputs were integrated with data from the available literature. Protein interaction networks were obtained with STRING v. 119, which was also used to provide information on known gene ontology categories. Venn diagrams were depicted using a web tool at http://bioinformatics.psb.ugent.be/webtools/Venn. The Protein Abundance Database (PAXdb) was also queried to evaluate quantitative levels of modified A. thaliana proteins at https://pax-db.org.

Statistical Analysis

Data concerning phenotypical determinations and qPCR-based gene expression analysis were analysed by one-way ANOVA coupled with the LSD post hoc test, when they followed a normal distribution and there was homogeneity of variances. When one of these two prerequisites was violated, a Kruskal-Wallis test for non-parametric statistic was performed and the significance letters were graphically assigned using a box-and-whisker plot with a median notch. Significant differences between the means (P < 0.05) are indicated by different letters in the figures and tables.

Results

Effects of Iodine on the Plant Phenotype

The effects of low amounts of KIO3 (0.20 or 10 μM) on hydroponically grown Arabidopsis plants, compared to plants cultivated on a control nutrient solution, were evaluated in terms of plant morphology, biomass, and seed production (exp. 1-phenotype). No phytotoxicity symptoms were observed on plants and the most significant phenotypical effect was a delay of flowering in the control plants, compared to KIO3 (either 0.20 or 10 μM) (Figures 1A,B). Twelve days after the opening of the first flower, plants treated with 0.20 or 10 μM KIO3 were close to complete flowering, as 87 and 96% of the plants had bloomed, respectively, vs. 69% of the control plants (Figure 1B). Control plants took about 18 days to complete blooming.

FIGURE 1
www.frontiersin.org

Figure 1. Impact of iodine on plant growth and development (exp. 1-phenotype). (A) Lateral view of plants after 4 or 6 weeks from the onset of KIO3 treatment. (B) Flowering time curve; the percentage of bloomed plants/tray was calculated every 3 days after the opening of the first flower (day 0). (C) Morphological data on plant FW, DW, dry matter content, rosette diameter, inflorescence length and number of produced siliques/plant, determined 1 month after the onset of KIO3 treatments. Values indicated by different letters significantly differ from each other (according with one-way ANOVA, LSD post hoc test, P ≤ 0.05). In particular, the statistical analysis of flowering (B) was performed by comparing the percentage of bloomed plants of each tray (considered as biological replicates) within each sampling point. When data followed a Normal distribution and there was homogeneity of variances, they were subjected to one-way ANOVA and values indicated by different letters significantly differ from each other (LSD post hoc test, P ≤ 0.05). When one of this two prerequisite was violated, a Kruskal-Wallis test was performed. Error bars (±SE) are shown in graphs.

Plant biomass, evaluated 1 month after the addition of KIO3 to the nutrient solution, was significantly lower in control plants, both in FW and DW (Figure 1C). When compared to the control, the plant FW increased by approximately 7.7 and 13% with addition of 0.20 and 10 μM KIO3 in the nutrient solution, respectively, whereas the DW increased by 13 and 22%, respectively. The effect on plant FW was mostly ascribable to the inflorescence, as no significant differences were evident in terms of the rosette FW values (Supplementary Table S3). The concentration of iodine in the nutrient solution had a marked effect on the inflorescence length, which was approximately 41 and 45% longer compared to the control in 0.20 and 10 μM KIO3 treated plants, respectively (Figure 1C), and a comparable effect was seen on the inflorescence FW and DW (Supplementary Table S3). Additionally, the rosette diameter in the control was smaller, and the application of 0.20 or 10 μM KIO3 increased it by approximately 5 and 9%, respectively (Figure 1C). The plant dry matter content positively correlated with the increased iodate concentrations (Figure 1C).

Seed production was determined in terms of total seed weight, seeds/silique and number of siliques/plant. The number of seeds contained in each silique was not affected by iodate treatments (Supplementary Table S3), whereas the number of siliques produced by each plant was lower in the control, compared to the addition of both 0.20 and 10 μM KIO3 (Figure 1C). This influenced the total seed production, which, 1 month after the addition of KIO3 to the nutrient solution, was much higher in plants treated with iodate (more than 50 and 35%, respectively, in 0.20 and 10 μM KIO3 treated plants in comparison with the control) (Supplementary Table S3).

Adding exogenous iodine in the form of KIO3 countered the delay in flowering of control plants (Figures 1A,B). This was confirmed in experiment 2 (exp. 2-phenotype), when iodine was added in the form of KI or NaI (Figures 2A,B). The possible effects of potassium or bromide, as an alternative halogen, were evaluated and then ruled out, as a similar behaviour was observed in plants treated with KI or NaI, but not with KBr (Figures 2A,B).

FIGURE 2
www.frontiersin.org

Figure 2. Impact of iodine on plant growth and development (exp. 2-phenotype). (A) Flowering curve; the percentage of bloomed plants/tray was calculated every 2 days after the opening of the first flower (day 0). (B) Representative control, and KI-, NaI- or KBr-treated plants (10 and 30 μM) after 15 days from the onset of the treatments. Pictures were taken after 4 days from the opening of the first flower on the main stem. (C) Morphological data on plant FW, DW, dry matter content determined 15 days after the onset of the treatments. Values indicated by different letters significantly differ from each other (according with one-way ANOVA, LSD post hoc test, P ≤ 0.05). In particular, the statistical analysis of flowering (A) was performed by comparing the percentage of bloomed plants of each tray (considered as biological replicates) within each sampling point. When data followed a Normal distribution and there was homogeneity of variances, they were subjected to one-way ANOVA and values indicated by different letters significantly differ from each other (LSD post hoc test, P ≤ 0.05). When one of this two prerequisite was violated, a Kruskal-Wallis test was performed. Error bars (±SE) are shown in graphs.

The application of 10 μM KI and NaI promoted flowering, without negatively impacting the plant biomass production (Figure 2C), whereas 30 μM KI or NaI reduced plant growth (Figure 2C), although the promoting effect of iodine on flowering was still present (Figures 2A,B). Four days after the opening of the first flower (day 0), more than 50% of KI- and NaI-treated plants had bloomed vs. 14% of the control plants and 10% and 14% of the 10 and 30 μM KBr-treated plants, respectively. Moreover, the floral transition was almost complete in 10 μM KI- and NaI-fed plants in the subsequent 6 days (10 days after day 0). Two and four more days were required for 30 μM KI- and NaI-fed plants, respectively (12 and 14 days after day 0), whereas the control and KBr-treated plants completed blooming in the subsequent 18 days (22 days after day 0) (Figure 2A).

Effects of Iodine on Gene Expression

The response of plants to iodine was analyzed at the transcriptomic level. To identify genes whose expression was specifically altered by iodine, Arabidopsis plants were treated by adding 10 μM of NaI, KI, or KBr to the nutrient solution, compared to the untreated control plants. The resulting RNAs were analyzed by hybridisation on ATH1 microarrays. To rule out the possible generic effects of halogens, we searched the microarray dataset for genes that responded to KI and NaI, but not to KBr. In addition, a comparison between KBr- and KI-treated plants enabled us to rule out the possible transcriptional regulation of genes exerted by potassium, as K+ ion was common to both salts.

Data visualisation with a Venn diagram showed that several genes were specifically regulated by iodine, as up- or down-regulated genes in both NaI- and KI- but not in KBr-treated plants were 33 (51% of DEGs) and 15 (33% of DEGs) in the shoot (Figure 3A), and 398 (95% of DEGs) and 133 (79% of DEGs) in the root (Figure 3B), respectively. The similarity and specificity in the expression pattern of KI- and NaI-treated plants were confirmed by the heatmaps generated from the analysis of the shoot (Figure 3C) and root (Figure 3D) expression data.

FIGURE 3
www.frontiersin.org

Figure 3. Transcriptional regulation of gene expression induced by iodine. Venn diagram showing the number of genes differentially regulated in shoot (A) or root (B) tissues of KBr-, NaI-, and KI-treated plants (10 μM—48 h), when compared with the control. Heatmap showing the pattern of expression of the genes analysed in the shoot (C) or root (D) tissues in response to NaI, KI or KBr treatments, when compared with the control. qPCR validation of selected genes up- or down-regulated by iodine treatments (commonly regulated by NaI and KI, but not KBr) in shoot (E) or root (F) tissues. qPCR data are mean ± SE of four biological replicates, each composed of a pool of three different rosettes. Values indicated by different letters significantly differ from each other (according with one-way ANOVA, LSD post hoc test, P ≤ 0.05).

To validate the microarray analysis, a subset of three I-induced and three I-repressed genes were analysed by qPCR, corroborating the specific regulation of iodine on their expression in both shoot (Figure 3E) and root (Figure 3F) samples.

The complete list of the KI and NaI commonly and not responding to KBr up- and down-regulated genes is reported in Supplementary Tables S4, S5 (shoot tissue), and Supplementary Tables S6, S7 (root tissue), respectively.

The polypeptides codified by the iodine-regulated genes did not show a preferential site of action in the cell, as their predicted localisations include cytoplasm, chloroplast, cell wall, nucleus, mitochondrion, vacuole and apoplast (Supplementary Tables S4S7).

The gene ontology (GO) analysis identified several functional categories regulated by iodine in the roots (Figure 4, Supplementary Figure S2 and Supplementary Tables S8, S9), whereas no statistically significant GO terms were identified by analysing the DEG data on the shoots.

FIGURE 4
www.frontiersin.org

Figure 4. Overview of the main biological processes affected by iodine based on the GO terms enrichment analysis carried out in root tissues. Only genes regulated in NaI- and KI-treated plants, and not in KBr-treated plants, when compared with the control, were analysed. The figure was extracted from GOrilla (http://cbl-gorilla.cs.technion.ac.il) and reproduced. In this analysis, DEGs with log2FC ≥ 2.5 or log2FC ≤ –2.5 were used.

The most representative biological processes affected by iodine in the roots were related to the response to stimulus (GO:0050896), and the downstream categories associated with response to abiotic (GO:0009628) and biotic stimulus (GO:0009607) (Figure 4 and Supplementary Table S8). The main molecular functions regulated by iodine in the roots were related to antioxidant (GO:0016209) and oxidoreductase activity (GO:0016491) and related child terms, in particular peroxidase activity (GO:0004601) and oxidoreductase activity, acting on peroxide as acceptor (GO:0016684) (Supplementary Figure S2 and Supplementary Table S9).

DEGs analysis performed with MapMan highlighted an over-representation of several genes in root samples that were related to calcium regulation and protein modification/degradation (Supplementary Figure S3A), together with genes encoding for the large enzyme families including peroxidases, oxidases, glutathione S-transferases, and cytochrome P450 (Supplementary Figure S3B).

The relatively low number of genes regulated by iodine in the shoots prevented a gene ontology analysis from being performed. However, in terms of the most well characterised genes specifically regulated by iodine treatments in the shoot, the main pathways affected were directly or indirectly involved in biotic (approximately 48 or 40% of up- or down-regulates genes, respectively) or abiotic (approximately 45 or 33% of up- or down-regulates genes, respectively) stress response pathways (Supplementary Tables S4, S5). Several genes playing a role in the transition to flowering (At4g19191 and At1g75750) and embryo and pollen development (i.e., At1g21310, At3g54150) are also worth mentioning.

The involvement of iodine in the defence response, highlighted by the previous analyses performed on root samples, was also suggested by querying all publicly available microarray datasets (see footnote) using the list of iodine-responsive genes of both shoot (Supplementary Figure S4) and root (Supplementary Figure S5) tissues. The majority of the up- or down-regulated genes were commonly modulated by the presence of fungal infection, salicylic acid (SA) or synthetic analogues of SA, such as benzothiadiazole (Kouzai et al., 2018).

Protein Iodination in Plants

Iodine can be found in plant tissues not only in a mineral form but also in organic compounds (Wang et al., 2014; Smoleń et al., 2020). To verify the possible in vivo incorporation of iodine into proteins, we carried out two different experiments by feeding hydroponically grown plants with 125I and carrying out the autoradiography of the SDS-PAGE of the relative protein extracts to detect possible radio-labelled proteins. The experiments were performed first with Arabidopsis plants, and then with other species, namely maize, tomato, wheat and lettuce.

The experiment carried out with Arabidopsis plants revealed the presence of at least six radio-labelled bands at different molecular mass values in the protein extracts from shoot tissues (Figure 5A; exp. 1-radioactive) and eleven radio-labelled bands from root tissues (Figure 5B; exp. 1-radioactive), indicating the presence of proteins likely containing iodo-amino acids. Iodinated proteins were preferentially present in root tissues, as the abundance and intensity of 125I-labelled bands were higher in the root than in the shoot extracts. No radioactive signals were observed in the shoot and root control samples (samples added with 125I solution during protein extraction).

FIGURE 5
www.frontiersin.org

Figure 5. Autoradiographies of the SDS-PAGE gels. Comparison between the position and relative intensities of 125I radiolabelled bands of representative shoot (A) and root (B) protein extracts from 125I treated Arabidopsis (exp. 1-radioactive), and maize, tomato, wheat and lettuce plants (exp. 2-radioactive). Sampling was performed after 48 h of 125I incubation. In both the experiments, autoradiographies were acquired after 72 h of gel exposition to the multipurpose phosphor storage screen. Representative pictures of total stained protein extracts (SDS-PAGE) and of the autoradiographies of control samples after 15 days of exposition are also shown. Controls consisted in protein extracts obtained from plants untreated with 125I during their growth, to which the radioactive solution containing 125I was added during the extraction process.

Several 125I -labelled bands were also observed in the leaf extracts of tomato, wheat and lettuce samples, whereas no 125I-containing bands were visible in the leaf protein extracts of maize (Figure 5A; exp. 2-radioactive). A clear radioactive signal was detected in several root proteins extracted from all the species analysed, including maize (Figure 5B; exp. 2-radioactive). Also in this case, the intensity of the radiolabelled bands was higher in root than in shoot extracts. A good degree of conservation of the molecular mass values of the putatively iodinated proteins was observed among the five plant species analysed (Figure 5).

Identification of Iodinated Proteins in Arabidopsis thaliana

The identification of the radiolabelled proteins described above was hampered by the presence of a radioactive isotope, which meant that our samples did not meet the safety rules for proteomic facilities. To maximise the probability of success in identifying targets of protein iodination, we then focused on the nano-LC-ESI-MS/MS raw data already acquired within the framework of experimental studies on different organs/subcellular districts of Arabidopsis thaliana, and released in the public repository PRoteomics IDEntification Database (PRIDE) Archive (Perez-Riverol et al., 2019).

The datasets considered for our analysis refer to many different experimental conditions in terms of plant growth, treatment, and cultivation regimen, as well as sample processing and fractionation performed before proteomic analysis. No experiments were explicitly related to iodination studies; the presence of iodine occurred accidentally, as a consequence of its natural presence in the cultivation environment (i.e., soil, air, irrigation water), or because it was conventionally present in the MS growing medium (Murashige and Skoog, 1962), which is widely used in studies based on in vitro plant tissue culture.

In proteins, iodination affects various amino acids, depending on the reaction conditions (Ramachandran, 1956), but generally following the reactivity order Tyr >> His ≥ Trp > Cys.

Mono-iodination at Tyr and His residues were thus considered in the searching parameters as variable modifications. The output of the database search, in terms of proteins iodinated at Tyr or His residues has been reported in Supplementary Table S10.

The iodinated peptides were identified in 16 out of the 21 datasets analysed. A total of 106 iodinated peptides, corresponding, respectively, to 42 and 40 protein accessions in the TAIR10 database of A. thaliana leaves (chloroplast, cauline, and rosette) (Table 1) and roots (Table 2), were identified. Iodinated sequences differently modified for deamidation, and/or protein N-terminal acetylation, and/or Met oxidation are reported as a unique iodinated peptide inventory.

TABLE 1
www.frontiersin.org

Table 1. Iodinated peptides identified in A. thaliana leaves (chloroplast, cauline, and rosette) by database searching of mass spectrometric data retrieved from PRIDE repository.

TABLE 2
www.frontiersin.org

Table 2. Iodinated peptides identified in A. thaliana roots by database searching of mass spectrometric data retrieved from PRIDE repository.

Most of the modified peptides were found to be iodinated at Tyr residues, while His iodination was identified in only five peptides. Representative MS/MS spectra of Tyr-iodinated peptides are reported in Figure 6A.

FIGURE 6
www.frontiersin.org

Figure 6. Iodination in A. thaliana proteins identified by database searching of nano-LC-ESI-MS/MS raw data from a public repository (PRIDE). (A) Unambiguous assignment of iodination sites by MS/MS analysis in two peptides from chloroplast (light harvesting complex of photosystem II 5, AT4G10340.1), upper panel, and roots (peroxidase superfamily protein, AT4G30170.1), lower panel. The peptides are identified by both y, and b ions. Red labels in the spectra evidence the mass shift corresponding to the iodinated tyrosine (i-Y). (B) Venn Diagram showing the iodinated peptide sequences identified in the datasets of chloroplasts (Chlor), cauline (Cau), rosette (Ros), and roots (Root).

To evaluate the entire output of iodinated peptides identified, iodinated sequences for chloroplasts, caulines, rosettes, and roots were processed and visualised in a Venn diagram (Figure 6B). This showed the presence of the common iodinated peptides for the chloroplast cauline and rosette subsets. The root subset was clearly distinct from the other three subsets that were all from the green parts of the plant.

Iodinated Proteins in A. thaliana Leaves

Iodinated peptides identified in 11 datasets of cauline, rosette, and leaf-isolated chloroplasts were assigned to 42 proteins (Supplementary Table S10). Most of the modified species were in the dataset of chloroplastic proteins (Figure 6B). STRING interaction analysis of the modified proteins revealed a single network of 40 iodinated proteins (PPI enrichment p< 1.0e–16) (Supplementary Figure S6A and Supplementary Table S11). A total of 31 of the 40 proteins in this network were involved in photosynthesis (GO:0015979) (Supplementary Table S12), as also attested by their functional analysis, according to MapMan categories (Supplementary Figure S6B and Supplementary Table S13). Moreover, the main molecular functions in which iodinated proteins are involved were related to chlorophyll binding (GO:0016168), protein domain specific binding (GO:0019904), tetrapyrrole binding (GO:0046906), and electron transfer activity (GO:0009055) (Supplementary Table S12). In particular, ten proteins were identified in two or more datasets, thus resulting in the most representative targets of the iodination reaction. Some of these proteins were constituents of the photosystem II (PSII), i.e., proteins of the reaction centre (PsbA, PsbB), oxygen evolving centre (PsbO, PsbP) and light harvesting complex II (LHCB1B1), or part of the photosystem I (PSI) (PsaF, PETE2). Three proteins showed ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCo) activity and were involved in the Calvin Cycle (RBCL, ORF110A, RBCS1A). These proteins were largely abundant in leaf tissues, especially the PSII component PsbB and RuBisCo large and small chains (RBCL, RBCS1A), according to the data reported in PAXdb.

Iodinated Proteins in A. thaliana Roots

Iodinated peptides identified in 5 datasets of roots were assigned to 40 proteins (Supplementary Table S10). The STRING interaction analysis recognised three networks containing 24 of the 40 iodinated proteins (Supplementary Figure S7A and Supplementary Table S11). The GO analysis for these proteins showed a significant over-representation of biological processes related to the response to stress (GO:0006950), response to oxidative stress (GO:0006979), response to toxic substances (GO:0009636), and response to stimulus (GO:0050896). For the molecular functions, the most enriched categories were copper ion binding (GO:0005507) and peroxidase activity (GO:0004601) (Supplementary Table S14).

The functional analysis of the iodinated proteins in roots, according to MapMan categories, highlighted a broad range of biological roles (Supplementary Figure S7B and Supplementary Table S13). In particular, among the iodinated proteins identified, 12 were found in two or more datasets. Five proteins belonged to the classical plant (class III) peroxidase subfamily (At4g30170, At1g05240, At2g37130, At3g01190, At5g17820). The alignment of the protein sequences of the peroxidases mentioned above showed that iodinated residues in all peptides preferentially corresponded to conserved Tyr residues, while only two iodinated tyrosines were unrelated (Supplementary Figure S8).

The other iodinated proteins included: (i) copper amine oxidase (CUAOy2), a cell-wall oxidase showing primary amine oxidase activity; (ii) beta-galactosidase 5 (BGAL5), a glycoside hydrolase involved in the modification of cell wall polysaccharides; (iii) glycosyl hydrolases family 32 protein (ATBFRUCT1), acting as a cell wall invertase; (iv) Pole Ole1 allergen/extension domain (IPR006041)-containing proline-rich protein-like 1 (PRPL1-MOP10) and root hair specific 13 protein (RHS13), which are cell-wall components; (v) D-mannose binding lectin protein (MBL1); and (vi) glyceraldehyde-3-phosphate dehydrogenase C sub 1 (GAPC1), a key enzyme in glycolysis. According to PAXdb, some of the above-reported proteins are abundant in the roots.

Discussion

Iodine Influences Plant Growth and Development and Can Modulate the Plant Transcriptome

Establishing whether iodine is important for a plant’s life is complex, as it is always present in variable amounts in the soil, water, and atmosphere. Plants can take up iodine from the soil solution through the root system, but they also assimilate it from the air or absorb it through the leaves if dissolved in salt solutions or in rain. All these processes occur naturally (Fuge and Johnson, 1986; Ashworth, 2009), thus a plant cannot be grown in the complete absence of iodine.

To identify whether iodine can act as a micro-nutrient we supplied it to plants at very low concentrations (in the micromolar range). These concentrations are typical of many mineral elements that are beneficial or essential when taken up in low doses, and phytotoxic when in excess (Welch and Shuman, 1995). We observed a difference in plant growth between the control and iodine-treated plants. Where these could be perceived as positive effects of addition of a beneficial compound, these can also be interpreted as a negative effect of removal of iodine from the plant’s nutrition. An increase in biomass and seed production, together with a very evident hastening of flowering, was observed by feeding plants with KIO3 (Figure 1 and Supplementary Table S1) or KI (Figure 2) at 0.2 and/or 10 μM. On addition of iodine, flowering was always early and appeared to be specific for iodine, since it was present in the KIO3, KI and NaI treatments but completely absent in KBr (Figures 1, 2). However, the positive effects of iodine on growth were lost at 30 μM. This suggests that a concentration of 30 μM applied as I may be above the toxicity threshold.

In the range of 1–10 μM iodine increases biomass in vegetables, e.g., spinach (Zhu et al., 2003), lettuce (Blasco et al., 2013), tomato (Lehr et al., 1958; Borst Pauwels, 1961), and strawberry (Li et al., 2016), or staple crops, e.g., barley (Borst Pauwels, 1961) and wheat (Cakmak et al., 2017). In tomato, Lehr et al. (1958) demonstrated that treatments with iodine accelerated plant growth, causing early flowering associated with an increase in yield. Similarly, Umaly and Poel (1970) found that the addition of 4 μM KI to the nutrient solution stimulated tomato plants to produce flowers 2–3 days earlier than the control, whereas the use of iodine at a higher concentration (80 μM KI) delayed flower formation and reduced the number of inflorescences. It must be noted that in most biofortification studies with iodine, its native occurrence in nutrient solution or soil of the control plants is not always reported; where iodine concentrations in leaf or root tissue in control plants are reported, these always indicate that iodine was available for uptake and accumulation, regardless of the exogenous administrations (e.g., Borst Pauwels, 1961; Cakmak et al., 2017).

Flowering is a complex physiological process affected by a multitude of internal and external factors, and its hastening may represent an evolutionary adaptive mechanism to guarantee species survival, by optimising the seed set in the case of biotic or abiotic stresses (Ionescu et al., 2017). In Arabidopsis, heat and drought stress are correlated with early flowering, which in turn is generally associated with a reduction in plant growth (Balasubramanian et al., 2006; Schmalenbach et al., 2014). In our study, 10 μM iodate or iodide promoted flowering without negatively impacting on biomass production (Figures 1, 2), which was increased by KIO3, thus suggesting the specific flowering-promoting role of iodine in the process.

Our transcriptomic analysis of plants treated with 10 μM KI, NaI, or KBr for 48 h showed that several genes were specifically regulated by iodine, as the large part of DEGs similarly responded to KI and NaI, but not to KBr (Figure 3). This was more evident in root than in shoot tissues, probably because iodine was added to the nutrient solution and was used first by roots before the green parts. Interestingly, transcripts specifically regulated by iodine in the roots were mostly involved in the plant response to biotic and/or abiotic stresses (Figure 4 and Supplementary Tables S6S8) and the selective regulation of iodine on these groups of genes was also observed in the shoots (Supplementary Tables S4, S5 and Supplementary Figure S4).

Although no previous data are available on the response of Arabidopsis to iodine at the transcriptomic level, the induction of HALIDE ION METHYLTRANSFERASE, SALICYLIC ACID CARBOXYL METHYLTRANSFERASE, and SALICYLIC ACID 3-HYDROXYLASE genes by aromatic iodine compounds, indicating a possible involvement of iodine in the SA metabolism, has already been described in tomato plants (Halka et al., 2019). SA is a signalling molecule involved in local defence reactions at infection sites and the induction of systemic resistance (Vlot et al., 2009).

Iodine likely has an indirect effect on plant resistance given that iodine treatments induce the biosynthesis of several enzymatic or non-enzymatic compounds involved in the plant response to environmental stresses (Leyva et al., 2011; Blasco et al., 2013; Gupta et al., 2015). The antioxidant response induced by increasing KI and KIO3 levels was found to be strongly associated with the synthesis of phenolic compounds (Incrocci et al., 2019; Kiferle et al., 2019), in agreement with our transcriptomic data (Figure 4 and Supplementary Table S8). Iodine treatments have also been associated with the modulation of the essential oil composition in basil plants (Kiferle et al., 2019), which plays a key role in defensive and attraction mechanisms in response to the environment (Cseke et al., 2007; Bakkali et al., 2008).

In our study, the relationship between iodine and plant resistance to stress was also suggested by the activation of GO terms associated with hypoxia (Figure 4 and Supplementary Table S8) and the over-representation of data points related to calcium (Ca) regulation (Supplementary Figure S3A). Ca plays a central role in the plant perception of stress by activating a general defence mechanism, which relies on a Ca spiking mechanism and thus on the battery of Ca-dependent proteins that sense Ca and transduce the signal to downstream targets (Dodd et al., 2010; Pucciariello et al., 2012). Our microarray results did not provide evidence that iodine regulated the expression of anaerobic core genes, such as PYRUVATE DECARBOXYLASE 1 or ALCOHOL DEHYDROGENASE (Mustroph et al., 2010). This suggests that the regulation of hypoxic genes by iodine was not associated with limiting O2 levels but with an unspecific plant response to an environmental/biotic stress, as it is also highlighted by the activation of enzyme families associated with antioxidant response and xenobiotic detoxification, such as peroxidases, cytochrome P450 or glutathione S-transferases (Supplementary Figure S3B; Mittler et al., 2004; Jouili et al., 2011; Pandian et al., 2020).

At very low doses, iodine thus activates a general defence response which takes place before any biotic or abiotic danger and thus may prepare the plant for a possible future attack or environmentally unfavourable conditions.

The impact of iodine on the transcriptome related to defence was impressive. For evaluation of iodine as a plant nutrient the possible influence on growth and development is also important, and this appeared to be sustained by more scattered elements, requiring further studies. For instance, we found iodine-driven modulation in the shoots of the expression of specific genes that are known to regulate flowering (Xing et al., 2013; Trapalis et al., 2017). This is line with our phenotypical data that indicate that iodine promotes early blooming (Figures 1B, 2A).

Iodine Is Incorporated Into Plant Proteins

Using radiolabelled iodine, we observed iodine incorporation into leaf and root proteins in various plant species, with a good level of conservation of the molecular mass values (Figure 5), suggesting that iodination plays a functional role in specific proteins. The absence of radioactive bands in maize leaf extracts may be due to a specific characteristic of the species: maize was the only C4 plant analysed in our study.

Under alkaline conditions, iodine is known to react in vitro with free amino acids, such as Tyr and His, possibly leading to the formation of several I-labelled proteins (Scott, 1954). In our experiments, however, no radioactive signals were present when 125I was added to shoot and root control samples after protein extraction, indicating that protein iodination occurred in vivo.

To the best of our knowledge, the presence of naturally occurring iodinated proteins in higher plants has never been described before, although it is well known in seaweed. For instance, the fraction of iodine bound to proteins in Sargassum kjellanianum accounts for 65.5% of the total element content of this organism (Hou et al., 2000). In addition, 1D and 2D gel-based electrophoresis combined with laser ablation ICP-MS highlighted several iodinated proteins in Nori seaweed, although no further analyses were conducted to identify them (Romarís-Hortas et al., 2014).

In this study, we identified several iodinated peptides from several high-quality proteomic datasets on plants grown without intentional enrichment of the growing media with iodine, but accumulating the iodine naturally present in soil and water or in MS media. Interestingly, these iodinated peptides belong to proteins, which appeared to be involved in well-defined biological contexts within the plant (Tables 1, 2 and Supplementary Table S10).

In roots, some of the iodinated protein networks that have been identified belong to various class III peroxidases (EC 1.11.17). These are a large molecular family of isozymes present in higher plants, which catalyse redox processes between H2O2 and reductants and are involved in growth, cell wall differentiation and the response to various biotic/abiotic stresses (Moerschbacher, 1992; Hiraga et al., 2001). They are abundant in root tissues, in which high concentrations of H2O2 occur during root development (Dunand et al., 2007). In the presence of low concentrations of iodide, various peroxidases catalytically degrade H2O2 at neutral pH values, generating hypoiodous acid (HIO) (Davies et al., 2008; Vlasova, 2018), a strong iodinating agent that can modify proteins in the environment surrounding the site of its generation, also inducing modification at Tyr. Our findings thus suggest that the class III peroxidases involved in the neutralisation of H2O2 present in root tissues, in the presence of iodide and resulting HIO, become themselves the targets of iodination, which in our study occurred at different Tyr residues in their structure. In fact we also found that Tyr iodination occurred in other proteins that directly/indirectly interact with or are functionally linked to class III peroxidases, and/or are present at high levels in root tissues (Supplementary Figure S7A and Supplementary Table S11).

Regarding leaves, we identified a number of iodinated peptides from proteins in proteomic datasets from Arabidopsis chloroplast, cauline and rosette extracts (Figure 6B). Most of these proteins are well-known constitutive subunits of molecular complexes (PSII, PSI, Cytb6f and ATPase) present in the plant photosynthetic machinery (Dekker and Boekema, 2005; Nevo et al., 2012). In in vitro labelling experiments (Machold and Aurich, 1981) demonstrated that some of them can be iodinated.

The above-mentioned macromolecular assemblies are involved in the generation and transfer of reactive electrons, from their early formation up to the coupling reactions, where their chemical potential allows the generation of plant ATP, NADP and carbohydrates (Foyer, 2018). In this context, high light intensity is one of the major stress factors in green plant tissues, which leads to the production of highly reactive oxygen species (ROS) (1O2, O2⋅–, H2O2), and hydroxyl radicals (Schmitt et al., 2014; Foyer, 2018). These reactions also occur during normal photosynthetic conditions and, whenever not controlled by dedicated endogenous antioxidants, can impact on photosynthetic efficiency.

The principal target of light stress is the chloroplast, which is the preferential site of iodine accumulation in the leaf (Weng et al., 2013), and PSI and PSII are the main sites of O2⋅– and 1O2 production, respectively. High light illumination and the corresponding generation of ROS cause photoinhibition of PSII as well as the modification of other photosynthetic complexes, and induce an accelerated turnover of components of these molecular machineries (Li et al., 2018). The latter phenomenon is generally accomplished through the rapid degradation of photo-damaged proteins and concomitant substitution of them with newly synthetised functional copies (Aro et al., 2005; Yamamoto et al., 2014). This oxidative modification is also known to affect redox-regulated enzymes involved in the Calvin cycle. Proteomic studies have already characterised the nature and the sites of various oxidative and nitrosative modifications at Tyr, Trp and His residues in components of PSII, PSI, Cytb6f and ATP synthase complexes in thylakoid membranes from plants exposed to intense illumination (Galetskiy et al., 2011). These modifications were induced by ROS and/or other reactive nitrogen species, e.g., peroxinitrite, which are formed after the reaction of ROS with nitric oxide and other plant nitrogenous species (Bachi et al., 2013; Lu and Yao, 2018).

Given that such oxidised and nitrated proteins coincide with those found iodinated in our study, ROS likely also react with iodo-containing ions present in the chloroplasts to generate iodinating species that affect Tyr and His residues. We found that iodination processes also affected other proteins functionally related to subunits of PSII, PSI, Cytb6f and ATP synthase complexes, such as those involved in the Calvin Cycle, which are present at high concentrations in the same subcellular district, and have already been reported to directly/indirectly interact with the above-mentioned photosynthetic assemblies (Supplementary Figure S6A and Supplementary Table S11).

The addition of 1O2 to I forms peroxyiodide, which decomposes into highly reactive iodine and iodo-containing radicals during dye-mediated photodynamic bacterial inactivation in the presence of potassium iodide (Wen et al., 2017). This process is fundamental for antimicrobial photodynamic therapy (Hamblin and Abrahamse, 2018). Thus, in the presence of iodide, the above-mentioned photoactivated reactions or their possible process variants may have contributed to generate highly reactive iodo-containing molecules and/or radicals, thereby leading to the formation of the iodinated proteins observed in this study. The possible functional meaning of these is certainly worth studying.

Conclusion

We demonstrated that very low amounts of iodine (between 0.20 and 10 μM, i.e., in the range of the concentrations required by plants for several micro-nutrients) improved plant growth and development thus promoting both biomass production and early flowering, and that this effect could not be achieved by another halogen most resembling iodine (Br). Secondly, we found that iodine was able to modulate gene expression in a specific way, activating multiple pathways, mostly involved in defence responses. Finally, we demonstrated that iodine can be a structural component of several different proteins, and conserved iodinated proteins are synthesised in both the roots and shoots of phylogenetically distant species.

These three lines of evidence highlight that iodine has a nutritional role in plants. This means that the influence of iodine on plants is not merely the result of an indirect priming effect by a potentially phytotoxic compound. Considering that plant nutrients are chemical elements that are components of biological molecules and/or influence essential metabolic functions, iodine matches at least the first part of this definition. Further studies on the importance of organification of iodine in proteins on their catalytic and/or regulatory function will help to complete the picture on the functional role of iodine as a plant nutrient.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/geo/, GSE157643.

Author Contributions

PP, CK, MM, KH, and HTH: conceived the project. CK, MM, and SB: experiments on plant phenotype and transcriptomic. PS, CK, and MM: experiments on radioactive feeding. AMS: bioinformatics analysis of mass spectrometry-based proteomics data sets. CK and AMS: data analysis and figures preparation. CK, AMS, AS, and SG: writing—original draft preparation. PP, AS, SG, PS, KH, and HTH: general discussion and revision of the article. All authors read and contributed to previous versions and approved the final version.

Funding

This work was supported by the Scuola Superiore Sant’Anna and SQM International N.V. The funder was not involved in the study design, collection, analysis, interpretation of data, the writing of this article or the decision to submit it for publication.

Conflict of Interest

KH and HH were employees of SQM International N.V., a company active in the sector of fertilisers.

The remaining 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.

Acknowledgments

This manuscript has been released as a pre-print at BiorXiv (Kiferle et al., 2020).

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.616868/full#supplementary-material

Footnotes

  1. ^ http://cbl-gorilla.cs.technion.ac.il
  2. ^ http://mapman.gabipd.org/mapman
  3. ^ https://genevestigator.com
  4. ^ http://medgen.ugent.be/j̃vdesomp/genorm
  5. ^ https://www.ebi.ac.uk/pride/archive
  6. ^ www.arabidopsis.org
  7. ^ https://www.maxquant.org/maxquant/
  8. ^ http://mapman.gabipd.org/web/guest/app/Mercator
  9. ^ http://string-db.org

References

Aro, E. M., Suorsa, M., Rokka, A., Allahverdiyeva, Y., Paakkarinen, V., Saleem, A., et al. (2005). Dynamics of photosystem II: A proteomic approach to thylakoid protein complexes. J. Exp. Bot. 56, 347–356. doi: 10.1093/jxb/eri&041

CrossRef Full Text | Google Scholar

Ashworth, D. J. (2009). “Transfers of Iodine in the soil-plant-air System: Solid-Liquid Partitioning, Migration, Plant Uptake and Volatilization,” in Comprehensive Handbook of Iodine, eds V. R. Preedy, G. N. Burrow, and R. Ross Watson (Cambridge: Academic Press), 107–118.

Google Scholar

Bachi, A., Dalle-Donne, I., and Scaloni, A. (2013). Redox proteomics: chemical principles, methodological approaches and biological/biomedical promises. Chem. Rev. 113, 596–698. doi: 10.1021/cr300073p

PubMed Abstract | CrossRef Full Text | Google Scholar

Bakkali, F., Averbeck, S., Averbeck, D., and Idaomar, M. (2008). Biological effects of essential oils – A review. Food Chem. Toxicol. 46, 446–475. doi: 10.1016/j.fct.2007.09.106

PubMed Abstract | CrossRef Full Text | Google Scholar

Balasubramanian, S., Sureshkumar, S., Lempe, J., and Weigel, D. (2006). Potent induction of Arabidopsis thaliana flowering by elevated growth temperature. PLoS Genet. 2:0980–0989. doi: 10.1371/journal.pgen.0020106

PubMed Abstract | CrossRef Full Text | Google Scholar

Blasco, B., Leyva, R., Romero, L., and Ruiz, J. M. (2013). Iodine effects on phenolic metabolism in lettuce plants under salt stress. J. Agric. Food. Chem. 61, 2591–2596. doi: 10.1021/jf303917n

PubMed Abstract | CrossRef Full Text | Google Scholar

Borst Pauwels, G. W. F. H. (1961). Iodine as a micronutrient for plants. Plant Soil 14, 377–392. doi: 10.1007/BF01666295

CrossRef Full Text | Google Scholar

Brahmachari, U., Guo, Z., Konecny, S. E., Obi, E. N. C., and Barry, B. A. (2018). Engineering proton transfer in photosynthetic oxygen evolution: Chloride, Nitrate, and Trehalose reorganize a hydrogen-bonding network. J. Phys. Chem. B 122, 6702–6711. doi: 10.1021/acs.jpcb.8b02856

PubMed Abstract | CrossRef Full Text | Google Scholar

Cakmak, I., Prom-u-thai, C., Guilherme, L. R. G., Rashid, A., Hora, K. H., Yazici, A., et al. (2017). Iodine biofortification of wheat, rice and maize through fertilizer strategy. Plant Soil 418, 319–335. doi: 10.1007/s11104-017-3295-9

CrossRef Full Text | Google Scholar

Cseke, L. J., Kaufman, P. B., and Kirakosyan, A. (2007). The biology of essential oils in the pollination of flowers. Nat. Prod. Commun. 2, 1317–1336. doi: 10.1177/1934578X0700201225

CrossRef Full Text | Google Scholar

Davies, M. J., Hawkins, C. L., Pattison, D. I., and Rees, M. D. (2008). Mammalian heme peroxidases: From molecular mechanisms to health implications. Antioxid. Redox Signal. 10, 1199–1234. doi: 10.1089/ars.2007.1927

PubMed Abstract | CrossRef Full Text | Google Scholar

Dekker, J. P., and Boekema, E. J. (2005). Supramolecular organization of thylakoid membrane proteins in green plants. Biochim. Biophys. Acta Bioen. 1706, 12–39. doi: 10.1016/j.bbabio.2004.09.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Dodd, A. N., Kudla, J., and Sanders, D. (2010). The language of calcium signaling. Annu. Rev. Plant Biol. 61, 593–620. doi: 10.1146/annurev-arplant-070109-104628

PubMed Abstract | CrossRef Full Text | Google Scholar

Dunand, C., Crèvecoeur, M., and Penel, C. (2007). Distribution of superoxide and hydrogen peroxide in Arabidopsis root and their influence on root development: possible interaction with peroxidases. N. Phytol. 174, 332–341. doi: 10.1111/j.1469-8137.2007.01995.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Eales, J. G. (1997). Iodine metabolism and thyroid-related functions in organisms lacking thyroid follicles: Are thyroid hormones also vitamins? Exp. Biol. Med. 214, 302–317. doi: 10.3181/00379727-214-44098

PubMed Abstract | CrossRef Full Text | Google Scholar

Foyer, C. H. (2018). Reactive oxygen species, oxidative signaling and the regulation of photosynthesis. Environ. Exp. Bot. 154, 134–142. doi: 10.1016/j.envexpbot.2018.05.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Fuge, R., and Johnson, C. C. (1986). The geochemistry of iodine - a review. Environ. Geochem. Health 8, 31–54. doi: 10.1007/BF02311063

PubMed Abstract | CrossRef Full Text | Google Scholar

Galetskiy, D., Lohscheider, J. N., Kononikhin, A. S., Popov, I. A., Nikolaev, E. N., and Adamska, I. (2011). Mass spectrometric characterization of photooxidative protein modifications in Arabidopsis thaliana thylakoid membranes. Rapid Commun. Mass Spectrom 25, 184–190. doi: 10.1002/rcm.4855

PubMed Abstract | CrossRef Full Text | Google Scholar

Gonzali, S., Kiferle, C., and Perata, P. (2017). Iodine biofortification of crops: agronomic biofortification, metabolic engineering and iodine bioavailability. Curr. Opin. Biotechnol. 44, 16–26. doi: 10.1016/j.copbio.2016.10.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Gupta, N., Bajpai, M., Majumdar, R., and Mishra, P. (2015). Response of iodine on antioxidant levels of Glycine max L. grown under Cadmium stress. Adv. Biol. Res. 9, 40–48. doi: 10.5829/idosi.abr.2015.9.1.9183

CrossRef Full Text | Google Scholar

Halka, M., Smoleń, S., Czernicka, M., Klimek-Chodacka, M., Pitala, J., and Tutaj, K. (2019). Iodine biofortification through expression of HMT, SAMT and S3H genes in Solanum lycopersicum L. Plant Physiol. Biochem. 144, 35–48. doi: 10.1016/j.plaphy.2019.09.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Hamblin, M. R., and Abrahamse, H. (2018). Inorganic salts and antimicrobial photodynamic therapy: mechanistic conundrums? Molecules 23, 1–18. doi: 10.3390/molecules23123190

PubMed Abstract | CrossRef Full Text | Google Scholar

Hiraga, S., Sasaki, K., Ito, H., Ohashi, Y., and Matsui, H. (2001). A large family of class III plant peroxidases. Plant Cell Physiol. 42, 462–468. doi: 10.1093/pcp/pce061

PubMed Abstract | CrossRef Full Text | Google Scholar

Hou, X., Yan, X., and Chai, C. (2000). Chemical species of iodine in some seaweeds II. Iodine-bound biological macromolecules. J. Radioanal. Nucl. Chem. 245, 461–467. doi: 10.1023/A:1006728420096

CrossRef Full Text | Google Scholar

Incrocci, L., Carmassi, G., Maggini, R., Poli, C., Saidov, D., Tamburini, C., et al. (2019). Iodine accumulation and tolerance in Sweet basil (Ocimum basilicum L.) with green or purple leaves grown in floating system technique. Front. Plant. Sci. 10:1–15. doi: 10.3389/fpls.2019.01494

PubMed Abstract | CrossRef Full Text | Google Scholar

Ionescu, I. A., Møller, B. L., and Sánchez-Pérez, R. (2017). Chemical control of flowering time. J. Experim. Bot. 68, 369–382. doi: 10.1093/jxb/erw427

PubMed Abstract | CrossRef Full Text | Google Scholar

Jouili, H., Bouazizi, H., and El Ferjani, E. (2011). Plant peroxidases: biomarkers of metallic stress. Acta Physiol. Plant. 33, 2075–2082. doi: 10.1007/s11738-011-0780-2

CrossRef Full Text | Google Scholar

Kiferle, C., Ascrizzi, R., Martinelli, M., Gonzali, S., Mariotti, L., Pistelli, L., et al. (2019). Effect of Iodine treatments on Ocimum basilicum L.: Biofortification, phenolics production and essential oil composition. PLoS One 14:e0226559. doi: 10.1371/journal.pone.0226559

PubMed Abstract | CrossRef Full Text | Google Scholar

Kiferle, C., Martinelli, M., Salzano, A. M., Gonzali, S., Beltrami, S., Salvadori, P. A., et al. (2020). Evidences for a nutritional role of iodine in plants. BiorXiv preprint.

Google Scholar

Kouzai, Y., Noutoshi, Y., Inoue, K., Shimizu, M., Onda, Y., and Mochida, K. (2018). Benzothiadiazole, a plant defense inducer, negatively regulates sheath blight resistance in Brachypodium distachyon. Sci. Rep. 8, 1–7. doi: 10.1038/s41598-018-35790-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Lehr, J. J., Wybenga, J. M., and Rosanow, M. (1958). Iodine as a micronutrient for tomatoes. Plant Physiol. 33, 421–427. doi: 10.1104/pp.33.6.421

PubMed Abstract | CrossRef Full Text | Google Scholar

Leyva, R., Sánchez-Rodríguez, E., Ríos, J. J., Rubio-Wilhelmi, M. M., Romero, L., Ruiz, J. M., et al. (2011). Beneficial effects of exogenous iodine in lettuce plants subjected to salinity stress. Plant Sci. 181, 195–202. doi: 10.1016/j.plantsci.2011.05.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, L., Aro, E. M., and Millar, A. H. (2018). Mechanisms of photodamage and protein turnover in photoinhibition. Trends Plant Sci. 23, 667–676. doi: 10.1016/j.tplants.2018.05.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, R., Wang, Z. J., Wang, L., Ma, B. C., Ghasimi, S., Lu, H., et al. (2016). Photocatalytic selective bromination of electron-rich aromatic compounds using microporous organic polymers with visible light. ACS Catal. 6, 1113–1121. doi: 10.1021/acscatal.5b02490

CrossRef Full Text | Google Scholar

Loreti, E., Poggi, A., Novi, G., Alpi, A., and Perata, P. (2005). A genome-wide analysis of the effects of sucrose on gene expression in arabidopsis seedlings under anoxia. Plant Physiol. 137, 1130–1138. doi: 10.1104/pp.104.057299

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, Y., and Yao, J. (2018). Chloroplasts at the crossroad of photosynthesis, pathogen infection and plant defense. Int. J. Mol. Sci. 19, 1–37. doi: 10.3390/ijms19123900

PubMed Abstract | CrossRef Full Text | Google Scholar

Machold, O., and Aurich, O. (1981). Thylakoid membrane proteins of higher plant chloroplasts II. Application of lactoperoxidase catalyzed iodination to study species specific differences and membrane topography. Biochem. Physiol. Pflanz. 176, 549–561. doi: 10.1016/s0015-3796(81)80028-5

CrossRef Full Text | Google Scholar

Medrano-Macías, J., Leija-Martínez, P., González-Morales, S., Juárez-Maldonado, A., and Benavides-Mendoza, A. (2016). Use of iodine to biofortify and promote growth and stress tolerance in crops. Front. Plant Sci. 7:1–20. doi: 10.3389/fpls.2016.01146

PubMed Abstract | CrossRef Full Text | Google Scholar

Mengel, K., and Kirkby, E. A. (2001). “Plant Nutrients,” in Principles of Plant Nutrition, eds K. Mengel, E. A. Kirkby, H. Kosegarten, and T. Appel (Dordrecht: Kluwer Academic Publishers), 1–13.

Google Scholar

Mittler, R., Vanderauwera, S., Gollery, M., and Van Breusegem, F. (2004). Reactive oxygen gene network of plants. Trends Plant Sci. 9, 490–498. doi: 10.1016/j.tplants.2004.08.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Moerschbacher, B. M. (1992). “Plant peroxidases involvement in response to pathogen,” in Plant Peroxidases 1980–1990, eds C. Panel, T. Gaspar, and H. Greppin (Geneva: University of Geneva Press), 91–99.

Google Scholar

Murashige, T., and Skoog, F. (1962). A revised medium for rapid growth and bio assays with tobacco tissue cultures. Physiol. Plant 15, 473–497. doi: 10.1111/j.1399-3054.1962.tb08052.x

CrossRef Full Text | Google Scholar

Mustroph, A., Lee, S. C., Oosumi, T., Zanetti, M. E., Yang, H., Ma, K., et al. (2010). Cross-kingdom comparison of transcriptomic adjustments to low-oxygen stress highlights conserved and plant-specific responses. Plant Physiol. 152, 1484–1500. doi: 10.1104/pp.109.151845

PubMed Abstract | CrossRef Full Text | Google Scholar

Nevo, R., Charuvi, D., Tsabari, O., and Reich, Z. (2012). Composition, architecture and dynamics of the photosynthetic apparatus in higher plants. Plant J. 70, 157–176. doi: 10.1111/j.1365-313X.2011.04876.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Pandian, B. A., Sathishraj, R., Djanaguiraman, M., Prasad, P. V. V., and Jugulam, M. (2020). Role of cytochrome P450 enzymes in plant stress response. Antioxid 9:454. doi: 10.3390/antiox9050454

CrossRef Full Text | Google Scholar

Perata, P., Matsukura, C., Vernieri, P., and Yamaguchi, J. (1997). Sugar repression of a gibberellin-dependent signaling pathway in barley embryos. Plant Cell 9, 2197–2208. doi: 10.1105/tpc.9.12.2197

PubMed Abstract | CrossRef Full Text | Google Scholar

Perez-Riverol, Y., Csordas, A., Bai, J., Bernal-Llinares, M., Hewapathirana, S., Kundu, D. J., et al. (2019). The PRIDE database and related tools and resources in 2019: Improving support for quantification data. Nucl. Acids Res. 47, D442–D450. doi: 10.1093/nar/gky1106

PubMed Abstract | CrossRef Full Text | Google Scholar

Pucciariello, C., Parlanti, S., Banti, V., Novi, G., and Perata, P. (2012). Reactive oxygen species-driven transcription in Arabidopsis under oxygen deprivation. Plant Physiol. 159, 184–196. doi: 10.1104/pp.111.191122

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramachandran, L. K. (1956). Protein-iodine interaction. Chem. Rev. 56, 199–218. doi: 10.1021/cr50008a001

CrossRef Full Text | Google Scholar

Romarís-Hortas, V., Bianga, J., Moreda-Piñeiro, A., Bermejo-Barrera, P., and Szpunar, J. (2014). Speciation of iodine-containing proteins in Nori seaweed by gel electrophoresis laser ablation ICP-MS. Talanta 127, 175–180. doi: 10.1016/j.talanta.2014.04.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmalenbach, I., Zhang, L., Reymond, M., and Jiménez-Gómez, J. M. (2014). The relationship between flowering time and growth responses to drought in the Arabidopsis Landsberg erecta x Antwerp-1 population. Front. Plant Sci. 5:1–9. doi: 10.3389/fpls.2014.00609

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmitt, F. J., Renger, G., Friedrich, T., Kreslavski, V. D., Zharmukhamedov, S. K., Los, D. A., et al. (2014). Reactive oxygen species: Re-evaluation of generation, monitoring and role in stress-signaling in phototrophic organisms. Biochim. Biophys. Acta Bioen. 1837, 835–848. doi: 10.1016/j.bbabio.2014.02.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Scott, R. (1954). Observations on the iodo-amino-acids of marine algae using iodine-131. Nature 173, 1098–1099. doi: 10.1038/1731098b0

CrossRef Full Text | Google Scholar

Smoleń, S., Kowalska, I., Halka, M., Ledwozyw-Smoleń, I., Grzanka, M., and Skoczylas, Ł et al. (2020). Selected aspects of iodate and iodosalicylate metabolism in lettuce including the activity of vanadium dependent haloperoxidases as affected by exogenous vanadium. Agronomy 10:10010001. doi: 10.3390/agronomy10010001

CrossRef Full Text | Google Scholar

Sonneveld, C. (2002). “Composition of nutrient solutions,” in Hydroponic Production of Vegetables and Ornamentals, eds D. Savvas and H. C. Passam (Athens: Embryo Publications), 179–210.

Google Scholar

Trapalis, M., Li, S. F., and Parish, R. W. (2017). The Arabidopsis GASA10 gene encodes a cell wall protein strongly expressed in developing anthers and seeds. Plant Sci. 260, 71–79. doi: 10.1016/j.plantsci.2017.04.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Umaly, R. C., and Poel, L. W. (1970). Effects of various concentrations of iodine as potassium iodide on the growth of barley, tomato and pea in nutrient solution culture. Ann. Bot. 34, 919–926. doi: 10.1093/oxfordjournals.aob.a084423

CrossRef Full Text | Google Scholar

Vlasova, I. I. (2018). Peroxidase activity of human hemoproteins: keeping the fire under control. Molecules 23, 1–27. doi: 10.3390/molecules23102561

PubMed Abstract | CrossRef Full Text | Google Scholar

Vlot, A. C., Dempsey, D. A., and Klessig, D. F. (2009). Salicylic acid, a multifaceted hormone to combat disease. Annu. Rev. Phytopathol. 47, 177–206. doi: 10.1146/annurev.phyto.050908.135202

PubMed Abstract | CrossRef Full Text | Google Scholar

Voogt, W., Holwerda, H. T., and Khodabaks, R. (2010). Biofortification of lettuce (Lactuca sativa L.) with iodine: the effect of iodine form and concentration in the nutrient solution on growth, development and iodine uptake of lettuce grown in water culture. J. Sci. Food Agric. 90, 906–913. doi: 10.1002/jsfa.3902

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Zhou, X., Fredimoses, M., Liao, S., and Liu, Y. (2014). Naturally occurring organoiodines. RSC Adv. 4, 57350–57376. doi: 10.1039/c4ra09833a

CrossRef Full Text | Google Scholar

Welch, R. M., and Shuman, L. (1995). Micronutrient Nutrition of Plants. CRC Crit. Rev. Plant Sci. 14, 49–82. doi: 10.1080/07352689509701922

CrossRef Full Text | Google Scholar

Wen, X., Zhang, X., Szewczyk, G., El-Hussein, A., Huang, Y. Y., Sarna, T., et al. (2017). Potassium iodide potentiates antimicrobial photodynamic inactivation mediated by Rose bengal in in vitro and in vivo studies. Antimicrob. Agents Chemother. 61:17. doi: 10.1128/AAC.00467-17

PubMed Abstract | CrossRef Full Text | Google Scholar

Weng, H. X., Hong, C. L., Xia, T. H., Bao, L. T., Liu, H. P., and Li, D. W. (2013). Iodine biofortification of vegetable plants-An innovative method for iodine supplementation. Chin. Sci. Bull. 58, 2066–2072. doi: 10.1007/s11434-013-5709-2

CrossRef Full Text | Google Scholar

Weng, H. X., Yan, A. L., Hong, C. L., Xie, L. L., Qin, Y. C., and Cheng, C. Q. (2008). Uptake of different species of iodine by water spinach and its effect to growth. Biol. Trace Elem. Res. 124, 184–194.

Google Scholar

Xing, J. H., Wang, F. R., Jia, J., Zhang, J., Li, L., Chen, Z., et al. (2013). Mutations in AP22.65 accelerate flowering in Arabidopsis thaliana. J. Plant Biol. 56, 49–58. doi: 10.1007/s12374-012-0379-1

CrossRef Full Text | Google Scholar

Yamamoto, Y., Kai, S., Ohnishi, A., Tsumura, N., Ishikawa, T., Hori, H., et al. (2014). Quality control of PSII: Behavior of PSII in the highly crowded grana thylakoids under excessive light. Plant Cell Physiol. 55, 1206–1215. doi: 10.1093/pcp/pcu043

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, Y. G., Huang, Y. Z., Hu, Y., and Liu, Y. X. (2003). Iodine uptake by spinach (Spinacia oleracea L.) plants grown in solution culture: Effects of iodine species and solution concentrations. Environ. Int. 29, 33–37. doi: 10.1016/S0160-4120(02)00129-0

CrossRef Full Text | Google Scholar

Zimmermann, M. B., Jooste, P. L., and Pandav, C. S. (2008). Iodine-deficiency disorders. Lancet 372, 1251–1262. doi: 10.1016/S0140-6736(08)61005-3

CrossRef Full Text | Google Scholar

Keywords: iodine, Arabidopsis, plant growth, plant nutrition, plant phenotype, transcriptomics, proteomics

Citation: Kiferle C, Martinelli M, Salzano AM, Gonzali S, Beltrami S, Salvadori PA, Hora K, Holwerda HT, Scaloni A and Perata P (2021) Evidences for a Nutritional Role of Iodine in Plants. Front. Plant Sci. 12:616868. doi: 10.3389/fpls.2021.616868

Received: 13 October 2020; Accepted: 04 January 2021;
Published: 17 February 2021.

Edited by:

Antonios Chrysargyris, Cyprus University of Technology, Cyprus

Reviewed by:

Artem Kasianov, Russian Academy of Sciences, Russia
Iwona Ledwożyw-Smoleń, University of Agriculture in Krakow, Poland

Copyright © 2021 Kiferle, Martinelli, Salzano, Gonzali, Beltrami, Salvadori, Hora, Holwerda, Scaloni and Perata. 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) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Claudia Kiferle, c.kiferle@santannapisa.it; c.kiferle@sssup.it; Pierdomenico Perata, p.perata@santannapisa.it

These authors have contributed equally to this work