Quantitative trait locus mapping and functional genomics of an organophosphate resistance trait in the western corn rootworm, Diabrotica virgifera virgifera.

The western corn rootworm, Diabrotica virgifera virgifera, is an insect pest of corn and population suppression with chemical insecticides is an important management tool. Traits conferring organophosphate insecticide resistance have increased in frequency amongst D. v. virgifera populations, resulting in the reduced efficacy in many corn-growing regions of the USA. We used comparative functional genomic and quantitative trait locus (QTL) mapping approaches to investigate the genetic basis of D. v. virgifera resistance to the organophosphate methyl-parathion. RNA from adult methyl-parathion resistant and susceptible adults was hybridized to 8331 microarray probes. The results predicted that 11 transcripts were significantly up-regulated in resistant phenotypes, with the most significant (fold increases ≥ 2.43) being an α-esterase-like transcript. Differential expression was validated only for the α-esterase (ST020027A20C03), with 11- to 13-fold greater expression in methyl-parathion resistant adults (P < 0.05). Progeny with a segregating methyl-parathion resistance trait were obtained from a reciprocal backcross design. QTL analyses of high-throughput single nucleotide polymorphism genotype data predicted involvement of a single genome interval. These data suggest that a specific carboyxesterase may function in field-evolved corn rootworm resistance to organophosphates, even though direct linkage between the QTL and this locus could not be established.

Odorant receptors of a primitive hymenopteran pest, the wheat stem sawfly.

The wheat stem sawfly, Cephus cinctus, is an herbivorous hymenopteran that feeds exclusively on members of the Graminae family. Synanthropically, it has become one of the most important insect pests of wheat grown in the northern Great Plains region of the USA and Canada. Insecticides are generally ineffective because of the wheat stem sawfly's extended adult flight period and its inaccessible larval stage, during which it feeds within the wheat stems, making it virtually intractable to most pest management strategies. While research towards integrated pest management strategies based on insect olfaction has proved promising, nothing is known about the molecular basis of olfaction in this important pest species. In this study we identified 28 unique odorant receptor (Or) transcripts from an antennal transcriptome. A phylogenetic analysis with the predicted Ors from the honey bee and jewel wasp genomes revealed at least four clades conserved amongst all three Hymenoptera species. Antennal expression levels were analysed using quantitative real-time PCR, and one male-biased and five female-biased Ors were identified. This study provides the basis for future functional analyses to identify behaviourally active odours that can be used to help develop olfactory-mediated pest management tools.

A functional DNA methylation system in the pea aphid, Acyrthosiphon pisum.

Methylation of cytosine is one of the main epigenetic mechanisms involved in controlling gene expression. Here we show that the pea aphid (Acyrthosiphon pisum) genome possesses homologues to all the DNA methyltransferases found in vertebrates, and that 0.69% (+/-0.25%) of all cytosines are methylated. Identified methylation sites are predominantly restricted to the coding sequence of genes at CpG sites. We identify twelve methylated genes, including genes that interact with juvenile hormone, a key endocrine signal in insects. Bioinformatic prediction using CpG ratios for all predicted genes suggest that a large proportion of genes are methylated within the pea aphid.

The insect chemoreceptor superfamily of the parasitoid jewel wasp Nasonia vitripennis.

Chemoreception is important for locating food, mates and other resources in many insects, including the parasitoid jewel wasp Nasonia vitripennis. In the insect chemoreceptor superfamily, Nasonia has 58 gustatory receptor (Gr) genes, of which 11 are pseudogenes, leaving 47 apparently intact proteins encoded. No carbon dioxide receptors, two candidate sugar receptors, a DmGr43a orthologue, and several additional Gr lineages were identified, including significant gene subfamily expansions related to the 10 Grs found in the honey bee Apis mellifera. Nasonia has a total of 301 odorant receptor (Or) genes, of which 76 are pseudogenes, leaving 225 apparently intact Ors. Phylogenetic comparison with the 174 honey bee Ors reveals differential gene subfamily expansion in each hymenopteran lineage, along with a few losses from each species. The only simple orthologous relationship is the expected single DmOr83b orthologue. The large number of Nasonia Ors is the result of several major subfamily expansions, including one of 55 genes. Nasonia does not have the elaborate social chemical communication of honey bees, nor the diversity of floral odours honey bees detect, however, Nasonia wasps might need to detect a diversity of odours to find potential mates and hosts or avoid harmful substances in its environment.

The gustatory receptor family in the silkworm moth Bombyx mori is characterized by a large expansion of a single lineage of putative bitter receptors.

The gustatory receptor (Gr) family of insect chemoreceptors includes receptors for sugars and bitter compounds, as well as cuticular hydrocarbons and odorants such as carbon dioxide. We have annotated a total of 65 Gr genes from the silkworm Bombyx mori genome. The Gr family in the silkworm moth includes putative carbon dioxide receptors and sugar receptors, as well as duplicated orthologues of the orphan DmGr43a receptor. Most prominent in this 65-gene family, however, is a single large expansion of 55 Grs that we propose are predominantly 'bitter' receptors involved in perception of the large variety of secondary plant chemicals that caterpillars and moths encounter. These Grs might therefore mediate food choice and avoidance as well as oviposition site preference.

Molecular characterization of the Aedes aegypti odorant receptor gene family.

The olfactory-driven blood-feeding behaviour of female Aedes aegypti mosquitoes is the primary transmission mechanism by which the arboviruses causing dengue and yellow fevers affect over 40 million individuals worldwide. Bioinformatics analysis has been used to identify 131 putative odourant receptors from the A. aegypti genome that are likely to function in chemosensory perception in this mosquito. Comparison with the Anopheles gambiae olfactory subgenome demonstrates significant divergence of the odourant receptors that reflects a high degree of evolutionary activity potentially resulting from their critical roles during the mosquito life cycle. Expression analyses in the larval and adult olfactory chemosensory organs reveal that the ratio of odourant receptors to antennal glomeruli is not necessarily one to one in mosquitoes.

Female-biased expression of odourant receptor genes in the adult antennae of the silkworm, Bombyx mori.

Olfaction plays an important role in the life history of insects, including key behaviours such as host selection, oviposition and mate recognition. Odour perception by insects is primarily mediated by the large diverse family of odourant receptors (Ors) that are expressed on the dendrites of olfactory neurones housed within chemosensilla. However, few Or sequences have been identified from the Lepidoptera, an insect order that includes some of the most important pest species worldwide. We have identified 41 Or gene sequences from the silkworm (Bombyx mori) genome, more than double the number of published Or sequences from the Lepidoptera. Many silkworm Ors appear to be orthologs of the 17 published tobacco budworm (Heliothis virescens) Ors indicating that many Or lineages may be conserved within the Lepidoptera. The majority of the Or genes are expressed in adult female and male antennae (determined by quantitative real-time PCR analysis), supporting their probable roles in adult olfaction. Several Or genes are expressed at high levels in both male and female antennae, suggesting they mediate the perception of common host or conspecific volatiles important to both sexes. BmOrs 45-47 group together in the same phylogenetic branch and all three are expressed at moderate female-biased ratios, six to eight times higher in female compared to male moth antennae. Interestingly, BmOrs19 and 30 appear to be expressed predominantly in female antennae, opposite to that of the published silkworm pheromone receptors BmOrs 1 and 3 that are specific to male antennae. These results suggest that BmOr19 and 30 may detect odours critical to female behaviour, such as oviposition cues or male-produced courtship pheromones.

Molecular evolutionary analysis of the widespread piggyBac transposon family and related "domesticated" sequences.

piggyBac is a short inverted-repeat-type DNA transposable element originally isolated from the genome of the moth Trichoplusia ni. It is currently the gene vector of choice for the transformation of various insect species. A few sequences with similarity to piggyBac have previously been identified from organisms such as humans ( Looper), the pufferfish Takifugu rubripes ( Pigibaku), Xenopus ( Tx), Daphnia ( Pokey), and the Oriental fruit fly Bactrocera dorsalis. We have now identified 50 piggyBac-like sequences from publicly available genome sequences and expressed sequence tags (ESTs). This survey allows the first comparative examination of the distinctive piggyBac transposase, suggesting that it might contain a highly divergent DDD domain, comparable to the widespread DDE domain found in many DNA transposases and retroviral integrases which consists of two absolutely conserved aspartic acids separated by about 70 amino acids with a highly conserved glutamic acid about 35 amino acids further away. Many piggyBac-like sequences were found in the genomes of a phylogenetically diverse range of organisms including fungi, plants, insects, crustaceans, urochordates, amphibians, fishes and mammals. Also, several instances of "domestication" of the piggyBac transposase sequence by the host genome for cellular functions were identified. Novel members of the piggyBac family may be useful in genetic engineering of many organisms.

Bmmar6, a second mori subfamily mariner transposon from the silkworm moth Bombyx mori.

A second member of the divergent mori subfamily of mariner transposons, Bmmar6, is described from the silkworm moth Bombyx mori genome. A confident consensus sequence for Bmmar6 was obtained from a single genomic copy, 17 EST sequences, and the direct sequencing of a 'family' sequence from an amplification of all full-length genomic copies. Bmmar6 is most similar to Bmmar1 in the mori subfamily, which now also includes several fly and nematode transposons. These might be viewed as a discrete family of transposons within the IS630-Tc1-mariner superfamily with a distinctive D,D37D catalytic motif, and another small divergent D,D41D clade is recognized as their sister group of transposons.

Calcium channel gamma subunits provide insights into the evolution of this gene family.

The gamma subunits of voltage-dependent calcium channels influence calcium current properties and may be involved in other physiological functions. Five distinct gamma subunits have been described from human and/or mouse. The first identified member of this group of proteins, gamma(1), is a component of the L-type calcium channel expressed in skeletal muscle. A second member, gamma(2), identified from the stargazer mouse regulates the targeting of alpha-amino-3-hydroxy-5-methylisoxazole-4-propionic acid (AMPA) receptors to the postsynaptic membrane. We report here the identification of three novel gamma subunits from rat and mouse as well as the unidentified rat, mouse and human orthologs of the previously described subunits. Phylogenetic analysis of the 24 mammalian gamma subunits suggests the following relationship ((((gamma(2), gamma(3)), (gamma(4), gamma(8))), (gamma(5), gamma(7))), (gamma(1), gamma(6))) that indicates that they evolved from a common ancestral gamma subunit via gene duplication. Our analysis reveals that the novel gamma subunit gamma(6) most closely resembles gamma(1) and shares with it the lack of a PSD-95/DLG/ZO-1 (PDZ)-binding motif that is characteristic of most other gamma subunits. Rat gamma subunit mRNAs are expressed in multiple tissues including brain, heart, lung, and testis. The expression of gamma(1) mRNA and the long isoform of gamma(6) mRNA is most robust in skeletal muscle, while gamma(6) is also highly expressed in cardiac muscle. Based on our analysis of the molecular evolution, primary structure, and tissue distribution of the gamma subunits, we propose that gamma(1) and gamma(6) may share common physiological functions distinct from the other homologous gamma subunits.

Candidate odorant receptors from the malaria vector mosquito Anopheles gambiae and evidence of down-regulation in response to blood feeding.

Olfaction plays a major role in host preference and blood feeding, integral behaviors for disease transmission by the malaria vector mosquito Anopheles gambiae sensu stricto (henceforth A. gambiae). We have identified four genes encoding candidate odorant receptors from A. gambiae that are selectively expressed in olfactory organs, contain approximately seven transmembrane domains, and show significant similarity to several putative odorant receptors in Drosophila melanogaster. Furthermore, one of the putative A. gambiae odorant receptors exhibits female-specific antennal expression and is down-regulated 12 h after blood feeding, a period during which substantial reduction in olfactory responses to human odorants has been observed. Taken together, these data suggest these genes encode a family of odorant receptors in A. gambiae, whose further study may aid in the design of novel antimalarial programs.

Taste: independent origins of chemoreception coding systems?

A large family of divergent candidate gustatory receptors has been identified in Drosophila. As with the odorant receptors, one receptor is expressed per sensory neuron, each class of which projects to discrete regions of the brain, allowing a combinatorial coding system for specific recognition of ligands.

Loss of transposase-DNA interaction may underlie the divergence of mariner family transposable elements and the ability of more than one mariner to occupy the same genome.

Mariners are a large family of eukaryotic DNA-mediated transposable elements that move via a cut-and-paste mechanism. Several features of the evolutionary history of mariners are unusual. First, they appear to undergo horizontal transfer commonly between species on an evolutionary timescale. They can do this because they are able to transpose using only their own self-encoded transposase and not host-specific factors. One consequence of this phenomenon is that more than one kind of mariner can be present in the same genome. We hypothesized that two mariners occupying the same genome would not interact. We tested the limits of mariner interactions using an in vitro transposition system, purified mariner transposases, and DNAse I footprinting. Only mariner elements that were very closely related to each other (ca. 84% identity) cross-mobilized, and then inefficiently. Because of the dramatic suppression of transposition between closely related elements, we propose that to isolate elements functionally, only minor changes might be necessary between elements, in both inverted terminal repeat and amino acid sequence. We further propose a mechanism to explain mariner diversification based on this phenomenon.

Updating the str and srj (stl) families of chemoreceptors in Caenorhabditis nematodes reveals frequent gene movement within and between chromosomes.

The seven transmembrane receptor (str) and srj (renamed from stl) families of chemoreceptors have been updated and the genes formally named following completion of the Caenorhabditis elegans genome sequencing project. Analysis of gene locations revealed that 84% of the 320 genes and pseudogenes in these two families reside on the large chromosome V. Movements to other chromosomes, especially chromosome IV, have nevertheless been relatively common, but only one has led to further gene family diversification. Comparisons with homologs in C. briggsae indicated that 22.5% of these genes have been newly formed by gene duplication since the species split, while also showing that four have been lost by large deletions. These patterns of gene evolution are similar to those revealed by analysis of the equally large srh family of chemoreceptors, and are likely to reflect general features of nematode genome dynamics. Thus large random deletions presumably balance the rapid proliferation of genes and their degeneration into pseudogenes, while gene movement within and between chromosomes keeps these nematode genomes in flux.

The tetraspanin superfamily in insects.

We describe four members of the tetraspanin/TM4SF superfamily of proteins that were identified in expressed sequence tag projects on the antennae of Manduca sexta moths and Apis mellifera honey bees. The three moth genes are expressed in the sensillar epithelium of male antennae, and some are expressed in female antennae, haemocytes, wing scale cell primordia and/or embryonic tissues. These proteins are probably involved in diverse cellular processes, much like their vertebrate homologues. A phylogenetic analysis of all known tetraspanins, including thirty-seven members of the superfamily revealed by the Drosophila melanogaster genome and twenty in the nematode Caenorhabditis elegans genome, reveals some phylum-specific gene amplification, in particular a contiguous array of eighteen genes in the D. melanogaster genome.

Phylogenomic analysis of the alpha proteasome gene family from early-diverging eukaryotes.

We employed a phylogenomic approach to study the evolution of alpha subunits of the proteasome gene family from early diverging eukaryotes. BLAST similarity searches of the Giardia lamblia genome identified all seven alpha proteasome genes characteristic of eukaryotes from the crown group. In addition, a PCR strategy for the amplification of multiple alpha subunit sequences generated single alpha proteasome products for representatives of the Kinetoplastida (Leishmania major), the Parabasalia (Trichomonas vaginalis), and the Microsporidia (Vairimorpha sp., Nosema sp., Endoreticulata sp., and Spraguea lophii). The kinetoplastid Trypanosoma cruzi and the eukaryote crown group Acanthamoeba castellanii yielded two distinct alpha proteasome genes each. The presence of seven distinct alpha proteasome genes in G. lamblia, one of the earliest-diverging eukaryotes, indicates that the alpha proteasome gene family evolved rapidly from a minimum of one gene in Archaea to seven or more in Eukarya. Results from the phylogenomic analysis are consistent with the idea that the Diplomonida (as represented by G. lamblia), the Kinetoplastida, the Parabasalia, and the Microsporidia diverged after the duplication events that originated the alpha proteasome gene family. A model for the early origin and evolution of the proteasome gene family is presented.

In vivo transposon mutagenesis of the methanogenic archaeon Methanosarcina acetivorans C2A using a modified version of the insect mariner-family transposable element Himar1.

We present here a method for in vivo transposon mutagenesis of a methanogenic archaeon, Methanosarcina acetivorans C2A, which because of its independence from host-specific factors may have broad application among many microorganisms. Because there are no known Methanosarcina transposons we modified the mariner transposable element Himar1, originally found in the insect Hematobia irritans, to allow its use in this organism. This element was chosen because, like other mariner elements, its transposition is independent of host factors, requiring only its cognate transposase. Modified mini-Himar1 elements were constructed that carry selectable markers that are functional in Methanosarcina species and that express the Himar1 transposase from known Methanosarcina promoters. These mini-mariner elements transpose at high frequency in M. acetivorans to random sites in the genome. The presence of an Escherichia coli selectable marker and plasmid origin of replication within the mini-mariner elements allows facile cloning of these transposon insertions to identify the mutated gene. In preliminary experiments, we have isolated numerous mini-mariner-induced M. acetivorans mutants, including ones with insertions that confer resistance to toxic analogs and in genes that encode proteins involved in heat shock, nitrogen fixation, and cell-wall structures.

Evolution of different subfamilies of mariner elements within the medfly genome inferred from abundance and chromosomal distribution.

The abundance and distribution pattern of eight mariner elements from three different subfamilies in the genome of the medfly Ceratitis capitata were determined. The copy numbers, as determined by slot-blot analysis, were very different for these elements. Their abundance did not change significantly within the native, the ancient or the newly derived populations, indicating that the rapid colonization process of the medfly had not affected the copy number of mariner elements. The distribution of the mariner elements was analyzed using fluorescent in situ hybridization (FISH) with charge-coupled device (CCD) camera analysis. The pattern of distribution in euchromatin and heterochromatin varied greatly and was distinctive and specific for each element. The implications of these findings are discussed and it is concluded that they generally support the hypothesis of a transposition/selection model in which the abundance and distribution patterns of these elements are regulated primarily by selection against deleterious effects due to meiotic ectopic recombination, while genetic drift would have played a minor role.

The large srh family of chemoreceptor genes in Caenorhabditis nematodes reveals processes of genome evolution involving large duplications and deletions and intron gains and losses.

The srh family of chemoreceptors in the nematode Caenorhabditis elegans is very large, containing 214 genes and 90 pseudogenes. It is related to the str, stl, and srd families of seven-transmembrane or serpentine receptors. Like these three families, most srh genes are concentrated on chromosome V, and mapping of their chromosomal locations on a phylogenetic tree reveals 27 different movements of genes to other chromosomes. Mapping of intron gains and losses onto the phylogenetic tree reveals that the last common ancestral gene of the family had five introns, which are inferred to have been lost 70 times independently during evolution of the family. In addition, seven intron gains are revealed, three of which are fairly recent. Comparisons with 20 family members in the C. briggsae genome confirms these patterns, including two intron losses in C. briggsae since the species split. There are 14 clear C. elegans orthologs for these 20 genes, whose average amino acid divergence of 68% allows estimation of 85 gene duplications in the C. elegans lineage since the species split. The absence of six orthologs in C. elegans also indicates that gene loss occurs; consideration of all deletions and terminal truncations of srh pseudogenes reveals that large deletions are common. Together these observations provide insight into the evolutionary dynamics of this compact animal genome.