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.

Pheromone biosynthesis-activating neuropeptide (PBAN), produced in subesophageal ganglion, is known to stimulate pheromone biosynthesis in Plutella xylostella. The pheromone production 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 analyze gene expression involving in pheromone biosynthesis for 24 hrs, we performed transcriptomes of pheromone glands which were isolated at every 4 h. Eleven pheromone biosynthesis-related genes, acetyl-CoA carboxylase (ACC), fatty acid synthase (FAS), acyl-CoA dehydrogenase, enoyl-CoA hydratase, 3-hydroxyacyl-CoA dehydrogenase, Δ9 desaturase, Δ11 desaturase, fatty acid reductase (FAR), alcohol oxidase, aldehyde oxidase, and aldehyde reductase were identified. Among these genes, ACC, FAS and chain shortening enzymes involving in early stage of pheromone biosynthesis exhibited their highest expression between AM9 and PM5. Desaturases, Δ9 and Δ11, showed the peak of expression at PM1 and AM5 or PM5, respectively. Interestingly, FAR expression was the highest at AM1, active reproductive period. These results suggest that genes involving in pheromone biosynthesis can be sequentially regulated for their biological roles.

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 stimulated a neuropeptide, pheromone biosynthesis activating neuropeptide which is produced in subesophageal ganglion. The pheromone production 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 analyze gene expression related to pheromone biosynthesis, we performed transcriptomes of pheromone glands which were isolated at every 4 h. Eleven pheromone biosynthesis-related genes, acetyl-CoA carboxylase, fatty acid synthase, acyl-CoA dehydrogenase, enoyl-CoA hydratase, 3-hydroxyacyl-CoA dehydrogenase, Δ9 desaturase, Δ11 desaturase, fatty acid reductase, alcohol oxidase, aldehyde oxidase, and aldehyde reductase were identified. Among these genes, the expression of aldehyde reductase and aldehyde oxidase were relatively higher in the scotophase than in photophase, which may affect increase of pheromone production in the scotophase. Expression of signal genes involving in pheromone biosynthesis such as acyl-CoA desaturase, FAR, PBAN receptor, fatty acid transporter and acyl-CoA binding protein did not exhibited any significant difference.

A rapid cold hardening (RCH) and supercooling capacity usually play crucial roles in survival during the overwintering period in the tobacco budworm, Helicoverpa assulta, a freeze-susceptible species. Cryoprotectants such as polyol or sugar affect RCH and maintain fluid status of hemolymph. This study is performed to identify cryoprotectants as a RCH factor in H. assulta. Pre-exposure of H. assulta larvae to 4℃ significantly increased survival at -10℃ in all developmental stages from egg to adult. RCH was dependent on the duration of the pre-exposure period. RCH also significantly enhanced the supercooling capacity. Analysis of cryoprotectants from the hemolymph revealed that the pre-exposure treatment allowed the larvae to accumulate glycerol and trehalose among polyols examined. In addition, unknown materials from 2 peaks were also increased. TIC analysis revealed 3 predicted formulas for unknown materials, C26H24O20S or C3H4N6OS and C22H20O21. Glycerol-3-phosphate dehydrogenase (GPDH) and glycerol 3-phosphatase (G3P) that involving in glycerol biosynthesis were identified from the transcriptome of H. assulta 4th instar larvae. Based on the expression level of transcripts, the expressions of GPDH and G3P were relatively increased when compared to that of the control, suggesting that these genes contribute to overwintering and biosynthesis of glycerol.

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.

Upon freezing temperatures, most insects should avoid cellular freezing by migration to warm hibernating sites, by becoming cold-hardy or by undergoing diapause development. However, a highly threatened butterfly, Parnassius bremeri, terminates egg diapause at early winter season and grows during entire winter and spring. Thus, the cold hardiness of P. bremeri needs to be explored to understand its cold tolerance limit and physiological factors. Supercooling points (SCPs) of P. bremeri vary from -10℃ to -48℃ among season. Especially, the young larvae during Jan – Mar kept SCPs at below -20℃. Larval plasma contained high level of glycerol (39.7 mM) at March, but it decreased the level (2.4 mM) at May. Transcriptome analysis indicated high levels of gene expressions associated with glycerol synthesis. Temporal expression patterns of polyol synthesis genes supported the change of glycerol. This study suggests that glycerol is a major cryoprotectant of P. bremeri to be cold-handy against freezing temperatures during winter.

The Asian honeybee, Apis cerana, is a native honeybee species in Korea which is important in agriculture for pollination and honey production. For better understanding of the physiology of A. cerana, high-throughput Illumina transcriptome sequencing was performed to analyze the gene expression profiles of queen, worker and larva. A total of 219,799,682 clean reads corresponding to 22.2Gb of nucleotide sequences was obtained from the whole body total RNA samples. The Apis mellifera reference mRNA sequence database was used to measure the gene expression level with Bowtie2 and eXpress software and the Illumina short reads were mapped to 11,459 out of 11,736 A. mellifera reference genes. Total of 9,221 genes with FPKM value greater than 5 of each sample group were subjected to evolutionary genealogy of genes: Non-supervised Orthologous Groups (eggNOG) with BLASTX for gene ontology analysis. The differential gene expression between queen and worker, and worker and larva were analyzed to screen the overexpressed genes in each sample group. in the queen and worker sample group, total of 1,766 genes were differentially expressed with 887 and 879 genes overexpressed over two folds in queen and worker, respectively. In the worker and larva sample group, total of 1,410 genes were differentially expressed with 1,009 and 401 genes overexpressed over two folds in worker and larva, respectively.

Rice stripe virus (RSV) is one of the serious plant pathogenic viruses for rice and mediated by small brown planthopper, Laodalphax striatellus. So far, the studies have been mainly focused on the interaction between the host plant and the virus. In this study, for better comprehension of the interactions among Rice stripe virus, rice and small brown planthopper, transcriptomes of the RSV-viruliferous (RVLS) and non-viruliferous L. striatellus (NVLS) were comparatively analysed. For this, non-viruliferous L. striatellus were collected from non-infected rice field and fed RSV-infected rice for 5 days. With the RNAs prepared from the RSV-viruliferous and the non-viruliferous small brown planthoppers, we conducted Illumina RNA sequencing (Hiseq 2000) and then two transcriptome databases were generated from RVLS and NVLS, respectively. The transcriptome of RVLS and NVLS were campared to figure out how the gene expression of the insects affected by Rice Stripe Virus. RSV-dependently regulated genes analysed from this study may have important functions in the transmission and replication of RSV.

Teratocytes (TCs) are the cells derived from the embryonic serosal membrane of some parasitic hymenopteran insects. As a parasitic factor, TCs are multifunctional in host regulation by inducing nutritional, immune, and developmental alterations. However, little is understood about their genetic constituents. This study reveals a comprehensive view of the genes expressed by TCs through a transcriptome analysis based on RNAseq technology. More than 6.29 Gb sequences were used to assemble 34,686 contigs (>200 bp) and annotated into different functional categories. The TC transcriptome profile was clearly distinct from those of hemocytes and the fat body. The TC transcriptome contained components of insulin signaling and biosynthesis of juvenile hormone and 20-hydroxyecdysone. TCs also expressed various groups of digestive enzymes, supporting its nutritional role for the growing parasitoid larvae in parasitism. Furthermore, this transcriptome analysis annotated two kinds of immunosuppressive serine protease inhibitors (serpins) and Rho GTPase-activating proteins (RhoGAPs). To determine the biological functions of these factors, we devised ex vivo RNA interference (RNAi) by conducting knockdown of gene expression in in vitro cultured TCs followed by injection of the treated TCs to test insects. Ex vivo RNAi revealed that some serpins and RhoGAPs expressed in TCs inhibited host cellular immunity. This study reports a transcriptome of the unique TC animal cell, and its immunosuppressive genetic factors using ex vivo RNAi technology.

Previously, we demonstrated the presence of a second copy of LPS myristoyl transferase in enterohemorrhagic Escherichia coli O157:H7, an important zoonotic diarrheagenic food-borne pathogen; the pO157-encoded ecf (an eae-conserved fragment) and the chromosomally-encoded lpxM (also referred to as msbM) genes. Although both genes share the same function as an LPS myristoyl transferase, the pO157-encoded ecf is thermoregulated via an intrinsically curved DNA while the chromosomal lpxM is regulated by the PhoP/Q two component regulatory system. However, it is unclear why E. coli O157: H7 carries two copies of LPS myristoyl transferase that are differentially regulated. In this study, a whole genome-scale transcriptome specific to E. coli O157:H7 was carried out for identification of the genes differentially expressed in the amyristoylated E. coli O157:H7. The results identified a total of 110 EHEC genes that were up- or down-regulated in the amyristoylated E. coli O157:H7 strain, including genes associated with virulence (26.36%), metabolism (20.91%), transport (10.91%), signal transduction (4.55%), genetic information processing (3.64%), stress response (2.73%), regulatory function (2.73%), motility/adherence (3.64%), cell envelope (2.73%), cell division (1.82%) and ORFs of unknown function (17.27%). Of particular interest, the expression of LEE pathogenicity island genes was significantly influenced by LPS structural defects.

Sacbrood virus (SBV) is one of the most fatal pathogens against Asian honeybee, Apis cerana. This virus cause failure of the insect larvae to pupate and death of the adult insects. This study has analyzed the host genes affected by viral infection, by comparing the expression level of host transcripts infected with or without SBV. As a first step, we sequenced the cDNA libraries of Asian honeybee by using illumina RNA sequencing. The sequences were de novo assembled to acquire honeybee transcriptome sequences. The transcriptome was annotated by the sequence comparison to known protein sequences by BLASTX and evolutionary genealogy of genes: Non-supervised Orthologous Groups (eggNOG) database with functional categories and description. By mapping the RNA-seq data to de novo assembled transcripts, we characterized the differentially expressed transcripts between SBV-infected and non-infected Asian honeybee.

Rice stripe virus disease (RSVD), one of the most serious disease of rice is mediated through the sucking by small brown planthopper, Laodalphax striatellus. So far, the studies have been mainly focused on the interaction between the host plant and the virus. In this study, for better comprehension of the interactions among the host plant, vector insect and plant-pathogenic virus, we investigated transcriptome of the vector insect and the differences between viruliferous and naïve L.striatellus. For this, naïve L. striatellus were collected from non-infected rice field and 50 L.striatellus of them were fed RSV-infected rice for 5 days. With the RSV-viruliferous and the naïve insects, we conducted Illumina RNA sequencing (Hiseq 2000) and obtained 175,243,488 and 146,031,348 reads from viruliferous and naïve L.striatellus, respectively. These reads were assembled into contigs and two transcriptome databases were generated. The transcriptome of naïve and RSV-viruliferous L. striatellus were campared to figure out up-regulated or down-regulated genes. These RSV-dependently regulated genes may have important function in the behavior of planthoppers or the transmission of RSV.

To investigate genes differentially expressed in the venom of social and solitary wasps, a comparative transcriptome analysis was conducted. Subtractive expressed sequence tag (EST) libraries specific to the venom gland and sac (gland/sac) of a social wasp species, Vespa tropica and a solitary hunting wasp species, Rhynchium brunneum, was constructed by suppression subtractive hybridization. In BLASTx analysis, 41% and 56% of the total ESTs showed statistically best-matched hits (E ≤ 10-4) in the libraries of V. tropica and R. brunneum, respectively. Although the functional category analysis did not show remarkable differences in the distribution of functional categories between the two venom gland/sac cDNA libraries, perhaps due to the lack of functional information on many of the venom components, there were groups of genes that are specific to either V. tropica or R. brunneum. Venom allergen 5 and serine protease were found to be social wasp-specific venom transcripts. In contrast, venom peptides, metalloendopeptidases, arginine kinase and dendrotoxin were observed in solitary wasp at much higher frequencies.

Cotesia plutellae has been known as a natural enemy against the Diamondback moth, Plutella xylostella via laying eggs into a larva. When the larva hatches from the egg, teratocytes also are released and expected to work as immune suppressor via secreting immune suppressive factors. In order to analyze the gene expression in teratocytes, total RNAs were isolated and genes expressed in the teratocyte were sequenced by Illumina HiSeq2000 RNASeq analysis. The information on RNA sequences was assembled by Trinity and contigs were annotated by Blast analysis. The levels of gene expression were calculated by FPKM. Approximately, 6.3 Gbs were obtained and 34,686 contigs were found and annotated. Forty two percent of contigs were homologous to previously reported genes and classified by gene ontologies: the highly abundant components are metabolic process, biological regulation and cellular process in biological function; binding, catalytic activity and transporter activity in molecular function; cell part, membrane part and organelle in cellular function, respectively. In addition, some teratocyte transcripts of C. plutellae are related to host regulation such as immunosuppression and nutrition: Ankyrin repeat proteins, Serpin, protease, lipase, chitinase and scavenger receptor.