Register      Login
Functional Plant Biology Functional Plant Biology Society
Plant function and evolutionary biology
RESEARCH ARTICLE (Open Access)

Physiological and cannabinoid responses of hemp (Cannabis sativa) to rock phosphate dust under tropical conditions

Luca De Prato https://orcid.org/0000-0001-5203-4331 A B * , Omid Ansari https://orcid.org/0000-0001-9969-7044 C , Giles E. S. J. Hardy D E , John Howieson A , Graham O’Hara A and Katinka X. Ruthrof D F
+ Author Affiliations
- Author Affiliations

A Murdoch University, Food Futures Institute, Murdoch, WA, Australia.

B Medicann Health Aust Pty Ltd, Osborne Park, WA, Australia.

C HempGenTech Pty Ltd, Kenmore, Qld, Australia.

D Murdoch University, Harry Butler Institute, Murdoch, WA, Australia.

E ArborCarbon, Murdoch University, Murdoch, WA, Australia.

F Department of Biodiversity, Conservation and Attractions, Kensington, WA 6151, Australia.

* Correspondence to: L.DePrato@murdoch.edu.au

Handling Editor: Tim Cavagnaro

Functional Plant Biology 50(5) 378-389 https://doi.org/10.1071/FP22264
Submitted: 1 February 2022  Accepted: 27 February 2023   Published: 28 March 2023

© 2023 The Author(s) (or their employer(s)). Published by CSIRO Publishing. This is an open access article distributed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License (CC BY-NC-ND)

Abstract

Growing a high-value crop such as industrial hemp (Cannabis sativa L.) in post-mining environments is economically and environmentally attractive but faces a range of biotic and abiotic challenges. An opportunity to investigate the cultivation of C. sativa presented itself as part of post-mining activities on Christmas Island (Australia) to profitably utilise disused phosphate (PS) quarries. Challenges to plant growth and cadmium (Cd) uptake were addressed in this study using potted plants under fully controlled conditions in a growth chamber. A complete nutritional spectrum, slow-release fertiliser was applied to all plants as a control treatment, and two levels of rock PS dust, a waste product of PS mining that contains 35% phosphorus (P) and 40 ppm of naturally occurring Cd, were applied at 54 and 162 g L−1. After 12 weeks, control plants (no PS dust) significantly differed in phenological development, with no flower production, lower aboveground biomass and reduced photosynthesis efficiency than those with P applied as rock dust. Compared with the controls, the 54 g L−1 level of P dust increased shoot biomass by 38%, while 162 g L−1 increased shoot biomass by 85%. The concentration of Δ9-tetrahydrocannabinol also increased with the higher P levels. Cd uptake from PS dust by C. sativa was substantial and warrants further investigation. However, there was no increase in Cd content between the 54 and 162 g L−1 application rates in seed and leaf. Results indicate that hemp could become a high-value crop on Christmas Island, with the readily available rock PS dust providing a source of P.

Keywords: agriculture post-mining, cadmium, Cannabis sativa, flowering, heavy metal, industrial hemp, nutrition, photosynthesis.


References

Allan RJ (1995) Impact of mining activities on the terrestrial and aquatic environment with emphasis on mitigation and remedial measures. In ‘Heavy metals’. (Ed. U Förstner, W Salomons, P Mader) pp. 119–140. (Springer: Berlin, Heidelberg)

Angelova V, Ivanova R, Delibaltova V, Ivanov K (2004) Bio-accumulation and distribution of heavy metals in fibre crops (flax, cotton and hemp). Industrial Crops and Products 19, 197–205.
Bio-accumulation and distribution of heavy metals in fibre crops (flax, cotton and hemp).Crossref | GoogleScholarGoogle Scholar |

Aubin M-P, Seguin P, Vanasse A, Tremblay GF, Mustafa AF, Charron J-B (2015) Industrial hemp response to nitrogen, phosphorus, and potassium fertilization. Crop, Forage & Turfgrass Management 1, 1–10.
Industrial hemp response to nitrogen, phosphorus, and potassium fertilization.Crossref | GoogleScholarGoogle Scholar |

Bates D, Mächler M, Bolker B, Walker S (2014) Fitting linear mixed-effects models using lme4. Journal of Statistical Software 67, 48
Fitting linear mixed-effects models using lme4.Crossref | GoogleScholarGoogle Scholar |

