# Sustainable laser-based technology for insect pest control

R. Gaetani, V. Lacotte, Vincent Dufour, A. Clavel, G. Duport, Karen Gaget, Federica Calevro, P. da Silva, A. Heddi, D. Vincent, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

R. Gaetani, V. Lacotte, Vincent Dufour, A. Clavel, G. Duport, et al. Sustainable laser-based technology for insect pest control. Scientific Reports, 2021, 11 (1), pp.11068. 10.1038/s41598-021-90782-7. hal-03239008

HAL Id: hal-03239008

https://hal.science/hal-03239008

Submitted on 22 Jul 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Check for updates

OPEN

# Sustainable laser-based technology for insect pest control

R. Gaetani<fnref n="1" />, V. Lacotte², V. Dufour³, A. Clavel², G. Duport², K. Gaget², F. Calevro², P. Da Silva², A. Heddi², D. Vincent<fnref n="1" /> & B. Masenelli<fnref n="1" />✉

Aphids damage directly or indirectly cultures by feeding and spreading diseases, leading to huge economical losses. So far, only the use of pesticides can mitigate their impact, causing severe health and environmental issues. Hence, innovative eco-friendly and low-cost solutions must be promoted apart from chemical control. Here, we have investigated the use of laser radiation as a reliable solution. We have analyzed the lethal dose required to kill 90% of a population for two major pest aphid species (Acyrthosiphon pisum and Rhopalosiphum padi). We showed that irradiating insects at an early stage (one-day old nymph) is crucial to lower the lethal dose without affecting plant growth and health. The laser is mostly lethal, but it can also cause insect stunting and a reduction of survivors' fecundity. Nevertheless, we did not notice any significant visible effect on the offspring of the surviving irradiated generation. The estimated energy cost and the harmless effect of laser radiation on host plants show that this physics-based strategy can be a promising alternative to chemical pesticides.

Damage caused by insect pests on human cultures is a major plague. This impact of insects has become particularly critical in the last half century as the need to feed the rapidly growing world population has required the development of industrial agriculture and the monoculture model, which reduces crop diversity and favors pest proliferation<fnref n="1" />. Insect pest damage results in annual product losses to the agriculture industry estimated from 2 to 43% in Brazil, 25 to 43% in Asia and 20% to 90% in Africa depending on the nature of cultures²⁻⁴. Because of its efficiency and low cost, the use of chemical pesticides has become the main solution to control insect populations since the Second World War, particularly in Europe and the USA⁵. However, it is now established that these toxins, either through polluted groundwater or direct air intake, have alarming consequences to public health and the environment⁶⁻¹¹, including a significantly increased risk of brain tumors¹² and Parkinson's disease¹³. In the "Ecophyto" plan, the French government announced the need to reduce pesticide use by half by 2018. Facing the failure to meet this commitment without ecologically sound alternatives, France has postponed the objective to 2025¹⁴,¹⁵.

Recently, several alternatives have been explored to reduce pesticide use but many of them rely either on non-specific (i.e. harmful to any insect without discrimination) or toxic chemistry-based strategies¹⁶⁻¹⁸. Bio-ecological methods have legitimately gained much interest and knowledge. They represent an important improvement strategy. However, some bio-ecological methods grant only a limited yield¹⁹⁻²¹. Besides these methods, complementary strategies stemming from other scientific areas than chemistry and biology are emerging. Among these, physics-based strategies are now being contemplated. A mechanical-based solution has been suggested using high-pressure water sprays²². This method can damage host plants and thus cannot be considered viable. Another uses electrical fields to repel or destroy pests but it is not specific to targeted insect pests and hardly applicable on large areas such as fields of agricultural crops²³. Moreover, it requires high voltages (a few kV), which makes its use energy consuming and hazardous.

Few pioneering studies have been conducted looking into the control of small insects, such as cockroaches, mosquitoes, fruit flies and floor beetles, using laser systems from UV to IR²⁴⁻²⁸. A laser can deliver a collimated long-range high-power beam. It is therefore adapted to open-air operation conditions. UV light is well-known to induce DNA damage at high irradiances²⁴,²⁵. However, such a method is not suited for pesticide substitution since it requires a long time before the DNA damage can develop and lead to a global lethal effect over a large population (e.g. an irradiation of 48 h²⁴, corresponding to a global delivered energy dose of 80 J cm⁻²). Moreover, provoking mutations to the DNA of insect population is not innocuous in the long term since it could lead to genetically modified and adapted species. In contrast, the use of pulsed (~ 25 ms) green laser light (532 nm

<footnote n="1">¹INSA Lyon, CNRS, Ecole Centrale de Lyon, Université Claude Bernard Lyon 1, CPE Lyon, INL, UMR5270, Univ Lyon, 69621 Villeurbanne, France. ²INSA Lyon, INRAE, BF2I, UMR 203, Univ Lyon, 69621 Villeurbanne, France. ³INSA Lyon, Ecole Centrale de Lyon, CNRS, Université Claude Bernard Lyon 1, AMPERE, UMR5005, Univ Lyon, 69621 Villeurbanne, France. ✉email: bruno.masenelli@insa-lyon.fr</footnote>
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Aphids life cycles and morphology. (a) Monoecious holocyclic cycle of aphids *A. pisum* LL01 and YR2 lines³². (b) Dioecious holocyclic cycle of the aphid *R. padi*³³. (c) Morphology of aphids *A. pisum* LL01 (green), YR2 (pink) and *R. padi* (black) at stages N1 (1d nymph) and A15 (15d adult). *R. padi* are 2 to 3 times smaller than *A. pisum* individuals. Details on aphid cycles are provided in the Supplementary Information section. Pictures taken by Virginie Lacotte (BF2i laboratory) except for sexual forms of the aphids and their eggs³⁴,³⁵ (from the EncyclopAphid website).</figcaption></figure>

wavelength) or CO₂ IR (10.6 µm wavelength) lasers turned out to be effective to kill mosquitoes (*Anopheles stephensi*) with moderate energy doses (about 1 J cm⁻²)²⁵.

