Peer-Review Report by Denise Crampton (Reviewer R) : https://bio.jmirx.org/2025/1/e85566/
Peer-Review Report by Osamu Kondo (Reviewer AE): https://bio.jmirx.org/2025/1/e85567/
Authors' Response to Peer-Review Reports: https://bio.jmirx.org/2025/1/e85565/
Abstract
Background: The Jōmon period (14,000 to 800 cal BC), when people on the Japanese archipelago engaged mainly in hunting and gathering and experienced relatively more severe climate changes, is an important prehistoric period for investigating how people reacted to environmental fluctuations in human evolutionary history. Anthropologists have extensively discussed the population history of the Jōmon period of the Japanese archipelago via their morphological variations. Some have supported the notion of relative morphological uniformity within the Jōmon population, which could be sustained by widespread population interactions, although others have claimed that spatiotemporal morphological differences (especially geographical clines) may exist to some extent.
Objective: The aim of this study is to examine the morphological interphase and interregion variations among the Jōmon populations to investigate the population interactions during this period (ie, how widely and continuously they interacted).
Methods: This study used geometric morphometrics of a much larger sample of 3D data of Jōmon human crania than previous studies (n=363 from 97 sites, including 146 females, 215 males, and 2 unknown-sex individuals). The configurations of landmarks were processed with generalized Procrustes analysis and principal component analysis. The principal scores were statistically analyzed with the Steel-Dwass test. We also compared the Jōmon crania with the Yayoi crania in the same way.
Results: Notably, the results of statistical tests on the Jōmon crania show that principal component (PC) 1 did not show any significant differences both in spatial and temporal comparisons. Regional differences did not show geographical clines in PC1. PC2 indicated that statistically significant differences were found in some regions (eg, the Tohoku and Tokai regions: Z=6.375, P<.001; the Kanto and Tokai regions: Z=4.880, P<.001), and a gradual geographical cline was found among the different regions and in phases (eg, the Early and Final phases: Z=3.118, P=.02; the Middle and Final phases: Z=4.233, P<.001). Comparative results between the Jōmon and the Yayoi populations also showed that the Jōmon populations were spatiotemporally less varied than the Yayoi populations and that individual variation within a site was more variable in the Jōmon site than that of the Yayoi site.
Conclusions: This observation is consistent with the possibility that the population interactions of the Jōmon people had been widespread and continuous, which has an important implication for their resilience against severe climate changes at that time. It is possible that the relative stability of the Jōmon society was sustained by their frequent interactions with various populations, as suggested by insights from relevant archeological, ethnographic, and genetic research.
doi:10.2196/72432
Keywords
Introduction
General Background
On the Japanese archipelago, the Jōmon period (14,000-800 cal BC; Incipient: 14,000-9300 cal BC; Initial: 9300-5200 cal BC; Early: 5200-3400 cal BC; Middle: 3400-2400 cal BC; Late: 2400-1200 cal BC; Final: 1200-800 cal BC []; the Incipient period was excluded in this research because no human skeletal remains from this period were found), when people engaged mainly in hunting and gathering, is prehistoric period that is particularly important for researchers investigating how people reacted to environmental fluctuations in human evolutionary history. This period experienced more severe climatic changes, including the Younger Dryas, than did subsequent Yayoi (800‐250 AD) and Kofun (250‐600 AD) societies [-]. Many have claimed that prehistoric societies undergoing radical environmental fluctuations tended to be accompanied by societal change [-]. The Jōmon society may be no exception. Research has suggested that climate change has affected population size in many regions [,], which may have influenced other aspects of the Jōmon society.
However, features such as the emergence of social hierarchies, shifts in subsistence patterns, and significant technological developments commonly seen in other periods were less radical in the Jōmon society. Jōmon people predominantly adhered to a hunter-gatherer lifestyle and persisted in using stone tools for over 10 millennia. Although recent findings at the Usu-Moshiri site in Hokkaido found evidence of violent injuries in more than 5 individuals, such occurrences remained relatively rare within the broader context of the Jōmon period. This stands in stark contrast to the following Yayoi period [,-] and can be distinguished from patterns observed in other relevant regions worldwide [-]. Although the notion of a “static view” of the Jōmon society has been challenged over the years [,], it remains clear that societal changes during this era were relatively less drastic, even in the face of more drastic climatic variations.
One potential strategy for mitigating the impact of external environmental stresses is sustained interaction between different populations. The continuous exchange of genetic material, knowledge, skills, and resources can enhance the resilience and adaptability of populations in the face of environmental disruptions. Many relevant studies of genetic diversity in organisms other than humans support this possibility, including plants, various animal species, and whole ecosystems [-]. The specific mechanisms behind the genetic diversity and resilience remain controversial, although it is generally argued that genetic diversity is somehow related to morphological and functional diversity, and that damage to populations due to ecological disturbance could be compensated for by such diversity.
This is also evident in modern hunter-gatherer societies, which exhibit higher genetic diversity, indicative of ongoing interactions and gene flow between groups [-]. Moreover, ethnographic research shows that mobile hunter-gatherer communities with lower population densities often seek mates from distant groups [-]. Archeological and anthropological research also suggest that higher mobility of prehistoric hunter-gatherers was an important factor for their societal stability [,,].
Previous Research on Jōmon Population Interactions
Anthropologists have extensively discussed the population interaction history during the Jōmon period of the Japanese archipelago. Typically, previous research used a traditional biodistance method using distances between specific measuring points[,] (see also ), and statistically significant differences between morphological biodistances among different populations have been regarded as populational and genetic differences. For example, as stated in the next paragraph, if geographical clines and statistically significant differences are found among morphological biodistances of different populations, it is commonly suggested that populations moved in a certain geographical direction and that they did not migrate widely.
Some research on human skeletal remains including crania and limbs has supported the notion of relative morphological uniformity within the Jōmon population, which could be sustained by widespread population interactions, although others have claimed that spatiotemporal differences may exist to some extent [-]. Dodo [,] argued that statistical differences found in cranial biodistances among Jōmon populations in the Tohoku, Kanto, and western regions were much smaller than among Jōmon populations and more recent populations. On the other hand, it has often been claimed that geographical clines are found from north to south because the Jōmon people were assumed to be influenced by both the southern and northern populations outside the Japanese archipelago. Fukase et al [] claimed that such geographical clines were found in the length of limbs. Hanihara and Ishida [] and Kondo [] argued that geographical clines were also found in cranial biodistances. Ogata [] and Yamaguchi [,] showed that the body size of the Incipient and Early Jōmon people was smaller than that of the other subperiods due to the dietary conditions.

Many archeological studies have also discussed population interactions via a widespread distribution and exchange of artifacts such as pottery and lithics over extensive geographic regions during this period [,,-]. It has often been pointed out that a certain type of pottery was widely distributed across regions. An example is the Ento-kaso type distributed in both the Hokuriku and the northern Tohoku regions of the early phase (see for regions in the Japanese archipelago). Even if the distribution of a type of pottery was regionally restricted, it may have been affected by a distant region. The Funamoto type of the Middle phase was distributed mainly in the Sanyo regions, and this type of pottery strongly influenced the pottery in the Kyushu region during that same phase. In the Late phase, when samples were obtained from the largest number of different regions in this study, it was found that a large quantity of the Horinouchi type I pottery originating from the Kanto region was excavated in the Kinki region; similarly, the Nakatsu type pottery from the Kinki region was transported to the Kanto region and influenced the origin of Shomyoji pottery [,]. Such diverse and widespread distribution of pottery suggests that the Jōmon people interacted widely, although it has often been argued that Jōmon society was regionally diverse and distinct, with differences in settlement styles, lifestyles, and types of archeological remains, such as clay figurines. This is why the Jōmon cultural complex has been referred as “Jōmon cultures” [,,].

The geometric morphometrics of 3D data, the approach used in this research, which explores shape variations of targeted objects typically using the coordination of landmarks, has gained traction in various fields, extending to objects such as lithics [-]. Although the morphological variations of the Jōmon people have been explored in previous studies using the traditional biodistance method, this study has the following novelties. The 3D data of Jōmon crania have rarely been examined using geometric morphometrics, with some exceptions [-], and these exceptions focused on smaller samples from restricted regions. Less biased larger samples are important for estimating macroscopic population interactions. Furthermore, when analyzing the resulting data, geometric morphometrics is able to track morphological variation among configurations of each measuring point at once, whereas the traditional biodistance method compares each measured distance independently. These differences could provide some new insights.
The Aim of This Research
This study aims to examine the hypothesis that the Jōmon populations interacted widely and continuously in an anthropological way (ie, through geometric-morphological analysis of the cranial features of the Jōmon people) using a much larger 3D cranial sample collected from broader regions than previous studies. If the populations interacted widely and continuously, there should be fewer morphological differences among different regions and phases (ie, fewer statistically significant differences between each region and phase are found). Moreover, when the Jōmon populations are compared with populations from a different period, the hypothesis would expect that interregion and interphase variations should be lower than in the populations from a different period.
Methods
Landmark Selection for Geometric Morphometric Examinations
We chose 31 representative and/or less complicated landmarks for this geometric morphometric study, which were placed on 3D models using the geomorph package in R (version 4.3.2; R Foundation for Statistical Computing; ). Our research focused on the human cranium, enabling us to use the landmark method, which was used in similar studies [-]. Landmark selection was based on those described by Martin [] and Caple and Stephan [] (see for the Bookstein type of landmark [] selected in this study).
3D Data Collection
We collected relatively well-preserved and less skewed 3D cranial data from the Jōmon period in collaboration with local archeological centers, museums, and universities. We used a range of scanning technologies, specifically Creaform HandySCAN BLACK, HandySCAN BLACK Elite, and Einscan HD, in addition to Structure from Motion and Multi-View Stereo techniques based on 2D photographs [,]. Previous studies have confirmed the consistency of the models generated by the aforementioned methods for human crania [,]. Our dataset comprised 363 Jōmon crania (including 146 females, 215 males, and 2 individuals with unknown sex), with specific information available in and . These 3D data were collected from various regions spanning the island of Honshu within the Japanese archipelago, encompassing regions from Tohoku to Kyushu (). Although buried crania from the prehistoric period tend to be distorted due to soil conditions, we chose to measure crania that are regarded as less skewed and better preserved based on their appearance and relevant excavation information.
| Phases | Initial | Early | Middle | Late | Final | Unknown | Yayoi | Total | ||||||||||||||||
| Area | F | M | U | F | M | U | F | M | U | F | M | U | F | M | U | F | M | U | F | M | U | F | M | U |
| Tohoku | 0 | 0 | 0 | 1 | 1 | 0 | 9 | 14 | 0 | 6 | 14 | 0 | 4 | 4 | 0 | 6 | 5 | 0 | 0 | 0 | 0 | 26 | 38 | 0 |
| Kanto | 0 | 1 | 0 | 3 | 1 | 0 | 10 | 20 | 0 | 14 | 35 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 28 | 57 | 0 |
| Hokuriku | 1 | 1 | 0 | 2 | 6 | 0 | 0 | 2 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 3 | 10 | 0 |
| Tokai | 3 | 3 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 4 | 11 | 0 | 42 | 48 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 49 | 63 | 1 |
| Kinki | 2 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 1 | 1 | 0 | 0 | 0 | 0 | 4 | 2 | 0 |
| Sanin | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 5 | 13 | 0 | 5 | 13 | 0 |
| Sanyo | 0 | 1 | 0 | 0 | 1 | 0 | 7 | 10 | 0 | 4 | 3 | 0 | 15 | 16 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 26 | 32 | 0 |
| Shikoku | 2 | 2 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 3 | 3 | 1 |
| Kyushu | 0 | 0 | 0 | 0 | 1 | 0 | 1 | 1 | 0 | 2 | 6 | 0 | 0 | 0 | 0 | 4 | 2 | 0 | 4 | 18 | 0 | 11 | 28 | 0 |
| Total | 8 | 8 | 1 | 6 | 11 | 0 | 27 | 48 | 0 | 32 | 71 | 0 | 62 | 68 | 1 | 11 | 9 | 0 | 0 | 0 | 0 | 155 | 246 | 2 |
aF: female.
bM: male.
cU: sex unknown.
To provide a basis for comparison, we included cranial data from the Kuma-Nishioda site of the Middle Yayoi period (350 cal BC to AD 30) in the Kyushu region and the Aoyakamijichi site of the Late Yayoi period (AD 30 to 250) in the Sanin region, comprising 22 and 18 crania, respectively. The choice of the Kuma-Nishioda site was motivated by the common assumption that Yayoi people, especially in the northern Kyushu region, might exhibit lower genetic and morphological diversity due to a genetic bottleneck caused by gene flow from migrating Korean populations [,,,]. The Aoyakamijichi site was selected because it is the only one with a relatively large number of 3D samples outside of the Kyushu region []. Since we could not gather good 3D data from the Early Yayoi period, we did not include this period.
Data Processing
The file size of the original 3D data was so large that it took a long time to read the data and locate the landmarks in R; the files ranged from 100 to 600 MB depending on the measurement methods. To address this issue, we reduced the size to 2.5%‐15% of the original using MeshLab []. Previous studies have confirmed that this reduction process does not result in a model that is significantly different from the original [,,].
Some prehistoric crania were incomplete, so we attempted to reconstruct lost areas in 3D models. This reconstruction involved mirroring remaining portions along the median line, passing through the nasion, prosthion, and bregma [,,,,]. Since this mirroring could include some subjective biases, mirroring was limited to cases with smaller distortions and clearer original positions of remaining parts. In cases where mirroring was not possible, particularly when corresponding parts were entirely lost, we applied the estimate missing function, using the thin-plate spline method within the geomorph package (version 4.0.4) of R (version 4.2.1; R Foundation for Statistical Computing) and R studio (2022.07.0 + 548) [-].
Sex and age estimation of skeletal remains depended mainly on published excavation reports, where anthropologists estimated the sex and age. If any relevant information was not given in the reports, we made educated guesses by investigating pelves and crania (especially postcranial features and sizes of the zygomatic arch and mastoid process) for sex [], while the pubic symphyses [] and cranial sutures were used to estimate age [,].
Statistical Analysis
After conducting generalized Procrustes analysis to rotate and adjust the landmark configurations for all Jōmon and Yayoi crania, we performed principal component analysis (PCA) to investigate the spatiotemporal variations of these crania. We did not perform PCA for each phase or region because no region had samples from all phases, and no phase had samples from all regions. Additionally, we also performed the Steel-Dwass test, which is a nonparametric multiple comparison test with an expansion of the Mann-Whitney U test, to evaluate the principal component (PC) scores of Jōmon crania based on each phase and region. All statistical analyses were performed using R and R Studio.
We also conducted the PCA and the same statistical tests on all the Jōmon and Yayoi populations to clarify how varied the Jōmon populations were when compared to populations from a different period. Next, we selected cranial data from Jōmon sites with relatively larger sample sizes to compare with the Kuma-Nishioda site. They were compared via PCA by each site. The latter analysis included the Ebishima (from the Tohoku region of the Middle phase), Nakazuma (from the Kanto region of the Late phase), Yoshigo (from the Tokai region of the Final phase), and Tsukumo sites (from the Sanyo region of the Final phase), and their locations are indicated in and . Such comparisons of the Jōmon populations from different phases and regions with the Yayoi population could reveal the degree of variation of individual Jōmon populations.
Ethical Considerations
Ethical approval was not required for this study according to the local legislation and institutional requirements because under the Law for the Protection of Cultural Properties in Japan, human skeletal remains are regarded as cultural properties []. Curators of all samples are described in .
Results
Spatiotemporal Variations Among the Jōmon Populations
The PCA results indicated cumulative contribution rates exceeding 75% up to the 20th PC (). To streamline the discussion, this study primarily focused on PCs with magnitudes larger than 5% (PC1-PC5). PC1 was positively correlated with narrow facial width, prognathism (protrusion of the jaws), and high and short occipital areas, while PC2 was positively linked to temporal length and negatively related to facial length. PC3 primarily captured temporal length, and PC4 reflected facial height. PC5, primarily representing cranial asymmetry, was not included in the analysis as it was influenced by ground pressure during burial (). We also showed deformation patterns of a cranium from the Fujizuka shell midden (ID 341; see ) according to the morphological change each PC captured ().
| Contribution rate, % | Cumulative proportion, % | |
| PC1 | 11.1 | 11.1 |
| PC2 | 9.0 | 20.1 |
| PC3 | 6.4 | 26.5 |
| PC4 | 6.1 | 32.6 |
| PC5 | 5.4 | 38.0 |
| PC6 | 4.9 | 42.9 |
| PC7 | 4.1 | 47.0 |
| PC8 | 3.6 | 50.6 |
| PC9 | 3.0 | 53.6 |
| PC10 | 3.0 | 56.6 |
| PC11 | 2.7 | 59.3 |
| PC12 | 2.5 | 61.8 |
| PC13 | 2.3 | 64.1 |
| PC14 | 2.1 | 66.2 |
| PC15 | 2.0 | 68.2 |
| PC16 | 1.8 | 70.1 |
| PC17 | 1.7 | 71.8 |
| PC18 | 1.6 | 73.4 |
| PC19 | 1.5 | 74.9 |
| PC20 | 1.4 | 76.4 |
aPC: principal component.


