Mining an Ostrinia nubilalis midgut expressed sequence tag (EST) library for candidate genes and single nucleotide polymorphisms (SNPs)

Identification of a metallothionein gene and the role of biological thiols in stress induced by short-term Cd exposure in Ostrinia nubilalis

Cadmium (Cd) is a non-essential metal that is highly toxic to all living forms, characterized by an extremely high affinity for thiol (SH) groups. The aim of this work was to identify and experimentally verify metallothionein gene and to analyze the role of biological thiols in stress induced by short-term Cd exposure in Ostrinia nubilalis, one of the most important corn pests. The coding region of a metallothionein (MT) gene in O. nubilalis was identified, encoding protein, OnMT1, which contains 46 amino acids, including 12 cysteine residues, and has no aromatic amino acids. Phylogenetic analysis revealed that OnMT1 clustered together with metallothionein from Bombyx mori. Structural bioinformatics analysis strongly suggests that OnMT1 is a metallothionein with affinity for multiple transition metals. Further, in order to elucidate the role of biological thiols, O. nubilalis L5 larvae were exposed to increasing Cd concentrations in diet (6.85, 41.71, 77.35 mg kg^-1) during a 48 h period, after which Cd concentration in larvae was measured (3.50, 12.02, 47.37 mg kg^-1, respectively). Due to short-term Cd exposure, concentration of free protein SH groups and relative expression of OnMT1 and thioredoxin (Trx) genes was elevated, while the reduced glutathione content remained unchanged. The presented results provide evidence that OnMT1 plays a role in Cd detoxification and homeostasis, and confirm the importance of biological thiols, especially OnMT1 and Trx, in the early response of O. nubilalis to Cd poisoning, indicating interaction between Cd and thiol-linked redox reactions. Insects provide valuable insight into molecular adaptations to metals.

Newly identified APN splice isoforms suggest novel splicing mechanisms may underlie circRNA circularization in moth

Circular RNA (circRNA) have long been considered by-products of splicing errors, but the coordination of RNA transcription and exon circularization events remains poorly understood. Here, we investigated this question using genes encoding aminopeptidases N (APNs), which are receptors of Bacillus thuringiensis toxins, in the cotton bollworm, Helicoverpa armigera. We cloned and sequenced the cDNA of ten APN genes (HaAPN1-10) located in the same APN gene cluster, and detected 20 and 14 novel splicing isoforms with exon skipping in HaAPN1 and HaAPN3, respectively, whereas no or very few variants were found in the remaining genes. Further study identified 14 and 6 circular RNA (circRNA) in HaAPN1 and HaAPN3, respectively. Neither novel splicing isoforms nor circRNA were detected in HaAPN2 and HaAPN5. Distinct from the conventional GT/AG splicing signal, short co-directional repeats were involved in the splicing of the linear and circular isoforms of HaAPN1 and HaAPN3. Identification of the splice sites revealed that the linear isoforms may be related in some way to the circularization. Moreover, phylogenetic analysis and detection of circRNA of the APN gene of the diamondback moth, Plutella xylostella (PxAPN3), showed that circRNA formation is relatively conserved during the lepidopteran evolutionary process. These results contribute to an improved understanding of lepidopteran APNs and this novel class of insect circRNA.

Down-regulation of aminopeptidase N and ABC transporter subfamily G transcripts in Cry1Ab and Cry1Ac resistant Asian corn borer, Ostrinia furnacalis (Lepidoptera: Crambidae)

The Asian corn borer (ACB), Ostrinia furnacalis (Lepidoptera: Crambidae), is a highly destructive pest of cultivated maize throughout East Asia. Bacillus thuringiensis (Bt) crystalline protein (Cry) toxins cause mortality by a mechanism involving pore formation or signal transduction following toxin binding to receptors along the midgut lumen of susceptible insects, but this mechanism and mutations therein that lead to resistance are not fully understood. In the current study, quantitative comparisons were made among midgut expressed transcripts from O. furnacalis susceptible (ACB-BtS) and laboratory selected strains resistant to Cry1Ab (ACB-AbR) and Cry1Ac toxins (ACB-AcR) when feeding on non-Bt diet. From a combined de novo transcriptome assembly of 83,370 transcripts, ORFs of ≥ 100 amino acids were predicted and annotated for 28,940 unique isoforms derived from 12,288 transcripts. Transcriptome-wide expression estimated from RNA-seq read depths predicted significant down-regulation of transcripts for previously known Bt resistance genes, aminopeptidase N1 (apn1) and apn3, as well as a putative ATP binding cassette transporter group G (abcg) gene in both ACB-AbR and -AcR (log2[fold-change] ≥ 1.36; P < 0.0001). The transcripts that were most highly differentially regulated in both ACB-AbR and -AcR compared to ACB-BtS (log2[fold-change] ≥ 2.0; P < 0.0001) included up- and down-regulation of serine proteases, storage proteins and cytochrome P450 monooxygenases, as well as up-regulation of genes with predicted transport function. This study predicted the significant down-regulation of transcripts for previously known Bt resistance genes, aminopeptidase N1 (apn1) and apn3, as well as abccg gene in both ACB-AbR and -AcR. These data are important for the understanding of systemic differences between Bt resistant and susceptible genotypes.

