The small brown planthopper (SBPH), Laodelphax striatellus, is one of the most serious pest insects of rice plants.Buprofezin has been used to control SBPH for more than a decade, however, the occurrence of buprofezin resistant SBPHwas reported recently. To develop an alternative pest control an alternative pest control strategy, RNA-seq of buprofezin-treatedSBPH was performed to screen the insecticidal target genes for RNA interference (RNAi). Six genes were selected fordsRNA synthesis, and applied to SBPH to assess the insecticidal efficacy. Two and three of those dsRNAs showed moderatedand substantial insecticidal activity up to 60% of mortality in one week, respectively. These results demonstrated the potentialof gene screening strategy for the development of RNAi-based pest management program.

To identify genes that commonly respond to the treatment of different insecticides and are responsible for the toleranceenhancement, transcriptomic profiles of larvae treated with sublethal doses of the five insecticides were compared withthat of untreated control. A total of 117,181 transcripts with a mean length of 662 bp were generated by de novo assembly,of which 35,329 transcripts were annotated. Among them, 125, 143, 182, 215 and 149 transcripts were determined tobe up-regulated whereas 67, 45, 60, 60 and 38 genes were down-regulated following treatments with these five insecticides.The most notable examples of commonly responding over-transcribed genes were two cytochrome P450 genes and ninecuticular protein genes. In contrast, several genes composing the mitochondrial energy generation system were significantlydown-regulated in all treated larvae. Considering the distinct structure and mode of action of the five insecticides tested,the differentially expressed genes identified in this study appear to be involved in general chemical defense at the initialstage of intoxication. Their possible roles in the tolerance/resistance development were discussed.

Worldwide, increasing numbers insecticide resistant insect is one of the main problem in agriculture not only in the field but also in the storage. The rusty grain beetle, Cryptolestes ferrugineus is one of the cosmopolitan insect that infests a wide range of stored cereals and related commodities. Until quite recently, phosphine (PH3) has been effective in controlling this species in worldwide including Korea. However, strongly resistant populations of RGB have been detected in Australia that could threaten market access of infested commodities. Resistant populations detected in Australia showed extremely high levels of resistance to phosphine, up to 1300 folds higher than that of susceptible strain. So here we tried to identify their phosphine resistance mechanism based on transcriptome analysis using RNaseq in adult stage. Over 10Gb were sequenced in each strains and some of specific P450 were over expressed in resistance strain.

The pheromone biosynthesis in Plutella xylostella is more active in the scotophase than in the photophase, indicating that there may be changes of gene expression in the pheromone glands. To identify genes contributing to change in pheromone production, we analyzed transcriptomes of pheromone glands from both decapitated females (PG-minus) in the photophase and normal ones (PG-plus) in the scotophase. Deep sequencing for mRNAs in the pheromone gland yielded approximately 7.5Gb and 6,671 transcripts showing positive FPKM value were analyzed. Differentially expressed gene analysis revealed that up- and down-regulated transcripts were 310 and 326 in the PG-plus transcriptome, respectively. Genes putatively involved in the pheromone biosynthesis pathway were identified such as acetyl-CoA carboxylase, acetyl-CoA dehydrogenase, fatty acid synthase (FAS), desaturases (Δ9 and Δ 11) and fatty acid reductases of pheromone gland (pgFAR), alcohol oxidase, aldehyde oxidase and aldehyde reductase, etc. Quantitative RT-PCR revealed that expressions of FAS, Δ11 desaturase and pgFAR were significantly higher in PG-plus, suggesting that they may have crucial roles in sex pheromone biosynthesis of P. xylostella

The aminoacyl-tRNA synthetases (ARSs) are ancient house-keeping enzymes that catalyze the ligation of tRNAs to their cognate amino acids in the first step of protein synthesis. During the evolution of higher eukaryotes, cytoplasmic ARSs have undergone significant changes including the addition of new domains that are not part of the enzymatic core. These additional regions have been found to be associated with a broad range of biological functions beyond protein synthesis. The non-translational functions of ARSs appear to be regulated by their presence within a cytoplasmic multi-tRNA synthetase complex (MSC), which is assembled through the appended domains. We recently reported that the MSC member glutamylprolyl- tRNA synthetase (EPRS) promotes antiviral gene expression through its infection-specific phosphorylation and release from the MSC. Here, we conducted transcriptome analysis of influenza A virusinfected cells. We particularly focused on the analysis of chemokine-related gene expression, in combination with chemokine array analysis against virus infection. Moreover, the correlation between chemokine expression pattern and EPRS function in response to different stimuli was assessed. The results showed that viral infection increases interferon-response and pro-inflammatory chemokine expression. In contrast, the level of chemokine expression was suppressed in interferon-γ treated cells. Thus, these results further demonstrate the previously reported stimulus-specific EPRS functions in immune responses.

Varroa destructor is a devastating ectoparasitic mite which attacks Honeybee, Apis mellifera. V. destructor feeds on honeybee hemolymph, and often harbors small RNA viruses such as the deformed wing virus to transmit these viruses in the infested bee hive. To survey the genes of V. destructor, total RNA was subjected to high-throughput transcriptome sequencing to construct in silico cDNA library by using the Illumina HiSeq 2000 platform. Total of 2×107,748,792 paired-end short reads were obtained and quality filtered reads were subjected to Trinity de novo assembler followed by TransDecoder, and CD-HIT program to make a V. destructor reference cDNA library containing 28,023 of clustered contigs with protein coding capacity. These cDNA sequences will help us to understand the molecular biology of V. destructor.

A cDNA isolated from female adult heads of Maruca vitrata encodes 197 amino acids including PBAN. Synthetic Mav-PBAN induced pheromone production in the pheromone gland, indicating that this synthetic peptide was biologically functional. Expression of Mav-PBAN cDNA was found in all examined body parts whereas PBAN receptor only in the pheromone gland. Transcriptomic analysis revealed that 191 contigs involved in the pheromone biosynthesis such as PBAN receptor, PBAN, fatty acid transport proteins, acetyl-CoA carboxylases, fatty acid synthases, desaturases (FAD), β-oxidation enzymes, and fatty acyl-CoA reductases (FARs) were identified.

A subtropical insect species, Maruca vitrata, is invasive to temperate zones. Supercooling capacity of M. vitrata was confirmed in all developmental stages, in which egg exhibited the lowest supercooling point (SCP) at –22.5oC. However, low temperature lethality was observed at much higher temperatures than SCPs in all developmental stages. In addition, this nonfreezing injury increased with incubation time at rearing condition. Uncontrolled phenoloxidase activation appeared to cause the nonfreezing injury. Preexposure to a cool temperature before lethal low temperature significantly increased survival of M. vitrata. This cold hardening accompanied with increase of glycerol content in hemolymph. Transcriptome of M. vitrata under cond-hardening was assessed.

Beauveria bassiana (Bb) is an entomopathogenic fungus with a wide host range, and is commonly used as an environment-friendly biopesticide. However, the molecular mechanisms of Bb-host interactions are not well understood. Here, RNA isolated from a highly virulent strain of B. bassiana (Bb JEF-007) and Riptortus pedestris (Hemiptera: Alydidae) (bean bug) infected with this strain were subjected to high throughput next generation sequencing (NGS) to analyze and compare transcriptomes. Differentially expressed gene (DEG) analysis showed that 2,381 genes were up-regulated and 2,303 genes were down-regulated upon infection. Most DEGs were classified into the categories of single-organism, cellular and metabolism processes by gene ontology (GO) analysis. Carbon metabolism-related enzymes in the glyoxylate cycle were significantly up-regulated, suggesting a possible role for them in Bb growth in the host. This work provides insight into how entomopathogenic B. bassiana occupies agriculturally harmful bean bug at the late stage, which might be essential during fungal infection.

The small brown planthopper (SBPH), Laodelphax striatellus Fallén (Hemiptera: Delphacidae) is one of the major insect pest against rice, Oryza sativa L. in Korea. High density of SBPH could cause severe damage on rice plant by directly sucking and indirectly transmitting viral pathogens, Rice stripe virus and Rice streaked dwarf virus. As a preliminary study for de novo whole-genome sequencing of SBPH, we investigated 6 transcriptomes isolated from different developmental stages, sex, and tissue (egg, 1st ~ 3rd nymphs, 4th ~ 5th nymphs, female and male adults, salivary gland). Clean-sequence data of 19.3 Gb were obtained from total 47.8 Gb raw data after adaptor and quality trimming (Q30) and overlapped reads joining. As a suitable assembler, Bridger was selected based on the results of reference mapping (93.45%) and CEGMA completeness (95.97%). Finally, we obtained 158,207 reads (size range: 201 ~ 22,162 bp; Mean size: 1,048.04 bp; N50: 2,417 bp) after clustering the assembly results by CD-HIT-EST (similarity threshold: 99%). Based on these results, we are conducting further studies such as transcript expression pattern among different developmental stages and gene annotation.

Tropilaelaps mercedesae is an ectoparasite of immature honey bees belonging to the genus Tropilaelaps (Acari: Laelapidae). T. mercedesae has become a major threat to the Western honey bee Apis mellifera in Asia, including Korea, and is expanding its geographical range to northern regions due to global warming. To establish gene resources of T. mercedesae, the whole transcriptome was analyzed by RNA sequencing. An mRNA-focused library was generated from total RNA extracted from the mixed stages using the TruSeq RNA Library Preparation kit and sequenced using the HiSeq 2000 platform. A total of 6.0 Gb reads were obtained with 85% Q30 value. Trimmed sequence data were de novo assembled using the CLC Assembly Cell v 4.2. A total of 64,868 non-duplicate contigs were finally obtained and annotated by the Blast2GO using the NCBI nr database. The most abundant species in the resulting 14,336 Blast hits (22.1%) was Metaseiulus occidentalis, a predatory mite, followed by Ixodes scapularis and Tribolium castaneum, suggesting that the T. mercedesae transcriptome matches well with closely related other arthropod species, including mites and ticks. In order to provide basic information for efficient control and monitoring of potential resistance in T. mercedesae, acaricide target genes were annotated and characterized. One voltage-sensitive sodium channel gene encoding the molecular target of fluvalinate, a pyrethroid acaricide most widely used for the control of T. mercedesae, was identified and its molecular properties were investigated. In addition, other acaricide target genes, including acetylcholinesterase and glutamate (or GABA)-gated chloride channel, were identified and characterized.

The ascomycete fungus Beauveria bassiana is a wide host range entomo- pathogenic fungus, which is commonly used as an environmental friendly biopesticide. However, the molecular mechanisms of host-pathogen interaction of B. bassiana are not well understood. Here, the high throughput next generation sequencing was performed to analyze the transcriptome of B. bassiana JEF-007 infected bean bug (Riptorus pedestris). Differentially expressed gene (DEG) analysis results showed that total 4,684 genes including 2,381 up and 2,303 down regulated genes were identified. Most of the DEGs were classified into single- organism, cellular and metabolism processes by gene ontology (GO) analysis. Metabolism pathway was the most abound category of DEGs via KEGG pathway mapping. Several possible candidates of virulence factors were dramatically expressed after infection, such as cytotoxic lectin, bacterial-like toxin, and proteins related to cell wall, hyphal growth, nutrient uptake and halogenated compounds synthesis. Furthermore, we also found the highest expression of a novel small RNA virus in the infected bean bug, but the relationship between fungal virulence and the RNA virus was under determination. The functional roles of these possible virulence factors are remained unclear, but this work provides a new insight for further fungal studies. Our results reflect systemic impacts of fungal pathogenesis and these findings represent a significant advance in the fungal functional genomics.

The pheromone biosynthesis in Plutella xylostella is more active in the scotophase than in the photophase, which suggests that there may be changes of gene expression in the pheromone glands. To identify genes contributing to change in pheromone production, we analyzed transcriptomes from pheromone glands of both decapitated females in the photophase and normal ones in the scotophase. Comparative analysis were performed with transcriptomes of pheromone glands from non-decapitated (PG) females and decapitated ones for identification and expression of putative genes associated with pheromone biosynthesis pathway. Deep sequencing for mRNAs in the pheromone gland yielded approximately 7.5Gb and totally 17265 transcript were constructed under a homology cutoff of 10-6 Evalue. Genes putatively involved in pheromone biosynthesis were identified such as acetyl-CoA carboxylase, acetyl-CoAdehydrogenase, fatty acid synthase (FAS), desaturases (Δ9 and Δ 11) and fatty acid reductases (FAR) including pgFAR, alcohol oxidase, aldehyde oxidase and aldehyde reductase, etc. Expression of 6 signal genes involving in pheromone biosynthesis such as acyl-CoA desaturase, FAR, PBAN receptor, fatty acid transporter, acyl-CoA binding protein did not exhibited ant significant different in both transcriptomes. Quantitative RT-PCR revealed that expressions of FAS, Δ11 desaturase and pgFAR were higher in PG than that in ΔPG. Based on results, Δ11 desaturase and pgFAR may have a crucial role in sex pheromone biosynthesis of P. xylostella.