107
Views
62
CrossRef citations to date
0
Altmetric
Original Research

Downregulation of hsa_circ_0011946 suppresses the migration and invasion of the breast cancer cell line MCF-7 by targeting RFC3

, , , , &
Pages 535-544 | Published online: 19 Mar 2018

Abstract

Introduction

Although some circRNAs have been found to regulate the progression of malignancies, their functions and coupled molecular mechanisms are still unclear. In our study, we sought to assess the underlying molecular mechanisms of circRNAs in breast cancer and therefore explored the differentially expressed circRNAs and co-expression networks, followed by in vitro experiments.

Materials and methods

High-throughput RNA sequencing was performed to obtain an unbiased profile of circRNA expression. CircRNA-miRNA-mRNA co-expression networks were predicted, and sequence analyses were carried out. The MTT, transwell migration and invasion assay was conducted in Michigan Cancer Foundation-7 cells that had been transfected with si-circRNA and si-negative control (si-NC).

Results

A total of 152 circRNAs were differentially expressed in breast cancer tissues, among which 85 were upregulated and 67 downregulated. Out of these, hsa_circ_0011946 was selected and the subsequent bioinformatics analysis predicted that hsa_circ_0011946 sponging miR-26a/b directly targeted replication factor C subunit 3 (RFC3) and that its knockdown could inhibit RFC3 mRNA and protein expression. Furthermore, hsa_circ_0011946 loss-of-function significantly suppressed the migration and invasion of Michigan Cancer Foundation-7 cells.

Conclusion

Together, these results indicate that hsa_circ_0011946 and RFC3 comprise a novel pathway involved in the progression of breast cancer.

Introduction

Breast cancer is the most frequently occurring female malignancy and has a gradually increasing incidence all over the world.Citation1 Genetic mutations and epigenetic alterations make the mechanisms of its progression difficult to understand. With the development of gene chip and high-throughput sequencing technology, many studies have been focused on the identification of differentially expressed non-coding RNAs, including microRNAs, long non-coding RNAs, and circular RNAs (circRNAs), which are considered to be promising biomarkers for the early detection of breast cancer, as well as potential therapeutic targets.Citation2Citation4

CircRNAs are a class of special non-coding RNA and are characterized by a highly stable, covalently closed continuous loop, without any 5′–3′ polarity or a polyadenylated tail, giving them the distinct ability to counteract RNA exonucleolytic digestion.Citation5 In recent years, circRNA research has experienced explosive growth, and more than 20,000 circRNAs have now been identified in eukaryocytes.Citation6 Functional studies indicate that circRNAs are involved in multiple physiological and pathological processes by acting as microRNA sponges to regulate gene and protein expression.Citation7,Citation8 An increasing number of studies have demonstrated that the dysregulated expression of circRNAs is closely related to the development and progression of tumors, including breast cancer.Citation9,Citation10 Intriguingly, blood-based circRNAs can serve as novel noninvasive biomarkers for the diagnosis of cancers.Citation11 Moreover, hsa_circ_006054, hsa_circ_100219, and hsa_ circ_406697 are upregulated in breast cancer tissues, and the combination of these circRNAs has a high diagnostic value for the detection of breast cancer.Citation12 Together, these findings suggest that circRNAs play a major role in tumorigenesis. It is important to highlight that the precise mechanisms behind the involvement of circRNAs in breast cancer progression have not been completely clarified.

Replication factor C (RFC) is composed of 5 subunits and plays an important role in eukaryotes in deoxyribonucleic acid (DNA) replication, DNA damage repair, and checkpoint control during cell cycle progression.Citation13,Citation14 RFC subunit 3 (RFC3) is one of the small subunits on the RFC complex and was originally purified from HeLa cells and shown to be essential for the in vitro replication of Simian virus 40.Citation14 In recent years, the various roles of RFC3 have been increasingly illuminated in studies of malignancies, including hepatocellular carcinoma,Citation15 ovarian carcinoma,Citation16 esophageal adenocarcinoma,Citation17 colorectal cancer,Citation18 and breast cancer,Citation19,Citation20 suggesting that RFC3 may be a potential oncogenic gene involved in tumorigenesis. However, it is not clear how RFC3 regulates the progression of breast cancer when combined with circRNAs.

In the present study, we utilized high-throughput sequencing and experimental validation in vitro to uncover the differentially expressed circRNAs in breast cancer tissues and cell lines and found that hsa_circ_0011946 was significantly upregulated. Moreover, a bioinformatics analysis showed RFC3 to be a target gene of hsa_circ_0011946. Furthermore, we downregulated hsa_circ_0011946 in the Michigan Cancer Foundation-7 (MCF-7) cell line and performed cell proliferation, migration, and invasion experiments, in order to elaborate on its function.

