Confluence and convergence of Dscam and Pcdh cell-recognition codes.

The ability of neurites of the same neuron to avoid each other (self-avoidance) is a conserved feature in both invertebrates and vertebrates. The key to self-avoidance is the generation of a unique subset of cell-surface proteins in individual neurons engaging in isoform-specific homophilic interactions that drive neurite repulsion rather than adhesion. Among these cell-surface proteins are fly Dscam1 and vertebrate clustered protocadherins (cPcdhs), as well as the recently characterized shortened Dscam (sDscam) in the Chelicerata. Herein, we review recent advances in our understanding of how cPcdh, Dscam, and sDscam cell-surface recognition codes are expressed and translated into cellular functions essential for neural wiring.

A systematic CRISPR screen reveals redundant and specific roles for Dscam1 isoform diversity in neuronal wiring.

Drosophila melanogaster Down syndrome cell adhesion molecule 1 (Dscam1) encodes 19,008 diverse ectodomain isoforms via the alternative splicing of exon 4, 6, and 9 clusters. However, whether individual isoforms or exon clusters have specific significance is unclear. Here, using phenotype-diversity correlation analysis, we reveal the redundant and specific roles of Dscam1 diversity in neuronal wiring. A series of deletion mutations were performed from the endogenous locus harboring exon 4, 6, or 9 clusters, reducing to 396 to 18,612 potential ectodomain isoforms. Of the 3 types of neurons assessed, dendrite self/non-self discrimination required a minimum number of isoforms (approximately 2,000), independent of exon clusters or isoforms. In contrast, normal axon patterning in the mushroom body and mechanosensory neurons requires many more isoforms that tend to associate with specific exon clusters or isoforms. We conclude that the role of the Dscam1 diversity in dendrite self/non-self discrimination is nonspecifically mediated by its isoform diversity. In contrast, a separate role requires variable domain- or isoform-related functions and is essential for other neurodevelopmental contexts, such as axonal growth and branching. Our findings shed new light on a general principle for the role of Dscam1 diversity in neuronal wiring.

Chelicerata sDscam isoforms combine homophilic specificities to define unique cell recognition.

Thousands of Down syndrome cell adhesion molecule (Dscam1) isoforms and ∼60 clustered protocadhrein (cPcdh) proteins are required for establishing neural circuits in insects and vertebrates, respectively. The strict homophilic specificity exhibited by these proteins has been extensively studied and is thought to be critical for their function in neuronal self-avoidance. In contrast, significantly less is known about the Dscam1-related family of ∼100 shortened Dscam (sDscam) proteins in Chelicerata. We report that Chelicerata sDscamα and some sDscamß protein trans interactions are strictly homophilic, and that the trans interaction is meditated via the first Ig domain through an antiparallel interface. Additionally, different sDscam isoforms interact promiscuously in cis via membrane proximate fibronectin-type III domains. We report that cell-cell interactions depend on the combined identity of all sDscam isoforms expressed. A single mismatched sDscam isoform can interfere with the interactions of cells that otherwise express an identical set of isoforms. Thus, our data support a model by which sDscam association in cis and trans generates a vast repertoire of combinatorial homophilic recognition specificities. We propose that in Chelicerata, sDscam combinatorial specificity is sufficient to provide each neuron with a unique identity for self-nonself discrimination. Surprisingly, while sDscams are related to Drosophila Dscam1, our results mirror the findings reported for the structurally unrelated vertebrate cPcdh. Thus, our findings suggest a remarkable example of convergent evolution for the process of neuronal self-avoidance and provide insight into the basic principles and evolution of metazoan self-avoidance and self-nonself discrimination.

RNA secondary structures in Dscam1 mutually exclusive splicing: unique evolutionary signature from the midge.

The Drosophila melanogaster gene Dscam1 potentially generates 38,016 distinct isoforms via mutually exclusive splicing, which are required for both nervous and immune functions. However, the mechanism underlying splicing regulation remains obscure. Here we show apparent evolutionary signatures characteristic of competing RNA secondary structures in exon clusters 6 and 9 of Dscam1 in the two midge species (Belgica antarctica and Clunio marinus). Surprisingly, midge Dscam1 encodes only ∼6000 different isoforms through mutually exclusive splicing. Strikingly, the docking site of the exon 6 cluster is conserved in almost all insects and crustaceans but is specific in the midge; however, the docking site-selector base-pairings are conserved. Moreover, the docking site is complementary to all predicted selector sequences downstream from every variable exon 9 of the midge Dscam1, which is in accordance with the broad spectrum of their isoform expression. This suggests that these cis-elements mainly function through the formation of long-range base-pairings. This study provides a vital insight into the evolution and mechanism of Dscam1 alternative splicing.

