A Hox Gene, Antennapedia, Regulates Expression of Multiple Major Silk Protein Genes in the Silkworm Bombyx mori

Differential expression of fibroin-related genes in middle silk glands is induced by dietary differences in a strain-dependent manner in Bombyx mori

The silkworm (Bombyx mori) is a model organism for lepidopteran insects. It is an oligophagous insect that primarily feeds on mulberry leaves and has industrial use for the production of raw silk. The development of artificial diets has provided an alternative nutrient source for silkworms; however, one significant issue is that the production of cocoons is lower in silkworms reared on artificial diets compared with those reared on mulberry leaves. The differences in the silk gland in the late-stage fifth instar silkworm larvae, when silk synthesis is most active, between those raised on artificial diets and mulberry leaves, are unknown. In this study, we identified differences in the transcriptomes of the middle and posterior silk glands of fifth instar day five silkworm larvae reared on artificial diets compared with those reared on mulberry leaves using three strains: Daizo, Nichi01, and J137 × C146. We found that the silk-related genes fibrohexamerin (fhx), fibroin-light-chain (fibL), and fibroin-heavy-chain (fibH) in the middle silk gland, and ser1 in the posterior silk gland, were differentially expressed in a strain-dependent manner. In silkworms reared on artificial diets, fhx, fibL, and fibH in the middle silk gland were upregulated in Nichi01 and downregulated in J137 × C146, whereas ser1 in the posterior silk gland was upregulated in J137 × C146 compared with silkworms reared on mulberry leaves. Our results demonstrate that the diet and strain of silkworm larvae affect the expression of genes related to silk production in their silk glands during the late fifth instar stage.

Dynamic chromatin conformation and accessibility changes mediate the spatial-specific gene regulatory network in Bombyx mori

Silk gland genes of Bombyx mori can have strict spatial expression patterns, which impact their functions and silk quality; however, our understanding of their regulation mechanisms is currently insufficient. To address this, the middle silk gland (MSG) and posterior silk gland (PSG) of the silkworm were investigated. Gene ontology annotation showed that spatially specific expressed genes were involved in the formation of H3k9me and chromatin topology. Chromatin conformation data generated by Hi-C showed that the topologically associated domain boundaries around FibL and Sericin1 genes were significantly different between MSG and PSG. Changes in chromatin conformation led to changes in chromatin activity, which significantly affected the expression of nearby genes in silkworm. Chromatin accessibility regions of MSG and PSG were analyzed using FAIRE-seq, and 1006 transcription factor motifs were identified in open chromatin regions. Furthermore, the spatial-specific expression patterns of silk gland genes were mainly associated with homeobox-contained transcription factors, such as POU-M2, which was specifically bound and relatively highly expressed in the MSG. The regulatory network mediated by POU-M2 regulated most of the spatial-specific expressed genes in MSG, such as ADH1. These results can aid in improving silk performance, optimizing silkworm breeding, and improving the gene spatial regulatory model research for insects.

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

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

Homeodomain proteins POU-M2, antennapedia and abdominal-B are involved in regulation of the segment-specific expression of the clip-domain serine protease gene CLIP13 in the silkworm, Bombyx mori

