Effects of CORO2A on Cell Migration and Proliferation and Its Potential Regulatory Network in Breast Cancer

Coronin 2A (CORO2A) is a novel component of the N-CoR (nuclear receptor co-repressor) complex. Abnormal CORO2A expression is associated with carcinogenesis. We used databases from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO), and analyzed CORO2A expression and gene regulation networks in breast cancer. Expression was analyzed using GEO and TCGA database and further validated in breast cancer samples collected in our clinic. The prognostic value of CORO2A was explored by using the Kaplan–Meier survival analysis and Cox proportional hazards regression analysis. LinkedOmics was used to identify coexpressed genes associated with CORO2A. After analyzing the intersection of coexpressed genes correlated with CORO2A and differentially expressed genes after CORO2A silencing, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses of the intersecting genes were conducted by using FunRich software. Transwell assays were performed in breast cancer cells to determine the effect of CORO2A on cell migration. MTS, colony formation, and cell cycle distribution assays were performed in breast cancer cells to determine the effect of CORO2A on cell proliferation. Gene enrichment analysis was employed to explore the target networks of transcription factors and miRNAs. We found that CORO2A was upregulated and that the elevated expression of CORO2A was associated with poor overall survival (OS) and relapse-free survival (RFS) in TNBC patients. Further bioinformatics analysis of public sequencing data and our own RNA-Seq data revealed that CORO2A was probably involved in the epithelial-to-mesenchymal transition process and might have a significant effect on the migration of breast cancer cells, which might be mediated via pathways involving several miRNAs and MYC transcription factors. Functionally, the knockdown of CORO2A inhibited cell migration, decreased viability, and colony formation and induced cell cycle arrest in the G0/G1 phase in breast cancer cells. These results demonstrate that bioinformatics-based analysis efficiently reveals information about CORO2A expression and its potential regulatory networks in breast cancer, laying a foundation for further mechanistic research on the role of CORO2A in carcinogenesis. Moreover, CORO2A promotes the migration and proliferation of breast cancer cells and may have an important function in breast cancer progression. CORO2A is a potential prognostic predictor for TNBC patients. Targeting CORO2A may provide promising therapy strategies for breast cancer treatment.

Coronin 2A (CORO2A) is a novel component of the N-CoR (nuclear receptor co-repressor) complex. Abnormal CORO2A expression is associated with carcinogenesis. We used databases from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO), and analyzed CORO2A expression and gene regulation networks in breast cancer. Expression was analyzed using GEO and TCGA database and further validated in breast cancer samples collected in our clinic. The prognostic value of CORO2A was explored by using the Kaplan-Meier survival analysis and Cox proportional hazards regression analysis. LinkedOmics was used to identify coexpressed genes associated with CORO2A. After analyzing the intersection of coexpressed genes correlated with CORO2A and differentially expressed genes after CORO2A silencing, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses of the intersecting genes were conducted by using FunRich software. Transwell assays were performed in breast cancer cells to determine the effect of CORO2A on cell migration. MTS, colony formation, and cell cycle distribution assays were performed in breast cancer cells to determine the effect of CORO2A on cell proliferation. Gene enrichment analysis was employed to explore the target networks of transcription factors and miRNAs. We found that CORO2A was upregulated and that the elevated expression of CORO2A was associated with poor overall survival (OS) and relapse-free survival (RFS) in TNBC patients. Further bioinformatics analysis of public sequencing data and our own RNA-Seq data revealed that CORO2A was probably involved in the epithelial-to-mesenchymal transition process and might have a significant effect on the migration of breast cancer cells, which might be mediated via pathways involving several miRNAs and MYC transcription factors. Functionally, the knockdown of CORO2A inhibited cell migration, decreased viability, and colony formation and induced cell cycle arrest in the G0/G1 phase in breast cancer cells. These results demonstrate that bioinformatics-based analysis efficiently reveals information about CORO2A expression and its potential regulatory networks in breast cancer, laying a foundation for further mechanistic research on the role of CORO2A in carcinogenesis. Moreover, CORO2A promotes the migration and proliferation of breast cancer cells and may have an important function in breast cancer progression. CORO2A is a potential prognostic predictor for TNBC patients. Targeting CORO2A may provide promising therapy strategies for breast cancer treatment.

