The desert woodrat (Neotoma lepida) induces a diversity of biotransformation genes in response to creosote bush resin.

Liver biotransformation enzymes have long been thought to enable animals to feed on diets rich in xenobiotic compounds. However, despite decades of pharmacological research in humans and rodents, little is known about hepatic gene expression in specialized mammalian herbivores feeding on toxic diets. Leveraging a recently identified population of the desert woodrat (Neotoma lepida) found to be highly tolerant to toxic creosote bush (Larrea tridentata), we explored the expression changes of suites of biotransformation genes in response to diets enriched with varying amounts of creosote resin. Analysis of hepatic RNA-seq data indicated a dose-dependent response to these compounds, including the upregulation of several genes encoding transcription factors and numerous phase I, II, and III biotransformation families. Notably, elevated expression of five biotransformation families - carboxylesterases, cytochromes P450, aldo-keto reductases, epoxide hydrolases, and UDP-glucuronosyltransferases - corresponded to species-specific duplication events in the genome, suggesting that these genes play a prominent role in N. lepida's adaptation to creosote bush. Building on pharmaceutical studies in model rodents, we propose a hypothesis for how the differentially expressed genes are involved in the biotransformation of creosote xenobiotics. Our results provide some of the first details about how these processes likely operate in the liver of a specialized mammalian herbivore.

Transcriptome Dynamics of Brassica juncea Leaves in Response to Omnivorous Beet Armyworm (Spodoptera exigua, Hübner).

Cruciferous plants manufacture glucosinolates (GSLs) as special and important defense compounds against insects. However, how insect feeding induces glucosinolates in Brassica to mediate insect resistance, and how plants regulate the strength of anti-insect defense response during insect feeding, remains unclear. Here, mustard (Brassica juncea), a widely cultivated Brassica plant, and beet armyworm (Spodoptera exigua), an economically important polyphagous pest of many crops, were used to analyze the changes in GSLs and transcriptome of Brassica during insect feeding, thereby revealing the plant-insect interaction in Brassica plants. The results showed that the content of GSLs began to significantly increase after 48 h of herbivory by S. exigua, with sinigrin as the main component. Transcriptome analysis showed that a total of 8940 DEGs were identified in mustard challenged with beet armyworm larvae. The functional enrichment results revealed that the pathways related to the biosynthesis of glucosinolate and jasmonic acid were significantly enriched by upregulated DEGs, suggesting that mustard might provide a defense against herbivory by inducing JA biosynthesis and then promoting GSL accumulation. Surprisingly, genes regulating JA catabolism and inactivation were also activated, and both JA signaling repressors (JAZs and JAMs) and activators (MYCs and NACs) were upregulated during herbivory. Taken together, our results indicate that the accumulation of GSLs regulated by JA signaling, and the regulation of active and inactive JA compound conversion, as well as the activation of JA signaling repressors and activators, collectively control the anti-insect defense response and avoid over-stunted growth in mustard during insect feeding.

A persistent major mutation in canonical jasmonate signaling is embedded in an herbivory-elicited gene network.

When insect herbivores attack plants, elicitors from oral secretions and regurgitants (OS) enter wounds during feeding, eliciting defense responses. These generally require plant jasmonate (JA) signaling, specifically, a jasmonoyl-L-isoleucine (JA-Ile) burst, for their activation and are well studied in the native tobacco Nicotiana attenuata. We used intraspecific diversity captured in a 26-parent MAGIC population planted in nature and an updated genome assembly to impute natural variation in the OS-elicited JA-Ile burst linked to a mutation in the JA-Ile biosynthetic gene NaJAR4. Experiments revealed that NaJAR4 variants were associated with higher fitness in the absence of herbivores but compromised foliar defenses, with two NaJAR homologues (4 and 6) complementing each other spatially and temporally. From decade-long seed collections of natural populations, we uncovered enzymatically inactive variants occurring at variable frequencies, consistent with a balancing selection regime maintaining variants. Integrative analyses of OS-induced transcriptomes and metabolomes of natural accessions revealed that NaJAR4 is embedded in a nonlinear complex gene coexpression network orchestrating responses to OS, which we tested by silencing four hub genes in two connected coexpressed networks and examining their OS-elicited metabolic responses. Lines silenced in two hub genes (NaGLR and NaFB67) co-occurring in the NaJAR4/6 module showed responses proportional to JA-Ile accumulations; two from an adjacent module (NaERF and NaFB61) had constitutively expressed defenses with high resistance. We infer that mutations with large fitness consequences can persist in natural populations due to compensatory responses from gene networks, which allow for diversification in conserved signaling pathways and are generally consistent with predictions of an omnigene model.

