Open Access

DNA hypomethylation of a transcription factor binding site within the promoter of a gout risk gene NRBP1 upregulates its expression by inhibition of TFAP2A binding

Clinical EpigeneticsThe official journal of the Clinical Epigenetics Society20179:99

https://doi.org/10.1186/s13148-017-0401-z

Received: 19 January 2017

Accepted: 5 September 2017

Published: 15 September 2017

Abstract

Background

Genome-wide association studies (GWASs) have identified dozens of loci associated with gout, but for most cases, the risk genes and the underlying molecular mechanisms contributing to these associations are unknown. This study sought to understand the molecular mechanism of a common genetic variant, rs780093, in the development of gout, both in vitro and in vivo.

Results

Nuclear receptor binding protein 1 (NRBP1), as a gout risk gene, and its regulatory region, 72 bp upstream of the transcription start site, designated as B1, were identified through integrative analyses of genome-wide genotype and DNA methylation data. We observed elevated NRBP1 expression in human peripheral blood mononuclear cells (PBMCs) from gout patients. In vitro luciferase reporter and protein pulldown assay results showed that DNA methylation could increase the binding of the transcription factor TFAP2A to B1, leading to suppressed gene expression. There results were further confirmed by in vivo bisulfite pyrosequencing showing that hypomethylation on B1 is associated with increased NRBP1 expression in gout patients.

Conclusions

Hypomethylation at the promoter region of NRBP1 reduces the binding of TFAP2A and thus leads to elevated NRBP1 expression, which might contribute to the development of gout.

Keywords

Gout Uric acid DNA methylation NRBP1 TFAP2A

Background

Gout is a complex disorder caused by deposition of monosodium urate (MSU) crystals within joints with a prevalence of 1.14% in adults in China [1]. The male to female sex ratio is generally around 3–4 to 1 [2]. Initial symptom includes severely painful episodes of peripheral joint synovitis, but it can eventually lead to joint damage and deformity, chronic usage-related pain, and subcutaneous tophus deposition. Central to the development of gout is elevated serum uric acid concentrations [3]. Gout is also associated with many conditions that affect longevity and well-being, such as hypertension, diabetes mellitus, metabolic syndrome, and renal and cardiovascular disease [49]. In particular, gout is increasingly recognized as an independent cardiovascular risk factor [2, 10]. Over the past decade, genome-wide association studies (GWASs) and subsequent meta-analyses have led to a considerable expansion in the knowledge of common genetic loci that are associated with gout and elevated serum urate level [1115]. However, apart from genes coding uric acid transporters, the mechanisms of how most of genetic variances identified from GWAS regulating urate level and developing into gout still remain poorly understood. As a common human disease, it is possible that both genetic predisposition and environmental exposure sharp its risk and contribute to its pathogenesis [16].

Epigenetics, such as histone modifications, DNA methylation, higher chromatin structure, and noncoding RNAs, can be influenced by genetic as well as environmental factors [17], and thus it can be a potential mechanism leading towards gout [1820]. DNA methylation, an ancient and crucial epigenetic modification, is the heritable methylation of cytosines in the context of cytosine-guanine dinucleotides (CpGs) without changing the DNA sequence itself [21]. It is commonly believed that DNA methylation can disrupt the interactions between transcription factors (TFs) and DNA either directly [22] or indirectly by recruiting methyl CpG-binding proteins, such as MeCP2, that occupy the methylated promoters and compete for the TF binding sites [22]. However, some recent findings show that TFs could specifically bind to methylated CpG sites with specific DNA sequences (DNA motifs) and promote gene expression in vitro [23].

In the past decade, it has been shown that DNA methylation may play a role in many human autoimmune diseases, such as systemic lupus erythematosus, rheumatoid arthritis, multiple sclerosis, systemic scleroderma, and Sjogren’s syndrome [24]. The identification of dysregulation of DNA methylation may inspire the discovery of other uncharacterized mechanisms. Understanding the molecular mechanisms involved in the pathophysiology of autoimmune diseases is essential for the discovering of diagnostic biomarkers, as well as the introduction of effective and target-directed therapies [25]. However, there is no research regarding the role of epigenetics in gout up to now.

Here, by integrating genotype with DNA methylation, we identified a gout risk gene, nuclear receptor binding protein 1 (NRBP1). We demonstrated that hypomethylation of its promoter region, 72 bp upstream of NRBP1 transcription start site (designated as B1 thereafter), is associated with increased gene expression both in vitro and in vivo. Moreover, gout-associated increased expression of NRBP1 is regulated through methylation-dependent TFAP2A binding to the B1 region.

Methods

Study design

The study design is shown in Additional file 1: Figure S1.

Clinical samples

