Gonnelli, Marsili-Libelli, Baker and Gabbrielli (2000) fitted a two-equation model — a Richards sigmoidal growth function for root elongation plus a first-order saturating nickel-uptake equation — to in-vitro time courses of the nickel hyperaccumulator Alyssum bertolonii, and proposed as a mechanistic hypothesis that Ni uptake occurs mainly during the active growth phase, with shoot saturation (about 12,600–19,600 micrograms per gram dry weight) set by influx acting against a declining growth rate and a practical harvest/toxicity threshold near 12,000 micrograms per gram shoot Ni. The model deliberately treated growth and uptake as uncoupled, and the authors identified the true growth–uptake interaction, experimental confirmation of the hypothesis, field translation of the laboratory-calibrated parameters, and generalization to other hyperaccumulators and multi-harvest cropping as directions for further work. No published work builds and validates a two-way coupled growth–uptake model with toxicity feedback for a Ni hyperaccumulator, experimentally discriminates the growth-phase-restricted uptake hypothesis from toxicity and dilution alternatives, or validates the laboratory-calibrated parameters and thresholds under soil or field conditions. The integrated objective remains open.

Background

Phytoextraction is the phytoremediation technique in which plants remove metals from contaminated soil by taking them up into harvestable aboveground tissue; phytomining, now often called agromining, applies the same idea to ultramafic (serpentine-derived, Ni-rich, Mg-rich) soils, treating the harvested metal-laden biomass as a bio-ore. Nickel hyperaccumulators are plants that concentrate Ni in shoots to extreme levels — by the field's convention, more than 1,000 micrograms of Ni per gram of dry shoot tissue, and in the best species more than 10,000 — whereas ordinary plants accumulate toxic symptoms above roughly a few hundred. The model species Alyssum bertolonii (frequently placed in the genus Odontarrhena in recent taxonomy) is a serpentine-endemic plant of central Italy and a canonical Ni hyperaccumulator.

Gonnelli et al. (2000) grew A. bertolonii plantlets in vitro in liquid culture medium for roughly 30 days under control conditions and 10, 50, and 100 micromolar NiSO4, and jointly fitted a two-equation model to the measured time courses. The first equation is a Richards growth function, an empirically flexible sigmoidal growth law dX/dt=rX(1−(X/K)n)dX/dt=rX(1-(X/K)^n) for root length XX with asymptote KK and shape exponent nn, generalizing logistic growth; the second is a first-order saturating uptake equation in which the tissue Ni concentration CC rises toward a saturation level CsatC_{sat} at a rate governed by a coefficient KcK_c. Five parameters were estimated per treatment by flexible polyhedron search, a Nelder–Mead-type direct-search optimization that requires no derivatives. The fits gave root saturation near 6,800 micrograms per gram, shoot saturations of 12,635–19,604 micrograms per gram depending on Ni supply, and approach to saturation within roughly two to three weeks. Growth curves of exposed plantlets diverged from controls at about day 16, when shoot Ni was near 12,000 micrograms per gram, and the authors took that value as a practical threshold marking toxicity onset and the latest sensible harvest time.

The essential limitation is that growth and uptake were modeled as uncoupled: tissue Ni burden does not feed back on growth, and growth phase does not gate uptake — the two processes interact only through treatment-specific fitted parameters. The authors explicitly advanced, as a plausible explanation rather than a demonstrated result, the mechanistic hypothesis that Ni uptake occurs mainly during the active growth phase and that shoot saturation arises from continued Ni influx interacting with a declining growth rate. Here root-to-shoot translocation is the movement of Ni from roots into shoots, largely carried in the transpiration stream (the water flow driven by evaporation from leaves) through the xylem, the plant's water-conducting tissue; and the competing explanations to consider alongside the growth-phase hypothesis are a pure toxicity effect (Ni impairs growth independently of growth phase), dilution by growth (expanding biomass lowers the measured tissue concentration), and passive transpiration-driven accumulation (Ni reaches shoots in the water stream regardless of growth phase). They left open the true growth–uptake interaction, the causal transport steps (root-to-shoot translocation kinetics and intracellular sequestration of Ni, largely in leaf-cell vacuoles), the translation of the in-vitro parameters and the 12,000 micrograms per gram threshold to soil-grown plants and the field, generalization beyond a single species, and optimization of multi-harvest cropping systems.

By 2000, hyperaccumulation had long been known, field demonstrations of hyperaccumulator phytoextraction existed, and soil-side models of metal uptake by roots had been developed, but no validated model coupled whole-plant growth dynamics with metal uptake and toxicity for a hyperaccumulator. Robust, transferable, integrated uptake modeling for phytoextraction remains an unmet goal of the field.

Problem Statement