Long Noncoding RNA GMAN, Up-regulated in Gastric Cancer Tissues, Is Associated With Metastasis in Patients and Promotes Translation of Ephrin A1 by Competitively Binding GMAN-AS.

BACKGROUND & AIMS: We aimed to identify long noncoding RNAs (lncRNAs) that are up-regulated in gastric cancer tissues from patients and study their function in gastric tumor metastasis. METHODS: We collected gastric tumor and nontumor tissues from patients in China and analyzed levels of lncRNAs by microarray analysis, proteins by immunohistochemistry, and RNAs by quantitative reverse-transcription polymerase chain reaction; we compared these with survival times of patients and tumor progression. RNA levels were knocked down or knocked out in BGC-823, SGC-7901, and MKN45 cell lines using small interfering or short hairpin RNAs or clustered regularly interspaced short palindromic repeats (ie, CRISPR)/CRISPR associated protein 9 (ie, Cas9) vectors. Genes were overexpressed from transfected plasmids in HGC-27 cells. Cells were analyzed by Northern blot and immunoblot, polysome profiling assay, and cell invasion assay. Cells were injected into the tail veins or spleens of nude mice or SCID mice; lung and liver tissues were collected, and metastases were counted. lncRNAs were cloned by using rapid amplification of complementary DNA ends. Their interactions with other genes were determined by RNA pulldown and mapping assays. RESULTS: In microarray analyses, we identified 151 lncRNAs expressed at significantly higher levels in gastric tumor vs nontumor tissues. Levels of an lncRNA that we called gastric cancer metastasis associated long noncoding RNA (GMAN) were increased in gastric tumor tissues, compared with nontumor tissues; its up-regulation was associated with tumor metastasis and shorter survival times of patients. The GMAN gene overlaps with the ephrin A1 gene (EFNA1) and was highly expressed in BGC-823 and MKN45 cells. Knockdown of GMAN in these cells did not affect proliferation, colony formation, or adhesion but did reduce their invasive activity in Transwell assays. Ectopic expression of GMAN increased the invasive activity of HGC-27 cells. BGC-823 and MKN45 cells with knockdown of GMAN formed fewer metastases after injection into tail veins of nude mice. Knockdown or knockout of GMAN also reduced levels of ephrin A1 protein in cells. We found that GMAN promoted translation of ephrin A1 messenger RNA into protein by binding to the antisense GMAN RNA (GMAN-AS)-this antisense sequence is also complementary to that of ephrin A1 mRNA. Levels of ephrin A1 protein were also increased in gastric tumors from patients with metastases than in those without metastases. Knockout of ephrin A1 in BGC-823 cells reduced their invasive activity in Transwell assays and ability to form metastases after injection into SCID mice. Ectopic expression of ephrin A1 in BGC-823 cells with knockdown or knockout of GMAN restored their invasive activities and ability form metastases in nude or SCID mice. A CRISPR/Cas9-based strategy to disrupt the GMAN gene significantly reduced the numbers of metastases formed from SGC-7901 cells in mice. CONCLUSIONS: We identified an lncRNA, which we call GMAN, that is increased in gastric tumors from patients and associated with survival and formation of metastases. It regulates translation of ephrin A1 mRNA by binding competitively to GMAN-AS. Knockdown or knockout of GMAN or ephrin A1 in gastric cancer cell lines reduces their invasive activity and ability to form metastases after injection into mice. These genes might be targeted to prevent or reduce gastric cancer metastasis.

Competing RNA pairings in complex alternative splicing of a 3' variable region.

Alternative pre-mRNA splicing remarkably expands protein diversity in eukaryotes. Drosophila PGRP-LC can generate three major 3' splice isoforms that exhibit distinct innate immune recognition and defenses against various microbial infections. However, the regulatory mechanisms underlying the uniquely biased splicing pattern at the 3' variable region remain unclear. Here we show that competing RNA pairings control the unique splicing of the 3' variable region of Drosophila PGRP-LC pre-mRNA. We reveal three roles by which these RNA pairings jointly regulate the 3' variant selection through activating the proximal 3' splice site and concurrently masking the intron-proximal 5' splice site, in combination with physical competition of RNA pairing. We also reveal that competing RNA pairings regulate alternative splicing of the highly complex 3' variable regions of Drosophila CG42235 and Pip Our findings will facilitate a better understanding of the regulatory mechanisms of highly complex alternative splicing as well as highly variable 3' processing.

Revisiting Dscam diversity: lessons from clustered protocadherins.

The complexity of neuronal wiring relies on the extraordinary recognition diversity of cell surface molecules. Drosophila Dscam1 and vertebrate clustered protocadherins (Pcdhs) are two classic examples of the striking diversity from a complex genomic locus, wherein the former encodes more than 10,000 distinct isoforms via alternative splicing, while the latter employs alternative promoters to attain isoform diversity. These structurally unrelated families show remarkably striking molecular parallels and even similar functions. Recent studies revealed a novel Dscam gene family with tandemly arrayed 5' cassettes in Chelicerata (e.g., the scorpion Mesobuthus martensii and the tick Ixodes scapularis), similar to vertebrate clustered Pcdhs. Likewise, octopus shows a more remarkable expansion of the Pcdh isoform repertoire than human. These discoveries of Dscam and Pcdh diversification reshape the evolutionary landscape of recognition molecule diversity and provide a greater understanding of convergent molecular strategies for isoform diversity. This article reviews new insights into the evolution, regulatory mechanisms, and functions of Dscam and Pcdh isoform diversity. In particular, the convergence of clustered Dscams and Pcdhs is highlighted.

A chelicerate-specific burst of nonclassical Dscam diversity.

BACKGROUND: The immunoglobulin (Ig) superfamily receptor Down syndrome cell adhesion molecule (Dscam) gene can generate tens of thousands of isoforms via alternative splicing, which is essential for both nervous and immune systems in insects. However, further information is required to develop a comprehensive view of Dscam diversification across the broad spectrum of Chelicerata clades, a basal branch of arthropods and the second largest group of terrestrial animals. RESULTS: In this study, a genome-wide comprehensive analysis of Dscam genes across Chelicerata species revealed a burst of nonclassical Dscams, categorised into four types-mDscam, sDscamα, sDscamß, and sDscamγ-based on their size and structure. Although the mDscam gene class includes the highest number of Dscam genes, the sDscam genes utilise alternative promoters to expand protein diversity. Furthermore, we indicated that the 5' cassette duplicate is inversely correlated with the sDscam gene duplicate. We showed differential and sDscam- biased expression of nonclassical Dscam isoforms. Thus, the Dscam isoform repertoire across Chelicerata is entirely dominated by the number and expression levels of nonclassical Dscams. Taken together, these data show that Chelicerata evolved a large conserved and lineage-specific repertoire of nonclassical Dscams. CONCLUSIONS: This study showed that arthropods have a large diversified Chelicerata-specific repertoire of nonclassical Dscam isoforms, which are structurally and mechanistically distinct from those of insects. These findings provide a global framework for the evolution of Dscam diversity in arthropods and offer mechanistic insights into the diversification of the clade-specific Ig superfamily repertoire.

Long-range RNA pairings contribute to mutually exclusive splicing.

Mutually exclusive splicing is an important means of increasing the protein repertoire, by which the Down's syndrome cell adhesion molecule (Dscam) gene potentially generates 38,016 different isoforms in Drosophila melanogaster. However, the regulatory mechanisms remain obscure due to the complexity of the Dscam exon cluster. Here, we reveal a molecular model for the regulation of the mutually exclusive splicing of the serpent pre-mRNA based on competition between upstream and downstream RNA pairings. Such dual RNA pairings confer fine tuning of the inclusion of alternative exons. Moreover, we demonstrate that the splicing outcome of alternative exons is mediated in relative pairing strength-correlated mode. Combined comparative genomics analysis and experimental evidence revealed similar bidirectional structural architectures in exon clusters 4 and 9 of the Dscam gene. Our findings provide a novel mechanistic framework for the regulation of mutually exclusive splicing and may offer potentially applicable insights into long-range RNA-RNA interactions in gene regulatory networks.

Predicted networks of protein-protein interactions in Stegodyphus mimosarum by cross-species comparisons.

BACKGROUND: Stegodyphus mimosarum is a candidate model organism belonging to the class Arachnida in the phylum Arthropoda. Studies on the biology of S. mimosarum over the past several decades have consisted of behavioral research and comparison of gene sequences based on the assembled genome sequence. Given the lack of systematic protein analyses and the rich source of information in the genome, we predicted the relationships of proteins in S. mimosarum by bioinformatics comparison with genome-wide proteins from select model organisms using gene mapping. RESULTS: The protein-protein interactions (PPIs) of 11 organisms were integrated from four databases (BioGrid, InAct, MINT, and DIP). Here, we present comprehensive prediction and analysis of 3810 proteins in S. mimosarum with regard to interactions between proteins using PPI data of organisms. Interestingly, a portion of the protein interactions conserved among Saccharomyces cerevisiae, Homo sapiens, Arabidopsis thaliana, and Drosophila melanogaster were found to be associated with RNA splicing. In addition, overlap of predicted PPIs in reference organisms, Gene Ontology, and topology models in S. mimosarum are also reported. CONCLUSIONS: Addition of Stegodyphus, a spider representative of interactomic research, provides the possibility of obtaining deeper insights into the evolution of PPI networks among different animal species. This work largely supports the utility of the "stratus clouds" model for predicted PPIs, providing a roadmap for integrative systems biology in S. mimosarum.

Role and convergent evolution of competing RNA secondary structures in mutually exclusive splicing.

Exon or cassette duplication is an important means of expanding protein and functional diversity through mutually exclusive splicing. However, the mechanistic basis of this process in non-arthropod species remains poorly understood. Here, we demonstrate that MRP1 genes underwent tandem exon duplication in Nematoda, Platyhelminthes, Annelida, Mollusca, Arthropoda, Echinodermata, and early-diverging Chordata but not in late-diverging vertebrates. Interestingly, these events were of independent origin in different phyla, suggesting convergent evolution of alternative splicing. Furthermore, we showed that multiple sets of clade-conserved RNA pairings evolved to guide species-specific mutually exclusive splicing in Arthropoda. Importantly, we also identified a similar structural code in MRP exon clusters of the annelid, Capitella teleta, and chordate, Branchiostoma belcheri, suggesting an evolutionarily conserved competing pairing-guided mechanism in bilaterians. Taken together, these data reveal the molecular determinants and RNA pairing-guided evolution of species-specific mutually exclusive splicing spanning more than 600 million years of bilaterian evolution. These findings have a significant impact on our understanding of the evolution of and mechanism underpinning isoform diversity and complex gene structure.

BIOINFORMATICS, TISSUE DISTRIBUTION, AND SUBCELLULAR LOCALIZATION ANALYSES OF FK506 BINDING PROTEIN 12B FROM SILKWORMS.

FK506 binding proteins (FKBPs) are intracellular receptors of the immunosuppressant FK506 and play important roles in the correct folding of new proteins and the self-assembly of biological macromolecules. FKBP12 is a member of the FKBP family that is widely expressed and highly conserved in many species. In this study, we identified the complete cDNA sequence encoding the FKBP12 ortholog in Bombyx mori, named Bm-FKBP12B (GenBank accession no. DQ443423). Multiple-sequence alignment among different species revealed a high similarity among FKBP12 paralogs and orthologs. Bioinformatics analysis of the Bm-FKBP12B gene showed that it is located on chromosome 20 and consists of three exons and two introns. We cloned, expressed, and purified the Bm-FKBP12B protein in Escherichia coli and generated a specific polyclonal antibody against Bm-FKBP12B. The real-time quantitative reverse-transcription (qRT) PCR and Western blotting results showed that Bm-FKBP12B was present throughout all of the development stages, but it was abundant in the adult and embryo stages. Bm-FKBP12B expression was higher in the silk gland and gut, suggesting that it might play important roles in regulating gene expression in the silk gland and during silk fiber formation. Bm-FKBP12B protein was distributed in the cytoplasm, nucleus, and nuclear membrane.

Treatment of type 1 myotonic dystrophy by engineering site-specific RNA endonucleases that target (CUG)(n) repeats.

Myotonic dystrophy type 1 (DM1) is caused by the expansion of (CTG)n in the 3' untranslated region of the dystrophia myotonica-protein kinase (DMPK) gene, which is transcribed as (CUG)n repeats that accumulate in the nucleus. The RNA repeats specifically sequester or change the expression levels of several RNA-binding proteins, leading to aberrant splicing of many target genes. In this study, we developed artificial site-specific RNA endonucleases (ASREs) that specifically bind and cleave (CUG)n repeats RNA. We have generated one ASRE that can target the expanded RNA repeats in DM1 patient cells and specifically degrade the pathogenic DMPK messenger RNAs with minimal effect on wild-type alleles. Such ASRE treatment significantly decreased the number of nuclear foci in DM1 patient cells and can reverse the missplicing of many genes affected in DM1 patients. Taken together, the application of ASRE provides a new route of gene therapy for DM1 treatment.

Origin and evolution of a new retained intron on the vulcan gene in Drosophila melanogaster subgroup species.

Although numerous intron gains have been discovered, the mechanisms of intron creation have proven to be elusive. Previous study revealed that the vulcan gene of Drosophila melanogaster contained four exons in its coding region. In the current study, a newly created intron (Intron L) was identified on exon 2 of vulcan in D. melanogaster by comparing expression sequence tags. The RT-PCR experiment revealed that Intron L was associated with intron retention, in which two alternative transcripts of the gene differ by the inclusion or removal of an intron. It was found that Intron L was created by intronization of exonic sequence, and its donor and acceptor splice sites were created by synonymous mutation, leading to the origin of a new vulcan protein that is 22 amino acids shorter than the previously reported vulcan protein. Moreover, to track the origin of Intron L, 36 orthologous genes of species of Drosophila were cloned or annotated, and phylogenetic analysis was carried out. It indicated that the common ancestor of D. melangaster subgroup species created Intron L about 15 million years ago.

Long non-coding RNAs: a novel endogenous source for the generation of Dicer-like 1-dependent small RNAs in Arabidopsis thaliana.

The biological relevance of long non-coding RNAs (lncRNAs) is emerging. Whether the lncRNAs could form structured precursors for small RNAs (sRNAs) production remains elusive. Here, 172 713 DCL1 (Dicer-like 1)-dependent sRNAs were identified in Arabidopsis. Except for the sRNAs mapped onto the microRNA precursors, the remaining ones led us to investigate their originations. Intriguingly, 65 006 sRNAs found their loci on 5891 lncRNAs. These sRNAs were sent to AGO (Argonaute) enrichment analysis. As a result, 1264 sRNAs were enriched in AGO1, which were then subjected to target prediction. Based on degradome sequencing data, 109 transcripts were validated to be targeted by 96 sRNAs. Besides, 44 lncRNAs were targeted by 23 sRNAs. To further support the origination of the DCL1-dependent sRNAs from lncRNAs, we searched for the degradome-based cleavage signals at either ends of the sRNA loci, which were supposed to be produced during DCL1-mediated processing of the long-stem structures. As a result, 63 612 loci were supported by degradome signatures. Among these loci, 6606 reside within the dsRNA-seq (double-stranded RNA sequencing) read-covered regions of 100 nt or longer. These regions were subjected to secondary structure prediction. And, 43 regions were identified to be capable of forming highly complementary long-stem structures. We proposed that these local long-stem structures could be recognized by DCL1 for cropping, thus serving as the sRNA precursors. We hope that our study could inspire more research efforts to study on the biological roles of the lncRNAs in plants.

A tobacco homolog of DCN1 is involved in pollen development and embryogenesis.

KEY MESSAGE: We show that DCN1 binds ubiquitin and RUB/NEDD8, associates with cullin, and is functionally conserved. DCN1 activity is required for pollen development transitions and embryogenesis, and for pollen tube growth. Plant proteomes show remarkable plasticity in reaction to environmental challenges and during developmental transitions. Some of this adaptability comes from ubiquitin-mediated protein degradation regulated by cullin-RING E3 ubiquitin ligases (CRLs). CRLs are activated through modification of the cullin subunit with the ubiquitin-like protein RUB/NEDD8 by an E3 ligase called defective in cullin neddylation 1 (DCN1). Here we show that tobacco DCN1 binds ubiquitin and RUB/NEDD8 and associates with cullin. When knocked down by RNAi, tobacco pollen formation was affected and zygotic embryogenesis was blocked around the globular stage. Additionally, we found that RNAi of DCN1 inhibited the stress-triggered reprogramming of cultured microspores from their intrinsic gametophytic mode of development to an embryogenic state. This stress-induced developmental switch is a known feature in many important crops and leads ultimately to the formation of haploid embryos and plants. Compensating the RNAi effect by re-transformation with a promoter-silencing construct restored pollen development and zygotic embryogenesis, as well as the ability for stress-induced formation of embryogenic microspores. Overexpression of DCN1 accelerated pollen tube growth and increased the potential for microspore reprogramming. These results demonstrate that the biochemical function of DCN1 is conserved in plants and that its activity is involved in transitions during pollen development and embryogenesis, and for pollen tube growth.

A potential link between aromatics-induced oviposition repellency behaviors and specific odorant receptor of Aedes albopictus.

BACKGROUND: The Asian tiger mosquito, Aedes albopictus, is a competent vector for the spread of several viral arboviruses including dengue, chikungunya, and Zika. Several vital mosquito behaviors linked to survival and reproduction are primarily dependent on a sophisticated olfactory system for semiochemical perception. However, a limited number of studies has hampered our understanding of the relationship between the A. albopictus acute olfactory system and the complex chemical world. RESULTS: Here, we performed a qRT-PCR assay on antennae from A. albopictus of differing sex, age and physiological states, and found that AalbOr10 was enriched in blood-fed female mosquitoes. We then undertook single sensillum recording to de-orphan AalbOr10 using a panel of physiologically and behaviorally relevant odorants in a Drosophila 'empty neuron' system. The results indicated that AalbOr10 was activated by seven aromatic compounds, all of which hampered egg-laying in blood-fed mosquitoes. Furthermore, using a post-RNA interference oviposition assay, we found that reducing the transcript level of AalbOr10 affected repellent activity mediated by 2-ethylphenol at low concentrations (10-4 vol/vol). Computational modeling and molecular docking studies suggested that hydrogen bonds to Y68 and Y150 mediated the interaction of 2-ethylphenol with AalbOr10. CONCLUSION: We reveal a potential link between aromatics-induced oviposition repellency behaviors and a specific odorant receptor in A. albopictus. Our findings provide a foundation for identifying active semiochemicals for the monitoring or controlling of mosquito populations. © 2024 Society of Chemical Industry.

Construction of small RNA-mediated gene regulatory networks in the roots of rice (Oryza sativa).

BACKGROUND: The root systems play essential roles for plants to anchorage to the soil, and to exploit the mineral and water resources. The molecular mechanisms underlying root development have been extensively studied to improve root system architecture, especially for the crops. Several microRNA (miRNA) families have been demonstrated to be involved in plant root development. However, whether the other small RNA (sRNA) species, which occupy a dominant portion of the plant endogenous sRNA population, possess potential roles in root development remains unclear. RESULTS: In this study, by using sRNA high-throughput sequencing data, we made a comparison of the sRNA accumulation levels between the rice root tips and the whole roots. The sRNAs highly accumulated in the root tips and in the whole roots were extracted respectively. After Argonaute 1 (AGO1) enrichment analysis, the sRNAs with great potential of performing target cleavages were included for target prediction and degradome sequencing data-based validation. As a result, lists of the targets regulated by the AGO1-enriched sRNAs were obtained for both the root tips and the whole roots. Further evidences were identified from microarray data of the target genes to support some of the sRNA-target interactions. Specifically, the expression patterns of certain target genes in the root tips and the whole roots were contrary to those of the regulating sRNAs. Besides, several targets were indicated to play important roles in root development based on literature mining. CONCLUSIONS: Taken together, the regulatory networks mediated by the sRNAs highly accumulated in the root tips or in the whole roots could advance our current understanding of the sRNA-involved molecular mechanisms underlying rice root development. And, the sRNA-target lists could serve as the basis for further functional investigations.

RIP-seq of BmAgo2-associated small RNAs reveal various types of small non-coding RNAs in the silkworm, Bombyx mori.

BACKGROUND: Small 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. RESULTS: Here, 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. CONCLUSIONS: Using 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.