Next Article in Journal
A Comprehensive Review of the Coffee Leaf Miner Leucoptera coffeella (Lepidoptera: Lyonetiidae)—A Major Pest for the Coffee Crop in Brazil and Others Neotropical Countries
Next Article in Special Issue
Changes in Alpine Butterfly Communities during the Last 40 Years
Previous Article in Journal
Diversity and Distribution of Forest Ants (Hymenoptera: Formicidae) in Nepal: Implications for Sustainable Forest Management
Previous Article in Special Issue
A Small-Scale Analysis of Elevational Species Richness and Beta Diversity Patterns of Arthropods on an Oceanic Island (Terceira, Azores)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Insect Herbivore Populations and Plant Damage Increase at Higher Elevations

1
Department of Entomology, The Pennsylvania State University, State College, PA 16802, USA
2
Microbial Solutions Team, AgResearch Ltd., Lincoln 7674, Christchurch 8140, New Zealand
3
Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, TN 37996, USA
4
Louisiana Department of Health, Northeast Delta Human Services Authority, Monroe, LA 71201, USA
5
Plantwise, Center for Agriculture and Bioscience International (CABI), South Asia Office, New Delhi 110012, India
*
Author to whom correspondence should be addressed.
Submission received: 12 October 2021 / Revised: 13 December 2021 / Accepted: 15 December 2021 / Published: 17 December 2021
(This article belongs to the Collection Insects in Mountain Ecosystems)

Abstract

:

Simple Summary

It is vitally important to understand the effects of climate change on insect pest populations and crop losses. Using elevation as a proxy for climate change, a field study was conducted in farmer’s fields in Nepal at various elevations in the Himalayan Mountains. At higher elevations, natural herbivore populations and plant damage from herbivory were significantly higher compared to their low-elevation counterparts. Temperature varied with elevation in the field and significantly affected both insect populations and plant damage. A geographical shift of insect pests towards higher elevations is predicated, so it is important to better understand how biotic and abiotic ecological factors and evolutionary processes will act together to affect ecosystem dynamics to reliably predict future pest problems.

Abstract

Elevation gradients are used as a proxy to simulate climate change effects. A field study was conducted along an elevational gradient in Nepal to understand the effects of abiotic conditions on agriculturally important insect herbivore populations (tobacco caterpillar: Spodoptera litura, tomato fruit worm: Helicoverpa armigera, and South American leaf miner, Tuta absoluta) and herbivory damage on tomatoes. Elevation ranged from 100 m to 1400 m above sea level, representing different climatic zones where tomatoes are grown. Contrary to our hypothesis, natural herbivore populations and herbivory damage significantly increased at higher elevations. Individual insect species responses were variable. Populations of S. litura and T. absoluta increased at higher elevations, whereas the H. armigera population was highest at the mid-elevational range. Temperature variations with elevation also affected insect catch numbers and the level of plant damage from herbivory. In the context of climate warming, our results demonstrate that the interactive effects of elevation and climatic factors (e.g., temperature) will play an important role in determining the changes in insect pest populations and the extent of crop losses.

1. Introduction

Elevational gradients are characterized by large variations in environmental conditions over relatively short distances [1]. Therefore, they are used in experiments to simulate changes expected under climate change scenarios [2]. Studies on how insect–plant interactions change with elevation are largely skewed towards woody and perennial herbaceous plants while impacts on agriculture crops are limited [3].
Crops are generally highly responsive to changes in both biotic and abiotic stressors [4,5]. Abiotic factors (e.g., temperature, CO2, UV radiation, air pressure) vary with elevation. In general, a low-level of herbivore abundance and damage is predicted at higher elevations due to harsher environmental conditions [6,7]. In addition, the level of plant secondary compounds (e.g., phenolics, glucosinolates, volatile organic compounds) and leaf morphological characteristics (e.g., leaf toughness and trichome density) can also vary with changing abiotic conditions affecting plant resistance to herbivores along elevation gradients [3,8,9]. A combination of herbivore pressure and plant defense responses can presage the level of plant damage and crop losses in the context of a warming climate [10]. Ambient temperature, which declines at higher elevations, is considered one of the most important factors affecting herbivore and plant traits [11,12].
There is limited knowledge of variation in insect populations and plant damage along elevational gradients in vegetable systems in Nepal. The present study evaluates the changes in natural herbivore populations and herbivory damage in tomatoes (Solanum lycopersicum) along an elevation gradient in Western Nepal. We tested the hypothesis that the extent of insect herbivory will decrease along an elevational gradient. Insect populations were monitored over time using pheromone traps, whereas data on plant damage were measured as a percentage of total leaves damaged from herbivory. Nepal is an ideal country for the research because it is highly diverse in terms of climate, ecosystems, and floral and faunal species; within a north–south trajectory of only 80 km, climatic regimes from the Florida Keys to the Arctic are found [13].

