BACKGROUND: Multiple factors, including interactions between genetic and environmental risks, are important in susceptibility to rheumatoid arthritis (RA). However, the underlying mechanism is not fully understood. This study was undertaken to evaluate whether DNA methylation can mediate the interaction between genotype and smoking in the development of anti-citrullinated peptide antibody (ACPA)-positive RA.
METHODS: We investigated the gene-smoking interactions in DNA methylation using 393 individuals from the Epidemiological Investigation of Rheumatoid Arthritis (EIRA). The interaction between rs6933349 and smoking in the risk of developing ACPA-positive RA was further evaluated in a larger portion of the EIRA (1119 controls and 944 ACPA-positive patients with RA), and in the Malaysian Epidemiological Investigation of Rheumatoid Arthritis (MyEIRA) (1556 controls and 792 ACPA-positive patients with RA). Finally, mediation analysis was performed to investigate whether DNA methylation of cg21325723 mediates this gene-environment interaction on the risk of developing of ACPA-positive RA.
RESULTS: We identified and replicated one significant gene-environment interaction between rs6933349 and smoking in DNA methylation of cg21325723. This gene-smoking interaction is a novel interaction in the risk of developing ACPA-positive in both Caucasian (multiplicative P value = 0.056; additive P value = 0.016) and Asian populations (multiplicative P value = 0.035; additive P value = 0.00027), and it is mediated through DNA methylation of cg21325723.
CONCLUSIONS: We showed that DNA methylation of cg21325723 can mediate the gene-environment interaction between rs6933349 and smoking, impacting the risk of developing ACPA-positive RA, thus being a potential regulator that integrates both internal genetic and external environmental risk factors.
Human populations have been shaped by catastrophes that may have left long-lasting signatures in their genomes. One notable example is the second plague pandemic that entered Europe in ca. 1,347 CE and repeatedly returned for over 300 years, with typical village and town mortality estimated at 10%-40%.1 It is assumed that this high mortality affected the gene pools of these populations. First, local population crashes reduced genetic diversity. Second, a change in frequency is expected for sequence variants that may have affected survival or susceptibility to the etiologic agent (Yersinia pestis).2 Third, mass mortality might alter the local gene pools through its impact on subsequent migration patterns. We explored these factors using the Norwegian city of Trondheim as a model, by sequencing 54 genomes spanning three time periods: (1) prior to the plague striking Trondheim in 1,349 CE, (2) the 17th-19th century, and (3) the present. We find that the pandemic period shaped the gene pool by reducing long distance immigration, in particular from the British Isles, and inducing a bottleneck that reduced genetic diversity. Although we also observe an excess of large FST values at multiple loci in the genome, these are shaped by reference biases introduced by mapping our relatively low genome coverage degraded DNA to the reference genome. This implies that attempts to detect selection using ancient DNA (aDNA) datasets that vary by read length and depth of sequencing coverage may be particularly challenging until methods have been developed to account for the impact of differential reference bias on test statistics.