Genome-wide identification of polycomb target genes reveals a functional association of Pho with Scm in Bombyx mori

Unique territorial and sub-chromosomal organization revealed in the holocentric moth Bombyx mori

The hallmarks of chromosome organization in multicellular eukaryotes are chromosome territories (CT), chromatin compartments, and different types of domains, including topologically associated domains (TADs). Yet, most of these concepts derive from analyses of organisms with monocentric chromosomes. Here we describe the 3D genome architecture of an organism with holocentric chromosomes, the silkworm Bombyx mori . At the genome-wide scale, B. mori chromosomes form highly separated territories and lack substantial trans contacts. As described in other eukaryotes, B. mori chromosomes segregate into an active A and an inactive B compartment. Remarkably, we also identify a third compartment, Secluded "S", with a unique contact pattern. Compartment S shows strong enrichment of short-range contacts and depletion of long-range contacts. It hosts a unique combination of genetic and epigenetic features, localizes at the periphery of CTs and shows developmental plasticity. Biophysical modeling shows that formation of such secluded domains requires a new mechanism - a high density of extruded loops within them along with low level of extrusion and compartmentalization of A and B. Together with other evidence of loop extrusion in interphase, this suggests SMC-mediated loop extrusion in this insect. Overall, our analyses highlight the evolutionary plasticity of 3D genome organization driven by a new combination of known processes.

Adaptation of Helicoverpa armigera to Soybean Peptidase Inhibitors Is Associated with the Transgenerational Upregulation of Serine Peptidases

Molecular phenotypes induced by environmental stimuli can be transmitted to offspring through epigenetic inheritance. Using transcriptome profiling, we show that the adaptation of Helicoverpa armigera larvae to soybean peptidase inhibitors (SPIs) is associated with large-scale gene expression changes including the upregulation of genes encoding serine peptidases in the digestive system. Furthermore, approximately 60% of the gene expression changes induced by SPIs persisted in the next generation of larvae fed on SPI-free diets including genes encoding regulatory, oxidoreductase, and protease functions. To investigate the role of epigenetic mechanisms in regulating SPI adaptation, the methylome of the digestive system of first-generation larvae (fed on a diet with and without SPIs) and of the progeny of larvae exposed to SPIs were characterized. A comparative analysis between RNA-seq and Methyl-seq data did not show a direct relationship between differentially methylated and differentially expressed genes, while trypsin and chymotrypsin genes were unmethylated in all treatments. Rather, DNA methylation potential epialleles were associated with transcriptional and translational controls; these may play a regulatory role in the adaptation of H. armigera to SPIs. Altogether, our findings provided insight into the mechanisms of insect adaptation to plant antiherbivore defense proteins and illustrated how large-scale transcriptional reprograming of insect genes can be transmitted across generations.

SUMOylation of Translationally Regulated Tumor Protein Modulates Its Immune Function

Translationally controlled tumor protein (TCTP) is a highly conserved protein possessing numerous biological functions and molecular interactions, ranging from cell growth to immune responses. However, the molecular mechanism by which TCTP regulates immune function is largely unknown. Here, we found that knockdown of Bombyx mori translationally controlled tumor protein (BmTCTP) led to the increased susceptibility of silkworm cells to virus infection, whereas overexpression of BmTCTP significantly decreased the virus replication. We further demonstrated that BmTCTP could be modified by SUMOylation molecular BmSMT3 at the lysine 164 via the conjugating enzyme BmUBC9, and the stable SUMOylation of BmTCTP by expressing BmTCTP-BmSMT3 fusion protein exhibited strong antiviral activity, which confirmed that the SUMOylation of BmTCTP would contribute to its immune responses. Further work indicated that BmTCTP is able to physically interact with interleukin enhancer binding factor (ILF), one immune molecular, involved in antivirus, and also induce the expression of BmILF in response to virus infection, which in turn enhanced antiviral activity of BmTCTP. Altogether, our present study has provided a novel insight into defending against virus via BmTCTP SUMOylation signaling pathway and interacting with key immune molecular in silkworm.

