2. Materials and Methods
SH-SY5Y cell culture and differentiation
Human SH-SY5Y neuroblastoma cells (ATCC® CRL-2266) were maintained in 1:1 EMEM/F12 (ATCC, Gibco), supplemented with 10% FBS and 1% penicillin/streptomycin 100X solution (Sigma Aldrich) at 37°C, 5% atmospheric CO2 in a humidified incubator. The ATCC CRL-2266 protocol was implemented for routine neuroblastoma subcultivation procedures. The cells were differentiated according to the protocol described by [
20,
21] with some modifications. Briefly, 2-3 X 10
4 cells/cm2 were seeded onto culture dishes and were allowed to be attached overnight in a standard medium with 10% FBS heat inactivated. The next day the FBS content of the culture medium was reduced to 5%, and the cells were exposed to 10µM all-trans retinoic acid (ATRA) (VWR). The cells were kept under these conditions for five days with a medium changes every second day. On the sixth day of differentiation, the cells were cultured in a standard medium without FBS, with addition of 10 µM ATRA and 50ng/ml brain-derived neurotropic factor (BDNF) (VWR) for an additional 3-4 days. Cells were ready to be used for further experiments.
Immunofluorescence microscopy
For immunofluorescence microscopy, the SH-SY5Y cultured cells were fixed with a fixing solution (4% paraformaldehyde, 4% sucrose in PBS, filtered) for 6-7 minutes at room temperature. After fixation, cells were washed three times in warmed PBS with 4% sucrose for 5 minutes at RT. Then SH-SY5Y cells were washed three times in PBS, permeabilized for 10 minutes in PBST (PBS with 0.3% Triton X-100; BioShop), and blocked for 2 h in 3% BSA (Sigma-Aldrich) in PBST at RT. The cells were then incubated with mouse monoclonal anti-MAP2 primary antibody (1:1000; M1406, Sigma-Aldrich) overnight at 4°C, washed three times in PBS, and thereafter incubated with the secondary polyclonal antibody conjugated with Alexa Fluor® 488 (goat anti-mouse, 1:500, Invitrogen) for 2h in RT. After washing, as indicated previously, coverslips with cells were mounted on glass slides with Fluoromount-G mounting medium, dried at room temperature for 5-6h and observed with ZEISS Spinning DISC confocal microscope with a 20x objective lens. Acquisition and analysis were performed with Java-based image processing program (ImageJ).
Preparation of nuclear extracts
Nuclear extracts from SH-SY5Y were prepared according to Bethyl Laboratories’ protocol with slight modifications. Briefly, cells were washed twice with cold PBS and resuspended in buffer A (1 M HEPES, pH 7.9, 1 M KCl, 0.5 M EDTA, 100% Nonindet-40 (Sigma), 1 M DTT (Sigma), and protease inhibitor cocktail (Sigma). Plate with cells was then incubated for 10 minutes at RT, scraped, and lysates were transferred to a microcentrifuge tube. Nuclei were pelleted at 4°C at 13 200 rpm for 3 minutes. Cytoplasmatic extracts were transferred to the new tube, and nuclei were resuspended in buffer B (1M HEPES, pH 7.9, 5M NaCl, O.5 M EDTA, glycerol, 1M DTT, and protease inhibitor cocktail). After pipetting, samples were vortexed every 15 minutes for 2 hours, centrifuged at 4°C at 13 200 rpm for 5 minutes, and supernatants were used as nuclear extracts. Protein concentration was estimated by The Pierce Detergent Compatible Bradford Assay Kit [Thermo Scientific #23246] using Synergy H4 (BioTek).
Transfection
Undifferentiated and differentiated SH-SY5Y cells were passaged at the same time and then cultured for five days. On the fifth day, cells were trypsinized, counted, and plated in 98-well plates at 4 x 104 cells per well so they were 70-90% confluent on the day of transfection. Undifferentiated SH-SY5Y cells were plated in a 100 µl standard medium containing serum with antibiotics. In contrast, differentiated SH-SY5Y cells were grown in a medium with 5% FBS, 10 µM ATRA and antibiotics (5th day of differentiation protocol –
Figure 1a). The 6th-day medium was changed to antibiotics-free medium, and 10 µl of the mixture containing 0,3 µl Lipofectamine TM 3000 Reagent (Thermo Fischer Scientific), 200 ng appropriate construct (Addgene), and 2µl/µg DNA P3000TM Reagent diluted in 1:1 DMEM/F12 was added to each well. After 6 hours after transfection medium was changed to standard medium in case of undifferentiated SH-SY5Y cells, and to serum-free medium containing 10 µM ATRA, and 50ng/ml BDNF in case of differentiated SH-SY5Y cells.
Figure 1.
SH-SY5Y differentiation into neurons. a. Timetable of SH-SY5Y human neuroblastoma cell line differentiation into neurons. b. The photographs of undifferentiated and differentiated SH-SY5Y neuroblastoma cells immunostained with the neuronal marker MAP-2 antibody. Alexa-488 was used to visualize MAP-2. Fluorescence images were obtained using a 20x objective. Scale bar: 100 µm. c. Expression of synaptophysin in SH-SY5Y cells and human neurons derived from SH-SY5Y cells analyzed by RT-qPCR. Data is presented as fold change in mRNA expression. The unpaired t-test test was followed by the D’Agostino and Pearson normality test. Values are given as mean ± SEM (****, p < 0.001; n=12). d. The overall expression level of the MMP-9 gene in SH-SY5Y cells and in human neurons. Data is presented as fold change in mRNA expression. The unpaired t-test test was followed by the D’Agostino and Pearson normality test. Values are given as mean ± SEM (n=9).
Figure 1.
SH-SY5Y differentiation into neurons. a. Timetable of SH-SY5Y human neuroblastoma cell line differentiation into neurons. b. The photographs of undifferentiated and differentiated SH-SY5Y neuroblastoma cells immunostained with the neuronal marker MAP-2 antibody. Alexa-488 was used to visualize MAP-2. Fluorescence images were obtained using a 20x objective. Scale bar: 100 µm. c. Expression of synaptophysin in SH-SY5Y cells and human neurons derived from SH-SY5Y cells analyzed by RT-qPCR. Data is presented as fold change in mRNA expression. The unpaired t-test test was followed by the D’Agostino and Pearson normality test. Values are given as mean ± SEM (****, p < 0.001; n=12). d. The overall expression level of the MMP-9 gene in SH-SY5Y cells and in human neurons. Data is presented as fold change in mRNA expression. The unpaired t-test test was followed by the D’Agostino and Pearson normality test. Values are given as mean ± SEM (n=9).
Luciferase reporter assay
DNA fragments corresponding to the region containing MMP-9-1562C/T SNP were amplified by PCR using heterozygous genomic DNA as a template and then cloned into KpnI and XhoI sites of the pGL3-promoter vector (Promega). DNA sequence analysis determined the insert's orientation and allele identity. The undifferentiated and differentiated SH-SY5Y cells were transfected with 200ng of appropriate construct using Lipofectamine TM 3000 Reagent. After 24h, the Dual-Glo® Luciferase Assay (Promega) was performed according to the manufacturer’s instructions.
Western blot
The nuclear or cytoplasmatic extracts (10 ug per well) were separated on 12.5% SDS-PAGE in Tris-glycine running buffer (25 mM Tris, 250 mM glycine (pH 8.3), 0.1% SDS) and electroblotted onto polyvinylidene difluoride membrane (Millipore) in transfer buffer (48 mM Tris, 39 mM glycine, 0.037% SDS, 20% methanol). Membranes were blocked in 5% nonfat milk/TBST (25 mm Tris-HCl (pH 8.0), 125 mm NaCl, 0.1% Tween 20) for 2 h at RT. Blots were then incubated with an appropriate primary antibody: TBP (TATA) Monoclonal Antibody (51841) 1: 1000 [Invitrogen #MA5-14739]; GAPDH 1:5000 [Millipore #MAB374] at 4 °C overnight. Peroxidase labeled Anti-mouse IgG secondary antibody [Vector Laboratories # PI-2000] was diluted 1:10000 in blocking buffer and incubated for over nihgt at 4°C. Proteins were detected by Western Bright Quantum HRP Substrate [Advansta #K-12042-D10] and chemiluminescence was measured using a UVITEC (Alliance) apparatus.
Gel shift assay
Gel shift assays were performed using LightShift® Chemiluminescent EMSA Kit (Thermo Scientific). Briefly, 5-10 µg of protein extracts prepared as described above were incubated with 20 fmol of 33-bp nucleotides labeled with biotin and 33-bp non-labeled competitor MMP-9-1562-C or –T (Supplementary Table1). Supershift assays were carried out by incubating 10 µg nuclear protein extracts with 1 µl of HDAC1 (D5C6U) XP® Rabbit mAb [Cell Signaling Technology #34589] or 0,5 µl of IgG normal Rabbit Ab [Santa Cruze sc-2027 X] (isotype control) for 30 min at RT. Then, DNA probes were added and incubated for 20 min at RT. Reactions with 5 µl ZNF384 Rabbit Ab [Invitrogen #PA-52044] or 0,5 µl of IgG normal Rabbit Ab (isotype control) were run at 4°C. The reactions with competitors (without Ab`s) were incubated for 20 min at RT before the electrophoretic run. DNA-protein and DNA protein-antibody complexes were resolved in 8% polyacrylamide gels in 0,5x TBE buffer (5xTBE: 450mM Tris, 450mM Boric Acid, 10 mM EDTA; pH = 8,3) for 1h with 100V at RT and transferred to Biodyne™ B Nylon Membrane (Thermo Scientific #77016) in 0,5x TBE buffer, on ice with 100V for 30 min. The complexes were visualized using UVITEC (Alliance).
Pull-down assay
60 pmol of 33 bp biotinylated nucleotides containing SNP of MMP-9-1562C/T (Supplementary Table 1.) were bound to 1mg of Dynabeads™ MyOne™ Streptavidin C1 (Thermo Fischer Scientific). 25 µl of SH-SY5Y nuclear extract was precleared in 1 ml of binding buffer (20 mM HEPES pH 7.6, 1mM EDTA, 10 mM (NH4)2SO4, 1mMDTT, 0,2% Tween-20, and 30mM KCl) in the presence of 1 mg of unbound beads and 288pmol of 33 bp non-labeled competitor MMP-9-1562-C or –T (Supplementary Table 1). Following magnetic separation, the supernatant was incubated with 1 mg of probe bound Dynabeads for 30 min at 4°C. The beads were then washed five times in 1 ml of binding buffer, frozen at -20°C and sent for mass spectrometry analysis.
Sample preparation for mass spectrometry analysis
Mass spectrometry experiments were performed in the Mass Spectrometry Laboratory at the Institute of Biochemistry and Biophysics PAS. Proteins were reduced by 1-hour incubation with 5 mM tris(2-carboxyethyl) phosphine (TCEP) at 60oC, followed by 10 min incubation at RT with 8 mM methyl methanethiosulfonate (MMTS). Proteins were digested on beads with 0.5 ugs of trypsin/Lys-C mix (Promega) overnight at 37 °C. Supernatants were transferred to a fresh tube and beads were washed with the extraction solution (2% acetonitrile, 0.1% formic acid) to assure maximum peptide recovery. Aliquots were dried on SpeedVac and reconstituted in 40 ul 10 mM HEPES buffer pH 8.0. Next steps were performed following SP3 protocol [Ultrasensitive proteome analysis using paramagnetic bead technology], with some modifications. Briefly, SP3 magnetic beads mix was prepared by combining equal parts of Sera-Mag Carboxyl hydrophilic and hydrophobic particles (09-981-121 and 09-981-123, GE Healthcare). The mix was washed three times with MS-grade water and resuspended in a working concentration of 10 ug/ul. Samples were mixed with 16 μl of the prepared bead mix and 1 ml of acetonitrile containing 0.1% formic acid to facilitate peptide binding. Beads were rinsed two times with 1 ml of isopropanol followed by acetonitrile. Peptides were eluted from the beads by subsequent incubation with MS-grade water and 2% DMSO with sonication during each step. Pulled aliquots were dried in SpeedVac and resuspended in 50 ul 2% acetonitrile and 0.1% formic acid. Peptide concentrations were measured using Pierce Quantitative Colorimetric Peptide Assay (Thermo Scientific).
Mass spectrometry and data analysis
1 ug of peptides from each sample was analyzed using nanoAcquity UPLC (Waters) directly coupled to a QExactive mass spectrometer (Thermo Scientific). Peptides were collected on a C18 trap column (180 µm x 20 mm, Waters) with 0.1% FA in water as a mobile phase and transferred to a nanocavity BEH C18 column (75 µm x 250 mm, 1.7 µm, Waters) using ACN gradient (0 – 35% ACN in 160 min) in the presence of 0.1% FA at a flow rate of 250 ml/min. Measurements were performed in the data-dependent mode with the top 12 precursors selected for MS2. Full MS scans covering the range of 300–1650 m/z were acquired at a resolution of 70,000 with a maximum injection time of 60 ms and an AGC target value of 1e6. MS2 scans were acquired at the resolution of 17,500 and the AGC target value of 5e5. Dynamic exclusion was set to 30 s. Obtained data were pre-processed with Mascot Distiller software (Matrixscience), and protein identification was performed using Mascot Server 2.5 (Matrixscience) against the Homo sapiens protein sequences (20496 sequences) deposited in the Swiss-Prot database (201904, 560118 sequences; 201292445 residues). The parameters settings were as follows: enzyme – Trypsin, missed cleavages – 2, fixed modifications – Methylthio (C), variable modifications –Oxidation(M), instrument – HCD. To reduce mass errors, peptide and fragment mass tolerance settings were established separately for each file after an off-line mass recalibration [
22]. The confidence assessment was based on the target/decoy database search strategy [
23], which provided q-value estimates for each peptide spectrum match. All queries with values> 0.01, subset proteins, and proteins identified with one peptide were discarded from further analysis. The mass recalibration, FDR computations, and data filtering were done with Mscan software and developed in-house (
http://proteom.ibb.waw.pl/mscan/)
HDAC1 and ZNF384 silencing with DsiRNAs
Dicer-Substrate Short Interfering RNAs (DsiRNAs) and TriFECTa® Kits (Integrated DNA Technologies) were used for RNA silencing. For transfection, SH-SY5Y cells were differentiated as described above. On the fifth day, cells were transferred to 12-well dishes at a density of 450,000 cells/well in the presence of ATRA and antibiotics. The next day, medium was changed to a fresh one without antibiotics. Cells were transfected using the X-tremeGENE siRNA Transfection Reagent (Roche) according to manufacturer’s protocol. 10 nM DsiRNA against HDAC1, ZNF384, positive and negative controls, and 0,5ug of the pGL3 vector with the promoter of the MMP-9 gene containing the C or T allele, respectively, were added to the medium. Cells were incubated in the transfection mixture for about 6.5 h, followed by a change of medium with the addition of ATRA, BDNF, and an antibio tics. On the seventh day the medium was changed again to the same medium as above (1:1 EMEM/F12 +ATRA+BDNF + antibiotics) On the eighth day, ca. 48h after transfection, the cells were flooded with TRIZOL, and the dishes were placed at -80ºC until the subsequent RNA isolation.
RNA isolation and reverse transcription quantitative PCR
Total RNA was extracted from SH-SY5Y cells using TRIzol® Reagent [Life Technologies # 15596018] according to the manufacturer's protocol. 500 ng of RNA samples were subjected to RT reaction using TaqMan Reverse Transcription Kit [Applied Biosystem #N8080234]. Real-time PCR analysis was performed with 5× HOT FIREPol EvaGreen qPCR Mix Plus [Solis BioDyne # 08-25-00020] using the LightCycler480 System (Roche). All experiments were performed in triplicate. CT values were chosen in the linear range of amplification, and the comparative CT method was used to calculate differences in gene expression. Primers and conditions for RT-qPCR are listed in Supplementary Table 2.
Luciferase assay after HDAC1 and ZNF384 silencing
48h after transfection, luciferase activity was measured using the Dual-Glo® Luciferase Assay System (Promega) according to manufactured protocol. Each luciferase construct was assayed in triplicate, and each transfection experiment was repeated at least three times.
4. Discussion
Earlier studies have revealed the implication of MMP-9-1562C/T polymorphism in the pathophysiology of various diseases related to the central nervous system including brain [
4]. Robert & Pelletier reported that the presence of SNPs in transcriptional regulatory elements influence the mRNA expression of this gene [
33]. Therefore, the cytosine (C)>thymidine (T) single nucleotide polymorphism (SNP) at position -1562 in the MMP-9 promoter may affect the expression level of MMP-9 mRNA and consequently susceptibility to various diseases related to the incorrect expression of this protease.
This is the first study to examine the role of MMP-9-1562C/T gene promoter polymorphism in the expression of MMP-9 in brain cells and related transcriptional regulation. Here we showed that MMP-9 promoter activity was almost two-fold higher in the presence of the C allele as compared to control (human neurons transfected with empty pGl3 vector). In case of the T allele the MMP-9 promoter activity was 1.4-fold higher as compared to control although the differences were not statistically significant. Similarly, the mRNA expression of MMP-9 varied depending on MMP-9-1562C/T polymorphism. MMP-9 mRNA level was significantly higher when the C allele was present in human neurons as compared to control. The expression of MMP-9 mRNA increased also in human neurons transfected with a pGl3 vector containing the T allele as compared to control, even though the differences were not statistically significant. Experiments performed by Zhang et al., however, have produced conflicting results. They showed that in MALU cells, the reporter activity of the MMP-9 promoter was almost 1.5-times greater under the control of the T allele than the reporter activity controlled by the C allele [
1]. Interestingly the allele-dependent effect of the MMP-9-1562C/T polymorphism was not shown in the studies conducted on primary amnion epithelial cell cultures, WISH, or THP-1 cells [
34]. The discrepancy between their studies and ours in this respect is puzzling Thus, we speculate that, that mechanism responsible for the differences in the allele-dependent MMP-9 promoter activity is perhaps cell type-dependent.
Next, we investigated differences in transcriptional regulatory proteins binding to the MMP-9 gene promoter in relation to the C or T allele. Using the electrophoretic mobility shift assay, we showed two bands in case of allele C, which was competed by the non-labeled allele C but not by the allele T. In comparison, experiment with the allele T showed one shifted band, competed by non-labeled allele T but not by the allele C in purified SH-SY5Y nuclear extracts. Similar conclusion was reached from the experiment with neuronal nuclear extract showing three bands in case of allele C, which were competed by non-labeled allele C and one shifted band in allele T. The results suggest that different sets of proteins regulate the transcriptional mechanism in the promoter of the MMP-9 gene in an allele-specific manner. The results are similar to those obtained by Zhang et al. who showed that the bands present on the C allele were less prominent in the T allele [
1].
We have extended our experiment further by identifying MM9-1562C/T binding proteins in human neurons using the pull-down assay followed by mass spectrometry technique. The analysis revealed 554 proteins which were bound to the C allele in the MMP-9-1562C/T polymorphism of the MMP-9 gene, and 586 proteins bound to the T allele, of which 542 were shared. After mass recalibration, FDR computations, and data filtering with Mscan software, we have selected proteins for further experiments. Additionally, we performed other bioinformatics analyzes that helped us find transcription factors that probably interact with the MMP-9 gene promoter, such as analysis in MatInpector and String software. For the C allele in the MMP-9 gene polymorphism, the Chromobox protein homolog 3 (CBX3), Nuclear factor 1 X-type (NFIX), and Histone deacetylase 1 (HDAC1) were chosen. These proteins as shown in Table 3b bound to the C allele in the presence of a competitor in the form of an MMP-9 promoter fragment with the T allele, so proteins that were not specific for the C allele were caught with the T allele. Reports have shown that CBX3 is involved in transcriptional silencing in heterochromatin-like complexes where recognizes and binds histone H3 tails methylated at 'Lys-9', leading to epigenetic repression [
35]. It is also a repressor in mammalian cells when tethered to DNA [
36]. In turn, NFIX has DNA-binding transcription factor activity [
37,
38] and regulates RNA polymerase II transcription [
39], while HDAC1 is responsible for transcriptional repression and gene silencing [
40,
41]. On the other hand, for the T variant in the -1562C / T polymorphism of the MMP-9 gene, the following proteins were selected RuvB-like 2 (RUVBL2), Histone-lysine N-methyltransferase 2A (KMT2A) and Zinc finger protein 384 (ZNF384). It has been reported that RUVBL2 is a component of the NuA4 histone acetyltransferase complex, which is involved in transcriptional activation of selected genes mainly by acetylation of nucleosomal histones H4 and H2A [
42]. Moreover, positively regulates transcription by RNA polymerase II [
43]. Further, KMT2A in the MLL1/MLL complex mediates H3K4me, a specific tag for epigenetic transcriptional activation [
44,
45,
46], positively regulates transcription [
47]. The last one from the selected genes, ZNF 384 is not well studied, however it is described as a transcription factor that binds sequence-specific double-stranded DNA [
48].
To investigate whether indeed selected transcription factors might specifically bind to the MMP-9 promoter, we performed EMSA supershift in the presence of RUVBL2, KMT2A, CBX3, NFIX, HDAC1 and ZNF384 specific antibodies. We have discovered that antibodies against RUVBL2, KMT2A, CBX3, NFIX do not bind to the promoter fragment of the MMP-9 gene in both the C and T allele polymorphism (data not shown). However, we observed complexes formed with antibodies directed to HDAC1 and ZNF384 in the C and T allele when human neurons were the source of nuclear proteins. Interestingly, in case of the T allele, the supershifts were much weaker as compared to the C allele suggesting that in the presence of the T allele, HDAC1 and ZNF384 may bind to the promoter of the MMP-9 to a lesser extent and/or in the presence of other transcription factors.
The results described above prompted us to evaluate the influence of HDAC1 and ZNF384 depletion onto MMP-9 promoter activity in human neurons. We depleted HDAC1 and ZNF384 expression by siRNA for HDAC1 and ZNF384 and measured the MMP-9 activity by luciferase assay. The MMP-9 promoter activity was higher under the control of the T allele after both HDAC1 and ZNF384 silencing as compared to control, i.e., cells transfected with empty pGL3 vector, although the increase in promoter activity was not statistically significant. As for the C allele, after HDAC1 knocking-down the MMP-9 activity was slightly higher compared to control. Opposite, after ZNF384 silencing, MMP-9 activity was lower compared to control. Previous data showed that overexpression of HDAC1 downregulated MMP-9 promoter activity in HT1080 fibrosarcoma cells in a dose-dependent manner [
49]. Surprisingly, after silencing of both HDAC1 and ZNF 384, we have shown that the MMP-9 mRNA expression level significantly increased when the T allele was present in the promoter as compared to the control, i.e., cells transfected with the empty pGL3 vector. Although, MMP-9 expression was higher in the case of C allele compared to control, this increase was not statistically significant. The above results strongly suggests that HDAC1 and ZNF384 play an important role in the regulation of MMP-9 gene in the presence of the T allele. Previously performed experiment showed that knockdown of HDAC1 suppressed the expression and protein level of MMP-9 in human breast cancer MDA-MB-321, MCF-7 cells and in U251, T98G glioblastoma cell line [
50,
51]. Our data show that after ZNF384 and HDAC1 silencing, both promoter activity and MMP-9 gene expression levels increase significantly, but only in the presence of the T allele. In the case of the C allele, after HDAC1 and ZNF384 silencing, there are no major differences in promoter activity and gene expression levels of MMP-9. Data obtained by other researchers showed that when HDAC1 was silenced, promoter activity and MMP-9 gene expression decreased [
49,
50,
51]. In the case of the T allele, as shown in the EMSA supershift image, HDAC1 and ZNF384 transcription factors bind to the MMP-9 gene promoter in smaller amounts, so differences in the activity of the MMP-9 gene promoter and its expression may depend strictly on a specific allele. Unfortunately, no other studies showing the effect of silencing ZNF384 on the activity of the MMP-9 promoter, or its expression have been published.