Hemiptera, the largest non-holometabolous order of insects, represents approximately 7% of metazoan diversity. With extraordinary life histories and highly specialized morphological adaptations, hemipterans have exploited diverse habitats and food sources through approximately 300 Myr of evolution. To elucidate the phylogeny and evolutionary history of Hemiptera, we carried out the most comprehensive mitogenomics analysis on the richest taxon sampling to date covering all the suborders and infraorders, including 34 newly sequenced and 94 published mitogenomes. With optimized branch length and sequence heterogeneity, Bayesian analyses using a site-heterogeneous mixture model resolved the higher-level hemipteran phylogeny as (Sternorrhyncha, (Auchenorrhyncha, (Coleorrhyncha, Heteroptera))). Ancestral character state reconstruction and divergence time estimation suggest that the success of true bugs (Heteroptera) is probably due to angiosperm coevolution, but key adaptive innovations (e.g. prognathous mouthpart, predatory behaviour, and haemelytron) facilitated multiple independent shifts among diverse feeding habits and multiple independent colonizations of aquatic habitats.

Quantitative real-time PCR (qRT-PCR), a primary tool in gene expression analysis, requires an appropriate normalization strategy to control for variation among samples. The best option is to compare the mRNA level of a target gene with that of reference gene(s) whose expression level is stable across various experimental conditions. In this study, expression profiles of eight candidate reference genes from the diamondback moth, Plutella xylostella, were evaluated under diverse experimental conditions. RefFinder, a web-based analysis tool, integrates four major computational programs including geNorm, Normfinder, BestKeeper, and the comparative ΔCt method to comprehensively rank the tested candidate genes. Elongation factor 1 (EF1) was the most suited reference gene for the biotic factors (development stage, tissue, and strain). In contrast, although appropriate reference gene(s) do exist for several abiotic factors (temperature, photoperiod, insecticide, and mechanical injury), we were not able to identify a single universal reference gene. Nevertheless, a suite of candidate reference genes were specifically recommended for selected experimental conditions. Our finding is the first step toward establishing a standardized qRT-PCR analysis of this agriculturally important insect pest.

Accurate evaluation of gene expression requires normalization relative to the expression of reliable reference genes. Expression levels of “classical” reference genes can differ, however, across experimental conditions. Although quantitative real-time PCR (qRT-PCR) has been used extensively to decipher gene function in the sweetpotato whitefly Bemisia tabaci, a world-wide pest in many agricultural systems, the stability of its reference genes has rarely been validated. In this study, 15 candidate reference genes from B. tabaci were evaluated using two Excel-based algorithms geNorm and Normfinder under a diverse set of biotic and abiotic conditions. At least two reference genes were selected to normalize gene expressions in B. tabaci under experimental conditions. Specifically, for biotic conditions including host plant, acquisition of a plant virus, developmental stage, tissue (body region of the adult), and whitefly biotype, ribosomal protein L29 was the most stable reference gene. In contrast, the expression of elongation factor 1 alpha, peptidylprolyl isomerase A, NADH dehydrogenase, succinate dehydrogenase complex subunit A and heat shock protein 40 were consistently stable across various abiotic conditions including photoperiod, temperature, and insecticide susceptibility. Our finding is the first step toward establishing a standardized quantitative real-time PCR procedure following the MIQE (Minimum Information for publication of Quantitative real time PCR Experiments) guideline in an agriculturally important insect pest, and provides a solid foundation for future RNA interference based functional study in B. tabaci.

Background: Termite lignocellulose digestion is achieved through a collaboration of host plus prokaryotic and eukaryotic symbionts. In the present work, we took a combined host and symbiont metatranscriptomic approach for investigating the digestive contributions of host and symbiont in the lower termiteReticulitermes flavipes. Our approach consisted of parallel high-throughput sequencing from (i) a host gut cDNA library and (ii) a hindgut symbiont cDNA library. Subsequently, we undertook functional analyses of newly identified phenoloxidases with potential importance as pretreatment enzymes in industrial lignocellulose processing.Results: Over 10,000 expressed sequence tags (ESTs) were sequenced from the 2 libraries that aligned into 6,555 putative transcripts, including 171 putative lignocellulase genes. Sequence analyses provided insights in two areas. First, a non-overlapping complement of host and symbiont (prokaryotic plus protist) glycohydrolase gene families known to participate in cellulose, hemicellulose, alpha carbohydrate, and chitin degradation were identified. Of these, cellulases are contributed by host plus symbiont genomes, whereas hemicellulases are contributed exclusively by symbiont genomes. Second, a diverse complement of previously unknown genes that encode proteins with homology to lignase, antioxidant, and detoxification enzymes were identified exclusively from the host library (laccase, catalase, peroxidase, superoxide dismutase, carboxylesterase, cytochrome P450). Subsequently, functional analyses of phenoloxidase activity provided results that were strongly consistent with patterns of laccase gene expression. In particular, phenoloxidase activity and laccase gene expression are mostly restricted to symbiont-free foregut plus salivary gland tissues, and phenoloxidase activity is inducible by lignin feeding.Conclusion: To our knowledge, this is the first time that a dual host-symbiont transcriptome sequencing effort has been conducted in a single termite species. This sequence database represents an important new genomic resource for use in further studies of collaborative host-symbiont termite digestion, as well as development of coevolved host and symbiont-derived biocatalysts for use in industrial biomass-to-bioethanol applications. Additionally, this study demonstrates that: (i) phenoloxidase activities are prominent in theR. flavipesgut and are not symbiont derived, (ii) expands the known number of host and symbiont glycosyl hydrolase families inReticulitermes, and (iii) supports previous models of lignin degradation and host-symbiont collaboration in cellulose/hemicellulose digestion in the termite gut. All sequences in this paper are available publicly with the accession numbers FL634956-FL640828 (Termite Gut library) and FL641015-FL645753 (Symbiont library).

For many insect-vectored plant viruses, the relationship between feeding behavior and vector competence may prove integral to an understanding of the epidemiology of the resulting plant disease. While plant-infecting viruses are well known to change host plant physiology in a way that makes them more attractive to vectors, viral manipulation of the vectors themselves has only recently been reported. Previous research suggested that the rapid spread of Tomato yellow leaf curl virus (TYLCV) throughout China has been facilitated by its primary vector, the whitefly Bemisia tabaci . We conducted two experiments testing the impact of TYLCV infection of the host plant (tomato) and vector ( B. tabaci biotypes B and Q) on whitefly feeding behavior. Whiteflies of biotypes B and Q both appeared to find TYLCV-infected plants more attractive, probing them more quickly and having a greater number of feeding bouts; this did not, however, alter the total time spent feeding. Viruliferous whiteflies fed more readily than uninfected whiteflies and spent more time salivating into sieve tube elements. Because vector salivation is essential for viral transmission, this virally mediated alteration of behavior should provide TYLCV a direct fitness benefit. This is the first report of such manipulation by a nonpropagative virus that belongs to an exclusively plant-infecting family of viruses ( Geminiviridae ). In the context of previous research showing that feeding on TYLCV-infected plants harms biotype B but helps biotype Q, the fact that both biotypes were equally affected by TYLCV also suggests that the virus may alter the biotype B-biotype Q competitive interaction in favor of biotype Q.

The arms race between entomopathogenic bacteria and their insect hosts is an excellent model for decoding the intricate coevolutionary processes of host-pathogen interaction. Here, we demonstrate that the MAPK signaling pathway is a general switch to trans-regulate differential expression of aminopeptidase N and other midgut genes in an insect host, diamondback moth (Plutella xylostella), thereby countering the virulence effect of Bacillus thuringiensis (Bt) toxins. Moreover, the MAPK cascade is activated and fine-tuned by the crosstalk between two major insect hormones, 20-hydroxyecdysone (20E) and juvenile hormone (JH) to elicit an important physiological response (i.e. Bt resistance) without incurring the significant fitness costs often associated with pathogen resistance. Hormones are well known to orchestrate physiological trade-offs in a wide variety of organisms, and our work decodes a hitherto undescribed function of these classic hormones and suggests that hormonal signaling plasticity is a general cross-kingdom strategy to fend off pathogens.