The study was approved by the Ethics Committees of Huashan Hospital, Fudan University. PBMCs were harvested from males with gout before any medical treatment and matched healthy volunteers. Gout patients were diagnosed and met the requirements of the 1977 American Rheumatism Association preliminary criteria [26]. Gout patients were also with primary hyperuricemia, defined as serum uric acid concentrations higher than 7 mg/dL [27]. The inclusion criteria for gout patients are (1) no other diseases except for gout, (2) no treatments in the last 3 months, (3) male, (4) age between 18 and 50, and (5) never smoker. Clinical characteristics of cases and controls are shown in Table 1. Additionally, no study subjects have ever suffered from cancer. PBMCs from 17 gout patients and 15 matched controls were used for measurement of NRBP1 RNA level and methylation of B1. PBMCs from an additional set of 7 gout patients and 7 matched controls were used for measurement of NRBP1 protein level.
Table 1

Clinical characteristics of study subjects

Characteristic

Cases (n = 15)

Controls (n = 17)

P value

Age (years)

39.33 ± 5.94

36.12 ± 7.66

0.12

Male (n)

15

17

1

Serum uric acid (μmol/L)

535.80 ± 38.31

320.65 ± 66.80

1.42E−11

Creatinine (μmol/L)

68.80 ± 28.18

56.94 ± 12.00

0.21

Body mass index (kg/m2)

24.81 ± 3.64

22.16 ± 3.52

0.04

Cholesterol (mg/dL)

197.20 ± 99.09

134.24 ± 49.02

0.03

Triglycerides (mg/dL)

149.73 ± 33.21

120.29 ± 43.33

0.04

Cancer (n)

0

0

1

Smoker (n)

0

0

1

Values are mean ± SD or numbers. P value for the difference between cases and controls

Integrative analyses of genotype and DNA methylation

Data collection and processing for integrative genotype and DNA methylation data were performed as described previously with minor modifications [17]. Briefly, DNA methylation data from human blood and genotype data from the same cohort was obtained from the Gene Expression Omnibus (GEO) with accession number GSE42861. Disease-associated SNPs were collected from the NHGRI Catalog of Published GWASs [1115], and a total of 105 SNPs, identified previously to be associated with gout or urate levels, were included for downstream genotype-methylation analyses. Associations between genotype and DNA methylation were evaluated using an additive minor-allele dosage model and SNP-CpG pairs with a stringent Bonferroni-adjusted P value less than 0.05 were determined to be significant, as suggested previously [17].

Quantitative real-time PCR (qPCR)

Total cellular RNA of PBMCs was isolated with TRIzol reagent (Invitrogen, Carlsbad, USA). cDNA was generated by reverse transcription (PrimeScript RT Reagent kit; Takara), and the expression of NRBP1 was quantified by qPCR using LightCycler 480 SYBR Green I Master Mix (Roche) normalized to GAPDH (more details in Additional file 2: File 1a). A non-template negative control (water) and a standard curve were conducted for each assay. The efficiency of amplification was determined based on the standard curve for mRNA quantification. All assays were performed in triplicates, and no obvious outlier was observed based on the SD.

Generation of methylated DNA and luciferase reporter

Single-stranded DNA (ssDNA) were synthetized in Sangon Biotech (Shanghai, China) and annealed to generate double-stranded DNA (ds-DNA), which were then purified using Illustra MicroSpin G-25 Columns (GE Healthcare) according to the manufacturer’s protocol (Additional file 2: File 1b) [28]. Methylation of ds-DNA or luciferase reporter constructs was generated by their treatment with M.SssI CpG methyltransferase (Zymo Research) as described previously (Additional file 2: File 1b) [29]. The DNA were then recovered by ethanol precipitation. The methylation status was confirmed by bisulfite pyrosequencing.

Protein pulldown assay

Biotinylated ssDNA were synthetized in Sangon Biotech (Shanghai, China) and annealed to generate dsDNA as described previously. Methylated or unmethylated biotinylated dsDNA probes with eight repeats of B1 were generated as above. GST-tagged human transcriptional factor TFAP2A was generated and purified from yeast as described previously [30]. The streptavidin magnetic beads (Roche, USA) were incubated with 1 pmol of the methylated or unmethylated biotinylated probes in the binding buffer (Additional file 2: File 1c) at room temperature for 10 min. The beads were then washed with washing buffer (Additional file 2: File 1c) for three times followed by incubation with 10 pmol of purified TFAP2A protein in PBS at 4 °C overnight. Beads were enriched on magnets and unbound supernatant were collected for Western blot. The beads, together with 10% input and 10% supernatant, were heated at 100 °C for 10 min and loaded on a 10% polyacrylamide gel by electrophoresis. The proteins were detected using antibody against TFAP2A (1:500; ab52222, Abcam) and visualized using luminescent image analyzer Image Quant Las 4000 mini (GE Healthcare Japan, Chiba, Japan).

293T cell culture

293T cells were maintained at 37 °C in 5% CO2 and grown in Dulbecco’s modified Eagle’s medium (GIBCO) supplemented with 10% fetal bovine serum (GIBCO). One day prior to transfection, 293T cells were seeded at a density of 100,000 cells per well on 48-well plates.

