key: cord-0839522-9tq0pr9z authors: Huang, Songqian; Nishiumi, Shinya; Asaduzzaman, Md; Pan, Yida; Liu, Guanting; Yoshitake, Kazutoshi; Maeyama, Kaoru; Kinoshita, Shigeharu; Nagai, Kiyohito; Watabe, Shugo; Yoshida, Tetsuhiko; Asakawa, Shuichi title: Exosome-derived small non-coding RNAs reveal immune response upon grafting transplantation in Pinctada fucata (Mollusca) date: 2022-05-04 journal: Open biology DOI: 10.1098/rsob.210317 sha: 19705c23d3bc5d9af09a7bd9b07ddb117e53561f doc_id: 839522 cord_uid: 9tq0pr9z Exosomes, a subset of small extracellular vesicles, carry various nucleic acids, proteins, lipids, amino acids and metabolites. They function as a mode of intercellular communication and molecular transfer. Exosome cargo molecules, including small non-coding RNAs (sncRNAs), are involved in the immune response in various organisms. However, the role of exosome-derived sncRNAs in immune responses in molluscs remains unclear. Here, we aimed to reveal the sncRNAs involved in the immune response during grafting transplantation by the pearl oyster Pinctada fucata. Exosomes were successfully extracted from the P. fucata haemolymph during graft transplantation. Abundant microRNAs (miRNAs) and PIWI-interacting RNAs (piRNAs) were simultaneously discovered in P. fucata exosomes by small RNA sequencing. The expression patterns of the miRNAs and piRNAs at the grafting and initial stages were not substantially different, but varied significantly between the initial and later stages. Target prediction and functional analysis indicate that these miRNAs and piRNAs are related to immune response upon grafting transplantation, whereas piRNAs may also be associated with transposon silencing by targeting with genome transposon elements. This work provides the basis for a functional understanding of exosome-derived sncRNAs and helps to gain further insight into the PIWI/piRNA pathway function outside of germline cells in molluscs. on model animals or higher organisms involved in immune response, health and disease [6, 7] . With the development of biotechnology and the continuous expansion of research scope, exosomes have been detected in various non-model animal body fluids, including arthropods [8] , crustaceans [9] , molluscs [10] [11] [12] and teleosts [13, 14] . MicroRNAs (miRNAs) are small non-coding RNA (sncRNA) molecules that are 20-23 nucleotides (nt) in length. They post-transcriptionally regulate the expression of target genes by binding to the 3 0 -untranslated regions (3 0 UTRs) of mRNA, which are involved in cell proliferation, tissue differentiation, embryonic development, metabolism and apoptosis [15, 16] . PIWI-interacting RNAs ( piRNAs) comprise a class of sncRNA molecules that are 24-31 nt in length. miRNAs are processed from long, single-stranded transcripts to form a short hairpin structure associated with AGO subfamily members of the Argonaute proteins [17] , whereas piRNAs are generated from long, single-stranded precursors associated with PIWI subfamily members [18] . piRNA biogenesis is classified into two categories, namely the primary and secondary biogenesis pathways. In the primary biogenesis pathway, long piRNA precursors are transcribed from specific genomic loci called piRNA clusters, cleaved and modified by complicated biogenesis factors in the cytoplasm [19, 20] . Subsequently, primary piRNAs are subjected to an amplification system to enforce the high expression of piRNAs, which is called the amplification loop or the ping-pong cycle [21] . Cells secrete various EVs, including exosomes, into the internal environment for various forms of cell-to-cell communication [22] . Exosomes act as a natural material transport capsule for intercellular communication to realize the gene expression regulatory activity (RNAi) of sncRNAs [23, 24] . miRNA-containing exosomes have been well studied in various animals and humans [25] . However, a few studies have also indicated that piRNAs are detected in exosomes and play key roles in disease and sex development [26] [27] [28] . In molluscs, shell and pearl formation is a complex biomineralization process that involves extensive participation of cells and secreted exosomes [28] . The epidermis, gills and feet of the mantle absorb Ca 2+ and HCO À 3 , which are transported to the mantle by haemocytes and exosomes and secreted by outer mantle epithelial (OME) cells. The shell matrix proteins (SMPs), mostly secreted from OMEs, are directly delivered to the mineralization site via either exosomes or the classical secretory pathway [29, 30] . To date, abundant SMPs, which play a vital role in shell and pearl formation, have been identified [31] [32] [33] . However, the process by which pearl formation occurs remains poorly understood. Currently, most studies on exosomes are primarily focused on their roles in disease and immune response [34, 35] , while exosomes are also involved in shell and pearl formation in molluscs [28] . In pearl farming, small pieces of mantle tissue from donor oysters are implanted into the gonad of host oysters along with an inorganic bead (termed the 'nucleus') for large-scale production of artificial pearls [33, 36] . Surgical implantation of the mantle can induce an immune reaction in response to transplantation for the survival of the host oyster [37] . Therefore, a vast number of immune-related genes are activated upon transplantation [33, [38] [39] [40] . However, the exosome-mediated immunological reaction that occurs in host oysters during the subsequent stages of pearl sac formation is still unclear. Accordingly, an improved understanding of exosome-mediated immune responses in host oysters upon acceptance of a transplant is required to further improve the effectiveness of pearl production. In this study, we aimed to reveal the function of exosome-derived sncRNAs upon mantle grafting using small RNA sequencing (sRNA-seq) profiling. Intriguingly, both miRNAs and piRNAs were discovered in the exosomes of Pinctada fucata, and their target prediction and functional analysis indicated that these small RNAs are related to the immune response upon grafting transplantation. Mantle pieces were sheared from donor oysters and implanted into the gonads of host oysters along with an inorganic bead for pearl culture. The haemolymph was collected from host oysters during grafting transplantation for exosome isolation and purification. Exosome morphologies were detected by transmission electron microscopy (TEM) and size by nanoparticle tracking analysis. As shown in figure 1a, typical spherical structures of 50-200 nm were observed in the purified exosome preparation. Nanoparticle tracking analysis indicated Figure 1 . The exosomes isolated from Pinctada fucata were detected by transmission electron microscopy (a) and nanosight particle tracking analysis (b). Scale bar = 1 µm. royalsocietypublishing.org/journal/rsob Open Biol. 12: 210317 that the concentration of exosomes was 4.72 ± 0.66 × 10 8 particles ml −1 with an average diameter of 169 ± 14.6 nm (figure 1b). All the measurement results suggested that the exosome preparation isolated from P. fucata haemolymph contains a heterogeneous mixture of exosomes and macrovesicles. (d ) Venn diagram of the overlap of known miRNA profiles and conserved miRNA family in P. fucata and related species. Thirty-two known miRNAs were identified from P. fucata exosomes in the present study, 186 known miRNAs from multiple somatic and gonadal tissues of P. fucata [41] , 104 known miRNAs from mantle tissues of P. fucata martensii using the same analysis pipeline [42] and 54 known miRNAs from Crassostrea gigas deposited in MirGeneDB 2.1 [43] . (e) Heatmap of identified miRNA expression patterns in all investigated samples of P. fucata. PA, just after grafting; PB, before grafting; P1D, 1 day post-grafting; P2D, 2 days postgrafting; P3D, 3 days post-grafting; P1W, 1 week post-grafting; P2W, 2 weeks post-grafting. royalsocietypublishing.org/journal/rsob Open Biol. 12: 210317 and known RNAs in Rfam, 37.24 million reads (43.88% of raw reads) remained for miRNA and piRNA identification and characterization. A total of 32 known and six novel miRNAs were identified, with a large variation in expression levels (electronic supplementary material, table S2). The miRNA counts were normalized according to the reads per million reads (RPM) method. The principal component analysis of miRNA expression showed that the miRNA expression patterns at the grafting and initial stages were not substantially different, but varied significantly between the initial and later stages (figure 2c). Among these miRNAs, seven of the 32 known miRNAs, namely miR-2a/b/c, miR-71, miR-317, miR-184 and miR-1990c, were simultaneously discovered in the exosomes and related species, while seven miRNAs, namely miR-745b, miR-1735, miR-3391, miR-4867, miR-6999, miR-9065 and miR-9885, were specifically expressed in P. fucata exosomes (figure 2d). Several immune-related miRNAs were observed in P. fucata exosomes, including miR-2a/b/c and miR-92. Meanwhile, miRNAs, namely miR-750, miR-1990c, miR-1991 and miR-1994, which were especially observed in molluscs, were also detected. Among these miRNAs, miR-67 and miR-750 were upregulated at the initial stage of mantle grafting (P1D, P2D and P3D), whereas miR-87 and novel-1 were downregulated at later stages (P1W and P2W) (figure 2e). Except for miRNAs detected in P. fucata exosomes, abundant piRNA-sized reads 24/25 nt and 29/31 nt in length were also detected by sRNA-seq (figure 3a). All putative reads were merged from all samples and processed for piRNA processing. A total of 277 139 putative unique piRNAs were detected in P. fucata exosomes (electronic supplementary royalsocietypublishing.org/journal/rsob Open Biol. 12: 210317 material, table S3). The total putative piRNA reads of each library were used for normalization to analyse the expression patterns of piRNAs among different samples according to the RPM method. The expression patterns of piRNAs were not substantially different between the grafting and initial stages, but varied significantly between the initial stages and later stages (electronic supplementary material, figure S1). A significant 10 nt overlap between the sense and antisense strands of putative reads was detected in the small RNAs derived from P. fucata exosomes, which is consistent with piRNAs generated by the ping-pong amplification mechanism (figure 3b) with significantly enriched pingpong pairs in each sample (figure 3c). Ping-pong matrices also illustrate that the ping-pong pairs combine piRNAs of homogeneous of 28/31 in length and 28/31 and 24/25 nt in length (figure 3d). Furthermore, the ping-pong pairs were heavily biased for a 5 0 uridine (U) for 28/31 nt piRNAs and stronger bias for an adenine (A) at position 10 for 24/25 nt piRNAs (figure 3e). All the evidence reveals the existence of piRNAs in P. fucata exosomes. To further parse piRNAs in P. fucata exosomes, the putative piRNA reads were merged for clustered locus detection on the genome. A total of 32 piRNA clusters were identified in the genome with a total of 366.94 kb in length. As the background, 539.3 Mb (52.1% of the total genome) of repetitive sequences were annotated in the genome (figure 4a). The identified piRNA clusters showed a lower enrichment for transposon sequences than the whole genome, which contained 39.6% transposons (figure 4b). Intriguingly, the composition of transposons in piRNA clusters does not at all reflect the transposon landscape of the whole genome. Instead, piRNA clusters were enriched for DNA/Crypton, LTR/Ngaro and LTR/ERVK, showing up to sixfold enrichment in piRNA clusters (figure 4c). However, thousands of distinct piRNAs do not cluster on the genome or map to the transposon in P. fucata (figure 4d ). The expression patterns of transposon elements show an opposite tendency of exosome-derived piRNAs, such as a higher expression density of exosome-derived piRNAs in P1D and lower in P1W, as well as a lower expression density of transposon elements in P1D and higher in P1W (figure 4e,f ). Within the 191 086 unigenes from our previous transcriptome analysis of gene expression performance upon grafting transplantation in P. fucata, 8469 and 12 326 unigenes were identified as target genes for all expressed miRNAs and piRNAs, respectively. Enrichment analysis of predicted targets indicated that a vast number of immune-related pathways were significantly ( p < 0.01) enriched, such as the Wnt signalling pathway (ko04310), cGMP-PKG signalling pathway (ko04022), platelet activation (ko04611), MAPK signalling pathway (ko04010) and various disease-related pathways, such as hepatocellular carcinoma (ko05225) and gastric cancer (ko05226) (table 1; electronic supplementary material, table S4 ). The Gene Ontology (GO) enrichment analysis also indicated that various immune-related pathways were enriched in the predicted targets of exosomal miRNAs, which is consistent with the results of the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis (electronic supplementary material, figure S2 ). In particular, the predicted targets of miR-750, which were upregulated at the initial stage of grafting transplantation, were enriched in the immune response (electronic supplementary material, figure S3 ). In addition, the identified piRNAs royalsocietypublishing.org/journal/rsob Open Biol. 12: 210317 were also used for target prediction to understand their potential function in gene regulation. A total of 12 326 unigenes were predicted to be targeted by exosomal piRNAs. Enrichment analysis of these predicted targets also indicated that a vast number of immune-related pathways were significantly ( p < 0.01) enriched, such as the cGMP-PKG signalling pathway (ko04022), calcium signalling pathway (ko04020), Wnt signalling pathway (ko04310), mTOR signalling pathway (ko04150), MAPK signalling pathway (ko04010) and various disease-related pathways (table 1; electronic supplementary material, table S5 ). Lastly, the expression patterns of miR-750 and three piRNAs ( piRNA-8149, piRNA-3419 and piRNA-7882) and their predicted targets were compared. miR-750 was upregulated at the initial stage and downregulated at a later stage upon grafting transplantation, whereas the genes MAF, HTH, GRIA3, APOB, TITIN and CCD39, the miR-750 targets, were downregulated to the initial stage, and upregulated at later stages ( figure 5 ). piR-8149 was immediately upregulated after grafting transplantation and downregulated to initial levels at later stages, whereas the targets of PRDX6, PCD16, RN213, TPPP2, S6A13 and DYHC were downregulated after grafting transplantation and upregulated at later stages. piR-3419 was upregulated at a later stage, and the targets of TBA2, MP, PIF, ZIG8 and THS78 were downregulated at later stages. piR-7882 was immediately upregulated after grafting transplantation and quickly returned to its initial level and downregulated at a later stage; the targets of DAXX, SIN3A and RPC2 showed opposite expression patterns of piR-7882 ( figure 6 ). Among these targets, most are related to the immune response in cellular processes. Nuclear transfer is often used to produce pearls on a large scale; however, the grafting transplantation process causes pearl oyster stress [33] . Exosomes, a cell-to-cell communication tool, are small membrane-enclosed vesicles actively released royalsocietypublishing.org/journal/rsob Open Biol. 12: 210317 by cells into the extracellular environment [44] . The functions of mollusc exosomes are mysterious and diverse, and they are currently known to play an important role in immune response and biomineralization [10] [11] [12] 28] . In this study, we successfully extracted exosomes from the haemolymph of P. fucata, a well-studied mollusc in biomineralization, and further explored the carrying small RNA molecules and their functions during artificial grafting transplantation. TEM showed typical exosome morphology with a spherical shape, and the size of exosomes was analysed using a NanoSight LM10 system with an average size of 170 nm in diameter, which was slightly larger than that of regular exosomes in size, but consistent with extracted exosomes in molluscs of Crassostrea gigas and Hyriopsis cumingii [10] [11] [12] . Exosomes carry cargo molecules that play a pivotal role in cell-to-cell communication, reflecting the physiological or pathological state of donor cells [22, 45] . As one of the most important types of small non-coding molecules in exosomes, miRNAs play an important role in intercellular communication by mediating the repression of critical mRNA targets in neighbouring or distant recipient cells [46] , while the function of specific miRNA in EVs is still debatable in many animals [47] . Small RNA molecules were isolated and sequenced using sRNA-seq from P. fucata exosomes. A total of 38 miRNA species were identified in the examined samples during grafting transplantation, which was fewer than those identified from sRNA-seq from somatic and gonadal tissues and genome-wide prediction [41, 48] . Furthermore, the expression levels of most miRNAs were very low in the exosomes. Among these identified miRNAs, miR-67 and miR-750 increased at the initial stage, and miR-87 and novel-1 were upregulated at the later stage of grafting transplantation. Prediction of miRNA targets is of central importance for a functional understanding of miRNAs in P. fucata exosomes. We performed a computational prediction of the targets of exosomal miRNAs based on the assembly unigenes of our previous transcriptome analysis of grafting transplantation in P. fucata. Thousands of targets were identified based on physical sequence interactions with these miRNAs. The computational prediction of miRNA targets always yields a large proportion of false positives [49] [50] [51] . A complex mRNA-miRNA network was predicted by computational methods in molluscs, which will require a huge amount of experimental work to confirm the predicted targets. However, it does royalsocietypublishing.org/journal/rsob Open Biol. 12: 210317 not affect the understanding of the functional landscape of the whole miRNAs in exosomes. In molluscs, the immune response that occurs during the grafting process plays a vital role in response to host oyster survival and regeneration, and immune-related genes are enriched at the initial stage of mantle grafting transplantation [33] . In this study, functional royalsocietypublishing.org/journal/rsob Open Biol. 12: 210317 annotations of GO and KEGG pathways also indicated that immune response-related pathways were enriched in the predicted target genes. The target prediction of miR-750 also indicated that this miRNA is related to the immune response. In fact, miR-2, an invertebrate-specific family member, plays a vital role in the immune response [52] , and miR-29 is involved in the immune response by regulating neuropeptide Y receptor 2 (Y2R) in Pinctada martensii [53] . In addition, several miRNAs, namely miR-750, miR-1175, miR-1990 and miR-1991, were also observed in P. fucata exosomes, which were specifically discovered in molluscs by genome-wide prediction [48] . Specifically, miR-1990 showed a mantle-specific expression pattern in P. fucata [41] , and miR-1991 was specifically expressed in the adductor muscle of C. gigas [54] and exosome-derived miR-15b negatively regulated hcApo to mediate shell colour formation in H. cumingii [11] . All the evidence indicates that exosomes play an important role in various immune responses, growth, development and biomineralization processes in molluscs. piRNAs are a widespread strategy to effectively suppress transposable element activity to safeguard the genome from detrimental insertion mutagenesis, a common issue in the germline of most animals and higher organisms [55] . Intriguingly, piRNAs have been discovered in the soma of arachnids [56] , cnidarians [57, 58] and molluscs [41, 59, 60] . In the present study, abundant exosome-derived piRNAs with canonical characteristics were first discovered in molluscs; these provide important insights into the evolutionary origins and functional understanding of piRNAs in bilaterians. Approximately 10% of unique piRNAs were clustered into 32 piRNA clusters. In line with the transposon element suppressive role of piRNAs, these piRNA clusters show a lower enrichment for transposon elements than the whole genome. Particularly, the expression level of piRNAs was decreased in P1W, whereas the expression levels of transposons increased accordingly. Exosome-derived piRNAs from the haemolymph may be transmitted to the gonadal tissue and act on the germline cells to suppress transposon activity; however, thousands of distinct piRNAs from exosomes do not map to any transposon sequences, indicating a potential additional target and function of piRNAs in P. fucata. Recently, sufficient evidence has supported the involvement of PIWI/piRNA pathways in protection against invading viruses in invertebrates and vertebrates [61, 62] . In flatfish, piRNAs were identified from skin mucus exosomes to establish an effective immune response to bacterial infection in Cynoglossus semilaevis [63] . Neural stem cell-derived exosomes produce piRNAs that are potentially enriched during the adaptive immunity-like reaction in mice [64] . miRanda, an RNA-RNA interacting prediction tool, was used to predict potential target genes of exosome-derived piRNAs [65] . Tens of thousands of genes were predicted to be targets by exosome-derived piRNAs. Enrichment analysis of these targets implied that exosome-derived piRNAs might also be related to the immune response in P. fucata during grafting transplantation. Despite the extensive studies on molecules and pathways of immune response in molluscs, the repertoire of exosomederived sncRNAs composing the mollusc immunity remains to be revealed. In this study, pathway analysis based on the GO and KEGG databases was performed to characterize the functional consequences of exosome-derived sncRNAs associated with the immune response upon grafting transplantation in P. fucata. Consequently, these immunerelated pathways, which have been comprehensively determined in numerous species, were characterized. These pathways included the cGMP-PKG, MAPK, Wnt and mTOR signalling pathways, which are both enriched by miRNA and piRNA targets, suggesting an important role of these pathways in immune response upon grafting transplantation. The mitogen-activated protein kinase (MAPK) plays a conserved role in the regulation of various molecular and physiological processes, including cell differentiation and proliferation, cell growth and death and immune response [66] . A MAPK homologue, MKK4, was found to be involved in the response to the nucleus insertion operation, indicating its role in host defence mechanisms, potentially in protecting P. fucata from injury caused by grafting or disease [67] . Moreover, the MAPK signalling pathway was also enriched at the early stage of pearl sac and pearl formation via grafting transplantation in P. fucata [33] . In conjunction with the activation of the Ras signalling pathway, MAPK activation induces the expression of multiple genes that regulate the inflammatory response [66] . mTOR is an evolutionarily conserved serine/threonine protein kinase from yeast to mammals that senses and integrates a variety of cellular physiological and environmental signals to regulate cell growth [68] . Dephosphorylation of mTOR contributes to induced autophagy responses to contaminants in marine mussels (Mytilus galloprovincialis) [69] . In addition, the potential target genes related to immune response were screened based on the expression level of transcriptome analysis and trends of miRNA and piRNAs, which could provide a basis for functional studies of these miRNAs and piRNAs in the future. The findings of small RNAs in exosomes contribute to the expanding knowledge base for miRNA and piRNA analysis in molluscs as well as for elucidating the ancestral state of somatic piRNAs in bilaterians. Approximately 2-year-old healthy P. fucata pearl oysters were used as donors and recipients for mantle grafting transplantation. Grafting transplantation was performed by a skilled expert at the Mikimoto Pearl Research Institute, Mie Prefecture, Japan. During the grafting transplantation, we collected seven groups of grafting samples: before grafting (PB), just after grafting (PA), 1 day (P1D), 2 days (P2D), 3 days (P3D), 1 week (P1W) and 2 weeks (P2W) post-grafting. The blood lymph fluid was collected from the adductor muscle of five individuals from each group and mixed for exosome isolation. The mixed haemolymph was centrifuged at 2000g for 30 min to yield the supernatant and stored at −80°C until use. Exosomes were isolated from the stored supernatant using a Total Exosome Isolation (from serum) Kit (Invitrogen, Carlsbad, CA, USA) following the manufacturer's instructions. The purified exosome solution (10 µl) was crossed with a hydrophilic collodion membrane for 15 min and treated with 1% uranyl acetate for 30 s. The grid was then placed royalsocietypublishing.org/journal/rsob Open Biol. 12: 210317 on a TEM sample table for observation and photographed using iTEM. The purified exosome solution was also quantified and sized using the NanoSight LM10 system (Malvern Panalytical, UK) following the manufacturer's instructions. Total RNA was extracted from the resuspended exosomes using the Total Exosome RNA and Protein Isolation Kit (Invitrogen) according to the manufacturer's protocol. The quantity and quality of exosomal total RNA were assessed using a Qubit RNA Assay Kit in a Qubit 3.0 Fluorometer (Life Technologies, Carlsbad, CA, USA) and Agilent 2200 TapeStation (Agilent Technologies, Waldbronn, Germany), respectively. Small RNA-seq libraries were constructed using the SMARTer smRNA-Seq Kit for Illumina (Takara, San Jose, CA, USA) following the manufacturer's recommendations. In brief, total RNA was first polyadenylated to provide a priming sequence for an oligo(dT) primer and then the first-strand cDNA molecules were synthesized by adding a second adapter sequence to the 3 0 end. Full-length Illumina adapters were subsequently linked during polymerase chain reaction amplification. The cDNA libraries were tested using the High Sensitivity D1000 ScreenTape Kit in Agilent 2200 TapeStation (Agilent Technologies). The clean-up products were purified and size selected using AMPure XP beads (Beckman Coulter, CA, USA) to obtain the appropriate size. The post-size-selection library was also validated by Agilent 2200 TapeStation (Agilent Technologies) with an expected size of 180 bp that combined the size of processed small RNAs plus adapters, followed by sequencing on an Illumina HiSeq 2500 platform with a 50 bp single-end module (Macrogen, Tokyo). Three independent small RNA libraries were constructed and sequenced for each time node sample. The raw sequencing data were saved as FASTQ files. We first filtered out the reads with low quality, reads containing a strong poly(A) tail and reads containing poly(N). The 3 0 end adapter sequences were trimmed from sequencing reads using fastx_toolkit (fastx_clipper -a AAAAAAAA -c -l 15 -d 0 -Q 33), followed by size filtration with a length range of 15-35 nt. The clean reads were then mapped with the reference genome sequences [70] and Rfam [71] to remove genome unmapped reads and rRNA, tRNA, snRNA and snoRNA sequences using bowtie (v1.2.1, -f -k 3 -v2). The filtered reads were pooled to perform miRNA prediction with miRBase 22.08 using miRDeep2 (v. 2.0.0.8) using default parameters [72] . The identified miRNAs in exosomes were also compared with previously discovered miRNAs in P. fucata [41, 42] and C. gigas in MirGeneDB 2.1 [43] . The miRNA family was annotated according to a previous study [73] . The total miRNA reads of each library were used for normalization to analyse differentially expressed miRNAs among different samples according to the RPM method. A read found more than five times was considered a potential miRNA in P. fucata exosomes. The remnant reads that were within 24-32 nt in length and did not produce a match to any known ncRNAs were considered putative piRNAs and subsequent piRNA discovery. PPmeter [59] , unitas [74] and proTRAC [75] were used to quantify and compare the amount of ongoing ping-pong amplification, search ping-pong signatures, calculate a Z-score for the enrichment of 10 bp overlaps and predict genomic piRNA clusters, respectively. For repeat annotation, we performed a de novo prediction of repetitive elements in the P. fucata genome using RepeatMasker (v. 4.0.7) [76] based on a de novo repeat library, constructed using RepeatModeler (v. 1.0.11) [77] . Analysis was conducted with the entire repeat dataset as well as with repeats localized in the predicted piRNA clusters. A gene set of transcriptome analysis from our previous study on gene expression performance during grafting transplantation in P. fucata was used for miRNA and piRNA target prediction. The assembly unigenes were aligned NR, Swis-s_Prot, COG, GO and KEGG databases using the Triotate pipeline [78] . The 3 0 UTR sequences were isolated based on the annotation results. Target predictions for miRNAs and piRNAs were performed by miRanda with restricted conditions of a score greater than 150 and a free energy ≤15 kcal mol −1 [79] . Enrichment analysis of predicted targets was performed using the OmicShare tools, a free online platform for omics analysis (http://www.omicshare.com/tools). In the present study, exosomes were successfully extracted from the haemolymph of P. fucata (Mollusca) after transplantation. A vast number of exosome-derived sncRNAs have been identified using sRNA-seq. Intriguingly, in addition to the detection of tens of miRNAs, abundant piRNAs were also discovered in the exosomes. The expression patterns of the miRNAs and piRNAs at the grafting and initial stages were not substantially different, but varied significantly between the initial and later stages. Target prediction and functional analysis indicate that these miRNAs and piRNAs are related to the immune response upon grafting transplantation, whereas piRNAs might also be associated with transposons but are not clear in detail. This work provides the basis for a functional understanding of exosome-derived sncRNAs in molluscs and helps to gain further insight into the PIWI/piRNA pathway function outside of germline cells in invertebrates, which might prove to be applicable to other species in the future. All authors gave final approval for publication and agreed to be held accountable for the work performed herein. Conflict of interest declaration. We declare we have no conflicts of interests with the contents of this article. Funding. This work was supported by the Japan Society for the Pro- royalsocietypublishing.org/journal/rsob Open Biol. 12: 210317 New insights into extracellular vesicle biogenesis and function 2020 Roles of exosomes in ocular diseases Exosomes maintain cellular homeostasis by excreting harmful DNA from cells Emerging roles of extracellular vesicles in cellular senescence and aging The biology, function, and biomedical applications of exosomes Circulating immune cells mediate a systemic RNAi-based adaptive antiviral response in Drosophila Challenges in exosome isolation and analysis in health and disease 2020 Arthropod exosomes as bubbles with message(s) to transmit vectorborne diseases Exosome-mediated apoptosis pathway during WSSV infection in crustacean mud crab Transcriptomic analysis of exosomal shuttle mRNA in Pacific oyster Crassostrea gigas during bacterial stimulation The mantle exosome and microRNAs of Hyriopsis cumingii involved in nacre color formation The mantle exosome proteins of Hyriopsis cumingii participate in shell and nacre color formation Seminal plasma exosomes: promising biomarkers for identification of male and pseudo-males in Cynoglossus semilaevis 2020 Identification and application of PIWIinteracting RNAs from seminal plasma exosomes in Cynoglossus semilaevis Biogenesis of small RNAs in animals MiRNA biogenesis and regulation of diseases: an overview Biogenesis of mammalian microRNAs: a global view PIWIinteracting RNA: its biogenesis and functions Diversity of the piRNA pathway for nonself silencing: worm-specific piRNA biogenesis factors PIWI proteins and PIWI-interacting RNAs in the soma Biology of PIWIinteracting RNAs: new insights into biogenesis and function inside and outside of germlines Specificities of secretion and uptake of exosomes and other extracellular vesicles for cell-to-cell communication Exosomes: from biogenesis and secretion to biological function The emerging role of exosome-derived non-coding RNAs in cancer biology Exosome and exosomal microRNA: trafficking, sorting, and function Exosomal piRNA sequencing reveals differences between heart failure and healthy patients 2020 Identification and comparison of piRNA expression profiles of exosomes derived from human stem cells from the apical papilla and bone marrow mesenchymal stem cells The oyster genome reveals stress adaptation and complexity of shell formation Hemocyte-mediated shell mineralization in the eastern oyster Recent advances of shell matrix proteins and cellular orchestration in marine molluscan shell biomineralization. Front Deep sequencing of ESTs from nacreous and prismatic layer producing tissues and a screen for novel shell formation-related genes in the pearl oyster The diversity of shell matrix proteins: genome-wide investigation of the pearl oyster, Pinctada fucata Gene expression profiles at different stages for formation of pearl sac and pearl in the pearl oyster Pinctada fucata Exosomes and their roles in immune regulation and cancer Exosomes in cancer development, metastasis, and immunity Pearl production Identification of genes potentially related to biomineralization and immunity by transcriptome analysis of pearl sac in pearl oyster Pinctada martensii Cloning and gene expression of signal transducers and activators of transcription (STAT) homologue provide new insights into the immune response and nucleus graft of the pearl oyster Pinctada fucata Deep transcriptome profiling sheds light on key players in nucleus implantation induced immune response in the pearl oyster Pinctada martensii Differentially expressed immune-related genes in hemocytes of the pearl oyster Pinctada fucata against allograft identified by transcriptome analysis Identification and characterization of microRNAs and their predicted functions in biomineralization in the pearl oyster royalsocietypublishing.org/journal/rsob Comprehensive analysis of microRNAs in the mantle central and mantle edge provide insights into shell formation in pearl oyster Pinctada fucata martensii MirGeneDB 2.1: toward a complete sampling of all major animal phyla Exosomes-vesicular carriers for intercellular communication The function and clinical application of extracellular vesicles in innate immune regulation Functional delivery of viral miRNAs via exosomes MicroRNAs are minor constituents of extracellular vesicles that are rarely delivered to target cells 2021 Discovery and functional understanding of miRNAs in molluscs: a genomewide profiling approach Issues in current microRNA target identification methods Too many false targets for microRNAs: challenges and pitfalls in prediction of miRNA targets and their gene ontology in model and non-model organisms Evolution and function of the extended miR-2 microRNA family 2015 miR-29a participated in nacre formation and immune response by targeting Y2R in Pinctada martensii Identification of conserved and novel microRNAs in the Pacific oyster Crassostrea gigas by deep sequencing The piRNA pathway guards the germline genome against transposable elements Pan-arthropod analysis reveals somatic piRNAs as an ancestral defence against transposable elements Characterization of the piRNA pathway during development of the sea anemone Nematostella vectensis Jellyfish genomes reveal distinct homeobox gene clusters and conservation of small RNA processing PIWI genes and piRNAs are ubiquitously expressed in mollusks and show patterns of lineage-specific adaptation Deep sequencing of small RNAs reveals the repertoire of miRNAs and piRNAs in Biomphalaria glabrata PIWI pathway against viruses in insects 2019 piRNA-guided CRISPR-like immunity in eukaryotes Extracellular vesicles Piwi-interacting RNAs from skin mucus for identification of infected Cynoglossus semilaevis with Vibrio harveyi Innate and adaptive immunity of murine neural stem cellderived piRNA exosomes/microvesicles against pseudotyped SARS-CoV-2 and HIV-based lentivirus Pachytene piRNAs instruct massive mRNA elimination during late spermiogenesis The complexity of mitogen-activated protein kinases (MAPKs) made simple Identification and analysis of an MKK4 homologue in response to the nucleus grafting operation and antigens in the pearl oyster, Pinctada fucata 2010 mTOR regulation of autophagy Role of mTOR in autophagic and lysosomal reactions to environmental stressors in molluscs Bivalve-specific gene expansion in the pearl oyster genome: implications of adaptation to a sessile lifestyle Rfam 14: expanded coverage of metagenomic, viral and microRNA families 2012 miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades 2021 Digging into bivalve miRNAomes: between conservation and innovation unitas: the universal tool for annotation of small RNAs 2012 proTRAC: a software for probabilistic piRNA cluster detection, visualization and analysis Using RepeatMasker to identify repetitive elements in genomic sequences RepeatModeler2 for automated genomic discovery of transposable element families A tissue-mapped axolotl de novo transcriptome enables identification of limb regeneration factors Exosome-derived small noncoding RNAs reveal immune response upon grafting transplantation in Pinctada fucata (Mollusca)