Skip to content

Advertisement

  • Research
  • Open Access

Identification of SUV39H2 as a potential oncogene in lung adenocarcinoma

Contributed equally
Clinical Epigenetics201810:129

https://doi.org/10.1186/s13148-018-0562-4

  • Received: 14 June 2018
  • Accepted: 9 October 2018
  • Published:

Abstract

Background

SUV39H2 (suppressor of variegation 3-9 homolog 2), which introduces H3K9me3 to induce transcriptional repression, has been reported to play critical roles in heterochromatin maintenance, DNA repair, and recently, carcinogenesis. Dysregulation of SUV39H2 expression has been observed in several types of cancers. However, neither the genomic landscape nor the clinical significance of SUV39H2 in lung adenocarcinoma has been probed comprehensively.

Methods

In this research, we conducted bioinformatics analysis to primarily sort out potential genes with dysregulated expressions. After we identified SUV39H2, RNA-seq was performed for a high-throughput evaluation of altered gene expression and dysregulated pathways, followed by a series of validations via RT-qPCR and bioinformatics analyses. Finally, to assess the potential oncogenic role of SUV39H2, we employed the invasion assay and clone formation assay in vitro and tumorigenesis assays in mouse models in vivo.

Results

Through bioinformatics analyses, we found that SUV39H2 underwent a severe upregulation in the tumor tissue, which was also confirmed in the surgically removed tissues. Overexpression of SUV39H2 was mainly associated with its amplification and with shorter patient overall survival. Then, the RNA-seq demonstrated that TPM4, STOM, and OPTN might be affected by the loss of function of SUV39H2. Finally, in vitro and in vivo experiments with SUV39H2 knockdown all suggested a potential role of SUV39H2 in both carcinogenesis and metastasis.

Conclusions

SUV39H2 expression was elevated in lung adenocarcinoma. TPM4, OPTN, and STOM were potentially regulated by SUV39H2. SUV39H2 might be a potential oncogene in lung adenocarcinoma, mediating tumorigenesis and metastasis.

Background

Chromatin modifications account for a major part of epigenetic regulation in mammalian cells, which are mainly composed of various posttranslational modifications (PTMs) [1], among which histone acetylation induces chromatin relaxation and therefore transcriptional activation. However, histone methylation can induce both transcriptional activation and repression. Generally, methylation of histone H3 lysine 9 (H3K9) and H3K27 is correlated with repression, and that of H3K4, with activation. Hence, under the complicated regulation processes through histone lysine methyltransferases (KMTs), lysine methylation has emerged as one of the most important PTMs in the processes from development to disease [2, 3].

One of the major KMT families that introduce transcriptional repression is the SUV39 sub-family, including SUV39H2 (KMT1B), G9a (EHMT2), G9a-like protein1, GLP (EHMT1), SETDB1 (KMT1E), and SETDB2 (KMT1F) [4], among which SUV39H1 and SUV39H2 can preferentially read H3K9me1 via their chromodomain and catalyze H3K9me3 [57]. Moreover, SUV39H1 and SUV39H2 exert mutually compensated expression throughout the embryonic development, implying a functional redundancy between the two enzymes [8]. However, complete knockout of both enzymes resulted in prenatal lethality, accompanied by a universal reduction of H3K9me3 levels, indicating that H3K9 methylation plays dominant roles in development [9]. Additionally, knockout of Suv39h1 and Suv39h2 could lead to abnormally long telomeres resulting from the absence of their enzymatic activities of targeting the heterochromatin [10]. Furthermore, in the dysregulation in cancer, SUV39H1 is reported to play a tumor-suppressor role due to its cell proliferation-suppressing activity [1114]. On one hand, SUV39H2 was reported to be associated with the maintenance of heterochromatin via the introduction of H3K9me3 [15]. SUV39H2 was known to be related to DNA repair by methylating the histone H2AX [16], indicating a tumor-promoting activity. Additionally, it was found to be upregulated in bladder cancer [16], hepatocellular carcinoma [17], and acute lymphoblastic leukemia [18]. Besides, dysregulation of SUV39H2 is also seen in several other diseases. A mutation of SUV39H2 was reported to be responsible for hereditary nasal parakeratosis [19]; this provided us insights into SUV39H2-mediated epigenetic regulation in keratinocytes. Additionally, dysregulation of SUV39H2 mediated by H3K9me3 is associated with autoimmune diabetes [20] and steatohepatitis [21] in mice. The activity of SUV39H2 is modulated in various manners. It was reported that automethylation [22] and alternative splicing [23] were responsible for its specificity and function.

As stated previously, to some degree, SUV39H1 and SUV39H2 possess the opposing functions in cancer, even though they have an overlapping enzymatic activity, which might be due to the specificity of their substrate preference [24]. Identification of the target genes regulated by SUV39H2 might further elucidate the molecular mechanisms underlying its overexpression in tumor tissues. Thus, we evaluated cancer genomic data in lung adenocarcinoma and functional small-interfering RNA (siRNA), to explore the clues of SUV39H2 overexpression and its pathological activity in tumor cell lines.

Methods

Antibodies and reagents

The source of the anti-SUV39H2 and anti-GAPDH antibodies was Abcam, Hong Kong, China. The siRNAs were purchased from Sigma-Aldrich.

Cell culture and transfection

All the cell lines were purchased from ATCC and maintained in the indicated media as instructed. Before culture, 10% fetal bovine serum was added to all media. Cells were cultured in a humidified incubator equilibrated with 5% CO2 at 37 °C. Transfection of siRNA was carried out utilizing the Lipofectamine® RNAiMAX Reagent (Invitrogen, Carlsbad, CA, USA), according to the manufacturer’s instructions.

