Abstract
Background Immune response changes have been reported in amyotrophic lateral sclerosis (ALS), but their clinical relevance remains undetermined. Therefore, we aim to evaluate the relationships between blood leukocyte subpopulations and prognosis of ALS.
Methods A longitudinal cohort of 288 ALS patients with up to 5 years of follow-up during 2015-2020 were recruited at the only tertiary referral center for ALS in Stockholm, Sweden. Routine differential leukocyte counts, and determination of lymphocyte subpopulations including an extended T cell panel with flow cytometry, collected at diagnosis and at regular intervals thereafter. The primary outcome was risk of death (alternatively use of invasive ventilation) after diagnosis of ALS. The secondary outcomes included repeatedly measured functional status - through Amyotrophic Lateral Sclerosis Functional Rating Scale-revised (ALSFRS-R) score and disease progression rate. Cox model was used to evaluate the associations between leukocytes and risk of death. Generalized estimating equation model (GEE) was used to assess the correlation between leukocytes and ALSFRS-R score and disease progression rate.
Results The counts of leukocytes, neutrophils and monocytes increased gradually over time since diagnosis and were negatively correlated with ALSFRS-R score, but not associated with risk of death or disease progression rate. Focusing on lymphocyte subpopulations, increasing counts of natural killer (NK) cells (HR=0.61, 95% CI= [0.42-0.88] per SD increase) and proportions of Th2-diffrentiated CD4+ central memory (CM) T cells (HR=0.64, 95% CI= [0.48-0.85] per SD increase) were correlated with a lower risk of death. Increasing proportions of CD4+ effector memory cells re-expressing CD45RA (EMRA) T cells (HR=1.39, 95% CI= [1.01-1.92] per SD increase) and CD8+ T cells (HR=1.38, 95% CI= [1.03-1.86] per SD increase) were associated with a higher risk of death. None of the lymphocyte subpopulations was correlated with ALSFRS-R score or disease progression rate.
Conclusion Our findings suggest a dual role of immune responses in ALS prognosis, where neutrophils and monocytes primarily reflect functional status whereas NK cells and different T lymphocyte populations act as prognostic markers for survival. The findings also provide insights for cell-specific treatment for ALS.
Funding This work was supported by the European Research Council (ERC) Starting Grant (MegaALS, No. 802091), the Swedish Research Council (No. 2019-01088), Karolinska Institutet (Strategic Research Area in Epidemiology and Senior Researcher Award), and China Scholarship Council.
Introduction
Amyotrophic lateral sclerosis (ALS) is a rare but devastating neurodegenerative disease. Although a genetic cause has been demonstrated for some cases of ALS, the etiology remains unknown for most of the patients with ALS. There is currently no cure or effective treatment available for ALS. A range of potential disease mechanisms have however been proposed, with potential for informing on novel therapeutic targets1.
Neuroinflammatory features, including local glial activation and T cell infiltration in the central nervous system (CNS), are well documented in ALS2. Animal studies have demonstrated that altering the function of microglia and infiltrating T cells to the CNS affects disease progression in experimental ALS3-5. Human studies have also documented signs of systemic immune activation in ALS patients, compared with healthy controls, suggesting that peripheral immune activation may play a role in human ALS through the peripheral-central neuroimmune crosstalk6,7. Most of the previous studies7-9, however, are not population-based and have limited sample size, raising the concern of potentially insufficient internal (e.g., selection bias and chance finding) and external (e.g., lack of generalizability) validity. Further, few studies have recorded immune cells longitudinally after ALS diagnosis, as most previous studies relied on a single cross-sectional measurement10,11.
The main purpose of this study was to determine cellular immune changes occurring over time since diagnosis of ALS in a large community-based sample. For this purpose, we enrolled a longitudinal cohort of ALS patients, representing a large proportion of the entire ALS population in Stockholm, Sweden, to 1) describe the temporal changes of peripheral leukocytes over time since ALS diagnosis, 2) assess the association of immune cell dynamics with the risk of death after ALS diagnosis, and 3) evaluate the correlations of different cell populations with the functional status and disease progression rate of ALS.
Methods
Study cohort
The Swedish Motor Neuron Disease (MND) Quality Registry was established in 2015, collecting information on clinical characteristics, biological test results, and quality of life outcomes from >80% of MND patients in Sweden12. Since 2017, the MND Quality Registry has included 99% of MND patients in the Stockholm area among whom 97.1% are diagnosed with ALS. All ALS diagnoses were made by a specialist in neurology and followed up by a neuromuscular specialist, and met the diagnostic requirement of definite, probable, probable laboratory-supported, or possible ALS according to the revised El Escorial criteria13,14. To ensure the accuracy of diagnosis, all patients in the registry are re-evaluated annually to update diagnosis, whenever needed.
Through the MND Quality Registry, we first identified 420 patients with ALS diagnosed from the start of the registry until October 7th, 2020, in Stockholm. We reviewed the medical records of these patients to identify information on peripheral leukocyte populations (i.e., differential leukocyte counts). During this process, we excluded 12 patients who were not diagnosed at the Karolinska University Hospital – the only tertiary referral center for ALS in Stockholm, three patients with unknown time of symptom onset, 82 patients lacking leukocyte counts, and 35 patients with counts outside of the stipulated observation period (i.e., from three months before date of diagnosis until October 7th, 2020). The final analysis cohort included 288 patients (68.6%), with at least one recorded differential leukocyte counts during the observation period. The included patients did not differ significantly in terms of demographic and clinical characteristics from the excluded patients (Table 1). About half of the patients had a single measurement of leukocytes whereas the other half had been sampled two or more times.
Characteristics of the 288 patients with amyotrophic lateral sclerosis (ALS) included in the study, compared with the entire population of ALS patients during the study period in Stockholm, Sweden
Among the ALS patients diagnosed at the Karolinska University Hospital and with a date of symptom onset (N=405), we further performed flow cytometry in 92 patients (“FlowC cohort”) to determine lymphocyte subpopulations (i.e., T, B and natural killer [NK] cells) as well as an extended T lymphocyte panel “FITMaN” - an internationally standardized panel reported by the Flow Immunophenotyping Technical Meeting at NIH15. Compared to the main study cohort, patients of the FlowC cohort were slightly younger and more likely to have a limb onset (eTable 1). In both the main and FlowC cohorts, we followed the ALS patients from date of diagnosis or first cell measurement (differential leukocyte counts or FlowC), whichever came later, until occurrence of the outcome of interest (i.e., death or use of invasive ventilation) or October 7th, 2020, whichever came first.
Blood samples and flow cytometric analysis
All sample processing and analyzing procedures were according to the validated protocol at the Departments of Clinical Chemistry (differential leukocyte counts) and Clinical Immunology and Transfusion Medicine (FlowC), Karolinska University Hospital. All analyses were performed during daytime, within 24 hours of sampling. Differential leukocyte counts were done on a Sysmex XN-9000 (Sysmex, Kobe, Japan). FlowC was implemented in clinical routine based on the standardized phenotyping panel by the Human Immunophenotyping Consortium with a set of defined 8-color antibody cocktails15. The experiments were performed on a triple-laser Beckman Coulter Gallios and analyzed by Kaluza Software (Beckman Coulter, Brea, CA).
Differential leukocyte counts included neutrophils, lymphocytes, monocytes, eosinophils, and basophils. We did not include eosinophils and basophils in the analysis as they in most cases were low to undetectable. The FITMaN panel included measures of 23 lymphocyte subpopulations, including 1) counts of B cells, NK cells, and T cells; 2) %s of CD4+ and CD8+ T cell subtypes (i.e., naive, central memory [CM], effector memory [EM], and effector memory cells re-expressing CD45RA [EMRA] T cells based on CCR7 and CD45RA expression, as well as Th1, Th2, and Th17 of CM and EM CD4+ T cells based on CXCR3 and CCR6 expression); and 3) subtypes of activated CD4+ and CD8+ T cells based on the expression of HLA-DR and CD38 (i.e., CD4+HLA-DR+CD38- cells, CD4+HLA-DR+CD38+ cells, CD8+HLA-DR+CD38- cells, and CD8+HLA-DR+CD38+ cells). T cells were gated from a lymphocyte (FCS/SSC) gate as cells expressing CD45 and CD3. The unit of cell count was 10^9/L whereas the %s were expressed as the proportions of the lower-level immune cell populations out of the upper-level immune cells. The reference values for FlowC were reported in the form of 5th to 95th percentiles using reference normal ranges obtained from 50 healthy adults. All cell counts and %s were retrieved from patients’ clinical files.
Outcomes of interest
The primary study outcome was risk of death or use of invasive ventilation after ALS diagnosis, identified from the MND Quality Registry. The secondary study outcomes included functional status measured through the Amyotrophic Lateral Sclerosis Functional Rating Scale – revised (ALSFRS-R) and disease progression rate. ALSFRS-R is a questionnaire-based scale that measures the motor function and disease severity of ALS patients and is considered the gold standard measure of disability progression16. Higher ALSFRS-R score indicates better functional status. We acquired information on all available ALSFRS-R scores for the ALS patients from the MND Quality Registry. Progression rate measures the rate of ALSFRS-R decline and was calculated by dividing the difference between 48 (the full score) and measured ALSFRS-R score at a specific time point by the time difference between time of symptom onset to the measurement time of ALSFRS-R (in months). Progression rate is an independent prognostic predictor for ALS17.
Statistical analysis
We first calculated the mean levels of measured cell populations among ALS patients, by summarizing all measurements from three months before diagnosis until end of follow-up. To visualize the temporal patterns of the cell populations, we drew a trajectory line of all measurements for each cell type and each patient. We then used the locally estimated scatterplot smoothing (LOESS) curves with 95% confidence intervals (CIs) to show the temporal pattern of the predicted median level of each cell type after ALS diagnosis. If the CIs did not overlap with the normal ranges of the cell populations, we considered the observed levels among ALS patients to be statistically deviant from normal ranges. We also used linear mixed model to assess the within-individual temporal changes of cell populations after ALS diagnosis. In this analysis, we included a random intercept to account for the initial differences between individuals and adjusted for age at diagnosis and sex. We first analyzed all ALS patients together, then separated by site of onset and presence of C9orf72 expansions. Patients with C9orf72 expansions have been suggested to demonstrate a different immune phenotype compared with ALS patients without such expansions18,19.
We next used Cox model to derive hazard ratio (HR) and 95% CI to assess the association of different cell populations with the risk of death, after adjustment for other prognostic indicators of ALS including age at diagnosis, sex, site of onset, diagnostic delay, ALSFRS-R score, time difference between the measure of ALSFRS-R score and diagnosis, body mass index (BMI), and time difference between the measure of BMI and diagnosis. We obtained information on BMI from the MND Quality Registry. We used time since diagnosis as the underlying time scale and cluster robust variance estimation to account for dependence of repeated measurements. In this analysis, we log-transformed the numbers of leukocytes, neutrophils, lymphocytes and monocytes, as well as %s of CD4+ EM cells, CD4+ EMRA cells, Th2 of CD4+ EM cells, Th17 of CD4+ EM cells, Th2 of CD4+ CM cells, CD8+ T cells, CD8+ CM cells, CD4+HLA-DR+CD38- cells, CD4+HLA-DR+CD38+ cells, CD8+HLA-DR+CD38- cells, and CD8+HLA-DR+CD38+ cells, to achieve a better normal distribution. Values of other cell populations were used as is (i.e., without transformation). For all markers, we estimated the effect size per standard deviation (SD) increase. Because some patients had their first cell measurements after diagnosis, we also took into account left truncation in all analyses.
The proportional hazards assumption of the Cox model was assessed using the Schoenfeld residual test. After stratifying the analysis by site of onset which deviated from the assumption, the assumption became satisfied for all other variables. As HRs obtained in the stratified analysis were nearly the same as those obtained without stratification, for simplicity and consistency, we reported all findings from the original models (without stratification by site of onset).
Among the 288 patients, 57 (19.8%) were diagnosed before the start of the MND Quality Registry. In a sensitivity analysis, we excluded these patients with the aim to see if the results obtained in the main analysis would pertain to a cohort of incident patients. This analysis was not performed for the FlowC cohort as all the 92 patients were incident patients. We also conducted a sensitivity analysis by only focusing on the first cell measure of each patient to examine whether the main results would differ after removal of repeated measurements. We then performed another sensitivity analysis by excluding patients with C9orf72 expansions.
We further used a generalized estimating equation (GEE) model to assess the correlations of cell populations with functional status (i.e., ALSFRS-R score) and disease progression rate measured at the same time as the cell markers. The GEE model considers the correlations between repeated measurements within the same individual. Like the survival analysis above, we used log-transformed values for some cell measures whereas original values for others and estimated the effect size per SD increase of the cell measures, after adjustment for age at diagnosis and sex.
Finally, to assess whether the results on cell measures and the study outcomes would be influenced by ongoing infections, we performed another sensitivity analysis after excluding all measurements taken with the presence of an ongoing infection. This was done by reviewing clinical records of all patients with a high count of leukocyte (>8.8×109/L) for diagnosis of any infectious disease, tests for infection, and patient-reported infectious symptoms.
All analyses were performed using R 3.6.2. A two-sided p value of <0.05 was considered statistically significant. To correct for multiple testing, we also computed the Benjamini-Hochberg false discovery rate (FDR)20.
Results
Leukocytes and lymphocyte subtypes in ALS
eTable 2 shows the distribution of leukocyte populations (N=288 patients) and lymphocyte subpopulations (N=92 patients) across all measures after ALS diagnosis. The vast majority of the cell populations were within the normal range, except for CD8+ CM cells, CD4+HLA-DR+CD38+ cells and CD8+HLA-DR+CD38+ cells, which were above the normal range.
In the main cohort, the levels of leukocytes, neutrophils and monocytes increased progressively over time, especially from 20 months after diagnosis onward (Figure 1). These trends were statistically significant, with or without adjustment for age and sex, and remained statistically significant after correction for multiple testing (Table 2). In contrast, no clear temporal trend was noted for lymphocytes. These results remained largely similar when stratifying the patients by site of onset or presence of C9orf72 expansions (eFigure 1).
The black lines show measured levels of leukocyte populations for each patient. The blue lines and shadow areas show the mean levels of leukocyte populations with 95% confidence intervals. Pink areas indicate normal range. P for trend shows the P value of within-individual temporal change of each cell population after taking into account the relatedness of repeated measurements.
Temporal changes of leukocyte populations after diagnosis of amyotrophic lateral sclerosis (ALS), a cohort study of 288 patients with ALS in Stockholm, Sweden*
In the FlowC cohort, no clear temporal trend was noted for any lymphocyte subpopulation, although ALS patients demonstrated persistently higher proportions of CD8+ CM, CD4+HLA-DR+CD38+ and CD8+HLA-DR+CD38+ cells than the reference ranges (eFigure 2). After adjustment for age and sex, there was a decreasing %s of naïve CD4+ T cells whereas increasing %s of CD4+ EMRA, CD4+HLA-DR+CD38- and CD8+HLA-DR+CD38- cells since ALS diagnosis (eTable 3). Patients with limb onset had lower levels of CD8+ CM and CD4+HLA-DR+CD38+ cells compared with patients with other anatomical disease onset, whereas carriers of C9orf72 expansions had higher levels of NK cells and T cells, but lower levels of naïve CD4+ T cells and CD8+ EM cells, than other patients (eFigure 3).
Survival
During a median follow-up of 1.1 years, we observed 163 deaths or use of invasive ventilation among the 288 patients of the main cohort. No association was noted between the level of leukocytes, neutrophils, lymphocytes, or monocytes with risk of death (eTable 4). This result did not change after excluding patients diagnosed before the start of the MND Quality Registry, focusing on first cell measure only, or excluding patients with C9orf72 expansions (eTable 5).
In the FlowC cohort, we found higher NK cell counts and %s of Th2-diffrentiated CD4+ CM cells to be associated with lower risk of death, whereas higher %s of CD4+ EMRA cells and CD8+ T cells were associated with higher risk of death (Figure 2). These results were largely similar after restricting the analysis to first measure of lymphocytes (data not shown) or after excluding patients with C9orf72 expansions (eFigure 4).
Functional status and disease progression rate
In the main cohort, a higher level of leukocytes, neutrophils, or monocytes was associated with a lower ALSFRS-R score measured at the time of sampling, whereas no such correlation was evident for lymphocyte counts (Table 3). There was, however, no association of leukocytes, neutrophils, lymphocytes, or monocytes with disease progression rate.
Cross-sectional correlations between leukocyte populations and ALS functional rating scale-revised (ALSFRS-R) score and disease progression rate, a cohort study of 288 ALS patients in Stockholm, Sweden*
In the FlowC cohort, none of the lymphocyte subtypes was associated with ALSFRS-R score or disease progression rate measured at the time of sampling (eTable 6).
Sensitivity analysis for ongoing infection
We excluded 15 patients in the main cohort who had been sampled with the presence of infection. The results on risk of death, ALSFRS-R score and disease progression rate remained similar although some results lost statistical significance (eTable 7 and eTable 8). There was no patient with sampling during ongoing infection in the FlowC cohort.
Discussion
We here report a longitudinal cohort study of temporal dynamics of white blood cell populations in 288 ALS patients in Stockholm, Sweden. We found higher counts of blood leukocytes, neutrophils and monocytes to be associated with a lower functional status, but not with the risk of death after ALS diagnosis or disease progression rate. In a subsample of 92 patients, we found that higher NK cell counts and proportions of Th2-diffrentiated CD4+ T cells were associated with a lower risk of death, whereas higher counts of CD8+ T cells and proportions of CD4+ EMRA T cells were associated with higher mortality risk.
ALS patients demonstrated gradually increasing counts of leukocytes, neutrophils and monocytes over time since diagnosis. Benjamin et al. also reported increasing counts of leukocytes and neutrophils in ALS patients21. Accumulating monocytes have also been shown in the cervical and lumbar spinal cord of ALS model over time22. The finding that peripheral leukocytes, neutrophils, and monocytes were correlated to ALSFRS-R, but not progression rate or survival, indicates that these leukocyte subtypes may serve better as markers for functional status than prognosis. A previous study also found CD16 expression on neutrophils and non-classical monocytes to correlate with ALS disease severity23. As neutrophils and monocytes have phagocytic function, increased levels of circulating neutrophils and monocytes might indicate enhanced muscle damage, explaining the correlation with deteriorating functional status. It has also been reported that circulating monocytes from ALS patients preferentially differentiate to a M1 proinflammatory phenotype and produce more Interleukin 6 (IL-6) and tumor necrosis factor α (TNF-α), compared with monocytes from healthy controls24. Regardless, although it has been proposed that peripheral neutrophils and monocytes are recruited to the CNS through a disrupted brain-blood barrier (BBB) and affect the disease progression of experimental ALS through secreting proinflammatory cytokines and influencing other cells25-27, the functional relevance in human ALS is still unclear22,28.
We found a higher level of blood-borne NK cells to be associated with a lower risk of death after ALS diagnosis. NK cells are considered a critical part of the innate immunity and function through lysing infected, oncogenic, apoptotic, and MHC class I-deficient cells29. Although NK cells are known to penetrate BBB and interact with microglia, astroglia, and neurons30, the role of NK cells in disease progression of ALS has rarely been addressed. NK cells may exert a cytotoxic role in the brain by its natural function, as suggested by a study showing that depletion of NK cells prolonged survival in ALS mouse models31. As infiltration of NK cells to CNS may lead to a decrease in peripheral NK populations, the protective role of higher peripheral NK cells against risk of death, as observed in the present study, might be partly attributable to a lower level of its infiltration to the CNS. On the other hand, NK cells may indeed also play a neuroprotective role. For instance, findings from experimental autoimmune encephalomyelitis, a model of multiple sclerosis, showed NK cells to suppress neuroinflammation, diminish tissue damage and protect motoneurons32,33. Patients with multiple sclerosis have also been shown to experience clinical improvement through NK cell expansions34. NK cells have also been shown to exert a protective role in brain by removing viral infection and activated microglia30.
We also found that higher Th2 differentiation of CD4+ T cells was associated with a lower risk of death after ALS diagnosis, corroborating previous findings of the neuroprotective role of CD4+ T cells in ALS5,35. Animal studies have suggested Th1 cells and Th17 cells to promote neuroinflammation by producing proinflammatory cytokines and enhancing microglia-mediated neurotoxic effects36 whereas Th2 cells and regulatory T cells (Tregs) to suppress neuroinflammation by producing anti-inflammatory cytokines and enhancing microglia-mediated neuroprotective effects37,38. Although not statistically significant, our study indeed found a trend for Th1-differentiatated CD4+ CM cells and Th17-differentiated CD4+ CM cells to be positively associated with the risk of death after ALS diagnosis.
A novel finding of the present study is that higher proportions of CD8+ T cells and CD4+ EMRA T cells were associated with a higher risk of death after ALS diagnosis. The precise underlying mechanisms linking together these cell types with ALS survival are unknown. Previously, CD8+ T cells have been found to be present in the spinal cord only at the end stage of ALS5. However, such cells have been suggested to have a pathogenic potential in multiple sclerosis and to communicate with mononuclear phagocytes39,40. Previous studies also indicated that CD4+ EMRA cells could demonstrate cytotoxic features and express the chemokine receptor CX3CR1 in the setting of Dengue virus infection41, which are associated with cytotoxic lymphocytes with cytoplasmic granules containing perforin and granzymes42. Further studies are clearly needed to study more in detail the functional relevance of CD8+ and CD4+ EMRA T cells in ALS.
Expansions in C9orf72, the most common genetic cause of ALS, have been shown to be associated with immune features43, including activation of microglia and elevated levels of peripheral inflammatory cytokines44. In our study, although patients with or without C9orf72 expansions did not differ greatly in terms of major leukocyte populations, patients with C9orf72 expansions appeared to display changes in certain lymphocyte subpopulations including NK, T, naïve CD4+ T and CD8+ EM cells compared with patients without such expansions. Because of the relatively small number of patients with C9orf72 expansions in the study, these results should however be interpreted with caution until validated further.
Our study is the first, to our knowledge, to report the role of a comprehensive list of immune cells, including neutrophils, lymphocytes, monocytes and detailed T cell phenotypes, on the prognosis of ALS. The strengths of the present study are the large number of ALS patients which were representative of all ALS patients in the source population, the rich information on disease characteristics including genetic causes, the long (up to 5 years after diagnosis) and complete follow-up (through the MND Quality Registry), the availability of both routine cell counts and detailed lymphocyte phenotyping, as well as the access to repeated measures over time. The study also has limitations. First, the cohort was heterogeneous in terms of the numbers of cell measurements and the time intervals between measurements. Some concern of indication bias, for instance, due to impact of infections, might exist. We therefore excluded cell measures taken at the time of infection in a sensitivity analysis. The largely similar results obtained in this analysis relieved this concern to some extent. Further, the FITMaN panel does not include B cell subtypes or Treg cells, which await to be studied further. Finally, causal inferences on the functional implications of the reported associations are not possible due to the observational study design.
In conclusion, our findings suggest a dual role of immune responses in ALS prognosis, where neutrophils and monocytes primarily reflect functional status whereas different T lymphocyte populations act as prognostic markers for survival, which also provide insights for cell-based therapy in prolonging survival in ALS.
Data Availability
The original data are held by the Swedish Motor Neuron Disease Quality Registry and are not publicly available due to Swedish laws. However, any researcher who is interested can access the data by obtaining an ethical approval from regional ethical review board. Detailed information of the Swedish Ethical Review Authority can be found at https://etikprovningsmyndigheten.se. Anonymized data that support the findings of this study are available on reasonable request from the corresponding authors, in agreement with European regulations and Regional Ethical Review Board in Sweden. One source file containing R code for table 2 and 3 has been provided.
Competing interests
The authors have declared that no conflict of interest exists.
Online Supplementary Content
Characteristics of the 92 patients with amyotrophic lateral sclerosis (ALS) included in the analysis of FlowC test, compared with the entire population of ALS patients during the study period in Stockholm, Sweden
Mean levels of leukocyte subpopulations (N=288 patients) and lymphocyte subpopulations (N=92 patients) across all measures
Temporal changes of lymphocyte populations after diagnosis of amyotrophic lateral sclerosis (ALS), analysis of 92 patients with FlowC test*
Associations of leukocyte populations with the risk of death after a diagnosis of amyotrophic lateral sclerosis (ALS), a cohort study of 288 patients with ALS in Stockholm, Sweden
Sensitivity analyses of the associations of leukocyte populations with risk of death after a diagnosis of amyotrophic lateral sclerosis (ALS)*
Cross-sectional correlations between lymphocyte populations and ALS functional rating scale-revised (ALSFRS-R) score and disease progression rate, a cohort study of 92 ALS patients in Stockholm, Sweden*
Sensitivity analyses of the associations of leukocyte populations with ALS functional rating scale-revised (ALSFRS-R) score and disease progression rate, a cohort study of 273 ALS patients in Stockholm, Sweden*
Sensitivity analyses of associations of leukocyte populations with the risk of death after a diagnosis of amyotrophic lateral sclerosis (ALS), a cohort study of 273 ALS patients in Stockholm, Sweden*
Acknowledgements
This work was supported by the European Research Council (ERC) Starting Grant (MegaALS, No. 802091), the Swedish Research Council (No. 2019-01088), Karolinska Institutet (Strategic Research Area in Epidemiology and Senior Researcher Award), and China Scholarship Council. The funders had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.