To confirm the presence of BDsEVs, three techniques were utilized: western blotting, NTA, and TEM (Figure 3). Western blot results (Figure 3A and Supplementary File 1) show the presence of all five positive markers (CD9, CD63, CD81, Flot-1, and TSG101) and the absence of Calnexin in sEVs (used as a negative control), confirming no contamination with cellular contents. As expected, the brain homogenates (BH) show more protein than observed in the sEVs for the markers tested. TEM images confirmed the morphology of BDsEVs and their relative size within the sample (Figure 3B), which was further confirmed by nanoparticle tracking analysis (Figure 3C).
Figure 3C (upper panel) shows the isolated particle sizes falling between 50-200 nm. From the NTA data, taking the scatter result as 100% of the represented particles, CMO represents around 52.35% of particles (Figure 3C, lower panel) to have a lipid bilayer. Of these CMO-stained particles, 34.66% contained CD9, 15.49% contained CD63, and 12.01% contained CD81. Refer to Supplementary File 2 for NTA control data.
Quantification data of isolated RNA (Table 2) show levels ranging from 3.24-8.76 (ng/µL), values which were normalized when generating the cDNA libraries. TapeStation assay results (Figure 4) show major peaks of around 120-160 bp, with smaller peaks below this of around 50 bp and other smaller peaks that are higher of around 484 bp.
Quality checks for the next-generation sequencing data show very high-quality reads after adapter trimming (Figure 5A). All sequence positions show a Phred score of >30, indicating error rates of less than 1 in 1000 nucleotides. Figure 5B shows the most sequences to be around 21-22 bp, matching the expected size range for microRNAs. Through further downstream mapping (using BowTie2), using the Samtools flagstat analysis biotool, out of 1,304,100 aligned sequences, 1,116 264 (85%) were mapped to the microRNA regions in the human genome (build hg38; Table 3).
Of the total miRNAs, 808 miRNAs were expressed through all 3 samples, with 344 miRNAs shared between all three samples (~43%), 5 shared between control 1 and control 2, 183 shared between control 1 and control 3, and 16 shared between control 2 and control 3 (Figure 6). From this data, the highest expressed miRNA across all three BDsEVs are hsa-let-7b-5p, hsa-miR-143-3p, hsa-miR-30a-5p, hsa-miR-221-3p and hsa-let7i-5p. Equally, the top 5 least differentially expressed miRNA across BDsEVs samples (reads >10) are hsa-miR-128-2-5p, hsa-miR-182-5p, hsa-miR-193b-5p, hsa-miR-448 and hsa-miR-505-3p. All raw and normalized read counts are given in Supplementary File 3.

Figure 1: Tissue processing workflow. The figure shows the stepwise process for the isolation of brain-derived sEVs (BDsEVs) from fresh-frozen brain sections using size exclusion chromatography. Please click here to view a larger version of this figure.

Figure 2: RNA isolation workflow from BDsEVs. The figure shows the stepwise process of isolating total RNA from the BDsEVs, including the pre-treatment step of removing any extracellular RNA using proteinase K and RNase A. Please click here to view a larger version of this figure.

Figure 3: BDsEV characterisation workflow. The figure shows various methods to characterize the BDsEVs. (A) Western blot to detect sEV-specific protein markers. BH Brain homogenate. Calnexin is a negative control for sEVs - absence in sEV and presence in BH confirms the sEV sample purity. For other markers, as expected, the protein quantity is higher in BH than in the sEV samples (20 µg of protein loaded per well). (B) Representative transmission electron microscopy (TEM) image for BDsEVs at 200 µM scale, confirming the expected morphology with a lipid bilayer membrane. (C) Size distribution plot for Nanoparticle Tracking Analysis (NTA) for all particles (upper panel). Each data point represents the variation across three replicates. Normalized proportions for biological particles (sEVs) and relative proportions for specific surface markers (lower panel) (data shown as mean ± SEM of triplicate wells (n = 3)). Please click here to view a larger version of this figure.

Figure 4: Quality check for cDNA library preparation. Tapestation analysis for individual samples after cDNA library preparation. The peaks displayed in each graph display the size of each preparation (horizontal axis, in bp) and the respective intensity of each peak (vertical axis, FU). Expected small RNA library peaks around 120-160 bp, as seen in the figure. Please click here to view a larger version of this figure.

Figure 5: Quality check for miRnome sequencing. (A) FASTQC plot for read quality for the individual samples after adapter trimming. The horizontal axis is for nucleotide position, and the vertical axis is for Phred scores. Box plots for each position show the distribution of the Phred score across all reads for each position. All distributions within the green zone confirm very high-quality data (>Q30) implicating high-quality sample preparation. (B) The sequence length distribution for each sample shows that for all three samples, the majority of the reads were around 21-22 nucleotides - the expected size range for microRNA. Please click here to view a larger version of this figure.

Figure 6: Overlap of MicroRNA expression across all samples. The number of miRNAs detected from each sample (read count ≥ 10) and specific microRNAs shared in more than one sample are shown in the Venn diagram. All raw and normalized reads are shown in Supplementary File 3. Please click here to view a larger version of this figure.
| Antibody Epitope | Dilution Factor of Antibody |
| Mouse Anti-CD9 | 1:2000 |
| Mouse Anti-CD63 | 1:2000 |
| Mouse Anti-CD81 | 1:2000 |
| Mouse Anti-Flotillin 1 | 1:2000 |
| Mouse Anti-TSG101 | 1:1000 |
| Mouse Anti-Calnexin | 1:10 000 |
| Anti-mouse Secondary Antibody | 1:3000 |
Table 1: Antibody epitopes used for the characterization of BDsEVs and their respective dilution factors.
| Sample | miRNA Concentration (ng/μL) |
| 1 | 8.76 |
| 2 | 3.24 |
| 3 | 5.33 |
Table 2: Quantification data of isolated RNA.
| Sample | Aligned Sequences | Mapped miRNA | Percentage (%) |
| Control Sample 1 | 1 444 665 | 1 133 437 | 78.46 |
| Control Sample 2 | 442 808 | 381 631 | 86.18 |
| Control Sample 3 | 2 024 828 | 1 983 726 | 97.97 |
| Average | 1 304 100.33 | 1 166 264.67 | 89.43 |
Table 3: Amount of aligned and successfully mapped miRNA sequences.
Supplementary File 1: Supplementary western blot analysis data. Please click here to download this File.
Supplementary File 2: Supplementary NTA control data. Please click here to download this File.
Supplementary File 3: Raw and normalized read counts of miRNA expression analysis. Please click here to download this File.
Supplementary File 4: Supplementary protein concentration data. Please click here to download this File.