Quantitative Anatomical StudiesView this Special Issue
Research Article | Open Access
Roqueline A. G. M. F. Aversi-Ferreira, Hisao Nishijo, Tales Alexandre Aversi-Ferreira, "Reexamination of Statistical Methods for Comparative Anatomy: Examples of Its Application and Comparisons with Other Parametric and Nonparametric Statistics", BioMed Research International, vol. 2015, Article ID 902534, 10 pages, 2015. https://doi.org/10.1155/2015/902534
Reexamination of Statistical Methods for Comparative Anatomy: Examples of Its Application and Comparisons with Other Parametric and Nonparametric Statistics
Various statistical methods have been published for comparative anatomy. However, few studies compared parametric and nonparametric statistical methods. Moreover, some previous studies using statistical method for comparative anatomy (SMCA) proposed the formula for comparison of groups of anatomical structures (multiple structures) among different species. The present paper described the usage of SMCA and compared the results by SMCA with those by parametric test (t-test) and nonparametric analyses (cladistics) of anatomical data. In conclusion, the SMCA can offer a more exact and precise way to compare single and multiple anatomical structures across different species, which requires analyses of nominal features in comparative anatomy.
Some biological sciences are dependent on subjective decision of scientists due to lack of appropriate numerical methods to determine observations. The description of biological structures, for instance, is sometimes exhaustive and depends on scientists’ subjective observations. It is well known that the ancient researchers mistook analyses of the anatomical structures including their numbers, for example, the number of the cranial nerves, probably because of the scarce conditions during the works at that time . Nevertheless, most structures in the body analyzed in the past were reported to contain the same number of structures nowadays.
In fact, gross anatomy requires descriptions of qualitative variables including innervation, vascularization, origin and insertion of muscles, and arteries and nerve branches. It also requires quantitative analyses of mass, area, volume, size, and dimensional measures of such biological structures, which can be characterized by parametric statistical methods. However, parametric statistical methods can be hardly applied to description of qualitative variables in a scope of gross anatomy . This is one of the reasons why the different authors with great experiences differently described structures of the same species without considering anatomical variations (for an illustrative example of confusion regarding the name and description of the tibial artery, see ).
Here, we propose the main methodology to characterize qualitative data in gross anatomy, which enables us to describe and compare objectively anatomical structures across different species. Specifically, there have been no adequate quantitative methods to compare discrete-nominal variables of anatomical structures. It is desirable to quantitatively assess discrete variables . Indeed, a nonparametric statistical method has been proposed for comparative anatomy  and actually used in some works [5, 6]. Although these studies compared their nonparametric statistical method with one of nonparametric methods (chi-square), the authors did not compare their nonparametric method with parametric ones, nor did they specify possible problems of interpretation of the data. Moreover, for statistical method for comparative anatomy (SMCA), the formula must be modified so that it accepts data of multiple structures across different species.
For objective descriptions of the anatomical structures, some authors used the chi-square comparison to analyze nominal variables [7, 8] by converting frequencies of anatomical characteristics to percentages. Thus, they calculated standard deviations of the data in percentages. However, standard deviations of nonparametric data derived from a kind of the discrete categorical variables usually generate a statistical error type I . Furthermore, the basis of the chi-square statistic is causality among the data, a hypothesis that is not consistent with a theory of evolution assuming the concept of the common ancestral . Anatomy of body structures such as innervation, origin, and/or insertion of muscles of the arm, for instance, seems to be not random in animals that evolved from a common ancestral animal, since the ancestral animal provided basic structures and could generate derivative features in descendant animals (for a detailed review, see ). Nevertheless, the chi-square statistic is an important tool among multivariate analyses of discrete variables that are considered to be independent in quantitative psychology . Another nonparametric method (cladistics) was applied to comparative anatomy to analyze primitive and derivative features in evolution [11, 12] but has not yet been applied to descriptive anatomy.
In this paper we compared the nonparametrical method for SMCA [5, 6] with another parametric method (-test) analyzing means of the samples, using the previously published data . The SMCA including Comparative Anatomy Index for groups of structures (GCAI) that enables comparison of multiple structures  (see Section 3.2 for details) was also compared with another nonparametric (cladistics) method.
This work reanalyzed the previously published data in comparative anatomy statistic (SMCA)  to verify this method in detail  and applied the new formula (GCAI) to compare groups of structures among different species. Furthermore, using the published data [5, 6] as examples, we show the steps to calculate the SMCA and compare between the SMCA and a parametric test (-test) and also between the SMCA and other nonparametric (cladistics) methods .
3.1. Methods to Compare Samples from Same Species
The first step of this statistical method for comparative anatomy (SMCA) is to analyze the frequency based on the anatomical concept of normality and variation. “A normal structure” means that it is observed in greater than 50% of cases within the same species; therefore, the variation can be observed in less than 50% of cases :where is the total number of analyzed structures of the samples, is the number of structures with variation, and is the number of normal structures . The subscript () indicates specific species such as humans, Cebus, and baboon, while the subscript () indicates specific structures (flexor pollicis longus, pronator quadratus, etc.), and the subscript () indicates parameters of the specific structures. For muscle studies, the parameters should include at least the following 4 parameters: innervation, origin, insertion, and vascularization of muscles. For example, in case of the flexor pollicis longus muscle in Cebus , the data analyses in this step should be performed in terms of the following 4 parameters: the innervation (), origin (), insertion (), and vascularization (). Furthermore, number of muscles () and shape () could be added for more detailed analyses. In addition, further detailed parameters (subscript ()) could be added.
The relative frequency () of normal structures in each parameter against the total number of structures is defined as follows:
When structures are pair, will be the number of individuals multiplied by 2. It is also possible to calculate in separate pieces of bodies, as well. Although any values can be used as , smaller number of will result in lower statistical power. Normal structure in each parameter means in practical terms. However, in mathematical ones with normality concept, can vary as follows: .
In the same species, is usually greater than 0.5. However, in comparison among different species, normality is different among the species; for instance, in the comparison of the dorsoepitrochlearis muscle among primates and modern humans, this muscle is rarely observed in humans and approximate is 0.05 , while, in nonhuman primates, the dorsoepitrochlearis muscle is a normal feature, and the is 1.00.
On the other hand, the palmaris longus could be defect in humans  and its prevalence is around 90% ; therefore, might be 90% of total individuals. Thus, to calculate innervation, vascularization, origin, or insertion of the palmaris longus, only 90% receives attention and data from the remaining 10% are discarded. Such case is common in comparative studies, where, usually, only data in specific species are studied. Furthermore, some muscles have more than one origin or insertion, as in the triceps brachii with 3 heads, and ultimately this muscle has 4 heads of origin in modern humans [13, 15, 16]. In this case, just 2 types of the origin are observed: type 1 with 3 heads that is the normal feature and type 2 with 4 heads that is a variation.
For a more detailed analysis, it is required to calculate by including such parameters in muscle, nerves, bones, arteries, and so forth. For instance, in muscle studies, the parameters have to be chosen according to the goal of researches; parameters for muscle studies should include, at least, (1) innervation (), (2) origin (), (3) insertion (), and (4) vascularization (). Furthermore, (5) number of muscles () and (6) shape () could be added for more detailed analyses. It is noted that small number of parameters means that the studied structure is less characterized. For instance, in case of contrahentes muscles, the number of the muscles must be analyzed because they show variation within the same species and in different species of primates as well .
The calculation of is the first step in the SMCA analysis, and single could be compared among different species. However, what we are seeking is comparison of multiple among different species. The next step is to specify pondered values for coefficients (i.e., weighted coefficients) () that are multiplied by . The coefficients must be determined based on anatomical perspective; a parameter for a specific feature with small is not important when we assess similarity of a given structure. That is, since the small is ascribed to greater number of variations, small must accompany small weighted coefficient, while large (i.e., small number of variations) must accompany greater weighted coefficients.
We gave the weighted coefficient 3 to innervation in case of muscles. When the muscles are formed during the development of animals, a specific nerve terminates on a specific muscle . Thus, variation in nerve innervation of muscles is small, and a variation of innervation is highly sensitive to differences among different individuals within the same species as well as among different species. Among the 4 parameters noted above (innervation , origin , insertion , and vascularization ), origin and insertion usually show similar variations. Thus, the both weighted coefficients should receive the same weight coefficient 2 ( for origin and for insertion). Finally, the parameter with greater variation, vascularization (), received the weighted coefficient 1 (). Indeed, vascularization can be different between the same muscles in bilateral sides within the same individuals .
Zero cannot be accepted as weighted coefficient (). Therefore, must be greater than zero; that is, . To make the calculation easier and to keep clear parameters, the best choice is to use only integer values; that is, . Based on the above inference, an important rule here is clear; when the weighted coefficients are defined, the values should depend on different degree of variation (highest value to the weighted coefficient for the parameter with the lowest degree of variation or the same value to weighted coefficients for the parameters with same degree of variation). The values also should be discrete since it is difficult to find proportional values that represent exact difference among nominal variables. Thus, the best way is to choose integer values according to variations of studied features.
After designation of pondered values for weighted coefficients and calculation of , the next step for SMCA is to calculate the Pondered Average of Frequencies , according to the following formula:where is the relative frequency and is the weighted coefficient attached to a given parameter. For example, in muscle 1 of species 1, is relative frequency of innervation, and weighted coefficient is 3; is relative frequency of the muscle origin, and is 2; is relative frequency of muscle insertion, and is 2; and is relative frequency of vascularization, and is 1 .
In practical terms, must be greater than 0.5 and less than or equal to 1; that is, . In fact, could be 1 if every has maximal value 1, and if every is minimum (), the will be 0.5, as well. In mathematical terms, again, regardless of concept of normality, can vary within the range of , since could be zero or less than 0.5 in case of the analyses among different species. The value of can be used to assess quantitative difference among studied structures in that equal values indicate high similarity and large difference in the values between two species indicates dissimilarities or less similarity.
3.2. Methods to Compare Different Structures of the Same Species and among Different Species
To compare structures among different species or different structures in the same species, has to be calculated in each structure in each species and the must be estimated in comparison with the data of reference species (control species). For instance, the corachobrachialis muscle has one or two heads of origin depending on species of primates . In case of this muscle, (relative number of heads) could be different according to the number of heads in the control species. Thus, before calculating , it is important to make sure that the must be consistently calculated in comparison with control species (see below).
For example, the maximum number of types of origin is 2 in the corachobrachialis ; type 1 has one origin and type 2 has two origins . could take different values according to the number of heads in the reference species (control species) . For example, for noncontrol species to be studied , the of type 1 (number of origin is 1) will be 1 in reference to the species with 1 head, and of type 1 will be 0.5 in reference to the species with 2 heads. In case of the muscle that has 1 to 3 heads of origins across different species, the value should be divided by maximum number (i.e., 3) of heads (), since should not be greater than 1. Thus, when control species with 3 heads of origin is reference, in species with 3 heads is 1.000, in species with 2 heads is 0.667 (2/3), and in species with 1 head will be 0.333 (1/3).
It is also important that the values of should be obtained firstly in the control species. If the control species have normally two heads of origin in the corachobrachialis () and if 100% of individuals in this species have two heads of origin, in this species will be 1. In the case wherein 90% of individuals in this species have two heads, will be 0.9. In other noncontrol species in which the normal is one head of origin of the corachobrachialis, if 100% of individuals of the studied sample have one head, the will be 0.5, and if 90% of the samples have one head, will be 0.45. These values are used to estimate , which will be applied to the CAI analyses (see below).
Although any species can be used as control species, the species studied in the first time or the species with much known data should be chosen as control species. To compare any single structure (e.g., muscle) between two different species , the data in any noncontrol species can be compared one by one with those in the control species using the Comparative Anatomy Index (CAI) defined by the following formula:The represents an absolute difference of weighted averages () of a single structure between the control () and other noncontrol () species. To compare one structure with one parameter between the control and noncontrol species, the formula can be modified as follows:
It is noted that the ranges from 0 to 1; that is, . This is because the maximum value of is 1 and the minimum is 0. Note that this equation permits only comparison of just one structure between the 2 species. However, the SMCA analysis of the muscles in the forearm  reported necessity to compare multiple muscles among different species, for example, to compare groups of the deep flexor muscles in the forearm among different species, because these muscles work together for a common function. Comparisons of them as a group would indicate similarities in relation to functions, phylogeny, and taxonomy. Thus, the authors  suggested the GCAI to compare a group of the muscles among species, one by one based on the sum of the , as follows:and is the number of studied structures (e.g., muscles); indeed, because the same number of structures is mostly studied in each species.
The GCAI, which represents difference in based on multiple muscle structures between the control () and other noncontrol () species, is defined by the following formula:or
Based on the above inferences, using SMCA, the values close to 0.000 suggest high similarity of the structures between the species, and the value 1.000 indicates that those are completely different structures. Thus, the GCAI is the absolute difference in mean weighted averages of for multiple muscles between the two species and is defined in Table 1.
|Note. : proportion of normal structured organs, where represents individual species, represents individual structures (i.e., muscles), and represent individual parameters of the muscles.|
4.1. Examples of Application and Calculation of the SMCA
We reanalyzed the previous data of the anatomical structures in comparative anatomy by application of the SMCA. The statistic for comparative anatomy was applied to the muscle extensor pollicis brevis in the forearm to compare the muscles among primates (Cebus (Ce) (now Sapajus), baboons (Ba), chimpanzees (Ch), and modern humans (Hu)), and Ce was used as the control species . According to the data by Aversi-Ferreira et al. , the characteristics of this muscle are shown in Table 2.
Firstly, the control species was Cebus (Ce), and the was calculated for the muscle extensor pollicis brevis using parameters of innervation, origin, insertion, and vascularization. Eight samples of Ce were used; therefore, 16 muscles were analyzed. All specimens contained the muscle extensor pollicis brevis, and all studied muscles showed same innervation and vascularization , identical origin (), and insertion ().
For the innervation of the muscle in Ce, relative frequency is as follows:
The same analysis must be done for each parameter in each muscle of different species. It is noted that the parameters in primates other than Ce were obtained from previous literatures [15, 17, 18]. The insertion in Hu does not match with that in Ce (Table 2). Then, the for the parameter insertion is zero. In Ba , every will be zero, because this muscle is absent in this species.
The next step is to calculate the in each species. The of Ce is as follows:
In Hu, the is zero, and then the is as follows:
The for Ch () is equal to . In Ba (), since the extensor pollicis brevis is absent, is zero. The conclusion is , , and .
The next step is to calculate the CAI between the control species and others, one by one, as follows:
The results by the CAI calculation represent degree of difference in the extensor pollicis brevis among the species. The CAI values suggest that Ce and Ch have identical extensor pollicis brevis (high similarity), while Ce and Hu show some differences (somewhat similar), and Ce and Ba are highly different. It is because the extensor pollicis brevis does not exist in Ba or that muscle have no match between Ce and Hu. These results provide quantitative measures with previous papers that described the differences in the extensor pollicis brevis among different species.
It is also important to verify usefulness of the GCAI. We reanalyzed the previous data  and applied CAI to individual deep flexor muscles in the forearm (Table 3) and GCAI to the group of the deep flexor muscles in the forearm (Table 4). Although CAIs indicate the difference in individual muscles, it is unknown how the deep flexor muscles in the forearm as a group are different between the species. The GCAI values suggest that the flexor deep muscles in the forearm of the Cebus (Ce) are more similar to, in order, chimpanzee, modern humans, and baboon (Table 4).
|Note. : proportion of normal structured organs, where represents individual species, represents individual structures (i.e., muscles), and represent individual parameters of the muscles.|
4.2. Comparison of the SMCA with Parametrical Method (-Test)
A previous study analyzed the palmaris longus among 9 nonhuman primates and humans  by calculating ratios of muscle length to its tendon length, and these data were submitted to -test. Additionally, the characteristics of these muscles were compared among the 9 species by the SMCA (Table 5) (data from Aotus and modern humans were not used; for detailed analysis, see ). Using CAI (SMCA), nominal variables of features in the palmaris longus including origin, insertion, innervation, muscle presence, and type of belly arrangement were analyzed.
The results by parametric analysis of the ratios (muscle length/tendon length) are consistent with the primate grouping (distinction of apes, old world primate, new world primates, and prosimii); within the same group, the ratios were not significantly different (Table 5). These data are further consistent with those by CAI (SMCA) analysis of the palmaris longus (Table 5). These parametric data are at least partially consistent with those by the nonparametric (CAI) comparison (Table 6).
4.3. Comparison of the SMCA with Other Nonparametric Methods
Another possibility to study nominal variables is to use the cladistics method. This method considers the binary characters and any other possibility could be an error, because these features are mutually excluded . Indeed this characteristic limits its application to morphological analyses of structures since it considers just two parameters: 0 for absent characteristic and 1 for its presence. Nevertheless, this method is important in evolutionary studies, since this method might provide a different concept; the cladistics analysis prioritizes the primitive and derivative features , while the morphological analysis prioritizes utmost characters observed in a structure.
CAI in Table 7 indicated that the pronator quadratus is similar among Ce, Ch, and Ba and different between Ce and Hu. The cladistics analysis also indicated the same results. For the flexor digitorum profundus, CAI indicated that this muscle is similar between Hu and Ce, somewhat similar between Ch and Ce, and completely different between Ba and Ce. The cladistics analysis also indicated the same characteristics; Hu, Ch, and Ce shared the same features of this muscle but were totally different from Ba. For the flexor pollicis longus, CAI indicated this muscle is identical in Ce and Ch, somewhat similar between Hu and Ce, and highly different between Ba and Ce. The cladistics analysis indicated that this muscle was similar between Ce and Ch and different between Hu and Ce and between Ba and Ce. Thus, the results by both methods are consistent. Nevertheless, CAI provides quantitative data to assess the features. Furthermore, GCAI provides further information demonstrating that the group of these muscles of Ce was more similar to Ch, Hu, and Ba in this order (Table 7).
Firstly, an accurate analysis is fundamental in comparative anatomy. Photographs of structures can be used to verify the new data mainly when the new data disagree with others. However, numerical methods can help to avoid subjectivities and mistakes and make it easier to assess group similarities and to observe differences among structures and specimens. The application of numerical methods to comparison of biological structures could avoid prolix texts and ensure more exact and precise conclusions on similarities or differences among biological structures. This is especially true for gross anatomy where the nominal variables are studied and described.
Nonparametric statistics are less exact than parametric ones . Numerical parametric analyses can be performed with observation of frequency of characteristics based on fundamental concepts of normality and variation in anatomy. However, the nonparametric methods take advantage when parametric methods cannot be used in situations where nominal variables are included in the analysis and the measured data do not match nonparametric one. That is, parametric statistics are preferable when continuous variables are available , while nonparametric methods could be applied to nominal variables. As shown in Table 5, both methods match partially when the same samples were analyzed by both methods .
The analyses by nonparametric methods (chi-square and cladistics) indicate that cladistics might be more applicable to anatomical analysis. However, this method cannot analyze frequency of features. Furthermore, chi-square method also has a problem in its basic premises; causality of data used in this method disagrees with evolution theory assuming a common ancestral , and it is not easy to obtain numerical nominal data from expected and observed data. Although the results by cladistics and CAI are consistent (Table 7), it is noted that CAI provides quantitative data to assess features. Furthermore, GCAI provides additional information of a group of multiple structures.
The objective of cladistics method is different from that of SMCA; cladistics method offers solution when data from large groups are analyzed in evolution analysis, while the goal of SMCA is to obtain utmost details from anatomical structures when multiple features are analyzed. The cladistics is used in evolution analysis, especially, to obtain genealogic tree. However, it can be applied to comparative anatomy, because observation of presence and absence of a specific feature provides a general idea of similarity or difference among the structures.
In conclusion, the SMCA can offer more exact and precise method to compare structures to assess same or different groups, which requires analyses of nominal features in comparative anatomy.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
The authors thank the anonymous reviewers for their wise, detailed, and very helpful comments on previous version of this paper. R. A. G. M. F. Aversi-Ferreira was recipient of a doctoral fellowship from CAPES/Brazil and T. A. Aversi-Ferreira is recipient of Scholarship Research Productivity from National Council of Technology and Development (CNPq/Brazil). This work was supported partly by the Japan Society for the Promotion of Science Asian Core Program and the Ministry of Education, Science, Sports and Culture, Grant-in-Aid for Scientific Research (B) (25290005). The authors declare no competing financial interests.
- G. D. A. Pfrimer, T. de Abreu, V. S. Vieira et al., “Historic and teaching aspects of anatomy and cebus genus role in contemporary anatomy,” International Journal of Morphology, vol. 30, no. 2, pp. 607–612, 2012.
- T. A. Aversi-Ferreira, “A new statistical method for comparative anatomy,” International Journal of Morphology, vol. 27, no. 4, pp. 1051–1058, 2009.
- R. A. G. M. F. Aversi-Ferreira, T. de Abreu, G. A. Pfrimer et al., “Comparative anatomy of the hind limb vessels of the bearded capuchins (Sapajus libidinosus) with apes, baboons, and Cebus capucinus: with comments on the vessels' role in bipedalism,” BioMed Research International, vol. 2013, Article ID 737358, 15 pages, 2013.
- W. Thomson, “Electrical units of measurement,” in The Practical Applications of Electricity, The Institution of Civil Engeneers, London, UK, 1884.
- T. A. Aversi-Ferreira, R. S. Maior, F. O. Carneiro-e-Silva et al., “Comparative anatomical analyses of the forearm muscles of Cebus libidinosus (Rylands et al. 2000): manipulatory behavior and tool use,” PLoS ONE, vol. 6, no. 7, Article ID e22165, 8 pages, 2011.
- R. A. G. M. F. Aversi-Ferreira, R. V. Bretas, R. S. Maior et al., “Morphometric and statistical analysis of the palmaris longus muscle in human and non-human primates,” BioMed Research International, vol. 2014, Article ID 178906, 6 pages, 2014.
- R. A. C. Barros, I. L. S. Prada, Z. Silva, A. R. Ribeiro, and D. C. O. Silva, “Constituição do plexo lombar do macaco Cebus apella,” Brazilian Journal of Veterinary Research and Animal Science, vol. 40, no. 5, pp. 373–381, 2003.
- A. R. Ribeiro, I. L. S. Prada, Z. Silva, R. A. C. Barros, and D. C. O. Silva, “Origem do plexo braquial do macaco Cebus apella,” Brazilian Journal of Veterinary Research and Animal Science, vol. 42, no. 2, pp. 143–149, 2005.
- C. Darwin, “Afinidades mútuas dos seres vivos; morfologia; embriologia; órgãos rudimentares,” in The Origin of Species, John Murray, London, UK, 6th edition, 1876.
- Y. Takane and L. Zhou, “Anatomy of Pearson's Chi-square statistic in three-way contingency tables,” in New Developments in Quantitative Psychology, vol. 66 of Springer Proceedings in Mathematics & Statistics, pp. 41–57, Springer, New York, NY, USA, 2013.
- R. A. Pimentcl and R. Riggins, “The nature of cladistic data,” Cladisticis, vol. 3, no. 3, pp. 201–209, 1987.
- J. A. Hawkins, C. E. Hughes, and R. W. Scotland, “Scotland primary homology assessment, characters and character states,” Cladistics, vol. 13, no. 3, pp. 275–283, 1997.
- K. L. Moore, A. F. Dalley II, and A. M. R. Agur, “Intoduction to clinically oriented anatomy,” in Clinically Oriented Anatomy, Wolters Kluwer Health, Lippincott Williams & Wilkins, Philadelphia, Pa, USA, 7th edition, 2014.
- S. Gibbs, Comparative soft tissue morphology of the extant Hominoidea, including man [Ph.D. thesis], The University of Liverpool, Liverpool, UK, 1999.
- S. Standring, “Pelvis girdle and lower limb,” in Gray's Anatomy: The Anatomical Basis of Clinical Practice, Churchill Livingstone, London, UK, 2008.
- H. Gray, “Miology,” in Anatomy of the Human Body, Lea & Febiger, Philadelphia, Pa, USA, 1918.
- K. L. Moore, T. V. N. Persaud, and M. G. Torchia, Muscular System in the Developing Human: Clinically Oriented Embryology, Saunders, Philadelphia, Pa, USA, 9th edition, 2013.
- D. R. Swindler and C. D. Wood, In an Atlas of Primates Gross Anatomy: Baboon, Chimpanzee and Man, University of Washington Press, Seattle, Wash, USA, 1973.
- C. R. Kitchen, “Nonparametric versus parametric tests of location in biomedical research,” American Journal of Ophthalmology, vol. 147, no. 4, pp. 571–572, 2009.
- D. J. Sheskin, Handbook of Parametric and Nonparametric Statistical Procedures, Chapman & Hall/CRC, Boca Raton, Fla, USA, 3rd edition, 2004.
Copyright © 2015 Roqueline A. G. M. F. Aversi-Ferreira et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.