Next Article in Journal
Pet Ownership in Aotearoa New Zealand: A National Survey of Cat and Dog Owner Practices
Next Article in Special Issue
Naked Mole-Rats Demonstrate Profound Tolerance to Low Oxygen, High Carbon Dioxide, and Chemical Pain
Previous Article in Journal
Are Currently Selected Laboratory Animals Useful in the Research of How Female Hormones Influence Orthodontic Biomechanics?
Previous Article in Special Issue
Laterality in the Damaraland Mole-Rat: Insights from a Eusocial Mammal
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pup Recruitment in a Eusocial Mammal—Which Factors Influence Early Pup Survival in Naked Mole-Rats?

by
Michaela Wetzel
1,*,
Alexandre Courtiol
1,
Heribert Hofer
1,2,
Susanne Holtze
1 and
Thomas B. Hildebrandt
1,2
1
Leibniz Institute for Zoo and Wildlife Research, Alfred-Kowalke-Str. 17, 10315 Berlin, Germany
2
Department of Veterinary Medicine, Freie Universität Berlin, Oertzenweg 19b, 14163 Berlin, Germany
*
Author to whom correspondence should be addressed.
Animals 2023, 13(4), 630; https://doi.org/10.3390/ani13040630
Submission received: 22 January 2023 / Revised: 8 February 2023 / Accepted: 9 February 2023 / Published: 11 February 2023
(This article belongs to the Special Issue Physiology and Behaviour of African Mole-Rats)

Abstract

:

Simple Summary

The naked mole-rat is a small, long-lived mammal that lives in large colonies in subterranean burrows. The lifestyle of naked mole-rats makes their study difficult, and many details about their biology remain enigmatic. For evolutionary biologists and ecologists, a special trait exhibited by this species, their eusocial lifestyle, is of particular interest. In this study, we explored the determinants of early offspring survival using data from observations of 14 captive colonies over a total period of seven years. Our study revealed that early pup survival was significantly improved by higher pup body mass and maternal number of mammae and significantly reduced by increased maternal body mass and colony size. The latter negative effect may, however, be an artifact of the captive conditions in which colonies were kept. We further discuss the implications of the similarities and differences of naked mole-rats to eusocial insects and cooperatively breeding mammals in these determinants to shed light on the origin and maintenance of eusociality in mammals.

Abstract

In eusocial insects, offspring survival strongly depends on the quality and quantity of non-breeders. In contrast, the influence of social factors on offspring survival is more variable in cooperatively breeding mammals since maternal traits also play an important role. This difference between cooperative insects and mammals is generally attributed to the difference in the level of sociality. Examining offspring survival in eusocial mammals should, therefore, clarify to what extent social organization and taxonomic differences determine the relative contribution of non-breeders and maternal effects to offspring survival. Here, we present the first in-depth and long-term study on the influence of individual, maternal, social and environmental characteristics on early offspring survival in a eusocial breeding mammal, the naked mole-rat (Heterocephalus glaber). Similarly to other mammals, pup birth mass and maternal characteristics such as body mass and the number of mammae significantly affected early pup survival. In this eusocial species, the number of non-breeders had a significant influence on early pup survival, but this influence was negative—potentially an artifact of captivity. By contrasting our findings with known determinants of survival in eusocial insects we contribute to a better understanding of the origin and maintenance of eusociality in mammals.

1. Introduction

Many factors can influence offspring survival in mammals (Table 1), including individual [1,2], maternal [3,4], social [5,6,7] and environmental factors [1,2,3,4,5,6,7,8]. In particular, maternal characteristics are known to have a strong impact on early juvenile survival (Table 1) [1,3,9,10] because of the prolonged period of offspring nutritional dependence during gestation and lactation [9]. In mammals, offspring survival is, therefore, likely to depend on a breeding female’s ability to lactate and deliver adequate nutrients and energy through milk [11,12]. Whether this situation also applies to eusocial mammals is, however, an open question. On the one hand, studies from non-eusocial cooperatively breeding mammals show that in many cases, reproduction is monopolized by a dominant breeding pair [13], and brood care is provided by both parents and subordinate group members [14,15]. On the other hand, studies from eusocial insects show that in these species, parental investment is almost exclusively limited to brood production (egg-laying) by the queen(s), whereas all brood care is provided by helpers termed workers [16].
The African mole-rat family Bathyergidae is unique among mammals with respect to their variable social organizations, ranging from solitary to eusocial lifestyles [17,18]. The Damaraland mole-rat (Fukomys damarensis) and some of its social relatives in the genera Cryptomys and Fukomys were classified as eusocial by some researchers, while others disagree [16,18,19,20]. In contrast, the naked mole-rat (Heterocephalus glaber), with the highest reproductive skew of all social mole-rats, is the species with a lifestyle that most resembles that of eusocial insects and, therefore, unambiguously meets the original definition of eusociality [21,22,23]. Similar to eusocial diploid insects such as termites (Isoptera), naked mole-rat colonies exhibit cooperative brood care, a reproductive division of labor with a single, dominant breeding female (the queen), overlapping generations of several subordinate and reproductively suppressed colony members (non-breeders) [19] and distinct morphological and behavioral castes [24,25,26]. Naked mole-rats are, therefore, an appropriate species to study maternal as well as non-breeder effects on early pup survival and determine whether eusocial mammals have a distribution of parental care that resembles either that of other cooperatively breeding mammals or that of other truly eusocial species.
Although factors influencing offspring survival have already been investigated in most cooperatively breeding mammals [15], they have not yet been studied in detail in any eusocial mammal. Braude [27] showed that wild colonies of naked mole-rats have high rates of both recruitment and loss. Wild colonies increased by as many as 104 pups per year (mean ± standard error = 34 ± 25, n = 42 cases in which a colony was captured in successive years) and lost between 12% and 91% of their non-breeders per year (46 ± 23). As pups younger than two weeks do not leave their subterranean nest [28], early postnatal survival rates of wild pups remain unexplored. In captivity, pup survival is known to be highly variable and differs during various life stages [28,29]. As in other mammals [4,30,31,32], naked mole-rat pup survival is reported to be lowest during the first postnatal days [33,34]. Moreover, pup survival in naked mole-rats may be influenced by maternal and social factors [28,35,36]. For instance, Jarvis [28] showed that neonates may die within the first few postnatal days caused by maternal neglect and problems with milk production. Other colony members in addition to the queen may also influence survival of littermates. Non-breeders provide both direct and indirect pup care, e.g., huddling in the nest, grooming, donating caecotrophs, carrying food to the nest, maintaining and defending the burrow system [35]. They also were seen performing detrimental behaviors such as pushing, displacing or consuming pups [28,35,37]. Payne [37] observed that young were sometimes ignored or actually abused by non-breeders. Jarvis [28] further showed that a rapid increase in the number of non-breeders and frequent nest disturbances may also negatively affect early pup survival in captive colonies. It is, therefore, unclear what the net effect of non-breeders on pup survival is and how it compares with the effect of maternal care.
The aim of this study was to quantify early pup survival and identify factors influencing pup survival in the first 10 days postpartum in a eusocial breeding mammal, using data from seven years of observing captive naked mole-rat colonies. While the demographic and environmental conditions in captivity may not necessarily mimic those in the wild, the subterranean lifestyle would make such a study difficult, if not impossible, to replicate in the wild. In particular, the detailed monitoring of our captive colonies enabled us to disentangle the effects of (1) individual pup characteristics (sex, birth body mass), (2) maternal characteristics (body mass, number of mammae), (3) social characteristics (colony size, litter size) and (4) environmental characteristics (ambient temperature, ambient humidity, number of nestbox changes).
Table 1. Effect of different variables on pup survival in mammals. The effect of variables on pup survival is denoted as positive (+), negative (−) and no significant effect (n.e.); studies were conducted under laboratory conditions (captive), natural conditions (wild) or fenced outdoor enclosures (semi-captive).
Table 1. Effect of different variables on pup survival in mammals. The effect of variables on pup survival is denoted as positive (+), negative (−) and no significant effect (n.e.); studies were conducted under laboratory conditions (captive), natural conditions (wild) or fenced outdoor enclosures (semi-captive).
EffectSpeciesConditionSources
Pup birth body mass+Fukomys mechowiicaptive[1]
+Ctenomys talarumcaptive[2]
+Peromyscus californicuscaptive[38]
+Oryctolagus cuniculussemi-captive[6]
Pup age+Leptonychotes weddelliiwild[39]
+Fukomys mechowiicaptive[1]
+Meriones unguiculatuscaptive[31]
+Mus musculus L.wild[40]
+Oryctolagus cuniculussemi-captive[6]
+Heterocephalus glabercaptive[41,42]
Maternal body mass+Clethrionomys glareoluscaptive[4]
+Oryctolagus cuniculussemi-captive[6]
n.e.Octodon deguscaptive[43]
n.e.Sigmodon hispiduswild/captive[44]
n.e.Peromyscus polionotuscaptive[10]
Maternal age+Leptonychotes weddelliiwild[39,45]
Heterocephalus glabercaptive[41]
Maternal behavior+Peromyscus polionotuscaptive[10]
+Crocuta crocutawild[3]
Colony size+Marmota marmotawild[46]
+Lycaon pictuswild[5]
+Suricata suricattawild[47]
Meriones unguiculatuscaptive[7]
Arctocephalus gazellawild[48]
Suricata suricattawild[47]
n.e.Fukomys mechowiiwild[1]
n.e.Meriones unguiculatuswild[49]
n.e.Canis lycaonwild[50]
n.e.Meles meleswild[51]
n.e.Myotis myotiswild[52]
n.e.Mirounga leoninawild[53]
Litter size+Rattus norvegicuscaptive[54]
Meriones unguiculatuscaptive[31]
n.e.Fukomys mechowiiwild[1]
n.e.Clethrionomys glareoluscaptive[4]
n.e.Octodon deguscaptive[43]
n.e.Mus musculuscaptive[55]
n.e.Peromyscus leucopuswild[56]
n.e.Marmota marmotawild[46]
n.e.Sus scrofa domesticuscaptive[57]
n.e.Canis lycaonwild[50]
Ambient temperature+Meriones unguiculatuscaptive[49]
n.e.Heterocephalus glabercaptive[58]
n.e.Rhabdomys pumiliocaptive[59]
Handling/disturbanceCryptomys hottentotuscaptive[60]
n.e.Heterocephalus glabercaptive[37]
n.e.Rattus norvegicuscaptive[61]
n.e.Mus musculuscaptive[55]

