Biological rhythms exert considerable regulatory control over an organism’s physiology, behavior, and development, which are particularly visible as daily metabolic fluctuations attributable to the central circadian clock [1]. Dental hard tissues (e.g., enamel and dentine) reflect these cycles, but also provide signatures of much longer periods (e.g., days to years), offering a unique window into the study of organismal life history. Bone, on the other hand, has not been recognized as a rhythmic tissue except for intermittent growth cycles commonly observed in seasonally growing animals [2].

Bone forming cells, the osteoblasts, secrete their organic matrix in discrete tissue patterns, which frequently include formation of the lamella, a fundamental microanatomical unit of bone observed in the postnatal tissues of many vertebrate taxa [3, 4]. The lamella is characterized in histological thin sections by a highly oriented band of collagen separated from adjacent bands by an interlamellar zone of less oriented collagen [5], resulting in the discrimination of one lamella from the next in polarized light. Lamellar bone is strikingly incremental in appearance, but the periodic formation of lamellae is the subject of only two preliminary reports [6, 7]. In this early research a 24-h rhythm was recorded for the formation of one lamella in three small mammals and a longer but uncertain rhythm in one larger mammal.

Mammalian teeth, on the other hand, are well documented to exhibit incremental structures representing successive forming fronts of enamel at varying time scales. Enamel forming cells, the ameloblasts, secrete their matrix with a circadian (ca. 24-h) periodicity, giving rise to daily increments referred to as cross striations, or alternating varicosities and constrictions, which are observed by a variety of microscopy techniques and widely used in studies of primate dental development and life history [810]. Enamel also exhibits a long-period increment called the stria of Retzius [11], which represents a temporary but profound periodic slowing of enamel formation in various multiples of whole days. The number of cross striations between striae is referred to as the repeat interval, which is identical for all teeth of an individual organism but is variable between and sometimes within species [12].

Cross striations are the result of daily variation in ameloblast secretory rate [10]. The etiology of striae of Retzius is unknown despite the exceedingly periodic nature of their expression. Equally enigmatic is the suggestion that the repeat interval scales with body size. Among primates it has been observed that smaller-bodied taxa have short repeat intervals while those of larger bodied primates are longer, suggesting a correlation between repeat interval and body size [13, 14], although no functional explanation for the scaling of striae of Retzius with body size has been asserted.

It is not unreasonable to assume that growth increments in bone would show a relationship with body size, insofar as the scaling of bone mass to body mass is an axiom of vertebrate hard tissue biology [15]. Because enamel and bone are each incrementally formed by cells secreting an organic matrix that is subsequently mineralized, we hypothesize that enamel tissue rhythms will have a correspondence with the incremental growth of lamellar bone, and that this association should be reflected during bone growth and bone remodeling processes that relate to an organism’s body mass.

To test this hypothesis, we investigate ‘ground’ sections of the bones and teeth of a large sample of mammalian taxa, using a variety of analytical microscopy techniques. In three steps we aim (i) to further validate a relationship between striae repeat interval and body mass across many mammalian taxa, including fossil hominins (extinct human ancestors); (ii) to determine whether mammalian osteocyte density—a surrogate for inferring rates of bone growth and remodeling—is related to body mass across various Primate and non-Primate taxa; and (iii) having independently established relationships between both stria repeat interval and bone formation rates with body mass, to logically extend the argument and establish if there is a relation between the enamel striae of Retzius and the formation rate of bone lamellae.

Finally, we also aim to contextualize our results with a literature on the central autonomic regulation of bone/body mass, thus reconciling our investigation of long-period rhythms with skeletal development and organismal life history research. This step, while speculative, is important for directing the focus of future experimental research onto physiological pathways implicated in the control of long-period rhythms.

Materials and Methods

Specimens

To support and expand earlier research on the relation between striae of Retzius repeat interval and body mass, tooth samples from 42 modern and fossil species/subspecies of the order Primates, representing the infraorders Platyrrhini and Catarrhini and the suborder Strepsirhini are investigated (Table 1). This relationship is as well examined in a compilation of striae of Retzius repeat intervals and body masses for 11 extinct and living members of the human lineage (Table 2). To determine whether this relationship might also be extended to another and disparate mammalian taxon, seven members of the order Proboscidea were investigated (Table 3).

Table 1 Primate striae of Retzius repeat interval (RI) and body mass
Table 2 Hominin striae of Retzius repeat interval (RI) and body mass
Table 3 Proboscidean striae of Retzius repeat interval (RI) and body mass

With respect to the relationship between osteocyte density and body mass, measurements were made on eight primates representing the infraorder Catarrhini and the suborder Strepsirhini and three nonprimate mammals (Table 4). Measurements for human were represented by five individual females aged 38–62 from people of known life history.

