Bombyx Mori Genome Research Articles

Bombyx mori UFL1 facilitates BmNPV proliferation by regulating host cell apoptosis through PERK.

Ubiquitin-fold modifier 1 (UFM1) is attached to protein substrates through the sequential activity of an E1 (UBA5)-E2 (UFC1)-E3 (UFL1) cascade. UFL1 is the E3 ligase for UFMylation in vertebrates. However, there have been no studies on UFL1 in silkworm to date. In this study, we identified a UFL1 ortholog in Bombyx mori genome. Spatio-temporal expression profiles showed that BmUFL1 expression was high in the midgut, epidermis, and testis and in the pupa-adult stage. BmUFL1 knockdown inhibited B. mori nucleopolyhedrovirus (BmNPV)proliferation, while BmUFL1 overexpression promoted BmNPV proliferation. Mechanically, protein kinase RNA-like endoplasmic reticulum kinase (PERK) signaling and cell apoptosis are involved in BmUFL1-regulatedBmNPV proliferation. Overall, these results suggest that BmUFL1 facilitates BmNPV proliferation in silkworm.

Bmpali, Bmb1 and Bmcap are necessary for uric acid granule formation in Bombyx mori

Uric acid is the end-product of nitrogen metabolism of the silkworm and other lepidopterans. The accumulation of uric acid particles in the epidermis causes the larval silkworm to appear white and opaque. However, the mechanism of uric acid granule formation is still unclear. Silkworm epidermis color is linked to the genes responsible for uric acid particle formation. We first identified two genes in the Bombyx mori genome that encode subunits of the Bloc-1 (Biogenesis of Lysosome-related Organelles Complex-1) by homology to these genes in other eukaryotes, Bmpali and Bmb1. Mutation in these genes caused a transparent phenotype in the silkworm larvae, and the loss of BmBloc-1 subunit gene Bmcap resulted in the same phenotype. These three genes are highly conserved between human and silkworm. We discovered that Bmpali, Bmcap, and Bmb1 localize in the cytoplasm of BmN cells. Yeast two-hybrid assays demonstrated that the Bmpali physically interacts with both Bmcap and Bmb1. Investigating the roles of Bmpali, Bmb1, and Bmcap is essential for uric acid granule formation understanding in Bombyx mori. These mutants present a valuable silkworm model for studying the biogenesis of lysosome-related organelles (LROs).

Glucose-regulated protein 94 facilitates the proliferation of the Bombyx mori nucleopolyhedrovirus via inhibiting apoptosis

Glucose regulatory protein 94 (GRP94) is an endoplasmic reticulum (ER)-resident member of the heat shock protein 90 (HSP90) family, that plays an important role in secreted protein folding. Bombyx mori nuclear polyhedrosis virus (BmNPV) is one of the main pathogens in sericulture, causing serious economic losses every year. Previous studies showed that HSP90 members promote BmNPV replication in silkworm, but the function of BmGRP94 in BmNPV infection and proliferation is still not understood. In this study, we investigated the interplay between BmGRP94 and BmNPV infection in silkworm. We first identified a single gene of BmGRP94 in the Bombyx mori genome, which encodes a polypeptide with 810 amino acids in length. Spatio-temporal expression profiles showed that BmGRP94 was highly expressed in hemocytes and midgut, and was significantly induced by BmNPV infection. Furthermore, overexpression of BmGRP94 facilitates viral proliferation, while BmGRP94 inhibition evidently decreased BmNPV proliferation in BmN cells and in silkworm midgut. Mechanistically, BmGRP94 inhibition triggers ER stress, as judged by increased expression of PERK/ATF4/ERO1, H2O2 production, and ER calcium efflux, which promotes cell apoptosis to restrict BmNPV replication in silkworm. These results suggest that BmGRP94 plays an important role in facilitating BmNPV proliferation, and provides a potential molecular target for BmNPV prevention.

Overexpression of BmJHBPd2 Repressed Silk Synthesis by Inhibiting the JH/Kr-h1 Signaling Pathway in Bombyx mori.