Bernstein N, Gorelick J, Zerahia R, Koch S (2019) Impact of N, P, K, and humic acid supplementation on the chemical profile of medical Cannabis (Cannabis sativa L). Frontiers in Plant Science 10, 736
Impact of N, P, K, and humic acid supplementation on the chemical profile of medical Cannabis (Cannabis sativa L).Crossref | GoogleScholarGoogle Scholar |

Bhargava A, Carmona FF, Bhargava M, Srivastava S (2012) Approaches for enhanced phytoextraction of heavy metals. Journal of Environmental Management 105, 103–120.
Approaches for enhanced phytoextraction of heavy metals.Crossref | GoogleScholarGoogle Scholar |

Bolland MDA (1985) Effects of phosphorus on seed yields of subterranean clover, serradella and annual medics. Australian Journal of Experimental Agriculture 25, 595–602.
Effects of phosphorus on seed yields of subterranean clover, serradella and annual medics.Crossref | GoogleScholarGoogle Scholar |

Bócsa I, Máthé P, Hangyel L (1997) Effect of nitrogen on tetrahydrocannabinol (THC) content in hemp (Cannabis sativa L.) leaves at different positions. Journal of the International Hemp Association 4, 78–79.

Böhm C, Quinkenstein A, Freese D, Hüttl RF (2011) Assessing the short rotation woody biomass production on marginal post-mining areas. Journal of Forest Science 57, 303–311.
Assessing the short rotation woody biomass production on marginal post-mining areas.Crossref | GoogleScholarGoogle Scholar |

Calzolari D, Magagnini G, Lucini L, Grassi G, Appendino GB, Amaducci S (2017) High added-value compounds from Cannabis threshing residues. Industrial Crops and Products 108, 558–563.
High added-value compounds from Cannabis threshing residues.Crossref | GoogleScholarGoogle Scholar |

Cen H, Weng H, Yao J, He M, Lv J, Hua S, Li H, He Y (2017) Chlorophyll fluorescence imaging uncovers photosynthetic fingerprint of citrus Huanglongbing. Frontiers in Plant Science 8, 1509
Chlorophyll fluorescence imaging uncovers photosynthetic fingerprint of citrus Huanglongbing.Crossref | GoogleScholarGoogle Scholar |

Chiou T-J, Liu H, Harrison MJ (2001) The spatial expression patterns of a phosphate transporter (MtPT1) from Medicago truncatula indicate a role in phosphate transport at the root/soil interface. The Plant Journal 25, 281–293.
The spatial expression patterns of a phosphate transporter (MtPT1) from Medicago truncatula indicate a role in phosphate transport at the root/soil interface.Crossref | GoogleScholarGoogle Scholar |

Chong J, Wishart DS, Xia J (2019) Using metaboanalyst 4.0 for comprehensive and integrative metabolomics data analysis. Current Protocols in Bioinformatics 68, e86
Using metaboanalyst 4.0 for comprehensive and integrative metabolomics data analysis.Crossref | GoogleScholarGoogle Scholar |

Citterio S, Santagostino A, Fumagalli P, Prato N, Ranalli P, Sgorbati S (2003) Heavy metal tolerance and accumulation of Cd, Cr and Ni by Cannabis sativa L. Plant and Soil 256, 243–252.
Heavy metal tolerance and accumulation of Cd, Cr and Ni by Cannabis sativa L.Crossref | GoogleScholarGoogle Scholar |

Cockson P, Schroeder-Moreno M, Veazie P, Barajas G, Logan D, Davis M, Whipker BE (2020) Impact of phosphorus on Cannabis sativa reproduction, cannabinoids, and terpenes. Applied Sciences 10, 7875
Impact of phosphorus on Cannabis sativa reproduction, cannabinoids, and terpenes.Crossref | GoogleScholarGoogle Scholar |

De Prato L (2021) Genotype by environment interactions of industrial hemp (Cannabis sativa L.) varieties under tropical conditions. PhD thesis, Murdoch University.

De Prato L, Ansari O, Hardy GESJ, Howieson J, O’Hara G, Ruthrof KX (2022a) The cannabinoid profile and growth of hemp (Cannabis sativa L.) is influenced by tropical daylengths and temperatures, genotype and nitrogen nutrition. Industrial Crops and Products 178, 114605
The cannabinoid profile and growth of hemp (Cannabis sativa L.) is influenced by tropical daylengths and temperatures, genotype and nitrogen nutrition.Crossref | GoogleScholarGoogle Scholar |

De Prato L, Timmins M, Ansari O, Ruthrof KX, Hardy GESJ, Howieson J, O’Hara G (2022b) Semi-quantitative analysis of cannabinoids in hemp (Cannabis sativa L.) using gas chromatography coupled to mass spectrometry. Journal of Cannabis Research 4, 51
Semi-quantitative analysis of cannabinoids in hemp (Cannabis sativa L.) using gas chromatography coupled to mass spectrometry.Crossref | GoogleScholarGoogle Scholar |

De Prato L, Ansari O, Hardy GESJ, Howieson J, O’Hara G, Ruthrof KX (2022c) Morpho-physiology and cannabinoid concentrations of hemp (Cannabis sativa L.) are affected by potassium fertilisers and microbes under tropical conditions. Industrial Crops and Products 182, 114907
Morpho-physiology and cannabinoid concentrations of hemp (Cannabis sativa L.) are affected by potassium fertilisers and microbes under tropical conditions.Crossref | GoogleScholarGoogle Scholar |

Deikman J (1997) Molecular mechanisms of ethylene regulation of gene transcription. Physiologia Plantarum 100, 561–566.
Molecular mechanisms of ethylene regulation of gene transcription.Crossref | GoogleScholarGoogle Scholar |

Diarra I, Prasad S (2021) The current state of heavy metal pollution in Pacific Island Countries: a review. Applied Spectroscopy Reviews 56, 27–51.
The current state of heavy metal pollution in Pacific Island Countries: a review.Crossref | GoogleScholarGoogle Scholar |

Duursma RA (2015) Plantecophys – an R package for analysing and modelling leaf gas exchange data. PLoS ONE 10, e0143346
Plantecophys – an R package for analysing and modelling leaf gas exchange data.Crossref | GoogleScholarGoogle Scholar |

Elisa B, Marsano F, Cavaletto M, Berta G (2007) Copper stress in Cannabis sativa roots: morphological and proteomic analysis. Caryologia 60, 96–101.
Copper stress in Cannabis sativa roots: morphological and proteomic analysis.Crossref | GoogleScholarGoogle Scholar |

Foyer C, Spencer C (1986) The relationship between phosphate status and photosynthesis in leaves. Planta 167, 369–375.
The relationship between phosphate status and photosynthesis in leaves.Crossref | GoogleScholarGoogle Scholar |

Gamalero E, Lingua G, Berta G, Glick BR (2009) Beneficial role of plant growth promoting bacteria and arbuscular mycorrhizal fungi on plant responses to heavy metal stress. Canadian Journal of Microbiology 55, 501–514.
Beneficial role of plant growth promoting bacteria and arbuscular mycorrhizal fungi on plant responses to heavy metal stress.Crossref | GoogleScholarGoogle Scholar |

Godfray HCJ, Beddington JR, Crute IR, Haddad L, Lawrence D, Muir JF, Pretty J, Robinson S, Thomas SM, Toulmin C (2010) Food security: the challenge of feeding 9 billion people. Science 327, 812–818.
Food security: the challenge of feeding 9 billion people.Crossref | GoogleScholarGoogle Scholar |

Howieson J, Calmy H, Ballard N, Skinner P, WO’Hara G, Skinner L, Ruthrof KX, Swift R, Ballard V, St Hardy GE, McHenry MP (2017) Bread from stones: post-mining land use change from phosphate mining to farmland. The Extractive Industries and Society 4, 290–299.
Bread from stones: post-mining land use change from phosphate mining to farmland.Crossref | GoogleScholarGoogle Scholar |

Jin JB, Hasegawa PM (2008) Flowering time regulation by the SUMO E3 ligase SIZ1. Plant Signaling & Behavior 3, 891–892.
Flowering time regulation by the SUMO E3 ligase SIZ1.Crossref | GoogleScholarGoogle Scholar |

Lauer MJ, Blevins DG, Sierzputowska-Gracz H (1989) 31P-nuclear magnetic resonance determination of phosphate compartmentation in leaves of reproductive soybeans (Glycine max L.) as affected by phosphate nutrition. Plant Physiology 89, 1331–1336.
31P-nuclear magnetic resonance determination of phosphate compartmentation in leaves of reproductive soybeans (Glycine max L.) as affected by phosphate nutrition.Crossref | GoogleScholarGoogle Scholar |

Linger P, Müssig J, Fischer H, Kobert J (2002) Industrial hemp (Cannabis sativa L.) growing on heavy metal contaminated soil: fibre quality and phytoremediation potential. Industrial Crops and Products 16, 33–42.
Industrial hemp (Cannabis sativa L.) growing on heavy metal contaminated soil: fibre quality and phytoremediation potential.Crossref | GoogleScholarGoogle Scholar |

Linger P, Ostwald A, Haensler J (2005) Cannabis sativa L. growing on heavy metal contaminated soil: growth, cadmium uptake and photosynthesis. Biologia Plantarum 49, 567–576.
Cannabis sativa L. growing on heavy metal contaminated soil: growth, cadmium uptake and photosynthesis.Crossref | GoogleScholarGoogle Scholar |

Liu J, Uhde-Stone C, Li A, Vance C, Allan D (2001) A phosphate transporter with enhanced expression in proteoid roots of white lupin (Lupinus albus L.). Plant and Soil 237, 257–266.
A phosphate transporter with enhanced expression in proteoid roots of white lupin (Lupinus albus L.).Crossref | GoogleScholarGoogle Scholar |

López-Arredondo DL, Leyva-González MA, González-Morales SI, López-Bucio J, Herrera-Estrella L (2014) Phosphate nutrition: improving low-phosphate tolerance in crops. Annual Review of Plant Biology 65, 95–123.
Phosphate nutrition: improving low-phosphate tolerance in crops.Crossref | GoogleScholarGoogle Scholar |

Magagnini G, Grassi G, Kotiranta S (2018) The effect of light spectrum on the morphology and cannabinoid content of Cannabis sativa L. Medical Cannabis and Cannabinoids 1, 19–27.
The effect of light spectrum on the morphology and cannabinoid content of Cannabis sativa L.Crossref | GoogleScholarGoogle Scholar |

Malceva M, Vikmane M, Stramkale V (2011) Changes of photosynthesis-related parameters and productivity of Cannabis sativa under different nitrogen supply. Environmental and Experimental Biology 9, 61–69.

Maxwell K, Johnson GN (2000) Chlorophyll fluorescence – a practical guide. Journal of Experimental Botany 51, 659–668.
Chlorophyll fluorescence – a practical guide.Crossref | GoogleScholarGoogle Scholar |

Mediavilla V, Jonquera M, Schmid-Slembrouck I, Soldati A (1998) Decimal code for growth stages of hemp (Cannabis sativa L.). Journal of the International Hemp Association 5, 68–74.

Mengistu T, Sterck FJ, Anten NPR, Bongers F (2012) Frankincense tapping reduced photosynthetic carbon gain in Boswellia papyrifera (Burseraceae) trees. Forest Ecology and Management 278, 1–8.
Frankincense tapping reduced photosynthetic carbon gain in Boswellia papyrifera (Burseraceae) trees.Crossref | GoogleScholarGoogle Scholar |

Ozanne PG, Keay J, Biddiscombe EF (1969) The comparative applied phosphate requirements of eight annual pasture species. Australian Journal of Agricultural Research 20, 809–818.
The comparative applied phosphate requirements of eight annual pasture species.Crossref | GoogleScholarGoogle Scholar |

Petrová S, Benešová D, Soudek P, Vaněk T (2012) Enhancement of metal(loid)s phytoextraction by Cannabis sativa L. Journal of Food, Agriculture and Environment 10, 631–641.

Plaxton WC, Tran HT (2011) Metabolic adaptations of phosphate-starved plants. Plant Physiology 156, 1006–1015.
Metabolic adaptations of phosphate-starved plants.Crossref | GoogleScholarGoogle Scholar |

Plénet D, Etchebest S, Mollier A, Pellerin S (2000a) Growth analysis of maize field crops under phosphorus deficiency. Plant and Soil 223, 119–132.
Growth analysis of maize field crops under phosphorus deficiency.Crossref | GoogleScholarGoogle Scholar |

Plénet D, Mollier A, Pellerin S (2000b) Growth analysis of maize field crops under phosphorus deficiency. II. Radiation-use efficiency, biomass accumulation and yield components. Plant and Soil 224, 259–272.
Growth analysis of maize field crops under phosphorus deficiency. II. Radiation-use efficiency, biomass accumulation and yield components.Crossref | GoogleScholarGoogle Scholar |

Prasad MNV, Strzalka K (2013) ‘Physiology and biochemistry of metal toxicity and tolerance in plants.’ pp. 432. (Springer Netherlands: Dordrecht, NL)

RStudio Team (2020) RStudio: Integrated Development for R. RStudio, PBC, Boston, MA. Available at http://www.rstudio.com/

Ruthrof KX, Fontaine JB, Hopkins AJM, McHenry MP, O’Hara G, McComb J, Hardy GESJ, Howieson J (2018a) Potassium amendment increases biomass and reduces heavy metal concentrations in Lablab purpureus after phosphate mining. Land Degradation & Development 29, 398–407.
Potassium amendment increases biomass and reduces heavy metal concentrations in Lablab purpureus after phosphate mining.Crossref | GoogleScholarGoogle Scholar |

Ruthrof KX, Steel E, Misra S, McComb J, O’Hara G, Hardy GESJ, Howieson J (2018b) Transitioning from phosphate mining to agriculture: responses to urea and slow release fertilizers for Sorghum bicolor. Science of The Total Environment 625, 1–7.
Transitioning from phosphate mining to agriculture: responses to urea and slow release fertilizers for Sorghum bicolor.Crossref | GoogleScholarGoogle Scholar |

Ryan PR, Delhaize E, Jones DL (2001) Function and mechanism of organic anion exudation from plant roots. Annual Review of Plant Physiology and Plant Molecular Biology 52, 527–560.
Function and mechanism of organic anion exudation from plant roots.Crossref | GoogleScholarGoogle Scholar |

Saloner A, Bernstein N (2021) Nitrogen supply affects cannabinoid and terpenoid profile in medical cannabis (Cannabis sativa L.). Industrial Crops and Products 167, 113516
Nitrogen supply affects cannabinoid and terpenoid profile in medical cannabis (Cannabis sativa L.).Crossref | GoogleScholarGoogle Scholar |

Saloner A, Sacks MM, Bernstein N (2019) Response of medical cannabis (Cannabis sativa L.) genotypes to K supply under long photoperiod. Frontiers in Plant Science 10, 1369
Response of medical cannabis (Cannabis sativa L.) genotypes to K supply under long photoperiod.Crossref | GoogleScholarGoogle Scholar |

Sera B, Sery M, Gavril B, Gajdova I (2017) Seed germination and early growth responses to seed pre-treatment by non-thermal plasma in hemp cultivars (Cannabis sativa L.). Plasma Chemistry and Plasma Processing 37, 207–221.
Seed germination and early growth responses to seed pre-treatment by non-thermal plasma in hemp cultivars (Cannabis sativa L.).Crossref | GoogleScholarGoogle Scholar |

Shiponi S, Bernstein N (2021a) Response of medical cannabis (Cannabis sativa L.) genotypes to P supply under long photoperiod: functional phenotyping and the ionome. Industrial Crops and Products 161, 113154
Response of medical cannabis (Cannabis sativa L.) genotypes to P supply under long photoperiod: functional phenotyping and the ionome.Crossref | GoogleScholarGoogle Scholar |

Shiponi S, Bernstein N (2021b) The highs and lows of P supply in medical Cannabis: effects on cannabinoids, the ionome, and morpho-physiology. Frontiers in Plant Science 12, 657323
The highs and lows of P supply in medical Cannabis: effects on cannabinoids, the ionome, and morpho-physiology.Crossref | GoogleScholarGoogle Scholar |

Small E, Pocock T, Cavers PB (2003) The biology of Canadian weeds. 119. Cannabis sativa L. Canadian Journal of Plant Science 83, 217–237.
The biology of Canadian weeds. 119. Cannabis sativa L.Crossref | GoogleScholarGoogle Scholar |

Stirbet A, Govindjee (2011) On the relation between the Kautsky effect (chlorophyll a fluorescence induction) and photosystem II: basics and applications of the OJIP fluorescence transient. Journal of Photochemistry and Photobiology B: Biology 104, 236–257.
On the relation between the Kautsky effect (chlorophyll a fluorescence induction) and photosystem II: basics and applications of the OJIP fluorescence transient.Crossref | GoogleScholarGoogle Scholar |

Strasser RJ, Srivastava A, Tsimilli-Michael M (2000) The fluorescence transient as a tool to characterize and screen photosynthetic samples. In ‘Probing photosynthesis: mechanisms, regulation and adaptation’. (Eds M Yunus, U Pathre, P Mohanty) pp. 445–483. (Taylor & Francis Group: Boca Raton, FL, USA)

Tang K, Struik PC, Yin X, Thouminot C, Bjelková M, Stramkale V, Amaducci S (2016) Comparing hemp (Cannabis sativa L.) cultivars for dual-purpose production under contrasting environments. Industrial Crops and Products 87, 33–44.
Comparing hemp (Cannabis sativa L.) cultivars for dual-purpose production under contrasting environments.Crossref | GoogleScholarGoogle Scholar |

Tang K, Struik PC, Amaducci S, Stomph T-J, Yin X (2017) Hemp (Cannabis sativa L.) leaf photosynthesis in relation to nitrogen content and temperature: implications for hemp as a bio-economically sustainable crop. GCB Bioenergy 9, 1573–1587.
Hemp (Cannabis sativa L.) leaf photosynthesis in relation to nitrogen content and temperature: implications for hemp as a bio-economically sustainable crop.Crossref | GoogleScholarGoogle Scholar |

Thorsteinsson B, Tillberg J-E (1987) Carbohydrate partitioning, photosynthesis and growth in Lemna gibba G3. II. Effects of phosphorus limitation. Physiologia Plantarum 71, 271–276.
Carbohydrate partitioning, photosynthesis and growth in Lemna gibba G3. II. Effects of phosphorus limitation.Crossref | GoogleScholarGoogle Scholar |

Valle JR, Vieira JE, Aucélio JG, Valio IF (1978) Influence of photoperiodism on cannabinoid content of Cannabis sativa L. Bulletin on Narcotics 30, 67–68.

Vance CP, Uhde-Stone C, Allan DL (2003) Phosphorus acquisition and use: critical adaptations by plants for securing a nonrenewable resource. New Phytologist 157, 423–447.
Phosphorus acquisition and use: critical adaptations by plants for securing a nonrenewable resource.Crossref | GoogleScholarGoogle Scholar |

Vera CL, Malhi SS, Raney JP, Wang ZH (2004) The effect of N and P fertilization on growth, seed yield and quality of industrial hemp in the Parkland region of Saskatchewan. Canadian Journal of Plant Science 84, 939–947.
The effect of N and P fertilization on growth, seed yield and quality of industrial hemp in the Parkland region of Saskatchewan.Crossref | GoogleScholarGoogle Scholar |

Vera CL, Malhi SS, Phelps SM, May WE, Johnson EN (2010) N, P, and S fertilization effects on industrial hemp in Saskatchewan. Canadian Journal of Plant Science 90, 179–184.
N, P, and S fertilization effects on industrial hemp in Saskatchewan.Crossref | GoogleScholarGoogle Scholar |

Wickham H (2011) ggplot2. Wiley Interdisciplinary Reviews: Computational Statistics 3, 180–185.
ggplot2.Crossref | GoogleScholarGoogle Scholar |

Xiong Y, DeFraia C, Williams D, Zhang X, Mou Z (2009) Deficiency in a cytosolic ribose-5-phosphate isomerase causes chloroplast dysfunction, late flowering and premature cell death in Arabidopsis. Physiologia Plantarum 137, 249–263.
Deficiency in a cytosolic ribose-5-phosphate isomerase causes chloroplast dysfunction, late flowering and premature cell death in Arabidopsis.Crossref | GoogleScholarGoogle Scholar |

Yu T-S, Lue W-L, Wang S-M, Chen J (2000) Mutation of Arabidopsis plastid phosphoglucose isomerase affects leaf starch synthesis and floral initiation. Plant Physiology 123, 319–326.
Mutation of Arabidopsis plastid phosphoglucose isomerase affects leaf starch synthesis and floral initiation.Crossref | GoogleScholarGoogle Scholar |

Zuur AF, Ieno EN, Walker NJ, Saveliev AA, Smith GM (2009) ‘Mixed effects models and extensions in ecology with R.’ (Eds M Gail, K Krickeberg, JM Samet, A Tsiatis, W Wong) (Springer: New York, NY, USA)