Efficient Delivery of dsRNA and DNA in Cultured Silkworm Cells for Gene Function Analysis Using PAMAM Dendrimers System

: Polyamidoamine (PAMAM) dendrimers are emerging as intriguing nanovectors for nucleic acid delivery because of their unique well-defined architecture and high binding capacity, which have been broadly applied in DNA- and RNA-based therapeutics. The low-cost and high-efficiency of PAMAM dendrimers relative to traditional liposomal transfection reagents also promote their application in gene function analysis. In this study, we first investigated the potential use of a PAMAM system in the silkworm model insect. We determined the binding property of G5-PAMAM using dsRNA and DNA in vitro, and substantially achieved the delivery of dsRNA and DNA from culture medium to both silkworm BmN and BmE cells, thus leading to efficient knockdown and expression of target genes. Under treatments with different concentrations of G5-PAMAM, we evaluated its cellular cytotoxicity on silkworm cells, and the results show that G5-PAMAM had no obvious toxicity to cells. The presence of serum in the culture medium did not affect the delivery performance of DNA and dsRNA by G5-PAMAM, revealing its convenient use for various purposes. In conclusion, our data demonstrate that the PAMAM system provides a promising strategy for delivering dsRNA and DNA in cultured silkworm cells and promote its further application in individuals.

Transcriptome-wide analysis of N6-methyladenosine uncovers its regulatory role in gene expression in the lepidopteran Bombyx mori

N6-methyladenosine (m6A) methylation is the most common form of RNA modification in eukaryotes and has been shown to act as an important epitranscriptomic marker, which can regulate gene expression, messenger RNA stability, alternative splicing and translation efficiency. Recent progress on determining insect m6A function has been limited to the dipteran Drosophila melanogaster, in which m6A is involved in neuronal functions and sex determination. Methylation and function of m6A in other insects, however, remain unknown. Here, we investigated a transcriptome-wide profile of m6A in the lepidopteran Bombyx mori and identified the methyltransferase subunits B. mori methyltransferase-like 3 (BmMETTL3) and BmMETTL14 in the m6A methylation pathway. Strikingly, loss of BmMETTL3 and BmMETTL14 in cultured B. mori cells led to arrest of cell cycle progression and caused deficiency of chromosome alignment and segregation. Specifically, we identified 2853 m6A peaks representing transcripts of 2043 genes, and the target genes with m6A methylation were shown to be involved in gene expression and translation. It was interesting that we found that the highly expressed genes tended to be methylated by m6A, and comparative analysis of RNA m6A and DNA N6-methyladenine (6mA) revealed two distinct regulatory mechanisms for gene expression. Overall, our work suggests RNA m6A and DNA 6mA play important roles in RNA and DNA epigenetic regulation in B. mori.

A functional polypeptide N-acetylgalactosaminyltransferase (PGANT) initiates O-glycosylation in cultured silkworm BmN4 cells

Mucin-type O-glycosylation is initiated by UDP-GalNAc:polypeptide N-acetylgalactosaminyltransferases (ppGalNAc-Ts or PGANTs), attaching GalNAc to serine or threonine residue of a protein substrate. In the insect model from Lepidoptera, silkworm (Bombyx mori), however, O-glycosylation pathway is totally unexplored and remains largely unknown. In this study, as the first report regarding protein O-glycosylation analysis in silkworms, we verified the O-glycan profile that a common core 1 Gal (β1-3) GalNAc disaccharide branch without terminally sialylated structure is mainly formed for a baculovirus-produced human proteoglycan 4 (PRG4) protein. Intriguingly, functional screenings in cultured silkworm BmN4 cells for nine Bmpgants reveal that Bmpgant2 is the solo functional BmPGANT for PRG4, implying that Bmpgants may have unique cell/tissue or protein substrate preferences. Furthermore, a recombinant BmPGANT2 protein was successfully purified from silkworm-BEVS and exhibited a high ability to transfer GalNAc for both peptide and protein substrates. Taken together, the present results clarified the functional BmPGANT2 in cultured silkworm cells, providing crucial fundamental insights for future studies dissecting the detailed silkworm O-glycosylation pathways and productions of glycoproteins with O-glycans.

Multiprotein bridging factor 2 regulates the expression of the fibroin heavy chain gene by interacting with Bmdimmed in the silkworm Bombyx mori

Multiprotein bridging factor 2 (MBF2) was first isolated from the posterior silk gland of Bombyx mori. However, its function in B. mori is still unknown. Herein, MBF2 transcripts were detected mainly in the posterior silk gland and Malpighian tubules of B. mori larvae via a quantitative PCR analysis. An analysis of temporal expression patterns showed that the expression pattern of MBF2 was the opposite of that of the fibroin heavy chain (fibH) gene, as its expression was high during the fourth-instar moulting stage, decreased gradually during the fifth-instar feeding stage and disappeared at the end of the fifth-instar phase. Furthermore, bimolecular fluorescent complementation and Far-Western blot assays showed that MBF2 interacted with the basic helix-loop-helix transcription factor Bmdimmed. Dual luciferase reporter assays showed that MBF2 down-regulated the promoter activity of fibH and inhibited the effect of Bmdimmed (Bmdimm) on fibH expression. MBF2 expression was induced in silk glands after treatment with 20-hydroxyecdysone in vivo and in vitro. These findings suggest that MBF2 is a transcriptional repressor that is involved in controlling the regulation of the fibH gene in the posterior silk gland by interacting with Bmdimm.

Ecdysone receptor (EcR) is involved in the transcription of cell cycle genes in the silkworm

EcR (ecdysone receptor)-mediated ecdysone signaling pathway contributes to regulate the transcription of genes involved in various processes during insect development. In this work, we detected the expression of EcR gene in silkworm ovary-derived BmN4 cells and found that EcR RNAi result in an alteration of cell shape, indicating that EcR may orchestrate cell cycle progression. EcR RNAi and EcR overexpression analysis revealed that in the cultured BmN4 cells, EcR respectively promoted and suppressed the transcription of E2F-1 and CycE, two genes controlling cell cycle progression. Further examination demonstrated that ecdysone application in BmN4 cells not only changed the transcription of these two cell cycle genes like that under EcR overexpression, but also induced cell cycle arrest at G2/M phase. In vivo analysis confirmed that E2F-1 expression was elevated in silk gland of silkworm larvae after ecdysone application, which is same as its response to ecdysone in BmN4 cells. However, ecdysone also promotes CycE transcription in silk gland, and this is converse with the observation in BmN4 cells. These results provide new insights into understanding the roles of EcR-mediated ecdysone signaling in the regulation of cell cycle.

Dynamics of polycomb proteins-mediated histone modifications during UV irradiation-induced DNA damage

Polycomb group (PcG) complexes are known to be chromatin modifiers and transcriptional repressors. In this work, we reported that the histone-modifying PcG complexes are able to participate in the repair process of ultraviolet (UV)-induced DNA lesions in the silkworm, Bombyx mori. The silkworm cells with depletion of PcG genes showed hypersensitive to UV-C irradiation and increased inhibition of cell proliferation. Interestingly, an SQ site in the silkworm-human chimeric H2A protein synthesized here was phosphorylated rapidly upon UV-C exposure, which could be used as a marker for monitoring the response to DNA damage in silkworm cells. Under these UV-C irradiated conditions, we found that PRC1-mediated ubiquitylation of H2AX, but not of H2AZ, were decreased and this deubiquitylation was independent of its phosphorylation event. In contrast, UV-C irradiation induced the increase of trimethylation of lysine 27 on histone H3 (H3K27me3), a mark of transcriptionally silent chromatin catalyzed by another PcG subcomplex, PRC2. Collectively, we provided the first evidence on chromatin remodeling in response to UV-C lesion in silkworm and revealed another layer role for PcG complexes-mediated histone modifications in contributing to creating an open chromatin structure for the efficient repair of DNA damages.

Nuclear import of transcription factor BR-C is mediated by its interaction with RACK1

The transcription factor Broad Complex (BR-C) is an early ecdysone response gene in insects and contains two types of domains: two zinc finger domains for the activation of gene transcription and a Bric-a-brac/Tramtrack/Broad complex (BTB) domain for protein-protein interaction. Although the mechanism of zinc finger-mediated gene transcription is well studied, the partners interacting with the BTB domain of BR-C has not been elucidated until now. Here, we performed a yeast two-hybrid screen using the BTB domain of silkworm BR-C as bait and identified the receptor for activated C-kinase 1 (RACK1), a scaffolding/anchoring protein, as the novel partner capable of interacting with BR-C. The interaction between BR-C and RACK1 was further confirmed by far-western blotting and pull-down assays. Importantly, the disruption of this interaction, via RNAi against the endogenous RACK1 gene or deletion of the BTB domain, abolished the nuclear import of BR-C in BmN4 cells. In addition, RNAi against the endogenous PKC gene as well as phosphorylation-deficient mutation of the predicted PKC phosphorylation sites at either Ser373 or Thr406 in BR-C phenocopied RACK1 RNAi and altered the nuclear localization of BR-C. However, when BTB domain was deleted, phosphorylation mimics of either Ser373 or Thr406 had no effect on the nuclear import of BR-C. Moreover, mutating the PKC phosphorylation sites at Ser373 and Thr406 or deleting the BTB domain significantly decreased the transcriptional activation of a BR-C target gene. Given that RACK1 is necessary for recruiting PKC to close and phosphorylate target proteins, we suggest that the PKC-mediated phosphorylation and nuclear import of BR-C is determined by its interaction with RACK1. This novel finding will be helpful for further deciphering the mechanism underlying the role of BR-C proteins during insect development.

A conserved SUMOylation signaling for cell cycle control in a holocentric species Bombyx mori

SUMOylation is an essential post-translational modification that regulates a variety of cellular processes including cell cycle progression. Although the SUMOylation pathway has been identified and investigated in many eukaryotes, the mechanisms of SUMOylation in regulating the functions of various substrates are still poorly understood. Here, we utilized a model species, the silkworm Bombyx mori that possesses holocentric chromosomes, to exploit the role of the SUMOylation system in cell cycle regulation. We identified all the components that are involved in the SUMOylation pathway in the silkworm genome. Our data revealed a cell cycle-dependent transcription of the SUMOylation genes, localization of the SUMOylation proteins, and abundance of the SUMOylation substrates in cultured silkworm cells. Importantly, the proliferation of the silkworm cells was strikingly inhibited by interference with SUMOylation genes expression, possibly due to an arrest of the SUMOylation-deficient cells at the G2/M phase. Furthermore, disruption of the SUMOylation genes induced the defects of holocentric chromosome congression and segregation during mitosis, which was consistent with high expressions of the SUMOylation genes and high enrichments of global SUMOylation at this stage, suggesting that the SUMOylation system in silkworm is essential for cell cycle regulation, with one particular role in mitosis.

RNAi silencing of the SoxE gene suppresses cell proliferation in silkworm BmN4 cells

The transcription factor SoxE is mainly expressed in the gonad and involved in the regulation of gonad development and sex determination in animals. Here, we used the silkworm ovary-derived BmN4-SID1 cell line to survey the roles of the silkworm SoxE protein (BmSoxE) and predict its candidate binding targets. RNAi-mediated silencing of BmSoxE expression suppressed cell proliferation in BmN4-SID1 cells. A further cell cycle analysis revealed that this inhibition of cell proliferation was largely due to cell cycle arrest in G1 phase when BmSoxE expression was blocked in BmN4-SID1 cells. Genome-wide microarray expression analyses demonstrated that the expression levels of a set of genes were significantly altered following BmSoxE RNAi. More than half of these genes contained conserved binding sites for HMG box domain of the Sox proteins and were predicted to be candidate binding targets for BmSoxE. Importantly, some of the candidate targets may be associated with the effect of BmSoxE on cell proliferation. Several candidate target genes showed gonad-specific expression in silkworm larvae. Taken together, these data demonstrate that BmSoxE is required for cell proliferation in silkworm BmN4-SID1 cells and provide valuable information for further investigations of the molecular control exerted by the BmSoxE protein over cell proliferation and gonad development in the silkworm.