Evolution of chemosensory and detoxification gene families across herbivorous Drosophilidae.

Herbivorous insects are exceptionally diverse, accounting for a quarter of all known eukaryotic species, but the genomic basis of adaptations that enabled this dietary transition remains poorly understood. Many studies have suggested that expansions and contractions of chemosensory and detoxification gene families-genes directly mediating interactions with plant chemical defenses-underlie successful plant colonization. However, this hypothesis has been challenging to test because the origins of herbivory in many insect lineages are ancient (>150 million years ago (mya)), obscuring genomic evolutionary patterns. Here, we characterized chemosensory and detoxification gene family evolution across Scaptomyza, a genus nested within Drosophila that includes a recently derived (<15 mya) herbivore lineage of mustard (Brassicales) specialists and carnation (Caryophyllaceae) specialists, and several nonherbivorous species. Comparative genomic analyses revealed that herbivorous Scaptomyza has among the smallest chemosensory and detoxification gene repertoires across 12 drosophilid species surveyed. Rates of gene turnover averaged across the herbivore clade were significantly higher than background rates in over half of the surveyed gene families. However, gene turnover was more limited along the ancestral herbivore branch, with only gustatory receptors and odorant-binding proteins experiencing strong losses. The genes most significantly impacted by gene loss, duplication, or changes in selective constraint were those involved in detecting compounds associated with feeding on living plants (bitter or electrophilic phytotoxins) or their ancestral diet (fermenting plant volatiles). These results provide insight into the molecular and evolutionary mechanisms of plant-feeding adaptations and highlight gene candidates that have also been linked to other dietary transitions in Drosophila.

Evolution and genomic basis of the plant-penetrating ovipositor: a key morphological trait in herbivorous Drosophilidae.

Herbivorous insects are extraordinarily diverse, yet are found in only one-third of insect orders. This skew may result from barriers to plant colonization, coupled with phylogenetic constraint on plant-colonizing adaptations. The plant-penetrating ovipositor, however, is one trait that surmounts host plant physical defences and may be evolutionarily labile. Ovipositors densely lined with hard bristles have evolved repeatedly in herbivorous lineages, including within the Drosophilidae. However, the evolution and genetic basis of this innovation has not been well studied. Here, we focused on the evolution of this trait in Scaptomyza, a genus sister to Hawaiian Drosophila, that contains a herbivorous clade. Our phylogenetic approach revealed that ovipositor bristle number increased as herbivory evolved in the Scaptomyza lineage. Through a genome-wide association study, we then dissected the genomic architecture of variation in ovipositor bristle number within S. flava. Top-associated variants were enriched for transcriptional repressors, and the strongest associations included genes contributing to peripheral nervous system development. Individual genotyping supported the association at a variant upstream of Gαi, a neural development gene, contributing to a gain of 0.58 bristles/major allele. These results suggest that regulatory variation involving conserved developmental genes contributes to this key morphological trait involved in plant colonization.

Common resistance mechanisms are deployed by plants against sap-feeding herbivorous insects: insights from a meta-analysis and systematic review.

Despite their abundance and economic importance, the mechanism of plant resistance to sap-feeding insects remains poorly understood. Here we deploy meta-analysis and data synthesis methods to evaluate the results from electrophysiological studies describing feeding behaviour experiments where resistance mechanisms were identified, focussing on studies describing host-plant resistance and non-host resistance mechanisms. Data were extracted from 108 studies, comprising 41 insect species across eight insect taxa and 12 host-plant families representing over 30 species. Results demonstrate that mechanisms deployed by resistant plants have common consequences on the feeding behaviour of diverse insect groups. We show that insects feeding on resistant plants take longer to establish a feeding site and have their feeding duration suppressed two-fold compared with insects feeding on susceptible plants. Our results reveal that traits contributing towards resistant phenotypes are conserved across plant families, deployed against taxonomically diverse insect groups, and that the underlying resistance mechanisms are conserved. These findings provide a new insight into plant-insect interaction and highlight the need for further mechanistic studies across diverse taxa.

Additive genetic effects in interacting species jointly determine the outcome of caterpillar herbivory.

Plant-insect interactions are common and important in basic and applied biology. Trait and genetic variation can affect the outcome and evolution of these interactions, but the relative contributions of plant and insect genetic variation and how these interact remain unclear and are rarely subject to assessment in the same experimental context. Here, we address this knowledge gap using a recent host-range expansion onto alfalfa by the Melissa blue butterfly. Common garden rearing experiments and genomic data show that caterpillar performance depends on plant and insect genetic variation, with insect genetics contributing to performance earlier in development and plant genetics later. Our models of performance based on caterpillar genetics retained predictive power when applied to a second common garden. Much of the plant genetic effect could be explained by heritable variation in plant phytochemicals, especially saponins, peptides, and phosphatidyl cholines, providing a possible mechanistic understanding of variation in the species interaction. We find evidence of polygenic, mostly additive effects within and between species, with consistent effects of plant genotype on growth and development across multiple butterfly species. Our results inform theories of plant-insect coevolution and the evolution of diet breadth in herbivorous insects and other host-specific parasites.

Whole genome assembly of the armored loricariid catfish Ancistrus triradiatus highlights herbivory signatures.

The catfish Ancistrus triradiatus belongs to the species-rich family Loricariidae. Loricariids display remarkable traits such as herbivory, a benthic lifestyle, the absence of scales but the presence of dermal bony plates. They are exported as ornamental fish worldwide, with escaped fishes becoming a threat locally. Although genetic and phylogenetic studies are continuously increasing and developmental genetic investigations are underway, no genome assembly has been formally proposed for Loricariidae yet. We report a high-quality genome assembly of Ancistrus triradiatus using long and short reads, and a newly assembled transcriptome. The genome assembly is composed of 9530 scaffolds, including 85.6% of ray-finned fish BUSCOs, and 26,885 predicted protein-coding genes. The genomic GC content is higher than in other catfishes, reflecting the higher metabolism associated with herbivory. The examination of the SCPP gene family indicates that the genes presumably triggering scale loss when absent, are present in the scaleless A. triradiatus, questioning their explanatory role. The analysis of the opsin gene repertoire revealed that gene losses associated to the nocturnal lifestyle of catfishes were not entirely found in A. triradiatus, as the UV-sensitive opsin 5 is present. Finally, most gene family expansions were related to immunity except the gamma crystallin gene family which controls pupil shape and sub-aquatic vision. Thus, the genome of A. triradiatus reveals that fish herbivory may be related to the photic zone habitat, conditions metabolism, photoreception and visual functions. This genome is the first for the catfish suborder Loricarioidei and will serve as backbone for future genetic, developmental and conservation studies.

Diet evolution of carnivorous and herbivorous mammals in Laurasiatheria.

BACKGROUND: Laurasiatheria contains taxa with diverse diets, while the molecular basis and evolutionary history underlying their dietary diversification are less clear. RESULTS: In this study, we used the recently developed molecular phyloecological approach to examine the adaptive evolution of digestive system-related genes across both carnivorous and herbivorous mammals within Laurasiatheria. Our results show an intensified selection of fat and/or protein utilization across all examined carnivorous lineages, which is consistent with their high-protein and high-fat diets. Intriguingly, for herbivorous lineages (ungulates), which have a high-carbohydrate diet, they show a similar selection pattern as that of carnivorous lineages. Our results suggest that for the ungulates, which have a specialized digestive system, the selection intensity of their digestive system-related genes does not necessarily reflect loads of the nutrient components in their diets but appears to be positively related to the loads of the nutrient components that are capable of being directly utilized by the herbivores themselves. Based on these findings, we reconstructed the dietary evolution within Laurasiatheria, and our results reveal the dominant carnivory during the early diversification of Laurasiatheria. In particular, our results suggest that the ancestral bats and the common ancestor of ruminants and cetaceans may be carnivorous as well. We also found evidence of the convergent evolution of one fat utilization-related gene, APOB, across carnivorous taxa. CONCLUSIONS: Our molecular phyloecological results suggest that digestive system-related genes can be used to determine the molecular basis of diet differentiations and to reconstruct ancestral diets.

Successes and limitations of quantitative diet metabarcoding in a small, herbivorous mammal.

