Skip to main content
  • Research article
  • Open access
  • Published:

Female behavior drives the formation of distinct social structures in C57BL/6J versus wild-derived outbred mice in field enclosures

A Correction to this article was published on 01 April 2024

This article has been updated

Abstract

Background

Social behavior and social organization have major influences on individual health and fitness. Yet, biomedical research focuses on studying a few genotypes under impoverished social conditions. Understanding how lab conditions have modified social organizations of model organisms, such as lab mice, relative to natural populations is a missing link between socioecology and biomedical science.

Results

Using a common garden design, we describe the formation of social structure in the well-studied laboratory mouse strain, C57BL/6J, in replicated mixed-sex populations over 10-day trials compared to control trials with wild-derived outbred house mice in outdoor field enclosures. We focus on three key features of mouse social systems: (i) territory establishment in males, (ii) female social relationships, and (iii) the social networks formed by the populations. Male territorial behaviors were similar but muted in C57 compared to wild-derived mice. Female C57 sharply differed from wild-derived females, showing little social bias toward cage mates and exploring substantially more of the enclosures compared to all other groups. Female behavior consistently generated denser social networks in C57 than in wild-derived mice.

Conclusions

C57 and wild-derived mice individually vary in their social and spatial behaviors which scale to shape overall social organization. The repeatable societies formed under field conditions highlights opportunities to experimentally study the interplay between society and individual biology using model organisms.

Background

Laboratory house mice are the premier model organism in biomedical research due to their small size, rapid breeding cycle, and the ready deployment of precise experimental manipulations using powerful genetic and neurobiological tools [1,2,3,4,5,6]. The wide availability of classical inbred mouse strains has allowed the scientific community to amass diverse physiological, genomic, neurobiological, and behavioral datasets on repeatable genotypes across labs and studies [2, 7,8,9,10,11,12,13,14,15]. While highly controlled conditions are required for many experiments, there is growing recognition that environmentally impoverished traditional laboratory approaches limit our ability to understand many complex biological processes [16,17,18,19,20,21]. For example, constrained lab environments inherently limit the study of patterns of space use or social behavior that contribute to social organization in natural populations, yet the consequences of a population’s social organization on individuals is increasingly recognized as a key factor shaping lifetime patterns of health and fitness [22,23,24,25]. At the same time, there have been repeated calls to study traditional model organisms, especially mice, under more natural contexts [26,27,28,29,30,31,32,33]. Indeed, studies of lab mice in outdoor enclosures have already revealed effects of more natural conditions relative to traditional laboratory housing on traits including foraging, hippocampal neurogenesis, immunity, microbiome, and cancer progression [34,35,36,37,38,39,40,41], though we do not yet know the social structure of lab mice under semi-natural field conditions.

Over the last decade, multiple groups have developed high-throughput lab assays using modestly-sized arenas or interconnected cages that allow for increased complexity of social interactions [42,43,44]. Yet, even relatively large and enriched lab settings [45,46,47] fail to capture many of the relevant features of social structures inferred by studies of wild mouse populations to be important to mouse natural history, such as territorial social organization and space use [48,49,50,51,52,53]. Accumulating evidence suggests that generations of inbreeding and artificial selection in lab mouse strains has impacted their social behavior and interactions in small groups [48]. Indeed, even simple lab behavioral assays can reveal differences in classical lab strain and wild-derived mouse behavior [54,55,56,57,58,59]. But we have little understanding of the ecological validity of such studies on lab mice because of the lack of studies addressing the socioecology of lab mice under natural or semi-natural field conditions. At stake is not only our understanding of the ecological validity of studies of social behavior in constrained conditions, but also how we interpret and understand the difference between domesticated lab and genetically wild mice more generally.

An immediate solution is to study the social and spatial behavior of lab mice in large naturalistic spaces. Providing ample space for individuals to interact or avoid each other is critical for assessing social structures because it allows animals to freely express their preferences. There is a long history of studies utilizing semi-natural indoor or large outdoor enclosures to study the population biology of house mice under free-range conditions [60,61,62,63,64,65,66,67,68,69,70,71]. These studies tend to use feral or wild-derived populations of outbred house mice and find that male mice establish and aggressively defend territories occupied by several females and their offspring. Fully adult males are most often associated with high quality territories, while juveniles and subadults typically aggregate in lower quality spaces within the environment [61, 72, 73]. Adult females may compete for nest sites but form strong associations and may even co-nest with close relatives [74, 75]. The competitive environments generated under free-range conditions can have a strong impact on social behavior relative to the lab environment [65, 76, 77], motivating our study of mice in field enclosures where they can freely compete in a semi-natural environment.

Understanding how lab mice behave under more natural field conditions will inform how their behaviors have evolved relative to wild mice, providing a crucial piece of biological data for the best studied mammalian genotype, C57BL/6J (hereafter “C57”). Here, we report the space use and social behavior of replicated mixed-sex populations of C57 and wild-derived outbred (hereafter “WD”) mice in large outdoor field enclosures. We aimed to address three empirical questions with this study. (1) We sought to determine whether male C57 would establish and defend territories in the field, as opposed to generating societies with a single integrated dominance hierarchy (a common outcome in the lab [47, 78,79,80]). We also compare the territorial behaviors of C57 mice and WD outbred mice. (2) We sought to test whether space use and the resulting social relationships that form differ between C57 and WD outbred females. Female fitness in wild mice depends on their abilities to compete for nesting sites [61, 74, 81, 82], and they show spatial and social biases toward other related females [75]. We hypothesized that lab mouse husbandry practices, which regularly lead to isogenic females living and occasionally breeding at very high density [1], may have led to changes in the magnitude of social biases towards familiar females. (3) Finally, we sought to describe the social networks that emerged from these individual behaviors and how differences in behavior between C57 and WD mice scaled up to shape the larger social structure. We hypothesized that, despite the myriad complexities and idiosyncrasies of individual decisions and environmental fluctuations that influence societies, repeated studies of societies that differed in a salient manipulation—in this case the genotype of the individuals—would result in consistently different outcomes across treatments. If this hypothesis is correct, populations with similar initial ecological and demographic conditions will reliably generate similar social structures, suggesting that the biological basis of social organization is amenable to study.

Results

Over a 3-month period (June 2020–August 2020), we examined the emergent social organization generated in enclosures stocked with adult male (n = 10 per trial) and adult female (n = 10 per trial) C57 (n = 4 trials) and WD outbred house mice (n = 3 trials) over the course of 10-day trials. The study employed a common garden design and compared the behavior of C57 mice to wild-derived house mice originating from upstate New York. All mice were bred and reared in the same lab colony prior to release. A detailed discussion of the study design, including the density of mice, length of trials, and genotypes tested is provided in the “Methods” section. Each field enclosure (38.1 m × 15.24 m; ~570 m2) contained eight weather protected resource zones arranged in a 2 × 4 grid pattern (Fig. 1; Additional file 1: Fig. S1A-B), resulting in a slight excess of mice of each sex relative to the available resource zones. Mice were implanted with PIT tags and activity at the resource zones was monitored via an RFID antenna (see the “Methods” section). We obtained high density sampling of mouse RFID reads for all trials (1,307,712 ± 135,646 reads per trial; mean ± SEM) and a mean of 6771 ± 275 RFID reads per mouse per day. Mice were able to quickly traverse the distance between the zones despite the ground vegetation (minimum inter-zone travel time = 10 s; Additional file 1: Fig. S1D). We estimated individual mouse location for a total of 5653 h across all trials (mean = 808 ± 51 h/trial). On average, we inferred individual mice spent 4.2 ± 0.1 h/day at resource zones (range: 1.1 s–18.7 h).

Fig. 1
figure 1

Experimental design for replicate populations of C57 and wild-derived mice in field enclosures. Photos demonstrate the layout of the field enclosures and the eight resource zones arranged in a 2 × 4 grid pattern. Resource zones had a single entrance tube and food and water towers provisioned ad libitum. We observed a variety of behaviors in the resource zones including co-feeding between females, mating and courtship, social investigation, and male-directed aggression towards intruders. Top left image created with BioRender

Territory establishment is slower in C57 compared to wild-derived males

Territorial behavior is reported from wild mouse populations and field enclosure studies at moderate density or indoor environments with high physical complexity [63, 72, 83,84,85]. Territorial males attempt to exclude other males from the spaces that they control, and the ability to compete for and maintain a territory is a key driver of male fitness in freely mating mouse populations [70]. In contrast, in laboratory and natural settings where there is low environmental complexity such that individual dominant mice can readily patrol most of the available space, an alternative social structure predominates where groups of male mice form a social dominance hierarchy [47, 78, 79, 86]. Mouse social hierarchies can vary in the extent of despotism or egalitarianism [47], though are generally characterized by a dominant alpha who is aggressive to all others and has broad access to space, followed by a linear hierarchy of subordinates [80]. Commensal mouse populations may show a mixture of these social forms, with hierarchies formed among the males residing within a single territory or deme but are nevertheless characterized by restricted movement of mice between territorial spaces.

We sought to determine whether C57 and WD males in our enclosures formed territories or dominance hierarchies. If males form territories in our enclosures, a subset of the males in each trial would each monopolize one or more resource zones to the relative exclusion of others. What is more, no male would regularly access all resource zones. In contrast, if males establish a dominance hierarchy within the enclosure, the top-ranking male or males are expected to have relatively free reign and regularly access all or nearly all resource zones. Consistent with previous field enclosure studies of wild-derived mice [72], we find evidence that males of both genotypes established and defended territories, though the intensity of territorial behaviors was muted in C57 relative to the WD outbred mice. Three related, but distinct, analyses support this conclusion.

First, males appear to establish monopolized spaces and rarely visit resource zones that they do not control (Fig. 2A; Additional file 2: Fig. S2A). By the third day, 93% ± 1 of all male-sourced RFID reads within each resource zone belonged to a single male, with zones in WD trials becoming monopolized by males more rapidly than in C57 trials over the 10-day period (t500 = −2.93, P = 0.004; Fig. 2B). No individual male consistently accessed all or nearly all resource zones (range of mean daily zones visited = 1.0–3.9, median = 1.6). Instead, males spent the vast majority of their time within their single most visited resource zone, with WD males developing this site fidelity more rapidly than C57 males (t67 = 2.31, P = 0.02; Additional file 2: Fig. S2B). Thus, male space use was consistent with territory defense, and inconsistent with a single integrated dominance hierarchy.

Fig. 2
figure 2

C57 and wild-derived male mice establish territories rather than dominance hierarchies. A Males generally visited only one or two zones, and no males visited all zones, as shown by representative space use data from three C57 and three WD males (rows, as indicated by three letter individual identification codes). Each 2 × 4 grid shows the eight resource zones on each of the 10 days (columns) of a trial, with warmer colors indicating a high concentration of RFID reads in a given zone. B WD males form territories more rapidly than C57 males. Here, the y-axis represents the percent male-sourced RFID reads per zone from the focal zone’s most present male, with values of 100% indicating complete monopolization. C Male-male dyadic spatiotemporal association bout events were shorter, less frequent, and deteriorated more quickly in WD mice compared to C57 mice. D To quantify territoriality, we calculated a priority access score, a cumulative metric of resource zone access across the 10 days of the trial. For both WD and C57 males, this territory metric deviated from unimodality, indicating a group of males that had consistent access to resource zones (territory holders) and those without (non-territorial). Here, scores near or greater than +10 indicate a male that controlled a resource zone for the duration of the trial. Negative scores are indicative of males that failed to capture a territory. E Consistent with territoriality, but not a dominance hierarchy, territory-holding males were much more likely to win a spatial dispute with another male, defined as displacing the other male, when the dispute occurred in their home territory as opposed to in a different male’s territory. F Territorial control (as categorized by having a day 10 PAS > 0) conferred benefits in the form of increased access to females. This benefit of territoriality was stronger among WD as compared to C57 males (notice that WD males without a territory essentially never spend any time with females, dashed blue line). Data are plotted as means ± SEM

Second, males spend limited time in spatiotemporal overlap with other males. In all trials, males were placed into the enclosure within a resource zone along with their cage mate (brother). Males therefore overlapped in space at the start of the trials, but episodes of male spatiotemporal overlap rapidly decreased in both frequency (t58 = −9.47, P < 0.0001; Additional file 2: Fig. S2C) and duration (t2116 = −13.15, P < 0.0001; Fig. 2C). The collapse of previously existing cage mate relationships was remarkably rapid in WD males. At least half of all the estimated male dyadic spatiotemporal overlap time in WD mice elapsed in the first 54 min of a trial on average (range = 32–93 min), compared to 590 min on average (range = 210–1501 min; P = 0.02) for C57 males. Male-male spatial overlaps shortened in duration over the course of the trial more quickly in WD males as compared to C57 males (t2116 = 2.77, P = 0.006; Fig. 2C) and were less frequent overall (t27 = −4.14, P = 0.0003; Additional file 2: Fig. S2C), suggesting that the intensity of male-male competition in C57 may be weaker than in genetically wild mice. Strikingly, on days 5–10, when territories were clearly established, the amount of temporal overlap between a dominant territorial male in a given zone and other males that visited that same zone was 86% less than expected if each males’ space use was independent of each other (range across trials = 74–97% less; t6 = −58, P < 0.0001). This pattern could be explained either by (a) territorial males quickly expelling intruders when they co-occur or (b) non-territorial males biasing their visitation behavior to times that they expect the territorial male to be absent.

To characterize the variation in males’ territorial behavior and resource access, we developed a priority access score (PAS) metric which tracked changes over time in the degree to which mice monopolized access to resource zones relative to same sex conspecifics (see the “Methods” section). Briefly, for the 10-day trials reported here, strongly positive final scores (near +10) indicate an individual consistently monopolized a single resource zone on each day of the trial while strongly negative scores (near -10) indicate an individual was consistently excluded from most spaces. High scores (>>10) indicate individuals who monopolized more than one zone over the course of the 10-day trial. Scores closer to zero indicate individuals that share spaces to some extent with other individuals of the same sex. Consistent with the hypothesized presence of territorial and non-territorial males, we find a bimodal distribution of males with high and low PAS values (excess mass test for unimodality: P < 0.02; Fig. 2D; Additional file 2: Fig. S2D; in contrast, females of both genotypes do not strongly monopolize resource zones at this density, see Additional file 3: Fig. S3B-C).

The third and ultimate indicator of territory formation rather than a single dominance hierarchy would be for contest outcomes to be predicated on spatial ownership. In other words, territorial males should win competitive interactions in their territory and lose in other males’ territories, such that encounter outcomes depend on location. In contrast, under a social dominance hierarchy, a given male should win or lose depending on their rank without respect to location. To assess this prediction, we investigated male-male dyadic interactions within territories by identifying when the territory-holding male spatiotemporally overlapped with an intruder. For each such inferred overlap, we identified which animal was the first to leave the interaction (i.e., was displaced and lost acute resource access) and assigned that male as the “loser” of the interaction event. We assigned the male remaining in the zone as the “winner”.

We identified 1380 two-male interactions which, as expected, overwhelmingly involved territory holders (n = 1290 events, 93.5% across all trials; mean = 25.6, median = 14, 95% C.I. [15.8, 35.2] per territory holding male). To compare win and loss rates for the same individuals, we examined 32 territorial males that were observed to engage in both “home” and “away” displacement events (n = 23 C57, n = 9 WD; n.b.- many territorial males were not observed in displacement events in an “away” context). We find that overall win rates of territory holders at home (0.82 ± 0.03) are dramatically higher than win rates away (0.18 ± 0.05; t31 = 12.19, P < 0.0001; Fig. 2E), with only one male showing a lower win rate at home than away. The one male that had an excess win rate away was observed to participate in only a single away displacement event, which he won. More than half of the males (59%, 19/32) lost all their contests when away from their territories. Winning away was most common when boundaries were still in flux; 72% (50/69) of displacement events won away occurred on the first day of the trial. Together, these results support the territory hypothesis and reject a hypothesis where male mice in the enclosures formed a single integrated dominance hierarchy.

Territorial control was associated with increased access to females. Males with positive (> 0) final PAS values spent more total time with females than males with negative (< 0) final PAS values over the course of the trial (t661 = −3.6, P = 0.0004; Fig. 2F). The effect of territoriality on males’ access to females was stronger among WD males. Among WD males, those with low PAS values essentially never spent any time with females (range of total time spent with females = 0.73–74.5 min), even on the first day of the experiment. In contrast, although C57 males with low PAS values spend dramatically less time with females over the course of the experiment as compared to high PAS males, they spend comparable amounts of time with females during the first day of the experiment, suggesting that the competitive exclusion of males from access to females was slower to develop in C57 mice.

Establishment of non-overlapping territories led to similar patterns of space use in both C57 and WD males. Examples of males’ spatial behavior are shown in spatial heatmaps in Fig. 2A highlighting the spatial separation and fidelity of males in a single trial (see Additional file 2: Fig. S2A for full example trials). The males in each trial that were not able to capture a resource zone (those with priority access scores less than 0) appeared to adopt an alternative strategy in which they briefly visited several zones each day. These non-territorial males, who were not tied to a specific space that they needed to defend, tended to visit more zones over time on average as compared to the territorial males (t65 = 2.61, P = 0.01; Additional file 2: Fig. S2F).

Distinct patterns of space use and social associations in C57 females

Studies of free-living house mice document competition among females for nesting space. Females tend to tolerate their relatives while avoiding or even showing aggression towards unrelated females [75, 81, 82, 87]. At the same time, male territorial structure may shape female behavior as novel males can lead to pregnancy termination through the Bruce effect [88, 89] and novel males represent an infanticide threat [88,89,90,91]. Competition for nesting space and avoidance of non-sire males is likely to be most acute among breeding females [71, 92]. Our data examine how females use space when first introduced into a large novel social environment and assess the extent to which their space-use behaviors are influenced by their social environment.

C57 females exhibited substantial differences in space and movement patterns across several measures compared to WD females as well as both male genotypes (Figs. 2A and 3A; Additional file 3: Fig. S3A; Additional file 2: Fig. S2A). Unlike WD females and males of both genotypes, C57 female space use was not limited to a few neighboring resource zones, but instead was widespread across the enclosure space. Across sexes and genotypes, mice visited an average of 2.32 ± 0.11 resource zones per day over the course of the trial, though patterns of zone visits varied over time and among individuals. On average, the number of resource zones that animals visited each day increased as the trials progressed (t133 = 13.3, P < 0.0001), but this increase was driven entirely by the behavior of C57 females (P = 0.52 for non-C57 females; Fig. 3B). Although all mice explored an equivalently low number of resource zones during the first several days in the enclosure, by the fourth day, C57 females had significantly increased exploration of the available zones compared to all other groups (P < 0.05 for daily LMM contrast estimates for days 4–10), which did not differ in their extent of space use.

Fig. 3
figure 3

C57 females show distinct patterns of space use and social associations. A C57 females explore more of the available resource zone spaces than C57 males and WD males and females, as shown by three representative females from each genotype (rows, indicated by three letter individual identification codes). Each 2 × 4 grid shows the eight resource zones on each of the 10 days (columns) of a trial, with warmer colors indicating a high frequency of RFID reads in a given zone. B All sex and genotype combinations visit similar numbers of zones (unique per day) on the first several days of the trial, but C57 females increase the numbers of zones visited relative to any other sex or genotype combination as the trial continues. C C57 females visit more cumulative resource zones over the course of the trial period compared to all other sex and genotype combinations. D WD females show a strong and sustained preference for their cage mates, whereas C57 female preferences for cage mates quickly fall to null expectation levels (as indicated by the dashed line). Note that mice were placed in resource zones with their cage mates on the first day of the trial. Data are plotted as means ± SEM

In addition to visiting more zones on average per day, C57 females visited a greater proportion of all possible zones over the course of the trial (Fig. 3C). By the final day of the trial, C57 females had visited 6.3 ± 0.4 of the 8 available zones, which is more than C57 males (4.3 ± 0.4; t125 = 5.5, P < 0.0001), WD females (4.1 ± 0.5; t7.4 = 3.3, P = 0.045), and WD males (3.5 ± 0.5; t7.6 = 4.3, P = 0.01). Many more C57 females (42%, 16/38) visited all eight resource zones as compared to C57 males (8%, 3/39), WD females (3%, 1/29), and WD males (4%, 1/28) (generalized LMM: P < 0.05 for all comparisons). Finally, C57 females spent less time in their most occupied zone compared to WD females (t18 = 2.3, P = 0.04; Additional file 3: Fig. S3E).

Female house mice exhibit social preferences towards familiar same-sex conspecifics under free-living conditions [75, 81], but the degree to which this happens in lab strains like C57 is unclear. Our results suggest that C57 females are more tolerant towards females in general and less biased towards familiar social partners. Overall, C57 mice engage in longer female-female spatiotemporal association bouts compared to WD mice over the course of the trial (t11910 = −2.8, P = 0.02; Additional file 3: Fig. S3F). At the same time, the female preference for spending time with cage mates on days 2–10 of the trials was much greater in WD females than in C57 females (t5 = 4.8, P = 0.005). WD females (t5 = 7.2, P = 0.0007), but not C57 females (t5 = 1.1, P = 0.33), differed significantly from the null value based on the number of cage mates present within a trial (Fig. 3D). Indeed, WD females spent nearly all their female-female social time with cage mates after day 2, whereas C57 females tended to spend less than half of their time with cage mates (Additional file 3: Fig. S3G). These results suggest that inbreeding and/or selection from lab mouse colony rearing practices may have a profound effect on both female space use and social behavior.

Repeatable differences in social structure between C57 and wild-derived mouse societies

Analyses of the behavior of each sex show generally similar patterns of space use for males but strongly divergent patterns for females. We next sought to understand how behavioral differences measured for individuals scaled up to larger patterns of social association and the emergent structure of societies.

We first examined how mice overlapped in space and time to determine to what extent individuals associated, as well as the range of spatiotemporal group compositions that arose. For each trial, we estimated the time spent in group compositions of differing numbers of males and females (Fig. 4A). For this analysis, we considered each mouse’s time separately such that a pair of mice in a zone together for 60 min is recorded as two mouse hours.

Fig. 4
figure 4

C57 social structures are more interconnected than those in wild-derived mouse societies. A C57 mice spend less time alone and have larger numbers of social participants in association events as compared to WD mice. The contour plot shows the average duration of observed time spent in different male and female group compositions across trials for C57 (left) and WD (right) mice. The scores reported here are the average number of hours estimated per trial. A white square indicates that a particular combination was not observed. B C57 females spend less than half of their total observed time alone throughout the trial, less than any other group. C All sex and genotype combinations combinations decrease the time spent in same-sex groups over the course of the trial, with males of both genotypes spending minimal amounts of their observation time in same-sex interactions. D, E Daily social networks from an example C57 trial (D) demonstrates a typical pattern of persistently high female interconnectivity while an example WD trial (E) demonstrates increasing network modularity over time. The size of connections between nodes represents the edge weight whereas the size of nodes reflects the node edge strength or the sum of all edge weights for a single node. F Number of network components increases in WD, but not C57, social networks over time. G Network edge density increases in C57, but not WD, social networks over time. H C57 females have consistently higher measures of node degree centrality, indicating that they are consistently meeting large numbers of social partners on each day of the trial, relative to other sex and genotype combinations. Data are plotted as means ± SEM

In each trial, mice were alone for most of the time that they were recorded at resource zones (range 50–64% for C57; 70–82% for WD mice, solitary mouse time per trial; P = 0.03 based on the probability of all four C57 trials having the lowest values; Fig. 4A), but the proportion of time that individuals spent alone was strongly predicted by sex and genotype. On average, males spent a greater proportion of recorded time at resource zones alone than females (t1339 = 5.4, P < 0.0001; Fig. 4B). Overall, WD males increased their time alone over the course of the trial (t1339 = −3.1, P = 0.002; Fig. 4B) and were especially likely to spend time alone compared to all other groups; all of them (29/29) spent more than 50% of their total recorded time alone. In comparison, 83% of WD females (25/30), 78% of C57 males (31/40), and only 18% of C57 females (7/40) spent most of their recorded time alone. Given the interest in the biology of social isolation in mice [93,94,95,96,97], it is notable that when given the opportunity to freely interact, many mice spent a significant portion of their observed time alone over the course of their trial.

Though individuals spend a large portion of their time at the resource zones by themselves, we estimated more than 2000 mouse hours of spatiotemporal associations across the seven trials, defined as time with two or more mice at the same zone simultaneously (Additional file 1: Fig. S1F). Dyadic interactions accounted for > 50% of estimated association time in both genotypes (62–80% in C57, 78–86% in WD), though larger aggregations of mice were detected in all trials (Fig. 4A). On average, females spent a greater portion of their recorded time in same-sex associations than males over the course of the trial (t1339 = −6.5, P < 0.0001; Additional file 4: Fig. S4A). Most mice (77%, n = 103/134) spent > 50% of their recorded association time in opposite-sex groups, with males increasing this metric over the course of the trial relative to females (t1165 = 6.6, P < 0.0001; Additional file 4: Fig. S4B).

To investigate the emergent group structure of both genotypes, we analyzed the total and daily social networks formed for each trial. Overall, C57 mice formed more connected networks than WD mice, a difference which was largely driven by high levels of C57 female sociability (Fig. 4D–E). WD networks increased in the number of graph components—the portions of the network disconnected from each other—over time (t61 = 4.86, P < 0.0001; Fig. 4F), echoing the demic structure reported for many wild mouse populations [52, 98, 99]. In contrast, C57 networks became more connected over the course of the trials. Over time, the network edge density—a measure of the proportion of edges actually observed out of all possible edges in the network—increased in C57 social networks, but not in WD networks (t61 = −5.4, P < 0.0001; Fig. 4G). Notably, the two genotypes show completely non-overlapping distributions for network edge density from days 4 to 10 (P = 0.03, non-parametric, each day). The differences between genotypes reported in this section are based on whole network level metrics indicating that the basic overall structure of social organization differs between the genotypes.

Females of both genotypes had high degree centrality measures compared to their respective males, indicating females encountered more unique social partners on each day of the trial. There was a significant three-way interaction between sex, genotype, and time, such that C57 females rapidly increased their network centrality measures compared to all other sex and genotype combinations (t134 = 2.6, P = 0.01; Fig. 4H). Thus, many of the differences we see in social networks between the genotypes is driven by the propensity of C57 females to co-occur with many distinct individuals. By the final day of the trial, C57 mice had overlapped with many more of the available social partners present in the enclosures as compared to WD mice (t5.48 = −3.9, P = 0.01; Additional file 4: Fig. S4C), who failed, on average, to ever meet more than 50% of the potential social partners, at least within or around resource zones. Intriguingly, females of both genotypes showed high levels of vertex page rank scores, indicating that females of both strains maintain higher levels of connectedness over time than their respective males (t1343 = 3.2, P = 0.001; Additional file 4: Fig. S4D).

Discussion

Our replicated semi-natural field enclosure experiments demonstrate that C57BL/6J lab mice broadly recapitulate many of the behaviors of wild-derived outbred mice in free-ranging conditions, including clear evidence of male territories, but have different emergent social structures. This difference is largely due to C57 females being more exploratory and showing less biased patterns of social association. The organization of mammal societies is influenced by ecological [100,101,102], demographic [103,104,105,106], and phylogenetic factors [107,108,109], each of which were controlled for in our trials. Thus, these data show that genotype can have a strong effect on social structures in mammals [110, 111]. In this case, it demonstrates key social behaviors that have been altered during the process of lab mouse domestication [48]. These data also highlight the flexibility of mouse social behaviors across diverse ecological and demographic conditions. For example, in contrast to lab studies at high density, which identify dominance hierarchies among both lab and wild-derived male mice [46,47,48, 112], the males in our lower density populations consistently formed and defended territories regardless of genotype (Fig. 2).

Although our experiment only examined one set of ecological and demographic conditions, it demonstrates a tractable field approach in which variables including food resources, defensibility of spaces, and demographic compositions are all easily tunable. There has been considerable interest in recent years in high-throughput measures of mouse social behavior in lab settings [42, 44,45,46, 113], and our study provides a blueprint for similar studies to be conducted under semi-natural outdoor settings using different socio-ecological conditions or genotypes. Studying reproducible genotypes of mice under similar settings in other locations and seasons provides an exciting avenue to examine the role of environmental variation in shaping individual phenotypes and population dynamics.

What drives the difference that we saw in female space use across our trials (Fig. 3)? Space use in female mammals is often predicted by intra-sexual competition for food resources and nest sites [114], but resource availability and population density were identical across trials in our study. This suggests two non-mutually exclusive possibilities. First, there may be a genetic difference shaping space-use behavior between C57 females and wild-derived outbred females. For example, this could be manifest as an increased tendency to explore irrespective of their social environment. Alternatively, female behavior may respond to the social conditions present in our trials, which differ in some key respects between the two genotypes. C57 females experience a society where everyone has very high genetic similarity, while WD females experience a world with variation in relatedness. It is likely that mouse behavior in general might be sensitive to these parameters, though studies of social interactions in the lab frequently utilize single inbred strains of mice [115,116,117]. Female mice respond to variation in perceived relatedness between themselves and males [68, 118, 119], and thus it may be the case that females explore more in conditions when all social partners are genetically homogeneous but show more restricted space use when social partners are genetically heterogenous. In wild house mice, infanticide risk from both male and female conspecifics is thought to be a major driver of social behavior in females [120,121,122]. As a result, wild female house mice will aggressively defend space from other females [71, 81, 82, 123, 124]. In contrast, C57 mice have been bred to live in cages at high densities, especially among females, and this is associated with lower female aggression compared to wild-derived mouse genotypes [125]. Differences in relative tolerance of other females may be a key driver of the observed differences in social organization between C57 and outbred females in this study. Understanding how innate behavioral differences among genotypes and emergent properties generated by social interactions shape mammalian societies is an exciting future direction, one that could be extended to other commonly commercially available inbred (such as BALB/c) and outbred (such as Swiss Webster and CD-1) mouse strains.

Male space use in rodents and other mammals is frequently linked to patterns of female space use [114, 122]. Yet, despite differing patterns of female space use between genotypes, the male spatial and social structures across genotypes were quite similar, highlighting that some aspects of social organization are relatively less sensitive to other features of a population’s socioecology. Perhaps one of the most striking features of our study is the speed with which male-male social interactions decrease in frequency, especially among wild-derived outbred males (Fig. 2C). Previous studies of wild mouse behavior have reported males will defend territories and attempt to exclude other males [60, 64], and our data show this behavior is retained in male C57 mice. The formation and physiological consequences of dominance hierarchies among male mice have been the subject of recent study in the lab [47, 48, 112, 126], but our results suggest that when given ample and defensible spaces male mice will tend to avoid interacting with others and form individual territories rather than a single integrated dominance hierarchy. Differences between territoriality and dominance behaviors remain poorly understood at the mechanistic level. Social dominance hierarchies and the establishment of territorial boundaries could be mediated by the same physiological and neurological mechanisms, or they might be mediated by distinct mechanisms. The process by which dominant males recognize known, tolerated subordinates may be fundamentally distinct from the process by which a territorial male characterizes intruders as well as neighbors versus strangers [127, 128].

Across repeated trials, we identified differences in the higher-level social organizations of C57 lab mice and their wild-derived outbred counterparts (Fig. 4). Studies of social structures tend to come from idiosyncratic populations living in the wild, meaning that studies of social behavior in natural conditions are rarely replicated [129,130,131]. Studies of free-living populations are critically important, but this non-replicability makes understanding the specific genetic, neurobiological, ecological, and demographic factors influencing complex behavior challenging. The repeatability of social organization demonstrated here suggests that future work manipulating aspects of physiology or neural function in free-range mice will offer a unique opportunity to study not just differences in individual behavior, but also how and whether specific behaviors reliably influence society.

Conclusions

The importance of social interactions in shaping health, physiology, and fitness is increasingly recognized, though traditional biomedical approaches study animals in highly impoverished social conditions. Thus, approaches for studying animal models in more natural and free-ranging settings are needed as a complement to traditional laboratory studies that examine animals in highly socially and spatially constrained conditions. Using a classic common garden study design, we investigated the spatial and social behavior of the best studied genotype of mice in comparison to wild-derived outbred counterparts. We found that while male C57 and WD mice behavior was largely similar in that males of both genotypes formed and defended exclusive territories, male C57 mice did so less rapidly. Importantly, territoriality in our field enclosures was the observed social organization, in contrast to a dominance hierarchy which is more commonly found in studies of mice in laboratory environments. Female C57 mice differed strikingly from their wild-derived counterparts in terms of their increased exploration of the field enclosure, their time spent in association with other mice, and their lack of preference for spending time with previously known same-sex cage mates. Moreover, we show that these differences scale to profoundly shape the larger social networks formed by populations of C57 and WD mice, with C57 networks being persistently highly connected across days. These results indicate that the process of domestication in laboratory mice and selection for social tolerance among females may have profound impacts on the expression of group social behaviors and broader patterns of lab mouse societies. Consistency in individual and overall society-level differences within and between genotypes suggests that the wealth of genetic resources available in mice make them well-suited to address fundamental questions about the genetic basis of higher-level social organization.

Methods

Ethical note

All experimental procedures adhered to guidelines established by the U.S. National Institutes of Health and the ASAB/ABS guidelines for the use of animals in research [132] and have been approved by the Cornell University Institutional Animal Care and Use committee (IACUC: Protocol #2015-0060). All animals were briefly anesthetized and implanted with dual PIT tags (12 mm) prior to introduction to the field enclosures. This procedure is minimally invasive and is consistent with recommendations by the veterinary and animal care staff at Cornell University. PIT tag weight is negligible, and anesthesia removed any discomfort associated with handling and implantation of the tags. We observed no changes in mouse behavior as a result of this procedure, as assessed by normal food and water intake and daily activity on home cage running wheels.

Animals

We examined two groups of M. m. domesticus, C57BL/6J (C57) and wild-derived (WD) outbred mice. Rather than merely qualitatively comparing our data to published data on mouse behavior in field enclosures or semi-natural indoor settings [51, 63, 72, 133, 134], we conducted identical, simultaneous trials using wild-derived outbred mice to allow for a direct quantitative comparison of data collected using the same methods and the same physical environmental conditions. We note that while lab mouse ancestry includes three house mouse subspecies, the C57 genome is > 90% Mus musculus domesticus [135]. More information on the origins of C57 mice can be found elsewhere [136,137,138,139]. Therefore, as a “control” wild-derived mouse population, we used outbred mice derived from wild house mice caught in upstate New York. Thus, the comparison group of mice were born and raised in the lab for multiple generations but have a wild-derived outbred genetic background, following other studies examining genetically “wild” house mouse populations [48, 64, 85, 125]. We initially obtained C57 (#000664) mice from The Jackson Laboratory (Bar Harbor, Maine, USA), which were then maintained and bred in our colony until the start of the experiment. We maintained a wild-derived outbred mouse stock from intercrossing six different wild-caught families generated through distinct initial pairings of wild mice from Saratoga Springs, NY, USA, trapped by MJS in 2013 [140]. Thus, the potential adverse stressful effects of transferring wild-caught mice to laboratory environments for use in our enclosure study is avoided. We maintained an outbred colony with > 16 breeding cages prior to the start of the experiment in 2020. These mice are descended from the same initial collection in Saratoga Springs that gave rise to the wild-derived inbred mouse strains SarA/NachJ (#035346), SarB/NachJ (#035347), and SarC/NachJ (#035348) available from the Jackson Laboratory [141]. However, our mice are not descended from the wild-derived inbred strains deposited at JAX, rather they are simply related to them and share alleles from the wild Saratoga Springs, NY house mouse population. All of the wild-derived outbred mice in our experiment had unique diploid genotypes representative of wild house mouse haplotypes found near Saratoga Springs, New York. All of the mice (both C57 and wild-derived outbred) used in this study were bred in our lab colony, under standard housing conditions contemporaneously in the same room. All mice were between 15 and 28 weeks of age when they were released into the field enclosures (see Additional file 6: Table S2 for detailed information on the ages and litter relationships among mice).

The population density in our enclosure was 0.034 mice/m2, which falls within the range of typical population densities reported for wild mice [142] (~0.0011–0.11 mice/m2) and is the same order of magnitude as initial populations in other field enclosure studies of wild-derived mouse behavior [63, 68].

Study design

Field studies of wild populations provide a powerful means to link aspects of organismal biology to selection but are typically hampered by a lack of replication [129, 143]. Enclosure studies done over short, but biologically relevant, time periods provide an opportunity to observe replicate populations across multiple trials. A fundamental feature of enclosure studies, as compared to wild populations, is that the experimenters must choose among a range of possible resource distributions, population densities, sex ratios, etc. as the starting conditions of the trials. House mice naturally live under a wide range of resource distributions and densities across non-commensal and commensal settings [61, 66, 73], meaning that there exists a large set of naturalistic conditions rather than a single most appropriate condition for studying mice. While studying social organization across a range of socio-ecological starting conditions will be informative, for this study, we chose a single common set of initial resource distribution, density, and sex ratio conditions for all trials. For our trials, we released 10 males and 10 females into semi-natural field enclosures with eight resource zones which we monitored using a commercially available RFID system. This experimental design allowed us to assess the consistency of social behaviors and structures formed within the enclosures across replicates and to isolate the effect of host genetic background on these outcomes.

To initiate each trial, we placed mice implanted with passive integrated transponder (PIT) tags (see RFID analyses section) into one of eight resource zones (see the “Field site” section) with their same-sex sibling cage mates in the evening shortly before sunset, meaning that all individuals started the trials at a resource zone with one or more familiar social partners (see the Additional file 6: Table S2 for details on cage mate and relatedness information). We allowed trials to proceed for 10 days. This trial length was chosen because it allowed us to avoid females giving birth in the enclosures and made it feasible to conduct multiple replicated trials in the same enclosures over the course of a single summer. House mouse social structure has been reported to vary seasonally and with shifts in demographic parameters in natural populations [49, 134, 144], so different conditions or studies conducted at other times of year may yield different results. The experiment reported here focuses on behavior during the initial formation and establishment of social structures within replicated mouse populations in our enclosures while controlling for seasonal and demographic variation. Over the course of 10 days, mice explored the enclosures and resource zones and engaged in a variety of social interactions with conspecifics including courtship, mating, co-nesting, and fighting (Fig. 1; Additional file 7: Video S1). As our goal was to test hypotheses regarding patterns of mouse space use and social structure, we focused our analyses on the large RFID dataset.

Field site

We conducted field work at Cornell University’s Liddell Laboratory Field Station in Ithaca, New York, USA using two adjacent and identically sized outdoor field enclosures. Our enclosures are approximately 9,000 times larger than the area of a standard laboratory mouse cage (Additional file 1: Fig. S1B). The walls of the enclosures were made from sheet metal and stood approximately 1.2 m tall and extended 1.2 m into the ground to prevent the mice from tunneling and moving between the enclosures. Each enclosure was covered with netting to prevent aerial predation, and gravel was spread along the interior perimeter of each enclosure to discourage digging near the walls. Three days prior to releasing mice into the enclosures, we trapped in and around the enclosures to capture and remove any small mammals or snakes from the enclosure. The enclosures contained a mixture of local perennial grasses and plant communities which were mowed to a height of approximately 5 cm at least one day prior to the start of each trial to maintain similar ecological starting conditions across trials.

We supplied all resource zones with food and water accessible by the mice ad libitum. Resource zones were covered with waterproof corrugated roofing material attached to a polyvinyl chloride (PVC, Home Depot, USA) frame to shade the zones during the day and provide protection from rain. Resource zones were comprised of two nested storage tubs (Rubbermaid, USA) and a single PVC entrance tube (50 mm diameter) through which the mice could freely enter or exit the zone. Each resource zone contained feeder towers containing food and water in excess (approximately 50 g of mixed sunflower and bird seed, and 2 L of water), which was replenished within a trial as necessary and replaced between trials. Several pieces of plastic lumber were added to provide edges and elevated locations for the mice to perch within the resource zones. All resource zone materials were thoroughly washed with soap and water and then wiped down with 70% ethanol at the start of the trial and again once all mice were removed from the enclosure to remove to prevent the transfer of odorants between trials in these spaces.

At least 24 h prior to release in the enclosures, all subjects were briefly anesthetized with isoflurane (3-5%) and placed into a stereotaxic frame (Kopf Instruments, Tuhunga, CA, USA). Mice were subcutaneously implanted with dual PIT tags (BioMark, Boise, ID, USA) in the dorsal flank and periscapular region. Each resource zone was equipped with a 15 cm RFID antenna connected to a centralized data acquisition unit (BioMark, Small Scale System, Boise, ID, USA). Antennae were placed directly beneath the floor adjacent to the PVC zone entrance tubes to increase the likelihood of capturing mouse activity during entrances and exits from the resource zone. Scanning for PIT tags within the antenna range occurred at approximately 2–3 Hz continuously for 10 days.

After the 10-day observation period, > 50 live-catch traps (H.B. Sherman, Tallahassee, FL, USA) were baited with nesting material, sunflower seeds, and a moistened cotton ball and placed in a grid pattern in the enclosures in the evening (1900–2200 h) and were checked for occupancy the following morning (0600–0900 h). Trapping continued until all the mice were recovered or identified as deceased or missing (a conclusion reached if there were no RFID reads in the enclosure for a 24-h period after 3 days of trapping; see the Additional file 6: Table S2). The trap locations were recorded, and the individual identities of the mice were confirmed using a handheld RFID reader (BioMark, HPR Lite). Mice were immediately transported to a clean cage with ample nesting material, sunflower seeds, and water. At the conclusion of the experiment, all mice were euthanized using carbon dioxide inhalation followed by decapitation for the collection of tissues for future analyses.

RFID analyses

We monitored the resource zones continuously over the trial period via a RFID antenna placed beneath the sole entrance into the zone (Fig. 1; Additional file 1: Fig. S1C). Since the RFID antenna detected tags both above and below the horizontal plane of the antenna and mice moved underneath the resource zones, we inferred RFID reads to indicate that mice were within a space inclusive of the inside and underside of the storage tubs.

To convert instantaneous RFID reads into estimates of how long mice spent at the resource zones, we grouped RFID reads into discreet resource zone visits with estimated durations (Additional file 1: Fig. S1E-F). As expected, the total number of visits to a zone strongly correlated with the total estimated duration of time spent at a zone (Spearman’s correlation, R > 0.84 for all genotype and sex combinations; Additional file 1: Fig. S1G).

We conducted two types of RFID analyses focusing on (1) when and at what zone each mouse was detected and (2) estimating how long it spent at each zone. In the first set of analyses, we examined whether or not mice visited a given zone during each day of the trial (Figs. 2A and 3A–C; Additional file 1: Fig. S1G; Additional file 2: Fig. S2A&F; Additional file 3: Fig. S3A). We estimated the minimum distance traveled using consecutive transitions in the RFID reads between antennae in different resource zones and the known spatial layout of the enclosure (Additional file 3: Fig. S3D). To assess the degree to which resource zones were exclusively accessed by individual mice, we examined the percent of same-sex sourced RFID reads in zones for each mouse (Fig. 2B).

In our second set of analyses estimating the durations of time mice spent in zones, we first examined the time elapsed between consecutive RFID read events for each mouse within each resource zone (the RFID inter-read interval). We found that the distribution of all RFID inter-read intervals was heavily skewed (min = 1 s, median = 1 s, mean = 16.4 s, max = 32,683 s; Additional file 1: Fig. S1E). We grouped RFID reads for each mouse for each zone into visitation bouts using a 139 s (the cut-off for capturing 99% of all the within-mouse inter-read interval values within a single zone) sliding window method such that an RFID read which occurred within 139 s of the previous read extended the visit bout duration up to that read. Transitions between zones automatically ended and started bouts in the first and second zones, respectively. We assigned isolated single RFID reads (with no other reads within 139 s before or after the focal read) a visit bout duration of 1 s. Using this dataset, we estimated the time mice spent in resource zones (Additional file 2: Fig. S2B). Mice were defined as participating in a spatiotemporal association bout when two or more mice had overlapping resource zone visit bouts within the same zone (see Additional file 1: Fig. S1F for a graphical schematic of the approach). We used the spatiotemporal association bout dataset to estimate the duration and frequency of male-male (Fig. 2C), female-male (Additional file 2: Fig. S2E), and female-female (Additional file 3: Fig. S3F) dyadic association bouts. We omitted a subset of animals from a subset of days for all spatial and social analyses when they received no reads over multiple days and were presumed dead (see Additional file 6: Table S2).

Priority access score calculation

Priority access scores were calculated separately for male and female mice within a trial. First, we calculated the time a given mouse (M) occupied a resource zone (Z) as a percentage of the total time that zone was occupied by same-sex conspecifics on a given day (D).

$${Occupancy}_{M, D, Z}= \frac{{time}_{M,D, Z}}{\sum_{m=1}^{10}{time}_{m,D,Z}}$$

Next, we calculated a daily Capture Score by summing the Occupancy values for all available zones. Mice that did not have an occupancy value of greater than 0.5 (in other words, a majority share of the time spent in any zone), were penalized by subtracting 1 from the final Capture Score. The penalty indicates that on a given day, a mouse failed to capture any of the zones that mouse visited.

$${Capture\, Score }_{M,D}= \left\{\begin{array}{ll}\sum\limits_{z=1}^{8}{Occupancy}_{M,D,z}& {\text{if}}\, {\exists\, Occupancy}_{M,D,z}>0.5 \,\forall \,z=1, \dots ,8 \\ \sum\limits_{z=1}^{8}{Occupancy}_{M,D,z}-1& otherwise\end{array}\right.$$

To see how access to zones changed over time, we took the cumulative sum of an individual’s Capture Score ordinally across each day of the trial to derive an evolving priority access score on each day of the trial.

$${Priority \,Access \,Score}_{M,D}= \sum\limits_{d=1}^{D}{Capture\, Score}_{M,d}$$

As an example, if one male (male A) occupied a single resource zone every day of the trial for 4 h a day, while another male (male B) accessed only that same zone for 1 h per day, and the zone was visited by no other mice, male A’s daily Capture Score would equal 0.8, (because he controlled 4 out of 5 h), while male B’s daily Capture score would equal − 0.8 (because he controlled 1 out of 5 h and received a one-point penalty for not controlling any zones). If this pattern of visitation remained unchanged for all 10 days, then male A’s final priority access score would equal 8, while male B’s priority access score would equal − 8. The priority access score value thus provides a temporally evolving measure that captures the dynamics of territory formation, maintenance, and collapse (Additional file 2: Fig. S2D; Additional file 3: Fig. S3B).

Male spatiotemporal overlap and win/loss at home and away

We used the zone visit bout and spatiotemporal association bout dataset to investigate the degree to which territorial and intruder males avoided spatiotemporal overlap. Across all zones for the last six days of the trial, we independently calculated the total time territorial and intruder males spent in zones and the amount of that time that was spent in spatiotemporal overlap. We derived null expected values of territorial and intruder male spatiotemporal overlap ((total territorial male hours in zones × total intruder male hours in zones)/(24 h × 6 days × 8 zones)) and examined the percent difference between the null expectation and the observed time territorial and intruder males overlapped. We chose the last six days of the trials for this analysis as territorial relationships were clearly established by this point in the trial.

Next, we used the RFID and spatiotemporal association bout data to estimate the dynamics of male spatiotemporal overlap and territorial intrusion. We inferred contests by identifying instances in the dataset where a single male within a zone was joined by an additional male, followed by one of the males leaving. We identified the male participants as either the territory holder of that zone or as an intruder based on whether the focal male had captured > 50% of all of the male-derived RFID reads within the zone during the last 24 h. We derived “win” and “loss” rates both within and away from a home territory for each male participant on each day of the trial by assigning the remaining male and the male that left as the “winner” and “loser” of the contest, respectively (Fig. 2).

Female spatiotemporal overlap and cagemate associations

Using the known cage mate relationships between individual females, we examined the total time female mice spent in spatiotemporal association with cage mate (sibling) versus non-cage mate females (Additional file 3: Fig. S3G). Based on the number of cage mates and non-cage mate female social partners available on each day of the trial (which fluctuated for some trials based on mice disappearing from the trials, see the Additional file 6: Table S2), we calculated a daily null value of the expected time spent with cage mates and the percent observed value deviation from the expected separately for each female (Fig. 3D).

Social network analyses

We used the spatiotemporal association bout data set arranged in a group by individual matrix to construct social networks. Daily weighted adjacency matrices were derived from a Simple Ratio Index calculation based on counts of binary participation in spatiotemporally overlapping mouse zone visit bout events (where zero indicates that individuals were never in the same spatiotemporal association bout and one indicates that individuals shared all spatiotemporal association bouts; Additional file 1: Fig. S1F) using the asnipe [145] package in R 4.1.2 (R Development Core Team). All networks were constructed using the weighted matrices and the igraph [146] package in R. Node sizes in our network reflect the node edge strength, or the sum of all connected edges to the node, while edge widths reflect the edge weight derived from the daily weighted adjacency matrices. If an animal died or went missing, we no longer included it in a given trial’s daily network [147, 148] (see Additional file 6: Table S2).

Statistical analyses

We built mixed effects models using R 4.1.2 (R Development Core Team) and the R packages lme4 [149], lmerTest [150], and emmeans [151] to examine relationships between predictor and response variables. We include the full statistical tests and model outputs for all analyses in the Additional file 6: Table S2. Most analyses are conducted with a repeated measures design where data are examined for each day per mouse or other relevant unit as appropriate using mixed effects models. Across models we considered random effects of trial, mouse, and/or time in the trial. We included relevant random intercepts and random slopes in our models as appropriate—random slopes of time were generally included, but we occasionally needed to simplify random effects structures when necessary to avoid singular fits. All model formulae, including random effects structures, are explicitly reported in Additional file 5: Table S1. For testing for unimodality, we used the multimode [152] package with “ACR” method [153]. Graphing was done in R using the package ggplot2 [154] and in GraphPad Prism 9.3 (www.graphpad.com). We report all means ± standard error (SEM), unless otherwise stated, and consider all values statistically significant when P < 0.05.

Availability of data and materials

All statistical outputs are included in Additional file 5: Table S1. Metadata and trial information is included in Additional file 6: Table S2. Summary data files are deposited on Zenodo (DOI: 10.5281/zenodo.7655040). Code and scripts for reproducing figures and analyses are available on Zenodo (DOI: 10.5281/zenodo.10402900).

Change history

Abbreviations

C57:

C57BL/6J

WD:

Wild-derived outbred

PAS:

Priority access score

LMM:

Linear mixed effects models

SEM:

Standard error measure

References

  1. Fox JG, Barthold S, Davisson M, Newcomer CE, Quimby FW, Smith A. The mouse in biomedical research: normative biology, husbandry, and models. Amsterdam: Elsevier; 2006.

  2. Rosenthal N, Brown S. The mouse ascending: perspectives for human-disease models. Nat Cell Biol. 2007;9:993–9.

    Article  CAS  PubMed  Google Scholar 

  3. Snell GD. Biology of the laboratory mouse. 1941.

    Google Scholar 

  4. Waterston RH, Pachter L. Initial sequencing and comparative analysis of the mouse genome. Nature. 2002;420:520–62.

    Article  CAS  PubMed  Google Scholar 

  5. Taylor K, Gordon N, Langley G, Higgins W. Estimates for worldwide laboratory animal use in 2005. Altern Lab Anim. 2008;36:327–42.

    Article  CAS  PubMed  Google Scholar 

  6. Phifer-Rixey M, Nachman MW. Insights into mammalian biology from the wild house mouse Mus musculus. eLife. 2015;4:e05959.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Paigen K. One hundred years of mouse genetics: an intellectual history. I. The classical period (1902–1980). Genetics. 2003;163:1–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Paigen K. One hundred years of mouse genetics: an intellectual history. II. The molecular revolution (1981–2002). Genetics. 2003;163:1227–35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Blake JA, Baldarelli R, Kadin JA, Richardson JE, Smith CL, Bult CJ, et al. Mouse Genome Database (MGD): knowledgebase for mouse–human comparative biology. Nucleic Acids Res. 2021;49:D981–7.

    Article  CAS  PubMed  Google Scholar 

  10. Luo Y, Hitz BC, Gabdank I, Hilton JA, Kagda MS, Lam B, et al. New developments on the Encyclopedia of DNA Elements (ENCODE) data portal. Nucleic Acids Res. 2020;48:D882–9.

    Article  CAS  PubMed  Google Scholar 

  11. Schaum N, Karkanias J, Neff NF, May AP, Quake SR, Wyss-Coray T, et al. Single-cell transcriptomics of 20 mouse organs creates a Tabula Muris. Nature. 2018;562:367–72.

    Article  PubMed Central  Google Scholar 

  12. Stamatoyannopoulos JA, Snyder M, Hardison R, Ren B, Gingeras T, Gilbert DM, et al. An encyclopedia of mouse DNA elements (Mouse ENCODE). Genome Biol. 2012;13:418.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Dickinson ME, Flenniken AM, Ji X, Teboul L, Wong MD, White JK, et al. High-throughput discovery of novel developmental phenotypes. Nature. 2016;537:508–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Groza T, Gomez FL, Mashhadi HH, Muñoz-Fuentes V, Gunes O, Wilson R, et al. The International Mouse Phenotyping Consortium: comprehensive knockout phenotyping underpinning the study of human disease. Nucleic Acids Res. 2023;51:D1038–45.

    Article  CAS  PubMed  Google Scholar 

  15. Bogue MA, Ball RL, Philip VM, Walton DO, Dunn MH, Kolishovski G, et al. Mouse Phenome Database: towards a more FAIR-compliant and TRUST-worthy data repository and tool suite for phenotypes and genotypes. Nucleic Acids Res. 2022;51:D1067–74.

    Article  PubMed Central  Google Scholar 

  16. Graham AL. Naturalizing mouse models for immunology. Nat Immunol. 2021;22:111–7.

    Article  CAS  PubMed  Google Scholar 

  17. Jaric I, Voelkl B, Clerc M, Schmid MW, Novak J, Rosso M, et al. The rearing environment persistently modulates mouse phenotypes from the molecular to the behavioural level. PLoS Biol. 2022;20:e3001837.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Lahvis G. Unbridle biomedical research from the laboratory cage. eLife. 2017;6:e27438.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Leonelli S, Ankeny RA. What makes a model organism? Endeavour. 2013;37:209–12.

    Article  PubMed  Google Scholar 

  20. Cait J, Cait A, Scott RW, Winder CB, Mason GJ. Conventional laboratory housing increases morbidity and mortality in research rodents: results of a meta-analysis. BMC Biol. 2022;20:15.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Bailoo JD, Murphy E, Boada-Saña M, Varholick JA, Hintze S, Baussière C, et al. Effects of cage enrichment on behavior, welfare and outcome variability in female mice. Front Behav Neurosci. 2018;12:232.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. House JS, Landis KR, Umberson D. Social relationships and health. Science. 1988;241:540–5.

    Article  CAS  PubMed  Google Scholar 

  23. Sapolsky RM. The influence of social hierarchy on primate health. Science. 2005;308:648–52.

    Article  CAS  PubMed  Google Scholar 

  24. Sapolsky RM. Social status and health in humans and other animals. Annu Rev Anthropol. 2004;33:393–418.

    Article  Google Scholar 

  25. Kikusui T, Winslow JT, Mori Y. Social buffering: relief from stress and anxiety. Philos Trans R Soc B Biol Sci. 2006;361:2215–28.

    Article  CAS  Google Scholar 

  26. Shah S. The case for free-range lab mice. New York: The New Yorker; 2023.

  27. Lipp H-P, Wolfer DP. Natural neurobiology and behavior of the mouse: relevance for behavioral studies in the laboratory. In: Behavioral genetics of the mouse, Vol. 1: Genetics of behavioral phenotypes. New York: Cambridge University Press; 2013. p. 5–16.

  28. Dennis EJ, Hady AE, Michaiel A, Clemens A, Tervo DRG, Voigts J, et al. Systems neuroscience of natural behaviors in rodents. J Neurosci. 2021;41:911–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Datta SR, Anderson DJ, Branson K, Perona P, Leifer A. Computational neuroethology: a call to action. Neuron. 2019;104:11–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Freund J, Brandmaier AM, Lewejohann L, Kirste I, Kritzler M, Krüger A, et al. Emergence of individuality in genetically identical mice. Science. 2013;340:756–9.

    Article  CAS  PubMed  Google Scholar 

  31. Marashi V, Barnekow A, Sachser N. Effects of environmental enrichment on males of a docile inbred strain of mice. Physiol Behav. 2004;82:765–76.

    Article  CAS  PubMed  Google Scholar 

  32. Richter SH, Garner JP, Würbel H. Environmental standardization: cure or cause of poor reproducibility in animal experiments? Nat Methods. 2009;6:257–61.

    Article  CAS  PubMed  Google Scholar 

  33. Weber EM, Dallaire JA, Gaskill BN, Pritchett-Corning KR, Garner JP. Aggression in group-housed laboratory mice: why can’t we solve the problem? Lab Anim. 2017;46:157–61.

    Article  Google Scholar 

  34. Cope EC, Opendak M, LaMarca EA, Murthy S, Park CY, Olson LB, et al. The effects of living in an outdoor enclosure on hippocampal plasticity and anxiety-like behavior in response to nematode infection. Hippocampus. 2019;29:366–77.

    Article  PubMed  Google Scholar 

  35. Daan S, Spoelstra K, Albrecht U, Schmutz I, Daan M, Daan B, et al. Lab mice in the field: unorthodox daily activity and effects of a dysfunctional circadian clock allele. J Biol Rhythms. 2011;26:118–29.

    Article  PubMed  Google Scholar 

  36. Leung JM, Budischak SA, The HC, Hansen C, Bowcutt R, Neill R, et al. Rapid environmental effects on gut nematode susceptibility in rewilded mice. PLoS Biol. 2018;16:e2004108.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Lin J-D, Devlin JC, Yeung F, McCauley C, Leung JM, Chen Y-H, et al. Rewilding Nod2 and Atg16l1 mutant mice uncovers genetic and environmental contributions to microbial responses and immune cell composition. Cell Host Microbe. 2020;27:830-840.e4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Yeung F, Chen Y-H, Lin J-D, Leung JM, McCauley C, Devlin JC, et al. Altered immunity of laboratory mice in the natural environment is associated with fungal colonization. Cell Host Microbe. 2020;27:809-822.e6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Makin DF, Agra E, Prasad M, Brown JS, Elkabets M, Menezes JFS, et al. Using free-range laboratory mice to explore foraging, lifestyle, and diet issues in cancer. Front Ecol Evol. 2021;9:741389.

    Article  Google Scholar 

  40. Dell’Omo G, Ricceri L, Wolfer DP, Poletaeva II, Lipp H-P. Temporal and spatial adaptation to food restriction in mice under naturalistic conditions. Behav Brain Res. 2000;115:1–8.

    Article  PubMed  Google Scholar 

  41. Vyssotski AL, Dell’Omo G, Poletaeva II, Vyssotski DL, Minichiello L, Klein R, et al. Long-term monitoring of hippocampus-dependent behavior in naturalistic settings: mutant mice lacking neurotrophin receptor TrkB in the forebrain show spatial learning but impaired behavioral flexibility. Hippocampus. 2002;12:27–38.

    Article  CAS  PubMed  Google Scholar 

  42. Puścian A, Łęski S, Kasprowicz G, Winiarski M, Borowska J, Nikolaev T, et al. Eco-HAB as a fully automated and ecologically relevant assessment of social impairments in mouse models of autism. eLife. 2016;5:e19532.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Körholz JC, Zocher S, Grzyb AN, Morisse B, Poetzsch A, Ehret F, et al. Selective increases in inter-individual variability in response to environmental enrichment in female mice. eLife. 2018;7:e35690.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Zocher S, Schilling S, Grzyb AN, Adusumilli VS, Bogado Lopes J, Günther S, et al. Early-life environmental enrichment generates persistent individualized behavior in mice. Sci Adv. 2020;6:eabb1478.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Torquet N, Marti F, Campart C, Tolu S, Nguyen C, Oberto V, et al. Social interactions impact on the dopaminergic system and drive individuality. Nat Commun. 2018;9:3081.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Weissbrod A, Shapiro A, Vasserman G, Edry L, Dayan M, Yitzhaky A, et al. Automated long-term tracking and social behavioural phenotyping of animal colonies within a semi-natural environment. Nat Commun. 2013;4:2018.

    Article  PubMed  Google Scholar 

  47. Williamson CM, Lee W, Curley JP. Temporal dynamics of social hierarchy formation and maintenance in male mice. Anim Behav. 2016;115:259–72.

    Article  Google Scholar 

  48. Zilkha N, Chuartzman SG, Sofer Y, Pen Y, Cum M, Mayo A, et al. Sex-dependent control of pheromones on social organization within groups of wild house mice. Curr Biol. 2023;33:1407-1420.e4.

  49. Bronson FH. Reproductive ecology of the house mouse. Q R Biol. 1979;54:265–99.

    Article  CAS  Google Scholar 

  50. Hurst JL. The functions of urine marking in a free-living population of house mice, Mus domesticus Rutty. Anim Behav. 1987;35:1433–42.

    Article  Google Scholar 

  51. König B, Lindholm AK, Lopes PC, Dobay A, Steinert S, Buschmann FJ-U. A system for automatic recording of social behavior in a free-living wild house mouse population. Anim Biotelemetry. 2015;3:39.

    Article  Google Scholar 

  52. Pocock MJO, Searle JB, White PCL. Adaptations of animals to commensal habitats: population dynamics of house mice Mus musculus domesticus on farms. J Anim Ecol. 2004;73:878–88.

    Article  Google Scholar 

  53. Pocock MJO, Hauffe HC, Searle JB. Dispersal in house mice. Biol J Lin Soc. 2005;84:565–83.

    Article  Google Scholar 

  54. Smith J, Hurst JL, Barnard CJ. Comparing behaviour in wild and laboratory strains of the house mouse: levels of comparison and functional inference. Behav Proc. 1994;32:79–86.

    Article  CAS  Google Scholar 

  55. Blanchard RJ, Hebert MA, Ferrari P, Palanza P, Figueira R, Blanchard DC, et al. Defensive behaviors in wild and laboratory (Swiss) mice: the mouse defense test battery. Physiol Behav. 1998;65:201–9.

    Article  CAS  PubMed  Google Scholar 

  56. Augustsson H. Ethoexperimental studies of behaviour in wild and laboratory mice: risk assessment, emotional reactivity and animal welfare. Uppsala: Dept. of Large Animal Clinical Sciences, Swedish Univ. of Agricultural Sciences; 2004.

  57. Tuttle AH, Philip VM, Chesler EJ, Mogil JS. Comparing phenotypic variation between inbred and outbred mice. Nat Methods. 2018;15:994–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Tuttle AH, Tansley S, Dossett K, Tohyama S, Khoutorsky A, Maldonado-Bouchard S, et al. Social propinquity in rodents as measured by tube cooccupancy differs between inbred and outbred genotypes. Proc Natl Acad Sci U S A. 2017;114:5515–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. 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.

    Article  Google Scholar 

  60. Crowcroft P. Mice all over. 1966.

    Google Scholar 

  61. Crowcroft P, Rowe FP. Social organisation and territorial behaviour in the wild house mouse (Mus musculus L.). Proc Zool Soc Lond. 1963;140:517–31.

    Article  Google Scholar 

  62. Gerlach G. Emigration mechanisms in fetal house mice - a laboratory investigation of the influence of social structure, population density, and aggression. Behav Ecol Sociobiol. 1996;39:159–70.

    Article  Google Scholar 

  63. Lidicker WZ. Social behaviour and density regulation in house mice living in large enclosures. J Anim Ecol. 1976;45:677–97.

    Article  Google Scholar 

  64. Luzynski KC, Nicolakis D, Marconi MA, Zala SM, Kwak J, Penn DJ. Pheromones that correlate with reproductive success in competitive conditions. Sci Rep. 2021;11:21970.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Meagher S, Penn DJ, Potts WK. Male–male competition magnifies inbreeding depression in wild house mice. PNAS. 2000;97:3324–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Noyes RF, Barrett GW, Taylor DH. Social structure of feral house mouse (Mus musculus L.) populations: effects of resource partitioning. Behav Ecol Sociobiol. 1982;10:157–63.

    Article  Google Scholar 

  67. Potts WK, Manning CJ, Wakeland EK. Mating patterns in seminatural populations of mice influenced by MHC genotype. Nature. 1991;352:619–21.

    Article  CAS  PubMed  Google Scholar 

  68. Sherborne AL, Thom MD, Paterson S, Jury F, Ollier WER, Stockley P, et al. The genetic basis of inbreeding avoidance in house mice. Curr Biol. 2007;17:2061–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Stockley P, Ramm SA, Sherborne AL, Thom MD, Paterson S, Hurst JL. Baculum morphology predicts reproductive success of male house mice under sexual selection. BMC Biol. 2013;11:1–6.

    Article  Google Scholar 

  70. Wolff RJ. Mating behaviour and female choice: their relation to social structure in wild caught House mice (Mus musculus) housed in a semi-natural environment. J Zool. 1985;207:43–51.

    Article  Google Scholar 

  71. Hurst JL. Behavioural variation in wild house mice Mus domesticus Rutty: a quantitative assessment of female social organization. Anim Behav. 1987;35:1846–57.

    Article  Google Scholar 

  72. Crowcroft P. Territoriality in wild house mice, Mus musculus L. J Mammal. 1955;36:299–301.

    Article  Google Scholar 

  73. Brown RZ. Social behavior, reproduction, and population changes in the house mouse (Mus musculus L.). Ecol Monogr. 1953;23:218–40.

    Article  Google Scholar 

  74. Rusu AS, Krackow S. Kin-preferential cooperation, dominance-dependent reproductive skew, and competition for mates in communally nesting female house mice. Behav Ecol Sociobiol. 2004;56:298–305.

    Article  Google Scholar 

  75. Harrison N, Lindholm AK, Dobay A, Halloran O, Manser A, König B. Female nursing partner choice in a population of wild house mice (Mus musculus domesticus). Front Zool. 2018;15:4.

    Article  PubMed  PubMed Central  Google Scholar 

  76. Cunningham CB, Ruff JS, Chase K, Potts WK, Carrier DR. Competitive ability in male house mice (Mus musculus): genetic influences. Behav Genet. 2013;43:151–60.

    Article  PubMed  PubMed Central  Google Scholar 

  77. Jo Manning C, Dewsbury DA, Wakeland EK, Potts WK. Communal nesting and communal nursing in house mice, Mus musculus domesticus. Anim Behav. 1995;50:741–51.

    Article  Google Scholar 

  78. Wang F, Kessels HW, Hu H. The mouse that roared: neural mechanisms of social hierarchy. Trends Neurosci. 2014. https://doi.org/10.1016/j.tins.2014.07.005.

    Article  PubMed  PubMed Central  Google Scholar 

  79. Varholick JA, Pontiggia A, Murphy E, Daniele V, Palme R, Voelkl B, et al. Social dominance hierarchy type and rank contribute to phenotypic variation within cages of laboratory mice. Sci Rep. 2019;9:13650.

    Article  PubMed  PubMed Central  Google Scholar 

  80. Lee W, Milewski TM, Dwortz MF, Young RL, Gaudet AD, Fonken LK, et al. Distinct immune and transcriptomic profiles in dominant versus subordinate males in mouse social hierarchies. Brain Behav Immun. 2022;103:130–44.

    Article  CAS  PubMed  Google Scholar 

  81. Palanza P, Della Seta D, Ferrari PF, Parmigiani S. Female competition in wild house mice depends upon timing of female/male settlement and kinship between females. Anim Behav. 2005;69:1259–71.

    Article  Google Scholar 

  82. Palanza P, Re L, Mainardi D, Brain PF, Parmigiani S. Male and female competitive strategies of wild house mice pairs (Mus musculus domesticus) confronted with intruders of different sex and age in artificial territories. Behaviour. 1996;133:863–82.

    Article  Google Scholar 

  83. Mackintosh JH. Territory formation by laboratory mice. Anim Behav. 1970;18 PART 1:177–83.

    Article  Google Scholar 

  84. Poole TB, Morgan HDR. Social and territorial behaviour of laboratory mice (Mus musculus L.) in small complex areas. Anim Behav. 1976;24:476–80.

    Article  Google Scholar 

  85. Morris JS, Ruff JS, Potts WK, Carrier DR. A disparity between locomotor economy and territory-holding ability in male house mice. J Exp Biol. 2017;220:2521–8.

    PubMed  PubMed Central  Google Scholar 

  86. Davis DE. The role of density in aggressive behaviour of house mice. Anim Behav. 1958;6:207–10.

    Article  Google Scholar 

  87. Fischer S, Duffield C, Davidson AJ, Bolton R, Hurst JL, Stockley P. Fitness costs of female competition linked to resource defense and relatedness of competitors. Am Nat. 2023;201:256–68.

    Article  PubMed  Google Scholar 

  88. Bruce HM. An exteroceptive block to pregnancy in the mouse. Nature. 1959;184:105–105.

    Article  CAS  PubMed  Google Scholar 

  89. Bruce HM. A block to pregnancy in the mouse caused by proximity of strange males. Reproduction. 1960;1:96–103.

    Article  CAS  Google Scholar 

  90. Yamazaki K, Beauchamp GK, Wysocki CJ, Bard J, Thomas L, Boyse EA. Recognition of H-2 types in relation to the blocking of pregnancy in mice. Science. 1983;221:186–8.

    Article  CAS  PubMed  Google Scholar 

  91. Zipple MN, Roberts EK, Alberts SC, Beehner JC. Male-mediated prenatal loss: functions and mechanisms. Evol Anthropol Issues News Rev. 2019;28:114–25.

    Article  Google Scholar 

  92. Miller CH, Reichard TM, Yang J, Carlson-Clarke B, Vogt CC, Warden MR, et al. Reproductive state switches the valence of male urinary pheromones in female mice. 2022:2022.08.22.504866.

  93. Guo M, Wu CF, Liu W, Yang JY, Chen D. Sex difference in psychological behavior changes induced by long-term social isolation in mice. Prog Neuropsychopharmacol Biol Psychiatry. 2004;28:115–21.

    Article  PubMed  Google Scholar 

  94. Koike H, Ibi D, Mizoguchi H, Nagai T, Nitta A, Takuma K, et al. Behavioral abnormality and pharmacologic response in social isolation-reared mice. Behav Brain Res. 2009;202:114–21.

    Article  CAS  PubMed  Google Scholar 

  95. Matthews GA, Nieh EH, Vander Weele CM, Halbert SA, Pradhan RV, Yosafat AS, et al. Dorsal raphe dopamine neurons represent the experience of social isolation. Cell. 2016;164(4):617–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  96. Nonogaki K, Nozue K, Oka Y. Social isolation affects the development of obesity and type 2 diabetes in mice. Endocrinology. 2007;148:4658–66.

    Article  CAS  PubMed  Google Scholar 

  97. Zelikowsky M, Hui M, Karigo T, Choe A, Yang B, Blanco MR, et al. The neuropeptide Tac2 controls a distributed brain state induced by chronic social isolation stress. Cell. 2018;173:1265-1279.e19.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  98. Chambers LK, Singleton GR, Krebs C. Movements and social organization of wild house mice (Mus domesticus) in the wheatlands of northwestern Victoria, Australia. J Mammal. 2000;81:59–69.

    Article  Google Scholar 

  99. Singleton GR. The social and genetic structure of a natural colony of house mice, Mus musculus, at Healesville WildlifeSanctuary. Aust J Zool. 1983;31:155–66.

    Article  Google Scholar 

  100. Barton RA, Byrne RW, Whiten A. Ecology, feeding competition and social structure in baboons. Behav Ecol Sociobiol. 1996;38:321–9.

    Article  Google Scholar 

  101. Faulkes CG, Bennett NC, Bruford MW, Obrien HP, Aguilar GH, Jarvis JUM. Ecological constraints drive social evolution in the African mole-rats. Proc R Soc Lond Ser B Biol Sci. 1997;264:1619–27.

    Article  CAS  Google Scholar 

  102. Lukas D, Clutton-Brock T. Climate and the distribution of cooperative breeding in mammals. R Soc Open Sci. 2017;4:160897.

  103. Elliser CR, Herzing DL. Social structure of Atlantic spotted dolphins, Stenella frontalis, following environmental disturbance and demographic changes. Mar Mamm Sci. 2014;30:329–47.

    Article  Google Scholar 

  104. Ostfeld RS, Lidicker WZ, Heske EJ. The relationship between habitat heterogeneity, space use, and demography in a population of California voles. Oikos. 1985;45:433–42.

    Article  Google Scholar 

  105. Palombit RA, Seyfarth RM, Cheney DL. The adaptive value of ‘friendships’ to female baboons: experimental and observational evidence. Anim Behav. 1997;54:599–614.

    Article  CAS  PubMed  Google Scholar 

  106. Evans JC, Liechti JI, Boatman B, König B. A natural catastrophic turnover event: individual sociality matters despite community resilience in wild house mice. Proc R Soc B Biol Sci. 2020;287:20192880.

    Article  Google Scholar 

  107. Lukas D, Clutton-Brock T. Cooperative breeding and monogamy in mammalian societies. Proc R Soc B Biol Sci. 2012;279:2151–6.

    Article  Google Scholar 

  108. Dobson FS, Way BM, Baudoin C. Spatial dynamics and the evolution of social monogamy in mammals. Behav Ecol. 2010;21:747–52.

    Article  Google Scholar 

  109. Perez-Barberia FJ, Shultz S, Dunbar RIM. Evidence for coevolution of sociality and relative brain size in three orders of mammals. Evolution. 2007;61:2811–21.

    Article  PubMed  Google Scholar 

  110. Schradin C. Intraspecific variation in social organization by genetic variation, developmental plasticity, social flexibility or entirely extrinsic factors. Philos Trans R Soc B Biol Sci. 2013;368:20120346.

    Article  Google Scholar 

  111. Schradin C, Hayes LD, Pillay N, Bertelsmeier C. The evolution of intraspecific variation in social organization. Ethology. 2018;124:527–36.

    Article  Google Scholar 

  112. Li SW, Zeliger O, Strahs L, Báez-Mendoza R, Johnson LM, McDonald Wojciechowski A, et al. Frontal neurons driving competitive behaviour and ecology of social groups. Nature. 2022;603:661–6.

    Article  CAS  PubMed  Google Scholar 

  113. de Chaumont F, Ey E, Torquet N, Lagache T, Dallongeville S, Imbert A, et al. Real-time analysis of the behaviour of groups of mice via a depth-sensing camera and machine learning. Nat Biomed Eng. 2019;3:930–42.

    Article  PubMed  Google Scholar 

  114. Clutton-Brock T. Social evolution in mammals. Science. 2021;373:eabc9699.

    Article  CAS  PubMed  Google Scholar 

  115. Macbeth AH, Lee H-J, Edds J, Young WS. Oxytocin and the oxytocin receptor underlie intrastrain, but not interstrain, social recognition. Genes Brain Behav. 2009;8:558–67.

    Article  CAS  PubMed  Google Scholar 

  116. Ferguson JN, Young LJ, Hearn EF, Matzuk MM, Insel TR, Winslow JT. Social amnesia in mice lacking the oxytocin gene. Nat Genet. 2000;25:284–8.

    Article  CAS  PubMed  Google Scholar 

  117. Oliva A, Fernández-Ruiz A, Leroy F, Siegelbaum SA. Hippocampal CA2 sharp-wave ripples reactivate and promote social memory. Nature. 2020;587:264–9.

    Article  CAS  PubMed  Google Scholar 

  118. Barnard CJ, Fitzsimons J. Kin recognition and mate choice in mice: fitness consequences of mating with kin. Anim Behav. 1989;38:35–40.

    Article  Google Scholar 

  119. Krackow S, Matuschak B. Mate choice for non-siblings in wild house mice: evidence from a choice test and a reproductive test. Ethology. 1991;88:99–108.

    Article  Google Scholar 

  120. König B, Lindholm AK, Macholan M, Baird SJ, Munclinger P, Pialek J. The complex social environment of female house mice (Mus domesticus). In: Cambridge studies in morphology and molecules: new paradigms in evolutionary bio. 2012. p. 114–34.

  121. Wolff JO. Why are female small mammals territorial? Oikos. 1993;68:364–70.

  122. Wolff JO. Social biology of rodents. Integr Zool. 2007;2:193–204.

    Article  PubMed  Google Scholar 

  123. Rowe FP, Redfern R. Aggressive behaviour in related and unrelated wild house mice (Mus musculus L.). Ann Appl Biol. 1969;64:425–31.

    Article  Google Scholar 

  124. Chovnick A, Yasukawa NJ, Monder H, Christian JJ. Female behavior in populations of mice in the presence and absence of male hierarchy. Aggressive Behav. 1987;13:367–75.

    Article  Google Scholar 

  125. Chalfin L, Dayan M, Levy DR, Austad SN, Miller RA, Iraqi FA, et al. Mapping ecologically relevant social behaviours by gene knockout in wild mice. Nat Commun. 2014;5:4569.

    Article  CAS  PubMed  Google Scholar 

  126. Lee W, Milewski TM, Dwortz MF, Young RL, Gaudet AD, Fonken LK, et al. Distinct inflammatory and transcriptomic profiles in dominant versus subordinate males in mouse social hierarchies. 2021:2021.09.04.458987.

  127. Tumulty JP. Dear enemy effect. In: Vonk J, Shackelford T, editors. Encyclopedia of animal cognition and behavior. Cham: Springer International Publishing; 2018. p. 1–4.

    Google Scholar 

  128. Miller CH, Hillock MF, Yang J, Carlson-Clarke B, Haxhillari K, Lee AY, et al. Dynamic changes to signal allocation rules in response to variable social environments in house mice. 2022:2022.01.28.478242.

  129. Filazzola A, Cahill JF Jr. Replication in field ecology: identifying challenges and proposing solutions. Methods Ecol Evol. 2021;12:1780–92.

    Article  Google Scholar 

  130. Hatchwell BJ. Replication in behavioural ecology: a comment on Ihle et al. Behav Ecol. 2017;28:360–360.

    Article  Google Scholar 

  131. Kelly CD. Replicating empirical research in behavioral ecology: how and why it should be done but rarely ever is. Q Rev Biol. 2006;81:221–36.

    Article  PubMed  Google Scholar 

  132. ASAB Ethical Committee/ABS Animal Care Committee. Guidelines for the ethical treatment of nonhuman animals in behavioural research and teaching. Anim Behav. 2023;195:I–XI.

    Article  Google Scholar 

  133. Southwick CH. Population characteristics of house mice living in english corn ricks: density relationships. Proc Zool Soc Lond. 1958;131:163–75.

    Article  Google Scholar 

  134. Evans JC, Lindholm AK, König B. Long-term overlap of social and genetic structure in free-ranging house mice reveals dynamic seasonal and group size effects. Curr Zool. 2021;67:59–69.

    Article  PubMed  Google Scholar 

  135. Yang H, Wang JR, Didion JP, Buus RJ, Bell TA, Welsh CE, et al. Subspecific origin and haplotype diversity in the laboratory mouse. Nat Genet. 2011;43:648–55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  136. Joan S. Chapter 1 - The Laboratory Mouse. In: Biology of the laboratory mouse. 1966.

  137. Nishioka Y. The origin of common laboratory mice. Genome. 1995;38:1–7.

    Article  CAS  PubMed  Google Scholar 

  138. Parker M, Mulder G. The history of Black 6 mice. Charles River; 2021. https://www.criver.com/eureka/history-black-6-mouse. Accessed 20 Dec 2023.

  139. Engber D. The trouble with Black-6. Slate Magazine; 2011. http://www.slate.com/articles/health_and_science/the_mouse_trap/2011/11/black_6_lab_mice_and_the_history_of_biomedical_research.html. Accessed 20 Dec 2023.

  140. Phifer-Rixey M, Bi K, Ferris KG, Sheehan MJ, Lin D, Mack KL, et al. The genomic basis of environmental adaptation in house mice. PLoS Genet. 2018;14:1–28.

    Article  Google Scholar 

  141. Dumont BL, Gatti DM, Ballinger MA, Lin D, Phifer-Rixey M, Sheehan MJ, et al. Into the wild: a novel wild-derived inbred strain resource expands the genomic and phenotypic diversity of laboratory mouse models. 2023:2023.09.21.558738.

  142. Anderson PK. Density, social structure, and nonsocial environment in house-mouse populations and the implications for regulation of numbers. Trans N Y Acad Sci. 1961;23:447–51.

    Article  CAS  PubMed  Google Scholar 

  143. Fraser H, Barnett A, Parker TH, Fidler F. The role of replication studies in ecology. Ecol Evol. 2020;10:5197–207.

    Article  PubMed  PubMed Central  Google Scholar 

  144. Berry RJ. The natural history of the house mouse. 1970.

    Google Scholar 

  145. Farine DR. Animal social network inference and permutations for ecologists in R using asnipe. Methods Ecol Evol. 2013;4:1187–94.

    Article  Google Scholar 

  146. Csardi G, Nepusz T. The Igraph software package for complex network research. InterJournal. 2005;Complex Systems:1695.

    Google Scholar 

  147. Aplin LM, Firth JA, Farine DR, Voelkl B, Crates RA, Culina A, et al. Consistent individual differences in the social phenotypes of wild great tits, Parus major. Anim Behav. 2015;108:117–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  148. Farine DR, Whitehead H. Constructing, conducting and interpreting animal social network analysis. J Anim Ecol. 2015;84:1144–63.

    Article  PubMed  PubMed Central  Google Scholar 

  149. Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. 2014.

    Google Scholar 

  150. Kuznetsova A, Brockhoff PB, Christensen RHB. lmerTest package: tests in linear mixed effects models. J Stat Softw. 2017;82:1–26.

    Article  Google Scholar 

  151. Lenth RV. Least-squares means: the R Package lsmeans. J Stat Softw. 2016;69:1–33.

    Article  Google Scholar 

  152. Ameijeiras-Alonso J, Crujeiras RM, Rodriguez-Casal A. multimode: an R package for mode assessment. J Stat Softw. 2021;97:1–32.

    Article  Google Scholar 

  153. Ameijeiras-Alonso J, Crujeiras RM, Rodríguez-Casal A. Mode testing, critical bandwidth and excess mass. TEST. 2019;28:900–19.

    Article  Google Scholar 

  154. Wickham H. Programming with ggplot2. In: Wickham H, editor. ggplot2: elegant graphics for data analysis. Cham: Springer International Publishing; 2016. p. 241–53.

    Chapter  Google Scholar 

Download references

Acknowledgements

We thank Melissa R. Warden for helpful discussions during the planning of the experiment. We thank Gary Olz and Russel L. Ligon for their assistance in preparing the field site.

Funding

Funding was provided by Cornell University (Neurobiology and Behavior Departmental Grant to CCV), the USDA (Hatch Grant, NYC-191428 to MJS), NIH Animal Models for the Social Dimensions of Health and Aging Research Network Pilot Grant (NIH/NIH R24 AG065172 to MJS), and the NIH (R35 GM138284 to AHM).

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization: CCV, MJS. Study design: CCV, MJS. Methodology design: CCV, MJS. Field work: CCV, DDS, CHM, MSC, AHM, MJS. Data curation and code: CCV. Data analysis: CCV, MNZ, SXH, MKA, AMG, LMM, MJS. Figure creation: CCV, MJS. Writing—original draft: CCV, MJS. Writing—review and editing: CCV, MNZ, AHM, MJS. Supervision: CCV, AHM, MJS. Funding acquisition: CCV, AHM, MJS.

Corresponding authors

Correspondence to Caleb C. Vogt or Michael J. Sheehan.

Ethics declarations

Ethics approval and consent to participate

All experimental procedures were approved by the Cornell University Institutional Animal Care and Use committee (IACUC: Protocol #2015-0060).

Consent for publication

All authors read and approved the final manuscript for publication.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1.

Field site setup and RFID duration bout window selection. (A) Satellite image of the field enclosures showing the position of the data sheds housing the computer for downloading RFID data from the central RFID sheds. (B) Graphical schematic of the Alpha and Bravo enclosures indicating the resource zone layouts. (C) RFID monitoring of the resource zones. Two storage tubs were nested with a RFID antenna placed between them, beneath the entrance tunnel to prevent mice from directly contacting the antenna and wire. (D) Histogram of the daily inter-zone travel times for all mice for all days. (E) Histogram of the within zone inter-RFID read intervals and the 139-s threshold capturing 99% of all inter-RFID read intervals within the same zone which was used to group RFID reads into resource zone visitation bouts (see the “Methods” section). (F) Schematic of the time window capture threshold grouping procedure to determine resource zone visit bouts and spatiotemporal association bouts using mock RFID data and the time window capture threshold (139 s). (G) Correlation of estimated duration spent in each zone and the number of visits to that zone for all sex and genotype categories.

Additional file 2: Figure S2.

Territory establishment in C57 and wild-derived male mice. (A) Schematic of the resource zone locations (colored boxes) within the field enclosures (2 x 4 grids) showing patterns of male resource zone usage (rows) for an example C57 (left) and WD (right) trial across 10 days of activity (columns). White boxes indicate resource zones that were not visited by the focal individual. (B) Percentage of the total time a mouse was observed across all zones spent in a mouse’s top occupied zone (resource zones rank ordered by mouse occupancy time). (C) Male-male social grouping events fell rapidly over time and were less frequent in WD trials compared to C57 trials. (D) Over time males tend to either gain priority access to resource zones (i.e., are territorial) or do not (non-territorial males). Y-axis shows the evolving daily priority access scores over 10 days of observation for males of both genotypes (see the “Methods” section for additional details on calculation of the daily PAS value). (E) C57 mice had longer male-female social grouping bout durations overall. For visualization purposes, the y-axis is cut off at 40 (n = 15,891 events shown out of 15,898 total events). (F) Daily zones visited were lower for high status males (Day 10 PAS > 0) than for low status males (Day 10 PAS < 0). Data are plotted as means ± SEM.

Additional file 3: Figure S3.

Distinct patterns of space utilization in C57 females. (A) C57 females extensively explore the available resource zones over the course of the trial compared to WD females. Example full trial data of female zone use in C57 (left) and WD (right) mice. Schematic of the resource zone locations (colored boxes) within the field enclosures (2 x 4 grids) showing patterns of zone usage for animals (rows) across 10 days of activity (columns). White boxes indicate resource zones that were not visited by the focal individual. (B) Daily priority access scores over 10 days of observation for female mice. (C) Distributions of Day 10 priority access scores for female mice are not multi-modal (excess mass test for unimodality from the multimode package), indicating decreased or inconsistent monopolization of resource zones amongst females. Higher scores indicate the extent to which a mouse maintained majority access over one or more resource zones relative to same-sex conspecific competitors (see the “Methods” section for details). (D) C57 female mice (n = 40) differed from C57 male (n = 40) and WD male (n = 29) and female (n = 30) mice in their estimated minimum distance travelled over the course of 10 days. (E) WD females spent more time in their most occupied zone than C57 females. (F) Female-female social grouping bout durations over time. For visualization purposes, the y-axis is cut off at 40 (n = 11,928 events shown out of 11,933 total events). (G) WD females spent nearly all of their female-female social time with cage mates after day 2, in contrast with C57 females who generally spent less than half of their female-female social time with cage mates. Data are plotted as means ± SEM.

Additional file 4: Figure S4.

Repeatable differences in C57 social structure and network level properties over time. (A) Percentage of total recorded observation time spent in opposite-sex associations is higher in C57 mice over the course of the trial. (B) Percentage association time spent in opposite sex associations is higher in males than females for both genotypes. (C) C57 mice met a majority of the available novel social partners by the final day of the trial, while WD mice did not. (D) Both C57 and WD females exhibited high page rank scores relative to males of either genotype, indicating that females serve as major social connections through the network in both genotypes. Data are plotted as means ± SEM.

Additional file 5: Table S1.

Statistical outputs. Full statistical outputs for all analyses and figures.

Additional file 6: Table S2.

Metadata and trial information. Individual metadata information for all mice (sheet 1) as well as trial specific metadata information (sheet 2).

Additional file 7: Video S1. Sample behavioral videos within resource zones in the field. Video demonstrating sample social and foraging related behaviors in both C57 and WD animals within resource zones.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Vogt, C.C., Zipple, M.N., Sprockett, D.D. et al. Female behavior drives the formation of distinct social structures in C57BL/6J versus wild-derived outbred mice in field enclosures. BMC Biol 22, 35 (2024). https://doi.org/10.1186/s12915-024-01809-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12915-024-01809-0

Keywords