Experimental RNomics and genomic comparative analysis reveal a large group of species-specific small non-message RNAs in the silkworm Bombyx mori.
Ontology highlight
ABSTRACT:
SUBMITTER: Li D
PROVIDER: S-EPMC3089462 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BackgroundSmall non-coding RNAs (ncRNAs) are important regulators of gene expression in eukaryotes. Previously, only microRNAs (miRNAs) and piRNAs have been identified in the silkworm, Bombyx mori. Furthermore, only ncRNAs (50-500nt) of intermediate size have been systematically identified in the silkworm.ResultsHere, we performed a systematic identification and analysis of small RNAs (18-50nt) associated with the Bombyx mori argonaute2 (BmAgo2) protein. Using RIP-seq, we identified various types of small ncRNAs associated with BmAGO2. These ncRNAs showed a multimodal length distribution, with three peaks at ~20nt, ~27nt and ~33nt, which included tRNA-, transposable element (TE)-, rRNA-, snoRNA- and snRNA-derived small RNAs as well as miRNAs and piRNAs. The tRNA-derived fragments (tRFs) were found at an extremely high abundance and accounted for 69.90% of the BmAgo2-associated small RNAs. Northern blotting confirmed that many tRFs were expressed or up-regulated only in the BmNPV-infected cells, implying that the tRFs play a prominent role by binding to BmAgo2 during BmNPV infection. Additional evidence suggested that there are potential cleavage sites on the D, anti-codon and TψC loops of the tRNAs. TE-derived small RNAs and piRNAs also accounted for a significant proportion of the BmAgo2-associated small RNAs, suggesting that BmAgo2 could be involved in the maintenance of genome stability by suppressing the activities of transposons guided by these small RNAs. Finally, Northern blotting was also used to confirm the Bombyx 5.8 s rRNA-derived small RNAs, demonstrating that various novel small RNAs exist in the silkworm.ConclusionsUsing an RIP-seq method in combination with Northern blotting, we identified various types of small RNAs associated with the BmAgo2 protein, including tRNA-, TE-, rRNA-, snoRNA- and snRNA-derived small RNAs as well as miRNAs and piRNAs. Our findings provide new clues for future functional studies of the role of small RNAs in insect development and evolution.
| S-EPMC3849828 | biostudies-literature
Project description: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.
| S-EPMC4714849 | biostudies-literature
Project description: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.
| S-EPMC8228281 | biostudies-literature
Project description:BackgroundThe domesticated silkworm, Bombyx mori, is the model insect for the order Lepidoptera, has economically important values, and has gained some representative behavioral characteristics compared to its wild ancestor. The genome of B. mori has been fully sequenced while function analysis of BmChi-h and BmSuc1 genes revealed that horizontal gene transfer (HGT) maybe bestow a clear selective advantage to B. mori. However, the role of HGT in the evolutionary history of B. mori is largely unexplored. In this study, we compare the whole genome of B. mori with those of 382 prokaryotic and eukaryotic species to investigate the potential HGTs.ResultsTen candidate HGT events were defined in B. mori by comprehensive sequence analysis using Maximum Likelihood and Bayesian method combining with EST checking. Phylogenetic analysis of the candidate HGT genes suggested that one HGT was plant-to- B. mori transfer while nine were bacteria-to- B. mori transfer. Furthermore, functional analysis based on expression, coexpression and related literature searching revealed that several HGT candidate genes have added important characters, such as resistance to pathogen, to B. mori.ConclusionsResults from this study clearly demonstrated that HGTs play an important role in the evolution of B. mori although the number of HGT events in B. mori is in general smaller than those of microbes and other insects. In particular, interdomain HGTs in B. mori may give rise to functional, persistent, and possibly evolutionarily significant new genes.
| S-EPMC3116507 | biostudies-literature
Project description:Gloverin is one of the glycine-rich antimicrobial peptide exclusively found in Lepidoptera insects. It is generally activated through the innate immune system in insects. In this study, recombinant Gloverin2 from Bombyx mori (BmGlv2) was synthesized using a prokaryotic expression system. Circular dichroism spectroscopy showed that the recombinant BmGlv2 has random coil structure, which is relatively stable at the temperatures ranging from 15 to 82.5 °C. Antimicrobial activity analysis revealed that BmGlv2 significantly inhibited the growth of gram-negative bacteria, Escherichia coli JM109 and Pseudomonas putida, by disrupting cell integrity. Western blotting and immunofluorescence analyses suggested that BmGlv2 absorbed on the cell surface after incubation, which might be the first step in the antibacterial process. Our results also proved that the cell wall component lipopolysaccharides (LPS) induce a conformational change in BmGlv2 from a random coil to ?-helix. Subsequently, ?-helical BmGlv2 would recruit more BmGlv2 and form higher aggregation state. Collectively, these findings expand our understanding of antibacterial mechanism of BmGlv2.
| S-EPMC6121877 | biostudies-literature
Project description:BackgroundIn our previous study, we identified four isoforms of the Bmovo gene, Bmovo-1, Bmovo-2, Bmovo-3 and Bmovo-4 from the silkworm ovary and verified that ovarian development was regulated by the BmOVO proteins.ResultsTo understand the regulatory mechanisms of ovarian development, the regulation of four BmOVO isoforms on the B. mori ovarian tumor (Bmotu) promoter activity was investigated with luciferase reporter assays. The results showed the Bmotu promoter activity was positively regulated by BmOVO-1, BmOVO-2, BmOVO-3 and BmOVO-4 in a dose-dependent manner, of which BmOVO-2 had the highest transcriptional activation. However, the first (A1) and third acidic domains (A3) at the N-terminus of BmOVO-1 are transcriptional repression domains, while the fourth (A4) and fifth acidic domains (A5) are transcriptional activation domains. A recombinant BmOVO zinc-finger domain was found to bind to the GTACCGTTGTA sequence located at the Bmotu promoter. Furthermore, the Bmotu promoter activity was negatively regulated by 'Tal-like' peptide, which can trigger BmOVO-1 degradation at the N-terminus.ConclusionsThese results will help us to further understand the regulatory mechanisms of BmOVO isoforms on Bmotu promoter activity and ovarian development in the silkworm.
| S-EPMC6503385 | biostudies-literature
Project description:To investigate the patterns of nucleotide diversity in domesticated silkworm, Bombyx mori L. (Lepidoptera: Bombycidae) and its wild relative, Chinese wild silkworm, Bombyx mandarina Moore, we sequenced nine nuclear genes. Neutrality test and coalescent simulation for these genes were performed to look at bottleneck intensity and selection signature; linkage disequilibrium (LD) within and between loci was employed to investigate allele association. As a result, B. mori lost 33-49% of nucleotide diversity relative to wild silkworm, which is similar to the loss levels found in major cultivated crops. Diversity of B. mori is significantly lower than that of B. mandarina measured as π(total) (0.01166 vs. 0.1741) or θ(W)(0.01124 vs. 0.02206). Bottleneck intensity of domesticated silkworm is 1.5 (in terms of k = N(b) /d, N(b) -bottleneck population size; d-bottleneck duration) with different durations. Gene DefA showed signature of artificial selection by all analysis methods and might experience strong artificial selection in B. mori during domestication. For nine loci, both curves of LD decay rapidly within 200 bp and drop slowly when distance is > 200 bp, although that of B. mori decays slower than B. mandarina at loci investigated. However, LD could not be estimated at DefA in B. mori and at ER in both silkworms. Elevated LD observed in B. mori may be indicator of selection and demographic events.
| S-EPMC3391917 | biostudies-literature
Project description:Bombyx mori densovirus 1 (BmDNV-1), a major pathogen of silkworms, causes significant losses to the silk industry. The structure of the recombinant BmDNV-1 virus-like particle has been determined at 3.1-Å resolution using X-ray crystallography. It is the first near-atomic-resolution structure of a virus-like particle within the genus Iteravirus. The particles consist of 60 copies of the 55-kDa VP3 coat protein. The capsid protein has a β-barrel "jelly roll" fold similar to that found in many diverse icosahedral viruses, including archaeal, bacterial, plant, and animal viruses, as well as other parvoviruses. Most of the surface loops have little structural resemblance to other known parvovirus capsid proteins. In contrast to vertebrate parvoviruses, the N-terminal β-strand of BmDNV-1 VP3 is positioned relative to the neighboring 2-fold related subunit in a "domain-swapped" conformation, similar to findings for other invertebrate parvoviruses, suggesting domain swapping is an evolutionarily conserved structural feature of the Densovirinae.
| S-EPMC3126189 | biostudies-literature
Project description:BackgroundInsect cuticle plays essential roles in many physiological functions. During molting and metamorphosis tremendous changes occur in silkworm cuticle where multiple proteins exist and genes encoding them constitute about 1.5% of all Bombyx mori genes.ResultsIn an effort to determine their expression profiles, a microarray-based investigation was carried out using mRNA collected from larvae to pupae. The results showed that a total of 6676 genes involved in various functions and physiological pathways were activated. The vast majority (93%) of cuticular protein genes were expressed in selected stages with varying expression patterns. There was no correlation between expression patterns and the presence of conserved motifs. Twenty-six RR genes distributed in chromosome 22 were co-expressed at the larval and wandering stages. The 2 kb upstream regions of these genes were further analyzed and three putative elements were identified.ConclusionsData from the present study provide, for the first time, a comprehensive expression profile of genes in silkworm epidermal tissues and evidence that putative elements exist to allow massive production of mRNAs from specific cuticular protein genes.
| S-EPMC2848646 | biostudies-literature
Project description:The pathogenesis of Bombyx mori cytoplasmic polyhedrosis virus (BmCPV) infection is unclear, although accumulating evidence indicates that circular RNAs (circRNAs), which act as competing endogenous RNAs or positive regulators, play important roles in regulating gene expression in eukaryotes and, thus, may play a role in BmCPV infections. To explore the expression and biological functions of circRNAs in the silkworm midgut following BmCPV infection, silkworm circRNA expression profiles of normal midgut tissue (control) and BmCPV-infected midgut tissue (test) were determined using high-through sequencing. A total of 9,753 and 7,475circRNAs were detected from the control and test samples, respectively. The two samples shared 6,085 circRNAs, while 646 and 737 circRNAs were expressed specifically in the control and test samples, respectively. A total of 3,638 circRNAs were shown to be differentially expressed, and 400 circRNAs were substantially differentially expressed with a fold-change ≥ 2.0 (p< 0.05 and a false discover rate < 0.05), of which 294 were up-regulated and 106 were down-regulated following infection. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were conducted to determine the principal functions of the substantially differentially regulated genes. circRNA-miRNA interaction networks were constructed based on a correlation analysis between the differentially expressed circRNAs and the nature of their microRNA (miRNA) binding sites. The network inferred that 13 miRNAs interacting with 193 circRNAs were among the 300 most abundant relationships. bmo-miR-3389-5p, bmo-miR-745-3p, and bmo-miR-3262 were related to 30, 34, and 34 circRNAs, respectively. circRNA_8115, circRNA_9444, circRNA_4553, circRNA_0827, and circRNA_6649 contained six, five, four, four, and four miRNA binding sites, respectively. We further found that alternative circularization of circRNAs is a common feature in silkworms and that the junction sites of many silkworm circRNAs are flanked by canonical GT/AG splicing signals. Our study is the first to show the circRNA response to virus infection. Thus, it provides a novel perspective on circRNA-miRNA interactions during BmCPV pathogenesis, and it lays the foundation for future research of the potential roles of circRNAs in BmCPV pathogenesis.
| S-EPMC5798952 | biostudies-literature