Skip to main content

MACE-Seq-based coding RNA and TrueQuant-based small RNA profile in breast cancer: tumor-suppressive miRNA-1275 identified as a novel marker

Abstract

Introduction

Disruption of cellular processes in the breast by abnormally expressed miRNA is characterized to develop cancer. We aimed to identify the differential expression of small RNAs (sRNAs) and mRNAs in formalin-fixed paraffin-embedded (FFPE) tissue of the breast cancer (BC) and normal adjacent tissue (NAT). Another aim is to determine the differential expression of miR-1275 as a novel biomarker for BC and also identify its target genes.

Methods

TrueQuant method for analysis of sRNA expression and MACE-sequencing method for analysis of gene expression were used analyzing. The RT-qPCR technique was used to confirm miR-1275 down expression. Target genes of miR-1275 were computationally identified using target prediction sites and also the expression level of them was experimentally determined among the expressed genes.

Results

TrueQuant findings showed that 1400 sRNAs were differentially expressed in the FFPE tissue of two Kurdish cases with BC, as compared to NAT. Among the sRNAs, 29 small RNAs were shown to be significantly downregulated in BC cells. The RT-qPCR results confirmed that miR-1275 was significantly down-expressed in 20 Kurdish cases with BC compared to NAT. However, Overall survival (OS) analysis revealed that the correlation between the expression level of miR-1275 and clinical significance was highly corrected in cases with BC (OS rate: P = 0.0401). The MACE-seq results revealed that 26,843 genes were differentially expressed in the BC tissue compared to NAT, but 7041 genes were displayed in a scatter plot. Furthermore, putative target genes (DVL3, PPP2R2D, THSD4, CREB1, SYT7, and PRKACA) were computationally identified as direct targets of miR-1275 in several target predicted sites. The MACE-seq results revealed that the expression level of these targets was increased in BC tissue compared to NAT. The level of these targets was negatively associated with miR-1275 expression. Finally, the role of down-regulated miR-1275 on its targets in biological mechanisms of BC cells was identified; including cell growth, proliferation, movement, invasion, metastasis, and apoptosis.

Conclusion

Down-expressed miR-1275, a tumor suppressor, is a novel biomarker for early detection of BC. DVL3, PPP2R2D, THSD4, CREB1, SYT7, and PRKACA are newly identified to be targeted by miR-1275.

Peer Review reports

Introduction

Malignant breast tumor is a prominent type of cancers mostly diagnosed in females, and the second most frequent malignancy-associated deaths worldwide, especially in the US and Asian countries [1, 2]. Approximately two million females are annually diagnosed and more than 620,000 deaths are newly recorded every year [3, 4]. Frequently, breast cancer (BC) is developed as a result of a genomic mutation. However, about 10% of BCs is inheritably come down from parents to their generations; whereas, more than 85% of BCs is developed in their lifetime [4, 5]. Inherited abnormalities in TP53 and PTEN genes were studied to result in the high risk of malignant breast tumor progression [6, 7].

Gene expression profiling has recently played a critical role in medicinal selection for BC subtypes. The analysis of BC gene expression can be used for the molecular category of BC subtypes [8,9,10]. Two studies reported that this classification facilitates the determination of cure doses. The molecular subtypes of BC can be categorized into luminal-A, luminal-B (including HER2+ /−), HER2+, and triple-negative (TN) [11, 12]. These subtypes are pivotal for cure choice and are correlated to the biological characteristics of BC.

MicroRNA (miRNA or miR), which is a type of small non-coding RNAs (sncRNAs), is synthesized from eukaryotic genomes, consisting of a single-stranded RNA of about 19-22 nt in length [13, 14]. These short RNAs are described as regulators of coding (cRNAs) and sncRNAs in eukaryotic cells because they are involved in silencing and stabilizing targeted cRNAs [14, 15]. MiRNAs play also several regulatory roles in several cellular processes; cell development, proliferation, migration, invasion, and death [16]. Because more than 50% of RNA molecules has been detected to be controlled by miRNAs, these mRNAs were damaged because of the effect of aberrant miRNA in malignant cells.

MiRNAs in human malignant tumors were found to act as oncogenes or tumor suppressors. Oncogenic miRNAs in tumor progress play a negative role in stimulating genes that regulate cell development and the apoptosis process. Tumor suppressive miRNAs have a key role in silencing genes that modulate cell development and apoptosis [17, 18]. When normal cells do not undergo normal growth and apoptosis process, they normally cause tumor creation. Numerous recent experiments show that numerous miRNAs are directly implicated in modulating cell growth, proliferation as well as apoptosis [19, 20]. They play major roles in the pathogenesis of several human malignancies; such as breast, colorectal cancer, lung, leukemia liver, and brain [16, 18]. The miRNA expression level may be either down- or up-regulated in these cancers. Several molecular techniques, such as RNA sequencing, miRNA microarray, RT-PCR, and northern blot are applied to determine the expression level of them.

Numerous miRNAs have been recognized to be implicated in the pathogenesis of human breast cancer. It was found that the expression patterns of miR-145, −125b, − 155, and -21were significantly downregulated. In breast malignant cells, these miRNAs were observed to be associated with pathologic properties; cell proliferation, expression of progesterone, and estrogen receptors [21]. A recent study revealed that tumor-suppressive miRNA-204-5p plays a key role in targeting several oncogenic genes that are closely connected to BC pathogenesis [20]. Complete information on miR-1275 expression level and its targets in BC have not been available; whereas, the expression profile of which has been analyzed in some human cancers. A study reported that down-expressed miR-1275 leads to overexpression claudin11 in cancer stem cells or tumor-initiating cells (CSCs/TICs) [22]. According to a study carried out on young women with BC, 6 miRNAs; including miR-1275, miR-1228*, miR-139, miR-92b, miR-1207, and miR-3196, were involved in the processes of cell movement, proliferation, and invasion [23]. It was also found that this miRNA was found to be downregulated in gastric cancer [24]. Another study found that miR-1275 expression level was significantly abnormally deregulated in Ewing’s Sarcoma (ES) [25]. The significance of a large number of miRNAs has been reported to become an appropriate biomarker for human cancer diagnostics. However, the significance of miR-1275 in BC is not reported. The objective of this study was to determine the expression level of miR-1275 as a biomarker for BC diagnostics. Another objective is to identify the potential targeted genes of this miRNA.

