Immunocyte phenotypes underlying the causal autoimmune features in narcolepsy type 1
Original Article

Immunocyte phenotypes underlying the causal autoimmune features in narcolepsy type 1

Tianlong Li1#, Yabang Chen1#, Pingan Zhang2, Qingqian Zhu1, Zhuoshen Lin1, Baoxin Peng1, Rundong Qin2, Nicolas Steenbergen3, Rulong Hu1, Hua Qin1,4, Xiaowen Zhang1

1State Key Laboratory of Respiratory Disease, Department of Otolaryngology, Head & Neck Surgery, The First Affiliated Hospital of Guangzhou Medical University, Guangzhou, China; 2State Key Laboratory of Respiratory Disease, National Clinical Research Center for Respiratory Disease, Guangzhou Institute of Respiratory Health;Department of Allergy and Clinical Immunology, The First Affiliated Hospital of Guangzhou Medical University, Guangzhou, China; 3Department of Bioengineering, Imperial College London, London, UK; 4Interdisciplinary Center of Sleep Medicine, Charité-Universitätsmedizin Berlin, Berlin, Germany

Contributions: (I) Conception and design: X Zhang, H Qin; (II) Administrative support: R Qin, B Peng; (III) Provision of study materials or patients: P Zhang, Q Zhu; (IV) Collection and assembly of data: T Li, Y Chen; (V) Data analysis and interpretation: Z Lin, R Hu; (VI) Manuscript writing: All authors; (VII) Final approval of manuscript: All authors.

#These authors contributed equally to this work.

Correspondence to: Xiaowen Zhang. State Key Laboratory of Respiratory Disease, Department of Otolaryngology, Head & Neck Surgery, The First Affiliated Hospital of Guangzhou Medical University, Guangzhou, China. Email: entxiaowen@VIP.163.com; Hua Qin. State Key Laboratory of Respiratory Disease, Department of Otolaryngology, Head & Neck Surgery, The First Affiliated Hospital of Guangzhou Medical University, Guangzhou, China; Interdisciplinary Center of Sleep Medicine, Charité-Universitätsmedizin Berlin, Berlin, Germany. Email: hua.qin12@hotmail.com.

Background: Narcolepsy is a neurological sleep disorder associated with immune response. However, the autoimmune basis for narcolepsy remains unclear. This study aimed to evaluate the causal relationship between immune cells and narcolepsy type 1 (NT1).

Methods: We used a two-sample Mendelian randomization (MR) method to investigate the associations between 731 immune cell traits and NT1 based on a genome-wide association study (GWAS) database from the FinnGen consortium. The inverse-variance weighted (IVW) method was used as the primary method, followed by sensitivity analyses, including the MR-Egger intercept test, Cochran’s Q test, and MR pleiotropy residual sum and outlier (MR-PRESSO). Additional mediation analysis was conducted to investigate the mediating effect of 91 cytokines on immune cells to facilitate the immune processes in NT1.

Results: Immune cell traits showed significant causal associations with NT1. Risk-associated traits mainly involved human leukocyte antigen-DR (HLA-DR)-related monocyte phenotypes, T-cell-related traits, natural killer (NK) cell-related traits, and natural killer T (NKT) cell-related traits, whereas protective traits mainly involved CD4+ T-cell-related, plasmacytoid dendritic cell-related, monocyte-related, and B-cell-related phenotypes. Overall, ten immune cell traits were associated with an increased risk of narcolepsy, whereas five were associated with a reduced risk. Mediation analysis further indicated that interleukin-6 (IL-6) mediates the immune-inflammatory pathway linking monocyte-related traits to NT1. These findings remained consistent in all sensitivity analyses.

Conclusions: Our study identified immunophenotypes that are related to the development of NT1, providing insight into the autoimmune pathogenesis of NT1 and subsequent immunotherapy.

Keywords: Narcolepsy; immunocyte phenotypes; Mendelian randomization (MR); immune cell traits; cytokines


Submitted Nov 20, 2025. Accepted for publication Apr 06, 2026. Published online Jun 29, 2026.

doi: 10.21037/jtd-2025-aw-2380


Highlight box

Key findings

• Multiple immune cell phenotypes, particularly monocyte-, T-cell-, NK/NKT-cell-, and myeloid-related traits, were causally associated with narcolepsy type 1 (NT1).

What is known and what is new?

• NT1 is widely considered to have an autoimmune basis, and our Mendelian randomization study provides genetic evidence that specific immune cell phenotypes and IL-6-related inflammatory pathways contribute to disease susceptibility.

What is the implication, and what should change now?

• These findings strengthen the autoimmune framework of NT1 and support further investigation of targeted immunomodulatory strategies.


Introduction

Narcolepsy is a neurological sleep disorder with a prevalence of 0.02–0.05% of the population worldwide (1). It is thought to be caused by selective loss or dysfunction of hypocretin neurons in the lateral hypothalamus, as hypocretin (also known as orexin) is essential for sleep regulation (2-4). It is characterized by excessive daytime sleepiness, cataplexy, hypnagogic hallucinations, sleep paralysis, and rapid eye movement (REM) sleep intrusion during wakefulness (5). Narcolepsy is divided into two distinct phenotypes based on hypocretin deficiency, narcolepsy type 1 (NT1) and narcolepsy type 2 (NT2). They are classified based on hypocretin-1 (hcrt-1) levels in the cerebrospinal fluid (CSF), presence of cataplexy, and human leukocyte antigen (HLA) gene susceptibility (6). Specifically, NT1 experiences cataplexy and shows a low or undetectable CSF hcrt-1, whereas NT2 presents with less severe clinical symptoms and has a normal CSF hcrt-1 (7). There is also a strong genetic association between narcolepsy type and HLA locus particularly human leukocyte antigen-DQ (HLA-DQ) and human leukocyte antigen-DR (HLA-DR) (8), with 89–98% of NT1 patients positive for HLA-DQB1*06:02 (9). Although these markers can be used to distinguish NT1 from NT2, approximately 10% of NT2 can progress to NT1 over the course of the disease, suggesting that hypocretin cells undergo a progressive destructive process during this transition (4,7,10,11).

