Comparison of EST libraries from seven beetle species: towards a framework for phylogenomics of the Coleoptera

Validation of reference genes for quantitative real-time PCR in tiger beetles across sexes, body parts, sexual maturity and immune challenge

Reference genes are frequently used as normalizers for expression studies despite not being previously verified to present suitable stabilities. Considering the interest that tiger beetles have generated in the past years, resulting in a variety of studies, it is crucial to dispose of a validated reference gene panel for expression studies. Nine candidate genes were tested in Cicindela campestris and Calomera littoralis across several conditions and their transcription levels were assessed with geNorm, NormFinder, BestKeeper and ΔC_Tmethod algorithms. Results showed high stabilities across sexes, immune challenge and gonad developmental stages for all genes tested, while body parts comparison presented less constant expression values. Only two genes are sufficient to perform a proper normalization for most of the conditions tested, except for the body parts comparison in C. littoralis, which requires the use of at least three reference genes. On the whole, no universal gene is found to be suitable for all situations, but according to the acceptable range of values, NADH, B-t, Vatpase and ArgKin seem to present the most constant expression stability, indicating their suitability as reference genes in most of the conditions. This is the first report evaluating the stability of housekeeping genes in adephagan beetles.

Characterisation and expression analysis of UBC9 and UBS27 genes in developing gonads of cicindelids (Coleoptera: Cicindelidae)

Ubiquitin and small ubiquitin-like modifiers (SUMO) are post-translational modifiers essential in a variety of cellular processes, including gametogenesis. SUMO-conjugating enzyme (UBC9) and the ubiquitin ribosomal fusion protein UBS27 have been characterised in several model species. However, their expression in coleopteran remains unstudied. In this study, UBC9 and UBS27 genes have been characterised in the tiger beetle Cicindela campestris for the first time. Bioinformatic analysis showed that the Cc-UBC9 gene encoded a 159 amino acid protein with a predicted molecular weight of 18.18kDa, and the Cc-UBS27 gene encoded a 156 amino acid protein with a predicted molecular weight of 17.71kDa. Selection analyses carried out in several cicindelid species revealed that both genes were affected by purifying selection. Real time quantitative PCR analysis demonstrated that Cc-UBC9 and Cc-UBS27 were expressed in different tissues. The highest expression on both genes was found in the ovary and testis, and there were differential expression levels between immature and mature stages of testis development. The expression patterns of Cc-UBC9 and Cc-UBS27 suggest that these genes play important roles in gametogenesis in C. campestris. This information is relevant to better understand the reproductive process in cicindelids and the function of ubiquitin and small ubiquitin-related modifier genes in the Coleoptera.

Identification, structural characterisation and expression analysis of a defensin gene from the tiger beetle Calomera littoralis (Coleoptera: Cicindelidae)

In this study, a defensin gene (Clit-Def) has been characterised in the tiger beetle Calomera littoralis for the first time. Bioinformatic analysis showed that the gene has an open reading frame of 246bp that contains a 46 amino acid mature peptide. The phylogenetic analysis showed a high variability in the coleopteran defensins analysed. The Clit-Def mature peptide has the features to be involved in the antimicrobial function: a predicted cationic isoelectric point of 8.94, six cysteine residues that form three disulfide bonds, and the typical cysteine-stabilized α-helix β-sheet (CSαβ) structural fold. Real time quantitative PCR analysis showed that Clit-Def was upregulated in the different body parts analysed after infection with lipopolysaccharides of Escherichia coli, and also indicated that has an expression peak at 12h post infection. The expression patterns of Clit-Def suggest that this gene plays important roles in the humoral system in the adephagan beetle Calomera littoralis.

Insect phylogenomics

Phylogenomics, the integration of phylogenetics with genome data, has emerged as a powerful approach to study the evolution and systematics of species. Recently, several studies employing phylogenomic tools have provided better insights into insect evolution. Next-generation sequencing methods are now increasingly used by entomologists to generate genomic and transcript sequences of various insect species and strains. These data provide opportunities for comparative genomics and large-scale multigene phylogenies of diverse lineages of insects. Phy-logenomic investigations help us to better understand systematic and evolutionary relationships of insect species that play important roles as herbivores, predators, detritivores, pollinators and disease vectors. It is important that we critically assess the prospects and limitations of phylogenomic methods. In this review, I describe the current status, outline the major challenges and remark on potential future applications of phylogenomic tools in studying insect systematics and evolution.

