Transcriptome profiling of a spirodiclofen susceptible and resistant strain of the European red mite Panonychus ulmi using strand-specific RNA-seq

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

High-resolution QTL mapping in Tetranychus urticae reveals acaricide-specific responses and common target-site resistance after selection by different METI-I acaricides

Arthropod herbivores cause dramatic crop losses, and frequent pesticide use has led to widespread resistance in numerous species. One such species, the two-spotted spider mite, Tetranychus urticae, is an extreme generalist herbivore and a major worldwide crop pest with a history of rapidly developing resistance to acaricides. Mitochondrial Electron Transport Inhibitors of complex I (METI-Is) have been used extensively in the last 25 years to control T. urticae around the globe, and widespread resistance to each has been documented. METI-I resistance mechanisms in T. urticae are likely complex, as increased metabolism by cytochrome P450 monooxygenases as well as a target-site mutation have been linked with resistance. To identify loci underlying resistance to the METI-I acaricides fenpyroximate, pyridaben and tebufenpyrad without prior hypotheses, we crossed a highly METI-I-resistant strain of T. urticae to a susceptible one, propagated many replicated populations over multiple generations with and without selection by each compound, and performed bulked segregant analysis genetic mapping. Our results showed that while the known H92R target-site mutation was associated with resistance to each compound, a genomic region that included cytochrome P450-reductase (CPR) was associated with resistance to pyridaben and tebufenpyrad. Within CPR, a single nonsynonymous variant distinguished the resistant strain from the sensitive one. Furthermore, a genomic region linked with tebufenpyrad resistance harbored a non-canonical member of the nuclear hormone receptor 96 (NHR96) gene family. This NHR96 gene does not encode a DNA-binding domain (DBD), an uncommon feature in arthropods, and belongs to an expanded family of 47 NHR96 proteins lacking DBDs in T. urticae. Our findings suggest that although cross-resistance to METI-Is involves known detoxification pathways, structural differences in METI-I acaricides have also resulted in resistance mechanisms that are compound-specific

Point mutations in the voltage-gated sodium channel gene associated with pyrethroid resistance in Iranian populations of the European red mite Panonychus ulmi

The European red mite Panonychus ulmi (Koch) is a major pest of apple trees worldwide and causes significant damage to apple orchards in Iran. Pyrethroid insecticides/acaricides, such as fenpropathrin and fenvalerate, are widely used to control P. ulmi, but their long-term use may lead to low efficacy. Earlier studies investigating pyrethroid resistance in closely related mites such as Tetranychus urticae revealed that pyrethroid resistance was associated with point mutations in the voltage-gated sodium channel gene (vgsc). The aim of this study was to investigate the biochemical and molecular mechanisms of fenpropathrin and fenvalerate resistance in Iranian populations of P. ulmi. Pyrethroid toxicity bioassays were carried out on different P. ulmi field populations. Marand (resistance ratio, RR = 149), Maracieh (RR = 90) and Mianeh2 (RR = 71) populations exhibited high levels of resistance to fenpropathrin, compared to a susceptible field population (Shahin Dej). Resistance was also observed for fenvalerate with resistance ratio's ranging from 2- to 20-fold. Synergism experiments and enzyme activity assays predicted a minor role for classical detoxification enzymes. In contrast, two amino acid substitutions in the VGSC, L1024V and F15381, that were previously shown to confer pyrethroid resistance, were detected in all three resistant P. ulmi populations and point towards target-site insensitivity as the most likely resistance mechanism. Furthermore, sequencing after cloning of vgsc fragments from single haploid males revealed the presence of multiple copies of vgsc in a highly resistant strain. The link between resistance mutations and vgsc copy number variation should be the subject of future study, as this might be used to develop molecular markers for monitoring pyrethroid resistance of P. ulmi in the field

The relative contribution of target-site mutations in complex acaricide resistant phenotypes as assessed by marker assisted backcrossing in Tetranychus urticae

The mechanisms underlying insecticide and acaricide resistance in insects and mites are often complex, including additive effects of target-site insensitivity, increased metabolism and transport. The extent to which target-site resistance mutations contribute to the resistance phenotype is, however, not well studied. Here, we used marker-assisted backcrossing to create 30 congenic lines carrying nine mutations (alone, or in combination in a few cases) associated with resistance to avermectins, pyrethroids, mite growth inhibitors and mitochondrial complex III inhibitors (QoI) in a polyphagous arthropod pest, the spider mite Tetranychus urticae. Toxicity tests revealed that mutations in the voltage-gated sodium channel, chitin synthase 1 and cytochrome b confer high levels of resistance and, when fixed in a population, these mutations alone can result in field failure of acaricide treatment. In contrast, although we confirmed the implication of mutations in glutamate-gated chloride channels in abamectin and milbemectin insensitivity, these mutations do not lead to the high resistance levels that are often reported in abamectin resistant strains of T. urticae. Overall, this study functionally validates reported target-site resistance mutations in T. urticae, by uncoupling them from additional mechanisms, allowing to finally investigate the strength of the conferred phenotype in vivo

Complex evolutionary dynamics of massively expanded chemosensory receptor families in an extreme generalist chelicerate herbivore

While mechanisms to detoxify plant produced, anti-herbivore compounds have been associated with plant host use by herbivores, less is known about the role of chemosensory perception in their life histories. This is especially true for generalists, including chelicerate herbivores that evolved herbivory independently from the more studied insect lineages. To shed light on chemosensory perception in a generalist herbivore, we characterized the chemosensory receptors (CRs) of the chelicerate two-spotted spider mite, Tetranychus urticae, an extreme generalist. Strikingly, T. urticae has more CRs than reported in any other arthropod to date. Including pseudogenes, 689 gustatory receptors were identified, as were 136 degenerin/Epithelial Na+ Channels (ENaCs) that have also been implicated as CRs in insects. The genomic distribution of T. urticae gustatory receptors indicates recurring bursts of lineage-specific proliferations, with the extent of receptor clusters reminiscent of those observed in the CR-rich genomes of vertebrates or C. elegans. Although pseudogenization of many gustatory receptors within clusters suggests relaxed selection, a subset of receptors is expressed. Consistent with functions as CRs, the genomic distribution and expression of ENaCs in lineage-specific T. urticae expansions mirrors that observed for gustatory receptors. The expansion of ENaCs in T. urticae to > 3-fold that reported in other animals was unexpected, raising the possibility that ENaCs in T. urticae have been co-opted to fulfill a major role performed by unrelated CRs in other animals. More broadly, our findings suggest an elaborate role for chemosensory perception in generalist herbivores that are of key ecological and agricultural importance

Long-term population studies uncover the genome structure and genetic basis of xenobiotic and host plant adaptation in the herbivore Tetranychus urticae

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

Molecular and genetic analysis of resistance to METI-I acaricides in Iranian populations of the citrus red mite Panonychus citri

The citrus red mite, Panonychus citri, is a major pest on citrus all around the world. Mitochondrial Electron Transport Inhibitors of complex I (METI-I) acaricides such as fenpyroximate have been used extensively to control P. citri populations, which resulted in multiple reports of METI-I resistant populations in the field. In this study, biochemical and molecular mechanisms of fenpyroximate resistance were investigated in P. citri. Seven populations were collected from Northern provinces of Iran. Resistance ratios were determined and reached up to 75-fold in comparison to a fenpyroximate susceptible population. Cross-resistance to two additional METI-I acaricides, pyridaben and tebufenpyrad, was detected. PBO synergism experiments, in vivo enzyme assays and gene expression analysis suggest a minor involvement of cytochrome P450 monooxygenases in fenpyroximate resistance, which is in contrast with many reported cases for the closely related Tetranychus urticae. Next, we determined the frequency of a well-known mutation in the target-site of METI-Is, the PSST subunit, associated with METI-I resistance. Indeed, the H92R substitution was detected in a highly fenpyroximate resistant P. citri population. Additionally, a new amino acid substitution at a conserved site in the PSST subunit was detected, A94V, with higher allele frequencies in a moderately resistant population. Marker-assisted back-crossing in a susceptible background confirmed the potential involvement of the newly discovered A94V mutation in fenpyroximate resistance. However, introduction of the A94V mutation in the PSST homologue of D. melanogaster using CRISPR-Cas9 did not result in fenpyroximate resistant flies. In addition, differences in binding curves between METI-Is and complex I measured directly, in isolated transgenic and wildtype mitochondria preparations, could not be found

Genetics and fitness costs of acaricide resistance in spider mites

The identification of the complex molecular mechanisms that govern the development of acaricide resistance in mite pests allows for more efficient pest management strategies. In chapter 2 we present the first genome-wide sequence data set for the spider mite Panonychus ulmi, and investigate the molecular mechanism underlying spirodiclofen resistance. In chapter 3, we shed more light on the contribution of target-site mutations to resistance to Mitochondrial Electron Transport Inhibitors I (METI-Is) in the spider mite pest Tetranychus urticae. We gather a body of experimental evidence that supports a causal link between the H92R mutation and METI-I resistance. In chapter 4, we introduce a marker-assisted backcrossing approach to create T. urticae lines that share a common pesticide-susceptible genomic background, excluding 9 loci that carry mutations that are associated with resistance to avermectins, pyrethroids, mite growth inhibitors, and mitochondrial complex III inhibitors. By this approach, this study functionally validates previously reported target-site resistance mutations in T. urticae. Subsequently, in chapter 5, we further take advantage of these near-isogenic lines to investigate possible fitness costs that are associated with the well-characterized acaricide target- site resistance mutations. Findings reported in this thesis are important considering the global spread of acaricide resistance, and within an Integrated Pest Management context

Selectivity and molecular stress responses to classical and botanical acaricides in the predatory mite Phytoseiulus persimilis Athias-Henriot (Acari: Phytoseiidae)

BACKGROUND Acaricide application remains an integral component of integrated pest management (IPM) for the two-spotted spider mite Tetranychus urticae. Species and strains of phytoseiid predatory mites vary significantly in their response to acaricides. For the success of IPM, it is imperative to identify the determinants of selectivity and molecular stress responses of acaricides in predatory mites. RESULTS The three classical acaricides bifenazate, cyflumetofen, and fenbutatin oxide did not affect the survival and fecundity of Phytoseiulus persimilis regardless of the route of exposure. Selectivity of the orange oil and terpenoid blend-based botanical acaricides was low via a combination of direct exposure, acaricide-laced diet, and residual exposure but improved when limiting exposure only to diet. To gain insights into the molecular stress responses, the transcriptome of P. persimilis was assembled. Subsequent gene expression analysis of predatory mites orally exposed to fenbutatin oxide and orange oil yielded only a limited xenobiotic stress response. In contrast, P. persimilis exhibited target-site resistance mutations, including I260M in SdhB, I1017M in CHS1, and kdr and super-kdr in VGSC. Extending the screen using available Phytoseiidae sequences uncovered I136T, S141F in cytb, G119S in AChE, and A2083V in ACC, well-known target-sites of acaricides. CONCLUSION Selectivity of the tested botanical acaricides to P. persimilis was low but could be enhanced by restricting exposure to a single route. Differential gene expression analysis did not show a robust induced stress response after sublethal exposure. In contrast, this study uncovered target-site mutations that may help to explain the physiological selectivity of several classical acaricides to phytoseiid predators