Materials and methods

Patients and specimens

Three pairs of breast cancer and corresponding adjacent non-cancerous tissues were acquired from patients who underwent surgical operations in Sun Yat-sen University Cancer Center in 2016. The commonly used clinicopathologic features of breast cancer are given in . All patients were diagnosed with invasive ductal breast cancer and the tumor-node-metastasis (TNM) classification was in accordance with the American Joint Committee on Cancer TNM Staging. Samples were rapidly stored in liquid nitrogen for high-throughput RNA sequencing. Written informed consent was obtained from all participants prior to sample collection. The study was approved by the Ethics Committee of the First Affiliated Hospital of Xiamen University and Sun Yat-sen University Cancer Center.

Cell culture

Six breast cancer cell lines (HS-578T, T47D, MCF-7, BT549, MDA-MB-231, and SKBR-3) were purchased from the Cell Bank of China Academy of Sciences, Shanghai, China. All were cultured in Roswell Park Memorial Institute medium 1640 (Invitrogen, Carlsbad, CA, USA) containing 10% fetal bovine serum (FBS; Gibco, Grand Island, NY, USA) and 1% penicillin–streptomycin at 37°C in 5% CO2 and were plated in 6-well plates at a density of 2×105 well. After 2 days of incubation, the cells were collected via RNA isolation, and used in proliferation, migration, and invasion assays.

High-throughput RNA sequencing of circRNA

In order to build a library, high-throughput RNA sequencing was performed as previously described.Citation21 The clean reads were aligned to the reference genome with Bowtie2 (http://bowtie-bio.sourceforge.net/bowtie2/manual.shtml).Citation22 The junctions of the unmapped reads were then chosen using a back splice algorithm. Finally, circRNAs were verified using circbase (http://www.circbase.org). “Mapped backsplicing junction reads per million mapped reads” (RPM) was used to measure the expression level of each circRNA.

Kyoto encyclopedia of genes and genomes (KEGG) biological pathway and co-expression analysis

The KEGG pathway analysis was performed to determine the involvement of linear transcripts in different biological pathways, as previously described.Citation21 Differentially expressed circRNAs were used as an alignment to enrich the circRNA-miRNA-mRNA network according to the miRNA target prediction, which was based on TargetScan and miRanda, as previously described.Citation23,Citation24 Cytoscape software (V. 3.2.1) was employed to build a circRNA-miRNA-mRNA network.

Validation of circRNA and head-to-tail splicing

The complementary DNA (cDNA) was acquired as described below and genomic DNA (gDNA) was isolated using the DNeasy Blood & Tissue Kit (Qiagen, Inc., Valencia, CA, USA). Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) was performed with the SYBR Select Master Mix (Applied Biosystems) using ABI7300 System (Applied Biosystems, Foster City, CA, USA). We designed divergent primers that bound to the circRNA transcript (chr1: 41578954-41618413) formed in a 3′–5′ fashion. After gel purification using the QIAquick Gel Extraction Kit (Qiagen), the RT-qPCR product was sequenced using the Sanger method, in order to confirm the head-to-tail splicing, as previously described.Citation21,Citation25

MTT assay

The proliferation of each breast cancer cell line was monitored using an MTT Cell Proliferation/Viability Assay kit (R&D Systems, Inc., Minneapolis, MN, USA), according to the manufacturer’s protocol. Briefly, following transfection with si-hsa_circ_0011946 or si-NC, 15 µL of MTT was added to each well (5×104) and incubated for 4 h at 37°C. Then, the stop solution was added and the solution was kept overnight at room temperature. The optical density value was measured at 492 nm using a SpectraMax M5 ELISA plate reader (Molecular Devices, LLC, Sunnyvale, CA, USA).

Cell transfection

The small interfering RNAs (siRNA) utilized for cell transfection were synthesized by GenePharma Co., Ltd. (Shanghai, China) and the sequence of the functional si-hsa_circ_0011946 was ACCAAAGCATCTAGTGCTTTT. MCF-7 cells were transfected with si-hsa_circ_0011946 for 48 h at 37°C using Lipofectamine 2000 (Invitrogen), according to the manufacturer’s protocol.

RT-qPCR