Dual luciferase transcriptional reporter assay

The CpG-free luciferase reporter vector, pCpGL, is used for detecting promoter’s activity in dual luciferase transcriptional reporter assay [23]. Ds-DNA containing eight repeats of B1 (GGCGCAAGG) were synthesized and subcloned into pCpGL promoter region to generate pCpGL-8X-B1. Methylation of the luciferase reporter was carried out with M.SssI as described previously. The coding sequence (CDS) of TFAP2A was cloned into the expression vector, pCAGIG, with DNA recombination by using pEASY-Uni Seamless Cloning and Assembly Kit (Transgen, Shanghai, China). The following primers were used for the construct cloning; forwards: 5′-GTCTCATCATTTTGGCAAAGAATTCATCACAAGTTTGTACAAAAAAGCAG-3′ and reverse: 5′-ACGTAGCGGCCGCGATATCCTCGAGGCTTTCACCACTTTGTACAAGAAAG-3′. 293T cells were co-transfected with three constructs: pCpGL-8X-B1, pCAGIG expressing the TFAP2A at various concentrations, and pTK-RL (Promega, Madison, WI) using Lipofectamine 2000 (Life Technologies, USA). Cells were harvested 48 h post-transfection for luciferase reporter assay using the dual-luciferase reporter assay system (Promega), and signals were recorded by calculating firefly luciferase activity normalized to renilla activity. All assays were performed in triplicates and no obvious outlier was observed based on the SD.

Sodium bisulfite pyrosequencing

Genomic DNA of PBMCs were extracted using the DNeasy Blood and Tissue Kit (Qiagen, Hilden, Germany). A total of 400 ng of genomic DNA were treated with sodium bisulfite using the EZ DNA Methylation-Gold Kit (Zymo Research). Twenty-five nanograms of bisulfite-converted genomic DNA was amplified with unbiased nested PCR reaction using ExTaq DNA Polymerase (Takara, Kyoto, Japan). The primer sequences to amplify B1 region at the promoter of NRBP1 were as follows: outer primers forward: 5′-TTATTATTGAATGATAATTTTAATGAGTT-3′ and reverse: 5′-CCTAAACTACTAAATAAACAAAACC-3′; and inner primers forward: 5′-GTAGAATTATTTGGGGTATTTGGAT-3′ and reverse: 5′ biotin-AACCCTCTTTTCCCTAAAC-3′. To quantify the percentage of methylated cytosine in each CpG site, amplified DNA were sequenced using a pyrosequencing system (PyroMark Q96, Qiagen) [31]. This method treats each individual CpG site as a C/T polymorphism and generates quantitative data for the relative proportion of the methylated versus the unmethylated allele. The sequencing primer 5′-GGTGGGGTGGATAGAGA-3′ was designed as reverse run, and the percentage of DNA methylation on each measured CpG site was generated by Pyro Q-CpG methylation software (Biotage). A non-template negative control (water) and a standard curve were included for each run.

Statistical analysis

All results were expressed as means ± SEM of at least three independent experiments. Comparisons between groups were evaluated by using an unpaired two-tailed Student’s t test. The association between DNA methylation and gene expression was calculated by a linear regression model. A P value less than 0.05 was considered statistically significant.

Results

Identification of the gout risk gene NRBP1, by integrative analyses of GWAS with DNA methylation data

Our previous work has showed that DNA methylation can be influenced by genotype and integrative genotype-methylation analyses can help us to fine map the epigenetic variants that might be responsible for the disease phenotype [17]. Given this, we want to investigate how altered DNA methylation level may play a role in disease etiology, specifically in gout. To do this, we first examined the genotype-methylation associations for 105 SNPs, identified previously to be associated with gout or urate level, from the NHGRI Catalog of Published GWASs [1115]. We identified six disease-associated SNPs that were significantly associated with altered DNA methylation level (Table 2). Among them, only one CpG site, cg05102552, is located in the promoter region of a gene. Since it is well known that promoter methylation can regulate the gene expression, we thus decided to focus on this region for the downstream study.
Table 2

CpG sites whose methylation are controlled by gout-associated GWAS SNPs

GWAS SNP

CpG

Chromosome

SNP position (hg19)

CpG position (hg19)

Gene close to the CpG site

P value (geno vs. meth)

rs780093

cg05102552

chr2

27742603

27650867

NRBP1

3.13E−12

rs13129697

cg00071950

chr4

9926967

10020882

SLC2A9

3.87E−12

rs3775948

cg00071950

chr4

9997112

10020882

SLC2A9

< 2E−16

rs717615

cg00071950

chr4

10104670

10020882

SLC2A9

< 2E−16

rs1165205

cg15691649

chr6

25870542

25882328

SLC17A3

< 2E−16

rs505802

cg19131476

chr11

64357072

64387923

NRXN2

< 2E−16