Table 4 Selected primate and nonprimate mammal osteocyte density and body mass

To investigate striae of Retzius repeat intervals and their relationships with lamellar formation rates, we acquired the vitally labeled immature bone of five mammals for which the times of administration of the labels were known. Table 5 provides labeling schedules for rat [16], macaque [17], patas monkey (this study), sheep [18], and human (inferred) (see published sources for details of the labeling protocols). Teeth were obtained from the same species/variety or, in the case of the patas monkey, from the same individual.

Table 5 Striae of Retzius repeat interval (RI) and lamellar formation rate

Sample Preparation

All modern bone and tooth samples directly investigated by us were cleaned of nonmineralized organic components with daily changes of fresh 1% Terg-a-Zyme (Alconox, NY) aqueous enzyme solution at 50°C until clean [19]. Following this, samples were subject to graded ethanol substitution and then 50:50 isopropanol:heptane reflux in a Soxhlet apparatus for 7 to 14 days.

Upon removal from the Soxhlet apparatus, all modern and fossil samples, excepting early hominin specimens, were transferred and subject to polymethylmethacrylate (PMMA) substitution and embedding. Each cured block was sawn through midshaft with a Buehler (Lake Bluff, IL) Isomet low-speed saw, hand ground through graded carbide papers to 1200 grit on a Buehler Handimet II, mounted on a glass slide, and sectioned, and the section was polished on a Buehler Ecomet III with a 1-μm diamond suspension until a uniform thickness of 100 ± 4 μm was achieved, in most cases according to a published protocol [20].

Histological sections were typically prepared to a section thickness standard of 100 ± 4 μm. All ± 4-μm section thickness monitoring was determined using an Edge R400 3D Transmitted Light Microscope [21, 22] (Edge R400: Micro Science Technologies, Marina del Ray, CA) fitted with calibrated 1/10-μm staging in the Z direction and using a Nikon PlanApo 40/0.95 objective lens.

No early hominin specimen examined for this research was sectioned, these being examined by confocal microscopy (see Microscopy and Analysis below), thus requiring no preparation.

Microscopy and Analysis

Striae of Retzius and lamellae

Histological ‘ground’ thin sections were coverslipped with ethylene glycol and imaged in conventional transmitted light (LM) and circularly polarized transmitted light (CPL) using a Leica-Leitz DMRX/E Universal Microscope configured with a Marzhauser motorized stage and CPL filters and using Leica PL Fluotar 40/0.70 (enamel) and PL Fluotar 20/0.50 (lamellae) objective lenses. LM and CPL images were acquired with a JVC KYF55B color videocamera using Syncroscopy Montage Explorer software (Synoptics Inc., Frederick, MD).

Evaluations of the number of daily increments between adjacent striae of Retzius, or the repeat interval, were performed on images compiled into large XY montages using an automated Montage Explorer capture. In addition, averages of cross-striation widths in micrometers were divided into averages of the distance between striae of Retzius for two regions of the same tooth. This allowed a second estimate of the repeat interval; a value was recorded only when both sets of numbers agreed.

For vitally labeled bone, fluorescence microscopy (FM) was also performed on the same LM and CPL fields of view. When fluorescence was weak or diffuse, FM was performed with a Perkin Elmer (Waltham, MA) UltraView ERS 3E situated on an AxioImager M1 Zeiss microscope (Carl Zeiss MicroImaging, Inc.) and using a Zeiss LD Plan-Neofluor 40/0.60 objective lens. The 488- and 568-nm laser lines were used to excite vital labels. Automated Z-axis image capture on the Axiovert was performed with a piezo-z-stepper. Images were acquired with a cooled 12-bit CCD Orca ER camera at magnifications and fields of view matching Leitz DMRX/E images. Digital image overlays between DMRX/E CPL and UltraView ERS 3E FM images were performed using crop, rotation, and screen functions in Adobe Photoshop CS2 (San Jose, CA).

Early hominin images of enamel and bone were obtained by Portable Confocal Scanning Optical Microscopy (PCSOM) [23], employing a Technical Instrument Co. K2S-BIO confocal module (Zygo Corp., Sunnyvale, CA) with conventional white light, a Mitutoyo M Plan Apo 20/0.42 objective lens, and a JVC KYF1030 camera. Using Montage Explorer software we obtained high-resolution images of early hominin internal structures, such as cross striations and striae of Retzius from naturally occurring occlusocervical fractures in enamel, and lamellae from intact and fractured bones.

To demonstrate that lamellae develop incrementally, when specimens allowed (human and sheep) we imaged both primary and secondarily remodeled lamellar bone by FM.

Osteocyte Density

All imaging, excepting that of the early hominin specimen, was performed with the Edge R400 as described above, configured to provide transmitted oblique circularly polarized light contrast as useful. Sections were coverslipped with ethylene glycol, and XY fields of view conforming to a 270-μm field width by 200-μm field height, respectively, and containing no dehydration cracks or vascular canals, were located. The periosteal cortex was imaged in all individuals, including, for purposes of control, specifically the lateral periosteal cortex from five modern humans, the average of which is given in Table 4.

An XY field of view selected, the focal plane was situated at the top surface of the section and then lowered in the Z-axis by a virtual translation of ca. 5 μm through the section. The through series of images from each ca. 5-μm virtual focal plane was acquired using a JVC KYF75U camera and Syncroscopy Automontage (Synoptics Inc., Frederick, MD) until the final in-focus bottom surface image plane had been reached. A red-green anaglyph image with 20-deg angular separation was constructed from each through series with Automontage.

Each through series was presented for analysis on a PC configured with dual monitors. The Automontage through series was placed on one monitor, while the analglyph image was opened separately in Adobe Photoshop CS2 (San Jose, CA) on the other; the anaglyph was scaled to the same screen size as the Automontage source images. This side-by-side arrangement allowed the viewing of lacunae with red-green anaglyph glasses through the full 3D thickness of the histological section while also scrolling through and cross-checking source image details in Automontage. All lacunae were counted for each XYZ volume according to standard stereological procedure; those intersecting each of one X, Y, and Z surface were excluded, while those intersecting their contralateral X, Y, and Z surfaces were included. After all lacunae had been labeled, they were manually counted. For reasons mentioned in the Results, no attempt was made to subtract from the total count those lacunae whose cells may have been apoptotic. All measurements were standardized by extrapolation of the number of lacunae per measured bone volume, typically from a 5,400,000-μm cubic measurement volume to a 1-mm3 unit value.

Because of satisfactory image contrast obtained with the Edge R400, we are only rarely unsure about the existence of a lacuna. This was verified by intraobserver reinspection and recounting of five fields of view, in which error was limited to 1%. Interobserver error remained below 5%, most of the differences between observers being attributable to the detection of top and bottom surface intersecting lacunae.

Measurement of the early hominin was obtained by the PCSOM as described above [23]. In this case, a 270-μm field width by 200-μm field height field of view was imaged through 50 μm from below the first surface.

Enumerating Lamellae

Bone lamellae in histological thin section appear as microanatomical lineations, but frequently they are domains laid down in patches, superficial to but not strictly overlapping the other, conforming to the three-dimensional shape of the forming front, which itself may not be flat [24]. Lamellae are thus typically anisotropic in two dimensions, appearing and disappearing across the field of view. This anisotropy may pose difficulties for the objective determination of what lineation is or is not a lamella to be counted. For this reason, novel methods for the study of incremental growth rate variability have been developed [25]. In the present study, fields of view containing contiguous lamellae were sought and evaluated; ambiguities associated with the visualization of lamellae are presented in the Discussion.

There are other methodological caveats that must be considered when enumerating lamellae with respect to developmental time. Instances of growth reversal due to cortical drift and transient growth cessation in the bones of seasonally growing animals are two examples of contiguous lamellae not representing continuous growth. In these cases, the micromorphology of a previously resorbed surface and/or a cement line typically demarcates the interruption of growth. While any numbers of observed lamellae have temporal meaning, to use lamellae for representations of continuous time, one must have reason to expect that their development was uninterrupted.

Results

Striae of Retzius Repeat Interval Correlates with Body Mass

The relationship between striae of Retzius and body mass in a range of living and fossil ape taxa has been recently explored [13]. We broaden and focus our examination of this relationship using a large sample of 23 New World primates of the infraorder Platyrrhini (Table 1). Repeat intervals obtained on platyrrhines, together with published data on other primates of the infraorder Catarrhini and suborder Strepsirhini, were regressed against the log of body mass estimates and show a statistically high relationship described by r = 0.79, p < 0.01, r 2 = 0.62, y = 3.136 + 1.04 * x (Fig. 1, Table 1). Employing data from the most comprehensive study and compilation of Plio-Pleistocene hominin repeat intervals to date [26], we also tested this relationship for early hominins, which are found to have repeat intervals ranging from 6 to 9 and which do indeed significantly correlate with estimated body mass and described by the relationship r = 0.87, p < 0.01, r 2 = 0.75, y = 5.262 + 0.0470 * x (Fig. 2, Table 2). To further test the broader applicability of this relationship in mammals, we examined living and extinct and macrofaunal and insular dwarfs of the order Proboscidea, which also reveal a significant correlation between body mass and striae of Retzius repeat interval described by r = 0.79, p = 0.03, r 2 = 0.62, y = 5.649 + 2.098 * x (Fig. 3, Table 3). When all primates and proboscideans examined here are treated together, the relationship is described by r = 0.87, p < 0.01, r 2 = 0.75, y = 3.078 + 1.057 * x.

