Prediction of the effect of in-seed SNPs in 31 miRNAs on BC development.
Abstract
Breast cancer (BC) is the leading type of cancer among women. Findings have revolutionized current knowledge of microRNA (miRNA) in breast tumorigenesis. The seed region of miRNA regulates the process of gene expression negatively. The presence of SNPs in the seed regions of miRNA dramatically alters the mature miRNA function. Additionally, SNPs in the out-seed region of miRNAs have a significant impact on miRNA targeting. This study focuses on the in silico analysis procedure of mature miRNA SNPs and their impact on BC risk. The database annotated SNPs on mature miRNAs was used. Also, target gene alterations, miRNAs function in BC, and the interaction of miRNAs with targets were predicted. A list of 101 SNPs in 100 miRNAs with functional targets in BC was indicated. Under the SNPs allele variation, 10 miRNAs changed function, 6 miRNAs lost targets, 15 miRNAs gained targets, 48 onco-miRNAs remained unchanged, and 21 tumor suppressor miRNAs remained unchanged. At last, a list of 89 SNPs, which alter miRNA function and miRNA-mRNA interaction, were shown to be potentially associated with BC risk. This research theoretically generated a list of possible causative SNPs in the mature miRNA gene that might be used in future BC management studies.
Keywords
- mature microRNA
- SNP
- breast cancer
- bioinformatics
- in silico
1. Introduction
Breast cancer (BC) is the most prevalent cancer among women across the world. The malignant growth begins in the ducts (85%) or lobules epithelium (15%) (“
MicroRNAs (miRNAs), short non-coding RNAs of approximately 18–24 nucleotides in length, affect gene expression negatively by directly binding the 3′-untranslated region (UTR) of the target messenger RNA (mRNA) and diminishing its stability and translatability. Numerous miRNAs have roles in cell signaling, including proliferation, death, differentiation, and immunity [3]. It has been proposed that they have a role in the formation and development of human malignancies, making them important markers of cancer [4]. The involvement of microRNAs in complex diseases like BC has focused on recent studies [5]. Various miRNAs have been related to regulating genes that have a role in the development of BC. MiR-21, miR-26a, miR-155, miR-221/miR-222, and miR-495 are some of the onco-miRs involved in tumor proliferation and angiogenesis [6, 7, 8, 9, 10, 11, 12]. Various miRNAs (miR-100, miR-125b, miR-126, miR-145, miR-200c, miR-298, and miR-335) have been involved in cell cycle regulation, hypoxia and stress response, and apoptotic induction [13, 14, 15, 16, 17, 18, 19, 20, 21, 22]. In addition, McAnena et al. [23] found that circulating miR-332 and miR-195 may be utilized to discriminate between local and metastatic BC. Meanwhile, Sathipati et al. [24] indicated that 34 miRNAs could be employed to classify the early and late stages of BC development. These findings support the concept that a small selection of miRNAs can be considered biomarkers for BC risk prediction or prognosis [25].
The seed region of miRNA, located between the second and eighth positions in miRNA, regulates gene expression negatively [26]. As miRNAs are small, even a single alteration in the mature sequence of miRNAs might impact the development of miRNAs, which leads to the production of new miRNA [27]. A novel repertoire of target genes is produced by the presence of SNPs in the seed regions of miRNA, which dramatically alters the biological activity of the miRNA [28]. Additionally, SNPs in the out-seed region of miRNAs significantly impact miRNA targeting enhancements [29].
Therefore, in this study, the functional effect of SNPs in miRNAs regulating the genes implicated in BC has been predicted by bioinformatics approaches, including (1) screening miRNAs that reveal SNPs in mature sequences, (2) identifying miRNAs’ target genes related to BC, (3) predicting miRNAs function in BC, and (4) estimating the degree of interaction between miRNAs and target mRNAs. The results may be beneficial in determining potential SNPs to analyze further and investigate SNPs that are causative to BC. Furthermore, the findings may potentially help generate theories and evaluate therapies for BC.
2. Materials and methods
2.1 SNPs in miRNA mature sequences database
In this study, the list of SNPs in mature miRNA sequences was extracted from the “SNPs in pre-miRNAs” database provided by miRNASNP-v3 [30] (http://bioinfo.life.hust.edu.cn/miRNASNP#W/). Mature miRNA location, positioning, and sequences were obtained from miRbase (www.mirbase.org), release 22 [31]. SNP position were acquired from dbSNP [32] (MAF > 1%) (https://www.ncbi.nlm.nih.gov/projects/SNP/snp_summary.cgi).
2.2 In silico target prediction in BC pathways
The miRDB (http://mirdb.org/) was adapted to input both miRNA reference and SNP altered sequences for target prediction [33]. Target prediction was performed for all miRNAs containing reference and SNP alternative sequences. All predicted targets with prediction scores ≥80 that are most likely to be accurate [34], were collected and applied to conduct Kyoto encyclopedia of genes and genomes (KEGG) enrichment analyses for the identified target genes in the BC pathway (hsa05224). KEGG is a knowledge base for systematic analysis of gene functions [35] (https://www.genome.jp/pathway/hsa05224). MiRNAs that target genes in the BC pathway are considered to be involved in BC formation. The target with the highest predicted target score from miRDB is considered the gene most likely to affect BC via the miRNA regulation.
2.3 Function prediction of target genes and miRNAs in BC
The target genes with the highest target score from miRDB of each miRNA were further analyzed. The published expression level (Log2FC value) of genes in BC was collected from GENT2. GENT2 is a platform for exploring Gene Expression patterns across Normal and Tumor tissues [36] (http://gent2.appex.kr). From Log2FC value, target genes were considered as oncogenes (Log2FC > 0, P < 0.05) or tumor suppressors of BC (Log2FC < 0, P < 0.05).
MiRNAs may contribute to (onco-miRNAs) or repress (tumor suppressor miRNAs) the cancer phenotype by inhibiting the expression of tumor suppressor genes or oncogenes, respectively. Therefore, the role of miRNAs in BC can be predicted based on the role of target genes in BC. MiRNAs are considered oncogenic miRNAs when their target genes act as tumor suppressors. Furthermore, conversely, miRNAs were considered tumor-suppressive miRNAs when they target oncogenes. The presence of SNPs on the seed region can alter the target, thereby changing the role of miRNAs in BC.
2.4 Predicting the interaction between miRNAs and targets
The interaction between mature miRNA sequences, in both reference and alternative allele versions of the SNP, and their target genes were analyzed based on the number of seed binding sites of a miRNA at the 3’UTR by TargetScan 8.0 [37] and the Minimum Free Energy (MFE) of the whole mature miRNA sequence with the 3’UTR by RNAhybrid. The target sequences extracted from TargetScan 8.0 (http://www.targetscan.org/vert_80/). The miRNA sequences were further included in RNAhybrid (https://bibiserv.cebitec.uni-bielefeld.de/rnahybrid), a tool for finding the MFE hybridization of a long and a short RNA [38]. The MFE values should be between −10 and − 30 kcal/mole. The more negative the MFE value (closer to −30), the stronger the interaction of the mature miRNA sequence with the 3’UTR. A miRNA can bind at multiple sites on the 3’UTR of its target, and each interaction site has a corresponding MFE value. Therefore, the interaction between the 3’UTR and the whole mature sequence of miRNA was determined via the mean value of MFEs from all miRNA binding sites on 3’UTR (Mean MFE).
The level of miRNA inhibition on target genes could be influenced by the presence of an alternative allele SNP, which could be shown through the difference between the mean MFE value of the alternative and reference allele version (ΔG). A non-zero ΔG value means a difference in the interaction degree with the same target between miRNA reference and alteration, indicating the SNP is potentially associated with BC risk. A ΔG value of 0 indicates no change in the level of interaction between the miRNA reference and the miRNA variation with the same target, suggesting that the SNP is unlikely to relate to BC risk.
3. Results
3.1 SNPs in miRNA mature sequences
This study collected a list of 169 SNPs on mature sequences of 161 miRNAs from miRNASNP-v3. The mature reference sequences of 161 miRNAs were obtained from miRBase. The 169 alternative mature sequences of miRNAs corresponding to positional changes of 169 SNPs identified from dbSNP are described in Supplement 1.
3.2 In silico target prediction in BC pathways
After predicting the target genes of the reference and alternative sequences by miRDB (score ≥ 80) and comparing them in the list of genes on the BC pathway of KEGG, out of 161 miRNAs, there are 100 miRNAs with functional targets in BC (Supplement 2). Thirty-eight miRNAs with 38 in-seed SNPs had altered targets, and 62 with 63 out-seed SNPs did not change targets. The most likely target gene, which has the highest target score, of each miRNA was selected to analyze the possible effect on BC further.
3.3 Function prediction of target genes and miRNAs in BC
The functional target genes in BC were determined from the GENT2 database. Hence, the function of miRNAs in BC was also predicted. Of the 38 miRNAs with altered targets, the targets of 31 miRNAs were functionally altered in BC under the influence of SNPs (Table 1). Of these 31 miRNAs, 4 predicted onco-miRNAs were predicted to change to tumor suppressor miRNAs because of their new target working as oncogenes (hsa-miR-4257, hsa-miR-499a-3p, hsa-miR-501-3p, hsa-miR-593-5p), 5 predicted onco-miRNAs were predicted to lost targets in BC (hsa-miR-1178-5p, hsa-miR-4482-5p, hsa-miR-4661-3p, hsa-miR-5589-3p, hsa-miR-7854-3p), 6 predicted tumor suppressor miRNAs were predicted to change to onco-miRNAs because of their new target working as tumor suppressor genes (hsa-miR-1302, hsa-miR-2682-3p, hsa-miR-3117-3p, hsa-miR-4695-5p, hsa-miR-5692b, hsa-miR-6810-5p), 1 tumor suppressor miRNA was predicted to lost targets in BC (hsa-miR-627-5p), 10 miRNAs were predicted to gain onco-miRNAs (hsa-miR-4513, hsa-miR-4707-3p, hsa-miR-4741, hsa-miR-4781-3p, hsa-miR-4804-5p, hsa-miR-5090, hsa-miR-662, hsa-miR-6796-3p, hsa-miR-6826-5p, hsa-miR-6879-3p), and 5 miRNAs were predicted to gain tumor suppressor miRNAs (hsa-miR-1269b, hsa-miR-4467, hsa-miR-6717-5p, hsa-miR-6763-3p, hsa-miR-6777-5p). By altering the function of miRNAs, 31 in-seed SNPs on these miRNAs are likely to be associated with BC risk (Table 1).
miRNA | SNP | Reference | Alternative | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
ID | Position | Allele | Target | Role of miRNA | Allele | Target | Role of miRNA | |||||||
Gene | Log2FC | Expression trend | Category | Gene | Log2FC | Expression trend | Category | |||||||
hsa-miR-1178-5p | rs7311975 | chr12:119713688 | T | FGF9 | −1.22 | Down | Tumor suppressor genes | Onco-miRNA | C | |||||
hsa-miR-1269b | rs7210937 | chr17:12917329 | G | C | DVL3 | 0.28 | Up | Oncogenes | Tumor suppressor miRNA | |||||
hsa-miR-1302 | rs74647838 | chr12:112695096 | G | ESR1 | 0.95 | Up | Oncogenes | Tumor suppressor miRNA | A | CTNNB1 | −0.40 | Down | Tumor suppressor genes | Onco-miRNA |
hsa-miR-2682-3p | rs74904371 | chr1:98045291 | C | TNFSF11 | 1.80 | Up | Oncogenes | Tumor suppressor miRNA | T | SOS1 | −0.43 | Down | Tumor suppressor genes | Onco-miRNA |
hsa-miR-3117-3p | rs12402181 | chr1:66628488 | G | KRAS | 0.31 | Up | Oncogenes | Tumor suppressor miRNA | A | SOS2 | −0.37 | Down | Tumor suppressor genes | Onco-miRNA |
hsa-miR-4257 | rs74743733 | chr1:150551992 | G | MAPK1 | −0.09 | Down | Tumor suppressor genes | Onco-miRNA | A | E2F2 | 1.20 | Up | Oncogenes | Tumor suppressor miRNA |
hsa-miR-4467 | rs115101071 | chr7:102471476 | G | A | IGF1R | 0.41 | Up | Oncogenes | Tumor suppressor miRNA | |||||
hsa-miR-4482-5p | rs45596840 | chr10:104268396 | G | TCF7 | −0.47 | Down | Tumor suppressor genes | Onco-miRNA | A | |||||
hsa-miR-4513 | rs2168518 | chr15:74788737 | G | A | NCOA1 | −0.21 | Down | Tumor suppressor genes | Onco-miRNA | |||||
hsa-miR-4661-3p | rs12335005 | chr8:91205534 | G | NOTCH2 | −0.60 | Down | Tumor suppressor genes | Onco-miRNA | T | |||||
hsa-miR-4695-5p | rs79637190 | chr1:18883265 | C | SHC1 | 0.21 | Up | Oncogenes | Tumor suppressor miRNA | T | MAPK1 | −0.09 | Down | Tumor suppressor genes | Onco-miRNA |
hsa-miR-4707-3p | rs2273626 | chr14:22956973 | C | A | FGF9 | −1.22 | Down | Tumor suppressor genes | Onco-miRNA | |||||
hsa-miR-4741 | rs7227168 | chr18:22933411 | C | T | FZD8 | −0.70 | Down | Tumor suppressor genes | Onco-miRNA | |||||
hsa-miR-4781-3p | rs74085143 | chr1:54054127 | G | A | FGF2 | −1.44 | Down | Tumor suppressor genes | Onco-miRNA | |||||
hsa-miR-4804-5p | rs266435 | chr5:72878605 | C | G | WNT2B | −0.96 | Down | Tumor suppressor genes | Onco-miRNA | |||||
hsa-miR-499a-3p | rs3746444 | chr20:34990448 | A | TCF7L2 | −0.68 | Down | Tumor suppressor genes | Onco-miRNA | G | WNT4 | 0.38 | Up | Oncogenes | Tumor suppressor miRNA |
hsa-miR-501-3p | rs149912461 | chrX:50009773 | A | HEY1 | −0.39 | Down | Tumor suppressor genes | Onco-miRNA | G | FGF4 | 0.30 | Up | Oncogenes | Tumor suppressor miRNA |
hsa-miR-5090 | rs3823658 | chr7:102465754 | G | A | FZD4 | −1.31 | Down | Tumor suppressor genes | Onco-miRNA | |||||
hsa-miR-5589-3p | rs116796353 | chr19:10038396 | A | MAP2K1 | −0.16 | Down | Tumor suppressor genes | Onco-miRNA | G | |||||
hsa-miR-5692b | rs451887 | chr21:42951004 | T | FGF18 | 0.34 | Up | Oncogenes | Tumor suppressor miRNA | C | EGF | −0.71 | Down | Tumor suppressor genes | Onco-miRNA |
hsa-miR-593-5p | rs73721294 | chr7:128081882 | C | SP1 | −0.29 | Down | Tumor suppressor genes | Onco-miRNA | T | KRAS | 0.31 | Up | Oncogenes | Tumor suppressor miRNA |
hsa-miR-627-5p | rs2620381 | chr15:42199650 | A | DVL3 | 0.28 | Up | Oncogenes | Tumor suppressor miRNA | C | |||||
hsa-miR-662 | rs9745376 | chr16:770249 | G | A | IGF1 | −1.13 | Down | Tumor suppressor genes | Onco-miRNA | |||||
hsa-miR-6717-5p | rs117650137 | chr14:21023373 | G | A | SHC1 | 0.21 | Up | Oncogenes | Tumor suppressor miRNA | |||||
hsa-miR-6763-3p | rs3751304 | chr12:132582046 | C | T | WNT4 | 0.38 | Up | Oncogenes | Tumor suppressor miRNA | |||||
hsa-miR-6777-5p | rs56155608 | chr17:17813539 | G | A | CSNK1A1 | 0.13 | Up | Oncogenes | Tumor suppressor miRNA | |||||
hsa-miR-6796-3p | rs3745198 | chr19:40369893 | C | G | CDK6 | −0.67 | Down | Tumor suppressor genes | Onco-miRNA | |||||
hsa-miR-6810-5p | rs62182086 | chr2:218341922 | A | CSNK1A1 | 0.13 | Up | Oncogenes | Tumor suppressor miRNA | G | FGF19 | −1.16 | Down | Tumor suppressor genes | Onco-miRNA |
hsa-miR-6826-5p | rs6771809 | chr3:129272155 | T | C | MTOR | −0.22 | Down | Tumor suppressor genes | Onco-miRNA | |||||
hsa-miR-6879-3p | rs74814065 | chr11:65018557 | C | T | FGF9 | −1.22 | Down | Tumor suppressor genes | Onco-miRNA | |||||
hsa-miR-7854-3p | rs2925980 | chr16:81533949 | A | AKT3 | −0.41 | Down | Tumor suppressor genes | Onco-miRNA | G |
The remaining 7 of 38 miRNAs with altered targets were predicted to remain functional despite the allele change of the SNPs (Table 2). Four (hsa-miR-146a-3p, hsa-miR-4284, hsa-miR-4731-3p, hsa-miR-557) and three miRNAs (hsa-miR-3622a-5p, hsa-miR-449c-3p, hsa-miR-548 t-3p) were predicted to be onco-miRNAs and tumor suppressor miRNAs, respectively.
MiRNA | ID | Position | Reference | Alternative | ΔG | ||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Allele | Target | Role of miRNA | Mean MFE | Allele | Target | Role of miRNA | Mean MFE | ||||||||||
Gene | Log2FC | Expression trend | Category | Gene | Log2FC | Expression trend | Category | ||||||||||
has-miR-146a-3p | rs2910164 | chr5:160485411 | C | FZD4 | −1.31 | Down | Tumor suppressor genes | Onco-miRNA | −16.55 | G | WNT9B | −0.40 | Down | Tumor suppressor genes | Onco-miRNA | −18.2 | −1.7 |
hsa-miR-3622a-5p | rs66683138 | chr8:27701697 | G | E2F3 | 0.25 | Up | Oncogenes | Tumor suppressor miRNA | −22.6 | A | SHC1 | 0.21 | Up | Oncogenes | Tumor suppressor miRNA | −19.4 | 3.2 |
hsa-miR-4284 | rs11973069 | chr7:73711334 | C | FZD4 | −1.31 | Down | Tumor suppressor genes | Onco-miRNA | −23.7 | T | DLL1 | −0.64 | Down | Tumor suppressor genes | Onco-miRNA | −20.5 | 3.2 |
hsa-miR-449c-3p | rs35770269 | chr5:55172296 | A | KRAS | 0.31 | Up | Oncogenes | Tumor suppressor miRNA | −18.2 | T | AKT2 | 0.05 | Up | Oncogenes | Tumor suppressor miRNA | −17.6 | 0.6 |
hsa-miR-4731-3p | rs66507245 | chr17:15251649 | T | NOTCH2 | −0.60 | Down | Tumor suppressor genes | Onco-miRNA | −14.7 | A | CTNNB1 | −0.39 | Down | Tumor suppressor genes | Onco-miRNA | −19.7 | −5.0 |
hsa-miR-548 t-3p | rs73872515 | chr4:173268209 | A | NRAS | 0.21 | Up | Oncogenes | Tumor suppressor miRNA | −15.2 | C | IGF1R | 0.41 | Up | Oncogenes | Tumor suppressor miRNA | −16.7 | −1.5 |
hsa-miR-557 | rs78825966 | chr1:168375591 | C | RPS6KB1 | −0.05 | Down | Tumor suppressor genes | Onco-miRNA | −18.8 | T | PGR | −0.85 | Down | Tumor suppressor genes | Onco-miRNA | −18.5 | 0.3 |
For the 62 miRNAs with 63 out-seed SNPs, the targets were not changed because the seed sequence was preserved between the two alleles of the SNPs (Table 3). The result showed that 44 and 18 miRNAs were predicted to function as oncogene and tumor suppressors, respectively. The predictive results do not sufficiently demonstrate an association of these 70 SNPs with BC risk. Therefore, the next step is further to predict the interaction of 69 miRNAs and targets to elucidate the potential association of these 70 SNPs on BC risk.
MiRNA | ID | Position | Allele Ref./Alt. | Target | Role of miRNA | Mean MFE | ΔG | ||||
---|---|---|---|---|---|---|---|---|---|---|---|
Gene | Log2FC | Expression trend | Category | Ref. | Alt. | ||||||
hsa-miR-449b-5p | rs10061133 | chr5:55170716 | A/G | DLL1 | −0.64 | Down | Tumor suppressor genes | Onco-miRNA | −25.8 | −24.6 | 1.2 |
hsa-miR-6801-3p | rs10412196 | chr19:52222085 | T/C | AKT3 | −0.41 | Down | Tumor suppressor genes | Onco-miRNA | −20.8 | −19.8 | 1.1 |
hsa-miR-548ae-5p | rs10461441 | chr5:58530093 | A/G | CDK6 | −0.67 | Down | Tumor suppressor genes | Onco-miRNA | −14.2 | −12.9 | 1.3 |
hsa-miR-4700-3p | rs1055070 | chr12:120723245 | T/G | NCOA3 | −0.01 | Down | Tumor suppressor genes | Onco-miRNA | −18.1 | −18.6 | −0.5 |
hsa-miR-4302 | rs11048315 | chr12:25874055 | G/A | CDK6 | −0.67 | Down | Tumor suppressor genes | Onco-miRNA | −19.4 | −19.0 | 0.4 |
hsa-miR-5579-3p | rs11237828 | chr11:79422176 | T/C | JAG1 | −0.21 | Down | Tumor suppressor genes | Onco-miRNA | −18.6 | −18.6 | 0.0 |
hsa-miR-501-3p | rs112489955 | chrX:50009781 | G/A | CDK6 | −0.67 | Down | Tumor suppressor genes | Onco-miRNA | −18.3 | −18.3 | 0.0 |
hsa-miR-5691 | rs112511786 | chr11:9090358 | G/C | MAPK1 | −0.09 | Down | Tumor suppressor genes | Onco-miRNA | −19.1 | −19.7 | −0.6 |
hsa-miR-4786-5p | rs115063401 | chr2:239943064 | G/A | AKT3 | −0.41 | Down | Tumor suppressor genes | Onco-miRNA | −23.7 | −23.7 | 0.0 |
hsa-miR-580-5p | rs115089112 | chr5:36147955 | T/C | FZD3 | 0.31 | Up | Oncogenes | Tumor suppressor miRNA | −13.4 | −13.9 | −0.5 |
hsa-miR-3124-3p | rs115160731 | chr1:248826432 | C/A | PGR | −0.86 | Down | Tumor suppressor genes | Onco-miRNA | −16.0 | −17.0 | −1.0 |
hsa-miR-942-3p | rs115372145 | chr1:117094703 | C/T | SOS2 | −0.37 | Down | Tumor suppressor genes | Onco-miRNA | −18.3 | −18.3 | 0.0 |
hsa-miR-4457 | rs115769169 | chr5:1309319 | C/T | APC | −0.22 | Down | Tumor suppressor genes | Onco-miRNA | −16.7 | −16.1 | 0.6 |
hsa-miR-3130-3p | rs115772313 | chr2:206783285 | G/A | E2F3 | 0.25 | Up | Oncogenes | Tumor suppressor miRNA | −17.3 | −17.3 | 0.0 |
hsa-miR-4514 | rs116034786 | chr15:80997457 | A/G | ESR1 | 0.95 | Up | Oncogenes | Tumor suppressor miRNA | −22.1 | −20.4 | 1.7 |
hsa-miR-196a-3p | rs11614913 | chr12:53991815 | C/T | JAG1 | −0.21 | Down | Tumor suppressor genes | Onco-miRNA | −19.9 | −18.9 | 1.0 |
hsa-miR-548at-5p | rs11651671 | chr17:42494785 | G/A | PGR | −0.86 | Down | Tumor suppressor genes | Onco-miRNA | −15.3 | −13.8 | 1.5 |
hsa-miR-3192-3p | rs11907020 | chr20:18470681 | T/C | NCOA3 | −0.01 | Down | Tumor suppressor genes | Onco-miRNA | −14.0 | −14.0 | 0.0 |
hsa-miR-5700 | rs12314280 | chr12:94561809 | T/C | AXIN2 | −0.51 | Down | Tumor suppressor genes | Onco-miRNA | −12.9 | −14.4 | −1.5 |
hsa-miR-4433a-5p | rs12473206 | chr2:64340782 | C/G | KRAS | 0.31 | Up | Oncogenes | Tumor suppressor miRNA | −17.0 | −16.6 | 0.4 |
hsa-miR-6077 | rs1280409926 | chr1:148388307 | C/T | SP1 | −0.29 | Down | Tumor suppressor genes | Onco-miRNA | −21.7 | −21.9 | −0.2 |
hsa-miR-548 l | rs13447640 | chr11:94466555 | G/A | SP1 | −0.29 | Down | Tumor suppressor genes | Onco-miRNA | −14.5 | −15.9 | −1.5 |
hsa-miR-548ar-3p | rs141659366 | chr13:114244551 | G/A | FZD4 | −1.31 | Down | Tumor suppressor genes | Onco-miRNA | −14.4 | −12.9 | 1.5 |
hsa-miR-4444 | rs142357696 | chr3:75214531 | A/G | E2F1 | 1.18 | Up | Oncogenes | Tumor suppressor miRNA | −24.7 | −22.1 | 2.6 |
hsa-miR-888-5p | rs143634721 | chrX:145994837 | C/A | PGR | −0.86 | Down | Tumor suppressor genes | Onco-miRNA | −16.2 | −15.3 | 0.9 |
hsa-miR-892b | rs146806052 | chrX:145997215 | A/G | DLL1 | −0.64 | Down | Tumor suppressor genes | Onco-miRNA | −21.0 | −24.2 | −3.2 |
hsa-miR-8060 | rs1514422 | chr3:96360020 | G/A | KRAS | 0.31 | Up | Oncogenes | Tumor suppressor miRNA | −18.3 | −18.5 | −0.2 |
hsa-miR-6887-5p | rs1688017 | chr19:35122719 | G/A | PIK3R3 | 0.33 | Up | Oncogenes | Tumor suppressor miRNA | −28.7 | −25.2 | 3.5 |
hsa-miR-1304-3p | rs2155248 | chr11:93733700 | G/T | NOTCH2 | −0.60 | Down | Tumor suppressor genes | Onco-miRNA | −18.3 | −18.6 | −0.3 |
hsa-miR-3130-5p | rs2241347 | chr2:206783257 | C/T | TCF7 | −0.47 | Down | Tumor suppressor genes | Onco-miRNA | −22.3 | −22.6 | −0.4 |
hsa-miR-1343-5p | rs2986407 | chr11:34941869 | T/C | CSNK1A1 | 0.13 | Up | Oncogenes | Tumor suppressor miRNA | −21.7 | −20.9 | 0.8 |
hsa-miR-5087 | rs2992458 | chr1:148334528 | A/G | GSK3B | 0.14 | Up | Oncogenes | Tumor suppressor miRNA | −17.2 | −16.2 | 1.0 |
hsa-miR-5189-3p | rs35613341 | chr16:88468999 | C/G | JAG1 | −0.21 | Down | Tumor suppressor genes | Onco-miRNA | −25.6 | −28.3 | −2.7 |
hsa-miR-8084 | rs404337 | chr8:93029770 | G/A | NOTCH1 | −0.28 | Down | Tumor suppressor genes | Onco-miRNA | −12.7 | −13.4 | −0.7 |
hsa-miR-548ap-5p | rs4414449 | chr15:85825667 | G/A | CDK6 | −0.67 | Down | Tumor suppressor genes | Onco-miRNA | −14.0 | −13.0 | 0.9 |
hsa-miR-608 | rs4919510 | chr10:100975021 | C/G | CDKN1A | −0.36 | Down | Tumor suppressor genes | Onco-miRNA | −30.7 | −31.2 | −0.5 |
hsa-miR-7157-3p | rs56148568 | chr2:140586631 | T/C | EGFR | −1.44 | Down | Tumor suppressor genes | Onco-miRNA | −19.1 | −19.1 | 0.0 |
hsa-miR-6744-3p | rs56310773 | chr11:1256664 | C/T | PIK3R1 | −0.81 | Down | Tumor suppressor genes | Onco-miRNA | −32.5 | −30.0 | 2.5 |
hsa-miR-6805-3p | rs56312243 | chr19:55388234 | C/T | MAPK1 | −0.09 | Down | Tumor suppressor genes | Onco-miRNA | −19.5 | −19.7 | −0.3 |
hsa-miR-6071 | rs56790095 | chr2:85783659 | C/G | WNT3 | 0.34 | Up | Oncogenes | Tumor suppressor miRNA | −17.0 | −17.0 | 0.0 |
hsa-miR-548ab | rs59323834 | chr3:103524093 | C/T | CDK6 | −0.67 | Down | Tumor suppressor genes | Onco-miRNA | −14.6 | −12.1 | 2.5 |
hsa-miR-3928-5p | rs5997893 | chr22:31160117 | A/G | ESR1 | 0.95 | Up | Oncogenes | Tumor suppressor miRNA | −17.7 | −15.6 | 2.1 |
hsa-miR-596 | rs61388742 | chr8:1817259 | T/C | PTEN | −0.26 | Down | Tumor suppressor genes | Onco-miRNA | −19.7 | −17.5 | 2.2 |
hsa-miR-3922-5p | rs61938575 | chr12:104591665 | G/A | SHC1 | 0.21 | Up | Oncogenes | Tumor suppressor miRNA | −28.6 | −27.3 | 1.4 |
hsa-miR-4772-5p | rs62154973 | chr2:102432320 | C/T | FOS | −1.17 | Down | Tumor suppressor genes | Onco-miRNA | −20.3 | −20.3 | 0.0 |
hsa-miR-646 | rs6513497 | chr20:60308547 | T/G | PIK3R1 | −0.81 | Down | Tumor suppressor genes | Onco-miRNA | −23.4 | −20.9 | 2.6 |
hsa-miR-8063 | rs7162033 | chr15:36972836 | C/G | IGF1 | −1.13 | Down | Tumor suppressor genes | Onco-miRNA | −17.6 | −15.9 | 1.7 |
hsa-miR-8063 | rs7183051 | chr15:36972838 | G/A | IGF1 | −1.13 | Down | Tumor suppressor genes | Onco-miRNA | −17.6 | −16.9 | 0.7 |
hsa-miR-6868-3p | rs7208391 | chr17:76098024 | C/G | AXIN2 | −0.51 | Down | Tumor suppressor genes | Onco-miRNA | −18.8 | −18.6 | 0.2 |
hsa-miR-4679 | rs72810954 | chr10:89063382 | G/A | WNT16 | 0.12 | Up | Oncogenes | Tumor suppressor miRNA | −17.0 | −17.0 | 0.0 |
hsa-miR-4799-5p | rs72955519 | chr4:147782619 | G/A | MAPK1 | −0.09 | Down | Tumor suppressor genes | Onco-miRNA | −18.8 | −16.3 | 2.5 |
hsa-miR-4999-5p | rs72996752 | chr19:8389352 | A/G | CTNNB1 | −0.40 | Down | Tumor suppressor genes | Onco-miRNA | −18.5 | −17.4 | 1.1 |
hsa-miR-624-3p | rs73251987 | chr14:31014677 | C/G | GSK3B | 0.14 | Up | Oncogenes | Tumor suppressor miRNA | −15.5 | −13.7 | 1.9 |
hsa-miR-4727-5p | rs73295187 | chr17:38825855 | A/C | CDK6 | −0.67 | Down | Tumor suppressor genes | Onco-miRNA | −21.4 | −19.6 | 1.8 |
hsa-miR-4739 | rs73410309 | chr17:79707227 | G/C | BAK1 | 0.49 | Up | Oncogenes | Tumor suppressor miRNA | −32.5 | −31.6 | 0.9 |
hsa-miR-6504-5p | rs74469188 | chr16:81611365 | T/C | FZD5 | −0.56 | Down | Tumor suppressor genes | Onco-miRNA | −21.9 | −21.5 | 0.4 |
hsa-miR-323b-5p | rs75330474 | chr14:101056252 | C/T | PIK3R3 | 0.33 | Up | Oncogenes | Tumor suppressor miRNA | −18.5 | −16.5 | 2.0 |
hsa-miR-6841-3p | rs76347846 | chr8:24953808 | A/G | GADD45A | −0.30 | Down | Tumor suppressor genes | Onco-miRNA | −19.1 | −15.2 | 3.9 |
hsa-miR-4704-3p | rs76595065 | chr13:66218307 | T/C | GSK3B | 0.14 | Up | Oncogenes | Tumor suppressor miRNA | −23.2 | −23.4 | −0.2 |
hsa-miR-6839-5p | rs7804972 | chr7:64679085 | G/A | CDK6 | −0.67 | Down | Tumor suppressor genes | Onco-miRNA | −19.4 | −18.1 | 1.3 |
hsa-miR-6885-5p | rs78293125 | chr19:6389688 | A/G | TCF7L2 | −0.68 | Down | Tumor suppressor genes | Onco-miRNA | −24.4 | −24.4 | 0.0 |
hsa-miR-4520-3p | rs8078913 | chr17:6655449 | C/G | GRB2 | 0.35 | Up | Oncogenes | Tumor suppressor miRNA | −24.2 | −21.2 | 3.0 |
hsa-miR-548 h-5p | rs9913045 | chr17:13543607 | G/A | CDK6 | −0.67 | Down | Tumor suppressor genes | Onco-miRNA | −13.8 | −13.8 | 0.0 |
3.4 Predicting the interaction between miRNAs and targets
Among the 69 miRNAs, the interaction of 12 miRNAs with their targets did not differ under the allele variation of the SNPs (ΔG = 0) (Table 3), indicating 12 SNPs in these 12 miRNAs are not likely to be associated with BC risk. The remaining 57 miRNAs have different levels of interaction with their targets under the allele variation of the SNPs (ΔG ≠ 0) (Tables 2 and 3). It suggested that 58 SNPs in these miRNAs are potentially associated with BC risk (Tables 2 and 3).
4. Discussion
In knowledge and improvements in bioinformatics, computational predictions of causative factors are being used as a supplementary approach to support the practical assessment of multifactorial disorders. Although miRNA can modulate up to 92% of mammalian genes, only a few target pairings of miRNAs have been experimentally verified [39]. Numerous technical challenges, such as tissue selectivity, poor expression, 3’ UTR selection, and miRNA preservation, make existing methodologies difficult for experimental confirmation of interactions between miRNAs and their mRNA targets [40]. Recognizing functional SNPs in genes and examining their impacts on phenotypes may allow for more in-depth knowledge of the possible consequences of making such changes. Biogenesis, expression level, and biological function are all influenced by SNPs in human miRNA genes. For identifying the possible impacts of SNPs, several useful bioinformatics tools have been developed. All mature microRNAs with SNPs implicated in BC and their target genes were obtained. In addition, bioinformatics strategies for predicting these functional SNPs were presented. All steps and results are summarized in Figure 1.
The results showed that 12 SNPs did not appear to be associated with BC risk. Thirty-one in-seed SNPs are likely to be strongly associated with disease risk by altering miRNA function in BC (Table 1). 58 SNPs appear to be moderately or weakly associated with BC risk by altering the degree of interaction between miRNAs and their targets (Tables 2 and 3). According to the most significant difference between Mean MFE values in the hsa-miR-4731-3p, the effect of T/A replacement can significantly increase the change of interaction (ΔG = −5); thus, it increases the effect of SNP rs66507245 on the risk of BC (Table 2). It seems that the association of rs66507245 with BC risk was moderate among the group of 70 SNPs predicted based on the ΔG value. The most negligible difference occurs in hsa-miR-6077, hsa-miR-8060, hsa-miR-4704-3p, and hsa-miR-6868-3p, indicating the effect of the SNPs on BC risk were weak (Table 3).
These SNPs were not found in GWASs since GWAS is a whole-genome sequencing approach that does not determine SNPs in non-coding regions. Only four SNPs (rs3746444, rs2910164, rs11614913, and rs4919510) have been investigated for the association with BC risk in case–control studies. The presence of the C allele of rs2910164, an in-seed SNP, in hsa-miR-146a-3p was associated with an increased risk of BC (C vs. G: OR = 1.4, 95% CI = 1.03–1.85, p = 0.03) [41]. Allele C of rs11614913 in hsa-miR-196a-3p was found to be significantly associated with decreased risk of BC (C vs. T: OR = 0.64, 95%CI = 0.49–0.85, p = 0.0019) [42]. Allele G of rs3746444 in seed region of hsa-miR-499a-3p significantly increased BC susceptibility only among Asians (G vs. A: OR = 1.12, 95% CI = 1.00–1.26, P = 0.04) [43]. The G allele rs4919510 in hsa-miR-608 decreased the risk of BC (G vs. C: OR = 0.53, 95%CI 0.30–0.92, p = 0.024) [44]. The remaining 97 SNPs have not been identified for an association with BC risk in any association studies. However, 7 of these (rs7210937, rs12402181, rs2273626, rs2620381, rs35770269, rs10061133, and rs6513497) are associated with the risk of other cancers, including oral and pharyngeal squamous cell carcinoma [45], acute lymphoblastic leukemia [46], leukopenia [47], gastric cancer [48], colon cancer [49], esophageal squamous cell carcinoma [50], and hepatocellular carcinoma [51]. These shreds of evidence suggest that candidate SNPs in this study can be associated with risk in BC and other types of cancer.
At the level of investigating the role of miRNAs in BC, 31 miRNAs and 69 miRNAs were predicted to change and remain unchanged function in the presence of SNP allele changes, respectively. Among them, 7 miRNAs (miR-4513, miR-501-3p, miR-580-5p, miR-3130-3p, miR-196a-3p, miR-8084, and miR-3922-5p) had predicted roles in BC consistent with the results of previous functional studies. Hsa-miR-4513, has-miR-196a-3p, and has-miR-8084 suggested functions as an oncogene in the progression of BC [52, 53, 54] were also predicted as onco-miRNAs in this study (Tables 1 and 3). Hsa-miR-501-3p, has-miR-580-5p, has-miR-3130-3p, and has-miR-3922-5p indicated tumor suppressor function [55, 56, 57, 58] were also predicted as tumor suppressor miRNAs (Tables 1 and 3).
Our work gives vital insights into the pathophysiology and progression of BC by presenting essential information on the likely influence of SNPs and distinct regulation patterns on miRNA production and function. Consequently, genetic variations appear to be the appropriate criterion for early diagnosis of BCs in future.
5. Conclusion
In brief, numerous resources were applied following a comprehensive screening of mature miRNAs, including SNPs that play a deciding role on BC to identify SNP’s functional effect in the miRNA gene. We conducted a thorough investigation into the influence of the mined SNPs on miRNA function, including target prediction, miRNA-target interaction, and target expression level. Theoretically, this work revealed a list of possible causative SNPs in the mature miRNA gene that may be addressed for further practical research in BC management.
Acknowledgments
This research was funded by Vietnam National University, Ho Chi Minh City (VNU-HCM), under grant number 562-2020-2118-02.
Credit authorship contribution statement
Thanh Thi Ngoc Nguyen was involved in conceptualization, methodology, formal analysis, data curation, software, writing—original draft, writing—editing, visualization, project administration, and funding acquisition. Thu Huynh Ngoc Nguyen was involved in investigation, validation, and writing—review. Luan Huu Huynh and Hoang Ngo Phan were involved in investigation and writing—review. Hue Thi Nguyen was involved in resources, conceptualization, methodology, supervision, writing—review, project administration, and funding acquisition.
References
- 1.
Liu S, Clouthier SG, Wicha MS. Role of microRNAs in the regulation of breast cancer stem cells. Journal of Mammary Gland Biology and Neoplasia. 2012; 17 (1):15-21 - 2.
Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA: a Cancer Journal for Clinicians. 2021; 71 (3):209-249 - 3.
Wahid F, Shehzad A, Khan T, Kim YY. MicroRNAs: Synthesis, mechanism, function, and recent clinical trials. Biochimica et Biophysica Acta. 2010; 1803 (11):1231-1243 - 4.
Galvão-Lima LJ, Morais AHF, Valentim RAM, Barreto EJSS. miRNAs as biomarkers for early cancer detection and their application in the development of new diagnostic tools. BioMedical Engineering OnLine. 2021; 20 (1):21 - 5.
Loh HY, Norman BP, Lai KS, Rahman N, Alitheen NBM, Osman MA. The regulatory role of MicroRNAs in breast cancer. International Journal of Molecular Sciences. 2019; 20 (19):4940 - 6.
Wang H, Tan Z, Hu H, Liu H, Wu T, Zheng C, et al. microRNA-21 promotes breast cancer proliferation and metastasis by targeting LZTFL1. BMC Cancer. 2019; 19 (1):738 - 7.
Li P, Guo Y, Bledsoe G, Yang Z, Chao L, Chao J. Kallistatin induces breast cancer cell apoptosis and autophagy by modulating Wnt signaling and microRNA synthesis. Experimental Cell Research. 2016; 340 (2):305-314 - 8.
Liu J, Li X, Wang M, Xiao G, Yang G, Wang H, et al. A miR-26a/E2F7 feedback loop contributes to tamoxifen resistance in ER-positive breast cancer. International Journal of Oncology. 2018; 53 (4):1601-1612 - 9.
Tormo E, Adam-Artigues A, Ballester S, Pineda B, Zazo S, Gonzalez-Alonso P, et al. The role of miR-26a and miR-30b in HER2+ breast cancer trastuzumab resistance and regulation of the CCNE2 gene. Scientific Reports. 2017; 7 :41309 - 10.
Dinami R, Ercolani C, Petti E, Piazza S, Ciani Y, Sestito R, et al. miR-155 drives telomere fragility in human breast cancer by targeting TRF1. Cancer Research. 2014; 74 (15):4145-4156 - 11.
Rao X, Di Leva G, Li M, Fang F, Devlin C, Hartman-Frey C, et al. MicroRNA-221/222 confers breast cancer fulvestrant resistance by regulating multiple signaling pathways. Oncogene. 2011; 30 (9):1082-1097 - 12.
Mishra S, Srivastava AK, Suman S, Kumar V, Shukla Y. Circulating miRNAs revealed as surrogate molecular signatures for the early detection of breast cancer. Cancer Letters. 2015; 369 (1):67-75 - 13.
Gong Y, He T, Yang L, Yang G, Chen Y, Zhang X. The role of miR-100 in regulating apoptosis of breast cancer cells. Scientific Reports. 2015; 5 :11650 - 14.
Pakravan K, Babashah S, Sadeghizadeh M, Mowla SJ, Mossahebi-Mohammadi M, Ataei F, et al. MicroRNA-100 shuttled by mesenchymal stem cell-derived exosomes suppresses in vitro angiogenesis through modulating the mTOR/HIF-1alpha/VEGF signaling axis in breast cancer cells. Cellular Oncology (Dordrecht). 2017; 40 (5):457-470 - 15.
Akhavantabasi S, Sapmaz A, Tuna S, Erson-Bensan AE. miR-125b targets ARID3B in breast cancer cells. Cell Structure and Function. 2012; 37 (1):27-38 - 16.
Chen J, Tian W, Cai H, He H, Deng Y. Down-regulation of microRNA-200c is associated with drug resistance in human breast cancer. Medical Oncology. 2012; 29 (4):2527-2534 - 17.
Alhasan L. MiR-126 modulates angiogenesis in breast cancer by targeting VEGF-A-mRNA. Asian Pacific Journal of Cancer Prevention. 2019; 20 (1):193-197 - 18.
Zou C, Xu Q, Mao F, Li D, Bian C, Liu LZ, et al. MiR-145 inhibits tumor angiogenesis and growth by N-RAS and VEGF. Cell Cycle. 2012; 11 (11):2137-2145 - 19.
Shi J, Zhang Y, Zheng W, Michailidou K, Ghoussaini M, Bolla MK, et al. Fine-scale mapping of 8q24 locus identifies multiple independent risk variants for breast cancer. International Journal of Cancer. 2016; 139 (6):1303-1317 - 20.
Bao L, Hazari S, Mehra S, Kaushal D, Moroz K, Dash S. Increased expression of P-glycoprotein and doxorubicin chemoresistance of metastatic breast cancer is regulated by miR-298. The American Journal of Pathology. 2012; 180 (6):2490-2503 - 21.
Meng Y, Zou Q, Liu T, Cai X, Huang Y, Pan J. microRNA-335 inhibits proliferation, cell-cycle progression, colony formation, and invasion via targeting PAX6 in breast cancer cells. Molecular Medicine Reports. 2015; 11 (1):379-385 - 22.
Gao Y, Zeng F, Wu J-Y, Li H-Y, Fan J-J, Mai L, et al. MiR-335 inhibits migration of breast cancer cells through targeting oncoprotein c-met. Tumor Biology. 2015; 36 (4):2875-2883 - 23.
McAnena P, Tanriverdi K, Curran C, Gilligan K, Freedman JE, Brown JAL, et al. Circulating microRNAs miR-331 and miR-195 differentiate local luminal a from metastatic breast cancer. BMC Cancer. 2019; 19 (1):436 - 24.
Yerukala Sathipati S, Ho SY. Identifying a miRNA signature for predicting the stage of breast cancer. Scientific Reports. 2018; 8 (1):16138 - 25.
Farina NH, Ramsey JE, Cuke ME, Ahern TP, Shirley DJ, Stein JL, et al. Development of a predictive miRNA signature for breast cancer risk among high-risk women. Oncotarget. 2017; 8 (68):112170 - 26.
John B, Enright AJ, Aravin A, Tuschl T, Sander C, Marks DS. Human MicroRNA targets. PLoS Biology. 2004; 2 (11):e363 - 27.
Dong H, Lei J, Ding L, Wen Y, Ju H, Zhang X. MicroRNA: Function, detection, and bioanalysis. Chemical Reviews. 2013; 113 (8):6207-6233 - 28.
Ryan BM, Robles AI, Harris CC. Genetic variation in microRNA networks: The implications for cancer research. Nature Reviews. Cancer. 2010; 10 (6):389-402 - 29.
Grimson A, Farh KK, Johnston WK, Garrett-Engele P, Lim LP, Bartel DP. MicroRNA targeting specificity in mammals: Determinants beyond seed pairing. Molecular Cell. 2007; 27 (1):91-105 - 30.
Liu CJ, Fu X, Xia M, Zhang Q, Gu Z, Guo AY. miRNASNP-v3: A comprehensive database for SNPs and disease-related variations in miRNAs and miRNA targets. Nucleic Acids Research. 2021; 49 (D1):D1276-D1D81 - 31.
Kozomara A, Birgaoanu M, Griffiths-Jones S. miRBase: From microRNA sequences to function. Nucleic Acids Research. 2019; 47 (D1):D155-DD62 - 32.
Sherry ST, Ward MH, Kholodov M, Baker J, Phan L, Smigielski EM, et al. dbSNP: The NCBI database of genetic variation. Nucleic Acids Research. 2001; 29 (1):308-311 - 33.
Chen Y, Wang X. miRDB: An online database for prediction of functional microRNA targets. Nucleic Acids Research. 2020; 48 (D1):D127-DD31 - 34.
Wong N, Wang X. miRDB: An online resource for microRNA target prediction and functional annotations. Nucleic Acids Research. 2015; 43 (Database issue):D146-D152 - 35.
Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K. KEGG: New perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Research. 2017; 45 (D1):D353-DD61 - 36.
Park S-J, Yoon B-H, Kim S-K, Kim S-Y. GENT2: An updated gene expression database for normal and tumor tissues. BMC Medical Genomics. 2019; 12 (5):101 - 37.
McGeary SE, Lin KS, Shi CY, Pham TM, Bisaria N, Kelley GM, et al. The biochemical basis of microRNA targeting efficacy. Science (New York, NY). 2019; 366 (6472):eaav1741 - 38.
Rehmsmeier M, Steffen P, Hochsmann M, Giegerich R. Fast and effective prediction of microRNA/target duplexes. RNA. 2004; 10 (10):1507-1517 - 39.
Boissonneault V, Plante I, Rivest S, Provost P. MicroRNA-298 and microRNA-328 regulate expression of mouse beta-amyloid precursor protein-converting enzyme 1. The Journal of biological chemistry. 2009; 284 (4):1971-81 - 40.
Andrés-León E, Cases I, Alonso S, Rojas AM. Novel miRNA-mRNA interactions conserved in essential cancer pathways. Scientific Reports. 2017; 7 (1):46101. - 41.
Brincas HM, Augusto DG, Mathias C, Cavalli IJ, Lima RS, Kuroda F, et al. A genetic variant in microRNA-146a is associated with sporadic breast cancer in a southern Brazilian population. Genetics and Molecular Biology. 2020; 42 (4):e20190278 - 42.
Ahmad M, Shah AA. Predictive role of single nucleotide polymorphism (rs11614913) in the development of breast cancer in Pakistani population. Personalized Medicine. 2020; 17 (3):213-227 - 43.
Tan SC, Lim PY, Fang J, Mokhtar MFM, Hanif EAM, Jamal R. Association between MIR499A rs3746444 polymorphism and breast cancer susceptibility: A meta-analysis. Scientific Reports. 2020; 10 (1):3508 - 44.
Hashemi M, Sanaei S, Rezaei M, Bahari G, Hashemi SM, Mashhadi MA, et al. miR-608 rs4919510 C>G polymorphism decreased the risk of breast cancer in an Iranian subpopulation. Experimental Oncology. 2016; 38 (1):57-59 - 45.
Chen HC, Tseng YK, Chi CC, Chen YH, Yang CM, Huang SJ, et al. Genetic variants in microRNA-146a (C>G) and microRNA-1269b (G>C) are associated with the decreased risk of oral premalignant lesions, oral cancer, and pharyngeal cancer. Archives of Oral Biology. 2016; 72 :21-32 - 46.
Gutierrez-Camino A, Martin-Guerrero I, Dolzan V, Jazbec J, Carbone-Baneres A, Garcia de Andoin N, et al. Involvement of SNPs in miR-3117 and miR-3689d2 in childhood acute lymphoblastic leukemia risk. Oncotarget. 2018; 9 (33):22907-22914 - 47.
Huang Y, Feng Y, Ren H, Zhang M, Li H, Qiao Y, et al. Associations of genetic variations in MicroRNA seed regions with acute adverse events and survival in patients with rectal cancer receiving postoperative Chemoradiation therapy. International Journal of Radiation Oncology, Biology, Physics. 2018; 100 (4):1026-1033 - 48.
Raad M, Salehi Z, Habibzaadeh Baalsini M, Mashayekhi F, Saeidi SH. Association of rs2620381 polymorphism in miR-627 and gastric cancer. British Journal of Biomedical Science. 2020; 77 (2):76-80 - 49.
Mao Y, Zou C, Meng F, Kong J, Wang W, Hua D. The SNPs in pre-miRNA are related to the response of capecitabine-based therapy in advanced colon cancer patients. Oncotarget. 2018; 9 (6):6793-6799 - 50.
Zhang P, Wang J, Lu T, Wang X, Zheng Y, Guo S, et al. miR-449b rs10061133 and miR-4293 rs12220909 polymorphisms are associated with decreased esophageal squamous cell carcinoma in a Chinese population. Tumour Biology. 2015; 36 (11):8789-8795 - 51.
Wang R, Zhang J, Jiang W, Ma Y, Li W, Jin B, et al. Association between a variant in microRNA-646 and the susceptibility to hepatocellular carcinoma in a large-scale population. ScientificWorldJournal. 2014; 2014 :312704 - 52.
Li Y, Zhu H, Wang J, Qian X, Li N. miR-4513 promotes breast cancer progression through targeting TRIM3. American Journal of Translational Research. 2019; 11 (4):2431-2438 - 53.
Xiong M, Wang P, Pan B, Nie J, Wang S, He B. The diagnostic and prognostic values of microRNA-196a in cancer. Bioscience Reports. 2021; 41 (1):BSR20203559 - 54.
Gao Y, Ma H, Gao C, Lv Y, Chen X, Xu R, et al. Tumor-promoting properties of miR-8084 in breast cancer through enhancing proliferation, suppressing apoptosis and inducing epithelial-mesenchymal transition. Journal of Translational Medicine. 2018; 16 (1):38 - 55.
Hadavi R, Mohammadi-Yeganeh S, Razaviyan J, Koochaki A, Kokhaei P, Bandegi A. Expression of Bioinformatically candidate miRNAs including, miR-576-5p, miR-501-3p and miR-3143, targeting PI3K pathway in triple-negative breast cancer. Galen Medical Journal. 2019; 8 :e1646 - 56.
Zhan Y, Du J, Min Z, Ma L, Zhang W, Zhu W, et al. Carcinoma-associated fibroblasts derived exosomes modulate breast cancer cell stemness through exonic circHIF1A by miR-580-5p in hypoxic stress. Cell Death Discovery. 2021; 7 (1):141 - 57.
Wang GP, Mou ZL, Xu YY, Liu GX, Wang DM, Zhang HP. LINC01096 knockdown inhibits progression of triple-negative breast cancer by increasing miR-3130-3p. European Review for Medical and Pharmacological Sciences. 2019; 23 (17):7445-7456 - 58.
Shi SH, Jiang J, Zhang W, Sun L, Li XJ, Li C, et al. A novel lncRNA HOXC-AS3 acts as a miR-3922-5p sponge to promote breast cancer metastasis. Cancer Investigation. 2020; 38 (1):1-12