INTRODUCTION
Breast cancer is reported to be the most common malignancy in women, with an estimated 268,600 new cases and 41,760 deaths in 2019, which represents 30% of all new cancer cases and 15% of cancer-related deaths among American women (1,2). Breast cancer is a complex disease with cumulative genetic and epigenetic aberrations as well as cancer cell heterogeneity (3). Understanding these pathological molecular alterations has significant implications in cancer initiation and progression as well as in cancer prevention and treatment. Therefore, the identification of the potential key molecules involved in the development and progression of breast cancer, which may function as biomarkers for the diagnosis and prognosis of breast cancer, is urgently needed.
In our previous study, by using four public datasets from the Gene Expression Omnibus (GEO) and TCGA, we found that CORO2A was overexpressed in malignant breast tissues compared to normal breast tissues (4). Coronin 2A (CORO2A), also known as CLIPINB, IR10, or WDR2, is a member of the "short" coronin subfamily containing a single WD40repeat domain, which adopts the fold of a seven-bladed βpropeller (5). Compared to other well-characterized coronins, less information is available about CORO2A. In 2003, Yoon et al. (6) revealed that CORO2A is a novel component of the N-CoR (nuclear receptor co-repressor) complex. In addition, Huang et al. (7) uncovered that CORO2A mediates actindependent derepression of inflammatory response genes, which is regulated by phosphorylation and by nuclear receptor signaling pathways. Rastetter et al. (8) illuminated that the expression level of CORO2A is related to the malignant progression of colon carcinoma and that the overexpression of CORO2A increases tumor cell migration. Recently, Fridley et al. (9) observed shorter progression-free survival (PFS) in endometrioid carcinoma and clear cell carcinoma patients with high gene expression for CORO2A. Based on these studies, we studied the function and prognostic significance of CORO2A in breast cancer.
In the current study, we found that patients with elevated expression of CORO2A had poorer overall survival (OS) and relapse-free survival (RFS) in TNBC via Kaplan-Meier survival analysis. Loss-of-function experiments further indicated that knockdown of CORO2A reduces cell migration, cell viability, and colony-formation and induces cell cycle arrest in the G0/G1 phase of breast cancer cells. In addition, we also studied CORO2A expression in data from patients with breast cancer in The Cancer Genome Atlas (TCGA) and various public databases. Using multidimensional analysis methods, we analyzed the functional networks and validated the function of CORO2A in breast cancer. These findings suggest that CORO2A is a novel oncogene that may be involved in breast cancer initiation and progression, and functions as a potential biomarker for the diagnosis and prognosis of breast cancer.

MATERIALS AND METHODS
Expression Analysis of CORO2A in UALCAN UALCAN (http://ualcan.path.uab.edu.) is a user-friendly, interactive web resource for analyzing cancer transcriptome data (TCGA and MET500 transcriptome sequencing) from 31 cancer types (10). UALCAN allows users to identify biomarkers or to perform in silico validation of potential genes of interest. One of the portal's user-friendly features is that it allows analysis of relative expression of a query gene(s) across tumor and normal samples, as well as in various tumor molecular subtypes such as individual age of diagnosis, gender, stages, or other clinicopathological features, which allow us to explore the relative expression of CORO2A in breast cancer.

Bioinformatic Data Mining of TCGA Database
The breast cancer expression profile data and corresponding clinical data were downloaded from the TCGA database (http:// tcga-data.nci.nih.gov). And the clinical characteristics of the patients were shown in Table 1.

LinkedOmics Analysis
LinkedOmics is a publicly available platform for analyzing multiomics data from all 32 TCGA cancer types and mass spectrometry-based proteomics data, especially for TCGA breast, colorectal, and ovarian tumors (11). LinkedOmics is publicly available at http://www.linkedomics.org/login.php. Since the LinkFinder module of LinkedOmics allows users to search for attributes that are associated with a query attribute, such as mRNA, we investigated genes that were differentially expressed in correlation with CORO2A in the TCGA breast cancer cohort. The results were analyzed by Pearson's correlation coefficient and visualized by heat maps, volcano plots, or scatter plots. In addition, data from the LinkFinder module were ranked, and GSEA was adopted to conduct analyses of transcription factor-target enrichment and miRNA-target enrichment. The two network analyses were based on the Molecular Signatures Database (MSigDB) (12). The rank criterion was set as FDR < 0.05, and 1,000 simulations were performed.

Enrichment Analysis of Intersecting Genes
The intersecting genes of genes coexpressed with CORO2A and differentially expressed genes after CORO2A silencing were analyzed in human breast cancer. Briefly, first, the intersecting genes from downregulated genes after CORO2A silencing from RNA-seq analysis and genes positively correlated with CORO2A expression from LinkedOmics analysis were obtained. Second, the intersecting genes from upregulated genes after CORO2A silencing from RNA-seq analysis and genes negatively correlated with CORO2A expression from LinkedOmics analysis were also obtained. Next, FunRich software (FunRich_V3) was applied to perform Gene Ontology (BP and MF) and biological pathway enrichment analysis of these intersecting genes.

GeneMANIA Analysis
GeneMANIA (http://www.genemania.org) is a flexible, userfriendly web tool for the construction of a protein-protein interaction (PPI) network, for generating hypotheses about gene function, analyzing gene lists and prioritizing genes for functional assays (13). Therefore, GeneMANIA was used to visualize the gene networks of genes that were identified as being enriched in breast cancer: transcription factor MYC ( Figure 10) and miRNA-493 (Supplementary Figure 2).

Western Blotting
Western blot analysis was conducted as previously described (14). Briefly, protein extracts from harvested experimental cells were prepared using lysis buffer from Beyotime Biotechnology (Shanghai, China). Thirty micrograms of protein extracts were loaded onto a 12% SDS-PAGE gel, electrophoretically separated, and transferred to polyvinylidene membranes (Beyotime Biotechnology, Shanghai, China). After incubating the membrane with 5% non-fat milk blocking buffer at room temperature for 1 h, the membranes were incubated with primary antibody overnight at 4 • C. After incubation with a horseradish peroxidase-conjugated secondary antibody (1:2,500, Promega, Madison, WI, USA) at room temperature for 1 h, the signal was obtained by using an ECL Western blotting system (Promega, Madison, WI, USA) and visualized and quantified using the Bio-Rad ChemiDoc MP system. The primary antibodies against CORO2A (Sigma, catalog # HPA041302, 1:1,000) were obtained from Sigma (Sigma Chemical Co., St. Louis, MO, USA). β-actin, also obtained from Sigma, was used as the loading control.

Immunostaining
Immunofluorescence experiments were performed according to the instructions of the Fast ImmunoFluorescence Staining Kit (Protein Biotechnologies). The primary antibody was anti-CORO2A (Sigma, catalog # HPA041302, 1:400). The secondary antibody was Alexa Fluor 488 (Jackson ImmunoResearch). The nucleus was detected using DAPI Fluoromount-G R (Beyotime Biotechnology). The images were collected by using DeltaVision Elite Imaging System (GE Healthcare Life Sciences).

Scratch-Wound Assay
After transfection with CORO2A plasmids for 24 h, cells were seeded into a 12-well-culture plate until a monolayer of cells was formed, and then a straight line was scratched on the monolayer of cells with a sterile 1,000-µl pipette tip. Then, the culture medium was removed, and the wells were washed thrice with the culture medium to remove the debris. The width of the scratch area was measured 24 h later to estimate the migration capacity.

Transwell Assay
Cellular migration was evaluated using Transwell migration chambers precoated with a layer of Matrigel R (Corning Cat# 3422). Cells suspended in 100 µl medium without FBS were seeded in the upper chamber (3 × 10 4 cells/well) following transfection with CORO2A siRNA or control for 24 h. To the lower chamber, 600 µl medium supplemented with 10% FBS was added. After 24 h incubation at 37 • C, cells on the top of the membrane were removed with cotton swabs. The cells that migrated through the membrane were washed with PBS, fixed in 4% methanol for 20 min and stained with 0.1% crystal violet for 10 min at room temperature. The number of migratory cells was counted in five randomly selected fields under a light microscope (Leica Microsystems GmbH, Wetzlar, Germany).

Cell Viability Assay
The MTS assay was performed to detect viable cell numbers using a cell proliferation assay kit from Promega (Madison, WI, USA). For the MTS assay, MDA-MB-231 and BT549 cells were seeded in 96-well-plates at 2 × 10 3 cells/well, cultured for 24 h in growth media, and then treated with CORO2A siRNA for 24, 48, 72, and 96 h. The absorbance (A) at 490 nm was detected using a microplate reader.

Cell Colony Formation Assay
The cell colony formation assay was conducted as previously described with minor modifications (15). Briefly, 1 × 10 3 cells/well transfected with CORO2A siRNA or control for 24 h were initially seeded in a six-well-plate. The medium was refreshed every 2 days. At day 10 of treatment, when the colonies were visible by the naked eye, the cells were fixed with 4% paraformaldehyde (Sigma-Aldrich) for 30 min, stained with 1% crystal violet for 20 min and then counted. The colony numbers were quantified using ImageJ software.

Statistical Analysis
Each experiment was performed at least three times, and data are expressed as the mean ± SEM. Statistical analyses were performed using SPSS 23.0 software (SPSS, Inc., Chicago, IL, USA). One-way ANOVA was applied to determine the difference among multiple groups, and Student's t-test was adopted to compare the statistical significance between two groups. Kaplan-Meier method and log-rank test were applied to compare the survival rate. Multivariate Cox proportional hazard models were used to analyze the effect of clinical variables on patients' OS and RFS. A p < 0.05 was considered to be statistically significant.

CORO2A Expression in Breast Cancer
In a previous study, we determined that the CORO2A transcription level was consistently upregulated in multiple breast cancer cohorts from the four Gene Expression Omnibus (GEO) and TCGA databases (4). These four datasets were generated from the same detecting microarray platform: [HG-U133_Plus_2] Affymetrix Human Genome U133 Plus 2.0 Array. Therefore, there is a smaller difference between batches. The expression level of CORO2A in those Gene Expression Omnibus (GEO) datasets was visualized and shown in Figure 1A in    this study. Moreover, in the current study, further subgroup analysis of multiple clinicopathological features of breast cancer samples in TCGA consistently showed elevated level of CORO2A transcript. The level of CORO2A transcript was significantly upregulated in breast cancer tissue compared with normal tissue in subgroup analyses based on age, ethnicity, tumor stage, nodal metastatic status, disease subclass, and menopause status (Figures 1B-H). Significant overexpression of CORO2A was also validated in breast cancer samples (n = 22) collected in our clinic ( Figure 1I). Moreover, the prognostic value of CORO2A was explored by the analysis of TCGA database. Although no significant association was observed between CORO2A mRNA expression and overall survival (OS) or relapse-free survival (RFS) in breast cancer patients (Figures 2A,B). However, Kaplan-Meier survival curves of CORO2A in breast cancer patients with differential tumor subclasses revealed that elevated expression of CORO2A indicated shorter OS in triple-negative breast cancer (TNBC) patients ( Figure 2C). Moreover, the median RFS time of TNBC patients with high CORO2A mRNA expression was significantly shorter than those with low expression (Figure 2D). Unfortunately, these findings were not observed in luminal breast cancer patients (Figures 2E,F) (Tables 2, 3). Thus, these data suggest that high CORO2A expression may function as a potential diagnostic indicator in TNBC.

Analyses of the Intersection of CORO2A-coexpressed Genes and Differentially Expressed Genes After CORO2A Silencing in Human Breast Cancer
The Function module of LinkedOmics was applied to analyze mRNA sequencing data from 1,093 breast cancer patients in the TCGA. As shown in the volcano plot ( Figure 3A) To further gain insight into the functional roles of CORO2A, RNA sequencing analysis in CORO2A-knockdown BT549 cells vs. BT549 control cells was performed. We identified 765 genes that were significantly downregulated and 1,139 genes that were significantly upregulated (P < 0.05; fold change > 2) between BT549/si-CORO2A and BT549/NC cells (Figures 3D-H). After overlapping genes coexpressed with CORO2A, a total of 119 genes were identified to be downregulated after CORO2A silencing and positively correlated with CORO2A expression, while a total of 316 genes were identified to be upregulated after CORO2A silencing and negatively correlated with CORO2A expression (Figures 3G,H). This result indicated that these 435 differentially expressed genes (DEGs) were extremely likely to be involved in the functional networks of CORO2A.

Enrichment Analysis of CORO2A Functional Networks in Breast Cancer
Significant GO term analysis by FunRich showed that genes downregulated after CORO2A silencing and positively correlated with CORO2A expression participated primarily in skeletal development, electron transport, reproduction, protein transport, and cell-cell adhesion. They act as transcription factors (Figures 4A,B and Supplementary Table 3). Genes upregulated after CORO2A silencing and negatively correlated with CORO2A expression primarily participated in cell growth and/or maintenance, muscle development, signal transduction, and cell communication. They act as structural constituents of the cytoskeleton and extracellular matrix structural constituents (Figures 4C,D and  Supplementary Table 4).
It is interesting that KEGG pathway analysis of genes that were downregulated after CORO2A silencing and positively correlated with CORO2A expression showed significant enrichment in the following terms: notch receptor binds with a ligand, receptorligand binding initiates the second proteolytic cleavage of Notch receptor, and mesenchymal-to-epithelial transition pathways ( Figure 5A and Supplementary Table 3); KEGG pathway analysis of genes that were upregulated after CORO2A silencing and negatively correlated with CORO2A expression showed significant enrichment in PERK-regulated gene expression, epithelial-to-mesenchymal transition, and potassium channel pathways ( Figure 5B and Supplementary Table 4). This result indicated that CORO2A was extremely likely to be involved in the epithelial-to-mesenchymal transition process and might have a significant effect on the migration of breast cancer cells.

The Effect of CORO2A on Cell Migration in Breast Cancer Cells
As previously mentioned, we have revealed that CORO2A is probably involved in the epithelial-to-mesenchymal transition process and might have a significant effect on the cell migration of breast cancer (Figures 5A,B and Supplementary Tables 3,  4). Besides, CORO2A was upregulated and that the elevated expression of CORO2A was associated with poor overall survival (OS) and relapse-free survival (RFS) in TNBC patients. And we assessed the expression and subcellular localization of CORO2A in breast cancer cells and found that the CORO2A had a relative higher expression in TNBC than in luminal breast cancer cell lines (Figure 6A). Immunofluorescence assay showed that CORO2A was mainly accumulated in the cytoplasm of breast cancer cells (Figure 6B). TNBC is associated with poor prognosis and belongs to aggressive tumor subtype, due to aggressive tumor phenotype(s), early metastasis to brain or visceral organ, and lack of clinically established targeted therapies (16,17). Therefore, we first assessed the effect of CORO2A on cell migration in TNBC cell lines. Transwell assays were performed in MDA-MB-231 and BT549 TNBC cells. The knockdown of CORO2A at the mRNA and protein levels (Figures 7A,B) produced a decrease in migration ability ( Figure 7C) in both breast cancer cell lines. In contrast, overexpression of CORO2A at the protein level ( Figure 7B) produced an opposite effect on cell migration in these breast cancer cells (Figures 7D,E). These results indicate that CORO2A could promote the migration of breast cancer cells.

The Effect of CORO2A on Cell Growth in Breast Cancer Cells
Furthermore, we also determined the effect of CORO2A on cell growth in human breast cancer cells. MTS, colony formation, and cell cycle distribution assays were performed in MDA-MB-231 and BT549 breast cancer cells, respectively. The knockdown of CORO2A produced a decrease in viability (Figures 8A,B) and colony formation (Figures 8C,D) in both breast cancer cell lines. The cell cycle distribution analysis showed that the decreased level of CORO2A could noticeably induce cell cycle arrest in the G0/G1 phase (Figures 8E,F). MTS, colony formation, and cell cycle distribution assays (Figures 9A-C) were also performed in luminal MCF7 breast cancer cells, and the observed effects were the same as those in MDA-MB-231 and BT549 breast cancer cells. However, the overexpression of CORO2A did not produce a significantly increased effect on cell viability in MDA-MB-231 and BT549 breast cancer cells (these data were not showed). These results indicate that knockdown CORO2A could inhibit the growth of both TNBC and luminal breast cancer cells, suggesting that the decrease of CORO2A level may have "inhibitor" function in breast cancer progress.

CORO2A Networks of Transcription Factors or miRNA Targets in Breast Cancer
To further explore the targets of CORO2A in breast cancer, we analyzed the transcription factor and miRNA target networks of positively correlated gene sets generated by GSEA. The top four most significant transcription factor-target networks were related mainly to MYC, ROAZ, MYCMAX, and USF ( Table 4 and Supplementary Tables 5, 6). The miRNA-target network was associated with (ATGTACA) miR-493, (TTGCACT) miR-130a, miR-301, miR-130b, and (TTTGCAC) miR-19a, miR-19b. In addition, the protein-protein interaction network constructed by GeneMANIA further showed a correlation among genes for the TF MYC_Q2 (Figure 10) and miRNA-493 (Supplementary Figure 2). Different colors of the network edge represent different bioinformatics methods applied, such as coexpression, website prediction, colocalization, shared protein domains, and physical interactions.

DISCUSSION
CORO2A, a novel component of the N-CoR (nuclear receptor co-repressor) complex, is associated with actin-dependent derepression of inflammatory response genes (6,7). In the current study, we found that CORO2A was overexpressed in malignant breast tissues compared to normal breast tissues. Moreover, breast cancer patients with elevated expression of CORO2A had a poorer overall survival (OS) and relapse-free survival (RFS) in TNBC. Further bioinformatics analysis of public sequencing data and our own RNA-Seq data revealed that CORO2A is probably involved in the epithelial-tomesenchymal transition process and might have a significant effect on the cell migration of breast cancer (Figures 5A,B and Supplementary Tables 3, 4). In addition, we also demonstrated in vitro that knockdown of CORO2A reduces cell migration and proliferation and induces cell cycle arrest in the G0/G1 phase in breast cancer cells (Figures 7-9). To gain more detailed insights into the potential functions of CORO2A in breast cancer and its regulatory network, bioinformatics analysis of public sequencing data was performed to guide future research in breast cancer.
Our study confirmed that the CORO2A mRNA level is elevated in breast cancer compared with normal breast tissue via analysis of transcriptional sequencing data from more than 1,000 clinical samples of breast cancer in the GEO and TCGA (Figure 1). The fold change of CORO2A was similar across FIGURE 10 | Protein-protein interaction network of transcription factor MYC-target networks construction (GeneMANIA). Protein-protein interaction (PPI) network illustrates the gene set that was enriched in the target network of MYC. Different colors of the network edge represent the bioinformatics methods applied including co-expression, website prediction, co-localization, shared protein domains, and physical interactions.
the various breast cancer databases, and these results were also validated in the breast cancer samples (n = 22) collected in our clinic ( Figure 1I). This finding suggests that the overexpression of CORO2A occurs in many cases of breast cancer and deserves further clinical support as a potential diagnostic and prognostic marker.
It is interesting that KEGG pathway analysis of genes that were downregulated after CORO2A silencing and positively correlated with CORO2A expression showed significant enrichment in the following terms: notch receptor binds with a ligand, receptor-ligand binding initiates the second proteolytic cleavage of the Notch receptor, and mesenchymal-to-epithelial transition pathways ( Figure 5A and Supplementary Table 3); however, KEGG pathway analysis of genes that were upregulated after CORO2A silencing and negatively correlated with CORO2A expression showed significant enrichment in the PERK regulated gene expression, epithelial-to-mesenchymal transition, and potassium channel pathways (Figure 5B and Supplementary Table 4). This result indicated that CORO2A was probably involved in the epithelial-to-mesenchymal transition process and might have a significant effect on the migration of breast cancer cells. Thus, we validated the effects of CORO2A on cell migration and proliferation as well as on the cell cycle (Figures 7-9).
CORO2A participates in several important physiological functions, and its alteration may lead to changes of various downstream signaling pathways. GSEA enrichment analysis of target gene sets can help uncover important networks of transcription factors and miRNAs. Abnormal MYC expression is frequently involved in the initiation and development of breast cancer. MYC expression is elevated in TNBC and has been reported as one of the key features driving TNBC (18,19). Oncogenic MYC directly transcriptionally activates the IRE1/XBP1 pathway by binding to the IRE1 promoter and enhancer (20). Wang et al. (21) revealed that overexpression of lncRNA EPIC1 is associated with poor prognosis in luminal B breast cancer patients and enhances tumor growth in vitro and in vivo. Mechanistically, EPIC1 promotes cell cycle progression by interacting with MYC, and EPIC1 knockdown reduces the occupancy of MYC on its target genes. Furthermore, Lee et al. (22) demonstrated that MYC and MCL1 confer resistance to chemotherapy by expanding CSCs via mitochondrial oxidative phosphorylation (mtOXPHOS) and that targeting mitochondrial respiration and HIF-1α may reverse chemotherapy resistance in TNBC. Therefore, our analyses suggest that MYC is an important target of CORO2A and that CORO2A acts through this factor to regulate the cell migration, cell cycle, and proliferation capacity of breast cancer. Further studies should evaluate this hypothesis.
Our study identified several miRNAs that were associated with CORO2A. MicroRNAs are mainly involved in the posttranscriptional regulation of gene expression and have potential applications for the clinical assessment of patient outcome in cancer, as well as in cancer monitoring and therapy (23). The particular microRNAs in our study were associated with tumor proliferation, migration, invasion, metastasis, chemoresistance, and angiogenesis (24)(25)(26)(27). In fact, microRNA-493-5p was identified as most prominently upregulated in male breast cancer (28). MicroRNA-493 is a prognostic factor in TNBC and patients with high microRNA-493 expression have better disease-free survival than patients with low microRNA-493 expression (29). MicroRNA-493 and microRNA-19a can be used as diagnostic and prognostic markers of breast cancer (29,30). MicroRNA-493-5p, microRNA-130a, microRNA-130b-3p, microRNA-301, and microRNA-19a-3p suppress the migration and invasion of breast cancer cells (31)(32)(33)(34)(35). Dysregulation of these miRNAs is consistent with the phenotype of CORO2A upregulation in breast cancer.
This study provides multilevel evidence for the importance of CORO2A in tumorigenesis and its potential as a marker in breast cancer. CORO2A is specifically related to several tumorassociated miRNAs (such as miRNA-493) and transcription factors (such as MYC). This study uses online tools and RNAseq assays to perform target gene analyses on tumor data. Compared with traditional chip screening, this method has the advantages of large sample size, low cost, and simplicity. This enables large-scale genomics research as well as subsequent mechanistic studies of CORO2A in breast cancer. Nevertheless, whether breast cancer patients overexpressing CORO2A could benefit from the suppression of CORO2A or whether CORO2A is a promising therapeutic target still needs further experimental support, including preclinical and prospective clinical studies.

ETHICS STATEMENT
The studies involving human participants were reviewed and approved by the Ethics Committee of Xiangya Hospital. The patients/participants provided their written informed consent to participate in this study.

AUTHOR CONTRIBUTIONS
J-LD conceived, designed the research, conducted the experiments, analyzed and interpreted data, and wrote the manuscript. H-BZ, YZ, Y-HX, and YH participated in analyzing and interpreting the data. GW participated in the experimental design and provided financial support. All authors read and approved the final manuscript.