Differential expression of chemosensory-protein genes in midguts in response to diet of Spodoptera litura

While it has been well characterized that chemosensory receptors in guts of mammals have great influence on food preference, much remains elusive in insects. Insect chemosensory proteins (CSPs) are soluble proteins that could deliver chemicals to olfactory and gustatory receptors. Recent studies have identified a number of CSPs expressed in midgut in Lepidoptera insects, which started to reveal their roles in chemical recognition and stimulating appetite in midgut. In this study, we examined expression patterns in midgut of 21 Spodoptera litura CSPs (SlitCSPs) characterized from a previously reported transcriptome, and three CSPs were identified to be expressed highly in midgut. The orthologous relationships between midgut expressed CSPs in S. litura and those in Bombyx mori and Plutella xylostella also suggest a conserved pattern of CSP expression in midgut. We further demonstrated that the expression of midgut-CSPs may change in response to different host plants, and SlitCSPs could bind typical chemicals from host plant in vitro. Overall, our results suggested midgut expressed SlitCSPs may have functional roles, likely contributing to specialization and adaption to different ecosystems. Better knowledge of this critical component of the chemsensation signaling pathways in midguts may improve our understanding of food preference processes in a new perspective.

Microarray-based annotation of the gut transcriptome of the migratory locust, Locusta migratoria

The migratory locust, Locusta migratoria, is a serious agricultural pest and important insect model in the study of insect digestion and feeding behaviour. The gut is one of the primary interfaces between the insect and its environment. Nevertheless, knowledge on the gut transcriptome of L. migratoria is still very limited. Here, 48 802 expressed sequence tags were extracted from publicly available databases and their expression in larval gut and/or brain tissue was determined using microarray hybridization. Our data show 2765 transcripts predominantly or exclusively expressed in the gut. Many transcripts had putative functions closely related to the physiological functions of the gut as a muscular digestive organ and as the first barrier against microorganisms and a wide range of toxins. By means of a ranking procedure based on the relative signal intensity, we estimated 15% of the transcripts to show high expression levels, the highest belonging to diverse digestive enzymes and muscle-related proteins. We also found evidence for very high expression of an allergen protein, which could have important implications, as locusts form a traditional food source in various parts of the world, and were also recently added to the list of insects fit for human consumption in Europe. Interestingly, many highly expressed sequences have as yet unknown functions. Taken together, the present data provide significant insight into locust larval gut physiology, and will be valuable for future studies on the insect gut.

Transcriptome Comparison Analysis of Ostrinia furnacalis in Four Developmental Stages

The Asian corn borer, Ostrinia furnacalis, is one of the most destructive pests of maize and causes huge losses in maize yield each year. In order to characterize the different developmental stages, a high-throughput sequencing platform was employed to perform de novo transcriptome assembly and gene expression analysis for the egg, larva, pupa and adult stages. Approximately 185 million reads were obtained, trimmed, and assembled into 42,638 unigenes with an average length of 801.94 bp and an N50 length of 1,152 bp. These unigene sequences were annotated and classified by performing Gene Ontology (GO), Cluster of Orthologous Groups (KOG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional classifications. Comparison of the gene expression profiles of the two transitional stages revealed dramatic differences. Some differentially expressed genes are associated with digestion, cuticularization olfactory recognition and wing formation as well as growth and development. In total, 12 putative insect development-related genes were identified. Real-time quantitative PCR (RT-qPCR) results and sequencing based on relative expression levels of randomly selected genes confirmed these expression patterns. These data represent the most comprehensive transcriptomic resource currently available for O. furnacalis and will facilitate the study of developmental pathways, cuticularization, wing formation and olfactory recognition.

Development and Characterization of Novel Microsatellite Markers for the Peach Fruit Moth Carposina sasakii (Lepidoptera: Carposinidae) Using Next-Generation Sequencing

The peach fruit moth Carposina sasakii is an economically important pest on dozens of fruits from Rosaceae and Rhamnaceae in Northeast Asia. We developed novel microsatellite markers for C. sasakii from randomly sequenced regions of the genome using next-generation sequencing. In total, 95,153 microsatellite markers were isolated from 4.70 GB genomic sequences. Thirty-five polymorphic markers were developed by assessing in 63 individuals from two geographical populations. The allele numbers ranged from 2 to 9 with an average value of 4.60 per locus, while the polymorphism information content ranged from 0.075 to 0.696 with an average value of 0.407. Furthermore, the observed and expected heterozygosity varied from 0.000 to 0.677 and 0.062 to 0.771, respectively. The microsatellites developed provide abundant molecular markers for investigating genetic structure, genetic diversity, and existence of host-plant associated biotypes of C. sasakii.

Unlinked genetic loci control the reduced transcription of aminopeptidase N 1 and 3 in the European corn borer and determine tolerance to Bacillus thuringiensis Cry1Ab toxin

Transgenic expression of Bacillus thuringiensis (Bt) crystalline (Cry) toxins by crop plants result in reduced insect feeding damage, but sustainability is threatened by the development of resistance traits in target insect populations. We investigated Bt toxin resistance trait in a laboratory colony of the European corn borer, Ostrinia nubilalis, selected for increased survival when exposed to Cry1Ab and correlated survival on Cry1Ab toxin with a constitutive ∼146.2 ± 17.3-fold reduction in midgut aminopeptidase N1 (apn1) transcript levels. A 7.1 ± 1.9-fold reduction apn3 transcript level was also correlated with Cry1Ab resistance. Quantitative trait locus (QTL) mapping identified a single major genome region controlling Cry1Ab resistance on linkage group 24 (LG24), and a minor QTL on LG27. Both QTL were independent of apn1 and apn3 loci on LG02. Positional mapping identified genetic markers that may assist in the identification of causal gene(s) within QTL intervals. This study indicates that genetic factor(s) may act in trans to reduce both apn1 and apn3 expression in Cry1Ab resistant O. nubilalis larvae, and suggest that gene regulatory pathways can influence Bt resistance traits. These findings show that gene interactions (epistasis) may influence Bt resistance in target insect populations.

Development of reference transcriptomes for the major field insect pests of cowpea: a toolbox for insect pest management approaches in west Africa

Cowpea is a widely cultivated and major nutritional source of protein for many people that live in West Africa. Annual yields and longevity of grain storage is greatly reduced by feeding damage caused by a complex of insect pests that include the pod sucking bugs, Anoplocnemis curvipes Fabricius (Hemiptera: Coreidae) and Clavigralla tomentosicollis Stål (Hemiptera: Coreidae); as well as phloem-feeding cowpea aphids, Aphis craccivora Koch (Hemiptera: Aphididae) and flower thrips, Megalurothrips sjostedti Trybom (Thysanoptera: Thripidae). Efforts to control these pests remain a challenge and there is a need to understand the structure and movement of these pest populations in order to facilitate the development of integrated pest management strategies (IPM). Molecular tools have the potential to help facilitate a better understanding of pest populations. Towards this goal, we used 454 pyrosequencing technology to generate 319,126, 176,262, 320,722 and 227,882 raw reads from A. curvipes, A. craccivora, C. tomentosicollis and M. sjostedti, respectively. The reads were de novo assembled into 11,687, 7,647, 10,652 and 7,348 transcripts for A. curvipes, A. craccivora, C. tomentosicollis and M. sjostedti, respectively. Functional annotation of the resulting transcripts identified genes putatively involved in insecticide resistance, pathogen defense and immunity. Additionally, sequences that matched the primary aphid endosymbiont, Buchnera aphidicola, were identified among A. craccivora transcripts. Furthermore, 742, 97, 607 and 180 single nucleotide polymorphisms (SNPs) were respectively predicted among A. curvipes, A. craccivora, C. tomentosicollis and M. sjostedti transcripts, and will likely be valuable tools for future molecular genetic marker development. These results demonstrate that Roche 454-based transcriptome sequencing could be useful for the development of genomic resources for cowpea pest insects in West Africa.

Application of Cydia pomonella expressed sequence tags: Identification and expression of three general odorant binding proteins in codling moth

The codling moth, Cydia pomonella, is one of the most important pests of pome fruits in the world, yet the molecular genetics and the physiology of this insect remain poorly understood. A combined assembly of 8 341 expressed sequence tags was generated from Roche 454 GS-FLX sequencing of eight tissue-specific cDNA libraries. Putative chemosensory proteins (12) and odorant binding proteins (OBPs) (18) were annotated, which included three putative general OBP (GOBP), one more than typically reported for other Lepidoptera. To further characterize CpomGOBPs, we cloned cDNA copies of their transcripts and determined their expression patterns in various tissues. Cloning and sequencing of the 698 nt transcript for CpomGOBP1 resulted in the prediction of a 163 amino acid coding region, and subsequent RT-PCR indicated that the transcripts were mainly expressed in antennae and mouthparts. The 1 289 nt (160 amino acid) CpomGOBP2 and the novel 702 nt (169 amino acid) CpomGOBP3 transcripts are mainly expressed in antennae, mouthparts, and female abdomen tips. These results indicate that next generation sequencing is useful for the identification of novel transcripts of interest, and that codling moth expresses a transcript encoding for a new member of the GOBP subfamily.

Frequency of hybridization between Ostrinia nubilalis E-and Z-pheromone races in regions of sympatry within the United States

Female European corn borer, Ostrinia nubilalis, produce and males respond to sex pheromone blends with either E- or Z-Δ11-tetradecenyl acetate as the major component. E- and Z-race populations are sympatric in the Eastern United States, Southeastern Canada, and the Mediterranean region of Europe. The E- and Z-pheromone races of O. nubilalis are models for incipient species formation, but hybridization frequencies within natural populations remain obscure due to lack of a high-throughput phenotyping method. Lassance et al. previously identified a pheromone gland-expressed fatty-acyl reductase gene (pgfar) that controls the ratio of Δ11-tetradecenyl acetate stereoisomers. We identified three single nucleotide polymorphism (SNP) markers within pgfar that are differentially fixed between E- and Z-race females, and that are ≥98.2% correlated with female pheromone ratios measured by gas chromatography. Genotypic data from locations in the United States demonstrated that pgfar-z alleles were fixed within historically allopatric Z-pheromone race populations in the Midwest, and that hybrid frequency ranged from 0.00 to 0.42 within 11 sympatric sites where the two races co-occur in the Eastern United States (mean hybridization frequency or heterozygosity (H O) = 0.226 ± 0.279). Estimates of hybridization between the E- and Z-races are important for understanding the dynamics involved in maintaining race integrity, and are consistent with previous estimates of low levels of genetic divergence between E- and Z-races and the presence of weak prezygotic mating barriers. This work describes the development of new single nucleotide polymorphism (SNP) markers within the pheromone gland expressed fatty acyl reductase (pgfar) gene of Ostrinia nubilalis. These SNPs were shown to segregate based upon female pheromone production, and thus provide the first description of an assay for genetic determination of O. nubilalis pheromone strain from field-collected samples. These assays were applied to estimate hybridization within field populations, and represent valuable tools for future population genetic studies of this species.

De novo transcriptomic resources for two sibling species of moths: Ostrinia nubilalis and O. scapulalis

Background

This study aimed at enhancing the transcriptomic resources for two sibling species of moths, Ostrinia scapulalis (Adzuki bean borer) and Ostrinia nubilalis (European corn borer), as a foundation for future researches on their divergence history. Previous works on these species had shown that their genetic divergence was low, while they were reproductively isolated in natura and specialized on different host plants. Comparative genomic resources will help facilitate the understanding of the mechanisms involved in this isolation and adaptation to the host plants. Despite their fundamental interest, these species still lack the genomic resources to thoroughly identify candidate genes for functions of interest. We present here a high throughput sequencing and de novo transcriptome assembly for these two sibling species in line with this objective of comparative genomics.

Results

Based on 322,504 and 307,622 reads of 454 sequencing for O. scapulalis and O. nubilalis respectively, we reconstructed 11,231 and 10,773 transcripts, of which 40% were functionally annotated by BLAST analyzes. We determined the level of completeness of both assemblies as well as the recovery level of published Ostrinia genomic resources. Gene ontology (GO) of common and species-specific de novo transcripts did not reveal GO terms significantly enriched in one or the other species. By applying stringent homology searches on transcripts common to O. scapulalis and O. nubilalis, we identified a set of homologous transcripts, with a mean nucleotide identity value of 98.1%. In this set, the most divergent transcripts revealed candidate genes involved in developmental, sensorial and pathogen defense processes.

Conclusions

This data greatly increases the genomic resources of Ostrinia species and constitute a solid skeleton for future comparative analyzes of expression or diversity, despite we show that the transcriptomes for both species have not been assembled at full completion. In addition, we provide a set of homologous transcripts together with their annotation as a source of candidate genes for comparative analyzes.

Tissue-specific transcriptome profiling of Plutella xylostella third instar larval midgut

The larval midgut of diamondback moth, Plutella xylostella, is a dynamic tissue that interfaces with a diverse array of physiological and toxicological processes, including nutrient digestion and allocation, xenobiotic detoxification, innate and adaptive immune response, and pathogen defense. Despite its enormous agricultural importance, the genomic resources for P. xylostella are surprisingly scarce. In this study, a Bt resistant P. xylostella strain was subjected to the in-depth transcriptome analysis to identify genes and gene networks putatively involved in various physiological and toxicological processes in the P. xylostella larval midgut. Using Illumina deep sequencing, we obtained roughly 40 million reads containing approximately 3.6 gigabases of sequence data. De novo assembly generated 63,312 ESTs with an average read length of 416 bp, and approximately half of the P. xylostella sequences (45.4%, 28,768) showed similarity to the non-redundant database in GenBank with a cut-off E-value below 10(-5). Among them, 11,092 unigenes were assigned to one or multiple GO terms and 16,732 unigenes were assigned to 226 specific pathways. In-depth analysis identified genes putatively involved in insecticide resistance, nutrient digestion, and innate immune defense. Besides conventional detoxification enzymes and insecticide targets, novel genes, including 28 chymotrypsins and 53 ABC transporters, have been uncovered in the P. xylostella larval midgut transcriptome; which are potentially linked to the Bt toxicity and resistance. Furthermore, an unexpectedly high number of ESTs, including 46 serpins and 7 lysozymes, were predicted to be involved in the immune defense.As the first tissue-specific transcriptome analysis of P. xylostella, this study sheds light on the molecular understanding of insecticide resistance, especially Bt resistance in an agriculturally important insect pest, and lays the foundation for future functional genomics research. In addition, current sequencing effort greatly enriched the existing P. xylostella EST database, and makes RNAseq a viable option in the future genomic analysis.

Characterization of cDNAs encoding serine proteases and their transcriptional responses to Cry1Ab protoxin in the gut of Ostrinia nubilalis larvae

Serine proteases, such as trypsin and chymotrypsin, are the primary digestive enzymes in lepidopteran larvae, and are also involved in Bacillus thuringiensis (Bt) protoxin activation and protoxin/toxin degradation. We isolated and sequenced 34 cDNAs putatively encoding trypsins, chymotrypsins and their homologs from the European corn borer (Ostrinia nubilalis) larval gut. Our analyses of the cDNA-deduced amino acid sequences indicated that 12 were putative trypsins, 12 were putative chymotrypsins, and the remaining 10 were trypsin and chymotrypsin homologs that lack one or more conserved residues of typical trypsins and chymotrypsins. Reverse transcription PCR analysis indicated that all genes were highly expressed in gut tissues, but one group of phylogenetically-related trypsin genes, OnTry-G2, was highly expressed in larval foregut and midgut, whereas another group, OnTry-G3, was highly expressed in the midgut and hindgut. Real-time quantitative PCR analysis indicated that several trypsin genes (OnTry5 and OnTry6) were significantly up-regulated in the gut of third-instar larvae after feeding on Cry1Ab protoxin from 2 to 24 h, whereas one trypsin (OnTry2) was down-regulated at all time points. Four chymotrypsin and chymotrypsin homolog genes (OnCTP2, OnCTP5, OnCTP12 and OnCTP13) were up-regulated at least 2-fold in the gut of the larvae after feeding on Cry1Ab protoxin for 24 h. Our data represent the first in-depth study of gut transcripts encoding expanded families of protease genes in O. nubilalis larvae and demonstrate differential expression of protease genes that may be related to Cry1Ab intoxication and/or resistance.

Application of high-resolution DNA melting for genotyping in lepidopteran non-model species: Ostrinia furnacalis (Crambidae)

Development of an ideal marker system facilitates a better understanding of the genetic diversity in lepidopteran non-model organisms, which have abundant species, but relatively limited genomic resources. Single nucleotide polymorphisms (SNPs) discovered within single-copy genes have proved to be desired markers, but SNP genotyping by current techniques remain laborious and expensive. High resolution melting (HRM) curve analysis represents a simple, rapid and inexpensive genotyping method that is primarily confined to clinical and diagnostic studies. In this study, we evaluated the potential of HRM analysis for SNP genotyping in the lepidopteran non-model species Ostrinia furnacalis (Crambidae). Small amplicon and unlabeled probe assays were developed for the SNPs, which were identified in 30 females of O. furnacalis from 3 different populations by our direct sequencing. Both assays were then applied to genotype 90 unknown female DNA by prior mixing with known wild-type DNA. The genotyping results were compared with those that were obtained using bi-directional sequencing analysis. Our results demonstrated the efficiency and reliability of the HRM assays. HRM has the potential to provide simple, cost-effective genotyping assays and facilitates genotyping studies in any non-model lepidopteran species of interest.

Identification of a novel aminopeptidase P-like gene (OnAPP) possibly involved in Bt toxicity and resistance in a major corn pest (Ostrinia nubilalis)

Studies to understand the Bacillus thuringiensis (Bt) resistance mechanism in European corn borer (ECB, Ostrinia nubilalis) suggest that resistance may be due to changes in the midgut-specific Bt toxin receptor. In this study, we identified 10 aminopeptidase-like genes, which have previously been identified as putative Bt toxin receptors in other insects and examined their expression in relation to Cry1Ab toxicity and resistance. Expression analysis for the 10 aminopeptidase-like genes revealed that most of these genes were expressed predominantly in the larval midgut, but there was no difference in the expression of these genes in Cry1Ab resistant and susceptible strains. This suggested that altered expression of these genes was unlikely to be responsible for resistance in these ECB strains. However, we found that there were changes in two amino acid residues of the aminopeptidase-P like gene (OnAPP) involving Glu³⁰⁵ to Lys³⁰⁵ and Arg³⁰⁷ to Leu³⁰⁷ in the two Cry1Ab-resistant strains as compared with three Cry1Ab-susceptible strains. The mature OnAPP contains 682 amino acid residues and has a putative signal peptide at the N-terminus, a predicted glycosylphosphatidyl-inositol (GPI)-anchor signal at the C-terminal, three predicted N-glycosylation sites at residues N178, N278 and N417, and an O-glycosylation site at residue T653. We used a feeding based-RNA interference assay to examine the role of the OnAPP gene in Cry1Ab toxicity and resistance. Bioassays of Cry1Ab in larvae fed diet containing OnAPP dsRNA resulted in a 38% reduction in the transcript level of OnAPP and a 25% reduction in the susceptibility to Cry1Ab as compared with larvae fed GFP dsRNA or water. These results strongly suggest that the OnAPP gene could be involved in binding the Cry1Ab toxin in the ECB larval midgut and that mutations in this gene may be associated with Bt resistance in these two ECB strains.

Transcript analysis and comparative evaluation of shaker and slowmo gene homologues from the European corn borer, Ostrinia nubilalis

The movement and dispersal of larval Lepidoptera impact their survival and distribution within the natural landscape. Homologues of the Drosophila behaviour-linked genes shaker (shkr) and slowmo (slmo) were identified from Ostrinia nubilalis (Lepidoptera: Crambidae). Onshkr was isolated as a 1610-nucleotide (nt) constitutively expressed transcript encoding a membrane-localized 469-amino-acid (aa) protein with a conserved tetramerization domain and the six-domain architecture necessary for the molecule to fold into an active K(+) channel. Three expressed splice variants of 682, 970 and 1604 nt were identified for the Onslmo gene, and encode predicted 141 and 228 aa proteins with a conserved protein of relevant evolutionary and lymphoid interest (PRELI) domain that may function in mitochondrial protein sorting and perinuclear protein localization. Onshkr and Onslmo protein sequences aligned within monophyletic lepidopteran groups.

Transcriptome sequencing, and rapid development and application of SNP markers for the legume pod borer Maruca vitrata (Lepidoptera: Crambidae)

The legume pod borer, Maruca vitrata (Lepidoptera: Crambidae), is an insect pest species of crops grown by subsistence farmers in tropical regions of Africa. We present the de novo assembly of 3729 contigs from 454- and Sanger-derived sequencing reads for midgut, salivary, and whole adult tissues of this non-model species. Functional annotation predicted that 1320 M. vitrata protein coding genes are present, of which 631 have orthologs within the Bombyx mori gene model. A homology-based analysis assigned M. vitrata genes into a group of paralogs, but these were subsequently partitioned into putative orthologs following phylogenetic analyses. Following sequence quality filtering, a total of 1542 putative single nucleotide polymorphisms (SNPs) were predicted within M. vitrata contig assemblies. Seventy one of 1078 designed molecular genetic markers were used to screen M. vitrata samples from five collection sites in West Africa. Population substructure may be present with significant implications in the insect resistance management recommendations pertaining to the release of biological control agents or transgenic cowpea that express Bacillus thuringiensis crystal toxins. Mutation data derived from transcriptome sequencing is an expeditious and economical source for genetic markers that allow evaluation of ecological differentiation.

The application and performance of single nucleotide polymorphism markers for population genetic analyses of lepidoptera

Microsatellite markers are difficult to apply within lepidopteran studies due to the lack of locus-specific PCR amplification and the high proportion of “null” alleles, such that erroneous estimations of population genetic parameters often result. Herein single nucleotide polymorphism (SNP) markers are developed from Ostrinia nubilalis (Lepidoptera: Crambidae) using next generation expressed sequence tag (EST) data. A total of 2742 SNPs were predicted within a reference assembly of 7414 EST contigs, and a subset of 763 were incorporated into 24 multiplex PCR reactions. To validate this pipeline, 5 European and North American sample sites were genotyped at 178 SNP loci, which indicated 84 (47.2%) were in Hardy–Weinberg equilibrium. Locus-by-locus F_ST, analysis of molecular variance, and STRUCTURE analyses indicate significant genetic differentiation may exist between European and North American O. nubilalis. The observed genetic diversity was significantly lower among European sites, which may result from genetic drift, natural selection, a genetic bottleneck, or ascertainment bias due to North American origin of EST sequence data. SNPs are an abundant source of mutation data for molecular genetic marker development in non-model species, with shared ancestral SNPs showing application within closely related species. These markers offer advantages over microsatellite markers for genetic and genomic analyses of Lepidoptera, but the source of mutation data may affect the estimation of population parameters and likely need to be considered in the interpretation of empirical data.

Single nucleotide polymorphism discovery from expressed sequence tags in the waterflea Daphnia magna

BACKGROUND

Daphnia (Crustacea: Cladocera) plays a central role in standing aquatic ecosystems, has a well known ecology and is widely used in population studies and environmental risk assessments. Daphnia magna is, especially in Europe, intensively used to study stress responses of natural populations to pollutants, climate change, and antagonistic interactions with predators and parasites, which have all been demonstrated to induce micro-evolutionary and adaptive responses. Although its ecology and evolutionary biology is intensively studied, little is known on the functional genomics underpinning of phenotypic responses to environmental stressors. The aim of the present study was to find genes expressed in presence of environmental stressors, and target such genes for single nucleotide polymorphic (SNP) marker development.

RESULTS

We developed three expressed sequence tag (EST) libraries using clonal lineages of D. magna exposed to ecological stressors, namely fish predation, parasite infection and pesticide exposure. We used these newly developed ESTs and other Daphnia ESTs retrieved from NCBI GeneBank to mine for SNP markers targeting synonymous as well as non synonymous genetic variation. We validate the developed SNPs in six natural populations of D. magna distributed at regional scale.

CONCLUSIONS

A large proportion (47%) of the produced ESTs are Daphnia lineage specific genes, which are potentially involved in responses to environmental stress rather than to general cellular functions and metabolic activities, or reflect the arthropod's aquatic lifestyle. The characterization of genes expressed under stress and the validation of their SNPs for population genetic study is important for identifying ecologically responsive genes in D. magna.

Immunity in lepidopteran insects

Lepidopteran insects provide important model systems for innate immunity of insects, particularly for cell biology of hemocytes and biochemical analyses of plasma proteins. Caterpillars are also among the most serious agricultural pests, and understanding of their immune systems has potential practical significance. An early response to infection in lepidopteran larvae is the activation of hemocyte adhesion, leading to phagocytosis, nodule formation, or encapsulation. Plasmatocytes and granular cells are the hemocyte types involved in these responses. Infectious microorganisms are recognized by binding of hemolymph plasma proteins to microbial surface components. This "pattern recognition" triggers phagocytosis and nodule formation, activation of prophenoloxidase and melanization and the synthesis of antimicrobial proteins that are secreted into the hemolymph. Many hemolymph proteins that function in such innate immune responses of insects were first discovered in lepidopterans. Microbial proteinases and nucleic acids released from lysed host cells may also activate lepidopteran immune responses. Hemolymph antimicrobial peptides and proteins can reach high concentrations and may have activity against a broad spectrum of microorganisms, contributing significantly to clearing of infections. Serine proteinase cascade pathways triggered by microbial components interacting with pattern recognition proteins stimulate activation of the cytokine Spätzle, which initiates the Toll pathway for expression of antimicrobial peptides. A proteinase cascade also results inproteolytic activation of phenoloxidase and production of melanin coatings that trap and kill parasites and pathogens. The proteinases in hemolymph are regulated by specific inhibitors, including members of the serpin superfamily. New developments in lepidopteran functional genomics should lead to much more complete understanding of the immune systems of this insect group.

A novel class of miniature inverted repeat transposable elements (MITEs) that contain hitchhiking (GTCY)(n) microsatellites

The movement of miniature inverted repeat transposable elements (MITEs) modifies genome structure and function. We describe the microsatellite-associated interspersed nuclear element 2 (MINE-2), that integrates at consensus WTTTT target sites, creates dinucleotide TT target site duplications (TSDs), and forms predicted MITE-like secondary structures; a 5' subterminal inverted repeat (SIR; AGGGTTCCGTAG) that is partially complementary to a 5' inverted repeat (IR; ACGAAGCCCT) and 3'-SIRs (TTACGGAACCCT). A (GTCY)(n) microsatellite is hitchhiking downstream of conserved 5'MINE-2 secondary structures, causing flanking sequence similarity amongst mobile microsatellite loci. Transfection of insect cell lines indicates that MITE-like secondary structures are sufficient to mediate genome integration, and provides insight into the transposition mechanism used by MINE-2s.

Isolation of BAC clones containing conserved genes from libraries of three distantly related moths: a useful resource for comparative genomics of Lepidoptera

Lepidoptera, butterflies and moths, is the second largest animal order and includes numerous agricultural pests. To facilitate comparative genomics in Lepidoptera, we isolated BAC clones containing conserved and putative single-copy genes from libraries of three pests, Heliothis virescens, Ostrinia nubilalis, and Plutella xylostella, harboring the haploid chromosome number, n = 31, which are not closely related with each other or with the silkworm, Bombyx mori, (n = 28), the sequenced model lepidopteran. A total of 108–184 clones representing 101–182 conserved genes were isolated for each species. For 79 genes, clones were isolated from more than two species, which will be useful as common markers for analysis using fluorescence in situ hybridization (FISH), as well as for comparison of genome sequence among multiple species. The PCR-based clone isolation method presented here is applicable to species which lack a sequenced genome but have a significant collection of cDNA or EST sequences.

The gut transcriptome of a gall midge, Mayetiola destructor

The Hessian fly, Mayetiola destructor, is a serious pest of wheat and an experimental organism for the study of gall midge-plant interactions. In addition to food digestion and detoxification, the gut of Hessian fly larvae is also an important interface for insect-host interactions. Analysis of the genes expressed in the Hessian fly larval gut will enhance our understanding of the overall gut physiology and may also lead to the identification of critical molecules for Hessian fly-host plant interactions. Over 10,000 Expressed Sequence Tags (ESTs) were generated and assembled into 2007 clusters. The most striking feature of the Hessian fly larval transcriptome is the existence of a large number of transcripts coding for so-called small secretory proteins (SSP) with amino acids less than 250. Eleven of the 30 largest clusters were SSP transcripts with the largest cluster containing 11.3% of total ESTs. Transcripts coding for diverse digestive enzymes and detoxification proteins were also identified. Putative digestive enzymes included trypsins, chymotrypsins, cysteine proteases, aspartic protease, endo-oligopeptidase, aminopeptidases, carboxypeptidases, and alpha-amylases. Putative detoxification proteins included cytochrome P450s, glutathione S-transferases, peroxidases, ferritins, a catalase, peroxiredoxins, and others. This study represents the first global analysis of gut transcripts from a gall midge. The identification of a large number of transcripts coding for SSPs, digestive enzymes, detoxification proteins in the Hessian fly larval gut provides a foundation for future studies on the functions of these genes.

Geographic distribution of phylogenetically-distinct legume pod borer, Maruca vitrata (Lepidoptera: Pyraloidea: Crambidae)

Maruca vitrata Fabricius is a pantropical lepidopteran pest of legumes. Phylogenetic analysis of a mitochondrial cytochrome c oxidase-I gene (cox1) fragment indicates that three Maruca sp. mitochondrial lineages have unique geographic distributions [lineages 1 and 2: Australia, Taiwan, and West Africa (Niger, Nigeria, and Burkina Faso), and lineage 3: Puerto Rico]. The haplotype (T30, T114) is specific to lineages 1&2 and was assayed by NsiI and SacI polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) within population samples; it was not observed in the Puerto Rican samples, but was nearly fixed among samples from West Africa, Australia and Taiwan (85.5-100%). Re-sequencing and phylogenetic analyses of PCR-RFLP defined cox1 haplotypes indicate that nucleotide diversity is highest among samples from West Africa. Phylogenetic reconstruction based upon ribosomal DNA (rDNA) internal transcribed spacer-2 (ITS-2) sequences provided additional evidence for three Maruca sp. clades. These data suggest that multiple unique Maruca species or subspecies are present worldwide, which has implications for the management of this pest species-complex.

Pyrosequencing the Manduca sexta larval midgut transcriptome: messages for digestion, detoxification and defence

The tobacco hornworm Manduca sexta is an important model for insect physiology but genomic and transcriptomic data are currently lacking. Following a recent pyrosequencing study generating immune related expressed sequence tags (ESTs), here we use this new technology to define the M. sexta larval midgut transcriptome. We generated over 387,000 midgut ESTs, using a combination of Sanger and 454 sequencing, and classified predicted proteins into those involved in digestion, detoxification and immunity. In many cases the depth of 454 pyrosequencing coverage allowed us to define the entire cDNA sequence of a particular gene. Many new M. sexta genes are described including up to 36 new cytochrome P450s, some of which have been implicated in the metabolism of host plant-derived nicotine. New lepidopteran gene families such as the beta-fructofuranosidases, previously thought to be restricted to Bombyx mori, are also described. An unexpectedly high number of ESTs were involved in immunity, for example 39 contigs encoding serpins, and the increasingly appreciated role of the midgut in insect immunity is discussed. Similar studies of other tissues will allow for a tissue by tissue description of the M. sexta transcriptome and will form an essential complimentary step on the road to genome sequencing and annotation.

Measurements of Cry1F binding and activity of luminal gut proteases in susceptible and Cry1F resistant Ostrinia nubilalis larvae (Lepidoptera: Crambidae)

The biochemical mechanism of resistance to the Bacillus thuringiensis Cry1F toxin was studied in a laboratory-selected strain of Ostrinia nubilalis (Hübner) (Lepidoptera: Crambidae) showing more than 3000-fold resistance to Cry1F and limited cross resistance to other Cry toxins. Analyses of Cry1F binding to brush border membrane vesicles of midgut epithelia from susceptible and resistant larvae using ligand immunoblotting and Surface Plasmon Resonance (SPR) suggested that reduced binding of Cry1F to insect receptors was not associated with resistance. Additionally, no differences in activity of luminal gut proteases or altered proteolytic processing of the toxin were observed in the resistant strain. Considering these results along with previous evidence of relatively narrow spectrum of cross resistance and monogenic inheritance, the resistance mechanism in this Cry1F selected strain of O. nubilalis appears to be specific and may be distinct from previously identified resistance mechanisms reported in other Lepidoptera.

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.

Comparative performance of single nucleotide polymorphism and microsatellite markers for population genetic analysis

Microsatellite loci are standard genetic markers for population genetic analysis, whereas single nucleotide polymorphisms (SNPs) are more recent tools that require assessment of neutrality and appropriate use in population genetics. Twelve SNP markers were used to describe the genetic structure of Diabrotica virgifera virgifera (LeConte; Coleoptera: Chrysomelidae) in the United States of America and revealed a high mean observed heterozygosity (0.40 +/- 0.059) and low global F(ST) (0.029). Pairwise F(ST) estimates ranged from 0.007 to 0.045, and all but 2 populations showed significant levels of genetic differentiation (P < or = 0.008). Population parameters and conclusions based on SNP markers were analogous to that obtained by use of microsatellite markers from the identical population samples. SNP-based F(ST) estimates were 3-fold higher than corresponding estimates from microsatellites, wherein lower microsatellite F(ST) estimates likely resulted from an overestimate of migration rates between subpopulations due to convergence of allele size (homoplasy). No significant difference was observed in the proportion of SNP or microsatellite markers loci that were nonneutral within populations. SNP markers provided estimates of population genetic parameters consistent with those from microsatellite data, and their low back mutation rates may result in reduced propensity for error in estimation of population parameters.