Drosophila melanogaster Response to Feeding with Neomycin-Based Medium Expressed in Fluctuating Asymmetry

The fruit fly Drosophila melanogaster is a model species used for a wide range of studies. Contamination of Drosophila cultures with bacterial infection is common and is readily eradicated by antibiotics. Neomycin antibiotics can cause stress to D. melanogaster’s larvae and imagoes, which may affect the interpretation of the results of research using culture from neomycin-based medium. In the present study, fluctuating asymmetry (FA), one of the important bioindicators of stress, was measured. Larvae and imagoes of a wild-type D. melanogaster strain were exposed to various concentrations of neomycin. The size of anal papillae and selected wing veins were measured using scanning electron and light microscopy, respectively. Next, the FA was checked. The values obtained for larval anal papillae appeared to be concentration-dependant; the FA indices increased with the concentration of neomycin. The wing FA presented a large but variable correlation, depending on the measured vein. However, the mean length of veins was the highest for the control group, with neomycin-exposed groups showing lower values. The research showed that neomycin may cause sublethal stress in D. melanogaster, which manifests in increased FA indices. This suggests that neomycin can cause physiological and developmental stress in insects, which should be taken into account when interpreting the results of studies using these model organisms.


Introduction
The fruit fly Drosophila melanogaster (Meigen) (Diptera: Drosophilidae) is a commonly used model for a wide range of studies, because of its short life cycle, adaptation to different ecological conditions, and ease of rearing in the laboratory. Previous experiments have demonstrated that different concentrations of some organophosphorus insecticides, plant secondary metabolites, and food dyes have toxic effects on survivorship of D. melanogaster [1][2][3][4]. Antibiotics that possess antibacterial, antifungal, have shown differences in the morphology of wings reflected in the FA are sensitive indicators of environmental stress [36][37][38]. Moreover, Debat et al. [36] indicated that wing shape can vary in more flexible ways and be more informative than simple traits such as the size of other body parts when investigating deteriorated larval survival, development, and adult fitness of insects exposed to environmental stressors (temperature, insecticides, and other chemical pollutants, dietary additives such as antimicrobial agents). Since the method of FA analysis is easy to apply, fast, and relatively cheap, it can be routinely used to compare numerous potential stressors on various species, worldwide.
This method was therefore applied to determine the possible stressful effects of sublethal concentrations of neomycin on D. melanogaster after rearing larvae on diets amended with the aminoglycoside antibiotic neomycin.

Drosophila melanogaster Culture
Wild-type (Oregon R strain) D. melanogaster from a laboratory culture was used in this study. The experiments were carried out under laboratory conditions. The stock culture was xenically maintained by rearing first instars to adults on an artificial diet [39,40]. The colony was maintained and the tests were conducted at 25 ± 2 • C, 60-70% relative humidity, and on a photoperiod of 12:12 (L:D) h. Newly emerged adults were used to maintain the stock culture. Insects were reared in 15 mL vials (26.5 × 58.7 mm, Yıldız Chemistry Co., Ankara, Turkey) with~5 mL of artificial diet. The standard diet was composed of 8 g of agar-agar ultrapure (Merck, Darmstadt, Germany), 20 g of sucrose (BioUltra, 99 %, Sigma Chemical Co., St. Louis, MO, USA), 11.78 g of dried powder yeast (Dr. Oetker Food Co., Torbalı-İzmir, Turkey), 0.8 g of ascorbic acid (BioUltra, 99%, Sigma), 7.72 mL of nipagine (SigmaUltra, p-hydroxybenzoic acid methyl ester, crystal), 36 g of mashed potatoes (Knorr, Unilever Co., Ümraniye, Istanbul, Turkey), and 1000 mL of distilled water. Diet in liquid form was poured into the individual vials, and the vials were maintained at room temperature for 30 min to allow the diet to solidify. Ten or 15 males and females of the same age were placed in individual vials and maintained for 24 h, the time necessary for fertilization and oviposition. Vials were plugged with cotton to prevent escape of the flies and to ensure adequate ventilation. After 24 h, adults were transferred into new vials with the control and amended diets according to treatment group as per Reis [41]. Transfers occurred every 3 days. Methods used to prepare and dispense diets and the rearing of instars are described in detail in Roberts [42].

Feeding Experiments
Neomycin sulfate (white powder form, neomycin O-2,6-diamino-2,6-dideoksi-α-d-glukopiranozil- , sulfate (trisulfate salt, aqueous) water-soluble in white powder form 90-95%, Sigma Chemical Co., St. Louis, MO, USA) was directly incorporated into diets at four different concentrations of 150, 300, 600, and 900 mg/L diet while the diets were in the liquid state. These concentrations were based on the results of previous experiments with sodium tetraborate and boron derivatives studied using D. melanogaster [43,44] and other insects [7,45], because it allowed fruit flies to complete their development to the adult stage in pre-feeding tests. Control diets had no neomycin sulfate. Neonates were placed into separate vials using a fine brush and reared until adult emergence on the artificial diets amended with given concentrations of neomycin, which took about 10 days. Newly emerged third instar larvae and adults from each assay were selected for fluctuating asymmetry (FA) observations. Due to the size of the larval anal papillae and their location, scanning electron microscopy (SEM) was used to measure them. Larvae were fixed in 2% glutaraldehyde, buffered with 0.1 M sodium cacodylate (pH 7.2) for 2 h, and then postfixed in 1% osmium tetroxide and dehydrated through a series of ethanol-water solutions of increasing concentration. Samples were critical-point dried, coated with gold, and observed using a Zeiss Evo 40 SEM (Figure 1). Due to the size of the larval anal papillae and their location, scanning electron microscopy (SEM) was used to measure them. Larvae were fixed in 2% glutaraldehyde, buffered with 0.1 M sodium cacodylate (pH 7.2) for 2 h, and then postfixed in 1% osmium tetroxide and dehydrated through a series of ethanol-water solutions of increasing concentration. Samples were critical-point dried, coated with gold, and observed using a Zeiss Evo 40 SEM (Figure 1). Figure 1. SEM picture of measured anal papillae. Drosophila melanogaster larvae exposed to 300 mg neomycin/L. For larvae, the length of the two anal papillae was measured twice (left and right appendix were distinguished) by two independent people, and the absolute difference in the length of two appendixes generated similar results (Table 1), indicating high repeatability. The FA index was defined as an absolute difference between the right and left appendixes, standardized by the mean. The presence of asymmetry was confirmed using an ANOVA test (Table 2), which tested the differences in length between the longer and shorter appendix in a pair, as previously described [46,47].  For larvae, the length of the two anal papillae was measured twice (left and right appendix were distinguished) by two independent people, and the absolute difference in the length of two appendixes generated similar results (Table 1), indicating high repeatability. The FA index was defined as an absolute difference between the right and left appendixes, standardized by the mean. The presence of asymmetry was confirmed using an ANOVA test (Table 2), which tested the differences in length between the longer and shorter appendix in a pair, as previously described [46,47]. Table 1. Correlation between two independent measurements (calculated from absolute values of differences between right and left anal papillae of D. melanogaster larvae), *** p < 0.001.

Imagoes
Wild-type Drosophila melanogaster imagoes from each assay were selected. FA index was defined as the absolute difference between the right wing and the left wing veins standardized by the mean. The length of six veins was measured under the microscope (4× magnification), using the CellB program (Olympus Co.) ( Figure 2) as the base to compute the FA indices.

Imagoes
Wild-type Drosophila melanogaster imagoes from each assay were selected. FA index was defined as the absolute difference between the right wing and the left wing veins standardized by the mean. The length of six veins was measured under the microscope (4× magnification), using the CellB program (Olympus Co.) ( Figure 2) as the base to compute the FA indices. The length of veins (from the left and the right wings) was measured twice by two independent people, and the absolute difference in the length of two appendixes generated similar results (Table  3), indicating high repeatability. A-F-the measured veins, for their location see Figure 2; N-sample size.
The presence of asymmetry was confirmed in three veins, B, C, and D, by the ANOVA test (Table 4), which tested the differences in length between longer and shorter veins in the pairs [46,47]. Therefore Veins A, E, and F were not included in further analyses.

Biometric and Statistical Analyses
The type of asymmetry of the two anal papillae was revealed by testing the normality of the right-minus-left value distribution using the Shapiro-Wilk test, which detected the absence of anti-symmetry (1: p = 0.56, 2: p = 0.71). Anti-symmetry and trait size were checked according to the The length of veins (from the left and the right wings) was measured twice by two independent people, and the absolute difference in the length of two appendixes generated similar results (Table 3), indicating high repeatability. A-F-the measured veins, for their location see Figure 2; N-sample size.
The presence of asymmetry was confirmed in three veins, B, C, and D, by the ANOVA test (Table 4), which tested the differences in length between longer and shorter veins in the pairs [46,47]. Therefore Veins A, E, and F were not included in further analyses.

Biometric and Statistical Analyses
The type of asymmetry of the two anal papillae was revealed by testing the normality of the right-minus-left value distribution using the Shapiro-Wilk test, which detected the absence of anti-symmetry (1: p = 0.56, 2: p = 0.71). Anti-symmetry and trait size were checked according to the Insects 2020, 11, 378 6 of 17 protocol described by Kozlov and Niemelä [46]. Absence of directional asymmetry was confirmed by the zero mean value (1: p = 0.46, 2: p = 0.58). Thus, observed asymmetry was classified as FA.
For imagoes, the type of asymmetry was revealed by testing the normality of right-minus-left values distribution using the Shapiro-Wilk test, which detected the absence of anti-symmetry (B: p = 0.80, C: p = 0.62, D: p = 0.97). Anti-symmetry and trait size were checked according to the protocol described by Kozlov and Niemelä [46]. Absence of directional asymmetry was confirmed by the zero mean value (B: p = 0.46, C: p = 0.35, D: p = 0.58). Thus, observed asymmetry was classified as FA. The Shapiro-Wilk test showed a normal distribution of the FA indices. The descriptive statistics of arithmetic mean, standard deviation, range, and the coefficient of variability were calculated. Next, the Student t-test was conducted to check statistically significant differences in terms of the fluctuating asymmetry factor in the studied dose of neomycin [48]. To provide a graphical presentation of results, a principal-component analysis (PCA) and cluster analysis were conducted.

Results of Larvae Fluctuation Asymmetry Measurements
The value of FA ranged from 0.055 to 0.198 for the first anal papillae and from 0.061 to 0.223 for the second ones. For both appendixes, the highest FA value was in the group treated with a concentration of 900 mg/L, and the lowest was in the control group (Table 5, Figure 3). The mean length of both appendixes was the highest in the control group (CA). The lowest value for the first anal papillae was in the groups treated with concentrations of 300 mg/L and 900 mg/L of neomycin, and for the second anal papillae the lowest value was observed in the group treated with a 600 mg/L dose of antibiotic.
The highest values of differences in length of anal papillae for both appendixes were in the 900 mg/L neomycin group, and the lowest value was for the control group ( Figure 3). Principal-component analysis (PCA) (Figure 4) and unweighted pair group method with arithmetic mean (UPGMA) dendrogram ( Figure 5) made using FA indices based on the two anal papillae showed that the 150 mg/L and 300 mg/L treatments formed one group, the second group included doses of neomycin 600 mg/L and 900 mg/L, and the control group was distinct. The Student t-test was conducted to check statistically significant differences in terms of the fluctuating asymmetry index in the studied dose of neomycin. In both appendixes, the control group differed statistically from 600 mg/L and 900 mg/L of neomycin treatments. For the second anal papillae, the groups treated with concentrations of 150 mg/L and 300 mg/L of neomycin differed statistically from groups treated with 600 mg/L and 900 mg/L of neomycin ( Figure 6). Principal-component analysis (PCA) (Figure 4) and unweighted pair group method with arithmetic mean (UPGMA) dendrogram ( Figure 5) made using FA indices based on the two anal papillae showed that the 150 mg/L and 300 mg/L treatments formed one group, the second group included doses of neomycin 600 mg/L and 900 mg/L, and the control group was distinct. The Student t-test was conducted to check statistically significant differences in terms of the fluctuating asymmetry index in the studied dose of neomycin. In both appendixes, the control group differed statistically from 600 mg/L and 900 mg/L of neomycin treatments. For the second anal papillae, the groups treated with concentrations of 150 mg/L and 300 mg/L of neomycin differed statistically from groups treated with 600 mg/L and 900 mg/L of neomycin ( Figure 6).

Results of Imago Fluctuation Asymmetry Measurements
The value of FA ranged from 0.063 to 0.224 for Vein B, from 0.044 to 0.237 for Vein C, and from 0.017 to 0.099 for Vein D. For Vein B, the highest FA value was for the 150 mg/L neomycin treatment and the lowest was for the 300 mg/L treatment. For Vein C, the highest FA value was for the 900 mg/L treatment and the lowest was for the 300 mg/L treatment. For Vein D, the highest FA value was for the 150 mg/L treatment and the lowest was for the 900 mg/L treatment (Table 6, Figure 7). The mean lengths of Veins B and C were the highest in the control group (CA) and the lowest in the 150 mg/L neomycin treatment group. The mean length of Vein D was also the highest in the control group (CA) but the lowest in the 900 mg/L neomycin treatment group. For all veins, the CA group had the longest veins, and the shortest veins were in the 150 mg/L treatment group for Veins B and C, whereas for Vein D the shortest veins were in the 900 mg/L treatment group. The highest values of differences in vein length for Veins B and D were for the concentration of 150 mg/L neomycin, and for Vein C were in the 900 mg/L treatment. The lowest value of differences was in the 300 mg/L group for Veins B and C, and in the 900 mg/L group for Vein D. Table 6. Descriptive statistics for the three studied veins of D. melanogaster wings. Figure 6. Student t-test applied between studied doses of neomycin (150-900 mg/L) and control group (CA) for two studied appendixes. AP1-the first pair of anal papillae, AP2-the second pair of anal papillae, * p < 0.05, ** p < 0.01, *** p < 0.001, ns-p > 0.05.

Results of Imago Fluctuation Asymmetry Measurements
The value of FA ranged from 0.063 to 0.224 for Vein B, from 0.044 to 0.237 for Vein C, and from 0.017 to 0.099 for Vein D. For Vein B, the highest FA value was for the 150 mg/L neomycin treatment and the lowest was for the 300 mg/L treatment. For Vein C, the highest FA value was for the 900 mg/L treatment and the lowest was for the 300 mg/L treatment. For Vein D, the highest FA value was for the 150 mg/L treatment and the lowest was for the 900 mg/L treatment (Table 6, Figure 7). The mean lengths of Veins B and C were the highest in the control group (CA) and the lowest in the 150 mg/L neomycin treatment group. The mean length of Vein D was also the highest in the control group (CA) but the lowest in the 900 mg/L neomycin treatment group. For all veins, the CA group had the longest veins, and the shortest veins were in the 150 mg/L treatment group for Veins B and C, whereas for Vein D the shortest veins were in the 900 mg/L treatment group. The highest values of differences in vein length for Veins B and D were for the concentration of 150 mg/L neomycin, and for Vein C were in the Insects 2020, 11, 378 9 of 17 900 mg/L treatment. The lowest value of differences was in the 300 mg/L group for Veins B and C, and in the 900 mg/L group for Vein D.   Principal-component analysis (PCA) (Figure 8) and UPGMA dendrogram (Figure 9) made using the FA indices based on the three veins showed that the CA, 300 mg/L, and 600 mg/L doses formed one group, and the 150 mg/L and 900 mg/L doses of neomycin were distinct from this group and from each other.
(b), and differences in vein length (c) for three studied veins.
Principal-component analysis (PCA) (Figure 8) and UPGMA dendrogram (Figure 9) made using the FA indices based on the three veins showed that the CA, 300 mg/L, and 600 mg/L doses formed one group, and the 150 mg/L and 900 mg/L doses of neomycin were distinct from this group and from each other.  The Student t-test was conducted to check statistically significant differences in terms of the fluctuating asymmetry factor in the studied concentrations of neomycin. The 150 mg/L concentration neomycin treatment differed from the 600 mg/L treatment with regard to Vein B. The 900 mg/L concentration differed from the 300 mg/L, 600 mg/L, and CA treatments with regard to Vein B, and because Vein D in the 900 mg/L treatment differed from that in CA ( Figure 10). The Student t-test was conducted to check statistically significant differences in terms of the fluctuating asymmetry factor in the studied concentrations of neomycin. The 150 mg/L concentration neomycin treatment differed from the 600 mg/L treatment with regard to Vein B. The 900 mg/L concentration differed from the 300 mg/L, 600 mg/L, and CA treatments with regard to Vein B, and because Vein D in the 900 mg/L treatment differed from that in CA ( Figure 10).
The Student t-test was conducted to check statistically significant differences in terms of the fluctuating asymmetry factor in the studied concentrations of neomycin. The 150 mg/L concentration neomycin treatment differed from the 600 mg/L treatment with regard to Vein B. The 900 mg/L concentration differed from the 300 mg/L, 600 mg/L, and CA treatments with regard to Vein B, and because Vein D in the 900 mg/L treatment differed from that in CA ( Figure 10).

Discussion
Fluctuating asymmetry (FA) is a measure of developmental instability expressed as a difference in accurate bilateral symmetry, caused by environmental stresses, and genetic problems during development. Environmental stressors that may induce FA in insects, usually include temperature extremes, oxygen deficiency, poor diet, and chemical pollution [49][50][51][52]. There is no information on the effects of antimicrobial substances such as antibacterials, antivirals, and antihelminthics on the fluctuating asymmetry (FA) of insects, including D. melanogaster. In many studies, the effect of Figure 10. Student t-test between studied doses of neomycin (150-900 mg/L) and control group (CA) for three studied veins; * p < 0.05, ** p < 0.01, ns-p > 0.05.

Discussion
Fluctuating asymmetry (FA) is a measure of developmental instability expressed as a difference in accurate bilateral symmetry, caused by environmental stresses, and genetic problems during development. Environmental stressors that may induce FA in insects, usually include temperature extremes, oxygen deficiency, poor diet, and chemical pollution [49][50][51][52]. There is no information on the effects of antimicrobial substances such as antibacterials, antivirals, and antihelminthics on the fluctuating asymmetry (FA) of insects, including D. melanogaster. In many studies, the effect of various pesticides on the FA of numerous insect groups has been investigated [3,4,[49][50][51][52]. Antipin and Imasheva [52] examined the effects of a chlorine-organic insecticide, endosulfane (thiodan), on phenotypic and genetic variation in wing length, thorax length, the number of orbital bristles, and the number of sternopleural bristles of D. melanogaster. Experiments carried out on Ceriagrion sp. larvae suggest that FA may be regarded as an indicator of pesticide stress. [50]. The usefulness of FA has also been demonstrated by investigating the effect of temperature and pesticide exposure on larval development in Copera annulata (Selys), among other studies. Insecticide treatment did not significantly affect the mortality of the larvae of damselfly, but the FA values of three traits decreased at lower concentrations, and then increased slowly with increased insecticide doses [49]. The results of these studies indicated that traits related to body size, expressed by the FA index, could be used as a suitable universal indicator of the environmental stress caused by chemical pollutants in populations of insects. A previous study reported the effects of chemicals such as arsenic and lead on FA in Drosophila [53]. It is known that nutrients available in the larval period are effective in all developmental stages of the insect [54]. For this purpose, the content and quality of artificial diets used for growing insects in the laboratory environment must be balanced [55][56][57][58][59]. Various antifungal, antibacterial, and antihelminthic agents are frequently used to improve diet quality. In this study, we used an asymmetry indicator to determine whether neomycin impaired the nutritional quality of the diet and led to nutritional stress. We checked whether the relationship between asymmetry parameters and the quality of the artificial diet could be determined. The effect of nutritional stress during development on FA has been addressed in a range of organisms [60,61]. Vijendravarma and coworkers (2011b) observed that larval malnutrition of Drosophila melanogaster had a significant effect on all wing traits at both the plastic and evolutionary levels. When raised on poor larval food, both control and selected populations had smaller wings than when raised on standard food [62]. On the other hand, experiments with Drosophila ananassae showed that the level of fluctuating asymmetry was similar in flies reared on poor and standard media [54]. This means that drawing general conclusions based on the use of fluctuating asymmetry as an indicator of nutritional stress should be done with caution.
Symmetry is important for the movement of an animal, reveals its physiological state, and it also plays an important role in reproduction, courtship, and reproductive behavior [63][64][65]. The results of our research indicated that neomycin causes significant stress for larvae and imagoes of D. melanogaster. Both examined life stages revealed significant differences in the FA indices of measured parameters. The results obtained for larvae appeared to be concentration-dependent; the FA index increased with the concentration of neomycin. Interestingly, there was also a distinct difference between effects caused by lower (150 and 300 mg/L) and higher (600 and 900 mg/L) concentrations of the antibiotic. Although the lowest FA values were observed in the control group, low levels of antibiotics did not lead to a significant increase in asymmetry. The highest FA index, observed in the group exposed to 900 mg/L neomycin, showed that the larvae had a pronounced nutritional stress response. The Pearson correlation coefficients for anal papillae (0.65 and 0.71) indicated a correlation between neomycin concentration and asymmetry. The strong effect for larvae may be due to their softer cuticle, and the inability to avoid the exposed environment (the larvae cannot emigrate from the substratum). Following this stage, they undergo processes that lead to fast growth, and development of tissues and organs that prepare them for the pupal stage. Therefore, any stress may significantly affect their morphology. For measurements of the length of anal papillae, the results were inconclusive. Although the addition of neomycin reduced their size compared to the controls, which was consistent with our previous observations [3], there was no clear tendency to decrease in length as the concentration of the antibiotic increased. Interestingly, the lowest and highest concentrations did not cause the most distinct effects. It is possible that the low concentrations of the antibiotic did not stimulate the insects' detoxifying pathways intensively enough, and neomycin remained in the body. On the other hand, the highest concentration might have been too high to be neutralized completely and some amounts of the antibiotic were still present in the environment.
The main role of the anal papillae is to regulate osmotic conditions. Neomycin has been described as a factor that affects concentration of potassium ions and blocks ryanodine RyR2 receptor [66][67][68]. The ryanodine receptors are also responsible for calcium storage within cells [69]. Therefore, due to altered concentrations of ions, this antibiotic may disturb osmotic conditions and, in consequence, lead to morphological malformations of anal papillae, which manifested in the measured FA. The ryanodine receptors are necessary for proper muscle physiology in fruit flies [70]. Hence, neomycin may significantly affect the development of insects and cause lethal and sublethal effects. These receptors are also present in insects' photoreceptor cells, including in Drosophila [71]. However, to our best knowledge, the effect of the antibiotic on the physiology of D. melanogaster cuticle-forming epithelium and anal papillae, in particular, have not been reported to date. Therefore, we plan to carry out further studies on the ultrastructure and morphology of organs and tissues exposed to antibiotics in the future.
The wing FA indices indicated a large but variable correlation. The Pearson correlation coefficients varied from 0.52 to 0.87. However, there was no strict rule and the veins did not show the same pattern of the FA value alterations. On the other hand, the mean length of veins was the highest for the control group, with neomycin-exposed groups showing lower values. That suggests that the antibiotic may decrease the size of some morphological/anatomical organs. This result was in line with results obtained by Ventrella et al. [3] and Chowański et al. [4], where decreased size and shortened time of development of D. melanogaster exposed to some natural substances was observed. It is possible that the insects speed up their development to the imaginal stage in order to decrease the time of exposure to toxic substances. In consequence, some processes and metabolic pathways may malfunction; the growth may be decreased, which results in asymmetry and the limited size of morphological features. However, due to the longer period of development, the presence of the pupal stage (when the insect does not feed and reorganizes its internal structures), the ability to move away from the neomycin-containing substratum, and a stronger cuticle, the malformations may be less significant than in the case of larvae. Additionally, the process of detoxification may be more intensive in imagoes than in larvae. In consequence, FA indices were not as significant as in the case of larvae. However, asymmetrical veins may significantly affect flight ability in imagoes, since veins play a crucial role in wing strengthening and in conditioning the surface of wings. Therefore, the increased asymmetry may affect imagoes in their behavior, locomotion, and perhaps reproductive success.
Although not lethal in the tested range of concentrations, neomycin may cause significant sublethal changes that manifest in FA. In consequence, it may lead to lower vitality of flies exposed to the antibiotic and limit the population size. It is possible that similar results would be observed not only in our model species, but also in others. The harmful activity of neomycin on insects was previously reported for Pimpla turionellae [72]. Most probably, the effect is caused by reactive oxygen species (ROS), which are produced during exposure to neomycin and other aminoglycosides [73,74]. Additionally, an increased level of O 2 − and H 2 O 2 , together with a decreased level of Cu, Zn superoxide dismutase was reported in cells exposed to gentamicin [75]. An imbalance between the level of ROS and antioxidant enzymes leads to oxidative stress within exposed organisms. If the cells are unable to achieve homeostasis, the imbalance will cause physiological abnormalities and morphological defects, which will lead to developmental errors and imperfections, which correlates with an increase in FA [76]. We plan to continue this research in the future, focusing on other antibiotics added to the media used to rear fruit flies, especially since Drosophila melanogaster is becoming increasingly recognized as a model organism in food and nutrition research.

Conclusions
Fluctuating asymmetry has been used as a bioinidicator of stress. In this paper, we measured FA in D. melanogaster larvae and imagoes reared on a neomycin-containing medium. Our findings showed that antibiotic may cause sublethal stress in D. melanogaster larvae and adults, and contribute significantly to morphological malformations of the exposed insects. This, in consequence, may lead to limited longevity, lower reproductive success, and decreased population size. The results suggest that neomycin, and perhaps other antibiotics, may cause significant physiological and developmental stress in insects and may possibly be used as insecticides in the future. However, such activity demands further research.
Author Contributions: Z.A. was responsible for investigation, data curation (scanning electron microscopy, larval measurements), methodology and supervision (took part in planning the research), analysis of results, and writing of the manuscript. E.C. was responsible for methodology, supervision, data curation, took part in analysis of results, and writing-original draft preparation. A.L.-W. was responsible for investigation and methodology (measurements and analysis of data), and writing-editing manuscript. R.G. formal analysis, writing-review and editing. E.B. was responsible for data curation, supervision, methodology, a part of writing manuscript. K.B. was responsible for supervision, data curation, maintaining insect culture and feeding experiments, methodology, a part of writing manuscript. All authors have read and agreed to the published version of the manuscript.