The methylation level of cg05102552 is strongly associated with a SNP, rs780093 (Fig. 1a), which has been shown to associate with gout in both Caucasian [12] and Han Chinese populations [32]. Cg05102552 site is 552 bp upstream of NRBP1 transcription start site and is in a region with DNase I hypersensitive sites (DHSs) [33, 34] and high level of active histone modifications, such as H3K4me1, H3K4me3, and H3K27Ac (Fig. 1b) [35], suggesting that it may be the regulatory region for NRBP1 expression. Furthermore, both normalized RNA level (gout versus control 1.425 ± 0.0656 versus 1.185 ± 0.0432 (mean ± SEM); P value = 0.004) and protein level of NRBP1 in PBMCs were found to be significantly elevated in gout patients (Fig. 2), indicating that NRBP1 may be a gout risk gene.
Fig. 1

Identification of a gout risk gene and a potential TFAP2A binding site at its promoter region. a The association between a gout-associated GWAS SNP, rs780093, and DNA methylation level of cg05102552 from human peripheral blood mononuclear cells (PBMCs). Each dot represents an individual and average methylation level for each genotype group is indicated by red bar. Statistical significance for the association was evaluated with additive minor-allele dosage model and P value was indicated in the bottom. b Cg05102552 site is 552 bp upstream of NRBP1 transcription start site (TSS) and is in a region with DNase I hypersensitive sites (DHSs) and high level of active histone modifications, such as H3K4me1, H3K4me3, and H3K27Ac. The potential TFAP2A binding region, B1, was highlighted in the bottom

Fig. 2

NRBP1 is overexpressed in PBMCs from gout patients. a Relative NRBP1 mRNA expression in freshly isolated PBMCs from 15 gout patients and 17 healthy controls detected by quantitative real-time PCR (qPCR). The expression level of NRBP1 is normalized to GAPDH and is represented as mean ± SEM. (P value = 0.004, Student’s t test, unpaired, two-sided). b NRBP1 protein expression in freshly isolated PMBCs from 7 gout patients and healthy controls detected by Western blotting

Identification of a potential TFAP2A binding site at the promoter region of NRBP1

We next investigated how methylation at the promoter region of NRBP1 may regulate its gene expression. We overlapped this region with the Encyclopedia of DNA Elements (ENCODE) database and observed that, in addition to epigenetic marks associated with accessible chromatin, this region is also occupied by transcription factors, such as TFAP2A (Fig. 1b) [36]. Moreover, the DNA sequence in this region contains a motif highly similar to a known TFAP2A binding site (Fig. 1b), identified previously with protein microarrays in vitro [23]. We thereafter designated this potential TFAP2A binding site, 72 bp upstream of NRBP1 transcription start site, as B1.

NRBP1 expression is regulated by methylation-dependent TFAP2A binding to its promoter region in vitro

To determine whether B1 is a potential TFAP2A binding site important for NRBP1 expression, we performed a cellular-based dual-luciferase reporter assay. We observed significant reduction of reporter gene expression in the presence of TFAP2A, and this reduction was TFAP2A dose dependent (Fig. 3a), indicating that B1 is indeed a TFAP2A regulatory motif for gene expression. Moreover, the reduction of gene expression was much more enhanced after we methylated B1 region on the luciferase reporter (Fig. 3a), suggesting that B1 can regulate gene expression in a methylation-dependent manner in the presence of TFAP2A. To further interrogate whether this dose- and methylation-dependent inhibition of gene expression is mediated through direct TFAP2A binding to the B1 motif, we conducted a pulldown experiment and observed that purified TFAP2A protein could specifically bind to the B1 motif and this binding was much stronger if the CpG dinucleotide within B1 was methylated (Fig. 3b). Taken together, these results suggest that the transcription factor, TFAP2A, can specifically bind to the B1 region at the promoter of NRBP1 and inhibit gene expression in vitro, and DNA methylation on B1 can further enhance the binding of TFAP2A and the inhibition of gene expression.
Fig. 3

Increased DNA-binding of TFAP2A to methylated B1 inhibits gene transcription in vitro. a TFAP2A exhibited methylation-dependent inhibition of luciferase activity on B1. Dual-luciferase reporter assays were performed to evaluate the transcriptional activity of the methylated or unmethylated B1 by TFAP2A in 293T cells, with firefly luciferase activity normalized to renilla activity. Data are represented as mean ± SEM (**P value < 0.001, ****P value < 0.00001, Student’s t test, unpaired, two-sided). b DNA methylation on B1 can increase its direct interaction with TFAP2A. Purified TFAP2A was incubated with streptavidin bead-bound, biotinylated B1 with or without DNA methylation modification. Bead-bound proteins were fractionated by SDS-PAGE, and TFAP2A was detected by immunoblotting. Portions (10%) of the purified TFAP2A (input) or unbound fractions were assayed in parallel

Hypomethylation of B1 is associated with increased NRBP1 expression in gout patients

In order to test that increased NRBP1 expression we observed in gout patients is regulated through methylation-dependent TFAP2A binding to NRBP1 promoter region, we evaluated the DNA methylation level of the CpG dinucleotide on B1 using PBMCs of gout patients (n = 15) and healthy controls (n = 17). We observed that DNA methylation on B1 was significantly lower from gout patients compared to healthy controls (gout versus control 3.067 ± 0.182% versus 4.118 ± 0.2829%; P value = 0.005) (Fig. 4a). Moreover, there is a significant negative association between NRBP1 gene expression and DNA methylation on B1 (P value = 0.011) (Fig. 4b), which is consistent with our in vitro findings that hypomethylation on B1 is associated with increased gene expression. In addition to gout disease status, we also observed a marginal significant association between DNA methylation of B1 and serum uric acid (P value = 0.08) (Additional file 1: Figure S2a) and a significant association between NRBP1 expression and serum uric acid (P value = 0.03) (Additional file 1: Figure S2b). Thus, these results support the idea that hypomethylation of B1 leads to increased NRBP1 expression in gout patients.
Fig. 4

Hypomethylation of B1 is associated with increased NRBP1 expression in gout patients. a Bisulfite pyrosequencing showed DNA methylation on B1 at the promoter of NRBP1 is significantly lower in gout patients (n = 15) comparing to healthy controls (n = 17). Data are represented as mean ± SEM. (P value = 0.005, Student’s t test, unpaired, two-sided). b A significant negative association between NRBP1 gene expression and DNA methylation on B1 in the tested subjects (P value = 0.0113)

Discussion

In summary, we have identified a gout risk gene, NRBP1, using integrative analyses of genotype and DNA methylation. Through experiments both in vitro and in vivo, we demonstrated that increased NRBP1 expression in gout patients may be regulated through methylation-dependent binding of TFAP2A to the B1 region, 72 bp upstream of NRBP1 transcription start site (Fig. 5).
Fig. 5

Proposed model for methylation-dependent TFAP2A binding and regulation of NRBP1 expression in gout. a In healthy subjects, TFAP2A binds to methylated B1, which suppresses the transcription of NRBP1. b In contrast, hypomethylation of B1 at the promoter region of NRBP1 abrogates TFAP2A binding, which results in elevated NRBP1 expression in gout

Our data showed that even though TFAP2A can bind to unmethylated B1 (Fig. 3b) and inhibit gene expression (Fig. 3a) when B1 is not methylated, the binding between DNA and TFAP2A was much enhanced when B1 is methylated, and this leads to greater reduction of gene expression. However, since we only investigated the role of B1 methylation in vitro, whether hypomethylation of B1 can lead to reduced gene expression in vivo is still unknown. The direct functionality of cg05102552 methylation and its role in regulating NRBP1 gene expression still needs to be investigated in the future.

In our current work, we have identified and shown that NRBP1 may be a risk gene in gout. However, the understanding of downstream pathways explaining detailed molecular mechanisms awaits further study. NRBP1, previously described to play a role in tumor suppression, cellular homoeostasis, and protein regulation [37], was recently reported to be associated with an autoimmune disease, Takayasu’s arteritis [38]. Studies have shown that NRBP1 can regulate the Wnt/β-catenin signaling pathway [39, 40] and affect the expression of ATP binding cassette subfamily G member 2 (ABCG2) [41], a critical urate efflux transporter in kidney proximal tubule and intestine [42]. Considering the important role of ABCG2 in gout, it is worth further investigating the molecular mechanism of NRBP1 in gout development.

Even though we showed that there is a statistically significant difference for DNA methylation of B1 between gout patients and healthy controls, the effect size for the difference is small. This is actually not uncommon to observe relatively small methylation difference associated with human common diseases [4345]. Moreover, methylation level for the CpG sites around the B1 region is also decreased in gout patients (Additional file 1: Figure S3), suggesting that this is a differentially methylated region (DMR). Considering methylation level in the neighboring CpGs is usually correlated [17], this provides an additional evidence supporting hypomethylation of B1 in gout. Additionally, since DNA methylation level analyzed in current study was from PBMCs, which may not be the directly relevant tissue for gout etiology, methylation difference may be much severer in tissues directly involved in gout, such as leucocytes and monocytes/macrophages in affected joints or tissues in kidney and intestine that excrete uric acid. Nevertheless, methylation difference in blood can still reflect the epigenetic regulation involving in the disease mechanism and may be used as a potential biomarker for early diagnosis or intervention before gout attack and severe complication [24].

In spite of these limitations above, to the best of our knowledge, this is the first report showing how DNA methylation may affect the pathogenesis of gout. Contradictory to the conventional wisdom that methylated CpG dinucleotide at the promoter can abrogate transcription factor binding and lead to transcriptional silencing, we demonstrated that DNA methylation can also inhibit gene expression by facilitating transcription factor binding to the promoter. Given DNA methylation is reversible, this work provides a promising target for potential gout therapy.

