Pauline Brochet1
Barbara Ianni2
João P. S. Nunes1,2,3,4
Amanda F. Frade2,3,4
Priscila C. Teixeira2,3,4
Charles Mady5
Ludmila R. P. Ferreira6
Andreia Kuramoto2
Cristina W. Pissetti7
Bruno Saba8Darlan D. S. Cândido2,3,4
Fabrício Dias9
Marcelo Sampaio8José A. Marin-Neto9Abílio Fragata8
Ricardo C .F. Zaniratto2
Sergio Siqueira10Giselle D. L. Peixoto10
Vagner O. C. Rigaud2Paula Buck11
Rafael R. Almeida2,3,4
Hui Tzu Lin-Wang8
André Schmidt9Martino Martinelli10
Mario H. Hirata12
Eduardo Donadi9
Virmondes Rodrigues Junior7
Alexandre C. Pereira11
Jorge Kalil2,3,4
Lionel Spinelli1*
Edecio Cunha-Neto2,3,4*†
Christophe Chevillard1*†- 1Aix Marseille Univ, TAGC Theories and Approaches of Genomic Complexity, Institut MarMaRa, Marseille, France
- 2Laboratory of Immunology, Heart Institute Instituto do Coração(InCor), School of Medicine, University of São Paulo, São Paulo, Brazil
- 3Division of Clinical Immunology and Allergy, School of Medicine, University of São Paulo, São Paulo, Brazil
- 4Instituto Nacional de Ciência e Tecnologia, INCT, iii- Institute for Investigation in Immunology, São Paulo, Brazil
- 5Myocardiopathies and Aortic Diseases Unit, Heart Institute, Instituto do Coração (InCor), School of Medicine, University of São Paulo, São Paulo, Brazil
- 6RNA Systems Biology Laboratory (RSBL), Departamento de Morfologia, Instituto de Ciências Biológicas, Universidade Federal de Minas Gerais, Belo Horizonte, Minas Gerais, Brazil
- 7Laboratory of Immunology, Universidade Federal Do Triângulo Mineiro (UFTM), Uberaba, Brazil
- 8Laboratório de Investigação Molecular em Cardiologia, Instituto de Cardiologia Dante Pazzanese (IDPC), São Paulo, Brazil
- 9School of Medicine of Ribeirão Preto (FMRP), University of São Paulo, Ribeirão Preto, Brazil
- 10Pacemaker Clinic, Heart Institute Instituto do Coração (InCor), School of Medicine, University of São Paulo, São Paulo, Brazil
- 11Heart Institute Instituto do Coração (InCor), School of Medicine, University of São Paulo, São Paulo, São Paulo, Brazil
- 12Department of Clinical and Toxicological Analyses, Faculty of Pharmaceutical Sciences, University of São Paulo (USP), São Paulo, Brazil
Chagas disease is a parasitic disease from South America, affecting around 7 million people worldwide. Decades after the infection, 30% of people develop chronic forms, including Chronic Chagas Cardiomyopathy (CCC), for which no treatment exists. Two stages characterized this form: the moderate form, characterized by a heart ejection fraction (EF) ≥ 0.4, and the severe form, associated to an EF < 0.4. We propose two sets of DNA methylation biomarkers which can predict in blood CCC occurrence, and CCC stage. This analysis, based on machine learning algorithms, makes predictions with more than 95% accuracy in a test cohort. Beyond their predictive capacity, these CpGs are located near genes involved in the immune response, the nervous system, ion transport or ATP synthesis, pathways known to be deregulated in CCCs. Among these genes, some are also differentially expressed in heart tissues. Interestingly, the CpGs of interest are tagged to genes mainly involved in nervous and ionic processes. Given the close link between methylation and gene expression, these lists of CpGs promise to be not only good biomarkers, but also good indicators of key elements in the development of this pathology.
Introduction
Chagas disease is an endemic disease from South America, caused by a parasite, Trypanosoma cruzi, and affecting around 7 million people. With migration flow, this disease can now be found in non-endemic country, notably in North America (1) (n > 300,000), Europe (2) (n > 100,000), Japan (3) (n > 4,000) or Australia (4) (n > 1,000). After the infection, patients present an acute stage which is mostly asymptomatic (ASY). Then comes the chronic forms, where 70% of them remains asymptomatic, with no end organ damage (the so-called indeterminate stage). However, 30% develop Chagas disease Cardiomyopathy (CCC) (5). CCC had been divided in two stages based on heart ejection fraction: moderate CCC (EF ‗ 0.4) and severe CCC (EF < 0.4) (6–8). Some drugs are effective on T. cruzi, but does not cure the CCC, reducing the parasitemia, without having any effect on heart damage (9). The only way out for CCC patients is the placement of a pacemaker, or a heart transplant. The early diagnosis of Chagas disease is therefore essential.
During the acute stage, disease diagnosis is commonly made by microscopy, considering the limited sensitivity of the direct test (10). However, in chronic stage, parasitemia is very low, or even null. The Pan American Health Organization (PAHO) recommends using two serological tests (two techniques based on different antigens) in parallel and, in case of discordant results, to perform these tests again on a new sample (11). If the results remain unclear, a confirmation test should be achieved (12). For CCC especially, an ECG and/or an echocardiogram is made to confirm cardiac involvement (13). BNP and NT-proBNP, well-known markers be associated with cardiac dysfunction (14), have been associated to Chagas cardiomyopathy (15, 16), but are not specific to this pathology. Others markers, including miRNAs (17, 18), cytokines (19) or metalloproteinases (20) have been proposed as biomarkers for CCC, but no confirmation has been made in a test cohort at this time. The only diagnosis of CCC currently in place is a clinical diagnosis, which is difficult to access for the most remote populations.
A previous analysis (21) has highlighted differences of DNA methylation in blood of asymptomatic and CCC patients. Moreover, some differences have also been demonstrated between moderate and severe CCC. Blood DNA methylation has already been proposed as biomarker for several diseases (22–24). Here, we used machine learning methods on both asymptomatic and CCC blood DNA methylation data to predict Chagas disease, as well as Chagas disease stage.
Methods
Ethical considerations
The protocol was approved by the institutional review boards of the University of São Paulo School of Medicine and INSERM (French National Institute of Health and Medical Research). Written informed consent was obtained from all patients. All experimental methods comply with the Helsinki Declaration.
Blood DNA collection and DNA methylation analysis
Blood samples (5 to 15 ml of blood) from CCC patients were collected in EDTA tubes. Genomic DNA was isolated using standard salted methods and the methylation analysis was done using the same protocol as tissue DNAs.
Blood DNA methylation data
138 patients were selected randomly from our Chagas bank. It included 48 asymptomatic subjects, 46 moderate CCC patients and 44 severe CCC patients (Supplementary Table 1). The age and sex ratio were not significantly different between the 3 groups (age mean and ratio female/male for all phenotypes: asymptomatic: age: 57.63, ratio =1; moderate: age: 56.89, ratio=1.14; severe: age: 59.59, ratio=0.95). In a second time, these 138 samples were randomly distributed between the training (70%) and validation (30%) cohorts. This random distribution was done in such a way that the age and sex ratio was still not so different between the groups in the two sub-cohorts (Training cohort (age mean and ratio female/male for all phenotypes: asymptomatic: age: 62.45, ratio=1; moderate: age: 60.13, ratio=1.14; severe: age: 57.18, ratio=0.86), validation cohort (age and ratio female/male for all phenotypes: asymptomatic: age: 52.82, ratio=1; moderate: age: 53.67, ratio=1.13; severe: age: 62, ratio=1)). The methylation data are available under the reference: (GEO accession: GSE191082).
Biomarker identification for disease forms
Since data contains a lot of features (736,661), feature selection was performed in two steps, on the training group only. The scripts used for the following steps are available on Github (https://github.com/TAGC-ComplexDisease/biomarkersChagas). First, the delta beta (difference of beta means) was computed between the two phenotypes of interest. Only the CpGs having at least 10% methylation differences were retained. Then, a machine learning (ML) analysis was done in Python with Scikit-learn library. Four supervised ML methods were considered: decision tree, random forest, logistic regression and linear SVM (Support-Vector Machine, a linear classificator). For each method, recursive feature elimination (RFE) was performed, and the best model (best accuracy) with the minimal set of feature was selected using 10-time cross-validation. Finally, model parameters were optimized with a grid search to obtain the final prediction on the validation group.
Results
Symptomatic cardiac form prediction
After feature selection based on delta beta values, 86 CpGs were selected. Among all the tested models, linear SVM seems to have the better prediction on training dataset with the minimal number of features (Figure 2A). According to this analysis, the SVM was trained with 35 features (Supplementary Table 2). The model parameters optimization was performed using a grid search where the L2 penalty varies between 0.01 and 10. Finally, with a L2 penalty of 1, 42 of 44 patients phenotype of the validation dataset were correctly predicted (accuracy = 0.95), with a sensitivity of 0.96 and a specificity of 0.94 (area under the curve: 0.996) (Figure 1A). Those 35 features are mainly located in the body of genes (n=20), or in intergenic regions (n=11). Particularly, 3 CpGs are located in LHX6, and 3 in POU6F2. All those genes are involved in biological process associated to Chagas disease: nervous system (LHX6, POU6F2, MDGA1, DISC1, PCSK9), immune system (ZMIZ, HLA-DRB1), Wnt pathway (DISC1), ion transport (KCNK15, PCSK9), striated muscles (SMYD3) or ATP metabolic process (ATP5S).
Figure 1 Receiver Operating Characteristic (ROC) curves produced by (A) linear SVM to predict CCC on 44 patients and (B) random forest to predict CCC stage on 28 patients.
Chagas cardiomyopathy stage prediction
After feature selection based on delta beta, 108 CpGs were selected. Among all the tested models, random forest seems to have the better prediction on training dataset with the minimal number of features (Figure 2B). According to this analysis, the SVM was trained with 33 features (Supplementary Table 3). The model parameters optimization was performed using a grid search where the number of estimators varies between 50 and 200. Finally, 150 estimators, 27 of 28 patients phenotype (accuracy = 0.96) of the validation dataset were correctly predicted, with a sensitivity of 1 and a specificity of 0.93 (area under the curve: 1) (Figure 1B). Those 33 features are mainly located in 18 intergenic regions. Other CpGs are located in 15 genes, and more precisely in 6 gene body and 7 promoter regions. Here, genes are involved in various biological processes, from ion transport (KCNC1, MFI2), actin filament (PACSIN1), generation of neurons (TNN, PACSIN1) or MAPK cascade (DUSP22). 2 CpGs are in common with those used as biomarker between ASY and CCC: cg24000535 (LOC101928909) and cg21873524 (intergenic).
Figure 2 Evolution of the accuracy obtained with different machine learning models according to the number of top-explicative features chosen for (A) predict CCC or (B) predict CCC stage.
Discussion
Alterations of heart tissue DNA methylation profiles have been associated to the development of dilated cardiomyopathies (25, 26) and chagas diseases (21, 27). Recently, we have studied the DNA methylation in the blood of asymptomatic, moderate and severe CCC by hypothesizing that the blood data reflect the phenotype. We had found 12624 DMPs (Differentially Methylated Position) between asymptomatic and severe CCC blood samples and 6735 CpGs were found as DMPs between moderate and severe CCC.
In our study, based on machine learning approaches, we have identified 35 CCC-specific methylation markers. Those CpGs could distinguish controls (asymptomatic) to CCC from blood samples with 96% of sensitivity and 94% of specificity in independent validation sets. 3 of them (cg02872767, cg24540763, cg25134647) are also differentially methylated between asymptomatic and severe CCC in heart tissue (21). Similarly, 33 CpGs have been identified, allowing to predict the progression of this pathology (from moderate to severe CCC), with a sensitivity of 100% and a specificity of 93%. In conclusion, we identified two set of methylation markers potentially useful for Chagas disease diagnostic. The first one permit to discriminate patients with Chagas cardiomyopathy from asymptomatic patients, with 95% of precision. The second allows to predict Chagas cardiomyopathy severity stage, according to the heart ejection fraction rate, with 96% of precision.
Interestingly, most of these markers were not differentially in heart tissue of patients. The main message of this report is the finding that peripheral blood epigenetic marks are good markers of clinical form, implying that epigenetic events are closely related to CCC progression. These markers are highlighting the same biological processes that have been associated to the disease development such as ion transport, ATP metabolic process, immune system, Wnt system, nervous systems, striated muscles and actin filament. These findings are important as these are under-lighting biological pathways that will have to be targeted in drug design.
To propose large-scale reproducible biomarkers, a consensus Target Product Profile (TPP) has been developed for Chagas disease (28) stipulating that marker should be able to detect the effects of drug treatments, be detectable with limited resources and not vary according to the strain of the parasite. Given the high specificity of these assays, these methylation sites appear to be good candidates to decipher the pathogenic process, or to be used as blood biomarkers, and further studies will be necessary to potentially validate their possible use in the clinic, in accordance with the TPP consensus.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
The studies involving human participants were reviewed and approved by INSERM IRB University of Sao Paulo. The patients/participants provided their written informed consent to participate in this study.
Author contributions
Study design: PBr, JK, LS, ECN, CC. Phenotype characterization: BI, CM, SS, CWP, BS, FD, MS, JAMN, AF, GDLP, PBu, HTL-W, AS, MM, MHH, ED, ACP, VRJ Experimental analysis: PBr, AFF, JPSN, PCT, LRPF, AK, DDSC, RCFZ, VOCR, RRA. Statistical analysis: PBr, LS, ECN, CC. Manuscript preparation: PBr, LS, ECN, CC.
Funding
This work was supported by the Institut National de la Santé et de la Recherche Médicale (INSERM); the Aix-Marseille University (grant number: AMIDEX “International_2018” MITOMUTCHAGAS); the French Agency for Research (Agence Nationale de la Recherche-ANR (grant numbers: “Br-Fr-Chagas”, “landscardio”); the CNPq (Brazilian Council for Scientific and Technological Development); and the FAPESP (São Paulo State Research Funding Agency Brazil (grant numbers: 2013/50302-3, 2014/50890-5); the National Institutes of Health/USA (grant numbers: 2 P50 AI098461-02 and 2U19AI098461-06). This work was founded by the Inserm Cross-Cutting Project GOLD. This project has received funding from the Excellence Initiative of Aix-Marseille University - A*Midex a French “Investissements d’Avenir programme”- Institute MarMaRa AMX-19-IET-007. JN was a recipient of a MarMaRa fellowship. EC-N, JK, ALR and ECS are recipients of productivity awards by CNPq. The funders did not play any role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Acknowledgments
Center de Calcul Intensif d'Aix-Marseille is acknowledged for granting access to its high performance computing resources.
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/fimmu.2022.1020572/full#supplementary-material
References
1. Bern C, Messenger LA, Whitman JD, Maguire JH. Chagas Disease in the United States: a Public Health Approach. Clin Microbiol Rev (2019) 33. doi: 10.1128/CMR.00023-19
2. Antinori S, Galimberti L, Bianco R, Grande R, Galli M, Corbellino M. Chagas disease in Europe: A review for the internist in the globalized world. Eur J Intern Med (2017) 43:6–15. doi: 10.1016/j.ejim.2017.05.001
3. Imai K, Maeda T, Sayama Y, Osa M, Mikita K, Kurane I, et al. Chronic Chagas disease with advanced cardiac complications in Japan: Case report and literature review. Parasitol Int (2015) 64:240–2. doi: 10.1016/j.parint.2015.02.005
4. Jackson Y, Pinto A, Pett S. Chagas disease in Australia and New Zealand: risks and needs for public health interventions. Trop Med Int Health (2014) 19:212–8. doi: 10.1111/tmi.12235
5. Pérez-Molina JA, Molina I. Chagas disease. Lancet (2018) 391:82–94. doi: 10.1016/s0140-6736(17)31612-4
6. Frade AF, Teixeira PC, Ianni BM, Pissetti CW, Saba B, Wang LH, et al. Polymorphism in the alpha cardiac muscle actin 1 gene is associated to susceptibility to chronic inflammatory cardiomyopathy. PloS One (2013) 8:e83446. doi: 10.1371/journal.pone.0083446
7. Nogueira LG, Santos RH, Ianni BM, Fiorelli AI, Mairena EC, Benvenuti LA, et al. Myocardial chemokine expression and intensity of myocarditis in Chagas cardiomyopathy are controlled by polymorphisms in CXCL9 and CXCL10. PloS Negl Trop Dis (2012) 6:e1867. doi: 10.1371/journal.pntd.0001867
8. van Veldhuisen DJ, Linssen GC, Jaarsma T, van Gilst WH, Hoes AW, Tijssen JG, et al. B-type natriuretic peptide and prognosis in heart failure patients with preserved and reduced ejection fraction. J Am Coll Cardiol (2013) 61:1498–506. doi: 10.1016/j.jacc.2012.12.044
9. Morillo CA, Marin-Neto JA, Avezum A, Sosa-Estani S, Rassi A Jr., Rosas F, et al. Randomized Trial of Benznidazole for Chronic Chagas' Cardiomyopathy. N Engl J Med (2015) 373:1295–306. doi: 10.1056/NEJMoa1507574
10. Balouz V, Aguero F, Buscaglia CA. Chagas Disease Diagnostic Applications: Present Knowledge and Future Steps. Adv Parasitol (2017) 97:1–45. doi: 10.1016/bs.apar.2016.10.001
11. Organization PAH. Guidelines for the diagnosis and treatment of Chagas disease. Washington, D.C: Pan American Health Organization (PAHO (2019).
12. Lapa JS, Saraiva RM, Hasslocher-Moreno AM, Georg I, Souza AS, Xavier SS, et al. Dealing with initial inconclusive serological results for chronic Chagas disease in clinical practice. Eur J Clin Microbiol Infect Dis (2012) 31:965–74. doi: 10.1007/s10096-011-1393-9
13. Benck L, Kransdorf E, Patel J. Diagnosis and Management of Chagas Cardiomyopathy in the United States. Curr Cardiol Rep (2018) 20:131. doi: 10.1007/s11886-018-1077-5
14. Chow SL, Maisel AS, Anand I, Bozkurt B, de Boer RA, Felker GM, et al. Role of Biomarkers for the Prevention, Assessment, and Management of Heart Failure: A Scientific Statement From the American Heart Association. Circulation (2017) 135:e1054–e91. doi: 10.1161/CIR.0000000000000490
15. Echeverria LE, Rojas LZ, Gomez-Ochoa SA, Rueda-Ochoa OL, Sosa-Vesga CD, Muka T, et al. Cardiovascular biomarkers as predictors of adverse outcomes in chronic Chagas cardiomyopathy. PloS One (2021) 16:e0258622. doi: 10.1371/journal.pone.0258622
16. Brito BOF, Pinto-Filho MM, Cardoso CS, Di Lorenzo Oliveira C, Ferreira AM, de Oliveira LC, et al. Association between typical electrocardiographic abnormalities and NT-proBNP elevation in a large cohort of patients with Chagas disease from endemic area. J Electrocardiol (2018) 51:1039–43. doi: 10.1016/j.jelectrocard.2018.08.031
17. Gomez-Ochoa SA, Bautista-Nino PK, Rojas LZ, Hunziker L, Muka T, Echeverria LE. Circulating MicroRNAs and myocardial involvement severity in chronic Chagas cardiomyopathy. Front Cell Infect Microbiol (2022) 12:922189. doi: 10.3389/fcimb.2022.922189
18. Linhares-Lacerda L, Granato A, Gomes-Neto JF, Conde L, Freire-de-Lima L, de Freitas EO, et al. Circulating Plasma MicroRNA-208a as Potential Biomarker of Chronic Indeterminate Phase of Chagas Disease. Front Microbiol (2018) 9:269. doi: 10.3389/fmicb.2018.00269
19. De Alba-Alvarado M, Salazar-Schettino PM, Jimenez-Alvarez L, Cabrera-Bravo M, Garcia-Sancho C, Zenteno E, et al. Th-17 cytokines are associated with severity of Trypanosoma cruzi chronic infection in pediatric patients from endemic areas of Mexico. Acta Trop (2018) 178:134–41. doi: 10.1016/j.actatropica.2017.11.009
20. Medeiros NI, Gomes JAS, Fiuza JA, Sousa GR, Almeida EF, Novaes RO, et al. MMP-2 and MMP-9 plasma levels are potential biomarkers for indeterminate and cardiac clinical forms progression in chronic Chagas disease. Sci Rep (2019) 9:14170. doi: 10.1038/s41598-019-50791-z
21. Brochet P, Ianni B, Laugier L, Frade AF, Nunes JPS, Teixeira P, et al. Epigenetic regulation of transcription factor binding motifs promotes Th1 response in Chagas disease Cardiomyopathy. Front Immunol (2022). doi: 10.3389/fimmu.2022.958200
22. Fransquet PD, Lacaze P, Saffery R, Phung J, Parker E, Shah R, et al. Blood DNA methylation signatures to detect dementia prior to overt clinical symptoms. Alzheimers Dement (Amst) (2020) 12:e12056. doi: 10.1002/dad2.12056
23. Lange CP, Campan M, Hinoue T, Schmitz RF, van der Meulen-de Jong AE, Slingerland H, et al. Genome-scale discovery of DNA-methylation biomarkers for blood-based detection of colorectal cancer. PloS One (2012) 7:e50266. doi: 10.1371/journal.pone.0050266
24. Shu C, Justice AC, Zhang X, Marconi VC, Hancock DB, Johnson EO, et al. DNA methylation biomarker selected by an ensemble machine learning approach predicts mortality risk in an HIV-positive veteran population. Epigenetics (2021) 16:741–53. doi: 10.1080/15592294.2020.1824097
25. Jo BS, Koh IU, Bae JB, Yu HY, Jeon ES, Lee HY, et al. Methylome analysis reveals alterations in DNA methylation in the regulatory regions of left ventricle development genes in human dilated cardiomyopathy. Genomics (2016) 108:84–92. doi: 10.1016/j.ygeno.2016.07.001
26. Haas J, Frese KS, Park YJ, Keller A, Vogel B, Lindroth AM, et al. Alterations in cardiac DNA methylation in human dilated cardiomyopathy. EMBO Mol Med (2013) 5:413–29. doi: 10.1002/emmm.201201553
27. Laugier L, Frade AF, Ferreira FM, Baron MA, Teixeira PC, Cabantous S, et al. Whole-Genome Cardiac DNA Methylation Fingerprint and Gene Expression Analysis Provide New Insights in the Pathogenesis of Chronic Chagas Disease Cardiomyopathy. Clin Infect Dis (2017) 65:1103–11. doi: 10.1093/cid/cix506
Keywords: chagas disease, cardiomyopathy, blood, biomarkers, methylation
Citation: Brochet P, Ianni B, Nunes JPS, Frade AF, Teixeira PC, Mady C, Ferreira LRP, Kuramoto A, Pissetti CW, Saba B, Cândido DDS, Dias F, Sampaio M, Marin-Neto JA, Fragata A, Zaniratto RCF, Siqueira S, Peixoto GDL, Rigaud VOC, Buck P, Almeida RR, Lin-Wang HT, Schmidt A, Martinelli M, Hirata MH, Donadi E, Rodrigues Junior V, Pereira AC, Kalil J, Spinelli L, Cunha-Neto E and Chevillard C (2022) Blood DNA methylation marks discriminate Chagas cardiomyopathy disease clinical forms. Front. Immunol. 13:1020572. doi: 10.3389/fimmu.2022.1020572
Received: 16 August 2022; Accepted: 15 September 2022;
Published: 29 September 2022.
Edited by:
Ulrike Kemmerling, University of Chile, ChileReviewed by:
Juan Diego Maya, University of Chile, ChileJulio A. Urbina, Instituto Venezolano de Investigaciones Cientificas, Venezuela
Israel Molina, Vall d’Hebron University Hospital, Spain
Copyright © 2022 Brochet, Ianni, Nunes, Frade, Teixeira, Mady, Ferreira, Kuramoto, Pissetti, Saba, Cândido, Dias, Sampaio, Marin-Neto, Fragata, Zaniratto, Siqueira, Peixoto, Rigaud, Buck, Almeida, Lin-Wang, Schmidt, Martinelli, Hirata, Donadi, Rodrigues Junior, Pereira, Kalil, Spinelli, Cunha-Neto and Chevillard. 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: Christophe Chevillard, christophe.chevillard@univ-amu.fr; Edecio Cunha-Neto, edecunha@gmail.com; Lionel Spinelli, lionel.spinelli@univ-amu.fr
†These authors have contributed equally to this work