Skip to main content

ORIGINAL RESEARCH article

Front. Neurol., 28 November 2022
Sec. Movement Disorders

Data-driven clustering of combined Functional Motor Disorders based on the Italian registry

  • 1Section of Neurosciences, Department “G.F. Ingrassia”, University of Catania, Catania, Italy
  • 2Oasi Research Institute—IRCCS, Troina, Italy
  • 3Neurology Unit, Movement Disorders Division, Department of Neurosciences, Biomedicine and Movement Sciences, University of Verona, Verona, Italy
  • 4Center for Neurodegenerative Diseases (CEMAND), Department of Medicine, Surgery and Dentistry—Scuola Medica Salernitana, University of Salerno, Baronissi, Italy
  • 5Neurology Unit, Department of Clinical and Experimental Medicine, University of Pisa, Pisa, Italy
  • 6Department of Clinical and Experimental Sciences, University of Brescia, Brescia, Italy
  • 7FERB Onlus, Ospedale S. Isidoro, Bergamo, Italy
  • 8Parkinson and Movement Disorders Unit, Fondazione IRCCS Istituto Neurologico Carlo Besta, Milan, Italy
  • 9Department of Medical Area, Neurology Unit, ASST Pavia, Pavia, Italy
  • 10Botulinum Toxin Center, Neurology Unit A.O.U. Mater Domini, Catanzaro, Italy
  • 11Aldo Ravelli Research Center for Neurotechnology and Experimental Brain Therapeutics, Department of Health Sciences, University of Milan, Milan, Italy
  • 12IRCCS Neuromed, Pozzilli, Italy
  • 13Department of Medicine and Aging Sciences, University G. d'Annunzio, Chieti, Italy
  • 14Department of Medicine and Aging Sciences, University G. d'Annunzio, Pescara, Italy
  • 15IRCCS Humanitas Research Hospital, Milan, Italy
  • 16Department of Human Neurosciences, Università La Sapienza, Rome, Italy
  • 17Department of Advanced Medical and Surgery Sciences, University of Campania—Luigi Vanvitelli, Naples, Italy
  • 18Department of Neuroscience—Rita Levi Montalcini, University of Turin, Turin, Italy
  • 19Department of Biomedical and Neuromotor Sciences, University of Bologna, Bologna, Italy
  • 20IRCCS, Institute of Neurological Sciences of Bologna, Bologna, Italy
  • 21Movement Disorder Unit, Fondazione Policlinico Universitario A. Gemelli IRCCS, Rome, Italy
  • 22Neurosciences Research Centre, Molecular and Clinical Sciences Neurosciences Research Centre, Molecular and Clinical Sciences Research Institute, St George's University of London, London, United Kingdom
  • 23Department of Experimental and Clinical Medicine, University of Messina, Messina, Italy
  • 24Clinical Neurophysiology Unit, Cardarelli Hospital, Naples, Italy
  • 25Department of Neurosciences, Reproductive and Odontostomatological Sciences, University of Naples—Federico II, Naples, Italy
  • 26IRCCS Mondino Foundation, Pavia, Italy
  • 27Department of Brain and Behavioral Sciences, University of Pavia, Pavia, Italy
  • 28Clinical Neurology Unit, Department of Medical, Surgical and Health Services, University of Trieste, Trieste, Italy
  • 29University and Institute of Research and Medical Care San Raffaele, Rome, Italy
  • 30Ospedale degli Infermi, Biella, Italy
  • 31Psichiatria, Ospedale San Maurizio, Bolzano, Italy
  • 32Department of Medical Sciences and Public Health, University of Cagliari, Cagliari, Italy

Introduction: Functional Motor Disorders (FMDs) represent nosological entities with no clear phenotypic characterization, especially in patients with multiple (combined FMDs) motor manifestations. A data-driven approach using cluster analysis of clinical data has been proposed as an analytic method to obtain non-hierarchical unbiased classifications. The study aimed to identify clinical subtypes of combined FMDs using a data-driven approach to overcome possible limits related to “a priori” classifications and clinical overlapping.

Methods: Data were obtained by the Italian Registry of Functional Motor Disorders. Patients identified with multiple or “combined” FMDs by standardized clinical assessments were selected to be analyzed. Non-hierarchical cluster analysis was performed based on FMDs phenomenology. Multivariate analysis was then performed after adjustment for principal confounding variables.

Results: From a study population of n = 410 subjects with FMDs, we selected n = 188 subjects [women: 133 (70.7%); age: 47.9 ± 14.4 years; disease duration: 6.4 ± 7.7 years] presenting combined FMDs to be analyzed. Based on motor phenotype, two independent clusters were identified: Cluster C1 (n = 82; 43.6%) and Cluster C2 (n = 106; 56.4%). Cluster C1 was characterized by functional tremor plus parkinsonism as the main clinical phenotype. Cluster C2 mainly included subjects with functional weakness. Cluster C1 included older subjects suffering from anxiety who were more treated with botulinum toxin and antiepileptics. Cluster C2 included younger subjects referring to different associated symptoms, such as pain, headache, and visual disturbances, who were more treated with antidepressants.

Conclusion: Using a data-driven approach of clinical data from the Italian registry, we differentiated clinical subtypes among combined FMDs to be validated by prospective studies.

Introduction

Functional Motor Disorders (FMDs) still represent nosological entities with no clear phenotypic characterization (1, 2). To be defined, they require evidence of clinical inconsistency and incongruence in the context of functional weakness or hypo-/hyperkinesia. Once defined, FMD characterization is still principally based on main clinical features (1, 2). However, in most cases, patients may present a combined phenomenology of two or more FMDs, leading to a clinical classification that is usually based on the predominant symptom judged by the referring clinicians and then potentially biased (2).

We have recently observed that individuals with single, isolated FMD without comorbid neurological disorders have largely overlapping demographic and clinical non-motor characteristics (2). On the other hand, the characterization of multiple or “combined” FMDs remains still difficult, considering the variability related to the clinical judgment, as well as the temporal inconsistency of the disorder.

Data-driven approach using clustering analysis of clinical data has been proposed as an analytic method to obtain non-hierarchical, unbiased classifications to be tested in clinical practice. It has been applied to detect possible disease phenotypes in different fields of movement disorders, including Parkinson's disease (3) and essential tremor (4).

In this study, we aimed to identify possible clinical subtypes of combined FMDs using a data-driven approach, to overcome possible limits related to “a priori” classifications and clinical overlapping.

Materials and methods

Study population and data collection

Data were obtained by the Italian Registry of Functional Motor Disorders (IRFMD), managed by the University of Verona and the Italian Academy for the Study of Parkinson's Disease and Other Movement Disorders (Accademia LIMPE-DISMOV RADAC project) and Fondazione LIMPE. Details about subjects' recruitment and data collection have been described elsewhere (1, 2).

In brief, consecutive outpatients with FMDs were recruited from 25 tertiary movement disorders centers (11 in northern, five in central, six in southern Italy, and three in Sardinia/Sicily) between 1 September 2018 and 31 August 2019. Inclusion criteria were: age ≥10 years; the presence of one or more FMDs; a clinically definite diagnosis of FMDs based on Gupta and Lang diagnostic criteria (1, 2). Each patient with FMDs underwent a detailed clinical assessment, including screening for possible associated comorbidities. Patients with comorbid neurological diseases were included in the study (1, 5). FMD phenotypes were defined based on their specific phenomenological features in (a) tremor, (b) weakness, (c) dystonia; (d) myoclonus-like jerks; (e) facial motor disorders; (f) parkinsonism; and (g) gait disorders (1, 2). For functional weakness, we applied clues for the functional cause suggested by the referred Gupta and Lang diagnostic criteria, including the presence of distractibility maneuvers and the demonstration of positive signs (1). Since we focused on the symptom of functional weakness without selecting cases based on any etiological assumptions, no specific DSM criteria were applied for the diagnosis. Exclusion criteria were the presence of cognitive or physical impairment that precluded signing the informed consent form for participation in the study, as well as data acquisition from the structured interview (1).

Patients identified with one (isolated FMDs) or multiple FMDs (combined FMDs; e.g., dystonia plus tremor) underwent standardized clinical assessments. Patients were assessed at each center in a single session by a neurologist specialized in movement disorders. Demographic and clinical data were systematically collected, including associated symptoms, comorbidities, and practiced treatment (1, 2).

Approval was obtained by the Institutional Ethics Committee of the Coordinator Center (University of Verona, Azienda Ospedaliera Universitaria Integrata Verona, Prog. 1757CESC) and confirmed by the Committees of each participating center. All patients (or their guardians) gave their written consent to participate.

Statistical analysis

Data are presented as mean ± standard deviation for continuous variables and frequency (percentage) for categorical variables.

Non-hierarchical cluster analysis using the k-means method was performed based on the main clinical phenomenology of FMDs: tremor, weakness, dystonia, myoclonus-like jerks, facial movement disorder, parkinsonism, and gait disorder (n = 7 variables). The Calinski-Harabasz pseudo-F stopping-rule criterion was used to evaluate the optimal number of clusters to be analyzed. We performed a cluster analysis by selecting the subgroup of patients with combined FMDs from the study population. The χ2 goodness of fit test was applied to determine if multicategorical variables distribution in the selected sample followed the expected distribution of the entire study population.

Differences in scalar variables among groups were tested using the independent-samples t-test, after testing for normality. Categorical variables were tested using the χ2-test. Unconditional logistic regression analysis was performed and for each study variable, we calculated Odds Ratios (OR), 95% Confidence Interval (CI), and p-value (two-tailed). Multivariate analysis was also performed after adjustment for principal confounding variables (age, sex, and disease duration). The significance level was set at p = 0.1.

Results

We analyzed data from a total sample size of n = 410 enrolled subjects [women: 291 (71%); age 46.6 ± 15.8 years; disease duration: 5.6 ± 6.8 years)] presenting isolated and combined FMDs. FMDs phenomenology distribution in the study population was: tremor [n = 167 (40.7%)], weakness [n = 180 (43.9%)], dystonia [n = 119 (29%)], myoclonus [n = 53 (12.9%)], facial movement disorder [n = 47 (11.5%)], parkinsonism [n = 24 (5.8%)], and gait disorder [n = 109 (26.6%)]. Detailed descriptive analysis of demographic and clinical data of the study population has been already reported elsewhere (1).

A study subgroup of n = 188 (45.8%) subjects [women: 133 (70.7%); age: 47.9 ± 14.4 years; disease duration: 6.4 ± 7.7 years] with combined FMD was selected to be analyzed with data clustering. FMDs phenomenology distribution in this study subgroup was: tremor [n = 109 (58%)], weakness [n = 106 (56.4%)], dystonia [n= 78 (41.5%)], myoclonus [n = 35 (18.6%)], facial movement disorder [n = 35 (18.6%)], parkinsonism [n = 22 (11.7%)], and gait disorder [n = 92 (48.9%)]. There were no differences between the selected sample and the entire population in terms of age (t = 0.96; p = 0.337), sex (χ2 = 0.003; p = 0.954), and disease duration (t = 1.28; p = 0.201). The frequency distribution of overall FMDs differed from those in the entire study population (goodness of fit χ2 = 76.6; p < 0.001).

Two independent groups of patients based on clustering parameters were identified (for n = 2 to six clusters: n = 2 pseudo F = 45.34; n = 3 pseudo F = 40.16; n = 4 pseudo F = 36.68; n = 5 pseudo F = 31.27; n = 6 pseudo F = 34.66): Cluster C1 (n = 82; 43.6%) and Cluster C2 (n = 106; 56.4%). The distribution of FMDs among groups is shown in Table 1. Statistically significant differences in evaluated clinical variables between the two clusters are shown in Table 2, including multivariate analysis adjusted for age, sex, and disease duration.

TABLE 1
www.frontiersin.org

Table 1. Distribution of FMDs among identified clusters.

TABLE 2
www.frontiersin.org

Table 2. Differences in evaluated clinical variables between the identified clusters.

As for the total population, the first detected Cluster C1 principally included subjects with functional tremor plus parkinsonism as the main FMDs clinical phenotype. The tremor was present in n = 59 (71.9%) patients: in n = 49 (83%) patients, it was distributed in the upper limb, in n = 18 (30.5%) patients in the lower limb, and in n = 15 (25.4%) patients in body parts other than limbs. Parkinsonism was instead present in n = 18 (21.9%) patients: in n = 15 (83.3%) patients, it was distributed in the upper limb, in n = 8 (44.4%) patients in the lower limb, and in n = 10 (55.6%) patients in body parts other than limbs.

Functional weakness was instead more represented in Cluster C2, being present in all n = 106 patients. In n = 42 (39.6%) patients, it was distributed in the upper limb, in n = 86 (81.1%) patients in the lower limb, and in n = 12 (11.3%) patients in body parts other than limbs. In Cluster C2, parkinsonism was present in n = 4 (3.8%) patients: in n = 3 (75%) patients, it was distributed in the upper limb, in n = 3 (75%) patients in the lower limb, and in n = 3 (75%) patients in body parts other than limbs. In all patients affected by parkinsonism, functional weakness principally affecting limbs was observed. Other hyperkinetic functional movement disorders were equally distributed in the two clusters.

Cluster C1 included older subjects with associated anxiety. They presented a diagnosis of parkinsonism as neurological comorbidity, even if the association was lacking in the multivariate analysis, and they were more treated with BTX and antiepileptics. Cluster C2 included instead younger subjects referring to more associated symptoms, such as pain, headache, visual disturbances, and feelings of depersonalization/derealization, even though this last factor lacked a significant association in the multivariate analysis. They were treated with antidepressants with respect to Cluster C1.

Discussion

In this study, we aimed to identify possible clinical subtypes of combined FMDs by clustering registry-based data. We used patients' main clinical features as clustering variables, then we looked at differences in terms of comorbidities associated with FMDs (5), precipitating and associated factors, as well as response to therapies among the identified groups.

When cluster analysis was applied to main clinical features with no “a priori” distinction on combined FMDs, two separate clusters were identified with different characteristics. The first cluster (Cluster C1) principally included subjects with functional movement disorders, specifically functional tremor plus parkinsonism. It also mainly included older subjects experiencing anxiety. They presented a diagnosis of parkinsonism as neurological comorbidity and they were treated with BTX and antiepileptics. The second cluster (Cluster C2) was characterized by weakness as the main clinical phenotype. The group principally included younger subjects referring to different associated symptoms, including headache, debilitating pain, visual symptoms, and feelings of depersonalization/derealization. Subjects in this group were more treated with antidepressants.

The high frequency of functional parkinsonism and tremor in Cluster C1 may reflect previous results showing that functional tremor represents one of the most frequent FMDs observable in elderly subjects (6). No specific differences in terms of gender distribution were obtained between the two groups, even though Cluster C1 included overall more men than Cluster C2. Previous reports highlighted a higher prevalence of conversion symptoms, including functional weakness and sensory symptoms in women with respect to men, probably explained by a higher rate of consultation for the referred condition in this group (79).

Concerning differences in associated conditions between groups, highlighted characteristics for Cluster C2 seemed in agreement with the previously reported association between functional weakness and sensory symptoms with respect to other phenotypes (1, 7, 8). By a pathophysiological state point, it has been argued that such results may be explained by the hypoactivity of subcortical circuits controlling sensorimotor function and driving both voluntary motor behavior and sensory loss, differently from what is hypothesized for functional symptoms characterized by excessive motor output (1, 10).

We also found a high prevalence of feelings of depersonalization/derealization in Cluster C2, in which functional weakness is more represented, with respect to Cluster C1. This may reflect a previous hypothesis postulating that functional weakness could be a form of “hemi-depersonalization” in which the relationship with the functional weakness could be explained in terms of immobility and/or increased attention paid to a specific body part (11).

Coexisting neurological comorbidities in FMDs, including parkinsonism, were already described in our study population (5). Patients with comorbid FMDs were usually older and more frequently had tremors (5), characteristics observed in Cluster C1. Systematic review already demonstrated that tremor is the most common FMD observed in Parkinson's disease, usually occurring in the most affected side (5, 12). It has been suggested that the underlying pathological process in Parkinson's disease predisposes patients to abnormal attention-to-movement production (12). The same abnormal attention toward movement has been demonstrated in FMD, a notion supported by the clinical evidence of functional movement normalization by shifting attention away from movement (12). Based on this assumption, it could be argued that Parkinson's disease pathology may produce vulnerability to develop functional symptoms. In the case of functional symptoms in the context of Parkinson's disease, they would be mainly motor symptoms linked phenomenologically to the physical deficit caused by the disease (12).

Our study has some limitations that should be considered. Selection bias cannot be excluded since the enrollment was carried out in tertiary, referral centers and collected data may be also biased by recall bias. Furthermore, data clustering was applied just to the combined FMDs, thus affecting the results generalizability. However, it should be underlined that the exclusion of the “isolated” FMDs from the data clustering was due to methodological limits. Indeed, the stopping-rule method could not be applied to identify an optimal number of clusters in the subpopulation of isolated FMDs because the optimal number of identifiable clusters overlaps with the number of the identified FMDs phenotypes used for the motor classification, as previously defined.

There are also strengths of the present study. This is the first study attempting an unbiased, data-driven classification of FMDs when combined, to overcome possible evaluation bias due to the subjective “a priori” classification. It also includes a large multicenter sample of FMDs patients which were systematically evaluated by a standardized protocol.

In conclusion, using a data-driven approach of clinical data from the Italian registry, we were able to detect clinical subtypes of combined FMDs to be validated by prospective studies. In future studies, it will be crucial to expand the ambition of diagnostic quality to postmortem studies.

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 authors.

Ethics statement

The studies involving human participants were reviewed and approved by University of Verona, Azienda Ospedaliera Universitaria Integrata Verona, Prog. 1757CESC. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.

Author contributions

GM: study design, data collection, data and statistical analysis, writing of the first draft, editing, and final draft. CG, AL, MZa, MT, and REr: study design, data collection, review, critique, editing, and final draft. AN: study design, data collection, data and statistical analysis, writing of the first draft, review, critique, editing, and final draft. EM, PBa, RC, SM, APil, APa, LR, REl, CD, CA, FB, PAB, BD, OG, NM, EO, LB, AA, GF, RD, MZi, GC-B, MP, FM, ME, APis, PM, FS, MCM, IDV, LT, FDB, TE, and GD: study design, data collection, review, and critique. All authors contributed to the article and approved the submitted version.

Acknowledgments

List of study co-investigators has been provided as supplementary material (Appendix - Table A1).

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.

References

1. Tinazzi M, Morgante F, Marcuzzo E, Erro R, Barone P, Ceravolo R, et al. Clinical correlates of functional motor disorders: an Italian multicenter study. Mov Disord Clin Pract. (2020) 7:920–9. doi: 10.1002/mdc3.13077

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Tinazzi M, Geroin C, Marcuzzo E, Cuoco S, Ceravolo R, Mazzucchi S, et al. Functional motor phenotypes: to lump or to split? J Neurol. (2021) 268:4737–43. doi: 10.1007/s00415-021-10583-w

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Erro R, Vitale C, Amboni M, Picillo M, Moccia M, Longo K, et al. The heterogeneity of early Parkinson's disease: a cluster analysis on newly diagnosed untreated patients. PLoS ONE. (2013) 8:e70244. doi: 10.1371/journal.pone.0070244

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Mostile G, Terranova R, Rascunà C, Terravecchia C, Cicero CE, Giuliano L, et al. Clinical-Instrumental patterns of neurodegeneration in Essential Tremor: a data-driven approach. Parkinsonism Relat Disord. (2021) 87:124–9. doi: 10.1016/j.parkreldis.2021.05.011

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Tinazzi M, Geroin C, Erro R, Marcuzzo E, Cuoco S, Ceravolo R, et al. Functional motor disorders associated with other neurological diseases: beyond the boundaries of “organic” neurology. Eur J Neurol. (2021) 28:1752–8. doi: 10.1111/ene.14674

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Batla A, Stamelou M, Edwards MJ, Pareés I, Saifee TA, Fox Z, et al. Functional movement disorders are not uncommon in the elderly. Mov Disord. (2013) 28:540–3. doi: 10.1002/mds.25350

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Stone J, Carson A, Duncan R, Coleman R, Roberts R, Warlow C, et al. Symptoms 'unexplained by organic disease' in 1144 new neurology out-patients: how often does the diagnosis change at follow-up? Brain. (2009) 132:2878–88. doi: 10.1093/brain/awp220

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Stone J, Carson A, Duncan R, Coleman R, Roberts R, Warlow C, et al. Who is referred to neurology clinics?–the diagnoses made in 3,781 new patients. Clin Neurol Neurosurg. (2010) 112:747–51. doi: 10.1016/j.clineuro.2010.05.011

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Baizabal-Carvallo JF, Jankovic J. Gender differences in functional movement disorders. Mov Disord Clin Pract. (2019) 7:182–7. doi: 10.1002/mdc3.12864

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Baizabal-Carvallo JF, Hallett MJ. Jankovic Pathogenesis and pathophysiology of functional (psychogenic) movement disorders. Neurobiol Dis. (2019) 127:32–44. doi: 10.1016/j.nbd.2019.02.013

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Stone J, Warlow CM. Sharpe functional weakness: clues to mechanism from the nature of onset. J Neurol Neurosurg Psychiatry. (2012) 83:67–9. doi: 10.1136/jnnp-2011-300125

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Pareés I, Saifee TA, Kojovic M, Kassavetis P, Rubio-Agusti IA, Sadnicka KP, et al. Functional (psychogenic) symptoms in Parkinson's disease. Mov Disord. (2013) 28:1622–7. doi: 10.1002/mds.25544

PubMed Abstract | CrossRef Full Text | Google Scholar

Appendix

TABLE A1
www.frontiersin.org

Table A1. Co-investigators.

Keywords: cluster analysis, clinical phenotypes, Functional Motor Disorders, data-driven phenotyping, functional neurological disorder

Citation: Mostile G, Geroin C, Erro R, Luca A, Marcuzzo E, Barone P, Ceravolo R, Mazzucchi S, Pilotto A, Padovani A, Romito LM, Eleopra R, Dallocchio C, Arbasino C, Bono F, Bruno PA, Demartini B, Gambini O, Modugno N, Olivola E, Bonanni L, Albanese A, Ferrazzano G, De Micco R, Zibetti M, Calandra-Buonaura G, Petracca M, Morgante F, Esposito M, Pisani A, Manganotti P, Stocchi F, Coletti Moja M, Di Vico IA, Tesolin L, De Bertoldi F, Ercoli T, Defazio G, Zappia M, Nicoletti A and Tinazzi M (2022) Data-driven clustering of combined Functional Motor Disorders based on the Italian registry. Front. Neurol. 13:987593. doi: 10.3389/fneur.2022.987593

Received: 06 July 2022; Accepted: 17 October 2022;
Published: 28 November 2022.

Edited by:

Jesús De Pedro-Cuesta, Carlos III Health Institute (ISCIII), Spain

Reviewed by:

Mohammad Sami Al-Wardat, Jordan University of Science and Technology, Jordan
Robert Fekete, New York Medical College, United States

Copyright © 2022 Mostile, Geroin, Erro, Luca, Marcuzzo, Barone, Ceravolo, Mazzucchi, Pilotto, Padovani, Romito, Eleopra, Dallocchio, Arbasino, Bono, Bruno, Demartini, Gambini, Modugno, Olivola, Bonanni, Albanese, Ferrazzano, De Micco, Zibetti, Calandra-Buonaura, Petracca, Morgante, Esposito, Pisani, Manganotti, Stocchi, Coletti Moja, Di Vico, Tesolin, De Bertoldi, Ercoli, Defazio, Zappia, Nicoletti and Tinazzi. 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: Michele Tinazzi, michele.tinazzi@univr.it; Christian Geroin, christian.geroin@univr.it

These authors have contributed equally to this work

Disclaimer: 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.