Background: X chromosome-linked interleukin-1 receptor-associated kinase (IRAK1) polymorphisms have been demonstrated to be associated with the risks of several autoimmune diseases, such as systemic lupus erythematosus, systemic sclerosis, rheumatoid arthritis, and autoimmune thyroid diseases. However, no studies have investigated the association of IRAK1 polymorphisms with neuromyelitis optica spectrum disorder (NMOSD). This case-control study was performed to determine the correlation between IRAK1 polymorphisms and the risk of NMOSD. Methods: Two single nucleotide polymorphisms (SNPs) rs1059703G > A and rs3027898C > A of IRAK1 were selected and genotyped using SNPscan in a Chinese cohort, including 332 patients with NMOSD and 520 healthy controls. Chi-square tests and logistic regression analyses were used to determine the associations between IRAK1 polymorphisms and the risk of NMOSD. Results: Patients with NMOSD showed a lower frequency of the minor allele A of rs1059703 than did controls (Odds ratio [OR] = 0.68; 95% confidence intervals [CI], 0.52-0.88; P-co(rr) = 0.007). Compared with wild genotype GG of rs1059703, homozygous mutation AA and heterozygous mutation GA were significantly associated with the decreased risk of NMOSD after adjusting for sex and age (adjusted OR = 0.64; 95%CI, 0.49-0.84; P-co(rr) = 0.002). Similar associations were also observed for IRAK1 rs3027898C > A. Stratification analysis according to sex revealed that the significantly different allele distributions of the two SNPs were mainly found in females. However, IRAK1 polymorphisms were not correlated with aquaporin-4-IgG, onset symptoms, or age at onset. Conclusions: This study is first to demonstrate that X-chromosome-linked IRAK1 polymorphisms are associated with the risk of NMOSD and provide novel insights into the underlying mechanisms of this disease. Further studies are needed to elucidate the function of IRAK1 variants in the pathogenesis of NMOSD and the underlying molecular mechanisms.