Abstract

Although the development of acaricide resistance in spider mites is a long-standing issue in agricultural fields, recent problems with acaricide resistance may be characterized by the development of complex- and/or multiresistance to acaricides in distinct classes. Such complexity of resistance is not likely to be a single mechanism. Pesticide resistance involves the microevolution of arthropod pests, and population genetics underlies the evolution. In this review, we address the genetic mechanisms of acaricide resistance evolution. We discuss genetic diversity and linkage of resistance genes, relationships between mite habitat and dispersal, and the effect of dispersal on population genetic structure and the dynamics of resistance genes. Finally, we attempt to present a comprehensive view of acaricide resistance evolution and suggest risks under globalization as well as possible approaches to managing acaricide resistance evolution or emergence.

1. Introduction

Safe and sufficient agricultural production is one of the most important global issues of the 21st century. Along the demands of food safety and minimizing environmental impacts, the idea of pest management is rapidly moving from chemical control toward sustainable technologies involving the exploitation of ecological functions. Nevertheless, insect pests and diseases cause huge yield losses in crops production, and pesticides have played an important role in mitigating yield loss of agricultural produce. Considering that the occurrence and geographical distribution of pests may change with global warming and economic globalization, pesticides will continue to be important in future agricultural production.

However, repeated pesticide use results in microevolution due to the selection pressure ultimately leading to resistance development. Marked development of acaricide resistance has been found in spider mites (Acari: Tetranychidae), which have short generation and high propagation potential. Since the 1990s, many populations of the two-spotted spider mite, Tetranychus urticae Koch, and the European red mite, Panonychus ulmi (Koch), have developed global resistance against many new acaricides. Among the top 20 resistant agricultural and medical arthropod pests, ranked by number of unique compounds in Arthropod Pesticide Resistance Database (http://www.pesticideresistance.org/), the most resistance has developed in T. urticae, whereas P. ulmi is the ninth most resistant [1].

Recent acaricide resistance problems are characterized by the development of complex- and/or multiresistance to acaricides in distinct classes. Pleiotropic effects may explain multiresistance in that a single key factor, for example, trans-acting factors, confers resistance to multiple insecticides involving distinct classes via a simultaneous increase in the action of several metabolic mechanisms involved in resistance [2]. Similarly, acaricide selection frequently confers cross resistance to acaricides in the distinct classes on spider mites, for example, T. urticae [3, 4] and P. ulmi [5]. However, the mechanisms of such resistance are largely unknown.

The complexity of acaricide resistance found today may not be explained by a single mechanism. Tetranychus urticae is a ubiquitous species that vigorously attacks an enormous number of plant species. Cross-resistance and the inheritance of resistance are not necessarily common to localities in this mite. These findings suggest that each selection process possibly leads to a different type of resistance mechanism, depending on the location and past selection history, and may yield different cross-resistance patterns.

Pesticide-resistance development is an example of microevolution in arthropod pests. Population genetics, which is affected by many factors including biological and ecological traits of arthropod pests, properties of pesticides, and pesticide application patterns [68], underlies the evolution. In this review, we address genetic mechanisms during the process of acaricide resistance evolution and discuss features of resistance genes and their genetic linkage, relationships between habitat and mite dispersal, and the effect of dispersal on population genetic structure and the dynamics of resistance genes to gain insight into spider mite resistance management in the global physical distribution environment. Finally, we comprehensively evaluate the evolution of acaricide resistance and suggest the global risks and management prospects.

2. Aspects and Features of Acaricide Resistance

Today’s acaricide resistances assume a new complicated aspect rather than the past; that is, genetic mode of resistance is regionally different. Etoxazole resistance is inherited maternally in Korean T. urticae populations [9], implying that it is caused by a mutation in cytoplasmic mitochondrial DNA. In contrast, in Japanese populations, the etoxazole resistance locus is located on a nuclear chromosome [10]. Hexythiazox resistance is under monogenic control in Australian populations [11], whereas more than one locus is involved in Japanese populations [10]. Similarly, chlorfenapyr resistance in Belgian populations is under polygenic control [12, 13], but Japanese populations might have a resistant allele at only one of the possible chlorfenapyr resistance loci (monogenic) [14]. Moreover, while clofentezine resistance confers a high level of cross-resistance to hexythiazox [15] but not to etoxazole in Australia, hexythiazox resistance probably confers etoxazole resistance in Japanese local populations [10, 16]. Cross-resistance has not been reported between etoxazole and hexythiazox [17, 18].

Monogenic resistance is favored under field selection regimes [1921], and acaricide resistance is typically monogenic [11, 2224], although some instances of polygenic resistance have been reported [2527]. However, each polygenic allele for chlorphenapyr and hexythiazox resistance may be capable of conferring substantial resistance [10, 12, 13]. Previous standard polygenic inheritance descriptions likely cannot adequately predict the evolutionary process of this type of inheritance in which several or a few resistance genes have additive and major effects on the resistant phenotype. In this regard, the LC50 for hexythiazox is notably higher for polygenic resistance in Japan ( 10,000 mg L-1) [10] than for monogenic resistance in Australia (48 mg L-1) [11], whereas similar LC50 values for chlorfenapyr were reported in a Japanese resistant population (2,130 mg L-1) [14] and a Belgian resistant population (2,939 mg L-1) [12]. Moreover, the potential for each resistance gene is individually associated with other distinctive resistance genes.

3. Genetic Linkage among Acaricide Resistance Genes

Most spider mite species have arrhenotokous parthenogenesis; males develop from unfertilized eggs (haploid), and females develop from fertilized eggs (diploid) [28]. This type of parthenogenesis may favor the fixation of fewer but more advantageous alleles [29, 30] but see [31, 32] such as acaricide resistance alleles under acaricide application. Fast fixation of advantageous alleles in male haploids may result in rapid fixation of linked alleles [33], as in the hitchhiking effect [34, 35].

To date, little is known about the genetic linkages among acaricide resistance genes. Uesugi et al. [14] revealed that T. urticae etoxazole and chlorfenapyr resistance genes are genetically linked to each other with a recombination rate of 14.8% (15.3 cM) using the three-point cross with malate-dehydrogenase allozymes [36] and located at 55.8 and 34.3 cM from loci accompanying etoxazole and chlorfenapyr resistance, respectively (Figure 1). They predicted the probability of hitchhiking effects between two closely linked acaricide resistance genes under acaricide selection, calling this phenomenon “apparent cross-resistance.’’ Although hexythiazox resistance is under monogenic control in Australian T. urticae populations, Asahara et al. [10] demonstrated that more than one locus is involved in hexythiazox resistance in a Japanese T. urticae population. Moreover, they found that one of the loci associated with hexythiazox resistance was tightly or completely linked to the etoxazole resistance locus (Figure 1). This suggests that the development of hexythiazox resistance possibly conferred etoxazole resistance in a Japanese local population of T. urticae that had never been exposed to etoxazole application [16].

Construction of a genetic linkage map is helpful to determine the effects of genetic relationships among loci associated with distinct acaricides, including the hitchhiking effects during the development of acaricide resistance, as Yan et al. [37] analyzed in yellow fever mosquitoes, Aedes aegypti (L.) (Diptera: Culicidae). Genetic markers, especially microsatellites, should be useful tools for genetic mapping of resistance genes. A substantial number of microsatellites have been cloned and analyzed in spider mites, including the citrus red mite Panonychus citri (McGregor) [38], T. urticae [3941], and the Kanzawa spider mite Tetranychus kanzawai Kishida [42]. Microsatellites have also been applied to kin species such as Tetranychus turkestani (Ugarov and Nikolskii) [43]. We conduct a linkage analysis of these microsatellites and have revealed two (12 loci) and three (14 loci) linkage groups in T. urticae and T. kanzawai, respectively (Uesugi et al. unpublished data, Matsuo et al. unpublished data). A linkage map with higher resolution should be complete in the near future.

4. Breeding Patch and Genetic Differentiation

Incidental factors such as gene flow (migration rate) and founder effects affect the dynamics of the local frequency of resistance genes in a metapopulation. Follett and Roderick [44] demonstrated that several resistance foci formed in a metapopulation through model analyses in which migration followed a two-dimensional stepping-stone model under the assumption that individuals move only to neighboring subpopulations in a latticed distribution and a percentage of the populations go extinct at each generation. Spider mites are tiny, wingless arthropods. They exhibit two major dispersal behaviors, crawling and aerial dispersal [4547]. Crawling is assumed to be used over short distances, such as between leaves within a plant, whereas aerial dispersal may be used for longer distance migration.

We addressed spider mite movement and breeding patches on a host plant because there is some evidence of sibmating in Tetranychus spider mites: a fertilized female deposits many eggs on the same leaf and very little migration occurs during development under uncrowded conditions [4852]. Hinomoto and Takafuji [53, 54] studied the genetic structure of a T. urticae population on strawberries based on phosphoglucoisomerase allozyme frequency [36]. Interestingly, they found the highest level of genetic differentiation, as measured by the fixation index (FST) [55], at the smallest habitat (leaflet) level. FST decreased with increasing the habitat size, that is, leaflet   leaf   plant [53, 54]. In contrast, the value of the inbreeding coefficient (FIS) [55] was the highest at the largest habitat level (subdivided section in the strawberry beds) and the lowest at the smallest habitat level (leaflets) [53]. In this regard, the Wahlund effect probably became larger in larger habitat sizes pooled for F-statistics analyses. Therefore, such genetic structure implies that although T. urticae females and males randomly mate on breeding patches (leaflets), gene flow among the patches is limited. The tendency to form breeding patches probably increases when population densities are low [50, 52] and such genetic structure may often cause genetic differentiation among breeding patches through random genetic drift and/or founder effects.

Adult P. citri females are rather mobile; their microhabitat is greenish twigs with foliage in citrus orchards and the females deposit their eggs on plural leaves [56]. Osakabe and Komazaki [57] focused on the genetic structure in a limited habitat area and performed several laboratory experiments with small citrus seedlings using esterase locus (α-Est1) alleles [58]. They introduced females from two P. citri strains, in which the α-Est1 locus was fixed to distinct alleles on the small citrus seedlings, and analyzed phenotype frequencies. As a result, although the inbreeding coefficients within a seedling decreased with population increase, the inbreeding coefficients remained substantially higher than zero even after the mite population had expanded enough to overlap the introduced colonies. Panonychus citri may couple randomly by nature, at least within the smallest breeding patch of a single leaf [57]. Despite their expected mobility [56], adult P. citri females oviposited on very few leaves, at least when the population density was low, and most individuals stayed on the leaves where they had been oviposited during the juvenile developmental stages in the experiments by Osakabe and Komazaki [57].

Adult males usually guard quiescent deutonymphal females and copulate with newly emerged females immediately after their last molts [59]. Once the first copulation for a virgin female has been completed, subsequent mating with other males is ineffective [60, 61]. Many offspring mature on the same leaf, at least during low population density. These spider mite habits might facilitate founder effect in each breeding patch and influence preservation (or fixation) of rare genes, such as acaricide resistance genes.

Under a condition of no acaicide selection, Uesugi et al. [62] demonstrated that genetic drift also works on the dynamics of an acaricide resistance allele; they found a remarkable decrease in mortality following an etoxazole application in two of 32 experimental subpopulations isolated for 13 generations, whereas the resistance did not develop in another 64 subpopulations that experienced limited gene flow.

5. Dispersal Depending on Habitat Structures and Environments

Dispersal behavior, colonization, and the resulting population structure of spider mites should strongly affect the distribution and spread of acaricide resistance genes [63]. Because of vigorous reproduction, spider mites cause leaf and plant deterioration, resulting in both food shortage and desiccation [64, 65]. Such unfavorable conditions motivate dispersal behaviors.

Dispersal, which is an important event for settlement of new habitats [63, 6668], may be affected by the surroundings of microhabitats that may vary depending on the agricultural system, for example, between woody and herbaceous plants and between greenhouse and field situations. This should result in genetic structure diversity of each local population and a consequent geographic distribution of alleles. Goka and Takafuji [69] found that Japanese T. urticae populations on fruit trees and roses could be subdivided into three geographical groups based on allozyme loci allele frequencies, whereas no particular geographical patterns were detectable in populations that inhabited annually herbaceous crops.

5.1. Dispersal and Genetic Differentiation in Greenhouses

Most local populations on herbaceous crops occur only temporarily. Thus a population increase may begin by invasion of a small number of mites, facilitating founder effects and genetic differentiation between local populations even if their host plants are adjacently located [69, 70]. Particularly under greenhouse conditions, population growth on herbaceous crops usually starts with a small number of immigrants. They may often experience an intensive bottleneck due to acaricide application in a closed place. This is also the case for greenhouse-grown perennial plants such as roses, although spider mites may maintain their populations on host plants throughout the year. Spider mites are unlikely to use aerial dispersal by means of wind-borne in a greenhouse due to the calm conditions. Windborne dispersal requires adequate wind velocity of 1.5 m/s or more [7173]. Instead, dispersal by crawling or passive transportation by human movements through agricultural practices is expected in greenhouses [40].

Under such restricted dispersal conditions, genetic differentiation is likely to occur among greenhouse colonies, especially at low population densities [50, 52]. Hinomoto and Takafuji [53] found significant genetic differentiation in a T. urticae population on greenhouse strawberries ( , ). Tsagkarakou et al. [70] investigated population structures in greenhouses in France, where roses, carnations, egg plants, and tomatoes were cultivated. They also found significant differentiation within a greenhouse ( ). In contrast, Navajas et al. [40] obtained low estimates of FST (0.008-0.09). Nevertheless, the latter findings may not contradict the former examples. Navajas et al. [40] sampled mites from areas incontrovertibly wider than the breeding colony ranges; for example, 15 eggplants were sampled from each site (25 m long transects), and two females per plant were used for the analysis. Such sampling strategies might obscure the effects of genetic differentiation among breeding patches on the population genetic structure data. However, their data demonstrated a general trend toward significant heterozygote deficiency [40], suggesting a Wahlund effect due to substantial isolation of the mite population into smaller breeding patches.

Uesugi et al. [74] analyzed the fine-scaled genetic structure of T. urticae exposed to more than one acaricide heavily applied to rose trees in a greenhouse. They set 18 consecutive quadrats (length: 1.2 m) in two beds (lengths: 22 and 27 m) in a greenhouse and analyzed the genetic structure using microsatellites. As a result, the level of genetic differentiation ( ) [75] among quadrats was higher in the bed with lower mite density (0.221–0.34) than in that with higher mite density (0.134–0.14) [74]. An analysis of molecular variance (AMOVA) revealed a temporal change in allelic frequencies in the bed with lower mite density ( ), but not in that with higher mite density [74]. These differences likely reflect the larger genetic drift at lower mite density; small population size could cause strong genetic drift, thereby decreasing genetic variation, which is indexed by allelic richness and heterozygosity, and increasing genetic differentiation among quadrat populations [76, 77]. Uesugi et al. [74] found that isolation-by-distance [78], calculated by comparing the FST/(1-FST) values with the spatial distance of each quadrat population [79], was significant in the bed with higher mite density (Mantel test, ) but not in that with lower mite density. A significant positive autocorrelation was detected by spatial autocorrelation analysis [80] only within a short range (2.4–3.6 m) of the bed with higher mite density, suggesting that migration was limited to a short range, whereas larger-scale dispersal, such as aerial dispersal, did not appear to contribute to the genetic structure [74].

Spider mites within a breeding patch are likely to mate randomly rather than assortatively, and an increase in population density effectively reduces FST values among patches [54, 57]. Nevertheless, fine-scale genetic structure is preserved under stable environments such as a perennial host plant in a greenhouse [74]. Yano [81] demonstrated collective dispersal and resulting aggregation at a new feeding site in T. urticae; crawling dispersers tended to follow the trails of former dispersers, which may be associated with the threads spun by spider mites. Therefore, the tendency for high FIS values in greenhouse spider mites may be attributed to the traits of spider mites forming breeding patches on small habitats, such as an individual leaf, and the limitation of immigration among patches. Consequently, gene flow between breeding patches is delayed even if mites are located on the same plant. Instead, founder effects and genetic drift promote genetic differentiation between breeding patches via fixation to specific alleles in each breeding patch.

5.2. Aerial Dispersal Mode and Gene Flow in Orchards and Citrus Groves

Spider mites aerially disperse by means of wind-borne or ballooning. Young gravid adult females are the dominant dispersers in both aerial dispersal means [82, 83]. However, species that exploit the wind-borne and ballooning probably achieve different arrival distances. In a practical investigation, dispersal of only 10–20 m occurred from blackberry into corn [66, 71]. Model analyses have suggested 5.5–12 m and 16–48 m dispersal from corn and apple, respectively [84]. In peanut field experiments by Boykin and Campbell [71], regardless of release height for aerial dispersal, most T. urticae caught by sticky traps were only 60 cm above the ground (the lowest traps) and catch numbers decreased as the traps were set higher (the highest traps were set 2.1 m above the ground). Therefore, as a general rule for the passive aerial ground dispersal, the number of arrivals per unit area probably declines steeply with increasing distance from the point of departure [85]. In contrast, in the field experiments by Fleschner et al. [82], ballooning O. punicae were captured at the top 60 cm of a 6 m (  cm) pole in an avocado grove, suggesting that they drifted upward after takeoff and that mites dispersing by ballooning migrate longer distances than do windborne mites. In fact, P. citri has been observed aerially migrating between adjacent citrus groves [82] and from distant citrus groves to Japanese pear orchards [86]. Almost imperceptible air currents are sufficient for mites to balloon [82].

Uesugi et al. [87] investigated the structures of T. urticae populations in two apple orchards using five microsatellite loci. In the first orchard, apple trees were planted in lines at intervals of 5–8 m, and their branches did not contact neighboring trees. In the second orchard, apple trees were planted at intervals of 2 and 7.5 m within and between lines, respectively. To exclude effects on gene flow from direct crawling dispersal between apple trees, Uesugi et al. sampled mites from apple trees along a diagonal transect in the second orchard. Ultimately, mites were collected from apple trees separated by distances of 10–24 m within about 100 m of the lines in these orchards. As a result, heterozygote deficiency (f) [75] and among trees were positive in these orchards ( and in the first orchard, and and in the second orchard). This suggests that there was a certain level of limitation to mite migration, and complete genetic homogenization was prevented within a scale of the studied sites. Nevertheless, analyses of genetic similarity and spatial location using dendrograms, isolation-by-distance, and multiallelic correlograms revealed no significant genetic differentiation among trees in these orchards. This may be explained that the level of gene flow is relatively high over long distances of 100 m, resulting in no genetic structure among mite populations in distant trees. Similarly, Grafton-Cardwell et al. [63] found that the frequencies of T. urticae and Tetranychus pacificus McGregor resistant to cyhexatin and dicofol were equivalent in on neighboring almonds and cotton, although cyhexatin had been applied only on almonds while dicofol had been applied only on cotton, indicating the occurrence of substantial gene flow between mite populations on these plants by aerial dispersal. Notwithstanding, significant isolation-by-distance was found in an analysis of wider geographic scales (populations were 10 m to more than 100 km apart) using allozymes in Greece [88]. Those suggest limitations for this mite to regionally disperse using the wind.

Genetic differentiation among geographically distant populations was also found in P. citri. Osakabe et al. [89] analyzed P. citri genetic structure hierarchically from the level between patches within a tree to between regions (Honshu and Kyushu islands of Japan), using the coefficient of gene differentiation (GST) [90, 91]. They found that GST tended to be smaller with decreasing area and distance between populations [89]. This geographic trend was consistent with the results of a protein analysis by Osakabe and Sakagami [92]. However, the diversity of P. citri populations within a local population and/or a microhabitat was different from that of T. urticae, which maintained frequent gene flow in a local population [87, 93]. The divergence within a grove and/or a tree clearly accounted for the total diversity of P. citri populations; GST values between trees within a grove and between patches within a tree were larger than those between groves. These results suggest that P. citri frequently migrates between citrus groves within a locality but not across wider areas. Moreover, higher diversity among trees (or patches) within groves than between groves implies that this mite migrates mainly by aerial drift (ballooning) between groves [82], whereas crawling dispersal is limited. Thus interbreeding between patches is infrequent and solely dependent on interpatch distance [57].

Consequently, the two modes of aerial dispersal in spider mites, wind-borne and ballooning, likely result in widely different population structures. Dispersion patterns of species exploiting the wind-borne might be explained using an aerodynamic model, such as the seed flux model [94, 95], as described by Jung and Croft [84]. These patterns may represent a practical decline in the number of arrivals per unit area with increasing distance from departure point [85]. The significance of microhabitat for P. citri genetic population structures was also demonstrated in marine organisms with reduced dispersal, such as polychaetes, which disperse only during the planktonic larval stage [96].

6. Effects of Dispersal and Habitat Stability on Acaricide Resistance Development

Patch size and persistence of host plants may largely influence population biology as well as gene flow between patches (colonies) and/or local populations. Additionally, variation in the initial frequency of pesticide resistance genes and intensity of selection force by pesticide application may result in different pesticide resistance development rates among local pest populations [97], which may need to be considered to determine the velocity of acaricide resistance development and prevalence [19, 98100].

As described above, T. urticae populations inhabiting stable woody plant habitats (fruit trees and roses) were geographically differentiated into genetically structured groups in Japan, but such differentiation did not occur in populations inhabiting annual herbaceous crops [69]. Goka [101] assessed resistance to two acaricides, pyridaben and fenpyroximate (used in Japan since 1991), in T. urticae populations collected from 48 deciduous fruit tree orchards and rose gardens and from 42 patches of herbaceous plants (some of which had been grown in greenhouses) from spring 1993 to autumn 1994. They found a rapid development of acaricide resistance in the woody plant populations. Only 7 and 8 of 48 populations showed 100% mortality after the application of pyridaben (200 mg l-1) and fenpyroximate (50 mg l-1), respectively. In contrast, resistant individuals were found in more than 80% of the populations; more than half of the tested individuals were alive in 24 and 10 populations after the application of pyridaben and fenpyroximate, respectively. Remarkably, no individual died by either acaricide in six populations (Figures 2(a) and 2(b)). In contrast, resistant individuals were found in a few populations after application of the acaricides, and more than 50% of individuals survived in an identical population (Figures 2(c) and 2(d); arrows). Different velocities of resistance development between these habitats were not caused by differences in the frequency of acaricide application (Figure 3). Instead, these phenomena may have occurred due to differences in mite population size which can be harbored in the habitat and habitat stability. Mite populations on woody plants are stable residents such that gene flow between them results in the strong likelihood of selection and accumulation of rare genes for resistance, whereas those on herbaceous crops are unstable transients such that rare genes would be unlikely to be retained for long periods [101]. Yet, resistance to the acaricides was obviously developed in a single population on herbaceous plants (Figure 3; arrows), which could have resulted from founder effects.

In the assessment by Osakabe et al. [89], P. citri showed higher genetic diversity in acaricide resistance genes within a grove and within a tree than between groves; a two-level nested analysis of variance revealed significant difference in fenpyroximate resistance among trees within groves ( ) but not among groves ( ). This may reflect the population structure revealed using gene markers, as described above. Similar trends also occurred in the distribution of etoxazole resistance [89]. A variation in fenpyroximate susceptibility was found even among patches within a tree. In a field study by Yamamoto et al. [102], strong spatial aggregation of the hexythiazox resistance gene in P. citri was identified both between trees and within a single tree in a citrus grove. This result suggests that acaricide resistance genes are spread based on the population structure, that is, limited crawling dispersal within foliage [57] resulting in between-patch diversity and frequent longer aerial dispersal by ballooning resulting in between-grove homogeneity.

A study of genetic differentiation of a T. urticae population on roses in a greenhouse estimated short dispersal distances and thus low spread velocity of acaricide resistance genes [74]. Nevertheless, milbemectin resistance genes, which had been distributed only partially in the population in July 2006, spread over the bed after only 2 months (Uesugi et al. unpublished data). During those 2 months, the roses were sprayed twice with milbemectin. The mechanism of such quick resistance development is unclear. Given that gene flow among patches within a bed did not occur frequently, an increase in the gene frequencies as effects of selection likely functioned to develop resistance. As Navajas et al. [40] noted, human-promoted passive transportation through agricultural practices might facilitate the quick spread of resistance genes.

7. Comprehensive Aspects of Acaricide Resistance Evolution

The initial frequency of resistance genes and the gene dominance are influential factors that determine the rate of resistance development [44]. In general, given a nonstructured population, mutated alleles associated with pesticide resistance likely exist at very low frequencies, and so in most of the resistance alleles are expected to be heterozygotes. Nevertheless, theoretical model experiments have revealed that in a metapopulation condition even lower-frequency alleles form concentrated foci in a set of subpopulations with certain qualifications for selection of pesticides, extinction of subpopulations, and migration [44]. With these qualifications, subpopulation extinction and subsequent migration and reconstruction possibly occur based on incidental changes in gene frequency through founder effects. This sometimes results in the construction of subpopulations in which the resistance allele frequencies are relatively higher, making selection by pesticides more effective for resistance development.

Data concerning genetic structure suggest that spider mite populations are metapopulations. Although Uesugi et al. [87] did not find differentiation within an apple orchard, there was significant heterozygote deficiency (significantly positive FIS value) in all trees, suggesting that Wahlund effects may have been caused by the division of the population into small breeding patches. Young gravid females are the main dispersers of spider mites when they experience high density [83]. Whether they disperse aerially or crawl, the dispersers may arrive alone on a new feeding site. These dispersal behaviors and restricted dispersal ability by crawling imply that founder effects frequently act on spider mite colonies, promoting the development of resistance foci through selection for acaricides. The resistance genes will then spread limitedly to neighboring colonies by gene flow or colony expansion by propagation. Moreover, the resistance genes may be spread from foci to other distant habitats by aerial dispersal in each species.

We have no data about the probability of new-linkage formation among resistance genes controlled by distinctive loci on the same chromosome by recombination in field populations. The frequent occurrence of resistance foci may increase the mating chances between mites bearing different resistance genes and might be followed by the occurrence of hitchhiking effects. Precise genetic experiments will be necessary to elucidate this issue.

8. Acaricide Resistance Management and Globalization

Species identification is a basic issue for biosecurity as injurious species must be identified to prevent their invasion and establishment in new areas [103105]. Eight spider mite species including T. urticae, T. kanzawai, P. citri, and P. ulmi are serious cosmopolitan pests and are now listed as nonquarantine pests in the Plant Protection Law of Japan. Therefore, the risk of introducing exotic resistance alleles to inland areas is undoubtedly increasing. Even if there are inland populations resistant to the same acaricide, an invasion of an exotic resistance gene may be a potential risk over the long term because an introduced resistance gene would not be identical to that of inland populations. Therefore, more studies on the diversity of acaricide resistance genes should be performed using genetic marker linkage maps. Moreover, methods to control the spread of such exotic resistance genes may be important to prevent the development of more complex resistance, and mechanisms of regional dispersal and resulting gene flow may have to be addressed.

9. Perspectives

Identifying a mechanism for the development of pesticide resistance is important for advancing pesticide resistance management for arthropod pests. In spider mites, past genetic and ecological studies have comprehensively suggested that the local concentration of resistance genes (increasing gene frequency in breeding patches) resulting from genetic diversity within habitats based on their biological traits and selection by acaricides and gene flow from selection sites to surroundings (local and/or regional spread of resistance) are the processes of acaricide resistance evolution.

Our study has revealed the consequences of habitat stability to acaricide resistance development [101] and implies that disturbing the expansion of incipient colonies is effective for restraining acaricide resistance development and avoiding complications with spider mites. Although natural enemies such as phytoseiid mites (Acari: Phytoseiidae) play the role of the disturber, their settlement and predation activity are important. We expect the development of new technologies for the management of natural enemies, such as artificial diet provisions [106]. Other factors possibly disturbing incipient spider mite colonies such as ultraviolet rays, which are lethal to spider mites [107], should also be exploited. Overall, the establishment of an integrated mite management strategy aimed at lowering spider mite population density may be the best way to manage acaricide resistance.