# 2000 years of agriculture in the Atacama desert lead to changes in the distribution and concentration of iron in maize

Ale Elgeta Vidal, Nathalia Navarro, Mauricio Uribe, Kevin Robe, Frédéric Gaymard, Christian Dubos, María Fernanda Pérez, Hannetz Roschzttardtz

► To cite this version:

Ale Elgeta Vidal, Nathalia Navarro, Mauricio Uribe, Kevin Robe, Frédéric Gaymard, et al.. 2000 years of agriculture in the Atacama desert lead to changes in the distribution and concentration of iron in maize. Scientific Reports, 2021, 11 (1), 10.1038/s41598-021-96819-1. hal-03358782

HAL Id: hal-03358782

https://hal.inrae.fr/hal-03358782

Submitted on 29 Sep 2021

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_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

# 2000 years of agriculture in the Atacama desert lead to changes in the distribution and concentration of iron in maize

Ale Vidal Elgueta<fnref n="1" />,²✉, Nathalia Navarro³, Mauricio Uribe⁴, Kevin Robe⁵, Frédéric Gaymard⁵, Christian Dubos⁵, María Fernanda Pérez<fnref n="1" /> & Hannetz Roschzttardtz³✉

We performed a histological and quantitative study of iron in archaeological maize seeds from prehispanic times recovered from Tarapacá, Atacama Desert. Also, we examined iron distribution changes at the cell level in embryos from ancient versus new varieties of maize. Our results show a progressive decrease in iron concentration from the oldest maize to modern specimens. We interpret the results as an effect of prehispanic agriculture over the micronutrient composition of maize.

There is no consensus among agronomic, nutritional, and chemical sciences on the changes suffered by micro-nutrients such as iron, zinc, copper, and manganese in cereal crops<fnref n="1" />. The unequal conditions of the experiments, difficult comparisons², discrepancies in statistical methods³, and the lack of historical data⁴ causes these contradictory positions. The controversy deepens when it comes to the iron accumulation in maize because iron is quite variable (9.6 to 63.2 mg kg⁻<fnref n="1" />) in modern maize cultivars⁵. Moreover, despite the nutritional importance of maize today, with more than 713 million tons of annual production⁶, there are limited studies on the physiological changes associated with the absence or presence of micronutrients throughout the history of maize's agronomic production.

Intensive agriculture for thousands of years among prehispanic agricultural communities in Tarapacá, Atacama Desert (lat. 19°00′ 21°00′S), led to changes in maize morphology. Thus, our previous studies⁷ indicated that Tarapacá prehispanic farmers modified phenotypically their ancient maize. These modifications included the significant augmentation of cobs and kernels size from a mean length of 5 cm to more than 15 cm length for cobs. Kernel size augmentation ranges from 5 mm to more than 10 mm in length. Accordingly, Tarapacá maize is clearly differentiated between two phenotypic groups: a first group dated ca. 400 BC–500 AD and a later group dated ca. 500 AD to modern times (Supplementary Fig. 1). Tarapacá maize also shows a low genetic diversity caused by human selection and genetic similarities between the prehispanic and modern maize⁷.

A more recent process of plant selection has been studied in terms of nutritional quality. The Green Revolution had a negative impact in the total amount micronutrients in cereals⁸⁻¹⁰. Until date, the nature of these changes is still unknown. Then, the opportunity to evaluate the micronutrient distributions and accumulation directed over ancient maize could determine if all agricultural selection, indistinctive of time and techniques, results in the loss of micronutrients in cereals.

Additionally, different methods and plant models have been used to determine where micronutrients accumulate in seeds. Focusing on iron as a model, it has been shown, by a histological approach, that this essential micronutrient accumulates in vacuoles of a specific cell layer in the *Arabidopsis thaliana* mature dry embryos¹¹. Iron distribution in maize has not been studied extensively. Recently, using Perls staining on maize showed that iron accumulates principally in the embryo close to the endosperm face¹².