Conclusions

In conclusion, we identified a gout risk gene, NRBP1, by integrating genome-wide genotype with DNA methylation data. We demonstrated that hypomethylation of its promoter region, B1, 72 bp upstream of NRBP1 transcription start site, is associated with increased gene expression both in vitro and in vivo. Moreover, gout-associated increased NRBP1 expression is regulated through methylation-dependent TFAP2A binding to the B1 region, which might be involved in the pathogenesis of gout.

Abbreviations

ABCG2

ATP binding cassette subfamily G member 2

CpGs: 

Cytosine-guanine dinucleotides

DHSs: 

DNase I hypersensitive sites

DMR: 

Differentially methylated region

ENCODE: 

Encyclopedia of DNA Elements

GEO: 

Gene Expression Omnibus

GWASs: 

Genome-wide association studies

MSU: 

Monosodium urate

NRBP1

Nuclear receptor binding protein 1

PBMCs: 

Peripheral blood mononuclear cells

TF: 

Transcription factor

Declarations

Acknowledgements

The CpG-free luciferase reporter vector, pCpGL, was a gift from Dr. Zhu at the Johns Hopkins University School of Medicine. We are grateful to all the participants involved in this study.

Funding

This work was supported by the National Natural Science Foundation of China (Grant Nos. 81701594, 31471212, 81373213, and 81671588) and the Ministry of Science and Technology 973 (Grant No. 2015CB910401).

Availability of data and materials

The datasets supporting the conclusions of this article are included within the article (and its Additional files 1 and 2).

Authors’ contributions

ZHZ collected clinical samples and performed the experiments, generated the figures and table, and wrote the manuscript. WDM and XXZ participated in the revision of the manuscript. PRL performed bisulfite pyrosequencing and the statistical analyses. YL and HJZ designed and supervised the study and revised the manuscript. All authors read and approved the final manuscript.

Ethics approval and consent to participate

The study was approved by the Ethics Committees of Huashan Hospital, Fudan University.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

(1)
Division of Rheumatology and Immunology, Huashan Hospital, Fudan University
(2)
Institute of Rheumatology, Immunology and Allergy, Fudan University
(3)
The Ministry of Education Key Laboratory of Metabolism and Molecular Medicine, Department of Biochemistry and Molecular Biology, School of Basic Medical Sciences, Fudan University

References

  1. Miao Z, Li C, Chen Y, Zhao S, Wang Y, Wang Z, et al. Dietary and lifestyle changes associated with high prevalence of hyperuricemia and gout in the Shandong coastal cities of Eastern China. J Rheumatol. 2008;35:1859–64.PubMedGoogle Scholar
  2. Kuo C-F, Grainge MJ, Zhang W, Doherty M. Global epidemiology of gout: prevalence, incidence, and risk factors. Nat Rev Rheumatol. 2015;11:649–62.View ArticlePubMedGoogle Scholar
  3. Dalbeth N, Merriman TR, Stamp LK. Gout. Lancet. 2016;388:2039–52.View ArticlePubMedGoogle Scholar
  4. Keenan T, Zhao W, Rasheed A, Ho WK, Malik R, Felix JF, et al. Causal assessment of serum urate levels in cardiometabolic diseases through a mendelian randomization study. J Am Coll Cardiol. 2016;67:407–16.View ArticlePubMedPubMed CentralGoogle Scholar
  5. Palmer TM, Nordestgaard BG, Benn M, Tybjaerg-Hansen A, Davey Smith G, Lawlor D a, et al. Association of plasma uric acid with ischaemic heart disease and blood pressure: mendelian randomisation analysis of two large cohorts. BMJ. 2013;347:f4262.View ArticlePubMedPubMed CentralGoogle Scholar
  6. White J, Sofat R, Hemani G, Shah T, Engmann J, Dale C, et al. Plasma urate concentration and risk of coronary heart disease: a mendelian randomisation analysis. Lancet Diabetes Endocrinol. 2016;4:327–36.View ArticlePubMedPubMed CentralGoogle Scholar
  7. Hughes K, Flynn T, de Zoysa J, Dalbeth N, Merriman TR. Mendelian randomization analysis associates increased serum urate, due to genetic variation in uric acid transporters, with improved renal function. Kidney Int. 2014;85:344–51.View ArticlePubMedGoogle Scholar
  8. Sluijs I, Holmes MV, Van Der Schouw YT, Beulens JWJ, Asselbergs FW, Huerta JM, et al. A Mendelian randomization study of circulating uric acid and type 2 diabetes. Diabetes. 2015;64:3028–36.View ArticlePubMedGoogle Scholar
  9. Kleber ME, Delgado G, Grammer TB, Silbernagel G, Huang J, Kramer BK, et al. Uric acid and cardiovascular events: a Mendelian randomization study. J Am Soc Nephrol. 2015;26:2831–8.View ArticlePubMedPubMed CentralGoogle Scholar
  10. So A. Epidemiology: gout—bad for the heart as well as the joint. Nat Publ Gr. 2010;6:386–7.Google Scholar
  11. Dehghan A, Köttgen A, Yang Q, Hwang SJ, Kao WL, Rivadeneira F, et al. Association of three genetic loci with uric acid concentration and risk of gout: a genome-wide association study. Lancet. 2008;372:1953–61.View ArticlePubMedPubMed CentralGoogle Scholar
  12. Yang Q, Köttgen A, Dehghan A, Smith AV, Glazer NL, Chen H, et al. Multiple genetic loci influence serum urate and their relationship with gout and cardiovascular disease risk factors. Circ Cardiovasc Genet. 2010;3:523–30.View ArticlePubMedPubMed CentralGoogle Scholar
  13. Sulem P, Gudbjartsson DF, Walters GB, Helgadottir HT, Helgason A, Gudjonsson S a, et al. Identification of low-frequency variants associated with gout and serum uric acid levels. Nat Genet. 2011;43:1127–30.View ArticlePubMedGoogle Scholar
  14. Köttgen A, Albrecht E, Teumer A, Vitart V, Krumsiek J, Hundertmark C, et al. Genome-wide association analyses identify 18 new loci associated with serum urate concentrations. Nat Genet. 2013;45:145–54.View ArticlePubMedGoogle Scholar
  15. Li C, Li Z, Liu S, Wang C, Han L, Cui L, et al. Genome-wide association analysis identifies three new risk loci for gout arthritis in Han Chinese. Nat Commun. 2015;6:7041.View ArticlePubMedPubMed CentralGoogle Scholar
  16. Feinberg AP, Fallin MD. Epigenetics at the crossroads of genes and the environment. JAMA. 2015;314:1129–30.View ArticlePubMedGoogle Scholar
  17. Liu Y, Li X, Aryee MJ, Ekström TJ, Padyukov L, Klareskog L, et al. GeMes, clusters of DNA methylation under genetic control, can inform genetic and epigenetic analysis of disease. Am J Hum Genet. 2014;94:485–95.View ArticlePubMedPubMed CentralGoogle Scholar
  18. Gee Teng G, Pan A, Yuan JM, Koh WP. Cigarette smoking and the risk of incident gout in a prospective cohort study. Arthritis Care Res. 2016;68:1135–42.View ArticleGoogle Scholar
  19. Krishnan E, Lessov-Schlaggar CN, Krasnow RE, Swan GE. Nature versus nurture in gout: a twin study. Am J Med. 2012;125:499–504.View ArticlePubMedGoogle Scholar
  20. Sun B, Hu L, Luo ZY, Chen XP, Zhou HH, Zhang W. DNA methylation perspectives in the pathogenesis of autoimmune diseases. Clin Immunol. 2016;164:21–7.View ArticlePubMedGoogle Scholar
  21. Jones PA, Liang G. Rethinking how DNA methylation patterns are maintained. Nat Rev Genet. 2009;10:805–11.View ArticlePubMedPubMed CentralGoogle Scholar
  22. Boyes J, Bird A. DNA methylation inhibits transcription indirectly via a methyl-CpG binding protein. Cell. 1991;64:1123–34.View ArticlePubMedGoogle Scholar
  23. Hu S, Wan J, Su Y, Song Q, Zeng Y, Nguyen HN, et al. DNA methylation presents distinct binding sites for human transcription factors. elife. 2013;2:e00726.PubMedPubMed CentralGoogle Scholar
  24. Zhang Z, Zhang R. Epigenetics in autoimmune diseases: pathogenesis and prospects for therapy. Autoimmun Rev. 2015;14:854–63.View ArticlePubMedGoogle Scholar
  25. Hedrich CM, Tsokos GC. Epigenetic mechanisms in systemic lupus erythematosus and other autoimmune diseases. Trends Mol Med. 2011;17:714–24.View ArticlePubMedPubMed CentralGoogle Scholar
  26. Wallace SL, Robinson H, Masi AT, Decker JL, Mccarty DJ, Yü TF. Preliminary criteria for the classification of the acute arthritis of primary gout. Arthritis Rheum. 1977;20:895–900.View ArticlePubMedGoogle Scholar
  27. Zhang W, Doherty M, Pascual E, Bardin T, Barskova V, Conaghan P, et al. EULAR evidence based recommendations for gout. Part I: diagnosis. Report of a task force of the standing committee for international clinical studies including therapeutics (ESCISIT). Ann Rheum Dis. 2006;65:1301–11.View ArticlePubMedPubMed CentralGoogle Scholar
  28. Hu S, Xie Z, Blackshaw S, Qian J, Zhu H. Characterization of protein-DNA interactions using protein microarrays. Cold Spring Harb Protoc. 2011;6:499–508.Google Scholar
  29. Guo JU, Su Y, Zhong C, Ming GL, Song H. Hydroxylation of 5-methylcytosine by TET1 promotes active DNA demethylation in the adult brain. Cell. 2011;145:423–34.View ArticlePubMedPubMed CentralGoogle Scholar
  30. Hu S, Xie Z, Onishi A, Yu X, Jiang L, Lin J, et al. Profiling the human protein-DNA interactome reveals ERK2 as a transcriptional repressor of interferon signaling. Cell. 2009;139:610–22.View ArticlePubMedPubMed CentralGoogle Scholar
  31. Vaissière T, Cuenin C, Paliwal A, Vineis P, Hainaut P, Herceg Z. Quantitative analysis of DNA methylation after whole bisulfitome amplification of a minute amount of DNA from body fluids. Epigenetics. 2009;4:221–30.View ArticlePubMedGoogle Scholar
  32. Wang J, Liu S, Wang B, Miao Z, Han L, Chu N, et al. Association between gout and polymorphisms in GCKR in male Han Chinese. Hum Genet. 2012;131:1261–5.View ArticlePubMedGoogle Scholar
  33. Thurman RE, Rynes E, Humbert R, Vierstra J, Maurano MT, Haugen E, et al. The accessible chromatin landscape of the human genome. Nature. 2012;489:75–82.View ArticlePubMedPubMed CentralGoogle Scholar
  34. Neph S, Vierstra J, Stergachis AB, Reynolds AP, Haugen E, Vernot B, et al. An expansive human regulatory lexicon encoded in transcription factor footprints. Nature. 2012;489:83–90.View ArticlePubMedPubMed CentralGoogle Scholar
  35. ENCODE Project Consortium. An integrated encyclopedia of DNA elements in the human genome. Nature. 2012;489:57–74.View ArticleGoogle Scholar
  36. Gerstein MB, Kundaje A, Hariharan M, Landt SG, Yan KK, Cheng C, et al. Architecture of the human regulatory network derived from ENCODE data. Nature. 2012;489:91–100.View ArticlePubMedPubMed CentralGoogle Scholar
  37. Kerr JS, Wilson CH. Nuclear receptor-binding protein 1: a novel tumour suppressor and pseudokinase. Biochem Soc Trans. 2013;41:1055–60.View ArticlePubMedGoogle Scholar
  38. Renauer PA, Saruhan-Direskeneli G, Coit P, Adler A, Aksu K, Keser G, et al. Identification of susceptibility loci in IL6, RPS9/LILRB3, and an intergenic locus on chromosome 21q22 in Takayasu arteritis in a genome-wide association study. Arthritis Rheumatol. 2015;67:1361–8.View ArticlePubMedPubMed CentralGoogle Scholar
  39. Wilson CH, Crombie C, van der Weyden L, Poulogiannis G, Rust AG, Pardo M, et al. Nuclear receptor binding protein 1 regulates intestinal progenitor cell homeostasis and tumour formation. EMBO J. 2012;31:2486–97.View ArticlePubMedPubMed CentralGoogle Scholar
  40. Wei H, Wang H, Ji Q, Sun J, Tao L, Zhou X. NRBP1 is downregulated in breast cancer and NRBP1 overexpression inhibits cancer cell proliferation through Wnt/β-catenin signaling pathway. Onco Targets Ther. 2015;8:3721–30.PubMedPubMed CentralGoogle Scholar
  41. Dharmapuri G, Doneti R, Philip GH, Kalle AM. Celecoxib sensitizes imatinib-resistant K562 cells to imatinib by inhibiting MRP1-5, ABCA2 and ABCG2 transporters via Wnt and Ras signaling pathways. Leuk Res. 2014;39:696–701.View ArticleGoogle Scholar
  42. Merriman TR. An update on the genetic architecture of hyperuricemia and gout. Arthritis Res Ther. 2015;17:98.View ArticlePubMedPubMed CentralGoogle Scholar
  43. Hannon E, Dempster E, Viana J, Burrage J, Smith AR, Macdonald R, et al. An integrated genetic-epigenetic analysis of schizophrenia: evidence for co-localization of genetic associations and differential DNA methylation. Genome Biol. 2016;17:176.View ArticlePubMedPubMed CentralGoogle Scholar
  44. Dick KJ, Nelson CP, Tsaprouni L, Sandling JK, Aïssi D, Wahl S, et al. DNA methylation and body-mass index: a genome-wide analysis. Lancet. 2014;383:1990–8.View ArticlePubMedGoogle Scholar
  45. Liu Y, Aryee MJ, Padyukov L, Fallin MD, Hesselberg E, Runarsson A, et al. Epigenome-wide association data implicate DNA methylation as an intermediary of genetic risk in rheumatoid arthritis. Nat Biotechnol. 2013;31:142–7.View ArticlePubMedPubMed CentralGoogle Scholar

Copyright

© The Author(s). 2017