Abstract
BACKGROUND: Up to 20% of Lyme disease cases experience post-treatment Lyme disease syndrome (PTLDS). The biological basis for PTLDS is poorly understood and no evidence-based treatment has been identified. Genetic studies have the potential to elucidate PTLDS pathophysiology and identify treatment targets. METHODS: We used electronic health record data (EHR) and genetic data from a linked biorepository to conduct a genome-wide association study (GWAS) for PTLDS among patients from a Pennsylvania health system. We evaluated the validity of the GWAS results in two separate conditions that have hypothesized overlapping pathophysiology, fibromyalgia and myalgic encephalomyelitis - chronic fatigue syndrome (ME/CFS). GWAS analyses were performed using logistic regression in SUGEN, assuming an additive genetic model, and adjusting for age, sex, array, and the first 10 principal components calculated from whole genome genotyping to adjust for ancestry, and accounting for relatedness including all 1st degree relationships. The functional mapping and annotation analysis (FUMA) tool was used to explore top findings from our GWAS. RESULTS: Among the 161,875 eligible MyCode participants with genotyping, there were 3,585 who met the criteria for treated Lyme disease. A subset of 695 (19.4%) of these patients met the criteria for PTLDS and the remaining 2890 were classified as controls. We identified two PTLDS loci that reached the suggestive significance threshold (P < 5 × 10(- 7)), with lead variants rs77857587, near IRX1, and rs10833979, near GAS2. Our top index single nucleotide polymorphism (SNP), rs77857587, is in high linkage disequilibrium with a long-range protein quantitative locus SNP, rs111774530, for the MARC2 (Mitochondrial Amidoxime Reducing Component 2) protein. We identified 5,041 cases of fibromyalgia (150,599 controls) and 2,268 cases of ME/CFS (151,594 controls) among the MyCode participants. Neither of the two suggestively significant loci were associated with fibromyalgia or ME/CFS. CONCLUSION: We identified two PTLDS loci that reached a suggestive significance threshold. Our top index SNP is associated with the MARC2 protein, a protein that has been linked to multiple immune checkpoints. Further study is needed in a larger population to evaluate whether there is genetic evidence of the role of immune response in the occurrence of PTLDS.