2. Materials and Methods

2.1. Study Species

Naked mole-rats are eusocial mammals that live in extensive, subterranean burrows with a mean colony size of 74 ± 49 individuals (n = 30 colonies) [62]. A colony usually contains one breeding female, the queen, who monopolizes mating by reproductively suppressing other colony members and is solely responsible for producing the colony offspring [19,63,64]. Naked mole-rats are long-lived, with a maximum recorded lifespan of 37 years [65]. Queens form long-lasting pair bonds with their breeding male(s) and are capable of sustaining reproduction for at least 27 years [66]. After a gestation of 72–77 days [35] the queen gives birth to 12 ± 5.7 pups (n = 84 litters, 1–29 pups) in the colony nest [28,41], in up to five litters per year [67]. Pups are born blind, with little locomotor activity and are kept warm in the communal nest by other mole-rats (Figure 1) [34,68]. Pups are exclusively nursed by the queen until weaning at ~40 days [28]. One week before they start to eat solid foods, pups begin to beg for caecotrophs from all colony members [34,37]. Naked mole-rats do not exhibit pronounced incest avoidance and may inbreed to a high degree (r = 0.63–0.88) [69,70] although they have a preference for mating with unfamiliar males [71]. Non-breeders usually are close relatives of the queen and her offspring and may, thus, increase their inclusive fitness by staying in the colony throughout their life and helping to rear the offspring [72].

2.2. Housing Conditions

