Next Article in Journal
Ameliorative Effects of Phytochemical Ingestion on Viral Infection in Honey Bees
Previous Article in Journal
Bacillus thuringiensis Bioinsecticides Induce Developmental Defects in Non-Target Drosophila melanogaster Larvae
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Insects
Volume 11
Issue 10
10.3390/insects11100695
insects-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Aphid Density and Plant Taxa on Predatory Ladybeetle Abundance at Field and Landscape Scales

by
Hongsheng Pan
1,2,
Bing Liu
1,
Coline C. Jaworski
3,
Long Yang
1,
Yongqiang Liu
1,
Nicolas Desneux
4,
Eva Thomine
4 and
Yanhui Lu
1,*
1
State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing 100193, China
2
Scientific Observing and Experimental Station of Crop Pests in Korla, Ministry of Agriculture and Rural Affairs, Institute of Plant Protection, Xinjiang Academy of Agricultural Sciences, Urumqi 830091, China
3
Department of Zoology, University of Oxford, Oxford OX1 3PS, UK
4
Université Côte D’Azur, INRAE, CNRS, UMR ISA, 06000 Nice, France
*
Author to whom correspondence should be addressed.
Insects 2020, 11(10), 695; https://doi.org/10.3390/insects11100695
Submission received: 5 September 2020 / Revised: 7 October 2020 / Accepted: 9 October 2020 / Published: 13 October 2020
Browse Figures
Versions Notes

Abstract

:

Simple Summary

In agroecosystems, predatory ladybeetles play an important role in suppressing aphid populations. How ladybeetles make use of host plant diversity in multicropping landscapes has rarely been documented in China. In this study, we examined the relationship between aphid densities and ladybeetle densities at both the local field and landscape scales. Overall, we found that there was a positive correlation between aphid densities and ladybeetle densities. However, plant taxa had no significant influence on predatory ladybeetle abundance at the local field scale. In addition, the effect of aphids on ladybeetles abundance was influenced by the crop type and growing season at the regional landscape scale. There was a significant positive correlation between aphid and ladybeetle populations on cotton only in July and August, whereas the correlation was significant for maize throughout the whole growing season. The δ13C value indicated that most prey aphids for ladybeetles originated from crops where aphids are abundant (cotton in June and July; both maize and cotton in August). These findings improved our understanding of the migration and dispersal of ladybeetles among different habitats and plant species and provided insight into the promotion of regional conservation and pest control of natural enemies in Northern China.

Abstract

In agroecosystems, predatory ladybeetles play an important role in restraining aphid population growth and suppressing aphid populations. They can adapt to various habitats and make use of various aphid species associated with multiple host plants during their life cycle. Agricultural landscapes in China are composed of a mosaic of small fields with a diverse range of crops, and how ladybeetles make use of host plant diversity in such landscapes has rarely been documented. In this study, we examined the relationship between aphid densities and ladybeetle densities in two different settings: (i) on the majority of plant species (including crops, trees, and weeds) at a local field scale in 2013 and 2014, and (ii) in paired cotton and maize crop fields at a regional landscape scale in 2013. Overall, we found that aphid abundance determined predatory ladybeetle abundance at both the local field and landscape scales, and there was a positive correlation between aphid densities and ladybeetle densities. However, plant taxa had no significant influence on the predatory ladybeetle abundance at the local field scale. In addition, the effect of aphids on ladybeetles abundance was influenced by the crop type and growing season at the regional landscape scale. There was a significant positive correlation between aphids and ladybeetles populations on cotton only in July and August, whereas the correlation was significant for maize throughout the whole growing season. We also conducted an analysis of the stable carbon isotope ratios of the adult ladybeetles caught in cotton and maize fields (C3 and C4 crops, respectively) in a regional landscape-scale survey in 2013. The δ13Cvalue indicated that most prey aphids for ladybeetles originated from crops where aphids are abundant (cotton in June and July; both maize and cotton in August).These findings improved our understanding of the migration and dispersal of ladybeetles among different habitats and plant species and provided insight into the promotion of the regional conservation and pest control of natural enemies in northern China.

1. Introduction

Both adult and larvalladybeetles are important predators of various aphid speciesin agroecosystems [1,2]. Combinations of laboratory, greenhouse, and field studies carried out in multiple agricultural systems have reported the contributions and importance of predatory ladybeetles to both decrease aphid population growth and lower peak aphid densities, making them essential biological control agents [3,4,5,6]. For example, they can cause strong, season-long suppression of aphid populations in soybean fields [7]. Many factors could affect the population occurrence of predatory ladybeetles and their potential ecosystem services through pest suppression, but prey abundance is the most important and direct component [1,8,9]. Herbivores are highly dependent on their hosts—various types of plants, such as crops, forests, grass, or weeds, or other phytogroups; plant diversity and functionality could directly affect the diversity and abundance of herbivores and usually influence the higher trophic level-predatory enemies through interactions between plant–herbivore–predator [10,11,12,13].
Most of generalist predatory ladybeetle adults and larvae rely on aphids as a large part of their diets, although they may occasionally use alternative floral resources such as plant pollen and nectar at low aphid density [14,15,16,17]. Adult ladybeetles usually disperse to various habitats and lay eggs on diversifyinghost plant species to utilize aphid preys during the same season [9], and their mobility is critical in enabling them to cope with transient and unpredictable food resources that are scattered in both space and time [9,18,19,20,21]. Aphidophagous ladybeetle adults are highly mobile but become less active and lay eggs at sites with high aphid densities to provide adequate food resources to larvae [1,22,23,24,25]. However, landscape composition and plant diversity were likely to affect their foraging behavior and capacity to locate food resources [26,27,28,29]. The quality of landscape habitat, usually determined by the availability of food resources, is a major driver for the movement of natural enemies among crops [30,31]. Compared with the monoculture landscape, the incidence for natural enemies’ access to high quality habitat appear much higher within a complex one, which is beneficial for the biocontrol services contributed by natural enemies [32,33].
Carbon isotope analyses have provided insights into ladybeetle movements across fields with different plant types. Indeed, plants using a C3 versus a C4 photosynthetic pathway have a different ratio of 13C over 12C (i.e., δ13C), and this signature is transferred to the body of insect herbivores and their own insect predators [34,35,36,37,38]. Because of their low mobility, apterous aphids’ carbon isotope signature is entirely determined by their host plant [39]. Therefore, the combined information on ladybeetle and aphid locations and on ladybeetle carbon isotope signatures makes it possible to estimate ladybeetle movement between habitats composed of C3 and C4 plants [19] and thus to obtain insights into their host plant use.
In the multicropping system of northern China, Harmonia axyridis Pallas, Propylea japonica Thunberg, and Coccinellaseptempunctata L. are commongeneralist predatory ladybeetle species that play a major role in suppressing pest populations, especially those of various key aphid species [5,40,41]. For example, ladybeetles can effectively delay the establishment and subsequent population growth of the cotton aphid Aphis gossypii Glover during the cotton growing season [42], and P. japonica adults can aggregate on plants with high aphid density in agricultural ecosystems composed of cotton and maize [38]. A better knowledge of the responses of ladybeetles to aphid densities and of their temporal use of different host plant species will help promote their conservation and improve biocontrol services by preserving preferred habitats and enhancing their movements between crop fields.
We hypothesized that the aphid density and plant taxa is a key determinant of the population abundanceof predatory ladybeetles in the agroecosystems. We first surveyed the relationship between aphid densities and ladybeetle densities in habitats with the majority of plant species (including crops, trees, and weeds) locally present at a local field scale in 2013 and 2014, and the survey was simultaneously conducted in paired cotton and maize fields (two important crops with overlapping growth periods in northern China) at a regional landscape scale in 2013. In addition, we measured the movement of the ladybeetles by calculating the stable carbon isotope ratios of the adult ladybeetles collected in cotton and maize fields from our regional landscape-scale survey. Our study provides insight into how aphids and plant taxa drive ladybeetle occurrence or migration in various habitats and at various spatiotemporal scales in agricultural ecosystems.

2. Materials and Methods

2.1. Relationship between Aphid and Ladybeetle Abundance at the Field Scale

We examined the correlation between aphid densities and ladybeetle abundance in a mosaic, multicropping local farmland (approx. 50 ha.) near the Langfang Experimental Station, Chinese Academy of Agricultural Sciences (CAAS) (39.53° N, 116.70° E), in Hebei Province, China, from early May to early September in 2013 and 2014. This farmland corresponded to the typical local small-scale multicropping pattern, with multiple small fields (0.2 to 0.5 ha.) of each crop. Overall, wheat was the main crop in spring and early summer, and maize was the dominant crop in summer, along with many minor vegetable crops (e.g., tomato, cucumber, Chinese cabbage), fruit trees (e.g., Chinese date, grape, peach), and timber trees (e.g., poplar). A total of 96 common plant species in the farmland were sampled in 2013 (78 species) and 2014 (80 species) (Table S1). They were classified into three plant functional groups (called “Plant_FG”), i.e., crop plants (31 species), trees (12 species), and weeds (53 species) (Table S1). Plant species were identified using regional weed guides [43] or with the expertise of CAAS plant taxonomists. Aphids and ladybeetles were counted every 10 days and, in total, 13 times each year, on 10 to 20 samples per plant species following Pan et al. [44]. Briefly, sampling consisted of visual plant inspections to count the number ofaphids, as well as adult and larval ladybeetles on the same plant. The identity of adult and larval ladybeetles collected was assessed based on morphological features [45]. For abundant herbaceous species, a sample consisted of a total area of 2 to 20 m2, while for scattered herbaceous species with small distributions, all plants in a patch ranging from 0.02 m2 to 0.50 m2 were sampled. For trees, 10 young branches were randomly sampled (length of 30 cm). The population densities of the aphids or ladybeetles were unified by upscaling or downscaling the sampling area (the crops and weeds were mapped to 1 m2, and 10 branches per tree weresampled) of each plant species on each sampling date.

2.2. Relationship between Aphid and Ladybeetle Abundancein Paired Cotton and Maize Fields at the Regional Landscape Scale

We assessed the effects of aphids on the occurrence of predatory ladybeetle metapopulations on multiple scales. At the regional landscape scale, we surveyed 83 study sites distributed aroundLangfang city and Xiongxian County (belonging to Baoding city) in Hebei Province and around the cities of Wuqing and Jinghai in Tianjin Province in northern China in 2013 (see Figure S1). The distance between the two study sites was averagely 3–4 km. In each study site, two adjacent crop fields were cultivated with cotton (a C3 plant, crop variety “GuoXin”) and maize (a C4 plant, crop variety “ZhengDan958”), and each had an area greater than 5000 m2.The aphid and ladybeetle densities were measured three times: mid-June, mid-July, and mid-August, which corresponded to the budding, flowering, and bolling stages of cotton, and the seedling, whorl, and silking stages of maize. No insecticide was sprayed during the trials or study periods (preventing any side effects of such chemicals) [46]. In each crop field and for each of the three sampling times, 50 plants were sampled over five random locations per field (10 plants per location), and we counted the number of adult and larval ladybeetles of each species and mixed aphid populations (i.e., Rhopalosiphummaidis (Fitch), Macrosiphum miscanthi (Takahashi), Schizaphisgraminum (Rondani), R. padi L.) on maize plants and A. gossypii on cotton plants through visual inspection of 50 plants in the fields [40]. For each species, all insects counted on the 50 plants were summed and doubled to obtain a single value per 100 plants per crop field and sampling date, and then the mean number of insects per 100 plants with the standard deviation was calculated for 83 study sites.
During each field survey, adult ladybeetles (at most 10 individuals per site) were collected from neighboring cotton and maize fields (300 m2 for each crop, at least 20 m away from survey locations to decrease the disturbance to the ladybeetle population abundance) to characterize their diet composition and host plant use, which were assessed through the measurement of their carbon isotope ratio. Only P. japonica (the most abundant species at the landscape scale, see results) was used. The collected individuals were placed into a 1.5 mL centrifuge tube containing 95% ethanol and stored at −20 °C for further stable isotope analysis following the method of Ouyang et al. [38] to assess the movement of this ladybeetle species between these two adjacent crops.

2.3. Statistical Analysis

2.3.1. Field Scale

First, for each study year, we deployed a nonparametric test (Kruskal–Wallis test using “proc npar1way” of the SAS 9.4 software) to compare the difference in non-normality data of aphid density (AD) and ladybeetle metapopulation density (LMD)among these three plant taxonomic groups (crops, trees, weeds) across the whole sampling season. Then, we performed generalized linear mixed-effects models (GLMMs) to clarify the linear or nonlinear relationship between the response variable (LMD) and various explanatory variables [47], the model was fitted with negative binomial distributions due to the data overdispersionderived from the population fluctuation of aphids or ladybeetles on different plant species. The fixed effects included AD, a categorical variable of plant functional groups(Plant_FG with three levels: crops, trees, weeds) and their interaction, the sampling number of plant species (Plant_num) nested in each taxonomic group was the random effect. We first ran the full model which included all fixed effects, then, we removed nonsignificant effects by stepwise model selection, and the final model only contained one fixed effect of AD. Before the analysis, the abundance datawere log-transformed to decrease the model error estimated. GLMMs were fitted with the “glmer.nb” function in the “MASS” package [48] of the R 3.5.3 software [49].

2.3.2. Regional Landscape Scale

At the regional landscape scale, we first ran a repeated-measures analysis with a mixed effects linear model (“proc mixed” within the SAS 9.4 software) to assess the variationin aphid density (AD) and ladybeetle metapopulation density (LMD) on cotton and maize in three investigative periods from June to August 2013. Month (with three levels: June, July, and August) was the repeated factor; crop (with two levels: cotton and maize), month and their interaction (crop*month) were the fixed effects; and the study site was the random effect. Then, in order to assess the response relationship between LMD and AD at the regional landscape scale, GLMM analysis with the Poisson distribution of non-overdispersion data using the R package “lme4” [50] was deployed. The fixed effects included AD, crop, study month, and the interactions of different explanatory variables, the study site was the random effect. Finally, because we were interested in how aphid density affected ladybeetle density, and the fixed effects were confirmed to be significant, all these main effects are shown in the tendency diagram. Prior to all analyses, the abundance data were transformed by log10(x + 1) to decrease the model error estimated.
Moreover, to study the temporal movement of ladybeetles between cotton and maize, stable carbon isotope analysis was used to detect the feeding activity of the dominating species, P. japonica adults, on aphids from cotton and/or maize fields, and the δ13 values were determined for the field-collected P. japonica samples (per crop and per sampling date). The criteria for estimating the proportion of aphidsin the diets of P. japonica adults that were from C3 or C4 plants (i.e., cotton or maize) followed the existing linear equation Y = 0.12X − 22.72, where Y is the δ13 value of P. japonica and X is the proportion of aphids from C3 plants and C4 plants [38].
At both field and landscape scales, we did a Chi-square test (“proc freq” within the SAS 9.4 software) to clarify the species composition proportion of aphidophagousladybeetles in different plant functional groups and crops, respectively.

3. Results

3.1. Species Composition of Aphidophagous Ladybeetles

There werethreeaphidophagous ladybeetle species (C.septempunctata, P. japonica, H. axyridis) found at both the local field and landscape scales (Figure S2). H. axyridis was the dominant ladybeetle specieson all plants at the local field scale, with proportions of 68.6% and 70.2% in 2013 and 2014, respectively; however, no significant difference was found for the species composition of aphidophagous ladybeetles in the different plant functional groups at the field scale in 2013 (χ2 = 0.65, df = 2, p = 0.723, Figure S2a) and 2014 (χ2 = 1.52, df = 2, p = 0.468, Figure S2b). More interesting, we found P. japonica was the dominant ladybeetle specieson these two main crops—cotton (97.6%) and maize (86.2%)—at the regional landscape scale in 2013, and significant differences were found for the species compositions of aphidophagous ladybeetles in different sampling months on cotton (χ2 = 8.96, df = 2, p = 0.011, Figure S2c) and maize (χ2 = 22.26, df = 2, p < 0.001, Figure S2d).

