Genome-wide DNA methylation profiling in rheumatoid arthritis identifies disease-associated methylation changes that are distinct to individual T- and B-lymphocyte populations

Changes to the DNA methylome have been described in patients with rheumatoid arthritis (RA). In previous work, we reported genome-wide methylation differences in T-lymphocyte and B-lymphocyte populations from healthy individuals. Now, using HumanMethylation450 BeadChips to interrogate genome-wide DNA methylation, we have determined disease-associated methylation changes in blood-derived T- and B-lymphocyte populations from 12 female patients with seropositive established RA, relative to 12 matched healthy individuals. Array data were analyzed using NIMBL software and bisulfite pyrosequencing was used to validate array candidates. Genome-wide DNA methylation, determined by analysis of LINE-1 sequences, revealed higher methylation in B-lymphocytes compared with T-lymphocytes (P ≤ 0.01), which is consistent with our findings in healthy individuals. Moreover, loci-specific methylation differences that distinguished T-lymphocytes from B-lymphocytes in healthy individuals were also apparent in RA patients. However, disease-associated methylation differences were also identified in RA. In these cases, we identified 509 and 252 CpGs in RA-derived T- and B-lymphocytes, respectively, that showed significant changes in methylation compared with their cognate healthy counterparts. Moreover, this included a restricted set of 32 CpGs in T-lymphocytes and 20 CpGs in B-lymphocytes (representing 15 and 10 genes, respectively, and including two, MGMT and CCS, that were common to both cell types) that displayed more substantial changes in methylation. These changes, apparent as hyper- or hypo-methylation, were independently confirmed by pyrosequencing analysis. Validation by pyrosequencing also revealed additional sites in some candidate genes that also displayed altered methylation in RA. In this first study of genome-wide DNA methylation in individual T- and B-lymphocyte populations in RA patients, we report disease-associated methylation changes that are distinct to each cell type and which support a role for discrete epigenetic regulation in this disease.     

Early- and late-onset preeclampsia and the DNA methylation of circadian clock and clock-controlled genes in placental and newborn tissues

The placenta is important in providing a healthy environment for the fetus and plays a central role in the pathophysiology of preeclampsia (PE). Fetal and placental developments are influenced by epigenetic programming. There is some evidence that PE is controlled to an altered circadian homeostasis. In a nested case–control study embedded in the Rotterdam Periconceptional Cohort, we obtained placental tissue, umbilical cord leukocytes (UCL), and human umbilical venous endothelial cells of 13 early-onset PE, 16 late-onset PE and 83 controls comprising 36 uncomplicated and 47 complicated pregnancies, i.e. 27 fetal growth restricted and 20 spontaneous preterm birth. To investigate the associations between PE and the epigenetics of circadian clock and clock-controlled genes in placental and newborn tissues, genome-wide DNA methylation analysis was performed using the Illumina HumanMethylation450K BeadChip and a candidate-gene approach using ANCOVA was applied on 939 CpGs of 39 circadian clock and clock-controlled genes. DNA methylation significantly differed in early-onset PE compared with spontaneous preterm birth at 6 CpGs in placental tissue (3.73^E-5 ≤ p ≤ 0.016) and at 21 CpGs in UCL (1.09^E-5≤ p ≤ 0.024). In early-onset PE compared with fetal growth restriction 2 CpGs in placental tissue (p < 0.05) and 8 CpGs in uncomplicated controls (4.78^E-5≤ p ≤ 0.049) were significantly different. Moreover, significantly different DNA methylation in early-onset PE compared with uncomplicated controls was shown at 6 CpGs in placental tissue (1.36^E-4≤ p ≤ 0.045) and 11 CpGs in uncomplicated controls (2.52^E-6≤ p ≤ 0.009). No significant associations were shown with late-onset PE between study groups or tissues. The most differentially methylated CpGs showed hypomethylation in placental tissue and hypermethylation in uncomplicated controls. In conclusion, DNA methylation of circadian clock and clock-controlled genes demonstrated most differences in UCL of early-onset PE compared with spontaneous preterm birth. Implications of the tissue-specific variations in epigenetic programming for circadian performance and long-term health need further investigation.     