2. Materials and Methods

The study was conducted in Western Nepal, on a south–north, low- to high-land transect that bridges the country’s physiographic regions, from Naubasta, Banke (157 m: 28.2505° N, 81.6729° E) to Gadi, Surkhet (1389 m: 28.6410° N, 81.6136° E), passing through Jhilmile, Surkhet (605 m: 28.4398° N, 81.7363° E) (Figure 1). The gradient spanned an elevation difference of 1232 m with vegetable farms throughout. At each elevation, three replicate tomato plots (50 m2 area) within a 500 m radius were established (Figure S1) on working smallholder farms. The field study was conducted from June to December 2018. Temperature was continuously monitored throughout the experiment using a temperature logger (HOBO MX2300 Temp, Bourne, MA, USA) in each of the three elevations (one per elevational site). Temperature loggers were fitted onto a 5 m-tall stick at a height of 3 m and left in the open sun in an appropriate location at the center of one of the farms. Temperature was recorded automatically every 30 min and weekly mean was used for data analysis.
Tomato is one of the most important vegetables in Nepal, the consumption of which has tremendously increased in recent times [14,15]. Insect herbivore species are key pests of tomatoes in Nepal [16,17]. We assessed leaf herbivory damage as well as populations of tomato fruit worm (Helicoverpa armigera), tobacco cutworm (Spodoptera litura), and South American tomato leaf miner (Tuta absoluta), weekly by monitoring them along an elevational gradient using pheromone traps and crop damage estimates. H. armigera and S. litura are polyphagous pests native to Asia and surrounding regions, whereas T. absoluta is one of the recent invasive pests into the region.
Tomato var. Srijana (Nakkhu Seeds, Kathmandu, Nepal) was used for the experiment, which is a locally obtained popular tomato variety and is recommended for regions up to 2200–2400 m.a.s.l in Nepal. Seedlings (n = 36/site/elevation) were transplanted in black plastic bags (12″ × 20″) with a common potting soilless mixture (vermicompost and sand) to mitigate soil effects. Seedlings were planted during June–July (spacing- Rows × Plants: 75 × 60 cm). Urea, diammonium phosphate (DAP) and murate of potash (MOP) were used at the rate of @ 200:140:140 kg/ha; 50% of urea and 100% of DAP and MOP were mixed with the common potting mixture and two split doses of remaining urea after 15 and 30 days of transplanting. Organic manure prepared from local farm wastes and crop residues was used @ 30,000 kg/ha.
Populations of tomato fruit worm (H. armigera), tobacco caterpillar (S. litura) and South American tomato leaf miner (T. absoluta) were monitored weekly using pheromone traps targeted for these insects (Pest control India, Bangalore, India, https://www.pestcontrolindia.com/products/ (accessed on 21 March 2018)). Funnel traps were used for H. armigera and S. litura and water traps for T. absoluta. The farm settlements were isolated, and the villages where the experiments were conducted were within a 1–2 km area; therefore, one trap per insect species per elevation provided sufficient coverage. Traps were hung in the middle of the research plot at about 2 m above the ground level. There was at least 5 m between traps. Synthetic pheromone lures (Pest control India, Bangalore, India https://www.pestcontrolindia.com/products/ (accessed on 7 April 2018)) for each individual species were replaced every two months as recommended. Traps were checked weekly. Individual target species were identified, and the number trapped was expressed as the average number of insects captured/species/week/ elevation. Trap counts for the three insect species were analyzed separately as well as grouped together to calculate total insects trapped/week/elevation. There were 16 observation dates.
Leaves with characteristic damage by leaf-feeding herbivores and leaf miners (e.g., holes, mines, chewing damage) were considered damaged leaves. Two weeks after transplantation, 10 plants/site/elevation were randomly selected and checked weekly for leaf damage. Two different approaches were used. First, binary data on whether a plant had been damaged or not was recorded by checking individual leaves from a plant for any kind of herbivore damage. Secondly, the percentage (%) leaf damage from herbivory was calculated by counting the number of damaged leaves relative to the total number of leaves from a plant [18,19]. Using the percentage leaf damage approach helps address potential variation due to plant size across sites.
Statistical analysis: The total number of insects trapped (n) per week in different elevations for individual species and the sum of all three species was analyzed with ‘elevation’ as categorical predictors via generalized linear model (GLM) with Poisson distribution [20]. Pair-wise comparisons between elevational sites were carried out using Tukey mean comparisons.
The correlation between the temperature and insect count data was analyzed first by calculating the cross-correlation of the paired time series. The cross-correlation was calculated simultaneously by assuming a non-delayed insect number response to the temperature. This was followed by shifting the insect count series by one lag (one time interval) at a time behind the temperature series to see if there was any delayed response, as indicated by the cross-correlation value larger than the simultaneous cross-correlation value. Once the maximum cross-correlation was achieved, either with or without the time shifting of the insect count series, the identified maximum cross-correlation between each paired series was analyzed in time series regression. The insect count number was the dependent variable, and the temperature was the explanatory variable.
The binary data (plants damaged vs. undamaged) were compared among elevation sites using a GLM with binomial distributions through a logit link function with ‘elevation’ as a single factor. Percentage leaf damage from herbivore data was arcsine transformed to ensure that the distribution of analyzed data was approximately normal. The data from the three elevations were analyzed using analysis of variance (ANOVA) with ‘elevation’ as a single factor and pair-wise comparisons were made using the Fisher’s least significant difference method.
Later, the arcsine-transformed damage data were analyzed by linear regression: Arcsine (Damage) = Elevation + Temperature + Elevation × Temperature, where ‘Elevation’ was a categorical variable representing three elevations while ‘Temperature’ was a numerical variable, and Elevation × Temperature was the interaction effect of the two. Consequently, three different regression lines were fitted in a single analysis, with each line representing the correlation between temperature and plant damage at each elevation. These single regression analyses also allowed statistical comparison between the three regression lines; thus, we tested the elevation-and-temperature interaction effects by comparing the slopes between the lines. Analyses were carried out using Minitab software [21].

3. Results

3.1. Temperature Variation with Elevation

As expected, the temperature was reduced at higher elevations (Figure S1). The mean air temperature from June to December 2018 declined as the elevation increased (i.e., from 25.7 °C at 157 m to 20.2 °C at 1389 m). The lowest temperature recorded at Gadi (1389 m) was 12.9 °C while the highest was 26.0 °C with a median temperature of 21.6 °C. In Jhilmile (605 m), the lowest temperature recorded was 15.0 °C while the highest was 29.7 °C with a median temperature of 25.8 °C. Similarly, the lowest temperature recorded in Naubasta (157 m) was 17.2 °C while the highest was 32.0 °C with a median temperature of 24.8 °C. Despite their average temperature of between 26.7 °C and 32.2 °C, the maximum temperature in Jhilmile and Naubasta consistently exceeded 37.8 °C during daytime (Figure S2).

3.2. Insect Herbivore Populations

The sum of three individual insect species (H. armigera, S. litura, and T. absoluta) trapped per week significantly varied with elevation (χ2 (2, 30) = 54.56, p < 0.01). The overall insect catch number in Gadi (1389 m) was significantly higher than the other two elevations (p < 0.000), but Jhilmile (605 m) and Naubasta (157 m) were not different (Figure 2a). The independent effect of elevation was also significant for S. litura2 (2, 30) = 31.63, p < 0.01) and T. absoluta2 (2, 30) = 17.53, p < 0.01) but not for H. armigera2 (2, 30) = 3.00, p = 0.22). Average populations of S. litura and T. absoluta were highest at Gadi (1389 m) followed by Jhilmile (605 m) and Naubasta (157 m) (Figure 2b). The H. armigera population, however, was highest at Jhilmile (605 m) followed by Gadi (1389 m) and Naubasta.
The number of S. litura caught in pheromone traps and the sum of all three insect species were positively correlated with temperature in Gadi (1389 m) (Table 1). Other correlations were not significant.

