ORIGINAL RESEARCH article

Front. Oncol., 10 December 2020

Sec. Molecular and Cellular Oncology

Volume 10 - 2020 | https://doi.org/10.3389/fonc.2020.576362

Identification of Novel MeCP2 Cancer-Associated Target Genes and Post-Translational Modifications

  • 1. Department of Immunology and Molecular Microbiology, Texas Tech University Health Sciences Center, Lubbock, TX, United States

  • 2. Center for Biotechnology & Genomics, Texas Tech University, Lubbock, TX, United States

  • 3. Department of Internal Medicine, Texas Tech University, Lubbock, TX, United States

  • 4. Department of Surgery, Texas Tech University, Lubbock, TX, United States

Abstract

Abnormal regulation of DNA methylation and its readers has been associated with a wide range of cellular dysfunction. Disruption of the normal function of DNA methylation readers contributes to cancer progression, neurodevelopmental disorders, autoimmune disease and other pathologies. One reader of DNA methylation known to be especially important is MeCP2. It acts a bridge and connects DNA methylation with histone modifications and regulates many gene targets contributing to various diseases; however, much remains unknown about how it contributes to cancer malignancy. We and others previously described novel MeCP2 post-translational regulation. We set out to test the hypothesis that MeCP2 would regulate novel genes linked with tumorigenesis and that MeCP2 is subject to additional post-translational regulation not previously identified. Herein we report novel genes bound and regulated by MeCP2 through MeCP2 ChIP-seq and RNA-seq analyses in two breast cancer cell lines representing different breast cancer subtypes. Through genomics analyses, we localize MeCP2 to novel gene targets and further define the full range of gene targets within breast cancer cell lines. We also further examine the scope of clinical and pre-clinical lysine deacetylase inhibitors (KDACi) that regulate MeCP2 post-translationally. Through proteomics analyses, we identify many additional novel acetylation sites, nine of which are mutated in Rett Syndrome. Our study provides important new insight into downstream targets of MeCP2 and provide the first comprehensive map of novel sites of acetylation associated with both pre-clinical and FDA-approved KDACi used in the clinic. This report examines a critical reader of DNA methylation and has important implications for understanding MeCP2 regulation in cancer models and identifying novel molecular targets associated with epigenetic therapies.

Introduction

Epigenetic dysregulation involving mutations or abnormal expression of DNA methylation readers has been associated with a broad spectrum of disorders that range from Rett Syndrome to human cancers (17), and alterations in both the writing and reading of epigenetic marks have been linked with tumor progression at every stage (812). Aberrant DNA methylation not only promotes disease progression but is targeted via therapeutics applied in the clinic (1315). Because of the prevalence of abnormal epigenetic changes in tumor progression (1622), exploitation of this property led to FDA approved “epigenetic” therapies (23, 24). Interestingly, DNA methylation readers, such as methyl-CpG-binding protein 2 (MeCP2), bridge DNA methylation and histone modifications by binding to methylated DNA and recruiting co-repressor proteins (2528). While both normal and abnormal DNA methylation is read by MeCP2, much remains unknown about its role and regulation in cancer-associated pathologies. MeCP2 was shown early on to have an affinity for 5-methylcytosine in the context of methylated CG dinucleotides (mCG) (29, 30) and methylated CH (mCH), where H = A/C/T. MeCP2 binds methylated cytosine (3133) and shows selectivity for mCG sequences with adjacent A/T sequences (34). However, it also binds to hydoxymethylated cytosine (31, 32, 3537). While more investigation is needed, MeCP2 binding to mCH has been primarily noted on mCA (3133, 35, 36, 38). Studies have also shown that MeCP2 binding in mouse brain is proportional to mCAC + mCG density wherein transcription is sensitive to MeCP2 occupancy (38). Additionally, MeCP2 regulates tumor suppressor genes (TSG) silencing, and serves as a critical bridge for histone methyltransferases (HMTs) (25), histone deacetylases (HDACs) (26, 28, 39), and other proteins that bind modified histones or that mediate nucleosome remodeling (27, 40, 41). Moreover, MeCP2 has been reported to be amplified in diverse cancer including human triple-negative breast cancers (TNBC), and it activates growth factor pathways targeted by activated Ras, MAPK and PI3K pathways (42). Novel interacting protein partners and gene targets in brain tissue have also been identified (43). These are the types of enigmatic and versatile properties of MeCP2 that have contributed to long-standing knowledge gaps. We previously reported that inhibition of SIRT1 triggers acetylation of endogenous MeCP2 at lysine (K171), a site that regulates MeCP2 interaction with HDAC1 and ATRX (44). These findings demonstrated that MeCP2 post-translational modifications (PTMs) can critically impact its function, yet few PTMs have been mapped despite the potential that they might affect substrate specificity (35, 45). This knowledge gap is especially important given reports demonstrating unique characteristics of MeCP2 domains in determining binding specificity (46, 47) and the impact of MeCP2 on chromatin-dependent regulation of epigenetic writers (48). In the present study, we have identified additional novel PTMs across the length of MeCP2 and target genes in cancer models. Our findings provide new insight on the versatile role of MeCP2 which is known to be critical in regulating gene imprinting (49), transcriptional activation and repression (50) in disparate conditions that range from autism to cancer (4, 51, 52).

Materials and Methods

Cell Lines

MDA-MB-468 (HTB-132), MCF7 (HTB-22), MCF10A (CRL-10317), MCF12F (CRL-10783), PC3 (CRL-1435), T47D (HTB-133), BT549 (HTB-122), and MDA-MB 231 (HTB-26) cell lines used in this manuscript were purchased from ATCC which utilizes STR technology for cell authentication. Cells were used at a low passage (<20) within 6 months or less after receipt or resuscitation. MDA-MB-468, T47D, and BT549 cells were cultured in RPMI 1640 (Gibco). MCF10A and MCF12F were cultured in HuMEC medium supplemented with HuMEC supplement kit (Gibco). PC3 cells were cultured in ATCC formulated F-12K media (ATCC). MCF7 cells were propagated in MEM while MDA-MB-231 cells were cultured in DMEM (Gibco). T47D and MCF7 cells were cultured in media supplemented with 0.1% insulin (Sigma). All cells were grown in culture media supplemented with 1% pen-strep and 10% fetal bovine serum from GIBCO at 37 ° in 5% CO2.

Plasmids

pCDNA3.1 (−) was used as the backbone and Hemagglutinin (HA)-tagged-MeCP2-WT-pCDNA3.1 (−) (encoding MeCP2 e2 isoform), HA-tagged-K135Q-MeCP2-pCDNA3.1 (−) and HA-tagged-K135R-MeCP2-pCDNA3.1 (−) were generated using outward PCR method.

Bioluminescent MDA-MB-468 Cells

The pGL4.50[luc2/CMV/Hygro] plasmid (E1310) which encodes the luciferase reporter gene luc2 (Photinus pyralis) was purchased from Promega. MDA-MB-468 cells were plated in a 6-well plate (Genesee) at the seeding density of 2 × 105 cells in order to reach 60% confluency at the time of transfection. Cells were transfected with 1 µg of the pGL4.50[luc2/CMV/Hygro] plasmid for 48 h. Stable transfectants were selected with 0.5 mg/ml hygromycin (Sigma H3274-100MG)-containing media which was replaced every 3–4 days until total selection was achieved. Bioluminescence was confirmed by In Vivo Imaging System (IVIS) in the presence of luciferin substrate (Promega VivoGlo™ Luciferin, In Vivo Grade P1041).

MeCP2 Stable Knock-Down and Clonal Selection

MDA-MB-468 cells stably expressing pGL4.50[luc2/CMV/Hygro] plasmid (Promega E1310) were plated at the seeding density 2 × 105 cells in order to reach 60% confluency at the time of transduction 48 h prior to infection and then infected with pLKO.1-puro based shRNA MISSION lentiviral transduction particles purchased from Sigma for MeCP2 (TRCN0000330971, TRCN0000330972) and Non-Targeting shRN control transduction particles (SHC002V). The transduction was enhanced with 5 μg/ml polybrene (Sigma Millipore) and 2× multiplicity of infection (MOI) viral particles was added to the media. After 24 h, culture media was replaced with fresh media for 2 days. Stable clones were selected with 6 μg/ml puromycin-containing media which was replaced every 3-4 days until selection was achieved and knockdown confirmed by Western blots and qPCR.

MeCP2 Stable Overexpression and Clonal Selection

MDA-MB-468 cells stably expressing pGL4.50[luc2/CMV/Hygro] plasmid with >90% knocked down of endogenous MeCP2 were plated at the seeding density 2 × 105 cells in order to reach 60% confluency at the time of transfection, 48 h prior to transfection and then transfected with 1 ug of the pCDNA3.1 (−) backbone, Hemagglutinin (HA)-tagged-MeCP2-WT (encoding MeCP2 e2 isoform), HA-tagged-K135Q-MeCP2 and HA-tagged-K135R-MeCP2 plasmids. The G418 disulfate salt solution (Sigma G8168-10ML) selection was started 48 h after transfection at a concentration of 0.4 mg/ml, and the G418 containing media was replaced every 3–4 days until total selection was achieved and overexpression confirmed by Western blots.

RT-PCR and qPCR

Total RNA was isolated using the Aurum™ Total RNA Mini Kit (Bio Rad) and 2 µg of RNA was used to produce cDNA via the SuperScript® III First-Strand Synthesis System for RT-PCR (Invitrogen). Intron-spanning primers designed for gene expression analysis are summarized in Table 1. All primers were validated by end-point PCR (RT-PCR), a minus reverse transcription control (−RT control) was included in all RT-PCR experiments. Equal amount of synthesized cDNA was used for qPCR using the Power UP SYBR Green (Thermofisher Scientific, A25778) and the CFX96 Real-Time System C1000 Touch Thermal Cycler (Bio Rad). β-actin gene expression was used as endogenous control for mRNA quantification, as is not a MeCp2 target gene in both cell lines studied and its expression didn’t change after the depletion of MeCP2 in RNA-seq analysis.

Table 1

Forward primerReverse Primer
qPCR
CDH1TGC CCA GAA AAT GAA AAA GGGTG TAT GTG GCA ATG CGT TC
EGFRCTT CCT CCC AGT GCC TGA ATACTC CGT GGT CAT GCT CCA AT
HDAC1CGA GAC GGG ATT GAT GAC GAACT TGG CGT GTC CTT TGA TAG T
HIPK3AAA AGG ACG ATC TGC CCC TGTCC AAA GTG CTG AAC CTG ACT
IL6CCA GGA GAA GAT TCC AAA GATGGA AGG TTC AGG TTG TTT TCT G
KDM3ATCA CAG GAG CCA CAG TAG GAAGA TCA TCA AAC CTG GAA GGC A
KDM3BCCC ACA CCA GGT TCA CAA TCT AGCC AAC CGC ATC TTT CAC TG
KMT2BCTC CGG AAG TGC ACC TTT GACCG TGG ATG GCT GAT CTG TAG
LANCL2GGC AGC AAA AGT GGA CCA AGCTC GCT GCC AAA TCA CAT CG
MALAT1GCC TGG AAG CTG AAA AAC GGTGG AA AAC GCC TCA ATC CCA
MAP1BGAG ATG CTG CCA ATG CCt CTCAG GGT CAT TCC CAC TCA CC
MECP2CCC ATC AAC ACG GAA GAA AAG TGCA GGG TGG GGT CAT CAT AC
MRPS17TTG GCG GAG GTG ACC AAAACA TGC TTT GCT CGT GGA AC
NEAT1CCC TGT GCT TCC GAC TTC ATCCC TGG CCT AGT GGA AAT GG
NUPR2AGC TTT ACG ACT GCC TGG ACCTT CGA ACA GGT CCT TCG GT
PSPHCGG CAT AAG GGA GCT GGT AAGGC TGC GTC TCA TCA AAA CC
SRA1AGC CCA CAA GTT TCC CAG TCGGC TTG AAA GCT CTT GCA CC

qPCR Primers used in the current study.

Western Blots

Protein extracts were generated using RIPA lysis buffer supplemented with protease inhibitor cocktail (Thermo Fisher Scientific). The protein concentration was measured with the BCA method. Approximately 50 μg of protein from each sample was loaded on NuPAGE™ 4–12% Bis-Tris Protein Gels (Thermo Fisher Scientific) and run at 175 V constant voltage. A constant voltage of 30 V was used for protein transfer onto polyvinylidene fluoride (PVDF) membranes (Millipore-Sigma). Blots were probed with rabbit anti-MeCP2 antibodies (1:1,000; Cell Signaling) and mouse anti-HA antibodies (1:2,000; Santa Cruz) overnight at 4°C. After three washes with tris-buffered saline and polysorbate 20 (TBST; Fisher scientific), blots were then incubated with anti-rabbit HRP conjugate secondary antibody (1:5,000) and anti-mouse HRP conjugate at room temperature for 1 h. After washing three times, chemiluminescence (Pierce ECL Western blotting substrate: Thermofisher Scientific, A25778) was then used to visualize protein bands. β-actin antibody (1:10,000; Santa Cruz) was used as control.

Immunofluorescence

About 1 × 105 cells were plated on coverslips 48 h prior and they were washed with PBS and then fixed with 4% paraformaldehyde for 15 min. After washing with PBS they were then permeabilized with 0.2% Triton X-100 for 20 min and blocked with 5% BSA for 30 min. Following that, cells were incubated with primary antibodies MeCP2 (Cell Signaling, 3456S), or HA (Cell Signaling, 3724S) for 1 h at room temperature, washed with PBS three times and then incubated in the dark with Phalloidin 568 and secondary antibodies ALEXA-488 goat anti rabbit conjugate for 1 h at room temperature. After three PBS washes the coverslips were mounted with ProLong® Gold Antifade Mountant with DAPI (Thermo Fisher Scientific), the slides were allowed to cure for 48 h and then examined under the Nikon T-1E scanning confocal microscope, with a 60× objective, and analyzed with NIS software.

Chromatin Immunoprecipitation (ChIP)

Cells were plated for 72 h and the media was changed 24 h prior to the experiment. The cells were subjected to 1% formaldehyde cross-linking (Sigma) for 8 min at room temperature. The cross-linking reaction was quenched by adding glycine (Sigma) to a final concentration of 0.125 M for 5 min at room temperature. The medium was then removed and cells were washed twice with cold PBS containing protease inhibitor cocktail (Thermo Fisher Scientific). Cells were scraped in PBS and pelleted. Pellet was resuspended in SDS Lysis buffer (50 mM Tris–HCl pH 8.0, 10mM 0.5M EDTA, and 1% SDS) with protease inhibitor cocktail and sonicated in a Diagenode Bioruptor 300 sonicator. Sonication conditions involved 20 cycles (30” ON/30” OFF) for MDA-MB-468 and 25 cycles for MCF7 cells. Sonication was evaluated for every experiment in a 1% agarose gel and chromatin fragments showed a distribution in the 100–800 bp range. The soluble chromatin fraction was collected and incubated for 2 h at 4°C with either MeCP2 (Sigma, M9317), or rabbit IgG (Sigma I5006). Dynabeads Protein A (Invitrogen) were washed and then added to the chromatin-antibody mixture and mutated for 2 h at 4°C. Beads were washed with a low salt wash buffer (0.1% SDS, 1% Triton X-100, 2 mM EDTA, 20 mM Tris–HCl pH 8.1, and 150 mM NaCl), a high salt wash buffer (0.1% SDS, 1% Triton X-100, 2 mM EDTA, 20 mM Tris–HCl pH 8.1, and 500 mM NaCl), and TE (1 mM EDTA and 10 mM Tris–HCl pH 8). Reverse-crosslinking was performed at 65°C overnight followed by treatment with RNaseA (Invitrogen) for 2 h at 37°C and proteinase K (Invitrogen) at 55°C for 2 h. The chromatin was eluted and purified using Qiaquick PCR purification kit (Qiagen) and subjected to ChIP-PCR to evaluate the occupancy of MeCP2 along the promoters of target genes. The primers listed in Table 2 were used for analysis.

Table 2

Forward primerReverse Primer
ChIP-PCR
ACOT2GGT TTT GCT GTG ATG GCT CTAAG GTA GTG TGT GTT GGG GTA G
AP1M2ACA GAG ATG TGC GGT GCA ATTC CAC CCT CAG CCA TTG AT
AQP1ACT TCA GCA ACC ACT GGG TAGATT TCC TGT CCT CTG GCT GTC
ATG4DTGT ACC GTG GGC TTC TAT GCCAC CTT TTC AGG GGT GGA CA
CCT6ATTG GTT CAT TGC GAC TAC CATCA CTT GAG GCC AGG AGT TT
CHCHD2AGT AAT GGC GTG ACC CAA TGTTGG TTG GAA TTG GGA ACT TGA TG
DKK1ATT GGC AGG AAC AGG ATG TGTGAG TGG AAT GAG GAA GGA TTT GT
DLX1GGT CTT CAT TTG TTG CTC GCAAT CCT TCC TGC GCC CTA AT
DNMT1CAA AAG GGG AAC CTT GTT CACCT GGG AGG AAG AAA TAG GG
DUX4TTG CCA TTT GTT CAC TCT GCTGA TCC TGG GGT TGA AAG TC
EGFRCCT TGG CAC CTT TCT ACT GCTGG AAA ATC GGC TTC AAA AC
EIF3GATT GGC TCA GGG CTA TTG GTGGC TTT GCC TTC ATC AGC TT
EYA2ATG AAG AAC AAG CCC CCG AGGGG GCT GAG GGT GCA TAA AT
FBXL7ATC CAG AGA GGT GCG TTT GGGAG TGG GTG CAG ACA GAC AA
GBASCCT CTT GGG AGG TCA ACA AAGCA TAC GTC ATC TTG CAT GG
GHSRTTC AGA GTG GAG AGC TTT CTC TGGTGA GCT GAC TAT CTC TCT CCC A
GRK7AAT GGT TGC TGG ACG AAC ACATC CCT GTG GAC AAT ACT GGT G
H2AFB2AGA TAG CAC ACT CAA CGC CCCTT GGC CGT CAG GTA CTC AA
HDAC1CTG AGC TAA ATC AGC ACC CGCCT CCC ACT GCC CTA CAT AGA
HIPK3TCC TTC CCG ACC TCA CAC AACT GGC ACT CTT TCA TGG TGG
HIST1H4FTGT TTG TCT TCG ATC ATG TCT GGGGC GTC CCG TAT CAC ATT CT
ICAM1ATT GTC CGG GAA ACT GGA CGACA ACA GGC GGT GAG GAT TG
ICAM3CAT GGT CCA GTG GGA AAG GTATA GGC TTG TAC GCC ATC CC
ICAM5ACT AGA CGG AAG TGG GAC AGAGTC AAG TTT CCA CCA CGC AG
ID2GAG CTG TGC GTG AAA TTG CTACC GCT TAT TCA GCC ACA CA
IL6AGC ATC CCT CCA CTG CAA AGGTG CCC ATG CTA CAT TTG CC
ILF3AGG CAG CTA CTC CTA CTC GAACAC CAC TTG TCC TCC TCC TAA
JARID2ACG AGT GTA TGG GTG AGT GCCCA TTG CAG CCA TTT GTC CC
KDM1AAAG CCA ACG GAC AAG CTG TAACA TCA CAT CAT CTC TAC CCT CA
KDM1BAGT TTG GAA AAC CTG CAA CAC TAGA GTA GGT GAT TTC GCT GGG
KDM2ATGC TTC TCA ATG TGC TCT CCAGCC AGG CTG AAA ACA CTT ACT T
KDM3BAAC TCC TTT GCT CTC AGC GTTCC AAA TCT TAC CTC CCC GTC
KDM4AGGG TCA AAG CAC TTG GGG ATGCT TCA CAG AGC AAC AAG GC
KMT2AACC ACC ATG TGA CTA TTG GAC TTACA GCT CTT ACA GCG AAC ACA
KMT2BAAC CCC ACC CAT TTC CCT GTTTGG GAG GCC AGG AAG TTG AA
KRI1TGA TAA ATG CGG GGG TCC TTTCC ATC CTA ATC CCT ACG CTG A
LANCL2CCA TTA ACT TGG GAG GCT GAGGA CTG CAA TGT CAC CAA TG
METTL7AGCT CTG TGG ATG TGG TGG TCCTC ACA CCC TTT CAC TCA CCG
MRPS17TAG GTG CCA AGG ATG GTT TCCTC CCA AAG GTC AAG GAT CA
NUPR1LAAA GCC TGC GGA ACT TCA TAGGA TGG TTT CGA TCT CCT GA
OXCT2TTG ATG TCG TCC ACC GTC AGGAC CTG GCG AAC TGG ATG AT
P2RY11CTG GTG GTT GAG TTC CTG GTGCT GAT GCA GGT GAT GAA GA
PPAN-P2RY11GAC ACT GTC TCT CCC CAC AGACCC AAT CTG GGG CGT TCA ATC
PSPHGTG CTT GAA GGT GGG TAG GACTG TGG ATT CTG CAA GAG CA
SEC61GCCT GGG TTC AAG CAA TTC TGCAC CCG AGG TCA GGA GTT TA
SIRT1AAG AAA GGC AGT CGG ACC ATGCT GAC CTA CAG TAA GCA CTC A
SIRT5AAC CAC AGA CCT GCC TGA GTTCC CTC TCC CAT CAG GGT AT
SLC44A2ACC CTA CTT CAT GTC GCC CTCA TGC ACC CCC AGT CTA CAT
SUMF2CAA CAC AGA CCC CCA TCT CTCAT GGC TCA CTA CAG CCT CA
TCEB3CCTC AGA AAT CGC CTC CTG TCGAG AGT GCT TCT GGG TTT GC
TCEB3CL2CTC AGA AAT CGC CTC CTG TCGGA GAG TGC TTC TGG GTT TG
VOPP1TTC CAC AGC ACT CCT CAC AGGAA TGA GGC AGC AGA AGT CC
VSTM2AAAG GTT GGA TGG GTT TTT CCACA CTG GCA AAT TCC GTT TC
WNT3AGGC ATG GGG AGG TAT GCA ATGAA TCT GGC CGT GTG CTT TG
WT1GCT TGA ATG AGT GGT TGG GGAACC GCT GAC ACT GTG CTT CTC
ZNF154TCT CCA GCA TCA CAT CAC GGTCC TCT CAG TTG GGG AGC TT
ZNF713AGT CAC AAA AAT CCA GAG CCC AAGA ACA GGC AGG AAT CCA TGA

ChIP-PCR Primers used in the current study.

ChIP Sequencing

For ChIP-Seq experiments, ChIP DNA was prepared as described above library preparation was followed by high throughput sequencing with Illumina Hi-seq 2000 at GENEWIZ Corporation.

RNA Sequencing

RNA was prepared as described above, and library preparation and sequencing were performed at Center for Biotechnology & Genomics of Texas Tech University. RNA quality was determined using RNA Screen Tpe (Agilent). Ribosomal RNA depletion was achieved using NEB Next rRNA Depletion Kit (Human/Mouse/Rat) (NEB # E6310X). RNA fragmentation, double stranded cDNA and adaptor ligation was generated using NEBNext Ultra II Directional RNA Library Prep according to the manufacturer’s protocol (NEB # E7760L). PCR enriched libraries were quantified by Qubit and equimolar indexed libraries (different samples had different indexes for multiplexing) were pooled. Pooled libraries were quantitatively checked using the Agilent Tapestation 2200 and quantified using Qubit. The libraries were then diluted to 200 pM and spiked with 2% phiX libraries (Illumina control). The transcriptome sequencing was performed on the barcoded stranded RNA-Seq libraries using Illumina NovaSeq 6000 SP flow cell, paired-end reads (2 × 50 bp).

ChIP Sequencing and RNA Sequencing Data Analysis

For ChIP-Seq analysis, the FASTQ files were analyzed using DNASTAR’s Laser Gene software. MEME-ChIP was used to analyze MeCP2 binding motifs and TOMTOM to identify if those motifs were similar to known consensus sequences using the MEME Suite Programs http://meme-suite.org/index.html (53). We downloaded the FASTQ data sets of RRBS for MCF7 cells from the ENCODE portal (54) (https://www.encodeproject.org/) with the following identifiers: ENCSR943EFS, and ENCSR939RXT; then avisualized with Integrative Genomics Viewer (IGV). Venn diagrams to identify the overlapping genes were generated using the Venny tool https://bioinfogp.cnb.csic.es/tools/venny/index.html. For RNA-Seq analysis, the RNA-Seq reads were normalized by RPKM and assembled by mapping reads directly to the annotated human reference genome using the DNASTAR SeqMan software (DNASTAR, Inc., Madison, WI). Differential gene expression levels were quantified using Fisher’s Exact Test Signal Search in the DNASTAR ArrayStar software package (DNASTAR, Inc., Madison, WI). Differentially expressed genes were filtered if they met the criterion for a two-fold change, a p-value that was less than.05 at a 95% confidence interval. For each comparison, genes were sorted based on fold change, from low to high. The results were ported into Excel spreadsheets where the log2 of the fold change for each gene was calculated.

RNA Analysis In Silico

Relative RNA expression of 20 selected genes in breast cancer and normal adjacent tissue was downloaded from UCSC Xena platform on 11th of April 2020 (1,092 breast cancer primary tumors and 114 normal tissues).

Liquid Chromatography/Mass Spectrometry (LC–MS/MS)

PC3 and MDA-MB-468 cells were cultured and seeded in p150 mm dishes at 37°C under atmospheric oxygen conditions. Once 70% confluent, cells were treated with DMSO, 2µM panobinostat, 10 µM Inhibitor-IV, 10 µM Inhibitor-VII, and 10 µM pracinostat for 45 min to 1.5 h and harvested in RIPA buffer (with complete protease inhibitor cocktail, 1 µM Trichostatin A and 1 mM nicotinamide). Protein concentration was quantified by the BCA method. Immunoprecipitation was performed using 4 μg of anti-MeCP2 antibody (Cell Signaling) and incubated for 2 h at 4°C. Protein A dynabeads (Invitrogen) were added to the immune-complex and incubated for 2 h at 4°C. IP protocol was followed as mentioned above. Beads were washed with RIPA buffer (four times) and autoclaved water (two times). Dry beads were shipped to Applied Biomics Inc. (Hayward, CA) for acetylation site identification by LC–MS/MS mass spectrometry on a fee-based service. The specific lysine residues that were acetylated, exhibited ion peaks at mass/charge (m/z) ratio of ~126 as summarized in Figure 4A.

Statistical Analysis

Statistical analysis was performed using unpaired student’s t tests (Graph Pad Prism software) to assess whether differences observed in the various experiments were significant. All results are expressed as mean ± SEM and considered significant at *p <0.05, **p <0.01 and ***p <0.001.

Results

MeCP2 Binds Novel Genes in Breast Cancer Cells Associated With Diverse Biological Functions

Since the discovery that MeCP2 regulates transcription and mutations in the gene cause Rett Syndrome, there has been considerable interest in what regulates its function and what downstream genes are targeted (55, 56). DNA methylation and its readers influence transcription activation and repression in a context-dependent manner depending on the genomic location of binding (57, 58). While this process is known to be frequently altered in cancers (5961), many unknowns remain regarding the role of MeCP2 in regulating gene expression. Given abnormal DNA methylation in breast cancers (10, 12, 62, 63) and MeCP2 amplification in cancers (41), we wanted to identify new MeCP2 target genes and map novel sites of MeCP2 post-translational acetylation in breast cancer cells. We first examined MeCP2 protein expression in breast cancer cells and noted a range of expression across all lines with higher expression (64, 65) in MDA-MB-468 and BT-549 cells (Figure 1A). Both of these lines are derived from triple negative breast cancer (TNBC) which lack the expression of hormone receptors (ER and PR) and do not overexpress the growth factor receptor, HER2. To identify novel genomic targets of endogenous MeCP2, we performed MeCP2 ChIP-Seq analyses across two breast cancer cell lines (MCF7 and MDA-MB-468). These cells were chosen because they represent two different breast cancer subtypes and show relatively different MeCP2 protein expression levels. Also, inclusion of MCF7 in the ENCODE Project enabled comparison of our ChIP-Seq data with other publicly available data for epigenetic marks mapped in this cell line.

Figure 1

MeCP2 ChIP-Seq had not been done in MDA-MB-468 cells and our analysis revealed that MeCP2 binds to a wide spectrum of target genes (~20,000 in MDA-MB-468 and ~1,337 in MCF7 cells) ranging from miRNA, lncRNA, snRNA, processed and unprocessed pseudogenes, antisense and protein-coding genes. These genes are associated with a diverse range of cellular processes like gene expression, organization of the extracellular matrix, transport, or signal transduction, as shown in Figures 1B, C. In MeCP2 ChIP-Seq in MDA-MB-468 cells, we found that MeCP2 binds to multiple novel targets not previously associated with MeCP2 function in the context of breast cancer. Some of these included the following genes: a) SSU72 Homolog, RNA Polymerase II CTD Phosphatase (SSU72), a protein phosphatase that catalyzes the dephosphorylation of the C-terminal domain of RNA polymerase II (66); b) CAPN2 (Calpain 2), a calcium-sensitive cysteine protease (67); c) Plexin B2 (PLXNB2), a class B transmembrane receptor that participates in axon guidance and cell migration in response to semaphorins (68); d) Zinc Finger SWIM-Type Containing 4 (ZSWIM4); e) RUNX Family Transcription Factor 3 (RUNX3) a transcription factor that functions as a tumor suppressor and is frequently deleted or transcriptionally silenced in cancer (69, 70), and f) Solute Carrier Family 45 Member 4 (SLC45A4) (Figure 1B). Additionally, in MCF7 some of the notable genes included a) Ubiquitin Specific Peptidase 34 (USP34), a ubiquitin hydrolase that removes conjugated ubiquitin from AXIN1 and AXIN2, acting as a regulator of Wnt signaling pathway (71); b) Maltase–Glucoamylase (MGAM), an enzyme that plays a role in the digestion of starch (72); c) GDP-Mannose 4,6-Dehydratase (GMDS), an enzyme that participates in the synthesis of GDP-fucose from GDP-mannose (73); d) Solute Carrier Family 45 Member 4 (SLC45A4); e) CCDC26 Long Non-Coding RNA (CCDC26), a lncRNA class associated with Malignant Glioma and Astrocytoma (74, 75); and f) Sidekick Cell Adhesion Molecule 1 (SDK1) (Figure 1C). Moreover, 60% of the MCF7 loci (806 of 1,336) overlapped with MDA-MB-468 loci, including gene such as USP34, MGAM, GMDS, SLC45A4, CCDC26, and SDK1 (Figure 1D). We further analyzed the methylation status for genes in MCF7 cells for which publically available Reduced Representation Bisulfite Sequencing (RRBS) data was available (Figure S1). We found that MeCP2 binds to genes in MCF7 cells in regions where CpG methylation had been mapped such as SDK1, a cell adhesion molecule; Jagged 2 (JAG2), a Notch ligand; glycogenin 2 (GYG2), an enzyme involved in glycogen synthesis (Figure 1C, Figure S1C). These novel MeCP2 targets as well as others in Figure S1 had not previously been linked with MeCP2, but have been linked with pathobiology associated with cancer (7684) or genetic disorders such as Leigh syndrome (85) and Raine syndrome (86, 87). We also found that MeCP2 binds to genomic regions devoid of CpG methylation such as for USP34, MGAM, GMDS, SLC45A4, SSU72, CAPN2 and PLXN2 (Figures S1B–C). Similarly, while these are novel targets of MeCP2, many have been implicated in diverse cancers (67, 71, 8892). This further shows the complexity of MeCP2 binding across the genome. To identify the DNA motifs associated with MeCP2 genomic binding, we analyzed the genomic fragments sequenced in our MeCP2 ChIP-Seq analyses performed in triplicate. The MEME-ChIP analysis revealed a motif consistent across three independent experiments for both MDA-MB-468 and MCF7 cells (Figure 1E and Table 3).

Table 3

CONSENSUSWidthFragmentsE-value
MDA-MB-468
Replicate #1GTGATGTGTGYRTTCAACTCACAGAGTTGA3032328.1e-3515
Replicate #1CTAGACAGAAKMATTCTCAGAAACTT2635681.1e-2663
Replicate #1TTWCAYAGAGCAGWTTTGAAACACTCTTT2929236.5e-2645
Replicate #2TTGTGATGTGTGYRTTCAACTCACAGAGT2931021.4e-3473
Replicate #2AARCTAGACAGAAKVATTCTCAGAAA2632101.3e-2530
Replicate #2AACVTTYCTTTTCAYAGAGCAGWTTKGAAA3026855.0e-2508
Replicate #3TGWTGTGTGYDTTCAACTCACAGAGTTKAA3017629.7e-1782
Replicate #3TAGACAGAAKMATTCTCAGWAACTTCYTTG3013255.0e-1682
Replicate #3TYCTTTWCAYAGAGCAGWTTKGAAACACTC3010853.3e-1262
MCF7
Replicate #1CTTTGTGATGTGTGYRTTCAACTCACAGA308157.2e-2719
Replicate #1YTAGACAGARBARTTCTSARAMACTY2615302.0e-2663
Replicate #1GCAAGTGAKATTTVRACCKCTTTGAGGYC306584.7e-1616
Replicate #2CTTTGTGATGTGTGYRTTCAACTCACAGA296985.1e-2488
Replicate #2YTAGACAGARBARTTCTSARAMAC2412733.1e-2310
Replicate #2GMWTKGARKSSAATGGWRKRR2114287.5e-1387
Replicate #3TGTAATCCCAGCWMYTYGGGAGGCYG2668852.8e-1157
Replicate #3TTTTTKTWTTTTTWKTWGAGAC2273141.3e-709
Replicate #3GCCACCTCGCCCGGC1579151.6e-485

MeCP2 binding motifs.

MeCP2 Localizes to Novel Genes and Regulates Their Expression

We further determined the global occupancy of MeCP2 with respect to cellular functions and performed pathway analysis to identify the core pathways associated with the newly identified target genes in MDA-MB-468 and MCF7. We observed an enrichment of the gene expression, immune system, metabolism, metabolism of proteins, and signal transduction pathways (Figure 2A). We further randomly chose more than 100 genes identified in the triplicate analysis of MeCP2 ChIPseq in MDA-MB-468 cells and validated MeCP2 binding via MeCP2 ChIP-PCR, some of which are shown in Figure 2B. Consistent with our MeCP2 ChIPseq analyses, we found via MeCP2 ChIP-PCR that MeCP2 localizes to various gene promoters involved in diverse biological processes such as immune system regulation (IL6, ICAM3, and ICAM5), signal transduction (EGFR, WNT3A, and DKK1), transcription (KMT2A, SIRT1, HDAC1, DNMT1), developmental biology (DUX4) and lncRNAs (MALAT-1 and NEAT1) in MDA-MB-468 cells (Figure 2B). Several of the lncRNA targets identified are poorly studied, so we examined transcript expression patterns of some associated with MeCP2 MALAT1 and NEAT1 and established expression patterns across a panel of breast cancer cells. To determine whether MeCP2 depletion would lead to a change in expression of novel gene targets, MeCP2 was stably depleted with two different shRNA (sh1 and sh3) in MDA-MB-468 cells (Figure S2A). We also observed by quantitative RT-qPCR a change in mRNA expression of novel targets in which were validated for knockdown (Figure S2A). A minimum of three independent experiments showed that depletion of MeCP2 caused a change in the expression of several of the genes whose promoter it bound. We found that knockdown of MeCP2 in MDA-MB-468 cells caused an increase in some genes such as IL6, KDM3B, HIPK3, KDM3A, EGFR, and KMT2B and a reduction others such NUPR1L (also known as NUPR2), METTL7A, PSPH, LANCL2, MRPS17 and HDAC1. (Figure 2C and Figure S2B). Together these results show the complexity of MeCP2-mediated regulation of gene expression.