In this article, we used histological and quantitative methods to observe changes over iron distribution and iron concentration in archaeological and modern maize in continuous agricultural production in Tarapacá, Atacama Desert. Thus, we asked whether the loss of iron in crops could also be observed in archaeological maize, which was under a constant agriculture selection in prehispanic times.

<footnote n="1">¹Departamento de Ecología, Pontificia Universidad Católica de Chile, 8331150 Santiago, Chile. ²Escuela de Antropología, Pontificia Universidad Católica de Chile, Santiago, Chile. ³Departamento de Genética Molecular y Microbiología, Pontificia Universidad Católica de Chile, 8331150 Santiago, Chile. ⁴Departamento de Antropología, Universidad de Chile, 7800284 Santiago, Chile. ⁵BPMP, Univ Montpellier, CNRS, INRAE, Institut Agro, Montpellier, France. ✉email: aevidal@uc.cl; hroschzttardtz@bio.puc.cl</footnote>
## Results

**Dating of archaeological kernels.** The kernels used in this analysis came from seven archaeological sites. All the locations are in the same area of the Atacama Desert, close to each other, in the Tarapacá region (see Material and Methods and supplementary Fig. 1.).

To ensure the dates of Tarapacá maize we performed 19 radiocarbon dates directly on maize kernels and cobs from various sites of Tarapacá, including 10 dates for the samples reported in this study. The distribution of these calibrated age ranges is from 376 BC to 1799 AD (95.4% of probability) (Fig. 1 and Supplementary Table 1). This range covers a 2000 years' time span including five different periods of the Tarapacá culture (from the Early Formative Period; ca. 1000 BC-500 AD, to Colonial Period; ca.1541–1810 AD; Fig. 1).

The most ancient sample recovered for study is a popcorn kernel from Ramaditas site, dated at 376-206 cal. BC (95.4% of probability; Fig. 1). These results show for certain that maize was already present at the southern part of Tarapacá at least from ca. 400 BC and its production continued without interruption until ca. 600 AD. Then, we noticed a disruption in maize around 600 AD to 900 AD, which could be due to a sample methodology bias, but almost for sure it is due to the temporary discontinuance of agriculture. Thus, other organic material also lacks dates for these moments. Radiocarbon dates indicate a return to agricultural maize practices around 1100-1200 AD that continues intermittently to present times (Supplementary Fig. 1).

**Iron quantification and distribution in the archaeological kernels.** Quantitative analysis of iron content in individual kernels from the six archaeological sites and the modern variety is shown in Fig. 1. The kernels from Ramaditas site (376-206 cal. BC, Fig. 1) reached a high mean value of total iron concentration 36.6 ppm, while kernels from Alero Cerro Colorado 7 site (dated for seventeenth century, Fig. 1), do not exceed a mean of 13.05 ppm, corresponding to less than 36% of total iron concentration compared with the kernels from Ramaditas (Fig. 1). Our results indicate that there is a progressive decline of iron total content in the old kernels versus new varieties indicated by the relation between iron concentration in kernels and radiocarbon dates of kernels ($R^2$ = 0.1297, p value < 0.043).

Histological sections for kernels from the six archeological sites and the modern variety were obtained and the cellular structures and iron distribution were evaluated (Fig. 2). The analysis of toluidine blue stained sections indicated that cellular structures are well preserved in all the samples, including cell wall, cytoplasmic contents and organelles like nuclei are easily identifiable. No nuclei were detected in the endosperm cells indicating that the analyzed kernels correspond to mature grains (Fig. 2A,C,E,G,I,K,M). Regarding iron localization, Perls/DAB iron staining was performed on thin sections (Fig. 2). Iron is detected in the cytoplasm and cytoplasmic structures of embryo cells in all samples, while endosperm does not accumulate detectable iron by this histological method. No staining is observed in apoplast, indicating that iron pools are intracellular in embryo. Interestingly, nuclei are strongly stained in embryo cells from older kernels (ca. 400 BC to ca. 600 AD; Fig. 2B,D,F). However, nuclear iron pools are no longer detected in embryo cells in samples from Tarapacá 13 ca. 1500 AD to actual kernels (Fig. 2H,J,L). It should be noted that the no detection of nuclear iron in the embryo cells in samples from ca. 1500 AD to actual kernels is not due to differences in the samples preservation, as indicated by the toluidine blue staining where nuclei are perfectly observable (Fig. 2G,I,K). Additionally, protodermis is less stained by Perls/DAB staining in the samples from ca. 1500 AD to actual kernels compared to the older samples (Figs. 1 and 2).

## Discussion

In this study, we showed that iron accumulates in the cytoplasm or cytoplasmic structures of maize embryo cells (Fig. 2). Interestingly, iron is not detected in nuclei of embryo cells in the kernels from ca. 1200 AD (Tarapacá 13, Tarapacá Viejo, Alero Cerro Colorado 7) to modern kernels, while iron pool is clearly observed in embryo cells from kernels older than ca. 600 AD (Ramaditas, Pintados, and Tarapacá 40; Fig. 2). Iron accumulated in nuclei has been described in the literature, but its chemical speciation and function in this organelle remain unknown¹³. It has been suggested that nuclei may be a reservoir of iron during seed development¹⁴. A dynamic of intracellular iron accumulation has been observed during embryo maturation stages in brassicaceae plants, however, further studies should be performed to address these fundamental biological questions about the physiological role of the intracellular iron pools¹³,¹⁴.

During the Early Formative Period (ca. 1000 BC–500 AD), the first agricultural display produced maize kernels with high contents of iron (Figs. 1 and 2). During this period, maize had a low degree of manipulation⁷,¹⁵. Thus, there is a low amount of macrobotanical maize remains in archaeological sites that do not exceed 3% of the total macrobotanicals recovered in the period's settlements¹⁶. Also, isotope studies (δ¹³C and δ¹⁵N) in human remains from the Tarapacá 40 cemetery indicate a gradual transition to the consumption of C⁴ plants among the Tarapacá groups¹⁷. It is then stated that maize would not have been the main crop during the Early Formative period. Accordingly, agricultural production is extensive but not intensive. We suggest the mechanisms of artificial selection are not fully operating yet. At the end of the Late Formative Period (ca. 500–1200 AD), cultivation fields were extended to the north and south of Tarapacá¹⁸, kernels and cobs increased in size⁷, people are consuming maize¹⁷, and we observed iron loss in maize embryos (this study). The agricultural production is clearly intensified, and thus we find storage pits and debris of cobs and kernels in different prehispanic villages such as Caserones, Lluga Túmulos, among others. We postulate that the manipulation of maize is beginning to operate, with no clear consequences yet to iron distribution, but with an important decrease in iron accumulation. However, the full disclosure of human selection is recognizable during a third moment of agriculture (Late Intermediate and Late period). By then, most of agricultural production is relegated to artificial terraces, over 2500 m.a.s.l., implementing a system of terraces that depended almost exclusively on rainfall. A few settlements remain in the Pampa, such as Tarapacá Viejo, Pica 8, Lluga Túmulos, and minor sites. The presence of maize in these villages remains steady, and in some cases, there is an intensification of production. The reason for this
IRON (FE) QUANTITIES AND RADIOCARBON DATES IN MAIZE

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Iron (Fe) quantities and radiocarbon dates in maize. Quantifications of iron concentration in kernels are shown in color dots. Each site has been assigned to a different color and a time period. Linear regression is expressed in the dotted line. Linear regression results are: F-statistic = 4.57(1.23), R² = 0.166, R² adjusted = 0.129, p-value = 0.0432. The regression was calculated based on the quantities of iron in kernels and we assigned an archaeological date to each archaeological site indicated in the X bar. Date was assigned as years before present. Thus, we used quantitative dates. Radiocarbon dates are indicated in color bars which are the sum of probabilities for 19 dates. Dates in color are the maize kernels for the same archaeological site used in Fe analysis. Dates in grey are reported for maize kernels from other sites for Tarapacá region. Dates range from ca. 400 BC to ca. 1799 AD.</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Iron distribution in kernels from archaeological sites and modern samples from Tarapacá revealed by Perls/DAB staining. (A–B,O) pictures from Ramaditas, (C,D) pictures from Pintados, (E,F) pictures from Tarapacá 40, (G,H) pictures from Tarapacá 13, (I,J) pictures from Tarapacá Viejo, (K,L) pictures from Alero Cerro Colorado 7, (M,N) pictures from modern sample. (A,C,E,G,I,K,M) correspond to toluidine blue stained sections, used to verify integrity of cells and subcellular structures in the section. (B,D,F,H,J,L,N) correspond to sections stained by Perls/DAB in order to detect iron, the black labeling corresponds to iron. (O) is a negative control of the Perls/DAB staining (without K-ferrocyanide). For all pictures, embryo and endosperm correspond to Em and End, respectively. Arrows indicate some nuclei (for simplicity only in A). Bars: 50 µm.</figcaption></figure>
shift in production is not yet fully known, though environmental fluctuation could be related. But despite the
environmental challenges, cobs and kernels reached their major size probably due to the prolonged human
selection and the use of seabird guano fertilizer¹⁵.

Concomitant to the human selection responsible for increase in maize size (7, Supplementary Fig. 1), our
results show that iron is rarely observed in nuclei of embryo cells, and we observed low iron total amount. There-
fore, the presence of iron in nuclei of maize and the large quantities of iron reported for this period, reflect a
moment in which maize is not yet being manipulated, probably the first entrance of maize into the region 2400 yr.
BP. Previously, it has been suggested that a "dilution effect" could be the cause of micronutrient declination, this
is the mechanism by which the increase in size of dry matter in cereals is not accompanied by the increment of
micronutrient²¹. Several authors agree that similar effects accompanied the development of modern cultivars²².
Though we cannot establish the exact mechanism by which maize lost their micronutrient qualities in prehispanic
times, our results suggest that a process of ancient agriculture in Tarapacá desert, which included the progressive
intensification of agricultural production and selection of grain size affected maize kernels nutritional quality.
It is also remarkable that a similar process has been observed as a consequence of the Green Revolution in the
last 50 years¹⁰. It has been reported for wheat, that the introduction of semi-dwarf, high-yielding cultivars in
the mid 1960's, coincides with a progressive decrease in micronutrients contents in the grain⁸. Nonetheless, no
information is available about how intensification of agricultural production has affected the micronutrients
accumulation in maize kernels during the Green Revolution, and our results suggest, like for other crops, that
human selection could negatively affect the nutritional quality of grains.

There are other consequences of this study. The corroboration of ancient maize varieties which naturally have
large amounts of iron and the potential deciphering of Tarapacá maize genetics linked to iron presence. Then,
changes in iron distribution at subcellular levels is an interesting focus of research to understand the causes of
micronutrients loss in kernels maize by agriculture practices and selection. It is in this regard that the use of
archaeological seeds could help to understand the collateral effects on nutritional quality of crops under long
term agricultural practices, which is a crucial trait and a global challenge for worldwide agriculture. From an
agronomic point of view, it seems essential to consider in the improvement programs, in addition to the total
content of nutrients, the eventual changes regarding where micronutrients are accumulated in the grains, includ-
ing at the subcellular level.

**Material and methods**

**Ethics statement.** All permissions or licenses were obtained in order to collect and to use archaeological
maize samples. All methods were carried out in accordance with relevant guidelines and regulations.

