First national-scale evaluation of temephos resistance in Aedes aegypti in Peru

Background The development of resistance against insecticides in Aedes aegypti can lead to operational failures in control programs. Knowledge of the spatial and temporal trends of this resistance is needed to drive effective monitoring campaigns, which in turn provide data on which vector control decision-making should be based. Methods Third-stage larvae (L3) from the F1 and F2 generations of 39 Peruvian field populations of Ae. aegypti mosquitoes from established laboratory colonies were evaluated for resistance against the organophosphate insecticide temephos. The 39 populations were originally established from eggs collected in the field with ovitraps in eight departments of Peru during 2018 and 2019. Dose–response bioassays, at 11 concentrations of the insecticide, were performed following WHO recommendations. Results Of the 39 field populations of Ae. aegypti tested for resistance to temephos , 11 showed high levels of resistance (resistance ratio [RR] > 10), 16 showed moderate levels of resistance (defined as RR values between 5 and 10) and only 12 were susceptible (RR < 5). The results segregated the study populations into two geographic groups. Most of the populations in the first geographic group, the coastal region, were resistant to temephos, with three populations (AG, CR and LO) showing RR values > 20 (AG 21.5, CR 23.1, LO 39.4). The populations in the second geographic group, the Amazon jungle and the high jungle, showed moderate levels of resistance, with values ranging between 5.1 (JN) and 7.1 (PU). The exception in this geographic group was the population from PM, which showed a RR value of 28.8 to this insecticide. Conclusions The results of this study demonstrate that Ae. aegypti populations in Peru present different resistance intensities to temephos, 3 years after temephos use was discontinued. Resistance to this larvicide should continue to be monitored because it is possible that resistance to temephos could decrease in the absence of routine selection pressures. Graphical Abstract Supplementary Information The online version contains supplementary material available at 10.1186/s13071-022-05310-x.

which has an arid temperate climate with an annual rainfall of only 2 mm [6,7]. These adaptations have led to a wide dispersal of the vector, which in turn has led to the presence of arboviruses that often follow the same pattern of dispersal.
It is estimated that about 2.5 billion people, representing 40% of the human population worldwide, live in areas at risk of dengue transmission [8] and that 390 million cases of dengue occur per year in tropical and subtropical areas [3]. In addition, in the last 5 years, Ae. aegypti has been responsible for the spread of chikungunya and Zika to regions of the Americas, placing a significant burden on healthcare systems and also causing social and economic disruption [9]. While Ae. albopictus is widely distributed in the Americas (present in 21 of 44 countries), it has not been detected in Peru. However, Peru shares borders with countries reporting widespread distribution of Ae. albopictus, such as Brazil, Colombia, Bolivia and Ecuador (where Ae. albopictus was first detected in 2017) [10].
In Peru, the vector of dengue, chikungunya and Zika is Ae. aegypti. Aedes aegypti was eradicated locally in 1958, but it subsequently recolonized the country in 1984 [11]. The first outbreak of dengue occurred 6 years later, attributed to dengue virus serotype 1 (DENV-1) [12]. Low-incidence outbreaks occurred thereafter up to 2001, when a major outbreak occurred, with 23,304 cases. This outbreak confirmed the circulation of the four serotypes of DENV [13]. Significant outbreaks occurred between 2017 and 2020, with 68,290 cases and 52,826 cases, respectively [14]. The first reported cases of chikungunya and Zika occurred in 2015 and 2016, respectively [15,16].
In Peru, Aedes-borne arboviruses are present in three ecological regions: the coast, the Amazon jungle and the Andes Mountains [17,18]. The latter region (Andes Mountains) presents such a diversity of altitudes that is differentiated into two areas: the Andean highlands (> 2300 m a.s.l.) and high jungle (400-1400 m a.s.l.; located on the eastern flank of the Andes) [19]. Aedes-borne arboviruses only affect areas of lower altitude. Aedes aegypti are widely distributed in 21 of Peru's 24 departments and the constitutional province of Callao, and have been identified in 527 districts [20], where approximately 22 million people live, putting 70.4% of the Peruvian population at risk of contracting arboviruses transmitted by this vector.
The main dengue control strategy implemented in Peru until 2016 was focal treatment of larval habitats with the organophosphate (OP) insecticide temephos due to its easy dosage, application and acceptability by the community. For the same reasons, the adulticides used in the 1990s were the OPs fenitrothion and malathion [21]. Since the beginning of this century, pyrethroid (PY) insecticides have also been used (cyfluthrin, deltamethrin, alpha-cypermethrin and cypermethrin) [17,22], applied using thermal and cold fog equipment that can be manually carried or mounted on trucks. However, in 2015, Pinto et al. [23] detected resistance against PYs in field-caught Ae. aegypti populations in association with knockdown resistance (kdr) mutations (Phe1534Cys, Val1016Ile). This led the Peruvian Ministry of Health (MINSA) to implement a control strategy in which PYs were rotated with the OP malathion. For larval control, the OP temephos was switched with the insect growth regulator (IGR) pyriproxyfen because temephos shares the same mode of action as the OP malathion. This switch away from using temephos, however, was made without any knowledge of the susceptibility status of Ae. aegypti populations to this insecticide.
Surveillance of insecticide resistance in arthropod vectors in Peru is an activity that falls under the Instituto Nacional de Salud (INS; National Institute of Health) and Regional Reference Laboratories (LRRs) which, in the context of a decentralized healthcare system, requires coordinated work. However, many of the LRRs do not have entomology laboratories or insectaries and thus do not have the facilities or knowledge to routinely perform the various tasks/responsibilities placed on them [22,24]. This situation results in a weakening of the vector surveillance and control programs, despite ongoing increases in the numbers of cases of Aedes-borne arboviruses. The continuous use of insecticides results in selective pressures that in turn drive physiological and/or behavioral adaptation, a phenomenon known as resistance [25]. Insecticide resistance can lead to operational failures, ultimately requiring rotations of insecticides with different modes of action or mosaic treatments to manage insecticide resistance [26]. Ideally, these strategies should be carried out preemptively to preserve insecticide efficacy, but also reactively to mitigate or reverse resistance [25].
The aim of the present study was to determine the levels of resistance to temephos in 39 Peruvian field populations of Ae. aegypti from the three ecological regions of Peru after having used temephos for > 25 years, following an interruption in its use of 3 years.

Sampling and study area
The mosquito population tested comprised the F1 or F2 generations of Ae. aegypti from colonies maintained at INS that were established from eggs collected in 39 localities in the Peruvian departments of Tumbes, Piura, La Libertad, San Martín, Loreto, Ucayali, Junin and Madre de Dios during April 2018 and January 2019 (Fig. 1). The colonies were obtained from eggs collected with ovitraps, which were distributed every 200 linear meters covering the urban area of the locality, according to the parameters established by MINSA [27]. The ovitraps were installed in intra-and peri-domestic areas and contained a strip of paper towel as oviposition substrate [27,28] and 10% hay infusion as attractant [29]. The localities sampled are located in three ecological regions showing distinct weather patterns [6], variations in the incidence of dengue [17,18] and Ae. aegypti populations with differing insecticide resistance profiles [23].

Mosquito collection and laboratory rearing
The paper strips containing eggs were left to dry for 1 week and then soaked in water to hatch the F0 generation [30]. The F1 and F2 generations were subsequently established in the insectary of the Laboratorio de Referencia National de Entomología (LRNE), INS, Lima, Peru (Table 1). Humidity and temperature conditions in the insectary were maintained at 70 ± 10% relative humidify and 26 ± 2 °C, respectively. Aedes aegypti of the susceptible Rockefeller reference strain were used as a control for the resistance tests [31].

Larvicide susceptibility testing
Dose-response susceptibility tests were performed with temephos following WHO recommendations [32]. The larvae were exposed to a wide range of concentrations of the insecticide with the aim to evaluate larvicidal activity and thereby determine the 50% and 95% lethal concentration values (LC 50 and LC 95 ) for each study population. Four replicates per concentration, and 11 concentrations were evaluated with 20 third-stage larvae (L3) in 100 ml of insecticide solution per replicate. Insecticide solutions were prepared with ethanol solvent and temephos active ingredient (Chem Service, West Chester, PA, USA), using a concentration range of 0.004 to 0.324 mg/ ml for the field populations, and 0.002 to 0.012 mg/ml for the control Rockefeller strain. Simultaneously, a control group with four replicates exposed only to 600 µl of ethanol solvent was evaluated [32]. Each test was carried out on three different days to ensure the reproducibility of the method and consistency of the results [32,33].

Data analysis
The values for LC 50 and LC 95 were calculated from a log dosage-probit mortality regression line using probit analysis (Polo-PC statistics package [34]) for each population. Resistance ratios (RRs; i.e. LC of field population/LC of susceptible population) were also calculated to define the intensity of resistance in the field populations. Specifically, Eq. 1 was used to calculate the RR 50 (where i = 50) and the RR 95 (where i = 95 as shown:: The population was considered to be susceptible when the RR < 5; when the RR was between 5 and 10, the population was considered to have moderate resistance; and when the RR ≥ 10, the population was to be considered highly resistant.

Results
Third-stage larvae from the F1 and F2 generations of 39 Peruvian populations of Ae. aegypti mosquitoes from localities with different incidences of arboviruses and located in three ecological regions of Peru were tested. A total of 37,440 L3 were tested across all bioassays. Table 2 and Fig. 2 show the values of the LC 50 and LC 95 , the resistance ratios (RR 50 and RR 95 ) and the slope of the probit regression lines for the insecticide temephos in the 39 populations of Ae. aegypti. The differences between the RR 95 of all populations were analyzed according to the criteria of Mazzarri and Georghiou [35] and Sá et al. [33].The results showed that the three regions could be grouped into two groups: coastal and jungle (the latter including the Amazon and high jungles). Most of the population in the coastal group showed high levels of resistance to temephos (RR > 10), with three populations showing values > 20 (AG 21.5, CR 23.1, LO 39.4; for abbreviations of localities/populations, see Table 1 and Fig. 1). Some populations in th coastal group (PG, MI, TA, BE, CS, ES and LA) showed moderate levels of resistance (5 ≤ RR < 10) and only four populations (JO, CH, CA and VI) were considered to be susceptible (RR < 5). In the jungle group, most of the populations showed moderate levels of resistance (5 ≤ RR < 10), with RR values ranging from 5.1 in JN to 7.1 in PU, while some populations were considered to be susceptible (RR < 5), with values ranging from 2.2 in MO to 4.9 in YA and PA; the notable exception in the jungle group was the PM population, which was considered to be highly resistant (RR = 28.8).
However, the RRs of the populations were heterogeneous, even within the same department. In Tumbes, for example, most populations were highly resistant to temephos (RR > 10), except for the PG population (RR = 8.9) which showed moderate resistance to this insecticide. Piura, in comparison, had two populations (MA and LO) that exhibited high RR values (RR > 19), but also had two susceptible populations with low values (RR < 5). In La Libertad, the RRs varied between 2.1 in CA and 11.3 in FL, showing that the mosquito populations in this department also showed heterogeneous resistance RRs; however, there were also two populations in which the lowest RRs were detected (RR < 2.5). The departments of San Martin, Junín, Loreto and Ucayali were more homogeneous, with RR values ranging from 2. were considered to be more susceptible (RR < 5) in both the RR 50 and RR 95 . In general, jungle localities were more homogeneous.
The slope values of the probit regression lines of the field populations were lower than those obtained with the Rockefeller lineage, except in the CA population (slope = 7.18) (Table 2; Fig. 3). This result confirms the heterogeneity of temephos resistance in the field strains in relation to the reference strain.

Discussion
As part of a national strategic plan to eliminate dengue, chikungunya and Zika, a national survey to ascertain the susceptibility of Ae. aegypti to temephos was undertaken. The Peruvian INS, in coordination with Regional Health Directorates, surveyed 39 dengue endemic localities from April 2018 to January 2019. Overall, Peruvian populations of Ae. aegypti showed variable levels of resistance to temephos. Based on the results, these populations were clustered into two groups: (i) coastal populations with moderate to high levels of resistance to temephos, and (ii) jungle populations, which showed a moderate level of resistance to temephos, and some susceptibility, except for the Puerto Maldonado population which showed a very high level of resistance. Interestingly, from 2015 to 2017, the dengue cases reported in the coastal departments accounted for between 54 and 86% of the total number of dengue cases reported in Peru, while from 2018 to 2021, dengue cases reported in the jungle regions accounted for between 55 and 87% of the total number of dengue cases reported [36,37] (Additional file 1: Table S1). Higher levels of resistance to temephos were detected in coastal populations than in jungle populations, despite the northern coast having been recolonized by Ae. aegypti 10 years later than the Peruvian Amazon [38,39]. The reasons underlying this difference in susceptibility are likely diverse. One possible explanation is based on the notion that Ae. aegypti entered the country along its border with Ecuador [11]. It is known that Ae. aegypti populations in Ecuador have variable levels of resistance to temephos, with populations from Huaquillas (bordering Peru) and Arenas showing susceptibility, while the populations from San Lorenzo and Nueva Loja show high levels of resistance. This variability could be due to variations in the intensity and frequency of temephos use across Ecuador [40]. Another possible explanation is climate differences. The departments along the northern coast of Peru were the most affected by the El Niño climatic phenomenon (1997)(1998), with heavy rains that caused an increase in malaria and dengue cases [41,42]. This situation led to the establishment of prevention and control activities, including insecticide space spraying and indoor residual spraying, mapping and treatment of larval habitats (including chemical control with temephos) and campaigns to eliminate potential larval habitats. These activities were carried out periodically in 333 localities in the departments of Tumbes, Piura, Lambayeque and La Libertad, with the goal of eliminating larval habitats in urban, peri-urban and rural areas.
In addition, the urban population on the coast is growing, but there is a lack of planning and organization and thus insufficient basic sanitation facilities due to, among other reasons, constant migration from other regions due to violence, lack of job opportunities, education and access to technological resources, as well as poverty [42]. This rapid, unplanned urbanization favors the introduction and establishment of Ae. aegypti. However, the west coast of Peru is a long desert strip interrupted by valleys, and the main characteristic of the region is a scarcity of rainfall [6], which should serve as a limiting factor for the transmission of arboviruses in this region, with the one important exception being the north coast, which experiences high temperatures and rainfall in the summer due to its proximity to the equator [18]. Overall, therefore, the potential benefit of the mainly arid climate is partially or totally negated by the presence of larval habitats inside homes as a consequence of poor water storage facilities due to an inadequate piped water supply in urban centers [17,26]. Evidence of this can be seen in the types of larval habitats that predominate in and around homes, including water storage drums, cylinders, wells and flower vases [7]. As a result, the lack of access to and availability of water, as well as the poor quality of water, create conditions for the proliferation of vectors and transmission of arboviruses [43], as recognized by the Ministry of Health of Peru, which reported in 2016 that inadequate access to water was associated with 41.2% of cases of dengue [44]. The ecosytems of the Peruvian Amazon differ from those on the coast. These jungle areas typically include humid and rainy tropical forests [6]. In this region, due to the constant rainfall, there is a greater abundance of potential larval habitats, especially those known as "los inservibles" (discarded objects, passively filled with rainwater). Unlike along the coast, these larval habitats are typically peri-domestic and are not used intentionally to hold water [45]. In contrast to the urbanization along the coast, the jungle is the least populated region of the country, accounting for only 10% of the national population, and the infrastructure of Amazonian villages often reflects inadequate access to basic services, education and health [46]. The city of Puerto Maldonado is the capital of the Amazonian department of Madre de Dios, located in the southeastern part of the country. The Ae. aegypti population in this city was the only one in the Peruvian Amazon that presented a high level of resistance to temephos. The vector was introduced into this city in 1999 [47], 15 years after it first being recorded in the Peruvian jungle in 1984 [11]. Dengue cases were sporadic in the period 2000-2016, with only one major outbreak (2952 recorded cases) reported in 2010 [48]. This low burden of dengue is supported by the findings of Salmón-Mulanovich et al. [49], who found low seroprevalence to DENV in a retrospective study conducted in Puerto Maldonado in 2018. This low disease burden suggests that vector control activities were not routine or intense and that the local Ae. aegypti population has not experienced strong selection pressure from the larvicide temephos. The explanation for the high levels of temephos resistance detected in the present study could be due to the vector having spread across the border from Brazil and Bolivia. On the Brazilian side, the Ae. aegypti population from Rio Branco (Acre State) is highly resistant to temephos [50] and similarly, on the Bolivian side, the Ae. aegypti population in the border city of Cobija (department of Pando) shows moderate levels of resistance to temephos [51].
It is also important to consider that the use of pesticides in the agricultural sector in Peru is intense at the national level, with the most widely used insecticides being OPs [52]. In a historical context, organochlorine (OC) pesticides were used between 1940 and 1950 (dichlorodiphenyltrichloroethane [DDT], benzene hexachloride [BHC] and toxaphene), followed by the introduction of the OPs parathion and methamidophos in 1950 and their use for several decades, with parathion used until 2005 [52] and methamidophos until 2018-2019 [53]. In the 2000s, pyrethroids (deltamethrin, cypermethrin and alphacypermethrin) and carbamates [CA] (carbofuran, methomyl, carbosulfan and carbaryl) were introduced [52,53]. A study conducted in 2012 reported that 43% of farmers preferred to use OPs because they had a broad spectrum of action (contact, ingestion and fumigant effects) [54]. It is important to note that these pesticides were widely used on cotton, corn and potato crops in the northern and central coastal valleys, as well as in the Andes Mountains. Another important consideration is the illegal trade in pesticides through smuggling (mainly along the northern border), street sales and adulteration and counterfeiting of products, especially in the northern coast, central and southern parts of the Andes Mountains [55]. The continuous use of pesticides in agriculture has led to pest resistance to OCs and OPs due to misuse of pesticides and lack of pesticide management [52]. The agricultural sector is using the same classes of insecticides that are being used for public health, so mosquitoes and other non-target insects may experience selection pressure from insecticides used in agriculture, resulting in the selection of populations that exhibit resistance to multiple insecticides [56].
The results of the present study demonstrate that Peruvian Ae. aegypti populations show diverse levels of Palomino et al. Parasites & Vectors (2022) 15:254 resistance to the OP temephos, and are consistent with resistance patterns observed in other field populations that have been subjected to intense selection pressure from temephos [57]. In Peru, temephos was used for > 25 years; it is thus not surprising that resistance to this insecticide has reached very high levels in some areas. However, it is not possible to estimate or quantify the evolution of this resistance due to the lack of baseline information; the absence of data from Peru is also found in the review of resistance prepared by Moyes et al. [58]. Resistance to temephos has been reported worldwide, with high levels of resistance reported in Tamil Nadu (India) [59], Caldas (Colombia) [60], Pernambuco (Brazil) [61], Martinique [62] and Bahia (Brazil) [63]; moderate levels of resistance reported in Tocantins (Brazil) [33], Laos (Asia) [64], Paraná (Brazil) [65], Quindío (Colombia) [66], Delhi (India) [67], Sao Paulo and Northeast Region (Brazil) [68]; and susceptibility reported in Malaysia [69], Santiago Island (Cape Verde) [70] and Phitsanulok province in Thailand [56].
Chemical control of Ae. aegypti with temephos in Peru was continuous until 2015 when it was replaced by pyriproxyfen [71]; this replacement could explain the moderate and low levels of resistance found in the populations of Pampa Grande, Bellavista and Comunidad Saludable, among others. This finding suggests that resistance to temephos is unstable in the absence of continuous selection pressure. The notion of unstable resistance is supported by the findings of various researchers. Wirth and Georghiou [57] who detected significant decreases in the levels of resistance to temephos in the British Virgin Islands (Tortola) following an interruption of > 10 years in its use, from a RR of 46.8 in 1985 [24] to a RR of 12.1 in 1992-93 [72] and then to a RR of 6.3 in 1995-96 [73]. Similarly, in Colombia (Caldas), Conde et al. [60] observed a reduction in the level of resistance to temephos at > 4 years after discontinuation of its use, from RRs of 13.27 and 11.48 in 2007 to RRs of 4.75 and 5.61 in 2011. Similarly, in Brazil (Juazeiro do Norte) [74], a decrease in resistance to temephos was observed, from a RR of 10.4 to a RR of 7.2, 7 years after this larvicide was replaced by Bacillus thuringiensis israelensis (Bti). Also in Brazil, Rahman et al. [75] observed a reduction in the level of resistance to temephos in Rio de Janeiro, in the municipality of Campos dos Goytacazes, 15 years after discontinuationof its use, from a RR of 7.8 in 2001 to a RR of 2.6 in 2016. These authors also detected a significant decrease in the levels of resistance to temephos in the municipality of Itaperona, from a RR of 25.6 in 2011 to a RR of 7.3 in 2016, after this larvicide was substituted by IGR [75]. Consequently, rotating to a new insecticide with a different mode of action could be advantageous in terms of temephos resistance management. The WHO recommends the following compounds as alternative larvicides: Bti, diflubenzuron, methoprene, novaluron, pyriproxyfen and spinosad.
Following the observed resistance of Peruvian Ae. aegypti to pyrethroids, the OP malathion is being reintroduced [71]. It is uncertain how long this insecticide will remain effective if resistance to temephos has already been demonstrated or if there is cross-resistance between temephos and malathion. Wirth and Georghiou [57] suggested that resistance to malathion did not increase significantly under selection pressure with temephos and that adulticides exerted lower selection pressure than larvicides [28]. If this is indeed the case, it is reasonable to believe that malathion may still be effective in Peru, which is important considering the few chemical options available in Peru for vector control.
An important limitation to this study was the lack of additional evaluations that would allow us to better understand the evolution of resistance to temephos, as well as the characterization of the observed resistance by molecular and enzymatic methods. On this last point, Rodriguez et al. [76], in a study with Latin American populations, found that a Peruvian Ae. aegypti population presented variations in the intensity of resistance to different OPs (temephos, RR = 30; malathion, 1.5; fenthion, 6.6; pirimiphos-methyl, 10; fenitrothion, 1.1; chlorpyrifos, 4.3), with elevated activity of esterases related to resistance to temephos, while mono-oxygenases were associated with resistance to pirimiphosmethyl and chlorpyrifos.