The efficient production of silkworm silk is crucial to the silk industry. Silk protein synthesis is regulated by the juvenile hormone (JH) and 20-Hydroxyecdysone (20E). Therefore, the genetic regulation of silk production is a priority. JH binding protein (JHBP) transports JH from the hemolymph to target organs and cells and protects it. In a previous study, we identified 41 genes containing a JHBP domain in the Bombyx mori genome. Only one JHBP gene, BmJHBPd2, is highly expressed in the posterior silk gland (PSG), and its function remains unknown. In the present study, we investigated the expression levels of BmJHBPd2 and the major silk protein genes in the high-silk-producing practical strain 872 (S872) and the low-silk-producing local strain Dazao. We found that BmJHBPd2 was more highly expressed in S872 than in the Dazao strain, which is consistent with the expression pattern of fibroin genes. A subcellular localization assay indicated that BmJHBPd2 is located in the cytoplasm. In vitro hormone induction experiments showed that BmJHBPd2 was upregulated by juvenile hormone analogue (JHA) treatment. BmKr-h1 upregulation was significantly inhibited by the overexpression of BmJHBPd2 (BmJHBPd2OE) at the cell level when induced by JHA. However, overexpression of BmJHBPd2 in the PSG by transgenic methods led to the inhibition of silk fibroin gene expression, resulting in a reduction in silk yield. Further investigation showed that in the transgenic BmJHBPd2OE silkworm, the key transcription factor of the JH signaling pathway, Krüppel homolog 1 (Kr-h1), was inhibited, and 20E signaling pathway genes, such as broad complex (Brc), E74A, and ultraspiracle protein (USP), were upregulated. Our results indicate that BmJHBPd2 plays an important role in the JH signaling pathway and is important for silk protein synthesis. Furthermore, our findings help to elucidate the mechanisms by which JH regulates silk protein synthesis.

Comparison of Long-Read Methods for Sequencing and Assembly of Lepidopteran Pest Genomes

Lepidopteran species are mostly pests, causing serious annual economic losses. High-quality genome sequencing and assembly uncover the genetic foundation of pest occurrence and provide guidance for pest control measures. Long-read sequencing technology and assembly algorithm advances have improved the ability to timeously produce high-quality genomes. Lepidoptera includes a wide variety of insects with high genetic diversity and heterozygosity. Therefore, the selection of an appropriate sequencing and assembly strategy to obtain high-quality genomic information is urgently needed. This research used silkworm as a model to test genome sequencing and assembly through high-coverage datasets by de novo assemblies. We report the first nearly complete telomere-to-telomere reference genome of silkworm Bombyx mori (P50T strain) produced by Pacific Biosciences (PacBio) HiFi sequencing, and highly contiguous and complete genome assemblies of two other silkworm strains by Oxford Nanopore Technologies (ONT) or PacBio continuous long-reads (CLR) that were unrepresented in the public database. Assembly quality was evaluated by use of BUSCO, Inspector, and EagleC. It is necessary to choose an appropriate assembler for draft genome construction, especially for low-depth datasets. For PacBio CLR and ONT sequencing, NextDenovo is superior. For PacBio HiFi sequencing, hifiasm is better. Quality assessment is essential for genome assembly and can provide better and more accurate results. For chromosome-level high-quality genome construction, we recommend using 3D-DNA with EagleC evaluation. Our study references how to obtain and evaluate high-quality genome assemblies, and is a resource for biological control, comparative genomics, and evolutionary studies of Lepidopteran pests and related species.

Genetic Code Expansion System for Tight Control of Gene Expression in Bombyx mori Cell Lines.

Simple SummaryBombyx mori is a lepidopteran insect with economic value. Its genetic background is clear, and genome sequence is relatively complete, but the function of many genes has not been determined. The genetic code expansion system has become an important means of gene function research. In this study, a genetic code expansion system suitable for B. mori cells was established. This system included a modified tRNAPyl/Pyrrolysyl-tRNA synthetase (PylRS) pair from Methanosarcina mazei, the reporter gene D[TAG]G formed by DsRed and EGFP through amber stop codon TAG connection and the unnatural amino acid H-Lys(Boc)-OH. In silkworm BmE and BmNs cell lines, the reporter gene expression was strictly controlled by H-Lys(Boc)-OH in the presence of both PylRS and tRNAPyl. The silkworm genetic code expansion system established here is another useful controllable gene expression system besides tetracycline induced expression system.Inducible gene expression systems are important tools for studying gene function and to control protein synthesis. With the completion of the detailed map of the silkworm (Bombyx mori) genome, the study of Bombyx mori has entered the post-genome era. While the functions of many genes have been described in detail, many coding genes remain unidentified. Except for the available tetracycline induction system, there is currently a dearth of other effective induction systems for B. mori. A genetic code expansion system can be used for protein labeling and to regulate gene expression. Here, we have established a genetic code expansion system for B. mori based on the well-researched tRNAPyl/PylRS pair from Methanosarcina mazei. We used H-Lys(Boc)-OH, which is a lysine derivative to efficiently and tightly control the expression of the reporter gene DsRed[TAG]EGFP (D[TAG]G), which encoded a H-Lys(Boc)-OH-bearing protein fused with DsRed and EGFP (here regarded as D[Boc]G) in B. mori cell lines BmE and BmNs. In D[TAG]G, the amber stop codon is recognized as the orthogonal tRNAPyl. Successful application of genetic code expansion system in silkworm cell lines will support the research into the function of silkworm genes and paves the way for the identification of new genes and protein markers in silkworm.

The Landscapes of Full-Length Transcripts and Splice Isoforms as Well as Transposons Exonization in the Lepidopteran Model System, Bombyx mori.

The domesticated silkworm, Bombyx mori, is an important model system for the order Lepidoptera. Currently, based on third-generation sequencing, the chromosome-level genome of Bombyx mori has been released. However, its transcripts were mainly assembled by using short reads of second-generation sequencing and expressed sequence tags which cannot explain the transcript profile accurately. Here, we used PacBio Iso-Seq technology to investigate the transcripts from 45 developmental stages of Bombyx mori. We obtained 25,970 non-redundant high-quality consensus isoforms capturing ∼60% of previous reported RNAs, 15,431 (∼47%) novel transcripts, and identified 7,253 long non-coding RNA (lncRNA) with a large proportion of novel lncRNA (∼56%). In addition, we found that transposable elements (TEs) exonization account for 11,671 (∼45%) transcripts including 5,980 protein-coding transcripts (∼32%) and 5,691 lncRNAs (∼79%). Overall, our results expand the silkworm transcripts and have general implications to understand the interaction between TEs and their host genes. These transcripts resource will promote functional studies of genes and lncRNAs as well as TEs in the silkworm.

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.

CRISPR/Cas9 Initiated Transgenic Silkworms as a Natural Spinner of Spider Silk.

Using transgenic silkworms with their natural spinning apparatus has proven to be a promising way to spin spider silk-like fibers. The challenges are incorporating native-size spider silk proteins and achieving an inheritable transgenic silkworm strain. In this study, a CRISPR/Cas9 initiated fixed-point strategy was used to successfully incorporate spider silk protein genes into the Bombyx mori genome. Native-size spider silk genes (up to 10 kb) were inserted into an intron of the fibroin heavy or light chain (FibH or FibL) ensuring that any sequence changes induced by the CRISPR/Cas9 would not impact protein production. The resulting fibers are as strong as native spider silks (1.2 GPa tensile strength). The transgenic silkworms have been tracked for several generations with normal inheritance of the transgenes. This strategy demonstrates the feasibility of using silkworms as a natural spider silk spinner for industrial production of high-performance fibers.

Insights into the Structure of the Spruce Budworm (Choristoneura fumiferana) Genome, as Revealed by Molecular Cytogenetic Analyses and a High-Density Linkage Map.

Genome structure characterization can contribute to a better understanding of processes such as adaptation, speciation, and karyotype evolution, and can provide useful information for refining genome assemblies. We studied the genome of an important North American boreal forest pest, the spruce budworm, Choristoneura fumiferana, through a combination of molecular cytogenetic analyses and construction of a high-density linkage map based on single nucleotide polymorphism (SNP) markers obtained through a genotyping-by-sequencing (GBS) approach. Cytogenetic analyses using fluorescence in situ hybridization methods confirmed the haploid chromosome number of n = 30 in both sexes of C. fumiferana and showed, for the first time, that this species has a WZ/ZZ sex chromosome system. Synteny analysis based on a comparison of the Bombyx mori genome and the C. fumiferana linkage map revealed the presence of a neo-Z chromosome in the latter species, as previously reported for other tortricid moths. In this neo-Z chromosome, we detected an ABC transporter C2 (ABCC2) gene that has been associated with insecticide resistance. Sex-linkage of the ABCC2 gene provides a genomic context favorable to selection and rapid spread of resistance against Bacillus thuringiensis serotype kurstaki (Btk), the main insecticide used in Canada to control spruce budworm populations. Ultimately, the linkage map we developed, which comprises 3586 SNP markers distributed over 30 linkage groups for a total length of 1720.41 cM, will be a valuable tool for refining our draft assembly of the spruce budworm genome.

An integrated CRISPR Bombyx mori genome editing system with improved efficiency and expanded target sites

Genome editing enabled unprecedented new opportunities for targeted genomic engineering of a wide variety of organisms ranging from microbes, plants, animals and even human embryos. The serial establishing and rapid applications of genome editing tools significantly accelerated Bombyx mori (B. mori) research during the past years. However, the only CRISPR system in B. mori was the commonly used SpCas9, which only recognize target sites containing NGG PAM sequence. In the present study, we first improve the efficiency of our previous established SpCas9 system by 3.5 folds. The improved high efficiency was also observed at several loci in both BmNs cells and B. mori embryos. Then to expand the target sites, we showed that two newly discovered CRISPR system, SaCas9 and AsCpf1, could also induce highly efficient site-specific genome editing in BmNs cells, and constructed an integrated CRISPR system. Genome-wide analysis of targetable sites was further conducted and showed that the integrated system cover 69,144,399 sites in B. mori genome, and one site could be found in every 6.5 bp. The efficiency and resolution of this CRISPR platform will probably accelerate both fundamental researches and applicable studies in B. mori, and perhaps other insects.

Genome-Wide Identification, Characterization and Expression Analysis of the Solute Carrier 6 Gene Family in Silkworm (Bombyx mori).

The solute carrier 6 (SLC6) gene family, initially known as the neurotransmitter transporters, plays vital roles in the regulation of neurotransmitter signaling, nutrient absorption and motor behavior. In this study, a total of 16 candidate genes were identified as SLC6 family gene homologs in the silkworm (Bombyx mori) genome. Spatio-temporal expression patterns of silkworm SLC6 gene transcripts indicated that these genes were highly and specifically expressed in midgut, brain and gonads; moreover, these genes were expressed primarily at the feeding stage or adult stage. Levels of expression for most midgut-specific and midgut-enriched gene transcripts were down-regulated after starvation but up-regulated after re-feeding. In addition, we observed that expression levels of these genes except for BmSLC6-15 and BmGT1 were markedly up-regulated by a juvenile hormone analog. Moreover, brain-enriched genes showed differential expression patterns during wandering and mating processes, suggesting that these genes may be involved in modulating wandering and mating behaviors. Our results improve our understanding of the expression patterns and potential physiological functions of the SLC6 gene family, and provide valuable information for the comprehensive functional analysis of the SLC6 gene family.

Development of a Genomic Resource and Quantitative Trait Loci Mapping of Male Calling Traits in the Lesser Wax Moth, Achroia grisella.

In the study of sexual selection among insects, the Lesser Waxmoth, Achroia grisella (Lepidoptera: Pyralidae), has been one of the more intensively studied species over the past 20 years. Studies have focused on how the male calling song functions in pair formation and on the quantitative genetics of male song characters and female preference for the song. Recent QTL studies have attempted to elucidate the genetic architecture of male song and female preference traits using AFLP markers. We continued these QTL studies using SNP markers derived from an EST library that allowed us to measure both DNA sequence variation and map loci with respect to the lepidopteran genome. We report that the level of sequence variation within A. grisella is typical among other Lepidoptera that have been examined, and that comparison with the Bombyx mori genome shows that macrosynteny is conserved. Our QTL map shows that a QTL for a male song trait, pulse-pair rate, is situated on the Z chromosome, a prediction for sexually selected traits in Lepidoptera. Our findings will be useful for future studies of genetic architecture of this model species and may help identify the genetics associated with the evolution of its novel acoustic communication.

Identification and molecular structure analysis of a new noncoding RNA, a sbRNA homolog, in the silkworm Bombyx mori genome.

