High Doses of Pesticides Induce mtDNA Damage in Intact Mitochondria of Potato In Vitro and Do Not Impact on mtDNA Integrity of Mitochondria of Shoots and Tubers under In Vivo Exposure

It is well known that pesticides are toxic for mitochondria of animals. The effect of pesticides on plant mitochondria has not been widely studied. The goal of this research is to study the impact of metribuzin and imidacloprid on the amount of damage in the mtDNA of potato (Solanum tuberosum L.) in various conditions. We developed a set of primers to estimate mtDNA damage for the fragments in three chromosomes of potato mitogenome. We showed that both metribuzin and imidacloprid considerably damage mtDNA in vitro. Imidacloprid reduces the rate of seed germination, but does not impact the rate of the growth and number of mtDNA damage in the potato shoots. Field experiments show that pesticide exposure does not induce change in aconitate hydratase activity, and can cause a decrease in the rate of H2O2 production. We can assume that the mechanism of pesticide-induced mtDNA damage in vitro is not associated with H2O2 production, and pesticides as electrophilic substances directly interact with mtDNA. The effect of pesticides on the integrity of mtDNA in green parts of plants and in crop tubers is insignificant. In general, plant mtDNA is resistant to pesticide exposure in vivo, probably due to the presence of non-coupled respiratory systems in plant mitochondria.


Introduction
Archeological and genetic research indicates that potato domestication (Solanum tuberosum L.) started about 8000 years ago on the border of present-day Peru and Bolivia [1]. According to FAOSTAT to 2019, the largest quantities of potatoes are produced in China . The widespread use of pesticides was the main reason for the increase in potato yields [3].
However, significant progress in yield protection was halted by the evolution of resistance against pesticides in pest populations. Synthetic and natural pesticides have thus become crucial; however, resources for pest control are non-renewable. The development of resistance led to the development and wider use of more and more pesticides [4]. In are more sensitive for pesticide exposure compared to intact mitochondria [36]. We can assume that in vivo plant mitochondria may be significantly more resistant to intact plant mitochondria and animal mitochondria under pesticides exposure.
The goal of this research is to develop a method for identifying mtDNA regions that are most susceptible to lesions induced by exogenous compounds; studying the effect of imidacloprid and metribuzin on mtDNA damage in vitro and in vivo; and evaluating this pesticide effect on the potato germination, shoot growth, rate of reactive oxygen species (ROS) production, and the activity of mitochondrial enzymes in the potato tubers. The results obtained will make it possible to judge about plant mitochondria susceptible to pesticide-induced damage in vitro compared to in vivo condition. mitochondria have a number of features which may mediate resistance to pesticides. Recently, we showed that alternative NADH dehydrogenases is a factor responsible for plant resistance to xenobiotics, such as mitochondria-targeted pesticides. The permeabilized mitochondria are more sensitive for pesticide exposure compared to intact mitochondria [36]. We can assume that in vivo plant mitochondria may be significantly more resistant to intact plant mitochondria and animal mitochondria under pesticides exposure. The goal of this research is to develop a method for identifying mtDNA regions that are most susceptible to lesions induced by exogenous compounds; studying the effect of imidacloprid and metribuzin on mtDNA damage in vitro and in vivo; and evaluating this pesticide effect on the potato germination, shoot growth, rate of reactive oxygen species (ROS) production, and the activity of mitochondrial enzymes in the potato tubers. The results obtained will make it possible to judge about plant mitochondria susceptible to pesticide-induced damage in vitro compared to in vivo condition.   We did not find a correlation between the number of H 2 O 2 -induced mtDNA of potato and frequency of occurrence of GTGR (r s = −0.6) or RTGR (r s = −0.31) sequences (where R is A or G). It should be noted that H 2 O 2 -induced mtDNA was statistically significant (p < 0.05) or close to statistically significant (the p value equals somewhere between 0.05 and 0.06) in all studied mtDNA fragments, except for ccmB (p = 0.58) (Figure 1). The amount of mtDNA damage induced by H 2 O 2 is characterized by an extremely high coefficient of variation (CV) (CV = 177 for ccmB, while CV ranged from 14 to 46 in other fragments).

Pesticides Effect on the Amount of mtDNA Damage In Vitro
Both metribuzin and imidacloprid damaged the potato mitochondrial genome in vitro for all studied mtDNA regions ( Figure 2A). Typically, more lesions were induced by the metribuzin. For the region encoding the cob gene (F (2,17) = 6.1623, p = 0.0097), more mtDNA lesions were induced by the imidacloprid (the median amount of damage per 10 kb of mtDNA was 4.39 (4.15; 4.58), p = 0.002). Metribuzin-induced damage of the cob gene-encoded fragment was not statistically significant (the median amount of damage per 10 kb of mtDNA was 4.01 (2.65; 4.28), p = 0.12 compared with control) ( Table 1).
There were changes in the amount of damage in the region encoding the cox1 gene (F (2,17) = 39.88, p = 0.00001). The median number of metribuzin-induced damage per 10 kb of mtDNA was 4.8 (4.28; 5.03), p = 0.0006, and the median number of imidacloprid-induced damage per 10 kb of mtDNA was 4.55 (4.03; 4.78), p = 0.002. A similar trend was observed for the ccmB and sdh3 genes located on the third chromosome. For ccmB (F (2,17) = 17.0728, p = 0.00009), the median number of metribuzin-induced damage per 10 kb of mtDNA was 4.72 (3.89; 5.23), p = 0.001, and the median number of imidacloprid-induced damage per 10 kb of mtDNA was 4.22 (3.63; 4.38), p = 0.005. For sdh3 (F (2,17) = 6.5538, p = 0.0078), the median number of metribuzin-induced damage per 10 kb of mtDNA was 4.96 (3.84; 5.10), p = 0.015, and the the median number of imidacloprid-induced damage per 10 kb of mtDNA was 4.24 (3.80; 4.51), p = 0.005 (Table 1). The p value corresponds to Tukey's post-hoc test. 1 p < 0.05 offers a comparison between the metribuzin and control groups. 11 p < 0.01 offers a comparison between the metribuzin and control groups. 111 p < 0.001 offers a comparison between the metribuzin and control groups. 22 p < 0.01 offers a comparison between the imidacloprid and control groups. The amount of mtDNA damage in potato shoots after treatment by metribuzin and imidacloprid, and co-treatment by metribuzin and imidacloprid. (C) The amount of mtDNA damage in potato tubers obtained at harvest in a field experiment. The control group did not treat any pesticides. The metribuzin group: seeds were treated herbicide Lazurit before planting only. Imidacloprid−/+ group: seeds were not treated before planting, and treated with Tanrek during the vegetation. The imidacloprid+/− group: seeds were treated with Tanrek before planting only. The imidacloprid+/+ group: seeds were treated with Tanrek both before planting and during the vegetation. The p value corresponds to Tukey's post-hoc test. p < 0.05 was considered to be statistically significant.

Effect of Pesticides on Seed Germination and Potato Growth under Controlled Conditions
We found that the control (1 and 2) and imidacloprid-treated (3 and 4) groups had the same technical germination-82.1%. Germination was recorded only on the 16th day in the imidacloprid-treated group, which indicated a delay in comparison with the plants of control group, which germinated even on the 2nd day. The start of germination curves The amount of mtDNA damage, induced by incubation of intact mitochondria for 30 min with respiratory substrate and pesticides (metribuzin or imidacloprid). The control mitochondria were incubated with respiratory substrate only. (B) The amount of mtDNA damage in potato shoots after treatment by metribuzin and imidacloprid, and co-treatment by metribuzin and imidacloprid. (C) The amount of mtDNA damage in potato tubers obtained at harvest in a field experiment. The control group did not treat any pesticides. The metribuzin group: seeds were treated herbicide Lazurit before planting only. Imidacloprid−/+ group: seeds were not treated before planting, and treated with Tanrek during the vegetation. The imidacloprid+/− group: seeds were treated with Tanrek before planting only. The imidacloprid+/+ group: seeds were treated with Tanrek both before planting and during the vegetation. The p value corresponds to Tukey's post-hoc test. p < 0.05 was considered to be statistically significant.

Effect of Pesticides on Seed Germination and Potato Growth under Controlled Conditions
We found that the control (1 and 2) and imidacloprid-treated (3 and 4) groups had the same technical germination-82.1%. Germination was recorded only on the 16th day in the imidacloprid-treated group, which indicated a delay in comparison with the plants of control group, which germinated even on the 2nd day. The start of germination curves in both groups is spasmodic. Then, the curves of germination reach a plateau. The germination maximum was equal after 35 days ( Figure 3A). The germination energy in the control group was 71.4%, in the imidacloprid-treated group-78.6%. lings were observed (F (1,54) = 1.5644, p = 0.21642) ( Figure 3C). Metribuzin treatment also did not affect the growth dynamics of potato seedlings. In particular, there were no differences in the length of seedlings on the 52nd day of the experiment (F (3,52) = 0.04674, p = 0.98645) ( Figure 3D). On the 52nd day of the experiment, there were no statistically significant differences in the mass of the root system of the seedlings (F (3,52) = 1.9549, p = 0.1322) (Figure 4).  Imidacloprid treatment affected seed germination. So, in the control group, 50% of the tubers sprouted on the fourth day. In the group treated with imidacloprid, 50% of the tubers sprouted on the 20th day ( Figure 3A). Differences in the size of potato seedlings were observed in the first weeks of the experiment, in particular, on the 20th day of the experiment (F (1,54) = 8.2663, p = 0.00577), when the median of the control seedlings was 0.5 cm (0.15; 0.9), and the average length of potato seedlings treated with imidacloprid was 0 cm (0; 0.4) ( Figure 3B). However, on the 38th day of the experiment (the day of metribuzin treatment of some experimental groups), no differences in the growth of seedlings were observed (F (1,54) = 1.5644, p = 0.21642) ( Figure 3C). Metribuzin treatment also did not affect the growth dynamics of potato seedlings. In particular, there were no differences in the length of seedlings on the 52nd day of the experiment (F (3,52) = 0.04674, p = 0.98645) ( Figure 3D). On the 52nd day of the experiment, there were no statistically significant differences in the mass of the root system of the seedlings (F (3,52) = 1.9549, p = 0.1322) ( Figure 4).

Effect of Pesticides on the Amount of mtDNA Damage in Potato Shoots
There was a statistically significant difference in the mean number of mtDNA damage in the potato shoots (F (3,652) = 7.6758, p = 0.00005). Tuckey's post-hoc test showed a statistically significant decrease in the number of lesions per 10 kb of mtDNA of shoots treated with imidacloprid (median 2.16 (2.33; 4.0)), compared with shoots treated with metribuzin (median 3.53 (1.06; 4.72), p = 0.038) ( Figure 2B).
For the mtDNA region encoding the 26s rRNA gene, a pesticide-induced decrease in the amount of mtDNA lesions was observed (F (3,120) = 8.38, p = 0.00004). There was a significant decrease in the amount of damage per 10 kb of mtDNA in metribuzin (

Effect of Pesticides on the Amount of mtDNA Damage in Potato Shoots
There was a statistically significant difference in the mean number of mtDNA damage in the potato shoots (F (3,652) = 7.6758, p = 0.00005). Tuckey's post-hoc test showed a statistically significant decrease in the number of lesions per 10 kb of mtDNA of shoots treated with imidacloprid (median 2.16 (2.33; 4.0)), compared with shoots treated with metribuzin (median 3.53 (1.06; 4.72), p = 0.038) ( Figure 2B).
For the mtDNA region encoding the 26s rRNA gene, a pesticide-induced decrease in the amount of mtDNA lesions was observed (F (3,120) = 8.38, p = 0.00004). There was a significant decrease in the amount of damage per 10 kb of mtDNA in metribuzin (  Table 2). The p value corresponds to Tukey's post-hoc test. 1 p < 0.05 offers a comparison between the metribuzin and control groups. 111 p < 0.001 offers a comparison between the metribuzin and control groups. 22 p < 0.01 offers a comparison between the imidacloprid and control groups. 333 p < 0.001 offers a comparison between the metribuzin-imidacloprid and control groups. 44 p < 0.01 offers a comparison between the imidacloprid and metribuzin groups. 444 p < 0.001 offers a comparison between the imidacloprid and metribuzin groups. 5 p < 0.05 offers a comparison between the imidacloprid-metribuzin and metribuzin groups. 555 p < 0.001 offers a comparison between the imidacloprid and metribuzin groups. 666 p < 0.001 offers a comparison between the imidacloprid-metribuzin and imidacloprid groups.

Effect of Pesticides on the Potato Growth and the Weight of the Crop in the Field Experiment
The median weight of the crop in the control group was 95 g (47.5; 149.5), as well as in the metribuzin group (median 100 g (69; 164)). There was a tendency to increase the crop weight in the group of potatoes, which was treated with imidacloprid. The median weight of crop was 222 g in the imidacloprid+/− group (42; 285), 120 g in the imidacloprid−/+ (109; 187.5), and 133.5 g in the imidacloprid+/+ (98.5; 229). However, there was no statistically significant difference (F (4, 31) = 0.9549, p = 0.4459) ( Figure 5B).

Figure 5.
The effect of pesticides on potato growth and crop weight in the field experiment. The control group was not treated with any pesticides. The metribuzin group: seeds were treated with herbicide Lazurit before planting only. The imidacloprid−/+ group: seeds were not treated before planting, and were treated with Tanrek during the vegetation. The imidacloprid+/− group: seeds were treated with Tanrek before planting only. The imidacloprid+/+ group: seeds were treated with Tanrek both before planting and during the vegetation. (A) The length of sprouts at day 52. (B) The weight of potato harvest at the end of experiment. There were no statistically significant differences.  Table 3).  Table 3). The p value corresponds to Tukey's post-hoc test. 22 p < 0.01 offers a comparison between the imidacloprid−/+ and control groups. 33 p < 0.01 offers a comparison between the imidacloprid +/− and control groups. 4 p < 0.05 offers a comparison between the imidacloprid+/+ and control groups. 7 p < 0.05 offers a comparison between the imidacloprid+/+ and metribuzin groups. 777 p < 0.001 offers a comparison between the imidacloprid+/+ and metribuzin groups. 888 p < 0.001 offers a comparison between the imidacloprid−/+ and imidacloprid+/− groups. 99 p < 0.01 offers a comparison between the imidacloprid−/+ and imidacloprid+/+ groups. 999 p < 0.001 offers a comparison between the imidacloprid−/+ and imidacloprid+/+ groups. *** p < 0.001 offers a comparison between the imidacloprid+/− and imidacloprid+/+ groups.  (Table 3).

Effect of Pesticides on Enzyme Activity in Potato Tubers
Statistically significant differences in activity values were found for succinate dehydrogenase (F (4,26) = 2.6392; p = 0.0566) ( Figure 6A). The enzyme activity in the imidacloprid−/+ group increased (median 0.00051 mmol/min/g of FW (0.000265; 0.000794)) compared to the imidacloprid+/− group (median 0.00178 mmol/min/g of FW (0.001305; 0.002381); p = 0.02512). No statistically significant changes were found for aconitate hydratase (F (4,26) = 0.5365; p = 0.7101) ( Figure 6B). Figure 6. Activity of succinate dehydrogenase (A) and aconitate hydratase (B) in mitochondria from potato tubers. The control group was not treated with any pesticides. The metribuzin group: seeds were treated with herbicide Lazurit before planting only. The imidacloprid−/+ group: seeds were not treated before planting, and were treated with Tanrek during the vegetation. The imidaclo-prid+/− group: seeds were treated with Tanrek before planting only. The imidacloprid+/+ group: seeds were treated with Tanrek both before planting and during the vegetation. The p value corresponds to Tukey's post-hoc test. p < 0.05 was considered to be statistically significant.

Effect of Pesticides on the Rate of H2O2 Production in Potato Tubers
Analysis of the data obtained during the assessment of the rate of H2O2 production without oligomycin showed statistically significant differences between groups (F (4,39) = 38.0792; p = 0.00001) ( Figure 7A). Mitochondria of the control group had maximal values of the rate of the H2O2 production (median 181 pmol/min/mg of protein (149; 196)). The minimal rate of the H2O2 production was in mitochondria from imidacloprid+/+ tubers (median 49 pmol/min/mg of protein (44; 55), which was significantly less compared with the control group (p = 0.000126). The metribuzin and imidacloprid−/+ groups also had a lower rate of H2O2 production compared with the control group (median 121 pmol/min/mg of protein (104; 129) (p = 0.000136) and 135 pmol/min/mg of protein (129; 145) (p = 0.002271), respectively). The rate of the H2O2 production in the imidacloprid+/− group was almost equal to the control group (median pmol/min/mg of protein 174 (158; 193) (p= 0.966632)) ( Figure 7A).
Similar results were obtained for mitochondria, which were inhibited by oligomycin (F (4,39) = 47.1607; p = 0.00001) ( Figure 7B). The maximal rate of the H2O2 production was observed for the control (median 171 pmol/min/mg of protein (159; 186)) and imidaclo-prid+/− (median 180 pmol/min/mg of protein (165; 194)) groups. For other groups, we found a statistically significant decrease in the rate of the H2O2 production (median for metribuzin group was pmol/min/mg of protein (104; 123), p = 0.000133. The median for imidacloprid−/+ was 140 pmol/min/mg of protein (131; 153), p = 0.006312, and the median for imidacloprid+/+ was 62 pmol/min/mg of protein (53; 71), p = 0.000126) ( Figure 7B). Figure 6. Activity of succinate dehydrogenase (A) and aconitate hydratase (B) in mitochondria from potato tubers. The control group was not treated with any pesticides. The metribuzin group: seeds were treated with herbicide Lazurit before planting only. The imidacloprid−/+ group: seeds were not treated before planting, and were treated with Tanrek during the vegetation. The imidacloprid+/− group: seeds were treated with Tanrek before planting only. The imidacloprid+/+ group: seeds were treated with Tanrek both before planting and during the vegetation. The p value corresponds to Tukey's post-hoc test. p < 0.05 was considered to be statistically significant.

Effect of Pesticides on the Rate of H 2 O 2 Production in Potato Tubers
Analysis of the data obtained during the assessment of the rate of H 2 O 2 production without oligomycin showed statistically significant differences between groups (F (4,39) = 38.0792; p = 0.00001) ( Figure 7A). Mitochondria of the control group had maximal values of the rate of the H 2 O 2 production (median 181 pmol/min/mg of protein (149; 196)). The minimal rate of the H 2 O 2 production was in mitochondria from imidacloprid+/+ tubers (median 49 pmol/min/mg of protein (44; 55), which was significantly less compared with the control group (p = 0.000126). The metribuzin and imidacloprid−/+ groups also had a lower rate of H 2 O 2 production compared with the control group (median 121 pmol/min/mg of protein (104; 129) (p = 0.000136) and 135 pmol/min/mg of protein (129; 145) (p = 0.002271), respectively). The rate of the H 2 O 2 production in the imidacloprid+/− group was almost equal to the control group (median pmol/min/mg of protein 174 (158; 193) (p= 0.966632)) ( Figure 7A).
Similar results were obtained for mitochondria, which were inhibited by oligomycin (F (4,39) = 47.1607; p = 0.00001) ( Figure 7B). The maximal rate of the H 2 O 2 production was observed for the control (median 171 pmol/min/mg of protein (159; 186)) and imidacloprid+/− (median 180 pmol/min/mg of protein (165; 194)) groups. For other groups, we found a statistically significant decrease in the rate of the H 2 O 2 production (median for metribuzin group was pmol/min/mg of protein (104; 123), p = 0.000133. The median for imidacloprid−/+ was 140 pmol/min/mg of protein (131; 153), p = 0.006312, and the median for imidacloprid+/+ was 62 pmol/min/mg of protein (53; 71), p = 0.000126) ( Figure 7B). Rate of H2O2 production in oligomycin-inhibited mitochondria of potato tubers. The control group was not treated with any pesticides. The metribuzin group: seeds were treated with herbicide Lazurit before planting only. The imidacloprid−/+ group: seeds were not treated before planting, and were treated with Tanrek during the vegetation. The imidacloprid+/− group: seeds were treated with Tanrek before planting only. The imidacloprid+/+ group: seeds were treated with Tanrek both before planting and during the vegetation. The p value corresponds to Tukey's post-hoc test. p < 0.05 was considered to be statistically significant.

Identification of mtDNA Regions That are More Prone to Oxidative Damage
It is well known that mtDNA is more prone to oxidative damage than nuclear DNA [37]. The proximity of mtDNA to the sites of oxidative phosphorylation and ROS production and the lack of the protection by histones are the main reasons why mtDNA is susceptible to damage [38]. For this reason, mtDNA damage is a popular marker of oxidative stress and an indicator of the level of degenerative processes during pathogenesis of various diseases in animal models [39][40][41] and humans [42]. There are data which show that pesticides can induce mtDNA damage in insects [10], mammals [43], and human's cell cultures [16]. At the same time, the ability of plant mtDNA to accumulate oxidative damage is practically not studied in comparison with animals.
There are a number of methods to estimate DNA damage. The classical methods, such as Southern blot and high-performance liquid chromatography, have a number of limitations. Long-range PCR has a number of advantages in measuring mtDNA damage, including its high sensitivity and ability to measure mtDNA damage in specific regions based on primer design [44]. We have previously shown that H2O2 induced significant damage to mtDNA; however, the distribution of damage is unequal along the molecule [45]. D-loop, the non-coding regulatory region of animal's mtDNA, is the most susceptibility region for somatic mutations [46] and oxidative damage [45,47]. Plants most likely employ multiple mechanisms to replicate the mtDNA due to the complex structure of the mitochondrial genome; nevertheless, such a structure in the D-loop appears to be missing [48]. Among the three studied potato chromosomes, we did not identify a region that would be most susceptible to H2O2-induced damage relative to other regions (Figure 1). It was previously shown that regions of DNA that are enriched in RTGR (where R is A or G) sequences are more sensitive to H2O2-induced oxidative damage due to the specific interactions of Fe 2+ ions with thymine in the purine-T-G-purine context which make the RTGR sequence more prone to the Fenton reaction [49,50]. Furthermore, the correlation production in oligomycin-inhibited mitochondria of potato tubers. The control group was not treated with any pesticides. The metribuzin group: seeds were treated with herbicide Lazurit before planting only. The imidacloprid−/+ group: seeds were not treated before planting, and were treated with Tanrek during the vegetation. The imidacloprid+/− group: seeds were treated with Tanrek before planting only. The imidacloprid+/+ group: seeds were treated with Tanrek both before planting and during the vegetation. The p value corresponds to Tukey's post-hoc test. p < 0.05 was considered to be statistically significant.

Identification of mtDNA Regions That Are More Prone to Oxidative Damage
It is well known that mtDNA is more prone to oxidative damage than nuclear DNA [37]. The proximity of mtDNA to the sites of oxidative phosphorylation and ROS production and the lack of the protection by histones are the main reasons why mtDNA is susceptible to damage [38]. For this reason, mtDNA damage is a popular marker of oxidative stress and an indicator of the level of degenerative processes during pathogenesis of various diseases in animal models [39][40][41] and humans [42]. There are data which show that pesticides can induce mtDNA damage in insects [10], mammals [43], and human's cell cultures [16]. At the same time, the ability of plant mtDNA to accumulate oxidative damage is practically not studied in comparison with animals.
There are a number of methods to estimate DNA damage. The classical methods, such as Southern blot and high-performance liquid chromatography, have a number of limitations. Long-range PCR has a number of advantages in measuring mtDNA damage, including its high sensitivity and ability to measure mtDNA damage in specific regions based on primer design [44]. We have previously shown that H 2 O 2 induced significant damage to mtDNA; however, the distribution of damage is unequal along the molecule [45]. D-loop, the non-coding regulatory region of animal's mtDNA, is the most susceptibility region for somatic mutations [46] and oxidative damage [45,47]. Plants most likely employ multiple mechanisms to replicate the mtDNA due to the complex structure of the mitochondrial genome; nevertheless, such a structure in the D-loop appears to be missing [48]. Among the three studied potato chromosomes, we did not identify a region that would be most susceptible to H 2 O 2 -induced damage relative to other regions (Figure 1). It was previously shown that regions of DNA that are enriched in RTGR (where R is A or G) sequences are more sensitive to H 2 O 2 -induced oxidative damage due to the specific interactions of Fe 2+ ions with thymine in the purine-T-G-purine context which make the RTGR sequence more prone to the Fenton reaction [49,50]. Furthermore, the correlation between ROS-induced DNA damage and the frequency of occurrence of GTGR (i.e., GTGG and GTGA) has been shown [51]. Then, confirmed these data was confirmed in mtDNA of mice and has been shown strong correlations between the H 2 O 2 -induced number of mtDNA damage and the frequency of occurrence of GTGR sequences in the mtDNA of brain (rs = 0.78) and liver (rs = 0.9) [45].In this research, we did not find a correlation between the number of H 2 O 2 -induced mtDNA damage of potato and the frequency of occurrence of GTGR (r s = −0.6) or RTGR (r s = −0.31) sequences. It should be noted that H 2 O 2 induced significant damage in all mtDNA fragments, except for ccmB (Figure 1). In this fragment, the amount of mtDNA damage induced by H 2 O 2 is characterized by an extremely high coefficient of variation compared with the other fragments. The ccmB gene in potato mtDNA is characterized by the fact that about 50% of ccmB transcripts cannot be translated into functional proteins [30]. It is likely that this region of mtDNA contains a component that affects the efficiency of not only RNA polymerase, but also DNA polymerase. For this reason, a high variation coefficient of the amount of mtDNA damage is observed when studying this mtDNA fragment.

Effect of Metribuzin and Imidacloprid on Potato mtDNA In Vitro
When evaluating the results of an in vitro experiment, it was found that the damage to mtDNA obtained after treatment with metribuzin and imidacloprid at a concentration of 200 µM is equivalent to the damage induced by treatment with H 2 O 2 at a concentration of 500 µM (Figure 2A, Table 1). We have previously shown similar results when treating mouse brain mitochondria with 500 µM H 2 O 2 and 1 µM rotenone pesticide, where rotenone caused even slightly more damage than H 2 O 2 [52]. It is believed that mtDNA damage is mainly initiated by ROS [53]. It is known that rotenone, which inhibits ETC complex I of animal mitochondria, causes overproduction of H 2 O 2 [54]. However, we have previously shown that imidacloprid and metribuzin did not increase the rate of ROS production in vitro [31], but they initiate serious damage. It is likely that the mechanism of DNArelated toxicity of pesticides is not only associated with the induction of oxidative stress. It is known that most pesticides are electrophilic substances, which ensures their high reactivity towards DNA. Pesticides induced adduct formation due to the binding with DNA. Reactions of DNA bases with pesticides or their metabolites alter its structure and prevent DNA replication [55]. So, we can assume that the mechanism of pesticide toxicity in vitro is not associated with ROS production, and it is most likely that pesticides directly interact with mtDNA and damage it.

Effect of Metribuzin and Imidacloprid on Potato mtDNA In Vivo
However, we have shown that pesticide treatment does not cause damage to mtDNA in the green parts of plant seedlings ( Figure 2B). In contrast, a pesticide-induced reduction in the number of lesions was observed in the mtDNA region encoding the 26s RNA gene. Imidacloprid treatment caused a decrease in the amount of mtDNA damage in the mtDNA region encoding the ccmB gene. However, the imidacloprid has an effect on the germination rate of potato seeds. In the control group, 50% of the tubers germinated on day 4, while in the imidacloprid-treated group they did so on day 20 ( Figure 3A). In the first days of the experiment, the growth rate of seedlings in the imidacloprid-treated group was lower than in the control; however, by the middle of the experiment, the experimental seedlings "caught up" to control seedlings. Previously it has been shown that the average imidacloprid-induced decline in normal germination of rice seed was 18.4% in 12 of the 15 cultivars or rice. However, there were generally no significant differences in shoot lengths and root system dry weights between control and pesticide-treated plants [56].These data are consistent with our results. We did not find statistically significant differences in shoot lengths on the 38th day of the experiment and root weight at the end of the experiment (Figures 3C and 4). It can be assumed that the treatment of seeds with an imidacloprid causes some decline, which affects the rate of germination. However, in green seedlings, these effects are leveled due to the presence of non-coupled respiration systems in plant mitochondria, primarily alternative oxidases (AOX). AOX can probably participate in the detoxification of pesticides by reducing the level of oxidative stress [57]. It is well known that AOX acts as a metabolic and signaling regulator and can maintain homeostasis, which is most important under stress conditions [58]. Pesticide exposure can be considered as severe abiotic stress for plants [59]. Nevertheless, the role of AOX in the regulation of plant response on pesticide toxicity needs further study.
However, there was research which showed the significant impact of the pesticide on potato germination and growth. Bocharov and Fomin showed that insecticides Kruiser (thiamethoxam) and Prestige (imidacloprid, pencycuron), in a mixture with the fungicide Maxim (fludioxonil), inhibited growth processes at the initial stage of ontogenesis (germination-seedlings) [60]. Pesticides reduced the length of regeneration seedlings and their roots (total and working surface of the root system) [60]. The high concentration (2-4 times higher than recommended dose) of some commonly used pesticides (including imidacloprid) decreased seed germination, biomass production, and photosynthetic pigments in tomato. However, at lower doses (quarter and half of the recommended dose), the pesticides had some stimulatory effects on growth as compared to the control [61]. These data do not fundamentally contradict our results, which showed that high doses of imidacloprid (three times higher than the recommended dose) negatively affect seed germination ( Figure 3A); however, we have not shown the effect on biomass production ( Figure 3C,D, Figures 4 and 5). In general, imidacloprid seems to have a positive effect on yield. Imidacloprid in a 90% and 30% average increase in the total grain yield of the wheat crop raised from seeds treated with the mixture relative to the corresponding untreated control. Imidacloprid seed treatment has important advantages over a conventional spray treatment as it reduces hazards and has a longer period of protection on the irrigated crop [62].
Potato tubers, unlike aerial parts of plants, do not have chloroplasts. Consequently, more energy is required from the mitochondria, and alternative respiratory pathways would interfere with the coupled respiration and ATP synthesis. Indeed, the low activity of alternative oxidase pathways in potato tuber mitochondria was previously shown [63,64]. Consequently, the mitochondria of potato tubers should be less protected from pesticideinduced damage. Indeed, we showed that the imidacloprid+/− and imidacloprid−/+ groups of crops had a minor increased number of mtDNA damage compared with control ( Figure 2C). However, the imidacloprid-induced amount of damage in the tubers of the potato crop is significantly less than that induced by damage in vitro ( Figure 2C, Table 1). It is noteworthy that the smallest amount of damage was observed in the imidacloprid+/+ group. The mitochondria of tubers of this group also produced the least amount of H 2 O 2 in comparison with other groups (Figure 7). This looks rather strange, since imidacloprid is an ATPase inhibitor [65]. It is known that the inhibition of ATPase, for example with oligomycin, leads to the hyperpolarization of mitochondria [66] which can lead to an increase in the rate of H 2 O 2 production [67]. Nevertheless, we observed an inverse relationship. The group that received the most imidacloprid treatment had the lowest rate of H 2 O 2 production ( Figure 7) and the smallest amount of mtDNA damage ( Figure 2C). The activity of aconitate hydratase, which may serve as a biosensor for oxidants [68], was also minimal in the imidacloprid+/+ group; however, the differences with other groups were statistically insignificant ( Figure 6B). This generally supports data that imidacloprid seed treatment reduces hazards and has a long period of protection on the irrigated crop [62,69,70], which manifested in the reduced level of oxidative stress (Figure 7) and the amount of mtDNA damage ( Figure 2C).
So, the effect of pesticides on the integrity of mtDNA in green parts of plants and in crop tubers is insignificant ( Figure 2B,C), despite the extremely high toxicity of the studied pesticides in vitro (Figure 2A). Given the effect of pesticides on the integrity of insect [9,10] and mammalian mtDNA [45,52], we can assume that plant mtDNA is more resistant to pesticides than animal mtDNA. This is consistent with the data of This is consistent with the data that herbicide-induced membrane damage (lipid peroxidation) in potato mitochondria was significantly lower than in rat mitochondria. Authors assumed different susceptibility of both kinds of mitochondria with different activities of protective antioxidant agents [71].In general, plants can detoxify the pesticides through multi-specific pathways. The reactions are primarily divided into three canonical phases, including phase I metabolism (pesticide metabolism by cytochrome P 450s , hydrolases, laccases), phase II conjugation (the metabolic pathway of S-thiols conjugates, detoxification by glutathione S-transferases and glycosyltransferases), and phase III transport (metabolism of pesticides by ATP binding cassette transporters) [72]. We recently showed that plant mitochondria, which oxidized NADH by complex I (not alternative NAD(P)H dehydrogenases), was more prone to pesticide-induced inhibition of the rate of O 2 consumption [36] than intact mitochondria, which oxidized of succinate [31]. Thus, respiration of plant mitochondria in "animal-like manner" (electron flow from complex I to complex IV) is a factor that makes mitochondria more sensitive to pesticides [36]. In sum, with current research, we can assume that alternative respiratory pathways can also mediate the plant mitochondria resistance to pesticide exposure.

Optimization of the Method for Measuring of the Number of mtDNA Damage
The primers panel was developed for six genes of potato mitochondrial DNA (Solanum tuberosum L., 1753). The mitochondrial genome of potato was formed by three chromosomes. We examined lesions for two genes from each chromosome: the sites of 26s rRNA and cox2 genes for the first chromosome (MN114537.1), cob and cox1 for the second (MN114538.1), and ccmB and sdh3 for the third chromosome (MN114539.1). The primer sequences are shown in Table 4.
The amount of mtDNA damage was assessed by long-range PCR using Encyclo polymerase. The method used was based on the assumption that DNA lesions, such as single-strand breaks, slow down accumulation of the PCR product due to the inhibition of elongation process. Therefore, the rate of PCR product accumulation would be inversely proportional to the number of damaged DNA molecules [44]. The primers were selected based on our earlier results, showing that the optimal length of the PCR product is about 2 kb [45]. To assess the amount of lesions, we compared the ∆Cq values of the control and experimental long fragments with the ∆Cq values of the control and experimental short fragments, which were used to normalize the amount of mtDNA. To optimize the panel of primers, mitochondria from potato tubers were isolated according to the method described below. Intact mitochondria were divided into two tubes. The first part of intact mitochondria (control) was incubated with a respiration substrate (10 mM succinate) only. The second part was additionally treated with 500 µM of H 2 O 2 for 30 min according to our previous research [45]. Then, mtDNA was isolated from both tubes with mitochondria for subsequent measurement of the amount of mtDNA by qPCR. The specificity of the PCR product was assessed by agarose gel electrophoresis in the 2% agarose gel with TAE buffer.

Designs of Experiments
Experiment 1. In vitro experiment ( Figure 8A). Intact mitochondria were obtained from potato tubers. The resulting final pellet was purified from supernatant and divided into three tubes. All tubes contained the respiration substrate (10 mM succinate). The isolation medium was added to the first tube (control). A medium containing 200 µM of metribuzin was added to the second tube. A medium containing 200 µM of imidacloprid was added to the third tube. It was previously shown that 200 µM of these pesticides caused the maximal in vitro effect for mitochondrial parameters [28]. The tubes were incubated at room temperature for 30 min, after which mtDNA was isolated from mitochondria for subsequent measurement of the amount of mtDNA damage by long-range PCR. Experiment 2. The effect of pesticides on the vegetative parts of potatoes. The potatoes were grown in 250 mL peat pots under controlled conditions at t = 24 • C, cycles of 16 h of light/8 h of darkness, and a relative humidity of at least 40%. The potato seeds were divided into four groups ( Figure 8B). The first group (control) (n = 14) was not exposed to pesticides either before planting or during the growing season.
The seeds of the second group (metribuzin) (n = 14) were not treated with pesticides before sowing; however, on the 38th day of the experiment, the vegetative parts of the plants were treated with the herbicide Lazurit. The recommended work solution of Lazurit was 40 L per 100 kg of tubers. The concentration of metribuzin was 700 g/kg. We treated 2 kg of tubers. Thus, 800 mL of Lazurit solution was needed, or 1.82 mg of metribuzin per 2 kg of tubers. We used Lazurit which was three times higher than the manufacturer's recommended concentration, or 5.46 mg of metribuzin per 2 kg, or 2.73 mg/kg of metribuzin.
The seeds of the third group (imidacloprid) (n = 14) were pretreated with the Tanrek insecticide. The recommended work solution of Tanrek was 1 L per 100 kg of tubers. The concentration of imidacloprid was 200 g/L. We treated of 2 kg of tubers. Thus, 20 mL of Tanrek solution was needed, or 4 mg of imidacloprid per 2 kg of tubers. We used Tanrek which was three times higher than the manufacturer's recommended concentration, or 12 mg of imidacloprid per 2 kg, or 6 mg/kg of imidacloprid. Vegetative parts of plants were not treated during the growing season.
The seeds of the fourth group (metribuzin-imidacloprid) (n = 14) were pretreated with a solution of the insecticide Tanrek at the rate of 6 mg of imidaclopride per 1 kg of seeds. Then, on the 38th day of the experiment, the vegetative parts of the plants were treated with the herbicide Lazurit at a concentration of 2.73 g/kg per 1 kg of seeds. The duration of the experiment from the day of landing was 52 days. The length of the seedlings was analyzed daily.
The technical germination of potato seedlings was calculated using the Formula (1) where G is germination, n is the number of germinated tubers, and N is the total number of tubers. The germination energy was calculated using the following Formula (2): where E is the germination energy, n is the number of tubers germinated after 1 /2 the term, and N is the total number of tubers.
On the 52nd day, seedlings were excised and frozen at t = −80 • C for subsequent DNA isolation. The mass of the root system was estimated. Experiment 3. Field experiment. The field experiment was carried out in 2021 at the town Biryuch, Belgorod region, Russia (50 • 39 53.2 N; 38 • 23 50.4 E). The mean monthly temperature was +20 • C and +11 • C (daytime and nighttime temperatures, respectively) in May, +25 • C and +16 • C in June, +29 • C and +20 • C in July, and +29 • C and +18 • C in August. The mean precipitation in May was 92 mm, 108 mm in June, 80 mm in July, and 39 mm in August. Potato tubers, which were used as seed tubers, were planted to a depth of 15 cm on 3 May. Landing was carried out in five groups, each of which has its own processing scheme ( Figure 8C,D).

3.
Imidacloprid−/+: seeds were not treated with an insecticide before planting, and treated during the vegetation. 4. Imidacloprid+/+: seeds were treated with an insecticide both before planting and during the vegetation.

5.
Metribuzin: seeds were treated with a metribuzin before planting.
The imidacloprid+/− and imidacloprid+/+ groups were treated with Tanrek before planting. The recommended work solution of Tanrek was 1 L per 100 kg of tubers. The concentration of imidacloprid was 200 g/L. We treated of 2 kg of tubers. Thus, 20 mL of Tanrek solution was needed, or 4 mg of imidacloprid per 2 kg of tubers. We used Tanrek which was three times higher than the manufacturer's recommended concentration, or 12 mg of imidacloprid per 2 kg, or 6 mg/kg of imidacloprid.
The metribuzin group was treated with Lazurit before planting. The recommended dose of Lazurit for preparing work solution was 3.33 g/L. The concentration of metribuzin was 700 g/kg of Lazurit. The final metribuzin concentration in work solution was 2.33 g/L. The recommended work solution of Lazurit was 3 L per 100 m 2 of soil. The area occupied by the group was 9 m 2 . We used Lazurit which was three times higher than the manufacturer's recommended concentration. Accordingly, it took 2.7 g of Lazurit containing 1.89 g of metribuzin, which was diluted in 500 mL of water. Moreover, the test area was sprayed. The control group did not treat any pesticides. The metribuzin group: seeds were treated herbicide Lazurit before planting only. The imidacloprid−/+ group: seeds were not treated before planting, and treated with Tanrek during the vegetation. The imidacloprid+/− group: seeds were treated with Tanrek before planting only. The imidacloprid+/+ group: seeds were treated with Tanrek both before planting and during the vegetation. p values correspond to Tukey's post-hoc test. p < 0.05 was considered to be statistically significant.
The imidacloprid+/− and imidacloprid+/+ groups were treated with Tanrek before planting. The recommended work solution of Tanrek was 1 L per 100 kg of tubers. The concentration of imidacloprid was 200 g/L. We treated of 2 kg of tubers. Thus, 20 mL of The control group did not treat any pesticides. The metribuzin group: seeds were treated herbicide Lazurit before planting only. The imidacloprid−/+ group: seeds were not treated before planting, and treated with Tanrek during the vegetation. The imidacloprid+/− group: seeds were treated with Tanrek before planting only. The imidacloprid+/+ group: seeds were treated with Tanrek both before planting and during the vegetation. p values correspond to Tukey's post-hoc test. p < 0.05 was considered to be statistically significant.
The imidacloprid−/+ and imidacloprid+/+ groups were treated with Tanrek during the vegetation on 18 July. The recommended work solution of Tanrek was 1 mL per 100 m 2 of soil or 0.2 g of imidacloprid per 100 m 2 . The area occupied by the group was 9 m 2 . We used Tanrek which was three times higher than the manufacturer's recommended concentration, or 3 mL/100 m 2 . Additionally, 270 µL of Tanrek (54 mg of imidacloprid) was diluted in 500 mL of water. Furthermore, the test area was sprayed.
The potato crop was harvested on 29 August. The mass of the harvested potato tubers was measured. Moreover, these tubers were used to isolate the mitochondria and to subsequently measure the H 2 O 2 production rate, enzyme activity, and mtDNA isolation.

Mitochondria Isolation
Mitochondria were isolated from potato tubers in a medium containing 225 mM of sucrose, 100 mM of mannitol, 20 mM of HEPES, 1 mM of EDTA, 1 mM of DTT, and1.5 mg/L of fatty acids free BSA, pH 7.4 (KOH). The tubers of potatoes were homogenized with Dounce homogenizer (Thermo Fisher Scientific, USA) in the medium in a 1:1 ratio at +4 • C. The homogenate was centrifuged for 5 min at 1000× g. The resulting supernatant was transferred into a clean centrifuge tube and centrifuged for 5 min at 4000× g. The resulting supernatant was separated from the pellet, transferred to the new centrifuge tube, and centrifuged for 15 min at 15,000× g. Then, mitochondria were purified by centrifugation in a 17% Percoll gradient for 30 min at 30,000× g with a centrifuge acceleration of 6. After centrifugation, the mitochondrial fraction located 1 cm from the bottom (3-4 mL) was collected, diluted with the medium, and centrifuged for 15 min at 18,000× g. The pellet was transferred to a 1.5 mL Eppendorf tube and centrifuged for 10 min at 16,000× g. The supernatant was removed and the pellet was resuspended in 150 µL of isolation medium.

DNA Isolation
DNA was isolated using the PROBA-GS kit (DNA-Technologies, Moscow, Russia). DNA was isolated from the green part of seedlings and intact potato mitochondria. Next, 1-2 g of potato seedlings were homogenized in 500 µL of "Lysis buffer" from the kit. For DNA isolation from intact mitochondria, the final pellet, which contained mitochondria, was dissolved in the 500 µL of "Lysis buffer" from the kit. Then, 20 µL of the "Sorbent" was added to the tube with the homogenate (for isolation DNA from seedling) or resuspended pellet (for isolation DNA from intact mitochondria). The tube was mixed on a vortex. The resulting mixture was incubated for 20 min at 50 • C, and centrifuged at 16,000× g for one minute. The supernatant was removed. Then, 200 µL of "Wash Solution #1" was added to the sorbent, the mixture was vortexed, and then centrifuged at 16,000× g for one minute. Similar steps were carried out with "Wash Solution #2" and " Wash Solution #3". After removing the supernatant, the pellet in an open tube was dried by incubation at 50 • C for 5 min. Next, 100 µL of the "Eluent solution" was added to the pellet, mixed on a vortex, and incubated for 5 min at 50 • C. The mixture was centrifuged at 16,000× g for one minute, and the supernatant from the eluted DNA was collected. The quality and purity of isolated DNA were estimated by electrophoresis using a 2% agarose gel in a TAE buffer.

Measuring of the mtDNA Damage Amount
To measure the amount of mtDNA lesions, the real-time qPCR method was used with a CFX96 TM thermocycler (Bio-Rad, Hercules, CA, USA). PCR reaction protocol: preliminary denaturation-3 min, 95 • C; denaturation-95 • C, 10 s; annealing-59 • C, 30 s; elongation-72 • C, 4 min 30 s; 40 cycles. The PCR reaction mixture per sample included 1X qPCRmix-HS SYBR mix, 0.5 µM of forward and reverse primers, 5 ng of DNA, and 14 µL of deionized water. To assess the amount of lesions, we compared the ∆Cq values of the control and experimental long fragments with the ∆Cq values of the control and experimental short fragments, which were used to normalize the amount of mtDNA. The amount of mtDNA damage was calculated at 10 kb according to Formula (3): where D 10 -the value of damage per 10,000 base pairs; E-the value of the PCR efficiency; ∆Sh = Cq of the short target fragment-Cq of the short control fragment; ∆L = Cq of the long target fragment-Cq of the long control fragment; and L f -the length of the fragment, bp.

Assays of Enzyme Activity
The medium for measuring the activity of aconitate hydratase contained the following components: 20 mM of sodium citrate; 50 mM of Tris-HCl (pH 7.8). The reaction was initiated by the addition of mitochondria. The increase in the optical density was detected at 240 nm due to the formation of the double bond in cis-aconitate [73]. The medium for measuring the activity of succinate dehydrogenase contained the following components: 10 mM of HEPES, 50 mM of Tris-HCl (pH 7.4); 40 mM of succinate; 50 µM of 2,6-DCPIP; 50 µM of PMS; and 20 µM of EDTA. The reaction was initiated by the addition of succinate. The succinate dehydrogenase activity was determined by a decrease in the absorption at 600 nm, caused by reduction in the artificial electron acceptor DCPIP [74]. The measurements were carried out on a Hitachi U-2900 spectrophotometer at a wavelength of 240 nm for aconitate hydratase and 600 nm for succinate dehydrogenase. The measurements were carried out in 1 mL of medium with the addition of 15 µM of mitochondria for 3 min. The enzyme activity was calculated on the basis that the coefficient of molar extinction for aconitate hydratase is 3.08 cm −1 * mM −1 and that the succinate dehydrogenase is 21 cm −1 * mM −1 . The enzyme activity is presented in mmol/min/g of fresh weight (FW).

Assays of H 2 O 2 Rate Production
The assay of H 2 O 2 production rate was recorded using Amplex Red Ultra fluorescent dye. The measurements were carried out on a Hitachi F-7000 spectrofluorometer at an extinction wavelength of 568 nm and an emission wavelength of 581 nm. The assay of H 2 O 2 production rate was carried out in an acrylic cuvette in the medium contained 225 mM of sucrose, 100 mM of mannitol, 20 mM of HEPES; 1 mM E of DTA; 1.5 mg/L of BSA free from fatty acids, pH 7.4; 1 mL of medium; 10 mM of succinate as a respiration substrate; 4 mM of KH 2 PO 4 ; four units of horseradish peroxidase; and 1 µM of Amplex Ultra Red. Then, 15 µL of mitochondria was added to the cuvette. Moreover, for an assay of maximal rate of H 2 O 2 production, 10 µM of oligomycin was added. Values of the rate of H 2 O 2 production were recorded for mitochondria, before oligomycin adding and after oligomycin adding. The concentration of H 2 O 2 was measured as the fluorescence intensity of the resorufin formed during the reaction [75]. The protein concentration was measured using the PierceTM BCA Protein Assay Kit (Thermo Fisher Scientific, USA).

Statistical Analysis
Statistical analysis was performed using the Statistica 12 software package (StatSoft. Inc, Tulsa, OK, USA). DNADamageCalculator (Voronezh, Russia) was used for the quantification of the number of mtDNA damage. The distribution normality was determined using the Shapiro-Wilk test. The data were represented as the median (Q1, Q3). The results were analyzed by one-way analysis of variance (ANOVA). Tukey's post-hoc test was used to determine the significance level.

Conclusions
Thus, our studies have shown that the use of pesticides, such as metribuzin and imidacloprid, does not damage potato mitochondria. At the same time, it should be taken into account that pesticides can accumulate in water and soil, and can harm animals and humans. Agricultural workers who mix, load, transport, and directly use pesticides are generally considered to be the group that will receive the most exposure and, therefore, are most at risk [76]. Expanding understanding of the mechanism of plant protection from pesticide-induced toxicity can help in the development of pharmacological methods for protecting humans and animals from the negative effects of pesticides. Together, this and our previous research [35] suggest that alternative respiratory pathways of plant mitochondria, which provide uncoupled respiration, may be responsible for plant resistance to pesticides. The pharmacological pathway for the uncoupling of mitochondrial respiration in animal mitochondria may be a perspective approach to help develop methods for protecting agricultural workers from the harmful effects of pesticides.