Codon Optimization-based Whole-gene Scanning Identifies Hidden Nucleotides Essential for Bombyx mori Nucleopolyhedrovirus polyhedrin Hyperexpression

During the late stage of infection, alphabaculoviruses produce many occlusion bodies (OBs) in the nuclei of the insect host's cells through the hyperexpression of polyhedrin (POLH), a major OB component encoded by polh. The strong polh promoter has been used to develop a baculovirus expression vector system for recombinant protein expression in cultured insect cells and larvae. However, the relationship between POLH accumulation and the polh coding sequence remains largely unelucidated. This study aimed to assess the importance of polh codon usage and/or nucleotide sequences in POLH accumulation by generating a baculovirus Bombyx mori nucleopolyhedrovirus (BmNPV) expressing mutant polh (co-polh) optimized according to the codon preference of its host insect. Although the deduced amino acid sequence of CO-POLH was the same as that of wild-type POLH, POLH accumulation was significantly lower in cells infected with the co-polh mutant. This reduction was due to decreased polh mRNA levels rather than translational repression. Analysis of mutant viruses with chimeric polh revealed that a 30 base-pair (bp) 5' proximal polh coding region was necessary for maintaining high polh mRNA levels. Sequence comparison of wild-type polh and co-polh identified five nucleotide differences in this region, indicating that these nucleotides were critical for polh hyperexpression. Furthermore, luciferase reporter assays showed that the 30 bp 5' coding region was sufficient for maintaining the polh promoter-driven high level of polh mRNA. Thus, our whole-gene scanning by codon optimization identified important hidden nucleotides for polh hyperexpression in alphabaculoviruses.

Expression of the Wolbachia male-killing factor Oscar impairs dosage compensation in lepidopteran embryos

Wolbachia are intracellular bacteria in insects that can manipulate the sexual development and reproduction by male killing or other methods. We have recently identified a Wolbachia protein named Oscar that acts as a male-killing factor for lepidopteran insects. Oscar interacts with the Masculinizer (Masc) protein, which is required for both masculinization and dosage compensation (DC) in lepidopteran insects. Embryonic expression of Oscar inhibits masculinization and causes male killing in two lepidopteran species, Ostrinia furnacalis and Bombyx mori. However, it remains unknown whether Oscar-induced male killing is caused by a failure of DC. Here, we performed a transcriptome analysis of Oscar complementary RNA-injected O. furnacalis and B. mori embryos, and found that Oscar primarily targets the Masc protein, resulting in male killing by interfering with DC in lepidopteran insects.

A seamless connection from the burst sequence to the start codon is essential for polyhedrin hyperexpression in alphabaculoviruses

Alphabaculoviruses produce a large number of occlusion bodies (OBs) in host cells during the late stage of infection. OBs are mainly composed of polyhedrin (POLH), and high-level transcription of the polh gene has been exploited to express foreign proteins in insect cells. While making Bombyx mori nucleopolyhedrovirus (BmNPV) polh mutants using a conventional transfer vector-based method, we noticed that a virus with a short sequence insertion just before the polh start codon produces fewer very small OBs. Detailed analysis of several BmNPV mutants revealed that insertions between the burst sequence and start codon markedly decrease POLH accumulation and polh transcription. We further confirmed this decrease using recombinant viruses expressing a reporter gene driven by the polh promoter. These findings underscore the critical importance of a seamless connection from the burst sequence to the start codon for baculovirus polh hyperexpression.

Two Complete Genomes of Male-Killing Wolbachia Infecting Ostrinia Moth Species Illuminate Their Evolutionary Dynamics and Association with Hosts

Wolbachia is an extremely widespread intracellular symbiont which causes reproductive manipulation on various arthropod hosts. Male progenies are killed in Wolbachia-infected lineages of the Japanese Ostrinia moth population. While the mechanism of male killing and the evolutionary interaction between host and symbiont are significant concerns for this system, the absence of Wolbachia genomic information has limited approaches to these issues. We determined the complete genome sequences of wFur and wSca, the male-killing Wolbachia of Ostrinia furnacalis and Ostrinia scapulalis. The two genomes shared an extremely high degree of homology, with over 95% of the predicted protein sequences being identical. A comparison of these two genomes revealed nearly minimal genome evolution, with a strong emphasis on the frequent genome rearrangements and the rapid evolution of ankyrin repeat-containing proteins. Additionally, we determined the mitochondrial genomes of both species' infected lineages and performed phylogenetic analyses to deduce the evolutionary dynamics of Wolbachia infection in the Ostrinia clade. According to the inferred phylogenetic relationship, two possible scenarios were proposed: (1) Wolbachia infection was established in the Ostrinia clade prior to the speciation of related species such as O. furnacalis and O. scapulalis, or (2) Wolbachia infection in these species was introgressively transferred from a currently unidentified relative. Simultaneously, the relatively high homology of mitochondrial genomes suggested recent Wolbachia introgression between infected Ostrinia species. The findings of this study collectively shed light on the host-symbiont interaction from an evolutionary standpoint.

A method for screening the suppressor genes of siRNA and piRNA pathways using cultured silkworm cells

The BmN-4 cell line originates from the ovaries of silkworm, Bombyx mori , and possesses endogenous small interfering RNA (siRNA) and PIWI-interacting RNA (piRNA) pathways. BmN-4 cells are latently infected with Bombyx mori latent virus (BmLV), an RNA virus whose replication is strictly controlled by both siRNA and piRNA pathways. Knockdown or knockout of the core factors of these two small RNA pathways increases BmLV RNA amount, which in turn inhibits cell growth. Here, we used the known RNAi suppressor CrPV-1A to assess whether the BmN-4 cell line can be used for screening the suppressors of siRNA and piRNA pathways.

Non-gonadal somatic piRNA pathways ensure sexual differentiation, larval growth, and wing development in silkworms

PIWI-interacting RNAs (piRNAs) guide PIWI proteins to target transposons in germline cells, thereby suppressing transposon activity to preserve genome integrity in metazoans' gonadal tissues. Piwi, one of three Drosophila PIWI proteins, is expressed in the nucleus and suppresses transposon activity by forming heterochromatin in an RNA cleavage-independent manner. Recently, Piwi was reported to control cell metabolism in Drosophila fat body, providing an example of piRNAs acting in non-gonadal somatic tissues. However, mutant flies of the other two PIWI proteins, Aubergine (Aub) and Argonaute3 (Ago3), show no apparent phenotype except for infertility, blurring the importance of the piRNA pathway in non-gonadal somatic tissues. The silkworm, Bombyx mori, possesses two PIWI proteins, Siwi (Aub homolog) and BmAgo3 (Ago3 homolog), whereas B. mori does not have a Piwi homolog. Siwi and BmAgo3 are mainly expressed in gonadal tissues and play a role in repressing transposon activity by cleaving transposon RNA in the cytoplasm. Here, we generated Siwi and BmAgo3 loss-of-function mutants of B. mori and found that they both showed delayed larval growth and failed to become adult moths. They also exhibited defects in wing development and sexual differentiation. Transcriptome analysis revealed that loss of somatic piRNA biogenesis pathways results in abnormal expression of not only transposons but also host genes, presumably causing severe growth defects. Our results highlight the roles of non-gonadal somatic piRNAs in B. mori development.

Masculinizer is not post-transcriptionally regulated by female-specific piRNAs during sex determination in the Asian corn borer, Ostrinia furnacalis

Lepidopteran insects are heterogametic in females, although most insect species are heterogametic in males. In a lepidopteran model species, the silkworm Bombyx mori (Bombycoidea), the uppermost sex determinant Feminizer (Fem) has been identified on the female-specific W chromosome. Fem is a precursor of PIWI-interacting small RNA (piRNA). Fem piRNA forms a complex with Siwi, one of the two B. mori PIWI-clade Argonaute proteins. In female embryos, Fem piRNA-Siwi complex cleaves the mRNA of the male-determining gene Masculinizer (Masc), directing the female-determining pathway. In male embryos, Masc activates the male-determining pathway in the absence of Fem piRNA. Recently, W chromosome-derived piRNAs complementary to Masc mRNA have also been identified in the diamondback moth Plutella xylostella (Yponomeutoidea), indicating the convergent evolution of piRNA-dependent sex determination in Lepidoptera. Here, we show that this is not the case in the Asian corn borer, Ostrinia furnacalis (Pyraloidea). Although our previous studies demonstrated that O. furnacalis Masc (OfMasc) has a masculinizing function in the embryonic stage, the expression level of OfMasc was indistinguishable between the sexes at the timing of sex determination. Deep sequencing analysis identified no female-specific small RNAs mapped onto OfMasc mRNA. Embryonic knockdown of two PIWI genes did not affect the expression level of OfMasc in either sex. These results demonstrated that piRNA-dependent reduction of Masc mRNA in female embryos is not a common strategy of sex determination, which suggests the possibility of divergent evolution of sex determinants across the order Lepidoptera.

Recessive embryonic lethal mutations uncovered in heterozygous condition in silkworm semiconsomic strains

In this study, we found two embryonic lethal mutations, t04 lethal (l-t04) and m04 lethal (l-m04), in semiconsomic strains T04 and M04, respectively. In these semiconsomic strains, the entire diploid genome, except for one chromosome 4 of the wild silkworm Bombyx mandarina, is substituted with chromosomes of the domesticated silkworm B. mori, and l-t04 and l-m04 mutations are located on B. mandarina-derived chromosome 4. To clarify the cause of the lethalities and the genes responsible for these mutations, positional cloning and CRISPR/Cas9 mediated knockout screening were performed. Finally, genetic complementation tests identified the mutations responsible for the l-t04 and l-m04 as the Bombyx homolog of imaginal discs arrested (Bmida) and TATA box binding protein-associated factor 5 (BmTaf5), respectively. Lethal stages of each knockout mutant indicated the importance of these genes in B. mori late embryogenesis. The lethal mutations responsible for l-t04 and l-m04 were not found in parental strains or wild B. mandarina collected from 39 distinct locations in Japan, indicating that both mutations were independently introduced during or after the development of the semiconsomic strains. We conclude that the recessive embryonic lethality in the T04 and M04 strains is due to deleterious mutations produced in B. mandarina-derived chromosome 4.

Expression and localization of Bombyx mori nucleopolyhedrovirus GP37

Mitochondria play an essential role in intracellular energy metabolism. This study described the involvement of Bombyx mori nucleopolyhedrovirus (BmNPV) GP37 (BmGP37) in host mitochondria. Herein, the proteins associated with host mitochondria isolated from BmNPV-infected or mock-infected cells by two-dimensional gel electrophoresis were compared. One mitochondria-associated protein in virus-infected cells was identified as BmGP37 by liquid chromatography-mass spectrometry analysis. Furthermore, the BmGP37 antibodies were generated, which could react specifically with BmGP37 in the BmNPV-infected BmN cells. Western blot experiments showed that BmGP37 was expressed at 18 h post-infection and was verified as a mitochondria-associated protein. Immunofluorescence analysis demonstrated that BmGP37 localized to the host mitochondria during BmNPV infection. Furthermore, western blot analysis revealed that BmGP37 is a novel component protein of the occlusion-derived virus (ODV) of BmNPV. The present results indicated that BmGP37 is one of the ODV-associated proteins and may have important roles in host mitochondria during BmNPV infection.

The piRNA cluster torimochi is an expanding transposon in cultured silkworm cells

PIWI proteins and PIWI-interacting RNAs (piRNAs) play a central role in repressing transposable elements in animal germ cells. It is thought that piRNAs are mainly produced from discrete genomic loci named piRNA clusters, which often contain many "dead" transposon remnants from past invasions and have heterochromatic features. In the genome of silkworm ovary-derived cultured cells called BmN4, a well-established model for piRNA research, torimochi was previously annotated as a unique and specialized genomic region that can capture transgenes and produce new piRNAs bearing a trans-silencing activity. However, the sequence identity of torimochi has remained elusive. Here, we carefully characterized torimochi by utilizing the updated silkworm genome sequence and the long-read sequencer MinION. We found that torimochi is in fact a full-length gypsy-like LTR retrotransposon, which is exceptionally active and has massively expanded its copy number in BmN4 cells. Many copies of torimochi in BmN4 cells have features of open chromatin and the ability to produce piRNAs. Therefore, torimochi may represent a young, growing piRNA cluster, which is still "alive" and active in transposition yet capable of trapping other transposable elements to produce de novo piRNAs.

The two Gtsf paralogs in silkworms orthogonally activate their partner PIWI proteins for target cleavage

The PIWI-interacting RNA (piRNA) pathway is a protection mechanism against transposons in animal germ cells. Most PIWI proteins possess piRNA-guided endonuclease activity, which is critical for silencing transposons and producing new piRNAs. Gametocyte-specific factor 1 (Gtsf1), an evolutionarily conserved zinc finger protein, promotes catalysis by PIWI proteins. Many animals have multiple Gtsf1 paralogs; however, their respective roles in the piRNA pathway are not fully understood. Here, we dissected the roles of Gtsf1 and its paralog Gtsf1-like (Gtsf1L) in the silkworm piRNA pathway. We found that Gtsf1 and Gtsf1L preferentially bind the two silkworm PIWI paralogs, Siwi and BmAgo3, respectively, and facilitate the endonuclease activity of each PIWI protein. This orthogonal activation effect was further supported by specific reduction of BmAgo3-bound Masculinizer piRNA and Siwi-bound Feminizer piRNA, the unique piRNA pair required for silkworm feminization, upon depletion of Gtsf1 and Gtsf1L, respectively. Our results indicate that the two Gtsf paralogs in silkworms activate their respective PIWI partners, thereby facilitating the amplification of piRNAs.

Masculinizer-induced dosage compensation is achieved by transcriptional downregulation of both copies of Z-linked genes in the silkworm, Bombyx mori

