Figure 4. miRNAs that arise from snoRNA transcripts. The most thermodynamically stable secondary structures of snoRNA transcripts were generated by mfold [
20]. Annotated miRNA sequences embedded in these transcripts are highlighted in blue (
A–
E). (
A) miR-664a is embedded within the SNORA36B transcript and has been shown to function as a tumor-suppressing RNA in HCC and some types of cervical cancer while also functioning as a tumor-promoting RNA in cervical squamous cell carcinoma [
56]. (
B) miR-1291 is located within the SNORA2C transcript and is characterized as a tumor-suppressing RNA in pancreatic, prostate, breast and renal cell carcinoma [
59,
60,
61,
62,
63]. (
C) miR-1248 is derived from SNORA81 and has been determined to be upregulated in prostate cancer [
64]. (
D) miR-3651 is embedded within SNORA84 and is upregulated in colorectal cancer but downregulated in esophageal cancer [
65,
66]. (
E) miR-768 resides within the SNORD71 transcript and has been shown to be downregulated in lung, breast, ovary, liver, parotid gland, thyroid gland cancers and melanoma [
67,
68,
69]. (
F–
J) Mature miRNA sequences (highlighted in blue) were aligned to their corresponding snoRNA precursor transcripts and snoRNAs were aligned to their respective genomic regions. (
F) Alignment between human genome (GRCh38: chr1:220200526-220200696:-1) (
top), SNORA36B (ENSG00000222370) (
middle), and miR664a (ENSG00000281696) (
bottom). (
G) Alignment between human genome (GRCh38: chr12:48654362-48654538:-1) (
top), SNORA2C (ENSG00000221491) (
middle) and miR-1291 (ENSG00000281842) (
bottom). (
H) Alignment between human genome (GRCh38: chr 3:186786655-186786872:1) (
top), SNORA81 (ENSG00000221420) (
middle), and miR-1248 (ENSG00000283958) (
bottom). (
I) Alignment between human genome (GRCh38: chr 9:92292441-92292613:-1) (
top) SNORA84 (ENSG00000239183) (
middle) and miR-3651 (ENSG00000281156) (
bottom). (
J) Alignment between human genome (GRCh38: chr 16:71758382-71758507:-1) (
top), SNORD71 (ENSG00000223224) (
middle) and miR-768 (ENSG00000223224) (
bottom).
2.5. sd/miR-1248
There is a clear need for additional reliable biomarkers to aid care providers in their decision on which treatment to apply to prostate cancer patients following radical prostatectomy. Aimed at meeting this need, Pudova et al. published a paper in 2020 focused on identifying differentially expressed miRNAs between prostate cancer patient tumors with or without lymphatic dissemination [
64]. Among the nine miRNAs found to be significantly upregulated in prostate cancer with lymphatic dissemination, sd/miR-1248 was identified which arises completely from SNORA81.
2.6. hsa-sno-HBII-296B and hsa-sno-HBII-85-29
In 2015, Müller et al. performed small RNA sequencing on six PDAC tissue samples and five normal pancreas tissue samples to assess the expression of ncRNAs not typically found via microarray analysis [
70]. They considered 45 noncoding RNAs identified as significantly down-regulated in PDAC, of which there were fourteen sdRNAs and a single sno-derived piRNA. The most downregulated sdRNA in PDAC, as determined by log
2 fold change, was hsa-sno-HBII-85-29. Four sdRNAs of a total 78 ncRNAs were significantly upregulated in PDAC. The most upregulated sdRNA in PDAC, as determined by log
2 fold change, was hsa-sno-HBII-296B. Subsequent analyses were focused on a miRNA identified in the study, however this work constitutes the only publication that characterizes differentially expressed sdRNAs focused specifically in PDAC. PDAC is a particularly difficult cancer to detect and treat, with a reported 5-year survival rate of just 11% [
82].
4.7. Sno-miR-28
By measuring snoRNA expression following induced P53 activation, Yu et al. found significant downregulation of all snoRNAs associated with the SNGH1 snoRNA gene [
71]. One sdRNA arising from SNGH1-associated sno28, sno-miR-28, was identified in complex with AGO via HITS-CLIP data and was found to be abundantly expressed in patient breast cancer tissue. Elevating TP53 expression in vitro consequently reduced sno-mir-28, and in silico and in vitro target engagement experiments revealed the P53-stabilizing protein TAF9B as a likely sno-mir28 target. Taken together, the authors proposed that sno-mir28 directly regulates TAF9B to bring about indirect repression of P53, forming a loop as P53 overexpression decreases sno-mir28 levels. Taqman expression analysis of matched breast cancer patient tumor and normal tissues revealed that SNHG1, SNORD28 and sno-miR-28 were all significantly upregulated in tumors. Furthermore, using the MCF10A cell line as a model of undifferentiated breast epithelium, the authors found that sno-miR-28 overexpression enhanced breast cancer cells’ proliferative capacity and colony formation. In sum, this study defined a role for sno-miR-28 as an oncogenic sdRNA heavily involved in suppressing the P53 pathway.
2.8. sdRNA-93
In 2017 our lab published a study in which we investigated the role of sdRNAs in the aggressive breast cancer phenotype [
72]. RNAseq analysis identified 13 total full length snoRNAs differentially expressed (>7.5×) in MDA-MB-231 compared to MCF7, 10 of which consistently give rise to sdRNAs that associate with AGO in publicly available HITS-CLIP data. We elected to focus on the sdRNA arising from sno93 (sdRNA-93) due to it displaying the highest differential expression in MDA-MB-231 (≥75×) and previous publications implicating sdRNA-93 with miRNA-like silencing capabilities in luciferase assays [
83]. SdRNA-93 silencing in MDA-MB-231 reduced cell invasion by >90% while sdRNA-93 overexpression conversely enhanced cell invasion by >100% at the same time point. By contrast, sdRNA-93 silencing in MCF7 had no significant effect on invasion while overexpression resulted in a substantial ~80% increase in cell invasion. We then employed multiple in silico miRNA target-prediction algorithms which predicted with consensus that Pipox, a gene involved in sarcosine metabolism which has been implicated in breast cancer progression, is a likely candidate for regulation by sdRNA-93.
2.9. sdRNA-D19b and sdRNA-A24
A study published in 2022 by our lab identified 38 specifically-excised and differentially-expressed sdRNAs in prostate cancer [
73]. We began by querying PCa patient TCGA datasets alongside TCGA-normal prostate using an in-house web-based search algorithm SURFR (Short Uncharacterized RNA Fragment Recognition). SURFR aligns next generation sequencing (NGS) datasets to a frequently updated database of all human ncRNAs, performs a wavelet analysis to specifically determine the location and expression of ncRNA-derived fragments (ndRNAs) and then conducts an expression analysis to identify significantly differentially expressed ndRNAs. Two sdRNAs, sdRNA-D19b and sdRNA-A24 (
Figure 3), were among the most overexpressed in PCa patient tumors and were identified as AGO-associated in publicly available datasets so they were selected for further scrutiny. In vitro phenotypic assays in PC3 cells, a model cell line for castration-resistant prostate cancer (CRPC), revealed that both sdRNAs markedly increased PC3 cell proliferation and that sdRNA-D19b, in particular, greatly enhanced cell migration. When the sdRNAs were overexpressed alongside treatment with chemotherapeutic agents, both sdRNAs provided drug-specific resistances with sdRNA-D19b levels correlating with paclitaxel resistance and sdRNA-A24 conferring dasatinib resistance. Multiple in silico target-prediction algorithms provided a consensus prediction of the CD44 and CDK12 3′ UTRs as targets for sdRNA-D19b and sdRNA-A24 respectively, which we then confirmed in vitro via the Renilla luciferase assay. Taken together, this work outlines a biologically coherent mechanism by which sdRNAs downregulate tumor suppressors in CRPC to enhance proliferative/metastatic capabilities and to encourage chemotherapeutic resistance.
2.10. Sd/miR-140
Based on prior investigations linking high expression of the long-noncoding RNA MALAT1 with poor prostate cancer patient prognosis, Hao et al.’s 2020 paper focused on elucidating the mechanism of action by which MALAT1 brings about this effect [
41]. In silico target prediction and in vitro confirmation via the Renilla luciferase assay identified miR-140 as a target of MALAT1, indicating that MALAT1 functions as a miRNA sponge to reduce sd/miR-140 bioavailability. This was further supported by RIP-seq analysis where direct engagement of MALAT1 and sd/miR-140 was confirmed. Target prediction indicated that sd/miR-140’s target is the 3′UTR of BIRC, and this was confirmed via the Renilla luciferase assay. Knockdown of MALAT1 inhibited the growth of prostate cancer both in vitro and in vivo, an effect that the authors suggest is brought about through sd/miR-140’s release from repression by MALAT1. This outlines a pathway where sd/miR-140 functions as a tumor-suppressing RNA and is tightly regulated by the lncRNA MALAT1.
2.11. Sd/miR-151
A study in 2020 by Chen et al. found that sd/miR-151 is downregulated in human prostate cancer cell lines [
44]. In prostate cancer cells, Chen et al. demonstrated that sd/miR-151 overexpression inhibited cell proliferation, migration, and invasion; enhanced apoptosis; and sensitized cells to treatment with 5-FU, an antimetabolite chemotherapy. While no target prediction was performed, the authors did find that overexpression of sd/miR-151 reduces phosphorylation of PI3K/AKT. Even without a precisely mapped mechanism of action, this publication links sd/miR-151 with prostate cancer as a tumor-suppressive sdRNA.
2.12. Sd/miR-215
A 2015 study by Ge et al. found that sd/miR-215 expression was reduced in both epithelial ovarian cancer cell lines and patient tissue [
42]. Sd/miR-215 decreased cell proliferation, enhanced apoptosis, and enhanced sensitivity to paclitaxel treatment. In addition, increased sd/miR-205 resulted in a decrease of cellular X-linked inhibitor of apoptosis (XIAP) mRNA albeit without confirmation of target engagement. Vychytilova-Faltejskova et al. found in 2017 that sd/miR-215-5p is significantly reduced in CRC patient samples and that low sd/miR-215-5p expression is associated with late clinical stages of CRC as well as poor overall survival for CRC patients [
43]. Overexpression of sd/miR-215-5p reduced cell proliferation, viability, colony formation, invasion, migration, in vivo tumor volume, and weight while enhancing apoptosis. Epiregulin (EREG) and HOXB9 mRNA were confirmed to be two targets of sd/miR-215-5p. Both genes are involved in epithelial growth factor receptor (EGFR) signaling, a critical pathway exploited by CRC to promote tumor growth.
2.13. Sd/miR-605
With their 2011 publication, Xiao et al. found that sd/mir-605 overexpression significantly reduced cell survival and enhanced apoptosis in breast and colon cancer cell lines [
47]. The MDM2 mRNA was predicted as a target of sd/mir-605, and sd/mir-605 overexpression was found to cause a reduction in MDM2 protein. MDM2 is a regulator of P53 tumor suppressor, and sd/mir-605 transfection was shown to increase P53 activity without increasing P53 protein levels, suggesting the sdRNA relieves MDM2-mediated repression of P53. In 2014 a study by Huang et al. found that single nucleotide polymorphisms (SNPs) in the sd/miR-605 precursor were significantly correlated with a shorter prostate cancer biochemical recurrence for patients [
48].
2.14. Sd/miR-16-1
S/dmiR-16-1 is one of the miRNAs first implicated in human cancer [
1]. In 2002, a paper by Calin et al. found that the sd/miR-16-1 gene was deleted in more than 65% of chronic lymphocytotic leukemia (CLL). Northern blot analysis of sd/miR-16-1 expression in CLL patient cancer and normal tissue revealed that miR-16-1 is downregulated in cancer. Sd/miR-16-1 has since been verified as a tumor suppressor in a multitude of cancers, several of which have been detailed in this 2009 review by Aqeilan et al. [
84]. In 2017, the sdRNA was demonstrated to function as a tumor-suppressing RNA in gastric cancer by Wang et al. [
21]. The authors found sd/miR-16-1 to be downregulated in gastric cancer patient samples. Additionally, overexpression decreased gastric cancer cell migration, invasion, and colony formation in vitro while reducing tumor volume in vivo. Target engagement with the 3′-UTR of TWIST1 was confirmed via the Renilla luciferase assay. Interestingly, a 2018 paper by a separate group, Feng et al., discovered that sd/miR-16-1 also targets TWIST1 in non-small cell lung cancer (NSCLC) [
22]. Overexpression inhibited NSCLC proliferation, migration, and invasion in vitro. Together these publications add to the numerous examples of sd/miR-16-1 functioning as a tumor-suppressing RNA.
2.15. Sd/miR-27b
In 2012, Ishteiwy et al. found that sd/miR-27b is downregulated in castration-resistant prostate tumors compared to primary prostate cancer and normal tissue [
25]. Increasing sd/miR-27b expression in CRPC cell lines reduced invasion, metastasis, and colony formation. While there was no confirmation of target interactions, the authors found that modulating sd/miR-27b expression had a subsequent impact on Rac1 activity and E-cadherin expression. These findings indicate a tumor-suppressive role for sd/miR-27b in CRPC whose mechanism may involve the tumor suppressor gene E-cadherin.
2.16. Sd/miR-31
The diverse role of miR-31 in cancer has been well studied and is the subject of two excellent reviews both in 2013 and more recently in 2018 [
28,
29]. Some examples from the aforementioned reviews include sd/miR-31 acting as a tumor-promoting RNA in CRC, HNSCC, and lung cancer while acting as a tumor-suppressing RNA in glioblastoma, melanoma, and prostate cancer [
30,
31,
32,
33,
74,
75].
2.17. Sd/let-7g
The let-7g miRNA belongs to the let-7 family of miRNAs, all of which have been well studied in the context of human cancer. Reviews published in 2010, 2012, and 2017 detail the let-7 family as well as let-7g specifically [
86,
87,
88]. One example from this includes an early publication from 2008 by Kumar et al. where let-7g was found to suppress the NSCLC phenotype both in vitro and in vivo through the RAS pathway [
34]. Many publications implicating sd/let-7g as a tumor-suppressing RNA in cancer have emerged since the aforementioned review, including a 2019 publication by Chang et al. in colorectal cancer and a 2019 paper by Biamonte et al. in ovarian cancer [
35,
36].
2.18. Sd/miR-28
Sd/miR-28 has been identified as a tumor-suppressing RNA in several cancers, including a 2014 publication by Schneider et al. focused on B-cell lymphoma (BCL) [
37]. The sdRNA was downregulated in patient BCL samples (n = 25) compared to normal B cells (n = 4), and sd/miR-28 overexpression in BCL cell lines suppressed proliferation and colony formation while enhancing apoptosis. Target engagement was confirmed for sd/miR-28 with the 3′-UTRs of MAD2L1, BAG1, RAP1B, and RAB23 mRNA. Taken together, sd/miR-28 suppresses the BCL phenotype by regulating the expression of four oncogenes involved in cell cycle progression and apoptosis.
2.19. pi-sno74, pi-sno75, pi-sno44, pi-sno78, and pi-sno81
Small RNA-seq (smRNAseq) of patient breast cancer tissue identified five differentially expressed sno-derived piRNAs (pi-snos) within the long noncoding RNA (lncRNA) GAS5 locus [
76]. These pi-sno’s included pi-sno74, pi-sno75, pi-sno44, pi-sno78, and pi-sno81. Each of these were found to be downregulated in breast cancer compared to adjacent normal tissue. Microarray and PCR revealed that pi-sno75 upregulates expression of tumor necrosis factor (TNF)-related apoptosis-inducing ligand (TRAIL), a proapoptotic protein. In silico target prediction identified a locus of 169 basepairs (bp) upstream of the TRAIL transcription start site that was highly complementary to pi-sno75 with predicted thermodynamically stable binding. To further investigate this mechanism of action, knockout (KO) and Co-IPs were performed to identify PIWIL1/4 binding partners during pi-sno75 overexpression. This revealed that pi-sno75/PIWIL can specifically interact with the methyltransferase complex protein WDR5, thereby increasing TRAIL expression. In MCF7 cells, pi-sno75 overexpression greatly increased TRAIL expression and, in combination with doxorubicin treatment, enhanced apoptosis. Treatment with pi-sno75 alone resulted in marked reductions of tumor volume in MCF7 and MDA-MB-231 mouse xenografts. In summary, the sno-derived piRNA pi-sno75 recruits epigenetic machinery to specifically upregulate TRAIL, and thus functions as a tumor-suppressive sdRNA.
2.20. pi-sno78 (Sd78-3′)
In 2011, a study by Martens-Uzunova et al. that focused on identifying a miRNA expression signature of prostate cancer also identified significantly differentially expressed sdRNAs [
89]. While the majority of reads aligned to defined microRNAs, sdRNAs were found to be predominantly upregulated in the metastatic LN-PCa compared to local prostate cancer. The authors then focused on canonical miRNAs, leaving these enticing putative sdRNA drivers of metastatic LN-PCa uncharacterized until their follow-up publication in 2015. In their 2015 publication, Martens-Uzunova et al. specifically focused on investigating the role of sdRNAs in prostate cancer progression [
77]. Consistent with their 2011 study, they identified an sdRNA from the 3′ end of snoRNA78 “sd78-3′” that was upregulated in the more aggressive LN-PCa patient samples. By examining expression at multiple stages of PCa, the authors concluded that globally overexpressed sdRNAs, including sd78-3′, are already present at early stages of cancer but exhibit striking overexpression concordantly with malignant transformation. While the authors did not investigate the mechanism of action of sd78-3′, this sdRNA is actually the pi-sno78 characterized in the above study by He et al. [
76].
2.21. piR-017061
Revisiting the 2015 publication by Muller et al. which queried small RNA expression in PDAC patient samples and normal pancreas controls, a total of 123 ncRNAs were found to be significantly differentially expressed [
70]. This approach considered 45 noncoding RNAs identified as significantly down-regulated in PDAC, of which there were fourteen sdRNAs and a single sno-derived piRNA, piR-017061.