The PCA results, graphically presented in and , suggested that the differences between phases are relatively small, in line with the statistical tests on PCs that indicated statistical differences only in the Early and Final (Z=3.118, P=.02), Middle and Final (Z=4.233, P<.001), Late and Final (Z=4.040, P=.001) in PC2, and Late and Final phase (Z=2.946, P=.038) in PC4 (see ). Notably, the scatter plots in and the box plots in showed more prominent regional differences, consistent with the statistical tests on PCs showing significant distinctions between various regions. These differences were particularly notable between the Tohoku and Tokai regions (Z=6.375, P<.001), the Tohoku and Sanyo regions (Z=5.852 P<.001), the Kanto and Tokai regions (Z=4.880, P<.001), the Kanto and Sanyo regions (Z=4.180, P<.001), the Hokuriku and Tokai regions (Z=3.773, P=.004), and the Hokuriku and Sanyo regions (Z=3.829, P=.003) in PC2. Additionally, there were significant differences in PC4 between the Kanto and Shikoku regions (Z=3.302, P=.02) and the Shikoku and Kyushu regions (Z=3.088, P=.04). It should be noted that there were no significant differences observed in PC1 and PC3, and the regional differences remained relatively limited. Furthermore, there were few sexual differences (see ), with some exceptions, specifically PC3 and PC4 in the early phase of the Tokai region and PC4 in the middle phase of the Kanto region, which was supported by statistical tests between the sexes in each PC score (U=15,926, P=.81 in PC1, U=13,925, P=.07 in PC2, U=12,693, P=.002 in PC3, and U=13,695, P=.004 in PC4 comparisons).


