# Bone loss markers in the earliest Pacific Islanders

Justyna Miszkiewicz, Frédérique Valentin, Christina Vrahnas, Natalie Sims, Jitraporn Vongsvivut, Mark Tobin, Geoffrey Clark

<figure><img src="image_1.png" /></figure>

## To cite this version:

Justyna Miszkiewicz, Frédérique Valentin, Christina Vrahnas, Natalie Sims, Jitraporn Vongsvivut, et al. Bone loss markers in the earliest Pacific Islanders. Scientific Reports, 2021, 11 (1), 10.1038/s41598-021-83264-3. [DOI] [arXiv]

HAL Id: hal-03322149

https://hal.science/hal-03322149

Submitted on 16 Dec 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

OPEN

# Bone loss markers in the earliest Pacific Islanders

Justyna J. Miszkiewicz<fnref n="1" />✉, Frédérique Valentin²,³, Christina Vrahnas⁴,⁵,⁶, Natalie A. Sims⁵,⁶, Jitraporn Vongsvivut⁷, Mark J. Tobin⁷ & Geoffrey Clark³

Kingdom of Tonga in Polynesia is one of the most obese nations where metabolic conditions, sedentary lifestyles, and poor quality diet are widespread. These factors can lead to poor musculoskeletal health. However, whether metabolic abnormalities such as osteoporosis occurred in archaeological populations of Tonga is unknown. We employed a microscopic investigation of femur samples to establish whether bone loss afflicted humans in this Pacific region approximately 3000 years ago. Histology, laser confocal microscopy, and synchrotron Fourier-transform infrared microspectroscopy were used to measure bone vascular canal densities, bone porosity, and carbonate and phosphate content of bone composition in eight samples extracted from adult Talasiu males and females dated to 2650 BP. Compared to males, samples from females had fewer vascular canals, lower carbonate and phosphate content, and higher bone porosity. Although both sexes showed evidence of trabecularised cortical bone, it was more widespread in females (35.5%) than males (15.8%). Our data suggest experiences of advanced bone resorption, possibly as a result of osteoporosis. This provides first evidence for microscopic bone loss in a sample of archaeological humans from a Pacific population widely afflicted by metabolic conditions today.

The global obesity pandemic is increasing at an alarming rate posing growing health concerns and an economic burden for human generations today and in the future<fnref n="1" />. In 2019, the cost of treating obesity was estimated at USD $2 trillion (~2.8% of gross domestic product globally), with 9% of world population considered obese, and a further 392 million of youth predicted to become overweight or obese by 2025¹,². Human populations in the Pacific Islands are consistently classified as some of the world's most obese nations³. The Kingdom of Tonga in Polynesia, in particular, has one of the largest prevalence in type 2 diabetes⁴. The reasons behind such high prevalence of metabolic conditions across the Pacific are complex. Studies have identified hereditary factors³,⁵ and environmental constraints that interplay with genetics, individual lifestyle, and access to resources⁶,⁷. While modern sedentary lifestyle is understood to underlie many metabolic conditions⁸, shedding light on the roots of these conditions can be predominantly achieved by examining health and disease markers left on human remains in the archaeological record⁹,¹⁰. Prior research examining the surviving skeletons of archaeological Pacific Islanders reported possible evidence of diseases such as gout¹¹ and Forestier's disease (diffuse idiopathic skeletal hyperostosis—DISH)¹² in ca. 3,000-year-old Vanuatu. Both diseases occur in the same populations today¹³, which suggests that the origins of these cannot be completely modern. Therefore, the study of archaeological human skeletal remains has the potential to fundamentally shift our understanding of the origins, or at least continuity in occurrence, of modern health problems¹⁴,¹⁵.

Obesity is typically associated with physical inactivity and nutrition of poor quality¹⁶. Both factors can have detrimental effects on bone physiology and structure because of a minimised mechanical loading regime¹⁷ and restricted dietary calcium and vitamin consumption¹⁸,¹⁹, which can lead to an increase in bone turnover where bone loss dominates bone gain¹⁸. This advances bone tissue age and can result in greater fragility, susceptibility to fractures, and a diagnosis of osteopenia or osteoporosis¹⁹. The deteriorating bone microarchitecture is caused by a shift in the activity balance between cells that normally renew the skeleton by bone remodelling such that osteoclast-mediated bone resorption exceeds bone deposition by osteoblasts²⁰,²¹. Continuing bone loss in such

<footnote n="1">¹School of Archaeology and Anthropology, Australian National University, 44 Linnaeus Way, Canberra, ACT 2601, Australia. ²CNRS, UMR 7041, ArScAn, Ethnologie préhistorique, Maison René-Ginouvès, Archéologie et Ethnologie, 21 Allée de l'Université, 92023 Nanterre Cedex, France. ³Archaeology and Natural History, School of Culture History and Language, College of Asia and the Pacific, Australian National University, Canberra, ACT 2601, Australia. ⁴Bone Biology and Disease Unit, St. Vincent's Institute of Medical Research, 9 Princes Street, Fitzroy, Melbourne, VIC 3065, Australia. ⁵Department of Medicine, St. Vincent's Hospital, The University of Melbourne, Melbourne, VIC 3065, Australia. ⁶MRC Protein Phosphorylation and Ubiquitylation Unit, James Black Centre, University of Dundee, Dundee DD1 5EH, UK. ⁷Infrared Microspectroscopy Beamline, ANSTO - Australian Synchrotron, 800 Blackburn Road, Clayton, VIC 3168, Australia. ⁸email: Justyna.Miszkiewicz@anu.edu.au</footnote>
fragile bones of individuals, particularly the elderly, with untreated osteoporosis impairs lifestyle and mobility as fractures occur repeatedly or heal poorly. Osteoporosis is common in modern aging populations, and as such poses a large economic burden. The condition is known as the “silent disease” manifesting only once fracture has occurred²²⁻²⁴. The difficulty in completely preventing osteoporosis lies in its multi-factorial and inter-linked aetiology. Osteoporosis, meaning ‘porous bones’, can take different forms (e.g. primary or secondary) and as such develops as a result of various factors including ageing²⁵, female biological sex²⁶, lack of physical activity²⁷, genetic predisposition²⁸, hormone deficiencies²⁹, dietary restrictions¹⁸,¹⁹, including those resulting from psychiatric conditions such as anorexia nervosa³⁰, and a combination of these and other variables that can relate to socio-economic status¹⁴,¹⁵. The effect of sex on bone health has been particularly well studied in post-menopausal women, where loss of oestrogen, a hormone essential in inhibiting bone resorption, is a causative factor³¹,³².

There has been much discussion about the relationship between obesity and osteoporosis³³⁻³⁶. Not only does a lack of exercise negatively impact bone density, a cell based link between fat and bone tissue has been proposed³⁷. As bone loss increases with age, the number of bone marrow fat cells (adipocytes) rises³⁷. Bone-forming osteoblasts and fat cells differentiate from the same bone-marrow precursor cells, yet whether bone loss is a precursor to fat infiltration or vice versa remains a point of contention³⁷. To complicate matters further, a heavier frame can have a positive effect on bone formation due to the higher weight-bearing from increased body mass—a relationship which has been termed the “obesity paradox”³⁵,³⁶. The development of osteoporotic bone in obese patients has been demonstrated in clinical studies, though data have not been consistent³⁸⁻⁴⁰. For example, some sites of the skeleton, such as vertebrae, but not the hip, have been classified as osteoporotic in obese or overweight individuals³⁹. Obesity incidence has been reported in Tonga⁴, but osteoporosis data are extremely limited. Reports from the past two decades paint a conflicting picture. In 2006, Tonga (included amongst other countries from the “Western Pacific”) was ranked second highest, after Europe, for prevalence of osteoporotic fractures⁴¹. In contrast, “the Polynesian bone phenotype” has been characterised as having higher bone density than that of European counterparts⁴². The International Osteoporosis Foundation is yet to map fracture incidence across the Pacific Islands comprehensively⁴³. The efficiency of data collection there is possibly hindered by limited radiography infrastructure⁴².

Clinically, osteoporosis in modern populations is diagnosed from measures of bone mineral density (BMD) and occurrence of fractures⁴⁴⁻⁴⁶. Whether an individual suffers from abnormal bone fragility defined as osteopenia or osteoporosis is diagnosed from BMD T-scores that are lower than those of healthy age-matched cohorts⁴⁷. A gold standard methodology to collect BMD in living people is dual energy x-ray absorptiometry (bone densitometry, DXA)⁴⁸. Diagnosis can be supplemented through bone turnover markers from bone histomorphometry⁴⁹, or quantitative ultrasound techniques⁵⁰, and various other types of microscopy that measure bone porosity two- (2D) or three- (3D) dimensionally⁵¹.

Similar means have been implemented to confirm the incidence of possible osteoporosis in various archaeological European⁵²⁻⁵⁷, North American⁵⁸, and African⁵⁹⁻⁶¹ osteological collections (see reviews¹⁵,⁶²⁻⁶⁴). Where access to large and well-preserved specimens has been possible, DXA, cortical and/or trabecular bone volume, and bone turnover assessment methods have indicated significant bone loss⁵²⁻⁵⁵,⁶⁰,⁶⁵⁻⁶⁷. Analyses of fracture types and their patterns (such as classic Colles' fractures) have also provided insights into archaeological experiences of osteoporosis⁶⁶⁻⁶⁹. Microscopic methods reconstructing bone turnover, mineral composition, and trabecular and cortical bone microarchitecture in archaeological human remains have also been successful⁷⁰⁻⁷². Non-invasive approaches such as 3D micro-computed tomography (micro-CT) can help in the calculation of trabecular or cortical bone thickness, separation, and the visualisation of porosity network⁷³,⁷⁴. Invasive methods such as histology rely on an extraction of samples to describe more localised bone changes⁷²,⁷⁵. Therefore, invasive techniques are less often applied to archaeological samples because of the irreplaceable nature of excavated remains⁷⁶, and so typically rely on smaller sample sizes⁷⁷. Additionally, preservation issues (such as degradation of bone microstructure post-mortem⁷⁸) can further limit the application of modern clinical osteoporosis diagnostic criteria to archaeological remains⁷¹,⁷⁹. Therefore, while the examination of disease in archaeological remains can provide fruitful results, response of living bone to disease cannot be observed. However, a combination of different methods can provide multiple lines of data for differential diagnosis and interpretation.

To the best of our knowledge, bone loss as a result of possible osteoporosis has not yet been investigated in archaeological samples in the Pacific. Obesity specifically as a factor leading to osteoporosis cannot be directly reconstructed for archaeological specimens as the available methods of body mass estimation still provide large error⁸⁰. However, given the widespread health issues in contemporary Pacific nations, bone loss markers in archaeological human remains in this region could be an indicator of musculoskeletal health and metabolic issues similar to those seen there today. The 2,650 BP archaeological site of Talasiu in the Kingdom of Tonga is one of the most significant sites ever uncovered in the Pacific⁸¹,⁸². This site has previously produced the first evidence for Early Polynesian mortuary behaviour⁸¹,⁸², archaeological Polynesian diet relying on marine resources⁸³, and insights into the abandonment of ceramics by early Polynesian people⁸⁴. We set out to test whether archaeological Tongan human bone would show evidence of abnormal bone loss marked by sex-specific differences. We employed an invasive microscopic examination to characterise bone metabolic fluctuations in eight samples extracted from femora in the Talasiu individuals, sub-divided into groups of estimated sex (four males, four females). All individuals were classified as 'adults' who likely lived beyond 30 years of age, which means their bone tissue had entered an age-driven degeneration phase following a peak bone mass accrual point of lifespan. One female (ID: Sk3.1) was possibly an elderly individual (Table S1), who might have survived into a sixth life decade and possibly beyond. We used several microscopic methods to provide multiple lines of data. We recorded sexually dimorphic femur size differences, and then combined 2D histology⁷⁵,⁸⁵, 3D laser scanning confocal microscopy⁵¹, and synchrotron particle accelerator methods (Fourier-transform infrared microspectroscopy, sFTIRM)⁸⁶,⁸⁷ (Fig. 1). This allowed us to examine compact femur bone Haversian canal densities as a proxy for remodelling (Fig. 2), abnormal porosity as a marker of prolonged bone resorption (Figs. 2, 3), and bone matrix
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Summary of methods and key findings in the present study. Posterior view of two right archaeological femora from Talasiu individuals estimated as female (♀, ID: Sk3.1) and male (♂, ID: Sk3.2) shows the sectioning location (black dashed line) and approximate sample removed (red dashed box). Methodological steps included synchrotron sourced infrared microspectroscopy to measure bone mineral composition, histomorphometry to estimate Haversian canal densities (red dots), and laser confocal scanning of bone topography to provide a qualitative illustration of intra-cortical bone porosity producing trabecularisation.</figcaption></figure>

composition through carbonate and phosphate content (Fig. 1), averages of which we hypothesised to differ between the estimated female and male sex groups.

## Results

We found evidence for bone loss, which differed between the estimated sex groups in our samples. Porosity which manifested as cortical bone trabecularisation was particularly advanced in the females. As expected, due to sexual dimorphism, males in our sample had larger femora when compared to the females (Table 1, Fig. 1). The average midshaft circumference was 92.25 mm (standard deviation (SD) = 4.79) in males, but 79.75 mm (SD = 7.09) in females (U = 15.5; *p* = 0.029). The diameter of femur midshaft was also larger and more pronounced in males in both the antero-posterior (A-P), and medio-lateral (M-L) planes. The A-P diameter was 30.75 mm (SD = 2.87) in males, but 27.26 mm (SD = 2.15) in females, albeit this difference was not statistically significant.

The sexual dimorphism was likely also reflected in the 2D histological measurements of bone Haversian canal density collected from the cortical area of bone unaffected by porosity producing trabecularisation, which was higher in males (17.445/mm²) than in females (15.603/mm²) (Table 1, S2; Figs. 1, 2, 4). However, the area of bone samples showing porosity producing trabecularisation which represented as 'giant' coalescing pores was more evident in females than in males (Fig. 1). The average porosity producing trabecularisation was 35.5% (SD = 16.40) in females, but only 15.8% (SD = 24.17) in males. In fact, all four female samples were affected by a high degree of porosity producing trabecularisation that was widespread from the most inner bone layer (endo-cortical) through the mid-section region (intra-cortical), reaching the bone areas immediate to the most outer bone layer (periosteum) (Figs. 1, 2). This led to cortical wall thinning⁸⁸,⁸⁹ and was a consistent pattern throughout the sample as supported by a statistically significant negative correlation between percent areas of trabecularisation and cortical wall area ($Rho = -0.762, p = 0.028$; Fig. 5). However, comparing only the porosity and Haversian canal variables between the two sex groups did not result in statistical significance as males were clearly impacted by intra-cortical trabecularisation as well. In fact, one of the male samples (ID: Sk9.3, Figure S1) showed an estimated 51.8% of trabecularisation of intra-cortical bone. This data point is maximum across our entire sample of eight bones. The remaining three male samples showed an overall consistent appearance of porosity whereby the bone surfaces were dense, showing a limited range of variably sized porosity regions and their occurrence (one complete absence, Figs. 1, 3). On the contrary, the female samples exhibited much more
<figure><img src="image_5.png" /><figcaption>**Figure 2.** A summary of histomorphometric procedures for estimating densities of Haversian canals and abnormal intra-cortical bone porosity. (a, b): localised histology (60×) viewed under linearly polarised light showing secondary osteon lamellae (white arrows) surrounding Haversian canals (oval features), which were obscured by patches of diagenesis (grey arrows). Cement lines (black arrow in b) were inconsistently preserved (examples shown are from individuals BG2 (a) and Sk9.3 (b)). Haversian canals were counted manually from thresholded images (c–f): (c) shows numbered canals, (d) shows manual segmentation of cortical wall from abnormal intra-cortical porosity (red dashed line), while the solid line marks whole section area in individual BG2; (e) illustrates final counts of Haversian canals, while (f) shows 'extracted' abnormal bone porosity.</figcaption></figure>

variability in porosity producing trabecularisation (Table 1, S2; Figs. 1, 2), ranging from 13.9% minimum to 50.2% maximum of cortical bone space. The 3D confocal laser scans of bone surfaces from two cases showing extremes in bone micro-organisation (one female Sk3.1, one male Sk3.2) (Figs. 1, 3, S10, S11) demonstrated an almost entirely limited topography in the male sample, but more variable topography in the female sample which was punctuated with relatively deep and 'giant pores' in the bone cortex.

The bone carbonate and phosphate content measured through sFTIRM was consistently lower in females (Figs. 4, 5; Table 2, S3–S6) ($p$ < 0.05 in all sFTIRM variables except for integrated phosphate area). The average peak height of phosphate (from 848 spectra; $U$ = 1,118,955.000; $p$ < 0.013), and $v_2$ (from 1676 spectra; $U$ = 3,767,893.000; $p$ < 0.0001) and $v_3$ carbonates (from 2159 spectra; $U$ = 4,978,269.000; $p$ < 0.0001) were lower in females. While the phosphate difference of 2.6% in the means seems marginal, it was more substantial when considering $v_2$ (30.3%) and $v_3$ (32%) carbonates. We note that more spectra were successfully measured in the samples from males (2496 phosphate spectra, minimum 3329 $v_2$ and 3264 $v_3$ spectra) (Table 2, S3–S5). Compared to the male samples, the average integrated areas of phosphate and two carbonate peaks were also lower in females. The phosphate area in females was only lower by 1.1% (Table S3). However, the $v_2$ ($U$ = 3,765,291.500;
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Intra-cortical surfaces scanned using a three-dimensional laser confocal topography microscope Olympus OLS5000 with applied false-coloured heat map to illustrate extreme differences in porosity in archaeological posterior midshaft femur samples from Talasiu individuals estimated as female (♀, ID: Sk3.1) and male (♂, ID: Sk3.2). (a) shows topography viewed superiorly from the bone surface, whereas (b) is reverted upside down to show the inferior view of the same scan as (a).</figcaption></figure>

<table><caption>**Table 1.** Sample details and raw gross morphometric data of the femora representing adult archaeological males and females from Talasiu, Tongatapu, Kingdom of Tonga. Cortical width and midshaft circumference are reported in mm. A-P: antero-posterior, M-L: medio-lateral, H.Dn: Haversian canal density per mm², %Po.Ar: percent of sample impacted by porosity producing trabecularisation. *left femur, **probable sex.</caption><thead><tr><th rowspan="2">Individual ID</th><th rowspan="2">Estimated sex</th><th colspan="2">Cortical width</th><th rowspan="2">Midshaft circumference</th><th colspan="2">Quantitative histology</th></tr><tr><th>A-P</th><th>M-L</th><th>H.Dn</th><th>%Po.Ar</th></tr></thead><tbody><tr><td>BG2</td><td>Female</td><td>25.28</td><td>23.35</td><td>76</td><td>22.223</td><td>31.890</td></tr><tr><td>Sk3.1</td><td>Female</td><td>28.93</td><td>24.87</td><td>86</td><td>13.159</td><td>50.232</td></tr><tr><td>Sk3.2</td><td>Male</td><td>33.87</td><td>29.27</td><td>101</td><td>15.627</td><td>7.261</td></tr><tr><td>Sk5</td><td>Male</td><td>27.72</td><td>22.17</td><td>86</td><td>20.798</td><td>0.000</td></tr><tr><td>Sk9.2*</td><td>Female**</td><td>29.30</td><td>26.54</td><td>81</td><td>14.815</td><td>46.001</td></tr><tr><td>Sk9.3</td><td>Male**</td><td>29.00</td><td>27.34</td><td>87</td><td>16.530</td><td>51.793</td></tr><tr><td>Sk12</td><td>Female**</td><td>25.54</td><td>23.26</td><td>76</td><td>12.216</td><td>13.908</td></tr><tr><td>Sk14</td><td>Male</td><td>32.41</td><td>26.92</td><td>95</td><td>16.826</td><td>4.197</td></tr></tbody></table>

$p < 0.0001$); and $v_3$ ($U = 4,847,614.000; p < 0.0001$) carbonates were 35.3% and 44.12% greater in the samples from males (Table S4, S5). The carbonate $v_2$:phosphate integrated area ratios from the female samples were lower by 5% when compared to males ($U = 1,421,109.000; p < 0.0001$; Table S6). This was despite females recording a maximum ratio spectrum of 0.143 contrasted with a much lower maximum spectrum of 0.038 in males, and the higher total number of spectra collected in males (2495) contrasted with 846 spectra in females.

Discussion

We report the first microscopic record of bone loss characteristics in archaeological femora representing some of the earliest Pacific Islanders. On the basis of how commonplace contemporary metabolic issues, including obesity and diabetes, are in Tonga, we interpret our evidence to indicate that similar problems that might have led to osteoporosis occurred there ca. 3000 years ago. Two key implications are presented. Firstly, our small sample of archaeological humans from the Pacific appears to be afflicted by cortical bone porosity producing trabecularisation despite the notion that osteoporosis and related conditions are modern-day diseases. This confirms that the occurrence of widespread metabolic disease in this region today may have roots in the past. It agrees with bone loss experiences reported for archaeological collections from around the world (see reviews$^{15,62–64}$),
BONE MORPHOLOGY AND HISTOMORPHOMETRY

<figure><img src="image_7.png" /></figure>

PHOSPHATE

<figure><img src="image_8.png" /></figure>

CARBONATE, AND CARBONATE:PHOSPHATE RATIOS

<figure><img src="image_9.png" /></figure>

**Figure 4.** Graphical summary of results illustrating descriptive comparisons of data between archaeological Tongan individuals of estimated male (♂) and female (♀) sex. The boxplots indicate that females present with: (a) smaller femoral midshaft (dm: diameter), with lower Haversian canal density (per mm²), but higher intra-cortical porosity producing trabecularisation (%); (b, c) lower phosphate content; and (d–f) lower carbonate, and carbonate:phosphate ratios. Outliers are marked by asterisks and circles. \*p < 0.05; \*\*p < 0.0001.

and archaeological metabolic bone disease indicators reported in other parts of the Pacific¹¹,¹². Secondly, the extensive cortical bone porosity leading to trabecularisation in our sample of archaeological Tongan females suggests experiences of oestrogen loss, mirroring modern day female osteoporosis incidence. This may further suggest at least the females represented in our sample lived to menopause age. The occurrence of trabecularisation in some of the male samples further points to advanced age despite widely held assumptions that past human longevity was shorter than today. We are unable to ascertain that archaeological Tongans specifically suffered from obesity, but we explain our findings through a differential diagnosis considering various aetiologies that might have impacted these individuals—complementarily or in an alternate fashion.

**Occurrence of metabolic abnormalities as a result of lifestyle and environmental factors.** The examination of human remains from archaeological backgrounds, such as tracing the origins of tuberculosis or leprosy⁹⁰, has revolutionised our current understanding of disease in the modern world. Osteoporosis is often considered a modern disease because it mostly manifests in the elderly and is strongly linked to modern sedentary lifestyles. However, its incidence was previously confirmed in various archaeological European⁵²⁻⁵⁷, North American⁵⁸, and African⁵⁹⁻⁶¹ osteological assemblages. We now report the first microscopic record of abnormally porous human bones from the Pacific region, adding to the growing body of evidence that modern metabolic conditions in the Pacific may be rooted in the past¹¹,¹².

Talasiu people are some of the earliest representatives of the Polynesian society emerging after the first settlement of archipelago at c. 2850 cal BP⁸⁴. The estimated Tongan skeletal fragility presented here is consistent with a historic basis for other current highly prevalent widespread metabolic problems such as obesity and type
## PHOSPHATE AND CARBONATES, AND THEIR RATIOS

<figure><img src="image_10.png" /></figure>

## CORTICAL BONE TRABECULARISATION AND CORTICAL WALL THINNING

<figure><img src="image_11.png" /><figcaption>**Figure 5.** Statistically significant Mann Whitney *U* ranking comparisons of the sFTIRM data for phosphate peaks (a), carbonate peaks and integrated area (b–e), and ratios computed using carbonate $v_2$: phosphate area (f). A strong negative correlation between cortical area and trabecularised cortical bone in the whole sample is shown (g).</figcaption></figure>

2 diabetes in the region³⁻⁵. Published evidence also exists for other types of skeletal abnormalities reported from different archaeological sites in Tonga⁹¹⁻⁹⁴, and Vanuatu¹¹,¹²,⁹⁵. For example, scurvy and hypervitaminosis A, and iron-deficiency anaemia were discussed amongst several possible conditions arising in archaeological Tongan communities possibly as a result of poor nutrition⁹¹. This came from skeletal evidence in 17 children aged 6 months to 3 years old suggesting a combination of infections and metabolic bone disease at the 'Atele pre-European burial mounds dated to AD 1100–1250 from Western Tonga⁹¹. Further, poorly mineralised 21 archaeological deciduous teeth obtained from infant skeletons excavated at 'Atele in Tonga also indicated developmental disturbances in dental enamel linked to physiological stress arising from possible infections and nutritional deficiencies⁹². Stable isotope data extracted from archaeological human Tongan bone suggested diets to be predominantly based on starchy and marine food sources⁸³,⁹³. Selection of a wider variety of food sources would have been limited by the island size in addition to norms arising from social organisation⁸³. Another site in Tonga, ca. 750–150 BP To-At-36 at Ha 'ateiho, yielded examples of dental developmental disturbance and growth
<table><caption>**Table 2.** Histomorphometry and synchrotron source Fourier-transform infrared microspectroscopy descriptive data per estimated sex (four males and four females). SD: standard deviation, min.: minimum data, max.: maximum data, A-P: antero-posterior, M-L: medio-lateral, H.Dn: Haversian canal density per mm², %Po.Ar: percent area of sample impacted by porosity producing trabecularisation.</caption><thead><tr><th>Estimated sex</th><th>N</th><th>Min</th><th>Max</th><th>Mean</th><th>SD</th></tr></thead><tbody><tr><td colspan="6"><strong>Femur size and quantitative histology</strong></td></tr><tr><td colspan="6"><strong>Female</strong></td></tr><tr><td>A-P midshaft diameter (mm)</td><td>4</td><td>25.28</td><td>29.30</td><td>27.26</td><td>2.15</td></tr><tr><td>M-L midshaft diameter (mm)</td><td>4</td><td>23.26</td><td>26.54</td><td>24.51</td><td>1.54</td></tr><tr><td>Midshaft circumference (mm)</td><td>4</td><td>76.00</td><td>86.00</td><td>79.75</td><td>4.79</td></tr><tr><td>H.Dn (/mm²)</td><td>4</td><td>12.216</td><td>22.223</td><td>15.60</td><td>4.54</td></tr><tr><td>%Po.Ar</td><td>4</td><td>13.908</td><td>50.232</td><td>35.508</td><td>16.397</td></tr><tr><td colspan="6"><strong>Male</strong></td></tr><tr><td>A-P midshaft diameter (mm)</td><td>4</td><td>27.72</td><td>33.87</td><td>30.75</td><td>2.87</td></tr><tr><td>M-L midshaft diameter (mm)</td><td>4</td><td>22.17</td><td>29.27</td><td>26.42</td><td>3.02</td></tr><tr><td>Midshaft circumference (mm)</td><td>4</td><td>86.00</td><td>101.00</td><td>92.25</td><td>7.09</td></tr><tr><td>H.Dn (/mm²)</td><td>4</td><td>15.63</td><td>20.80</td><td>17.45</td><td>2.295</td></tr><tr><td>%Po.Ar</td><td>4</td><td>0.000</td><td>51.793</td><td>15.813</td><td>24.171</td></tr><tr><td colspan="6"><strong>Fourier-transform infrared microspectroscopy</strong></td></tr><tr><td colspan="6"><strong>Female</strong></td></tr><tr><td>Peak phosphate height</td><td>848</td><td>0.201</td><td>1.139</td><td>0.481</td><td>0.189</td></tr><tr><td>Phosphate integrated area</td><td>848</td><td>15.087</td><td>71.381</td><td>34.518</td><td>11.673</td></tr><tr><td>Carbonate v<sub>2</sub> height</td><td>1676</td><td>0.001</td><td>0.216</td><td>0.023</td><td>0.018</td></tr><tr><td>Carbonate v<sub>2</sub> integrated area</td><td>2159</td><td>0.001</td><td>2.899</td><td>0.392</td><td>0.327</td></tr><tr><td>Carbonate v<sub>3</sub> height</td><td>1676</td><td>0.001</td><td>0.050</td><td>0.011</td><td>0.008</td></tr><tr><td>Carbonate v<sub>3</sub> integrated area</td><td>2071</td><td>0.001</td><td>2.873</td><td>0.475</td><td>0.498</td></tr><tr><td>Carbonate v<sub>2</sub>, phosphate ratio</td><td>846</td><td>0.001</td><td>0.143</td><td>0.019</td><td>0.007</td></tr><tr><td colspan="6"><strong>Male</strong></td></tr><tr><td>Peak phosphate height</td><td>2496</td><td>0.200</td><td>1.156</td><td>0.494</td><td>0.171</td></tr><tr><td>Phosphate integrated area</td><td>2496</td><td>14.937</td><td>72.184</td><td>34.899</td><td>10.346</td></tr><tr><td>Carbonate v<sub>2</sub> height</td><td>3339</td><td>0.001</td><td>0.104</td><td>0.033</td><td>0.017</td></tr><tr><td>Carbonate v<sub>2</sub> integrated area</td><td>3339</td><td>0.003</td><td>1.625</td><td>0.576</td><td>0.303</td></tr><tr><td>Carbonate v<sub>3</sub> height</td><td>3329</td><td>0.001</td><td>0.352</td><td>0.017</td><td>0.011</td></tr><tr><td>Carbonate v<sub>3</sub> integrated area</td><td>3264</td><td>0.001</td><td>8.510</td><td>0.850</td><td>0.526</td></tr><tr><td>Carbonate v<sub>2</sub>, phosphate ratio</td><td>2495</td><td>0.001</td><td>0.038</td><td>0.020</td><td>0.003</td></tr></tbody></table>

disruption indicators in 33 adult and 11 juveniles⁹⁴. Lower spine degeneration (spondylosis) was also observed on archaeological Tongan lumbar vertebrae⁹⁴. While spondylosis occurring as a result of bone degeneration due to age can be accelerated through weight bearing from active lifestyles⁹⁴ (combat sports at archaeological Tonga mounds were previously also discussed⁹⁶, and see our next discussion section about physical activity links to age-related bone loss), clinical studies have also linked it to increased body weight in obesity⁹⁷. The experiences of degenerative bone conditions in our sample of the Tongan individuals can be inferred from joint lesions indicating osteoarthritis (Table S1).

In Vanuatu, an archipelago western to Tonga, previous reports of most likely diagnoses of scurvy due to vitamin C deficiencies⁹⁵, gouty arthritis and DISH¹¹,¹², point to significant metabolic and nutritional issues afflicting Pacific Islanders thousands of years ago. While we cannot provide one aetiological explanation for the occurrence of bone metabolic conditions in our sample, we propose that a combination of lifestyle factors such as nutrition, and genetic predisposition played a role in the poor Tongan bone metabolism. Theoretical attempts to explain the high incidence of metabolic disease in the Pacific have predominantly focused on genetics (see⁹⁸,⁹⁹ for discussion). “The Thrifty Genotype” hypothesis proposes that thrifty genes, which would have been advantageous during famine events in the human past, drive diabetes and obesity in modern environments⁹⁹. However, considering our data, and the aforementioned archaeological skeletal data from other sites in the region¹¹,¹²,⁹¹⁻⁹⁵, it is clear that genes alone cannot elucidate the human obesity trends in the Pacific. Additionally, as noted previously by Gosling et al⁹⁸, thrifty genes cannot alone account for a range of selective pressures that characterised the Pacific islands. Alternative, or complementary, interpretations to consider include early childhood exposure to other diseases⁹¹,⁹²,¹⁰⁰, compromising adult immunity which can equally result in altered bone turnover and poor skeletal mineralisation¹⁰⁰,¹⁰¹. Ongoing macroscopic pathology examination of the Talasiu sample will provide further data and help shed light on possible incidence of conditions such as DISH or gout, expanding our reported microscopic bone loss markers.
**Sex-specific and ageing driven bone fragility.** Osteoporosis has a long history of afflicting females more than males due to menopause-driven loss of estrogen³¹,³². Oestrogen plays a key role in inhibiting osteoclast-mediated bone resorption, by regulating osteoclast apoptosis³¹. After menopause, prolonged bone loss occurs leading to increasingly porous and weak bones. Females are biologically disadvantaged because of this phenomenon, with modern postmenopausal women experiencing four times the level of osteoporosis than men¹⁰². Our data for a sample of archaeological Tongan females appear to match this sex-specific difference in bone fragility, providing another line of interpretation as an alternative or complementary aetiology. All our female samples showed a trabecularisation effect whereby the compact bone exhibited trabeculae-like architecture accompanied by cortical wall thinning (Figs. 1, 2). Female samples exhibiting higher cortical bone loss than male samples as a result of sex differences is consistent with prior archaeological⁵⁸ and clinical reports¹⁰³⁻¹⁰⁵. The cortical wall bone not impacted by abnormal/trabecularised porosity in the female samples showed a lower density of Haversian canals, which can be interpreted as a proxy for the amount of remodelled bone⁵⁸,⁸⁵. This can be explained by males possibly experiencing higher mechanically stimulated bone remodelling than females, body size differences, and/or within-sample age differences as noted in prior studies⁷³,¹⁰³,¹⁰⁴,¹⁰⁶.

Inter-woven with the sex-specific differences in bone microstructure is ageing⁶⁵,⁷²,¹⁰³⁻¹⁰⁵. It is well understood that the variability of age impacts on bone structure with sex manifests substantially in the femoral cross-section (and in other lower limb bones)¹⁰⁷. Ageing of bone tissue results in a greater resorption on the endocortical surface, leading to modified long bone cross-sectional geometric properties¹⁰⁷. In a classic study examining this in a sample of US human cadavers¹⁰⁷, age-related endocortical resorption manifested both in males and females. However, males appeared to exhibit a simultaneous formation of bone sub-periosteally and resorption endo-cortically, which did not drastically impact bone strength¹⁰⁷. In females, while the medullary cavity expanded, there was no associated expansion of the sub-periosteal bone, meaning that their cross-sectional properties weakened with age¹⁰⁷. When considering preindustrial humans such as the Pecos who undertook high levels of physical activity¹⁰⁸, sex-differences in bone loss while present were not as extensive as those seen in the aforementioned cadavers¹⁰⁷. In the Pecos, sub-periosteal bone and biomechanical properties of the femur (and tibia) followed a general pattern of increase with age in both sexes¹⁰⁸. We were not able to undertake a biomechanical analysis of cross-sectional geometry, but the substantial thinning of cortical wall in the Tongan female samples, along with lower Haversian canal densities, could suggest not enough mechanically stimulated bone remodelling around the time of their peak bone mass accrual phase of the life-course (third life decade)¹⁰⁹. Indeed, sedentary lifestyles are so widespread in Tonga today that several modern intervention efforts targeting adolescents have been unsuccessful¹¹⁰. This raises an important consideration for future studies of archaeological bone loss comparing populations from different parts of the world, whereby environmental and lifestyle factors impacting bone building in the early ontogeny will differ in accordance to genetic and cultural determinants¹⁰⁹.

One of the negative repercussions of overall 'improved' longevity and mortality across contemporary human populations is the advanced age of soft and hard tissues that results in multiple degenerative diseases experienced by the elderly¹¹¹. Bone fragility, increased fracture incidence, and hip-replacement surgeries, are some of the most common issues for the elderly of significantly deteriorated skeletal health. One of our female individuals (ID: Sk3.1) showed highly advanced skeletal deterioration characteristics (e.g. loss of all mandibular teeth antemortem, Figure S2) so she was classified as the oldest, and likely the only elderly, individual in our sample. This implies she might have survived into her sixth life decade (and possibly beyond). The combined age-at-death data and her abnormal intra-cortical porosity, are evidence that, at least some, archaeological Tongan females might have well surpassed the commonly assumed short longevity of humans in the distant past¹¹²,¹¹³. The trabecularisation of the cortical bone in our samples is remarkably similar to histology and microradiography images of midshafts examined in modern Australians aged 89 year old⁸⁸, and 67, 78, 90 year old¹¹⁴. This encourages further research combining bone biology, skeletal anatomy, archaeology, and social structures of Tonga to elucidate aspects of care of elderly in past Polynesia¹¹⁵.

Archaeological context-specific and methodological limitations of our study mean we cannot provide a single diagnosis and aetiology of the bone loss characteristics reported. Access to hundreds of well-preserved skeletons, as has been the case in some prior archaeological osteoporosis studies⁵⁴,⁵⁵, is not possible at Talasiu as the site is on a remote island constrained by its land and population size. Given the archaeological age of the Talasiu samples, the preservation of bone is not comparable to modern or post-mortem tissue, and thus cannot be experimentally examined to the same level. Future microscopic research applied to archaeological human remains from the Pacific will hopefully generate more comparative data, ideally using complementary 3D and 2D methods where feasible. While our study is limited by the broad age-at-death ranges, this is an issue impossible to overcome in biological anthropology as a narrow and exact chronological age for archaeological individuals cannot be ascertained from gross osteological analysis alone¹¹⁶. The assignment of biological sex is also a probability estimate based on well-established sexually dimorphic features of the human skeleton, which, without future aDNA validation, will be the best sex estimate possible.

## Conclusions

Our results are the first microscopic record of cortical bone loss in a sample of archaeological humans from Tonga. Given the small sample size, we suggest caution in the generalisation of our results in regards to the wider archaeological Tongan societies. With larger archaeological samples from across the Pacific islands, a pattern in bone loss may be shown in the future. Nevertheless, our results are evidence that a possible occurrence of bone metabolic conditions in a sample of archaeological Pacific individuals can be detected microscopically. We think this is an important step forward for discussions about metabolic diseases in the past and present Pacific. We discussed several explanations for the observed bone loss markers including: sex-specific driven bone micro-architectural deterioration with the possibility of our sample of Tongan females experiencing hormonal changes
due to menopause; age-related bone loss similar to that seen in modern aging populations; and lifestyle factors such as poor physical activity and nutrition. Equally, other diseases in the region compromising immunity might have contributed to abnormal bone physiology in our sample. Given how commonplace contemporary metabolic issues, including obesity and diabetes, are in Tonga, our key interpretation is that similar problems might have led to osteoporosis in our sample there ca. 3,000 years ago.

## Methods

We implemented an invasive methodology and so were restricted to eight femora (seven right and one left) (Table 1). The choice of femur side was determined on the basis of macroscopic preservation and availability of midshaft bone for extraction. Brief summarises of gross anatomically visible skeletal lesions per individual are reported in Table S1, with several individuals (6/8) showing evidence for osteoarthritis. While limited, this sample size is comparable to previous research implementing microscopic methods of bone assessment in archaeological samples<sup>77,85</sup>. The biological sex of the individuals represented by the femora had been previously estimated<sup>82,83</sup> following established methods<sup>117</sup>. Four males and four females (including one probable male and two probable females, Table 1) were estimated. The inclusion of 'probable' sexes is standard practice where archaeological human remains that do not show a directional consistency in sexually dimorphic skeletal features<sup>119</sup>. The adult age-at-death was assigned based on morphological features of the pelvic auricular surface<sup>118</sup>. The true chronological age of these individuals cannot be ascertained from bone morphology alone.

The skeletal remains represented by the specimens reported here were excavated and analysed with the permission from the Lapaha Community and Nobles, and the Ministry of Internal Affairs (Government of Tonga). The samples are curated at the Australian National University, Canberra until further repatriation notice. Talasiu on Tongatapu in the Kingdom of Tonga is one of the most archaeologically significant sites in Asia-Pacific. It represents some of the earliest occupations of the Neolithic people migrating into the Pacific alongside the Lapita culture some 3200-2850 cal BP<sup>120</sup>. It is a shoreline site located to the north of Lapaha, and is mainly comprised of shell middens and ceramic deposits<sup>115</sup>. The first excavation of Talasiu had taken place in 1957 which was subsequently followed by dating of the recovered shellfish to 2800 +/- 70 BP<sup>120</sup>. Excavations from 2008 onwards revealed burials containing human remains, including burnt remains that provided an insight into early mortuary behaviours in the region<sup>81</sup>. At least 19 late-Lapita/ immediately post-Lapita single to multiple burials have been reported, detailing a possible total number of 45 individuals represented and complex mortuary practices<sup>81,82</sup>. Shell midden samples from Talasiu indicate a sedentary human occupation dated to approximately 2700-2650 cal BP<sup>120</sup>, with more recent calibrations of 2650 BP for the burials<sup>82</sup>. Twenty-one of these individuals were recently examined as part of a multi-isotopic analysis reconstructing the first Polynesian diet<sup>83</sup>. A sub-sample of these Talasiu human individuals with well-preserved femora was studied here to conduct the analysis of bone loss markers. Applying the microscopic methods reported here to bones previously sampled for isotopic analysis, and limiting sample extraction to a cortical quadrant (rather than a complete long bone cross-section) ensures minimal invasion of the archaeological human material as per ethical recommendations<sup>76</sup>.

Samples were collected from the midshaft femur due to its biomechanical versatility, and prior published data reporting abnormal porosity occurring there<sup>103,104</sup>. Prior to sampling, each femur was photographed and measured at midshaft using standard WorkZone calipers and a measure tape to obtain A-P and M-L diameters (in mm), and the midshaft circumference (in mm) (Table 1). As the femora were fragmented on the distal and proximal ends, we could not measure maximum femoral lengths (Fig. 1). However, the determination of the midshaft location was easily achieved by locating the *linea aspera*, which is also the region from which samples were extracted (Fig. 1). This followed prior methods examining femur bone histology in archaeological specimens<sup>75,85</sup>, and reports of increased porosity on the posterior aspect of the femur in modern cadavers<sup>103</sup>. Approximately 1 cm thick cortical quadrants were removed using a Dremel 200/2-30 Two-Speed Rotary Tool equipped with a rotary blade. Cutting was performed as per prior methodologies<sup>121</sup>. The samples detached loosely after the blade had reached the medullary cavity having made two parallel longitudinal and transverse cuts on the bone exterior.

The preparation for three different microscopic analyses (2D histology, 3D confocal laser topography scanning, sFTIRM, Fig. 1) included sequential partitioning of bone samples into slices. Using the Dremel blade, each sample was further cut in half to designate an approximately 0.5 cm thick bone section for confocal topography scanning. The remainder of the sample was embedded in Buehler epoxy resin to impregnate the internal bone structure for histology and sFTIRM. Using a Kemet MICRACUT precision cutter equipped with a 150 mm diamond blade, a ~ 150 µm thin slice of bone was removed from the embedded block and designated for sFTIRM. The remaining portion of the sample was set aside for histological preparation.

The preparation of thin sections followed standard methods applicable to archaeological human bone<sup>75</sup>. The histology surface revealed by cutting on the low speed cutter was glued to a microscope glass slide using Stuk epoxy glue. It was further trimmed on the saw and ground using a series of silicon carbide pads until ~ 100 µm thickness was reached. The sections were then polished using Buehler MicroPolish alumina powder and cleared in an ultrasonic bath. This was followed by dehydration in ethanol, clearing using xylene, and mounting with cover slips. Imaging of the thin sections was undertaken using a high powered Olympus BX53 microscope, a DP74 camera, and associated Olympus cellSens Life Science Technology software. Images were scanned under a 40× total magnification (17 mm working distance) and 'auto-stitched' using the "Multiple Image Alignment—MIA" tool available from the Olympus cellSens software. Each image was then imported into Adobe Photoshop CC 2014 (replicated recently in Photoshop CC 2020), greyscaled (Black & White), converted to type 16 bit, and thresholded so that Haversian canals throughout cortical bone were 'enhanced'. This essentially converted cortical bone areas into black pixels and porous bone space into empty/white areas of the image (Fig. 2). Thresholding greyscale images for microscopic quantitative analysis is common practice<sup>122,123</sup> and can be applied to threshold
out Haversian canals¹²⁴. The images were imported into the open access ImageJ vol. 1.52 software for measurements of quantifiable components of the thin sections.

In order to differentiate between porosity producing trabecularisation and Haversian canal density the sections had to be manually segmented. We followed descriptions in<sup>88,89</sup> where a delimitation between 'dense' cortical bone and porosity producing trabecularisation cortical space can be estimated by visually separating the two bone matrices (see red dashed line in Fig. 2). As such, five variables as reported in previous bone porosity studies<sup>123-126</sup> (some acronyms were modified following the nomenclature standards by Dempster et al<sup>127</sup>) were collected. Total Bone Area (T.B.Ar in mm²)<sup>123</sup> was measured using the ImageJ vol. 1.52 "Polygon" tool by tracing the outer outline endosteal and periosteal borders of the section. Porosity Area (Po.Ar in mm²)<sup>124</sup> was measured using the same tools, but by manually selecting intra-cortical bone regions characterised with abnormal 'giant pores' that originated from within the endo-cortical section areas, which we here define as porosity producing trabecularisation. Cortical area (Ct.Ar) was the cortical bone region comprising cortical walls where no porosity producing trabecularisation (Po.Ar) was noted (Ct.Ar = T.B.Ar - Po.Ar). From these measurements, we calculated porosity producing trabecularisation (%Po.Ar = P.Ar/T.B.Ar × 100). Therefore, if no trabecularised regions are consistently observed, this method allows for a null result of Po.Ar and %Po.Ar. From within the Ct.Ar, Haversian canal number (H.N) was counted manually<sup>125,126</sup> (using the "Multi-Point" tool in ImageJ vol. 1.52 see red dots in Fig. 1). To estimate Haversian canal density (H.Dn)<sup>125,126</sup>, H.N was divided by Ct.Ar to obtain a value per mm². Therefore, we took into consideration two types of bone porosity measures—H.Dn as a proxy for the amount of remodelled bone present in cortical bone, and %Po.Ar which estimates the area of intra-cortical bone affected by trabecularised cortex. Through H.Dn, we worked with an assumption that one Haversian canal represents one secondary osteon. Because cement lines of secondary osteons were not consistently preserved in these archaeological samples, this technique does not account for fragmentary osteons (Fig. 2). The canals that were counted were identified as Haversian canals to the best of our expertise. This excluded micro-features that resembled other pores, which might have occurred as a result of diagenesis. That way, %Po.Ar is exclusively composed of 'giant' and coalescing pores, occurring consistently intra-cortically, originating on the endo-cortical part of the sample. These regions should otherwise be filled with dense cortical bone if the individual did not experience significant bone loss. Both measures can provide an insight into a Bone Multi Cellular unit (BMU) activity tuningling through cortical bone<sup>128</sup>, whereby H.Dn approximates the number of BMUs that once existed per mm², and %Po.Ar indicates prolonged bone resorption through osteoclast-mediated activity.

To provide a qualitative illustration of bone surface topography in relation to the porosity producing trabecularisation, we scanned two contrasting male and female samples (IDs: Sk3.1 and Sk3.2) using an OLS5000 3D laser confocal microscope. Confocal laser scanning microscopy is a recommended technique for characterising porous structures such as bones<sup>129</sup>. The associated OLS5000 2017 LEXT data acquisition and data analysis application software (Olympus LEXT, Japan) (Fig. 1) was then used to apply a heat map of false colours that ranged from red to blue indicating highest to lowest depth, respectively (Fig. 3, S10). This resulted in yellow to green colours indicating low topography (Fig. 3, S11) in the male, and red to blue marking high topography in the female (Fig. 3, S10). We used a 5×LEXT short working distance objective (20 mm) with a 405 nm violet laser that scans 4,096 pixels along the x-axis, with the zoom as at 1.0x. We used the LEXT automated stitching tool to collectively scan six regions of bone located on the mid-line of the cortical area of the posterior bone quadrant. The scanned area of the female sample was 7192.648 µm long (y-axis) and 4872.699 µm wide (x-axis). The area scanned on the male sample was 7185.926 µm long (y-axis) and 4876.763 µm wide (x-axis). The z-plane depth (height of laser reaching the bone surface) was approximately 3741.925 µm.

Each sample was examined for phosphate and carbonate content using sFTIRM<sup>86,130-136</sup>. Phosphate was selected because it is necessary for bone metabolism, and occurs in skeletal tissue as part of hydroxyapatite crystals<sup>134</sup>. Carbonates ($v_2$ and $v_3$) substitute calcium apatite and influence the properties of crystal in bone, and thus influence bone function at the macroscopic scale<sup>133</sup>. Carbonate in bone is mainly A-type whereby it substitutes for phosphate, but it can also be naturally accompanied by B-type (substituting for hydroxide)<sup>133</sup>. Assessing total (A- and B-type combined) carbonate content, in addition to phosphate, and $v_2$ carbonate and phosphate ratio of integrated areas under trace, can thus provide an insight into the extent to which calcium phosphate has been resorbed through osteoclast-mediated activity,<sup>86,135,136</sup>

The samples were scanned at the IRM beamline at the Australian Synchrotron facility in Melbourne (Victoria). The scanning technique followed settings reported by Vrahnas et al.<sup>86</sup> with the exception that we used an attenuated total reflectance (ATR) attachment that allowed a direct contact with bone surface<sup>85,87</sup>. The synchrotron light source provides a highly intense infrared beam that was used to analyse mineral content of bone in situ and the ATR attachment allows the collection of infrared (IR) data from sample sections that are otherwise too thick for conventional IR transmission analysis<sup>87</sup>. The sFTIRM measurement was performed using a Bruker V80v FTIR spectrometer and a Hyperion 3000 IR microscope, which produced high quality spectra in terms of signal-to-noise ratios at 1–2 µm spatial resolution when the synchrotron IR beam coupled to the ATR crystal<sup>87</sup>.

Each sample had four regions of interest (ROIs) identified on the sub-periosteal area of bone, avoiding the endocortical surfaces particularly that they are so extensively affected by abnormal porosity in the Talasius samples (Fig. 1). Before placing the samples on the microscope stage, the ROIs were identified from images captured at 2×magnification under a basic dissecting microscope (Fig. 1). We measured carbonates $v_2$ (890–850 cm⁻¹) and $v_3$ (1500–1400 cm⁻¹), and phosphate (1180–916 cm⁻¹), which were then used to calculate carbonate($v_2$):phosphate ratios (890–850 cm⁻¹:1180–916 cm⁻¹)<sup>85,86,132,133</sup>. Data analysis was undertaken in OPUS 7.2 and 8.0.19 (Bruker Optik, Germany) by creating integration files and extracting peak height (absorbance/AU) and area under the trace values in each spectrum. Each ROI was scanned for 220 spectra, which totalled 880 spectra per sample, totalling 6820 spectra in the entire sample (this excludes the second region of interest, which was unsuccessfully scanned, therefore deducting 220 spectra in data for individual Sk 3.1).
Once data were inspected, it became apparent that not all spectra were suitable for analysis. This might have occurred as a result of the ATR attachment not making contact with all scanned bone regions. Additionally, as our samples derive from an archaeological context, and are thus impacted by diagenetic processes, diagenetic calcite deposited post-mortem can contribute to carbonate sites of the spectra. To account for these issues, the spectra were inspected and those of poor quality were removed. Our supplemental file (Dataset S1) with raw data reports all spectra collected during the scanning, as well as the datasets following inspection. The inspection involved removing phosphate peak values (and corresponding integrated areas under the peak) of <0.2. All carbonate spectra values in the negative range to 0.001 were also removed. To further assess for possible diagenesis impacting the carbonate data (sites where calcite can occur), statistical correlations were performed to check how well carbonate $v_2$ and $v_3$ data aligned within each specimen (Table S6, S7). Where no strong correlations were identified, this was taken as an indication that diagenesis might have impacted the data as $v_3$ and $v_2$ values are scattered randomly. We also visually examined the sites of $v_3$ for a strong peak at $1427 cm^{-1}$ indicating presence of calcite, which should otherwise be replaced by a trough in apatite spectra (creating a doublet)¹³⁷ (Figure S3–S9). Overall, the inspection reduced the number of suitable phosphate spectra to 3344 in males and 844 in females. Carbonate spectra were minimum 1676 ($v_2$) and 2159 ($v_3$) in females, and 3339 ($v_2$) and 3329 ($v_3$) in males (Table 2). The carbonate $v_2$:phosphate ratios calculated from integrated areas under the trace were possible to compute for 846 in females and 2495 in males. We note that all data for Sk12 appeared unsuitable, and so this specimen was excluded from data comparisons.

All data obtained in our study were firstly summarised descriptively. Descriptive comparisons were made within-sample by contrasting results between the sexes, and identifying whether either of the groups showed relatively higher or lower values of data. Next, as our sample size is small (n = 8), we were limited in the choice of statistical analyses to conduct inferential tests. On this basis, we applied non-parametric tests when comparing the sex groups (using a Mann Whitney U test¹³⁸) and a Spearman's Rho correlation to test if %Po.Ar was significantly associated with cortical thinning. All inferential statistical testing was conducted in IBM Statistical Package for Social Sciences 26 (SPSS) software. For the sFTIRM spectra inspection through correlations, data normality Kolmogorov-Smirnov tests were performed to check whether non-parametric or parametric correlations tests should be applied (Table S7). This step informed the use of Spearman's Rho correlations as the spectra values were not normally distributed (Table S7, S8). Strong correlations ($Rho = 0.68–1.00$)¹³⁹ were taken as an indication that data can be interpreted to draw conclusions for our research question.

## Data availability

All data generated or analysed during this study are included in this published article (and its Supplementary Information files).

Received: 26 June 2020; Accepted: 25 January 2021

Published online: 17 February 2021

## References

1. Swinburn, B. *et al.* The global syndemic of obesity, undernutrition, and climate change: the Lancet Commission report. *The Lancet* **23**, 791–846. https://doi.org/10.1016/S0140-6736(18)32822-8 (2019).

2. NCD Risk Factor Collaboration (NCD-RisC), Trends in adult body-mass index in 200 countries from 1975 to 2014: a pooled analysis of 1698 population-based measurement studies with 19.2 million participants. *The Lancet* **387**(1377–96), 2016. https://doi.org/10.1016/S0140-6736(16)30054-X (2016).

3. Hawley, N. L. & McGarvey, S. T. Obesity and diabetes in Pacific Islanders: the current burden and the need for urgent action. *Curr. Diab. Rep.* **15**, 29. https://doi.org/10.1007/s11892-015-0594-5 (2015).

4. Lin, S. *et al.* Diabetes and obesity trends in Tonga over 40 years. *Asia Pac. J. Public Health* **28**, 475–485. https://doi.org/10.1177/1010539516645156 (2016).

5. Furusawa, T. *et al.* The Q223R polymorphism in LEPR is associated with obesity in Pacific Islanders. *Hum. Genet.* **127**, 287–294. https://doi.org/10.1007/s00439-009-0768-9 (2010).

6. Hetherington, M. M. & Cecil, J. E. Gene-environment interactions in obesity. In *Frontiers in Eating and Weight Regulation* (eds Langhans, W. & Geary, N.) 195–203 (Krager Publishers, Basel, 2010).

7. Snowdon, W., Lawrence, M., Schultz, J., Vivili, P. & Swinburn, B. Evidence-informed process to identify policies that will promote a healthy food environment in the Pacific Islands. *Public Health Nutr.* **13**, 886–892. https://doi.org/10.1017/S136898001000011X (2010).

8. Freese, J., Klement, R. J., Ruiz-Núñez, B., Schwarz, S. & Lötzerich, H. The sedentary (r)evolution: Have we lost our metabolic flexibility. *F1000Res.* **6**, 1787. https://doi.org/10.12688/f1000research.12724.2 (2017).

9. Waldron, T. *Palaeopathology* (Cambridge University Press, Cambridge, 2008).

10. Brickley, M. & Ives, R. *The Bioarchaeology of Metabolic Bone Disease* (Elsevier Academic Press, Amsterdam, 2008).

11. Buckley, H. R. Possible gouty arthritis in Lapita-associated skeletons from Teouma, Efate Island, Central Vanuatu. *Curr. Anthropol.* **48**, 741–749. https://doi.org/10.1086/520967 (2007).

12. Foster, A. *et al.* Possible diffuse idiopathic skeletal hyperostosis (DISH) in a 3000-year-old Pacific Island skeletal assemblage. *J. Archaeol. Sci. Rep.* **18**, 408–419. https://doi.org/10.1016/j.jasrep.2018.01.002 (2018).

13. Bateman, M., Hapuarachchi, K., Pinto, C. & Doyle, A. J. Diffuse idiopathic skeletal hyperostosis (DISH): increased prevalence in Pacific Islanders. *J. Med. Imag. Rad. Oncol.* **6**, 188–193. https://doi.org/10.1111/1754-9485.12679 (2018).

14. Miszkiewicz, J. J., Brennan-Olsen, S. & Riancho, J. A. *Bone Health: A Reflection of the Social Mosaic* (Springer Nature Medicine, Singapore, 2019).

15. Miszkiewicz, J. J. & Cooke, K. M. Socio-economic determinants of bone health from past to present. *Clin. Rev. Bone Miner. Metab.* **17**, 109–122. https://doi.org/10.1007/s12018-019-09263-1 (2019).

16. Eaton, S. B. & Eaton, S. B. Physical inactivity, obesity, and type 2 diabetes: an evolutionary perspective. *Res. Q. Exercise Sport* **88**, 1–8. https://doi.org/10.1080/02701367.2016.1268519 (2017).

17. Pagnotti, G. M. *et al.* Combating osteoporosis and obesity with exercise: leveraging cell mechanosensitivity. *Nat. Rev. Endocrinol* **15**, 339–355. https://doi.org/10.1038/s41574-019-0170-1 (2019).

18. Patsch, J. M. *et al.* Increased bone resorption and impaired bone microarchitecture in short-term and extended high-fat diet-induced obesity. *Metabolism* **60**, 243–249. https://doi.org/10.1016/j.metabol.2009.11.023 (2011).
19. Cândido, F. G. & Bressan, J. Vitamin D: link between osteoporosis, obesity, and diabetes?. *Int. J. Mol. Sci.* **15**, 6569–6591. https://doi.org/10.3390/ijms15046569 (2014).

20. Peel, N. Orthopaedics I: general principles. Bone remodelling and disorders of bone metabolism. *Surgery (Oxford)* **27**, 70–74. https://doi.org/10.1016/j.mpsur.2008.12.007 (2009).

21. Martin, T. J. & Sims, N. A. Osteoclast-derived activity in the coupling of bone formation to resorption. *Trends Mol. Med.* **11**, 76–81. https://doi.org/10.1016/j.molmed.2004.12.004 (2005).

22. Jepsen, K. J., Bigelow, E. M., Ramcharan, M., Schlecht, S. H. & Karvonen-Gutierrez, C. A. Moving toward a prevention strategy for osteoporosis by giving a voice to a silent disease. *Women's Midlife Health* **2**, 1–5. https://doi.org/10.1186/s40695-016-0016-0 (2016).

23. Lewiecki, E. M. Osteoporotic fracture trends in a population of US managed care enrollees from 2007 to 2017. *Osteoporos. Int.* **31**, 1299–1304. https://doi.org/10.1007/s00198-020-05334-y (2020).

24. Nguyen, T. V. Lifestyle factors and bone density in the elderly: implications for osteoporosis prevention. *J. Bone Miner. Res.* **9**, 339–1346. https://doi.org/10.1002/jbmr.5650090904 (1994).

25. Hannan, M. T. Risk factors for longitudinal bone loss in elderly men and women: the Framingham Osteoporosis Study. *J. Bone Miner. Res.* **15**, 710–720. https://doi.org/10.1359/jbmr.2000.15.4.710 (2000).

26. Karasik, D. & Ferrari, S. L. Contribution of gender-specific genetic factors to osteoporosis risk. *Ann. Hum. Genet.* **72**, 696–714. https://doi.org/10.1016/S0140-6736(02)08657-9 (2008).

27. Harding, A. T. & Beck, B. R. Exercise, osteoporosis, and bone geometry. *Sports https://doi.org/10.3390/sports5020029 (2017).*

28. Hsu, Y. H., Xu, X. & Jeong, S. Genetic determinants and pharmacogenetics of osteoporosis and osteoporotic fracture. In *Osteoporosis Pathophysiology and Clinical Management* (eds Leder, B. Z. & Wein, M. N.) 485–506 (Cham, Springer Nature Switzerland, 2020). https://doi.org/10.1007/978-3-319-69287-6_25.

29. Bassett, J. D. et al. Thyroid hormone excess rather than thyrotropin deficiency induces osteoporosis in hyperthyroidism. *Mol. Endocrinol.* **21**, 1095–1107. https://doi.org/10.1210/me.2007-0033 (2007).

30. Anand, P. & Mehler, P. S. Osteoporosis recovery in severe anorexia nervosa: a case report. *J. Eat. Disord.* **7**, 38. https://doi.org/10.1186/s40337-019-0269-8 (2019).

31. Klein-Nulend, J., van Oers, R. F., Bakker, A. D. & Bacabac, R. G. Bone cell mechanosensitivity, estrogen deficiency, and osteoporosis. *J. Biomech.* **48**, 855–865. https://doi.org/10.1016/j.jbiomech.2014.12.007 (2015).

32. Shieh, A., Greendale, G. A., Cauley, J. A. & Karlamangla, A. S. The association between fast increase in bone turnover during the menopause transition and subsequent fracture. *J. Clin. Endocrinol. Metab.* **105**, e1440–e1448. https://doi.org/10.1210/cline/m/dgz281 (2020).

33. Zhao, L. J. et al. Relationship of obesity with osteoporosis. *J. Clin. Endocrinol. Metab.* **92**, 1640–1646. https://doi.org/10.1210/jc.2006-0572 (2007).

34. Rosen, C. J. & Bouxsein, M. L. Mechanisms of disease: is osteoporosis the obesity of bone?. *Nat. Clin. Pract. Rheumatol.* **2**, 35–43. https://doi.org/10.1038/npr rheum0070 (2006).

35. Greco, E. A. et al. Is obesity protective for osteoporosis? Evaluation of bone mineral density in individuals with high body mass index. *Int. J. Clin. Pract.* **64**, 817–820. https://doi.org/10.1111/j.1742-1241.2009.02301.x (2010).

36. Fassio, A. et al. The obesity paradox and osteoporosis. *Eat. Weight Disord. St.* **23**, 293–302. https://doi.org/10.1007/s40519-018-0505-2 (2018).

37. Murugananandan, S., Roman, A. A. & Sinal, C. J. Adipocyte differentiation of bone marrow-derived mesenchymal stem cells: cross talk with the osteoblastogenic program. *Cell. Mol. Life Sci.* **66**, 236–253. https://doi.org/10.1007/s00018-008-8429-z (2009).

38. Poiana, C., Carsote, M., Radoi, V., Mihai, A. & Capatina, C. Prevalent osteoporotic fractures in 622 obese and non-obese menopausal women. *J. Med. Life* **8**, 462 (2015).

39. Sun, H. B. et al. Preliminary study of obese patients with chronic obstructive pulmonary disease suffering from painful osteoporotic vertebral compression fracture treated by percutaneous vertebroplasty in improved prone position and right lateral position. *World Neurosurg.* **130**, e933–e940. https://doi.org/10.1016/j.wneu.2019.07.040 (2019).

40. Tanaka, S., Kuroda, T., Saito, M. & Shiraki, M. Overweight/obesity and underweight are both risk factors for osteoporotic fractures at different sites in Japanese postmenopausal women. *Osteoporos. Int.* **24**, 69–76. https://doi.org/10.1007/s00198-012-2209-1 (2013).

41. Johnell, O. & Kanis, J. A. An estimate of the worldwide prevalence and disability associated with osteoporotic fractures. *Osteoporosis Int.* **17**, 1726–1733. https://doi.org/10.1007/s00198-006-0172-4 (2006).

42. Stride, P. Polynesian bones. *J. Adv. Med. Med. Res.* **16**, 1–9. https://doi.org/10.9734/BJMMR/2016/25651 (2016).

43. Ebeling, P. R. et al. Secondary prevention of fragility fractures in Asia Pacific: an educational initiative. *Osteoporosis Int.* **31**, 805–826. https://doi.org/10.1007/s00198-019-05197-y (2019).

44. Elaine, W. Y. Osteoporosis screening and diagnosis. In *Osteoporosis. Contemporary Endocrinology* (eds Leder, B. & Wein, M.) 139–149 (Springer Nature Switzerland, Cham, 2020). https://doi.org/10.1007/978-3-319-69287-6_7.

45. Fuggle, N. R. et al. Fracture prediction, imaging and screening in osteoporosis. *Nat. Rev. Endocrinol.* **15**, 535–547. https://doi.org/10.1038/s41574-019-0220-8 (2019).

46. Kanis, J. A. Diagnosis of osteoporosis and assessment of fracture risk. *The Lancet* **359**, 1929–1936. https://doi.org/10.1016/S0140-6736(02)08761-5 (2002).

47. Link, T. M. Metabolic bone disease. In *Measurements in Musculoskeletal Radiology* (eds Cassar-Pullicino, V. & Davies, A.) 785–807 (Springer, Berlin, 2020). https://doi.org/10.1007/978-3-540-68897-6_19.

48. Lorente-Ramos, R. et al. Dual-energy x-ray absorptiometry in the diagnosis of osteoporosis: a practical guide. *Am. J. Roentgenol.* **196**, 897–904. https://doi.org/10.2214/AJR.10.5416 (2011).

49. Kulak, C. A. M. & Dempster, D. W. Bone histomorphometry: a concise review for endocrinologists and clinicians. *Arq. Bras. Endocrinol. Metabol.* **54**, 87–98. https://doi.org/10.1590/S0004-27302010000200002 (2010).

50. dos Santos Floter, M., Bittar, C. K., Zabeu, J. L. A. & Carneiro, A. C. R. Review of comparative studies between bone densitometry and quantitative ultrasound of the calcaneus in osteoporosis. *Acta. Reumatol. Port.* **36**, 327–335. https://doi.org/10.1590/S0004-27302010000200002 (2011).

51. Cardoso, L., Fritton, S. P., Gailani, G., Benalla, M. & Cowin, S. C. Advances in assessment of bone porosity, permeability and interstitial fluid flow. *J. Biomech.* **46**, 253–265. https://doi.org/10.1016/j.jbiomech.2012.10.025 (2013).

52. Mays, S., Turner-Walker, G. & Syversen, U. Osteoporosis in a population from medieval Norway. *Am. J. Phys. Anthropol.* **131**, 343–351. https://doi.org/10.1002/aipa.20445 (2006).

53. Mays, S. Age-dependent cortical bone loss in women from 18th and early 19th century London. *Am. J. Phys. Anthropol.* **112**, 349–361. https://doi.org/10.1002/1096-8644(200007)112:3%3c349::AID-AJPA6%3e3.0.CO;2-0 (2000).

54. Drusini, A. G., Bredariol, S., Carrara, N. & Rippa-Bonati, M. Cortical bone dynamics and age-related osteopenia in a Longobard archaeological sample from three graveyards of the Veneto region (northeast Italy). *Int. J. Ostearchaeol.* **10**, 268–279. https://doi.org/10.1002/1099-1212(200007)10:4%3c268::AID-OA529%3e3.0.CO;2-Q (2000).

55. Brickley, M. & Howell, P. G. T. Measurement of changes in trabecular bone structure with age in an archaeological population. *J. Archaeol. Sci.* **26**, 151–157. https://doi.org/10.1006/jasc.1998.0313 (1999).

56. Ekenman, I., Eriksson, S. A. V. & Lindgren, J. U. Bone density in medieval skeletons. *Calcif. Tissue Int.* **56**, 355–358. https://doi.org/10.1007/BF00301601 (1995).
57. González-Reimers, E. et al. Double-energy X-ray absorptiometry in the diagnosis of osteopenia in ancient skeletal remains. *Am. J. Phys. Anthropol.* **118**, 134–145. https://doi.org/10.1002/ajpa.10076 (2002).

58. Ericksen, M. F. Cortical bone loss with age in three Native American populations. *Am. J. Phys. Anthropol.* **45**, 443–452. https://doi.org/10.1002/ajpa.1330450306 (1976).

59. Dewey, J. R., Armelagos, G. J. & Bartley, M. H. Femoral cortical involution in three Nubian archaeological populations. *Hum. Biol.* **41**, 13–28 (1969).

60. Martin, D. L. & Armelagos, G. J. Morphometrics of compact bone: an example from Sudanese Nubia. *Am. J. Phys. Anthropol.* **51**, 571–577. https://doi.org/10.1002/ajpa.1330510409 (1979).

61. Zaki, M. E., Hussien, F. H. & El Banna, R. A. E. S. Osteoporosis among ancient Egyptians. *Int. J. Osteoarchaeol.* **19**, 78–89. https://doi.org/10.1002/oa.978 (2009).

62. Agarwal, S. C. & Stout, S. D. *Bone Loss and Osteoporosis: An Anthropological Perspective* (Springer, Boston, 2003).

63. Curate, F. Osteoporosis and paleopathology: a review. *J. Anthropol. Sci.* **92**, 119–146. https://doi.org/10.4436/jass.92003 (2014).

64. Agarwal, S. C. & Grynpas, M. D. Bone quantity and quality in past populations. *Anat. Rec.* **246**, 423–432. https://doi.org/10.1002/(SICI)1097-0185(199612)246:4%3c423::AID-AR1%3e3.0.CO;2-W (1996).

65. Agarwal, S. C., Dumitriu, M., Tomlinson, G. A. & Grynpas, M. D. Medieval trabecular bone architecture: the influence of age, sex, and lifestyle. *Am. J. Phys. Anthropol.* **124**, 33–44. https://doi.org/10.1002/ajpa.10335 (2004).

66. Curate, F. et al. A glimpse from the past: osteoporosis and osteoporotic fractures in a Portuguese identified skeletal sample. *Acta Reumatol. Port.* **38**, 20–27 (2013).

67. Brickley, M. An investigation of historical and archaeological evidence for age-related bone loss and osteoporosis. *Int. J. Osteoarchaeol.* **12**, 364–371. https://doi.org/10.1002/oa.635 (2002).

68. Morgan, B., Mant, M., de la Cova, C. & Brickley, M. B. Osteoporosis, osteomalacia, and hip fracture: a case study from the Terry collection. *Int. J. Paleopathol.* **30**, 17–21. https://doi.org/10.1016/j.ijpp.2020.03.004 (2020).

69. Curate, F., Lopes, C. & Cunha, E. A 14th–17th century osteoporotic hip fracture from the Santa Clara-a-Velha Convent in Coimbra (Portugal). *Int. J. Osteoarchaeol.* **20**, 591–596. https://doi.org/10.1002/oa.1076 (2010).

70. Roberts, C. & Wakely, J. Microscopical findings associated with the diagnosis of osteoporosis in palaeopathology. *Int. J. Osteoarchaeol.* **2**, 23–30. https://doi.org/10.1002/oa.1390020105 (1992).

71. Brickley, M. B. & Agarwal, S. C. Techniques for the investigation of age-related bone loss and osteoporosis in archaeological bone. In *Bone Loss and Osteoporosis: An Anthropological Perspective* (eds Agarwal, S. C. & Stout, S. D.) 157–172 (Springer, Boston, 2003). https://doi.org/10.1007/978-1-4419-8891-1_10.

72. Cho, H. & Stout, S. D. Bone remodeling and age-associated bone loss in the past: a histomorphometric analysis of the Imperial Roman skeletal population of Isola Sacra. In *Bone Loss and Osteoporosis: An Anthropological Perspective* (eds Agarwal, S. C. & Stout, S. D.) 207–228 (Springer, Boston, 2003). https://doi.org/10.1007/978-1-4419-8891-1_13.

73. Cooper, D. M. L., Matyas, J. R., Katzenberg, M. A. & Hallgrimsson, B. Comparison of microcomputed tomographic and micro-radiographic measurements of cortical bone porosity. *Calcif. Tissue Int.* **74**, 437–447. https://doi.org/10.1007/s00223-003-0071-z (2004).

74. Faccia, K. et al. Bone quality in prehistoric, cis-baikal forager femora: A micro-CT analysis of cortical canal microstructure. *Am. J. Phys. Anthropol.* **154**, 486–497. https://doi.org/10.1002/ajpa.22533 (2014).

75. Miszkiewicz, J. J. & Mahoney, P. Human bone and dental histology in an archaeological context. In *Human Remains: Another Dimension* (eds Thompson, T. & Errickson, D.) 29–43 (Elsevier Academic Press, London, The Application of Imaging to the Study of Human Remains, 2017). https://doi.org/10.1016/B978-0-12-804602-9.00004-7.

76. Mays, S., Elders, J., Humphrey, L., White, W., & Marshall, P. *Science and the Dead: A Guideline for the Destructive Sampling of Archaeological Human Remains for Scientific Analysis* (English Heritage Publishing with the Advisory Panel on the Archaeology of Burials in England, 2013).

77. Bennike, P., Bohr, H. & Toft, T. Determination of mineral content and organic matrix in bone samples using dual photon absorptiometry. *Int. J. Anthropol.* **8**, 111–116. https://doi.org/10.1007/BF02447620 (1993).

78. Booth, T. The rot sets. In: low-powered microscopic investigation of taphonomic changes to bone microstructure and its application to funerary contexts. In *Human Remains: Another Dimension. The Application of Imaging to the Study of Human Remains* (eds Thompson, T. & Errickson, D.) 7–28 (Elsevier Academic Press, London, 2017). https://doi.org/10.1016/B978-0-12-804602-9.00003-5.

79. Caruso, V. et al. Bone diagenesis in archaeological and contemporary human remains: an investigation of bone 3D microstructure and minero-chemical assessment. *Archaeol. Anthropol. Sci.* **12**, 162. https://doi.org/10.1007/s12520-020-01090-6 (2020).

80. Jeanson, A. L. et al. Body mass estimation from the skeleton: an evaluation of 11 methods. *Forensic Sci. Int.* **281**(183), e1-183. e8. https://doi.org/10.1016/j.forsciint.2017.10.026 (2017).

81. Valentin, F. & Clark, G. Early Polynesian mortuary behaviour at the Talasiu site, Kingdom of Tonga. *J. Pac. Archaeol.* **4**, 1–14 (2013).

82. Valentin, F., Clark, G., Parton, P. & Reepmeyer, C. The oldest Polynesian burial site (Talasiu site, Kingdom of Tonga): Implications for the reconstruction of the Ancestral Polynesian Society. *Antiquity* **94**, 999–1014 (2020).

83. Herrscher, E. et al. Multi-isotopic analysis of first Polynesian diet (Talasiu, Tongatapu, Kingdom of Tonga). *J. Archaeol. Sci. Rep.* **18**, 308–317. https://doi.org/10.1016/j.jasrep.2018.01.012 (2018).

84. Burley, D. V., Connaughton, S. P. & Clark, G. Early cessation of ceramic production for ancestral Polynesian society in Tonga. *PLoS ONE* **13**, e0193166. https://doi.org/10.1371/journal.pone.0193166 (2018).

85. Miszkiewicz, J. J. et al. Asymmetric midshaft femur remodeling in an adult male with left sided hip joint ankylosis, Metal Period Nagsabaran, Philippines. *Int. J. Paleopath.* **31**, 14–22. https://doi.org/10.1016/j.ijpp.2020.07.003 (2020).

86. Vrahnas, C. et al. Anabolic action of parathyroid hormone (PTH) does not compromise bone matrix mineral composition or maturation. *Bone* **93**, 146–154. https://doi.org/10.1016/j.bone.2016.09.022 (2016).

87. Vongsvivut, J. et al. Synchrotron macro ATR-FTIR microspectroscopy for high-resolution chemical mapping of single cells. *Analyst* **144**, 3226–3238. https://doi.org/10.1039/C8AN01543K3161-3452 (2019).

88. Lerebours, C., Buenzli, P. R., Scheiner, S. & Pivonka, P. A multiscale mechanobiological model of bone remodelling predicts site-specific bone loss in the femur during osteoporosis and mechanical disuse. *Biomech. Model. Mechanobiol.* **15**, 43–67. https://doi.org/10.1007/s10237-015-0705-x (2016).

89. Buenzli, P. R., Thomas, C. D. L., Clement, J. G. & Pivonka, P. Endocortical bone loss in osteoporosis: the role of bone surface availability. *Int. J. Numer. Method. Biomed. Eng.* **29**, 1307–1322. https://doi.org/10.1002/cnm.2567 (2013).

90. Expanding the dialogue between immunology and paleopathology. Crespo, F., White, J., & Roberts, C. Revisiting the tuberculosis and leprosy cross-immunity hypothesis. *Int. J. Paleopath.* **26**, 37–47. https://doi.org/10.1016/j.ijpp.2019.05.005 (2019).

91. Buckley, H. R. Subadult health and disease in prehistoric Tonga, Polynesia. *Am. J. Phys. Anthropol.* **113**, 481–505. https://doi.org/10.1002/1096-8644(200012)113:4%3c481::AID-AIPA3%3e3.0.CO;2-1 (2000).

92. Farah, R. A. et al. Physiological insult or the burial environment: Differentiating developmental defects from post-mortem stained enamel in deciduous dentition from the Chiefdom Period of Tonga, Polynesia. *Matters Sel.* https://doi.org/10.19185/matters.20160500005 (2016).

93. Stantis, C., Kinaston, R. L., Richards, M. P., Davidson, J. M. & Buckley, H. R. Assessing human diet and movement in the Tongan maritime chiefdom using isotopic analyses. *PLoS ONE* https://doi.org/10.1371/journal.pone.0123156 (2015).
94. Pietrusewsky, M., Douglas, M. T., Ikehara-Quebral, R. M. & Kadohiro Lauer, K. Skeletal and dental health of early Tongans: The bioarchaeology of the human skeletons from the To-At-36 Site, Ha 'ateiho, Tongatapu, Tonga. *J. Isl. Coast. Archaeol.* **15**, 204–243. https://doi.org/10.1080/15564894.2018.1564711 (2019).

95. Buckley, H. R. et al. Scurvy in a tropical paradise? Evaluating the possibility of infant and adult vitamin C deficiency in the Lapita skeletal sample of Teouma, Vanuatu, Pacific islands. *Int. J. Paleopath.* **5**, 72–85. https://doi.org/10.1016/j.ijpp.2014.03.001 (2014).

96. Scott, R. M. & Buckley, H. R. Exploring prehistoric violence in Tonga: Understanding skeletal trauma from a biocultural perspective. *Curr. Anthropol.* **55**, 335–347. https://doi.org/10.1086/676477 (2014).

97. Sheng, B., Feng, C., Zhang, D., Spitler, H. & Shi, L. Associations between obesity and spinal diseases: a medical expenditure panel study analysis. *Int. J. Environ. Res. Public Health* **14**, 183. https://doi.org/10.3390/ijerph14020183 (2017).

98. Gosling, A. L., Buckley, H. R., Matisoo-Smith, E. & Merriman, T. R. Pacific populations, metabolic disease and 'Just-So Stories': a critique of the 'Thrifty Genotype' hypothesis in Oceania. *Ann. Hum. Genet.* **79**, 470–480. https://doi.org/10.1111/ahg.12132 (2015).

99. Qasim, A. et al. On the origin of obesity: identifying the biological, environmental and cultural drivers of genetic risk among human populations. *Obes. Rev.* **19**, 121–149. https://doi.org/10.1111/obr.12625 (2018).

100. Buckley, H. R. & Tayles, N. G. The functional cost of tertiary yaws (*Treponema pertenue*) in a prehistoric Pacific Island skeletal sample. *J. Archaeol. Sci.* **30**, 1301–1314. https://doi.org/10.1016/S0305-4403(03)00020-7 (2003).

101. Weitzmann, M. N. The role of inflammatory cytokines, the RANKL/OPG axis, and the immunoskeletal interface in physiological bone turnover and osteoporosis. *Scientifica* https://doi.org/10.1155/2013/125705 (2013).

102. Alswat, K. A. Gender disparities in osteoporosis. *J. Clin. Med. Res.* **9**, 382–387. https://doi.org/10.14740/jocmr2970w (2017).

103. Thomas, C. D. L., Feik, S. A. & Clement, J. G. Regional variation of intracortical porosity in the midshaft of the human femur: age and sex differences. *J. Anat.* **206**, 115–125. https://doi.org/10.1111/j.1469-7580.2005.00384.x (2005).

104. Bousson, V. et al. Distribution of intracortical porosity in human midfemoral cortex by age and gender. *J. Bone Miner. Res.* **16**, 1308–1317. https://doi.org/10.1359/jbmr.2001.16.7.1308 (2001).

105. Bjørnerem, Å. The clinical contribution of cortical porosity to fragility fractures. *BoneKey Rep.* https://doi.org/10.1038/bonekey.2016.77 (2016).

106. Cooper, D. M. et al. Age-dependent change in the 3D structure of cortical porosity at the human femoral midshaft. *Bone* **40**, 957–965. https://doi.org/10.1016/j.bone.2006.11.011 (2007).

107. Ruff, C. B. & Hayes, W. C. Sex differences in age-related remodeling of the femur and tibia. *J. Orthop. Res.* **6**, 886–896. https://doi.org/10.1002/jor.1100060613 (1988).

108. Ruff, C. B. & Hayes, W. C. Cross-sectional geometry of Pecos Pueblo femora and tibiae—a biomechanical investigation: II. Sex, age, and side differences. *Am. J. Phys. Anthropol.* **60**, 383–400. https://doi.org/10.1002/ajpa.1330600309 (1983).

109. Heaney, R. P. et al. Peak bone mass. *Osteoporos. Int.* **11**, 985–1009. https://doi.org/10.1007/s001980070020 (2000).

110. Fotu, K. F. et al. Outcome results for the M'a'lahi Youth Project, a Tongan community-based obesity prevention programme for adolescents. *Obes. Rev.* **12**, 41–50. https://doi.org/10.1111/j.1467-789X.2011.00923.x (2011).

111. Chien, K. R. & Karsenty, G. Longevity and lineages: toward the integrative biology of degenerative diseases in heart, muscle, and bone. *Cell* **120**, 533–544. https://doi.org/10.1016/j.cell.2005.02.006 (2005).

112. Appleby, J. E. Why we need an archaeology of old age, and a suggested approach. *Nor. Archaeol. Rev.* **43**, 145–168. https://doi.org/10.1080/00293652.2010.531582 (2010).

113. Crews, D. E. Biological anthropology and human aging: some current directions in aging research. *Annu. Rev. Anthropol.* **22**, 395–423 (1993).

114. Zebaze, R. M. et al. Intracortical remodelling and porosity in the distal radius and post-mortem femurs of women: a cross-sectional study. *The Lancet* **375**, 1729–1736. https://doi.org/10.1016/S0140-6736(10)60320-0 (2010).

115. Clark, G., Burley, D. & Murray, T. Monumentality and the development of the Tongan maritime chiefdom. *Antiquity* **82**, 994–1008. https://doi.org/10.1017/S0003598X00097738 (2008).

116. Buikstra, J. E., & Ubelaker, D. H. *Standards for Data Collection from Human Skeletal Remains* (Arkansas Archaeological Survey Research Series, 1994).

117. Bruzek, J. A method for visual determination of sex, using the human hip bone. *Am. J. Phys. Anthropol.* **117**, 157–168. https://doi.org/10.1002/ajpa.10012 (2002).

118. Schmitt, A. Une nouvelle méthode pour estimer l'âge au décès des adultes à partir de la surface sacro-pelvienne iliaque. *Bull. Mem. Soc. Anthropol. Paris.* **17**, 89–101 (2005).

119. Sofaer, J. R. 10. Gender, bioarchaeology and human ontogeny. In *The Social Archaeology of Funerary Remains* (eds Gowland, E. & Knüsel, C.) 155–167 (Alden Press, Oxford, 2006).

120. Clark, G., Grono, E., Ussher, E. & Reepmeyer, C. Early settlement and subsistence on Tongatapu, Kingdom of Tonga: Insights from a 2700–2650 cal BP midden deposit. *J. Archaeol. Sci. Rep.* **3**, 513–524. https://doi.org/10.1016/j.jasrep.2015.08.005 (2015).

121. Miszkiewicz, J. J. & Mahoney, P. Ancient human bone microstructure in medieval England: comparisons between two socio-economic groups. *Anat. Rec.* **299**, 42–59. https://doi.org/10.1002/ar.23285 (2016).

122. Kaczmarek, E., Gorna, A. & Majewski, P. Techniques of image analysis for quantitative immunohistochemistry. *Rozk. Akad. Med. Bialymst. Supl.* **49**, 155–158 (2004).

123. Ashique, A. M. et al. Lacunar-canalicular network in femoral cortical bone is reduced in aged women and is predominantly due to a loss of canalicular porosity. *Bone Rep.* **7**, 9–16. https://doi.org/10.1016/j.bonr.2017.06.002 (2017).

124. Hunter, R. L. & Agnew, A. M. Intraskeletal variation in human cortical osteocyte lacunar density: implications for bone quality assessment. *Bone Rep.* **5**, 252–261. https://doi.org/10.1016/j.bonr.2016.09.002 (2016).

125. Bell, K. L. et al. Regional differences in cortical porosity in the fractured femoral neck. *Bone* **24**(57–64), 1999. https://doi.org/10.1016/S8756-3282(98)00143-457-64 (1999).

126. Thompson, D. D., Salter, E. M. & Laughlin, W. S. Bone core analysis of Baffin Island skeletons. *Arctic Anthropol.* **18**, 87–96 (1981).

127. Dempster, D. W. et al. Standardized nomenclature, symbols, and units for bone histomorphometry: a 2012 update of the report of the ASBMR Histomorphometry Nomenclature Committee. *J. Bone Miner. Res.* **28**(1), 2–17. https://doi.org/10.1002/jbmr.1805 (2013).

128. Lassen, N. E. et al. Coupling of bone resorption and formation in real time: new knowledge gained from human Haversian BMUs. *J. Bone Miner. Res.* **32**, 1395–1405. https://doi.org/10.1002/jbmr.3091 (2017).

129. Ren, F., Smith, I. O., Baumann, M. J. & Case, E. D. Three-dimensional microstructural characterization of porous hydroxyapatite using confocal laser scanning microscopy. *Int. J. Appl. Ceram. Technol.* **2**(200–211), 2005. https://doi.org/10.1111/j.1744-7402.2005.02022.x (2005).

130. Heaney, R. P., Recker, R. R., Watson, P. & Lappe, J. M. Phosphate and carbonate salts of calcium support robust bone building in osteoporosis. *Am. J. Clin. Nutr.* **92**, 101–105. https://doi.org/10.3945/ajcn.2009.29085 (2010).

131. Bertrand, L. et al. Development and trends in synchrotron studies of ancient and historical materials. *Phys. Rep.* **519**, 51–96. https://doi.org/10.1016/j.physrep.2012.03.003 (2012).

132. Ou-Yang, H., Paschalis, E. P., Mayo, W. E., Boskey, A. L. & Mendelsohn, R. Infrared microscopic imaging of bone: spatial distribution of CO32-. *J. Bone Miner Res.* **16**, 893–900. https://doi.org/10.1359/jbmr.2001.16.5.893 (2001).
133. Madupalli, H., Pavan, B. & Tecklenburg, M. M. Carbonate substitution in the mineral component of bone: Discriminating the structural changes, simultaneously imposed by carbonate in A and B sites of apatite. *J. Solid State Chem.* https://doi.org/10.1016/j.jssc.2017.07.02527-35 (2017).

134. Penido, M. G. M. & Alon, U. S. Phosphate homeostasis and its role in bone health. *Pediatr. Nephrol.* **27**(2039–2048), 2012. https://doi.org/10.1007/s00467-012-2175-z2039-2048 (2012).

135. Doi, Y., Iwanaga, H., Shibutani, T., Moriwaki, Y. & Iwayama, Y. Osteoclastic responses to various calcium phosphates in cell cultures. *J. Biomed. Mater. Res.* **47**, 424–433. https://doi.org/10.1002/(SICI)1097-4636(19991205)47:3%3c424::AID-JBM19%3e3.0.CO;2-0 (1999).

136. Paschalis, E. P. et al. FTIR microspectroscopic analysis of human osteonal bone. *Calcif. Tissue Int.* **59**, 480–487. https://doi.org/10.1007/BF00369214 (1996).

137. Fleet, M. E. Infrared spectra of carbonate apatites: v₂-Region bands. *Biomaterials* **30**, 1473–1481. https://doi.org/10.1016/j.bioma-terials.2008.12.007 (2009).

138. Nachar, N. The Mann-Whitney U: A test for assessing whether two independent samples come from the same distribution. *Tutor. Quant. Methods Psychol.* **4**, 13–20. https://doi.org/10.20982/ttmp.04.1.p013 (2008).

139. Taylor, R. Interpretation of the correlation coefficient: a basic review. *J. Diagn. Med. Sonogr.* **6**, 35–39. https://doi.org/10.1177/875647939000600106 (1990).

## Acknowledgments

The Lapaha Community and Nobles (His Serene Highness Prince Kalaniuvalu Fotofili, Her Serene Highness Princess Marcella Taumoepeau Tupou Kalaniuvalu Fotofili, and Her Royal Highness Mele Siu'ilikutapu Kalaniuvalu Fotofili) and the Ministry of Internal Affairs (Government of Tonga) permitted the excavation and analysis of the skeletal material. This study was undertaken as part of projects funded by the Australian Research Council (ARC): "Ancient human bone metabolism change with lifestyle in Asia-Pacific populations" (DE190100068 to Miszkiewicz) and "The First Polynesians: Their Origins, Lifeways and Environmental Challenges" (DP200102872 to Valentin and Clark). Earlier associated archaeological excavations were supported by the MAEDI and CNRS funding allocated to Valentin and an ARC grant to Clark (FT0990591). The Research School of Humanities and the Arts at the Australian National University (ANU) funded stay on-site at the ANSTO - Australian Synchrotron facility (Clayton, Victoria), where the sFTIRM data were collected using the IRM beamline. The ANU College of Arts and Social Sciences, and the ANU Major Equipment Grant (18MEC26 to Miszkiewicz) funded microscopy equipment. David McGregor offered continuous technical assistance on the project.

## Author contributions

J.J.M. led the project, secured funding, carried out lab work, data analysis, wrote the draft of the manuscript; C.V., N.A.S., M.J.T., J.V. were part of the synchrotron lab work and data analysis; F.V., G.C. co-designed the study, excavated specimens, collected data in the field, secured funding. All edited the manuscript and gave approval for publication.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-021-83264-3.

**Correspondence** and requests for materials should be addressed to J.J.M.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_12.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021