Skip to main content

KAZN as a diagnostic marker in ovarian cancer: a comprehensive analysis based on microarray, mRNA-sequencing, and methylation data

Abstract

Background

Ovarian cancer (OC) is among the deadliest malignancies in women and the lack of appropriate markers for early diagnosis leads to poor prognosis in most cases. Previous studies have shown that KAZN is involved in multiple biological processes during development, such as cell proliferation, differentiation, and apoptosis, so defects or aberrant expression of KAZN might cause queer cell behaviors such as malignancy. Here we evaluated the KAZN expression and methylation levels for possible use as an early diagnosis marker for OC.

Methods

We used data from Gene Expression Omnibus (GEO) microarrays, The Cancer Genome Atlas (TCGA), and Clinical Proteomic Tumor Analysis Consortium (CPTAC) to investigate the correlations between KAZN expression and clinical characteristics of OC by comparing methylation levels of normal and OC samples. The relationships among differentially methylated sites in the KAZN gene, corresponding KAZN mRNA expression levels and prognosis were analyzed.

Results

KAZN was up-regulated in ovarian epithelial tumors and the expression of KAZN was correlated with the patients’ survival time. KAZN CpG site cg17657618 was positively correlated with the expression of mRNA and the methylation levels were significantly differential between the group of stage “I and II” and the group of stage “III and IV”. This study also presents a new method to classify tumor and normal tissue in OC using DNA methylation pattern in the KAZN gene body region.

Conclusions

KAZN was involved in ovarian cancer pathogenesis. Our results demonstrate a new direction for ovarian cancer research and provide a potential diagnostic biomarker as well as a novel therapeutic target for clinical application.

Peer Review reports

Background

Ovarian cancer is one of the most lethal gynecological malignancies in women. Due to the absence of symptoms at early stages and late detection, ovarian cancer is usually fatal. Data from SEER 182010–2016 show that 5-years relative survival for all patients with ovarian cancer is 46.8%, but 5-years relative survival in patients at early-stages is significantly better than in patients at advanced-stages (92.6% vs30.2%) [1]. The most common histological types are epithelial tumors, including high-grade serous carcinoma (HGSOC, accounting for ~ 60%), endometrioid, clear cell, mucinous, and low-grade serous carcinomas [2]. Recent studies report that the different histological subtypes of ovarian cancer may have distinct origins. Defects in BRCA1 and BRCA2 are well known genetic risks [3, 4], and mutations in many other genes are also associated with ovarian cancer, such as MSH2, MSH6, and TP53 [5,6,7,8,9], but in many cases the involved genetic factors remain unknown. Because of the genetic diversity and histological heterogeneity of ovarian cancer, the application of these studies for clinical diagnosis and effective therapy is rare and further investigation is urgently required [10].

Epigenetic changes are integral to all aspects of cancer genomics, such as DNA methylation, which is among the common epigenetic mechanisms involved in the formation and development of cancer [11, 12]. A large body of studies has shown that aberrant methylation of global DNA or specific genes may affect the progression [13] and prognosis of ovarian cancer [14], but the mechanisms for the aberrant DNA methylations to be involved in ovarian cancer remain largely unknown.

KAZN is an evolutionarily conserved gene initially identified in keratinocytes and then found to be widely expressed across different human tissues [15]. A study with mouse eggs demonstrates that the cellular localization of KAZN changes dynamically during the development [16]. KAZN is partially co-localized with desmoplakin and periplakin at desmosome and involved in the interplay between adherens junctions and desmosomes [15]. By binding to actin and intermediate filament, KAZN can affect the cell shape and remodeling of cytoskeletal networks [17]. In the nucleus, KAZN was found to be associated with the cell cycle, gene regulation, and matrix stability. KAZN interacts with apoptotic regulators, ARC and Bax, and plays an important role in apoptosis and cell growth [18]. KAZN is up-regulated during keratinocyte terminal differentiation [17] and is dispensable for murine epidermal morphogenesis and homeostasis [19]. Since KAZN is involved virtually in all aspects of cell development, its roles in carcinogenesis can be speculated [20, 21]. Therefore, according to the dynamic changes of KAZN expression during cell development and the regulation of gene expression by epigenetics in time and space [22, 23], we postulate that the expression and methylation of KAZN are crucial in the occurrence and development of ovarian cancer.

In this study, we conducted a comprehensive analysis on KAZN, using the data from GEO, TCGA, GTEx, and CPTAC datasets. We compared KAZN expression at levels of mRNA, protein, and DNA methylation, and detected the correlation between KAZN expression and survival time of the patients. Our results indicated that the expression and methylation patterns of KAZN were closely associated with the oncogenesis of ovarian cancer. These results can contribute to understanding the molecular mechanisms of tumor occurrence and progression and can be used to develop new diagnostic as well as treatment strategies.

Materials and methods

Gene expression omnibus (GEO) database

The GEO database (https://www.ncbi.nlm.nih.gov/geo/) is a public storage repository of microarray, methylation, and next-generation sequencing data. We searched the GEO datasets with the keywords “ovarian cancer” and “mRNA” and then filtered the documents by the following criteria: 1, organisms from the Homo species; and 2, datasets that contain both the tumor and non-tumor tissues (normal, benign, or para-carcinomas tissue). R package ‘GEOquery’ was used for data download and preprocessing. R package ‘Limma’ was used to normalize and analyze the data. A t-test is a statistical test used to compare the means of two groups. R package ‘ggpubr’ was used to plot the boxplot.

Comprehensive meta-analysis

R package ‘meta’ was used for performing a comprehensive meta-analysis of GEO data. The analysis of KAZN expression in the normal and tumor group was displayed on forest plots that illustrate the standardized mean difference (SMD) and the 95% confidential interval (CI). The chi-squared test of Q and the I2 statistic were calculated to assess heterogeneity across the studies and to determine the appropriateness of applying either a random-effects model or fixed-effects model to the pooling process. Influence analyses were conducted to investigate the relative influence of each individual study on the pooled effect size using R package ‘meta’ metainf function. To measure publication bias, Egger’s and Begg’s tests and a funnel plot, for which significance was p < 0.05, were performed.

TCGA, GTEx, and GEO gene expression data integration and differentially expressed genes (DEGs) identification

RNA-Seq data of 356 ovarian cancer and 180 normal tissue samples were downloaded from Genotype-Tissue Expression (GTEx) (https://www.gtexportal.org/home/) and TCGA via NCI’s Genomic Data Commons (GDC) portal (https://portal.gdc.cancer.gov/). Besides, five gene expression profiles - Datasets GSE137238, GSE101108, GSE101948, GSE143897, and GSE132107 - were downloaded from GEO, with a sample size of 323 ovarian cancer and 8 normal tissue samples. We integrated the data and ran ‘ComBat-seq’ in the R package ‘SVAseq’ (Version 3.38; http://bioconductor.org/packages/release/bioc/html/sva.html) to correct the batch effect. ‘DESeq2’ (Version 1.30. 1[24]; http://bioconductor.org/packages/release/bioc/html/DESeq2.html) was used to detect the differential gene expression between ovarian cancer and normal ovarian tissue samples.

Proteomic analysis

Proteomic data were downloaded from The National Cancer Institute’s Clinical Proteomic Tumor Analysis Consortium (CPTAC) (https://proteomics.cancer.gov/programs/cptac) [25]. After quality evaluation, a total of 84 ovarian cancer tissue and 22 normal tissue samples from the CPTAC Ovarian Cancer Confirmatory Study were retained for further proteomic study. Differentially expressed proteins (DEPs) between normal and tumor tissues were identified by R package ‘Limma’ with a FDR < 0.05 cut-off criteria.

Survival analysis

Overall survival was computed as the number of years between the year of diagnosis and the year of death from all causes, the date of the last follow-up, or 5-year censored survival data. Kaplan–Meier curves comparing overall survival according to subgroups divided by expression of KAZN and a log-rank test were used to assess the survival distributions across the subgroups. Samples were divided into either “high” or “low” groups with the cutoff at the lower quartile and upper quartile of KAZN gene expression. In log-rank, p < 0.05 was considered statistically significant. The survival curve was plotted by R package ‘survival’ and ‘survminer’.

Cox regressions analysis and ROC curve

The relationship between KAZN mRNA and patients’ overall survival was analyzed by univariate Cox regression using R package ‘survival’ and ‘survminer’. The forest plot for the KAZN hazard ratios and confidence intervals was created by ggplot2. To assess the performance of the gene risk model and compare the prognostic value, the time-dependent receiver operating characteristic (ROC) curve for this model was plotted using R package ‘survivalROC’.

Methylation level and mRNA expression level correlation analysis

The DNA methylation data (Illumina Human Methylation 27 k) of ovarian cancer (n = 601) and normal tissues (n = 12) and corresponding clinical information were obtained from TCGA (https://portal.gdc.cancer.gov/). The obtained DNA methylation data were further analyzed using Perl script and R package ‘ggpubr’. The Spearman rank correlation coefficient test was used to examine statistical significance in differences between DNA methylation and the expression of KAZN. R package ‘ggplot2’ was used to plot.

Tissue preparation

For the quantitative real-time PCR (qPCR) analyses, 6 ovarian cancer tissues (HGSOC) and 8 normal ovarian tissues were obtained from ovarian cancer patients in operation from The Third Affiliated Hospital, Harbin Medical University (Harbin, China).

RNA extraction and qRT-PCR analyses

RNA isolation of ovarian tissue samples were conducted through TRizol reagent (Invitrogen) according to the manufacturer’s instructions. Total RNAs were reversely transcribed into cDNAs and then used to perform qRT-PCR with Biosystems (ABI) 7500 platform. Homo sapiens β-actin (beta ACTB) was selected as the internal reference gene. The primer sequences were as follows: KAZN forward 5′-GGCAGATGAAGGAGATGTTGGCGAAGG-3′; KAZN reverse 5′-CTCTCCTTGCGGTGCTGCTCATAGTTG-3′; β-actin forward 5′-GGGAAATCGTGCGTGACATT-3′; β-actin reverse 5′-GGAACCGCTCATTGCCAAT-3′. The KAZN gene expression was determined by the subtracting their threshold cycle values (CT) to CT of β-actin gene.

KAZN gene methylation pattern analysis

GEO datasets Illumina Infinium 450 K BeadChips were used to detect the methylation sites. R package ‘minfi’ was used to check the different CpG sites. Online wANNOVAR webserver was used to annotate the CpG sites [26]. R package ‘pheatmap’ was used to generate the figure.

Results

Assessment of KAZN mRNA level in ovarian cancer, based on gene expression omnibus (GEO) datasets

The expression data of KAZN in ovarian cancer were obtained through the GEO database. A total of 11 microarrays from the GEO database met the entry criteria. The features of the selected GEO datasets are depicted in Table 1. Clinical features of samples in each dataset were listed in supplementary material Table S1. Expression of KAZN was significantly increased in ovarian cancer tissues in GSE105437, GSE18520, GSE27651, GSE36668, GSE38666, GSE40595, GSE66957, GSE69428, GSE54388, and GSE14407 (p = 0.0017, p = 0.0012, p < 0.0001, p = 0.0075, p < 0.0001, p < 0.0001, p < 0.0001, p = 0.0177, p = 0.0133, and p = 0.0049, respectively) (Fig. 1); no statistical difference was detected in GSE29450 (p = 0.0549, supplemental file: Fig. S1).

Table 1 Features of the enrolled Gene Expression Omnibus datasets
Fig. 1
figure 1

Expression of KAZN in ovarian cancer tissues and normal tissues based on Gene Expression Omnibus datasets. The expressions levels of KAZN are up-regulated in GEO datasets GSE105437, GSE18520, GSE27651, GSE36668, GSE38666, GSE40595, GSE66957, GSE69428, GSE54388, and GSE14407

Meta-analysis of GEO datasets

Meta-analysis was conducted based on 11 included microarrays from the GEO database (Fig. 2A). Given the apparent heterogeneity (p = 0.01, I2 = 56%), a random-effects model was applied, and remarkable up-regulation (SMD = 1.18, 95% CI: 0.76, 1.61) of KAZN mRNA was found in ovarian cancer group.

Fig. 2
figure 2

Meta-analysis of Gene Expression Omnibus (GEO) data. A Forest plot of GEO chips. The standard mean deviation is 0.99 (95% CI: 0.76, 1.61) with great heterogeneity (I2 = 56%, p = 0.01) showing that KAZN expression was markedly up-regulated in the ovarian cancer tissues. B Sensitivity analysis of GEO chips (p < 0.001). C A funnel plot for evaluating the publication bias of GEO chips (z = − 0.86, p = 0.3918)

Sensitivity analysis was performed to explore whether a particular microarray played a vital role in significant heterogeneity. By removing an individual study per time of meta-analysis to assess the influence of each study, the result showed that no study was found to have played a crucial role in any of the enrolled studies (overall effect p < 0.001, Fig. 2B). A funnel plot showed that no evidence of publication bias was observed for this analysis (p = 0.3918, Fig. 2C).

Validation of KAZN expression by qRT-PCR

We used qRT-PCR to detect the expression levels of the KAZN in 14 human samples (8 normal ovarian tissues and 6 tumor tissues). The results showed that KAZN expression was significantly upregulated in ovarian cancer tissue compared with normal control (p = 0.0087, Fig. 3), which is. consistent with the results of GEO database.

Fig. 3
figure 3

Quantitative real-time PCR analysis the expression of KAZN in human ovarian cancer tissues and normal ovarian tissues. ΔCt was used to determine the relative amounts of mRNA. There is Significant difference in 8 normal ovarian tissues and 6 tumor tissues by ANOVA statistical test, p = 0.0087

Up-regulation of KAZN affects overall survival in ovarian cancer

TCGA ovarian cancer datasets have only cancer tissue samples, no normal tissue as a control, so we compared them to GTEx samples, which have expression data from normal ovary tissue of GTEx donors who did not have cancer. To eliminate the batch effect, we conducted ComBat-seq to integrate different sourced datasets [24]. The comprehensive analysis based on TCGA, GTEx, and GEO datasets showed that the KAZN mRNA expression was significantly differential, and the expression was up-regulated in the TCGA group compared to GTEx group (p < 0.0001, Fig. 4A).

Fig. 4
figure 4

KAZN expression in ovarian cancer and survival analysis based on TCGA database. A Box plot, illustrating median expression levels of KAZN in normal ovarian tissue (Left, from GTEx) and ovarian tumors (Right, from TCGA). B Kaplan-Meier survival curve, illustrating overall survival for patients who had a tumor with up-regulated or down-regulated KAZN

To further study the clinical effects of KAZN in ovarian cancer, we divided the cases from TCGA OC datasets into four groups by the quantiles of KAZN counts, use quantile 1 and 4 as the KAZN high expression group and KAZN low expression group, and analyzed the survival status between the two groups. A Kaplan-Meier curve was used to identify the effects of the expression of KAZN on survival time and showed that the KAZN low expression group had significantly longer survival time than the KAZN high expression group (Fig. 4B). Univariate Cox regression analysis showed the association of KAZN mRNA with overall survival, demonstratig KAZN as a risk factor for OC (HR, 1.19, 95% CI, 1.03-1.37, p = 0.015, Fig. S2). Time-dependent ROC analysis indicated the prognostic accuracies were 0.652 at 6.8 years (Fig. S3).

KAZN protein was up-regulated in ovarian cancer

The Clinical Proteomic Tumor Analysis Consortium (CPTAC) has produced huge amounts of cancer proteomics data providing unprecedented research opportunities. The dataset from CPTAC Ovarian Cancer Confirmatory Study contains 41 normal participants and 169 tumor participants, which was used to validate the KAZN protein expression in Ovarian cancer. The result showed that in the tumor group, the KAZN protein expression was significantly higher than in the normal group (p = 0.0027, Fig. 5). KAZN protein expression was consistent with the mRNA level in this study.

Fig. 5
figure 5

KAZN protein level in ovarian cancer. The box plot illustrates the median of the KAZN protein in normal ovarian tissue (Left) and ovarian tumors (Right) based on CPTAC datasets

Correlation of KAZN methylation level with mRNA expression in ovarian cancer

Methylated CpG sites have a moderate to strong associations with gene expression changes across the phases in cancer-involved genes with specific functions. To study the correlation between the expression of KAZN and DNA methylation, we detected the DNA methylation level of CpG sites in the KAZN gene body region in TCGA datasets. From the TCGA database, we obtained 27 K DNA methylation array data of ovarian cancer, which contains 601 tumors and 12 normal samples. We detected the correlation between the expression of KAZN and cg17657618 and found the methylated CpG site cg17657618 positively correlated with the rising expression of the KAZN gene (Fig. 6A).

Fig. 6
figure 6

Correlation between KAZN methylation and mRNA expression and Correlation between KAZN methylation and clinical stage. A This plot illustrates that KAZN mRNA is positively correlated to cg17957618 methylation level (R = 0.29, p = 2.1e-08). B Boxplot for the cg17657618 methylation level between “stage I and II” group and “stage III and IV” group

To detect the correlation between KAZN DNA methylation and the progress of OC, we divided all cases into two groups by stage of clinical traits. We compared the methylation level between the two groups. The result showed that the cg17657618 methylation level of the stage “I and II” group was significantly lower than the stage “III and IV” group (Fig. 6B). Survival analysis of cg17657618 showed a trend of longer survival time for the low CpG beta value group, though not statistically significant (p = 0.089) (Fig. S4), suggesting that KAZN methylation may have an important role in the development and outcome of ovarian cancer.

The diagnostic value of KAZN methylation status in ovarian cancer

Illumina Infinium 450 K BeadChip covers 134 CpG sites for KAZN. The overview of all CpG sites for KAZN is shown in Fig. 7. In the GSE146552 and GSE81224 datasets, we found 13 and 12 differentially methylated CpG sites in the KAZN gene body region, respectively (details in supplemental Table S2). An unsupervised hierarchical clustering based on the 9 both different methylated CpG sites (cg00763594, cg00866953, cg01567509, cg02927252, cg08468082, cg13502395, cg14976342, cg21581845, cg27538859) in the two datasets in KAZN was constructed. The 9 CpG sites methylation pattern almost perfectly divided the samples of GSE146552 and GSE81224 into two groups, the tumor cluster, and the non-tumor cluster (Fig. 8), suggesting that the CpG pattern may potentially be used as a new biomarker for the diagnosis of ovarian cancer.

Fig. 7
figure 7

All CpG sites in KAZN gene body based on Illumina Infinium 450 K BeadChip. The upward ball-bar shape represents the CpG site in the forward strand. The downward ball-bar represents the CpG site in the reverse strand. The length of the bar represents the average methylation beta value. The red CpG sites are differentially methylated in the GSE146552 and GSE81224 datasets

Fig. 8
figure 8

Pheatmap of 9 KAZN methylation CpG sites. The tumor (epithelial ovarian cancer) group and non-tumor (non-tumor, includes epithelial layer of normal ovary and fimbriae fallopian tubes) group have different methylation patterns

Discussions

Previous studies have shown that KAZN is involved in multiple biological processes during development, such as cell proliferation and differentiation [17], as well as apoptosis [18], which prompted us to attempt identifying any possible associations of KAZN defects or abnormal expression levels with queer behaviors of cells, e.g., malignancies. To date, most of KAZN related studies have been focused on keratinocytes [15, 17, 27]. However, as KAZN is also expressed in many other tissues, its roles in health and diseases in general need to be investigated.

In the present study, we demonstrated that the expression of KAZN was significantly associated with OC. We compared 11 GEO microarray datasets to detect KAZN expression at the mRNA level and also conducted meta-analysis. KAZN mRNA was significantly up-regulated in OC in 10 GEO datasets, including 9 serous ovarian cancers and one unspecified ovarian cancer. In the remaining one dataset, which is clear cell cancer (supplemental file: Fig. S1), we did not find significant differences in KAZN expression, suggesting that the expression of KAZN may be correlated to the histological subtype of OC.

One of the merit of this work is the use of multiple data sources, which all indicated that KAZN was differentially expressed both on the mRNA level and on the protein level between ovarian cancer and normal tissues, with KAZN mRNA expression negatively correlated with survival time of the patients. KAZN can dually regulate proliferation and differentiation by Rho-dependent and –independent mechanisms [17], and also plays an important role in regulating cellular apoptosis by interacting with ARC and Bax [18]. During invasion and metastasis, cancer cells undergo changes in morphology and disruptions of cell connections. Overexpression of KAZN in keratinocytes can cause changes in cell shape and impair the assembly of intercellular junctions. As a component of the desmosome, KAZN also participates in the formation of cell connections [15, 17]. A recent study suggests that KAZN F is highly expressed in human cervical cancer tissues and could promote cell proliferation, migration and invasion in vitro by inhibiting apoptosis and facilitating epithelial-to-mesenchymal transition (EMT) [18]. It was further proved that KAZN F was directly regulated by miR-186, which influences the sensitivity of ovarian cancer cells to paclitaxel and cisplatin [28]. As such, it is reasonable to speculate that KAZN gene may participate in tumorigenesis, invasion, and metastasis of OC by affecting these key processes.

Growing evidence has shown that epigenetic changes are involved in cancer development and progression [22]. Indeed, DNA hypermethylation can lead to the silencing of tumor suppressor genes, whilst DNA hypomethylation can induce genomic instability, increasing transcription and facilitating protein activation [29]. A large number of putative suppressor genes that are silenced or activated by aberrant methylation have been identified in ovarian cancer, including hypermethylation of OPCML, RASSF1, CDKN2B as well as classical tumor suppressors BRCA1, p16 and MLH1, and hypomethylation of LINE-1, SLC6A12 and PRAME [30,31,32,33,34,35,36,37,38]. A recent study reported the association of CpG site cg17657618 with Endometriosis [39], but whether this CpG site may be associated with ovarian cancer was unknown. Our results based on TCGA human methylation 27 K BeadChips showed that cg17657618 was hypermethylated and positively correlated with the expression of the KAZN in OC. Analysis of the 450 K array dataset revealed 9 methylation sites that were differentially expressed in both datasets. The differentially expressed CpG sites, including 5 hypermethylated and 4 hypomethylated CpG sites in OC samples, may affect KAZN expression in different ways, whereby participating in the formation and progression of ovarian cancer. These findings all support the critical roles of KAZN gene expression and its methylation in OC occurrence and progress.

It is a common practice to use cluster analysis to divide patients into different groups based on the expression values of multiple genes or the top-ranked methylation CpG sites ß-value. A recent study discriminated EOC from normal ovarian tissues with high specificity and sensitivity based on the methylation of RASSF1A, OPCML and HOXA9 [40]. In another study, methylation of OPCML, together with methylation of RUNX3 and TFPI2, was demonstrated to be an early diagnostic marker for OC with higher sensitivity and specificity than classical CA125 [41]. However, our results based on single gene methylation CpG site ß-value were peculiar: the CpG methylation pattern of KAZN could be a new biomarker for predicting and diagnosing OC.

Previous studies have shown that KAZN has six splice variants, which are differentially expressed in a wide variety of cell types [15, 18]. As the transcripts are mostly tissue-specific and expressed in a temporal sequence, multiple methylation sites in the single gene are expected to play different roles at different stages in the development of ovarian cancer. In contrast to genetic changes, the relatively reversible character of the epigenetic alterations like DNA methylation determines that it has the potential to be artificially regulated. DNA methyltransferases (DNMTs) are responsible for the establishment and maintenance of the DNA methylation patterns on human genome. Hence, with continuous discovery of various global and specific DNMT inhibitors, CpG methylation sites of KAZN could be a new potential therapeutic target for ovarian cancer treatment. Further understanding of the function of the KAZN gene and the mechanistic relationship between KAZN expression and methylation in cancer will open a new horizon for the control of malignant diseases.

Availability of data and materials

The dataset supporting the conclusions of this article is included within the article and its supplemental tables. Methylation and mRNA expression of ovarian cancer tissues were mined from The Cancer Genome Atlas (TCGA)(https://cancergenome.nih.gov/) and Gene expression omnibus (GEO, https://www.ncbi.nlm.nih.gov/geo/. The accession number is GSE105437, GSE18520, GSE27651, GSE36668, GSE38666, GSE40595, GSE66957, GSE69428, GSE54388, GSE14407, GSE137238, GSE101108, GSE101948, GSE143897, GSE132107, GSE146552 and GSE81224, respectively).

References

  1. Guo Q, He Y, Sun L, Kong C, Cheng Y, Zhang G. In silico detection of potential prognostic circRNAs through a re-annotation strategy in ovarian cancer. Oncol Lett. 2019;17(4):3677–86.

    CAS  PubMed  PubMed Central  Google Scholar 

  2. Wang Q, Lu Z, Ma J, Zhang Q, Wang N, Qian L, et al. Six-mRNA risk score system and nomogram constructed for patients with ovarian cancer. Oncol Lett. 2019;18(2):1235–45.

    PubMed  PubMed Central  Google Scholar 

  3. Hou Y, Wang Y, Xu S, Qi G, Wu X. Bioinformatics identification of microRNAs involved in polycystic ovary syndrome based on microarray data. Mol Med Rep. 2019;20(1):281–91.

    CAS  PubMed  PubMed Central  Google Scholar 

  4. Zhuang H, Cheng L, Wang Y, Zhang YK, Zhao MF, Liang GD, et al. Dysbiosis of the gut microbiome in lung Cancer. Front Cell Infect Microbiol. 2019;9:112.

    CAS  Article  Google Scholar 

  5. Kulbe H, Otto R, Darb-Esfahani S, Lammert H, Abobaker S, Welsch G, et al. Discovery and validation of novel biomarkers for detection of epithelial ovarian Cancer. Cells Basel. 2019;8(7):713.

  6. Jiang Y, Lyu T, Che X, Jia N, Li Q, Feng W. Overexpression of SMYD3 in ovarian Cancer is associated with ovarian Cancer proliferation and apoptosis via Methylating H3K4 and H4K20. J Cancer. 2019;10(17):4072–84.

    CAS  Article  Google Scholar 

  7. Zheng J, Zhou Y, Li XJ, Hu JM. MiR-574-3p exerts as a tumor suppressor in ovarian cancer through inhibiting MMP3 expression. Eur Rev Med Pharmacol Sci. 2019;23(16):6839–48.

    CAS  PubMed  Google Scholar 

  8. Harrington BS, Annunziata CM. NF-kappaB Signaling in Ovarian Cancer. Cancers (Basel). 2019;11:8.

    Article  Google Scholar 

  9. Mohamed NE, Hay T, Reed KR, Smalley MJ, Clarke AR. APC2 is critical for ovarian WNT signalling control, fertility and tumour suppression. BMC Cancer. 2019;19(1):677.

    Article  Google Scholar 

  10. Li L, Gu H, Chen L, Zhu P, Zhao L, Wang Y, et al. Integrative network analysis reveals a MicroRNA-based signature for prognosis prediction of epithelial ovarian Cancer. Biomed Res Int. 2019;2019:1056431.

    PubMed  PubMed Central  Google Scholar 

  11. Zhao J, Huang J, Geng X, Chu W, Li S, Chen ZJ, et al. Polycystic ovary syndrome: novel and hub lncRNAs in the insulin resistance-associated lncRNA-mRNA network. Front Genet. 2019;10:772.

    CAS  Article  Google Scholar 

  12. Gao YC, Zhou XH, Zhang W. An ensemble strategy to predict prognosis in ovarian Cancer based on gene modules. Front Genet. 2019;10:366.

    CAS  Article  Google Scholar 

  13. Wang L, Ni S, Du Z, Li X. A six-CpG-based methylation markers for the diagnosis of ovarian cancer in blood. J Cell Biochem. 2020;121(2):1409–19.

  14. Mase S, Shinjo K, Totani H, Katsushima K, Arakawa A, Takahashi S, et al. ZNF671 DNA methylation as a molecular predictor for the early recurrence of serous ovarian cancer. Cancer Sci. 2019;110(3):1105–16.

    CAS  Article  Google Scholar 

  15. Groot KR, Sevilla LM, Nishi K, DiColandrea T, Watt FM. Kazrin, a novel periplakin-interacting protein associated with desmosomes and the keratinocyte plasma membrane. J Cell Biol. 2004;166(5):653–9.

    CAS  Article  Google Scholar 

  16. Mu XQ, Liu BB, Hui E, Huang W, Yao LC, Duo LB, et al. A rapid loop-mediated isothermal amplification (LAMP) method for detection of the macrolide-streptogramin type B resistance gene msrA in Staphylococcus aureus. J Glob Antimicrob Resist. 2016;7:53–8.

    Article  Google Scholar 

  17. Sevilla LM, Nachat R, Groot KR, Watt FM. Kazrin regulates keratinocyte cytoskeletal networks, intercellular junctions and differentiation. J Cell Sci. 2008;121(Pt 21):3561–9.

    CAS  Article  Google Scholar 

  18. Zhou YJ, Zhu S, Yang DH, Zhao DD, Li JJ, Liu SL. Characterization of Klebsiella sp. strain S1: a bacterial producer of secoisolariciresinol through biotransformation. Can J Microbiol. 2017;63(1):1–10.

    CAS  Article  Google Scholar 

  19. Chhatriwala MK, Cipolat S, Sevilla LM, Nachat R, Watt FM. Exons 5-15 of kazrin are dispensable for murine epidermal morphogenesis and homeostasis. J Invest Dermatol. 2012;132(8):1977–87.

    CAS  Article  Google Scholar 

  20. Li QH, Jin G, Wang JY, Li HN, Liu H, Chang XY, et al. Live attenuated Salmonella displaying HIV-1 10E8 epitope on fimbriae: systemic and mucosal immune responses in BALB/c mice by mucosal administration. Sci Rep. 2016;6:29556.

    CAS  Article  Google Scholar 

  21. Li FF, Han Y, Shi S, Li X, Zhu XD, Zhou J, et al. Characterization of transcriptional repressor gene MSX1 variations for possible associations with congenital heart diseases. PLoS One. 2015;10(11):e0142666.

    Article  Google Scholar 

  22. Dong J, Xu M. A 19miRNA support vector machine classifier and a 6miRNA risk score system designed for ovarian cancer patients. Oncol Rep. 2019;41(6):3233–43.

    CAS  PubMed  PubMed Central  Google Scholar 

  23. Xie H, Xu H, Hou Y, Cai Y, Rong Z, Song W, et al. Integrative prognostic subtype discovery in high-grade serous ovarian cancer. J Cell Biochem. 2019;120(11):18659–66.

    CAS  Article  Google Scholar 

  24. Luo M, Zhou L, Zhan SJ, Cheng LJ, Li RN, Wang H, et al. ALPL regulates the aggressive potential of high grade serous ovarian cancer cells via a non-canonical WNT pathway. Biochem Biophys Res Commun. 2019;513(2):528–33.

    CAS  Article  Google Scholar 

  25. Zhao H, Sun Q, Li L, Zhou J, Zhang C, Hu T, et al. High expression levels of AGGF1 and MFAP4 predict primary platinum-based Chemoresistance and are associated with adverse prognosis in patients with serous ovarian Cancer. J Cancer. 2019;10(2):397–407.

    CAS  Article  Google Scholar 

  26. Feng Y, Hang W, Sang Z, Li S, Xu W, Miao Y, et al. Identification of exosomal and nonexosomal microRNAs associated with the drug resistance of ovarian cancer. Mol Med Rep. 2019;19(5):3376–92.

    CAS  PubMed  PubMed Central  Google Scholar 

  27. Goicoechea SM, Zinn A, Awadia SS, Snyder K, Garcia-Mata R. A RhoG-mediated signaling pathway that modulates invadopodia dynamics in breast cancer cells. J Cell Sci. 2017;130(6):1064–77.

    CAS  PubMed  PubMed Central  Google Scholar 

  28. Li T, Ferns K, Yan ZQ, Yin SY, Kou JJ, Li D, et al. Acanthopanax senticosus: photochemistry and anticancer potential. Am J Chin Med. 2016;44(8):1543–58.

    CAS  Article  Google Scholar 

  29. Chen Y, Bi F, An Y, Yang Q. Identification of pathological grade and prognosis-associated lncRNA for ovarian cancer. J Cell Biochem. 2019;120(9):14444–54.

    CAS  Article  Google Scholar 

  30. Tian X, Guan W, Zhang L, Sun W, Zhou D, Lin Q, et al. Physical interaction of STAT1 isoforms with TGF-beta receptors leads to functional crosstalk between two signaling pathways in epithelial ovarian cancer. J Exp Clin Cancer Res. 2018;37(1):103.

    Article  Google Scholar 

  31. Ferrero H, Corachan A, Aguilar A, Quinonero A, Carbajo-Garcia MC, Alama P, et al. Single-cell RNA sequencing of oocytes from ovarian endometriosis patients reveals a differential transcriptomic profile associated with lower quality. Hum Reprod. 2019;34(7):1302–12.

    CAS  Article  Google Scholar 

  32. Bi X, Zhai Z, Wang S. Identification of the key pathways and genes related to polycystic ovary syndrome using bioinformatics analysis. Gen Physiol Biophys. 2019;38(3):205–14.

    CAS  Article  Google Scholar 

  33. Li J, Yue H, Yu H, Lu X, Xue X. Development and validation of SIRT3-related nomogram predictive of overall survival in patients with serous ovarian cancer. J Ovarian Res. 2019;12(1):47.

    Article  Google Scholar 

  34. Liu S, Matsuzaki J, Wei L, Tsuji T, Battaglia S, Hu Q, et al. Efficient identification of neoantigen-specific T-cell responses in advanced human ovarian cancer. J Immunother Cancer. 2019;7(1):156.

    Article  Google Scholar 

  35. Um H, Tixier F, Bermudez D, Deasy JO, Young RJ, Veeraraghavan H. Impact of image preprocessing on the scanner dependence of multi-parametric MRI radiomic features and covariate shift in multi-institutional glioblastoma datasets. Phys Med Biol. 2019;64(16):165011.

    Article  Google Scholar 

  36. Josahkian JA, Saggioro FP, Vidotto T, Ventura HT. Candido dos Reis FJ, de Sousa CB, Tiezzi DG, de Andrade JM, Koti M, squire JA: increased STAT1 expression in high grade serous ovarian Cancer is associated with a better outcome. Int J Gynecol Cancer. 2018;28(3):459–65.

    Article  Google Scholar 

  37. Li N, Zhan X. Signaling pathway network alterations in human ovarian cancers identified with quantitative mitochondrial proteomics. EPMA J. 2019;10(2):153–72.

    Article  Google Scholar 

  38. Ma SY, Wei P, Qu F. KCNMA1-AS1 attenuates apoptosis of epithelial ovarian cancer cells and serves AS a risk factor for poor prognosis of epithelial ovarian cancer. Eur Rev Med Pharmacol Sci. 2019;23(11):4629–41.

    PubMed  Google Scholar 

  39. Wang X, Han L, Zhou L, Wang L, Zhang LM. Prediction of candidate RNA signatures for recurrent ovarian cancer prognosis by the construction of an integrated competing endogenous RNA network. Oncol Rep. 2018;40(5):2659–73.

    CAS  PubMed  PubMed Central  Google Scholar 

  40. Hinchcliff E, Paquette C, Roszik J, Kelting S, Stoler MH, Mok SC, et al. Lymphocyte-specific kinase expression is a prognostic indicator in ovarian cancer and correlates with a prominent B cell transcriptional signature. Cancer Immunol Immunother. 2019;68(9):1515–26.

    CAS  Article  Google Scholar 

  41. Wu L, Shang W, Zhao H, Rong G, Zhang Y, Xu T, et al. In Silico screening of circulating MicroRNAs as potential biomarkers for the diagnosis of ovarian Cancer. Dis Markers. 2019;2019:7541857.

    PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We highly appreciate the time and efforts of the researchers who provided invaluable data to the public databases and who maintain the public databases.

Funding

This work was supported by the grants from National Natural Science Foundation of China (NSFC81871623, 82020108022).

Author information

Authors and Affiliations

Authors

Contributions

Zhu S and Liu SL conceived and designed the study. Zhu S and Zhang M contributed to the computational analyses and confirmed the results. Zhu S and Bao H wrote the main manuscript text. Zhao X gave advice on software applications. All authors have read and approved the final version of this manuscript.

Corresponding author

Correspondence to Shu-Lin Liu.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Ethics Committee of Harbin medical university.

Consent for publication

Not applicable.

Competing interests

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.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1: Supplementary Table 1.

Clinical feature information of samples in each dataset.

Additional file 2: Supplementary Table 2.

Differentially methylated CpG sites in the KAZN gene body region.

Additional file 3: Supplementary Figure 1.

Expression of KAZN in ovarian cancer tissues and normal tissues based on Gene Expression Omnibus datasets. The expressions of KAZN are not significantly differential in GSE29450.

Additional file 4: Supplementary Figure 2.

Forest plots of hazard ratios and their confidence intervals for the KAZN mRNA.

Additional file 5: Supplementary Figure 3.

The time-dependent receiver operating characteristic (ROC) curve. The area under ROC curve (AUC) reached 0.652 at 6.8 year.

Additional file 6: Supplementary Figure 4.

The Kaplan-Meier survival curve illustrate overall survival for patients who had tumor with hypermethylated or hypomethylated cg17657618.

Rights and permissions

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

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zhu, S., Bao, H., Zhang, MC. et al. KAZN as a diagnostic marker in ovarian cancer: a comprehensive analysis based on microarray, mRNA-sequencing, and methylation data. BMC Cancer 22, 662 (2022). https://doi.org/10.1186/s12885-022-09747-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-022-09747-2

Keywords

  • KAZN
  • Ovarian cancer
  • Methylation
  • Bioinformatics
  • Survival