A Genome-Wide Survey for Host Response of Silkworm, Bombyx mori during Pathogen Bacillus bombyseptieus Infection

Host-pathogen interactions are complex relationships, and a central challenge is to reveal the interactions between pathogens and their hosts. Bacillus bombysepticus (Bb) which can produces spores and parasporal crystals was firstly separated from the corpses of the infected silkworms (Bombyx mori). Bb naturally infects the silkworm can cause an acute fuliginosa septicaemia and kill the silkworm larvae generally within one day in the hot and humid season. Bb pathogen of the silkworm can be used for investigating the host responses after the infection. Gene expression profiling during four time-points of silkworm whole larvae after Bb infection was performed to gain insight into the mechanism of Bb-associated host whole body effect. Genome-wide survey of the host genes demonstrated many genes and pathways modulated after the infection. GO analysis of the induced genes indicated that their functions could be divided into 14 categories. KEGG pathway analysis identified that six types of basal metabolic pathway were regulated, including genetic information processing and transcription, carbohydrate metabolism, amino acid and nitrogen metabolism, nucleotide metabolism, metabolism of cofactors and vitamins, and xenobiotic biodegradation and metabolism. Similar to Bacillus thuringiensis (Bt), Bb can also induce a silkworm poisoning-related response. In this process, genes encoding midgut peritrophic membrane proteins, aminopeptidase N receptors and sodium/calcium exchange protein showed modulation. For the first time, we found that Bb induced a lot of genes involved in juvenile hormone synthesis and metabolism pathway upregulated. Bb also triggered the host immune responses, including cellular immune response and serine protease cascade melanization response. Real time PCR analysis showed that Bb can induce the silkworm systemic immune response, mainly by the Toll pathway. Anti-microorganism peptides (AMPs), including of Attacin, Lebocin, Enbocin, Gloverin and Moricin families, were upregulated at 24 hours post the infection

Microarray-based gene expression profiles in multiple tissues of the domesticated silkworm, Bombyx mori

Using a genome-wide oligonucleotide microarray, gene expression was surveyed in multiple silkworm tissues on day 3 of the fifth instar, providing a new resource for annotating the silkworm genome

Genome-wide annotation and comparative analysis of cuticular protein genes in the noctuid pest \u3cem\u3eSpodoptera litura\u3c/em\u3e

Insect cuticle is considered an adaptable and versatile building material with roles in the construction and function of exoskeleton. Its physical properties are varied, as the biological requirements differ among diverse structures and change during the life cycle of the insect. Although the bulk of cuticle consists basically of cuticular proteins (CPs) associated with chitin, the degree of cuticular sclerotization is an important factor in determining its physical properties. Spodoptera litura, the tobacco cutworm, is an important agricultural pest in Asia. Compared to the domestic silkworm, Bombyx mori, another lepidopteran whose CP genes have been well annotated, S. litura has a shorter life cycle, hides in soil during daytime beginning in the 5th instar and is exposed to soil in the pupal stage without the protection of a cocoon. In order to understand how the CP genes may have been adapted to support the characteristic life style of S. litura, we searched its genome and found 287 putative cuticular proteins that can be classified into 9 CP families (CPR with three groups (RR-1, RR-2, RR-3), CPAP1, CPAP3, CPF, CPFL, CPT, CPG, CPCFC and CPLCA), and a collection of unclassified CPs named CPH. There were also 112 cuticular proteins enriched in Histidine residues with content varying from 6% to 30%, comprising many more His-rich cuticular proteins than B. mori. A phylogenetic analysis between S. litura, M. sexta and B. mori uncovered large expansions of RR-1 and RR-2 CPs, forming large gene clusters in different regions of S. liturachromosome 9. We used RNA-seq analysis to document the expression profiles of CPs in different developmental stages and tissues of S. litura. The comparative genomic analysis of CPs between S. litura and B. moriintegrated with the unique behavior and life cycle of the two species offers new insights into their contrasting ecological adaptations

Expression map of a complete set of gustatory receptor genes in chemosensory organs of \u3cem\u3eBombyx mori\u3c/em\u3e

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

SilkDB v2.0: a platform for silkworm (Bombyx mori ) genome biology

The SilkDB is an open-access database for genome biology of the silkworm (Bombyx mori). Since the draft sequence was completed and the SilkDB was first released 5 years ago, we have collaborated with other groups to make much remarkable progress on silkworm genome research, such as the completion of a new high-quality assembly of the silkworm genome sequence as well as the construction of a genome-wide microarray to survey gene expression profiles. To accommodate these new genomic data and house more comprehensive genomic information, we have reconstructed SilkDB database with new web interfaces. In the new version (v2.0) of SilkDB, we updated the genomic data, including genome assembly, gene annotation, chromosomal mapping, orthologous relationship and experiment data, such as microarray expression data, Expressed Sequence Tags (ESTs) and corresponding references. Several new tools, including SilkMap, Silkworm Chromosome Browser (SCB) and BmArray, are developed to access silkworm genomic data conveniently. SilkDB is publicly available at the new URL of http://www.silkdb.org

Functional Divergence among Silkworm Antimicrobial Peptide Paralogs by the Activities of Recombinant Proteins and the Induced Expression Profiles

Antimicrobial peptides are small-molecule proteins that are usually encoded by multiple-gene families. They play crucial roles in the innate immune response, but reports on the functional divergence of antimicrobial peptide gene families are rare. In this study, 14 paralogs of antimicrobial peptides belonging to cecropin, moricin and gloverin families were recombinantly expressed in pET expression systems. By antimicrobial activity tests, peptides representing paralogs in the same family of cecropin and moricin families, displayed remarkable differences against 10 tested bacteria. The evolutionary rates were relatively fast in the two families, which presented obvious functional divergence among paralogs of each family. Four peptides of gloverin family had similar antimicrobial spectrum and activity against tested bacteria. The gloverin family showed similar antimicrobial function and slow evolutionary rates. By induced transcriptional activity, genes encoding active antimicrobial peptides were upregulated at obviously different levels when silkworm pupae were infected by three types of microbes. Association analysis of antimicrobial activities and induced transcriptional activities indicated that the antimicrobial activities might be positively correlated with induced transcriptional activities in the cecropin and moricin families. These results suggest that representative BmcecB6, BmcecD and Bmmor as the major effector genes have broad antimicrobial spectrum, strong antimicrobial activity and high microbe-induced expression among each family and maybe play crucial roles in eliminating microbial infection

Gene Expression Pattern and Regulatory Network of α-Toxin Treatment in Bombyx mori

Bacillus bombyseptieus is a pathogen of Bombyx mori; it can cause bacterial septicemia in silkworm. One of the components of the parasporal crystal toxin of B. bombyseptieus, α-toxin, plays an important role in the process of infection in silkworm. In this study, we investigated the immune response of silkworm induced by α-toxin by using RNA-seq. We compared the changes in gene expression in the midgut, fatbody, and hemocytes of silkworm and in the B. mori embryonic cell line (BmE) after treatment with α-toxin and identified 952 differentially expressed genes and 353 differentially expressed long noncoding RNAs (lncRNAs). These regulated genes in different tissues were found to be enriched in different pathways. The upregulated genes in the midgut were mainly involved in peptidoglycan catabolic process and tyrosine kinase signaling pathway, whereas the downregulated genes were mainly involved in chitin metabolic pathways. The upregulated genes in fatbody were also involved in peptidoglycan catabolic process, but they were for a different peptidoglycan subtype. Further, genes encoding cecropins were enriched in the fatbody. The downregulated genes were mainly involved in the metabolic pathways of fundamental substances such as cellular protein metabolic process and nucleobase-containing compound metabolic process. These results suggest that α-toxin can induce various immune responses in silkworm, and further studies are warranted to understand the mechanism of α-toxin action in silkworm. Further, lncRNAs and differentially expressed genes were correlated using coexpression network analysis. Our findings revealed potential candidate genes and lncRNAs that might play important physiological functions in the immune response to α-toxins in silkworm

Effects of Osiris9a on Silk Properties in Bombyx mori Determined by Transgenic Overexpression

Osiris is an insect-specific gene family with multiple biological roles in development, phenotypic polymorphism, and protection. In the silkworm, we have previously identified twenty-five Osiris genes with high evolutionary conservation and remarkable synteny among several insects. Bombxy mori Osiris9a (BmOsi9a) is expressed only in the silk gland, particularly in the middle silk gland (MSG). However, the biological function of BmOsi9a is still unknown. In this study, we overexpressed BmOsi9a in the silk gland by germline transgene expression. BmOsi9a was overexpressed not only in the MSG but also in the posterior silk gland (PSG). Interestingly, BmOsi9a could be secreted into the lumen in the MSG but not in the PSG. In the silk fiber, overexpressed BmOsi9a interacted with Sericin1 in the MSG, as confirmed by a co-immunoprecipitation assay. The overexpression of BmOsi9a altered the secondary structure and crystallinity of the silk fiber, thereby changing the mechanical properties. These results provide insight into the mechanisms underlying silk proteins secretion and silk fiber formation

Circular transcriptome sequencing of the middle silk gland and posterior silk gland in the Bombyx mori

Circular transcriptome sequencing of the middle silk gland (MSG) and posterior silk gland (PSG) in the Bombyx mori are presented. The middle silk gland and posterior silk gland were collected from the third day of fifth-instar B.mori larvae. The circular RNAs enriched by using RNase R to degrade the linear RNA molecules, and circular RNA sequencing (circRNA-seq) was performed using an Illumina Hiseq. 2500 sequencing platform. Samples are described in the SRA portal (SRP100385) and FASTQ files have been deposited in Sequence Read Archive (accession numbers: SRX2577343 and SRX2577342). The interpretation of these data is presented in the following research article: âIdentification of circular RNA in the Bombyx mori silk glandâ [1] (Gan et al., 2017). Keywords: CircRNA-seq, Bombyx mori, MSG, PS