DNA metabarcoding is widely used to determine wild animal diets, but whether this technique provides accurate, quantitative measurements is still under debate. To test our ability to accurately estimate the abundance of dietary items using metabarcoding, we fed wild-caught desert woodrats (Neotoma lepida) diets consisting of constant amounts of juniper (Juniperus osteosperma, 15%) and varying amounts of creosote (Larrea tridentata, 1%-60%), cactus (Opuntia sp., 0%-100%) and commercial chow (0%-85%). Using metabarcoding, we compared the representation of items in the original diet samples to that in the faecal samples to test the sensitivity and accuracy of diet metabarcoding, the performance of different bioinformatic pipelines and our ability to correct sequence counts. Metabarcoding, using standard trnL primers, detected creosote, juniper and chow. Different pipelines for assigning taxonomy performed similarly. While creosote was detectable at dietary proportions as low as 1%, we failed to detect cactus in most samples, probably due to a primer mismatch. Creosote read counts increased as its proportion in the diet increased, and we could differentiate when creosote was a minor and major component of the diet. However, we found that estimates of juniper and creosote varied. Using previously suggested methods to correct these errors did not improve accuracy estimates of creosote, but did reduce error for juniper and chow. Our results indicate that metabarcoding can provide quantitative information on dietary composition, but may be limited. We suggest that researchers use caution when quantitatively interpreting diet metabarcoding results unless they first experimentally determine the extent of possible biases.

Evolution of Olfactory Receptors Tuned to Mustard Oils in Herbivorous Drosophilidae.

The diversity of herbivorous insects is attributed to their propensity to specialize on toxic plants. In an evolutionary twist, toxins betray the identity of their bearers when herbivores coopt them as cues for host-plant finding, but the evolutionary mechanisms underlying this phenomenon are poorly understood. We focused on Scaptomyza flava, an herbivorous drosophilid specialized on isothiocyanate (ITC)-producing (Brassicales) plants, and identified Or67b paralogs that were triplicated as mustard-specific herbivory evolved. Using in vivo heterologous systems for the expression of olfactory receptors, we found that S. flava Or67bs, but not the homologs from microbe-feeding relatives, responded selectively to ITCs, each paralog detecting different ITC subsets. Consistent with this, S. flava was attracted to ITCs, as was Drosophila melanogaster expressing S. flava Or67b3 in the homologous Or67b olfactory circuit. ITCs were likely coopted as olfactory attractants through gene duplication and functional specialization (neofunctionalization and subfunctionalization) in S. flava, a recently derived herbivore.

Evolutionary and dietary relationships of wild mammals based on the gut microbiome.

Gut microbiome influence the health and evolution of mammals and multiple factors modulate the structure and function of gut microbiome. However, the specific changes of the diets and phylogeny on the gut microbiome were unclear. Here, we compared the gut microbiome of 16 rare wild mammals. All data (>200G 16S rRNA gene sequences) were generated using a high-throughput sequencing platform. Firmicutes and Bacteroidetes were the most predominant phyla in all mammals. However, Proteobacteria was an additionally dominant phylum specifically detected in the microbiome of carnivores and omnivores. Moreover, the dominant phyla in canids were Firmicutes, Bacteroidetes, Proteobacteria, and Fusobacteria. Phylogenetic reconstructions based on the gut microbiome and mitochondrial genome of these mammals were similar. The impact of the host on the microbiome community composition was most evident when considering conspecific and congeneric relationships. Similarity clustering showed that the gut microbiome of herbivores was clustered together, and the other clade comprised both omnivores and carnivores. Collectively, these results revealed that phylogenetic relationships and diet have an important impact on the gut microbiome, and thus the gut microbiome community composition may reflect both the phylogenetic relationships and diets. This study provides valuable basic data to facilitate future efforts related to animal conservation and health.

Rapid PCR-based method for herbivore dietary evaluation using plant-specific primers.

Polyphagous pests cause significant economic loss worldwide through feeding damage on various cash crops. However, their diets in agricultural landscapes remain largely unexplored. Pest dietary evaluation in agricultural fields is a challenging task currently approached through visual observation of plant feeding and microscopic identification of semi-digested plant material in pest's guts. While molecular gut content analysis using metabarcoding approaches using universal primers (e.g., rbcl and trnL) have been successful in evaluating polyphagous pest diet, this method is relatively costly and time-consuming. Hence, there is a need for a rapid, specific, sensitive, and cost-effective method to screen for crops in the gut of pests. This is the first study to develop plant-specific primers that target various regions of their genomes, designed using a whole plant genome sequence. We selected Verticillium wilt disease resistance protein (VE-1) and pathogenesis related protein-coding genes 1-5 (PR-1-5) as our targets and designed species-specific primers for 14 important crops in the agroecosystems. Using amplicon sizes ranging from 115 to 407 bp, we developed two multiplex primer mixes that can separate nine and five plant species per PCR reaction, respectively. These two designed primer mixes provide a rapid, sensitive and specific route for polyphagous pest dietary evaluation in agroecosystems. This work will enable future research to rapidly expand our knowledge on the diet preference and range of crops that pests consume in various agroecosystems, which will help in the redesign and development of new crop rotation regimes to minimize polyphagous pest pressure and damage on crops.