The first naked mole-rat breeding stock (KI) used in this study was purchased from the ABQ BioPark Zoo (Albuquerque, NM, USA) in 2008. Based on this colony, nine more colonies (KII, W, I, H, L, J, K, M, G) were founded at our institute. In 2012, we purchased two other colonies—one colony (C) from the Osnabrück Zoo (Osnabrück, Germany) and one colony (N) from the Schönbrunn Zoo (Vienna, Austria). The breeding stock of colony (D) was composed of individuals from the colonies (KI and N), and colony (E) was founded by animals from the colonies (KI and C). This study presents data from 14 captive colonies closely observed for a total of seven years (September 2008–December 2015), each comprising a breeding pair as well as 29 ± 18 non-breeders. Observation periods of the colonies were as follows: 7–18 months (colonies: L, I, G, W, K, J, H, E, M), 34–39 months (colonies: D, C, N) and 76–87 months (colonies: KI, KII). Each colony was separately kept in an artificial burrow system inside a climatized box consisting of 5–8 acrylic glass boxes of 25 cm diameter and interconnecting tunnels. The ambient temperature was adjusted to 28.0 ± 1.25 °C via a heating coil around the burrow system and to a relative ambient humidity of 61.0 ± 13.8% using wet tissue paper inside the box lids, simulating environmental conditions similar to their native habitat [73,74]. The chambers were organized by the colonies themselves as nestbox, toilet(s) or food-store and contained wood shavings for bedding, unbleached paper tissue as nesting material and wooden branches and blocks for enrichment. Depending on the degree of soiling, boxes were cleaned and partly refreshed every day or every few days. Naked mole-rats were fed ad libitum with a mixed diet of vegetables and fruits, supplemented with high protein and vitamin enriched cereals. No free water was provided. All naked mole-rats were kept and bred at the Leibniz Institute for Zoo and Wildlife Research (IZW) in Berlin, with the approval by the local ethics committee of the ‘Landesamt für Gesundheit und Soziales’, Berlin, Germany (reference no. #ZH 156; 23 September 2008).

2.3. Data Collection

By inspecting all colonies for signs of delivery and newborn litters several times each day, we could record the birth of all litters within 12 h of postnatal life (Day 0). From the day of birth until an age of 80 days, all newborns were daily taken out of the nest, counted and weighed using a digital top-loading balance with a precision of 0.1 g (Kern 822–37, Gottl. Kern and Sohn, Albstadt, Germany). The colonies appeared relatively undisturbed by the regular checking procedure used to monitor the pups since no queen abandoned any of her neonates. Naked mole-rats remain, however, extremely sensitive to vibrations in and around their burrow system and, in our holding facility, responded to disturbances by hastily leaving the nest and ceasing their current behavior. If disturbances are too frequent, colonies usually move their communal nest to another place in the burrow. Thus, we used nestbox changing behavior within a colony as an indicator of stressful environment. To do so, we recorded every day (before the colonies were checked and cleaned) the location of the nestbox and toilet(s). We also recorded the ambient temperature and humidity in each colony system as proxies of the environmental conditions.
All newborns were marked within 12 h after birth, as previously described [28,75], and individual tissue samples were stored for later DNA extraction. Offspring sex was determined by a modified PCR protocol [76] following Szafranski et al. [77]. For permanent marking, all surviving pups received a transponder microchip at the age of three months (size 7 × 1 mm). The gestation period was 73.0 ± 3.6 days (n = 11 observed copulations) with a 6–11 day postpartum estrus (this study) [35]. If conception was not directly observed (n = 68 out of 79 litters) during behavioral estrus (2–24 h [28]), the conception date was defined as being 73 days before the birth of a litter. At conception, queens were weighed and their mammae numbers counted.

2.4. Statistical Analysis

The main goal of this study was to assess potential predictors of early pup survival in naked mole-rats (Table 2). The first days after birth appear to be the most critical days for pup survival (Figure 2B) [28,34,36]. We, therefore, summarized pup survival with a binary variable indicating whether neonates survived past their 10th day, including day 10, or not (i.e., scoring 1 when the pup survived from day 0 till 9, including day 9, and 0 otherwise). We used this variable as the response variable in a logistic regression model. We examined four categories of effects encompassing nine variables as predictors in the model. These predictors were organized as follows: (1) individual characteristics (pup sex and birth body mass), (2) maternal characteristics (body mass and number of mammae), (3) social characteristics (colony size, litter size) and (4) environmental characteristics (ambient temperature, ambient humidity, number of nestbox changes). All variables except pup sex were measured as quantitative variables. Other predictors were initially considered, including litter order, maternal body length, colony age, burrow length or number of related naked mole-rats to the queen, but they turned out to be highly collinear with maternal body mass and colony size and less revealing than the finally chosen predictors. The fitting procedure estimated a total of 10 parameters for the fixed effects (1 for the intercept and 9 for modeling the effect of the predictor variables). The analysis was performed at the level of the pup. After removing missing values from the entire dataset (16 queens and the 79 litters), 585 out of 869 pups of 15 queens from 57 litters were considered in these analyses. To account for these dependencies in the analysis, we considered both the identity of the mother and the identity of the litter as two random effects with a Gaussian distribution, which led to the estimation of two additional parameters, the two corresponding variances. We fitted the model using the function ‘fitme’ from the package ‘spaMM’ (version 4.1.2.) [78] in the program R (version 4.2.2.) [79]. To assess the statistical significance of each predictor, we carried out a likelihood ratio test comparing the goodness of fit between the full model and a model refitted after excluding the variable under consideration. The significance of the likelihood ratio tests was assessed by parametric bootstraps, which represent a more robust alternative to the traditional asymptotic likelihood ratio test. As estimates for fixed effects are expressed on the logit scale in logistic regressions, we converted the estimates into odds ratios in the text to facilitate interpretation. P-values are computed for 2-tailed tests. Correlations were calculated accordingly as Spearman’s rank correlations.

3. Results

3.1. Description of Early Pup Survival

During the seven years of study, 869 pups (79 litters) were born to 16 queens in 14 colonies (Figure 2A). The proportion of pups that lived more than 10 days was 57.0% (491 of 862 pups for which survival at this age was known). Survival did not change much after this point within the duration of the study, dropping to 50.6% (436 of 862 pups) by the age of 60 days (Figure 2B) and, therefore, corresponding to a survival rate of 88.8% (436/491) between day 11 and day 60, or an average 10-day survival rate of 97.6% (=((436/491)^(1/(60−11)))^10). Until day 10, all pups survived in 30.8% of the litters (24 of 78 litters for which survival was known for all pups) and all pups died in 14.1% of the litters (11 of 78 litters; Figure 2A). Seven queens successfully raised 100% of their litters and pups (n = 14 litters). The remaining nine queens had together at least one or more surviving pups in 83.1% of their litters (54 out of 65 litters). Four of these nine queens had 11 single litters in which all pups died (Figure 2A).

3.2. Predictors of Early Pup Survival

The statistical analysis of survival to the age of 10 days showed that individual, maternal and social characteristics all significantly influenced early pup survival (Table 2). Environmental characteristics had no significant effects (Table 2). The predictors with the largest effect sizes were pup birth body mass, maternal body mass, number of mammae and colony size (Table 2). Recomputing the regression slopes after standardizing these quantitative variables (with mean = 0, sd = 1) indicated that pup birth body mass showed the largest absolute standardized effect size (β = 2.06), followed by number of mammae (β = 1.56), colony size (β = −1.24) and maternal body mass (β = −1.00).

3.2.1. Effect of Individual Characteristics

Individual pup characteristics had a highly significant effect on early pup survival (Table 2). Amongst these, pup birth body mass was the best predictor of early pup survival. Newborns (n = 843) weighed on average 1.97 ± 0.31 g (range 1.03 to 3.07 g). Heavier pups had a significantly higher chance to survive the first 10 days (Table 2, Figure 3C). Surviving pups had a birth body mass of at least 1.34 g. The model predicted that if birth body mass of pups would be increase by 0.5 g, the probability to survive would improve 27-fold. Pup sex did not significantly influence pup survival (Table 2).

3.2.2. Effect of Maternal Characteristics

The survival of pups was significantly affected by characteristics of the mother (Table 2). The body masses of queens at conception varied widely from 27 to 65 g (50.0 ± 7.9 g). All breeding females showed fluctuating body masses at conception, with mass changes from −5 to +16 g (4.6 ± 4.2 g). Queens with a high body mass at conception were significantly less successful at raising pups than lighter queens (Table 2, Figure 3B). The model predicts that if maternal body mass would be decrease by 10 g, the probability for the pups to survive would improve 5-fold. Maternal body mass increased with the number of litters produced, a possible proxy for maternal experience (Spearman’s rank correlation, ρ = 0.375, n = 58, p < 0.001). Maternal body mass also increased with age (ρ = 0.486, n = 58, p < 0.001). All mammae were functional and their number remained constant for each queen. The number of mammae ranged from 11 to 13 (11.7 ± 0.9). For the 15 queens for which we knew the number of mammae, 13 had an odd number of mammae (86.70%) and 2 had an even number (13.3%). In 41.0% of the litters (n = 31 litters), the number of littermates exceeded or was equal (9.0%; n = 7 litters) to the number of mammae. Number of mammae positively influenced pup survival (Table 2, Figure 3A). The model predicts that with each additional mammae the probability for the pups to survive improve 9-fold.

3.2.3. Effect of Social Characteristics

Pup survival significantly decreased with colony size (Table 2, Figure 3D). In our 14 established colonies, colony size ranged from 2 to 55 individuals. An increase in colony size and density resulted in decreased pup survival over the first 10 days. For each 10 additional mole-rats within the pre-defined size of the burrow system, our model predicts that pup survival would decrease 2-fold. The model further forecasts that pup survival dropped below 50.0% in colonies with more than 37 individuals (Figure 3D). Litter size ranged from 1 to 22 pups (11.0 ± 4.4), but we found no significant relationship between litter size and early pup survival (Table 2).

3.2.4. Effect of Environmental Characteristics

Depending on the number of individuals in a colony, ambient temperature ranged from 24 °C to 30 °C (28 ± 1.25 °C). Ambient humidity varied between 10 and 90% (61.0 ± 13.8%) across colonies and during this study. The frequency of nestbox changes varied between colonies from 1 to 6 (2.1 ± 0.7) changes per week during the first 10 days postpartum. Despite such variation in these environmental variables, we were unable to detect any significant effect on early pup survival (Table 2).

4. Discussion

Our study of the determinants of early pup survival in naked mole-rats show that survival is variable and lowest in the first postnatal days (Figure 2B), similar to what was reported for small mammals [4,30,31,32], mole-rats [1,80,81] and other naked mole-rat colonies [28,34,36,82]. Body mass at birth showed the strongest effect on early pup survival (Figure 3C), with heavier pups exhibiting a higher chance to survive the first postnatal days than individuals with lower birth body mass. Comparable results (Table 1) were also reported for many other small mammals [1,2,6]. In naked mole-rats, all colony members tend to sleep and huddle together in the nest, with the breeding female typically assuming the nursing position at the top of the pile (Figure 1). Heavier naked mole-rat pups were observed to reach and remain on top of the pile more easily, therefore, close to the queen and her teats (Figure 1). This advantageous position may reduce their likelihood of being crushed or suffocated by other mole-rats in the nest.
Maternal characteristics also appear to be important predictors of early pup survival in naked mole-rats (Table 2). For eusociality to evolve in mammals, queens must allocate large amounts of energy and nutrients to offspring production during pregnancy and lactation. Hood et al. [12] showed that the milk of naked mole-rat queens is relatively dilute and low in energy compared to that of other rodents. Therefore, given the high reproductive output of breeding naked mole-rat queens (litter size: 1–29 pups) [83] and the absence of food provisioning by non-breeders before weaning, queens must produce an usually large amount of milk to meet nutritive requirements of their pups. Hood et al. [12] estimated, on the basis of pup masses and average litter sizes, a daily milk yield for lactating naked mole-rat queens of 30.5 g, equivalent to 58.0% of maternal mass. This is a large maternal investment, as a comparison with other rodents shows. For instance, a lactating mouse only produces an approximate daily milk yield of 14 g, equivalent to 35.0% of its body mass [84]. The view that maternal investment in milk production is critical in naked mole-rats is further supported by our finding that lactating queens with a higher number of mammae also had higher pup survival (Figure 3A). Among mammals in general, and rodents in particular, mean litter sizes usually are approximately one-half the number of mammae, and maximum litter sizes approximate mammary numbers [83,85]. Naked mole-rats appear to be an exception to this rule, as breeding females exhibit a maximum number of 14 mammae but can have litters of up to 29 pups [83]. In this study, in 50.7% (n = 38) of litters, the number of littermates was equal to or above the maternal number of mammae. So, as in other mammals but unlike eusocial insects, the presence of maternal care provided to pups in early life seems critical in naked mole-rats.
A high body mass in lactating queens may increase the amount of milk produced and, thus, could be advantageous for offspring survival, as seen in other mammalian species [4,6]. This idea is not supported by our results, which instead showed that queens with a higher body mass were less successful at rearing pups than light mothers (Figure 3B). The negative effect of maternal body mass may be real, but it could also be indirectly caused by the effect of other maternal characteristics.
A prime suspect for such a confounding variable is age. For age to qualify as a confounding variable, variation in age should be present between queens, larger queens should be older, and older queens should provide less maternal care to their pups. These three conditions seem to be satisfied in naked mole-rats. First, although we know the exact age of only some of our queens, estimating queen ages on the basis of morphological characteristics [41,86] revealed that their age ranged from 1 to 20 years. Second, we observed that the body mass of queens was positively correlated to the age estimate, with almost all queens experiencing a gain in body mass, averaging 12.0% over the course of the study. Naked mole-rat queens are expected to gain mass as they age because they generally undergo fecundity-enhancing vertebral elongation, another unique characteristic of their biology [24,25,26]. Our study contributes little information about the third condition—that maternal care may be lower in older queens. We did not discover obvious lactation problems in any of our queens, but it remains possible that adequate milk supply decreased with maternal age. Buffenstein [41] reported that elderly naked mole-rat queens tended to have very large litters and that pup survival was extremely low, possibly because of inadequate maternal milk production. In sum, it is likely that aging is at least partly responsible for the observed negative effect of maternal body mass on early pup survival. Future studies could clarify whether age is the sole factor explaining why heavier queens produce pups that are less likely to survive than lighter queens. We recommend future studies to also record aggressive behaviors since heavier breeding females in rodents are sometimes more aggressive [87], and queen aggressiveness may detrimentally affect pup survival in naked mole-rats [28].
The effect of alloparental behavior and the effect of the number of non-breeders on offspring survival are quite variable across cooperatively breeding species. Some studies show a positive link [5,46,88], others no correlation [1,47] or even a negative effect [7,89]. In naked mole-rats, the effect of the number of non-breeders (i.e., colony size) on pup survival had not been quantified before. Alloparental care is not essential for early pup survival in naked mole-rats, since a breeding pair alone can rear pups to weaning in both the lab (this study: pup survival rate when raised by pairs = 89.69%; [28]) and the field [62].
Nevertheless, since naked mole-rats are highly inbred [70], wild colonies contain a large number of individuals up to 295 [90] and non-breeders are reproductively suppressed, often philopatric throughout their life and while performing pup-care, it seems reasonable to expect the number of non-breeders to have a positive effect on pup survival as demonstrated for the closely related Damaraland mole-rat (Fukomys damarensis: Jarvis and Bennett, 1993 unpubl.).
In our study, however, an increase in the number of non-breeders significantly reduced early pup survival (Figure 3D). Our results may be influenced by the laboratory settings. As also reported by other researchers, captive naked mole-rat colonies often cease pup recruitment after an initial phase of high pup survival, and litter survivorship remains higher when some mole-rats are regularly removed from the colony [28,62]. Perhaps negative density dependence emerges because naked mole-rats tend to huddle all together in a single nest with the lactating queen on the top of the pile (Figure 1). Regardless of the actual colony size, all our colonies had the same nestbox size, which implies that the operational density at the nest necessarily increased with colony size. With an increasing number of mole-rats, pups were observed more frequently to be displaced to the bottom of the pile, and thus, away from the preferred nursing position, or crushed in the nest by other colony members. As studies of radio-tracking wild naked mole-rat demonstrated that colony nests vary considerably in shape and size and that mole-rats use several nest-sites at the same time [90], this phenomenon is likely to be an artifact of captivity.
A second reason for the higher pup mortality rate with an increasing number of non-breeders may be the fact that non-breeding mole-rats are also frequently observed to carry out negative pup-care behaviors [34,35,91]. These behaviors cause a substantial number of offspring deaths and often result in the failure of the whole litter (this study) [28]. In unsuccessful litters, pups are constantly dumped, kicked and dragged out of the nest until they die from physical trauma [28,35]. Colony members were also seen to kill and consume pups (this study) [28,34,91]. The motivation or purpose of these negative pup care behaviors by older colony members in a eusocial species are still unexplored but do not appear to be driven by the environmental stress factors that we recorded [91]. Whether these behaviors occur in the wild or only in captivity remains to be determined.
A third explanation for the negative effect of the number of non-breeders on pup survival could be that the positive effects of these individuals are lower in captivity than in the wild. In wild colonies, the breeding pair and their newborns may benefit from behaviors that occur little in captivity because they are irrelevant or have little or no effect in captivity. In particular, in the wild, non-breeders perform the high risk tasks of foraging, defending the colony and maintaining the tunnel system [27]—all behaviors that are of no relevance in an artificial burrow system.
A fourth explanation could be that the effect of non-breeders may vary during different offspring life stages. Clutton-Brock et al. [47], for example, showed that in meerkats (Suricata suricatta), pup mortality was not significantly related to group size between birth and emergence but was important for the time between emergence and five months of age. In our study, the survival of naked mole-rat pups was negatively correlated with colony size in the first 10 days postpartum, so it remains possible that positive effects of non-breeders may be received by pups at a later age.
A fifth and final possible explanation is that the lower rate of offspring survival at a higher density of naked mole-rats in the burrow system may also be of adaptive value. While small and recently founded colonies will likely need to acquire a large workforce as fast as possible to be able to defend and expand the burrow system to find sufficient food, a large and established colony in an environment with scarce resources may benefit from a high pup mortality that retains only the fittest pups and decreases competition.

5. Conclusions

Our findings highlight the diversity of factors influencing early pup survival in a eusocial mammal. Naked mole-rats combine the characteristics of cooperatively breeding mammals and eusocial insects. Therefore, as expected, both maternal and social predictors exhibited a significant effect on early pup survival. Surprisingly, maternal body mass and colony size had a significant negative effect on early pup survival. Our results thus reveal that studies performed in captivity may underestimate the benefits provided by non-breeders and overestimate the costs of increasing colony size unless variable nestbox sizes are available.

Author Contributions

Conceptualization, M.W., A.C. and H.H.; methodology, M.W. and A.C.; formal analysis, M.W. and A.C.; investigation, M.W. and S.H.; data curation, M.W.; writing—original draft preparation, M.W.; writing—review and editing, A.C.; funding acquisition, T.B.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Leibniz-Institute for Zoo and Wildlife Research and the Leibniz Association (SAW-2012-FLI-2). The publication of this article was funded by the Open Access Fund of the Leibniz Association and by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation)—Project number 491292795.

