KAIKObase: an integrated silkworm genome database and data mining tool

fruitless is sex-differentially spliced and is important for the courtship behavior and development of silkmoth Bombyx mori

Sexual dimorphisms of the brain play essential roles in successful reproduction. Silkmoth Bombyx mori exhibits extensive sexual differences in sexual behavior, as well as their morphology. Although the neural circuits that transmit information about sex pheromone in the male brain are extensively analyzed, the molecular mechanisms that regulate their development are still elusive. In the present study, we focused on the silkmoth ortholog of fruitless (fru) as a candidate gene that regulates sexual dimorphisms of the brain. fru transcripts were expressed from multiple promoters in various tissues, and brain-specific transcripts were sex-specifically spliced, in a manner similar to Drosophila. Interestingly, fru was highly expressed in the adult female brain and the male larval testis. Analysis of CRISPR/Cas9-mediated fru knockout strains revealed that fru plays important roles in survival during late larval and pupal stages, testis development, and adult sexual behavior. fru mutant males exhibited highly reduced levels of courtship and low copulation rate, indicating that fru plays significant roles in the sexual behavior of silkmoths, although it is not absolutely necessary for copulation. In the fru mutant males, sexually dimorphic pattern of the odorant receptor expression was impaired, possibly causing the defects in courtship behavior. These results provide important clues to elucidate the development of sexual dimorphisms of silkmoth brains, as well as the evolution of fruitless gene in insects.

Bombyx Vasa sequesters transposon mRNAs in nuage via phase separation requiring RNA binding and self-association

Bombyx Vasa (BmVasa) assembles non-membranous organelle, nuage or Vasa bodies, in germ cells, known as the center for Siwi-dependent transposon silencing and concomitant Ago3-piRISC biogenesis. However, details of the body assembly remain unclear. Here, we show that the N-terminal intrinsically disordered region (N-IDR) and RNA helicase domain of BmVasa are responsible for self-association and RNA binding, respectively, but N-IDR is also required for full RNA-binding activity. Both domains are essential for Vasa body assembly in vivo and droplet formation in vitro via phase separation. FAST-iCLIP reveals that BmVasa preferentially binds transposon mRNAs. Loss of Siwi function derepresses transposons but has marginal effects on BmVasa-RNA binding. This study shows that BmVasa assembles nuage by phase separation via its ability to self-associate and bind newly exported transposon mRNAs. This unique property of BmVasa allows transposon mRNAs to be sequestered and enriched in nuage, resulting in effective Siwi-dependent transposon repression and Ago3-piRISC biogenesis.

Global Profiling of Genes Expressed in the Silk Glands of Philippine-Reared Mulberry Silkworms (Bombyx mori)

RNA sequencing was used to assemble transcriptome data for Philippine-reared silkworm and compare gene expression profiles of strains reared in high- and low-temperature environments. RNA was isolated from the silk glands of fifth instar larvae and mRNA-enriched libraries were sequenced using Illumina NextSeq 500. Transcriptome reads were assembled using reference-based and de novo assemblers, and assemblies were evaluated using different metrics for transcriptome quality, including the read mapping rate, E90N50, RSEM-eval, and the presence of single-copy orthologs. All transcriptome assemblies were able to reconstruct >40,000 transcripts. Differential expression analysis found 476 differentially expressed genes (DEGs; 222 upregulated, 254 downregulated) in strains reared in different temperatures. Among the top DEGs were myrosinase, heat shock proteins, serine protease inhibitors, dehydrogenases, and regulators of the juvenile hormone. Validation of some of the top DEGs using qPCR supported the findings of the in silico analysis. GO term enrichment analysis reveals an overrepresentation of GO terms related to nucleotide metabolism and biosynthesis, lipid and carbohydrate metabolic processes, regulation of transcription, nucleotide binding, protein binding, metal binding, catalytic activity, oxidoreductase activity, and hydrolase activity. The data provided here will serve as a resource for improving local strains and increasing silk production of Philippine-reared B. mori strains.

Bibliometric Analysis of Trends in Mulberry and Silkworm Research on the Production of Silk and Its By-Products

Simple Summary

Over the past two decades scientific research on sericulture, the agricultural activity of silk production, generated a great number of outputs in the form of articles reported and classified by one of the most well-known and used database dealing with scientific literature. This occurrence demonstrates an increasing interest in this sector especially starting from 2000s; results presented in relevant papers showed their applicability even in fields apparently not related to silk production as commonly meant, like medicine, cosmetics, and engineering. To understand how sericulture has been transcending its usual boundaries, which are its current “hotspots”, and links with other fields of study, the authors propose a text-mining based analysis of the outputs of scientific research on sericulture and silk; the final goal is to establish “quantitative” indicators for researchers, entrepreneurs, and scholars.

Abstract

Traditionally, sericulture is meant as the agricultural activity of silk production, from mulberry (Morus sp.pl.) cultivation to silkworm (Bombyx mori L.) rearing. The aim of the present work is to analyze the trends and outputs of scientific research on sericulture-related topics during the last two decades, from 2000 to 2020. In this work the authors propose a text-mining analysis of the titles, abstracts and keywords of scientific articles focused on sericulture and available in the SCOPUS database considering the above-mentioned period of time; from this article collection, the 100 most recurrent terms were extracted and studied in detail. The number of publications per year in sericulture-related topics increased from 87 in 2000 to 363 in 2020 (+317%). The 100 most recurrent terms were then aggregated in clusters. The analysis shows how in the last period scientific research, besides the traditional themes of sericulture, also focused on alternative products obtainable from the sericultural practice, as fruits of mulberry trees (increment of +134% of the occurrences in the last five years) and chemical compounds as antioxidants (+233% of occurrences), phenolics (+330% of occurrences) and flavonoids (+274% of occurrences). From these considerations, the authors can state how sericulture is an active and multidisciplinary research field.

Whole genome sequencing of spotted stem borer, Chilo partellus, reveals multiple genes encoding enzymes for detoxification of insecticides

Spotted stem borer, Chilo partellus, is the most important constraint for increasing the production and productivity of maize and sorghum, the two major coarse cereals in Asia and Africa. The levels of resistance to this pest in the cultivated germplasm are low to moderate, and hence, farmers have to use insecticides for effective control of this pest. However, there is no information on the detoxification mechanisms in C. partellus, which is one of the constraints for deployment of appropriate insecticides to control this pest. The ability to detoxify insecticides varies across insect populations, and hence, we sequenced different populations of C. partellus to identify and understand detoxification mechanisms to devise appropriate strategies for deployment of different insecticides for controlling this pest. Larval samples were sequenced from three different cohorts of C. partellus using the Illumina HiSeq 2500 platform. The data were subjected to identify putative genes that are involved in detoxification on insecticides in our cohort insect species. These studies resulted in identification of 64 cytochrome P450 genes (CYP450s), and 36 glutathione S-transferases genes (GSTs) encoding metabolic detoxification enzymes, primarily responsible for xenobiotic metabolism in insects. A total of 183 circadian genes with > 80% homolog and 11 olfactory receptor genes that mediate chemical cues were found in the C. partellus genome. Also, target receptors related to insecticide action, 4 acetylcholinesterase (AChE), 14 γ-aminobutyric acid (GABA), and 15 nicotinic acetylcholine (nAChR) receptors were detected. This is the first report of whole genome sequencing of C. partellus useful for understanding mode of action of different insecticides, and mechanisms of detoxification and designing target-specific insecticides to develop appropriate strategies to control C. partellus for sustainable crop production.

An Expanded Survey of the Moth PBP/GOBP Clade in Bombyx mori: New Insight into Expression and Functional Roles

We studied the expression profile and ontogeny (from the egg stage through the larval stages and pupal stages, to the elderly adult age) of four OBPs from the silkworm moth Bombyx mori. We first showed that male responsiveness to female sex pheromone in the silkworm moth B. mori does not depend on age variation; whereas the expression of BmorPBP1, BmorPBP2, BmorGOBP1, and BmorGOBP2 varies with age. The expression profile analysis revealed that the studied OBPs are expressed in non-olfactory tissues at different developmental stages. In addition, we tested the effect of insecticide exposure on the expression of the four OBPs studied. Exposure to a toxic macrolide insecticide endectocide molecule (abamectin) led to the modulated expression of all four genes in different tissues. The higher expression of OBPs was detected in metabolic tissues, such as the thorax, gut, and fat body. All these data strongly suggest some alternative functions for these proteins other than olfaction. Finally, we carried out ligand docking studies and reported that PBP1 and GOBP2 have the capacity of binding vitamin K1 and multiple different vitamins.

Identification and function of ETH receptor networks in the silkworm Bombyx mori

Insect ecdysis triggering hormones (ETHs) released from endocrine Inka cells act on specific neurons in the central nervous system (CNS) to activate the ecdysis sequence. These primary target neurons express distinct splicing variants of ETH receptor (ETHR-A or ETHR-B). Here, we characterized both ETHR subtypes in the moth Bombyx mori in vitro and mapped spatial and temporal distribution of their expression within the CNS and peripheral organs. In the CNS, we detected non-overlapping expression patterns of each receptor isoform which showed dramatic changes during metamorphosis. Most ETHR-A and a few ETHR-B neurons produce multiple neuropeptides which are downstream signals for the initiation or termination of various phases during the ecdysis sequence. We also described novel roles of different neuropeptides during these processes. Careful examination of peripheral organs revealed ETHRs expression in specific cells of the frontal ganglion (FG), corpora allata (CA), H-organ and Malpighian tubules prior to each ecdysis. These data indicate that PETH and ETH are multifunctional hormones that act via ETHR-A and ETHR-B to control various functions during the entire development—the ecdysis sequence and associated behaviors by the CNS and FG, JH synthesis by the CA, and possible activity of the H-organ and Malpighian tubules.

Host Response against Virus Infection in an Insect: Bidensovirus Infection Effect on Silkworm (Bombyx mori)

Silk cocoons obtained from silkworms are the primary source of commercial silk, making the silkworm an economically important insect. However, the silk industry suffers significant losses due to various virus infections. Bombyx mori bidensovirus (BmBDV) is one of the pathogens that cause flacherie disease in silkworms. Most silkworm strains die after BmBDV infection. However, certain silkworm strains show resistance to the virus, which is determined by a single recessive gene, nsd-2. The +^nsd-2 gene (allele of nsd-2; the susceptibility gene) encodes a putative amino acid transporter expressed only in the insect’s midgut, where BmBDV can infect, suggesting that this membrane protein may function as a receptor for BmBDV. Interestingly, the expression analysis revealed no changes in the +^nsd-2 gene expression levels in virus-uninfected silkworms, whereas the gene expression drastically decreased in the virus-infected silkworm. This condition indicates that the host factor’s expression, the putative virus receptor, is affected by BmBDV infection. It has recently been reported that the expression levels of some host genes encoding cuticle, antioxidant, and immune response-related proteins were significantly regulated by BmBDV infection. In this review, we discuss the host response against virus infection based on our knowledge and long-term research experience in this field.

Bombyx mori β1,4-N-acetylgalactosaminyltransferase possesses relaxed donor substrate specificity in N-glycan synthesis

N-Glycosylation is one of the most important post-translational protein modifications in eukaryotic cells. Although more than 200 N-glycogenes contributing to N-glycan biosynthesis have been identified and characterized, the information on insect N-glycosylation is still limited. Here, focusing on insect N-glycosylation, we characterized Bombyx mori N-acetylgalactosaminyltransferase (BmGalNAcT) participating in complex N-glycan biosynthesis in mammals. BmGalNAcT localized at the Golgi and was ubiquitously expressed in every organ and in the developmental stage of the middle silk gland of fifth instar larvae. Analysis of recombinant BmGalNAcT expressed in Sf9 cells showed that BmGalNAcT transferred GalNAc to non-reducing terminals of GlcNAcβ1,2-R with β1,4-linkage. In addition, BmGalNAcT mediated transfer of galactose and N-acetylglucosamine residues but not transfer of either glucose or glucuronic acid from the UDP-sugar donor substrate to the N-glycan. Despite this tri-functional sugar transfer activity, however, most of the endogenous glycoproteins of insect cells were present without GalNAc, Gal, or GlcNAc residues at the non-reducing terminal of β1,2-GlcNAc residue(s). Moreover, overexpression of BmGalNAcT in insect cells had no effect on N-acetylgalactosaminylation, galactosylation, or N-acetylglucosaminylation of the major N-glycan during biosynthesis. These results suggested that B. mori has a novel multifunctional glycosyltransferase, but the N-glycosylation is highly and strictly regulated by the endogenous N-glycosylation machineries.

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.

SilkDB 3.0: visualizing and exploring multiple levels of data for silkworm

SilkDB is an open-accessibility database and powerful platform that provides comprehensive information on the silkworm (Bombyx mori) genome. Since SilkDB 2.0 was released 10 years ago, vast quantities of data about multiple aspects of the silkworm have been generated, including genome, transcriptome, Hi-C and pangenome. To visualize data at these different biological levels, we present SilkDB 3.0 (https://silkdb.bioinfotoolkits.net), a visual analytic tool for exploring silkworm data through an interactive user interface. The database contains a high-quality chromosome-level assembly of the silkworm genome, and its coding sequences and gene sets are more accurate than those in the previous version. SilkDB 3.0 provides a view of the information for each gene at the levels of sequence, protein structure, gene family, orthology, synteny, genome organization and gives access to gene expression information, genetic variation and genome interaction map. A set of visualization tools are available to display the abundant information in the above datasets. With an improved interactive user interface for the integration of large data sets, the updated SilkDB 3.0 database will be a valuable resource for the silkworm and insect research community.

Insect genomes: progress and challenges

In the wake of constant improvements in sequencing technologies, numerous insect genomes have been sequenced. Currently, 1219 insect genome-sequencing projects have been registered with the National Center for Biotechnology Information, including 401 that have genome assemblies and 155 with an official gene set of annotated protein-coding genes. Comparative genomics analysis showed that the expansion or contraction of gene families was associated with well-studied physiological traits such as immune system, metabolic detoxification, parasitism and polyphagy in insects. Here, we summarize the progress of insect genome sequencing, with an emphasis on how this impacts research on pest control. We begin with a brief introduction to the basic concepts of genome assembly, annotation and metrics for evaluating the quality of draft assemblies. We then provide an overview of genome information for numerous insect species, highlighting examples from prominent model organisms, agricultural pests and disease vectors. We also introduce the major insect genome databases. The increasing availability of insect genomic resources is beneficial for developing alternative pest control methods. However, many opportunities remain for developing data-mining tools that make maximal use of the available insect genome resources. Although rapid progress has been achieved, many challenges remain in the field of insect genomics.

BmBlimp-1 gene encoding a C2H2 zinc finger protein is required for wing development in the silkworm Bombyx mori

Cys2-His2 zinc finger (C2H2-ZF) proteins represent the most common class of transcription factors. These factors have great potential for the management of developmental progression by regulating the specific spatiotemporal expression of genes. In this study, we cloned one C2H2-ZF protein gene of Bombyx mori, BGIBMGA000319, that is orthologous to B-lymphocyte-induced maturation protein-1 (Blimp-1); we thus named it as Bombyx mori Blimp-1 (BmBlimp-1). In the silkworm, the BmBlimp-1 gene is specifically upregulated during day 2 of the pupal to adult stage and is highly expressed in wing discs on day 3 of the pupa. Knockdown of its expression level in the pupal stage results in a crumpled-winged silkworm moth. Using the predicted DNA-binding sequences of BmBlimp-1 to search the silkworm genome to screen target genes of BmBlimp-1, 7049 genes were identified to have at least one binding site of BmBlimp-1 on their 1 kb upstream and downstream genome regions. Comparisons of those genes with a reported pupal wing disc transcriptome data resulted in 4065 overlapping genes being retrieved. GO enrichment analysis of the overlapping genes showed that most of the genes were enriched in the binding term. Combining functional annotation and real-time quantitative PCR, 15 genes were identified as the candidate target genes of BmBlimp-1, including several wing cuticular protein genes, chitin synthase A, and wing disc development genes, such as Wnt1, cubitus interruptus (ci) and engrailed (en). Moreover, the amino acid sequence of the zinc finger motif of Blimp-1 gene was highly conserved among the 15 insect species. We propose that BmBlimp-1 is an important regulatory factor in silkworm wing development.

Influence of host plant, geography and pheromone strain on genomic differentiation in sympatric populations of Ostrinia nubilalis

Patterns of mating for the European corn borer (Ostrinia nubilalis) moth depend in part on variation in sex-pheromone blend. The ratio of (E)-11- and (Z)-11-tetradecenyl acetate (E11- and Z11-14:OAc) in the pheromone blend that females produce and males respond to differs between strains of O. nubilalis. Populations also vary in female oviposition preference for and larval performance on maize (C4) and nonmaize (C3) host plants. The relative contributions of sexual and ecological trait variation to the genetic structure of O. nubilalis remains unknown. Host-plant use (¹³ C/¹⁴ C ratios) and genetic differentiation were estimated among sympatric E and Z pheromone strain O. nubilalis males collected in sex-pheromone baited traps at 12 locations in Pennsylvania and New York between 2007 and 2010. Among genotypes at 65 single nucleotide polymorphism marker loci, variance at a position in the pheromone gland fatty acyl-reductase (pgfar) gene at the locus responsible for determining female pheromone ratio (Pher) explained 64% of the total genetic differentiation between males attracted to different pheromones (male response, Resp), providing evidence of sexual inter-selection at these unlinked loci. Principal coordinate, Bayesian clustering, and distance-based redundancy analysis (dbRDA) demonstrate that host plant history or geography does not significantly contribute to population variation or differentiation among males. In contrast, these analyses indicate that pheromone response and pgfar-defined strain contribute significantly to population genetic differentiation. This study suggests that behavioural divergence probably plays a larger role in driving genetic variation compared to host plant-defined ecological adaptation.

Genome-wide identification and expression profiling of the C2H2-type zinc finger protein genes in the silkworm Bombyx mori

Cys2-His2 zinc finger (C2H2-ZF) proteins comprise the largest class of putative eukaryotic transcription factors. The zinc finger motif array is highly divergent, indicating that most proteins will have distinctive binding sites and perform different functions. However, the binding sites and functions of the majority of C2H2-ZF proteins remain unknown. In this study, we identified 327 C2H2-ZF protein genes in the silkworm, 290 in the monarch butterfly, 243 in the fruit fly, 107 in elegans, 673 in mouse, and 1,082 in human. The C2H2-ZF protein genes of the silkworm were classified into three main grouping clades according to a phylogenetic classification, and 312 of these genes could be mapped onto 27 chromosomes. Most silkworm C2H2-ZF protein genes exhibited specific expression in larval tissues. Furthermore, several C2H2-ZF protein genes had sex-specific expression during metamorphosis. In addition, we found that some C2H2-ZF protein genes are involved in metamorphosis and female reproduction by using expression clustering and gene annotation analysis. Among them, five genes were selected, BGIBMGA002091 (CTCF), BGIBMGA006492 (fru), BGIBMGA006230 (wor), BGIBMGA004640 (lola), and BIGBMGA004569, for quantitative real-time PCR analysis from larvae to adult ovaries. The results showed that the five genes had different expression patterns in ovaries, among which BGIBMGA002091 (CTCF) gene expression level was the highest, and its expression level increased rapidly in late pupae and adult stages. These findings provide a basis for further investigation of the functions of C2H2-ZF protein genes in the silkworm, and the results offer clues for further research into the development of metamorphosis and female reproduction in the silkworm.

Unprecedented reorganization of holocentric chromosomes provides insights into the enigma of lepidopteran chromosome evolution

Chromosome evolution presents an enigma in the mega-diverse Lepidoptera. Most species exhibit constrained chromosome evolution with nearly identical haploid chromosome counts and chromosome-level gene collinearity among species more than 140 million years divergent. However, a few species possess radically inflated chromosomal counts due to extensive fission and fusion events. To address this enigma of constraint in the face of an exceptional ability to change, we investigated an unprecedented reorganization of the standard lepidopteran chromosome structure in the green-veined white butterfly (Pieris napi). We find that gene content in P. napi has been extensively rearranged in large collinear blocks, which until now have been masked by a haploid chromosome number close to the lepidopteran average. We observe that ancient chromosome ends have been maintained and collinear blocks are enriched for functionally related genes suggesting both a mechanism and a possible role for selection in determining the boundaries of these genome-wide rearrangements.

Genome-Wide Analysis and Hormone Regulation of Chitin Deacetylases in Silkworm

Chitin deacetylases (CDAs) are a group of enzymes involved in chitin metabolism in insects; they play a critical role in molting, pupation, and the modification of chitin. In this study, we identified several CDAs in the silkworm, Bombyx mori (BmCDA), and investigated the effect of various hormones on their expression in B. mori larvae and embryo cell lines (BmE). Eight genes encoding BmCDAs were identified in the silkworm genome. They showed different expression patterns in different tissues, and were classified into three types based on where they were expressed: the exoskeleton, digestive organs, and genital organs. Moreover, we found that some BmCDAs showed upregulated expression during the molting period, especially during the fourth molting period in larvae. We also verified that the expression of BmCDA1-6 was upregulated by treatment with 20-hydroxyecdysone not only in larvae, but also in BmE cells. Interestingly, juvenile hormone analog treatment also upregulated the expression of some BmCDAs. The overexpression of several transcription factors revealed that the POU transcription factor POUM2 may play a major role in the regulation of BmCDA expression. Finally, the silencing of BmCDA1 and BmCDA2 did not lead to abnormal phenotypes or death, but may have led to delays in silkworm pupation. These results provide important information about lepidopteran insects in terms of chitin deacetylases and the regulation of their expression.

Silkworm, Bombyx mori, as an alternative model organism in toxicological research

Silkworms, Bombyx mori, are a promising model animal in health safety and environmental pollution assessment due to their sensitivity to chemical compounds like pesticides, drugs, and heavy metals, in addition to other features like their low cost and body characteristics and their full genome sequencing. In this review, we summarize the silkworm advantages as a model organism in toxicological research. Graphical abstract ᅟ.

Genome editing reveals the function of Yorkie during the embryonic and early larval development in silkworm, Bombyx mori

As a transcriptional coactivator, Yorkie (Yki) is a major downstream target of the Hippo signalling pathway to regulate the organ size during animal development and regeneration. Previous microarray analysis in the silkworm, Bombyx mori, has shown that genes associated with the Hippo pathway were primarily expressed in gonads and imaginal discs. The RNA-interference-mediated silencing of Yki at the early wandering stage delayed B. mori development and ovary maturation, whereas baculovirus-mediated overexpression at the late larval instar facilitated organ growth and accelerated metamorphosis. Here, we employed CRISPR/Cas9-mediated mutagenesis to investigate the function of Yki in B. mori (BmYki) at the embryonic and early larval stages. Knocking out of BmYki led to reduced body size, moulting defects and, eventually, larval lethality. Sequence analysis of CRISPR/Cas9 mutants exhibited an array of deletions in BmYki. As a critical downstream effector of the Hippo kinase cassette, silencing of BmYki at the embryonic stage is indispensable and the consequence is lethal. Given that the Hippo signalling pathway is evolutionarily conserved, Yki has the potential to be a novel molecular target for genetic-based pest management practices.

Decrease in the expression level of the gene encoding the putative Bombyx mori bidensovirus receptor during virus infection

Bombyx mori bidensovirus (BmBDV) is a pathogen that replicates only in the midgut columnar cells of silkworms, causing fatal disease. Resistance to BmBDV, which does not depend on the viral dose, is determined by a single gene, nsd-2 (resistance gene). Previously, we identified nsd-2 by positional cloning using B. mori genome information and found that this gene encodes a putative amino acid transporter that may function as a receptor for BmBDV. In this study, to understand the relationship between BmBDV and the putative virus receptor, we performed expression analysis of +^nsd-2 (allele of nsd-2; susceptibility gene) after virus infection. Quantitative RT-PCR analysis using total RNA isolated from the midgut of an uninfected and a virus-infected silkworm revealed no change in the expression levels of +^nsd-2 in the uninfected silkworm, whereas the expression levels of +^nsd-2 drastically decreased in the virus-infected silkworm. Moreover, comparison of the expression pattern between the BmBDV-derived transcript and +^nsd-2 revealed that the expression level of +^nsd-2 decreased with an increase in the virus-derived transcript. In addition, expression analysis of 26 genes encoding other transporters in the midgut demonstrated that the expression levels of three other genes also decreased similarly to the decrease of the expression levels of +^nsd-2 after virus infection. Thus, our results suggest that some transporters, including +^nsd-2, are affected by BmBDV infection.

Non-Pleiotropic Coupling of Daily and Seasonal Temporal Isolation in the European Corn Borer

Speciation often involves the coupling of multiple isolating barriers to produce reproductive isolation, but how coupling is generated among different premating barriers is unknown. We measure the degree of coupling between the daily mating time and seasonal mating time between strains of European corn borer (Ostrinia nubilalis) and evaluate the hypothesis that the coupling of different forms of allochrony is due to a shared genetic architecture, involving genes with pleiotropic effects on both timing phenotypes. We measure differences in gene expression at peak mating times and compare these genes to previously identified candidates that are associated with changes in seasonal mating time between the corn borer strains. We find that the E strain, which mates earlier in the season, also mates 2.7 h earlier in the night than the Z strain. Earlier daily mating is correlated with the differences in expression of the circadian clock genes cycle, slimb, and vrille. However, different circadian clock genes associate with daily and seasonal timing, suggesting that the coupling of timing traits is maintained by natural selection rather than pleiotropy. Juvenile hormone gene expression was associated with both types of timing, suggesting that circadian genes activate common downstream modules that may impose constraint on future evolution of these traits.

Deficiency of a pyrroline-5-carboxylate reductase produces the yellowish green cocoon 'Ryokuken' of the silkworm, Bombyx mori

The silkworm cocoon colour has attracted researchers involved in genetics, physiology and ecology for a long time. 'Ryokuken' cocoons are yellowish green in colour due to unusual flavonoids, prolinylflavonols, while 'Sasamayu' cocoons are light green and contain only simple flavonol glucosides. We found a novel gene associated with the cocoon colour change resulting from a change in flavonoid composition and named it Lg (light green cocoon). In the middle silk glands of the + ^Lg /+ ^Lg larvae, 1-pyrroline-5-carboxylic acid (P5C) was found to accumulate due to a decrease in the activity of pyrroline-5-carboxylate reductase (P5CR), an enzyme reducing P5C to proline. Sequence analysis of BmP5CR1, the candidate gene for Lg, revealed a 1.9 kb insertion and a 4 bp deletion within the 1st intron, a 97 bp deletion within the 4th intron, and a > 300 bp insertion within the 3'-UTR, in addition to two amino acid changes on exons 3 and 4 in + ^Lg /+ ^Lg compared to Lg/Lg. Decreased expression of BmP5CR1 was observed in all of the investigated tissues, including the middle silk glands in + ^Lg /+ ^Lg , which was probably caused by structural changes in the intronic regions of BmP5CR1. Furthermore, a BmP5CR1 knockout strain exhibited a yellowish green cocoon with the formation of prolinylflavonols. These results indicate that the yellowish green cocoon is produced by a BmP5CR1 deficiency. To our knowledge, this is the first report showing that the defect of an enzyme associated with intermediate metabolism promotes the conjugation of phytochemicals derived from foods with endogenously accumulating metabolites in animal tissues.

Intrinsic antimicrobial properties of silk spun by genetically modified silkworm strains

The domesticated silkworm, Bombyx mori, is a fundamental insect for silk industry. Silk is obtained from cocoons, protective envelopes produced during pupation and composed of single raw silk filaments secreted by the insect silk glands. Currently, silk is used as a textile fibre and to produce new materials for technical and biomedical applications. To enhance the use of both fabrics and silk-based materials, great efforts to obtain silk with antimicrobial properties have been made. In particular, a convincing approach is represented by the enrichment of the textile fibre with antimicrobial peptides, the main effectors of the innate immunity. To this aim, silkworm-based transgenic techniques appear to be cost-effective strategies to obtain cocoons in which antimicrobial peptides are integrated among the silk proteins. Recently, cocoons transgenic for a recombinant silk protein conjugated to the silkworm Cecropin B antimicrobial peptide were obtained and showed enhanced antibacterial properties (Li et al. in Mol Biol Rep 42:19-25, https://doi.org/10.1007/s11033-014-3735-z , 2015a). In this work we used the piggyBac-mediated germline transformation to generate several transgenic B. mori lines able to overexpress Cecropin B or Moricin antimicrobial peptides at the level of the silk gland. The derived cocoons were characterised by increased antimicrobial properties and the resulting silk fibre was able to inhibit the bacterial growth of the Gram-negative Escherichia coli. Our results suggest that the generation of silkworm overexpressing unconjugated antimicrobial peptides in the silk gland might represent an additional strategy to obtain antimicrobial peptide-enriched silk, for the production of new silk-based materials.

Identification of circular RNA in the Bombyx mori silk gland

Bombyx mori is an economically important holometabolous lepidopteran insect. In B. mori endogenous noncoding RNAs such as microRNAs (miRNAs) and Piwi-interacting RNAs play crucial biological functions in metamorphosis and sex determination. In addition, circular RNAs (circRNAs) have been recently identified as noncoding RNAs in most common model organisms and show potential as gene regulators. However, to date, there have been few studies on the circRNAs present in the B. mori genome conducted to date. Here, we identified 3916 circRNAs by deep circular transcriptome sequencing using the silk gland of B. mori. 3155 circRNAs were found to be derived from 1727 parental genes. The circRNAs displayed tissue-specific expression between the middle silk gland (MSG) and posterior silk gland (PSG), with 2532 and 880 being upregulated circRNAs in the MSG and PSG, respectively. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses showed that the parental genes from the MSG and PSG were generally annotated to similar categories and pathways. The interaction network of circRNAs and miRNAs showed that circRNAs might act as miRNA sponges or interact with miRNAs in some other way. Overall, the results revealed the complicated patterns of circRNAs in the B. mori silk gland providing a new angle from which to explore the mechanisms of complex gene regulation and efficient silk protein synthesis.

Metabolomics differences between silkworms (Bombyx mori) reared on fresh mulberry (Morus) leaves or artificial diets

Silkworms (Bombyx mori) reared on artificial diets have great potential applications in sericulture. However, the mechanisms underlying the enhancement of metabolic utilization by altering silkworm nutrition are unclear. The aim of this study was to investigate the mechanisms responsible for the poor development and low silk protein synthesis efficiency of silkworms fed artificial diets. After multi-generational selection of the ingestive behavior of silkworms to artificial diets, we obtained two strains, one of which developed well and another in which almost all its larvae starved to death on the artificial diets. Subsequently, we analyzed the metabolomics of larval hemolymph by gas chromatography/liquid chromatography–mass spectrometry, and the results showed that vitamins were in critically short supply, whereas the nitrogen metabolic end product of urea and uric acid were enriched substantially, in the hemolymph of the silkworms reared on the artificial diets. Meanwhile, amino acid metabolic disorders, as well as downregulation of carbohydrate metabolism, energy metabolism, and lipid metabolism, co-occurred. Furthermore, 10 male-dominant metabolites and 27 diet-related metabolites that differed between male and female silkworms were identified. These findings provide important insights into the regulation of silkworm metabolism and silk protein synthesis when silkworms adapt to an artificial diet.

RNA sequencing reveals a complete but an unconventional type of dosage compensation in the domestic silkworm Bombyx mori

Sex chromosomal dose difference between sexes is often normalized by a gene regulatory mechanism called dosage compensation (DC). Studies indicate that DC mechanisms are generally effective in XY rather than ZW systems. However, DC studies in lepidopterans (ZW system) gave bewildering results. In Manduca sexta, DC was complete and in Plodia interpunctella, it was incomplete. In Heliconius species, dosage was found to be partly incomplete. In domesticated silkmoth Bombyx mori, DC studies have yielded contradictory results thus far, showing incomplete DC based on microarray data and a possible existence of DC based on recent reanalysis of same data. In this study, analysis of B. mori sexed embryos (78, 96 and 120 h) and larval heads using RNA sequencing suggest an onset of DC at 120 h. The average Z-linked expression is substantially less than autosomes, and the male-biased Z-linked expression observed at initial stages (78 and 96 h) gets almost compensated at 120 h embryonic stage and perfectly compensated in heads. Based on these findings, we suggest a complete but an unconventional type of DC, which may be achieved by reduced Z-linked expression in males (ZZ). To our knowledge, this is the first next-generation sequencing report showing DC in B. mori, clarifying the previous contradictions.

Comparative transcriptomes analysis of the wing disc between two silkworm strains with different size of wings

Wings of Bombyx mori (B. mori) develop from the primordium, and different B. mori strains have different wing types. In order to identify the key factors influencing B. mori wing development, we chose strains P50 and U11, which are typical for normal wing and minute wing phenotypes, respectively. We dissected the wing disc on the 1^st-day of wandering stage (P50D1 and U11D1), 2^nd-day of wandering stage (P50D2 and U11D2), and 3^rd-day of wandering stage (P50D3 and U11D3). Subsequently, RNA-sequencing (RNA-Seq) was performed on both strains in order to construct their gene expression profiles. P50 exhibited 628 genes differentially expressed to U11, 324 up-regulated genes, and 304 down-regulated genes. Five enriched gene ontology (GO) terms were identified by GO enrichment analysis based on these differentially expressed genes (DEGs). KEGG enrichment analysis results showed that the DEGs were enriched in five pathways; of these, we identified three pathways related to the development of wings. The three pathways include amino sugar and nucleotide sugar metabolism pathway, proteasome signaling pathway, and the Hippo signaling pathway. The representative genes in the enrichment pathways were further verified by quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR). The RNA-Seq and qRT-PCR results were largely consistent with each other. Our results also revealed that the significantly different genes obtained in our study might be involved in the development of the size of B. mori wings. In addition, several KEGG enriched pathways might be involved in the regulation of the pathways of wing formation. These results provide a basis for further research of wing development in B. mori.

Genetic basis of allochronic differentiation in the fall armyworm

BACKGROUND

Very little is known on how changes in circadian rhythms evolve. The noctuid moth Spodoptera frugiperda (Lepidoptera: Noctuidae) consists of two strains that exhibit allochronic differentiation in their mating time, which acts as a premating isolation barrier between the strains. We investigated the genetic basis of the strain-specific timing differences to identify the molecular mechanisms of differentiation in circadian rhythms.

RESULTS

Through QTL analyses we identified one major Quantitative trait chromosome (QTC) underlying differentiation in circadian timing of mating activity. Using RADtags, we identified this QTC to be homologous to Bombyx mori C27, on which the clock gene vrille is located, which thus became the major candidate gene. In S. frugiperda, vrille showed strain-specific polymorphisms. Also, vrille expression differed significantly between the strains, with the rice-strain showing higher expression levels than the corn-strain. In addition, RT-qPCR experiments with the other main clock genes showed that pdp1, antagonist of vrille in the modulatory feedback loop of the circadian clock, showed higher expression levels in the rice-strain than in the corn-strain.

CONCLUSIONS

Together, our results indicate that the allochronic differentiation in the two strains of S. frugiperda is associated with differential transcription of vrille or a cis-acting gene close to vrille, which contributes to the evolution of prezygotic isolation in S. frugiperda.

Expression map of a complete set of gustatory receptor genes in chemosensory organs of Bombyx mori

Most lepidopteran species are herbivores, and interaction with host plants affects their gene expression and behavior as well as their genome evolution. Gustatory receptors (Grs) are expected to mediate host plant selection, feeding, oviposition and courtship behavior. However, due to their high diversity, sequence divergence and extremely low level of expression it has been difficult to identify precisely a complete set of Grs in Lepidoptera. By manual annotation and BAC sequencing, we improved annotation of 43 gene sequences compared with previously reported Grs in the most studied lepidopteran model, the silkworm, Bombyx mori, and identified 7 new tandem copies of BmGr30 on chromosome 7, bringing the total number of BmGrs to 76. Among these, we mapped 68 genes to chromosomes in a newly constructed chromosome distribution map and 8 genes to scaffolds; we also found new evidence for large clusters of BmGrs, especially from the bitter receptor family. RNA-seq analysis of diverse BmGr expression patterns in chemosensory organs of larvae and adults enabled us to draw a precise organ specific map of BmGr expression. Interestingly, most of the clustered genes were expressed in the same tissues and more than half of the genes were expressed in larval maxillae, larval thoracic legs and adult legs. For example, BmGr63 showed high expression levels in all organs in both larval and adult stages. By contrast, some genes showed expression limited to specific developmental stages or organs and tissues. BmGr19 was highly expressed in larval chemosensory organs (especially antennae and thoracic legs), the single exon genes BmGr53 and BmGr67 were expressed exclusively in larval tissues, the BmGr27-BmGr31 gene cluster on chr7 displayed a high expression level limited to adult legs and the candidate CO₂ receptor BmGr2 was highly expressed in adult antennae, where few other Grs were expressed. Transcriptional analysis of the Grs in B. mori provides a valuable new reference for finding genes involved in plant-insect interactions in Lepidoptera and establishing correlations between these genes and vital insect behaviors like host plant selection and courtship for mating.

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.

Gene expression analysis in the larval silk gland of the eri silkworm Samia ricini

Insects produce silk for a range of purposes. In the Lepidoptera, silk is utilized as a material for cocoon production and serves to protect larvae from adverse environmental conditions or predators. Species in the Saturniidae family produce an especially wide variety of cocoons, for example, large, golden colored cocoons and those with many small holes. Although gene expression in the silk gland of the domestic silkworm (Bombyx mori L.) has been extensively studied, considerably fewer investigations have focused on members of the saturniid family. Here, we established expression sequence tags from the silk gland of the eri silkworm (Samia ricini), a saturniid species, and used these to analyze gene expression. Although we identified the fibroin heavy chain gene in the established library, genes for other major silk proteins, such as fibroin light chain and fibrohexamerin, were absent. This finding is consistent with previous reports that these latter proteins are lacking in saturniid silk. Recently, a series of fibrohexamerin-like genes were identified in the Bombyx genome. We used this information to conduct a detailed analysis of the library established here. This analysis identified putative homologues of these genes. We also found several genes encoding small silk protein molecules that are also present in the silk of other Lepidoptera. Gene expression patterns were compared between eri and domestic silkworm, and both conserved and nonconserved expression patterns were identified for the tested genes. Such differential gene expression might be one of the major causes of the differences in silk properties between these species. We believe that our study can be of value as a basic catalogue for silk gland gene expression, which will yield to the further understanding of silk evolution.

Emerging developmental genetic model systems in holometabolous insects

The number of insect species that are amenable to functional genetic studies is growing rapidly and provides many new research opportunities in developmental and evolutionary biology. The holometabolous insects represent a disproportionate percentage of animal diversity and are thus well positioned to provide model species for a wide variety of developmental processes. Here we discuss emerging holometabolous models, and review some recent breakthroughs. For example, flies and midges were found to use structurally unrelated long-range pattern organizers, butterflies and moths revealed extensive pattern formation during oogenesis, new imaging possibilities in the flour beetle Tribolium castaneum showed how embryos break free of their extraembryonic membranes, and the complex genetics governing interspecies difference in head shape were revealed in Nasonia wasps.

A comprehensive view of the web-resources related to sericulture

Recent progress in the field of sequencing and analysis has led to a tremendous spike in data and the development of data science tools. One of the outcomes of this scientific progress is development of numerous databases which are gaining popularity in all disciplines of biology including sericulture. As economically important organism, silkworms are studied extensively for their numerous applications in the field of textiles, biomaterials, biomimetics, etc. Similarly, host plants, pests, pathogens, etc. are also being probed to understand the seri-resources more efficiently. These studies have led to the generation of numerous seri-related databases which are extremely helpful for the scientific community. In this article, we have reviewed all the available online resources on silkworm and its related organisms, including databases as well as informative websites. We have studied their basic features and impact on research through citation count analysis, finally discussing the role of emerging sequencing and analysis technologies in the field of seri-data science. As an outcome of this review, a web portal named SeriPort, has been created which will act as an index for the various sericulture-related databases and web resources available in cyberspace.Database URL: http://www.seriport.in/.

The potential role of sexual conflict and sexual selection in shaping the genomic distribution of Mito-nuclear genes

Mitochondrial interactions with the nuclear genome represent one of life's most important co-evolved mutualisms. In many organisms, mitochondria are maternally inherited, and in these cases, co-transmission between the mitochondrial and nuclear genes differs across different parts of the nuclear genome, with genes on the X chromosome having two-third probability of co-transmission, compared with one-half for genes on autosomes. These asymmetrical inheritance patterns of mitochondria and different parts of the nuclear genome have the potential to put certain gene combinations in inter-genomic co-adaptation or conflict. Previous work in mammals found strong evidence that the X chromosome has a dearth of genes that interact with the mitochondria (mito-nuclear genes), suggesting that inter-genomic conflict might drive genes off the X onto the autosomes for their male-beneficial effects. Here, we developed this idea to test coadaptation and conflict between mito-nuclear gene combinations across phylogenetically independent sex chromosomes on a far broader scale. We found that, in addition to therian mammals, only Caenorhabditis elegans showed an under-representation of mito-nuclear genes on the sex chromosomes. The remaining species studied showed no overall bias in their distribution of mito-nuclear genes. We discuss possible factors other than inter-genomic conflict that might drive the genomic distribution of mito-nuclear genes.

Morphological characterization and molecular mapping of an irradiation-induced Speckled mutant in the silkworm, Bombyx mori

Speckled (Spc), an X-ray-induced lethal mutant of Bombyx mori, exhibits a mosaic dark-brown-spotted larval epidermis in both sexes and egg-laying problems only in females. Here, we report the morphological characterization and molecular mapping of the Spc mutant. Morphological investigations revealed that the epidermal ultrastructure of the small, dark-brown spots was more dense than that of the white regions in both Spc/+ mutants and wild type, and that the lethality of the Spc/Spc mutants occurred during early embryogenesis. Furthermore, the ovarioles and ovipositor were disconnected in approximately 85.5% of Spc/+ females, a further 2.5% had a connection between the ovarioles and ovipositor that was too narrow to lay eggs. The remaining females showed a normal connection similar to that of the wild type. We successfully narrowed down the location of the Spc mutation to a region on chromosome 4 that was ∼1041 kb long. Gene-prediction analysis identified 25 candidate genes in this region. Chromosome structure analysis indicated that a ∼305 kb deletion was included in the mapping region. Temporal and spatial reverse transcription PCR (RT-PCR) analysis showed that several genes in the mapped region are associated with the Spc mutant. Although the genes responsible for the Spc mutation were not definitively identified, our results further the current understanding of the complex mechanism underlying the multiple morphological defects in Spc mutants.

Advanced technologies for genetically manipulating the silkworm Bombyx mori, a model Lepidopteran insect

Genetic technologies based on transposon-mediated transgenesis along with several recently developed genome-editing technologies have become the preferred methods of choice for genetically manipulating many organisms. The silkworm, Bombyx mori, is a Lepidopteran insect of great economic importance because of its use in silk production and because it is a valuable model insect that has greatly enhanced our understanding of the biology of insects, including many agricultural pests. In the past 10 years, great advances have been achieved in the development of genetic technologies in B. mori, including transposon-based technologies that rely on piggyBac-mediated transgenesis and genome-editing technologies that rely on protein- or RNA-guided modification of chromosomes. The successful development and application of these technologies has not only facilitated a better understanding of B. mori and its use as a silk production system, but also provided valuable experiences that have contributed to the development of similar technologies in non-model insects. This review summarizes the technologies currently available for use in B. mori, their application to the study of gene function and their use in genetically modifying B. mori for biotechnology applications. The challenges, solutions and future prospects associated with the development and application of genetic technologies in B. mori are also discussed.

Insect Resistance to Bacillus thuringiensis Toxin Cry2Ab Is Conferred by Mutations in an ABC Transporter Subfamily A Protein

The use of conventional chemical insecticides and bacterial toxins to control lepidopteran pests of global agriculture has imposed significant selection pressure leading to the rapid evolution of insecticide resistance. Transgenic crops (e.g., cotton) expressing the Bt Cry toxins are now used world wide to control these pests, including the highly polyphagous and invasive cotton bollworm Helicoverpa armigera. Since 2004, the Cry2Ab toxin has become widely used for controlling H. armigera, often used in combination with Cry1Ac to delay resistance evolution. Isolation of H. armigera and H. punctigera individuals heterozygous for Cry2Ab resistance in 2002 and 2004, respectively, allowed aspects of Cry2Ab resistance (level, fitness costs, genetic dominance, complementation tests) to be characterised in both species. However, the gene identity and genetic changes conferring this resistance were unknown, as was the detailed Cry2Ab mode of action. No cross-resistance to Cry1Ac was observed in mutant lines. Biphasic linkage analysis of a Cry2Ab-resistant H. armigera family followed by exon-primed intron-crossing (EPIC) marker mapping and candidate gene sequencing identified three independent resistance-associated INDEL mutations in an ATP-Binding Cassette (ABC) transporter gene we named HaABCA2. A deletion mutation was also identified in the H. punctigera homolog from the resistant line. All mutations truncate the ABCA2 protein. Isolation of further Cry2Ab resistance alleles in the same gene from field H. armigera populations indicates unequal resistance allele frequencies and the potential for Bt resistance evolution. Identification of the gene involved in resistance as an ABC transporter of the A subfamily adds to the body of evidence on the crucial role this gene family plays in the mode of action of the Bt Cry toxins. The structural differences between the ABCA2, and that of the C subfamily required for Cry1Ac toxicity, indicate differences in the detailed mode-of-action of the two Bt Cry toxins.

Transgenic crops expressing the insecticidal protein Cry2Ab from Bacillus thuringiensis (Bt) are used worldwide to suppress damage by lepidopteran pests, often used in combination with Cry1Ac toxin to delay resistance evolution. Until now, the Cry2Ab mode of action and the mechanism of resistance were unknown, with field-isolated Cry2Ab resistant Helicoverpa armigera showing no cross-resistance to Cry1Ac. In this study, biphasic linkage analysis of a Cry2Ab-resistant H. armigera family followed by EPIC marker mapping and candidate gene sequencing identified three independent INDEL mutations in an ATP-Binding Cassette transporter subfamily A gene (ABCA2). A deletion mutation was identified in the same gene of resistant H. punctigera. All four mutations are predicted to truncate the ABCA2 protein. This is the first molecular genetic characterization of insect resistance to the Cry2Ab toxin, and detection of diverse Cry2Ab resistance alleles will contribute to understanding the micro-evolutionary processes that underpinned lepidopteran Bt-resistance.

InsectBase: a resource for insect genomes and transcriptomes

The genomes and transcriptomes of hundreds of insects have been sequenced. However, insect community lacks an integrated, up-to-date collection of insect gene data. Here, we introduce the first release of InsectBase, available online at http://www.insect-genome.com. The database encompasses 138 insect genomes, 116 insect transcriptomes, 61 insect gene sets, 36 gene families of 60 insects, 7544 miRNAs of 69 insects, 96 925 piRNAs of Drosophila melanogaster and Chilo suppressalis, 2439 lncRNA of Nilaparvata lugens, 22 536 pathways of 78 insects, 678 881 untranslated regions (UTR) of 84 insects and 160 905 coding sequences (CDS) of 70 insects. This release contains over 12 million sequences and provides search functionality, a BLAST server, GBrowse, insect pathway construction, a Facebook-like network for the insect community (iFacebook), and phylogenetic analysis of selected genes.

A comprehensive analysis of the chorion locus in silkmoth

Despite more than 40 years of intense study, essential features of the silkmoth chorion (eggshell) are still not fully understood. To determine the precise structure of the chorion locus, we performed extensive EST analysis, constructed a bacterial artificial chromosome (BAC) contig, and obtained a continuous genomic sequence of 871,711 base pairs. We annotated 127 chorion genes in two segments interrupted by a 164 kb region with 5 non-chorion genes, orthologs of which were on chorion bearing scaffolds in 4 ditrysian families. Detailed transcriptome analysis revealed expression throughout choriogenesis of most chorion genes originally categorized as “middle”, and evidence for diverse regulatory mechanisms including cis-elements, alternative splicing and promoter utilization, and antisense RNA. Phylogenetic analysis revealed multigene family associations and faster evolution of early chorion genes and transcriptionally active pseudogenes. Proteomics analysis identified 99 chorion proteins in the eggshell and micropyle localization of 1 early and 6 Hc chorion proteins.

Identification of two juvenile hormone inducible transcription factors from the silkworm, Bombyx mori

Juvenile hormone (JH) regulates many physiological processes in insects. However, the signal cascades in which JH is active have not yet been fully elucidated, particularly in comparison to another major hormone ecdysteroid. Here we identified two JH inducible transcription factors as candidate components of JH signaling pathways in the silkworm, Bombyx mori. DNA microarray analysis showed that expression of two transcription factor genes, E75 and Enhancer of split mβ (E(spl)mβ), was induced by juvenile hormone I (JH I) in NIAS-Bm-aff3 cells. Real time RT-PCR analysis confirmed that expression of four E75 isoforms (E75A, E75B, E75C and E75D) and E(spl)mβ was 3-8 times greater after JH I addition. Addition of the protein synthesis inhibitor cycloheximide did not suppress JH-induced expression of the genes, indicating that they were directly induced by JH. JH-induced expression of E75 and E(spl)mβ was also observed in four other B. mori cell lines and in larval hemocytes of final instar larvae. Notably, E75A expression was induced very strongly in larval hemocytes by topical application of the JH analog fenoxycarb; the level of induced expression was comparable to that produced by feeding larvae with 20-hydroxyecdysone. These results suggest that E75 and E(spl)mβ are general and direct target genes of JH and that the transcription factors encoded by these genes play important roles in JH signaling.

Sialylation potentials of the silkworm, Bombyx mori; B. mori possesses an active α2,6-sialyltransferase

N-Glycosylation is an important post-translational modification in most secreted and membrane-bound proteins in eukaryotic cells. However, the insect N-glycosylation pathway and the potentials contributing to the N-glycan synthesis are still unclear because most of the studies on these subjects have focused on mammals and plants. Here, we identified Bombyx mori sialyltransferase (BmST), which is a Golgi-localized glycosyltransferase and which can modify N-glycans. BmST was ubiquitously expressed in different organs and in various stages of development and localized at the Golgi. Biochemical analysis using Sf9-expressed BmST revealed that BmST encoded α2,6-sialyltransferase and transferred N-acetylneuraminic acid (NeuAc) to the nonreducing terminus of Galβ1-R, but exhibited the highest activity toward GalNAcβ1,4-GlcNAc-R. Unlike human α2,6-sialyltransferase, BmST required the post-translational modification, especially N-glycosylation, for its full activity. N-Glycoprotein analysis of B. mori fifth instar larvae revealed that high-mannose-type structure was predominant and GlcNAc-linked and fucosylated structures were observed but endogenous galactosyl-, N-acetylgalactosaminyl- and sialyl-N-glycoproteins were undetectable under the standard analytical approach. These results indicate that B. mori genome encodes an α2,6-sialyltransferase, but further investigations of the sialylation potentials are necessary.

A FISH-based chromosome map for the European corn borer yields insights into ancient chromosomal fusions in the silkworm

A significant feature of the genomes of Lepidoptera, butterflies and moths, is the high conservation of chromosome organization. Recent remarkable progress in genome sequencing of Lepidoptera has revealed that syntenic gene order is extensively conserved across phylogenetically distant species. The ancestral karyotype of Lepidoptera is thought to be n=31; however, that of the most well-studied moth, Bombyx mori, is n=28, and diverse studies suggest that three chromosomal fusion events occurred in this lineage. To identify the boundaries between predicted ancient fusions involving B. mori chromosomes 11, 23 and 24, we constructed fluorescence in situ hybridization (FISH)-based chromosome maps of the European corn borer, Ostrinia nubilalis (n=31). We first determined a 511 Mb genomic sequence of the Asian corn borer, O. furnacalis, a congener of O. nubilalis, and isolated bacterial artificial chromosomes and fosmid clones that were expected to localize in candidate regions for the boundaries using these sequences. Combined with FISH and genetic analysis, we narrowed down the candidate regions to 40 kb-1.5 Mb, in strong agreement with a previous estimate based on the genome of a butterfly, Melitaea cinxia. The significant difference in the lengths of the candidate regions where no functional genes were observed may reflect the evolutionary time after fusion events.

Analyses of interactions among pair-rule genes and the gap gene Krüppel in Bombyx segmentation

In the short-germ insect Tribolium, a pair-rule gene circuit consisting of the Tribolium homologs of even-skipped, runt, and odd-skipped (Tc-eve, Tc-run and Tc-odd, respectively) has been implicated in segment formation. To examine the application of the model to other taxa, I studied the expression and function of pair-rule genes in Bombyx mori, together with a Bombyx homolog of Krüppel (Bm-Kr), a known gap gene. Knockdown embryos of Bombyx homologs of eve, run and odd (Bm-eve, Bm-run and Bm-odd) exhibited asegmental phenotypes similar to those of Tribolium knockdowns. However, pair-rule gene interactions were similar to those of both Tribolium and Drosophila, which, different from Tribolium, shows a hierarchical segmentation mode. Additionally, the Bm-odd expression pattern shares characteristics with those of Drosophila pair-rule genes that receive upstream regulatory input. On the other hand, Bm-Kr knockdowns exhibited a large posterior segment deletion as observed in short-germ insects. However, a detailed analysis of these embryos indicated that Bm-Kr modulates expression of pair-rule genes like in Drosophila, although the mechanisms appear to be different. This suggested hierarchical interactions between Bm-Kr and pair-rule genes. Based on these results, I concluded that the pair-rule gene circuit model that describes Tribolium development is not applicable to Bombyx.

Transcriptome analysis of Bombyx mori larval midgut during persistent and pathogenic cytoplasmic polyhedrosis virus infection

Many insects can be persistently infected with viruses but do not show any obvious adverse effects with respect to physiology, development or reproduction. Here, Bombyx mori strain Daizo, persistently infected with cytoplasmic polyhedrosis virus (BmCPV), was used to study the host's transcriptional response after pathogenic infection with the same virus in midgut tissue of larvae persistently and pathogenically infected as 2nd and 4th instars. Next generation sequencing revealed that from 13,769 expressed genes, 167 were upregulated and 141 downregulated in both larval instars following pathogenic infection. Several genes that could possibly be involved in B. mori immune response against BmCPV or that may be induced by the virus in order to increase infectivity were identified, whereas classification of differentially expressed transcripts (confirmed by qRT-PCR) resulted in gene categories related to physical barriers, immune responses, proteolytic/metabolic enzymes, heat-shock proteins, hormonal signaling and uncharacterized proteins. Comparison of our data with the available literature (pathogenic infection of persistently vs. non-persistently infected larvae) unveiled various similarities of response in both cases, which suggests that pre-existing persistent infection does not affect in a major way the transcriptome response against pathogenic infection. To investigate the possible host's RNAi response against BmCPV challenge, the differential expression of RNAi-related genes and the accumulation of viral small RNAs (vsRNAs) were studied. During pathogenic infection, siRNA-like traces like the 2-fold up-regulation of the core RNAi genes Ago-2 and Dcr-2 as well as a peak of 20 nt small RNAs were observed. Interestingly, vsRNAs of the same size were detected at lower rates in persistently infected larvae. Collectively, our data provide an initial assessment of the relative significance of persistent infection of silkworm larvae on the host response following pathogenic infection with CPV, while they also highlight the relative importance of RNAi as an antiviral mechanism.