The age of tomato plants affects the development of Macrosiphum euphorbiae (Thomas, 1878) (Hemiptera) colonies Received for publication: 24 July, 2020. Accepted for publication: 25 February, 2021. Doi: 10.15446/agron.colomb.v39n1.89301 1 Scuola di Scienze Agrarie, Forestali, Alimentari ed Ambientali, Università degli Studi della Basilicata, Potenza (Italy). 2 Istituto Comprensivo “Via Pirotta”, Rome (Italy). 3 Dipartimento di Scienze, Università degli Studi della Basilicata, Potenza (Italy). 4 Jardín Botánico de Bogotá “José Celestino Mutis”, Bogota (Colombia). * Corresponding author: jduran@jbb.gov.co Agronomía Colombiana 39(1), 108-112, 2021 ABSTRACT RESUMEN We tested the hypothesis that the intensity and duration of Macrosiphum euphorbiae infestations in tomato depend on both the age (phenological stage) of the host plant and the initial number of aphids present in the colony. We compared the effects of three initial levels of infestation and two pheno- logical stages of the plant (pre-f lowering and f lowering stages) on infestation curves. The position of the infestation peak over time was significantly affected by the plant phenological phase. Populations of M. euphorbiae reached the highest peak of abun- dance on plants infested at the pre-f lowering stage compared to those subsequently infested. Within a phenological phase, the maximum abundance also varied according to the initial aphid density on the plant. The implications concerning the management of the pest in the field are brief ly discussed. Se planteó la hipótesis que la intensidad y duración de las in- festaciones de Macrosiphum euphorbiae en el tomate dependen simultáneamente de la edad (etapa fenológica) de la planta hospedera y del número inicial de áfidos de la colonia. Se compararon los efectos de tres niveles iniciales de infestación y dos fases fenológicas de la planta (fases de pre-f loración y f loración) sobre las curvas de infestación. La posición del pico de infestación a lo largo del tiempo estuvo afectada significa- tivamente por la fase fenológica de la planta. Las poblaciones de M. euphorbiae alcanzaron su máximo pico de abundancia en las plantas infestadas en la fase de pre-f loración compa- radas con las plantas infestadas en fases posteriores. Dentro de la fase fenológica, la máxima abundancia varió también de acuerdo con la densidad inicial de áfidos en la planta. Las implicaciones concernientes al manejo de la plaga en campo se discuten brevemente. Key words: aphids, curve of infestation, plant phenological stages, Solanum lycopersicum, trophic interactions. Palabras clave: áfidos, curva de infestación, estados fenológicos de la planta, Solanum lycopersicum, interacciones tróficas. The age of tomato plants affects the development of Macrosiphum euphorbiae (Thomas, 1878) (Hemiptera) colonies La edad de las plantas de tomate afecta el desarrollo de las colonias de Macrosiphum euphorbiae (Thomas, 1878) (Hemiptera) Vincenzo Trotta1, Irene Toma2, Pierluigi Forlano1, Paolo Fanti3, Juliana Durán Prieto4*, and Donatella Battaglia1 Introduction Aphids are r-strategist (Gadgil & Solbrig, 1972) insects characterized by their capacity for extremely rapid popula- tion growth together with a transient relationship with the host plant (Powell et al., 2006). Aphid colonies typically decline after an initial period of rapid growth. This last phase is largely due to the production of winged morphs and is stimulated by crowding, the presence of natural enemies, and the decrease in plant quality (Müller et al., 2001; Irwin et al., 2007). The curve of infestation for the same aphid species can be dramatically affected by the host plant despite keeping all the other factors equal (Larocca et al., 2011). Macrosiphum euphorbiae (Thomas, 1878) (Hemiptera: Aphididae) is an important aphid pest that causes the most significant direct damage among the aphid species that attack tomato (Solanum lycopersicum L. 1753) (Perring et al., 2018). In a previous study concerning the development of M. euphorbiae colonies on tomato plants, we observed that cultivars and water stress affect the peak but do not interfere with the length of the infestation (Rivelli et al., 2013). This result was confirmed by field observations (Colella et al., 2014). Several models of aphid population dynamics, which exclude predator inf licted mortality as a regulating fac- tor, identify the initial number of aphids as the main SCIENTIFIC NOTE https://doi.org/10.15446/agron.colomb.v39n1.89301 109Trotta, Toma, Forlano, Fanti, Durán Prieto, and Battaglia: The age of tomato plants affects the development of Macrosiphum euphorbiae (Thomas, 1878) (Hemiptera) colonies determinant of the width of the density curve and/or size and the timing of the peak of maximum density (Kindl- mann et al., 2007). In addition, some of these models take into account variations in the “aphid-carrying capacity” due to changes in host plant quality over time. Plant qual- ity may largely depend on its phenology, which modifies the physiological priorities for resource allocation (Boege & Marquis, 2005). In this study, we hypothesized that the initial number of aphids and/or the age of the plant are significant factors that can inf luence the duration and, consequently, the harmfulness of M. euphorbiae infestations in tomato. Materials and methods Tomato plants of the cultivar Rio Grande (pear-shaped processing tomato for paste and concentrated juice) were used for both M. euphorbiae rearing and the experiment. We bought tomato plants in polystyrene trays from a nurs- ery that used seeds produced by OLTER (Piacenza, Italy). Seedlings were transplanted and grown individually in plastic pots (18 cm diameter, 19 cm height) with commercial soil (COMPO SANA® universal potting soil, International Kingenta Group, Italy). Macrosiphum euphorbiae individuals were originally col- lected in a tomato crop in Scafati (Salerno, Italy) and reared on tomato plants at 20 ± 1°C, 65 ± 5% relative humidity, and a photoperiod of 18 h light:6 h dark. The experiment was carried out at the University of Ba- silicata, Italy (40°36’ N, 15°48’ E), in the summertime, in a naturally lit and temperature-controlled greenhouse maintained at 20°C (with an oscillation between 15°C at night and 28°C at 12-14 pm). We compared three initial levels of infestation and two phenological phases of the plant. Nine plants were infested 7 d after transplant (pre- f lowering stage; 18.0 ± 2.6 cm height; root biomass 0.10 ± 0.04 g; leaf biomass 0.43 ± 0.08 g) and 9 plants were infested 26 d after transplant (beginning of f lowering; 39.5 ± 4.6 cm height; root biomass 1.5 ± 0.3 g; leaf biomass 6.4 ± 3.0 g). Three levels of initial infestation were obtained by placing 6, 10, and 15 M. euphorbiae adults (three replicates for each level of initial infestation: R1, R2, and R3), respectively, on plants. We selected this range of initial infestation to simulate a variable but incipient attack of the pest. We used apterous females for simplicity, even if it is more plausible that the infestation starts with winged females under natural conditions. All plants were infested on the same day. The plants in the f lower transition stage had been transplanted 19 d before the plants in pre-f lowering. Plants were placed in the greenhouse in a completely randomized experimental design and regularly checked. The number of aphids in the whole plant was counted using a magnifying glass. The infestation curves for each plant stage and each level of initial infestation were adjusted to second degree polynomi- als passing through the origin according to a “cumulative density model” (Kindlmann et al., 2007). In the present study, we used the polynomial parameters to calculate the axes of the vertex of the second-degree polynomials: the theoretical values of maximum abundance of aphids (the Y axis) and of the day needed to reach the maximum abundance (the X axis). The theoretical values of the maxi- mum number of aphids and of the days needed to reach it were then analyzed with factorial ANOVAs, with “stage” (7 and 26 d after transplant, i.e., pre-f lowering and f lowering stages) and the initial level of infestation as fixed effects. The same analyses were also performed on the observed values of the maximum number of aphids and of the days needed to reach it. All the analyses in this study were per- formed using the R 3.2.3 software (R Core Team, 2014). Results and discussion The response curves for the three replicates of the pre- f lowering and f lowering stages are shown separately for each level of initial infestation in Figure 1. In all cases, values of adjusted R2 are significant and the adjustments to second degree polynomial provide plausible values for the vertex. Positions of the theoretical and observed vertices of all the curves (mean values of three replicates) with the respective standard errors are shown in Figure 2. The mean value of the peak and time when the peak occurred for each plant stage are shown in Tables 1 (theoretical values) and 2 (observed values). Aphid population reached a much higher peak on plants infested after 7 d from transplant (pre-f lowering) com- pared to those subsequently infested (f lowering) (theo- retical values: F1,12 = 48.18, P<0.001; observed values: F1,12 = 83.28, P<0.001). Significant differences in the maximum abundance were also detected for the level of infestation (theoretical values: F2,12 = 4.45, P<0.05; observed values: F2,12 = 5.49, P<0.05) but not for the interaction between the phenological stage and the level of infestation. These differences are mainly due to the peak reached on plants initially infested with 15 aphids on the 7th d (Figs. 1C and 2). The vertices of the other curves are all very close. 110 Agron. Colomb. 39(1) 2021 7065605550454035302520151050 1200 1000 800 600 400 200 0 1400 B - Initial infestation: 10 aphids Days after initial infestation N um be r of a ph id s pe r pl an t 65 70605550454035302520151050 1200 1000 800 600 400 200 0 1400 C - Initial infestation: 15 aphids Days after initial infestation Theor. curve Pre-flowering R3 Theor. curve Pre-flowering R2 Theor. curve Pre-flowering R1 Theor. curve Flowering R3 Theor. curve Flowering R2 Theor. curve Flowering R1 N um be r of a ph id s pe r pl an t 7065 75605550454035302520151050 1200 1000 800 600 400 200 0 1400 A - Initial infestation: 6 aphids Days after initial infestation N um be r of a ph id s pe r pl an t 454035302520151050 1200 1000 800 600 400 200 0 1400 Days after initial infestation Observed flowering 6 aphids Observed flowering 15 aphids Observed flowering 10 aphids Observed pre-flowering 6 aphids Observed pre-flowering 15 aphids Observed pre-flowering 10 aphids Estimated flowering 6 aphids Estimated flowering 15 aphids Estimated flowering 10 aphids Estimated pre-flowering 6 aphids N um be r of a ph id s pe r pl an t Estimated pre-flowering 15 aphids Estimated pre-flowering 10 aphids FIGURE 1. Curves of Macrosiphum euphorbiae infestation estimated for plants initially infested at the pre-flowering stage or at the beginning of flowering. A) Initial infestation: 6 aphids; B) initial infestation: 10 aphids; C) initial infestation: 15 aphids. Continuous lines connect the experi- mental points; dotted lines (theor.) represent the theoretical adjustments to second degree polynomial; R1, R2 and R3 are the three replicates for each level of initial infestation and stage. FIGURE 2. Maximum abundance of Macrosiphum euphorbiae colonies according to the initial infestation level (6, 10, or 15 aphids per plant) and the stage of the attacked plant (pre-flowering and flowering stages). TABLE 1. Theoretical maximum abundance and time of maximum abun- dance according to the plant stage at the beginning of infestation and the initial number of aphids (mean ± standard error). Plant stage at the beginning of infestation Initial number of aphids Time of maximum abundance (d) Maximum abundance Pre-flowering 6 35.05 ± 2.4 364.3 ± 120 Pre-flowering 10 35.63 ± 0.5 357.7 ± 36 Pre-flowering 15 38.16 ± 6.3 689.4 ± 102 Flowering 6 15.93 ± 2.4 71.7 ± 25 Flowering 10 18.65 ± 1.7 73.4 ± 26 Flowering 15 17.51 ± 1.4 99.0 ± 30 TABLE 2. Observed maximum abundance and time of maximum abun- dance according to the plant stage at the beginning of the infestation and the initial number of aphids (mean ± standard error). Plant stage at the beginning of infestation Initial number of aphids Time of maximum abundance (d) Maximum abundance Pre-flowering 6 34.67 ± 1.3 745.3 ± 209 Pre-flowering 10 36.00 ± 0.0 617.7 ± 95 Pre-flowering 15 33.33 ± 1.3 1186.3 ± 26 Flowering 6 15.00 ± 3.0 94.3 ± 35 Flowering 10 13.67 ± 2.6 116.0 ± 39 Flowering 15 11.33 ± 1.5 154.0 ± 35 111Trotta, Toma, Forlano, Fanti, Durán Prieto, and Battaglia: The age of tomato plants affects the development of Macrosiphum euphorbiae (Thomas, 1878) (Hemiptera) colonies The position of the peaks of infestation over time (the day after infestation when the aphid population reached its maximum abundance) was significantly affected only by the plant stage at the start of infestation (theoretical values: F1,12 = 57.35, P<0.001; observed values: F1,12 = 191, P<0.001). Our results show that the duration of M. euphorbiae in- festation in tomato mainly depends on the phenological stage of the host plant (it is higher and persists longer on early infested plants), while the degree of infestation seems to be irrelevant considering the initial levels used in this experiment. Plant defense significantly changes with plant development from the seedling to juvenile to mature and senescent stages (Barton & Boege, 2017). Ontogenetic changes in plant de- fenses observed in many species do not only concern the intensity of the defensive response but also the mechanisms involved. This can be partially explained by the fact that allocation costs and resource allocation priorities vary as plants develop (Boege & Marquis, 2005). Two important changes that can impact the development of insect infestations occur during the f lower transition: the negative regulation of herbivory-induced jasmonic acid (JA) signaling (Gaquerel & Stitz, 2017) and the increase of the C:N ratio in the phloem sap (Corbesier et al., 2002). The repression of the JA signaling-based induction ex- plains why the ability of tomato plants to induce defenses against Manduca sexta (Linnaeus, 1763) is lost when the reproductive stage is reached (Wolfson & Murdock,1990). The down-regulation of the JA signaling pathway mainly favors chewing insects. On the other hand, the increase of the C:N ratio in the phloem sap is unfavorable for aphids since amino acid availability significantly affects their growth and reproduction (Ponder et al., 2000; Karley et al., 2002). In the case of aphids, the nutritional quality of the host plant can have an even higher impact than that of the induced defenses (Battaglia et al., 2013). In addition to the nutritional aspects, we must consider that plants have evolved multiple defense strategies against aphids, including constitutive as well as inducible factors (Nalam et al., 2019) that may change during ontogenesis. For example, expanded leaves have greater trichome den- sity and resistance in tomato plants in reproductive stages than in vegetative ones (Mymko & Avila-Sakar, 2019). Furthermore, herbaceous plants usually show a significant increase in secondary chemistry across the entire ontoge- netic trajectory (Barton & Koricheva, 2010). The ontogenetic changes in the nutritional quality of the plant and in the defense strategies may explain the greater development of aphid colonies we observed when we in- fested tomato plants in the pre-f lowering stage compared to the f lowering stage. Our results confirm the better per- formance of aphids on young plants previously reported for other aphid-plant systems, as in the case of Diuraphis noxia (Kurdjumov 1913) on barley (Ma & Bechinski, 2009), and Myzus persicae (Sulzer, 1776) and M. euphorbiae on potato (Karley et al., 2002). The best performance of aphids on young plants has practi- cal consequences on the development of aphid infestations in the field and, therefore, on the control strategies that can be implemented. In fact, the extent of aphid infestations not only depends on the number of winged forms that move from the winter host to the herbaceous crops but also, to a large extent, on the phenological stage in which crop plants are at the time of the aphid invasion. The anticipation of the sowing time, when possible, allows plants to reach a less suitable phenological stage before aphid colonization takes place. This could be a strategy for the control of aphids, as several field studies also suggest (Perring et al., 1988; John et al., 2017; Alam et al., 2020). Acknowledgments We would like to thank the Università degli studi della Basillicata and the Italian Ministry of Education, University and Scientific Research (MIUR) for funding. Author’s contributions DB, PF, and VT were involved in planning and supervised the work; IT, PF, and JD performed the experiments; VT processed the experimental data and performed the analy- sis, and DB drafted the manuscript. All authors discussed the results and commented on the manuscript. Literature cited Alam, J., Nahar, K., Khatun, K., Rashid, H., & Ahmed, K. S. (2020). Impact assessment of different sowing dates on maize aphid, Rhopalosiphum madis infestation in Bangladesh. Sustainability in Food and Agriculture, 1(2), 99–105. https://doi.org/10.26480/ sfna.02.2020.99.105 Barton, K. E., & Boege, K. (2017). Future directions in the ontogeny of plant defence: understanding the evolutionary causes and consequences. Ecology Letters, 20(4), 403–411. https://doi. org/10.1111/ele.12744 Barton, K. E., & Koricheva, J. (2010). The ontogeny of plant defense and herbivory: characterizing general patterns using meta- analysis. The American Naturalist, 175(4), 481–493. https:// doi.org/10.1086/650722 https://doi.org/10.26480/sfna.02.2020.99.105 https://doi.org/10.26480/sfna.02.2020.99.105 https://doi.org/10.1111/ele.12744 https://doi.org/10.1111/ele.12744 https://doi.org/10.1086/650722 https://doi.org/10.1086/650722 112 Agron. Colomb. 39(1) 2021 Battaglia, D., Bossi, S., Cascone, P., Digilio, M. C., Duran Prieto, J., Fanti, P., Guerrieri, E., Iodice, L., Lingua, G., Lorito, M., Maffei, M. E., Massa, N., Ruocco, M., Sasso, R., & Trotta, V. (2013). Tomato below ground-above ground interactions: Trichoderma longibrachiatum affects the performance of Mac- rosiphum euphorbiae and its natural antagonists. Molecular Plant-Microbe Interactions, 26(10), 1249–1256. https://doi. org/10.1094/MPMI-02-13-0059-R Boege, K., & Marquis, R. J. (2005). Facing herbivory as you grow up: the ontogeny of resistance in plants. Trends in Ecol- ogy and Evolution, 20(8), 441–448. https://doi.org/10.1016/j. tree.2005.05.001 Colella, T., Candido, V., Campanelli, G., Camele, I., & Battaglia, D. (2014). Effect of irrigation regimes and artificial mycor- rhization on insect pest infestations and yield in tomato crop. Phytoparasitica, 42, 235–246. https://doi.org/10.1007/ s12600-013-0356-3 Corbesier, L., Bernier, G., & Périlleux, C. (2002). C:N ratio increases in the phloem sap during f loral transition of the long-day plants Sinapis alba and Arabidopsis thaliana. Plant and Cell Physiology, 43(6), 684–688. https://doi.org/10.1093/pcp/pcf071 Gadgil, M., & Solbrig, O. T. (1972). The concept of r- and K-selection: evidence from wild f lowers and some theoretical consider- ations. The American Naturalist, 106(947), 14–31. https://doi. org/10.1086/282748 Gaquerel, E., & Stitz, M. (2017). Insect resistance: an emerging molecular framework linking plant age and JA signaling. Molecular Plant, 10(4), 537–539. https://doi.org/10.1016/j. molp.2017.02.006 Karley, A. J., Douglas, A. E., & Parker, W. E. (2002). Amino acid composition and nutritional quality of potato leaf phloem sap for aphids. Journal of Experimental Biology, 205, 3009–3018. Kindlmann, P., Jarošík, V., & Dixon, A. F. G. (2007). Population dynamics. In H. F. van Emden, & R. Harrington (Eds.), Aphids as crop pests (pp. 311–329). CAB International. https://doi. org/10.1079/9780851998190.0311 Irwin, M. E., Kampmeier, G. E., & Weisser, W. W. (2007). Aphid movement: process and consequences. In H. F. van Emden, & R. Harrington (Eds.), Aphids as crop pests (pp. 153–186). CAB International. https://doi.org/10.1079/9780851998190.0153 John, F., Saeed, N. A., Nadeem, S., & Hamed, M. (2017). Integration of planting time and insecticide to manage aphid infestations in wheat for better crop productivity. Pakistan Journal of Zoology, 49(4), 1343–1351. https://doi.org/10.17582/journal. pjz/2017.49.4.1343.1351 Larocca, A., Fanti, P., Molinaro, A., Mattia, M. F., & Battaglia, D. (2011). Aphid performance on Vicia faba and two southern Italy Phaseolus vulgaris landraces. Bulletin of Insectology, 64(1), 101–106. Ma, Z. S., & Bechinski, E. J. (2009). Life tables and demographic sta- tistics of Russian wheat aphid (Hemiptera: Aphididae) reared at different temperatures and on different host plant growth stages. European Journal of Entomology, 106(2), 205–210. https://doi.org/10.14411/eje.2009.026 Mymko, D., & Avila‐Sakar, G. (2019). The inf luence of leaf onto- genetic stage and plant reproductive phenology on trichome density and constitutive resistance in six tomato variet- ies. Arthropod-Plant Interactions, 13, 797–803. https://doi. org/10.1007/s11829-019-09690-3 Müller, C. B., Williams, I. S., & Hardie, J. (2001). The role of nutrition, crowding and interspecific interactions in the development of winged aphids. Ecological Entomology, 26(3), 330–340. https:// doi.org/10.1046/j.1365-2311.2001.00321.x Nalam, V., Louis, J., & Shah, J. (2019). Plant defense against aphids, the pest extraordinaire. Plant Science, 279, 96–107. https://doi. org/10.1016/j.plantsci.2018.04.027 Perring, T. M., Battaglia, D., Walling, L. L., Toma, I., & Fanti, P. (2018). Aphids: biology, ecology, and management. In W. Wakil, G. E. Brust, & T. M. Perring (Eds.), Sustainable man- agement of arthropod pests of tomato (pp. 15–48). Academic Press. https://doi.org/10.1016/B978-0-12-802441-6.00002-4 Perring, T. M., Farrar, C. A., & Toscano, N. C. (1988). Relationships among tomato planting date, potato aphids (Homoptera: Aphi- didae) and natural enemies. Journal of Economic Entomology, 81(4), 1107–1112. https://doi.org/10.1093/jee/81.4.1107 Ponder, K. L., Pritchard, J., Harrington, R., & Bale, J. S. (2000). Dif- ficulties in location and acceptance of phloem sap combined with reduced concentration of phloem amino acids explain lowered performance of the aphid Rhopalosiphum padi on nitrogen deficient barley (Hordeum vulgare) seedlings. Ento- mologia Experimentalis et Applicata, 97(2), 203–210. https:// doi.org/10.1046/j.1570-7458.2000.00731.x Powell, G., Tosh, C. R., & Hardie, J. (2006). Host plant selection by aphids: behavioral, evolutionary, and applied perspectives. Annual Review of Entomology, 51(1), 309–330. https://doi. org/10.1146/annurev.ento.51.110104.151107 R Core Team (2014). R: a language and environment for statistical computing. R Foundation for Statistical Computing. http:// www.R-project.org/ Rivelli, A. R., Trotta, V., Toma, I., Fanti, P., & Battaglia, D. (2013). Relation between plant water status and Macrosiphum euphor- biae (Hemiptera: Aphididae) population dynamics on three cultivars of tomato. European Journal of Entomology, 110(4), 617–625. https://doi.org/10.14411/eje.2013.084 Wolfson, J. L., & Murdock, L. L. (1990). Growth of Manduca sexta on wounded tomato plants: role of induced proteinase inhibi- tors. Entomologia Experimentalis et Applicata, 54(3), 257–264. https://doi.org/10.1111/j.1570-7458.1990.tb01336.x https://doi.org/10.1094/MPMI-02-13-0059-R https://doi.org/10.1094/MPMI-02-13-0059-R https://doi.org/10.1094/MPMI-02-13-0059-R https://doi.org/10.1094/MPMI-02-13-0059-R https://doi.org/10.1094/MPMI-02-13-0059-R https://doi.org/10.1094/MPMI-02-13-0059-R https://doi.org/10.1016/j.tree.2005.05.001 https://doi.org/10.1016/j.tree.2005.05.001 https://doi.org/10.1007/s12600-013-0356-3 https://doi.org/10.1007/s12600-013-0356-3 https://doi.org/10.1093/pcp/pcf071 https://doi.org/10.1086/282748 https://doi.org/10.1086/282748 https://doi.org/10.1016/j.molp.2017.02.006 https://doi.org/10.1016/j.molp.2017.02.006 https://doi.org/10.1079/9780851998190.0311 https://doi.org/10.1079/9780851998190.0311 https://doi.org/10.1079/9780851998190.0153 https://doi.org/10.17582/journal.pjz/2017.49.4.1343.1351 https://doi.org/10.17582/journal.pjz/2017.49.4.1343.1351 https://doi.org/10.14411/eje.2009.026 https://doi.org/10.1007/s11829-019-09690-3 https://doi.org/10.1007/s11829-019-09690-3 https://doi.org/10.1046/j.1365-2311.2001.00321.x https://doi.org/10.1046/j.1365-2311.2001.00321.x https://doi.org/10.1016/j.plantsci.2018.04.027 https://doi.org/10.1016/j.plantsci.2018.04.027 https://doi.org/10.1016/B978-0-12-802441-6.00002-4 https://doi.org/10.1093/jee/81.4.1107 https://doi.org/10.1046/j.1570-7458.2000.00731.x https://doi.org/10.1046/j.1570-7458.2000.00731.x https://doi.org/10.1146/annurev.ento.51.110104.151107 https://doi.org/10.1146/annurev.ento.51.110104.151107 http://www.R-project.org/ http://www.R-project.org/ https://doi.org/10.14411/eje.2013.084 https://doi.org/10.1111/j.1570-7458.1990.tb01336.x