Roles of Piwi proteins in transcriptional regulation mediated by HP1s in cultured silkworm cells

Piwi proteins are part of a superfamily of Argonaute proteins, which are one of the core components of the RNA silencing pathway in many eukaryotes. Piwi proteins are thought to repress the transposon expression both transcriptionally and post-transcriptionally. Recently, Drosophila melanogaster Piwi was recently reported to associate with chromatin and to interact directly with the Heterochromatin Protein 1 (HP1a). However, similar interactions have not been reported in other higher eukaryotes. Here we show that silkworm Piwi proteins interact with HP1s in the nucleus. The silkworm, Bombyx mori, has two Piwi proteins, Ago3 and Siwi, and two typical HP1 proteins, HP1a and HP1b. We found that HP1a plays an important role in the interaction between Ago3/Siwi and HP1b in the ovary-derived BmN4 cell line. We also found that Ago3/Siwi regulates the transcription in an HP1-dependent manner. These results suggest that silkworm Piwi proteins function as a chromatin regulator in collaboration with HP1a and HP1b.

Characterization of Tudor-sn-containing granules in the silkworm, Bombyx mori

The Tudor-sn protein, which contains four staphylococcal nuclease domains and a Tudor domain, is a ubiquitous protein found in almost all organisms. It has been reported that Tudor-sn in mammals participates in various cellular pathways involved in gene regulation, cell growth, and development. In insects, we have previously identified a Tudor-sn ortholog in the silkworm, Bombyx mori, and detected its interactions between with Argonaute proteins. The role of Tudor-sn in silkworm, however, still remains largely unknown. In this study, we demonstrated that silkworm Tudor-sn is a stress granule (SG) protein, and determined its interactions with other SG proteins using Bimolecular Fluorescence Complementation assay and Insect Two-Hybrid method. Depletions of Argonaute proteins and SG-marker protein Tia1 by RNAi impaired the involvement of Tudor-sn in the SG formation. Protein domain deletion analysis of Tudor-sn demonstrated that SN2 is the key domain required for the aggregation of Tudor-sn in SGs.

Cell cycle-dependent recruitment of polycomb proteins to the ASNS promoter counteracts C/ebp-mediated transcriptional activation in Bombyx mori

Epigenetic modifiers and transcription factors contribute to developmentally programmed gene expression. Here, we establish a functional link between epigenetic regulation by Polycomb group (PcG) proteins and transcriptional regulation by C/ebp that orchestrates the correct expression of Bombyx mori asparagine synthetase (BmASNS), a gene involved in the biosynthesis of asparagine. We show that the cis-regulatory elements of YY1-binding motifs and the CpG island present on the BmASNS promoter are required for the recruitment of PcG proteins and the subsequent deposition of the epigenetic repression mark H3K27me3. RNAi-mediated knockdown of PcG genes leads to derepression of the BmASNS gene via the recruitment of activators, including BmC/ebp, to the promoter. Intriguingly, we find that PcG proteins and BmC/ebp can dynamically modulate the transcriptional output of the BmASNS target in a cell cycle-dependent manner. It will be essential to suppress BmASNS expression by PcG proteins at the G2/M phase of the cell cycle in the presence of BmC/ebp activator. Thus, our results provide a novel insight into the molecular mechanism underlying the recruitment and regulation of the PcG system at a discrete gene locus in Bombyx mori.