Philip Batterham

Copy number variation and transposable elements feature in recent, ongoing adaptation at the Cyp6g1 locus.

The increased transcription of the Cyp6g1 gene of Drosophila melanogaster, and consequent resistance to insecticides such as DDT, is a widely cited example of adaptation mediated by cis-regulatory change. A fragment of an Accord transposable element inserted upstream of the Cyp6g1 gene is causally associated with resistance and has spread to high frequencies in populations around the world since the 1940s. Here we report the existence of a natural allelic series at this locus of D. melanogaster, involving copy number variation of Cyp6g1, and two additional transposable element insertions (a P and an HMS-Beagle). We provide evidence that this genetic variation underpins phenotypic variation, as the more derived the allele, the greater the level of DDT resistance. Tracking the spatial and temporal patterns of allele frequency changes indicates that the multiple steps of the allelic series are adaptive. Further, a DDT association study shows that the most resistant allele, Cyp6g1-[BP], is greatly enriched in the top 5% of the phenotypic distribution and accounts for ∼16% of the underlying phenotypic variation in resistance to DDT. In contrast, copy number variation for another candidate resistance gene, Cyp12d1, is not associated with resistance. Thus the Cyp6g1 locus is a major contributor to DDT resistance in field populations, and evolution at this locus features multiple adaptive steps occurring in rapid succession.

Characterization of Drosophila melanogaster cytochrome P450 genes

Cytochrome P450s form a large and diverse family of heme-containing proteins capable of carrying out many different enzymatic reactions. In both mammals and plants, some P450s are known to carry out reactions essential for processes such as hormone synthesis, while other P450s are involved in the detoxification of environmental compounds. In general, functions of insect P450s are less well understood. We characterized Drosophila melanogaster P450 expression patterns in embryos and 2 stages of third instar larvae. We identified numerous P450s expressed in the fat body, Malpighian (renal) tubules, and in distinct regions of the midgut, consistent with hypothesized roles in detoxification processes, and other P450s expressed in organs such as the gonads, corpora allata, oenocytes, hindgut, and brain. Combining expression pattern data with an RNA interference lethality screen of individual P450s, we identify candidate P450s essential for developmental processes and distinguish them from P450s with potential functions in detoxification.

The genetic, molecular and phenotypic consequences of selection for insecticide resistance

Studies of insecticide resistance allow theories of the adaptive process to be tested where the selective agent, the insecticide, is unambiguously defined. Thus, the consequences of selection of phenotypic variation can be investigated in genetic, biochemical, molecular, population biological and, most recently, developmental contexts. Are the options limited biochemically and molecularly? Is the genetic mechanism monogenic or polygenic, general or population/species specific? Are fitness and developmental patterns associated? These questions of general evolutionary significance can be considered with experimental approaches to determine how insecticide resistance evolves.

Quantitative trait symmetry independent of Hsp90 buffering: Distinct modes of genetic canalization and developmental stability

The Hsp90 chaperone buffers development against a wide range of morphological changes in many organisms and in Drosophila masks the effects of hidden genetic variation. Theory predicts that genetic and nongenetic buffering will share common mechanisms. For example, it is argued that Hsp90 genetic buffering evolved solely as a by-product of environmental buffering, and that Hsp90 should mask morphological deviations from any source. To test this idea, we examined the effect of Hsp90 on purely nongenetic variation in phenotype, measured as differences between the left and right sides of several bilaterally symmetrical bristle and wing traits in individual flies. Consistent with previous reports, Hsp90 buffered the expression of rare morphogenic variants specific to particular genetic backgrounds. However, neither trait-by-trait nor global asymmetry was affected in outbred flies treated with an Hsp90 inhibitor or across a series of inbred genetic backgrounds from a wild population tested in isogenic F1 heterozygotes carrying either (i) a dominant negative Hsp90 allele on a mutant 3rd chromosome or (ii) a null P-insertion mutation, which was introgressed into the control genetic background on all chromosomes. By contrast, Hsp90-regulated trait means and significant effects of sex, temperature, and genetic background on trait symmetry were clearly detected. We conclude that, by maintaining the function of signaling proteins, Hsp90 masks variation affecting target pathways and traits in populations independent of purely nongenetic sources of variation, refuting the idea that a single Hsp90-dependent process generally controls genetic canalization and developmental stability.

The acetylcholinesterase gene and organophosphorus resistance in the Australian sheep blowfly, Lucilia cuprina

Acetylcholinesterase (AChE), encoded by the Ace gene, is the primary target of organophosphorous (OP) and carbamate insecticides. Ace mutations have been identified in OP resistants strains of Drosophila melanogaster. However, in the Australian sheep blowfly, Lucilia cuprina, resistance in field and laboratory generated strains is determined by point mutations in the Rop-1 gene, which encodes a carboxylesterase, E3. To investigate the apparent bias for the Rop-1/E3 mechanism in the evolution of OP resistance in L. cuprina, we have cloned the Ace gene from this species and characterized its product. Southern hybridization indicates the existence of a single Ace gene in L. cuprina. The amino acid sequence of L. cuprina AChE shares 85.3% identity with D. melanogaster and 92.4% with Musca domestica AChE. Five point mutations in Ace associated with reduced sensitivity to OP insecticides have been previously detected in resistant strains of D. melanogaster. These residues are identical in susceptible strains of D. melanogaster and L. cuprina, although different codons are used. Each of the amino acid substitutions that confer OP resistance in D. melanogaster could also occur in L. cuprina by a single non-synonymous substitution. These data suggest that the resistance mechanism used in L. cuprina is determined by factors other than codon bias. The same point mutations, singly and in combination, were introduced into the Ace gene of L. cuprina by site-directed mutagenesis and the resulting AChE enzymes expressed using a baculovirus system to characterise their kinetic properties and interactions with OP insecticides. The K(m) of wild type AChE for acetylthiocholine (ASCh) is 23.13 microM and the point mutations change the affinity to the substrate. The turnover number of Lucilia AChE for ASCh was estimated to be 1.27x10(3) min(-1), similar to Drosophila or housefly AChE. The single amino acid replacements reduce the affinities of the AChE for OPs and give up to 8.7-fold OP insensitivity, while combined mutations give up to 35-fold insensitivity. However, other published studies indicate these same mutations yield higher levels of OP insensitivity in D. melanogaster and A. aegypti. The inhibition data indicate that the wild type form of AChE of L. cuprina is 12.4-fold less sensitive to OP inhibition than the susceptible form of E3, suggesting that the carboxylesterases may have a role in the protection of AChE via a sequestration mechanism. This provides a possible explanation for the bias towards the evolution of resistance via the Rop-1/E3 mechanism in L. cuprina.

Mitochondrial DNA analysis of field populations of Helicoverpa armigera (Lepidoptera: Noctuidae) and of its relationship to H. zea

BackgroundHelicoverpa armigera and H. zea are amongst the most significant polyphagous pest lepidopteran species in the Old and New Worlds respectively. Separation of H. armigera and H. zea is difficult and is usually only achieved through morphological differences in the genitalia. They are capable of interbreeding to produce fertile offspring. The single species status of H. armigera has been doubted, due to its wide distribution and plant host range across the Old World. This study explores the global genetic diversity of H. armigera and its evolutionary relationship to H zea.ResultsWe obtained partial (511 bp) mitochondrial DNA (mtDNA) Cytochrome Oxidase-I (COI) sequences for 249 individuals of H. armigera sampled from Australia, Burkina Faso, Uganda, China, India and Pakistan which were associated with various host plants. Single nucleotide polymorphisms (SNPs) within the partial COI gene differentiated H. armigera populations into 33 mtDNA haplotypes. Shared haplotypes between continents, low F-statistic values and low nucleotide diversity between countries (0.0017 – 0.0038) suggests high mobility in this pest. Phylogenetic analysis of four major Helicoverpa pest species indicates that H. punctigera is basal to H. assulta, which is in turn basal to H. armigera and H. zea. Samples from North and South America suggest that H. zea is also a single species across its distribution. Our data reveal short genetic distances between H. armigera and H. zea which seem to have been established via a founder event from H. armigera stock at around 1.5 million years ago.ConclusionOur mitochondrial DNA sequence data supports the single species status of H. armigera across Africa, Asia and Australia. The evidence for inter-continental gene flow observed in this study is consistent with published evidence of the capacity of this species to migrate over long distances. The finding of high genetic similarity between Old World H. armigera and New World H. zea emphasises the need to consider work on both pests when building pest management strategies for either.

Molecular Evolution of Glutathione S-Transferases in the Genus Drosophila

As classical phase II detoxification enzymes, glutathione S-transferases (GSTs) have been implicated in insecticide resistance and may have evolved in response to toxins in the niche-defining feeding substrates of Drosophila species. We have annotated the GST genes of the 12 Drosophila species with recently sequenced genomes and analyzed their molecular evolution. Gene copy number variation is attributable mainly to unequal crossing-over events in the large δ and ε clusters. Within these gene clusters there are also GST genes with slowly diverging orthologs. This implies that they have their own unique functions or have spatial/temporal expression patterns that impose significant selective constraints. Searches for positively selected sites within the GSTs identified G171K in GSTD1, a protein that has previously been shown to be capable of metabolizing the insecticide DDT. We find that the same radical substitution (G171K) in the substrate-binding domain has occurred at least three times in the Drosophila radiation. Homology-modeling places site 171 distant from the active site but adjacent to an alternative DDT-binding site. We propose that the parallel evolution observed at this site is an adaptive response to an environmental toxin and that sequencing of historical alleles suggests that this toxin was not a synthetic insecticide.

Beyond barcoding : a mitochondrial genomics approach to molecular phylogenetics and diagnostics of blowflies (Diptera: Calliphoridae)

Members of the Calliphoridae (blowflies) are significant for medical and veterinary management, due to the ability of some species to consume living flesh as larvae, and for forensic investigations due to the ability of others to develop in corpses. Due to the difficulty of accurately identifying larval blowflies to species there is a need for DNA-based diagnostics for this family, however the widely used DNA-barcoding marker, cox1, has been shown to fail for several groups within this family. Additionally, many phylogenetic relationships within the Calliphoridae are still unresolved, particularly deeper level relationships. Sequencing whole mt genomes has been demonstrated both as an effective method for identifying the most informative diagnostic markers and for resolving phylogenetic relationships. Twenty-seven complete, or nearly so, mt genomes were sequenced representing 13 species, seven genera and four calliphorid subfamilies and a member of the related family Tachinidae. PCR and sequencing primers developed for sequencing one calliphorid species could be reused to sequence related species within the same superfamily with success rates ranging from 61% to 100%, demonstrating the speed and efficiency with which an mt genome dataset can be assembled. Comparison of molecular divergences for each of the 13 protein-coding genes and 2 ribosomal RNA genes, at a range of taxonomic scales identified novel targets for developing as diagnostic markers which were 117-200% more variable than the markers which have been used previously in calliphorids. Phylogenetic analysis of whole mt genome sequences resulted in much stronger support for family and subfamily-level relationships. The Calliphoridae are polyphyletic, with the Polleninae more closely related to the Tachinidae, and the Sarcophagidae are the sister group of the remaining calliphorids. Within the Calliphoridae, there was strong support for the monophyly of the Chrysomyinae and Luciliinae and for the sister-grouping of Luciliinae with Calliphorinae. Relationships within Chrysomya were not well resolved. Whole mt genome data, supported the previously demonstrated paraphyly of Lucilia cuprina with respect to L. sericata and allowed us to conclude that it is due to hybrid introgression prior to the last common ancestor of modern sericata populations, rather than due to recent hybridisation, nuclear pseudogenes or incomplete lineage sorting.

Copper homoeostasis in Drosophila melanogaster S2 cells

Copper homoeostasis was investigated in the Drosophila melanogaster S2 cell line to develop an insect model for the study of copper regulation. Real-time PCR studies have demonstrated expression in S2 cells of putative orthologues of human Cu regulatory genes involved in the uptake, transport, sequestration and efflux of Cu. Drosophila orthologues of the mammalian Cu chaperones, ATOX1 (a human orthologue of yeast ATX1), CCS (copper chaperone for superoxide dismutase), COX17 (a human orthologue of yeast COX17), and SCO1 and SCO2, did not significantly respond transcriptionally to increased Cu levels, whereas MtnA, MtnB and MtnD (Drosophila orthologues of human metallothioneins) were up-regulated by Cu in a time- and dose-dependent manner. To examine the effect on Cu homoeostasis, expression of several key copper homoeostasis genes was suppressed using double-stranded RNA interference. Suppression of the MTF-1 (metal-regulatory transcription factor 1), reduced both basal and Cu-induced gene expressions of MtnA, MtnB and MtnD, significantly reducing the tolerance of these cells to increased Cu. Suppression of either Ctr1A (a Drosophila orthologue of yeast CTR1) or Ctr1B significantly reduced Cu uptake from media, demonstrating that both these proteins function to transport Cu into S2 cells. Significantly, Cu induced Ctr1B gene expression, and this could be prevented by suppressing MTF-1, suggesting that Ctr1B might be involved in Cu detoxification. Suppression of DmATP7, the putative homologue of human Cu transporter genes ATP7A and ATP7B, significantly increased Cu accumulation, demonstrating that DmATP7 is essential for efflux of excess Cu. This work is consistent with previous studies in mammalian cells, validating S2 cells as a model system for studying Cu transport and identifying novel Cu regulatory mechanisms.

Developing compound eye in lozenge mutants of Drosophila: lozenge expression in the R7 equivalence group

Abstract The lozenge locus is genetically complex, containing two functionally distinct units, cistrons A and B, that influence the structure of the compound eye. Extreme mutations of either cistron produce adult phenotypes that share similarities and that have striking differences. We have analyzed the expression of several developmentally important eye genes including boss, scabrous, rhomboid, seven-up, and Bar in lozenge mutant backgrounds representing both cistrons. This analysis follows the progressive recruitment of photoreceptor neurons during eye development and has confirmed that the initial development of photoreceptors is normal up to the five cell precluster stage (R8, R2/5 and R3/4). However, when lozenge is mutant, further eye development is perturbed. As cells R1, R6 and R7 are recruited, patterns of gene expression for seven-up and Bar become abnormal. We have also characterized the expression of two different enhancer trap alleles of lozenge. The lozenge product(s) appear to be first expressed in the eye disc in undifferentiated cells shortly after the five cell precluster forms. Then, as distinct cells are recruited to a fate, lozenge expression persists and is refined in those cells. Our data suggests that lozenge functions in cone cells and pigment cells as well as in specific glia. With respect to photoreceptor neurons, lozenge biases the developmental potential of cells R1, R6 and R7, by directly influencing the expression of genes important for establishing cell fate.