- 1Department of Neurology, University Hospital Dresden, Dresden, Germany
- 2Department Klinische Forschung, University Hospital Basel, Basel, Switzerland
- 3Department of Neurology, University Hospital Basel, Basel, Switzerland
Purpose: Neurofilament light chain in serum (sNfL) has been suggested as a biomarker for the assessment of neuroaxonal damage. Since NfL are not expressed in muscle, elevated sNfL in patients with primary myopathies suggest additional nervous system involvement. To verify this hypothesis, we measured sNfL in a series of patients with myopathies.
Methods: sNfL were determined in 62 patients with molecular proven primary myopathies in whom some nervous system involvement may be predicted: myotonic dystrophy type I and II (DM I, II) and mitochondrial disease. In addition, sNfL were measured in 8 patients with facioscapulohumeral muscular dystrophy (FSHD) and in a disease control group caused by genetic defects exclusively expressed in muscle.
Results: sNfL values were significantly elevated in the DM I, the DM II and the mitochondrial group, with FSHD patients showing the lowest sNfL elevations. sNfL levels in the disease control group were not different from the healthy controls. A significant correlation between repeat length and sNfL levels was found in the DM I patients, but not in the DM II patients. Mitochondrial patients with encephalopathy showed significantly higher sNfL concentrations compared to patients with only muscular symptoms.
Conclusion: sNfL levels are elevated in myopathies with, based on the underlying molecular defect or clinical features, established nervous system involvement, i.e., myotonic dystrophies and mitochondrial disorders. sNfL were also raised in FSHD, where involvement of the nervous system is not usually clinically apparent. Thus, sNfL concentrations may serve as a biomarker for additional neuronal damage in primary myopathies.
Introduction
Primary myopathies are a diverse group of genetic and acquired conditions, which may present with a plethora of clinical symptoms, and it is becoming increasingly evident that many hereditary primary myopathies can be ascribed to underlying multisystem disorders. MRI studies have uncovered central nervous system (CNS) involvement, with evidence of cerebral atrophy, gray and white matter lesions and ventricular enlargement in a number of muscle diseases (Angelini and Pinzan, 2019). These diseases include certain muscular dystrophies (e.g., Duchenne), the dystrophic myotonias and mitochondrial myopathies.
The myotonic disorders are classified into dystrophic (myotonic dystrophy Type I/DM1 and myotonic dystrophy Type II/DM 2/PROMM) and non-dystrophic myotonias. The underlying genetic defect in the myotonic dystrophies affects transcription and translation of multiple cellular genes, thus also causing extramuscular pathology, with cardiac, endocrinological and CNS-symptoms.
The mitochondrial diseases represent the archetype of multisystem disorders because of the ubiquitous occurrence of mitochondria, but they also display a high degree of heterogeneity due to the involvement of either the nuclear or mitochondrial genome and, in the case of mitochondrial DNA mutations, frequent disparity in the level of mutations between the different tissues. Not surprisingly, because of their high energy requirements, both the peripheral and the central nervous system are frequently involved in mitochondrial diseases. Our cohort included patients with proven defects of nuclear genes (POLG, encoding the mitochondrial polymerase gamma and DNAJC30, encoding DnaJ heat shock protein family (Hsp40) member C30) and of mtDNA encoded mitochondrial genes (MT-ATP6, encoding ATPase 6, MT-TL1, encoding tRNA Leucine 1, MT-TK encoding tRNA Lys, MT-T1 encoding tRNA Ile) and patients with a single, large scale deletion in mtDNA.
Clinically evident CNS involvement is rare in facioscapulohumeral muscular dystrophy (FSHD), but detailed neurophysiological testing has shown subclinical involvement in a number of cases (Stübgen, 2007).
Until now, the assessment of CNS involvement in myopathies has mainly been based upon clinical parameters such as neuropsychological testing or cerebral imaging by MRI. Therefore, there is a need for reliable and easily accessible biomarkers for the detection of relevant nervous system damage, possibly even at a pre-symptomatic stage.
Neurofilaments constitute an important part of the neuronal cytoskeleton and are relevant for axonal growth and transport. Elevated levels of serum neurofilament light chain have been shown to specifically reflect neuroaxonal damage in a variety of neurodegenerative conditions such as amyotrophic lateral sclerosis (Steinacker et al., 2016), Alzheimer’s dementia (Mattsson et al., 2017) and multiple sclerosis (Kuhle et al., 2017). Plasma or serum NfL have also been measured in peripheral nervous system disorders, including Charcot-Marie-Tooth disease (Sandelius et al., 2018), the acquired polyneuropathies (Mariotto et al., 2018; Van Lieverloo et al., 2019) and Guillain-Barré syndrome (Martín-Aguilar et al., 2021). In all cases, increased NfL levels were associated with greater disease activity and severity, thus monitoring of serum NfL levels may also be used to assess disease progression.
The aim of this study was to ascertain whether, using sNfL as a non-invasive, easily accessible biomarker, nervous system involvement can be demonstrated in patients with primary myopathies. We chose three multisystem disorders, DM1, DM2/PROMM, and mitochondrial disease, with frequent nervous system involvement and one myopathy without overt nervous system involvement (FSHD) to test this hypothesis. In addition, we measured sNfL concentrations in a group of disease controls, consisting of monogenic primary myopathies caused by pathogenic mutations in genes which are expressed in muscle tissue only. These comprise structural muscle proteins, channel proteins of the muscle membrane and proteins related to muscle energy metabolism.
Methods
Clinical data and blood samples were prospectively collected in the Neuromuscular Outpatient Clinic of the University Hospital Dresden, Germany. Molecular genetic confirmation of the diagnosis was available for all patients. Patients with known acquired diseases or injuries of the CNS were excluded from the study. The NfL study was approved by the local Ethics committee (ID: EK394102018), conforms with World Medical Association Declaration of Helsinki and the patients or their authorized representatives consented to the use and storage of the biosamples and publication of their data.
All serum samples were frozen immediately and stored at −20°C. sNfL concentrations were determined using the single-molecule array (SIMOA) analysis as previously described (Akgün et al., 2019; Sutter et al., 2020). The mean intra-assay coefficient of variation of duplicates and the mean inter-assay coefficient of variation was<10%. The control group comprised 485 samples from healthy volunteers, which were collected in the University Hospital Basel as previously described (Sutter et al., 2020).
To account for the known positive correlation of NfL with age, age-corrected sNfL percentile values have been derived based on healthy controls using a Generalized Additive Model for Location, Scale and Shape (GAMLSS) as described previously (Sutter et al., 2020). The percentile value was then calculated for each data point. Unadjusted group differences in sNfL levels were visualized using boxplots and were assessed using the Wilcoxon rank sum test. To additionally account for the minor age differences between groups, a multivariable linear regression model with log (sNfL) as the dependent variable was built with group and age as predictors. Estimates were back-transformed (exponentiated) and therefore represent multiplicative effects on the geometric mean of sNfL. Correlation of sNfL with age and repeat length within a patient group were assessed using Pearson correlation coefficients. Differences between different groups of mitochondrial patients were assessed using the unpaired t-test (GraphPad Prism). All remaining analyses were done with the statistical Software R (version 4.1.0). We used the STARD reporting guidelines (Bossuyt et al., 2015).
Results
Between December 2019 and September 2021, 62 patients with either DM1, DM 2/PROMM, mitochondrial disease or FSHD and 13 patients for the disease control group were recruited for the study. Clinical characteristics and demographics of the patient cohorts are shown in Table 1.
Concomitant neurological diseases (depression, dissociative seizures, fibromyalgia syndrome, migraine, chronic/neuropathic pain, restless legs syndrome, essential tremor) in some patients with DM I/DM II and the medication taken by the patients are not, according to current knowledge, associated with structural neuronal damage and would therefore not affect sNfL levels. Several patients in the mitochondrial group showed, in addition to the myopathy, clinically apparent CNS-symptoms (Table 1). In these patients with mitochondrial myopathy, the CNS-symptoms are a feature of the primary mitochondrial disorder and therefore clinically confirm neuronal damage underlying the sNfL elevation.
All disease groups showed significantly higher median serum NfL levels than the healthy controls, whilst the smallest elevation of sNfL concentrations was found in the FSHD group (Figure 1). After adjustment for age, serum NfL was on average 2.1 times higher (CI: 1.6;2.7, p-value: < 0.0001) in DM I, compared to levels in healthy controls, 1.4 times (CI: 1.2;1.7, p-value: < 0.0001) in DM 2/PROMM, 1.4 times (CI: 1.1;1.9, p-value: 0.0227) in FSHD and 3.2 times (CI: 2.7;3.8, p-value: < 0.0001) in the Mito-patient cohort.
Figure 1. Boxplot diagram of serum neurofilament light chain (sNfL) concentrations, comparing the disease groups against healthy control samples. Boxes represent median and interquartile range (IQR) and whiskers the extreme value within 1.5 × IQR above and below the median. The p-values listed were determined using a Wilcoxon rank-sum test.
Compared to the controls, serum NfL levels were above the 90th percentile in 6/9 (66%) patients with DM I and in 9/22 (41%) patients with PROMM (Figure 2). In the FSHD group, only 3/8 (38%) patients were above the 90th percentile. The mitochondrial patients had the highest NfL levels of all groups, and 16/23 (70%) patients had values above the 90th percentile (Figure 2).
Figure 2. Serum neurofilament light chain (sNfL) concentrations in the four disease groups are shown in comparison to 485 control samples. The distribution of sNfL is shown as a function of age and expressed as percentile curves (coloured lines), based on the control samples (gray circles). Patient 18 from the mitochondrial group is not included in this figure as the sNfL value is outside the range of the axes.
Within each group, only the DM I and DM 2/PROMM patients showed a significant correlation between age and sNfL levels (r = 0.69, p = 0.04 and r = 0.83, p < 0.0001, respectively) (Supplementary Figures 1A,B), whilst no significant correlation was seen within the other groups (Supplementary Figures 1C–E).
No significant correlation was found between sNFL levels and repeat length in the DM2/PROMM patients (r = 0.5; p = 0.07), but there was a significant correlation in the DM I patients (r = 0.86; p = 0.014). P7 was omitted from the regression analysis as an outlier value, because the clinical symptoms (gait ataxia out of proportion to the myopathy), noted at the last examination when the sample was obtained, suggested additional neuronal pathology.
The mitochondrial cohort is a more heterogeneous group than the other groups, in terms of both clinical symptoms and underlying genetic cause. Within the mitochondrial group, those patients with only clinically apparent muscle involvement had lower sNFL levels than patients with additional non-muscle symptoms. Thus, patients with a single deletion in mtDNA and manifesting CPEO only, had significantly lower sNFL levels than patients with CPEO + and a single or multiple deletions (due to a mutation in POLG) in mtDNA in muscle (Table 1 and Supplementary Figure 2; 11.3 ± 2 vs. 24.6 ± 5.3; mean ± SD; p < 0.05). The greatest range in sNFL levels occurred in the group of patients with a mutation in mtDNA, where the levels ranged from 11.9 (patient 17 with a heteroplasmic m.8344A > G mutation in tRNA lysine) to 1,250 pg/ml (patient 18, with a homoplasmic m.4290T > C mutation in tRNA Ile) (Table 1 and Supplementary Figure 2). Four patients from the mitochondrial group harbored a heteroplasmic m.3243A > G mutation in tRNA leucine 1 (the so-called MELAS mutation, Table 1). Although this group is small, the sNFL levels in these 4 patients showed no correlation with the mutation load in either blood or urine (Pearson correlation test, p = 0.7 and 0.6, respectively). Interestingly, the two patients with less severe neuronal involvement (P11 and P14 with polyneuropathy only) had similar sNFL levels 13.7 and 12.7 pg/ml), despite the difference in the underlying genetic cause. P11 harbors a heterozygous mutation in the nuclear gene POLG, whilst P14 harbors a m.9185T > C in the mitochondrial gene MT-ATP6 which encodes a subunit of complex V.
Discussion
This study demonstrates that sNfL levels can be used as a sensitive biomarker of ongoing neuronal damage in primary myopathies. As proof of concept, sNFL levels were significantly elevated in patients suffering from myopathies in which CNS-manifestations have been established, either as a clinically apparent symptom of the underlying genetic defect or by auxiliary clinical investigations (MRI, neurophysiology). Thus, we found sNfL to be high in patients with DM1, PROMM and mitochondrial disease, but also—even though to a lesser degree—in patients with FSHD. Our concept of sNfL as a biomarker for nervous system involvement in primary myopathies is further validated by the lack of sNfL elevations in patients from the disease control group (myopathies due to mutations in genes exclusively expressed in muscle tissue).
The severity of neuronal damage appeared to be loosely correlated with the degree of sNfL elevation: disease groups with clinically common and more severe CNS-symptoms (DM1, mitochondrial myopathy) showed on average higher median sNfL levels than disorders in which signs of CNS-involvement are mild and less common (PROMM/FSHD).
Within the Mito-group, patients with a single deletion in mtDNA and only CPEO had on average, a significantly lower sNfL level than patients with CPEO plus (e.g., P6 with ataxia and P9 with sensorineural hearing loss and cognitive impairment). In addition, the highest sNfL levels were observed in patients with the most severe involvement, even when the underlying mutation resides in the same mtDNA gene. Thus, in the two patients with a homoplasmic mutation in the mtDNA gene MT-ATP6 (P14 and P15), the highest sNFL levels were detected in P15, with ataxia, dysarthria, epilepsy, optic atrophy, cognitive impairment and polyneuropathy as compared to P14 with polyneuropathy. The highest sNfL levels were observed in patients who experienced symptoms such as epilepsy, ataxia or stroke-like episodes (Table 1: P15, P16, P18, P19, P20). Patient 18 (described in Limongelli et al., 2004) had the highest sNFL levels of all of the mitochondrial patients (1,250 pg/ml), and she has the most severe neurological involvement of the group with Leigh disease with onset at the age of 16 years and necrotizing encephalopathy of basal ganglia and brain stem.
One of the greatest problems associated with reaching a prognosis in individuals who harbor a mutation in mtDNA is that these mutations are often heteroplasmic, i.e., both mutant and wild type mtDNA are present in an affected individual, and the mutation load may vary from tissue to tissue. The mutation level in urinary epithelial cells can be a useful predictor of the mutation load in muscle (Whittaker et al., 2009), but the extent to which the level in this tissue mirrors that in nervous tissue is unknown. Three of our four patients with the m.3243A > G mutation (tRNA leucine 1) harbor similar mutation loads in urinary epithelial cells, but show different degrees of neuronal involvement which suggests that the sNFL levels may be a better indicator of neuronal involvement than the degree of heteroplasmy in urinary epithelial cells. Thus, regular measurement of sNFL levels may provide an important prognostic tool for treatment of these patients. Indeed, higher sNfL concentrations in patients with MELAS-syndrome during acute stroke-like episodes than in the attack-free interval suggest the applicability of sNfL for the assessment of neuronal damage in mitochondrial diseases (Zheng et al., 2021). Also in our cohort, P20 had only manifested sensorineural hearing loss until the age of 61, when she developed stroke-like episodes with a very high sNFL level of 70.4 pg/ml.
Interestingly, only patients with DM I and DM2/PROMM showed a significant correlation of sNfL with age (Supplementary Figures 1A,B) within their cohort. This suggests that additional neuronal damage, above that which occurs with normal aging, accumulates in the myotonic dystrophies. Thus, our observations could support the concept of myotonic dystrophies as progeroid disorders with premature aging, as suggested previously (Meinke et al., 2018).
Regarding repeat length, sNfL correlated well to repeat length for the DM I patients, which is supported by previous observations of a negative correlation between brain volume on neuroimaging and CTG repeat length (van der Plas et al., 2019). However, one patient (P7) presented sNfL greatly above the range of all other DM I patients; this patient, as mentioned above, also exhibited clinical symptoms which indicated additional neuronal involvement separate from DM I.
In contrast, no correlation could be demonstrated between sNfL and repeat length in the DM2/PROMM patients. This is not surprising, however, considering the lack of a correlation between repeat length and the severity of the myopathy in these patients (Schoser, 2006).
Interestingly, even patients with FSHD demonstrated elevated sNfL compared to age matched controls, although CNS-involvement is not generally regarded to be clinically relevant in FSHD. However, neurophysiological (Stübgen, 2007) and cognitive (Sistiaga et al., 2009) assessment revealed alterations of CNS function in up to 70% of patients with FSHD.
Owing to the rarity of the myopathies examined in this study, the size of the disease groups was small, but despite this drawback, the estimated differences in sNfL levels were very high. Other confounding variables affecting NfL concentrations, such as a subclinical age-dependent vascular encephalopathy (Thebault et al., 2020), were accounted for by comparison with NfL values in a large group of healthy controls of a similar age range. Thus, NfL elevations due to an age-dependent encephalopathy in the patient groups would be matched in the age-adjusted controls. Other concomitant CNS-disease which might influence NfL levels, was excluded on clinical grounds.
We propose from this pilot study that sNfL could be used as simple and non-invasive biomarker to detect and monitor neuronal damage in myopathies, even in those without clinically evident CNS-involvement. sNfL might also be helpful to assess treatment effects in view of upcoming new gene therapies for myopathies with nervous system involvement, such as the dystrophinopathies, the myotonic dystrophies or mitochondrial disorders. Larger longitudinal studies comparing sNfL levels with imaging or neuropsychological course parameters, are required to elucidate the value of sNfL as an outcome parameter in clinical practice.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Ethics Statement
The studies involving human participants were reviewed and approved by the Local Ethics committee of the Technische Universität Dresden (ID: EK394102018). The patients/participants provided their written informed consent to participate in this study.
Author Contributions
AS, JS, and SJ designed the concept of the study. AS, JS, JK, PB, SJ, EW, TZ, and KA contributed to data acquisition and analysis. KA and TZ developed ethics protocol. AS, JS, SJ, and PB drafted the manuscript and figures. All authors contributed to the article and approved the submitted version.
Funding
Funding for the open access publication fee is granted by the Technische Universität Dresden.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2021.791670/full#supplementary-material
References
Akgün, K., Kretschmann, N., Haase, R., Proschmann, U., Kitzler, H. H., Reichmann, H., et al. (2019). Profiling individual clinical responses by high-frequency serum neurofilament assessment in MS. Neurol. Neuroimmunol. Neuroinflamm. 6:e555. doi: 10.1212/NXI.0000000000000555
Angelini, C., and Pinzan, E. (2019). Advances in imaging of brain abnormalities in neuromuscular diseases. Ther. Adv. Neurol. Disord. 12, 1–24.
Bossuyt, P. M., Reitsma, J. B., Bruns, D. E., Gatsonis, C. A., Glasziou, P. P., Irwig, L., et al. (2015). for the STARD Group, STARD 2015: an updated list of essential items for reporting diagnostic accuracy studies. Clin. Chem. 61, 1446–1452. doi: 10.1373/clinchem.2015.246280
Kuhle, J., Nourbakhsh, B., Grant, D., Morant, S., Barro, C., Yaldizli, O., et al. (2017). Serum neurofilament is associated with progression of brain atrophy and disability in early MS. Neurology 88, 826–831. doi: 10.1212/WNL.0000000000003653
Limongelli, A., Schaefer, J., Jackson, S., Invernizzi, F., Kirino, Y., Suzuki, T., et al. (2004). Variable penetrance of a familial progressive necrotising encephalopathy due to a novel tRNA(Ile) homoplasmic mutation in the mitochondrial genome. J. Med. Genet. 41, 342–349. doi: 10.1136/jmg.2003.016048
Mariotto, S., Farinazzo, A., Magliozzi, R., Alberti, D., Monaco, S., and Ferrari, S. (2018). Serum and cerebrospinal neurfilament light chain levels in patients with acquired peripheral neuropathies. J. Peripher. Nerv. Syst. 23, 174–177. doi: 10.1111/jns.12279
Martín-Aguilar, L., Camps-Renom, P., Lleixà, C., Pascual-Goñi, E., Díaz-Manera, J., Rojas-García, R., et al. (2021). Serum neurofilament light chain predicts long-term prognosis in Guillain-Barré syndrome patients. J. Neurol. Neurosurg. Psychiatry 92, 70–77. doi: 10.1136/jnnp-2020-323899
Mattsson, N., Andreasson, U., Zetterberg, H., and Blennow, K. Alzheimer’s Disease Neuroimaging Initiative (2017). Association of plasma neurofilament light with neurodegeneration in patients with alzheimer disease. JAMA Neurol. 74, 557–566. doi: 10.1001/jamaneurol.2016.6117
Meinke, P., Hintze, S., Limmer, S., and Schoser, B. (2018). Myotonic dystrophy - a progeroid disease? Front. Neurol. 9:601. doi: 10.3389/fneur.2018.00601
Sandelius, A., Zetterberg, H., Blennow, K., Adiutori, R., Malaspina, A., Laura, M., et al. (2018). Plasma neurofilament light chain concentration in the inherited peripheral neuropathies. Neurology 90, e518–e524. doi: 10.1212/WNL.0000000000004932et
Schoser, B. (2006). “Myotonic Dystrophy Type 2,” in GeneReviews, eds M. P. Adam, H. H. Ardinger, R. A. Pagon, S. E. Wallace, L. J. H. Bean, and G. Mirzaa (Seattle: University of Washington).
Sistiaga, A., Camaño, P., Otaegui, D., Ibáñez, B., Ruiz-Martinez, J., Martí-Massó, J. F., et al. (2009). Cognitive function in facioscapulohumeral dystrophy correlates with the molecular defect. Genes Brain Behav. 8, 53–59. doi: 10.1111/j.1601-183X.2008.00442.x
Steinacker, P., Feneberg, E., Weishaupt, J., Brettschneider, J., Tumani, H., Andersen, P. M., et al. (2016). Neurofilaments in the diagnosis of motoneuron diseases: a prospective study on 455 patients. J. Neurol. Neurosurg. Psychiatry 87, 12–20. doi: 10.1136/jnnp-2015-311387
Stübgen, J. P. (2007). Facioscapulohumeral muscular dystrophy. Multimodal evoked potentials and electroretinogram. Electromyogr. Clin. Neurophysiol. 47, 233–241.
Sutter, R., Hert, L., De Marchis, G., Twerenbold, R., Kappos, L., Naegelin, Y., et al. (2020). Serum neurofilament light chain levels in the intensive care unit: comparison between severely ill patients with and without covid-19. Ann. Neurol. 89, 610–616. doi: 10.1002/ana.26004
Thebault, S., Booth, R. A., and Freedman, M. S. (2020). Blood neurofilament light chain: the neurologist’s troponin? Biomedicines 8:523.
van der Plas, E., Hamilton, M. J., Miller, J. N., Koscik, T. R., Long, J. D., Cumming, S., et al. (2019). Brain structural features of myotonic dystrophy type 1 and their relationship with CTG Repeats. J. Neuromuscul. Dis. 6, 321–332. doi: 10.3233/JND-190397
Van Lieverloo, G. G. A., Wieske, L., Verhamme, C., Vrancken, A. F. J., van Doorn, P. A., Michalak, Z., et al. (2019). Serum neurofilament light chain in chronic inflammatory demyelinating polyneuropathy. J. Peripher. Nerv. Syst. 24, 187–194. doi: 10.1111/jns.12319
Whittaker, R. G., Blackwood, J. K., Alston, C. L., Blakely, E. L., Elson, J. L., McFarland, R., et al. (2009). Urine heteroplasmy is the best predictor of clinical outcome in the m.3243A>G mtDNA mutation. Neurology 72, 568–569.
Keywords: myopathy, serum neurofilament, nervous system involvement, mitochondriopathies, myotonic dystrophies (DM1 and DM2), facio scapulo humeral dystrophy
Citation: Saak A, Benkert P, Akgün K, Willemse E, Kuhle J, Ziemssen T, Jackson S and Schaefer J (2021) Serum Neurofilament Light Chain: A Marker of Nervous System Damage in Myopathies. Front. Neurosci. 15:791670. doi: 10.3389/fnins.2021.791670
Received: 08 October 2021; Accepted: 29 November 2021;
Published: 17 December 2021.
Edited by:
Helene Blasco, INSERM U1253 Imagerie et Cerveau (iBrain), FranceReviewed by:
Dario Ronchi, University of Milan, ItalyMukesh Gautam, Northwestern University, United States
Copyright © 2021 Saak, Benkert, Akgün, Willemse, Kuhle, Ziemssen, Jackson and Schaefer. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Annika Saak, QW5uaWthLlNhYWtAdWtkZC5kZQ==
†These authors share senior authorship