The RNA was extracted using TRIzol (Invitrogen; Thermo Fisher Scientifc, Inc., Waltham, MA, USA). Moloney Murine Leukemia Virus reverse transcriptase (Promega Corporation, Madison, WI, USA) and oligo dT 15 primers (Thermo Fisher Scientific, Inc.) were used to synthesize cDNA. Divergent primers were designed to ensure the amplification of the head-to-tail splicing of circRNA using the ABI7300 System (Applied Biosystems) with the SYBR Select Master Mix (Applied Biosystems), as previously described.Citation26 Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and U6 levels were utilized to normalize the expressions of the circRNAs and RFC3 mRNA, respectively. The PCR primers used in this study were as follows: hsa_circ_0011946, forward 5′-GCTGGTGTTCCTTGACTGGA-3′ and reverse 5′-CACTGTAGCAAACCAGCATTTCT-3′; hsa_circ_0007534, forward 5′-GTGACGGAAATCCAATTGCACC-3′ and forward 5′-ATGGAATTGCTGGCGAGTTG-3′; hsa_circ_0093859, forward 5′-ACCCAATGGAGTTCTCAGCAG-3′ and forward 5′-TGGAAGGTCTGACAGGAATGTG-3′; RFC3, forward 5′-TCCCTGCTTCTGATTTCCTTTACC-3′ and forward 5′-GGCTTCCCTGACCACCCTATTTA-3′; GAPDH, forward 5′-GCACCGTCAAGCTGAGAAC-3′ and reverse 5′-TGGTGAAGACGCCAGTGGA-3′; U6, forward 5′-AGCCCGCACTCAGAACATC-3′ and reverse 5′-GCCACCAAGACAATCATCC-3′. The reaction conditions were 95°C for 2 min, 36 cycles of denaturation at 95°C for 15 s, annealing at 60°C for 30 s, and extension at 60°C for 30 s. The relative expression levels of circRNAs and RFC3 were calculated using the 2−ΔΔCq method.Citation27

Western blotting

The protein was extracted using NP-40 buffer (Beyotime Institute of Biotechnology, Haimen, China), and the concentration was determined using the Bicinchoninic Acid Kit for Protein Determination (Sigma-Aldrich, St. Louis, MO, USA). Samples containing 50 µg of protein were separated on a 10% sodium dodecyl sulfate polyacrylamide gel electrophoresis and transferred to nitrocellulose membranes (Bio-Rad Laboratories, Inc., Hercules, CA, USA). The primary antibody anti-RFC3 was purchased from Santa Cruz Biotechnology, Dallas, TX, USA (Cat. no: sc-390293; dilution: 1:1,000). The membranes were incubated with the appropriate horseradish peroxidase-conjugated secondary antibody (Cat. no: sc-516102; dilution: 1:10,000; Santa Cruz Biotechnology), following visualization with chemiluminescence (Thermo Fisher Scientific, Inc.). GAPDH primary antibody (Cat. no: 2118; dilution: 1:2,000; Cell Signaling Technology, Inc., Danvers, MA, USA) was used as the control antibody.

Migration and invasion assays

The MCF-7 cells were transfected with si-circRNA and si-NC for 48 h and cells were resuspended in 5% FBS medium to achieve a density of 1×106 cells/mL. A total of 100 µL cell suspension medium with 5% FBS was added to the upper chamber, whereupon 600 µL complete medium was added to the bottom chamber and incubated at 37°C with 5% CO2. For the transwell migration assays, the cells on the upper surface of the membrane were removed with cotton swabs, and those on the lower surface were the migrated cells. After being fixed with 4% paraformaldehyde and stained with 0.1% crystal violet solution, the cells that passed through the filter were photographed by inverted fluorescence microscope (Leica Microsystems GmbH, Wetzlar, Germany). The transwell invasion assay was carried out as described above and as previously described.Citation28

Statistical analysis

SPSS (version 15.0; SPSS, Inc., Chicago, IL, USA) and GraphPad Prism (version 6.0; GraphPad Software, Inc., La Jolla, CA, USA) were used for the general statistical analysis. The student’s t-test was used to analyze 2-group differences. Inter-group differences were analyzed by one-way analysis of variance, followed by a post hoc Tukey test for multiple comparisons. P<0.05 was considered statistically significant.

Results

Differential expression of circRNAs in breast tumor and normal tissues

Differentially expressed circRNAs were selected by false discovery rate ≤0.001 and fold change ≥2 or fold change ≤0.5. We found that 152 circRNAs were differentially expressed between the breast cancer tissues and the corresponding adjacent non-cancerous tissues, among which 85 were upregulated and 67 downregulated (). In order to analyze the features and function of differentially expressed circRNAs and their host genes, KEGG signaling pathway analysis was performed using the host genes, which showed a strong correlation with PI3K/AKT, focal adhesion, and protein digestion, and the absorption signaling pathway. The PI3K/AKT signaling pathway was the top pathway for the host genes (). Our findings suggest that these pathways might contribute to the pathogenesis and development of breast cancer.