**Plant samples and radiocarbon dating.** Tarapacá region is situated in the northern part of Chile, and the hinterland is part of the Atacama Desert characterized as a hyper arid environment crossed by several basins. Today only two basins have permanent runoff (Tarapacá basin and Camiña basin). The archaeological sites selected for this study are all situated in the same cultural area. Thus, representing different moments of the material, social and cultural life of prehispanic communities of Tarapacá.

A total of 65 samples were analyzed. Fifty-nine archaeological maize samples were used as follows: 19 kernels were dated, 18 kernels were used for histological analysis, and 22 kernels were used for iron quantification. Additionally, 3 modern kernels were used for histological analysis, and other 3 more modern kernels were added for iron quantification. The criteria to choose these samples were the preservation of the material and the number of kernels available for each site. In some cases, samples were limited. Then, we prioritized radiocarbon dating samples. The sample belonged to the following periods and archeological sites: Early Formative Period (ca. 2000 BC-500 AD, Ramaditas), Late Formative Period (ca. 500-1200 AD, Pintados and Tarapacá 40), Late Intermediate Period (ca. 1200-1450 AD, Tarapacá 13), Late Period (ca. 1450-1541 AD, Tarapacá Viejo), Colonial Period (ca. 1520-1650 AD, Alero Cerro Colorado 7). The archaeological sites are four domestic contexts (Ramaditas, Pintados, Tarapacá 13, and Tarapacá Viejo), one funerary context (Tarapacá 40) and one ceremonial deposit (Alero Cerro Colorado 7). Archaeological kernels were cultivated locally near the sites mentioned and deposited as debris in the domestic contexts. For Tarapacá 40 and Alero Cerro Colorado 7 samples were deposited as offerings. The archaeological samples are corn types and floury types (Supplementary Fig. 1). Modern samples were collected from Camiña basin (lat. 19°18′S, lon. 69°25′W) in traditional cultivars in 2014. The modern landraces were identified as Capio Chileno and Harinoso Tarapaqueño, and one European cultivar from INRAE-France.

Radiocarbon dating was made by DirectAMS radiocarbon dating service laboratory, Seattle, Washington,
USA. The samples were calibrated using the version OXCAL v4.3.2²³ and SCHCal13 atmospheric curve²⁴. All the
technical parameters of the radiocarbon dating analysis are included in the Supplementary Table 1.

**Sample preparations and Perls/DAB staining.** Maize kernels were longitudinally sectioned using a blade before the fixation step. The fixation and Perls/DAB staining were performed following previous studies¹³,¹⁹,²⁰. Briefly, all samples were vacuum infiltrated with a fixation solution (2% w/v paraformaldehyde in 1 mM phosphate buffer pH 7.0) for 45 min and incubated for 16 h in the same solution. The fixated material was dehydrated with different concentration of ethanol (50%, 70%, 80%, 90%, 95% and 100%), then incubated 12 h with a solution of butanol/ethanol 1:1 (v/v), and finally incubated several hours with 100% butanol. The kernels were embedded in the Technovit 7100 resin (Kulzer) according to the manufacturer's instructions and thin Sects. (3 µm) were obtained. For Perls/DAB staining, slides were incubated with 2% (v/v) HCl and 2% (w/v) K-ferrocyanide (Sigma Aldrich) for 45 min. For the DAB intensification, each glass slide was washed with distilled water, and incubated in a methanol solution containing 0.01 M NaN₃ (Sigma Aldrich) and 0.3% (v/v) H₂O₂ (Merck) for 1 h, and then washed with 0.1 M phosphate buffer (pH 7.4). For the intensification reaction, the
glass slides were incubated in 0.1 M phosphate buffer (pH 7.4) solution containing 0.025% (w/v) DAB (Sigma), 0.005% (v/v) H₂O₂, and 0.005% (w/v) CoCl₂·6H₂O. To stop the reaction the slides were rinsed with distilled water. The samples were observed and photographed with a Nikon Eclipse 80i microscope.

**Iron quantification.** Iron quantification was performed according to Gao et al., 2020²⁵. Fifteen milligrams of ground material (dry weight) per sample (one kernel) was mixed with 750 mL of nitric oxide (65% [v/v]) and 250 mL of hydrogen peroxide (30% [v/v]) prior to homogenization. Following 10 min of incubation at room temperature, the samples were mineralized using the microwave digestion systems (Berghof). Once mineralized, the nitric oxide proportion present in the samples was adjusted to 5 to 10% of the final volume by adding ultrapure water. Iron content present in the samples was then measured by microwave plasma atomic emission spectroscopy (Agilent Technologies).

Received: 26 February 2021; Accepted: 17 August 2021

Published online: 27 August 2021

**References**

1. Marles, R. Mineral nutrient composition of vegetables, fruits and grains: The context of reports of apparent historical declines. *J. Food Compos. Anal.* **56**, 93–103 (2017).

2. Davis, D. Declines in iron content of foods. *Br. J. Nutr.* **109**, 2111 (2013).

3. Davis, D. Commentary on: "Historical variation in the mineral composition of edible horticultural products" (White, P. J and Broadley, M.R (2005) Journal of Horticultural Science & Biotechnology, 80, 660-667). *J. Horticult. Sci. Biotechnol.* **81**(3), 553–554 (2006).

4. Broadley, M. R., Mead, A. & White, P. J. Replay to Davis (2006) Commentary. *J. Horticult. Sci. Technol.* **81**(3), 554–555 (2006).

5. Teklic, T., Loncaric, Z., Kovacevic, V. & Singh, B. R. Metallic trace elements in cereal grain—a review: How much metal do we eat? *Food Energy Secur.* **2**(2), 81–95 (2013).

6. Ranum, P., Peña-Rosas, J. P. & García-Casal, M. N. Global maize production, utilization and consumption. *Ann N Y Acad Sci.* **1312**, 105–112 (2014).

7. Vidal Elgueta, A., Hinojosa, L. F., Pérez, M. F., Peralta, G. & Rodríguez, M. U. Genetic and phenotypic diversity in 2000 years old maize (*Zea mays* L.) samples from the Tarapacá region, Atacama Desert, Chile. *PLoS ONE* **14**(1), e0210369. https://doi.org/10.1371/journal.pone.0210369 (2019).

8. Fan, M. S., Fairweather, S., Polton, P., Dunham, S. & Mcrath, S. Evidence of decreasing mineral density in wheat grain over the last 160 years. *J. Trace Elem. Med. Biol.* **22**, 315–324 (2008).

9. McGrath, S. The effects of increasing yields on the macro- and microelement concentrations and offtakes in the grain of winter wheat. *J. Sci. Food Agric.* **36**, 1073–1083 (1985).

10. De Fries, R., Fanzo, J., Remans, R., Palm, C. & Wood, S. Metrics for land-scarce agriculture: Nutrient content must be better integrated into planning. *Science* **349**(6245), 238–240 (2015).

11. Roschztardtz, H., Conéjéro, G., Curie, C. & Mari, S. Identification of the endotermal vacuole as the iron storage compartment in the arabidopsis embryo. *Plant Physiol.* **151**, 1329–1338 (2009).

12. Zang, J. et al. Maize YSL2 is required for iron distribution and development in kernels. *J. Exp. Bot.* **71**, 5896–5910 (2020).

13. Roschztardtz, H. et al. Plant cell nucleolus as a hot spot for iron. *J. Biol. Chem.* **286**, 27863–27866 (2011).

14. Ibeas, M., Grant-Grant, S., Navarro, N., Perez, F. & Roschztardtz, H. Dynamic subcellular localization of iron during embryo development in Brassicaceae seeds. *Front. Plant Sci.* **8**, 2186 (2017).

15. Santana-Sagredo, F. et al. 'White gold' guano fertilizer drove agricultural intensification in the Atacama Desert from AD 1000. *Nat. Plants.* **7**, 152–158 (2021).

16. García, M. et al. Alimentos, tecnologías vegetales y paleoambiente en las aldeas formativas de la pampa del Tamarugal (ca. 900 a.C.–800 d.C.). *Estudios Atacameños.* **47**, 33–58 (2014).

17. Santana-Sagredo, F., Uribe, M., Herrera, M. J., Retamal, R. & Flores, S. Brief communication: Dietary practices in ancient populations from northern chile during the transition to agriculture (Tarapaca Region, 1000 BC-AD 900). *Am. J. Phys. Anthropol.* **158**(4), 751–758 (2014).

18. Santoro, C. M. et al. Continuities and discontinuities in the socio-environmental systems of the Atacama Desert during the last 13,000 years. *J. Anthropol. Archaeol.* **46**, 28–39 (2017).

19. Roschztardtz, H., Conejero, G., Curie, C. & Mari, S. Identification of the endodermal vacuole as the iron storage compartment in the arabidopsis embryo. *Plant Physiol.* **151**, 1329–1338 (2009).

20. Ibeas, M. et al. The diverse iron distribution in Eudicotyledoneae seeds: From Arabidopsis to Quinoa. *Front. Plant Sci.* **15**, 1985 (2019).

21. Davis, D., Epp, M. & Riordan, H. Changes in USDA food composition data for 43 Garden crops, 1950 to 1990. *J. Am. Coll. Nutr.* **23**(6), 669–682 (2004).

22. White, P. J. & Broadley, M. R. Historical variation in the mineral composition of edible horticultural products. *J. Horticult. Sci. Biotchnol.* **80**(6), 660–667 (2005).

23. Bronk Ramsey, C. Methods for summarizing radiocarbon datasets. *Radiocarbon* **59**(2), 1809–1833 (2017).

24. Hogg, A. G. et al. SHCali3 southern hemisphere calibration, 0–50,000 years cal BP. *Radiocarbon* **55**(4), 1889–1903 (2013).

25. Gao, F., Robe, K., Bettembourg, M., Navarro, N., Rofidal, V., Santoni, V., Gaymard, F., Vignols, F., Roschztardtz, H., Izquierdo, E., & Dubos, C. The transcription factor bHLH121 interacts with bHLH105 (ILR3) and its closest homologs to regulate iron homeostasis in arabidopsis. *Plant Cell*, **32**, 508–524 (2020).

**Acknowledgements**

FONDECYT 1160334, FONDECYT 1181829, FONDECYT 1171369, ANID-ECOS Project C18B04, Pontificia Universidad Católica de Chile, Museo Regional de Iquique, Facultad de Ciencias Sociales Universidad de Chile, Universidad Arturo Prat (Estación Canchones), Instituto de Investigaciones Arqueológicas, Museo Gustavo Le Paige. We are grateful to Marcela Vidal for the graphic design of Figure 1 and to Dr. Evandro Ferrada for critically reading the manuscript. We also thank Sandrine Chay (BPMP, SAME platform, Montpellier, France) for technical support with plant Fe determination.
**Author contributions**

Conceived and designed the experiments: A.V.E., P.M.F., H.R. Performed the experiments: A.V.E., N.N., R.K., H.R. Analyzed the data: A.V.E., U.M., G.F., D.C., P.M.F., H.R. Contributed reagents/materials/ analysis tools: A.V.E., R.K., G.F., D.C., P.M.F., R.H. Wrote the paper: A.V.E., H.R.

**Competing interests**

The authors declare no competing interests.

**Additional information**

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-021-96819-1.

**Correspondence** and requests for materials should be addressed to A.V.E. or H.R.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

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

© The Author(s) 2021, corrected publication 2021