The Cancer Genome Atlas (TCGA) data for lung cancer

The mutation, mRNA expression, DNA copy number, and clinical data used in the research were obtained from the cBio Cancer Genomics Portal [25, 26]. The copy numbers for SUV39H1 and SUV39H2 were generated from the copy number analysis algorithm Genomic Identification of Significant Targets in Cancer (GISTIC) and classified as copy number per gene; “− 2” stands for a deep loss, “− 1” stands for a heterozygous deletion, “0” stands for neutral or diploid, “1” is a low-level gain, and “2” indicates a high-level amplification. The expression of mRNA was represented by the average reads from the RNA-seq of tumor tissues, which were obtained from TCGA [25, 26]. The value of the Z-scores suggests the relative deviations from the diploid samples. Somatic mutation data were obtained from exosome sequencing [25, 26].

The Lung Adenocarcinoma (TCGA, Provisional) dataset

The Lung Adenocarcinoma (TCGA, Provisional) dataset contains approximately 520 primary lung adenocarcinoma cases with fully detailed clinical follow-up. All data were available without restrictions and limitations. For the expression and mutation analysis, the Entrez-IDs of the indicated genes were selected as the “perfect” evidence of annotation.

The non-small cell lung carcinoma (NSCLC) in female nonsmoker (“GSE19804”) dataset

The “GSE19804” dataset (Additional file 1: Table S1) includes 120 paired tissue samples from 60 NSCLC female non-smoker patients. A detailed description of the dataset can be obtained from the original research [27]. The specimens were analyzed utilizing Affymetrix U133plus2.0 expression arrays. The analysis of the expression of the indicated genes was conducted via Illumina probes as the annotation augments.

RT-PCR and quantitative RT-PCR (qPCR)

Total cellular RNA was obtained utilizing the Trizol reagent, following the manufacturer’s instructions (Invitrogen). cDNA was synthesized using the MMLV Reverse Transcriptase (Promega). Relative quantitation was determined utilizing the ABI PRISM 7500 sequence detection system (Applied Biosystems, Foster City, CA, USA) through the measurement of real-time SYBR green fluorescence, and the results were analyzed by means of the comparative Ct method (2−ΔΔCt) with GAPDH as an internal control. Each sample with the indicated primers was submitted to the PCR analyses in three replicates, and the experiment was performed in triplicate. The primers used are listed in Additional file 2: Table S2.

ChIP-qPCR

ChIP was performed in MCF-7 cells as described previously [28, 29]. Briefly, cells were cross-linked with 1% formaldehyde, sonicated, pre-cleared, and incubated with 5–10 μg proper antibody, followed by addition of protein A/G Sepharose CL-4B beads. The beads were then washed in buffer of high and low salt concentration, and then, DNA was eluted for qPCR assays.

Cell invasion assay

Transwell chamber filters (Becton Dickinson) were coated with Matrigel. After being transfected with siRNAs (Additional file 3: Table S3), A549 cells were suspended in serum-free 1640 medium, and then, 5 × 104 cells were seeded into the upper chamber in a volume of 500 μl. The chamber was then cultured in a well containing 500 μl of 1640 media with 10% fetal bovine serum at 37 °C for 18 h. Cells on the upper side of the membrane were removed by cotton swabs and those on the other side were stained and counted. Four high-powered fields were counted for each membrane. For each complete experiment, three independent samples in the indicated groups were subjected to analysis, and every experiment was performed in triplicate.

Tumor tissue samples

Tumor tissues were obtained from eight first-time admitted lung adenocarcinoma patients in the Tianjin Medical University Cancer Institute and Hospital (Tianjin, China), after receiving written informed consent. The utilized procedure was approved by the Ethics Committee of the Tianjin Medical University Cancer Institute and Hospital. None of the patients ever received surgery, radiotherapy, chemotherapy, or other medical intervention before the sample collection.

Luciferase assay

The sequence of OPTN promoter was obtained from UCSC and cloned into pGL3-Basic for transfection. 3xFLAG-SUV39H2 and vector were transfected respectively into A549 cells, and pGL3-OPTN and Renilla were co-transfected in both groups of cells simultaneously. Both groups were harvested and examined for luciferase activity according to the manufacturer’s instruction from Promega. Each group of cells contained six independent cell cultures and each experiment was conducted in a triplicate manner.

Xenograft assay

A549 cells stably expressing luciferase were infected with the indicated lentivirus and seeded either subcutaneously in nude athymic BALB/c mice or intravenously in immunocompromised severe combined immunodeficiency (SCID) mice, and the bioluminescence of the neoplasia was measured via the IVIS imaging system (Xenogen). Six mice were randomly assigned into each group.

Statistical analysis