Comparisons Between the Jōmon Populations and the Yayoi Population
The two comparative results between the Jōmon and Yayoi populations were summarized in and and and . In the second site-to-site comparison, to streamline the presentation, these figures showed only PC1 and PC2, as all sites displayed substantial overlap in PC3 and PC4 except the results between the Yoshigo and Kuma-Nishioda sites (see ).


The statistical results of the first comparison indicated that statistically significant results were found only in PC2 and PC3. PC2 focused on the interphase and interregion differences among all populations, and PC2 and PC3 captured the difference between the Jōmon and Yayoi periods (). showed that the Jōmon populations substantially overlapped, while the Yayoi populations were more scattered spatiotemporally in PC2 and PC3, which was supported by the fact that the effect sizes of statistical tests between the Yayoi populations (Z/n=.113 in spatial and temporal comparisons in PC2) were higher than the ones between the Jōmon populations (Z/n<.055 in spatial comparisons and <.029 in temporal comparisons in PC2; ).
showed that the most pronounced differences of the second comparison were observed in PC1, with one notable exception seen when comparing the Ebishima and Kuma-Nishioda sites, where the differences were more conspicuous in PC2. Overall, the Jōmon individual populations were more varied than the Yayoi population. shows a visual representation of the configurational changes in landmarks for each comparison. The figure was constructed using the plotRefToTarget function of the geomorph package in R. Straight bars indicate the degree of variation in each landmark. These overall differences are primarily related to facial height and tooth position, specifically highlighting that individuals from the Kuma-Nishioda site tend to exhibit taller facial features and a more anterior tooth placement.

Discussion
Principal Findings and Comparison With Previous Research
The results of the PCA presented herein were not straightforward. PC2 showed statistically significant spatiotemporal differences between some regions and phases, which was consistent with previous biometric studies claiming interphase variation [-] and geographical clines from north to south [-,], although it is difficult to conclude that their differences are pronounced. This is because PC1 did not show any statistical differences, and a similar study of the Kofun period, examining a larger set of 3D data of human crania using geometric morphometrics, exhibited a geological cline in PC1 []. Moreover, as mentioned in the Introduction, given the archeological evidence for the spatiotemporal distinctiveness of the Jōmon material cultures, it would be expected to find statistically significant differences in PC1 as well. Thus, it is possible that the spatiotemporal differences throughout the Jōmon period were more nuanced or relatively small.
Comparisons between the Jōmon and Yayoi populations are consistent with the above interpretations. It is suggested that morphological variations among Jōmon populations are relatively less varied spatiotemporally than in the Yayoi populations, and individual populations are more diverse than the Yayoi population from the Kuma-Nishioda site. It has been widely proposed that Yayoi populations, particularly in the Kyushu region, are descended from migratory groups originating from the Korean peninsula [,,]. Suppose that a genetic bottleneck resulting from this migration contributed to the reduced variation observed in the Yayoi people of northern Kyushu. It suggests that the morphological and genetic diversity among the Jōmon populations were not relatively limited.
Examining population interactions in the Jōmon period could provide important insights into their resilience to environmental fluctuations at the time. As stated in the Introduction, the Jōmon period experienced more severe climate changes than subsequent periods such as the Yayoi and Kofun periods, which sometimes showed a radical societal change. The reason why less evidence of societal changes—such as the occurrence of warfare and emergence of social hierarchies—was found in the Jōmon period may be because the Jōmon populations interacted widely and continuously and frequently exchanged knowledge, skills, and resources against environmental disruptions.
Limitations and Future Directions
First, even if our results suggest that Jōmon people moved or interacted widely, it is still possible that such movements or interactions were not continuous but rather discontinuous with some drastic environmental changes such as the Kikai-Akahoya eruption and the 4.2 ka event [-]. The frequency and degree of interactions naturally depended on the region and phase, as suggested by research on pottery distributions []. It should be noted, however, that certain regional interactions must be maintained to sustain homogeneity across different regions. How frequent and continuous interactions could maintain the morphological variations revealed in this study should be explored in future work.
Second, our discussions have focused on PCs with a contribution rate of higher than 5 percentage points, which covers less than 40% of total cumulative proportions. Although the results of PCA on 3D data typically tend to be dispersed across more PCs [,,], we should be careful not to overestimate the present results. The question of how frequently the Jōmon people interacted should be explored from additional aspects, including quantitative investigations of relevant archeological remains and mathematical simulations or modeling of the Jōmon population.
Third, some results from genetic or isotopic research on Jōmon skeletal remains, suggesting that Jōmon people interacted with populations from different regions, are consistent with our interpretation [,]. Cooke et al [] argued that Jōmon people on the mainland of the Japanese archipelago were isolated from other islands, which might have contributed to the morphological homogeneity of Jōmon people. Nevertheless, to strengthen or verify this hypothesis, further genetic data and morphological data from other relevant regions overseas should be collected and combined [,,]. For instance, our sample size is much larger than that in a study by Buck et al [], which claimed that dietary factors influenced the shape of the Jōmon neurocranium, while our results should also be combined with isotopic dietary data in future work.
Fourth, while the dataset comprising larger sample sizes is acknowledged, potential sampling bias may exist, particularly due to restrictions on obtaining samples, notably in the Kyushu region. However, previous biometric studies have suggested that individuals from the northern Kyushu region do not exhibit significant differences compared to other regions []. This contributes to the overall reliability of the findings presented in this study. Our study exclusively focuses on cranial 3D data, although previous investigations have examined other parts of skeletal remains, such as limbs [,], warranting a comprehensive examination of their 3D data.
Finally, the selection of landmarks in this study is a critical aspect that warrants consideration. Although the chosen landmarks were selected for their clarity and suitability, the effectiveness of the landmarks in capturing morphologically significant changes may vary depending on various factors such as geographic regions, time periods, age groups, and other contextual factors.
Acknowledgments
We thank for the following cultural property or archeological centers, institutions, museums, and universities (see for detailed information): Iwate Prefectural Museum, Osaki City Board of Education, Ofunato City Museum, Kamikuroiwa Archaeological Museum, Shiura Historical Museum, Tohoku History Museum, Tome City History Museum, Niigata University of Health and Welfare, Yokohama City Archaeological Center, Hamamatsu City Museum, Matsudo City Board of Education, Anjo City Museum of History, Okayama University of Science, Kasaoka City Board of Education, Kyoto University, Takayama Village History Museum, National Museum of Nature and Science, Sakaki Town Archaeological Center, Sakiyama shell midden Museum, Shiga Prefectural Archaeological Center, Kagoshima Prefecture Archaeological Center, Toride City Archaeological Center, Shirokawa History Museum, Aomori Prefectural Museum, Chiba Prefectural Board of Education, Osaka Metropolitan University, Tahara Municipal Museum, the Tohoku University Museum, Minamichita Town Board of Education, Toyama Prefecture Archaeological Center, Hirajo Public Hall, Nagoya City Archaeological Center, Nagoya City Museum, Kariya City History Museum, Kasukabe City History Museum, Mizuko shell midden Museum, Saitama City Urawa Museum, Kaizu History Museum, Tobinodai Historic Site Park Museum, Iwata City Board of Education, and Kazuhiko Tanaka (Nagano Nishi High School). We also appreciate the following researchers for their highly valuable suggestions, comments, information, and cooperation: Kazuhiro Sakaue (National Museum of Science), Atsushi Fujisawa (The Tohoku University Museum), Takafumi Nara (Niigata University of Medicine and Welfare), Mikiko Abe, Masatake Kai (Osaka Metropolitan University), Naoto Tomioka (Okayama University of Science), Isao Yumura (Aoya Kamijichi Historical Park), and Ai Takeuchi (Nanzan University).
Funding
The present research was supported by JSPS (No. 24H02201 and 20K00019) and Nanzan University Pache Research Subsidy II-B for the 2025 academic year.
Authors' Contributions
All authors designed the research. HN, TN, MY, and KT gathered the 3D data. HN analyzed the data and wrote the original draft. All authors edited the paper and approved the final manuscript.
Conflicts of Interest
None declared.
Metadata of the samples.
XLSX File, 293 KBSummary of the statistical tests.
XLSX File, 38 KBPlots of sexual differences.
PDF File, 7731 KBComparisons of the Jōmon and Yayoi populations.
PDF File, 2661 KBReferences
- Taniguchi Y. Intoroduction of Archaeology of Jōmon Period. Douseisha; 2019. ISBN: 4886217915
- Habu J. Ancient Jōmon of Japan. Cambridge University Press; 2014. ISBN: 0521772133
- Imamura K. Prehistoric Japan: New Perspective on Insular East Asia. Routledge; 1996. ISBN: 1857286170
- Matsugi T. What Are Kofuns? From a Cognitive Archaeology Perspective. Kadokawa; 2011. ISBN: 4044007632
- Mizoguchi K. The Archaeology of Japan: From the Earliest Rice Farming Villages to the Rise of the State. Cambridge University Press; 2013. ISBN: 0521711886
- Burroughs WJ. Climate Change in Prehistory: The End of the Reign of Chaos. Cambridge University Press; 2005. [CrossRef]
- Maher LA, Banning EB, Chazan M. Oasis or mirage? Assessing the role of abrupt climate change in the prehistory of the Southern Levant. CAJ. Feb 2011;21(1):1-30. [CrossRef]
- Munoz SE, Gajewski K, Peros MC. Synchronous environmental and cultural change in the prehistory of the northeastern United States. Proc Natl Acad Sci U S A. Dec 21, 2010;107(51):22008-22013. [CrossRef] [Medline]
- Schug GR, editor. The Routledge Handbook of the Bioarchaeology of Climate and Environmental Change. Routledge; 2020. [CrossRef] ISBN: 0367560380
- Koyama S. The Jōmon Period. Chuo Koronsha; 1984. ISBN: 4121007336
- Kodama D. Komakino stone circle and its significance for the study of Jomon social structure. In: Hongo H, editor. Hunter-Gatherers of the North Pacific Rim (Senri Ethnological Studies 63). National Museum of Ethnology; 2003:235-261. [CrossRef]
- Nakao H, Tamura K, Arimatsu Y, Nakagawa T, Matsumoto N, Matsugi T. Violence in the prehistoric period of Japan: the spatio-temporal pattern of skeletal evidence for violence in the Jomon period. Biol Lett. Mar 2016;12(3):20160028. [CrossRef] [Medline]
- Nakao H, Nakagawa T, Tamura K, Yamaguchi Y, Matsumoto N, Matsugi T. Violence and climate change in the Jōmon period. In: Schug GR, editor. The Routledge Handbook of the Bioarchaeology of Climate and Environmental Change. Routledge; 2020:364-376. [CrossRef]
- Nakagawa T, Nakao H. A critical evaluation of recent violence and warfare arguments with reference to data from human skeletal remains in Mesolithic Europe [Website in Japanese]. Nihon Kokogaku. 2017;44:65-77. URL: https://ndlsearch.ndl.go.jp/books/R000000004-I028677126 [Accessed 2025-11-05]
- Nakagawa T, Nakao H, Tamura K, Arimatsu Y, Matsumoto N, Matsugi T. Violence and warfare in prehistoric Japan. Lett Evol Behav Sci. 2017;8(1):8-11. [CrossRef]
- Nakagawa T, Tamura K, Yamaguchi Y, Matsumoto N, Matsugi T, Nakao H. Population pressure and prehistoric violence in the Yayoi period of Japan. J Archaeol Sci. Aug 2021;132:105420. [CrossRef]
- Frayer DW. Ofnet: evidence for a mesolithic massacre. In: Martin DL, Frayer DW, editors. Troubled Times: Violence and Warfare in the Past. Gordon & Breach; 1997:181-216. ISBN: 9056995340
- Lahr MM, Rivera F, Power RK, et al. Inter-group violence among early Holocene hunter-gatherers of West Turkana, Kenya. Nature New Biol. Jan 21, 2016;529(7586):394-398. [CrossRef]
- Wendorf F, editor. The Prehistory of Nubia. Southern Methodist University Press; 1968. ISBN: 087074125X
- Kosugi Y, Taniguchi Y, Nishida Y, Mizunoe K, Yano K, editors. The Outline of the Jōmon Culture: A Comparative Approach. Doseisha; 2010. ISBN: 9784886215086
- Booy G, Hendriks RJJ, Smulders MJM, Van Groenendael JM, Vosman B. Genetic diversity and the survival of populations. Plant Biol (Stuttg). Jul 2000;2(4):379-395. [CrossRef]
- Hughes AR, Stachowicz JJ. Genetic diversity enhances the resistance of a seagrass ecosystem to disturbance. Proc Natl Acad Sci U S A. Jun 15, 2004;101(24):8998-9002. [CrossRef] [Medline]
- Hughes AR, Inouye BD, Johnson MTJ, Underwood N, Vellend M. Ecological consequences of genetic diversity. Ecol Lett. Jun 2008;11(6):609-623. [CrossRef] [Medline]
- Verdu P, Austerlitz F, Estoup A, et al. Origins and genetic diversity of pygmy hunter-gatherers from Western Central Africa. Curr Biol. Feb 24, 2009;19(4):312-318. [CrossRef] [Medline]
- Walsh SJ, Mitchell RJ, Watson N, Buckleton JS. A comprehensive analysis of microsatellite diversity in Aboriginal Australians. J Hum Genet. 2007;52(9):712-728. [CrossRef] [Medline]
- Wang S, Lewis CM, Jakobsson M, et al. Genetic variation and population structure in native Americans. PLoS Genet. Nov 2007;3(11):e185. [CrossRef] [Medline]
- Harpending H. Regional variation in !Kung populations. In: DeVore I, Lee RB, editors. Kalahari Hunter-Gatherers: Studies of the !Kung San and Their Neighbors. Harvard University Press; 1976:152-165. [CrossRef] ISBN: 1583481257
- MacDonald DH, Hewlett BS. Reproductive interests and forager mobility. Curr Anthropol. Aug 1999;40(4):501-524. [CrossRef]
- Marlowe FW. The Hadza: Hunter-Gatherers of Tanzania. University of California Press; 2010. ISBN: 0520253418
- Morgan C. Climate change, uncertainty and prehistoric hunter–gatherer mobility. J Anthropol Archaeol. Dec 2009;28(4):382-396. [CrossRef]
- Ritchison BT, Thompson VD, Lulewicz I, Tucker B, Turck JA. Climate change, resilience, and the Native American Fisher-hunter-gatherers of the late Holocene on the Georgia coast, USA. Quat Int. May 2021;584:82-92. [CrossRef]
- Martin R. Lehrbuch Der Anthropologie in Systematischer Darstellung [Book in German]. 2nd ed. Gustav Fischer; 1928. [CrossRef]
- Caple J, Stephan CN. A standardized nomenclature for craniofacial and facial anthropometry. Int J Legal Med. May 2016;130(3):863-879. [CrossRef] [Medline]
- Dodo Y. A metrical analysis of Jomon crania from the Tohoku district. The Journal of Anthropological Society of Nippon. 1982;90(Supplement):119-128. [CrossRef]
- Dodo Y. Metrical and non‒metrical analyses of Jōmon crania from eastern Japan. Prehistoric Hunter-Gatherers in Japan. 1986;27:137-161. URL: https://cir.nii.ac.jp/crid/1570291224292584064 [Accessed 2025-11-05]
- Fukase H, Wakebe T, Tsurumoto T, Saiki K, Fujita M, Ishida H. Geographic variation in body form of prehistoric Jomon males in the Japanese archipelago: its ecogeographic implications. American J Phys Anthropol. Sep 2012;149(1):125-135. [CrossRef]
- Hanihara T, Ishida H. Regional differences in craniofacial diversity and the population history of Jomon Japan. American J Phys Anthropol. Jul 2009;139(3):311-322. [CrossRef]
- Kondo O. Regional variation of Jōmon people as viewed from cranial metrics. Bulletin of the National Museum of Japanese History. 2018;208:249-267. [CrossRef]
- Mouri T. Incidences of cranial nonmetric characters in five Jomon populations from west Japan. The Journal of Anthropological Society of Nippon. 1988;96(3):319-337. [CrossRef]
- Ogata T, editor. Japanese I: Anthropology. Vol 5. Yuzankaku shuppann. Japanese; 1981. ISBN: 9784639025092
- Yamaguchi B. A review of the osteological characteristics of the Jomon population in prehistoric Japan. The Journal of Anthropological Society of Nippon. 1982;90(Supplement):77-90. [CrossRef]
- Yamaguchi B. Human skeletal remains in the Jōmon period. In: Kato S, Kobayashi T, Fujimoto T, editors. Research on the Jōmon Culture. Yuzankaku; 1982:15-88. ISBN: 9784639008798
- Abe Y. Population diffusion and return of nomadic hunters in the initial Jōmon period. In: Abe Y, editor. Social History of Migration and Distribution in the Jōmon Period. Yuzankaku; 2010:233-253. ISBN: 9784639021445
- Fukunaga M. East and West in the Jōmon Society: Social Structure of the Late Jōmon Period. Yuzankaku; 2020. ISBN: 9784639026952
- Yano K. The Jōmon Society of the West Japan from the Point of View of Pottery Types. Doseisha; 2016. ISBN: 9784886217240
- Suzuki N. Coexistence of pottery from different lineages in Shomyoji type: process of transmission and change of pottery from different lineages. In: Imamura K, editor. Interactions between Different Types of Pottery. Doseisha; 2011:49-71. ISBN: 9784886215802
- Matsumoto N. Jōmon Community and Society. Iwanami shoten; 2005. ISBN: 9784000067881
- Yamada Y. The History of Jōmon Period. Koudansha; 2019. ISBN: 9784065143681
- García-Medrano P, Maldonado-Garrido E, Ashton N, Ollé A. Objectifying processes: the use of geometric morphometrics and multivariate analyses on Acheulean tools. JLS. 2020;7(1). [CrossRef]
- Hashemi SM, Vahdati Nasab H, Berillon G, Oryat M. An investigation of the flake-based lithic tool morphology using 3D geometric morphometrics: a case study from the Mirak Paleolithic Site, Iran. Journal of Archaeological Science: Reports. Jun 2021;37:102948. [CrossRef]
- Shott MJ, Trail BW. Exploring new approaches to lithic analysis: laser scanning and geometric morphometrics. Lithic Technology. Sep 2010;35(2):195-220. [CrossRef]
- Buck LT, De Groote I, Hamada Y, Hassett BR, Ito T, Stock JT. Evidence of different climatic adaptation strategies in humans and non-human primates. Sci Rep. Jul 30, 2019;9(1):11025. [CrossRef] [Medline]
- Buck LT, Menéndez LP, De Groote I, Hassett BR, Matsumura H, Stock JT. Factors influencing cranial variation between prehistoric Japanese forager populations. Archaeol Anthropol Sci. 2024;16(1):3. [CrossRef] [Medline]
- Makishima H, Ogihara N. Three-dimensional geometric morphometric study of craniofacial variations in Jomon populations. AS (J Ser). 2009;117(1):11-21. [CrossRef]
- Nakao H, Nakagawa T, Kaneda A, Tamura K, Noshita K. Demic diffusion of the Yayoi people in the Japanese archipelago. Lett Evol Behav Sci. 2023;14(2):58. [CrossRef]
- Nakao H, Kaneda A, Tamura K, Noshita K, Nakagawa T. Macro-scale population patterns in the Kofun period of the Japanese archipelago: quantitative analysis of a larger sample of three-dimensional data from ancient human crania. Humans. 2024;4(2):131-147. [CrossRef]
- Kuzminsky SC, Reyes Báez O, Arriaza B, et al. Investigating cranial morphological variation of early human skeletal remains from Chile: a 3D geometric morphometric approach. American J Phys Anthropol. Feb 2018;165(2):223-237. [CrossRef]
- Zelditch ML. Geometric Morphometrics for Biologists: A Primer. Academic Press; 2012. ISBN: 9780123869036
- Kaneda A, Nakagawa T, Tamura K, Noshita K, Nakao H. A proposal of a new automated method for SfM/MVS 3D reconstruction through comparisons of 3D data by SfM/MVS and handheld laser scanners. PLOS One. 2022;17(7):e0270660. [CrossRef] [Medline]
- Nakagawa T, Kaneda A, Tamura K, Nakao H. A comparative study of measurement methods digitizing human skeletal remains: SfM and laser scanning. Nabunken Ronso. 2022;3:39-64. URL: https://www.nabunken.go.jp/english/paper/ronso-no3.html [Accessed 2025-11-05]
- Nakao H, Nakagawa T, Yoshida M. 3D data of human skeletal remains acquired by two kinds of laser scanners: Einscan Pro HD and Creafrom HandySCAN BLACK. Journal of Nanzan Academic Society Humanities and Natural Sciences. 2022;24:309-314. URL: https://nanzan-u.repo.nii.ac.jp/records/4127 [Accessed 2025-11-05]
- Naito A. Human skeletal remains in the Yayoi period. In: Ogata T, editor. Japanese I: Anthropology. 1981:57-99. ISBN: 9784639000778
- Nakahashi T. Physical traits of the Yayoi people. In: Nagai M, Nasu T, Kanaseki H, Sahara M, editors. Research on the Yayoi Culture. 1989:23-51. ISBN: 9784639008590
- Nakagawa T. Violence in the Aoyakamijichi site. Mater Cult; 2021:105-121. URL: https://ndlsearch.ndl.go.jp/books/R000000004-I031586567 [Accessed 2025-11-05]
- Cignoni P, Callieri M, Corsini M, Dellepiane M, Ganovelli F, Ranzuglia G. MeshLab: an open-source mesh processing tool. 2008. Presented at: Sixth Eurographics Italian Chapter Conference; Jul 2-4, 2008:129-136; Salerno, Italy. [CrossRef]
- Noshita K, Kaneda A, Tamura K, Nakagawa T, Nakao H. A mathematical study on outline morphologies of Yayoi pottery: a case study focusing on the Ongagawa-style pottery from the Yamura, Yano, and Ayaragi-go sites. Nabunken Ronso. 2022;3:65-82. URL: https://www.nabunken.go.jp/english/paper/ronso-no3.html [Accessed 2025-11-05]
- Noshita K, Kaneda A, Tamura K, Nakagawa T, Nakao H. A comparative analysis of two- and three-dimensional data using Ongagawa pottery as an example. J Comput Archaeolo. 2022;27(1/2):1-10. URL: https://ndlsearch.ndl.go.jp/books/R000000004-I032174522 [Accessed 2025-11-05]
- Fantini M, de Crescenzio F, Persiani F, Benazzi S, Gruppioni G. 3D restitution, restoration and prototyping of a medieval damaged skull. Rapid Prototyp J. Sep 26, 2008;14(5):318-324. [CrossRef]
- Nakagawa T, Yoshida M, Nakao H. A geometric morphometric analysis of human craniums in the Kofun period from the Okayama, Hiroshima, and Hyogo prefecture. Kodai Kibi; 2022:43-60. URL: https://ndlsearch.ndl.go.jp/books/R000000004-I032847725 [Accessed 2025-11-05]
- Adams DC, Collyer ML. Geomorph: Software for geometric morphometric analyses R package version 404. 2022. URL: https://cran.r-project.org/package=geomorph [Accessed 2025-10-22]
- R Core Team. The R Project for Statistical Computing. 2020. URL: https://www.R-project.org [Accessed 2025-10-22]
- RStudio: Integrated Development for R RStudio, PBC. 2020. URL: https://posit.co/download/rstudio-desktop/ [Accessed 2025-10-22]
- Buikstra JE, Ubelaker D. Standards for data collection from human skeletal remains. Proceedings of a seminar at the Field Museum of Natural History. Internet Archive. 1994. URL: https://archive.org/details/isbn_1563490757_t0j9 [Accessed 2025-11-05]
- Brooks S, Suchey JM. Skeletal age determination based on the os pubis: a comparison of the Acsádi-Nemeskéri and Suchey-Brooks methods. Hum Evol. Jun 1990;5(3):227-238. [CrossRef]
- Meindl RS, Lovejoy CO. Ectocranial suture closure: a revised method for the determination of skeletal age at death based on the lateral‐anterior sutures. American J Phys Anthropol. Sep 1985;68(1):57-66. [CrossRef]
- Sakaue K. A Bayesian approach to age estimation from cranial suture closure in Japanese people. Bulletin of National Museum of Natural Science Series; 2015:1-11. URL: https://www.kahaku.go.jp/english/research/publication/anthropology/v41.html [Accessed 2025-11-05]
- Cultural properties protection act. E-gov. 1950. URL: https://laws.e-gov.go.jp/law/325AC0100000214 [Accessed 2025-10-22]
- Kajita H, Isaji Y, Kato R, et al. Climatic change around the 4.2 ka event in coastal areas of the East China Sea and its potential influence on prehistoric Japanese people. Palaeogeogr Palaeoclimatol Palaeoecol. Jan 2023;609:111310. [CrossRef]
- Kawahata H. Climatic reconstruction at the Sannai-Maruyama site between Bond events 4 and 3—implication for the collapse of the society at 4.2 ka event. Prog Earth Planet Sci. Dec 2019;6(1):63. [CrossRef]
- Maeno F, Taniguchi H. Spatiotemporal evolution of a marine caldera-forming eruption, generating a low-aspect ratio pyroclastic flow, 7.3 ka, Kikai caldera, Japan: implication from near-vent eruptive deposits. Journal of Volcanology and Geothermal Research. Nov 2007;167(1-4):212-238. [CrossRef]
- Noshita K, Nakagawa T, Kaneda A, Tamura K, Nakao H. The cultural transmission of Ongagawa style pottery in the prehistoric Japan: quantitative analysis on three-dimensional data of archaeological pottery in the early Yayoi period. J R Soc Interface. Feb 2025;22(223):20240889. [CrossRef] [Medline]
- Kawanishi K. A preliminary report on the odake shell midden project in the reiwa 4. Toyama Prefecture Cultural Property Center; Feb 19, 2023.
- Kusaka S, Saito T, Ishimaru E, Yamada Y. Strontium isotope analysis on human skeletal remains from the Hobi and Ikawazu shell-mounds in Aichi Prefecture, Japan. AS. 2022;130(1):25-32. [CrossRef]
- Cooke NP, Mattiangeli V, Cassidy LM, et al. Ancient genomics reveals tripartite origins of Japanese populations. Sci Adv. Sep 17, 2021;7(38):eabh2419. [CrossRef] [Medline]
- Adachi N, Kanzawa-kiriyama H, Nara T, Kakuda T, Nishida I, Shinoda KI. Ancient genomes from the initial Jomon period: new insights into the genetic history of the Japanese archipelago. AS. 2021;129(1):13-22. [CrossRef]
- Wang T, Wang W, Xie G, et al. Human population history at the crossroads of East and Southeast Asia since 11,000 years ago. Cell. Jul 8, 2021;184(14):3829-3841. [CrossRef] [Medline]
Abbreviations
| PC: principal component |
| PCA: principal component analysis |
Edited by Amy Schwartz; submitted 10.Feb.2025; peer-reviewed by Denise Crampton, Osamu Kondo; final revised version received 08.Aug.2025; accepted 24.Sep.2025; published 12.Nov.2025.
Copyright© Hisashi Nakao, Akihiro Kaneda, Kohei Tamura, Koji Noshita, Mayu Yoshida, Tomomi Nakagawa. Originally published in JMIRx Bio (https://bio.jmirx.org), 12.Nov.2025.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work, first published in JMIRx Bio, is properly cited. The complete bibliographic information, a link to the original publication on https://bio.jmirx.org/, as well as this copyright and license information must be included.

