Improved behavioral indices of welfare in continuous compared to intermittent pair‐housing in adult female rhesus macaques (Macaca mulatta)

Limiting opportunities for captive nonhuman primates (NHPs) to express species‐specific social behaviors may disrupt the adaptive drive for social companionship and may lead to increases in coping behaviors and inactivity. While captive NHPs show improved welfare when moving to pair‐housing from single‐housing, the impact of daily separation of pair‐mates, as is implemented in intermittent pair‐housing, is not fully understood. We compared behavioral indices of welfare exhibited by adult female rhesus macaques (Macaca mulatta) in two conditions: (1) intermittent pair‐housing, involving daily overnight separation of pair‐mates, and (2) continuous pair‐housing, involving little separation of pair‐mates. A within‐subjects study design tested two groups of females experiencing both pairing conditions in an alternate order, switching either from continuous to intermittent pair‐housing, or from intermittent to continuous pair‐housing. Behavioral observations, recording activity state, self‐directed, abnormal, and social behaviors, were conducted at midday when all females were paired, and in the afternoon when intermittent pairs were separated. Females exhibited higher levels of inactivity and self‐directed behavior when separated due to intermittent pair‐housing in comparison to continuous pair‐housing. In addition, intermittently paired females showed higher levels of grooming and other types of affiliation when paired, than during the same time frame when they were continuously paired. These results suggest that females in the continuous presence of a social partner experience improved levels of activity and do not need to elevate levels of behavioral coping mechanisms (e.g., self‐scratching, increased affiliation) as they receive the benefits associated with social companionship consistently throughout the day. Overall, this study provides the first evidence that continuous pair‐housing affords better welfare than intermittent pair‐housing in adult female rhesus macaques. Pair‐housing options, such as continuous pairing, that reduce reliance on behavioral coping mechanisms and promote adaptive social behavior throughout the entirety of the day should be prioritized over husbandry care scheduled for convenience.


| INTRODUCTION
Social interactions and relationships are fundamentally important to primates, providing individuals with a range of fitness and health benefits (Majolo & Huang, 2018;Ostner & Schülke, 2018;Silk, 2007;Snyder-Mackler et al., 2020). For captive nonhuman primates (NHPs), social housing is the most effective form of environmental enrichment (Lutz & Novak, 2005), improving physiological, behavioral, and psychological measures of welfare (Olsson & Westlund, 2007). Indoor-housed NHPs, particularly macaque species, are often pairhoused where two individuals inhabit connected adjacent cages (Baker & Crockett et al., 2012). Ideally, pair-mates would be continuously pair-housed, having complete physical access to their pairmate with very little separation, but this presents a problem for some research objectives (e.g., urine sample collection). Intermittent pairhousing, involving the temporary daily or weekly separations of pairmates that can last 12 or more hours (often occurring overnight), has been developed to simultaneously allow some social contact and facilitate research protocols.
Temporary separations due to intermittent pair-housing, however, may still negatively impact NHP welfare due to the reduced opportunity to receive the benefits social contact provides. In a study of adult female rhesus macaques (Macaca mulatta), higher urinary cortisol concentration (a measure of stress physiology) was related to the interaction of intermittent pair-housing and a poor pair relationship quality (Hannibal et al., 2018). By contrast, in the same species, no differences in activity, abnormal, or anxiety-related behaviors were seen between continuous and intermittent pair-housing conditions while pair-mates were together (Baker & Bloomsmith et al., 2012;Baker & Crockett et al., 2012); however, these behavioral measures were not assessed while animals were separated.
Pair-housed adult female rhesus macaques strongly prefer to be in close proximity with their pair-mate, particularly overnight (Eaton, Kelley, Axthelm, Iliff-Sizemore, & Shiigi, 1994). Co-sleeping behavior is suggested to have evolved in response to ecological pressures such as predation  and low ambient temperature (Gilbert et al., 2009). While captivity obviously removes these pressures, NHPs instinctively seek out social partners to sleep with during the night (e.g., Mochida & Nishikawa, 2014) and the prevention of such behavior could impair welfare. Therefore, evaluating measures of welfare during the time pair-mates spend apart is critical, as intermittently paired NHPs are separated for up to three-quarters of the day, including overnight.
Here, we investigate behavioral indices of welfare in intermittent and continuous pair-housed adult female rhesus macaques (M. mulatta), quantifying time spent inactive, exhibiting abnormal behavior, exhibiting affiliative behavior, and bouts of SDB. We predicted that inactivity, abnormal behavior, and SDB would be highest during the intermittent pair-housing phase, when pair-mates were separated. Additionally, affiliative behaviors were predicted to be higher during the intermittent pair-housing phase.

| Study subjects
Subjects were 24 adult female rhesus macaques, born and reared at the CNPRC. Females were raised in either 0.2-hectare outdoor enclosures (i.e., field cage) each containing up to 180 monkeys of all age and sex classes (N = 17 individuals), or in ∼43.7-m 2 outdoor enclosures (i.e., corn crib) each containing up to 30 monkeys of all age and sex classes (N = 7 individuals). The age of the subjects ranged from 4.9 to 10.9 years (mean ± standard deviation: 6.5 ± 1.7 years). Females weighed 4.7 to 8.7 kg (6.5 ± 1.2 kg) and were not pregnant. Before study participation, animals were housed indoors for at least 4 months (18.2 ± 18.1 months) and paired together for more than 3 months (9.9 ± 6.1 months) in their initial pairing condition. Study subjects were housed with one other female either initially as a continuous (N = 7 pairs, 12 study animals) or intermittent (N = 6 pairs, 12 study animals) pair; two females were paired with pair-mates not included in the study as they did not meet our subject selection criteria for rearing history (i.e., mother-reared in an outdoor social group during the first year of life). Two of six pairs that were initially intermittently paired had a previous relationship with each other (e.g., outdoor cage mates); five of seven pairs that were initially continuously paired had previous relationships, where four of these pairs were from the same matriline as their pair-mate.

| Housing
The CNPRC and other primate centers throughout the United States utilize intermittent and continuous pair-housing to provide indoorhoused NHPs with contact socialization with conspecifics. During the current study, intermittent pair-housed rhesus macaques were separated before afternoon feeding (about 14:00) and reunited the next day after morning feeding (by 08:00); these animals thus experienced about 6 h of physical contact per day with their pair-mate.
In contrast, continuously pair-housed animals were paired for over 18 h daily and experienced separation only during feeding times to eliminate the likelihood of food aggression and facilitate urine collection for another project (Hannibal et al., 2018). Overnight separation for subjects in the intermittent contact phase, while the continuous pairs remained together overnight, allowed the behavioral effects of this difference to be tested.
Each subject was housed in a stainless-steel cage (0.80 m high, floor space of 0.4 m 2 ). These cages sat on a two-tiered rack so that two pairs of cages were stacked one on top of the other. The cages of paired subjects were adjoined by a square opening (0.3 × 0.3 m) that was equipped with a sliding solid steel partition which was used to separate pair-mates when necessary and prevented physical and direct visual contact. Each cage was sanitized daily and was equipped with a metal reflective mirror, chew toy, forage board, and a metal perch. Animals were fed nutritionally complete biscuits and a forage mixture (i.e., rice, split peas, and oats) twice a day by CNPRC animal care staff, before 08:00 and after 14:00. Fresh water was available ad libitum. Additional facility enrichment (e.g., puzzle tubes filled with fruit, coconuts, novel forage, and alfalfa cubes) was provided once every 2 weeks. Animals experienced a 12-h light:dark cycle, with light onset at 06:00. Room temperature was maintained between 18 and 29°C and relative humidity between 30% and 70%. Subjects were housed in rooms containing 25 to 77 animals (48.5 ± 20.6 animals) at the beginning of the study period with visual, auditory, and olfactory access to other monkeys from a variety of age/sex classes.

| Experimental study design
Subjects were studied for 5 weeks, over which they were observed in their initial pairing condition at the beginning of the study for 2 weeks and then in their alternative pairing condition for three weeks ( Figure 1). Changes to pair-housing condition occurred due to the experimental manipulation of the present study as females were not enrolled in a concurrent project. Each pair of females was classified into one of two experimental groups: females in the CI experimental group were initially continuously (C) housed and then intermittently (I) housed. Conversely, IC experimental group animals began in intermittent pair-housing and then were switched to continuous pair-housing ( Figure 1; Table 1). Pairs were split into two cohorts of 12 subjects each, balanced by the experimental group, to accommodate data collection for all animals. Cohort 1 was studied from March 23 to April 24, 2015, and Cohort 2 was studied from CASSIDY ET AL.

| 3 of 14
April 27 to May 29, 2015. Behavior during the initial pairing condition was observed for each pair on weekdays for 9-10 data collection days and for the alternative pairing condition for 14-15 days.

| Data collection
Behavioral data were collected by Lauren C. Cassidy during observations of both animals in a pair during the two different pairhousing conditions. Each pair was observed for 8-min behavioral observations that took place four times a day during two observation periods, twice around midday (between 11:15 to 13:45; i.e., midday observation period), and twice in the late afternoon (between 15:30 and 18:00; i.e., afternoon observation period) after intermittent pairs were separated for the day. During the initial pairing condition, 3.6-4.3 h (4.0 ± 0.3 h) per animal of behavioral data were collected, and during the alternative pairing condition 7.3-8.0 h (7.7 ± 0.3 h) of data were collected; in total, 303.6 h of behavioral data were collected from all animals observed. Pair observation order was randomized to account for variation in behavior depending on time of F I G U R E 1 Experimental design. Each female in the study experienced intermittent (IP) and continuous (CP) pair-housing throughout the study. Females were assigned to an experimental group (IC, CI) depending on which pair-housing condition they began in T A B L E 1 Behavioral and model variable categories analyzed, with subcategories where applicable.

Category
Subcategory Definition

Model variable definitions
Predictor variables Current condition (factor) Current pair-housing condition (continuous or intermittent) Observation period (factor) Time of day observation occurred (midday or afternoon) Experimental group (factor) Pair began as intermittent and then experimentally changed to continuous (IC), or began as continuous and then experimentally changed to intermittent (CI)

Control variables Weight (covariate)
Weight in kilograms (kg) determined at most recent health assessment before the study.
Age (covariate) Age in years at the beginning of the study.

Dominant (factor)
Whether the animal received the greatest proportion of status signaling behaviors (move away, turn away, silent bared teeth) displayed between pair-mates over the entire study (yes or no).

Time spent indoors (covariate)
Total time in months that the subject was living in indoor housing before the study.

Previous relationship (factor)
Whether the pair-mates were familiar with each other before pairing (yes or no).

Duration paired (covariate)
Total time in months that the subject was living with current pair-mate before the study.
Note: The sampling method for each behavior (1/0, event) and type of model variable (factor, covariate) is indicated in parentheses following the category or subcategory where relevant. Definitions were derived from Hannibal et al., (2018).
day and to minimize the potential for disruptions due to CNPRC daily activity.
Behavioral data were recorded using fields and forms created Behaviors of interest (Table 1) were coded using either one-zero sampling (1/0) or all occurrences (i.e., event) sampling during a 20-s sample interval (Altmann, 1974). Proximity was coded concurrently with other behaviors. Other behaviors (e.g., locomotion, foraging) relating to activity budget were recorded, but not included in analyses. All behaviors except for SDB were recorded using 1/0 sampling. All bouts of SDB were recorded using event sampling and a new entry was made if separate bouts of the behavior occurred within the same interval. A break in a bout was defined if three or more seconds passed without the behavior occurring (Gilbert & Baker, 2011). Disruptions were recorded if the behavior of the focal subjects was interrupted due to noise outside of the room or CNPRC personnel entering the room.

| Data analyses
One pair from the IC experimental group required temporary separation due to injuries received due to an intrapair conflict part way through the study; therefore, the behavioral data of To assess whether the housing condition affected NHP welfare, we fit four generalized linear mixed models (GLMMs) with inactivity, SDB, abnormal behavior, and affiliation as the response variables with a negative binomial error structure using maximum likelihood (Bolker et al., 2009 Analyses were conducted in R (version 3.6.0; R Core Team, 2019), and GLMMs were fit using the "glmmTMB" package with the family argument set to "nbinom1" (version 0.2.3; Brooks et al., 2017).
Before fitting the models, we checked the distributions of control covariates and log-transformed age, duration paired, and inverse transformed time spent indoors to obtain more normal distributions.
Covariates were z-transformed to a mean of zero and a standard deviation of one to provide more comparable estimates and facilitate easier interpretation regarding interactions when scaled (Aiken, West, & Reno, 1991;Schielzeth, 2010).
Model stability was assessed by excluding subjects one at a time and comparing the subset model estimates with those from the full dataset. Using the "vif" function within the "car" package (version 3.0-2; (Fox & Weisberg, 2018)

| Inactivity
Inactivity occurred on average in 43.83% of intervals subjects were observed for. The model comparison between our full and null models for inactivity behavior was significant (LRT: χ 2 = 27.36, df = 5, p < .001). The model included significant interactions between current pair-housing condition and experimental group (LRT: χ 2 = 6.00, df = 1, p = .014), and current pairing housing condition and observation period (LRT: χ 2 = 11.27, df = 1, p = .001; Table 2). Model estimates indicate that inactivity levels were relatively stable across both observation periods when animals were continuously paired, whereas intermittently paired subjects were inactive more often in the afternoon when they were separated from their pair-mate than at midday when they were paired (Figure 2a). A comparison of model estimates across pair-housing conditions within each observation period indicates that inactivity was marginally lower for intermittent pairs than continuous pairs at midday but substantially higher for the former over the latter during the afternoon observation period when intermittently paired subjects were separated from their pair-mate ( Figure 2a). Additionally, model estimates for animals in the CI experimental group weakly indicate that inactivity was lower when they were intermittently paired than when they were continuously paired; the IC experimental group exhibited comparatively higher levels of inactivity when they were intermittently paired than the CI experimental group and when they were continuously paired ( Figure 2b).

| Self-directed behavior
Bouts of SDB occurred on average 4.5 times per minute of observation. The full model for total bouts of SDB was significantly different than its null model (LRT: χ 2 = 16.93, df = 5, p = .013). The interaction of the current pair-housing condition and observation period was significant (LRT: χ 2 = 7.45, df = 5, p = .006; Table 3). Model estimates indicate that bouts of SDB exhibited by continuously paired study subjects stayed relatively stable over observation periods, whereas bouts of SDBs exhibited when subjects were intermittently paired were lower at midday (when subjects had physical access to one another) than in the afternoon when subjects were separated from one another (Figure 3a). In comparison to continuous pair-housing levels, SDB was lower during the midday observation period, but higher in the afternoon when pair-mate separation occurred ( Figure 3a). Additionally, the interaction of current pairhousing condition and experimental group was significant (LRT: χ 2 = 4.26, df = 5, p = .039; Table 3). Animals in the CI experimental group exhibited more bouts of SDB when intermittently paired than when continuously paired, whereas animals in the IC experimental group exhibited less SDB when intermittently paired than when continuously paired (Figure 3b).

| Abnormal behavior
Abnormal behavior occurred on average in 3.47% of intervals subjects were observed for. The full-null model comparison for abnormal behavior was not significant (LRT: χ 2 = 5.11, df = 5, p = .402), indicating that the interactions or main effects of current pair-housing condition with the experimental group, and current condition with observation period did not explain the occurrence of abnormal behavior in our study population.

| DISCUSSION
Our study provides the first evidence that continuous pair-housing confers behavioral benefits over intermittent pair-housing.  Levels of afternoon inactivity were reduced when females had access to their pair-mate (continuous pair-housing condition) compared to when pair-mates were separated (intermittent pair-housing condition). Furthermore, females exhibited SDB less frequently during afternoon observations when they were continuously paired than when they were intermittently paired. Decreases in these measures are widely accepted as indicators of an improved behavioral repertoire in captive NHPs (Fureix & Meagher, 2015;Leeds & Lukas, 2019;Maestripieri et al., 1992).
Other pair-housing studies of female rhesus macaques found that levels of inactivity and bouts of SDB decreased when subjects in single-housing and protected contact switched to either intermittent or continuous pair-housing (Baker et al., 2014;Baker & Bloomsmith et al., 2012). However, they found no differences in levels of inactivity and SDB between intermittent and continuous housing conditions, leading to the conclusion that temporary separations due to intermittent pair-housing did not impede the overall benefits reaped from limited contact with a pair-mate (Baker et al., 2014).
It is important to note that Baker et al. (2014) did not assess behavioral measures of well-being while pair-mates were separated. The presence of conspecifics in captive environments allows captive NHPs to exhibit a wide range of species-typical patterns of behavior (Novak & Suomi, 1988). Physical separation, even if only temporary, completely removes the opportunity for tactile, species-specific behaviors (Baker & Bloomsmith et al., 2012;Novak et al., 1988), particularly affiliation and co-sleep which are key adaptive behaviors for group-living NHPs (Anderson, 2000;Dunbar, 1991). Preventing such behaviors from occurring may consequently result in elevated levels of inactivity or SDB, as seen in this study. Furthermore, increased levels of inactivity may have detrimental physiological (e.g., weight gain; Katzmarzyk, 2010) and psychological outcomes (e.g., boredom, lack of cognitive stimulation; Meagher, 2019). A pairhousing option that consistently minimizes abnormal behavior, inactivity, and SDB, and that promotes positive adaptive behaviors throughout the entirety of the day, would, therefore, provide superior well-being for laboratory NHPs.  Interestingly, the levels of grooming and other types of affiliative behaviors measured in this study were elevated when females were intermittently paired, whereas proximity and coalition behaviors were not related to the current pair-housing condition or experimental group. This result is similar to the trend Baker et al. (2014) found in their study, where adult female rhesus macaques tended to socialize more often in intermittent contact than when in continuous contact with their pair-mate. There could be several reasons why our females elevated their affiliative behavior when in intermittent contact. First, females in intermittent contact may increase rates of affiliative behavior within the window of time that they are paired, to mitigate daily stressors induced by a captive environment. Furthermore, it seems that behaviors involving physiological touch (i.e., grooming, other types of affiliation)-and not those such as proximity and coalitionary behaviors-may be the salient means for adult female rhesus to cope with the loss of physical access to their pair-mate. The social buffering hypothesis proposes that social partners temper behavioral and physiological stress responses (Balasubramaniam, Beisner, Vandeleest, Atwill, & McCowan, 2016;Cohen & Wills, 1985;Hennessy, Kaiser, & Sachser, 2009;Kikusui, Winslow, & Mori, 2006 Gilbert & Baker, 2011;husbandry: Boccia, Laudenslager, & Reite, 1995). The "tend-and-befriend" hypothesis further proposes that when levels of affiliation are inadequate, biological signaling systems (i.e., release of oxytocin) stimulate purposeful social behavior with conspecifics (Taylor et al., 2000). If social partners are supportive and comforting, these stress responses decline (Taylor, 2006).
The higher rates of affiliative behaviors exhibited by females in intermittent pair-housing suggest that these females are actively trying to cope with the daily temporary separation; because continuous housing conditions provide the constant benefit of social companionship, elevated levels of these behaviors during the middle of the day may be unnecessary. This interpretation is also supported by the results found for individual levels of inactivity and SDB during the midday observation period, where these behaviors were slightly lower than levels exhibited at midday when they were continuously housed.
The second explanation for this difference in affiliative behavior between pair-housing conditions is that females may have redistributed the time they spent engaging in affiliative behavior throughout the day when in continuous pair-housing. This redistribution could occur due to the reduced pressure to affiliate during a tight time window. Furthermore, continuous pair-housing allows pair-mates to huddle overnight, during which they can receive the benefits of physical contact without needing to make up a deficit the following day.
Alternatively, lower levels of active affiliative behavior during the continuous pair-housing phase may occur to minimize the length of physical contact bouts with a pair-mate who may not be very compatible. Although we did not explore measures of pair compatibility in this study (e.g., Hannibal et al., 2018), this alternative explanation is unlikely, due to the thorough pairing process conducted by trained staff at the CNPRC at the time of our study. Despite having met CNPRC's compatibility criteria, however, one IC pair did inflict minor wounds (scratches to face and arms) to each other during their continuous pair-housing phase, were consequently temporarily separated, and therefore not included in analyses. The change in pair-housing condition from intermittent to continuous may have triggered this incident to occur between these two pairmates due to a sudden increase in access time. The implications of changes to pair-housing condition are not very well understood and require longer term studies (Hannibal et al., 2017).
Although we did not find differences between experimental groups for affiliative behavior, we did find opposing differences in inactivity and SDB between the groups. Over both observation periods, females in the CI experimental group had slightly higher levels of SDB and lower inactivity when intermittently paired versus than when they were continuously paired, whereas IC females exhibited lower SDB and higher inactivity in the former pair-housing condition versus the latter. One distinct group difference is that most females assigned to the CI experimental group had pair-mates they previously knew from a social group (including kin and non-kin; five of seven pairs previously knew each other), whereas most IC females did not previously know each other (one of five pairs previously knew each other). We did account for previous relationship status as a control variable in our affiliation models and these analyses were explored at the pair level, finding no effect. However, this variable may partially explain these experimental group effects seen for inactivity and SDB. At the time of the study, CNPRC pair-housing staff used this knowledge of animals' previous relationships to match pairs as these individuals likely had previously established dominance relationships that would likely continue and therefore expedite the pair-housing process. Consequently, CI paired females may be more compatibly matched overall and may have been more agitated with the change in pair-housing condition (and reduced access to their pair-mate) as they went from being continuous pair-housed to intermittently pair-housed. Conversely, the IC experimental group experienced the reverse pair-housing change, which improved their overall activity but elevated SDB behavior, suggesting these females needed to boost these behaviors to cope with the new pair-housing arrangement.
Another explanation for experimental group effects could be that our results may have captured short-term adjustments and potential transient effects to changes in pair-housing (Baker et al., 2014). Although these effects appear to be moderate, researchers should consider incorporating a period of habituation into research plans as the impacts of pair-housing changes are relatively unknown (Capitanio, Kyes, & Fairbanks, 2006;Hannibal et al., 2017). Therefore, further long-term exploration of behavioral, as well as physiological, responses to intermittent and continuous pair-housing is necessary to understand if intermittent contact is a less favorable housing arrangement. Factors intrinsic to the individual (e.g., temperament, developmental history, early experience), the non-social (e.g., feeding regime), and social environment (e.g., who they are paired with) may also shape an individual's response (Boccia et al., 1995;Coleman, 2012), but exploring these was not within the scope our study. These numerous factors highlight the importance of assessing welfare in NHPs at the individual level. Nonetheless, understanding which captive NHP management practices, such as type of pair-housing, enhance welfare overall, is a vital first step forward.
The occurrence of abnormal behavior was not related to current pair-housing condition, observation period, or experimental group; however, this may be because these behaviors occurred infrequently.
Likely the low prevalence of abnormal behavior in our study can be attributed to our animal selection criteria (mother-reared in outdoor social groups). Previous studies have found that one of the main predictors of abnormal behavior development is restricted rearing background, such as nursery-rearing (Gottlieb, Capitanio, & McCowan, 2013;Novak et al., 2006;Rommeck, Anderson, Heagerty, Cameron, & McCowan, 2009;Vandeleest, McCowan, & Capitanio, 2011). None of our study subjects were nursery-reared, thus plausibly accounting for the low frequency of abnormal behavior observed. Recording these behaviors in the context of housing management practices, however, is critical to evaluate, particularly for vulnerable NHP populations, as these behaviors are known to occur due to experimental procedures (e.g., blood draws) and other housing management practices (i.e., cage location), and could develop into deleterious self-injurious behaviors in some cases (Gottlieb et al., 2013;Gottlieb, Maier, & Coleman, 2015).
Differences of the model estimates of SDB and other types of affiliative behavior were seemingly small (Figures 3 and 4), but statistically meaningful. Average SDB bouts increased by 69% from midday to afternoon observations for intermittently paired females and the average number of intervals other affiliation occurred increased by 44% during intermittent pair-housing in comparison to continuous pair-housing. From a practical standpoint, subtle shifts of these behaviors alone may not provide strong evidence that one pairhousing condition is superior. However, we assessed these behaviors alongside other behavioral indices of welfare and found the results largely in agreement that continuous housing confers greater behavioral welfare benefits over intermittent housing.
In the laboratory setting, intermittent pair-housing remains an alternative to single-housing when some separation is required as it offers NHPs some social contact with conspecifics while facilitating laboratory research procedures. However, modifications for carrying out common experimental procedures should be strongly considered before defaulting to intermittent pair-housing. For example, health monitoring for diarrhea or fecal sample collection can be conducted by offering pair-mates food items with differently colored dyes to distinguish feces. Macaques with a more limited social history, particularly while young, may have different outcomes under intermittent pairing. Overall, social-housing practices should prioritize the biological and psychological needs of the animals and not merely be determined by the ease with which they can be implemented by husbandry staff such as scheduling activities for convenience (Brando & Buchanan-Smith, 2018).
The feasibility of using integrated social housing systems to simultaneously maximize welfare and facilitate research should be assessed in future pair-housing studies. Both welfare and research objectives may be optimized by utilizing protected pair-housing in place of complete temporary separation. This method would, for example, prevent feeding aggression and facilitate the collection of separate biological specimens, such as urine, from individuals, while maximizing the social benefits of protected pair-housing during these times.

| CONCLUSIONS
Monitoring changes in a captive animal's behavioral repertoire provides a valuable means for assessing welfare noninvasively (Honess, Johnson, & Wolfensohn, 2004). One of the main goals of captive animal welfare research is to continually refine the management practices implemented in animal care, especially in the biomedical research environment. Refining husbandry practices not only improve animal welfare but also improve the scientific integrity of the research being carried out (Jennings & Prescott, 2009;Poole, 1997).
The present study provides evidence that continuous pair-housing confers welfare benefits over intermittent pair-housing for adult female rhesus macaques. These results enhance the current literature on pair-housing and will help guide further efforts to improve NHP welfare in the biomedical research setting.

DATA AVAILABILITY STATEMENT
The data that support the findings of this study are available from the corresponding author upon reasonable request.