Evaluation of dental morphometrics during the orthodontic treatment
© Kašparováet al.; licensee BioMed Central Ltd. 2014
Received: 6 February 2014
Accepted: 23 May 2014
Published: 3 June 2014
Diagnostic orthodontic and prosthetic procedures commence with an initial examination, during which a number of individual findings on occlusion or malocclusion are clarified. Nowadays we try to replace standard plaster casts by scanned objects and digital models.
Geometrically calibrated images aid in the comparison of several different steps of the treatment and show the variation of selected features belonging to individual biomedical objects. The methods used are based on geometric morphometrics, making a new approach to the evaluation of the variability of features. The study presents two different methods of measurement and shows their accuracy and reliability.
The experimental part of the present paper is devoted to the analysis of the dental arch objects of 24 patients before and after the treatment using the distances between the canines and premolars as the features important for diagnostic purposes. Our work proved the advantage of measuring digitalized orthodontic models over manual measuring of plaster casts, with statistically significant results and accuracy sufficient for dental practice.
A new method of computer imaging and measurements of a dental stone cast provides information with the precision required for orthodontic treatment. The results obtained point to the reduction in the variance of the distances between the premolars and canines during the treatment, with a regression coefficient R C=0.7 and confidence intervals close enough for dental practice. The ratio of these distances pointed to the nearly constant value of this measure close to 0.84 for the given set of 24 individuals.
KeywordsOrthodontic modelling Geometric morphometrics Digital models Dental arch features Digital signal processing Regression analysis Computational intelligence
It is becoming increasingly evident in orthodontics and dentofacial orthopedics that the timing of the onset of treatment may be as critical as the selection of a specific treatment protocol. By beginning at the patient’s optimal maturational stage, the most favorable response with the lowest potential morbidity can be anticipated. The issue of optimal timing for dentofacial orthopedics is closely linked to the identification of growth periods that can contribute significantly to the correction of skeletal imbalances in the individual patient. Cephalometric investigations on longitudinal samples have identified a pubertal spurt in mandibular growth that possesses wide individual variations in onset, duration, and rate [1–4].
The multidisciplinary dental care of patients is not a simple matter. New prosthodontic methods, including implant insertion, can be instituted using therapy by the surgeon, orthodontist and prosthodontist . Dental casts play an important role in the diagnosis and treatment planning [6–8] in prosthodontics and orthodontics. Digitalization [9, 10] is an important part of medicine using paperless patient information systems as virtual charts, digital photographs, and digital dental cross.
Although plaster casts are the gold standard in treatment planning, not only in orthodontics, the replacement of plaster dental casts by their digital models [11–14] can be advantageous, especially in saving space in storage areas, the efficiency of having patient records accessible through a computer, the possibility of sharing the models with other specialists needed during the therapy, and the possibility of accurate measurements and the use of diagnosis setups.
The analysis of a study cast, consisting of the three-dimensional assessment of the maxillary and mandibular dental arches and the maximal intercuspal relationship, is one of the basic tools of diagnosis and treatment planning. The arch form, dimensions, and variations obtained from orthodontic and prosthodontic treatment have been the subject of study for many years now . Some authors have tried to identify a geometric curve that would facilitate the accurate definition of arch forms . The first studies showing the importance of 3D dental arch analysis were conducted in the previous century. Bonwill and Hawley described the alignment of the upper anterior teeth as a circumference arch; whilst MacConaill and Scher maintained that the dental arch resembled a catenary curve [16, 17]. Izard, in trying to relate the dimension of the dental arch to the facial dimensions, found that the arch form could be accurately represented by an elliptical curve , as can be also found in our group of plaster casts.
Both methods, manual measurement and measurement using 3D models of dental casts [19–23], have advantages and disadvantages which can influence the result of the measurement and therefore also the treatment planning, treatment evaluation or face reconstruction .
The aim of this study is to validate the accuracy and use of digitized dental plaster casts in dental practice [35, 36], using a group of 48 plaster casts. It forms a contribution to the study of dental arch parameters using segmentation, classification and registration [37, 38] of orthodontic data.
We created a set of 24 patient’s plaster casts, randomly selected from the patients who underwent an orthodontic treatment at the Department of Orthodontics, 2nd Medical Faculty, Charles University, Prague, Czech Republic. Informed consents were obtained from all patients or their legal representatives. No ethical approval was required for this study.
All patients suffered from Class II malocclusion and had already finished orthodontic treatment by a fixed orthodontic appliance and extraction of the first premolars. There were no orthodontic appliances present. Two study plaster casts of the upper jaw (pre- and post treatment) were made of each patient, which means there were 48 plaster casts. All of the plaster casts were of high quality with no fractures or other damage.
A Roland LPX-250 scanner was used to scan all of the plaster casts. This is non-contact laser scanner which picks up several points to create coordinates and converts them into 3D data [39, 40] with a lateral resolution of 200 μ m. Multiple scanning can be done in up to six planes. The laser automatically works in conjunction with the Dr. Picza program, in which scanned data is entered and from which rough models are exported to other programs. The subsequent adjustment of the models involves the projection of surface re-scans onto the individual rotational scans, adjustment of all scans and re-scans (e.g. identification and removal of abnormal surfaces, or deletion of excessively long surfaces), and the registration and merging of individual re-scans with the rotational scan. The imperfect surface is further adjusted using a smoothing function; the small apertures that develop during adjustment of the scans are filled in according to the original model. The resulting models may then be further analysed using existing software, i.e., sections may be made according to previously defined planes, and the shape of any structure may be measured.
The Pixform software package was used for measuring the digitalized plaster casts. All models were displayed in the same colour and zoom and rotation functions were used to find the ideal point of measurement.
For the purposes of this study, the measurements of the dental arch width were made twice on each of the casts and twice on each of the digital models of dental casts. All measurements were made by a single examiner in the same conditions and were repeated in 72 hours.
In the first part of the study, the measurements of the width of the dental arch between the canines (3-3) and the second premolars (5-5) of the 24 sets of plaster casts were made manually. A sliding scale was used for this purpose. While measuring the plaster casts, the casts were held in the hand and manual rotation was used to find the tips in the cusps of the measured teeth.
In the second part of the study, the same measurement as in the first part was made on digitalized plaster casts of the study group. The casts were scanned on a rotating plate from a position perpendicular to the occlusal plane. The raw scan data were processed using Pixform reverse engineering software (Roland DG). This procedure included cleaning, merging of multiple scans, hole-filling, decimating, smoothing, and global remeshing.
Distances [mm] between corresponding canines (3-3) and premolars (5-5) measured twice both for digital models and plaster casts for 24 patients (i) before and (ii) after the treatment
Digital model measurement
Plaster cast measurement
for the t and χ2 distributions, respectively, and selected confidence level p=100 (1−α)%.
The comparison of the results obtained from the computational digital models and from the plaster casts  was studied through histograms and the distribution of the measured values. The differences in the distances obtained from the digital models and those from the plaster casts were tested by the pair t-test. The same test was used for the differences in the distances before and after the treatment. The precisions of the distance measurement of the digital models and the plaster casts were tested using the standard F-test  for the null hypothesis that the two normal sets have the same variance.
All the tests were performed at the selected significance level α=0.05 separately for the distances measured between the canines (3-3) and those between the premolars (5-5). Selected statistical tests were performed in the Matlab Statistical Toolbox [49, 50].
The plaster casts were evaluated using classical methods and geometrical morphometrics. Changes in the transverse direction of growth were determined from the distance between the tips of the permanent canines and the premolars, which we denote by 3-3 and 5-5.
Statistical characteristics including confidence intervals (CI) of digital model (DM) and plaster cast (PC) before and after the treatment
DM measurement differences
PC measurement differences
DM 3-3 before the treatment
PC 3-3 before the treatment
DM 5-5 before the treatment
PC 5-5 before the treatment
DM 3-3 after the treatment
PC 3-3 after the treatment
DM 5-5 after the treatment
PC 5-5 after the treatment
Digital measurements provide slightly better results, with a smaller standard deviation. The distances between the canines (3-3) and the premolars (5-5) before and after the treatment obtained by digital and manual measurements are presented in Table 2 together with their 95% confidence intervals (CI) as well. These results correspond to registering digital models obtained after and before the treatment .
The t-test used for comparison of the distances measured (i) on digital models and (ii) plaster casts proved that data in both vectors are random samples from normal distributions with equal means and variances at the 5% significance levels with the 95% confidence interval of the difference between the population means for distances 3-3 in the range 〈−0.55;1.78〉, and for distances 5-5 in the range 〈−1.25;0.92〉. The F-test applied to the two different kinds of measurements proved that the null hypothesis of the variances’ being equal can be accepted at the significance level of 5% with value p=0.58.
With present day knowledge, the analysis of study casts as a whole is often considered to have limited diagnostic value . A new method of computer imaging and measurements on a dental stone cast is a ubiquitous tool in dentistry and aids in the recording of the precise information [55, 56] required for treatment, does not destroy hard dental tissues, and provides stable results. The present paper monitored this process, step-by-step. In both methods we have to count on the possibility of subjective mistakes, which means both methods are influenced by the skills of the measuring person. Because of this fact, we can not reliably say which of the two methods mentioned above is better.
The statistical analysis of the geometric measurements obtained from the digital models points to a decrease in the variances of the distances between the premolars (5-5) related to the distances between the canines (3-3) after the treatment, and to the regression coefficient being close to 0.7 for the given set of individuals, with a high correlation and close confidence intervals. The mean value of ratio of distances between the canines (3-3) and the premolars (5-5) after the treatment for 24 individuals is 0.84 with the regression coefficient R C=0.0008 close to zero. This result corresponds with typical dental arch parameters.
This contribution compared the accuracy of two different methods used for the measurement of distances in orthodontics and used for treatment planning and its evaluation. The resulting digital model was then used for the analysis of selected dental morphometrics.
As a result of our measurements, we can say that measurements of digital models have the same accuracy as measurements of plaster casts, and they moreover allow more convenient way of data processing in clinical practice.
This conclusion is in accordance with modern trends to move to digital formats for orthodontic models, as suggested by the American Board of Orthodontics [57–60]. Physical casts, either plaster casts or models printed by 3D digital printers, can be replaced by digital models in many cases.
A specific following study devoted to distances between the canines (3-3) and the premolars (5-5) before and after the treatment for 24 individuals performed by digital models pointed to the nearly constant value of their ratio close to 0.84. The mean value of the relative difference between values evaluated from the digital model and plaster cast was 3%.
Our further research will be devoted to more precise 3-D modelling and to an algorithmic approach to the evaluation of the orthodontic parameters used for the treatment and follow up care of patients. The digital detection of reference points on three dimensional models and the volume registration  of models in different stages of the treatment will be studied as well.
Real orthodontic data were kindly provided by the Faculty Hospital Motol.
- Lewis A, Roche AF, Wagner B: Pubertal spurts in cranial base and mandible: comparisons within individuals. Angle Orthod 1985, 55: 17–30.Google Scholar
- Hagg U, Pancherz H, Taranger J: Pubertal growth and orthodontic treatment. In Craniofacial Growth During Adolescence. Craniofacial Growth Series, Volume 20. Edited by: Carlson DS, Ribbens KA. Ann Arbor, MI: Center for Human Growth and Development, University of Michigan; 1987:87–115.Google Scholar
- Opheij DG, Opdebeeck H, van Steenberghe D, Quirynen M: Age as compromising factor for implant insertion. Periodontol 2000 2003, 33: 172–184. 10.1046/j.0906-6713.2003.03314.xView ArticleGoogle Scholar
- Thilander B, Odman J, Jemt T: Single implants in the upper incisor region and their relationship to the adjacent teeth. an 8-year follow-up study. Clin Oral Implants Res 1999, 10: 346–355. 10.1034/j.1600-0501.1999.100502.xView ArticleGoogle Scholar
- Dostalova T, Racek J, Lozekova E, Rerchova M: Composite veneers, crowns, and inlay bridges after orthodontic therapy - a three-year prospective study. Gen Dent 2003, 51: 129–132.Google Scholar
- Conti MF, Filho MV, Vedovello SAS, Valdrighi HC, Kuramae M: Longitudinal evaluation of dental arches individualized by the WALA ridge method. J Orthod (Dent Press) 2011,16(2):65–74. 10.1590/S2176-94512011000200009View ArticleGoogle Scholar
- Rheudea B, Sadowsky P. L, Ferrierac A, Jacobson A: An evaluation of the use of digital study models in orthodontic diagnosis and treatment planning. Angle Orthod 2005,75(3):300–304.Google Scholar
- Ghislanzoni LTH, Lineberger M, Cevidanes LHS, Mapelli A, Sforza C, McNamara JA: Evaluation of tip and torque on virtual study models: a validation study. Prog Orthod 2013,14(19):1–6.Google Scholar
- Peluso MJ, Josell SD, Levine SW, Lorei BJ: Digital models: an introduction. Elsevier: Semin Orthod 2004,10(3):226–238. 10.1053/j.sodo.2004.05.007Google Scholar
- Santoro M, Galkin S, Teredesai M, Nicolay OF, Cangialosi TJ: Comparison of measurements made on digital and plaster models. Am J Orthod Dentofacial Orthop 2003,124(1):101–105. 10.1016/S0889-5406(03)00152-5View ArticleGoogle Scholar
- Costalos PA, Sarraf K, Cangialosi TJ, Efstratiadis S: Evaluation of the accuracy of digital model analysis for the American board of rthodontics objective grading system for dental casts. Am J Orthod Dentofacial Orthop 2005,128(5):624–628. 10.1016/j.ajodo.2004.08.017View ArticleGoogle Scholar
- Grauer D, Proffit WR: Accuracy in tooth positioning with a fully customized lingual orthodontic appliance. Am J Orthod Dentofacial Orthop 2011,140(2):433–444.View ArticleGoogle Scholar
- Hildebrand JC, Palomo JM, Palomo L, Sivik M, Hans M: Evaluation of a software program for applying the American board of orthodontics objective grading system to digital casts. Am J Orthod Dentofacial Orthop 2008,133(2):283–289. 10.1016/j.ajodo.2006.03.035View ArticleGoogle Scholar
- Leifert MF, Leifert MM, Efstratiadis SS, Cangialosi TJ: Comparison of space analysis evaluations with digital models and plaster dental casts. Am J Orthod Dentofacial Orthop 2009,136(1):16–1164. 10.1016/j.ajodo.2009.03.002View ArticleGoogle Scholar
- Murad SM, Al-Mulla AA: Accuracy of measurements made on digital and study models (A comparative study). Malays Dental J (MDJ) 2010,7(1):71–82.Google Scholar
- Tweed CH: The frankfort mandibular incisor angle in orthodontic diagnosis, treatment planning and prognosis. Angle Orthod 1954, 15: 1212–1269.Google Scholar
- MacConaill MA, Scher E: The ideal form of the human dental arcade, with some prosthetic application. Dent Rec 1949, 69: 285–302.Google Scholar
- Izard G: New method for the determination of the normal arch by the function of the face. Int J Orthod 1927,13(7):582–595.Google Scholar
- Adaskevicius R, Vasiliauskas A: Evaluation of dental arch form using 3D dental cast scanning technology. Electron Electrical Eng: Med Technol 2009,93(5):99–102.Google Scholar
- Adaskevicius R, Vasiliauskas A: Three-dimensional determination of dental occlusion and facial structures using soft tissue cephalometric analysis. Electron Electrical Eng: Syst Eng Comput Technol 2012,121(5):93–96.Google Scholar
- Grauer D: Three-dimensional applications in orthodontics. PhD thesis. University of North Carolina, Department of Oral Biology at the School of Dentistry; 2010Google Scholar
- Al-Khatib AR, Rajion ZA, Masudi SM, Hassan R, Townsen GC: Validity and reliability of tooth size and dental arch measurements a stereo photogrammetric study. Aust Orthod J 2012,28(1):22–29.Google Scholar
- El-Zanaty HM, El-Beialy AR, El-Ezz AMA, Attia KH, El-Bialy AR, Mostafa YA: Three-dimensional dental measurements: an alternative to plaster models. Am J Orthod Dentofacial Orthop 2010,137(2):259–265. 10.1016/j.ajodo.2008.04.030View ArticleGoogle Scholar
- Rosati R, DeMenezes M, Rossetti A, Sforza C, Ferrario V. F: Digital dental cast placement in 3-dimensional, full-face reconstruction: a technical evaluation. Am J Orthod Dentofacial Orthop 2010,138(1):84–88. 10.1016/j.ajodo.2009.10.035View ArticleGoogle Scholar
- Yamamoto K, Hayashi S, Nishikawa H, Nakamura S, Mikami T: Measurements of dental cast profile and three-dimensional tooth movement during orthodontic treatment. IEEE Trans on Biomed Eng 1991,38(4):360–365. 10.1109/10.133232View ArticleGoogle Scholar
- Kondo T, Ong SH, Foong KWC: Tooth segmentation of dental study models using range images. IEEE Trans on Med Imaging 2004,23(3):350–362. 10.1109/TMI.2004.824235View ArticleGoogle Scholar
- Chapuis J, Schramm A, Pappas I, Hallermann W, Schwenzer-Zimmerer K, Langlotz F, Caversaccio M: A new system for computer-aided preoperative planning and intraoperative navigation during corrective jaw surgery. IEEE Trans on Inf Technol Biomed 2007,11(3):274–287.View ArticleGoogle Scholar
- Yaqi M, Zhongke L: Computer aided orthodontics treatment by virtual segmentation and adjustment. In 2010 International Conference on Image Analysis and Signal Processing (IASP). IEEE; 2010:336–339.View ArticleGoogle Scholar
- Chang YB, Xia JJ, Gateno J, Xiong Z, Zhou X, Wong STC: An automatic and robust algorithm of reestablishment of digital dental occlusion. IEEE Trans on Med Imaging 2010,29(9):1652–1663.View ArticleGoogle Scholar
- Liang W, Yang L, Wang S, Wang B: Three-dimensional finite element analysis of maxillary first molar orthodontics. In 2010 3rd International Conference on Biomedical Engineering and Informatics (BMEI 2010). IEEE; 2010:1287–1291.View ArticleGoogle Scholar
- Zelditch ML, Swiderski DL, Sheets HD, Fink WL: Geometric Morphometrics for Biologists. Elsevier, Academic Press, London; 2012.Google Scholar
- Slice DE: Geometric morphometrics. Annu Rev Anthropol 2007, 36: 261–281. 10.1146/annurev.anthro.34.081804.120613View ArticleGoogle Scholar
- Dostalova T, Racek J, Tauferova E, Seydlova M, Smutny V, Bartonova M: Composite veneers, crowns, and inlay bridges after orthodontic therapy - a three-year prospective study. Methods Inf Med 2006, 45: 191–194.Google Scholar
- Stevens DR, Flores-Mir C, Nebbe B, Raboud DW, Heo G, Major PW: Validity, reliability, and reproducibility of plaster vs digital study models: comparison of peer assessment rating and Bolton analysis and their constituent measurements. Am J Orthod Dentofacial Orthop 2006,129(6):794–803. 10.1016/j.ajodo.2004.08.023View ArticleGoogle Scholar
- Quimby ML, Vig KW, Rshid RG, Firestone AR: The accuracy and reliability of measurements made on computer-based digital models. Angle Orthod 2004,74(3):298–303.Google Scholar
- Bootvong K, Liu Z, McGrath C, Hagg U, Wong RW, Bendeus M, Yeung S: Virtual model analysis as an alternative approach to plaster model analysis: reliability and validity. Eur J Orthod 2010,32(5):589–595. 10.1093/ejo/cjp159View ArticleGoogle Scholar
- Grafova L, Kasparova M, Kakawand S, Prochazka A, Dostalova T: Study of edge detection task in dental panoramic x-ray images. Dentomaxillofacial Radiol 2012,42(7):0391/1–20120391/12. http://www.ncbi.nlm.nih.gov/pubmed/23640989 Google Scholar
- Grauer D, Cevidanes LH, Tyndall D, Styner MA, Flood PM, Proffit WR: Registration of Orthodontic Digital Models. In Effective and efficient orthodontic tooth movement. Craniofacial Growth Series, Volume Monograph 48. Edited by: McNamara JA, Hatch N, Kapila SD. USA: Needham Press, Needham, MA 02492; 2011:377–392.Google Scholar
- Dugelay J. L, Baskurt A, Daoudi M: 3D Object Processing. Chichester, UK: John Wily & Sons; 2008.View ArticleGoogle Scholar
- Dvorak P, Prochazka A, Kasparova M, Dostalova T: Orthodontic Data Acquisition and Visualization. In Technical Computing. Humusoft, ICT Prague; 2011:1–6.Google Scholar
- Jacquet W, Nyssen E, Ibel G, Vannet B. V: On the augmented reproducibility in measurements on 3D orthodontic digital dental models and the definition of feature points. Aust Orthod J 2013,29(1):28–33.Google Scholar
- Freedman D. A: Statistical Models. New York: Cambridge University Press; 2005.Google Scholar
- Gardiner WP: Statistics for the Biosciences. New York: Prentice Hall; 1997.Google Scholar
- Goulden CH: Methods of Statistical Analysis. New York: Wiley; 1956.Google Scholar
- Armitage P, Berry G, Matthews JNS: Statistical Methods in Medical Research. Oxford: Blackwell Science; 2002.View ArticleGoogle Scholar
- Bury K: Statistical Distributions in Engineering. Cambridge, UK: Cambridge University Press; 1999.View ArticleGoogle Scholar
- Fleming PS, Marinho V, Johal A: Orthodontic measurements on digital study models compared with plaster models: a systematic review. Orthod Craniofac Res 2011,14(1):1–16. 10.1111/j.1601-6343.2010.01503.xView ArticleGoogle Scholar
- Brandt S: Data Analysis: Statistical and Computational Methods for Scientists and Engineers. USA: Springer-Verlag New York Inc; 1999.View ArticleGoogle Scholar
- Martinez WL, Martinez AR: Computational Statistics Handbook with MATLAB. Florida, USA: Chapman & Hall/CRC; 2008.Google Scholar
- Statistics Toolbox User’s Guide. 3 Apple Hill Drive, Natick, MA: The Mathworks. Inc; 2014.Google Scholar
- Creed B, Kau CH, English JD, Xia JJ, Lee RP: A comparison of the accuracy of linear measurements obtained from cone beam computerized tomography images and digital models. Elsevier: Semin Orthod 2011,17(1):49–56. 10.1053/j.sodo.2010.08.010Google Scholar
- Luu NS, Nikolcheva LG, Retrouveyc JM, Flores-Mird C, El-Bialye T, Careyf JP, Majorg PW: Linear measurements using virtual study models: a systematic review. Angle Orthod 2012,82(6):1098–1106. 10.2319/110311-681.1View ArticleGoogle Scholar
- Mayers M, Firestone AR, Rashid R, Vigd KWL: Comparison of peer assessment rating (PAR) index scores of plaster and computer-based digital models. Am J Orthod Dentofacial Orthop 2005,128(4):431–434. 10.1016/j.ajodo.2004.04.035View ArticleGoogle Scholar
- Graber TM, Vanarsdall RL, Vig KWL: Orthodontics: Current Principles and Techniques. Philadelphia PA 19103–2899; USA: Elsevier; 2012.Google Scholar
- Keating AP, Knox J, Bibb R, Zhurov AI: A comparison of plaster, digital and reconstructed study model accuracy. J Orthod 2008,35(3):191–201. 10.1179/146531207225022626View ArticleGoogle Scholar
- Kasparova M, Grafova L, Dvorak P, Dostalova T, Prochazka A, Eliasova H, Prusa J, Kakawand S: Possibility of reconstruction of dental plaster cast from 3d digital study models. BioMed Eng OnLine 2013,12(49):1–11.Google Scholar
- ABO: The American Board of Orthodontics (ABO) Digital Model Requirements. 2014.Google Scholar
- Kau CH, Littlefield J, Rainy N, Nguyen JT, Creed B: Evaluation of CBCT digital models and traditional models using the Little’s index. Angle Orthod 2010, 80: 435–439. 10.2319/083109-491.1View ArticleGoogle Scholar
- Cuperus AM, Harms MC, Rangel FA, Bronkhorst EM, Schols JG, Breuning KH: Dental models made with an intraoral scanner: a validation study. Am J Orthod Dentofacial Orthop 2012, 142: 308–313. 10.1016/j.ajodo.2012.03.031View ArticleGoogle Scholar
- Tarazona B, Llamas JM, Cibrian R, Gandia JL, Paredes V: A comparison between dental measurements taken from CBCT models and those taken from a digital method. Eur J Orthod 2013, 35: 1–6. 10.1093/ejo/cjr005View ArticleGoogle Scholar
- Ramalingam S, Taguchi Y: A theory of minimal 3D point to 3D plane registration and its generalization. Int J Comput Vis 2013,102(1–3):73–90. 10.1007/s11263-012-0576-xMathSciNetView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.