Bioactivities of rose-scented geranium nanoemulsions against the larvae of Anopheles stephensi and their gut bacteria

Anopheles stephensi with three different biotypes is a major vector of malaria in Asia. It breeds in a wide range of habitats. Therefore, safer and more sustainable methods are needed to control its immature stages rather than chemical pesticides. The larvicidal and antibacterial properties of the Pelargonium roseum essential oil (PREO) formulations were investigated against mysorensis and intermediate forms of An. stephensi in laboratory conditions. A series of nanoemulsions containing different amounts of PREO, equivalent to the calculated LC50 values for each An. stephensi form, and various quantities of surfactants and co-surfactants were developed. The physical and morphological properties of the most lethal formulations were also determined. PREO and its major components, i.e. citronellol (21.34%), L-menthone (6.41%), linalool (4.214%), and geraniol (2.19%), showed potent larvicidal activity against the studied mosquitoes. The LC50/90 values for mysorensis and intermediate forms were computed as 11.44/42.42 ppm and 12.55/47.69 ppm, respectively. The F48/F44 nanoformulations with 94% and 88% lethality for the mysorensis and intermediate forms were designated as optimized formulations. The droplet size, polydispersity index, and zeta-potential for F48/F44 were determined as 172.8/90.95 nm, 0.123/0.183, and -1.08/-2.08 mV, respectively. These results were also confirmed by TEM analysis. Prepared formulations displayed antibacterial activity against larval gut bacteria in the following order of decreasing inhibitory: LC90, optimized nanoemulsions, and LC50. PREO-based formulations were more effective against mysorensis than intermediate. Compared to the crude PREO, the overall larvicidal activity of all nanoformulations boosted by 20% and the optimized formulations by 50%. The sensitivity of insect gut bacteria may be a crucial factor in determining the outcome of the effect of toxins on target insects. The formulations designed in the present study may be a good option as a potent and selective larvicide for An. stephensi.

Introduction alcohols, ethers or oxides, aldehydes, ketones, esters, amines, amides, phenols, heterocycles, and chiefly the terpenes [39]. In recent years, nanoemulsion-based delivery systems have been the focus of many studies. Thanks to their subcellular size, nanocarriers have potential to boost the bioactivity of EOs since they allow a wide tissue penetration and an effortless cellular uptake. Additionally, they render possible to adjust the release of active ingredients at the target site. The nanoemulsion-based delivery systems have been proposed to ameliorate the EOs' physico-chemical properties by decreasing their volatility, improving the stability, and enhancing water solubility, as well as by protecting them from the interaction with the environment [40,41]. These formulations may correspondingly enhance other characteristics, namely phytotoxicity (when applied on vegetation [42]) and spraying improvements (i.e. fog treatments [43]). Nanoemulsion of different plant EOs has hitherto been prepared, and their mosquito larvicidal properties have been investigated [44][45][46][47][48][49][50][51][52][53].
The idea that insects are colonized by numerous microorganisms, particularly bacteria, has widely been acknowledged by life scientists. The external and internal parts of the insect body, especially digestive tract, offer conditions and resources needed to support beneficial microbiota [54]. Several bacterial symbionts promote host fitness by contributing to nutrition [55], reproduction [56,57], speciation [58], immunity/defense [59][60][61][62], ecological communication [63], and pathogen transmission ability [64]. More recent studies have pinpointed that there is a potential linkage between insect gut microbiota and their susceptibility to insecticides [65][66][67]; therefore, the sensitivity of gut bacteria needs to be taken into account prior to any intervention measures.
In spite of its supreme biological potential, the Pelargonium roseum essential oil (PREO) remains, to our knowledge, unexplored concerning the development of an innovative nanolarvicide. Therefore, the present study aimed to assess the larvicidal activity of PREO and its components, to prepare various nanoemulsions of PREO, to select optimized formulations in terms of physical and biological properties and to evaluate the interaction of larval gut bacteria with PREO formulations as a mechanism for the attribution of the oil insecticidal activity.

Biological materials
Late third and fourth instars larvae of two BFs of An. stephensi (mysorensis and intermediate) were purchased from the National Insectarium of Iran, Malaria and Vector Research Group. These strains were maintained separately in the laboratory conditions according to the MR4 standard protocols [68]. The origins of An. stephensi mysorensis and An. stephensi intermediate forms were from Chabahar City in Sistan and Baluchestan Province in Southeastern Iran and from Bandar Abbas City in Hormozgan Province in Southern Iran, respectively. Prior to any experiment, the BFs of the An. stephensi specimens were determined by counting the number of egg ridges, as described in the literature [69].