Figure 1 The changes in the circRNA expression profile in breast cancer.

Notes: Differentially expressed circRNAs were selected by false discovery rate ≤0.001 and fold change ≥2 or fold change ≤0.5 and that 152 circRNAs were differentially expressed between breast cancer tissues and the corresponding adjacent non-cancerous tissues. Eighty-five circRNAs were upregulated and 67 were downregulated using high-throughput RNA sequencing technology (A). A KEGG pathway analysis was performed to determine the correlation between differentially expressed circRNA and pathways in breast cancer (B).

Abbreviations: DEG, differentially expressed gene; CR, cancer tissue; KEGG, Kyoto encyclopedia of genes and genomes; MAPK, mitogen activated protein kinase; NOD, nucleotide oligomerization domain; PR, paraneplastic normal tissue; TGF, transforming growth factor; TNF, tumor necrosis factor; TPM, transcripts per Kilobase Million.

Figure 1 The changes in the circRNA expression profile in breast cancer.Notes: Differentially expressed circRNAs were selected by false discovery rate ≤0.001 and fold change ≥2 or fold change ≤0.5 and that 152 circRNAs were differentially expressed between breast cancer tissues and the corresponding adjacent non-cancerous tissues. Eighty-five circRNAs were upregulated and 67 were downregulated using high-throughput RNA sequencing technology (A). A KEGG pathway analysis was performed to determine the correlation between differentially expressed circRNA and pathways in breast cancer (B).Abbreviations: DEG, differentially expressed gene; CR, cancer tissue; KEGG, Kyoto encyclopedia of genes and genomes; MAPK, mitogen activated protein kinase; NOD, nucleotide oligomerization domain; PR, paraneplastic normal tissue; TGF, transforming growth factor; TNF, tumor necrosis factor; TPM, transcripts per Kilobase Million.

Construction of the circRNA-miRNA-mRNA interaction network

We found a network consisting of 13 circRNAs and its target genes. All the circRNAs are capable of indirectly regulating RFC3 by interacting with miR-26a, miR-26b, miR-200b, miR-200c, or miR-429 (). This indicates the tight correlation between RFC3 and breast cancer progression.

Figure 2 Three circRNAs (hsa_circ_0093859, hsa_circ_0007534, and hsa_circ_0011946) were up-regulated in both cells and breast tissues.

Notes: The circRNA-miRNA-mRNA networks were predicted based on online databases (A), the sexangle denotes miRNA, the inverted triangle represents circRNA, the rhombus is mRNA, and their relationship is depicted by the linear connections. Divergent primers were used to amplify circRNAs in cDNA but not gDNA (B). Sanger sequencing was performed to confirm the head-to-tail back-spliced site (red arrow) (C). The levels of hsa_circ_0093859, hsa_circ_0007534, and hsa_circ_0011946 expression were carried out using RT-qPCR among 6 different breast cancer cell lines (HS-578T, T47D, MCF-7, BT549, MDA-MB-231, and SKBR-3) (D) and 12 pairs of breast cancer tissues (E).

Abbreviations: cDNA, complementary DNA, gDNA, genome DNA; RFC, replication factor C; RFC3, RFC subunit 3; RT-qPCR, reverse transcription-quantitative polymerase chain reaction.

Figure 2 Three circRNAs (hsa_circ_0093859, hsa_circ_0007534, and hsa_circ_0011946) were up-regulated in both cells and breast tissues.Notes: The circRNA-miRNA-mRNA networks were predicted based on online databases (A), the sexangle denotes miRNA, the inverted triangle represents circRNA, the rhombus is mRNA, and their relationship is depicted by the linear connections. Divergent primers were used to amplify circRNAs in cDNA but not gDNA (B). Sanger sequencing was performed to confirm the head-to-tail back-spliced site (red arrow) (C). The levels of hsa_circ_0093859, hsa_circ_0007534, and hsa_circ_0011946 expression were carried out using RT-qPCR among 6 different breast cancer cell lines (HS-578T, T47D, MCF-7, BT549, MDA-MB-231, and SKBR-3) (D) and 12 pairs of breast cancer tissues (E).Abbreviations: cDNA, complementary DNA, gDNA, genome DNA; RFC, replication factor C; RFC3, RFC subunit 3; RT-qPCR, reverse transcription-quantitative polymerase chain reaction.

Verification of circRNA expression by RT-qPCR

We selected 3 circRNAs (hsa_circ_0093859, hsa_circ_ 0007534, and hsa_circ_0011946) on which to focus in our study. Divergent primers were designed for use in detecting the circRNAs. Our results indicated that the circRNAs were only amplified using the divergent primers and cDNA as templates; however, the amplified bands showed no signs of gDNA (). This confirmed that a head-to-tail splicing existed in these circRNAs. In order to further investigate the head-to-tail back-spliced site, Sanger sequencing was performed and subsequently validated the back-spliced junction sites labeled by red arrows, as shown in . Based on the above observations, the hsa_circ_0093859, hsa_circ_0007534, and hsa_circ_0011946 expression levels were measured in 6 different breast cancer cell lines (HS-578T, T47D, MCF-7, BT549, MDA-MB-231, and SKBR-3). We noted that hsa_circ_0093859, hsa_circ_0007534, and hsa_circ_0011946 had an obviously higher expression in the MCF-7 cell line than in other cell lines. More importantly, hsa_circ_0011946 was stably expressed in the majority of cell lines, except for MDA-MB-231 (). Therefore, we finally focused on hsa_circ_0011946 and the MCF-7 cell line in our study. In addition, we also measured the relative expressions of hsa_circ_0093859, hsa_circ_0007534, and hsa_circ_0011946 circRNAs in 12 pairs of breast cancer patient specimens (). The results showed that these 3 circRNAs were upregulated in breast cancer tissues.

Downregulated hsa_circ_0011946 inhibited RFC3 expression and migration, and invasion in MCF-7 cells

To distinguish the effects of hsa_circ_0011946 on RFC3 expression, we examined the effects of its downregulation on RFC3 expression in MCF-7 cells transfected with hsa_circ_0011946 inhibitors. As expected, the downregulation of hsa_circ_0011946 significantly inhibited the mRNA expression of RFC3 in MCF-7 cells compared with the NC group (). In addition, the protein expression of RFC3 was also suppressed in MCF-7 cells transfected with hsa_circ_0011946 inhibitors, when compared with the control group ().

Figure 3 Downregulated hsa_circ_0011946 inhibited RFC3 expression in MCF-7 cells.

Notes: After transfection with si-hsa_circ_0011946 for 48 h, the expression levels of hsa_circ_0011946 and RFC3 were assessed with a RT-qPCR assay (A). The protein expression of RFC3 was measured by Western blotting (B). n=3 in each group. *P<0.05, **P<0.001 compared with the NC group.

Abbreviations: GAPDH, glyceraldehyde-3-phosphate dehydrogenase; NC, negative control; RFC, replication factor C; RFC3, RFC subunit 3; RT-qPCR, reverse transcription-quantitative polymerase chain reaction.

Figure 3 Downregulated hsa_circ_0011946 inhibited RFC3 expression in MCF-7 cells.Notes: After transfection with si-hsa_circ_0011946 for 48 h, the expression levels of hsa_circ_0011946 and RFC3 were assessed with a RT-qPCR assay (A). The protein expression of RFC3 was measured by Western blotting (B). n=3 in each group. *P<0.05, **P<0.001 compared with the NC group.Abbreviations: GAPDH, glyceraldehyde-3-phosphate dehydrogenase; NC, negative control; RFC, replication factor C; RFC3, RFC subunit 3; RT-qPCR, reverse transcription-quantitative polymerase chain reaction.

In order to assess the potential effects of hsa_circ_0011946 on cell proliferation, an MTT assay was performed after si-hsa_circ_0011946 transfection at 24, 48, or 72 h. Nevertheless, no obvious differences were found in proliferation between the siRNA-NC and si-hsa_circ_0011946 groups at any time point (). To investigate whether the migration and invasion of MCF-7 cells was regulated by hsa_circ_0011946, si-hsa_circ_0011946, or siRNA-NC, transfected MCF-7 cells were cultured in a transwell plate. After 24-h incubation, the number of migrated cells in the si-hsa_circ_0011946 group was significantly lower than that in the siRNA-NC group. By using a Boyden chamber coated with matrigel, we then measured the effect of si-hsa_circ_0011946 on MCF-7 cell invasion after 24-h incubation. Compared with the siRNA-NC group, the si-hsa_circ_0011946 transfected cells showed a significant reduction in the number of invasive cells (). These findings suggest that hsa_circ_0011946 may be an oncogene and that its loss-of-function can inhibit the migration and invasion of MCF-7 cells.

