Edward J. Osborne

The genome of Tetranychus urticae reveals herbivorous pest adaptations

The spider mite Tetranychus urticae is a cosmopolitan agricultural pest with an extensive host plant range and an extreme record of pesticide resistance. Here we present the completely sequenced and annotated spider mite genome, representing the first complete chelicerate genome. At 90 megabases T. urticae has the smallest sequenced arthropod genome. Compared with other arthropods, the spider mite genome shows unique changes in the hormonal environment and organization of the Hox complex, and also reveals evolutionary innovation of silk production. We find strong signatures of polyphagy and detoxification in gene families associated with feeding on different hosts and in new gene families acquired by lateral gene transfer. Deep transcriptome analysis of mites feeding on different plants shows how this pest responds to a changing host environment. The T. urticae genome thus offers new insights into arthropod evolution and plant–herbivore interactions, and provides unique opportunities for developing novel plant protection strategies.

Multiple reference genomes and transcriptomes for Arabidopsis thaliana

Genetic differences between Arabidopsis thaliana accessions underlie the plant’s extensive phenotypic variation, and until now these have been interpreted largely in the context of the annotated reference accession Col-0. Here we report the sequencing, assembly and annotation of the genomes of 18 natural A. thaliana accessions, and their transcriptomes. When assessed on the basis of the reference annotation, one-third of protein-coding genes are predicted to be disrupted in at least one accession. However, re-annotation of each genome revealed that alternative gene models often restore coding potential. Gene expression in seedlings differed for nearly half of expressed genes and was frequently associated with cis variants within 5 kilobases, as were intron retention alternative splicing events. Sequence and expression variation is most pronounced in genes that respond to the biotic environment. Our data further promote evolutionary and functional studies in A. thaliana, especially the MAGIC genetic reference population descended from these accessions.

DNA methylation in Arabidopsis has a genetic basis and shows evidence of local adaptation

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. DOI: http://dx.doi.org/10.7554/eLife.05255.001

Population bulk segregant mapping uncovers resistance mutations and the mode of action of a chitin synthesis inhibitor in arthropods

Because of its importance to the arthropod exoskeleton, chitin biogenesis is an attractive target for pest control. This point is demonstrated by the economically important benzoylurea compounds that are in wide use as highly specific agents to control insect populations. Nevertheless, the target sites of compounds that inhibit chitin biogenesis have remained elusive, likely preventing the full exploitation of the underlying mode of action in pest management. Here, we show that the acaricide etoxazole inhibits chitin biogenesis in Tetranychus urticae (the two-spotted spider mite), an economically important pest. We then developed a population-level bulk segregant mapping method, based on high-throughput genome sequencing, to identify a locus for monogenic, recessive resistance to etoxazole in a field-collected population. As supported by additional genetic studies, including sequencing across multiple resistant strains and genetic complementation tests, we associated a nonsynonymous mutation in the major T. urticae chitin synthase (CHS1) with resistance. The change is in a C-terminal transmembrane domain of CHS1 in a highly conserved region that may serve a noncatalytic but essential function. Our finding of a target-site resistance mutation in CHS1 shows that at least one highly specific chitin biosynthesis inhibitor acts directly to inhibit chitin synthase. Our work also raises the possibility that other chitin biogenesis inhibitors, such as the benzoylurea compounds, may also act by inhibition of chitin synthases. More generally, our genetic mapping approach should be powerful for high-resolution mapping of simple traits (resistance or otherwise) in arthropods.

A burst of ABC genes in the genome of the polyphagous spider mite Tetranychus urticae

BackgroundThe ABC (ATP-binding cassette) gene superfamily is widespread across all living species. The majority of ABC genes encode ABC transporters, which are membrane-spanning proteins capable of transferring substrates across biological membranes by hydrolyzing ATP. Although ABC transporters have often been associated with resistance to drugs and toxic compounds, within the Arthropoda ABC gene families have only been characterized in detail in several insects and a crustacean. In this study, we report a genome-wide survey and expression analysis of the ABC gene superfamily in the spider mite, Tetranychus urticae, a chelicerate ~ 450 million years diverged from other Arthropod lineages. T. urticae is a major agricultural pest, and is among of the most polyphagous arthropod herbivores known. The species resists a staggering array of toxic plant secondary metabolites, and has developed resistance to all major classes of pesticides in use for its control.ResultsWe identified 103 ABC genes in the T. urticae genome, the highest number discovered in a metazoan species to date. Within the T. urticae ABC gene set, all members of the eight currently described subfamilies (A to H) were detected. A phylogenetic analysis revealed that the high number of ABC genes in T. urticae is due primarily to lineage-specific expansions of ABC genes within the ABCC, ABCG and ABCH subfamilies. In particular, the ABCC subfamily harbors the highest number of T. urticae ABC genes (39). In a comparative genomic analysis, we found clear orthologous relationships between a subset of T. urticae ABC proteins and ABC proteins in both vertebrates and invertebrates known to be involved in fundamental cellular processes. These included members of the ABCB-half transporters, and the ABCD, ABCE and ABCF families. Furthermore, one-to-one orthologues could be distinguished between T. urticae proteins and human ABCC10, ABCG5 and ABCG8, the Drosophila melanogaster sulfonylurea receptor and ecdysone-regulated transporter E23. Finally, expression profiling revealed that ABC genes in the ABCC, ABCG ABCH subfamilies were differentially expressed in multi-pesticide resistant mite strains and/or in mites transferred to challenging (toxic) host plants.ConclusionsIn this study we present the first comprehensive analysis of ABC genes in a polyphagous arthropod herbivore. We demonstrate that the broad plant host range and high levels of pesticide resistance in T. urticae are associated with lineage-specific expansions of ABC genes, many of which respond transcriptionally to xenobiotic exposure. This ABC catalogue will serve as a basis for future biochemical and toxicological studies. Obtaining functional evidence that these ABC subfamilies contribute to xenobiotic tolerance should be the priority of future research.

Reciprocal Responses in the Interaction between Arabidopsis and the Cell-Content-Feeding Chelicerate Herbivore Spider Mite

Analysis of genome-wide Arabidopsis responses to spider mite identifies the secondary metabolites that limit the mites ability to use Arabidopsis as a host. Most molecular-genetic studies of plant defense responses to arthropod herbivores have focused on insects. However, plant-feeding mites are also pests of diverse plants, and mites induce different patterns of damage to plant tissues than do well-studied insects (e.g. lepidopteran larvae or aphids). The two-spotted spider mite (Tetranychus urticae) is among the most significant mite pests in agriculture, feeding on a staggering number of plant hosts. To understand the interactions between spider mite and a plant at the molecular level, we examined reciprocal genome-wide responses of mites and its host Arabidopsis (Arabidopsis thaliana). Despite differences in feeding guilds, we found that transcriptional responses of Arabidopsis to mite herbivory resembled those observed for lepidopteran herbivores. Mutant analysis of induced plant defense pathways showed functionally that only a subset of induced programs, including jasmonic acid signaling and biosynthesis of indole glucosinolates, are central to Arabidopsis’s defense to mite herbivory. On the herbivore side, indole glucosinolates dramatically increased mite mortality and development times. We identified an indole glucosinolate dose-dependent increase in the number of differentially expressed mite genes belonging to pathways associated with detoxification of xenobiotics. This demonstrates that spider mite is sensitive to Arabidopsis defenses that have also been associated with the deterrence of insect herbivores that are very distantly related to chelicerates. Our findings provide molecular insights into the nature of, and response to, herbivory for a representative of a major class of arthropod herbivores.

Divergence, convergence, and the ancestry of feral populations in the domestic rock pigeon.

Domestic pigeons are spectacularly diverse and exhibit variation in more traits than any other bird species [1]. In The Origin of Species, Charles Darwin repeatedly calls attention to the striking variation among domestic pigeon breeds-generated by thousands of years of artificial selection on a single species by human breeders-as a model for the process of natural divergence among wild populations and species [2]. Darwin proposed a morphology-based classification of domestic pigeon breeds [3], but the relationships among major groups of breeds and their geographic origins remain poorly understood [4, 5]. We used a large, geographically diverse sample of 361 individuals from 70 domestic pigeon breeds and two free-living populations to determine genetic relationships within this species. We found unexpected relationships among phenotypically divergent breeds as well as convergent evolution of derived traits among several breed groups. Our findings also illuminate the geographic origins of breed groups in India and the Middle East and suggest that racing breeds have made substantial contributions to feral pigeon populations.

Limited Contribution of DNA Methylation Variation to Expression Regulation in Arabidopsis thaliana

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.

Wham: Identifying Structural Variants of Biological Consequence.

Existing methods for identifying structural variants (SVs) from short read datasets are inaccurate. This complicates disease-gene identification and efforts to understand the consequences of genetic variation. In response, we have created Wham (Whole-genome Alignment Metrics) to provide a single, integrated framework for both structural variant calling and association testing, thereby bypassing many of the difficulties that currently frustrate attempts to employ SVs in association testing. Here we describe Wham, benchmark it against three other widely used SV identification tools–Lumpy, Delly and SoftSearch–and demonstrate Wham’s ability to identify and associate SVs with phenotypes using data from humans, domestic pigeons, and vaccinia virus. Wham and all associated software are covered under the MIT License and can be freely downloaded from github (https://github.com/zeeev/wham), with documentation on a wiki (http://zeeev.github.io/wham/). For community support please post questions to https://www.biostars.org/.

Genomic Rearrangements in Arabidopsis Considered as Quantitative Traits

Structural Rearrangements can have unexpected effects on quantitative phenotypes. Surprisingly, these rearrangements can also be considered as... 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.