Clip-domain serine proteases (CLIPs) play important roles in insect innate immunity and development. Our previous studies indicated that CLIP13, an epidermis-specific gene, was involved in cuticle remodeling during molting and metamorphosis in the silkworm, Bombyx mori. However, the transcriptional regulatory mechanism and regulatory pathways of CLIP13 remained unclear. In the present study, we investigated CLIP13 expression and the regulation pathway controlled by 20-hydroxyecdysone (20E) in the silkworm. At the transcriptional level, expression of CLIP13 exhibited pronounced spatial and temporal specificity in different regions of the epidermis; homeodomain transcription factors POU-M2, antennapedia (Antp), and abdominal-B (Abd-B) showed similar expression change trends as CLIP13 in the head capsule, thorax, and abdomen, respectively. Furthermore, results of cell transfection assays, electrophoretic mobility shift assays, and chromatin immunoprecipitation demonstrated that POU-M2, Antp, and Abd-B were involved in the transcriptional regulation of CLIP13 by directly binding to their cis-response elements in CLIP13 promoter. RNA interference-mediated silencing of POU-M2, Antp, and Abd-B led to a decrease of CLIP13 expression in the head capsule, the epidermis of the 1st to 3rd thoracic segments and the 7th to 10th abdominal segments, respectively. Consistent with CLIP13, 20E treatment significantly upregulated expression of POU-M2, Antp, and Abd-B in the silkworm epidermis. Taken together, these data suggest that 20E positively regulates transcription of CLIP13 via homeodomain proteins POU-M2, Antp, and Abd-B in different regions of the silkworm epidermis during metamorphosis, thus affecting the molting process. Our findings provide new insight into the functions of homeodomain transcription factors in insect molting.

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

Simple Summary

Previous studies have reported how the silk production ability of Bombyx mori can be enhanced, but the mechanism that regulates silk protein genes remains unclear. We performed co-expression network analysis using networkz, an in-house program, which led to the identification of 91 transcription factors were co-expressed with silk protein genes. Of them, 13 transcripts were identified to be novel regulatory factors by time-course expression analysis during the fifth instar larvae stage. Their expression patterns were highly relevant to those of silk protein genes. Our results suggest that the two-step expression screening was robust and highly sensitive to screen relative genes, and a complex mechanism regulates silk protein production in B. mori. The novel candidates that were identified herein can serve as key genes to develop methods to enhance the silk protein production ability of B. mori.

Abstract

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

Transcriptomic analysis of the bagworm moth silk gland reveals a number of silk genes conserved within Lepidoptera

Lepidopteran insects produce cocoons with unique properties. The cocoons are made of silk produced in the larval tissue silk gland and our understanding of the silk genes is still very limited. Here, we investigated silk genes in the bagworm moth Eumeta variegata, a species that has recently been found to produce extraordinarily strong and tough silk. Using short-read transcriptomic analysis, we identified a partial sequence of the fibroin heavy chain gene and its product was found to have a C-terminal structure that is conserved within nonsaturniid species. This is in accordance with the presence of fibroin light chain/fibrohexamerin genes and it is suggested that the bagworm moth is producing silk composed of fibroin ternary complex. This indicates that the fibroin structure has been evolutionarily conserved longer than previously thought. Other than fibroins we identified candidates for sericin genes, expressed strongly in the middle region of the silk gland and encoding serine-rich proteins, and other silk genes, that are structurally conserved with other lepidopteran homologues. The bagworm moth is thus considered to be producing conventional lepidopteran type of silk. We further found a number of genes expressed in a specific region of the silk gland and some genes showed conserved expression with Bombyx mori counterparts. This is the first study allowing comprehensive silk gene identification and expression analysis in the lepidopteran Psychidae family and should contribute to the understanding of silk gene evolution as well as to the development of novel types of silk.

A cis-Acting Mutation in the PxABCG1 Promoter Is Associated with Cry1Ac Resistance in Plutella xylostella (L.)

The molecular mechanisms of insect resistance to Cry toxins generated from the bacterium Bacillus thuringiensis (Bt) urgently need to be elucidated to enable the improvement and sustainability of Bt-based products. Although downregulation of the expression of midgut receptor genes is a pivotal mechanism of insect resistance to Bt Cry toxins, the underlying transcriptional regulation of these genes remains elusive. Herein, we unraveled the regulatory mechanism of the downregulation of the ABC transporter gene PxABCG1 (also called Pxwhite), a functional midgut receptor of the Bt Cry1Ac toxin in Plutella xylostella. The PxABCG1 promoters of Cry1Ac-susceptible and Cry1Ac-resistant strains were cloned and analyzed, and they showed clear differences in activity. Subsequently, a dual-luciferase reporter assay, a yeast one-hybrid (Y1H) assay, and RNA interference (RNAi) experiments demonstrated that a cis-mutation in a binding site of the Hox transcription factor Antennapedia (Antp) decreased the promoter activity of the resistant strain and eliminated the binding and regulation of Antp, thereby enhancing the resistance of P. xylostella to the Cry1Ac toxin. These results advance our knowledge of the roles of cis- and trans-regulatory variations in the regulation of midgut Cry receptor genes and the evolution of Bt resistance, contributing to a more complete understanding of the Bt resistance mechanism.