The evolution of dosage compensation produces similar expression of sex-linked and autosomal genes in the heterogametic sex. The silkworm (Bombyx mori), a lepidopteran insect, has a female heterogametic WZ sex determination system. A Z-linked gene, Masculinizer (Masc), is the primary determinant of maleness and dosage compensation in B. mori. However, it remains unknown whether one of the two Z chromosomes is inactivated or both Z chromosomes are suppressed in B. mori males. Hence, we performed transcriptome analysis using hybrids between two B. mori strains and analysed allele-specific expression to distinguish these alternatives. Our analysis revealed that genes on both the maternal and paternal Z chromosomes are transcriptionally upregulated in Masc knocked down males. We therefore conclude that both Z chromosomes are transcriptionally downregulated in B. mori males, similar to the system in Caenorhabditis elegans.

Mutations in the polyhedrin NLS affect the assembly and polyhedral shape of alphabaculovirus occlusion bodies

Alphabaculoviruses produce occlusion bodies (OBs) in the nucleus of the infected cells at the late stage of infection. OBs are mainly composed of a single viral protein called polyhedrin (POLH). Autographa californica multiple nucleopolyhedrovirus (AcMNPV) POLH possesses a monopartite nuclear localization signal sequence (NLS), KRKK, from 32nd to 35th residues. However, the functions of POLH NLS of other alphabaculoviruses remain unknown. Here, POLH NLS mutants of Bombyx mori nucleopolyhedrovirus (BmNPV) were generated and NLS function as well as the relationship between NLS and OB localization or morphology was investigated. Deletion or mutation of BmNPV POLH NLS severely affected POLH and OB intracellular localization. Additionally, viruses in which the arginine residue at the 33rd position of POLH was mutated produced a lower number of OBs, which was presumably due to decreased POLH accumulation in the infected cells. Furthermore, cytoplasmic OBs were morphologically aberrant, even though nuclear OB morphology was normal in the same cell. These results indicate that NLS is required for nuclear localization and efficient accumulation of BmNPV POLH, which heavily affect the number and morphology of OBs.

Resistance mechanism of Nid-1, a dominant non-susceptibility gene, against Bombyx mori densovirus 1 infection

Bombyx mori densovirus 1 (BmDV1) is a pathogen that causes flacherie disease in mulberry silkworms (B. mori). The absolute resistance (non-susceptibility) to BmDV1 of certain silkworm strains is determined independently by two genes, nsd-1 and Nid-1. Previously, we investigated the expression of viral transcript in virus-inoculated silkworms carrying different nsd-1 and Nid-1 genotypes, and observed that nsd-1 and Nid-1 expression blocked the early and late steps of BmDV1 infection, respectively. In addition, we found that nsd-1 encoded a Bombyx-specific mucin-like membrane protein only present on the surface of the midgut, where BmDV1 could infect. In this study, we dissected the resistance mechanism by Nid-1 against BmDV1 infection by investigating the sequential changes in the accumulation of viral DNA, transcripts, and proteins derived from BmDV1 in susceptible strain (pxj) and Nid-1-carrying resistant strain (No. 908) after inoculation with BmDV1. Genomic PCR results showed that the BmDV1 DNA was detected immediately after the infection in both strains but rapidly decreased in the Nid-1-carrying strain No. 908 compared with the susceptible strain pxj. RT-PCR results also showed that the BmDV1 transcripts of Nid-1-carrying strain No. 908 were rapidly decreased after the infection. Moreover, BmDV1-derived proteins were not detected in No. 908 throughout the infection. These results suggest that Nid-1 expression might inhibit the accumulation of viral DNA and transcripts. As Nid-1 has not been molecularly characterized, its identification will contribute to the elucidation of the interactions between the silkworm and BmDV1.

SilkBase: an integrated transcriptomic and genomic database for Bombyx mori and related species

We introduce SilkBase as an integrated database for transcriptomic and genomic resources of the domesticated silkworm Bombyx mori and related species. SilkBase is the oldest B. mori database that was originally established as the expressed sequence tag database since 1999. Here, we upgraded the database by including the datasets of the newly assembled B. mori complete genome sequence, predicted gene models, bacterial artificial chromosome (BAC)-end and fosmid-end sequences, complementary DNA (cDNA) reads from 69 libraries, RNA-seq data from 10 libraries, PIWI-interacting RNAs (piRNAs) from 13 libraries, ChIP-seq data of 9 histone modifications and HP1 proteins and transcriptome and/or genome data of four B. mori-related species, i.e. Bombyx mandarina, Trilocha varians, Ernolatia moorei and Samia ricini. Our new integrated genome browser easily provides a snapshot of tissue- and stage-specific gene expression, alternative splicing, production of piRNAs and histone modifications at the gene locus of interest. Moreover, SilkBase is useful for performing comparative studies among five closely related lepidopteran insects.

Database URL: https://silkbase.ab.a.u-tokyo.ac.jp

Transcriptome analysis in the silkworm Bombyx mori overexpressing piRNA-resistant Masculinizer gene

Dosage compensation is a process that produces a similar expression of sex-linked and autosomal genes. In the silkworm Bombyx mori with a WZ sex-determination system, the expression from the single Z in WZ females matches that of ZZ males due to the suppression of Z-linked genes in males. A primary maleness determinant gene, Masculinizer (Masc), is also required for dosage compensation. In females, silkworm Piwi is complexed with the W chromosome-derived female-specific Feminizer (Fem) PIWI-interacting RNA (piRNA) and cleaves Masc mRNA. When Fem piRNA-resistant Masc cDNA (Masc-R) is overexpressed in both sexes, only female larvae are dead during the larval stage. In this study, transcriptome analysis was performed in neonate larvae to examine the effects of Masc-R overexpression on a global gene expression profile. Z-linked genes were globally repressed in Masc-R-overexpressing females due to force-driven dosage compensation. In contrast, Masc-R overexpression had little effect on the expression of Z-linked genes and the male-specific isoform of B. mori insulin-like growth factor II mRNA-binding protein in males, indicating that excessive Masc expression strengthens neither dosage compensation nor maleness in males. Fourteen genes were differentially expressed between Masc-R-overexpressing and control neonate larvae in both sexes, suggesting Masc functions other than dosage compensation and masculinization.

Mutations in a β-group of solute carrier gene are responsible for egg and eye coloration of the brown egg 4 (b-4) mutant in the silkworm, Bombyx mori

The brown egg 4 (b-4) is a recessive mutant in the silkworm (Bombyx mori), whose egg and adult compound eyes exhibit a reddish-brown color instead of normal purple and black, respectively. By double digest restriction-site associated DNA sequencing (ddRAD-seq) analysis, we narrowed down a region linked to the b-4 phenotype to approximately 1.1 Mb that contains 69 predicted gene models. RNA-seq analysis in a b-4 strain indicated that one of the candidate genes had a different transcription start site, which generates a short open reading frame. We also found that exon skipping was induced in the same gene due to an insertion of a transposable element in other two b-4 mutant strains. This gene encoded a putative amino acid transporter that belongs to the β-group of solute carrier (SLC) family and is orthologous to Drosophila eye color mutant gene, mahogany (mah). Accordingly, we named this gene Bmmah. We performed CRISPR/Cas9-mediated gene knockout targeting Bmmah. Several adult moths in generation 0 (G₀) had totally or partially reddish-brown compound eyes. We also established three Bmmah knockout strains, all of which exhibit reddish-brown eggs and adult compound eyes. Furthermore, eggs from complementation crosses between the b-4 mutants and the Bmmah knockout mutants also exhibited reddish-brown color, which was similar to the b-4 mutant eggs, indicating that Bmmah is responsible for the b-4 phenotypes.

Functional Characterization of Silkworm PIWI Proteins by Embryonic RNAi

The molecular mechanisms of sex-determination systems among insect orders and species are diverse. Therefore, genes involved in sex determination are strong candidates for insect pest management. Even though lepidopterans are major agricultural insect pests that cause widespread economic damage to various crops, their sex-determination systems have not been fully elucidated, even in the silkworm (Bombyx mori), a model lepidopteran insect. In 2014, we found that a female-specific W chromosome-derived PIWI-interacting RNA (piRNA) determines femaleness in silkworms. To analyze the function of two core silkworm piRNA biogenesis pathway genes, Siwi and BmAgo3, in the sex-determination system, we developed a genomic DNA and total RNA extraction strategy for a siRNA-injected single embryo. The siRNA-injected embryo can be molecularly sexed by W chromosome-specific DNA markers. Using complementary DNA (cDNA) reverse transcribed from the sexed RNA, we evaluated the knockdown effect of the PIWI protein-coding genes on a sexual development-related gene, Bombyx mori doublesex.

RNase κ promotes robust piRNA production by generating 2',3'-cyclic phosphate-containing precursors

In animal germlines, PIWI proteins and the associated PIWI-interacting RNAs (piRNAs) protect genome integrity by silencing transposons. Here we report the extensive sequence and quantitative correlations between 2',3'-cyclic phosphate-containing RNAs (cP-RNAs), identified using cP-RNA-seq, and piRNAs in the Bombyx germ cell line and mouse testes. The cP-RNAs containing 5'-phosphate (P-cP-RNAs) identified by P-cP-RNA-seq harbor highly consistent 5'-end positions as the piRNAs and are loaded onto PIWI protein, suggesting their direct utilization as piRNA precursors. We identified Bombyx RNase Kappa (BmRNase κ) as a mitochondria-associated endoribonuclease which produces cP-RNAs during piRNA biogenesis. BmRNase κ-depletion elevated transposon levels and disrupted a piRNA-mediated sex determination in Bombyx embryos, indicating the crucial roles of BmRNase κ in piRNA biogenesis and embryonic development. Our results reveal a BmRNase κ-engaged piRNA biogenesis pathway, in which the generation of cP-RNAs promotes robust piRNA production.

Bombyx mori nucleopolyhedrovirus ptp and egt genes are dispensable for triggering enhanced locomotory activity and climbing behavior in Bombyx mandarina larvae

Baculoviruses are classic pathogens that alter host behavior to enhance their dispersal and transmission. While viral protein tyrosine phosphatase (ptp) has been considered as a critical factor for inducing enhanced locomotory activity, preceding investigations have reported that viral ecdysteroid UDP-glucosyltransferase (egt) contributes to triggering climbing behavior in some virus and host species. Here we found that both egt and ptp were dispensable for these abnormal behaviors in Bombyx mandarina larvae induced by Bombyx mori nucleopolyhedrovirus, thus implying that there is an unknown core mechanism of baculovirus-induced alteration of host behaviors.

Loss of p24 from the Bombyx mori nucleopolyhedrovirus genome results in the formation of cuboidal occlusion bodies

Some insect viruses produce the occlusion body (OB), a large crystalline particle comprising a viral protein that occludes virions to protect them from harsh environments. The shapes and sizes of OBs are diverse depending on baculovirus species, but the detailed molecular mechanism determining them has yet to be totally clarified yet. Here we generated Bombyx mori nucleopolyhedrovirus (BmNPV) mutants of the p24 gene that encodes a viral capsid protein and found that p24-mutated BmNPVs produced cuboidal OBs with a slightly larger size than typical truncated octahedral OBs produced by wild-type BmNPVs. Meanwhile, p24 disruption has no significant impact on progeny virus production and viral pathogenicity. In addition, we experimentally demonstrated that a single amino acid substitution found in the P24 protein of the BmNPV Cubic isolate caused cuboidal OB production. These results suggest that p24 has a crucial role in generating the typical shape of OBs.

Horizontal Gene Transfer and Gene Duplication of β-Fructofuranosidase Confer Lepidopteran Insects Metabolic Benefits

Horizontal gene transfer (HGT) is a potentially critical source of material for ecological adaptation and the evolution of novel genetic traits. However, reports on posttransfer duplication in organism genomes are lacking, and the evolutionary advantages conferred on the recipient are generally poorly understood. Sucrase plays an important role in insect physiological growth and development. Here, we performed a comprehensive analysis of the evolution of insect β-fructofuranosidase transferred from bacteria via HGT. We found that posttransfer duplications of β-fructofuranosidase were widespread in Lepidoptera and sporadic occurrences of β-fructofuranosidase were found in Coleoptera and Hymenoptera. β-fructofuranosidase genes often undergo modifications, such as gene duplication, differential gene loss, and changes in mutation rates. Lepidopteran β-fructofuranosidase gene (SUC) clusters showed marked divergence in gene expression patterns and enzymatic properties in Bombyx mori (moth) and Papilio xuthus (butterfly). We generated SUC1 mutations in B. mori using CRISPR/Cas9 to thoroughly examine the physiological function of SUC. BmSUC1 mutant larvae were viable but displayed delayed growth and reduced sucrase activities that included susceptibility to the sugar mimic alkaloid found in high concentrations in mulberry. BmSUC1 served as a critical sucrase and supported metabolic homeostasis in the larval midgut and silk gland, suggesting that gene transfer of β-fructofuranosidase enhanced the digestive and metabolic adaptation of lepidopteran insects. These findings highlight not only the universal function of β-fructofuranosidase with a link to the maintenance of carbohydrate metabolism but also an underexplored function in the silk gland. This study expands our knowledge of posttransfer duplication and subsequent functional diversification in the adaptive evolution and lineage-specific adaptation of organisms.

Characterization of nuclear localization signal in Ostrinia furnacalis Masculinizer protein

Bombyx mori Masculinizer protein (BmMasc) is essential for both masculinization and dosage compensation in B. mori. We previously identified a bipartite nuclear localization signal (NLS) of BmMasc and two essential residues (lysine at 274 [K274] and arginine at 275 [R275]) implicated in its function. Sequence comparison showed the presence of putative NLSs in lepidopteran Masc proteins, but their functional properties and critical residues are unknown. Here we characterized a putative NLS of Ostrinia furnacalis Masc (OfMasc) using B. mori ovary-derived BmN-4 cell line. Deletion and alanine scanning mutagenesis revealed that a putative NLS is required for nuclear localization of OfMasc. However, mutations at both K227 and R228, which correspond to K274 and R275 of BmMasc, respectively, do not greatly abolish the NLS activity. Additional mutagenesis analysis revealed that triple mutations at K227, R228, and K240 almost completely inhibited OfMasc nuclear localization. These results suggest that lepidopteran Masc proteins possess a common functional NLS, but the critical residues for its activity are different. Moreover, we examined the masculinizing activity of OfMasc derivatives and found that nuclear localization is not required for the masculinizing activity of OfMasc. The results from our studies indicate that lepidopteran Masc proteins function in the cytoplasm to drive masculinizing cascade.