Fig. 1
figure 1

Striae of Retzius repeat interval (RI) linearly regressed against log body mass for primates: (●) platyrrhine; (□) catarrhine; (◯, Δ) strepsirrhine primates (least squares model). The number of daily cross striations between striae of Retzius is referred to as the RI. Raw data provided in Table 1. Summary statistics for primates combined are as follows: r = 0.79, p < 0.01, r 2 = 0.62, y = 3.136 + 1.04 * x. For platyrrhines, r = 0.92, p < 0.01, r 2 = 0.85, y = 3.451 + 1.416 * x. For catarrhines, r = 0.88, p < 0.01, r 2 = 0.77, y = 1.381 + 1.708 * x. The X-axis represents the log of body mass measured in kilograms

Fig. 2
figure 2

Average striae of Retzius repeat interval (RI) linearly regressed against log body mass for humans and their fossil ancestors (least squares model). The number of daily cross striations between striae of Retzius is referred to as the RI. Raw data are provided in Table 2. Summary statistics are as follows: r = 0.87, p < 0.01, r 2 = 0.75, y = 5.262 + 0.0470 * x. The X-axis represents body mass measured in kilograms

Fig. 3
figure 3

Striae of Retzius repeat interval (RI) linearly regressed against log body mass for proboscideans (least squares model). The number of daily cross striations between striae of Retzius is referred to as the RI. Raw data provided in Table 3. Summary statistics are as follows: r = 0.79, p = 0.03, r 2 = 0.62, y = 5.649 + 2.098 * x. The X-axis represents the log of body mass measured in kilograms

Osteocyte Lacuna Density Correlates with Body Mass

Using the novel three-dimensional tallying method described under Materials and Methods, we compared osteocyte densities from the midshaft femur of eight primate and three nonprimate taxa (Table 4). Primary lamellar or parallel-fibered bone sites could be measured for most taxa. However, only secondary osteonal lamellae were available for measurement on the modern human lateral periosteal cortex selected and on the modern and fossil hippo material. Results indicate that osteocyte density is inversely proportional to body mass for mammals, the relationship being described by r = −0.85, p < 0.01, r 2 = 0.72, y = 10.500 – 0.122 * x (Fig. 4). In contrast, the relationship between osteocyte lacuna density and body mass within the modern human sample is positive and described by r = 0.97, p < 0.01, r 2 = 0.99, y = 8.040 + 0.458 * x (regression not shown; see Table 4 for raw data).

Fig. 4
figure 4

Log osteocyte density linearly regressed against log body mass for primates ● and selected nonprimate mammals (♦) (least squares model). Raw data are provided in Table 4. When all primates and mammals examined here are treated together, the relationship is described by r = −0.85, p < 0.01, r 2 = 0.72, y = 10.500 – 0.122 * x. Summary statistics for primates are r = −0.84, = 0.01, r 2 = 0.71, y = 10.500 – 0.149 * x; for selected mammals, r = −0.99, = 0.01, r 2 = 0.99, y = 10.804 – 0.144 * x; and for the sample of five humans, r = −0.97, p < 0.01, r 2 = 0.99, y = 8.040 + 0.458 * x. The X-axis represents the log of body mass measured in kilograms

Striae of Retzius Repeat Interval Compares to the Lamellar Formation Rate

Primary Lamellar Bone

We imaged primary lamellar bone separately from secondarily remodeled lamellar bone tissue, recording the number of lamellae formed between vital labels when the number of days between the administration of vital labels was known or could be inferred (Table 5). In each case we find that one lamella is formed in the period matching the species-specific striae of Retzius repeat interval (e.g., Fig. 5, Table 5). Of note, patas monkey label fluorescence is saturated over three lamellae, with only a partial lamella following the last label, which is consistent with its repeat interval. Also, vital labels in the human (Fig. 5) are presumed to be 7- to 10-day treatment courses of tetracycline antibiotic; this duration would label only one lamella, which is consistent with the range of human repeat intervals.

Fig. 5
figure 5

Circularly polarized light image of primary lamellar bone from the midshaft femur of a growing macaque, Macaca nemestrina (FW = 237 μm; Table 5), with growth direction from right to left. Vital labels were administered at 2-week intervals, labeling the forming front with minocycline (yellow), xylenol orange, and, finally, DCAF (green). The animal was killed 18 days following the last label. The image of vital labels is truncated at the top and bottom to permit visualization of the association between label and lamellae

Secondarily Remodeled Lamellar Bone

Having confirmed that the lamellar formation rate in primary lamellar bone corresponds with the striae of Retzius repeat interval, we were also concerned to know if lamellae formed during a secondary remodeling event and contained within an osteon would provide the same result. Fluorescence microscopy of our prepared ‘ground’ thin sections of the midshaft femur of postadolescent humans revealed vital labels inferred to be therapeutic courses of tetracycline antibiotic. In these individuals, multiple labels could be identified both in primary lamellar bone laid down during growth and in secondary remodeled osteonal bone laid down during bone maintenance (e.g., Fig. 6). There are two important observations to make about the images of human lamellae presented in Fig. 6. First, there exists a mix of reasonably clear primary and secondary osteonal lamellae, which, in both tissue types around the 10th or 11th lamella, becomes more difficult to discern. This suggests that the physiological milieu responsible for lamella formation is similar between the two circumstances. Second is a caveat described under Materials and Methods, which is because of drift in the course of the forming secondary osteon, the cross section does not represent lamellae equally when counted on one side and then the other; this must be taken into account when enumerating lamellae.

Fig. 6
figure 6

Primary and secondary osteonal lamellae from the midshaft femur of a 38-year-old human female (FW = 237 μm for each image). The same numbers of lamellae are formed during bone growth and maintenance

We also observed primary lamellar and secondary osteonal lamellae in sheep labeled at known intervals (Fig. 7, Table 5), and as observed for human, the number of lamellae between labels is identical between the two lamellar tissue types.

Fig. 7
figure 7

Primary and secondary osteonal lamellae from the midshaft femur of a female sheep (FW = 237 μm for each image; Table 5). Two lamellae are contained within the 11-day labeling protocol, consistent with a 5-day striae of Retzius repeat interval. The same numbers of lamellae are formed during bone growth and maintenance

Discussion

Our aim has been to conduct a series of tests to determine if lamellar bone is an incremental tissue having relevance to interpretations of organismal life history, as has been widely recognized for enamel. The first test, using modern and fossil primates and proboscideans, was designed to test the relationship between enamel striae of Retzius repeat intervals and body mass. The second test, employing primates and a few nonprimate mammals, examined the relationship between rates of lamellar bone growth and body mass, where osteocyte density was used as a surrogate for rates of bone cell proliferation in lamellar bone. Having established these relationships with body mass, the conditions were set for testing the correspondence between enamel striae of Retzius and lamellar bone development.

Striae of Retzius Repeat Interval Correlates with Body Mass

Our analysis of the relationship between striae of Retzius repeat interval and body mass was performed on the largest and most taxonomically diverse sample of mammalian species to date. Significant correlations were found between these two variables, which hold true not only for the order Primates generally, including fossil hominins, but also for members of the Proboscidea, which include insular-dwarf forms. There is thus a clear association between mammalian long-period tooth growth rhythms and body mass, with a noted exception in one group of primates, the strepsirrhines (subfossil giant lemurs) having very low repeat intervals for their body mass [27]. When this is considered in light of our results on the relationship between bone and tooth development, this interesting exception may be explained (see below).

In addition to the interspecific comparisons examined here, the relationships between striae of Retzius repeat interval and body mass may also help to explain long-period variability within a species. Repeat intervals in Homo sapiens are highly variable, but modern humans also have an unusually high body mass variability for a primate. Preliminary research by us on six humans of known life history provides a statistically significant negative relationship between striae of Retzius repeat interval and body height (weight was extremely variable in the sample), suggesting that long period rhythms influence body size in humans (data not shown). Interspecifically, the relationship between repeat interval and body mass is positive, indicating that body mass increases are accomplished by slowing developmental rates over longer growth periods; i.e., larger-bodied species are characterized by relatively protracted life histories. Within a species and life history gamut, such as suggested by our preliminary research on humans, increases in body mass are at least partly accomplished by increasing the developmental rate within the same growth period and life history; i.e., the long-period rhythm is shortened. Long-period rhythm variation relative to body mass thus offers a potential, if not partial, solution to questions of life history variation, while remaining cognizant of other genetic factors responsible for growth factors operating within what are otherwise restricted genetic/historical programs [28, 29].

Osteocyte Lacuna Density Correlates with Body Mass

