Background: Quantitative real-time reverse transcriptase PCR (RT-qPCR) has been widely used for quantification of mRNA as a way to determine key genes involved in different biological processes. For accurate gene quantification analysis, normalization of RT-qPCR data is absolutely essential. To date, normalization is most frequently achieved by the use of internal controls, often referred to as reference genes. However, several studies have shown that the reference genes used for the quantification of mRNA expression can be affected by the experimental set-up or cell type resulting in variation of the expression level of these key genes. Therefore, the evaluation of reference genes is critical for gene expression profiling, which is often neglected in gene expression studies of insects. For this purpose, ten candidate reference genes were investigated in three different tissues (midgut, Malpighian tubules, and fat body) of the oriental fruit fly, Bactrocera dorsalis (Hendel). Results: Two different programs, geNorm and Normfinder, were used to analyze the data. According to geNorm, α-TUB + ACT5 are the most appropriate reference genes for gene expression profiling across the three different tissues in the female flies, while ACT3 + α-TUB are considered as the best for males. Furthermore, we evaluated the stability of the candidate reference genes to determine the sexual differences in the same tissue. In the midgut and Malpighian tubules, ACT2 + α-TUB are the best choice for both males and females. However, α-TUB + ACT1 are the best pair for fat body. Meanwhile, the results calculated by Normfinder are quite the same as the results with geNorm; α-TUB is always one of the most stable genes in each sample validated by the two programs. Conclusions: In this study, we validated the suitable reference genes for gene expression profiling in different tissues of B. dorsalis. Moreover, appropriate reference genes were selected out for gene expression profiling of the same tissues taking the sexual differences into consideration. This work not only formed a solid basis for future gene expression study in B. dorsalis, but also will serve as a resource to screen reference genes for gene expression studies in any other insects.

The oriental fruit fly, Bactrocera dorsalis (Hendel), is one of the most economically important pests in the world, causing serious damage to fruit production. However, lack of genetic information on this organism is an obstacle to understanding the mechanisms behind its development and its ability to resist insecticides. Analysis of the B. dorsalis transcriptome and its expression profile data is essential to extending the genetic information resources on this species, providing a shortcut that will support studies on B. dorsalis. We performed de novo assembly of a transcriptome using short read sequencing technology (Illumina). The results generated 484,628 contigs, 70,640 scaffolds, and 49,804 unigenes. Of those unigenes, 27,455 (55.13%) matched known proteins in the NCBI database, as determined by BLAST search. Clusters of orthologous groups (COG), gene orthology (GO), and the Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations were performed to better understand the functions of these unigenes. Genes related to insecticide resistance were analyzed in additional detail. Digital gene expression (DGE) libraries showed differences in gene expression profiles at different developmental stages (eggs, third-instar larvae, pupae, and adults). To confirm the DGE results, the expression profiles of six randomly selected genes were analyzed. This transcriptome greatly improves our genetic understanding of B. dorsalis and makes a huge number of gene sequences available for further study, including both genes of known importance and genes of unknown function. The DGE data provide comprehensive insight into gene expression profiles at different developmental stages. This facilitates the study of the role of each gene in the developmental process and in insecticide resistance.

Wing dimorphism is a phenomenon of phenotypic plasticity in aphid dispersal. However, the signal transduction for perceiving environmental cues (e.g., crowding) and the regulation mechanism remain elusive. Here, we found that aci-miR-9b was the only down-regulated microRNA (miRNA) in both crowding-induced wing dimorphism and during wing development in the brown citrus aphid Aphis citricidus . We determined a targeted regulatory relationship between aci-miR-9b and an ABC transporter ( AcABCG4 ). Inhibition of aci-miR-9b increased the proportion of winged offspring under normal conditions. Overexpression of aci-miR-9b resulted in decline of the proportion of winged offspring under crowding conditions. In addition, overexpression of aci-miR-9b also resulted in malformed wings during wing development. This role of aci-miR-9b mediating wing dimorphism and development was also confirmed in the pea aphid Acyrthosiphon pisum . The downstream action of aci-miR-9b- AcABCG4 was based on the interaction with the insulin and insulin-like signaling pathway. A model for aphid wing dimorphism and development was demonstrated as the following: maternal aphids experience crowding, which results in the decrease of aci-miR-9b. This is followed by the increase of ABCG4 , which then activates the insulin and insulin-like signaling pathway, thereby causing a high proportion of winged offspring. Later, the same cascade, “ miR-9b-ABCG4-insulin signaling ,” is again involved in wing development. Taken together, our results reveal that a signal transduction cascade mediates both wing dimorphism and development in aphids via miRNA. These findings would be useful in developing potential strategies for blocking the aphid dispersal and reducing viral transmission.

The oriental fruit fly, Bactrocera dorsalis (Hendel), is a serious insect pest that causes large losses to orchards globally. In this study, we conducted experi

Booklice (order Psocoptera) in the genus Liposcelis are major pests to stored grains worldwide and are closely related to parasitic lice (order Phthiraptera). We sequenced the mitochondrial (mt) genome of Liposcelis bostrychophila and found that the typical single mt chromosome of bilateral animals has fragmented into and been replaced by two medium-sized chromosomes in this booklouse; each of these chromosomes has about half of the genes of the typical mt chromosome of bilateral animals. These mt chromosomes are 8,530 bp (mt chromosome I) and 7,933 bp (mt chromosome II) in size. Intriguingly, mt chromosome I is twice as abundant as chromosome II. It appears that the selection pressure for compact mt genomes in bilateral animals favors small mt chromosomes when small mt chromosomes co-exist with the typical large mt chromosomes. Thus, small mt chromosomes may have selective advantages over large mt chromosomes in bilateral animals. Phylogenetic analyses of mt genome sequences of Psocodea (i.e. Psocoptera plus Phthiraptera) indicate that: 1) the order Psocoptera (booklice and barklice) is paraphyletic; and 2) the order Phthiraptera (the parasitic lice) is monophyletic. Within parasitic lice, however, the suborder Ischnocera is paraphyletic; this differs from the traditional view that each suborder of parasitic lice is monophyletic.