Potential for small RNA production against Bombyx mori latent virus in Bombyx mori ovaries

Bombyx mori latent virus (BmLV) is a positive, single-stranded insect RNA virus closely related to plant maculaviruses. BmLV was first isolated from Bombyx mori ovary-derived cell line BmN-4, and this virus has already infected most B. mori-derived cultured cell lines. We previously reported that small interfering RNA (siRNA) and PIWI-interacting RNA (piRNA) pathways function cooperatively to maintain the amount of BmLV RNA for normal BmN-4 cell growth. On the other hand, BmLV does not propagate in B. mori larvae. Here we conducted BmLV injection into the larval body cavities of B. mori, and examined BmLV accumulation in larval ovaries where siRNA and piRNA pathways are both active, to investigate whether this in vivo resistance is governed by small RNA pathways. Expression levels of RNA-dependent RNA polymerase, coat protein, and p15 genes in BmLV-injected larval ovaries were extremely low compared with those in B. mori cultured cells, indicating that B. mori larval ovaries are more resistant to BmLV than B. mori cultured cells. We also sequenced small RNAs prepared from BmLV-injected larval ovaries and mapped them onto the BmLV genome. Although their amounts were very small, we were able to detect BmLV-derived small RNAs in the ovaries. According to their length distribution and nucleotide bias, they were likely to be siRNAs and piRNAs. These results suggest that B. mori ovaries can potentially produce small RNAs against BmLV, but the resistance of larval ovaries against BmLV is not dependent on RNA silencing pathways.

Formation of the CenH3-Deficient Holocentromere in Lepidoptera Avoids Active Chromatin

Despite the essentiality for faithful chromosome segregation, centromere architectures are diverse among eukaryotes^1,2 and embody two main configurations: mono- and holocentromeres, referring, respectively, to localized or unrestricted distribution of centromeric activity. Of the two, some holocentromeres offer the curious condition of having arisen independently in multiple insects, most of which have lost the otherwise essential centromere-specifying factor CenH3³ (first described as CENP-A in humans).^4-7 The loss of CenH3 raises intuitive questions about how holocentromeres are organized and regulated in CenH3-lacking insects. Here, we report the first chromatin-level description of CenH3-deficient holocentromeres by leveraging recently identified centromere components^6,7 and genomics approaches to map and characterize the holocentromeres of the silk moth Bombyx mori, a representative lepidopteran insect lacking CenH3. This uncovered a robust correlation between the distribution of centromere sites and regions of low chromatin activity along B. mori chromosomes. Transcriptional perturbation experiments recapitulated the exclusion of B. mori centromeres from active chromatin. Based on reciprocal centromere occupancy patterns observed along differentially expressed orthologous genes of Lepidoptera, we further found that holocentromere formation in a manner that is recessive to chromatin dynamics is evolutionarily conserved. Our results help us discuss the plasticity of centromeres in the context of a role for the chromosome-wide chromatin landscape in conferring centromere identity rather than the presence of CenH3. Given the co-occurrence of CenH3 loss and holocentricity in insects,⁷ we further propose that the evolutionary establishment of holocentromeres in insects was facilitated through the loss of a CenH3-specified centromere.

Hsp90 function is required for stable transcription of the baculovirus transactivator ie-1 gene

A molecular chaperone heat shock protein 90 (Hsp90) is required for efficient infection by several viruses. Hsp90 has been recently implicated in baculovirus infection, but its exact role remains obscure. This study investigated the effect of 17-N-allylamino-17-demethoxygeldanamycin (17-AAG), an Hsp90-specific inhibitor, on Bombyx mori nucleopolyhedrovirus (BmNPV) infection. The 17-AAG treatment significantly decreased the production of budded viruses and occlusion bodies in BmNPV-infected Bombyx mori cultured cells. Immunoblot and SDS-PAGE analyses showed that the expression of early and delayed early gene products, DBP and BRO, was delayed and dysregulated, and the very late gene product POLH was almost completely diminished. RT-qPCR experiments revealed that 17-AAG treatment did not affect initiation of the immediate early gene ie-1 expression, but the expression decreased by ∼50 % during the late stage of infection. 17-AAG treatment also decreased ie-1 promoter activity by ∼50 %. In addition, the expression of delayed early and late genes was dysregulated and inhibited, respectively. These results indicated that Hsp90 function is required for stable ie-1 transcription. Inhibiting Hsp90 function negatively affects ie-1 expression, resulting in dysregulation of delayed early genes and a severe decrease in late and very late gene expression.

The genome sequence of Samia ricini, a new model species of lepidopteran insect

Samia ricini, a gigantic saturniid moth, has the potential to be a novel lepidopteran model species. Samia ricini is far more resistant to diseases than the current model species Bombyx mori, and therefore can be more easily reared. In addition, genetic resources available for S. ricini rival those for B. mori: at least 26 ecoraces of S. ricini are reported and S. ricini can hybridize with wild Samia species, which are distributed throughout Asian countries, and produce fertile progenies. Physiological traits such as food preference, integument colour and larval spot pattern differ among S. ricini strains and wild Samia species so that those traits can be targeted in forward genetic analyses. To facilitate genetic research in S. ricini, we determined its whole genome sequence. The assembled genome of S. ricini was 458 Mb with 155 scaffolds, and the scaffold N50 length of the assembly was ~ 21 Mb. In total, 16,702 protein coding genes were predicted. While the S. ricini genome was mostly collinear with that of B. mori with some rearrangements and few S. ricini-specific genes were discovered, chorion genes and fibroin genes seemed to have expanded in the S. ricini lineage. As the first step of genetic analyses, causal genes for "Blue," "Yellow," "Spot," and "Red cocoon" phenotypes were mapped to chromosomes.

Pheromonal activities of the bombykol isomer, (10E,12E)-10,12-hexadecadien-1-ol, in the pheromone gland of the silkmoth Bombyx mori

Bombykol (EZ) is the single component of the female sex pheromone in the silkmoth Bombyx mori. EZ alone evokes full courtship behaviors from conspecific males; however, its geometric isomer (EE) was consistently detected in the pheromone glands (PG) of 16 B. mori strains and a field population of the wild silkmoth Bombyx mandarina, which also uses EZ as the single pheromone component. We investigated the pheromonal activities of EE using a commercial hybrid strain of B. mori, Kinshu × Showa. The behavioral assay demonstrated that a 10⁴-10⁵-fold higher dose of EE than EZ was able to elicit behavioral responses from males. To elucidate whether the trace contaminant of EZ in the EE standard is responsible for these responses, we examined the responses of male antennae to EE using a gas chromatograph-electroantennographic detector system (GC-EAD). The EE, at high doses elicited marginal responses from the male antennae. We next examined antennal and behavioral responses of B. mori whose BmOR1 gene, which is responsible for the reception of bombykol, was knocked out. The knockout of BmOR1 resulted in the complete loss of antennal and behavioral responses to EE and EZ, demonstrating that if EE itself is active, it induces these responses via the incidental stimulation of BmOR1, not via the stimulation of EE-specific receptors. The existence of EE in the PG of B. mori and B. mandarina is discussed from the viewpoints of pheromone biosynthesis and the evolution of pheromone communication systems.

Duplication and diversification of trehalase confers evolutionary advantages on lepidopteran insects

Gene duplication provides a major source of new genes for evolutionary novelty and ecological adaptation. However, the maintenance of duplicated genes and their relevance to adaptive evolution has long been debated. Insect trehalase (Treh) plays key roles in energy metabolism, growth, and stress recovery. Here, we show that the duplication of Treh in Lepidoptera (butterflies and moths) is linked with their adaptation to various environmental stresses. Generally, two Treh genes are present in insects: Treh1 and Treh2. We report three distinct forms of Treh in lepidopteran insects, where Treh1 was duplicated into two gene clusters (Treh1a and Treh1b). These gene clusters differ in gene expression patterns, enzymatic properties, and subcellular localizations, suggesting that the enzymes probably underwent sub- and/or neofunctionalization in the lepidopteran insects. Interestingly, selective pressure analysis provided significant evidence of positive selection on duplicate Treh1b gene in lepidopteran insect lineages. Most positively selected sites were located in the alpha-helical region, and several sites were close to the trehalose binding and catalytic sites. Subcellular adaptation of duplicate Treh1b driven by positive selection appears to have occurred as a result of selected changes in specific sequences, allowing for rapid reprogramming of duplicated Treh during evolution. Our results suggest that gene duplication of Treh and subsequent functional diversification could increase the survival rate of lepidopteran insects through various regulations of intracellular trehalose levels, facilitating their adaptation to diverse habitats. This study provides evidence regarding the mechanism by which gene family expansion can contribute to species adaptation through gene duplication and subsequent functional diversification.

Masc-induced dosage compensation in silkworm cultured cells

The Masculinizer (Masc) gene encodes a CCCH‐tandem zinc finger protein that controls both masculinization and dosage compensation in the silkworm Bombyx mori. We previously measured the masculinizing activity of the lepidopteran Masc proteins using B. mori ovary‐derived cell line BmN‐4. Here, we established an RNA‐seq data‐based assay system in which the level of B. mori Masc (BmMasc)‐induced dosage compensation can be estimated in BmN‐4 cells. Using this system, we found that a cysteine residue at position 301, which was shown to be essential for the masculinizing activity of BmMasc, is also required for dosage compensation. We further investigated the relationships between Masc‐induced cell growth inhibition, masculinizing activity, and the level of dosage compensation, using Masc genes from three lepidopteran insects. In summary, we have established a cell‐based system to monitor levels of Masc‐induced dosage compensation.

High-quality genome assembly of the silkworm, Bombyx mori