Identification of the PREO constituents
Chemical constituents of PREO were characterized by gas chromatography-mass spectrometry (GC-MS), which was achieved on a GC (HP 6890, Agilent, USA) equipped with a mass spectrometer detector (HP 5973, Agilent). The MS was operated in the electron ionization mode (70 eV). The MS ion source temperature and the MS quadrupole temperature were kept at 230˚C and 150˚C, respectively. Afterwards, 1 μL of the diluted sample (10 μL in 1 mL of heptane) was injected by autosampler using a 100:1 split ratio and analyzed on a capillary column (TRB-5MS, 30 m, 250 μm, and 0.25 μm). Helium functioned as the carrier gas (99.9995% pure), and its flow rate in the column was adjusted to 1 mL/min -1 . The sample was assessed under the following settings: initial oven temperature at 36˚C for 5 min, ramp-up at 4˚C/min to 200˚C and continued for 8 min, then increased up to 280˚C with a ramp-up of 40˚C/min for 10 min and overall run time of 66 min.
Compounds were determined by comparing their respective mass spectra, retention time, and relative abundance of acceptance match criteria with those of standards and also by comparing with the NIST05 (National Institute of Standards and Technology, Gaithersburg, MD, USA) and the Wiley Registry of Mass Spectral Libraries. The retention were indices calculated for each component using a mixture of n-alkanes (C9-C24) dissolved in n-hexane based on the following formula: Retention Indices = 100 × [n + (Tu-Tn)/(TN-Tn)] Where n = the number of carbons in the alkane preceding compound; N = the number of carbons in the alkane following compound; Tu = the retention time of the unknown compound; Tn = the retention time of the preceding alkane; TN = the retention time of the following alkane.

Larvicidal activity of PREO and its main constituents
Larvae of the two BFs of An. stephensi were exposed to seven (1.56 to 100 ppm) serially diluted concentrations of the PREO for 24 h, according to the methods described previously [35, 70,71]. Due to the hydrophobic behavior of PREO, the oil was initially dissolved in ETH 96% as stock (100 ppm = 100 μg/ml), and subsequent concentrations were prepared by the stock dilution. Also, larvicidal activities of the main components of the PREO (citronellol, L-menthone, linalool, and geraniol) were evaluated both individually and in combination at two concentrations equal to LC 50/90 values in the same way. The components were evaluated individually based on the LC 50 (11.44 and 12.55 μl) and LC 90  For all treatments, 1 mL of each concentration of the PREO solution was supplemented with 99 mL of the dechlorinated water containing 0.001% TW 20, to make up 100 mL of test solution in a 200-mL glass beaker. A plastic rod was applied to stir the oil-ethanol-water solution for 30 s. Four different beakers, one containing 0.125 mg/L of temephos (13.2 μL in 1mL ETH), one comprising of 1% ETH, one holding of 0.001% TW, and the other including untreated dechlorinated water, were set as control solutions. Batches of~20 healthy larvae were gently transferred to the beakers by a fine strainer. The bioassay was carried out in a room with the conditions of 24 ± 1˚C, 50 ± 5% relative humidity, and 12:12 light and dark photoperiodicity. After 24 h of exposure, larval mortality was observed, and four replicates was used for Probit analysis and LC 50/90 calculation. If the mortality in the untreated group was between 5% and 20%, the bioassay was then corrected via Abbott's formula [72].

Preparation of PREO nanoemulsions
To reduce the volatility and improve aqueous solubility of PREO, 16 oil-in-water (O/W) nanoemulsions of the oil were prepared, as described previously [50]. Briefly, different concentrations of TW 20 (2-16%) and ETH 96% (2-32%) were mixed at 600 rpm at room temperature for 5 minutes. Amounts of 11.44 and 12.55 ppm of PREO, equivalent to the calculated percentage of LC 50 values for An. stephensi mysorensis and An. stephensi intermediate, respectively, were added to the TW (as surfactant) and ETH (as co-surfactant) mixture and stirred at the same conditions for 15 min. Subsequently, deionized water was gradually supplemented up to 10 mL and stirred further for 30 min.

Characterization of PREO nanoemulsions
The larvicidal activity of 16 PREO nanoemulsions (F22-F1632) was investigated against BFs of An. stephensi, and the physical properties of the most lethal ones were determined. The oil particle size (PS), polydispersity index (PDI), and zeta-potential (ZP) were measured by dynamic light scattering (DLS) with a Zetasizer Nano ZS (Malvern, ZEN3600, UK) working at 633 nm at 25˚C and equipped with a backscatter detector at 173˚. The droplets' morphology of optimized formulations was also determined by transition electron microscopy (TEM) (RASTAK Lab; Tehran, Iran). The appearance and ZP of the optimum nanoemulsions stored at 4˚C were checked after 40 days.

Comparison of the larvicidal activity of optimized nanoformulations with PREO
Larvicidal activities of optimized nanoemulsions, F48 for An. stephensi mysorensis and F44 for An. stephensi intermediate, were precisely performed as mentioned before. The mortality caused by each formulation was compared with their corresponding PREO LC 50 , 11.5 ppm for An. stephensi mysorensis and 12.5 ppm for An. stephensi intermediate.

Interaction of larval aerobic gut bacteria with PREO and optimized nanoemulsions
Five fourth instar larvae of each An. stephensi form were randomly selected for dissection and identification of gut-associated bacteria. Before dissection, specimens were washed twice thoroughly with sterile PBS (phosphate buffer saline) 1× and then surface sterilized in 70% ETH for 2 min. Guts were removed aseptically within a sterile PBS drop on a microscopic slide under a microbiological laboratory hood. Then guts were mechanically grinded in 1.5-ml Eppendorf tubes containing 100 μl of PBS. The homogenates were inoculated into BHI broth medium and subsequently plated on the BHI agar medium to obtain pure colonies. Colonies with different phenotypes were isolated and sub-cultured successively. Individual bacteria were identified using 16S rDNA sequencing method, as described in the experiments [31]. Based on the frequency of colony-forming units (CFU) on the BHI agar medium, two bacterial isolates, one Gram-positive bacterium and one Gram-negative bacterium, were chosen from each form of An. stephensi to explore the interaction of bacteria with the PREO formulations.
In each interaction test, 32.650 ml of BHI broth, 330 μl of each PREO dilution (equal volume of the LC 50/90 for each biological form as well as their corresponding optimized nanoemulsions), and 20 μl of suspension of each target microorganism (adjusted to 0.5 McFarland standard turbidity corresponding to 10 8 CFU/mL) were mixed together. The interacting mixtures were incubated at 37˚C and shaken at 200 rpm for 24 h. The growth of bacteria was quantified by measuring the OD at 600 nm (Eppendorf BioPhotometer plus, Germany) at different time points, i.e. 2,4,8,12, and 24 hours post interaction (hpi). Experiments were performed in triplicate. Four culture media, comprising of one medium without any additives and two media containing TW and ETH (as negative controls), and a medium including only the desired bacterium (as positive control), were used for the clarification of the interaction results. The percentage inhibition of the growth of bacteria exposed to treatments during 24 h was calculated by the formula: % inhibition = |Test OD-Control OD/Control OD| × 100, where the test group OD corresponded to bacteria treated with antimicrobial compound concentrations, and the control group corresponded to the untreated bacteria grown under normal conditions.

Statistical analysis
The design was a comparative study on the efficacy of the various formulations of PREO in the mortality of the larvae of two BFs of An. stephensi. The LC 50 and LC 90 values for the respective BFs were calculated using probit analysis according to the method described previously [71]. The LC 50 and LC 90 values were chosen based on chi-squared values and degrees of freedom. The significance of the slope probit-log (dose) regression was assessed using the z test (z = b sðbÞ ), if the P values less than 0.05 represented correlations between PREO doses and mortalities. The heterogeneity factor h of the regression equation was calculated to adjust for large χ2. h was defined as h = w2 df If h < 1, the model provided a good fit to the data [73]. SPSS V22.0 software was used for data analysis. Graphs were drawn using GraphPad Prism1 v.5.00 (Graph-Pad Software Inc, San Diego, CA, USA). Mean ± standard deviation (SD) was used to describe quantitative variables. Percentage and frequency were also applied for qualitative variables. The assumptions of normality and homogeneity of variance were tested by Shapiro-Wilks and Brown-Forsythe tests, respectively. The data did not violate these assumptions. The two-tailed Student's t test was used to comparison of larvicidal activity of optimized nanoformulations with corresponding PREO LC 50 . The One-way ANOVA analysis was used to compare the larval mortalities of two BFs exposed to components of PREO, as well as to evaluate growth inhibition of gut bacteria exposed to diverse PREO formulations. If ANOVA test was significant, Tukey's post-hoc analysis was applied. Significance level was considered less than 0.05.

BFs of the An. stephensi
The BFs of An. stephensi strains of Chabahar and Bandar Abbas were confirmed as mysorensis and intermediate, respectively, based on egg ridges counts (detailed results are not shown).

Larvicidal activity of PREO and its main constituents
After 24 h exposure, the PREO showed potent larvicidal activity against BFs of An. stephensi. The mortality in the control groups did not exceed 5% in all concentrations; thus, there was no need for corrections. Temephos caused~90% mortality in both mysorensis and intermediate forms. Mortality of larvae exposed to PREO increased in a dose-dependent manner ( Table 2). Minimum and maximum larval mortality for both BFs were determined to be at 1.56 ppm and 100 ppm concentrations, respectively (Fig 1). The LC 50 and LC 90 values against the larvae of An. stephensi mysorensis were 11.44 and 42.42 ppm, while they were 12.55 and 47.69 ppm for An. stephensi intermediate, respectively (Table 2).
The four major components identified in PREO showed the larvicidal activity against BFs of An. stephensi, as well. The statistical analysis of the larval mortality for each component of the PREO is shown individually in the corresponding columns in Table 3. There were no statistically significant differences in the larval mortality between the PREO components at the dose of 11.44 ppm [F(4,15) = 1.29, P = 0.31]. However, at the concentration of 42.42 ppm, there was a significant difference between the PREO components [F(4,15) = 842.8, P�0.001]. The results of post hoc tests (Tukey's test) revealed a significant difference between citronellol, geraniol, and the mixture of four components (P�0.05) and between the first two mentioned components and L-menthone (P�0.05) with linalool. There was also a significant difference between the studied components at a dose of 12.55 ppm [F(4,15) = 9.60, P = 0.005], as well. The larvicidal activity of L-menthone and geraniol (P�0.05), the mixture of four compounds (P�0.01), and citronellol (P�0.001) were significantly more than linalool. Statistically The highest and the lowest larvicidal activity in both forms of An. stephensi were related to the citronellol and linalool, respectively. Details on the larvicidal activity of PREO and its main constituents are offered in Tables 2 and 3.

Larvicidal activity of PREO nanoemulsions and selection of the optimum formulations
Sixteen O/W formulations of PREO were freshly prepared, and their larvicidal activity was determined as explained before. Among the formulations, F48 and F44 showed the highest

Comparison of larvicidal activity of optimized nanoformulations with PREO
The larvicidal activity of the optimal nanoformulations affecting both BFs of An. stephensi was found to be twice more than that of PREO LC 50 equivalent (P�0.01; Table 4).

Characterization of PREO nanoemulsions
The average PS of F48 and F44 nanoformulations was determined as 172.8 nm and 90.95 nm, and their PDI values were also found to be 0.123 and 0.183, respectively. The ZPs of the two particles were -1.08 mV and -2.08 mV, respectively (Fig 3A, 3B, 3D and 3E). TEM results revealed that the size of the sphere-shaped droplets was about 100 nm (Fig 3C and 3F), accordant with the DLS measures. After 40 days, the opacity of the F48 and F44 nanoemulsions did not change in the storage conditions; however, the ZPs were correspondingly raised to -1.79 and 2.52 (Fig 4). These data indicate that the stability of the prepared nanoemulsions have been maintained during this time. In general, the PREO and optimized nanoemulsions displayed potent antibacterial properties against the gut bacteria of the larvae of An. stephensi (Fig 5). In treatments with only one bacterium, the growth of bacteria was well suppressed for up to 4 hpi, then increased exponentially up to 8 hpi, and eventually became constant to some extents by the end of the interaction. However, in the presence of both Gram-positive and Gram-negative bacteria, the growth pattern occurred relatively in a linear manner (Fig 5).

Interaction of larval aerobic gut bacteria with PREO and optimized nanoemulsions
PREO-based formulations indicated antibacterial activity against four bacteria species in the following order of decreasing inhibitory: LC 90 , optimized nanoemulsions, and LC 50 (Table 5). Overall, the optimized formulations showed better antibacterial activity in the

PLOS ONE
individual interaction mode compared to the LC 50 and in the combined mode compared to the LC 50 and temephos (P�0.01). Also, in the combined interactions, the performance of the optimized formulations was similar to that of LC 90 . Details on the growth inhibition of bacteria exposed to different concentrations of PREO formulations are presented in the Table 5.

Discussion
Herbal EOs, as green pesticides, have practically displayed insecticidal, fumigant, antifeedant, attractive, repellent, and growth-reducing effects on various arthropod pests, including insects of medical and veterinary significance [34, 74,75]. Botanical pesticides made from EOs exploit the toxicity of aromatic hydrocarbons included in the oils [76]. This is an important advantage

PLOS ONE
because these bioinsecticides affect only the target insects, do not damage beneficial natural enemies and provide residue-free food and safe environment [77]. As concluded in the literature, most of the EOs and their active substances are nontoxic to mammals. In general, their LD 50 for rats have been reported to be 800-3,000 mg kg -1 for raw compounds and �5,000 mg kg -1 for formulated insecticides [78]. Thence most EOs and their active constituents can be used without toxicological and ecotoxicological studies [78]. Pelargoniums, with a vast variety

PLOS ONE
of growth habits and habitats, are cultivated both for their beauty as ornamental plants and for their scent, as an important factor, in perfume, food, and beverages industries [79,80]. Previous studies have shown the excito-repellency effects of PREO on adults An. stephensi [81] and on larvicidal activities against Aedes aegypti [82], Culex pipiens [83], and Cx. quinquefasciatus [84]. Nevertheless, its larvicidal impact on An. stephensi, which breeds in varied habitats, has not yet been investigated. The results of the current study disclosed the larvicidal effect of PREO and its major components on two BFs of An. stephensi. These findings are significant from two perspectives; first, the larvicidal activity of pure PREO is acceptable in comparison with the temephos, and second, the susceptibility of the two forms of An. stephensi, mysorensis and intermediate, was found to differ in terms of exposure to the oil. In this study, the PREO, even its major components, was more effective against An. stephensi mysorensis than An. stephensi intermediate (Fig 1 and Table 2). Moreover, the LC 50 achieved an average of 12 ppm for both BFs studied. According to the categories suggested for the larvicidal activity of plant EOs [71], PREO can be categorized into the third class of active compounds, demanding for further consideration and investigation.
In agreement with other studies [83,[85][86][87], the major profiles of PREO were determined as citronellol, geraniol, linalool, and L-menthone by GC-MS analysis (Table 1). Citronellol, well-known as dihydro-geraniol, is the most abundant component of the oil. Geraniol, an acyclic monoterpene alcohol, is actually the main ingredient of citronellol synthesis in many plants. This compound is synthesized via ionization-dependent reaction or converted to citronellol by microbiological reduction [88]. Linalool is also a structural isomer of geraniol. Lmenthone is a monoterpene with a minty flavor that occurs naturally in a number of EOs. With the available libraries, it was possible for us to accurately identify~82% of the components in PREO, while the remaining 18%, which may contain the PREO metabolites, was not distinguished. Although the expected bioactivity was achieved in the identified components, inclusive libraries are required to determine unknown chemical compounds.
The bioactivities of the oxygenated monoterpenes, known in PREO, have been demonstrated in various surveys [88][89][90][91][92]. In our work, all the major components of PREO displayed the larvicidal activity; however, the highest larvicidal activity was related to citronellol and the  [83,85]. Lipophilic citronellol can disrupt the membrane integrity by inducing free radical generation [93]. Geraniol has also been considered as a penetration enhancer in transdermal drug delivery [94]. Likewise, geraniol can restore susceptibility to drugs in numerous Gram-negative bacteria by targeting efflux pumps [95].
Owing to their volatility, EOs have an easier degradation and cause less spread of environmental pollution. However, their solubility and stability problems need to be resolved in some ways, for example through the development of nanoformulations [96]. Accordingly, we developed 16 O/W nanoemulsions. Compared to the crude PREO, the overall larvicidal activity of all nanoformulations was boosted by 20% and the optimized formulations by 50% (Fig 2 and Table 4). F48 with~94% and F44 with~88% larval mortalities were nominated as the optimum formulations for An. stephensi mysorensis and An. stephensi intermediate, respectively (Fig 2). The physical properties of the optimized nanoformulations provided acceptable results under preparation conditions. Although there were no significant changes in the physical characteristics of the formulations kept for up to 40 days, future studies may examine their bioactivities, which were not the purpose of this study.
PREO is stable under normal temperature conditions, and the larvicidal assays of this study were performed under constant laboratory conditions; nonetheless, the effectiveness of EOs could be significantly influenced by post use temperature. This phenomenon has been noticed in the Pavela and Sedlak's study [97] in which the effects of temperature on the insecticidal efficacy of an essential oil from Thymus vulgaris, in terms of acute toxicity against the larvae of Spodoptera littoralis and Cx. quinquefasciatus, have been explored. They reported that the lethal doses significantly decreased with rising/decreasing temperature against S. littoralis/Cx. quinquefasciatus larvae, respectively [97]. Therefore, the effects of ambient temperature on the larvicidal activity of PREO against An. stephensi forms, which are often found in sub-tropical regions, call for further investigation.
The results of this study showed the optimization of the PREO nanoformulations by applying different amounts of surfactant and co-surfactant. Dissimilarity between the optimal formulations affecting the two BFs may reflect the physiological and even genetic variations in the populations of the two studied BFs. Anopheles stephensi has not been regarded as a species complex in the literature; however, morphological differences in the number of egg ridges [69,98], spiracular index [99], genetic variations in the intron I sequences of odorant-binding protein 1 [100], cytogenetic characteristics [101,102], cuticular hydrocarbon profiles [103], and disparities in ecological, behavioral and mating characteristics [104], together with the findings of this study call into question this hypothesis.
The contributions of insect gut microbiota in the susceptibility and resistance to antimicrobials have recently been the focus of researchers [67,105]. The function of gut bacteria against toxins can be protective or synergistic. Basically, insect gut bacteria are resistant to antimicrobials and to toxic ingredients of the food, or at least actively participate in detoxifying such compounds [106]. In this context, the antibacterial effects of PREO and its optimized nanoemulsions were investigated on representative gut bacteria of two BFs of An. stephensi. The whole PREO-based formulations displayed potent antibacterial activities against the symbiotic bacteria of the mysorensis and intermediate forms. However, the bacterial growth inhibition results varied depending on the species and Gram staining (Fig 5). Surprisingly, in the treatments with both Gram-positive and Gram-negative bacteria, all the PREO-based formulations uncovered a better inhibitory effect than temephos. This behavior may reflect the interaction of the bacteria under study, which encourages future studies to consider the sensitivity or resistance of the whole microbial community in the evaluation of antimicrobials. It should also be clarified in future investigations whether the antibacterial consequence of PREO inhibits the bacterial growth (bacteriostatic) or destroys bacterial cells (bactericidal), as well as determines the mechanisms underlying such actions.
EOs are active against a broad variety of organisms with multiple mechanisms, though their anti-insect mode of action can be categorized as behavioral and physiological. Behavioral mode of action will be true in adult insects in which volatile components of the oil can disrupt the communication behavior of insects by blocking the function of antennal sensilla [107]. The activities of some EOs appear to be the result of effects on the insect nervous system, either by the inhibition of acetylcholinesterase [108] or by the antagonism of the octopamine receptors [109]. The rapid action against some pests is indicative of a neurotoxic mode of action. The lack of octopamine receptors in vertebrates provides the mammalian selectivity of EOs as insecticides [110,111]. Therefore, plant EOs and their components can influence the physiological functions of different insect species and modulate their gut microbiota, as the results of the present study indicated.
In the present research, the antibacterial performance of PREO nanoemulsions was comparable with temephos, and this type of performance may be attributable to their antibacterial activity; however, further (proteomic and enzymatic) studies are needed to identify the inhibitors of enzymes involved in susceptibility/resistance to toxic compounds.
As a non-systemic organophosphorus insecticide, temephos has extensively been used to treat water infested with disease-carrying insects, including mosquitoes, midges, and black fly larvae. It has also been recommended for the treatment of drinking-water sources and containers [112,113]. Regrettably, resistance (and tolerance in this study) to temephos has been detected in the Culex spp. [114,115], Aedes spp. [116][117][118], and An. stephensi [119,120]. This issue sounds the alarm for using temephos in Iran and around the world.

Conclusion
The results of this study approved the larvicidal and antibacterial effects of pure PREO and its nanoformulations on two BFs of An. stephensi and their intestinal bacteria. Difference between the PREO-based formulations affecting the two BFs may reflect the physiological and even genetic variations in the populations of An. stephensi studied. The larvicidal activity of PREObased formulations can ascribed to their antibacterial activity, which calls for further investigation. The susceptibility or resistance of the intestinal bacterial flora may determine the outcomes of the pesticides effects on the target insects. This situation affirms the potential of the bacteria as a major contributor to the survival or inexistence of insects. The formulations developed herein are cost-effective and possibly will have the least damage to humans, environment, and non-target organisms and would be comparable to industrial larvicides, e.g. temephos, if they pass semi-field and field trials successfully.