Search

Bombyx mori is an economically important holometabolous lepidopteran insect. In B. mori endogenous noncoding RNAs such as microRNAs (miRNAs) and Piwi-interacting RNAs play crucial biological functions in metamorphosis and sex determination. In addition, circular RNAs (circRNAs) have been recently identified as noncoding RNAs in most common model organisms and show potential as gene regulators. However, to date, there have been few studies on the circRNAs present in the B. mori genome conducted to date. Here, we identified 3916 circRNAs by deep circular transcriptome sequencing using the silk gland of B. mori. 3155 circRNAs were found to be derived from 1727 parental genes. The circRNAs displayed tissue-specific expression between the middle silk gland (MSG) and posterior silk gland (PSG), with 2532 and 880 being upregulated circRNAs in the MSG and PSG, respectively. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses showed that the parental genes from the MSG and PSG were generally annotated to similar categories and pathways. The interaction network of circRNAs and miRNAs showed that circRNAs might act as miRNA sponges or interact with miRNAs in some other way. Overall, the results revealed the complicated patterns of circRNAs in the B. mori silk gland providing a new angle from which to explore the mechanisms of complex gene regulation and efficient silk protein synthesis.

We established a large-scale separation and purification platform to obtain kilogram amounts of natural compounds from the extraction of the fruiting bodies of C. militaris. Seven monomeric compounds, N6-(2-hydroxyethyl) adenosine (HEA), ergosterol (E), ergosta-7,22-diene-3,5,6-triol (EI), 5α,8α-epidioxy-(22E,24R)-ergosta-6,22-dien-3β-ol (ED),ergosta-7,22-dien-3β,5α-dihydroxy-6-one (EO), (20S,22E,24R)-Eegosta-7,22-dien-3β,5α,6β,9α-tetraol (ET), and (24S)-5,22-stigmastadien-3β-ol (SE), were harvested using different solvents, and the structure of each compound was identified. The activities and functions of the isolated compounds were tested by label-free, real-time cell analysis methods at the cellular level, and their antitumor effects were verified using mouse models of Lewis and H22 tumors. The anti-insomnia effect of HEA was tested in an anti-insomnia mouse model. The interactions between E and 8 A549 cell proteins were determined. The biosynthetic pathways of HEA and E, which possess pharmacologically active monomers, were determined. This platform can provide a theoretical basis for the further development and discovery of novel natural medicines.

Tarantula venoms provide a model system for studying toxin selectivity, structure-activity relationships and molecular evolution of peptide toxins. Previous studies have identified a large number of peptide toxins in the venom of the Chinese bird spider Haplopelma hainanum, generally regarded as a highly venomous spider. However, the lack of available RNA-seq transcriptomic and genomic data is an obstacle to understanding its venom at the molecular level. In this study, we investigated the venom gland transcriptome of H. hainanum by RNA-seq, in the absence of an available genomic sequence. We identified 201 potential toxins among 57 181 de novo assembled transcripts, including knottins, Kunitz-type toxins, enzymes and other proteins. We systematically identified most of the knottins and Kunitz-type toxins, some of which showed strongly biased expression in the venom gland, including members of the huwentoxin-1, huwentoxin-2 and magi-1 families. We also discovered several novel potential toxins. These data demonstrate the high molecular and structural diversity in the venom toxins of H. hainanum. This study offers a useful strategy for exploring the complex components of spider venoms.

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, PSG

The complete genome of Cordyceps militaris was sequenced using single-molecule real-time (SMRT) sequencing technology at a coverage over 300×. The genome size was 32.57 Mb, and 14 contigs ranging from 0.35 to 4.58 Mb with an N50 of 2.86 Mb were assembled, including 4 contigs with telomeric sequences on both ends and an additional 8 contigs with telomeric sequences on either the 5′ or 3′ end. A methylome database of the genome was constructed using SMRT and m4C and m6A methylated nucleotides, and many unknown modification types were identified. The major m6A methylation motif is GA and GGAG, and the major m4C methylation motif is GC or CG/GC. In the C. militaris genome DNA, there were four types of methylated nucleotides that we confirmed using high-resolution LCMS-IT-TOF. Using PacBio Iso-Seq, a total of 31,133 complete cDNA sequences were obtained in the fruiting body. The conserved domains of the nontranscribed regions of the genome include TATA boxes, which are the initial regions of genome replication. There were 406 structural variants between the HN and CM01 strains, and there were 1,114 structural variants between the HN and ATCC strains.

