Elemental composition of primary lamellar bone differs between parous and nulliparous rhesus macaque females

Paola Cerrito; Bin Hu; Justin Z. Goldstein; Rachel Kalisher; Shara E. Bailey; Timothy G. Bromage

doi:10.1371/journal.pone.0276866

Abstract

Extracting life history information from mineralized hard tissues of extant and extinct species is an ongoing challenge in evolutionary and conservation studies. Primary lamellar bone is a mineralized tissue with multidien periodicity that begins deposition prenatally and continues until adulthood albeit with concurrent resorption, thus maintaining a record spanning several years of an individual’s life. Here, we use field-emission scanning electron microscopy and energy-dispersive X-ray analysis to measure the relative concentrations of calcium, phosphorous, oxygen, magnesium and sodium in the femora of seven rhesus macaque with known medical and life-history information. We find that the concentration of these elements distinguishes parous from nulliparous females; that in females calcium and phosphorus are lower in bone formed during reproductive events; and that significant differences in relative magnesium concentration correlate with breastfeeding in infants.

Citation: Cerrito P, Hu B, Goldstein JZ, Kalisher R, Bailey SE, Bromage TG (2022) Elemental composition of primary lamellar bone differs between parous and nulliparous rhesus macaque females. PLoS ONE 17(11): e0276866. https://doi.org/10.1371/journal.pone.0276866

Editor: Inge Roggen, Universitair Kinderziekenhuis Koningin Fabiola: Hopital Universitaire des Enfants Reine Fabiola, BELGIUM

Received: March 14, 2022; Accepted: October 16, 2022; Published: November 1, 2022

Copyright: © 2022 Cerrito et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: The complete dataset (both normalized and in percentage) and the EDS spectra for all the segments are available on the Open Science Framework at the following URL: https://osf.io/6w24g/.

Funding: Support was provided by the NSF DDRIG (grant number: 2018357); by the Program in Early Cultures at Brown University; by the 2010 Max Planck Research Award to TGB in respect of the Hard Tissue Research Program in Human Paleobiomics. The Zeiss Gemini 300 FE-SEM was provided by the NIH S10 Shared Instrumentation Program (grant number: 1S10OD026989-01). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

The reconstruction of demographic patterns and life-history scheduling of both extant and extinct mammalian species is an ongoing challenge in several disciplines, including evolutionary biology, conservation studies and biological anthropology. The reconstruction of life-history scheduling has been routinely addressed via the histological analysis of mineralized tissues. These tissues act as recoding structures [1], keeping track of the physiological changes experienced by an organism in the form of accentuated growth lines and changes in elemental composition. A large body of work spanning several decades has shown that a variety of both intrinsic and extrinsic events, including birth [2, 3], reproduction [4], weaning [5], menopause [6], illnesses [7], dietary changes [8] and climate variation [9–11], leave a permanent record in calcified tissues which is preserved even in fossil specimens [12]. The microstructural organization of teeth is not altered during life [13]. As such, teeth have been the structure of choice for the identification and chronological reconstruction of past physiological changes in mammals. Fortunately, primary lamellar bone is not significantly remodeled in shorter-lived species such a rhesus macaques [14] and other catarrhine primates [15]. Large segments of this tissue are also maintained well into adulthood in humans [9], with some variation across geographic groups.

Lamellar bone is an incremental tissue formed according to a coordinated internal biological rhythm. It is the same rhythm visible in the Retzius Periodicity (RP) of enamel [16] and it remains constant throughout the entire life of the individual. Its recording sensitivity is high. In rhesus macaques one lamella has a thickness in the order of 3.5–4μm and a period of four days (rarely five, as per personal observation of TGB). While its period of registration spans prenatal life to adulthood, the duration of the record is limited due to resorption of earlier-formed bone.

The skeleton participates to whole organism physiology and reflects alterations in endogenous and exogenous conditions by undergoing changes in both growth rate [9] and elemental composition [17]. Research on humans [18] and rodents [19] has shown that several adaptations such as increased absorption, decreased renal excretion and increased skeletal resorption allow for the calcium (Ca) and phosphorous (P) needs of the forming fetuses to be met without necessitating a concurrent increase in maternal Ca and P intake. Especially during the last term of gestation, maternal serum Ca and P, but not magnesium (Mg) [17] or sodium (Na) [20] are decreased. This is remarkable considering that, in rodents, up to ~80% of the Ca present in maternal circulation is transferred, each hour, to the developing litter [17].

Elemental analysis employing laser ablation inductively coupled plasma mass spectrometry on skeletal remains associated with known ages at life-history milestones has shown that weaning is recorded in the teeth of walruses [5], humans [21] and macaques [22]. Similarly, Changes in zinc concentrations in teeth have been reported to reflect the onset of sexual maturity in walruses [23], whereas declining fertility has been shown to correlate with increased copper concentrations in rat teeth [24]. Other research on both primates and cetaceans, have employed scanning electron microscopy (SEM) equipped with energy dispersive X-ray spectroscopy (EDS) to assess the changes in major and minor element concentrations over time [25, 26]. However, the chemical signature of reproductive events in hard tissues is yet to be characterized in any mammalian species.

In this study we use field emission (FE) SEM-EDS to measure changes in five major and minor element concentrations (oxygen, calcium, magnesium, sodium and phosphorus) in the femora of seven male (n = 3) and female (n = 4) rhesus macaques with known medical and life-history information. Given the participation of the skeleton to the metabolic demands of reproduction, we predict the following:

The elemental composition across all bone of parous females will be significantly different from that of males and nulliparous females.
P and Ca concentrations will be lower in the tissues formed during pregnancy and lactation.
Maternal relative percent concentrations of Mg will be higher in the bone formed during pregnancy and lactation since other elements, such as P and Ca, will be transferred to a proportionally greater extent than Mg.
Relative percent concentrations of Mg in infants will be lower in the tissues formed during breastfeed than after weaning is completed, since other elements, such as P and Ca, will be transferred to a proportionally greater extent than Mg.

Materials and methods

Specimen information

The present work did not require IACUC approval as it neither used live subjects, nor did it require sacrificing any animals. All the individuals in this study died from causes unrelated to this research. The research adhered to the American Society of Primatologists’ Principles for the Ethical Treatment of Nonhuman Primates. The skeletonized specimens were acquired from the Caribbean Primate Research Center (CPRC). Prior to their death the animals were housed at the Sabana Seca field station of the University of Puerto Rico where they were provisioned with monkey chow. Veterinarians at the Field Station monitored and recorded information on their health and reproductive history. The IDs of the seven femora used in this study, together with each individuals’ sex, parity and dates of birth and death are reported in Table 1. The medical and life history information is reported in S1 Table. In order to calculate the growth rate of the lamellar bone [16], we also acquired the value of the RP from the individual dental sections of the same individuals [27]. This value (four days) was the same for all seven specimens.

Download:

Table 1. Specimen information.

https://doi.org/10.1371/journal.pone.0276866.t001

Block preparation and imaging

We selected the right femora of each of the seven individuals. We removed a ~1.0 cm block from the midshaft of each femur. Each block was cleaned, dehydrated using progressively higher-concentration solutions of ethanol, and then embedded in poly methyl methacrylate (PMMA) following procedures described in detail by Goldman and colleagues [28]. Both sides of the block were then ground to 1200 grit finish using progressively graded carbide papers and polished with a 0.1-μm diamond suspension. The blocks were mounted on an EXAKT (Oklahoma City, OK, USA) plastic slide using cyanoacrylate adhesive.

Each block was imaged in reflected circularly polarized light (CPL) using a Leica DMRXE microscope fitted with an automated X, Y and Z stage. We collected a montaged micrograph of each block at a resolution of 1.3 μm / pixel. Given the high resolution with which these micrographs were acquired, the identification of contiguous patches of lamellae was unambiguous. On each block, we used the reflected CPL montages to identify two regions fulfilling our criteria of contiguity across which we collected the data, from here onwards identified as Transect 1 and Transect 2. The CPL montages were used as references for the SEM imaging, which was performed before proceeding with the EDS data collection. Once we identified on the EDS the same field of view that we had chosen with the CPL images, we verified that the patches of lamellae were indeed contiguous, and only then proceeded to acquire EDS data.

Elemental data collection

To be suitable for high-vacuum SEM imaging and EDS analysis, we uniformly carbon-coated the blocks using a Denton Desk V (Denton Vacuum, Moorestown, NJ, USA). The thickness of the carbon coating is ~20nm. It is also worth noting that Monte Carlo simulations [29, 30] clearly demonstrate that the depth of penetration of the electron beam means that signal emanates from a thickness of bone that is about 2 orders of magnitude deeper than the carbon coating thickness.

All data were collected using a Bruker Quantax 200 XFlash 6160 EDS detector (Bruker, Billerica, MA, USA) coupled with a Zeiss Gemini-300 FE-SEM (Carl Zeiss Microscopy, White Plains, NY, USA). Before collecting the data, we confirmed by imaging in the SEM that the surface was free of any topographical relief. We acquired the data with the following SEM settings: in high vacuum mode, with an accelerating voltage of 15kV, with an electron beam aperture of 60μm, in high current mode and with a working distance of 8.5mm, which is the analytical working distance of this particular microscope and detector setup. We acquired the data in line scan mode and collected one line scan per segment. The spot size was set each time to match the radius of the interaction volume of the beam, which was consistently ~1.0μm since all the samples had similar composition and were analyzed with the same experimental conditions. The scan dwell time was set at 64μs, and the spectrometer throughput at 275Kcps. To ensure that we had at least 10,000 counts for each peak, the scan time was manually determined for each specimen since they varied in dimension (length of the segment). We used the zero kV peak to make sure that each spectrum was properly aligned. After acquisition was complete, we used the P/B-ZAF standardless analysis quantification model [31] in interactive mode which allowed us to manually validate each peak.

Each transect (two per specimen, 14 in total) was acquired in several separate segments (56 in total), moving from the periosteal surface towards the endosteal one. The spectra corresponding to each segment is available online as supplementary data (see the Data Availability Statement). Each segment was chosen to avoid osteocyte lacunae and dehydration cracks as best we could. Details regarding the segments collected for each transect are available in Table 2. The elements that were detected with sufficient counts per peak in significantly above-background concentrations in all 56 segments were: C, O, P, Ca, Na, Mg. We saved our data as.csv files in which for each point across the femoral transect we recorded the number of counts for each of the elements detected in above-threshold concentrations.

Download:

Table 2. Dimensional information regarding each segment and transect for all seven femora.

https://doi.org/10.1371/journal.pone.0276866.t002

Thin section preparation and imaging

The carbon-coated side of each block was mounted on an EXAKT (Oklahoma City, OK, USA) plastic slide using cyanoacrylate adhesive. Using a Buehler (Lake Bluff, IL, USA) Isomet 1000 saw we sectioned the block so that a ~150-μm thick section remained attached to the slide. Each section was then ground to 100 ± 10 μm and polished with a 0.1-μm diamond suspension. The resulting thin sections were imaged using the same Leica DMRXE microscope but in transmitted circularly polarized light (CPL). We acquired a montaged image of each section at a resolution of 0.325 μm / pixel. Since we mounted the carbon-coated surface directly on the plastic slide, without polishing or otherwise touching the surface, these micrographs clearly show the segments that were targeted by the electron beam of the SEM.

Chronological age assignation

To analyze the possible correlations between the elemental data and the differing physiological states of the animals we assigned a chronological age range to each transect (Table 2). First, we used the montages acquired in CPL mode to measure transect-specific growth rates of the bone. Along each transect of each individual we measured the thickness of the lamellae five times and calculated the average value (S2 Table). By knowing the RP values (four days), we were able to infer the transect-specific average yearly growth rate. For the four younger individuals (MA024, CK85, AV27 and BD24) the femoral diameter was still increasing at the time of death. Therefore, we aged the segments of primary lamellar bone from the periosteal margin towards the endosteal one, assuming that the age at the periosteal one corresponded to the age at death. For the older individuals M361, M526 and M591 primary lamellar growth had ceased before their death. We aged the segments of these individuals from the endosteal margin of the primary lamellar bone towards the periosteal one and derived the age of the endosteal one from corresponding diameters of the immature individuals (Fig 1 and S1–S7 Figs). The only exception is Transect 2 of M526 for which, due to extensive remodeling and cortical drift it was not possible to match the endosteal margin with corresponding primary bone of younger female individuals. Hence, this transect was aged from the periosteal towards the endosteal region by matching its periosteal margin to that of Transect 1 of the same individual. Primary lamellar bone growth ceased between seven and eight years for females, and up to 12 years for males, likely relating to their increase in body mass well into adulthood [32, 33].

Download:

Fig 1. Montages of micrographs acquired in reflected light of the blocks of the seven femora.

For each femur we report the sex and age of each. The micrographs are oriented with the posterior portion towards the bottom and the medial one towards the left.

https://doi.org/10.1371/journal.pone.0276866.g001

Data preparation

First, we discarded the data points collected over any cracks. Then, within each data file and corresponding to each acquisition point across the femoral transect, we added the following variables: age; category and specific type of event relating to reproduction and weaning; type of individual (male, nulliparous female, parous female).

Illnesses and traumas were detailed in the life histories, so we were able to track in time the progression of an affliction. For gestation, we began at the point in bone corresponding to an offspring’s date of birth and counted back 5.5 months (c. 165 days) prior [34–37]. Lactation periods may vary by individual, so for its coding we began at the point in bone corresponding with an offspring’s date of birth and counted ahead 5 months, as this was found to be the minimum duration of breastfeeding in various studies [34, 38, 39]. Ages at weaning were available for 4/7 macaques, and thus this period was averaged (9.5 months) and applied to the two individuals for which data was lacking but which experienced breastfeeding during the periods of cementum formation that we analyzed. The averaged period of 9.5 months falls in line with previous scholarship on macaque lactation duration [38]. All cementum deposited before the age at weaning was considered “suckling”. We use suckling to indicate the period from birth to first introduction of solid foods and weaning as the period between the first introduction of solid foods and the transition to al 100% non-milk diet.

In most cases, there was only one Specific Event for a given period, but for some individuals, two Specific Events may have been experienced simultaneously. These cases were mostly related to females who experienced illness while pregnant or lactating. In these scenarios, the reproductive specific event took precedent and thus we coded it accordingly.

Carbon coating our specimens artificially increased concentrations of C, rendering absolute elemental values problematic. Instead, we analyzed the data in a way that they would reveal any changes in elemental composition over time (and physiological events). To this end, we pre-processed our original data in two different ways, each one aimed at addressing a different question:

Within each data file we scaled the data for each element, which was originally recorded as net total (NETTO) counts, such that the minimum value would be 0 and the maximum would be 1. This type of processing provides information regarding the behavior of a single element over time. We use this data to address prediction 2. This dataset is referred to as “normalized data”.
Within each sample, we eliminated the values for C and then calculated the percent concentration for each element at each acquisition point. For example, at t₁%Mg = NETTO_Mg / (NETTO_Mg + NETTO_Na + NETTO_P + NETTO_Ca + NETTO_O)*100. This implies that the resulting relative percent concentrations reveal how changes in the concentration of one element are correlated with changes in the concentration of the other ones. This type of data allows us to analyze of the simultaneous behavior of several elements as their numerical relationship is maintained, even if all the values are artificially increased by the absence of C. We used these data to address predictions 1, 3 and 4. This dataset is referred to as “percentage data”.

Data analysis

All data analysis was carried out using R version 3.6.3 (R Core Team, 2020). For prediction 1, we used proprietary scripts (S1 File, functions pair.test.all, remove.ns and final) to identify individuals who were significantly different from all other individuals belonging to different groups, defined as either male, nulliparous or parous female. Additionally, we performed a principal component analysis (PCA) to visually appreciate the separation between the different types of individuals. For predictions 2, 3 and 4 we ran Bonferroni-corrected two-tailed pairwise t-tests to assess whether there were significant differences in elemental composition associated with either maternal reproductive events or infant breastfeeding. Primary lamellar bone stopped forming before M591 experienced her single reproductive event. Since the age range of primary lamellar bone for M591 (5.3 to 7.5 years) did not encompass the age at which this female gave birth (7.8 years), data from this individual was not included in the analyses for predictions 2 and 3, which test the elemental composition of bone formed during reproductive activity.

Results

Prediction 1

Our results support the prediction that the elemental composition of reproductive females is different from that of nulliparous females and young and old males. The average elemental percentage of Ca, O, P, Mg and Na, across both segments is significantly different (p < 0.05) for the two parous females (M591 and M526, S3 Table). This is readily appreciable in the PCA (Fig 2), where there is almost no overlap between nulliparous and parous females. We also conducted a post-hoc PCA to assess whether this difference is a function of the age of each data point but found that the difference between parous females and other individuals is present regardless of the specific chronological age of data points (S8 Fig).

Download:

Fig 2. PCA on percentage data grouped by individual.

Both segments from the same individual are pooled together. M526 and M591 are the only parous females, while the other adult female in the dataset (BD24) is nulliparous.

https://doi.org/10.1371/journal.pone.0276866.g002

Prediction 2

Our results fully support our prediction (Fig 3). Normalized Ca values are lower during reproductive events in Transect 1 (p = 0.025) and Transect 2 (p < 0.00001) of M526. Normalized P values are also lower during reproductive events in both Transect 1 and Transect 2 (p < 0.00001) of M526. Both Ca and P normalized values are significantly lower, in M526, during lactation than during gestation in Transect 2 (p < 0.00001), but not in Transect 1, where values are instead lower during gestation (p < 0.00001) (S10 and S11 Figs).

Download:

Fig 3.

Boxplot showing the difference in normalized concentrations of calcium (left) and phosphorous (right) in the two transects (Trans1 and Trans2) of the bone of M526 formed during reproductive events (gestation and lactation) or during times of no known physiologically challenging events (NKE). Calcium and phosphorus concentrations during reproductive events are significantly lower (p<0.05) in both transects.

https://doi.org/10.1371/journal.pone.0276866.g003

Prediction 3

We do not find conclusive support for prediction 3 (S12 Fig), as Mg relative percent concentration is lower in the bone formed during reproductive events only in Transect 2 (p < 0.00001), but not in Transect 1 (p = 0.88) of M526.

Prediction 4

Prediction 4 is fully supported by our results. Mg relative percent concentration is significantly lower during breastfeeding (p < 0.00001) in both transects of MA024 (Fig 4) and CK85 (S13 Fig). We find also that Mg relative percent concentration is significantly lower in prenatal cortical bone than in the postnatal cortical bone formed both during and after breastfeeding. This is true in both transects of MA024, which is the only individual for which prenatal bone is preserved (Fig 4).

Download:

Fig 4. Boxplot showing the difference in relative percent concentrations of magnesium in the bone of MA024 formed prenatally, during breastfeeding (suckling and weaning) and during times of no known physiologically challenging events (NKE).

The difference is significant in both transects (p<0.00001).

https://doi.org/10.1371/journal.pone.0276866.g004

Discussion

As per prediction 1, we find that there is a significant difference (S3 Table) in relative percent concentration of the five elements considered, between parous females and all other individuals, regardless of the age at which the bone was formed. It is important to notice that one of the parous females included in our sample (M591) had not yet experienced her single reproductive event during bone formation. Therefore, all the bone sampled for this female was formed while she was still nulliparous. One could therefore hypothesize that the difference detected between parous and nulliparous females is instead reflecting a difference between younger and older females, but this is not the case. The observed difference between parous and nulliparous females is independent of age (S8 and S9 Figs). Our results are in agreement with the recent work [40] showing that stable isotopes record sex-specific, but not age-specific changes in rhesus macaques. Work investigating the effects of parity on bone mineral content and density in rhesus macaques has also found that bone characteristics differ between age-matched parous and nulliparous females [41]. These combined findings indicate a retroactive effect of reproduction on bone elemental composition, as previously deposited bone likely contributes to meeting the elemental demands of the developing fetuses [42]. It is also possible that our records fail to report any unsuccessful pregnancy that M591 might have experienced before the first one reported at 7.8 years. Indeed, the mean age at first birth in rhesus macaques is 3.8 years [43], which is four years before the first reproductive event recorded for M591.

In agreement with the vast body of literature documenting the physiological response to the mineral needs of the developing fetus and infant [17–19, 44], bone formed during gestation and lactation has lower concentrations of both Ca and P.

The combined findings of prediction 1 and prediction 2 indicate that: overall Ca concentrations are lower in parous females compared to other individuals, and that they are even significantly lower in the bone of these females formed during reproductive events as bone resorption increases [45]. Conversely, we found no changes in the concentrations of O (S14 Fig).

In agreement with prediction 3 and partially with prediction 4, Mg relative percent concentrations are significantly different in the maternal bone formed during gestation and lactation, and in the bone formed while an infant is breastfeeding. Since the data pre-processing that we performed returns relative percent concentration of each element at each acquisition point across the femoral transects, our results indicate that during reproduction Mg might be lost by the mother proportionally less than the other elements, while concurrently received proportionally less than other elements by the breastfeeding infant. This is in agreement with research on humans showing that serum Mg is minimally altered during reproductive activity [46–48], while serum Ca and P decrease during the same time period [44]. Furthermore, in MA024, for which the endosteal-most portion of the transects is of prenatal formation, we observe an even lower concentration of Mg than during the breastfeeding (suckling and weaning) phases (Fig 4). This suggests a delay in Mg transfer from the maternal organism until after the birth of the fetus. This delay in Mg transfer has also been directly observed in the body fluids of dam and fetal rabbits [49].

Due to the margin of error likely associated with the chronological age assigned to each bone transect, we did not test for elemental signatures of more punctual events, such as physical trauma.

The findings regarding the elemental changes associated with reproduction are of relevance as the detection of parturitions from mineralized tissues is still a vastly unexplored area of research with significant implications for evolutionary, conservation and archaeological studies. However, given the small size of our sample and limited taxonomic breadth, further research on more animals, preferably from wild populations, is necessary to determine the scalability of these results. It would be important to use femora from wild populations as the present specimens were both provisioned and sheltered from extreme climatic events. It is therefore possible that the signal of reproductive events and weaning that we detected could be masked, in wild populations, by physiological responses to changing diets and environments. From an analytical perspective, studies on animals with lower lamellar growth rate would benefit from an increase in the temporal resolution of the data, while those on species with higher growth rate (e.g., humans) would likely have to employ sampling methods with higher spatial resolution (e.g. laser-ablation inductively coupled mass-spectrometry). Indeed, in our sample rhesus macaque primary lamellar bone growth rate is the range of 0.75 to 1.25 μm/day, while in humans a range of 0.63 to 0.88 μm/day has been reported [50]. Hence, at a given spatial resolution (e.g., 1μm, as in the present study) the temporal resolution for macaques is slightly higher than for humans. Finally, the present work has the caveat of not providing absolute quantitative values of elemental concentrations that can be associated with either reproduction or suckling. Our results are relative either to other animals (e.g., parous females vs. nulliparous ones), or to different periods of a same animal’s life (e.g. reproductive events vs. no physiological stressors) and therefore need reference data to be interpreted. Future work employing absolute quantification methods (such as laser ablation inductively coupled mass spectrometry) are necessary in order to provide absolute quantitative values corresponding to specific physiological events. Quantitative data, especially of elements present in low concentrations, can also be successfully collected at below micrometer spatial resolution using wavelength dispersive spectroscopy (WDX) [51]. EDX technique has the advantage of being much faster and allowing whole spectrum acquisition, however it does have high detection limits, which make it not suitable for trace element identification. Conversely, WDX is much more time consuming as the entire range of wavelengths must be scanned one at the time, but this allows for much lower detection limits. An additional advantage of WDX comparing to EDX is that it is suitable also for the quantification of irregularly shaped surfaces [52]. It therefore appears to be a trade-off between the two techniques, with EDX being more suitable for exploratory analysis of a wide range of elements, and WDX best suited for hypothesis-driven experiments targeting a small number of elements. The present work should therefore serve as foundational to future research employing different analytical methods to target elements present in low concentrations.

Supporting information

S1 Fig. Montaged micrograph of the femur of MA024 obtained in reflected light.

The red lines indicate the locations of transects 1 and 2. For each transect we report the length in micrometers; the yearly growth rate; the age range covered by the transect. For a description of the method used to derive the ages, see the Methods section of the main manuscript.

https://doi.org/10.1371/journal.pone.0276866.s001

(TIF)

S2 Fig. Montaged micrograph of the femur of M591 obtained in reflected light.

The red lines indicate the locations of transects 1 and 2. For each transect we report the length in micrometers; the yearly growth rate; the age range covered by the transect. For a description of the method used to derive the ages, see the Methods section of the main manuscript.

https://doi.org/10.1371/journal.pone.0276866.s002

(TIF)

S3 Fig. Montaged micrograph of the femur of M526 obtained in reflected light.

Montaged micrograph of the femur of M526 obtained in reflected light. The red lines indicate the locations of transects 1 and 2. For each transect we report the length in micrometers; the yearly growth rate; the age range covered by the transect. For a description of the method used to derive the ages, see the Methods section of the main manuscript.

https://doi.org/10.1371/journal.pone.0276866.s003

(TIF)

S4 Fig. Montaged micrograph of the femur of M361 obtained in reflected light.

The red lines indicate the locations of transects 1 and 2. For each transect we report the length in micrometers; the yearly growth rate; the age range covered by the transect. For a description of the method used to derive the ages, see the Methods section of the main manuscript.

https://doi.org/10.1371/journal.pone.0276866.s004

(TIF)

S5 Fig. Montaged micrograph of the femur of CK85 obtained in reflected light.

The red lines indicate the locations of transects 1 and 2. For each transect we report the length in micrometers; the yearly growth rate; the age range covered by the transect. For a description of the method used to derive the ages, see the Methods section of the main manuscript.

https://doi.org/10.1371/journal.pone.0276866.s005

(TIF)

S6 Fig. Montaged micrograph of the femur of BD24 obtained in reflected light.

The red lines indicate the locations of transects 1 and 2. For each transect we report the length in micrometers; the yearly growth rate; the age range covered by the transect. For a description of the method used to derive the ages, see the Methods section of the main manuscript.

https://doi.org/10.1371/journal.pone.0276866.s006

(TIF)

S7 Fig. The red lines indicate the locations of transects 1 and 2.

For each transect we report the length in micrometers; the yearly growth rate; the age range covered by the transect. For a description of the method used to derive the ages, see the Methods section of the main manuscript.

https://doi.org/10.1371/journal.pone.0276866.s007

(TIF)

S8 Fig. PCA of relative percentage data grouped by individual type (both segments from the same individual are pooled together) and colored according to the age associated with each datapoint of each individual.

As visible by the superimposition between blue and yellow datapoints, age is not a discriminating factor.

https://doi.org/10.1371/journal.pone.0276866.s008

(TIF)

S9 Fig. PCA of relative percentage data of only females, colored by parity state (blue = nulliparous, yellow = parous).

The data from both segments of the same individual are pooled together. There is clear separation between parous and nulliparous females. Dot size is proportional to age. It is clearly visible that the separation between parous and nulliparous females is independent of age (dot size).

https://doi.org/10.1371/journal.pone.0276866.s009

(TIF)

S10 Fig. Boxplot showing the difference in relative percent concentrations of calcium in the bone of M526 formed during gestation and during lactation.

Ca values are significantly lower only in Transect 2 (p<0.00001).

https://doi.org/10.1371/journal.pone.0276866.s010

(TIF)

S11 Fig. Boxplot showing the difference in relative percent concentrations of phosphorous in the bone of M526 formed during gestation and during lactation.

P values are significantly lower only in Transect 2 (p<0.00001).

https://doi.org/10.1371/journal.pone.0276866.s011

(TIF)

S12 Fig. Boxplot showing the difference in relative percent concentrations of magnesium in the bone of M526 formed during reproduction (gestation and lactation) or during times of no known physiologically challenging events (NKE).

There is a significant difference only in Transect 2 (p<0.001).

https://doi.org/10.1371/journal.pone.0276866.s012

(TIF)

S13 Fig. Boxplot showing the difference in relative percent concentrations of magnesium in the bone of CK85 formed during breastfeeding (suckling and weaning) or during times of no known physiologically challenging events (NKE).

The difference is significant in both transect 1 and transect 2 (p<0.00001).

https://doi.org/10.1371/journal.pone.0276866.s013

(TIF)

S14 Fig. Boxplot showing the difference in relative percent concentrations of oxygen in the bone of M526 formed during reproduction (gestation and lactation) or during times of no known physiologically challenging events (NKE).

https://doi.org/10.1371/journal.pone.0276866.s014

(TIF)

S1 Table. List of age at known impactful physiological events for each individual.

Ages are reported in years. The duration of each medical event was inferred from the duration of the medical treatment received by the animal, while the duration of gestation was assumed to be the same for all individuals and taken from the published literature, as reported in the Supplementary Methods section.

https://doi.org/10.1371/journal.pone.0276866.s015

(XLSX)

S2 Table. Detailed measurements for each transect.

For each transect: five measures of lamella thickness and corresponding yearly growth rate; average (avg) value of lamella thickness and corresponding growth rate, together with the standard deviation (sd).

https://doi.org/10.1371/journal.pone.0276866.s016

(XLSX)

S3 Table. Results of the raw (p) and Bonferroni-corrected (p.adj) pairwise t-tests between individuals of different types (males, parous and nulliparous females).

Values reported as 0 are actually extremely small, but below the reporting threshold of the pairwise.t.test function in R.

https://doi.org/10.1371/journal.pone.0276866.s017

(XLSX)

S1 File. Supplementary methods and scripts.

https://doi.org/10.1371/journal.pone.0276866.s018

(DOCX)

Acknowledgments

We thank Hosanna Schroeder of Bruker for guidance in data collection.

References

1. Klevezal G. A., Recording structures of mammals. Boca Raton, Florida: CRC Press, 1995.
2. Mahoney P. et al., “Growth of Neanderthal infants from Krapina (120–130 ka), Croatia,” Proc. R. Soc. B, vol. 288, no. 1963, p. 20212079, 2021. pmid:34814754
- View Article
- PubMed/NCBI
- Google Scholar
3. Müller W., Nava A., Evans D., Rossi P. F., Alt K. W., and Bondioli L., “Enamel mineralization and compositional time-resolution in human teeth evaluated via histologically-defined LA-ICPMS profiles,” Geochim. Cosmochim. Acta, vol. 255, pp. 105–126, 2019.
- View Article
- Google Scholar
4. Dean M. C. and Elamin F., “Parturition lines in modern human wisdom tooth roots: do they exist, can they be characterized and are they useful for retrospective determination of age at first reproduction and/or inter-birth intervals?,” Ann. Hum. Biol., vol. 41, no. 4, pp. 358–367, Jul. 2014, pmid:24932749
- View Article
- PubMed/NCBI
- Google Scholar
5. Clark C. T., Horstmann L., and Misarti N., “Evaluating tooth strontium and barium as indicators of weaning age in Pacific walruses,” Methods Ecol. Evol., 2020. pmid:33381293
- View Article
- PubMed/NCBI
- Google Scholar
6. Cerrito P. et al., “Dental cementum virtual histology of Neanderthal teeth from Krapina (Croatia, 130–120 kyr): an informed estimate of age, sex and adult stressors,” J. R. Soc. Interface, vol. 19, no. 187, p. 20210820, 2022. pmid:35193386
- View Article
- PubMed/NCBI
- Google Scholar
7. Kagerer P. and Grupe G., “Age-at-death diagnosis and determination of life-history parameters by incremental lines in human dental cementum as an identification aid,” Forensic Sci. Int., vol. 118, no. 1, pp. 75–82, 2001. pmid:11343858
- View Article
- PubMed/NCBI
- Google Scholar
8. Lieberman D. E., “Life history variables preserved in dental cementum microstructure,” Science, vol. 261, no. 5125, pp. 1162–1164, 1993. pmid:8356448
- View Article
- PubMed/NCBI
- Google Scholar
9. Bromage T. G., Juwayeyi Y. M., Smolyar I., Hu B., Gomez S., and Chisi J., “Enamel-calibrated lamellar bone reveals long period growth rate variability in humans,” Cells Tissues Organs, vol. 194, no. 2–4, pp. 124–130, 2011. pmid:21525718
- View Article
- PubMed/NCBI
- Google Scholar
10. Cipriano A., “Cold stress in captive great apes recorded in incremental lines of dental cementum,” Folia Primatol. (Basel), vol. 73, no. 1, pp. 21–31, 2002. pmid:12065938
- View Article
- PubMed/NCBI
- Google Scholar
11. Marín-Moratalla N., Jordana X., and Köhler M., “Bone histology as an approach to providing data on certain key life history traits in mammals: implications for conservation biology,” Mamm. Biol., vol. 78, no. 6, pp. 422–429, 2013.
- View Article
- Google Scholar
12. Smith T. M. et al., “Wintertime stress, nursing, and lead exposure in Neanderthal children,” Sci. Adv., vol. 4, no. 10, p. eaau9483, 2018. pmid:30402544
- View Article
- PubMed/NCBI
- Google Scholar
13. Boyde A. and Bromage T. G., “Histology of Dental Hard Tissues,” in Vertebrate Skeletal Histology and Paleohistology, CRC Press, 2021, pp. 259–288.
14. Burr D. B., “Estimated intracortical bone turnover in the femur of growing macaques: implications for their use as models in skeletal pathology,” Anat. Rec., vol. 232, no. 2, pp. 180–189, 1992. pmid:1546797
- View Article
- PubMed/NCBI
- Google Scholar
15. McFarlin S. C., Terranova C. J., Zihlman A. L., and Bromage T. G., “Primary bone microanatomy records developmental aspects of life history in catarrhine primates,” J. Hum. Evol., vol. 92, pp. 60–79, 2016. pmid:26989017
- View Article
- PubMed/NCBI
- Google Scholar
16. Bromage T. G. et al., “Lamellar bone is an incremental tissue reconciling enamel rhythms, body size, and organismal life history,” Calcif. Tissue Int., vol. 84, no. 5, pp. 388–404, 2009. pmid:19234658
- View Article
- PubMed/NCBI
- Google Scholar
17. Kovacs C. S., “Maternal mineral and bone metabolism during pregnancy, lactation, and post-weaning recovery,” Physiol. Rev., 2016.
- View Article
- Google Scholar
18. Widdowson E., “Metabolic relationship of calcium, magnesium and phosphorus in the foetus and newly born.,” Voeding, vol. 23, pp. 62–71, 1962. pmid:14006729
- View Article
- PubMed/NCBI
- Google Scholar
19. Comar C. L., “Radiocalcium studies in pregnancy,” Ann. N. Y. Acad. Sci., vol. 64, no. 3, pp. 281–298, 1956.
- View Article
- Google Scholar
20. de Flamingh J. and van der Merve J., “A serum biochemical profile of normal pregnancy,” S. Afr. Med. J., vol. 65, no. 14, pp. 552–555, 1984. pmid:6369570
- View Article
- PubMed/NCBI
- Google Scholar
21. Li Q., Nava A., Reynard L. M., Thirlwall M., Bondioli L., and Müller W., “Spatially-resolved Ca isotopic and trace element variations in human deciduous teeth record diet and physiological change,” Environ. Archaeol., pp. 1–10, 2020.
- View Article
- Google Scholar
22. Austin C. et al., “Barium distributions in teeth reveal early-life dietary transitions in primates,” Nature, vol. 498, no. 7453, pp. 216–219, Jun. 2013, pmid:23698370
- View Article
- PubMed/NCBI
- Google Scholar
23. Clark C. T., Horstmann L., and Misarti N., “Zinc concentrations in teeth of female walruses reflect the onset of reproductive maturity,” Conserv. Physiol., vol. 8, no. 1, 2020. pmid:32308984
- View Article
- PubMed/NCBI
- Google Scholar
24. Rahnama M., “Influence of estrogen deficiency on the copper level in rat teeth and mandible.,” in Annales Universitatis Mariae Curie-Sklodowska. Sectio D: Medicina, 2002, vol. 57, no. 1, pp. 352–356. pmid:12898943
- View Article
- PubMed/NCBI
- Google Scholar
25. Arnold W. H. and Gaengler P., “Quantitative analysis of the calcium and phosphorus content of developing and permanent human teeth,” Ann. Anat.-Anat. Anz., vol. 189, no. 2, pp. 183–190, 2007. pmid:17419551
- View Article
- PubMed/NCBI
- Google Scholar
26. McCormack M. A., McFee W. E., Whitehead H. R., Piwetz S., and Dutton J., “Exploring the Use of SEM–EDS Analysis to Measure the Distribution of Major, Minor, and Trace Elements in Bottlenose Dolphin (Tursiops truncatus) Teeth,” Biol. Trace Elem. Res., pp. 1–13, 2021.
- View Article
- Google Scholar
27. Cerrito P., Cerrito L., Hu B., Bailey S. E., Kalisher R., and Bromage T. G., “Weaning, parturitions and illnesses are recorded in rhesus macaque (Macaca mulatta) dental cementum microstructure,” Am. J. Primatol., vol. 83, no. 3, p. e23235, 2021. pmid:33522634
- View Article
- PubMed/NCBI
- Google Scholar
28. Goldman H. M., Kindsvater J., and Bromage T. G., “Correlative light and backscattered electron microscopy of bone—Part I: Specimen preparation methods,” Scanning, vol. 21, no. 1, pp. 40–43, 1999. pmid:10070782
- View Article
- PubMed/NCBI
- Google Scholar
29. Boyde A., Elliott J. C., and Jones S. J., “Stereology and histogram analysis of backscattered electron images: age changes in bone,” Bone, vol. 14, no. 3, pp. 205–210, 1993. pmid:8363858
- View Article
- PubMed/NCBI
- Google Scholar
30. Howell P. G. T. and Boyde A., “Monte Carlo simulations of electron scattering in bone,” Bone, vol. 15, no. 3, pp. 285–291, 1994. pmid:8068449
- View Article
- PubMed/NCBI
- Google Scholar
31. Goldstein J. I. et al., “Quantitative X-ray analysis: the basics,” in Scanning Electron Microscopy and X-ray Microanalysis, Springer, 2003, pp. 391–451.
32. Schwartz S. M. and Kemnitz J. W., “Age-and gender-related changes in body size, adiposity, and endocrine and metabolic parameters in free-ranging rhesus macaques,” Am. J. Phys. Anthropol., vol. 89, no. 1, pp. 109–121, 1992. pmid:1530057
- View Article
- PubMed/NCBI
- Google Scholar
33. Turcotte C. M. et al., “The ontogeny of sexual dimorphism in free‐ranging rhesus macaques,” Am. J. Phys. Anthropol., vol. 177, no. 2, pp. 314–327, 2022. pmid:35571460
- View Article
- PubMed/NCBI
- Google Scholar
34. Gomendio M., “Differences in fertility and suckling patterns between primiparous and multiparous rhesus mothers (Macaca mulatta),” Reproduction, vol. 87, no. 2, pp. 529–542, 1989.
- View Article
- Google Scholar
35. Kessler M. J. and Rawlins R. G., Cayo Santiago Macaques, The: History, Behavior, and Biology. SUNY Press, 1986.
36. Ross C., “Life history patterns and ecology of macaque species,” Primates, vol. 33, no. 2, pp. 207–215, 1992.
- View Article
- Google Scholar
37. Silk J., Short J., Roberts J., and Kusnitz J., “Gestation length in rhesus macaques (Macaca mulatta),” Int. J. Primatol., vol. 14, no. 1, pp. 95–104, 1993.
- View Article
- Google Scholar
38. Reitsema L. J., Partrick K. A., and Muir A. B., “Inter‐individual variation in weaning among rhesus macaques (Macaca mulatta): Serum stable isotope indicators of suckling duration and lactation,” Am. J. Primatol., vol. 78, no. 10, pp. 1113–1134, 2016. pmid:26284697
- View Article
- PubMed/NCBI
- Google Scholar
39. Tanaka I., “Three phases of lactation in free-ranging Japanese macaques,” Anim. Behav., vol. 44, no. 1, pp. 129–139, 1992.
- View Article
- Google Scholar
40. Boucher R. D., Alavi S. E., de Jong H. N., Godfrey L. V., and Vogel E. R., “Stable isotope evidence (Fe, Cu) suggests that sex, but not aging is recorded in rhesus macaque (Macaca mulatta) bone,” Am. J. Phys. Anthropol., vol. 176, no. 1, pp. 80–92, 2021. pmid:33973647
- View Article
- PubMed/NCBI
- Google Scholar
41. Cerroni A. M., Tomlinson G. A., Turnquist J. E., and Grynpas M. D., “Effect of parity on bone mineral density in female rhesus macaques from Cayo Santiago,” Am. J. Phys. Anthropol. Off. Publ. Am. Assoc. Phys. Anthropol., vol. 121, no. 3, pp. 252–269, 2003.
- View Article
- Google Scholar
42. Lujano-Negrete A. Y. et al., “Bone metabolism and osteoporosis during pregnancy and lactation,” Arch. Osteoporos., vol. 17, no. 1, pp. 1–6, 2022. pmid:35182221
- View Article
- PubMed/NCBI
- Google Scholar
43. Pittet F., Johnson C., and Hinde K., “Age at reproductive debut: Developmental predictors and consequences for lactation, infant mass, and subsequent reproduction in rhesus macaques (Macaca mulatta),” Am. J. Phys. Anthropol., vol. 164, no. 3, pp. 457–476, 2017.
- View Article
- Google Scholar
44. Gonen E., Sahin I., Ozbek M., Kovalak E., Yologlu S., and Ates Y., “Effects of pregnancy and lactation on bone mineral density, and their relation to the serum calcium, phosphorus, calcitonin and parathyroid hormone levels in rats,” J. Endocrinol. Invest., vol. 28, no. 6, pp. 322–326, 2005.
- View Article
- Google Scholar
45. Ott S. M., Lipkin E. W., and Newell‐Morris L., “Bone physiology during pregnancy and lactation in young macaques,” J. Bone Miner. Res., vol. 14, no. 10, pp. 1779–1788, 1999. pmid:10491226
- View Article
- PubMed/NCBI
- Google Scholar
46. Edward A., Obeagu E. I., Okorie H. M., Vincent C. C. N., and Bot Y. S., “Studies of Serum Calcium, Inorganic Phosphate and Magnesium Levels in Lactating Mothers in Owerri,” Studies, vol. 33, no. 41B, 2021.
- View Article
- Google Scholar
47. Kumar A. N. and Kumar R., “Correlation of serum magnesium level with increase in parity and comparison among rural and urban pregnant women,” Int J Cur Bio Med Sci, vol. 2, no. 1, pp. 185–187, 2012.
- View Article
- Google Scholar
48. Pathak P., Kapoor S. K., Kapil U., and Dwivedi S. N., “Serum magnesium level among pregnant women in a rural community of Haryana State, India,” Eur. J. Clin. Nutr., vol. 57, no. 11, pp. 1504–1506, 2003. pmid:14576766
- View Article
- PubMed/NCBI
- Google Scholar
49. Kriesten K., Palawinskas R., and Sommer H., “Magnesium levels in the body fluids of dam and fetal rabbits during the reproductive phase.,” Comp. Biochem. Physiol. A, vol. 79, no. 1, pp. 197–200, 1984. pmid:6148182
- View Article
- PubMed/NCBI
- Google Scholar
50. Bromage T. G. et al., “Signposts ahead: Hard tissue signals on rue Armand de Ricqlès,” Comptes Rendus Palevol, vol. 10, no. 5–6, pp. 499–507, 2011.
- View Article
- Google Scholar
51. Mitić Ž. et al., “Instrumental methods and techniques for structural and physicochemical characterization of biomaterials and bone tissue: A review,” Mater. Sci. Eng. C, vol. 79, pp. 930–949, 2017. pmid:28629099
- View Article
- PubMed/NCBI
- Google Scholar
52. Weinbruch S., Wentzel M., Kluckner M., Hoffmann P., and Ortner H. M., “Characterization of individual atmospheric particles by element mapping in electron probe microanalysis,” Microchim. Acta, vol. 125, no. 1, pp. 137–141, 1997.
- View Article
- Google Scholar

[ref1] 1. Klevezal G. A., Recording structures of mammals. Boca Raton, Florida: CRC Press, 1995.

[ref2] 2. Mahoney P. et al., “Growth of Neanderthal infants from Krapina (120–130 ka), Croatia,” Proc. R. Soc. B, vol. 288, no. 1963, p. 20212079, 2021. pmid:34814754
View Article
PubMed/NCBI
Google Scholar

[3] View Article

[4] PubMed/NCBI

[5] Google Scholar

[ref3] 3. Müller W., Nava A., Evans D., Rossi P. F., Alt K. W., and Bondioli L., “Enamel mineralization and compositional time-resolution in human teeth evaluated via histologically-defined LA-ICPMS profiles,” Geochim. Cosmochim. Acta, vol. 255, pp. 105–126, 2019.
View Article
Google Scholar

[7] View Article

[8] Google Scholar

[ref4] 4. Dean M. C. and Elamin F., “Parturition lines in modern human wisdom tooth roots: do they exist, can they be characterized and are they useful for retrospective determination of age at first reproduction and/or inter-birth intervals?,” Ann. Hum. Biol., vol. 41, no. 4, pp. 358–367, Jul. 2014, pmid:24932749
View Article
PubMed/NCBI
Google Scholar

[10] View Article

[11] PubMed/NCBI

[12] Google Scholar

[ref5] 5. Clark C. T., Horstmann L., and Misarti N., “Evaluating tooth strontium and barium as indicators of weaning age in Pacific walruses,” Methods Ecol. Evol., 2020. pmid:33381293
View Article
PubMed/NCBI
Google Scholar

[14] View Article

[15] PubMed/NCBI

[16] Google Scholar

[ref6] 6. Cerrito P. et al., “Dental cementum virtual histology of Neanderthal teeth from Krapina (Croatia, 130–120 kyr): an informed estimate of age, sex and adult stressors,” J. R. Soc. Interface, vol. 19, no. 187, p. 20210820, 2022. pmid:35193386
View Article
PubMed/NCBI
Google Scholar

[18] View Article

[19] PubMed/NCBI

[20] Google Scholar

[ref7] 7. Kagerer P. and Grupe G., “Age-at-death diagnosis and determination of life-history parameters by incremental lines in human dental cementum as an identification aid,” Forensic Sci. Int., vol. 118, no. 1, pp. 75–82, 2001. pmid:11343858
View Article
PubMed/NCBI
Google Scholar

[22] View Article

[23] PubMed/NCBI

[24] Google Scholar

[ref8] 8. Lieberman D. E., “Life history variables preserved in dental cementum microstructure,” Science, vol. 261, no. 5125, pp. 1162–1164, 1993. pmid:8356448
View Article
PubMed/NCBI
Google Scholar

[26] View Article

[27] PubMed/NCBI

[28] Google Scholar

[ref9] 9. Bromage T. G., Juwayeyi Y. M., Smolyar I., Hu B., Gomez S., and Chisi J., “Enamel-calibrated lamellar bone reveals long period growth rate variability in humans,” Cells Tissues Organs, vol. 194, no. 2–4, pp. 124–130, 2011. pmid:21525718
View Article
PubMed/NCBI
Google Scholar

[30] View Article

[31] PubMed/NCBI

[32] Google Scholar

[ref10] 10. Cipriano A., “Cold stress in captive great apes recorded in incremental lines of dental cementum,” Folia Primatol. (Basel), vol. 73, no. 1, pp. 21–31, 2002. pmid:12065938
View Article
PubMed/NCBI
Google Scholar

[34] View Article

[35] PubMed/NCBI

[36] Google Scholar

[ref11] 11. Marín-Moratalla N., Jordana X., and Köhler M., “Bone histology as an approach to providing data on certain key life history traits in mammals: implications for conservation biology,” Mamm. Biol., vol. 78, no. 6, pp. 422–429, 2013.
View Article
Google Scholar

[38] View Article

[39] Google Scholar

[ref12] 12. Smith T. M. et al., “Wintertime stress, nursing, and lead exposure in Neanderthal children,” Sci. Adv., vol. 4, no. 10, p. eaau9483, 2018. pmid:30402544
View Article
PubMed/NCBI
Google Scholar

[41] View Article

[42] PubMed/NCBI

[43] Google Scholar

[ref13] 13. Boyde A. and Bromage T. G., “Histology of Dental Hard Tissues,” in Vertebrate Skeletal Histology and Paleohistology, CRC Press, 2021, pp. 259–288.

[ref14] 14. Burr D. B., “Estimated intracortical bone turnover in the femur of growing macaques: implications for their use as models in skeletal pathology,” Anat. Rec., vol. 232, no. 2, pp. 180–189, 1992. pmid:1546797
View Article
PubMed/NCBI
Google Scholar

[46] View Article

[47] PubMed/NCBI

[48] Google Scholar

[ref15] 15. McFarlin S. C., Terranova C. J., Zihlman A. L., and Bromage T. G., “Primary bone microanatomy records developmental aspects of life history in catarrhine primates,” J. Hum. Evol., vol. 92, pp. 60–79, 2016. pmid:26989017
View Article
PubMed/NCBI
Google Scholar

[50] View Article

[51] PubMed/NCBI

[52] Google Scholar

[ref16] 16. Bromage T. G. et al., “Lamellar bone is an incremental tissue reconciling enamel rhythms, body size, and organismal life history,” Calcif. Tissue Int., vol. 84, no. 5, pp. 388–404, 2009. pmid:19234658
View Article
PubMed/NCBI
Google Scholar

[54] View Article

[55] PubMed/NCBI

[56] Google Scholar

[ref17] 17. Kovacs C. S., “Maternal mineral and bone metabolism during pregnancy, lactation, and post-weaning recovery,” Physiol. Rev., 2016.
View Article
Google Scholar

[58] View Article

[59] Google Scholar

[ref18] 18. Widdowson E., “Metabolic relationship of calcium, magnesium and phosphorus in the foetus and newly born.,” Voeding, vol. 23, pp. 62–71, 1962. pmid:14006729
View Article
PubMed/NCBI
Google Scholar

[61] View Article

[62] PubMed/NCBI

[63] Google Scholar

[ref19] 19. Comar C. L., “Radiocalcium studies in pregnancy,” Ann. N. Y. Acad. Sci., vol. 64, no. 3, pp. 281–298, 1956.
View Article
Google Scholar

[65] View Article

[66] Google Scholar

[ref20] 20. de Flamingh J. and van der Merve J., “A serum biochemical profile of normal pregnancy,” S. Afr. Med. J., vol. 65, no. 14, pp. 552–555, 1984. pmid:6369570
View Article
PubMed/NCBI
Google Scholar

[68] View Article

[69] PubMed/NCBI

[70] Google Scholar

[ref21] 21. Li Q., Nava A., Reynard L. M., Thirlwall M., Bondioli L., and Müller W., “Spatially-resolved Ca isotopic and trace element variations in human deciduous teeth record diet and physiological change,” Environ. Archaeol., pp. 1–10, 2020.
View Article
Google Scholar

[72] View Article

[73] Google Scholar

[ref22] 22. Austin C. et al., “Barium distributions in teeth reveal early-life dietary transitions in primates,” Nature, vol. 498, no. 7453, pp. 216–219, Jun. 2013, pmid:23698370
View Article
PubMed/NCBI
Google Scholar

[75] View Article

[76] PubMed/NCBI

[77] Google Scholar

[ref23] 23. Clark C. T., Horstmann L., and Misarti N., “Zinc concentrations in teeth of female walruses reflect the onset of reproductive maturity,” Conserv. Physiol., vol. 8, no. 1, 2020. pmid:32308984
View Article
PubMed/NCBI
Google Scholar

[79] View Article

[80] PubMed/NCBI

[81] Google Scholar

[ref24] 24. Rahnama M., “Influence of estrogen deficiency on the copper level in rat teeth and mandible.,” in Annales Universitatis Mariae Curie-Sklodowska. Sectio D: Medicina, 2002, vol. 57, no. 1, pp. 352–356. pmid:12898943
View Article
PubMed/NCBI
Google Scholar

[83] View Article

[84] PubMed/NCBI

[85] Google Scholar

[ref25] 25. Arnold W. H. and Gaengler P., “Quantitative analysis of the calcium and phosphorus content of developing and permanent human teeth,” Ann. Anat.-Anat. Anz., vol. 189, no. 2, pp. 183–190, 2007. pmid:17419551
View Article
PubMed/NCBI
Google Scholar

[87] View Article

[88] PubMed/NCBI

[89] Google Scholar

[ref26] 26. McCormack M. A., McFee W. E., Whitehead H. R., Piwetz S., and Dutton J., “Exploring the Use of SEM–EDS Analysis to Measure the Distribution of Major, Minor, and Trace Elements in Bottlenose Dolphin (Tursiops truncatus) Teeth,” Biol. Trace Elem. Res., pp. 1–13, 2021.
View Article
Google Scholar

[91] View Article

[92] Google Scholar

[ref27] 27. Cerrito P., Cerrito L., Hu B., Bailey S. E., Kalisher R., and Bromage T. G., “Weaning, parturitions and illnesses are recorded in rhesus macaque (Macaca mulatta) dental cementum microstructure,” Am. J. Primatol., vol. 83, no. 3, p. e23235, 2021. pmid:33522634
View Article
PubMed/NCBI
Google Scholar

[94] View Article

[95] PubMed/NCBI

[96] Google Scholar

[ref28] 28. Goldman H. M., Kindsvater J., and Bromage T. G., “Correlative light and backscattered electron microscopy of bone—Part I: Specimen preparation methods,” Scanning, vol. 21, no. 1, pp. 40–43, 1999. pmid:10070782
View Article
PubMed/NCBI
Google Scholar

[98] View Article

[99] PubMed/NCBI

[100] Google Scholar

[ref29] 29. Boyde A., Elliott J. C., and Jones S. J., “Stereology and histogram analysis of backscattered electron images: age changes in bone,” Bone, vol. 14, no. 3, pp. 205–210, 1993. pmid:8363858
View Article
PubMed/NCBI
Google Scholar

[102] View Article

[103] PubMed/NCBI

[104] Google Scholar

[ref30] 30. Howell P. G. T. and Boyde A., “Monte Carlo simulations of electron scattering in bone,” Bone, vol. 15, no. 3, pp. 285–291, 1994. pmid:8068449
View Article
PubMed/NCBI
Google Scholar

[106] View Article

[107] PubMed/NCBI

[108] Google Scholar

[ref31] 31. Goldstein J. I. et al., “Quantitative X-ray analysis: the basics,” in Scanning Electron Microscopy and X-ray Microanalysis, Springer, 2003, pp. 391–451.

[ref32] 32. Schwartz S. M. and Kemnitz J. W., “Age-and gender-related changes in body size, adiposity, and endocrine and metabolic parameters in free-ranging rhesus macaques,” Am. J. Phys. Anthropol., vol. 89, no. 1, pp. 109–121, 1992. pmid:1530057
View Article
PubMed/NCBI
Google Scholar

[111] View Article

[112] PubMed/NCBI

[113] Google Scholar

[ref33] 33. Turcotte C. M. et al., “The ontogeny of sexual dimorphism in free‐ranging rhesus macaques,” Am. J. Phys. Anthropol., vol. 177, no. 2, pp. 314–327, 2022. pmid:35571460
View Article
PubMed/NCBI
Google Scholar

[115] View Article

[116] PubMed/NCBI

[117] Google Scholar

[ref34] 34. Gomendio M., “Differences in fertility and suckling patterns between primiparous and multiparous rhesus mothers (Macaca mulatta),” Reproduction, vol. 87, no. 2, pp. 529–542, 1989.
View Article
Google Scholar

[119] View Article

[120] Google Scholar

[ref35] 35. Kessler M. J. and Rawlins R. G., Cayo Santiago Macaques, The: History, Behavior, and Biology. SUNY Press, 1986.

[ref36] 36. Ross C., “Life history patterns and ecology of macaque species,” Primates, vol. 33, no. 2, pp. 207–215, 1992.
View Article
Google Scholar

[123] View Article

[124] Google Scholar

[ref37] 37. Silk J., Short J., Roberts J., and Kusnitz J., “Gestation length in rhesus macaques (Macaca mulatta),” Int. J. Primatol., vol. 14, no. 1, pp. 95–104, 1993.
View Article
Google Scholar

[126] View Article

[127] Google Scholar

[ref38] 38. Reitsema L. J., Partrick K. A., and Muir A. B., “Inter‐individual variation in weaning among rhesus macaques (Macaca mulatta): Serum stable isotope indicators of suckling duration and lactation,” Am. J. Primatol., vol. 78, no. 10, pp. 1113–1134, 2016. pmid:26284697
View Article
PubMed/NCBI
Google Scholar

[129] View Article

[130] PubMed/NCBI

[131] Google Scholar

[ref39] 39. Tanaka I., “Three phases of lactation in free-ranging Japanese macaques,” Anim. Behav., vol. 44, no. 1, pp. 129–139, 1992.
View Article
Google Scholar

[133] View Article

[134] Google Scholar

[ref40] 40. Boucher R. D., Alavi S. E., de Jong H. N., Godfrey L. V., and Vogel E. R., “Stable isotope evidence (Fe, Cu) suggests that sex, but not aging is recorded in rhesus macaque (Macaca mulatta) bone,” Am. J. Phys. Anthropol., vol. 176, no. 1, pp. 80–92, 2021. pmid:33973647
View Article
PubMed/NCBI
Google Scholar

[136] View Article

[137] PubMed/NCBI

[138] Google Scholar

[ref41] 41. Cerroni A. M., Tomlinson G. A., Turnquist J. E., and Grynpas M. D., “Effect of parity on bone mineral density in female rhesus macaques from Cayo Santiago,” Am. J. Phys. Anthropol. Off. Publ. Am. Assoc. Phys. Anthropol., vol. 121, no. 3, pp. 252–269, 2003.
View Article
Google Scholar

[140] View Article

[141] Google Scholar

[ref42] 42. Lujano-Negrete A. Y. et al., “Bone metabolism and osteoporosis during pregnancy and lactation,” Arch. Osteoporos., vol. 17, no. 1, pp. 1–6, 2022. pmid:35182221
View Article
PubMed/NCBI
Google Scholar

[143] View Article

[144] PubMed/NCBI

[145] Google Scholar

[ref43] 43. Pittet F., Johnson C., and Hinde K., “Age at reproductive debut: Developmental predictors and consequences for lactation, infant mass, and subsequent reproduction in rhesus macaques (Macaca mulatta),” Am. J. Phys. Anthropol., vol. 164, no. 3, pp. 457–476, 2017.
View Article
Google Scholar

[147] View Article

[148] Google Scholar

[ref44] 44. Gonen E., Sahin I., Ozbek M., Kovalak E., Yologlu S., and Ates Y., “Effects of pregnancy and lactation on bone mineral density, and their relation to the serum calcium, phosphorus, calcitonin and parathyroid hormone levels in rats,” J. Endocrinol. Invest., vol. 28, no. 6, pp. 322–326, 2005.
View Article
Google Scholar

[150] View Article

[151] Google Scholar

[ref45] 45. Ott S. M., Lipkin E. W., and Newell‐Morris L., “Bone physiology during pregnancy and lactation in young macaques,” J. Bone Miner. Res., vol. 14, no. 10, pp. 1779–1788, 1999. pmid:10491226
View Article
PubMed/NCBI
Google Scholar

[153] View Article

[154] PubMed/NCBI

[155] Google Scholar

[ref46] 46. Edward A., Obeagu E. I., Okorie H. M., Vincent C. C. N., and Bot Y. S., “Studies of Serum Calcium, Inorganic Phosphate and Magnesium Levels in Lactating Mothers in Owerri,” Studies, vol. 33, no. 41B, 2021.
View Article
Google Scholar

[157] View Article

[158] Google Scholar

[ref47] 47. Kumar A. N. and Kumar R., “Correlation of serum magnesium level with increase in parity and comparison among rural and urban pregnant women,” Int J Cur Bio Med Sci, vol. 2, no. 1, pp. 185–187, 2012.
View Article
Google Scholar

[160] View Article

[161] Google Scholar

[ref48] 48. Pathak P., Kapoor S. K., Kapil U., and Dwivedi S. N., “Serum magnesium level among pregnant women in a rural community of Haryana State, India,” Eur. J. Clin. Nutr., vol. 57, no. 11, pp. 1504–1506, 2003. pmid:14576766
View Article
PubMed/NCBI
Google Scholar

[163] View Article

[164] PubMed/NCBI

[165] Google Scholar

[ref49] 49. Kriesten K., Palawinskas R., and Sommer H., “Magnesium levels in the body fluids of dam and fetal rabbits during the reproductive phase.,” Comp. Biochem. Physiol. A, vol. 79, no. 1, pp. 197–200, 1984. pmid:6148182
View Article
PubMed/NCBI
Google Scholar

[167] View Article

[168] PubMed/NCBI

[169] Google Scholar

[ref50] 50. Bromage T. G. et al., “Signposts ahead: Hard tissue signals on rue Armand de Ricqlès,” Comptes Rendus Palevol, vol. 10, no. 5–6, pp. 499–507, 2011.
View Article
Google Scholar

[171] View Article

[172] Google Scholar

[ref51] 51. Mitić Ž. et al., “Instrumental methods and techniques for structural and physicochemical characterization of biomaterials and bone tissue: A review,” Mater. Sci. Eng. C, vol. 79, pp. 930–949, 2017. pmid:28629099
View Article
PubMed/NCBI
Google Scholar

[174] View Article

[175] PubMed/NCBI

[176] Google Scholar

[ref52] 52. Weinbruch S., Wentzel M., Kluckner M., Hoffmann P., and Ortner H. M., “Characterization of individual atmospheric particles by element mapping in electron probe microanalysis,” Microchim. Acta, vol. 125, no. 1, pp. 137–141, 1997.
View Article
Google Scholar

[178] View Article

[179] Google Scholar

Figures

Abstract

Introduction

Materials and methods

Specimen information

Block preparation and imaging

Elemental data collection

Thin section preparation and imaging

Chronological age assignation

Data preparation

Data analysis

Results

Prediction 1

Prediction 2

Prediction 3

Prediction 4

Discussion

Supporting information

S1 Fig. Montaged micrograph of the femur of MA024 obtained in reflected light.

S2 Fig. Montaged micrograph of the femur of M591 obtained in reflected light.

S3 Fig. Montaged micrograph of the femur of M526 obtained in reflected light.

S4 Fig. Montaged micrograph of the femur of M361 obtained in reflected light.

S5 Fig. Montaged micrograph of the femur of CK85 obtained in reflected light.

S6 Fig. Montaged micrograph of the femur of BD24 obtained in reflected light.

S7 Fig. The red lines indicate the locations of transects 1 and 2.

S8 Fig. PCA of relative percentage data grouped by individual type (both segments from the same individual are pooled together) and colored according to the age associated with each datapoint of each individual.

S9 Fig. PCA of relative percentage data of only females, colored by parity state (blue = nulliparous, yellow = parous).

S10 Fig. Boxplot showing the difference in relative percent concentrations of calcium in the bone of M526 formed during gestation and during lactation.

S11 Fig. Boxplot showing the difference in relative percent concentrations of phosphorous in the bone of M526 formed during gestation and during lactation.

S12 Fig. Boxplot showing the difference in relative percent concentrations of magnesium in the bone of M526 formed during reproduction (gestation and lactation) or during times of no known physiologically challenging events (NKE).

S13 Fig. Boxplot showing the difference in relative percent concentrations of magnesium in the bone of CK85 formed during breastfeeding (suckling and weaning) or during times of no known physiologically challenging events (NKE).

S14 Fig. Boxplot showing the difference in relative percent concentrations of oxygen in the bone of M526 formed during reproduction (gestation and lactation) or during times of no known physiologically challenging events (NKE).

S1 Table. List of age at known impactful physiological events for each individual.

S2 Table. Detailed measurements for each transect.

S3 Table. Results of the raw (p) and Bonferroni-corrected (p.adj) pairwise t-tests between individuals of different types (males, parous and nulliparous females).

S1 File. Supplementary methods and scripts.

Acknowledgments

References