Substantiating microCT for diagnosing bioerosion in archaeological bone using a new Virtual Histological Index (VHI)

Recent technological advances have broadened the application of palaeoradiology for non-destructive investigation of ancient remains. X-ray microtomography (microCT) in particular is increasingly used as an alternative to histological bone sections for interpreting pathological alterations, trauma, microstructure, and, more recently, bioerosion with direct or ancillary use of histological indices. However, no systematic attempt has been made to confirm the reliability of microCT for histotaphonomic analysis of archaeological bone. The objective of this study is therefore to compare thin sections of human femora rated with the Oxford Histological Index to microCT sections using the newly developed Virtual Histological Index in order to provide an accessible methodology for the evaluation of bioerosion in archaeological bone. We provide detailed descriptions of virtual sections and assess the efficacy of the method on cranial and postcranial elements, cremated long bones, and faunal samples. The traditional histological and virtual methods showed a strong correlation, providing the first systematic data substantiating lab-based microCT as a suitable alternative tool for reconstructing post-mortem history in the archaeological record, and for the reliable, non-destructive screening of samples for further analyses.


Introduction
As human remains offer unique perspectives on mortuary behaviour in the past, the continued development and testing of bioarchaeological methods and techniques play an important role in refining our understanding of long-term taphonomic processes. In recent years, numerous histotaphonomic studies of archaeological human and faunal bone have looked at bioerosion to understand the mechanisms by which a corpse turns into a skeleton (e.g. Smith et al. 2002;Turner-Walker and Peacock 2008;Fernández-Jalvo et al. 2010;Hollund et al. 2012Hollund et al. , 2018Booth 2016Booth , 2017Brönnimann et al. 2018). Although there is disagreement concerning the origin of bioeroding microorganisms and what histotaphonomic analysis can truly tell us about early taphonomic events (Parker Pearson et al. 2005;Kontopoulos et al. 2016;Turner-Walker 2019), such analyses have proven to be highly effective for assessing the presence and type of bioerosion in archaeological bone. To improve the accessibility of histotaphonomic research in archaeology, this paper tests the utility of lab-based microCT imaging and presents the Virtual Histological Index (VHI), which together provide a rapid and minimally destructive method for the reconstruction of the early post-mortem history of archaeological bone.
The restructuring of bone microanatomy by microorganisms such as fungi and bacteria is referred to as biogenic bone degradation, or bioerosion (Piepenbrink 1986;Turner-Walker 2012). Upon death, at the outset of the post-mortem interval, bone is beset by osteolytic microorganismal activity induced by the host's enteric microbiome and/or by microbiota in the burial environment, which affect its surface and microstructure, eroding the histological anatomy (Child 1995;Janaway 1996;Jans 2008;White and Booth 2014;Kontopoulos et al. 2016;Metcalf et al. 2016;Damann and Jans 2017;Brönnimann et al. 2018). This not only influences chemical and radiological analyses but is also-along with other factors such as the pH value of the soil-responsible for the fragile nature of archaeological skeletons (Millard 2001;Nielsen-Marsh et al. 2007). Differences in environmental and soil conditions may affect the proliferation of osteolytic bacteria or saprophagous organisms (Wedl 1865;Kendall et al. 2018).

Background
Studies in the fields of palaeontology (Higgs et al. 2011) and anthropology (Le Garff et al. 2017) applied microCT to bone to identify different types of diagenesis, and in combination with histological indices (Dal Sasso et al. 2014;Booth et al. 2016), while others have used scanning electron microscopy (SEM) (Turner-Walker and Syversen 2002) or synchrotron radiation microCT (SR-microCT) (Caruso et al. 2020, 2021). Dal Sasso et al. (2014 and Booth et al. (2016) were amongst the first to see the value of lab-based microCT scanning for the evaluation of diagenetic alteration to archaeological bone, stating that virtual cross-sections are comparable to thin section micrographs as rated via the Oxford Histological Index (OHI). However, there has been no systematic attempt to confirm this, nor the broader utility of this technology for investigating microscopic diagenetic alterations in comparison to the well-established method of histotaphonomy.
MicroCT with commercial lab-based systems can produce detailed volume images of intact objects with microscopic resolutions. This obviously requires samples small enough to fit into the tomography setup, and in most commercial systems, the largest sample that can be imaged with ca. 5 µm voxels is 5 to 10 mm in diameter and 100 to 200 mm in length. Thus, the lab-based imaging we describe applies to small bones, bone fragments, or bone samples previously prepared for sectioning. For objects with sufficient inherent X-ray contrast, such as bone, no preparation is required apart from stable mounting of the specimen.
Importantly, studies comparing microCT images with histology to evaluate skeletal trauma (Baier et al. 2019), pathological alterations (Rühli et al. 2007), bone morphometry (Uchiyama et al. 1997;Müller et al. 1998), and cortical bone microstructure (Particelli et al. 2012), identify strong correlations between conventional histological section methods and virtual assessments. Therefore, comparing bioerosion in thin sections with virtual sections is a crucial first step towards confirming the reliability of microCT for identifying microscopic taphonomic alterations to bone. The aim of this research is thus to systematically compare transverse thin sections rated with the OHI to microCT scan images rated with the newly developed VHI. In doing so, we contribute an accessible, minimally destructive investigative tool for reconstructing post-mortem history in the archaeological record, and for the reliable screening of samples for further analysis. Importantly, the method also permits assessment of the entire virtual sample, including individual virtual slices through the 3D image, chosen locations within the sample, or larger virtual sections composed of several virtual slices in any orientation.

Dataset and sampling
The dataset is based on 28 samples selected from five archaeological sites in the temperate European environment of Lower Austria dating from the Early Neolithic to Late Iron Age (Table 1). Studies have demonstrated intra-skeletal variation for diagenetic alteration, and that the cortex of long bones tends to best describe bioerosion (e.g. Jans et al. 2004;Dal Sasso et al. 2014;Booth 2017), which is likely due to the longitudinal alignment of osteons and collagen axes, which simplify the histology and direction of tunnelling. To facilitate comparison with previous research (e.g. Jans et al. 2004;White and Booth 2014;Booth 2017), samples 1-2 cm thick were preferentially extracted from the anterior aspect of the proximal femoral diaphysis (n = 15) inferior to the surgical neck. Rib Fig. 1 Light micrographs showing the described histological features. a Inhibited bioerosion (MFD) and exogenous mineral precipitation in canals in femoral sample S01. The blue and yellow visualize collagen birefringence. The micrograph was taken with a lambda plate and 10(2 ×) magnification. b Very good collagen preservation in femoral sample GÖ01. The micrograph was taken with cross polarised light and 10(1.5 ×) magnification. c Inhibited bioerosion around Haversian canals (H) with Wedl type 2 tunnelling (W2) identified by enlarged canaliculi in femoral sample GÖ01. The micrograph was taken with polarised light and 10(2 ×) magnification (n = 2), mandible (n = 1), skull (n = 2), and humerus (n = 1) fragments were also selected to evaluate the efficacy of the method on non-femoral samples (Supplementary Materials S1). Five cremated long bones and two faunal mammalian long bones were also sampled. Samples were extracted with a diamond wheel blade attached to a Dremel® 3000 electric drill. The histological investigation was conducted using the well-established OHI (Hedges et al. 1995;Millard 2001), while the virtual assessment was completed using the newly developed VHI. Hedges et al. (1995) note the subjectivity inherent in rating samples with the OHI but found general agreement between observers within a single rating unit. This approach is applied here where we performed inter-rater reliability (IRR) tests between the two OHIs, between the two VHIs, and between the averaged OHIs and the averaged VHIs using the Pearson correlation coefficient (PCC), joint-probability of agreement (JPA), and limits of agreement (LOA). Overall IRR analysis, also known as inter-rater agreement or inter-observer reliability, is the degree to which there is agreement amongst independent observers rating the same phenomenon. The PCC is the linear correlation between two sets of data where it measures covariances, resulting in values ranging between − 1 (highly uncorrelated) and + 1 (highly correlated). JPA is the percentage of times that raters agree, and is the least robust measurement of IRR as it does not account for agreement based on chance. To compensate for this weakness, LOA was applied. This approach calculates the difference between each pair of observations with the mean value acting as the bias value ± 2*standard deviation (Bland and Altman 1986). This type of analysis provides insight into how much random variation influences the ratings, as depicted in the Bland-Altman plots. Furthermore, to assess the true scope of variation of agreement, the JPA calculations were evaluated for an exact match (identical observations), an approximate match (≤ 0.5 difference), and a general match (≤ 1.0 difference). This approach captures the degree to which observers entirely agree (disregarding the principle of equifinality), and with applied thresholds for when observers somewhat and generally agree.

Statistics
All statistical analyses were performed in R (v. 4.0.2). All figures were formatted using GIMP software (v. 2.99.6), and the ScientiFig plugin (Aigouy and Mirouse 2013) for Fiji software (Schindelin et al. 2012).

Light microscopy (LM) and scanning electron microscopy (SEM)
Twenty-four samples were selected for embedding in a two-compound resin (Biodur® E1/E2) following a protocol for undecalcified bone (Schultz 2001). The blocks were sectioned in the transverse plane with a Leica SP 1600 saw-blade microtome (Leica Microsystems) to a thickness ranging from 30 to 60 μm, and then mounted to a glass slide with Histokitt (ROTI®) and covered with a glass cover slip. A Nikon Eclipse Ni light microscope was used with 4 × , 10 × , and 20 × magnifications under plane light, polarised light (XPL), and differential phase contrast light. Additionally, a lambda plate was used with XPL to visualize collagen. Micrographs were taken with a DS-Ri2 camera (Nikon) mounted to the microscope. All visualizations were conducted with NIA Elements BR Software (Nikon). Six samples were mounted on carbon planchets and examined with an IT 300 LAB6 (JEOL Ltd.) scanning electron microscope at low vacuum at 25.0 kV under high magnification. Images were taken in backscatter electron mode.

Oxford Histological Index (OHI)
The OHI (Table 2) was developed to determine the percentage of bone unaffected by bioerosion within a single cross-section using thick, polished sections as evaluated through reflected light microscopy by the application of an ordinal scale from zero, which indicates poor preservation without original features save for Haversian canals, to five, which indicates good preservation, Table 1 Provenance and number of 28 cortical bone samples by period: femora (n = 15 + 5 cremated), ribs (n = 2), parietal (n = 2), humerus (n = 1), mandible (n = 1), faunal long bones (n = 2). For further details, see Supplementary Materials S1 like unaltered bone ( Fig. 2) (Hedges et al. 1995;Millard 2001). Jans et al. (2002) made two modifications to the OHI to provide a more detailed examination by (1) using thin sections, and (2) dividing the sample transversely, according to anatomy, into periosteal, endosteal, and midline sections, and then assigning each an OHI (Jans 2008). In the present study, Jans et al.'s (2002) modifications to the OHI were used to quantify bioerosion in samples via LM and SEM. Two observers separately assessed the thin sections using the OHI (Supplementary Materials S2).  (Hedges et al., 1995;Millard, 2001)

X-ray microtomography (microCT)
The same samples from which thin sections were taken for LM and SEM were placed in plastic cylinders stabilized with synthetic foam prior to microCT scanning. Scans were made at the University of Vienna Theoretical Biology imaging lab using a Zeiss/Xradia MicroXCT-200 or a Bruker/SkyScan 1272 microCT system, with image resolutions ranging from 4 to 10.9 μm (isotropic voxel size), sufficient to image osteolytic bacteria-induced non-Wedl tunnelling. Further scanning specifications and parameters are outlined in Supplementary Materials S3.

Results
Eight of the twenty-eight samples were selected to report what can be visualized in virtual cross-sections in comparison to what was identified in the thin section micrographs as an introduction to how bioerosion presents at each level of the VHI. This subset consists of samples from the femur and mandible, and includes one faunal sample, a cremated sample, and a severely osteoporotic sample. Selected images are included in the text while the remaining ones can be found in the Supplementary Materials (S4-S28).

SB01/mandibular ramus (OHI 5/VHI 5) (Fig. 3)
OHI with LM No MFD can be found in the entire thin section. There is reduced collagen content along the periosteal surface, but overall, collagen concentration is high.

GE15/femur (OHI 4/VHI 4.5) (Fig. 4)
OHI with LM The bacterial attack is arrested. The highest concentration of MFD is within the midline and endosteum; the trabeculae are affected and the periosteum mildly so.
Microfissures follow the microanatomy along the periosteal surface. Collagen preservation is very good throughout the entire section. There is exogenous mineral precipitation within the Haversian systems, primarily along the endosteal surface with a milder manifestation along the midline.
Virtual cross-sections Bioerosion is minimal. Grey values are mostly consistent throughout and between virtual slices with minor density loss (darker grey values), particularly around the midline and towards the endosteal surface. Bioerosion is concentrated around Haversian systems with numerous MFDs present as small clusters of perforations

S01/femur (OHI 2.75/VHI 3) (Fig. 5)
OHI with LM Bioerosion is inhibited. MFD is equally dispersed throughout the entire section. Collagen preservation is mediocre though it is better within the unaffected lamellae. There are many enlarged canaliculi and lacunae, evidence of Wedl type 2 tunnelling. Most of the canals are filled with exogenous mineral precipitation. The endosteal surface is not eroded, though the periosteal surface is damaged.
Virtual cross-sections A moderate level of bioerosion is present. Canals can be distinguished throughout much of the bone; when not visible, it is not due to bioerosion; the sample is an early adult and much of the endosteal lamellae remain where at this age, fewer secondary osteons have formed with cement lines. The periosteal lamellae are more heavily eroded. There are MFDs and patches of reduced density within the circumferential lamellae of the endosteal surface that connect with bioerosion surrounding Haversian systems. Lamellae of the Haversian canals are distinct within areas of relatively unaffected bone; the bioerosion, which presents as both MFD and patches of dark grey values, has primarily attacked the periosteum and the midline towards the periosteum. A large embedding artefact traverses the sample from the endosteum to the midline. Small patches of high density (white) material fill a few canals towards the endosteal surface and along the midline; the most external lamellae of the endosteum are also a higher density.
R03/femur (OHI 1/VHI 1.75) (Fig. 6) OHI with LM Bioerosion is present throughout the entire section. MFD is strongest along the midline towards the endosteum. A large area along the midline towards the endosteum remains mostly unaltered with unaffected lamellae still visible, though there are MFDs around the Haversian canals. Within this area, collagen content is very good and the lacunae remain a normal size. All canals and microfissures are filled with calcite. A large microfissure runs from the periosteum of the entire section through the histological anatomy, which may be the result of shrinkage as the bone dries in areas where there is no tunnelling. Microfissures following the microanatomy travel from the periosteum into the midline. Sediment is attached to the periosteal surface. GÖ03/femur (OHI 0/VHI 0) (Fig. 7) OHI with SEM The bacterial attack is extensive. MFDs are identifiable along the entire section with only Haversian canals remaining. A very small area of well-preserved bone remains along the endosteum; however, this surface is heavily eroded. Most of the canals are filled with an exogenous matrix, possibly bone fragments, soil, and/or calcite. The periosteum is also heavily eroded. Virtual cross-sections The sample has undergone extensive bioerosion with little unaltered bone remaining. Canals remain distinct but bioerosion is extensive throughout and between virtual slices; however, the bone is also porous and demineralized in part due to an osteoporosis-induced reduction in bone mineral density (Marcus and Bouxsein, 2010

IN01/cremated femur (OHI 4.75/VHI 4.5) (Fig. 9)
OHI with LM Macroscopic examination revealed the sample was cremated as confirmed by histology. It appears that there is no bioerosion, though discoloration makes the assessment difficult. Osteocyte lacunae are numerous and slightly enlarged, which is a common feature of cremated bone.
Howship's lacunae are also visible (Supplementary Material S7). There is exogenous mineral precipitation within the Haversian canals. The collagen signal is weak (Supplementary Material S7). Overall, the classic characteristics of concentric lamellae are lost as the matrix becomes homogeneous due to cremation (Hunger and Leopold 1978;Schultz 1986).

Virtual cross-sections
The consistently medium-dark grey grey-values of each virtual slice within the virtual stack (save for high density (bright white) exogenous mineral precipitation) confirm cremation; however, the evaluation was hindered by the inability to visualize certain anatomical features due to burning. The periosteal lamellae are missing, though some remain along the endosteum. Canals can be distinguished throughout, though they appear smaller than usual by visual inspection. However, a morphometric assessment is required to verify this. It is demonstrated that burning at medium-to-high/high temperatures results in shrinkage, cracking, and crystallisation (Hanson and Cain 2007;Boschin et al. 2015;Ellingham and Sandholzer 2020). Several long, narrow microfissures travel from both the periosteum and endosteum towards the midline; all are dense (whitish), which may be calcite. The visually undifferentiated lamellae surrounding Haversian canals are obscured by burning but are more clearly visible towards the endosteum and along the midline than the periosteum. That they remain visible may be due to the slightly denser cement lines that are visible around some Haversian systems.  . 10) OHI with LM This bone was macroscopically identified as non-human. There is a clear pattern of bioerosion with MFD visible in patches throughout the entire section; however, bioerosion is inhibited. Under polarised light, a thick, highdensity band can be seen along the periosteum and a thinner band along the endosteum. There is brownish staining along both surfaces. Haversian canals are visible in the histology, and the canaliculi are slightly enlarged (seen at higher magnifications, Supplementary Material S8). In areas without MFD, the lacunae are typical in size. Collagen preservation is poor even in the islands of bone that are free of MFD.

Virtual cross-sections
The bone has undergone moderate bioerosion. Some circumferential lamellae, individually visually undifferentiated, remain along the endosteum and periosteum; both surfaces have light grey, higher density bands as noted in the histology. Bioerosion manifests in patches relatively evenly throughout and between virtual slices, with a slightly higher concentration along the midline. Although MFDs have spread throughout the sample, they remain distinct and do not generally present as undifferentiated dark grey patches. There are several shrinkage artefacts that primarily travel longitudinally from the periosteum. There are few canals, as expected for plexiform animal bone, and the few Haversian lamellae that are visible are difficult to distinguish. There is potentially also osteon banding between two longitudinal microfissures below the transverse crack along the midline.

Statistical results
The following results (Table 3) show that exact interobserver matching within each histological index is at least 54% but diminishes across all observations. However, it is useful to also assess approximate and general observational matches, as the histological indices are categorical across several criteria. This approach elucidates the rate of approximate agreement of minimal variation and minor variation.
Assessing for these minimal disagreements (≤ 0.5) demonstrates a significant increase in observation agreement for both OHI and VHI ratings. Furthermore, a general matching with a threshold of ≤ 1.0 between observations demonstrates a significant increase of all categories, with a minimum of 93% agreement. The results of our PCC analyses (Table 4) demonstrate a strong positive (0.94) correlation between the OHI as rated by two observers, and a stronger positive (0.96) correlation between the VHI as rated by two observers, with a strong correlation overall (0.97) between the mean values of OHI and VHI observations. To account for observational agreement by chance, LOA analysis was performed. As histological indices are based on multiple criteria with multiple descriptive indicators, the problem of equifinality arises. This approach helps address the chance to which observers agree and assesses agreement based on genuinely similar observations using the same rating scale. The LOA analysis (Table 5) demonstrates that for the OHI, observer 1  Table 3 Joint-probability of agreement calculations between OHI observations, VHI observations, and across all observations. Each set of calculations was analysed for an exact match, approximate match (within ≤ 0.5), and general match (within ≤ 1.0) Joint-probability of agreement (n = 28)

Exact match
Approx. match (≤ 0.5) measures slightly biased (0.25), whereas for VHI observer 2 measures negligibly biased (− 0.02). These negligible differences indicate that both observers produce relatively similar results, with a trend of observer 1 rating marginally higher observations (4 s and 5 s) and observer 2 rating slightly lower observations (1 s and 2 s) (Fig. 11a, b).
A 95% confidence interval is used to determine the limits of agreement across observations (see Bland and Altman 1986). For the OHI, the systematic difference is greater than the VHI observations (Fig. 11a, b) as indicated by the upper and lower limit ranges. One observation (GÖ25) falls outside the LOA for the OHI assessment (Fig. 12a); however, this is the only difference where the observers rated a specimen significantly different. Similarly, VHI observations were consistently matching, and observer bias is distributed evenly with 1 observation falling outside the limits of agreement: GE08 where observer 1 rated the sample higher than observer 2 (Fig. 12b). Furthermore, the trend across observations is similar for both the OHI and VHI where observer 1 tends to rate more specimens higher and observer 2 tends to rate more specimens lower (Fig. 12c). However, the OHI illustrates a

Fig. 11
Frequency of a OHI ratings per observer and b VHI ratings per observer slightly more accentuated slope (y = 0.0639x + 0.0753) compared to the VHI (y = 0.0537x − 0.1671), demonstrating the minor differences in observations between the extreme ratings. This further shows that the observers are more prone to rate at either end of the scale and with minimal inherent bias in one observer compared to the other.

Discussion
The IRR results demonstrate that the VHI evaluation of virtual sections is congruent with the OHI rating of histological thin sections. As noted, certain structures taken into consideration by the OHI, such as individual concentric lamellae, cannot be differentiated, nor can structures smaller than 1 μm (e.g. osteocyte lacunae) be visualized without synchrotron-based microCT (Andronowski et al. 2017a,b), but undifferentiated concentric lamellae, and features such as shallow pits and tunnels can be seen with lab-based microCT. Moreover, our inability to visualize osteocyte lacunae may result from a contrast-to-resolution issue during the scan. Thus, the VHI is able to describe the level of bioerosion that a sample has undergone, specifically that which results from bacterial attack. It may be possible that non-MFD bioerosion can be detected and differentiated with microCT images; however, the samples included in this study almost exclusively exhibited MFD, and not type 1 Wedl cyanobacterial-induced tunnelling, which is primarily found in aquatic environments, nor type 2. Importantly, in other respects the virtual images are also consistent with what has been qualitatively described in conventional thin sections. Jans et al. (2002) specify five categories of diagenetic alterations as visualized via histology: (1) presence and type of MFD;

Qualitative comparisons with histology
(2) presence of inclusions (e.g. sand, fungi, framboids); (3) presence of infiltrations, such as stains exogenous to the bone matrix; (4) presence of microfissures; and (5) birefringence intensity (see Fig. 1b). We found that MFD can be clearly visualized in virtual sections; these can, moreover, be virtually measured (as can anatomical features) to refine their identification by size (Fig. 13). Such measurements are possible due to the pixel-size calibration inherent in microCT scans (Withers et al. 2021). In badly affected samples, it is difficult to ascertain which type of tunnelling is present, as the bioerosion presents as large dark grey Fig. 12 Bland-Altman plot. a Differences between OHI observations. b Differences between VHI observations. c Differences between all observations patches. However, this problem is also encountered with thin section micrographs of badly preserved samples, and, moreover, though morphologically distinct, different MFDs are very likely, "all aspects of the same type of bacterial attack, differing only in tissue microarchitecture and local hydrology" (Turner-Walker 2019, p. 35). We were also unable to visualize the enlarged or stained canaliculi diagnosed via histology. This is likely due to the small size of the tunnels, which, although they may be long, are not thick enough to be visualized at the resolution we used. We agree with Booth et al. (2016) that although inclusions can be observed based on differences in density (grey levels) and texture, no further identification of their composition can be made with benchtop microCT. Infiltrations such as stains are more difficult to distinguish in microCT images. While it may be possible to detect a stain after having been recognized through histology, it is not immediately apparent using virtual sections alone that a change in grey values results from, for example, the impregnation of inorganic substances as manifested macro-or microscopically as a stain, or if changes to the grey values result from demineralization caused by bacterial attack. Furthermore, while dissolution and recrystallization of HAp can be visualized as changes in grey values, it is not possible to determine the uptake or exchange of specific materials like uranium or fluorine, or nitrogen content, nor can microCT visualize the hypermineralized cuffs surrounding tunnels that result from mineral reprecipitation as is specifically achievable with SEM (Turner-Walker and Syversen 2002; Kendall et al. 2018). Microfissures recognized in histological sections can be clearly detected in virtual sections; it is also possible to determine if the microfissures are postdepositional, embedding artefacts, or result from excavation damage. Optical birefringence intensity cannot be assessed, though collagen content can be visualized with staining procedures not employed in this study (Handschuh et al. 2017). Lastly, virtual evaluation of sub-micron spongiform porosity (0.1-1 μm in diameter), a common form of tunnelling in archaeological bone Turner-Walker 2012), requires SR-microCT (see Caruso et al. 2020Caruso et al. , 2021 or synchrotron-based microscopy (Biswas et al., 2022).
As concerns qualitative descriptions of how bacteria invade bone, we once again find congruence between conventional microscopy and virtual imaging. Yoshino et al. (1991) describe regions heavily affected by bacteria and fungi as potentially corresponding to low X-ray density. As noted, areas of low density (dark grey values) that result from mineral dissolution and redistribution (Hackett 1981) correspond to MFD in the virtual sections. Yoshino et al. (1991) also report 5 to 10 μm in diameter vacuoles surrounding circumferential and Haversian lamellae, while Jackes et al. (2001) describe the orientation of MFDs along circumferential lamellae and the axis of Haversian systems. These vacuoles, or MFDs, were also visible in virtual sections within and surrounding the lamellae of Haversian systems and in circumferential lamellae. It has also been suggested that bacterial attack and chemical demineralization are impeded by cement lines Turner-Walker 2012;Kendall et al. 2018). We may be visualizing this in samples such as GE08 (Supplementary Materials S9) where bioerosion appears to specifically attack the undifferentiated concentric lamellae, while in some regions of S01 (see Fig. 5) bioerosion seems to affect the interstitial lamellae, surrounding but not penetrating the concentric lamellae. However, there were no observable cement lines in sample S01. Moreover, after a certain point, when bioerosion has become severe, even if only within patches as with SB02 (Supplementary Materials S10), it is not possible to distinguish where bacterial attack begins or ends. Further research is thus required to confirm that bacterial attack is curbed by cement lines.

Effects of sampling
According to both Dal Sasso et al. (2014) and Booth et al. (2016), individual slices within a stack are representative of the sample, or entire femoral midshaft, and thus, the OHI applied to a single transverse cross-section may be "considered representative of the sample itself" (Dal Sasso et al. 2014, p. 37). We assessed the complete stack of virtual sections from microCT images and confirmed that within a given stack individual slices are generally homogeneously affected by bioerosion. While it appears probable that this holds true for the entire midshaft, further research should assess different regions of the femoral diaphysis to confirm homogeneity within the entire shaft, particularly when staining is present in patches. Using microCT, Dal Sasso et al. (2014) found that skull, rib, and femur fragments present with different levels of bioerosion with femoral samples being the best preserved, and further suggest that the poorer preservation of rib and skull bone microstructure results from the higher volume of trabecular bone and porosity associated with the latter, which can be more heavily affected by inclusions and MFD. In the current study, we had both femoral and rib samples from only two individuals (GÖ01/GÖ24 (S11/S12) and GÖ04/GÖ25 (S13/S14)), none of which was badly affected by bioerosion; thus, we can neither confirm nor refute their findings. Nevertheless, our research suggests that, although the femoral diaphysis, which has a thick cortical layer, is preferred for visualizing canal structure and Haversian systems, the VHI can also be applied to mandibular, cranial, humerus, and rib samples.
As concerns burnt remains, Lemmers et al. (2020) conducted experiments to ascertain the visibility of bioerosion in histological sections of bone that had been cremated, concluding that it can be identified in burnt remains, though further research should investigate this to confirm these findings. Thus, the classical histological analysis provides good results; however, further experiments are required assessing the effects of burning on bioeroded versus uneroded remains in microCT images to better understand how bioerosion may be visualized in cremated remains. However, given the high concordance between the VHIs and OHIs of cremated samples, as demonstrated by the IRR, we included them in the dataset.

Conclusion
A strong correlation between histotaphonomy and virtual methods provides the first steps towards substantiating microCT for assessing microtaphonomy in archaeological bone. Though the evaluation of thin sections via light microscopy or SEM remains the best means to identify the presence and type of chemical and microbial diagenesis (e.g. collagen leaching using lambda filters, Wedl MFDs, and hypermineralized cuffs) affecting a sample, visualization with these methods entails time-consuming, labour-intensive, multi-step processes during any step of which artefacts may be accrued. Numerous 30 to 55-μm-thick slices must be sectioned with a microtome from an embedded sample with small areas of each section assessed individually. Moreover, transverse histological thin sections permit visualization in only a single direction. Although SR-microCT (Caruso et al. 2020(Caruso et al. , 2021, which can provide superior resolution, has also been used, it is not as readily accessible as lab-based microCT scanners, the costs of which are within the reach of start-up grants. This study has demonstrated that microCT can be used as a complementary method to histology for the rapid, accessible evaluation of bioerosion in archaeological bone, and for the screening of samples for further, more destructive analyses. The VHI method, like the OHI, remains subjective; however, the strong IRR results suggest that this is not a major problem. Further research may amalgamate the quantitative methods of Turner-Walker and Syversen (2002) with virtual techniques to improve accuracy and reproducibility. We also recommend that future analyses explore the effects of (1) cremation on the visualization of bioerosion, which has rarely been examined histologically (but see Hanson and Cain, 2007;Lemmers et al., 2020;Végh et al., 2021) as the presence or absence of bacterial attack may affect indices of cremation temperature, (2) the visibility of Wedl type 2 and cyanobacterial tunnelling in virtual sections, and (3) the applicability of other histological indices to microCT scan images. Lastly, the use of non-destructive and minimally invasive techniques extends the range of potential specimens to include very rare or fragile samples, and should be considered, where possible, before more destructive methods are employed to mitigate the ethical ramifications of sampling.
Funding Open access funding provided by University of Vienna. This paper is part of a dissertation project and sponsored by the Vienna Doctoral School of Ecology and Evolution (VDSEE) 2021-Completion Grant of the University of Vienna awarded to K.M. Thin section preparation was funded by the DEEPDEAD project (Deploying the Dead: Artefacts and Human Bodies in Socio-Cultural Transformations) under the HERA Joint Research Programme "Uses of the Past" (UP) (Collaborative Research Project (CRP) no. 15.055) and the European Union's Horizon 2020 research and innovation programme under grant agreement No. 649307.
Data availability All original data are available in the online supplementary information.
Code availability Not applicable.

Declarations
Ethics approval Not applicable.

Conflict of interest The authors declare no competing interests.
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:// creat iveco mmons. org/ licen ses/ by/4. 0/.