Epithelial ovarian cancer (EOC) is a highly lethal malignancy due to a lack of early detection approaches coupled with poor outcomes for patients with clinically advanced disease. Cancer-testis (CT) or cancer-germline genes encode antigens known to generate spontaneous anti-tumor immunity in cancer patients. CT45 genes are a recently discovered 6-member family of X-linked CT genes with oncogenic function. Here, we determined CT45 expression in EOC and fully defined its epigenetic regulation by DNA methylation. CT45 was silent and hypermethylated in normal control tissues, but a large subset of EOC samples showed increased CT45 expression in conjunction with promoter DNA hypomethylation. In contrast, copy number status did not correlate with CT45 expression in the TCGA database for EOC. CT45 promoter methylation inversely correlated with both CT45 mRNA and protein expression, the latter determined using IHC staining of an EOC TMA. CT45 expression was increased and CT45 promoter methylation was decreased in late-stage and high-grade EOC, and both measures were associated with poor survival. CT45 hypomethylation was directly associated with LINE-1 hypomethylation, and CT45 was frequently co-expressed with other CT antigen genes in EOC. Decitabine treatment induced CT45 mRNA and protein expression in EOC cells, and promoter transgene analyses indicated that DNA methylation directly represses CT45 promoter activity. These data verify CT45 expression and promoter hypomethylation as possible prognostic biomarkers, and suggest CT45 as an immunological or therapeutic target in EOC. Treatment with decitabine or other epigenetic modulators could provide a means for more effective immunological targeting of CT45.
Epithelial ovarian cancer (EOC) is a highly lethal malignancy due to a lack of early detection approaches coupled with poor outcomes for patients with clinically advanced disease. Cancer-testis (CT) or cancer-germline genes encode antigens known to generate spontaneous anti-tumor immunity in cancerpatients. CT45 genes are a recently discovered 6-member family of X-linked CT genes with oncogenic function. Here, we determined CT45 expression in EOC and fully defined its epigenetic regulation by DNA methylation. CT45 was silent and hypermethylated in normal control tissues, but a large subset of EOC samples showed increased CT45 expression in conjunction with promoter DNA hypomethylation. In contrast, copy number status did not correlate with CT45 expression in the TCGA database for EOC. CT45 promoter methylation inversely correlated with both CT45 mRNA and protein expression, the latter determined using IHC staining of an EOC TMA. CT45 expression was increased and CT45 promoter methylation was decreased in late-stage and high-grade EOC, and both measures were associated with poor survival. CT45 hypomethylation was directly associated with LINE-1 hypomethylation, and CT45 was frequently co-expressed with other CT antigen genes in EOC. Decitabine treatment induced CT45 mRNA and protein expression in EOC cells, and promoter transgene analyses indicated that DNA methylation directly represses CT45 promoter activity. These data verify CT45 expression and promoter hypomethylation as possible prognostic biomarkers, and suggest CT45 as an immunological or therapeutic target in EOC. Treatment with decitabine or other epigenetic modulators could provide a means for more effective immunological targeting of CT45.
Entities:
Keywords:
CNA, copy number alteration; CT antigen genes, cancer-testis or cancer-germline antigen genes; CT45; DAC, decitabine, 5-Aza-2′-deoxycytidine; DFS, disease-free survival; DNA methylation; DNMT, DNA methyltransferase; EOC, epithelial ovarian cancer; FTE, normal fallopian tube epithelia; HGSOC, high-grade serous ovarian cancer; IHC, immunohistochemistry; NO, bulk normal ovary; OS, overall survival; OSE, normal ovary surface epithelia; RLM-RACE, 5′ RNA ligase-mediated rapid amplification of cDNA ends; RNA-seq, RNA sequencing; TCGA, The Cancer Genome Atlas; TMA, tissue microarray; TSS, transcription start site; cancer germline genes; cancer testis antigen genes; decitabine; epithelial ovarian cancer; tumor antigens
copy number alterationcancer-testis or cancer-germline antigen genesdecitabine, 5-Aza-2′-deoxycytidinedisease-free survivalDNA methyltransferaseepithelial ovarian cancernormal fallopian tube epitheliahigh-grade serous ovarian cancerimmunohistochemistrybulk normal ovaryoverall survivalnormal ovary surface epithelia5′ RNA ligase-mediated rapid amplification of cDNA endsRNA sequencingThe Cancer Genome Atlastissue microarraytranscription start site
Introduction
Ovarian cancer is the seventh most common cancer among women worldwide, and the fifth leading cause of cancer death in US women, with an estimated 14,180 deaths in 2015. Ninety percent of ovarian cancer is epithelial in origin (EOC). The 5 year survival rate for women diagnosed with localized EOC exceeds 90%; however, more than 60% of cases are diagnosed with advanced disease, and the 5 year survival rate for these patients is less than 30%. Current standard of care therapy for EOC is cytoreductive surgery followed by combination platinum-taxane chemotherapy. Disease recurrence occurs in greater than 70% of patients initially responsive to this therapy, and there is no durable cure for recurrent ovarian cancer. Achieving an improved understanding of the molecular pathogenesis of EOC is required to develop new biomarkers for early detection as well as novel therapeutic approaches.Cancer testis (CT) antigen genes, also known as cancer germline genes, are normally expressed only in germ cells, but are aberrantly expressed in a wide range of humancancers. CT genes frequently encode antigens that can elicit spontaneous immune responses in cancerpatients, and this, coupled with their highly restricted expression in normal somatic tissues, makes them ideal immunotherapy targets. CT antigen genes are often regulated at the transcriptional level by promoter DNA methylation and other epigenetic mechanisms, including global DNA methylation status, histone acetylation, and nucleosome occupancy.CT45 is a recently discovered sub-family of CT genes composed of 6 highly related paralogs: CT45A1 (GenBank ID: 541466), CT45A2 (GenBank ID: 728911), CT45A3 (GenBank ID: 441519), CT45A4 (GenBank ID: 441520), CT45A5 (GenBank ID: 441521), and CT45A6 (GenBank ID: 541465). CT45 genes are located in a tandem arrangement at chromosome Xq26.3 and encode highly similar mRNAs and proteins, with 98% conservation at the amino acid level, suggesting evolutionary emergence via gene duplications. CT45 is localized in the nucleus, and 2 nuclear localization signals have been defined. In normal tissues, CT45 mRNA expression was detected with a strong signal in testis and weak signal in lung, heart, spinal cord, and brain, while the CT45 protein was detected only in primary spermatocytes and gonocytes. Unlike its highly restricted expression in normal tissues, CT45 is expressed widely in humancancer. In some tumor types, CT45 is expressed in a low proportion of cases, e.g., myeloma (16%), lung (13%), colon (11%), and breast (5–10%), while, in other tumor types, CT45 exhibits widespread expression, e.g., seminoma (>90%), Hodgkin's lymphoma (>50%) , EOC (37%), and diffuse large B-cell lymphoma (22%). In addition to its mere expression in cancer, CT45 expression has also been correlated with more aggressive disease. In Hodgkin's lymphoma, CT45 expression correlated with more aggressive histological subtypes, disease symptoms, and advanced stages. In multiple myeloma, CT45 expression correlated with poor survival. In breast cancer, CT45 expression correlated with high grade, and ER negative, HER2 negative status, suggesting it could be used as a potential immunotherapeutic target in triple negative breast cancer.Widespread expression in specific humancancers coupled with increased expression in more aggressive disease states have led to speculation that CT45 may play direct roles in oncogenesis. In support of this hypothesis, 2 recent studies of CT45 provide initial evidence for an oncogenic role for CT45 proteins. Koop et al. demonstrated that CT45 knockdown in Hodgkin's lymphoma and multiple myeloma cell lines reduced cell adhesion, migration, and invasion. Shang et al. similarly found that CT45A1 overexpression increased migration and invasion, promoted epithelial to mesenchymal transition (EMT), and increased xenograft growth and metastasis, using the MCF7 breast cancer cell line as a model. In both studies, CT45 was not observed to play a role in cell growth or cell cycle progression, but only in phenotypic characteristics related to metastatic tumor growth.The mechanism leading to CT45 expression in cancer has not been well defined. One study reported that treatment with decitabine, a DNA hypomethylating agent, induced CT45 expression in HeLa cells and stimulated peripheral blood lymphocytes. However, decitabine has both on-target and off-target effects, and the promoter methylation status of CT45 genes has not been directly determined, either in samples expressing these genes or in those in which it is silent. Furthermore, demonstration that DNA methylation directly impacts CT45 promoter activity has not been provided. In the current study, we conducted a thorough investigation of the expression of CT45 genes in EOC and normal control tissues, and the role of DNA methylation (both promoter-specific and global methylation) in its regulation. In addition, we analyzed the relationship between CT45 expression, DNA methylation, and EOC disease progression and patient survival.
Results
CT45 mRNA expression in EOC
Due to the high sequence conservation of the 6 CT45 paralogs, we were unable to accurately measure individual CT45 gene expression by RT-PCR (data not shown). Therefore, we used pan-CT45 amplification primers to measure CT45 expression by reverse transcriptase quantitative PCR (RT-qPCR) (Table S1). We first determined CT45 expression in a large set (>100) of EOC samples as well as bulk normal ovary (NO) samples as controls (EOC clinical information shown in Table S2). As shown in , CT45 was not expressed in NO but was activated in a significant proportion of EOC (25%). The majority of EOC, including those in our sample collection, have serous histology (Table S2). We observed that serous, but not non-serous EOC, expressed significantly higher CT45 as compared to NO (; Fig. S1). We examined the influence of disease stage and grade on CT45 expression and found that CT45 is activated only in late stage and high-grade disease (). In agreement, CT45 protein levels, as determined by immunohistochemistry (IHC) staining of EOC tissue microarrays (TMA), indicated that CT45 protein expression was increased in late stage, high-grade disease (Table S3). Increased CT45 expression correlated with reduced overall survival (OS) in EOC (). A similar trend was observed for disease-free survival (DFS) but the difference was not statistically significant, likely due to the small set number of patients available for this analysis ().
Figure 1.
CT45 mRNA expression in EOC and its relation to clinicopathology. (A-D). CT45 mRNA expression relative to 18s rRNA expression is shown for: (A) Bulk normal ovary (NO) and EOC, (B) EOC separated into non-serous and serous histotypes, (C) EOC separated by disease stage, and (D) EOC separated by pathological grade. Gene expression was determined by RT-qPCR. The two-tailed Mann-Whitney test p-value is indicated (*: P<0.05). For A-B, the proportion of EOC and serous samples with elevated CT45 mRNA expression, respectively, is indicated. (E) Overall survival (OS) and (F) disease-free survival (DFS) of EOC patients as a function of CT45 expression. Patients were dichotomized at the median CT45 mRNA expression value. Logrank P-values are indicated.
CT45 mRNA expression in EOC and its relation to clinicopathology. (A-D). CT45 mRNA expression relative to 18s rRNA expression is shown for: (A) Bulk normal ovary (NO) and EOC, (B) EOC separated into non-serous and serous histotypes, (C) EOC separated by disease stage, and (D) EOC separated by pathological grade. Gene expression was determined by RT-qPCR. The two-tailed Mann-Whitney test p-value is indicated (*: P<0.05). For A-B, the proportion of EOC and serous samples with elevated CT45 mRNA expression, respectively, is indicated. (E) Overall survival (OS) and (F) disease-free survival (DFS) of EOC patients as a function of CT45 expression. Patients were dichotomized at the median CT45 mRNA expression value. Logrank P-values are indicated.
CT45 copy number alterations do not correlate with CT45 mRNA expression in EOC
Copy number alterations are very common in EOC, and have the capacity to deregulate gene expression. As mentioned earlier, CT45 genes are located in a tight arrangement on Xq26.3, and coordinated chromosomal gains or amplifications at this region could potentially account for CT45 activation in EOC. To investigate this possibility, we retrieved The Cancer Genome Atlas (TCGA) data, including RNA sequencing (RNA-seq) and copy number alteration (CNA) data for high-grade serous ovarian cancer (HGSOC), the most common EOC subtype. Of 265 tumors with RNA-seq data (microarray data was sparse for CT45), 37 (14%) showed activation of at least one CT45 gene (). This was a lower frequency than we observed by RT-qPCR in our data set, potentially due to the increased sensitivity of RT-qPCR or to our measurement of Pan-CT45 expression (). As expected, CT45 family members showed overlapping expression (). Among 579 tumors with CNA data, 70% of cases showed a CT45 alteration, and there was complete overlap among the 6 CT45 genes (Fig. S2). The CNAs included roughly equivalent amounts of gains and deletions (). Notably, CT45 expression did not correlate with copy number status (; representative data shown for CT45A1). Thus, alternative mechanisms are likely to account for the expression pattern of CT45 observed in EOC.
Figure 2.
CT45 mRNA expression and copy number alterations (CNA) in EOC. Data were obtained from the TCGA and analyzed using cBioPortal as described in Materials and Methods. (A)
CT45 family mRNA expression in 265 HGSOC EOC tumors, as determined by RNA-seq. Upregulation, downregulation, and the percentage of cases altered are indicated. (B) Synopsis of CT45 CNA in 579 HGSOC samples. (C)
CT45A1 expression plotted as a function of CT45A1 copy number status. Mann-Whitney t-tests indicated no significant differences in any pairwise group analysis; similar results were obtained for other CT45 family members (data not shown).
CT45 mRNA expression and copy number alterations (CNA) in EOC. Data were obtained from the TCGA and analyzed using cBioPortal as described in Materials and Methods. (A)
CT45 family mRNA expression in 265 HGSOC EOC tumors, as determined by RNA-seq. Upregulation, downregulation, and the percentage of cases altered are indicated. (B) Synopsis of CT45 CNA in 579 HGSOC samples. (C)
CT45A1 expression plotted as a function of CT45A1 copy number status. Mann-Whitney t-tests indicated no significant differences in any pairwise group analysis; similar results were obtained for other CT45 family members (data not shown).
5′ end CT45 transcript mapping in normal testis and EOC
The 5′ end and promoter regions of CT45 genes do not contain classical CpG islands; however, the predicted CT45 transcriptional start sites (TSS) are flanked by several CpG dinucleotides, suggesting that DNA methylation could play a transcriptional regulatory role for these genes (Fig. S3). As a prelude to understanding the potential role of DNA methylation in CT45 regulation, we first mapped the TSS of CT45 genes. Specifically, we used 5′ RNA ligase-mediated rapid amplification of cDNA ends (RLM-RACE) to determine the TSS in both normal human testis and an admixture of 2 EOC samples with high CT45 expression. RLM-RACE primer sequences are shown in Table S1. We observed several alternative TSSs for CT45, all but one of which was downstream of the NCBI-predicted TSS (Fig. S3). The variable nature of the CT45 TSS is consistent with the absence of a TATA box in its promoter. The TSS detected in normal testis and EOC overlapped, suggesting relative conservation in normal tissues and cancer. Unexpectedly, RLM-RACE additionally detected fusion sequences in EOC, in which the CT45 sequence contained an insert derived from a nearby locus on chromosome X (data not shown). A previous study reported MLL-CT45A2 fusion transcripts in pediatric biphenotypic acute leukemia. These observations suggest that the CT45 locus could be a hotspot for cancer-related translocations.
CT45 promoter hypomethylation in EOC
Based on the NCBI-predicted TSS for CT45, as well as that determined by RLM-RACE, we designed bisulfite clonal sequencing primers covering 2 relevant regions of the CT45 promoter region (Table S1; ). We analyzed 5 samples by bisulfite sequencing, including one NO, 2 EOC with high CT45 expression, and 2 EOC with low CT45 expression. We observed highest CT45 methylation of both regions in NO, some hypomethylation in the CT45 negative EOC samples, and the highest degree of hypomethylation in the CT45 positive EOC samples (), suggesting an indirect association between CT45 expression and promoter methylation. To confirm the relationship between promoter hypomethylation and CT45 expression in EOC, we designed high-throughput bisulfite pyrosequencing assays to interrogate 4 CpG sites showing correlation between hypomethylation and expression in bisulfite sequencing (Table S1; ). In these experiments, in addition to using bulk normal ovary (NO) as a control, we also measured CT45 methylation in primary ovarian surface epithelia (OSE) and fallopian tube epithelia (FTE), 2 presumed tissues of origin of EOC. As shown in , OSE and FTE had high CT45 promoter methylation, NO also had high levels of CT45 methylation (however a small degree of hypomethylation was apparent in Region 2), and EOC samples showed heterogeneous methylation, with roughly half (Region 1) or 2-thirds (Region 2) of the samples showing hypomethylation as compared to control tissues. Reduced methylation at Region 2 in bulk NO, as compared to OSE and FTE, may be due to the presence of other cell types, including germ cells, in this tissue (). Methylation at the 2 different CT45 promoter regions was highly correlated in EOC ().
Figure 3.
Sodium bisulfite clonal sequencing of the CT45 promoter. (A) Upper: diagram of the CT45 promoter region, indicating CpG sites (black hash marks), predicted or mapped TSS (right broken arrows, see Fig. S3), regions analyzed by pyrosequencing (rectangles; see subsequent figures), and a CpG site differentially present in CT45 family members (inverted triangle). Lower: bisulfite clonal sequencing results for 2 regions of the CT45 promoter. Filled and open circles indicate methylated and unmethylated CpG sites, and each row represents one sequenced allele. Data are shown for a NO sample (NO1), 2 EOC samples with low CT45 expression (EOC99 and EOC83), and 2 EOC samples with high CT45 expression (EOC100 and EOC89). (B) Graphical summary of the data presented in A.
Figure 4.
Sodium bisulfite pyrosequencing of the CT45 promoter. The methylation status of 2 different regions of the CT45 promoter (indicated by rectangles in ), were analyzed by sodium bisulfite pyrosequencing in an expanded set of clinical samples including: OSE: normal ovary surface epithelia, FTE: normal fallopian tube epithelia, NO: normal bulk ovary tissue, and EOC: epithelial ovarian cancer. (A) Pyrosequencing data for Region 1. (B) Pyrosequencing data for Region 2. The proportion of EOC samples with CT45 promoter hypomethylation as compared to OSE and FTE are indicated. The two-tailed Mann-Whitney test P-value is indicated (*: P <0.05; ***: P<0.001; ****: P<0.0001). (C) Correlation between CT45 Region 1 and Region 2 methylation in EOC samples. Spearman test results are shown.
Sodium bisulfite clonal sequencing of the CT45 promoter. (A) Upper: diagram of the CT45 promoter region, indicating CpG sites (black hash marks), predicted or mapped TSS (right broken arrows, see Fig. S3), regions analyzed by pyrosequencing (rectangles; see subsequent figures), and a CpG site differentially present in CT45 family members (inverted triangle). Lower: bisulfite clonal sequencing results for 2 regions of the CT45 promoter. Filled and open circles indicate methylated and unmethylated CpG sites, and each row represents one sequenced allele. Data are shown for a NO sample (NO1), 2 EOC samples with low CT45 expression (EOC99 and EOC83), and 2 EOC samples with high CT45 expression (EOC100 and EOC89). (B) Graphical summary of the data presented in A.Sodium bisulfite pyrosequencing of the CT45 promoter. The methylation status of 2 different regions of the CT45 promoter (indicated by rectangles in ), were analyzed by sodium bisulfite pyrosequencing in an expanded set of clinical samples including: OSE: normal ovary surface epithelia, FTE: normal fallopian tube epithelia, NO: normal bulk ovary tissue, and EOC: epithelial ovarian cancer. (A) Pyrosequencing data for Region 1. (B) Pyrosequencing data for Region 2. The proportion of EOC samples with CT45 promoter hypomethylation as compared to OSE and FTE are indicated. The two-tailed Mann-Whitney test P-value is indicated (*: P <0.05; ***: P<0.001; ****: P<0.0001). (C) Correlation between CT45 Region 1 and Region 2 methylation in EOC samples. Spearman test results are shown.
CT45 promoter hypomethylation correlates with CT45 mRNA and protein expression in EOC
To determine if CT45 hypomethylation correlates with CT45 expression, we first compared CT45 methylation as determined by pyrosequencing to CT45 mRNA expression as determined by RT-qPCR. For both pyrosequencing regions, we observed a significant inverse correlation with CT45 mRNA expression in EOC (). Second, by comparison to a previously published TMA IHC data set, we found that CT45 promoter methylation was significantly lower in EOC samples positive for CT45 protein expression as compared to those negative for CT45 (). Representative examples of negative, weak, and strong CT45 IHC-staining HGSOC tumors are shown in . The staining pattern observed was nuclear ().
Figure 5.
CT45 expression and promoter methylation in EOC. (A-B)
CT45 mRNA expression was measured by RT-qPCR and was normalized to 18s rRNA expression. CT45 Region 1 (A) and Region 2 (B) methylation were measured by pyrosequencing. Samples with undetectable CT45 mRNA are plotted on the x-axis. Spearman test results are shown. (C-D) CT45 protein expression was measured by IHC using an EOC TMA as described in Materials and Methods. EOC samples with weak (W), moderate (M), or strong (S) CT45 IHC signal (see Table S1) were classified as CT45-positive, while samples with no CT45 IHC signal were classified as CT45-negative. CT45 Region 1 (C) and Region 2 (D) methylation were measured by pyrosequencing. The two-tailed Mann-Whitney test p-value is indicated (*: P < 0.05). (E) CT45 IHC TMA staining of HGSOC samples, showing representative images of negative, weak, and strong staining. Upper and lower rows show low and high magnification images of 0.6 mm TMA cores, respectively.
CT45 expression and promoter methylation in EOC. (A-B)
CT45 mRNA expression was measured by RT-qPCR and was normalized to 18s rRNA expression. CT45 Region 1 (A) and Region 2 (B) methylation were measured by pyrosequencing. Samples with undetectable CT45 mRNA are plotted on the x-axis. Spearman test results are shown. (C-D) CT45 protein expression was measured by IHC using an EOC TMA as described in Materials and Methods. EOC samples with weak (W), moderate (M), or strong (S) CT45 IHC signal (see Table S1) were classified as CT45-positive, while samples with no CT45 IHC signal were classified as CT45-negative. CT45 Region 1 (C) and Region 2 (D) methylation were measured by pyrosequencing. The two-tailed Mann-Whitney test p-value is indicated (*: P < 0.05). (E) CT45 IHC TMA staining of HGSOC samples, showing representative images of negative, weak, and strong staining. Upper and lower rows show low and high magnification images of 0.6 mm TMA cores, respectively.
CT45 promoter methylation is decreased during disease progression and is associated with reduced survival in EOC patients
We next compared CT45 promoter methylation levels, at both pyrosequencing regions, with clinical stage and histological grade in EOC. While both pyrosequencing regions showed hypomethylation in advanced stage disease, the trend in CT45 promoter Region 1 also met statistical significance (). Similar results were observed for tumor grade (). Thus, similar to CT45 expression (), CT45 promoter hypomethylation correlates with disease progression. With regards to survival, OS was reduced in patients with CT45 promoter hypomethylation, and this difference was highly significant for Region 1 (). Although fewer patients were available for DFS analysis, we observed a trend toward reduced DFS with hypomethylation of CT45 Region 1, and a statistically significant difference for CT45 Region 2 ().
Figure 6.
CT45 promoter methylation and EOC disease progression. (A-B)
CT45 promoter Region 1 (A) and Region 2 (B) methylation was determined by pyrosequencing and is plotted as a function of EOC clinical stage. (C-D)
CT45 promoter Region 1 (C) and Region 2 (D) methylation was determined by pyrosequencing and is plotted as a function of EOC histological grade. The two-tailed Mann-Whitney test p-value is indicated (*: P < 0.05).
Figure 7.
CT45 promoter methylation and EOC patient survival. (A-B)
CT45 promoter Region 1 (A) and Region 2 (B) methylation was determined by pyrosequencing. Patients were dichotomized at the median methylation value and OS was compared using Kaplan Meier analysis and Logrank test. (C-D) Same as A-B except DFS was compared.
CT45 promoter methylation and EOC disease progression. (A-B)
CT45 promoter Region 1 (A) and Region 2 (B) methylation was determined by pyrosequencing and is plotted as a function of EOC clinical stage. (C-D)
CT45 promoter Region 1 (C) and Region 2 (D) methylation was determined by pyrosequencing and is plotted as a function of EOC histological grade. The two-tailed Mann-Whitney test p-value is indicated (*: P < 0.05).CT45 promoter methylation and EOC patient survival. (A-B)
CT45 promoter Region 1 (A) and Region 2 (B) methylation was determined by pyrosequencing. Patients were dichotomized at the median methylation value and OS was compared using Kaplan Meier analysis and Logrank test. (C-D) Same as A-B except DFS was compared.
CT45 hypomethylation and expression are associated with LINE-1 hypomethylation, and CT45 is co-expressed with other CT antigen genes, in EOC
Prior data from the study of other CT antigen genes indicates connections between global loss of DNA methylation in cancer, promoter demethylation, and gene expression. To address this critical question in the context of CT45 in EOC, we measured LINE-1 methylation as a surrogate for global methylation status, and performed correlation analysis with CT45 promoter methylation as determined by pyrosequencing. As shown in , there was significant direct correlation between CT45 promoter methylation and LINE-1 methylation in EOC. Consistently, there was significant indirect correlation between CT45 expression and LINE-1 methylation (). These data suggest that transcriptional regulation of CT45 genes in EOC may in part reflect the global DNA methylation status of these tumors. If true, we would expect CT45 expression to correlate with the expression of other CT antigen genes in EOC. Consistent with this idea, analysis of Oncomine datasets confirmed that CT45A3 is co-expressed with other CT antigen genes in EOC, including several known to be regulated by DNA methylation, including MAGE, GAGE, and XAGE genes ().
Figure 8.
CT45 methylation and expression, LINE-1 methylation, and CT gene co-expression in EOC. (A-B)
CT45 promoter Region 1 (A) and Region 2 (B) methylation as compared to LINE-1 methylation in EOC. All methylation data were determined by pyrosequencing. (C)
CT45 mRNA expression as compared to LINE-1 methylation in EOC. CT45 mRNA expression was measured by RT-qPCR and was normalized to 18s rRNA expression. Samples with undetectable CT45 mRNA levels were plotted on the x-axis. (D)
Oncomine analysis of the genes most closely correlated with CT45A3 in the Bittner ovarian cancer dataset (see Materials and Methods for details). Correlation coefficient and gene name are shown at left. All genes besides C18orf2 are known, or likely (GTSF1), CT genes (see CTDatabase, http://www.cta.lncc.br/).
CT45 methylation and expression, LINE-1 methylation, and CT gene co-expression in EOC. (A-B)
CT45 promoter Region 1 (A) and Region 2 (B) methylation as compared to LINE-1 methylation in EOC. All methylation data were determined by pyrosequencing. (C)
CT45 mRNA expression as compared to LINE-1 methylation in EOC. CT45 mRNA expression was measured by RT-qPCR and was normalized to 18s rRNA expression. Samples with undetectable CT45 mRNA levels were plotted on the x-axis. (D)
Oncomine analysis of the genes most closely correlated with CT45A3 in the Bittner ovarian cancer dataset (see Materials and Methods for details). Correlation coefficient and gene name are shown at left. All genes besides C18orf2 are known, or likely (GTSF1), CT genes (see CTDatabase, http://www.cta.lncc.br/).
Decitabine treatment induces CT45 mRNA and protein expression in EOC cell lines
The primary tumor data presented above establishes a correlation between CT45 promoter hypomethylation and CT45 expression, but does not demonstrate causation. To establish a mechanistic connection, we first treated EOC cell lines with varying levels of baseline CT45 expression with the DNA methyltransferase (DNMT) inhibitor decitabine (DAC), a cytosine analog that specifically inhibits DNA methylation by covalently trapping DNMTs to DNA. Following DAC treatment, the CT45 promoter was hypomethylated in all treated cell lines (). Notably, CT45 expression was robustly induced in each of the low expressing cell lines (i.e., IOSE121, OVCAR3, SNU119), with less or no induction in 2 lines with higher basal CT45 expression (A2780, OVSAHO) (). Importantly, CT45 protein was also induced, albeit to varying levels, in EOC cell lines treated with DAC ().
Figure 9.
Decitabine-mediated induction of CT45 expression in EOC cells. The indicated cell lines were treated with 1µM decitabine (DAC) for 5 days as described in Materials and Methods. (A)
CT45 mRNA expression and CT45 promoter methylation after vehicle (PBS) or DAC treatment. CT45 mRNA expression was determined by RT-qPCR and CT45 promoter methylation was determined by pyrosequencing (average of Region 1 + Region 2 methylation). Data represent mean +/– SD. The two-tailed paired t-test p-value is indicated (*:P < 0.05; **: P < 0.01; ****: P < 0.0001; NS: not significant). (B) Western blot analysis of CT45 protein expression after vehicle (PBS, –) or DAC (+) treatment. β-actin and FastGreen staining are shown as loading controls.
Decitabine-mediated induction of CT45 expression in EOC cells. The indicated cell lines were treated with 1µM decitabine (DAC) for 5 days as described in Materials and Methods. (A)
CT45 mRNA expression and CT45 promoter methylation after vehicle (PBS) or DAC treatment. CT45 mRNA expression was determined by RT-qPCR and CT45 promoter methylation was determined by pyrosequencing (average of Region 1 + Region 2 methylation). Data represent mean +/– SD. The two-tailed paired t-test p-value is indicated (*:P < 0.05; **: P < 0.01; ****: P < 0.0001; NS: not significant). (B) Western blot analysis of CT45 protein expression after vehicle (PBS, –) or DAC (+) treatment. β-actin and FastGreen staining are shown as loading controls.
DNA methylation directly represses CT45 promoter activity
To determine whether the inhibitory action of DNA methylation on CT45 expression is direct or indirect, we used a promoter-luciferase reporter approach. We cloned 2 overlapping regions of the CT45 promoter into the CpG-free vector pCpGL-basic, which allowed specific methylation of CT45 promoter inserts. A schematic of the CT45A2 promoter (the gene chosen for analysis) along with the position of different methyltransferase acceptor sites is shown in Figure 10A. We tested the activity of both CT45 promoter constructs, in the absence or presence of methylation, using 2 different EOC cell lines, OVCAR3 and A2780. The constructs were either mock-methylated or methylated with one of 3 different CpG methylases that targeted a different number of CpG sites in each construct (). As shown in Figs. 10B-C, highly similar results were obtained in each EOC cell line, with repression of each promoter construct by DNA methylation. Interestingly, construct 2, but not construct 1, showed maximal repression with either HhaI or HpaII methylation, despite the fact that less CpG sites are methylated by these enzymes as compared to M.SssI (Figs. 10B-C). The in vivo relevance of this observation remains to be determined. Regardless, our data establish DNA methylation as a direct repressor of CT45 promoter activity in EOC cells.
Figure 10.
Effect of DNA methylation on CT45 promoter activity. (A) Diagram of the CT45A2 promoter region and 2 constructs used for promoter activity assay. Black hash marks represent CpG sites. NCBI-predicted and RLM-RACE mapped TSS are indicated by the right broken arrows as in . The two regions analyzed by pyrosequencing are indicated by black lines. The different methyltransferase acceptor sites are labeled in the legend. (B-C) Luciferase activity data for CT45 promoter constructs in (B) OVCAR3 and (C) A2780 EOC cell lines. Constructs were either mock-methylated or methylated with HhaI, HpaII, or M.SssI prior to transfection as described in Materials and Methods. CT45 firefly luciferase/renilla luciferase (control) enzymatic activity data are shown. Data represent mean ± SD.
Effect of DNA methylation on CT45 promoter activity. (A) Diagram of the CT45A2 promoter region and 2 constructs used for promoter activity assay. Black hash marks represent CpG sites. NCBI-predicted and RLM-RACE mapped TSS are indicated by the right broken arrows as in . The two regions analyzed by pyrosequencing are indicated by black lines. The different methyltransferase acceptor sites are labeled in the legend. (B-C) Luciferase activity data for CT45 promoter constructs in (B) OVCAR3 and (C) A2780 EOC cell lines. Constructs were either mock-methylated or methylated with HhaI, HpaII, or M.SssI prior to transfection as described in Materials and Methods. CT45 firefly luciferase/renilla luciferase (control) enzymatic activity data are shown. Data represent mean ± SD.
Discussion
We examined the role of DNA methylation in the regulation of CT45 expression in EOC. CT45 was expressed in ∼25% of EOC patient samples, and its expression was significantly increased in serous tumors, advanced stage, and high-grade lesions. Promoter hypomethylation of CT45 genes was observed in >50% of EOC samples, and this was exacerbated in advanced disease. Both CT45 expression and promoter hypomethylation were associated with reduced OS and DFS, consistent with their alteration during disease progression, and suggesting potential utility as prognostic biomarkers for EOC.We validated a link between DNA methylation alterations and CT45 expression in EOC both in clinical samples and EOC cell models. Importantly, CT45 promoter hypomethylation was predictive of both CT45 mRNA and CT45 protein expression. CT45 hypomethylation also showed strong correlation with LINE-1 hypomethylation, providing evidence that CT45 expression is linked to the global hypomethylation process that frequently impacts the cancer genome. While this link could be used to argue that CT45 (and more generally CT antigen) gene expression could be a passenger effect, the emerging oncogenic functions of CT antigen genes, including CT45, alternatively suggests that CT antigen gene expression is under positive selection in cancer. In the case of CT45, recent studies in Hodgkin's disease and breast cancer have established its function in promoting cell migration and invasion, as well as tumor growth and metastasis in vivo. The mechanism by which CT45 promotes these activities is not well understood, but may include regulation of the mRNA splicing machinery and/or transcriptional activation of critical pro-oncogenic genes such as TWIST. Our future studies will include investigation of the oncogenic role of CT45 proteins in EOC.Pharmacological inhibition of DNA methylation induced both CT45 mRNA and protein expression, and site-specific methylation of the CT45 promoter significantly downregulated its activity. Further work is needed to define additional epigenetic and genetic mechanisms that regulate CT45 expression in EOC and other malignancies, and to compare these to the mechanisms regulating other CT antigen genes. The more frequent occurrence of CT45 hypomethylation than CT45 expression in EOC additionally suggests that there may be a threshold of hypomethylation necessary to elicit CT45 expression. This threshold level could be distinct in each tumor, depending upon the genomic and epigenomic context.Treatment approaches for advanced EOC are challenged by the emergence of chemotherapy-resistant disease. The immune system may be uniquely suited to combat this problem, by virtue of its recognition of antigens specific to cancer tissues and due to its adaptive nature. In this context, the CT antigen gene family is a potentially fruitful immunological target in EOC and other cancers due to their highly restricted pattern of expression and inherent immunogenicity. Although spontaneous immune responses to CT45 in cancerpatients have not been commonly observed, this does not exclude its potential utility as a cancer vaccination target. In support of this idea, we observed a high titer antibody response to CT45 in a recent phase I clinical trial of doxil + decitabine + NY-ESO-1 vaccine in EOC. In addition, a recent study of immune modulation using anti-CTLA4 antibody in melanomapatients demonstrated seroconversion to CT45. However, our data suggest that CT45 based vaccination approaches might be most beneficial when occurring in the context of co-treatments that induce CT45 expression in vivo. Epigenetic modulators such as decitabine, SGI-110, or similar agents, may be well suited for utilization in this context. Finally, we point out that the emerging oncogenic functions of CT45 suggest that therapeutics that directly target this protein family could also have utility in cancers that commonly express these proteins, including ER-negative breast cancer and HGSOC. A key advantage of such an approach would be a high level of specificity for tumor cells over normal cells.
Materials and Methods
Human tissues
NO, OSE, FTE, and EOC samples were obtained under IRB-approved protocols at Roswell Park Cancer Institute as described previously. Cell extracts of clinical samples were prepared as described previously. EOC TMAs were constructed at Roswell Park Cancer Institute as described previously. EOC clinical characteristics, and the molecular assay data for each sample, are shown in Table S2.
Human cell lines
HumanOSE cells immortalized with SV40 large T antigen (IOSE121), and EOC cell lines (A2780, OVCAR3), and their culturing methods were described previously. SNU119 HGSOC cells were obtained from the Korean Cell Line Bank and cultivated in RPMI1640 media supplemented with 10% FBS and 1% Pen-Strep. OVSAHO HGSOC cells were obtained from the Japanese Collection of Research Bioresources, and cultivated in RPMI1640 media supplemented with 10% FBS and 1% Pen-Strep.
DNA, RNA, and protein extractions
Genomic DNA was isolated using the Puregene Tissue Kit (Qiagen), total RNA was purified using TRIzol® (Invitrogen), and total protein was extracted using RIPA buffer containing protease inhibitor cocktail (Sigma), according to standard methods.
Decitabine (DAC) treatment
DAC was obtained from Sigma. Cells at ∼20% confluence were treated with 1 μM DAC (day 0), passaged at day 2, re-treated with 1 μM DAC at day 3, and harvested for cell extracts at day 5. PBS was used as the vehicle control. We observed slight to moderate toxicity in EOC cells using this regimen.
RT-qPCR
RNA was processed using the DNA-free kit (Ambion), and cDNA synthesis was performed using the iScript cDNA Synthesis Kit (BioRad). Gene expression was determined by qPCR using a BioRad CFX Connect system and the SYBR green method. Standard curves were prepared using a mixture of cDNA from EOC cell lines. Pan-CT45 mRNA expression was normalized to 18s rRNA. Primer sequences are provided in Table S1.
Western blotting
Equal amounts of protein extracts (20 µg) were loaded per lane, separated using NuPAGE 4–12% Bis-Tris Gels (Invitrogen), and transferred to PVDF membranes (Roche) using the Trans-Blot SD Semi-Dry Transfer Cell (BioRad). Membranes were stained with FastGreen (Fisher) as a loading control. Western blotting was performed according to standard methods using humanCT45-specific mouse polyclonal antibody ab169438 (1:100 dilution) (Abcam), and human β-Actin mouse monoclonal antibody sc-47778 (1:2000 dilution) (Santa Cruz).
CT45 protein expression on EOC TMAs
EOC TMAs and CT45 protein staining by IHC was described previously. Representative CT45 IHC staining of 0.6 mm TMA cores is shown in .
RLM-RACE
The TSS of CT45 were determined using the FirstChoice RLM-RACE Kit (Ambion). RNA samples included human normal testis (BioChain Institute), and total RNA from a mixture of 2 different EOC samples (EOC 89, EOC 100). Primer sequences are provided in Table S1. PCR products were analyzed as described previously.
DNA methylation analyses
Bisulfite clonal sequencing and bisulfite pyrosequencing were used to determine the DNA methylation status of 5′ regions of CT45 genes. Sodium bisulfite conversions utilized the EZ DNA Methylation Kit (Zymo Research). CT45bisulfite sequencing primers were designed using MethPrimer, and targeted 4 CT45 family members: CT45A2, CT45A3, CT45A4 and CT45A5. Bisulfite pyrosequencing assays were designed based on the CT45A2 sequence, and assay 1 targets CT45A2, A3, A4 and A5; assay 2 targets CT45A2, A3, A4, A5 and A6. LINE-1 methylation was determined using bisulfite pyrosequencing as described previously. Primer sequences are provided in Table S1.
CT45 co-expression with other CT antigen genes in EOC
Oncomine (Compendia Bioscience, Ann Arbor, MI) was used to analyze CT45 expression in publically available human microarray data. Data are shown in for the analysis of CT45A3 co-expression in the Bittner humanovarian cancer data set, which used Affymetrix HG U133 Plus 2.0 arrays to profile 19,574 genes in 241 ovarian cancer samples of diverse histological origin (GEO accession number GSE2109).
CT45 promoter activity
Two CT45 promoter fragments (see ) were amplified using standard methods and cloned into pCR2.1 TOPO vector (Invitrogen). Inserts were then sub-cloned into the CpG-free vector, pCpGL-basic. Purified plasmids were methylated in vitro with M.SssI, HpaII, HhaI, or no enzyme (mock), in the presence of S-adenosyl methionine (all enzymes and co-factor obtained from NEB). Methylation efficiency was confirmed by digesting the methylated constructs using endonucleases MspI, McrBC, HhaI, and HpaII (NEB) (data not shown). Cell lines were transiently transfected with 500 ng CT45 firefly luciferase construct + 250 ng renilla luciferase control construct (pRL-TK) (Promega), using Lipofectamine 2000 (Invitrogen). Transfections were carried out in duplicate wells of 24 well plates. Cell extracts were harvested 24–48 hrs post-transfection, and firefly and renilla luciferase activity were measured using the Dual-Luciferase Reporter Assay System (Promega).
Survival analyses
OS was defined as the time between the date of diagnosis and death. Patients who were alive at the time of analysis were censored at the date of last follow-up. DFS was defined as the time between the date of surgery and disease recurrence. Patients who were alive and disease-free were censored at the date of last follow-up. For each molecular marker, samples were dichotomized at the median marker value, and the difference between the 2 groups was compared using Kaplan Meier analysis. The null hypothesis of no difference in the survival distributions was assessed using the Logrank test.
Analysis of CT45 mRNA expression and CNA in TGCA data
TCGA genomic data for HGSOC (Ovarian Serous Cystadenocarcinoma-TCGA, Provisional) was retrieved and analyzed using cBioPortal (data upload March 9, 2015). The following genomic profiles were analyzed: mutations, putative copy-number alterations from GISTIC using Onco Query Language (OQL), and mRNA expression data (RNA seq V2 RSEM). All parameters were set at default.
Authors: Hans-Jürgen Heidebrecht; Alexander Claviez; Marie Luise Kruse; Marc Pollmann; Friedrich Buck; Sönke Harder; Markus Tiemann; Wolfgang Dörffel; Reza Parwaresch Journal: Clin Cancer Res Date: 2006-08-15 Impact factor: 12.531
Authors: Michael A Postow; Jason J Luke; Mark J Bluth; Nikhil Ramaiya; Katherine S Panageas; Donald P Lawrence; Nageatte Ibrahim; Keith T Flaherty; Ryan J Sullivan; Patrick A Ott; Margaret K Callahan; James J Harding; Sandra P D'Angelo; Mark A Dickson; Gary K Schwartz; Paul B Chapman; Sacha Gnjatic; Jedd D Wolchok; F Stephen Hodi; Richard D Carvajal Journal: Oncologist Date: 2013-05-28
Authors: Ashok Sharma; Mustafa Albahrani; Wa Zhang; Christina N Kufel; Smitha R James; Kunle Odunsi; David Klinkebiel; Adam R Karpf Journal: Epigenetics Date: 2019-03-04 Impact factor: 4.528
Authors: Jianhang Zhao; Zhaoxu Xu; Yan Liu; Xiaobin Wang; Xinli Liu; Yanan Gao; Ying Jin Journal: Am J Cancer Res Date: 2022-02-15 Impact factor: 6.166
Authors: Amit M Oza; Ursula A Matulonis; Angeles Alvarez Secord; John Nemunaitis; Lynda D Roman; Sarah P Blagden; Susana Banerjee; William P McGuire; Sharad Ghamande; Michael J Birrer; Gini F Fleming; Merry Jennifer Markham; Hal W Hirte; Diane M Provencher; Bristi Basu; Rebecca Kristeleit; Deborah K Armstrong; Benjamin Schwartz; Patricia Braly; Geoff D Hall; Kenneth P Nephew; Simone Jueliger; Aram Oganesian; Sue Naim; Yong Hao; Harold Keer; Mohammad Azab; Daniela Matei Journal: Clin Cancer Res Date: 2019-12-12 Impact factor: 12.531
Authors: Charly Ryan Good; Shoghag Panjarian; Andrew D Kelly; Jozef Madzo; Bela Patel; Jaroslav Jelinek; Jean-Pierre J Issa Journal: Cancer Res Date: 2018-06-11 Impact factor: 12.701
Authors: Fabian Coscia; Ernst Lengyel; Jaikumar Duraiswamy; Bradley Ashcroft; Michal Bassani-Sternberg; Michael Wierer; Alyssa Johnson; Kristen Wroblewski; Anthony Montag; S Diane Yamada; Blanca López-Méndez; Jakob Nilsson; Andreas Mund; Matthias Mann; Marion Curtis Journal: Cell Date: 2018-09-20 Impact factor: 41.582