Behavioral and Physiological Plasticity Provides Insights into Molecular Based Adaptation Mechanism to Strain Shift in Spodoptera frugiperda.

How herbivorous insects adapt to host plants is a key question in ecological and evolutionary biology. The fall armyworm, (FAW) Spodoptera frugiperda (J.E. Smith), although polyphagous and a major pest on various crops, has been reported to have a rice and corn (maize) feeding strain in its native range in the Americas. The species is highly invasive and has recently established in China. We compared behavioral changes in larvae and adults of a corn population (Corn) when selected on rice (Rice) and the molecular basis of these adaptational changes in midgut and antennae based on a comparative transcriptome analysis. Larvae of S. frugiperda reared on rice plants continuously for 20 generations exhibited strong feeding preference for with higher larval performance and pupal weight on rice than on maize plants. Similarly, females from the rice selected population laid significantly more eggs on rice as compared to females from maize population. The most highly expressed DEGs were shown in the midgut of Rice vs. Corn. A total of 6430 DEGs were identified between the populations mostly in genes related to digestion and detoxification. These results suggest that potential adaptations for feeding on rice crops, may contribute to the current rapid spread of fall armyworm on rice crops in China and potentially elsewhere. Consistently, highly expressed DEGs were also shown in antennae; a total of 5125 differentially expressed genes (DEGs) s were identified related to the expansions of major chemosensory genes family in Rice compared to the Corn feeding population. These results not only provide valuable insight into the molecular mechanisms in host plants adaptation of S. frugiperda but may provide new gene targets for the management of this pest.

Late Quaternary dynamics of Arctic biota from ancient environmental genomics.

During the last glacial-interglacial cycle, Arctic biotas experienced substantial climatic changes, yet the nature, extent and rate of their responses are not fully understood1-8. Here we report a large-scale environmental DNA metagenomic study of ancient plant and mammal communities, analysing 535 permafrost and lake sediment samples from across the Arctic spanning the past 50,000 years. Furthermore, we present 1,541 contemporary plant genome assemblies that were generated as reference sequences. Our study provides several insights into the long-term dynamics of the Arctic biota at the circumpolar and regional scales. Our key findings include: (1) a relatively homogeneous steppe-tundra flora dominated the Arctic during the Last Glacial Maximum, followed by regional divergence of vegetation during the Holocene epoch; (2) certain grazing animals consistently co-occurred in space and time; (3) humans appear to have been a minor factor in driving animal distributions; (4) higher effective precipitation, as well as an increase in the proportion of wetland plants, show negative effects on animal diversity; (5) the persistence of the steppe-tundra vegetation in northern Siberia enabled the late survival of several now-extinct megafauna species, including the woolly mammoth until 3.9 ± 0.2 thousand years ago (ka) and the woolly rhinoceros until 9.8 ± 0.2 ka; and (6) phylogenetic analysis of mammoth environmental DNA reveals a previously unsampled mitochondrial lineage. Our findings highlight the power of ancient environmental metagenomics analyses to advance understanding of population histories and long-term ecological dynamics.

Reference Gene Selection for Analyzing the Transcription Patterns of Two Fatty Acyl-CoA Reductase Genes From Paracoccus marginatus (Hemiptera: Pseudococcidae).