Here, we have investigated the use of laser radiation for insect control as a part of a global program that includes an initial step to specifically detect insect pests by their spectral profiles. This strategy has been proposed and patented recently²⁹, and is attracting interest on an European scale³⁰,³¹. For this strategy to be ecologically sustainable, two criteria must be met. First, the chosen laser radiation must be efficiently absorbed by the pests, leading to a lethal effect for minimal supplied energy. Second, the energy must be delivered on a short-time scale (< 100 ms) so that a wide portion of a crop field can be treated in a reasonable time (one day).

In the present study, we have focused on three distinct laser radiations; at 532 nm, 1070 nm and 10.6 µm. This is meant to probe the efficiency of two distinct mechanisms, namely absorption by the insect exoskeleton pigments and by the cell water, respectively. The first mechanism is species-dependent while the second is not. The goal was to determine the minimum energy required to kill 90% of a pest population in the minimum time. Aphids were chosen as insect pests as they are the most widespread and present a dangerous threat to cultures. Apart from being directly and indirectly extremely harmful to crops by feeding and spreading diseases, aphids are also a model system representative of other sap-sucking insect pests, including leafhoppers, mealybugs and whiteflies. We used three distinct lines belonging to two different aphid species: *Acyrthosiphon pisum* LL01 (green pea aphids), *Acyrthosiphon pisum* YR2 (pink pea aphids) and *Rhopalosiphum padi* (dark-green bird cherry-oat aphids). They offered the advantage of testing the laser on insects with different colors (and presumably with different exoskeleton properties) and of different sizes. Their life cycles are extremely complex and alternate, in the field, sexual and clonal parthenogenetic reproductions (Fig. 1), but all can be reared under obligate parthenogenesis in laboratory conditions. To be efficient and less energy consuming, any annihilation strategy thus requires treating the plants as early as possible and specifically in the early aphid development before they start reproducing clonally. Obviously, as mentioned before, a detection system that identifies the presence of pests in crop fields as soon as possible must be exploited before a selective lethal treatment is dispensed.

## Results and discussion

To assess a global strategy using a laser, the lethal dose to ensure the annihilation of 90% of the population (LD 90) after one day was first determined. To accurately retrieve this value, non-synchronized adult aphid (A. *pisum* LL01, A. *pisum* YR2 and *R. padi*) lines were used. The choice of the non-synchronized population is justified as being representative of what would be found in a field if no early detection system was coupled to the lethal system. Aphids were exposed to three distinct laser radiations, namely in the visible (532 nm wavelength), in the near infra-red (1070 nm) and in the far infra-red (10.6 µm). Based on the results obtained in this first experiment, we next focused on synchronized one-day-old nymph (N1) populations irradiated by a CO₂ laser radiation of 10.6 µm wavelength. The corresponding LD 90 was measured at different times after irradiation. The effects on the irradiated aphids, as well as on their progeny, were investigated. This allowed us to establish the efficiency and energy effectiveness of the proposed strategy.

**Lethal dose (LD90) of non-synchronized adult aphids for visible and infrared wavelengths.** Mortality rates measured one day after irradiation as a function of the fluence at 532 nm and 10.6 µm for samples of 48 aphids (see “Methods” section) are displayed in Fig. 2. Curve fitting was achieved using a logistic regression. All of the statistical tests conducted showed that the logistic regression is the good
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Mortality rates as a function of the fluence for 532 nm and 10.6 µm lasers at Day +1. Error bars are plus and minus exact binomial 95% confidence intervals. Control samples showed that 96% or more aphids survived. (a) *A. pisum* LL01 (green aphids). (b) *A. pisum* YR2 (pink aphids). (c) *R. padi* (dark-green aphids). Blue dots and Gray squares refer to experiments done using the 532 nm and 10.6 µm lasers respectively. The purple and red solid lines are logistic regressions to the data obtained with the 532 nm and 10.6 µm lasers respectively. Statistical tests are available in Supplementary Table S1.</figcaption></figure>

model to fit the data (Supplementary Table S1). Lethal doses for 90% death rate in the populations (LD 90 values) extracted from logistic fits (Fig. 2) for the three aphid lines are gathered in Supplementary Table S2.

Both the green and CO₂ lasers are efficient at killing the three aphid lines, with LD90 ranging from 1 to more than 50 J cm⁻². However, the CO₂ laser reduces the LD90 value by about 75% for *A. pisum* LL01 and *A. pisum* YR2, whereas it is nearly equivalent for *R. padi* with both lasers (reduction by only 21%). In addition, by selecting one wavelength, we noticed that *A. pisum* LL01 and *A. pisum* YR2 exhibit similar behavior, while *R. padi* remains the most sensitive to radiation, showing the smallest LD90. While *A. pisum* LL01 and YR2 both have the obligatory symbiont *Buchnera aphidicola* which participates in their nutrition, only YR2 has a secondary symbiont *Regiella insecticola* that can provide the host with more resistance to biotic and abiotic stresses (temperature, parasitoid attack, etc.)³⁶. The similar results of *A. pisum* LL01 and YR2 aphids indicate that despite the presence of the secondary symbiont, aphids are equally sensitive to the new stress caused by the CO₂ laser. Otherwise, the
high sensitivity of *R. padi* outlines a better absorption in the visible and far-infrared range. Besides, the 1070 nm wavelength is inefficient for killing *R. padi* aphids, requiring lethal doses four orders of magnitude higher than the two other wavelengths (Supplementary Table S3 and Supplementary Discussion).

It was reported on other insects that laser induced damages are due to their optical absorption properties²⁵. This implies that a better absorption coefficient leads to a higher temperature in the sample. The fundamental effect of lethality would thus be a rapid increase of the body temperature. Considering water contained in the cells to be the main cause of laser absorption, we can state that this is inefficient at 532 nm, somewhat efficient at 1070 nm, and very efficient at 10.6 µm. The absorption coefficient is a few $10^{-3}$ cm⁻¹ at 532 nm, 10 cm⁻¹ at 1070 nm and 1000 cm⁻¹ at 10.6 µm³⁷. Considering a 1 mm thick aphid (typical size of an *R. padi* specimen), much less than 1% of the 532 nm laser radiation entering the aphid is absorbed by the cells. Nonetheless, this fraction amounts to 60% and 100% at 1070 nm and 10.6 µm, respectively. Moreover, in the latter case, the radiation is absorbed over a typical depth of 50 µm. These estimations do not take into account the different reflectance for each wavelength due to the exoskeleton. Nevertheless, they partly explain why the 10.6 µm radiation is the most efficient.

If absorption by water is most likely the lethal cause at 10.6 µm, it is obviously not the case at 532 nm and 1070 nm. At both these wavelengths, absorption by aphid pigments is a more probable mechanism. Besides, if we consider the reflectance of YR2 and LL01 *A. pisum* aphids as obtained from hyperspectral measurements (Supplementary Fig. S1), we note that reflectance at 532 nm is about 0.35 or 0.4 for aphids, depending on their color. For wavelengths over 1000 nm, the reflectance is twice as large (about 0.75). This means that only half as much of the incident energy is available for pigment absorption at 1070 nm as compared to 532 nm. This is in qualitative agreement with the observed inefficiency of the 1070 nm laser to destroy aphids.

So far, we did not consider the time of energy delivery as a parameter. However, using ultra-short pulses with identical energy, very large peak powers may result in significant changes, as has already been proven at smaller length scales (scale of bacteria) ³⁸. When dealing with pigment absorption, it might be beneficial. On the contrary, when the absorption occurs in water, as the heat penetration formula is proportional to the square root of the time $t$, reducing $t$ by a factor of two results in a decrease in penetration depth of more than 30%. It will consequently be less harmful.

**Mortality rate evolution and transgenerational effects in synchronized green aphids.** *N1 mortality rate evolution.* LD90 values of about 15 J cm⁻² previously found for non-synchronized populations of *A. pisum* might seem already satisfying. However, in order to make the pest control system faster, more energy efficient and secure, any reduction of the minimum lethal dose has to be sought. This can be expected if one considers a synchronized population. In addition, if we relax the constraint of killing 90% of a population after only one day and impose the same mortality rate after a longer period, the minimum lethal dose could be lowered even further.

N1 nymphs were selected because they are the youngest and thus presumably the weakest population among the synchronized ones. The *A. pisum* LL01 nymphs (green aphids) were chosen as representative aphids because they are as sensitive as *A. pisum* YR2, and less sensitive than *R. padi*. Furthermore, the life history traits of the parthenogenetic morphs of this aphid line in the laboratory conditions used in this work are well documented in the literature³⁹. This will help to avoid any under-estimation of the LD90 for N1 populations. We considered the LD90 dose seven days after irradiation (Day + 7-timescale), and compared it to the LD90 dose one day after irradiation (Day + 1-timescale). This is justified by the development time of aphid nymphs.

Figure 3a displays the mortality rate as a function of the fluence using the CO₂ laser at Day + 1 and Day + 7. The retrieved LD90 value at Day + 1 is four times lower ($3.29$ J cm⁻²), compared to the value found for adults of the same species (12.91 J cm⁻²). Moreover, this same value is lowered by a supplemental factor of 2.6 when going from a Day + 1 to a Day + 7-timescale, reaching the minimum value of $1.15$ J·cm⁻². Finally, going from adults at Day + 1 to N1 at Day + 7 lowers the LD90 value by more than ten times. Nymphs appear to be more sensitive than adults, so treating the aphids at this development stage and at Day + 7 is of great interest if one wants to build an energy-efficient pest control system based on this method.

Moreover, a closer analysis of the mortality rate evolution (Fig. 3b) revealed that irradiation effects are not linear. First, we can roughly distinguish three categories of samples. Those irradiated with a fluence less than or equal to $0.57$ J cm⁻² are hardly affected and their mortality rate remains insignificant. For the samples irradiated with a fluence between 0.57 and $0.95$ J cm⁻², the laser radiation produces significant defects but is not enough to achieve a mortality rate of 90%. Over a fluence of $0.95$ J cm⁻², the 90% mortality rate is achieved. Second, for the latter samples, a clear threshold is observed at Day + 3 timescale. Three days after the irradiation, the change of the mortality rate is minor. This shows that, on the one hand, the minimal energy necessary to destroy N1 *A. pisum* LL01 is $1.15$ J cm⁻² and, on the other hand, the shortest time after which one can perform a second treatment is after a delay of three days.

**F0 and F1 generation effects.** Once the optimal conditions for total pest destruction have been characterized, we have wondered whether a non-lethal dose can affect the physiology of aphids and influence their fecundity and intrinsic rate of population increase, as several biotic and abiotic stresses can influence these parameters⁴⁰. Even in the case where irradiated aphids can reproduce and give birth to nymphs of second generation (F1), the laser irradiation could also affect this F1 generation. To this end, we have examined biological parameters indicative of the health status of aphids (detailed in the “Methods” section) on F0 and F1 aphid generations (born the 1st, 5th, 10th and 15th day during the F0 reproduction stage) and on their offspring (Fig. 4).

A LD50 fluence of $0.76$ J·cm⁻² corresponding to approximately 50% of mortality (Fig. 3b) was selected. Our assumption was that 50% of N1 population surviving the chosen dose could be affected in such a way that their
<figure><img src="image_6.png" /><figcaption>**Figure 3.** (a) Mortality rates as a function of the fluence using a 10.6 µm laser at Day + 1 (Blue triangles) and Day + 7 (Red squares). Control samples showed that 96% or more aphids survived. Error bars are plus and minus exact binomial 95% confidence intervals. The gray and red solid lines are logistic regressions obtained at Day + 1 and Day + 7 respectively. Statistical tests are available in Supplementary Table S4. (b) Evolution of the mortality rate over seven days at different fluences using the CO₂ laser on synchronized N1 *A. pisum* LL01 nymphs (green aphids). Letters in the right refer to the results of the 95% confidence interval GLM Gaussian test.</figcaption></figure>

development could be impaired and their offspring reduced. We observed that, when irradiated at that fluence, N1 individuals could still generate nymphs via parthenogenesis. Only the F0 nymphs that received the treatment showed significant mortality while their F1 offspring exhibited no mortality (Fig. 4a). At day seven of development (end of the nymphal development), a significant weight reduction was observed for the irradiated sample (F0 L). However, their offspring F1 showed no difference with the control, attesting that fluence of 0.76 J cm⁻² does not influence offspring's weight.

The reduction of weight on the irradiated F0 aphids comes with a developmental delay (Fig. 4c). While all the untreated individuals of the first generation (F0 C as Control) reach the adult stage, only 63% do so in the irradiated first-generation population (F0 L as Laser), the rest remaining at the end of the nymphal period (N4). The development of the F1 L nymphs is no different from the control, except those laid on D1 and D5 of which 6 and 14% were not yet adults. It is possible that the first offspring of irradiated F0 aphids have a slight delay in development.

Finally, even if the irradiated F0 aphids can still reproduce, their fecundity is also impacted (Fig. 4d and Supplementary Fig. S2). Again, only the irradiated (F0 L) population is affected with a slight decrease of fecundity and intrinsic rate of increase ($r_m$). This reduced fecundity is independent of the stunting observed previously since all the F0 aphids chosen to lay the F1 generation were adults at D0. However, the F0 L adults were smaller than the control adults (30% reduced weight), which could influence their production of embryos. The cumulated fecundity for the irradiated population (F0 L) after fifteen days is 11% lower than the control (F0 C) and its intrinsic rate of increase is 0.327 compared with 0.348 for the control.
<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 4.** (a) Mortality rates over time for the first (F0) and second (F1) aphid generations. F0 C refers to the untreated first generation, while F0 L refers to the irradiated first generation. F1 JX refers to the second generation born the 1st, 5th, 10th and 15th day during the F0 reproduction period (95% interval confidence GLM Gaussian test). (b) Aphids' mass (mg) at day seven (95% interval confidence ANOVA test; sample size F1J15 too small for statistics). (c) Proportion of adults (Ad) aphids in samples at day seven. N4 refers to the fourth stage of nymph, which precedes the adult stage (95% interval confidence χ² and Fisher tests). (d) Evolution of cumulated fecundity (95% interval confidence GLM Poisson test).</figcaption></figure>
In view of these results, we can theorize that laser treatment only influences irradiated aphid populations and not the next generations. Indeed, their offspring only shows a slight delay in reaching their maturity and their survival is not compromised.

**Global strategy against pests.** In order to be a reliable solution, a laser-based method must meet several criteria. First, the laser radiation must be efficient to kill the targeted insect. We assess this in "Introducution" and "Results and discussion" sections. The use of CO₂ laser radiation proves to be effective. The underlying mechanism, based on heating water contained in the cells, is a generic one. Indeed, prior research has found a LD90 value for mosquitoes of about 1 J cm⁻² therefore quite close to ours using the same wavelength²⁵. Targeting the water of the cells is more robust than targeting pigments of insects.

Second, the laser-based destruction of pests must be energy efficient. In order to estimate the actual energy delivered by the laser while shooting, we have considered the average area for a French sugar beet field of 17.6 ha with a row every 50 cm and a plant every 25 cm representing 13.6 million of plants per field⁴¹. Furthermore, we suppose that the early detection of N1 aphids allows one to treat only 5% of plants with one hundred aphids killed per plant. Based on the aforementioned LD90 values, this gives a total delivered energy of 200Wh per field. Obviously, this is a minor part of the energy cost estimation. Actually, the laser consumes 200 W while shooting and 150 W in stand-by. If, as assumed, the early detection allows to treat a small fraction of the field, the laser will be in stand-by mode most of the time. Consequently, the energy used will ultimately depend on the velocity of the chosen mobility vector. This point is out of the scope of the present study. However, it has been addressed first in a patent²⁹ held by the GreenShield Technology start up and more recently investigated by the consortium of academic and industrial laboratories in the framework of the EU project WeLASER³¹. Both studies highlight that a strategy combining a sensitive early detection method and a subsequent laser-based destruction technique can be economically and energetically sustainable.

Third, the strategy has to be ecologically sustainable. It must be robust with respect to a false positive event. We must thus ensure that laser radiation hitting the plant host and not an aphid does not damage the plant so as to threaten its life or growth. For that purpose, we have shot several plant organs of *Vicia faba* and *Triticum aestivum* that are our aphid host plants, using LD90 values for both aphid species at different stages (N1 or Adult) and days (D + 1 or D + 7). We used samples of 20 plants per targeted plant organ and per fluence. The results show that the fluences we use are harmless for their growth (Methods, Supplementary Discussion, and Supplementary Fig. S3, S4 and S5).

However, since the CO₂ laser acts on the water contained in the cells, the fluences used are probably harmful to non-target insect community (like aphid predators for instance). To avoid this collateral damage effect, the hyperspectral identification of insects present in the fields must be as accurate as possible, prior to the laser treatment.

## Conclusion

We have analyzed the LD90 for three aphid species using a visible and an IR laser. The CO₂ IR laser, giving rise to the lowest LD90, proves to be the most suited. Its mechanism of interaction with living matter makes it also adapted to all pests. Our study also shows that the minimum LD90 is obtained first by irradiating nymphs instead of adults and second by considering a lethal date seven days after the irradiation instead of one day. However, the laser is only effective on directly irradiated aphids and not on the offspring issued from the irradiated generation. Hence, the laser is an effective solution against aphid proliferation, particularly when coupled with an early treatment of aphids. We also have estimated the energy cost and showed that the laser effect on host plants in case of a false positive detection is negligible. Thus, a global strategy based on a laser destruction of pests can be considered sound.

## Methods

**Insects and plants.** *Acythosiphon pisum* individuals are from parthenogenetic females from two field-originated lines, LL01 and YR2 (Supplementary Fig. S6). The two lines were maintained on young broad bean plants (*Vicia faba*, L. cv. Aquadulce), in a phytotron at 21 °C ± 1 °C and 65% RH ± 5% with a 16 h photoperiod. Sampling of synchronized N1 aphids was performed as described in³⁹.

*Rhopalosiphum padi* individuals LR01 (Supplementary Fig. S7) have been reared on wheat seedling (*Triticum aestivum* [Linnaeus] cv. Orvantis) in the same controlled conditions than *A. pisum*.

In order to get close to field conditions, non-synchronized parthenogenetic adults were first used for the radiation experiments.

Organic seeds of *Vicia faba*, L. cv. Aquadulce and *Triticum aestivum* [Linnaeus] cv. Orvantis were supplied by Graines-Voltz® (Loire-Authion, France) and Nature et Progrès® (Alès, France), respectively. The seeds were germinated in commercial peat substrate Florentaise® TRH400 (Saint-Mars-Du-Désert, France) in a phytotron using the same controlled conditions as for aphid rearing. The collection of plant material in this manuscript complies with relevant institutional, national and international guidelines and can be obtained from various seed companies worldwide.

**Optical set up.** The experiment used three different lasers. A 2.4-mm-beam-diameter CO₂ pulsed laser manufactured by Access Laser, a 1-mm-beam-diameter DPSS CW 532 nm laser manufactured by CNI Lasers and an 80-µm-beam-diameter pulsed laser at 1070 nm manufactured by IPG. A red pointer was associated with them in order to aim at the targets more easily, especially for CO₂ laser pulses. Laser beams were brought to aphids using a two mirrors-movable system operated by a Jetson GPU from NVIDIA linked to an Arduino Uno card to control laser triggering (Supplementary Fig. S7). The command system triggered a time-regulated laser
impulsion through the optical system that reached the aphids' abdomen. Upstream, we monitored command
signals on an oscilloscope to check its shape and frequency.

Apart from the nature of the aphids (*A. pisum* LL01, *A. pisum* YR2 and *R. padi*), the variables were laser flu-
ence controlled by changing the impulsion time, the Pulse Width Modulation (PWM) for the CW laser (through
external command), the number of shots for pulsed lasers. The 532 nm laser cannot reach impulsion time shorter
than 5 ms due to its technical characteristics making the beam inconsistent under this value. Changing the PWM
allowed exploring the fluence under this impulsion time. An assumption made here is that changing the PWM is
equivalent to reducing the laser impulsion time and that only the average power matters. We choose impulsion
time and frequency for the CO₂ lasers as 400 µs and 750 Hz, respectively, in order to have the maximum average
power. These values were not changed during experiments each time using this laser. Thus, the overall energy
deposited by this laser was monitored through the number of delivered pulses.

Laser beam power is measured a posteriori thanks to S425-C sensor and the PM100D console manufactured by ThorLabs with a more or less 5% error announced by the manufacturer⁴².

A FX10 hyperspectral camera from Specim® have been used to measure the aphids light reflectance spectrum between 400 and 1 000 nm. Hyperspectral measurements were made on 15 days-old adults of the three aphid lines and repeated three times.

**Mortality counts.** Aphids were collected from plants, treated with the laser and then put into a translucent box with leaves changed after one day to avoid death by food insufficiency or directly on *Vicia faba* plant depending of the duration of the monitoring. We observed each aphid after one day or every day for longer period, depending on the experiment timescale (see Supplementary Fig. S9). Examinations allowed aphids discrimination in three categories: alive, sick and dead. Examination was tricky because aphids often retracted their feet and immobilized themselves when stressed. In addition, under the heat generated by the laser beam, the abdomen could darken, probably because of melanization. It does not necessarily mean that the aphid is dead. Hence, examinations must be done very meticulously, especially when performed with the smallest aphids (N1): when N1 aphids looked sick, their status were reexamined one day latter to ensure whether the insect was dead or still alive and able to recover. Pushing forward aphids by pressing on their abdomen to make them move or stack two or three aphids were useful maneuvers because they resulted in a prompt or quasi-prompt reaction from alive or sick aphids. As long as an aphid moved, even only slightly its feet, it was considered alive or sick depending on its behavior. Alive aphids can move quite quickly and can overturn when set on their back. On the contrary, sick aphids barely move and hardly or cannot overturn properly when set on their back. Furthermore, dead aphids are those with a much-darkened abdomen, spread stiff feet and do not react to any stimulus.

**N1 *A. pisum* LL01 population study protocol.** This protocol aimed at retrieving the energy inhibiting the most aphid development using the minimum energy. To this end, we first have irradiated samples of 48 N1 *A. pisum* LL01 nymphs with the CO₂ laser. Fluences applied came from those applied on *R. padi* adults because of their size similitude. For each sample, we sowed *V. faba* eleven days before depositing aphids. Therefore, these were big enough to receive the sample. *A. pisum* specimens were deposited on plants after laser irradiation. We reviewed mortality daily for seven days, using the same method described in "Mortality counts" section. This seven-day supervision resulted in the time evolution of the mortality rate (Fig. 3b).

Secondly, we have investigated laser-induced effects on first- and second-generation aphids during a biologi-
cal monitoring (Supplementary Fig. S8 and S9). For each modality (F0 C, F0 L ...), the samples were made of 50
synchronized N1 *A. pisum* LL01 nymphs splitted into 5 repetitions of 10 aphids per cage containing a plant. We
chose to apply a 50%-lethality fluence (LD50). Deposited N1 nymphs achieved their development cycle (seven
days) and reproduction cycle (fifteen days) during the biological monitoring (Fig. 1). Mortality was assessed
every day. After seven days or more, aphids reached their adult stage depending on laser damages. We, therefore,
measured their mass and count developmentally delayed aphids. *V. faba* were changed once they had reached
their maximum development stage.

Then, one adult aphid per repetition corresponding to five aphids among the samples of fifty aphids were randomly taken and deposited on a separated plant. Thus, we measured every day during fifteen days their fecundity by counting the number of offspring and their health state (mortality and deformity). The nymphs laid at day one (F1J1), five (F1J5), ten (F1J10) and fifteen (F1J15) were again deposited on a new plant. As before, we monitored their development during seven days followed by the fecundity of five randomly chosen aphids during fifteen days. We also calculated the intrinsic rate of increase of aphids for every treatment using the equation r_m = 0.74 (ln M_d)/T, with T as the time from birth to the start of reproduction, M_d as the reproductive output per aphid during a period equal to T, and 0.74 as a correction factor⁴⁰ (Supplementary Fig. S2). Hence, the biological monitoring was performed on two generations (Fig. 4).

**Laser induced effects on infested plants protocol.** This protocol aimed at reproducing a shot irra-
diating directly the plant due to a false positive detection or a shot imprecision. We monitored *Vicia faba* and
*Triticum aestivum* whose are host plants for *A. pisum* and *R. padi* respectively during seven days. Each plant was
kept in the phytotron with the same experimental conditions applied before (21 °C ± 1 °C, 65% RH ± 5%).

Thus, each sample contained 20 seedlings per fluence and per organ targeted. According to aphids on their host plants, we irradiated the adaxial (up) and abaxial (under) surfaces of the leaf, the apex and the stem. We targeted four times the same organ. To avoid the underestimation of plant damage, we have applied fluences corresponding to LD90 values for adults or N1 at 1 or 7 days (Supplementary Table S5) with the CO₂ laser.

The plants were pictured every day and biological parameters, including height, weight and leaf surface, were
measured on the seventh day (Supplementary Table S6).
**Statistical tests.** Logistic fits have been already used in similar studies²⁵. Statistic tests were done on every fit to ensure that the logistic model corresponded to the experimental results. The first one is the log-likelihood ratio and the pseudo-R². For instance, for *A. pisum* YR2 aphids irradiated with the CO₂ laser, the obtained values are 0.40 and 0.57. Moreover, the Hosmer–Lemeshow test (HLT), the area under the curve (AUC) of the Receiver Operating Characteristic (ROC) curve are other tests performed. The HLT is a statistical test for goodness of logistic fit. ROC curve shows the true positive rate over the false positive rate of the model and the closer this curve is to one, the better the fit. Using the same previous example, the HLT is passed and the AUC value is 0.91, which shows that the fit is accurate and correct.

Statistical test results for the other wavelengths and species are available (Supplementary Table S1, 4). For the 532 nm laser experiment done with adults *A. pisum* LL01, we add “Before readjustment” and “After readjustment” labels because the Hosmer–Lemeshow test was failed and the sample at 50 J cm⁻² seemed to be abnormal. This sample was removed to check whether it distorted the test. The fit was barely changed but when the HLT was done, it was used to prove that the fit was good but disturbed by this experimental sample. This is something to watch out when using the HLT because a single aberrant experimental point can make the test failed.

During the N1 green aphids (*A. pisum* LL01) population study, all the treatments were repeated five times. Statistical tests were performed using R Studio software. All statistical tests were carried out with a significance level of 5% (0.05). The daily measurements, which showed a linear change, were tested using the GLM linear regression model with Gaussian or Poisson regression after observation of the residuals (daily mortality and fecundity). The continuous quantitative measurements taken at day 7 on large samples (aphid mass, plant mass, plant height) have been evaluated with a Shapiro–Wilk normality test followed by a Levene homogeneity of variance test and an ANOVA comparison of means followed by a test of Tukey HSD. Finally, measurements relating to a proportion such as the adult aphid rate were analysed with the χ² ratio comparison test 2 to 2 or with Fisher's exact test for smaller samples.

Received: 8 March 2021; Accepted: 5 May 2021

Published online: 26 May 2021

## References

1. Andow, D. The extent of monoculture and its effects on insect pest populations with particular reference to wheat and cotton. *Agric. Ecosyst. Environ.* 9, 25–35 (1983).

2. Savary, S., Ficke, A., Aubertot, J.-N. & Hollier, C. Crop losses due to diseases and their implications for global food production losses and food security. *Food Sec.* 4, 519–537 (2012).

3. Badii, K. B., Billah, M. K., Afreh Nuamah, K., Obeng Ofori, D. & Nyarko, G. Review of the pest status, economic impact and management of fruit-infesting flies (Diptera: Tephritidae) in Africa. *Afr. J. Agric. Res.* 10, 1488–1498 (2015).

4. Oliveira, C. M., Auad, A. M., Mendes, S. M. & Frizzas, M. R. Crop losses and the economic impact of insect pests on Brazilian agriculture. *Crop Prot.* 56, 50–54 (2014).

5. Eurostat numbers from http://appsso.eurostat.ec.europa.eu/nui/submitViewTableAction.do. Accessed 9 Sept 2020.

6. Upadhayay, J., Rana, M., Juyal, V., Bisht, S. S. & Joshi, R. Impact of pesticide exposure and associated health effects. In *Pesticides in Crop Production* 69–88 (Wiley, 2020). https://doi.org/10.1002/9781119432241.ch5.

7. Kim, K.-H., Kabir, E. & Jahan, S. A. Exposure to pesticides and the associated human health effects. *Sci. Total Environ.* 575, 525–535 (2017).

8. Ye, M., Beach, J., Martin, J. W. & Senthilselvan, A. Pesticide exposures and respiratory health in general populations. *J. Environ. Sci.* 51, 361–370 (2017).

9. Fantke, P. & Jolliet, O. Life cycle human health impacts of 875 pesticides. *Int. J. Life Cycle Assess.* 21, 722–733 (2016).

10. Mrema, E. J., Ngowi, A. V., Kishinhi, S. S. & Mamuya, S. H. Pesticide exposure and health problems among female horticulture workers in Tanzania. *Environ. Health Insights* 11, 117863021771523 (2017).

11. Mahmood, I., Imadi, S. R., Shazadi, K., Gul, A. & Hakeem, K. R. Effects of pesticides on environment. In *Plant, Soil and Microbes: Volume 1: Implications in Crop Science* (eds. Hakeem, K. R., Akhtar, M. S. & Abdullah, S. N. A.) 253–269 (Springer International Publishing, 2016). https://doi.org/10.1007/978-3-319-27455-3_13.

12. Provost, D. et al. Brain tumours and exposure to pesticides: A case-control study in southwestern France. *Occup. Environ. Med.* 64, 509–514 (2007).

13. James, K. A. & Hall, D. A. Groundwater pesticide levels and the association with Parkinson disease. *Int J Toxicol* 34, 266–273 (2015).

14. Butault, J.-P. et al. Synthèse du rapport d'étude Écophyto R&D: quelles voies pour réduire l'usage des pesticides? *INRA Editions* 93 (2010).

15. Guichard, L. et al. Le plan Ecophyto de réduction d'usage des pesticides en France : décryptage d'un échec et raisons d'espérer. *Cah. Agric.* 26, 14002 (2017).

16. Pimentel, D. *Techniques for Reducing Pesticide Use: Economic and Environmental Benefits*. (Wiley, 1997).

17. Gill, H. K. & Garg, H. Pesticides: Environmental impacts and management strategies. *Pestic. Toxic Asp.* https://doi.org/10.5772/57399 (2014).

18. Epstein, L. & Bassein, S. Patterns of pesticide use in California and the implications for strategies for reduction of pesticides. *Annu. Rev. Phytopathol.* 41, 351–375 (2003).

19. Flint, M. L., Dreistadt, S. H., Clark, J. K. & Program, U. of C. I. P. M. *Natural Enemies Handbook: The Illustrated Guide to Biological Pest Control*. (UC Division of Agriculture and Natural Sciences, 1998).

20. Crowder, D. W., Northfield, T. D., Strand, M. R. & Snyder, W. E. Organic agriculture promotes evenness and natural pest control. *Nature* 466, 109–112 (2010).

21. Luff, M. L. The potential of predators for pest control. *Agric. Ecosyst. Environ.* 10, 159–181 (1983).

22. Drees, B. M. & Leroy, T. R. Evaluation of alternative methods for suppression of crape myrtle aphids. In *Upper Coast 1990–1991 Entomological Result Demonstration Handbook* 21–22 (1991).

23. Kusakari, S., Okada, K., Shibao, M. & Toyoda, H. High voltage electric fields have potential to create new physical pest control systems. *Insects* 11, 447 (2020).

24. Hori, M., Shibuya, K., Sato, M. & Saito, Y. Lethal effects of short-wavelength visible light on insects. *Sci. Rep.* 4, 7383 (2015).

25. Keller, M. D. et al. Laser induced mortality of Anopheles stephensi mosquitoes. *Sci. Rep.* 6, 1–11 (2016).

26. Wilde, W. H. A. Laser effects on some phytophagous arthropods and their hosts1. *Ann. Entomol. Soc. Am.* 60, 204–207 (1967).
27. Dáder, B., Gwynn-Jones, D., Moreno, A., Winters, A. & Fereres, A. Impact of UV-A radiation on the performance of aphids and whiteflies and on the leaf chemistry of their host plants. *J. Photochem. Photobiol. B* **138**, 307–316 (2014).

28. Beard, R. L. Lethal action of UV irradiation on insects. *J. Econ. Entomol.* **65**, 650–654 (1972).

29. Feugier, F. G. Pest control system, pest control method and pest control program. WO 2013/069059 A1. (2013).

30. Sustainable Weed Management in Agriculture With Laser-Based Autonomous Tools | WeLASER Project | H2020 | CORDIS | European Commission. https://cordis.europa.eu/project/id/101000256.

31. Rakhmatulin, I. & Andreasen, C. A concept of a compact and inexpensive device for controlling weeds with laser beams. *Agronomy* **10**, 1616 (2020).

32. Ogawa, K. & Miura, T. Aphid polyphenisms: Trans-generational developmental regulation through viviparity. *Front. Physiol.* **5**, 1–11. https://doi.org/10.3389/fphys.2014.00001 (2014).

33. Peng, X., Qiao, X. & Chen, M. Responses of holocyclic and anholocyclic *Rhopalosiphum padi* populations to low-temperature and short-photoperiod induction. *Ecol. Evol.* **7**, 1030–1042. https://doi.org/10.1002/ece3.2720 (2017).

34. https://www6.inrae.fr/encyclopedie-pucerons_eng.

35. Hullé, M., Chaubet, B., Turpeau, E. & Simon, J. C. Encyclop'Aphid: A website on aphids and their natural enemies. *Entomol. Gen. https://doi.org/10.1127/entomologia/2019/0867* (2020).

36. Luo, C. et al. Hosting certain facultative symbionts modulate the phenoloxidase activity and immune response of the pea aphid *Acyrthosiphon pisum*. *Insect Sci.* https://doi.org/10.1111/1744-7917.12888 (2020).

37. Oh, H.-S. & Kim, J.-S. Clinical Application of CO₂ Laser in "CO₂ Laser" IntechOpen Ed. Dan C. Dumitras 2012. https://doi.org/10.5772/2496.

38. de Miguel, I. et al. Plasmon-based biofilm inhibition on surgical implants. *Nano Lett.* **19**, 2524–2529 (2019).

39. Simonet, P. et al. Direct flow cytometry measurements reveal a fine-tuning of symbiotic cell dynamics according to the host developmental needs in aphid symbiosis. *Sci. Rep.* **6**, 19967 (2016).

40. Legrand, A. & Barbosa, P. Pea aphid (Homoptera: Aphididae) fecundity, rate of increase, and within-plant distribution unaffected by plant morphology. *Environ. Entomol.* **29**, 987–993. https://doi.org/10.1603/0046-225X-29.5.987 (2000).

41. Confédération Générale des Planteurs de Betteraves. Rapport Annuel : Faits & Chiffres 2019. http://www.cgb-france.fr/wp-content/uploads/2019/12/RA2019-Complet-web-bd.pdf (2019).

42. Thermal Power Sensors (C-Series). https://www.thorlabs.com/newgrouppage9.cfm?objectgroup_id=3333&pn=S425C-L.

## Acknowledgements

This work was supported INSA, INRAE, and by the French "Agence Nationale de la Recherche", Project N° ANR-17-CE34-0012-03 (GREENSHIELD). The authors are also grateful to GreenShield Technology for discussion. They also acknowledge in the invaluable help of Rainer Kling (business unit manager), Mathieu Valentin (R&D laser process engineer) from ALPhANOV Company (Bordeaux, France) and Christophe Arnaud (R&D engineer) from ALPhANOV at the time of our study, currently at Solutions NOVIKA (Quebec, Canada). The *A. pisum* YR2 line was a kind donation of Jean-Christophe Simon to F.C. The authors are grateful to Dr. Tom Wood for careful editing of the manuscript.

## Author contributions

R.G., V.L., V.D., P.S, F.C., A.H, B.M. and D.V. set up the experiment hardware. R.G., V.L., A.C., G.D., K.G., P.S. and V.D. carried out the experiment. R.G., V.L., P.S. A.H. F.C. and B.M. analyzed the data. R.G., V.L., P.S, A.H and B.M. wrote the manuscript. V.D., D.V. reviewed the manuscript.

## Competing interests

All authors declare no competing interests, either financial or non-financial, except R.G. who was partially employed by GreenShield Technology company while performing this work.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-021-90782-7.

**Correspondence** and requests for materials should be addressed to B.M.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_11.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021