- Open Access
SFRP1 is a possible candidate for epigenetic therapy in non-small cell lung cancer
BMC Medical Genomics volume 9, Article number: 28 (2016)
Non-small cell lung cancer (NSCLC) remains a lethal disease despite many proposed treatments. Recent studies have indicated that epigenetic therapy, which targets epigenetic effects, might be a new therapeutic methodology for NSCLC. However, it is not clear which objects (e.g., genes) this treatment specifically targets. Secreted frizzled-related proteins (SFRPs) are promising candidates for epigenetic therapy in many cancers, but there have been no reports of SFRPs targeted by epigenetic therapy for NSCLC.
This study performed a meta-analysis of reprogrammed NSCLC cell lines instead of the direct examination of epigenetic therapy treatment to identify epigenetic therapy targets. In addition, mRNA expression/promoter methylation profiles were processed by recently proposed principal component analysis based unsupervised feature extraction and categorical regression analysis based feature extraction.
The Wnt/β-catenin signalling pathway was extensively enriched among 32 genes identified by feature extraction. Among the genes identified, SFRP1 was specifically indicated to target β-catenin, and thus might be targeted by epigenetic therapy in NSCLC cell lines. A histone deacetylase inhibitor might reactivate SFRP1 based upon the re-analysis of a public domain data set. Numerical computation validated the binding of SFRP1 to WNT1 to suppress Wnt signalling pathway activation in NSCLC.
The meta-analysis of reprogrammed NSCLC cell lines identified SFRP1 as a promising target of epigenetic therapy for NSCLC.
Non-small cell lung cancer (NSCLC) is still lethal despite many proposed therapeutic strategies. Among the many alternative strategies, epigenetic therapy is regarded as a promising method , and a histone deacetylase (HDAC) inhibitor  or DNA methyltransferase inhibitor  were shown to be promising NSCLC treatments, especially when combined . There has been extensive research regarding the clinical usefulness of epigenetic therapy for NSCLC; however, studies investigating the target genes of these treatments are limited, although some promising candidates have been proposed . The potential reasons for the small number of epigenetic therapy target gene reports might be the difficulty of in vitro studies . Compared with many clinical studies regarding the efficiency of epigenetic therapy, there have been few in vitro studies of epigenetic therapy [6, 7]. Thus, alternative strategies to direct in vitro experiments for epigenetic therapy such as the investigation of reprogrammed cancer cell lines are required to investigate the effect of epigenetic therapy in NSCLC.
It is thought that epigenetic therapy targets epigenetic effects, e.g., DNA methylation and/or histone modification, which might be affected by reprogramming. Thus, a detailed and extensive comparative study might indirectly identify the effect of epigenetic therapy in NSCLC cell lines.
This study performed a meta-analysis of reprogrammed NSCLC cell lines to identify genes associated with epigenetic alterations and expression changes during reprogramming and to identify promising candidate genes for targets of epigenetic therapy. Among those identified, secreted frizzled-related protein (SFRP)1 was of interest. Using in vitro epigenetic therapy experiments, we confirmed that SFRP1 mRNA expression and its histone modification were altered. Furthermore, SFRP1 might suppress the Wnt signalling pathway by binding to Wnt genes. An in silico study indicated the potential binding of SFRP1 with WNT1; thus, the reactivation of SFRP1 suppressed in NSCLC might be a promising candidate target for the epigenetic therapy of NSCLC.
Identification of biologically significant genes
To identify genes targeted by epigenetic therapy in NSCLC, we analysed gene expression and promoter methylation in reprogrammed NSCLC cell lines . Although it is useful to consider histone modification and promoter methylation together because epigenetic therapies targets both, suitable data sets were not publically available for histone modification; therefore, as promoter methylation often reflects the effect of histone modification , a data set containing gene expression and promoter methylation information was analysed. The primary aim of this analysis was to identify genes associated with aberrant gene expression and promoter methylation during reprogramming because associated genes are most likely targeted by epigenetic therapy.
Although promoter methylation was generally expected to be negatively correlated with gene expression, this was not always observed, especially when histone modification was also considered . Because this study aimed to identify targets of epigenetic therapy including both DNA methylation and histone modification, we did not restrict candidate biologically significant genes such as those associated with negative correlations between promoter methylation and gene expression, but considered all genes associated with significant correlations between promoter methylation and gene expression independent of the direction.
To select biologically significant genes, we used principal component analysis (PCA) based unsupervised feature extraction (FE) [11–24]. PCA based unsupervised FE is useful when there is no information regarding how to order multiple classes. It also allows us to restrict number of pairs whose correlations must be computed, which can reduce the possibility that selected genes are rejected because of P-values adjustments based on multiple comparison correction criteria. Therefore, because many cell lines, including those that were reprogrammed and differentiated, were used in this study, PCA based unsupervised FE was a suitable method for analysis. To select principal components (PCs) with a significant correlation between gene expression and promoter methylation for FE, we performed hierarchical clustering (see Methods) to identify a pair of PCs associated with a high correlation between promoter methylation and gene expression. PC3 and PC4 were the most suitable candidate pairs (Fig. 1).
One may wonder why H1, IMR90, and iPSIMR90 must be included in the analysis. As can be seen in Additional file 1, reprogrammed cell lines have similar values to H1 or iPSIMR90 that are pluripotent; this suggested that the inclusion of H1 and iPSIMR90 guaranteed that cell lines were correctly reprogrammed. Similarly, inclusion of IMR90 guaranteed that reprogrammed cell lines are often distinct from IMR90 that is differentiated. Thus, H1, IMR90, and iPSIMR90 are worthwhile being included.
To determine the stability of pairs of PC3 and PC4 between gene expression and promoter methylation we constructed hierarchical clustering with only 23 samples (because there are 24 samples, there was a sequential removal of one specific sample from the 24 samples; see Methods and full results are Additional file 2). Pairs of PC3 and PC4 between gene expression and promoter methylation were conserved for 22 hierarchical clusters among 24 samples. Thus, the pairs of PC3 and PC4 in Fig. 1 are not accidental but robust.
Although we previously  considered only PC1 and PC2 when analysing the same data set for another study, gene expression and promoter methylation showed a lower correlation when compared with PC3 and PC4 in this study. This might explain why promising candidate genes were not identified in our previous study . Therefore, we used PC3 and PC4 for gene selection in this study. In addition to PCA based unsupervised FE, we used another FE that is also suitable for multiclasses that lack a pre-decided order, FE based upon categorical regression (see Methods).
Table 1 summarizes the genes selected by PCA based unsupervised FE and categorical regression based FE. In total, we identified 32 unique gene candidates (three genes were identified by more than one method). Gene expression and promoter methylation of specific PCs and genes and their correlation information is summarized in Table 1 and detailed in Additional files 1, 3, 4 and 5.
Biological significance of selected genes
Disease association of genes
To validate the biological significance of the selected genes, we used the Gendoo server  to search the literature for genes associated with diseases. For most of the genes examined (excluding LAD1, KIF1A, SLC16A12, SCG3 and IGSF21), there were significant associations with cancer-related diseases as summarized in Table 1A (and Additional file 6). Many oncogenes and tumor suppressors are not unique for specific cancers but are related to several cancers. This suggested that the genes selected in this study have the potential to be involved with NSCLC tumourigenesis.
Pathway/Gene Ontology (GO) term analysis
Because disease association is not always informative regarding how the genes are involved in tumourigenesis, we uploaded a list including gene IDs to two gene annotation servers [26, 27] (Table 2). The selection of “extracellular region” was reasonable, because this is a reprogramming study, thus cell surface receptors should be activated to initiate differentiation, which is related to another selected GO term, “cell proliferation”. When mapping selected genes to the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, most were at the cell surface (see Additional file 7).
To determine whether the selected genes were specifically related to NSCLC and pluripotency because the data analysed was from reprogramming experiments, we performed an extensive literature search. Studies regarding the relationship between NSCLC and proteins reported to bind to any of the 32 genes listed in Table 1 were collected and analysed with BioGrid , which reports literature-based protein-protein interactions. Most of the genes identified were specifically related to NSCLC tumourigenesis and some were also related to pluripotency (Table 1C and Additional file 8). Thus, this methodology is promising.
The Wnt/β-catenin signalling pathway as an NSCLC therapy target
Although we found that the genes identified in this study were mostly related to NSCLC tumourigenesis, genes should be selected according to their potential for epigenetic therapy of NSCLC. However, as can be seen in Table 2, no significant pathway enrichments were detected. In order to investigate pathway enrichment, we performed literature searches manually. Then we have found that multiple genes selected in this study were related to the Wnt/β-catenin signalling pathway (Table 1D) that was recently reported to be a major pathway in NSCLC tumourigenesis . One may wonder why Wnt/β-catenin signalling pathway was not detected in enrichment analyses in Table 2. First of all, even if no significant enrichments were detected, it does not always mean the lack of enrichment, but often simply means the lack of ability of the specific statistical tests. Second, as can be seen in the following, some genes detected by literature searches, e.g., EPCAM and TACSTD2, are not included into KEGG pathway. This means that we need more sophisticated investigations than simple enrichment analyses. This is the reason why we additionally performed literature searches.
SALL4 is part of the Wnt signalling pathway  and regulates the stemness of EPCAM-positive hepatocellular carcinomas [31, 32]. EPCAM was recently reported to be an endoderm-specific Wnt derepressor . ANGPT1 was reported to be upregulated via the overexpression of β-catenin that is a key factor of the Wnt signalling pathway . TACSTD2 was proposed to be a Wnt target  identified through consistent gene expression changes in APC-mutant intestinal adenomas from humans and mice; EFNB1 :Eph-related receptor is a Wnt signalling target gene in colorectal cancer  that binds to the EFNB1 ligand. MEST inhibits Wnt signalling through the regulation of LRP6 glycosylation . F2R(PAR1) stabilizes β-catenin in mammary gland tissues . DKK3 binds to LRP5/6 and inhibits the initiation of Wnt signalling . SFRP1 binds to FZD and WNT to suppress the activation of Wnt signalling . HOXA5 expression increased the retention of β-catenin in adherens junctions and reduced permeability . KIF1A binds to at least two β-catenin binding proteins , ESR1  and AR . TM4SF1 might have a role in coordinating Wnt signalling and migration during endocrine pancreas specification , and TM4SF1 and TM4SF4 belong to the tetraspanin L6 domain family. GPR56: The Wnt/β-catenin signalling pathway regulates genes involved in cell proliferation, survival, migration and invasion through by the regulation of T-cell factor (TCF)-4 transcription factor proteins that activate GPR56 in HCC . S100P: Increased expression of S100P promoted cellular proliferation by increasing the nuclear translocation of β-catenin in endometrial cancer . SPINT2: The epigenetic silencing of SPINT2 promoted cancer cell motility via HGF-MET pathway activation in melanoma , and β-catenin formed a complex with c-Met (HGF receptor) . CDH1 (E-cadherin) is involved in the inactivation of Wnt/β-catenin signalling in urothelial carcinoma and normal urothelial cells . LAMC2 (Laminin γ2) mediated the Wnt5a-induced invasion of gastric cancer cells . HMGA1 interacted with β-catenin to positively regulate Wnt/β-catenin signaling in colorectal cancer cells . PFKFB3: The altered expression of PFKFB3 is associated with Wnt signalling pathway genes . UCHL1 is a colorectal cancer oncogene that activated the β-catenin/TCF pathway through its deubiquitinating activity . ALDH3A1 is overexpressed in a subset of hepatocellular carcinoma characterized by activation of the Wnt/β-catenin pathway . EPB41L3 (DAL1) binds to YWHAZ , and the YWHAZ/β-catenin axis promoted epithelial-mesenchymal transition and lung cancer metastasis . LAMA1 (laminin): Overexpression of the Wnt antagonist FRZB1 decreased RNA levels of the essential basement membrane genes fibronectin and laminin .
β-catenin is often reported to be related to NSCLC
Although β-catenin is extensively related to the selected genes in this study, β-catenin was overexpressed in NSCLC . β-catenin expression was also prognostic for improved NSCLC survival . Nuclear β-catenin accumulation was associated with the increased expression of NANOG protein and predicted a poor prognosis of NSCLC . Promoter methylation-mediated silencing of β-catenin enhanced the invasiveness of NSCLC and predicted an adverse prognosis . All of these studies strongly suggest the importance of β-catenin in NSCLC.
These genes are also related to epigenetic therapy
The following genes were also suggested to be related to epigenetic therapy. Recently, the combined usage of two drugs, romidepsin and decitabine, restored SFRP1 activity in four cancer cell lines, A498, KIJ265T, MDA-231, and BT-20 . This strategy might be useful for NSCLC therapy because an HDAC inhibitor, MPT0E028, enhanced erlotinib-induced cell death in epidermal growth factor receptor tyrosine kinase inhibitor (EGFR-TKI)-resistant NSCLC cells  and SAHA, a HDAC inhibitor, had profound anti-growth activity against NSCLC cells . Other evidence includes an organosulfur derivative of the HDAC inhibitor, valproic acid, which sensitised human lung cancer cell lines to apoptosis and to cisplatin cytotoxicity . EGFR-TKI resistance by BIM polymorphism was circumvented when combined with HDAC inhibition , and antitumour activity of histone deacetylase inhibitors was observed in NSCLC cells . The effect of HDAC inhibitors can be improved by in silico drug screening . In addition, SALL4 was reported to form a protein complex with HDAC [67–69] (Fig. 4a in ), which suggests that epigenetic regulation of the Wnt signalling pathway is a key factor in the tumourigenesis of NSCLC. Interestingly, promoters of SALL4 and SFRP1 were methylated simultaneously [70–72]. Although there have been no reports to suggest a direct relationship between HOXA5 and the Wnt signalling pathway in NSCLC, HOX is related to the Wnt signalling pathway, which controls HOX gene expression in C. elegans , while HOX genes control Wnt signalling . Furthermore, WNT7A has a strong relationship with HOX genes . In addition, from an evolutionary point of view, HOX and Wnt might be related . Thus, HOXA5 might be involved in Wnt signalling in NSCLC and might also be influenced by HDAC .
SFRP1 is a potential epigenetic therapy target
Overall, we concluded that the Wnt signalling pathway is a likely target of epigenetic therapy in NSCLC cell lines. A previous study suggested that the reactivation of Wnt antagonists, including DKK3 and SFRP1, in NSCLC might be a beneficial epigenetic therapy . Among the genes potentially related to Wnt signalling, we considered SFRP1 to be the most promising candidate gene targeted by epigenetic treatment. A previous study reported that treatment with 5-aza-2’-deoxycytidine, a DNA methyltransferase inhibitor, enhanced SFRP1 expression in NSCLC . Transcriptional silencing of the gene was also due to hypermethylation of its promoter region in NSCLC . SFRP1 has been reported as a marker for NSCLC [81, 82]. In addition, SFRP1 was also reported to be related to β-catenin. For example, SFRP1, SFRP2, and SFRP5 regulate Wnt/β-catenin and planar cell polarity pathways during early trunk formation in mice . Loss of SFRP1 expression was associated with aberrant β-catenin distribution and tumor progression in mucoepidermoid carcinoma of salivary glands .
To confirm whether the HDAC inhibitor affected SFRP1 regulation in NSCLC, we analysed a public domain data set. Miyanaga et al.  compared various cell lines to determine whether they were resistant to HDAC inhibitors. We investigated SFRP1 expression between HDAC inhibitor-resistant cell lines and non-resistant cell lines for adenocarcinoma and squamous cell carcinoma and found different levels of SFRP1 expression (Table 3). SFRP1 expression was upregulated in non-resistant cell lines compared with resistant cell lines because cell lines with downregulated SFRP1 required greater levels of HDAC suppression to increase SFRP1 expression. In addition, histone acetylation of SFRP1 in NSCLC was enhanced by HDAC inhibitors compared with DKK3 and TACSTD1 , but not in CL1-1 generated from the cervix. These results are in accord with the hypothesis that the therapeutic effect of HDAC in NSCLC is caused by the reactivation of SFRP1. Interestingly, the histone acetylation of SALL4 in NSCLC was also enhanced by the HDAC inhibitor in an A549 cell line (Table 3, P-values for the CL1-1 cell line were very small, but because of the deacetylation, this effect is not likely to be caused by the HDAC inhibitor). Unfortunately, the microarray analysis by Miyagawa et al. , did not include SALL4, thus we cannot check whether SALL4 expression was coincident with HDAC inhibitor resistance.
Potential of SFRP1 binding to WNT1
Next, we validated the ability of SFRP1 to bind to WNT1 as it was the most promising target from our study that affected the Wnt signalling pathway. Although Wnt8 and FZ8 share a cysteine-rich domain (CRD) that forms a protein complex with SFRP1  and FRZB1 that is a similar protein to SFRP1  was reported to bind to WNT1 in Xenopus , there have been no direct reports investigating the binding of SFRP1 to WNT1 . Therefore, we examined the formation of a SFRP1-WNT1 protein complex using numerical simulation. The tertiary structures of WNT1 and SFRP1 were inferred by RaptorX . Then, the obtained structures were uploaded to the ZDOCK server , a rigid body based protein complex predictor. The 10 top ranked protein complex structures obtained were further uploaded to Fiberdock  that refines the protein complex structures obtained by other methods by considering the flexible structures of the proteins. Finally, the best candidate (with the minimum energy) reported by Fiberdock was identified as the most reliable candidate for the WNT1-SFRP1 protein complex. Figure 2a shows the structure of the protein complex obtained using this procedure. This structure is very similar to the WNT8-FZ8 complex (Fig. 2d) because the CRDs of SFRP1 and FZ8 were similarly sandwiched by two arms of the Wnt protein. This suggests that SFRP1 can bind to WNT1 and suppress the Wnt signalling pathway.
To confirm the reliability of this structure, we performed two tests. The first was to upload WNT8 and the CRD of FZ8 separately to ZDOCK and Fiberdock, as for the protein tertiary structures of SFRP1 and WNT1 obtained by RaptorX, to determine whether the correct structure was identified as that with the minimum energy. Figure 2b shows a protein complex identified by the combined usage of ZDOCK and Fiberdock. The expected structure contains an FZ8 CRD sandwiched by two arms of WNT8, thus, this confirms the use of our strategy using the combined ZDOCK and Fiberdock systems.
The second test was a molecular dynamics (MD) simulation to test the stability of the protein complex predicted by the combined usage of ZDOCK and Fiberdock. The protein complex inferred by ZDOCK and Fiberdock was used as the initial state and MD simulation was performed by GROMACS . The obtained structure was modified to have minimum energy and was simulated under NVT (constant number of molecules, volume and temperature) and NPT (constant number of molecules, pressure and temperature) conditions over 0.1 ns, respectively. Then, a 2 ns equilibration MD was performed. Figure 3 shows the time developments of the root mean square deviation (RMSD) during the first and second 1 ns in equilibration MD. Although the structure heavily fluctuates because RMSD increased with time, the SFRP1 CRD structure was sandwiched by two arms of WNT1 and was maintained even after 2 ns equilibration MD (Fig. 2c) indicating this structure was stable. Thus, SFRP1 binds to WNT1 to suppress the Wnt signalling pathway.
Possibilities that proteins other than SFRP1 are epigenetic therapy targets
Although this study focused on SFRP1 as a promising candidate for epigenetic therapy of NSCLC, 31 other proteins selected in Table 1 might also have potential for NSCLC epigenetic therapy. SALL4, DKK3 and HOXA5 were suggested to be related to the Wnt signalling pathway. DKK3 suppresses the Wnt signalling path-way by binding to LPR5/6 proteins. Other genes also contribute to many other pathways involved in tumourigenesis (Additional files 6 and 8). Thus, although no strong evidence that proteins other than SFRP1 were epigenetic therapy targets for NSCLC was obtained, it is likely that other proteins identified in this study are epigenetic therapy targets.
In conclusion, this meta-analysis of reprogrammed NSCLC cell lines indicated that SFRP1 was a promising candidate for NSCLC epigenetic therapy. The reliability of SFRP1 binding to WNT1 to suppress the Wnt signalling pathway was confirmed using numerical investigations. Thus, we propose that SFRP1 is an epigenetic therapy target for NSCLC.
Gene expression and promoter methylation profiles during reprogramming of NSCLC cell lines
Gene expression and promoter methylation profiles were downloaded from the Gene Expression Omnibus (GEO) using GEO ID: GSE35913. Files including gene expression, GSE35911_SampleProbeProfile.txt.gz, were provided as a supplementary file in the subseries GEO ID: GSE35911. Columns annotated as “AVG Signal” were used. Promoter methylation profiles were obtained from “Series Matrix File(s)” in the subseries GEO ID: GSE35912. They consisted of eight cell lines, H1 (ES cell), H358 and H460 (NSCLC), IMR90 (Human Caucasian foetal lung fibroblast), iPCH358, iPCH460, iPSIMR90 (reprogrammed cell lines), and piPCH358 (re-differentiated iPCH358) with three biological replicates. In total, there were 3 replicates × 8 cell lines × 2 properties (gene expression and promoter methylation) = 48 samples. No further normalization processes were applied. Multiple probes attributed to the same gene were not integrated before feature selection and if genes with more than one attributable probes were selected, the genes were regarded as extracted features.
PCA-based unsupervised FE
PCA based unsupervised FE was performed as previously described . Briefly, x ij is ith feature that represents gene expression/promoter methylation/histone modification of the jth sample. In contrast to the standard usage of PCA, features but not samples were embedded into low dimensional space. Then, kth PC score attributed to ith feature, x ik, was expressed as
where c jk is the contribution of jth sample to kth PC (PC loadings) and 〈···〉 i is the mean over i. After identifying biologically critical PCs (in this study, PC3 and PC4, based on hierarchical clustering, Fig. 1), features that are outliers along the specified PC were extracted, because outliers were expected to reflect biological significance (in this study, there was a high correlation between gene expression and promoter methylation) that specified PCs represent. The number of features to be extracted as outliers was decided empirically. In this study, the top 300 outliers, those with larger or smaller (larger absolute values of) x ik of kth PC selected for FE, were selected for gene expression and promoter methylation profiles, respectively. Genes listed in Table 1 were those commonly selected as the top ranked outliers for gene expression and methylation profiles in the following four combinations of rankings: larger x ik for gene expression and larger x ik for promoter methylation; larger x ik for gene expression and smaller (negatively larger) x ik for promoter methylation; smaller (negatively larger) x ik for gene expression and larger x ik for promoter methylation; and smaller (negatively larger) x ik for gene expression and smaller (negatively larger) x ik for promoter methylation.
Identification of biologically meaningful PCs
Although there are several ways to identify biologically meaningful PCs, this study used the coincidence of c jk between gene expression and promoter methylation. To identify mostly coincident PCs between gene expression and promoter methylation, we computed the correlation coefficients between kth and k’th PCs of gene expression or promoter methylation as follows: (all pairs were considered among PCs of gene expression and promoter methylation)
where 〈···〉 j is the mean over j. Then, the negative signed absolute value of ρ kk′ , − |ρ kk′ |, was used as the distance for hierarchical clustering (Unweighted Pair Group Method with Arithmetic mean, [UPGMA], was employed). If a pair of PCs for gene expression and promoter methylation were clustered together with a smaller distance, i.e., with a larger absolute value of correlation coefficient, we employed the pair of PCs of gene expression and promoter methylation, for FE.
Hierarchical clustering with 23 samples
PCs were computed with 23 samples and hierarchical clustering was performed. As a result, we had 24 hierarchical clusters. Because the correspondence between PCs with a distinct set of 23 samples is incomplete, we re-labelled the first five PCs (PC1, PC2, PC3, PC4 and PC5) and compared them with PCs obtained using 24 samples.
Categorical regression-based FE
Categorical regression-based FE was defined as follows: x ij reflects the ith feature of the jth samples as described above; therefore x ij can be represented as:
with δ jl = 1 only when the jth sample belongs to the lth category (in this study, each category corresponds to the type of cell line), otherwise it is 0. Category summation was determined and a jl s were the fitting parameters. Because independent variables are categorical, the above regression equation belongs to a category of equations often named categorical regression. For each ith feature, P-values were computed using the lm function implemented in R  (this can be easily performed if factors corresponding to the types of cell lines are used as independent variables in lm). Genes were ranked based upon obtained P-values and the top 300 ranked significant genes (those with smaller P-values) were extracted. Genes listed in Table 1 were those commonly selected as top ranked genes for gene expression and methylation profiles.
Disease associations with genes
Disease associations with genes were investigated by Gendoo , a literature-based disease-gene association database.
For gene expression, probe annotations were based on the “Accession” column (for RefSeq gene ID) or “Symbol” column (for gene symbol) in the GSE35911_SampleProbeProfile.txt.gz file. For promoter methylation, GPL849065.txt available from the GEO ID: GSE35912 file was used and the “Accession” column was used to assign a Refseq gene ID to each probe.
Gene expression for comparison of resistant and non-resistant cell lines
Gene expression profiles used for the comparison between resistant and non-resistant cell lines were obtained from GSE4127 . The data set included in the “Series Matrix File(s)” was used for analysis without further normalization. RERF-LC-MS, PC14, PC9, A549, and RERF-LC-KJ were regarded as non-resistant cell lines and CP7, ABC-1, PC3 and LC2/ad were regarded as resistant cell lines.
Histone modification during treatment with an HDAC inhibitor
Histone modification profiles used for analysis during treatment with an HDAC inhibitor were obtained from GSE20304 . Data sets included in the “Series Matrix File(s)” were used for analysis without further normalization.
Inference of protein tertiary structures
Amino acid sequences extracted from Uniprot (Q8N474.1 for SFRP1 HUMAN and P04628.1 for WNT1 HUMAN) were uploaded to the RaptorX server and inferred protein structures were used for further analyses, i.e., uploading to ZDOCK, Fiberdock and MD by GROMACS.
MD by GROMACS
GROMACS 5.0.4 compiled with enabling mpi (Message Passing Interface) and GPU usage was used for MD. The protein complex of SFRP1 and WNT1 inferred by the combined usage of ZDOCK and Fiberdock (Fig. 2a) was used as the initial structure of the protein complex. Force field used was OPLS-AA/L (all-atom force field) and group 13 “SOL” was employed for embedding ions. At first, energy minimization was performed and NVT and NPT conditions of the simulation followed. Finally, an equilibration run was executed over 2 ns. For all the procedures, we followed the instructions shown in .
Mehta A et al. Epigenetics in lung cancer diagnosis and therapy. Cancer Metastasis Rev. 2015;34(2):229–41.
Huffman K, Martinez ED. Pre-clinical studies of epigenetic therapies targeting histone modifiers in lung cancer. Front Oncol. 2013;3:235.
Momparler RL. Epigenetic therapy of non-small cell lung cancer using decitabine (5-aza-2’-deoxycytidine). Front Oncol. 2013;3:188.
Jakopovic M et al. Targeting the epigenome in lung cancer: expanding approaches to epigenetic therapy. Front Oncol. 2013;3:261.
Vendetti FP et al. Evaluation of azacitidine and entinostat as sensitization agents to cytotoxic chemotherapy in preclinical models of non-small cell lung cancer. Oncotarget. 2015;6(1):56–70.
Kovalev RA et al. Possibilities of epigenetic anti-tumor therapy in in-vitro models. Vopr Onkol. 2012;58(6):800–7.
Gore SD. In vitro basis for treatment with hypomethylating agents and histone deacetylase inhibitors: can epigenetic changes be used to monitor treatment? Leuk. Res. 2009;33 Suppl 2:2–6.
Mahalingam D et al. Reversal of aberrant cancer methylome and transcriptome upon direct reprogramming of lung cancer cells. Sci Rep. 2012;2:592.
Hohmann CF, Ebner FF. Basal forebrain lesions facilitate adult host fiber ingrowth into neocorticaltransplants. Brain Res. 1988;448(1):53–66.
Wagner JR et al. The relationship between DNA methylation, genetic and expression inter-individual variation in untransformed human fibroblasts. Genome Biol. 2014;15(2):37.
Taguchi YH. Identification of aberrant gene expression associated with aberrant promoter methylation in primordial germ cells between E13 and E16 rat F3 generation vinclozolin lineage. BMC Bioinformatics. 2015;16 Suppl 18:S16.
Murakami Y et al. Comprehensive analysis of transcriptome and metabolome analysis in Intrahepatic Cholangiocarcinoma and Hepatocellular Carcinoma. Scientific Reports. 2015;5:16294.
Taguchi YH, et al. Heuristic principal component analysis-based unsupervised feature extraction and its application to gene expression analysis of amyotrophic lateral sclerosis data sets. IEEE Conf Comput Intell Bioinf Comput Biol. 2015. http://dx.doi.org/10.1109/CIBCB.2015.7300274.
Taguchi YH et al. Principal component analysis-based unsupervised feature extraction applied to in silico drug discovery for posttraumatic stress disorder-mediated heart disease. BMC Bioinformatics. 2015;16(1):139.
Taguchi Y-h et al. Heuristic principal component analysis-based unsupervised feature extraction and its application to bioinformatics. In: Wang B, Li R, Perrizo W, editors. Big Data Analytics in Bioinformatics and Healthcare. 2015. p. 138–62.
Taguchi Y-h, Okamoto A. Principal component analysis for bacterial proteomic analysis. In: Shibuya T, Kashima H, Sese J, Ahmad S, editors. Pattern Recognition in Bioinformatics. LNCS, vol. 7632. Heidelberg: Springer; 2012. p. 141–52.
Murakami Y et al. Comprehensive miRNA expression analysis in peripheral blood can diagnose liver disease. PLoS ONE. 2012;7(10):48366.
Ishida S et al. Bioinformatic Screening of Autoimmune Disease Genes and Protein Structure Prediction with FAMS for Drug Discovery. Protein Pept Lett. 2014;21(8):828–39.
Taguchi YH, Murakami Y. Principal component analysis based feature extraction approach to identify circulating microRNA biomarkers. PLoS ONE. 2013;8(6):66714.
Kinoshita R et al. Genes associated with genotype-specific DNA methylation in squamous cell carcinoma as candidate drug targets. BMC Syst Biol. 2014;8 Suppl 1:S4.
Taguchi YH, Murakami Y. Universal disease biomarker: can a fixed set of blood microRNAs diagnose multiple diseases? BMC Res Notes. 2014;7:581.
Murakami Y et al. Comparison of Hepatocellular Carcinoma miRNA Expression Profiling as Evaluated by Next Generation Sequencing and Microarray. PLoS ONE. 2014;9(9):106314.
Umeyama H, Iwadate M, Taguchi YH. TINAGL1 and B3GALNT1 are potential therapy target genes to suppress metastasis in non-small cell lung cancer. BMC Genomics. 2014;15 Suppl 9:S2.
Taguchi Y-h. Integrative analysis of gene expression and promoter methylation during reprogramming of a non-small-cell lung cancer cell line using principal component analysis-based unsupervised feature extraction. In: Huang D-S, Han K, Gromiha M, editors. Intelligent Computing in Bioinformatics. LNCS, vol. 8590. Heidelberg: Springer; 2014. p. 445–55.
Nakazato T, Bono H, Matsuda H, Takagi T. Gendoo: functional profiling of gene and disease features using MeSH vocabulary. Nucleic Acids Res. 2009;37(Web Server issue):166–9.
Chen YA, Tripathi LP, Mizuguchi K. TargetMine, an integrated data warehouse for candidate gene prioritisation and target discovery. PLoS ONE. 2011;6(3):17844.
Reimand J, Arak T, Vilo J. g:Profiler–a web server for functional interpretation of gene lists (2011 update). Nucleic Acids Res. 2011;39(Web Server issue):307–15.
Chatr-Aryamontri A et al. The BioGRID interaction database: 2015 update. Nucleic Acids Res. 2015;43(Database issue):470–8.
Stewart DJ. Wnt signaling pathway in non-small cell lung cancer. J Natl Cancer Inst. 2014;106(1):356.
Bohm J et al. SALL4 is directly activated by TCF/LEF in the canonical Wnt signaling pathway. Biochem Biophys Res Commun. 2006;48(3):898–907.
Zeng SS et al. The transcription factor SALL4 regulates stemness of EpCAM-positive hepatocellular carcinoma. J Hepatol. 2014;60(1):127–34.
Lu H et al. EpCAM is an endoderm-specific Wnt derepressor that licenses hepatic development. Dev Cell. 2013;24(5):543–53.
Zhang Y et al. Aberrant expression of a beta-catenin gain-of-function mutant induces hyperplastic transformation in the mouse cornea. J Cell Sci. 2010;123(Pt 8):1285–94.
Segditsas S et al. Putative direct and indirect Wnt targets identified through consistent gene expression changes in APC-mutant intestinal adenomas from humans and mice. Hum Mol Genet. 2008;17(24):3864–75.
Clevers H, Batlle E. EphB/EphrinB receptors and Wnt signaling in colorectal cancer. Cancer Res. 2006;66(1):2–5.
Jung H, Lee SK, Jho EH. Mest/Peg1 inhibits Wnt signalling through regulation of LRP6 glycosylation. Biochem J. 2011;436(2):263–9.
Yin YJ et al. Mammary gland tissue targeted overexpression of human protease-activated receptor 1 reveals a novel link to beta-catenin stabilization. Cancer Res. 2006;66(10):5224–33.
Arderiu G et al. HoxA5 stabilizes adherens junctions via increased Akt1. Cell Adh Migr. 2007;1(4):185–95.
Albers M et al. Automated yeast two-hybrid screening for nuclear receptor-interacting proteins. Mol Cell Proteomics. 2005;4(2):205–13.
Kouzmenko AP et al. Wnt/beta-catenin and estrogen signaling converge in vivo. J Biol Chem. 2004;279(39):40255–8.
Lee E et al. Inhibition of androgen receptor and β-catenin activity in prostate cancer. Proc Natl Acad Sci USA. 2013;110(39):15710–5.
Anderson KR et al. The L6 domain tetraspanin Tm4sf4 regulates endocrine pancreas differentiation and directed cell migration. Development. 2011;138(15):3213–24.
Stubbe P, von zur Muhlen A, Heidemann P. LH-RH Test in prepuberal children (author’s transl). Monatsschr Kinderheilkd. 1975;123(3):100–5.
Guo L et al. The expression of S100P increases and promotes cellular proliferation by increasing nuclear translocation of β-catenin in endometrial cancer. Int J Clin Exp Pathol. 2014;7(5):2102–12.
Hwang S et al. Epigenetic Silencing of SPINT2 Promotes Cancer Cell Motility via HGF-MET Pathway Activation in Melanoma. J Invest Dermatol. 2015;135(9):2283–91.
Apte U et al. Activation of Wnt/beta-catenin pathway during hepatocyte growth factor-induced hepatomegaly in mice. Hepatology. 2006;44(4):992–1002.
Thievessen I et al. E-cadherin involved in inactivation of WNT/beta-catenin signalling in urothelial carcinoma and normal urothelial cells. Br J Cancer. 2003;88(12):1932–8.
Yamamoto H et al. Laminin gamma2 mediates Wnt5a-induced invasion of gastric cancer cells. Gastroenterology. 2009;137(1):242–52.
Xing J, Cao G, Fu C. HMGA1 interacts with β-catenin to positively regulate Wnt/β-catenin signaling in colorectal cancer cells. Pathol Oncol Res. 2014;20(4):847–51.
Parker GE, Pederson BA, Obayashi M, Schroeder JM, Harris RA, Roach PJ. Gene expression profiling of mice with genetically modified muscle glycogen content. Biochem J. 2006;395(1):137–45.
Zhong J et al. UCHL1 acts as a colorectal cancer oncogene via activation of the β-catenin/TCF pathway through its deubiquitinating activity. Int J Mol Med. 2012;30(2):430–6.
Calderaro J et al. ALDH3A1 is overexpressed in a subset of hepatocellular carcinoma characterised by activation of the Wnt/β-catenin pathway. Virchows Arch. 2014;464(1):53–60.
Ewing RM et al. Large-scale mapping of human protein-protein interactions by mass spectrometry. Mol Syst Biol. 2007;3:89.
Chen CH, Chuang SM, Yang MF, Liao JW, Yu SL, Chen JJ. A novel function of YWHAZ/β-catenin axis in promoting epithelial-mesenchymal transition and lung cancer metastasis. Mol Cancer Res. 2012;10(10):1319–31.
Dickinson AJ, Sive HL. The Wnt antagonists Frzb-1 and Crescent locally regulate basement membrane dissolution in the developing primary mouth. Development. 2009;136(7):1071–81.
Retera JM, Leers MP, Sulzer MA, Theunissen PH. The expression of beta-catenin in non-small-cell lung cancer: a clinicopathological study. J Clin Pathol. 1998;51(12):891–4.
Chiu CG et al. Beta-catenin expression is prognostic of improved non-small cell lung cancer survival. Am J Surg. 2012;203(5):654–9.
Li XQ et al. Nuclear β-catenin accumulation is associated with increased expression of Nanog protein and predicts poor prognosis of non-small cell lung cancer. J Transl Med. 2013;11:114.
Miao Y et al. Promoter methylation-mediated silencing of β-catenin enhances invasiveness of non-small cell lung cancer and predicts adverse prognosis. PLoS ONE. 2014;9(11):112258.
Cooper SJ et al. Reexpression of tumor suppressor, sFRP1, leads to antitumor synergy of combined HDAC and methyltransferase inhibitors in chemoresistant cancers. Mol Cancer Ther. 2012;11(10):2105–15.
Chen MC et al. The HDAC inhibitor, MPT0E028, enhances erlotinib-induced cell death in EGFR-TKI-resistant NSCLC cells. Cell Death Dis. 2013;4:810.
Komatsu N et al. SAHA, a HDAC inhibitor, has profound anti-growth activity against non-small cell lung cancer cells. Oncol Rep. 2006;15(1):187–91.
Tesei A et al. Organosulfur derivatives of the HDAC inhibitor valproic acid sensitize human lung cancer cell lines to apoptosis and to cisplatin cytotoxicity. J Cell Physiol. 2012;227(10):3389–96.
Nakagawa T et al. EGFR-TKI resistance due to BIM polymorphism can be circumvented in combination with HDAC inhibition. Cancer Res. 2013;73(8):2428–34.
Miyanaga A et al. Antitumor activity of histone deacetylase inhibitors in non-small cell lung cancer cells: development of a molecular predictive model. Mol Cancer Ther. 2008;7(7):1923–30.
Tambunan US et al. In silico modification of suberoylanilide hydroxamic acid (SAHA) as potential inhibitor for class II histone deacetylase (HDAC). BMC Bioinformatics. 2011;12 Suppl 13:S23.
Lu J et al. Stem cell factor sall4 represses the transcriptions of pten and sall1 through an epigenetic repressor complex. PLoS ONE. 2009;4:5577.
Gao C et al. Targeting transcription factor SALL4 in acute myeloid leukemia by interrupting its interaction with an epigenetic complex. Blood. 2013;121(8):1413–21.
Wang F et al. The embryonic stem cell gene sall4. Mol Cell Oncol. 2014;1(4):969169.
Sugai T et al. Molecular analysis of gastric differentiated-type intramucosal and submucosal cancers. Int J Cancer. 2010;127(11):2500–9.
Vincent A et al. Genome-wide analysis of promoter methylation associated with gene expression profile in pancreatic adenocarcinoma. Clin Cancer Res. 2011;17(13):4341–54.
McGarvey KM et al. Defining a chromatin pattern that characterizes DNA-hypermethylated genes in colon cancer cells. Cancer Res. 2008;68(14):5753–9.
Maloof JN et al. A Wnt signaling pathway controls hox gene expression and neuroblast migration in C. elegans. Development. 1999;126(1):37–49.
Breau MA et al. A Hox gene controls lateral line cell migration by regulating chemokine receptor expression downstream of Wnt signaling. Proc Natl Acad Sci USA. 2013;110(42):16892–7.
Calvo R et al. Altered HOX and WNT7A expression in human lung cancer. Proc Natl Acad Sci USA. 2000;97(23):12776–81.
Ryan JF, Baxevanis AD. Hox, Wnt, and the evolution of the primary body axis: insights from the early-divergent phyla. Biol Direct. 2007;2:37.
Yan W et al. Histone deacetylase inhibitors suppress mutant p53 transcription via histone deacetylase 8. Oncogene. 2013;32(5):599–609.
Licchesi JD et al. Epigenetic alteration of Wnt pathway antagonists in progressive glandular neoplasia of the lung. Carcinogenesis. 2008;29(5):895–904.
Fang HL, Yu ZC, Zhu HB, Jin YT. Effects of 5-Aza-2-deoxycytidine on DNA methylation of anti-oncogenes in non-small cell lung cancer cells. Zhonghua Zhong Liu Za Zhi. 2012;34(9):658–63.
Fukui T et al. Transcriptional silencing of secreted frizzled related protein 1 (SFRP 1) by promoter hypermethylation in non-small-cell lung cancer. Oncogene. 2005;24(41):6323–7.
Zhang YW et al. Transcriptional inactivation of secreted frizzled-related protein 1 by promoter hypermethylation as a potential biomarker for non-small cell lung cancer. Neoplasma. 2010;57(3):228–33.
Zhang Y et al. Methylation of multiple genes as a candidate biomarker in non-small cell lung cancer. Cancer lett. 2011;303(1):21–8.
Satoh W et al. Sfrp1, Sfrp2, and Sfrp5 regulate the Wnt/beta-catenin and the planar cell polarity pathways during early trunk formation in mouse. Genesis. 2008;46(2):92–103.
Lee CH et al. Loss of SFRP1 expression is associated with aberrant beta-catenin distribution and tumor progression in mucoepidermoid carcinoma of salivary glands. Ann Surg Oncol. 2010;17(8):2237–46.
Tang YA et al. A novel histone deacetylase inhibitor exhibits antitumor activity via apoptosis induction, F-actin disruption and gene acetylation in lung cancer. PLoS ONE. 2010;5(9):12417.
Janda CY et al. Structural basis of Wnt recognition by Frizzled. Science. 2012;337(6090):59–64.
Hoang BH, Thomas JT, Abdul-Karim FW, Correia KM, Conlon RA, Luyten FP, Ballock RT. Expression pattern of two Frizzled-related genes, Frzb-1 and Sfrp-1, during mouse embryogenesis suggests a role for modulating action of Wnt family members. Dev Dyn. 1998;212(3):364–72.
Bafico A et al. Interaction of frizzled related protein with Wnt ligands and the frizzled receptor suggests alternative mechanisms for FRP inhibition of Wnt signaling. J Biol Chem. 1999;274(23):16180–7.
Surana R et al. Secreted frizzled related proteins: Implications in cancers. Biochim Biophys Acta. 2014;1845(1):53–65.
Kallberg M et al. Template-based protein structure modeling using the RaptorX web server. Nat Protoc. 2012;7(8):1511–22.
Pierce BG et al. ZDOCK server: interactive docking prediction of protein-protein complexes and symmetric multimers. Bioinformatics. 2014;30(12):1771–3.
Mashiach E, Nussinov R, Wolfson HJ. FiberDock: a web server for flexible induced-fit backbone refinement in molecular docking. Nucleic Acids Res. 2010;38(Web Server issue):457–61.
Pronk S et al. GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit. Bioinformatics. 2013;29(7):845–54.
R Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing. Vienna: R Foundation for Statistical Computing; 2014. http://www.R-project.org/. Accessed 1 Aug 2015.
Lemkul JA. Lysozyme in Water. http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/lysozyme/index.html. Accessed 1 Aug 2015.
Publication of this article has been funded by KAKENHI 23300357 and 2610528 and a Chuo University Joint Research Grant. The authors confirm that the funder had no influence over the study design, content of the article, or selection of this journal.
This article has been published as part of BMC Medical Genomics Volume 9 Supplement 1, 2016. Selected articles from the 5th Translational Bioinformatics Conference (TBC 2015): medical genomics. The full contents of the supplement are available online https://bmcmedgenomics.biomedcentral.com/articles/supplements/volume-9-supplement-1
Availability of data and materials
All data analyzed in this manuscript was taken from GEO.
YHT, MI and UM designed researches and wrote the paper. YHT performed analyses and numerical simulations. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
c jk of PC3 and PC4 employed for feature extraction. The left column corresponds to gene expression and the right column corresponds to promoter methylation. PC3 and PC4 show distinct sample dependence. PC3 represents sample dependence that distinguishes between two NSCLC cell lines and non-NSCLC cell lines, while PC4 represents the distinction between two NSCLC cell lines in addition to that between two NSCLC cell lines and non-NSCLC cell lines. (PDF 7 kb)
Hierarchical clustering with 23 samples. Hierarchical clustering between PCs with 23 samples. Left (right) column corresponds to those before (after) re-labelling. (PDF 77 kb)
FE based upon categorical regression. Gene expression/promoter methylation of genes selected by categorical regression. FE based upon categorical regression does not always guarantee a high correlation between gene expression and promoter methylation, because it simply filters those with distinct expression/methylation among samples. However, among eight genes selected, four had a significant (P < 0.05) negative correlation coefficient between gene expression and promoter methylation. This suggested the feasibility of FE based upon categorical regression. The correlation between gene expression and promoter methylation is very high. However, only one gene had a significant positive correlation. (PDF 19 kb)
Genes selected based upon PC3. Gene expression/promoter methylation of genes selected by PCA based unsupervised FE employing PC3. Among eleven selected genes, eight genes had a significant (P < 0.05) negative correlation between gene expression and promoter methylation. Because we did not restrict the selection of genes to those with negative correlations, large numbers of genes with negative correlations demonstrate the feasibility of our methodology. (PDF 26 kb)
Genes selected based upon PC4. Gene expression/promoter methylation of genes selected by PCA based unsupervised FE employing PC4. Among 16 selected genes, 12 genes had a significant (P < 0.05) negative correlation between gene expression and promoter methylation. Because we did not restrict the selection of genes to those with negative correlations, large numbers of genes with negative correlations indicate the feasibility of our methodology. (PDF 36 kb)
Disease associations of genes in Table 1. Disease associations listed with Gendoo server and associated P-values. (PDF 99 kb)
KEGG pathways associated with genes in Table 1. Pathway image files downloaded from KEGG and the html file is linked to these images. (ZIP 963 kb)
About this article
Cite this article
Taguchi, Yh., Iwadate, M. & Umeyama, H. SFRP1 is a possible candidate for epigenetic therapy in non-small cell lung cancer. BMC Med Genomics 9, 28 (2016). https://doi.org/10.1186/s12920-016-0196-3
- Histone Modification
- Promoter Methylation
- Root Mean Square Deviation
- HDAC Inhibitor
- NSCLC Cell Line