Methylation of the genes ROD1, NLRC5, and HKR1 is associated with aging in Hainan centenarians

Background Human aging is a hot topic in biology, and it has been associated with DNA methylation changes at specific genomic sites. We aimed to study the changes of DNA methylation at a single-CpG-site resolution using peripheral blood samples from centenarians. Methods Using Illumina 450 K Methylation BeadChip microarray assays, we carried out a pool-based, epigenome-wide investigation of DNA methylation of blood samples from 12 centenarians and 12 healthy controls. Differentially methylated cytosine-phosphate-guanosine (CpG) sites were selected for further pyrosequencing analysis of blood samples from 30 centenarians and 30 healthy controls. Result We identified a total of 31 high-confidence CpG sites with differential methylation profiles between the groups: 9 (29%) were hypermethylated and 22 (71%) were hypomethylated in centenarians. It was also found that hypermethylation of HKR1 and hypomethylation of ROD1 and NLRC5 genes strongly correlated with age in centenarians. Conclusion Our results indicate that the methylation profile combination of HKR1, ROD1, and NLRC5 could be a promising biomarker for aging in Hainan centenarians. Electronic supplementary material The online version of this article (10.1186/s12920-018-0334-1) contains supplementary material, which is available to authorized users.


Background
Improvements in living conditions have meant that the average life expectancy of human beings is rising. Indeed, it has increased by approximately 10 years over the last century [1]. On the other hand, aging has attracted more and more attention. Aging is a very complex trait, involving environmental, genetic, and stochastic factors [2], and, in recent decades, has been a focus for scientists worldwide. Centenarian cohort is one of the most valuable models to study the mechanisms involved in human aging [3]. Centenarians are considered to have reached the extreme limits of the human life span, but still remain in relatively good health maintaining physiological functions and avoid succumbing to common fatal diseases [4,5]. Hainan, with a centenarian density of 21.46 per 100,000 people, is one of the most longevous regions in China. We managed to recruit 42 centenarians and 42 healthy adult controls from Hainan Province, and performed complete physical examinations on all subjects.
Human aging has largely been driven by extrinsic improvements to living conditions, including improved diet and reduced disease prevalence, as well as by both genetic and epigenetic factors [6][7][8]. Here, we evaluated the crucial role of epigenetic modification in gene regulation, with a particular focus on DNA methylation [9][10][11]. As early as 1975, it was proposed that DNA methylation could regulate gene activity. Methylation occurs primarily at CpG dinucleotides, which can form CpG islands containing an above average CpG density [12,13]. DNA methylation has been shown to regulate gene expression and plays a role in alternative silencing [14,15]. Furthermore, CpG islands overlap the transcription start sites (TSSs) of a majority of human genes, and TSS methylation is known to prevent transcriptional initiation [16,17]. Yet, the role of methylation in the gene body is less clear. Recent studies have demonstrated the presence of aging-associated CpG sites, which were either hypermethylated or hypomethylated [18][19][20][21]. Human genomic methylation has recently been advanced by utilization of the Illumina Human Methylation 450 K BeadChip Array (450 K array), which includes probes spanning 99% of human genes, evenly distributed across the genome. This system can detect the shores and shelves of CpG islands and non-CpG islands, as well as methylation sites in promoter regions, gene bodies, and untranslated regions (UTRs) [22][23][24].
In the present study, we aimed to identify DNA methylation sites by exploring differences in genomic methylation between centenarians and healthy controls. Using microarray technology and pyrosequencing and a pool-based strategy, we performed an epigenome-wide investigation of DNA methylation in a total of 42 centenarians from Hainan Province. We found that ROD1, NLRC5, and HKR1 genes were differentially methylated between centenarians and healthy controls, suggesting they might be associated with aging in Hainan centenarians.

Subjects
This study formed part of the "Health Investigation and Longevous Mechanism Study on Centenarians in Hainan Province" project, conducted between 2014 and 2016. The subject consisted of 42 centenarians (male: 15, female: 27; Han nationality) and 42 healthy adult controls (male: 18, female: 24; Han nationality). Characteristic information on centenarians and healthy controls is described in Table 1. Significant differences in characteristic information between centenarians and healthy controls are age-related. Healthy controls were excluded if they were suffering from known infectious diseases, or if they had a history of cancer, cardiovascular disease, or autoimmune disease. Neither the centenarians nor the healthy controls included in this study had any infections, and neither group had received any vaccinations in the 15 days prior to the collection of blood samples. Peripheral blood samples were collected by morning fasting. Additional informed consent was obtained from all individual participants for whom identifying information is included in this article.

Genomic DNA extraction
Peripheral blood samples were collected in EDTAtreated collection tubes by experienced nurses during a home visit, and then transferred to the Hainan branch of PLA General Hospital laboratory for processing within 4 h. Genomic DNA was extracted using a QIAamp DNA Mini kit (Qiagen, Germany), according to the manufacturer's instructions for the Spin Procedure. DNA was eluted in 60 μL of AE elution buffer, and then stored at − 20°C until use. The concentration and purity of DNA samples were assessed using a Qubit dsDNA HS Assay kit (Invitrogen, Eugene, OR, USA) according to the manufacturer's instructions, and 1% agarose gel electrophoresis.

Microarray analysis
Infinium Human Methylation 450BeadChips (Illumina) were used to assay genome-wide DNA methylation. The microarray chip contained 485,577 locusspecific CpG sites. Bisulfite conversion of DNA samples was done using the EZ DNA Methylation Kit (Zymo Research). Microarrays were normalized by methylumi Bioconductor [25][26][27]. Differentially methylated CpGs were selected using an algorithm in IMA Bioconductor. Probes were used in the analysis, after excluding those with a P value higher than 0.01 in more than 25% of the samples, those in chromosome Y and SNP sites. In this study, we assessed the meandifference β-value (Δβ) between the two sample groups for each CpG site. The Illumina platform shows a significant dye bias (probe-type design differences) in the two channels which will lead to bias in the estimates of Beta on the GoldenGate platform. Therefore, some normalization was required. Basically, it looks at the median intensities in the methylated and unmethylated channels (each measured in one color on the GoldenGate platform) at very low and very high beta values and sets these medians equal. And then a quantile approach was used to eliminate the batch effects among every microarray. In microarray data analysis, false positives are loci that are found to be statistically different between conditions. However, clinical samples may have relatively high individual heterogeneity, and the amount of samples is not enough. If we use p-values for multiple test corrections, it eliminated the false positive, but may have a certain false negative at the same time. Therefore, we considered a probe as differentially methylated if the absolute Δβ was higher than 0.2 and the statistical test was significant (p value < 0.01) with wilcox [28][29][30].

Bioinformatic analysis
As described above, genome-wide DNA methylation was assessed using a pool-based approach, in which the β-value for each probe represents the average methylation level of all samples in the pool. To identify differentially methylated probes (DMPs), it was assumed that the genome-wide methylation conformed to a Gaussian model [31], and then differences in methylation between the two data sets were calculated. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were conducted using the DAVID database. GO terms and KEGG information were downloaded from the Bioconductor site.

Bisulfite pyrosequencing analysis
Based on the results of the methylation microarrays, three CpG sites were selected for further pyrosequencing analysis, using blood samples from 30 centenarians and 30 healthy controls. Bisulfite modification of 1-2 μg DNA was performed using an EpiTect Bisulfite kit (Qiagen) according to the manufacturer's instructions. Next, bisulfite pyrosequencing primers were designed using PyroMark Assay Design software version 2.0 (Qiagen). Sequencing reactions and methylation level quantification were then performed using the PyroMark Q96 ID system and Fig. 1 Overview of the average genome-wide DNA methylation level of genomic loci at each chromosome in Hainan centenarians and healthy controls, analyzed by Circos software. The inner and outer tracks represent the average methylation level for centenarians and healthy controls, respectively. All chromosomes are presented as 10Mbp-wide windows. The average methylation level in each region is represented as the average β-value (0-1) for all of the probes in that region software (Qiagen), respectively. Pyrosequencing primer sequences are presented in additional file 1.

Statistical analyses
SPSS software version 17.0 (SPSS Inc., Chicago, IL, USA) was only used for the pyrosequencing analyses. Differences between two independent groups were compared using Student's t-tests and P < 0.01 was used to define statistical significance.

Results
Differential analysis of the DNA methylomes of Hainan centenarians and healthy controls Using microarray and pyrosequencing methods, we performed an analysis of DNA methylation patterns in An overview of the genome-wide methylation profiles for both groups is shown in Fig. 1. Different patterns of methylation between centenarians and healthy controls were observed at the majority of methylation loci, indicating that DNA methylation-based epigenetic profiles might be significant in aging. Having discovered these differences in DNA methylation, we then searched for individual DMPs. To do this, we performed a pooledtest to analyze the average methylation level at all of the CpG loci in the centenarian and healthy control groups.
In total, 31 DMPs were identified using the criteria P < 0.01 and |Δβ| > 0.2, as shown in Fig. 2a and b. In this analysis, Δβ > 0.2 and Δβ < − 0.2 refer to loci that were hypermethylated or hypomethylated in centenarians relative to healthy controls, respectively. Among the 31 DMPs identified, 9 were hypermethylated CpG sites and 22 were hypomethylated CpG sites (Table 2).

GO and KEGG functional enrichment analysis of differentially methylated genes
Our results showed that the 31 identified significantly differentially methylated sites described above could be mapped to 23 different genes. We performed GO term and KEGG pathway enrichment analyses to elucidate the biological functions of these genes, using the R package GOstats described previously [32]. The biological processes, molecular functions, and cellular components that were enriched in the GO analysis of longevity-associated DMP-containing genes (DMGs) are summarized in Table 3. In total, 3 biological processes, 6 cellular components, and 2 molecular function terms were enriched in the 23 DMGs. The enriched biological functions, such as antigen processing and presentation of peptide antigen (GO: 0048002) are associated with immune system function. Furthermore, enriched cellular components include major histocompatibility complex (MHC) protein complex and insoluble fraction (GO: 0042611 and GO: 0005626, respectively). MHC receptor activity, class I and II (GO: 0032393 and GO: 0032395, respectively), are enriched molecular functions in the DMG subset. KEGG analysis revealed the enrichment of only several pathways in the DMG subset (Table 4), with the most significantly enriched biological pathways being graft-versus-host disease (KEGG: 05332) and allograft rejection (KEGG: 05330).

Changes in ROD1, NLRC5 and HKR1 DNA methylation in centenarians
The CpGs associated with aging were not uniformly distributed across chromosomes, CpG islands, or genes; significantly differentially methylated sites were overrepresented on chromosomes 3, 6, 9, and 19 ( Table 2).
The gene locations of the DMPs were assessed, and defined as follows: 1500 bp upstream of the TSS (TSS1500), 200 bp upstream of the TSS (TSS200), the 5'-UTR, exon 1, the gene body (encompassing all exons except exon 1), and the 3'-UTR. Some of the enriched DMPs were located within the gene body, but the role of methylation within the gene body is less clear. DNA methylation of promoter regions is a repressive epigenetic mark that regulates gene expression. The three most significantly hypermethylated genes in the centenarian group were ELOVL fatty acid elongase 2 (ELOVL2; ID cg24724428), histidine kinase rhodopsins 1 (HKR1; ID cg23448505; ID cg14166009; ID Table 3 Analysis of significantly hypermethylated differentially methylated probes (DMPs) GO  The CpG islands in the ELOVL2 gene was located within the gene promoter (TSS1500), HKR1 gene was located within TSS1500, TSS200 and 5'-UTR, while RPTOR genes were located within the gene body. Previous studies demonstrated that methylation of a CpG island within the ELOVL2 gene strongly correlates with age; indeed such methylation was declared a novel epigenetic marker for aging [33,34]. Furthermore, seven gene sites (ID cg23448505; ID cg14166009; ID cg13687570; ID cg26734888; ID cg12948621; ID cg24834889; IDcg05280698) belong to the same gene(HKR1), therefore, we chose one of HKR (ID cg26734888) to pyrosequence in a large sample set. The three most significantly hypomethylated genes in the centenarians were NLR family CARD domain containing 5 (NLRC5; ID cg07839457), regulator of differentiation 1 (ROD1; ID cg00008629) and LOC153328 (ID cg04233669). The CpG islands in the NLRC5 gene was located TSS1500; ROD1 gene was located 5'UTR, while LOC153328 gene was located within the gene body. It has been reported that the methylation level of CpGs within the coding region can regulate gene transcription [35]. To confirm these observations, we chose NLRC5 (ID cg07839457) and ROD1 (ID cg00008629) genes to pyrosequence in a large sample set. In total, 60 peripheral blood samples (30 samples from each group) underwent pyrosequencing. As expected, significant differences in the level of methylation at the sites cg00008629 (ROD1), cg07839457 (NLRC5), and cg26734888 (HKR1) were observed between Hainan centenarians and healthy controls ( Table 5 and Fig. 3). Using a regression-based approach, we further analyzed the DNA methylation data based on cell type compositions [36]. Our further analyses indicated that the differences in the level of methylation of cg00008629 (ROD1), cg07839457 (NLRC5), and cg26734888 (HKR1) sites were significant (p value < 0.05) between centenarians and healthy controls ( Table 6).

Discussion
In our current study, we analyzed theDNA methylation of blood samples collected from centenarians and healthy controls in Hainan province, China. In summary, we found 31 DMPs using a 450 K array. Furthermore, our pyrosequencing analysis showed that methylation level at three of CpGs, ROD1, NLRC5, and HKR1 strongly correlates with aging in centenarians.
Our results also indicate that 31 high-confidence CpG sites with differential methylation profiles between the groups; 9 (29%) were hypermethylated and 22 (71%) were hypomethylated in centenarians. KEGG pathway analysis showed that the most significantly  . 3 The level of methylation at the indicated CpG sites in centenarians (blue boxes) and healthy controls (green boxes). The boxes represent the 25th to the 75th centiles, and solid lines represent the median value for each sample. Error bars represent 95% confidence intervals. * p < 0.05 indicated the significant different between healthy controls and centenarians enriched biological pathways being graft-versus-host disease and allograft rejection. Further studies are needed to study how these pathways relate to the aging exhibited by centenarians. The ROD1 gene was discovered to be hypomethylated in centenarians. This gene encodes an RNAbinding protein that binds preferentially to poly(G) and poly(U) sequences in vitro, and was initially regarded as an inhibitor of differentiation [37]. Several ROD1 transcriptional variants, encoding different isoforms, have been identified. In addition, ROD1 was recently found to be a member of the heterogeneous nuclear ribonucleoprotein family, which regulate the alternative splicing of pre-RNA, and can thus modulate gene expression by repressing the expression of exons in a tissue specific manner [38,39]. Further, ROD1 binding to C-terminal peptide was strongly inhibited by synthetic RNA and weakly inhibited by a synthetic phosphorylated peptide designed to mimic the C-terminal domain of RNA polymerase II. Interestingly, ROD1 expression and knockdown experiments suggest a role in the regulation of C-terminal peptide mitogenic activity [40], suggesting a potential correlation between hypomethylation of ROD1 gene and aging in centenarians.
NLRC5, a member of the NOD-like receptor gene family and widely expressed in tumor tissues, is a key transcriptional regulator of MHC class I molecules [41,42]. Cancers can escape immune surveillance by reducing the expression of MHC class I molecules. NLRC5, an important regulator of innate immune responses, may affect patient survival by regulating this immune evasion via MHC class I molecules. NLRC5 could therefore be exploited to restore tumor immunogenicity and to stimulate protective antitumor immunity [43]. While the role of NLRC5 in both innate and adaptive immune signaling is well established, its function in inflammation remains poorly understood. The regulation of NF-κB activation is crucial for many biological functions, and NF-κB dysregulation is associated with immune deficiency, infectious disease, inflammation, and cancer [44,45]. NLRC5 has been reported to inhibit NF-κB activation by interacting with IKK-α/−β and blocking its phosphorylation [46]. To date, there have been no reports which implicate NLRC5 in aging, however, we found that NLRC5 was hypomethylated in centenarians. Methylation of the NLRC5 is negatively correlated with NLRC5 expression in many cells [47,48]. It is possible that this NLRC5 hypomethylation is involved in aging through the regulation of various biological pathways.
In contrast to ROD1 and NLRC5, HKR1 was hypermethylated in centenarians. HKR1 encodes Krueppelrelated zinc finger protein 1, a member of the GLI-Kruppel zinc finger family that is important for gene regulation [49]. The HKR1 protein possesses a serine/threonine-rich domain, and a Ca 2+ -binding consensus sequence within its cytoplasmic domain. Recently, HKR1 mRNA expressions were found to be higher in lung cancer tissues than in normal lung tissues, and HKR1 may be involved in the regulation of a signaling pathway involved in lung cancer progression [50]. Although the physiological roles of HKR1 remain unclear, HKR1 knockout is lethal, indicating that this gene is essential for survival. Further research is needed to fully understand the significance of HKR1 hypermethylation in centenarians and its effect on aging.
However, it is worth mentioning that peripheral blood is a heterogeneous collection of different cell types, each with a very different DNA methylation profile. Differences in the relative proportions of these components may contribute to differences in DNA methylation [51][52][53]. To identify the potential effect of cellular proportions in the DNA methylation differentiation, we used a regression-based approach to adjusting for cell type composition in DNA methylation data. Our results indicated that the methylation cg00008629 (ROD1), cg07839457 (NLRC5), and cg26734888 (HKR1) were also found significant differences between centenarians and healthy controls.
To our knowledge, the association between the genes ROD1, NLRC5, and HKR1 and aging have not been previously reported. We here found significantly different levels of DNA methylation at the ROD1, NLRC5, and HKR1 genes in centenarians and healthy controls. Our findings implicate that these genes may play important roles in aging via the regulation of various biological pathways and could be used as biomarker genes for aging.

Conclusions
In summary, we present a comprehensive comparison of genomic methylation in centenarians and healthy controls residing in Hainan Province. We report an epigenetic signature consisting of three genes, which can be used to estimate aging from blood samples. Changes to the methylation of CpG islands, namely, hypomethylation within the ROD1 and NLRC5 and hypermethylation within HKR1 genes, are possibly correlated with aging in centenarians. While further studies are required to understand the mechanisms and signaling pathways through which ROD1, NLRC5, and HKR1 methylation affect aging, the methylation of these genes constitutes a signature that could be used to evaluate aging. The combination of ROD1, NLRC5, and HKR1 gene methylation is therefore a potential novel biomarker for aging.

Additional file
Additional file 1: "Pyrosequencing primer sequences of three sites". This file contains detailed description of pyrosequencing primer sequences.

Acknowledgements
We appreciate all the participants of the Hainan Centenarian Study for their continued cooperation. We thank Hanqi Yin (Guangdong Longsee biomedical corporation) for guidance on data analysis.

Funding
This study was funded by Key research funding of Hainan (grant nos. 16A200143 and 14A110069), and Key Research and Development Project of Hainan (grant nos. ZDYF2016135 and ZDYF2017095). There is no role of the funding body in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Availability of data and materials
The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.
Authors' contributions QZ designed the study and drafted the manuscript. XPC and QZH conducted experiments. CXN performed statistical analyses. YY prepared the samples. YLZ and FXL conceived the study, participated in its design and coordination, and helped to draft the manuscript. The authors declare that they have no conflict of interest. All authors have read and approved the final manuscript.
Ethics approval and consent to participate Prior to investigation and physical examination, we have got the consent from the participants or from their family members such as sons and daughter-in-laws. If the participants could not sign their names because of illiteracy or being cognitive-impaired, their names were signed by their direct relatives after the consent of the participants. The Ethics Committee of the Hainan branch of the Chinese People's Liberation Army General Hospital (Sanya, Hainan) approved the study protocol (No. of serial: 301hn11201601).

Consent for publication
Not applicable.

Competing interests
On behalf of all authors, the corresponding author states that there is no conflict of interest.