De novo transcriptome assembly and SNP discovery in the wing polymorphic salt marsh beetle Pogonus chalceus (Coleoptera, Carabidae)

Background

The salt marsh beetle Pogonus chalceus represents a unique opportunity to understand and study the origin and evolution of dispersal polymorphisms as remarkable inter-population divergence in dispersal related traits (e.g. wing development, body size and metabolism) has been shown to persist in face of strong homogenizing gene flow. Sequencing and assembling the transcriptome of P. chalceus is a first step in developing large scale genetic information that will allow us to further study the recurrent phenotypic evolution in dispersal traits in these natural populations.

Methodology/Results

We used the Illumina HiSeq2000 to sequence 37 Gbases of the transcriptome and performed de novo transcriptome assembly with the Trinity short read assembler. This resulted in 65,766 contigs, clustering into 39,393 unique transcripts (unigenes). A subset of 12,987 show similarity (BLAST) to known proteins in the NCBI database and 7,589 are assigned Gene Ontology (GO). Using homology searches we identified all reported genes involved in wing development, juvenile- and ecdysteroid hormone pathways in Tribolium castaneum. About half (56.7%) of the unique assembled genes are shared among three life stages (third-instar larva, pupa, and imago). We identified 38,141 single nucleotide polymorphisms (SNPs) in these unigenes. Of these SNPs, 26,823 (70.3%) were found in a predicted open reading frame (ORF) and 6,998 (18.3%) were nonsynonymous.

Conclusions

The assembled transcriptome and SNP data are essential genomic resources for further study of the developmental pathways, genetic mechanisms and metabolic consequences of adaptive divergence in dispersal power in natural populations.

Mosquito genomics: progress and challenges

The whole-genome sequencing of mosquitoes has facilitated our understanding of fundamental biological processes at their basic molecular levels and holds potential for application to mosquito control and prevention of mosquito-borne disease transmission. Draft genome sequences are available for Anopheles gambiae, Aedes aegypti, and Culex quinquefasciatus. Collectively, these represent the major vectors of African malaria, dengue fever and yellow fever viruses, and lymphatic filariasis, respectively. Rapid advances in genome technologies have revealed detailed information on genome architecture as well as phenotype-specific transcriptomics and proteomics. These resources allow for detailed comparative analyses within and across populations as well as species. Next-generation sequencing technologies will likely promote a proliferation of genome sequences for additional mosquito species as well as for individual insects. Here we review the current status of genome research in mosquitoes and identify potential areas for further investigations.

Ribosomal protein genes of holometabolan insects reject the Halteria, instead revealing a close affinity of Strepsiptera with Coleoptera

The phylogenetic relationships among holometabolan insect orders remain poorly known, despite a wealth of previous studies. In particular, past attempts to clarify the sister-group of the enigmatic order Strepsiptera with rRNA genes have led to intense debate about long-branch attraction (the 'Strepsiptera problem'), without resolving the taxonomic question at hand. Here, we appealed to alternative nuclear sequences of 27 ribosomal proteins (RPs) to generate a data matrix of 10,731 nucleotides for 22 holometabolan taxa, including two strepsipteran species. Phylogenetic relationships among holometabolan insects were analyzed under several nucleotide-coding schemes to explore differences in signal and systematic biases. Saturation and compositional bias particularly affected third positions, which greatly differed in AT content (18-72%). Such confounding factors were best reduced by R-Y coding and removal of third codon positions, resulting in more strongly supported topologies, whereas amino acid coding gave poor resolution. The placement of Strepsiptera with Coleoptera (the Coleopterida) was recovered under most coding schemes and analytical methods, if often with modest support and ambiguity. In contrast, an alternative sister-group with Diptera (the Halteria) was only found in one analysis using parsimony, and weakly supported. The topologies here generally support a Coleoptera+Strepsiptera as sister-group to Mecopterida (Siphonaptera+Mecoptera+Diptera+Lepidoptera+Trichoptera), while Hymenoptera were always recovered as sister-group to the remaining Holometabola.

