Acetylcholinesterase 1 (AmAChE1) has low catalytic activity and is abundantly expressed in both neuronal and non-neuronal tissues. In previous experiments, we observed that AmAChE1 is rarely expressed in summer while highly expressed in winter. Through additional experiments, the expression of AmAChE1 was suggested to be associated with brood rearing status. Under the assumption that abnormal suppression of brood rearing activity may result in stressful condition in honey bee social community, it was further suggested that AmAChE1 is likely involved in stress management particularly during winter. We hypothesized that the increased docility usually observed in overwintering bees is likely an outcome of stress management in colony, which is mediated by AmAChE1 expression. To verify this, worker bees expressing abundant AmAChE1 were collected in early winter and injected with Amace1 dsRNA to knockdown Amace1. Then, the behavioral activity of the bees was investigated using the EthoVison video tracking system. Honey bees injected with Amace1 dsRNA showed significantly increased motility, which was strongly correlated with the suppressed expression level of AmAChE1 in the abdomen. No apparent reduced expression of AmAChE1 in the head was observed perhaps due to the limited efficacy of RNA interference in the blood-brain-barrier. Our finding suggests that behavioral activity can be regulated, at least, by AmAChE1 expression level in non-neuronal tissue (i.e., fatbody) perhaps via metabolic alteration.

Human body and head lice are obligatory human ectoparasites. Although both body and head lice belong to a single species, Pediculus humanus, only body lice are known to be a vector of several bacterial diseases. The higher vector competence of body lice is assumed to be due to their weaker immune response than that of head lice. To test this hypothesis, immune reactions were compared between body and head lice following infections by two model bacteria, Staphylococcus aureus and Escherichia coli, and a human pathogen, Bartonella quintana. Following dermal or oral challenge, the number of these bacteria increased both in hemocoel and alimentary tract of body lice but not in head lice and the viability of the B. quintana was significantly higher in body louse feces, the major route of infection to human. In addition, body lice showed the lower basal/induced transcription level of major immune genes, cytotoxic reactive oxygen species and phagocytosis activity compared with head lice. These findings suggest that a reduced immune response may be responsible, in part, for the increased proliferation and excretion of viable bacteria which are associated with the high level of human infectivity seen in body versus head lice.

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 acetylcholinesterase 1 (AmAChE1) of the honey bee is known to be abundantly expressed both in the central and peripheral nervous systems. AmAChE1 exists mostly in the soluble form with little catalytic activity and has non-neuronal functions. Our preliminary observation showed that AmAChE1 expression fluctuated between the forages and nurses. A more systematic expression profiling of AmAChE1 over a year cycle on a monthly basis revealed that AmAChE1 was predominantly expressed during the winter months with being moderately expressed during the rainy summer time. However, no significant difference in AmAChE1 expression was noticed between the nurse and forager workers. Interestingly, AmAChE1 expression was inhibited when bees were allowed for brooding by placing overwintering bee hives in strawberry green houses with the supplement of pollen diets whereas it was resumed when the bee hives were removed from the green houses, thereby suppressed brooding. To confirm whether brooding status is a main determining factor for the suppression of AmAChE1 expression, active bee hives were placed in a screen tent, thereby hindering foraging, until brooding was completely suppressed, and then allowed to restore brooding by removing the screen. The AmAChE1 expression in the head was up-regulated when brooding was suppressed whereas its expression was down-regulated when brooding was resumed. These finding demonstrates that AmAChE1 expression in the central nervous system (i.e., head) is related with brooding status of honey bee. To understand the connection between the AmAChE1 expression and other pathways related with brooding, currently in progress are the analyses of head transcriptomes of honey bee workers with or without their brooding suppressed.