Paracoccus marginatus (Hemiptera: Pseudococcidae), known as the papaya mealybug, could cause considerable yield loss of several plants. To date, there is no molecular-based study of P. marginatus. Fatty acyl-CoA reductases (FARs) are key enzymes involved in wax synthesis. In the present study, we cloned and characterized coding sequences (CDS) of two FAR genes from P. marginatus. The results showed that PmFAR1 and PmFAR2 CDS were 1,590 and 1,497 bp in length, respectively, and sequence analysis indicated that these two genes both had the conservative motifs belonging to FAR_C superfamily. Furthermore, seven candidate reference genes were analyzed for their expression stability by using common algorithms including comparative ΔCq method, geNorm, NormFinder, BestKeeper, and RefFinder. Eventually, ß-actin and GAPDH were the best reference genes in evaluating the expression of those two FAR genes. We found that PmFAR1 and PmFAR2 showed distinct expression patterns in different life stages. Moreover, the transcription of PmFAR1 and PmFAR2 in P. marginatus fed on resistant cassava cultivars was significantly lower compared with those fed on susceptible ones, indicating the potential function of FAR genes in cassava resistance to P. marginatus. The present study might help in better understanding the molecular mechanism of cassava resistance to mealybug.

Large-scale identification of extracellular plant miRNAs in mammals implicates their dietary intake.

Extracellular microRNAs (miRNAs) have been proposed to function in cross-kingdom gene regulation. Among these, plant-derived miRNAs of dietary origin have been reported to survive the harsh conditions of the human digestive system, enter the circulatory system, and regulate gene expression and metabolic function. However, definitive evidence supporting the presence of plant-derived miRNAs of dietary origin in mammals has been difficult to obtain due to limited sample sizes. We have developed a bioinformatics pipeline (ePmiRNA_finder) that provides strident miRNA classification and applied it to analyze 421 small RNA sequencing data sets from 10 types of human body fluids and tissues and comparative samples from carnivores and herbivores. A total of 35 miRNAs were identified that map to plants typically found in the human diet and these miRNAs were found in at least one human blood sample and their abundance was significantly different when compared to samples from human microbiome or cow. The plant-derived miRNA profiles were body fluid/tissue-specific and highly abundant in the brain and the breast milk samples, indicating selective absorption and/or the ability to be transported across tissue/organ barriers. Our data provide conclusive evidence for the presence of plant-derived miRNAs as a consequence of dietary intake and their cross-kingdom regulatory function within human circulating system.

The draft genome of the Asian corn borer yields insights into ecological adaptation of a devastating maize pest.

The Asian corn borer (ACB) is the most devastating pest on maize in the western Pacific region of Asia. Despite broad interests in insecticide resistance, seasonal adaptation, and larval color mimicry regarding the ACB system, lacking of reference genomic information and a powerful gene editing approach have hindered the in-depth studies of these aspects. Here we present a 455.7 Mb draft genome of ACB with 98.4% completeness. Comparative genomics analysis showed an evident expansion in gene families of gustatory receptors (105), which is related to polyphagous characteristics. Based on the comparative transcriptome analysis of resistant and susceptible ACB against Bt Cry1Ab toxin, we identified 26 genes related to Cry1Ab resistance. Additionally, transcriptomics of insects exposed to conditions of low temperature and diapause (LT) vs. room temperature and diapause (RT) provided insights into the genetic mechanisms of cold adaptation. We also successfully developed an efficient CRISPR/Cas9-based genome editing system and applied it to explore the role of color pattern genes in the ecological adaptation of ACB. Taken together, our study provides a fully annotated high-quality reference genome and efficient gene editing system to realize the potential of ACB as a study system to address important biological questions such as insecticide resistance, seasonal adaptation, and coloration. These valuable genomic resources will also benefit the development of novel strategies for maize pest management.

Adaptive divergence and post-zygotic barriers to gene flow between sympatric populations of a herbivorous mite.

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.

Comparative analysis of mite genomes reveals positive selection for diet adaptation.

Diet is a powerful evolutionary force for species adaptation and diversification. Acari is one of the most abundant clades of Arachnida, exhibiting diverse dietary types, while the underlying genetic adaptive mechanisms are not fully understood. Based on comparative analyses of 15 Acari genomes, we found genetic bases for three specialized diets. Herbivores experienced stronger selection pressure than other groups; the olfactory genes and gene families involving metabolizing toxins showed strong adaptive signals. Genes and gene families related to anticoagulation, detoxification, and haemoglobin digestion were found to be under strong selection pressure or significantly expanded in the blood-feeding species. Lipid metabolism genes have a faster evolutionary rate and been subjected to greater selection pressures in fat-feeding species; one positively selected site in the fatty-acid amide hydrolases 2 gene was identified. Our research provides a new perspective for the evolution of Acari and offers potential target loci for novel pesticide development.