Abstract
Potential roles of euchromatic histone methyltransferase 2 (EHMT2 or G9a) in invasion and metastasis are not well understood in non–small cell lung cancer (NSCLC). Here, we investigated the effect and underlying mechanisms of G9a and therapeutic implications of targeting G9a in the invasion and metastasis of NSCLC. Overexpression of G9a significantly enhanced in vitro proliferation and invasion, while knockdown of G9a drastically suppressed in vivo growth and metastasis of A549 and H1299 NSCLC cells. Knockdown or inhibition of G9a significantly decreased the expression of focal adhesion kinase (FAK) protein and activation of FAK pathway. In addition, defactinib, a potent FAK inhibitor, partially abolished the G9a-enhanced invasion in these NSCLC cells. Furthermore, targeting G9a was found to suppress NF-κB transcriptional activity in NSCLC cells through stabilizing NF-κB inhibitor alpha (IκBα), while an NF-κB inhibitor Parthenilide partially abolished the G9a-enhanced FAK activation, which suggests that G9a-enhanced invasion and activation of FAK is mediated by elevated NF-κB activity. Notably, a strong positive correlation between the IHC staining of G9a and phosphorylated FAK proteins was identified in H1299 xenografts and 159 cases of NSCLC tissues (R = 0.408).
The findings of this study strongly demonstrate that G9a may promote invasion and metastasis of NSCLC cells by enhancing FAK signaling pathway via elevating NF-κB transcriptional activity, indicating potential significance and therapeutic implications of these pathways in the invasion and metastasis of NSCLCs that overexpress G9a protein.
Introduction
Lung cancer is a major global health problem and is the leading cause of cancer-related mortality in both men and women (1). Although advances including lung cancer screening, targeted therapies, and immunotherapy have improved outcomes in patients with lung cancer (2), the overall 5-year survival of patients with lung cancer is still just 5% for patients with advanced lung cancer (3). Therefore, novel therapeutic approaches for lung cancer are urgently needed.
G9a histone methyltransferase encoded by the euchromatic histone-lysine N-methyltransferase 2 (EHMT2) is responsible for catalyzing monomethylation and dimethylation of Histone H3 Lysine 9 (H3K9me1 and H3K9me2) and plays an important role in regulating gene expression (4). G9a epigenetically blocks tumor suppressors and activates oncogenes leading to carcinogenesis and cancer cell growth. Previous studies have shown that G9a expression is elevated in many types of human cancers, including pancreatic, bladder, and breast cancers (5–7). Overexpression of G9a is associated with enhanced proliferation and metastasis in many cancer types. However, the underlying mechanisms about how G9a participates in lung cancer metastasis are still poorly understood.
Tumor metastasis is the major cause of cancer-related death, yet precise mechanisms of metastasis remain incompletely understood (8). Focal adhesion kinase (FAK) has been shown to play a central role in metastasis. FAK is a nonreceptor cytoplasmic tyrosine kinase, which is transcriptionally regulated by P53 and NF-κB (9). It contributes to almost every aspect of tumor metastasis, as well as cancer cell adhesion, growth, migration, and invasion (10, 11). Enhanced phosphorylation of FAK at specific sites, especially at Y397, has been reported in a number of cancer types (12). Studies have shown that FAK is upregulated in non–small cell lung cancers (NSCLC; ref. 13) and is associated with an aggressive clinical course (14). In patients with KRAS mutation, FAK is a potential druggable target as it is a downstream effector of KRAS signaling (15). In KRAS-driven lung adenocarcinomas, FAK inhibitors showed potent antitumor effects in Kras-G12V-INK4a/ARF-deficient lung cancers in mice (16). Preclinical data in KRAS-mutated lung cancers showed that inhibition of FAK resulted in sustained DNA damage by suppressing DNA repair mechanisms and enhancing radiation sensitivity (17). In addition, FAK is considered to be a potential therapeutic target to prevent tumor metastasis in a number of solid tumors (18).
Recently, we reported that G9a was frequently overexpressed in NSCLC tissues and that targeting G9a potently suppressed the growth of NSCLC cells, suggesting that G9a may be a therapeutic target for NSCLC (19). In this study, we extended our investigation into the role of G9a in NSCLC cellular migration, invasion, and metastasis. We studied potential mechanisms of G9a-mediated invasion and its impact on FAK signaling pathway in NSCLC. In addition, we explored the effect of targeting G9a on invasive potential and the activation of FAK signaling pathway in NSCLC cells.
Materials and Methods
Cultured cell lines and treated inhibitors
Two human RAS-mutated non–small cell lung adenocarcinoma cell lines were used in this study. Both cell lines were purchased from ATCC. All of these cell lines were authenticated by Integrative Genomics Core by using short-tandem repeat polymorphism analysis. H1299 (p53-Null, NRAS mutation) was cultured in RPMI1640 (ATCC) medium and A549 (p53-wild-type, KRAS mutation) was cultured with DMEM/F12 (Corning) medium. For all cell lines, cell authentic growth medium was supplemented with 10% FBS, 100 U/mL penicillin, and 100 mg/mL streptomycin. Cell lines were tested Mycoplasma-free using colorimetric Mycoplasma detection assay with HEK-Blue-2 cells. Cells used for experiments were within 20 passages from thawing. UNC0638, a selective inhibitor for G9a and GLP methyltransferase, was purchased from Cayman Chemical Company and used at a concentration of 2.5 μmol/L and 5 μmol/L. Defactinib hydrochloride, a FAK inhibitor, which inhibits FAK phosphorylation at Tyr397, was purchased from MedChemExpress and used at a concentration of 1 μmol/L. Parthenilide, an NF-κB inhibitor, was purchased from Abcam and used at a concentration of 20 μmol/L.
Plasmids, cell transfection, and qRT-PCR
G9a mRNA (isoform a) was cloned into pcDNA3.1 vector with Hind III and EcoR I (New England BioLabs). Two independent G9a siRNAs: catalog no. 10630318 (5′-GGCAUCUCAGGAGGAUGCCAAUGAA-3′) and catalog no. 10620319 (5′-UUCAUUGGCAUCCUCCUGAGAUGCC-3′), purchased from Thermo Fisher Scientific Corporation, were used to silence G9a expression by using RNAi Maxi Reagent (Invitrogen) as described previously (20). A FAK siRNA (catalog no.: sc-29310) was purchased from Santa Cruz Biotechnology. Cells were transfected with G9a-pcDNA3.1 plasmid and pcDNA3.1 plasmid as control by using Lipofectamine 3000 Reagent (Invitrogen) according to the manufacturer’s protocol. tRNAs from cells were extracted with RNeasy Mini Kit (catalog no.: 74104) and cDNA was synthesized from 1.0 μg RNA with QuantiTech Reverse Transcription Kit (catalog no.: 205313), which were purchased from Qiagen. EMHT2 (forward: GGACATCATCACTCATGCGGAAA, reverse: GCAAAACCATGTCCAAACCAGG) PTK2 (forward: GCAATGGAGCGAGTATTAAAGGTC, reverse: TGGCCACATGCTTTACTTTGTGAC) GAPDH (forward: GGGAAGGTGAAGGTCGGAGTC, reverse: GAGTTAAAAGCAGCCCTGGTGAC) and NFKBIA (forward: GCCATCATCCATGAAGAAAAG, reverse: GGAGTCTGCAGGTTCTT) primers were used in the qRT-PCR assay to quantitatively measure the mRNA expression. Power SYBR Green PCR Master Mix was purchased from Life Technology Ltd (catalog no.: 4367695). Data were presented as ratio to control to show the relative expression of target, normalized to internal control, and relative to the calibration sample.
Cell proliferation, migration, and invasion assays
For cell proliferation, cells were seeded as 1 × 105 per well in 6-well plates. After 72 hours, cells were harvested and counted with Cello-meter Spectrum (Nexcellom Bioscience). For migration assay, after being transfected with siRNA, plasmids, or being treated with inhibitors for 48 hours, cells were detached from tissue culture plate by using 0.25% trypsin–EDTA solution, and pelleted by centrifugation. The cells were resuspended in serum-free culture medium and seeded as 5 × 105 per well in a 6-well plate. After 24 hours, cells were scratched with the tip of 20 μL pipette and washed twice with PBS. Cells were then cultured with serum-free medium and treated with inhibitors. Images were captured immediately after scratching (0 hours), and then at 24 hours for H1299 and 36 hours for A549. Five locations were imaged for each well. Migration distances were measured by Image J software. Three measurements were performed for each image, and average migratory distance was determined. Cell invasion assay was performed using 8 μm Matrigel Invasion Chambers (Corning). Cells were detached after being transfected or treated with inhibitors for 48 hours as described previously (21). Then, 2.5 × 104 cells were seeded on top of the insert with serum-free medium. Lower chambers were filled with growth medium with 10% FBS. The culture medium was continuously supplemented with respective inhibitor in the inhibitor-treated group. After 22 hours, media containing the remaining cells that did not migrate from the top of the membrane was carefully removed. Cells were fixed by placing the insert in 500 μL 10% formalin for 15 minutes, air-drying for 10 minutes at room temperature, and staining with 1% crystal violet solution for 10 minutes. The fixed inserts were then washed with distilled water. Transwells were imaged by microscope, and cells were counted using Qu path software (22). Data were analyzed and presented as the ratio of the average of total number in treated group divided by the average of total number in the paired control group.
Western blot analysis
Total cellular protein samples were extracted with SDS lysis buffer and heated in 95°C for 7 minutes. Proteins were separated by electrophoresis and transferred onto Polyvinylidenefluoride (PVDF) membranes (Thermo Fisher Scientific). Membrane was blocked with 5% nonfat milk in PBST for 40 minutes at room temperature and blotted with primary and secondary antibodies, respectively. A primary antibody against β-actin (1:3000) was purchased from Santa Cruz Biotechnology. FAK (1:1000), pFAK (Tyr397; 1:1000), IκBα (1:1500), H3K9-me2 (1: 2000) antibodies were purchased from Cell Signaling Technology. Antibody against G9a (1:1000) was purchased from GeneTex, Inc.
Luciferase reporter assay of NF-κB
The pSI-Check2-hRluc-NFkB-firefly reporter plasmid was a gift from Qing Deng, Purdue University, West Lafayette, IN (Addgene plasmid # 106979; http://n2t.net/addgene:106979; RRID: Addgene_106979; ref. 23). For each well, 50 ng NF-κB reporter plasmids were transfected into cells at a concentration of 3 × 103 per well for H1299 and 5 × 103 per well for A549 in 96-well plates using lipofectamine 3000 reagent, with either siRNAs or the overexpression plasmids. For inhibitor experiments, cells were treated with UNC0638 12 hours after transfection with the NF-κB plasmid. At 72 hours posttransfection, cells were washed with PBS and lysed with 40 μL 1× passive lysis buffer. Firefly and Renilla luciferase activity were detected by using Dual-Luciferase Reporter Assay System (Promega). For each group, five replicate experiments were performed. Firefly and Renilla luciferase activities were determined and calculated as described previously (24).
In vivo xenograft and metastasis model
All animal protocols were approved by the institutional animal care and use committee (IACUC, 16005) of City of Hope and performed in the animal facility at City of Hope in accordance with federal, local, and institutional guidelines. H1299 cell line was transfected with G9a-shRNA-GFP plasmid, and transfected cells were selected with puromycin. NOD/SCID/IL2R gamma null mice (NSG, 24–27 g, 8–10 weeks of age, female:male = 1:1) from The Jackson Laboratory were used for xenograft experiment. A suspension of 5 × 106 tumor cells (H1299 shRNA-control and H1299 G9a-KD) in 0.1 mL RPMI1640 was injected into the subcutaneous dorsa of mice at the proximal midline (25). About 3 weeks after injection, mice were euthanized by CO2 inhalation, and tumors were excised and fixed for IHC staining. To investigate the impact of G9a knockdown on in vivo metastasis formation, a suspension of 3 × 106 GFP-expressing parental H1299 cells transfected with a scramble control shRNA (shCtrl) and G9a-attenuated H1299 cells (G9a-KD) in which G9a is stably knocked down by a G9a specific shRNA in 0.1 mL PBS was injected into the tail vein of mice, respectively. About 4 weeks after injection, mice were euthanized by CO2 inhalation; livers and lungs were collected and frozen or fixed immediately. The metastatic nodules of each cancer cell in the frozen liver and lung tissue slides were analyzed and imaged by GFP imaging using a Fluorescence Microscope (Nikon).
IHC analysis
This study was reviewed and approved by the Institutional Review Board (IRB#: 13240) of City of Hope National Medical Center (Duarte, CA). All subjects gave written informed consent. A total of 250 patients with lung adenocarcinoma and squamous carcinoma who underwent surgical resection for curative intent between 2002 and 2014 without preoperative chemotherapy or radiation therapy were included, as previously reported (19). Formaldehyde-fixed paraffin embedded (FFPE) tumor samples were sectioned and stained for G9a, H3K9me2, and pFAK (Tyr397) at the pathology core laboratory, as described previously (26). Antibody against pFAK (Tyr397; 1:2000) was purchased from Invitrogen (catalog no.: 700255). Antibody against G9a (1:300) was purchased from GeneTex (catalog no.: 129153). Antibody against H3K9me2 (1:200) was purchased from Abcam (catalog no.: 5327). G9a and pFAK (Tyr397) IHC staining was scored according to different expression positive percentage and intensity as 0 negative, + minimal, ++ moderate, and +++ strong.
Data analysis
All experiments were performed in duplicates or triplicates and repeated at least two times. For group statistics analysis, data were analyzed for variation and significance using Student t test. All data are shown as mean ± SD. Statistical significance was set at P < 0.05. Pearson’s correlation coefficient was used to measure correlation of G9a and pFAK (Tyr397) gene expression.
Results
G9a regulates PTK2 gene expression in NSCLC cells
Through RNA-Seq analysis (see GSE113493 on NCBI GEO website) and gene set enrichment analysis, we found that knockdown of G9a significantly suppressed gene sets of cell motility and cell adhesion signaling pathways of NSCLC cells, which are critical for cancer invasion and metastasis (19); notably, PTK2 that encodes FAK protein is among these significantly downregulated genes in the gene sets (Fig. 1A), suggesting that PTK2 may represent an important G9a target. To further validate the finding and explore the association between G9a and FAK expression, we first knocked down G9a in H1299 and A549 NSCLC cells by two independent siRNAs. As shown in Fig. 1B, in both A549 and H1299 cell lines, upon knockdown of G9a, the level of H3K9me2 was decreased, and FAK protein was also dramatically decreased in both H1299 and A549 cell lines. qRT-PCR analysis showed that G9a mRNA expression level was significantly silenced in both H1299 and A549 cell lines transfected with specific G9a siRNAs (Fig. 1C; P < 0.01). Simultaneously, mRNA expression of PTK2 gene was downregulated in these two cells (Fig. 1D). These data indicate that FAK expression is associated with G9a expression in NSCLC.
G9a regulates PTK2 gene expression in NSCLC cells. A, RNA-seq data of mRNA levels of EMHT2 (G9a) gene and PTK2 (FAK) gene in G9a-attenuated A549 and H1299 cells by gene-specific siRNAs. B, Western blot analysis shows H3K9me2 and FAK was decreased upon silence of G9a in A549 and H1299 cells. C and D, qRT-PCR analysis of mRNAs. C, EMHT2 gene and D, PTK2 gene in G9a-attenuated A549 and H1299 cells. Data are presented as percentages related to the control cells transfected with scramble siRNAs, and is for averages ± standard division based on data of three independent experiments (**, P < 0.01; ***, P < 0.001, compared with control cells).
G9a regulates PTK2 gene expression in NSCLC cells. A, RNA-seq data of mRNA levels of EMHT2 (G9a) gene and PTK2 (FAK) gene in G9a-attenuated A549 and H1299 cells by gene-specific siRNAs. B, Western blot analysis shows H3K9me2 and FAK was decreased upon silence of G9a in A549 and H1299 cells. C and D, qRT-PCR analysis of mRNAs. C, EMHT2 gene and D, PTK2 gene in G9a-attenuated A549 and H1299 cells. Data are presented as percentages related to the control cells transfected with scramble siRNAs, and is for averages ± standard division based on data of three independent experiments (**, P < 0.01; ***, P < 0.001, compared with control cells).
Knockdown and inhibition of G9a suppresses cell migratory and invasive potential of NSCLC cells
To investigate the potential roles of G9a in cell migration and invasion of NSCLCs, cancer cells were first transfected with two different G9a-specific siRNAs for in vitro migration and invasion assays. We observed that cell proliferation was also significantly suppressed (P < 0.01) upon G9a knockdown in these two lung cancer cells (Supplementary Fig. S1A and S1B). As shown in Fig. 2A, compared with the control group, cell migratory potential was significantly suppressed upon knockdown of G9a in both H1299 (top panel) and A549 (bottom panel) cell lines. A significantly decreased migration distance (P < 0.001) was observed in both H1299 and A549 cells upon G9a knockdown (Fig. 2B).
Knockdown and inhibition of G9a suppresses cell migratory and invasive potential of NSCLC cells. A, Representative images for scratch assays of A549 and H1299 cells transfected with control siRNA (Ctrl siRNA) and G9 siRNA 1 and 2 (G9a si1/2; scale bar, 400 μm). B, Quantitative analysis of migration distances in G9a-attenuated NSCLC cells. C, Representative images for Matrigel transwell invasion assays of A549 and H1299 cells transfected with control siRNA (Ctrl siRNA) and G9a siRNA 1 and 2 (G9a si1/2), (scale bar, 200 μm). D, Quantitative analysis of invasive cells in G9a-attenuated NSCLC cells. E, Western blots shows H3K9me2 was decreased by G9a inhibitor UNC0638 in A549 and H1299 cells. F, Quantitative analysis of decreased cell migration (left panel) and invasion (right panel) in A549 and H1299 cells treated with UNC0638. Data are normalized to the untreated control cells and are presented as percentages related to the control cells (*, P < 0.05; **, P < 0.01; ***, P < 0.001, compared with control cells).
Knockdown and inhibition of G9a suppresses cell migratory and invasive potential of NSCLC cells. A, Representative images for scratch assays of A549 and H1299 cells transfected with control siRNA (Ctrl siRNA) and G9 siRNA 1 and 2 (G9a si1/2; scale bar, 400 μm). B, Quantitative analysis of migration distances in G9a-attenuated NSCLC cells. C, Representative images for Matrigel transwell invasion assays of A549 and H1299 cells transfected with control siRNA (Ctrl siRNA) and G9a siRNA 1 and 2 (G9a si1/2), (scale bar, 200 μm). D, Quantitative analysis of invasive cells in G9a-attenuated NSCLC cells. E, Western blots shows H3K9me2 was decreased by G9a inhibitor UNC0638 in A549 and H1299 cells. F, Quantitative analysis of decreased cell migration (left panel) and invasion (right panel) in A549 and H1299 cells treated with UNC0638. Data are normalized to the untreated control cells and are presented as percentages related to the control cells (*, P < 0.05; **, P < 0.01; ***, P < 0.001, compared with control cells).
Invasion assay was carried out in G9a-attenuated cells. As shown in Fig. 2C, cells transfected with G9a siRNA showed lower invasive potential compared with cell transfected with the control siRNA in both H1299 (top panel) and A549 (bottom panel) cells. Statistical analysis showed that compared with the controls, the invasive potential of the G9a-attenuated groups decreased significantly (P < 0.05) in the two NSCLC cells (Fig. 2D).
To examine whether pharmacologic inhibition on G9a activity will also suppress the invasive potential of lung cancer, UNC0638, a selected G9a methyltransferase inhibitor (27), was used to suppress the G9a activity. As shown in Fig. 2E, a drastic decrease of H3K9me2 protein was found in cancer cells treated with UNC0638, indicating the methyltransferase activity of G9a was inhibited significantly. Similar to the data of G9a knockdown experiment, cell migration in these two cells was also suppressed by UNC0638 treatment (Supplementary Fig. S2A). Compared with the control group, a significant decrease was observed in UNC0638-treated cells (Fig. 2F, left panel) in these two cells (P < 0.05). Besides, the number of invasive cells in A549 and H1299 cell lines was also significantly decreased by UNC0638 treatment (Supplementary Fig. S2B). Quantitative analysis of invasive data was presented in Fig. 2F, right panel (P < 0.05). Therefore, above data suggest G9a plays an important role in migratory and invasive potential of NSCLC cell lines.
Inhibition of G9a suppresses the activation of FAK signaling pathway in NCSLC
Considering the critical role of FAK in cancer migration and invasion, we hypothesized that G9a may regulate cell invasion and migration through FAK signaling pathway. To investigate this underlying mechanism, total protein was extracted from cells either transfected with G9a siRNA or treated with UNC0638. Western blot analysis showed that, compared with control siRNA group, FAK and phosphorylation of FAK (pFAK at Tyr397, an autophosphorylation site on the activated FAK that is used as an indicator for FAK activation; Fig. 3A) were drastically decreased in both A549 and H1299 cells upon G9a knockdown. Similarly, after being treated with UNC0638 for 72 hours, total FAK and phospho-FAK proteins were decreased in the both two cell lines (Fig. 3B). In contrast, overexpression of G9a significantly elevated the levels of FAK and pFAK (Tyr397; Fig. 3C). These data indicated that G9a positively regulates the expression of PTK2 gene and activation of FAK signaling pathway in NSCLC cells.
Inhibition of G9a suppresses the activation of FAK signaling pathway in NCSLC. A, Western blot analysis of G9a, total FAK, and pFAK proteins in G9a-attenuated NSCLC cells. B, Western blot analysis of decreased G9a, H3K9me2, total FAK, and pFAK proteins in the cells treated with UNC0638. C, Western blot analysis of G9a, total FAK, and pFAK proteins in A549 and H1299 cells, and the relative protein levels in the cancer cells overexpressing G9a were presented as the normalized folds (under each band) to the control cells transfected with empty plasmid (set as 1).
Inhibition of G9a suppresses the activation of FAK signaling pathway in NCSLC. A, Western blot analysis of G9a, total FAK, and pFAK proteins in G9a-attenuated NSCLC cells. B, Western blot analysis of decreased G9a, H3K9me2, total FAK, and pFAK proteins in the cells treated with UNC0638. C, Western blot analysis of G9a, total FAK, and pFAK proteins in A549 and H1299 cells, and the relative protein levels in the cancer cells overexpressing G9a were presented as the normalized folds (under each band) to the control cells transfected with empty plasmid (set as 1).
Overexpression of G9a enhances cell migratory and invasive potential of NCSLC cells
To further validate the role of G9a in migration, invasion, and activation of FAK signaling pathway of NSCLC cells, we inserted G9a gene into the pcDNA 3.1 vector and examined whether G9a overexpression will enhance the invasion and activation of FAK signaling pathway. Consistently, compared with control cells, we observed overexpression of G9a significantly increased lung cancer cell proliferation (Supplementary Fig. S3A and S3B; P < 0.05). Meanwhile, migration of NSCLC cells (Fig. 4A) was significantly increased with overexpressed G9a protein (P < 0.001; Fig. 4B and C). Furthermore, compared with control group, the invasion of NSCLC cells was also increased significantly with overexpression of G9a protein (Fig. 4D and E; P < 0.05). These data further demonstrated a key role of G9a in cell invasion, migration, as well as activation of FAK signaling pathway in NSCLC.
Overexpression of G9a enhances cell migratory and invasive potential of NCSLC cells. A, Western blot analysis of overexpressed G9a in H1299 and A549 cells. B, Representative images for the scratch assays of H1299 cells overexpressing G9a and A549 cells overexpressing G9a (at 0 hours and 24 hours for H1299/36 hours for A549). C, Quantitative analysis of migration distances in the cells overexpressing G9a. D, Representative images for Matrigel transwell invasion assays of G9a-overexpressed A549 and H1299 cells. E, Quantitative analysis for invasive cells in G9a-overexpressed NSCLC cell lines. Data are presented as percentage of the control cells (*, P < 0.05; **, P < 0.01; ***, P < 0.001, compared with the control cells).
Overexpression of G9a enhances cell migratory and invasive potential of NCSLC cells. A, Western blot analysis of overexpressed G9a in H1299 and A549 cells. B, Representative images for the scratch assays of H1299 cells overexpressing G9a and A549 cells overexpressing G9a (at 0 hours and 24 hours for H1299/36 hours for A549). C, Quantitative analysis of migration distances in the cells overexpressing G9a. D, Representative images for Matrigel transwell invasion assays of G9a-overexpressed A549 and H1299 cells. E, Quantitative analysis for invasive cells in G9a-overexpressed NSCLC cell lines. Data are presented as percentage of the control cells (*, P < 0.05; **, P < 0.01; ***, P < 0.001, compared with the control cells).
G9a regulates cell migration and invasion through FAK signaling pathway
To investigate whether G9a promotes cell invasion and migration directly or indirectly through activating FAK signaling pathway, a rescue experiment was designed and performed. In this experiment, we investigated whether knockdown of FAK by a FAK-specific siRNA can suppress G9a-enhanced FAK activation and cell invasion in NSCLC cells. As shown in Fig. 5A, overexpression of G9a enhanced the phosphorylation of FAK, while knockdown of FAK significantly attenuated or even completely abolished the elevated phosphorylated FAK. While cell migration (Fig. 5B and C) and invasion (Fig. 5D and E) increased after overexpression of G9a, similar to the change in phosphorylated FAK, the enhanced migration and invasion was reversed by knockdown of FAK siRNA in these two cell lines. In addition, the FAK inhibitor also suppressed G9a-enhanced FAK activation and cellular invasion in these lung cancer cells (Supplementary Fig. S4). Therefore, the above data suggest that the increased cell invasion by overexpression of G9a was partially abolished by FAK siRNA as well as FAK inhibitor, indicating that G9a promotes cell invasion and migration via activating of FAK signaling pathway.
FAK knockdown attenuates the enhanced cell migration and invasion by G9a overexpression. A, Western blot analysis showing that activation of FAK signaling pathway was attenuated by knockdown of FAK by siRNA in G9a-overexpressed NSCLC cells. B, Representative images of scratch assays. 0 hours and 24 hours/36 hours for the controls (pcDNA 3.1) and G9a-overexpressed (G9a-pcDNA 3.1) H1299/A549 upon knockdown of FAK by its specific siRNA, respectively. C, Quantitative analysis of migration distances in the control and G9a-overexpressed A549 and H1299 cells after FAK is attenuated by its specific siRNA. D, Representative images for Matrigel transwell invasion assays of the controls (pcDNA 3.1) and G9a-overexpressed (G9a-pcDNA 3.1) H1299 and A549 transfected with control and FAK siRNA. E, Quantitative analysis for invasive cells in the controls and G9a-overexpressed A549 and H1299 cancer cells upon knockdown of FAK. Data are presented as percentages of the control cells transfected with the empty plasmid pcDNA 3.1 (*, P < 0.05; **, P < 0.01; ***, P < 0.001, NS represents no significance; compared with the controls cells).
FAK knockdown attenuates the enhanced cell migration and invasion by G9a overexpression. A, Western blot analysis showing that activation of FAK signaling pathway was attenuated by knockdown of FAK by siRNA in G9a-overexpressed NSCLC cells. B, Representative images of scratch assays. 0 hours and 24 hours/36 hours for the controls (pcDNA 3.1) and G9a-overexpressed (G9a-pcDNA 3.1) H1299/A549 upon knockdown of FAK by its specific siRNA, respectively. C, Quantitative analysis of migration distances in the control and G9a-overexpressed A549 and H1299 cells after FAK is attenuated by its specific siRNA. D, Representative images for Matrigel transwell invasion assays of the controls (pcDNA 3.1) and G9a-overexpressed (G9a-pcDNA 3.1) H1299 and A549 transfected with control and FAK siRNA. E, Quantitative analysis for invasive cells in the controls and G9a-overexpressed A549 and H1299 cancer cells upon knockdown of FAK. Data are presented as percentages of the control cells transfected with the empty plasmid pcDNA 3.1 (*, P < 0.05; **, P < 0.01; ***, P < 0.001, NS represents no significance; compared with the controls cells).
G9a activates FAK signaling pathway by elevating NF-kB transcriptional activity
Studies have already demonstrated that FAK was transcriptionally regulated by P53 and NF-κB transcription factors (9). We observed the regulation of FAK expression by G9a in both p53-wild-type A549 and p53-null H1299 cells. Considering that we hypothesized that G9a might activate FAK expression through NF-κB signaling pathway. Therefore, we investigated the effect of knockdown and overexpression of G9a on transcriptional activity of the NF-κB–controlled luciferase reporter. Dual-luciferase assays showed that silencing of G9a significantly suppressed NF-κB luciferase activity (Fig. 6A; P < 0.05), while overexpression of G9a significantly increased NF-κB activity (Fig. 6B; P < 0.01) in these two NSCLC cell lines. UNC0638 also significantly suppressed NF-κB luciferase activity in both the cell lines (Fig. 6C; P < 0.05). Interestingly, the level of IκBα protein, an inhibitor of NF-κB signaling pathway, was found to be upregulated by knockdown of G9a and downregulated by overexpression of G9a (Fig. 6D), indicating that suppression of G9a on IκBα expression may contribute to the overactivation of NF-κB signaling pathway.
G9a activates FAK signaling pathway by elevating NF-kB transcriptional activity. A, Knockdown of G9a suppressed NF-κB luciferase reporter activity in H1299 and A549 cells (Ctrl siR for control siRNA, G9a si1/2 for G9a siRNA1/2; *, P < 0.05; **, P < 0.01, compared with the control siRNA). B, Increased NF-κB luciferase reporter activity in G9a-overexpressed (G9+) A549 and H1299 cells; (**, P < 0.01, compared with the control cells transfected with pcDNA 3.1). C, NF-κB luciferase reporter activity was suppressed in these cancer cells treated with UNC0638 (*, P < 0.05; **, P < 0.01, compared with the control cells treated with DMSO). D, Western blot analysis of upregulated IκBα (top panel) in G9a-attenuated cancer cells and decreased IκBα (bottom panel) in G9-overexpressed cancer cells. E, Western blots shows NF-κB inhibitor Parthenilide partially abolished the elevated pFAK (Tyr397) by overexpression of G9a. F, qRT-PCR analysis of EMHT2 and NFKBIA mRNAs in G9a-attenuated (left panel) and G9a-overexpressed (right panel) H1299 and A549 cells. G, Western blot analysis of G9a and IκBα proteins in the controls and G9a-attenuated H1299 cells treated with cycloheximide (CHX). Cells were first transfected with either control siRNA or G9a siRNA, at 48 hours posttransfection, cells were then treated with 25 μg/mL CHX. Data are presented as percentages of the control cells (*, P < 0.05; **, P < 0.01; ***, P < 0.001, NS represents no significance; compared with the controls cells).
G9a activates FAK signaling pathway by elevating NF-kB transcriptional activity. A, Knockdown of G9a suppressed NF-κB luciferase reporter activity in H1299 and A549 cells (Ctrl siR for control siRNA, G9a si1/2 for G9a siRNA1/2; *, P < 0.05; **, P < 0.01, compared with the control siRNA). B, Increased NF-κB luciferase reporter activity in G9a-overexpressed (G9+) A549 and H1299 cells; (**, P < 0.01, compared with the control cells transfected with pcDNA 3.1). C, NF-κB luciferase reporter activity was suppressed in these cancer cells treated with UNC0638 (*, P < 0.05; **, P < 0.01, compared with the control cells treated with DMSO). D, Western blot analysis of upregulated IκBα (top panel) in G9a-attenuated cancer cells and decreased IκBα (bottom panel) in G9-overexpressed cancer cells. E, Western blots shows NF-κB inhibitor Parthenilide partially abolished the elevated pFAK (Tyr397) by overexpression of G9a. F, qRT-PCR analysis of EMHT2 and NFKBIA mRNAs in G9a-attenuated (left panel) and G9a-overexpressed (right panel) H1299 and A549 cells. G, Western blot analysis of G9a and IκBα proteins in the controls and G9a-attenuated H1299 cells treated with cycloheximide (CHX). Cells were first transfected with either control siRNA or G9a siRNA, at 48 hours posttransfection, cells were then treated with 25 μg/mL CHX. Data are presented as percentages of the control cells (*, P < 0.05; **, P < 0.01; ***, P < 0.001, NS represents no significance; compared with the controls cells).
To further investigate whether G9a activates FAK signaling pathway through NF-κB, we treated the G9a-overexpressed H1299 cells with a NF-κB inhibitor, Parthenilide. Western blot analysis showed that NF-κB inhibitor could partially abolish the elevated phosphorylation of FAK (Tyr397) that was caused by overexpression of G9a (Fig. 6E).
To further address how G9a regulates IκBα protein. qRT-PCR was used to measure the mRNA expression level of NFKBIA gene, which encodes IκBα protein. Expression of NFKBIA was not significantly modulated by either knockdown of G9a (Fig. 6F, left panel) nor overexpression of G9a in H1299 cells (Fig. 6F, right panel), indicating the regulation may be at the posttranscription level. Next, we investigated whether G9a impacts the stability of IκBα. Cycloheximide was used to suppress protein synthesis to check the stability of IκBα protein. As shown in Fig. 6G, knockdown of G9a stabilized IκBα protein, leading to inhibition of NF-κB signaling pathway. Taken together, these data indicate that G9a knockdown inhibits FAK signaling pathway partially through inhibiting NF-κB signal pathway via stabilizing IκBα protein, although the exact mechanisms for IκBα stabilized by G9a remains to be further explored.
The impact of knockdown of G9a on in vivo growth, metastasis, FAK activation, and the correlation of G9a to pFAK proteins in NSCLC tissues
We first investigated the correlation between G9a expression with in vivo growth and FAK activation in xenografts using stable G9a-KD H1299 cells. As shown in Fig. 7A, the tumor sizes of G9a-KD group were significantly smaller than the control group. Compared with the xenograft tissues of the controls, the IHC intensities of G9a and H3K9me2 were strongly decreased in G9a-attenuated xenograft tissues and were consistent with the in vitro data, and pFAK (Tyr397) levels were downregulated dramatically in G9a-attenuated xenograft tissues as expected (Fig. 7B).
The impact of knockdown of G9a on in vivo growth, metastasis, FAK activation, and the correlation of G9a to pFAK proteins in NSCLC tissues. A, The images (left panel) and weighs (right panel) of xenografts generated by the control cells and G9a-attenuated H1299 (compared with the control xenografts, * P < 0.05. B, Representative images for hematoxylin and eosin (H&E) and G9a, H3K9me2, and pFAK IHC staining in the xenograft tissues of the control and G9a-attenuated H1299 cancer cells (scale bar, 200 μm). C, Left panel, Representative GFP images for liver metastasis (the top of left panel) and lung metastasis (the bottom of left panel) formed by the GFP-labeled control and G9a-attenuated H1299 cells. Right panel, Statistical analysis of the metastasized cells and the number of nodules in liver (the top of right panel) and lung (the right of bottom panel; ** P < 0.01, compared with the control cells). D and E, Representative images for G9a (E) proteins and pFAK (F) proteins scored as 0, 1+, 2+, 3+ in NSCLC tissues (Scale bar, 200 μm). F, Cross-distribution and correlation of G9a IHC scores to pFAK IHC scores in 159 NSCLC tissues. X-axis is for G9a IHC scores, y-axis is for pFAK scores. Number in dot is the sample size with various G9a and pFAK. R is the Spearman's correlation coefficient for G9a and pFAK. G, Schematic diagram of the potential mechanisms for G9a-enhanced metastasis, modulated FAK and NF-κB signaling pathways, and related therapeutic implications in NSCLC.
The impact of knockdown of G9a on in vivo growth, metastasis, FAK activation, and the correlation of G9a to pFAK proteins in NSCLC tissues. A, The images (left panel) and weighs (right panel) of xenografts generated by the control cells and G9a-attenuated H1299 (compared with the control xenografts, * P < 0.05. B, Representative images for hematoxylin and eosin (H&E) and G9a, H3K9me2, and pFAK IHC staining in the xenograft tissues of the control and G9a-attenuated H1299 cancer cells (scale bar, 200 μm). C, Left panel, Representative GFP images for liver metastasis (the top of left panel) and lung metastasis (the bottom of left panel) formed by the GFP-labeled control and G9a-attenuated H1299 cells. Right panel, Statistical analysis of the metastasized cells and the number of nodules in liver (the top of right panel) and lung (the right of bottom panel; ** P < 0.01, compared with the control cells). D and E, Representative images for G9a (E) proteins and pFAK (F) proteins scored as 0, 1+, 2+, 3+ in NSCLC tissues (Scale bar, 200 μm). F, Cross-distribution and correlation of G9a IHC scores to pFAK IHC scores in 159 NSCLC tissues. X-axis is for G9a IHC scores, y-axis is for pFAK scores. Number in dot is the sample size with various G9a and pFAK. R is the Spearman's correlation coefficient for G9a and pFAK. G, Schematic diagram of the potential mechanisms for G9a-enhanced metastasis, modulated FAK and NF-κB signaling pathways, and related therapeutic implications in NSCLC.
To further examine the potential role of G9a in metastasis formation, parental and G9a-attenuated H1299 cells were injected into the tail veins of mice to establish in vivo metastasis. We observed that, at about 4 weeks of postinjection, H1299 cells developed visible metastasis in both liver (Fig. 7C, the left of top panel) and lung (Fig. 7C, the left of bottom panel). Metastatic tumors were more numerous and larger in livers than in lungs. Compared with control H1299 cells, G9a-KD H1299 formed fewer and smaller nodules in both liver tissues (Fig. 7C, the right of top panel; Supplementary Fig. S5A) and lung tissues (Fig. 7C, the right of bottom panel; Supplementary Fig. S5B). Similarly, pFAK (Tyr397) levels were also dramatically decreased in metastasis of G9a-attenuated H1299 cells (data not shown), suggesting the activation of FAK signaling pathway was suppressed upon G9a knockdown.
Furthermore, we analyzed the correlation between G9a and pFAK expression after IHC staining of pFAK (Tyr397). IHC analysis demonstrated nuclear G9a staining (Fig. 7D), and quantitative analysis of G9a on IHC slides was conducted as previously reported (19). IHC analysis showed cytoplasmic pFAK (Tyr397) staining in tumors. pFAK (Tyr397) IHC staining in the same tissue arrays was quantitatively scored using the same scoring method, and representative images of IHC scoring are shown in Fig. 7E. Spearman correlation analysis was used to examine the correlation between G9a and pFAK (Tyr397) IHC staining in these tumor tissues, and the analysis revealed that pFAK (Tyr397) staining was significantly correlated with G9a staining (Fig. 7F; R = 0.408; P < 0.001).
These data indicated that overexpression of G9a may enhance activation of FAK signaling pathway and then promote invasion and metastasis in NSCLC.
Discussion
This study demonstrated for the first time that G9a may promote invasion and metastasis of NSCLC by activating FAK signaling pathway. Our data suggest that in NSCLC, G9a promotes the degradation of IκBα protein, and then may activate NF-κB signaling pathway, which then promotes FAK expression and enhances cancer invasion and metastasis. Given the important role of both G9a and FAK in cancer cell growth and metastasis, this study underscores the potential role of G9a inhibition as a therapeutic target for lung cancer.
Functionally, G9a is responsible for mono or dimethylation of H3K9 and its activity leads to silencing of a number of tumor suppressor genes including CDH1, DUSP5, and SPRY4, etc. (28). Overexpression of G9a is closely associated with tumorigenesis and metastasis in many cancers. Recently, G9a has been demonstrated to promote tumor metastasis by upregulating ITGB3 in gastric cancer (29). It is also responsible for liver cancer development and metastasis by silencing tumor suppressor RARRES2 (30). In lung cancer, we previously reported that G9a promotes tumor growth via enhancing Wnt signaling pathway and HP1α and silencing APC2 expression (19). A previous study demonstrated that G9a promotes lung cancer metastasis via silencing Ep-CAM (31) and CASP1 (32). G9a inhibition significantly suppressed the migration and invasion enhanced by overexpression of SNAIL2 in lung cancer cells (33). In contrast, a recent study demonstrated that depletion of G9a during tumorigenesis enriches tumors in tumor-propagating cells and accelerates disease progression metastasis (34). Recently, we also found that G9a is involved in lung cancer stemness through epigenetic mechanisms of maintaining DNA methylation of multiple lung cancer stem cell genes and their expression (35). Consistent with findings of many other studies, we also found in this study that G9a promotes cancer invasion and metastasis through a novel mechanism, and targeting G9a suppressed invasion and metastasis in lung cancer. Studies revealed that G9a inhibitor UNC0638 reduced growth and metastasis of breast cancer and significantly suppressed epithelial-mesenchymal transition (EMT)-mediated cellular migration and invasion (36). Another G9a methyltransferase activity inhibitor, UNC0646, was shown to suppress hepatocellular tumor growth but activate cellular apoptosis and metastasis (37).
One of the novel findings of the study is that FAK signaling pathway may represent a downstream target of G9a. FAK plays a pivotal role in transducing signals from the plasma membrane to the nucleus, and phosphorylation of FAK promotes cancer cell growth, invasion, and angiogenesis (38). Studies have demonstrated that high level of pFAK (Tyr397) is associated with aggressive behavior of lung cancer, including rapid growth, frequent metastasis, and poor overall survival (39). Several phase I and phase II clinical trials are currently undergoing for FAK inhibitors in cancers (40). And we further showed that FAK inhibition can significantly attenuate G9a-enhanced invasion. Considering the advanced improvements on FAK inhibitor, FAK inhibition may represent a potential therapy for patients with lung cancer with overexpressed G9a.
Furthermore, this study has demonstrated for the first time that G9a may regulate FAK through NF-κB signaling pathway. In the study, we have found that G9a activates NF-κB signaling pathway in lung cancer. NF-κB activity not only promotes tumor cells proliferation, suppresses apoptosis, and attracts angiogenesis, but it also induces EMT, which facilitates distant metastasis (41). In lung cancer cells expressing oncogenic KRAS, NF-κB signaling pathway has been shown to be activated in a KRAS–dependent manner and represents a potential therapeutic target for KRAS–induced lung cancer (42). In many cancers, the primary regulation of the NF-κB pathway is through the association of NF-κB complexes with IκB proteins. IκBα is a major inhibitor of NF-κB signaling pathway that forms a complex with NF-kB subunits p50 and p65, with the degradation of IκBα, P50 and P65 translocate to the nucleus to regulate downstream gene expression (41). For example, a previous study demonstrated that Tetraspanin 15 (TSPAN15) promoted esophageal cancer metastasis via activating NF-κB through promoting the degradation of IκBα (43). Previous study reported that inhibition of G9a epigenetically increases activation of NF-κB signaling pathway in breast cancer (44). In this study, we revealed that the elevated expression of G9a promoted degradation of IκBα protein, and then activated NF-κB. NF-κB signaling pathway has previously been shown to be associated with FAK activation and increased cell proliferation in lung cancer (45). Therefore, although the exact underlying association between NF-κB and FAK signaling pathways should be further explored, our data have indicated that the NF-κB signaling pathway may positively contribute to the activation of FAK signaling pathway and cancer metastasis. We have also demonstrated that NF-κB inhibition suppresses G9a-enhanced invasion, indicating it may also represent a potential druggable target for lung cancer overexpressing G9a protein.
In our previous study, we found G9a activated Wnt signaling pathway. Wnt signaling usually negatively regulates NF-κB signaling pathway in liver, colon, and breast cancer (46–48). But such correlation was not found in head and neck cancer (49), suggesting a context-specific mechanism for the interaction between two signaling pathways. In lung cancer, a study revealed that overexpression of ING4 gene induced growth inhibition of A549 cells by inactivation of Wnt signaling pathway partially through suppressing NF-κB signaling pathway (50), indicating a potential positive interaction between the two pathways. In our studies, we have found that G9a may enhance activation of Wnt and NF-κB signaling pathways in NSCLC cells. And we have also shown overactivation of Wnt signaling pathway does not affect NF-κB signaling pathway in NSCLC cells (Supplementary Fig. S6).
To summarize, we have revealed a novel mechanism by which G9a mediates invasion and metastasis in NSCLC. We have demonstrated that the role of G9a in lung cancer invasion and metastasis is mediated, in part, by enhanced activation of FAK that appears to be mediated by NF-κB signaling pathway (Fig. 7G), suggesting a potential positive interaction between G9a, NF-κB pathway, and FAK pathway in NSCLC cells.
Authors’ Disclosures
No disclosures were reported.
Authors' Contributions
T. Sun: Data curation, formal analysis, validation, investigation, methodology, writing–original draft. K. Zhang: Conceptualization, data curation, supervision, investigation, methodology, project administration, writing–review and editing. R.P. Pangeni: Resources, investigation, methodology, writing–review and editing. J. Wu: Investigation, methodology. W. Li: Resources, methodology, writing–review and editing. Y. Du: Supervision, project administration, writing–review and editing. Y. Guo: Investigation, methodology. S. Chaurasiya: Investigation, methodology, writing–review and editing. L. Arvanitis: Investigation, writing–review and editing. D.J. Raz: Conceptualization, resources, supervision, funding acquisition, project administration, writing–review and editing.
Acknowledgments
The authors thank the generous support of the Baum Family Foundation in support of this laboratory. The authors thank Michael Lewallen at the pathology core for IHC analysis and Michael Nelson at the light microscope core for transwell imaging and cell counting. This study is supported by the V Foundation (D.J. Raz), the Doris Duke Charitable Foundation (D.J. Raz), and the NCI of the NIH (NIH P30CA33572) using several core facilities.
The costs of publication of this article were defrayed, in part, by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.






