- 1Department of Gynecology, Gansu Provincial Maternity and Child-Care Hospital, Lanzhou, China
- 2Department of Pathology, Gansu Provincial Maternity and Child-Care Hospital, Lanzhou, China
- 3Department of Gynecological Oncology, Gansu Provincial Maternity and Child-Care Hospital, Lanzhou, China
Background: Microchromosome maintenance protein 10 (MCM10) is required for DNA replication in all eukaryotes, and it plays a key role in the development of many types of malignancies. However, we currently still do not know the relationship between MCM10 and ovarian cancer (OV) prognosis and immune checkpoints.
Methods: The Gene Expression Profiling Interactive Analysis and Tumor Immunology Estimation Resource (TIMER) databases were used to investigate MCM10 expression in Fan cancer. The Kaplan-Meier Plotter and PrognoScan were used to assess the relationship between MCM10 and OV prognosis. The LinkedOmics database was used to analyze the MCM10 co-expression network and explore GO term annotation and the KEGG pathway. The relationship between MCM10 expression and immune infiltration in OV was investigated using the Tumor Immunology Estimation Resource database. cBioPortal database was used to explore the relationship between MCM10 expression and 25 immune checkpoints. Finally, quantitative real-time polymerase chain reaction (qRT-PCR) was performed to detect MCM10 expression. The prognosis was also analyzed by distinguishing between high and low expression groups based on median expression values.
Results: The results of the three data sets (220,651_s_at, 222,962_s_at and 223,570_at) in KM Plotter all indicated that the overall survivalof the high MCM10 expression group was lower than that of the low expression group OV, and the results of GSE9891 also reached the same conclusion. The expression level of MCM10 was negatively correlated with B cells and CD8+T cells, and positively correlated with CD4+T Cells and Macrophages. GO term annotation and KEGG pathway analysis showed that the co-expressed genes of MCM10 were mainly enriched in cell cycle and DNA replication. The alterations in MCM10 coexisted statistically with the immune checkpoints CTLA4, TNFSF4, TNFSF18, CD80, ICOSLG, LILRB1 and CD200. PCR results displayed that MCM10 was highly expressed in OV tissues, and the increased expression of MCM10 was significantly associated with poor overall survival.
Conclusion: These results demonstrated that high expression of MCM10 was associated with poor prognosis in OV and correlated with immune checkpoints.
Introduction
Ovarian cancer (OV) is a common but fatal gynaecological malignancy. Although the mortality rate from OV has declined over the past 40 years as medical care has improved, it still remains the second leading cause of death from gynaecological cancers in women and the eighth leading cause of death in women (Irusta, 2021). The current treatment for OV is mainly surgical resection and platinum-based chemotherapy, the combination of which usually brings good outcomes, with the addition of anti-angiogenic agents usually for poorly operated and stage IV patients (Lheureux et al., 2019). Despite many efforts to treat OV, the prognosis still remains poor due to the recurrence and metastasis of OV. Thus, identifying novel regulators as diagnostic and therapeutic targets for OV is still urgently required.
Minichromosome maintenance (MCM) proteins are essential for the initiation of DNA replication, and it has been detected to be overexpressed in various cancer tissues, including lung squamous cell carcinoma (Wu et al., 2018), breast cancer (Juríková et al., 2016), glioma (Cai et al., 2018), hepatocellular carcinoma (Liu et al., 2018), etc. As an important player in the initiation pathway of DNA replication, Minichromosome maintenance 10 (MCM10) was first identified in a yeast genetic screen and only presents in eukaryotes (Aves et al., 1998). At the same time, the involvement of MCM10 has also been found in DNA elongation, bolstering the activity of the CMG helicase on bypassing replication blocks (Langston et al., 2017; Lõoke et al., 2017) and promotion of replication fork progression and stability (Baxley and Bielinsky, 2017). In addition to the above group roles, MCM10, like other MCM family proteins, is abnormally expressed in various tumors and associated with prognosis. The overexpression of MCM10 is thought to promote the abnormal proliferation of prostate cancer (PC) cells and associated with poor prognosis of PC (Cui et al., 2018). Meanwhile, it is positively related to poor prognosis in breast cancer (Yang and Wang, 2019). In glioma, the knockdown of MCM10 in glioma cells resulted in decreased cell proliferation, migration and invasion (Kang et al., 2020). However, we have not found many articles on the relationship between MCM10 and OV. Based on the close relationship between MCM10 and various malignant tumors, we have reasons to believe that MCM10 is a potential prognostic marker.
Increasing research results proved that immunotherapy is a very promising therapeutic method in the treatment of malignant tumors, among which the blockade of immune checkpoints has displayed significant efficacy in various types of tumors (Topalian et al., 2016). Interfering with Cytotoxic T lymphocyte (CTLA-4) and ProgrammedDeath-1 (PD-1) reportedly has clinical benefits in several human cancers (Odunsi, 2017), so the characterizing associations between MCM10 and immune checkpoints will potentially enhance OV treatment.
In this study, the online tools TIMER and GEPIA were used to explore the expression of MCM10 in various malignancies. The prognostic value of MCM10 expression in OV was determined using the Kaplan-Meier Plotter and PrognoScan databases. LinkedOmics database was used to view genes and pathways associated with MCM10. The cBioPortal database was used to visualize and compare genetic alterations and explore the association between MCM10 and 25 immune checkpoints. Finally, qRT-PCR was used to detect the expression of MCM10 and analyze the relationship between its expression and prognosis. Our findings revealed an significant role for MCM10 in OV expression and prognosis, and also elucidated the relationship between MCM10 and multiple immune checkpoints.
Materials and Methods
Tissue Samples
A total of 22 cancerous and 50 paracancerous ovarian tissue samples were obtained during surgery. The study was approved by the jurisdictional clinical research ethics committees. All patients consented to the study.
qRT-PCR
Frozen tissues (100 mg) were ground into powder in liquid nitrogen, and then suspended in 1 ml TRIZOL Reagent (Invitrogen, United States of America), Total RNA was extracted using TRIZOL reagent. RNA was quantified using a spectrophotometer (Beckman, United States of America). RNA template and random primers were incubated at 70°C for 10 min to melt the secondary structure within the template, and cooled on ice for more than 2 min. Then the complete reaction mixture was incubated at 30°C for 10 min, 42°C for 60 min and 70°C for 15 min. PCR was performed in a total volume of 25 μl containing 1 μl of reverse-transcribed cDNA. After an initial incubation at 94°C for 5 min, the reaction mixtures were subjected to 35 cycles of amplification using the following protocols: 94°C for 45 s, 55°C for 45 s and 72°C for 45 s, followed by a final extension step at 72°C for 7 min. PCR products were analyzed by 1.2% agarose gel electrophoresis and stained with GoldView nucleic acid dye. Real-time RT-PCR was performed using ABI PRISM 7500 Sequence Detection System instrument and software (Applied Biosystems, United States). The relative expression level of MCM10 was measured using SYBR Green I dye-based method. The results were normalized to the expression of glyceraldehyde-3-phosphate dehydrogenase (GAPDH). The Ct values of the amplified products were used in conjunction with the 2−ΔΔCt method to analyze the data (Schmittgen and Livak, 2008; Zhang et al., 2010). The primers used were GAPDH: 5′-AGAAGGCTGGGGCTCATTTG-3′ (F), 5′-AGGGGCCATCCACAGTCTTC-3′ (R); MCM10:5′-CACAGAAATGAACAAGAA-3′(F),5′-AATAAGAACAAGGACACA-3′(R); Primers were synthesized by BGI Company.
GEPIA Database Analysis
GEPIA (http://gepia.cancer-pku.cn/index.html) is a recently developed bioinformatics platform that incorporates genotype tissue expression data from 9,736 tumors and 8,587 normal samples. In the current study, the “Expression DIY” component was used to analyze the EFNA1 expression levels in a variety of cancers and adjacent tissues, and p < 0.05 was used as the screening threshold significance level (Tang et al., 2017).
TIMER Database
TIMER (https://cistrome.shinyapps.io/timer/) is a tumor immunity database, including 10,897 cancer samples from the TCGA database, together with an abundance of tumor-infiltrating immune cells (TIICs) based on a deconvolution method from gene expression profiles. In the present study, the “Gene” module was applied to analyze the correlations between EFNA1 expression and immune cell infiltration. The immune cells analyzed included CD4+ T cells, CD8+ T cells, B cells, neutrophils, macrophages, and dendritic cells (DCs). The “Correlation” module of TIMER was used to analyze the associations between EFNA1 and other prognosis-related immune cell markers, including CD8+ T cells, all T cells collectively, B cells, monocytes, tumor-associated macrophages (TAMs), M1 and M2 macrophages, neutrophils, natural killer (NK) cells, and DCs (Li et al., 2017).
Survival Analysis and Prognostic Evaluation
Kaplan–Meier Plotter (http://kmplot.com/analysis/) and PrognoScan (http://dna00.bio.kyutech.ac.jp/PrognoScan/index.html) were used for prognostic analyses. Kaplan–Meier Plotter evaluated the prognostic significance of MCM10 mRNA expression in OV. Patients’ samples were divided into two groups based on the median MCM10 expression level, and the overall survival (OS) of patients with OV was analyzed. The examination probe ID was used for MCM10 was 220,651_s_at, 222,962_s_at and 223,570_at. The log-rank p-value and hazard ratio (HR) with 95% confidence intervals were calculated. The PrognoScan database mainly collects clinical prognostic information derived from 14 cancers from GEO (Gene Expression Omnibus) and various laboratories, and then applies a minimum p value approach in analyses. In the current study, it was used to analyze the prognostic value of MCM10 in OV and adjust the threshold to a Cox p value (Mizuno et al., 2009; Lánczky and Győrffy, 2021).
cBioPortal Database
cBioPortal for Cancer Genomics (http://cbioportal.org) is a large repository of genomics datasets. In the present study, cBioPortal was used to visualize and compare the changes in EFNA1 and immune checkpoints in OV. The correlations between MCM10 and immune checkpoints were also investigated. The immune checkpoints analyzed included PD-L1 (CD274), PD-L2 (PDCD1LG2), CD80, CD86, VTCN1, VSIR, HHLA2, TNFRSF14, PVR, CTL4, CD112 (NECTIN2), CD200, LGALS9, ICOSLG, TNFSF9, TNFSF4, CD70, TNFSF18, CD48, CTLA4, CD276, LILRB1, LILRB2, HAVCR2, CD47 and TNFRSF9(CD137) (Wu et al., 2019).
COSMIC Database Analysis of MCM10 Mutations in OV
The Catalogue of Somatic Mutations in Cancer (COSMIC) (https://cancer.sanger.ac.uk) is the most detailed and comprehensive resource to explore the effects of somatic mutations in human cancer. COSMIC database Contains 6 million coding mutations (Tate et al., 2019). Ovary in the“tissue distribution”and “mutation distribution”were chosen.
LinkedOmics Database Analysis of MCM10-Related Pathways
The LinkedOmics (http://www.linkedomics.org) database includes 32 cancer types from TCGA project and 11,158 patients with multiple omics and clinical data. It is also the first multi-omics database that integrates mass spectrometry–based global proteomics data generated by the Clinical Proteomics Cancer Analysis Alliance on selected TCGA tumor samples (Vasaikar et al., 2018). The differentially expressed genes related to MCM10 were screened from the TCGA OV cohort through the LinkFinder module in the database, and the correlations of the results were presented in volcano plots and heat maps, respectively, by Pearson correlation coefficient test. Gene Set Enrichment Analysis (GSEA) in the LinkInterpreter module performed functional module analysis of the Gene Ontology Biological Process (GO_BP), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. And 0.05 was considered as the p-value cutoff and the Spearman correlation test was conducted to analyze the results statistically.
Statistical Analysis
All data were expressed as means and standard deviations, and SPSS 19.0 software (SPSS Inc. Chicago, IL, United States) was the statistical analysis tool used. The student’s t-test and one-way analysis of variance were carried out to analyze the differences between groups. p < 0.05 was deemed to indicate the statistical significance. All experiments were conducted in triplicate as a minimum. The Kaplan–Meier Plotter and GEPIA results are presented as hazard ratios (HRs) and p values, and PrognoScan results are presented as Cox p values.
Results
mRNA Expression Levels of MCM10 in Different Types of Human Cancers
To evaluate the differences of MCM10 expression in tumor and normal tissues, the MCM10 mRNA levels in tumor and normal tissues of patients with multiple types of cancer were analyzed using the GEPIA and TIMER database. MCM10 expression was higher in BLCA, BRCA, CHOL, COAD, HNSC, KICH, KIRC, KIRP, LIHC, LUAD, PRAD, READ, STAD, THCA, and UCEC, compared with normal tissues in TIMER (Figure 1A). In GEPIA, high expression in BLCA, BRCA, CESC, COAD, DLBC, ESCA, GBM, HNSC, LUAD, LUSC, OV, READ, SKCM, STAD, TGCT, THYM, UCEC and UCS was observed. In addition, lower expression was observed in LAML (Figure 1B). MCM10 was highly expressed in OV (Figure 1C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/864578/fgene-13-864578-HTML/image_m/fgene-13-864578-g001.jpg)
FIGURE 1. MCM10 expression levels in different types of human cancers. (A) MCM10 expression in different tumor types in TIMER. (B) MCM10 expression in different tumor types in GEPIA. (C) Box plots comparing MCM10 expression in OV and unpaired normal tissues in GEPIA based on analysis of variance method (TCGA tumor versus TCGA normal + GTEx normal).*p < 0.05, **p < 0.01, ***p < 0.001.
Relationships Between MCM10 and the Prognosis of OV
The prognostic value of MCM10 expression in OV was evaluated using Kaplan-Meier plots and PrognoScan. The expression of MCM10 was significantly associated with the prognosis of OV patients. In analyses with the Kaplan-Meier plots, the OS 220651_s_at [HR = 1.14 (1 -1.3), p = 0.043], 222,962_s_at [HR = 1.59 (1.29 -1.97), p = 1.3e-05] and 223,570_at [HR = 1.43 91.17 -1.760, p = 0.00049] of OV patients with high MCM10 expression (Figures 2A–C) values were significantly lower than those of patients with low MCM10 expression. In PrognoScan database analysis of the prognostic potential of MCM10 in OV, high MCM10 expression in the GSE9891 cohort was associated with poor OS [HR = 1.36 (1.02 -1.810, p = 0.035626] (Figure 2D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/864578/fgene-13-864578-HTML/image_m/fgene-13-864578-g002.jpg)
FIGURE 2. Correlations between MCM10 expression and prognostic indicators in OV. (A–C) Correlations between MCM10 and OV prognoses in the Kaplan–Meier Plotter database from different expression callers (220,651_s_at,222,962_s_at, 223,570_at). (D) Survival curve from PrognoScan analysis for OS of patients with OV. HR = hazard ratio.
MCM10 Co‐Expression Network in OV
To understand the biological function of MCM10 in OV, the LinkFinder module in the LinkedOmics portal was used to examine the co-expression pattern of MCM10 in TCGA-OV. As shown in Figure 3A genes positively correlated with MCM10 were dark red dots, and 3,613 genes negatively correlated with MCM10 were dark green dots. Figures 3B–D represents the top 50 genes associated and negatively associated with the MCM10 signature, respectively. GO term annotation proved that the co‐expressed genes of MCM10 join mainly in chromosome segregation, spindle organization, DNA replication, cell cycle G2/M phase transition, mitotic cell cycle phase transition, cell cycle checkpoint, double-strand break repair, cytokinesis, negative regulation of mitotic cell cycle and protein localization to chromosome, etc (Figure 3C). KEGG pathway analysis indicated the enrichment in Cell cycle, DNA replication, Fanconi anemia pathway, Oocyte meiosis, Progesterone-mediated oocyte maturatio, Homologous recombination, Mismatch repair, Asthma, Graft-versus-host disease, and Staphylococcus aureus infection, etc (Figure 3C). It was found that GO terms and KEGG pathways were more concentrated in cell cycle and DNA replication.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/864578/fgene-13-864578-HTML/image_m/fgene-13-864578-g003.jpg)
FIGURE 3. The co‐expression genes with MCM10 from the LinkedOmics database in OV. (A) The whole significantly associated genes with MCM10 distinguished in OV cohort. (B–D) Top 50 genes positively and negatively related to MCM10 in OV (C,E) GO annotations and KEGG pathways of CLEC10A in LUAD cohort.
Mutation of MCM10 in OV
The cBioPortal was used to explore the mutation status of MCM families (MCM2, MCM3, MCM4, MCM5, MCM6, MCM7, MCM8, MCM9, MCM10 and MCMBP). Totally, 30.87% (96/311) of patients had genetic alterations (Figure 4A), of which MCM10 accounted for 6.43% (20/311), and the amplification was the most frequent mutation (Figure 4B). In COSMIC, we further assessed the mutation type of MCM10, Missense substitutions occurred in approximately 21.43% of the samples, synonymous substitutions occurred in 3.57% of the samples, and nonsense substitutions occurred in 3.57% of the samples (Figure 4C). The substitution mutations mainly occurred at A > G (25.00%) and G > T (25.00%), followed by A > T (12.50%), C > A (12.50%), C > G (12.50%), and G > A (12.50%) (Figure 4D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/864578/fgene-13-864578-HTML/image_m/fgene-13-864578-g004.jpg)
FIGURE 4. Mutation analysis of MCM10 in OV. (A) Mutation frequency of MCM10. (B) Mutation frequency of MCM10. (C–D) The mutation types of MCM10 in OV by Catalogue of Somatic Mutations in Cancer (COSMIC) database.
Correlation of MCM10 With Immune Infiltration and Immune Checkpoints
The presence of immune infiltration within tumors can generate important biomarkers to predict the prognosis of tumor patients, with impacts on radiotherapy, chemotherapy and therapy. Therefore, it is cardinal to study the relationship between MCM10 and immunity. We used the “Gene” module in TIMER for database search, entered the target gene MCM10, and selected OV. This module displays infiltration results, including TIMER, EPIC, MCP-COUNTER, CIBERSORT, CIBERSORT-ABS, QUANTISEQ, XCELL, QUANTISE and TIDE. The expression of MCM10 was positively correlated with CD4+T Cells (cor = 0.109, p = 0.0166) and Macrophages (cor = 0.101,p = 0.0275), and was positively correlated with B cells(cor = -0.1,p = 0.0288) and CD8+T Cells(cor = -0.139, p = 0.00223) (Figure 5A). The relationship between genetic changes in the MCM10 gene and 25 immune checkpoints was explored. Three datasets (MSK, TCGA and MSKCC), including 612 samples, were selected, and genomic studies revealed that MCM10 was involved in the alteration of OV immune checkpoints. The alterations of MCM10 and immune checkpoints in OV were visualized in a compact manner. In OV, the Genetic Alteration of MCM10 was 5% and mainly concentrated in Amplification. There was a part of Missense Mutation, second only to 6% of the immune checkpoint CD47, which was equal to HHLA2 and CD200. This indicated that MCM10 had a high mutation rate during the progression of OV (Figure 5B). Then, the association between MCM10 and each immune checkpoint was examined. Notably, the alterations in MCM10 showed statistically significant coexistence rather than rejection with the immune checkpoints CTLA4, TNFSF4, TNFSF18, CD80, ICOSLG, LILRB1 and CD200 (Table 1). These findings strongly suggested that MCM10 is a potential co-regulator of the OV immune checkpoint.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/864578/fgene-13-864578-HTML/image_m/fgene-13-864578-g005.jpg)
FIGURE 5. Correlation of MCM10 with immune infiltration and immune checkpoints. (A) the correlations between MCM10 and immune cell infiltrations from TIMER. (B) Landscape of MCM10 and immune checkpoint alteration in OV from cBioPortal.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/864578/fgene-13-864578-HTML/image_m/fgene-13-864578-t001.jpg)
TABLE 1. Mutual-exclusivity analysis between MCM10 and multiple-immune checkpoints in ovarian cancer.
Validation of MCM10 in OV Tissues
To explore the expression of MCM10 in OV, qRT-PCR validation was chosen. The results suggested that the expression of MCM10 was elevated in OV (Figure 6A). According to the median value of expression, the patients were divided into high and low expression groups to judge the prognosis of patients (Figure 6B). The results illustrated that OV patients with high MCM10 had a worse prognosis.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/864578/fgene-13-864578-HTML/image_m/fgene-13-864578-g006.jpg)
FIGURE 6. Expression of MCM10 in Independent OV Cohorts. (A) qRT-PCR showed that MCM10 expression was up-regulated in OV tissues. (B) High MCM10 expression have worse prognosis in OV patients.
Disscussion
One of the characteristics of tumor cells is unlimited proliferation. Therefore, many proteins related to DNA replication have been considered as potential cancer biomarkers, including MCM protein (Yu et al., 2020). MCM protein has been considered as a biomarker of dysplasia and tumor (Wang et al., 2020). There are several hypotheses about the mechanism of MCM dysregulation leading to tumorigenesis. The first is genomic instability (GIN), because the formation of cancer cells is caused by the accumulation of mutations in oncogenes and tumor suppressor genes. In many studies, GIN caused by MCM mutation has been proved to be related to the occurrence of malignant tumors. (Chuang et al., 2010). Secondly, damage to MCM induces replicative stress, a critical step in the initiation of the oncogenic process (Gaillard et al., 2015). Finally, the study found that MCM family proteins participate in the progression of cell cycle pathways. For example, knockdown of MCM2 reduces the expression of cyclinD1, cyclinA and CDK4, knockdown of MCM3 reduces the expression of cyclinA, knockdown of MCM6 causes CyclinA, CyclinB1,CyclinD1, silencing of MCM7 reduces cyclinD1, cyclinE2 and CDK2, and down-regulation of cyclinD1 in breast cancer cells with MCM10 knockdown (Zhang et al., 2015; Qiu et al., 2017; Liu et al., 2018; Yang and Wang, 2019). These results all suggest that mutations in MCM proteins cause changes in various cyclins, and ultimately aberrant cell cycle progression leads to tumorigenesis. There are few studies on MCM10 in the MCM protein family, especially its relationship with OV, so this study focused on this gene, which has not received much attention but is extremely significant. MCM10 acts as an vital scaffold for DNA replication and protection against replication stress under normal conditions. However, under pathological conditions, MCM10 is frequently deregulated, and gene amplification and overexpression are very common in cancer (Baxley and Bielinsky, 2017). Our results also demonstrated that the changes of the MCM10 gene in OV are mainly concentrated in the amplification, and it was discovered that most of the mutations in MCM10 are missense mutations (93%), and the rest are roughly divided into splicing mutations (3.7%) and nonsense mutations (3.2%) (Gao et al., 2013; Kang et al., 2013), which is consistent with our findings. The mechanism by which MCM10 causes OV, in addition to the aforementioned, may also be related to the specific relationship between MCM10 and female ovaries, where MCM10 is highly expressed in adult female ovaries (Graveley et al., 2011). Based on this, Reubens et al. conducted further research and believed that MCM10 plays a unique biological role in the development or maintenance of the female germline (Reubens et al., 2015). So the mutation of MCM10 may be another cause of OV progression. These data clearly indicated that MCM10 changes in the cancer genome, but whether these changes are the causes or the results of OV still needs further studies to confirm.
Our results proved that the expression of MCM10 in various malignant tumors is different from that in normal tissues, except for the low expression of LAML in the GEPIA database, because the data in the TIMER database are all from TCGA, and normal control samples of some malignant tumors are insufficient, but MCM10 is highly expressed in malignant tumors with differential expression. Since we only obtained high expression of MCM10 in OV from GEPIA, we further verified the expression of MCM10 in OV by qRT-PCR, and it was found that the expression of MCM10 in OV was higher than that in normal ovarian tissue. Subsequently, the relationship between the expression of MCM10 and the prognosis of OV patients was examined in four databases, and it was obvious that the prognosis of OV patients with high expression of MCM10 was worse. After being divided into high and low expression, the patients in the high expression group also had a poor prognosis. These findings all demonstrated that MCM10, as a promising prognostic biomarker, is increased in OV.
Having verified the differential expression and prognostic potential of MCM10 in OV, in order to better serve the clinic, we next explored the possible pathogenic mechanism of MCM10. In order to test our previous speculation about the pathogenic mechanism of MCM10, the effects of MCM10 changes in OV on the transcriptome were explored, and it was found that 4,297 positively correlated genes and 3,613 negatively correlated genes were changed accordingly, which suggested that the alterations in MCM10 have broad impacts on the transcriptome. The analysis results of GO and KEGG both illustrated that the pathway enrichment of other gene changes caused by MCM10 changes mainly concentrated in the related pathways of cell cycle and DNA replication. Since MCM10 itself participates in DNA replication, it is not surprising that the changes are mainly concentrated in DNA replication. However, the changes of MCM10 are closely related to the cell cycle, which may be closely related to the cell cycle, because the DNA replication process depends on the regulation of the cell cycle (Tachibana et al., 2005). Previous studies suggested that MCM10 may be part of a high-priority group of genes that may promote cell cycle-related processes in cancer cells (Cerami et al., 2012). A series of cell cycle checkpoint inhibitors targeting OV are already under development or clinical trials (Pujade-Lauraine, 2017). For example, cells initiate multiple responses to protect the genome and ensure survival against DNA damage, and unsuccessful DNA damage repair can lead to mitotic abnormalities and cell death (Lin et al., 2017). High-grade serous ovarian cancer (HGSOC) relies heavily on G2 checkpoint blockade to promote DNA damage repair, which is a process that opens up a new perspective for the treatment of OV (Haynes et al., 2018). If blocking MCM10 can affect both DNA replication and cell cycle, why not? Therefore, the combined treatment of DNA replication and cell cycle intervention may benefit OV patients.
People who care about OV treatment know that OV treatment cannot be cured by single-agent therapy because the results of single-agent studies on OV so far have been disappointing (Disis et al., 2019). The results of the drug combination may have some benefit, but the results are not better than the historical control (Lee et al., 2019). Therefore, it is necessary to explore a new combination therapy method. Tumor immunotherapy is considered as a new and potential tumor treatment method. The infiltration state of tumor-associated immune cells in vivo together constitutes the immune microenvironment of tumor cells, and these immune cells may have tumor antagonism or tumor promotion (Jain, 2021). Our results demonstrated that the expression level of MCM10 is negatively correlated with B cells and CD8+T Cells, and positively correlated with CD4+T Cells and Macrophages. This suggests that MCM10 plays a particular role in the immune infiltration of OV. In addition to the critical role of the immune microenvironment in anticancer immunity, another most popular approach in immunotherapy is immune checkpoint blockade (Huang et al., 2020). PD-1/PD-L1 and CTLA-4 are considered as the most principal immune checkpoints at present. In the past decade, great progress has been made in the field of immune checkpoint-related researches. immune checkpoint blockade (ICB) has been very successful in this type of cancer (Havel et al., 2019). However, since the current FAD-approved immune checkpoint inhibitors (ICI) are all monoclonal antibodies (mAbs), there are many shortcomings. Therefore, the therapeutic effects of ICB on OV are still limited. The study on small molecule inhibitors to eliminate the limitations of mAbs is a new direction for ICB therapy (Zhang et al., 2020; Zhang et al., 2021). More and more evidence showed that small molecule inhibitors that target oncogenic play a role far beyond the biological behavior of tumors. Some studies have found that some small molecule inhibitors directly participate in mediating the tumor microenvironment and promoting tumor cell death (Chen et al., 2019; Ziogas et al., 2021). For example, the inhibitors of CDK4/6 can synergize with PD-1 blockade and benefit the treatment of OV (Zhang et al., 2020). It has been revealed that the matrix metalloproteinases (MMPs) inhibitor, SB-3CT, can enhance the effects of PD-1 and CTLA-4 blockade in primary and metastatic tumors in studies (Ye et al., 2020). Small-molecule inhibitor JQ1 targeting BET bromodomains reduces PD-L1 expression, while attenuating progression in PC models (Mao et al., 2019). All of the above studies have proved that ICI combined with small molecule inhibitors is an effective way to address the shortcomings of current ICIs such as low oral availability, long tissue retention time and poor membrane permeability. Changes in seven ICIs occurred simultaneously, and considering that MCM10 is a member of a high-priority gene, the combination of its small-molecule inhibitor and ICI greatly benefits OV patients.
In our study, the expression and prognosis of MCM10 in OV were analyzed by bioinformatics methods, and the related pathways of MCM10 that were related to immune infiltration and immune checkpoints were analyzed. However, inevitably, our tests have certain limitations. Firstly, the number of patients enrolled in our validation experiment is relatively small, and we will update the number of patients in the later stage. Secondly, we only propose possible pathogenic pathways, and further experimental verification is needed, for example, we could test the effects of an artificial knockdown of MCM10 expression on tumor size and progression in cell models, validate qPCR results from tumor samples with Western blots, test for differences between MCM10 expression in CD8 mutant backgrounds and check if MCM10 expression differences cause differences in replication timing, genome stability or cell cycle defects. Finally, we need animal experiments and long-term clinical experiments for MCM10 small-molecule inhibitors, and then apply them to patients and use them in combination with ICI. Most importantly, we will solve these problems slowly, and our efforts have also achieved certain results. Therefore, future researches on Mcm10’s relationship to cancer development and progression may lead to discoveries with momentous prognostic and even therapeutic value.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics Statement
The studies involving human participants were reviewed and approved by the Ethics Committee of Gansu Provincial Maternity and Child-Care Hospital. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
JZ and ZW designed the study. YW and JL analyzed the data and wrote the manuscript. HW and XT revised manuscript. All authors have read and approved the final manuscript.
Funding
This study was supported by the Gansu Province Health Industry Research Plan (grant number GSWSKY 2020-16).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Aves, S. J., Tongue, N., Foster, A. J., and Hart, E. A. (1998). The Essential Schizosaccharomyces pombe Cdc23 DNA Replication Gene Shares Structural and Functional Homology with the Saccharomyces cerevisiae DNA43 ( MCM10 ) Gene. Curr. Genet. 34 (3), 164–171. doi:10.1007/s002940050382
Baxley, R., and Bielinsky, A.-K. (2017). Mcm10: A Dynamic Scaffold at Eukaryotic Replication Forks. Genes. 8 (2), 73. doi:10.3390/genes8020073
Cai, H.-Q., Cheng, Z.-J., Zhang, H.-P., Wang, P.-F., Zhang, Y., Hao, J.-J., et al. (2018). Overexpression of MCM6 Predicts Poor Survival in Patients with Glioma. Hum. Pathol. 78, 182–187. doi:10.1016/j.humpath.2018.04.024
Cerami, E., Gao, J., Dogrusoz, U., Gross, B. E., Sumer, S. O., Aksoy, B. A., et al. (2012). The cBio Cancer Genomics Portal: An Open Platform for Exploring Multidimensional Cancer Genomics Data. Cancer Discov. 2 (5), 401–404. doi:10.1158/2159-8290.cd-12-0095
Chen, S., Song, Z., and Zhang, A. (2019). Small-Molecule Immuno-Oncology Therapy: Advances, Challenges and New Directions. Ctmc 19 (3), 180–185. doi:10.2174/1568026619666190308131805
Chuang, C.-H., Wallace, M. D., Abratte, C., Southard, T., and Schimenti, J. C. (2010). Incremental Genetic Perturbations to MCM2-7 Expression and Subcellular Distribution Reveal Exquisite Sensitivity of Mice to DNA Replication Stress. PLoS Genet. 6 (9), e1001110. doi:10.1371/journal.pgen.1001110
Cui, F., Hu, J., Ning, S., Tan, J., and Tang, H. (2018). Overexpression of MCM10 Promotes Cell Proliferation and Predicts Poor Prognosis in Prostate Cancer. Prostate 78 (16), 1299–1310. doi:10.1002/pros.23703
Disis, M. L., Taylor, M. H., Kelly, K., Beck, J. T., Gordon, M., Moore, K. M., et al. (2019). Efficacy and Safety of Avelumab for Patients with Recurrent or Refractory Ovarian Cancer: Phase 1b Results From The JAVELIN Solid Tumor Trial. JAMA Oncol. 5 (3), 393–401. doi:10.1001/jamaoncol.2018.6258
Gaillard, H., García-Muse, T., and Aguilera, A. (2015). Replication Stress and Cancer. Nat. Rev. Cancer 15 (5), 276–289. doi:10.1038/nrc3916
Gao, J., Aksoy, B. A., Dogrusoz, U., Dresdner, G., Gross, B., Sumer, S. O., et al. (2013). Integrative Analysis of Complex Cancer Genomics and Clinical Profiles Using the cBioPortal. Sci. Signal. 6 (269), pl1. doi:10.1126/scisignal.2004088
Graveley, B. R., Brooks, A. N., Carlson, J. W., Duff, M. O., Landolin, J. M., Yang, L., et al. (2011). The Developmental Transcriptome of Drosophila melanogaster. Nature 471 (7339), 473–479. doi:10.1038/nature09715
Havel, J. J., Chowell, D., and Chan, T. A. (2019). The Evolving Landscape of Biomarkers for Checkpoint Inhibitor Immunotherapy. Nat. Rev. Cancer 19 (3), 133–150. doi:10.1038/s41568-019-0116-x
Haynes, B., Murai, J., and Lee, J.-M. (2018). Restored Replication Fork Stabilization, a Mechanism of PARP Inhibitor Resistance, Can Be Overcome by Cell Cycle Checkpoint Inhibition. Cancer Treat. Rev. 71, 1–7. doi:10.1016/j.ctrv.2018.09.003
Huang, X., Tang, T., Zhang, G., Hong, Z., Xu, J., Yadav, D. K., et al. (2020). Genomic Investigation of Co-targeting Tumor Immune Microenvironment and Immune Checkpoints in Pan-Cancer Immunotherapy. npj Precis. Onc. 4 (1), 29. doi:10.1038/s41698-020-00136-1
Irusta, G. (2021). Roads to the Strategic Targeting of Ovarian Cancer Treatment. Reprod. Camb. Engl. 161 (1), R1–R11. doi:10.1530/REP-19-0593
Jain, K. K. (2021). Personalized Immuno-Oncology. Med Princ Pract Int. J. the Kuwait Univ. Health Sci. Centre 30 (1), 1–16. doi:10.1159/000511107
Juríková, M., Danihel, Ľ., Polák, Š., and Varga, I. (2016). Ki67, PCNA, and MCM Proteins: Markers of Proliferation in the Diagnosis of Breast Cancer. Acta Histochem. 118 (5), 544–552. doi:10.1016/j.acthis.2016.05.002
Kang, G., Hwang, W. C., Do, I.-G., Wang, K., Kang, S. Y., Lee, J., et al. (2013). Exome Sequencing Identifies Early Gastric Carcinoma as an Early Stage of Advanced Gastric Cancer. PloS one 8 (12), e82770. doi:10.1371/journal.pone.0082770
Kang, P., Han, Z., Liao, Z., Zhang, H., Jia, W., and Tian, Y. (2020). Knockdown of MCM10 Gene Impairs Glioblastoma Cell Proliferation, Migration and Invasion and the Implications for the Regulation of Tumorigenesis. J. Mol. Neurosci. 70 (5), 759–768. doi:10.1007/s12031-020-01486-y
Lánczky, A., and Győrffy, B. (2021). Web-Based Survival Analysis Tool Tailored for Medical Research (KMplot): Development and Implementation. J. Med. Internet Res. 23 (7), e27633. doi:10.2196/27633
Langston, L. D., Mayle, R., Schauer, G. D., Yurieva, O., Zhang, D., Yao, N. Y., et al. (2017). Mcm10 Promotes Rapid Isomerization of CMG-DNA for Replisome Bypass of Lagging Strand DNA Blocks. eLife 6, e29118. doi:10.7554/eLife.29118
Lee, J. M., Minasian, L., and Kohn, E. C. (2019). New Strategies in Ovarian Cancer Treatment. Cancer 125 (Suppl. 24), 4623–4629. doi:10.1002/cncr.32544
Lheureux, S., Braunstein, M., and Oza, A. M. (2019). Epithelial Ovarian Cancer: Evolution of Management in the Era of Precision Medicine. CA Cancer J. Clin. 69 (4), 280–304. doi:10.3322/caac.21559
Li, T., Fan, J., Wang, B., Traugh, N., Chen, Q., Liu, J. S., et al. (2017). TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res. 77 (21), e108–e110. doi:10.1158/0008-5472.can-17-0307
Lin, A. B., McNeely, S. C., and Beckmann, R. P. (2017). Achieving Precision Death with Cell-Cycle Inhibitors that Target DNA Replication and Repair. Clin. Cancer Res. 23 (13), 3232–3240. doi:10.1158/1078-0432.ccr-16-0083
Liu, Z., Li, J., Chen, J., Shan, Q., Dai, H., Xie, H., et al. (2018). MCM Family in HCC: MCM6 Indicates Adverse Tumor Features and Poor Outcomes and Promotes S/G2 Cell Cycle Progression. BMC cancer 18 (1), 200. doi:10.1186/s12885-018-4056-8
Lõoke, M., Maloney, M. F., and Bell, S. P. (2017). Mcm10 Regulates DNA Replication Elongation by Stimulating the CMG Replicative Helicase. Genes. Dev. 31 (3), 291–305. doi:10.1101/gad.291336.116
Mao, W., Ghasemzadeh, A., Freeman, Z. T., Obradovic, A., Chaimowitz, M. G., Nirschl, T. R., et al. (2019). Immunogenicity of Prostate Cancer Is Augmented by BET Bromodomain Inhibition. J. Immunother. cancer 7 (1), 277. doi:10.1186/s40425-019-0758-y
Mizuno, H., Kitada, K., Nakai, K., and Sarai, A. (2009). PrognoScan: a New Database for Meta-Analysis of the Prognostic Value of Genes. BMC Med. Genomics 2, 18. doi:10.1186/1755-8794-2-18
Odunsi, K. (2017). Immunotherapy in Ovarian Cancer. Ann. Oncol. 28 (Suppl. l_8), viii1–viii7. doi:10.1093/annonc/mdx444
Pujade-Lauraine, E. (2017). New Treatments in Ovarian Cancer. Ann. Oncol. 28 (Suppl. l_8), viii57–viii60. doi:10.1093/annonc/mdx442
Qiu, Y.-T., Wang, W.-J., Zhang, B., Mei, L.-L., and Shi, Z.-Z. (2017). MCM7 Amplification and Overexpression Promote Cell Proliferation, Colony Formation and Migration in Esophageal Squamous Cell Carcinoma by Activating the AKT1/mTOR Signaling Pathway. Oncol. Rep. 37 (6), 3590–3596. doi:10.3892/or.2017.5614
Reubens, M. C., Biller, M. D., Bedsole, S. E., Hopkins, L. T., Ables, E. T., and Christensen, T. W. (2015). Mcm10 Is Required for Oogenesis and Early Embryogenesis in Drosophila. Mech. Dev. 138, 291–299. doi:10.1016/j.mod.2015.09.002
Schmittgen, T. D., and Livak, K. J. (2008). Analyzing Real-Time PCR Data by the Comparative C(T) Method. Nat. Protoc. 3 (6), 1101–1108. doi:10.1038/nprot.2008.73
Tachibana, K.-e. K., Gonzalez, M. A., and Coleman, N. (2005). Cell-cycle-dependent Regulation of DNA Replication and its Relevance to Cancer Pathology. J. Pathol. 205 (2), 123–129. doi:10.1002/path.1708
Tang, Z., Li, C., Kang, B., Gao, G., Li, C., and Zhang, Z. (2017). GEPIA: a Web Server for Cancer and Normal Gene Expression Profiling and Interactive Analyses. Nucleic acids Res. 45 (W1), W98–W102. doi:10.1093/nar/gkx247
Tate, J. G., Bamford, S., Jubb, H. C., Sondka, Z., Beare, D. M., Bindal, N., et al. (2019). COSMIC: the Catalogue of Somatic Mutations in Cancer. Nucleic acids Res. 47 (D1), D941–D947. doi:10.1093/nar/gky1015
Topalian, S. L., Taube, J. M., Anders, R. A., and Pardoll, D. M. (2016). Mechanism-driven Biomarkers to Guide Immune Checkpoint Blockade in Cancer Therapy. Nat. Rev. Cancer 16 (5), 275–287. doi:10.1038/nrc.2016.36
Vasaikar, S. V., Straub, P., Wang, J., and Zhang, B. (2018). LinkedOmics: Analyzing Multi-Omics Data within and across 32 Cancer Types. Nucleic acids Res. 46 (D1), D956–D963. doi:10.1093/nar/gkx1090
Wang, Y., Chen, H., Zhang, J., Cheng, A. S. L., Yu, J., To, K. F., et al. (2020). MCM Family in Gastrointestinal Cancer and Other Malignancies: From Functional Characterization to Clinical Implication. Biochimica Biophysica Acta (BBA) - Rev. Cancer 1874 (2), 188415. doi:10.1016/j.bbcan.2020.188415
Wu, P., Heins, Z. J., Muller, J. T., Katsnelson, L., de Bruijn, I., Abeshouse, A. A., et al. (2019). Integration and Analysis of CPTAC Proteomics Data in the Context of Cancer Genomics in the cBioPortal. Mol. & Cell. Proteomics 18 (9), 1893–1898. doi:10.1074/mcp.tir119.001673
Wu, W., Wang, X., Shan, C., Li, Y., and Li, F. (2018). Minichromosome Maintenance Protein 2 Correlates with the Malignant Status and Regulates Proliferation and Cell Cycle in Lung Squamous Cell Carcinoma. Ott 11, 5025–5034. doi:10.2147/ott.s169002
Yang, W. D., and Wang, L. (2019). MCM10 Facilitates the Invaded/migrated Potentials of Breast Cancer Cells via Wnt/β-Catenin Signaling and Is Positively Interlinked with Poor Prognosis in Breast Carcinoma. J. Biochem. Mol. Toxicol. 33 (7), e22330. doi:10.1002/jbt.22330
Ye, Y., Kuang, X., Xie, Z., Liang, L., Zhang, Z., Zhang, Y., et al. (2020). Small-molecule MMP2/MMP9 Inhibitor SB-3CT Modulates Tumor Immune Surveillance by Regulating PD-L1. Genome Med. 12 (1), 83. doi:10.1186/s13073-020-00780-z
Yu, S., Wang, G., Shi, Y., Xu, H., Zheng, Y., and Chen, Y. (2020). MCMs in Cancer: Prognostic Potential and Mechanisms. Anal. Cell. Pathol. Amst. 2020, 3750294. doi:10.1155/2020/3750294
Zhang, F., Wang, Z.-m., Liu, H.-y., Bai, Y., Wei, S., Li, Y., et al. (2010). Application of RT-PCR in Formalin-Fixed and Paraffin-Embedded Lung Cancer Tissues. Acta Pharmacol. Sin. 31 (1), 111–117. doi:10.1038/aps.2009.178
Zhang, J., Zhang, Y., Qu, B., Yang, H., Hu, S., and Dong, X. (2021). If Small Molecules Immunotherapy Comes, Can the Prime Be Far behind? Eur. J. Med. Chem. 218, 113356. doi:10.1016/j.ejmech.2021.113356
Zhang, Q.-F., Li, J., Jiang, K., Wang, R., Ge, J.-l., Yang, H., et al. (2020). CDK4/6 Inhibition Promotes Immune Infiltration in Ovarian Cancer and Synergizes with PD-1 Blockade in a B Cell-dependent Manner. Theranostics 10 (23), 10619–10633. doi:10.7150/thno.44871
Zhang, X., Teng, Y., Yang, F., Wang, M., Hong, X., Ye, L.-G., et al. (2015). MCM2 Is a Therapeutic Target of Lovastatin in Human Non-small Cell Lung Carcinomas. Oncol. Rep. 33 (5), 2599–2605. doi:10.3892/or.2015.3822
Keywords: ovarian cancer, Mcm10, bioinformatics, prognosis biomarkers, immune checkpoints
Citation: Wu Z, Wang Y, Li J, Wang H, Tuo X and Zheng J (2022) MCM10 is a Prognostic Biomarker and Correlated With Immune Checkpoints in Ovarian Cancer. Front. Genet. 13:864578. doi: 10.3389/fgene.2022.864578
Received: 28 January 2022; Accepted: 02 May 2022;
Published: 19 May 2022.
Edited by:
Alip Ghosh, University of Maryland, United StatesReviewed by:
Sriram Vijayraghavan, MUSC Health, United StatesDebanjali Dasgupta, Mayo Clinic, United States
Copyright © 2022 Wu, Wang, Li, Wang, Tuo and Zheng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jing Zheng, 394146806@qq.com
†These authors have contributed equally to this work