Reference Transcriptome Data in Silkworm Bombyx mori

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

Cell type innovation at the tips of the animal tree

Understanding how organs originate is challenging due to the twin problems of explaining how new cell types evolve and how collective interactions between cell types arise and become selectively advantageous. Animals are assemblages of organs and cell types of different antiquities, and among the most rapidly and convergently evolving are exocrine glands and their constituent secretory cell types. Such structures have arisen independently thousands of times across the Metazoa, impacting how animals chemically interact with their environments. The recurrent evolution of exocrine systems provides a paradigm for examining how qualitative phenotypic novelties arise from variation at the cellular level. Here, we take a hierarchical perspective, focusing on the evolutionary assembly of novel biosynthetic pathways and secretory cell types, and how both selection and non-adaptive molecular processes may combine to build the complex, modular architectures of many animal glands.

Molecular evolution of gland cell types and chemical interactions in animals

Across the Metazoa, the emergence of new ecological interactions has been enabled by the repeated evolution of exocrine glands. Specialized glands have arisen recurrently and with great frequency, even in single genera or species, transforming how animals interact with their environment through trophic resource exploitation, pheromonal communication, chemical defense and parental care. The widespread convergent evolution of animal glands implies that exocrine secretory cells are a hotspot of metazoan cell type innovation. Each evolutionary origin of a novel gland involves a process of 'gland cell type assembly': the stitching together of unique biosynthesis pathways; coordinated changes in secretory systems to enable efficient chemical release; and transcriptional deployment of these machineries into cells constituting the gland. This molecular evolutionary process influences what types of compound a given species is capable of secreting, and, consequently, the kinds of ecological interactions that species can display. Here, we discuss what is known about the evolutionary assembly of gland cell types and propose a framework for how it may happen. We posit the existence of 'terminal selector' transcription factors that program gland function via regulatory recruitment of biosynthetic enzymes and secretory proteins. We suggest ancestral enzymes are initially co-opted into the novel gland, fostering pleiotropic conflict that drives enzyme duplication. This process has yielded the observed pattern of modular, gland-specific biosynthesis pathways optimized for manufacturing specific secretions. We anticipate that single-cell technologies and gene editing methods applicable in diverse species will transform the study of animal chemical interactions, revealing how gland cell types are assembled and functionally configured at a molecular level.

Synergism of open chromatin regions involved in regulating genes in Bombyx mori

The dynamic variability of transcription factors (TFs) and their binding sites makes it challenging to conduct genome-wide transcription regulation research. The silkworm Bombyx mori, which produces silk, is one of the most valuable model insects in the order Lepidoptera. The "opening" and "closing" of chromatin in different silk yield strains is associated with changes in silk production, making this insect a good model for studying the transcriptional regulation of genes. However, few studies have examined the open chromatin regions (OCRs) of silkworms, and studying OCR synergism and their function in silk production remains challenging. Here, we performed formaldehyde-assisted isolation of regulatory elements (FAIRE) to isolate OCRs from the silk glands of fifth-instar larvae of the DaZao and D872 strains. In total, 128,908 high confidence OCRs were identified and approximately 80% of OCRs were located in non-coding regions. OCRs upregulated adjacent genes and showed signal-dependent vulnerability to single-nucleotide polymorphisms. Mid- and low-signal OCRs were more likely to have single-nucleotide polymorphisms (SNP). Further, OCRs interacted with each other within a distance of 5 kb. We named the OCR interaction complex as the "cluster of related regions" (COREs). The functions of the CORE and its harbored OCRs showed some differences. Additionally, COREs enriched many silk protein synthesis-associated genes, some of which were upregulated. This study identified numerous high confidence regulation sites and synergistic regulatory modes of OCRs that affect adjacent genes. These results provide insight into silkworm transcriptional regulation and improve our understanding of cis-element cooperation.

