http://orcid.org/0000-0003-3561-0776
Figures
Abstract
The abundance of skeletal remains of cave bears in Pleistocene deposits can offer crucial information on the biology and life history of this megafaunal element. The histological study of 62 femora from 23 different European localities and comparisons with specimens of five extant ursid species revealed novel data on tissue types and growth patterns. Cave bear’s femoral bone microstructure is characterized by a fibrolamellar complex with increasing amounts of parallel-fibered and lamellar bone towards the outer cortex. Remodelling of the primary bone tissue initially occurs close to the perimedullary margin of the bone cortex around the linea aspera. Although similar histological traits can be observed in many extant bear species, the composition of the fibrolamellar complex can vary greatly. Cave bears reached skeletal maturity between the ages of 10 and 14, which is late compared to other bear species. There is a significant correlation between altitude and growth, which reflects the different body sizes of cave bears from different altitudes.
Citation: Veitschegger K, Kolb C, Amson E, Scheyer TM, Sánchez-Villagra MR (2018) Palaeohistology and life history evolution in cave bears, Ursus spelaeus sensu lato. PLoS ONE 13(11): e0206791. https://doi.org/10.1371/journal.pone.0206791
Editor: Carlo Meloro, Liverpool John Moores University, UNITED KINGDOM
Received: August 5, 2017; Accepted: October 19, 2018; Published: November 21, 2018
Copyright: © 2018 Veitschegger 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: All relevant data are within the paper and its Supporting Information files.
Funding: This research was mostly funded by the Swiss National Science Foundation (http://www.snf.ch) grant 31003A_149605 and 31003A_169395 to MRS-V. Additional funding for materials was provided by the Swiss National Science Foundation (http://www.snf.ch) grant 149506 to TMS. The Alexander von Humboldt Foundation (https://www.humboldt-foundation.de) and Deutsche Forschungsgemeinschaft (DFG AM 517/1-1) provided funding for the research of EA. 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
During the Pleistocene, Eurasia was home to many large sized mammals, commonly referred to as megafauna [1, 2]. Among them are cave bears, Ursus spelaeus, one of the most commonly found species, recovered in numerous Eurasian localities ranging from Spain to Russia [3, 4]. The closest living relatives of U. spelaeus are the brown bear, U. arctos and the polar bear, U. maritimus [5, 6]. The lineage of cave bears split between 2.75 and 1.4 Ma [5, 6] from the aforementioned two taxa. The hypothesized ancestral species of cave bears, U. deningeri is according to some authors, directly on the evolutionary line of U. spelaeus and, thus, the anagenetic ancestor of the latter [3, 7]. Cave bears occupied a wide range of habitats and remains of this species were found from about sea level up to 3000 m in altitude [8]. Molecular studies have shown that several distinct haplotype groups of cave bears are recognisable [9–12]. Some authors suggest giving species and subspecies status to some of those [8, 13] but it remains unclear if species status can be assigned to those cave bear haplotypes [4]. As this discussion is not completely resolved yet, we use U. spelaeus sensu lato, including U. ingressus, U. s. spelaeus, U. s. eremus, and U. s. ladinicus in the text but include analyses considering these haplotypes and possible species. High altitude specimens of U. spelaeus s.l. are often smaller [8]. High altitudes, in general, represent a challenging environment for mammals due to hypoxic conditions, which impair growth [14, 15].
Bone histology of extinct mammals has received increased interest in recent years, as details found in thin-sections provide insights into growth and thus life history [16–21]. These works focus on deposition of different bone tissue types and/or by analyzing growth cycles by measuring the distance between lines of arrested growth (LAGs) or growth annuli [22]. Cave bears went into dormancy during winters and therefore spent a considerable amount of time inactive during the year [23]. Non-hibernating mammals risk losing balance of bone resorption and formation in favour of resorption when being inactive for extended periods, which leads to extensive loss of bone minerals in a short time [24]. Studies on extant black bears, U. americanus, however, show that the balance between resorption and formation is not affected by extended inactive periods during winters but the cortical bone turnover rate is lower [25, 26]. The cortical bone geometry and strength of the femoral midshaft remain unchanged [26].
Much of our work is based on previous knowledge of bone formation and change, including Amprino’s rule, which postulates that different bone tissue types are produced in variable rates [27]. The fastest bone tissue to be produced is woven-fibered bone, followed by parallel-fibered and lamellar bone [28, 29]. Besides these three, a fibrolamellar complex is also found in mammals. During the formation of this complex initially highly vascularised woven-fibered scaffolding is deposited. In a later stage, the mostly reticular or plexiform vascular canals are filled with parallel-fibered and lamellar bone [29]. Recently, the term fibrolamellar complex has come under scrutiny as it does not provide detailed information about formation or structure of primary bone tissue [30]. Nonetheless, we continue using fibrolamellar complex to ensure comparability to previous studies [21]. A LAG is a thin opaque or translucent, circumferential line of compact bone and indicates a cessation of growth [22, 29]. In recent years, evidence accumulated supporting the yearly formation of LAGs and that this pattern is independent from climate or metabolic rate [18, 29, 31–34]. Here, we aim to study several aspects of the bone histology of cave bears: (1) we describe the general histology of the femora of cave bears compared to extant relatives, (2) we aim to reconstruct growth of cave bears by using average inter-LAG distance (AILD) as a proxy and compare it to extant relatives, (3) we present a hypothesis on some life history variables for cave bears based on the investigated growth, and (4) we examine the intraspecific variation in growth in the context of the altitude of the excavation site in order to elucidate if the investigated growth is different in the smaller cave bear populations from high altitudes.
Material and methods
No permits were required for the described study, which complied with all relevant regulations. All sampled specimens are part of museum collections. The providing institutions and curators/collection managers are named in the material and methods section and acknowledgments.
We investigated 62 femora of cave bears (U. spelaeus s.l.) from 23 different localities (Fig 1, S1 Table), as well as specimens of one other extinct (U. deningeri) and five extant species of Ursidae. U. deningeri is represented by three femora. From the 23 localities, twelve preserved the haplotype U. ingressus and seven U. spelaeus s.s. For the remaining four it remains unclear which haplotype was preserved at their respective caves (S2 Table).
The following institutions provided material for sampling from their collections: Aristotle University of Thessaloniki—School of Geology (AUTH); Biologiezentrum Linz (BZL); Canadian Museum of Nature (CMN); Finnish Museum of Natural History (MZH); Institut für Paläontologie, Universität Erlangen (IPUE); Institut für Paläontologie, Universität Wien (PIUW); Museum für Naturkunde Berlin (MfN); Muséum national d'Histoire naturelle, Paris (MNHN); Naturhistorisches Museum Basel (NMB); Naturhistorisches Museum Wien (NHM); Naturmuseum Südtirol Bozen (PZO); Naturmuseum St. Gallen (NMSG); Paläontologische Sammlung der Universität Tübingen (GPIT); and Swedish Museum of Natural History (SMNH). The thin-sections prepared and investigated for this study are deposited in the corresponding collections with the inventory number of their femora (S1 Table).
The sampled femora of extant bear species were those of Helarctos malayanus (one), Melursus ursinus (one), U. americanus (two), U. arctos (six), and U. maritimus (two), For U. spelaeus s.l. and U. arctos, ontogenetic series were available and in case of U. arctos the approximate age of three individuals was known (S1 Table). The preparation of the thin-sections followed standard procedures [17, 20, 35]. The midshaft of the femur was chosen, as it preserves the most complete growth record [17, 36]. We investigated bone growth in accordance to Kolb et al. [17] by measuring the distance between consecutive LAGs (S1 Table). Thus, more precisely, we investigated the amount of bone deposited in a growth cycle as proxy for growth in the respective cycle. To differentiate this form of growth studies [17, 20] from studies more concerned with average growth rate [37], we refered to our proxy as “average inter-LAG distance” (AILD) The histology of the bones was observed under normal transmitted or cross-polarized light using a Leica DM 2500 M composite microscope. In many cases, histological details were more pronounced using a lambda compensator. Pictures were taken using a Leica DFC 420 C digital camera. Growth zones were defined as the distance between two LAGs and measured using Leica IM 50 Image Manager software.
Measurements were acquired in the medial quadrant of the thin-sections of the femora because in this area LAGs were least obscured by remodelling (Figures A, B, I in S1 File). Pre-OCL inter-LAG distance was averaged (AILD) to compare it among species (S2 Table). The outer circumferential layer (OCL) denotes skeletal maturity [29]. Thus, we avoided bias resulting from mixing pre-maturity AILD with post-maturity AILD. For the comparison of AILD among different cave bear localities, growth zones 2 to 7 were averaged to include also samples lacking a full growth record (S2 Table). These averaged variables were also used as basis to compare the slopes and intercepts between U. spelaeus s.s. (sensu stricto) (incl. U. s. eremus and U. s. ladinicus) and U. ingressus. Comparative methods were employed to identify and eventually take into account the phylogenetic relationships of the sampled taxa. We investigated the species/haplotype averaged data with a PGLS (phylogenetic generalized least squares) as implemented in the R packages ape [38] and caper [39] using maximum likelihood (see also Losos [40]). We used a model defined by a variance-covariance matrix that reflects the phylogenetic signal observed in the data [41, 42]. It retrieves the phylogenetic signal index lambda (λ; [43]), of which a value of 0 reflects that data are independent from phylogeny (all nodes of the tree are collapsed resulting in a star phylogeny data structur), a value of 1 denoting that data distribution is consistent with a Brownian motion model of evolution (BM; no nodes are collapsed; [43]). We retrieved the phylogenetic information from Bon et al. [5], a tree based on mitochondrial genome and that included all sampled species. Knapp et al. [4] was used to place U. ingressus in the phylogeny. Branch lengths representing time of divergence among our investigated taxa were retrieved from Nyakatura and Bininda-Emonds [44], Bon et al. [5], and Knapp et al. [4]. The analyses excluded U. deningeri because of the position of this species in the cave bear phylogeny [3, 7]. AILD was mapped onto the built timetree (function contMap, package phytools), for which maximum likelihood reconstruction of the ancestral states is performed [45]. The phylogenetic framework is shown in Fig 2.
Log10-transformed AILD was mapped on the time tree using a gradient from blue (low values) to red (high values).
We investigated the possible relationship between AILD and the altitude of different localities. For this, we used OLS regressions and Kendall’s tau to statistically assess the correlation of the variables. In this study, the greatest latero-medial diameter (LMD) of the midshaft measured was used as proxy for body size for a species or locality and compared to AILD. Altitude was compared to AILD and midshaft LMD. To further study if the effect of altitude was size-dependent, we used size-corrected AILD by investigating the residuals of an OLS regression of the AILD against midshaft LMD and compared these residuals with altitude. AILD and midshaft LMD were log10-transformed. Statistical analyses were performed in R, version 3.2.3 [46] and with the additional package Kendall [47]. Plots were produced with ggplot2 [48] and Adobe Illustrator CS5.
Results
Skeletochronology
An overview of the histological features identified in this study is depicted on Fig 3.
Overview of histology of bears exemplified with Ursus spelaeus s.l. MB. Ma. 10881 (A) and Ursus arctos MNHN 1904–244 (B). The thin-sections of the medial part of the femoral midshaft are shown in cross-polarised light with lambda compensator. White arrows indicate LAGs, H = Haversian tissue, FLC = fibrolamellar complex, PF = parallel-fibered bone (mainly), SO = secondary osteon, WO = woven-fibered bone (mainly). (Scale bars: 2 mm).
Femora of neonate cave bears exhibit a fibrolamellar complex with high amounts of woven-fibered bone tissue. In the posterior part of the cortex, high amounts of parallel-fibered bone are deposited. The vascularization is mostly longitudinal as well as reticular and the parallel-fibered bone located posteriorly is less vascularised. With increasing individual age, the amounts of parallel-fibered and lamellar bone within a fibrolamellar complex increase while vascularization changes to predominantly laminar and plexiform arrangement. In the outer cortex of adult animals, an avascular OCL is present. A narrow lamellar endosteal layer and/or trabecular bone can be found in many individuals. In senile animals, a scarcely vascularised mixture of woven-fibered and parallel-fibered bone is present. The few vascular canals are usually longitudinally arranged (Fig 3). A varying amount of parallel-fibered and lamellar bone is the most frequent histological variation among U. spelaeus s.l. femora (Figures A, B, I in S1 File). Juvenile U. arctos exhibits a fibrolamellar complex dominated by woven-fibered bone with plexiform and laminar vascularisation. In the innermost part of the cortex, a thin endosteal layer comprised of avascular, laminar bone can be distinguished. Older animals exhibit a fibrolamellar complex with a higher amount of parallel-fibered bone in the outer parts of the cortex, where the vascularisation is laminar and plexiform. The overall change in histology from the inner to the outer part of the cortex is not very distinct. An endosteal layer and OCL are present (Fig 4, Figure F in S1 File).
White arrows indicate LAGs. Juvenile Ursus spelaeus s.l. PZO 5136 (Scale bar: 0.5 mm), subadult Ursus spelaeus s.l. MB.Ma.10881 (MfN) (Scale bar: 2 mm), adult Ursus spelaeus s.l. PEC 1183 (AUTH) (Scale bar: 2 mm), senile Ursus spelaeus s.l. LAC 6277b (AUTH) (Scale bar: 2 mm), juvenile Ursus arctos SMNH 2016–5132 (Scale bar: 1 mm), subadult Ursus arctos SMNH 2016–5131 (Scale bar: 2 mm), adult Ursus arctos SMNH 2016–5025 (Scale bar: 2 mm), senile Ursus arctos MNHN 1904–244 (Scale bar: 2 mm). Note that lamellar bone in the outer cortex of adult and senile specimens is indicated as OCL. H = Haversian tissue, FLC = fibrolamellar complex, L = lamellar bone, SO = secondary osteon, WO = woven-fibered bone (mainly).
In accordance with the already described histology of cave bears and brown bears, the femora of H. malayanus, U. americanus, U. deningeri, and U. maritimus exhibit a fibrolamellar complex with increasing amounts of parallel-fibered bone towards the outer surface of the cortex. The vascularization is primarily laminar and plexiform. In the outer cortex of the femora of H. malayanus, U. americanus, and U. maritimus, an OCL is present (Figures C, E, G, H in S1 File).
The bone of M. ursinus exhibits a very narrow medullary cavity, as also depicted by Mátyás [49]. On the whole, histology in M. ursinus varies among the anterior, medial, posterior, and lateral quadrants of the femoral cortex. The anterior part of the bone exhibits a matrix of woven-fibered and parallel-fibered bone with mostly longitudinal and reticular vascularization in the inner part of the cortex (fibrolamellar complex). The outer part of the cortex of the anterior quadrant of the bone shows a higher amount of parallel-fibered bone and a distinct avascular OCL. Medially, the inner cortex of the bone is comprised of woven-fibered bone with mostly longitudinal and reticular vascularization whereas the outer cortex exhibits parallel-fibered bone with low amounts of laminar vascularization and layers of avascular lamellar bone. Medially, the lamellar bone of the OCL is interlayered by woven-fibered bone with longitudinal and reticular vascularization (Fig 5). The lateral and posterior parts of the thin-section exhibit little primary bone tissue. Laterally, the outer cortex exhibits avascular lamellar bone interrupted by woven-fibered bone with laminar and reticular vascularization (Figure D in S1 File).
The thin-section is shown in cross-polarised light with lambda compensator. White arrows indicate LAGs, OCL = outer circumferential layer (mostly formed by lamellar bone), PF = parallel-fibered bone, WO = woven-fibered bone. Note the woven-fibered bone within the OCL. (Scale bars: 2 mm).
Bone remodelling in U. spelaeus s.l. starts in the posterior inner cortex with scattered secondary osteons. In the inner part of the cortex, especially posteriorly, dense Haversian bone is present already in young individuals. In older individuals, scattered secondary osteons are found in many parts of the cortex, but remodelling is especially prevalent anteriorly. Posteriorly, around the linea aspera—a ridge on the posterior section of femora—dense Haversian bone reaches the outer surface of the cortex (Fig 6). In their first year, the femora of U. arctos exhibit no remodelling of the bone. First signs of remodelling start in the posterior part of the bone around the linea aspera. Over time, dense Haversian bone tissue forms with several generations of secondary osteons, which reach the outer surface of the cortex. In the other regions of the bone, remodelling starts in the inner part of the cortex with scattered secondary osteons. With increasing age more scattered secondary osteons appear also in outer parts of the cortex. In old animals, a dense Haversian bone tissue forms in the inner cortex of the bone (Fig 4).
The thin-section is shown in cross-polarised light with lambda compensator. Note the dense Haversian tissue, comprised of several generations of secondary osteons, which extends to the outer cortex (Scale bar: 2 mm).
In M. ursinus, U. americanus, U. deningeri, and U. maritimus bone remodelling is also most pronounced around the linea aspera in the posterior part of the bone. In this part, Haversian bone is dense and reaches the outer part of the cortex. In other parts of the bone, especially in the inner part of the cortex, scattered secondary osteons are present. The anterior part of the femora often exhibits a higher amount of scattered secondary osteons.
Growth marks and growth in cave bears and closely related species
All examined bear species and cave bear populations with preserved histology exhibit LAGs except for a few early juveniles. In cave bears, an ontogenetic series of thin-sections revealed that the first LAG is resorbed around the time the sixth one is deposited (Fig 4). The ontogenetic series of U. arctos shows that the number of LAGs is consistently below the minimum age documented for the individual, evidencing the loss of at least one LAG during ontogeny. Thus, for comparison, the loss of the first LAG and the first growth zone was inferred for all adult bear specimens within the study. Growth zones of cave bears commonly exhibit a regular pattern. After the termination by a previous LAG, the growth zone starts with the deposition of woven-fibered bone. Over time, the amount of parallel-fibered bone increases and finally the growth zone terminates with the deposition of the subsequent LAG. In contrast, in U. arctos and U. americanus the growth zone is mostly a layer of a uniform fibrolamellar complex and only close to the cessation of growth, marked by a LAG, a thin layer of parallel-fibered bone is produced, also called annulus (Fig 7). The growth zones of H. malayanus, M. ursinus U. deningeri, and U. maritimus are less uniform and exhibit somewhat alternating parallel-fibered and woven-fibered bone in the fibrolamellar complex (Figure J in S1 File).
White arrows indicate LAGs. Note that inner part of the growth zone exhibits more WO bone whereas the external part exhibits more PF bone. WO = woven-fibered bone, PF = parallel-fibered bone. (Scale bars: 0.5 mm).
The timing of the cessation of circumferential bone growth is indicated by the OCL [29]. Usually, an OCL is a layer of avascular lamellar bone potentially comprising closely spaced LAGs. The OCL of cave bears appears between LAG 9 and 13 and in U. arctos it appears between LAG 6 and 7. Skeletal maturity was reached in U. americanus after LAG 7, in H. malayanus after LAG 5, in U. maritimus after LAG 10, and in M. ursinus after LAG 6. No OCL was recorded in the available specimens of U. deningeri, suggesting that this animal did not reach skeletal maturity before the age of eleven (ten LAGs plus one resorbed, S1 Table).
The overall growth, as presented in Fig 2 and Fig 8A (using AILD as a proxy), shows that cave bears were growing in comparable speed to U. maritimus. The second best sampled species U. arctos exhibits slower growth rates. In our sampling, the fast growing species were U. deningeri, U. maritimus, and U. spelaeus s.l. (U. ingressus, U spelaeus s.s., the slow growing species were U. americanus, U. arctos, and M. ursinus (Fig 2, Fig 8A). For its size, according to our data, the fastest growing species was H. malayanus. However, the placement of this species in the scatterplot (Fig 8A) is unique. The phylogenetically informed analysis showed that the relationship among bear species is not significantly correlated with the scaling of their AILD (λ = 0) (Table 1, the coefficients of the OLS regression are very similar). Accordingly, mapping of AILD onto the ursid timetree (Fig 2) does not reveal clade-specific patterns.
Scatter plot of: a) AILD of different bear species compared to their femoral latero-medial diameterat midshaft (LMD) as a proxy for body size (open circle not included in phylogenetic analyses), b) AILD of different cave bear localities compared to their LMD, c) AILD of different cave bear localities compared to the altitude of the locality, d) AILD/LMD-residuals compared to the altitude of the locality (results of OLS regression in the left box, results of Kendall’s tau in the right box).
The standardized growth among different cave bear localities showed a considerable amount of variation in AILD (Fig 8B). Within cave bears, growth was significantly correlated with size (Fig 7B, Table 1). The slowest growth was recorded for Arzberghöhle, Zoolithenhöhle, and Brettsteinbärenhöhle whereas the fastest growth was documented for Loutraki and Gondenans-les-Moulins (Fig 8B). Kendall’s tau shows that AILD was correlated with altitude of the locality (Fig 8C, Table 1). However, this correlation is due to the effect of size as shown by the absence of a correlation between AILD/midshaft LMD-residuals and altitude (Fig 8D, Table 1). Indeed, in our sample, the proxy for body size (midshaft LMD) was correlated with altitude (Table 1), thus, smaller bears are found at higher altitudes.
The comparison of the slopes and intercepts of the cave bear haplotypes U. spelaeus s.s. and U. ingressus (Table 2) revealed a significant difference in the slope of these two groups (as indicated by the significant interaction between midshaft LMD and haplotype; Table 3). The slope of U. spelaeus s.s. is slightly negative (Table 2). However, the correlation is not significant. In contrast, U. ingressus has a significant positive correlation between growth rate and latero-medial diameter of the femora, suggesting that bigger animals grew faster than smaller ones. The intercept of both groups is not significantly different (Table 3) evidencing that there is no difference in AILD between these two haplotypes.
Femoral latero-medial diameter at midshaft (LMD) was used as a proxy for body size.
Comparison based on OLS regressions for the two different cave bear haplotypes U.spelaeus (incl. U. s. eremus and U. s. ladinicus) and U. ingressus with femoral latero-medial diameter at midshaft (LMD) as a covariate, including or excluding its interaction with the category of haplotype (upper panel and lower panel, respectively).
Discussion
Using AILD as a proxy for growth, our data suggest that extant and extinct bear species grew rather uniformly, with the exception of H. malayanus. It remains unclear, if the growth documented in our sample of H. malayanus is representative for the species or not. The origin of the bone is unknown; thus, it might be possible that the animal was raised under protected conditions in a zoo. The histology itself does not show clear differences compared to the Ursus species, despite scattered secondary osteons. However, the amount of secondary osteons was very variable in bear species with higher sample sizes.
The observed intraspecific variation within U. spelaeus s.l. was high (Fig 8A and 8B, Table 1). The histology of U. americanus and U. arctos, however, suggests a fast growth because of the amounts of woven-fibered bone. On the other hand, U. maritimus exhibits a great amount of slowly deposited parallel-fibered and lamellar bone, which would suggest a slower growth than found based on AILD. However, in contrast to artiodactyls [18], growth periods during the year are not known for bear species. As aforementioned, bears slow down bone deposition and resorption during hibernation [26]. Thus the extent of the cessation of growth between species could be different. Additionally, the length of hibernation varies among populations of one species [50]. The observed variation in AILD among cave bears from different localities was best explained by their size, which in turn is correlated to the altitude of cave entrance. Smaller high alpine cave bear populations [13] have AILDs that are not significantly smaller than the ones from their lowland relatives. In contrast, experiments mimicking high altitude conditions have shown that growth rate is usually smaller in elevated environments due to hypoxia as well as reduced appetite [14, 51, 52]. This effect is positively correlated with altitude [14, 53]. However, findings by Elia et al. [14] suggest that this effect is not found under 1850 m altitude. Our data show that size dependent AILD drops already around 500 m above sea level, excluding the outlying caves Arzberghöhle and Zoolithenhöhle (Fig 8C).
The presence of a fibrolamellar complex in bear species follows previously described thin-sections for the genus Ursus [54, 55]. However, notable differences were recorded in the composition of the fibrolamellar complex. Whereas, U. spelaeus s.l., U. arctos, and U americanus exhibit a bone histological pattern similar to the one described by Enlow and Brown [54], the one in other bear species is characterized by higher amounts of parallel-fibered and lamellar bone. A notable difference was found in M. ursinus, with considerable amounts of longitudinal vascular canals only present in this species. One of the first histological examinations of U. spelaeus s.l., from the Drachenhöhle of Mixnitz, described lamellar bone structure within trabecular bone [56]. Later studies on cave bear histology focused on mineralogical examinations of remains to investigate taphonomic processes [57] or pathological changes to the bone [58].
The ontogenetic change of bone tissue in U. arctos and U. spelaeus s.l. mirrors findings in deer femora with the exception that remodelling in bears starts in the inner cortex [17]. The cessation of bone growth as marked by an OCL has been associated with either sexual maturity [59, 60] or skeletal maturity [16]. In U. arctos, femora are skeletally mature between 5 and 8 years [61] and animals reach sexual maturity between the age of 3.5 and 5.5 years [50]. The OCL in U. arctos is produced between 7 and 8 years (6–7 LAGs plus one resorbed) suggesting that it indicates skeletal maturity. For the other examined bear species, the appearance of the OCL cannot be linked to sexual maturity because the amount of observed LAGs usually exceeds the known age of sexual maturity [50]. In cave bears the OCL appears between LAG 9 and 13, which suggest attaining femoral skeletal maturity at an advanced age. Evidence from histological thin-sections of Cervidae also suggest that the OCL is linked to skeletal maturity [17]. The best sampled species, U. spelaeus s.l. exhibits a gradual decrease in growth rate over time. In contrast, in deer the growth rate drops after growth zone 3 and remains fairly constant afterwards [17].
The remodelling of the femora of all examined bear species generally follows a very similar spatial pattern. Haversian tissue is dense around the linea aspera and the inner part of the cortex, whereas in the rest of the bone more scattered secondary osteons are present. A similar pattern of remodelling was described by Mátyás [49]. In deer, remodelling is also strongest close to the linea aspera [17]. The linea aspera serves as point of attachments for several muscles [62]. Thus, the remodelling might be the result of increased mechanical stress in this area.
As for now, the effect of hibernation on the composition of the fibrolamellar-complex remains elusive. Hibernation is a possible explanation for the formation of zonal bone in dinosaurs [55] but mammals express zonal bone independent of hibernation [this study, 21]. Our sampling does not conclusively add further information. Although bone histology of the non-hibernating bear M. ursinus is different from the other sampled bears, the same is not true for H. malayanus, which does not hibernate either [50].
Conclusion
Cave bear offspring had about the same size as the ones of U. arctos at birth [63–65] and thus grew faster [65] and longer (as indicated by the OCL) to attain their adult stature of 500 kg [66] or probably even up to 1,500 kg [3]. The maximum known age for cave bears is 30–32 years [67, 68]. However, most individuals did not reach that age, as mortality was often high among juveniles [67–70].
The investigated haplotypes of cave bears exhibit different growth patterns indicating variation in life history strategy during early phases of speciation. There is evidence that cave bears had a faster life history compared to close relatives because of their smaller relative brain size [71]. However, the overall pace of growth reconstructed for cave bears in this study is similar to those of their close relatives.
Supporting information
S1 Table. List of investigated specimens (with their ontogenetic stage and locality) and related histomorphometrical measurements.
https://doi.org/10.1371/journal.pone.0206791.s001
(XLSX)
S2 Table. Averaged inter-LAG distances by species or locality for U. spelaeus s.l.
https://doi.org/10.1371/journal.pone.0206791.s002
(XLSX)
S1 File. Additional figures depicting the observed variation in histology among different cave bear localities and the different thin-section quadrants of H. malayanus, M. ursinus, U. americanus, U. arctos, U. deningeri, and U. maritimus.
https://doi.org/10.1371/journal.pone.0206791.s003
(PDF)
Acknowledgments
We thank the following colleagues for kindly providing us with samples for this study: Christine Argot (MNHN), Björn Berning (BZL), Toni Bürgin (NMSG), Loïc Costeur (NMB), Stephanie Fassl (PIUW), Christine Frischauf (PIUW), Ursula B. Göhlich (NHM), Oliver Hampe (MfN), Voitto Haukisalmi (MZH), Brigitte Hilpert (IPUE), Daniela Kalthoff (SMNH), Kamal Khidas (CMN), Evelyn Kustatscher (PZO), Gernot Rabeder (PIUW), Thomas Schossleitner (MfN), Laura Smyk (CMN), Martin Studeny (BZL), Evangelia Tsoukala (AUTH), Géraldine Veron (MNHN), and Ingmar Werneburg (GPIT). Additionally, we want to thank Alexandra Houssaye, Gernot Rabeder, and other anonymous reviewers for their helpful suggestions. John Finarelli (Dublin) and Lukas Keller (Zurich) provided comments to earlier versions of this manuscript and the taxonomy for the multi-species concept of cave bear evolution was provided by Gernot Rabeder. This research was funded by the Swiss National Science Foundation (SNSF) grant 31003A_169395 to Marcelo R. Sánchez-Villagra and partly supported by funds from the SNF grant 149506 to Torsten M. Scheyer. Eli Amson was funded by the Alexander von Humboldt Foundation and the German Research Council (DFG AM 517/1-1).
References
- 1. Barnosky AD, Koch PL, Feranec RS, Wing SL, Shabel AB (2004) Assessing the causes of Late Pleistocene extinctions on the continents. Science 306(5693):70–75. pmid:15459379
- 2. Lorenzen ED, Nogues-Bravo D, Orlando L, Weinstock J, Binladen J, Marske KA, et al. (2011) Species-specific responses of Late Quaternary megafauna to climate and humans. Nature 479(7373):359–364. pmid:22048313
- 3.
Rabeder G, Nagel D, Pacher M (2000) Der Höhlenbär. Stuttgart: Thorbecke Verlag.
- 4. Knapp M, Rohland N, Weinstock J, Baryshnikov G, Sher A, Nagel D, et al. (2009) First DNA sequences from Asian cave bear fossils reveal deep divergences and complex phylogeographic patterns. Mol Ecol 18(6):1225–1238. pmid:19226321
- 5. Bon C, Caudy N, de Dieuleveult M, Fosse P, Philippe M, Maksud F, et al. (2008) Deciphering the complete mitochondrial genome and phylogeny of the extinct cave bear in the Paleolithic painted cave of Chauvet. PNAS 105(45):17447–17452. pmid:18955696
- 6. Krause J, Unger T, Noçon A, Malaspinas A-S, Kolokotronis S-O, Stiller M, et al. (2008) Mitochondrial genomes reveal an explosive radiation of extinct and extant bears near the Miocene-Pliocene boundary. BMC Evol Biol 8(1):1–12.
- 7. Rabeder G, Hofreiter M (2004) Der neue Stammbaum der alpinen Höhlenbären. Die Höhle 5558–77.
- 8. Rabeder G, Debeljak I, Hofreiter M, Withalm G (2008) Morphological responses of cave bears (Ursus spelaeus group) to high-alpine habitats. Die Höhle 5959–72.
- 9. Dabney J, Knapp M, Glocke I, Gansauge M-T, Weihmann A, Nickel B, et al. (2013) Complete mitochondrial genome sequence of a Middle Pleistocene cave bear reconstructed from ultrashort DNA fragments. PNAS 110(39):15758–15763. pmid:24019490
- 10. Stiller M, Molak M, Prost S, Rabeder G, Baryshnikov G, Rosendahl W, et al. (2014) Mitochondrial DNA diversity and evolution of the Pleistocene cave bear complex. Quat Int 339-340224-231.
- 11. Baca M, Mackiewicz P, Stankovic A, Popović D, Stefaniak K, Czarnogórska K, et al. (2014) Ancient DNA and dating of cave bear remains from Niedźwiedzia Cave suggest early appearance of Ursus ingressus in Sudetes. Quat Int 339-340217-223.
- 12. Knapp M (2014) From a molecules’ perspective–contributions of ancient DNA research to understanding cave bear biology. Hist Biol Online First ().
- 13. Rabeder G, Hofreiter M, Nagel D, Withalm G (2004) New taxa of Alpine cave bears (Ursidae, Carnivora). Cah Sci Hors série 249–67.
- 14. Elia R, Elgoyhen A, Bugallo G, Río M, Bozzini C (1985) Effect of acute exposure to reduced atmospheric pressures on body weight, food intake and body composition of growing rats. Acta Physiol Pharmacol Latinoam 35(3):311–318. pmid:2938414
- 15. Delaquerrière-Richardson L, Forbes S, Valdivia E (1965) Effects of stimulated high altitude on the growth rate of albino guinea pigs. J Appl Physiol 20(5):1022–1025. pmid:5891210
- 16. Nacarino-Meneses C, Jordana X, Köhler M (2016) First approach to bone histology and skeletochronology of Equus hemionus. C R Palevol 15(1–2):267–277.
- 17. Kolb C, Scheyer TM, Lister AM, Azorit C, de Vos J, Schlingemann MA, et al. (2015) Growth in fossil and extant deer and implications for body size and life history evolution. BMC Evol Biol 15(1):19.
- 18. Köhler M, Marín-Moratalla N, Jordana X, Aanes R (2012) Seasonal bone growth and physiology in endotherms shed light on dinosaur physiology. Nature 487358–361.
- 19. Köhler M, Moyà-Solà S (2009) Physiological and life history strategies of a fossil large mammal in a resource-limited environment. PNAS 106(48):20354–20358. pmid:19918076
- 20. Amson E, Kolb C, Scheyer TM, Sánchez-Villagra MR (2015) Growth and life history of Middle Miocene deer (Mammalia, Cervidae) based on bone histology. C R Palevol 14(8):637–645.
- 21. Kolb C, Scheyer TM, Veitschegger K, Forasiepi AM, Amson E, Van der Geer AAE, et al. (2015) Mammalian bone palaeohistology: a survey and new data with emphasis on island forms. PeerJ 3e1358.
- 22.
Castanet J, Francillon-Vieillot H, Meunier FJ, de Ricqlès A. Bone and individual aging. In: Hall BK, editor. Bone Volume 7: Bone growth—B. Boca Raton: CRC Press; 1993. p. 245–283.
- 23. Fortes GG, Grandal-d'Anglade A, Kolbe B, Fernandes D, Meleg IN, García-Vázquez A, et al. (2016) Ancient DNA reveals differences in behaviour and sociality between brown bears and extinct cave bears. Mol Ecol 25(19):4907–4918. pmid:27506329
- 24. Giangregorio L, Blimkie CJR (2002) Skeletal adaptations to alterations in weight-bearing activity. Sports Med 32(7):459–476. pmid:12015807
- 25. McGee-Lawrence M, Buckendahl P, Carpenter C, Henriksen K, Vaughan M, Donahue S (2015) Suppressed bone remodeling in black bears conserves energy and bone mass during hibernation. J Exp Biol 218(13):2067–2074.
- 26. McGee-Lawrence ME, Carey HV, Donahue SW (2008) Mammalian hibernation as a model of disuse osteoporosis: the effects of physical inactivity on bone metabolism, structure, and strength. Am J Physiol Regul Integr Comp Physiol 295(6):R1999–R2014. pmid:18843088
- 27. Amprino R (1947) La structure du tissu osseux envisagee comme expression de differences dans la vitesse de l'accroissement. Arch Biol 58315–330.
- 28. de Margerie E, Robin JP, Verrier D, Cubo J, Groscolas R, Castanet J (2004) Assessing a relationship between bone microstructure and growth rate: a fluorescent labelling study in the king penguin chick (Aptenodytes patagonicus). J Exp Biol 207(5):869–879.
- 29.
Huttenlocker AK, Woodward HN, Hall BK. The biology of bone. In: Padian K, Lamm E-T, editors. Bone histology of fossil tetrapods—advancing methods, analysis, and interpretation. Berkeley and Los Angeles: University of California Press; 2013. p. 13–34.
- 30. Stein K, Prondvai E (2014) Rethinking the nature of fibrolamellar bone: an integrative biological revision of sauropod plexiform bone formation. Biol Rev 89(1):24–47. pmid:23647662
- 31.
Hinrichs J (2016). Are lines of arrested growth in bone indicative of seasonal metabolic supression in bears. Fort Collins, Colorado: Colorado State University.
- 32. Marangoni F, Schaefer E, Cajade R, Tejedo M (2009) Growth-mark formation and chronology of two neotropical anuran species. J Herpetol 43(3):546–550.
- 33. Buffrénil V, Castanet J (2000) Age estimation by skeletochronology in the nile monitor (Varanus niloticus), a highly exploited species. J Herpetol 34(3):414–424.
- 34. Castanet J, Croci S, Aujard F, Perret M, Cubo J, Margerie E (2004) Lines of arrested growth in bone and age estimation in a small primate: Microcebus murinus. J Zool 263
- 35. Chinsamy A, Raath MA (1992) Preparation of fossil bone for histological examination. Palaeontol afr 2939–44.
- 36. Sander PM, Andrassy P (2006) Lines of arrested growth and long bone histology in Pleistocene large mammals from Germany: what do they tell us about dinosaur physiology? Palaeontogr Abt A 277143–159.
- 37. Cubo J, Le Roy N, Martinez-Maza C, Montes L (2012) Paleohistological estimation of bone growth rate in extinct archosaurs. Paleobiology 38(2):335–349.
- 38. Paradis E, Claude J, Strimmer K (2004) APE: analyses of phylogenetics and evolution in R language. Bioinformatics 20:289–90. pmid:14734327
- 39. Orme D, Freckleton R, Thomas G, Petzoldt T, Fritz S, Isaac N, et al. (2013) caper: comparative analyses of phylogenetics and evolution in R. R package version 0.5.2. https://CRAN.R-project.org/package=caper. Accessed 25 Jan 2016.
- 40. Losos JB (2011) Seeing the forest for the trees: the limitations of phylogenies in comparative biology. Am Nat 177(6):709–727. pmid:21597249
- 41. Revell LJ (2010) Phylogenetic signal and linear regression on species data. Methods Ecol Evol 1(4):319–329.
- 42. Clauss M, Dittmann MT, Müller DWH, Meloro C, Codron D (2013) Bergmann′ s rule in mammals: a cross-species interspecific pattern. Oikos 122(10):1465–1472.
- 43. Pagel M. (1999) Inferring the historical patterns of biological evolution. Nature, 401,877–884. pmid:10553904
- 44. Nyakatura K, Bininda-Emonds OR (2012). Updating the evolutionary history of Carnivora (Mammalia): a new species-level supertree complete with divergence time estimates. BMC Biol 10(1):1–31.
- 45. Revell LJ (2012) phytools: An R package for phylogenetic comparative biology (and other things). Methods Ecol Evol 3 217–223.
- 46. R Development Core Team (2015) R: A language and environment for statistical computing. Vienna, Austria Version 3.2.3
- 47. McLeod AI (2011) Kendall: Kendall rank correlation and Mann-Kendall trend test. R package Version 2.2
- 48.
Wickham H (2009) ggplot2: Elegant graphics for data analysis. New York: Springer-Verlag.
- 49.
Mátyás J (1955) Mikroskopische Untersuchungen der biologischen Resorptionen in den Röhrenknochen. Budapest: Akadémiai Kiadó.
- 50.
Hunter L (2011) Carnivores of the world. Princeton (New Jersey): Princeton University Press.
- 51. Westerterp-Plantenga MS, Westerterp KR, Rubbens M, Verwegen CRT, Richelet J-P, Gardette B (1999) Appetite at “high altitude” [Operation Everest III (Comex-’97)]: a simulated ascent of Mount Everest. J Appl Physiol 87(1):391–399. pmid:10409600
- 52. Morel OE, Aubert R, Richalet J-P, Chapelot D (2005) Simulated high altitude selectively decreases protein intake and lean mass gain in rats. Physiol Behav 86(1):145–153.
- 53. Bozzini C, Champin GM, Alippi RM, Bozzini CE (2013) Static biomechanics in bone from growing rats exposed chronically to simulated high altitudes. High Alt Med Biol 14(4):367–374. pmid:24377344
- 54. Enlow DH, Brown SO (1958) A comparative histological study of fossil and recent bone tissues. Part III. Tex J Sci 10187–230
- 55. Chinsamy A, Rich T, Vickers-Rich P (1998) Polar dinosaur bone histology. J Vert Paleontol 18(2):385–390.
- 56.
Breuer R. Zur Anatomie, Pathologie und Histologie der Zähne und der Kiefer von Ursus spelaeus. In: Abel O, Kyrle G, editors. Die Drachenhöhle bei Mixnitz. Speläologische Monographien 7/8. Wien: Österreichische Staatsdruckerei; 1931. p. 581–623.
- 57. Rogóż A, Sawłowicz Z, Socha P, Stefaniak K (2009) Mineralization of teeth and bones of the cave bear (Ursus spelaeus) from the Biśnik Cave, Southern Poland. Mineralogia 40(1–4):65–84.
- 58. Nowakowski D, Stefaniak K (2015) Pathological changes of the cranium of a young female cave bear (Ursus spelaeus R.)-a case study (the Sudety Mts, Poland). Int J Osteoarchaeol 25(1):119–125.
- 59. Jordana X, Marín-Moratalla N, Moncunill-Solè B, Nacarino-Meneses C, Köhler M (2016) Ontogenetic changes in the histological features of zonal bone tissue of ruminants: a quantitative approach. C R Palevol 15(1–2):255–266.
- 60. Marin-Moratalla N, Jordana X, Köhler M (2013) Bone histology as an approach to providing data on certain key life history traits in mammals: implications for conservation biology. Mamm Biol 78
- 61. Weinstock J (2009) Epiphyseal fusion in brown bears: a population study of grizzlies (Ursus arctos horribilis) from Montana and Wyoming. Int J Osteoarchaeol 19(3):416–423.
- 62.
Drake RL, Wayne Vogl A, Mitchell AWM (2015) Gray's anatomy for students, 3rd edition. Philadelphia: Churchill Livingstone/Elsevier.
- 63.
Ehrenberg K. Die Variabilität der Backenzähne beim Höhlenbären. In: Abel O, Kyrle G, editors. Die Drachenhöhle bei Mixnitz. Speläologische Monographien 7/8. Wien: Österreichische Staatsdruckerei; 1931. p. 537–573.
- 64. Ehrenberg K (1973) Ein fast vollständiges Höhlenbärenneonatenskelett aus der Salzofenhöhle im Toten Gebirge. Ann Nathist Mus Wien 7769–113.
- 65. Fosse P, Cregut-Bonnoure E (2014) Ontogeny/growth of (sub)modern brown bear (Ursus arctos) skeleton: A guideline to appraise seasonality for cave bear (Ursus spelaeus) sites? Quat Int 339275–288.
- 66. Christiansen P (1999) What size were Arctodus simus and Ursus spelaeus (Carnivora: Ursidae)? Ann Zool Fenn 36(2):93–102.
- 67. Debeljak I (2007) Fossil population structure and mortality of the cave bear from the Mokrica cave (North Slovenia). Acta Carsologica 36(3):475–484.
- 68. Veitschegger K, Kolb C, Amson E, Sánchez-Villagra MR (2018) Longevity and life history of cave bears–a review and novel data from tooth cementum and relative emergence of permanent dentition. Hist Biol Online First ().
- 69. Grandal-D'Anglade A, Vidal Romaní JR (1997) A population study on the cave bear (Ursus spelaeus Ros.-Hein.)) from Cova Eirós (Triacastela, Galicia, Spain). Geobios 30(5):723–731.
- 70.
Debeljak I. Determination of individual age and season at death in cave bear from Ajdovska jama near Krško (Slovenia). In: Pacher M, Pohar V, Rabeder G, editors. Ajdovska Jama—Palaeontology, zoology and archaeology of Ajdovska Jama near Krško in Slovenia. Mitteilungen der Kommission für Quartärforschung der östereichischen Akademie der Wissenschaften, Volume 20. Wien: Österreichische Akademie der Wissenschaften; 2011. p. 51–63.
- 71. Veitschegger K (2017) The effect of body size evolution and ecology on encephalization in cave bears and extant relatives. BMC Evol Biol 17(1):124. pmid:28583080