1. Introduction {#sec1-animals-10-01443}
===============

Excessive body mass, i.e., being overweight or obese, is a health concern. Such conditions are associated with fat-levels that may impair health \[[@B1-animals-10-01443]\], and problems such as diabetes, heart disease, impaired reproduction, orthopaedic disorders, and cancers \[[@B2-animals-10-01443],[@B3-animals-10-01443],[@B4-animals-10-01443],[@B5-animals-10-01443],[@B6-animals-10-01443],[@B7-animals-10-01443]\]. Methods of determining body condition include weight-for-height (kg/m^2^) indices in humans ("body mass index": 1) and other primates (e.g., \[[@B8-animals-10-01443]\]), skinfold thickness (e.g., \[[@B9-animals-10-01443]\]), and visual body condition scoring systems (e.g., \[[@B10-animals-10-01443],[@B11-animals-10-01443],[@B12-animals-10-01443]\]). Body mass is also used as a proxy for body fat: individuals ≥20--25% over the ideal are considered obese (e.g., in humans, dogs, *Canis lupus familaris*, and cats, *Felis catus*: \[[@B7-animals-10-01443],[@B13-animals-10-01443]\]). Positive energy balance, i.e., calorific intake greater than expenditure, is central to weight gain \[[@B13-animals-10-01443],[@B14-animals-10-01443]\] but is not the sole factor. Energy balance is affected by environmental, e.g., diet and physical activity levels, and genetic factors \[[@B13-animals-10-01443]\], e.g., obesity-risk varies between human \[[@B15-animals-10-01443]\] and non-human animal populations \[[@B6-animals-10-01443],[@B14-animals-10-01443],[@B16-animals-10-01443]\], and is heritable (e.g., vervet monkeys, *Chlorocebus pygerythrus*: \[[@B17-animals-10-01443]\]).

Lemuriformes ("lemurs") are a primate infraorder native to Madagascar, and these species vary in susceptibility to weight gain in captivity. For instance, a healthy body condition is typical of greater bamboo, *Prolemur simus*, and red-bellied lemurs, *Eulemur rubriventer* \[[@B18-animals-10-01443],[@B19-animals-10-01443]\]. Conversely, ring-tailed, *Lemur catta*, and blue-eyed black lemurs, *E. flavifrons*, are prone to elevated body masses \[[@B18-animals-10-01443],[@B19-animals-10-01443]\], and for some species problems such as obesity are prevalent. For example, using a weight-based method of inferring body condition, 46.5% black-and-white ruffed and red ruffed lemurs, *Varecia variegata* and *V. rubra*, housed in European zoos were deemed obese (defined by the authors as ≥two standard deviations from the wild mean weight: \[[@B20-animals-10-01443]\]). Similarly, of 902 zoo lemurs from 14 species 54% were considered overweight or obese (likewise defined as ≥ two and ≥ four standard deviations from species-typical wild mean weights, respectively: \[[@B18-animals-10-01443]\]). Furthermore, obesity can hinder breeding programmes, and obese animals are unsuitable for wild reintroduction \[[@B21-animals-10-01443]\]---concerning for a taxonomic group as threatened as lemurs (cf. \[[@B22-animals-10-01443]\]). Therefore, to improve health and promote conservation, understanding the basis for some species' susceptibility to weight gain is fundamental.

Three aspects of species-typical ecology might explain species differences in susceptibility to weight gain. The first relates to lemurs' physiological and behavioural adaptations regarding fat-storage \[[@B23-animals-10-01443]\], to buffer against future resource-restriction (i.e., "thrifty genotypes" sensu: \[[@B24-animals-10-01443]\]). Madagascar has relatively poor soil quality and plant productivity, with a harsh and unpredictable climate, the severity of which varies geographically \[[@B25-animals-10-01443]\]. Rainfall affects plant growth \[[@B26-animals-10-01443],[@B27-animals-10-01443]\] and annual rainfall is commonly used to infer primary productivity and thus food availability \[[@B28-animals-10-01443],[@B29-animals-10-01443],[@B30-animals-10-01443]\]. Annual rainfall varies across Madagascar, e.g., northern areas have a very long dry season and the south receives very little rainfall, unlike eastern areas \[[@B31-animals-10-01443],[@B32-animals-10-01443],[@B33-animals-10-01443]\]. Additionally, there are regional differences in between-year annual rainfall, e.g., the north and south are especially variable; hypothesised to explain unpredictable fruit availability \[[@B34-animals-10-01443]\]. The frequency and severity of droughts, cyclones and El Niño events also vary across the island and between years---events which also restrict food resources \[[@B25-animals-10-01443],[@B35-animals-10-01443],[@B36-animals-10-01443],[@B37-animals-10-01443]\]. Therefore, while all lemur species naturally experience relatively poor and unpredictable food resources, this will be more severe for some than others, according to where they reside on the island. Additionally, captive and wild food quality likely differ. Commercial fruits especially are more energy- and sugar-rich than wild equivalents---a mismatch cited in the excessive body masses observed in some captive lemurs (e.g., \[[@B38-animals-10-01443],[@B39-animals-10-01443]\]). In captivity, this mismatch combined with thrifty adaptations to naturally poor and/or unpredictable food resources might leave certain species susceptible to weight gain.

The second aspect of species-typical ecology that might explain species variation in captive weight gain, is arboreality. Ring-tailed lemurs, for instance, are deemed "semi-terrestrial" \[[@B40-animals-10-01443]\], aye ayes, *Daubentonia madagascariensis*, also spend time at ground-level \[[@B41-animals-10-01443]\], and grey mouse lemurs, *Microcebus murinus*, hunt and occasionally nest on the ground \[[@B42-animals-10-01443]\]. However, other species are arboreal specialists, such as ruffed lemurs \[[@B43-animals-10-01443]\], only coming to ground if absolutely required, e.g., to drink \[[@B44-animals-10-01443]\] or to retrieve fallen young \[[@B45-animals-10-01443]\]. Across mammals, including lemurs, terrestrial species carry higher fat stores than arboreal ones \[[@B46-animals-10-01443]\]. Arboreality imposes upper limits on how much body fat species develop, as carrying excess weight imposes locomotive and agility costs \[[@B46-animals-10-01443],[@B47-animals-10-01443],[@B48-animals-10-01443],[@B49-animals-10-01443]\]. Therefore, being naturally arboreal could be *protective* against captive weight gain, meaning such species should be less prone to excessive body masses than those more terrestrial.

Finally, predation risk could drive species differences in susceptibility to weight gain. High predation pressure is associated with relatively smaller fat stores in birds, reptiles and mammals, because carrying excess weight impedes escape \[[@B48-animals-10-01443],[@B49-animals-10-01443],[@B50-animals-10-01443],[@B51-animals-10-01443],[@B52-animals-10-01443],[@B53-animals-10-01443],[@B54-animals-10-01443]\]. To the best of our knowledge, there is no evidence for such patterns in lemurs. However, lifespan in wild grey mouse lemurs was negatively correlated with body mass, suggested to result from selective predation on larger individuals \[[@B55-animals-10-01443]\]. High wild predation risk, therefore, could explain why some species are less prone to weight gain even when in captivity.

In addition to species-typical ecology, an individual's demographic profile and environment, e.g., diet and physical activity, also affect body condition (e.g., \[[@B13-animals-10-01443]\]). Epidemiological approaches \[[@B56-animals-10-01443]\] are useful in exploring such individual-level effects, having been successfully used to examine risk factors for health and welfare problems across zoo \[[@B57-animals-10-01443]\], farm \[[@B58-animals-10-01443]\], and companion animals \[[@B59-animals-10-01443],[@B60-animals-10-01443]\]. For captive lemurs, usable space and exercise opportunities (sensu \[[@B13-animals-10-01443]\]), the need to maintain body temperature if housed outdoors (sensu \[[@B61-animals-10-01443]\]), use of contraceptives/neutering \[[@B14-animals-10-01443],[@B19-animals-10-01443],[@B62-animals-10-01443]\], amount of fruit in the diet (cf. \[[@B38-animals-10-01443]\]), and activity levels (sensu \[[@B61-animals-10-01443]\]) are biologically relevant aspects that could contribute to susceptibility to weight gain.

Here, we use two complementary methods (cf. \[[@B63-animals-10-01443]\]) to better understand susceptibility to weight gain in captive lemurs. We use phylogenetic comparative methods to test the hypotheses relating species-typical ecology to captive weight gain described above, aiming to identify *species*-level effects, i.e., ecological risk factors. We also use an epidemiological approach, and by comparing across individuals we aim to identify *individual*-level effects, e.g., demographic and/or environmental risk factors for captive weight gain.

2. Materials and Methods {#sec2-animals-10-01443}
========================

Ethical approval for outcome data collection was granted by the Faculty of Health Sciences Research Ethics Committee of the University of Bristol (University Investigation Number: 37201).

2.1. Outcome Variable {#sec2dot1-animals-10-01443}
---------------------

As reviewed above, there are several ways to determine body condition in primates, but most require direct access and/or handling of animals (e.g., waist circumference measurements: \[[@B64-animals-10-01443]\]), and none to-date have been tested and validated for use in lemurs. A practical alternative is to use deviations from a species' normal mass to infer body condition. Many zoos routinely weigh their animals \[[@B65-animals-10-01443]\], and published values of lemur species-typical mean wild adult weights are available \[[@B66-animals-10-01443]\]. Comparisons between species-typical wild mean and captive representatives' body masses have already been used to infer body condition across several lemur species \[[@B18-animals-10-01443]\], and similar weight-based methods of inferring body condition are used in other species (e.g., in humans, dogs, and cats \[[@B7-animals-10-01443],[@B13-animals-10-01443],[@B67-animals-10-01443]\]. We acknowledge, however, that using mass alone to infer body condition is limited because frame is not accounted for (unlike, for example, body mass index: \[[@B1-animals-10-01443]\]); however, frame measurements are rarely readily available for captive wild animals.

Here, we use "relative body mass" as a proxy for body condition, i.e., the ratio between a captive animal's body mass and its species-typical wild mean (method after: \[[@B18-animals-10-01443]\], wild means from: \[[@B66-animals-10-01443]\]). Relative body mass values correspond with the following body conditions (after \[[@B18-animals-10-01443]\]): \<0.75 = *underweight*, i.e., the captive animal is \<75% of its wild mean; 0.75--1.25 = *healthy,* the captive animal is 75--125% of its wild mean; 1.25--1.5 = *overweight,* the captive animal is 125--150% of its wild mean; 1.5--2 = *obese*, the captive animal is 150--200% of its wild mean; and ≥2 = *morbidly obese*, the animal is 200% (or more) of its wild mean.

Our outcome for comparative analyses was species-typical median relative body mass, calculated from values from individuals across 13 species (medians being the preferred species summary statistic to reduce the effects of outliers and skew in the raw data: \[[@B68-animals-10-01443]\]). For epidemiological analyses of the four best-sampled species, our outcome was individuals' relative body masses.

2.2. Outcome Data Collection and Processing {#sec2dot2-animals-10-01443}
-------------------------------------------

Data on body masses and corresponding living conditions of captive lemurs were collected between August 2016--January 2018 mainly using an online survey made in Google Forms (see [Supplementary Materials File S1](#app1-animals-10-01443){ref-type="app"} for the questions). Participants could either provide body mass values within the survey or provide Zoological Information Systems (ZIMS) Specimen Reports for their lemurs. As the online survey was suitable for collections with ≤10 enclosures, larger collections completed a tailored Excel spreadsheet containing the same questions. On request, the Duke Lemur Center ([www.lemur.duke.edu](www.lemur.duke.edu)) kindly provided detailed copies of their most recent records of body masses (an update of the dataset of \[[@B69-animals-10-01443]\]), and housing, feeding and enrichment routines of their lemurs.

Using 'species holding' information from Species360 (then, the International Species Information System \[[@B70-animals-10-01443]\]) we emailed or used zoo website contact forms to recruit zoos known to hold lemurs. After initial contact, up to two further reminder requests were made. Zoos were classed as non-participatory if no response was received after the third email, or if the zoo declined to participate. The British and Irish Association of Zoos and Aquariums (BIAZA) Research Committee, the Association of Zoos and Aquariums (AZA, North America) Prosimian Taxon Advisory Group (TAG), and the Zoo and Aquarium Association (ZAA, Australasia) Primate TAG all supported our project.

Raw online survey responses were downloaded for processing, and responses contained in tailored Excel spreadsheets and the Duke Lemur Center's data were added to this, resulting in records from 1386 animals representing 22 species from 133 collections worldwide. Initially, we processed these raw data into individual-level responses, combining responses from specific questions where necessary to calculate our variables of interest, ready for calculation of species summary statistics for the comparative analyses (see [Table 1](#animals-10-01443-t001){ref-type="table"}). A subset of these individual-level responses from the four best-sampled species ([Section 2.2.2](#sec2dot2dot2-animals-10-01443){ref-type="sec"}) were used for epidemiological analyses. Note that the resulting datasets contain missing data for some variables resulting in sample size differences.

We excluded entries without body mass values, those from species-hybrids (because we were interested in species effects), and females known to be pregnant (as their weights might over-estimate their body condition). We calculated age at weighing from ZIMS Specimen reports or reported dates of birth to determine which lemurs were adult; for the minority of animals lacking this information, we accepted the respondents' judgement of "adult". Lemurs have a late "near-adult" growth period, in which animals are sexually mature but not yet fully-grown \[[@B69-animals-10-01443]\]. After Zehr et al., \[[@B69-animals-10-01443]\] we therefore excluded records from all but adult animals (defined by \[[@B69-animals-10-01443]\] as ≥twice the minimum dam age of reproduction for lemurs housed at the Duke Lemur Center; see [Appendix A](#app2-animals-10-01443){ref-type="app"} [Table A1](#animals-10-01443-t0A1){ref-type="table"} for these). Wild fat-tailed dwarf lemurs, *Cheirogaleus medius*, and grey mouse lemurs undergo seasonal hibernation or torpor, and naturally undergo programmed fattening (sensu \[[@B71-animals-10-01443]\]) prior to this \[[@B72-animals-10-01443],[@B73-animals-10-01443]\], and captive animals likewise demonstrate a similar pattern of annual weight fluctuations \[[@B74-animals-10-01443]\]. For these two species, most of their representatives' weights were taken immediately before their inactive periods (i.e., when they were at their heaviest). Because their relative body masses would be biased upwards, we excluded values from these two species from analyses. After these exclusions we were left with records, including some on living conditions, for 691 animals from 20 species. Our last exclusion was to remove species represented by \<5 animals (after \[[@B75-animals-10-01443],[@B76-animals-10-01443]\]), resulting in a final dataset of 675 animals from 13 species. For each animal in this dataset we calculated the ratio of its most recently taken body mass to its species-typical wild mean, yielding its relative body mass value.

### 2.2.1. Further Survey Processing for Comparative Analyses {#sec2dot2dot1-animals-10-01443}

For comparative analyses we calculated species summary statistics from the final individual-level dataset (described in [Table 1](#animals-10-01443-t001){ref-type="table"}). As some species did not have ≥5 animals per sex (see [Section 2.5.1](#sec2dot5dot1-animals-10-01443){ref-type="sec"}), we were unable to examine sex differences in comparative analyses (obesity and its related effects do differ between the sexes across primates: \[[@B79-animals-10-01443],[@B80-animals-10-01443],[@B81-animals-10-01443]\]). For our outcome variable, we calculated the median relative body mass of species' representatives, this yielding species-typical median relative body mass.

Environment and living conditions also contribute to weight gain \[[@B13-animals-10-01443]\] and, should these covary with our wild ecology predictors, could be confounds \[[@B63-animals-10-01443]\]. Therefore, using the information on living conditions we collected in our survey, we calculated six husbandry variables (see [Table 1](#animals-10-01443-t001){ref-type="table"}) believed *a priori* to likely influence body weight so we could later statistically control for their effects if necessary. For these we calculated an appropriate summary statistic (a median or proportion; see [Table 1](#animals-10-01443-t001){ref-type="table"}).

### 2.2.2. Further Survey Processing for Epidemiological Analyses {#sec2dot2dot2-animals-10-01443}

Early exploratory analyses \[[@B74-animals-10-01443]\] revealed that four of our final 13 species had complete cases across all variables enabling epidemiological analyses (n): ring-tailed lemurs (351), black-and-white ruffed lemurs (89), red ruffed lemurs (75), and mongoose lemurs, *Eulemur mongoz* (29). We excluded records from all other species. Because we were interested in possible sex differences, we also excluded any animal of unknown sex (one animal). Maximum recorded captive lifespans of these four species are similar and thus ages are directly comparable (i.e., mongoose lemur = 36.2 years; ring-tailed lemur = 37.3 years; red ruffed lemur = 37.6 years; black-and-white ruffed lemur = 39.4 years \[[@B82-animals-10-01443]\]). To boost sample sizes, we estimated some dates and corresponding seasons of weighing recorded as "unknown" for the comparative analyses, by assigning the former as the 1st of the month that we received weight records and using this to calculate the corresponding hemisphere-specific season of weighing. For two ring-tailed lemurs with dates of birth but without dates of weighing, we also used these estimated dates to calculate estimated age at weighing. Doing so resulted in them being deemed sub-adult rather than adult as per the comparative analyses, thus excluding them from the epidemiological analyses. Ring-tailed lemurs are adult from 978 days of age (cf. \[[@B69-animals-10-01443]\]). Being very close to adulthood (929 and 936 days), including these individuals in the comparative analyses where they just contributed to the species median is therefore likely to have minimal impact.

2.3. Comparative Analyses Predictor Variable Data Collation {#sec2dot3-animals-10-01443}
-----------------------------------------------------------

For the 13 species featuring in the comparative analyses, we collated data on six wild ecological predictor variables from published literature to test our hypotheses. See [Table 2](#animals-10-01443-t002){ref-type="table"} for details and rationale of wild ecological predictors and their calculations, and [Table 3](#animals-10-01443-t003){ref-type="table"} for species-typical values.

2.4. Predictors for Epidemiological Analyses {#sec2dot4-animals-10-01443}
--------------------------------------------

Details of husbandry and demographic predictor variables calculated from our survey and their rationale are shown in [Table 1](#animals-10-01443-t001){ref-type="table"}.

2.5. Statistical Procedures {#sec2dot5-animals-10-01443}
---------------------------

All analyses were performed in R version 3.6.1 \[[@B87-animals-10-01443]\].

### 2.5.1. Comparative Analyses {#sec2dot5dot1-animals-10-01443}

Prior to hypothesis-testing we made two types of confound check. First, we investigated correlations between wild ecology predictor variables belonging to different hypotheses. Based on these results, we made further checks during hypothesis-testing by including correlated predictors as extra terms, to assess whether they altered interpretation of the focal predictor (in practice they never did, so hypothesis-testing models are reported without them). Our second confound check was to assess relationships between wild ecology predictor variables and the six species-typical husbandry variables shown in [Table 1](#animals-10-01443-t001){ref-type="table"}. Correlated husbandry variable(s) were included in final hypothesis-testing models.

To control for species non-independence due to shared ancestry, we used phylogenetic generalised least squares (PGLS) regressions \[[@B96-animals-10-01443],[@B97-animals-10-01443]\] for all continuous outcomes \[[@B98-animals-10-01443]\], and phylogenetic logistic regression models when ground-use was the outcome during between-predictor confound checks ("phylolm" package: \[[@B99-animals-10-01443]\]). For comparative confound checks, we used a consensus lemur phylogenetic tree (from \[[@B100-animals-10-01443]\]). Models were only run if data were available for ≥5 species, or ≥5 species per level for categorical variables. Potential outliers were assessed on graphs and PGLS diagnostic plots \[[@B98-animals-10-01443]\], and results reported with outliers if model assumptions were not violated, and they did not have a studentised phylogenetic residual \>±3 (e.g., \[[@B101-animals-10-01443]\]). Homogeneity of residuals was checked on diagnostic plots, and normality assessed using a Shapiro-Wilk normality test---where necessary transformations were applied to satisfy these. Pagel's Lambda (λ), a measure of phylogenetic signal \[[@B97-animals-10-01443],[@B102-animals-10-01443],[@B103-animals-10-01443]\], was estimated in PGLS models using maximum likelihood. The phylogenetic logistic regression equivalent of λ is alpha (α), which is reported in relevant models.

To account for phylogenetic uncertainty, all final hypothesis-testing models were performed over a tree block of 1000 alternative lemur phylogenetic trees (from \[[@B100-animals-10-01443]\]), and summarised as medians and 95% CIs. Results were considered significant at *p* \< 0.05 and all *p* values are two-tailed.

### 2.5.2. Epidemiological Analyses {#sec2dot5dot2-animals-10-01443}

For epidemiological analyses, we used mixed effects models in 'lmer4′ \[[@B104-animals-10-01443]\] and 'ordinal' \[[@B105-animals-10-01443]\]. To account for variance explained by non-independence between animals sharing enclosures, and enclosures within the same zoo, 'enclosure' was nested in 'zoo' as a random effect. Relative body mass was natural log-transformed to better meet the requirements of linear modelling, as were age and enclosure area because of their skew and hence leverage of high values. During hypothesis-testing, outliers with studentised residuals \>±3 (e.g., after: \[[@B101-animals-10-01443]\]) were removed (four animals) as their inclusion meant that models' residuals did not pass a Shapiro--Wilks normality test.

Being non-experimental data, the explanatory power of predictor variables (cf. [Table 2](#animals-10-01443-t002){ref-type="table"} in \[[@B106-animals-10-01443]\]) was assessed by comparing Akaike information criterion (AIC) scores of models with and without a focal predictor \[[@B107-animals-10-01443]\]. Models with lower AIC scores by ≥2 were judged as improvements \[[@B108-animals-10-01443]\]. If a focal predictor yielded a better fit to the data than a model without it, the model's coefficients were examined to interpret the relationship. For categorical predictor terms with \>2 factors, we spilt the dataset by that factor and re-ran the analysis to assess the relationship between the focal predictor and relative body mass.

We first checked for relationships between predictors using all data available for each variable. In all cases we used AIC scores to initially assess whether varying intercepts only or intercepts and slopes provided the better fit \[[@B109-animals-10-01443]\], with the random effects structure mentioned above. Linear mixed models were used when age and enclosure size were analysed as outcomes; generalised linear mixed models (binomial family and logit link function) were used for contraceptive use/neutered status and sex, and for when % fruit (gamma family and inverse link function) and enrichment score (Poisson family and log link function) were analysed as outcomes; and cumulative link mixed models (probit link function) \[[@B105-animals-10-01443]\] were used when the ordinally ranked enclosure type was an outcome. Relationships between variables were confirmed if the full model with the predictor term was a better fit to the data than null without it \[[@B110-animals-10-01443]\]. We ran additional models at the univariable stage of hypothesis-testing (below) to assess if correlated predictors affected interpretation of focal predictors, by including them as extra predictor terms (in practice they never did).

To enable model comparisons for hypothesis-testing, next we reduced our dataset to animals with complete data across all variables. We followed a similar model-building approach to KilBride et al. \[[@B58-animals-10-01443]\], although we evaluated interactions. Initially, we assessed univariable associations between predictors and our outcome using the same linear mixed models and procedure already described (here, varying intercepts and slopes never gave a better fit to the data than varying intercepts only). Predictors found to have a univariable relationship with relative body mass, as determined by AIC, were taken forward to the multivariable model-building stage.

Multivariable model building began with a model including all predictors with a univariable association with relative body mass. We used a forwards multiple regression technique and by making comparisons to the simplest version of this model (i.e., without any interactions), we checked for improvements in model fit by progressively including interactions, stopping when additions did not yield further improvements. This became our minimal adequate (baseline) model. Using this baseline model, we reassessed predictor variables without a univariable association with relative body mass, by sequentially added these into the baseline model, again checking for interactions, to see if they now improved model fit (after: \[[@B58-animals-10-01443],[@B111-animals-10-01443]\]). We continued until none of the remaining predictors resulted in improvements, yielding our final minimal adequate model.

3. Results {#sec3-animals-10-01443}
==========

3.1. Survey Response Rate and Descriptive Statistics {#sec3dot1-animals-10-01443}
----------------------------------------------------

Of 359 zoos contacted 135 responded, representing a response rate of 38%. Of the 13 species in our final comparative dataset, their sample sizes and species-typical relative body masses are shown in [Table 3](#animals-10-01443-t003){ref-type="table"}. [Figure 1](#animals-10-01443-f001){ref-type="fig"} provides a visual depiction of ring-tailed lemurs in different body conditions (note we did not have images of those in underweight and morbidly obese conditions).

Boxplots of individuals' relative body masses are shown in [Figure 2](#animals-10-01443-f002){ref-type="fig"}, grouped by species. Across our 13 species, 0.30% of animals were underweight (species-typical relative body mass \<0.75), 57.63% were healthy (0.75--1.25), 28.30% were overweight (1.25--1.5), 13.33% were obese (1.5--2) and 0.44% were morbidly obese (\>2).

Species-typical median relative body mass values (indicated by vertical lines on the boxes of [Figure 2](#animals-10-01443-f002){ref-type="fig"}) of eight species corresponded with the healthy category: aye ayes, red-collared lemurs, *E. collaris*, mongoose lemurs, red-bellied lemurs, Alaotran gentle lemurs, *Hapalemur alaotrensis*, Coquerel's sifaka, *Propithecus coquereli*, and black-and-white and red ruffed lemurs. Four species were classed as being typically overweight: crowned lemurs, blue-eyed black lemurs, black lemurs, *E. macaco*, and ring-tailed lemurs. One species, brown lemur, *E. fulvus*, was classed as obese.

3.2. Comparative Analyses: Results {#sec3dot2-animals-10-01443}
----------------------------------

[Figure 3](#animals-10-01443-f003){ref-type="fig"} shows the phylogenetic tree of the 13 species featured in our comparative analyses, with tip-points coloured according to species-typical median relative body mass values.

### 3.2.1. Comparative Analyses: Results of Confound Checks {#sec3dot2dot1-animals-10-01443}

All between-predictor confound checks are shown in [Appendix A](#app2-animals-10-01443){ref-type="app"} [Table A2](#animals-10-01443-t0A2){ref-type="table"}. Two sets of between-predictor check models were significant. Species that spend more time on the ground live in ranges with low annual rainfall (*t*~6~ = −3.64, R^2^ = 0.69, λ = 0.51, *p* = 0.01), and vice versa (*t*~6~ = −9.31, R^2^ = 0.94, λ = 0, *p* \< 0.001). Species more heavily predated also live in ranges with low annual rainfall (*t*~11~ = −3.16, R^2^ = 0.48, λ = 0.77, *p* = 0.01). Results of predictor-husbandry checks are also shown in [Appendix A](#app2-animals-10-01443){ref-type="app"} [Table A3](#animals-10-01443-t0A3){ref-type="table"}---no significant relationships between wild ecology predictors and species-typical husbandry were found.

### 3.2.2. Comparative Analyses: Results of Hypothesis-Testing {#sec3dot2dot2-animals-10-01443}

No aspect of wild ecology significantly predicted species-typical median relative body mass ([Table 4](#animals-10-01443-t004){ref-type="table"}); however, we did uncover one trend. Species that experience large between-year variation in rainfall, and thus greater food resource unpredictability, tend to have larger species-typical relative body masses ([Figure 4](#animals-10-01443-f004){ref-type="fig"}; as the 95% CI values are identical \[see [Table 4](#animals-10-01443-t004){ref-type="table"}\], here we report median values for each parameter: *t*~11~ = 2.04, R^2^ = 0.27, λ \< 0.01, *p* = 0.07). To illustrate, black lemurs, crowned lemurs, and ring-tailed lemurs are examples of species that experience relatively large between-year variation in rainfall and, correspondingly, whose representatives are typically overweight in captivity.

3.3. Epidemiological Analyses: Results {#sec3dot3-animals-10-01443}
--------------------------------------

### 3.3.1. Epidemiological Analyses: Between-Predictor Checks Results {#sec3dot3dot1-animals-10-01443}

Results of predictor-predictor checks featuring the four best-sampled species included in the comparative analyses are shown and described in [Appendix A](#app2-animals-10-01443){ref-type="app"} [Table A4](#animals-10-01443-t0A4){ref-type="table"}.

### 3.3.2. Epidemiological Analyses: Univariable Results {#sec3dot3dot2-animals-10-01443}

The reduced (complete cases) epidemiological dataset contained data for 256 animals. Details of predictors with a univariable association with relative body mass are shown in [Appendix A](#app2-animals-10-01443){ref-type="app"} [Table A5](#animals-10-01443-t0A5){ref-type="table"}. In sum, models including species, sex, and age as predictors were found to individually yield better fits to the data than the null model. Including correlated predictors ([Section 2.5.2](#sec2dot5dot2-animals-10-01443){ref-type="sec"} and [Section 3.3.1](#sec3dot3dot1-animals-10-01443){ref-type="sec"}) did not change the interpretation of focal predictors in relevant models.

### 3.3.3. Epidemiological Analyses: Multivariable Results {#sec3dot3dot3-animals-10-01443}

The baseline model with the three predictors found to have a univariable association with relative body mass, is shown in [Appendix A](#app2-animals-10-01443){ref-type="app"} [Table A6](#animals-10-01443-t0A6){ref-type="table"}. Species × sex was a better fit to the data than a simpler model without it (AIC = −211.93 v −216.09). Male mongoose lemurs have smaller mean relative body masses, corresponding with healthy, than females, which are typically overweight (*t*~6~ = −2.76, *p* = 0.03); whereas male ring-tailed lemurs have larger mean values than females, but both sexes are typically overweight (*t*~127~ = 3.15, *p* \< 0.01). Species and age both had main effects. Black-and-white ruffed lemurs had smaller mean relative body mass values than ring-tailed lemurs (*t*~146~ = −5.11, *p* \< 0.001) and red ruffed lemurs (*t*~4~ = −3.94, *p* = 0.02), and relative body mass increased with age (*t*~157~ = 2.09, *p* = 0.04).

Sequentially adding predictors without a univariable relationship with the outcome into the baseline model yielded the following (see [Table 5](#animals-10-01443-t005){ref-type="table"}). Males in enclosures featuring only fixed climbing structures were on average obese, whereas those with flexible and fixed structures were typically overweight (*t*~157~ = 3.04, *p* \< 0.01). Females weighed in the winter had larger mean relative body masses, and were normally overweight, than those weighed in spring (*t*~178~ = 3.30, *p* \< 0.01) and summer (*t*~25~ = 3.70, *p* \< 0.01) which fell within the healthy range. See [Figure 5](#animals-10-01443-f005){ref-type="fig"} for visual representations of the multivariable results.

4. Discussion {#sec4-animals-10-01443}
=============

Our comparative study found no support for our hypotheses relating susceptibility to captive weight gain to arboreality, predation risk, or low productivity. However, given our sample sizes were small, likely reducing statistical power (*n* \< 20; cf. \[[@B112-animals-10-01443]\]), we discuss the trend consistent with one hypothesis. Lemurs for which native ranges have large between-year variation in annual rainfall, and who are thus assumed to be adapted to greater food resource unpredictability, tended to have larger relative body masses. They may therefore have "thrifty" adaptations that represent an ecological risk factor for susceptibly to captive weight gain, pre-disposing such species to obesity and its related health problems.

Our epidemiological analyses revealed four demographic and one environmental risk factor for large relative body masses across four species. We confirmed species differences in average relative body masses; found relative body mass to increase with age; and that for two species there was a sex effect. Climbing structure provision influenced relative body masses across males of all species. Males whose climbing structures were fixed had larger relative body masses than those housed with some flexible climbing structures too. Additionally, for females there was a seasonal effect. Next, we discuss findings from each study in turn. Based on these, we then make recommendations for zoos managing these species, describe limitations of our study, and detail ideas for future research.

4.1. Comparative Study {#sec4dot1-animals-10-01443}
----------------------

Regarding species differences our results show some agreement with two previous multi-species studies (\[[@B18-animals-10-01443]\]; note we used Terranova and Coffman's \[[@B19-animals-10-01443]\] body mass values to calculate the typical body condition of their lemurs ourselves), and our epidemiological analyses likewise confirmed species differences ([Table 5](#animals-10-01443-t005){ref-type="table"}). Thus, of species common to all three studies, red-bellied lemurs and black-and-white ruffed lemurs emerged as usually healthy, whereas blue-eyed black lemurs and black lemurs are typically overweight. Mongoose lemurs provide the one instance of disagreement, being typically healthy in our study but overweight in the previous two. This between-study difference might reflect differences in data sources and measurement methods, or recent improvements in captive diet and husbandry for this species (e.g., based on recommendations in: \[[@B38-animals-10-01443]\]).

How might "thrifty" adaptations to unpredictable wild food resources potentially increase susceptibility to captive weight gain in lemurs? Lemurs have various attributes proposed to facilitate survival in their native, and unpredictable, environment. These include increased resting, food-switching (e.g., from fruit to leaves and flowers), reduced basal metabolic rates and hibernation/torpor during lean periods, as well as seasonal breeding and weaning synchrony \[[@B25-animals-10-01443]\]. Therefore, lemurs are likely adapted to take advantage of plentiful food when available (cf. \[[@B25-animals-10-01443]\]) as, presumably, is the case in captive environments. Possibly, if food is delivered in two or three bouts (meals) rather than encountered slowly over the course of a day as per wild food, this might promote binge-eating. If captive food also tends to be energy-rich (cf. \[[@B38-animals-10-01443]\]) then our result might simply reflect a constant state of positive energy imbalance and subsequent weight gain (sensu \[[@B13-animals-10-01443]\]).

As unpredictability is key to the potential ecological risk factor identified here, paradoxically, for some species perhaps food unpredictability is being inadvertently signalled by captive feeding regimes, leading to animals going into "thrifty" mode (e.g., increased resting, reduced basal metabolic rate: \[[@B25-animals-10-01443]\]) ultimately resulting in weight gain. Food unpredictability, rather than greater access to energy-rich food *per se*, is proposed to help explain why people from lower socioeconomic positions in wealthy (but not poorer) countries are at most risk of being overweight or obese \[[@B113-animals-10-01443],[@B114-animals-10-01443]\]. For these people, it is the combination of food being sometimes perceived as unpredictable and, when accessible, being energy-rich that might explain this effect \[[@B114-animals-10-01443]\]. A similar scenario could be relevant for some captive lemurs too. Even if regularly provisioned daily, food provided in a small number of discrete meals rather than being available *ad libitum* might signal unpredictability. This may be especially so for animals in which eating just two or three times per day is likely at odds---a mismatch (cf. \[[@B115-animals-10-01443],[@B116-animals-10-01443]\])---with their foraging and eating behaviour during times of plenty in the wild. In humans, attempts to lose weight via episodes of restrictive dieting are usually counterproductive, ultimately resulting in weight gain because such dieting likely cues food insecurity, triggering fat-storage mechanisms when food becomes available \[[@B114-animals-10-01443],[@B117-animals-10-01443],[@B118-animals-10-01443],[@B119-animals-10-01443],[@B120-animals-10-01443]\]. For some captive lemurs, then, "thrifty" adaptations combined with being provisioned with energy-rich food (cf. \[[@B38-animals-10-01443],[@B39-animals-10-01443]\]) in a small number of daily bouts could lead to weight gain.

4.2. Epidemiological Study {#sec4dot2-animals-10-01443}
--------------------------

We found a species × sex effect in our epidemiological analyses. Because sex-specific wild means were not available for every species in our comparative analyses (see \[[@B66-animals-10-01443]\]), we used species means for relative body mass calculations. However, the sexes of some species do differ in size, and while these differences have not been found to be statistically significant \[[@B18-animals-10-01443]\] this could explain our result. Another, more speculative explanation relates to dominance hierarchy. In other taxa, subordinate individuals maintain higher fat reserves than dominants, under scenarios of relatively low energy cost and high food availability (e.g., \[[@B121-animals-10-01443]\]) like, presumably, captivity. Under such conditions, the trade-off between the risks of death from starvation and death from predation \[[@B122-animals-10-01443]\] differ according to an animal's place in the hierarchy \[[@B123-animals-10-01443],[@B124-animals-10-01443]\]. Because they do not have feeding priority, subordinates carry greater fat stores as their starvation risk outweighs the cost of carrying extra weight, e.g., increased risk of death by predation \[[@B49-animals-10-01443],[@B123-animals-10-01443]\]. Logically, as most lemurs, including our four epidemiological study species, are female-dominant \[[@B25-animals-10-01443],[@B125-animals-10-01443],[@B126-animals-10-01443],[@B127-animals-10-01443]\] we might expect captive males to be the heavier sex, as we found here for ring-tailed lemurs. But, then, why did we find the *opposite* for mongoose lemurs? A key assumption of the scenario described above, is that subordinate animals should only carry greater fat stores under specific conditions. Under other scenarios, like high starvation risk, dominants are expected to carry greater fat stores \[[@B124-animals-10-01443]\]. Potentially then, if species differ in their perception of starvation or predation risk, or if food presentation and/or quality is signalled as being unpredictable, one would expect dominant animals, i.e., females, to carry greater fat stores, especially if they get first choice of likely energy-rich captive food (sensu \[[@B123-animals-10-01443],[@B124-animals-10-01443],[@B128-animals-10-01443]\]).

We also revealed an environmental risk factor for large relative body masses in males: being housed with fixed climbing structures only (as opposed to having some flexible structures too). Opportunities for regular exercise/activity may affect weight via an animal's energy balance, and one can reasonably assume climbing fixed structures requires less physical effort, and therefore less energy expenditure, than does climbing flexible structures (sensu \[[@B13-animals-10-01443]\]). Exercise-induced weight loss is reported to vary between the sexes in humans: males expend relatively more energy than females, being larger-bodied with larger total daily energy expenditure \[[@B129-animals-10-01443]\]. A similar effect could explain why male lemurs, and not females, seem more affected by exercise opportunities (but see: \[[@B18-animals-10-01443]\]). Alternatively, assuming climbing structures are more preferred than the ground, dominant (i.e., female) (cf. \[[@B127-animals-10-01443]\]) exclusion of subordinates (i.e., males) from the preferred substrate might explain our result, if fixed climbing structures are easier for dominant animals to monopolise. A similar effect might also explain why the lowest ranking squirrel monkeys, *Saimiri sciureus*, in Marriott and Meyers \[[@B130-animals-10-01443]\] mainly used the ground and feeding areas, instead of suspended logs like their higher-ranking cage mates. Alternatively still, our result could reflect other sex differences in climbing structure use and preference (e.g., female chimpanzees, *Pan troglodytes*, used elevated areas more than males, perhaps to avoid the males' displays: \[[@B131-animals-10-01443]\]).

Females' relative body masses, but not males', differed seasonally. While known pregnant females were excluded from analyses (see [Section 2.2](#sec2dot2-animals-10-01443){ref-type="sec"}), females with undiagnosed pregnancy could be in our final dataset. Breeding is highly seasonal and often synchronised for wild lemurs (e.g., \[[@B12-animals-10-01443],[@B25-animals-10-01443]\]), and captive lemurs also maintain strongly seasonal breeding patterns \[[@B132-animals-10-01443],[@B133-animals-10-01443]\]. Thus, the seasonal pattern in relative body mass reported here could reflect pregnancies. Breeding is usually carefully managed in captivity, however. Alternatively, then, our seasonal effect might represent the close association between body condition and (preparation) for reproduction \[[@B134-animals-10-01443],[@B135-animals-10-01443]\]. Even if unable to breed, fat stores of females from seasonal breeding species might nevertheless fluctuate in preparation for reproduction (sensu "capital breeding" in: \[[@B136-animals-10-01443]\]). Our result might therefore indicate key times of the year during which females might be at increased risk of weight gain.

Finally, older animals had larger relative body mass values, in good agreement with other studies of weight-gain and age in non-human primates \[[@B137-animals-10-01443]\] and humans (inferred using body mass index: \[[@B138-animals-10-01443]\]). Although very old animals might instead be lighter than those younger, due to effects of chronic illness and other age-related complications (sensu \[[@B69-animals-10-01443]\]), there was no hint of a non-linear relationship in our data. Examining [Figure 5](#animals-10-01443-f005){ref-type="fig"} suggests the older half of our sample are more likely to be overweight, but that the regression line does not cross the obesity threshold. In other words, increasing age does not seem to pose a strong risk for potential obesity. However, note that most of our animals were relatively young: only 14% (37/256) were ≥50% of their species-specific maximum recorded captive lifespan, limiting our ability to extrapolate beyond the current dataset. Therefore, further research is required to determine the health relevance, if any, of the age effect observed in our study.

4.3. Recommendations for Zoos {#sec4dot3-animals-10-01443}
-----------------------------

Our results tentatively identify which species---the "thriftier" ones---might benefit most from weight management recommendations (e.g., in: \[[@B18-animals-10-01443],[@B38-animals-10-01443],[@B39-animals-10-01443]\]). Captive diets should be carefully designed to avoid overfeeding calorie-rich diets \[[@B38-animals-10-01443],[@B39-animals-10-01443]\], especially for species of concern in [Figure 3](#animals-10-01443-f003){ref-type="fig"}, e.g., brown lemurs and blue-eyed black lemurs. Limiting/avoiding commercial fruit also seems sensible (also see: \[[@B18-animals-10-01443],[@B38-animals-10-01443],[@B39-animals-10-01443]\]), especially as many wild lemurs readily switch from a primarily fruit-based diet to e.g., leaves and flowers \[[@B25-animals-10-01443]\]. Factoring appropriate exercise opportunities into enclosures could also help manage weight (also see: \[[@B18-animals-10-01443]\]), e.g., by providing a mixture of flexible and fixed climbing structures (see [Section 4.2](#sec4dot2-animals-10-01443){ref-type="sec"}); encouraging movement by placing food and other resources such as enrichment at various heights and distances around the enclosure; and necessitating climbing and distances to be covered into enclosures at the initial design stage.

Implicit in our comparative result is that the "thriftiness" enabling some species to succeed in the wild renders the same species prone to captive weight gain. Based on this, we predict which species outside our dataset may likewise be at risk of captive weight gain. Thus, white-fronted brown lemurs, *Eulemur albifrons*, might be at risk because their geographic range is nearby---and so is likely similarly unpredictable---to some of our typically prone species (e.g., crowned lemurs; brown lemurs). Furthermore, one might predict that other taxa with similarly unpredictable wild environments may also be prone to weight gain in captivity.

4.4. Limitations {#sec4dot4-animals-10-01443}
----------------

Our sampling method of captive animals is non-random, as we recruited Species360 member zoos. Using survey responses from different people likely introduced noise into our data (see \[[@B63-animals-10-01443]\]). Being observational rather experimental also means there may be other unmeasured variables not considered here, e.g., calorie intake (sensu \[[@B13-animals-10-01443]\]), which could affect our outcome and/or be confounded with our predictors. As lemur taxonomy is ever being updated, we cannot rule out potential species (mis)identification especially for *Eulemur* animals \[[@B41-animals-10-01443]\] within both our captive sample and animals used to calculate species-typical wild means \[[@B66-animals-10-01443]\]. Our study has unequal representation of families (just one Indriiae and no Cheirogaleidae representatives, the latter for reasons given in [Section 2.2](#sec2dot2-animals-10-01443){ref-type="sec"}), although it is likely representative of lemur species common in captivity (cf. \[[@B139-animals-10-01443]\]). Seasonal variation was likely represented in some species-typical wild mean weights (cf. \[[@B66-animals-10-01443]\]), but not all (e.g., \[[@B140-animals-10-01443],[@B141-animals-10-01443]\]), which might reduce the representativeness of these as 'gold standards' for body mass. We also assumed our wild ecology predictor variables accurately reflect evolved adaptations, e.g., that between-year annual rainfall CV indirectly reflects adaptations to wild food resource unpredictability: an assumption that needs testing. Finally, as a minimum of 20 species is recommended for comparative analyses for acceptable power and Type I error rates \[[@B112-animals-10-01443]\], these analyses are likely underpowered and should be treated cautiously. This limitation is unlikely to be readily solved; further species are unlikely to be taken into captivity and, if they are (e.g., to avoid extinction), will probably be in small numbers. Nonetheless, these limitations aside, the nature of our study yields results with predictive values as indicators of where problems arise from and how they may be addressed.

4.5. Future Research {#sec4dot5-animals-10-01443}
--------------------

Relative body mass as a proxy for body condition and the accuracy of its obesity threshold now requires validation, e.g., using lifespan (e.g., from \[[@B69-animals-10-01443]\]) and/or the prevalence of reproductive problems and other comorbidities of obesity (e.g., from records in ZIMS \[[@B142-animals-10-01443]\] or from \[[@B143-animals-10-01443]\]; but see Mellor et al. \[[@B63-animals-10-01443]\]). If the current threshold of ≥1.5 does correspond with obesity, then such animals should have reduced lifespan and higher prevalence of comorbidities. Alternatively, it might be that "thrifty" species are able to carry relatively greater fat-stores safely, without determinant to health---indeed, this would be adaptive, for wild animals at least. If this were the case, then using obesity thresholds from other taxa (20--25% over the ideal in dogs, cats and humans: \[[@B7-animals-10-01443],[@B13-animals-10-01443]\]) would be inappropriate, and even Taylor et al. \[[@B18-animals-10-01443]\]'s ≥ 50% of species-typical wild mean might prove inaccurate.

Future research might also examine relationships between the potential ecological risk factor identified here, weight gain and the captive feeding environment. Are animals from "thrifty" species simply easier to overfeed? Do some animals binge-eat? Or is food unpredictability accidentally being signalled by food provision schemes (see [Section 4.1](#sec4dot1-animals-10-01443){ref-type="sec"})? Lastly, a future comparative study might examine whether dominance hierarchy effects do promote fat storage in subordinate individuals, and whether this interacts with the potential ecological risk factor uncovered here (i.e., to assess whether "thriftiness" is more advantageous for subordinates than dominants). If so, one would predict the most despotic species should have the greatest weight difference between subordinates and dominants. This would have management implications; for example, in such species, providing food in multiple sites so it cannot be readily monopolised by dominant individuals.

Finally, the cross-species approach used here clearly highlights species whose typical captive body conditions are concerning, because they deviate the most from their wild norms. Brown lemurs (but see [Section 4.4](#sec4dot4-animals-10-01443){ref-type="sec"}), blue-eyed black lemurs, crowned lemurs, black lemurs, and ring-tailed lemurs (see [Figure 3](#animals-10-01443-f003){ref-type="fig"}) might therefore be priority species for research effort.

5. Conclusions {#sec5-animals-10-01443}
==============

By using two complementary methods to address our research question, we uncovered one potential ecological, four demographic, and one environmental risk factor for susceptibility to weight gain in captive lemurs. Using phylogenetic comparative methods, we found species that experience large between-year variation in annual rainfall, and thus greater food resource unpredictability, tended to have larger species-typical relative body masses. Using an epidemiological approach, we also report that for some species the sexes were unequally affected: female mongoose lemurs were heavier than males, but male ring-tailed lemurs were heavier than females. Relative body mass also increases with age, though it is unclear if increasing age *per se* is a risk factor for being overweight or obese. Climbing structure provision also affected males: those housed with only fixed climbing structures had larger relative body masses than those housed with a mixture of flexible and fixed structures. Finally, for females there was a seasonal effect, as those weighed in winter had larger relative body masses than those weighed in spring or summer. Based on these results, we were able to make targeted practical recommendations to help address weight concerns. The validation of relative body mass as a proxy for body condition is now required.

Thanks to Christie Eddie and the rest of the AZA Prosimian TAG; Amy Plowman and the rest of the BIAZA Research Committee; and Sarah Zehr. For contributing data, many thanks to: Adelaide Zoo, Al Bustan Zoological Centre, United Arab Emirates, Alameda Wildlife Park, Alligator Farm Zoological Park, Aquazoo Friesland, Audubon Zoo, Barcelona Zoo, Birmingham Zoo, Cameron Park, Cango Wildlife Ranch, Central Florida Zoo, Cheyenne Mountain Zoo, Cleveland Metroparks Zoo, Duke Lemur Center, Furuvik Zoo, Givskud Zoo---Zootopia, Greater Vancouver Zoo, Hannover Zoo, Hemsley Conservation Centre, Honolulu Zoo, Kaliningrad Zoo, Kosice Zoo, Lake Superior Zoological Society, Lakeland Wildlife Oasis, Lehigh Valley Zoo, Lemur Conservation Foundation, Little Rock Zoo, Minsk Zoo, Nashville Zoo, National Zoological Park, Nikolaev Zoo, Omaha's Henry Doorly Zoo and Aquarium, Paignton Zoo, Parc Animalier D'auvergne, Parc Merveilleux, Parco Faunistico Valcobra, Parken Zoo, Peak Wildlife Park, Perm Zoo, Point Defiance Zoo & Aquarium, Randers Regnskov, Rolling Hills Zoo, Rotterdam Zoo, Santa Fe College Teaching Zoo, Śląski Ogród Zoologiczny in Chorzów, The Cincinnati Zoo and Botanical Garden, Twycross Zoo and the East Midlands Zoological Society, Wellington Zoo Trust, Zoo Bratislava, Zoo de Lyon, Zoo Landau, Zoo Vienna (Tiergarten Schönbrunn), Zoo Wuppertal, Zoologischer Garten Schwerin, and all others that contributed.

The following are available online at <https://www.mdpi.com/2076-2615/10/8/1443/s1>, File S1: copy of survey questions used to collect data on captive lemurs' body masses and corresponding living conditions.

###### 

Click here for additional data file.

Original concept: G.J.M.; planning: all authors; analysis and first draft: E.L.M.; review and editing: E.L.M.; I.C.C.; C.S.; G.J.M. and M.M. All authors have read and agreed to the published version of the manuscript.

E.L.M. was funded by a University of Bristol PhD Scholarship for this work.

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

animals-10-01443-t0A1_Table A1

###### 

Ages from which each species is considered adult (see [Section 2.2](#sec2dot2-animals-10-01443){ref-type="sec"} for rationale). Values are rounded to the nearest whole number.

  ------------------------------------------------------------------------------------------------------------------------------
  Species Name                     Common Name                    Dam Minimum Age at Conception (Days)   Adult From: (Days) \*
  -------------------------------- ------------------------------ -------------------------------------- -----------------------
  *Daubentonia madagascariensis*   Aye aye                        1540                                   3081

  *Eulemur collaris*               Red-collared lemur             599                                    1197

  *Eulemur coronatus*              Crowned lemur                  624                                    1248

  *Eulemur flavifrons*             Blue-eyed black lemur          580                                    1161

  *Eulemur fulvus*                 Brown lemur                    507                                    1015

  *Eulemur macaco*                 Black lemur                    540                                    1080

  *Eulemur mongoz*                 Mongoose lemur                 650                                    1299

  *Eulemur rubriventer*            Red-bellied lemur              650                                    1299

  *Hapalemur alaotrensis*\         Alaotran gentle lemur          558                                    1117
  *(H. griseus)*                                                                                         

  *Lemur catta*                    Ring-tailed lemur              489                                    978

  *Propithecus coquereli*          Coquerel's sifaka              964                                    1927

  *Varecia rubra*                  Red ruffed lemur               610                                    1219

  *Varecia variegata*              Black-and-white ruffed lemur   588                                    1175
  ------------------------------------------------------------------------------------------------------------------------------

\* Twice the minimum dam age at conception (following: \[[@B69-animals-10-01443]\]). NB: As the Duke Lemur Center only holds one of the black-and-white ruffed lemur sub-species, and the sub-species are considered very similar, we used this (1175 days) as the cut off for both sub-species in our dataset. One species, Alaotran gentle lemur, *Hapalemur alaotrensis,* is not held by the Duke Lemur Center. We therefore used values from its closest relative, *Hapalemur griseus* (1117 days), as indicated in parentheses.

animals-10-01443-t0A2_Table A2

###### 

Results of models investigating potential relationships among predictor variables belonging to different hypotheses. PGLS models were used for continuous outcomes, and phylogenetic logistic regression was used to analyse ground use as an outcome. CV = coefficient of variation. Dashes indicate that the analysis could not be performed (because neither PGLS nor phylogenetic logistic regression models are appropriate for ordinal outcome variables like predation risk). Results in bold were considered significant at *p* \< 0.05, and trends (*p* \< 0.10) are italicised. All *p* values are two-tailed.

  ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
  Hypothesis:                    \(i\) Low Productivity                        \(ii\) Unpredictability           \(iv\) Arboreality      \(v\) Predation Risk                                                 
  ------------------------------ --------------------------------------------- --------------------------------- ----------------------- ------------------------------------ ------------------------------- ----
  **Annual rainfall (mm)**                                                     *t*~11~ = 0.25, R^2^ = 0.01,\     Z = −1.25, α = 0.04,\   ***t*~6~ = −9.31, R^2^ = 0.94,**\    *t*~11~ = 0.34, R^2^ = 0.01,\   \-
                                                                               λ = 0, *p* = 0.81                 *p* = 0.21              **λ = 0, *p* \< 0.001**              λ = 0, *p* = 0.74               

  **Between-year rainfall CV**   *t*~11~ = 0.40, R^2^ = 0.01,\                                                   Z= 0.04, α = 0.08,\     *t*~6~ = −0.23, R^2^ = 0.01,\        *t*~10~ = 1.11, R^2^ = 0.11,\   \-
                                 λ = 0.90, *p* = 0.70                                                            *p* = 0.97              λ = 0, *p* = 0.83                    λ = 0, *p* = 0.29 ^a^           

  **Ground use**\                *t*~10~ = −1.66, R^2^ = 0.22,\                *t*~10~ = 0.02, R^2^ \< 0.001,\                                                                                                \-
  **(rare v some)**              λ = 0.65, *p* = 0.13                          λ = 0, *p* = 0.99                                                                                                              

  **Ground time (% scans)**      ***t*** **~6~** **= −3.64, R^2^ = 0.69,**\    *t*~6~ = −0.23, R^2^ = 0.01,\                                                                                                  \-
                                 **λ** **= 0.51,** ***p*** **= 0.01**          λ = 0, *p* = 0.83                                                                                                              

  **Relative tail length \***    *t*~11~ = 0.29, R^2^ = 0.01,\                 *t*~10~ = 1.47, R^2^ = 0.18,\                                                                                                  \-
                                 λ = 0.91, *p* = 0.78                          λ = 0, *p* = 0.17 ^a^                                                                                                          

  **Predation score (0--3)**     ***t*** **~11~** **= −3.16, R^2^ = 0.48,**\   *t*~11~ = 0.52, R^2^ = 0.02,\     Z \< 0.01, α = 0.10,\   *t* ~6~ *= 2.19, R* ^2^ *= 0.45,*\   *t*~11~ = −0.01,\               
                                 **λ** **= 0.77,** ***p*** **= 0.01**          λ = 0, *p* = 0.61                 *p* = 0.99              *λ = 0, p = 0.07*                    R^2^ \<0.001, λ = 0,\           
                                                                                                                                                                              *p* = 0.99                      
  ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

\* This is the length of the tail as a proportion of the entire length of the animal, from head to tail tip. Outlier removed: ^a^ *Eulemur fulvus.*

animals-10-01443-t0A3_Table A3

###### 

Results of PGLS models investigating potential relationships between species-typical husbandry and wild ecology predictor variables. In all models the husbandry variables are outcome variables. Dashes indicate an analysis could not be performed due to too few species per level for ground use (see [Section 2.5.1](#sec2dot5dot1-animals-10-01443){ref-type="sec"}). CV = coefficient of variation. Results in bold were considered significant at *p* \< 0.05, and trends (*p* \< 0.10) are italicised. All *p* values are two-tailed.

  -----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
  Hypothesis:                   Wild Ecology Predictor                          Median Enclosure Area (m^2^)                    Median Enrichment Score (0--25)   Median Proportion of Diet Fruit   Proportion Given Contraception   Proportion Housed Indoors        Proportion with Flexible Climbing Materials
  ----------------------------- ----------------------------------------------- ----------------------------------------------- --------------------------------- --------------------------------- -------------------------------- -------------------------------- ---------------------------------------------
  \(i\) low productivity        **Total annual rainfall (mm)**                  *t*~11~ = 0.23,\                                *t*~11~ = −1.27,\                 *t*~11~ = 0.68,\                  *t*~10~ = −0.61,\                *t*~11~ = 0.73,\                 *t*~11~ = −0.32,\
                                                                                R^2^ = 0.01, λ \< 0.01, *p* = 0.82              R^2^ = 0.13, λ = 0,\              R^2^ = 0.04,\                     R^2^ = 0.04, λ = 0,\             R^2^ = 0.05, λ = 1,\             R^2^ = 0.01,\
                                                                                                                                *p* = 0.23                        λ = 0, *p* = 0.51                 *p* = 0.56                       *p* = 0.48                       λ = 0, *p* = 0.76

  \(ii\) unpredictability       **Between-year rainfall CV**                    *t*~11~ = 1.64,\                                *t*~11~ \< 0.001,\                *t*~11~ = 0.26,\                  *t*~10~ = 0.39,\                 *t*~11~ = 1.37,\                 *t*~11~ = 1.27, R^2^ = 0.13,\
                                                                                R^2^ = 0.20, λ = 1,\                            R^2^ = 0.08, λ = 0,\              R^2^ = 0.01,\                     R^2^ = 0.02,\                    R^2^ = 0.15,\                    λ = 0, *p* = 0.23
                                                                                *p* = 0.13                                      *p* = 0.36                        λ = 0, *p* = 0.80                 λ = 0, *p* = 0.70                λ = 1, *p* = 0.20                

  \(iii\) arboreality           **Ground use**\                                 *t*~10~ = 0.35,\                                *t*~10~ = −0.34,\                 *t~10~ = −2.16,*\                 \-                               *t*~10~ = 0.28,\                 *t*~10~ = 0.33, R^2^ = 0.01,\
                                **(some v rare)**                               R^2^ = 0.01, λ = 1,\                            R^2^ = 0.01, λ = 0,\              *R^2^ = 0.32,* *λ = 0,*\                                           R^2^ = 0.01,\                    λ = 0, *p* = 0.75
                                                                                *p* = 0.73                                      *p* = 0.74 ^a^                    *p = 0.06*                                                         λ = 0.97, *p* = 0.79             

  **Ground time**\              *t*~6~ = 0.35, R^2^ = 0.01, λ = 1, *p* = 0.73   *t*~6~ = 0.46, R^2^ = 0.03, λ = 0, *p* = 0.66   *t*~6~ = −1.05,\                  *t*~6~ = 0.26,\                   *t*~6~ = −0.15,\                 *t*~6~ = −0.01, R^2^ \< 0.01,\   
  **(% scans)**                                                                                                                 R^2^ = 0.15,\                     R^2^ = 0.01,\                     R^2^ \< 0.01,\                   λ = 0, *p* = 0.99                
                                                                                                                                λ = 0.06, *p* = 0.34              λ = 0, *p* = 0.81                 λ = 0, *p* = 0.88                                                 

  **Relative tail length \***   *t*~11~ = −0.46,\                               *t*~11~ = −1.64,\                               *t*~11~ = 0.28,\                  *t*~10~ = −1.27,\                 *t*~11~ = 0.38,\                 *t*~11~ = −0.28,\                
                                R^2^ = 0.02, λ = 0,\                            R^2^ = 0.20, λ = 0,\                            R^2^ = 0.01,\                     R^2^ = 0.14, λ = 0,\              R^2^ = 0.01,\                    R^2^ = 0.01,\                    
                                *p* = 0.66                                      *p* = 0.13                                      λ = 0, *p* = 0.79                 *p* = 0.23                        λ = 0.98, *p* = 0.71             λ = 0, *p* = 0.79                

  \(iv\) predation risk         **Predation score**\                            *t*~11~ = −0.47,\                               *t*~11~ = 1.36,\                  *t*~11~ = 0.27,\                  *t*~10~ = 0.65,\                 *t*~11~ = −0.22,\                *t*~11~ = −0.26,\
                                **(0--3)**                                      R^2^ = 0.02, λ = 0,\                            R^2^ = 0.15, λ = 0,\              R^2^ = 0.01,\                     R^2^ = 0.04,\                    R^2^ \< 0.01,\                   R^2^ = 0.01,\
                                                                                *p* = 0.66                                      *p* = 0.20                        λ = 0, *p* = 0.80                 λ = 0, *p* = 0.53                λ = 0.98, *p* = 0.83             λ = 0, *p* = 0.80
  -----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

\* This is the proportion of the length of the tail compared to the entire length of the animal, from head to tail tip. ^a^ Note that the residuals from this model did not pass a Shapiro-Wilks normality test, despite numerous transformation attempts, so this result should be treated cautiously.

animals-10-01443-t0A4_Table A4

###### 

Results of between-predictor checks for the epidemiological analyses. Note that species and season, being categorical with \>2 levels and nominal, were not analysable as outcome variables. In all cases, we first checked if varying intercepts only or intercepts and slopes were the best fit. In practice, varying intercepts and slopes were rarely a better fit than intercepts only, and these are indicated with an asterisk. Model comparisons were then made between the full model including the focal predictor and the null, using AIC (see [Section 2.5.2](#sec2dot5dot2-animals-10-01443){ref-type="sec"}). Associations between variables were confirmed if the full model yielded an improved fit, i.e., a lower AIC score by ≥ 2. The AIC of the null model is reported first. Full models found to be a better fit to the data than the null are shown in bold. For these, we assessed the nature of the relationship between the predictor and the outcome using the model's coefficients, and these are described in the text immediately below the table.

  ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
  Outcome:\                         Sex                           Age                        Contraceptive Use/Neutered    Enclosure Area                Enclosure Type           Climbing Structures   Enrichment Score            \% Fruit
  Predictor:                                                                                                                                                                                                                        
  --------------------------------- ----------------------------- -------------------------- ----------------------------- ----------------------------- ------------------------ --------------------- --------------------------- ------------------------
  **Species**\                      AIC: 708.84 v 711.28          **AIC: 640.79 v 637.73**   AIC: 423.59 v 425.24          **AIC: −405.81 v −434.91**    **AIC: 94.61 v 86.60**   AIC: 40.63 v 45.88    AIC: 2775.40 v 2780.60      AIC: 3376.50 v 3377.90
  Ref: Mongoose lemur                                                                                                                                                                                                               

  **Sex**\                                                        AIC: 640.79 v 641.64       AIC: 300.71 v 302.57 \*       AIC: −413.46 v −412.53 \*     AIC: 94.61 v 96.18       AIC: 40.63 v 46.59    AIC: 2775.40 v 2777.30      AIC: 3376.50 v 3377.90
  Ref: Female                                                                                                                                                                                                                       

  **Age**\                          AIC: 484.64 v 485.31                                     AIC: 296.92 v 298.27          AIC: −190.50 v −191.78        AIC: 48.33 v 50.22       AIC: 22.46 v 24.46    AIC: 1870.70 v 1872.70      AIC: 2468.80 v 2467.60
  (days)                                                                                                                                                                                                                            

  **Season**\                       **AIC: 708.84 v 699.01**      AIC: 640.79 v 640.31       AIC: 423.59 v 426.35          **AIC: −405.81 v −445.78**    AIC: 94.61 v 99.45       AIC: 40.63 v 46.58    AIC: 2775.40 v 2777.70      AIC: 3376.50 v 3379.90
  Ref: Autumn                                                                                                                                                                                                                       

  **Contraceptive use/neutered**\   **AIC: 557.99 v 551.91 \***   AIC: 531.78 v 532.93                                     AIC: −451.96 v −452.54        AIC: 94.55 v 96.55       AIC: 40.58 v 42.23    AIC: 2532.60 v 2534.40      AIC: 2800.20 v 2800.40
  Ref: No                                                                                                                                                                                                                           

  **Enclosure area** (m^2^)         AIC: 620.58 v 622.55          AIC: 569.46 v 571.32       AIC: 358.18 v 359.13                                        AIC: 61.79 v 62.44       AIC: 36.79 v 35.71    AIC: 2385.00 v 2386.80      AIC: 2759.60 v 2761.40

  **Enclosure type**\               AIC: 708.84 v 710.82          AIC: 640.79 v 642.23       AIC: 423.59 v 424.83          **AIC: −405.81 v −408.53**                             AIC: 40.63 v 44.60    **AIC: 2775.40 v 2773.3**   AIC: 3376.50 v 3377.80
  Ref: Indoor and outdoor                                                                                                                                                                                                           

  **Climbing structures**\          AIC: 708.84 v 710.61          AIC: 640.79 v 642.79       AIC: 423.59 v 425.58          **AIC: −405.81 v −408.61**    AIC: 94.61 v 96.60                             AIC: 2775.40 v 2777.40      AIC: 3376.50 v 3378.50
  Ref: Flexible and fixed                                                                                                                                                                                                           

  **Enrichment score**\             AIC: 704.60 v 706.06 \*       AIC: 640.79 v 642.28       **AIC: 405.16 v 423.59 \***   **AIC: −405.81 v −2832.36**   AIC: 95.22 v 96.96 \*    AIC: 40.63 v 42.63                                AIC: 3376.50 v 3376.10
  (0--25)                                                                                                                                                                                                                           

  **% fruit**\                      AIC: 690.20 v 690.38          AIC: 620.37 v 620.90       **AIC: 414.34 v 407.11**      AIC: −372.03 v −370.67        AIC: 96.30 v 98.07 \*    AIC: 40.62 v 42.35    AIC: 2697.80 v 2696.90      
  (0--100%)                                                                                                                                                                                                                         
  ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

When sex was the outcome variable, the models including contraceptive use/neutered and season as predictors were better fits to the data than the null. Animals given contraception were less likely to be male (Z = −2.91, *p* \< 0.01); and animals weighed in spring are more likely to be male (autumn v: spring: Z = 2.98, *p* \< 0.01; summer: Z = −1.06, *p* = 0.29; winter: Z = 0.80, *p* = 0.42).

For the age models, the model with species as a predictor was a better fit to the data than the null. Ring-tailed lemurs were younger than mongoose lemurs (*t*~241~ = −2.34, *p* = 0.02), as were, to a lesser extent, black-and-white ruffed lemurs (*t*~241~ = −1.73, *p* = 0.08; red ruffed v mongoose lemurs: *t*~241~ = −0.65, *p* = 0.52).

For the contraceptive use/neutered models, the models with enrichment score and percentage fruit as predictors were better fits to the data than the null. Animals in more enriched enclosures are more likely to be given contraception/neutered (Z = 3.39, *p* \< 0.01); and animals fed less fruit were more likely to be given contraception/neutered (Z = −2.71, *p* = 0.01).

For the enclosure area models, models including species, season, enclosure type, climbing structures, and enrichment score as predictors were all better fits than the null. Compared with mongoose lemurs, red ruffed lemurs are housed in smaller enclosures (mongoose lemur v: red ruffed lemurs: t~315~ = −1.86, *p* = 0.06; ring-tailed lemurs: *t*~315~ = 0.64, *p* = 0.52; black-and-white ruffed lemurs: *t*~315~ = 0.55, *p* = 0.58). Animals weighed in the spring (*t*~315~ = −3.67, *p* \< 0.001), summer (*t*~315~ = −3.05, *p* \< 0.01), and winter (*t*~315~ = −7.00, *p* \< 0.001) had smaller enclosures than those weighed in the autumn. Indoor only enclosures were smaller than indoor and outdoor enclosures (indoor v indoor and outdoor: t~67~ = −2.19, *p* = 0.03; outdoor v indoor and outdoor: *t*~67~ = −1.61, *p* = 0.11). Enclosures featuring fixed climbing structures only were smaller (*t*~68~ = −2.21, *p* = 0.03), and higher enrichment scores were also associated with smaller enclosures (*t*~317~ = −1037.32, *p* \< 0.001).

For the enclosure type models, the model including species as a predictor was a better fit than the null. Compared with mongoose lemurs, both black-and-white ruffed and red ruffed lemurs are more often housed outdoors (Z = 3.88, *p* \< 0.001; Z = 4.78, *p* \< 0.001).

For the enrichment scores models, the model with enclosure type as a predictor was a better fit to the data than the null. Outdoor enclosures have higher enrichment scores compared with indoor and outdoor enclosures (outdoor v indoor and outdoor: Z = 2.43, *p* = 0.01; indoor v indoor and outdoor: Z \< 0.01, *p* = 0.99).

Thus, predictors identified here as being related had their possible effects on focal predictors assessed during the univariable stage, by including these correlated terms in the relevant models (see also [Section 2.5.2](#sec2dot5dot2-animals-10-01443){ref-type="sec"}).

animals-10-01443-t0A5_Table A5

###### 

Details of predictor variables found to have a univariable relationship with the epidemiological outcome, relative body mass. Random effects were 'enclosure' nested in 'zoo'. For each predictor varying intercepts, but common slopes, were found to be the best fit to the data. n = number of animals. AIC ^−185.23^ = AIC value of the null model (i.e., without predictor terms). Coef. = coefficient. SE = standard error. df = degrees of freedom.

  Predictor Details      Comparisons to Null            Model Coefficients                                            
  ---------------------- ------------------------------ -------------------- --------- ------- ------- ------ ------- ------
  Species                Black-and-white ruffed lemur   36                   −206.70   −0.17   0.06    169    −2.80   0.01
  Red ruffed lemur       22                             −0.05                0.06      169     −0.78   0.44           
  Ring-tailed lemur      183                            0.03                 0.05      169     0.49    0.62           
  Mongoose lemur (ref)   15                             \-                   \-        \-      \-      \-             
  Sex                    Male                           146                  −189.18   0.05    0.02    171    2.45    0.02
  Female (ref)           110                            \-                   \-        \-      \-      \-             
  Age                    Days                           256                  −187.89   0.04    0.02    171    2.16    0.03

animals-10-01443-t0A6_Table A6

###### 

Details of baseline model featuring three predictors, and one interaction, all previously identified as having a univariable relationship with relative body mass ([Table A5](#animals-10-01443-t0A5){ref-type="table"}). The species × sex interaction was found to be a better fit to the data than a simpler model without it (AIC = −211.93 v −216.09). Random effects were 'enclosure' nested in 'zoo'. For each predictor varying intercepts but common slopes were found to be the best fit to the data. n = number of animals. Coef. = coefficient. SE = standard error. df = degrees of freedom.

  -----------------------------------------------------------------------------------------------------------------
  Predictor Details            Model Coefficients                                                          
  ---------------------------- ------------------------------------ ------- ------- ------ ------- ------- --------
  Species                      Black-and-white ruffed lemur         36      −0.22   0.08   164     −2.91   \<0.01

  Red ruffed lemur             22                                   −0.10   0.08    164    −1.20   0.23    

  Ring-tailed lemur            183                                  −0.05   0.07    164    −0.68   0.50    

  Mongoose lemur (ref)         15                                   \-      \-      \-     \-      \-      

  Sex                          Male                                 146     −0.10   0.07   164     −1.47   0.14

  Female (ref)                 110                                  \-      \-      \-     \-      \-      

  Age                          Days                                 256     0.05    0.02   164     2.60    0.01

  Species × Sex\               Male: Black-and-white ruffed lemur   22      0.14    0.09   164     1.60    0.11
  Sex ref: female                                                                                          

  Male: Red ruffed lemur       13                                   0.12    0.09    164    1.32    0.19    

  Male: Ring-tailed lemur      103                                  0.17    0.07    164    2.39    0.02    

  Male: Mongoose lemur (ref)   8                                    \-      \-      \-     \-      \-      
  -----------------------------------------------------------------------------------------------------------------

![Images of different adult lemurs in one of three body conditions as deemed by their relative body masses: healthy, overweight, or obese. The number at the bottom of each image is that animal's relative body mass at the time the photograph was taken. Sex of lemurs (credits). **Far left**, **top**: female (Lakeland Wildlife Oasis); **middle**: male (Nancy Nill, Palm Beach Zoo); **bottom**: female (Ashley Ashcraft). **Middle**, **top**: male (Valerie Schultz, Smithsonian's National Zoo); **middle**: female (Śląski Ogród Zoologiczny); **bottom**: male (Heidi Beal). **Right**, **top**: male (Nancy Nill, Palm Beach Zoo); **second from top**: female (Śląski Ogród Zoologiczny); **second from bottom**: male (Debbie Fenton); **bottom**: male (Rebecca Lambert, Taronga Conservation Society).](animals-10-01443-g001){#animals-10-01443-f001}

![Box-and-whisker plots of relative body mass values of individual adult lemurs, from the 13 species featuring in our comparative analyses. Ring-tailed lemurs, mongoose lemurs, black-and-white ruffed lemurs, and red ruffed lemurs also featured in our epidemiological analyses. Sample sizes are shown in parenthesis next to species' names, and are also are indicated by the relative width of their respective boxplot. Species' medians are indicated by vertical lines on the boxes; the extent of the boxes indicate their interquartile ranges; whiskers represent values within 1.5 times the interquartile range; and outliers outside this are depicted as points. The healthy range (0.75--1.25) is shown by yellow shading; overweight by orange (1.25--1.5); obese by red (1.5--2); and morbidly obese (\>2) by the darker red.](animals-10-01443-g002){#animals-10-01443-f002}

![Phylogenetic tree of 13 lemur species in our comparative analyses, with their species-typical median relative body mass values shown as tip-points. The colour of a species' tip-point represents its species-typical body condition: the healthy range (species-typical relative body mass 0.75--1.25) is shown by yellow shading; overweight by orange (1.25--1.5); obese by red (1.5--2); and morbidly obese (\>2) by the darker red.](animals-10-01443-g003){#animals-10-01443-f003}

![Relationship between between-year rainfall coefficient of variation and species-typical median relative body mass across 13 lemur species (*t*~11~ = 2.04, R^2^ = 0.27, λ \< 0.001, *p* = 0.07). The shaded area shows the 95% confidence region. *Labels*: AGL: Alaotran gentle lemur; AA: aye aye; BEBL: blue-eyed black lemur; BL: black lemur; BrL: brown lemur; BWRL: black-and-white ruffed lemur; CL: crowned lemur; RCL: red-collared lemur; RCL: Coquerel's sifaka; ML: mongoose lemur; RBL: red-bellied lemur; RRL: red ruffed lemur; RTL: ring-tailed lemur.](animals-10-01443-g004){#animals-10-01443-f004}

![Relationships between relative body mass and species × sex, age, climbing structure provision × sex, and season × sex. Dotted horizontal line indicates the point at which relative body mass = 1 (i.e., the captive animal is the same weight as its species-typical wild mean; 0 on the log-scale); the short dashed line indicates the overweight threshold (≥1.25, or 0.22 on the log-scale); and the long dashed line indicates the obese threshold (≥1.5, or 0.41 on the log-scale). On the three plots with error bars, points indicate the mean relative body mass value, and the whiskers the upper and lower 95% confidence intervals; differences between groups are indicated with asterisks; and females are indicated by black, males by grey. **Top Left**: BWRL: black-and-white ruffed lemur. ML: mongoose lemur. RRL: red ruffed lemur. RTL: ring-tailed lemur. Male mongoose lemurs have smaller relative body masses than females (*t*~6~ = −2.76, *p* = 0.03), whereas male ring-tailed lemurs have larger values than females (*t*~127~ = 3.15, *p* \< 0.01). **Top Right**: Relative body mass increases with age (*t*~157~ = 2.09, *p* = 0.04) across four lemur species (both log-transformed): black-and-white ruffed lemur = filled square; mongoose lemur = filled triangle; red ruffed lemur = filled circle; ring-tailed lemur = unfilled circle. The shaded area shows the 95% confidence region. **Bottom Left**: Males housed in enclosures featuring fixed climbing structures only (versus flexible and fixed structures) have larger relative body mass values (*t*~157~ = 3.04, *p* \< 0.01). **Bottom Right**: Females weighed in the winter had larger relative body masses than those weighed in spring (*t*~178~ = 3.30, *p* \< 0.01) and summer (*t*~25~ = 3.70, *p* \< 0.01).](animals-10-01443-g005){#animals-10-01443-f005}

animals-10-01443-t001_Table 1

###### 

Details of survey data processing. Individual-level variables are those used in the epidemiology analyses (*n* = number of animals with data), and species-level are those used for comparative analyses (*n* = number of species with data). Note that a subset of the total available data, from the four best-sampled species, was used for epidemiological analyses (see [Section 2.2.2](#sec2dot2dot2-animals-10-01443){ref-type="sec"}).

  ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
  Definition                                                                                                                                                                                                    Levels or Type                  Rationale                                                                                                                                                                                                 Individual-Level Variable (*n*)                                Species-Level Variable (*n*)                   Comments Regarding Species-Level Variable
  ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- ------------------------------- --------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- -------------------------------------------------------------- ---------------------------------------------- --------------------------------------------------------------------
  ***Outcome***                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         

  For adults only, the ratio of the most recent mass recorded (grams) to its species-typical wild mean (grams)                                                                                                  Continuous                      See [Section 2.1](#sec2dot1-animals-10-01443){ref-type="sec"}                                                                                                                                             Relative body mass (544)                                       Median relative body mass (13)                 Median across individuals

  ***Husbandry***                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       

  Type of climbing structures within the enclosure(s)\                                                                                                                                                          Fixed only\                     We assumed fixed climbing structures require less physical effort to climb, and thus might contribute to large body masses, than flexible structures (sensu \[[@B13-animals-10-01443]\])                  Climbing structures (544)                                      Proportion with some flexible climbing\        Proportion across enclosures
  *Fixed:* climbing structures that are rigid and fixed into place, e.g., platforms, bolted down logs\                                                                                                          Flexible and fixed ^1^                                                                                                                                                                                                                                                                                   structures ^2^ (13)                            
  *Flexible:* climbing structures that are unstable and flexible, e.g., ropes, branches on a living tree                                                                                                                                                                                                                                                                                                                                                                                                                                                                

  Contraceptive use/neuter status                                                                                                                                                                               Yes\                            Contraceptive usage is associated with weight gain in primates \[[@B19-animals-10-01443],[@B62-animals-10-01443]\] and castration can cause obesity in other taxa \[[@B14-animals-10-01443]\]             Contraceptive use/neuter status (excluding "Unknown"s) (491)   Proportion given contraception/neutered (12)   Proportion calculated across animals of known contraceptive status
                                                                                                                                                                                                                No\                                                                                                                                                                                                                                                                                                                                                     
                                                                                                                                                                                                                Unknown                                                                                                                                                                                                                                                                                                                                                 

  Total area (m^2^) of the enclosure ^3^                                                                                                                                                                        Continuous                      Used to infer quantity of exercise space, assuming that less exercise, and increased risk of weight gain, occurs in smaller enclosures (sensu \[[@B13-animals-10-01443]\])                                Enclosure area (469)                                           Median enclosure area (13)                     Median across enclosures

  Type of enclosure the animal is housed in.\                                                                                                                                                                   *Ordinally ranked:*\            Levels were ranked according to an assumed increasing need for thermoregulation, and therefore energy expenditure (sensu \[[@B61-animals-10-01443]\]) indoor only \< indoor and outdoor \< outdoor only   Enclosure type (544)                                           Proportion of enclosures which are indoors     Proportion calculated across enclosures
  *Note:* Animals with access to indoor and outdoor enclosures for all or part of the year, or any combination of these were scored as 'indoor and outdoor'                                                     Indoor only\                                                                                                                                                                                                                                                                                                                                            
                                                                                                                                                                                                                Indoor and outdoor\                                                                                                                                                                                                                                                                                                                                     
                                                                                                                                                                                                                Outdoor only                                                                                                                                                                                                                                                                                                                                            

  Types of enrichment scored according to their provision: daily (5), two/three times per week (4), weekly (3), monthly (2), less than monthly (1), never (0). Summed across for a maximum total score of 25\   Counts                          We assumed a less enriched environment represents fewer activity opportunities, which could contribute to weight gain (sensu \[[@B13-animals-10-01443]\])                                                 Enrichment score (544)                                         Median enrichment score (13)                   Median across enclosures
  *Enrichment types:* Foraging-related, manipulable, olfactory, audio/visual, and training                                                                                                                                                                                                                                                                                                                                                                                                                                                                              

  Proportion of the diet fed that is fruit                                                                                                                                                                      Continuous                      Commercial, domesticated fruit are more energy-rich than the fruits wild lemurs consume, and if fed in large amounts could contribute to weight gain \[[@B38-animals-10-01443]\]                          Proportion of the diet that is fruit (529)                     Median proportion of fruit in the diet (13)    Median across enclosures

  ***Demographics (epidemiological analyses only)***                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    

  Age at weighing (days)                                                                                                                                                                                        Continuous                      Older animals are at higher risk of weight gain in other primates \[[@B5-animals-10-01443]\]                                                                                                              Age (365)                                                      \-                                             \-

  Season of weighing                                                                                                                                                                                            Spring\                         In the wild, some species display seasonal variation in body masses, e.g., \[[@B77-animals-10-01443],[@B78-animals-10-01443]\]                                                                            Season of weighing (544)                                       \-                                             \-
                                                                                                                                                                                                                Summer\                                                                                                                                                                                                                                                                                                                                                 
                                                                                                                                                                                                                Autumn\                                                                                                                                                                                                                                                                                                                                                 
                                                                                                                                                                                                                Winter                                                                                                                                                                                                                                                                                                                                                  

  Sex                                                                                                                                                                                                           Female\                         Risk of obesity and its related health problems varies between the sexes across primates \[[@B79-animals-10-01443],[@B80-animals-10-01443],[@B81-animals-10-01443]\]                                      Sex (544)                                                      \-                                             \-
                                                                                                                                                                                                                Male                                                                                                                                                                                                                                                                                                                                                    

  Species name                                                                                                                                                                                                  Ring-tailed lemur\              Species differ in their susceptibility to weight gain in captivity \[[@B18-animals-10-01443],[@B19-animals-10-01443]\]                                                                                    Species name (544)                                             \-                                             \-
                                                                                                                                                                                                                Black-and-white ruffed lemur\                                                                                                                                                                                                                                                                                                                           
                                                                                                                                                                                                                Red ruffed lemur\                                                                                                                                                                                                                                                                                                                                       
                                                                                                                                                                                                                Mongoose lemur                                                                                                                                                                                                                                                                                                                                          
  ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

^1^ Only one enclosure featured flexible only climbing structures, so this was pooled into the 'Flexible and fixed' level for analysis. ^2^ For animals with both indoor and outdoor enclosures, if one had flexible climbing structures but not the other they were recorded as 'flexible'. ^3^ Summed across indoor and outdoor enclosures.

animals-10-01443-t002_Table 2

###### 

Table describing predictor variables used to test our comparative hypotheses, their rationale, calculation, and predicted relationship with species-typical median relative body mass. *n* = number of species with data.

  -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
  Hypothesis:\                                          Predictor Variable (*n*)                                                                                                                                                                                                                                                                                                   Rationale                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             Predicted Relationship with Species-Typical Relative Body Mass
  Susceptibility to Captive Weight Gain Relates to...                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
  ----------------------------------------------------- -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- ----------------------------------------------------------------
  \(i\) low productivity                                Annual rainfall (13)                                                                                                                                                                                                                                                                                                       Annual rainfall is commonly used in the primate literature to infer primary productivity \[[@B28-animals-10-01443],[@B29-animals-10-01443],[@B30-animals-10-01443]\], because it directly affects plant growth \[[@B26-animals-10-01443],[@B27-animals-10-01443]\] and thus food availability. We extracted the mean annual rainfall within each species' native geographic range (from \[[@B83-animals-10-01443]\]) for years 1970--2000 (from \[[@B84-animals-10-01443]\]) at 1 km^2^ spatial resolution using QGIS (version: Maderia 3.4; QGIS Development \[[@B85-animals-10-01443]\]), taking a median across range fragments for species with a fragmented range ^a^.                                                                                                                           Negative

  \(ii\) unpredictable food resources                   Between-year annual rainfall coefficient of variation (13)                                                                                                                                                                                                                                                                 We assumed that greater between-year variation in rainfall represents greater unpredictability in primary productivity and, thus, food resource availability (cf. \[[@B26-animals-10-01443],[@B27-animals-10-01443]\]). Using 0.5 × 0.5 gridded values, within each species' geographic range we calculated the mean total rainfall for every month between 1901--2016 (CRU version 4.01 \[[@B86-animals-10-01443]\]) in R (version 3.6.1 \[[@B87-animals-10-01443]\]) using the packages "maptools" \[[@B88-animals-10-01443]\], "raster" \[[@B89-animals-10-01443]\] and "rgdal" \[[@B90-animals-10-01443]\]. By summing a given year's monthly means we calculated annual total rainfall values, then calculated the coefficient of variation of these $CV\  = \ \frac{\sigma}{\mu} \times 100$\   Positive
                                                                                                                                                                                                                                                                                                                                                                                   Where $\sigma$ is the standard deviation of annual total rainfall values, and $\mu$ is the mean of these values.                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      

  \(iii\) natural arboreality                           Ground use, common v rare (11)                                                                                                                                                                                                                                                                                             From descriptive terms in published literature, we categorised species' ground use as "common" (e.g., often, sometimes, regularly forages and/or nests on the ground) or "rare" (e.g., rarely, almost never). Sources were: Mittermeier et al. \[[@B41-animals-10-01443]\], Garbutt \[[@B43-animals-10-01443]\], Macdonald \[[@B91-animals-10-01443]\]; and field studies of native, free-living wild-born lemur populations identified during systematic literature searches in Web of Science and British Library's EThOS Thesis Repository \[[@B92-animals-10-01443]\], using species' scientific and common names as terms. Using a species % scans at ground level (below), we also arbitrarily classed species spending ≥10% scans at ground-level as 'common' ^b^.                             Common \> rare

  Ground time, % scans (8)                              We used data from field studies identified in our systematic literature searches to calculate time spent at ground-level (% scans). For species with values from \>1 study, we calculated the median value across studies.                                                                                                 Positive                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              

  Relative tail length (13)                             The proportion of a species' total length (head to tail tip) that is the tail (from \[[@B41-animals-10-01443],[@B43-animals-10-01443]\]). Larger values correspond to relatively longer tails: an adaptation to arboreality across primates \[[@B93-animals-10-01443],[@B94-animals-10-01443],[@B95-animals-10-01443]\].   Negative                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              

  \(iv\) predation risk                                 Predation score,\                                                                                                                                                                                                                                                                                                          Data on direct predation pressure are rare, but broad taxonomic descriptions of predators (e.g., "raptors") are described \[[@B41-animals-10-01443],[@B43-animals-10-01443]\]. We assumed that being predated by an increasing number of different classes of species (i.e., mammals, birds, reptiles) likely represents increasing evolutionary investment in different predator-avoidance strategies by lemurs. We scored the maximum of 1 each for reported mammalian, avian, and/or reptilian predation on a species, and summed across these for a maximum possible score of 3 (from \[[@B41-animals-10-01443],[@B43-animals-10-01443]\]).                                                                                                                                                       Negative
                                                        0--3 (13)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
  -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

^a^ Note that while we could have also calculated annual rainfall values using data from CRU version 4.01 \[[@B86-animals-10-01443]\], our source here \[[@B84-animals-10-01443]\] provided these values at a finer, and thus preferred, scale. ^b^ Descriptive terms and behavioural data showed good agreement in all but one case. As the study providing behavioural data of crowned lemurs', *Eulemur coronatus*, ground-use (1% of scans) was only 4 months long, we judged the field-guides to be more reliable and categorised its ground-used as 'common'.

animals-10-01443-t003_Table 3

###### 

Details of species-typical variables used during comparative analyses. Acronyms are explained at the foot of the table. Species-typical outcome and husbandry data processing and calculations are described in [Section 2.2.1](#sec2dot2dot1-animals-10-01443){ref-type="sec"} and [Table 1](#animals-10-01443-t001){ref-type="table"}. Justification and calculation of species-typical wild ecology predictors are detailed in [Table 2](#animals-10-01443-t002){ref-type="table"}. Dashes indicate data were not available for that species.

  Species Name                     Common Name                    Species-Typical Outcome   Species-Typical Wild Ecology   Species-Typical Husbandry                                                                                        
  -------------------------------- ------------------------------ ------------------------- ------------------------------ --------------------------- --------- ------- -------- ------ ------ --- ----------- ------ ------ ------ ------ ------
  *Daubentonia madagascariensis*   Aye aye                        1.07                      12                             0.07                        2472.25   16.05   Common   \-     0.59   0   26.13       14     0.13   0.9    0      1
  *Eulemur collaris*               Red-collared lemur             1.11                      17                             0.09                        1445.59   19.59   Rare     2.5    0.57   2   38.05       18     0.29   1      0.15   0.20
  *Eulemur coronatus*              Crowned lemur                  1.28                      21                             0.24                        1438.58   26.0    Common   1      0.56   2   72          15     0.24   1      0      0.17
  *Eulemur flavifrons*             Blue-eyed black lemur          1.41                      13                             0.19                        1706.97   19.86   Rare     0.6    0.61   2   49.10       14     0.29   0.89   0      0
  *Eulemur fulvus*                 Brown lemur                    1.53                      11                             0.16                        1482.72   14.27   \-       \-     0.46   2   12.95       18     0.25   0.86   0.45   0.14
  *Eulemur macaco*                 Black lemur                    1.29                      14                             0.13                        1961.93   23.78   Rare     0.58   0.58   2   72.59       14.5   0.05   1      0.5    0.17
  *Eulemur mongoz*                 Mongoose lemur                 1.16                      29                             0.17                        1481.21   15.60   Rare     1      0.59   2   36          18     0.33   1      0.14   0.24
  *Eulemur rubriventer*            Red-bellied lemur              1.01                      14                             0.15                        1718.98   14.55   Rare     \-     0.56   2   102.25      13.5   0.29   0.83   0      0
  *Hapalemur alaotrensis*          Alaotran gentle lemur          0.97                      12                             0.09                        1151.84   14.20   Common   \-     0.50   2   20          17     0      1      0.08   0
  *Lemur catta*                    Ring-tailed lemur              1.25                      351                            0.24                        701.21    20.03   Common   27     0.58   3   142.91      14     0.27   0.90   0.33   0.08
  *Propithecus coquereli*          Coquerel's sifaka              1.08                      17                             0.14                        1527.47   17.23   Common   \-     0.54   0   14,188.39   14     0.01   1      \-     0
  *Varecia rubra*                  Red ruffed lemur               1.24                      75                             0.22                        2873.58   20.56   Rare     0.01   0.57   0   227.60      14     0.40   0.93   0.19   0.10
  *Varecia variegata*              Black-and white ruffed lemur   1.05                      89                             0.15                        1869.95   16.75   Rare     1.1    0.54   0   144.79      14     0.29   0.97   0.12   0.05

Acronyms: *Outcome variable:* RBM: species-typical median relative body mass; n: number of individuals with outcome data; SD: standard deviation of species-typical median relative body mass. *Species-typical wild ecology predictor variables:* AR: annual rainfall (mm); ARCV: between-year rainfall coefficient of variance; GU: ground use (common v rare); GT: ground time (% scans); RTL: relative tail length (proportion of tail length to whole body length); PS: predation score (0--3). *Species-typical husbandry variables:* MEA: median enclosure area; MEE: median enrichment score; MPF: median proportion of the captive diet that is fruit; PFC: proportion of enclosures with flexible climbing materials; PGC: proportion of captive animals known to be given contraception/neutered. PHI: proportion of enclosures which were indoors.

animals-10-01443-t004_Table 4

###### 

Results of hypothesis-testing model results performed over a tree block of 1000 alternative Lemuriform phylogenetic trees. Because values for upper and lower 95% confidence intervals were identical to medians in all cases, here we provide median values for each parameter. CV = coefficient of variation. Results were considered significant at *p* \< 0.05, and trends (*p* \< 0.10) are italicised. All *p* values are two-tailed.

  Hypothesis                Wild Ecology Predictor                                 Model Output \[95% CIs\]
  ------------------------- ------------------------------------------------------ -------------------------------------------------------
  \(i\) low productivity    Total annual rainfall (mm)                             *t*~11~ = −0.10, R^2^ \< 0.01, λ \< 0.01, *p* = 0.92
  \(ii\) unpredictability   Between-year rainfall CV                               *t*~11~*= 2.04, R^2^ = 0.27*, *λ \< 0.01*, *p = 0.07*
  \(iii\) arboreality       Ground use (some v rare)                               *t*~10~ = −0.59, R^2^ = 0.03, λ \< 0.01, *p* = 0.57
  Ground time (% scans)     *t*~6~ = −0.31, R^2^ = 0.02, λ \< 0.01, *p* = 0.76     
  Relative tail length \*   *t*~11~ = −0.15, R^2^ \< 0.01, λ \< 0.01, *p* = 0.89   
  \(iv\) predation risk     Predation score (0--3)                                 *t*~11~ = 1.18, R^2^ = 0.11, λ \< 0.01, *p* = 0.26

\* This is the proportion of the length of the tail compared to the entire length of the animal, from head to tail tip.

animals-10-01443-t005_Table 5

###### 

Final model of predictors that explain relative body mass. Random effects were 'enclosure' nested in 'zoo'. For each predictor, varying intercepts but common slopes were found to be the best fit to the data. n = number of animals. AIC ^−216.08^ = AIC value of the baseline model (the first four predictors are the baseline model itself, hence they do not have AIC values). Coef. = coefficient. SE = standard error. df = degrees of freedom. 95% CI^LL^ = lower 95% confidence interval of the coefficient; 95% CI^UL^ = upper 95% confidence interval of the coefficient.

  -------------------------------------------------------------------------------------------------------------------------------------------------------------
  Predictor Details                Comparisons to Baseline Model        Model Coefficients                                                             
  -------------------------------- ------------------------------------ -------------------- --------- ------- -------- ------- ------- ------ ------- --------
  Species                          Black-and-white ruffed lemur         36                   \-        −0.22   −0.37    −0.08   0.08    157    −2.94   \<0.01

  Red ruffed lemur                 22                                   −0.10                −0.25     0.06    0.08     157     −1.20   0.23           

  Ring-tailed lemur                183                                  −0.06                −0.19     0.06    0.07     157     −0.92   0.36           

  Mongoose lemur (ref)             15                                   \-                   \-        \-      \-       \-      \-      \-             

  Sex                              Male                                 146                  \-        −0.17   −0.32    −0.01   0.08    157    −2.06   0.04

  Female (ref)                     110                                  \-                   \-        \-      \-       \-      \-      \-             

  Age                              Days                                 256                  \-        0.04    \<0.01   0.07    0.02    157    2.09    0.04

  Species × Sex\                   Male: Black-and-white ruffed lemur   22                   \-        0.16    0.00     0.32    0.08    157    1.86    0.07
  Sex ref: female                                                                                                                                      

  Male: Red ruffed lemur           13                                   \-                   0.13      −0.04   0.30     0.09    157     1.43   0.15    

  Male: Ring-tailed lemur          103                                                       0.21      0.07    0.35     0.07    157     2.78   0.01    

  Male: Mongoose lemur (ref)       8                                                         \-        \-      \-       \-      \-      \-     \-      

  Climbing structures              Fixed                                19                   −214.66   −0.93   −0.24    0.06    0.08    27     −1.15   0.26

  Flexible and fixed (ref)         239                                                       \-        \-      \-       \-      \-      \-     \-      

  Season                           Spring                               37                   −220.20   −0.11   −0.24    0.01    0.06    157    −1.76   0.08

  Summer                           108                                  −0.04                −0.12     0.05    0.04     157     −0.83   0.41           

  Winter                           37                                   0.13                 0.02      0.25    0.06     157     2.18    0.03           

  Autumn (ref)                     74                                   \-                   \-        \-      \-       \-      \-      \-             

  Sex × Climbing structures\       Male: Fixed                          9                    −221.01   0.24    0.09     0.39    0.08    157    3.04    \<0.01
  Sex ref: female                                                                                                                                      

  Male: Flexible and fixed (ref)   137                                  \-                   \-        \-      \-       \-      \-      \-             

  Sex × Season\                    Spring                               26                   −226.70   0.09    −0.05    0.22    0.07    157    1.22    0.23
  Sex ref: female                                                                                                                                      

  Summer                           56                                   0.07                 −0.02     0.17    0.05     157     1.46    0.15           

  Winter                           18                                   −0.15                −0.29     −0.01   0.07     157     −2.06   0.04           

  Autumn (ref)                     46                                   \-                   \-        \-      \-       \-      \-      \-             
  -------------------------------------------------------------------------------------------------------------------------------------------------------------