Development of bacterial artificial chromosome library resources for parasitoid Hymenoptera (Nasonia vitripennis and Nasonia giraulti: Pteromalidae)

The species of the genus Nasonia possess qualities that make them excellent candidates for genetic and genomic studies. To increase the wealth of genomic resources for the genus we constructed publicly available bacterial artificial chromosome (BAC) libraries for Nasonia vitripennis and Nasonia giraulti. Libraries have 36 864 clones each, empty-vector contents of approximately 2% and average insert sizes of 113.1 and 97.7 Kb, respectively, representing 12 and 11 genome equivalents. The N. vitripennis library was used for genome sequence assembly and in efforts at positional cloning of a developmental gene. The genome assembly of N. vitripennis is currently composed on 6181 un-joined scaffolds. These BAC libraries can be used to identify and close regions between scaffolds of the genome assemblies of both species.

Leveraging skewed transcript abundance by RNA-Seq to increase the genomic depth of the tree of life

Assembling the tree of life is a major goal of biology, but progress has been hindered by the difficulty and expense of obtaining the orthologous DNA required for accurate and fully resolved phylogenies. Next-generation DNA sequencing technologies promise to accelerate progress, but sequencing the genomes of hundreds of thousands of eukaryotic species remains impractical. Eukaryotic transcriptomes, which are smaller than genomes and biased toward highly expressed genes that tend to be conserved, could potentially provide a rich set of phylogenetic characters. We sampled the transcriptomes of 10 mosquito species by assembling 36-bp sequence reads into phylogenomic data matrices containing hundreds of thousands of orthologous nucleotides from hundreds of genes. Analysis of these data matrices yielded robust phylogenetic inferences, even with data matrices constructed from surprisingly few sequence reads. This approach is more efficient, data-rich, and economical than traditional PCR-based and EST-based methods and provides a scalable strategy for generating phylogenomic data matrices to infer the branches and twigs of the tree of life.

Phylogenetic molecular function annotation

It is now easier to discover thousands of protein sequences in a new microbial genome than it is to biochemically characterize the specific activity of a single protein of unknown function. The molecular functions of protein sequences have typically been predicted using homology-based computational methods, which rely on the principle that homologous proteins share a similar function. However, some protein families include groups of proteins with different molecular functions. A phylogenetic approach for predicting molecular function (sometimes called "phylogenomics") is an effective means to predict protein molecular function. These methods incorporate functional evidence from all members of a family that have functional characterizations using the evolutionary history of the protein family to make robust predictions for the uncharacterized proteins. However, they are often difficult to apply on a genome-wide scale because of the time-consuming step of reconstructing the phylogenies of each protein to be annotated. Our automated approach for function annotation using phylogeny, the SIFTER (Statistical Inference of Function Through Evolutionary Relationships) methodology, uses a statistical graphical model to compute the probabilities of molecular functions for unannotated proteins. Our benchmark tests showed that SIFTER provides accurate functional predictions on various protein families, outperforming other available methods.

Phylogenomic analysis reveals bees and wasps (Hymenoptera) at the base of the radiation of Holometabolous insects

Comparative studies require knowledge of the evolutionary relationships between taxa. However, neither morphological nor paleontological data have been able to unequivocally resolve the major groups of holometabolous insects so far. Here, we utilize emerging genome projects to assemble and analyze a data set of 185 nuclear genes, resulting in a fully resolved phylogeny of the major insect model species. Contrary to the most widely accepted phylogenetic hypothesis, bees and wasps (Hymenoptera) are basal to the other major holometabolous orders, beetles (Coleoptera), moths (Lepidoptera), and flies (Diptera). We validate our results by meticulous examination of potential confounding factors. Phylogenomic approaches are thus able to resolve long-standing questions about the phylogeny of insects.

Multiple nuclear gene sequences identify phylogenetic species boundaries in the rapidly radiating clade of Mexican ambystomatid salamanders

