Next Article in Journal
FGF23 and Fetuin-A Interaction and Mesenchymal Osteogenic Transformation
Previous Article in Journal
Copy-Paste Mutagenesis: A Method for Large-Scale Alteration of Viral Genomes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Effect of Single CpG Demethylation on the Pattern of DNA-Protein Binding

Laboratory of Calcium Binding Proteins, Department of Molecular and Cellular Neurobiology, Nencki Institute of Experimental Biology of the Polish Academy of Sciences, 3 Pasteur Street, 02-093 Warsaw, Poland
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(4), 914; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20040914
Submission received: 25 January 2019 / Revised: 15 February 2019 / Accepted: 18 February 2019 / Published: 20 February 2019
(This article belongs to the Section Molecular Genetics and Genomics)

Abstract

:
Epidermal differentiation is a complex process and its regulation may involve epigenetic factors. Analysis of DNA methylation in 20 selected regions within the epidermal differentiation complex (EDC) gene cluster by targeted next-generation sequencing (NGS) detected no or only minor changes in methylation, mostly slight demethylation, occurring during the course of keratinocyte differentiation. However, a single CpG pair within the exon of the PGLYRP3 gene underwent a pronounced demethylation concomitant with an increase in PGLYRP3 expression. We have employed a DNA-affinity precipitation assay (DAPA) and mass spectrometry to examine changes in the composition of proteins that bind to DNA containing either methylated or unmethylated CpG. We found that the unmethylated probe attracted mostly RNA binding proteins, including splicing factors, which suggests that demethylation of this particular CpG may facilitate PGLYRP3 transcription and/or pre-mRNA splicing.

1. Introduction

Keratinocyte differentiation is a highly coordinated process that ensures the proper function of the epidermis. In the course of differentiation, keratinocytes produce numerous specific proteins that serve to build a physical barrier against environmental hazards, infections and organismal water loss [1]. The epidermal differentiation complex (EDC) is a 1.7 Mb gene cluster on human chromosome 1q21 that houses many epidermis-specific genes [2]. Namely, EDC comprises genes of four protein families: the S100 proteins, the S100-fused type proteins (SFTP), the small proline-rich region (SPRR) proteins, late cornified envelope (LCE) proteins, and several other genes [3,4,5]. The EDC genes have a different onset of expression during fetal and postnatal development [6]. Furthermore, their expression changes during keratinocyte differentiation. The regulatory mechanisms that control EDC gene expression have not been elucidated yet.
DNA methylation is an epigenetic mechanism of gene expression regulation predominantly associated with gene silencing, as methylated cytosine residues attract components of chromatin repressive complexes [7]. DNA methylation pattern is acquired during early embryonic development and is characteristic for a given cell type [8]. Marked changes in the relatively stable methylation landscape are usually observed in cancer cells and each tumor type is associated with a unique set of methylated genes [9]. However, increasing evidence indicates that changes in cytosine methylation may occur in response to various environmental stimuli or diet [10], during cell or organismal aging [11], and during cell differentiation [12]. Involvement of DNA methylation in regulation of EDC gene expression upon keratinocyte differentiation may be considered based on studies showing differential expression of EDC genes upon treatment with DNMT inhibitor 5-azacytidine [13] or following DNMT1 knock-out [14]. Furthermore, changes in cytosine methylation in EDC genes occur in various skin pathologies, which are associated with disturbed epidermal differentiation [15]. We have recently compared the extent of DNA methylation in promoters and other regions of selected EDC genes in undifferentiated and differentiated primary human keratinocytes [16]. This analysis did not detect major changes in cytosine methylation that would correspond to the robust changes in expression of the examined genes. In the current analysis, using targeted next-generation sequencing (NGS), we have examined DNA methylation in 20 selected EDC areas located within DNase I hypersensitive sites and H3K27ac/H3K4me1 abundant sites, that is in presumptive regulatory regions, as a function of keratinocyte differentiation. We have detected marked demethylation of a single cytosine residue located in the 5th exon of the PGLYRP3 gene, which associated with an increase in PGLYRP3 mRNA level. Using a DNA-affinity precipitation assay (DAPA) and mass spectrometry we have investigated differences in the composition of proteins that bound to probes containing unmethylated or methylated CpG. This analysis revealed a marked switch in the set of bound proteins. In particular, significantly more RNA-binding proteins, including splicing factors, appeared to bind to the unmethylated probe suggesting that demethylation of a single cytosine residue may promote PGLYRP3 transcription and/or pre-mRNA processing.

