- 1Department of Preventive Medicine, Institute of Nutrition and Diseases, Wenzhou Medical University, Wenzhou, China
- 2Department of Neurology, Institute of Geriatric Neurology, the Second Affiliated Hospital and Yuying Children’s Hospital, Wenzhou Medical University, Wenzhou, China
CLEC16A is a membrane-associated endosomal protein implicated in regulating autophagy and antigen presentation. Its genetic variants are broadly associated with multiple autoimmune diseases. Parkinson’s disease (PD), which undergoes autophagy disruption and neuroinflammation, has been clinically observed, for an extensive amount of time, to be associated with autoimmune diseases. In this study, we aimed to understand whether the autoimmune disease associated CLEC16A variants pleiotropically modulate PD risk. Five of such CLEC16A variants, including rs6498169, rs12708716, rs12917716, rs7200786, and rs2903692, were selected and analyzed in a Han Chinese cohort comprising 515 sporadic PD patients and 504 controls. Results showed that rs6498169 and rs7200786 were significantly associated with PD susceptibility (p = 0.005 and 0.004, respectively; recessive model, p = 0.002 and 0.001, respectively). Rs6498169 was also associated with the PD subtype of postural instability/gait difficulty (p = 0.002). Haplotype analysis showed that the AAG module in order of rs6498169, rs12708716, and rs2903692 was associated with the highest risk for PD (p = 0.0047, OR = 1.42, 95% CI = 1.11–1.82). Functional annotation analyses suggested that rs6498169 had high probability to affect transcription factor binding and target gene expression. In summary, the current study demonstrates that the autoimmune disease associated CLEC16A variants convey risk of PD in Han Chinese. Our findings suggest a pleiotropic role of CLEC16A and strengthen the link between PD and autoimmune diseases.
Introduction
CLEC16A (C-type lectin domain family 16, member A; Chr16; Gene ID: 23274; MIM: 611303), encoding a large protein of 1,053 amino acids with a putative C-type lectin-like domain, is a membrane-associated endosomal protein (Soleimanpour et al., 2014). CLEC16A genetic variations have been associated with multiple autoimmune diseases including multiple sclerosis, type 1 diabetes, Crohn’s disease, Addison’s disease, rheumatoid arthritis, and juvenile idiopathic arthritis (Vitale et al., 2002; Hakonarson et al., 2007; International Multiple Sclerosis Genetics et al., 2007; Marquez et al., 2009; Martinez et al., 2010; Skinningsrud et al., 2010; Nischwitz et al., 2011; International Multiple Sclerosis Genetics et al., 2013). The broad association with autoimmune disorders suggests that CLEC16A may functionally link to autoimmunity by certain common pathogenic pathways. Indeed, CLEC16A has been implicated in regulating autophagy/mitophagy and antigen presentation. CLEC16A interacts with E3 ubiquitin ligase NRDP1 (neuregulin receptor degradation protein 1), controlling the volume of Parkin and its mastered mitophagy and thereby mediates murine β cell function and diabetogenesis (Soleimanpour et al., 2014). In addition, Clec16a knockdown in the nonobese diabetic mice protects against autoimmunity by altering the T cell selection, possibly through the inhibition of thymic epithelial cell autophagy (Schuster et al., 2015). Clec16a mutant mice generated by gene-trap insertion develop neurodegeneration featured by impaired motor behaviors, loss of Purkinje cells, and abnormal autophagy homeostasis (Redmann et al., 2016). CLEC16A also plays a role in multiple sclerosis via controlling the late endosome biogenesis-mediated human leukocyte antigen class II antigen presentation (van Luijn et al., 2015).
Parkinson’s disease (PD) is the second most common neurodegenerative disorder characterized by the loss of dopaminergic neurons in the substantia nigra pars compacta. The majority of PD is sporadic and its etiology is complex and remains largely unknown. Autophagy dysfunction is one of the key mechanisms involved in PD, leading to the abnormal degradation and aggregation of proteins such as α-synuclein (Michel et al., 2016). Mutations in PRKN (Parkin encoding gene) are the common cause of the early-onset PD (Kalia and Lang, 2015). Recent studies have also revealed the presence of α-synuclein-reactive T cells in PD patients and an association of its reactivity with preclinical and early PD, indicating a link of PD to autoimmunity (Sulzer et al., 2017; Lindestam Arlehamn et al., 2020).
Indeed, PD risk has been observed, for an extensive amount of time, as associated with autoimmune diseases. For instance, a Swedish epidemiological study showed that subsequent risks of PD are increased in patients with 6 autoimmune diseases, including amyotrophic lateral sclerosis, Graves’s disease/hyperthyroidism, Hashimoto’s disease/hypothyroidism, multiple sclerosis, pernicious anemia, and polymyalgia rheumatica (Li et al., 2012), whereas being a casual or causal association between PD risk and multiple sclerosis remains in debate (Pedemonte et al., 2013; Nielsen et al., 2014). Nonetheless, it is hypothesized that there are common genetic risk variants between PD and autoimmune diseases (Witoelar et al., 2017). Given the functional overlap of CLEC16A with PD pathogenic mechanisms and its broad association with autoimmune diseases, we aimed, in this study, to investigate whether the autoimmune disease associated CLEC16A variants convey risk of PD in a Han Chinese population.
Materials and Methods
Subjects
A total of 1,019 subjects of Han Chinese ethnicity from eastern China were recruited in this study, comprising 515 sporadic PD patients (263 males and 252 females) and 504 controls (265 males and 239 females). The median age of the patients and controls was 66 (interquartile range: 60–73) and 66 (interquartile range: 57–72) years old, respectively. All PD patients were diagnosed by two movement disorder neurologists, according to the UK Parkinson’s disease Society Brain Bank Criteria (Hughes et al., 1992). Patients with a family history of PD or with secondary and atypical parkinsonism were excluded. The controls had no neurological disorders according to medical history, physical, and laboratory examinations. All subjects participating in the study signed written informed consents. The study was performed under the approval No. LCKY 2020-66 by the Ethics Committee of The Second Affiliated Hospital and Yuying Children’s Hospital, Wenzhou Medical University.
Selection of CLEC16A Risk Variants
Five autoimmune disease-associated risk variants of CLEC16A were selected, including rs6498169, rs12708716, rs7200786, rs2903692 and rs12917716. Among these, rs6498169, rs12708716 and rs7200786 were identified as susceptibility markers for multiple sclerosis (Vitale et al., 2002; International Multiple Sclerosis Genetics et al., 2007; International Multiple Sclerosis Genetics et al., 2013); rs2903692 and rs12917716 were found to be associated with type 1 diabetes (Hakonarson et al., 2007; Skinningsrud et al., 2010); rs6498169 was additionally associated with juvenile idiopathic arthritis and rheumatoid arthritis (Martinez et al., 2010; Skinningsrud et al., 2010); rs2903692 was additionally associated with Crohn’s disease and multiple sclerosis (Marquez et al., 2009; Martinez et al., 2010); rs12917716 was additionally associated with Addison’s disease, multiple sclerosis, and primary adrenal insufficiency (Skinningsrud et al., 2010; Nischwitz et al., 2011).
Genotyping
Genomic DNA was extracted from peripheral blood using TIANamp Genomic DNA kit (Tiangen, Beijing, China) according to the manufacturer’s instruction. Four SNPs, including rs6498169, rs12708716, rs12917716, and rs7200786, were genotyped using SNaPshot at BGI Technology (Wuhan, China) as described previously (Zou et al., 2018). Rs2903692 was not successfully genotyped by SNaPshot, and therefore was genotyped by conventional PCR and sequencing. The PCR condition was initial denaturation at 95°C for 3 min, followed by 32 cycles of 95°C for 30 s, 60°C for 30 s, and 72°C for 60 s, and a final extension at 72°C for 5 min. The PCR products were sequenced at BGI Technology (Wuhan, China). All relevant primers were listed in Supplementary Table S1.
Statistical Analysis
Statistical analyses were carried out using the Statistical Package for Social Science program (SPSS for Windows, version 23.0). Hardy-Weinberg equilibrium in genotype distribution was assessed using χ2 test. Following Kolmogorov-Smirnov test for normality, Mann-Whitney U test was used to evaluate age difference. The χ2 test was also used to assess differences in gender, genotype, allele, and haplotype frequencies between the PD patients and controls. The haplotype construction and association analysis were performed using SNPStats Online Version (https://www.snpstats.net/start.htm) (Sole et al., 2006). A backward elimination method as reported previously was used to identify the highest-risk haplotype for PD (Francis et al., 2007). Statistical power was calculated by the QUANTO version 1.2.4. A p value < 0.05 was considered statistically significant unless otherwise indicated.
Function Annotation
Function annotations of the variants were obtained from HaploReg v4.1 (http://pubs.broadinstitute.org/mammals/haploreg/haploreg.php) and RegulomeDB (http://www.regulomedb.org/). HaploReg was used to annotate potential causal links to disease pathogenesis (Ward and Kellis, 2016). RegulomeDB was used to annotate known and predicted genetic variations in regulatory elements in intergenic regions of the human genome. A score ranging from 1 to 6 was provided to indicate the potential function and the lower score indicates a higher probability that a variant affects binding and gene expression (Boyle et al., 2012). The expression quantitative trait locus (eQTL) and splicing quantitative trait locus (sQTL) of the variants were analyzed by the GTEx Portal (https://gtexportal.org/).
Results
Association of the CLEC16A Variants with PD Susceptibility
Genotype distributions of the 5 variants in controls were in accordance with Hardy-Weinberg equilibrium (p > 0.05). The PD cases and controls were comparable in both gender and age (p > 0.05). Difference was considered only after Bonferroni correction (threshold for significance = 0.01). Results showed that significant difference between PD cases and control cases was found only in genotypes of rs6498169 and rs7200786 (p = 0.005 and p = 0.004, respectively; Table 1). No difference was found within rs12708716, rs12917716, and rs2903692.
We further analyzed the rs6498169 and rs7200786 by three genetic models (additive, dominant, and recessive). Both SNPs were significantly associated with PD in the recessive model (p = 0.002, OR = 1.629, 95% CI = 1.203–2.205 for rs6498169; p = 0.001, OR = 1.895, 95% CI = 1.285–2.795 for rs7200786). The risk genotypes for PD were the AA of rs6498169 and the GG of rs7200786 (Table 2). The statistical power was 92% and 96%, respectively, for the recessive model of rs6498169 and rs7200786.
We also extracted a total of 265 PD patients recorded with PD subtypes to understand their association with the five CLEC16A SNPs. The subtypes were classified into postural instability/gait difficulty (PIGD), tremor dominant, and indeterminate. Results showed that the genotype distribution of rs6498169 was significantly different (p = 0.002) between the PIGD patients and the controls. No other difference was found in the variants between the PD subtypes and the controls (Supplementary Table S2).
Haplotype Analyses of the CLEC16A Variants
We analyzed whether haplotypes of these CLEC16A variants were associated with PD. Haplotypes were constructed in the following order: rs6498169, rs12708716, rs12917716, rs7200786, and rs2903692. As shown in Table 3, five haplotypes were listed as those with a frequency < 3.0% in both PD patients and controls were excluded. A significant difference (p = 0.034) in overall haplotype distribution was observed between the PD patients and controls. Two haplotypes, AAGAG (p = 0.02, OR = 1.47, 95% CI = 1.06–2.04) and AACGG (p = 0.03, OR = 1.47, 95% CI = 1.04–2.07), showed a significant difference between the cases and controls.
To characterize the highest-risk haplotype of CLEC16A towards PD, a backward elimination model was employed. Results showed that the strongest PD-associated haplotype was AAG in the order of rs6498169, rs12708716, and rs2903692 from the best 3-variant model (p = 0.0047, OR = 1.42, 95% CI = 1.11–1.82; Table 4).
Function Annotations of the PD-Associated CLEC16A Variants
Function annotations of the PD-associated variants were performed using the HaploReg, RegulomeDB, and GTEx. Based on the HaploRreg (Tables 5 and Supplementary Table S3), rs6498169 was located within enhancer histone marks in 13 tissues, DNase hypersensitivity regions in 7 tissues, and the region of GATA2 (GATA binding protein 2) binding site, and was predicted to significantly alter the binding motif of the HES1 (hes family bHLH transcription factor 1) transcription factor. Rs7200786 was located within enhancer histone marks in 13 tissues, and DNase hypersensitivity regions in 3 tissues and was predicted to significantly change the binding motifs of the CCNT2 (cyclin T2) and EBF (early B cell factor) transcription factors.
In RegulomeDB (Table 5), where the score <3 indicates a relatively high possibility of potential regulatory function (Luciano et al., 2018), rs6498169 was predicted with a score of 1b, representing that this variant is highly likely to affect transcription factor binding, certain motifs, DNase footprint and peaks, and potentially affect the expression of target genes. The score of rs7200786 was 4, representing that this variant may affect transcription factor binding and DNase peaks.
By using GTEx (Table 5; Supplementary Figures S1–S3), rs6498169 was suggested to be significant eQTL of the RMI2 (RecQ mediated genome instability 2) expression in brain tissues such as cortex (p = 0.0015), cerebellum (p = 0.008) and nucleus accumbens (p = 0.008). Rs7200786 was as significant eQTL of the RMI2 expression in brain tissues such as nucleus accumbens (p = 0.00071), cortex (p = 0.0028) and putamen (p = 0.0033), as well as the CLEC16A expression in brain tissues such as cerebellar hemisphere (p = 0.0012) and amygdala (p = 0.008). However, both rs6498169 and rs7200786 were not associated with the expression of RMI2 and CLEC16A in the substantia nigra. Results of the sQTL analysis showed that the risk alleles of both rs6498169 and rs7200786 (A and G, respectively) were significantly associated with an increased intron-excision ratio, which potentially leads to elevation of certain variant expression of CLEC16A (Supplementary Figure S4).
Discussion
Immunity disturbance is increasingly considered to be important in PD pathogenesis. The clinical observations of risk association between PD and autoimmune diseases indicate the possibility of having common genetic risk loci in between. By studying the five autoimmune disease associated CLEC16A variants in a Chinese cohort, we demonstrate that CLEC16A is pleiotropic for modulating PD risk. Two loci, rs6498169 and rs7200786, are recessively associated with PD susceptibility.
CLEC16A is involved in regulation of autophagy, T cell selection, antigen presentation, and neurodegeneration as suggested earlier (Soleimanpour et al., 2014; Schuster et al., 2015; van Luijn et al., 2015; Redmann et al., 2016). Autophagy plays a role, not only for substance clearance and recycling, but also in the presentation of antigenic peptides to the receptor of T cells in the context of antigen-presenting cells and major histocompatibility complex class II (Bonam and Muller, 2020). While disruption of the autophagy pathway serves as one of the key mechanisms in PD (Michel et al., 2016), PD patients are indeed observed with dysregulated innate and adaptive immune responses, particularly in those carrying autophagy-related gene mutations such as in Parkin and LRRK2 (leucine rich repeat kinase 2). Besides the previously mentioned α-synuclein reactive T cell responses (Sulzer et al., 2017; Lindestam Arlehamn et al., 2020), B cell-produced autoantibodies against α-synuclein antigen, GM1 gangliosides, and neuronal antigens are also shown to be elevated in blood and/or cerebral spinal fluid of PD patients (Zappia et al., 2002; van de Warrenburg et al., 2008; Scott et al., 2018). Hence, besides the identified genetic connection, CLEC16A may be functionally possible to participate in PD.
Our results show that rs6498169 and rs7200786 are associated with PD susceptibility. The modulation of rs7200786 on PD risk is consistent with a previous report in the Italian population (Strafella et al., 2021). The A of rs6498169 and the G of rs7200786 are recessive risk alleles. However, these two alleles are protective against multiple sclerosis and rheumatoid arthritis (Vitale et al., 2002; International Multiple Sclerosis Genetics et al., 2007; Martinez et al., 2010; Skinningsrud et al., 2010; International Multiple Sclerosis Genetics et al., 2013). When the relationship between PD and multiple sclerosis remains in debate (Li et al., 2012; Pedemonte et al., 2013; Nielsen et al., 2014), PD risk, indeed, appears to be negatively correlated with rheumatoid arthritis (Bacelis et al., 2021; Li et al., 2021). These results suggest that the risk variants of CLEC16A may have differential actual impact on PD and certain autoimmune diseases. Interestingly, a recent GWAS pooling study of European ancestry identified 17 shared susceptibility loci between autoimmune diseases and PD (Witoelar et al., 2017). Similar to our case, these loci have not been previously reported by PD GWAS studies (Obeso et al., 2017), but partially reported by individual polymorphism studies. Thus, a variety in research designs and ethnicities may still be valuable in searching for new genetic risks. Conversely, these findings may need further validation in additional populations. Being noted, the above European study identified none of the CLEC16A variants (Witoelar et al., 2017). To understand the discordance, we examined allele frequencies of the five variants in Europeans and East Asians in the gnomAD (https://gnomad.broadinstitute.org/). Interestingly, the two PD-associated CLEC16A variants happen to be with inverted major and minor allele distributions in these two ethnicities. In detail, the G allele frequency of rs6498169 is at 0.358 and 0.554, and the A allele frequency of rs7200786 is at 0.461 and 0.660, respectively, in Europeans and East Asians (Supplementary Table S4). In this case, the pleiotropic discordance may partially attribute to the ethnicity-associated evolutional divergency in CLEC16A genetic variation.
It is known that haplotypes are more powerful for the detection of susceptibility alleles than individual variants (Gabriel et al., 2002). By analyzing the five variants, we identified two haplotypes, AAGAG and AACGG (in the order of rs6498169, rs12708716, rs12917716, rs7200786, and rs2903692), serving as risk factor. Further analysis of the effect polymorphisms suggests that the AAG of rs6498169, rs12708716 and rs2903692 represents the core haplotype associated with PD. Results of the functional annotation analyses suggest that rs6498169 is highly probable to affect transcription factor binding and target gene expression, such as through the modulation of GATA2 and HES1 binding. In contrast, rs7200786 is with relatively less probability to be functional. These results appear to be in line with the above highest-risk haplotype analysis, which shows that the rs6498169 locus, but not rs7200786, is within the core in association with PD. Based on the quantitative trait locus analysis, the CLEC16A expression may be affected by rs6498169 in a splicing-regulating way and by rs7200786 in both expression- and splicing-regulating ways. In addition, rs6498169 and rs7200786 may also affect the RMI2 expression in an expression-regulating way. RMI2 is a component of the BLM (Bloom syndrome RecQ like helicase) complex and is essential for genome stability (Hudson et al., 2016). However, the effects of these two variants on gene expression appear mainly in brain regions other than the substantia nigra. Indeed, although the substantia nigra is most profoundly affected in PD, other brain regions are also important for PD pathogenesis and clinical manifestations, such as putamen, cortex and amygdala (Kalia and Lang, 2015).
In conclusion, the present study demonstrates that the autoimmune disease associated CLEC16A genetic variants are associated with PD susceptibility in Han Chinese. Specifically, the A of rs6498169 and the G of rs7200786 serve as recessive risk alleles towards PD. These findings provide genetic insights into the pleiotropic role of CLEC16A and strengthen the link between PD and autoimmune diseases.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The name of the repository and link can be found below: Dryad; https://datadryad.org/stash/share/_yc2e6e8jRgT3dCBI-WN2w4NZNY1CqRkqiSnD4q1h7w.
Ethics Statement
The studies involving human participants were reviewed and approved by the Second Affiliated Hospital and Yuying Children’s Hospital, Wenzhou Medical University. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
H-HF and LC designed the research; H-HF, LC, X-XJ, S-SL, K-YW, H-JD, and MM conducted experiments; LC, X-XJ, Y-DS, and H-MW analyzed data; XZ contributed samples; BO helped data analysis and edited the manuscript; LC and J-HZ wrote the manuscript; H-MW, J-HZ, and XZ supervised the study. All the authors have read and approved the final manuscript.
Funding
This work was supported by the funding from Zhejiang Provincial Natural Science Foundation (LD19H090001 and LZ19H090002), National Natural Science Foundation of China (81771380, 81771510, and 82071585), Zhejiang Provincial Medical Technology Program (2019KY445), and Wenzhou Medical University (89218014).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors are thankful to their colleagues for clinical data collection and management of patients. The authors appreciate all the subjects participating in this study.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.856493/full#supplementary-material
References
Bacelis, J., Compagno, M., George, S., Pospisilik, J. A., Brundin, P., Naluai, Å. T., et al. (2021). Decreased Risk of Parkinson's Disease after Rheumatoid Arthritis Diagnosis: A Nested Case-Control Study with Matched Cases and Controls. Jpd 11 (2), 821–832. doi:10.3233/JPD-202418
Bonam, S. R., and Muller, S. (2020). Parkinson's Disease Is an Autoimmune Disease: A Reappraisal. Autoimmun. Rev. 19 (12), 102684. doi:10.1016/j.autrev.2020.102684
Boyle, A. P., Hong, E. L., Hariharan, M., Cheng, Y., Schaub, M. A., Kasowski, M., et al. (2012). Annotation of Functional Variation in Personal Genomes Using RegulomeDB. Genome Res. 22 (9), 1790–1797. doi:10.1101/gr.137323.112
Francis, P. J., Schultz, D. W., Hamon, S., Ott, J., Weleber, R. G., and Klein, M. L. (2007). Haplotypes in the Complement Factor H (CFH) Gene: Associations with Drusen and Advanced Age-Related Macular Degeneration. PLoS One 2 (11), e1197. doi:10.1371/journal.pone.0001197
Gabriel, S. B., Schaffner, S. F., Nguyen, H., Moore, J. M., Roy, J., Blumenstiel, B., et al. (2002). The Structure of Haplotype Blocks in the Human Genome. Science 296 (5576), 2225–2229. doi:10.1126/science.1069424
Hakonarson, H., Grant, S. F. A., Bradfield, J. P., Marchand, L., Kim, C. E., Glessner, J. T., et al. (2007). A Genome-wide Association Study Identifies KIAA0350 as a Type 1 Diabetes Gene. Nature 448 (7153), 591–594. doi:10.1038/nature06010
Hudson, D. F., Amor, D. J., Boys, A., Butler, K., Williams, L., Zhang, T., et al. (2016). Loss of RMI2 Increases Genome Instability and Causes a Bloom-like Syndrome. Plos Genet. 12 (12), e1006483. doi:10.1371/journal.pgen.1006483
Hughes, A. J., Daniel, S. E., Kilford, L., and Lees, A. J. (1992). Accuracy of Clinical Diagnosis of Idiopathic Parkinson's Disease: a Clinico-Pathological Study of 100 Cases. J. Neurol. Neurosurg. Psychiatry 55 (3), 181–184. doi:10.1136/jnnp.55.3.181
International Multiple Sclerosis Genetics, C., Beecham, A. H., Patsopoulos, N. A., Xifara, D. K., Davis, M. F., Kemppinen, A., et al. (2013). Analysis of Immune-Related Loci Identifies 48 New Susceptibility Variants for Multiple Sclerosis. Nat. Genet. 45 (11), 1353–1360. doi:10.1038/ng.2770
International Multiple Sclerosis Genetics, C., Hafler, D. A., Compston, A., Sawcer, S., Lander, E. S., Daly, M. J., et al. (2007). Risk Alleles for Multiple Sclerosis Identified by a Genomewide Study. N. Engl. J. Med. 357 (9), 851–862. doi:10.1056/NEJMoa073493
Kalia, L. V., and Lang, A. E. (2015). Parkinson's Disease. The Lancet 386 (9996), 896–912. doi:10.1016/S0140-6736(14)61393-3
Li, C., Ou, R., and Shang, H. (2021). Rheumatoid Arthritis Decreases Risk for Parkinson's Disease: a Mendelian Randomization Study. Npj Parkinsons Dis. 7 (1), 17. doi:10.1038/s41531-021-00166-x
Li, X., Sundquist, J., and Sundquist, K. (2012). Subsequent Risks of Parkinson Disease in Patients with Autoimmune and Related Disorders: a Nationwide Epidemiological Study from Sweden. Neurodegener. Dis. 10 (1-4), 277–284. doi:10.1159/000333222
Lindestam Arlehamn, C. S., Dhanwani, R., Pham, J., Kuan, R., Frazier, A., Rezende Dutra, J., et al. (2020). α-Synuclein-specific T Cell Reactivity Is Associated with Preclinical and Early Parkinson's Disease. Nat. Commun. 11 (1), 1875. doi:10.1038/s41467-020-15626-w
Luciano, M., Hagenaars, S. P., Davies, G., Hill, W. D., Clarke, T.-K., Shirali, M., et al. (2018). Association Analysis in over 329,000 Individuals Identifies 116 Independent Variants Influencing Neuroticism. Nat. Genet. 50 (1), 6–11. doi:10.1038/s41588-017-0013-8
Márquez, A., Varadé, J., Robledo, G., Martínez, A., Mendoza, J. L., Taxonera, C., et al. (2009). Specific Association of a CLEC16A/KIAA0350 Polymorphism with NOD2/CARD15− Crohn's Disease Patients. Eur. J. Hum. Genet. 17 (10), 1304–1308. doi:10.1038/ejhg.2009.50
Martínez, A., Perdigones, N., Cénit, M. C., Espino, L., Varadé, J., Lamas, J. R., et al. (2010). Chromosomal Region 16p13: Further Evidence of Increased Predisposition to Immune Diseases. Ann. Rheum. Dis. 69 (1), 309–311. doi:10.1136/ard.2008.098376
Michel, P. P., Hirsch, E. C., and Hunot, S. (2016). Understanding Dopaminergic Cell Death Pathways in Parkinson Disease. Neuron 90 (4), 675–691. doi:10.1016/j.neuron.2016.03.038
Nielsen, N. M., Pasternak, B., Stenager, E., Koch-Henriksen, N., and Frisch, M. (2014). Multiple Sclerosis and Risk of Parkinson's Disease: a Danish Nationwide Cohort Study. Eur. J. Neurol. 21 (1), 107–111. doi:10.1111/ene.12255
Nischwitz, S., Cepok, S., Kroner, A., Wolf, C., Knop, M., Müller-Sarnowski, F., et al. (2011). More CLEC16A Gene Variants Associated with Multiple Sclerosis. Acta Neurol. Scand. 123 (6), 400–406. doi:10.1111/j.1600-0404.2010.01421.x
Obeso, J. A., Stamelou, M., Goetz, C. G., Poewe, W., Lang, A. E., Weintraub, D., et al. (2017). Past, Present, and Future of Parkinson's Disease: A Special Essay on the 200th Anniversary of the Shaking Palsy. Mov Disord. 32 (9), 1264–1310. doi:10.1002/mds.27115
Pedemonte, E., Trabucco, E., Cella, M., and Solaro, C. (2013). Parkinsonism in Multiple Sclerosis Patients: a Casual or Causal Association? Parkinsonism Relat. Disord. 19 (4), 492–493. doi:10.1016/j.parkreldis.2012.11.021
Redmann, V., Lamb, C. A., Hwang, S., Orchard, R. C., Kim, S., Razi, M., et al. (2016). Clec16a Is Critical for Autolysosome Function and Purkinje Cell Survival. Sci. Rep. 6, 23326. doi:10.1038/srep23326
Schuster, C., Gerold, K. D., Schober, K., Probst, L., Boerner, K., Kim, M.-J., et al. (2015). The Autoimmunity-Associated Gene CLEC16A Modulates Thymic Epithelial Cell Autophagy and Alters T Cell Selection. Immunity 42 (5), 942–952. doi:10.1016/j.immuni.2015.04.011
Scott, K. M., Kouli, A., Yeoh, S. L., Clatworthy, M. R., and Williams-Gray, C. H. (2018). A Systematic Review and Meta-Analysis of Alpha Synuclein Auto-Antibodies in Parkinson's Disease. Front. Neurol. 9, 815. doi:10.3389/fneur.2018.00815
Skinningsrud, B., Lie, B. A., Husebye, E. S., Kvien, T. K., Forre, O., Flato, B., et al. (2010). A CLEC16A Variant Confers Risk for Juvenile Idiopathic Arthritis and Anti-cyclic Citrullinated Peptide Antibody Negative Rheumatoid Arthritis. Ann. Rheum. Dis. 69 (8), 1471–1474. doi:10.1136/ard.2009.114934
Sole, X., Guino, E., Valls, J., Iniesta, R., and Moreno, V. (2006). SNPStats: a Web Tool for the Analysis of Association Studies. Bioinformatics 22 (15), 1928–1929. doi:10.1093/bioinformatics/btl268
Soleimanpour, S. A., Gupta, A., Bakay, M., Ferrari, A. M., Groff, D. N., Fadista, J., et al. (2014). The Diabetes Susceptibility Gene Clec16a Regulates Mitophagy. Cell 157 (7), 1577–1590. doi:10.1016/j.cell.2014.05.016
Strafella, C., Caputo, V., Termine, A., Assogna, F., Pellicano, C., Pontieri, F. E., et al. (2021). Immune System and Neuroinflammation in Idiopathic Parkinson's Disease: Association Analysis of Genetic Variants and miRNAs Interactions. Front. Genet. 12, 651971. doi:10.3389/fgene.2021.651971
Sulzer, D., Alcalay, R. N., Garretti, F., Cote, L., Kanter, E., Agin-Liebes, J., et al. (2017). T Cells from Patients with Parkinson's Disease Recognize α-synuclein Peptides. Nature 546 (7660), 656–661. doi:10.1038/nature22815
van de Warrenburg, B. P. C., Church, A. J., Martino, D., Candler, P. M., Bhatia, K. P., Giovannoni, G., et al. (2008). Antineuronal Antibodies in Parkinson's Disease. Mov. Disord. 23 (7), 958–963. doi:10.1002/mds.21929
van Luijn, M. M., Kreft, K. L., Jongsma, M. L., Mes, S. W., Wierenga-Wolf, A. F., van Meurs, M., et al. (2015). Multiple Sclerosis-Associated CLEC16A Controls HLA Class II Expression via Late Endosome Biogenesis. Brain 138 (Pt 6), 1531–1547. doi:10.1093/brain/awv080
Vitale, E., Cook, S., Sun, R., Specchia, C., Subramanian, K., Rocchi, M., et al. (2002). Linkage Analysis Conditional on HLA Status in a Large North American Pedigree Supports the Presence of a Multiple Sclerosis Susceptibility Locus on Chromosome 12p12. Hum. Mol. Genet. 11 (3), 295–300. doi:10.1093/hmg/11.3.295
Ward, L. D., and Kellis, M. (2016). HaploReg V4: Systematic Mining of Putative Causal Variants, Cell Types, Regulators and Target Genes for Human Complex Traits and Disease. Nucleic Acids Res. 44 (D1), D877–D881. doi:10.1093/nar/gkv1340
Witoelar, A., Jansen, I. E., Wang, Y., Desikan, R. S., Gibbs, J. R., Blauwendraat, C., et al. (2017). Genome-wide Pleiotropy between Parkinson Disease and Autoimmune Diseases. JAMA Neurol. 74 (7), 780–792. doi:10.1001/jamaneurol.2017.0469
Zappia, M., Crescibene, L., Bosco, D., Arabia, G., Nicoletti, G., Bagalà, A., et al. (2002). Anti-GM1 Ganglioside Antibodies in Parkinson's Disease. Acta Neurol. Scand. 106 (1), 54–57. doi:10.1034/j.1600-0404.2002.01240.x
Keywords: CLEC16A, Parkinson’s disease, autoimmune diseases, genetic variation, association
Citation: Fan H-H, Cui L, Jiang X-X, Song Y-D, Liu S-S, Wu K-Y, Dong H-J, Mao M, Ovlyakulov B, Wu H-M, Zhu J-H and Zhang X (2022) Autoimmune Disease Associated CLEC16A Variants Convey Risk of Parkinson’s Disease in Han Chinese. Front. Genet. 13:856493. doi: 10.3389/fgene.2022.856493
Received: 17 January 2022; Accepted: 09 February 2022;
Published: 30 March 2022.
Edited by:
Corrado Romano, University of Catania, ItalyReviewed by:
Félix Javier Jiménez-Jiménez, Hospital Universitario del Sureste, SpainPaul Lasko, McGill University, Canada
Copyright © 2022 Fan, Cui, Jiang, Song, Liu, Wu, Dong, Mao, Ovlyakulov, Wu, Zhu and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hong-Mei Wu, aG13dUB3bXUuZWR1LmNu; Jian-Hong Zhu, amh6aHVAd211LmVkdS5jbg==; Xiong Zhang, emhhbmd4aW9uZzk4QGdtYWlsLmNvbQ==
†These authors have contributed equally to this work