The role, magnitude, and molecular nature of trans -driven expression variation underlying the upregulation of detoxification genes in pesticide resistant arthropod populations has remained enigmatic. In this study, we performed expression quantitative trait locus (eQTL) mapping ( n = 458) between a pesticide resistant and a susceptible strain of the generalist herbivore and crop pest Tetranychus urticae . We found that a single trans eQTL hotspot controlled large differences in the expression of a subset of genes in different detoxification gene families, as well as other genes associated with host plant use. As established by additional genetic approaches including RNAi gene knockdown, a duplicated gene with a nuclear hormone receptor HR96-related ligand-binding domain was identified as causal for the expression differences between strains. The presence of a large family of HR96-related genes in T. urticae may enable modular control of detoxification and host plant use genes, facilitating this species’ known and rapid evolution to diverse pesticides and host plants.

The extreme adaptation potential of the generalist herbivore Tetranychus urticae (the two-spotted spider mite) to pesticides as well as diverse host plants has been associated with clade-specific gene expansions in known detoxifying enzyme families, and with extensive and rapid transcriptional responses. However, how this broad transcriptional potential is regulated remains largely unknown. Using a parental/F1 design in which four inbred strains were crossed to a common inbred strain, we assessed the genetic basis and inheritance of gene expression variation in T. urticae. Mirroring known phenotypic variation in the progenitor strains of the inbreds, we confirmed that the inbred strains we created were genetically distinct, varied markedly in pesticide resistance, and also captured variation in host plant fitness as is commonly observed in this species. By examining differences in gene expression between parents and allele-specific expression in F1s, we found that variation in RNA abundance was more often explained in trans as compared to cis, with the former associated with dominance in inheritance. Strikingly, in a gene ontology analysis, detoxification genes of the cytochrome P450 monooxygenase (CYP) family, as well as dioxygenases (DOGs) acquired from horizontal gene transfer from fungi, were specifically enriched at the extremes of trans-driven up- and downregulation. In particular, multiple CYPs and DOGs with broad substrate-specificities for pesticides or plant specialized compounds were exceptionally highly upregulated as a result of trans-regulatory variation, or in some cases synergism of cis and trans, in the most multi-pesticide resistant strains. Collectively, our findings highlight the potential importance of trans-driven expression variation in genes associated with xenobiotic metabolism and host plant use for rapid adaptation in T. urticae, and also suggests modular control of these genes, a regulatory architecture that might ameliorate negative pleiotropic effects.

Maize (Zea mays subsp. mays) yield loss from arthropod herbivory is substantial. While the basis of resistance to major insect herbivores has been comparatively well-studied in maize, less is known about resistance to spider mite herbivores, which are distantly related to insects and feed by a different mechanism. Two spider mites, the generalist Tetranychus urticae, and the grass-specialist Oligonychus pratensis, are notable pests of maize, especially during drought conditions. We assessed resistance (antibiosis) to both mites of 38 highly diverse maize lines, including several previously reported to be resistant to one or the other mite species. We found that line B96, as well as its derivatives B49 and B75, were highly resistant to T. urticae. In contrast, neither these three lines, nor any others included in our study, were notably resistant to the specialist O. pratensis. Quantitative trait locus (QTL) mapping with replicate populations from crosses of B49, B75, and B96 to susceptible B73 identified a QTL in the same genomic interval on chromosome 6 for T. urticae resistance in each of the three resistant lines, and an additional resistance QTL on chromosome 1 was unique to B96. Single-locus genotyping with a marker coincident with the chromosome 6 QTL in crosses of both B49 and B75 to B73 revealed that the respective QTL was large-effect; it explained ∼70% of the variance in resistance, and resistance alleles from B49 and B75 acted recessively as compared to B73. Finally, a genome-wide haplotype analysis using genome sequence data generated for B49, B75, and B96 identified an identical haplotype, likely of initial origin from B96, as the source of T. urticae resistance on chromosome 6 in each of the B49, B75, and B96 lines. Our findings uncover the relationship between intraspecific variation in maize defenses and resistance to its major generalist and specialist spider mite herbivores, and we identified loci for use in breeding programs and for genetic studies of resistance to T. urticae, the most widespread spider mite pest of maize.

Chemical control strategies are driving the evolution of pesticide resistance in pest populations. Understanding the genetic mechanisms of these evolutionary processes is of crucial importance to develop sustainable resistance management strategies. The acaricide pyflubumide is one of the most recently developed mitochondrial complex II inhibitors with a new mode of action that specifically targets spider mite pests. In this study, we characterize the molecular basis of pyflubumide resistance in a highly resistant population of the spider mite Tetranychus urticae. Classical genetic crosses indicated that pyflubumide resistance was incompletely recessive and controlled by more than one gene. To identify resistance loci, we crossed the resistant population to a highly susceptible T. urticae inbred strain and propagated resulting populations with and without pyflubumide exposure for multiple generations in an experimental evolution set-up. High-resolution genetic mapping by a bulked segregant analysis approach led to the identification of three quantitative trait loci (QTL) linked to pyflubumide resistance. Two QTLs were found on the first chromosome and centered on the cytochrome P450 CYP392A16 and a cluster of CYP392E6-8 genes. Comparative transcriptomics revealed a consistent overexpression of CYP392A16 and CYP392E8 in the experimental populations that were selected for pyflubumide resistance. We further corroborated the involvement of CYP392A16 in resistance by in vitro functional expression and metabolism studies. Collectively, these experiments uncovered that CYP392A16 N-demethylates the toxic carboxamide form of pyflubumide to a non-toxic compound. A third QTL coincided with cytochrome P450 reductase (CPR), a vital component of cytochrome P450 metabolism. We show here that the resistant population harbors three gene copies of CPR and that this copy number variation is associated with higher mRNA abundance. Together, we provide evidence for detoxification of pyflubumide by cytochrome P450s that is likely synergized by gene amplification of CPR.

Plant-herbivore interactions promote the generation and maintenance of both plant and herbivore biodiversity. The antagonistic interactions between plants and herbivores lead to host race formation: the evolution of herbivore types specializing on different plant species, with restricted gene flow between them. Understanding how ecological specialization promotes host race formation usually depends on artificial approaches, using laboratory experiments on populations associated with agricultural crops. However, evidence on how host races are formed and maintained in a natural setting remains scarce. Here, we take a multidisciplinary approach to understand whether populations of the generalist spider mite Tetranychus urticae form host races in nature. We demonstrate that a host race co-occurs among generalist conspecifics in the dune ecosystem of The Netherlands. Extensive field sampling and genotyping of individuals over three consecutive years showed a clear pattern of host associations. Genome-wide differences between the host race and generalist conspecifics were found using a dense set of SNPs on field-derived iso-female lines and previously sequenced genomes of T. urticae. Hybridization between lines of the host race and sympatric generalist lines is restricted by post-zygotic breakdown, and selection negatively impacts the survival of generalists on the native host of the host race. Our description of a host race among conspecifics with a larger diet breadth shows how ecological and reproductive isolation aid in maintaining intra-specific variation in sympatry, despite the opportunity for homogenization through gene flow. Our findings highlight the importance of explicitly considering the spatial and temporal scale on which plant-herbivore interactions occur in order to identify herbivore populations associated with different plant species in nature. This system can be used to study the underlying genetic architecture and mechanisms that facilitate the use of a large range of host plant taxa by extreme generalist herbivores. In addition, it offers the chance to investigate the prevalence and mechanisms of ecological specialization in nature.

The tomato russet mite, Aculops lycopersici, is among the smallest animals on earth. It is a worldwide pest on tomato and can potently suppress the host’s natural resistance. We sequenced its genome, the first of an eriophyoid, and explored whether there are genomic features associated with the mite’s minute size and lifestyle. At only 32.5 Mb, the genome is the smallest yet reported for any arthropod and, reminiscent of microbial eukaryotes, exceptionally streamlined. It has few transposable elements, tiny intergenic regions, and is remarkably intron-poor, as more than 80% of coding genes are intronless. Furthermore, in accordance with ecological specialization theory, this defense-suppressing herbivore has extremely reduced environmental response gene families such as those involved in chemoreception and detoxification. Other losses associate with this species’ highly derived body plan. Our findings accelerate the understanding of evolutionary forces underpinning metazoan life at the limits of small physical and genome size.

Water-stress commonly affects crops grown in arid and semi-arid regions. Apart from the direct impact of this abiotic stress on yield, a diverse community of herbivores can outbreak under these conditions, and plant responses to abiotic stress may alter plant defense responses that deter herbivores. Outbreaks of both generalist and specialist spider mites are strongly associated with hot and dry conditions in the field settings. To understand how water-stress impacts maize responses to spider mites, we conducted parallel greenhouse and field experiments with maize plants (B73 inbred line) subjected to optimal irrigation and water-stress conditions (50-60% and 5-10% volumetric water content (VWC), and 25-32% and 10-15% VWC, in the greenhouse and field, respectively). In addition to recording B73 responses to water-stress alone, we measured the population growth of the generalist twospotted spider mite (Tetranychus urticae, TSM) and the specialist Banks grass mite (Oligonychus pratensis, BGM) on optimally watered and water-stressed plants. We also measured plant defense protein activities (peroxidase (POD), polyphenol oxidase (PPO), chitinase (CHI) and trypsin inhibitor (TI)) at 1, 3, and 7 days post mite introduction for each irrigation treatment. For B73 plants exposed to water-stress, we observed increases in leaf temperature, leaf water potential, POD activity, as well as decreases in stomatal conductance and stem height. Populations of both mite species increased more rapidly on water-stressed B73 plants. While optimally irrigated B73 plants responded with similar plant defense activity to both mite species, combinations of plant water-stress and TSM herbivory resulted in increases in CHI and TI activity that were not observed for the respective treatments with BGM. Collectively, our results highlight a role for species-specific factors, possibly associated with herbivore host plant breadth, in impacting plant responses to herbivory in combination with an abiotic stress.

BACKGROUND:Vector control is the main intervention in malaria control and elimination strategies. However, the development of insecticide resistance is one of the major challenges for controlling malaria vectors. Anopheles arabiensis populations in Ethiopia showed resistance against both DDT and the pyrethroid deltamethrin. Although an L1014F target-site resistance mutation was present in the voltage gated sodium channel of investigated populations, the levels of resistance indicated the presence of additional resistance mechanisms. In this study, we used genome-wide transcriptome profiling by RNAseq to assess differentially expressed genes between three deltamethrin and DDT resistant An. arabiensis field populations - Asendabo, Chewaka and Tolay - and two susceptible strains - Sekoru and Mozambique.RESULTS:Both RNAseq analysis and RT-qPCR showed that a glutathione-S-transferase, gstd3, and a cytochrome P450 monooxygenase, cyp6p4, were significantly overexpressed in the group of resistant populations compared to the susceptible strains, suggesting that the enzymes they encode play a key role in metabolic resistance against deltamethrin or DDT. Furthermore, a gene ontology enrichment analysis showed that expression changes of cuticle related genes were strongly associated with insecticide resistance. Although this did not translate in increased thickness of the procuticle, a higher cuticular hydrocarbon content was observed in a resistant population.CONCLUSION:Our transcriptome sequencing of deltamethrin and DDT resistant An. arabiensis populations from Ethiopia suggests non-target site resistance mechanisms and paves the way for further investigation of the role of cuticle composition in insecticide resistance of malaria vectors. © 2019 Society of Chemical Industry.

Bulked segregant analysis (BSA) is a cross-based method for genetic mapping in sexually reproducing organisms. The method's use of bulked (pooled) samples markedly reduces the genotyping effort associated with traditional linkage mapping studies. Further, it can be applied to species with life histories or physical attributes (as for micro-insects) that render genetic mapping with other methods impractical. Recent studies in both insects and mites have revealed that advanced BSA experimental designs can resolve causal loci to narrow genomic intervals, facilitating follow-up investigations. As high-quality genomes become more widely available, BSA methods are poised to become an increasingly important tool for the rapid mapping of both monogenic and polygenic traits in diverse arthropod species.

Pesticide resistance arises rapidly in arthropod herbivores, as can host plant adaptation, and both are significant problems in agriculture. These traits have been challenging to study as both are often polygenic and many arthropods are genetically intractable. Here, we examined the genetic architecture of pesticide resistance and host plant adaptation in the two-spotted spider mite, Tetranychus urticae, a global agricultural pest. We show that the short generation time and high fecundity of T. urticae can be readily exploited in experimental evolution designs for high-resolution mapping of quantitative traits. As revealed by selection with spirodiclofen, an acetyl-CoA carboxylase inhibitor, in populations from a cross between a spirodiclofen-resistant and a spirodiclofen-susceptible strain, and which also differed in performance on tomato, we found that a limited number of loci could explain quantitative resistance to this compound. These were resolved to narrow genomic intervals, suggesting specific candidate genes, including acetyl-CoA carboxylase itself, clustered and copy variable cytochrome P450 genes, and NADPH cytochrome P450 reductase, which encodes a redox partner for cytochrome P450s. For performance on tomato, candidate genomic regions for response to selection were distinct from those responding to the synthetic compound and were consistent with a more polygenic architecture. In accomplishing this work, we exploited the continuous nature of allele frequency changes across experimental populations to resolve the existing fragmented T. urticae draft genome to pseudochromosomes. This improved assembly was indispensable for our analyses, as it will be for future research with this model herbivore that is exceptionally amenable to genetic studies.

While substantial progress has been made in understanding defense responses of cereals to insect herbivores, comparatively little is known about responses to feeding by spider mites. Nevertheless, several spider mite species, including the generalist Tetranychus urticae and the grass specialist Oligonychus pratensis, cause damage on cereals such as maize and wheat, especially during drought stress. To understand defense responses of cereals to spider mites, we characterized the transcriptomic responses of maize and barley to herbivory by both mite species, and included a wounding control against which modulation of defenses could be tested. T. urticae and O. pratensis induced highly correlated changes in gene expression on both maize and barley. Within 2 h, hundreds of genes were upregulated, and thousands of genes were up- or downregulated after 24 h. In general, expression changes were similar to those induced by wounding, including for genes associated with jasmonic acid biosynthesis and signaling. Many genes encoding proteins involved in direct defenses, or those required for herbivore-induced plant volatiles, were strongly upregulated in response to mite herbivory. Further, biosynthesis genes for benzoxazinoids, which are specialized compounds of Poaceae with known roles in deterring insect herbivores, were induced in maize. Compared to chewing insects, spider mites are cell content feeders and cause grossly different patterns of tissue damage. Nonetheless, the gene expression responses of maize to both mite herbivores, including for phytohormone signaling pathways and for the synthesis of the benzoxazinoid 2-hydroxy-4,7-dimethoxy-1,4-benzoxazin-3-one glucoside, a known defensive metabolite against caterpillars, resembled those reported for a generalist chewing insect, Spodoptera exigua. On maize plants harboring mutations in several benzoxazinoid biosynthesis genes, T. urticae performance dramatically increased compared to wild-type plants. In contrast, no difference in performance was observed between mutant and wild-type plants for the specialist O. pratensis. Collectively, our data provide little evidence that maize and barley defense responses differentiate herbivory between T. urticae and O. pratensis. Further, our work suggests that the likely route to specialization for O. pratensis involved the evolution of a robust mechanism to cope with the benzoxazinoid defenses of its cereal hosts.

Synergists can counteract metabolic insecticide resistance by inhibiting detoxification enzymes or transporters. They are used in commercial formulations of insecticides, but are also frequently used in the elucidation of resistance mechanisms. However, the effect of synergists on genome-wide transcription in arthropods is poorly understood. In this study we used Illumina RNA-sequencing to investigate genome-wide transcriptional responses in an acaricide resistant strain of the spider mite Tetranychus urticae upon exposure to synergists such as S, S, S-tributyl phosphorotrithioate (DEF), diethyl maleate (DEM), piperonyl butoxide (PBO) and cyclosporin A (CsA). Exposure to PBO and DEF resulted in a broad transcriptional response and about one third of the differentially expressed genes (DEGs), including cytochrome P450 monooxygenases and UDP-glycosyltransferases, was shared between both treatments, suggesting common transcriptional regulation. Moreover, both DEF and PBO induced genes that are strongly implicated in acaricide resistance in the respective strain. In contrast, CsA treatment mainly resulted in downregulation of Major Facilitator Superfamily (MFS) genes, while DEGs of the DEM treatment were not significantly enriched for any GO-terms.

Carotenoids underlie many of the vibrant yellow, orange, and red colors in animals, and are involved in processes ranging from vision to protection from stresses. Most animals acquire carotenoids from their diets because de novo synthesis of carotenoids is primarily limited to plants and some bacteria and fungi. Recently, sequencing projects in aphids and adelgids, spider mites, and gall midges identified genes with homology to fungal sequences encoding de novo carotenoid biosynthetic proteins like phytoene desaturase. The finding of horizontal gene transfers of carotenoid biosynthetic genes to three arthropod lineages was unprecedented; however, the relevance of the transfers for the arthropods that acquired them has remained largely speculative, which is especially true for spider mites that feed on plant cell contents, a known source of carotenoids. Pigmentation in spider mites results solely from carotenoids. Using a combination of genetic approaches, we show that mutations in a single horizontally transferred phytoene desaturase result in complete albinism in the two-spotted spider mite, Tetranychus urticae, as well as in the citrus red mite, Panonychus citri. Further, we show that phytoene desaturase activity is essential for photoperiodic induction of diapause in an overwintering strain of T. urticae, consistent with a role for this enzyme in provisioning provitamin A carotenoids required for light perception. Carotenoid biosynthetic genes of fungal origin have therefore enabled some mites to forgo dietary carotenoids, with endogenous synthesis underlying their intense pigmentation and ability to enter diapause, a key to the global distribution of major spider mite pests of agriculture.

The extent to which epigenetic variation affects complex traits in natural populations is not known. We addressed this question using transcriptome and DNA methylation data from a sample of 135 sequenced A. thaliana accessions. Across individuals, expression was significantly associated with cis-methylation for hundreds of genes, and many of these associations remained significant after taking SNP effects into account. The pattern of correlations differed markedly between gene body methylation and transposable element methylation. The former was usually positively correlated with expression, and the latter usually negatively correlated, although exceptions were found in both cases. Finally, we developed graphical models of causality that adapt to a sample with heavy population structure, and used them to show that while methylation appears to affect gene expression more often than expression affects methylation, there is also strong support for both being independently controlled. In conclusion, although we find clear evidence for epigenetic regulation, both the number of loci affected and the magnitude of the effects appear to be small compared to the effect of SNPs.

Richard Mott, Simon Myers and colleagues present a new imputation method, STITCH, which does not require genotyping arrays or high-quality reference panels. They use STITCH to accurately impute genotypes in both outbred laboratory mice and a sample human population directly from low-coverage (<2×) sequencing data. Inexpensive genotyping methods are essential for genetic studies requiring large sample sizes. In human studies, array-based microarrays and high-density haplotype reference panels allow efficient genotype imputation for this purpose. However, these resources are typically unavailable in non-human settings. Here we describe a method (STITCH) for imputation based only on sequencing read data, without requiring additional reference panels or array data. We demonstrate its applicability even in settings of extremely low sequencing coverage, by accurately imputing 5.7 million SNPs at a mean r2 value of 0.98 in 2,073 outbred laboratory mice (0.15× sequencing coverage). In a sample of 11,670 Han Chinese (1.7× coverage), we achieve accuracy similar to that of alternative approaches that require a reference panel, demonstrating that our approach can work for genetically diverse populations. Our method enables straightforward progression from low-coverage sequence to imputed genotypes, overcoming barriers that at present restrict the application of genome-wide association study technology outside humans.

The two-spotted spider mite Tetranychus urticae is an extremely polyphagous crop pest. Alongside an unparalleled detoxification potential for plant secondary metabolites, it has recently been shown that spider mites can attenuate or even suppress plant defenses. Salivary constituents, notably effectors, have been proposed to play an important role in manipulating plant defenses and might determine the outcome of plant-mite interactions. Here, the proteomic composition of saliva from T. urticae lines adapted to various host plants—bean, maize, soy, and tomato—was analyzed using a custom-developed feeding assay coupled with nano-LC tandem mass spectrometry. About 90 putative T. urticae salivary proteins were identified. Many are of unknown function, and in numerous cases belonging to multimembered gene families. RNAseq expression analysis revealed that many genes coding for these salivary proteins were highly expressed in the proterosoma, the mite body region that includes the salivary glands. A subset of genes encoding putative salivary proteins was selected for whole-mount in situ hybridization, and were found to be expressed in the anterior and dorsal podocephalic glands. Strikingly, host plant dependent expression was evident for putative salivary proteins, and was further studied in detail by micro-array based genome-wide expression profiling. This meta-analysis revealed for the first time the salivary protein repertoire of a phytophagous chelicerate. The availability of this salivary proteome will assist in unraveling the molecular interface between phytophagous mites and their host plants, and may ultimately facilitate the development of mite-resistant crops. Furthermore, the technique used in this study is a time- and resource-efficient method to examine the salivary protein composition of other small arthropods for which saliva or salivary glands cannot be isolated easily.

To understand the population genetics of structural variants (SVs), and their effects on phenotypes, we developed an approach to mapping SVs, particularly transpositions, segregating in a sequenced population, and which avoids calling SVs directly. The evidence for a potential SV at a locus is indicated by variation in the counts of short-reads that map anomalously to the locus. These SV traits are treated as quantitative traits and mapped genetically, analogously to a gene expression study. Association between an SV trait at one locus and genotypes at a distant locus indicate the origin and target of a transposition. Using ultra-low-coverage (0.3x) population sequence data from 488 recombinant inbred Arabidopsis genomes, we identified 6,502 segregating SVs. Remarkably, 25% of these were transpositions. Whilst many SVs cannot be delineated precisely, PCR validated 83% of 44 predicted transposition breakpoints. We show that specific SVs may be causative for quantitative trait loci for germination, fungal disease resistance and other phenotypes. Further we show that the phenotypic heritability attributable to sequence anomalies differs from, and in the case of time to germination and bolting, exceeds that due to standard genetic variation. Gene expression within SVs is also more likely to be silenced or dysregulated. This approach is generally applicable to large populations sequenced at low-coverage, and complements the prevalent strategy of SV discovery in fewer individuals sequenced at high coverage.

Background: The European red mite, Panonychus ulmi, is among the most important mite pests in fruit orchards, where it is controlled primarily by acaricide application. However, the species rapidly develops pesticide resistance, and the elucidation of resistance mechanisms for P. ulmi has not kept pace with insects or with the closely related spider mite Tetranychus urticae. The main reason for this lack of knowledge has been the absence of genomic resources needed to investigate the molecular biology of resistance mechanisms.Results: Here, we provide a comprehensive strand-specific RNA-seq based transcriptome resource for P. ulmi derived from strains susceptible and resistant to the widely used acaricide spirodiclofen. From a de novo assembly of the P. ulmi transcriptome, we manually annotated detoxification enzyme families, target-sites of commonly used acaricides, and horizontally transferred genes implicated in plant-mite interactions and pesticide resistance. In a comparative analysis that incorporated sequences available for Panonychus citri, T. urticae, and insects, we identified radiations for detoxification gene families following the divergence of Panonychus and Tetranychus genera. Finally, we used the replicated RNA-seq data from the spirodiclofen susceptible and resistant strains to describe gene expression changes associated with resistance. A cytochrome P450 monooxygenase, as well as multiple carboxylcholinesterases, were differentially expressed between the susceptible and resistant strains, and provide a molecular entry point for understanding resistance to spirodiclofen, widely used to control P. ulmi populations.Conclusions: The new genomic resources and data that we present in this study for P. ulmi will substantially facilitate molecular studies of underlying mechanisms involved in acaricide resistance.

Epigenome modulation potentially provides a mechanism for organisms to adapt, within and between generations. However, neither the extent to which this occurs, nor the mechanisms involved are known. Here we investigate DNA methylation variation in Swedish Arabidopsis thaliana accessions grown at two different temperatures. Environmental effects were limited to transposons, where CHH methylation was found to increase with temperature. Genome-wide association studies (GWAS) revealed that the extensive CHH methylation variation was strongly associated with genetic variants in both cis and trans, including a major trans-association close to the DNA methyltransferase CMT2. Unlike CHH methylation, CpG gene body methylation (GBM) was not affected by growth temperature, but was instead correlated with the latitude of origin. Accessions from colder regions had higher levels of GBM for a significant fraction of the genome, and this was associated with increased transcription for the genes affected. GWAS revealed that this effect was largely due to trans-acting loci, many of which showed evidence of local adaptation.

The acaricides clofentezine, hexythiazox and etoxazole are commonly referred to as 'mite growth inhibitors', and clofentezine and hexythiazox have been used successfully for the integrated control of plant mite pests for decades. Although they are still important today, their mode of action has remained elusive. Recently, a mutation in chitin synthase 1 (CHS1) was linked to etoxazole resistance. In this study, we identified and investigated a Tetranychus urticae strain (HexR) harboring recessive, monogenic resistance to each of hexythiazox, clofentezine, and etoxazole. To elucidate if there is a common genetic basis for the observed cross-resistance, we adapted a previously developed bulk segregant analysis method to map with high resolution a single, shared resistance locus for all three compounds. This finding indicates that the underlying molecular basis for resistance to all three compounds is identical. This locus is centered on the CHS1 gene, and as supported by additional genetic and biochemical studies, a non-synonymous variant (I1017F) in CHS1 associates with resistance to each of the tested acaricides in HexR. Our findings thus demonstrate a shared molecular mode of action for the chemically diverse mite growth inhibitors clofentezine, hexythiazox and etoxazole as inhibitors of an essential, non-catalytic activity of CHS1. Given the previously documented cross-resistance between clofentezine, hexythiazox and the benzyolphenylurea (BPU) compounds flufenoxuron and cycloxuron, CHS1 should be also considered as a potential target-site of insecticidal BPUs.