2. Results

2.1. DNA Methylation Analysis

In order to search for possible changes in DNA methylation that might occur during keratinocyte differentiation, we have selected 20 sequences mostly located within putative enhancer regions in EDC of normal human epidermal keratinocytes (NHEK), defined based on the presence of H3K27ac/H3K4m1 enriched and DNase I hypersensitive sites (USCS Genome Browser, GRCh37/hg19). Additional criteria for sequence selection were the preferable location in intergenic regions separating gene clusters or individual genes, number of CpGs, the density of transcription factor (TF) binding sites, and the presence of binding sites for methylation-sensitive TFs [17] and/or for TFs that bind to keratinocyte-specific enhancers [18]. The analyzed sequences covered 4781 bp and 103 CpG dinucleotides. The size and genomic location of these sequences together with the number of reads obtained in targeted NGS is presented in Table S1. The first analysis, performed on DNA from NHEK derived from an individual donor (Figure 1A), covered 17 sequences. It showed that, in general, the analyzed regions tend to be either highly methylated (75–100% methylation of individual CpGs) or unmethylated (0–10% methylation) with only one sequence (no. 99) showing a mixed methylation pattern. Upon differentiation, the level of methylation of individual CpGs remained unchanged or changed only slightly with most cytosine residues undergoing minor demethylation (by about 5%). There were only several cytosine residues that underwent marked demethylation. To verify these changes in cytosine methylation, the highly methylated regions were analyzed again in DNA derived from a heterogeneous keratinocyte population (NHEK obtained from many donors). The results (Figure 1B) confirmed that the overall methylation pattern of the studied regions is the same as in the material derived from an individual donor. The only cytosine residue that underwent marked demethylation according to both analyses was contained in sequence 95 and corresponded to the 5th exon of the PGLYRP3 gene (Figure S1). To check whether this demethylation event coincided with altered gene expression, we have performed PCR analysis, which showed a concomitant increase in PGLYRP3 mRNA level in differentiated versus undifferentiated NHEK and HaCaT cells (Figure 2).

2.2. Analysis of Protein Binding to a Sequence Containing the Variably Methylated CpG Pair

To explore whether the prominent change in methylation of one of the examined CpGs, which occurred during differentiation of NHEK, induces a switch in the amount/composition of protein-DNA complexes formed in its vicinity, we have performed EMSA assays using a 26 bp long oligonucleotide containing the examined CpG pair in either unmethylated or methylated form. As shown in Figure 3, prominent protein complexes of identical mobility bound to both probes. The same mobility of the complexes and the fact that the binding of each complex could be competed by the excess of both the unmethylated and methylated probes suggested that their binding, although specific, was not critically dependent on the CpG methylation state. Therefore, to look for differently bound proteins, we have performed a DNA-affinity precipitation assay (DAPA) followed by mass spectrometry analysis. To focus our analysis on proteins whose binding was dependent on the presence of the CpG dinucleotide in the examined sequence, we have preincubated the nuclear extract from HaCaT cells with a mutated sequence, containing TA instead of CG. The unbound fraction was then incubated with either the unmethylated or methylated sequence. Altogether, four experiments were performed. Among all the bound proteins, we searched for those that bound uniquely to either of the two sequences. Depending on the experiment, we identified 44–149 proteins that bound only to the unmethylated sequence and 7–144 that bound to the methylated probe. As could be expected, proteins classified by the PANTHER program as nucleic acid binding proteins and transcription factors constituted the largest class of proteins bound to each sequence. Namely, those proteins accounted for 54.27% of all proteins bound to the unmethylated sequence and for 38.57% of proteins bound to the methylated one (Figure 4). As to the unmethylated probe, 10–37% of these proteins, depending on the experiment, were identified in at least two analyses (and were not identified as bound to the methylated probe in any of the analyses). For the methylated probe, the set of bound proteins was more heterogeneous across experiments and the respective percentage amounted to 11–28%. Most interestingly, when nucleic acid binding proteins were further classified into DNA- and RNA-binding proteins, the proportion was definitely different for the two sequences (Figure 4). RNA-binding proteins constituted, on average, 87.25% of proteins bound to the unmethylated sequence, that is significantly more than in the case of the methylated one (57.15%). Of those, splicing factors of the SRSF family, RPLs or SMARC proteins were found, rather consistently, among proteins bound to the unmethylated probe (Table S2). On the other hand, there were, on average, more proteins belonging to the DNA-binding protein category bound to the methylated probe (Figure 4). Of note, MBD2, a protein that binds methylated cytosine, was identified in two analyses.

