Environmental Risk Assessment (ERA) is a strategic conceptual framework to characterize the nature and magnitude of risks, to humans and biodiversity, of the release of chemical contaminants in the environment. Several measures have been suggested to enhance the science and application of ERA, including the identification and acknowledgment of uncertainties that potentially influence the outcome of risk assessments, and the appropriate consideration of temporal scale and its linkage to assessment endpoints [1].
Baudrot & Charles [2] proposed to approach these questions by coupling toxicokinetics-toxicodynamics models, which describe the time-course of processes leading to the adverse effects of a toxicant, with Bayesian inference. TKTD models separate processes influencing an organismal internal exposure (´toxicokinetics´, i.e., the uptake, bioaccumulation, distribution, biotransformation and elimination of a toxicant) from processes leading to adverse effects and ultimately its death (´toxicodynamics´) [3]. Although species and substance specific, the mechanistic nature of TKTD models facilitates the comparison of different toxicants, species, life stages, environmental conditions and endpoints [4].
Baudrot & Charles [2] investigated the use of a Bayesian framework to assess the uncertainties surrounding the calibration of General Unified Threshold Models of Survival (a category of TKTD) with data from standard toxicity tests, and their propagation to predictions of regulatory toxicity endpoints such as LC(x,t) [the lethal concentration affecting any x% of the population at any given exposure duration of time t] and MF(x,t) [an exposure multiplication factor leading to any x% effect reduction due to the contaminant at any time t].
Once calibrated with empirical data, GUTS models were used to explore individual survival over time, and under untested exposure conditions. Lethal concentrations displayed a strong curvilinear decline with time of exposure. For a given total amount of contaminant, pulses separated by short time intervals yielded higher mortality than pulses separated by long time intervals, as did few pulses of high amplitude when compared to multiple pulses of low amplitude. The response to a pulsed contaminant exposure was strongly influenced by contaminant depuration times. These findings highlight one important contribution of TKTD modelling in ecotoxicology: they represent just a few of the hundreds of exposure scenarios that could be mathematically explored, and that would be unfeasible or even unethical to conduct experimentally.
GUTS models were also used for interpolations or extrapolations of assessment endpoints, and their marginal distributions. A case in point is the incipient lethal concentration. The responses of model organisms to contaminants in standard toxicity tests are typically assessed at fixed times of exposure (e.g. 24h or 48h in the Daphnia magna acute toxicity test). However, because lethal concentrations are strongly time-dependent, it has been suggested that a more meaningful endpoint would be the incipient (i.e. asymptotic) lethal concentration when time of exposure increases to infinity. The authors present a mathematical solution for calculating the marginal distribution of such incipient lethal concentration, thereby providing both more relevant information and a way of comparing experiments, compounds or species tested for different periods of time.
Uncertainties were found to change drastically with time of exposure, being maximal at extreme values of x for both LC(x,t) and MF(x,t). In practice this means that assessment endpoints estimated when the effects of the contaminant are weak (such as LC10, the contaminant concentration resulting in the mortality of 10% of the experimental population), a commonly used assessment value in ERA, are prone to be highly variable.
The authors end with recommendations for improved experimental design, including (i) using assessment endpoints at intermediate values of x (e.g., LC50 instead of LC10) (ii) prolonging exposure and recording mortality over the course of the experiment (iii) experimenting one or few peaks of high amplitude close to each other when assessing pulsed exposure. Whereas these recommendations are not that different from current practices, they are based on a more coherent mechanistic grounding.
Overall, this and other contributions from Charles, Baudrot and their research group contribute to turn TKTD models into a real tool for Environmental Risk Assessment. Further enhancement of ERA´s science and application could be achieved by extending the use of TKTD models to sublethal rather than lethal effects, and to chronic rather than acute exposure, as these are more controversial issues in decision-making regarding contaminated sites.

References

[1] Dale, V. H., Biddinger, G. R., Newman, M. C., Oris, J. T., Suter, G. W., Thompson, T., ... & Chapman, P. M. (2008). Enhancing the ecological risk assessment process. Integrated environmental assessment and management, 4(3), 306-313. doi: 10.1897/IEAM_2007-066.1
[2] Baudrot, V., & Charles, S. (2018). Recommendations to address uncertainties in environmental risk assessment using toxicokinetics-toxicodynamics models. bioRxiv, 356469, ver. 3 peer-reviewed and recommended by PCI Ecol. doi: 10.1101/356469
[3] EFSA Panel on Plant Protection Products and their Residues (PPR), Ockleford, C., Adriaanse, P., Berny, P., Brock, T., Duquesne, S., ... & Kuhl, T. (2018). Scientific Opinion on the state of the art of Toxicokinetic/Toxicodynamic (TKTD) effect models for regulatory risk assessment of pesticides for aquatic organisms. EFSA Journal, 16(8), e05377. doi: 10.2903/j.efsa.2018.5377
[4] Jager, T., Albert, C., Preuss, T. G., & Ashauer, R. (2011). General unified threshold model of survival-a toxicokinetic-toxicodynamic framework for ecotoxicology. Environmental science & technology, 45(7), 2529-2540. doi: 10.1021/es103092a

Functional traits are “morpho-physio-phenological traits which impact fitness indirectly via their effects on growth, reproduction and survival” [1]. Most functional traits are defined at organ level, e.g. for leaves, roots and stems, and reflect key aspects of resource acquisition and resource use by organisms for their development and reproduction [2]. More rarely, some functional traits can be related to spatial development, such as vegetative height and lateral spread in plants.
Organ-level traits are especially popular because they can be measured in a standard way and easily compared over many plants. But these traits can broadly vary during the life of an organism. For instance, Roggy et al. [3] found that Leaf Mass Area can vary from 30 to 140 g.m^(-2) between seedling and adult stages for the canopy tree Dicorynia guianensis in French Guiana. Fortunel et al. [4] have also showed that developmental stages much contribute to functional trait variation within several Micropholis tree species in lowland Amazonia.
The way plants grow and invest resources into organs is variable during life and allows defining specific developmental sequences and architectural models [5,6]. There is clear ontogenic variation in leaf number, leaf properties and ramification patterns. Ontogenic variations reflect changing adaptation of an individual over its life, depending on the changing environmental conditions.
In this regard, measuring a single functional trait at organ level in adult trees should miss the variation of resource acquisition and use strategies over time. Thus we should built a more integrative approach of ecological development, also called “eco-devo” approach [7].
Although the ecological significance of ontogeny and developmental strategies is now well known, the extent to which it contributes to explain species survival and coexistence in communities is still broadly ignored in functional ecology. Levionnois et al. [8] investigated intraspecific variation of functional traits and growth trajectories in a typical, early-successional tree species in French Guiana, Amazonia. This species, Cecropia obtusa, is generalist regarding soil type and can be found on both white sand and ferralitic soil. The study examines whether there in intraspecific variation in functional traits and growth trajectories of C. obtusa in response to the contrasted soil types.
The tree communities observed on the two types of soils include species with distinctive functional trait values, that is, there are changes in species composition related to different species strategies along the classical wood and leaf economic spectra. The populations of C. obtusa found on the two soils showed some difference in functional traits, but it did not concern traits related to the main economic spectra. Conversely, the populations showed different growth strategies, in terms of spatial and temporal development.
The major lessons we can learn from the study are:
(i) Functional traits measured at organ level cannot reflect well how long-lived plants collect and invest resources during their life. The results show the potential of considering architectural and developmental traits together with organ-level functional traits, to better acknowledge the variation in ecological strategies over plant life, and thus to better understand community assembly processes.
(ii) What makes functional changes between communities differs when considering interspecific and intraspecific variation. Species turnover should encompass different corteges of soil specialists. These specialists are sorted along economic spectra, as shown in tropical rainforests and globally [2]. Conversely, a generalist species such as C. obtusa does occur on contrasted soil, which entails that it can accommodate the contrasted ecological conditions. However, the phenotypic adjustment is not related to how leaves and wood ensure photosynthesis, water and nutrient acquisition, but regards the way the resources are allocated to growth and reproduction over time.
The results of the study stress the need to better integrate growth strategies and ontogeny in the research agenda of functional ecology. We can anticipate that organ-level functional traits and growth trajectories will be more often considered together in ecological studies. The integration should help better understand the temporal niche of organisms, and how organisms can coexist in space and time with other organisms during their life. Recently, Klimešová et al. [9] have proposed standardized protocols for collecting plant modularity traits. Such effort to propose easy-to-measure traits representing plant development and ontogeny, with clear functional roles, should foster the awaited development of an “eco-devo” approach.

References

[1] Violle, C., Navas, M. L., Vile, D., Kazakou, E., Fortunel, C., Hummel, I., & Garnier, E. (2007). Let the concept of trait be functional!. Oikos, 116(5), 882-892. doi: 10.1111/j.0030-1299.2007.15559.x
[2] Díaz, S. et al. (2016). The global spectrum of plant form and function. Nature, 529(7585), 167-171. doi: 10.1038/nature16489
[3] Roggy, J. C., Nicolini, E., Imbert, P., Caraglio, Y., Bosc, A., & Heuret, P. (2005). Links between tree structure and functional leaf traits in the tropical forest tree Dicorynia guianensis Amshoff (Caesalpiniaceae). Annals of forest science, 62(6), 553-564. doi: 10.1051/forest:2005048
[4] Fortunel, C., Stahl, C., Heuret, P., Nicolini, E. & Baraloto, C. (2020). Disentangling the effects of environment and ontogeny on tree functional dimensions for congeneric species in tropical forests. New Phytologist. doi: 10.1111/nph.16393
[5] Barthélémy, D., & Caraglio, Y. (2007). Plant architecture: a dynamic, multilevel and comprehensive approach to plant form, structure and ontogeny. Annals of botany, 99(3), 375-407. doi: 10.1093/aob/mcl260
[6] Hallé, F., & Oldeman, R. A. (1975). An essay on the architecture and dynamics of growth of tropical trees. Kuala Lumpur: Penerbit Universiti Malaya.
[7] Sultan, S. E. (2007). Development in context: the timely emergence of eco-devo. Trends in Ecology & Evolution, 22(11), 575-582. doi: 10.1016/j.tree.2007.06.014
[8] Levionnois, S., Tysklind, N., Nicolini, E., Ferry, B., Troispoux, V., Le Moguedec, G., Morel, H., Stahl, C., Coste, S., Caron, H. & Heuret, P. (2020). Soil variation response is mediated by growth trajectories rather than functional traits in a widespread pioneer Neotropical tree. bioRxiv, 351197, ver. 4 peer-reviewed and recommended by PCI Ecology. doi: 10.1101/351197
[9] Klimešová, J. et al. (2019). Handbook of standardized protocols for collecting plant modularity traits. Perspectives in Plant Ecology, Evolution and Systematics, 40, 125485. doi: 10.1016/j.ppees.2019.125485

Urban areas are expanding worldwide, and have become a dominant part of the landscape for many species. Urbanization can fragment pre-existing populations of vulnerable species leading to population declines and the loss of connectivity. On the other hand, expansion of urban areas can also facilitate the spread of human commensals including pests. Knowledge of the features of cityscapes that facilitate gene flow and maintain diversity of pests is thus key to their management and eradication.
Cities are complex mosaics of natural and manmade surfaces, and habitat quality is not only influenced by physical aspects of the cityscape but also by socioeconomic factors and human behaviour. Constant development means that cities also change rapidly in time; contemporary urban life reflects only a snapshot of the environmental conditions faced by populations. It thus remains a challenge to identify the features that actually drive ecology and evolution of populations in cities [1]. While several studies have highlighted strong urban clines in genetic structure and adaption [2], few have considered the influence of factors beyond physical aspects of the cityscape or historical processes.
In this paper, Stragier et al. [3] sought to identify the current and past features of the cityscape and socioeconomic factors that shape genetic structure and diversity of the house mouse (Mus musculus domesticus) in Dakar, Senegal. The authors painstakingly digitized historical maps of Dakar from the time of European settlement in 1862 to present. The authors found that the main spatial genetic cline was best explained by historical cityscape features, with higher apparent gene flow and genetic diversity in areas that were connected earlier to initial European settlements. Beyond the main trend of spatial genetic structure, they found further evidence that current features of the cityscape were important. Specifically, areas with low vegetation and poor housing conditions were found to support large, genetically diverse populations. The authors demonstrate that their results are reproducible using several statistical approaches, including modeling that explicitly accounts for spatial autocorrelation.
The work of Stragier et al. [3] thus highlights that populations of city-dwelling species are the product of both past and present cityscapes. Going forward, urban evolutionary ecologists should consider that despite the potential for rapid evolution in urban landscapes, the signal of a species’ colonization can remain for generations.

References

[1] Rivkin, L. R., Santangelo, J. S., Alberti, M. et al. (2019). A roadmap for urban evolutionary ecology. Evolutionary Applications, 12(3), 384-398. doi: 10.1111/eva.12734
[2] Miles, L. S., Rivkin, L. R., Johnson, M. T., Munshi‐South, J. and Verrelli, B. C. (2019). Gene flow and genetic drift in urban environments. Molecular ecology, 28(18), 4138-4151. doi: 10.1111/mec.15221
[3] Stragier, C., Piry, S., Loiseau, A., Kane, M., Sow, A., Niang, Y., Diallo, M., Ndiaye, A., Gauthier, P., Borderon, M., Granjon, L., Brouat, C. and Berthier, K. (2020). Interplay between historical and current features of the cityscape in shaping the genetic structure of the house mouse (Mus musculus domesticus) in Dakar (Senegal, West Africa). bioRxiv, 557066, ver. 4 peer-reviewed and recommended by PCI Ecology. doi: 10.1101/557066

Are biological invasions driven by a few pioneers, running ahead of their conspecifics? Or are these pioneers constantly being caught up by, and folded into, the larger flux of propagules from the established populations behind them?
In ecology and beyond, these two scenarios are known as "pulled" and "pushed" fronts, and they come with different expectations. In a pushed front, invasion speed is not just a matter of how good individuals are at dispersing and settling new locations. It becomes a collective, density-dependent property of population fluxes. And in particular, it can depend on the equilibrium abundance of the established populations inside the range, i.e. the species’ carrying capacity K, factoring in its abiotic environment and biotic interactions.
This realization is especially important because it can flip around our expectations about which species expand fast, and how to manage them. We tend to think of initial colonization and long-term abundance as two independent axes of variation among species or indeed as two ends of a spectrum, in the classic competition-colonization tradeoff [1]. When both play into invasion speed, good dispersers might not outrun good competitors. This is useful knowledge, whether we want to contain an invasion or secure a reintroduction.
In their study "When higher carrying capacities lead to faster propagation", Haond et al [2] combine mathematical analysis, Individual-Based simulations and experiments to show that various mechanisms can cause pushed fronts, whose speed increases with the carrying capacity K of the species. Rather than focus on one particular angle, the authors endeavor to demonstrate that this qualitative effect appears again and again in a variety of settings.
It is perhaps surprising that this notable and general connection between K and invasion speed has managed to garner so little fame in ecology. A large fraction of the literature employs the venerable Fisher-KPP reaction-diffusion model, which combines local logistic growth with linear diffusion in space. This model has prompted both considerable mathematical developments [3] and many applications to modelling real invasions [4]. But it only allows pulled fronts, driven by the small populations at the edge of a species range, with a speed that depends only on their initial growth rate r.
This classic setup is, however, singular in many ways. Haond et al [2] use it as a null model, and introduce three mechanisms or factors that each ensure a role of K in invasion speed, while giving less importance to the pioneers at the border.
Two factors, the Allee effect and demographic stochasticity, make small edge populations slower to grow or less likely to survive. These two factors are studied theoretically, and to make their claims stronger, the authors stack the deck against K. When generalizing equations or simulations beyond the null case, it is easy to obtain functional forms where the parameter K does not only play the role of equilibrium carrying capacity, but also affects dynamical properties such as the maximum or mean growth rate. In that case, it can trivially change the propagation speed, without it meaning anything about the role of established populations behind the front. Haond et al [2] avoid this pitfall by disentangling these effects, at the cost of slightly more peculiar expressions, and show that varying essentially nothing but the carrying capacity can still impact the speed of the invasion front.
The third factor, density-dependent dispersal, makes small populations less prone to disperse. It is well established empirically and theoretically that various biological mechanisms, from collective organization to behavioral switches, can prompt organisms in denser populations to disperse more, e.g. in such a way as to escape competition [5]. The authors demonstrate how this effect induces a link between carrying capacity and invasion speed, both theoretically and in a dispersal experiment on the parasitoid wasp, Trichogramma chilonis.
Overall, this study carries a simple and clear message, supported by valuable contributions from different angles. Although some sections are clearly written for the theoretical ecology crowd, this article has something for everyone, from the stray physicist to the open-minded manager. The collaboration between theoreticians and experimentalists, while not central, is worthy of note. Because the narrative of this study is the variety of mechanisms that can lead to the same qualitative effect, the inclusion of various approaches is not a gimmick, but helps drive home its main message. The work is fairly self-contained, although one could always wish for further developments, especially in the direction of more quantitative testing of these mechanisms.
In conclusion, Haond et al [2] effectively convey the widely relevant message that, for some species, invading is not just about the destination, it is about the many offspring one makes along the way.

References

[1] Levins, R., & Culver, D. (1971). Regional Coexistence of Species and Competition between Rare Species. Proceedings of the National Academy of Sciences, 68(6), 1246–1248. doi: 10.1073/pnas.68.6.1246
[2] Haond, M., Morel-Journel, T., Lombaert, E., Vercken, E., Mailleret, L., & Roques, L. (2018). When higher carrying capacities lead to faster propagation. BioRxiv, 307322. doi: 10.1101/307322
[3] Crooks, E. C. M., Dancer, E. N., Hilhorst, D., Mimura, M., & Ninomiya, H. (2004). Spatial segregation limit of a competition-diffusion system with Dirichlet boundary conditions. Nonlinear Analysis: Real World Applications, 5(4), 645–665. doi: 10.1016/j.nonrwa.2004.01.004
[4] Shigesada, N., & Kawasaki, K. (1997). Biological Invasions: Theory and Practice. Oxford University Press, UK.
[5] Matthysen, E. (2005). Density-dependent dispersal in birds and mammals. Ecography, 28(3), 403–416. doi: 10.1111/j.0906-7590.2005.04073.x

Invasive plant species are widely studied by the ecologist community, especially in wetlands. Indeed, alien plants are considered one of the major threats to wetland biodiversity (Reid et al., 2019). Ludwigia grandiflora subsp. hexapetala (Hook. & Arn.) G.L.Nesom & Kartesz, 2000 (Lgh) is one of them and has received particular attention for a long time (Hieda et al., 2020; Thouvenot, Haury, & Thiebaut, 2013). The ecology of this invasive species and its effect on its biotic and abiotic environment has been studied in previous works. Different processes were demonstrated to explain their invasibility such as allelopathic interference (Dandelot et al., 2008), resource competition (Gérard et al., 2014), and high phenotypic plasticity (Thouvenot, Haury, & Thiébaut, 2013), to cite a few of them. However, although vegetative reproduction is a well-known invasive process for alien plants like Lgh (Glover et al., 2015), the sexual reproduction of this species is still unclear and may help to understand the Lgh population dynamics.

Portillo Lemus et al. (2021) showed that two floral morphs of Lgh co-exist in natura, involving self-compatibility for short-styled phenotype and self-incompatibility for long-styled phenotype processes. This new article (Portillo Lemus et al., 2022) goes further and details the underlying mechanisms of the sexual reproduction of the two floral morphs.

Complementing their previous study, the authors have described a late self-incompatible process associated with the long-styled morph, which authorized a small proportion of autogamy. Although this represents a small fraction of the L-morph reproduction, it may have a considerable impact on the L-morph population dynamics. Indeed, authors report that “floral morphs are mostly found in allopatric monomorphic populations (i.e., exclusively S-morph or exclusively L-morph populations)” with a large proportion of L-morph populations compared to S-morph populations in the field. It may seem counterintuitive as L-morph mainly relies on cross-fecundation. 

Results show that L-morph autogamy mainly occurs in the fall, late in the reproduction season. Therefore, the reproduction may be ensured if no exogenous pollen reaches the stigma of L-morph individuals. It partly explains the large proportion of L-morph populations in the field. 

Beyond the description of late-acting self-incompatibility, which makes the Onagraceae a third family of Myrtales with this reproductive adaptation, the study raises several ecological questions linked to the results presented in the article. First, it seems that even if autogamy is possible, Lgh would favour allogamy, even in S-morph, through the faster development of pollen tubes from other individuals. This may confer an adaptative and evolutive advantage for the Lgh, increasing its invasive potential. The article shows this faster pollen tube development in S-morph but does not test the evolutive consequences. It is an interesting perspective for future research. It would also be interesting to describe cellular processes which recognize and then influence the speed of the pollen tube. Second, the importance of sexual reproduction vs vegetative reproduction would also provide information on the benefits of sexual dimorphism within populations. For instance, how fruit production increases the dispersal potential of Lgh would help to understand Lgh population dynamics and to propose adapted management practices (Delbart et al., 2013; Meisler, 2009).

To conclude, the study proposes a morphological, reproductive and physiological description of the Lgh sexual reproduction process. However, underlying ecological questions are well included in the article and the ecophysiological results enlighten some questions about the role of sexual reproduction in the invasiveness of Lgh. I advise the reader to pay attention to the reviewers’ comments; the debates were very constructive and, thanks to the great collaboration with the authorship, lead to an interesting paper about Lgh reproduction and with promising perspectives in ecology and invasion ecology.

References

Dandelot S, Robles C, Pech N, Cazaubon A, Verlaque R (2008) Allelopathic potential of two invasive alien Ludwigia spp. Aquatic Botany, 88, 311–316. https://doi.org/10.1016/j.aquabot.2007.12.004

Delbart E, Mahy G, Monty A (2013) Efficacité des méthodes de lutte contre le développement de cinq espèces de plantes invasives amphibies : Crassula helmsii, Hydrocotyle ranunculoides, Ludwigia grandiflora, Ludwigia peploides et Myriophyllum aquaticum (synthèse bibliographique). BASE, 17, 87–102. https://popups.uliege.be/1780-4507/index.php?id=9586

Gérard J, Brion N, Triest L (2014) Effect of water column phosphorus reduction on competitive outcome and traits of Ludwigia grandiflora and L. peploides, invasive species in Europe. Aquatic Invasions, 9, 157–166. https://doi.org/10.3391/ai.2014.9.2.04

Glover R, Drenovsky RE, Futrell CJ, Grewell BJ (2015) Clonal integration in Ludwigia hexapetala under different light regimes. Aquatic Botany, 122, 40–46. https://doi.org/10.1016/j.aquabot.2015.01.004

Hieda S, Kaneko Y, Nakagawa M, Noma N (2020) Ludwigia grandiflora (Michx.) Greuter & Burdet subsp. hexapetala (Hook. & Arn.) G. L. Nesom & Kartesz, an Invasive Aquatic Plant in Lake Biwa, the Largest Lake in Japan. Acta Phytotaxonomica et Geobotanica, 71, 65–71. https://doi.org/10.18942/apg.201911

Meisler J (2009) Controlling Ludwigia hexaplata in Northern California. Wetland Science and Practice, 26, 15–19. https://doi.org/10.1672/055.026.0404

Portillo Lemus LO, Harang M, Bozec M, Haury J, Stoeckel S, Barloy D (2022) Late-acting self-incompatible system, preferential allogamy and delayed selfing in the heteromorphic invasive populations of Ludwigia grandiflora subsp. hexapetala. bioRxiv, 2021.07.15.452457, ver. 4 peer-reviewed and recommended by Peer Community in Ecology. https://doi.org/10.1101/2021.07.15.452457

Portillo Lemus LO, Bozec M, Harang M, Coudreuse J, Haury J, Stoeckel S, Barloy D (2021) Self-incompatibility limits sexual reproduction rather than environmental conditions in an invasive water primrose. Plant-Environment Interactions, 2, 74–86. https://doi.org/10.1002/pei3.10042

Reid AJ, Carlson AK, Creed IF, Eliason EJ, Gell PA, Johnson PTJ, Kidd KA, MacCormack TJ, Olden JD, Ormerod SJ, Smol JP, Taylor WW, Tockner K, Vermaire JC, Dudgeon D, Cooke SJ (2019) Emerging threats and persistent conservation challenges for freshwater biodiversity. Biological Reviews, 94, 849–873. https://doi.org/10.1111/brv.12480

Thouvenot L, Haury J, Thiebaut G (2013) A success story: water primroses, aquatic plant pests. Aquatic Conservation: Marine and Freshwater Ecosystems, 23, 790–803. https://doi.org/10.1002/aqc.2387

Thouvenot L, Haury J, Thiébaut G (2013) Seasonal plasticity of Ludwigia grandiflora under light and water depth gradients: An outdoor mesocosm experiment. Flora - Morphology, Distribution, Functional Ecology of Plants, 208, 430–437. https://doi.org/10.1016/j.flora.2013.07.004