3.2. Effects of Aphids and Plant Taxa on Ladybeetles at the Field Scale

At the field scale, both AD (aphid density) and LMD (ladybeetle metapopulation density) showed obvious seasonal dynamics from 2013 to 2014 (Figure 1). The ADs on three plant functional groups (i.e., crops, trees, and weeds) wereall higher from mid-May to early June in 2013 (Figure 1a), but the AD was only higheron weeds in mid-May in 2014 (Figure 1b). Meanwhile, the LMD on trees was higher in mid-May, and the peak value on crop plants was in early June in 2013 (Figure 1c), whereas the peak values on trees and weeds appeared in early May and late May in 2014, respectively(Figure 1d). However, there wereno significant differences inAD and LMD among the threeplant functional groups throughout the whole sampling period at the local field scale in 2013 (AD: χ2 = 0.75, df = 2, p = 0.688; LMD, χ2 = 4.63, df = 2, p = 0.099) and 2014 (AD: χ2 = 2.12, df = 2, p = 0.347; LMD, χ2 = 2.22, df = 2, p = 0.330) (small plotsnested in Figure 1), respectively.
For the GLMMs, no significant impacts were found for the interaction of aphid density (AD) and the three plant functional groups (Plant_FG) in 2013 (LRT chi-square test = 1.68, df = 2, p = 0.432) or 2014 (LRT chi-square test = 1.47, df = 2, p = 0.480) (AD: Plant_FG in Table S2, Figure S3a,b). However, the LMD wassignificantly positively related to AD across all plants in both years (2013: Wald Z = 4.14, p < 0.001, Figure 2a; 2014: Z = 5.62, p < 0.001, Figure 2b). These results indicated that aphid density was the main determinant of aphidophagouspredatory ladybeetle abundances at the local field scale.
Based upon the occurrence and distribution of plant species as well asaphid density on crops in northern China, the host plant shifts of aphidophagousladybeetles were determinedduring three successive periods atlocal farmland. From early May to early June, ladybeetles mainly remained on species including Triticum aestivum L., Prunus persica L., Malus pumila Mill., Descurainiasophia (L.) Webb. ex Prantl, Hemisteptalyrata Bunge, and Populustomentosa Carr. (seedling). From mid-June to early August Gossypium hirsutum L., Vitis vinifera L., Artemisia lavandulaefolia DC., A. annua L., P. tomentosa (seedling), etc., were the dominant habitat plants of ladybeetles. In the last periods, i.e., from mid-August to early September, a large number of predatory ladybeetles moved to Zea mays L., A. lavandulaefolia, Chenopodium album L., Ulmuspumila L., etc. (Figure S4).

3.3. Effects of Aphids on Ladybeetles at the Regional Landscape Scale

At the regional landscape scale, the AD varied throughout the months (F = 81.48, df = 2, 328, p < 0.001) on both cotton and maize and was higher on cotton than maize throughout the whole sampling period (F = 941.03, df = 1, 82, p < 0.001). The AD on cotton decreased from June to July and August but that on maize increased from June and July to August (Figure 3a). The mean LMDabundance was higher on cotton than on maize throughout the whole sampling period (F = 10.97, df = 1, 82, p = 0.001). The LMD on maizewas also higher in August than in June (t = 16.18, p < 0.001) and July (t = 10.2, p < 0.001), but we did not find a significant decrease in the LMD on cotton (August vs June: t = 0.39, p = 0.694; August vs July: t = 1.38, p = 0.168) (Figure 3b). There were also significant differences in the AD or LMD between cotton and maize in each sampling month (AD, June: t = 28.77, p < 0.001; July: t = 22.24, p < 0.001; August: t = 4.52, p < 0.001. LMD, June: t = 9.72, p < 0.001; July: t = 2.39, p = 0.017; August: t = −6.00, p < 0.001). Moreover, LMD was enhanced following the increasing AD on maize, which indicated a following relationship between the ladybeetle andaphid population.
Furthermore, the GLMM analysis also revealed the relationship between the occurrence of ladybeetles and aphids at the regional landscape scale (Table 1). LMD was significantly positively related to AD (Wald Z = 3.93, p < 0.001), and the interactions such as AD:crop (LRT chi-square test = 100.18, df = 1, p < 0.001) and AD:month (LRT chi-square test = 223.38, df = 2, p < 0.001) were both significant, which indicated the positive effects of AD on LMD were also significant between different crops and different sampling periods. For individual crop, the positive effect of the AD on LMD was only significant on maize (Z = 91.37, p < 0.001, AD:crop-maize interaction) rather than cotton among different sampling periods (from June to August) (Figure 4). The positive relationship between LMD and AD was also significant in July (Z = 6.01, p < 0.001, AD:month-July interaction) and August (Z = −87.56, p < 0.001, AD:month-August interaction) on both cotton (Figure 4a) and maize (Figure 4b).
In the stable carbon isotope analysis, the δ13Cvalues of the P. japonica adults from the cotton and maize fields showed that they were nearly all positive (~100%) for consumption of aphids originated from C3plants in June and July. In August, ladybeetles in cotton and maize are mostly resident in the crop but there is a slight migration between the two crops, i.e., 13.43% of ladybeetles in cotton came from maize (C4plants) and 28.29% in maize came from cotton (C3plants) (Table 2). These results indicated thatmost prey aphids for ladybeetles originated from crops where aphids areabundant (cotton in June and July; both maize and cotton in August), and ladybeetle density had a positive response to aphid density on these two crops following their growing period.

4. Discussion