DNA methylation pattern of CALCA in preterm neonates with bacterial sepsis as a putative epigenetic biomarker

Diagnosis of bacterial sepsis in preterm neonates can be difficult when using serum markers that rely on physiological changes because these changes may not necessarily be the result of bacterial infections alone. This retrospective investigation explores the potential use of the DNA methylation pattern of CpG sites in the promoter region of the calcitonin-related polypeptide α (CALCA) gene as an epigenetic biomarker for bacterial sepsis in preterm newborns. Four novel changes in the DNA methylation of eight CpG sites were detected in this gene and are present only in neonates with bacterial sepsis: (1) partial methylation at -769 CpG in gram-negative or gram-positive early onset sepsis (EOS) and late onset sepsis (LOS) episodes; (2) demethylation of 8 CpGs in gram-negative EOS followed by LOS (ELS) and in gram-negative EOS; (3) demethylation of 7 CpGs in gram-positive ELS and gram-positive EOS; (4) -771 C:G > T:A; 5′ de novo -778 CpG mutation on both alleles in EOS. These changes were not detected in birth weight and gestational age matched controls or in newborns with isolated infections. Our results indicate that the DNA methylation pattern of the promoter region of the CALCA gene varies in different types of bacterial preterm sepsis, thus suggesting a potential use as an epigenetic biomarker. A prospective confirmation of these results is essential.     

CLOCK,PER2 and BMAL1 DNA Methylation: Association with Obesity and Metabolic Syndrome Characteristics and Monounsaturated Fat Intake

The circadian clock system instructs 24-h rhythmicity on gene expression in essentially all cells, including adipocytes, and epigenetic mechanisms may participate in this regulation. The aim of this research was to investigate the influence of obesity and metabolic syndrome (MetS) features in clock gene methylation and the involvement of these epigenetic modifications in the outcomes. Sixty normal-weight, overweight and obese women followed a 16-weeks weight reduction program. DNA methylation levels at different CpG sites of CLOCK, BMAL1 and PER2 genes were analyzed by Sequenom's MassARRAY in white blood cells obtained before the treatment. Statistical differences between normal-weight and overweight?+?obese subjects were found in the methylation status of different CpG sites of CLOCK (CpGs 1, 5-6, 8 and 11-14) and, with lower statistical significance, in BMAL1 (CpGs 6-7, 8, 15 and 16-17). The methylation pattern of different CpG sites of the three genes showed significant associations with anthropometric parameters such as body mass index and adiposity, and with a MetS score. Moreover, the baseline methylation levels of CLOCK CpG 1 and PER2 CpGs 2-3 and 25 correlated with the magnitude of weight loss. Interestingly, the percentage of methylation of CLOCK CpGs 1 and 8 showed associations with the intake of monounsaturated and polyunsaturated fatty acids. This study demonstrates for the first time an association between methylation status of CpG sites located in clock genes (CLOCK, BMAL1 and PER2) with obesity, MetS and weight loss. Moreover, the methylation status of different CpG sites in CLOCK and PER2 could be used as biomarkers of weight-loss success, particularly CLOCK CPGs 5-6. (Author correspondence: garaulet@um.es)     

DNA methylation analysis reveals the effect of arsenic on gestational diabetes mellitus

BackgroundArsenic (As) exposure is one of the risk factors for gestational diabetes mellitus (GDM). This study aimed to explore the effect of As-exposure on DNA methylation in GDM and to establish a risk assessment model of GDM in As exposed pregnant women.MethodWe collected elbow vein blood of pregnant women before delivery to measure As concentration and DNA methylation data. Then compared the DNA methylation data and established a nomogram.ResultWe identified a total of 10 key differentially methylated CpGs (DMCs) and found 6 corresponding genes. Functions were enriched in Hippo signaling pathway, cell tight junction, prophetic acid metabolism, ketone body metabolic process, and antigen processing and presentation. A nomogram was established that can predict GDM risks (c-index = 0.595, s:p = 0.973).ConclusionWe found 6 genes associated with GDM with high As exposure. The prediction of the nomograms has been proven to be effective.     

Maternal cigarette smoking during pregnancy and offspring DNA methylation in midlife

Maternal smoking in pregnancy (MSP) has been associated with DNA methylation in specific CpG sites (CpGs) in infants and children. We investigated whether MSP, independent of own personal active smoking, was associated with midlife DNA methylation in CpGs that were previously identified in studies of MSP-DNA methylation in children. We used data on MSP collected from pregnant mothers of 89 adult women born in 1959–1964 and measured DNA methylation in blood (granulocytes) collected in 2001–2007 (mean age: 43 years). Seventeen CpGs were differentially methylated by MSP, with multiple CpGs mapping to CYP1A1, MYO1G, AHRR, and GFI1. These associations were consistent in direction with prior studies (e.g., MSP associated with more and less methylation in AHRR and CYP1A1, respectively) and, with the exception of AHRR CpGs, were not substantially altered by adjustment for active smoking. These preliminary results confirm prior prospective reports that MSP influences the offspring DNA methylation, and extends the timeframe to midlife, and suggest that these effects may persist into adulthood, independently of active smoking.     

Tobacco exposure-related alterations in DNA methylation and gene expression in human monocytes: the Multi-Ethnic Study of Atherosclerosis (MESA)

Alterations in DNA methylation and gene expression in blood leukocytes are potential biomarkers of harm and mediators of the deleterious effects of tobacco exposure. However, methodological issues, including the use of self-reported smoking status and mixed cell types have made previously identified alterations in DNA methylation and gene expression difficult to interpret. In this study, we examined associations of tobacco exposure with DNA methylation and gene expression, utilizing a biomarker of tobacco exposure (urine cotinine) and CD14+ purified monocyte samples from 934 participants of the community-based Multi-Ethnic Study of Atherosclerosis (MESA). Urine cotinine levels were measured using an immunoassay. DNA methylation and gene expression were measured with microarrays. Multivariate linear regression was used to test for associations adjusting for age, sex, race/ethnicity, education, and study site. Urine cotinine levels were associated with methylation of 176 CpGs [false discovery rate (FDR)<0.01]. Four CpGs not previously identified by studies of non-purified blood samples nominally replicated (P value<0.05) with plasma cotinine-associated methylation in 128 independent monocyte samples. Urine cotinine levels associated with expression of 12 genes (FDR<0.01), including increased expression of P2RY6 (Beta ± standard error = 0.078 ± 0.008, P = 1.99 × 10^?22), a gene previously identified to be involved in the release of pro-inflammatory cytokines. No cotinine-associated (FDR<0.01) methylation profiles significantly (FDR<0.01) correlated with cotinine-associated (FDR<0.01) gene expression profiles. In conclusion, our findings i) identify potential monocyte-specific smoking-associated methylation patterns and ii) suggest that alterations in methylation may not be a main mechanism regulating gene expression in monocytes in response to cigarette smoking.     

Genome-wide DNA methylation profiling in rheumatoid arthritis identifies disease-associated methylation changes that are distinct to individual T- and B-lymphocyte populations

Changes to the DNA methylome have been described in patients with rheumatoid arthritis (RA). In previous work, we reported genome-wide methylation differences in T-lymphocyte and B-lymphocyte populations from healthy individuals. Now, using HumanMethylation450 BeadChips to interrogate genome-wide DNA methylation, we have determined disease-associated methylation changes in blood-derived T- and B-lymphocyte populations from 12 female patients with seropositive established RA, relative to 12 matched healthy individuals. Array data were analyzed using NIMBL software and bisulfite pyrosequencing was used to validate array candidates. Genome-wide DNA methylation, determined by analysis of LINE-1 sequences, revealed higher methylation in B-lymphocytes compared with T-lymphocytes (P ≤ 0.01), which is consistent with our findings in healthy individuals. Moreover, loci-specific methylation differences that distinguished T-lymphocytes from B-lymphocytes in healthy individuals were also apparent in RA patients. However, disease-associated methylation differences were also identified in RA. In these cases, we identified 509 and 252 CpGs in RA-derived T- and B-lymphocytes, respectively, that showed significant changes in methylation compared with their cognate healthy counterparts. Moreover, this included a restricted set of 32 CpGs in T-lymphocytes and 20 CpGs in B-lymphocytes (representing 15 and 10 genes, respectively, and including two, MGMT and CCS, that were common to both cell types) that displayed more substantial changes in methylation. These changes, apparent as hyper- or hypo-methylation, were independently confirmed by pyrosequencing analysis. Validation by pyrosequencing also revealed additional sites in some candidate genes that also displayed altered methylation in RA. In this first study of genome-wide DNA methylation in individual T- and B-lymphocyte populations in RA patients, we report disease-associated methylation changes that are distinct to each cell type and which support a role for discrete epigenetic regulation in this disease.     

Deep Bisulfite Sequencing of Aberrantly Methylated Loci in a Patient with Multiple Methylation Defects

NLRP7 is a maternal effect gene as maternal mutations in this gene cause recurrent hydatidiform moles, spontaneous abortions and stillbirths, whereas live births are very rare. We have studied a patient with multiple anomalies born to a mother with a heterozygous NLRP7 mutation. By array-based CpG methylation analysis of blood DNA from the patient, his parents and 18 normal controls on Illumina Infinium HumanMethylation27 BeadChips we found that the patient had methylation changes (delta ß ≥ 0.3) at many imprinted loci as well as at 87 CpGs associated with 85 genes of unknown imprinting status. Using a pseudoproband (permutation) approach, we found methylation changes at only 7-24 CpGs (mean 15; standard deviation 4.84) in the controls. Thus, the number of abberantly methylated CpGs in the patient is more than 14 standard deviations higher. In order to identify novel imprinted genes among the 85 conspicuous genes in the patient, we selected 19 (mainly hypomethylated) genes for deep bisulfite amplicon sequencing on the ROCHE/454 Genome Sequencer in the patient and at least two additional controls. These controls had not been included in the array analysis and were heterozygous for a single nucleotide polymorphism at the test locus, so that allele-specific DNA methylation patterns could be determined. Apart from FAM50B, which we proved to be imprinted in blood, we did not observe allele-specific DNA methylation at the other 18 loci. We conclude that the patient does not only have methylation defects at imprinted loci but (at least in blood) also an excess of methylation changes at apparently non-imprinted loci.     

Key epigenetic changes associated with lung cancer development

Epigenetic alterations are a common event in lung cancer and their identification can serve to inform on the carcinogenic process and provide clinically relevant biomarkers. Using paired tumor and non-tumor lung tissues from 146 individuals from three independent populations we sought to identify common changes in DNA methylation associated with the development of non-small cell lung cancer. Pathologically normal lung tissue taken at the time of cancer resection was matched to tumorous lung tissue and together were probed for methylation using Illumina GoldenGate arrays in the discovery set (n = 47 pairs) followed by bisulfite pyrosequencing for validation sets (n = 99 pairs). For each matched pair the change in methylation at each CpG was calculated (the odds ratio), and these ratios were averaged across individuals and ranked by magnitude to identify the CpGs with the greatest change in methylation associated with tumor development. We identified the top gene-loci representing an increase in methylation (HOXA9, 10.3-fold and SOX1, 5.9-fold) and decrease in methylation (DDR1, 8.1-fold). In replication testing sets, methylation was higher in tumors for HOXA9 (p < 2.2 × 10^?16) and SOX1 (p < 2.2 × 10^?16) and lower for DDR1 (p < 2.2 × 10^?16). The magnitude and strength of these changes were consistent across squamous cell and adenocarcinoma tumors. Our data indicate that the identified genes consistently have altered methylation in lung tumors. Our identified genes should be included in translational studies that aim to develop screening for early disease detection.     

Age-related epigenome-wide DNA methylation and hydroxymethylation in longitudinal mouse blood

DNA methylation at cytosine-phosphate-guanine (CpG) dinucleotides changes as a function of age in humans and animal models, a process that may contribute to chronic disease development. Recent studies have investigated the role of an oxidized form of DNA methylation – 5-hydroxymethylcytosine (5hmC) – in the epigenome, but its contribution to age-related DNA methylation remains unclear. We tested the hypothesis that 5hmC changes with age, but in a direction opposite to 5-methylcytosine (5mC), potentially playing a distinct role in aging. To characterize epigenetic aging, genome-wide 5mC and 5hmC were measured in longitudinal blood samples (2, 4, and 10 months of age) from isogenic mice using two sequencing methods – enhanced reduced representation bisulfite sequencing and hydroxymethylated DNA immunoprecipitation sequencing. Examining the epigenome by age, we identified 28,196 unique differentially methylated CpGs (DMCs) and 8,613 differentially hydroxymethylated regions (DHMRs). Mouse blood showed a general pattern of epigenome-wide hypermethylation and hypo-hydroxymethylation with age. Comparing age-related DMCs and DHMRs, 1,854 annotated genes showed both differential 5mC and 5hmC, including one gene – Nfic – at five CpGs in the same 250 bp chromosomal region. At this region, 5mC and 5hmC levels both decreased with age. Reflecting these age-related epigenetic changes, Nfic RNA expression in blood decreased with age, suggesting that age-related regulation of this gene may be driven by 5hmC, not canonical DNA methylation. Combined, our genome-wide results show age-related differential 5mC and 5hmC, as well as some evidence that changes in 5hmC may drive age-related DNA methylation and gene expression.     

Genome-Wide DNA Methylation Analysis of Systemic Lupus Erythematosus Reveals Persistent Hypomethylation of Interferon Genes and Compositional Changes to CD4+ T-cell Populations

Systemic lupus erythematosus (SLE) is an autoimmune disease with known genetic, epigenetic, and environmental risk factors. To assess the role of DNA methylation in SLE, we collected CD4+ T-cells, CD19+ B-cells, and CD14+ monocytes from 49 SLE patients and 58 controls, and performed genome-wide DNA methylation analysis with Illumina Methylation450 microarrays. We identified 166 CpGs in B-cells, 97 CpGs in monocytes, and 1,033 CpGs in T-cells with highly significant changes in DNA methylation levels (p<1×10⁻⁸) among SLE patients. Common to all three cell-types were widespread and severe hypomethylation events near genes involved in interferon signaling (type I). These interferon-related changes were apparent in patients collected during active and quiescent stages of the disease, suggesting that epigenetically-mediated hypersensitivity to interferon persists beyond acute stages of the disease and is independent of circulating interferon levels. This interferon hypersensitivity was apparent in memory, naïve and regulatory T-cells, suggesting that this epigenetic state in lupus patients is established in progenitor cell populations. We also identified a widespread, but lower amplitude shift in methylation in CD4+ T-cells (>16,000 CpGs at FDR<1%) near genes involved in cell division and MAPK signaling. These cell type-specific effects are consistent with disease-specific changes in the composition of the CD4+ population and suggest that shifts in the proportion of CD4+ subtypes can be monitored at CpGs with subtype-specific DNA methylation patterns.     

Regions of variable DNA methylation in human placenta associated with newborn neurobehavior

The placenta regulates the in utero environment and functionally impacts fetal development. Candidate gene studies identified variation in placental DNA methylation is associated with newborn neurologic and behavioral outcomes including movement quality, lethargic behavior, attention, and arousal. We sought to identify novel regions of variable DNA methylation associated with newborn attention, lethargy, quality of movement, and arousal by performing an epigenome-wide association study in 335 infants from a US birth cohort. Methylation status was quantified using the Illumina HumanMethylation450 BeadChip array and associations to newborn outcomes assessed by the NICU Network Neurobehavioral Scales (NNNS) were identified while incorporating established bioinformatics algorithms to control for confounding by cell type composition. Methylation of CpGs within FHIT (cg15970800) and ANKRD11 (cg16710656) demonstrated genome-wide significance (P < 1.8 × 10^?7) in specific associations with infant attention. CpGs whose differential methylation was associated with all 4 neurobehavioral outcomes were common to 50 genes involved in biological processes relating to cellular adhesion and nervous system development. Comprehensive methylation profiling identified relationships between methylation of FHIT and ANKRD11, which have been previously linked to neurodevelopment and behavioral outcomes in genetic association studies. Subtle changes in DNA methylation of these genes within the placenta may impact normal variation of a newborn's ability to alter and track visual and auditory stimuli. Gene ontology analysis suggested that those genes with variable methylation related to these outcomes are over-represented in biological pathways involved in brain development and placental physiology, supportive of our hypothesis for a key role of the placenta in neurobehavioral outcomes.     

Cumulative lifetime maternal stress and epigenome-wide placental DNA methylation in the PRISM cohort

Evolving evidence links maternal stress exposure to changes in placental DNA methylation of specific genes regulating placental function that may have implications for the programming of a host of chronic disorders. Few studies have implemented an epigenome-wide approach. Using the Infinium HumanMethylation450 BeadChip (450K), we investigated epigenome-wide placental DNA methylation in relation to maternal experiences of traumatic and non-traumatic stressors over her lifetime assessed using the Life Stressor Checklist-Revised (LSC-R) survey (n = 207). We found differential DNA methylation at epigenome-wide statistical significance (FDR = 0.05) for 112 CpGs. Additionally, we observed three clusters that exhibited differential methylation in response to high maternal lifetime stress. Enrichment analyses, conducted at an FDR = 0.20, revealed lysine degradation to be the most significant pathway associated with maternal lifetimes stress exposure. Targeted enrichment analyses of the three largest clusters of probes, identified using the gap statistic, were enriched for genes associated with endocytosis (i.e., SMAP1, ANKFY1), tight junctions (i.e., EPB41L4B), and metabolic pathways (i.e., INPP5E, EEF1B2). These pathways, also identified in the top 10 KEGG pathways associated with maternal lifetime stress exposure, play important roles in multiple physiological functions necessary for proper fetal development. Further, two genes were identified to exhibit multiple probes associated with maternal lifetime stress (i.e., ANKFY1, TM6SF1). The methylation status of the probes belonging to each cluster and/or genes exhibiting multiple hits, may play a role in the pathogenesis of adverse health outcomes in children born to mothers with increased lifetime stress exposure.     

Epigenome-wide differential DNA methylation between HIV-infected and uninfected individuals

Epigenetic control of human immunodeficiency virus-1 (HIV-1) genes is critical for viral integration and latency. However, epigenetic changes in the HIV-1-infected host genome have not been well characterized. Here, we report the first large-scale epigenome-wide association study of DNA methylation for HIV-1 infection. We recruited HIV-infected (n = 261) and uninfected (n = 117) patients from the Veteran Aging Cohort Study (VACS) and all samples were profiled for 485,521 CpG sites in DNA extracted from the blood. After adjusting for cell type and clinical confounders, we identified 20 epigenome-wide significant CpGs for HIV-1 infection. Importantly, 2 CpGs in the promoter of the NLR family, CARD domain containing gene 5 (NLRC5), a key regulator of major histocompatibility complex class I gene expression, showed significantly lower methylation in HIV-infected subjects than in uninfected subjects (cg07839457: t = ?6.03, P_nominal = 4.96 × 10^?9; cg16411857: t = ?7.63, P_nominal = 3.07 × 10^?13). Hypomethylation of these 2 CpGs was replicated in an independent sample (GSE67705: cg07839457: t = ?4.44, P_nominal = 1.61 × 10^?5; cg16411857: t = ?5.90; P = 1.99 × 10^?8). Methylation of these 2 CpGs in NLRC5 was negatively correlated with viral load in the 2 HIV-infected samples (cg07839457: P = 1.8 × 10^?4; cg16411857: P = 0.03 in the VACS; and cg07839457: P = 0.04; cg164111857: P = 0.01 in GSE53840). Our findings demonstrate that differential DNA methylation is associated with HIV infection and suggest the involvement of a novel host gene, NLRC5, in HIV pathogenesis.     

Global DNA methylation profiling in peripheral blood cells of South African women with gestational diabetes mellitus

Background/Objective: Recently, several studies have reported that DNA methylation changes in tissue are reflected in blood, sparking interest in the potential use of global DNA methylation as a biomarker for gestational diabetes mellitus (GDM). This study investigated whether global DNA methylation is associated with GDM in South African women.

Methods: Global DNA methylation was quantified in peripheral blood cells of women with (n?=?63) or without (n?=?138) GDM using the MDQ1 Imprint® DNA Quantification Kit.

Results: Global DNA methylation levels were not different between women with or without GDM and were not associated with fasting glucose or insulin concentrations. However, levels were 18% (p?=?0.012) higher in obese compared to non-obese pregnant women and inversely correlated with serum adiponectin concentrations (p?=?0.005).

Discussion: Contrary to our hypothesis, global DNA methylation was not associated with GDM in our population. These preliminary findings suggest that despite being a robust marker of overall genomic methylation that offers opportunities as a biomarker, global DNA methylation profiling may not offer the resolution required to detect methylation differences in the peripheral blood cells of women with GDM. Moreover, global DNA methylation in peripheral blood cells may not reflect changes in placental tissue. Further studies in a larger sample are required to explore the candidacy of a more targeted approach using gene-specific methylation as a biomarker for GDM in our population.     

Epigenome-wide profiling identifies significant differences in DNA methylation between matched-pairs of T- and B-lymphocytes from healthy individuals

Multiple reports now describe changes to the DNA methylome in rheumatoid arthritis and in many cases have analyzed methylation in mixed cell populations from whole blood. However, these approaches may preclude the identification of cell type-specific methylation, which may subsequently bias identification of disease-specific changes. To address this possibility, we conducted genome-wide DNA methylation profiling using HumanMethylation450 BeadChips to identify differences within matched pairs of T-lymphocytes and B-lymphocytes isolated from the peripheral blood of 10 healthy females. Array data were processed and differential methylation identified using NIMBL software. Validation of array data was performed by bisulfite pyrosequencing. Genome-wide DNA methylation was initially determined by analysis of LINE-1 sequences and was higher in B-lymphocytes than matched T-lymphocytes (69.8% vs. 65.2%, P ≤ 0.01). Pairwise analysis identified 679 CpGs, representing 250 genes, which were differentially methylated between T-lymphocytes and B-lymphocytes. The majority of sites (76.6%) were hypermethylated in B-lymphocytes. Pyrosequencing of selected candidates confirmed the array data in all cases. Hierarchical clustering revealed perfect segregation of samples into two distinct clusters based on cell type. Differentially methylated genes showed enrichment for biological functions/pathways associated with leukocytes and T-lymphocytes. Our work for the first time shows that T-lymphocytes and B-lymphocytes possess intrinsic differences in DNA methylation within a restricted set of functionally related genes. These data provide a foundation for investigating DNA methylation in diseases in which these cell types play important and distinct roles.