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.

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.

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.

Interspecific and intraspecific horizontal transfer of Wolbachia in Drosophila.

Cytoplasmic incompatibility (CI) in Drosophila simulans is related to infection of the germ line by a rickettsial endosymbiont (genus Wolbachia). Wolbachia were transferred by microinjection of egg cytoplasm into uninfected eggs of both D. simulans and D. melanogaster to generate infected populations. Transinfected strains of D. melanogaster with lower densities of Wolbachia than the naturally infected D. simulans strain did not express high levels of CI. However, transinfected D. melanogaster egg cytoplasm, transferred back into D. simulans, generated infected populations that expressed CI at levels near those of the naturally infected strain. A transinfected D. melanogaster line selected for increased levels of CI expression also displayed increased symbiont densities. These data suggest that a threshold level of infection is required for normal expression of CI and that host factors help determine the density of the symbiont in the host.

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.

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.

Modified P elements that mimic the P cytotype in Drosophila melanogaster.

Activity of the P family of transposable elements in Drosophila melanogaster is regulated primarily by a cellular condition known as P cytotype. It has been hypothesized that P cytotype depends on a P element-encoded repressor of transposition and excision. We provide evidence in support of this idea by showing that two modified P elements, each with lesions affecting the fourth transposase exon, mimic most of the P cytotype effects. These elements were identified by means of two sensitive assays capable of detecting repression by a single P element. One assay makes use of cytotype-dependent gene expression of certain P element insertion mutations at the singed bristle locus. The other measures suppression of transposase activity from the unusually stable genomic P element, delta 2-3(99B), that normally produces transposase in both germinal and somatic tissues. The P cytotype-like effects include suppression of snw germline hypermutability, snw somatic mosaicism, pupal lethality, and gonadal dysgenic sterility. Unlike P cytotype, however, there was no reciprocal cross effect in the inheritance of repression.

Wolbachia infections and the expression of cytoplasmic incompatibility in Drosophila sechellia and D. mauritiana.

Various stocks of Drosophila mauritiana and D. sechellia were found to be infected with Wolbachia, a Rickettsia-like bacterium that is known to cause cytoplasmic incompatibility and other reproductive abnormalities in arthropods. Testing for the expression of cytoplasmic incompatibility in these two species showed partial incompatibility in D. sechellia but no expression of incompatibility in D. mauritiana. To determine whether absence of cytoplasmic incompatibility in D. mauritiana was due to either the bacterial or host genome, we transferred bacteria from D. mauritiana into an uninfected strain of D. simulans, a host species known to express high levels of incompatibility with endogenous Wolbachia. We also performed the reciprocal transfer of the natural D. simulans Riverside infection into a tetracycline-treated stock of D. mauritiana. In each case, the ability to express incompatibility was unaltered by the different host genetic background. These experiments indicate that in D. simulans and D. mauritiana expression of the cytoplasmic incompatibility phenotype is determined by the bacterial strain and that D. mauritiana harbors a neutral strain of Wolbachia.

Factors affecting transposition of the Himar1 mariner transposon in vitro.

Mariner family transposable elements are widespread in animals, but their regulation is poorly understood, partly because only two are known to be functional. These are particular copies of the Dmmar1 element from Drosophila mauritiana, for example, Mos1, and the consensus sequence of the Himar1 element from the horn fly, Haematobia irritans. An in vitro transposition system was refined to investigate several parameters that influence the transposition of Himar1. Transposition products accumulated linearly over a period of 6 hr. Transposition frequency increased with temperature and was dependent on Mg2+ concentration. Transposition frequency peaked over a narrow range of transposase concentration. The decline at higher concentrations, a phenomenon observed in vivo with Mos1, supports the suggestion that mariners may be regulated in part by "overproduction inhibition." Transposition frequency decreased exponentially with increasing transposon size and was affected by the sequence of the flanking DNA of the donor site. A noticeable bias in target site usage suggests a preference for insertion into bent or bendable DNA sequences rather than any specific nucleotide sequences beyond the TA target site.

A stable genomic source of P element transposase in Drosophila melanogaster.