Leptotrombidium pallidum is the major vector mites for Orientia tsutsugamushi, the causative agent of scrub typhus. To understand the molecular mechanism of L. pallidum, we sequenced the whole genome using Illumina sequencing technology. Totally four genomic libraries with different insert sizes ranging from 280 bp to 8 kb were used to generate 45.1 Gb of genome in the combination of paired-end and mate-pairs sequencing reads. Quality filtering and correction of paired-end reads for very small and/or bad-quality sequences yielded 26.9 Gb of high-quality sequences, which are used to estimate the genome size as 175 Mbusing kmer methods and assembled into a 193.7 Mb genomic sequence scaffolds with N50 length of 92,945 bp. Furthermore, 94% of CEGMA completeness score were obtained from genome scaffold assembly. To facilitate gene annotation, we used a combination of de novo and homology based tools to predict gene models in the chigger mite genome. A combination of evidence-based and de novo approaches predicted 15,842 high-confidence protein-coding genes with an average transcript length of 1,511 bp and 2.4 exons per gene which corresponds to about 12.4% total gene length. Bacterial endosymbiosis are very common in mite species and can range from mutualistic to pathogenic associations. Henceforth, the endosymbionts in L. pallidum were predicted using the NCBI microbial draft genomes and mitochondrial genome. Besides, this L. pallidum draft genome can be used as a significant reference for comparative genomic studies across mite species.

Leptotrombidium pallidum is the major vector mite for Orientia tsutsugamushi, the causative agent of scrub typhus, in Asian countries, including Korea. The genome size of L. pallidum was previously estimated to be 191 ± 7 Mb (Kim et al., 2014). Genomic DNA (gDNA) was extracted from a single female from a 9-generation inbred L. pallidum colony and used for whole genome amplification (WGA). The resulting amplified gDNA was used for the construction of paired-end and mate-pair libraries and sequenced using Illumina platforms (HiSeq2000 and MiSeq). An unamplified gDNA sample extracted from 20 female mites was also used for sequencing in parallel. More than 45Gb sequence reads from both paired-end and mate-pair libraries of the WGA gDNA were trimmed and then de novo assembled using the CLC Asembly Cell v.4.0 for contig assembly and SSPACE for scaffolding. The assembly generated approximately 6,545 scaffolds with N50 value of 92,945 and total size of ~193Mb, which was in a good agreement with our previous estimation. Repeat analysis showed that about 30% of genome (~58Mb) was masked as repeats, most of which were unclassified novel elements. For gene predictions, generated were the PASA models based on genomic alignments of RNA-seq reads from 4 different chigger mite samples (i.e. male, female, larva, and protonymph) and the GeneWise models based on genomic alignments of protein sequences from 4 closely related species with chigger mite. Independently, ab initio gene predictions were performed with AUGUSTUS and FgeneSH with custom trained matrices optimized for L. pallidum and GENEID with pre-trained matrix for Acyrthopsiphon pisum. By combining all together, 15,842 genes were predicted finally. Manual curation is in progress for various groups of genes, including chemosensory receptor genes, immune-related genes, acaricide target genes, etc.

Recently, the expression of acetylcholinesterase1 (AChE1) in honeybee worker has been found to be seasonally fluctuated. Seasonal investigation on the AChE1 expression profiles revealed that it is abundantly expressed in January but its expression was completely abolished in February in both head and abdomen. In an attempt to predict the physiological function of seasonally expressed AChE1, proteomic analysis of honeybee worker was conducted using the samples collected in January and February. Total protein samples separately extracted from the head and abdomen of honeybee forager were compared by 2-D electrophoresis (2-DE). More than 2-fold differences in expression patterns between the two different samples were observed in 50 and 85 protein spots in the head and abdomen, respectively. Among them, 20 protein spots showing >17-fold differences in expression between the two different samples of the head were identified by mass spectrometry. Most of the proteins were identified to be the major royal jelly protein (MRJP) families (e.g., MRJP, MRJP2 and MRJP3), which are known to be expressed in nurse bees during brooding season, and their expression was significantly higher in January than in February. This result was unexpected because brooding usually began in the study site apiary during February and the worker bees used for analysis were assumed to be foragers (old workers). Thus, current findings suggest, though speculative, that the workers collected in January may function as nurses despite their old ages in January or that MRJPs may have other not-yet-characterized functions, which is apart from the conventionally known roles. Finally, possible association of MRJPs with AChE1 was discussed.

Leptotrombidium pallidum and Leptotrombidium scutellare are the major vector mites for Orientia tsutsugamushi, the causative agent of scrub typhus. Before these organisms can be subjected to whole-genome sequencing, the genome sizes of L. pallidum and L. scutellare were estimated by a method based on quantitative real-time PCR. In addition, k-mer analysis of the genome sequences obtained from Illumina sequencing was conducted to verify the mutual compatibility and reliability of results. The genome sizes estimated by qPCR were 191.3±7 Mb for L. pallidum and 262.1±13 Mb for L. scutellare. The estimated genome sizes based on k-mer analysis were 175.5 Mb for L. pallidum and 286.6 Mb for L. scutellare. The estimates from two independent methods were mutually complementary and in a similar range to those of other Acariform mites. The relatively small genome size would facilitate genome analysis, which could contribute to understanding Arachnida genome evolution and mite vector competence and provide key information for scrub typhus prevention.

The chigger mite, Leptotrombidium pallidum, is widely distributed throughout South Korea and is a major vector for Orientia tsutsugamushi, the causative agent of scrub typhus. In this study, the genome size of the chigger mite was estimated to determine the necessary coverage level prior to whole genome sequencing. Cloning of EF1α and RpS3 as putative single copy reference genes were conducted and their partial sequences were determined. Using the serially diluted reference genes with known amount as standard templates, the weight of a single copy of the genome was predicted by a method based on quantitative real time PCR. The average genome length estimated from the weight using two methods was 191 ± 7 Mb. When the genome size of other arthropods (Drosophila melanogster, Apis mellifera and Tetranychus urticae), with their genome analysis completed, were estimated using the same method and compared with actual values, the estimation accuracy was 79.8-98.9%, suggesting our current estimation of L. pallidum genome size is reliable. The estimated L. pallidum genome size is in a similar range to other Acariform mites, such as the dust mite and scabie mite, but appoximately 10-fold smaller compared to the deer tick, which belongs to Parasitiform. Our finding provides key information for further genome sequencing and understanding of mite genome evolution.

The body and head lice (Pediculus humanus humanus and Pediculus humanus capitis, respectively) are hematophagous ectoparasites of humans and only the body louse between two is known to transmit three bacterial diseases through its feces. The proliferation profiles of Bartonella quintana, the causative agent of trench fever, inside the louse body and its excretion patterns were investigated in the two louse subspecies following oral challenge with B. quintana-infected blood meal. The initial density of B. quintana was sustained inside head lice without any noticeable proliferation for the entire period after infection. In contrast, B. quintana proliferated rapidly inside body lice and the maximum density reached at 10 days post-infection. The numbers of bacteria detected in feces from infected lice were almost the same and steadily decreased over time in both body and head lice. Nevertheless, the viability of the bacteria, as determined by fluorescence, was significantly higher in body louse feces, especially at 1 day post-infection and this tendency lasted for 11 days. These findings suggest that excretion of feces containing more viable B. quintana that is proliferated inside body lice following ingestion of infected blood meal is responsible for the higher vector competence of body lice.

The screening of effective lethal genes was conducted via the systemic delivery of dsRNA for the RNAi-based management of Tetranychus urticae. Six candidate genes (COPI coatmer, T_COPI; ESCRT III_Snf7, T_SNF7; Ribosomal protein S4, T_RPS4; v-ATPase A subunit 2, T_V-ATPase; Aminopeptidase N, T_APN3; Acetylcholinesterase, T_AChE) and two reference genes (EGFP and T_AChEintron) were tested for the experiment. The permeated dsRNA to the leaf disc (ca. 30 mm diameter) was detected at 6 h after treatment, indicating that dsRNA could move through veins on the leaf. In the reference gene selection, the T_AChEintron was chosen for its low mortality compared with EGFP gene. In the evaluation of mortality, the T_COPI, T_V-ATPase and T_RPS4 exerted higher toxicities at 24 hour after treatment among six genes tested. Interestingly, T_APN3 showed toxicity after 72 hour. In summary, the dsRNA delivery via leaf disc was effective in screening lethal genes and some genes, such as COPI, V-ATPase and RPS4, can be applicable for establishing a RNAi-based control system against T. urticae.

The resistance levels against carbamates (CB) and organophosphates (OP) were determined through bioassay and quantitative sequencing (QS) methods in 16 field populations of Nilaparvata lugens. The resistance levels to CB and OP by bioassay were 1.3~47.5-fold and 1.4~14.4-fold higher than a susceptible strain, respectively. The QS protocol was established to determine the allele frequencies of eight point mutations on acetylcholinesterase putatively associated CB and OP resistance. The allele frequencies of four mutations in local populations (G119A, F/S330Y, F331H and I332L) ranged from ca. 0.0~51.7%, 1.0~44.3%, 8.5~57.3% and 7.12~56.6%, respectively. The average prediction limits were –9.6±5.1~7.7±2.9%. The F330Y, F331H and I332L were tightly linked each other, suggesting these mutations may occur simultaneously. In the correlation analysis, G119A was not well correlated with both insecticides (r2= less 0.25), whereas F/S330Y, F331H and I332L showed better correlation with the resistance levels of carbamate (r2=0.590) than organophosphate (r2=0.235). This finding indicates that F/S330Y, F331H and I332L mutation frequencies are suitable for detecting carbamate resistance in N. lugens. QS will be applicable for the rapid monitoring of resistance levels to CB insecticides in N. lugens.

The human body and head louse are ectoparasites of humans for thousands of years. Although both body and head lice belong to a single species, Pediculus humanus, only body lice are known to transmit several bacterial diseases. This different vector competence is assumed to be due to their different immune responses. Here, the immune reactions in the alimentary canal were investigated in both two louse subspecies following oral challenge of Escherichia coli as a model gram-negative bacteria. In propagation assay, head lice suppressed the proliferation of E. coli in their epithelial cells effectively at the early stage of infection, resulting in gradual reduction of E. coli number in gut tissues. In contrast, the number of E. coli steadily increased in gut tissues of body lice. No apparent alteration of transcription was observed following E. coli challenge in three important genes for the humoral immune responses, PGRP as a recognition gene and defensin1 and 2 as effector genes. Nevertheless, the basal transcription levels of these genes were higher in the gut tissues of head lice than body lice. Considering that there is no cellular immune reactions in gut tissues, these findings suggest that the higher constitutive transcription levels of major immune genes in head lice can contribute to their initial defense and immune capacity against intestinal bacterial infection.

Body and head lice (Pediculus humanus humanus and Pediculus humanus capitis, respectively) are typical ectoparasites of humans. They differ not only in the ecological habitat but also in the vector competence in spite of their conspecific nature. Only body lice transmit several bacterial pathogens to humans, including Bartonella quintana, Rickettsia prowazekii and Borrelia recurrentis. In this study, the proliferation rates of two model bacteria, a gram positive Staphylococcus aureus and a gram negative Escherichia coli, were determined following bacterial challenge by cuticular injection. Both bacteria proliferated rapidly in body lice at the early stage of bacterial challenge but not in head lice, suggesting that head lice have more sensitive immune responses to these bacteria. In vivo phagocytosis assay revealed that head lice have much higher phagocytic activity against E. coli than body lice whereas only slight differences in phagocytic activity against S. aureus were observed between the two lice species. Taken together, these findings suggest that the reduced phagocytosis activity of body lice contributes, at least in part, to their higher vector competence.

The differences in the immune response between body lice, Pediculus humanus humanus, and head lice, Pediculus humanus capitis, were regarded as primary factors determining their differential vector competence. To find any differences in genetic components in immune system between body and head lice, whole genome sequences of head lice were determined by both SBS [sequencing by synthesis, Illumina Genome Analyzer (Illumina-GA)] and pyrosequencing (Roche GS FLX), and compared with the reference genome sequences of body lice. The short DNA reads from Illumina-GA (an average mapping depth of 50-fold) were aligned first to the body louse reference genome, to which Roche GS FLX DNA reads (an average depth of 2.5-fold) were subsequently assembled to make up gaps between mapped consensus. Total consensus showed a size of 114 Mb and a coverage of 96% of the published body louse genome sequences. From this head louse genome sequences, a total of 12,651 genes were predicted and used for comparing with the 10,775 genes previously reported from the body louse genome. The homolog analysis identified 873 head louse-specific genes and 422 body lice-specific genes. Comparison of immune response genes between both louse species showed head lice have more number of immune-related genes than body lice. Head lice were determined to possess all of the 107 immune-related genes reported in the previous study (Kim et al., 2011), suggesting that there is no difference in genetic make-up in terms of the 107 immune-related genes between body and head lice.