SYSTEMATIC REVIEW article

Front. Cell Dev. Biol., 11 March 2021

Sec. Epigenomics and Epigenetics

Volume 9 - 2021 | https://doi.org/10.3389/fcell.2021.570898

OPCML Methylation and the Risk of Ovarian Cancer: A Meta and Bioinformatics Analysis

  • 1. Nanjing Maternity and Child Health Care Institute, Nanjing Maternity and Child Health Care Hospital, Women's Hospital of Nanjing Medical University, Nanjing, China

  • 2. The First People's Hospital of Zhangjiagang City, The Zhangjiagang Affiliated Hospital of Soochow University, Suzhou, China

  • 3. Department of Radiotherapy, The Affiliated Cancer Hospital of Nanjing Medical University, Nanjing, China

  • 4. Jiangsu Institute of Cancer Research, The Affiliated Cancer Hospital of Nanjing Medical University, Nanjing, China

Article metrics

View details

8

Citations

3,9k

Views

2k

Downloads

Abstract

Background: The association of opioid binding protein cell adhesion molecule-like (OPCML) gene methylation with ovarian cancer risk remains unclear.

Methods: We identified eligible studies by searching the PubMed, Web of Science, ScienceDirect, and Wanfang databases. Odds ratios (ORs) and 95% confidence intervals (95% CIs) were used to determine the association of OPCML methylation with ovarian cancer risk. Meta-regression and subgroup analysis were used to assess the sources of heterogeneity. Additionally, we analyzed the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) datasets to validate our findings.

Results: Our study included 476 ovarian cancer patients and 385 controls from eight eligible studies. The pooled OR was 33.47 (95% CI = 12.43–90.16) in the cancer group vs. the control group under the random-effects model. The association was still significant in subgroups according to sample type, control type, methods, and sample sizes (all P < 0.05). Sensitivity analysis showed that the finding was robust. No publication bias was observed in Begg's (P = 0.458) and Egger's tests (P = 0.261). We further found that OPCML methylation was related to III/IV (OR = 4.20, 95% CI = 1.59–11.14) and poorly differentiated grade (OR = 4.37; 95% CI = 1.14–16.78). Based on GSE146552 and GSE155760, we validated that three CpG sites (cg16639665, cg23236270, cg15964611) in OPCML promoter region were significantly higher in cancer tissues compared to normal tissues. However, we did not observe the associations of OPCML methylation with clinicopathological parameters and overall survival based on TCGA ovarian cancer data.

Conclusion: Our findings support that OPCML methylation is associated with an increased risk of ovarian cancer.

Introduction

Ovarian cancer is the most frequent cause of death in women with gynecological malignancies. In 2020, 313,959 new ovarian cancer cases were diagnosed and 207,252 new deaths occurred from ovarian cancer in the world (Sung et al., 2021). Most patients are diagnosed at an advanced or incurable stage because of asymptomatic development (Lheureux et al., 2019). Therefore, identifying novel biomarkers for early diagnosis and research on treatment strategies could reduce the incidence of ovarian cancer and improve the survival rate in advanced ovarian cancer patients.

Aberrant gene expression, including the activation of oncogenes and inactivation of tumor suppressor genes (TSGs), plays a vital role in tumorigenesis (Ehrlich, 2002). Emerging studies have shown that DNA methylation is a key mechanism of epigenetic variability in gene expression (Wilson et al., 2007). Compared with normal cells, cancer cells exhibit extensive changes in DNA methylation patterns (Vidal et al., 2017). Cancer cells undergo alterations in 5-methylcytosine distribution including global DNA hypomethylation and hypermethylation of CpG islands (Jones, 2012). Notably, aberrant DNA methylation is not only a well-known mechanism of TSG inactivation but also one of the earliest events in carcinogenesis (Rauscher et al., 2015; Widschwendter et al., 2017). DNA methylation is considered a promising biomarker for cancer diagnosis and the prediction of treatment and prognosis (Szyf, 2012; Widschwendter et al., 2017). Intriguingly, the Food and Drug Administration (FDA) has approved the first methylation-based assay in colorectal cancer screening and early detection (Song et al., 2017; Widschwendter et al., 2018). Therefore, evaluating the alteration of DNA methylation in ovarian cancer is urgent.

Opioid binding protein cell adhesion molecule-like (OPCML) gene at 11q25 is a glycosyl-phosphatidylinositol (GPI)-anchored cell adhesion-like molecule and belongs to the IgLON family (Wu and Sood, 2012). OPCML demonstrates tumor-suppressor function in epithelial ovarian cancer both in vitro and in vivo (Sellar et al., 2003). Subsequently, the same team described the mechanism underlying the phenotype of OPCML. It negatively regulates a specific receptor tyrosine kinase (RTK) repertoire comprising erythropoietin-producing hepatocellular receptor-2 (EPHA2), fibroblast growth factor receptor-1 (FGFR1), fibroblast growth factor receptor-3 (FGFR3), human epidermal growth factor receptor-2 (HER2), and human epidermal growth factor receptor-4 (HER4) receptors by binding to the extracellular domains of RTKs, thus promoting their degradation via a polyubiquitination-associated proteasomal mechanism and leading to growth inhibition (McKie et al., 2012). Additionally, they revealed that a recombinant protein based on OPCML is an active anticancer agent preclinical in vivo. Another study showed that OPCML restoration interfered with HER2-EGFR heterodimer formation and potentiated lapatinib and erlotinib therapy in ovarian and breast cancer cell lines overexpressing HER2 (Zanini et al., 2017). Recent researches on OPCML protein reveal that OPCML is a potential anti-cancer therapy (Birtley et al., 2019; Simovic et al., 2019). Considering the biological and clinical effects of OPCML, research on the mechanism of aberrant expression patterns may benefit ovarian cancer patients from OPCML-based therapy.

Many studies have indicated that OPCML DNA hypermethylation frequently occurs in ovarian cancers. However, the findings remain unclear regarding the small population and different methods in previous studies (Czekierdowski et al., 2006; Zhang et al., 2006, 2014; Chen et al., 2007; Liu et al., 2008; Zhou et al., 2011, 2014a,b; Wang et al., 2015, 2017; Xing et al., 2015). Therefore, we performed a meta-analysis to identify the association of OPCML methylation with ovarian cancer risk. We further assessed whether OPCML methylation was associated with other clinicopathological parameters in eligible studies, such as age, stage, histological type, and tumor differentiation. Additionally, we validated our results based on the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) ovarian cancer datasets.

Materials and Methods

Literature Search Strategy

A systematic literature search was performed in the PubMed, Web of Science, ScienceDirect, and Wanfang databases using the following keywords and search items: (ovarian OR ovary) AND (cancer OR carcinoma OR tumor) AND (OPCML methylation). The search was updated until November 17, 2020.

Study Selection

The inclusion criteria were as follows: (1) studies primarily evaluating the association between OPCML methylation and ovarian cancer, (2) studies including case and control populations, (3) the incidence of OPCML methylation in case and control groups, and (4) studies using human tissues or blood. Based on the inclusion criteria, the titles and abstracts from the preliminary search were evaluated. Subsequently, all related studies were evaluated as full-text papers. Exclusion criteria were (1) book section, conference abstracts, and reviews, (2) studies with insufficient data to provide the methylation frequencies in cancer and control groups, and (3) studies without human samples. In the case of duplicated publications from the shared cohort, we selected the most complete information to be included in the meta-analysis.

Data Extraction and Quality Assessment

The following information was extracted for each eligible study: the first author's name, publication year, country, age, number of methylated and unmethylated samples in cancer and control groups, sample type, control type, methylation detection methods, and other clinicopathological parameters. In our study, controls included adjacent non-cancerous ovarian tissues (AT), tissue and blood samples from benign ovarian tumor (BOT) and other cancer-free patients (NT). Additionally, we used the Newcastle-Ottawa Scale (NOS) to assess the quality. The NOS includes the selection of the research groups (four stars), comparison of the groups (two stars), and ascertainment of the outcome (three stars).

Bioinformatics Analysis

Two DNA methylation data were downloaded from the Gene Expression Omnibus (GEO) database: GSE146552 (20 high grade serous ovarian cancers, 16 epithelial layer by mechanical scraping of resected fimbriae fallopian tubes and ovaries) and GSE155760 (23 high-grade serous ovarian carcinoma, 11 fallopian tube mucosa from cancer-free normal control). The GEO2R online analysis tool (https://www.ncbi.nlm.nih.gov/geo/geo2r/) was used to identify the differential DNA methylated sites (P < 0.05) between cancer tissues and normal tissues. For the cancer genome atlas (TCGA) ovarian serous cystadenocarcinoma (OV) data, we downloaded the OPCML expression data (log2-transformed RSEM normalized count) based on RNA sequencing (RNA-Seq), OPCML-associated DNA methylation data (cg03923934, cg25853078) based on HumanMethylation27 and the clinical data including age, stage, grade, and overall survival from the website (https://xena.ucsc.edu/).

Statistical Analysis

We calculated the pooled odds ratios (ORs) and 95% confidence intervals (95% CIs) to evaluate the association between OPCML methylation and ovarian cancer risk. The heterogeneity of studies was evaluated by the χ2-based Q-test. When the P > 0.05, the fixed-effects model was used to combine the effect size; otherwise, the random-effects model was adopted. We conducted meta-regression to explore sources of heterogeneity and further performed a subgroup analysis to evaluate the source of the heterogeneity. The contribution of each study to our findings was examined according to sensitivity analysis. We used the Begg's rank correlation method (Begg and Mazumdar, 1994) and a funnel plot for Egger's test (Egger et al., 1997) to evaluate publication bias.

The correlation between OPCML methylation and its expression in TCGA was calculated by Spearman's rank correlation. Kruskal-Wallis and Wilcoxon-Mann-Whitney U tests were applied for comparisons, as appropriate. Survival analysis was performed using the log-rank test and a Cox regression model. Hazard ratios (HRs) and corresponding 95% CIs were calculated to assess the associations of factors with overall survival. All statistical tests were two-sided and performed using the R software (version 3.6). P-value < 0.05 were considered statistically significant.

Results

Study Characteristics

The detailed information of how we selected relevant articles is shown in Figure 1. We initially identified 126 articles from PubMed, Web of Science, ScienceDirect, and Wanfang databases and excluded 34 duplicated articles. After reviewing 92 articles based on titles and abstracts, 75 articles did not meet the selection criteria, leaving 17 studies for detailed full-text evaluation. Seven articles were excluded because they did not report the methylation frequencies in cancers or controls. Three articles were excluded because of overlapping cohorts (Zhou et al., 2011; Zhang et al., 2014; Wang et al., 2017). The study performed by Zhou et al. (2014a) included two sample types, tissue and serum. Thus, we considered this paper as two separate studies. Therefore, seven eligible articles with eight studies including 476 ovarian cancer patients and 385 controls were included in the meta-analysis.

Figure 1

Figure 1

Flowchart of study selection in the meta-analysis.

The characteristics of the eight studies are summarized in Table 1. Seven studies were performed in China and one study was performed in Poland. The cases comprised of cancer tissues and blood from ovarian cancer patients. Among the eight included studies, six studies clearly described that the cases were primary tumor patients whereas two studies (Wang et al., 2015; Xing et al., 2015) did not. In the study performed by Liu et al. (2008), there are 41 metastatic tissues of pelvic and abdomen cavity. We provided this in Table 1. However, we did not include these samples when we calculated the pooled OR. The controls in these studies were not matched to cases for any features. OPCML promoter methylation was detected in all eligible articles. There are four studies that used methylation-specific polymerase chain reaction (MSP), three used restriction enzyme-related analysis [methylation-sensitive restriction enzyme-polymerase chain reaction (MSRE-PCR) and restriction enzyme cut analysis], and one study used cationic conjugated polymer (CCP)-based fluorescence resonance energy transfer (FRET) to detect OPCML methylation in ovarian cancer and controls. According to the NOS criteria, the quality of these studies was six, indicating a relatively high quality. Detailed features are shown in Supplementary Table 1.

Table 1

ReferencesCountryAge (mean or median)Cancer groupControl groupSample typeControl typeMethodMethylation frequency of OPCML associated with other clinicopathological parameters
MUMU
Wang et al. (2015)ChinaCancer, 54.06 years. Benign, 33.95 years. Healthy, 46.17 years.64711a112aSerumBOT&NTMSPFIGO stage: I/II (33/39), III/IV (30/32).
Xing et al. (2015)ChinaNR27829TissuesNTCCP-based FRETAge: <60 years (23/28), ≥60 years (4/7). FIGO stage: I/II (5/11), III/IV (22/24). Histological type: serous (22/26), others (5/9)b. Grade of differentiationc: High-grade serous (19/21), others (8/14) Ascites: positive (21/24), Negative (6/11). CA125: >200 U/mL (21/25), <200 U/mL (6/10)
Zhou et al. (2014a)ChinaCancer and normal, 50 years.3960d40dTissuesNTMSRE-PCRTissues FIGO stage: I/II (5/10), III/IV (34/35). Histological type: serous (18/20), others (21/25)e. Tumor differentiation: Well and moderately differentiated (10/16), Poorly differentiated (29/29).
Healthy, 50 years.3690d20dSerumNTMSRE-PCRSerum FIGO stage: I/II (5/10), III/IV (31/35). Histological type: serous (17/20), others (19/25)e. Tumor differentiation: Well and moderately differentiated (9/16), Poorly differentiated (27/29).
Zhou et al. (2014b)ChinaCancer, 54 years. Benign, 53 years. Normal, 54 years.80220f30fTissuesNTMSPAge: <55 years (39/54), ≥55 years (41/48). Pathological stage*: I/II (31/44), III/IV (49/58). Histological typeb: serous (42/57), others (38/45).
––28f57fTissuesBOTMSP
Czekierdowski et al. (2006)PolandNRg202304TissuesNTMSPFIGO stage: I/II (4/5), III/IV (16/38). Histological type: serous (8/15), others (12/28)h. Histological Grading: G1+G2 (11/20), G3 (9/23)k
Zhang et al. (2006)ChinaCancer, 51.5 years. Benign, 39.0 years. Normal, 46.5 years.3240020TissuesNTRestriction enzyme cut
017TissuesBOTRestriction enzyme cut
Liu et al. (2008)iChinaCancer, 53.8 years. Normal, 53.2 years.3033020TissuesNTMSPFIGO stage: I/II (2/22), III/IV (16/41). Histological type: serous (20/34), others (10/29)j. Grading: Well and Moderately differentiated (6/36), Poorly differentiated (12/27)
213TissuesAT

Characteristics of included studies.

M, methylation; U, unmethylation; BOT, benign ovarian tissues or blood; NT, normal ovarian tissues or blood of cancer-free patients or healthy people; AT, Adjacent non-cancerous ovarian tissues; MSP, methylation-specific polymerase chain reaction; MSRE-PCR, methylation-sensitive restriction enzyme-polymerase chain reaction; CCP-based FRET, cationic conjugated polymer (CCP)-based fluorescence resonance energy transfer (FRET); FIGO, International Federation of Gynecology and Obstetrics. NR, not reported the median or mean age in all patients.

a

The control group include 43 benign ovarian carcinoma and 80 healthy women in this study.

b

Others included Mucinous, Clear cell, and Endometrioid carcinomas.

c

In further analysis, we regarded the high-grade serous as Poorly differentiated., others as well and moderately differentiated.

d

This study include 45 ovarian tumor tissues, 40 normal ovarian tissues, and 20 serum samples.

e

Others included Mucinous, Clear cell, Endometrioid, and undifferentiated carcinomas.

f

This study include 102 ovarian cancer tissues, 85 benign ovarian tumor, and 30 normal ovarian.

*

Pathological stage was classified according to the tumor, lymph node, and metastasis (TNM) classification of the American Joint Committee on Cancer (AJCC).

g

This study only reported the mean ages of different histological type, but not reported the mean age of cancer in tumor and normal tissue.

h

Others included Mucinous, Clear cell, Endometrioid, undifferentiated, metastatic, and other carcinomas.

i

This study include 63 epithelial ovarian cancer tissues, 41 Metastatic tissues of pelvic and abdomen cavity, 15 Adjacent non-cancerous ovarian tissues, 20 Normal ovairan tissues.

j

Others included Mucinous, Clear cell, and Endometrioid carcinomas.

k

In further analysis, we regarded the G1+G2 as well and moderately differentiated, G3 as poorly differentiated.

Relationship Between OPCML Methylation and Ovarian Cancer Risk

Because of significant heterogeneity among the included studies, a random-effects model was used to evaluate the effect size (P < 0.01; I2 = 65%; Figure 2). In the overall meta-analysis, the OPCML methylation frequency was significantly associated with an increased risk of ovarian cancer (summary OR: 33.47; 95% CI = 12.43–90.16; P < 0.001; Figure 2).

Figure 2

Figure 2

The summary estimates for OPCML methylation frequency associated with ovarian cancer in the meta-analysis.

Meta-Regression and Subgroup Analysis

To explore the potential sources of heterogeneity, we performed multiple regression model using six variables (publication year, country, sample type, control type, methods, and sample size). However, the source of heterogeneity was not observed among these factors (all P > 0.05; Table 2). We also conducted subgroup analysis to assess the source of the heterogeneity according to country, sample type, sample size, control type, and method (Table 3). In the subgroup analysis of country, the OR was 38.17 (95% CI = 13.21–110.28) in China under the random-effects model. The OR was 102.86 (95% CI = 39.37–268.74) under the fixed-effects model in the serum subgroup, and 18.22 (95% CI = 10.93–30.36) under the fixed-effects model in the tissues subgroup. In the subgroup analysis based on the control type, the OR was 52.77 (95% CI = 16.92–164.55) in the NT group under the fixed-effects model and 35.50 (95% CI = 6.83–197.47) in the NT and BOT groups under the random-effects model. In the subgroup analysis of methods, the OR was 21.78 (95% CI = 6.10–77.79) in the MSP group under the random-effects model and 139.45 (95% CI = 24.06–808.26) in the restriction enzyme-related analysis group under the fixed-effects model. The OR for the sample size subgroup was 34.19 (95% CI = 14.11–82.86) in the <100 subgroup under the fixed-effects model and 35.50 (95% CI = 6.38–197.47) in the ≥100 subgroup under the random-effects model.

Table 2

Heterogeneity sourcesCoefficient95% CIP
Publication year0.814(−0.751, 2.380)0.308
Country10.538(−9.721, 30.795)0.308
Sample type−0.716(−3.370, 1.939)0.597
Control type4.851(−4.042,13.744)0.285
Method3.748(−1.419, 8.916)0.155
Sample size−0.001(−0.020, 0.018)0.889

Meta-regression analysis.

Table 3

GroupStudiesCaseControlFixed effects modelRandom effects modelHeterogeneity
nMUMUOR (95% CI)OR (95% CI)I2 (%)P
Country
    China73081254333827.47 (17.52, 43.08)38.17 (13.21, 110.28)70<0.01
    Poland12023047.85 (0.40, 154.75)7.85 (0.40, 154.75)––
Sample type
     Serum21001611132102.86 (39.37, 268.74)98.42 (38.37, 252.45)00.73
     Tissues62281323221018.22 (10.93, 30.36)19.13 (7.71, 47.44)390.15
Control type
     NT41224627352.77 (16.92, 164.55)47.94 (7.38, 311.27)530.09
    NT&BOT3176693923623.41 (13.96, 39.25)35.50 (6.38, 197.47)84<0.01
     AT&NT1303323315.00 (3.31, 67.93)15.00 (3.31, 67.93)––
Methods
    MSP4194854123619.52 (11.94, 31.92)21.78 (6.10, 77.79)76<0.01
    CCP-based FRET12782915.19 (2.71, 85.10)15.19 (2.71, 85.10)––
    Restriction enzyme related methods310755097139.45 (24.06, 808.26)163.86 (31.13, 862.61)00.59
Sample size
    <100515279410634.19 (14.11, 82.86)32.63 (8.42, 126.53)450.12
    ≥1003176693923623.41 (13.96, 39.25)35.50 (6.38, 197.47)84<0.01

Subgroup analysis of associations between OPCML methylation and ovarian cancer risk.

M, methylation; U, unmethylation; BOT, benign ovarian tissues; NT, normal ovarian tissues of cancer-free patients or healthy people; AT, Adjacent non-cancerous ovarian tissues; MSP, methylation-specific polymerase chain reaction; CCP-based FRET, a cationic conjugated polymer (CCP)-based fluorescence resonance energy transfer (FRET); Restriction enzyme related methods include methylation-sensitive restriction enzyme-polymerase chain reaction and restriction enzyme cut analysis.

Sensitivity Analysis

As shown in Figure 3, leave-one-out sensitivity analysis was performed by removing a single study at a time under the random-effects model to assess the stability of the results. Our results showed that the ORs ranged from 24.04 (95% CI = 9.22–62.72) to 44.47 (95% CI = 16.61–119.08).

Figure 3

Figure 3

The sensitivity analysis of pooled OR for OPCML methylation and ovarian cancer risk under the random-effects model.

Publication Bias

Begg's and Egger's tests were performed to estimate the publication bias of the included studies. The funnel plot was approximately symmetric (Figure 4), indicating no publication bias in Begg's (P = 0.458) or Egger's tests (P = 0.261).

Figure 4

Figure 4

The funnel plot for assessment of publication bias in the meta-analysis.

Relationship Between OPCML Methylation and Clinicopathological Features

We divided the ovarian cancer patients into groups based on clinical features including stage, histological type, tumor differentiation, and age, and then further valuated the associations between other variables and OPCML methylation (Figure 5). The findings suggested that methylation was related to III/IV (OR = 4.20; 95% CI = 1.59–11.14) and poorly differentiated tumor (OR = 4.37; 95% CI = 1.14–16.78), while methylation was not related to serous histology (OR = 1.46; 95% CI = 0.88–2.41) or older age (OR = 1.42; 95% CI = 0.61–3.30).

Figure 5

Figure 5

Forest plots for OPCML methylation and clinicopathological features of ovarian cancer patients. (A) Forest plot for the relationship between OPCML methylation and III/IV stage. (B) Forest plot for the relationship between OPCML methylation and serous ovarian cancer. (C) Forest plot for the relationship between OPCML methylation and poorly differentiated tumor. (D) Forest plot for the relationship between OPCML methylation and older age.

Validation of the Results in Multiple Public Databases

Based on GSE146552 and GSE155760 datasets, we evaluated the DNA methylation profiling in ovarian cancer tissues and normal tissues (Figure 6A). We firstly analyzed the differential methylated sites between cancer tissues and normal tissues (Figure 6B). Then, we found that methylation levels of 55 differential methylated CpG sites were consistent in two datasets (Table 4). Among them, three CpG sites (cg16639665, cg23236270, cg15964611) were located in OPCML promoter regions (Figures 6C,D, all P < 0.05). In consistent with the results of meta-analysis, the methylation levels of three CpG sites were significantly higher in cancer tissues compared to normal tissues. Based on the TCGA ovarian cancer datasets, two OPCML-associated CpG sites (cg03923934, cg25853078, Figures 7A,B) were significantly associated with OPCML expression (rho = 0.261, P = 5.97 × 10−6; rho = 0.211, P = 2.75 × 10−4). However, we did not observe the significant associations for the clinicopathological parameters and overall survival time (Figures 7C–F).

Figure 6

Figure 6

The OPCML methylation status in ovarian cancer tissues and normal tissues. (A) The workflow of differential methylated sites in GSE146552 and GSE155760. (B) Volcano plots of differential methylated sites based on GSE146552 and GSE155760. (C) There are seven CpG sites in OPCML promoter regions. (D) The methylation status of seven CpG sites in GSE146552 and GSE155760. Among them, three differential methylated sites (cg16639665, cg23236270, cg15964611) were observed in both GSE146552 and GSE155760.

Table 4

IDchromchromStartchromEndGSE155760GSE146552Gene
logFCP-valuelogFCP-value
cg01875106chr11132287341132287343−0.63601.050 × 10−12−0.27401.500 × 10−4OPCML
cg11500790chr11132321480132321482−0.59801.400 × 10−7−0.34408.150 × 10−5OPCML
cg05806819chr11132432636132432638−0.01953.980 × 10−2−0.05233.340 × 10−2OPCML
cg03784083chr11132513894132513896−0.43004.610 × 10−4−0.21604.350 × 10−3OPCML
cg08080489chr11132526946132526948−0.57206.510 × 10−8−0.28001.240 × 10−4OPCML
cg17242596chr11132527295132527297−0.60003.930 × 10−8−0.39901.650 × 10−5OPCML
cg06967998chr11132527437132527439−0.62001.490 × 10−7−0.38802.850 × 10−5OPCML
cg10376306chr11132588617132588619−0.30607.050 × 10−4−0.15905.340 × 10−3OPCML
cg12156287chr11132662865132662867−0.40303.800 × 10−12−0.16002.330 × 10−2OPCML
cg18944144chr11132673498132673500−0.54701.730 × 10−6−0.26501.120 × 10−3OPCML
cg05035315chr11132691483132691485−0.61902.010 × 10−9−0.37307.990 × 10−6OPCML
cg13136241chr11132709827132709829−0.31901.080 × 10−3−0.17701.210 × 10−2OPCML
cg26824847chr11132728087132728089−0.51702.470 × 10−10−0.38101.270 × 10−6OPCML
cg19743254chr11132735813132735815−0.34101.090 × 10−3−0.19501.950 × 10−2OPCML
cg10966440chr11132745700132745702−0.65707.910 × 10−11−0.46708.810 × 10−7OPCML
cg23813681chr11132772322132772324−0.40601.880 × 10−4−0.18102.100 × 10−2OPCML
cg16813466chr11132811377132811379−0.40906.750 × 10−5−0.24401.190 × 10−3OPCML
cg07649562chr111328135031328135050.06624.840 × 10−20.04192.210 × 10−2OPCML
cg18710784chr111328135621328135640.25309.680 × 10−30.21001.020 × 10−3OPCML
cg15618210chr11132816802132816804−0.60109.760 × 10−8−0.30205.000 × 10−4OPCML
cg02993944chr11132831380132831382−0.32301.030 × 10−3−0.15201.250 × 10−2OPCML
cg16123848chr11132853100132853102−0.49609.470 × 10−6−0.33203.560 × 10−4OPCML
cg11747462chr11132869599132869601−0.26506.650 × 10−4−0.13202.100 × 10−2OPCML
cg03887438chr11132910641132910643−0.34103.090 × 10−13−0.15405.610 × 10−3OPCML
cg18413062chr11132931731132931733−0.42501.660 × 10−9−0.23505.720 × 10−4OPCML
cg13296579chr11132934961132934963−0.16303.160 × 10−2−0.08853.160 × 10−2OPCML, AP004782.1
cg20417723chr11132935119132935121−0.25704.250 × 10−3−0.14401.560 × 10−2OPCML, AP004782.1
cg26906285chr11132935462132935464−0.17201.820 × 10−2−0.15603.880 × 10−3OPCML, AP004782.1
cg05886537chr11132936427132936429−0.32502.030 × 10−3−0.17703.810 × 10−3OPCML, AP004782.1
cg08945802chr11132946454132946456−0.19801.710 × 10−2−0.08644.530 × 10−2OPCML
cg05807849chr11132946691132946693−0.41002.210 × 10−5−0.20907.820 × 10−3OPCML
cg11983942chr11132946933132946935−0.39108.510 × 10−8−0.19308.810 × 10−4OPCML
cg08357627chr11132949125132949127−0.51503.900 × 10−11−0.24803.460 × 10−4OPCML
cg11515095chr11132949436132949438−0.33201.500 × 10−3−0.18408.340 × 10−3OPCML
cg10335834chr11132949780132949782−0.28001.370 × 10−3−0.15701.030 × 10−2OPCML
cg15702349chr11132950103132950105−0.19704.550 × 10−3−0.09891.620 × 10−2OPCML
cg26237595chr11132955641132955643−0.29406.960 × 10−3−0.19106.130 × 10−3OPCML
cg01311955chr11132965249132965251−0.27801.140 × 10−4−0.26701.270 × 10−3OPCML
cg20972294chr11133005918133005920−0.37109.260 × 10−4−0.19401.420 × 10−2OPCML
cg12697833chr11133085432133085434−0.53709.690 × 10−8−0.28501.670 × 10−3OPCML
cg26311804chr11133155302133155304−0.37404.070 × 10−4−0.16806.610 × 10−3OPCML
cg04913291chr11133187955133187957−0.41509.420 × 10−6−0.28301.090 × 10−3OPCML
cg22209355chr11133208306133208308−0.36506.770 × 10−4−0.20506.410 × 10−4OPCML
cg05131064chr11133228365133228367−0.28202.300 × 10−3−0.13702.400 × 10−2OPCML
cg08872013chr11133229156133229158−0.23903.790 × 10−3−0.11002.460 × 10−2OPCML
cg20886183chr11133230390133230392−0.25304.390 × 10−3−0.19801.230 × 10−2OPCML, OPCML-IT1
cg01880176chr11133230949133230951−0.30707.980 × 10−5−0.15404.180 × 10−3OPCML, OPCML-IT1
cg06379368chr11133231146133231148−0.26501.750 × 10−3−0.16301.330 × 10−2OPCML, OPCML-IT1
cg13974487chr11133231880133231882−0.27201.790 × 10−3−0.17101.340 × 10−3OPCML, OPCML-IT1
cg18204273chr11133232926133232928−0.18502.260 × 10−2−0.15401.490 × 10−2OPCML, OPCML-IT1
cg01140660chr11133233199133233201−0.46601.190 × 10−4−0.21605.430 × 10−3OPCML, OPCML-IT1
cg09531376chr11133234072133234074−0.33701.850 × 10−3−0.15604.630 × 10−3OPCML, OPCML-IT1
cg16639665chr111334024961334024980.13803.450 × 10−20.11301.090 × 10−2OPCML
cg23236270chr111334024991334025010.16302.090 × 10−20.11504.540 × 10−3OPCML
cg15964611chr111334025441334025460.20702.630 × 10−30.10803.180 × 10−2OPCML

Consistently differential DNA methylated sites (P < 0.05 and the same direction of logFC in both datasets) between ovarian cancer tissues and normal tissues based on GEO DNA methylation datasets.

Bold values represent CpG sites in OPCML promoter region.

Figure 7

Figure 7

The methylation status of OPCML-associated CpG sites in the TCGA dataset. (A) There are only two CpG sites from HumanMethylation27k, and the location of two CpG sites in OPCML. (B) There are significantly positive associations between two CpG sites and OPCML expression. (C) The methylation levels were not significantly associated with clinical stage. (D) The methylation levels were not significantly associated with clinical grade. (E) No significant associations were observed in patient with different methylation levels using log-rank analysis. (F) No significant associations were observed in patient with different methylation levels using cox analysis.

Discussion

In the present study, we pooled eight studies together and compared the frequency of OPCML methylation in 476 ovarian cancers with 385 non-malignant tissues or blood samples. The results indicated that OPCML methylation was associated with ovarian cancer risk. The association was still significant in subgroups according to sample type, control type, methods, and sample size. Additionally, we observed significant associations between OPCML methylation and clinical stages and the grade of differentiation in ovarian cancer based on meta-analysis.

OPCML widely downregulated in many tumor types including brain tumors (Reed et al., 2007), non-small cell lung carcinoma (Tsou et al., 2007), bladder cancer (Duarte-Pereira et al., 2011), cholangiocarcinoma (Sriraksa et al., 2011), nasopharyngeal, esophageal, gastric, hepatocellular, colorectal, breast, cervical cancers, and ovarian cancer (Cui et al., 2008). Zanini et al. (2017) found that OPCML expression is associated with better progression free survival in HER2-positive ovarian cancer patients. Notably, exogenous recombinant OPCML protein inhibited ovarian cancer cell growth in vitro and in vivo (McKie et al., 2012; Wu and Sood, 2012). Recently, the advances in structure-function relationships of OPCML give rise to its potential as an anti-cancer therapy (Birtley et al., 2019). DNA methylation is the most common epigenetic modification and is widely reported for the transcriptional silencing of tumor suppressor genes in ovarian cancer. Previous studies have indicated that OPCML methylation is associated with the inactivation of OPCML, which has a tumor suppressor function in ovarian cancer (Simovic et al., 2019). However, studies with small sample sizes and different control types and methylation methods might produce inconsistent results. For example, in the study performed by Czekierdowski et al. (2006) OPCML methylation was observed in 20 out of 43 cases of ovarian cancer, but no methylation was found in four normal ovaries. No significant association was observed in this study. The control sample size in this study was small. Additionally, most previous studies reported the methylation frequency of OPCML in tumor and non-tumor tissues or serum, but the precise estimates for the associations were unclear. In the present study, we attempted to consolidate the available data and found that OPCML methylation was associated with an increased risk of ovarian cancer. The findings from the sensitivity analysis showed that our results are reliable and stable. Overall, given the TSG role and frequent methylation inactivation of OPCML in ovarian cancer, the restoration of OPCML expression could be a promising approach for ovarian cancer treatment.

In the subgroup analysis, the OR was 38.17 in China but 7.85 in Poland. The cause may be due to the limited studies in Poland. Therefore, more studies should be performed to confirm these observations in Europeans. In the serum subgroup, the OR value was highest (OR = 102.86) than that in the tissue subgroup (OR = 18.22). The results should be interpreted with caution because of the relatively small number of subjects included in the serum subgroup. Given that serum is a promising biomarker for non-invasive ovarian cancer diagnosis, further well-designed studies are warranted to explore the diagnostic performance of OPCML methylation in serum. Regarding the methods, the OR was 21.78 in the MSP group and 139.45 in the restriction enzyme-related analysis group. Technically, MSP amplify either methylated (M) or unmethylated (U) alleles after bisulfite conversion and evaluate methylation status in CpG regions (Herman et al., 1996). MSRE-PCR can recognize and degrade unmethylated DNA sequences whereas the methylated DNA sequences remain intact (Ramsahoye et al., 1997). Taken together, the results of subgroup analyses showed that OPCML methylation was significantly associated with an increased risk of ovarian cancer, regardless of the sample size, sample type, control type, and methods. Additionally, we found that the methylation status was associated with III/IV stage and poorly differentiated tumors. This finding suggests that OPCML is related to the malignant progression of ovarian cancer.

Based on publicly available datasets, we observed that the DNA methylation levels of cg16639665, cg23236270, and cg15964611 were consistently higher in cancer tissues than that in normal tissues, which validated some results of our meta-analysis. However, we did not confirm the difference of the OPCML methylation levels in different groups of clinical stage, grade, and so on. We found the positive association between DNA methylation level of cg25853078 in OPCML promoter region and OPCML expression. This positive association did not fit the hypothesis that DNA methylation levels in promoter region suppress the expression and contribute to ovarian cancer development. The DNA methylation levels in previous published papers were regions but not one CpG site. There are some CpG sites in OPCML promoter region, but we only assessed one CpG site in the promoter region based on the HumanMethylation27 platform. Therefore, other CpG sites were not calculated based on the TCGA data, which may have significant roles in the regulation of OPCML expression. Our study possessed several limitations. First, the total sample size in our study was relatively small and enrolled studies had low quality. Second, all the eligible studies were case-control designs and we could not explore the causation between OPCML methylation and ovarian cancer risk. Additionally, we used the incidence of OPCML methylation in the case and control groups to evaluate the pooled estimates because the adjusted ORs were not reported in the original studies. Further prospective studies with large samples are warranted to determine the role of OPCML in ovarian cancer risk and progression.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

Author contributions

YS, KX, and JW: conception and design. YS and JK: acquisition of data. HX, XW, YC, WL, JH, and DL: analysis and interpretation of data. YS, KX, and JW: manuscript drafting or revising. All authors contributed to the article and approved the submitted version.

Funding

This work was financially supported by the Nanjing Medical Science and Technique Development Foundation (JQX18009) and National Natural Science Foundation of China (81702569).

Acknowledgments

We are thankful to the publicly available datasets (TCGA, GEO).

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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.570898/full#supplementary-material

References

  • 1

    BeggC. B.MazumdarM. (1994). Operating characteristics of a rank correlation test for publication bias. Biometrics50, 1088–1101. 10.2307/2533446

  • 2

    BirtleyJ. R.AlomaryM.ZaniniE.AntonyJ.MabenZ.WeaverG. C.et al. (2019). Inactivating mutations and X-ray crystal structure of the tumor suppressor OPCML reveal cancer-associated functions. Nat. Commun. 10:3134. 10.1038/s41467-019-10966-8

  • 3

    ChenH.YeF.ZhangJ.LuW.ChengQ.XieX. (2007). Loss of OPCML expression and the correlation with CpG island methylation and LOH in ovarian serous carcinoma. Eur. J. Gynaecol. Oncol. 28, 464–467.

  • 4

    CuiY.YingY.van HasseltA.NgK. M.YuJ.ZhangQ.et al. (2008). OPCML is a broad tumor suppressor for multiple carcinomas and lymphomas with frequently epigenetic inactivation. PLoS ONE3:e2990. 10.1371/journal.pone.0002990

  • 5

    CzekierdowskiA.CzekierdowskaS.SzymanskiM.WielgosM.KaminskiP.KotarskiJ. (2006). Opioid-binding protein/cell adhesion molecule-like (OPCML) gene and promoter methylation status in women with ovarian cancer. Neuro Endocrinol. Lett. 27, 609–613.

  • 6

    Duarte-PereiraS.PaivaF.CostaV. L.Ramalho-CarvalhoJ.Savva-BordaloJ.RodriguesA.et al. (2011). Prognostic value of opioid binding protein/cell adhesion molecule-like promoter methylation in bladder carcinoma. Eur. J. Cancer. 47, 1106–1114. 10.1016/j.ejca.2010.12.025

  • 7

    EggerM.Davey SmithG.SchneiderM.MinderC. (1997). Bias in meta-analysis detected by a simple, graphical test. BMJ315, 629–634. 10.1136/bmj.315.7109.629

  • 8

    EhrlichM. (2002). DNA methylation in cancer: too much, but also too little. Oncogene21, 5400–5413. 10.1038/sj.onc.1205651

  • 9

    HermanJ. G.GraffJ. R.MyöhänenS.NelkinB. D.BaylinS. B. (1996). Methylation-specific PCR: a novel PCR assay for methylation status of CpG islands. Proc. Natl. Acad. Sci. U.S.A. 93, 9821–9826. 10.1073/pnas.93.18.9821

  • 10

    JonesP. A. (2012). Functions of DNA methylation: islands, start sites, gene bodies and beyond. Nat. Rev. Genet. 13, 484–492. 10.1038/nrg3230

  • 11

    LheureuxS.GourleyC.VergoteI.OzaA. M. (2019). Epithelial ovarian cancer. Lancet393, 1240–1253. 10.1016/S0140-6736(18)32552-2

  • 12

    LiuM.WangY.CaoX.ZhouF.QiaoH.LuX.et al. (2008). The relationship between the incidence of OPCML methylation and the development of epithelial ovarian cancer (in Chinese). China Oncol. 18, 950–952.

  • 13

    McKieA. B.VaughanS.ZaniniE.OkonI. S.LouisL.de SousaC.et al. (2012). The OPCML tumor suppressor functions as a cell surface repressor-adaptor, negatively regulating receptor tyrosine kinases in epithelial ovarian cancer. Cancer Discov. 2, 156–171. 10.1158/2159-8290.CD-11-0256

  • 14

    RamsahoyeB. H.BurnettA. K.TaylorC. (1997). Restriction endonuclease isoschizomers ItaI, BsoFI and Fsp4HI are characterised by differences in their sensitivities to CpG methylation. Nucl. Acids Res. 25, 3196–3198. 10.1093/nar/25.16.3196

  • 15

    RauscherG. H.KresovichJ. K.PoulinM.YanL.MaciasV.MahmoudA. M.et al. (2015). Exploring DNA methylation changes in promoter, intragenic, and intergenic regions as early and late events in breast cancer formation. BMC Cancer15:816. 10.1186/s12885-015-1777-9

  • 16

    ReedJ. E.DunnJ. R.du PlessisD. G.ShawE. J.ReevesP.GeeA. L.et al. (2007). Expression of cellular adhesion molecule 'OPCML' is down-regulated in gliomas and other brain tumours. Neuropathol. Appl. Neurobiol. 33, 77–85. 10.1111/j.1365-2990.2006.00786.x

  • 17

    SellarG. C.WattK. P.RabiaszG. J.StronachE. A.LiL.MillerE. P.et al. (2003). OPCML at 11q25 is epigenetically inactivated and has tumor-suppressor function in epithelial ovarian cancer. Nat. Genet. 34, 337–343. 10.1038/ng1183

  • 18

    SimovicI.Castaño-RodríguezN.KaakoushN. O. (2019). OPCML: a promising biomarker and therapeutic avenue. Trends Cancer5, 463–466. 10.1016/j.trecan.2019.06.002

  • 19

    SongL.JiaJ.PengX.XiaoW.LiY. (2017). The performance of the SEPT9 gene methylation assay and a comparison with other CRC screening tests: a meta-analysis. Sci. Rep. 7:3032. 10.1038/s41598-017-03321-8

  • 20

    SriraksaR.ZellerC.El-BahrawyM. A.DaiW.DaduangJ.JearanaikoonP.et al. (2011). CpG-island methylation study of liver fluke-related cholangiocarcinoma. Br. J. Cancer104, 1313–1318. 10.1038/bjc.2011.102

  • 21

    SungH.FerlayJ.SiegelR. L.LaversanneM.SoerjomataramI.JemalA.et al. (2021). Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 10.3322/caac.21660. [Epub ahead of print].

  • 22

    SzyfM. (2012). DNA methylation signatures for breast cancer classification and prognosis. Genome Med. 4:26. 10.1186/gm325

  • 23

    TsouJ. A.GallerJ. S.SiegmundK. D.LairdP. W.TurlaS.CozenW.et al. (2007). Identification of a panel of sensitive and specific DNA methylation markers for lung adenocarcinoma. Mol. Cancer6:70. 10.1186/1476-4598-6-70

  • 24

    VidalE.SayolsS.MoranS.Guillaumet-AdkinsA.SchroederM. P.RoyoR.et al. (2017). A DNA methylation map of human cancer at single base-pair resolution. Oncogene36, 5648–5657. 10.1038/onc.2017.176

  • 25

    WangB.YuL.LuoX.HuangL.LiQ. S.ShaoX. S.et al. (2017). Detection of OPCML methylation, a possible epigenetic marker, from free serum circulating DNA to improve the diagnosis of early-stage ovarian epithelial cancer. Oncol. Lett. 14, 217–223. 10.3892/ol.2017.6111

  • 26

    WangB.YuL.YangG. Z.LuoX.HuangL. (2015). Application of multiplex nested methylated specific PCR in early diagnosis of epithelial ovarian cancer. Asian Pa.c J. Cancer Prev. 16, 3003–3007. 10.7314/APJCP.2015.16.7.3003

  • 27

    WidschwendterM.JonesA.EvansI.ReiselD.DillnerJ.SundströmK.et al. (2018). Epigenome-based cancer risk prediction: rationale, opportunities and challenges. Nat. Rev. Clin. Oncol. 15, 292–309. 10.1038/nrclinonc.2018.30

  • 28

    WidschwendterM.ZikanM.WahlB.LempiäinenH.PaprotkaT.EvansI.et al. (2017). The potential of circulating tumor DNA methylation analysis for the early detection and management of ovarian cancer. Genome Med. 9:116. 10.1186/s13073-017-0500-7

  • 29

    WilsonA. S.PowerB. E.MolloyP. L. (2007). DNA hypomethylation and human diseases. Biochim. Biophys. Acta1775, 138–162. 10.1016/j.bbcan.2006.08.007

  • 30

    WuS. Y.SoodA. K. (2012). New roles opined for OPCML. Cancer Discov. 2, 115–116. 10.1158/2159-8290.CD-11-0356

  • 31

    XingB. L.LiT.TangZ. H.JiaoL.GeS. M.QiangX.et al. (2015). Cumulative methylation alternations of gene promoters and protein markers for diagnosis of epithelial ovarian cancer. Genet. Mol. Res. 14, 4532–4540. 10.4238/2015.May.4.11

  • 32

    ZaniniE.LouisL. S.AntonyJ.KaraliE.OkonI. S.McKieA. B.et al. (2017). The tumor-suppressor protein OPCML potentiates anti-EGFR- and anti-HER2-targeted therapy in HER2-positive ovarian and breast cancer. Mol Cancer Ther. 16, 2246–2256. 10.1158/1535-7163.MCT-17-0081

  • 33

    ZhangJ.XingB.SongJ.ZhangF.NieC.JiaoL.et al. (2014). Associated analysis of DNA methylation for cancer detection using CCP-based FRET technique. Anal. Chem. 86, 346–350. 10.1021/ac402720g

  • 34

    ZhangJ.YeF.ChenH. Z.YeD. F.LuW. G.XieX. (2006). Deletion of OPCML gene and promoter methylation in ovarian epithelial carcinoma (in Chinese). Zhongguo Yi Xue Ke Xue Yuan Xue Bao28, 173–177.

  • 35

    ZhouF.CaoX.MaM.LiuM.TaoG. (2011). Detection of OPCML methylation in tissues of ovarian cancer by MSRE-PCR. Chin. J. Clin. Lab. Sci. 29, 427–429.

  • 36

    ZhouF.MaM.TaoG.ChenX.XieW.WangY.et al. (2014a). Detection of circulating methylated opioid binding protein/cell adhesion molecule-like gene as a biomarker for ovarian carcinoma. Clin. Lab. 60, 759–765. 10.7754/Clin.Lab.2013.130446

  • 37

    ZhouF.TaoG.ChenX.XieW.LiuM.CaoX. (2014b). Methylation of OPCML promoter in ovarian cancer tissues predicts poor patient survival. Clin. Chem. Lab. Med. 52, 735–742. 10.1515/cclm-2013-0736

Summary

Keywords

OPCML, methylation, ovarian cancer, risk, progression

Citation

Shao Y, Kong J, Xu H, Wu X, Cao Y, Li W, Han J, Li D, Xie K and Wu J (2021) OPCML Methylation and the Risk of Ovarian Cancer: A Meta and Bioinformatics Analysis. Front. Cell Dev. Biol. 9:570898. doi: 10.3389/fcell.2021.570898

Received

01 July 2020

Accepted

08 February 2021

Published

11 March 2021

Volume

9 - 2021

Edited by

Xiao Zhu, Guangdong Medical University, China

Reviewed by

Abhijit Shukla, Memorial Sloan Kettering Cancer Center, United States; Mingzhi Han, University of Bergen, Norway

Updates

Copyright

*Correspondence: Jiangping Wu Kaipeng Xie

This article was submitted to Epigenomics and Epigenetics, a section of the journal Frontiers in Cell and Developmental Biology

†These authors have contributed equally to this work

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics