Research Article | Open Access
Shouhei Ueda, Taito Nozawa, Tetsuya Matsuzuki, Ryo-ichi Seki, Shinya Shimamoto, Takao Itino, "Phylogeny and Phylogeography of Myrmica rubra Complex (Myrmicinae) in the Japanese Alps", Psyche: A Journal of Entomology, vol. 2012, Article ID 319097, 7 pages, 2012. https://doi.org/10.1155/2012/319097
Phylogeny and Phylogeography of Myrmica rubra Complex (Myrmicinae) in the Japanese Alps
We investigated the genetic diversification of the mountain ant, Myrmica kotokui, in the Japanese Alps by using molecular phylogenetic analyses. Myrmica kotokui is widely distributed in Japan, and in the central Japanese Alps it is found only between elevations of approximately 1000 to 2000 m. We hypothesized that genetically distinct clades of this ant species might inhabit different mountain ranges in central Japan. To test this hypothesis, we reconstructed a molecular phylogeny using the DNA sequences of the mitochondrial cytochrome oxidase I gene and the nuclear long-wavelength rhodopsin gene of M. kotokui specimens collected from six mountain ranges in the Japanese Alps. The phylogeny showed four highly differentiated clades. However, the correspondence between the clades and morphological species was a little confusing. Two clades were composed only of M. kotokui specimens, whereas the other two clades were composed of multispecies, suggesting the possibility of multispecies composition of putative M. kotokui. The distribution pattern of these clades did not support our hypothesis of geographical differentiation, because two were distributed across all ranges, and a third was distributed in five of the six ranges. On the other hand, we found a pattern in the altitudinal distribution of the clades: one clade was distributed only at higher elevations, and the others were distributed at lower elevations. Thus, the ant clades do not show geographical segregation by mountain range, but they do show altitudinal differences.
Modern molecular phylogenetic techniques have revealed that some morphological species are composed of several genetically distinct cryptic species . To evaluate the biodiversity of a taxonomic group, it is essential to identify cryptic species, and biogeographical studies that determine the distributions of cryptic species are of crucial importance in conservation biology. Although many cryptic species have been found in several ant genera by using molecular phylogenetic techniques [2–5], it is difficult to identify cryptic species by using traditional morphological classification techniques alone, because intraspecific morphological variation may be greater than the interspecific variation and because convergently evolved characters may not show any morphological variation among species [3, 5]. For instance, Schlick-Steiner et al.  used mitochondrial DNA (mtDNA) to reconstruct the phylogeny of a Tetramorium ant species complex. Although the species of this complex are difficult to distinguish morphologically, by using their mtDNA phylogeny, morphology, and cuticular hydrocarbons, Schlick-Steiner et al.  succeeded in identifying seven cryptic species.
In the genus Myrmica, about 180 species are known in the Holarctic region . The taxonomy of Myrmica is well understood in Palearctic region, and the Palearctic species are morphologically classified into some species groups . In Japan, currently nine species are known, but taxonomic understanding of Myrmica species in Japan is not yet satisfactory . Myrmica kotokui, which belongs to Myrmica rubra species group (Figure 1), is known from the Russian Far East, Korea, and Japan [9, 10] and is distributed in Japan from Hokkaido in the north to the Yakushima Islands in the south . This species, which probably originated in Eurasia [11, 12], is adapted to a cold-temperate climate. As a result, in lower latitudes, M. kotokui is restricted to high elevations. For example, on Mt. Norikura in the Japanese Alps, this species is mainly found at elevations from approximately 1000 to 2000 m . Thus, in the mountains of Japan, the habitat of M. kotokui is fragmented into several patches of mountain ranges , which has led us to hypothesize that M. kotokui might be composed of several genetically distinct clades, or cryptic species, each inhabiting a different mountain range. Moreover, this possibility is supported by the dimorphism and two social systems of M. kotokui that have been observed in Hokkaido. Populations of M. kotokui in Hokkaido have been divided into two groups according to their social structure, monogynous, or polygynous, and these two groups can also be differentiated by morphological characters of queens, namely, head width and the ratio of wing length or thorax width to head width [15–17]. Although queen body size of M. kotokui was significantly larger in the monogynous population than in the polygynous population, the size overlapped extensively between the two populations. Therefore, social structure of M. kotokui cannot be distinguished simply by queen body size .
The taxonomic classification of M. kotokui is controversial. Forel  described this species as a variety of M. ruginodis, but Collingwood  reclassified M. kotokui as a separate species. Onoyama  pointed out, however, that M. kotokui should be considered a subspecies of M. ruginodis because there are no remarkable morphological differences between M. ruginodis and M. kotokui. Moreover, in Japan, M. kotokui is likely to be confused with M. rubra . Wetterer and Radchenko  also reported that in early published records from East Asia, some M. kotokui specimens appear to have been misidentified as M. rubra. M. ruginodis and M. rubra are distributed through the northern Palaearctic region, and M. ruginodis lives at higher latitude and altitude than M. rubra . M. rubra, which is known as an alien species spreading through temperate North America , is infrequently discovered in Hokkaido .
Recently, Jansen et al.  reconstructed a broad-scale phylogeny of the genus Myrmica and confirmed that M. kotokui forms a genetically distinct species from M. ruginodis and M. rubra and that M. ruginodis is a closest relative of M. kotokui (also see ). Their phylogeny, however, was based on only two specimens of M. kotokui from Hokkaido; no M. kotokui specimens from Honshu were included. Therefore, to clarify the taxonomy of M. kotokui in Japan, a molecular phylogeny should be reconstructed using samples from both Honshu and Hokkaido.
For this study, we collected 72 ant colonies of M. kotokui in the Japanese Alps, central Honshu, and reconstructed a molecular phylogeny on the basis of the mitochondrial cytochrome oxidase I gene and the nuclear long-wavelength rhodopsin gene. Our aim was twofold: first, to clarify the phylogenetic position of M. kotokui among related taxa and second, to test our hypothesis that M. kotokui is composed of several genetic clades, each inhabiting a different mountain range.
2. Materials and Methods
2.1. Sampling and Specimens
Between June 4, and October 10, 2009, we collected 72 ant colonies of M. kotokui from 36 locations in six mountainous ranges of Honshu (Joshinetsu, Kita-Alps, Hijiri, Yatsugatake, Chuo-Alps, and Minami Alps) at elevations ranging from 900 to 1800 m a.s.l. Voucher specimens are deposited in the Faculty of Science, Shinshu University, Matsumoto, Japan. For the phylogenetic study, we also used as ingroups two sequences of M. kotokui published by Jansen et al. . For outgroups, we used the COI and LwRh sequences of 11 species used by Savolainen and Vepsäläinen , Jansen et al.  and Leppänen et al. : three species (M. arisana, M. rubra, and M. ruginodis) in the rubra group, which also includes M. kotokui; seven species (M. formosae, M. serica, M. yamanei, M. schoedii, Myrmica n. sp. 2, M. indica, and M. weberi) in the ritae group, a sister clade of the rubra group; one species (M. rugosa) of the rugosa group, a clade distant from the rubra group. In addition, for outgroups, we added four species (M. excelsa, M. jessensis and M. taediosa in lobicornis group; M. luteola in luteola group), which are distributed in Japan. Collection locations (latitude, longitude, and elevation) of the specimens and their GenBank accession numbers are listed in Table S1 (see Supplementary materials available online at http://dx.doi.org/10.1155/2012/319097).
2.2. DNA Extraction, Polymerase Chain Reaction, and Sequencing
DNA was extracted from the whole body of each ant with a DNeasy Blood & Tissue Kit (Qiagen, Hilden, Germany) following the manufacturer’s protocols. A mitochondrial COI gene and a nuclear LwRh gene were amplified by polymerase chain reaction (PCR) using Takara Ex Taq (Takara Bio, Shiga, Japan). The PCR primers for COI were mtD-6 (5′-GGA TCA CCT GAT ATA GCA TTC CC-3′) and Nancy (5′-CCC GGT AAA ATT AAA ATA TAA ACT TC-3′) , and that for LwRh were LR143F (5′-GAC AAA GTK CCA CCR GAR ATG CT-3′)  and LR672R (5′-CCR CAM GCW GTC ATG TTR CCT TC-3′) . The PCR temperature profile used for COI was 30 cycles of 95°C for 30 s, 50°C for 30 s, and 72°C for 40 s and that used for LwRh was 30 cycles of 98°C for 10 s, 55°C for 30 s, and 72°C for 40 s. After amplification, the PCR products were purified with ExoSap-IT (USB, Cleveland, OH, USA). Cycle sequencing reactions for both strands were performed with a BigDye Terminator v1.1 Cycle Sequencing Kit (ABI, Weiterstadt, Germany) on an ABI 3130 Genetic Analyzer.
2.3. Sequence Alignment and Character Statistics
Mitochondrial COI and nuclear LwRh sequences were edited and aligned with SeqScape v. 2.5 (ABI, Weiterstadt, Germany). Base frequency homogeneity was tested by using the test in the PAUP* 4.0b10 . Parsimony-uninformative sites were excluded from the test. The test did not reject the hypothesis of homogeneity of nucleotide frequencies in each pair of taxa (; Table ). To test conflicts in phylogenetic signal among each dataset we conducted an incongruence length differences (ILDs) test (Farris et al., 1994) in PAUP* 4.0b10 with heuristic searches of 100 replicates with tree bisection and reconnection (TBR) and 10 random addition replicates. The ILD test revealed no conflict between COI and LwRh (). The degree of substitution saturation in the third codon position of the COI and LwRh sequences was assessed by plotting the transitions (Tis) and transversions (Tvs) ratios against genetic distance for each data set with the DAMBE software package by the method of Xia and Xie  (Figure ). In the saturation plot analysis, the simple JC69 substitution model  was used instead of the J2 substitution model (COI)  or the HKY substitution model (LwRh) , which was selected in the model selections described in the next paragraph, because DAMBE does not support the J2 and the HKY substitution model. Substitution saturation at the third codon position was not detected (; Figure ). In addition, we did not find the evidence of mitochondrial pseudogenes, called nuclear mitochondrial transfer (numts), in the mitochondrial COI sequences, which leads to the erroneous phylogeny : there are not indels and stop codons in COI sequences. Following the results of the character statistics, we used all data sets and all positions in each data set for the phylogenetic analyses.
2.4. Phylogenetic Analyses
Best-fit substitution models were selected for each codon position of each gene by using Bayesian information criterion 5 (BIC5) and the Kakusan4 program  (Table ). Maximum likelihood (ML) analysis was performed with TREEFINDER version October 2008  and the models selected by Kakusan4. Clade support was assessed by 1000 bootstrap replications in TREEFINDER. In addition, Bayesian posterior probabilities and maximum parsimony (MP) bootstrap support were obtained with MrBayes version 3.1.2  and PAUP* 4.0b10 software , respectively. The models selected by BIC5 were also used in the Bayesian analysis (Table ) with the default run settings, in which two independent analyses are performed, each with four chains (one cold and three heated). The Bayesian analysis was run for 1 million generations, with sampling every 1000 generations. We assessed the log likelihood of each sampling point against generation time to identify when the Markov chains reached a stationary distribution and accordingly discarded the initial 2000 trees as burn-in. The parsimony bootstrap support was assessed with 1000 bootstrap replicates by using heuristic searches with tree bisection and reconnection (TBR) and 100 random addition replicates for each.
2.5. Statistical Analysis
Differences in the elevation range among the distributions of the ant clades were analyzed by one-way analysis of variance (ANOVA) and Tukey’s multiple comparison test in the R software package .
3. Results and Discussion
3.1. Taxonomy of M. kotokui Based on DNA Phylogeny
Molecular phylogenies of M. kotokui were inferred from a total of 855 bp from the two genes (508 bp sequences of the mitochondrial COI gene and 347 bp sequences of the nuclear LwRh gene). Although there is no conflict of the phylogenetic signal between the mitochondrial COI gene and the nuclear LwRh gene, the further analyses of additional nuclear gene sequences will be required to reveal the conflict between mtDNA and nuclear DNA phylogenies due to introgression or incomplete lineage sorting [34, 35]. It is because the substitution rate of LwRh is much lower than that of COI and genetic divergences between L1 and L2 and between L3 and L4 were not detected in the analysis of the nuclear LwRh marker alone (data not shown).
Four independent clades were identified by the ML, MP, and Bayesian analyses (L1 to L4; Figure 2), suggesting that the single morphological species known as M. kotokui is composed of several putative cryptic species. However, the correspondence between the clades and morphological species was a little confusing. Two clades (L3 and L4) were composed only of M. kotokui specimens, whereas the other two clades composed of multispecies: L1 was composed of M. rubra and M. kotokui specimens, and L2 was composed of M. arisana and M. ruginodis as well as M. kotokui specimens (Figure 2).
The possibility of multispecies composition of M. kotokui may account for the confused taxonomy of M. kotokui in Japan as reported by several studies [10, 15, 19, 36]. The absence of morphological differences between M. kotokui and M. ruginodis reported by Onoyama , for example, might be because M. kotokui specimens of the L2 clade, which is closely related to M. ruginodis (Figure 2), were used for the morphological comparison of the two species. Moreover, the misidentification of M. kotokui specimens as M. rubra [10, 19] may have occurred because the L1 clade of M. kotokui is closely related to M. rubra (Figure 2).
Taken at face value, it is possible that the L1 and L2 of putative M. kotokui are misidentification of M. rubra and M. ruginodis, respectively, and that the L3 or L4 clades, or both, should be regarded as the true species M. kotokui. Here, the taxonomy of M. kotokui cannot be clarified because detailed morphological comparisons have not been made with type material. All specimens collected from the Japanese Alps for this study were identified as M. kotokui by using the taxonomic key for Japanese ant species of Morisita et al. , which was made based on the original species description of Collingwood , but further morphological study of the specimens is desirable. Further multidisciplinary studies of the morphology, social structure, and cuticular hydrocarbons of M. kotokui and its relatives [4, 5] are also necessary to solve these taxonomic problems.
3.2. Geographical Distribution and Altitudinal Distribution
We investigated the geographical distribution of each clade to determine whether they were allopatrically isolated. We found that L1 and L3 were distributed across all six mountain ranges, L2 was distributed in all mountain ranges except Yatsugatake, and L4 (present in only two samples) was found only in Yatsugatake (Figure 3). Thus, the distribution pattern of the ant clades did not support our hypothesis of geographical genetic differentiation by mountain range. Therefore, we investigated the altitudinal distributions of the ant clades in the Japanese Alps between 900 and 1800 m a.s.l. and found that the L3 clade was distributed at a significantly higher elevation range than L1 and L2 (Figure 4; one-way ANOVA, ; Tukey’s multiple comparisons: L1 × L3, ; L2 × L3, ). A significant difference was not detected between the L4 distribution and the other distributions because of the small L4 sample size.
These results raise the following question. Did the genetic differentiation of the ant clades occur in the Japanese Alps or did clades that had already diverged on the Eurasian continent migrate to the Japanese Alps? The differentiation of the clades may not have occurred in the Japanese Alps for two reasons: (1) no physical barrier separates altitudinal ranges in the Japanese Alps and (2) the genetic distances between the clades are too large (6.5%) to have been accumulated after the last glacial maximum. Thus, the altitudinally stratified distribution of the ant clades may reflect multiple migrations from the continent with clades adapted to different levels of coldness becoming segregated in different habitats because of the drastic environmental differences along the elevation gradient . Taken together, our results suggest that a cryptic highland clade of Myrmica exists. We also propose that further biodiversity may be hidden in ants that live at high elevations.
The authors thank T. Komatsu, Y. Kataoka, T. Yamamoto, and T. Ando for sampling support; H. Mohri, S. Hiraga for sampling advice. They thank the Chubu Regional Office for Nature Conservation for permission of conducting our survey in the Japanese Alps. This research was funded by the Nagano Society for the Promotion of Science. This study was supported by Research and Education Funding for Japanese Alps Inter-Universities Cooperative Project, MEXT, Japan.
Supplementary Table S1: Specimens used for the phylogenetic analyses.
Supplementary Table S2: Table S2. Character statistics of the data sets.
Supplementary Table S3: Substitution models selected for each data set
Supplementary Figure S1: Plots of the transition and transversion ratios in the third codon position against genetic distance.
- D. Bickford, D. J. Lohman, N. S. Sodhi et al., “Cryptic species as a window on diversity and conservation,” Trends in Ecology & Evolution, vol. 22, no. 3, pp. 148–155, 2007.
- C. Bernasconi, D. Cherix, B. Seifert, and P. Pamilo, “Molecular taxonomy of the Formica rufa group (red wood ants) (Hymenoptera: Formicidae): a new cryptic species in the Swiss Alps?” Myrmecological News, vol. 14, pp. 37–47, 2011.
- K. G. Ross, D. Gotzek, M. S. Ascunce, and D. D. Shoemaker, “Species delimitation: a case study in a problematic ant taxon,” Systematic Biology, vol. 59, no. 2, pp. 162–184, 2010.
- B. C. Schlick-Steiner, F. M. Steiner, K. Moder et al., “A multidisciplinary approach reveals cryptic diversity in Western Palearctic Tetramorium ants (Hymenoptera: Formicidae),” Molecular Phylogenetics and Evolution, vol. 40, no. 1, pp. 259–273, 2006.
- B. Seifert, “Cryptic species in ants (Hymenoptera: Formicidae) revisited: we need a change in the alpha-taxonomic approach,” Myrmecological News, vol. 12, pp. 149–166, 2009.
- B. Bolton, G. Alpert, P. S. Ward, and P. Nasrecki, Bolton's Catalogue of Ants of the World, Harvard University Press, Cambridge, Mass, USA, 2006.
- A. G. Radchenko and G. W. Elmes, Myrmica Ants (Hymenoptera: Formicidae) of the Old World, Natura Optima dux Foundation, Warsaw, Poland, 2010.
- Japanese Ant Database Group, Japanese Ant Image Database 2008, Japanese Ant Database Group, Sendai, Japan, 2008, (in Japanese).
- C. A. Collingwood, “Ants (Hymenoptera, Formicidae) from North Korea,” Annales Historico-Naturales Musei Nationalis Hungarici, vol. 68, pp. 295–309, 1976.
- J. K. Wetterer and A. G. Radchenko, “Worldwide spread of the ruby ant, Myrmica rubra (Hymenoptera: Formicidae),” Myrmecological News, vol. 14, pp. 87–96, 2011.
- G. Jansen and R. Savolainen, “Molecular phylogeny of the ant tribe Myrmicini (Hymenoptera: Formicidae),” Zoological Journal of the Linnean Society, vol. 160, no. 3, pp. 482–495, 2010.
- G. Jansen, R. Savolainen, and K. Vepsäläinen, “Phylogeny, divergence-time estimation, biogeography and social parasite-host relationships of the Holarctic ant genus Myrmica (Hymenoptera: Formicidae),” Molecular Phylogenetics and Evolution, vol. 56, no. 1, pp. 294–304, 2010.
- R. Togai, S. Ueda, M. Hattori, and T. Itino, “The vertical distribution of ants in Mt. Norikura,” Japanese Journal of Environmental Entomology and Zoology, vol. 23, no. 3, pp. 119–126, 2012 (Japanese).
- M. Morisita, “Ari-rui (ants),” in Nihonseibutsushi, H. Furukawa, Ed., pp. 1–56, Kenkyusha, Tokyo, 1945.
- T. Kikuchi, “Between- and within-population morphological comparisons of all castes between monogynous and polygynous colonies of the ant Myrmica kotokui,” Ecological Entomology, vol. 27, no. 4, pp. 505–508, 2002.
- T. Kikuchi, S. Higashi, and T. Murakami, “A morphological comparison of alates between monogynous and polygynous colonies of Myrmica kotokui in Northernmost Japan,” Insectes Sociaux, vol. 46, no. 3, pp. 250–255, 1999.
- A. Mizutani, “On the two forms of the ant Myrmica ruginodis Nylander (Hymenoptera, Formicidae) from Sapporo and its vicinity, Japan,” Japanese Journal of Ecology, vol. 31, pp. 131–137, 1981.
- A. Forel, “Die Ameisen des K. Zoologischen Museums in München,” Sitzungsberichte der Bayerischen Akademie der Wissenschaften zu Muenchen, vol. 1911, pp. 249–303, 1911.
- K. Onoyama, “Confirmation of the occurrence of Myrmica rubra in Japan, with taxonomic and ecological notes (Hymenoptera, Formicidae),” Japanese Journal of Ecology, vol. 57, no. 1, pp. 131–135, 1989.
- B. Seifert, “A taxonomic revision of the Myrmica species of Europe, Asia Minor, and Caucasus (Hymenoptera, Formicidae),” Abhandlungen und Berichte des Naturkundemuseums Görlitz, vol. 62, pp. 1–75, 1988.
- R. Savolainen and K. Vepsäläinen, “Sympatric speciation through intraspecific social parasitism,” Proceedings of the National Academy of Sciences of the United States of America, vol. 100, no. 12, pp. 7169–7174, 2003.
- J. Leppänen, K. Vepsäläinen, and R. Savolainen, “Phylogeography of the ant Myrmica rubra and its inquiline social parasite,” Ecology and Evolution, vol. 1, no. 1, pp. 46–62, 2011.
- C. Simon, F. Frati, A. Beckenbach, B. Crespi, H. Liu, and P. Flook, “Evolution, weighting, and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers,” Annals of the Entomological Society of America, vol. 87, no. 6, pp. 651–701, 1994.
- P. S. Ward and D. A. Downie, “The ant subfamily Pseudomyrmecinae (Hymenoptera: Formicidae): phylogeny and evolution of big-eyed arboreal ants,” Systematic Entomology, vol. 30, no. 2, pp. 310–335, 2005.
- D. L. Swofford, PAUP∗: Phylogenetic Analysis Using Parsimony. Ver. 4.0b10, Sinauer Associates, Sunderland, Mass, USA, 2002.
- X. Xia and Z. Xie, “DAMBE: software package for data analysis in molecular biology and evolution,” Journal of Heredity, vol. 92, no. 4, pp. 371–373, 2001.
- T. H. Jukes and C. R. Cantor, “Evolution of protein molecules,” in Mammalian Protein Metabolism, H. N. Munro, Ed., pp. 21–132, Academic Press, New York, NY, USA, 1969.
- G. Jobb, A. von Haeseler, and K. Strimmer, “TREEFINDER: a powerful graphical analysis environment for molecular phylogenetics,” BMC Evolutionary Biology, vol. 4, article 18, 2004.
- M. Hasegawa, H. Kishino, and T. A. Yano, “Dating of the human-ape splitting by a molecular clock of mitochondrial DNA,” Journal of Molecular Evolution, vol. 22, no. 2, pp. 160–174, 1985.
- J. E. Buhay, “‘Coi-like’ sequences are becoming problematic in molecular systematic and DNA barcodeng studies,” Journal of Crustacean Biology, vol. 29, no. 1, pp. 96–110, 2009.
- A. S. Tanabe, “KAKUSAN: a computer program to automate the selection of a nucleotide substitution model and the configuration of a mixed model on multilocus data,” Molecular Ecology Notes, vol. 7, no. 6, pp. 962–964, 2007.
- J. P. Huelsenbeck and F. Ronquist, “MRBAYES: bayesian inference of phylogenetic trees,” Bioinformatics, vol. 17, no. 8, pp. 754–755, 2001.
- R. Ihaka and R. Gentleman, “R: a language for data analysis and graphics,” Journal of Computational and Graphical Statistics, vol. 5, no. 3, pp. 299–314, 1996.
- C. R. Linnen and B. D. Farrell, “Mitonuclear discordance is caused by rampant mitochondrial introgression in Neodiprion (Hymenoptera: Diprionidae) sawflies,” Evolution, vol. 61, no. 6, pp. 1417–1438, 2007.
- T. Sota and A. P. Vogler, “Incongruence of mitochondrial and nuclear gene trees in the Carabid beetles Ohomopterus,” Systematic Biology, vol. 50, no. 1, pp. 39–59, 2001.
- T. Kikuchi, F. Tomizuka, and S. Higashi, “Reproductive strategy in orphaned colonies of Myrmica kotokui forel, the Japanese species of the M. ruginodis complex (Hymenoptera, Formicidae),” Insectes Sociaux, vol. 47, no. 4, pp. 343–347, 2000.
- M. Morisita, M. Kubota, K. Onoyama et al., A Guide for the Identification of Japanese Ants (III). Myrmicinae and Supplement to Leptanillinae. (Hymenoptera: Formicidae), The Myrmecological Society of Japan, Tokyo, Japan, 1992.
- A. S. Jump, C. Mátyás, and J. Peñuelas, “The altitude-for-latitude disparity in the range retractions of woody species,” Trends in Ecology & Evolution, vol. 24, no. 12, pp. 694–701, 2009.
Copyright © 2012 Shouhei Ueda 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.