The small noncoding group of RNAs called stem-bulge RNAs (sbRNAs), first reported in Caenorhabditis elegans, is described as molecules homologous to the Y RNAs, a specific class of noncoding RNAs that is present in vertebrates. This homology indicates the possibility of the existence of sbRNAs in other invertebrate organisms. In this work, we used bioinformatic tools and conserved sequences of sbRNAs from C. Elegans and Y RNAs to search for homologous sbRNA sequences in the Bombyx mori genome. This analysis led to the discovery of one noncoding gene, which was translated into RNA segments and comparatively analysed with segments from human and hamster Y RNAs and C. elegans sbRNAs in molecular dynamic simulations. This gene represents the first evidence for a new sbRNA-like noncoding RNA, the BmsbRNA gene, in this Lepidoptera genome.

In vivo site-specific integration of transgene in silkworm via PhiC31 integrase-mediated cassette exchange

Current techniques for genetic engineering of the silkworm Bombyx mori genome utilize transposable elements, which result in positional effects and insertional mutagenesis through random insertion of exogenous DNA. New methods for introducing transgenes at specific positions are therefore needed to overcome the limitations of transposon-based strategies. Although site-specific recombination systems have proven powerful tools for genome manipulation in many organisms, their use has not yet been well established for the integration of transgenes in the silkworm. We describe a method for integrating target genes at pre-defined chromosomal sites in the silkworm via phiC31/att site-specific recombination system-mediated cassette exchange. Successful recombinase-mediated cassette exchange (RMCE) was observed in the two transgenic target strains with an estimated transformation efficiency of 3.84–7.01%. Our results suggest that RMCE events between chromosomal attP/attP target sites and incoming attB/attB sites were more frequent than those in the reciprocal direction. This is the first report of in vivo RMCE via phiC31 integrase in the silkworm, and thus represents a key step toward establishing genome manipulation technologies in silkworms and other lepidopteran species.

Characterization of a novel chromodomain-containing gene from the silkworm, Bombyx mori

Heterochromatin protein 1 (HP1) is an evolutionarily conserved protein across different eukaryotic species, and is crucial in the establishment and maintenance of heterochromatin. HP1 proteins have two distinct functional domains, an N-terminal chromodomain (CD) and a C-terminal chromoshadow domain (CSD), which are required for the selective binding of HP1 proteins to modified histones. During our screen for HP1-like proteins in the Bombyx mori genome, we found a novel silkworm gene, Bombyx mori chromodomain protein 1 (BmCdp1), encoding a putative chromobox protein with only two CDs. The BmCdp1 family proteins are closely related to the HP1 proteins, and most of them belong to insect lineages. qRT-PCR analysis indicated that BmCdp1 mRNA was most abundantly expressed in early embryos, and relatively higher expression was observed in larval testes, hemocytes, and pupal ovaries. Western blot and immunostaining experiments showed that BmCdp1 was localized mainly in the nucleus of BmN4 cells. We searched BmCdp1-bound loci in the Bombyx genome by ChIP-seq analysis using Flag-tagged BmCdp1-expressing BmN4 cells. Combined with ChIP-qPCR experiments, we identified two reliable BmCdp1-bound loci in the genome. siRNA-mediated knockdown of BmCdp1 in BmN4 cells and early embryos did not affect the expression of the gene located close to the BmCdp1-bound locus.

FISH identification of Helicoverpa armigera and Mamestra brassicae chromosomes by BAC and fosmid probes

