Search

Abstract Background The primary plant cell wall is a complex mixture of polysaccharides and proteins encasing living plant cells. Among these polysaccharides, cellulose is the most abundant and useful biopolymer present on earth. These polysaccharides also represent a rich source of energy for organisms which have evolved the ability to degrade them. A growing body of evidence suggests that phytophagous beetles, mainly species from the superfamilies Chrysomeloidea and Curculionoidea, possess endogenous genes encoding complex and diverse families of so-called plant cell wall degrading enzymes (PCWDEs). The presence of these genes in phytophagous beetles may have been a key element in their success as herbivores. Here, we combined a proteomics approach and transcriptome sequencing to identify PCWDEs present in larval gut contents of the mustard leaf beetle, Phaedon cochleariae. Results Using a two-dimensional proteomics approach, we recovered 11 protein bands, isolated using activity assays targeting cellulose-, pectin- and xylan-degrading enzymes. After mass spectrometry analyses, a total of 13 proteins putatively responsible for degrading plant cell wall polysaccharides were identified; these proteins belong to three glycoside hydrolase (GH) families: GH11 (xylanases), GH28 (polygalacturonases or pectinases), and GH45 (β-1,4-glucanases or cellulases). Additionally, highly stable and proteolysis-resistant host plant-derived proteins from various pathogenesis-related protein (PRs) families as well as polygalacturonase-inhibiting proteins (PGIPs) were also identified from the gut contents proteome. In parallel, transcriptome sequencing revealed the presence of at least 19 putative PCWDE transcripts encoded by the P. cochleariae genome. All of these were specifically expressed in the insect gut rather than the rest of the body, and in adults as well as larvae. The discrepancy observed in the number of putative PCWDEs between transcriptome and proteome analyses could be partially explained by differences in transcriptional level. Conclusions Combining proteome and transcriptome sequencing analyses proved to be a powerful tool for the discovery of active PCWDEs in a non-model species. Our data represent the starting point of an in-depth functional and evolutionary characterization of PCWDE gene families in phytophagous beetles and their contribution to the adaptation of these highly successful herbivores to their host plants.

Abstract Background Hormesis is a biphasic biological response characterized by the stimulatory effect at relatively low amounts of chemical compounds which are otherwise detrimental at higher concentrations. A hormetic response in larval growth rates has been observed in cotton-feeding insects in response to increasing concentrations of gossypol, a toxic metabolite found in the pigment glands of some plants in the family Malvaceae. We investigated the developmental effect of gossypol in the cotton bollworm, Helicoverpa armigera, an important heliothine pest species, by exposing larvae to different doses of this metabolite in their diet. In addition, we sought to determine the underlying transcriptional responses to different gossypol doses. Results Larval weight gain, pupal weight and larval development time were measured in feeding experiments and a hormetic response was seen for the first two characters. On the basis of net larval weight gain responses to gossypol, three concentrations (0%, 0.016% and 0.16%) were selected for transcript profiling in the gut and the rest of the body in a two-color double reference design microarray experiment. Hormesis could be observed at the transcript level, since at the low gossypol dose, genes involved in energy acquisition such as β-fructofuranosidases were up-regulated in the gut, and genes involved in cell adhesion were down-regulated in the body. Genes with products predicted to be integral to the membrane or associated with the proteasome core complex were significantly affected by the detrimental dose treatment in the body. Oxidoreductase activity-related genes were observed to be significantly altered in both tissues at the highest gossypol dose. Conclusions This study represents the first transcriptional profiling approach investigating the effects of different concentrations of gossypol in a lepidopteran species. H. armigera's transcriptional response to gossypol feeding is tissue- and dose-dependent and involves diverse detoxifying mechanisms not only to alleviate direct effects of gossypol but also indirect damage such as pH disturbance and oxygen radical formation. Genes discovered through this transcriptional approach may be additional candidates for understanding gossypol detoxification and coping with gossypol-induced stress. In a generalist herbivore that has evolved transcriptionally-regulated responses to a variety of different plant compounds, hormesis may be due to a lower induction threshold of growth-promoting, stress-coping responses and a higher induction threshold of detoxification pathways that are costly and cause collateral damage to the cell.

Abstract Background The cell suicide pathway of apoptosis is a necessary event in the life of multicellular organisms. It is involved in many biological processes ranging from development to the immune response. Evolutionarily conserved proteases, called caspases, play a central role in regulating apoptosis. Reception of death stimuli triggers the activation of initiator caspases, which in turn activate the effector caspases. In Lepidoptera, apoptosis is crucial in processes such as metamorphosis or defending against baculovirus infection. The discovery of p35, a baculovirus protein inhibiting caspase activity, has led to the characterization of the first lepidopteran caspase, Sf-Caspase-1. Studies on Sf-Caspase-1 mode of activation suggested that apoptosis in Lepidoptera requires a cascade of caspase activation, as demonstrated in many other species. Results In order to get insights into this gene family in Lepidoptera, we performed an extensive survey of lepidopteran-derived EST datasets. We identified 66 sequences distributed among 27 species encoding putative caspases. Phylogenetic analyses showed that Lepidoptera possess at least 5 caspases, for which we propose a unified nomenclature. According to homology to their Drosophila counterparts and their primary structure, we determined that Lep-Caspase-1, -2 and -3 are putative effector caspases, whereas Lep-Caspase-5 and -6 are putative initiators. The likely function of Lep-Caspase-4 remains unclear. Lep-Caspase-2 is absent from the silkworm genome and appears to be noctuid-specific, and to have arisen from a tandem duplication of the Caspase-1 gene. In the tobacco hawkmoth, 3 distinct transcripts encoding putative Caspase-4 were identified, suggesting at least 2 duplication events in this species. Conclusions The basic repertoire of five major types of caspases shared among Lepidoptera seems to be smaller than for most other groups studied to date, but gene duplication still plays a role in lineage-specific increases in diversity, just as in Diptera and mammals.

Abstract Background The intertidal zone of seacoasts, being affected by the superimposed tidal, diurnal and lunar cycles, is temporally the most complex environment on earth. Many marine organisms exhibit lunar rhythms in reproductive behaviour and some show experimental evidence of endogenous control by a circalunar clock, the molecular and genetic basis of which is unexplored. We examined the genetic control of lunar and diurnal rhythmicity in the marine midge Clunio marinus (Chironomidae, Diptera), a species for which the correct timing of adult emergence is critical in natural populations. Results We crossed two strains of Clunio marinus that differ in the timing of the diurnal and lunar rhythms of emergence. The phenotype distribution of the segregating backcross progeny indicates polygenic control of the lunar emergence rhythm. Diurnal timing of emergence is also under genetic control, and is influenced by two unlinked genes with major effects. Furthermore, the lunar and diurnal timing of emergence is correlated in the backcross generation. We show that both the lunar emergence time and its correlation to the diurnal emergence time are adaptive for the species in its natural environment. Conclusions The correlation implies that the unlinked genes affecting lunar timing and the two unlinked genes affecting diurnal timing could be the same, providing an unexpectedly close interaction of the two clocks. Alternatively, the genes could be genetically linked in a two-by-two fashion, suggesting that evolution has shaped the genetic architecture to stabilize adaptive combinations of lunar and diurnal emergence times by tightening linkage. Our results, the first on genetic control of lunar rhythms, offer a new perspective to explore their molecular clockwork.

Abstract Background The European Corn Borer, Ostrinia nubilalis (Hubner), is a keystone model for studies on the evolution of sex pheromone diversity and its role in establishing reproductive isolation. This species consists of two sympatric races, each utilizing opposite isomers of the same compound as their major pheromone component. Female production and male response are congruent in each race, and males from each strain exhibit phenotypic differences in peripheral physiology. Both strains possess co-localized pheromone-sensitive olfactory sensory neurons characterized by a larger amplitude action potential (spike) responding to the major pheromone component, and a smaller spike amplitude cell responding to the minor component, i.e. the opposite isomer. These differences in amplitude correspond to differences in dendritic diameter between the two neurons. Previous studies showed that behavioral response to the pheromone blend was sex-linked, but spike amplitude response to pheromone components matched autosomal, not sex-linked inheritance. Results As part of a larger study to finely map the loci responsible for pheromone communication in this species, we have reanalyzed peripheral physiology among parental, and first and second generation hybrids between the two pheromone strains using tungsten electrode electrophysiology. Our results reveal that differences in spike amplitude ratio between male pheromone-sensitive sensory neurons in O. nubilalis races are controlled, at least partially, by sex-linked genes that exhibit E-strain dominance. Conclusions We propose that peripheral olfactory response in O. nubilalis may be affected both by autosomal and sex-linked genes exhibiting a cross-locus dominance effect, and suggest that the genetic basis for pheromone reception and response in the species is more closely linked than previously thought.

Abstract Background Developing lepidopteran microsatellite DNA markers can be problematical, as markers often exhibit multiple banding patterns and high frequencies of non-amplifying "null" alleles. Previous studies identified sequences flanking simple sequence repeat (SSR) units that are shared among many lepidopteran species and can be grouped into microsatellite-associated DNA families. These families are thought to be associated with unequal crossing-over during DNA recombination or with transposable elements (TEs). Results We identified full-length lepidopteran non-LTR retrotransposable elements of the RTE clade in Heliconius melpomene and Bombyx mori. These retroelements possess a single open reading frame encoding the Exonuclease/Endonuclease/Phosphatase and the Reverse Transcriptase/nLTR domains, a 5' UTR (untranslated region), and an extremely short 3' UTR that regularly consists of SSR units. Phylogenetic analysis supported previous suggestions of horizontal transfer among unrelated groups of organisms, but the diversity of lepidopteran RTE elements appears due to ancient divergence of ancestral elements rather than introgression by horizontal transfer. Similarity searches of lepidopteran genomic sequences in GenBank identified partial RTE elements, usually consisting of the 3' terminal region, in 29 species. Furthermore, we identified the C-terminal end of the Reverse Transcriptase/nLTR domain and the associated 3' UTR in over 190 microsatellite markers from 22 lepidopteran species, accounting for 10% of the lepidopteran microsatellites in GenBank. Occasional retrotransposition of autonomous elements, frequent retrotransposition of 3' partial elements, and DNA replication slippage during retrotransposition offers a mechanistic explanation for the association of SSRs with RTE elements in lepidopteran genomes. Conclusions Non-LTR retrotransposable elements of the RTE clade therefore join a diverse group of TEs as progenitors of SSR units in various organisms. When microsatellites are isolated using standard SSR enrichment protocols and primers designed at complementary repeated regions, amplification from multiple genomic sites can cause scoring difficulties that compromise their utility as markers. Screening against RTE elements in the isolation procedure provides one strategy for minimizing this problem.

Abstract Background A molecular population genetics understanding is central to the study of ecological and evolutionary functional genomics. Population genetics identifies genetic variation and its distribution within and among populations, it reveals the demographic history of the populations studied, and can provide indirect insights into historical selection dynamics. Here we use this approach to examine the demographic and selective dynamics acting of a candidate gene involved in plant-insect interactions. Previous work documents the macroevolutionary and historical ecological importance of the nitrile-specifier protein (Nsp), which facilitated the host shift of Pieridae butterflies onto Brassicales host plants ~80 Myr ago. Results Here we assess the microevolutionary dynamics of the Nsp gene by studying the within and among-population variation at Nsp and reference genes in the butterfly Pieris rapae (Small Cabbage White). Nsp exhibits unexpectedly high amounts of amino acid polymorphism, unequally distributed across the gene. The vast majority of genetic variation exists within populations, with little to no genetic differentiation among four populations on two continents. A comparison of synonymous and nonsynonymous substitutions in 70 randomly chosen genes among P. rapae and its close relative Pieris brassicae (Large Cabbage White) finds Nsp to have a significantly relaxed functional constraint compared to housekeeping genes. We find strong evidence for a recent population expansion and no role for strong purifying or directional selection upon the Nsp gene. Conclusions The microevolutionary dynamics of the Nsp gene in P. rapae are dominated by recent population expansion and variation in functional constraint across the repeated domains of the Nsp gene. While the high amounts of amino acid diversity suggest there may be significant functional differences among allelic variants segregating within populations, indirect tests of selection could not conclusively identify a signature of historical selection. The importance of using this information for planning future studies of potential performance and fitness consequences of the observed variation is discussed.

Abstract Background The chemical components of sex pheromones have been determined for more than a thousand moth species, but so far only a handful of genes encoding enzymes responsible for the biosynthesis of these compounds have been identified. For understanding the evolution of moth sexual communication, it is essential to know which genes are involved in the production of specific pheromone components and what controls the variation in their relative frequencies in the pheromone blend. We used a transcriptomic approach to characterize the pheromone gland of the Noctuid moth Heliothis virescens, an important agricultural pest, in order to obtain substantial general sequence information and to identify a range of candidate genes involved in the pheromone biosynthetic pathway. Results To facilitate identifying sets of genes involved in a broad range of processes and to capture rare transcripts, we developed our majority of ESTs from a normalized cDNA library of Heliothis virescens pheromone glands (PG). Combining these with a non-normalized library yielded a total of 17,233 ESTs, which assembled into 2,082 contigs and 6,228 singletons. Using BLAST searches of the NR and Swissprot databases we were able to identify a large number of putative unique gene elements (unigenes), which we compared to those derived from previous transcriptomic surveys of the larval stage of Heliothis virescens. The distribution of unigenes among GO Biological Process functional groups shows an overall similarity between PG and larval transcriptomes, but with distinct enrichment of specific pathways in the PG. In addition, we identified a large number of candidate genes in the pheromone biosynthetic pathways. Conclusion These data constitute one of the first large-scale EST-projects for Noctuidae, a much-needed resource for exploring these pest species. Our analysis shows a surprisingly complex transcriptome and we identified a large number of potential pheromone biosynthetic pathway and immune-related genes that can be applied to population and systematic studies of Heliothis virescens and other Noctuidae.

Abstract 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.

Abstract Background Poly- and oligophagous insects are able to feed on various host plants with a wide range of defense strategies. However, diverse food plants are also inhabited by microbiota differing in quality and quantity, posing a potential challenge for immune system mediated homeostasis in the herbivore. Recent studies highlight the complex interactions between environmentally encountered microorganisms and herbivorous insects, pointing to a potential adaptational alteration of the insects' physiology. We performed a differential gene expression analysis in whole larvae and eggs laid by parents grown on different diets to identify potential novel genes related to elevated microbial content in the caterpillars' food. Results We used GeneFishing, a novel differential display method, to study the effects of dietary bacteria on the general gene expression in different life stages and tissues of the cabbage looper (Trichoplusia ni). We were able to visualize several hundred transcripts on agarose gels, one fifth of which were differentially expressed between treatments. The largest number of differentially expressed genes was found in defense-related processes (13) and in recognition and metabolism (16). 21 genes were picked out and further tested for differential gene expression by an independent method (qRT-PCR) in various tissues of larvae grown on bacterial and bacteria-free diet, and also in adults. We detected a number of genes indicative of an altered physiological status of the insect, depending on the diet, developmental stage and tissue. Conclusion Changes in immune status are accompanied by specific changes in the transcript levels of genes connected to metabolism and homeostasis of the organism. Our findings show that larval feeding on bacteria-rich diet leads to substantial gene expression changes, potentially resulting in a reorganization of the insects' metabolism to maintain organismal homeostasis, not only in the larval but also in the adult stage. Furthermore, differences in gene expression levels can also be seen in the next generation, strongly influenced by parental diet.

Abstract Background The fall armyworm Spodoptera frugiperda (Lepidoptera; Noctuidae) consists of two distinct strains with different host plant preferences for corn and rice. To assess whether pheromonal-mediated behavioral isolation accompanies the habitat isolation on different host plants, we compared the sex pheromone composition among females of the two strains. Pheromone glands were extracted with or without injection of pheromone biosynthesis activating neuropeptide (PBAN). To assess the mode of inheritance of this variation, we also analyzed the pheromone composition of F1 hybrid females. Results Relative to intra-strain variation, the pheromone composition of the two strains differed significantly. Corn strain females contained significantly more of the second most abundant pheromone compound Z11-16:Ac (m), and significantly less of most other compounds, than rice strain females. When females were injected with PBAN before their glands were extracted, the differences between the strains were less pronounced but still statistically significant. The pheromone composition of hybrid females showed a maternal inheritance of the major component Z9-14:Ac (M) as well as of Z11-16:Ac (m). Most other compounds showed an inheritance indicating genetic dominance of the corn strain. The within-strain phenotypic correlations among the various components were consistent with their hypothesized biosynthetic pathway, and between-strain differences in the correlation structure suggested candidate genes that may explain the pheromone differences between the two strains. These include Δ9- and Δ11 desaturases, and possibly also a Δ7-desaturase, although the latter has not been identified in insects so far. Conclusion The two host strains of S. frugiperda produce systematically differing female sex pheromone blends. Previously-documented geographic variation in the sexual communication of this species did not take strain identity into account, and thus may be partly explained by different strain occurrence in different regions. The finding of pheromone differences reinforces the possibility of incipient reproductive isolation among these strains, previously shown to differ in the timing of nocturnal mating activity and host plant use. Finding the genetic basis of the pheromone differences, as well as these other biological traits, will help to elucidate the role of premating isolation in the continuing differentiation of these two strains that may eventually lead to speciation.

Abstract Background Insects helped pioneer, and persist as model organisms for, the study of specific aspects of immunity. Although they lack an adaptive immune system, insects possess an innate immune system that recognizes and destroys intruding microorganisms. Its operation under natural conditions has not been well studied, as most studies have introduced microbes to laboratory-reared insects via artificial mechanical wounding. One of the most common routes of natural exposure and infection, however, is via food; thus, the role of dietary microbial communities in herbivorous insect immune system evolution invites study. Here, we examine the immune system response and consequences of exposing a lepidopteran agricultural pest to non-infectious microorganisms via simple oral consumption. Results Immune system response was compared between Trichoplusia ni larvae reared on diets with or without non-pathogenic bacteria (Escherichia coli and Micrococcus luteus). Two major immune response-related enzymatic activities responded to diets differently – phenoloxidase activity was inhibited in the bacteria-fed larvae, whereas general antibacterial activity was enhanced. Eight proteins were highly expressed in the hemolymph of the bacteria fed larvae, among them immune response related proteins arylphorin, apolipophorin III and gloverin. Expression response among 25 putative immune response-related genes were assayed via RT-qPCR. Seven showed more than fivefold up regulation in the presence of bacterial diet, with 22 in total being differentially expressed, among them apolipophorin III, cecropin, gallerimycin, gloverin, lysozyme, and phenoloxidase inhibiting enzyme. Finally, potential life-history trade-offs were studied, with pupation time and pupal mass being negatively affected in bacteria fed larvae. Conclusion The presence of bacteria in food, even if non-pathogenic, can trigger an immune response cascade with life history tradeoffs. Trichoplusia ni larvae are able to detect and respond to environmental microbes encountered in the diet, possibly even using midgut epithelial tissue as a sensing organ. Potential benefits of this immune system priming may outweigh the observed tradeoffs, as priming based on environmentally sensed bacterial may decrease risk of serious infection. These results show that food plant microbial communities represent a dynamic and unstudied part of the coevolutionary interactions between plants and their insect herbivores.

Abstract Background Helicoverpa armigera and H. zea are amongst the most significant polyphagous pest lepidopteran species in the Old and New Worlds respectively. Separation of H. armigera and H. zea is difficult and is usually only achieved through morphological differences in the genitalia. They are capable of interbreeding to produce fertile offspring. The single species status of H. armigera has been doubted, due to its wide distribution and plant host range across the Old World. This study explores the global genetic diversity of H. armigera and its evolutionary relationship to H zea. Results We obtained partial (511 bp) mitochondrial DNA (mtDNA) Cytochrome Oxidase-I (COI) sequences for 249 individuals of H. armigera sampled from Australia, Burkina Faso, Uganda, China, India and Pakistan which were associated with various host plants. Single nucleotide polymorphisms (SNPs) within the partial COI gene differentiated H. armigera populations into 33 mtDNA haplotypes. Shared haplotypes between continents, low F-statistic values and low nucleotide diversity between countries (0.0017 – 0.0038) suggests high mobility in this pest. Phylogenetic analysis of four major Helicoverpa pest species indicates that H. punctigera is basal to H. assulta, which is in turn basal to H. armigera and H. zea. Samples from North and South America suggest that H. zea is also a single species across its distribution. Our data reveal short genetic distances between H. armigera and H. zea which seem to have been established via a founder event from H. armigera stock at around 1.5 million years ago. Conclusion Our mitochondrial DNA sequence data supports the single species status of H. armigera across Africa, Asia and Australia. The evidence for inter-continental gene flow observed in this study is consistent with published evidence of the capacity of this species to migrate over long distances. The finding of high genetic similarity between Old World H. armigera and New World H. zea emphasises the need to consider work on both pests when building pest management strategies for either.