Incidence, Prevalence, and Geographical Clustering of Motor Neuron Disease in the Netherlands
Citation Manager Formats
Make Comment
See Comments

Abstract
Objective To assess time trends in motor neuron disease (MND) incidence, prevalence, and mortality and to investigate geographic clustering of MND cases in the Netherlands from 1998 to 2017, we analyzed data from the Netherlands Personal Records database, the Netherlands MND Center, and the Netherlands Patient Association of Neuromuscular Diseases.
Methods In this prospective cohort study, Poisson regression was used to assess time trends in MND risk. We calculated age- and sex-standardized, observed, and expected cases for 1,694 areas. Bayesian smoothed risk mapping was used to investigate geographic MND risk.
Results We identified 7,992 MND cases, reflecting an incidence of 2.64 (95% confidence interval [CI] 2.62–2.67) per 100,000 person-years and a prevalence of 9.5 (95% CI 9.1–10.0) per 100,000 persons. Highest age-standardized prevalence and mortality rates occurred at a later age in men than in women (p < 0.001). Unadjusted mortality rates increased by 53.2% from 2.57 per 100,000 person-years in 1998 to 3.86 per 100,000 person-years in 2017. After adjustment for age and sex, an increase in MND mortality rate of 14.1% (95% CI 5.7%–23.2%, p < 0.001) remained. MND relative risk ranged from 0.78 to 1.43 between geographic areas; multiple urban and rural high-risk areas were identified.
Conclusions We found a significant national increase in MND mortality from 1998 through 2017, explained only partly by an aging Dutch population, and a geographic variability in MND risk, suggesting a role for environmental or demographic risk factors.
Glossary
- ALS=
- amyotrophic lateral sclerosis;
- CI=
- confidence interval;
- DPR=
- Netherlands Personal Records;
- ENCALS=
- European Network for the Cure of ALS;
- ICD-10=
- International Classification of Diseases, 10th revision;
- MND=
- motor neuron disease;
- RR=
- relative risk
Motor neuron disease (MND) is a fatal neurodegenerative disease caused by both genetic and environmental factors, many of which remain to be elucidated.1,2 Worldwide MND incidence and prevalence have been shown to be increasing in the last decades.3 Aging of the world population is seen as its primary driver but does not account entirely for the increasing incidence and prevalence.3,4 This indicates that important genetic or potentially preventable environmental risk factors are currently driving the increase in MND risk. Urgent identification of these risk factors is paramount given the debilitating and care-intensive nature of MND and its future effect on health care services.2,5
Although assessing whether MND risk is increasing may aid planning of health services, it does not directly identify underlying risk factors. Geographic studies, on the other hand, can stimulate etiologic research by identifying specific disease clusters. Examples are identification of patient clusters in Guam6 and the Kii Peninsula7 that were subsequently linked to neurotoxic plant consumption and the C9orf72 repeat expansion. Studies of time trends or clusters have, however, often been limited by methodologic issues such as incomplete case ascertainment or geographic risk assessment in a limited area.8
We therefore use multiple, independent, national data sources in the Netherlands, including a nationwide population-based study operational since 2006, to determine the change in MND incidence, prevalence, mortality rates and survival time in the Netherlands over the last 2 decades and to investigate geographic clustering of MND cases.
Methods
Data Sources
Three independent sources of prospectively collected data were used to estimate incidence, prevalence, and mortality of adult-onset MND (that is, amyotrophic lateral sclerosis [ALS], primary lateral sclerosis, and progressive muscular atrophy) in the Netherlands. The first source was the Netherlands Personal Records (DPR) database covering the period from 1998 to 2017. All Dutch residents are registered in the DPR database, and Dutch law requires that all deaths in the Netherlands be registered by a physician according to the ICD-10. MND (ICD-10 codes G12.21–G12.24) is registered as the cause of death even if it was an underlying cause (e.g., in the case of death due to pneumonia). We also selected Alzheimer cases (ICD-10 code G30) to compare the effect of age on Alzheimer and MND mortality rates. As a second source, we used data from the population-based Netherlands MND Center registry, described in more detail elsewhere.9 In brief, patients diagnosed with MND, according to the revised El Escorial criteria,10 have been prospectively registered centrally by the Netherlands MND Center since 2006. Patients were identified via annual screening of large hospital registries and specialized MND rehabilitation clinic registries and by contacting neurologists individually.9 Survival times after date of diagnosis (date of death or date last known to be alive) were obtained for patients in the Netherlands MND registry by checking the online municipal register at 3-month intervals. Third, patients were identified within the register of a national patient advocacy organization, the Netherlands Patients Association for Neuromuscular Diseases.11
Statistical Analysis
Comprehensive MND health surveillance estimates were calculated as follows. Age- and sex-adjusted incidence and mortality rates were calculated by dividing the number of observed cases by the person-years of observation. Cases were identified in the DPR database from 1998 through 2017 and in the Netherlands MND registry from 2006 through 2017. MND prevalence was the number of MND patients alive on December 31, 2017, expressed as a proportion of the total Dutch population. The number of unobserved cases in the MND registry was estimated via 2 methods. First, 2-source capture-recapture methodology was applied to each 5-year age and sex stratum, using data from the MND registry and the Netherlands Patients Association for Neuromuscular Diseases. Capture-recapture methodology aims to correct for underascertainment in health surveillance studies12 and has been successfully applied in other neuromuscular epidemiologic studies.13 The Chapman12 formulas were used to estimate the total number of patients with prevalent MND and the proportion identified by the MND registry. As a sensitivity analysis, we also estimated incidence and prevalence on the basis of the MND deaths per year recorded in the DPR database. Individual survival times were sampled from a Weibull distribution14 and subtracted from the date of death to calculate date of diagnosis. Simulations were repeated 50,000 times to obtain empirical 95% confidence intervals (CIs) for incidence and prevalence estimates.
Next, we assessed whether MND mortality rates changed between 1998 and 2017. We used a Poisson generalized linear model and data from the DPR database to estimate MND mortality rate per year. We incorporated age (quadratic), sex, their interaction, and years since 1998 as covariates with the natural logarithm of population size as an offset. To test both whether the effect of age differs across sexes and whether rates have increased since 1998, we tested each term using a likelihood ratio test. As a sensitivity analysis, we evaluated how loss of competing risks and longevity of a subpopulation susceptible to MND influence MND mortality rates as described elsewhere.15 We evaluated whether patient characteristics of the Netherlands MND registry changed between early periods (2006–2009) and late periods (2014–2017). Cox regression was used to assess whether survival recorded by the Netherlands MND registry changed from 2006 through 2017. The hazard ratio for time in years since 2006 was adjusted for the patient's individual risk profile according to the European Network for the Cure of ALS (ENCALS) survival prediction model.16 The ENCALS risk profile is a relative measure and indicates whether a patient's survival time is shorter or longer relative to other patients. Five prognostic subgroups were defined on the basis of the quintiles of the individual risk profiles.17 Missing values in baseline characteristics, necessary to determine the ENCALS risk profile, were accounted for by creating 100 imputed datasets as described earlier.17 Estimates were pooled across imputations with the Rubin rules.18
Finally, we evaluated the geographic distribution of MND cases in the Netherlands. To map MND risk and to identify both high- and low-risk areas, data from the DPR database were used. As per January 1, 2018, the Netherlands comprises 13,305 geographic areas (i.e., neighborhoods). The DPR database contains detailed individual residential history data, dating from approximately 1990. The address where a case had resided the longest since 1990 was used to calculate the number of observed cases per geographic area. Nationwide mortality rates per 5-year age and sex group were applied to the mean population per time period to calculate the number of expected MND cases per geographic area. To obtain sufficient sample size to conduct spatial smoothing of MND risk, areas with an expected incidence of <3 cases (i.e., ≤6,000 residents from 1998 through 2017) were merged iteratively with adjacent areas with the lowest expected incidence, until all areas had at least an expected incidence of 3. Next, crude standardized mortality ratios were calculated by dividing the observed by the expected number of cases. Standardized mortality ratios were subsequently smoothed to estimate MND relative risk (RR) per area by using a bayesian autoregression model with local random effects as described elsewhere.18 As an example, a RR of 1.3 means that the observed MND rate in that area is estimated to be 30% higher than expected (i.e., the national rate of MND).8 For each area, we determined the posterior probability that the local risk of MND would be either lower or higher than the national risk according to 10,000 draws from the posterior distribution (i.e., the posterior probability that RR <1.0 or RR >1.0, respectively). Bayesian analyses control for multiplicity by using a conservative prior that shrinks all estimates toward the null hypothesis (i.e., an RR of 1.0). Therefore, it is not necessary to adjust the posterior probability for multiplicity. In addition, the prior ensures that the MND RR can be estimated even in rural areas when data are sparse. As an exploratory analysis, we assessed the weighted association between RR and the average age, percentage of men, and population density per area. Population density was analyzed as number of residents per square kilometer. The following definitions from Statistics Netherlands were used: <500 (very low), 500 to 1,000 (low), 1,000 to 1,500 (average), 1,500 to 2,500 (high), and >2,500 (very high).
Standard Protocol Approvals, Registrations, and Patient Consents
This study is reported according to the Standards of Reporting of Neurological Disorders guideline.19 The medical ethics committee and institutional review board of the University Medical Center Utrecht approved this study.
Data Availability
All protocols, analyses, and anonymized data will be shared on request from any qualified investigator. Access to data from the DPR database can be requested from Statistics Netherlands (CBS) via cbs.nl/en-gb/our-services/customised-services-microdata/microdata-conducting-your-own-research/microdata-catalogue.
Results
MND Incidence, Prevalence, and Mortality
From 1998 through 2017, there were, on average, 16.4 million residents and a total of 3.2 million deaths recorded in the Netherlands, of which 7,992 (0.25%) were MND related. The total number of patients with incident MND in this time period was estimated to be 8,676 (95% CI 8,585–8,767), resulting in an average annual incidence of 2.64 (95% CI 2.62–2.67) per 100,000 person-years. Since the start of the MND registry in 2006, the number of incident patients was 5,764, 4,152 (72.0%) of whom were identified by the MND registry. Table 1 presents their patient characteristics; patients identified via the MND registry were, on average, younger (67.7 years vs 69.7 years) and more likely to be male (58.4% vs 54.8%) compared to patients in the DPR database.
Baseline Characteristics of the Netherlands MND Registry and Netherlands Personal Records Database
On December 31, 2017, we identified 1,215 prevalent patients in the MND registry (source 1) and the Netherlands Patient Association for Neuromuscular Diseases (source 2). A total of 953 patients were unique to source 1, 83 were unique to source 2, and 179 patients were identified by both sources. The total number of prevalent patients with MND was estimated at 1,654, resulting in a prevalence of 9.6 (95% CI 9.5–9.8) per 100,000 persons (supplementary table 1, doi.org/10.5061/dryad.4xgxd257d). The simulation-based sensitivity analysis based on data from the DPR database resulted in a similar prevalence estimate (n = 1,639 patients, 9.5 per 100,000 persons).
Overall MND incidence and mortality rates were highest in the 70- to 74-year and 75- to 79-year age groups, respectively (figure 1). Incidence and prevalence data are available in Dryad (figure 1, doi.org/10.5061/dryad.4xgxd257d). There was a differential effect of age for men and women (i.e., different shape of the age-standardized rates graphs in figure 1) on incidence (p < 0.01 for interaction term), prevalence, and mortality (both p < 0.001). For example, highest MND prevalence and mortality rates occurred at later ages for men compared to women (80–84 years vs 75–79 years, respectively). In contrast, highest Alzheimer mortality rates occurred in the oldest age group (>95 years) for both men and women (figure 1D).
This panel plot shows motor neuron disease (MND) (A) incidence, (B) prevalence, and (C) mortality and (D) Alzheimer mortality in the Netherlands from 1998 through 2017 determined via Poisson regression. Dashed lines indicate the age group with maximum rates. There was a differential effect of age for men and women on incidence (p < 0.01 for interaction term), prevalence, and mortality (both p < 0.001). For example, highest MND prevalence and mortality rates occurred at later ages for men than for women. In contrast, highest Alzheimer mortality rates occurred in the oldest age group (>95 years) for both men and women. This suggests that MND risk is not merely a result of aging. Age- and sex-stratified MND incidence was determined on the basis of data from the Netherlands MND registry. MND prevalence was determined by capture-recapture methodology using data from the Netherlands MND registry and Netherlands Patient Association for Neuromuscular Diseases. MND mortality and Alzheimer mortality were determined from mortality records in the Netherlands Personal Records database.
Time Trends
In 2017, 0.31% of all recorded deaths in the Netherlands were MND related, resulting in a lifetime risk of 1 in 323. The absolute number of MND deaths increased from 299 in 1998 to 466 in 2017 with an average annual increase of 2.7% (95% CI 2.4%–3.1%) (figure 2A). Unadjusted estimated mortality rates increased by 53.2% (95% CI 42.0%–65.4%, p < 0.001) from 2.57 in 1998 to 3.86 per 100,000 person-years in 2017 (figure 2B).
(A) Number of motor neuron disease (MND) deaths observed in the Netherlands between 1998 and 2017. Absolute number of MND deaths increased from 299 in 1998 to 466 in 2017 with an average annual increase of 2.6%. (B) Adjusted for population size, mortality rates increased by 53.2% or by 14.1% After additional age and sex adjustment (both p < 0.001) in the time period of 1998 to 2017, dashed lines indicate regression lines. (C) MND mortality rates are provided separately for men (blue) and women (red). Annual increase in MND mortality rates was not different for men and women (p = 0.62 for interaction term). (D) Similar to panel C, we provide the annual male:female ratio, revealing a stable sex ratio of 1.19 (dashed line) over time (p = 0.09 for time trend).
Overall, we found no evidence that the increase in MND mortality rates differed among age groups (p = 0.57 for interaction term). The observed increase did, however, seem to be driven largely by younger age groups; mortality rates increased by 60.3% (95% CI 40.8%–82.6%) in people <70 years of age and by 15.6% (95% CI 4.1%–28.4%) in age groups >70 years. After adjustment for age and sex, MND mortality rates increased by 14.1% (95% CI 5.7%–23.2%, p < 0.001) in this time period; this increase was similar for men and women (p = 0.62 for interaction term) (figure 2C). Consequently, the average male:female ratio among MND cases did not increase significantly over time (p = 0.09) and was, on average, 1.19 (95% CI 1.14–1.25) (figure 2D). Similarly, the mean age at death due to MND in the DPR database was, on average, 69.7 years and did not significantly change during the study period (change in age was 0.62 years, 95% CI −0.23 to 1.47, p = 0.15).
In contrast to the observed increase in MND mortality rates, age- and sex-adjusted mortality rates from competing (i.e., non-MND) causes of death decreased by 1.98% (95% CI 1.96%–2.00%) per year. This annual decrease was larger for men than for women (p < 0.001 for interaction term): 2.48% (95% CI 2.41%–2.54%) and 1.51% (95% CI 1.48%–1.54%) for men and women, respectively. As a sensitivity analysis, we used the Gompertz model15 to evaluate whether longevity of a subpopulation susceptible to MND, as a result of loss of competing risks, may be driving the increase in MND mortality rates. We estimated that the susceptible subpopulation decreased by −253 men per 100,000 (95% CI −529 to 23, p = 0.05) and by −210 women per 100,000 (95% CI −356 to −56, p = 0.006).
Figure 3A shows the annual distribution of the 5 prognostic subgroups defined by the ENCALS personalized prediction model of patients enrolled in the Netherlands MND registry since 2006. According to the ENCALS prediction model, the average risk profile at time of diagnosis has worsened since 2006 (p < 0.001), meaning that in recent years more patients have been diagnosed with a relatively poor risk profile. This finding is in line with table 1, which shows that, even though diagnostic delay had not significantly changed since 2006 (−0.76 months, 95% CI −4.4 to 2.9, p = 0.31), the average age of patients with MND had increased and more patients had been diagnosed with concomitant frontotemporal dementia since 2006. This suggests that the Netherlands MND registry has become more population based by improving recruitment strategies to ascertain patients with incident disease instead of a mixture of patients with incident and prevalent disease (i.e., long-surviving) in earlier register years. Figure 3B shows the adjusted annual probability of survival in the Netherlands. The adjusted risk of death during follow-up decreased by 16.2% during the time period of 2006 to 2017 (hazard ratio per 5 years 0.93, 95% CI 0.88–0.98, p = 0.006). This resulted in an adjusted median survival time that increased from 19.6 months (95% CI 18.7–20.9) to 22.4 months (95% CI 21.1–23.7) for patients diagnosed in 2006 and 2017, respectively.
(A) Annual distribution of the 5 prognostic subgroups defined by the European Network for the Cure of ALS (ENCALS) personalized prediction model. Ideally, each prognostic subgroup should have a prevalence of ≈20%/y. As can be seen, in 2006, the Netherlands motor neuron disease (MND) registry enrolled relatively more very long- and long-surviving patients, indicating that the registry recruited mainly prevalent cases. Fortunately, as recruitment strategies improved to register also short- and very short-surviving subgroups, the Netherlands MND registry has become more population based. (B) Effect of year of diagnosis on survival time was modeled with a Cox proportional hazards model adjusted for the ENCALS risk prediction. Risk of death during follow-up decreased by 16.2% between 2006 and 2017 (hazard ratio per 5 years 0.93, 95% confidence interval [CI] 0.88–0.98, p = 0.006). Adjusted median survival increased from 19.6 months (95% CI 18.7–20.9) in 2006 to 22.4 months (95% CI 21.1–23.7) in 2017.
Mapping of MND Risk
Finally, the probability that local MND risk is higher than expected (i.e., the national MND risk) per geographic area is shown in figure 4. The DPR database contained 1 patient with a missing address (0.01%). After merging of areas to obtain sufficient sample sizes, there were 1,694 areas with a median population size of 9,507 residents and 4.0 observed MND cases per area. We found geographic variation in MND risk with RRs ranging from 0.78 to 1.43, with the highest risk near Ten Boer in the province of Groningen and the lowest risk near Wijk bij Duurstede in the province of Utrecht. The table in figure 4 presents areas that have a probability of at least 0.90 that local MND risk is either higher or lower than national MND risk (i.e., RR >1.0 or RR <1.0). We identified 10 areas, both urban and rural, with a >0.90 probability of either high or low MND risk scattered across the Netherlands. As an exploratory analysis, we assessed whether population density, mean age, and proportion of men per area were related to MND risk. Median population density was 2,591 residents per 1 km2 (range, 82–14,452) and was not significantly related to MND risk: RR changed by −0.015% (95% CI −0.006 to 0.003) for every doubling in population density (p = 0.49). Similarly, mean age and proportion of men were also unrelated. The changes in RR for a 10-year increase in mean age and an absolute increase of 10% in the proportion of men were 0.003 (95% CI −0.024 to 0.0295, p = 0.85) and −0.79 (95% CI -5.94 to 4.35, p = 0.76), respectively.
Motor neuron disease (MND) relative risk for 1,694 areas was calculated and subsequently smoothed based on data from the Dutch Personal Records database from 1998 to 2017. This map shows the probability (Prob) that MND relative risk (RR) per area is greater than expected according to the national rate of MND (i.e., the probability of RR >1.0 based on 10,000 draws from the posterior distribution). Population (Pop) density was estimated as number of residents per 1 km2; the following definitions from Statistics Netherlands were used: <500 (very low), 500 to 1,000 (low), 1,000 to 1,500 (average), 1,500 to 2,500 (high), and >2,500 (very high). Wijk bij D. = Wijk bij Duurstede.
Discussion
In this study, we used multiple independent data sources to provide comprehensive health surveillance estimates of MND in the Netherlands from 1998 to 2017, resulting in an average estimated incidence of 2.64 per 100,000 person-years, a prevalence of 9.5 per 100,000 persons, and a lifetime MND risk of 1 in 323. Mean diagnostic delay did not significantly change while median survival time improved by 3 months from 2006 through 2017, possibly as a result of improved multidisciplinary care. This observed increase in median survival time could imply that the increase in MND incidence may be larger than the observed increase in MND mortality rates since 1998. After age and sex adjustment, MND mortality rates increased by 14.1% over the last 2 decades in both men and women. The 53.2% increase in unadjusted MND mortality rates could have significant consequences for planning of future health care services in that it may indicate a doubling of MND prevalence by the year 2050. Our results suggest that either genetic or potentially preventable environmental risk factors are driving MND risk, and urgent identification is needed, which may be aided by detailed geographic studies of MND risk.
We found a different effect of age when comparing MND and Alzheimer risk. The decline of MND risk in the oldest age groups contrasts with the risk of Alzheimer and suggests that there is a time period of maximal susceptibility and that MND is not only a result of aging. Short survival in cases of MND compared to cases of Alzheimer may also play a role. MND may be exclusive to a small susceptible subpopulation, the majority of whom are deceased by 70 to 80 years of age as a result of either MND or other unrelated causes. There is evidence that loss of competing risks and longevity of this susceptible subpopulation could also be driving the increase in MND mortality rates.15 Indeed, we found that the risk of death due to causes other than MND decreased during the study period. Our findings of a reduction of the size of a susceptible subpopulation and larger increase of MND mortality rates among younger age groups, however, contrast with this hypothesis. Nevertheless, we cannot rule out that longevity of a susceptible subpopulation might had affected the increase of MND mortality rates.
Aging of the general population is an important driver of the increasing incidence of neurodegenerative diseases.20 An earlier study found no statistically significant change in MND incidence in Minnesota in the United States in the period of 1925 to 1998.21 In contrast, several studies have indicated that MND incidence has increased in the last 2 decades.3,4,22 For example, in an epidemiologic study in Northern Italy, MND incidence increased by 14% from 1995 to 2014, mostly in women.4 A possible explanation for the larger increase in MND incidence in more recent years may be that those born during the baby boom after World War II reached ages with highest incidence in the period of 2005 to 2017. Most studies were based on registry data only while using capture-recapture methodology to correct for underascertainment, which can underestimate or overestimate the total number of cases when data sources are not independent.23 In our study, we have therefore performed both capture-recapture and simulation-based sensitivity analyses with 3 independent data sources, including unselected data from a national compulsory health database. In absolute terms, the number of cases per year increased from 298 in 1998 to 466 in 2017. If this increase continues, the lifetime risk of MND would increase to 1 in 205 by 2050. This observed increase is not driven solely by aging of the general population, which could explain only ≈75% of the increase. Increased awareness and timelier referral of patients suspected of having MND since the start of the Netherlands MND Center in 2003 may play a role, although several environmental risk factors that have been linked to MND could be involved such as exposure to organic dust, agricultural pesticides, air pollution, and other occupational hazards.24,–,27
To help develop hypotheses as to which risk factors are currently driving MND risk, we mapped geographic MND risk from 1998 through 2017 to identify high-risk areas. The use of national data is particularly advantageous because assessing geographic risk in a limited area is an approach known to exaggerate the likelihood of identifying high-risk areas.8 There was considerable spatial variation in MND risk, which indicates that associating spatial risk with characteristics of that particular area may be worthwhile (i.e., identify factors that explain spatial variation in risk). For example, by quantifying risk factors such as lifestyle, physical activity levels, and environmental exposures in well-defined geographic areas,27,–,29 we may be able to efficiently identify risk factors. We have illustrated its potential use by associating MND risk with population density, which may act as a proxy for other risk factors.30 Although we found no direct correlation of population density with MND risk, its simplicity could stimulate etiologic research in MND. For example, air pollution could be more severe in cities; by combining spatial risk estimates with exposure levels, on both an individual and a geographic level, one may increase power to find potential associations. A similar approach may hold true for exposure to agricultural pesticides, which are more common in rural areas. An illustrative example of the latter is the increased risk of Parkinson disease in the Netherlands associated with living in the vicinity of agricultural fields where pesticides are used.31 Furthermore, genetic risk factors may also cause geographic variability in MND risk.7 Previously, studies found that rare genetic variants appear to be important drivers of MND risk32 and that genetic variation can be localized geographically, meaning that there is relatively little migration in the Netherlands.33 As a result, local MND incidence may increase via a genetic founder effect when new genetic variations occur. We hypothesize that both environmental risk factors, unique to each area, and local genetic admixture may thus explain geographic variation in MND risk. Quantification of a broad spectrum of risk factors, both on an individual level as per geographic area, remains paramount.
Due to privacy regulations, our study was limited in that we could not confirm a clinical diagnosis in the DPR database for individual cases. Therefore, we cannot rule out some misclassification in the DPR database. Nevertheless, we considered the DPR database to be of a high standard for several reasons. Because registration is mandatory, the DPR database includes all Dutch residents. Furthermore, it is likely that patients with MND will visit a neurologist at least once because the Dutch public healthcare system ensures that there are no financial or physical hurdles to receiving health care. As a result, it is probable that if a patient receives an MND diagnosis, this will also be communicated to that individual’s general practitioner or nursing home physician. Dutch law requires that a physician (usually the general practitioner or nursing home physician) list all possible contributing causes of death on the death certificate. In addition, we found that both age at death and male:female ratio were stable, suggesting that case ascertainment has not changed during the study period. Moreover, the incidence and prevalence estimates that we found are in line with estimates found in other North European countries with compulsory national health databases. In Denmark, Norway, and Sweden, incidence and prevalence estimates ranged from 2.47 to 3.54 per 100,000 person-years and 8.0 to 8.7 per 100,000 persons, respectively.34,–,40 Last, we found similar prevalence estimates using capture-recapture methodology and a simulation-based approach using data from independent sources, thus further confirming the accurate case ascertainment of the data sources used.
Using the national DPR database enabled us to assess that 72.0% of all Dutch patients with MND were identified by the Netherlands MND registry. This coverage rate is comparable to or higher than other large population-based registries.41,–,43 The higher mean age of patients and higher proportion of female patients in the DPR database suggest that older female patients may be underrepresented in the Dutch population–based MND registry. This unintended selection of younger male patients with, on average, longer survival was also observed when trial participants were compared to the eligible MND population.17 It is important to recognize this unintended selection of patients because it can lead to biased effect sizes in population-based studies when investigating etiologic risk factors or genes.29,44 To resolve unintended selection in population-based studies, future studies that investigate reasons for nonparticipation in observational studies are needed to identify these as-yet unknown patient factors.
We found a national 53.2% increase in MND mortality rates from 1998 through 2017 in both men and women. This large increase in the last 2 decades underlines the impact that the care of patients with MND will have on health services in the future. There was geographic variation in MND risk, which could be related to local genetic admixture or complex environmental risk factors. Associating spatial risk with possible MND risk factors could stimulate etiologic research in MND.
Study Funding
This study was funded by the Netherlands ALS Foundation.
Disclosure
The authors report no disclosures relevant to the manuscript. Go to Neurology.org/N for full disclosures.
Appendix Authors

Footnotes
Go to Neurology.org/N for full disclosures. Funding information and disclosures deemed relevant by the authors, if any, are provided at the end of the article.
↵* These authors contributed equally to this work.
The Article Processing Charge was funded by The Netherlands ALS Foundation.
- Received July 29, 2020.
- Accepted in final form October 26, 2020.
- Copyright © 2021 The Author(s). Published by Wolters Kluwer Health, Inc. on behalf of the American Academy of Neurology.
This is an open access article distributed under the terms of the Creative Commons Attribution-NonCommercial-NoDerivatives License 4.0 (CC BY-NC-ND), which permits downloading and sharing the work provided it is properly cited. The work cannot be changed in any way or used commercially without permission from the journal.
References
- 1.↵
- Al-Chalabi A,
- Fang F,
- Hanby MF, et al
- 2.↵
- 3.↵
- 4.↵
- Chio A,
- Mora G,
- Moglia C, et al
- 5.↵
- Moore A,
- Young CA,
- Hughes DA
- 6.↵
- Spencer PS,
- Nunn PB,
- Hugon J, et al
- 7.↵
- 8.↵
- Olsen SF,
- Martuzzi M,
- Elliott P
- 9.↵
- Huisman MH,
- de Jong SW,
- van Doormaal PT, et al
- 10.↵
- 11.↵
- 12.↵
- Chapman DG
- 13.↵
- Deenen JC,
- Arnts H,
- van der Maarel SM, et al
- 14.↵
- van Eijk RPA,
- Nikolakopoulos S,
- Roes KCB, et al
- 15.↵
- 16.↵
- 17.↵
- van Eijk RPA,
- Westeneng HJ,
- Nikolakopoulos S, et al
- 18.↵
- 19.↵
- Bennett DA,
- Brayne C,
- Feigin VL, et al
- 20.↵
- 21.↵
- Sorenson EJ,
- Stalker AP,
- Kurland LT,
- Windebank AJ
- 22.↵
- 23.↵
- 24.↵
- 25.↵
- Visser AE,
- D'Ovidio F,
- Peters S, et al
- 26.↵
- 27.↵
- Huisman MH,
- Seelen M,
- van Doormaal PT, et al
- 28.↵
- Visser AE,
- Rooney JPK,
- D'Ovidio F, et al
- 29.↵
- 30.↵
- 31.↵
- Brouwer M,
- Huss A,
- van der Mark M, et al
- 32.↵
- 33.↵
- 34.↵
- Vieira VM,
- Hansen J,
- Gredal O,
- Weisskopf MG
- 35.↵
- Kioumourtzoglou MA,
- Seals RM,
- Himmerslev L,
- Gredal O,
- Hansen J,
- Weisskopf MG
- 36.↵
- 37.↵
- Nakken O,
- Lindstrom JC,
- Tysnes OB,
- Holmoy T
- 38.↵
- Longinetti E,
- Mariosa D,
- Larsson H, et al
- 39.↵
- Longinetti E,
- Regodon Wallin A,
- Samuelsson K, et al
- 40.↵
- 41.↵
- Zhou S,
- Qian S,
- Li X,
- Zheng L,
- Chang W,
- Wang L
- 42.↵
- Golby R,
- Poirier B,
- Fabros M,
- Cragg JJ,
- Yousefi M,
- Cashman N
- 43.↵
- Nelson LM,
- Topol B,
- Kaye W, et al
- 44.↵
Disputes & Debates: Rapid online correspondence
NOTE: All authors' disclosures must be entered and current in our database before comments can be posted. Enter and update disclosures at http://submit.neurology.org. Exception: replies to comments concerning an article you originally authored do not require updated disclosures.
- Stay timely. Submit only on articles published within the last 8 weeks.
- Do not be redundant. Read any comments already posted on the article prior to submission.
- 200 words maximum.
- 5 references maximum. Reference 1 must be the article on which you are commenting.
- 5 authors maximum. Exception: replies can include all original authors of the article.
- Submitted comments are subject to editing and editor review prior to posting.
You May Also be Interested in
Related Articles
- No related articles found.