CoChat
Sign Up
Cited by 226Epigenetics and DNA MethylationChronic Kidney Disease and DiabetesGenomics and Chromatin Dynamics

The DOI matches its Crossref record; 35 of the 35 references it cites were found in a registry or an index.

Cytosine methylation changes in enhancer regions of core pro-fibrotic genes characterize kidney fibrosis development

Open access
Ko, Yi-An, Mohtat, Davoud, Suzuki, Masako, Park, Ae Seo Deok, Izquierdo, Maria Concepcion, Han, Sang Youb, Kang, Hyun Mi, Si, Han, Hostetter, Thomas, Pullman, James M, Fazzari, Melissa, Verma, Amit, Zheng, Deyou, Greally, John M, Susztak, Katalin·2013·Genome Biology
DOI 10.1186/gb-2013-14-10-r108PMC PMC4053753
Open this paper in CoChat — free accountOpen PDF

Abstract

Background: One in eleven people is affected by chronic kidney disease, a condition characterized by kidney fibrosis and progressive loss of kidney function. Epidemiological studies indicate that adverse intrauterine and postnatal environments have a long-lasting role in chronic kidney disease development. Epigenetic information represents a plausible carrier for mediating this programming effect. Here we demonstrate that genome-wide cytosine methylation patterns of healthy and chronic kidney disease tubule samples obtained from patients show significant differences. Results: We identify differentially methylated regions and validate these in a large replication dataset. The differentially methylated regions are rarely observed on promoters, but mostly overlap with putative enhancer regions, and they are enriched in consensus binding sequences for important renal transcription factors. This indicates their importance in gene expression regulation. A core set of genes that are known to be related to kidney fibrosis, including genes encoding collagens, show cytosine methylation changes correlating with downstream transcript levels. Conclusions: Our report raises the possibility that epigenetic dysregulation plays a role in chronic kidney disease development via influencing core pro-fibrotic pathways and can aid the development of novel biomarkers and future therapeutics.

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1186/gb-2013-14-10-r108

Cites

35 of 35 references verified0 citation errors · 0 not checkable

Statistics

Figures

Figure 1 Statistically significant cytosine methylation differences in chronic kidney disease. (A) Volcano plot analysis of cytosine methylation differences. The x-axis represents the relative cytosine methylation difference of control (CTL) versus CKD samples, the y-axis represents the negative log2 of the P-value of that locus. The mean P-value and mean difference of 1.3 million loci present on the chips are plotted on the graph. The green and red lines represent the statistical criteria used for further analysis (P-value and fold change, respectively). (B) Hierarchical cluster analysis of the differentially methylated regions. Each column represents changes from one individual kidney sample; blue indicates hypermethylation in CKD, while red represents hypomethylation in CKD. The chart below shows the clinical parameters of the samples: glomerular filtration rate, diabetes status (DM, diabetes mellitus), sex, and age (aged >65 years or <65 years). (C) Gene Ontology analysis of the 1,535 DMRs mapped to unique genes using DAVID gene ontology annotation groups (biological process level 1 annotation).
Figure 1 Statistically significant cytosine methylation differences in chronic kidney disease. (A) Volcano plot analysis of cytosine methylation differences. The x-axis represents the relative cytosine methylation difference of control (CTL) versus CKD samples, the y-axis represents the negative log2 of the P-value of that locus. The mean P-value and mean difference of 1.3 million loci present on the chips are plotted on the graph. The green and red lines represent the statistical criteria used for further analysis (P-value and fold change, respectively). (B) Hierarchical cluster analysis of the differentially methylated regions. Each column represents changes from one individual kidney sample; blue indicates hypermethylation in CKD, while red represents hypomethylation in CKD. The chart below shows the clinical parameters of the samples: glomerular filtration rate, diabetes status (DM, diabetes mellitus), sex, and age (aged >65 years or <65 years). (C) Gene Ontology analysis of the 1,535 DMRs mapped to unique genes using DAVID gene ontology annotation groups (biological process level 1 annotation).CC BY · PubMed Central
Figure 2 External and internal validation of the observed changes. (A) Correlation of the DMRs identified in 26 samples using the HELP method and an external dataset containing 87 human kidney samples analyzed using Illumina Infinium 450K arrays. We found concordant regulation of 1,061 transcripts (98%) from the 1,092 mapped genes using this validation dataset. Of transcripts that showed both differential methylation and expression, 404 (97%) were confirmed. (B) An example of a DMR identified by the HELP assay and confirmed in the validation dataset. This DMR is localized in the intronic region of COLIVA1. (C) Methylation status and (D) gene expression of COLIVA1 in the original HELP dataset in control and CKD samples. (E) MassArray confirmation of methylation status of the COLIVA1 locus (blue is mean ± standard deviation of control samples, red is mean ± standard deviation of CKD samples). (F) Methylation status of the COLIVA1 locus in the validation dataset (Infinium 450K arrays). The data represent the mean differences in absolute methylation levels of individual cytosines at the COLIVA1 locus. (G) Immunohistochemistry of COLIVA1 expression in control (CTL), diabetic (DM), DKD and CKD kidneys.
Figure 2 External and internal validation of the observed changes. (A) Correlation of the DMRs identified in 26 samples using the HELP method and an external dataset containing 87 human kidney samples analyzed using Illumina Infinium 450K arrays. We found concordant regulation of 1,061 transcripts (98%) from the 1,092 mapped genes using this validation dataset. Of transcripts that showed both differential methylation and expression, 404 (97%) were confirmed. (B) An example of a DMR identified by the HELP assay and confirmed in the validation dataset. This DMR is localized in the intronic region of COLIVA1. (C) Methylation status and (D) gene expression of COLIVA1 in the original HELP dataset in control and CKD samples. (E) MassArray confirmation of methylation status of the COLIVA1 locus (blue is mean ± standard deviation of control samples, red is mean ± standard deviation of CKD samples). (F) Methylation status of the COLIVA1 locus in the validation dataset (Infinium 450K arrays). The data represent the mean differences in absolute methylation levels of individual cytosines at the COLIVA1 locus. (G) Immunohistochemistry of COLIVA1 expression in control (CTL), diabetic (DM), DKD and CKD kidneys.CC BY · PubMed Central
Figure 3 Chronic kidney disease differentially methylated regions are localized to kidney-specific enhancer regions. (A) RefSeq annotation of the DMRs. Relative enrichment ratio of the DMRs compared with the representation of the different elements on the methylation microarray. TSS, transcription start site; TTS, transcription termination site. (B) DMRs overlap with regulatory element (chromatin state) annotation maps of renal tubule epithelial cells (HKC8) and adult kidney cortex, indicating that most differentially methylated cytosines are localized to enhancer (yellow and orange) regions in kidney epithelial cells. The color code annotation of the chromatin state map is provided bottom right. (C) Chromatin annotation of the DMRs in five different ENCODE cell lines (H1, human embryonic stem cells; HepG2, hepatocytes; HUVEC, endothelial cells; K562, erythroid cells; NHLF, human lung fibroblasts).
Figure 3 Chronic kidney disease differentially methylated regions are localized to kidney-specific enhancer regions. (A) RefSeq annotation of the DMRs. Relative enrichment ratio of the DMRs compared with the representation of the different elements on the methylation microarray. TSS, transcription start site; TTS, transcription termination site. (B) DMRs overlap with regulatory element (chromatin state) annotation maps of renal tubule epithelial cells (HKC8) and adult kidney cortex, indicating that most differentially methylated cytosines are localized to enhancer (yellow and orange) regions in kidney epithelial cells. The color code annotation of the chromatin state map is provided bottom right. (C) Chromatin annotation of the DMRs in five different ENCODE cell lines (H1, human embryonic stem cells; HepG2, hepatocytes; HUVEC, endothelial cells; K562, erythroid cells; NHLF, human lung fibroblasts).CC BY · PubMed Central
Figure 4 Chronic kidney disease differentially methylated regions are enriched for kidney-specific transcription factor binding sites. (A) The DMR and DHS sites contain consensus sequences. The transcription factor binding site motifs and their statistical enrichment from the de novo searched consensus sequences in DMR and DHS sites. (B) A specific example of an intronic DMR (within the EZR gene). The genomic location of the DMR is at the top, followed by the RefSeq representation of EZR; fetal kidney (FK)-specific DHS tracks (in blue); HKC8 cell-specific H3K4me1 and H3K4me3 tracks; HKC8 cell specific ChromHMM annotation of the locus (yellow, enhancer; red, promoter; green, transcription-associated region; the full color coding key is shown bottom right) - the sequences contain consensus-binding sites for the key kidney transcription factor SIX2/3, with the SIX2/3 binding motif illustrated as a sequence logo plot below; and adult kidney (AK)-specific H3K4me1 (blue) and H3K4me3 (green) tracks [19].
Figure 4 Chronic kidney disease differentially methylated regions are enriched for kidney-specific transcription factor binding sites. (A) The DMR and DHS sites contain consensus sequences. The transcription factor binding site motifs and their statistical enrichment from the de novo searched consensus sequences in DMR and DHS sites. (B) A specific example of an intronic DMR (within the EZR gene). The genomic location of the DMR is at the top, followed by the RefSeq representation of EZR; fetal kidney (FK)-specific DHS tracks (in blue); HKC8 cell-specific H3K4me1 and H3K4me3 tracks; HKC8 cell specific ChromHMM annotation of the locus (yellow, enhancer; red, promoter; green, transcription-associated region; the full color coding key is shown bottom right) - the sequences contain consensus-binding sites for the key kidney transcription factor SIX2/3, with the SIX2/3 binding motif illustrated as a sequence logo plot below; and adult kidney (AK)-specific H3K4me1 (blue) and H3K4me3 (green) tracks [19].CC BY · PubMed Central
Figure 5 Differentially methylated regions correlate with transcript changes. (A) The 4,751 DMRs mapped to 1,092 unique genes that were present in the Affymetrix arrays. There were 415 transcripts that showed differences both in their methylation status and their expression in CKD samples. The RefSeq-based locations of the DMRs are also shown. While most differentially methylated regions localize to gene body regions, they also show correlation with the expression of many of those genes. Not only are the 415 transcripts differentially expressed, they also show differences in their cytosine methylation profiles as well. (B) DAVID-based gene ontology annotation of the 415 transcripts. (C) Network chart of the genes that are both differentially expressed and methylated). (D) Methylation and gene expression level of key molecules (RUNX3, RARB, SMAD6) identified by the network analysis in control (CTL) and CKD samples.
Figure 5 Differentially methylated regions correlate with transcript changes. (A) The 4,751 DMRs mapped to 1,092 unique genes that were present in the Affymetrix arrays. There were 415 transcripts that showed differences both in their methylation status and their expression in CKD samples. The RefSeq-based locations of the DMRs are also shown. While most differentially methylated regions localize to gene body regions, they also show correlation with the expression of many of those genes. Not only are the 415 transcripts differentially expressed, they also show differences in their cytosine methylation profiles as well. (B) DAVID-based gene ontology annotation of the 415 transcripts. (C) Network chart of the genes that are both differentially expressed and methylated). (D) Methylation and gene expression level of key molecules (RUNX3, RARB, SMAD6) identified by the network analysis in control (CTL) and CKD samples.CC BY · PubMed Central
Figure 6 Regulation of transcripts by a DNA methyltransferase inhibitor in in vitro cultured human tubular epithelial cells. Gene ontology terms of transcripts showing differential expression in the decitabine-treated cells. (A) Illustration of regions that showed differential methylation of cultured HKC8 cells treated with 0.5 μM decitabine (5'DAC). CTL, control. (B) The interconnected network analysis highlighted the differential expression of cell adhesion and developmental pathways. These genes are also differentially expressed and methylated in the original CKD dataset. GO, gene ontology.
Figure 6 Regulation of transcripts by a DNA methyltransferase inhibitor in in vitro cultured human tubular epithelial cells. Gene ontology terms of transcripts showing differential expression in the decitabine-treated cells. (A) Illustration of regions that showed differential methylation of cultured HKC8 cells treated with 0.5 μM decitabine (5'DAC). CTL, control. (B) The interconnected network analysis highlighted the differential expression of cell adhesion and developmental pathways. These genes are also differentially expressed and methylated in the original CKD dataset. GO, gene ontology.CC BY · PubMed Central
Figure 7 Gene body cytosine methylation changes drive gene expression differences. RUNX1 methylation and gene expression were examined only (A,B) In the original discovery dataset, the gene body region of RUNX1 was hypomethylated (A) and the corresponding transcript level was increased (B) in the CKD (discovery) dataset. CTL, control. (C,D) The differential methylation (C) and expression (D) of RUNX1 in the DKD replication dataset. (E,F) Transcript levels are increased (F) in vitro in cultured tubules after decreasing the methylation level of the locus following 0.5 μM decitabine (DAC) treatment (E). (G) Genomic representation of the RUNX1 locus showing DMRs in the DKD dataset and in the CKD dataset. Different tracks are shown for the RUNX1 locus, including RefSeq gene, DMRs in the DKD dataset, DMRs in the CKD cells, and histone ChIP-seq data for H3K4me1 and H3K4me3 for adult kidney cortex and DHS sites from fetal kidneys. In addition, ENCODE-based transcription factor binding sites are also shown.
Figure 7 Gene body cytosine methylation changes drive gene expression differences. RUNX1 methylation and gene expression were examined only (A,B) In the original discovery dataset, the gene body region of RUNX1 was hypomethylated (A) and the corresponding transcript level was increased (B) in the CKD (discovery) dataset. CTL, control. (C,D) The differential methylation (C) and expression (D) of RUNX1 in the DKD replication dataset. (E,F) Transcript levels are increased (F) in vitro in cultured tubules after decreasing the methylation level of the locus following 0.5 μM decitabine (DAC) treatment (E). (G) Genomic representation of the RUNX1 locus showing DMRs in the DKD dataset and in the CKD dataset. Different tracks are shown for the RUNX1 locus, including RefSeq gene, DMRs in the DKD dataset, DMRs in the CKD cells, and histone ChIP-seq data for H3K4me1 and H3K4me3 for adult kidney cortex and DHS sites from fetal kidneys. In addition, ENCODE-based transcription factor binding sites are also shown.CC BY · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Citations

Cites

Cited by

Sorted by citations

See all 226

Related

See all

Metadata

Authors
Ko, Y.-A., Mohtat, D., Suzuki, M., Park, A. S. D., Izquierdo, M. C., Han, S. Y., Kang, H. M., Si, H., Hostetter, T., Pullman, J. M., Fazzari, M., Verma, A., Zheng, D., Greally, J. M., & Susztak, K.
Affiliations
  • University of Pennsylvania (US)
  • Albert Einstein College of Medicine (US)
  • Case Western Reserve University (US)
  • Montefiore Medical Center (US)
Corresponding author
Greally, John M
ORCID
Venue
Genome Biology
Publisher
Springer Science and Business Media LLC
Published
2013-12-10
Volume
14
Issue
10
Pages
R108-R108
ISSN
  • 1465-6906 (ISSN-L)
  • 1474-760X
DOI
10.1186/gb-2013-14-10-r108
PMC
PMC4053753
OpenAlex
W2130816592
Cited by
225
References
35
Citations by year
2026: 9 · 2025: 10 · 2024: 19 · 2023: 9 · 2022: 10 · 2021: 14 · 2020: 16 · 2019: 24 · 2018: 29 · 2017: 20 · 2016: 16 · 2015: 30 · 2014: 16 · 2013: 3
FWCI
9.71
Type
journal article
Language
en
Open access
Licence
Locations
  • Genome biology (publishedVersion, cc-by)
  • Europe PMC (PubMed Central) (publishedVersion, cc-by)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
Funders
  • National Institutes of Health
  • National Institute of Diabetes and Digestive and Kidney Diseases (R01DK087635)
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Epigenetics and DNA Methylation
  • Health Sciences › Medicine › Nephrology › Chronic Kidney Disease and Diabetes
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Genomics and Chromatin Dynamics
Keywords
cytosine methylation, renal fibrosis, chronic kidney disease, differentially methylated regions, enhancers, epigenetic dysregulation, collagen genes, gene expression regulation, developmental programming, biomarkers
MeSH terms
Aged, Base Sequence, Binding Sites, Cytosine, Enhancer Elements, Genetic, Fibrosis, Gene Expression Regulation, Humans, Kidney, Middle Aged, Transcription Factors, Transcription, Genetic, Reproducibility of Results, Cluster Analysis, Case-Control Studies, DNA Methylation, Gene Expression Profiling, Renal Insufficiency, Chronic, Gene Regulatory Networks, Position-Specific Scoring Matrices, Nucleotide Motifs

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review