Genomic tools and cDNA derived markers for butterflies

A phylogenomics approach to characterizing sensory neuron membrane proteins (SNMPs) in Lepidoptera

Sensory neuron membrane proteins (SNMPs) play a critical role in the insect olfactory system but there is a deficit of functional studies beyond Drosophila. Here, we use a combination of available genome sequences, manual curation, genome and transcriptome data, phylogenetics, expression profiling and gene knockdown to investigate SNMP superfamily in various insect species with a focus on Lepidoptera. We curated 81 genes from 36 insect species and identified a novel lepidopteran SNMP gene family, SNMP3. Phylogenetic analysis shows that lepidopteran SNMP3, but not the previously annotated lepidopteran SNMP2, is the true homologue of the dipteran SNMP2. Digital expression, microarray and qPCR analyses show that the lepidopteran SNMP1 is specifically expressed in adult antennae. SNMP2 is widely expressed in multiple tissues while SNMP3 is specifically expressed in the larval midgut. Microarray analysis suggest SNMP3 may be involved in the silkworm immunity response to virus and bacterial infections. We functionally characterized SNMP1 in the silkworm using RNA interference (RNAi) and behavioral assays. Our results suggested that Bombyx mori SNMP1 is a functional orthologue of the Drosophila melanogaster SNMP1 and plays a critical role in pheromone detection. Split-ubiquitin yeast hybridization study shows that BmorSNMP1 has a protein-protein interaction with the pheromone receptor (BmorOR1), and the co-receptor (BmorOrco). Concluding, we propose a novel molecular model in which BmorOrco, BmorSNMP1 and BmorOR1 form a heteromer in the detection of the silkworm sex pheromone bombykol.

The functional basis of wing patterning in Heliconius butterflies: the molecules behind mimicry

Wing-pattern mimicry in butterflies has provided an important example of adaptation since Charles Darwin and Alfred Russell Wallace proposed evolution by natural selection >150 years ago. The neotropical butterfly genus Heliconius played a central role in the development of mimicry theory and has since been studied extensively in the context of ecology and population biology, behavior, and mimicry genetics. Heliconius species are notable for their diverse color patterns, and previous crossing experiments revealed that much of this variation is controlled by a small number of large-effect, Mendelian switch loci. Recent comparative analyses have shown that the same switch loci control wing-pattern diversity throughout the genus, and a number of these have now been positionally cloned. Using a combination of comparative genetic mapping, association tests, and gene expression analyses, variation in red wing patterning throughout Heliconius has been traced back to the action of the transcription factor optix. Similarly, the signaling ligand WntA has been shown to control variation in melanin patterning across Heliconius and other butterflies. Our understanding of the molecular basis of Heliconius mimicry is now providing important insights into a variety of additional evolutionary phenomena, including the origin of supergenes, the interplay between constraint and evolvability, the genetic basis of convergence, the potential for introgression to facilitate adaptation, the mechanisms of hybrid speciation in animals, and the process of ecological speciation.

Comprehensive microarray-based analysis for stage-specific larval camouflage pattern-associated genes in the swallowtail butterfly, Papilio xuthus

Background

Body coloration is an ecologically important trait that is often involved in prey-predator interactions through mimicry and crypsis. Although this subject has attracted the interest of biologists and the general public, our scientific knowledge on the subject remains fragmentary. In the caterpillar of the swallowtail butterfly Papilio xuthus, spectacular changes in the color pattern are observed; the insect mimics bird droppings (mimetic pattern) as a young larva, and switches to a green camouflage coloration (cryptic pattern) in the final instar. Despite the wide variety and significance of larval color patterns, few studies have been conducted at a molecular level compared with the number of studies on adult butterfly wing patterns.

Results

To obtain a catalog of genes involved in larval mimetic and cryptic pattern formation, we constructed expressed sequence tag (EST) libraries of larval epidermis for P. xuthus, and P. polytes that contained 20,736 and 5,376 clones, respectively, representing one of the largest collections available in butterflies. A comparison with silkworm epidermal EST information revealed the high expression of putative blue and yellow pigment-binding proteins in Papilio species. We also designed a microarray from the EST dataset information, analyzed more than five stages each for six markings, and confirmed spatial expression patterns by whole-mount in situ hybridization. Hence, we succeeded in elucidating many novel marking-specific genes for mimetic and cryptic pattern formation, including pigment-binding protein genes, the melanin-associated gene yellow-h3, the ecdysteroid synthesis enzyme gene 3-dehydroecdysone 3b-reductase, and Papilio-specific genes. We also found many cuticular protein genes with marking specificity that may be associated with the unique surface nanostructure of the markings. Furthermore, we identified two transcription factors, spalt and ecdysteroid signal-related E75, as genes expressed in larval eyespot markings. This finding suggests that E75 is a strong candidate mediator of the hormone-dependent coordination of larval pattern formation.

Conclusions

This study is one of the most comprehensive molecular analyses of complicated morphological features, and it will serve as a new resource for studying insect mimetic and cryptic pattern formation in general. The wide variety of marking-associated genes (both regulatory and structural genes) identified by our screening indicates that a similar strategy will be effective for understanding other complex traits.

Developing EPIC markers for chalcidoid Hymenoptera from EST and genomic data

Increasing numbers of phylogeographic studies make comparative inferences about the histories of co-distributed species. Although the aims of such studies are best achieved by jointly analysing sequences from multiple loci in a model-based framework, such data currently exist for few nonmodel systems. We used existing genomic data and expressed sequence tags (ESTs) for Hymenoptera and other insects to design intron-crossing primers for 40 loci, mainly ribosomal proteins, for chalcidoid parasitoids. Amplification success was scored on a range of taxa associated with two natural communities; oak galls and figs. Taxa were chosen at increasing distance from Nasonia, which was used for primer design, (i) within Pteromalids, (ii) within Chalcidoidea (Eupelmidae, Eulophidae, Eurytomidae, Ormyridae, Torymidae) and (iii) for a selection of distantly related gall and fig wasps (Cynipidae, Agaonidae). To assess the utility of these loci for phylogeographic and population genetic studies, we compared genetic diversity between Western Palaearctic refugia for two species. Our results show that it is feasible to design a large number of exon-primed-intron-crossing (EPIC) loci that may be informative about phylogeographic history within species but amplify across a large taxonomic range.

Isolation, annotation and applications of expressed sequence tags from the olive fly, Bactrocera oleae

The olive fruit fly, Bactrocera oleae, is the major pest of the olive tree. Despite its importance, very little genetic and molecular knowledge is available. The present study is a first attempt to identify and characterize B. oleae expressed sequence tags (ESTs). One hundred and ninety-five randomly selected cDNA clones were isolated and the obtained sequences were annotated through BLASTX similarity searches. A set of 159 unique putative transcripts were functionally assigned using Gene Ontology terms in broad categories of biological process, molecular function and cellular component based on D. melanogaster matches. Moreover, the cytogenetic location of 35 ESTs was determined by in situ hybridization to B. oleae polytene chromosomes. The resulting low-resolution EST map more than doubles the available entry points to the insect's genome and can assist syntenic comparisons with other distant species. The deduced codon usage of the isolated ESTs suggested a conserved pattern of B. oleae with its closest relatives. Additionally, the comparative analysis of B. oleae ESTs with the homologous D. melanogaster genes led to the development of 17 nuclear EPIC-PCR markers for the amplification of intron sequences of 11 Tephritidae species. Sequencing analysis of several cross-amplified intron sequences revealed a high degree of conservation among Bactrocera species and a varying transferability of the generated markers across the examined genera, suggesting that this method can provide a useful tool for the clarification of phylogenetic relationships among different species, particularly in cases of species complexes.

Population-level transcriptome sequencing of nonmodel organisms Erynnis propertius and Papilio zelicaon

BACKGROUND

Several recent studies have demonstrated the use of Roche 454 sequencing technology for de novo transcriptome analysis. Low error rates and high coverage also allow for effective SNP discovery and genetic diversity estimates. However, genetically diverse datasets, such as those sourced from natural populations, pose challenges for assembly programs and subsequent analysis. Further, estimating the effectiveness of transcript discovery using Roche 454 transcriptome data is still a difficult task.

RESULTS

Using the Roche 454 FLX Titanium platform, we sequenced and assembled larval transcriptomes for two butterfly species: the Propertius duskywing, Erynnis propertius (Lepidoptera: Hesperiidae) and the Anise swallowtail, Papilio zelicaon (Lepidoptera: Papilionidae). The Expressed Sequence Tags (ESTs) generated represent a diverse sample drawn from multiple populations, developmental stages, and stress treatments. Despite this diversity, > 95% of the ESTs assembled into long (> 714 bp on average) and highly covered (> 9.6x on average) contigs. To estimate the effectiveness of transcript discovery, we compared the number of bases in the hit region of unigenes (contigs and singletons) to the length of the best match silkworm (Bombyx mori) protein--this "ortholog hit ratio" gives a close estimate on the amount of the transcript discovered relative to a model lepidopteran genome. For each species, we tested two assembly programs and two parameter sets; although CAP3 is commonly used for such data, the assemblies produced by Celera Assembler with modified parameters were chosen over those produced by CAP3 based on contig and singleton counts as well as ortholog hit ratio analysis. In the final assemblies, 1,413 E. propertius and 1,940 P. zelicaon unigenes had a ratio > 0.8; 2,866 E. propertius and 4,015 P. zelicaon unigenes had a ratio > 0.5.

CONCLUSIONS

Ultimately, these assemblies and SNP data will be used to generate microarrays for ecoinformatics examining climate change tolerance of different natural populations. These studies will benefit from high quality assemblies with few singletons (less than 26% of bases for each assembled transcriptome are present in unassembled singleton ESTs) and effective transcript discovery (over 6,500 of our putative orthologs cover at least 50% of the corresponding model silkworm gene).

Combined EST and proteomic analysis identifies rapidly evolving seminal fluid proteins in Heliconius butterflies

Seminal fluid proteins (SFPs) directly influence a wide range of reproductive processes, including fertilization, sperm storage, egg production, and immune response. Like many other reproductive proteins, the molecular evolution of SFPs is generally characterized by rapid and frequently adaptive evolution. However, the evolutionary processes underlying this often-documented pattern have not yet been confidently determined. A robust understanding of the processes governing SFP evolution will ultimately require identifying SFPs and characterizing their evolution in many different taxa, often where only limited genomic resources are available. Here, we report the first comprehensive molecular genetic and evolutionary analysis of SFPs conducted in Lepidoptera (moths and butterflies). We have identified 51 novel SFPs from two species of Heliconius butterflies (Heliconius erato and Heliconius melpomene) by combining "indirect" bioinformatic and expression analyses of expressed sequence tags from male accessory gland and wing tissues with "direct" proteomic analyses of spermatophores. Proteomic analyses identified fewer SFPs than the indirect criteria but gave consistent results. Of 51 SFPs, 40 were identified in both species but fewer than half could be functionally annotated via similarity searches (Blast, IPRscan, etc.). The majority of annotated Heliconius SFPs were predicted to be chymotrypsins. Comparisons of Heliconius SFPs with those from fruit fly, mosquito, honeybee, and cricket suggest that gene turnover is high among these proteins and that SFPs are rarely conserved across insect orders. Pairwise estimates of evolutionary rates between SFPs and nonreproductive proteins show that, on average, Heliconius SFPs are evolving rapidly. At least one of these SFPs is evolving adaptively (dN/dS > 1), implicating a role for positive selection in this rapid evolution. This work establishes a strong precedent for future research on the causes and consequences of reproductive protein evolution in the Lepidoptera. Butterflies and moths have an extremely rich history of organismal research, which will provide an informative ecological context for further molecular evolutionary investigations.

Next generation transcriptomes for next generation genomes using est2assembly

Background

The decreasing costs of capillary-based Sanger sequencing and next generation technologies, such as 454 pyrosequencing, have prompted an explosion of transcriptome projects in non-model species, where even shallow sequencing of transcriptomes can now be used to examine a range of research questions. This rapid growth in data has outstripped the ability of researchers working on non-model species to analyze and mine transcriptome data efficiently.

Results

Here we present a semi-automated platform 'est2assembly' that processes raw sequence data from Sanger or 454 sequencing into a hybrid de-novo assembly, annotates it and produces GMOD compatible output, including a SeqFeature database suitable for GBrowse. Users are able to parameterize assembler variables, judge assembly quality and determine the optimal assembly for their specific needs. We used est2assembly to process Drosophila and Bicyclus public Sanger EST data and then compared them to published 454 data as well as eight new insect transcriptome collections.

Conclusions

Analysis of such a wide variety of data allows us to understand how these new technologies can assist EST project design. We determine that assembler parameterization is as essential as standardized methods to judge the output of ESTs projects. Further, even shallow sequencing using 454 produces sufficient data to be of wide use to the community. est2assembly is an important tool to assist manual curation for gene models, an important resource in their own right but especially for species which are due to acquire a genome project using Next Generation Sequencing.

454 pyrosequencing based transcriptome analysis of Zygaena filipendulae with focus on genes involved in biosynthesis of cyanogenic glucosides

Background

An essential driving component in the co-evolution of plants and insects is the ability to produce and handle bioactive compounds. Plants produce bioactive natural products for defense, but some insects detoxify and/or sequester the compounds, opening up for new niches with fewer competitors. To study the molecular mechanism behind the co-adaption in plant-insect interactions, we have investigated the interactions between Lotus corniculatus and Zygaena filipendulae. They both contain cyanogenic glucosides which liberate toxic hydrogen cyanide upon breakdown. Moths belonging to the Zygaena family are the only insects known, able to carry out both de novo biosynthesis and sequestration of the same cyanogenic glucosides as those from their feed plants. The biosynthetic pathway for cyanogenic glucoside biosynthesis in Z. filipendulae proceeds using the same intermediates as in the well known pathway from plants, but none of the enzymes responsible have been identified. A genomics strategy founded on 454 pyrosequencing of the Z. filipendulae transcriptome was undertaken to identify some of these enzymes in Z. filipendulae.

Results

Comparisons of the Z. filipendulae transcriptome with the sequenced genomes of Bombyx mori, Drosophila melanogaster, Tribolium castaneum, Apis mellifera and Anopheles gambiae indicate a high coverage of the Z. filipendulae transcriptome. 11% of the Z. filipendulae transcriptome sequences were assigned to Gene Ontology categories. Candidate genes for enzymes functioning in the biosynthesis of cyanogenic glucosides (cytochrome P450 and family 1 glycosyltransferases) were identified based on sequence length, number of copies and presence/absence of close homologs in D. melanogaster, B. mori and the cyanogenic butterfly Heliconius. Examination of biased codon usage, GC content and selection on gene candidates support the notion of cyanogenesis as an "old" trait within Ditrysia, as well as its origins being convergent between plants and insects.

Conclusion

Pyrosequencing is an attractive approach to gain access to genes in the biosynthesis of bio-active natural products from insects and other organisms, for which the genome sequence is not known. Based on analysis of the Z. filipendulae transcriptome, promising gene candidates for biosynthesis of cyanogenic glucosides was identified, and the suitability of Z. filipendulae as a model system for cyanogenesis in insects is evident.

Obtaining accurate translations from expressed sequence tags

The genomes of an increasing number of species are being investigated through the generation of expressed sequence tags (ESTs). However, ESTs are prone to sequencing errors and typically define incomplete transcripts, making downstream annotation difficult. Annotation would be greatly improved with robust polypeptide translations. Many current solutions for EST translation require a large number of full-length gene sequences for training purposes, a resource that is not available for the majority of EST projects. As part of our ongoing EST programs investigating these "neglected" genomes, we have developed a polypeptide prediction pipeline, prot4EST. It incorporates freely available software to produce final translations that are more accurate than those derived from any single method. We describe how this integrated approach goes a long way to overcoming the deficit in training data.

Expressed sequence tags from larval gut of the European corn borer (Ostrinia nubilalis): exploring candidate genes potentially involved in Bacillus thuringiensis toxicity and resistance

Background

Lepidoptera represents more than 160,000 insect species which include some of the most devastating pests of crops, forests, and stored products. However, the genomic information on lepidopteran insects is very limited. Only a few studies have focused on developing expressed sequence tag (EST) libraries from the guts of lepidopteran larvae. Knowledge of the genes that are expressed in the insect gut are crucial for understanding basic physiology of food digestion, their interactions with Bacillus thuringiensis (Bt) toxins, and for discovering new targets for novel toxins for use in pest management. This study analyzed the ESTs generated from the larval gut of the European corn borer (ECB, Ostrinia nubilalis), one of the most destructive pests of corn in North America and the western world. Our goals were to establish an ECB larval gut-specific EST database as a genomic resource for future research and to explore candidate genes potentially involved in insect-Bt interactions and Bt resistance in ECB.

Results

We constructed two cDNA libraries from the guts of the fifth-instar larvae of ECB and sequenced a total of 15,000 ESTs from these libraries. A total of 12,519 ESTs (83.4%) appeared to be high quality with an average length of 656 bp. These ESTs represented 2,895 unique sequences, including 1,738 singletons and 1,157 contigs. Among the unique sequences, 62.7% encoded putative proteins that shared significant sequence similarities (E-value ≤ 10^-3)with the sequences available in GenBank. Our EST analysis revealed 52 candidate genes that potentially have roles in Bt toxicity and resistance. These genes encode 18 trypsin-like proteases, 18 chymotrypsin-like proteases, 13 aminopeptidases, 2 alkaline phosphatases and 1 cadherin-like protein. Comparisons of expression profiles of 41 selected candidate genes between Cry1Ab-susceptible and resistant strains of ECB by RT-PCR showed apparently decreased expressions in 2 trypsin-like and 2 chymotrypsin-like protease genes, and 1 aminopeptidase genes in the resistant strain as compared with the susceptible strain. In contrast, the expression of 3 trypsin- like and 3 chymotrypsin-like protease genes, 2 aminopeptidase genes, and 2 alkaline phosphatase genes were increased in the resistant strain. Such differential expressions of the candidate genes may suggest their involvement in Cry1Ab resistance. Indeed, certain trypsin-like and chymotrypsin-like proteases have previously been found to activate or degrade Bt protoxins and toxins, whereas several aminopeptidases, cadherin-like proteins and alkaline phosphatases have been demonstrated to serve as Bt receptor proteins in other insect species.

Conclusion

We developed a relatively large EST database consisting of 12,519 high-quality sequences from a total of 15,000 cDNAs from the larval gut of ECB. To our knowledge, this database represents the largest gut-specific EST database from a lepidopteran pest. Our work provides a foundation for future research to develop an ECB gut-specific DNA microarray which can be used to analyze the global changes of gene expression in response to Bt protoxins/toxins and the genetic difference(s) between Bt- resistant and susceptible strains. Furthermore, we identified 52 candidate genes that may potentially be involved in Bt toxicity and resistance. Differential expressions of 15 out of the 41 selected candidate genes examined by RT-PCR, including 5 genes with apparently decreased expression and 10 with increased expression in Cry1Ab-resistant strain, may help us conclusively identify the candidate genes involved in Bt resistance and provide us with new insights into the mechanism of Cry1Ab resistance in ECB.

A gene-based linkage map for Bicyclus anynana butterflies allows for a comprehensive analysis of synteny with the lepidopteran reference genome

Lepidopterans (butterflies and moths) are a rich and diverse order of insects, which, despite their economic impact and unusual biological properties, are relatively underrepresented in terms of genomic resources. The genome of the silkworm Bombyx mori has been fully sequenced, but comparative lepidopteran genomics has been hampered by the scarcity of information for other species. This is especially striking for butterflies, even though they have diverse and derived phenotypes (such as color vision and wing color patterns) and are considered prime models for the evolutionary and developmental analysis of ecologically relevant, complex traits. We focus on Bicyclus anynana butterflies, a laboratory system for studying the diversification of novelties and serially repeated traits. With a panel of 12 small families and a biphasic mapping approach, we first assigned 508 expressed genes to segregation groups and then ordered 297 of them within individual linkage groups. We also coarsely mapped seven color pattern loci. This is the richest gene-based map available for any butterfly species and allowed for a broad-coverage analysis of synteny with the lepidopteran reference genome. Based on 462 pairs of mapped orthologous markers in Bi. anynana and Bo. mori, we observed strong conservation of gene assignment to chromosomes, but also evidence for numerous large- and small-scale chromosomal rearrangements. With gene collections growing for a variety of target organisms, the ability to place those genes in their proper genomic context is paramount. Methods to map expressed genes and to compare maps with relevant model systems are crucial to extend genomic-level analysis outside classical model species. Maps with gene-based markers are useful for comparative genomics and to resolve mapped genomic regions to a tractable number of candidate genes, especially if there is synteny with related model species. This is discussed in relation to the identification of the loci contributing to color pattern evolution in butterflies.

An expressed sequence tag (EST) library for Drosophila serrata, a model system for sexual selection and climatic adaptation studies

BACKGROUND

The native Australian fly Drosophila serrata belongs to the highly speciose montium subgroup of the melanogaster species group. It has recently emerged as an excellent model system with which to address a number of important questions, including the evolution of traits under sexual selection and traits involved in climatic adaptation along latitudinal gradients. Understanding the molecular genetic basis of such traits has been limited by a lack of genomic resources for this species. Here, we present the first expressed sequence tag (EST) collection for D. serrata that will enable the identification of genes underlying sexually-selected phenotypes and physiological responses to environmental change and may help resolve controversial phylogenetic relationships within the montium subgroup.

RESULTS

A normalized cDNA library was constructed from whole fly bodies at several developmental stages, including larvae and adults. Assembly of 11,616 clones sequenced from the 3' end allowed us to identify 6,607 unique contigs, of which at least 90% encoded peptides. Partial transcripts were discovered from a variety of genes of evolutionary interest by BLASTing contigs against the 12 Drosophila genomes currently sequenced. By incorporating into the cDNA library multiple individuals from populations spanning a large portion of the geographical range of D. serrata, we were able to identify 11,057 putative single nucleotide polymorphisms (SNPs), with 278 different contigs having at least one "double hit" SNP that is highly likely to be a real polymorphism. At least 394 EST-associated microsatellite markers, representing 355 different contigs, were also found, providing an additional set of genetic markers. The assembled EST library is available online at http://www.chenowethlab.org/serrata/index.cgi.

CONCLUSION

We have provided the first gene collection and largest set of polymorphic genetic markers, to date, for the fly D. serrata. The EST collection will provide much needed genomic resources for this model species and facilitate comparative evolutionary studies within the montium subgroup of the D. melanogaster lineage.

EST analysis of male accessory glands from Heliconius butterflies with divergent mating systems

Background

Heliconius butterflies possess a remarkable diversity of phenotypes, physiologies, and behaviors that has long distinguished this genus as a focal taxon in ecological and evolutionary research. Recently Heliconius has also emerged as a model system for using genomic methods to investigate the causes and consequences of biological diversity. One notable aspect of Heliconius diversity is a dichotomy in mating systems which provides an unusual opportunity to investigate the relationship between sexual selection and the evolution of reproductive proteins. As a first step in pursuing this research, we report the generation and analysis of expressed sequence tags (ESTs) from the male accessory gland of H. erato and H. melpomene, species representative of the two mating systems present in the genus Heliconius.

Results

We successfully sequenced 933 ESTs clustering into 371 unigenes from H. erato and 1033 ESTs clustering into 340 unigenes from H. melpomene. Results from the two species were very similar. Approximately a third of the unigenes showed no significant BLAST similarity (E-value <10^-5) to sequences in GenBank's non-redundant databases, indicating that a large proportion of novel genes are expressed in Heliconius male accessory glands. In both species only a third of accessory gland unigenes were also found among genes expressed in wing tissue. About 25% of unigenes from both species encoded secreted proteins. This includes three groups of highly abundant unigenes encoding repetitive proteins considered to be candidate seminal fluid proteins; proteins encoded by one of these groups were detected in H. erato spermatophores.

Conclusion

This collection of ESTs will serve as the foundation for the future identification and evolutionary analysis of male reproductive proteins in Heliconius butterflies. These data also represent a significant advance in the rapidly growing collection of genomic resources available in Heliconius butterflies. As such, they substantially enhance this taxon as a model system for investigating questions of ecological, phenotypic, and genomic diversity.

annot8r: GO, EC and KEGG annotation of EST datasets

Background

The expressed sequence tag (EST) methodology is an attractive option for the generation of sequence data for species for which no completely sequenced genome is available. The annotation and comparative analysis of such datasets poses a formidable challenge for research groups that do not have the bioinformatics infrastructure of major genome sequencing centres. Therefore, there is a need for user-friendly tools to facilitate the annotation of non-model species EST datasets with well-defined ontologies that enable meaningful cross-species comparisons. To address this, we have developed annot8r, a platform for the rapid annotation of EST datasets with GO-terms, EC-numbers and KEGG-pathways.

Results

annot8r automatically downloads all files relevant for the annotation process and generates a reference database that stores UniProt entries, their associated Gene Ontology (GO), Enzyme Commission (EC) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) annotation and additional relevant data. For each of GO, EC and KEGG, annot8r extracts a specific sequence subset from the UniProt dataset based on the information stored in the reference database. These three subsets are then formatted for BLAST searches. The user provides the protein or nucleotide sequences to be annotated and annot8r runs BLAST searches against these three subsets. The BLAST results are parsed and the corresponding annotations retrieved from the reference database. The annotations are saved both as flat files and also in a relational postgreSQL results database to facilitate more advanced searches within the results. annot8r is integrated with the PartiGene suite of EST analysis tools.

Conclusion

annot8r is a tool that assigns GO, EC and KEGG annotations for data sets resulting from EST sequencing projects both rapidly and efficiently. The benefits of an underlying relational database, flexibility and the ease of use of the program make it ideally suited for non-model species EST-sequencing projects.

Developmental and genetic mechanisms for evolutionary diversification of serial repeats: eyespot size in Bicyclus anynana butterflies

Serially repeated pattern elements on butterfly wings offer the opportunity for integrating genetic, developmental, and functional aspects towards understanding morphological diversification and the evolution of individuality. We use captive populations of Bicyclus anynana butterflies, an emerging model in evolutionary developmental biology, to explore the genetic and developmental basis of compartmentalized changes in eyespot patterns. There is much variation for different aspects of eyespot morphology, and knowledge about the genetic pathways and developmental processes involved in eyespot formation. Also, despite the strong correlations across all eyespots in one butterfly, B. anynana shows great potential for independent changes in the size of individual eyespots. It is, however, unclear to what extent the genetic and developmental processes underlying eyespot formation change in a localized manner to enable such individualization. We use micromanipulations of developing wings to dissect the contribution of different components of eyespot development to quantitative differences in eyespot size on one wing surface. Reciprocal transplants of presumptive eyespot foci between artificial selection lines and controls suggest that while localized antagonistic changes in eyespot size rely mostly on localized changes in focal signal strength, concerted changes depend greatly on epidermal response sensitivities. This potentially reflects differences between the signal-response components of eyespot formation in the degrees of compartmentalization and/or the temporal pattern of selection. We also report on the phenotypic analysis of a number of mutant stocks demonstrating how single alleles can affect different eyespots in concert or independently, and thus contribute to the individualization of serially repeated traits.

Rapid transcriptome characterization for a nonmodel organism using 454 pyrosequencing

We present a de novo assembly of a eukaryote transcriptome using 454 pyrosequencing data. The Glanville fritillary butterfly (Melitaea cinxia; Lepidoptera: Nymphalidae) is a prominent species in population biology but had no previous genomic data. Sequencing runs using two normalized complementary DNA collections from a genetically diverse pool of larvae, pupae, and adults yielded 608,053 expressed sequence tags (mean length = 110 nucleotides), which assembled into 48,354 contigs (sets of overlapping DNA segments) and 59,943 singletons. BLAST comparisons confirmed the accuracy of the sequencing and assembly, and indicated the presence of c. 9000 unique genes, along with > 6000 additional microarray-confirmed unannotated contigs. Average depth of coverage was 6.5-fold for the longest 4800 contigs (348-2849 bp in length), sufficient for detecting large numbers of single nucleotide polymorphisms. Oligonucleotide microarray probes designed from the assembled sequences showed highly repeatable hybridization intensity and revealed biological differences among individuals. We conclude that 454 sequencing, when performed to provide sufficient coverage depth, allows de novo transcriptome assembly and a fast, cost-effective, and reliable method for development of functional genomic tools for nonmodel species. This development narrows the gap between approaches based on model organisms with rich genetic resources vs. species that are most tractable for ecological and evolutionary studies.

Conserved synteny of genes between chromosome 15 of Bombyx mori and a chromosome of Manduca sexta shown by five-color BAC-FISH

The successful assignment of the existing genetic linkage groups (LGs) to individual chromosomes and the second-generation linkage map obtained by mapping a large number of bacterial artificial chromosome (BAC) contigs in the silkworm, Bombyx mori , together with public nucleotide sequence databases, offer a powerful tool for the study of synteny between karyotypes of B. mori and other lepidopteran species. Conserved synteny of genes between particular chromosomes can be identified by comparatively mapping orthologous genes of the corresponding linkage groups with the help of BAC-FISH (fluorescent in situ hybridization). This technique was established in B. mori for 2 differently labeled BAC probes simultaneously hybridized to pachytene bivalents. To achieve higher-throughput comparative mapping using BAC-FISH in Lepidoptera, we developed a protocol for five-color BAC-FISH, which allowed us to map simultaneously 6 different BAC probes to chromosome 15 in B. mori. We identified orthologs of 6 B. mori LG15 genes (RpP0, RpS8, eIF3, RpL7A, RpS23, and Hsc70) for the tobacco hornworm, Manduca sexta , and selected the ortholog-containing BAC clones from an M. sexta BAC library. All 6 M. sexta BAC clones hybridized to a single M. sexta bivalent in pachytene spermatocytes. Thus, we have confirmed the conserved synteny between the B. mori chromosome 15 and the corresponding M. sexta chromosome (hence provisionally termed chromosome 15).

ButterflyBase: a platform for lepidopteran genomics

With over 100 000 species and a large community of evolutionary biologists, population ecologists, pest biologists and genome researchers, the Lepidoptera are an important insect group. Genomic resources [expressed sequence tags (ESTs), genome sequence, genetic and physical maps, proteomic and microarray datasets] are growing, but there has up to now been no single access and analysis portal for this group. Here we present ButterflyBase (http://www.butterflybase.org), a unified resource for lepidopteran genomics. A total of 273 077 ESTs from more than 30 different species have been clustered to generate stable unigene sets, and robust protein translations derived from each unigene cluster. Clusters and their protein translations are annotated with BLAST-based similarity, gene ontology (GO), enzyme classification (EC) and Kyoto encyclopaedia of genes and genomes (KEGG) terms, and are also searchable using similarity tools such as BLAST and MS-BLAST. The database supports many needs of the lepidopteran research community, including molecular marker development, orthologue prediction for deep phylogenetics, and detection of rapidly evolving proteins likely involved in host–pathogen or other evolutionary processes. ButterflyBase is expanding to include additional genomic sequence, ecological and mapping data for key species.

Synteny and chromosome evolution in the lepidoptera: evidence from mapping in Heliconius melpomene

The extent of conservation of synteny and gene order in the Lepidoptera has been investigated previously only by comparing a small subset of linkage groups between the moth Bombyx mori and the butterfly Heliconius melpomene. Here we report the mapping of 64 additional conserved genes in H. melpomene, which contributed 47 markers to a comparative framework of 72 orthologous loci spanning all 21 H. melpomene chromosomes and 27 of the 28 B. mori chromosomes. Comparison of the maps revealed conserved synteny across all chromosomes for the 72 loci, as well as evidence for six cases of chromosome fusion in the Heliconius lineage that contributed to the derived 21-chromosome karyotype. Comparisons of gene order on these fused chromosomes revealed two instances of colinearity between H. melpomene and B. mori, but also one instance of likely chromosomal rearrangement. B. mori is the first lepidopteran species to have its genome sequenced, and the finding that there is conserved synteny and gene order among Lepidoptera indicates that the genomic tools developed in B. mori will be broadly useful in other species.

Butterfly genomics eclosing

Technological and conceptual advances of the last decade have led to an explosion of genomic data and the emergence of new research avenues. Evolutionary and ecological functional genomics, with its focus on the genes that affect ecological success and adaptation in natural populations, benefits immensely from a phylogenetically widespread sampling of biological patterns and processes. Among those organisms outside established model systems, butterflies offer exceptional opportunities for multidisciplinary research on the processes generating and maintaining variation in ecologically relevant traits. Here we highlight research on wing color pattern variation in two groups of Nymphalid butterflies, the African species Bicyclus anynana (subfamily Satyrinae) and species of the South American genus Heliconius (subfamily Heliconiinae), which are emerging as important systems for studying the nature and origins of functional diversity. Growing genomic resources including genomic and cDNA libraries, dense genetic maps, high-density gene arrays, and genetic transformation techniques are extending current gene mapping and expression profiling analysis and enabling the next generation of research questions linking genes, development, form, and fitness. Efforts to develop such resources in Bicyclus and Heliconius underscore the general challenges facing the larger research community and highlight the need for a community-wide effort to extend ongoing functional genomic research on butterflies.

Transgenic approaches to study wing color pattern development in Lepidoptera

The extremely diverse lepidopteran wing patterns make useful models to study the evolution of development and the molecular changes that enable it. Until now, the implication of candidate genes in the differentiation of color patterns has relied primarily on correlational evidence, i.e., gene expression patterns in a developing wing mapping closely to the adult color pattern. The use of transgenic techniques in the Lepidoptera, including the manipulation of gene expression, will finally allow researchers to test hypotheses of gene function at various levels of the patterning hierarchy, from signaling ligands and transcription factors to pigment enzymes. Here we present an overview of transgenic techniques employed in lepidopteran systems and highlight areas where current and future research will provide exciting opportunities to deepen our understanding of the mechanisms of morphological evolution.

A conserved supergene locus controls colour pattern diversity in Heliconius butterflies

We studied whether similar developmental genetic mechanisms are involved in both convergent and divergent evolution. Mimetic insects are known for their diversity of patterns as well as their remarkable evolutionary convergence, and they have played an important role in controversies over the respective roles of selection and constraints in adaptive evolution. Here we contrast three butterfly species, all classic examples of Müllerian mimicry. We used a genetic linkage map to show that a locus, Yb, which controls the presence of a yellow band in geographic races of Heliconius melpomene, maps precisely to the same location as the locus Cr, which has very similar phenotypic effects in its co-mimic H. erato. Furthermore, the same genomic location acts as a "supergene", determining multiple sympatric morphs in a third species, H. numata. H. numata is a species with a very different phenotypic appearance, whose many forms mimic different unrelated ithomiine butterflies in the genus Melinaea. Other unlinked colour pattern loci map to a homologous linkage group in the co-mimics H. melpomene and H. erato, but they are not involved in mimetic polymorphism in H. numata. Hence, a single region from the multilocus colour pattern architecture of H. melpomene and H. erato appears to have gained control of the entire wing-pattern variability in H. numata, presumably as a result of selection for mimetic "supergene" polymorphism without intermediates. Although we cannot at this stage confirm the homology of the loci segregating in the three species, our results imply that a conserved yet relatively unconstrained mechanism underlying pattern switching can affect mimicry in radically different ways. We also show that adaptive evolution, both convergent and diversifying, can occur by the repeated involvement of the same genomic regions.

Heliconius wing patterns: an evo-devo model for understanding phenotypic diversity

Evolutionary Developmental Biology aims for a mechanistic understanding of phenotypic diversity, and present knowledge is largely based on gene expression and interaction patterns from a small number of well-known model organisms. However, our understanding of biological diversification depends on our ability to pinpoint the causes of natural variation at a micro-evolutionary level, and therefore requires the isolation of genetic and developmental variation in a controlled genetic background. The colour patterns of Heliconius butterflies (Nymphalidae: Heliconiinae) provide a rich suite of naturally occurring variants with striking phenotypic diversity and multiple taxonomic levels of variation. Diversification in the genus is well known for its dramatic colour-pattern divergence between races or closely related species, and for Müllerian mimicry convergence between distantly related species, providing a unique system to study the development basis of colour-pattern evolution. A long history of genetic studies has showed that pattern variation is based on allelic combinations at a surprisingly small number of loci, and recent developmental evidence suggests that pattern development in Heliconius is different from the eyespot determination of other butterflies. Fine-scale genetic mapping studies have shown that a shared toolkit of genes is used to produce both convergent and divergent phenotypes. These exciting results and the development of new genomic resources make Heliconius a very promising evo-devo model for the study of adaptive change.

SPODOBASE: an EST database for the lepidopteran crop pest Spodoptera

Background

The Lepidoptera Spodoptera frugiperda is a pest which causes widespread economic damage on a variety of crop plants. It is also well known through its famous Sf9 cell line which is used for numerous heterologous protein productions. Species of the Spodoptera genus are used as model for pesticide resistance and to study virus host interactions. A genomic approach is now a critical step for further new developments in biology and pathology of these insects, and the results of ESTs sequencing efforts need to be structured into databases providing an integrated set of tools and informations.

Description

The ESTs from five independent cDNA libraries, prepared from three different S. frugiperda tissues (hemocytes, midgut and fat body) and from the Sf9 cell line, are deposited in the database. These tissues were chosen because of their importance in biological processes such as immune response, development and plant/insect interaction. So far, the SPODOBASE contains 29,325 ESTs, which are cleaned and clustered into non-redundant sets (2294 clusters and 6103 singletons). The SPODOBASE is constructed in such a way that other ESTs from S. frugiperda or other species may be added. User can retrieve information using text searches, pre-formatted queries, query assistant or blast searches. Annotation is provided against NCBI, UNIPROT or Bombyx mori ESTs databases, and with GO-Slim vocabulary.

Conclusion

The SPODOBASE database provides integrated access to expressed sequence tags (EST) from the lepidopteran insect Spodoptera frugiperda. It is a publicly available structured database with insect pest sequences which will allow identification of a number of genes and comprehensive cloning of gene families of interest for scientific community. SPODOBASE is available from URL:

A wing expressed sequence tag resource for Bicyclus anynana butterflies, an evo-devo model

BACKGROUND

Butterfly wing color patterns are a key model for integrating evolutionary developmental biology and the study of adaptive morphological evolution. Yet, despite the biological, economical and educational value of butterflies they are still relatively under-represented in terms of available genomic resources. Here, we describe an Expression Sequence Tag (EST) project for Bicyclus anynana that has identified the largest available collection to date of expressed genes for any butterfly.

RESULTS

By targeting cDNAs from developing wings at the stages when pattern is specified, we biased gene discovery towards genes potentially involved in pattern formation. Assembly of 9,903 ESTs from a subtracted library allowed us to identify 4,251 genes of which 2,461 were annotated based on BLAST analyses against relevant gene collections. Gene prediction software identified 2,202 peptides, of which 215 longer than 100 amino acids had no homology to any known proteins and, thus, potentially represent novel or highly diverged butterfly genes. We combined gene and Single Nucleotide Polymorphism (SNP) identification by constructing cDNA libraries from pools of outbred individuals, and by sequencing clones from the 3' end to maximize alignment depth. Alignments of multi-member contigs allowed us to identify over 14,000 putative SNPs, with 316 genes having at least one high confidence double-hit SNP. We furthermore identified 320 microsatellites in transcribed genes that can potentially be used as genetic markers.

CONCLUSION

Our project was designed to combine gene and sequence polymorphism discovery and has generated the largest gene collection available for any butterfly and many potential markers in expressed genes. These resources will be invaluable for exploring the potential of B. anynana in particular, and butterflies in general, as models in ecological, evolutionary, and developmental genetics.

Localization of Müllerian mimicry genes on a dense linkage map of Heliconius erato

We report a dense genetic linkage map of Heliconius erato, a neotropical butterfly that has undergone a remarkable adaptive radiation in warningly colored mimetic wing patterns. Our study exploited natural variation segregating in a cross between H. erato etylus and H. himera to localize wing color pattern loci on a dense linkage map containing amplified fragment length polymorphisms (AFLP), microsatellites, and single-copy nuclear loci. We unambiguously identified all 20 autosomal linkage groups and the sex chromosome (Z). The map spanned a total of 1430 Haldane cM and linkage groups varied in size from 26.3 to 97.8 cM. The average distance between markers was 5.1 cM. Within this framework, we localized two major color pattern loci to narrow regions of the genome. The first gene, D, responsible for red/orange elements, had a most likely placement in a 6.7-cM region flanked by two AFLP markers on the end of a large 87.5-cM linkage group. The second locus, Sd, affects the melanic pattern on the forewing and was found within a 6.3-cM interval between flanking AFLP loci. This study complements recent linkage analysis of H. erato's comimic, H. melpomene, and forms the basis for marker-assisted physical mapping and for studies into the comparative genetic architecture of wing-pattern mimicry in Heliconius.

Simple sequence repeat-based consensus linkage map of Bombyx mori

We established a genetic linkage map employing 518 simple sequence repeat (SSR, or microsatellite) markers for Bombyx mori (silkworm), the economically and culturally important lepidopteran insect, as part of an international genomics program. A survey of six representative silkworm strains using 2,500 (CA)n- and (CT)n-based SSR markers revealed 17-24% polymorphism, indicating a high degree of homozygosity resulting from a long history of inbreeding. Twenty-nine SSR linkage groups were established in well characterized Dazao and C108 strains based on genotyping of 189 backcross progeny derived from an F(1) male mated with a C108 female. The clustering was further focused to 28 groups by genotyping 22 backcross progeny derived from an F(1) female mated with a C108 male. This set of SSR linkage groups was further assigned to the 28 chromosomes (established linkage groups) of silkworm aided by visible mutations and cleaved amplified polymorphic sequence markers developed from previously mapped genes, cDNA sequences, and cloned random amplified polymorphic DNAs. By integrating a visible mutation p (plain, larval marking) and 29 well conserved genes of insects onto this SSR-based linkage map, a second generation consensus silkworm genetic map with a range of 7-40 markers per linkage group and a total map length of approximately 3431.9 cM was constructed and its high efficiency for genotyping and potential application for synteny studies of Lepidoptera and other insects was demonstrated.