HLA regions are closely associated with autoimmune diseases, signaling their potential role in narcolepsy pathophysiology (12,13). However, the lack of studies investigating hypocretin neuron-specific antibodies makes it difficult to define narcolepsy autoimmune disease (14,15). Despite this, data from genetic, immunologic, and clinical studies provide information on the immune-mediated pathophysiology of NT1 (9,16,17). Increasing evidence indicates that T cell is a key contributor to the attack of hypocretin neurons in narcolepsy, as high levels of certain T cell types have been found to target proteins normally expressed in neurons in the brain (18,19). However, the contribution of CD4+ T cells and CD8+ T cells to the autoimmune response against orexin-producing neurons remains unclear (20,21). Other immune cells such as B cells, natural killer T (NKT) cells, and monocytes are presumed to be involved in disease initiation and progression (10).

However, the mechanism through which hypocretin-producing cells activate immune cell destruction is underexplored. Epidemiologic studies show an increased incidence of narcolepsy after upper airway infections and H1N1 pandemic vaccination (22,23). Combined detection of HLA positivity and anti-streptococcal antibodies was identified close to the onset of narcolepsy (24). It is thus assumed that molecular mimicry or bystander activation towards infectious triggers and H1N1 vaccines are crucial immunological pathways for neuronal destruction (17,25,26). Furthermore, genome-wide association analysis demonstrated that genetic factors associated with autoimmunity, such as the T cell receptor alpha (TCRα) locus, tumor necrosis factor superfamily member 4 (TNFSF4/OX40L), Cathepsin H (CTSH) and P2RY11-DNMT1 loci, are present in narcolepsy (16,27,28). There is also evidence of immunoreactivity due to increased serum anti-Tribbles homolog 2 (TRIB2) antibodies in narcolepsy patients, correlated with the onset and severity of cataplexy and sleepiness (14,29,30). Thus, genetic predisposition and environmental triggers may contribute to the initiation of the autoimmune onset of narcolepsy with cataplexy. Increased levels of cytokines such as interleukin (IL)-4, IL-6, and tumor necrosis factor (TNF)-alpha tightly promote the ongoing disease process (31,32), but their mediation effect in the causal pathway towards disease risk is still controversial (33,34). In addition, narcolepsy is closely correlated with an elevated risk of other autoimmune diseases (e.g., multiple sclerosis and systemic lupus erythematosus), suggesting a shared onset (16).