Institutional Review Board Statement

The animal protocol was approved by the local ethics committee of the ‘Landesamt für Gesundheit und Soziales’, Berlin, Germany (reference no. #ZH 156; 23 September 2008).

Informed Consent Statement

Not applicable.

Data Availability Statement

All data collected are available in the present paper.

Acknowledgments

We are grateful for funding from the Leibniz-Institute for Zoo and Wildlife Research and the Leibniz Association (SAW-2012-FLI-2) and the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation). We would like to thank Jette Ziep, Nga Nguyen, Rafaela Fiúza and Melanie Schwochow for expert technical assistance during this project.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Scharff, A.; Begall, S.; Grütjen, O.; Burda, H. Reproductive characteristics and growth of Zambian giant mole-rats, Cryptomys mechowi (Rodentia: Bathyergidae). Mammalia 1999, 63, 217–230. [Google Scholar] [CrossRef]
  2. Zenuto, R.R.; Antinuchi, C.D.; Busch, C. Bioenergetics of Reproduction and Pup Development in a Subterranean Rodent (Ctenomys talarum). Physiol. Biochem. Zool. 2002, 75, 469–478. [Google Scholar] [CrossRef] [PubMed]
  3. Hofer, H.; East, M.L. Behavioral processes and costs of co-existence in female spotted hyenas: A life history perspective. Evol. Ecol. 2003, 17, 315–331. [Google Scholar] [CrossRef]
  4. Jonsson, P.; Agrell, J.; Koskela, E.; Mappes, T. Effects of litter size on pup defence and weaning success of neighbouring bank vole females. Can. J. Zool. 2002, 80, 1–5. [Google Scholar] [CrossRef]
  5. Creel, S.; Creel, N.M. Opposing effects of group size on reproduction and survival in African wild dogs. Behav. Ecol. 2015, 26, 1414–1422. [Google Scholar] [CrossRef]
  6. Rödel, H.; Starkloff, A.; Seltmann, M.W.; Prager, G.; Von Holst, D. Causes and predictors of nest mortality in a European rabbit population. Mamm. Biol. 2009, 74, 198–209. [Google Scholar] [CrossRef]
  7. Scheibler, E.; Weinandy, R.; Gattermann, R. Social factors affecting litters in families of Mongolian gerbils, Meriones unguiculatus. Folia Zool. 2005, 54, 61. [Google Scholar]
  8. Hofer, H.; East, M. Population dynamics, population size, and the commuting system of Serengeti spotted hyenas. In Serengeti II: Dynamics, Management, and Conservation of an Ecosystem; Sinclaire, A.R.E., Arcese, P., Eds.; University of Chicago Press: Chicago, IL, USA, 1995; p. 332. [Google Scholar]
  9. Maestripieri, D.; Mateo, J.M. Maternal Effects in Mammals; University of Chicago Press: Chicago, IL, USA, 2009. [Google Scholar]
  10. Margulis, S.W.; Nabong, M.; Alaks, G.; Walsh, A.; Lacy, R.C. Effects of early experience on subsequent parental behaviour and reproductive success in oldfield mice, Peromyscus polionotus. Anim. Behav. 2005, 69, 627–634. [Google Scholar] [CrossRef]
  11. Hofer, H.; Benhaiem, S.; Golla, W.; East, M.L. Trade-offs in lactation and milk intake by competing siblings in a fluctuating environment. Behav. Ecol. 2016, 27, 1567–1578. [Google Scholar] [CrossRef]
  12. Hood, W.R.; Kessler, D.S.; Oftedal, O.T. Milk composition and lactation strategy of a eusocial mammal, the naked mole-rat. J. Zool. 2014, 293, 108–118. [Google Scholar] [CrossRef]
  13. Jennions, M.D. Cooperative breeding in mammals. Trends Ecol. Evol. 1994, 9, 89–93. [Google Scholar] [CrossRef] [PubMed]
  14. Ciszek, D. Behavior of the Naked Mole-Rat (Heterocephalus glaber), with Emphasis on Factors Related to the Attainment of Breeding Status. Ph.D. Thesis, University of Michigan, Ann Arbor, MI, USA, 2001. [Google Scholar]
  15. Tardif, S.D.; Moehlman, P.D.; Hofer, H.; Asa, C.S.; Creel, S.R.; Waser, P.M.; Lucas, J.R.; Getz, L.L.; Carter, C.S.; Roberts, R.L.; et al. Cooperative Breeding in Mammals; Solomon, N.G., French, J.A., Eds.; Cambridge University Press: Cambridge, UK, 1997. [Google Scholar] [CrossRef]
  16. Crespi, B.J.; Yanega, D. The definition of eusociality. Behav. Ecol. 1995, 6, 109–115. [Google Scholar] [CrossRef]
  17. Bennett, N.C.; Faulkes, C.G. African Mole-Rats: Ecology and Eusociality; Cambridge University Press: Cambridge, UK, 2000. [Google Scholar]
  18. Burda, H.; Honeycutt, R.L.; Begall, S.; Locker-Grütjen, O.; Scharff, A. Are naked and common mole-rats eusocial and if so, why? Behav. Ecol. Sociobiol. 2000, 47, 293–303. [Google Scholar] [CrossRef]
  19. Jarvis, J.U.M. Eusociality in a Mammal: Cooperative Breeding in Naked Mole-Rat Colonies. Science 1981, 212, 571–573. [Google Scholar] [CrossRef]
  20. Jarvis, J.U.M.; Bennett, N. Eusociality has evolved independently in two genera of bathyergid mole-rats—But occurs in no other subterranean mammal. Behav. Ecol. Sociobiol. 1993, 33, 253–260. [Google Scholar] [CrossRef]
  21. Batra, S.W. Nests and social behavior of halictine bees of India (Hymenoptera: Halictidae). Indian J. Entomol. 1966, 28, 375. [Google Scholar]
  22. Michener, C.D. Comparative Social Behavior of Bees. Annu. Rev. Entomol. 1969, 14, 299–342. [Google Scholar] [CrossRef]
  23. Wilson, E.O. The Insect Societies; Harvard University Press: Cambridge, MA, USA, 1971. [Google Scholar]
  24. Dengler-Crish, C.M.; Catania, K.C. Phenotypic plasticity in female naked mole-rats after removal from reproductive suppression. J. Exp. Biol. 2007, 210, 4351–4358. [Google Scholar] [CrossRef]
  25. Henry, E.C.; Crish, C.; Catania, K.C. Growing out of a caste—Reproduction and the making of the queen mole-rat. J. Exp. Biol. 2007, 210, 261–268. [Google Scholar] [CrossRef]
  26. O’Riain, M.J.; Jarvis, J.U.M.; Alexander, R.; Buffenstein, R.; Peeters, C. Morphological castes in a vertebrate. Proc. Natl. Acad. Sci. USA 2000, 97, 13194–13197. [Google Scholar] [CrossRef]
  27. Braude, S. Dispersal and new colony formation in wild naked mole-rats: Evidence against inbreeding as the system of mating. Behav. Ecol. 2000, 11, 7–12. [Google Scholar] [CrossRef]
  28. Jarvis, J.U.M. Reproduction of Naked Mole-Rats. In The Biology of the Naked Mole-Rat; Sherman, P.W., Jarvis, J.U.M., Alexander, R.D., Eds.; Princeton University Press: Princeton, NJ, USA, 1991; pp. 384–425. [Google Scholar] [CrossRef]
  29. Buffenstein, R. Negligible senescence in the longest living rodent, the naked mole-rat: Insights from a successfully aging species. J. Comp. Physiol. B 2008, 178, 439–445. [Google Scholar] [CrossRef]
  30. Hallett, A.F.; Meester, J. Early Postnatal Development of the South African Hamster Mystromys Albicaudatus. Zool. Afr. 1971, 6, 221–228. [Google Scholar] [CrossRef]
  31. Norris, M.L.; Adams, C.E. Mortality from birth to weaning in the mongolian gerbil, Meriones unguiculatus. Lab. Anim. 1972, 6, 49–53. [Google Scholar] [CrossRef]
  32. Rödel, H.G.; Starkloff, A.; Bautista, A.; Friedrich, A.-C.; Von Holst, D. Infanticide and Maternal Offspring Defence in European Rabbits under Natural Breeding Conditions. Ethology 2008, 114, 22–31. [Google Scholar] [CrossRef]
  33. Edrey, Y.H.; Hanes, M.; Pinto, M.; Mele, J.; Buffenstein, R. Successful Aging and Sustained Good Health in the Naked Mole Rat: A Long-Lived Mammalian Model for Biogerontology and Biomedical Research. ILAR J. 2011, 52, 41–53. [Google Scholar] [CrossRef] [PubMed]
  34. O’Riain, J. Pup Ontogeny and Factors Influencing Behavioural and Morphological Variation in Naked Mole-Rats, Heterocephalus glaber (Rodentia, Bathyergidae). Ph.D. Thesis, University of Cape Town, Cape Town, South Africa, March 1996. [Google Scholar]
  35. Lacey, E.A.; Sherman, P.W. Social Organization of Naked Mole-Rat Colonies: Evidence for Divisions of Labor. In The Biology of the Naked Mole-Rat; Sherman, P.W., Jarvis, J.U.M., Alexander, R.D., Eds.; Princeton University Press: Princeton, NJ, USA, 1991; pp. 275–336. [Google Scholar] [CrossRef]
  36. Ke, Z.; Vaidya, A.; Ascher, J.; Seluanov, A.; Gorbunova, V. Novel husbandry techniques support survival of naked mole rat (Heterocephalus glaber) pups. J. Am. Assoc. Lab. Anim. Sci. 2014, 53, 89–91. [Google Scholar] [PubMed]
  37. Payne, S.F. Social Organization of the Naked Mole-Rat (Heterocephalus Glaber): Cooperation in Colony Labor and Reproduction. Master Thesis, University of California, Santa Cruz, CA, USA, December 1982. [Google Scholar]
  38. Dudley, D. Contributions of paternal care to the growth and development of the young in Peromyscus californicus. Behav. Biol. 1974, 11, 155–166. [Google Scholar] [CrossRef]
  39. Proffitt, K.M.; Rotella, J.J.; Garrott, R.A. Effects of pup age, maternal age, and birth date on pre-weaning survival rates of Weddell seals in Erebus Bay, Antarctica. Oikos 2010, 119, 1255–1264. [Google Scholar] [CrossRef]
  40. Brown, R.Z. Social Behavior, Reproduction, and Population Changes in the House Mouse (Mus musculus L.). Ecol. Monogr. 1953, 23, 217–240. [Google Scholar] [CrossRef]
  41. Buffenstein, R. The Naked Mole-Rat: A New Long-Living Model for Human Aging Research. J. Gerontol. Ser. A 2005, 60, 1369–1377. [Google Scholar] [CrossRef]
  42. O’Riain, M.J.; Jarvis, J.U.M.; Faulkes, C.G. A dispersive morph in the naked mole-rat. Nature 1996, 380, 619–621. [Google Scholar] [CrossRef]
  43. Ebensperger, L.A.; Hurtado, M.J.; León, C. An experimental examination of the consequences of communal versus solitary breeding on maternal condition and the early postnatal growth and survival of degu, Octodon degus, pups. Anim. Behav. 2007, 73, 185–194. [Google Scholar] [CrossRef]
  44. Campbell, M.T.; Slade, N.A. The effect of maternal mass on litter size and offspring survival in the hispid cotton rat (Sigmodon hispidus). Can. J. Zool. 1995, 73, 133–140. [Google Scholar] [CrossRef]
  45. Hastings, K.K.; Testa, J.W. Maternal and birth colony effects on survival of Weddell seal offspring from McMurdo Sound, Antarctica. J. Anim. Ecol. 1998, 67, 722–740. [Google Scholar] [CrossRef]
  46. Allainé, D.; Brondex, F.; Graziani, L.; Coulon, J.; Till-Bottraud, I. Male-biased sex ratio in litters of Alpine marmots supports the helper repayment hypothesis. Behav. Ecol. 2000, 11, 507–514. [Google Scholar] [CrossRef]
  47. Clutton-Brock, T.H.; Gaynor, D.; McIlrath, G.M.; Maccoll, A.D.C.; Kansky, R.; Chadwick, P.; Manser, M.; Skinner, J.D.; Brotherton, P.N.M. Predation, group size and mortality in a cooperative mongoose, Suricata suricatta. J. Anim. Ecol. 1999, 68, 672–683. [Google Scholar] [CrossRef]
  48. Doidge, D.W.; Croxall, J.P.; Baker, J.R. Density-dependent pup mortality in the Antarctic fur seal Arctocephalus gazella at South Georgia. J. Zool. 1984, 202, 449–460. [Google Scholar] [CrossRef]
  49. French, J.A. Alloparents in the Mongolian gerbil: Impact on long-term reproductive performance of breeders and opportunities for independent reproduction. Behav. Ecol. 1994, 5, 273–279. [Google Scholar] [CrossRef]
  50. Mills, K.J.; Patterson, B.R.; Murray, D.L. Direct Estimation of Early Survival and Movements in Eastern Wolf Pups. J. Wildl. Manag. 2008, 72, 949–954. [Google Scholar] [CrossRef]
  51. Woodroffe, R.; Macdonald, D.W. Helpers provide no detectable benefits in the European badger (Meles meles). J. Zool. 2000, 250, 113–119. [Google Scholar] [CrossRef]
  52. Zahn, A. Reproductive success, colony size and roost temperature in attic-dwelling bat Myotis myotis. J. Zool. 1999, 247, 275–280. [Google Scholar] [CrossRef]
  53. Pistorius, P.A.; Bester, M.N.; Kirkman, S.P.; Taylor, F.E. Pup mortality in southern elephant seals at Marion Island. Polar Biol. 2001, 24, 828–831. [Google Scholar] [CrossRef]
  54. Kahan, E.; Rosen, M. The influence of litter size and favourable conditions on mortality and weaning weights in rats. Lab. Anim. 1984, 18, 247–251. [Google Scholar] [CrossRef] [PubMed]
  55. Gates, W.H. Litter size, birth weight, and early growth rate of mice (Mus musculus). Anat. Rec. 1925, 29, 183–193. [Google Scholar] [CrossRef]
  56. Goundie, T.R.; Vessey, S.H. Survival and Dispersal of Young White-Footed Mice Born in Nest Boxes. J. Mammal. 1986, 67, 53–60. [Google Scholar] [CrossRef]
  57. Milligan, B.N.; Fraser, D.; Kramer, D.L. Birth weight variation in the domestic pig: Effects on offspring survival, weight gain and suckling behaviour. Appl. Anim. Behav. Sci. 2001, 73, 179–191. [Google Scholar] [CrossRef]
  58. Woodley, R.; Buffenstein, R. Thermogenic changes with chronic cold exposure in the naked mole-rat (Heterocephalus glaber). Comp. Biochem. Physiol. Part A: Mol. Integr. Physiol. 2002, 133, 827–834. [Google Scholar] [CrossRef]
  59. Schradin, C.; Pillay, N. The influence of the father on offspring development in the striped mouse. Behav. Ecol. 2004, 16, 450–455. [Google Scholar] [CrossRef]
  60. Malherbe, G.; Bennett, N.; Schoeman, A. The reproductive biology and postnatal development of the highveld mole-rat, Cryptomys hottentotus pretoriae (Rodentia: Bathyergidae). Afr. Zool. 2003, 38, 161–167. [Google Scholar] [CrossRef]
  61. Burn, C.C.; Mason, G.J. Effects of cage-cleaning frequency on laboratory rat reproduction, cannibalism, and welfare. Appl. Anim. Behav. Sci. 2008, 114, 235–247. [Google Scholar] [CrossRef]
  62. Braude, S. The Behavior and Demographics of the Naked Mole-Rat, Heterocephalus glaber. Ph.D. Thesis, University of Michigan, Ann Arbor, MI, USA, 1991. [Google Scholar]
  63. Edwards, P.D. For the physiology of cooperative breeding, it’s time to move beyond stress: A Comment on: ‘Stress in an underground empire’ by Medger. Biol. Lett. 2022, 18, 20220375. [Google Scholar] [CrossRef] [PubMed]
  64. Faulkes, C.G.; Abbott, D.H. Evidence that primer pheromones do not cause social suppression of reproduction in male and female naked mole-rats (Heterocephalus glaber). Reproduction 1993, 99, 225–230. [Google Scholar] [CrossRef] [PubMed]
  65. Lee, B.P.; Smith, M.; Buffenstein, R.; Harries, L.W. Negligible senescence in naked mole rats may be a consequence of well-maintained splicing regulation. Geroscience 2020, 42, 633–651. [Google Scholar] [CrossRef] [PubMed]
  66. O’Riain, M.J.; Faulkes, C.G. African Mole-Rats: Eusociality, Relatedness and Ecological Constraints. In Ecology of Social Evolution; Korb, J., Heinze, J., Eds.; Springer: Berlin/Heidelberg, Germany, 2008; pp. 207–223. [Google Scholar] [CrossRef]
  67. Brett, R.A. The Population Structure of Naked Mole-Rat Colonies. In The Biology of the Naked Mole-Rat; Sherman, P.W., Jarvis, J.U.M., Alexander, R.D., Eds.; Princeton University Press: Princeton, NJ, USA, 1991; pp. 97–136. [Google Scholar] [CrossRef]
  68. Yahav, S.; Buffenstein, R. Huddling Behavior Facilitates Homeothermy in the Naked Mole Rat Heterocephalus glaber. Physiol. Zool. 1991, 64, 871–884. [Google Scholar] [CrossRef]
  69. Ingram, C.M.; Troendle, N.J.; Gill, C.A.; Braude, S.; Honeycutt, R.L. Challenging the inbreeding hypothesis in a eusocial mammal: Population genetics of the naked mole-rat, Heterocephalus glaber. Mol. Ecol. 2015, 24, 4848–4865. [Google Scholar] [CrossRef]
  70. Reeve, H.K.; Westneat, D.F.; Noon, W.A.; Sherman, P.W.; Aquadro, C.F. DNA “fingerprinting” reveals high levels of inbreeding in colonies of the eusocial naked mole-rat. Proc. Natl. Acad. Sci. USA 1990, 87, 2496–2500. [Google Scholar] [CrossRef]
  71. Clarke, F.M.; Faulkes, C.G. Kin discrimination and female mate choice in the naked mole-rat Heterocephalus glaber. Proc. R. Soc. B Boil. Sci. 1999, 266, 1995–2002. [Google Scholar] [CrossRef]
  72. Hamilton, W. The genetical evolution of social behaviour. II. J. Theor. Biol. 1964, 7, 17–52. [Google Scholar] [CrossRef]
  73. McNab, B.K. The Metabolism of Fossorial Rodents: A Study of Convergence. Ecology 1966, 47, 712–733. [Google Scholar] [CrossRef]
  74. Withers, P.; Jarvis, J. The effect of huddling on thermoregulation and oxygen consumption for the naked mole-rat. Comp. Biochem. Physiol. Part A Physiol. 1980, 66, 215–219. [Google Scholar] [CrossRef]
  75. Röllig, K.; Drews, B.; Göritz, F.; Hildebrandt, T.B. The Long Gestation of the Small Naked Mole-Rat (Heterocephalus glaber RÜPPELL, 1842) Studied with Ultrasound Biomicroscopy and 3D-Ultrasonography. PLoS ONE 2011, 6, e17744. [Google Scholar] [CrossRef] [PubMed]
  76. Katsushima, K.; Nishida, C.; Yosida, S.; Kato, M.; Okanoya, K.; Matsuda, Y. A multiplex PCR assay for molecular sexing of the naked mole-rat (Heterocephalus glaber). Mol. Ecol. Resour. 2010, 10, 222–224. [Google Scholar] [CrossRef] [PubMed]
  77. Szafranski, K.; Wetzel, M.; Holtze, S.; Büntjen, I.; Lieckfeldt, D.; Ludwig, A.; Huse, K.; Platzer, M.; Hildebrandt, T. The Mating Pattern of Captive Naked Mole-Rats Is Best Described by a Monogamy Model. Front. Ecol. Evol. 2022, 10, 855688. [Google Scholar] [CrossRef]
  78. Rousset, F.; Ferdy, J.-B. Testing environmental and genetic effects in the presence of spatial autocorrelation. Ecography 2014, 37, 781–790. [Google Scholar] [CrossRef]
  79. R-Core-Team. R: A Language and Environment for Statistical Computing, 2014; R Foundation for Statistical Computing: Vienna, Austria, 2016. [Google Scholar]
  80. Burda, H. Constraints of pregnancy and evolution of sociality in mole-rats with special reference to reproductive and social patterns in Cryptomys hottentotus (Bathyergidae, Rodentia). J. Zool. Syst. Evol. Res. 1990, 28, 26–39. [Google Scholar] [CrossRef]
  81. Šumbera, R.; Burda, H.; Chitaukali, W.N. Reproductive Biology of a Solitary Subterranean Bathyergid Rodent, the Silvery Mole-Rat (Heliophobius Argenteocinereus). J. Mammal. 2003, 84, 278–287. [Google Scholar] [CrossRef]
  82. Yu, C.; Wang, S.; Yang, G.; Zhao, S.; Lin, L.; Yang, W.; Tang, Q.; Sun, W.; Cui, S. Breeding and Rearing Naked Mole-Rats (Heterocephalus glaber) under Laboratory Conditions. J. Am. Assoc. Lab. Anim. Sci. 2017, 56, 98–101. [Google Scholar]
  83. Sherman, P.W.; Braude, S.; Jarvis, J.U.M. Litter Sizes and Mammary Numbers of Naked Mole-Rats: Breaking the One-Half Rule. J. Mammal. 1999, 80, 720–733. [Google Scholar] [CrossRef]
  84. Knight, C.H.; Maltz, E.; Docherty, A.H. Milk yield and composition in mice: Effects of litter size and lactation number. Comp. Biochem. Physiol. Part A Physiol. 1986, 84, 127–133. [Google Scholar] [CrossRef]
  85. Stockley, P.; Parker, G.A. Life history consequences of mammal sibling rivalry. Proc. Natl. Acad. Sci. USA 2002, 99, 12932–12937. [Google Scholar] [CrossRef]
  86. Braude, S.; Ciszek, D.; Berg, N.E.; Shefferly, N. The ontogeny and distribution of countershading in colonies of the naked mole-rat (Heterocephalus glaber). J. Zool. 2001, 253, 351–357. [Google Scholar] [CrossRef]
  87. Hilakivi-Clarke, L.A.; Lister, R.G. The role of body weight in resident-intruder aggression. Aggress. Behav. 1992, 18, 281–287. [Google Scholar] [CrossRef]
  88. Heymann, E.W.; Soini, P. Offspring number in pygmy marmosets, Cebuella pygmaea *, in relation to group size and the number of adult males. Behav. Ecol. Sociobiol. 1999, 46, 400–404. [Google Scholar] [CrossRef]
  89. Harbo, J.R. Effect of Population Size on Brood Production, Worker Survival and Honey Gain in Colonies of Honeybees. J. Apic. Res. 1986, 25, 22–29. [Google Scholar] [CrossRef]
  90. Brett, R.A. The Ecology and Behaviour of the Naked Mole-Rat, Heterocephalus glaber Rüppell (Rodenti: Bathyergidae). Ph.D. Thesis, University College London (University of London), London, UK, 1986. [Google Scholar]
  91. Smith, M.; Buffenstein, R. Managed Care of Naked Mole-Rats. Adv. Exp. Med. Biol. 2021, 1319, 381–407. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Naked mole-rat colonies huddling together in the nest, with the breeding female lying on the top of the pile. Picture (A) shows a small naked mole-rat colony (six non-breeders + queen + her pups) with enough space in the nestbox for the newborn pups to easily crawl around the huddle. Picture (B) shows a large colony (60 mole-rats) with the breeding female lactating her pups on the top of the crowded huddle.
Figure 1. Naked mole-rat colonies huddling together in the nest, with the breeding female lying on the top of the pile. Picture (A) shows a small naked mole-rat colony (six non-breeders + queen + her pups) with enough space in the nestbox for the newborn pups to easily crawl around the huddle. Picture (B) shows a large colony (60 mole-rats) with the breeding female lactating her pups on the top of the crowded huddle.
Animals 13 00630 g001
Figure 2. (A) Pup survival status at day 10 for 869 pups from 79 litters born to 16 queens in 14 colonies. Each row of squares represents one litter, with each square representing a pup. The color of the square refers to the survival status of the pup at day 10 as indicated in the legend. Litters grouped together correspond to the different litters produced by a given queen (e.g., C, D, E…). Asterisks label the first (*) and second queen (**) of a colony. (B) Survival curve of male and female naked mole-rat pups until 60 days postpartum. The curves represent the proportion of male (black curve) and female (grey curve) pups born that remained alive as time passed (i.e., Kaplan–Meier estimator). The day 10, used to determine the survival status of pups in the survival analyses, is indicated by a vertical, dashed line.
Figure 2. (A) Pup survival status at day 10 for 869 pups from 79 litters born to 16 queens in 14 colonies. Each row of squares represents one litter, with each square representing a pup. The color of the square refers to the survival status of the pup at day 10 as indicated in the legend. Litters grouped together correspond to the different litters produced by a given queen (e.g., C, D, E…). Asterisks label the first (*) and second queen (**) of a colony. (B) Survival curve of male and female naked mole-rat pups until 60 days postpartum. The curves represent the proportion of male (black curve) and female (grey curve) pups born that remained alive as time passed (i.e., Kaplan–Meier estimator). The day 10, used to determine the survival status of pups in the survival analyses, is indicated by a vertical, dashed line.
Animals 13 00630 g002
Figure 3. Predictions (black line) and 95.0% confidence intervals (black dashed lines) displaying the probability for naked mole-rat pups to have lived for longer than 10 days, i.e., to survive from birth (day 0) to day 9, as a function of: (A) maternal number of mammae, (B) body mass of queens at conception, (C) pup birth body mass and (D) colony size (i.e., number of all colony members) with the grey dashed line representing a 50.0% probability of survival on pup survival.
Figure 3. Predictions (black line) and 95.0% confidence intervals (black dashed lines) displaying the probability for naked mole-rat pups to have lived for longer than 10 days, i.e., to survive from birth (day 0) to day 9, as a function of: (A) maternal number of mammae, (B) body mass of queens at conception, (C) pup birth body mass and (D) colony size (i.e., number of all colony members) with the grey dashed line representing a 50.0% probability of survival on pup survival.
Animals 13 00630 g003
Table 2. Summary table of the likelihood ratio tests used to compare offspring survival with respect to different predictors (n = 585 pups).
Table 2. Summary table of the likelihood ratio tests used to compare offspring survival with respect to different predictors (n = 585 pups).
Predictor VariablesBetaSEχ2dfp-Value
Individual characteristics 88.0220.001
Pup birth body mass6.550.9088.0110.001
Sex (male)−0.220.330.5310.453
Maternal characteristics 16.6620.015
Queen body mass−0.160.067.9610.015
Queen mammae2.250.876.3110.048
Social characteristics 6.1720.090
Colony size−0.070.036.1710.039
Litter size0.090.090.9110.423
Environmental characteristics 4.3130.316
Temperature0.260.580.2510.653
Humidity0.040.040.9410.404
Nestbox change−0.320.232.5310.182
Note: Results presented as unstandardized parameter estimates (beta), the standard error for the unstandardized parameter estimates (SE), the chi-squared test statistic (χ2), the degrees of freedom (df) and the p-value (p).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wetzel, M.; Courtiol, A.; Hofer, H.; Holtze, S.; Hildebrandt, T.B. Pup Recruitment in a Eusocial Mammal—Which Factors Influence Early Pup Survival in Naked Mole-Rats? Animals 2023, 13, 630. https://doi.org/10.3390/ani13040630

AMA Style

Wetzel M, Courtiol A, Hofer H, Holtze S, Hildebrandt TB. Pup Recruitment in a Eusocial Mammal—Which Factors Influence Early Pup Survival in Naked Mole-Rats? Animals. 2023; 13(4):630. https://doi.org/10.3390/ani13040630

Chicago/Turabian Style

Wetzel, Michaela, Alexandre Courtiol, Heribert Hofer, Susanne Holtze, and Thomas B. Hildebrandt. 2023. "Pup Recruitment in a Eusocial Mammal—Which Factors Influence Early Pup Survival in Naked Mole-Rats?" Animals 13, no. 4: 630. https://doi.org/10.3390/ani13040630

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop