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.

Hybridization in the absence of an ecotone favors hybrid success in woodrats (Neotoma spp.).

Hybridization is a common process that has broadly impacted the evolution of multicellular eukaryotes; however, how ecological factors influence this process remains poorly understood. Here, we report the findings of a 3-year recapture study of the Bryant's woodrat (Neotoma bryanti) and desert woodrat (Neotoma lepida), two species that hybridize within a creosote bush (Larrea tridentata) shrubland in Whitewater, CA, USA. We used a genotype-by-sequencing approach to characterize the ancestry distribution of individuals across this hybrid zone coupled with Cormack-Jolly-Seber modeling to describe demography. We identified a high frequency of hybridization at this site with ~40% of individuals possessing admixed ancestry, which is the result of multigenerational backcrossing and advanced hybrid-hybrid crossing. F1, F2, and advanced generation hybrids had apparent survival rates similar to parental N. bryanti, while parental and backcross N. lepida had lower apparent survival rates and were far less abundant. Compared to bimodal hybrid zones where hybrids are often rare and selected against, we find that hybrids at Whitewater are common and have comparable survival to the dominant parental species, N. bryanti. The frequency of hybridization at Whitewater is therefore likely limited by the abundance of the less common parental species, N. lepida, rather than selection against hybrids.

Trio-binned genomes of the woodrats Neotoma bryanti and Neotoma lepida reveal novel gene islands and rapid copy number evolution of xenobiotic metabolizing genes.

The genomic architecture underlying the origins and maintenance of biodiversity is an increasingly accessible feature of species, due in large part to third-generation sequencing and novel analytical toolsets. Applying these techniques to woodrats (Neotoma spp.) provides a unique opportunity to study how herbivores respond to environmental change. Neotoma bryanti and N. lepida independently achieved a major dietary feat in the aftermath of a natural climate change event: switching to the novel, toxic food source creosote bush (Larrea tridentata). To better understand the genetic mechanisms underlying this ability, we employed a trio binning sequencing approach with a N. bryanti × N. lepida F1 hybrid, allowing the simultaneous assembly of genomes representing each parental species. The resulting phased, chromosome-level, highly complete haploid references enabled us to explore the genomic architecture of several gene families-cytochromes P450, UDP-glucuronosyltransferases (UGTs), and ATP-binding cassette (ABC) transporters-known to play key roles in the metabolism of naturally occurring toxic dietary compounds. In addition to duplication events in the ABCG and UGT2B subfamilies, we found expansions in three P450 gene families (2A, 2B, 3A), including the evolution of multiple novel gene islands within the 2B and 3A subfamilies, which may have provided the crucial substrate for dietary adaptation. Our assemblies demonstrate that trio binning from an F1 hybrid rodent effectively recovers parental genomes from species that diverged more than a million years ago.

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.

Wild herbivorous mammals (genus Neotoma) host a diverse but transient assemblage of fungi.

Fungi are often overlooked in microbiome research and, as a result, little is known about the mammalian mycobiome. Although frequently detected in vertebrate guts and known to contribute to digestion in some herbivores, whether these eukaryotes are a persistent part of the mammalian gut microbiome remains contentious. To address this question, we sampled fungi from wild woodrats (Neotoma spp.) collected from 25 populations across the southwestern United States. For each animal, we collected a fecal sample in the wild, and then re-sampled the same individual after a month in captivity on a controlled diet. We characterized and quantified fungi using three techniques: ITS metabarcoding, shotgun metagenomics and qPCR. Wild individuals contained diverse fungal assemblages dominated by plant pathogens, widespread molds, and coprophilous taxa primarily in Ascomycota and Mucoromycota. Fungal abundance, diversity and composition differed between individuals, and was primarily influenced by animal geographic origin. Fungal abundance and diversity significantly declined in captivity, indicating that most fungi in wild hosts came from diet and environmental exposure. While this suggests that these mammals lack a persistent gut mycobiome, natural fungal exposure may still impact fungal dispersal and animal health.

Toxin tolerance across landscapes: Ecological exposure not a prerequisite.

Little is known about the tolerances of mammalian herbivores to plant specialized metabolites across landscapes.We investigated the tolerances of two species of herbivorous woodrats, Neotoma lepida (desert woodrat) and Neotoma bryanti (Bryant's woodrat) to creosote bush (Larrea tridentata), a widely distributed shrub with a highly toxic resin. Woodrats were sampled from 13 locations both with and without creosote bush across a 900 km transect in the US southwest. We tested whether these woodrat populations consume creosote bush using plant metabarcoding of feces and quantified their tolerance to creosote bush through feeding trials using chow amended with creosote resin.Toxin tolerance was analyzed in the context of population structure across collection sites with microsatellite analyses. Genetic differentiation among woodrats collected from different locations was minimal within either species. Tolerance differed substantially between the two species, with N. lepida persisting 20% longer than N. bryanti in feeding trials with creosote resin. Furthermore, in both species, tolerance to creosote resin was similar among woodrats near or within creosote bush habitat. In both species, woodrats collected greater than 25 km from creosote had markedly lower tolerances to creosote resin compared to animals from within the range of creosote bush.The results imply that mammalian herbivores are adapted to the specialized metabolites of plants in their diet, and that this tolerance can extend several kilometers outside of the range of dietary items. That is, direct ecological exposure to the specialized chemistry of particular plant species is not a prerequisite for tolerance to these compounds. These findings lay the groundwork for additional studies to investigate the genetic mechanisms underlying toxin tolerance and to identify how these mechanisms are maintained across landscape-level scales in mammalian herbivores.

Addressing nontarget amplification in DNA metabarcoding studies of arthropod-feeding rodents.

High-throughput sequencing approaches have revolutionized how we study animal diets by enabling the detection of dietary components from the metabarcoding of DNA in excrement. Mitochondrial cytochrome oxidase C subunit I (mtCOI) DNA metabarcoding is commonly used to study the diets of arthropod-feeding animals; however, this approach is susceptible to nontarget amplification of the consumer species mtCOI locus. Nontarget amplification is often an unforeseen complication that can drastically reduce the quality and utility of the results generated by high-throughput amplicon sequencing. By interrogating the diets of new world rodents in the genus Neotoma (woodrats) in both natural and captive settings, we demonstrate that nontarget amplification can drastically reduce the total read abundance of detected arthropod taxa in fecal samples and inhibit downstream analyses of dietary diversity and composition metrics. Using the results from these investigations, we offer a guide on how to identify concerns for nontarget amplification when selecting degenerate primers for DNA metabarcoding studies and recommend several approaches that can reduce or eliminate nontarget amplification. Lastly, for the community interested in investigating the diets of arthropod-feeding rodents, we generated a database containing the degree of mismatch between publicly available Rodentia mtCOI sequences and four common universal mtCOI primer sets to be used as a resource for inferring the relative risk of nontarget amplification when designing arthropod metabarcoding studies in rodent systems. This guide will be especially useful for researchers working with consumer species that have not previously been studied.

Microbiome stability and structure is governed by host phylogeny over diet and geography in woodrats (Neotoma spp.).

The microbiome is critical for host survival and fitness, but gaps remain in our understanding of how this symbiotic community is structured. Despite evidence that related hosts often harbor similar bacterial communities, it is unclear whether this pattern is due to genetic similarities between hosts or to common ecological selection pressures. Here, using herbivorous rodents in the genus Neotoma, we quantify how geography, diet, and host genetics, alongside neutral processes, influence microbiome structure and stability under natural and captive conditions. Using bacterial and plant metabarcoding, we first characterized dietary and microbiome compositions for animals from 25 populations, representing seven species from 19 sites across the southwestern United States. We then brought wild animals into captivity, reducing the influence of environmental variation. In nature, geography, diet, and phylogeny collectively explained ∼50% of observed microbiome variation. Diet and microbiome diversity were correlated, with different toxin-enriched diets selecting for distinct microbial symbionts. Although diet and geography influenced natural microbiome structure, the effects of host phylogeny were stronger for both wild and captive animals. In captivity, gut microbiomes were altered; however, responses were species specific, indicating again that host genetic background is the most significant predictor of microbiome composition and stability. In captivity, diet effects declined and the effects of host genetic similarity increased. By bridging a critical divide between studies in wild and captive animals, this work underscores the extent to which genetics shape microbiome structure and stability in closely related hosts.

Maize Inbred Line B96 Is the Source of Large-Effect Loci for Resistance to Generalist but Not Specialist Spider Mites.

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.

Genome streamlining in a minute herbivore that manipulates its host plant.

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.

Arthropods are a group of invertebrates that include insects ­ such as flies or beetles ­ arachnids ­ like spiders or scorpions ­ and crustaceans ­ including shrimp and woodlice. One of the tiniest species of arthropods, measuring less than 0.2 millimeters, is the tomato russet mite Aculops lycopersici. This arachnid is among the smallest animals on Earth, even smaller than some single-celled organisms, and only has four legs, unlike other arachnids. It is a major pest on tomato plants, which are toxic to many other animals, and it feeds on the top cell layer of the stems and leaves. Tomato growers need a way to identify and treat tomato russet mite infestations, but this tiny species remains something of a mystery. One way to tackle this pest may be to take a closer look at its genome, as this could reveal what genes the mite uses to detoxify its diet. Examining the mite's genome could also reveal information about how evolution handles creatures becoming smaller. An area of particular interest is the overall size of its genome. Not all of the DNA in a genome is part of genes that code for proteins; there are also sections of so-called 'non-coding' DNA. These sequences play important roles in controlling how and when cells use their genes. In the human genome, for example, just 1% of the DNA codes for protein. In fact, most human protein-coding genes are interrupted by sequences of non-coding DNA, called introns. Here, Greenhalgh, Dermauw et al. sequence the entire tomato russet mite genome and reveal that not only is the mite's body size miniature: these tiny animals have the smallest arthropod genome reported to date, almost a hundred times smaller than the human genome. Part of this genetic miniaturization seems to be down to massive loss of non-coding DNA. Around 40% of the mite genome codes for protein, and 80% of its protein coding genes contain no introns. The rest of the miniaturization involves loss of genes themselves. The mites have lost some of the genes that determine body structure, which could explain why they have fewer legs than other arachnids. Additionally, they only carry a small set of genes involved in sensing chemicals and clearing toxins, which could explain why they are mostly found on tomato plants. Greenhalgh, Dermauw et al.'s findings shed light on what may happen to the genome at the extremes of size evolution. Sequencing the genomes of other mites could reveal when in evolutionary history this genetic miniaturization occurred. Furthermore, a better understanding of the tomato russet mite genome could lead to the development of methods to detect the infestation of plants earlier and be highly beneficial for tomato agriculture.

Convergent evolution of cytochrome P450s underlies independent origins of keto-carotenoid pigmentation in animals.

Keto-carotenoids contribute to many important traits in animals, including vision and coloration. In a great number of animal species, keto-carotenoids are endogenously produced from carotenoids by carotenoid ketolases. Despite the ubiquity and functional importance of keto-carotenoids in animals, the underlying genetic architectures of their production have remained enigmatic. The body and eye colorations of spider mites (Arthropoda: Chelicerata) are determined by ß-carotene and keto-carotenoid derivatives. Here, we focus on a carotenoid pigment mutant of the spider mite Tetranychus kanzawai that, as shown by chromatography, lost the ability to produce keto-carotenoids. We employed bulked segregant analysis and linked the causal locus to a single narrow genomic interval. The causal mutation was fine-mapped to a minimal candidate region that held only one complete gene, the cytochrome P450 monooxygenase CYP384A1, of the CYP3 clan. Using a number of genomic approaches, we revealed that an inactivating deletion in the fourth exon of CYP384A1 caused the aberrant pigmentation. Phylogenetic analysis indicated that CYP384A1 is orthologous across mite species of the ancient Trombidiformes order where carotenoids typify eye and body coloration, suggesting a deeply conserved function of CYP384A1 as a carotenoid ketolase. Previously, CYP2J19, a cytochrome P450 of the CYP2 clan, has been identified as a carotenoid ketolase in birds and turtles. Our study shows that selection for endogenous production of keto-carotenoids led to convergent evolution, whereby cytochrome P450s were independently co-opted in vertebrate and invertebrate animal lineages.

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.

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.

Generalist and Specialist Mite Herbivores Induce Similar Defense Responses in Maize and Barley but Differ in Susceptibility to Benzoxazinoids.

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.

The effect of insecticide synergist treatment on genome-wide gene expression in a polyphagous pest.

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.

Disruption of a horizontally transferred phytoene desaturase abolishes carotenoid accumulation and diapause in Tetranychus urticae.

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.

Epistatic and allelic interactions control expression of ribosomal RNA gene clusters in Arabidopsis thaliana.

BACKGROUND: Ribosomal RNA (rRNA) accounts for the majority of the RNA in eukaryotic cells, and is encoded by hundreds to thousands of nearly identical gene copies, only a subset of which are active at any given time. In Arabidopsis thaliana, 45S rRNA genes are found in two large ribosomal DNA (rDNA) clusters and little is known about the contribution of each to the overall transcription pattern in the species. RESULTS: By taking advantage of genome sequencing data from the 1001 Genomes Consortium, we characterize rRNA gene sequence variation within and among accessions. Notably, variation is not restricted to the pre-rRNA sequences removed during processing, but it is also present within the highly conserved ribosomal subunits. Through linkage mapping we assign these variants to a particular rDNA cluster unambiguously and use them as reporters of rDNA cluster-specific expression. We demonstrate that rDNA cluster-usage varies greatly among accessions and that rDNA cluster-specific expression and silencing is controlled via genetic interactions between entire rDNA cluster haplotypes (alleles). CONCLUSIONS: We show that rRNA gene cluster expression is controlled via complex epistatic and allelic interactions between rDNA haplotypes that apparently regulate the entire rRNA gene cluster. Furthermore, the sequence polymorphism we discovered implies that the pool of rRNA in a cell may be heterogeneous, which could have functional consequences.

Genomic Rearrangements in Arabidopsis Considered as Quantitative Traits.

To understand the population genetics of structural variants and their effects on phenotypes, we developed an approach to mapping structural variants that segregate in a population sequenced at low coverage. We avoid calling structural variants directly. Instead, the evidence for a potential structural variant at a locus is indicated by variation in the counts of short-reads that map anomalously to that locus. These structural variant traits are treated as quantitative traits and mapped genetically, analogously to a gene expression study. Association between a structural variant trait at one locus, and genotypes at a distant locus indicate the origin and target of a transposition. Using ultra-low-coverage (0.3×) population sequence data from 488 recombinant inbred Arabidopsis thaliana genomes, we identified 6502 segregating structural variants. Remarkably, 25% of these were transpositions. While many structural variants cannot be delineated precisely, we validated 83% of 44 predicted transposition breakpoints by polymerase chain reaction. We show that specific structural variants may be causative for quantitative trait loci for germination and resistance to infection by the fungus Albugo laibachii, isolate Nc14. Further we show that the phenotypic heritability attributable to read-mapping anomalies differs from, and, in the case of time to germination and bolting, exceeds that due to standard genetic variation. Genes within structural variants are also more likely to be silenced or dysregulated. This approach complements the prevalent strategy of structural variant discovery in fewer individuals sequenced at high coverage. It is generally applicable to large populations sequenced at low-coverage, and is particularly suited to mapping transpositions.

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.