A single P element insert in Drosophila melanogaster, called P[ry+ delta 2-3](99B), is described that caused mobilization of other elements at unusually high frequencies, yet is itself remarkably stable. Its transposase activity is higher than that of an entire P strain, but it rarely undergoes internal deletion, excision or transposition. This element was constructed by F. Laski, D. Rio and G. Rubin for other purposes, but we have found it to be useful for experiments involving P elements. We demonstrate that together with a chromosome bearing numerous nonautonomous elements it can be used for P element mutagenesis. It can also substitute efficiently for "helper" plasmids in P element mediated transformation, and can be used to move transformed elements around the genome.

Type I repressors of P element mobility.

We describe here a family of P elements that we refer to as type I repressors. These elements are identified by their repressor functions and their lack of any deletion within the first two-thirds of the canonical P sequence. Elements belonging to this repressor class were isolated from P strains and were made in vitro. We found that type I repressor elements could strongly repress both a cytotype-dependent allele and P element mobility in somatic and germline tissues. These effects were very dependent on genomic position. Moreover, we observed that an element's ability to repress in one assay positively correlated with its ability to repress in either of the other two assays. The type I family of repressor elements includes both autonomous P elements and those lacking exon 3 of the P element. Fine structure deletion mapping showed that the minimal 3' boundary of a functional type I element lies between nucleotide position 1950 and 1956. None of 12 elements examined with more extreme deletions extending into exon 2 made repressor. We conclude that the type I repressors form a structurally distinct group that does not include more extensively deleted repressor elements such as the KP element described previously.

Somatic effects of P element activity in Drosophila melanogaster: pupal lethality.

Nonautonomous P elements normally excise and transpose only when a source of transposase is supplied, and only in the germline. The germline specificity depends on one of the introns of the transposase gene which is not spliced in somatic cells. To study the effects of somatic P activity, a modified P element (delta 2-3) lacking this intron was used as a source of transposase. Nonautonomous P elements from a strain called Birmingham, when mobilized in somatic cells by delta 2-3, were found to cause lethality, although neither component was lethal by itself. The three major Birmingham chromosomes acted approximately independently in producing the lethal effect. This lethality showed a strong dependence on temperature. Although temperature sensitivity was limited to larval stages, the actual deaths occurred at the pupal stage. Survivors, which could be recovered by decreasing the temperature or by reducing the proportion of the Birmingham genome present, often showed multiple developmental anomalies and reduced longevity reminiscent of the effects of cell death from radiation damage. Although the genetic damage occurred in dividing imaginal disc cells, the phenotypic manifestations--death and abnormalities--are not observed until later. The survivors also showed gonadal dysgenic (GD) sterility, a well-known characteristic of P-M hybrid dysgenesis. To explain these findings, we suggest that pupal lethality and GD sterility are both caused by massive chromosome breakage in larval cells, resulting from excision and transposition of genomic P elements acting as substrate for the transposase.

Molecular evolution of the second ancient human mariner transposon, Hsmar2, illustrates patterns of neutral evolution in the human genome lineage.

A consensus sequence for the second ancient mariner identified in the human genome, Hsmar2, was constructed by majority rule from full-length and partial sequences of 44 of the +/-1000 copies in the genome. This 1300 base pair (bp) consensus has 31 bp imperfect terminal repeats (ITRs) and encodes a 351 amino acid (aa) mariner transposase. The sequence of this transposase has allowed classification of Hsmar2 as a basal lineage of the irritans subfamily of mariners, sharing at most 38% aa identity with other members of the subfamily. The individual copies in the human genome are all highly mutated from the consensus, having suffered numerous small and some large insertions and deletions (indels), including many insertions of S and J subfamily Alu elements. The copies differ, on average, from the consensus by 11.6%, have suffered 11.8 indels per kilobase (kb), and only 3.7% of the 30 hypermutable CpG dinucleotide pairs in the consensus remain intact. This level of divergence indicates that the ancestrally active Hsmar2 element represented by the consensus was present in the human genome lineage about 80 million years (Myr) ago. Each copy has apparently evolved since then largely independently of the others, and with little constraint on its transposase coding capacity. This pattern of molecular evolution fits the current model for mariner transposon evolution. These copies provide multiple independent datasets for evaluating the pattern of neutral evolution in the human genome, for example, they confirm that most indels are very short and that deletions are twice as common as insertions.

Molecular evolution of an ancient mariner transposon, Hsmar1, in the human genome.

A confident consensus sequence for Hsmar1, the first mariner transposon recognized in the human genome, was generated using three genomic and 15 cDNA sequences. It is thought to represent the ancestrally active copy that invaded an early primate genome. The consensus is 1287 base pairs (bp) long, has 30 bp perfect inverted terminal repeats (ITRs), and encodes a 343 amino acid (aa) mariner transposase. Each copy has diverged from the consensus largely independently of the others and mostly neutrally, and most are now defective. They differ from the consensus by an average of 7.8% in DNA sequence and 7.5 indels per kilobase, both of which values indicate that the copies were formed about 50 Myr ago. On average, only 20% of the 73 surmised CpG hypermutable sites in the consensus remain. A remarkable exception to this loss of functionality is revealed by a set of ten cDNA clones derived from a particular genomic copy that has diverged only 2.4% from the consensus, retained 54% of its hypermutable CpG pairs, and which has a full-length transposase open reading frame. The complete sequence of one of these cDNAs (NIB1543) indicates that the transposase gene of this copy may have been conserved because it is spliced to a human cellular gene encoding a SET domain protein. A specific PCR assay was used to reveal the presence of Hsmar1 copies in all primates examined representing all major lineages, but not in close relatives of primates. PCR fragments cloned and sequenced from a representative sample of primates confirmed that Hsmar1 copies are present in all major lineages, and also revealed another cecropia subfamily mariner in prosimians only, and a third highly divergent mariner present in the greater slow loris Nycticebus coucang. There are about 200 copies of Hsmar1 in the human genome, as well as +/-2400 copies of a derived 80 bp paired ITR structure and +/-4600 copies of solo ITRs. Thus, this transposon had a considerable insertional mutagenic effect on past primate genomes.

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.

Bmmar1: a basal lineage of the mariner family of transposable elements in the silkworm moth, Bombyx mori.

We describe a transposable element, called Bmmar1, from the genome of the silkworm moth, Bombyx mori. This element has features of the Tc1-mariner superfamily of transposable elements. Bmmar1 was first detected as a fragment in the 5' region of the larval serum protein (BmLSP) gene. Six genomic clones characterized each differed from a consensus sequence by 3-5 insertions and deletions, as well as an average of 2.3% in nucleotide sequence. The genome contains approximately 2400 copies of Bmmar1. Maximum parsimony phylogenetic analysis of the relationship of Bmmar1 and other members of the Tc1-mariner superfamily, based on their encoded transposase amino acid sequences, indicates that it represents a basal lineage of the mariner family. In particular Bmmar1 encodes a D,D37D motif thought to be the catalytic domain of mariner transposases. Bmmar1 considerably increases the known diversity of this widespread family of transposons. A new naming system is proposed for members of the family.

Expression of lacunin, a large multidomain extracellular matrix protein, accompanies morphogenesis of epithelial monolayers in Manduca sexta.

Morphogenesis is a complex process operating at several levels of organization--organism, tissues, cells, and molecules. Complex interactions occur between and within these levels. Many of the molecules that mediate these interactions are predictably turning out to be large multidomain proteins. Here we describe one such novel protein associated with remodeling of epithelial monolayers in embryos and developing wings of the moth Manduca sexta. On the basis of its sequence and its expression pattern along lacunae of developing wings, we propose the name lacunin for this extracellular matrix protein that contains nine different types of domains, most of which are present in multiple copies. These include domains of various types: Kunitz proteinase inhibitors, thrombospondin type I, immunoglobulin-like, and several newly defined domains of unknown function (PAL, PLAC, and lagrin domains). This rich patchwork of distinct domains probably exerts multiple effects on a variety of cell behaviors associated with the complex phenomenon of epithelial morphogenesis.

Inherited copper toxicosis in Bedlington terriers.

Chronic hepatitis and increased hepatic copper concentrations, from 1,600 to 6,361 micrograms/g dry tissue were found in 4 related, Australian-bred Bedlington terriers. Two dogs were asymptomatic and 2 were clinically ill with signs referable to liver dysfunction. Two dogs were treated with d-penicillamine. After one year there was no improvement in the histopathological liver changes in either dog or significant lowering of hepatic copper level in one dog.