Material and method

Collection of FFPE-blocks of BC samples

All procedures done in this research study involving human participants, human material, or human data were followed and approved by local Human Research Ethics Committee (HREC) at Science College in Salahaddin University-Erbil (Reference no. 4c/132). All methods while performing the study were also performed in accordance with the 1964 Helsinki Declaration and the written informed-consent and Permission for the publication were obtained from all research participants. FFPE-Block of 22 Kurdish cases with BC were collected at clinicopathological laboratories, called Al-Mufti and Luay. Any chemotherapy or radiotherapy was taken by none of these cases. For each patient, two paraffin blocks (NAT and BC tissue) were collected. We diagnosed the malignant and non-malignant areas. The non-malignant tissues were taken nearly 2 cm away from the malignant area. Clinical features of 22 cases were obtained using a questionnaire. The features were displayed in Table 1.

Table 1 Clinicopathological features of 22 cases with BC

Total RNA extraction and library preparation for MACE-seq and small RNA seq

MACE-sequencing for mRNA expression and small RNA sequencing for small RNA expression were performed for two cases 1 and 2 (Table 1). The practical part of this task was performed in GenXPro GmbH laboratory. The protocol of total RNA extraction, small RNA separation, and library preparation was prepared from GenXPro GmbH in Germany. For each sample, total RNA was prepared from 10 μm of FFPE tissues employing the MACE-seq kit (GenXPro GmbH, Frankfurt, Germany). This kit contains all reagents to separate the molecule of the messenger RNAs (mRNAs). Then TrueQuant small RNA kit was used to separate small RNAs (sRNAs); including miRANs, small nucleolar RNAs (snoRNAs), small interfering RNAs (siRNAs), small nuclear RNAs (snRNAs), PIWI-interacting RNAs (piRNAs). Massive Analysis of cDNA Ends (MACE)-seq technique was used for the analysis of mRNA expression (GenXPro GmbH, Frankfurt, Germany). In MACE-seq, a specific molecular barcode, known as TrueQuant barcode, was used to barcode each mRNA molecule before PCR steps. The large molecules of the total RNA were utilized for the construction of the MACE-libraries and the small molecules for the construction of sncRNA libraries.

Four MACE-libraries (two mRNA libraries and two sRNA libraries) were generated for 2 BC and 2 NAT samples. Poly-adenylated mRNA was purified from the large fraction of total RNA using MACE-Kit (GenXPro GmbH, Frankfurt, Germany). First and second-strand cDNA was synthesized using the SuperScript® III First-Strand Synthesis System (GenXPro GmbH, Frankfurt, Germany). Then the cDNA products were attached to barcoded poly-T adapters. These barcodes are appropriate to attach to the MACE-seq flow well. The cDNA cap structures were specifically biotinylated. The full-length cDNAs which had biotinylated cap are attracted by streptavidin beads and that were randomly sheared to be about 300 bps, GenXPro’s TrueQuant barcodes were ligated. Four barcoded samples were sequenced simultaneously in one lane of an Illumina Hiseq2000 with 1 × 100 bps. To prepare miRNA libraries, small RNAs (16–25 bps) were chosen from the entire RNA molecules utilizing the flashPAGE™ Fractionator System (Life Technologies). TrueQuant Adapters were directly attached to the small RNAs, basically as mentioned by Hafner and his colleagues [26]. A Hiseq2000 machine (Illumina) was used to sequence the miRNA group comprising p5 and p7 adapters. The GenXPro GmbH prepared MACE and sRNA libraries. For removing flanking adapters, the organized reads were cut out for high-quality sequences. Bowtie 2 tool was later used for aligning the sorted reads to the nominated reference sequences and annotating with corresponding properties [27]. HT-seq tool was used to execute the quantification of mapped reads to each gene [28]. DESeq2 was used to execute differential expression analysis [29], and also is based on negative binomial generalized linear models. Results were compiled into a final table including significance parameters; including P.value, FDR and log2Fold-Changes. R-scripts were used to perfume Final Data visualization of the significantly expressed gene and the down/over-expressed genes.

Total RNA extraction and cDNA synthesis for miR-1275

From the TrueQuant results, miR-1275 was selected to confirm its expression from 20 Kurdish cases using Real Time-quantitative polymerase chain reaction (RT-qPCR) technique. The practical part of this task was performed in Salahuddin University Research Center (SURC). Forty block specimens (20 NAT and 20 BC tissues) were used for analysis of differential expression of miR-1275. For each experiment, 10 μm of the FFPE section was put in the microcentrifuge tube. Deparaffinization of the tissue was performed by incubation in 1 ml of xylene for 5 min at 50 °C. After incubation, the tissue was centrifugated to produce pellets at 14,000 rpm for 2 min. Then, the deparaffinized pellet was rehydrated with 1 ml of absolute ethanol. The tissue was centrifugated at 14,000 rpm for 2 min. The rehydration and centrifugation were repeated again. The tissue pellet was air-dried for 20 at room temperature. Then, the tissue pellet was homogenized using a PowerGen 125 Tissue Homogenizer. The total RNA molecules including miRNA were extracted using FFPE RNA/DNA Purification Plus kit (Cat. No. 54300, NORGEN BIOTEK CORP, Canada). Complementary DNA (cDNA) was synthesized using the miRNA All-In-One cDNA Synthesis Kit (Cat. No. G898, abmgood company, US).

Determination of differential expression of miR-1275 using RT-qPCR

All the items of RT-qPCR technique were ordered from abmgood company, US. RT-qPCR applications were performed to amplify the targeted miRNA among total cDNA molecules. For each solution well, the total volume was 20 μl which included 10 μl of BrightGreen miRNA qPCR MasterMix-ROX (Cat. No. MasterMix-mR), 0.5 μl for each forward and reverse primers (Cat. No. MPH01104), 1.4 μl of template cDNA, and 7.6 μl of nuclease-free water. In this study, two Universal/housekeeping miRNA primers were used, which represented the U6–2 primers (Cat. No. MPH0001) and SNORD44 primers (Cat. No. MPH0003). qPCR reaction was performed setting up the following 3-step cycling program. Enzyme activation was at 95 °C for 10 min, 40 cycles were set up for denaturing at 95 °C for 10 s, annealing at 63 °C for 20 s and extension at 72 °C for 20 s.

Most common putative targeted genes regulated by miR-1275

In this study, 11 databases were searched for finding the most common putative targets of miR-1275 (Table 2). Six putative targets (DVL3, PPP2R2D, THSD4, CREB1, SYT7, and PRKACA) were determined to possess binding sequence to miR-1275. Then differential expression of them was identified in the MACE-seq findings. GraphPad Prism, Version 8.0.1, was used to reveal the differential expression of these selected target genes.

Table 2 Brief information on target predicted databases was shown to find putative targets possessing binding sequence to miR-1275

Gene enrichment analysis of gene ontology for target genes of hsa-miR-1275

