ABSTRACT
Background REM sleep behavior disorder (RBD) is an early manifestation of alpha-synucleinopathy in many cases. Dream enactment behavior (DEB), the clinical hallmark of RBD, has many etiologies and cannot be used alone to predict underlying alpha-synucleinopathy. We compared the proportion of people with alpha-synucleinopathy, as measured by CSF alpha-synuclein seed amplification assay (CSFasynSAA), between people with polysomnographic-confirmed RBD (RBD-PSG) and people who reported DEB on a questionnaire and were further selected with smell testing and DAT-SPECT.
Methods Participants were enrolled in the Parkinson’s Progression Marker Initiative (PPMI) and ≥60 years old without a diagnosis of Parkinson’s disease. Participants had either RBD-PSG or self-reported DEB. Self-reported DEB participants had to have hyposmia (<10th percentile for age/sex) and at least mild DAT-SPECT abnormality (<100% age/sex-expected). We compared CSFasynSAA between RBD-PSG and self-reported DEB with hyposmia (DEB+Hypos). RBD-PSG participants also underwent smell testing and DAT-SPECT; we determined the predictive value of these tests in RBD-PSG with regards to CSFasynSAA.
Results CSFasynSAA was positive in 171/240 (71%) of RBD-PSG and %) 180/210 (86%) of DEB+Hypos participants. Among RBD-PSG, hyposmia strongly predicted CSFasynSAA+ (PPV: 92% [95% CI 87%-97%]). Smell identification was more accurate than DAT-SPECT in predicting CSFasynSAA+ in RBD-PSG (AUC for UPSIT: 0.89 [95% CI 0.84 – 0.94]; AUC for DAT-SPECT: 0.65 [95% CI 0.58 – 0.73]).
Conclusions Smell testing may be an effective and scalable method to identify people with alpha-synucleinopathy among those with self-reported DEB. Among individuals with RBD diagnosed by PSG, smell testing improved prediction of positive CSF alpha-synuclein biomarker.
INTRODUCTION
Identifying individuals with early alpha-synucleinopathy is of high interest as clinical trials aiming to test interventions to reduce disability from these diseases are in the planning stages.1 Several non-motor and motor features may be present in early stages of alpha-synucleinopathy, but many are non-specific. Isolated REM sleep behavior disorder (iRBD) is a parasomnia defined by clinical manifestations of dream enactment behavior (DEB) and by absence of atonia during REM sleep on polysomnogram (PSG) in individuals without a diagnosis of PD or other neurodegenerative disorders.2 As many as 90% of people with PSG confirmed iRBD have misfolded alpha-synuclein (asyn) in their cerebrospinal fluid (CSF),3–7 and iRBD is now widely accepted as early-stage alpha-synucleinopathy.8–11 Identifying people with iRBD is a high priority to establish clinical trial-ready cohorts.12 13
However, iRBD symptoms overlap with other disorders (e.g. sleep and psychiatric disorders) that are more prevalent and not due to neurodegeneration. The PSG required for definitive diagnosis of iRBD14 is not practical for large-scale screening. DEB is common – occurring in up to 19% of the general population15 – and can be identified in the community based on self-reported questionnaires. Yet DEB can result from other medical or psychiatric conditions and is only due to PSG confirmed RBD in a small subset.16 More accurate methods to identify underlying alpha-synucleinopathy in people with DEB are needed.
We therefore sought to investigate, via a tiered screening paradigm, how other assessments, specifically olfactory testing and dopamine transporter imaging, could be combined with DEB to accurately identify underlying alpha-synucleinopathy. Olfactory dysfunction is highly prevalent in neuronal alpha-synucleinopathy and is associated with disease progression in iRBD.3 17–20 Olfactory testing has been used to identify individuals with dysfunction in dopaminergic pathways as assessed by imaging.21 Among people with iRBD, dopamine imaging abnormalities predict development of parkinsonism.10 Determining how self-report of DEB, followed by smell testing and dopamine imaging, compares to PSG in identifying individuals with iRBD would inform large-scale screening strategies that seek to identify participants eligible for clinical trials.
METHODS
This was a cross-sectional analysis of individuals enrolled in Parkinson’s Progression Markers Initiative (PPMI), an international multisite longitudinal biomarker study.
Sample
Participants are recruited for this study either in-clinic or through remote recruitment strategies, including media campaigns.22 23 Inclusion and exclusion criteria are shown in Table 1. All participants had to be at least 60 years of age or older without a clinical diagnosis of Parkinson’s disease to be eligible.
In-clinic recruitment involved identification by study site investigators of individuals with RBD based on REM sleep without atonia on a PSG report (RBD-PSG). RBD-PSG were eligible to enroll in the study regardless of other features (iRBD). We additionally analyzed a small group of participants who were identified by investigators as having DEB thought to be likely due to RBD but without documented PSG confirmation (possible RBD, pRBD). People with pRBD also undergo smell testing with the UPSIT and dopamine transporter DAT-SPECT. People with pRBD who have an UPSIT < 10th %ile and DAT-SPECT binding less than 100% expected for age and sex (determined based on the lowest putamen specific binding ratio) are eligible for enrollment.
People recruited remotely go through a staged-screening pathway previously described.22 Briefly, participants were identified either through PPMI Online – an online longitudinal study involving participant reported outcomes in people with and without PD24 – or Smell Test (ST) Direct – a series of recruitment campaigns that sent out smell tests directly to participants. Participants were included in this analysis if they completed screening questionnaires (Supplementary Table 1) and reported a prior diagnosis of RBD made by a health professional or denied a diagnosis of RBD but reported DEB based on a single item question.25 Participants were mailed a University of Pennsylvania Smell Identification Test (UPSIT) to complete independently online. If the UPSIT score was < 10th percentile (%ile) compared to age- and sex-stratified normative values, participants were contacted for additional phone screening by a centralized team to confirm eligibility to undergo dopamine transporter imaging at a PPMI site using [123I]FP-CIT SPECT (DAT-SPECT).26–28 The lowest putamen specific binding ratio (SBR) was determined and compared to age- and sex-adjusted expected values. Participants with a lowest putamen SBR < 100% expected for age and sex were then invited to enroll. Because the threshold for UPSIT was much more stringent than the threshold for DAT-SPECT, where indeterminate deficits between 75% and 100% expected are less specific neurodegenerative synucleinopathy,29 30 we refer to this group collectively as DEB+Hypos.
Assessments
Enrolled participants underwent clinical assessments by PPMI site study personnel. Assessments included in this analysis include the Movement Disorders Society Unified Parkinson’s Disease Rating Scale (MDS-UPDRS),31 the REM sleep behavior disorder questionnaire (RBDSQ),32 and the Scales of Outcome of Parkinson’s Disease (SCOPA-Aut).33 Even though it was not used in selection criteria, participants with RBD-PSG also completed UPSIT and DAT-SPECT. Participants also underwent serum and cerebrospinal fluid (CSF) collection. CSF was sent for alpha-synuclein seed amplification assay (CSFasynSAA) through Amprion Inc.6 7 34 The technical document describing the assay is available for download in the PPMI database (www.ppmi-info.org/access-data-specimens/download-data). Briefly, the aSynSAA had a dual output: one for the detection or not-detection of synuclein seeds (positive, negative, inconclusive) and another for the type of seeds detected (Type 1, Type 2, and undetermined). Positive samples with Type 1 seeds presented high fluorescence values (≥45,000 RFU) and these seeds are predominantly found in participants with neuronal synuclein inclusions, while Type 2 seeds presented intermediate fluorescence values (≥3,000 RFU & <45,000 RFU) and these seeds are predominantly found in patients with underlying glial cytoplasmic inclusions and/or clinical presentation of multiple system atrophy (MSA).35
Statistical analysis
We included participants who had enrolled and completed a baseline visit in PPMI with CSFasynSAA available. We first used descriptive statistics to compare the demographic characteristics and prevalence of biomarkers (UPSIT and DAT-SPECT imaging result) of participants in each group (iRBD, pRBD, and DEB+Hypos). Univariate statistics were estimated using Chi-square or Fisher’s Exact tests for categorical variables and t-tests for continuous variables. We then determined the proportion of people with positive Type 1 CSFasynSAA in each group. We separately report the proportion with positive Type 2 CSFasynSAA. Among people with RBD-PSG, we determined the association between either an UPSIT < 10th %ile or a lowest putamen SBR < 75% expected for age and sex, based on the level of binding reduction observed in Parkinson’s disease,29 30 with positive Type 1 CSFasynSAA. We then used multivariate logistic regression to determine the association between age, sex, UPSIT, and DAT- SPECT and positive Type 1 CSFasynSAA. Finally, we constructed receiver operating characteristics (ROC) curves and calculated the area under the ROC curve (AUC) to determine the accuracy of UPSIT and DAT-SPECT in identifying aSyn-SAA.
RESULTS
Participant selection and characteristics
From in-clinic study site recruitment, we included 240 RBD-PSG and 61 pRBD participants. Through remote recruitment, 44,865 participants at least 60 years old completed an UPSIT remotely and answered screening questions related to RBD and DEB; 1,402 (3%) endorsed an RBD diagnosis while 5,345 (12%) denied an RBD diagnosis but endorsed DEB (Figure 1). Out of these 6,747 participants, 1,907 (28%) had an UPSIT < 10th %ile and 441 out of 557 (79%) who had completed DAT-SPECT had a lowest putamen SBR less than 100% expected for age and sex. While many participants are still in the process of enrolling, as of the time of data download, 210 DEB+Hypos participants had fully enrolled and had CSFasynSAA available for analysis.
Characteristics between RBD-PSG and DEB+Hypos are shown in Table 2; pRBD is shown in Supplementary Table 1. While age was similar in all groups, a greater proportion of the RBD-PSG group were male compared to the DEB+Hypos group.
Differences in UPSIT and DAT-SPECT across groups
As expected based on inclusion criteria, the average UPSIT percentile in the DEB+Hypos group was lower than the RBD-PSG group (Table 2). DAT-SPECT was lower for the same reason, though a similar proportion had only mildly reduced DAT-SPECT binding (between 75% and 100% expected for age and sex) across groups (Table 2).
Prevalence of CSFasynSAA positivity
Among the group recruited in-clinic, 171/240 (71%) of RBD-PSG had positive CSFasynSAA, whereas 54/61 (89%) of the pRBD group were positive. Among the remotely recruited DEB+Hypos group, 180/210 (86%) were positive (Table 2).
Determinants of CSFasynSAA in RBD-PSG
Given the broad range of smell test results in the RBD-PSG group, we stratified the group by a smell test cutoff of the 10th %ile, our inclusion criteria for other groups. Out of 240 RBD-PSG participants, 139 (58%) had an UPSIT < 10th %ile. DAT-SPECT did not significantly differ between the RBD-PSG participants with and without hyposmia (Table 3). The proportion of people with positive CSFasynSAA increased: among participants with RBD-PSG and hyposmia, 92% (128/139) had Type 1 positive CSFasynSAA compared to 43% (43/101) among participants with RBD-PSG without hyposmia.
We then determined the association between UPSIT, DAT-SPECT, and CSFasynSAA status (Figure 2; Supplementary Table 3). Among the 171 RBD-PSG participants with Type 1 positive CSFasynSAA, 128 (75%) had UPSIT < 10th %ile and 85 (50%) had lowest putamen SBR < 75% expected for age and sex on DAT-SPECT. An UPSIT < 10th %ile had a PPV of 92% (95% CI 88% – 97%) and 12.3 (95% CI 7.1 – 21.2, p<0.001) times increased adjusted odds of positive CSFasynSAA. DAT-SPECT < 75% expected for age and sex had a PPV of 82% (95% CI 74% - 89%) and 1.7 (95% CI 1.0 – 2.8, p = 0.05) times increased adjusted odds of having positive Type 1 CSFasynSAA (Table 4). When used as continuous measures among people with RBD-PSG to predict CSFasynSAA status, UPSIT showed an AUC of 0.89 (95% CI 0.84 – 0.94) compared to 0.65 (95% CI 0.58 – 0.73) for DAT-SPECT imaging results (Figure 3).
Clinical differences between participants with positive CSFasynSAA across groups
To determine whether different recruitment strategy and inclusion criteria may lead to selection for certain demographic or clinical features, we compared characteristics across participants from different groups with positive Type 1 CSFasynSAA. Demographics were similar except a higher proportion of females were seen in the DEB+Hypos group (Table 5). The MDS-UPDRS Part III was higher, though not by a large magnitude, in DEB+Hypos compared to RBD-PSG potentially due to the DAT-SPECT selection requirement. The RBDSQ was highest in the RBD-PSG with UPSIT < 10th %ile group. Otherwise, no substantial clinical differences were observed. Characteristics of CSFasynSAA+ pRBD participants are shown in Supplementary Table 4.
DISCUSSION
In this study we have investigated, among individuals with dream enactment behavior, different methods of identifying early alpha-synucleinopathy. A key finding was that a questionnaire and remotely completed smell test can be used to identify a high proportion of people with positive CSFasynSAA. In our study, the proportion of people with positive CSFasynSAA was higher than in participants with RBD-PSG who were not selected based on smell. PSG remains the gold standard for diagnosing RBD but is costly36 and requires access to specialist sleep centers for proper administration and interpretation, which is less available in underserved populations.37 Questionnaires about DEB alone are not sufficient: this strategy can identify iRBD when used in people presenting to sleep clinics,25 but have low positive predictive value when used in the community setting.16 In our study, a brief questionnaire followed by a smell test and identified a population with 86% positivity for misfolded CSF alpha-synuclein. While the individuals who were identified by the questionnaire strategy could not be said to meet current diagnostic criteria for iRBD, which depends on PSG confirmation, they do have clinical and biological features consistent with early alpha-synucleinopathy. Moreover, our goal is not to make a clinical diagnosis of iRBD, but rather to identify individuals with early alpha-synucleinopathy based on positive CSFasynSAA status who are at risk for disease progression. Within this framework, the questionnaire/olfactory testing strategy is both more efficient and likely less costly than PSG.
While our algorithm included DAT-SPECT prior to SAA testing, our data further suggest that the screening protocol we used will likely be just as effective without the DAT-SPECT imaging component, which is still costly and burdensome to participants. This hypothesis is supported by our second key finding: among people with RBD-PSG, hyposmia is a strong determinant of CSFasynSAA (AUC: 0.89) compared to DAT-SPECT (AUC: 0.65). This finding is consistent with prior studies of people with RBD-PSG demonstrating more prominent hyposmia in positive compared to negative CSFasynSAA38 and identifying hyposmia as a risk factor for future clinical diagnosis of PD.3 10 We therefore expect hyposmia alone to be a sufficient screen to identify early alpha-synucleinopathy among people with self-reported DEB or RBD, but future confirmatory studies are needed. Questionnaires in combination with smell testing, which can be remotely and broadly applied, would be an efficient and inclusive recruitment approach for clinical trials seeking to test therapies in people with early alpha-synucleinopathy.
Notably, the proportion of RBD-PSG participants with positive CSFasynSAA in our study (71%) was lower than what has been reported in some studies,3 6 where cohorts of PSG-confirmed RBD have been ∼90% positive for CSFasynSAA. This difference may be because of broader inclusion in our RBD-PSG cohort compared to other studies, leading this cohort to be more heterogenous than other cohorts. In PPMI, an investigator confirms RBD based on PSG report, but full criteria may vary from center to center. Formal PSG confirmation requires stringent criteria2 which may improve the specificity for underlying synucleinopathy. Arguably, the criteria in PPMI are more generalizable to clinical situations and easier to implement for trials, therefore showing what could be expected in terms of enrichment for CSFasynSAA. Other recent studies have reported similar proportions of CSFasynSAA in RBD-PSG.38 39 Importantly, even in our cohort of PSG- confirmed RBD, hyposmia significantly enriched for CSFasynSAA. UPSIT is a readily available inexpensive test that can be used broadly in sleep clinics as a proxy for likely underlying asyn pathology. However, therapeutic studies aiming to enroll iRBD participants specifically testing asyn targeting therapies will need to establish a paradigm for biomarker testing.
Our observation that DAT-SPECT was not as strongly associated with CSFasynSAA as smell supports the hypothesis that hyposmia comes prior to substantial striatal neurodegeneration. Most people with RBD-PSG and positive CSFasynSAA did not have substantial reduction in DAT-SPECT binding: ∼50% had normal or only mild reductions in DAT-SPECT binding (lowest putamen SBR > 75% expected for age and sex). Longitudinal studies will test our expectation that DAT-SPECT will continue to decrease in these participants. Still, in this early stage, hyposmia appears more useful than DAT-SPECT at identifying people with underlying asyn pathology who do not yet have neurodegeneration.
Recent validation of in vivo alpha-synuclein biomarkers4 7 40 41 has enabled a biologic definition of Neuronal Synuclein Disease (NSD),8 9 a new terminology encompassing clinical syndromes of Parkinson’s disease (PD) and dementia with Lewy Bodies (DLB) as well as pre-diagnostic syndromes, to describe individuals with neuronal a-synucleinopathy (n-asyn). A staging system, the NSD Integrated Staging System (NSD-ISS), was also put forth, defined based on a combination of biomarkers, clinical signs and symptoms, and their functional consequences.42
People with biomarkers of n-asyn and possibly dopaminergic transporter (DAT) dysfunction without signs or symptoms are NSD-ISS Stage 1. Stage 2 is comprised of people with positive n- asyn only (Stage 2A), additional evidence of dopaminergic dysfunction (Stage 2B), and presence of signs or symptoms but without meaningful functional impairment. Many individuals with iRBD are in stage 2 NSD,42 representing a desirable cohort for intervention to prevent the onset and progression of disabling symptoms. Since n-asyn positivity is a mandatory criterion for NSD and for application of the NSD-ISS, identifying strategies that efficiently detect individuals with n-asyn positivity is crucial, and individuals with clinical history or self-report of DEB offer the prime opportunity to develop these methods. As demonstrated in this study, combining DEB with hyposmia increased the yield for identifying individuals with NSD.
Limitations
A limitation of this study is the lack of available longitudinal data. Longitudinal data will determine how these groups may differ in the rate and type of symptoms that arise, further informing these recruitment strategies. These studies will be the focus of future analyses. Additionally, the exclusion of participants with UPSIT > 10th %ile and DAT-SPECT > 100% expected in the pRBD and DEB+Hypos groups prevented us from understanding the full test characteristics of these assessments in our groups. Further information on normosmic participants with a wide range of dopaminergic function, with and without self-reported RBD/DEB, would be useful to understand these associations more completely. This selection limits the interpretation of clinical differences between these cohorts in our report. However, the goal of this work was to demonstrate the predictive value of a recruitment protocol implemented in the real world, which can be easily replicated for ongoing clinical trials targeting people with early stage NSD.
In conclusion, our study demonstrates that self-reported DEB combined with screening for severe hyposmia can identify individuals with high probability of having n-asyn pathology. This offers a scalable inexpensive approach for population screening to identify individuals in early stages of NSD. Further, assessment of sense of smell even in individuals with PSG-confirmed RBD, the current diagnostic gold standard, may be more useful than DAT-SPECT to identify people with n-asyn pathology. Smell testing could therefore provide a scalable, safe, and low-cost strategy to identify people with early stage NSD for clinical trial recruitment.
Funding
Funding support for the data analysis was provided by The Michael J. Fox Foundation for Parkinson’s Research (MJFF). PPMI – a public-private partnership – is funded by MJFF and funding partners, including 4D Pharma, Abbvie, AcureX, Allergan, Amathus Therapeutics, Aligning Science Across Parkinson’s, AskBio, Avid Radiopharmaceuticals, BIAL, BioArctic, Biogen, Biohaven, BioLegend, BlueRock Therapeutics, BristolMyers Squibb, Calico Labs, Capsida Biotherapeutics, Celgene, Cerevel Therapeutics, Coave Therapeutics, DaCapo Brainscience, Denali, Edmond J. Safra Foundation, Eli Lilly, Gain Therapeutics, GE HealthCare, Genentech, GSK, Golub Capital, Handl Therapeutics, Insitro, Janssen Neuroscience, Jazz Pharmaceuticals, Lundbeck, Merck, Meso Scale Discovery, Mission Therapeutics, Neurocrine Biosciences, Neuropore, Pfizer, Piramal, Prevail Therapeutics, Roche, Sanofi, Servier, Sun Pharma Advanced Research Company, Takeda, Teva, UCB, Vanqua Bio, Verily, Voyager Therapeutics, the Weston Family Foundation and Yumanity Therapeutics.
Data Availability Statement
Data used in the preparation of this article were obtained on July 8th 2024 from the Parkinson’s Progression Markers Initiative (PPMI) database (www.ppmi-info.org/access-data-specimens/download-data), RRID:SCR_006431. For up-to-date information on the study, visit www.ppmi-info.org. This analysis was conducted by the PPMI Statistics Core and used actual dates of activity for participants, a restricted data element not available to public users of PPMI data. This analysis used DAT-SPECT and aSynSAA results for participants of the Prodromal Cohort, obtained from PPMI upon request after approval by the PPMI Data Access Committee. Protocol information for The Parkinson’s Progression Markers Initiative (PPMI) Clinical - Establishing a Deeply Phenotyped PD Cohort AM 3.2. can be found on protocols.io or by following this link: https://dx.doi.org/10.17504/protocols.io.n92ldmw6ol5b/v2.
Code Availability Statement
Statistical analysis codes used to perform the analyses in this article are shared on Zenodo (10.5281/zenodo.13951754)
Competing Interests
EB declares grants to his institution from MJFF, Biogen (clinical trial funding), National Institutes of Health, the Department of Defense, and the Gateway Institute for Brain Research. EB declares consultancies for Guidepoint Inc, Rune Labs Inc, and a member of the scientific advisory board of 153 Therapeutics. LC declares grants to her institution from Biogen (clinical trial funding), MJFF, UPMC Competitive Medical Research Fund, National Institutes of Health, and University of Pittsburgh; grant and travel support from MJFF; royalties from Wolters Kluwel (for authorship); and in-kind donation by Advanced Brain Monitoring of equipment for research study to her institution. AV has nothing to disclose. MM declares research funding to her institution from The Michael J. Fox Foundation. RK declares research funding to her institution from The Michael J. Fox Foundation. CC-G declares research funding to her institution from The Michael J. Fox Foundation. MB declares travel grants from The Michael J. Fox Foundation. CCo declares grants from The Michael J. Fox Foundation and NIH/NINDS. AS declares consultancies for Mitsubishi, GE healtcare, Capsida Therapeutics and Parkinson Study Group; grants from The Michael J. Fox Foundation (member of PPMI Steering Committee); and participation on DSMB boards at Spark Therapeutics, Cerevance. Alerity, Wave Life Sciences, Inhibikase, Prevai (Eli Lilly), Huntington Study Group and Massachusetts General Hospital. TS declares consultancies for AcureX, Adamas, AskBio, Amneal, Blue Rock Therapeutics, Critical Path for Parkinson’s Consortium, Denali, The Michael J. Fox Foundation, Neuroderm, Roche, Sanofi, Sinopia, Takeda, and Vanqua Bio; on advisory boards for AcureX, Adamas, AskBio, Biohaven, Denali, GAIN, Neuron23 and Roche; on scientific advisory boards for Koneksa, Neuroderm, Sanofi and UCB; and received research funding from Amneal, Biogen, Roche, Neuroderm, Sanofi, Prevail and UCB and an investigator for NINDS, MJFF, Parkinson’s Foundation. KM declares support to his institution (Institute for Neurodegenerative Disorders) from The Michael J. Fox Foundation. KM also declares consultancies for Invicro, The Michael J. Fox Foundation, Roche, Calico, Coave, Neuron23, Orbimed, Biohaven, Anofi, Koneksa, Merck, Lilly, Inhibikase, Neuramedy, IRLabs and Prothena and participates on DSMB at Biohaven. CT declares consultancies for CNS Ratings, Australian Parkinson’s Mission, Biogen, Evidera, Cadent (data safety monitoring board), Adamas (steering committee), Biogen (via the Parkinson Study Group steering committee), Praxis (via the International Parkinsons and Movement Disorder Society), Kyowa Kirin (advisory board), Lundbeck (advisory board), Jazz/Cavion (steering committee), Acorda (advisory board), Bial (DMC) and Genentech. CT also declares grant support to her institution from The Michael J. Fox Foundation, National Institute of Health, Gateway LLC, Department of Defense, Roche Genentech, Biogen, Parkinson Foundation and Marcus Program in Precision Medicine. CT declares membership on the NPJ Parkinson’s Disease Editorial Board.
AUTHOR CONTRIBUTIONS
EB, AV, LC, AS, TS, KM, and CT contributed to the conceptualization. EB, AV, LC, CCo, AS, TS, KM, and CT contributed to the investigations. AV, LC, MD, CCo, AS, TS, KM, and CT contributed to the supervision. EB, MM, RK, CC-G, MB, and CCo contributed to the data curation and formal analysis. EB, LC, MM, RK, CC-G, MB, and CCo contributed to the validation. EB wrote the original draft. EB, AV, LC, CCo, AS, TS, KM, and CT provided project administration. EB, CCo, AS, TS, KM, and CT acquired funding. CCo, AS, TS, KM, and CT provided resources. All authors contributed to the methodology, reviewed and edited the final version, and had direct access to and verified the underlying data.