Beauveria bassiana, Cordyceps militaris and Ophiocordyceps sinensis (Ascomycotina) are traditional Chinese medicines. Here, mitogenomes of these three Ascomycotina fungi were sequenced and de-novo assembled using single-molecule real-time sequencing. The results showed that their complete mitogenomes were 31,258, 31,854 and 157,584 bp, respectively, with sequencing depth approximately 278,760×, 326,283× and 69,385×. Types of repeat sequences were mainly (AA)n, (AAT)n, (TA)n and (TATT)n. DNA methylation motifs were revealed in DNA modifications of these three fungi. We discovered new models of RNA editing through analysis of transcriptomes from B. bassiana and C. militaris. These data lay a solid foundation for further genetic and biological studies about these three fungi, especially for elucidating the mitogenome evolution and exploring the regulatory mechanism of adapting environment.

Spodoptera litura is a polyphagous pest and can feed on more than 100 species of plants, causing great damage to agricultural production. The SNP results showed that there were gene exchanges between different regions. To explore the variations of larger segments in S. litura genome, we used genome resequencing samples from 14 regions of China, India, and Japan to study the copy number variations (CNVs). We identified 3976 CNV events and 1581 unique copy number variation regions (CNVRs) occupying the 108.5 Mb genome of S. litura. A total of 5527 genes that overlapped with CNVRs were detected. Selection signal analysis identified 19 shared CNVRs and 105 group-specific CNVRs, whose related genes were involved in various biological processes in S. litura. We constructed the first CNVs map in S. litura genome, and our findings will be valuable for understanding the genomic variations and population differences of S. litura.

Although the draft genome sequence of silkworm is available for a decade, its genetic variations, especially structural variations, are far from well explored. In this study, we identified 1,298,659 SNPs and 9,731 indels, of which 32 % of SNPs and 92.2 % of indels were novel compared to previous silkworm re-sequencing analysis. In addition, we applied a read depth-based approach to investigate copy number variations among 21 silkworm strains at genome-wide level. This effort resulted in 562 duplicated and 41 deleted CNV regions, and among them 442 CNV were newly identified. Functional annotation of genes affected by these genetic variations reveal that these genes include a wide spectrum of molecular functions, such as immunity and drug detoxification, which are important for the adaptive evolution of silkworms. We further validated the predicted CNV regions using q-PCR. 94.7 % (36/38) of the selected regions show divergent copy numbers compared to a single-copy gene OR2. In addition, potential presence/absence variations are also observed in our study: 11 genes are present in the reference genome, but absent in other strains. Overall, we draw an integrative map of silkworm genetic variation at genome-wide level. The identification of genetic variations in this study improves our understanding that these variants play important roles in shaping phenotypic variations between wild and domesticated silkworms.

Genome editing is one of the most powerful tools for revealing gene function and improving crop plants. Recently, RNA-guided genome editing using the type II clustered regularly interspaced short palindromic repeats (CRISPR)-associated protein (Cas) system has been used as a powerful and efficient tool for genome editing in various organisms. Here, we report genome editing in tobacco (Nicotiana tabacum) mediated by the CRISPR/Cas9 system. Two genes, NtPDS and NtPDR6, were used for targeted mutagenesis. First, we examined the transient genome editing activity of this system in tobacco protoplasts, insertion and deletion (indel) mutations were observed with frequencies of 16.2-20.3% after transfecting guide RNA (gRNA) and the nuclease Cas9 in tobacco protoplasts. The two genes were also mutated using multiplexing gRNA at a time. Additionally, targeted deletions and inversions of a 1.8-kb fragment between two target sites in the NtPDS locus were demonstrated, while indel mutations were also detected at both the sites. Second, we obtained transgenic tobacco plants with NtPDS and NtPDR6 mutations induced by Cas9/gRNA. The mutation percentage was 81.8% for NtPDS gRNA4 and 87.5% for NtPDR6 gRNA2. Obvious phenotypes were observed, etiolated leaves for the psd mutant and more branches for the pdr6 mutant, indicating that highly efficient biallelic mutations occurred in both transgenic lines. No significant off-target mutations were obtained. Our results show that the CRISPR/Cas9 system is a useful tool for targeted mutagenesis of the tobacco genome.

The tobacco cutworm, Spodoptera litura, is among the most widespread and destructive agricultural pests, feeding on over 100 crops throughout tropical and subtropical Asia. By genome sequencing, physical mapping and transcriptome analysis, we found that the gene families encoding receptors for bitter or toxic substances and detoxification enzymes, such as cytochrome P450, carboxylesterase and glutathione-S-transferase, were massively expanded in this polyphagous species, enabling its extraordinary ability to detect and detoxify many plant secondary compounds. Larval exposure to insecticidal toxins induced expression of detoxification genes, and knockdown of representative genes using short interfering RNA (siRNA) reduced larval survival, consistent with their contribution to the insect’s natural pesticide tolerance. A population genetics study indicated that this species expanded throughout southeast Asia by migrating along a South India–South China–Japan axis, adapting to wide-ranging ecological conditions with diverse host plants and insecticides, surviving and adapting with the aid of its expanded detoxification systems. The findings of this study will enable the development of new pest management strategies for the control of major agricultural pests such as S. litura.

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.

Insect pests have developed resistance to chemical insecticides, insecticidal toxins as bioinsecticides or genetic protection built into crops. Consequently, novel, orally active insecticidal toxins would be valuable biological alternatives for pest control. Here, we identified a novel insecticidal toxin, parasporal crystal toxin (PC), from Bacillus bombysepticus (Bb). PC shows oral pathogenic activity and lethality towards silkworms and Cry1Ac-resistant Helicoverpa armigera strains. In vitro assays, PC after activated by trypsin binds to BmAPN4 and BtR-175 by interacting with CR7 and CR12 fragments. Additionally, trypsin-activated PC demonstrates cytotoxicity against Sf9 cells expressing BmAPN4, revealing that BmAPN4 serves as a functional receptor that participates in Bb and PC pathogenicity. In vivo assay, knocking out BtR-175 increased the resistance of silkworms to PC. These data suggest that PC is the first protein with insecticidal activity identified in Bb that is capable of causing silkworm death via receptor interactions, representing an important advance in our understanding of the toxicity of Bb and the contributions of interactions between microbial pathogens and insects to disease pathology. Furthermore, the potency of PC as an insecticidal protein makes it a good candidate for inclusion in integrated agricultural pest management systems.

The wild silkworm Bombyx mandarina is widely believed to be an ancestor of the domesticated silkworm, Bombyx mori. Silkworms are often used as a model for studying the mechanism of species domestication. Here, we performed transcriptome sequencing of the wild silkworm using an Illumina HiSeq2000 platform. We produced 100,004,078 high-quality reads and assembled them into 50,773 contigs with an N50 length of 1764 bp and a mean length of 941.62 bp. A total of 33,759 unigenes were identified, with 12,805 annotated in the Nr database, 8273 in the Pfam database, and 9093 in the Swiss-Prot database. Expression profile analysis found significant differential expression of 1308 unigenes between the middle silk gland (MSG) and posterior silk gland (PSG). Three sericin genes (sericin 1, sericin 2, and sericin 3) were expressed specifically in the MSG and three fibroin genes (fibroin-H, fibroin-L, and fibroin/P25) were expressed specifically in the PSG. In addition, 32,297 Single-nucleotide polymorphisms (SNPs) and 361 insertion-deletions (INDELs) were detected. Comparison with the domesticated silkworm p50/Dazao identified 5,295 orthologous genes, among which 400 might have experienced or to be experiencing positive selection by Ka/Ks analysis. These data and analyses presented here provide insights into silkworm domestication and an invaluable resource for wild silkworm genomics research.

Bacillus bombysepticus is a Gram-positive spore-forming bacterium. Here, we announce the first complete genome sequence of this organism isolated from the cadavers of silkworm larvae that had been sick. The genome contains a single circular chromosome and a circular plasmid. Analyses of the B. bombysepticus genome will provide insights into its pathomechanisms and biology.

The ATP-binding cassette (ABC) superfamily is a large protein family with diverse physiological functions in all kingdoms of life. One distinguished subfamily, the pleiotropic drug resistance (PDR) transporters, has only been identified in plants and fungi. Here, we identified a Nicotiana tabacum PDR gene, NtPDR6, which is a homolog of Petunia hybrida PDR1. The full-length cDNA of NtPDR6 had a 4482-bp open reading frame encoding a full-size ABC transporter with 1493 amino acids. Sequence comparison showed that NtPDR6 had high homology with plant PDR proteins. NtPDR6 was strongly induced by phosphate starvation as well as by 1-naphthalene acetic acid. Tissue expression pattern analysis showed that NtPDR6 was detected in all surveyed tissues but preferentially in roots. We cloned the 1.3-kb NtPDR6 promoter and found that there was one phosphate starvation response-related element Pho-like and several root-specific expression-related elements rootmotiftapox1 in the NtPDR6 promoter. A tissue-specific pattern of NtPDR6 promoter-β-glucuronidase expression was dominantly observed in subepidermal cells and the elongation zone of lateral roots. RNA interference technology was used to knock down NtPDR6 expression, and there was a significantly increased branching phenotype in the NtPDR6 knockdown plants. These data suggest that NtPDR6 plays a key role in regulation of shoot branching processes.

In the silkworm Bombyx mori, pigment mutants with diverse body colors have been maintained throughout domestication for about 5000 years. The silkworm larval body color is formed through the mutual interaction of melanin, ommochromes, pteridines and uric acid. These pigments/compounds are synthesized by the cooperative action of various genes and enzymes. Previous reports showed that melanin, ommochrome and pteridine are increased in silkworm quail (q) mutants. To understand the pigment increase and alterations in pigment synthesis in q mutant, transcriptome profiles of the silkworm integument were investigated at 16 h after head capsule slippage in the fourth molt in q mutants and wild-type (Dazao). Compared to the wild-type, 1161 genes were differentially expressed in the q mutant. Of these modulated genes, 62.4% (725 genes) were upregulated and 37.6% (436 genes) were downregulated in the q mutant. The molecular function of differently expressed genes was analyzed by Blast2GO. The results showed that upregulated genes were mainly involved in protein binding, small molecule binding, transferase activity, nucleic acid binding, specific DNA-binding transcription factor activity and chromatin binding, while exclusively down-expressed genes functioned in oxidoreductase activity, cofactor binding, tetrapyrrole binding, peroxidase activity and pigment binding. We focused on genes related to melanin, pteridine and ommochrome biosynthesis; transport of uric acid; and juvenile hormone metabolism because of their importance in integument coloration during molting. This study identified differently expressed genes implicated in silkworm integument formation and pigmentation using silkworm q mutant. The results estimated the number and types of genes that drive new integument formation.

Long noncoding RNAs (lncRNAs) are emerging as important regulators in various biological processes. However, to date, no systematic characterization of lncRNAs has been reported in the silkworm Bombyx mori. In the present study, we generated eighteen RNA-seq datasets with relatively high depth. Using an in-house designed lncRNA identification pipeline, 11,810 lncRNAs were identified for 5,556 loci. Among these lncRNAs, 474 transcripts were intronic lncRNAs (ilncRNAs), 6,250 transcripts were intergenic lncRNAs (lincRNAs), and 5,086 were natural antisense lncRNAs (lncNATs). Compared with protein-coding mRNAs, silkworm lncRNAs are shorter in terms of full length but longer in terms of exon and intron length. In addition, lncRNAs exhibit a lower level of sequence conservation, more repeat sequences overlapped and higher tissue specificity than protein-coding mRNAs in the silkworm. We found that 69 lncRNA transcripts from 33 gene loci may function as miRNA precursors, and 104 lncRNA transcripts from 72 gene loci may act as competing endogenous RNAs (ceRNAs). In total, 49.47% of all gene loci (2,749/5,556) for which lncRNAs were identified showed sex-biased expression. Co-expression network analysis resulted in 19 modules, 12 of which revealed relatively high tissue specificity. The highlighted darkgoldenrod module was specifically associated with middle and posterior silk glands, and the hub lncRNAs within this module were co-expressed with proteins involved in translation, translocation, and secretory processes, suggesting that these hub lncRNAs may function as regulators of the biosynthesis, translocation, and secretion of silk proteins. This study presents the first comprehensive genome-wide analysis of silkworm lncRNAs and provides an invaluable resource for genetic, evolutionary, and genomic studies of B. mori.

Several pathogenic microorganisms have been used to investigate the genome-wide transcriptional responses of Bombyx mori to infection. However, studies have so far each focused on one microorganism, and systematic genome-wide comparison of transcriptional responses to different pathogenic microorganisms has not been undertaken. Here, we surveyed transcriptional responses of B. mori to its natural bacterial, viral, and fungal pathogens, Bacillus bombyseptieus, B. mori nucleopolyhedrovirus (BmNPV), and Beauveria bassiana, respectively, and to nonpathogenic Escherichia coli, by microarray analysis. In total, the expression of 2,436, 1,804, 1,743, and 912 B. mori genes was modulated by infection with B. bombyseptieus, BmNPV, B. bassiana, and E. coli, respectively. Notably, the expression of 620, 400, 177, or 165 of these genes was only modulated by infection with B. bombyseptieus, BmNPV, B. bassiana, or E. coli, respectively. In contrast to the expression of genes related to juvenile hormone synthesis and metabolism, that of genes encoding juvenile hormone binding proteins was microorganism-specific. Three basal metabolic pathways were modulated by infection with any of the four microorganisms, and 3, 14, 5, and 2 metabolic pathways were specifically modulated by infection with B. bombyseptieus, BmNPV, B. bassiana, and E. coli, respectively. Interestingly, BmNPV infection modulated the JAK/STAT signaling pathway, whereas both the Imd and Toll signaling pathways were modulated by infection with B. bombyseptieus, B. bassiana, or E. coli. These results elucidate potential molecular mechanisms of the host response to different microorganisms, and provide a foundation for further work on host–pathogen interaction.