3.3. Leaf Damage from Herbivory

There were 915 observations (Gadi: 345, Jhilmile: 345, Naubasta: 225) over the entire study period, of which 299 plants had observable damage from herbivory. Approximately 6% of plants died during the experiment. There was a significant effect of elevation on the percentage of leaves damaged (F (2,296) = 13.37, p < 0.001) but not for the number of plants damaged vs. undamaged. Based on pair-wise comparisons, plants at Gadi (1389 m) had a significantly higher proportion of leaf damage per plant than Naubasta (189 m) and Jhilmile (605 m) (p < 0.001) (Figure 3). Level of damage between Naubasta (189 m) and Jhilemile (605 m) was not significantly different.
There was a significantly positive relationship between temperature and level of plant damage in Gadi (1389m), but they were negatively correlated in Jhilmile (605 m) (Table 2). When compared among elevations, the temperature–plant-damage relationship in Gadi (1389 m) was significantly different (p < 0.001) from Jhilmile (605 m) and Naubasta (189 m). In addition, Jhilmile (605 m) was also statistically (p = 0.033) different from Naubasta (189 m). The elevation-and-temperature interaction effect was statistically significant (F (2,293) = 27.65, p < 0.001) (Figure 4 and Table 2).

4. Discussion

Insect herbivore populations and level of plant damage from herbivory increased with elevation. Comparable results were reported for Antesia Bug, Antestiopsis thunbergia, in Arabica coffee in Tanzania, where the population density of the insects increased with elevation, (1000–1700 m a.s.l) [22]. Herbivore populations and plant damage in Ficus spp. (200 m to 2700 m a.s.l) were also significantly higher at a mid-elevation range (700 to 1200 m a.s.l), which is comparable to our experimental elevation range (605 to 1389 m a.s.l) [23]. These results, however, contrast with the general prediction of greater herbivore damage in plants growing closer to sea level because of a warmer climate [7,24]. However, the generality of this prediction should be interpreted with caution. For example, a correlation between elevational gradients and herbivory was not found in non-woody species in a global analysis of insect herbivory with 1027 plant species [11]. Similarly, Moriera et al. in their review also questioned the hypothesis, providing evidence of several interactions including positive, negative and non-linear associations with elevation [3]. These discrepancies in the findings may have resulted from a large variation in biotic and abiotic factors along an elevational gradient and their independent and combined influences on insect–plant interactions [3,11,25]. Paudel et al. termed this phenomenon an ‘asymmetric response to climate warming’ [26].
Temperature variations within the elevation gradient affected natural herbivore levels and plant damage. This result has important ecological consequences in the context of climate change. An increase in herbivore abundance, as well as a geographical shift of insect pests towards higher elevations, is expected with rising temperature [24,27]. Temperatures in tropical regions for most insects are already in an optimum range; therefore, any further increases may significantly reduce their abundance and diversity [28]. At the same time, increasing temperatures at higher latitudes or temperate regions might benefit from a northward range expansion resulting in higher insect abundance and crop losses [29,30]. Therefore, the potential of insect pests to better adapt will determine species distribution and community interactions in the face of climate change.
We also observed asymmetric elevation-specific variations in individual insect herbivore populations. This may be due to the patterns of temperature changes in these three elevational sites as well as the cropping system. The temperature in Naubasta (157 m) and Jhilmile (605 m) went consistently over 37.8 °C, which is considered a critical temperature threshold for several insect activities: 35.9 °C to 38.5 °C for T. absoluta [31], 38.0 °C for S. litura [32], >37.0 °C for H. armigera [33]. In contrast, the high-elevation site (Gadi, 1389 m) experienced comparatively moderate conditions. Furthermore, tomato production in Nepal is staggered in time, depending on elevation. In the hills (500 to 2000 m), tomatoes are successfully produced year round under plastic houses with open sides [34], which may provide food and shelter to insect herbivores specializing on tomatoes. At lower elevations (<500 m), the production window is limited (Aug–Feb) and is constrained by hot temperatures, low fruit set/flowering and diseases.

