High laboratory mouse pre-weaning mortality associated with litter overlap, advanced dam age, small and large litters

Abstract

High and variable pre-weaning mortality is a persistent problem in laboratory mouse breeding. Assuming a modest 15% mortality rate across mouse strains, means that approximately 1 million more pups are produced yearly in the EU to compensate for those which die. This paper presents the first large study under practical husbandry conditions to determine the risk factors associated with mouse pre-weaning mortality. We analysed historical records from 219,975 pups from two breeding facilities, collected as part of their management routine and including information on number of pups born and weaned per litter, parents’ age and identification, and dates of birth and death of all animals. Pups were counted once in their first week of life and at weaning, and once every one or two weeks, depending on the need for cage cleaning. Dead pups were recorded as soon as these were found during the daily cage screening (without opening the cage). It was hypothesized that litter overlap (i.e. the presence of older siblings in the cage when new pups are born), a recurrent social configuration in trio-housed mice, is associated with increased newborn mortality, along with advanced dam age, large litter size, and a high number and age of older siblings in the cage. The estimated probability of pup death was two to seven percentage points higher in cages with litter overlap compared to those without. Litter overlap was associated with an increase in death of the entire litter of five and six percentage points, which represent an increase of 19% and 103% compared to non-overlapped litters in the two breeding facilities, respectively. Increased number and age of older siblings, advanced dam age, small litter size (less than four pups born) and large litter size (over 11 pups born) were associated with increased probability of pup death.

About the Authors:

Gabriela M. Morello

Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Supervision, Writing – original draft

* E-mail: [email protected]

Affiliations Laboratory Animal Science, IBMC-Institute of Molecular and Cellular Biology, University of Porto, Porto, Portugal, i3S –Institute for Investigation and Innovation in Health, University of Porto, Porto, Portugal

ORCID logo http://orcid.org/0000-0002-5439-969X

Jan Hultgren

Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Writing – review & editing

Affiliation: Department of Animal Environment and Health, Swedish University of Agricultural Sciences, Skara, Sweden

Sara Capas-Peneda

Roles Data curation, Writing – review & editing

Affiliations Laboratory Animal Science, IBMC-Institute of Molecular and Cellular Biology, University of Porto, Porto, Portugal, i3S –Institute for Investigation and Innovation in Health, University of Porto, Porto, Portugal

ORCID logo http://orcid.org/0000-0002-0167-8812

Marc Wiltshire

Roles Data curation, Resources, Writing – review & editing

Affiliation: Babraham Institute, Babraham, Cambridge, United Kingdom

Aurelie Thomas

Roles Data curation, Resources, Writing – review & editing

Affiliation: Wellcome Sanger Institute, Hinxton, Cambridgeshire, United Kingdom

Hannah Wardle-Jones

Roles Data curation, Resources, Writing – review & editing

Affiliation: Wellcome Sanger Institute, Hinxton, Cambridgeshire, United Kingdom

Sophie Brajon

Roles Investigation, Writing – review & editing

Affiliations Laboratory Animal Science, IBMC-Institute of Molecular and Cellular Biology, University of Porto, Porto, Portugal, i3S –Institute for Investigation and Innovation in Health, University of Porto, Porto, Portugal

Colin Gilbert

Roles Conceptualization, Data curation, Funding acquisition, Investigation, Resources, Writing – review & editing

Affiliation: Babraham Institute, Babraham, Cambridge, United Kingdom

ORCID logo http://orcid.org/0000-0003-2126-5828

I. Anna S. Olsson

Roles Conceptualization, Funding acquisition, Investigation, Project administration, Resources, Writing – review & editing

Affiliations Laboratory Animal Science, IBMC-Institute of Molecular and Cellular Biology, University of Porto, Porto, Portugal, i3S –Institute for Investigation and Innovation in Health, University of Porto, Porto, Portugal

Introduction

High pre-weaning mortality of laboratory mice is a major welfare and economic problem affecting mouse breeding at academic and industrial laboratories worldwide. Previous studies report pup mortalities from less than 10% [1, 2] to as high as 49% [3] for C57BL/6 mice, one of the most commonly used mouse strains. Despite the general ongoing effort to reduce the number of animals in research and improve their welfare according to the 3R principle for research [4], high pre-weaning mortality rates persist and very little systematic research has been done to identify causes of poor survival. Data from experimental and observational studies conducted by the authors of this work at different breeding facilities in three different countries revealed that 32% of 344 litters (retrospective analysis, Germany [5]), 33% of 55 litters (experimental data, U.K. [6]), and 18% of 510 litters (experimental data, Portugal [7]) were completely lost, with the overall mortality varying from 25% [7] to 52% in trio-bred mice [6] in the experimental studies. If a modest level of 15% mortality is assumed across all mouse strains, at least 1 million more mice must be produced every year just in the European Union (EU) to compensate for pups that die before they can be used in science (estimate based on the number of mice used yearly in research in the EU; European Commission 2020 [8]). Such losses are contrary to the 3R principle that is now explicit in EU legislation [9] and incur extra breeding costs of €5–8 million yearly. Several environmental, management and behavioural factors have been linked to pup mortality, such as thermal environment of the cage, level of parental care, dam age, litter size, provision of nest material, and cage manipulation [2, 6, 10–14]. Recently, we identified the presence of older litter mates in the cage when a new litter is born (litter overlap) as a major factor affecting pup survival [6]. In a study with 55 litters of C57BL/6 mice (n = 521 pups) housed in trios [6], a 2.3 fold increase was found in litter loss (death of the entire litter) in cages where older littermates were present, compared to trio cages with no older littermates. Litter overlap happens in both trio (two adult females and one male) and pair (one adult female and one male) housing, which are the most common configurations in mouse breeding. Although litter overlap is more frequent in trios due to the presence of two breeding females, trios wean more pups per cage [15] and findings are inconsistent as to whether the number of pups weaned per female is similar [15] or reduced [16] in trios compared to pairs. One possible reason for this is that litter overlap in pair cages affects pup mortality more severely as compared to trio cages. In pair cages, litter overlap occurs when the only female of the cage gives birth before weaning her previous litter. In these cases, the age gap between litters becomes large, which might be especially detrimental to pup survival.

Previous research into factors affecting laboratory mouse reproduction used primarily experimental study approaches, where the sample size was small and animal management and data collection differed from standard practice in a breeding facility. With the increasing use of breeding management software, it is now possible to use much larger datasets representing the reality of practical laboratory mouse breeding. In this study, a dataset of 219,975 pups was analysed from two different collaborating breeding facilities in the UK (58,692 and 161,283 pups), by modelling the risk of a newborn mouse dying as a function of litter overlap, age and number of older littermates, number of pups born, as well as of age of the dam. It was hypothesized that litter overlap is a recurrent social configuration and that the risk of pup mortality increases with litter overlap, advanced dam age, large litter size, as well as a high number and age of older siblings in the cage.

Material and methods

Data retrieval

Historical mouse breeding data was provided by two collaborating facilities. Therefore, this study did not involve any type of animal manipulation, observation, or use. Mouse breeding in the collaborating facilities was performed in line with the UK Animals (Scientific Procedures) Act of 1986.

Mouse breeding data were made available by two collaborating breeding facilities (C1, the Babraham Institute and C2, the Wellcome Sanger Institute). Historical production data were downloaded directly from their breeding management software (MCMS, Mouse Colony Management System, the Wellcome Sanger Institute Data Centre). C1 provided data from January 2014 to October 2018, and C2 provided data from January 2010 to March 2019. The datasets contained information on litter identity, breeding adults’ identities, date of birth, date of death, number of pups born, and number of pups weaned.

Animals, housing, and management

The studied dataset contained a total of 34,949 C57BL/6 litters and 219,975 pups. All mice were housed in trios (two females and one male) in individually ventilated cages (IVC). Details on animals, housing, and management are shown in Table 1.

[Figure omitted. See PDF.]

Table 1. Animal, housing, and management characteristics for the collaborating animal facilities.

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

Statistical analysis

The historical mouse breeding data that was analyzed represent 58,692 pups in 9,261 litters and 161,283 pups in 25,688 litters from C1 and C2, respectively. Required information was retrieved using Scilab (version 6.0.1, Scilab Enterprises, Rungis, France), resulting in one data line per pup. A total of 11% (C1) and 21% (C2) of the data originally provided was excluded (see S1 Fig), mainly due to incongruent data records, implausibly large litters (more than 13 pups, unless confirmed as correct), unreliable information on number and age of older pups in the cage, or missing information. Male pups at C2 euthanized at day 7 post-partum or later were coded as surviving. Litters with males euthanized before day 7 post-partum were excluded.

Pup mortality before weaning was coded as 0 (survived) or 1 (died) and used as the dependent variable. Environmental and social factors were considered as risk factors for pup death. Independent variables representing environmental factors considered for analysis were Collaborator (C1 or C2), Season (Winter, Spring, Summer, Fall), Month (as an alternative predictor to Season), Weekday, and Year, while independent variables for social risk factors included Dam Age (continuous), Father Age (continuous), Litter Size (number of pups born; continuous), litter Overlap (whether or not older siblings were present at the time of birth of the focus litter; no or yes), Sibling Number (number of older pups in the cage at the time of birth of the focus litter; continuous) and Sibling Age (age of the older siblings; continuous).

The risk of pup death was modelled by mixed logistic regression, using the GLIMMIX procedure in SAS (2018 University Edition, SAS Institute Inc., Cary, NC, USA). Multicollinearity among independent variables was checked by using the variance inflation factor (VIF) and regressing each independent variable on the others. As a consequence, Year, Month, and Father Age were excluded from the analysis.

Data of C1 and C2 were combined into one dataset for analysis. Two separate analyses were run: one on the complete dataset (Model 1) and one on a subset of the data containing only pups born in situations with litter Overlap (Model 2) (see S1 Fig). With the first model (Model 1), we aimed at comparing overlapped with non-overlapped litters, accounting for the effects of the remaining covariates and confounders. With the second model (Model 2), we further evaluated how pup probability of death varied with the variation in the age and number of older siblings, also considering the effects of the remaining covariates and confounders. In both models, litter identity was included as a random effect to account for clustering. The models were built by adding one independent variable of interest (Dam Age, Litter Size, Overlap (Model 1) or Sibling Number and Sibling Age (Model 2) at a time in a stepwise process with bidirectional elimination. Independent variables with p ≤ 0.05 were kept in the model. Weekday, Season, and Collaborator were then tested one at a time as confounders, followed by possible interactions and higher order terms. Pair-wise t-tests were performed to compare least-squares means of Overlap, Weekday, and Season, applying a Tukey-Kramer adjustment for multiple pairwise comparisons.

Results

The percentage of pups dying before weaning was 39% at C1 and 14% at C2, while the mean Litter Size was 7.6 pups born/litter in both collaborators. In 42% of the C1 litters and 78% of the C2 litters no pups died. The percentage of litters with at least 90% death rate was 28% at C1 and 9% at C2. Approximately 50% and 57% of the litters were born with the presence of older siblings in the cage (litter Overlap) in C1 and C2, respectively.

Model details are available in S1 and S2 Tables. In both models, pup probability of death was affected by Collaborator, Season, Weekday, Dam Age, and Litter Size in a quadratic fashion. Additionally, pup probability of death was affected by the variables of interest Overlap in Model 1 (all pups), as well as Sibling Number, Sibling Age, and their interaction in Model 2 (pups born in overlapped litters). Models’ outcomes are presented separately and by Collaborator, due to existent interactions among this and the remaining variables, as described below.

Outcomes from Model 1: Overlap vs. non-overlapped litters

The estimated probability of pup death was seven (C1) and two (C2) percentage points higher (p < 0.001) in cages with the presence of older siblings compared to cages without an older litter (Fig 1A and 1B). At C1, 31% of the overlapped and 26% of the non-overlapped litters had a total litter loss (all pups dying, Fig 1C), whereas at C2, the corresponding figures were 12% and 6% (Fig 1D). Mean Litter Size at birth was 6.3 ± 2.9 pups with no significant differences between collaborators (p = 0.087).

[Figure omitted. See PDF.]

Fig 1. Probability of pup death and litter mortality distribution.

Probability of pup death in litters without (NO) or with (YES) the presence of older siblings in the cage (litter Overlap) at (A) C1, the Babraham Institute, and (B) C2, the Wellcome Sanger Institute, based on least-square means. Percentage of litters by category of pre-weaning mortality in Overlapped and Non-Overlapped litters at (C) C1 and (D) C2, based on raw data. Numbers within brackets in the x-axis designate lower (left side) and upper (right side) endpoints of mortality range. An open bracket next to a number designates a non-inclusive endpoint.

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

Although Weekday and Season were not added to the models as variables of interest, these factors turned out to be confounders to the models, while interacting with Collaborator (p < 0.001, p = 0.010, respectively). The probability of death consistently decreased towards the end of the week at both collaborators, while the effect of Season lacked a consistent pattern between collaborators (available in S2 Fig). C1 was able to provide records on cage cleaning dates per Weekday for the period of 17 months (April 2018 to November 2019) for the colony. Cage change events also peaked in the beginning of the week and were lower towards the end of the week (available in S3 Fig).

Outcomes from Model 2: Effects of Sibling Number, Sibling Age, Dam Age, and Litter Size

The predicted probability for a pup to die as a function of Sibling Number and Age, Dam Age, and Litter Size, is illustrated by Figs 2 and 3. Increased Sibling Number, Sibling Age, and Dam Age were associated (p < 0.001) with an increase in the probability of pups dying at both C1 and C2 (Figs 2 and 3).

[Figure omitted. See PDF.]

Fig 2. Probability of pup death by Dam Age and Sibling Age.

Predicted probabilities (least-squares means) of a pup to die as a function of Dam Age for three distinct levels of Litter Size (number of pups born), at C1, the Babraham Institute and C2, the Wellcome Sanger Institute. Each line corresponds to predictions for a specific value of Sibling Age, as depicted in the legend next to the top left graph. Predictions were obtained while assuming six older pups in the cage, the most recurrent Weekday (Thursday) and the most common Season (Spring) in the combined dataset.

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

[Figure omitted. See PDF.]

Fig 3. Probability of pup death by Sibling Number and Sibling Age.

Predicted probabilities (least-squares means) of a pup to die as a function of Sibling Number (number of older pups in the cage at the birth of the focal litter) for three distinct levels of Litter Size (number of pups born), at C1, the Babraham Institute and C2, the Wellcome Sanger Institute. Each line corresponds to predictions for a specific value of Sibling Age, as depicted in the legend next to the top left graph. Predictions were obtained while assuming a mean Dam Age of 130d, the most recurrent Weekday (Thursday) and the most common Season (Spring) in the combined dataset.

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

Mean litter size at birth for overlapped data was slightly larger (p < 0.001) in C1 (6.1 ± 2.7 pups) compared to C2 (5.8 ± 2.9 pups). Pup death probability seem to be higher in either small or large litters due to the quadratic nature of this relationship (see S2 Table). Prediction estimates (least-square means) indicate that, generally, the risk of pup death increases as Litter Size decreases below four and increases above 11 pups per litter, for all combinations of Sibling Age, Sibling Number and Dam Age.

Discussion

Neonatal mortality is a large problem in laboratory mouse breeding and improving pup survival is a key to improve the efficiency and sustainability of producing laboratory mice while complying with the 3R principle [9]. Here, we present the first large scale study of management, environmental and animal factors affecting pup survival based on data from over 200,000 C57BL/6 pups born during a period of 5 to 10 years in two large mouse breeding facilities. The study indicates that litter overlap, a social configuration which frequently occurs in trio-breeding cages, results in a 30% to 60% increase in the probability of neonatal pup death (considering 7 and 2 percentage points increase from mean pup probability of deaths in non-overlapped litters of 21% and 3% in C1 and C2, respectively, Fig 1A and 1B). The higher the number and age of older siblings in the cage, the greater is the risk of neonatal pup death. The probability of pups to die in the presence of older siblings was also affected by dam age and litter size.

Pup probability of death

Irrespective of any mortality risks identified, the average litter mortality rates obtained both at C1 (39%) and C2 (14%) were higher than what was previously reported for C57BL/6J mice (8% pre-weaning mortality [17] with pup counting at weaning and 3% mortality at three days post-partum [18] with pup number obtained from video-records). Litter mortality was higher in C1 compared to C2 and no differences were found between collaborators in the number of pups born per litter for the whole dataset. Husbandry differences between C1 and C2 may contribute to the mortality difference between both institutes as cage temperature, as well as nest and bedding amount and quality affect rodents’ breeding performance [11–14]. An additional factor may be the practice of euthanizing male pups in C2. Euthanized males after seven days of age were considered as survivals, but this is an assumption and may have led to an underestimation of mortality. Finally, any differences in accuracy of data entry may have affected the results. C1 had a more consistent and early counting of pups than C2. Thus it is possible that C1 has a better accuracy in detecting the number of pups born compared to C2 where pups born could be underestimated, considering that most of deaths happen within 48h post-partum and dead pups often get cannibalized by the dam, thus not seen by the caretakers. Productivity differences in pup survivability due to the distinct mouse sub-strains between C1 and C2 could also underlie the differences in overall mortality found between C1 and C2.

We found higher mouse pre-weaning risk of death in trios with overlapped litters, i.e. litters born when an older litter was present, compared to non-overlapped litters, accounting for the variation in litter size, dam age, and confounders (weekday and season of birth), which is in agreement with previous experimental findings. In outbred mice derived from the C57BL/6J, BALB/c, and DBA/1J strains, Schmidt et al. [19] found pup mortality to increase with increasing age gap between the litters sharing a cage at a specific time and attributed this to infanticide. Understanding why being born into a cage with an older litter is so dangerous requires information about events around pup death and the condition of dead and dying pups. Pup mortality is often associated with infanticide [16, 20], assuming that cannibalized pups were killed before they were eaten. However, our previous behavior studies suggest that infanticide precedes less than 15% of the cannibalism events [6, 21] and that pups die primarily from other causes than direct killing. Litter asynchrony, which often leads to overlap, is likely to increase unequal competition for access to milk and parental care, trauma caused by trampling and stepping of newborns by the adults or the older siblings, and problems related with increased cage stocking density.

Early access to milk is essential for the survival of newborn pups. Measurements of pup energy losses due to metabolism between nursing bouts, extrapolated for a period of 24 hours, revealed that if pups did not receive milk during their first day of life, they would lose approximately 8% of their birth weight [22], which would likely reduce their chances of survival. The presence of older and consequently heavier, more developed and more mobile pups in the cage may have interfered with the access of the newborns to milk in general and specifically to steal the iron-rich milk that is present in higher concentrations during the first week of lactation [22]. Also, older and heavier siblings may be able to displace light newborn pups more easily from the dam's nipples, as compared to younger siblings. This could partly explain the interaction found in this work between the number and the age of older siblings in the cage, affecting pup death probability. Our results showed that the smaller the age gap between the newborn and the older litter, the lower was the pup probability of death. Thus, older pups which are less than a week older than the younger siblings may not be able to interfere with the youngest access to appropriate amount and quality of nutrition.

The presence of two litters in the cage has been demonstrated to increase parturition duration and affect parental behavior. Adults in trio cages with two litters were observed to care for their newborn pups a total of 20% less time (all the three adults together) than adults with one single litter in the cage [6], while parental investment is known to improve the chances of survival of young mice. In fact, C57BL/6 females which lost their litters entirely have been found to spend more time outside the nest and invested less time in building the nest prior to parturition [23], while the presence of males in cages with breeding females (CD-1) has been demonstrated to increase pup survival by facilitating maternal behavior [24]. Thus, reduced parental care in cages with more than one litter can be one of the mechanisms through which pup survivability is reduced in the presence of an older litter.

Most often, when there are two females sharing a cage, they also share the same nest, and younger lighter pups get clustered together with the older, heavier, and more mobile pups. Data on post-mortem inspection performed in 324 C57BL/6J pups found dead, by the authors of this study, revealed that 24% of the pups had some kind of traumatic lesion, including bite wounds and bruises [25].

Alloparenting, meaning that an individual contributes to parenting young which are not their offspring, is known to take place in both house mice [26] and laboratory mice [6, 16]. House mouse females are able to switch strategy during their reproductive lifetime: Ferrari et al. [27] found that half of the females switched and that young and inexperienced females were more prone to join another female, while older and more experienced females (and, presumably, with a larger body mass and more capable of producing larger amounts of milk) tended to raise their litters solitarily.

The reproductive consequences of communal breeding are complex. In their study of wild house mice populations, Ferrari et al. [27] found an increase in pup mortality (between days 1–3 and 13 days old of the older litter), aggravated with each additional litter that was found together with the focal litter, but another study of the same population found an increase in pup survival in litters reared communally [28]. In contrast, comparing pairs (where alloparenting is impossible) and trios (where alloparenting is common), Garner et al. [16] found no difference in mortality. However, Garner et al. [16] found number of pups weaned per female to be reduced in trios, whereas Wasson [15] found no difference between pairs and trios. The reduction in breeding performance found by Garner et al. [16] was attributed to a possible reproduction suppression in trios before birth. In the present work, however, the differences in breeding performance between overlapped and non-overlapped litters are likely not due to reproductive suppression before birth, because our models accounted for the variation in number of pups born per litter.

Higher stocking density leads to increased humidity and gas concentration in the air, with ammonia levels increasing in 30 ppm with litter growth from birth to weaning [29] and reaching over 100 ppm in individually ventilated cages with trios 4d post-partum with a maximum concentration of over 600 ppm at weaning [30]. The presence of a second litter in the cage substantially increased ammonia levels, which indicated an effect of housing density on ammonia concentration [30]. Whereas the impact of these ammonia levels on newborn mice has not been studied, ammonia levels from 25 to over 600 ppm have been demonstrated to increase pathological scores, destroy the surface layers of the trachea epithelium lining and increase the severity of rhinitis, otitis, tracheitis, and pneumonia in rats and mice [29–32]. The gas concentration problem may be aggravated by the fact that animal care-takers generally tend to avoid cleaning cages when litters were just born (to avoid pup disturbance).

Dam Age

Pup death probability increased as Dam Age increased in both collaborators. In this study, Dam Age and parity were confounded. Thus, it was not possible to distinguish effects on pup death probability of the dam's age and its birthing experience. Decreased productivity and increased mortality in first-parity litters have been reported for a few different species [33, 34], but for mice this subject remains controversial. While first-parity BALB/c and 129/Sv dams were reported to wean fewer pups per litter compared to later parity ones [35], we previously found an increase in pup survival with lower parities [6] in an experimental study conducted in C1, whereas another study did not find any significant differences in pup loss between first- and later parity C57BL/6 or BALB/c dams [5].

The results for Dam Age are in agreement with those from Tarín et al. [10], who found increased pup mortality and incidence of litters with at least one cannibalized pup with increased parity. Tarín et al. [10] also compared breeding performance between dams (F1 of C57BL/6JIco × CBA/JIco) who started their reproductive life at age 70d (young) and 357d (old). The authors found no differences between dam age group on pre-weaning mortality and litter size both at birth and at weaning, but reported that young dams produced F2 litters with higher expectation of survival and body weight than those of old dams.

The effect of maternal age on mortality may be related to the age-related hearing loss that is common to females of the strain studied [36]. Pup vocalizations are important signals of cold stress and distress situations and affect maternal behavior such as pup retrieval and maternal care [37, 38]. Thus, if older C57BL/6 females are not able to hear and communicate with their pups appropriately, pup chances of survival could be reduced through a reduction in maternal care, especially knowing how important maternal care is for pup survival and how delicate pups are when they are born (as reviewed by Latham and Mason [39]).

Number of pups born

The reduced survivability in small litters is in agreement with previous reports for C57BL/6 [6] and F1 hybrid (C57BL/6JIco × CBA/JIco) mice [10]. One possible explanation for this is that a reduced number of pups may reduce the huddle formed by the pups in nest, which is an important mechanism through which newborn rodents maintain their body heat [40].

Litter size may also affect parental care, which is essential for pup survival. Ehret and Bernecker [41] demonstrated that early pup vocalization, which gradually increases in frequency after birth, is essential to maintain maternal attention at high levels, which leads to improved pup weight gain, as compared to pups from dams which were unable to hear them. Therefore, it is possible that a small newly born litter does not emit sufficient vocal cues to ensure sufficient maternal care. Rat litters [42] with one single pup were found to perform only about 10% and 5% of the suckling stimuli performed by litters of 10 and 22 pups. As a consequence, the milk yield of dams (estimated based on adjusted measures of the pups' daily weight gain) raising single pups was only -0.4% to 7.0% of those raising 10 pups, which led one-pup litters to have the lowest growth rate. More than half of the one-pup litters did not show any weight gain in the first five days post-partum. From an evolutionary perspective, a small litter is less worth investing in than a larger litter: Maestripieri and Alleva [43] demonstrated that CD-1 dams of large litters (eight pups) spent more than twice as much time displaying litter defense behaviors against intruder males than dams of small litters (four pups). The increase in pup death probability found in litters of 12 pups and above, on the other hand, may be a result of increased sibling competition for access to milk, as discussed above, and also may represent a ceiling in milk production capacity by the dams [22, 44, 45].

Weekday and Season

In both collaborators, there seemed to be a decrease the probability of pup death towards the end of the week, possibly associated with the timing of cage changes, a management routine which affects the mice as well as the accuracy of mortality detection. In C1, which provided records on cage cleaning dates, these closely mirrored the pattern of pup death probability. To reliably count the number of pups, the cage must be opened and animals moved, something that often only happens at cage cleaning when manipulation is unavoidable. Mortality is, therefore, likely to be more accurately detected for litters born on cage changing days. For example, an eight-pup litter born on a Tuesday with cage cleaning schedule for the same day will be recorded as an eight-pup litter. If two of these pups die in the following 24 hours, this litter's pre-weaning mortality will be recorded as being 25% at weaning. A similar litter born on a Saturday with two pups dying on Sunday, and subsequently cannibalized, will be recorded at the Tuesday cage change as a litter with six pups born with no pre-weaning deaths.

Still, the mouse disturbance hypothesis cannot be disregarded. If pup mortality is affected by cage change, the same pattern would be expected in cage change frequency as in pup death probability per Weekday (of birth). Reeb-Whitaker et al. [1] found a higher pup mortality in cages with weekly changes than those changed once every two weeks. Cage change requires that mice are moved from the dirty to a clean cage, an event that triggers a stress response evidenced by increases in serum corticosterone [46] and general activity [47]. It is possible, therefore, that cage change interferes with parental behavior in breeding cages, which could aggravate pup mortality around those days.

Conclusions

The present study revealed that high pre-weaning mortality in laboratory mice (C57BL/6) is associated with advanced dam age, litter overlap, the presence of a high number and age of older siblings in the cage, and a small (less than four) or large (more than 11 pups) litter. The dynamics of parental care, sibling competition for access to milk, and issues related with the number of animals in the cage may underlie the effects found in pup mortality caused by the identified risks. Future studies should address sibling competition and parental behavior in asynchronized litters.

Supporting information

[Figure omitted. See PDF.]

S1 Fig. Data processing and modeling.

Schematic illustration of the data analysis process performed on the data provided by each of the two collaborators. (A) After original datasets were processed with the help of a data-cleaning algorithm, the data were split into two datasets: One with overlapped and non-overlapped litters and another with only overlapped litters. (B) Model 1 was performed to evaluate which environmental and social factors were a risk for pup mortality, while Model 2 was performed to provide in more details how the social factors affect pup mortality in overlapped litters. Model outcomes are listed in (C).

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

(PDF)

S2 Fig. Probability of pup death by Weekday and Season.

Predicted Probability (least-square means) of a pup to die as a function of (A) Weekday and (B) Season of birth at C1, the Babraham Institute, and C2, the Wellcome Sanger Institute. Data points with distinct labeled letters indicate statistical difference at 95% confidence level.

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

(TIF)

S3 Fig. Cage change frequency and probability of pup death by Weekday.

Cage change frequency and predicted probability (least-square means) of a pup to die as a function of Weekday at C1, the Babraham Institute. Data points with distinct labeled letters indicate statistical difference at 95% confidence level. Cage change frequency is depicted as the percentage per weekday of the 78 cage change episodes which happened from April 2018 to November 2019 (available data records), in the studied room of C1.

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

(TIF)

S1 Table. Solutions for fixed effects of the final model predicting the odds of pup death fitted in the whole processed dataset.

n.a. = not applicable. aVariable Litter Size was centered by its mean.

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

(PDF)

S2 Table. Solutions for fixed effects of the final model predicting the odds of pup death fitted in the dataset containing only overlapped litters.

n.a. = not applicable. aVariable Litter Size was centered by its mean.

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

(PDF)

Acknowledgments

Dr. Luiz Henrique Antunes Rodrigues for providing technical input for the statistical modeling process. The Babraham Institute and the Wellcome Sanger Institute, UK, for providing the historical data used in this manuscript. Dr. Sofia Lamas for providing technical and scientific input.

Citation: Morello GM, Hultgren J, Capas-Peneda S, Wiltshire M, Thomas A, Wardle-Jones H, et al. (2020) High laboratory mouse pre-weaning mortality associated with litter overlap, advanced dam age, small and large litters. PLoS ONE 15(8): e0236290. https://doi.org/10.1371/journal.pone.0236290

References

1. Reeb-Whitaker CK, Paigen B, Beamer WG, Bronson RT, Churshill GA, Schweitzer IB, et al. The impact of reduced frequency of cage changes on the health of mice housed in ventilated cages. Lab Anim. 2001;35:58–73. pmid:11201289

2. Gaskill BN, Pritchett-Corning KR, Gordon CJ, Pajor EA, Lucas JR, Davis JK, et al. Energy reallocation to breeding performance through improved nest building in laboratory mice. PLoS One. 2013;8(9):e74153. pmid:24040193

3. Inglis CA, Campbell ER, Auciello SL, Sarawar SR. Effects of Enrichment Devices on Stress-Related Problems in Mouse Breeding. Johns Hopkins Cent Altern to Anim Test Final Rep Anim Welf Enhanc Award 2004. 2004;(1):1–9.

4. Russell WMS, Burch RL. The principles of humane experimental technique. Methuen; 1959.

5. Weber EM, Algers B, Wurbel H, Hultgren J, Olsson IAS. Influence of Strain and Parity on the Risk of Litter Loss in Laboratory Mice. Reprod Dom Anim. 2013;48:292–6.

6. Brajon S, Morello GM, Teixeira MS, Hultgren J, Gilbert C, Olsson IAS. Social environment as a cause of litter loss in laboratory mouse: A behavioural study. Appl Anim Behav Sci. 2019;218:104827.

7. Morello GM, Brajon S, Capas-Peneda S, Hultgren J, Ferreira JM, Gilbert C, et al. Understanding pup mortality in laboratory mouse breeding: How the presence of an older litter in the cage aggravates pre-weaning mortality in mice housed in trios and in pairs. 2019 Jul 3 [cited 2020 February 23]. In UFAW International Animal Welfare Sc. Symp Available from https//www.ufaw.org.uk/downloads/bruges-2019—ufaw-symposium-programme-booklet-online-v1.pdf.

8. European Commission. 2019 report on the statistics on the use of animals for scientific purposes in the member states of the European Union in 2015–2017. 2020 Feb 5 [Cited 2020 Feb 23]. Brussels. Available from: https://eur-lex.europa.eu/legal-content/EN/.

9. The European Parliament and the Council of 22 September 2010. 2010 Sep 22 [Cited 2020 Feb 23]. Directive 2010/63/EU on the Protection of Animals Used for Scientific Purposes. Brussels. Available from: https://eur-lex.europa.eu/homepage.html.

10. Tarín JJ, Gómez-Piquer V, Rausell F, Navarro S, Hermenegildo C, Cano A. Delayed motherhood decreases life expectancy of mouse offspring. Biol Reprod. 2005;72(6):1336–43. pmid:15689534

11. Potgieter FJ, Wilke PI. Effect of different bedding materials on the reproductive performance of mice. J S Afr Vet Assoc. 1997;68(1):8–15. pmid:9186933

12. Eskola S, Kaliste-korhonen E. Nesting material and number of females per cage: effects on mouse productivity in BALB/c, C57BL/6J, DBA/2 and NIH/S mice. Lab Anim. 1999;33:122–8. pmid:10780814

13. Paul MJ, Tuthill C, Kauffman AS, Zucker I. Pelage insulation, litter size, and ambient temperature impact maternal energy intake and offspring development during lactation. Physiol Behav. 2010;100:128–34. pmid:20184907

14. Gaskill BN, Winnicker C, Garner JP, Pritchett-corning KR. The naked truth: Breeding performance in nude mice with and without nesting material. Appl Anim Behav Sci. 2013;143(2–4):110–6.

15. Wasson K. Retrospective analysis of reproductive performance of pair-bred compared with trio-bred mice. J Am Assoc Lab Anim Sci. 2017;56(2):190–3. pmid:28315650

16. Garner JP, Gaskill BN, Pritchett-Corning KR. Two of a kind or a full house? Reproductive suppression and alloparenting in laboratory mice. PLoS One. 2016;11(5):e0154966. pmid:27148872

17. Mouse phenome database at the Jackson Laboratory. Phenotype measure: Jax5 wean_born_ration. [Cited 2020 Feb 23]. Bar Harbor. Available from: https://phenome.jax.org/measures/31406.

18. Murray SA, Morgan JL, Kane C, Sharma Y, Heffner CS, Lake J, et al. Mouse gestation length is genetically determined. PLoS One. 2010;5(8):e12418. pmid:20811634

19. Schmidt J, Kosztolányi A, Tökölyi J, Hugyecz B, Illés I, Király R, et al. Reproductive asynchrony and infanticide in house mice breeding communally. Anim Behav. 2015;101:201–11.

20. Schmidt J, Kosztolányi A, Tökölyi J, Hugyecz B, Illés I, Király R, et al. Reproductive asynchrony and infanticide in house mice breeding communally. Anim Behav. 2015;101:201–11.

21. Weber EM, Algers B, Hultgren J, Olsson IAS. Pup mortality in laboratory mice—infanticide or not? Acta Vet Scand. 2013;55(1):83.

22. König B, Riester J, Markl H. Maternal care in house mice (Mus musculus): II. The energy cost of lactation as a function of litter size. J Zool. 1988;216:195–210.

23. Weber EM, Hultgren J, Algers B, Olsson IAS. Do laboratory mouse females that lose their litters behave differently around parturition? PLoS One. 2016;11(8):e0161238. pmid:27575720

24. Wright SL, Brown RE. Maternal behavior, paternal behavior, and pup survival in CD-1 albino mice (Mus musculus) in three different housing conditions. J Comp Psychol. 2000;114(2):183–92. pmid:10890590

25. Capas-Peneda S, Morello GM, Lamas S, Olsson IAS, Gilbert C. Causes of death in newborn C57BL/6 mice. BioRxiv [Preprint]. 2020 [cited 2020 Feb 23]. Available from: https://doi.org/10.1101/2020.02.25.964551.

26. König B. Cooperative care of young in mammals. Naturwissenschaften. 1997;84(3):95–104. pmid:9112240

27. Ferrari M, Lindholm AK, König B. Fitness consequences of female alternative reproductive tactics in house mice (Mus musculus domesticus). Am Nat. 2019;193(1):106–24. pmid:30624110

28. Auclair Y, König B, Lindholm AK. Socially mediated polyandry: A new benefit of communal nesting in mammals. Behav Ecol. 2014;25(6):1467–73. pmid:25419087

29. Eveleigh JR. Murine cage density: Cage ammonia levels during the reproductive performance of an inbred strain and two outbred stocks of monogamous breeding pairs of mice. Lab Anim. 1993;27(2):156–60. pmid:8501897

30. Di Vincenti LJ, Moorman-White D, Bavlov N, Garner M, Wyatt J. Effects of housing density on nasal pathology of breeding mice housed in individually ventilated cages. Lab Anim (NY). 2012;41(3):68–76.

31. Vogelweid CM, Zapien KA, Honigford MJ, Li L, Li H, Marshall H. Effects of a 28-day cage-change interval on intracage ammonia levels, nasal histology, and perceived welfare of CD1 mice. J Am Assoc Lab Anim Sci. 2011;50(6):868–78. pmid:22330779

32. Broderson JR, Lindsey JR, Crawford JE. The role of environmental ammonia in respiratory mycoplasmosis of rats. Am J Pathol. 1976;85(1):115–30. pmid:970435

33. Gourdine JL, Bidanel JP, Noblet J, Renaudeau D. Effects of season and breed on the feeding behavior of multiparous lactating sows in a tropical humid climate. J Anim Sci. 2006;84:469–80. pmid:16424276

34. Dhakal K, Maltecca C, Cassady JP, Baloche G, Williams CM, Washburn SP. Calf birth weight, gestation length, calving ease, and neonatal calf mortality in Holstein, Jersey, and crossbred cows in a pasture system. J Dairy Sci. 2013;96(1):690–8. pmid:23084888

35. Whitaker JW, Moy SS, Pritchett-Corning KR, Fletcher CA. Effects of enrichment and litter parity on reproductive performance and behavior in BALB/c and 129/Sv mice. J Am Assoc Lab Anim Sci. 2016;55(4):387–99. pmid:27423144

36. Henry KR. Sex- and age-related elevation of cochlear nerve envelope response (CNER) and auditory brainstem response (ABR) thresholds in C57BL/6 mice. Hear Res. 2002;170(1–2):107–15. pmid:12208545

37. Ehret G, Bernecker C. Low-frequency sound communication by mouse pups (Mus musculus): Wriggling calls release maternal behaviour. Anim Behav. 1986;34(3):821–30.

38. Hahn ME, Lavooy MJ. A review of the methods of studies on infant ultrasound production and maternal retrieval in small rodents. Behav Genet. 2005;35(1):31–52. pmid:15674531

39. Latham N, Mason G. From house mouse to mouse house: The behavioural biology of free-living Mus musculus and its implications in the laboratory. Appl Anim Behav Sci. 2004;86:261–89.

40. Alberts JR. Huddling by rat pups: Group behavioral mechanisms of temperature regulation and energy conservation. J Comp Physiol Psychol. 1978;92(2):231–45. pmid:670452

41. Ehret G, Bernecker C. Low-frequency sound communication by mouse pups (Mus musculs): Wriggling calls release maternal behaviour. Anim Behav. 1986;34:821–30.

42. Russell JA. Milk yield, suckling behavior and milk ejection in the lactating rat nursing litters of different sizes. J Physiol. 1980;303(1):403–15.

43. Maestripieri D, Alleva E. Litter defence and parental investment allocation in house mice. Behav Processes. 1991;23(3):223–30. pmid:24923517

44. Johnson MS, Thomson SC, Speakman JR. Limits to sustained energy intake: I. Lactation in the laboratory mouse Mus musculus. J Exp Biol. 2001;204:1925–35. pmid:11441034

45. Vaanholt LM, Duah OA, Balduci S, Mitchell SE, Hambly C, Speakman JR. Limits to sustained energy intake. XXVII. Trade-offs between first and second litters in lactating mice support the ecological context hypothesis. J Exp Biol. 2018;221(5).

46. Rasmussen S, Miller MM, Filipski SB, Tolwani RJ. Cage change influences serum corticosterone and anxiety-like behaviors in the mouse. J Am Assoc Lab Anim Sci JAALAS. 2011;50(4):479–83. pmid:21838975

47. Castelhano-Carlos MJ, Baumans V. The impact of light, noise, cage cleaning and in-house transport on welfare and stress of laboratory rats. Lab Anim. 2016;43:311–27.