A major endogenous glycoside hydrolase mediating quercetin uptake in Bombyx mori.

Quercetin is a common plant flavonoid which is involved in herbivore-plant interactions. Mulberry silkworms (domestic silkworm, Bombyx mori, and wild silkworm, Bombyx mandarina) take up quercetin from mulberry leaves and accumulate the metabolites in the cocoon, thereby improving its protective properties. Here we identified a glycoside hydrolase, named glycoside hydrolase family 1 group G 5 (GH1G5), which is expressed in the midgut and is involved in quercetin metabolism in the domestic silkworm. Our results suggest that this enzyme mediates quercetin uptake by deglycosylating the three primary quercetin glycosides present in mulberry leaf: rutin, quercetin-3-O-malonylglucoside, and quercetin-3-O-glucoside. Despite being located in an unstable genomic region that has undergone frequent structural changes in the evolution of Lepidoptera, GH1G5 has retained its hydrolytic activity, suggesting quercetin uptake has adaptive significance for mulberry silkworms. GH1G5 is also important in breeding: defective mutations which result in discoloration of the cocoon and increased silk yield are homozygously conserved in 27 of the 32 Japanese white-cocoon domestic silkworm strains and 12 of the 30 Chinese ones we investigated.

Combination of target site mutation and associated CYPs confers high-level resistance to pyridaben in Tetranychus urticae.

Pyridaben is a mitochondrial electron transport complex I inhibitor. The H110R mutation in the PSST subunit has been reported as a major factor in pyridaben resistance in the two-spotted spider mite, Tetranychus urticae. However, backcross experiments revealed that the mutant PSST alone conferred only moderate resistance. In contrast, inhibition of cytochrome P450 (CYP) markedly reduces resistance levels in a number of highly resistant strains. It was reported previously that maternal factors contributed to the inheritance of pyridaben resistance in the egg stage, but the underlying mechanisms have yet to be elucidated. Here, we studied the combined effects of the PSST H110R mutation and candidate CYPs, as metabolic resistance factors, on pyridaben resistance in T. urticae. We found that the maternal effects of inheritance of resistance in the egg stage were associated with CYP activity. Analysis of differential gene expression by RNA-seq identified CYP392A3 as a candidate causal factor for the high resistance level. Congenic strains, where the alleles of both PSST and CYP392A3 were derived from a resistant strain (RR_i; i = 1 or 2) and a susceptible strain (SS_i) in a common susceptible genetic background, were constructed by marker-assisted backcrossing. RR_i showed upregulation of CYP392A3 and high resistance levels (LC50 > 10,000 mg L-1), while SS_i had LC50 < 10 mg L-1. To disentangle the individual effects of PSST and CYP392A3 alleles, we also attempted to uncouple these genes in RR_i. We conclude that given the variation in LC50 values and expression levels of CYP392A3 in the congenic and uncoupled strains, it is likely that the high pyridaben resistance levels are due to a synergistic or cumulative effect of the combination of mutant PSST and associated CYPs, including CYP392A3, but other yet to be discovered factors cannot be excluded.

Molecular mechanism underlying juvenile hormone-mediated repression of precocious larval-adult metamorphosis.

Juvenile hormone (JH) represses precocious metamorphosis of larval to pupal and adult transitions in holometabolous insects. The early JH-inducible gene Krüppel homolog 1 (Kr-h1) plays a key role in the repression of metamorphosis as a mediator of JH action. Previous studies demonstrated that Kr-h1 inhibits precocious larval-pupal transition in immature larva via direct transcriptional repression of the pupal specifier Broad-Complex (BR-C). JH was recently reported to repress the adult specifier gene Ecdysone-induced protein 93F (E93); however, its mechanism of action remains unclear. Here, we found that JH suppressed ecdysone-inducible E93 expression in the epidermis of the silkworm Bombyx mori and in a B. mori cell line. Reporter assays in the cell line revealed that the JH-dependent suppression was mediated by Kr-h1. Genome-wide ChIP-seq analysis identified a consensus Kr-h1 binding site (KBS, 14 bp) located in the E93 promoter region, and EMSA confirmed that Kr-h1 directly binds to the KBS. Moreover, we identified a C-terminal conserved domain in Kr-h1 essential for the transcriptional repression of E93 Based on these results, we propose a mechanism in which JH-inducible Kr-h1 directly binds to the KBS site upstream of the E93 locus to repress its transcription in a cell-autonomous manner, thereby preventing larva from bypassing the pupal stage and progressing to precocious adult development. These findings help to elucidate the molecular mechanisms regulating the metamorphic genetic network, including the functional significance of Kr-h1, BR-C, and E93 in holometabolous insect metamorphosis.

Functional crosstalk across IMD and Toll pathways: insight into the evolution of incomplete immune cascades.

In insects, antimicrobial humoral immunity is governed by two distinct gene cascades, IMD pathway mainly targeting Gram-negative bacteria and Toll pathway preferentially targeting Gram-positive bacteria, which are widely conserved among diverse metazoans. However, recent genomic studies uncovered that IMD pathway is exceptionally absent in some hemipteran lineages like aphids and assassin bugs. How the apparently incomplete immune pathways have evolved with functionality is of interest. Here we report the discovery that, in the hemipteran stinkbug Plautia stali, both IMD and Toll pathways are present but their functional differentiation is blurred. Injection of Gram-negative bacteria and Gram-positive bacteria upregulated effector genes of both pathways. Notably, RNAi experiments unveiled significant functional permeation and crosstalk between IMD and Toll pathways: RNAi of IMD pathway genes suppressed upregulation of effector molecules of both pathways, where the suppression was more remarkable for IMD effectors; and RNAi of Toll pathway genes reduced upregulation of effector molecules of both pathways, where the suppression was more conspicuous for Toll effectors. These results suggest the possibility that, in hemipterans and other arthropods, IMD and Toll pathways are intertwined to target wider and overlapping arrays of microbes, which might have predisposed and facilitated the evolution of incomplete immune pathways.

T929I and K1774N mutation pair and M918L single mutation identified in the voltage-gated sodium channel gene of pyrethroid-resistant Thrips tabaci (Thysanoptera: Thripidae) in Japan.

Pyrethroid-resistance in onion thrips, Thrips tabaci, has been reported in many countries including Japan. Identifying factors of the resistance is important to correctly monitoring the resistance in field populations. To identify pyrethroid-resistance related genes in T. tabaci in Japan, we performed RNA-Seq analysis of seven T. tabaci strains including two pyrethroid-resistant and five pyrethroid-susceptible strains. We identified a pair of single point mutations, T929I and K1774N, introducing two amino acid mutations, in the voltage-gated sodium channel gene, a pyrethroid target gene, in the two resistant strains. The K1774N is a newly identified mutation located in the fourth repeat domain of the sodium channel. Genotyping analysis of field-collected populations showed that most of the T. tabaci individuals in resistant populations carried the mutation pair, indicating that the mutation pair is closely associated with pyrethroid-resistance in Japan. Another resistance-related mutation, M918L, was also identified in part of the resistant populations. Most of the individuals with the mutation pair were arrhenotokous while all individuals with the M918L single mutation were thelytokous. The result of differentially expressed gene analysis revealed a small number of up-regulated detoxification genes in each resistant strain which might be involved in resistance to pyrethroid. However, no up-regulated detoxification genes common to the two resistant strains were detected. Our results indicate that the mutation pair in the sodium channel gene is the most important target for monitoring pyrethroid-resistance in T. tabaci, and that pyrethroid-resistant arrhenotokous individuals with the mutation pair are likely to be widely distributed in Japan.

The R81T mutation in the nicotinic acetylcholine receptor of Aphis gossypii is associated with neonicotinoid insecticide resistance with differential effects for cyano- and nitro-substituted neonicotinoids.

The cotton aphid, Aphis gossypii Glover, is one of the most agriculturally important insect pests. Neonicotinoid insecticides and sulfoxaflor have generally shown excellent control of A. gossypii, but these aphids have recently developed resistance against neonicotinoid insecticides. We previously characterized a field-collected A. gossypii Kushima clone that showed higher resistance to nitro-substituted neonicotinoids, such as imidacloprid, than to cyano-substituted neonicotinoids, such as acetamiprid. This Kushima clone harbors the R81T mutation in the nicotinic acetylcholine receptor (nAChR) ß1 subunit; this mutation is the source of neonicotinoid insecticide resistance. In the present study, electrophysiological analyses and molecular modeling were employed to investigate the differential effects of the R81T mutation on cyano- and nitro-substituted neonicotinoids and sulfoxaflor. We isolated full-length coding sequences of A. gossypii nAChR α1, α2, and ß1 subunits. When co-expressed in Xenopus laevis oocytes with chicken ß2 nAChR, A. gossypii α1 evoked inward currents in a concentration-dependent manner in response to acetylcholine (ACh) and showed sensitivity to neonicotinoid and sulfoxaflor. Additionally, the chicken ß2 T77R+E79V (equivalent double mutant of R81T) mutation resulted in a lower effect to cyano-substituted neonicotinoids and sulfoxaflor than to nitro-substituted neonicotinoids. Electrophysiological data and nAChR homology modeling analysis suggested that the Kushima clone exhibited different levels of resistance to cyano- and nitro-substituted neonicotinoid insecticides.

Field survey of reproductive modes and sodium channel mutations associated with pyrethroid resistance in Thrips tabaci.

Using PCR-Restriction Fragment Length Polymorphism (RFLP) with mitochondrial cytochrome c oxidase subunit I sequences, we examined the reproductive modes of female adults of Thrips tabaci collected at 54 sites across Japan. Results showed the presence of heteroplasmic insects harboring mitochondria associated with arrhenotoky and thelytoky. Using the insects, we also applied PCR-RFLP to examine the genotypes for the amino acid mutation (T929I) site involved in pyrethroid resistance. Findings showed the presence of thelytokous heterozygotes under the circumstance that most arrhenotokous insects are resistant homozygotes, and many thelytokous insects are susceptible homozygotes. These results suggest that, in the field, genetic exchange occurs between insects through of both reproductive modes. A survey of the genotypes for the other amino acid mutations using nucleotide sequencing showed a decline of insects with an M918T and L1014F pair and an increase of insects with M918L. These results suggest the evolutional progression of amino acid mutations associated with pyrethroid resistance in T. tabaci.

Search for genes gained by horizontal gene transfer in an entomopoxvirus, with special reference to the analysis of the transfer of an ABC transporter gene.

Although it is generally believed that large DNA viruses capture genes by horizontal gene transfer (HGT), the detailed manner of such transfer has not been fully elucidated. Here, we searched for genes in the coleopteran entomopoxvirus (EV) Anomala cuprea entomopoxvirus (ACEV) that might have been gained by ACEV by HGT. We classified the potential source organisms for HGT into three categories: the host A. cuprea; other organisms, including viruses unrelated to EVs; and organisms with uncertain host attribution. Of the open reading frames (ORFs) of the ACEV genome, 2.1 % were suggested to have been gained from the host by ACEV or its recent ancestor via HGT; 8.7 % were possibly from organisms other than the host, and 3.7 % were possibly from the third category of organisms via HGT. The analysis showed that ACEV contains some interesting ORFs obtained by HGT, including a large ATP-binding cassette protein (ABC transporter) ORF and a tenascin ORF (IDs ACV025 and ACV123, respectively). We then performed a detailed analysis of the HGT of the ACEV large ABC transporter ORF-the largest of the ACEV ORFs. mRNA sequences obtained by RNA-seq from fat bodies-sites of ACEV replication-and midgut tissues-sites of initial infection-of the virus's host A. cuprea larvae were subjected to BLAST analysis. One type of ABC transporter ORF from the fat bodies and two types from the midgut tissues, one of which was identical to that in the fat bodies, had the greatest identity to the ABC transporter ORF of ACEV. The two types from the host had high levels of identity to each other (approximately 95 % nucleotide sequence identity), strongly suggesting that the host ABC transporter group consisting of the two types was the origin of ACV025. We then determined the sequence (12,381 bp) containing a full-length gene of the A. cuprea ABC transporter. It turned out to be a transcription template for the abovementioned mRNA found in both tissues. In addition, we determined a large part (ca. 6.9 kb) of the template sequence for the mRNA found only in the midgut tissues. The results showed that the ACEV ABC transporter ORF is missing parts corresponding to introns of the host ABC transporter genes, indicating that the ORF was likely acquired by HGT in the form of mRNA. The presence of definite duplicated sequences adjacent to the ACEV ABC transporter genes-a sign of LINE-1 retrotransposon-mediated HGT-was not observed. An approximately 2-month ACV025 transcription experiment suggested that the transporter sequence is presumed to be continuously functional. The amino acid sequence of ACV025 suggests that its product might function in the regulation of phosphatide in the host-cell membranes.

Time-course transcriptome data of silk glands in day 0-7 last-instar larvae of Bombyx mori (w1 pnd strain).

Time-course transcriptome expression data were constructed for four parts of the silk gland (anterior, middle, and posterior parts of the middle silk gland, along with the posterior silk gland) in the domestic silkworm, Bombyx mori, from days 0 to 7 of the last-instar larvae. For sample preparation, silk glands were extracted from one female and one male larva every 24 hours accurately after the fourth ecdysis. The reliability of these transcriptome data was confirmed by comparing the transcripts per million (TPM) values of the silk gene and quantitative reverse transcription PCR results. Hierarchical cluster analysis results supported the reliability of transcriptome data. These data are likely to contribute to the progress in molecular biology and genetic research using B. mori, such as elucidating the mechanism underlying the massive production of silk proteins, conducting entomological research using a meta-analysis as a model for lepidopteran insect species, and exploring medical research using B. mori as a model for disease species by utilising transcriptome data.

Evolutionary origin and distribution of amino acid mutations associated with resistance to sodium channel modulators in onion thrips, Thrips tabaci.

In onion thrips Thrips tabaci, reduced sensitivity of the sodium channel caused by several sodium channel mutations have been correlated with pyrethroid resistance. For this study, using mitochondrial cytochrome c oxidase subunit I gene sequences, we examined the phylogenetic relation among a total of 52 thelytokous and arrhenotokous strains with different genotypes of the sodium channel mutations. Then, we used flow cytometry to estimate their ploidy. Results showed that the strains are divisible into three groups: diploid thelytoky, triploid thelytoky, and diploid arrhenotoky. Using 23 whole genome resequencing data obtained from 20 strains out of 52, we examined their genetic relation further using principal component analysis, admixture analysis, and a fixation index. Results showed that diploid and triploid thelytokous groups are further classifiable into two based on the sodium channel mutations harbored by the respective group members (strains). The greatest genetic divergence was observed between thelytokous and arrhenotokous groups with a pair of T929I and K1774N. Nevertheless, they shared a genomic region with virtually no polymorphism around the sodium channel gene loci, suggesting a hard selective sweep. Based on these findings, we discuss the evolutionary origin and distribution of the sodium channel mutations in T. tabaci.

KONAGAbase: a genomic and transcriptomic database for the diamondback moth, Plutella xylostella.

BACKGROUND: The diamondback moth (DBM), Plutella xylostella, is one of the most harmful insect pests for crucifer crops worldwide. DBM has rapidly evolved high resistance to most conventional insecticides such as pyrethroids, organophosphates, fipronil, spinosad, Bacillus thuringiensis, and diamides. Therefore, it is important to develop genomic and transcriptomic DBM resources for analysis of genes related to insecticide resistance, both to clarify the mechanism of resistance of DBM and to facilitate the development of insecticides with a novel mode of action for more effective and environmentally less harmful insecticide rotation. To contribute to this goal, we developed KONAGAbase, a genomic and transcriptomic database for DBM (KONAGA is the Japanese word for DBM). DESCRIPTION: KONAGAbase provides (1) transcriptomic sequences of 37,340 ESTs/mRNAs and 147,370 RNA-seq contigs which were clustered and assembled into 84,570 unigenes (30,695 contigs, 50,548 pseudo singletons, and 3,327 singletons); and (2) genomic sequences of 88,530 WGS contigs with 246,244 degenerate contigs and 106,455 singletons from which 6,310 de novo identified repeat sequences and 34,890 predicted gene-coding sequences were extracted. The unigenes and predicted gene-coding sequences were clustered and 32,800 representative sequences were extracted as a comprehensive putative gene set. These sequences were annotated with BLAST descriptions, Gene Ontology (GO) terms, and Pfam descriptions, respectively. KONAGAbase contains rich graphical user interface (GUI)-based web interfaces for easy and efficient searching, browsing, and downloading sequences and annotation data. Five useful search interfaces consisting of BLAST search, keyword search, BLAST result-based search, GO tree-based search, and genome browser are provided. KONAGAbase is publicly available from our website (http://dbm.dna.affrc.go.jp/px/) through standard web browsers. CONCLUSIONS: KONAGAbase provides DBM comprehensive transcriptomic and draft genomic sequences with useful annotation information with easy-to-use web interfaces, which helps researchers to efficiently search for target sequences such as insect resistance-related genes. KONAGAbase will be continuously updated and additional genomic/transcriptomic resources and analysis tools will be provided for further efficient analysis of the mechanism of insecticide resistance and the development of effective insecticides with a novel mode of action for DBM.

Highly accurate genome assembly of an improved high-yielding silkworm strain, Nichi01.

The silkworm (Bombyx mori) is an important lepidopteran model insect and an industrial domestic animal traditionally used for silk production. Here, we report the genome assembly of an improved Japanese strain Nichi01, in which the cocoon yield is comparable to that of commercial silkworm strains. The integration of PacBio Sequel II long-read and ddRAD-seq-based high-density genetic linkage map achieved the highest quality genome assembly of silkworms to date; 22 of the 28 pseudomolecules contained telomeric repeats at both ends, and only four gaps were present in the assembly. A total of 452 Mbp of the assembly with an N50 of 16.614 Mbp covered 99.3% of the complete orthologs of the lepidopteran core genes. Although the genome sequence of Nichi01 and that of the previously reported low-yielding tropical strain p50T assured their accuracy in most regions, we corrected several regions, misassembled in p50T, in our assembly. A total of 18,397 proteins were predicted using over 95 Gb of mRNA-seq derived from 10 different organs, covering 96.9% of the complete orthologs of the lepidopteran core genes. The final assembly and annotation files are available in KAIKObase (https://kaikobase.dna.affrc.go.jp/index.html) along with a genome browser and BLAST searching service, which would facilitate further studies and the breeding of silkworms and other insects.

A male-killing gene encoded by a symbiotic virus of Drosophila.

In most eukaryotes, biparentally inherited nuclear genomes and maternally inherited cytoplasmic genomes have different evolutionary interests. Strongly female-biased sex ratios that are repeatedly observed in various arthropods often result from the male-specific lethality (male-killing) induced by maternally inherited symbiotic bacteria such as Spiroplasma and Wolbachia. However, despite some plausible case reports wherein viruses are raised as male-killers, it is not well understood how viruses, having much smaller genomes than bacteria, are capable of inducing male-killing. Here we show that a maternally inherited double-stranded RNA (dsRNA) virus belonging to the family Partitiviridae (designated DbMKPV1) induces male-killing in Drosophila. DbMKPV1 localizes in the cytoplasm and possesses only four genes, i.e., one gene in each of the four genomic segments (dsRNA1-dsRNA4), in contrast to ca. 1000 or more genes possessed by Spiroplasma or Wolbachia. We also show that a protein (designated PVMKp1; 330 amino acids in size), encoded by a gene on the dsRNA4 segment, is necessary and sufficient for inducing male-killing. Our results imply that male-killing genes can be easily acquired by symbiotic viruses through reassortment and that symbiotic viruses are hidden players in arthropod evolution. We anticipate that host-manipulating genes possessed by symbiotic viruses can be utilized for controlling arthropods.

ABC transporter subfamily B1 as a susceptibility determinant of Bombyx mori larvae to Cry1Ba, Cry1Ia and Cry9Da toxins.

ATP binding cassette (ABC) transporters are a diverse family of transmembrane proteins. Specific subfamily members expressed in the lepidopteran midgut can act as susceptibility determinants for several insecticidal Bt Cry proteins. However, the susceptibility determinants to many Cry toxins still remain unclear. Therefore, we knocked out a series of ABC transporters that are highly expressed in the midgut of Bombyx mori larvae by transcription activator-like effector nuclease (TALEN)-mediated gene editing, and the lineages that became resistant to Cry toxins were searched by toxin overlay bioassay. As a result, the B. mori ABC transporter subfamily B1 (BmABCB1) knockout lineage showed 19.17-fold resistance to Cry1Ba, 876.2-fold resistance to Cry1Ia, and 29.1-fold resistance to Cry9Da, suggesting that BmABCB1 is the determinant of susceptibility to these toxins. BmABCC2 and BmABCC3 have been shown to be susceptibility determinants based on their function as receptors. Therefore, we next heterologously expressed these ABC transporters in HEK293T cells and performed a cell swelling assay to examine whether these molecules could exert receptor functions. As a result, BmABCB1-expressing cells showed swelling response to Cry1Ia and Cry9Da, and cells expressing PxABCB1, which is the Plutella xylostella ortholog of BmABCB1, showed swelling for Cry1Ba, suggesting that ABCB1 is a susceptibility determinant by functioning as a receptor to these toxins. Furthermore, in order to clarify how high binding affinity is based on receptor function, we performed surface plasmon resonance analysis and found that each KD of Cry1Ba, Cry1Ia, and Cry9Da to BmABCB1 were 7.69 × 10-8 M, 2.19 × 10-9 M, and 4.17 × 10-6 M respectively.

Reference Genome Sequences of the Oriental Armyworm, Mythimna separata (Lepidoptera: Noctuidae).

Lepidopteran insects are an important group of animals, including those used as biochemical and physiological model species in the insect and silk industries as well as others that are major agricultural pests. Therefore, the genome sequences of several lepidopteran insects have been reported. The oriental armyworm, Mythimna separata, is an agricultural pest commonly used to study insect immune reactions and interactions with parasitoid wasps as hosts. To improve our understanding of these research topics, reference genome sequences were constructed in the present study. Using long-read and short-read sequence data, de novo assembly and polishing were performed and haplotigs were purged. Subsequently, gene predictions and functional annotations were performed. To search for orthologs of the Toll and Immune Deficiency (IMD) pathways and for C-type lectins, annotation data analysis, BLASTp, and Hummer scans were performed. The M. separata genome is 682 Mbp; its contig N50 was 2.7 Mbp, with 21,970 genes and 24,452 coding sites predicted. All orthologs of the core components of the Toll and IMD pathways and 105 C-type lectins were identified. These results suggest that the genome data were of sufficient quality for use as reference genome data and could contribute to promoting M. separata and lepidopteran research at the molecular and genome levels.

Co-Expression Network and Time-Course Expression Analyses to Identify Silk Protein Regulatory Factors in Bombyx mori.

Bombyx mori is an important economic insect and an animal model in pharmacomedical research. Although its physiology has been studied for many years, the mechanism via which silk protein genes are regulated remains unclear. In this study, we performed two-step expression screening, namely co-expression network and time-course expression analyses to screen silk protein regulation factors. A co-expression network analysis using RNA-seq data that were obtained from various tissues, including the silk glands of B. mori, was performed to identify novel silk protein regulatory factors. Overall, 91 transcription factors, including some known ones, were found to be co-expressed with silk protein genes. Furthermore, time-course expression analysis during the fifth instar larvae stage revealed that the expression pattern of 13 novel transcription factors was highly relevant to that of silk protein genes and their known regulatory factor genes. In particular, the expression peak of several transcription factors (TFs) was detected before the expression of silk protein genes peak. These results indicated that a larger number of genes than expected may be involved in silk protein regulation in B. mori. Functional analyses of function-unknown transcription factors should enhance our understanding of this system.

A feasibility trial of genomics-based diagnosis detecting insecticide resistance of the diamondback moth.

BACKGROUND: Insecticide resistance management has been key for crop protection for over 70 years and is increasingly important because the development of new active ingredients has decreased in recent years. By monitoring the development of resistance in a timely manner, we can effectively prolong insecticide efficacy. Genomic-based diagnosis can reliably predict resistance development if information on resistant mutations against major pesticides is available. Here, we developed a feasibility trial of genomics-based diagnosis of insecticide resistance in diamondback moth (Plutella xylostella) populations in Nagano Prefecture, Japan. Amplicon sequencing analyses using a next-generation sequencer (Illumina MiSeq) for major insecticides, including diamides, pyrethroids, Bacillus thuringiensis (Bt) toxin (Cry1Ac), organophosphates, and spinosyns, were conducted. RESULTS: Mutations related to the resistance of pyrethroids, organophosphates, and diamides (flubendiamide and chlorantraniliprole) prevailed, while those of a diamide (cyantraniliprole), Bt (Cry1Ac), and spinosyns were scanty, suggesting that they are still effective. The results of the genomics-based diagnosis were generally concordant with the results of bioassays. Resistance development tendencies were generally uniform across Nagano. CONCLUSION: An insecticide-resistance management campaign can be conducted in Nagano Prefecture with a quick genomic-based diagnosis in early spring while bioassay is the only option for monitoring resistances whose mutations are unavailable. Our study is the first step in the future management of insecticide resistance in all significant pests. © 2022 Society of Chemical Industry.

Characterization of E93 in neometabolous thrips Frankliniella occidentalis and Haplothrips brevitubus.

Insect metamorphosis into an adult occurs after the juvenile hormone (JH) titer decreases at the end of the juvenile stage. This generally coincides with decreased transcript levels of JH-response transcription factors Krüppel homolog 1 (Kr-h1) and broad (br), and increased transcript levels of the adult specifier E93. Thrips (Thysanoptera) develop through inactive and non-feeding stages referred to as "propupa" and "pupa", and this type of distinctive metamorphosis is called neometaboly. To understand the mechanisms of hormonal regulation in thrips metamorphosis, we previously analyzed the transcript levels of Kr-h1 and br in two thrips species, Frankliniella occidentalis (Thripidae) and Haplothrips brevitubus (Phlaeothripidae). In both species, the transcript levels of Kr-h1 and br decreased in the "propupal" and "pupal" stages, and their transcription was upregulated by exogenous JH mimic treatment. Here we analyzed the developmental profiles of E93 in these two thrips species. Quantitative RT-PCR revealed that E93 expression started to increase at the end of the larval stage in F. occidentalis and in the "propupal" stage of H. brevitubus, as Kr-h1 and br mRNA levels decreased. Treatment with an exogenous JH mimic at the onset of metamorphosis prevented pupal-adult transition and caused repression of E93. These results indicated that E93 is involved in adult differentiation after JH titer decreases at the end of the larval stage of thrips. By comparing the expression profiles of Kr-h1, br, and E93 among insect species, we propose that the "propupal" and "pupal" stages of thrips have some similarities with the holometabolous prepupal and pupal stages, respectively.

Reference Transcriptome Data in Silkworm Bombyx mori.

Herein, we performed RNA-seq analysis of ten major tissues/subparts of silkworm larvae. The sequences were mapped onto the reference genome assembly and the reference transcriptome data were successfully constructed. The reference data provided a nearly complete sequence for sericin-1, a major silk gene with a complex structure. We also markedly improved the gene model for other genes. The transcriptomic expression was investigated in each tissue and a number of transcripts were identified that were exclusively expressed in tissues such as the testis. Transcripts strongly expressed in the midgut formed tight genomic clusters, suggesting that they originated from tandem gene duplication. Transcriptional factor genes expressed in specific tissues or the silk gland subparts were also identified. We successfully constructed reference transcriptome data in the silkworm and found that a number of transcripts showed unique expression profiles. These results will facilitate basic studies on the silkworm and accelerate its applications, which will contribute to further advances in lepidopteran and entomological research as well as the practical use of these insects.

An update of KAIKObase, the silkworm genome database.

KAIKObase was established in 2009 as the genome database of the domesticated silkworm Bombyx mori. It provides several gene sets and genetic maps as well as genome annotation obtained from the sequencing project of the International Silkworm Genome Consortium in 2008. KAIKObase has been used widely for silkworm and insect studies even though there are some erroneous predicted genes due to misassembly and gaps in the genome. In 2019, we released a new silkworm genome assembly, showing improvements in gap closure and covering more and longer gene models. Therefore, there is a need to include new genome and new gene models to KAIKObase. In this article, we present the updated contents of KAIKObase and the methods to generate, integrate and analyze the data sets. Database URL: https://kaikobase.dna.affrc.go.jp.