Turn the uncoupled, in-vitro, single-species model of Gonnelli et al. (2000) into a predictive framework for hyperaccumulator nickel phytoextraction by coupling growth and uptake in both directions, with the validation of that coupling planned as an experimental and field-test program. The four items below jointly define this objective — a mechanistically coupled, transferable model of hyperaccumulator phytoextraction: (1) Explicitly model and calibrate the two-way interaction between plant growth and metal uptake: represent a possible bidirectional coupling, in which tissue Ni burden and the resulting toxicity may feed back on growth and growth phase may gate uptake, with the toxicity-feedback and growth-phase-gating coefficients left free and identifiable for estimation and testing, replacing the source model's treatment of growth and uptake as independent processes joined only by treatment-specific parameters; if item (2) refutes growth-phase gating, the corresponding coefficient is to be estimated to zero rather than asserted. (2) Test the mechanistic hypothesis that Ni uptake in the hyperaccumulator occurs mainly during the active growth phase and that shoot saturation is set by continued influx interacting with a declining growth rate, using the coupled model together with a program of targeted experiments, and specify the experiments that would establish or refute the causal links, including the transport steps of root-to-shoot translocation kinetics and intracellular sequestration. (3) Translate the laboratory-calibrated parameter values and the harvest/toxicity threshold of roughly 12,000 micrograms per gram shoot Ni to soil-grown plants and field conditions, where environmental heterogeneity, soil chemistry, and whole-plant physiology over longer lifecycles may alter saturation levels, uptake rates, and toxicity thresholds. (4) As a stated extension, generalize the framework to other hyperaccumulator species, soil-grown plants, and the optimization of multi-harvest cropping systems. The coupled model is to be developed and calibrated against the published in-vitro time-course summaries and thresholds given in Background (together with any openly available joint growth–uptake data), and the targeted experiments of item (2) and the soil- and field-translation of item (3) are to be specified as an experimental protocol and a field-validation plan rather than executed as new measurements.

