The 3p21.31 genetic locus promotes progression to type 1 diabetes through the CCR2/CCL2 pathway =============================================================================================== * Paul Tran * Sharad Purohit * Eileen Kim * Khaled bin Satter * Diane Hopkins * Kathleen Waugh * Fran Dong * Suna Onengut-Gumuscu * Stephen S. Rich * Marian Rewers * Jin-Xiong She ## Abstract **Objective** Multiple cross-sectional and longitudinal studies have shown that serum levels of the chemokine ligand 2 (CCL-2) are associated with type 1 diabetes (T1D), although the direction of effect differs. We assessed CCL-2 serum levels in a longitudinal cohort to clarify this association, combined with genetic data to elucidate the regulatory role of CCL-2 in T1D pathogenesis. **Research design and methods** The Diabetes Autoimmunity Study in the Young (DAISY) followed 310 subjects with high risk of developing T1D. Of these, 42 became persistently seropositive for islet autoantibodies but did not develop T1D (non-progressors); 48 did develop T1D (progressors). CCL-2 serum levels among the three study groups were compared using linear mixed models adjusting for age, sex, HLA genotype, and family history of T1D. Summary statistics were obtained from the CCL-2 protein quantitative trait loci (pQTL) and *CCR2* expression QTL (eQTL) studies. The T1D fine mapping association data were provided by the Type 1 Diabetes Genetics Consortium (T1DGC). **Results** Serum CCL-2 levels were significantly lower in both progressors (p=0.004) and non-progressors (p=0.005), compared to controls. Two SNPs (rs1799988 and rs746492) in the 3p21.31 genetic locus, which includes the CCL-2 receptor, *CCR2*, were associated with increased *CCR2* expression (p = 8.2e-5 and 5.2e-5, respectively), decreased CCL-2 serum level (p = 2.41e-9 and 6.21e-9, respectively), and increased risk of T1D (p = 7.9e-5 and 7.9e-5, respectively). **Conclusions** The 3p21.31 genetic region is associated with developing T1D through regulatory control of the CCR2/CCL2 immune pathway. Keywords * CCR2 * CCL2 * MCP1 * T1D * DAISY * autoimmune * fine mapping ## INTRODUCTION Type 1 diabetes (T1D) is an autoimmune disorder where some genetically predisposed individuals are exposed to an environmental trigger which leads to immune recognition and destruction of pancreatic islet cells1. T1D is associated with microvascular and macrovascular complications and a shortened lifespan2,3. A primary goal in T1D is primary prevention of the disease. Recent success in anti-CD3 inhibition delaying progression from islet autoantibody positivity to T1D has brought promise to immunomodulation for T1D prevention4. The use of genome wide association studies (GWAS) may help identify other immune targets5,6. GWAS have associated the 3p21.31 locus with T1D and pediatric autoimmune diseases (including T1D)6,7; however, the causal variants and mechanism of association have not been established for this locus. This region of the human genome is enriched for chemokine receptors that play an important role in autoimmunity8, including the C-C motif chemokine receptors (*CCR1, CCR2, CCR3, CCR5, CCR9, CCRL2, CXCR6*) and X-C motif chemokine receptor 1 (*XCR1*). Associated pathways include chromatin remodeling genes (*SETD2, ELP6*, and *SMARCC1*), putative transcriptional regulators (*CCDC12, CCDC51, CCDC36*, and *ZNF589*), and cell replication associated genes (*KIF9* and *CDC25A*)9, as well as the *CCRL2*-*LINC02009*6 and *DAG1*7 loci. CCR2 is a pro-inflammatory chemokine receptor which, upon binding with its ligand CCL-2, promotes effector T-cell activation10. Non-obese diabetic (NOD) *CCR2*-/-mice develop T1D at a slower rate than NOD wild-type mice11. Serum CCL2 levels are different between T1D and controls in case-control studies; however, the direction of difference is contested12-14. Lower serum levels of CCL-2 have been reported in children with persistent islet autoantibodies compared to controls15, suggesting a potential causative role of *CCR2* in T1D pathogenesis16,17. *CCR2* has also been implicated in other immune diseases, including rheumatoid arthritis (RA)18-20 and multiple sclerosis (MS)21-23. However, clinical trials of CCR2 antagonists did not show reduce disease progression in either RA or MS patients24,25. CCR2 inhibition is associated with increased serum levels of its ligand, CCL2 or MCP-1, in both mouse and human studies24,26,27. CCR2 activation by CCL2 has been shown to cellular uptake and internalization of the CCR2/CCL2 complex28. A recent GWAS identified SNPs near *CCR2* in the 3p21.31 genetic locus associated with serum CCL2 levels. These studies suggest that *cis*-acting genetic variants which alter CCR2 expressivity could paradoxically change serum CCL2 levels29. In this report, the direction of the association between CCL2 level and T1D progression is determined in a longitudinal cohort (DAISY) and potential causal variants in the 3p21.31 locus were identified. ## RESULTS ### Serum CCL-2 levels decrease over time in DAISY participants who develop T1D compared to controls In this study, we measured CCL-2 levels in longitudinally collected serum samples from 310 children at high risk for T1D participating in the Diabetes Autoimmunity Study in the Young (DAISY)30,31. At the time of analysis, 48 DAISY participants developed persistent islet autoantibodies and progressed to T1D (progressors) and an additional 42 participants developed persistent islet autoantibodies but have not progressed to T1D (non-progressors). Their longitudinal CCL-2 levels were compared to those of 220 islet autoantibody negative children (controls) similar to the cases in terms of age, sex, family history of T1D and HLA genotypes (**Table 1)**. View this table: [Table 1.](http://medrxiv.org/content/early/2021/08/05/2021.08.01.21261449/T1) Table 1. Characteristics of study participants by study group. Categorical variables were analyzed using Pearson χ2 tests. Continuous variables were tested using the t test for differences in means or the Wilcoxon rank sum test for differences in medians The subjects were followed for 16.80 years on average (range 1.98 to 31.00 years). On average, 9.34 (range 3 to 30) serial serum samples were collected from each subject. The mean (SD) age at the appearance of islet autoantibodies was 4.59± 3.25 years in progressors and 7.54± 4.24 years in non-progressors. The mean age of T1D diagnosis in progressors was 12.22+ 5.68 years. CCL-2 was detectable in all samples from the 310 DAISY subjects. CCL-2 levels decreased with age in all groups. We next tested whether the average CCL-2 levels differed among the three study groups, using linear mixed models adjusting for age, sex, HLA genotype, FDR status (**Table 2**). CCL-2 levels were significantly lower in both progressors (p=0.004) and non-progressors (p=0.005), compared to controls (**Supplemental Figure 1**). Subjects with high-risk HLA genotype (DR3/4 or DR4/4 or DR4/X) had lower MCP1 level compared to subjects with moderate-risk HLA genotype DR3/3 or DR3/X or DRX/X (p=0.009). However, CCL-2 levels did not differ by sex or by the presence of a first-degree relative. ![Supplemental Fig. 1.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/08/05/2021.08.01.21261449/F4.medium.gif) [Supplemental Fig. 1.](http://medrxiv.org/content/early/2021/08/05/2021.08.01.21261449/F4) Supplemental Fig. 1. The following graph is spaghetti plot of the CCL-2 level over age among different groups with total number of subjects=310. ![Supplemental Fig 2.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/08/05/2021.08.01.21261449/F5.medium.gif) [Supplemental Fig 2.](http://medrxiv.org/content/early/2021/08/05/2021.08.01.21261449/F5) Supplemental Fig 2. Modeled CCL-2 levels for controls (upper curve), progressors, and non-progressor from the mixed model with a quadratic curve over age. Fit computed for male and HLA group of DR3/3 or DR3/X or DRX/X and GP. View this table: [Table 2.](http://medrxiv.org/content/early/2021/08/05/2021.08.01.21261449/T2) Table 2. Multivariate mixed models analysis of the CCL-2 levels ![Figure 1.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/08/05/2021.08.01.21261449/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2021/08/05/2021.08.01.21261449/F1) Figure 1. GWAS for CCL-2 serum levels and fine mapping for type 1 diabetes. A) Manhattan plot of the serum CCL-2 protein quantitative trait loci (pQTL) from Folkersen et al29. Data were downloaded from Xenodo and plotted using the LocusZoom plot package. The genome wide significant SNPs are from the 3p21.31 region. B) LocusZoom plot of SNPs significance on T1D association testing from T1DGC fine mapping study. C) LocusZoom plot of SNPs significance on CCL-2 pQTL from Figure 1A. D) LocusZoom plot of SNPs significance on CCR2 eQTL from GTEx. The potential causal variants rs1799988 (chr3:46412259) and rs746492 (chr3:46417312) are indicated with arrows. ### The 3p21.31 genetic locus is associated with serum CCL-2 levels Genetic variants associated with serum CCL-2 levels were determined by downloading summary statistics from a study estimating cis- and trans-protein quantitative trait loci (pQTLs) for over 600 proteins using over 20,000 subjects 29. One locus contributing to serum CCL-2 levels at 3p21.31 was identified **(Figure 1A)**, centered around the *CCR2* gene **(Figure 1B)**. The SNPs on chromosome 3 accounted for ∼10% of the variation in CCL-2 serum level, which explains 70% of estimated heritability29. An additional study in the GWAS catalog32 validated the association between the 3p21.31 locus and serum CCL-2 levels33. ### Potential causal variants for serum CCL-2 level regulation and T1D pathogenesis Since the 3p21.31 locus is associated with both CCL-2 serum levels and T1D, and CCL-2 serum levels are associated with T1D, genetic variants associated with T1D through regulatory control of CCL-2 serum levels were targeted for evaluation. Multiple levels of data were integrated to identify a potential set of causal variants from the 3p21.31 locus. Summary statistics were obtained from the T1DGC ImmunoChip fine mapping data6, which profiled 135,870 custom selected SNPs on 18,892 subjects. After merging the summary statistics from CCL-2 protein quantitative trait loci (pQTL) and T1DGC data, a total of 8,804 SNPs on chromosome 3 overlapped both the CCL-2 and T1DGC studies, including 576 SNPs between position chr3:46,300,000 - 46,800,000 (**Figure 1C**). Using a Bonferroni-adjusted p-value cut-off, 225 SNPs were associated (p < 8.7 x10−5) with CCL-2 serum levels, with 63 SNPs associated with T1D. A total of 47 SNPs were associated with both T1D and CCL-2 serum levels. Expression QTLs (eQTLs) associated with *CCR2* were obtained from the Genotype-Tissue Expression (GTEx) portal for whole blood. Eighty-six SNPs overlapped the combined CCL-2 pQTL and T1DGC fine mapping data and the GTEx whole blood eQTL data at chr 3:46,300,000 - 46,800,000 (**Figure 1D**). Two SNPs, rs1799988 and rs746492, were significantly associated with T1D, CCL2 serum pQTL, and *CCR2* eQTLs (**Figure 2**), suggesting that the two SNPs could play a regulatory role in CCR2/CCL2 pathway activation and progression to T1D. We therefore propose a mechanism linking the two potential causal variants to T1D pathogenesis. Individuals with these variants have associated increased immune cell expression of *CCR2* and the associated increase in CCR2/CCL2 signaling, which causes immune-mediated islet destruction and depletion of serum CCL-2 levels **(Figure 3)**. ![Figure 2.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/08/05/2021.08.01.21261449/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2021/08/05/2021.08.01.21261449/F2) Figure 2. Forest plot of potential causal SNPs (rs1799988 and rs746492) *CCR2* eQTL, CCL2 pQTL, and T1D Immunochip Association. The axis represents effect size with bars showing 95% confidence interval. ![Figure 3.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/08/05/2021.08.01.21261449/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2021/08/05/2021.08.01.21261449/F3) Figure 3. Proposed model of 3p21.31 locus influence on progression to type 1 diabetes. Individuals with high risk *CCR2* genotypes are more prone to increased *CCR2* expression, increased CCR2/CCL2 axis signaling with subsequent depletion of serum CCL2, increased monocyte recruitment, and increased rate of pancreatic islet destruction, resulting in an increased rate of progression to type 1 diabetes. Created with BioRender.com. ## DISCUSSION While some studies have shown that serum CCL-2 levels are increased in T1D subjects compared to normal controls, our previous work and current findings do not agree with this12-14. By analyzing serum CCL-2 levels in a longitudinal cohort, we were able to demonstrate the decrease in serum CCL-2 levels as subjects progressed from autoantibody positivity to T1D **(Table 2)**. The decrease in serum CCL-2 with disease progression is in agreement with our proposed mechanism. *Cis*-acting genetic variants increase *CCR2* expression above the average population level which increases CCR2/CCL2 axis signaling and monocyte recruitment to the pancreatic islet. The increased signaling depletes the serum pool of CCL2 **(Figure 3)**. These findings contradict the simple logic that increased serum inflammatory markers indicate a pro-inflammatory state. A similar mechanism has been described where decreased IL-6 serum levels are associated with worse cardiovascular prognosis due to increased IL-6R/IL-6 signaling in the vascular wall34,35. Our findings provide an additional example of decreased serum levels of a pro-inflammatory molecule associated with an increased inflammatory state. Interestingly, serum CCL2 levels were more likely to be elevated at 3-4 months in infants who will develop islet autoimmunity in The Environmental Determinants of Diabetes in the Young36 (TEDDY) high risk genetic cohort (M. Rewers, personal communication). This may be due to early enterovirus infection, which has been shown to increase type 1 interferon and CCL2 expression37. A lack of statistical significance of decreased CCL2 levels in progression to islet autoimmunity and type 1 diabetes may be due to lack of statistical power in this cohort since 1) less samples were collected after seroconversion and 2) the cohort is a genetically pre-select group which may lack the traditional variation in CCL2 levels. We were able to narrow the list of potential causal variants by integrating the CCL-2 pQTL and T1D Immunochip fine mapping. Although this integration provides high confidence for these SNPs, they will require functional validation. Additionally, while the T1D fine mapping study used a custom SNP array for autoimmune diseases38, the CCL-2 GWAS and the GTeX eQTL study39 used a broader SNP array. Thus, SNPs which were not assessed in all SNP arrays could not be fairly evaluated and may still regulate the CCR2/CCL2 axis. Our findings support the potential of CCR2 inhibition in delaying the progression of islet autoantibody positive individuals to T1D. Several CCR2 inhibitors have been tested in clinical trials for cardiovascular disease, diabetic nephropathy, cancer, and autoimmunity40-44. Within autoimmunity, CCR2 inhibition did not change disease progression for subjects with rheumatoid arthritis or multiple sclerosis. One potential explanation is that CCR2 inhibition was targeted too late in the disease course and that earlier targeting based on serology and before clinical manifestations would modulate clinical progression45. Additionally, individuals with high CCL-2 serum levels, indicating a normal or low expression and activity of CCR2 would likely not benefit as much as individuals with high *CCR2* expression and activity with lower CCL-2 serum levels. Thus, potential future trials for T1D prevention should target autoantibody positive individuals with lower serum levels of CCL-2. In conclusion, we mapped the T1D association of the 3p21.31 locus to the *CCR2* gene. Increased *CCR2* gene expression decreases serum levels of CCL-2 and is associated with increased immune destruction of pancreatic islets and quicker progression to T1D. ## METHODS ### Human subjects The prospective cohort was part of the DAISY study30,31. The sex, age, family history of T1D, and HLA-DQB1 genotypes for all subjects were summarized in **Table 1**. Serum samples were collected in clot activator tubes and processed and stored within less than 1 hour. Repeated freezer/thawing cycles were avoided for all samples. Islet autoantibodies to insulin, GAD, IA-2 and ZnT8 were measured by radio-binding (RBA) assay. Persistence was defined by presence of at least one of the autoantibodies on at least two consecutive visits 3-6 months apart. ### Laboratory assays Serum levels of CCL-2 were measured using multiplex fluorescent immunoassays (Millipore, Billerica MA) according to the manufacturer’s recommendations. The coefficients of variation (CV, <10%) and the minimum detectable dose (MDD, CCL-2: 0.9) for the assays were provided by the manufacturer with the assay kits. Serum samples (5x diluted) were incubated with capture antibodies immobilized on polystyrene beads for one hour. The beads were then washed and further incubated with biotinylated detection antibody cocktail for one hour. Beads were washed twice to remove unbound detection antibody and then incubated with phycoerythrin-labeled streptavidin for thirty minutes. Prior to reading, the beads were washed 2 times and suspended in 60µl of wash buffer. The median fluorescence intensities (MFI) were measured on a FlexMAP 3D array reader (Millipore, Billerica, MA) using the following instrument settings: events/bead: 50, minimum events: 0, flow rate: 60µl/min, Sample size: 50ul and discriminator gate: 8000-13500. Before profiling, the serum dilutions were optimized by performing the assays at different serum dilutions to ensure that the majority of the data falls within the linear range of the standard curve. Protein concentrations were estimated using a regression fit to the standard curve with known concentration included on each plate using a serial dilution series. The data were subjected to several quality control steps before further analysis12. ### Data acquisition DAISY phenotype data was provided by MR. CCL-2 pQTL data were downloaded from Zenodo ([https://zenodo.org/record/2615265#.YGs9GEhKjYk](https://zenodo.org/record/2615265#.YGs9GEhKjYk)) on 05/01/2021. T1DGC summary statistics data were provided by SOG and SR6. CCR2 eQTL data were downloaded from GTEx portal ([https://www.gtexportal.org/home/gene/CCR2](https://www.gtexportal.org/home/gene/CCR2)) on 05/01/2021. ### Statistical analysis Statistical analyses of serum CCL2 levels were performed using SAS software version 9.4. A p-value less than 0.05 was considered significant. Because CCL-2 values are not normally distributed, Box-Cox transformation was performed to normalize the data. The optimal power parameter (lamda) was obtained using PROC TRANSREG by using a maximum likelihood criterion46. The optimal power parameter (Lamda) was found to be 0.75 and tCCL-2=((CCL-2)**0.75-1)/0.75. Linear mixed models (PROC mixed) are used to compare the difference in Box-Cox transformed CCL-2 levels among the three study groups, adjusting for age, gender (female vs male) and HLA genotype ([DR3/4 or DR4/4 or DR4/X] vs [DR3/3 or DR3/X or DRX/X]) and FDR status (FDR vs GP) with random intercept in the model47. Fit statistics between different curves over age were compared for a lower AIC/BIC and the model with quadratic curve over age was chosen as the final model. Additional statistical analyses were performed using the R language and environment for statistical computing (R version 4.0.5; R Foundation for Statistical Computing; [www.r-project.org](http://www.r-project.org)). PLINK1.948 was used for pre-processing the genotype data. LocusZoom49 was used to generate both Manhattan and LocusZoom plots. All p-values were two–tailed and a p< 0.05 after adjusting for multiple testing was considered statistically significant. ## Data Availability The data that support the findings from this study are openly available from Zenodo (https://zenodo.org/record/2615265#.YGs9GEhKjYk) and the GTEx portal (https://www.gtexportal.org/home/gene/CCR2). The T1DGC Immunochip summary statistics data and the CCL-2 serum DAISY data are available by request to the authors. ## Author contributions PMHT wrote the manuscript and researched data. SP acquired the Luminex data; SP, KBS, EK, PMHT, FD, KW and MR were responsible for data analysis. DH, KW, MR and JXS contributed to clinical samples. S O-G and SSR generated fine mapping data. All authors contributed to writing and editing of the manuscript. ## Acknowledgements Support for this work to JXS was provided by the National Institutes of Health ([R21HD050196](https://public.era.nih.gov/grantfolder/piAppDetails/genericStatus.do?encryptedParam=(v2)ETMsDgAAAWDgE1r-ABRBRVMvQ0JDL1BLQ1M1UGFkZGluZwCAABAAEMwuYtIFSDk0YfK-Nuphb90AAAAgAcuKccRAZX9Zf6uYMr6TpbwEibMk2IFnjXgZPZ6CbIsAFP1BxIWLl3WiLLbVVdQJJof8dBPG), R33HD050196, and 2RO1HD37800). SP was supported by Postdoctoral Fellowship (10-2006-792 and 3-2004-195) and Career Development Award (2-2011-153) from JDRF. PMHT was supported by NIH/NIDDK fellowship (F30DK121461). Diabetes Autoimmunity Study in the Young (MR, KW, and FD) is supported by the National Institutes of Health (R01 DK032493 and P30 DK116073) and Helmsley Charitable Trust (G-1901-03687). This research utilizes resources provided by the Type 1 Diabetes Genetics Consortium, a collaborative clinical study sponsored by the National Institute of Diabetes and Digestive and Kidney Diseases (NIDDK), National Institute of Allergy and Infectious Diseases (NIAID), National Human Genome Research Institute (NHGRI), National Institute of Child Health and Human Development (NICHD), and JDRF and supported by U01 [DK062418](https://www.ncbi.nlm.nih.gov/nuccore/DK062418). All authors declare no conflict of interest. * Received August 1, 2021. * Revision received August 1, 2021. * Accepted August 5, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. 1.Eisenbarth, G.S. Type I diabetes mellitus. A chronic autoimmune disease. N Engl J Med 314, 1360–1368 (1986). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJM198605223142106&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=3517648&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1986C377500006&link_type=ISI) 2. 2.DiMeglio, L.A., Evans-Molina, C. & Oram, R.A. Type 1 diabetes. Lancet 391, 2449–2462 (2018). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) 3. 3.Livingstone, S.J., et al. Estimated life expectancy in a Scottish cohort with type 1 diabetes, 2008-2010. Jama 313, 37–44 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2014.16425&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25562264&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) 4. 4.Herold, K.C., et al. An Anti-CD3 Antibody, Teplizumab, in Relatives at Risk for Type 1 Diabetes. New England Journal of Medicine 381, 603–613 (2019). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31180194&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) 5. 5.Barrett, J.C., et al. Genome-wide association study and meta-analysis find that over 40 loci affect risk of type 1 diabetes. Nat Genet 41, 703–707 (2009). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.381&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19430480&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000266411700019&link_type=ISI) 6. 6.Onengut-Gumuscu, S., et al. Fine mapping of type 1 diabetes susceptibility loci and evidence for colocalization of causal variants with lymphoid gene enhancers. Nature genetics 47, 381–386 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3245&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25751624&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) 7. 7.Li, Y.R., et al. Meta-analysis of shared genetic architecture across ten pediatric autoimmune diseases. Nat Med 21, 1018–1027 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nm.3933&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26301688&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) 8. 8.Godessart, N. & Kunkel, S.L. Chemokines in autoimmune disease. Curr Opin Immunol 13, 670–675 (2001). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0952-7915(01)00277-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11677088&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000171797600007&link_type=ISI) 9. 9.Kent, W.J., et al. The human genome browser at UCSC. Genome research 12, 996–1006 (2002). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ2Vub21lIjtzOjU6InJlc2lkIjtzOjg6IjEyLzYvOTk2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDgvMDUvMjAyMS4wOC4wMS4yMTI2MTQ0OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 10. 10.Bakos, E., et al. CCR2 Regulates the Immune Response by Modulating the Interconversion and Function of Effector and Regulatory T Cells. The Journal of Immunology 198, 4659–4671 (2017). 11. 11.Solomon, M., Balasa, B. & Sarvetnick, N. CCR2 and CCR5 chemokine receptors differentially influence the development of autoimmune diabetes in the NOD mouse. Autoimmunity 43, 156–163 (2010). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19824873&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) 12. 12.Purohit, S., et al. Large-Scale Discovery and Validation Studies Demonstrate Significant Reductions in Circulating Levels of IL8, IL-1Ra, MCP-1, and MIP-1β in Patients With Type 1 Diabetes. The Journal of Clinical Endocrinology & Metabolism 100, E1179–E1187 (2015). 13. 13.Grieb, G., Merk, M., Bernhagen, J. & Bucala, R. Macrophage migration inhibitory factor (MIF): a promising biomarker. Drug news & perspectives 23, 257 (2010). 14. 14.Ismail, N.A., et al. Monocyte chemoattractant protein 1 and macrophage migration inhibitory factor in children with type 1 diabetes. Journal of Pediatric Endocrinology and Metabolism 29, 641–645 (2016). 15. 15.Waugh, K., et al. Increased inflammation is associated with islet autoimmunity and type 1 diabetes in the Diabetes Autoimmunity Study in the Young (DAISY). PLoS One 12, e0174840 (2017). 16. 16.Panee, J. Monocyte Chemoattractant Protein 1 (MCP-1) in obesity and diabetes. Cytokine 60, 1–12 (2012). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cyto.2012.06.018&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22766373&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000310095000001&link_type=ISI) 17. 17.Piemonti, L., et al. Human Pancreatic Islets Produce and Secrete MCP-1/CCL2: Relevance in Human Islet Transplantation. Diabetes 51, 55–65 (2002). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZGlhYmV0ZXMiO3M6NToicmVzaWQiO3M6NzoiNTEvMS81NSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA4LzA1LzIwMjEuMDguMDEuMjEyNjE0NDkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 18. 18.Quinones, M.P., et al. The complex role of the chemokine receptor CCR2 in collagen-induced arthritis: implications for therapeutic targeting of CCR2 in rheumatoid arthritis. J Mol Med (Berl) 83, 672–681 (2005). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00109-005-0637-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15827759&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000231602200003&link_type=ISI) 19. 19.Koch, A.E., et al. Enhanced production of monocyte chemoattractant protein-1 in rheumatoid arthritis. J Clin Invest 90, 772–779 (1992). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1172/JCI115950&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1522232&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1992JN95900013&link_type=ISI) 20. 20.Haringman, J.J., Smeets, T.J., Reinders-Blankert, P. & Tak, P.P. Chemokine and chemokine receptor expression in paired peripheral blood mononuclear cells and synovial tissue of patients with rheumatoid arthritis, osteoarthritis, and reactive arthritis. Ann Rheum Dis 65, 294–300 (2006). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTE6ImFubnJoZXVtZGlzIjtzOjU6InJlc2lkIjtzOjg6IjY1LzMvMjk0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDgvMDUvMjAyMS4wOC4wMS4yMTI2MTQ0OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 21. 21.Bose, S. & Cho, J. Role of chemokine CCL2 and its receptor CCR2 in neurodegenerative diseases. Archives of pharmacal research 36, 1039–1050 (2013). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s12272-013-0161-z&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23771498&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) 22. 22.Novakova, L., et al. Cerebrospinal fluid biomarkers as a measure of disease activity and treatment efficacy in relapsing-remitting multiple sclerosis. Journal of neurochemistry 141, 296–304 (2017). 23. 23.Sindern, E., et al. Differential release of β-chemokines in serum and CSF of patients with relapsing–remitting multiple sclerosis. Acta neurologica scandinavica 104, 88–91 (2001). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1034/j.1600-0404.2001.104002088.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11493224&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000170274200006&link_type=ISI) 24. 24.Vergunst, C.E., et al. Modulation of CCR2 in rheumatoid arthritis: A double-blind, randomized, placebo-controlled clinical trial. Arthritis & Rheumatism 58, 1931–1939 (2008). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/art.23591&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18576354&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000257469800007&link_type=ISI) 25. 25.Cui, L.-Y., Chu, S.-F. & Chen, N.-H. The role of chemokines and chemokine receptors in multiple sclerosis. International Immunopharmacology 83, 106314 (2020). 26. 26.Gilbert, J., et al. Effect of CC chemokine receptor 2 CCR2 blockade on serum C-reactive protein in individuals at atherosclerotic risk and with a single nucleotide polymorphism of the monocyte chemoattractant protein-1 promoter region. The American journal of cardiology 107, 906–911 (2011). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.amjcard.2010.11.005&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21247529&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) 27. 27.Wang, Y., et al. CCR2 and CXCR4 regulate peripheral blood monocyte pharmacodynamics and link to efficacy in experimental autoimmune encephalomyelitis. J Inflamm (Lond) 6, 32 (2009). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1476-9255-6-32&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19906300&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) 28. 28.Zhao, B.N., et al. CCR2-Mediated Uptake of Constitutively Produced CCL2: A Mechanism for Regulating Chemokine Levels in the Blood. The Journal of Immunology 203, 3157–3165 (2019). 29. 29.Folkersen, L., et al. Genomic and drug target evaluation of 90 cardiovascular proteins in 30,931 individuals. Nat Metab 2, 1135–1148 (2020). 30. 30.Rewers, M., et al. Newborn screening for HLA markers associated with IDDM: diabetes autoimmunity study in the young (DAISY). Diabetologia 39, 807–812 (1996). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s001250050514&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8817105&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1996UV87800008&link_type=ISI) 31. 31.Norris, J.M., et al. Omega-3 polyunsaturated fatty acid intake and islet autoimmunity in children at increased risk for type 1 diabetes. Jama 298, 1420–1428 (2007). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.298.12.1420&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17895458&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) 32. 32.Buniello, A., et al. The NHGRI-EBI GWAS Catalog of published genome-wide association studies, targeted arrays and summary statistics 2019. Nucleic Acids Res 47, D1005–d1012 (2019). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gky1120&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30445434&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) 33. 33.Hillary, R.F., et al. Multi-method genome- and epigenome-wide studies of inflammatory protein levels in healthy older adults. Genome Med 12, 60 (2020). 34. 34.Sarwar, N., et al. IL6R Genetics Consortium Emerging Risk Factors Collaboration: Interleukin-6 receptor pathways in coronary heart disease: A collaborative meta-analysis of 82 studies. Lancet 379, 1205–1213 (2012). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(11)61931-4&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22421339&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000302230400033&link_type=ISI) 35. 35.Swerdlow, D.I., et al. The interleukin-6 receptor as a target for prevention of coronary heart disease: a mendelian randomisation analysis. Lancet 379, 1214–1224 (2012). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(12)60110-X&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22421340&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000302230400034&link_type=ISI) 36. 36.Hagopian, W.A., et al. The Environmental Determinants of Diabetes in the Young (TEDDY): genetic criteria and international diabetes risk screening of 421 000 infants. Pediatr Diabetes 12, 733–743 (2011). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1399-5448.2011.00774.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21564455&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000298170000010&link_type=ISI) 37. 37.Poma, A.M., et al. Immune Transcriptome of Cells Infected with Enterovirus Strains Obtained from Cases of Type 1 Diabetes. Microorganisms 8, 1031 (2020). 38. 38.Cortes, A. & Brown, M.A. Promise and pitfalls of the Immunochip. Arthritis Res Ther 13, 101–101 (2011). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/ar3382&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21345260&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) 39. 39.Lonsdale, J., et al. The Genotype-Tissue Expression (GTEx) project. Nature Genetics 45, 580–585 (2013). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2653&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23715323&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) 40. 40.Flores-Toro, J.A., et al. CCR2 inhibition reduces tumor myeloid cells and unmasks a checkpoint inhibitor effect to slow progression of resistant murine gliomas. Proceedings of the National Academy of Sciences 117, 1129–1138 (2020). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMDoiMTE3LzIvMTEyOSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA4LzA1LzIwMjEuMDguMDEuMjEyNjE0NDkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 41. 41.Linehan, D., et al. Overall survival in a trial of orally administered CCR2 inhibitor CCX872 in locally advanced/metastatic pancreatic cancer: Correlation with blood monocyte counts. Journal of Clinical Oncology 36, 92–92 (2018). 42. 42.Gale, J.D., et al. Effect of PF-04634817, an oral CCR2/5 chemokine receptor antagonist, on albuminuria in adults with overt diabetic nephropathy. Kidney international reports 3, 1316–1327 (2018). 43. 43.Vergunst, C.E., et al. Modulation of CCR2 in rheumatoid arthritis: a double-blind, randomized, placebo-controlled clinical trial. Arthritis & Rheumatism: Official Journal of the American College of Rheumatology 58, 1931–1939 (2008). 44. 44.Gilbert, J., et al. Effect of CC chemokine receptor 2 CCR2 blockade on serum C-reactive protein in individuals at atherosclerotic risk and with a single nucleotide polymorphism of the monocyte chemoattractant protein-1 promoter region. Am J Cardiol 107, 906–911 (2011). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.amjcard.2010.11.005&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21247529&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) 45. 45.Dayan, C.M., Korah, M., Tatovic, D., Bundy, B.N. & Herold, K.C. Changing the landscape for type 1 diabetes: the first step to prevention. Lancet 394, 1286–1296 (2019). 46. 46.Draper, N.R. & Smith, H. Applied regression analysis, (John Wiley & Sons, 1998). 47. 47.Laird, N.M. & Ware, J.H. Random-effects models for longitudinal data. Biometrics 38, 963–974 (1982). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2307/2529876&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7168798&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1982QF89700008&link_type=ISI) 48. 48.Purcell, S., et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81, 559–575 (2007). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/519795&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17701901&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) 49. 49.Pruim, R.J., et al. LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics 26, 2336–2337 (2010). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btq419&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20634204&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F05%2F2021.08.01.21261449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000281714100054&link_type=ISI)