Delimiting the boundaries of species involved in radiations is critical to understanding the tempo and mode of lineage formation. Single locus gene trees may or may not reflect the underlying pattern of population divergence and lineage formation, yet they constitute the vast majority of the empirical data in species radiations. In this study we make use of an expressed sequence tag (EST) database to perform nuclear (nDNA) and mitochondrial (mtDNA) genealogical tests of species boundaries in Ambystoma ordinarium, a member of an adaptive radiation of metamorphic and paedomorphic salamanders (the Ambystoma tigrinum complex) that have diversified across terrestrial and aquatic environments. Gene tree comparisons demonstrate extensive nonmonophyly in the mtDNA genealogy of A. ordinarium, while seven of eight independent nuclear loci resolve the species as monophyletic or nearly so, and diagnose it as a well-resolved genealogical species. A differential introgression hypothesis is supported by the observation that western A. ordinarium localities contain mtDNA haplotypes that are identical or minimally diverged from haplotypes sampled from a nearby paedomorphic species, Ambystoma dumerilii, while most nDNA trees place these species in distant phylogenetic positions. These results provide a strong example of how historical introgression can lead to radical differences between gene trees and species histories, even among currently allopatric species with divergent life history adaptations and morphologies. They also demonstrate how EST-based nuclear resources can be used to more fully resolve the phylogenetic history of species radiations.

Phylogenetic and structural relationships of the PR5 gene family reveal an ancient multigene family conserved in plants and select animal taxa

Pathogenesis-related group 5 (PR5) plant proteins include thaumatin, osmotin, and related proteins, many of which have antimicrobial activity. The recent discovery of PR5-like (PR5-L) sequences in nematodes and insects raises questions about their evolutionary relationships. Using complete plant genome data and discovery of multiple insect PR5-L sequences, phylogenetic comparisons among plants and animals were performed. All PR5/PR5-L protein sequences were mined from genome data of a member of each of two main angiosperm groups-the eudicots (Arabidoposis thaliana) and the monocots (Oryza sativa)-and from the Caenorhabditis nematode (C. elegans and C. briggsase). Insect PR5-L sequences were mined from EST databases and GenBank submissions from four insect orders: Coleoptera (Diaprepes abbreviatus and Biphyllus lunatus), Orthoptera (Schistocerca gregaria), Hymenoptera (Lysiphlebus testaceipes), and Hemiptera (Toxoptera citricida). Parsimony and Bayesian phylogenetic analyses showed that the PR5 family is paraphyletic in plants, likely arising from 10 genes in a common ancestor to monocots and eudicots. After evolutionary divergence of monocots and eudicots, PR5 genes increased asymmetrically among the 10 clades. Insects and nematodes contain multiple sequences (seven PR5-Ls in nematodes and at least three in some insects) all related to the same plant clade, with nematode and insect sequences separating as two clades. Protein structural homology modeling showed strong similarity among animal and plant PR5/PR5-Ls, with divergence only in surface-exposed loops. Sequence and structural conservation among PR5/PR5-Ls suggests an important and conserved role throughout the evolutionary divergence of the diverse organisms from which they reside.

Caste-specific cytochrome P450 in the damp-wood termite Hodotermopsis sjostedti (Isoptera, Termopsidae)

Termites are eusocial insects with a well-defined caste system, which is an example of polyphenism. This polyphenism is based on hormonally controlled differential gene expression. In the damp-wood termite Hodotermopsis sjostedti, we induced differentiation into the soldier caste by using juvenile hormone analogue treatment. We then investigated specific gene expression, which appeared during the hormonal response and triggered caste differentiation, using fluorescent differential display. A candidate cDNA sequence with similarity to cytochromes P450, CYP6AM1, was characterized and its transcript shown to be repressed between 1 and 3 days after hormone treatment. CYP6AM1 was specifically expressed in the fat body of pseudergates and soldiers. The putative function of this P450 is discussed with respect to the caste differentiation system.

OrthologID: automation of genome-scale ortholog identification within a parsimony framework

MOTIVATION

The determination of gene orthology is a prerequisite for mining and utilizing the rapidly increasing amount of sequence data for genome-scale phylogenetics and comparative genomic studies. Until now, most researchers use pairwise distance comparisons algorithms, such as BLAST, COG, RBH, RSD and INPARANOID, to determine gene orthology. In contrast, orthology determination within a character-based phylogenetic framework has not been utilized on a genomic scale owing to the lack of efficiency and automation.

RESULTS

We have developed OrthologID, a Web application that automates the labor-intensive procedures of gene orthology determination within a character-based phylogenetic framework, thus making character-based orthology determination on a genomic scale possible. In addition to generating gene family trees and determining orthologous gene sets for complete genomes, OrthologID can also identify diagnostic characters that define each orthologous gene set, as well as diagnostic characters that are responsible for classifying query sequences from other genomes into specific orthology groups. The OrthologID database currently includes several complete plant genomes, including Arabidopsis thaliana, Oryza sativa, Populus trichocarpa, as well as a unicellular outgroup, Chlamydomonas reinhardtii. To improve the general utility of OrthologID beyond plant species, we plan to expand our sequence database to include the fully sequenced genomes of prokaryotes and other non-plant eukaryotes.

AVAILABILITY

http://nypg.bio.nyu.edu/orthologid/

Dense taxonomic EST sampling and its applications for molecular systematics of the Coleoptera (beetles)

Expressed sequence tag (EST) sequences can provide a wealth of data for phylogenetic and genomic studies, but the utility of these resources is restricted by poor taxonomic sampling. Here, we use small EST libraries (<1,000 clones) to generate phylogenetic markers across a broad sample of insects, focusing on the species-rich Coleoptera (beetles). We sequenced over 23,000 ESTs from 34 taxa, which produced 8,728 unique sequences after clustering nonredundant sequences. Between taxa, the sequences could be grouped into 731 gene clusters, with the largest corresponding to mitochondrial DNA transcripts and gene families chymotrypsin, actin, troponin, and tubulin. While levels of paralogy were high in most gene clusters, several midsized clusters including many ribosomal protein (RP) genes appeared to be free of expressed paralogs. To evaluate the utility of EST data for molecular systematics, we curated available transcripts for 66 RP genes from representatives of the major groups of Coleoptera. Using supertree and supermatrix approaches for phylogenetic analysis, the results were consistent with the emerging phylogenetic conclusions about basal relationships in Coleoptera. Numerous small EST libraries from a taxonomically densely sampled lineage can provide a core set of genes that together act as a scaffold in phylogenetic reconstruction, comparative genomics, and studies of gene evolution.

Animal phylogenomics: multiple interspecific genome comparisons

The utility of DNA sequence information for phylogenetics and phylogeography is now well known. Rather than attempt to summarize studies addressing this well-demonstrated utility, this chapter focuses on fundamental approaches and techniques that implement the collection of DNA sequence data for comparative phylogenetic purposes in a genomic context (phylogenomics). Whole genome sequencing approaches have changed the way we think about phylogenetics and have opened the way for new perspectives on "old" phylogenetics concerns. Some of these concerns are which gene regions to use and how much sequence information is needed for robust phylogenetic inference. Whole genome sequences of a few animal model organisms have gone a long way to implement approaches to better understand these important phylogenetic concerns. This chapter also addresses how genomics has made it more important for a clear understanding of orthology of gene regions in comparative biology. Finally, genome-enabled technologies that are affecting comparative biology are also discussed.

Using exon and intron sequences of the gene Mp20 to resolve basal relationships in Cicindela (Coleoptera:Cicindelidae)

The genus Cicindela (Coleoptera: Cicindelidae) is a species-rich cosmopolitan group of tiger beetles useful for comparing clade diversification worldwide. Knowledge about relationships of major groups is important for this analysis but basal nodes in Cicindela have been difficult to resolve with standard mtDNA markers. Here we developed the Mp20 gene, a single-copy nuclear marker coding for a muscle-associated protein in insects, for phylogenetic analysis of basal groups of Cicindela. Nearly full-length sequences were obtained for 51 cicindelids, including major taxonomic groups from all continents. Sequences of Mp20 were between 1.2 and 1.7 kb and spanning three introns. Phylogenetic signal of exon and intron sequences was compared with that from four gene regions of mtDNA (COI, COIII, Cytb, 16S rRNA; 2.4 kb total). Because introns differed in length, sequence alignment was conducted using various procedures of phenetic and parsimony-based character coding of indels to assess their phylogenetic information content, but major nodes were recovered consistently. Mp20 sequences contributed two thirds of the total support of the combined analysis, with most signal from the introns. We found major clades of Cicindela to be geographically largely coincident with continental regions, confined to Australasia, the Holarctic, the Indian subcontinent, Africa, and South and Central America. Clock estimates using various maximum-likelihood (ML) branch length calculations resulted in roughly similar divergence times whether Mp20 exon, introns, or mtDNA were used, and they were not greatly affected by different procedures for coding and optimizing indel characters. Based on existing clock calibrations in Cicindela, basal splits of continental lineages occurred in the mid-Miocene, placing the radiation of basal groups of Cicindela to a period when their open-vegetation habitats expanded globally.

The phylogeny of acorn weevils (genus Curculio) from mitochondrial and nuclear DNA sequences: the problem of incomplete data

We considered the contribution of two mitochondrial and two nuclear data sets for the phylogenetic reconstruction of 22 species of seed beetles in the genus Curculio (Coleoptera: Cuculionidae). A phylogenetic tree from representatives found on various hosts was inferred from a combined data set of mitochondrial DNA cytochrome oxidase subunit I, mitochondrial cytochrome b, nuclear elongation factor 1alpha, and nuclear phosphoglycerate mutase, used for the first time as a molecular marker. Separate parsimony analyses of each data set showed that individual gene trees were mainly congruent and often complementary in the support of clades but the analysis was complicated by failure of PCR amplification of nuclear genes for many taxa and hence missing data entries. When the four gene partitions were combined in a simultaneous analysis despite the missing data, this increased the resolution and taxonomic coverage compared to the individual source trees. Alternative approaches of combining the information via supertree methodology produced a comparatively less resolved tree, and hence seem inferior to combining data matrices even in cases where numerous taxa are missing. The molecular data suggest a classification of the European species into two species groups that are in accordance with morphological characteristics but the data do no support any of the previously recognised American species groups.

Comparison of gene representation in midguts from two phytophagous insects, Bombyx mori and Ips pini, using expressed sequence tags

Midgut proteins may provide new molecular targets for insect control. This could be particularly important for some pests, such as pine bark beetles, which are difficult to control by conventional methods. Expressed sequence tags (ESTs) provide information about the activity of a particular tissue, and, in the case of pest insects, may quickly identify potential targets. We present here an EST project representing 574 tentative unique genes (TUGs) expressed in the midgut of the male pine engraver beetle, Ips pini. This tissue uses the mevalonate pathway to produce the monoterpenoid pheromone component, ipsdienol, de novo in response to juvenile hormone (JH) III. Comparison of our ESTs with those previously isolated from larval silkmoth (Bombyx mori) midguts revealed interesting similarities and differences in gene representation that correlate with the conserved and divergent functions of these two tissues. For example, seven mevalonate pathway genes were represented in the I. pini ESTs, while none were found from B. mori. This type of comparison may assist the identification of species-specific targets for future pest control strategies.

The construction of an EST database for Bombyx mori and its application

To build a foundation for the complete genome analysis of Bombyx mori, we have constructed an EST database. Because gene expression patterns deeply depend on tissues as well as developmental stages, we analyzed many cDNA libraries prepared from various tissues and different developmental stages to cover the entire set of Bombyx genes. So far, the Bombyx EST database contains 35,000 ESTs from 36 cDNA libraries, which are grouped into approximately 11,000 nonredundant ESTs with the average length of 1.25 kb. The comparison with FlyBase suggests that the present EST database, SilkBase, covers >55% of all genes of Bombyx. The fraction of library-specific ESTs in each cDNA library indicates that we have not yet reached saturation, showing the validity of our strategy for constructing an EST database to cover all genes. To tackle the coming saturation problem, we have checked two methods, subtraction and normalization, to increase coverage and decrease the number of housekeeping genes, resulting in a 5-11% increase of library-specific ESTs. The identification of a number of genes and comprehensive cloning of gene families have already emerged from the SilkBase search. Direct links of SilkBase with FlyBase and WormBase provide ready identification of candidate Lepidoptera-specific genes.

Transcriptome identification of putative genes involved in protein catabolism and innate immune response in human body louse (Pediculicidae: Pediculus humanus)

Genomics information relating to human body lice is surprisingly scarce, and this has constrained studies of their physiology, immunology and vector biology. To identify novel body louse genes, we used engorged adult lice to generate a cDNA library. Initially, 1152 clones were screened for inserts, edited for removal of vector sequences and base pairs of poor quality, and viewed for splicing variations, gene families and polymorphism. Computational methods identified 506 inferred open reading frames including the first predicted louse defensin. The inferred defensin aligns well with other insect defensins and has highly conserved cysteine residues, as are known for other defensin sequences. Two cysteine and five serine proteinases were categorized according to their inferred catalytic sites. We also discovered seven putative ubiquitin-pathway genes and four iron metabolizing deduced enzymes. Finally, glutathione-S-transferases and cytochrome P450 genes were among the detoxification enzymes found. Results from this first systematic effort to discover human body louse genes should promote further studies in Phthiraptera and lice.

Genomics, evolution and biological functions of the pacifastin peptide family: a conserved serine protease inhibitor family in arthropods

The last decade, a new serine protease inhibitor family has been described in arthropods. Eight members were purified from the locusts Locusta migratoria (LMPI-1-2 and HI) and Schistocerca gregaria (SGPI-1-5). The light chain of the heterodimeric protease inhibitor pacifastin, from the freshwater crayfish Pacifastacus leniusculus, was found to be composed of nine consecutive inhibitory domains (PLDs). These domains share a pattern of six conserved cysteine residues (Cys-Xaa(9-12)-Cys-Asn-Xaa-Cys-Xaa-Cys-Xaa(2-3)-Gly-Xaa(3-6)-Cys-Thr-Xaa(3)-Cys) with the locust inhibitors. Via cDNA cloning, eight pacifastin-related precursors have been identified in locusts. Interestingly, additional pacifastin-related precursors have been identified in Diptera, Lepidoptera and Coleoptera utilising an in silico data mining approach.

Transcriptome analysis of the cowpea weevil bruchid: identification of putative proteinases and alpha-amylases associated with food breakdown

We describe here the first systematic work to discover insect genes involved in food breakdown using a cDNA library enriched for gut-expressed transcripts from Callosobruchus maculatus. A total of 1056 clones were screened for cDNA insert-containing plasmids, and 503 nonredundant open reading frames were discovered. Twenty-three inferred genes potentially involved in digestive processes in cowpea weevil were identified, including proteinases and amylases. The predicted catalytic sites were identified in the inferred cysteine and aspartic acid proteinases, and in alpha-amylases. Transcriptome analysis of the cowpea bruchid will potentially permit gene discovery in other beetles, an insect order of major economic and ecological importance that is poorly represented in genomic databases.

Expressed sequence tags: alternative or complement to whole genome sequences?

Over three million sequences from approximately 200 plant species have been deposited in the publicly available plant expressed sequence tag (EST) sequence databases. Many of the ESTs have been sequenced as an alternative to complete genome sequencing or as a substrate for cDNA array-based expression analyses. This creates a formidable resource from both biodiversity and gene-discovery standpoints. Bioinformatics-based sequence analysis tools have extended the scope of EST analysis into the fields of proteomics, marker development and genome annotation. Although EST collections are certainly no substitute for a whole genome scaffold, this "poor man's genome" resource forms the core foundations for various genome-scale experiments within the as yet unsequenceable plant genomes.

Genomics in pure and applied entomology

Genomics is the study of the structure and function of the genome: the set of genetic information encoded in the DNA of the nucleus and organelles of an organism. It is a dynamic field that combines traditional paths of inquiry with new approaches that would have been impossible without recent technological developments. Much of the recent focus has been on obtaining the sequence of entire genomes, determining the order and organization of the genes, and developing libraries that provide immediate physical access to any desired DNA fragment. This has enabled functional studies on a genome-wide level, including analysis of the genetic basis of complex traits, quantification of global patterns of gene expression, and systematic gene disruption projects. The successful contribution of genomics to problems in applied entomology requires the cooperation of the private and public sectors to build upon the knowledge derived from the Drosophila genome and effectively develop models for other insect Orders.