Whole‐exome sequencing reveals the major genetic factors contributing to neuromyelitis optica spectrum disorder in Chinese patients with aquaporin 4‐IgG seropositivity

Neuromyelitis optica spectrum disorder (NMOSD) is an autoimmune disease. Although genetic factors are involved in its pathogenesis, limited evidence is available in this area. The aim of the present study was to identify the major genetic factors contributing to NMOSD in Chinese patients with aquaporin 4 (AQP4)-IgG seropositivity.Whole-exome sequencing (WES) was performed on 228 Chinese NMOSD patients seropositive for AQP4-IgG and 1400 healthy controls in Guangzhou, South China. Human leukocyte antigen (HLA) sequencing was also utilized. Genotype model and haplotype, gene burden, and enrichment analyses were conducted.A significant region of the HLA composition is on chromosome 6, and great variation was observed in DQB1, DQA2 and DQA1. HLA sequencing confirmed that the most significant allele was HLA-DQB1*05:02 (p < 0.01, odds ratio [OR] 3.73). The genotype model analysis revealed that HLA-DQB1*05:02 was significantly associated with NMOSD in the additive effect model and dominant effect model (p < 0.05). The proportion of haplotype "HLA-DQB1*05:02-DRB1*15:01" was significantly greater in the NMOSD patients than the controls, at 8.42% and 1.23%, respectively (p < 0.001, OR 7.39). The gene burden analysis demonstrated that loss-of-function mutations in NOP16 were more common in the NMOSD patients (11.84%) than the controls (5.71%; p < 0.001, OR 2.22). The IgG1-G390R variant was significantly more common in NMOSD, and the rate of the T allele was 0.605 in patients and 0.345 in the controls (p < 0.01, OR 2.92). The enrichment analysis indicated that most of the genetic factors were mainly correlated with nervous and immune processes.Human leukocyte antigen is highly correlated with NMOSD. NOP16 and IgG1-G390R play important roles in disease susceptibility.