References

  1. C. Gonnelli, S. Marsili-Libelli, A. Baker, R. Gabbrielli, 'Assessing Plant Phytoextraction Potential Through Mathematical Modeling', International Journal of Phytoremediation 2(4):343–351 (2000). DOI: 10.1080/15226510008500043, https://doi.org/10.1080/15226510008500043
  2. S. N. Whiting, M. R. Broadley, P. J. White, 'Applying a solute transfer model to phytoextraction: Zinc acquisition by Thlaspi caerulescens', Plant and Soil 249:45–56 (2003). DOI: 10.1023/A:1022542725880, https://doi.org/10.1023/A:1022542725880
  3. A. Schnepf, M. L. Himmelbauer, M. Puschenreiter, T. Schrefl, E. Lombi, W. J. Fitz, W. Loiskandl, W. W. Wenzel, 'Model development for simulating the bioavailability of Ni to the hyperaccumulator Thlaspi goesingense', in Biogeochemistry of Trace Elements in the Rhizosphere (Elsevier, 2005), pp. 391–418. DOI: 10.1016/B978-044451997-9/50015-5, https://doi.org/10.1016/B978-044451997-9/50015-5
  4. J. E. Low, M. L. Whitfield Åslund, A. Rutter, B. A. Zeeb, 'Effect of Plant Age on PCB Accumulation by Cucurbita pepo ssp. pepo', Journal of Environmental Quality 39:245–250 (2010). DOI: 10.2134/jeq2009.0169, https://doi.org/10.2134/jeq2009.0169
  5. L. Trakal, D. Martínez-Fernández, M. Vítková, M. Komárek, 'Phytoextraction of Metals: Modeling Root Metal Uptake and Associated Processes', in Phytoremediation (Springer, 2015), pp. 69–83. DOI: 10.1007/978-3-319-10395-2_6, https://doi.org/10.1007/978-3-319-10395-2_6
  6. F. J. Pereira, E. M. de Castro, 'A novel phytoremediation index for aquatic systems contaminated by toxic elements', Bioremediation Journal 27:84–91 (2021, online first; print issue 2023). DOI: 10.1080/10889868.2021.1988896, https://doi.org/10.1080/10889868.2021.1988896
  7. P. Drenning, A. Enell, D. Berggren Kleja, Y. Volchko, J. Norrman, 'Development of simplified probabilistic models for predicting phytoextraction timeframes of soil contaminants: demonstration at the DDX-contaminated Kolleberga tree nursery in Sweden', Environmental Science and Pollution Research 31:40925–40940 (2024). DOI: 10.1007/s11356-024-33858-x, https://doi.org/10.1007/s11356-024-33858-x
  8. D. Coinchelin, D. Stemmelen, F. Bartoli, 'A simple model for estimating Ni availability and leaf Ni accumulation for the Ni-hyperaccumulator Leptoplax emarginata', Plant and Soil 374:131–147 (2013, online first; print 2014). DOI: 10.1007/s11104-013-1873-z, https://doi.org/10.1007/s11104-013-1873-z
  9. B. Robinson, J.-E. Fernández, P. Madejón, T. Marañón, J. M. Murillo, S. Green, B. Clothier, 'Phytoextraction: an assessment of biogeochemical and economic viability', Plant and Soil 249:117–125 (2003). DOI: 10.1023/A:1022586524971, https://doi.org/10.1023/A:1022586524971
  10. S. Dehghani, K. R. Zupfer, L. Vasiluk, M. D. Dutton, M. Bellantino-Perco, B. A. Hale, 'Modeling phytoremediation of aged soil Ni from anthropogenic deposition using Alyssum murale', Chemosphere 267:128861 (2021). DOI: 10.1016/j.chemosphere.2020.128861, https://doi.org/10.1016/j.chemosphere.2020.128861
  11. G. Brunetti, R. Kodešová, J. Šimůnek, 'Modeling the Translocation and Transformation of Chemicals in the Soil–Plant Continuum: A Dynamic Plant Uptake Module for the HYDRUS Model', Water Resources Research 55:8967–8989 (2019). DOI: 10.1029/2019WR025432, https://doi.org/10.1029/2019WR025432
  12. P. Cârdei, C. Tudora, V. Vlăduț, M. A. Pruteanu, I. Găgeanu, D. Cujbescu, D.-M. Bordean, N. Ungureanu, G. Ipate, O. D. Cristea, 'Mathematical Model to Simulate the Transfer of Heavy Metals from Soil to Plant', Sustainability 13(11):6157 (2021). DOI: 10.3390/su13116157, https://doi.org/10.3390/su13116157
  13. L. Shi, J. Li, K. N. Palansooriya, Y. Chen, D. Hou, E. Meers, D. C. W. Tsang, X. Wang, Y. S. Ok, 'Modeling phytoremediation of heavy metal contaminated soils through machine learning', Journal of Hazardous Materials 441:129904 (2023). DOI: 10.1016/j.jhazmat.2022.129904, https://doi.org/10.1016/j.jhazmat.2022.129904
  14. G. Fellet, T. Centofanti, R. L. Chaney, C. E. Green, 'NiO(s) (bunsenite) is not available to Alyssum species', Plant and Soil 319:219–223 (2009). DOI: 10.1007/s11104-008-9863-2, https://doi.org/10.1007/s11104-008-9863-2
  15. H. Küpper, E. Lombi, F.-J. Zhao, G. Wieshammer, S. P. McGrath, 'Cellular compartmentation of nickel in the hyperaccumulators Alyssum lesbiacum, Alyssum bertolonii and Thlaspi goesingense', Journal of Experimental Botany 52(365):2291–2300 (2001). DOI: 10.1093/jexbot/52.365.2291, https://doi.org/10.1093/jexbot/52.365.2291
  16. S. Corem, A. Carpaneto, P. Soliani, L. Cornara, F. Gambale, J. Scholz-Starke, 'Response to cytosolic nickel of Slow Vacuolar channels in the hyperaccumulator plant Alyssum bertolonii', European Biophysics Journal 38:495–501 (2009). DOI: 10.1007/s00249-008-0400-2, https://doi.org/10.1007/s00249-008-0400-2
  17. Y.-M. Li, R. L. Chaney, E. P. Brewer, J. S. Angle, J. Nelkin, 'Phytoextraction of Nickel and Cobalt by Hyperaccumulator Alyssum Species Grown on Nickel-Contaminated Soils', Environmental Science & Technology 37(7):1463–1468 (2003). DOI: 10.1021/es0208963, https://doi.org/10.1021/es0208963
  18. A. Bani, G. Echevarria, S. Sulçe, J. L. Morel, A. Mullai, 'In-situ phytoextraction of Ni by a native population of Alyssum murale on an ultramafic site (Albania)', Plant and Soil 293:79–89 (2007). DOI: 10.1007/s11104-007-9245-1, https://doi.org/10.1007/s11104-007-9245-1
  19. R. F. Saad, G. Echevarria, B. Rodríguez-Garrido, P. Kidd, E. Benizri, 'A two-year field study of nickel-agromining using Odontarrhena chalcidica co-cropped with a legume on an ultramafic soil: temporal variation in plant biomass, nickel yields and taxonomic and bacterial functional diversity', Plant and Soil 461:471–488 (2021). DOI: 10.1007/s11104-021-04834-y, https://doi.org/10.1007/s11104-021-04834-y
  20. R. R. Canales-Pastrana, M. Paredes, 'Phytoremediation Dynamic Model as an Assessment Tool in the Environmental Management', Open Journal of Applied Sciences 3(2):208–217 (2013). DOI: 10.4236/ojapps.2013.32028, https://doi.org/10.4236/ojapps.2013.32028