In 2008, the genome assembly and gene models for the domestic silkworm, Bombyx mori, were published by a Japanese and Chinese collaboration group. However, the genome assembly contains a non-negligible number of misassembled and gap regions due to the presence of many repetitive sequences within the silkworm genome. The erroneous genome assembly occasionally causes incorrect gene prediction. Here we performed hybrid assembly based on 140 × deep sequencing of long (PacBio) and short (Illumina) reads. The remaining gaps in the initial genome assembly were closed using BAC and Fosmid sequences, giving a new total length of 460.3 Mb, with 30 gap regions and an N50 comprising 16.8 Mb in scaffolds and 12.2 Mb in contigs. More RNA-seq and piRNA-seq reads were mapped on the new genome assembly compared with the previous version, indicating that the new genome assembly covers more transcribed regions, including repetitive elements. We performed gene prediction based on the new genome assembly using available mRNA and protein sequence data. The number of gene models was 16,880 with an N50 of 2154 bp. The new gene models reflected more accurate coding sequences and gene sets than old ones. The proportion of repetitive elements was also reestimated using the new genome assembly, and was calculated to be 46.8% in the silkworm genome. The new genome assembly and gene models are provided in SilkBase (http://silkbase.ab.a.u-tokyo.ac.jp).

Two CCCH-type zinc finger domains in the Masc protein are dispensable for masculinization and dosage compensation in Bombyx mori

The Masculinizer (Masc) gene encodes a novel lepidopteran-specific protein that controls both masculinization and dosage compensation in the silkworm Bombyx mori. The Masc protein possesses two CCCH-type zinc finger domains (ZFs), a nuclear localization signal, and an 11-amino-acid region that is highly conserved among lepidopteran insects. Using a cell-based assay system, we revealed that two cysteine residues localized in the conserved region, but not ZFs, are required for masculinization. In addition, nuclear localization of the Masc protein is not associated with masculinizing activity. Because dosage compensation is considered to occur in the nucleus, we inferred that the two ZFs play a role in the establishment of dosage compensation. To investigate this hypothesis at the organism level, we utilized the CRISPR/Cas9 system and established three B. mori strains whose Masc is partially deleted at different regions. The strain lacking the 210 C-terminal amino acids of the Masc protein showed male-specific embryonic lethality due to its low abundance and/or instability. The male embryos of this strain expressed the female-type splice variants of B. mori doublesex and did not express the male-type mRNA of B. mori IGF-II mRNA-binding protein. Furthermore, mRNA levels of Z-linked genes were abnormally enhanced only in male embryos. In contrast, the strain lacking both ZFs grew normally and did not show any defective phenotypes including sexual differentiation and the expression of Z-linked genes, demonstrating that the two CCCH-type ZFs, which are conserved in lepidopteran Masc homologs, are dispensable for masculinization and dosage compensation.

A reexamination on the deficiency of riboflavin accumulation in Malpighian tubules in larval translucent mutants of the silkworm, Bombyx mori

A variety of insects accumulate high contents of riboflavin (vitamin B₂) in their Malpighian tubules (MTs). Although this process is known to be genetically controlled, the mechanism is not known. In the 1940s and the 1950s, several studies showed that riboflavin contents were low in the MTs of some Bombyx mori (silkworm) mutants with translucent larval skin mutations (e.g., w-3, od, oa, and otm) and that genes responsible for these translucent mutations also affected riboflavin accumulation in the MTs. Since the 2000s, it has been shown that the w-3 gene encodes an ABC transporter, whereas genes responsible for od, oa, and otm mutations encode for the biogenesis of lysosome-related organelles. These findings suggest that some genes of ABC transporters and biogenesis of lysosome-related organelles may control the accumulation of riboflavin in MTs. Therefore, we reexamined the effects that translucent mutations have on the accumulation of riboflavin in MTs by using the translucent and wild-type segregants in mutant strains to measure the specific effect that each gene has on riboflavin accumulation (independent of genomic background). We used nine translucent mutations (w-3^oe, oa, od, otm, Obs, oy, or, oh, and obt) even though the genes responsible for some of these mutations (Obs, oy, or, oh, and obt) have not yet been isolated. Through observation of larval MTs and measurements of riboflavin content using high-performance liquid chromatography, we found that the oa, od, otm, and or mutations were responsible for low contents of riboflavin in MTs, whereas the Obs and oy mutations did not affect riboflavin accumulation. This indicates that the molecular mechanism for riboflavin accumulation is similar but somewhat different than the mechanism responsible for uric acid accumulation in epidermal cells. We found that the genes responsible for oa, od, and otm mutations were consistent with those already established for uric acid accumulation in larval epidermis. This suggests that these three genes control riboflavin accumulation in MTs through a mechanism similar to that of uric acid accumulation, although we do not yet know why the or mutation also controls riboflavin accumulation.

Inhibitory role of the Bm8 protein in the propagation of Bombyx mori nucleopolyhedrovirus

Lepidopteran nucleopolyhedroviruses have distinct viral tissue tropisms in host larvae. We previously identified the Bm8 gene of Bombyx mori nucleopolyhedrovirus (BmNPV), the product of which inhibits viral propagation in the middle silk gland (MSG). However, it is unknown whether this inhibitory function of the Bm8 protein is specific to MSGs. Here we generated a Bm8-disrupted recombinant BmNPV expressing green fluorescent protein (GFP) and examined viral propagation in B. mori cultured cells and larvae. We found that Bm8-disrupted BmNPV produced fewer budded viruses and more occlusion bodies (OBs) than the wild-type virus in both cultured cells and larvae. Microscopic observation of OB production and GFP expression revealed that Bm8 disruption accelerated the progression of viral infection in various larval tissues. Furthermore, quantitative reverse transcription-polymerase chain reaction experiments showed that the loss of Bm8 enhanced viral gene expression in BmNPV-infected larval tissues. These results indicate that the Bm8 protein suppresses viral propagation to varying degrees in each larval tissue, which may establish BmNPV tissue tropisms in B. mori larvae.

Transcriptome profiling reveals infection strategy of an insect maculavirus

Bombyx mori macula-like virus (BmMLV) is a positive, single-stranded insect RNA virus that is closely related to plant maculaviruses. BmMLV is currently characterized as an unclassified maculavirus. BmMLV accumulates at extremely high levels in cell lines derived from the silkworm, Bombyx mori, but it does not lead to lethality and establishes persistent infections. It is unknown how this insect maculavirus replicates and establishes persistent infections in insect cells. Here, we showed that BmMLV p15, which is located on a subgenomic fragment and is not found in plant maculaviruses, is highly expressed in BmMLV-infected silkworm cells and that p15 protein is required to establish BmMLV infections in silkworm cells. We also showed that two distinct small RNA-mediated pathways maintain BmMLV levels in BmMLV-infected silkworm cells, thereby allowing the virus to establish persistent infection. Virus-derived siRNAs and piRNAs were both produced as the infection progressed. Knockdown experiments demonstrated that the exogenous RNAi pathway alone or RNAi and piRNA pathways function cooperatively to silence BmMLV RNA and that both pathways are important for normal growth of BmMLV-infected silkworm cells. On the basis of our study, we propose a mechanism of how a plant virus-like insect virus can establish persistent infections in insect cells.

Bombyx ortholog of the Drosophila eye color gene brown controls riboflavin transport in Malpighian tubules

The Drosophila eye color gene brown is known to control the transport of pteridine precursors in adult eyes. The Brown protein belongs to the ATP-binding cassette (ABC) transporter G family, which includes proteins encoded by the genes brown, scarlet, and white. These genes are responsible for pigmentation in Drosophila and the domestic silkworm Bombyx mori. Although orthologs of brown are conserved among insects, the function of this gene is only known in Drosophila. Here, we elucidated the function of the B. mori ortholog Bm-brown. We examined the spatial and temporal expression profiles of Bm-brown and found that this gene was specifically and continuously expressed in larval Malpighian tubules (MTs), indicating this gene has a special function in MTs. We then successfully obtained a Bm-brown knockout (KO) strain based on a wild-type (WT) strain using the clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated nuclease 9 (Cas9) system. We found that larval MTs of the KO strain were white, whereas those of WT were yellow. It is known that larval yellow MTs of WT are due to the accumulation of riboflavin. Therefore, we compared the riboflavin contents of MTs of KO and WT strains, and found that the riboflavin level in the KO strain was 20 fold less than that in WT during the 5th instar period. MTs are known to exhibit a similar milky color in w-3 mutant larvae due to a deficiency of riboflavin accumulation. The responsible gene for w-3 mutant is the Bmwh3 gene, which is orthologous to Drosophila white. Thus, we speculate that Bm-brown is heterodimerized with Bmwh3, similar to Brown/White in Drosophila, and acts as a riboflavin transporter in silkworm MTs.

Artificial "ping-pong" cascade of PIWI-interacting RNA in silkworm cells

PIWI-interacting RNAs (piRNAs) play essential roles in the defense system against selfish elements in animal germline cells by cooperating with PIWI proteins. A subset of piRNAs is predicted to be generated via the "ping-pong" cascade, which is mainly controlled by two different PIWI proteins. Here we established a cell-based artificial piRNA production system using a silkworm ovarian cultured cell line that is believed to possess a complete piRNA pathway. In addition, we took advantage of a unique silkworm sex-determining one-to-one ping-pong piRNA pair, which enabled us to precisely monitor the behavior of individual artificial piRNAs. With this novel strategy, we successfully generated artificial piRNAs against endogenous protein-coding genes via the expected back-and-forth traveling mechanism. Furthermore, we detected "primary" piRNAs from the upstream region of the artificial "ping-pong" site in the endogenous gene. This artificial piRNA production system experimentally confirms the existence of the "ping-pong" cascade of piRNAs. Also, this system will enable us to identify the factors involved in both, or each, of the "ping" and "pong" cascades and the sequence features that are required for efficient piRNA production.

Bombyx mori nucleopolyhedrovirus BM5 protein regulates progeny virus production and viral gene expression

Bombyx mori nucleopolyhedrovirus (BmNPV) orf5 (Bm5) is a core gene of lepidopteran baculoviruses and encodes the protein with the conserved amino acid residues (DUF3627) in its C-terminus. Here, we found that Bm5 disruption resulted in lower titers of budded viruses and fewer numbers of occlusion bodies (OBs) in B. mori cultured cells and larvae, although viral genome replication was not affected. Bm5 disruption also caused aberrant expression of various viral genes at the very late stage of infection. Immunocytochemical analysis revealed that BM5 localized to the nuclear membrane. We also found that DUF3627 is important for OB production, transcriptional regulation of viral genes, and subcellular localization of BM5. Compared with wild-type BmNPV infection, larval death was delayed when B. mori larvae were infected with Bm5 mutants. These results suggest that BM5 is involved in progeny virus production and regulation of viral gene expression at the very late stage of infection.

Transgenic Expression of the piRNA-Resistant Masculinizer Gene Induces Female-Specific Lethality and Partial Female-to-Male Sex Reversal in the Silkworm, Bombyx mori

In Bombyx mori (B. mori), Fem piRNA originates from the W chromosome and is responsible for femaleness. The Fem piRNA-PIWI complex targets and cleaves mRNAs transcribed from the Masc gene. Masc encodes a novel CCCH type zinc-finger protein and is required for male-specific splicing of B. mori doublesex (Bmdsx) transcripts. In the present study, several silkworm strains carrying a transgene, which encodes a Fem piRNA-resistant Masc mRNA (Masc-R), were generated. Forced expression of the Masc-R transgene caused female-specific lethality during the larval stages. One of the Masc-R strains weakly expressed Masc-R in various tissues. Females heterozygous for the transgene expressed male-specific isoform of the Bombyx homolog of insulin-like growth factor II mRNA-binding protein (ImpM) and Bmdsx. All examined females showed a lower inducibility of vitellogenin synthesis and exhibited abnormalities in the ovaries. Testis-like tissues were observed in abnormal ovaries and, notably, the tissues contained considerable numbers of sperm bundles. Homozygous expression of the transgene resulted in formation of the male-specific abdominal segment in adult females and caused partial male differentiation in female genitalia. These results strongly suggest that Masc is an important regulatory gene of maleness in B. mori.

Identification and Functional Analysis of the Pre-piRNA 3' Trimmer in Silkworms

PIWI-interacting RNAs (piRNAs) play a crucial role in transposon silencing in animal germ cells. In piRNA biogenesis, single-stranded piRNA intermediates are loaded into PIWI-clade proteins and cleaved by Zucchini/MitoPLD, yielding precursor piRNAs (pre-piRNAs). Pre-piRNAs that are longer than the mature piRNA length are then trimmed at their 3' ends. Although recent studies implicated the Tudor domain protein Papi/Tdrkh in pre-piRNA trimming, the identity of Trimmer and its relationship with Papi/Tdrkh remain unknown. Here, we identified PNLDC1, an uncharacterized 3'-5' exonuclease, as Trimmer in silkworms. Trimmer is enriched in the mitochondrial fraction and binds to Papi/Tdrkh. Depletion of Trimmer and Papi/Tdrkh additively inhibits trimming, causing accumulation of ∼35-40-nt pre-piRNAs that are impaired for target cleavage and prone to degradation. Our results highlight the cooperative action of Trimmer and Papi/Tdrkh in piRNA maturation.

Morphological and histomorphological structures of testes and ovaries in early developmental stages of the silkworm, Bombyx mori

The gonad develops as a testis in male or an ovary in female. In the silkworm, B. mori, little is known about testis and ovary in the embryonic stages and early larval stages. In this study, we performed morphological and histomorphological observations of ovaries and testes from the late embryonic stage to the 1st instar larval stage. Results obtained with lack of accurate information on sex of examined individuals may be misleading, thus we performed phenotypic observations of gonads by utilizing sex-limited strain that enables us to easily discriminate female embryos from male ones based on those egg colors. In testis, four testicular follicles were clearly observed in the testis at the first instar larval stage, and boundary layers were formed between the testicular follicles. At the late embryonic stage, the testis consisted of four testicular follicles, while the boundary layers were still obscure. In ovary, four ovarioles were easily recognizable in the ovary at the first instar larval stage, and boundary layers were formed between the ovarioles. However, in the late embryonic stage, it was quite difficult to identify four ovarioles. Morphological characteristics were almost similar between testis and ovary in early developmental stages. Our present study demonstrates that the most reliable difference between testis and ovary in early developmental stages is the attaching point of the duct. Formation and development of the duct may be sensitive to the sex-determining signal and display sexual dimorphism in early embryonic stages.

Identification of the silkworm quail gene reveals a crucial role of a receptor guanylyl cyclase in larval pigmentation

Diverse color patterns on the integument of lepidopteran larvae play important roles in their survival through camouflage, mimicry, sexual signaling, and aposematism. In the silkworm Bombyx mori, many color pattern variations have been preserved in inbred strains making them a good model for elucidating the molecular mechanisms that underlie color pattern formation. In this study, we focused on the silkworm quail (q) mutant, which exhibits abnormalities in multiple pigment biosynthesis pathways. Positional cloning of the q gene revealed that disruption of a guanylyl cyclase gene, BmGC-I, is responsible for its abnormal pigmentation. In q mutants, we identified a 16-bp deletion in the BmGC-I transcript, resulting in the production of a premature stop codon. Knockout of the BmGC-I gene resulted in the q-like abnormal pigmentation, thereby demonstrating that the BmGC-I gene is involved in the pigment biosynthesis pathway in the integument. Moreover, quantitative reverse transcription polymerase chain reaction showed that BmGC-I was strongly expressed in the fourth instar on day 2. Our results suggest that BmGC-I deficiency affects the pigment biosynthesis pathway, which supports the involvement of guanylyl cyclase in larval coloration.

Phosphatase activity of Bombyx mori nucleopolyhedrovirus PTP is dispensable for enhanced locomotory activity in B. mori larvae

Baculovirus-induced enhanced locomotory activity (ELA) is not induced in caterpillars infected with a mutant Bombyx mori nucleopolyhedrovirus (BmNPV) or Autographa californica multiple nucleopolyhedrovirus (AcMNPV) lacking a functional protein tyrosine phosphatase gene (ptp). Previous studies suggest that the PTP proteins from BmNPV and AcMNPV act in different ways to induce ELA, i.e., BmNPV PTP is utilized as a virion structural component, whereas AcMNPV PTP requires its phosphatase activity. Here, I generated and characterized two new BmNPV mutants expressing enzymatically inactive PTP proteins and confirmed that the phosphatase activity of PTP is not required for ELA induction in BmNPV-infected B. mori larvae.

Inactivation of Bombyx mori macula-like virus under physical conditions

The Bombyx mori macula-like virus (BmMLV) is a member of the genus Maculavirus, family Tymoviridae, and contains a positive-sense single-stranded RNA genome. Previously, we reported that almost all B. mori-derived cell lines have already been contaminated with BmMLV via an unknown infection route. Since B. mori-derived cell lines are used for the baculovirus expression vector system, the invasion of BmMLV will cause a serious safety risk in the production of recombinant proteins. In this study, to determine the inactivation effectiveness of BmMLV, viruses were treated with various temperatures as well as gamma and ultraviolet (UV) light radiation. After these treatments, the virus solutions were inoculated into BmMLV-free BmVF cells. At 7 days postinoculation, the amount of virus in cells was evaluated by real-time reverse transcription PCR. Regarding heat treatment, conditions under 56°C for 3 h were tolerated, whereas infectivity disappeared after treatment at 75°C for 1 h. Regarding gamma radiation treatment, viruses were relatively stable at 1 kGy; however, their infectivity was entirely eliminated at a dose of 10 kGy. With 254 nm UV-C treatment, viruses were still active at less than 120 mJ/cm(2); however, their infectivity was completely lost at greater than 140 mJ/cm(2) UV-C radiation. These results provide quantitative evidence of the potential for BmMLV inactivation under a variety of physical conditions.