3. Discussion

The targeted NGS analysis of 20 selected sequences within EDC localized in presumptive regulatory/enhancer regions revealed that, in these particular regions, DNA tends to be either highly (75–100%) or hardly (0–5%) methylated, concordant with the general pattern of genome methylation. Of note, as indicated in the Methods Section, bisulfite-based techniques cannot distinguish between methylated and hydroxymethylated cytosine so the presence of the latter modification in the examined regions cannot be excluded. Upon differentiation of primary human keratinocytes, the majority of the analyzed CpG residues were subject to only minor changes in the extent of methylation—with slight demethylation (1–5%) representing the prevailing trend. As an exception to this rule, we noticed a prominent decrease in methylation of a single CpG dinucleotide in primary keratinocytes derived both from an individual donor and in a keratinocyte population derived from many donors. The consistency and uniqueness of this demethylation event prompted us to examine it further.
The CpG pair undergoing demethylation is located in the 5th exon of the PGLYRP3 gene (Figure S1). The peptidoglycan recognition protein 3 (PGLYRP3) encoded by this gene is an innate immunity molecule, one of the four mammalian representatives of the PGRP protein family [19]. Two PGLYRP genes, encoding PGLYRP3 and PGLYRP4, are located in EDC. According to an interesting hypothesis, the ancestral PGRYLP and S100 genes, through multiple fusion and duplication events, could have given rise to all EDC gene families [20]. PGLYRP3 and PGLYRP4 are expressed mainly in the epidermis and epidermal appendages such as hair follicles, sebaceous glands and sweat glands [19]. Both proteins can be secreted and are bactericidal for many pathogenic and nonpathogenic bacteria. Changes in PGLYRP3 expression have been described in the context of microbial components and Toll-like receptor activation [21] or in skin pathologies [22], but less information is available concerning PGLYRP3 expression during epidermal differentiation. Our results indicate that the level of PGLYRP3 mRNA increased during differentiation of both primary human keratinocytes (NHEK) and keratinocyte-derived HaCaT cells. This increase occurs concomitantly with the loss of methylation of the CpG dinucleotide in the 5th exon of this gene, but a causative link between these two events could not be proved at the moment.
There are numerous reports that methylation of a single CpG pair, either in gene promoters [23] or introns and exons [24] can silence gene expression. Conversely, demethylation at single CpG residues should be able to trigger gene expression and such events have been observed in gene promoters [25,26]. It is not known how methylation or demethylation of a single CpG residue, especially within an otherwise unmethylated or methylated genomic area, is able to influence the structure of chromatin and DNA-protein interactions to execute a profound change in transcription rate. To get some insight into this issue, we have performed EMSA using double-stranded oligonucleotides, which matched the sequence containing the examined CpG pair, in either methylated or unmethylated form. The analysis detected a specific protein complex, which, however, bound to the examined sequence regardless of its methylation and could not, therefore, account for an altered expression rate. To further pursue this issue we employed DNA-affinity precipitation assay (DAPA) followed by mass spectrometry analysis to detect proteins that bound differentially to the methylated or unmethylated probes. This assay, followed by Western blot analysis, is routinely used to examine the binding of transcription factors and other proteins of interest to the analyzed sequence [27]. The combination of DAPA and mass spectrometry has been used earlier, for example, to compare protein binding to hydroxymethylcytosine and formylcytosine [28]. Using this experimental setup, we were able to discern a set of proteins that bound preferably to each probe. Interestingly, the unmethylated probe attracted mostly RNA binding proteins, including splicing factors, which were poorly represented or absent from the pool of proteins bound to the methylated probe. Given that the cytosine residue undergoing prominent demethylation during keratinocyte differentiation is located within an exon, this may indicate that the unmethylated sequence is more accessible to the transcriptional and/or spliceosomal machinery than the methylated one. Thus, loss of methylation could result in increased PGLYRP3 expression upon differentiation. DNA methylation in exons is implicated in both transcription and splicing. Methylation within exons is generally believed to favor transcription [12]. On the other hand, the exact role of cytosine methylation in the regulation of alternative splicing is still a subject of intensive studies (Maor et al., 2015) [29]. Our observation that splicing factors of the SRSF and U2AF family bind preferably to unmethylated DNA sequence may be a valuable contribution to solving this issue.

4. Methods

4.1. Cell Culture

Normal human epidermal keratinocytes (NHEK) were purchased from PromoCell and used at early passages. Cells were maintained in Keratinocyte Growth Medium 2 (PromoCell, Heidelberg, Germany) supplemented with 0.06 mM calcium chloride, human keratinocyte growth supplement (HKGS), 100 U/mL penicillin and 100 µg/mL streptomycin. Cell differentiation was induced by increasing calcium chloride concentration to 1.8 mM. HaCaT cells, spontaneously immortalized human keratinocytes, purchased from Cell line Service (Germany), were cultured and differentiated as described in Reference [30]. Cells were cultured in a cell incubator at 37 °C and 5% CO2.

4.2. DNA Methylation Analysis

DNA was isolated from undifferentiated and differentiated NHEK using a standard method and conditions as described in Sobiak et al., 2016. The DNA was then subjected to bisulfite treatment, which converts cytosine into uracil while methylated (or hydroxymethylated) cytosine remains intact. Twenty selected regions within EDC, 119–314 bp in length, were amplified and subjected to targeted next-generation sequencing (NGS) reaction performed by EpigenDX (Hopkinton, MA, USA). The minimal sequencing depth was set at 30. The analysis was performed twice: once on DNA from NHEK derived from an individual donor and the second time on DNA from cell population obtained from many donors.

4.3. Preparation of Nuclear Extracts

HaCaT cells grown in a 10 cm plate were scraped in PBS and pelleted by centrifugation at 600× g. Cells were then resuspended in 0.5 mL of hypotonic buffer (20 mM Hepes pH 8.0, 0.2% Igepal CA-630, 1 mM EDTA, 1 mM DTT) supplemented with protease inhibitor cocktail (Roche) and incubated for 10 min on ice. The lysate was centrifuged at 5000× g for 10 min and the resulting pellet containing cell nuclei was incubated for 30 min on ice in the extraction buffer (20 mM Hepes pH 8.0, 450 mM NaCl, 20% glycerol, 1 mM EDTA, 1 mM DTT and protease inhibitor cocktail) with occasional mixing. Following centrifugation at 12,000× g for 15 min at 4 °C the supernatant was collected and, after protein determination with protein assay dye reagent (Bio-Rad, Hercules, CA, USA), used immediately or portioned and kept at −20 °C.

4.4. RNA Isolation and PCR

Total RNA was isolated from undifferentiated and differentiated NHEK and HaCaT cells using RNeasy mini kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. Two microgram of RNA was used for cDNA synthesis catalyzed by M-MLV reverse transcriptase (Sigma Aldrich, St. Louis, MO, USA) in the presence of random nonamers (Sigma Aldrich). To assess the level of PGLYRP3 transcript, a PCR reaction was conducted with the following primers: 5′-CGCCCAGGTATATTCAGCCA-3′ and 5′-CCGTGTGTCCATGTGAAAGG-3′ for 35 cycles.

4.5. Electrophoretic Mobility Shift Assay (EMSA)

A 26 base long oligonucleotide, the sequence of which corresponded to the fragment of the probe no. 95 containing the differentially methylated cytosine residue, was synthesized at the DNA Sequencing and Oligonucleotides Synthesis Laboratory, Institute of Biochemistry and Biophysics, PAS, Warsaw, Poland, in a form containing either unmethylated (5′-TCCCAAGCAGATCGTTTGATGATGTT-3’) or methylated cytosine (5′-TCCCAAGCAGATmCGTTTGATGATGTT-3’). Both oligonucleotides were hybridized with complementary strands and radioactively labelled using T4 kinase and (γ32P) ATP. An amount of probe equivalent to 100,000 cpm was incubated with 10 µg of nuclear extract from HaCaT cells (grown in standard DMEM). The labeling and binding reactions and the subsequent gel electrophoresis were performed essentially as described earlier [31].

4.6. DNA Affinity Precipitation Assay (DAPA)

The two oligonucleotides used in EMSA (5′-TCCCAAGCAGATCGTTTGATGATGTT-3’ and 5′ -TCCCAAGCAGATmCGTTTGATGATGTT-3’) and a corresponding oligonucleotide with the CG pair mutated (5′-TCCCAAGCAGATATTTTGATGATGTT-3’) were synthesized and biotinylated at the 5′-termini by the same facility as oligonucleotides used for EMSA. They were then annealed with their complementary, non-biotinylated, strands. The DNA affinity precipitation assay was performed as follows: 600 μg of nuclear extract diluted 5 times in binding buffer (2 µg herring sperm DNA, 20 mM HEPES, pH 7.9, 2 mM MgCl2, 0.2 mM EDTA, 1 mM DTT, and 10% (v/v) glycerol) was first precleared by incubation with 10 μL of streptavidin-agarose beads for 1 h. After separation of the beads, the supernatant was first incubated with 2 μg of the mutated DNA probe and 40 μL of streptavidin-agarose for 1 h at 4 °C. The unbound fraction was divided into two equal portions, which were incubated with either the methylated or unmethylated probe as above. Beads were collected and washed three times with the binding buffer containing 0.5% Nonidet P-40. Proteins bound to the beads were eluted with elution buffer (5 mM Tris-HCl, pH 7.5, 0.5 mM EDTA, 0.5 M NaCl), and subjected to mass spectrometry analysis performed at the Mass Spectrometry Laboratory, Institute of Biochemistry and Biophysics, PAS, Warsaw, on an Orbitrap spectrometer. Data were analyzed by the Mascot program (www.matrixscience.com) based on the NCBI sequence database. The cut-off point was set at Mascot score 50.

4.7. Statistical Analysis

Analysis of proteins retrieved by Mascot was performed using PANTHER (protein analysis through evolutionary relationships) (http://www.pantherdb.org) and Fisher’s exact test with FDR correction. Statistical data analysis of data presented in Figure 2 and Figure 4 was performed using Student’s t-test.

5. Conclusions

Demethylation of a single CpG in the 5th exon of PGLYRP3 facilitates binding of proteins involved in RNA processing, including splicing, and may contribute to the enhanced expression of this protein during keratinocyte differentiation.

Supplementary Materials

Supplementary materials can be found at https://0-www-mdpi-com.brum.beds.ac.uk/1422-0067/20/4/914/s1.

Author Contributions

Conceptualization, W.L. and B.S.; methodology, W.L. and B.S.; validation, W.L. and B.S.; formal analysis, W.L.; investigation B.S.; writing—original draft preparation, W.L.; writing—review and editing, W.L. and B.S.; visualization, B.S.; supervision, W.L.; project administration, W.L. and B.S.; funding acquisition, W.L. and B.S.

Funding

This work was supported by a National Science Center grants: 2016/23/N/NZ2/02376 to B.S. and 2012/05/B/NZ2/02944 to W.L. and by the statutory funds from the Nencki Institute of Experimental Biology, PAS.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fuchs, E.; Raghavan, S. Getting under the skin of epidermal morphogenesis. Nat. Rev. Genet. 2002, 3, 199–209. [Google Scholar] [CrossRef] [PubMed]
  2. Mischke, D.; Korge, B.P.; Marenholz, I.; Volz, A.; Ziegler, A. Genes encoding structural proteins of epidermal cornification and S100 calcium-binding proteins form a gene complex “epidermal differentiation complex” on human chromosome 1q21. J. Investig. Dermatol. 1996, 106, 989–992. [Google Scholar] [CrossRef] [PubMed]
  3. Cabral, A.; Voskamp, P.; Cleton-Jansen, A.M.; South, A.; Nizetic, D.; Backendorf, C. Structural organization and regulation of the small proline-rich family of cornified envelope precursors suggest a role in adaptive barrier function. J. Biol. Chem. 2001, 276, 19231–19237. [Google Scholar] [CrossRef] [PubMed]
  4. Kypriotou, M.; Huber, M.; Hohl, D. The human epidermal differentiation complex: Cornified envelope precursors, S100 proteins and the “fused genes” family. Exp. Dermatol. 2012, 21, 643–649. [Google Scholar] [CrossRef] [PubMed]
  5. Leśniak, W.; Graczyk-Jarzynka, A. The S100 proteins in epidermis: Topology and function. Biochim. Biophys. Acta 2015, 1850, 2563–2572. [Google Scholar] [CrossRef] [PubMed]
  6. Bickenbach, J.R.; Greer, J.M.; Bundman, D.S.; Rothnagel, J.A.; Roop, D.R. Loricrin expression is coordinated with other epidermal proteins and the appearance of lipid lamellar granules in development. J. Investig. Dermatol. 1995, 104, 405–410. [Google Scholar] [CrossRef] [PubMed]
  7. Bird, A.P.; Wolffe, A.P. Methylation-induced repression-belts, braces and chromatin. Cell 1999, 99, 451–454. [Google Scholar] [CrossRef]
  8. Okano, M.; Bell, D.W.; Haber, D.A.; Li, E. DNA methyltransferases Dnmt3a and Dnmt3b are essential for de novo methylation and mammalian development. Cell 1999, 99, 247–257. [Google Scholar] [CrossRef]
  9. Das, P.M.; Singal, R. DNA Methylation and Cancer. J. Clin. Oncol. 2004, 22, 4632–4642. [Google Scholar] [CrossRef]
  10. Crider, K.S.; Yang, T.P.; Berry, R.J.; Bailey, L.B. Folate and DNA methylation: A review of molecular mechanisms and the evidence for folate’s role. Adv. Nutr. 2012, 3, 21–38. [Google Scholar] [CrossRef]
  11. Johnson, A.A.; Akman, K.; Calimport, S.R.G.; Wuttke, D.; Stolzing, A.; de Magalhaes, J.P. The role of DNA methylation in aging, rejuvenation, and age-related disease. Rejuv. Res. 2012, 15, 483–494. [Google Scholar] [CrossRef] [PubMed]
  12. Laurent, L.; Wong, E.; Li, G.; Huynh, T.; Tsirigos, A.; Ong, C.T.; Low, H.M.; Kin Sung, K.W.; Rigoutsos, I.; Loring, J.; et al. Dynamic changes in the human methylome during differentiation. Genome Res. 2010, 20, 320–331. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Elder, J.T.; Zhao, X. Evidence for local control of gene expression in the epidermal differentiation complex. Exp. Dermatol. 2002, 11, 406–412. [Google Scholar] [CrossRef] [PubMed]
  14. Sen, G.L.; Reuter, J.A.; Webster, D.E.; Zhu, L.; Khavari, P.A. DNMT1 maintains progenitor function in self-renewing somatic tissue. Nature 2010, 463, 563–567. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Roberson, E.D.O.; Liu, Y.; Ryan, C.; Joyce, C.E.; Duan, S.; Cao, L.; Martin, A.; Liao, W.; Menter, A.; Bowcock, A.M. A subset of methylated CpG sites differentiate psoriatic from normal skin. J. Investig. Dermatol. 2012, 132, 583–592. [Google Scholar] [CrossRef] [PubMed]
  16. Sobiak, B.; Graczyk-Jarzynka, A.; Leśniak, W. Comparison of DNA methylation and expression pattern of S100 and other Epidermal Differentiation Complex genes in differentiating keratinocytes. J. Cell. Biochem. 2016, 11, 1092–1098. [Google Scholar] [CrossRef] [PubMed]
  17. Zhu, H.; Wang, G.; Qian, J. Transcription factors as readers and effectors of DNA methylation. Nat. Rev. Genet. 2016, 17, 551–565. [Google Scholar] [CrossRef] [Green Version]
  18. Cavazza, A.; Miccio, A.; Romano, O.; Petiti, L.; Malagoli Tagliazucchi, G.; Peano, C.; Severgnini, M.; Rizzi, E.; De Bellis, G.; Bicciato, S.; et al. Dynamic transcriptional and epigenetic regulation of human epidermal keratinocyte differentiation. Stem Cell Rep. 2016, 6, 618–632. [Google Scholar] [CrossRef]
  19. Dziarski, R.; Gupta, D. The peptidoglycan recognition proteins (PGRPs). Genome Biol. 2006, 7, 232. [Google Scholar] [CrossRef]
  20. Strasser, B.; Mlitz, V.; Hermann, M.; Rice, R.H.; Eigenheer, R.A.; Alibardi, L.; Tschachler, E.; Eckhart, L. Evolutionary origin and diversification of epidermal barrier proteins in Amniotes. Mol. Biol. Evol. 2014, 31, 3194–3205. [Google Scholar] [CrossRef]
  21. Ma, P.; Wang, Z.; Pflugfelder, S.C.; Li, D.Q. Toll-like receptors mediate induction of peptidoglycan recognition proteins in human corneal epithelial cells. Exp. Eye Res. 2010, 90, 130–136. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Park, S.Y.; Gupta, D.; Kim, C.H.; Dziarski, R. Differential effects of peptidoglycan recognition proteins on experimental atopic and contact dermatitis mediated by Treg and Th17 cells. PLoS ONE 2011, 6, e24961. [Google Scholar] [CrossRef] [PubMed]
  23. Wang, T.; Chen, M.; Liu, L.; Cheng, H.; Yan, Y.E.; Feng, Y.H.; Wang, H. Nicotine induced CpG methylation of Pax6 binding motif in StAR promoter reduces the gene expression and cortisol production. Toxicol. Appl. Pharmacol. 2011, 257, 328–337. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Zhang, X.; Wu, M.; Xiao, H.; Lee, M.T.; Levin, L.; Leung, Y.K.; Ho, S.M. Methylation of a single intronic CpG mediates expression silencing of the PMP24 gene in prostate cancer. Prostate 2010, 70, 765–776. [Google Scholar] [CrossRef] [PubMed]
  25. Venza, I.; Visalli, M.; Fortunato, C.; Ruggeri, M.; Ratone, S.; Caffo, M.; Caruso, G.; Alafaci, C.; Tomasello, F.; Teti, D.; et al. PGE2 induces interleukin-8 derepression in human astrocytoma through coordinated DNA demethylation and histone hyperacetylation. Epigenetics 2012, 7, 1315–1330. [Google Scholar] [CrossRef]
  26. Ceccarelli, V.; Racanicchi, S.; Martelli, M.P.; Nocentini, G.; Fettucciari, K.; Riccardi, C.; Marconi, P.; Di Nardo, P.; Grignani, F.; Binaglia, L.; et al. Eicosapentaenoic acid demethylates a single CpG that mediates expression of tumor suppressor CCAAT/enhancer-binding protein delta in U937 leukemia cells. J. Biol. Chem. 2011, 286, 27092–27102. [Google Scholar] [CrossRef] [PubMed]
  27. Miller, D.E.; Patel, Z.H.; Lu, X.; Lynch, A.T.; Weirauch, M.T.; Kottyan, L.C. Screening for functional non-coding genetic variants using electrophoretic mobility shift assay (EMSA) and DNA-affinity precipitation assay (DAPA). J. Vis. Exp. 2016, 114, 54093. [Google Scholar] [CrossRef]
  28. Iurlaro, M.; Ficz, G.; Oxley, D.; Raiber, E.A.; Bachman, M.; Booth, M.J.; Andrews, S.; Balasubramanian, S.; Reik, W. A screen for hydroxymethylcytosine and formylcytosine binding proteins suggests functions in transcription and chromatin regulation. Genome Biol. 2013, 14, R119. [Google Scholar] [CrossRef]
  29. Maor, G.L.; Yearim, A.; Ast, G. The alternative role of DNA methylation in splicing regulation. Trends Genet. 2015, 31, 274–280. [Google Scholar] [CrossRef]
  30. Graczyk, A.; Lesniak, W. S100A6 expression in keratinocytes and its impact on epidermal differentiation. Int. J. Biochem. Cell Biol. 2014, 57, 135–141. [Google Scholar] [CrossRef]
  31. Leśniak, W.; Jezierska, A.; Kuznicki, J. Upstream stimulatory factor is involved in the regulation of the human calcyclin (S100A6) gene. Biochim. Biophys. Acta 2000, 1517, 73–81. [Google Scholar] [CrossRef]
Figure 1. Targeted next-generation sequencing NGS analysis of CpG methylation in selected regions of EDC in DNA from undifferentiated and differentiated primary keratinocytes (NHEK). (A) Analysis performed on DNA from cells derived from an individual donor. (B) Analysis of highly methylated sequences shown in (A) in DNA from cells obtained from many donors. Sequences 109 and 107 could be read only in analysis B, and sequences 98 and 100 in analysis A. Solid line—undifferentiated keratinocytes; dashed line—differentiated keratinocytes. Arrows indicate the demethylated cytosine residue in sequence no. 95.
Figure 1. Targeted next-generation sequencing NGS analysis of CpG methylation in selected regions of EDC in DNA from undifferentiated and differentiated primary keratinocytes (NHEK). (A) Analysis performed on DNA from cells derived from an individual donor. (B) Analysis of highly methylated sequences shown in (A) in DNA from cells obtained from many donors. Sequences 109 and 107 could be read only in analysis B, and sequences 98 and 100 in analysis A. Solid line—undifferentiated keratinocytes; dashed line—differentiated keratinocytes. Arrows indicate the demethylated cytosine residue in sequence no. 95.
Ijms 20 00914 g001
Figure 2. Analysis of PGLYRP3 mRNA level in undifferentiated (undiff) and differentiated (diff) primary keratinocytes (NHEK) and HaCaT cells. (A) representative PCR results. (B) statistical analysis of results from 3 experiments. PGLYRP3 mRNA level in undifferentiated NHEK or HaCaT cells is represented as 1.00 (white bar). The level in differentiated NHEK and HaCaT cells is represented by striped and gray bars, respectively. Data are presented as mean ± SEM.
Figure 2. Analysis of PGLYRP3 mRNA level in undifferentiated (undiff) and differentiated (diff) primary keratinocytes (NHEK) and HaCaT cells. (A) representative PCR results. (B) statistical analysis of results from 3 experiments. PGLYRP3 mRNA level in undifferentiated NHEK or HaCaT cells is represented as 1.00 (white bar). The level in differentiated NHEK and HaCaT cells is represented by striped and gray bars, respectively. Data are presented as mean ± SEM.
Ijms 20 00914 g002
Figure 3. EMSA analysis. Nuclear extract of HaCaT cells (10 µg protein) was incubated with radioactively labelled probes containing the analyzed CpG in either an unmethylated or methylated form. Protein binding to each probe was competed with 10, 30 and 50 times excess of the other probe. The DNA-protein complex is indicated by an arrow.
Figure 3. EMSA analysis. Nuclear extract of HaCaT cells (10 µg protein) was incubated with radioactively labelled probes containing the analyzed CpG in either an unmethylated or methylated form. Protein binding to each probe was competed with 10, 30 and 50 times excess of the other probe. The DNA-protein complex is indicated by an arrow.
Ijms 20 00914 g003
Figure 4. Percentage of proteins classified as nucleic acid binding and, within this category, as RNA- or DNA-binding proteins, bound to unmethylated (white bars) and methylated (striped bars). Data obtained in 4 experiments are presented as mean ± SEM.
Figure 4. Percentage of proteins classified as nucleic acid binding and, within this category, as RNA- or DNA-binding proteins, bound to unmethylated (white bars) and methylated (striped bars). Data obtained in 4 experiments are presented as mean ± SEM.
Ijms 20 00914 g004

Share and Cite

MDPI and ACS Style

Sobiak, B.; Leśniak, W. The Effect of Single CpG Demethylation on the Pattern of DNA-Protein Binding. Int. J. Mol. Sci. 2019, 20, 914. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20040914

AMA Style

Sobiak B, Leśniak W. The Effect of Single CpG Demethylation on the Pattern of DNA-Protein Binding. International Journal of Molecular Sciences. 2019; 20(4):914. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20040914

Chicago/Turabian Style

Sobiak, Barbara, and Wiesława Leśniak. 2019. "The Effect of Single CpG Demethylation on the Pattern of DNA-Protein Binding" International Journal of Molecular Sciences 20, no. 4: 914. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20040914

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop