A residual contact vial plus water (RCVpW) bioassay method was established to monitor insectiside resistance in field populations of the melon thrips, Thrips palmi. Resistance level against six major insecticides were evalutated in five regions to test applicability of RCVpW as an on-site resistance monitoring tool. Reduced mortality in response to six test insecticides were exhibited compared to the RDA susceptable strain showing 100 % mortality, indicating different degree of resistance. An apparently reduced mortality to emamectin benzoate and chlofenapyr was observed in some field populations, suggesting uneven distribution of resistance to these insecticides in field populations. In addition, spinosad resistance was high and widely distributed in the test regions. Synergistic bioassay revealed that cytochrome P450-mediated metabolic factor is involved in spinosad resistance in the Korean population.

One of major mechanisms of insecticide resistance is the reduced susceptibility caused by mutations on the main target sites, such as sodium channel and acetylcholinesterase. Bioassay is a useful method to diagnose insecticide resistance of mosquitoes; however, it is hard to establish regional/annual resistance database based on bioassay. Recently, various markers of mutation and copy number variation have been identified through insecticide mechanism studies. Thus, molecular detection of resistance based on the resistance marker is now feasible, which can be readily implemented as a novel resistance monitoring tool to complement or replace the conventional bioassay method. In Korea, the density of vector mosquitoes native to the subtropical areas has increased due to climate change. Therefore, it is required to establish an efficient resistance monitoring system based on molecular markers to facilitate the construction of a nation-wide resistance map for more effective control of mosquitos. In addition, alternative insecticides should be introduced to the areas where mosquitoes develop high levels of resistance to maximize control efficacy against resistant populations

Honey bee, Apis mellifera L., have been widely used as a model organism for biological science because of its highly developed sociality, specialized labor division and passive population management. In order to examine the expression patterns of genes putatively involved in social development in honey bee, quantitative real-time PCR (qRT-PCR) that has been widely used to investigate the expression level of target gene can be used in honey bee study. However, the selection and validation of optimal reference genes is a crucial step prior to running qRT-PCR. In the present study, therefore, the seasonal expression stability of five candidate reference genes in the abdomen of forager and nurse was investigated using three programs (geNorm, NormFinder and BestKeeper), and selected reference genes were validated by the normalization of expression level of vg encoding vitellogenin. Although three programs revealed slightly different gene stability values, overall the combination of two genes (rpS18 and gapdh encoding ribosomal protein S18 and glyceraldehyde-3-phosphate dehydrogenase, respectively) was resulted in the most suitable use for normalization of the target gene in forager. However, a single gene, either rpL32 or rpS18 in the nurse or either rpL32, rpS18, or gapdh in the comparison between foragers and nurses, were suggested to be applied for normalization of seasonal and labor-specific gene expression by qRT-PCR.

Molecular diagnostic markers are necessary for establishing highthroughput screening systems to support insecticide-resistant population management. Here, we identified single amino acid substitution mutations related to carbamate resistance in Laodelphax striatellus Fallén type-1 acetylcholinesterase (Lsace1) using carbofuran-selected strains. The phenotypic resistance profiles of the final selection strain (SEL9) compared to the susceptible strain revealed a 14-fold higher resistance ratio based on topical application, 1.2-fold higher general esterase activity, and 4.3- fold higher acetylcholinesterase insensitivity based on the 50% inhibitory concentration (I50), suggesting that insensitivity of the target site could occur as a resistance factor. Comparison of the nucleotide sequences of Lsace1 of five strains (SUS, SEL0, SEL3, SEL6, and SEL9) revealed two amino acid substitutions (F330Y and F331H). To understand the roles of these mutations, we determined the allele frequency of both point mutations in the selected strains using quantitative sequencing methods. In addition, several quantitative genotypic traits (e.g., gene copy numbers and transcript levels of Lsace1, Lsace2, and LS.CarE1) were assessed. A correlation analysis of genotypic and phenotypic traits revealed strong correlations between resistance level and I50 with F331H allele frequency. Interestingly, the F331H mutation was negatively correlated with transcript levels of Lsace1, suggesting that selection pressure might result in a reduction of the target gene. Overall, the F331H mutation and reduced mRNA are important factors in the development of carbamate resistance. Furthermore, the point mutation can be used to monitor rapid carbofuran resistance in conjunction with molecular diagnostic methods such as quantitative sequencing.

To identify the venom components and their expression patterns of some Aculeata bees/wasps, venom gland-specific transcriptome analysis was conducted. FPKM values were normalized with the average of the transcription level of reference gene (a-tubulin). Common components in both solitary and social wasp venoms include hyaluronidase, phospholipase A2, metalloendopeptidase, etc. Although it has been expected that more diverse bioactive components with the functions of prey inactivation and physiology manipulation are present in solitary wasps, the information on venom compositions of solitary wasps obtained in this study was not sufficient to generalizae this notion. Nevertheless, some neurotoxic peptides (e.g., pompilidotoxin and dendrotoxin-like peptide) and proteins (e.g., insuline-like peptide binding protein) appear to be specific to solitary wasp venom. In contrast, several proteins, such as venom allergen 5 protein, venom acid phosphatase, and various phospholipases, appear to be relatively more abundant in social wasp venom. In the venom gland trancsriptome of bumblebees, major allergens or pain producing factors were barely identified, implying that bumblebee venoms are relatively less toxic than those of social or solitary wasps.

Lately, CRISPR-Cas9 has become one of the most essential tools to understand gene’s function. In the honeybee, however, the application of CRISPR technology has been hindered by various factors leading to very few reports of success in genome editing. Among these, collection of honeybee embryos for microinjection has been a time-consuming procedure, mainly limiting the applicability of the genome editing technique to honeybees. To improve the drawbacks of the conventional plastic plug-based system, we have developed a film-assisted honeybee embryo collection system (FECS) using transparent film as a detachable bottom layer. In this new system, eggs are laid on the detachable film surface and collected in a batch, and thus no additional alignment is required for microinjectoin. As the film can be easily replaced with in a few seconds, embryo collection can be repeated continuously after a single caging of a queen. Also, unlike conventional plug-based systems, the new system utilizes 100% of the eggs laid by the queen, thereby increases the yield three times in theory. The main unit of the system can be printed with ordinary SLA/DLP type 3D printer and the stl file for 3D printing will be distributed online.

The honey bee soluble acetylcholinesterase 1 (AmAChE1) is overexpressed under the overwintering and brood rearing-suppressed conditions. To investigate the role of AmAChE1 in regulating acetylcholine (ACh) titer, ACh concentrations both in the head (neuronal) and abdomen (non-neuronal) were analyzed. ACh titer was significantly lower in both tissues of worker bees under the overwintering and brood rearing-suppressed conditions compared to control bees. The expression levels of another two factors that regulate ACh titer, choline acetyltransferase (AmAChT) and acetylcholinesterase 2 (AmAChE2), were not altered as judged by qPCR and native PAGE, suggesting that the lower ACh titer was mainly regulated by AmAChE1. For precise verification of AmAChE1 as an ACh titer regulator, honey bees were put under brood rearing-suppressed condition to induce AmAChE1 and injected AmAChE1 dsRNA to knock down the gene. The ACh titer of AmAChE1-knocked down honey bees was 1.9 and 2.6 folds higher than that of control bees in head and abdomen, respectively. Taken together, in spite of its extremely low catalytic activity, the overexpression of AmAChE1 is likely to be related with the low level of ACh homeostasis, perhaps via ACh sequestration, under brood rearingsuppressed condition, and likely induce metabolic changes through ACh receptors-related pathways.

Residual contact vial (RCV) method was used to monitor insecticide resistance in field populations of the melon thrips, Thrips palmi. Median lethal doses (LD50) at 8 h post-treatment of six insecticides (chlorfenapyr, cyantraniliprole, cypermethrin, dinotefuran, emamectin benzoate and spinosad), which are commonly used for T. palmi control, were determined at 8 h post-treatment using a susceptible RDA strain. The diagnostic doses for on-site resistance monitoring of the six insecticides, which were determined as two-fold higher doses of LD90 for the RDA strain, were in the range of 0.299 to 164,25 μg-1cm2. To test the applicability of RCV for T. palmi, insecticide resistance levels in three field populations (Gyeonggi; GG_AS, Chungbuk; CB_CJ, Jeonbuk; JB_KJ) were evaluated. Field populations showed reduced mortality (0-50% mortality) to spinosad, cypermethrin, cyantraniliprole and emamectin benzoate, that they have different degree of resistances to these insecticides. In particular, all test field populations exhibited 0% mortality to spinosad, suggesting wide spread of spinosad resistance in the field. Moreover, no detectable mortality to emamectin benzoate was observed in JB_KJ strain, suggesting uneven distribution of emamectin benzoate resistant population of T. palmi. To provide more precise information on resistance profiles and distribution in T. palmi populations, it would be necessary to conduct a large scale resistance mapping for broad geographical regions.

Laodelphax striatellus is an important pest of rice due to not only sucking rice seedlings, but also transmitting serious plant viruses. Among various kinds of insecticide groups (carbamates, organophosphorus, neonicotinoids, etc.), carbofuran, a systemic carbamate insecticide, has been most extensively used to control rice pests including L. striatellus, resulting in widespread carbamate resistance in Korea and other Northeast Asia countries. To identify the genes associated with carbofuran resistance, we obtained a 14-fold higher resistant strain (SEL9) from the mixed-field population (SEL0) by consecutive selection. A transcriptome-based analysis was conducted and differentially expressed genes (DEG) were compared between the SEL9 and SEL0 strains. A total of 96,185,150 reads were analyzed, of which 62,860,430 reads were mapped. From these reads, 15,356 transcripts were annotated. A total of 327 up-regulated and 275 down-regulated genes were identified in the resistant SEL9 strain compared to SEL0 strain by DEG analysis. Gene ontology (GO) analysis was performed using DEGs showing statistical significance (P < 0.05). The GO analysis identified 1,320 genes in biological process group, 1,100 genes in cellular component group and 428 genes in molecular_function group.

Frankliniella occidentalis is a notorious polyphagous crop pest causing tremendous economic loss. It damages flowers and leaves of host plants and also carries severe plant viruses. During last few decades, it has spread to all continents via transport of plant materials. Following extensive use of insecticides to control F. occidentalis, it has developed high level of resistance due to its short life cycle and high reproductive potential. In this study, RNA interference (RNAi)-based bioassay system was developed to find an alternative control measure for insecticide-resistant population of F. occidentalis. A variety of genes involved in various physiological mechanisms were selected for the test of dsRNA potency (tubulin, v-ATPase, amylase, aquaporin etc.). Each bioassay unit made by 3D printing has a leaf disc placed on 150 ㎕ of 50 ng/ul dsRNA solution and 20 thrips. The mortality was checked, and the dsRNA and leaf disc were replaced every 24 h for 72 h. Of the 20 genes tesetd, tubulin, v-ATPase, and aquaporin showed 31, 38, 38 and 45% of corrected mortality at 72 h post-treatment, respectively. This result suggests the potential of these genes as candidate lethal genes for RNAi-based F. occidentalis control system.

Small hive beetle (Aethina tumida) (SHB) is an invasive species to most northern hemisphere countries, including Korea. In an attempt to obtain basic information for efficient management of SHB, genes encoding conventional insecticide targets [voltage-sensitive sodium channel α-subunit (VSSC) and acetylcholinesterase (AChE)] were annotated and characterized following the analysis of whole transcriptomes of adults and larvae. A single VSSC gene was identified but no apparent mutations associated with pyrethroid resistance were detected. Genes encoding two AChEs (AtAChE1 and AtAChE2) were identified from the SHB transcriptome. AtAChE1 was determined to be the main catalytic enzyme, thereby being a toxicologically more relevant target. No apparent mutations associated with resistance to organophosphorus and carbamate insecticides was identified in the AtAChE1 gene, whereas the S238G mutation, originally identified from the Colorado potato beetle, was detected in the AtAChE2 gene.

The common bed bug, Cimex lectularius, possesses a cholinesterase expressed exclusively in the salivary gland (ClSChE). In this paper, we investigated the molecular structure, tissue distribution patterns, and biochemical properties of ClSChE and showed that ClSChE exists as a soluble monomeric form or a soluble dimeric form connected by a disulfide bridge. Immunohistochemical analysis confirmed that ClSChE was expressed in the epithelial cells of both the salivary gland and the duct. In addition, the secretion of monomeric ClSChE through the proboscis during feeding was detected by western blotting using a ClSChE-specific antibody. To predict the role of ClSChE injected into the tissue of an animal host, we analyzed the extent of sequestration and hydrolysis of acetylcholine (ACh)/choline (Ch) by ClSChE by ultra-performance liquid chromatography-tandem mass spectrometry. Kinetic analysis revealed that ClSChE possesses extremely low Km (high affinity to ACh) and Vmax values. These findings suggest that ClSChE functions as a sequestering enzyme specific to ACh (not to Ch) by having a very strong affinity to ACh but an extremely long turnover time.

The honey bee soluble acetylcholinesterase 1 (AmAChE1) is overexpressed under the overwintering and brood rearing-suppressed conditions. To investigate the role of AmAChE1 in regulating acetylcholine (ACh) titer, ACh concentrations in both the head (central nervous system) and abdomen (peripheral nervous system) were analyzed. ACh titer was significantly lower in both tissues of worker bees under the overwintering and brood rearing-suppressed conditions compared to control bees. Interestingly, the expression levels of choline acetyltransferase (AmChAT) and molecular marker genes of immune systems were significantly reduced in honey bee head under the same conditions. Taken together, ACh titer appears to be reduced via a cooperative interaction of the AmAChE1 overexpression and AmChAT underexpression and to be linked to reduced inmmune responses under the overwintering and brood rearing-suppressed conditions. The roles of AmAChE1 (with little catalytic activity) and AmChAT in the ACh homeostasis and signaling was discussed in the contexts of immune response and longevity regulation in honey bees.

It is necessary to understand of temporal and spatial dynamics by establishing a periodical monitoring system for theproper management in small brown planthopper (SBPH). A dataset is including the number of SBPHs by location, collectionmethod [aerial collection net (AeCN) or light trap (LT)] and period (May~Aug.) for five years (2011~2015), and missingvalues were imputed using multiple imputation methods. Of the 15,848 individuals collected, approximately 47% and 52.9%were collected using the AeCN and LT methods, respectively. A high incidence of migratory SBPHs was observed duringJulian days 144-166 using the AeCN method. Generally, the migratory SBPHs from China composed 39.4% of the totalpopulations of SBPHs. These results would provide valuable information to predict the incidence period of migratory SBPHsand establish a proactive management system against SBPH.

We compared the genetic structures of overwintered indigenous Korean and Chinese populations. The eight Koreanpopulations consisted of 33 haplotypes, and 16 haplotypes were newly identified. The genetic diversity of the Koreanpopulation revealed high haplotype diversity and low nucleotide diversity of 0.86 and 0.0024 on average, respectively,due to the high dispersal ability, which is similar to that of the Chinese population (Sun et al., 2015). Comparison with30 Chinese populations using a population tree showed that the Korean populations grouped with 12 Chinese populationsand that 67% were located near Jiangsu province. Moreover, the three frequent migration regions by migratory SBPH,including the Buan, Shinan and Taean counties, were grouped together with high supporting values. These results mightsupport the presence of gene flow between the Korean and Chinese populations by migratory SBPHs.

Among two different acetylcholinesterase (AmAChE1 and AmAChE2) of the western honey bee, the soluble AmAChE1might be related with a stress response as judged from its over-expression in honey bee workers when brood rearingwas suppressed. In this study, to ensure the nature of AmAChE1 responding to stress factors, the expression patternsof AmAChE1 were investigated following various treatments, including varroa mite infestation, bacterial challenge, broodrearing suppression, thermal stresses, chemical treatments, ultraviolet B irradiation, starvation, water restriction and crowdingstress. In addition, transcription profiles of four heat shock protein genes known as general stress markers and vitellogeningene, which is induced in several stress conditions, were tested as positive references. In every tested condition, onlybrood rearing suppression and heat shock were related with the expression of AmAChE1.

Biological activities of bombolitins from Bombus ardens, B. consobrinus, B. terrestris and B. ussurensis (bombolitinsA, C, T and U, respectively) were examined using hemolytic, anti-microbial, anti-fungal and anti-tumor activity assays.Among the four bombolitins tested, bombolitin T showed the highest hemolytic and anti-tumor activities. All bombolitinsexhibited strong anti-microbial and anti-fungal activities, and bombolitin A specifically possessed the highest anti-microbialactivity against the Gram-negative bacteria Escherichia coli. Circular dichroism spectrometry analysis revealed that allfour bombolitins had over 61.7% and 45.5% of α-helicity in 30 mM sodium dodecyl sulfate and 50% trifluoroethanolbuffers, respectively, which form lipid-membrane-mimicking environments. Bombolitin T showed the lowest IC50 valuesof 8.5 μM and 8.8 μM against SK-OV-3 and NIH-OVCAR-3 cell lines, respectively, after 72 h of treatment, but itsrelative hemolytic activity at a concentration of 200 μM was 2.3-fold higher than that of 0.1% Triton X-100. Thisstudy provides new information on the biological and molecular properties of venom peptides of bumblebees. However,further studies on reducing cytotoxicity of bombolitins are needed for designing selective anti-tumor peptides.

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.

Perturbation of normal behaviors (e.g., nursing and foraging) in honey bee colonies by any external factors would immediately reduce the colony’s capacity for brood rearing, which can eventually lead to collapse of entire colony. To investigate the effects of brood rearing suppression in the biology of honey bee workers (nurse and forager), the gene-set enrichment analysis (GSEA) was performed for the transcriptomes of worker bees with or without their brood rearing being suppressed, from which functional profiles of pathways under influences by each condition were identified. Blocking of normal labor (i.e., nursing or foraging) induced the over-representation of pathways related with reshaping of worker bee physiology, suggesting that transition of labor is physiologically reversible. In addition, brood rearing suppression appeared to result in the reduction of neuronal excitability and aggressiveness in both forager and nurse, which would be necessary to manage the in-hive stress under unfavorable conditions

Honey bee swarming is a natural phenomenon that occurs by changes of colony (i.e. population size and queen condition) and environment conditions. As cuticular hydrocarbons (CHCs) are known to be involved in the communication between honey bee nest-mates, we investigated and compared the CHC profiles of worker bees from individual colonies of 9-days before swarming (PPSC), a day before swarming (PSC), swarming (SG) and remaining (non-swarming) (RG). A total of 53 CHCs were identified by GC-MS, among which 11 compounds showed significantly differential expression patterns between swarming states. Before swarming (between PPSC and PSC), detection levels of 4 CHCs were significantly different, suggesting that production of some CHCs changed prior to swarming for swarming preparation. Six CHCs were deferentially produced between PSC and RG. The differential profiles of CHCs with respect to different swarming states are currently under investigation.