Similar to other autoimmune diseases, narcolepsy has a detrimental impact on health outcomes, leading to a low life quality, and increased cardiovascular risk and all-cause mortality (35-37). Additionally, the efficacy of immunotherapy in the prevention or reduction of the loss of the hypocretin neurons remains controversial. Inconsistent results obtained from immunotherapy trials aimed at improving autoimmune conditions in patients with narcolepsy make it a challenge to confirm treatment benefits (38). This is because specific risk immune system genes and relevant signaling pathways in previous studies have given rise to different potential immunocyte phenotypes that regulate immune reactions and autoimmunity in narcolepsy. Therefore, the establishment of a causal relationship between hypocretin-targeting immune cells and narcolepsy is critical for determining suitable therapeutic targets and improving the efficacy of immunotherapy. However, the identification of autoreactive immune cell subtypes underlying immune-mediated responses against neurons in narcolepsy has not been well documented. Mendelian randomization (MR) is a robust epidemiological method that estimates the genetic variants linked to risk factors to determine the causal relationship between exposure and outcome (39). This is an effective method because it prevents confounding factors and reverse causation bias from influencing results. MR analysis has been employed to identify the associations between sleep traits and various diseases (40). Therefore, we performed two-sample MR and mediation analyses to evaluate the complex genetic pattern of the interactions among the immune system, inflammatory factors, and NT1. We present this article in accordance with the STROBE-MR reporting checklist (available at https://jtd.amegroups.com/article/view/10.21037/jtd-2025-aw-2380/rc).


Methods

Study design

We assessed the potential causal relationship between immune cell traits and narcolepsy using a genome-wide association study (GWAS) dataset obtained from the FinnGen consortium (41). The application of an MR analysis that uses genetic variation as an instrumental variable (IV) to avoid unbiased causal effects is required to adhere to the assumptions of correlation, independence, and exclusion as follows: (I) the genetic variants are strongly correlated with exposure factors; (II) there was no significant association between the IVs and confounding factors; (III) the genetic variants only affect outcomes via exposure pathways. The study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments. The institutional review board of the FinnGen Consortium approved the GWAS protocol, and all participants provided signed informed consent. No additional ethical approval was required for this study, as it is a secondary analysis of the publicly available FinnGen GWAS dataset.

GWAS data sources

Immune cell traits as exposure data

The GWAS summary-level dataset of immune cells is publicly available from the FinnGen GWAS database directory (https://www.finngen.fi/fi). The initial dataset was sourced from the FinnGen Consortium for European ancestry. A total of 731 immune phenotypes (GCST90001391–GCST90002121) were incorporated into this study (42). Immune profiles included absolute cell (AC) counts (n=118), median fluorescence intensities (MFI) reflecting surface antigen levels (n=389), morphological parameters (MP) (n=32), and relative cell (RC) counts (n=192). The AC, MFI, and RC profiles encompassed classical dendritic cells (CDCs), T-cell maturation, bone marrow cells, monocytes, TBNK [T cells, B cells, and natural killer (NK) cells], and Treg (regulatory T cells) panels. The MR profiles were comprised of TBNK panels and CDC (cytotoxic lymphocytes). GWAS data of immune cells were studied using a sample of 313,00 individuals of European descent and approximately 21,303,723 genetic variants of single nucleotide polymorphisms (SNPs) affecting 731 immune cell traits.

Cytokines as exposure data

The GWAS summary-level dataset of cytokines was publicly accessed from GWAS Catalog (https://www.ebi.ac.uk/gwas/downloads). 91 cytokines (GCST90274758 to GCST90274848, N=14,824) were included in the study (43).

NT1 as outcome data

The summary results of the narcolepsy traits were publicly sourced from the directory of the FinnGen GWAS database. This dataset included 313,001 participants (198 patients with narcolepsy and 312,803 controls), comprising 86 males and 112 females. Genetic data for NT1 (GWAS ID, finn-bG6_NARCOCATA, narcolepsy, and cataplexy) encompassed 16,380,329 genetic variants of the SNPs. The mean age at disease onset of all participants with narcolepsy was 39.8 years with a mean 36.6-year-old males and mean 43.9-year-old females.

IV selection

IVs were selected from SNPs in the summary-level datasets of GWAS of immune cells and NT1. We used a significance level of 5×10−8 for each IV of immune cell traits to boost the pool of SNPs available for the MR analyses. Furthermore, we used a chain imbalance coefficient of r2 and region width of 10,000 kb to eliminate the influence of pleiotropy on the results. The F-statistic was applied to diminish the potential impact of weak genetic instrument variables and to exclude IVs with F<10. SNPs with a minor allele frequency <0.01 were excluded to ensure data robustness, and SNPs associated with potential confounding factors were removed to prevent their impact on the data.

Statistical analysis

MR results were expressed as odds ratios (ORs) with the 95% confidence intervals (CIs). The causal relationships between 731 immunophenotypes and NT1 were estimated using two-sample MR. The inverse-variance weighted (IVW) method with Bonferroni correction was conducted as the primary analysis. In addition, we conducted four other MR methods, including the MR-Egger, weighted median, weighted mode, and simple mode, to verify the robustness of our results. For exploratory analysis, immune cells and cytokines with significant causal effects on NT1 were included in the mediation analysis (step 1 and step 2 by a two-step MR approach in Figure 1) (44). Multiple MR analyses were performed to explore the mediating effect of cytokines in the pathway from immune cells to narcolepsy, to determine if immune cells were causally associated with cytokines (step 3 in Figure 1). For sensitivity analysis, we performed Cochran’s Q test and MR Egger intercept test to address heterogeneity and pleiotropy within the summary estimates and selected IVs (45). To further verify the robustness of our results further, we used the MR pleiotropy residual sum and outlier (MR-PRESSO) test and leave-one-out analysis to identify and exclude outliers. All aforementioned MR methods and sensitivity analyses were applied to evaluate the causal associations of 91 cytokines with NT1 and immunophenotypes with cytokines. A p-value of <0.05, with a consistent direction of IVW and MR-Egger, was considered statistically significant. Bonferroni correction was applied to further validate the MR estimates of immunocyte phenotypes (P<0.008 for B cells, maturation stages of T cells, and monocyte panel; P<0.01 for cDC panel; P<0.01 for myeloid cell panel; and P<0.005 for TBNK and Treg panel). All analyses were conducted in R 4.2.1 with TwoSampleMR package and MRPRESSO package.

Figure 1 Flowchart of Mendelian randomization and mediation analysis on causal associations among immune cells, cytokines and narcolepsy type 1.

Results

Immunocyte phenotypes

Risk factors of immune phenotype for NT1

The summarized results of the IVW method used to describe immunophenotypic patterns significantly associated with narcolepsy were shown in Figures 2,3. After performing the Bonferroni correction, we identified 10 immune phenotypes as risk factors for the development of NT1 (Figure 2) including.

Figure 2 Forest plots shows causal associations between immune cell traits as risk factors and narcolepsy type 1. AC, absolute cell count; CI, confidence interval; HLA DR, human leukocyte antigen-DR; IVW, inverse variance weighting; MR, Mendelian randomization; NK, natural killer cell; NKT, natural killer T cell; nSNP, number of single-nucleotide polymorphisms; OR, odds ratio; SSC-A, side scatter area; TBNK, T cells; Treg, regulatory T cell.
Figure 3 Forest plots shows causal associations between immune cell traits as protective factors and narcolepsy type 1. BAFF-R, B-cell activating factor receptor; cDC, conventional dendritic cell; CI, confidence interval; DC, dendritic cell; IVW, inverse variance weighting; MR, Mendelian randomization; nSNP, number of single-nucleotide polymorphisms; OR, odds ratio; SSC-A, side scatter area; TBNK, T cells.

(I) TBNK panel: HLA DR++ monocyte %leukocyte (OR =3.39, 95% CI: 2.15–5.35, P=1.44E-07), CD45 on CD4+ (OR =2.04, 95% CI: 1.3–3.19, P=0.002), SSC-A on NKT (OR =1.83, 95% CI: 1.22–2.75, P=0.003); HLA DR on HLA DR+ NK (OR =1.63, 95% CI: 1.25–2.11, P<0.001); (II) Treg panel: CD25 on CD45RA+ CD4 not Treg (OR =1.78, 95% CI: 1.21–2.62, P=0.003); (III) monocyte panel: HLA DR on monocyte (OR =1.74, 95% CI: 1.39–2.19, P=1.89E-06), HLA DR on CD14+ CD16 monocyte (OR =1.64, 95% CI: 1.2–2.24, P=0.002), HLA DR on CD14+ monocyte (OR =1.64, 95% CI: 1.18–2.28, P=0.003); (IV) maturation stages of T cell panel : naïve CD4+ AC (OR =1.66, 95% CI: 1.15–2.39, P=0.007); (V) myeloid cell panel: HLA DR on CD33br HLA DR+ CD14dim (OR =1.63, 95% CI: 1.3–2.05, P=2.30E-05).

Protector factors of immune phenotype for NT1

We observed 5 immune phenotypes as protective factors for NT1 after implementing Bonferroni correction, as shown in Figure 3 including:

(I) TBNK panel: SSC-A on CD4+ (OR =0.6, 95% CI: 0.44–0.81, P<0.001); (II) cDC panel: CD80 on plasmacytoid DC (OR =0.59, 95% CI: 0.41–0.83, P=0.003), CD80 on CD62L+ plasmacytoid DC (OR =0.59, 95% CI: 0.41–0.83, P=0.003); (III) monocyte panel: HLA DR on CD14+ CD16+ monocyte (OR =0.54, 95% CI: 0.42–0.7, P=1.85E-06); (IV) B cell panel: B-cell activating factor receptor on CD20 (OR =0.38, 95% CI: 0.24–0.58, P=1.39E-05).

Cytokine

The results of the IVW method for the association between cytokines and NT1 without Bonferroni correction were shown in Figure 3.

The 4 following inflammatory factors are correlated with an elevated incidence of NT1:

IL-6 levels (OR =2.53, 95% CI: 1.17–5.45, P=0.02), Eotaxin levels (OR =2.03, 95% CI: 1.19–3.46, P<0.001), Neurturin levels (OR =1.73, 95% CI: 1.01–2.95, P=0.046), TNF-beta levels (OR =1.42, 95% CI: 1.04–1.93, P=0.03).

Of these, 4 cytokines reduce the risk of developing NT1:

TNF ligand superfamily member 14 levels (OR =0.59, 95% CI: 0.38–0.92, P=0.02), eukaryotic translation initiation factor 4E-binding protein 1 levels (OR =0.49, 95% CI: 0.24–0.99, P=0.046), IL-15 receptor subunit alpha levels (OR =0.48, 95% CI: 0.26–0.86, P=0.02), and macrophage colony-stimulating factor 1 levels (OR =0.47, 95% CI: 0.26–0.84, P=0.01).

Mediation analysis of cytokines on immune cells and NT1

Two-step MR analysis was used to determine the impact of immune traits on NT1 via cytokines. (Figure 1, Figure S1, and Table S1), We found that IL-6 serves as a mediator in the immune-inflammatory pathway. Both HLA-DR on CD14+ monocytes and HLA-DR on monocytes exhibited weak correlations with IL-6 levels (OR =1.03, 95% CI: 1.00–1.06, P=0.049 and OR =1.02, 95% CI: 1.00–1.05, P=0.04, respectively). No evidence of pleiotropy or weak IVs were found (Table S2).

Sensitivity analyses

The sensitivity analyses shown in Tables 1,2 generally supported the robustness and reliability of the observed associations between the identified immunocyte phenotypes, cytokines, and NT1. Although Cochran’s Q test indicated heterogeneity among the IVs for certain immune cell traits, neither the MR-Egger intercept test nor the MR-PRESSO global test detected significant directional horizontal pleiotropy or outlier-driven bias (all pleiotropy P>0.05). These results indicate that the observed heterogeneity is attributable to differences in SNP-specific effect estimates rather than clear directional instrumental bias. Accordingly, the absence of directional pleiotropy confirms the validity and robustness of our primary IVW estimates.

Table 1

The results of sensitivity analyses on immune cell traits to the development of narcolepsy type 1

Exposure Heterogeneity Horizontal pleiotropy
Cochran’s Q test MR-egger intercept test MR-PRESSO
Q Q df P Egger intercept SE P Outliers Global test P
Risk factors
   HLA DR++ monocyte %leukocyte 5.568 5 0.35 −0.208 0.259 0.467 NA 0.45
   CD45 on CD4+ 9.498 11 0.58 −0.061 0.09 0.511 NA 0.61
   SSC-A on NKT 19.528 15 0.19 0.094 0.154 0.552 NA 0.19
   CD25 on CD45RA+ CD4 not Treg 16.401 17 0.50 −0.063 0.091 0.5 NA 0.56
   HLA DR on monocyte 17.645 12 0.13 −0.195 0.086 0.046 NA 0.17
   Naive CD4+ AC 23.61 19 0.21 −0.116 0.1 0.261 NA 0.22
   HLA DR on CD14+ CD16 monocyte 30.686 18 0.03 −0.177 0.089 0.062 NA 0.08
   HLA DR on CD14+ monocyte 30.76 17 0.02 −0.182 0.092 0.067 NA 0.07
   HLA DR on CD33br HLA DR+ CD14dim 26.168 17 0.07 −0.228 0.096 0.031 NA 0.054
   HLA DR on HLA DR+ NK 33.707 22 0.053 0.164 0.084 0.064 NA 0.06
Protective factors
   SSC-A on CD4+ 24.21 27 0.62 0.1 0.06 0.107 NA 0.68
   CD80 on plasmacytoid DC 7.448 16 0.96 −0.046 0.101 0.655 NA 0.98
   CD80 on CD62L+ plasmacytoid DC 7.469 16 0.96 −0.047 0.101 0.65 NA 0.98
   HLA DR on CD14+ CD16+ monocyte 17.209 14 0.25 0.12 0.092 0.211 NA 0.20
   BAFF-R on CD20 9.845 11 0.54 0.049 0.122 0.698 NA 0.64

AC, absolute cell count; BAFF-R, B-cell activating factor receptor; DC, dendritic cell; df, degree of freedom; MR-PRESSO, Mendelian randomization pleiotropy residual sum and outlier; NK, natural killer cell; NKT, natural killer T cell; SE, standard error; SSC-A, side scatter area; Treg, regulatory T cell.

Table 2

The results of sensitivity analyses on associations between cytokine and narcolepsy type 1

Exposure Heterogeneity Horizontal pleiotropy
Cochran’s Q test MR-egger intercept test MR-PRESSO
Q Q df P Egger intercept SE P Outliers Global test P
Risk factors
   IL-6 17.703 12 0.13 0.07 0.111 0.54 NA 0.12
   Eotaxin 25.253 26 0.51 −0.03 0.06 0.62 NA 0.54
   Neurturin 26.005 24 0.35 −0.046 0.059 0.45 NA 0.35
   TNF-beta 41.421 41 0.45 0.025 0.049 0.62 NA 0.50
Protective factors
   TNFSF14 29.09 36 0.79 −0.051 0.044 0.25 NA 0.82
   EIF4EBP1 12.712 17 0.76 0.06 0.087 0.50 NA 0.72
   IL-15RA 14.168 19 0.77 0.039 0.076 0.61 rs2836464, NA 0.049
   M-CSF1 12.903 26 0.99 −0.036 0.066 0.59 NA 0.99

df, degree of freedom; EIF4EBP1, eukaryotic translation initiation factor 4E-binding protein 1; IL-15RA, interleukin-15 receptor subunit alpha; IL-6, interleukin-6; M-CSF1, macrophage colony-stimulating factor 1; MR-PRESSO, Mendelian randomization pleiotropy residual sum and outlier; Q, Cochran’s Q statistic; SE, standard error; TNF, tumor necrosis factor; TNFSF14, tumor necrosis factor ligand superfamily member 14.


Discussion

To the best of our knowledge, this is the first MR analysis to investigate the causal relationship between multiple immunophenotypes and NT1 expression. We observed that 10 immunocyte traits in TBNK cell, Treg, monocyte cells, maturation stages of T cell, and myeloid cells were positively linked to the development of narcolepsy. Four inflammatory factors, including IL-6, eotaxin, neurturin, and TNF-beta, increase inflammatory responses to facilitate disease development. Additional mediation analysis showing the immune-inflammatory pathways suggested that IL-6 produced by HLA-DR on CD14+ monocytes and HLA-DR on monocytes is a pivotal mediator of NT1. These results extend the current findings on the basis of neurological autoimmunity in the pathophysiology of narcolepsy.

The genetic association of the HLA allele and immune pathway analysis in narcolepsy assumes an autoimmune origin for the disease (46). This is supported by many attempts to explore multiple immune cell subtypes involved in cell-specific autoimmune attack to hypocretin-producing neurons (47,48). Single-cell analysis of the immune signature of narcolepsy demonstrated that the activation of CD4+ and CD8+ T cells and cytokine production by B cells catalyze the susceptibility to narcolepsy (47). However, no differences were observed between the levels of CD4+ and CD8+ T cells, B cells, NK cells, NKT cells, and monocytes in peripheral blood mononuclear cells (PBMCs) samples (47). In contrast, our results not only confirmed the contribution of T cells, but also showed that NK cells, NKT cells, monocytes, and myeloid cells contribute to the immune reaction, which has not been fully investigated in narcolepsy (Figure 2).

The T cell-driven autoimmune response against hypocretin neurons causes the initiation and late-stage progression of narcolepsy. Importantly, questions regarding the contribution of CD4+ and CD8+ T cells as well as other immune cells or antibodies remain unanswered and require further investigation. CD4+ T cells indirectly target hypocretin-producing neurons because they bind to HLA class II, including HLA-DQ6 and HLA-DR, which are normally absent in neurons (8,49). CD8+ T cells can directly kill hypocretin-producing neurons via a pathway that recognizes antigens bound to HLA class I proteins that are expressed in neurons (49). Latorre et al. observed that high levels of T cells (a high frequency of hypocretin-specific CD8+T cells) in blood samples recognized and targeted hypocretin proteins present in the hypothalamus (50). Pedersen et al. observed that patients with NT1 had increased hypocretin-specific CD8+ T cell levels compared to HLA-DQB1*06:02-positive controls, possibly because of the activation of pre-existing CD8+ T cells facilitated by exogenous antigens (51). Bernard-Valnet et al. showed that cytotoxic CD8+ T cells potentially act as effector cells that cause hypocretin neuron destruction owing to T cell infiltration and cytolytic granule polarization in mouse models (52).

Some subsets of CD4+ T cells have been reported to be associated with the development of autoimmune diseases (e.g., type 1 diabetes, multiple sclerosis, and rheumatoid arthritis) (53-56). However, Latorre et al. showed that a large proportion of autoreactive CD4+ T cells failed to identify hypocretin peptides bound to HLA-DQB1*06:02-restricted HLA-DQ6 proteins but not HLA-DR proteins, suggesting a more nuanced mechanism (50). Despite efforts to identify T cells that target neuronal proteins, there is no strong recognition of narcolepsy-related T cell subtypes that initiate disease onset, as the blood and CSF samples used in previous studies were from a wide range of disease progressions (57). However, we found that three phenotypes of CD4+ T cells (CD45 on CD4+, naïve CD4+ AC, and CD25 on CD45RA+ CD4, not Tregs) were associated with an increased risk of narcolepsy. CD45 is a receptor-type protein tyrosine phosphatase that is essential for T cell activation via TCR signaling. CD45RA is a phenotypic marker of naïve T cells. CD25 expressed on activated T cells is the alpha chain of interleukin-2 receptor (IL-2R). CD25 expression on CD45RA+ CD4+ T cells, but not on Tregs, is likely to serve as an activated effector T cell. These findings support the hypothesis that genetic variants within TCR genes regulate autoimmunity, and that these immune responses are highly predisposed to narcolepsy susceptibility. Subsequent protein production is likely responsible for the observed symptoms. CD4+ and CD8+ T cells induce high levels of P2Y11 protein associated with increased sleep fragmentation as well as increased levels of IL-2, IL-4, IL-13, IL-21, and TNF, which affect B cell differentiation and are involved in age-dependent progression. Additionally, cytokine expression patterns are potentially associated with disease duration, which further elucidates the inflammatory process and autoimmune etiology. For example, CD69+ T cells are associated with short disease duration, whereas IL-4 producing T cell are associated with long disease duration (58).

In addition to a T cell-mediated immune response, other autoimmune mediators also increase the risk of developing narcolepsy via potential immunological pathways responsible for hypocretin neuronal loss. Upregulation of HLA-DR expression in monocytes/macrophages tends to be activated toward infection or inflammation. In this study, we found that elevated expression of HLA-DR++ on leukocytes (HLA-DR++ monocyte % leukocytes) had the highest predisposition to narcolepsy (OR =3.39). As antigen-presenting cells (APCs), they interact with CD4+ T-cell for recognition and activation, by processing and presenting hypocretin-specific antigens. Moreover, the expression of HLA-DR in monocytes and their subsets (CD14+ CD16 and CD14+) was significantly correlated with a high risk of narcolepsy. The four immune cell traits (HLA DR++ monocyte %leukocyte, HLA DR on monocytes, HLA DR on CD14+ CD16 monocytes, and HLA DR on CD14+ monocytes) observed in this study attained causal significance for antigen processing, antigen presentation, and downstream T cell stimulatory effects. These results are consistent with the APC mechanisms of narcolepsy.

Moreover, an immunochip study implicated an association between narcolepsy and polymorphisms in CTSH (27). CTSH is an endosomal lysosomal enzyme highly expressed in monocytes, dendritic cells, and B cells. Activated CTSH can modify the presentation of major histocompatibility complex class II peptide repertoire to T cells, inducing disease development. Moreover, monocytes in patients with NT1 secrete increased levels of IL-6, which is related to sleepiness, compared to those in controls that are positive for HLA-DQB1*06:02. Mediation analysis showed that IL-6 is associated with HLA DR on CD14+ monocytes and HLA DR on monocytes, which further indicated the process of antigen presentation to T cells by HLA class II-restricted APCs in the pathogenesis of narcolepsy. Although there are inconsistent results of etiological effect of circulating cytokines on narcolepsy (18,33,59), our results support the causal link between monocytes and IL-6 in immune-inflammatory processes. From a therapeutic perspective, our findings provide insights beyond orexin receptor agonist-based symptomatic therapies (46,60,61). Whereas orexin receptor agonists primarily compensate for downstream hypocretin deficiency, our results highlight upstream immune mechanisms, particularly monocyte-related antigen presentation and IL-6-mediated inflammation, that causally contribute to disease susceptibility. These observations indicate that the interplay between monocytes and IL-6 signaling is highly relevant to the immune pathophysiology of NT1 and highlights promising targets for future mechanistic and immunomodulatory studies (46,61). Our findings also help explain why broad immunosuppressive strategies such as systemic steroid administration cannot be universally recommended in NT1. In our analysis, not all cytokines were associated with increased NT1 risk, demonstrating that non-selective immunosuppression can inadvertently suppress both detrimental and potentially protective immune pathways. In addition, NT1 is often recognized after substantial hypocretin neuronal loss has already occurred, which may further limit the potential benefit of late-stage, non-specific immunosuppressive treatment (46,62).

The regulatory functional mechanisms of HLA-DR on CD33br HLA-DR+ CD14dim from a subset of myeloid cells, dendritic cells, or other myeloid-derived suppressor cells in narcolepsy have not yet been reported. However, we found that the presence of NK cells, such as HLA-DR on HLA-DR+ NK cells, plays a role in the immunogenetic aspects of narcolepsy. HLA-DR expression in cytotoxic lymphocytes suggests their activation or immunoregulatory function. HLA-DR+ NK cells share phenotypic features of both NK cells and dendritic cells, which produce cytokines such as interferon-gamma (IFN-γ) and present certain antigens to CD4+ and CD8+ T cells, leading to activation and proliferation. CD4+ T cells produce IFN-γ, which stimulates CD8+ T cells and drives a T helper 1 (Th1) cell-mediated cellular response. Corroborating this mechanism, high plasma levels of IFN-γ have been observed in patients with narcolepsy. Vuorela et al. demonstrated the upregulation of IFN-γ stimulated by T-cell reactivity to single peptides from hemagglutinin, neuraminidase, or nucleoprotein, either in vivo or NT1 patients (26). Similarly, SSC-A in NKT cells is also associated with IFN-γ production. Finally, this study found that IL-6, eotaxin, neurturin, and TNF-α have a tendency to increase the risk of narcolepsy. In this context, it still remains debated that changes in levels of serum, plasma and CSF cytokine can be indicative of an ongoing disease process or are consequences of hypocretin neuronal loss in narcolepsy (18,34,63). Further studies on the role of cytokines as biomarkers for hypocretin-deficient narcolepsy are warranted. At the same time, the cytokines identified in this study are pleiotropic rather than exclusive to NT1. Inflammatory mediators such as IL-6 and TNF-related factors have also been implicated in other autoimmune diseases associated with narcolepsy, including systemic lupus erythematosus and multiple sclerosis (62,64). These findings reflect shared autoimmune-inflammatory pathways across immune-mediated disorders rather than an NT1-specific cytokine signature. In NT1, their biological relevance is better interpreted in conjunction with the accompanying immune-cell context rather than as disease-specific biomarkers in isolation. Nevertheless, in NT1 these cytokines operate within a specific context of monocyte activation and HLA-DR-related antigen presentation, as supported by our mediation findings.


Conclusions

Our study demonstrated a causal relationship between ten risk immunophenotypes and NT1 via bidirectional MR analysis, highlighting the involvement of various immune cells in immune modulation and autoimmunity in narcolepsy. IL-6 expression mediated by monocyte subsets plays a role in the pathogenesis of narcolepsy. These findings reveal the underlying mechanism of autoimmune and inflammatory dysregulation in narcolepsy, providing evidence for immunological intervention and the prevention of narcolepsy.

Figure 4 Mendelian randomization results of the relationship between cytokines and narcolepsy type 1. Beta, effect estimate; CI, confidence interval; EIF4EBP1, eukaryotic translation initiation factor 4E-binding protein 1 levels; IL-15RA, interleukin-15 receptor subunit alpha levels; IVW, inverse variance weighting; M-CSF1, macrophage colony-stimulating factor 1 levels; MR, Mendelian randomization; nSNP, number of single-nucleotide polymorphisms; OR, odds ratio; TNF, tumor necrosis factor; TNFSF14, tumor necrosis factor ligand superfamily member 14 levels.

Acknowledgments

We thank the participants and investigators of the FinnGen study and GWAS Resource for the public data.


Footnote

Reporting Checklist: The authors have completed the STROBE-MR reporting checklist. Available at https://jtd.amegroups.com/article/view/10.21037/jtd-2025-aw-2380/rc

Peer Review File: Available at https://jtd.amegroups.com/article/view/10.21037/jtd-2025-aw-2380/prf

Funding: This work was supported by the National Youth Foundation of China (No. 82301290, to H.Q., No. 82300040, to R.Q.) and State Key Laboratory of Respiratory Disease (No. SKLRD-Z-202509, to H.Q.).

Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at https://jtd.amegroups.com/article/view/10.21037/jtd-2025-aw-2380/coif). The authors have no conflicts of interest to declare.

Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. The study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments. The institutional review board of the FinnGen Consortium approved the GWAS protocol, and all participants provided signed informed consent. No additional ethical approval was required for this study, as it is a secondary analysis of the publicly available FinnGen GWAS dataset.

Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.


References

  1. Longstreth WT Jr, Koepsell TD, Ton TG, et al. The epidemiology of narcolepsy. Sleep 2007;30:13-26. [Crossref] [PubMed]
  2. Nishino S, Ripley B, Overeem S, et al. Hypocretin (orexin) deficiency in human narcolepsy. Lancet 2000;355:39-40. [Crossref] [PubMed]
  3. Thannickal TC, Moore RY, Nienhuis R, et al. Reduced number of hypocretin neurons in human narcolepsy. Neuron 2000;27:469-74. [Crossref] [PubMed]
  4. Lopez R, Barateau L, Evangelista E, et al. Temporal Changes in the Cerebrospinal Fluid Level of Hypocretin-1 and Histamine in Narcolepsy. Sleep 2017;40:zsw010. [Crossref] [PubMed]
  5. Scammell TE. Narcolepsy. N Engl J Med 2015;373:2654-62. [Crossref] [PubMed]
  6. AASM. International classification of sleep disorders. 3rd ed. Darien (IL): American Academy of Sleep Medicine; 2014.
  7. Bassetti CLA, Adamantidis A, Burdakov D, et al. Narcolepsy - clinical spectrum, aetiopathophysiology, diagnosis and treatment. Nat Rev Neurol 2019;15:519-39. [Crossref] [PubMed]
  8. Mignot E, Lin L, Rogers W, et al. Complex HLA-DR and -DQ interactions confer risk of narcolepsy-cataplexy in three ethnic groups. Am J Hum Genet 2001;68:686-99. [Crossref] [PubMed]
  9. Mignot E, Hayduk R, Black J, et al. HLA DQB1*0602 is associated with cataplexy in 509 narcoleptic patients. Sleep 1997;20:1012-20.
  10. Andlauer O, Moore H 4th, Hong SC, et al. Predictors of hypocretin (orexin) deficiency in narcolepsy without cataplexy. Sleep 2012;35:1247-55F. [Crossref] [PubMed]
  11. Baumann CR, Mignot E, Lammers GJ, et al. Challenges in diagnosing narcolepsy without cataplexy: a consensus statement. Sleep 2014;37:1035-42. [Crossref] [PubMed]
  12. Dendrou CA, Petersen J, Rossjohn J, et al. HLA variation and disease. Nat Rev Immunol 2018;18:325-39. [Crossref] [PubMed]
  13. Gough SC, Simmonds MJ. The HLA Region and Autoimmune Disease: Associations and Mechanisms of Action. Curr Genomics 2007;8:453-65. [Crossref] [PubMed]
  14. Smith AJ, Jackson MW, Neufing P, et al. A functional autoantibody in narcolepsy. Lancet 2004;364:2122-4. [Crossref] [PubMed]
  15. Kornum BR. Narcolepsy type 1: what have we learned from immunology? Sleep 2020;43:zsaa055. [Crossref] [PubMed]
  16. Ollila HM, Sharon E, Lin L, et al. Narcolepsy risk loci outline role of T cell autoimmunity and infectious triggers in narcolepsy. Nat Commun 2023;14:2709. [Crossref] [PubMed]
  17. Tesoriero C, Codita A, Zhang MD, et al. H1N1 influenza virus induces narcolepsy-like sleep disruption and targets sleep-wake regulatory neurons in mice. Proc Natl Acad Sci U S A 2016;113:E368-77. [Crossref] [PubMed]
  18. Liblau RS, Latorre D, Kornum BR, et al. The immunopathogenesis of narcolepsy type 1. Nat Rev Immunol 2024;24:33-48. [Crossref] [PubMed]
  19. Hallmayer J, Faraco J, Lin L, et al. Narcolepsy is strongly associated with the T-cell receptor alpha locus. Nat Genet 2009;41:708-11. [Crossref] [PubMed]
  20. Tafti M, Lammers GJ, Dauvilliers Y, et al. Narcolepsy-Associated HLA Class I Alleles Implicate Cell-Mediated Cytotoxicity. Sleep 2016;39:581-7. [Crossref] [PubMed]
  21. Jiang W, Birtley JR, Hung SC, et al. In vivo clonal expansion and phenotypes of hypocretin-specific CD4(+) T cells in narcolepsy patients and controls. Nat Commun 2019;10:5247. [Crossref] [PubMed]
  22. Sarkanen TO, Alakuijala APE, Dauvilliers YA, et al. Incidence of narcolepsy after H1N1 influenza and vaccinations: Systematic review and meta-analysis. Sleep Med Rev 2018;38:177-86. [Crossref] [PubMed]
  23. Han F, Lin L, Warby SC, et al. Narcolepsy onset is seasonal and increased following the 2009 H1N1 pandemic in China. Ann Neurol 2011;70:410-7. [Crossref] [PubMed]
  24. Aran A, Lin L, Nevsimalova S, et al. Elevated anti-streptococcal antibodies in patients with recent narcolepsy onset. Sleep 2009;32:979-83. [Crossref] [PubMed]
  25. Luo G, Ambati A, Lin L, et al. Autoimmunity to hypocretin and molecular mimicry to flu in type 1 narcolepsy. Proc Natl Acad Sci U S A 2018;115:E12323-32. [Crossref] [PubMed]
  26. Vuorela A, Freitag TL, Leskinen K, et al. Enhanced influenza A H1N1 T cell epitope recognition and cross-reactivity to protein-O-mannosyltransferase 1 in Pandemrix-associated narcolepsy type 1. Nat Commun 2021;12:2283. [Crossref] [PubMed]
  27. Faraco J, Lin L, Kornum BR, et al. ImmunoChip study implicates antigen presentation to T cells in narcolepsy. PLoS Genet 2013;9:e1003270. [Crossref] [PubMed]
  28. Kornum BR, Kawashima M, Faraco J, et al. Common variants in P2RY11 are associated with narcolepsy. Nat Genet 2011;43:66-71. [Crossref] [PubMed]
  29. Cvetkovic-Lopes V, Bayer L, Dorsaz S, et al. Elevated Tribbles homolog 2-specific antibody levels in narcolepsy patients. J Clin Invest 2010;120:713-9. [Crossref] [PubMed]
  30. Kawashima M, Lin L, Tanaka S, et al. Anti-Tribbles homolog 2 (TRIB2) autoantibodies in narcolepsy are associated with recent onset of cataplexy. Sleep 2010;33:869-74. [Crossref] [PubMed]
  31. Dauvilliers Y, Jaussent I, Lecendreux M, et al. Cerebrospinal fluid and serum cytokine profiles in narcolepsy with cataplexy: a case-control study. Brain Behav Immun 2014;37:260-6. [Crossref] [PubMed]
  32. Tanaka S, Honda Y, Honda M. Identification of differentially expressed genes in blood cells of narcolepsy patients. Sleep 2007;30:974-9. [Crossref] [PubMed]
  33. Mohammadi S, Mayeli M, Saghazadeh A, et al. Cytokines in narcolepsy: A systematic review and meta-analysis. Cytokine 2020;131:155103. [Crossref] [PubMed]
  34. Kornum BR, Pizza F, Knudsen S, et al. Cerebrospinal fluid cytokine levels in type 1 narcolepsy patients very close to onset. Brain Behav Immun 2015;49:54-8. [Crossref] [PubMed]
  35. Tadrous R, O’Rourke D, Mockler D, et al. Health-related quality of life in narcolepsy: A systematic review and meta-analysis. J Sleep Res 2021;30:e13383. [Crossref] [PubMed]
  36. Jennum PJ, Plazzi G, Silvani A, et al. Cardiovascular disorders in narcolepsy: Review of associations and determinants. Sleep Med Rev 2021;58:101440. [Crossref] [PubMed]
  37. Ohayon MM, Black J, Lai C, et al. Increased mortality in narcolepsy. Sleep 2014;37:439-44. [Crossref] [PubMed]
  38. Szabo ST, Thorpy MJ, Mayer G, et al. Neurobiological and immunogenetic aspects of narcolepsy: Implications for pharmacotherapy. Sleep Med Rev 2019;43:23-36. [Crossref] [PubMed]
  39. Emdin CA, Khera AV, Kathiresan S. Mendelian Randomization. JAMA 2017;318:1925-6. [Crossref] [PubMed]
  40. Sonti S, Grant SFA. Leveraging genetic discoveries for sleep to determine causal relationships with common complex traits. Sleep 2022;45:zsac180. [Crossref] [PubMed]
  41. Kurki MI, Karjalainen J, Palta P, et al. FinnGen provides genetic insights from a well-phenotyped isolated population. Nature 2023;613:508-18. [Crossref] [PubMed]
  42. Orrù V, Steri M, Sidore C, et al. Complex genetic signatures in immune cells underlie autoimmunity and inform therapy. Nat Genet 2020;52:1036-45. [Crossref] [PubMed]
  43. Zhao JH, Stacey D, Eriksson N, et al. Genetics of circulating inflammatory proteins identifies drivers of immune-mediated disease risk and therapeutic targets. Nat Immunol 2023;24:1540-51. [Crossref] [PubMed]
  44. Carter AR, Sanderson E, Hammerton G, et al. Mendelian randomisation for mediation analysis: current methods and challenges for implementation. Eur J Epidemiol 2021;36:465-78. [Crossref] [PubMed]
  45. Verbanck M, Chen CY, Neale B, et al. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet 2018;50:693-8. [Crossref] [PubMed]
  46. Valizadeh P, Momtazmanesh S, Plazzi G, et al. Connecting the dots: An updated review of the role of autoimmunity in narcolepsy and emerging immunotherapeutic approaches. Sleep Med 2024;113:378-96. [Crossref] [PubMed]
  47. Hartmann FJ, Bernard-Valnet R, Quériault C, et al. High-dimensional single-cell analysis reveals the immune signature of narcolepsy. J Exp Med 2016;213:2621-33. [Crossref] [PubMed]
  48. Khajavi L, Nguyen XH, Queriault C, et al. The transcriptomics profiling of blood CD4 and CD8 T-cells in narcolepsy type I. Front Immunol 2023;14:1249405. [Crossref] [PubMed]
  49. Mahoney CE, Cogswell A, Koralnik IJ, et al. The neurobiological basis of narcolepsy. Nat Rev Neurosci 2019;20:83-93. [Crossref] [PubMed]
  50. Latorre D, Kallweit U, Armentani E, et al. T cells in patients with narcolepsy target self-antigens of hypocretin neurons. Nature 2018;562:63-8. [Crossref] [PubMed]
  51. Pedersen NW, Holm A, Kristensen NP, et al. CD8(+) T cells from patients with narcolepsy and healthy controls recognize hypocretin neuron-specific antigens. Nat Commun 2019;10:837. [Crossref] [PubMed]
  52. Bernard-Valnet R, Yshii L, Quériault C, et al. CD8 T cell-mediated killing of orexinergic neurons induces a narcolepsy-like phenotype in mice. Proc Natl Acad Sci U S A 2016;113:10956-61. [Crossref] [PubMed]
  53. Sun L, Su Y, Jiao A, et al. T cells in health and disease. Signal Transduct Target Ther 2023;8:235. [Crossref] [PubMed]
  54. Dittel BN. CD4 T cells: Balancing the coming and going of autoimmune-mediated inflammation in the CNS. Brain Behav Immun 2008;22:421-30. [Crossref] [PubMed]
  55. Orrù V, Steri M, Sole G, et al. Genetic variants regulating immune cell levels in health and disease. Cell 2013;155:242-56. [Crossref] [PubMed]
  56. Todd JA, Wicker LS. Genetic protection from the inflammatory disease type 1 diabetes in humans and animal models. Immunity 2001;15:387-95. [Crossref] [PubMed]
  57. Ramberger M, Högl B, Stefani A, et al. CD4+ T-Cell Reactivity to Orexin/Hypocretin in Patients With Narcolepsy Type 1. Sleep 2017;40:zsw070. [Crossref] [PubMed]
  58. Lecendreux M, Churlaud G, Pitoiset F, et al. Narcolepsy Type 1 Is Associated with a Systemic Increase and Activation of Regulatory T Cells and with a Systemic Activation of Global T Cells. PLoS One 2017;12:e0169836. [Crossref] [PubMed]
  59. Tanaka S, Honda M, Toyoda H, et al. Increased plasma IL-6, IL-8, TNF-alpha, and G-CSF in Japanese narcolepsy. Hum Immunol 2014;75:940-4. [Crossref] [PubMed]
  60. Dauvilliers Y, Plazzi G, Mignot E, et al. Oveporexton, an Oral Orexin Receptor 2-Selective Agonist, in Narcolepsy Type 1. N Engl J Med 2025;392:1905-16. [Crossref] [PubMed]
  61. Vringer M, Zhou J, Gool JK, et al. Recent insights into the pathophysiology of narcolepsy type 1. Sleep Med Rev 2024;78:101993. [Crossref] [PubMed]
  62. Su X, Yu H, Lei Q, et al. Systemic lupus erythematosus: pathogenesis and targeted therapy. Mol Biomed 2024;5:54. [Crossref] [PubMed]
  63. Hinze-Selch D, Wetter TC, Zhang Y, et al. In vivo and in vitro immune variables in patients with narcolepsy and HLA-DR2 matched controls. Neurology 1998;50:1149-52. [Crossref] [PubMed]
  64. Mazziotti V, Crescenzo F, Turano E, et al. The contribution of tumor necrosis factor to multiple sclerosis: a possible role in progression independent of relapse? J Neuroinflammation 2024;21:209. [Crossref] [PubMed]
Cite this article as: Li T, Chen Y, Zhang P, Zhu Q, Lin Z, Peng B, Qin R, Steenbergen N, Hu R, Qin H, Zhang X. Immunocyte phenotypes underlying the causal autoimmune features in narcolepsy type 1. J Thorac Dis 2026;18(7):714. doi: 10.21037/jtd-2025-aw-2380

Download Citation