Since the Bombyx mori genome sequence was published, conserved synteny between B. mori and some other lepidopteran species has been revealed by either FISH (fluorescence in situ hybridization) with BAC (bacterial artificial chromosome) probes or linkage analysis. However, no species belonging to the Noctuidae, the largest lepidopteran family which includes serious polyphagous pests, has been analyzed so far with respect to genome-wide conserved synteny and gene order. For that purpose, we selected the noctuid species Helicoverpa armigera and Mamestra brassicae, both with n = 31 chromosomes. Gene-defined fosmid clones from M. brassicae and BAC clones from a closely related species of H. armigera, Heliothis virescens, were used for a FISH analysis on pachytene chromosomes. We recognized all H. armigera chromosomes from specific cross-hybridization signals of 146 BAC probes. With 100 fosmid clones we identified and characterized all 31 bivalents of M. brassicae. Synteny and gene order were well conserved between the two noctuid species. The comparison with the model species B. mori (n = 28) showed the same phenomenon for 25 of the 28 chromosomes. Three chromosomes (#11, #23 and #24) had two counterparts each in H. armigera and M. brassicae. Since n = 31 is the modal chromosome number in Lepidoptera, the noctuid chromosomes probably represent an ancestral genome organization of Lepidoptera. This is the first identification of a full karyotype in Lepidoptera by means of BAC cross-hybridization between species. The technique shows the potential to expand the range of analyzed species efficiently.

Paralytic Peptide Binding Protein (PP-BP) Gene Expression During Egg Diapause and Its Multi-Gene Organization in the Silkworm Bombyx mori

Paralytic peptide binding proteins (PP-BP) are 30KP proteins that show similarity to ENF binding proteins. The ENF-BP act as active regulators of ENF peptides. ENF peptides are multifunctional insect cytokines. The comparison of gene expression in diapause induced and non-diapause eggs at different time intervals after oviposition showed an upregulation of PP at 18h as well as PP-BP at 12 and 18h after oviposition along with few other genes. The current study has been taken up to investigate the role of PP as well as PP-BP in diapause induction in polyvoltine silkworms and to study the multigene organization of PP-BP in the Bombyx mori genome. The tissue specific expression analysis revealed that, PP-BP is highly expressed in fat body followed by egg and brain while no expression was observed in midgut. The expression levels of PP and PP-BP in diapause and non-diapause eggs from 0h to 48h after oviposition, validated through realtime PCR revealed that PP is highly expressed at 18 and 24h while PP-BP expression is higher at 12 and 18h time intervals suggesting their possible role in diapause induction. The whole genome survey of the PP-BP paralogous sequences revealed a total of 46 B. mori PP-BP homologs that are classified into 3 categories viz., ENF-BP, Typical 30KPs and serine/threonine rich 30KPs. These paralogous sequences are distributed on chromosomes 7, 20, 22 and 24, all 30KP and S/T rich 30KP proteins are present in the same locus of chromosome 20.

Characterization of an Ancient Lepidopteran Lateral Gene Transfer

Bacteria to eukaryote lateral gene transfers (LGT) are an important potential source of material for the evolution of novel genetic traits. The explosion in the number of newly sequenced genomes provides opportunities to identify and characterize examples of these lateral gene transfer events, and to assess their role in the evolution of new genes. In this paper, we describe an ancient lepidopteran LGT of a glycosyl hydrolase family 31 gene (GH31) from an Enterococcus bacteria. PCR amplification between the LGT and a flanking insect gene confirmed that the GH31 was integrated into the Bombyx mori genome and was not a result of an assembly error. Database searches in combination with degenerate PCR on a panel of 7 lepidopteran families confirmed that the GH31 LGT event occurred deep within the Order approximately 65–145 million years ago. The most basal species in which the LGT was found is Plutella xylostella (superfamily: Yponomeutoidea). Array data from Bombyx mori shows that GH31 is expressed, and low dN/dS ratios indicates the LGT coding sequence is under strong stabilizing selection. These findings provide further support for the proposition that bacterial LGTs are relatively common in insects and likely to be an underappreciated source of adaptive genetic material.

FLP Recombinase-Mediated Site-Specific Recombination in Silkworm, Bombyx mori

A comprehensive understanding of gene function and the production of site-specific genetically modified mutants are two major goals of genetic engineering in the post-genomic era. Although site-specific recombination systems have been powerful tools for genome manipulation of many organisms, they have not yet been established for use in the manipulation of the silkworm Bombyx mori genome. In this study, we achieved site-specific excision of a target gene at predefined chromosomal sites in the silkworm using a FLP/FRT site-specific recombination system. We first constructed two stable transgenic target silkworm strains that both contain a single copy of the transgene construct comprising a target gene expression cassette flanked by FRT sites. Using pre-blastoderm microinjection of a FLP recombinase helper expression vector, 32 G3 site-specific recombinant transgenic individuals were isolated from five of 143 broods. The average frequency of FLP recombinase-mediated site-specific excision in the two target strains genome was approximately 3.5%. This study shows that it is feasible to achieve site-specific recombination in silkworms using the FLP/FRT system. We conclude that the FLP/FRT system is a useful tool for genome manipulation in the silkworm. Furthermore, this is the first reported use of the FLP/FRT system for the genetic manipulation of a lepidopteran genome and thus provides a useful reference for the establishment of genome manipulation technologies in other lepidopteran species.