Based on the observed relationship between striae of Retzius repeat interval and body mass, and our appreciation for the scaling of bone mass to body mass, we advanced the working hypothesis that proliferation rates of osteoblasts, which become embedded in bone as osteocytes, vary with body mass. Osteocyte density is of special significance in this connection because, while osteoblast fates and functions are complex [30] (e.g., only a subpopulation are destined to become osteocytes), the positive relationship between osteoblast density and the proportion of osteoid surface covered by osteoblasts [31] permits specific predictions to be made concerning osteocyte density and life history. Osteocyte lacuna densities, as they reflect the rate of osteoblast proliferation, transformation, and incorporation into bone as osteocytes during growth, should be higher in organisms with rapid growth, with a small body mass, and whose osteoblast proliferation rates would lead to higher osteocyte lacuna densities; i.e., an inverse relationship between osteocyte density and body mass. No analysis has been performed to examine such a relationship within a mammalian order. We thus opted to test this hypothesis with a selection of very small to large primate and nonprimate mammal midshaft femurs (Table 4) and found that, indeed, osteocyte density is inversely proportional to body mass (Fig. 4). We ignored osteocyte apoptosis for the main reason that our interest in the development of bone mass as it relates to body mass concerns only the rate of osteoblast proliferation and the production of bone matrix by actively secreting cells. While osteocyte apoptosis most certainly occurs in some small percentage of cells, this is not a principal concern for analyses of the rate at which osteoblasts becomes embedded in the bone matrix; our interest here is independent, at least relatively so as far as we know presently, of the living bone cell population.

Mullender and colleagues [32] found a similar relationship for osteocyte lacuna densities within the cancellous bone tissues of the proximal femur in five mammals, speculating that osteocyte density is related to rates of osteoblast proliferation and obeys scaling laws relating to body mass, metabolic rate, life span, and growth rate. The well-recognized reduction in the capacity for cell proliferation during postreproductive aging and disease notwithstanding [33], the inverse relationship of osteocyte number per unit lamellar bone volume with body size thus appears to be an indicator of organismal life history; the faster the life history and skeletal development, the faster are osteoblast mitotic rates (the extent to which osteocyte density in woven bone tissue may be similarly interpreted is an area worth further investigation [34]).

Though the sample is small and the analysis preliminary, it is interesting that the relationship between osteocyte lacuna density and body mass in the modern human sample is positive. Interspecifically, the relationship is negative, indicating that body mass increases are accomplished by slowing developmental rates over longer growth periods; i.e., larger-bodied species are characterized by relatively protracted life histories. However, within the human species, body mass increases are accomplished by increasing the developmental rate. These inter- and intraspecific interpretations match those described for the relationship between repeat interval and body mass above.

Striae of Retzius Repeat Interval Correlates with Lamellar Bone Formation Rate

The bone lamellae of all mammals are strikingly similar in character, typical widths falling in the range of 3–6 μm, though thinner and much wider widths have been reported (to date there is no formal comparative mammalian study of lamellar widths). The lamellar design is one commonly associated with incremental patterns comprised of a layered structure, such as that represented by enamel cross striations, each layer being developed during one cycle of growth [25]. The temporal association of enamel formation being known, it was our aim to examine lamellar formation rates in mammals for which the times of administration of the labels were known or could be inferred (see Materials and Methods and Table 5).

In each case we found that one lamella is formed in the period matching the species-specific striae of Retzius repeat interval (e.g., Fig. 5, Table 5). For example, in the pigtailed macaque (Fig. 5), one vital label to the next was observed to contain between three and four lamellae, followed by four or five lamellae from the last label to the end of growth, which is consistent with a lamellar formation rate of 4 days/lamella and the striae of Retzius repeat interval for this taxon [8]. The lamellar structure in Fig. 5 contains some anisotropy, however, which is increased in the lower portion of the image. Concerning this anisotropy, there are two points to note. First, where one lamella is observed in the upper portion of this image, two may arise below. Second, the fluorescence image derives from a single confocal plane overlain onto a circularly polarized light image of the full section thickness, thus vital labels are observed to waver slightly around lamellae; i.e., this is a result arising because of our optical methods. Both of these phenomena are explained by the anisotropy discussed under Materials and Methods, which is informative and helps us to appreciate the three-dimensional complexity of forming bone surfaces.

In both humans and sheep it was also found that the lamellar formation rate was identical in both primary lamellar bone laid down during growth and secondary remodeled osteonal bone laid down during bone maintenance. Evidence for primary and secondary lamellar bone formation combined demonstrates that lamellar bone is an incremental tissue with a consistent time-dependent developmental morphology, similar to that described for enamel.

Long-Period Rhythms and Organismal Life History

Because the bone organ provides mechanical support for body mass, developmental rates for body mass should reflect developmental rates for bone mass. We thus hypothesize that the metabolic rhythm affecting osteoblast behavior and lamellar bone formation sets the pace of skeletal development; the effect this rhythm has on ameloblast function and enamel formation is secondary. Specifically, our morphological research supports the hypothesis that lamellar bone is an incremental tissue and that one lamella, mirrored by the stria of Retzius repeat interval in enamel, represents a centrally regulated long-period rhythm controlling bone and body mass. The functional significance of this is central to the evolution of vertebrate life histories. This hypothesized period, we suggest, is deeply integrated with the neuroendocrine functions of leptin in its control of bone and body mass, and most probably the role of osteocalcin in metabolism [3538] (see theoretical remarks below). The posited function of this long period clock is to regulate osteoblast proliferation and matrix secretion and whole-body metabolism [35, 36]. Provided that taxon-specific striae of Retzius repeat intervals are available or that long-period metabolic signatures will be identified in future research, bone lamellae may be calibrated in time. This will permit analyses of growth rate variability as a primary source of information about the duration and amplitude of growth, periodic phenomena, and other natural history events occurring during hard tissue formation. Lamellar bone can be evaluated over known time intervals for the durations of experimental protocols or in studies of metabolic dysfunction or disease affecting hard tissue formation.

Information about biological periodicities, environmental and/or physiological cycles, and perturbations to the developing organism are all potentially contained in growth rate variability studies of lamellar incremental patterns. For instance, given that lamellae are formed within defined periods of time, quantitative measures of widths of individual lamellae may provide time-resolved growth rate variability and reveal rhythms heretofore unknown. An image of fossil hominin (Paranthropus robustus) lamellar bone, calibrated at 7 days/lamella, which is the striae of Retzius repeat interval for this taxon (Table 2), reveals a long-period rhythm resembling the lunar cycle, not previously observed in bone (Fig. 8).

Fig. 8
figure 8

Lamellar growth rate variability of Paranthropus robustus metacarpal (SK 3646) lamellae is based on a method developed for assessment of growth rate of anisotropic incremental structures [25]. The Z-axis represents lamellar width, or growth rate in micrometers, according to the color-coded scale; the X-axis is a relative time scale of the growth rate in two dimensions, where X = 0 is at the beginning and X = 100 is at the end of the lamellar growth sequence; the Y-axis denotes different levels of structural lamellar detail, where Y = 0 refers to the most basic elements of the lamellar structure and Y = 9 means that all structural elements associated with nine measurement transects were taken into account. Based on a 7- day period for Paranthropus robustus (Table 2) and a number of lamellae near to 35 in the original image, lamellar growth rate illustrates a periodicity of roughly 4 weeks

As noted above, the subfossil lemur, Palaeopropithecus ingens, is an exception, having a lower repeat interval for its body size than observed for primates generally. Exceptions of this nature have the potential of offering wonderful insight into life history processes, and in this case, the correspondence between repeat intervals and lamellar formation enables us to reconcile the enigmatic relationship between their rapid dental development and their estimated large body size of 45 kg. In the absence of any knowledge of the incremental nature of their skeletal development, one plausible interpretation would be that the pace of dental and somatic growth was decoupled, allowing the dentition to form quickly in relation to a slower-growing body. However, we hypothesize that dental and skeletal development are in a manner coupled, and that the very short, 2-day, striae of Retzius repeat interval in these primates relates specifically to their large size, achieved by retaining a historically constrained rapid lamellar formation rate over a somewhat prolonged period of somatic development.

The hypothesized long-period rhythm is thus worthy of future study both as a potential analytical tool and as a paradigm for analyses of organismal life history variation and evolution. Quite apart from its effects on hard tissues, the hypothesized rhythm is also related to the overall pace of development and life at high taxonomic levels. We suggest that this long period likely regulates size and rate variation, which coordinate the organismal life history gamut [39]. The lower height of the primate osteocyte production function (Fig. 4), for instance, agrees with growth law predictions indicating that primate production energy is significantly less compared with other mammals of the same body size [40]. These predictions are part of a larger model, which explains that growth energy determines also how much energy is available for reproduction later. This suggests that production rate, bone and body size, and reproduction are centrally linked to this long-period rhythm, which is a mediating mechanism in the regulation of organismal life history according to metabolic constraints [41]. These mechanisms may also regulate the cell proliferation rates of some soft tissues and thus of increases in body organ masses during postnatal development. Cell division rates in vivo appear to regulate longevity [42] and thus the integrative study of long period rhythms and local aging mechanisms (e.g., neuroendocrine imbalance, DNA repair and metabolism, and related genetic disorders) is expected to provide insight on their combined influence on longevity and the aging process. We expect future neuroendocrine and physiological research to identify the role of long-period developmental rhythms and address such suggestions.

A Hypothesis for Central Control of Long-Period Rhythms and Bone/Body Mass in Mammals

Research on biological rhythms in mammals has revealed a suite of well-known circadian oscillations that regulate fundamental day-to-day metabolic requirements, growth, and maintenance of the organism [1]. However, concentrations on circadian biological rhythms, particularly those established using “gold standards” based on small-bodied mice, rats, or hamsters, whose rhythmic growth is fast, would miss fundamental long-period physiological rhythms regulating developmental rates in organisms having slower and longer life histories. Slow-growing and large-bodied organisms have protracted life histories that suggest the possibility of long-period physiological and developmental regulation. A long-period rhythm might explain how the slow growth rate and protracted life history of a large-bodied organism can be achieved by its cells proliferating over a longer period than that of a circadian rhythm.

One long-period rhythm has been noticed; a circaseptan (near-7-day) periodicity in heart rate and pressure has been demonstrated in humans [4345]. Furthermore, it has been specifically suggested that this centrally regulated autonomic, sympathetic-mediated heart rate and pressure oscillation in humans relates to striae of Retzius repeat intervals [43], which have a modal frequency of 8 days [13]. We conjecture that the resemblance between circaseptan variability in sympathetic drive (tone) and striae of Retzius formation in humans is no coincidence, and that they are physiological manifestations of a centrally regulated biological rhythm. Such periodic rhythms, we believe, are likely to be driven by an oscillator rather than an hourglass mechanism, the latter of which, for instance, explain single life events in an organism (e.g., menopause) [46].

As described above, the relationship between osteocyte lacuna density and body mass suggests that many cells, a proportion of which are subsequently incorporated as osteocytes, are required to form a single lamella during a relatively brief period in a small mammal, but that fewer mitoses would be required in the development of a large-bodied organism, their osteoblasts forming one lamella over a longer period of time. An important body of work has recently demonstrated that leptin-bound receptors in the hypothalamus mediate sympathetic inhibition of osteoblast proliferation through the β2-adrenergic receptor, thus exerting control of bone mass [3638, 47]. This research has shown that leptin-mediated control of bone mass is complementary to leptin’s other hypothalamic actions relating to appetite, satiety, and regulation of body mass. It has also been shown that integrity of the peripheral (i.e., osteoblast population) molecular clock is necessary for this regulatory system to operate. A necessary coadaptive feature of body mass variation is its relationship with bone mass: i.e., a bone mass consistent with the support of body mass. Takeda and colleagues [48] have highlighted the importance of several discrete hypothalamic nuclei for integrated leptin-mediated autonomic control of bone and body mass and reproduction.

Recently the endoluminal release of serotonin by enterochromaffin cells of the duodenum has also been demonstrated to regulate (inhibit) osteoblast proliferation and bone formation through the cAMP response element binding (CREB) transcription mediator [49]. CREB is as well implicated in sympathetic inhibition of osteoblast proliferation through the β2-adrenergic receptor, as mentioned above. While the parasympathetic nervous system may be one link in the autonomic regulation of bone mass via vagal pathways, efferent stimulation of splanchnic nerves of the sympathetic nervous system also cause release of duodenal serotonin through a β2-adrenergic mechanism [50], raising the possibility of a two-pronged central sympathetic regulatory system for controlling bone mass.

A hypothesized autonomic long-period rhythm controlling osteoblast proliferation and striae of Retzius repeat intervals agrees with data presented on mammalian osteocyte density. Thus we believe that a strong circumstantial case exists for central autonomic control of long-period rhythms over aspects of life history such as body mass. We term this hypothesized period the Havers-Halberg Oscillation, with reference to Clopton Havers [51], the first to observe and describe both the lamellae in bone and the striae of Retzius in enamel, and Franz Halberg, a long-time explorer of long-period rhythms [52, 53].

We wish to add, however, that it remains to be determined whether there are circumstances concerning continuously forming bone surfaces in which the correspondence between striae of Retzius repeat intervals and the formation of lamellae falters. In a 38-year-old (Fig. 6) and a 41-year-old human female (data not shown), both primary and secondary osteonal lamellar tissues were simultaneously labeled during a time in these women’s lives when they were still growing, the two bone tissue types in each individual exhibiting an identical number of lamellae between labels. However, while we believe it to be unlikely that a centrally regulated long-period rhythm would change with age, we have not tested for its maintenance in the aged. Further, it must be determined whether there are any pathological conditions that might alter long-period and lamellar formation rhythms, operable through the SCN and/or in combination with other hypothalamic nuclei and peripheral clocks.

In sum, and in wait for further experimental testing, we propose that vertebrate long-period rhythms are a function of natural selection operating on a hypothalamic autonomic clock responsible for the pace of development and bone/body mass, the study of which we believe holds much promise for the integrative science of organismal life history.