Charles Robin

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.

Universal primers for fluorescent labelling of PCR fragments—an efficient and cost-effective approach to genotyping by fluorescence

Directly labelling locus‐specific primers for microsatellite analysis is expensive and a common limitation to small‐budget molecular ecology projects. More cost‐effective end‐labelling of PCR products can be achieved through a three primer PCR approach, involving a fluorescently labelled universal primer in combination with modified locus‐specific primers with 5′ universal primer sequence tails. This technique has been widely used but has been limited largely due to a lack of available universal primers suitable for co‐amplifying large numbers of size overlapping loci and without requiring locus‐specific PCR conditions to be modified. In this study, we report a suite of four high‐performance universal primers that can be employed in a three primer PCR approach for efficient and cost‐effective fluorescent end‐labelling of PCR fragments. Amplification efficiency is maximized owing to high universal primer Tm values (approximately 60+ °C) that enhance primer versatility and enable higher annealing temperatures to be employed compared with commonly used universal primers such as M13. We demonstrate that these universal primers can be combined with multiple fluorophores to co‐amplify multiple loci efficiently via multiplex PCR. This method provides a level of multiplexing and PCR efficiency similar to microsatellite fluorescent detection assays using directly labelled primers while dramatically reducing project costs. Primer performance is tested using several alternative PCR strategies that involve both single and multiple fluorophores in single and multiplex PCR across a wide range of taxa.

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.

A next-generation sequencing method for overcoming the multiple gene copy problem in polyploid phylogenetics, applied to Poa grasses

BackgroundPolyploidy is important from a phylogenetic perspective because of its immense past impact on evolution and its potential future impact on diversification, survival and adaptation, especially in plants. Molecular population genetics studies of polyploid organisms have been difficult because of problems in sequencing multiple-copy nuclear genes using Sanger sequencing. This paper describes a method for sequencing a barcoded mixture of targeted gene regions using next-generation sequencing methods to overcome these problems.ResultsUsing 64 3-bp barcodes, we successfully sequenced three chloroplast and two nuclear gene regions (each of which contained two gene copies with up to two alleles per individual) in a total of 60 individuals across 11 species of Australian Poa grasses. This method had high replicability, a low sequencing error rate (after appropriate quality control) and a low rate of missing data. Eighty-eight percent of the 320 gene/individual combinations produced sequence reads, and >80% of individuals produced sufficient reads to detect all four possible nuclear alleles of the homeologous nuclear loci with 95% probability.We applied this method to a group of sympatric Australian alpine Poa species, which we discovered to share an allopolyploid ancestor with a group of American Poa species. All markers revealed extensive allele sharing among the Australian species and so we recommend that the current taxonomy be re-examined. We also detected hypermutation in the trnH-psbA marker, suggesting it should not be used as a land plant barcode region. Some markers indicated differentiation between Tasmanian and mainland samples. Significant positive spatial genetic structure was detected at <100 km with chloroplast but not nuclear markers, which may be a result of restricted seed flow and long-distance pollen flow in this wind-pollinated group.ConclusionsOur results demonstrate that 454 sequencing of barcoded amplicon mixtures can be used to reliably sample all alleles of homeologous loci in polyploid species and successfully investigate phylogenetic relationships among species, as well as to investigate phylogeographic hypotheses. This next-generation sequencing method is more affordable than and at least as reliable as bacterial cloning. It could be applied to any experiment involving sequencing of amplicon mixtures.

Gene identification and proteomic analysis of the esterases of the cotton bollworm, Helicoverpa armigera.

Some of the resistance of Helicoverpa armigera to conventional insecticides such as organophosphates and synthetic pyrethroids appears to be due to metabolic detoxification by carboxylesterases. To investigate the H. armigera carboxyl/cholinesterases, we created a data set of 39 putative paralogous H. armigera carboxyl/cholinesterase sequences from cDNA libraries and other sources. Phylogenetic analysis revealed a close relationship between these sequences and 70 carboxyl/cholinesterases from the recently sequenced genome of the silkworm, Bombyx mori, including several conserved clades of non-catalytic proteins. A juvenile hormone esterase candidate from H. armigera was identified, and B. mori orthologues were proposed for 31% of the sequences examined, however low similarity was found between lepidopteran sequences and esterases previously associated with insecticide resistance from other insect orders. A proteomic analysis of larval esterases then enabled us to match seven of the H. armigera carboxyl/cholinesterase sequences to specific esterase isozymes. All identified sequences were predicted to encode catalytically active carboxylesterases, including six proteins with N-terminal signal peptides and N-glycans, with two also containing C-terminal signals for glycosylphosphatidylinositol anchor attachment. Five of these sequences were matched to zones of activity on native PAGE at relative mobility values previously associated with insecticide resistance in this species.

Duplication and divergence of the genes of the alpha-esterase cluster of Drosophila melanogaster.

The α-esterase cluster ofD. melanogaster contains 11 esterase genes dispersed over 60 kb. Embedded in the cluster are two unrelated open reading frames that have sequence similarity with genes encoding ubiquitin-conjugating enzyme and tropomyosin. The esterase amino acid sequences show 37–66% identity with one another and all but one have all the motifs characteristic of functional members of the carboxyl/cholinesterase multigene family. The exception has several frameshift mutations and appears to be a pseudogene. Patterns of amino acid differences among cluster members in relation to generic models of carboxyl/cholinesterase protein structure are broadly similar to those among other carboxyl/cholinesterases sequenced to date. However the α-esterases differ from most other members of the family in: their lack of a signal peptide; the lack of conservation in cysteines involved in disulfide bridges; and in four indels, two of which occur in or adjacent to regions that align with proposed substrate-binding sites of other carboxyl/cholinesterases. Phylogenetic analyses clearly identify three simple gene duplication events within the cluster. The most recent event involved the pseudogene which is located in an intron of another esterase gene. However, relative rate tests suggest that the pseudogene remained functional after the duplication event and has become inactive relatively recently. The distribution of indels also suggests a deeper node in the gene phylogeny that separates six genes at the two ends of the cluster from a block of five in the middle.

Gene duplication in the major insecticide target site, Rdl, in Drosophila melanogaster.

The Resistance to Dieldrin gene, Rdl, encodes a GABA-gated chloride channel subunit that is targeted by cyclodiene and phenylpyrazole insecticides. The gene was first characterized in Drosophila melanogaster by genetic mapping of resistance to the cyclodiene dieldrin. The 4,000-fold resistance observed was due to a single amino acid replacement, Ala301 to Ser. The equivalent change was subsequently identified in Rdl orthologs of a large range of resistant insect species. Here, we report identification of a duplication at the Rdl locus in D. melanogaster. The 113-kb duplication contains one WT copy of Rdl and a second copy with two point mutations: an Ala301 to Ser resistance mutation and Met360 to Ile replacement. Individuals with this duplication exhibit intermediate dieldrin resistance compared with single copy Ser301 homozygotes, reduced temperature sensitivity, and altered RNA editing associated with the resistant allele. Ectopic recombination between Roo transposable elements is involved in generating this genomic rearrangement. The duplication phenotypes were confirmed by construction of a transgenic, artificial duplication integrating the 55.7-kb Rdl locus with a Ser301 change into an Ala301 background. Gene duplications can contribute significantly to the evolution of insecticide resistance, most commonly by increasing the amount of gene product produced. Here however, duplication of the Rdl target site creates permanent heterozygosity, providing unique potential for adaptive mutations to accrue in one copy, without abolishing the endogenous role of an essential gene.

Identification of candidate genes for fenvalerate resistance in Helicoverpa armigera using cDNA-AFLP

A cDNA‐amplified fragment length polymorphisms approach was undertaken to screen for candidate genes associated with fenvalerate resistance in the AN02 strain of Helicoverpa armigera. Larvae and adults of this strain manifest ~50‐fold resistance, which is suppressible by piperonyl butoxide and controlled by the semidominant gene RFen1 previously mapped to AFLP Linkage Group 13. Two cytochrome P450s (CYP337B1 and CYP4S1), one carboxylesterase‐like protein and one glutathione transferase were found to be constitutively upregulated in resistant insects. Mapping of these potential detoxification genes showed that one of them, the novel P450 CYP337B1, was tightly linked to the resistance locus. This suggests that the RFen1R allele has a cis‐acting effect on CYP337B1 expression, and possible trans‐acting effects on expression of other genes.

Recognition and Detoxification of the Insecticide DDT by Drosophila melanogaster Glutathione S-Transferase D1

GSTD1 is one of several insect glutathione S-transferases capable of metabolizing the insecticide DDT. Here we use crystallography and NMR to elucidate the binding of DDT and glutathione to GSTD1. The crystal structure of Drosophila melanogaster GSTD1 has been determined to 1.1 A resolution, which reveals that the enzyme adopts the canonical GST fold but with a partially occluded active site caused by the packing of a C-terminal helix against one wall of the binding site for substrates. This helix would need to unwind or be displaced to enable catalysis. When the C-terminal helix is removed from the model of the crystal structure, DDT can be computationally docked into the active site in an orientation favoring catalysis. Two-dimensional (1)H,(15)N heteronuclear single-quantum coherence NMR experiments of GSTD1 indicate that conformational changes occur upon glutathione and DDT binding and the residues that broaden upon DDT binding support the predicted binding site. We also show that the ancestral GSTD1 is likely to have possessed DDT dehydrochlorinase activity because both GSTD1 from D. melanogaster and its sibling species, Drosophila simulans, have this activity.

Dissecting the Insect Metabolic Machinery Using Twin Ion Mass Spectrometry: A Single P450 Enzyme Metabolizing the Insecticide Imidacloprid in Vivo

Insecticide resistance is one of the most prevalent examples of anthropogenic genetic change, yet our understanding of metabolic-based resistance remains limited by the analytical challenges associated with rapidly tracking the in vivo metabolites of insecticides at nonlethal doses. Here, using twin ion mass spectrometry analysis of the extracts of whole Drosophila larvae and excreta, we show that (i) eight metabolites of the neonicotinoid insecticide, imidacloprid, can be detected when formed by susceptible larval genotypes and (ii) the specific overtranscription of a single gene product, Cyp6g1, associated with the metabolic resistance to neonicotinoids, results in a significant increase in the formation of three imidacloprid metabolites that are formed in C-H bond activation reactions; that is, Cyp6g1 is directly linked to the enhanced metabolism of imidacloprid in vivo. These results establish a rapid and sensitive method for dissecting the metabolic machinery of insects by directly linking single gene products to insecticide metabolism.