5. Conclusions

Overall, the present study suggests that natural herbivore populations and plant damage increase at higher elevations and are influenced by abiotic factors (e.g., temperature). While this study was focused on a single elevational gradient, the results are indicative of some of the changes in insect pest populations and plant damage with global warming. The asymmetric effect of site-specific temperature variations on individual insect pests is indicative of the complexity involved in predicting the climate-change-mediated effects on insect–plant interactions. Future studies including broader and multiple elevational ranges with a complete set of insect pests associated with a crop are necessary to bolster the results and conclusions drawn from this study.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/insects12121129/s1, Figure S1: Average air temperatures from June to December 2018 in three different elevations, Naubasta (157 m), Jhilmile (605 m) and Gadi (1389 m); Figure S2: Variation in air temperatures from May to Nov 2018 in three different elevations, Gadi (1389 m), Jhilmile (605 m) and Naubasta (157 m) as recorded by the temperature logger.

Author Contributions

Conception and design: S.P., G.W.F. and E.G.R.; Data recording, analysis and interpretation: P.K., D.B., V.P., S.P., G.W.F. and E.G.R.; Preparation of manuscript and review: S.P., V.P., G.W.F. and E.G.R.; Obtaining of funding: E.G.R., G.W.F. and S.P. All authors have read and agreed to the published version of the manuscript.

Funding

S.P. was funded through the Integrated Pest Management Innovation Lab (IPM IL), United States Agency for International Development (USAID) Agreement No. AID-OAA-L-15-00001. This research was also partially supported by the US Department of Agriculture Grant AFRI 2017- 67013-26596, National Science Foundation Grant IOS- 1645548, CABI Plantwise, and International Program, College of Agriculture Sciences, Pennsylvania State University.

Data Availability Statement

Data is contained within the article and supplementary material.

Acknowledgments

