The laboratory mouse is commonly used to study aging and test putative aging interventions (Ackert-Bicknell et al., 2015; Yuan et al., 2011). While many studies have used lifespan extension as an endpoint, lifespan in mice is not necessarily the best criterion for assessing the efficacy of aging interventions or for evaluating their potential translation to humans; measuring health is an important additional source of information (Hansen and Kennedy, 2016; Fischer et al., 2016). Health is a multi-parameter, organism-level state, and as multiple aspects of health decline with age, a sophisticated assessment of aging and putative aging interventions should evaluate as many organism-level aspects of health as possible (Freund, 2019).

Recognizing this, the field has developed assays to evaluate different aspects of animal behavior, physiology, and function. Particular emphasis has been given to assays designed to measure cognitive, metabolic, and neuromuscular metrics, as well as body size and composition (Huffman et al., 2016; Bellantuono et al., 2020; Sukoff Rizzo et al., 2018; Richardson et al., 2016). However, many of these assays are labor- and time-intensive, thus the number that any one study can include is limited. Further, these assays often utilize a challenge, for example, performance on a difficult task, and due to the stress this places on the animals as well as changes in performance due to training, challenge-based assays have limited repeatability, reducing their utility for longitudinal analysis (Huffman et al., 2016). Additionally, investigators must pay attention to several confounding variables and potential artifacts such as motivation to execute the task, order of assays, and time allowed for recovery between assays (Crabbe et al., 1999; Chesler et al., 2002b; McIlwain et al., 2001). Lastly, many assays are highly operator-dependent, adversely impacting reproducibility (Brown et al., 2018; Voelkl et al., 2020). Thus, although the use of multiple challenge-based assays provides insight into fitness and health, this phenotyping strategy has serious limitations.

To develop an alternative approach, we utilized commercially available, automated monitoring cages that simultaneously and continuously record physiological and behavioral parameters (voluntary animal position, respiration, weight, food and water intake, and wheel activity) (Brown et al., 2018). This type of hardware is not new, but the data from these systems has been underutilized. Each of the sensors in the cage provides dense time series data, with time windows ranging from <1 s (photobeam breaks) to 3 min (gas measurements). The most common analytical approach is to average this data into 12 hr time bins, that is, dark and light phase. Although useful for detecting large effects, more granular analyses provide more information. We hypothesized that, by generating hundreds of features from these raw data streams, insights into behavior and physiology could be uncovered via tools such as network modeling, which have previously been restricted to lower levels of biological organization. Additionally, because these measurements are acquired non-invasively during the course of daily living, observation can be continued for multiple days and repeated multiple times throughout an animal’s life. Observation across multiple days provides important information about circadian stability and captures animal behavior during both the light and dark phases; this is difficult to accomplish using any other method (Brown et al., 2018). Repeated, longitudinal phenotyping provides major advantages for aging studies by allowing for baseline normalization (which increases statistical power), direct measurement of rates of change, predictive modeling of future outcomes, and parsing of survivorship bias (Bellantuono et al., 2020; Zhang and Pincus, 2016). Further, automated phenotyping removes confounders such as operator-induced variability, increasing data robustness (Chesler et al., 2002a). Lastly, the system can be readily scaled without compromising data continuity and without the need for training highly specialized staff.

We used this platform to assess aging in the diversity outbred (DO) population, a heterogeneous stock derived from the intercrossing of eight inbred founder strains, in which each mouse is genetically unique (Churchill et al., 2012). Widespread reliance on a small number of inbred mouse strains for preclinical studies raises questions about generalizability. Phenotype and intervention response are highly dependent on genetic background, meaning that the findings in one mouse strain may not generalize to other mouse strains, let alone to humans (Voelkl et al., 2020; Sittig et al., 2016; Mandillo et al., 2008; Kafkafi et al., 2018; Tuttle et al., 2018). This is a concern across many disease areas, including aging (Liao et al., 2010). In recent years, mouse resources with high genetic diversity and phenotypic variation have become available, and this diversity provides important advantages. First, it allows for genetic analysis of phenotypes, which can be used to validate phenotypic data based on measurable heritability, to better understand the correlation between traits, and, with sufficient sample size, to genetically map traits (Svenson et al., 2012; Churchill et al., 2012). Second, genetic diversity introduces inter-individual variation, which is useful when building networks based on phenotypic correlations. Third, and perhaps most importantly, genetic diversity introduces deliberate biological variation in order to increase the external validity of study findings (Voelkl et al., 2020).

We generated ~60 years of data across 415 female DO mice, from an age range of 6 weeks to 40.5 months (3.4 years). We utilized a staggered enrollment design in which mice were enrolled at different ages, between 3 and 24 months, to ensure that the dataset contains hundreds of runs for all age groups up to 30+ months of age, providing consistent statistical power across an average mouse lifespan. To our knowledge, this represents the highest-resolution assessment of murine physiological aging to date: we examined animals monthly, whereas previous assessments of physiological aging sampled animals in ≥6-month intervals (Ladiges et al., 2017; Fischer et al., 2016; Petr et al., 2021). High-resolution phenotyping was manageable because our passive monitoring platform required only 20 hr of hands-on time per week from a single operator; as we and others have experienced, challenge-based phenotyping pipelines require so much hands-on time and often place such stress on the animals that they cannot be repeated frequently.

The platform detected hundreds of aging-related changes starting early in life, demonstrating that this passive monitoring approach has the sensitivity to quantify aging even in a highly variable strain background. Network analysis demonstrated that different behavioral and physiological clusters exhibited distinct aging dynamics. By employing analytical methods to measure overall network connectivity, we quantified resilience, an emergent property, and found that it exhibited an accelerating decline that was distinct from any individual phenotype. Finally, we develop a new method for determining biological age: a Combined Age and Survival Prediction of Aging Rate (CASPAR) model that is trained to simultaneously predict both chronological age and survival time. In summary, we demonstrate that automated phenotyping of murine behavior and physiology, when combined with high-dimensional analytics, identifies aspects of aging and resilience that have not been previously described and provides important advantages over existing approaches in terms of throughput, temporal resolution, and physiological scope. Our study builds upon pre-existing literature from other model organisms, particularly nematodes, demonstrating that passive, automated monitoring can be used to quantify multi-dimensional, organism-level aging (Zhang et al., 2016; Le et al., 2020; Martineau et al., 2020). This combined hardware/software platform can be replicated in any vivarium and used for in-depth characterization of behavior and physiology as well as intervention testing.

Many assessments of aging utilize a challenge, for example, a maze for cognitive function or a rotating rod for coordination, on the premise that measuring the limits of an animal’s ability is the optimal way to detect aging-related changes (Bellantuono et al., 2020). We hypothesized that aging-related changes could also be detected without an overt challenge, by employing sufficiently sensitive and temporally dense measurements during the course of normal daily living, a concept analogous to measuring one’s breathing rate after a brisk walk instead of after sprinting 100 m. To test this idea, we constructed a system of 64 semi-automated phenotyping cages that continuously measured multiple aspects of mouse physiology and behavior (Figure 1A). Cages record oxygen and carbon dioxide concentration, water vapor, weight of food and water hoppers, mouse body mass (when inside the sleeping chamber), running wheel rotations, and mouse position via a three-dimensional array of infrared beams. Gas measurements are acquired every 3 min because the time constant of the cage (rate at which air is replaced) only allows for meaningful changes in gas concentrations every ~3 min. Other measurements are acquired every second. Because food and water are readily available, mice can be monitored continuously, without disruption, for up to weeks at a time. We chose to monitor each mouse for 1 week per month (Figure 1B), as we reasoned that (1) a 7-day run would provide enough data for a robust and stable estimation of physiological/behavioral status and (2) monthly rotations would be sufficiently granular to capture the pace of aging in mice. This schedule allowed us to regularly monitor fourfold as many animals as we had cages, that is, 256. We placed the same animal in different cages for different runs, to mitigate any cage-specific measurement effects.

Figure 1 with 1 supplement see all

Download asset Open asset

We initially enrolled 256 female DO mice (Churchill et al., 2012) split among four age groups (n = 64 at 7, 14, 21, and 25 months), and we enrolled younger (~3 months) animals in semi-regular waves as the older cohorts expired. Enrolling animals at different ages ensured that different ages were similarly represented and mitigated the loss of power that occurs with time due to mortality. Enrolling animals at different calendar dates ensured that age and calendar date were not perfectly correlated, attenuating the effect of any unknown, but potentially confounding, environmental effects such as seasonal changes or temperature drift over time. We utilized DO mice in order to avoid strain-specific results that would be less likely to generalize; however, as a practical necessity, we only examined female mice in this study. In total, we enrolled 415 animals, with an age distribution that provided large numbers of animals at both young and old ages (Figure 1—figure supplement 1A). The study was run for approximately 1.5 years, during which 212 mice died (51%) and 203 mice were censored at the end of the study (49%). This design is similar to human longitudinal studies in which subjects are enrolled at different ages and at different dates, but followed for similar lengths of time, leading to left-truncated and right-censored data. Accounting for truncation and censoring, DO mice had a median lifespan of 25 months (~2 years) and a maximal lifespan of 41 months (~3.4 years) (Figure 1C, Figure 1—figure supplement 1B). In total, after quality control filtering, we obtained data for 3143 7-day runs, that is, ~60 years of dense time series data, from an age range of 6 weeks to 40.5 months (3.4 years), with virtually all ages containing data from multiple cohorts (Figure 1D). For analysis of aging in the following sections, we grouped runs into 3-month age bins, providing a large number of runs (>200) in most bins (Figure 1E). Even with this staggered enrollment design, animals older than 33 months were rare; so for the oldest animals, we combined all runs for animals 33–41 months, for a total of 87 runs in the oldest age bin.

Because of their genetic diversity, DO mice are generally more variable than inbred strains. Indeed, we found a broad body mass distribution in this cohort, ranging from 15 to 60 g (Figure 1—figure supplement 1C). Given this heterogeneity, we sought evidence that the phenotyping cages, which were optimized for less diverse strains, were functioning properly and capturing physiologically meaningful information. In a large, multi-location analysis of indirect calorimetry data, the most consistent relationship was a positive correlation between body mass and energy expenditure (Corrigan et al., 2020). This positive relationship has been proposed as one indication of the technical success of an experiment (Tschöp et al., 2011). We compared body mass to energy expenditure (kcal/hr) in this dataset and saw a highly significant correlation (Figure 1—figure supplement 1D). As we already had body mass information, we were primarily interested in changes to energy expenditure and related parameters that were body mass independent, therefore for all gas-derived measurements except RQ (VO₂, VCO₂, VH₂O, and energy expenditure), because it is a ratio, we normalized for body mass via linear regression. This removed the positive correlation with body mass (Figure 1—figure supplement 1E). This normalization had the potential to induce spurious correlations between gas-derived measurements, so we compared the correlation coefficients between these measurements before and after body weight correction (Figure 1—figure supplement 1F). As expected, there was some correlation between most gas measurements prior to correction. This correlation increased slightly, but not appreciably, upon body mass correction (up to a 0.04 increase in R). We believe that the benefit of focusing on body mass-independent information outweighs the slight increase in correlation between gas measurements that we observed, so we chose to use body mass-corrected gas measurements for all subsequent analysis steps.

Automated phenotyping identifies robust physiological changes at each stage of life

We developed a data processing pipeline that provides aggregated, quality-controlled data at various levels of resolution (Figure 2A). Initially, data are processed into 14 ‘base features’ in 3 min bins: BodyMass (weight), Food (food intake), Water (water intake), WheelMeters (distance traveled on wheel), PedMeters (distance traveled in cage), AllMeters (PedMeters + fine motor movements), XBreak (number of photobeam breaks in the X axis), YBreak (number of photobeam breaks in the Y axis), ZBreak (number of photobeam breaks in the vertical axis), VO₂ (oxygen consumption), VCO₂ (carbon dioxide production), energy expenditure (kcal/hr), RQ (respiratory quotient − VCO₂/VO₂), VH₂O (water vapor). After quality control, each run was split into 24 hr windows and those windows were averaged to generate the average 24 hr trace for each feature for each run. We then grouped runs into the 3-month age bins described above and examined the 14 base features for signs of aging-related change (Figure 2—figure supplement 1A). As expected, animals exhibited circadian patterns of behavior and physiology, with increased activity during the dark phase. Clear aging-related change were observed for most of the 14 base features, starting from the youngest ages we examined. One of the most striking changes we observed was a concave body mass trajectory: animals gained weight until 21–24 months of age, then progressively lost weight (Figure 2B). Other striking changes included a significant decline in wheel running, food intake, water intake, and a decline in mass-adjusted energy expenditure and rate of whole-body water loss (VH₂O) (Figure 2B). Consistent with these observations, it has been previously reported that body weight declines in older mice despite a decrease in energy expenditure (Petr et al., 2021). Not all base features declined with age – ambulatory activity (PedMeters, AllMeters, and YBreaks) tended to increase with age (Figure 2—figure supplement 1A). None of these changes mirrored the concave trajectory of body mass, demonstrating that multiple aging-related changes to physiology and behavior are not driven solely by weight changes.

Figure 2 with 3 supplements see all

Download asset Open asset

Aging-related changes could arise from altered behavior and/or physiology. For example, reduced distance run on the wheel could arise from reduced running speed or reduced time spent running, which might arise from different underlying mechanisms. To quantify behavior, we fit a robust hidden Markov model (HMM) that uses base features of the phenotyping cages, excluding body mass, to classify each 3 min time window into one of six states. A mouse can perform many different behaviors in 3 min, meaning that each ‘state’ is a combination of individual behaviors. Each state was characterized by a distinct combination of base features (Figure 2—figure supplement 1B), and we named the states Sleep, Rest, Eat, Eat&Drink, Run, and Active based on our qualitative interpretation of those patterns. Active, Eat, and Eat&Drink in particular are mixed states, characterized by different combinations of eating, drinking, movement around the cage, and running on the wheel. As measured by energy expenditure, Sleep, Rest, and Eat are lower activity states, whereas Eat&Drink, Run, and Active are higher activity states. As expected, mice organized their behavior around the light cycle, with sharp transitions at lights on/off. They spent the vast majority of the light phase in lower activity states, whereas high activity state occupancy increased in the dark phase (Figure 2—figure supplement 1C). There were differences in overall state occupancy between young (3–6 months) and old (30–33 months) mice; compared to young mice, old mice spent more time in low activity states (Figure 2C). This was true during both light and dark phases, though it was particularly obvious during the dark phase, where old mice spent far less time Running or Active than young mice and more time in the other four, lower activity, states.

Circadian averages are a coarse description of animal behavior, and the temporal resolution of the data allows for more thorough examination. At the hourly (Figure 2D) and 4-hr level (Figure 2—figure supplement 1D), we observed a number of interesting patterns of mouse behavior. Animals exhibited a sharp increase in activity when the lights turned off, reaching their highest activity in the early part of the dark phase, about 2 hr after lights off. During this time, young animals spent nearly 70% of their time Running or Active. Older animals also exhibited their highest activity during this time, but their peak activity was much lower than young mice, and their transition to high activity when the lights initially turned off was slower. Running and Active then generally declined in frequency throughout the rest of the dark phase. Mid-level activity states, that is, Eat and Eat&Drink, exhibited a different pattern: animals spent most of their time in these states near the light/dark transitions. At the end of the dark phase, in particular, this corresponded to a large increase in food and water intake (see Figure 2—figure supplement 1A), likely indicating that the animals were ‘stocking up’ before the low-activity light phase. Prior to this large increase in food and water intake, animals reached a dark phase minimum for food/water intake and a maximum for Sleep and Rest, that is, they took a nap. Napping behavior exhibited different aging dynamics than the decline in high activity noted above; whereas Run and Active declined until ~24 months of age, napping increased in frequency until ~18 months and then stabilized (see Rest during dark phase time windows, Figure 2—figure supplement 1D).

We also saw evidence of circadian phase shifts with age, though not in the manner we expected. We assumed older animals would show less entrainment to the light cycle as they aged, indicative of circadian deregulation. However, we found that old animals were actually more entrained to the light cycle than young animals: young animals reached their peak Eat and Eat&Drink occupancy 1 hr after the lights turned on, during what is generally considered the low activity, ‘sleep’ phase of the circadian cycle (Figure 2D). In contrast, older animals reached their peak Eat and Eat&Drink occupancy in the hour prior to lights on, then quickly transitioned to Sleep and Rest once the lights turned on. The reason for this is unclear – it may be that young mice organize their circadian behavior around wheel running instead of eating, that is, they spend so much of the dark phase running that eating necessarily occurs during the light phase (Edgar and Dement, 1991; Yamanaka et al., 2013; Yasumoto et al., 2015). As wheel running declines with age, there is more time during the dark phase for eating, so peak Eat&Drink occupancy shifts (Valentinuzzi et al., 1997). Whatever the causal chain, this increased entrainment of feeding behavior to the light cycle with age represents an aging-related circadian shift that, to our knowledge, has not been previously described.

Irrespective of age, all mice reached a peak Sleep occupancy of ~50% in the middle of the light phase, and overall sleep duration did not change with age (Supplementary file 1). However, the temporal granularity of the data allowed for an assessment of sleep fragmentation, that is, the duration of continuous periods of sleep, for which we did observed an aging effect. In young mice, periods of continuous sleep averaged 1.1 hr in duration; this declined quickly with age, plateauing at ~0.6 hr by 9–12 months of age, representing a 30 min reduction in the duration of an average sleep period (Figure 2E). At the same time, number of sleeping periods increased with age, representing increased sleep fragmentation. Increased behavioral fragmentation was not observed for all states: in contrast, the average duration of feeding periods increased with age while number of feeding periods decreased (Figure 2F), representing a consolidation of feeding behavior. Increased food spilling has been reported in aging mice (Starr and Saito, 2012); this may result in old mice requiring more time than young mice to ingest an equal quantity of food and may explain the extended duration of feeding periods with age.

The above analyses showed that murine aging is characterized by marked changes in behavioral patterns. The aging-related changes in base features noted above (Figure 2—figure supplement 1A) are thus affected by changes in both behavior and underlying physiology. To decouple these contributions, we computed averages of each of the 14 base features conditioned on the HMM state and compared how these changed with age (Figure 2—figure supplement 2A). We included two additional features – WheelSpeed and PedSpeed – which are the speed of wheel running and pedestrian locomotion, respectively, for a total of 16 base features. Needless to say, these two rate metrics are the most physiologically informative when the animal is actively running or walking, and the plots show that calculated running and walking speed during other states (e.g., Sleep) is under-estimated because the animal spends considerable time not running or walking during these states.

Examining energy expenditure while mice were resting demonstrated that the decline in overall mass-adjusted energy expenditure was, at least partly, a decline in resting metabolic rate – mice exhibited a decline in low-activity-state (Rest and Sleep) energy expenditure with age (Figure 2—figure supplement 2B). Energy expenditure during high activity states (Run and Active) was also decreased, indicating reduced activity intensity. Consistent with this, animals ran more slowly as they aged throughout most of life (Figure 2—figure supplement 2C). Animals also exhibited a aging-related decline in mass-adjusted VH₂O, that is, whole-body water loss rate, during all states, including Rest and Sleep (Figure 2—figure supplement 2D), suggesting that this is a consequence of altered internal physiology rather than altered behavior. Interestingly, the decline in water loss rate precedes the decline in overall water intake (see Figure 2A), so reduced water intake does not explain this effect; it may instead represent aging-related changes in body composition, as lean and fat mass retain different amounts of water (Vu et al., 2017).