Figure 4 The downregulation of hsa_circ_0011946 inhibited the migration and invasion of MCF-7 cells.

Notes: Cell proliferation was monitored using an MTT assay after transfection with si-hsa_circ_0011946 for 24, 48, and 72 h (A). After transfection with si-hsa_circ_0011946 for 48 h, MCF-7 cells showed reduced migration and invasion abilities (B). n=3 in each group. *P<0.05, **P<0.001 compared with the NC group.

Abbreviations: MTT, 3-(4, 5-dimethylthiazol-2-yl)-2, 5-diphenyltetrazolium bromide; NC, negative control; OD, optical density.

Figure 4 The downregulation of hsa_circ_0011946 inhibited the migration and invasion of MCF-7 cells.Notes: Cell proliferation was monitored using an MTT assay after transfection with si-hsa_circ_0011946 for 24, 48, and 72 h (A). After transfection with si-hsa_circ_0011946 for 48 h, MCF-7 cells showed reduced migration and invasion abilities (B). n=3 in each group. *P<0.05, **P<0.001 compared with the NC group.Abbreviations: MTT, 3-(4, 5-dimethylthiazol-2-yl)-2, 5-diphenyltetrazolium bromide; NC, negative control; OD, optical density.

Discussion

CircRNAs perform a wide variety of biological functions in eukaryotic cells by competing with endogenous RNAs (ceRNAs) or acting as miRNA sponges, interacting with RNA binding proteins, modulating the stability of mRNAs, regulating gene transcription, or translating proteins.Citation8 As a novel gene regulator, circRNAs are dysregulated in multiple pathological processes, including malignancies.Citation11,Citation12 Currently, the clinical diagnostic value of circRNAs has begun to emerge for breast cancer.Citation12,Citation25 Moreover, the overexpression of circ-ABCB10 and circDENND4C has been shown to promote the proliferation of breast cancer cells in vitro.Citation9,Citation10 In the present study, we performed the circRNA expression profiling analysis in breast cancer tissues and in the corresponding adjacent non-cancerous tissues using high-throughput sequencing and found that 85 of 152 differentially expressed circRNAs were significantly upregulated and 67 significantly downregulated in breast cancer tissues. Specifically, the bioinformatics-based circRNA-miRNA-mRNA interaction network revealed that the hsa_circ_0011946/RFC3 signaling pathway might be associated with the progression of breast cancer. Downregulated hsa_circ_0011946 could inhibit the migration and invasion of breast cancer cells, but the underlying mechanism was mediated, at least partially, through the suppression of RFC3 mRNA and protein expression.

Hsa_circ_0011946 is located in chr1: 41578954-41618413, the spliced sequence length is 782 ribonucleotides, and its associated-gene symbol is scm polycomb group protein homolog 1 (SCMH1; circBase, http://www.circbase.org/). SCMH1 encodes a protein with several characteristic domains, including the malignant brain tumor domains, which directly interact with histone methylation, H3K4, H3K9, H3K27, H3K36, and H4K20.Citation29,Citation30 Previous research has confirmed that H3K4 methylation is increased in breast tumors and is driven by an activated PI3K/AKT signaling pathway and inactivated demethylase KDM5A.Citation31 Surprisingly, the KEGG signaling pathway analysis suggested that PI3K/AKT signaling was most closely linked to the abnormal expression of circRNAs in breast cancer tissues. All these findings support the hypothesis that SCMH1-originated hsa_circ_0011946 regulates the development of breast cancer by PI3K/AKT-mediated histone methylation.

Accumulating evidence has demonstrated that the main function of circRNA is regulating miRNA expression by acting as a ceRNA.Citation32Citation34 For example, hsa_circ_0006528 can inhibit miR-7-5p expression in adriamycin-resistant MCF-7 breast cancer cells.Citation35 Moreover, circ-ABCB10 acts as a sponge for miR-1271 in breast cancer cells.Citation10 Although the fact that circRNA harbors miRNA-binding sites has not been experimentally verified in our study, we did establish circRNA-miRNA-mRNA co-expression networks and found that the top 3 circRNAs (hsa_circ_0011946, hsa_circ_0007534, and hsa_circ_0093859) were related to 5 miRNAs (miR-26a, miR-26b, miR-200b, miR-200c, and miR-429). In MDA-MB-231 cells, the overexpression of miR-26a/b led to a marked decrease in the proliferative, migratory, and invasive capacity of cells.Citation36 MiR-200b/c was also shown to inhibit tumor growth and metastasis in MCF-7 and MDA-MB-231 breast cancer cells.Citation37,Citation38 MiR-429 is a member of the miR-200 family, which is associated with lymph node metastasis in breast cancer patients and the inhibition of migration and invasion in breast cancer cells.Citation39,Citation40 From these studies, we deduced that hsa_circ_0011946, as a miRNA sponge, might regulate the progression of breast cancer by binding up miR-26a/b, miR-200b/c, and miR-429, which would otherwise suppress their expression and reverse their post-translational regulatory mechanism. Unexpectedly, all these miRNAs had a mutual target gene, RFC3, as confirmed by the online prediction software miRanda (http://www.microRNA.org) and TargetScan (http://www.targetscan.org).

Previous studies indicated that the overexpression of RFC3 promotes breast cancer metastasis and MCF-7 cell proliferation, and is associated with poor prognosis via the epithelial–mesenchymal transition, while the knockdown of RFC3 expression can inhibit MCF-7 cell proliferation and invasion.Citation19,Citation20 Consistent with these conclusions, our findings showed that the downregulation of RFC3 levels with small interfering hsa_circ_0011946 resulted in the inhibition of migration and invasion in MCF-7 cells. Breast cancer patients with high RFC3 levels have a higher risk of distant metastases,Citation19 thus RFC3 might be an oncogene and might serve as a therapeutic target for breast cancer.

Conclusion

Taken together, these circRNA expression profile findings from high-throughput RNA sequencing technology showed 152 circRNAs changing remarkably in breast cancer tissues and the preliminarily determined hsa_circ_0011946 being a key regulator of breast cancer. The outcome of the biomathematical prediction and in vitro experiments manifested that the inactivation of the hsa_circ_0011946/RFC3 signaling pathway could inhibit the migration and invasion capacities of MCF-7 cells.

Acknowledgments

This work was supported by grants from the National Natural Science Foundation of China (No. 81402527), the Natural Science Foundation of Fujian Province (No. 2015J01550, 2016J01635), the Science and Technology Planning Projects of Xiamen Science and Technology Bureau (No. 3502Z20174070), and the Guangdong Medical Research Foundation (No. A2017023).

Supplementary material

Table S1 Clinical characteristics of the study population

Disclosure

The authors report no conflicts of interest in this work.

References

  • DeSantisCMaJBryanLJemalABreast cancer statistics, 2013CA Cancer J Clin2014641526224114568
  • LuanTZhangXWangSLong non-coding RNA MIAT promotes breast cancer progression and functions as ceRNA to regulate DUSP7 expression by sponging miR-155-5pOncotarget2017844761537616429100300
  • MuluhngwiPAlizadeh-RadNVittitowSLKalbfleischTSKlingeCMThe miR-29 transcriptome in endocrine-sensitive and resistant breast cancer cellsSci Rep201771520528701793
  • ZhangHDJiangLHSunDWHouJCJiZLCircRNA: a novel type of biomarker for cancerBreast Cancer20182511728721656
  • HentzeMWPreissTCircular RNAs: splicing’s enigma variationsEMBO J201332792392523463100
  • GlažarPPapavasileiouPRajewskyNcircBase: a database for circular RNAsRNA201420111666167025234927
  • HansenTBJensenTIClausenBHNatural RNA circles function as efficient microRNA spongesNature2013495744138438823446346
  • MengSZhouHFengZCircRNA: functions and properties of a novel potential biomarker for cancerMol Cancer20171619428535767
  • LiangGLiuZTanLSuANJiangWGGongCHIF1α-associated circDENND4C promotes proliferation of breast cancer cells in hypoxic environmentAnticancer Res20173784337434328739726
  • LiangHFZhangXZLiuBGJiaGTLiWLCircular RNA circ-ABCB10 promotes breast cancer proliferation and progression through sponging miR-1271Am J Cancer Res2017771566157628744405
  • ChenSLiTZhaoQXiaoBGuoJUsing circular RNA hsa_ circ_0000190 as a new biomarker in the diagnosis of gastric cancerClin Chim Acta201746616717128130019
  • LüLSunJShiPIdentification of circular RNAs as a promising new class of diagnostic biomarkers for human breast cancerOncotarget2017827440964410728484086
  • ShimadaMOkuzakiDTanakaSReplication factor C3 of Schizosaccharomyces pombe, a small subunit of replication factor C complex, plays a role in both replication and damage checkpointsMol Biol Cell199910123991400310588638
  • XiaSXiaoLGannonPLiXRFC3 regulates cell proliferation and pathogen resistance in ArabidopsisPlant Signal Behav20105216817020023430
  • YaoZHuKHuangHShRNA-mediated silencing of the RFC3 gene suppresses hepatocellular carcinoma cell proliferationInt J Mol Med20153651393139926397132
  • ShenHXuJZhaoSShiHYaoSJiangNShRNA-mediated silencing of the RFC3 gene suppress ovarian tumor cells proliferationInt J Clin Exp Pathol2015888968897526464638
  • LockwoodWWThuKLLinLIntegrative genomics identified RFC3 as an amplified candidate oncogene in esophageal adenocarcinomaClin Cancer Res20121871936194622328562
  • KimYRSongSYKimSSAnCHLeeSHYooNJMutational and expressional analysis of RFC3, a clamp loader in DNA replication, in gastric and colorectal cancersHum Pathol201041101431143720573375
  • HeZYWuSGPengFUp-regulation of RFC3 promotes triple negative breast cancer metastasis and is associated with poor prognosis Via EMTTransl Oncol20171011927888707
  • MaengSKimGJChoiEJYangHOLeeDSSohnYC9-Cis-retinoic acid induces growth inhibition in retinoid-sensitive breast cancer and sea urchinembryonic cells via retinoid X receptor α and replication factor C3Mol Endocrinol201226111821183522949521
  • LiLGuoJChenYChangCXuCComprehensive CircRNA expression profile and selection of key CircRNAs during priming phase of rat liver regenerationBMC Genomics20171818028086788
  • MemczakSJensMElefsiniotiACircular RNAs are a large class of animal RNAs with regulatory potencyNature2013495744133333823446348
  • EnrightAJJohnBGaulUTuschlTSanderCMarksDSMicroRNA targets in DrosophilaGenome Biol200351R114709173
  • PasquinelliAEMicroRNAs and their targets: recognition, regulation and an emerging reciprocal relationshipNat Rev Genet201213427128222411466
  • NairAANiuNTangXCircular RNAs and their associations with breast cancer subtypesOncotarget2016749809678097927829232
  • XiaWQiuMChenRCircular RNA has_circ_0067934 is upregulated in esophageal squamous cell carcinoma and promoted proliferationSci Rep201663557627752108
  • LivakKJSchmittgenTDAnalysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) methodMethods200125440240811846609
  • LiJGuoYDuanLAKR1B10 promotes breast cancer cell migration and invasion via activation of ERK signalingOncotarget2017820336943370328402270
  • YasunagaSOhtsuboMOhnoYScmh1 has E3 ubiquitin ligase activity for geminin and histone H2A and regulates gemininstability directly or indirectly via transcriptional repression of Hoxa9 and Hoxb4Mol Cell Biol201333464466023207902
  • GrimmCde Ayala AlonsoAGRybinVStructural and functional analyses of methyl-lysine binding by the malignant brain tumour repeatprotein sex comb on midlegEMBO Rep20078111031103717932512
  • SpangleJMDreijerinkKMGronerACI3K/AKT signaling regulates H3K4 methylation in breast cancerCell Rep201615122692270427292631
  • HuangXYHuangZLXuYHComprehensive circular RNA profiling reveals the regulatory role of the circRNA-100338/miR-141-3p pathway in hepatitis B-related hepatocellular carcinomaSci Rep201771542828710406
  • YangZGAwanFMDuWWThe circular RNA interacts with STAT3, increasing its nuclear translocation and wound repair by modulating Dnmt3a and miR-17 functionMol Ther20172592062207428676341
  • HanDLiJWangHCircular RNA circMTO1 acts as the sponge of microRNA-9 to suppress hepatocellular carcinoma progressionHepatology20176641151116428520103
  • GaoDZhangXLiuBScreening circular RNA related to chemotherapeutic resistance in breast cancerEpigenomics2017991175118828803498
  • MaXDongWSuZFunctional roles of sialylation in breast cancer progression through miR-26a/26b targetingST8SIA4Cell Death Dis2016712e256128032858
  • ZhengQCuiXZhangDmiR-200b inhibits proliferation and metastasis of breast cancer by targeting fucosyltransferaseIV and α1, 3-fucosylated glycansOncogenesis201767e35828692034
  • SongCLiuLZPeiXQmiR-200c inhibits breast cancer proliferation by targeting KRASOncotarget2015633349683497826392416
  • XuFHeHHuangWDecreased expression of MicroRNA-200 family in human breast cancer is associated with lymph node metastasisClin Transl Oncol201618328328826201425
  • YeZBMaGZhaoYHmiR-429 inhibits migration and invasion of breast cancer cells in vitroInt J Oncol201546253153825405387