To study the biological role of hsa-miR-1275 in the BC, the online tool GenXpro (http://tools.genxpro.net/modules/GO Enrichment Tool v2) was applied to execute for the predicted target genes of miR-1275 to determine its function. The gene enrichment analyses of Gene Ontology (GO) include Cellular Component (CC), Biology Process (BP), and Molecular Function (MF).

Analysis of clinicopathological data associated with BC

Association between miR-1275 and its target genes was computationally analyzed to determine the clinical significance using databases of cBioPortal (http://www.cbioportal.org/) and OncoLnc (http://www.oncolnc.org/). Clinical data and expression levels of the miR-1275 and its target genes were gained from these sites and then were downloaded on 10 September 2020. Then overall survival (OS) curve per day was designed using R software. The data and R scripts will be stored for future studies.

Result

Construction of expression profile of sRNAs for BC by TrueQuant technique

By comparing sRNA expression profiles of BC and NAT, 1400 sRNAs were filtered out by SAM software. The raw data were then standardized and log2-transformed to show on a scatter plot (Fig. 1a). Among 1400 sRNAs, 723 non-coding RNAs were downregulated, which are located on the right and lower side of the trumpet plot but 677 sRNAs were upregulated, which are on the left and upper side of the trumpet. Each dot on the scatter plot represents the sRNA. Among 1400 sRNAs, 520 miRNAs were differentially expressed. 185 miRNAs were down-expressed, but 335 miRNAs were overexpressed. The x-axis denotes the data of the NAT and the y-axis denotes the data of the BC. A correlation plot was then constructed to show expression levels of sRNAs between BC and NAT (Fig. 1b). Blue color denoted the correlation of sRNAs between the BC and NAT but red color denoted the misrelation of sRNAs. However, a Heat map was designed to show 29 miRNAs which were markedly downregulated in BC compared to NAT (Fig. 1c). Moreover, the Limma package in R language was utilized to searching differentially expressed ncRNAs (DEncRNAs) between BC and NAT. Bonferroni in the multtest package was applied to adjust the P. value into the FDR. The FDR < 0.05 and |log2 FC| > 0.5 were applied as the cutoff criteria for the DEncRNAs. Table 3 shows the information on the accession number, P. value, genome loci, Log2 fold change (Log2FC), and False Discovering Rate (FDR) of these 29 miRNAs. In this study, has-miR-1275, which was underlined with red color in the heat map, was focused to identify expression level and sequence.

Fig. 1
figure1

Differential expression analysis of non-coding RNAs by MACE-sequencing. a the expression profile of 1400 sRNAs in BC compared to NAT is plotted. Red dots represent the sample sRNAs. A light green dot represents miR-1275 expression level in BC tissue compared to NAT. b Heat-map based clustering of several downregulated miRNAs in BC compared to NAT. MiR-1275 is underlined. c correlation between BC and NAT in differential expression of sRNAs

Table 3 Comparison of marked down-expressed miRNAs in BC with NAT

Confirmation of miR-1275 expression level by RT-qPCR

In small RNA sequencing findings, the miR-1275 expression level in 2 cases with BC was observed to be significantly downregulated compared to NAT. The p.value of which was 0.614 (Fig. 2a). Then, this result was confirmed using an RT-qPCR machine in 20 other cases with BC. The miR-1275 expression level was detected to be markedly decreased in BC compared to NAT and the P-value of this was 0.001 (Fig. 2b). The mature sequence of miR-1275 in the BC and NAT was sequenced and made up of 17 nucleotides (10G, 3 T, 2C, and 2A). The mature sequence in both is also identical (Fig. 2c). However, Kaplan–Meier overall survival curve was designed to show the effect of the miR-1275 expression on the prognosis of cases with BC. Data of overall survival (OS) was gained from The Cancer Genome Atlas (TCGA) database and analyzed by the R program. OS curve displayed that cases were separated into two classes according to miR-1275 expression. The decreased miR-1275 (P-value = 0.0401) was related to overall survival in cases with BC (Fig. 2d).

Fig. 2
figure2

a Comparison of differential expression of miR-1275 in BC compared to NAT using MACE-seq method. b Comparison of differential expression of miR-1275 in BC compared to NAT using RT-qPCR method. c mature sequence of miR-1275 in NAT and BC was shown to be identical. d Kaplan–Meier overall survival curve designed to show the differential expression of miR-1275 related to overall survival in the patients with the BC

Comparison of differential expression of genes in BC with NAT tissues

In MACE-sequencing results, differentially expressed genes in two cases with BC were compared with NAT tissues (Fig. 3a). Twenty-six thousand eight hundred forty-three differentially expressed genes (P ≤ 0.05) were filtered out by SAM software. In order to show genes that were more significantly different in their expressions, 7041 genes were standardized and log2-transformed to show on a scatter plot. Three thousand six hundred twenty-four genes of which were significantly overexpressed and 3417 genes were significantly down-expressed in BC compared to NAT. The P. value for that was ranged from smaller (Blue) to greater (Red). Each point on the scatter plot represents the gene. The x-axis denotes the data of the NAT and the y-axis denotes the data of the cancerous tissue.

Fig. 3
figure3

a Scatter plot analysis of gene expression profile displays up- or down-regulation of genes in BC tissue, compared to NAT. Each point denotes the average value of one transcript in the experiment. The expression difference is taken account of significance for a P. value (0.05). b Outlined points and names denote the selected target genes for miR-1275

Candidate target genes regulated in BC by miR-1275

More than ten computational prediction programs were utilized for discovering the strongest candidate genes possessed miR-1275 binding sites in the 3′ -UTR. Six predicted genes (DVL3, PPP2R2D, THSD4, CREB1, SYT7, and PRKACA) were selected to have the binding sites to miR-1275 (red bonds and bps between targets and miR-1275 in Fig. 4). The information on these six predicted genes was summarized in Table 4. Eleven databases showed that DVL3 and PPP2R2D possessed the binding site to miR-1275; whereas, THSD4, and CREB1 were confirmed in ten prediction programs to be targeted by miR-1275, but SYT7 and PRKACA were confirmed by six tools to be predicted targets. These putative target genes are important for biological analysis of the BC tissues because the over-or down-expression of which can play a damaging role in several cellular processes and contribute to cancer progression and tumorigenesis.

Fig. 4
figure4

Target genes and miR-1275 are combined in the seed region, including 6 to 8 nucleotides in the 5′ end, showing in red color. Differential expression level of target genes in two cases with BC compared to NAT

Table 4 Candidate target genes possessing binding sequence to miR-1275 was shown

Determination of expression level of candidate target genes in MACE-seq findings

In MACE-seq finding, the six predicted genes (DVL3, PPP2R2D, THSD4, CREB1, SYT7, and PRKACA) were pointed and named in the BC cells as compared to NAT (Fig. 3b). Then the differential expression of which and their binding sites to has-miR-1275 were shown in Fig. 4. Among the targets, the SYT7 gene was more overexpressed in BC, as compared to the PRKACA gene. The overexpression level of the THSD4 gene was also higher than the up-regulation of PPP2R2D and DVL3 genes. CREB1 was upregulated but more over-expressed than the ST73 gene. Table 5 demonstrates the information on P.value, False Discovering Rate (FDR), and Fold Change (FC) of these six predicted genes possessed miR-1275 binding sites in the 3′ -UTR. They were identified as potentially modulated by miR-1275 using computational prediction databases and TCGA algorithm.

Table 5 Experimentally validated target genes of miR-1275 in BC

However, the relationship between the expression level of these 6 candidates and histopathological significance was analyzed based on data from the TCGA database. Among 204 target genes, DVL3: P = 6.98E-04, PPP2R2D: P = 1.53E-03, THSD4: P = 2.28E-12, CREB1: P = 3.12E-01, SYT7: P = 1.28E-48, and PRKACA: P = 5.68E-02 were markedly relationship with worse prognosis in cases with BC. For all targets, the Kaplan Meier curve was designed to show the correlation between OS and gene expression in cases (Fig. 5).

Fig. 5
figure5

An association between the expression levels of six genes (DVL3, PPP2R2D, TSHD4, CREB1, SYT7, and PRKACA) and histopathological significance was shown using data from TCGA database. The Kaplan–Meier overall survival curves show that patients with BC were separated into 2 classes according to their expression levels

Enrichment analysis of gene ontology (GO)

Based on the GenXpro databases, GO enrichment analysis of differentially expressed genes (DEGs) was carried out for determination of the biological role of the target genes of hsa-miR-1275 in the BC. The up-regulated DEG of DVL3, PPP2R2D, THSD4, CREB1, SYT7, and PRKACA was enriched in the GO_terms, such as the extracellular exosome (GO:0070062; P = 0.014527), cytosol (GO:0005829; P = 0.041234), phosphoprotein (GO:0004721; P = 0.036901) (Fig. 6).

Fig. 6
figure6

The gene ontology functional enrichment analyses of differentially expressed genes (DEGs) for targeted genes of hsa-miR-1275

The role of miR-1275 by targeting selected putative genes in BC

MiRNAs are implicated in silencing mRNA transcripts through matching or mismatching with target mRNAs. As hypothesis of microRNA biogenesis, significant strands of miR-1275 come from miRNA duplex are joined into the RISC protein and modulate the mRNA transcripts, but minor strands are broken down and cannot modulate gene expression. Figure 7 shows the relationship between miR-1275 and these target genes. miR-1275 can play an essential role in regulating several biological mechanisms; including cell growth, migration, differentiation, proliferation, and apoptosis. In this study, miR-1275 was observed to regulate six genes related to tumor development. In MACE-seq. Findings, these six genes (DVL3, PPP2R2D, THSD4, CREB1, SYT7, and PRKACA) were detected to be over-expressed in BC cells as compared with NAT. Downregulation of miR-1275 in breast cancer promotes cancer cell proliferation, cell differentiation, tumor growth, invasion and migration and also inhibits apoptosis through these six genes. All targeted genes were negatively regulated by miR-1275. PPP2R2D acts as a tumor suppressor in the signaling pathway in BC. The overexpression of which decreases AKT and RACK1 abilities. Then these regulators increase cell survival and migration. DVL3 is implicated in the breast cancer pathways and negatively controlled by miR-1275. The up-regulation of which increases cancer cell proliferation, migration, and invasion. The cancer cell proliferation ability is also increased when miR-1275 becomes overexpressed. Moreover, CREB1 was found to reduce the apoptosis process and increase cell proliferation in breast cancer; whereas, PRKACA plays a key role in tumorigenesis and development of BC. However, the function of THSD4 and SYT7, currently unidentified, may boost tumor growth in breast cancer.

Fig. 7
figure7

miR-1275 putative targets and their roles in the BC

Discussion

There is evidence that a single microRNA modulates multiple protein-coding and non-coding genes in different ordinary cells and cancerous cells. In human cancer cells, new RNA groups can be detected using the specific nature of microRNA from relevant microRNA analysis. Multiple high-throughput approaches, such as DNA microarrays, MACE-sequencing, PCR-based arrays, and RNA-sequencing, are now available and have made microRNA expression profiles of BC, showing the irregular expression of numerous miRNAs [3, 30,31,32]. One approach to detect the most essential miRNA from numerous miRNAs is to detect differential expression of miRNAs which have been shown in numerous experiments. Several types of research have reported that multiple miRNAs, such as miR-100, miR-107, miR-205-3p, miR-122 and miR-99a-5p, are continuously down-expressed and act as tumor-suppressive miRNA in BC cells [4, 30,31,32,33,34,35]. In this study, these miRNAs were found to be downregulated in tumor cells but only miR-1275 was focused and their putative target genes were newly explored in paraffin-embedded BC tissues. Several types research have reported that this miRNA is down-expressed in gastric and nasopharyngeal carcinoma and function as a tumor suppressor [36,37,38]; whereas, this miRNA is overexpressed in non-small lung cell cancer, squamous carcinoma, and chronic myelogenous leukemia [39,40,41]. Some recent studies revealed that this downregulated miRNA was detected to have an essential effect on cancer cell proliferation, migration, invasion, metastasis, and angiogenesis through targeting multiple oncogenic genes HOXB5, WNT7B and LncRNA HAND2-AS1 [36, 38, 41]. One previous study showed that miR-1275 regulates IGF1, NFIX, Claudin11 in very young women with BC [23]. Whereas, down-expression of miR-1275 in all subtypes of paraffin BC tissues was not fully investigated. In this study, down-expression of miRNA was observed in all subtypes of paraffin-embedded BC tissues of 21 cases of different ages.

After that, the miR-1275-modulated putative targets and their pathways were aimed to explain in the cells of BC. The expression of DVL3, PPP2R2D, TSHD4, CREB1, SYT7, and PRKACA were experimentally observed to be high in the cells of BC. Based on the databases of miRNA target prediction, they were selected and closely correlated with poor prognosis. The biological role of these genes was determined according to the GO terms enrichment with DEGs. These genes were involved in the GO terms of the cytosol, extracellular exosome, protein binding, phosphoprotein and biological adhesion. Two previous studies reported that the overexpression of which could induce a significant impact on the biological processes of the BC progression [42, 43]. Among these candidate genes, three genes (PPP2R2D, DVL3, and CREB1) were shown to be strongly targeted by the miR-1275 in the BC cells. Studies showed that these regulators were found to reduce cell survival and migration in cancer cells [42, 43]. DVL3 is observed to be implicated in the breast cancer pathways [36] and negatively regulated by miR-1275. The up-regulation of this gene can increase cancer cell proliferation, migration, and invasion in BC cells. The overexpression of which increases cancer cell proliferation ability [44]. Another target gene, CREB1, and PRKACA show also a negative correlation with miR-1275 level. Whereas CREB1 was found to reduce the apoptosis process and increase cell proliferation in breast cancer [45], PRKACA plays a key role in tumorigenesis and development of BC [46]. However, the role of THSD4 and SYT7, currently unidentified, may enhance tumor growth in a variety of cancers, especially breast cancer [47, 48].

Conclusion

In the present study, differential expression profiles of mRNA transcripts and sRNAs were analyzed in two cases with BC compared to NAT. Down expression of miR-1275 was measured using the RT-qPCR technique. Down-expression of which develops breast cancer by increasing the activity of biological processes; such as growth, migration, invasion, and metastasis. Upregulated miR-1275 prevented BC development by modulating the direct expression of DVL3, PPP2R2D, TSHD4, CREB1, SYT7, and PRKACA. This is the first study revealing that miR-1275 function as a tumor-suppressive miRNA in BC cells, regulating numerous targets which were closely related with BC pathogenesis and oncogenesis.

Availability of data and materials

Although row data of MACE-sequencing and sRNA may be available in the database of GenXpro, at (https://genxpro.net/, available project: 20190618_Majed_MACE SMajed 2019-06-18 18:44:18.841903 + 00:00), These data will be further studied for another research in the future. The findings described in this manuscript were provided by the Co-author.

Abbreviations

MiRNA:

MicroRNA

FFPE:

Formalin fixed paraffin embedded

BC:

Breast cancer

NAT:

Normal adjacent tissue

ncRNA:

Non-coding RNA

sRNA:

Small RNA

HREC:

Human Research Ethics Committee

MACE:

Massive Analysis of cDNA Ends

E.R.:

Estrogen receptor

Pg. R.:

Progesterone receptor

HER2:

Human epidermal growth factor receptor 2

OS:

Overall survival

TCGA:

The Cancer Genome Atlas

FDR:

False Discovering Rate

FC:

Fold change

DEncRNA:

Differentially expressed ncRNA

DGE:

Differentially expressed gene

GO:

Gene ontology

CC:

Cellular compartment

MF:

Molecular function

BP:

Biological process

References

  1. 1.

    Control, C.f.D. and Prevention. Cancer survivors--United States, 2007. MMWR Morb Mortal Wkly Rep. 2011;60(9):269.

    Google Scholar 

  2. 2.

    Jiang X, Tang H, Chen T. Epidemiology of gynecologic cancers in China. J Gynecol Oncol. 2018;29(1):e7.

    Article  Google Scholar 

  3. 3.

    Adhami M, Haghdoost AA, Sadeghi B, Malekpour Afshar R. Candidate miRNAs in human breast cancer biomarkers: a systematic review. Breast Cancer. 2018;25(2):198–205. https://doi.org/10.1007/s12282-017-0814-8.

    Article  PubMed  Google Scholar 

  4. 4.

    Toda H, Seki N, Kurozumi S, Shinden Y, Yamada Y, Nohata N, et al. RNA-sequence-based microRNA expression signature in breast cancer: tumor-suppressive miR-101-5p regulates molecular pathogenesis. Mol Oncol. 2020;14(2):426–46. https://doi.org/10.1002/1878-0261.12602.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    DeSantis C, Howlader N, Cronin KA, Jemal A. Breast cancer incidence rates in US women are no longer declining. Cancer Epidemiol Prev Biomark. 2011;20(5):733–9. https://doi.org/10.1158/1055-9965.EPI-11-0061.

    Article  Google Scholar 

  6. 6.

    Campeau PM, Foulkes WD, Tischkowitz MD. Hereditary breast cancer: new genetic developments, new therapeutic avenues. Hum Genet. 2008;124(1):31–42. https://doi.org/10.1007/s00439-008-0529-1.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Economopoulou P, Dimitriadis G, Psyrri A. Beyond BRCA: new hereditary breast cancer susceptibility genes. Cancer Treat Rev. 2015;41(1):1–8. https://doi.org/10.1016/j.ctrv.2014.10.008.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Perou CM, et al. Molecular portraits of human breast tumours. Nature. 2000;406(6797):747–52.

    CAS  Article  Google Scholar 

  9. 9.

    Rouzier R, Perou CM, Symmans WF, Ibrahim N, Cristofanilli M, Anderson K, et al. Breast cancer molecular subtypes respond differently to preoperative chemotherapy. Clin Cancer Res. 2005;11(16):5678–85. https://doi.org/10.1158/1078-0432.CCR-04-2421.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Eroles P, Bosch A, Alejandro Pérez-Fidalgo J, Lluch A. Molecular biology in breast cancer: intrinsic subtypes and signaling pathways. Cancer Treat Rev. 2012;38(6):698–707. https://doi.org/10.1016/j.ctrv.2011.11.005.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Park S, Koo JS, Kim MS, Park HS, Lee JS, Lee JS, et al. Characteristics and outcomes according to molecular subtypes of breast cancer as classified by a panel of four biomarkers using immunohistochemistry. Breast. 2012;21(1):50–7. https://doi.org/10.1016/j.breast.2011.07.008.

    Article  PubMed  Google Scholar 

  12. 12.

    Prat A, Pineda E, Adamo B, Galván P, Fernández A, Gaba L, et al. Clinical implications of the intrinsic molecular subtypes of breast cancer. Breast. 2015;24:S26–35. https://doi.org/10.1016/j.breast.2015.07.008.

    Article  PubMed  Google Scholar 

  13. 13.

    Duttagupta R, Jiang R, Gollub J, Getts RC, Jones KW. Impact of cellular miRNAs on circulating miRNA biomarker signatures. PLoS One. 2011;6(6):e20769. https://doi.org/10.1371/journal.pone.0020769.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  14. 14.

    Tomimaru Y, Eguchi H, Nagano H, Wada H, Kobayashi S, Marubashi S, et al. Circulating microRNA-21 as a novel biomarker for hepatocellular carcinoma. J Hepatol. 2012;56(1):167–75. https://doi.org/10.1016/j.jhep.2011.04.026.

    CAS  Article  PubMed  Google Scholar 

  15. 15.

    Gebert LF, MacRae IJ. Regulation of microRNA function in animals. Nat Rev Mol Cell Biol. 2019;20(1):21–37. https://doi.org/10.1038/s41580-018-0045-7.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Treiber T, Treiber N, Meister G. Regulation of microRNA biogenesis and its crosstalk with other cellular pathways. Nat Rev Mol Cell Biol. 2019;20(1):5–20. https://doi.org/10.1038/s41580-018-0059-1.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Zhang B, Pan X, Cobb GP, Anderson TA. microRNAs as oncogenes and tumor suppressors. Dev Biol. 2007;302(1):1–12. https://doi.org/10.1016/j.ydbio.2006.08.028.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Rufino-Palomares EE, et al. MicroRNAs as oncogenes and tumor suppressors. MicroRNAs Med. 2014:223–43.

  19. 19.

    Cheng AM, Byrom MW, Shelton J, Ford LP. Antisense inhibition of human miRNAs and indications for an involvement of miRNA in cell growth and apoptosis. Nucleic Acids Res. 2005;33(4):1290–7. https://doi.org/10.1093/nar/gki200.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Toda H, Kurozumi S, Kijima Y, Idichi T, Shinden Y, Yamada Y, et al. Molecular pathogenesis of triple-negative breast cancer based on microRNA expression signatures: antitumor miR-204-5p targets AP1S3. J Hum Genet. 2018;63(12):1197–210. https://doi.org/10.1038/s10038-018-0510-3.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Iorio MV, Ferracin M, Liu CG, Veronese A, Spizzo R, Sabbioni S, et al. MicroRNA gene expression deregulation in human breast cancer. Cancer Res. 2005;65(16):7065–70. https://doi.org/10.1158/0008-5472.CAN-05-1783.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    Kwon MJ. Emerging roles of claudins in human cancer. Int J Mol Sci. 2013;14(9):18148–80. https://doi.org/10.3390/ijms140918148.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Peña-Chilet M, Martínez MT, Pérez-Fidalgo JA, Peiró-Chova L, Oltra SS, Tormo E, et al. MicroRNA profile in very young women with breast cancer. BMC Cancer. 2014;14(1):529. https://doi.org/10.1186/1471-2407-14-529.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  24. 24.

    Liu D, Hu X, Zhou H, Shi G, Wu J. Identification of aberrantly expressed miRNAs in gastric cancer. Gastroenterol Res Pract. 2014;2014:1–9. https://doi.org/10.1155/2014/473817.

    CAS  Article  Google Scholar 

  25. 25.

    Parafioriti A, et al. Ewing’s sarcoma: an analysis of miRNA expression profiles and target genes in paraffin-embedded primary tumor tissue. Int J Mol Sci. 2016;17(5):656.

    Article  Google Scholar 

  26. 26.

    Hafner M, Landgraf P, Ludwig J, Rice A, Ojo T, Lin C, et al. Identification of microRNAs and other small regulatory RNAs using cDNA library sequencing. Methods. 2008;44(1):3–12. https://doi.org/10.1016/j.ymeth.2007.09.009.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9(4):357–9. https://doi.org/10.1038/nmeth.1923.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Anders S, Pyl PT, Huber W. HTSeq--a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31(2):166–9. https://doi.org/10.1093/bioinformatics/btu638.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550. https://doi.org/10.1186/s13059-014-0550-8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  30. 30.

    Kurozumi S, Yamaguchi Y, Kurosumi M, Ohira M, Matsumoto H, Horiguchi J. Recent trends in microRNA research into breast cancer with particular focus on the associations between microRNAs and intrinsic subtypes. J Hum Genet. 2017;62(1):15–24. https://doi.org/10.1038/jhg.2016.89.

    CAS  Article  PubMed  Google Scholar 

  31. 31.

    Gupta I, Sareyeldin R, al-Hashimi I, al-Thawadi HA, al Farsi H, Vranic S, et al. Triple negative breast cancer profile, from gene to microRNA, in relation to ethnicity. Cancers. 2019;11(3):363. https://doi.org/10.3390/cancers11030363.

    CAS  Article  PubMed Central  Google Scholar 

  32. 32.

    Khordadmehr M, Shahbazi R, Ezzati H, Jigari-Asl F, Sadreddini S, Baradaran B. Key microRNAs in the biology of breast cancer; emerging evidence in the last decade. J Cell Physiol. 2019;234(6):8316–26. https://doi.org/10.1002/jcp.27716.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Gebeshuber CA, Martinez J. miR-100 suppresses IGF2 and inhibits breast tumorigenesis by interfering with proliferation and survival signaling. Oncogene. 2013;32(27):3306–10. https://doi.org/10.1038/onc.2012.372.

    CAS  Article  PubMed  Google Scholar 

  34. 34.

    Li XY, Luo QF, Wei CK, Li DF, Li J, Fang L. MiRNA-107 inhibits proliferation and migration by targeting CDK8 in breast cancer. Int J Clin Exp Med. 2014;7(1):32–40.

    PubMed  PubMed Central  Google Scholar 

  35. 35.

    Wang B, Wang H, Yang Z. MiR-122 inhibits cell proliferation and tumorigenesis of breast cancer by targeting IGF1R. PLoS One. 2012;7(10):e47053. https://doi.org/10.1371/journal.pone.0047053.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  36. 36.

    Wang Y-W, Zhang W, Ma R. Bioinformatic identification of chemoresistance-associated microRNAs in breast cancer based on microarray data. Oncol Rep. 2018;39(3):1003–10. https://doi.org/10.3892/or.2018.6205.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  37. 37.

    Mei JW, Yang ZY, Xiang HG, Bao R, Ye YY, Ren T, et al. MicroRNA-1275 inhibits cell migration and invasion in gastric cancer by regulating vimentin and E-cadherin via JAZF1. BMC Cancer. 2019;19(1):740. https://doi.org/10.1186/s12885-019-5929-1.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  38. 38.

    Sun KY, Peng T, Chen Z, Huang J, Zhou XH. MicroRNA-1275 suppresses cell growth, and retards G1/S transition in human nasopharyngeal carcinoma by down-regulation of HOXB5. J Cell Commun Signal. 2016;10(4):305–14. https://doi.org/10.1007/s12079-016-0351-9.

    Article  PubMed  PubMed Central  Google Scholar 

  39. 39.

    He J, Yu L, Wang CM, Zhou XF. MiR-1275 promotes non-small cell lung cancer cell proliferation and metastasis by regulating LZTS3 expression. Eur Rev Med Pharmacol Sci. 2018;22(9):2680–7. https://doi.org/10.26355/eurrev_201805_14964.

    CAS  Article  PubMed  Google Scholar 

  40. 40.

    Liu MD, Wu H, Wang S, Pang P, Jin S, Sun CF, et al. MiR-1275 promotes cell migration, invasion and proliferation in squamous cell carcinoma of head and neck via up-regulating IGF-1R and CCR7. Gene. 2018;646:1–7. https://doi.org/10.1016/j.gene.2017.12.049.

    CAS  Article  PubMed  Google Scholar 

  41. 41.

    Yang JR, Shi MX, Zeng Y. LncRNA HAND2-AS1 inhibits proliferation and promotes apoptosis of chronic myeloid leukemia cells by sponging with micRNA-1275. Eur Rev Med Pharmacol Sci. 2019;23(5):2103–11. https://doi.org/10.26355/eurrev_201903_17254.

    Article  PubMed  Google Scholar 

  42. 42.

    Kuo Y-C, Huang KY, Yang CH, Yang YS, Lee WY, Chiang CW. Regulation of phosphorylation of Thr-308 of Akt, cell proliferation, and survival by the B55α regulatory subunit targeting of the protein phosphatase 2A holoenzyme to Akt. J Biol Chem. 2008;283(4):1882–92. https://doi.org/10.1074/jbc.M709585200.

    CAS  Article  PubMed  Google Scholar 

  43. 43.

    Kiely PA, O'Gorman D, Luong K, Ron D, O'Connor R. Insulin-like growth factor I controls a mutually exclusive association of RACK1 with protein phosphatase 2A and β1 integrin to promote cell migration. Mol Cell Biol. 2006;26(11):4041–51. https://doi.org/10.1128/MCB.01868-05.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  44. 44.

    Jiang H, Cheng L, Hu P, Liu R. MicroRNA-663b mediates TAM resistance in breast cancer by modulating TP73 expression. Mol Med Rep. 2018;18(1):1120–6. https://doi.org/10.3892/mmr.2018.9064.

    CAS  Article  PubMed  Google Scholar 

  45. 45.

    Friedrich M, Heimer N, Stoehr C, Steven A, Wach S, Taubert H, et al. CREB1 is affected by the microRNAs miR-22-3p, miR-26a-5p, miR-27a-3p, and miR-221-3p and correlates with adverse clinicopathological features in renal cell carcinoma. Sci Rep. 2020;10(1):6499. https://doi.org/10.1038/s41598-020-63403-y.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  46. 46.

    Tang W, Li GS, Li JD, Pan WY, Shi Q, Xiong DD, et al. The role of upregulated miR-375 expression in breast cancer: an in vitro and in silico study. Pathol Res Pract. 2020;216(1):152754. https://doi.org/10.1016/j.prp.2019.152754.

    CAS  Article  PubMed  Google Scholar 

  47. 47.

    Cohen H, Ben-Hamo R, Gidoni M, Yitzhaki I, Kozol R, Zilberberg A, et al. Shift in GATA3 functions, and GATA3 mutations, control progression and clinical presentation in breast cancer. Breast Cancer Res. 2014;16(6):464. https://doi.org/10.1186/s13058-014-0464-0.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  48. 48.

    Liu X, Li C, Yang Y, Liu X, Li R, Zhang M, et al. Synaptotagmin 7 in twist-related protein 1-mediated epithelial – Mesenchymal transition of non-small cell lung cancer. EBioMedicine. 2019;46:42–53. https://doi.org/10.1016/j.ebiom.2019.07.071.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

I wish to thank the patients and their families. I would like to express my specially thanks of gratitude to my supervisor for her frequent support and encouragement. I also offer my sincere appreciation for Al mufti and Luay laboratories which helped and supported me through collecting FFPE blocks of BC and clinical data.

Funding

This study was supported by Sevan Majed who is a Ph.D. student.

Author information

Affiliations

Authors

Contributions

Suhad Mustafa and Sevan Majed were responsible for the experimental design. They led to the execution of the experiments. Although analyses of MACE-sequencing and sRNA sequencing were done in Genxpro company, in Germany. Data analysis and bioinformatic tasks were done by Sevan Majed. Sevan Majed also discussed and interpreted the data and also did the manuscript mapping and submission, but Suhad Mustafa supervised the project. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Sevan Omer Majed.

Ethics declarations

Ethics approval and consent to participate

Prior to the study, all methods done in this study including human participants and human materials were approved by Human Research Ethics Committee (HREC) at Science College in Salahuddin University-Erbil (Approval no.4c/132). All methods while performing the study were also performed in accordance with the 1964 Helsinki Declaration and the written informed-consent and Permission for the publication were obtained from all research participants.

Consent for publication

Not applicable.

Competing interests

The author announces no conflict of interest. Mr. Sevan is a teacher of Salahaddin University-Erbil, a subsidiary of ministry of higher education in Kurdistan region government (KRG).

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Majed, S.O., Mustafa, S.A. MACE-Seq-based coding RNA and TrueQuant-based small RNA profile in breast cancer: tumor-suppressive miRNA-1275 identified as a novel marker. BMC Cancer 21, 473 (2021). https://doi.org/10.1186/s12885-021-08218-4

Download citation

Keywords

  • Breast cancer
  • miRNA
  • Small RNA and gene expression
  • miR-1275 and its target genes
  • Pathogenesis
  • And tumor suppressor