Many other features changed with age as well, so in order to obtain a holistic view of aging patterns, we generated 309 features from the raw data using a number of different approaches, including those mentioned above: (1) the per-run average of each of the 16 base features noted above + average occupancy of the 6 HMM states ($16+6=22$ features), (2) the 16 base features conditioned on each state ($16\times 6=96$ features), (3) average values of the 16 base features and 6 HMM states during each of six 4-hr time bins (${\displaystyle [(16+6)]\times 6=132}$ features), (4) the ratio of values from before and after a light transition for each of the 16 base features and 6 HMM states ($(16+6)\times 2=44$ features), and (5) period number, average and maximum period length, average and maximum time interval between periods for eating, sleeping, and high-activity ($5\times 3=15$ features) (see Figure 2A, green box). We chose these different aggregation approaches because we reasoned they had the highest chance of capturing unique aspects of physiology, though additional features may be informative. For example, the length of sleeping periods might show a different aging pattern than overall sleep amount, and the ratio of pedestrian locomotion before and after a light transition might provide different physiological information than pedestrian locomotion at a particular time of day.

After multiple hypothesis correction, 244 features (79%) significantly changed with age (Supplementary file 1). Features exhibited a variety of trajectories and ages of onset (Figure 2G), many of which were distinct from the base features, demonstrating that feature engineering uncovered additional physiological information. We chose a selection of three such features to highlight the different ages of onset we observe (Figure 2H). This diversity of aging patterns shows that (1) passive phenotyping has the sensitivity to detect aging across the entire lifespan and (2) physiological aging is more complex than simple linear trajectories and thus both dense time sampling and multi-dimensional phenotyping are needed to accurately capture the many different domains of physiological aging.

Quantifying relationships between physiological and behavioral features demonstrates a decline in resilience with age

Thus far we had examined how individual features change with age, but because the interplay between physiological domains is hypothesized to be a primary determinant of system resilience (Cohen, 2016), we were also interested in the connectivity between these features and how that connectivity changes with age. The strong genetic and phenotypic correlations noted above motivated us to explore this in more detail. A unique feature of this dataset is that all features were measured simultaneously in all runs, which allowed us to quantify the covariance between features, that is, how each feature is correlated with every other feature across runs. More precisely, we fit a graphical model (a semi-parametric Gaussian copula model) using sparse precision matrix estimation methods and visualized the model as a force-directed network (Figure 3A). We determined the number of clusters and clustering of features via consensus clustering, which resulted in 22 clusters. Some of these clusters were defined by the physical sensor used to generate the data, for example, VH₂O and CO₂ production, but most were comprised of features from multiple sensors. For example, beam break features did not form a single cluster, nor did all food, water, or gas features. Instead, clusters were defined by combinations of behaviors and aggregation methods: cluster 16 represents overall metabolic rate; cluster 12 represents food and water intake in a specific time window (the light phase); and cluster 13 represents the ratio of activity before/after a light transition, etc. (Supplementary file 2). Ratio-based features, state-conditioned features, and time-of-day conditioned features often clustered separately from one another; this demonstrates that different data aggregation approaches convey different information. Furthermore, this suggests that nodes within a cluster share similar physiological underpinnings rather than simply similar sensor hardware.

Figure 3 with 1 supplement see all

Download asset Open asset

Some clusters were tightly integrated with other clusters, whereas others exhibited sparse connectivity, demonstrating their independence from other physiological domains. To better visualize this, we generated a chord diagram in which each cluster is placed around a circle, with chords demonstrating connectivity between nodes from different clusters (Figure 3B). The radius of each cluster arc is proportional to its overall connectivity to other clusters, thus larger clusters indicate more connections (not the number of nodes in the cluster). One of the most well-defined and least-connected clusters was body mass, as demonstrated by its small arc radius: all body mass metrics formed a tight cluster that was largely independent of other clusters. The tight clustering between mass metrics arises because body mass is mostly invariant with respect to behavior and time of day, so different body mass features from the same run are highly correlated to one another. The independence from other clusters arises partly because, as noted above, gas-derived measurements other than RQ (VO₂, VCO₂, kcal/hr, VH₂O) were adjusted for body mass. However, body mass was also largely independent of features that were not mass-adjusted, such as wheel, photobeam, food, and water features. The low connectivity of body mass to other phenotype clusters demonstrates that variance in body mass does not drive all, or even most, of the variance in other features in this dataset. This conclusion is further supported by the observation that body mass exhibits a clear concave trajectory with age, whereas most other features do not. In contrast, metabolic rate and high activity pattern were the most densely connected clusters, indicating that these physiological characteristics impact (or are impacted by) many other physiological domains.

We next sought to understand whether network connectivity changed with age, as a readout of system resilience. Resilience refers to the ability of a system to maintain function in the face of change. This is a broad concept that is unlikely to be fully captured by a single number, and multiple approaches to measuring organism-level resilience have been proposed (Huffman et al., 2016). Here, we propose a metric that is based on the relationship between physiological features. In a resilient system, a small disruption in one physiological domain will have limited system-wide effect, whereas in a less resilient system, the same perturbation will affect more of the network. Thus, in less resilient systems, individual subsystems (nodes) are sensitive to the fluctuations of other subsystems, leading to a rising correlation between subsystems (Scheffer et al., 2018). For a graphical network in which correlation between nodes is visualized by edge weight, this manifests as increased network connectivity. Conversely, in a resilient system, because a small disruption in one node only affects a small number of other nodes, the average correlation between nodes is low, leading to lower network connectivity. Keeping node positions constant, we calculated network edges for each age bin (Figure 3C, Figure 3—figure supplement 1A). Visual inspection of these per-age networks from 3 to 6 months to 33+ months suggested increasing connectivity with time, suggesting decreasing system resilience.

To quantify network connectivity directly, we calculated the copula entropy/negative multivariate mutual information (MI) (Singh and Póczos, 2017) of the graphical model for each age bin. Multivariate MI is a measurement of how much information can be determined about the network as a whole from a subset of nodes. In one extreme case where all nodes are perfectly correlated (low resilience), knowing the value of a single node allows one to perfectly predict the value of all other nodes. Such a system has zero degrees of freedom and thus high MI. In the other extreme, nodes are completely independent (high resilience), provide no information about one another, and all must be measured to achieve a full description of the system – low MI. Thus, the negative of MI directly corresponds to a more rigid network and is a quantitative measurement of system resilience. Resilience (negative MI) exhibited a striking decrease with age, with an accelerating rate of decline; this pattern qualitatively reflects the accelerating failure rate associated with aging better than the semi-linear or early-life dropoff we observed with many individual parameters in this study (Figure 3D). Thus, resilience may be more accurate representations of physiological aging than are individual features.

To better understand what network-level changes were driving this decline in resilience with age, we examined both inter-cluster connectivity as well as intra-cluster connectivity, that is, we examined the relationship between physiological domains as well as the relationship between related features within physiological domains. To quantify inter-cluster connectivity, we represented each cluster by the exemplar feature that most reflected the behavior of that cluster (as measured by median consensus with other features in the same cluster). We then re-built networks using only those exemplars and calculated resilience. Like overall resilience, inter-cluster resilience exhibited a striking decrease with age (Figure 3E), suggesting the connectivity between subsystems increases with age. This effect was not merely driven by a few outlying clusters; increased connectivity with age was observed for the majority of clusters (Figure 3F). We also examined intra-cluster resilience, by measuring the negative MI between nodes within each cluster (Figure 3—figure supplement 1B). Most clusters showed a decline in resilience with age (increasing connectivity), but trajectories varied: some clusters exhibited a nearly linear decline in resilience, whereas others exhibited a late-life acceleration, perhaps indicating that certain domains of physiology are more prone to a loss of resilience with age than others.

Age and time to death can be simultaneously predicted from physiological data

An increasingly common approach when dealing with multi-dimensional aging data is to train models that predict relevant endpoints (so-called aging ‘clocks’). This approach was popularized using DNA methylation data (Horvath, 2013), and more recent reports use other data modalities including physiological measurements (Schultz et al., 2020). These studies have developed models that predict either chronological age or mortality risk separately, but these separate approaches are limited because neither chronological age nor mortality risk fully encapsulates an individual’s multi-system, aging-related health status, that is, biological age.

To address this, we developed an aging rate regression model in which biological age is determined from a combination of chronological age and health status (in this case, time to death, though other health proxies such as a frailty score could be used). The model includes a hyperparameter that allows for tuning of the relative weighting of chronological age and time to death, allowing us to generate models with different behaviors. More specifically, this hyperparameter (denoted ${\sigma }_{\beta }$) quantifies our belief that different individuals age at different rates. If a ground truth measurement of individual aging rates existed, this hyperparameter could be measured empirically. Unfortunately, there remains no agreed-upon definition of biological age and no such ground truth is available. Therefore, here we explore model behavior under several different values of ${\sigma }_{\beta }$. A low value of ${\sigma }_{\beta }$ causes the model to assume that all individuals age at similar rates, meaning that the biological age of individuals of the same chronological age should be similar. In this case, model training heavily weighs chronological age, and the resulting model approximates a standard age clock model. Conversely, a high value of ${\sigma }_{\beta }$ causes the model to assume that individuals can age at different rates, and thus model training disregards chronological age, instead emphasizing health status (time to death), and the resultant model approximates a standard accelerated failure time model. Neither chronological age nor time to death are perfect representations of aging rate, and they are not particularly well correlated with one another (Figure 4—figure supplement 1A), thus optimizing the prediction of one necessarily reduces performance for the other, resulting in a tunable tension in model behavior and the ability to explore intermediate states that may avoid overfitting to either of these imperfect biological age surrogates. Because this framework utilizes both chronological age and survival time as outcome variables, we name this approach the “Combined Age and Survival Prediction of Aging Rate”, or CASPAR.

To understand how this trade-off behavior between chronological age and time to death manifests, and to determine if there is an acceptable middle ground, we trained models using different values of ${\sigma }_{\beta }$. For each model, we trained on the same 15 different random splits of the data and report the average performance of the model on the held-out mice. In each split, we allocated 90% of mice to training and 10% of mice to testing, with full animal holdouts (i.e., all runs of the same animal were either entirely in the training set or entirely in the test set). Note that because the model utilizes time to death, we only included the mice that died during the course of the study, approximately half the mice in the overall dataset. Leveraging the longitudinal nature of our dataset to improve the robustness of the model to run-specific effects, during model training we assume that the aging rate of an animal is constant over time, that is, that the observed aging rate of the same animal at different ages is modeled by a latent common aging rate plus a standard normal noise term.

At the lowest values of ${\sigma }_{\beta }$ we tested (e.g., ${\sigma }_{\beta }=0.1$, heavy weighting of chronological age), when assessed against the test sets, the model’s scores were highly correlated to chronological age, with a mean correlation across the 15 data splits of 0.77 (Figure 4A). However, as expected, this model performed poorly on time to death prediction, with a mean correlation between model scores and time to death of 0.33 (11% variance explained). Not surprisingly, in this case the model barely performs better than true chronological age, which explains about 8% of the variance in time to death (${\displaystyle \mathrm{p}<1\mathrm{e}-4}$) (Figure 4—figure supplement 1A). This behavior represented one extreme of model performance, as the correlations for age and time to death were stable for all values of ${\sigma }_{\beta }\le 1$ (Figure 4B).

Figure 4 with 1 supplement see all

Download asset Open asset

In contrast, at the highest value of ${\sigma }_{\beta }$ (e.g., ${\sigma }_{\beta }=100$, heavy weighting of time to death), the model’s scores were better correlated to actual time to death with a mean correlation of 0.47 (Figure 4A). As expected, correlation with the chronological age was lower in this model, dropping to a correlation of 0.58. This behavior represented the other extreme of model performance, as the correlations for age and time to death were stable for all values of ${\sigma }_{\beta }\ge 30$ (Figure 4B). All 15 training/test splits of the data led to similar model performance, demonstrating that model performance is robust to the specific animals allocated to the training and test sets (Figure 4B).

Interestingly, the trade-off between age and time to death prediction was not linear – intermediate values of ${\sigma }_{\beta }$ improved prediction in one dimension without an equivalent loss of performance in the other dimension (Figure 4A). To better visualize this, we chose three values of ${\sigma }_{\beta }$ to represent three different regimes of the CASPAR model: the chronological age model (${\sigma }_{\beta }=0.1$), the time to death model (${\sigma }_{\beta }=100$), and the hybrid model, which represents a roughly equivalent trade-off between age and time to death (${\sigma }_{\beta }=3$). For each model class, we plotted actual age and actual time to death versus predicted age for each run in the test set, for each of the 15 models trained on different data splits (Figure 4C). As noted earlier, the chronological age model predicted age quite well ($R^2=0.60$, ${\displaystyle \mathrm{p}<1\mathrm{e}-4}$) but did not predict time to death much better than chronological age itself ($R^2=0.11$, $p=0.001$), whereas the time to death model predicted time to death relatively well ($R^2=0.22$, ${\displaystyle p<1e-4}$) but chronological age relatively poorly ($R^2=0.34$, ${\displaystyle p<1e-4}$). Interestingly, the hybrid model performed relatively well on both tasks, predicting chronological age with an $R^2=0.50$ (${\displaystyle p<1e-4}$) and predicting time to death with an $R^2=0.20$ (${\displaystyle p<1e-4}$), which was substantially better than the true chronological age of the animal and nearly as good as the time to death model itself. This balanced prediction of both dimensions suggests that the hybrid model may represent a qualitatively different approach to quantifying health-relevant biological age than either age or survival prediction alone.

We examined the models to understand which features were most informative. For all three model classes (age, time to death, and hybrid), the top 20 most informative features came from multiple clusters and sensors (gas, photobeams, wheel, body mass) and aggregation methods (state-conditioned, time of day, lights on/off ratio) (Figure 4—figure supplement 1B, Figure 4—figure supplement 1C). No single cluster or sensor dominated, demonstrating the importance of a high-dimensional feature set. The features that informed each of the model classes were partially overlapping (Figure 4D). Of the top 20 features from each model, 7 were shared across all three (Figure 4—figure supplement 1C). Perhaps unsurprisingly, these include major aspects of physiology – overall size (BodyMass), activity (wheel running and pedestrian locomotion), and a surrogate of basal metabolic rate (VCO₂ while Sleeping). Outside of this overlap, quite a few features were distinct to either the age model or time to death model. Features more associated with time to death than with age (‘TTD and Hybrid only’ and ‘TTD only’) include measures of activity (wheel, food, and water intake) during the light phase or low activity states, perhaps indicating that a breakdown in circadian stability is an indicator of impending demise. Conversely, features more associated with age than time to death (‘Age only’) include the extent of activity during the dark phase, perhaps indicating that the best predictor of chronological age is a decline in high intensity activity. Limited overlap between age and mortality prediction models has been reported elsewhere using different measures of physiological age (Schultz et al., 2020), thus this may represent a more general truth about organism-level aging: aging-related changes in physiological function, though of clear health relevance to the individual, are only partially related to mortality risk, demonstrating the need for a thoughtful balance between these conceptualizations of biological age.

Assessments of behavior and physiology are essential aspects of many preclinical studies. However, while technological advances such as sequencing have allowed researchers to explore molecular and cellular phenotypes in high-dimensional space via systems-level analyses, organism-level phenotyping has not benefited from a similar advance. To remedy this, we combined automated phenotyping cages with a sophisticated analysis pipeline to create a a platform for high-dimensional assessment of physiology and behavior in mice. This platform could be utilized to study multiple organism-level processes and diseases, for example, cognitive and mood disorders, neuromuscular deficits, or metabolic disease. We chose to focus on aging because (1) there is increasing interest in therapeutically modulating aging, (2) aging affects multiple physiological domains, so we expected broad and clear effects, and (3) current approaches to quantifying organism-level aging are extremely labor- and time-intensive (Bellantuono et al., 2020; Sukoff Rizzo et al., 2018). We were able to measure over 200 aging phenotypes across 22 physiological clusters, as well as network effects and related emergent phenomena (resilience) that are difficult to quantify via other methods, in >400 outbred mice. Data generation required ~20 hr of hands-on operator time per week for 1.5 years, that is, 0.75 FTE-years total. This study advances the state of the art for healthspan assessment in terms of throughput, resolution, and physiological scope.

At the outset of this study, it was not clear that monitoring of mice in a normal living environment would provide sufficient sensitivity to detect age-related changes. The rationale behind utilizing challenge-based assays is that animals must be pushed to the limit of their abilities in order to quantify functional decline (Sukoff Rizzo et al., 2018; Bellantuono et al., 2020). Although some aspects of aging undoubtedly require such assays, we found that our data show a plethora of aging-related changes in all age groups, including prior to 6 months of age, demonstrating that automated phenotyping of voluntary activity can detect even early aging-related changes. Mice reach sexual maturity at ~6 weeks of age, but are arguably still undergoing developmental changes for weeks or months after that. This shifts the question from when we can detect aging-related changes to when we should. Deciding whether aging-related changes in animals younger than 6 months of age represents aging, development, or a combination of the two is a complex issue that cannot be definitively resolved here. However, we did notice that phenotypes exhibited a diversity of trajectories across life – unchanged, parabolic, linear, logarithmic, and nearly exponential – and we propose that phenotypes which change monotonically, particularly when the trajectory is linear, can be considered part of aging even when that change begins early. For example, wheel running declines near-monotonically with age starting at 3 months; it seems reasonable to propose that this is an aspect of aging that begins quite early. Changes that are non-monotonic, for example, body weight, are more difficult to interpret through this lens.

One difficulty that arises when interpreting automated phenotyping data is distinguishing changes in physiology from changes in behavior. Although this is a somewhat arbitrary distinction, it is meaningful – a decline in overall energy expenditure because an animal runs less is different from a decline in energy expenditure due to reduced basal metabolic rate. To address this, we developed a robust HMM to assign a behavioral state to each 3 min time bin. We were then able to examine features conditioned on the state of the animal, for example, VO₂ while running. This turned out to be an informative approach, as state-conditioned features arising from the same sensor often clustered separately from one another, indicating that they contained complementary physiological information (e.g., pedestrian locomotion while eating versus while resting). More qualitatively, it allowed us to assess specific aspects of physiology that have long been considered the domain of specialized procedures. For example, energy expenditure while sleeping provides a reasonable surrogate of resting metabolic rate whereas VO₂ while running provides a reasonable surrogate of maximal oxygen consumption, both of which decline with age. We propose that behavioral inference based on automated phenotyping data is a useful technique that can be applied in a number of preclinical contexts.

A powerful use of multi-dimensional data is the application of network analysis to better understand the wiring of the system. We built a network of physiological and behavioral phenotypes using sparse precision matrix estimation methods; in other words, we determined every pairwise correlation between features after accounting for all other features. In this context, two features connected by a strong edge (high covariance) are more likely to be mechanistically connected than two features with a weaker edge, and thus features that form a cluster are likely driven by the same mechanism(s), whereas features from a different cluster are likely driven by distinct mechanism(s). Aging involves a multitude of detrimental changes to health and well-being, but it is unknown how many distinct causal mechanisms drive these changes (Freund, 2019). Analyzing aging as a phenotypic network informs how many clusters exist, therefore how many independent mechanisms likely exist. The concept is analogous to identifying clusters of coordinately expressed genes; with sufficient data, one can conclude that coordinately expressed genes are regulated by a similar transcriptional program, and the number of gene clusters provides a reasonable estimate of the number of transcriptional programs. In this study, we identified 22 distinct organism-level clusters, though this number will undoubtedly be refined as additional studies are incorporated into the network analysis framework.

Network modeling also allows for quantification of overall network connectivity, which serves as a useful proxy for resilience. We uncovered striking changes with age in this dimension: more than virtually any individual feature, resilience smoothly and monotonically declined with age. The decline in resilience accelerated with age, decreasing more rapidly in old animals than in young animals, a pattern that is qualitatively consistent with the accelerating decline in health and corresponding exponential increase in mortality with age. This emergent property of the system was only detectable by analyzing the relationship between phenotypes, rather than the individual values of the phenotypes themselves – an analytical approach that is infeasible for data from challenge-based assays. Increasing network connectivity indicates that, with age, runs increasingly resemble one another, that is, old animals occupy a lower diversity of phenotypic states than younger animals. A similar phenomenon has been reported for human frailty – inter-individual variation in frailty scores decreases with age (Rockwood et al., 2004). It may be that individual animals become more physiologically inflexible with age; alternatively, this observation may be a consequence of survival bias, as DO mice reach 50% mortality by ~2 years of age, and only a small percentage of animals survive to >30 months. This latter explanation would imply that there are a limited number of viable aging trajectories, and animals that do not follow those trajectories die early, removing their contribution to phenotypic diversity. Longitudinal assessment of resilience in individuals would help address this question, and future work could develop measurements of individual animal resilience, rather than the population-level resilience we have calculated here. In particular, the second-by-second and minute-by-minute data streams from phenotyping cages could be used to examine how similar an animal is to itself at a later time. This ‘temporal autocorrelation’ analysis is only possible with time series data and has been proposed as a useful measure of individual resilience (Scheffer et al., 2018).

High-dimensional data is information-rich but also difficult to interpret. As such, there is value in developing summary statistics that capture meaningful aspects of the data. In the case of aging, this leads to the concept of biological age – a single number that is meant to reflect the aging-related health status of an animal better than does its chronological age. An increasingly popular method to quantify biological age from multi-dimensional data is to train a model to predict chronological age, then treat the error in that model (i.e., the difference between actual age and predicted age) as biologically meaningful. In many cases, this so-called ‘delta age’ correlates to meaningful outcomes, such as survival. Despite its popularity, the conceptual foundation of this approach is questionable. Consider that the computational goal of an age prediction model is to predict chronological age perfectly, but perfect performance would destroy the utility of the model. A model that predicts all 70-year-old people, healthy or unhealthy, to be exactly 70 provides no information other than chronological age itself; $R^2$ = 1 is not useful. However, although it is clear that perfect performance is not desirable, it is unclear how much error is optimal; that is, it is unclear whether the error should be 1%, 5%, 20%, or more. Further, the desire for error is not explicitly stated during model training – algorithms push toward perfect performance whenever possible and do not aim to capture any other type of information. At some point, improved age prediction begins to ignore the biological information we care about, but we do not know when that point is reached because we have not defined a ‘ground truth’ measurement of the variable we actually care about, that is, biological age. In short, if you don’t know what you want, you’re unlikely to get it. The non-circular approach is to define a measurement of biological age that is independent of the training data and use that as the outcome variable when training models. However, because we don’t fully understand aging, and because health is a difficult concept to define, let alone measure, it is difficult to confidently create such a ‘ground truth’ biological age metric.

A common alternative is to use lifespan, or survival time, as the outcome variable. Survival time is objectively measurable and fulfills the criteria of being independent of the training data, but it is still not exactly what we care about from an aging perspective: not all individuals with equal remaining lifespan are equally healthy or equally ‘aged’. Lifespan is influenced by a large number of external factors, and mortality risk is often dominated by particular pathologies (e.g., cancer, in the case of many mouse strains). An anti-cancer drug or any other single-cause-of-death modulator does not meet most people’s criteria of an aging intervention, but would be viewed as such if survival time were the single outcome variable assessed.

As neither chronological age prediction nor survival time prediction is an optimal method to evaluate biological age, we developed a compromise: CASPAR, which incorporates both chronological age and survival time prediction. CASPAR explicitly captures two key notions: (1) biological age should be correlated to chronological age and (2) deviations of biological age from chronological age should be informative of some proxy of health status (e.g., time to death). These two assumptions are often used as secondary validation metrics for mortality-only or age-only models, respectively: biological age predictions from mortality-trained models are tested for their correlation with chronological age, and biological age predictions from chronological-age models are tested for their ability to predict mortality. However, our approach explicitly quantifies these assumptions, allowing for smooth interpolation between chronological age regression and survival regression. We trained multiple models, varying the relative importance of age versus survival time prediction, and evaluated their performance on held-out animals. Age prediction was always more accurate than survival time prediction, but this is not surprising – the former is an exercise in estimating a concurrent variable (chronological age), whereas the latter predicts a future outcome that is dependent on a number of stochastic external influences that have not yet occurred and may have no relationship to an animal’s current state. One extreme, the ‘chronological age model’, predicted chronological age well but had poor survival time prediction. Conversely, the ‘time to death’ model predicted survival time well, but showed substantially reduced age prediction performance. Therefore, both age and mortality information are present in the data, but the feature weighting that predicts them is largely distinct. A similar dichotomy between age and survival time prediction has been reported in other physiological datasets (Schultz et al., 2020; Fischer et al., 2016; Levine et al., 2018). Thus, considerable thought must go into model design, as different outcome variables are likely to emphasize different aspects of biology. We also demonstrate the existence of a hybrid model that gives some weight to both age and survival time prediction. This hybrid model predicted survival time nearly as well as a standard survival model, but gained substantial improvement in chronological age prediction. Notably, predicted biological age from the hybrid model was more strongly correlated to time to death than chronological age itself. It is beyond the scope of this manuscript to dictate exactly where the age versus survival prediction trade-off should lie, but we hypothesize that such hybrid models are likely to be more valuable and aging-relevant than either extreme.

The CASPAR approach also allowed us to take advantage of longitudinal measurements of the same mouse. Versions of the model in which we assumed that aging rates are relatively constant throughout life, thus allowing us to use multiple runs of the same animal to estimate a single age rate for that animal, performed significantly better than models without this additional constraint on held-out validation splits.

The model framework we have developed, with the ability to tune the relative weighting of age and survival time, can be applied more broadly. First, the underlying data does not need to be physiological. A similar approach can be used with any other high-dimensional, per-individual dataset, including methylation data, transcriptomic data, and blood biomarkers. For the purposes of preclinical intervention testing, we favor physiological data because it is, almost by definition, health-relevant. Our models relied on features like body mass, running speed, energy expenditure, and sleeping behavior – this is intuitively sensible and provides confidence that the resulting model truly reflects animal function. In contrast, although molecular data is often easier to acquire and contains more features, there is limited prior information on each individual feature, meaning that molecular data-based models are more difficult to interpret and sanity check. A second way to broaden the application of this modeling framework is by using an outcome variable other than death. This may be particularly useful in human datasets, where follow-up mortality data is often limited. Age could be combined with a moribidity endpoint, or a pre-determined ‘health score’, and these outcome variables could be used to train the model. This may allow biological age assessment in cohorts with limited follow-up time.

Our study has several important limitations. First, it is not a fully longitudinal design. Mice were enrolled starting at different ages and followed for a fixed length of time (~1.5 years), thus our data are both left-truncated and right-censored. Having full lifespan curves for all animals would provide additional insight. Human clinical trial data has a similar left-censored, right-truncated data structure, so our ability to quantify health-relevant aging changes despite incomplete life histories bodes well for translationally relevant intervention testing using this platform. Second, the cohort was entirely female. This was a practical necessity; due to their aggressive tendencies, DO males are generally singly housed at all times, leading to unsustainable financial and vivarium space requirements. Nevertheless, a mixed dataset would be preferable, not only in terms of sex, but also in terms of strain and environmental variables, because mixed datasets increase the generalizability of conclusions (Voelkl et al., 2020, Webster and Rutz, 2020). It is quite likely that the specific physiological changes that develop with age will differ in male mice, in other strains, and in animals housed under different conditions; however, we suspect that broad patterns such as reduced wheel running, metabolic activity, and resilience, will generalize. In this manuscript, we have primarily focused on the analytical methodology and tools we developed to understand and summarize the physiological changes we observed because we believe those will be more generalizable, and thus useful, to the field than the particular set of physiological changes that we identify.

Additionally, although the automated monitoring cages we used provided a tremendous amount of data, they are blind to many health-relevant phenotypes such as body posture and coat condition. Additional high content, automated phenotyping modalities, such as video monitoring paired with machine vision feature extraction, would increase the value of this platform. Fourth, our most granular analyses used 3 min time windows. However, for many of the cage sensors, the data are acquired every second. This provides the opportunity for more sophisticated time series analyses that track individual behaviors and rapid fluctuations in physiology.

Changes with age occur at all levels of biological organization (molecules, cells, tissues, etc.). We chose to focus on physiology and behavior because changes to this level of biological organization are most proximal to changes in health and quality of life. We created a platform that can measure physiological and behavioral aging at any stage of life in outbred mice using automated phenotyping. This system provides a number of advances that allow organism-level aging to be studied with improved throughput, resolution, and physiological scope while reducing the activation energy that comes with highly specialized assay procedures. We encourage more widespread adoption of automated phenotyping and high-dimensional analysis in order to study aging and putative aging interventions.

All data and associated code are available on Github (https://github.com/calico/catnap, (copy archived at swh:1:rev:2f18eea0d02c23501bdd36558822a6974f99f640)).

1. 1. Chen Z
   2. Raj A
   3. Prateek GV
   4. Francesco AD
   5. Liu J
   6. Keyes BE
   7. Kolumam G
   8. Jojic V
   9. Freund A

   (2022) Github

   ID Github. Automated, high-dimensional evaluation of physiological aging and resilience in outbred mice.

   https://github.com/calico/catnap

1. 1. Ackert-Bicknell CL
   2. Anderson LC
   3. Sheehan S
   4. Hill WG
   5. Chang B
   6. Churchill GA
   7. Chesler EJ
   8. Korstanje R
   9. Peters LL

   (2015) Aging Research Using Mouse Models

   Current Protocols in Mouse Biology 5:95–133.

   https://doi.org/10.1002/9780470942390.mo140195

   • PubMed
   • Google Scholar
2. 1. Bellantuono I
   2. de Cabo R
   3. Ehninger D
   4. Di Germanio C
   5. Lawrie A
   6. Miller J
   7. Mitchell SJ
   8. Navas-Enamorado I
   9. Potter PK
   10. Tchkonia T
   11. Trejo JL
   12. Lamming DW

   (2020) A toolbox for the longitudinal assessment of healthspan in aging mice

   Nature Protocols 15:540–574.

   https://doi.org/10.1038/s41596-019-0256-1

   • PubMed
   • Google Scholar
3. 1. Broman KW
   2. Gatti DM
   3. Simecek P
   4. Furlotte NA
   5. Prins P
   6. Sen Ś
   7. Yandell BS
   8. Churchill GA

   (2019) R/qtl2: Software for Mapping Quantitative Trait Loci with High-Dimensional Data and Multiparent Populations

   Genetics 211:495–502.

   https://doi.org/10.1534/genetics.118.301595

   • PubMed
   • Google Scholar
4. 1. Brown SDM
   2. Holmes CC
   3. Mallon AM
   4. Meehan TF
   5. Smedley D
   6. Wells S

   (2018) High-throughput mouse phenomics for characterizing mammalian gene function

   Nature Reviews. Genetics 19:357–370.

   https://doi.org/10.1038/s41576-018-0005-2

   • PubMed
   • Google Scholar
5. Software

   1. Chen Z

   (2002) calico/catnap, version swh:1:rev:2f18eea0d02c23501bdd36558822a6974f99f640

   GitHub.

   https://archive.softwareheritage.org/swh:1:dir:21cb0b936f54e859757b50758d337a9a2164c993;origin=https://github.com/calico/catnap;visit=swh:1:snp:0e7844388acd7f01f38f968b6f84048165ddc4a7;anchor=swh:1:rev:2f18eea0d02c23501bdd36558822a6974f99f640
6. Conference

   1. Chen T
   2. Guestrin C

   (2016) Xgboost: A scalable tree boosting system

   Proceedings of the 22nd acm sigkdd international conference on knowledge discovery and data mining. pp. 785–794.

   https://doi.org/10.1145/2939672.2939785

   • Google Scholar
7. Preprint

   1. Chen Z
   2. Raj A
   3. Prateek GV
   4. Di Francesco A
   5. Liu J
   6. Keyes BE
   7. Kolumam G
   8. Jojic V
   9. Freund A

   (2021) Automated, High-Dimensional Evaluation of Physiological Aging and Resilience in Outbred Mice

   bioRxiv.

   https://doi.org/10.1101/2021.08.02.454830

   • Google Scholar
8. 1. Chesler EJ
   2. Wilson SG
   3. Lariviere WR
   4. Rodriguez-Zas SL
   5. Mogil JS

   (2002a) Identification and ranking of genetic and laboratory environment factors influencing a behavioral trait, thermal nociception, via computational analysis of a large data archive

   Neuroscience and Biobehavioral Reviews 26:907–923.

   https://doi.org/10.1016/s0149-7634(02)00103-3

   • PubMed
   • Google Scholar
9. 1. Chesler EJ
   2. Wilson SG
   3. Lariviere WR
   4. Rodriguez-Zas SL
   5. Mogil JS

   (2002b) Influences of laboratory environment on behavior

   Nature Neuroscience 5:1101–1102.

   https://doi.org/10.1038/nn1102-1101

   • PubMed
   • Google Scholar
10. 1. Churchill GA
    2. Gatti DM
    3. Munger SC
    4. Svenson KL

    (2012) The Diversity Outbred mouse population

    Mammalian Genome 23:713–718.

    https://doi.org/10.1007/s00335-012-9414-2

    • PubMed
    • Google Scholar
11. 1. Cohen AA

    (2016) Complex systems dynamics in aging: new evidence, continuing questions

    Biogerontology 17:205–220.

    https://doi.org/10.1007/s10522-015-9584-x

    • PubMed
    • Google Scholar
12. 1. Corrigan JK
    2. Ramachandran D
    3. He Y
    4. Palmer CJ
    5. Jurczak MJ
    6. Chen R
    7. Li B
    8. Friedline RH
    9. Kim JK
    10. Ramsey JJ
    11. Lantier L
    12. McGuinness OP
    13. Banks AS
    14. Mouse Metabolic Phenotyping Center Energy Balance Working Group

    (2020) A big-data approach to understanding metabolic rate and response to obesity in laboratory mice

    eLife 9:e53560.

    https://doi.org/10.7554/eLife.53560

    • PubMed
    • Google Scholar
13. 1. Crabbe JC
    2. Wahlsten D
    3. Dudek BC

    (1999) Genetics of mouse behavior: interactions with laboratory environment

    Science (New York, N.Y.) 284:1670–1672.

    https://doi.org/10.1126/science.284.5420.1670

    • PubMed
    • Google Scholar
14. 1. Edgar DM
    2. Dement WC

    (1991) Regularly scheduled voluntary exercise synchronizes the mouse circadian clock

    The American Journal of Physiology 261:R928–R933.

    https://doi.org/10.1152/ajpregu.1991.261.4.R928

    • PubMed
    • Google Scholar
15. 1. Fischer KE
    2. Hoffman JM
    3. Sloane LB
    4. Gelfond JAL
    5. Soto VY
    6. Richardson AG
    7. Austad SN

    (2016) A cross-sectional study of male and female C57BL/6Nia mice suggests lifespan and healthspan are not necessarily correlated

    Aging 8:2370–2391.

    https://doi.org/10.18632/aging.101059

    • PubMed
    • Google Scholar
16. 1. Freund A

    (2019) Untangling Aging Using Dynamic, Organism-Level Phenotypic Networks

    Cell Systems 8:172–181.

    https://doi.org/10.1016/j.cels.2019.02.005

    • PubMed
    • Google Scholar
17. 1. Furlotte NA
    2. Eskin E

    (2015) Efficient multiple-trait association and estimation of genetic correlation using the matrix-variate linear mixed model

    Genetics 200:59–68.

    https://doi.org/10.1534/genetics.114.171447

    • PubMed
    • Google Scholar
18. 1. Gatti DM
    2. Svenson KL
    3. Shabalin A
    4. Wu LY
    5. Valdar W
    6. Simecek P
    7. Goodwin N
    8. Cheng R
    9. Pomp D
    10. Palmer A
    11. Chesler EJ
    12. Broman KW
    13. Churchill GA

    (2014) Quantitative trait locus mapping methods for diversity outbred mice

    G3: Genes, Genomes, Genetics 4:1623–1633.

    https://doi.org/10.1534/g3.114.013748

    • PubMed
    • Google Scholar
19. Conference

    1. Hallac D
    2. Park Y
    3. Boyd S
    4. Leskovec J

    (2017) Network Inference via the Time-Varying Graphical Lasso

    KDD: Proceedings. International Conference on Knowledge Discovery & Data Mining. pp. 205–213.

    https://doi.org/10.1145/3097983.3098037

    • Google Scholar
20. 1. Hansen M
    2. Kennedy BK

    (2016) Does Longer Lifespan Mean Longer Healthspan?

    Trends in Cell Biology 26:565–568.

    https://doi.org/10.1016/j.tcb.2016.05.002

    • PubMed
    • Google Scholar
21. 1. Horvath S

    (2013) DNA methylation age of human tissues and cell types

    Genome Biology 14:R115.

    https://doi.org/10.1186/gb-2013-14-10-r115

    • PubMed
    • Google Scholar
22. 1. Huffman DM
    2. Justice JN
    3. Stout MB
    4. Kirkland JL
    5. Barzilai N
    6. Austad SN

    (2016) Evaluating Health Span in Preclinical Models of Aging and Disease: Guidelines, Challenges, and Opportunities for Geroscience

    The Journals of Gerontology. Series A, Biological Sciences and Medical Sciences 71:1395–1406.

    https://doi.org/10.1093/gerona/glw106

    • PubMed
    • Google Scholar
23. 1. Kafkafi N
    2. Agassi J
    3. Chesler EJ
    4. Crabbe JC
    5. Crusio WE
    6. Eilam D
    7. Gerlai R
    8. Golani I
    9. Gomez-Marin A
    10. Heller R
    11. Iraqi F
    12. Jaljuli I
    13. Karp NA
    14. Morgan H
    15. Nicholson G
    16. Pfaff DW
    17. Richter SH
    18. Stark PB
    19. Stiedl O
    20. Stodden V
    21. Tarantino LM
    22. Tucci V
    23. Valdar W
    24. Williams RW
    25. Würbel H
    26. Benjamini Y

    (2018) Reproducibility and replicability of rodent phenotyping in preclinical studies

    Neuroscience and Biobehavioral Reviews 87:218–232.

    https://doi.org/10.1016/j.neubiorev.2018.01.003

    • PubMed
    • Google Scholar
24. 1. Kang HM
    2. Zaitlen NA
    3. Wade CM
    4. Kirby A
    5. Heckerman D
    6. Daly MJ
    7. Eskin E

    (2008) Efficient control of population structure in model organism association mapping

    Genetics 178:1709–1723.

    https://doi.org/10.1534/genetics.107.080101

    • PubMed
    • Google Scholar
25. 1. Kim SJ
    2. Koh K
    3. Boyd S
    4. Gorinevsky D

    (2009) $\ell_1$ Trend Filtering

    SIAM Review 51:339–360.

    https://doi.org/10.1137/070690274

    • Google Scholar
26. 1. Ladiges W
    2. Snyder JM
    3. Wilkinson E
    4. Imai DM
    5. Snider T
    6. Ge X
    7. Ciol M
    8. Pettan-Brewer C
    9. Pillai SPS
    10. Morton J
    11. Quarles E
    12. Rabinovitch P
    13. Niedernhofer L
    14. Liggitt D

    (2017) A New Preclinical Paradigm for Testing Anti-Aging Therapeutics

    The Journals of Gerontology. Series A, Biological Sciences and Medical Sciences 72:760–762.

    https://doi.org/10.1093/gerona/glx019

    • PubMed
    • Google Scholar
27. 1. Le KN
    2. Zhan M
    3. Cho Y
    4. Wan J
    5. Patel DS
    6. Lu H

    (2020) An automated platform to monitor long-term behavior and healthspan in Caenorhabditis elegans under precise environmental control

    Communications Biology 3:297.

    https://doi.org/10.1038/s42003-020-1013-2

    • PubMed
    • Google Scholar
28. 1. Levine ME
    2. Lu AT
    3. Quach A
    4. Chen BH
    5. Assimes TL
    6. Bandinelli S
    7. Hou L
    8. Baccarelli AA
    9. Stewart JD
    10. Li Y
    11. Whitsel EA
    12. Wilson JG
    13. Reiner AP
    14. Aviv A
    15. Lohman K
    16. Liu Y
    17. Ferrucci L
    18. Horvath S

    (2018) An epigenetic biomarker of aging for lifespan and healthspan

    Aging 10:573–591.

    https://doi.org/10.18632/aging.101414

    • PubMed
    • Google Scholar
29. 1. Liao CY
    2. Rikke BA
    3. Johnson TE
    4. Diaz V
    5. Nelson JF

    (2010) Genetic variation in the murine lifespan response to dietary restriction: from life extension to life shortening

    Aging Cell 9:92–95.

    https://doi.org/10.1111/j.1474-9726.2009.00533.x

    • PubMed
    • Google Scholar
30. 1. Liu H
    2. Lafferty J
    3. Wasserman L

    (2009)

    The nonparanormal: Semiparametric estimation of high dimensional undirected graphs

    Journal of Machine Learning Research 10:2295–2328.

    • Google Scholar
31. 1. Mandillo S
    2. Tucci V
    3. Hölter SM
    4. Meziane H
    5. Banchaabouchi MA
    6. Kallnik M
    7. Lad HV
    8. Nolan PM
    9. Ouagazzal AM
    10. Coghill EL
    11. Gale K
    12. Golini E
    13. Jacquot S
    14. Krezel W
    15. Parker A
    16. Riet F
    17. Schneider I
    18. Marazziti D
    19. Auwerx J
    20. Brown SDM
    21. Chambon P
    22. Rosenthal N
    23. Tocchini-Valentini G
    24. Wurst W

    (2008) Reliability, robustness, and reproducibility in mouse behavioral phenotyping: a cross-laboratory study

    Physiological Genomics 34:243–255.

    https://doi.org/10.1152/physiolgenomics.90207.2008

    • PubMed
    • Google Scholar
32. 1. Martineau CN
    2. Brown AEX
    3. Laurent P

    (2020) Multidimensional phenotyping predicts lifespan and quantifies health in Caenorhabditis elegans

    PLOS Computational Biology 16:e1008002.

    https://doi.org/10.1371/journal.pcbi.1008002

    • PubMed
    • Google Scholar
33. 1. McIlwain KL
    2. Merriweather MY
    3. Yuva-Paylor LA
    4. Paylor R

    (2001) The use of behavioral test batteries: effects of training history

    Physiology & Behavior 73:705–717.

    https://doi.org/10.1016/s0031-9384(01)00528-5

    • PubMed
    • Google Scholar
34. Conference

    1. Mohan K
    2. Chung MJY
    3. Han S
    4. Witten D
    5. Lee SI
    6. Fazel M

    (2012)

    Structured Learning of Gaussian Graphical Models

    Advances in Neural Information Processing Systems. pp. 629–637.

    • Google Scholar
35. 1. Monti S
    2. Tamayo P
    3. Mesirov J
    4. Golub T

    (2003) Consensus clustering: a resampling-based method for class discovery and visualization of gene expression microarray data

    Machine Learning 52:91–118.

    https://doi.org/10.1023/A:1023949509487

    • Google Scholar
36. 1. Morgan AP
    2. Fu CP
    3. Kao CY
    4. Welsh CE
    5. Didion JP
    6. Yadgary L
    7. Hyacinth L
    8. Ferris MT
    9. Bell TA
    10. Miller DR
    11. Giusti-Rodriguez P
    12. Nonneman RJ
    13. Cook KD
    14. Whitmire JK
    15. Gralinski LE
    16. Keller M
    17. Attie AD
    18. Churchill GA
    19. Petkov P
    20. Sullivan PF
    21. Brennan JR
    22. McMillan L
    23. Pardo-Manuel de Villena F

    (2015) The Mouse Universal Genotyping Array: From Substrains to Subspecies

    G3: Genes, Genomes, Genetics 6:263–279.

    https://doi.org/10.1534/g3.115.022087

    • PubMed
    • Google Scholar
37. 1. Petr MA
    2. Alfaras I
    3. Krawcyzk M
    4. Bair W-N
    5. Mitchell SJ
    6. Morrell CH
    7. Studenski SA
    8. Price NL
    9. Fishbein KW
    10. Spencer RG
    11. Scheibye-Knudsen M
    12. Lakatta EG
    13. Ferrucci L
    14. Aon MA
    15. Bernier M
    16. de Cabo R

    (2021) A cross-sectional study of functional and metabolic changes during aging through the lifespan in male mice

    eLife 10:e62952.

    https://doi.org/10.7554/eLife.62952

    • PubMed
    • Google Scholar
38. 1. Richardson A
    2. Fischer KE
    3. Speakman JR
    4. de Cabo R
    5. Mitchell SJ
    6. Peterson CA
    7. Rabinovitch P
    8. Chiao YA
    9. Taffet G
    10. Miller RA
    11. Rentería RC
    12. Bower J
    13. Ingram DK
    14. Ladiges WC
    15. Ikeno Y
    16. Sierra F
    17. Austad SN

    (2016) Measures of Healthspan as Indices of Aging in Mice-A Recommendation

    The Journals of Gerontology. Series A, Biological Sciences and Medical Sciences 71:427–430.

    https://doi.org/10.1093/gerona/glv080

    • PubMed
    • Google Scholar
39. 1. Rockwood K
    2. Mogilner A
    3. Mitnitski A

    (2004) Changes with age in the distribution of a frailty index

    Mechanisms of Ageing and Development 125:517–519.

    https://doi.org/10.1016/j.mad.2004.05.003

    • PubMed
    • Google Scholar
40. Book

    1. Romano JP
    2. DiCiccio C

    (2019)

    Multiple Data Splitting for Testing

    Stanford, United States: Stanford University.

    • Google Scholar
41. 1. Scheffer M
    2. Bolhuis JE
    3. Borsboom D
    4. Buchman TG
    5. Gijzel SMW
    6. Goulson D
    7. Kammenga JE
    8. Kemp B
    9. van de Leemput IA
    10. Levin S
    11. Martin CM
    12. Melis RJF
    13. van Nes EH
    14. Romero LM
    15. Olde Rikkert MGM

    (2018) Quantifying resilience of humans and other animals

    PNAS 115:11883–11890.

    https://doi.org/10.1073/pnas.1810630115

    • PubMed
    • Google Scholar
42. 1. Schultz MB
    2. Kane AE
    3. Mitchell SJ
    4. MacArthur MR
    5. Warner E
    6. Vogel DS
    7. Mitchell JR
    8. Howlett SE
    9. Bonkowski MS
    10. Sinclair DA

    (2020) Publisher Correction: Age and life expectancy clocks based on machine learning analysis of mouse frailty

    Nature Communications 11:5143.

    https://doi.org/10.1038/s41467-020-19046-8

    • PubMed
    • Google Scholar
43. Conference

    1. Singh S
    2. Póczos B

    (2017)

    Nonparanormal Information Estimation

    Proceedings of the 34th International Conference on Machine Learning. pp. 3210–3219.

    • Google Scholar
44. 1. Sittig LJ
    2. Carbonetto P
    3. Engel KA
    4. Krauss KS
    5. Barrios-Camacho CM
    6. Palmer AA

    (2016) Genetic Background Limits Generalizability of Genotype-Phenotype Relationships

    Neuron 91:1253–1259.

    https://doi.org/10.1016/j.neuron.2016.08.013

    • PubMed
    • Google Scholar
45. 1. Starr ME
    2. Saito H

    (2012) Age-related increase in food spilling by laboratory mice may lead to significant overestimation of actual food consumption: implications for studies on dietary restriction, metabolism, and dose calculations

    The Journals of Gerontology. Series A, Biological Sciences and Medical Sciences 67:1043–1048.

    https://doi.org/10.1093/gerona/gls009

    • PubMed
    • Google Scholar
46. 1. Sukoff Rizzo SJ
    2. Anderson LC
    3. Green TL
    4. McGarr T
    5. Wells G
    6. Winter SS

    (2018) Assessing Healthspan and Lifespan Measures in Aging Mice: Optimization of Testing Protocols, Replicability, and Rater Reliability

    Current Protocols in Mouse Biology 8:e45.

    https://doi.org/10.1002/cpmo.45

    • PubMed
    • Google Scholar
47. 1. Svenson KL
    2. Gatti DM
    3. Valdar W
    4. Welsh CE
    5. Cheng R
    6. Chesler EJ
    7. Palmer AA
    8. McMillan L
    9. Churchill GA

    (2012) High-resolution genetic mapping using the Mouse Diversity outbred population

    Genetics 190:437–447.

    https://doi.org/10.1534/genetics.111.132597

    • PubMed
    • Google Scholar
48. 1. Tschöp MH
    2. Speakman JR
    3. Arch JRS
    4. Auwerx J
    5. Brüning JC
    6. Chan L
    7. Eckel RH
    8. Farese RV
    9. Galgani JE
    10. Hambly C
    11. Herman MA
    12. Horvath TL
    13. Kahn BB
    14. Kozma SC
    15. Maratos-Flier E
    16. Müller TD
    17. Münzberg H
    18. Pfluger PT
    19. Plum L
    20. Reitman ML
    21. Rahmouni K
    22. Shulman GI
    23. Thomas G
    24. Kahn CR
    25. Ravussin E

    (2011) A guide to analysis of mouse energy metabolism

    Nature Methods 9:57–63.

    https://doi.org/10.1038/nmeth.1806

    • PubMed
    • Google Scholar
49. 1. Tuttle AH
    2. Philip VM
    3. Chesler EJ
    4. Mogil JS

    (2018) Comparing phenotypic variation between inbred and outbred mice

    Nature Methods 15:994–996.

    https://doi.org/10.1038/s41592-018-0224-7

    • PubMed
    • Google Scholar
50. 1. Valentinuzzi VS
    2. Scarbrough K
    3. Takahashi JS
    4. Turek FW

    (1997) Effects of aging on the circadian rhythm of wheel-running activity in C57BL/6 mice

    The American Journal of Physiology 273:R1957–R1964.

    https://doi.org/10.1152/ajpregu.1997.273.6.R1957

    • PubMed
    • Google Scholar
51. 1. Voelkl B
    2. Altman NS
    3. Forsman A
    4. Forstmeier W
    5. Gurevitch J
    6. Jaric I
    7. Karp NA
    8. Kas MJ
    9. Schielzeth H
    10. Van de Casteele T
    11. Würbel H

    (2020) Reproducibility of animal research in light of biological variation

    Nature Reviews. Neuroscience 21:384–393.

    https://doi.org/10.1038/s41583-020-0313-3

    • PubMed
    • Google Scholar
52. 1. Vu JP
    2. Luong L
    3. Parsons WF
    4. Oh S
    5. Sanford D
    6. Gabalski A
    7. Lighton JR
    8. Pisegna JR
    9. Germano PM

    (2017) Long-Term Intake of a High-Protein Diet Affects Body Phenotype, Metabolism, and Plasma Hormones in Mice

    The Journal of Nutrition 147:2243–2251.

    https://doi.org/10.3945/jn.117.257873

    • PubMed
    • Google Scholar
53. 1. Webster MM
    2. Rutz C

    (2020) How STRANGE are your study animals?

    Nature 582:337–340.

    https://doi.org/10.1038/d41586-020-01751-5

    • Google Scholar
54. 1. Wei LJ

    (1992) The accelerated failure time model: a useful alternative to the Cox regression model in survival analysis

    Statistics in Medicine 11:1871–1879.

    https://doi.org/10.1002/sim.4780111409

    • PubMed
    • Google Scholar
55. Conference

    1. Wen Q
    2. Gao J
    3. Song X
    4. Sun L
    5. Xu H
    6. Zhu S

    (2019) RobustSTL: A Robust Seasonal-Trend Decomposition Algorithm for Long Time Series

    Proceedings of the AAAI Conference on Artificial Intelligence. pp. 5409–5416.

    https://doi.org/10.1609/aaai.v33i01.33015409

    • Google Scholar
56. Preprint

    1. Wright KM
    2. Deighan A
    3. Di Francesco A
    4. Freund A
    5. Jojic V
    6. Churchill G
    7. Raj A

    (2020) Age and Diet Shape the Genetic Architecture of Body Weight in Diversity Outbred Mice

    bioRxiv.

    https://doi.org/10.1101/2020.11.04.364398

    • Google Scholar
57. 1. Yamanaka Y
    2. Honma S
    3. Honma K -i

    (2013) Daily exposure to a running wheel entrains circadian rhythms in mice in parallel with development of an increase in spontaneous movement prior to running-wheel access

    American Journal of Physiology. Regulatory, Integrative and Comparative Physiology 305:R1367–R1375.

    https://doi.org/10.1152/ajpregu.00389.2013

    • PubMed
    • Google Scholar
58. 1. Yasumoto Y
    2. Nakao R
    3. Oishi K

    (2015) Free access to a running-wheel advances the phase of behavioral and physiological circadian rhythms and peripheral molecular clocks in mice

    PLOS ONE 10:e0116476.

    https://doi.org/10.1371/journal.pone.0116476

    • PubMed
    • Google Scholar
59. 1. Yuan R
    2. Peters LL
    3. Paigen B

    (2011) Mice as a mammalian model for research on the genetics of aging

    ILAR Journal 52:4–15.

    https://doi.org/10.1093/ilar.52.1.4

    • PubMed
    • Google Scholar
60. 1. Zhang WB
    2. Pincus Z

    (2016) Predicting all-cause mortality from basic physiology in the Framingham Heart Study

    Aging Cell 15:39–48.

    https://doi.org/10.1111/acel.12408

    • PubMed
    • Google Scholar
61. 1. Zhang WB
    2. Sinha DB
    3. Pittman WE
    4. Hvatum E
    5. Stroustrup N
    6. Pincus Z

    (2016) Extended Twilight among Isogenic C. elegans Causes a Disproportionate Scaling between Lifespan and Health

    Cell Systems 3:333–345.

    https://doi.org/10.1016/j.cels.2016.09.003

    • PubMed
    • Google Scholar

Author details

1. Zhenghao Chen

   Calico Life Sciences LLC, South San Francisco, South San Francisco, United States

   Contribution

   Data curation, Formal analysis, Methodology, Software, Visualization, Writing – original draft, Writing – review and editing

   Competing interests

   is affiliated with Calico Life Sciences LLC. The author has no financial interests to declare

2. Anil Raj

   Calico Life Sciences LLC, South San Francisco, South San Francisco, United States

   Contribution

   Data curation, Formal analysis, Methodology, Software, Visualization, Writing – original draft, Writing – review and editing

   Competing interests

   is affiliated with Calico Life Sciences LLC. The author has no financial interests to declare

   "This ORCID iD identifies the author of this article:" 0000-0003-4412-0883

3. GV Prateek

   Calico Life Sciences LLC, South San Francisco, South San Francisco, United States

   Contribution

   Data curation, Formal analysis, Methodology, Software, Writing – original draft, Writing – review and editing

   Competing interests

   is affiliated with Calico Life Sciences LLC. The author has no financial interests to declare

4. Andrea Di Francesco

   Calico Life Sciences LLC, South San Francisco, South San Francisco, United States

   Contribution

   Investigation, Methodology, Resources, Writing – original draft, Writing – review and editing

   Competing interests

   is affiliated with Calico Life Sciences LLC. The author has no financial interests to declare

   "This ORCID iD identifies the author of this article:" 0000-0001-6867-8203

5. Justin Liu

   Calico Life Sciences LLC, South San Francisco, South San Francisco, United States

   Contribution

   Investigation, Methodology, Resources

   Competing interests

   is affiliated with Calico Life Sciences LLC. The author has no financial interests to declare

   "This ORCID iD identifies the author of this article:" 0000-0002-5338-6491

6. Brice E Keyes

   Calico Life Sciences LLC, South San Francisco, South San Francisco, United States

   Contribution

   Investigation, Methodology, Resources

   Competing interests

   was affiliated with Calico Life Sciences LLC during the time of his contribution to the study. The author has no financial interests to declare

7. Ganesh Kolumam

   Calico Life Sciences LLC, South San Francisco, South San Francisco, United States

   Contribution

   Supervision

   Competing interests

   is affiliated with Calico Life Sciences LLC. The author has no financial interests to declare

8. Vladimir Jojic

   Calico Life Sciences LLC, South San Francisco, South San Francisco, United States

   Contribution

   Data curation, Formal analysis, Methodology, Software, Supervision, Writing – original draft, Writing – review and editing

   Competing interests

   was affiliated with Calico Life Sciences LLC during the time of his contribution to the study. The author has no financial interests to declare

9. Adam Freund

   Calico Life Sciences LLC, South San Francisco, South San Francisco, United States

   Contribution

   Conceptualization, Formal analysis, Funding acquisition, Project administration, Resources, Supervision, Writing – original draft, Writing – review and editing

   For correspondence

   adam@afbio.co

   Competing interests

   was affiliated with Calico Life Sciences LLC during the time of his contribution to the study. The author has no financial interests to declare

   "This ORCID iD identifies the author of this article:" 0000-0002-7956-5332

• 2,832

  views

• 274

  downloads

• 13

  citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.