Objectives The main objective of this study was to determine whether the DNA methylation profile of children born to mothers with rheumatoid arthritis (RA) is different from that of children born to mothers from the general population. In addition, we aimed to determine whether any differences in methylation are associated with maternal RA disease activity or medication use during pregnancy.
Methods For this study, genome-wide DNA methylation was measured at cytosine-phosphate-guanine (CpG) sites, using the Infinium Illumina HumanMethylation 450K BeadChip, in 80 blood samples from children (mean age=6.8 years) born to mothers with RA. As controls, blood samples from 354 children (mean age=6.0 years) from the population-based Generation R Study were used. Linear mixed models were performed to investigate differential methylation between the groups, corrected for relevant confounders.
Results A total of 147 CpGs were differentially methylated between blood samples of children born to mothers with RA and the control blood samples. The five most significantly associated CpGs were cg06642177, cg08867893, cg06778273, cg07786668 and cg20116574. The differences in methylation were not associated with maternal RA disease activity or medication use during pregnancy.
Conclusions DNA methylation at 147 CpGs differed between children born to mothers with RA and children born to mothers from the general population. It remains unknown whether the identified associations are causal, and if so whether they are caused by the disease or treatment. More research, including replication of these results, is necessary in order to strengthen the relevance of our findings for the later-life health of children born to mothers with RA.
- rheumatoid arthritis
- dna methylation
This is an open access article distributed in accordance with the Creative Commons Attribution 4.0 Unported (CC BY 4.0) license, which permits others to copy, redistribute, remix, transform and build upon this work for any purpose, provided the original work is properly cited, a link to the licence is given, and indication of whether changes were made. See: https://creativecommons.org/licenses/by/4.0/.
Statistics from Altmetric.com
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.
What is already known about this subject?
Adverse exposures in early life are associated with later-life health.
Epigenetic changes are thought to be one of the underlying mechanisms.
There is not much known about the consequences of maternal rheumatoid arthritis (RA) on the offsprings’ long-term health.
What does this study add?
DNA methylation is different in children born to mothers with RA compared with mothers from the general population.
How might this impact on clinical practice or future developments?
Maternal RA disease during pregnancy might have lifelong consequences for the offspring.
More research in this particular field must be undertaken.
Adverse exposures in early life are associated with later-life health, which is referred to as the developmental origins of health and disease hypothesis.1–5 Epigenetic processes are thought to be one of the mechanisms underlying the associations of early-life exposures and later-life health outcomes.6 7 DNA methylation is the best studied and understood epigenetic modification.8 9 Factors that have been demonstrated to be associated with fetal DNA methylation include maternal disease,7 malnutrition,10–13 smoking,14 placental insufficiency,15 corticosteroids,16 folate depletion17 and cytokines.18 DNA methylation usually occurs at cytosine-phosphate-guanine (CpG) sites.8 9 The effect of hypermethylation and hypomethylation on gene expression depends on the CpG location.19 20 The most pronounced changes in DNA methylation occur during early pregnancy.7 11
During embryogenesis, there are three germ layers that form in the developing fetus (ectoderm, mesoderm and endoderm). When DNA methylation is altered in early pregnancy, all germ layers are affected.21 22
Rheumatoid arthritis (RA) may be considered as an adverse exposure during pregnancy.23 Therefore, it is plausible that maternal RA may induce changes in fetal DNA methylation, and that it is related with the later-life health of the offspring. Interleukin-6 is known to influence DNA methylation.18 RA treatment during pregnancy includes among others sulfasalazine (SSZ) and corticosteroids such as prednisone. SSZ is a known folate antagonist that crosses the placenta and could influence DNA methylation in this respect.17 24 Furthermore, corticosteroids might influence DNA methylation.16 Especially during early pregnancy, when the placenta is not completely developed, prednisone passively diffuses to the fetus.25–28
In the current study we investigated whether the DNA methylation profile of children born to mothers with RA was different from that of children born to mothers from the general population. Furthermore, we investigated whether any differentially methylated CpGs were associated with RA disease activity or medication use during pregnancy, or with indicators of future metabolic and cardiovascular diseases. In addition, we examined whether these CpGs were associated with the expression of genes using expression quantitative trait methylation (eQTM) analysis.
This study is embedded in the Pregnancy-induced Amelioration of Rheumatoid Arthritis (PARA) study, a prospective cohort study on pregnancy and RA.29 From 2002 to 2008, 369 female patients with RA who had a wish to conceive (or already pregnant) were enrolled.30 31 After participation in the PARA study, 196 children and their parents were invited to participate in a follow-up study, the FEtal Programming in Rheumatoid Arthritis (FEPRA) study. For this study, 108 children with a mean age of 6.8 years (SD=1.3) visited Erasmus Medical Centre in Rotterdam, and the parents of 85 children (all of European ancestry) gave informed consent for studies on DNA methylation of their children. Furthermore, the parents of 71 children provided cheek swabs from their children. There were no statistical differences in baseline characteristics between the participating and non-participating group.
Generation R Study
The control group consisted of children with a mean age of 6.0 years (SD=0.4), included in the Generation R Study, a population-based prospective cohort study from pregnancy onwards in Rotterdam, the Netherlands.32 In this study, all pregnant women living in Rotterdam with a delivery date between April 2002 and January 2006 were invited to participate, and 9778 mothers were enrolled.32 At the age of 6 years, DNA methylation was measured in a subgroup of 493 children of European ancestry.
In the PARA study, data on mother (eg, disease activity (with the Disease Activity Score in 28 joints using C reactive protein levels, DAS28-CRP3)) and child were collected.31 For the FEPRA study, data on blood pressure, growth and body composition of the children were measured. Also, blood, which is a mesoderm germ layer derivate, was drawn for DNA methylation analysis.33 Cheek swabs were collected for the analysis of DNA methylation in buccal epithelial cells, which is an ectoderm germ layer derivate.
Generation R Study
In the Generation R Study, mothers were seen three times during pregnancy. The children were followed from birth until childhood. Data collection in children and their mothers included questionnaires, detailed physical examinations and blood sampling.32
DNA methylation analysis
Genomic DNA was extracted from whole peripheral blood samples and from the cheek swab samples. Bisulfite conversion of 500 ng of genomic DNA was performed using the Zymo EZ-96 DNA Methylation Kit (Zymo Research, Irvine, California, USA) according to the manufacturer’s protocol.
Genomic methylation profiling was performed using the Infinium Illumina HumanMethylation 450K BeadChip arrays (Illumina, San Diego, USA) according to the manufacturer’s protocol. The Illumina array measures methylation status of 485 512 CpG sites in the gene and non-gene regions across the entire human genome. To prevent possible batch effects, blood and cheek swab samples were measured in the same run.
Quality control and normalisation
The data were preprocessed using the minfi package in R V.3.4.1 (www.r-project.org). Samples with incomplete or poor bisulfite conversion, extension, hybridisation or specificity were excluded.34 In addition, samples with sex mismatch and samples with a call rate <95% were removed. This quality control (QC) was done separately for blood samples and for cheek swab samples. During QC, 5 blood and 14 cheek swab samples from the FEPRA study were excluded, resulting in 80 and 57 samples, respectively. From the Generation R blood samples, 27 were excluded due to corticosteroid use or RA disease in the mother, and 32 were excluded during QC, resulting in 441 blood samples. In addition, 87 cases with missing data from the Generation R Study were excluded, leaving 354 samples to analyse. The intensity values were then quantile normalised using the incorporating Control Probe Adjustment and reduction of global CORrelation (CPACOR) workflow.34 Methylation at each CpG was calculated as the beta value (beta=intensity of the methylated allele (M)/(intensity of the unmethylated allele (U)+M+100)), containing values from 0 to 1. Blood cell composition of the samples was estimated using the Houseman method with the Reinius reference set.35 36 The Reinius reference set is not yet validated in children. However, it is the best method available, and it has been used in other epigenetic studies in children.37 38 Probes with single nucleotide polymorphisms (SNPs) at single base extension, probes with improper binding, and CpGs on the X and Y chromosome were removed from the analysis.37–40 From the initial 485 512 CpGs, 32 057 were excluded during QC, leaving 453 456 CpGs for analysis.
For all subjects, descriptive statistics were calculated using Stata V.14.1 (https://www.stata.com/stata14/). Student’s t-tests and χ2 tests were used to compare the baseline characteristics. For these analyses, p values <0.05 were considered statistically significant.
Linear mixed models were performed to analyse differential methylation between the groups, using R. The first model was created to compare the blood samples from the FEPRA study with the blood samples of the Generation R Study to determine whether the DNA methylation profile of children born to mothers with RA was different from that of children born to mothers from the general population. This model was corrected for biological covariates (age, body mass index (BMI) SD scores (SDS), adjusted for age and sex according to the Dutch reference values, using the Growth Analyser (V.4.0; Growth Analyser, Rotterdam, the Netherlands, www.growthanalyser.org,) sex, gestational age at delivery, maternal age, folic acid use during pregnancy, socioeconomic status (SES), maternal smoking and white blood cell subtypes), technical covariates (technical batch effects (array identifier (ID) and position on array)) and five hidden confounders. Technical covariates were added as random effects in the models. The hidden confounders were calculated using the CATE package41 42 while correcting for the group (RA vs non-RA offsprings), all biological covariates and technical covariates. This resulted in hidden confounders that were independent of all included covariates. The BACON package42 was used to correct for unobserved covariates in order to reduce test statistic bias and inflation. The genomic inflation factor (λ)43 was calculated. After QC, 453 456 CpGs remained for analysis. Therefore, a Bonferroni-adjusted p value of 0.05/453 456=1.10×10− 7 was used.44
CpGs were annotated for nearby genes with the Genomic Regions Enrichment of Annotations Tool,45 a general method to present epigenetic data.37 38 40 Significant CpGs from the first analysis were then further analysed within the blood samples of the FEPRA study to explore if RA disease activity, prednisone use or SSZ use during pregnancy would explain the differences in methylation found in the first model. After that, two linear mixed models were created with the significant CpGs to identify whether these were associated with the BMI SDS or the fat percentage SDS of the child, as indicators of future metabolic disease. These models were analysed in the 80 blood samples of the FEPRA study. CpGs with a p value below 0.05/147=3.4×10−4 were considered statistically significant.
The significant CpGs found in the first analysis were also analysed in DNA derived from buccal epithelial cells, obtained by cheek swabs from the FEPRA study, to explore if the differentially methylated CpGs were also differentially methylated in that germ layer derivate, as a kind of validation of the results. For this analysis, CpGs with a p value below significance level of 3.4×10−4 were considered statistically significant.
eQTMs are sites at which DNA methylation is known to influence the expression of one or more genes.46 To analyse whether any of the significant CpGs were linked to the expression of nearby genes, an eQTM analysis was performed. For these analyses the online BIOS QTL browser (https://molgenis26.target.rug.nl/downloads/biosqtlbrowser/) was used.39
The flow chart of the study population is shown in figure 1. A total of 80 blood and 57 cheek swab samples from the FEPRA study (children born to mothers with RA) and 354 blood samples from the Generation R Study (children born to mothers from the general population) remained for analysis.
Descriptive statistics of the study population are presented in table 1. Overall, children from the FEPRA study were slightly older compared with children in the Generation R Study (p<0.001). In the FEPRA study, 45% of the women did not start using folic acid before or in early pregnancy compared with 92.7% in the Generation R Study (p<0.001). Approximately half of the women (47.5%) from the FEPRA study had a high SES compared with 68.4% of the women from the Generation R Study (p<0.001). One woman (1.3%) from the FEPRA study and 86 women (24.3%) from the Generation R Study smoked periconceptionally or at any time during pregnancy (p<0.001).
DNA methylation analysis
In the first linear mixed model, blood samples from the FEPRA study were compared with blood samples from the Generation R Study, corrected for the covariates mentioned in the Methods section. In total, 147 CpGs were significantly differentially methylated between children in the FEPRA study and children in the Generation R Study (figure 2). The QQ plot is shown in figure 3. The genomic inflation factor (λ) was 1.06.
In table 2 the five most significant CpGs and the five CpGs with the largest effect sizes (within the significant CpGs) are described in detail. The complete list of significantly associated CpGs is provided as an online supplementary file. A positive beta represents higher methylation levels in children born to mothers with RA. Methylation was higher in children born to mothers with RA than in those from the general population at all five CpGs with the lowest p values. Methylation at the five CpGs with the largest effect sizes was lower in children born to mothers with RA, with the exception of cg06656994.
Subsequent analysis of DNA methylation
The mean DAS28-CRP(3) in the first, second and third trimesters (3.6, 3.5 and 3.3, respectively) were highly correlated (>0.6). The DAS28-CRP(3) in the third trimester was available in all patients. Therefore, this timepoint was chosen for the analysis. None of the 147 CpGs were significantly associated with maternal RA disease activity (DAS28-CRP(3)) in the third trimester or medication (prednisone or SSZ) use. In addition, none of the CpGs were associated with BMI SDS or fat percentage SDS in the children.
Analysis in buccal epithelial cells
A total of 10 out of the 147 CpGs significantly associated with maternal RA in blood were also associated in buccal epithelial cells. From these, four were in the same direction as in blood (table 3). CpG cg11336323 was located in a promoter region.
Two CpGs, cg21384971 and cg11220663, were associated with expression of the COPZ2 and ADD2 genes, respectively (table 4). These two genes were also the nearest genes to those CpGs.45 Both CpGs were hypermethylated in the children born to mothers with RA and were associated with decreased expression of COPZ2 and ADD2 in the BIOS eQTM lookup browser.
This is the first study investigating the differences in DNA methylation of children born to mothers with RA compared with children born to mothers from the general population. In this unique study, all participants were followed prospectively from pregnancy onwards. Our study showed differential DNA methylation between the two groups. The differentially methylated CpG sites were not associated with disease activity and/or medication use, nor to BMI SDS and fat percentage SDS.
In total, 147 CpGs were significantly associated with maternal RA after adjustment for multiple biological and technical covariates and hidden confounders. Of the five most significant CpGs, interestingly, two (cg06642177 and cg07786668) have been associated with myocardial infarction.47 The most significant CpG, cg06642177, is located on chromosome 6 near the SLC2A12 gene,48 associated with insulin sensitivity,48 49 heart failure and diabetes50 in animal models. cg07786668, located on chromosome 16, is located in the ZFHX3 gene. ZFHX3 has been associated in multiple human studies with atrial fibrillation,51–55 coronary heart disease56 and obesity in a Korean population.57 cg20116574 was annotated to the NCOA5, a protein coding gene, which has been associated with diabetes mellitus type 2 in animal models.58 59
From the remaining 142 CpGs, 1 (cg17218495), annotated to the SMARCA4 gene, has independently and significantly been associated with myocardial infarction.47 The other significant CpGs from our study have not been associated with disease phenotypes.
Out of the initial 147 significant CpGs, 10 were also significantly differentially methylated in buccal epithelial cells obtained by cheek swabs. From these, four were in the same direction as in blood. When DNA methylation is altered in more than one germ layer derivate, it is likely that these alterations occurred in early development.21 22 These four CpGs have not been associated with disease phenotypes in humans. Unfortunately, in the Generation R Study, DNA methylation in buccal epithelial cells was not available. There were no publicly available epigenetic data sets on buccal epithelial cells in healthy children with a similar age and background.
As mentioned before, CpGs that are associated with eQTMs can influence the expression levels of genes.46 In our study, from the 147 significant CpGs, 2 were present in the BIOS QTL browser. cg21384971 is associated with less expression of the COPZ2 gene, which has been studied as a therapeutic opportunity for proliferation-independent selective killing of tumour cells.60 cg11220663 is associated with less expression of the ADD2 gene, also known as beta-adducin. ADD2 is involved in multiple pathogenic processes with a wide range of diseases.61 ADD2 gene variants are associated with hypertension,62 cancer61 and systemic lupus erythematosus.63
Pathway analysis of the genes annotated to the 147 significant CpGs, using the WebGestalt (WEB-based Gene Set AnaLysis Toolkit),64 did not result in significant pathways (data not shown).
Thus, some of the associated CpGs (cg06642177, cg07786668 and cg17218495) have been associated with cardiovascular disease in previous studies, while others are located in or near genes that are associated with cardiovascular or metabolic disease. Maternal RA during pregnancy, a chronic inflammatory disease, might be associated with later-life health and disease risk in the offspring.
None of the significant CpGs were associated with RA disease activity or medication use during pregnancy. However, this may have been due to a lack of power, since these analyses were performed in the 80 samples of the children born to women with RA. The same power problem also applied to the analysis of the CpGs with indicators for future metabolic and cardiovascular disease (BMI SDS and fat percentage).
Remarkably, a large percentage of women with RA did not use folic acid before or during pregnancy. Even though this is outside the scope of our study, this requires additional attention during preconceptional counselling.
Our study has some limitations. First, although in its kind it is a large study, a study on DNA methylation including 80 subjects and 354 controls is still relatively small. Despite this, a large number of CpGs reached significance. Correcting for biological and technical covariates, as well as hidden confounders, and using BACON resulted in a λ near 1, reflecting that there was no inflation. Second, we were not able to collect a new independent cohort of children born to mothers with RA to replicate the results. At the time our study was performed, there were no other comparable prospective studies available. Currently, European research groups are conducting new prospective cohort studies on the impact of RA on pregnancy and offspring. We encourage these research groups, possibly with international collaborations, to replicate our study.
Thus, the results of this study may support follow-up research of children born to mothers with RA. Based on our data, we recommend that at least indicators for future cardiovascular and metabolic disease should be considered. The effects of RA disease activity and medication use on DNA methylation should be investigated in studies with larger sample sizes. Furthermore, in the last years the use of tumour necrosis factor (TNF) inhibitors during pregnancy in patients with RA has increased. This often results in lower RA disease activity during pregnancy. Future research should also cover the effects of the use of TNF inhibitors on the differentially methylated CpGs in children born to mothers with RA.
In addition, it would be interesting for a future study to measure DNA methylation in the mothers of the children from our study and compare that with the methylation of their offspring. Since mothers and children live in a shared environment (generally), there may be an overlap in differentially methylated CpGs.
In conclusion, maternal RA during pregnancy is associated with differential DNA methylation in offspring. It remains unknown whether the identified associations are causal, and if so whether they are caused by the disease or treatment. Some of the differentially methylated CpGs or their nearby genes were associated with cardiovascular or metabolic disease. Maternal RA disease might have lifelong consequences for the offspring. However, more research in this particular field must be undertaken in order to strengthen the relevance of our findings.
We are grateful to all patients participating in the PARA study, FEPRA study and the Generation R Study. We thank all research assistants and laboratory workers for their help with data collection and processing of the samples. We gratefully acknowledge the contribution of Dr Florentien de Steenwinkel for her work on the FEPRA study. We thank Dr Bas Heijmans for his advice on the statistical analysis of the data. The Generation R Study is conducted by the Erasmus Medical Centre in close collaboration with the School of Law and Faculty of Social Sciences of the Erasmus University Rotterdam, the Municipal Health Service Rotterdam area, Rotterdam, the Rotterdam Homecare Foundation, Rotterdam, and the Stichting Trombosedienst & Artsenlaboratorium Rijnmond (STAR-MDC), Rotterdam. We gratefully acknowledge the contribution of children and parents, general practitioners, hospitals, midwives and pharmacies in Rotterdam. The generation and management of the Illumina 450K methylation array data (EWAS data) for the Generation R Study were executed by the Human Genotyping Facility of the Genetic Laboratory of the Department of Internal Medicine, Erasmus MC, the Netherlands. We thank Mr Michael Verbiest, Ms Mila Jhamai, Ms Sarah Higgins, Mr Marijn Verkerk and Dr Lisette Stolk for their help in creating the EWAS database. We thank Dr A Teumer for his work on the quality control and normalisation scripts.
Handling editor Josef S Smolen
Correction notice This article has been corrected since it first published online. The open access licence type has been amended.
Contributors JMH, RJEMD, JBvM, JF, LD, HI-A: study design and analysis plan. JMH, RJEMD, JF: data acquisition. HI-A, PM, JF: cleaning of data. PM, HI-A: data analysis. HI-A, RJEMD: draft paper. JMH, RJEMD, JF, JBvM, PM, LD, HI-A: revision of the paper and final approval for publication. The authors confirm that they all meet the Uniform Requirements for Manuscripts Submitted to Biomedical Journals criteria for authorship.
Funding This study was funded by the Dutch Arthritis Society (ReumaNederland), a non-commercial fund-raising organisation. The general design of the Generation R Study is made possible by financial support from the Erasmus Medical Centre, Rotterdam, the Erasmus University Rotterdam, the Netherlands Organisation for Health Research and Development, and the Ministry of Health, Welfare and Sport. The EWAS data was funded by a grant to Vincent W Jaddoe (Generation R Study group) from the Netherlands Genomics Initiative (NGI)/Netherlands Organisation for Scientific Research (NWO), Netherlands Consortium for Healthy Ageing (NCHA; project nr 050-060-810), by funds from the Genetic Laboratory of the Department of Internal Medicine, Erasmus MC, and by a grant from the National Institute of Child and Human Development (R01HD068437). LD received funding from the European Union's Horizon 2020 cofunded programme ERA-Net on Biomarkers for Nutrition and Health (ERA HDHL) (ALPHABET project (no 696295; 2017), ZonMW The Netherlands (no 529051014; 2017)). The Generation R Study was supported by funding from the European Union’s Horizon 2020 research and innovation programme (733206, LifeCycle Project).
Competing interests None declared.
Patient consent for publication Not required.
Ethics approval This study is in compliance with the Helsinki Declaration. Informed consent was obtained for all participants. The Medical Ethics Committee of Erasmus MC, University Medical Center Rotterdam, the Netherlands, approved the PARA study (MEC-214.320/2002/117), FEPRA study (MEC-2009-441) and the Generation R Study (MEC198.782.2001.31).
Provenance and peer review Not commissioned; externally peer reviewed.
Data availability statement Data may be obtained from a third party and are not publicly available.