Statistical analysis was conducted via R software (http://www.r-project.org) and SPSS Statistics 22. For the data process of the GSE19084 dataset, we employed the function “affycoretools” in the R software for organizing the expression matrix from the raw data and function “limma” to determine the dysregulated genes between the paired specimens in the empirical Bayes (eBayes) manner. For the TCGA dataset, the expressions and mutations of indicated genes were paired with the clinical data in the R software. The expression sets of genes were obtained from the RNA sequencing dataset (data_RNA_Seq_v2_expression_median), and the copy number variations were from the exosome sequencing data (data_CNA). When we compared the mRNA expression of SUV39H2 and overall patient survival, the samples were divided into lower and higher expression groups based on the cutoff value obtained from the ROC curve in SPSS 22. The univariate and multivariate analyses were also performed using SPSS22 with SUV39H2 expression grouping and other indicators in the TCGA dataset.

Results

Overexpression of SUV39H2 in lung adenocarcinoma

To determine the expression status of SUV39H2 in lung adenocarcinoma, we analyzed the SUV39H2 expression in lung adenocarcinoma (GSE19084). The specimens were divided into the “adjacent normal tissue” and “tumor” groups, followed by statistical analysis to identify the significantly dysregulated genes in the eBayes manner (Additional file 1: Table S1). Among all the identified genes, SUV39H2 and several other classical oncogenes were prominent: FOXA1 [30, 31], PCNA [32, 33], and EZH2 [34]; additionally, we also found OPTN, which is reported to be involved in the autophagy-related apoptosis of lung cancer cells [35], and STOM, which is a membrane protein whose downregulation indicates poor prognosis and metastasis in non-small cell lung cancer [36] and HER2-positive breast cancer [37] patients. Furthermore, the expression of SUV39H1 did not differ in the two groups (Fig. 1a and Additional file 1: Table S1). DNA copy number alterations (CNAs) can result in the activation of oncogenes and silencing of tumor suppressors in human cancers [38], and CNA-derived gene dysregulation is becoming an emerging issue in lung cancer examination and treatment [39, 40]. Therefore, we hypothesized that the overexpression of SUV39H2 originated from the copy number increase of the gene. Then, we identified that the amplification percentages of SUV39H2 in different tumor stages were consistently higher than those of SUV39H1 (Fig. 1b). Moreover, according to the nonparametric test of the data, we found that the expression level of SUV39H2 was correlated with the CNA of the gene (Fig. 1c), which implies that the gain of copy number in tumor cells might be the underlying mechanism that causes SUV39H2 overexpression. Furthermore, we employed the TCGA data for survival analysis to reveal the potential prognostic value underlying the SUV39H2 overexpression in lung cancer tissue. The results indicated that patients with higher SUV39H2 expression tended to suffer from poorer survival (Fig. 1d). We then employed univariate analysis to examine the relationship between the indicated factors and OS (Table 1), which exhibited that the SUV39H2 expression level, TNM staging, T staging, and N staging were univariately correlated with OS. Furthermore, we performed a multivariate analysis, which suggested that a higher SUV39H2 expression level (p = 0.016, HR = 1.564) and later TNM staging (p < 0.001, HR = 1.527) were associated with poorer prognosis.
Fig. 1
Fig. 1

SUV39H2 was upregulated in lung adenocarcinoma. a The raw data of “GSE19084” chip assay was downloaded from GEO datasets and analyzed via R software. The samples were grouped into adjacent normal tissues and tumor tissues in a pair-wise manner, and the expressions of the indicated genes were analyzed in the eBayes manner. The heatmap profiled the mRNA expression of the indicated genes, and the genes with higher expression (p < 0.001) in the specimen are highlighted in red, and those with lower expression are highlighted in blue. b The CNA data of SUV39H1 and SUV39H2 were sorted from the Lung Adenocarcinoma (TCGA, Provisional) dataset of TCGA. The amplification percentages of the indicated genes in each tumor stage were analyzed according to the clinical data. c Different CNA statuses of SUV39H2 were plotted against the corresponding mRNA expressions of the gene, demonstrating a gradual rise as the CNA status leveled up. Correlation of each indicated pair was assessed via nonparametric test. **p < 0.01, *p < 0.05. d Kaplan-Meier plots of overall survival associated with the mRNA expression level of SUV39H2 in TCGA lung adenocarcinomas

Table 1

Univariate and multivariate analysis of variables correlated to OS

 

Univariate analysis

Multivariate analysis

p value

HR

95% CI

p value

HR

95% CI

SUV39H2 expression

0.004

1.632

1.165–1.846

0.016

1.564

1.088–2.248

Gender

0.758

1.047

0.783–1.400

   

Age

0.208

1.208

0.900–1.621

   

Smoking status

0.547

0.919

0.699–1.209

   

TNM staging

< 0.001

1.621

1.424–1.846

< 0.001

1.527

1.327–1.758

T staging

< 0.001

1.562

1.300–1.878

   

N staging

< 0.001

1.721

1.451–2.040

   

SUV39H2 overexpression in lung cancer tissues and lung cancer cell lines potentially influenced the expression of a variety of genes

Since the overexpression of SUV39H2 was identified to be of importance in the online datasets, we then evaluated the expression of SUV39H2 in surgically removed NSCLC tissues via Western blotting, the results of which showed high consistency with those of the bioinformatics analyses (Fig. 2a). Given that SUV39H2 introduces H3K9me3 to the histone and mediates transcriptional repression [4], what the underlying mechanisms or the specific target genes of SUV39H2 could be is more intriguing. Hence, the expression spectrum analysis of SUV39H2 in lung cancer cell lines was performed, and SUV39H2 overexpression was seen in cells of A549, LTEP, and GLC-82, all of which are lung adenocarcinoma cell lines (Fig. 2b). Then, we performed immunohistochemistry in a tissue micro-array (TMA), which was composed of 32 lung adenocarcinoma and 12 normal adjacent tissue samples. Results demonstrated that SUV39H2 expression was significantly elevated in tumor tissue samples (Fig. 2c). Next, we utilized small interfering RNA to specifically knock down SUV39H2 in A549 cells, followed by the validation of knockdown efficiency through RT-qPCR and Western blotting (Fig. 2d, e). Then, the mRNA of the indicated A549 cells was extracted and subjected to next-generation sequencing (NGS); the representative upregulated and downregulated genes were plotted in the heatmap (Fig. 2f), from which we may conclude that the knockdown of SUV39H2 resulted in a significant transcriptome alteration.
Fig. 2
Fig. 2

SUV39H2 overexpression in lung cancer tissues and lung cancer cell lines potentially influenced the expression of a variety of genes. a The globular proteins of eight surgically removed tumor tissues and the adjacent normal tissues were extracted, and the expression of SUV39H2 was assessed via Western blotting, with GAPDH serving as the internal control. b The whole-cell lysates of the indicated cell lines were used in the Western blot analysis to examine the expression of SUV39H2; GAPDH served as the control. c The representative image of SUV39H2 stain in indicated tissue sample of TMA and the percentage of SUV39H2-positive cells in each group are presented. d, e qRT-PCR and Western blotting were utilized to verify the knockdown of SUV39H2 at both the mRNA and protein levels. Each sample with the indicated primers was used for the PCR analyses in three replicates, and the experiment was performed in triplicate. f Heatmap of the representative genes identified in the RNA-seq. The samples were from three independent transfections with three different siRNA sequences

Identification of potential target genes of SUV39H2

A pathway analysis was performed against the NGS data. The results suggested that knockdown of SUV39H2 might have an impact on a series of pathways including metabolism in cancer, TGF-β signaling pathway, Hippo signaling pathway, and pathways in cancer (Fig. 3a). Among all the genes with altered expressions, we identified OPTN and STOM, which were also downregulated in the previous analysis and participated in tumorigenesis. Moreover, we found TPM4 [41, 42] and CCDC80 [43, 44], which were reported to exert a tumor repression activity, to be upregulated due to SUV39H2 knockdown. Then, we conducted RT-qPCR analysis to validate the potential regulation of the four genes by SUV39H2. Results showed that all these four genes were significantly upregulated in SUV39H2-knockdown A549 cells (Fig. 3b). Furthermore, we tested the protein change of STOM and TPM4 in SUV39H2 steadily knockdown A549 cells, and the results indicated the same as that from the RT-qPCR assay (Fig. 3c). Moreover, we collected four tissue samples of surgically removed lung adenocarcinoma and their paired normal adjacent tissue to evaluate the protein level of STOM and TPM4 (Fig. 3d). To further consolidate the potential relationship, we again utilized the TCGA dataset and analyzed the relationship between SUV39H2 and OPTN or STOM via the Spearman correlation analysis. The results show that the expression of SUV39H2 was significantly negatively correlated with OPTN (r = − 0.176, p < 0.001) and STOM (r = − 0.325, p < 0.001) (Fig. 3e, f). As stated previously, OPTN could suppress the growth and tumorigenicity of lung cancer cells [35]. As the negative correlation between OPTN and SUV39H2 was strongly indicated, we employed luciferase assay against the promoter of OPTN, which suggested that overexpression of SUV39H2 significantly repressed the luciferase activity via interacting with the promoter of OPTN (Fig. 3g). Then, we performed chromatin immunoprecipitation (ChIP) to further validate the transcriptional regulation of SUV39H2 against OPTN. The results demonstrated that loss of SUV39H2 not only diminished the enrichment of SUV39H2 but also the enrichment of H3K9me3 on the promoter of OPTN (Fig. 3h), so we primarily identified the transcriptional repression of OPTN mediated by SUV39H2. Finally, we grouped the patients according to the expressions of both SUV39H2 and OPTN as indicated and determined the overall survival (OS) in the two groups (Fig. 3i). The Kaplan-Meier plot demonstrated that patients with a high level of SUV39H2 expression and low level of OPTN expression tended to suffer from shorter OS than patients with a low level of SUV39H2 expression (p = 0.001).
Fig. 3
Fig. 3

Identification of the potential target genes of SUV39H2. a The pathway enrichment assay was conducted according to the KEGG pathway database, and the top 20 pathways in the assay are enlisted in the figure. b mRNA expression of the indicated genes in control or SUV39H2 knockdown A549 cells. The error bars represent the mean ± SD of three independent experiments. *p < 0.05, **p < 0.01 (two-tailed unpaired t test). Each sample with the indicated primers was used for the PCR analyses in three replicates, and the experiment was performed in triplicate. c The protein level of STOM and TPM4 and H3K9me3 level in SUV39H2 knockdown A549 cells. d The protein level of STOM and TPM4 in surgically removed lung adenocarcinoma tissue samples and their paired normal adjacent tissue samples. e, f The correlation analysis of SUV39H2 with OPTN or STOM with regard to mRNA expression. g The OPTN promoter-mediated luciferase activity in indicated cells. Each group was composed of six samples and each experiment was repeated three times. h The ChIP-qPCR assay conducted with indicated antibodies in each group of cells. **p < 0.01. i Kaplan-Meier plots of overall survival associated with mRNA expression levels of SUV39H2 and OPTN in TCGA lung adenocarcinomas

Potential function of SUV39H2 in lung adenocarcinoma cell lines

As shown in the pathway analysis, the cell junction pathway was identified to be dysregulated, and we performed the invasion assay in A549 cells with or without SUV39H2 knockdown. Results showed that the invasion of the SUV39H2 knockdown cells was severely inhibited, suggesting a prominent role of SUV39H2 in promoting A549 cell invasion (Fig. 3a, b). Due to the potential regulation of SUV39H2 over OPTN and the tumor-repressive activity of OPTN, we next conducted clone formation assays in the A549 cells with SUV39H2 knockdown. As demonstrated in Fig. 4c, the clone formation activity of A549 cells with SUV39H2 knockdown was severely impaired. Next, we investigated the role of SUV39H2 in tumor progression in vivo. For this purpose, we performed a bioluminescence assay to measure tumor growth in situ in nude athymic BALB/c mice or in immunocompromised SCID mice. A549 cells stably expressing luciferase were infected with the indicated lentiviruses and seeded either subcutaneously or intravenously in nude athymic BALB/c mice or in immunocompromised SCID mice, respectively, and the bioluminescence of the neoplasia was measured via the IVIS imaging system (Xenogen). As shown in Fig. 4d, the knockdown of SUV39H2 severely impaired the tumor growth in situ. Additionally, in the intravenous injection group, liver metastasis was significantly abrogated in the SUV39H2 knockdown group (Fig. 4e). Consistently, both the in vitro and in vivo assays demonstrated that the loss of function of SUV39H2 resulted in suppressed tumor growth and invasion. The validation of SUV39H2 knockdown is shown in Fig. 4f.
Fig. 4
Fig. 4

Potential function of SUV39H2 in lung adenocarcinoma cell lines. a, b A549 cells were transfected with siControl or siRNA targeting SUV39H2, before being seeded into the transwell chambers. The invaded cells were stained and counted, and the images represent one field under the microscope (× 10 magnification). The error bars represent the mean ± SD of three independent experiments. *p < 0.05, **p < 0.01 (two-tailed unpaired t test). For each complete experiment, three independent samples were subjected to analysis in the indicated group, and every experiment was conducted in triplicate. c A549 cells infected with the indicated lentiviruses were maintained in culture media for 5 days prior to being stained with crystal violet. Representative photos are shown on the left; they were statistically analyzed as shown on the right. For each complete experiment, three independent samples were analyzed in the indicated group, and every experiment was conducted in triplicate. d A549 cells infected with lentiviruses expressing shSCR or shSUV39H2 were inoculated subcutaneously in 6-week-old female nude mice (n = 6). Six tumors were quantified using bioluminescence imaging 1 week after the initial implantation. The error bars indicate the mean ± SD. *p < 0.05, **p < 0.01 (two-tailed t test). e A549 cells infected with lentiviruses expressing shSCR of shSUV39H2 were injected intravenously in 6-week-old female SCID mice (n = 3). The tumors were detected via bioluminescence imaging. f The knockdown efficiency was confirmed by RT-qPCR and Western blotting

Discussion

In this research, we analyzed the chip assay of lung adenocarcinoma samples and identified SUV39H2 as a potential oncogene due to its overexpression in the tumor tissues. Because of its emerging role in cancers [18, 22, 45], the underlying mechanisms are becoming the focus of studies about SUV39H2. Since it plays a role in regulating the high-order chromatin dynamics [8], the transcriptional repression induced by SUV39H2 might be a major molecular basis for its overexpression. Therefore, we mainly focused on the downregulated genes in the RNA-seq results and identified OPTN, STOM, TPM4, and CCDC80 as potential target genes (Fig. 5). Moreover, we failed to identify SUV39H1 as a differentially expressed gene in the chip analysis, indicating that it might have a unique role in carcinogenesis. We also probed the molecular basis that might cause the specific overexpression of SUV39H2 in tumor cells. Compared to SUV39H1, the amplification percentage of SUV39H2 was much higher throughout the various stages of lung adenocarcinoma; guided by which, we found that the CNA and the expression of SUV39H2 shared a positive relation, leading to an explanation of its overexpression [46]. Furthermore, we also uncovered the overexpression of some classical oncogenes (FOXA1, PCNA, EZH2) that were reported to mediate carcinogenesis in lung adenocarcinoma along with SUV39H2.
Fig. 5
Fig. 5

The work flow and experiment design of the research

We selected the most case-abundant dataset [25, 26], and the survival analysis strongly suggested the oncogenic role of SUV39H2. Furthermore, the statistical analysis again verified its correlation with the overall survival. The cell line assay of SUV39H2 expression showed that SUV39H2 was overexpressed to a greater degree in lung adenocarcinoma than in squamous cell carcinoma, which might also explain the fact that we did not find a notable benefit from the low expression of SUV39H2 in squamous cell carcinoma (TCGA, Provisional) (data not shown).

Moreover, after the primary analysis of the online datasets, we managed to further clarify the oncogenic potential of SUV39H2 via RNA-seq. Through pathway enrichment, we found that several important pathways including metabolism in cancer [47], TGF-β signaling [48], and cell junction [49], all of which are involved in tumorigenesis and progression, are druggable targets during cancer treatment. We also identified several potential target genes of SUV39H2 due to their elevated expressions in the context of SUV39H2 knockdown; however, the underlying transcriptional regulation is still unclear and worth probing.

In addition, further ChIP assays of H3K9me3 indicated that the knockdown of SUV39H2 resulted in the impairment of H3K9me3 enrichment on the promoter of the OPTN, implicating that SUV39H2 could still have an impact on the transcriptional regulation of OPTN. However, the detailed molecular mechanism is still unclear. Analogously, SUV39H2 was reported to be recruited on the promoter of SLIT1 and promote the proliferation and metastasis of colorectal cancer [50], suggesting that the epigenetic modulation of SUV39H2 was of importance. Moreover, histone demethylase LSD1 was found to be stabilized by SUV39H2 via methylation [45], which raised another hypothesis that SUV39H2 correlates with LSD1 during the transcription regulation.

Since the utilization of methyltransferase-targeted therapy is becoming increasingly important in cancer treatment [51, 52], this research illustrated the potential oncogenic activity of SUV39H2 in lung adenocarcinoma and the possible target genes, which provided a deeper insight into its role in carcinogenesis. The inhibitors targeting EZH2, which is co-overexpressed with SUV39H2 in lung cancer, are being thoroughly researched [5355]. This research provides strong clues that SUV39H2 can be a potential drug target in lung cancer and a candidate indicator of cancer prognosis.

Conclusions

In this study, we discovered that SUV39H2 is a potential oncogene in lung adenocarcinoma, whose expression was elevated in tumor tissues. SUV39H2 could potentiate the tumorigenesis and invasion of lung adenocarcinoma cells, probably by repressing OPTN and STOM.

Notes

Abbreviations

EHMT1: 

Euchromatic histone lysine methyltransferase 1

EHMT2: 

Euchromatic histone lysine methyltransferase 2

EZH2: 

Enhancer of zeste 2 polycomb repressive complex 2 subunit

FOXA1: 

Forkhead box A1

KMT: 

Histone lysine methyltransferase

NSCLC: 

Non-small cell lung carcinoma

OPTN: 

Optineurin

PCNA: 

Proliferating cell nuclear antigen

SETDB1: 

SET domain bifurcated 1

SETDB2: 

SET domain bifurcated 2

STOM: 

Stomatin

SUV39H2: 

Suppressor of variegation 3-9 homolog 2

TCGA: 

The Cancer Genome Atlas

TPM4: 

Tropomyosin 4

Declarations

Acknowledgements

Not applicable

Funding

This work was supported by a grant from the National Key Technology R&D Program (No. 2015BAI12B12) and National Natural Science Foundation of China (81572265). This work was supported by a grant from Natural Science Foundation of Tianjin (18JCQNJC10400). This work was partially supported by a grant (B1513 to Y. Z.) from the institutional fund of the Tianjin Medical University Cancer Institute and Hospital.

Availability of data and materials

The dataset “GSE19804” analyzed during the current study is available in the GEO dataset repository, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE19804 [27].

Authors’ contributions

YZ carried out the bioinformatics analysis. YZ and BHL conducted the animal experiments. YZ and JW performed the biochemical assays. All authors read and approved the final manuscript.

Ethics approval and consent to participate

The study was approved by the ethical review board of the Tianjin Medical University Cancer Institute and Hospital, Tianjin, China.

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

(1)
Department of Immunology, Tianjin Medical University Cancer Institute and Hospital, Tianjin, 300060, China
(2)
National Clinical Research Center for Cancer, Tianjin, China
(3)
Key Laboratory of Cancer Prevention and Therapy, Tianjin, China
(4)
Tianjin’s Clinical Research Center for Cancer, Tianjin, China
(5)
Key Laboratory of Cancer Immunology and Biotherapy, Tianjin, China

References

  1. Bannister AJ, Kouzarides T. Regulation of chromatin by histone modifications. Cell Res. 2011;21(3):381–95.View ArticleGoogle Scholar
  2. Greer EL, Shi Y. Histone methylation: a dynamic mark in health, disease and inheritance. Nat Rev Genet. 2012;13(5):343–57.View ArticleGoogle Scholar
  3. Zhang Y, Reinberg D. Transcription regulation by histone methylation: interplay between different covalent modifications of the core histone tails. Genes Dev. 2001;15(18):2343–60.View ArticleGoogle Scholar
  4. Rao VK, Pal A, Taneja R. A drive in SUVs: from development to disease. Epigenetics. 2017;12(3):177–86.View ArticleGoogle Scholar
  5. Wang T, Xu C, Liu Y, Fan K, Li Z, Sun X, et al. Crystal structure of the human SUV39H1 chromodomain and its recognition of histone H3K9me2/3. PLoS One. 2012;7(12):e52977.View ArticleGoogle Scholar
  6. Rice JC, Briggs SD, Ueberheide B, Barber CM, Shabanowitz J, Hunt DF, et al. Histone methyltransferases direct different degrees of methylation to define distinct chromatin domains. Mol Cell. 2003;12(6):1591–8.View ArticleGoogle Scholar
  7. Loyola A, Bonaldi T, Roche D, Imhof A, Almouzni G. PTMs on H3 variants before chromatin assembly potentiate their final epigenetic state. Mol Cell. 2006;24(2):309–16.View ArticleGoogle Scholar
  8. O'Carroll D, Scherthan H, Peters AH, Opravil S, Haynes AR, Laible G, et al. Isolation and characterization of Suv39h2, a second histone H3 methyltransferase gene that displays testis-specific expression. Mol Cell Biol. 2000;20(24):9423–33.View ArticleGoogle Scholar
  9. Peters AH, O'Carroll D, Scherthan H, Mechtler K, Sauer S, Schofer C, et al. Loss of the Suv39h histone methyltransferases impairs mammalian heterochromatin and genome stability. Cell. 2001;107(3):323–37.View ArticleGoogle Scholar
  10. Garcia-Cao M, O'Sullivan R, Peters AH, Jenuwein T, Blasco MA. Epigenetic regulation of telomere length in mammalian cells by the Suv39h1 and Suv39h2 histone methyltransferases. Nat Genet. 2004;36(1):94–9.View ArticleGoogle Scholar
  11. Vandel L, Nicolas E, Vaute O, Ferreira R, Ait-Si-Ali S, Trouche D. Transcriptional repression by the retinoblastoma protein through the recruitment of a histone methyltransferase. Mol Cell Biol. 2001;21(19):6484–94.View ArticleGoogle Scholar
  12. Dick FA, Rubin SM. Molecular mechanisms underlying RB protein function. Nat Rev Mol Cell Biol. 2013;14(5):297–306.View ArticleGoogle Scholar
  13. Nielsen SJ, Schneider R, Bauer UM, Bannister AJ, Morrison A, O'Carroll D, et al. Rb targets histone H3 methylation and HP1 to promoters. Nature. 2001;412(6846):561–5.View ArticleGoogle Scholar
  14. Albacker CE, Storer NY, Langdon EM, Dibiase A, Zhou Y, Langenau DM, et al. The histone methyltransferase SUV39H1 suppresses embryonal rhabdomyosarcoma formation in zebrafish. PLoS One. 2013;8(5):e64969.View ArticleGoogle Scholar
  15. Puschendorf M, Terranova R, Boutsma E, Mao X, Isono K, Brykczynska U, et al. PRC1 and Suv39h specify parental asymmetry at constitutive heterochromatin in early mouse embryos. Nat Genet. 2008;40(4):411–20.View ArticleGoogle Scholar
  16. Sone K, Piao L, Nakakido M, Ueda K, Jenuwein T, Nakamura Y, et al. Critical role of lysine 134 methylation on histone H2AX for gamma-H2AX production and DNA repair. Nat Commun. 2014;5:5691.View ArticleGoogle Scholar
  17. Hung SY, Lin HH, Yeh KT, Chang JG. Histone-modifying genes as biomarkers in hepatocellular carcinoma. Int J Clin Exp Pathol. 2014;7(5):2496–507.PubMedPubMed CentralGoogle Scholar
  18. Mutonga M, Tamura K, Malnassy G, Fulton N, de Albuquerque A, Hamamoto R, et al. Targeting suppressor of variegation 3-9 homologue 2 (SUV39H2) in acute lymphoblastic leukemia (ALL). Transl Oncol. 2015;8(5):368–75.View ArticleGoogle Scholar
  19. Jagannathan V, Bannoehr J, Plattet P, Hauswirth R, Drogemuller C, Drogemuller M, et al. A mutation in the SUV39H2 gene in labrador retrievers with hereditary nasal parakeratosis (HNPK) provides insights into the epigenetics of keratinocyte differentiation. PLoS Genet. 2013;9(10):e1003848.View ArticleGoogle Scholar
  20. Liu XY, Li H. Reduced histone H3 lysine 9 methylation contributes to the pathogenesis of latent autoimmune diabetes in adults via regulation of SUV39H2 and KDM4C. Journal of diabetes research. 2017;2017:8365762.PubMedPubMed CentralGoogle Scholar
  21. Fan Z, Li L, Li M, Zhang X, Hao C, Yu L, et al. The histone methyltransferase Suv39h2 contributes to nonalcoholic steatohepatitis in mice. Hepatology. 2017;65(6):1904–19.View ArticleGoogle Scholar
  22. Piao L, Nakakido M, Suzuki T, Dohmae N, Nakamura Y, Hamamoto R. Automethylation of SUV39H2, an oncogenic histone lysine methyltransferase, regulates its binding affinity to substrate proteins. Oncotarget. 2016;7(16):22846–56.View ArticleGoogle Scholar
  23. Mauger O, Klinck R, Chabot B, Muchardt C, Allemand E, Batsche E. Alternative splicing regulates the expression of G9A and SUV39H2 methyltransferases, and dramatically changes SUV39H2 functions. Nucleic Acids Res. 2015;43(3):1869–82.View ArticleGoogle Scholar
  24. Schuhmacher MK, Kudithipudi S, Kusevic D, Weirich S, Jeltsch A. Activity and specificity of the human SUV39H2 protein lysine methyltransferase. Biochim Biophys Acta. 2015;1849(1):55–63.View ArticleGoogle Scholar
  25. Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer discovery. 2012;2(5):401–4.View ArticleGoogle Scholar
  26. Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013;6(269):pl1.View ArticleGoogle Scholar
  27. Lu TP, Tsai MH, Lee JM, Hsu CP, Chen PC, Lin CW, et al. Identification of a novel biomarker, SEMA5A, for non-small cell lung carcinoma in nonsmoking women. Cancer epidemiology, biomarkers & prevention : a publication of the American Association for Cancer Research, cosponsored by the American Society of Preventive Oncology. 2010;19(10):2590–7.View ArticleGoogle Scholar
  28. Wang Y, Zhang H, Chen Y, Sun Y, Yang F, Yu W, et al. LSD1 is a subunit of the NuRD complex and targets the metastasis programs in breast cancer. Cell. 2009;138(4):660–72.View ArticleGoogle Scholar
  29. Yang Y, Liu R, Qiu R, Zheng Y, Huang W, Hu H, et al. CRL4B promotes tumorigenesis by coordinating with SUV39H1/HP1/DNMT3A in DNA methylation-based epigenetic silencing. Oncogene. 2013;34:104–18.View ArticleGoogle Scholar
  30. Zhu X, Li D, Yu F, Jia C, Xie J, Ma Y, et al. miR-194 inhibits the proliferation, invasion, migration, and enhances the chemosensitivity of non-small cell lung cancer cells by targeting forkhead box A1 protein. Oncotarget. 2016;7(11):13139–52.View ArticleGoogle Scholar
  31. Jamal-Hanjani M, Wilson GA, McGranahan N, Birkbak NJ, Watkins TBK, Veeriah S, et al. Tracking the evolution of non-small-cell lung cancer. N Engl J Med. 2017;376(22):2109–21.View ArticleGoogle Scholar
  32. Sotgia F, Lisanti MP. Mitochondrial markers predict survival and progression in non-small cell lung cancer (NSCLC) patients: use as companion diagnostics. Oncotarget. 2017;8:68095–107.PubMedPubMed CentralGoogle Scholar
  33. Wang Y, Chen T, Huang H, Jiang Y, Yang L, Lin Z, et al. miR-363-3p inhibits tumor growth by targeting PCNA in lung adenocarcinoma. Oncotarget. 2017;8(12):20133–44.PubMedPubMed CentralGoogle Scholar
  34. Kim NY, Pyo JS. Clinicopathological significance and prognostic role of EZH2 expression in non-small cell lung cancer. Pathol Res Pract. 2017;213(7):778–82.View ArticleGoogle Scholar
  35. Liu Z, Chen P, Gao H, Gu Y, Yang J, Peng H, et al. Ubiquitylation of autophagy receptor optineurin by HACE1 activates selective autophagy for tumor suppression. Cancer Cell. 2014;26(1):106–20.View ArticleGoogle Scholar
  36. Arkhipova KA, Sheyderman AN, Laktionov KK, Mochalnikova VV, Zborovskaya IB. Simultaneous expression of flotillin-1, flotillin-2, stomatin and caveolin-1 in non-small cell lung cancer and soft tissue sarcomas. BMC Cancer. 2014;14:100.View ArticleGoogle Scholar
  37. Chen CY, Yang CY, Chen YC, Shih CW, Lo SS, Lin CH. Decreased expression of stomatin predicts poor prognosis in HER2-positive breast cancer. BMC Cancer. 2016;16:697.View ArticleGoogle Scholar
  38. Albertson DG. Gene amplification in cancer. Trends in genetics : TIG. 2006;22(8):447–55.View ArticleGoogle Scholar
  39. Tian HX, Zhang XC, Wang Z, Chen JG, Chen SL, Guo WB, et al. Establishment and application of a multiplex genetic mutation-detection method of lung cancer based on MassARRAY platform. Cancer biology & medicine. 2016;13(1):68–76.View ArticleGoogle Scholar
  40. Desai A, Menon SP, Dy GK. Alterations in genes other than EGFR/ALK/ROS1 in non-small cell lung cancer: trials and treatment options. Cancer biology & medicine. 2016;13(1):77–86.View ArticleGoogle Scholar
  41. Yang R, Zheng G, Ren D, Chen C, Zeng C, Lu W, et al. The clinical significance and biological function of tropomyosin 4 in colon cancer. Biomed Pharmacother. 2018;101:1–7.View ArticleGoogle Scholar
  42. Jeong S, Lim S, Schevzov G, Gunning PW, Helfman DM. Loss of Tpm4.1 leads to disruption of cell-cell adhesions and invasive behavior in breast epithelial cells via increased Rac1 signaling. Oncotarget. 2017;8(20):33544–59.View ArticleGoogle Scholar
  43. Ferraro A, Schepis F, Leone V, Federico A, Borbone E, Pallante P, et al. Tumor suppressor role of the CL2/DRO1/CCDC80 gene in thyroid carcinogenesis. J Clin Endocrinol Metab. 2013;98(7):2834–43.View ArticleGoogle Scholar
  44. Leone V, Ferraro A, Schepis F, Federico A, Sepe R, Arra C, et al. The cl2/dro1/ccdc80 null mice develop thyroid and ovarian neoplasias. Cancer Lett. 2015;357(2):535–41.View ArticleGoogle Scholar
  45. Piao L, Suzuki T, Dohmae N, Nakamura Y, Hamamoto R. SUV39H2 methylates and stabilizes LSD1 by inhibiting polyubiquitination in human cancer cells. Oncotarget. 2015;6(19):16939–50.View ArticleGoogle Scholar
  46. Gu DL, Chen YH, Shih JH, Lin CH, Jou YS, Chen CF. Target genes discovery through copy number alteration analysis in human hepatocellular carcinoma. World J Gastroenterol. 2013;19(47):8873–9.View ArticleGoogle Scholar
  47. Tang L, Wei F, Wu Y, He Y, Shi L, Xiong F, et al. Role of metabolism in cancer cell radioresistance and radiosensitization methods. J Exp Clin Cancer Res. 2018;37(1):87.View ArticleGoogle Scholar
  48. Colak S, Ten Dijke P. Targeting TGF-beta signaling in cancer. Trends Cancer. 2017;3(1):56–71.View ArticleGoogle Scholar
  49. Cerutti C, Ridley AJ. Endothelial cell-cell adhesion and signaling. Exp Cell Res. 2017;358(1):31–8.View ArticleGoogle Scholar
  50. Shuai W, Wu J, Chen S, Liu R, Ye Z, Kuang C, et al. SUV39H2 promotes colorectal cancer proliferation and metastasis via tri-methylation of the SLIT1 promoter. Cancer Lett. 2018;422:56–69.View ArticleGoogle Scholar
  51. Katoh M. Mutation spectra of histone methyltransferases with canonical SET domains and EZH2-targeted therapy. Epigenomics. 2016;8(2):285–305.View ArticleGoogle Scholar
  52. Yan W, Herman JG, Guo M. Epigenome-based personalized medicine in human cancer. Epigenomics. 2016;8(1):119–33.View ArticleGoogle Scholar
  53. Zhang H, Qi J, Reyes JM, Li L, Rao PK, Li F, et al. Oncogenic deregulation of EZH2 as an opportunity for targeted therapy in lung cancer. Cancer discovery. 2016;6(9):1006–21.View ArticleGoogle Scholar
  54. Chen YT, Zhu F, Lin WR, Ying RB, Yang YP, Zeng LH. The novel EZH2 inhibitor, GSK126, suppresses cell migration and angiogenesis via down-regulating VEGF-A. Cancer Chemother Pharmacol. 2016;77(4):757–65.View ArticleGoogle Scholar
  55. Frankel AE, Liu X, Minna JD. Developing EZH2-targeted therapy for lung cancer. Cancer discovery. 2016;6(9):949–52.View ArticleGoogle Scholar

Copyright

© The Author(s). 2018

Advertisement