Figure 2

MeCP2 Targets Genes With Differential Expression Between Breast Cancer and Normal Samples

To evaluate the global effect of MeCP2 on RNA expression we performed RNA-Seq in MDA-MB-468 cells. We analyzed three independent experiments of non-targeting control (NTC) versus sh1 MeCP2 and NTC versus sh3 MeCP2 and found changes in 899 genes and 875 genes, respectively (Figure 3A). Overlap of ChIP-Seq hits and RNA-Seq hits showed 175 potential transcriptional targets of MeCP2 (Figure 3B). A pathway enrichment analysis of these potential targets showed their participation in the immune system, metabolism, metabolism of proteins, and signal transduction, among other pathways (Table 4). Moreover, these genes were differentially expressed in normal vs. breast cancer tissue (Figure 3C), and several of these target genes have been previously reported to be tumor suppressors (9399) while others were reported to be oncogenes (100103) (Figure 3D).

Figure 3

Table 4

Pathway identifierPathway nameGenes
R-HSA-9612973AutophagyRRAGD
R-HSA-1640170Cell CycleDMGDH;MECP2;DACH1;PIM1;L1CAM;PHLDA1
R-HSA-8953897Cellular responses to external stimuliDMGDH;DACH1;NCF1;RRAGD;ADRA1B;CRYAB
R-HSA-1266738Developmental BiologyDMGDH;NRP2;ADGRV1;PLXND1;SEMA3E;IRS2;L1CAM;MECP2;DACH1;COL4A1;ITGA10;COL5A2;PI3;KRTAP3-1;EPHA2
R-HSA-1643685DiseaseJAG2;DMGDH;PDGFRA;FLII;PCDH8;IRS2;ISG15;FMO3;L1CAM;MECP2;NT5E;DACH1;TBXAS1;PIM1;GPC2;LRAT;GGT1;TMEM45A;MUC6
R-HSA-73894DNA RepairTNNT2;ISG15
R-HSA-1474244Extracellular matrix organizationMECP2;MMP25;COL4A1;ITGB4;COL7A1;ITGA10;COL5A2;P3H2;LTBP1
R-HSA-74160Gene expression (Transcription)DMGDH;MECP2;GRIN2A;DACH1;CAV1;IGFBP3;RRAGD;PCDH8;NALCN;ANKRD6
R-HSA-109582HemostasisPDGFRA;CAV1;TFPI2;L1CAM;F3;PSG9;CEACAM3;PSG8;ISLR;PSG5;PSG4;ITGA10;HABP4;TIMP3;SLC16A8;KIF1A;RAPGEF3;CD177
R-HSA-168256Immune SystemDMGDH;NCF1;MAST4;ITGB4;FBXO27;IRS2;CST6;MECP2;C4B;MMP25;DACH1;PIM1;SNCG;PI3;MUC6;CD177;PDGFRA;OSBPL6;IL34;LIF;ISG15;OLFM4;CEACAM3;BST2;RAPGEF3
R-HSA-1430728MetabolismARL9;DMGDH;CLIC5;NRP2;MOGAT2;WNK4;LTBP1;NALCN;NT5E;DACH1;GPC2;LRAT;ST3GAL6;SLC16A8;GGT1;UGT1A7;GGT2;OSBPL6;CAV1;CYP4B1;FMO3;KMO;TIAM2;AKR1B10;TBXAS1;RAPGEF3
R-HSA-392499Metabolism of proteinsDMGDH;OSBPL6;NCF1;FBXO27;IGFBP3;PCDH8;P3H2;LTBP1;NALCN;ZDBF2;MECP2;NT5E;LGALS1;DACH1;COL4A1;COL7A1;TNNT2;ATP13A4;CPE;ST3GAL6;FOLR1;TMEM45A;MUC6;RABL6
R-HSA-1852241Organelle biogenesis and maintenanceRABL6
R-HSA-9609507Protein localizationP3H2
R-HSA-162582Signal TransductionDMGDH;CLIC5;NRP2;NCF1;PLXND1;FLII;STK39;IRS2;GRIK2;ADRA1B;NALCN;MECP2;NT5E;DACH1;GPC2;LRAT;ST3GAL6;NTSR1;JAG2;PDGFRA;EDN1;CAV2;PDE6G;CAV1;GPR55;L1CAM;F3;TIAM2;FOSL1;AKR1B10;COL4A1;COL5A2;PPP1R1B;RRAGD;RAPGEF3;RABL6
R-HSA-382551Transport of small moleculesDMGDH;SLC15A1;DACH1;TTYH1;ATP10D;WNK4;ATP13A4;ATP13A5;SLC16A8;SLC12A7;NALCN
R-HSA-5653656Vesicle-mediated transportTBC1D2;MAST4;COL4A1;COL7A1;P3H2;HSF4;KIF1A;FOLR1;KIAA0319;RABL6

Pathway analysis of the intersecting genes of MeCP2 ChIP-Seq and RNA-Seq.

Endogenous MeCP2 Is Acetylated at Key Lysine Residues and KDI Further Influence Its Acetylation Patterns

We previously reported that MeCP2 undergoes acetylation on Lys-171 in both MCF7 and RKO cells. We further demonstrated that a K171 acetylation mimetic did not perturb binding to select gene targets, but it diminished interaction of MeCP2 with binding partners such as ATRX and HDAC1 in colorectal cancer cells (44). In vivo and in vitro studies have demonstrated the importance of MeCP2 post-translational regulation (45, 104107), yet little has been done to comprehensively map novel MeCP2 PTMs. In the current study we wanted to extend our analyses and provide a comprehensive map of post-translational acetylation in other cancer cell line models. In order to further understand how MeCP2 is post-translationally regulated in TNBC breast and prostate cancer cell lines, we systematically identified the specific lysines on endogenous MeCP2 where acetylation was induced upon lysine deacetylase inhibition (KDACi). We inhibited SIRT1, a class III lysine deacetylase, using 10 µM Inhibitor-IV or 10 µM Inhibitor-VII, as well as the class I/II/IV lysine deacetylases using 2 µM panobinostat and 10 µM pracinostat. Given the links between DNA methylation and/or aberrant expression of DNA methylation readers in prostate cancer (4, 8) and TNBC (12, 42), we focused on two model lines representing each cancer, PC3 and MDA-MB-468, respectively. Next, we performed immunoprecipitation of endogenous MeCP2 and analyzed the samples using LC–MS/MS. Figure 4A summarizes the specific lysine residues that were acetylated and exhibited ion peaks at mass/charge (m/z) ratio of ~126 under basal (vehicle control) and KDI-induced conditions (i.e., cells treated with panobinostat, Inhibitor-IV, Inhibitor-VII, and pracinostat) (also see Figure S3). The mass spectrometry analyses showed that endogenous MeCP2 was acetylated at eight lysine residues under basal conditions (i.e., vehicle control) with induction in acetylation on K417 with 2 µM panobinostat; K364, K417, K431, K435 with 10 µM Inhibitor-IV; K22, K24, K27, K210 with 10 µM Inhibitor-VII; and K12, K135, K144, K171, K233 with 10 µM pracinostat. We found changes in acetylation patterns induced by exposure to both pre-clinical KDIs such as SIRT1 inhibitors and pracinostat as well as an FDA-approved inhibitor, panobinostat, which is used in the clinic to treat leukemias and lymphomas (23, 108). Interestingly, some of the lysine residues detected as acetylation sites (K22 and K135) were also sites mutated in Rett Syndrome. Moreover, some of the lysine residues detected as acetylation sites (such as K135), have been previously reported as sites linked with ubiquitination (4). We found acetylated lysine residues across the length of the protein, including at the N-terminus, in the methyl-binding domain (MBD), in the intermediate domain (ID) and the transcriptional repression domain (TRD) as well at the C-terminus region (Figure 4B). Together, these results indicate that MeCP2 is acetylated under basal and KDI-induced conditions in multiple cancer cell lines.

Figure 4

Next, we wanted to determine the impact of K135 acetylation on MeCP2 subcellular localization. We chose to study this site since it is situated in a highly conserved MBD domain and is a residue mutated in Rett syndrome patients. In order to probe the functional significance of MeCP2 acetylation, we generated HA-tagged wild-type MeCP2, HA-tagged deacetylation mimetics (K135R), HA-tagged acetylation mimetics (K135Q). Once the mutations were confirmed by sequencing, we then transfected and selected MDA-MB-468 cells with the plasmids for stable expression. Overexpression of HA-tagged MeCP2 constructs was confirmed by protein expression of WT and point mutants in MDA-MB-468 cells (Figure S4B). Using immunofluorescence assays, we detected that HA-tagged wild-type MeCP2, deacetylation mutants (K135R), and acetylation mutants (K135Q) were mostly in the nucleus of stably expressing MDA-MB-468 cells (Figure 5). These data demonstrate that post-translational acetylation on K135 lysine residue does not alter MeCP2 sub-cellular localization and calls for future studies to examine the role of acetylation at this residue as well as others identified in this report.

Figure 5

Discussion

The present study provides valuable insight on two important fronts. First, we identify novel genes that are subject to MeCP2-mediated regulation. Second, we provide a comprehensive identification of novel sites of post-translational acetylation associated with different cancer types and in response to multiple classes of deacetylase inhibitors. Concerning genomic analyses, these findings are important because we identify novel MeCP2 target genes linked with tumor progression which were not previously linked with MeCP2. While global DNA hypomethylation frequently occurs during tumorigenesis (60, 61), the promoters of TSGs may undergo hypermethylation (109111) and these aberrant changes in both the marks and the enzymes that modify them are being intensively examined for novel therapies (112116). These epigenomic changes may instigate genomic instability or generate a heritable molecular signature, which enables tumor progression, so identification of novel genomic targets of MeCP2 is very important (117119).

Previous reports linked MeCP2 expression with ER status (3) and with BRCA1 promoter silencing (120), which provided further rationale for assessing genome-wide MeCP2 profiling in both MCF7 and MDA-MB-468 cells, which represent two subtypes of breast cancer. We found that MeCP2 binds to multiple regions of genes, including promoters, exons, and introns. These novel targets have been associated with a wide range of regulatory and signaling pathways. We found that there was an overlap of around 800 genes between the two cell lines, and there were distinct MeCP2 binding motif enrichments between both cell lines. We observed that not only did MeCP2 bind many novel gene targets, but its depletion also led to both increases and decreases in their corresponding RNA transcripts. This is especially important given that studies demonstrate MeCP2 binds to methylated cytosines and hydoxymethylated cytosines in mCH dinucleotides, a property wherein many unknowns remain (3133, 35, 36, 38, 121). We discovered that MeCP2 localizes at various gene promoters involved in diverse processes such as autophagy (ATG4D), immune cell regulation (IL6), chromatin organization (KDM3B, KDM2A, KMT2B, KMT2A, KDM1A, HDAC1, HIST1H4F), circadian clock (SIRT1), developmental biology (EGFR, DKK1, SUMF2), extracellular matrix organization (ICAM5, ICAM3, ICAM1) and metabolism (EIF3G, SLC44A2, SUMF2, OXCT2, ACOT2, PSPH).

Recently, MeCP2 was shown to be amplified in human tumors and can mimic the function of activated Ras in cancer models (42), and also acts as a critical bridge linking information encoded in methylated DNA to epigenetic regulators (40, 122). Although MeCP2 binds methyltransferases (25), co-repressors (123) histone deacetylases (26), chromatin modulators (41), and long noncoding RNAs (lncRNAs) (124, 125) and other epigenetic regulators (4), much remains unknown about what regulates these interactions and what regulates binding to mCG vs. mCA dinucleotides as well as methylation-independent binding (126). However, our previous report provided some of the first insight into the role of post-translational regulation of MeCP2 binding to co-repressor proteins. We found that K171 acetylation regulates MeCP2 interaction with HDAC1 and ATRX (44). It is worth noting that the severity of Rett syndrome is influenced by the location of the mutation. For example, the site of the mutation can strongly impact cognitive and psychomotor skills as well as neonatal encephalopathy and death. Some mutations have been shown to disrupt conserved AT-hook regions and cause differential localization of α-Thalassemia/Mental Retardation Syndrome X-Linked (ATRX) (23). ATRX is a critical SWI/SNF-like chromatin remodeler (127) and our previous report demonstrates that a novel PTM on MeCP2 regulates MeCP2-ATRX binding, which is a critical aspect of MeCP2 function (43). Our present study identified an additional 17 novel sites of post-translational MeCP2 acetylation in triple-negative breast cancer and prostate cancer cell lines. Notably, nine of these lysines (K12, K135, K144, K177, K210, K233, K289, K364, K352, and K417) have been shown to be mutated in patients with Rett Syndrome. We identified four of these sites in our previous study mapping MeCP2 acetylation in MCF7 breast cancer cells and RKO colon cancer cells (K22, K135, K171, and K289) (43). Based on previous findings one may reason that one or more of these novel PTMs may be influencing MeCP2 function in cancer progression. Another example of the impact of post-translational regulation comes from transgenic models involving single MeCP2 serine residues that undergo post-translational regulation which show distinct neurological defects (106, 128), and phosphorylation of specific serine residues is enriched at specific gene promoters (104). However, much less is known about the role of MeCP2 acetylation as a regulatory switch in any context. Our more thorough mapping of novel MeCP2 acetylation PTMs performed here is a first step in defining their functional significance which is beyond the scope of the present study. Based on MeCP2 acetylation patterns induced by the various pre-clinical or clinical lysine deacetylase inhibitors, it is likely that KDACi’s that target class I/II vs. class III HDACs will influence MeCP2 function in both common and distinct ways. Based on acetylation mapping one can also reason that MeCP2 interaction with different KDACs may lead to important role in cell-type-specific biology driven by unique acetylation patterns. We previously demonstrated that lysine acetylation serves as a regulatory switch in Wnt pathway signaling (129, 130) and cancer-associated steroidogenesis (131, 132). The current study provides yet another example of the scope of post-translational acetylation and may help explain how SIRT1 preferentially targets active (133, 134) vs. repressed genes (135) depending on its deacetylation of specific non-histone partners (136, 137). Future work may identify more factors involved in this SIRT1-MeCP2 regulatory network, and through such work, our understanding of the key molecular relationships in cancer may lead to deeper understanding of the mechanism of action of epigenetic therapies and KDAC inhibitors.

Funding

This work was supported by the National Institute of Health (CA155223) to KP and a Cancer Prevention and Research Institute of Texas (CPRIT), Recruitment of Rising Stars Award (RR410008) to KP.

Statements

Data availability statement

Sequences and processed ChIP-Seq and RNA-Seq data files were deposited in the NCBI Gene Expression Omnibus (GEO) database under accession number GSE160150 and the BioProject: PRJNA667107.

Author contributions

IC-P, DV, MS, SP, LC, DM, MF, and KP discussed and designed the experiments. IC-P, DV, MS, SP, LC, DM, and MF performed experiments and/or analyzed data. IC-P, MS, DM, and KP wrote and edited the paper with input from all authors. IC-P, MS, SP, LC, DM, MF, JN, SA, RL, FR, and KP reviewed and revised the paper. All authors contributed to the article and approved the submitted version.

Acknowledgments

Immunofluorescence images were generated in the Image Analysis Core Facility, supported in part by TTUHSC.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2020.576362/full#supplementary-material

Data Sheet 1

MeCP2 acetylation site identification by LC-MS/MS mass spectrometry in MDA-MB-468 cells.

Data Sheet 2

MeCP2 acetylation site identification by LC-MS/MS mass spectrometry in PC3 cells.

Supplementary Figure 1

MeCP2 binds to unmethylated DNA (A) An assembly of IgG (first row), MeCP2 (second row) ChIP-Seq, and RRBS (third row) data in MCF7 for USP34, MGAM, GMDS, SLC45A4, CCDC26 and SDK1 genes, visualized by IGV. Each column is 4,000 bp wide. The third rows show the gene nearest to the ChIP-Seq alignment including its location, and orientation. The medium thick dark lines are the UTRs of the gene and the thicker dark regions are exons followed by thin lines with arrows which are the introns. (B) Schematic representations of the top ChIP-Seq hits where the red box indicates the MeCp2 binding, the white circles indicated methylated regions identified by RRBS and the medium thick blue lines are the UTRs of the gene and the thicker blue regions are exons followed by thin lines which are the introns. (C) Table of MeCp2 ChIP-Seq hits summarizing the presence/absence of CG methylation identified by RRBS in MCF7 cells.

Supplementary Figure 2

Validation of MeCP2 knockdown. (A) Left panel. Western blot to evaluate the protein levels of MeCP2 in MDA-MB-468 (NTC, sh1 MeCP2, and sh3 MeCP2) cells. Right panel. RT-qPCR analysis to evaluate expression of MeCP2 in MDA-MB-468 (NTC, sh1 MeCP2, and sh3 MeCP2) cells. Transcript levels were normalized to actin transcript levels. (B) Representative of two-independent RT-qPCR-based analysis to evaluate expression changes of NUPR2, PSPH, LANCL2, MRPS17, HDAC1, KDM3B, HIPK3, KDM3A, EGFR and KMT2B genes in MDA-MB-468 (NTC and sh MeCP2) cells. Transcript levels were normalized to actin transcript levels.

Supplementary Figure 3

Pharmacological inhibition of lysine deacetylases and key lysine residues acetylated on endogenous MeCP2. Acetylation of MeCP2 detected by Western blotting. (A) PC3 cells were treated with deacetylase inhibitors: DMSO as vehicle control, and SIRT1/2 Inhibitor-IV for a short time period range from 10 min to 1:15 h. (B) MDA-MB-468 were treated with deacetylase inhibitors: DMSO as vehicle control, and SIRT1/2 Inhibitor-IV for a short time period range from 10 to 120 min. (C) MDA-MB-468 were treated with deacetylase inhibitors: DMSO as vehicle control, and with various doses of SIRT1/2 Inhibitor-IV. For all immunoprecipitations equal amount of protein were loaded for each immunoprecipitation set up using acetyl-lysine (Ac-K) antibody as per protocol. Acetylation of MeCP2 was detected by Western blotting along with positive control, whole cell extract (WCE) using MeCP2 specific antibody. Species-matched IgG was used as a negative control. IgG heavy chain (IgG Hc) was blotted for as a control for equal antibody loading for immunoprecipitation and GAPDH for WCE. (D) The table indicates putative lysine residues that were found to be acetylated on MeCP2 under basal condition (DMSO) and upon deacetylase inhibition using 2 µM panobinostat (PANO), 10 µM Inhibitor-IV (IV), 10 µM Inhibitor-VII (VII), and 10 µM pracinostat (PRAC) and showed ion peaks at mass/charge (m/z) ratio of ∼126 in PC3 and MDA-MB-468 cells.

Supplementary Figure 4

Expression profile of lncRNA across normal and breast cancer cell lines. (A) RNA samples were extracted and converted to cDNA by reverse transcriptase enzyme. RT-PCR was performed to determine the expression of MALAT-1, MEG3, NEAT-1, CDKN2B, GAS5, SRA1, MIR31HG LncRNAs, and Beta actin as positive control in MCF12F normal breast cells and MCF7, BT549, MDA-MB-468, MDA-MB-231 and T47D breast cancer cell lines. (B) Stable expression of empty vector (EV), HA-epitope tagged MeCP2 wild type (WT), HA-tagged deacetylation mutants (K to R), HA-tagged acetylation mutants (K to Q) on K135 lysine residues in MDA-MB-468 cells.

Abbreviations

ChIP, Chromatin immunoprecipitation; ChIP-Seq, Chromatin Immunoprecipitation Sequencing; ENCODE, Encyclopedia of DNA Elements; FDA, U.S. Food and Drug Administration; HDAC, histone deacetylase; HMT, histone methyltransferase; KDI and KDACi, lysine deacetylase inhibitors; lncRNA, long non-coding RNA; MBD, methyl-CpG binding domain; MeCP2, methyl-CpG-binding protein 2; NTC, non-targeting control; PDX, Patient-Derived Xenograft; PTM, post-translational modifications; RNA-Seq, RNA-Sequencing; sh RNA, short hairpin RNA; SIRT1, Sirtuin 1; TNBC, triple negative breast cancer; TRD, transcriptional repression domain.

References

  • 1

    AmirREVan den VeyverIBWanMTranCQFranckeUZoghbiHY. Rett syndrome is caused by mutations in X-linked MECP2, encoding methyl-CpG-binding protein 2. Nat Genet (1999) 23:185–8. doi: 10.1038/13810

  • 2

    BaubecTDefossezPA. Reading DNA Modifications. J Mol Biol (2020) S0022-2836(20):30096‐6. doi: 10.1016/j.jmb.2020.02.001

  • 3

    MullerHMFieglHGoebelGHubalekMMWidschwendterAMuller-HolznerEet al. MeCP2 and MBD2 expression in human neoplastic and non-neoplastic breast tissue and its association with oestrogen receptor status. Br J Cancer (2003) 89:1934–9. doi: 10.1038/sj.bjc.6601392

  • 4

    PandeySPruittK. Functional assessment of MeCP2 in Rett syndrome and cancers of breast, colon, and prostate. Biochem Cell Biol (2017) 95:368–78. doi: 10.1139/bcb-2016-0154

  • 5

    LaveryLAUreKWanYWLuoCTrostleAJWangWet al. Losing Dnmt3a dependent methylation in inhibitory neurons impairs neural function by a mechanism impacting Rett syndrome. Elife (2020) 9:e52981. doi: 10.7554/eLife.52981

  • 6

    MannersMTErtelATianYAjitSK. Genome-wide redistribution of MeCP2 in dorsal root ganglia after peripheral nerve injury. Epigenetics Chromatin (2016) 9:23. doi: 10.1186/s13072-016-0073-5

  • 7

    ZhaoLLiuYTongDQinYYangJXueMet al. MeCP2 Promotes Gastric Cancer Progression Through Regulating FOXF1/Wnt5a/beta-Catenin and MYOD1/Caspase-3 Signaling Pathways. EBioMedicine (2017) 16:87100. doi: 10.1016/j.ebiom.2017.01.021

  • 8

    BernardDGilJDumontPRizzoSMonteDQuatannensBet al. The methyl-CpG-binding protein MECP2 is required for prostate cancer cell growth. Oncogene (2006) 25:1358–66. doi: 10.1038/sj.onc.1209179

  • 9

    DuQLuuPLStirzakerCClarkSJ. Methyl-CpG-binding domain proteins: readers of the epigenome. Epigenomics (2015) 7:1051–73. doi: 10.2217/epi.15.39

  • 10

    MatheAWong-BrownMLockeWJStirzakerCBrayeSGForbesJFet al. DNA methylation profile of triple negative breast cancer-specific genes comparing lymph node positive patients to lymph node negative patients. Sci Rep (2016) 6:33435. doi: 10.1038/srep33435

  • 11

    NoberiniRRestelliniCSavoiaEORaimondiFGhianiLJodiceMGet al. Profiling of Epigenetic Features in Clinical Samples Reveals Novel Widespread Changes in Cancer. Cancers (Basel) (2019) 11:723. doi: 10.3390/cancers11050723

  • 12

    StirzakerCZotenkoESongJZQuWNairSSLockeWJet al. Methylome sequencing in triple-negative breast cancer reveals distinct methylation clusters with prognostic value. Nat Commun (2015) 6:5899. doi: 10.1038/ncomms6899

  • 13

    FandyTEJiemjitAThakarMRhodenPSuarezLGoreSD. Decitabine induces delayed reactive oxygen species (ROS) accumulation in leukemia cells and induces the expression of ROS generating enzymes. Clin Cancer Res (2014) 20:1249–58. doi: 10.1158/1078-0432.CCR-13-1453

  • 14

    O’ConnorOAFalchiLLueJKMarchiEKinahanCSawasAet al. Oral 5-azacytidine and romidepsin exhibit marked activity in patients with PTCL: a multicenter phase 1 study. Blood (2019) 134:1395–405. doi: 10.1182/blood.2019001285

  • 15

    SiJBoumberYAShuJQinTAhmedSHeRet al. Chromatin remodeling is required for gene reactivation after decitabine-mediated DNA hypomethylation. Cancer Res (2010) 70:6968–77. doi: 10.1158/0008-5472.CAN-09-4474

  • 16

    BallestarEPazMFValleLWeiSFragaMFEspadaJet al. Methyl-CpG binding proteins identify novel sites of epigenetic inactivation in human cancer. EMBO J (2003) 22:6335–45. doi: 10.1093/emboj/cdg604

  • 17

    CouronneLBastardCBernardOA. TET2 and DNMT3A mutations in human T-cell lymphoma. N Engl J Med (2012) 366:95–6. doi: 10.1056/NEJMc1111708

  • 18

    LeyTJDingLWalterMJMcLellanMDLamprechtTLarsonDEet al. DNMT3A mutations in acute myeloid leukemia. N Engl J Med (2010) 363:2424–33. doi: 10.1056/NEJMoa1005143

  • 19

    MunozDPLeeELTakayamaSCoppeJPHeoSJBoffelliDet al. Activation-induced cytidine deaminase (AID) is necessary for the epithelial-mesenchymal transition in mammary epithelial cells. Proc Natl Acad Sci USA (2013) 110:E2977–86. doi: 10.1073/pnas.1301021110

  • 20

    PruittK. Molecular and Cellular Changes During Cancer Progression Resulting From Genetic and Epigenetic Alterations. Prog Mol Biol Transl Sci (2016) 144:347. doi: 10.1016/bs.pmbts.2016.09.001

  • 21

    RasmussenKDHelinK. Role of TET enzymes in DNA methylation, development, and cancer. Genes Dev (2016) 30:733–50. doi: 10.1101/gad.276568.115

  • 22

    XuJWangYYDaiYJZhangWZhangWNXiongSMet al. DNMT3A Arg882 mutation drives chronic myelomonocytic leukemia through disturbing gene expression/DNA methylation in hematopoietic cells. Proc Natl Acad Sci USA (2014) 111:2620–5. doi: 10.1073/pnas.1400150111

  • 23

    Garnock-JonesKP. Panobinostat: first global approval. Drugs (2015) 75:695704. doi: 10.1007/s40265-015-0388-8

  • 24

    KaminskasEFarrellATWangYCSridharaRPazdurR. FDA drug approval summary: azacitidine (5-azacytidine, Vidaza) for injectable suspension. Oncologist (2005) 10:176–82. doi: 10.1634/theoncologist.10-3-176

  • 25

    FuksFHurdPJWolfDNanXBirdAPKouzaridesT. The methyl-CpG-binding protein MeCP2 links DNA methylation to histone methylation. J Biol Chem (2003) 278:4035–40. doi: 10.1074/jbc.M210256200

  • 26

    JonesPLVeenstraGJWadePAVermaakDKassSULandsbergerNet al. Methylated DNA and MeCP2 recruit histone deacetylase to repress transcription. Nat Genet (1998) 19:187–91. doi: 10.1038/561

  • 27

    NanXHouJMacleanANasirJLafuenteMJShuXet al. Interaction between chromatin proteins MECP2 and ATRX is disrupted by mutations that cause inherited mental retardation. Proc Natl Acad Sci USA (2007) 104:2709–14. doi: 10.1073/pnas.0608056104

  • 28

    SuzukiMYamadaTKihara-NegishiFSakuraiTOikawaT. Direct association between PU.1 and MeCP2 that recruits mSin3A-HDAC complex for PU.1-mediated transcriptional repression. Oncogene (2003) 22:8688–98. doi: 10.1038/sj.onc.1207182

  • 29

    MeehanRRLewisJDBirdAP. Characterization of MeCP2, a vertebrate DNA binding protein with affinity for methylated DNA. Nucleic Acids Res (1992) 20:5085–92. doi: 10.1093/nar/20.19.5085

  • 30

    NanXMeehanRRBirdA. Dissection of the methyl-CpG binding domain from the chromosomal protein MeCP2. Nucleic Acids Res (1993) 21:4886–92. doi: 10.1093/nar/21.21.4886

  • 31

    ChenLChenKLaveryLABakerSAShawCALiWet al. MeCP2 binds to non-CG methylated DNA as neurons mature, influencing transcription and the timing of onset for Rett syndrome. Proc Natl Acad Sci USA (2015) 112:5509–14. doi: 10.1073/pnas.1505909112

  • 32

    GabelHWKindeBStroudHGilbertCSHarminDAKastanNRet al. Disruption of DNA-methylation-dependent long gene repression in Rett syndrome. Nature (2015) 522:8993. doi: 10.1038/nature14319

  • 33

    GuoJUSuYShinJHShinJLiHXieBet al. Distribution, recognition and regulation of non-CpG methylation in the adult mammalian brain. Nat Neurosci (2014) 17:215–22. doi: 10.1038/nn.3607

  • 34

    KloseRJSarrafSASchmiedebergLMcDermottSMStanchevaIBirdAP. DNA binding selectivity of MeCP2 due to a requirement for A/T sequences adjacent to methyl-CpG. Mol Cell (2005) 19:667–78. doi: 10.1016/j.molcel.2005.07.021

  • 35

    MellenMAyataPDewellSKriaucionisSHeintzN. MeCP2 binds to 5hmC enriched within active genes and accessible chromatin in the nervous system. Cell (2012) 151:1417–30. doi: 10.1016/j.cell.2012.11.022

  • 36

    SperlazzaMJBilinovichSMSinananLMJavierFRWilliamsDCJr.Structural Basis of MeCP2 Distribution on Non-CpG Methylated and Hydroxymethylated DNA. J Mol Biol (2017) 429:1581–94. doi: 10.1016/j.jmb.2017.04.009

  • 37

    SpruijtCGGnerlichFSmitsAHPfaffenederTJansenPWBauerCet al. Dynamic readers for 5-(hydroxy)methylcytosine and its oxidized derivatives. Cell (2013) 152:1146–59. doi: 10.1016/j.cell.2013.02.004

  • 38

    LaggerSConnellyJCSchweikertGWebbSSelfridgeJRamsahoyeBHet al. MeCP2 recognizes cytosine methylated tri-nucleotide and di-nucleotide sequences to tune transcription in the mammalian brain. PloS Genet (2017) 13:e1006793. doi: 10.1371/journal.pgen.1006793

  • 39

    NanXNgHHJohnsonCALahertyCDTurnerBMEisenmanRNet al. Transcriptional repression by the methyl-CpG-binding protein MeCP2 involves a histone deacetylase complex. Nature (1998) 393:386–9. doi: 10.1038/30764

  • 40

    AgarwalNHardtTBreroANowakDRothbauerUBeckerAet al. MeCP2 interacts with HP1 and modulates its heterochromatin association during myogenic differentiation. Nucleic Acids Res (2007) 35:5402–8. doi: 10.1093/nar/gkm599

  • 41

    LeightonGWilliamsDCJr.The Methyl-CpG-Binding Domain 2 and 3 Proteins and Formation of the Nucleosome Remodeling and Deacetylase Complex. J Mol Biol (2019) 432:1624–39. doi: 10.1016/j.jmb.2019.10.007

  • 42

    NeupaneMClarkAPLandiniSBirkbakNJEklundACLimEet al. MECP2 Is a Frequently Amplified Oncogene with a Novel Epigenetic Mechanism That Mimics the Role of Activated RAS in Malignancy. Cancer Discovery (2016) 6:4558. doi: 10.1158/2159-8290.CD-15-0341

  • 43

    Martinez dePAKhajaviLMartinHClaveria-GimenoRTomDSCheemaMSet al. MeCP2-E1 isoform is a dynamically expressed, weakly DNA-bound protein with different protein and DNA interactions compared to MeCP2-E2. Epigenetics Chromatin (2019) 12:63. doi: 10.1186/s13072-019-0298-1

  • 44

    PandeySSimmonsGEJr.MalyarchukSCalhounTNPruittK. A novel MeCP2 acetylation site regulates interaction with ATRX and HDAC1. Genes Cancer (2015) 6:408–21. doi: 10.18632/genesandcancer.84

  • 45

    BelliniEPavesiGBarbieroIBergoAChandolaCNawazMSet al. MeCP2 post-translational modifications: a mechanism to control its involvement in synaptic plasticity and homeostasis? Front Cell Neurosci (2014) 8:236. doi: 10.3389/fncel.2014.00236

  • 46

    KhrapunovSTaoYChengHPadlanCHarrisRGalanopoulouASet al. MeCP2 Binding Cooperativity Inhibits DNA Modification-Specific Recognition. Biochemistry (2016) 55:4275–85. doi: 10.1021/acs.biochem.6b00451

  • 47

    KhrapunovSWarrenCChengHBerkoERGreallyJMBrenowitzM. Unusual characteristics of the DNA binding domain of epigenetic regulatory protein MeCP2 determine its binding specificity. Biochemistry (2014) 53:3379–91. doi: 10.1021/bi500424z

  • 48

    RajaveluALunguCEmperleMDukatzMBrohmABrocheJet al. Chromatin-dependent allosteric regulation of DNMT3A activity by MeCP2. Nucleic Acids Res (2018) 46:9044–56. doi: 10.1093/nar/gky715

  • 49

    DrewellRAGoddardCJThomasJOSuraniMA. Methylation-dependent silencing at the H19 imprinting control region by MeCP2. Nucleic Acids Res (2002) 30:1139–44. doi: 10.1093/nar/30.5.1139

  • 50

    ChahrourMJungSYShawCZhouXWongSTQinJet al. MeCP2, a key contributor to neurological disease, activates and represses transcription. Science (2008) 320:1224–9. doi: 10.1126/science.1153252

  • 51

    CrespiB. Autism and cancer risk. Autism Res (2011) 4:302–10. doi: 10.1002/aur.208

  • 52

    KhwajaOSHoEBarnesKVO’LearyHMPereiraLMFinkelsteinYet al. Safety, pharmacokinetics, and preliminary assessment of efficacy of mecasermin (recombinant human IGF-1) for the treatment of Rett syndrome. Proc Natl Acad Sci USA (2014) 111:4596–601. doi: 10.1073/pnas.1311141111

  • 53

    BaileyTLBodenMBuskeFAFrithMGrantCEClementiLet al. MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res (2009) 37:W202–8. doi: 10.1093/nar/gkp335

  • 54

    DavisCAHitzBCSloanCAChanETDavidsonJMGabdankIet al. The Encyclopedia of DNA elements (ENCODE): data portal update. Nucleic Acids Res (2018) 46:D794–801. doi: 10.1093/nar/gkx1081

  • 55

    HiteKCAdamsVHHansenJC. Recent advances in MeCP2 structure and function. Biochem Cell Biol (2009) 87:219–27. doi: 10.1139/O08-115

  • 56

    LiCJiangSLiuSQLykkenEZhaoLTSevillaJet al. MeCP2 enforces Foxp3 expression to promote regulatory T cells’ resilience to inflammation. Proc Natl Acad Sci U S A (2014) 111:E2807–16. doi: 10.1073/pnas.1401505111

  • 57

    SchubelerD. Function and information content of DNA methylation. Nature (2015) 517:321–6. doi: 10.1038/nature14192

  • 58

    SuraniMA. Imprinting and the initiation of gene silencing in the germ line. Cell (1998) 93:309–12. doi: 10.1016/S0092-8674(00)81156-3

  • 59

    BaylinSBHermanJG. DNA hypermethylation in tumorigenesis: epigenetics joins genetics. Trends Genet (2000) 16:168–74. doi: 10.1016/S0168-9525(99)01971-X

  • 60

    FeinbergAPGehrkeCWKuoKCEhrlichM. Reduced genomic 5-methylcytosine content in human colonic neoplasia. Cancer Res (1988) 48:1159–61.

  • 61

    GoelzSEVogelsteinBHamiltonSRFeinbergAP. Hypomethylation of DNA from benign and malignant human colon neoplasms. Science (1985) 228:187–90. doi: 10.1126/science.2579435

  • 62

    ChristensenBCKelseyKTZhengSHousemanEAMarsitCJWrenschMRet al. Breast cancer DNA methylation profiles are associated with tumor size and alcohol and folate intake. PloS Genet (2010) 6:e1001043. doi: 10.1371/journal.pgen.1001043

  • 63

    SkvortsovaKMasle-FarquharELuuPLSongJZQuWZotenkoEet al. DNA Hypermethylation Encroachment at CpG Island Borders in Cancer Is Predisposed by H3K4 Monomethylation Patterns. Cancer Cell (2019) 35:297314. doi: 10.1016/j.ccell.2019.01.004

  • 64

    Ito-IshidaAYamalanchiliHKShaoYBakerSAHeckmanLDLaveryLAet al. Genome-wide distribution of linker histone H1.0 is independent of MeCP2. Nat Neurosci (2018) 21:794–8. doi: 10.1038/s41593-018-0155-8

  • 65

    Lopez-SerraLEstellerM. Proteins that bind methylated DNA and human cancer: reading the wrong words. Br J Cancer (2008) 98:1881–5. doi: 10.1038/sj.bjc.6604374

  • 66

    SpectorBMTurekMEPriceDH. Functional interaction of human Ssu72 with RNA polymerase II complexes. PloS One (2019) 14:e0213598. doi: 10.1371/journal.pone.0213598

  • 67

    MiaoCLiangCTianYXuAZhuJZhaoKet al. Overexpression of CAPN2 promotes cell metastasis and proliferation via AKT/mTOR signaling in renal cell carcinoma. Oncotarget (2017) 8:97811–21. doi: 10.18632/oncotarget.22083

  • 68

    PerrotVVazquez-PradoJGutkindJS. Plexin B regulates Rho through the guanine nucleotide exchange factors leukemia-associated Rho GEF (LARG) and PDZ-RhoGEF. J Biol Chem (2002) 277:43115–20. doi: 10.1074/jbc.M206005200

  • 69

    KangKAZhangRKimGYBaeSCHyunJW. Epigenetic changes induced by oxidative stress in colorectal cancer cells: methylation of tumor suppressor RUNX3. Tumour Biol (2012) 33:403–12. doi: 10.1007/s13277-012-0322-6

  • 70

    ManandharSLeeYM. Emerging role of RUNX3 in the regulation of tumor microenvironment. BMB Rep (2018) 51:174–81. doi: 10.5483/BMBRep.2018.51.4.033

  • 71

    LuiTTLacroixCAhmedSMGoldenbergSJLeachCADaulatAMet al. The ubiquitin-specific protease USP34 regulates axin stability and Wnt/beta-catenin signaling. Mol Cell Biol (2011) 31:2053–65. doi: 10.1128/MCB.01094-10

  • 72

    Diaz-SotomayorMQuezada-CalvilloRAverySEChackoSKYanLKLinAHet al. Maltase-glucoamylase modulates gluconeogenesis and sucrase-isomaltase dominates starch digestion glucogenesis. J Pediatr Gastroenterol Nutr (2013) 57:704–12. doi: 10.1097/MPG.0b013e3182a27438

  • 73

    PfeifferMJohanssonCKrojerTKavanaghKLOppermannUNidetzkyB. A Parsimonious Mechanism of Sugar Dehydration by Human GDP-Mannose-4,6-dehydratase. ACS Catal (2019) 9:2962–8. doi: 10.1021/acscatal.9b00064

  • 74

    JenkinsRBWrenschMRJohnsonDFridleyBLDeckerPAXiaoYet al. Distinct germ line polymorphisms underlie glioma morphologic heterogeneity. Cancer Genet (2011) 204:13–8. doi: 10.1016/j.cancergencyto.2010.10.002

  • 75

    WangXLuoTRuanMLiuPWangSZhuW. Association of the CCDC26 rs4295627 polymorphism with the risk of glioma: Evidence from 7,290 cases and 11,630 controls. Mol Clin Oncol (2016) 4:878–82. doi: 10.3892/mco.2016.813

  • 76

    BaggerSOHopkinsonBMPandeyDPBakMBrydholmAVVilladsenRet al. Aggressiveness of non-EMT breast cancer cells relies on FBXO11 activity. Mol Cancer (2018) 17:171. doi: 10.1186/s12943-018-0918-6

  • 77

    BenvenutoGTodeschiniPParacchiniLCaluraEFruscioRRomaniCet al. Expression profiles of PRKG1, SDF2L1 and PPP1R12A are predictive and prognostic factors for therapy response and survival in high-grade serous ovarian cancer. Int J Cancer (2020) 147:565–74. doi: 10.1002/ijc.32935

  • 78

    ChungCYSunZMullokandovGBoschAQadeerZACihanEet al. Cbx8 Acts Non-canonically with Wdr5 to Promote Mammary Tumorigenesis. Cell Rep (2016) 16:472–86. doi: 10.1016/j.celrep.2016.06.002

  • 79

    EngqvistHParrisTZKovacsANemesSWernerREDeLSet al. Immunohistochemical validation of COL3A1, GPR158 and PITHD1 as prognostic biomarkers in early-stage ovarian carcinomas. BMC Cancer (2019) 19:928. doi: 10.1186/s12885-019-6084-4

  • 80

    KunduANamHShelarSChandrashekarDSBrinkleyGKarkiSet al. PRDM16 suppresses HIF-targeted gene expression in kidney cancer. J Exp Med (2020) 217(6):e20191005. doi: 10.1084/jem.20191005

  • 81

    TanJJonesMKosekiHNakayamaMMunteanAGMaillardIet al. CBX8, a polycomb group protein, is essential for MLL-AF9-induced leukemogenesis. Cancer Cell (2011) 20:563–75. doi: 10.1016/j.ccr.2011.09.008

  • 82

    ThurnerLHartmannSFadleNKemeleMBockTBewarderMet al. LRPAP1 is a frequent proliferation-inducing antigen of BCRs of mantle cell lymphomas and can be used for specific therapeutic targeting. Leukemia (2019) 33:148–58. doi: 10.1038/s41375-018-0182-1

  • 83

    VeroneARDuncanKGodoyAYadavNBakinAKoochekpourSet al. Androgen-responsive serum response factor target genes regulate prostate cancer cell migration. Carcinogenesis (2013) 34:1737–46. doi: 10.1093/carcin/bgt126

  • 84

    YokoyamaSHigashiMKitamotoSOeldorfMKnippschildUKornmannMet al. Aberrant methylation of MUC1 and MUC4 promoters are potential prognostic biomarkers for pancreatic ductal adenocarcinomas. Oncotarget (2016) 7:42553–65. doi: 10.18632/oncotarget.9924

  • 85

    ImagawaEOsakaHYamashitaAShiinaMTakahashiESugieHet al. A hemizygous GYG2 mutation and Leigh syndrome: a possible link? Hum Genet (2014) 133:225–34. doi: 10.1007/s00439-013-1372-6

  • 86

    Hernandez-ZavalaACortes-CamachoFPalmaLGodinez-AguilarREspinosa-GarciaAMPerez-DuranJet al. Two Novel FAM20C Variants in A Family with Raine Syndrome. Genes (Basel) (2020) 11(2):222. doi: 10.3390/genes11020222

  • 87

    HungCYRodriguezMRobertsABauerMMihalekIBodamerO. A novel FAM20C mutation causes a rare form of neonatal lethal Raine syndrome. Am J Med Genet A (2019) 179:1866–71. doi: 10.1002/ajmg.a.61291

  • 88

    LinCXiaJGuZMengYGaoDWeiS. Downregulation of USP34 Inhibits the Growth and Migration of Pancreatic Cancer Cells via Inhibiting the PRR11. Onco Targets Ther (2020) 13:1471–80. doi: 10.2147/OTT.S228857

  • 89

    Vincent-ChongVKAnwarAKaren-NgLPCheongSCYangYHPradeepPJet al. Genome wide analysis of chromosomal alterations in oral squamous cell carcinomas revealed over expression of MGAM and ADAM9. PloS One (2013) 8:e54705. doi: 10.1371/journal.pone.0054705

  • 90

    WeiXZhangKQinHZhuJQinQYuYet al. GMDS knockdown impairs cell proliferation and survival in human lung adenocarcinoma. BMC Cancer (2018) 18:600. doi: 10.1186/s12885-018-4524-1

  • 91

    XiangGChengY. MiR-126-3p inhibits ovarian cancer proliferation and invasion via targeting PLXNB2. Reprod Biol (2018) 18:218–24. doi: 10.1016/j.repbio.2018.07.005

  • 92

    XieXQWangMJLiYLeiLPWangNLvZYet al. miR-124 Intensified Oxaliplatin-Based Chemotherapy by Targeting CAPN2 in Colorectal Cancer. Mol Ther Oncolytics (2020) 17:320–31. doi: 10.1016/j.omto.2020.04.003

  • 93

    DiFeoANarlaGCamacho-VanegasONishioHRoseSLBullerREet al. E-cadherin is a novel transcriptional target of the KLF6 tumor suppressor. Oncogene (2006) 25:6026–31. doi: 10.1038/sj.onc.1209611

  • 94

    DuCChenLZhangHWangZLiuWXieXet al. Caveolin-1 limits the contribution of BKCa channel to MCF-7 breast cancer cell proliferation and invasion. Int J Mol Sci (2014) 15:20706–22. doi: 10.3390/ijms151120706

  • 95

    GrantFMYangJNasrallahRClarkeJSadiyahFWhitesideSKet al. BACH2 drives quiescence and maintenance of resting Treg cells to promote homeostasis and cancer immunosuppression. J Exp Med (2020) 217(9):e20190711. doi: 10.1084/jem.20190711

  • 96

    KangJParkJHLeeHJJoUParkJKSeoJHet al. Caveolin-1 Modulates Docetaxel-Induced Cell Death in Breast Cancer Cell Subtypes through Different Mechanisms. Cancer Res Treat (2016) 48:715–26. doi: 10.4143/crt.2015.227

  • 97

    KimSKJangHRKimJHNohSMSongKSKimMRet al. The epigenetic silencing of LIMS2 in gastric cancer and its inhibitory effect on cell migration. Biochem Biophys Res Commun (2006) 349:1032–40. doi: 10.1016/j.bbrc.2006.08.128

  • 98

    LiuGHouGLiLLiYZhouWLiuL. Potential diagnostic and prognostic marker dimethylglycine dehydrogenase (DMGDH) suppresses hepatocellular carcinoma metastasis in vitro and in vivo. Oncotarget (2016) 7:32607–16. doi: 10.18632/oncotarget.8927

  • 99

    SlavinDAKoritschonerNPPrietoCCLopez-DiazFJChattonBBoccoJL. A new role for the Kruppel-like transcription factor KLF6 as an inhibitor of c-Jun proto-oncoprotein function. Oncogene (2004) 23:8196–205. doi: 10.1038/sj.onc.1208020

  • 100

    DesaiSDReedREBurksJWoodLMPullikuthAKHaasALet al. ISG15 disrupts cytoskeletal architecture and promotes motility in human breast cancer cells. Exp Biol Med (Maywood) (2012) 237:3849. doi: 10.1258/ebm.2011.011236

  • 101

    GillMKChristovaTZhangYYGregorieffAZhangLNarimatsuMet al. A feed forward loop enforces YAP/TAZ signaling during tumorigenesis. Nat Commun (2018) 9:3510. doi: 10.1038/s41467-018-05939-2

  • 102

    SongWHwangYYoungbloodVMCookRSBalkoJMChenJet al. Targeting EphA2 impairs cell cycle progression and growth of basal-like/triple-negative breast cancers. Oncogene (2017) 36:5620–30. doi: 10.1038/onc.2017.170

  • 103

    ZhouHTelonisAGJingYXiaNLBiedermanLJimboMet al. GPRC5A is a potential oncogene in pancreatic ductal adenocarcinoma cells that is upregulated by gemcitabine with help from HuR. Cell Death Dis (2016) 7:e2294. doi: 10.1038/cddis.2016.169

  • 104

    GonzalesMLAdamsSDunawayKWLasalleJM. Phosphorylation of distinct sites in MeCP2 modifies cofactor associations and the dynamics of transcriptional regulation. Mol Cell Biol (2012) 32:2894–903. doi: 10.1128/MCB.06728-11

  • 105

    StefanelliGGandagliaACostaMCheemaMSDiMDBarbieroIet al. Brain phosphorylation of MeCP2 at serine 164 is developmentally regulated and globally alters its chromatin association. Sci Rep (2016) 6:28295. doi: 10.1038/srep28295

  • 106

    TaoJHuKChangQWuHShermanNEMartinowichKet al. Phosphorylation of MeCP2 at Serine 80 regulates its chromatin association and neurological function. Proc Natl Acad Sci USA (2009) 106:4882–7. doi: 10.1073/pnas.0811648106

  • 107

    ZocchiLSassone-CorsiP. SIRT1-mediated deacetylation of MeCP2 contributes to BDNF expression. Epigenetics (2012) 7:695700. doi: 10.4161/epi.20733

  • 108

    BaertschMAHillengassJBlockaJSchonlandSHegenbartUGoldschmidtHet al. Efficacy and tolerability of the histone deacetylase inhibitor panobinostat in clinical practice. Hematol Oncol (2018) 36:210–6. doi: 10.1002/hon.2462

  • 109

    CostelloJFFruhwaldMCSmiragliaDJRushLJRobertsonGPGaoXet al. Aberrant CpG-island methylation has non-random and tumour-type-specific patterns. Nat Genet (2000) 24:132–8. doi: 10.1038/72785

  • 110

    JonesPAIssaJPBaylinS. Targeting the cancer epigenome for therapy. Nat Rev Genet (2016) 17:630–41. doi: 10.1038/nrg.2016.93

  • 111

    Simo-RiudalbasLEstellerM. Targeting the histone orthography of cancer: drugs for writers, erasers and readers. Br J Pharmacol (2015) 172:2716–32. doi: 10.1111/bph.12844

  • 112

    Gonzalez-GiraldoYForeroDAEcheverriaVGarcia-SeguraLMBarretoGE. Tibolone attenuates inflammatory response by palmitic acid and preserves mitochondrial membrane potential in astrocytic cells through estrogen receptor beta. Mol Cell Endocrinol (2019) 486:6578. doi: 10.1016/j.mce.2019.02.017

  • 113

    MahmoodNRabbaniSA. DNA Methylation Readers and Cancer: Mechanistic and Therapeutic Applications. Front Oncol (2019) 9:489. doi: 10.3389/fonc.2019.00489

  • 114

    MinucciSPelicciPG. Histone deacetylase inhibitors and the promise of epigenetic (and more) treatments for cancer. Nat Rev Cancer (2006) 6:3851. doi: 10.1038/nrc1779

  • 115

    SinghAKBishayeeAPandeyAK. Targeting Histone Deacetylases with Natural and Synthetic Agents: An Emerging Anticancer Strategy. Nutrients (2018) 10(6):731. doi: 10.3390/nu10060731

  • 116

    TemianDCPopLAIrimieAIBerindan-NeagoeI. The Epigenetics of Triple-Negative and Basal-Like Breast Cancer: Current Knowledge. J Breast Cancer (2018) 21:233–43. doi: 10.4048/jbc.2018.21.e41

  • 117

    GelatoKAShaikhibrahimZOckerMHaendlerB. Targeting epigenetic regulators for cancer therapy: modulation of bromodomain proteins, methyltransferases, demethylases, and microRNAs. Expert Opin Ther Targets (2016) 20:783–99. doi: 10.1517/14728222.2016.1134490

  • 118

    LengauerCKinzlerKWVogelsteinB. DNA methylation and genetic instability in colorectal cancer cells. Proc Natl Acad Sci USA (1997) 94:2545–50. doi: 10.1073/pnas.94.6.2545

  • 119

    MomparlerRL. Cancer epigenetics. Oncogene (2003) 22:6479–83. doi: 10.1038/sj.onc.1206774

  • 120

    MagdinierFBillardLMWittmannGFrappartLBenchaibMLenoirGMet al. Regional methylation of the 5’ end CpG island of BRCA1 is associated with reduced gene expression in human somatic cells. FASEB J (2000) 14:1585–94. doi: 10.1096/fj.99-0817com

  • 121

    ListerRPelizzolaMDowenRHHawkinsRDHonGTonti-FilippiniJet al. Human DNA methylomes at base resolution show widespread epigenomic differences. Nature (2009) 462:315–22. doi: 10.1038/nature08514

  • 122

    BallasNGrunseichCLuDDSpehJCMandelG. REST and its corepressors mediate plasticity of neuronal gene chromatin throughout neurogenesis. Cell (2005) 121:645–57. doi: 10.1016/j.cell.2005.03.013

  • 123

    KokuraKKaulSCWadhwaRNomuraTKhanMMShinagawaTet al. The Ski protein family is required for MeCP2-mediated transcriptional repression. J Biol Chem (2001) 276:34115–21. doi: 10.1074/jbc.M105747200

  • 124

    FilarskyMZillnerKArayaIVillar-GareaAMerklRLangstGet al. The extended AT-hook is a novel RNA binding motif. RNA Biol (2015) 12:864–76. doi: 10.1080/15476286.2015.1060394

  • 125

    MaxwellSSPelkaGJTamPPEl-OstaA. Chromatin context and ncRNA highlight targets of MeCP2 in brain. RNA Biol (2013) 10:1741–57. doi: 10.4161/rna.26921

  • 126

    NikitinaTShiXGhoshRPHorowitz-SchererRAHansenJCWoodcockCL. Multiple modes of interaction between the methylated DNA binding protein MeCP2 and chromatin. Mol Cell Biol (2007) 27:864–77. doi: 10.1128/MCB.01593-06

  • 127

    RatnakumarKBernsteinE. ATRX: the case of a peculiar chromatin remodeler. Epigenetics (2013) 8:39. doi: 10.4161/epi.23271

  • 128

    CohenSGabelHWHembergMHutchinsonANSadaccaLAEbertDHet al. Genome-wide activity-dependent MeCP2 phosphorylation regulates nervous system development and function. Neuron (2011) 72:7285. doi: 10.1016/j.neuron.2011.08.022

  • 129

    SharmaMCastro-PiedrasISimmonsGEJr.PruittK. Dishevelled: A masterful conductor of complex Wnt signals. Cell Signal (2018) 47:5264. doi: 10.1016/j.cellsig.2018.03.004

  • 130

    SharmaMMolehinDCastro-PiedrasIMartinezEGPruittK. Acetylation of conserved DVL-1 lysines regulates its nuclear translocation and binding to gene promoters in triple-negative breast cancer. Sci Rep (2019) 9:16257. doi: 10.1038/s41598-019-52723-3

  • 131

    MannaPRAhmedAUVartakDMolehinDPruittK. Overexpression of the steroidogenic acute regulatory protein in breast cancer: Regulation by histone deacetylase inhibition. Biochem Biophys Res Commun (2019) 509:476–82. doi: 10.1016/j.bbrc.2018.12.145

  • 132

    MolehinDCastro-PiedrasISharmaMSennouneSRArenaDMannaPRet al. Aromatase Acetylation Patterns and Altered Activity in Response to Sirtuin Inhibition. Mol Cancer Res (2018) 16:1530–42. doi: 10.1158/1541-7786.MCR-18-0047

  • 133

    HollowayKRBarbieriAMalyarchukSSaxenaMNedeljkovic-KurepaACameronMMet al. SIRT1 positively regulates breast cancer associated human aromatase (CYP19A1) expression. Mol Endocrinol (2013) 27:480–90. doi: 10.1210/me.2012-1347

  • 134

    SimmonsGEJrPandeySNedeljkovic-KurepaASaxenaMWangAPruittK. Frizzled 7 expression is positively regulated by SIRT1 and beta-catenin in breast cancer cells. PLoS ONE (2014) 9:e98861. doi: 10.1371/journal.pone.0098861

  • 135

    PruittKZinnRLOhmJEMcGarveyKMKangSHWatkinsDNet al. Inhibition of SIRT1 reactivates silenced cancer genes without loss of promoter DNA hypermethylation. PloS Genet (2006) 2:e40. doi: 10.1371/journal.pgen.0020040

  • 136

    Castro-PiedrasISharmaMdenBMMolehinDMartinezEGVartakDet al. DVL1 and DVL3 differentially localize to CYP19A1 promoters and regulate aromatase mRNA in breast cancer cells. Oncotarget (2018) 9:35639–54. doi: 10.18632/oncotarget.26257

  • 137

    SimmonsGEJrPruittWMPruittK. Diverse roles of SIRT1 in cancer biology and lipid metabolism. Int J Mol Sci (2015) 16:950–65. doi: 10.3390/ijms16010950

Summary

Keywords

MeCP2, breast cancer, post-translational modification, ChIP-Seq analysis, transcriptional regulation

Citation

Castro-Piedras I, Vartak D, Sharma M, Pandey S, Casas L, Molehin D, Rasha F, Fokar M, Nichols J, Almodovar S, Rahman RL and Pruitt K (2020) Identification of Novel MeCP2 Cancer-Associated Target Genes and Post-Translational Modifications. Front. Oncol. 10:576362. doi: 10.3389/fonc.2020.576362

Received

25 June 2020

Accepted

26 October 2020

Published

10 December 2020

Volume

10 - 2020

Edited by

Monther Al-alwan, King Faisal Specialist Hospital & Research Centre, Saudi Arabia

Reviewed by

Barry Zee, Boston Children's Hospital and Harvard Medical School, United States; Ola Hermanson, Karolinska Institutet (KI), Sweden

Updates

Copyright

*Correspondence: Kevin Pruitt,

This article was submitted to Molecular and Cellular Oncology, a section of the journal Frontiers in Oncology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics