Osteohistological variation in growth marks and osteocyte lacunar density in a theropod dinosaur (Coelurosauria: Ornithomimidae)
BMC Evolutionary Biology volume 14, Article number: 231 (2014)
Osteohistological examinations of fossil vertebrates have utilized a number of proxies, such as counts and spacing of lines of arrested growth (LAGs) and osteocyte lacunar densities (OLD), in order to make inferences related to skeletochronology and mass-specific growth rates. However, many of these studies rely on samplings of isolated bones from single individuals. These analyses do not take individual variation into account, and as a result may lead to misleading inferences of the physiology of extinct organisms. This study uses a multi-element, multi-individual sampling of ornithomimid dinosaurs to test the amount of individual variation in the aforementioned osteohistological indicators. Based on these results we also assess the conclusions of previous studies that tested paleohistological hypotheses using isolated elements.
LAG number was found to be consistent within the hind limb bones of each individual, with the exception of the fibula, which preserves one additional LAG. Considerable differences in LAG spacing were found between elements of the sampled individuals, with larger variation found in elements of the foot compared with the femur, fibula, and tibia. Osteocyte lacunar density ranged between 29000 and 42000 osteocyte lacunae per mm3, and was found to vary more between hind limb bones of an individual and within bones, than between the average values of individuals.
The variation between hind limb elements in LAG number and LAG spacing suggests that direct comparisons of these elements may be misleading, and that LAG spacing is not a reliable proxy for mass-specific growth rates of an individual. Sampling of multiple bones should be performed as an internal check of model-based LAG retro-calculation and growth equations. The observation that osteocyte lacunar density varies more between individual bone elements than between average individual values suggests that the choice of sampled element can greatly influence the result, and care should be taken to not bias interpretations of the physiology of fossil tetrapods.
Osteohistological examinations of fossil vertebrates have provided important insights into their ontogeny and physiology -. In addition to the description of bone tissue types, considerable attention has been paid to skeletochronological indicators, such as the number and relative spacing of lines of arrested growth (LAGs) and other cyclical growth marks within the bone cortex -. The presence of tightly spaced growth marks within relatively avascular tissues in the outer margin of the bone cortex (the external fundamental system, or EFS, of some authors, e.g. Horner et al. , or the outer circumferential layer, or OCL, following Cormack ) has been used to infer somatic maturity in sampled individuals ,,. In the absence of an EFS, patterns in the relative spacing of LAGs has been used to infer the ontogenetic stage of an individual, in some cases using a single sectioned element ,,-. On a finer scale, the density of osteocyte lacunae (osteocyte lacunar density, OLD) in the cortex of limb bones of tetrapods has been suggested to be related to mechanical forces and differential loadings, basal metabolic and growth rates, and element structural differences -, and osteocyte lacuna morphology itself has also been shown to be related to these properties ,,. Recently, the densities of osteocyte lacunae have been used to predict basal metabolic rates and relative growth rates in fossil tetrapods, including dinosaurs .
Few studies to date have focused on variation in LAG spacing between different limb bones of a single individual; individual variation in osteocyte lacunar morphology , and density , in the long bones of tetrapods suggests that element choice and sampling location may influence studies using these histological parameters to infer life history traits in fossil taxa. Here, we describe bone microstructure and assess intra- and inter-skeletal variation in skeletochronological indicators (LAGs) and osteocyte lacunar density (OLD) in a sample from the Horseshoe Canyon Formation (Maastrichtian, Alberta) of ornithomimid theropod skeletons, three of which were found closely associated in the same stratigraphic horizon and likely represent individuals dervied from the same population . Specifically, we use the intra-specific, multi-element osteohistological data to test two separate hypotheses: 1) that patterns of lines of arrested growth (LAG) spacing in different elements of a single individual provides a consistent signal for evaluating its relative maturity, and, 2) osteocyte lacunar densities (OLDs) are generally consistent between limb elements of a single individual, and between individuals of a given taxon due to their shared physiology and loading regimes. Growth marks, such as LAGs, are particularly well documented in theropod dinosaurs ,,,,-, and this detailed case study of individual variation in a single theropod taxon may provide insights that can be extended to other theropod dinosaurs, as well as tetrapods more generally.
In each sampled element of the bonebed individuals (Figure 1), the bone matrix is consistent in being predominantly a woven-parallel complex, showing a combination of laminar, plexiform, and reticular patterns of vascularization. An external fundamental system (EFS) is absent in all elements, with all periosteal margins showing active primary bone deposition and well-vascularized tissue ,,. However, some differences exist in the relative distributions of vascularization and the degree of secondary remodelling in elements of each individual, and in the number of lines of arrested growth (LAGs) preserved.
The bone tissue of the femur of the largest individual (CMN 12068) (Figure 2A) ranges from laminar to plexiform, with the density of radial canals varying throughout the section. Secondary remodelling is minimal in most areas of the femur; however, there is a heavily remodelled section of the medial cortex that is associated with a slight external eminence of the periosteal margin where the cortex is thickened relative to the rest of the section. Three LAGs are visible in the femur, with the innermost LAG partially lost due to remodelling and cortical drift. The outermost LAG appears as a double LAG in some limited sections of the cortex. The tibia of CMN 12068 (Figure 2B) is comprised primarily of laminar bone, with some plexiform regions due to localized increases in the density of radial canals. Secondary remodelling occurs in the inner cortex, although secondary osteon development is not extensive, nor is there evidence of considerable outward expansion of the medullary cavity along the endosteal margin. Cortical thickness is inconsistent, with the medial side of the element being considerably thicker. Three LAGs that completely extend around the cortex are visible in the tibia. The fibula (Figure 2C) of this individual shows plexiform and reticular bone, with greater vascularization towards the periosteal margin. There is extensive secondary remodelling on the medial side of the cortex, including several large, well-developed, secondary osteons, and a poorly developed medullary cavity. Unlike the other sectioned elements of this individual, the fibula preserves four visible LAGs. The third pedal phalanx (Figure 2D) exhibits primarily laminar vascularization, with some regions appearing semi-plexiform to plexiform. There is secondary remodelling concentrated along the inner cortex, particularly along the ventral side of the phalanx. This area of the cortex is also slightly thicker than the lateral or dorsal sides. Like the femur and tibia, three LAGs are preserved in the phalanx.
CMN 12069 (Figure 3A) is smaller than CMN 12068, but approximately the same size as CMN 12070. The tibia of CMN 12069 has primarily laminar to plexiform vascularization, with areas of reticular bone associated with increased density of radial canals. There is little to no secondary remodelling, nor is there evidence of erosion of the cortex via expansions of the medullary cavity. The cortex thickness is relatively consistent in the incompletely preserved section, and two LAGs are present. The fibula (Figure 3B) of this individual shows primarily reticular vascularization, with more extensive secondary remodelling on the medial side of the cortex. However, secondary remodelling is less extensive than in the larger fibula of CMN 12068. A well-defined medullary cavity is present, and the fibula preserves three LAGs. Metatarsal IV (Figure 3C) has reticular and plexiform vascularization, concentrated towards the inner cortex, and laminar vascularization concentrated towards the outer cortex. Much of the inner cortex is secondarily remodelled, although there are no outward expansions of the medullary cavity along the endosteal margin of the cortex. As in the tibia, but unlike the fibula, there are two LAGs present in the metatarsal. The pedal phalanx (Figure 3D) exhibits primarily laminar vascularization, with two LAGs present in the cortex.
The tibia of CMN 12070 (Figure 4A) has woven-parallel complex bone tissue with primarily laminar to plexiform vascularization, with little to no secondary remodelling . The cortex thickness appears relatively consistent, and two LAGs that extend around the cortex are visible. Like CMN 12069, the fourth metatarsal (Figure 4B) is comprised of plexiform and reticular vascularization concentrated towards the inner cortex, and laminar vascularization concentrated towards the outer cortex. There is secondary remodelling concentrated in the inner cortex and through sections of the medial and lateral outer cortex As in its associated tibia, there are two LAGs visible in the metatarsal of CMN 12070.
The fibula of ROM 852 (Figure 5) exhibits primarily plexiform vascularization with reticular vascularization in some sections towards the inner cortex, and less vascularized, more laminar tissue towards the periosteum. There is no definitive presence of an EFS/OCL. Secondary remodelling is relatively extensive on the medial side of the cortex. Five LAGs are present (indicated in Figure 5 and contained inset [I]), with inter-LAG spacing increasing between the first and second LAG, and decreasing considerably between the subsequent LAGs approaching the periosteum. As in CMN 12068, this specimen has double LAGs, most notably near the middle of the cortex in the second LAG, and additionally in variable positions of the fifth LAG. While most of the LAGs visibly extend around the cortex (the first and second most obviously, the fourth and fifth visible, but less distinct), the third LAG is visible in the posterior cortex, but becomes indistinct as it grades into the more heavily vascularized mid cortex and more remodelled anterior cortex.
Variability in growth line spacing
Among the three ornithomimids from the Horseshoe Canyon Formation bonebed (Figure 6A), there is a similarity in the pattern of LAG spacing in the femur, fibula, and tibia, with all three elements showing a relative increase in spacing from the periosteal margin to the innermost LAG. However, the lowermost hind limb elements (MT IV, pedal phalanx) did not show the same signal within a given individual; the metatarsals and pedal phalanges show a more consistent spacing between LAGs, and do not exhibit a decrease between the innermost LAG and endosteal margin. Although the absolute distance between LAGs varies within individual bones, the relative pattern of LAG spacing (increases versus decreases) within a single bone is consistent between different measurement transects (Figure 5). However, the relative proportion of cortex defined by each cycle can vary in different measurement transects.
Osteocyte lacunar density
Osteocyte lacunar density (OLD) in the total sample of four ornithomimids from the Horseshoe Canyon Formation ranged from ~29000 to ~42800 osteocyte lacunae (OL) per mm3; standard deviations (SD) ranged from ~600 to ~6000 OL/mm3, depending on the individual or element sampled (Table 1). The individual specimens averaged 34955 OL/mm3 (SD of 1855 OL/mm3) for CMN 12070, 35002 OL/mm3 (SD of 5943 OL/mm3) for CMN 12069, 37176 OL/mm3 (SD of 2982 OL/mm3) for CMN 12068, and 37314 OL/mm3 (SD of 1408 OL/mm3) for ROM 852; the global individual average OLD is ~36112 OL/mm3 (SD of 1310 OL/mm3). Sampling of the medial and lateral sides of each sectioned bone of CMN 12068 yielded 38681 OL/mm3 and 41533 OL/mm3 for the femur, (mean of 40107 OL/mm3, SD of 2017 OL/mm3) 34052 OL/mm3 and 33180 OL/mm3 for the tibia, (mean of 33616 OL/mm3, SD of 617 OL/mm3) 36840 OL/mm3 and 40000 OL/mm3 for the fibula, (mean of 38420 OL/mm3, SD of 2234 OL/mm3) and 34769 OL/mm3 and 38354 OL/mm3 for the pedal phalanx (mean of 33561 OL/mm3, SD of 2535 OL/mm3). Sampling of the inner, middle, and outer cortex of the fibula of ROM 852, yielded results of ~36004 OL/mm3 (SD of 1239 OL/mm3), ~37137 OL/mm3 (SD of 1202 OL/mm3), and ~38802 OL/mm3 (SD of 1508 OL/mm3) respectively.
For an estimated ornithomimid body size range of 90–125 kg , the minimum and maximum, average OLD values of the Horseshoe Canyon Formation ornithomimids plot midway between birds/mammals and non-avian theropod dinosaurs , but show considerable variation when individual bone measurements are plotted (Figure 7). This variation is illustrated further through a comparison of the range of OLD values for each element of the three bonebed individuals against the mean OLD values for the species derived from the individual averages (Figure 8A), a comparison of the variation in OLD within (via medial and lateral sampling) and between elements of CMN 12068 (Figure 8B), and a detailed comparison of the OLD values across a transect through the inner, middle, and outer cortex of the fibula of ROM 852 (Figure 8C). The latter comparison was also assessed via a series of ANOVAs, which found no significant differences in the OLD between the inner and middle cortex (p = 0.79), between the inner and outer cortex (p = 0.06), or between the middle and outer cortex (p = 0.97).
All examined ornithomimid sections are characterized by woven-parallel complex tissue, with a range of laminar to reticular vascularization (primarily plexiform), and lack evidence of an EFS/OCL . Each element sectioned has at least two LAGs. Within each bonebed individual, there is consistency in the number of LAGs within each element, with the exception of the fibula, which consistently has one more LAG. Additionally, localized double LAGs are present in the femur of CMN 12068 and the fibula of ROM 852, which suggests differential growth and localized suspension of bone deposition . The femur of CMN 12068, the fibulae of CMN 12068, CMN 12069, and ROM 852, as well as the fourth metatarsal of CMN 12070, show localized secondary remodelling across the cortex that is likely related to adductor muscle attachment . CMN 12069 and CMN 12070, the smaller individuals within the sample (metatarsal length ~29.3 cm & ~29.0 cm, respectively), show consistency in the number of LAGs in corresponding elements, which supports the hypothesis that this bonebed represents a mass death assemblage where the preserved individuals are part of the same cohort . In the largest bonebed individual, all corresponding elements preserve one additional LAG compared to those of the smaller individuals, suggesting that it is at least one year older ,,,. ROM 852, the largest individual overall (Additional file 1), had one more LAG in the fibula than CMN 12068. Although there is no evidence of an EFS/OCL in ROM 852, there is a decrease in tissue vascularity in the outermost cortex that suggests a reduction in bone growth rate ,. As no EFS/OCL is preserved in any of the sectioned ornithomimid elements from this study, none of these individuals had yet reached skeletal maturity, and were likely still growing in body size. These specimens are of similar overall size to many other known North American ornithomimids (Additional file 1), including TMP 1995.110.1, which is reported to have an EFS/OCL . This suggests that the maximum size of North American ornithomimids is not yet fully understood, and may show a large degree of individual or taxonomic variation.
This multi-element, multiple individual sampling approach permits the first in-depth, quantitative assessment of individual variation in osteohistological traits of a theropod dinosaur species. Differences in LAG counts in bones from the same individual is a pattern observed in other dinosaur histology studies ,,,, and is not unexpected due the possibility of different rates of bone growth and remodelling in different parts of the skeleton. However, in this study of hind limb elements, only the fibula differs from the other sampled elements in that it preserves one additional LAG among the multiple ornithomimid individuals. This differs from other studies that have sampled a wider representation of bones within the skeleton and have shown more variation within multiple elements of single individual ,,,. The presence of an additional LAG in the fibula suggests that it may be a preferred element for skeletochronology in theropods, because it preserves the longest growth record due to a small medullary cavity and less remodelling, as suggested in previous studies ,. A potential drawback of using the fibula is that its irregular shape, coupled with strong allometric growth in diaphysis shape (e.g., Figure 5), complicates retro-calculation-based age reconstructions and growth-curve model-fitting. Femora and tibiae do not have these growth related problems thus making them more favourable elements in equation-based retro-calculation growth reconstructions . Therefore, the diaphysis of these primary weight-bearing bones may be more effective for reconstructing growth curves than fibulae, despite the fact that the fibula appears to consistently preserve more LAGs in theropods. Therefore, we propose that a within-individual, multi-element sampling approach can act as an internal check of retro-calculation-based age estimations, and that this will be the most powerful means of reconstructing growth curves of extinct vertebrates going forward.
Spacing between lines of arrested growth is often used to qualitatively assess the relative growth rates and maturity in isolated bones, e.g. Tyrannosaurus ,, Alioramus , Raptorex ,,, and a putative Alaskan ornithomimosaur . This method assumes that a decrease in LAG spacing from the inner cortex to the periosteum corresponds with an individual approaching somatic maturity. In our sample, LAG spacing is variable between different elements within a single individual (Figure 6A). The femur, tibia, and fibula do show a similar trend of decreasing distance between LAGs, regardless of where within the cortex the spacing is measured. However, the metatarsal and pedal phalanges do not show the same pattern as their corresponding upper limb bones. Here, distances between LAGs in the outer cortex, which presumably represent the same periods of growth as record in the crus and femur, are more variable and do not show constant decreases in spacing towards the periosteal surface. Therefore, hind limb and foot bones of a single individual would yield different assessments of relative maturity if analysed in isolation. For example, the inference that an Alaskan ornithomimosaur  was approaching somatic maturity at the time of death based on a decreasing pattern of LAG spacing in a single, isolated metatarsal, is poorly substantiated.
This type of variability in LAG spacing between elements of a single individual is also evident in other fossil studies (Figure 6B). Differing LAG spacing signals also appear to be present within different elements of an individual adult Hypacrosaurus stebingeri (MOR 549), although the presence of an EFS/OCL makes this difficult to show definitively; the fibula shows a pronounced increase in LAG spacing from the periosteal margin to the inner cortex, whereas the rib has more variable spacing between LAGs. In addition, alternating patterns of increase vs. decrease in LAG spacing through the cortex are also present in both the tibia and fibula of the Asian ornithomimid Sinornithomimus dongi , as well as in a femur of Tenontosaurus . Variability in LAG spacing from the inner cortex to the periosteum is consistent with recent extant experimental work that shows that different long bones of an individual grow at different rates , and may also relate to annual variation in resource availiblity ,, or functional responses to ontogenetic gait shifts , among other factors. These problems could be further compounded if spacing is not assessed at comparable locations in the cortex, as they could lead to different interpretations of growth or ontogenetic status (such as in Raptorex, which has been suggested to represent either a small-bodied nearly-mature animal  or a juvenile Tarbosaurus , with both arguments made partly using LAG spacing). This demonstrates the unreliability of simple LAG spacing patterns for inferring a particular ontogenetic stage or relative maturity across a wide range of taxa or elements in somatically immature specimens where an EFS/OCL is not present ,.
Osteocyte lacunar density
Recently, Stein and Werner  quantified and compared osteocyte lacunar density (OLD) in a sample of adult limb elements from a variety of extant and extinct taxa, which included mammal OLD data from a study by Bromage et al. . They found significant relationships between OLD and body mass, and body mass-specific OLD and relative growth rates. Whereas the Stein and Werner study did not explicitly address how intra- and/or inter-individual variation in OLD may limit the interpretive power of their results, care was taken to sample from the outer cortex of each element, sample the femur when available (68% of sampled elements in the sample), and to use the largest individual per taxon in regression analyses when more than one sample was taken (as opposed to a mean taxon value), in order to avoid variation in OLD values resulting from ontogeny or sampling location. Although 41 percent of taxa included in Stein and Werner  were represented by more than one individual or sample (e.g., cf. Trilophosaurus, Hyperodapedon, Rhamphorhynchus, Albertosaurus, and Tyrannosaurus, in which more than two individuals were sampled), intra-specific, intra-individual, and intra-element variation in OLD was not explicitly quantified or discussed in context of their results. Given the variation in lacunar morphology (size and shape) between elements ,,, variation in OLD within individual human femora  and between fore-limb elements of modern taxa , as well as quantification issues to due sectioning plane , intra-specific variation is potentially problematic when trying to predict aspects of organismal palaeobiology using osteocyte lacunar density.
Measured ornithomimid osteocyte lacunar densities varied within individual elements, and between elements of the same individual (Table 1). All OLD values calculated for each individual are larger than those of larger tyrannosaurid theropods sampled by Stein and Werner , as predicted by the inverse relationship between OLD and body mass. OLD values (all raw counts and global average) are relatively larger than similarly sized mammals, but plot between mammals and birds of similar body size  (Figure 7). This result is not unexpected, given the status of ornithomimids as small- to medium-bodied theropods ,. However, our analysis indicates that there is considerably more variation in OLD between elements of an individual than in different positions within a single element (Figure 8B&C), or in the average OLD between different individuals (Figure 8A). This level of variation between elements of single individuals suggests that estimating metabolic , or genome size  characteristics of extinct organisms may be more complex than previously considered, and OLD counts may be affected by individual variation and allometry of different skeletal elements, as noted in studies of extant mammals -,. Despite this high intra-skeletal variation, the mean OLDs of each specimen, calculated from the lacunar densities of multiple sampled elements in its skeleton, are very similar, despite each specimen not always being sampled from the same set of elements. Given this consistency, multi-element sampling and the use of individual/species mean OLDs is recommended in order to capture and control for the high variation present in sampling isolated individual elements.
Stein and Werner  postulated a relationship between ontogenetic stage and OLD, and suggested sampling of OLD in the outer cortex of mature animals to control for changing OLD throughout ontogeny. Our data suggest that, when sampling element and location is kept consistent, ontogenetic variation in OLD appears negligible. Over the ontogenetic size range investigated here (Figure 8C), no significant differences were found in sampling sites in the inner, middle, or outer cortex. However, it should be noted that none of our sampled series represents a skeletally-mature animal, given the lack of EFS/OCL, and there appeared to be a weak, albeit non-significant, positive trend towards increasing OLD through growth.
Histological thin sectioning of multiple hind limb elements of three semi-articulated ornithomimid individuals, and a fibular thin section from a larger, stratigraphically contemporaneous specimen of Ornithomimus edmontonicus show a consistent pattern in LAG numbers. Each sectioned element within an individual has the same LAG number, with the exception of the fibula, which possesses one additional LAG. None of the sampled individuals appear to be skeletally mature, despite being of similar size to previously published mature individuals, suggesting a greater degree of individual variation in growth within North American ornithomimids than has been previously reported. Spacing of LAGs shows two distinct patterns, with femur/tibia/fibula LAG spacing decreasing relatively consistently from the inner to outer cortex, while the metatarsal/phalanx spacing remains more variable in spacing throughout the cortex. LAG spacing appears to be variable, and too easily affected by other factors of bone growth and modification, to be of repeatable use in assessing growth throughout ontogeny of individual animals or species (not counting the presence of an EFS/OCL). Similarly, osteocyte lacunar density appears more highly variable when making comparisons between bones of a single individual than when comparing the average values of different individuals. This suggests that, at least within closely related groups, element choice may have a greater impact on the measured result than ontogenetic stage of the individual. Localized factors relating to differential limb growth may play a larger role in the distribution and density of osteocyte lacunae than previously thought, and are potentially biasing the interpretations of analyses based on single samplings of isolated specimens.
Three articulated partial ornithomimid skeletons from the Horseshoe Canyon Formation (HCF) were histologically sampled with permission from the Canadian Museum of Nature (CMN) for this study. The skeletal completeness, anatomy, and taphonomy of CMN 12068, 12069, and 12070 (cf. Ornithomimus edmontonicus, originally Dromiceiomimus brevitertius) was described by Cullen et al. . One additional specimen sampled from the HCF was Royal Ontario Museum (ROM) 852 (O. edmontonicus, originally Struthiomimus ingens). The elements sampled were: CMN 12068 (femur, tibia, fibula, phalanx III-3), CMN 12069 (tibia, fibula, metatarsal IV, phalanx III-3), CMN 12070 (tibia, metatarsal IV), ROM 852 (fibula).
All thin-sections were created at the Royal Ontario Museum. Samples were initially embedded in Castolite AP polymer resin and cut as transverse diaphyseal sections using a Buehler Isomet 1000 wafer blade low-speed saw. These sections were then mounted to glass (CMN 12068, 12069, 12070) or plexiglass (ROM 852) slides using Scotch-Weld SF-100 cyanoacrylate and ground down to approximately 50–80 μm thick (for specific section thicknesses, see Table 1) using a Hilquist grinding cup, being ground progressively finer by hand using silicon carbide powder. Thin-sections were examined and photographed under plane-polarized light using a Nikon DS-Fi1 camera mounted to a Nikon AZ-100 microscope fitted with cross-polarizing and lambda filters. Specimen images were processed using Nikon NIS-Elements (Basic Research) v. 3.13.
Spacing between LAGs (=zone thickness) was measured in micrometres for each sectioned element of CMN 12068, CMN 12069, and CMN 12070. Measurements of LAG spacing were taken from the periosteal margin to the first visible LAG, then between each LAG, and from the innermost visible LAG to the endosteal margin of the cortex. Measurements were taken where the cortex was the thickest without extensive remodelling, or in the case of the tibia and metatarsal of CMN 12069, where the section of the cortex was most complete. This was done in an attempt to maximize comparability between sections, while minimizing subjectivity in space differences due to intra-element allometry. These values were then converted into percentages of the total cortex to allow for comparisons of relative spacing between elements. This was repeated for previously published histological sections of Hypacrosaurus , Tenontosaurus , Sinornithomimus , Tyrannosaurus , and Raptorex , and plotted in Figure 6B.
Osteocyte lacunar densities (OLD) were measured following the method (with slight modification) of . Sections were re-imaged at approximately 32× magnification, and a 250 μm × 250 μm area of the outer cortex was selected from the resulting images. Images were taken at each focal plane of this square and combined into a single image via z-stacking. In the selected area of this z-stacked image, osteocyte lacunae of each element were counted three times and averaged. The lacunae counts were then converted to lacunar density per 1 mm3 of cortex (Table 1). Additional density counts per individual were obtained through averaging of individual element values. Each sectioned bone of CMN 12068 was resampled at locations on opposite sides of the outer cortex to additionally assess variation within an individual element. Osteocyte lacunar densities of ROM 852 were obtained from positions in the inner, middle, and outer bone cortex. Each position was sampled five times, and the inner/middle/outer measurements were compared using ANOVA in the ‘stats’ package in R . Osteocyte lacunar measurements were plotted with those of Stein & Werner  in order to further illustrate the range of individual and inter-element variation in a single taxon of theropod (Figure 7), with the body mass estimated for the largest ornithomimid specimen (ROM 852) using the ‘MASSTIMATE’ package in R and the methods described in Campione et al. , with body masses of the other sub-adult individuals estimated using Developmental Mass Extrapolation  based on known/estimated femur lengths (Additional file 1). Osteocyte lacunar densities for sample, and the resulting means and standard deviations were also plotted as box plots in Figure 8, using the ‘ggplot2’ package in R .
Availability of supporting data
Full-section images and images used for OLD calculations are available in the MorphoBank repository, project 1218 at http://morphobank.org.
TMC, DCE, and MJR conceived of the project. TMC performed the thin-sectioning, imaged the resulting slides, and created the figures. TMC, DCE, and MJR conducted the general histological examination. TMC and DCE performed the LAG spacing and lacunar density analyses. TMC performed all statistical tests. TMC, DCE, MJR, PJC, and YK interpreted the results and contributed to the writing of the manuscript. All authors read and approved the final manuscript.
Burns ME, Vickaryous MK, Currie PJ: Histological variability in fossil and recent alligatoroid osteoderms: systematic and functional implications. J Morphol. 2013, 274 (6): 676-686. 10.1002/jmor.20125.
Bybee PJ, Lee AH, Lamm ET: Sizing the Jurassic theropod dinosaur Allosaurus: assessing growth strategy and evolution of ontogenetic scaling of limbs. J Morphol. 2006, 267 (3): 347-359. 10.1002/jmor.10406.
de Ricqlès A, Padian K, Knoll F, Horner JR: On the origin of high growth rates in archosaurs and their ancient relatives: complementary histological studies on Triassic archosauriforms and the problem of a “phylogenetic signal” in bone histology. Ann Paléontologie. 2008, 94 (2): 57-76. 10.1016/j.annpal.2008.03.002.
Erickson GM, Currie PJ, Inouye BD, Winn AA: A revised life table and survivorship curve for Albertosaurus sarcophagus based on the Dry Island mass death assemblage. Can J Earth Sci. 2010, 47 (9): 1269-1275. 10.1139/E10-051.
Erickson GM, Makovicky PJ, Currie PJ, Norell MA, Yerby SA, Brochu CA: Gigantism and comparative life-history parameters of tyrannosaurid dinosaurs. Nature. 2004, 430 (7001): 772-775. 10.1038/nature02699.
Horner JR, De Ricqlès A, Padian K: Long bone histology of the hadrosaurid dinosaur Maiasaura peeblesorum: growth dynamics and physiology based on an ontogenetic series of skeletal elements. J Vertebr Paleontol. 2000, 20 (1): 115-129. 10.1671/0272-4634(2000)020[0115:LBHOTH]2.0.CO;2.
Horner JR, Padian K: Age and growth dynamics of Tyrannosaurus rex . Proc Biol Sci. 2004, 271 (1551): 1875-1880. 10.1098/rspb.2004.2829.
Klein N, Sander M, Suteethorn V: Bone histology and its implications for the life history and growth of the Early Cretaceous titanosaur Phuwiangosaurus sirindhornae . Geol Soc Lond, Spec Publ. 2009, 315 (1): 217-228. 10.1144/SP315.15.
Kohler M, Marin-Moratalla N, Jordana X, Aanes R: Seasonal bone growth and physiology in endotherms shed light on dinosaur physiology. Nature. 2012, 487 (7407): 358-361. 10.1038/nature11264.
Prondvai E, Stein KHW, De Ricqlès AJ, Cubo J: Development-based revision of bone tissue classification: the importance of semantics for science. Biol J Linn Soc. 2014, 112: 799-816. 10.1111/bij.12323.
Varricchio D, Sereno P, Zhao X, Tan L, Wilson J, Lyon G: Mud − trapped herd captures evidence of distinctive dinosaur sociality. Acta Palaeontol Pol. 2008, 53: 567-578. 10.4202/app.2008.0402.
Werning S: The ontogenetic osteohistology of Tenontosaurus tilletti . PLoS One. 2012, 7 (3): e33539-10.1371/journal.pone.0033539.
Erickson GM: Assessing dinosaur growth patterns: a microscopic revolution. Trends Ecol Evol. 2005, 20 (12): 677-684. 10.1016/j.tree.2005.08.012.
Woodward H, Padian K, Lee AH: Skeletochronology. Bone Histology of Fossil Tetrapods: Advancing Methods, Analysis, and Interpretation. Edited by: Padian K, Lamm ET. 2013, University of California Press, Berkeley, 187-207.
Lee AH, Huttenlocker AK, Padian K, Woodward H: Analysis of growth rates. Bone Histology of Fossil Tetrapods: Advancing Methods, Analysis, and Interpretation. Edited by: Padian K, Lamm ET. 2013, University of California Press, Berkeley, 217-251.
Horner J, de Ricqlès A, Padian K: Variation in skeletochronological indicators of the hadrosaurid dinosaur Hypacrosaurus: implications for age assessment of dinosaurs. Paleobiology. 1999, 25 (3): 49-78.
Cormack DH, Ham AW: Ham's histology. 1987, Lippincott, Philadelphia
Brusatte SL, Carr TD, Erickson GM, Bever GS, Norell MA: A long-snouted, multihorned tyrannosaurid from the Late Cretaceous of Mongolia. Proc Natl Acad Sci U S A. 2009, 106 (41): 17261-17266. 10.1073/pnas.0906911106.
Castanet J: Time recording in bone microstructures of endothermic animals; functional relationships. Compt Rendus Palevol. 2006, 5 (3–4): 629-636. 10.1016/j.crpv.2005.10.006.
Castanet J, Croci S, Aujard F, Perret M, Cubo J, de Margerie E: Lines of arrested growth in bone and age estimation in a small primate: microcebus murinus. J Zool. 2004, 263 (1): 31-39. 10.1017/S0952836904004844.
Castanet J, Francillon-Vieillot H, Meunier F, De Ricqles A: Bone and lndividual Aging. Bone: Bone Growth-B. 1993, 7: 245-
de Margerie E: Assessing a relationship between bone microstructure and growth rate: a fluorescent labelling study in the king penguin chick (Aptenodytes patagonicus). J Exp Biol. 2004, 207 (5): 869-879. 10.1242/jeb.00841.
Fowler DW, Woodward HN, Freedman EA, Larson PL, Horner JR: Reanalysis of Raptorex kriegsteini: A Juvenile Tyrannosaurid Dinosaur from Mongolia. PLoS One. 2011, 6 (6): e21376-10.1371/journal.pone.0021376.
Hubner TR: Bone histology in Dysalotosaurus lettowvorbecki (Ornithischia: Iguanodontia)–variation, growth, and implications. PLoS One. 2012, 7 (1): e29958-10.1371/journal.pone.0029958.
Sereno PC, Tan L, Brusatte SL, Kriegstein HJ, Zhao X, Cloward K: Tyrannosaurid skeletal design first evolved at small body size. Science. 2009, 326 (5951): 418-422. 10.1126/science.1177428.
Watanabe A, Erickson GM, Druckenmiller PS: An ornithomimosaurian from the upper cretaceous prince creek formation of Alaska. J Vertebr Paleontol. 2013, 33 (5): 1169-1175. 10.1080/02724634.2013.770750.
Carter Y, Thomas CD, Clement JG, Peele AG, Hannah K, Cooper DM: Variation in osteocyte lacunar morphology and density in the human femur–a synchrotron radiation micro-CT study. Bone. 2013, 52 (1): 126-132. 10.1016/j.bone.2012.09.010.
Skedros JG, Grunander TR, Hamrick MW: Spatial distribution of osteocyte lacunae in equine radii and third metacarpals: considerations for cellular communication, microdamage detection and metabolism. Cells Tissues Organs. 2005, 180 (4): 215-236. 10.1159/000088938.
Stein KWH, Werner J: Preliminary analysis of osteocyte lacunar density in long bones of tetrapods: all measures are bigger in sauropod dinosaurs. PLoS One. 2013, 8 (10): e77109-10.1371/journal.pone.0077109.
D'Emic MD, Benson RB: Measurement, variation, and scaling of osteocyte lacunae: a case study in birds. Bone. 2013, 57 (1): 300-310. 10.1016/j.bone.2013.08.010.
Montanari S, Brusatte SL, De Wolf W, Norell MA: Variation of osteocyte lacunae size within the tetrapod skeleton: implications for palaeogenomics. Biol Lett. 2011, 7 (5): 751-754. 10.1098/rsbl.2011.0173.
Cullen T, Ryan M, Currie P, Kobayashi Y, Schröder-Adams C: An Ornithomimid (Dinosauria) bonebed from the late cretaceous of Alberta, with implications for the behavior, classification, and stratigraphy of North American Ornithomimids. PLoS One. 2013, 8 (3): e58853-10.1371/journal.pone.0058853.
Brusatte SL, Carr TD, Norell MA: The Osteology of Alioramus, a gracile and long-snouted tyrannosaurid (Dinosauria: Theropoda) from the late cretaceous of Mongolia. Bull Am Mus Nat Hist. 2012, 366: 1-197. 10.1206/770.1.
Ibiricu LM, Martínez RD, Casal GA, Cerda IA: The behavioral implications of a multi-individual bonebed of a small Theropod Dinosaur. PLoS One. 2013, 8 (5): e64253-10.1371/journal.pone.0064253.
Woodward HN, Horner JR, Farlow JO: Osteohistological evidence for determinate growth in the American Alligator. J Herpetol. 2011, 45 (3): 339-342. 10.1670/10-274.1.
Evans DC, Barrett PM, Brink KS, Carrano MT: Osteology and bone microstructure of new, small theropod dinosaur material from the early Late Cretaceous of Morocco. Gondwana Research In Press.,
Campione NE, Evans DC, Brown CM, Carrano MT: Body mass estimation in non-avian bipeds using a theoretical conversion to quadruped stylopodial proportions. Methods in Ecology and Evolution In Press.,
Cooper LN, Lee AH, Taper ML, Horner JR: Relative growth rates of predator and prey dinosaurs reflect effects of predation. Proc Biol Sci. 2008, 275 (1651): 2609-2615. 10.1098/rspb.2008.0912.
Padian K, Horner JR, De Ricqlès A: Growth in small dinosaurs and pterosaurs: the evolution of archosaurian growth strategies. J Vertebr Paleontol. 2004, 24 (3): 555-571. 10.1671/0272-4634(2004)024[0555:GISDAP]2.0.CO;2.
Zelenitsky DK, Therrien F, Erickson GM, DeBuhr CL, Kobayashi Y, Eberth DA, Hadfield F: Feathered non-avian dinosaurs from North America provide insight into wing origins. Science. 2012, 338 (6106): 510-514. 10.1126/science.1225376.
Erickson GM, Tumanova TA: Growth curve of Psittacosaurus mongoliensis Osborn (Ceratopsia: Psittacosauridae) inferred from long bone histology. Zool J Linn Soc. 2000, 130 (4): 551-566. 10.1111/j.1096-3642.2000.tb02201.x.
Woodward HN, Horner JR, Farlow JO: Quantification of intraskeletal histovariability in Alligator mississippiensis and implications for vertebrate osteohistology. Peer J. 2014, 2: e422-10.7717/peerj.422.
Zhao Q, Benton MJ, Sullivan C, Martin Sander P, Xu X: Histology and postural change during the growth of the ceratopsian dinosaur Psittacosaurus lujiatunensis . Nat Commun. 2013, 4: 2079-
Bromage TG, Lacruz RS, Hogg R, Goldman HM, McFarlin SC, Warshaw J, Dirks W, Perez-Ochoa A, Smolyar I, Enlow DH, Boyde A: Lamellar bone is an incremental tissue reconciling enamel rhythms, body size, and organismal life history. Calcif Tissue Int. 2009, 84: 388-404. 10.1007/s00223-009-9221-2.
Stein K, Prondvai E: Rethinking the nature of fibrolamellar bone: an integrative biological revision of sauropod plexiform bone formation. Biol Rev. 2014, 89: 24-47. 10.1111/brv.12041.
Brown CM, Evans DC, Campione NE, O'Brien LJ, Eberth DA: Evidence for taphonomic size bias in the Dinosaur Park Formation (Campanian, Alberta), a model Mesozoic terrestrial alluvial, paralic system. Palaeogeogr Palaeoclimatol Palaeoecol. 2013, 372: 108-122. 10.1016/j.palaeo.2012.06.027.
Russell D: Ostrich dinosaurs from the late cretaceous of Western Canada. Can J Earth Sci. 1972, 9: 375-402. 10.1139/e72-031.
R Core Team: ‘Stats Package’ - R version 2.15.2. In: Vienna, Austria: R Foundation for Statistical Computing
Wickham H: ggplot2: elegant graphics for data analysis. In: New York: Springer; 2009.
We thank the staff of the Canadian Museum of Nature (K. Shepherd, C. Kennedy, M. Currie, A. Macdonald) for access to materials and permission to perform histological sampling, the staff of the Royal Ontario Museum (Kevin Seymour and Brian Iwama) for access to materials and use of histological thin sectioning equipment, J. Horner, D. Varricchio, S. Werning, D. Zelenitsky, and F. Therrien for prompt access to thin section images for comparison to sampled material, Koen Stein and one anonymous reviewer for their helpful reviews, and K. Brink, A. LeBlanc, K. Chiba, C. Brown, N. Campione, C. Van Buren, and D. Larson for useful discussions. K. Brink and A. LeBlanc are additionally thanked for providing helpful comments on the manuscript. We also thank D. Dufault and S. Hartman for permission to use skeletal and silhouette images. This study was supported by a Natural Sciences and Engineering Research Council (NSERC) of Canada Discovery Grant to D.C.E., and NSERC Alexander Graham Bell Canada Graduate Scholarship to T.M.C.
The authors declare that they have no competing interests.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Cullen, T.M., Evans, D.C., Ryan, M.J. et al. Osteohistological variation in growth marks and osteocyte lacunar density in a theropod dinosaur (Coelurosauria: Ornithomimidae). BMC Evol Biol 14, 231 (2014). https://doi.org/10.1186/s12862-014-0231-y