Genome-wide open chromatin regions and their effects on the regulation of silk protein genes in Bombyx mori

Nucleosome-depleted open chromatin regions (OCRs) often harbor transcription factor (TF) binding sites that are associated with active DNA regulatory elements. To investigate the regulation of silk-protein genes, DNA molecules isolated from the silk glands of third-day fifth-instar silkworm larvae and embryo-derived (BmE) cells were subjected to formal dehyde-assisted isolation of regulatory elements (FAIRE) and high-throughput sequencing. In total, 68,000 OCRs were identified, and a number of TF-binding motifs were predicted. In particular, OCRs located near silk-protein genes contained potential binding sites for functional TFs. Moreover, many TFs were found to bind to clusters of OCRs upstream of silk-protein genes, and to regulate the expression of these genes. The expression of silk protein genes may be related not only to regulating TFs (such as fkh, Bmdimm, and Bmsage), but also to developmental and hormone-induced TFs (such as zen, eve, Br, and eip74ef). Elucidation of genome-wide OCRs and their regulatory motifs in silk protein genes will provide valuable data and clues for characterizing the mechanisms of transcriptional control of silk protein genes.

Histomorphometric and transcriptomic features characterize silk glands' development during the molt to intermolt transition process in silkworm

The molt-intermolt cycle is an essential feature in holometabolous and hemimetabolous insects' development. In the silkworm, silk glands are under dramatic morphological and functional changes with fibroin genes' transcription being repeatedly turned off and on during the molt-intermolt cycles. However, the molecular mechanisms controlling it are still unknown. Here, silk gland's histomorphology and transcriptome analysis were used to characterize changes in its structure and gene expression patterns from molt to intermolt stages. By using section staining and transmission electron microscope, a renewable cell damage was detected in the silk gland at the molt stage, and an increased number of autophagosomes and lysosomes were found in silk gland cells' cytoplasm. Next, by using RNA sequencing, 54,578,413 reads were obtained, of which 85% were mapped to the silkworm reference genome. The expression level analysis of silk protein genes and silk gland transcription factors revealed that fibroin heavy chain, fibroin light chain, P25/fhx, sericin1, sericin3 and Dimm had consistent alteration trends in temporal expression. In addition, differentially expressed genes (DEGs) were identified, and most of the DEGs associated with ecdysone signal transduction, mRNA degradation, protein proteolysis, and autophagy were significantly down-regulated in the transition from molt to intermolt, suggesting that these pathways were activated for the silk gland renewal. These findings provide insights into the molecular mechanisms of silk gland development and silk protein genes transcriptional regulation during the molt to intermolt transition process.

Regulation of Silk Genes by Hox and Homeodomain Proteins in the Terminal Differentiated Silk Gland of the Silkworm Bombyx mori

The silk gland of the silkworm Bombyx mori is a long tubular organ that is divided into several subparts along its anteroposterior (AP) axis. As a trait of terminal differentiation of the silk gland, several silk protein genes are expressed with unique regional specificities. Most of the Hox and some of the homeobox genes are also expressed in the differentiated silk gland with regional specificities. The expression patterns of Hox genes in the silk gland roughly correspond to those in embryogenesis showing "colinearity". The central Hox class protein Antennapedia (Antp) directly regulates the expression of several middle silk gland-specific silk genes, whereas the Lin-1/Isl-1/Mec3 (LIM)-homeodomain transcriptional factor Arrowhead (Awh) regulates the expression of posterior silk gland-specific genes for silk fiber proteins. We summarize our results and discuss the usefulness of the silk gland of Bombyx mori for analyzing the function of Hox genes. Further analyses of the regulatory mechanisms underlying the region-specific expression of silk genes will provide novel insights into the molecular bases for target-gene selection and regulation by Hox and homeodomain proteins.