NR2E1 and LSD1 are essential for the proliferation of BTICs
To study the role of NR2E1 and LSD1 in BTICs, two BTIC lines, BTIC-1 and BTIC-2 derived from Nestin-TV-a mice were employed for the following experiments. These BTICs can efficiently grow mouse brain tumors after transplantation in recipient mice(Zhu et al., 2014). Like neural stem cells (NSCs), BTICs could be maintained in monolayer or non-adherent suspension culture (Figure 1a) (Zhu et al., 2014). The capacity of BTICs to form tumor spheres suggests their stem cell properties. NR2E1 and LSD1 are highly expressed in NSCs and promote the self-renewal of NSCs (Sun et al., 2010). Real-time PCR and western blot assays showed that both NR2E1 and LSD1 were more highly expressed in BTICs than NSCs (Figure 1b and 1c). We therefore set to examine whether NR2E1 and LSD1 coordinate with each other to regulate the proliferation of BTICs. To knock down Nr2e1 and Lsd1, we designed two different shRNAs for each gene and cloned them into pSuper-puro vector. We transfected these shRNAs in to BTICs respectively by electroporation and selected the transfected cells with puromycin(Sun et al., 2010). Three days after puromycin selection, Nr2e1 and Lsd1 shRNA knockdown led to a significantly lower amount of BTICs in culture (Figure 1d and Supplementary Figure 1a). Real-time PCR revealed that the expression of Nr2e1 and Lsd1 was downregulated to about 20% to 40% of their original levels by their respective shRNAs (Figure 1e and Supplementary Figure 1b). We then examined cell viability by the MTT cell proliferation assay, which revealed that the knockdown of Nr2e1 and Lsd1 resulted in reduced cell viability compared to the control knockdown BTICs (Figure 1f and Supplementary Figure 1c). These results suggest that Nr2e1 and Lsd1 are required for the proliferation of BTICs.
NR2E1 and LSD1 synergistically repress PTEN to promote BTIC proliferation
To examine whether NR2E1 and LSD1 play the same role in BTICs as in NSCs, we performed co-immunoprecipitation assay using whole cell lysate. An anti-NR2E1 antibody could pull down endogenous LSD1, but control IgG could not, suggesting that NR2E1 and LSD1 form a complex in BTICs (Figure 2a). To investigate whether NR2E1 and LSD1 also regulate Pten in BTICs, we knocked down the expression of Nr2e1 and Lsd1 in BTICs by shRNAs. The downregulation of NR2E1 and LSD1 led to the upregulation of PTEN at both the mRNA and protein levels (Figure 2b and 2c, Supplementary Figure 2a). To examine whether PTEN is one of the major effectors of NR2E1 and LSD1, we generated shRNA that could efficiently downregulate Pten expression (Figure 2d). Downregulation of Pten expression in Nr2e1 knockdown BTICs or Lsd1 knockdown BTICs could rescue the cell death phenotype caused by only knockdown of Nr2e1 or Lsd1 (Figure 2e- 2f, Supplementary Figure 2b-2c). To further examine whether BTICs employed the same NR2E1-LSD1 regulatory mechanism as NSCs, we performed chromatin immunoprecipitation (ChIP) assay to investigate NR2E1 and LSD1 binding profile. The ChIP assay revealed that NR2E1 and LSD1 both bind to the promoter of Pten (Figure 2g). To examine whether NR2E1 and LSD1 are functional at the promoter of Pten, we performed ChIP assay with antibodies against H3K4me1 and H3K4me2 using chromatin extracted from Nr2e1 or Lsd1 knockdown BTICs. It turned out that both downregulation of NR2E1 and LSD1 led to upregulated enrichment of H3K4me1 and H3K4me2 at the Pten promoter, suggesting that NR2E1 and LSD1 indeed directly repress Pten in BTICs by demethylating H3K4me and H3K4me2 at its promoter (Figure 2h).
Prediction of Lsd1 peptides involved in the NR2E1-LSD1 interaction
To understand how NR2E1 and LSD1 synergistically function, we employed Amide Hydrogen/Deuterium Exchange and Mass Spectrometry (HDX-MS) to investigate the interaction between NR2E1 and LSD1. A total of 56 pepsin digested fragments covering about 80% of the LSD1 primary sequence were identified and analyzed. The difference in deuterium uptake for all the fragments between LSD1 alone and LSD1:NR2E1 complex was measured at 30 seconds, 1, 2, 5 and 10 minutes. As the difference in deuterium exchange was maximum at 1 minute, the deuterium uptake for each peptide in the 1 minute samples was used to monitor the effects of NR2E1 binding with LSD1. A number of regions in LSD1 showed decreased exchange upon interactions with NR2E1, but the maximum difference occurred at the regions within the AO domain. Peptides 333-350, 333-353 and 354-377 from the AO domain showed a significant decrease of 2.2, 2.5 and 4.0 deuterons, respectively (Figure 3a). Also, in the LSD1:NR2E1 complex, peptides 196-211, 197-211, 320-332, 333-344, 378-385, 419-441, 481-501, 498-511, 500-510, 537-546, 601-614 and 623-650 from the SWIRM, Tower and AO domains showed a decrease of about 1.6 deuterons. Mass spectral isotope envelopes for four peptides, 197-211, 354-377, 481-501 and 537-546, from 1 minute HDX samples showed the most significant difference between LSD1 alone and NR2E1-LSD1 complex after deuterium uptake, suggesting these LSD1 peptides may be involved in forming a complex with NR2E1 (Figure 3b). The difference in deuterium uptake for each peptide was calculated and the results from 1 minute samples were mapped onto the crystal structure of LSD1 (PDB ID: 2Z3Y) (Figure 3c).
The crystal structure of the NR2E1 ligand binding domain (LBD) has been solved (Zhi et al., 2015). NR2E1 LBD interacts with both the SWIRM and AO domains of LSD1 (Yokoyama et al., 2008). We have predicted the interaction between NR2E1 LBD (PDB code: 4XAI) and LSD1 (PDB code: 3ZMU) using the ZDOCK program (Chen et al., 2003). LSD1 residues 256-333, 410-435, 575-626, 712-744 and 792-819 that are not supposed to interact with NR2E1 were set as block residues. As a result, 3600 docking complexes were generated and clustered into 650 groups using the MMSTB clustering method (Feig et al., 2004) with a root-mean-square deviation (RMSD) cutoff at 8 Å. The buried solvent-accessible area, which is deemed as the possible binding surface of NR2E1 LBD by LSD1 peptides 197-211, 354-377, 481-501 and 537-546, was calculated by Naccess (Hubbard and Thornton, 1993) with a probe radius of 1.4 Å. LSD1 peptides 197-211, 354-377 and 481-501 individually showed varying predicted binding surface on NR2E1. However, LSD1-537-546 showed almost no predicted binding surface, suggesting that LSD1-537-546 is not likely involved in the interaction between LSD1 and NR2E1 LBD (Figure 3d).
Role of LSD1 peptides in NR2E1-LSD1 interaction
To further characterize the role of the LSD1 peptides 197-211, 354-377, 481-501 and 537-546 identified by HDX-MS, we generated Flag-tagged Lsd1 mutant clones by deleting the peptide encoding regions. These clone plasmids were then co-transfected with a plasmid expressing HA-NR2E1 into 293T cells for co-immunoprecipitation assay. Western blot revealed that while Flag-LSD1-∆197-211, Flag-LSD1-∆354-377 and Flag-LSD1-∆537-546 could be stably expressed, deletion of residues 481-501 of LSD1 however led to no detectable protein, suggesting this region is critical for the stable expression of LSD1 (Figure 4a). Then we performed immunoprecipitation with an anti-HA antibody using whole cell lysate and followed by immunoblotting with an anti-Flag antibody. It turned out that HA-NR2E1could pull down Flag-LSD1, Flag-LSD1-∆354-377 and Flag-LSD1-∆537-546, but not Flag-LSD1-∆197-211. This suggests that deletion of peptide LSD1-354-377 or LSD1 537-546 did not disturb the interaction between LSD1 and NR2E1, only deletion of LSD1-197-211 destroyed the NR2E-LSD1 complex (Figure 4a). Therefore, LSD1-197-211 is essential for NR2E1-LSD1 complex formation.
We next investigated the function of overexpression of these peptides in BTICs. We cloned the four peptides into the pCAG-puro plasmids and expressed them as Flag-tag peptides. Twelve hours after transfection, puromycin was added to select the transfected cells. At day 3 after puromycin selection, cells were harvested to check protein expression. Western blot with an anti-Flag antibody confirmed that all four Flag-tagged LSD1 peptides were expressed. Interestingly, the level of PTEN was increased in Flag-LSD1-197-211 overexpressed BTICs compared to GFP overexpressed BTICs (Figure 4b). Immunostaining BTICs at this stage with an anti-Ki67 antibody revealed that overexpression of LSD1-197-211 led to fewer Ki67 positive cells (Figure 4c). Further MTT assay confirmed that LSD1-197-211 transfected BTICs showed the most drastic reduction of viable cells, while LSD1-354-377 only slightly decreased BTICs and both LSD1-481-501 and LSD1-537-546 showed almost no effect (Figure 4d).
As transient express LSD1-197-211 could inhibit the self-renewal of BTICs, we next generated doxycycline inducible lentivirus to stably express LSD1-197-211 peptide-P2A-EGFP in BTICs. LSD1-197-211 peptide and EGFP are connected with P2A, a self-cleavage peptide(Supplementary Figure 3a). After translation, LSD1-197-211 and EGFP are cut apart at P2A site. Hence,GFP protein level can reflect LSD1-197-211 level in transduced cells. We first infected 293T cells with these lentiviruse. After doxycycline induction for three days, it is clearly that both GFP lentivirus and LSD1-197-211-GFP lentivirus express similar level of GFP. Therefore, LSD1-197-211 did not inhibit the growth of 293T cells (Supplementary Figure 3b). Next, we infected BTICs with the same batch of lentivirus and purified the transduced cells by puromycin selection. We then induced the peptide expression with doxycycline. 60 hours after induction, we performed immunostaining with antibody against cleaved Caspase-3, the marker of apoptosis. It turned out that BTICs that expressed LSD1-197-211 and GFP were often cleaved Caspase-3 positive, while it was not the case for the control GFP overexpressed BTICs (Supplementary Figure 4a). This result suggests that the expression of LSD1-197-211 in BTICs leads to apoptosis of the cells. Further extending the cell culture time, we observed that GFP positive LSD1-197-211 BTICs gradually lost their shape and detached from the plate. Eventually very rare GFP positive cells in LSD1-197-211 lentivirus transduced BTIC culture were observed (Supplementary Figure 4b). Analysis of GFP positive cell ratio with flow cytometry revealed that in contrast to the 78.6% GFP positive cells in the GFP lentivirus transduced BTICs, only 0.71% of the LSD1-197-211 lentivirus transduced BTICs showed weak GFP expression, suggesting that LSD1-197-211 peptide expressed BTICs were dead (Supplementary Figure. 4c). Soft agar colony formation assay also revealed that LSD1-197-211 overexpression drastically reduced the colony formation capacity of BTICs (Figure 4e). Furthermore, transwell assay revealed that overexpression of LSD1-197-211 drastically reduced the migration capacity of BTICs (Figure 4f). As expected, the size of tumor spheres formed by LSD1-197-211 overexpressed BTICs was smaller than GFP overexpressed BTICs (Figure 4g). Besides, the overall sphere number was also much less in LSD1-197-211 overexpressed BTICs than GFP overexpressed BTICs (Figure 4h), suggesting an inhibitory role of LSD1-197-211 on BTIC in vitro.
To examine whether the inhibitory effect of LSD1-197-211 on BTIC proliferation is due to the interference of NR2E1 and LSD1 synergistical function, we performed ChIP assay with antibodies against H3K4me1 and H3K4me2 with chromatin extracted from GFP, LSD1-197-211 and LSD1-354-377 overexpressed BTICs respectively. H3K4me1 and H3K4me2 modification at Pten promoter was significantly increased in LSD1-197-211 overexpressed BTICs, compared to the GFP overexpressed BTICs and LSD1-354-377 overexpressed BTICs (Figure 4i). These data suggest that LSD1-197-211 can inhibit BTIC proliferation by interfering the demethylation function of NR2E1-LSD1 complex on H3K4 methylation at Pten promoter.
Specificity of peptide Lsd1-197-211
Although LSD1-197-211 could interfere with the synergistic function of NR2E1 and LSD1 and inhibit the proliferation of BTICs, the specificity of this peptide is unclear. Both NR2E1 and LSD1 are highly expressed in 293T cells (Figure 5a). Knockdown of Nr2e1 by shRNA did not, however, lead to the upregulation of Pten at the mRNA level, suggesting that the NR2E1-LSD1 mechanism is not involved in the proliferation of 293T cells (Figure 5b). To test whether LSD1-197-211 had any effect on the cells that do not rely on the NR2E1-LSD1 based cell proliferation, we overexpressed GFP and LSD1 peptides in 293T cells separately with puromycin selection for three days. Western blot showed that the PTEN protein level was similar in the GFP and LSD1 peptide overexpressed BTICs (Figure 5c). MTT assay was further performed to check the cell viability. Apart from LSD1-354-377, which exhibited slight inhibition of 293T cell growth, other peptides showed no obvious inhibitory effect (Figure 5d). This result suggests that the LSD1-197-211 peptide shows relatively specific inhibition on BTICs.
To further characterize the specificity of LSD1-197-211, we determined the crystal structure of the human LSD1 SWIRM domain, residues 183-267 (Supp. Table 1). The SWIRM structure mainly contains a long central helix separating two smaller helix-loop-helix motifs at both sides (Figure 5e). The SWIRM domain highly resembles the SWIRM domain of previously determined human LSD1 crystal structure (PDB ID: 2Z3Y) (Supplementary Figure 5). Out of 687 aligned atoms, 528 atoms of the two structures can be well aligned with a root mean square of 0.418. Unlike the SWIRM domain of SWI3 and ADA2, which binds to DNA, the SWIRM domain of LSD1 can neither bind to DNA nor does it contain any typical DNA-binding patch (Da et al., 2006). One major difference between LSD1 SWIRM domain and SWI3 and ADA2 SWIRM domains is mainly confined to the N-terminus. LSD1-197-211 is located at this region of the SWIRM domain, comprising a part of helix H2, a loop and a part of the long central helix H3 (Figure 5e). This region protrudes away from the hydrophobic core formed by H6 and the AO domain, which is involved in substrate demethylation (Chen et al., 2006). It forms a stable binding pocket with the N terminal loop and provides enough space to interact with other proteins (Figure 5e). Unlike LSD1 SWIRM domain, the corresponding region of ADA2 SWIRM is blocked by the N-terminal loop, making it impossible to interact with other proteins. For the SWI3 SWIRM domain, a big loop connects the H2 and H3 helices and no binding pocket is formed in this region (Figure 5e). The structural difference among LSD1, ADA2 and SWI3 SWIRM domains substantiates their different ligand binding properties.
LSD1 197-211 inhibits the brain tumor formation of BTIC
Human and mouse NR2E1 protein sequences share more than 97% similarity,so does LSD1. In addition, human and mouse LSD1-197-211 peptides are exactly the same. Therefore, we deduced that LSD1-197-211 peptide should be able to repress human BTICs (hBTICs) as well. To test this hypothesis, we derived two hBTIC lines from glioblastoma patients (Supplementary Figure 6a) and generated doxycycline inducible LSD1-197-211-GFP and control GFP transduced hBTICs by abovementioned lentivirus system. Next, we tested the effect of LSD1-197-211 on tumor sphere formation by adding doxcycycline to the culture medium. 3 days after doxycycline induction, LSD1-197-211-GFP transduced hBTICs formed slightly smaller tumor sphere than GFP transduced hBTICs. GFP signal in LSD1-197-211-GFP transduced hBTICs was also weaker than GFP transduced hBTICs. 6 days later, the difference was more drastic (Figure 6a and Supplementary Figure 6b). Not only the number of tumor spheres formed by LSD1-197-211-GFP transduced hBTICs was less than GFP transduced hBTICs, but also the average tumor sphere size of LSD1-197-211-GFP transduced hBTICs was smaller than GFP transduced hBTICs (Figure 6b). These results confirmed that LSD1-197-211 inhibits hBTICs.
To further investigate the function of LSD1-197-211 in vivo, we transplanted the transduced hBTICs into the brains of nude mice. The next day after the transplantation, we started to feed the mice with water containing doxycycline to induce transduced gene expression. Mice transplanted with GFP transduced hBTICs showed obvious brain tumor growth as revealed by GFP positive signals, while almost no GFP positive signals were detected in mice transplanted with LSD1-197-211-GFP transduced hBTICs (Figure 6c-6d). Further Hematoxylin-eosin staining revealed obvious tumor infiltration in the brains of GFP transduced hBTIC transplanted mice, while most brains of mice transplanted with LSD1-197-211-GFP transduced hBTICs showed no or very small tumor infiltration (Figure 6e). By examining the genotype of the BTIC transplanted mouse brain via PCR, we confirmed that 11 mice were GFP hBTIC positive and 13 mice were LSD1-197-211-GFP hBTIC positive (Supplementary Figure 6c). The average brain tumor volume of LSD1-197-211-GFP hBTIC transplanted mice was obviously smaller than control BTIC transplanted mice (Supplementary Figure 6d). In addition, mice transplanted with GFP transduced hBTICs were much weaker and skinner than mice transplanted with LSD1-197-211-GFP transduced hBTICs (Supplementary Figure 6e-6f). Out of 11 GFP transduced hBTIC transplanted mice, 3 mice died before week seven after transplantation, while none of LSD1-197-211-GFP transduced hBTIC transplanted mouse died before week seven. Hence, LSD1-197-211 can suppress the tumor formed by hBTICs and improve the survival rate of hBITC transplanted mice.