Prey densities can affect the host plant selection and population growth of predatory ladybeetles [1,8]. In our survey, we provided further evidence that aphid densities are a strong determinant of ladybeetledensities in the farmlands of northern China, and that the host plant type also affected ladybeetle densities at a regional landscape scale. Especially at the local field scale, the trees hosted high ladybeetle densities in the middle of spring, while the crops hosted high ladybeetle densities at the end of spring. The stable carbon isotope analysis at the regional landscape survey confirmed that ladybeetles moved between crops following aphid population fluctuations.
In our two progressive surveys (local field scale and regional landscape scale), the ladybeetle densities increased with increasing aphid densities, as expected and consistent with earlier literatures [9,51,52,53,54]. The ladybeetles responded differently depending on the host plant type: they responded to high aphid densities in maize fields but had a less obvious response to those in cotton fields at the regional landscape scale, which was similar to the results of Ouyang et al. [38], and this response possibly occurred through a phenology effect/host shift, which is discussed below. At the local field scale in 2013, the ladybeetles responded more strongly to aphid densities in trees and crops than weeds. This may be due to the difference in the compensation of plant functional diversity [55]. In addition to prey, other plant nutrients can also benefit predators. For instance, Bertrand et al. [56] reported that two generalist predators of crop aphids (a ladybeetle species: H. axyridis, and a lacewing species: Chrysoperlacarnea) could utilize pollen sources from trees (e.g., Salix, Prunus, Quercus, and Acer) to survive during seasonal shifts in resources in agricultural landscapes. Plant-provided food supplements can affect diversity and biological pest control by omnivorous predators [57]. Some non-crop plants could attract and conserve aphid predators in specific crop fields [58], and plant morphology may affect ladybeetle mobility, foraging behavior, and access to prey, encouraging them to forage on preferred plant types [59,60]. For instance, Reynolds and Cuddington [61] showed that H. axyridis foraged more thoroughly and had a higher foraging ability on plants with more branches compared with leafier plants. H. axyridis also preferred to aggregate in the more sunny section of tree crowns [62], and they were of a higher abundance on trees than on herbaceous plants and cereals in central Europe [63,64]. Hence, trees would be a preferred habitat for them. Different plant species attacked by aphids may also be differentially attractive to ladybeetles due to differences in susceptibility to herbivorous attacks and in the emission of herbivore-induced plant volatiles [65,66,67,68,69]. Finally, the area occupied by the focal plant species in the local landscape may also affect ladybeetles’ foraging behavior via a concentration effect, i.e., plant species occupying large areas such as crops and trees may represent large, uniform, and easily detected prey reservoirs and hence be preferred by ladybeetles [70,71].In addition, a more complex landscape not only can increase plant functional diversity, but can also create the microclimate to influence the thermal tolerance of aphid, beetles, and parasitoids, which in turn could impact biological control strategies [72,73,74]. As the global climate warms, insectpests and their biological control agents maybe a mismatch, as suggested by Tougeron et al. [73] on aphid parasitoids, landscape complexity however can to some extent alleviate the influence of climatechange on the interaction of prey–predator.
In the regional landscape-scale survey, the stable carbon isotope analysis showed that adult ladybeetles moved from aphid-depleted patches to patches with high aphid density: the P. japonica adults collected in maize fields in June and July mostly preyed upon aphids originating from C3plants (likely cotton, the most abundant C3 plant in the environment). This suggested that they hunted for prey in the fields with the highest aphid densities, i.e., cotton fields in June and July. In August, the percentage of aphids consumed by P. japonica adults that had developed on C4plants (likely maize, the most abundant C4 plant in the environment) increased when compared with June and July. This result was consistent with the population growth of aphids in maize fields and the concomitant decrease in aphid densities in cotton fields from July to August. Similar prey-motivated changes in habitat selection by ladybeetles tracking high aphid densities have been observed for other ladybeetle species [1,21,75]. Additionally, crop phenology or temperature shifts can cause the population movement of natural enemies, for instance, maizeplants at whorl stage offer suitable micro-climatic conditions forforaging ladybeetlesand could act as a refuge during seasonally hot or dry conditions in agro-ecosystems of northern China [76]. Thus, global landscape management such as increasing plant diversification could improve conservation biological control in agro-ecosystems [77,78].
Evidence for temporal host plant shifts was also found at the local field scale, with ladybeetle densities decreasing in trees over time in 2013 and 2014 and increasing versus decreasing in weeds over time in 2013 and 2014, respectively, compared with in crops. Hence, trees may represent early season reservoirs forladybeetles before colonization of crops or weed patches [62]. However, the aphid densities in crops were already high in early May when many ladybeetles were in trees, which also supported high aphid densities; therefore, their shift from trees to crops may happen too late for efficient biocontrol of aphid populations in crop field [79,80]. There were also likely temporal shifts within plant types (e.g., from cotton to maize) that could not be tested with our models. Seasonal migrations of ladybeetles across host plants have often been described and proposed as a means to enhance aphid biocontrol via habitat enhancement and the conservation of reservoirs [3,8].
Interspecific competition between ladybeetle species determines their foraging behavior, especially in patch selection for oviposition [53,81,82,83]. In the present study, ladybeetle densities were significantly positively related toaphid densities at the local field and landscapescales, hinting at the absence of differential specialization to aphid densitiesby different ladybeetle species. However, we found that H. axyridis was the dominant ladybeetle species at the local field scale, whereas P. japonica was the dominant ladybeetle species at the landscape scale. This may be because we only surveyed two crops (i.e., cotton and maize) at the landscape scale, while there are 96 common plant species(crops, trees, and weeds) at the local field scale, which increase the diversity of host plant of their prey, especially H. axyridis more preferred to settle on trees compared with crop plants. Meanwhile, P. japonica maybe more competitive and more rapidly builds populations on cotton and maize crops, although H. axyridis has been found to be a strong intraguild competitor [2,84]. In Chinese agroecosystems, H. axyridis has greater predatory capacity on eggs of C. septempunctata and P. japonica in laboratory choice and no-choice trials, and all predator–prey combinations in the field yieldedpositive documentation of intraguild predation by using DNA-based gut-content detection [85]. Additionally, the two species may have slightly shifted phenologies [86], allowing decreased interspecific competition between them.

5. Conclusions

In summary, we reported that aphid density was a key determinant of predatory ladybeetle abundances, and plant taxa also to a certain extent affected ladybeetle densities. Host plant species likely influenced ladybeetle foraging behavior through temporal host shifts. Hence, it would be helpful to better characterize host plant reservoirs and temporal host shifts in more complex landscapes, as well as the spatial cues related to food resources, since non-crop habitats are also heavily used by predatory ladybeetles. By improving non-crop and crop habitats in synchrony with crop development and aphid population growth, it would help to maximize predatory ladybeetle movements from non-crop plants into crop fields, there by promote conservation biocontrol of aphid populations by predatory ladybeetles.

Supplementary Materials

The following are available online at https://www.mdpi.com/2075-4450/11/10/695/s1. Additional Supplementary Information may be found in the online version of this article: Table S1: Plant species and plant functional groups surveyed at the local field scale in 2013 and 2014. Table S2: Effects of aphids and plant taxa on ladybeetle metapopulation density (LMD) at the local field scale. Figure S1: The distribution of 83 study sites at the regional landscape scale in northern China in 2013. Figure S2: The species composition of aphidophagous ladybeetles at the local field and landscape scales. Figure S3: The GLMM analysis results for the effects of aphid density on ladybeetle metapopulation density in different plant functional groups at the local field scale. Figure S4: Dominant host plant species of ladybeetles and corresponding aphid densities at the local farmland scale from the surveys in 2013 and 2014.

Author Contributions

Conceptualization, H.P. and Y.L. (Yanhui Lu); data curation, H.P., B.L., and Y.L. (Yanhui Lu); formal analysis, H.P., B.L., N.D., C.C.J., and Y.L. (Yanhui Lu); funding acquisition, Y.L. (Yanhui Lu); investigation, H.P., B.L., L.Y., and Y.L. (Yongqiang Liu); methodology, Y.L. (Yanhui Lu); project administration, Y.L. (Yanhui Lu); resources, Y.L. (Yanhui Lu); software, Y.L. (Yanhui Lu); supervision, Y.L. (Yanhui Lu); validation, H.P., B.L., and Y.L. (Yanhui Lu); visualization, H.P., B.L., and Y.L. (Yanhui Lu);writing—original draft, H.P., B.L., and C.C.J.; writing—review & editing, E.T., N.D., and Y.L. (Yanhui Lu). All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Funds of China (Nos. 31621064, 31222046), Tianshan Innovation Team Plan (No. 2020D14004), and by the H2020 project EUCLID project (No. 633999).

Acknowledgments

The authors thank the graduate trainees who were at the Langfang Experimental Station, CAAS, from 2013 to 2014, for their assistance with the field surveys.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hodek, I.; Honĕk, A.; van Emden, H.F. Ecology and Behaviour of the Ladybeetle Beetles (Coccinellidae); John Wiley & Sons: Hoboken, NJ, USA, 2012. [Google Scholar]
  2. Dixon, A.F.G. Insect Predator-Prey Dynamics: Ladybeetle Beetles and Biological Control; Cambridge University Press: Cambridge, UK, 2000. [Google Scholar]
  3. Obrycki, J.J.; Harwood, J.D.; Kring, T.J.; O’Neil, R.J. Aphidophagy by Coccinellidae: Application of biological control in agroecosystems. Biol. Control 2009, 51, 244–254. [Google Scholar] [CrossRef]
  4. Fox, T.B.; Landis, D.A.; Cardoso, F.F.; Difonzo, C.D. Predators suppress Aphis glycines Matsumura population growth in soybean. Environ. Entomol. 2004, 33, 608–618. [Google Scholar] [CrossRef] [Green Version]
  5. Ali, A.; Desneux, N.; Lu, Y.H.; Wu, K.M. Key aphid natural enemies showing positive effects on wheat yield through biocontrol services in northern China. Agric. Ecosyst. Environ. 2018, 226, 1–9. [Google Scholar] [CrossRef]
  6. Karp, D.S.; Chaplin-Kramer, R.; Meehan, T.D.; Martin, E.A.; DeClerck, F.; Grab, H.; Gratton, C.; Hunt, L.; Larsen, A.E.; Martínez-Salinas, A.; et al. Crop pests and predators exhibit inconsistent responses to surrounding landscape composition. Proc. Nat. Acad. Sci. USA 2018, 115, E7863–E7870. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Costamagna, A.C.; Landis, D.A.; Difonzo, C.D. Suppression of soybean aphid by generalist predators results in a trophic cascade in soybeans. Ecol. Appl. 2007, 17, 441–451. [Google Scholar] [CrossRef]
  8. Obrycki, J.J.; Kring, T.J. Predaceous Coccinellidae in biological control. Ann. Rev. Entomol. 1998, 43, 295–321. [Google Scholar] [CrossRef] [Green Version]
  9. Schellhorn, N.A.; Bianchi, F.J.J.A.; Hsu, C.L. Movement of entomophagous arthropods in agricultural landscapes: Links to pest suppression. Ann. Rev. Entomol. 2014, 59, 559–581. [Google Scholar] [CrossRef]
  10. Haddad, N.M.; Crutsinger, G.M.; Gross, K.; Haarstad, J.; Knops, J.M.H.; Tilman, D. Plant species loss decreases arthropod diversity and shifts trophic structure. Ecol. Lett. 2009, 12, 1029–1039. [Google Scholar] [CrossRef]
  11. Lundgren, J.G.; Wyckhuys, K.A.G.; Desneux, N. Population responses by Oriusinsidiosus to vegetational diversity. Biocontrol 2009, 54, 135–142. [Google Scholar] [CrossRef]
  12. Scherber, C.; Eisenhauer, N.; Weisser, W.W.; Schmid, B.; Voigt, W.; Fischer, M.; Schulze, E.D.; Roscher, C.; Weigelt, A.; Allan, E.; et al. Bottom-up effects of plant diversity on multitrophic interactions in a biodiversity experiment. Nature 2010, 468, 553–556. [Google Scholar] [CrossRef]
  13. Thomine, E.; Jeavons, E.; Rusch, A.; Bearez, P.; Desneux, N. Effect of crop diversity on predation activity and population dynamics of the mirid predator Nesidiocoris tenuis. J. Pest Sci. 2020. [Google Scholar] [CrossRef]
  14. Wäckers, F.L. Suitability of Extra-Floral Nectar, Pollen, and Honeydew as Insect Food Sources. In Plant-Provided Food for Carnivorous Insects; Wäckers, F.L., van Rijn, P.C.J., Bruin, J., Eds.; Cambridge University Press: Cambridge, UK, 2005. [Google Scholar]
  15. Wäckers, F.L.; van Rijn, P.C.J. Pick and Mix: Selecting Flowering Plants to Meet the Requirements of Target Biological Control Insects. In Biodiversity and Insect Pests: Key Issues for Sustainable Management; Gurr, S.D., Snyder, W.E., Read, D.M.Y., Eds.; John Wiley and Sons Ltd.: Chichester, UK, 2012. [Google Scholar]
  16. Thomine, E.; Rusch, A.; Supplisson, C.; Monticelli, L.S.; Amiens-Desneux, E.; Lavoir, A.V.; Desneux, N. Highly diversified crop systems can promote the dispersal and foraging activity of the generalist predator Harmonia axyridis. Entomol. Gen. 2020, 40, 133–145. [Google Scholar] [CrossRef] [Green Version]
  17. Wolf, S.; Romeis, J.; Collatz, J. Utilization of plant-derived food sources from annual flower strips by the invasive harlequin ladybird Harmonia axyridis. Biol. Control 2018, 122, 118–126. [Google Scholar] [CrossRef]
  18. Leather, S.R.; Cooke, R.C.A.; Fellowes, M.D.E.; Rombe, R. Distribution and density of ladybeetles (Coleoptera: Coccinelldiae) in non-crop habitats. Eur. J. Entomol. 1999, 96, 23–27. [Google Scholar]
  19. Prasifka, J.R.; Heinz, K.M.; Minzenmayer, R.R. Timing, magnitude, rates, and putative causes of predator movement between cotton and grain sorghum fields. Environ. Entomol. 2004, 33, 282–290. [Google Scholar] [CrossRef] [Green Version]
  20. Elliott, N.C.; Kieckhefer, R.W.; Phoofolo, M.W. Prey foraging movements by Hippodamia convergens 1 in wheat are influenced by hunger and aphids. Southwest Entomol. 2013, 38, 163–172. [Google Scholar] [CrossRef]
  21. Villegas, C.M.; Verdugo, J.A.; Grez, A.A.; Tapia, J.; Lavandero, B. Movement between crops and weeds: Temporal refuges for aphidophagous insects in Central Chile. Cienc. Investig. Agrar. 2013, 40, 317–326. [Google Scholar] [CrossRef] [Green Version]
  22. Evans, E.W.; Dixon, A.F.G. Cues for oviposition by ladybeetle beetles (Coccinellidae): Response to aphids. J. Anim. Ecol. 1986, 55, 1027–1034. [Google Scholar] [CrossRef]
  23. Minoretti, N.; Weisser, W.W. The impact of individual ladybeetles (Coccinellaseptempunctata, Coleoptera: Coccinellidae) on aphids colonies. Eur. J. Entomol. 2000, 97, 475–479. [Google Scholar] [CrossRef] [Green Version]
  24. Van der Werf, W.; Evans, E.W.; Powell, J. Measuring and modelling the dispersal of Coccinellaseptempunctata (Coleoptera: Coccinellidae) in alfalfa fields. Eur. J. Entomol. 2000, 97, 487–493. [Google Scholar] [CrossRef] [Green Version]
  25. Evans, E.W. Searching and reproductive behaviour of female aphidophagous ladybeetles (Coleoptera: Coccinellidae): A review. Eur. J. Entomol. 2003, 100, 1–10. [Google Scholar] [CrossRef]
  26. Bianchi, F.J.J.A.; van der Werf, W. Model evaluation of the function of prey in non-crop habitats for biological control by ladybeetles in agricultural landscapes. Ecol. Model 2004, 171, 177–193. [Google Scholar] [CrossRef]
  27. Gardiner, M.M.; Landis, D.A.; Gratton, C.; Schmidt, N.; O’Neal, M.; Mueller, E.; Chacon, J.; Heimpel, G.E.; DiFonzo, C.D. Landscape composition influences patterns of native and exotic lady beetle abundance. Divers. Distrib. 2009, 15, 554–564. [Google Scholar] [CrossRef]
  28. Woltz, M.J.; Landis, D.A. Coccinellid response to landscape composition and configuration. Agric. For. Entomol. 2014, 16, 341–349. [Google Scholar] [CrossRef]
  29. Schellhorn, N.A.; Parry, H.R.; Macfadyen, S.; Wang, Y.; Zalucki, M.P. Connecting scales: Achieving in-field pest control from areawide and landscape ecology studies. Insect Sci. 2015, 22, 35–51. [Google Scholar] [CrossRef]
  30. Marrec, R.; Badenhausser, I.; Bretagnolle, V.; Börger, L.; Roncoroni, M.; Guillon, N.; Gauffre, B. Crop succession and habitat preferences drive the distribution and abundance of carabid beetles in an agricultural landscape. Agric. Ecosyst. Environ. 2015, 199, 282–289. [Google Scholar] [CrossRef] [Green Version]
  31. Djoudi, E.A.; Plantegenest, M.; Aviron, S.; Pétillon, J. Local vs. landscape characteristics differentially shape emerging and circulating assemblages of carabid beetles in agroecosystems. Agric. Ecosyst. Environ. 2019, 270, 149–158. [Google Scholar] [CrossRef]
  32. Gardiner, M.; Landis, D.; Gratton, C.; DiFonzo, C.; O’neal, M.; Chacon, J.; Wayo, M.; Schmidt, N.; Mueller, E.; Heimpel, G. Landscape diversity enhances biological control of an introduced crop pest in the north-central USA. Ecol. Appl. 2009, 19, 143–154. [Google Scholar] [CrossRef]
  33. Chaplin-Kramer, R.; Kremen, C. Pest control experiments show benefits of complexity at landscape and local scales. Ecol. Appl. 2012, 22, 1936–1948. [Google Scholar] [CrossRef]
  34. Ostrom, P.H.; Colunga-Garcia, M.; Gage, S.H. Establishing pathways of energy flow for insect predators using stable isotope ratios: Field and laboratory evidence. Oecologia 1997, 109, 108–113. [Google Scholar] [CrossRef]
  35. Oelbermann, K.; Scheu, S. Stable isotope enrichment (delta N-15 and delta C-13) in a generalist predator (Pardosalugubris, Araneae: Lycosidae): Effects of prey quality. Oecologia 2002, 130, 337–344. [Google Scholar] [CrossRef]
  36. Schallhart, N.; Wallinger, C.; Juen, A.; Traugott, M. Dispersal abilities of adult click beetles in arable land revealed by analysis of carbon stable isotopes. Agric. For. Entomol. 2009, 11, 333–339. [Google Scholar] [CrossRef]
  37. Forbes, K.J.; Gratton, C. Stable isotopes reveal different patterns of inter-crop dispersal in two ladybeetle species. Ecol. Entomol. 2011, 36, 396–400. [Google Scholar] [CrossRef]
  38. Ouyang, F.; Men, X.Y.; Yang, B.; Su, J.W.; Zhang, Y.S.; Zhao, Z.H.; Ge, F. Maize benefits the predatory beetle, Propylea japonica (Thunberg), to provide potential to enhance biological control for aphids in cotton. PLoS ONE 2012, 7, e44379. [Google Scholar] [CrossRef] [PubMed]
  39. Madeira, F.; di Lascio, A.; Carlino, P.; Costantini, M.L.; Rossi, L.; Pons, X. Stable carbon and nitrogen isotope signatures to determine predatordispersal between alfalfa and maize. Biol. Control 2014, 77, 66–75. [Google Scholar] [CrossRef]
  40. Lu, Y.H.; Wu, K.M.; Jiang, Y.Y.; Guo, Y.Y.; Desneux, N. Widespread adoption of Bt cotton and insecticide decrease promotes biocontrol services. Nature 2012, 487, 362–365. [Google Scholar] [CrossRef]
  41. Ali, A.; Desneux, N.; Lu, Y.H.; Liu, B.; Wu, K.M. Characterization of the natural enemy community attacking cotton aphid in the Bt cotton ecosystem in northern China. Sci. Rep. 2016, 6, 24273. [Google Scholar] [CrossRef]
  42. Han, P.; Niu, C.Y.; Desneux, N. Identification of top-down forces regulating cotton aphid population growth in transgenic Bt cotton in central China. PLoS ONE 2014, 9, e102980. [Google Scholar] [CrossRef]
  43. Wang, Z.R. Farmland Weeds in China: A Collection of Colored Illustrative Plates; Agricultural Publishing House: Beijing, China, 1990. [Google Scholar]
  44. Pan, H.S.; Lu, Y.H.; Wyckhuys, K.A.G. Early-Season host switching in Adelphocoris spp. (Hemiptera: Miridae) of differing host breadth. PLoS ONE 2013, 8, e59000. [Google Scholar] [CrossRef] [Green Version]
  45. Ren, S.X.; Wang, X.M.; Pang, H.; Peng, Z.Q.; Zeng, T. Colored Pictorial Handbook of Ladybeetle Beetles in China; Science Press: Beijing, China, 2009. [Google Scholar]
  46. Desneux, N.; Decourtye, A.; Delpuech, J.M. The sublethal effects of pesticides on beneficial arthropods. Ann. Rev. Entomol. 2007, 52, 81–106. [Google Scholar] [CrossRef]
  47. Bolker, B.M.; Brooks, M.E.; Clark, C.J.; Geange, S.W.; Poulsen, J.R.; Stevens, M.H.H.; White, J.S.S. Generalized linear mixed models: A practical guide for ecology and evolution. Trends Ecol. Evol. 2009, 24, 127–135. [Google Scholar] [CrossRef] [PubMed]
  48. Venables, W.N.; Ripley, B.D. Modern Applied Statistics with S. Fourth Edition; Springer: New York, NY, USA, 2002. [Google Scholar]
  49. R Core Team R. A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2019. [Google Scholar]
  50. Bates, D.; Mächler, M.; Bolker, B.; Walker, S. Fitting linear mixed-effects models using lme4. J. Stat. Soft. 2015, 67, 133–199. [Google Scholar] [CrossRef]
  51. Kareiva, P.; Odell, G. Swarms of predators exhibit “preytatxis” if individual predators use area-restricted search. Am. Nat. 1987, 130, 233–270. [Google Scholar] [CrossRef]
  52. Schellhorn, N.A.; Andow, D.A. Response of coccinellids to their aphid prey at different spatial scales. Popul. Ecol. 2005, 47, 71–76. [Google Scholar] [CrossRef]
  53. Seagraves, M.P. Lady beetle oviposition behavior in response to the trophic environment. Biol. Control 2009, 51, 312–322. [Google Scholar] [CrossRef]
  54. Honĕk, A.; Martinkova, Z.; Štrobach, J. Effect of aphid abundance and urbanization on the abundance of Harmonia axyridis (Coleoptera: Coccinellidae). Eur. J. Entomol. 2018, 115, 703–707. [Google Scholar] [CrossRef]
  55. Dassou, A.G.; Tixier, P. Response of pest control by generalist predators to local-scale plant diversity: A meta-analysis. Ecol. Evol. 2016, 6, 1143–1153. [Google Scholar] [CrossRef] [Green Version]
  56. Bertrand, C.; Eckerter, P.W.; Ammann, L.; Entling, M.H.; Gobet, E.; Herzog, F.; Mestre, L.; Tinner, W.; Albrecht, M. Seasonal shifts and complementary use of pollen sources by two bees, a lacewing and a ladybeetle species in European agricultural landscapes. J. Appl. Ecol. 2019, 56, 2431–2442. [Google Scholar] [CrossRef] [Green Version]
  57. Schuldiner-Harpaz, T.; Coll, M. Estimating the effect of plant-provided food supplements on pest consumption by omnivorous predators: Lessons from two coccinellid beetles. Pest Manag. Sci. 2017, 73, 976–983. [Google Scholar] [CrossRef]
  58. Fonseca, M.M.; Lima, E.; Lemos, F.; Venzon, M.; Janssen, A. Non-crop plant to attract and conserve an aphid predator (Coleoptera: Coccinellidae) in tomato. Biol. Control 2017, 115, 129–134. [Google Scholar] [CrossRef]
  59. Kareiva, P.; Sahakian, R. Tritrophic effects of a simple architectural mutation in pea plants. Nature 1990, 345, 433–434. [Google Scholar] [CrossRef]
  60. Grevstad, F.S.; Klepetka, B.W. The influence of plant architecture on the foraging efficiencies of a suite of ladybird beetles feeding on aphids. Oecologia 1992, 92, 399–404. [Google Scholar] [CrossRef] [PubMed]
  61. Reynolds, P.G.; Cuddington, K. Effects of plant gross morphology on predator searching behaviour. Environ. Entomol. 2012, 41, 516–522. [Google Scholar] [CrossRef] [PubMed]
  62. Honĕk, A.; Martinkova, Z.; Dixon, A.F.G.; Skuhrovec, J.; Roy, H.E.; Brabec, M.; Pekar, S. Life cycle of Harmonia axyridis in central Europe. BioControl 2018, 63, 241–252. [Google Scholar] [CrossRef] [Green Version]
  63. Honĕk, A.; Martinkova, Z.; Roy, H.E.; Dixon, A.F.G.; Skuhrovec, J.; Pekár, S.; Brabec, M. Differences in thephenology of Harmonia axyridis (Coleoptera: Coccinellidae) and nativecoccinellids in Central Europe. Environ. Entomol. 2019, 48, 80–87. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Honĕk, A.; Brabec, M.; Martinkova, Z.; Dixon, A.F.G.; Pekár, S.; Skuhrovec, J. Factors determining local and seasonal variation in abundance of Harmonia axyridis in Central Europe. Eur. J. Entomol. 2019, 116, 93–103. [Google Scholar] [CrossRef]
  65. Zhu, J.; Park, K.C. Methyl salicylate, a soybean aphid-induced plant volatile attractive to the predator CoccinellaSeptempunctata. J. Chem. Ecol. 2005, 31, 1733–1746. [Google Scholar] [CrossRef]
  66. Girling, R.D.; Hassal, M. Behavioural responses of the seven-spot ladybeetle Coccinellaseptempunctata to plant headspace chemicals collected from four crop Brassicas and Arabidopsis thaliana, infested with Myzuspersicae. Agric. For. Entomol. 2008, 10, 297–306. [Google Scholar] [CrossRef]
  67. Leroy, P.D.; Schillings, T.; Farmakidis, J.; Heuskin, S.; Lognay, G.; Verheggen, F.J.; Brostaux, Y.; Haubruge, E.; Francis, F. Testing semiochemicals from aphid, plant and conspecific: Attraction of Harmonia axyridis. Insect Sci. 2012, 19, 372–382. [Google Scholar] [CrossRef]
  68. Turlings, T.; Erb, M. Tritrophic interactions mediated by herbivore-induced plant volatiles: Mechanisms, ecological relevance, and application potential. Ann. Rev. Entomol. 2018, 63, 433–452. [Google Scholar] [CrossRef]
  69. Xiu, C.L.; Xu, B.; Pan, H.S.; Zhang, W.; Yang, Y.Z.; Lu, Y.H. Volatiles from Sophora japonica flowers attract Harmonia axyridis adults (Coleoptera: Coccinellidae). J. Integr. Agr. 2019, 18, 873–883. [Google Scholar] [CrossRef] [Green Version]
  70. Risch, S.J.; Wrubel, R.; Andow, D. Foraging by a predaceous beetle, Coleomegillamaculata (Coleoptera: Coccinellidae), in a polyculture: Effects of plant density and diversity. Env. Entomol. 1982, 11, 949–950. [Google Scholar] [CrossRef] [Green Version]
  71. Andow, D.A.; Risch, S.J. Predation in diversified agroecosystems: Relations between a coccinellid predator Coleomegillamaculata and its food. J. Appl. Ecol. 1985, 2, 357–372. [Google Scholar] [CrossRef]
  72. Vermunt, B.; Cuddington, K.; Sobek-Swant, S.; Crosthwaite, J.C.; Lyons, D.B.; Sinclair, B.J. Temperatures experienced by wood-boring beetles in the under-bark microclimate. Forest Ecol. Manag. 2012, 269, 149–157. [Google Scholar] [CrossRef]
  73. Tougeron, K.; van Baaren, J.; Burel, F.; Alford, L. Comparing thermal tolerance across contrasting landscapes: First steps towards understanding how landscape management could modify ectotherm thermal tolerance. Insect Conserv. Diver. 2016, 9, 171–180. [Google Scholar] [CrossRef] [Green Version]
  74. Alford, L.; Tougeron, K.; Pierre, J.S.; Burel, F.; Van Baaren, J. The effect of landscape complexity and microclimate on the thermal tolerance of a pest insect. Insect Sci. 2017, 25, 905–915. [Google Scholar] [CrossRef] [PubMed]
  75. Grez, A.A.; Torres, C.; Zaviezo, T.; Lavander, B.; Ramírez, M. Migration of coccinellids to alfalfa fields with varying adjacent vegetation in Central Chile. Cienc. Investig. Agrar. 2010, 37, 111–121. [Google Scholar] [CrossRef] [Green Version]
  76. Pan, H.S.; Xiu, C.L.; Liu, B.; Wyckhuys, K.A.G.; Lu, Y.H. Whorl-stage maize provides a microclimate refuge for predatory ladybeetles. Biol. Control 2020, 142, 104162. [Google Scholar] [CrossRef]
  77. Wyckhuys, K.A.G.; Lu, Y.H.; Morales, H.; Vazquez, L.L.; Legaspi, J.C.; Eliopoulos, P.A.; Hernandez, L.M. Current status and potential of conservation biological control for agriculture in the developing world. Biol. Control 2013, 65, 152–167. [Google Scholar] [CrossRef]
  78. Rayl, R.J.; Shields, M.W.; Tiwari, S.; Wratten, S.D. Conservation biological control of insect pests. Sust. Agric. Rev. 2018, 28, 103–124. [Google Scholar]
  79. Brewer, M.J.; Gordy, J.W.; Kerns, D.L.; Woolley, J.B.; Rooney, W.L.; Bowling, R.D. Sugarcane aphid population growth, plant injury, and natural enemies on selected grain sorghum hybrids in Texas and Louisiana. J. Econ. Entomol. 2017, 110, 2109–2118. [Google Scholar] [CrossRef] [PubMed]
  80. Jaworski, C.C.; Xiao, D.; Xu, Q.X.; Ramirez-Romero, R.; Guo, X.J.; Wang, S.; Desneux, N. Varying the spatial arrangement of synthetic herbivore-induced plant volatiles and companion plants to improve conservation biological control. J. Appl. Ecol. 2019, 56, 1176–1188. [Google Scholar] [CrossRef]
  81. Yasuda, H.; Tagaki, T.; Kogi, K. Effects of conspecific and heterospecific larval tracks on the oviposition behaviour of the predatory ladybeetle, Harmonia axyridis (Coleoptera: Coccinellidae). Eur. J. Entomol. 2000, 97, 551–553. [Google Scholar] [CrossRef] [Green Version]
  82. Hoogendoorn, M.; Heimpel, G.E. Competitive interactions between an exotic and a native ladybeetle: A field cage study. Entomol. Exp. Appl. 2004, 111, 19–28. [Google Scholar] [CrossRef]
  83. Oliver, T.H.; Timms, J.E.L.; Taylor, A. Oviposition responses to patch quality in the larch ladybeetle Aphidectaolbiterata (Coleoptera: Coccinellidae): Effects of aphid density, and con- and heterospecific tracks. B. Entomol. Res. 2006, 96, 25–34. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Ovchinnikov, A.N.; Belyakova, N.A.; Ovchinnikova, A.A.; Reznik, S.Y. Factors determining larval cannibalistic behavior in invasive and native populations of the multicolored Asian ladybird, Harmonia axyridis. Entomol. Gen. 2019, 38, 243–254. [Google Scholar] [CrossRef]
  85. Yang, F.; Wang, Q.; Wang, D.M.; Xu, B.; Xu, J.X.; Lu, Y.H.; Harwood, J.D. Intraguild predation among three common Coccinellids (Coleoptera: Coccinellidae) in China: Detection using DNA-based gut-content analysis. Environ. Entomol. 2017, 46, 1–10. [Google Scholar] [CrossRef]
  86. Honěk, A.; Dixon, A.F.G.; Soares, A.O.; Skuhrovec, J.; Martinkova, Z. Spatial and temporal changes in the abundance and 16omposition of ladybird (Coleoptera: Coccinellidae) communities. Curr. Opin. Insect Sci. 2017, 20, 61–67. [Google Scholar] [CrossRef]
Insects 11 00695 g001 550
Figure 1. Population dynamics of aphids and ladybeetles on different plant functional groupsatthe field scale in 2013 and 2014. Aphid density (AD) in 2013 (a) and 2014 (b) and ladybeetle metapopulation density (LMD) in 2013 (c) and 2014 (d) are the mean values on each sampling date in three plant functional groups. Boxplots nested in the top right corner are the mean values of AD and LMD throughout the whole sampling period. Diamond points with black color in the box are the mean value for each functional plant taxon, whereas the line in each box is the median, and box edges represent the lower (Q1) and upper (Q3) quartiles. The number of plant species of plant functional groups (crops, trees, weeds) was 21, 8, and 49 species in 2013 and 26, 8, and 46 species in 2014, respectively. The p value indicated no significant difference (p > 0.05) between these plant taxa within the Kruskal–Wallis test.
Figure 1. Population dynamics of aphids and ladybeetles on different plant functional groupsatthe field scale in 2013 and 2014. Aphid density (AD) in 2013 (a) and 2014 (b) and ladybeetle metapopulation density (LMD) in 2013 (c) and 2014 (d) are the mean values on each sampling date in three plant functional groups. Boxplots nested in the top right corner are the mean values of AD and LMD throughout the whole sampling period. Diamond points with black color in the box are the mean value for each functional plant taxon, whereas the line in each box is the median, and box edges represent the lower (Q1) and upper (Q3) quartiles. The number of plant species of plant functional groups (crops, trees, weeds) was 21, 8, and 49 species in 2013 and 26, 8, and 46 species in 2014, respectively. The p value indicated no significant difference (p > 0.05) between these plant taxa within the Kruskal–Wallis test.
Insects 11 00695 g001
Insects 11 00695 g002 550
Figure 2. The relationship between aphid density (AD) and ladybeetle metapopulation density (LMD) atthe local field scale in 2013 and 2014. The statistical results of a negative binomial generalized linear mixed-effects model (GLMM model using the “glmer.nb” function in the “MASS” package of the R software) are shown in Table S2. Because theeffects of plant functional groups were not significant, we only provided the aphid effects here. The data points (with marginal rug) of different colors are for different plant functional groups, and the red smooth curve with 95% confidence interval band (shaded areas) is across all plant species in 2013 (a) and 2014 (b) due to no difference in population density in these functional plant taxa. Data on AD shown in the scatter plotswere transformed by log10(x + 1) before analysis.
Figure 2. The relationship between aphid density (AD) and ladybeetle metapopulation density (LMD) atthe local field scale in 2013 and 2014. The statistical results of a negative binomial generalized linear mixed-effects model (GLMM model using the “glmer.nb” function in the “MASS” package of the R software) are shown in Table S2. Because theeffects of plant functional groups were not significant, we only provided the aphid effects here. The data points (with marginal rug) of different colors are for different plant functional groups, and the red smooth curve with 95% confidence interval band (shaded areas) is across all plant species in 2013 (a) and 2014 (b) due to no difference in population density in these functional plant taxa. Data on AD shown in the scatter plotswere transformed by log10(x + 1) before analysis.
Insects 11 00695 g002
Insects 11 00695 g003 550
Figure 3. The population density of aphids and ladybeetles in different sampling periods (months) on cotton and maize crops at the regional landscape scale in 2013. (a): aphid density (AD); (b): ladybeetle metapopulation density (LMD). The data are the average ± S.E. with log10(x + 1) transformed before using repeated-measures analysis with a mixed effect linear model (“proc mixed” in the SAS software). Asterisk (*) indicates the significant difference between cotton and maize in different sampling months, * p < 0.05, *** p < 0.001.
Figure 3. The population density of aphids and ladybeetles in different sampling periods (months) on cotton and maize crops at the regional landscape scale in 2013. (a): aphid density (AD); (b): ladybeetle metapopulation density (LMD). The data are the average ± S.E. with log10(x + 1) transformed before using repeated-measures analysis with a mixed effect linear model (“proc mixed” in the SAS software). Asterisk (*) indicates the significant difference between cotton and maize in different sampling months, * p < 0.05, *** p < 0.001.
Insects 11 00695 g003
Insects 11 00695 g004 550
Figure 4. The relationship between aphid density (AD) and ladybeetle metapopulation density (LMD) in different sampling periods (months) on cotton and maize crops at the regional landscape scale in 2013. The data are transformed by log10(x + 1) before using GLMM analysis with Poisson distribution (“lme4” package of the R software). The scatter plot with different colors shows the data points (with marginal rug) and the regression trend lines (smooth curves) and 95% confidence interval band (shaded areas) in different months (June: steel blue color; July: orange color; August: pale green color). The smooth regression curve in red indicates the trend across the entire sampling period (from June to August). Circle points and triangle points indicate the data from cotton (a) and maize (b), respectively. The statistical results of GLMMs are shown in Table 1.
Figure 4. The relationship between aphid density (AD) and ladybeetle metapopulation density (LMD) in different sampling periods (months) on cotton and maize crops at the regional landscape scale in 2013. The data are transformed by log10(x + 1) before using GLMM analysis with Poisson distribution (“lme4” package of the R software). The scatter plot with different colors shows the data points (with marginal rug) and the regression trend lines (smooth curves) and 95% confidence interval band (shaded areas) in different months (June: steel blue color; July: orange color; August: pale green color). The smooth regression curve in red indicates the trend across the entire sampling period (from June to August). Circle points and triangle points indicate the data from cotton (a) and maize (b), respectively. The statistical results of GLMMs are shown in Table 1.
Insects 11 00695 g004
Table
Table 1. Generalized linear mixed models (GLMMs) related ladybeetle metapopulation density (LMD) with aphid density (AD) on cotton and maize crops in different months at the regional landscape scale.The fixed effects included AD, crop, month, and interactions of different explanatory variables. Estimates and S.E. was the coefficient and standard error of the fixed effect. The bold p value indicates a significant effect (p < 0.05).
Table 1. Generalized linear mixed models (GLMMs) related ladybeetle metapopulation density (LMD) with aphid density (AD) on cotton and maize crops in different months at the regional landscape scale.The fixed effects included AD, crop, month, and interactions of different explanatory variables. Estimates and S.E. was the coefficient and standard error of the fixed effect. The bold p value indicates a significant effect (p < 0.05).
Fixed EffectsEstimatesS.E.Wald Zp
(Intercept)1.856000.07092026.17<0.001
AD0.000010.0000033.93< 0.001
Crop_maize−0.160500.029880−5.37<0.001
Month_August1.062000.03583029.65<0.001
AD:crop_maize0.000350.00000491.37<0.001
AD:month_July0.000180.0000316.01<0.001
AD:month_August−0.000340.000004−87.56<0.001
Table
Table 2. Estimated proportion of diet for Propylea japonica adults originating from C3 (cotton) and C4 (maize) plants at a regional landscape composed of cotton and maize in 2013. The carbon isotope ratioδ13Cdata are shown as the mean ±SE.
Table 2. Estimated proportion of diet for Propylea japonica adults originating from C3 (cotton) and C4 (maize) plants at a regional landscape composed of cotton and maize in 2013. The carbon isotope ratioδ13Cdata are shown as the mean ±SE.
CropsJuneJulyAugust
No. SamplesCarbon Isotope Ratioδ13CC3 PlantC4 PlantNo. SamplesCarbon Isotope Ratioδ13cC3 PlantC4 PlantNo. SamplesCarbon Isotope Ratioδ13cC3 PlantC4 Plant
Cotton119−26.70 ± 1.01~100.00%(119)~0.00%(0)101−25.46 ± 1.78~100.00%(101)~0.00%(0)458−21.11 ± 5.92~86.57%(396)~13.43%(62)
Maize199−26.62 ± 1.51~100.00%(199)~0.00%(0)141−23.72 ± 3.91~100.00%(141)~0.00%(0)852−14.12 ± 5.01~28.29%(241)~71.71%(611)

Share and Cite

MDPI and ACS Style

Pan, H.; Liu, B.; Jaworski, C.C.; Yang, L.; Liu, Y.; Desneux, N.; Thomine, E.; Lu, Y. Effects of Aphid Density and Plant Taxa on Predatory Ladybeetle Abundance at Field and Landscape Scales. Insects 2020, 11, 695. https://doi.org/10.3390/insects11100695

AMA Style

Pan H, Liu B, Jaworski CC, Yang L, Liu Y, Desneux N, Thomine E, Lu Y. Effects of Aphid Density and Plant Taxa on Predatory Ladybeetle Abundance at Field and Landscape Scales. Insects. 2020; 11(10):695. https://doi.org/10.3390/insects11100695

Chicago/Turabian Style

Pan, Hongsheng, Bing Liu, Coline C. Jaworski, Long Yang, Yongqiang Liu, Nicolas Desneux, Eva Thomine, and Yanhui Lu. 2020. "Effects of Aphid Density and Plant Taxa on Predatory Ladybeetle Abundance at Field and Landscape Scales" Insects 11, no. 10: 695. https://doi.org/10.3390/insects11100695

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