- Open Access
Tri-methylation of H3K79 is decreased in TGF-β1-induced epithelial-to-mesenchymal transition in lung cancer
Clinical Epigenetics volume 9, Article number: 80 (2017)
The epithelial-to-mesenchymal transition (EMT) enables epithelial cancer cells to acquire mesenchymal features and contributes to metastasis and resistance to treatment. This process involves epigenetic reprogramming for gene expression. We explored global histone modifications during TGF-β1-induced EMT in two non-small cell lung cancer (NSCLC) cell lines and tested different epigenetic treatment to modulate or partially reverse EMT.
Loss of classical epithelial markers and gain of mesenchymal markers were verified in A549 and H358 cell lines during TGF-β1-induced EMT. In addition, we noticed increased expression of the axonal guidance protein semaphorin 3C (SEMA3C) and PD-L1 (programmed death-ligand 1) involved in the inhibition of the immune system, suggesting that both SEMA3C and PD-L1 could be the new markers of TGF-β1-induced EMT. H3K79me3 and H2BK120me1 were decreased in A549 and H358 cell lines after a 48-h TGF-β1 treatment, as well as H2BK120ac in A549 cells. However, decreased H3K79me3 was not associated with expression of the histone methyltransferase DOT1L. Furthermore, H3K79me3 was decreased in tumors compared in normal tissues and not associated with cell proliferation. Associations of histone deacetylase inhibitor (SAHA) with DOT1L inhibitors (EPZ5676 or SGC0946) or BET bromodomain inhibitor (PFI-1) were efficient to partially reverse TGF-β1 effects by decreasing expression of PD-L1, SEMA3C, and its receptor neuropilin-2 (NRP2) and by increasing epithelial markers such as E-cadherin.
Histone methylation was modified during EMT, and combination of epigenetic compounds with conventional or targeted chemotherapy might contribute to reduce metastasis and to enhance clinical responses.
The epithelial-to-mesenchymal transition (EMT) is a highly dynamic and reversible mechanism by which epithelial cells can convert into a mesenchymal phenotype, allowing a loss of cellular adhesion, cellular polarity, and an improvement in migratory and invasive properties. This process occurs during embryonic development, wound healing, and metastatic expansion [1,2,3]. It also plays a major role in resistance to cancer treatment [4, 5]. A major inducer of EMT is TGF-β1, along with cytokines and growth factors secreted by the tumor microenvironment. Switch in gene expression during EMT is characterized by repression of epithelial genes and induction of mesenchymal genes.
The induction of EMT is associated with reprogramming of the epigenome characterized by chromatin remodeling, changes in DNA methylation, post-translational histone modifications or insertion of histone variants, and modifications of non-coding RNA expression [2, 6, 7]. The basic unit of chromatin, the nucleosome, is formed by 180–200 bp of DNA wrapped around a histone protein complex, composed by an octamer of two copies of each histone H2A, H2B, H3, and H4, and fixed by histone H1 . Each histone can be affected by post-translational modifications (PTMs), including acetylation, methylation, phosphorylation, ubiquitination or sumoylation, that define a complex “histone code.” PTMs are regulated by histone-modifying enzymes which add or erase these modifications. They are recognized and read by protein partners to control the accessibility of the transcriptional machinery to nearby genes . All these processes can be considered as epigenetic-based therapeutic strategy to treat cancer .
Lung carcinomas cause one fifth of cancer deaths worldwide . Large-scale genomic studies have characterized frequent modifications affecting epigenetic mechanisms [12,13,14]. Recurrent mutations of epigenetic modifying genes affect the SWI/SNF chromatin remodeling components (ARID1a and BRG1) and the H3K36 methyltransferase SETD2 in about 20% of lung adenocarcinomas [15, 16]. The H3K79 methyltransferase DOT1L mutations are less frequent and are described in 3% of lung adenocarcinomas . In squamous lung cancers, mutations of the H3K4 methyltransferase MLL2 are described in 20% . Abnormal expression of epigenetic enzymes includes overexpression of EZH2 (enhancer of zeste homolog 2), the catalytic subunit of the PRC2 repressive complex that methylates H3K27, associated with tumor progression and poor prognosis in lung cancer [19,20,21,22]. Abnormal epigenetic marks are noticed in lung cancer, mainly DNA methylation, histone acetylation, and methylation for the most studied marks . DNA methylation and miRNAs have emerged as potential biomarkers in body fluids for lung cancer [13, 23].
During EMT, the transcription factors SNAIL and ZEB1 recruit several epigenetic players, including histone methyltransferases, the lysine demethylase LSD1, HATs (histone acetyltransferases) and HDACs (histone deacetylases), sirtuins, and BRG1 . We reported a global decrease in H3K27 acetylation in a ZEB1-induced EMT lung cancer cell model and suggested that ZEB1 would recruit EZH2 . EMT is also associated with repression of the miR-200 family. However, EMT-related epigenetic reprogramming is still poorly understood in lung tumors.
In the present study, we focused on epigenetic modifications during TGF-β1-induced EMT in non-small cell lung cancer (NSCLC) cells and tested different epigenetic treatment to modulate or partially reverse EMT.
Cell lines, inhibitors, and antibodies
NSCLC cell lines, A549, and NCI-H358 (hereafter H358) were obtained from ATCC in 2014. The cells were grown in RPMI-1640 medium with 10% fetal bovine serum (FBS) and antibiotics-antimycotics (#15240-062, Invitrogen, Carlsbad, CA, USA), at 37 °C and 5% CO2, and controlled every month for mycoplasma contamination. Recombinant human TGF-β1 was from R&D Systems (Minneapolis, MN, USA). The inhibitors EPZ5676, SGC0946, PFI-1, and SAHA were purchased from Cayman (Ann Arbor, Michigan, USA). The antibodies are listed in Additional file 1.
RNA extraction and quantification by real-time quantitative PCR
Total RNA was extracted using the RNeasy Mini kit (Qiagen, Hilden, Germany) following the manufacturer’s instructions. RNA quality was controlled by electrophoresis on 0.8% agarose gel. Five hundred nanograms of total RNA was reverse-transcribed using the iScript Reverse Transcription Kit (Biorad, Hercules, CA, USA). Real-time quantitative PCR (qPCR) was performed using SsoAdvanced Universal SYBR Green Supermix (Bio-Rad) on a CFX384 system (Bio-Rad). Data were analyzed using the Bio-Rad CFX Manager software and normalized to GAPDH messenger RNA (mRNA) level using the 2−ΔCq method. The primer sequences are listed in Additional file 2.
Protein detection by immunoblotting
Cells were directly lysed in Laemmli loading buffer and sonicated. Whole-cell extracts were separated on SDS-PAGE gels (4–15% acrylamide gradient, #4568084, Bio-Rad) and transferred to 0.2 μm nitrocellulose membrane (#1704159, Bio-Rad) with the Trans-Blot Turbo Transfer System (Bio-Rad) for 7 min. The membranes were blocked with 5% BSA in PBS 0.05% Tween (PBS-T) for 90 min. They were then incubated with primary antibodies at the indicated dilutions (Additional file 1) in 5% BSA in PBS-T, at 4 °C overnight. Removal of excess primary antibodies was carried out by washing the membranes in PBS-T (3 × 5 min each). Secondary antibodies were incubated with the membranes in 5% BSA-PBS-T for 1 h at room temperature. The membranes were washed in PBS-T (3 × 5 min each) before exposition to Clarity Western ECL substrate (#1705060, Bio-Rad). Blot images were acquired with the ChemiDoc MP System (Bio-Rad) and quantified with the Image Lab software (Bio-Rad). Protein normalization was carried out with anti-actin or anti-total H3 histone staining.
Cells grown on 8-well Ibidi plate (Martinsried, Germany) were fixed for 15 min with 4% paraformaldehyde. After rinsing with PBS (3 × 5 min), the cells were permeabilized for 20 min with 0.2% Triton X-100 in PBS. After rinsing (3 × 5 min with PBS), blocking was performed with 3% normal chicken serum (#ab7477, Abcam, Cambrigde, UK) in PBS during 1 h. Primary and secondary antibodies were sequentially applied for 60 min at room temperature in 3% serum in PBS at the indicated dilutions (Additional file 1). DRAQ5 at 1:250 dilution (Thermo Scientific, Rockford, IL, USA) was applied with the secondary antibody for DNA staining. Stained slides were rinsed in PBS and mounted in Dako fluorescent media (Santa Clara, CA, USA). Images were captured with a confocal microscope (Olympus FluoView FV1000) at ×60 immersion oil objective.
5 × 105 cells were requested per condition and were diluted at 5 × 106 cells/ml in PBS 2% FBS. Near-IR staining (#L10119, Thermo Fisher, Waltham, MA, USA) was performed during 20 min at 4 °C to detect dead cells. The cells were washed once in PBS 2% FBS and stained for 1 h at 4 °C with PD-L1 antibody or BB515 isotype control at the indicated dilutions (Additional file 1). BB515 IgG1 isotype control is used as a negative control and binds specifically to KLH antigen which is not expressed in human cells. Flow cytometry analyses were conducted on a BD FACSVerse flow cytometer (BD Biosciences, San Jose, CA, USA). Data were analyzed using the FlowJo software (Tree Star, Ashland, OR, USA).
Cancer tissue microarray (TMA) and immunohistochemistry
Three commercial formalin-fixed, paraffin-embedded TMA slides were obtained from USBiomax (Rockville, MD, USA). LC121a is a lung cancer TMA (120 samples), MC6163 a multi-tissue TMA (616 samples) that includes normal and cancer tissues, and BCN962 a combined multiple normal and cancer tissue microarray (96 samples), with 17 types of common organs. Both TMAs LC121a and MC613 contain 110 and 48 lung cancer tissues, respectively, with evaluation of TNM disease stages. TMA slides were deparaffinized by incubations (2 × 3 min) in Histosol solution (National Diagnostics, Inc., Charlotte, NC, USA) and rehydrated by sequential immersions in 100% ethanol (2 × 1 min) and 70% ethanol (1 min). After washing 2 min in water, the slides were incubated for 2 min in a pressure cooker with the antigen unmasking solution (10 mM citrate buffer, pH = 6) and left for cooling for additional 20 min in the de-pressurized cooker. After washing 5 min in TBS-Tween (Dako Wash Buffer #S3006, Glostrup, Denmark), endogenous peroxidases were inhibited with 3% H2O2-PBS during 5 min at room temperature. Following washing in TBS-Tween (1 × 5 min), the tissues were permeabilized with 0.2% Triton in PBS for 20 min. Blocking was performed with 3% BSA-PBS during 20 min. Both TMAs LC121a and MC613 were incubated for 1 h at room temperature in a wet chamber with rabbit anti-H3K79me3 primary antibody (1:200; #pAb-068-050, Diagenode, Denville, NJ, USA,). For BCN962 slides, rabbit anti-H3K79me3 primary antibody was incubated for 1 h at room temperature or with mouse anti-Ki67 primary antibody (1:100; #M7240, Dako) in Emerald diluent (#936b-08, Cell Marque, Rocklin, CA, USA) overnight at 4 °C. Following washing in TBS-Tween (1 × 5 min), the slides were incubated for 30 min with the HRP-labeled polymer conjugated with secondary antibodies (#K4065, EnVision® + Dual Link System-HRP kit, Dako). Slides were washed 5 min in TBS-Tween and further incubated in 3,3-diaminobenzidine (DAB) for 5 min. After washing in H2O, tissues were incubated in hematoxylin for 2 min. TMA slides were dehydrated by sequential immersions in 70% ethanol (1 min) and 100% ethanol (2 × 1 min) and finally mounted using mount medium (#10046430011, Leica, Wetzlar, Germany). For TMAs LC121a and MC613, each sample was scored by multiplying the percentage of positive (0 to 100%) cells in the tumor compartment by the average level of staining intensity (0 to 3). For TMA BCN962, Ki67 expression was determined by the percentage of positive (0 to 100%) epithelial cells in the tumor compartment and H3K79me3 expression was analyzed by the level of staining (0 to 3) for epithelial cells only.
Data were summarized by median and range for quantitative variables, percentage, and confidence intervals when appropriate for qualitative variables. Percentage of positive cells staining intensity and expression score distributions in tumor, stroma, and normal tissues were compared using Wilcoxon rank-sum test and the Kruskal-Wallis test. Relationship between co-staining status and qualitative parameters (histology, nodes) was analyzed by Fisher’s exact test or chi-square as appropriate. Correlations between score expression and other quantitative parameters were determined with the Spearman rank correlation method. Potential relationships with baseline characteristics were also explored with the use of non-parametric test, as appropriate.
TGF-β1 exposure induces EMT in NSCLC cell lines
A549 and H358 NSCLC cell lines were selected for EMT induction by TGF-β1 because of their previous characterization in a ZEB1-induced EMT model . We first verified target gene expression over time during a 72-h period to validate our cellular models. As expected, we could verify the more epithelial status of H358 cells compared to that of A549 cells. Upon TGF-β1 treatment, the progressive loss of the classical epithelial marker like E-cadherin was noticed in addition to the gain of the mesenchymal markers N-cadherin and vimentin (Fig. 1a, b). Expression of other epithelial markers ST14, ESPR1, EpCAM, and Rab25 were also decreased like in the ZEB1-induced EMT model [24, 25]. ZEB1, one of the transcription factors involved in EMT induction, was increased as expected (Fig. 1a). Since EMT is often partial and NSCLC cell lines display intermediate/hybrid states with mixed epithelial and mesenchymal characteristics , we checked the expression of additional genes coding for the guidance proteins, semaphorin 3C (SEMA3C) and semaphorin 3F (SEMA3F), and one of their common receptors, neuropilin-2 (NRP2), because of their involvement in lung cancer progression and EMT [27, 28]. In both cell lines, SEMA3F expression did not change (result not shown), but SEMA3C and NRP2 were increased (Fig. 1a) suggesting a more aggressive phenotype upon TGF-β1 treatment. From these collected data, a 48-h TGF-β1 treatment was chosen for EMT induction for further analysis.
TGF-β1-induced EMT increases PD-L1 expression in H358 cells
For a better characterization of TGF-β1-induced EMT, we analyzed the programmed death-ligand 1 (PD-L1) expression in A549 and H358 cells. PD-L1 plays a critical role by associating programmed death 1 receptor (PD-1) on tumor-infiltrating T cells thus inhibiting the immune response. Indeed, in 1070 surgically resected NSCLC specimens, PD-L1 was expressed in 44% of them and strong PD-L1 staining correlated with poor prognostic . Interestingly, EMT is associated with an inflammatory tumor microenvironment in lung adenocarcinoma .
Untreated H358 cells were found positive for PD-L1, whereas PD-L1 was not detectable in A549 cells previously described as a negative control in several studies [30, 31] (Fig. 2a, b). After a 48-h TGF-β1 treatment, PD-L1 was still not detected in A549 cells but was increased at the mRNA and protein levels in H358 cells (Fig. 2a, b). Immunofluorescence confirmed increased membranous PD-L1 staining (Fig. 2c), and a 60% increase in intensity was estimated by flow cytometry (Fig. 2d). Therefore, H358 cells would become more aggressive upon TGF-β1 treatment.
H3K79 methylation is decreased in TGF-β1-induced EMT
We investigated global histone modifications during EMT induction in both A549 and H358 cell lines. Of interest, untreated A549 and H358 cells present differences for some of the tested histone marks, such as H3K9me1, -me2, H3K36me1, -me2, H3K79me1, -me2, -me3, and H2BK120ac. Whereas most global histone marks did not change, up to 72-h TGF-β1 treatment, some of them including H3K79me3 and H2BK120me1, decreased in both cell lines (Fig. 3a). Of note, A549 cells have a higher basal level of H3K79me1, -me2, and -me3 than H358 cells. Immunocytochemistry confirmed H3K79me3 decrease in A549 cells after a 48-h treatment, but variation of H3K79me3 in H358 cells was difficult to estimate because of the low basal level (Fig. 3b). The different methylation states of H3K79 are generated in a distributive manner by the histone methyltransferase DOT1L , and the di-methylation of H3K79 is activated by ubiquitination of H2BK120 , suggesting that H2BK120ub1 may associate with DOT1L to promote H3K79 methylation . However, H2BK120ub1 level was not modified in A549 neither in H358 cells during EMT induction, but H2BK120me1 was decreased (Fig. 3a). H2BK120ac was also decreased in A549 cells. This result suggests that H2BK120ub1 global level may not be directly involved in the reduction of H3K79me3 but that other modifications like acetylation and methylation of H2BK120 could be.
Next, we investigated DOT1L expression. After a 48-h treatment, DOT1L expression was decreased at the mRNA and protein levels in both cell lines (Fig. 4a, b). Of note and surprisingly, H358 cells express more DOT1L than A549 cells but H3K79me3 is less, suggesting a complex mechanism for H3K79 methylation.
Partial inhibition of EMT induction by epigenetic compounds
Epigenetic compounds were tested to reverse EMT in untreated cells and to inhibit EMT induced with TGF-β1 (Fig. 5). Results are summarized in Table 1. First, in order to evaluate the involvement of H3K79me3 in TGF-β1-induced EMT, A549 and H358 cells were treated with DOT1L inhibitors, EPZ5676 and SGC0946, simultaneously with TGF-β1 (Fig. 5a). The efficacy of these compounds was verified by H3K79me3 decrease (Fig. 5c). In absence of TGF-β1, these compounds did not modify expression of EMT-related genes except for an increase of ESRP1 with SGC0946 (Fig. 5b). In combination with TGF-β1, both DOT1L inhibitors did not reverse expression induced by TGF-β1 for the selected genes except for PD-L1 that was reduced at the protein level in H358 cells (Fig. 5c).
Histone acetylation is involved in DOT1L activation through H2B ubiquitination  (Fig. 6). Combined treatments with inhibitors of BET family proteins that are readers for acetylation, and inhibitors of histone deacetylases, induced anti-cancer effects in mouse models [36, 37]. For these reasons, we examined the impact of PFI-1 (a BRD2 and BRD4 inhibitor) and SAHA (a HDAC inhibitor). On control H358 cells, SAHA increased E-cadherin, Rab25, and decreased ST14 expression suggesting a partial gain of an epithelial phenotype (Fig. 5b, c). When added to TGF-β1, SAHA increased E-cadherin, Rab25, N-cadherin, and PD-L1 expression but decreased NRP2, SEMA3C, ST14, and ZEB1 expression compared to the TGF-β1 control.
PFI-1 treatment in H358 control cells, by itself, did not induce gross modifications but interestingly decreased PD-L1 and SEMA3C expression. When combined with TGF-β1, the major effect of PFI-1 was to decrease NRP2, PD-L1, SEMA3C, vimentin, and ZEB1 expression (ratio treated/untreated, 0.15, 0.27, 0.03, 0.58, and 0.26, respectively) to levels similar to control cells. However, expression of epithelial genes such as ESRP1 and ST14 was decreased. Next, DOT1L inhibitors or PFI-1 were combined with SAHA in presence or absence of TGF-β1 treatment. Interestingly, all tested combinations were able to reduce NRP2, SEMA3C, and ZEB1 in presence of TGF-β1. The combination between DOT1L inhibitors and SAHA was sufficient to increase E-cadherin expression in H358 cells (Fig. 5b). The addition of SAHA to PFI-1 improved the response compared to PFI-1 alone to increase E-cadherin expression and to decrease NRP2 expression in control cells. This combination did not improve the PFI-1 effect in TGF-β1-treated H358 cells.
In TFG-β1-treated A549 cells, SAHA restored H3K79me3 and induced epithelial gene expression but reduced NRP2 and SEMA3C. In combination with DOT1L inhibitors or PFI-1, NRP2 and ZEB1 expression was reduced (Additional file 3 a, b).
These results show that each treatment has benefits and disadvantages to modulate the EMT status or EMT induction. They suggest that SAHA and both combinations of either SAHA/DOT1L or SAHA/PFI-1 inhibitors are potential epigenetic-based therapies to partially reduce or reverse EMT (Table 1).
H3K79me3 in human normal tissues and cancers
Since H3K79me3 was decreased during EMT in cancer cell lines suggesting a relation with aggressiveness, we investigated the presence of H3K79me3 in lung cancers by immunohistochemistry on two commercial TMAs. Because of sample quality, some tumors were excluded, leaving for analysis 17 squamous cell carcinomas, 32 large carcinomas, 48 adenocarcinomas, 3 papillary adenocarcinomas on TMA LC121a, and 11 squamous cell carcinomas, 10 large carcinomas, 12 adenocarcinomas, and 10 small cell carcinomas on TMA MC6163. Scores were obtained by multiplying the percentage of positive (0 to 100%) cells in the tumor compartment by the average level of staining intensity (0 to 3) (Fig. 7a). Histology, tumor grade, and metastatic lymph node did not affect H3K79me3 scores. However, H3K79me3 intensity was statistically higher in large cell carcinomas than in adenocarcinomas in TMA LC121a (median = 2, range 0–3 versus median = 1, range 0–3, respectively, p = 0.0149), but the percentage of stained cells was not significantly different suggesting more H3K79me3 in large cell carcinomas. Because of the low number of corresponding samples, this difference was not found on TMA MC6163. In addition, a tendency for less staining was observed for small cell carcinomas. When lung tumors were split into two groups regarding the presence or absence of metastatic lymph nodes, scores were not significantly different.
When the multi-tissue TMA BCN962 was stained for H3K79me3 and Ki67 on two serial slides, both staining did not show the same repartition in normal tissues (n = 38) compared to malignant tissues (n = 42). H3K79me3 staining was lower in tumors compared to normal tissues (Fig. 7b; p = 9.994 × 10−7), and in tumors, a negative correlation was found between Ki67 and H3K79me3 (n = 42, r = −0.52984, p = 0.0003). In support of these observations, examples of staining are given for different tumors and normal tissues (Fig. 7c). Interestingly, in the normal skin, strong Ki67-positive cells are localized in the basal cell layer, whereas H3K79me3 staining is faint and extends beyond that zone. These results suggest that H3K79me3 is decreased in tumors compared to that in normal tissues and not associated with cell proliferation.
Epigenetics is involved in EMT to repress epithelial gene expression and to stimulate mesenchymal marker expression and, as a consequence, contributes to cancer invasion, metastasis, immune surveillance escape, and resistance to treatments [6, 7]. In this study using a model of TGF-β1-induced EMT in two NSCLC cancer cell lines, we found a decrease in global H3K79me3 staining with a reduced expression of the corresponding histone methyltransferase DOT1L upon TGF-β1 treatment. Expression of several genes of interest was increased such as PD-L1 involved in the inhibition of immune checkpoints and guidance molecules like SEMA3C and its receptor NRP2. We also found that a co-treatment with DOT1L inhibitors associated with HDAC or a bromodomain inhibitor was efficient to partially reverse modifications in gene expression induced by TGF-β1. In addition, in different human tissues and corresponding tumors, H3K79me3 was generally less in tumors compared to normal tissues, and an inverse correlation was found between H3K79me3 and Ki67, a marker of cell proliferation, suggesting that H3K79me3 is not requested for cell proliferation.
Among genes that showed modified expression upon TGF-β1 treatment, guidance molecules are of particular interest and represent potential new targets in cancers. We showed that both the cell surface receptor NRP2 and SEMA3C (one of NRP2 ligands) were overexpressed during TGF-β1 treatment. NRP2, one of the two members of the neuropilin family, plays an essential role in EMT induction through non-canonical TGF-β1 signaling involving ERK . Neuropilins are pleiotropic cell surface co-receptors for some secreted members of class 3 semaphorins (SEMA3) including SEMA3F and SEMA3C, integrins, and other ligands like VEGF and growth factors such as TGF-β1 . Neuropilins (NRP1 and NRP2) are overexpressed in several cancers, and their expressions correlate with increased invasion and poor prognostic in lung cancer [28, 39]. Recently, it was demonstrated that NRP2 regulates mTOR signaling , β-catenin signaling , endosome maturation, and EGFR trafficking sustaining cancer development . Therefore, increased NRP2 expression upon TGF-β1 treatment strongly supports a function of NRP2 for TGF-β1 response in EMT induction. Semaphorins were initially described as guidance molecules involved in growth cone migration but were further involved in developmental and pathologic processes including cancer [43,44,45,46], and SEMA3C is generally described as a tumor promoting semaphorin  Thus, both NRP2 and SEMA3C overexpression might facilitate TGF-β1 signaling and tumorigenesis, and targeting NRP2 and SEMA3C would be useful to reduce EMT.
We also found a significant overexpression of PD-L1 at the membrane of H358 cells during TGF-β1-induced EMT. PD-L1 (also known as B7-H1 or CD274) is one of the ligands of PD-1, an immune checkpoint which prevents T cell activation and limit autoimmunity leading to self-tolerance . We verified that PD-1 was not expressed on H358 cells (data not shown), suggesting absence of autocrine function of PD-L1/PD1 and cell-intrinsic PD-1 pathway as shown in melanoma cells . PD-L1 is overexpressed in cancer and is involved in tumor immune escape, leading to cancer development and metastasis [30, 49, 50]. It was associated with poor clinical outcomes in several types of cancer, including NSCLC [29, 51,52,53]. By an integrated analysis of three independent large datasets, PD-L1 was associated with lung adenocarcinomas displaying a “mesenchymal” phenotype . Indeed, a molecular link was found between PD-L1 and EMT where PD-L1 is negatively regulated by miR-200, and this inhibition is relieved by ZEB1, an EMT activator induced by TGF-β1 . PD-L1 upregulation would be a consequence of EMT induction. Therefore, compounds that partially reverse EMT would target PD-L1 as well and reduce tumor progression and metastasis by restoring the immune response. In fact, pembrolizumab and nivolumab targeting PD-1 were recently approved for NSCLC to prevent PD-L1 binding .
Histone methylation is a well-balanced mechanism guided by histone methyltransferases and demethylases. Disruptor of telomeric silencing 1-like (DOT1L) is a histone methyltransferase for H3K79, and its recruitment is ubiquitously coupled with transcription. Strong similarities were noticed between H3K79 and H3K4 methylation patterns, and H3K79me3, although present at very low levels, peaks just behind the transcriptional start site and gradually decreases on the gene body [56, 57]. However, the role of DOT1L in gene expression is still controversial, but reports of H3K79me3 association with gene repression might be linked to technical issues raised by the antibody quality and absence of SDS during chromatin preparation . DOT1L is also involved in chromosome integrity and heterochromatin formation. DOT1L is well-known to be associated with pathological functions where interactors recruit DOT1L to specific gene regions, like in MLL-rearranged leukemias which present abnormal H3K79 methylation pattern at HOX loci [58,59,60,61]. In breast cancers, DOT1L is associated with poorer survival and aggressiveness, and DOT1L can cooperate with c-MYC and histone acetyltransferases to activate EMT and enhance cancer stem cell-like properties . In lung adenocarcinomas, DOT1L has been recently described to be mutated in 3% of tumors, suggesting abnormal H3K79me in a subset of samples . In contrast to breast cancers, we found a decrease of global H3K79me3 in TFG-β1-induced EMT in NSCLC cells. We suggest that H3K79me3 decrease is a mark of TGF-β1 treatment and not a general mark of EMT since additional NSCLC cell lines with a mesenchymal phenotype like H661 and H460 show higher H3K79me3 than the more epithelial H358 cells (Additional file 4). In addition, H3K79me3 level does not correlate with DOT1L expression (Additional file 4). Therefore, DOT1L activity and association with specific partners need to be considered. Also, RE-IIBP expression and activity should be studied since it was recently identified as a histone methyltransferase for H3K79 . In addition, H3K79me3 demethylation could be involved in H3K79me3 decrease, but no specific demethylase has been described to our knowledge yet. Cross-talks with other histone marks should be considered as well. Indeed, H2B ubiquitination is a prerequisite for DOT1L enzymatic activity . Although we did not find gross difference of H2BK120ub level in H358 and A549 cells after TGF-β1 treatment, H2BK120me1 and HBK120ac were decreased suggesting more complex cross-talks between histone marks.
In lung tumors, H3K79me3 intensity was statistically higher in large cell carcinomas than in adenocarcinomas. Would this difference be related to poor differentiation? The answer cannot be given because of lack of information with commercial TMAs about the status of BRG1, MYC, and EGFR. Of note, loss of function of BRG1 in about 12% of lung adenocarcinomas was shown to be associated with poor differentiation with EMT features and poor survival [64, 65]. Considering that up to one third or more of large cell carcinomas have alteration of BRG1 or other SWI/SNF subunits, it seems reasonable to consider that large cell carcinomas show poor differentiation. However, to associate DOT1L with poor differentiation would be in contradiction with the fact that DOT1L reduced cell reprogramming from a somatic differentiated state to an undifferentiated pluripotent state in several models . In addition, we found H3K79me3 decrease in EMT known to generate gain of stem cell properties in models including immortalized mammary epithelial cells . Therefore, better characterization of these large cell carcinomas is requested since, before the recent 2015 WHO classification , they were a heterogeneous group of tumors that makes interpretation difficult with the samples we analyzed.
The presence or absence of metastatic lymph nodes did not affect H3K79me3 scores in primary tumors. Interestingly, in a panel of human normal tissues and cancers, H3K79me3 level seems to be less present in tumors in comparison to normal tissues and not associated with cell proliferation. However, in both H358 and A549 cell lines, cell proliferation is not the cause of H3K79me3 decrease as it was not affected by TGF-β1 for 48 h.
We tried epigenetic treatments to modulate EMT. DOT1L inhibitors did not induce EMT as might be expected if H3K79me3 decrease was one of the causes of EMT. They were rather not effective on EMT target genes. Of interest, DOT1L was recently shown to methylate non-histone proteins like the androgen receptor through recruitment of the PRNCR1 long non-coding RNA for subsequent activation of this receptor in prostate cancer after looping between enhancer and promoter sequences . Therefore, complex responses are expected with DOT1L inhibitors. Targeting selective histone methyltransferase is a promising approach, and DOT1L inhibitors are in clinical trials for leukemias with MLL rearrangements .
The best combination to partially reverse TGF-β1-induced EMT to a more epithelial phenotype was SAHA (a HDAC inhibitor) associated with a bromodomain inhibitor. With this treatment, NRP2, SEMA3C, and PD-L1 expression were reduced. For other classical EMT genes, the response was gene-specific. More experiments would be necessary with dose response for each compound alone or in combination, associated with functional tests for cell migration and invasion. In addition, this first screening would be improved by use of a 3D cell culture model. Several studies highlight the benefits of epigenetic-based therapeutic strategy in mouse models. Combining the bromodomain inhibitor JQ1 with the histone deacetylase inhibitor SAHA in pancreatic cancer inhibits both MYC activity and inflammatory signals as well as in an established adenocarcinoma lung cancer model with KRasG12D mutation and p53 loss . In a neuroblastoma mouse model, JQ1 and panabinostat (another HDAC inhibitor) in combination, synergistically reduced N-MYC protein and tumor progression . As MYC is often overexpressed in tumors including lung cancers and directly binds PD-L1 promoter , these co-treatments are expected to reduce PD-L1 expression. From a mechanistic point, bromodomains might be efficient by targeting BRD4 binding to acetylated H4K5 at proximity to superenhancers at highly transcribed genes. Di-methylation of H3K79 mediated by DOT1L is involved in this process, as described in MLL leukemia, suggesting an unrecognized functional interplay between BRD4 and DOT1L .
In summary, we described epigenetic changes that decrease global methylation of H3K79 during TGF-β1-induced EMT, and we suggest a potential therapy with HDAC inhibitors associated to a bromodomain inhibitor to partially reverse EMT. However, clinical applications need careful characterization of tumors, and caution must be exerted with these treatments not to induce MET (mesenchymal-to-epithelial transition), particularly in cases in which tumor cells have already disseminated [71, 72].
Non-small lung cancer cell
Programmed death-ligand 1
Craene BD, Berx G. Regulatory networks defining EMT during cancer initiation and progression. Nat Rev Cancer. 2013;13:97–110.
Nieto MA, Huang RY-J, Jackson RA, Thiery JP. EMT: 2016. Cell. 2016;166:21–45.
Thiery JP, Acloque H, Huang RYJ, Nieto MA. Epithelial-mesenchymal transitions in development and disease. Cell. 2009;139:871–90.
Fischer KR, Durrans A, Lee S, Sheng J, Li F, Wong STC, et al. Epithelial-to-mesenchymal transition is not required for lung metastasis but contributes to chemoresistance. Nature. 2015;527:472–6.
Zheng X, Carstens JL, Kim J, Scheible M, Kaye J, Sugimoto H, et al. Epithelial-to-mesenchymal transition is dispensable for metastasis but induces chemoresistance in pancreatic cancer. Nature. 2015;527:525–30.
Cieslik M, Hoang SA, Baranova N, Chodaparambil S, Kumar M, Allison DF, et al. Epigenetic coordination of signaling pathways during the epithelial-mesenchymal transition. Epigenetics Chromatin. 2013;6:28.
Tam WL, Weinberg RA. The epigenetics of epithelial-mesenchymal plasticity in cancer. Nat Med. 2013;19:1438–49.
Thoma F, Koller T, Klug A. Involvement of histone H1 in the organization of the nucleosome and of the salt-dependent superstructures of chromatin. J Cell Biol. 1979;83:403–27.
Jenuwein T, Allis CD. Translating the histone code. Science. 2001;293:1074–80.
Arrowsmith CH, Bountra C, Fish PV, Lee K, Schapira M. Epigenetic protein families: a new frontier for drug discovery. Nat Rev Drug Discov. 2012;11:384–400.
Torre LA, Siegel RL, Ward EM, Jemal A. Global Cancer Incidence and Mortality Rates and Trends—an update. Cancer Epidemiol Biomarkers Prev [Internet]. 2015. [cited 2016 Aug 23]; Available from: http://cebp.aacrjournals.org/content/early/2015/12/10/1055-9965.EPI-15-0578
Balgkouranidou I, Liloglou T, Lianidou ES. Lung cancer epigenetics: emerging biomarkers. Biomarkers Med. 2013;7:49–58.
Liloglou T, Bediaga NG, Brown BRB, Field JK, Davies MPA. Epigenetic biomarkers in lung cancer. Cancer Lett. 2014;342:200–12.
Van Den Broeck A, Ozenne P, Eymin B, Gazzeri S. Lung cancer. Cell Adhes Migr. 2010;4:107–13.
Cancer Genome Atlas Research Network. Comprehensive molecular profiling of lung adenocarcinoma. Nature. 2014;511:543–50.
Collisson EA, Campbell JD, Brooks AN, Berger AH, Lee W, Chmielecki J, et al. Comprehensive molecular profiling of lung adenocarcinoma. Nature. 2014;511:543–50.
Campbell JD, Alexandrov A, Kim J, Wala J, Berger AH, Pedamallu CS, et al. Distinct patterns of somatic genome alterations in lung adenocarcinomas and squamous cell carcinomas. Nat Genet. 2016;48:607–16.
Cancer Genome Atlas Research Network. Comprehensive genomic characterization of squamous cell lung cancers. Nature. 2012;489:519–25.
Behrens C, Soto LMS, Lin HY, Yuan P, Tang X, Kadara H, et al. EZH2 Protein Expression Associates With the Early Pathogenesis, Tumor Progression and Prognosis of Non-small Cell Lung Carcinoma. Clin Cancer Res. 2013;19(23):6556–65.
Huqun, Ishikawa R, Zhang J, Miyazawa H, Goto Y, Shimizu Y, et al. Enhancer of zeste homolog 2 is a novel prognostic biomarker in nonsmall cell lung cancer. Cancer. 2012;118:1599–606.
Kikuchi J, Kinoshita I, Shimizu Y, Kikuchi E, Konishi J, Oizumi S, et al. Distinctive expression of the polycomb group proteins Bmi1 polycomb ring finger oncogene and enhancer of zeste homolog 2 in nonsmall cell lung cancers and their clinical and clinicopathologic significance. Cancer. 2010;116:3015–24.
Takawa M, Masuda K, Kunizaki M, Daigo Y, Takagi K, Iwai Y, et al. Validation of the histone methyltransferase EZH2 as a therapeutic target for various types of human cancer and as a prognostic marker. Cancer Sci. 2011;102:1298–305.
Diaz-Lagares A, Mendez-Gonzalez J, Hervas D, Saigi M, Pajares MJ, Garcia D, et al. A novel epigenetic signature for early diagnosis in lung cancer. Clin Cancer Res. 2016;22:3361–71.
Roche J, Nasarre P, Gemmill R, Baldys A, Pontis J, Korch C, et al. Global decrease of histone H3K27 acetylation in ZEB1-induced epithelial to mesenchymal transition in lung cancer cells. Cancers (Basel). 2013;5:334–56.
Gemmill RM, Roche J, Potiron VA, Nasarre P, Mitas M, Coldren CD, et al. ZEB1-responsive genes in non-small cell lung cancer. Cancer Lett. 2011;300:66–78.
Schliekelman MJ, Taguchi A, Zhu J, Dai X, Rodriguez J, Celiktas M, et al. Molecular portraits of epithelial, mesenchymal, and hybrid states in lung adenocarcinoma and their relevance to survival. Cancer Res. 2015;75:1789–800.
Nasarre P, Potiron V, Drabkin H, Roche J. Guidance molecules in lung cancer. Cell Adhes Migr. 2010;4:130–45.
Nasarre P, Gemmill RM, Potiron VA, Roche J, Lu X, Barón AE, et al. Neuropilin-2 is upregulated in lung cancer cells during TGF-β1-induced epithelial-mesenchymal transition. Cancer Res. 2013;73:7111–21.
Sun J-M, Zhou W, Choi Y-L, Choi S-J, Kim SE, Wang Z, et al. Prognostic significance of PD-L1 in patients with non-small cell lung cancer: a large cohort study of surgically resected cases. J Thorac Oncol. 2016;11:1003–11.
Lou Y, Diao L, Cuentas ERP, Denning WL, Chen L, Fan YH, et al. Epithelial-mesenchymal transition is associated with a distinct tumor microenvironment including elevation of inflammatory signals and multiple immune checkpoints in lung adenocarcinoma. Clin Cancer Res. 2016;22(14):3630–42.
Fang W, Zhang J, Hong S, Zhan J, Chen N, Qin T, et al. EBV-driven LMP1 and IFN-γ up-regulate PD-L1 in nasopharyngeal carcinoma: Implications for oncotargeted therapy. Oncotarget. 2014;5:12189–202.
Frederiks F, Tzouros M, Oudgenoeg G, van Welsem T, Fornerod M, Krijgsveld J, et al. Nonprocessive methylation by Dot1 leads to functional redundancy of histone H3K79 methylation states. Nat Struct Mol Biol. 2008;15:550–7.
Zhu B, Zheng Y, Pham A-D, Mandal SS, Erdjument-Bromage H, Tempst P, et al. Monoubiquitination of human histone H2B: the factors involved and their roles in HOX gene regulation. Mol Cell. 2005;20:601–11.
Shahbazian MD, Zhang K, Grunstein M. Histone H2B ubiquitylation controls processive methylation but not monomethylation by Dot1 and Set1. Mol Cell. 2005;19:271–7.
Johnsen SA. The enigmatic role of H2Bub1 in cancer. FEBS Lett. 2012;586:1592–601.
Mazur PK, Herner A, Mello SS, Wirth M, Hausmann S, Sánchez-Rivera FJ, et al. Combined inhibition of BET family proteins and histone deacetylases as a potential epigenetics-based therapy for pancreatic ductal adenocarcinoma. Nat Med. 2015;21:1163–71.
Shahbazi J, Liu PY, Atmadibrata B, Bradner JE, Marshall GM, Lock RB, et al. The bromodomain Inhibitor JQ1 and the histone deacetylase inhibitor panobinostat synergistically reduce N-Myc expression and induce anticancer effects. Clin Cancer Res. 2016;22:2534–44.
Prud’homme GJ, Glinka Y. Neuropilins are multifunctional coreceptors involved in tumor initiation, growth, metastasis and immunity. Oncotarget. 2012;3:921–39.
Lantuéjoul S, Constantin B, Drabkin H, Brambilla C, Roche J, Brambilla E. Expression of VEGF, semaphorin SEMA3F, and their common receptors neuropilins NP1 and NP2 in preinvasive bronchial lesions, lung tumours, and cell lines. J Pathol. 2003;200:336–47.
Nakayama H, Bruneau S, Kochupurakkal N, Coma S, Briscoe DM, Klagsbrun M. Regulation of mTOR signaling by semaphorin 3F-neuropilin 2 interactions in vitro and in vivo. Sci Rep. 2015;5:11789.
Samuel S, Gaur P, Fan F, Xia L, Gray MJ, Dallas NA, et al. Neuropilin-2 mediated β-catenin signaling and survival in human gastro-intestinal cancer cell lines. PLoS One. 2011;6:e23208.
Dutta S, Roy S, Polavaram NS, Stanton MJ, Zhang H, Bhola T, et al. Neuropilin-2 regulates endosome maturation and EGFR trafficking to support cancer cell pathobiology. Cancer Res. 2016;76:418–28.
Gu C, Giraudo E. The role of semaphorins and their receptors in vascular development and cancer. Exp Cell Res. 2013;319:1306–16.
Nasarre P, Gemmill RM, Drabkin HA. The emerging role of class-3 semaphorins and their neuropilin receptors in oncology. Onco Targets Ther. 2014;7:1663–87.
Neufeld G, Mumblat Y, Smolkin T, Toledano S, Nir-Zvi I, Ziv K, et al. The role of the semaphorins in cancer. Cell Adhes Migr. 2016:10(6):652–674.
Rehman M, Tamagnone L. Semaphorins in cancer: biological mechanisms and therapeutic approaches. Semin Cell Dev Biol. 2013;24:179–89.
Francisco LM, Sage PT, Sharpe AH. The PD-1 pathway in tolerance and autoimmunity. Immunol Rev. 2010;236:219–42.
Kleffel S, Posch C, Barthel SR, Mueller H, Schlapbach C, Guenova E, et al. Melanoma cell-intrinsic PD-1 receptor functions promote tumor growth. Cell. 2015;162:1242–56.
Chen J, Jiang CC, Jin L, Zhang XD. Regulation of PD-L1: a novel role of pro-survival signalling in cancer. Ann Oncol. 2016;27:409–16.
He J, Hu Y, Hu M, Li B. Development of PD-1/PD-L1 pathway in tumor immune microenvironment and treatment for non-small cell lung cancer. Sci Rep. 2015;5:13110.
Boland JM, Kwon ED, Harrington SM, Wampfler JA, Tang H, Yang P, et al. Tumor B7-H1 and B7-H3 expression in squamous cell carcinoma of the lung. Clinical Lung Cancer. 2013;14:157–63.
Velcheti V, Schalper KA, Carvajal DE, Anagnostou VK, Syrigos KN, Sznol M, et al. Programmed death ligand-1 expression in non-small cell lung cancer. Lab Investig. 2014;94:107–16.
Zou W, Chen L. Inhibitory B7-family molecules in the tumour microenvironment. Nat Rev Immunol. 2008;8:467–77.
Chen L, Gibbons DL, Goswami S, Cortez MA, Ahn Y-H, Byers LA, et al. Metastasis is regulated via microRNA-200/ZEB1 axis control of tumour cell PD-L1 expression and intratumoral immunosuppression. Nat Commun. 2014;5:5241.
Viteri S, González-Cao M, Barrón F, Riso A, Rosell R. Results of clinical trials with anti-programmed death 1/programmed death ligand 1 inhibitors in lung cancer. Transl Lung Cancer Res. 2015;4:756–62.
Steger DJ, Lefterova MI, Ying L, Stonestrom AJ, Schupp M, Zhuo D, et al. DOT1L/KMT4 recruitment and H3K79 methylation are ubiquitously coupled with gene transcription in mammalian cells. Mol Cell Biol. 2008;28:2825–39.
Vlaming H, van Leeuwen F. The upstreams and downstreams of H3K79 methylation by DOT1L. Chromosoma. 2016;125:593–605.
Daigle SR, Olhava EJ, Therkelsen CA, Basavapathruni A, Jin L, Boriack-Sjodin PA, et al. Potent inhibition of DOT1L as treatment of MLL-fusion leukemia. Blood. 2013;122:1017–25.
Deshpande AJ, Bradner J, Armstrong SA. Chromatin modifications as therapeutic targets in MLL-rearranged leukemia. Trends Immunol. 2012;33:563–70.
Gilan O, Lam EYN, Becher I, Lugo D, Cannizzaro E, Joberty G, et al. Functional interdependence of BRD4 and DOT1L in MLL leukemia. Nat Struct Mol Biol. 2016;23:673–81.
McLean CM, Karemaker ID, van Leeuwen F. The emerging roles of DOT1L in leukemia and normal development. Leukemia. 2014;28:2131–8.
Cho M-H, Park J-H, Choi H-J, Park M-K, Won H-Y, Park Y-J, et al. DOT1L cooperates with the c-Myc-p300 complex to epigenetically derepress CDH1 transcription factors in breast cancer progression. Nat Commun. 2015;6:7821.
Woo Park J, Kim K-B, Kim J-Y, Chae Y-C, Jeong O-S, Seo S-B. RE-IIBP methylates H3K79 and induces MEIS1-mediated apoptosis via H2BK120 ubiquitination by RNF20. Sci Rep. 2015;5:12485.
Matsubara D, Kishaba Y, Ishikawa S, Sakatani T, Oguni S, Tamura T, et al. Lung cancer with loss of BRG1/BRM, shows epithelial mesenchymal transition phenotype and distinct histologic and genetic features. Cancer Sci. 2013;104:266–73.
Oike T, Ogiwara H, Tominaga Y, Ito K, Ando O, Tsuta K, et al. A synthetic lethality-based strategy to treat cancers harboring a genetic deficiency in the chromatin remodeling factor BRG1. Cancer Res. 2013;73:5508–18.
Mani SA, Guo W, Liao M-J, Eaton EN, Ayyanan A, Zhou AY, et al. The epithelial-mesenchymal transition generates cells with properties of stem cells. Cell. 2008;133:704–15.
Travis WD, Brambilla E, Nicholson AG, Yatabe Y, Austin JHM, Beasley MB, et al. The 2015 World Health Organization classification of lung tumors. J Thorac Oncol. 2015;10:1243–60.
Yang L, Lin C, Jin C, Yang JC, Tanasa B, Li W, et al. lncRNA-dependent mechanisms of androgen-receptor-regulated gene activation programs. Nature. 2013;500:598–602.
Morera L, Lübbert M, Jung M. Targeting histone methyltransferases and demethylases in clinical trials for cancer therapy. Clin Epigenetics. 2016;8:57.
Casey SC, Tong L, Li Y, Do R, Walz S, Fitzgerald KN, et al. MYC regulates the antitumor immune response through CD47 and PD-L1. Science. 2016;352:227–31.
Marcucci F, Stassi G, De Maria R. Epithelial-mesenchymal transition: a new target in anticancer drug discovery. Nat Rev Drug Discov. 2016;15:311–25.
Mishra VK, Johnsen SA. Targeted therapy of epigenomic regulatory mechanisms controlling the epithelial to mesenchymal transition during tumor progression. Cell Tissue Res. 2014;356:617–30.
This work benefited from the core facilities of the ImageUP platform (STIM-ERL7368, University of Poitiers, France). The authors thank Dr. A. Cantereau and Dr. A. Delwail for their expertise for the confocal microscopy and flow cytometry, respectively. They acknowledge Dr. J. Imbert (Luminy University, Marseille, France) and Dr. M. Ramnath (Eurofins Cerep, France) for the fruitful discussions and S-M. Huruguen (Eurofins Cerep, France) for the technical assistance. They are grateful to Dr. R. Hamamoto (University of Chicago, IL) for providing the H2BK120me1 antibody.
This work was supported by FEDER 38345-2013 and Région Poitou-Charentes (France).
Availability of data and materials
All data generated and analyzed during this study are included in this published article and in supplementary information files.
Ethics approval and consent to participate
Not applicable: Human tissue arrays were commercially available and purchased from USBiomax (Rockville, MD, USA).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Evanno, E., Godet, J., Piccirilli, N. et al. Tri-methylation of H3K79 is decreased in TGF-β1-induced epithelial-to-mesenchymal transition in lung cancer. Clin Epigenet 9, 80 (2017). https://doi.org/10.1186/s13148-017-0380-0
- Lung cancer
- Epigenetic treatment