We would like to dedicate this paper to the late Padam Khadka, who was instrumental in setting up field visits and experiments. He unfortunately lost his life this year. We thank Komal Pradhan, Lalit Sah, Chakra Rai and Muna Basel for their continuous assistance in carrying out the study in Nepal. Statistical help from Chikako van Koten is also highly appreciated.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Körner, C. The use of ‘altitude’in ecological research. Trends Ecol. Evol. 2007, 22, 569–574. [Google Scholar] [CrossRef]
  2. Buckley, J.; Widmer, A.; Mescher, M.C.; De Moraes, C.M. Variation in growth and defence traits among plant populations at different elevations: Implications for adaptation to climate change. J. Ecol. 2019, 107, 2478–2492. [Google Scholar] [CrossRef]
  3. Moreira, X.; Petry, W.K.; Mooney, K.A.; Rasmann, S.; Abdala-Roberts, L. Elevational gradients in plant defences and insect herbivory: Recent advances in the field and prospects for future research. Ecography 2018, 41, 1485–1496. [Google Scholar] [CrossRef] [Green Version]
  4. Smith, S.A.; Beaulieu, J.M. Life history influences rates of climatic niche evolution in flowering plants. Proc. R. Soc. B Biol. Sci. 2009, 276, 4345–4352. [Google Scholar] [CrossRef] [Green Version]
  5. Smilanich, A.M.; Fincher, R.M.; Dyer, L.A. Does plant apparency matter? Thirty years of data provide limited support but reveal clear patterns of the effects of plant chemistry on herbivores. New Phytol. 2016, 210, 1044–1057. [Google Scholar] [CrossRef] [Green Version]
  6. Pellissier, L.; Fiedler, K.; Ndribe, C.; Dubuis, A.; Pradervand, J.-N.; Guisan, A.; Rasmann, S. Shifts in species richness, herbivore specialization, and plant resistance along elevation gradients. Ecol. Evol. 2012, 2, 1818–1825. [Google Scholar] [CrossRef]
  7. Salgado, A.L.; Suchan, T.; Pellissier, L.; Rasmann, S.; Ducrest, A.-L.; Alvarez, N. Differential phenotypic and genetic expression of defence compounds in a plant–herbivore interaction along elevation. R. Soc. Open Sci. 2016, 3, 160226. [Google Scholar] [CrossRef] [Green Version]
  8. Dostdlek, T.; Rokaya, M.B.; Marsik, P.; Rezek, J.; Skuhrovec, J.; Pavela, R.; Munzbergova, Z. Trade-off among different anti-herbivore defence strategies along an altitudinal gradient. AoB Plants 2016, 8, plw026. [Google Scholar] [CrossRef] [Green Version]
  9. Talebi, S.M.; Nohooji, M.G.; Yarmohammadi, M.; Khani, M.; Matsyura, A. Effect of altitude on essential oil composition and on glandular trichome density in three Nepeta species (N. sessilifolia, N. heliotropifolia and N. fissa). Mediterr. Bot. 2019, 40, 81–93. [Google Scholar] [CrossRef] [Green Version]
  10. Pratt, J.D.; Mooney, K.A. Clinal adaptation and adaptive plasticity in Artemisia californica: Implications for the response of a foundation species to predicted climate change. Glob. Chang. Biol. 2013, 19, 2454–2466. [Google Scholar] [CrossRef]
  11. Wang, Q.; Fan, X.; Wang, M. Recent warming amplification over high elevation regions across the globe. Clim. Dyn. 2014, 43, 87–101. [Google Scholar] [CrossRef] [Green Version]
  12. Galmán, A.; Petry, W.K.; Abdala-Roberts, L.; Butron, A.; De La Fuente, M.; Francisco, M.; Kergunteuil, A.; Rasmann, S.; Moreira, X. Inducibility of chemical defences in young oak trees is stronger in species with high elevational ranges. Tree Physiol. 2018, 39, 606–614. [Google Scholar] [CrossRef]
  13. Dyoulgerov, M.; Bucher, A.; Zermoglio, F.; Forner, C. Climate Risk and Adaptation Country Profile. Nepal: Vulnerability, Risk Reduction and Adaptation to Climate Change. Global Facility for Disaster Reduction and Recovery; World Bank: Washington, DC, USA, 2011; Available online: https://www.gfdrr.org/en/publication/climate-risk-and-adaptation-country-profile-nepal (accessed on 15 January 2018).
  14. Ghimire, D.; Lamsal, G.; Paudel, B.; Khatri, S.; Bhusal, B. Analysis of trend in area, production and yield of major vegetables of Nepal. Trends Hortic. 2018. [Google Scholar] [CrossRef]
  15. Acharya, R.; Paudel, S.; Rajbhandari, B.P. Growth, yield formation, fruit quality and economics of tomato production under IPM package vs. conventional practice in Surkhet district. Nepal. J. Agric. Sci. 2015, 13, 96–102. [Google Scholar]
  16. Paudel, S.; Sah, L.P.; Pradhan, K.; Colavito, L.A.; Upadhyay, B.P.; Rajotte, E.G.; Muniappan, R. Development and Dissemination of Vegetable IPM Practices and Packages in Nepal. In Integrated Pest Management of Tropical Vegetable Crops; Springer: Dordrecht, The Netherlands, 2016; pp. 251–269. [Google Scholar]
  17. Venkatramanan, S.; Wu, S.; Shi, B.; Marathe, A.; Eubank, S.; Sah, L.; Giri, A.; Colavito, L.; Nitin, K.; Sridhar, V.; et al. Modeling commodity flow in the context of invasive species spread: Study of Tuta absoluta in Nepal. Crop. Prot. 2020, 135, 104736. [Google Scholar] [CrossRef] [Green Version]
  18. Kanhar, K.A.; Sahito, H.A.; Kanher, F.M.; Tunio, S.A.; Awan, R.R.H. Damage percent and biological parameters of leaf miner Acrocercops syngramma (Meyrick) on different mango varieties. J. Entomol. Zool. Stud. 2016, 4, 541–546. [Google Scholar]
  19. Paudel, S.; Sah, L.; Devkota, M.; Poudyal, V.; Prasad, P.; Reyes, M. Conservation Agriculture and Integrated Pest Management Practices Improve Yield and Income while Reducing Labor, Pests, Diseases and Chemical Pesticide Use in Smallholder Vegetable Farms in Nepal. Sustainability 2020, 12, 6418. [Google Scholar] [CrossRef]
  20. Demétrio, C.G.; Hinde, J.; Moral, R.A. Models for overdispersed data in entomology. In Ecological Modelling Applied to Entomology; Springer: Cham, Germany, 2014; pp. 219–259. [Google Scholar]
  21. Minitab (Version 18). Available online: http://www.minitab.com/en-US/products/minitab/default.aspx (accessed on 22 June 2020).
  22. Azrag, A.G.A.; Pirk, C.W.W.; Yusuf, A.A.; Pinard, F.; Niassy, S.; Mosomtai, G.; Babin, R. Prediction of insect pest distribution as influenced by elevation: Combining field observations and temperature-dependent development models for the coffee stink bug, Antestiopsis thunbergii (Gmelin). PLoS ONE 2018, 13, e0199569. [Google Scholar] [CrossRef]
  23. Sam, K.; Koane, B.; Sam, L.; Mrazova, A.; Segar, S.; Volf, M.; Moos, M.; Simek, P.; Sisol, M.; Novotny, V. Insect herbivory and herbivores of Ficus species along a rain forest elevational gradient in Papua New Guinea. Biotropica 2020, 52, 263–276. [Google Scholar] [CrossRef]
  24. Rasmann, S.; Pellissier, L.; Defossez, E.; Jactel, H.; Kunstler, G. Climate-driven change in plant-insect interactions along elevation gradients. Funct. Ecol. 2014, 28, 46–54. [Google Scholar] [CrossRef]
  25. Hodkinson, I.D. Terrestrial insects along elevation gradients: Species and community responses to altitude. Biol. Rev. 2005, 80, 489–513. [Google Scholar] [CrossRef] [Green Version]
  26. Paudel, S.; Lin, P.-A.; Hoover, K.; Felton, G.W.; Rajotte, E.G. Asymmetric Responses to Climate Change: Temperature Differentially Alters Herbivore Salivary Elicitor and Host Plant Responses to Herbivory. J. Chem. Ecol. 2020, 46, 891–905. [Google Scholar] [CrossRef]
  27. Pateman, R.M.; Hill, J.K.; Roy, D.B.; Fox, R.; Thomas, C.D. Temperature-Dependent Alterations in Host Use Drive Rapid Range Expansion in a Butterfly. Science 2012, 336, 1028–1030. [Google Scholar] [CrossRef]
  28. Johansson, F.; Orizaola, G.; Nilsson-Örtman, V. Temperate insects with narrow seasonal activity periods can be as vulnerable to climate change as tropical insect species. Sci. Rep. 2020, 1, 1–8. [Google Scholar] [CrossRef]
  29. Parmesan, C.; Ryrholm, N.; Stefanescu, C.; Hill, J.K.; Thomas, C.; Descimon, H.; Huntley, B.; Kaila, L.; Kullberg, J.; Tammaru, T.; et al. Poleward shifts in geographical ranges of butterfly species associated with regional warming. Nat. Cell Biol. 1999, 399, 579–583. [Google Scholar] [CrossRef]
  30. Konvicka, M.; Maradova, M.; Benes, J.; Fric, Z.; Kepka, P. Uphill shifts in distribution of butterflies in the Czech Republic: Effects of changing climate detected on a regional scale. Glob. Ecol. Biogeogr. 2003, 12, 403–410. [Google Scholar] [CrossRef] [Green Version]
  31. Özgökçe, M.S.; Bayındır, A.; Karaca, I.; Çobanoğlu, S.; Kumral, N.A. Temperature-dependent development of the tomato leaf miner, Tuta absoluta (Meyrick) (Lepidoptera: Gelechiidae) on tomato plant Lycopersicon esculentum Mill. (Solanaceae). Turk. J. Èntomol. 2016, 40, 51–59. [Google Scholar] [CrossRef] [Green Version]
  32. Fand, B.B.; Sul, N.T.; Bal, S.K.; Minhas, P.S. Temperature Impacts the Development and Survival of Common Cutworm (Spodoptera litura): Simulation and Visualization of Potential Population Growth in India under Warmer Temperatures through Life Cycle Modelling and Spatial Mapping. PLoS ONE 2015, 10, e0124682. [Google Scholar] [CrossRef] [Green Version]
  33. Kriticos, D.J.; Ota, N.; Hutchison, W.D.; Beddow, J.; Walsh, T.; Tay, W.T.; Borchert, D.M.; Paula-Moraes, S.V.; Czepak, C.; Zalucki, M. Correction: The Potential Distribution of Invading Helicoverpa armigera in North America: Is It Just a Matter of Time? PLoS ONE 2015, 10, e0119618. [Google Scholar] [CrossRef]
  34. Pandey, Y.; Pun, A.; Upadhyuy, K. Participatory Varietal Evaluation of Rainy Season Tomato under Plastic House Condition. Nepal Agric. Res. J. 2006, 7, 11–15. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Research sites in Nepal—Naubasta, Banke (157 m; 28.2505° N, 81.6729° E), Jhilmile, Surkhet (605 m; 28.4398° N, 81.7363° E), and Gadi, Surkhet (1389 m; 28.6410° N, 81.6136° E).
Figure 1. Research sites in Nepal—Naubasta, Banke (157 m; 28.2505° N, 81.6729° E), Jhilmile, Surkhet (605 m; 28.4398° N, 81.7363° E), and Gadi, Surkhet (1389 m; 28.6410° N, 81.6136° E).
Insects 12 01129 g001
Figure 2. Insects trapped in pheromone traps at three different elevations: Naubasta (157 m), Jhilmile (605 m) and Gadi (1389 m); (a) box-plot diagram showing variation in average number of total insect trapped including H. armigera, S. litura and T. absoluta (total insect trapped (n)/week) and (b) box-plot diagram showing variation in the average number of insects trapped/week/elevation delineated by each individual herbivore species. The segment inside the rectangle represents the median insect trapped and “whiskers” above and below the box correspond to the minimum and maximum number of insects trapped. Circles (〇) indicate outliers.
Figure 2. Insects trapped in pheromone traps at three different elevations: Naubasta (157 m), Jhilmile (605 m) and Gadi (1389 m); (a) box-plot diagram showing variation in average number of total insect trapped including H. armigera, S. litura and T. absoluta (total insect trapped (n)/week) and (b) box-plot diagram showing variation in the average number of insects trapped/week/elevation delineated by each individual herbivore species. The segment inside the rectangle represents the median insect trapped and “whiskers” above and below the box correspond to the minimum and maximum number of insects trapped. Circles (〇) indicate outliers.
Insects 12 01129 g002aInsects 12 01129 g002b
Figure 3. Box-plot diagram showing leaf damage by herbivores (%) in tomato plants from three different elevations: Naubasta (157 m), Jhilmile (605 m) and Gadi (1389 m). The segment inside the rectangle shows the median and “whiskers” above and below the box show the minimum and maximum damage (%), respectively. Circles (〇) symbolize outliers.
Figure 3. Box-plot diagram showing leaf damage by herbivores (%) in tomato plants from three different elevations: Naubasta (157 m), Jhilmile (605 m) and Gadi (1389 m). The segment inside the rectangle shows the median and “whiskers” above and below the box show the minimum and maximum damage (%), respectively. Circles (〇) symbolize outliers.
Insects 12 01129 g003
Figure 4. Scatterplot showing correlations between temperature and plant damage (leaf damage by herbivores/plant (%) in three different elevations: Gadi (1389 m), Jhilmile (605 m) and Naubasta (157 m). r2 coefficient and associated p-values are shown. Each circle within the plot represents % of leaves damaged by herbivores per plant per week. Different colors represent different elevation sites (blue: Gadi; red: Jhilmile, green: Naubasta).
Figure 4. Scatterplot showing correlations between temperature and plant damage (leaf damage by herbivores/plant (%) in three different elevations: Gadi (1389 m), Jhilmile (605 m) and Naubasta (157 m). r2 coefficient and associated p-values are shown. Each circle within the plot represents % of leaves damaged by herbivores per plant per week. Different colors represent different elevation sites (blue: Gadi; red: Jhilmile, green: Naubasta).
Insects 12 01129 g004
Table 1. Regression analysis of cross-correlation between temperature and insect catch number in different elevational sties.
Table 1. Regression analysis of cross-correlation between temperature and insect catch number in different elevational sties.
Elevation SiteInsect SpeciesTime Series Regression EquationSlope (Mean ± SE) Interceptr2 ValueF-Statisticsp-Value
Gadi (1389)H. armigera + S. litura + T. absolutaTrap catch number = −357.40 + 19.331 × Temperature19.331 ± 7.274−357.40 ± 167.100.414F (1,10) = 7.060.024
Jhilmile (605 m)H. armigera + S. litura + T. absolutaTrap catch number = 190.90 − 5.105 × Temperature−5.105 ± 6.453190.90 ± 176.000.073F (1,8) = 0.630.452
Naubasta (189 m)H. armigera + S. litura + T. absolutaTrap catch number = −63.88 + 3.621 × Temperature3.621 ± 1.730−63.88 ± 49.130.354F (1,8) = 4.380.070
Gadi (1389)H. armigeraTrap catch number = −9.36 + 0.982 × Temperature0.982 ± 0.818−9.36 ± 18.800.126F (1,10) = 1.440.258
Jhilmile (605 m)H. armigeraTrap catch number = 27.58 − 0.516 × Temperature−0.516 ± 0.684 27.58 ± 18.480.054F (1,10) = 0.570.469
Naubasta (189 m)H. armigeraTrap catch number = 30.83 − 0.668 × Temperature−0.668 ± 0.43230.83 ± 12.330.193F (1,10) = 2.390.153
Gadi (1389)S. lituraTrap catch number = −167.86 + 8.681 × Temperature8.681 ± 3.063−167.86 ± 70.380.445F (1,10) = 8.030.0018
Jhilmile (605 m)S. lituraTrap catch number = −31.45 + 1.987 × Temperature1.987 ± 3.647−31.45 ± 98.860.032F (1,10) = 0.300.599
Naubasta (189 m)S. lituraTrap catch number = −45.45 + 2.262 × Temperature2.262 ± 1.082−45.45 ± 30.730.353F (1,8) = 4.370.070
Gadi (1389)T. absolutaTrap catch number = −180.20 + 9.669 × Temperature9.669 ± 5.520−180.20 ± 126.800.235F (1,9) = 3.070.110
Jhilmile (605 m)T. absolutaTrap catch number = 98.44 − 3.040 × Temperature−3.040 ± 2.75198.44 ± 75.010.132F (1,10) = 0.300.599
Naubasta (189 m)T. absolutaTrap catch number = −20.10 + 1.004 × Temperature1.004 ± 0.533−20.10 ± 15.120.308F (1,8) = 3.560.096
Table 2. Elevation-level-specific relationships between temperature (°C) and arcsine-transformed leaf damage percentages estimated from single regression analysis with r2 value.
Table 2. Elevation-level-specific relationships between temperature (°C) and arcsine-transformed leaf damage percentages estimated from single regression analysis with r2 value.
Gadi (1389 m)
EquationArcsine(Damage) = − 120.323 + 6.005 × Temperature
Slope6.005 ± 0.769 (mean ± SE); F = 60.91, df = 1, 293, p < 0.001
Intercept−120.323 ± 17.532 (mean ± SE)
Jhilmile (605 m)
EquationArcsine(Damage) = 47.173 − 1.330 × Temperature
Slope−1.330 ± 0.675 (mean ± SE); F = 3.88, df = 1, 293, p = 0.050
Intercept47.173 ± 18.399 (mean ± SE)
Naubasta (189 m)
EquationArcsine(Damage) = −1.956 + 0.457 × Temperature
Slope0.457 ± 0.489 (mean ± SE); F = 0.87, df = 1, 293, p = 0.351
Intercept−1.956 ± 13.148 (mean ± SE)
r2 = 0.251
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Paudel, S.; Kandel, P.; Bhatta, D.; Pandit, V.; Felton, G.W.; Rajotte, E.G. Insect Herbivore Populations and Plant Damage Increase at Higher Elevations. Insects 2021, 12, 1129. https://0-doi-org.brum.beds.ac.uk/10.3390/insects12121129

AMA Style

Paudel S, Kandel P, Bhatta D, Pandit V, Felton GW, Rajotte EG. Insect Herbivore Populations and Plant Damage Increase at Higher Elevations. Insects. 2021; 12(12):1129. https://0-doi-org.brum.beds.ac.uk/10.3390/insects12121129

Chicago/Turabian Style

Paudel, Sulav, Pragya Kandel, Dependra Bhatta, Vinod Pandit, Gary W. Felton, and Edwin G. Rajotte. 2021. "Insect Herbivore Populations and Plant Damage Increase at Higher Elevations" Insects 12, no. 12: 1129. https://0-doi-org.brum.beds.ac.uk/10.3390/insects12121129

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop