Transition to siblinghood causes a substantial and long-lasting increase in urinary cortisol levels in wild bonobos

  1. Verena Behringer  Is a corresponding author
  2. Andreas Berghänel
  3. Tobias Deschner
  4. Sean M Lee
  5. Barbara Fruth
  6. Gottfried Hohmann
  1. Endocrinology Laboratory, German Primate Center, Leibniz Institute for Primate Research, Germany
  2. Max Planck Institute for Evolutionary Anthropology, Germany
  3. Domestication Lab, Konrad Lorenz Institute of Ethology, Department of Interdisciplinary Life Sciences, University of Veterinary Medicine Vienna, Austria
  4. Comparative BioCognition, Institute of Cognitive Science, University of Osnabrück, Germany
  5. Center for the Advanced Study of Human Paleobiology, Department of Anthropology, George Washington University, United States
  6. Max Planck Institute of Animal Behavior, Germany
  7. Centre for Research and Conservation, Royal Zoological Society of Antwerp, Belgium
2 figures, 3 tables and 3 additional files

Figures

Figure 1 with 5 supplements
Physiological changes in cortisol (A–C), neopterin (D–F), and total T3 (triiodothyronine) (G–I) levels in the older offspring, 7 years before and after sibling birth (sibling birth at 0).

Data points are physiological measures corrected for specific gravity (SG). All smooths are not controlled for age to show cumulative pattern. Axes for physiological variables are log-transformed. 95% confidence intervals are plotted. Left-hand plots (A, D, G): sex-specific trajectories around sibling birth (blue: males; red: females). Middle plots (B, E, H): age-specific trajectories around sibling birth for offspring that were older (purple) or younger (yellow) than the median value of 5.1 years at sibling birth. Right-hand plots (C, F, I): interaction plots visualizing how trajectories around sibling birth change with increasing offspring age at sibling birth (scale from dark green [lowest levels] to brown [highest levels]; white space: extrapolation would be unreliable due to lacking data) for the respective perspective plots, see Figure 1—figure supplement 5. (A) Urinary cortisol levels showed a significant, sudden rise to fivefold values at sibling birth (dotted line); no sex differences or age effects. (B, C) The sudden rise in cortisol levels was independent of the age of the older offspring at sibling birth. (D) Urinary neopterin levels decreased by 1/3 at sibling birth (dotted line; no sex differences or age effects). (E, F) The sudden decrease in neopterin levels was independent of the age of the older offspring at sibling birth. (G–I). Urinary total T3 levels increased around sibling birth, but this effect was indistinguishable from a general age effect. There was no significant sudden change at sibling birth in total T3 levels (G), and there was no significant effect of the age at sibling birth (H, I).

Figure 1—figure supplement 1
Explanation of the principal concept of the applied series of statistical models using the example of our main analysis on offspring urinary cortisol and neopterin levels 7 years before and after sibling birth (0 marks the time of sibling birth).

Confidence intervals: 95% (dark gray) and 99.9% (light gray). Dotted lines: sibling birth (left) and potential second discontinuity in data at 7 (cortisol) and 4.5 months (neopterin) after sibling birth. Data points are original values of physiological measures corrected for specific gravity (SG) and are depicted with females in red and males in blue. All smooths are not controlled for age to show cumulative pattern. Axes for physiological variables are log-transformed. Real models were additionally allowed for sex-specific trajectories (see main text). (A–C) Response variable: urinary cortisol. (A) For all our physiological and behavioral response variables, we first ran a model with a continuous smooth only that allows for nonlinear modeling of changes around sibling birth but not for a discontinuity ( = sudden change) in values. (B) Next, we ran an identical model as in (A) but additionally allowing for an intercept difference between before and after the birth of a younger sibling, thereby allowing for a discontinuity and thus an abrupt change of response values at sibling birth. Significance of this discontinuity was estimated through model comparison between (A) and (B). (C) The same as (B) but additionally with separate smooths ( = trajectories) for offspring that were younger or older than the median age (5.11 years) at sibling birth, allowing for a change in pattern with increasing offspring age. Significance of this age difference was estimated through model comparison between (B) and (C). To allow for visual inspection, we further provide continuous interaction plots showing how trajectories in the response variables around sibling birth change continuously with the age of the offspring at sibling birth. (D–F) Same as (A–C) but for urinary neopterin.

Figure 1—figure supplement 2
Different models of continuous smooths of (A) cortisol and (B) neopterin levels around sibling birth that are allowed for high levels of wiggliness and thus overfitting.

We tested whether the additional allowance for one or two discontinuities still provides a better model fit in a model comparison if compared with such continuous but highly flexible smooths that could theoretically also sufficiently fit the sudden changes. k = number of basis functions (here set to 50 for the predictor variables but kept at 6 for the random smooths to allow for high wiggliness and also for model comparison with the discontinuous models). sp = smoothing penalty, set to low values (and deactivating the default, automatic smoothing penalty estimation of the generalized additive mixed models [GAMM]).

Figure 1—figure supplement 3
Scatter plot of the older sibling urinary cortisol data (blue: males; red: females) in relation to sibling birth (A) with a vertical dotted line at sibling birth (sibling birth is at 0) and (B) with two vertical dotted lines, one at sibling birth and the second one at the end of a 7-month period.

Between the two vertical dotted lines, we did not find any low cortisol levels. Data points are physiological measures corrected for specific gravity (SG).

Figure 1—figure supplement 4
Physiological changes in cortisol (A–C) and neopterin (D–F) levels in the older offspring 7 years before and after sibling birth (sibling birth is at 0) with a sudden change at sibling birth and a second sudden change after a 7-month period (cortisol) or 4.5-month period (neopterin).

Data points are physiological measures corrected for specific gravity (SG). All smooths are not controlled for age to show cumulative pattern. Axes for physiological variables are log-transformed. 95% confidence intervals are plotted. Left-hand plots (A, D): sex-specific trajectories around sibling birth (blue: males; red: females). Middle plots (B, E): age-specific trajectories around sibling birth, for offspring that were older (purple) or younger (yellow) than the median value of 5.1 years at sibling birth. Right-hand plots (C, F): interaction plots visualizing how trajectories around sibling birth change with increasing offspring age at sibling birth (scale from dark green [lowest levels] to brown [highest levels]; white space: extrapolation would be unreliable due to lacking data).

Figure 1—figure supplement 5
Perspective plots of cortisol, neopterin and total T3, showing how the trajectories of the physiological measures change with increasing age of the offspring at sibling birth.

The plots are identical to the plots in Figure 1 and represent another type of visualization of the contour plots shown in Figure 1C, F, and I. Z-axes represent levels of physiological measures (log-transformed and corrected for specific gravity). All smooths are not controlled for age to show cumulative pattern. (A) Cortisol levels before (left) and after (right; both on same scale) sibling birth. (B) Neopterin levels before (left) and after (right; both on same scale) sibling birth. (C) Total T3 levels before and after sibling birth. Red line = time of sibling birth.

Figure 2 with 2 supplements
Behavioral changes in nipple contact (A–C), riding (D–F), body contact (G–I) and 5 m proximity (J–L) with the mother of the older sibling in relation to sibling birth (sibling birth is set to 0).

Vertical dotted lines = time of putative conception (left dotted line) and sibling birth (right dotted line). Data points represent the proportion of time and circle size the underlying sample size (square-rooted; ranges: riding 3–44, all other behaviors 3–303). All smooths are not controlled for age to show cumulative pattern. 95% confidence intervals are plotted. Left-hand plots (A, D, G, J): sex-specific trajectories around sibling birth (blue: males; red: females). Middle plots (B, E, H, K): age-specific trajectories around sibling birth for offspring that were older (purple) or younger (yellow) than the median value of 5.1 years at sibling birth. Right-hand plots (C, F, I, L): interaction plots visualizing how trajectories around sibling birth change with increasing offspring age at sibling birth (scale from dark green [lowest levels] to brown [highest levels]; white space: extrapolation would be unreliable due to lacking data) for the respective perspective plots, see Figure 2—figure supplement 1. (A–C) Proportion of time spent suckling decreased to zero already before sibling birth (A) and was largely absent in older offspring (B, C), without a sudden change at sibling birth. (D–F) The proportion of time riding on the mother showed a significant and sudden decline at sibling birth (D), but this cut was evident only in offspring younger than 5 years old at sibling birth and not anymore in older offspring (E, F). (G–I) The proportion of time spent in body contact with the mother showed a significant and sudden increase at sibling birth, irrespective of the sex or age of the offspring. (J–L) The proportion of time in 5 m proximity to the mother decreased around sibling birth, but this effect was indiscernible from a general age effect. There was no significant sudden change at sibling birth (J), and there was no significant effect of offspring age at sibling birth (K, L).

Figure 2—figure supplement 1
Perspective plots of niplle contact, riding on the mother, body contact and proximity to the mother, showing how the trajectories of the behavioral measures change with increasing age of the offspring at sibling birth.

The plots are identical with the plots in Figure 2 and represent another type of visualization of the contour plots (Figure 2C, F, I, and L). Z-axes represent proportion of time spent (A) nipple contact, (B) riding on mother, (C) in body contact with the mother, and (D) in 5 m proximity to the mother. All smooths are not controlled for age to show cumulative pattern. Red line = time of sibling birth.

Figure 2—figure supplement 2
Behavioral changes in the proportion of time spent foraging independently while the mother is foraging (to control for foraging opportunity).

All smooths are not controlled for age to show cumulative pattern. 95% confidence intervals are plotted. (A, B) Vertical dotted lines = time of putative conception (left dotted line) and sibling birth (right dotted line). Data points represent the proportion of time and circle size the underlying sample size (square-rooted; range 1–182). (A) Sex-specific trajectories around sibling birth (blue: males; red: females). (B) Age-specific trajectories around sibling birth for offspring that were older (purple) or younger (yellow) than the median value of 5.1 years at sibling birth. (C, D) Interaction plots visualizing how trajectories around sibling birth change with increasing offspring age at sibling birth. (C) Contour plot: scale from dark green (lowest levels) to brown (highest levels). White space: extrapolation would be unreliable due to lacking data. Vertical dotted lines = time of putative conception and sibling birth. (D) The respective perspective plot. The z-axis represents the proportion of time spent foraging independently. Red line = time of sibling birth.

Tables

Table 1
Summary of the main findings of analyses of physiological markers and scores of older offspring behavioral during the transition to siblinghood (TTS).
CortisolNeopterinTotal T3NursingRiding5m-proximity with motherBody contact with motherIndependent foraging
Sudden change at sibling birthYes, increaseYes, decreaseNoNoYesNoYes, but increaseNo
Effect of TTS attenuates with offspring ageNoNoNoNo: all changes before sibling birthYes, effect exists only up to 5 years oldNoNoNo: all changes finished before sibling birth
Table 2
General additive mixed model results for physiological changes (urinary cortisol, urinary neopterin, and urinary total T3 levels; all log-transformed) in the older offspring 7 years before and after sibling birth.
ReferenceLog cortisolLog neopterinLog total T3
Factor variables:CategoryEst.SEtpEst.SEtpEst.SEtp
(intercept)0.850.0516.232.410.0458.680.890.0517.26
MalesFemales0.110.052.090.037–0.030.04–0.690.488–0.110.05–2.190.030
After S-birth*Before0.430.085.27<0.001–0.190.06–3.010.0020.130.081.610.114
Smooth term variables:edfRef. dfFpedfRef. dfFpedfRef. dfFp
Time-S-birth: males2.653.171.170.2621.001.000.530.4691.001.003.760.054
Time-S-birth: females1.772.120.690.4331.501.830.560.6031.001.000.010.922
Time-S-birth: malesFemales1.001.000.050.8181.001.000.230.5851.001.002.850.093
Age: males1.001.003.220.0743.193.744.250.002----
Age: females1.001.001.370.2431.001.000.030.874----
Age: malesFemales1.001.000.430.5132.172.631.080.243----
Daytime1.201.3729.27<0.0011.001.004.820.0292.102.561.770.142
Seasonal effect2.383.0011.18<0.0010.513.000.220.2780.003.000.000.723
Random effects:
Time-S-birth per ID (smooth)0.00111.000.2380.00112.000.8507148.000.015
Age per ID (smooth)0.00109.000.2910.00108.000.737----
Mother ID (intercept)0.0013.000.1750.0013.000.313013.000.230
Year (intercept)0.001.000.0120.001.000.27701.000.850
R2adj (deviance explained)0.311 (33.8%)0.169 (19.6%)0.117 (15.3%)
N (p-value, full/null comp)319 (<0.001)314 (<0.001)319 (0.020)
  1. Green indicates classic interaction term derived from a separate model calculation (see ‘Methods’). Data points are physiological measures corrected for specific gravity (SG). All smooths are not controlled for age to show cumulative pattern.

  2. ID = individual; T3 = total triiodothyronine; S-birth = sibling birth. .

  3. *

    After sibling birth.

  4. Before sibling birth.

Table 3
Generalized additive mixed model (GAMM) results of behavioral changes (nipple contact, riding, and body contact and 5 m proximity with the mother) in the older offspring around sibling birth (±2 years).

Binomial GAMMs on proportions of time per day and individual.

ReferenceNipple contactRidingProximityBody contact with mother
Factor variables:CategoryEst.SEzpEst.SEzpEst.SEzpEst.SEzp
(intercept)–7.340.84–8.73–1.360.52–2.650.210.161.30–2.790.19–14.50
MalesFemales0.760.611.240.211.180.363.220.001–0.080.13–0.660.510.220.112.050.040
After YS-birth*Before1.390.881.550.12–2.000.55–3.64<0.0010.020.100.180.850.470.114.17<0.001
Year----------------
Smooth term Variables:edfRef. dfChi²pedfRef. dfChi²pedfRef. dfChi²pedfRef. dfChi²p
T-S-birth: males4.454.9512.380.0021.001.000.030.861.001.006.910.0091.001.0012.48<0.001
T-S-birth: females3.724.1912.380.0171.001.005.240.0223.283.378.030.0321.001.0028.70<0.001
T-S-birth: malesFemales1.001.000.200.661.001.002.780.0951.001.000.050.831.001.001.040.31
Age: males1.001.000.550.461.001.0019.89<0.001--------
Age: females1.001.000.120.711.001.004.390.036--------
Age: malesFemales1.001.000.000.991.001.004.380.036--------
Daytime3.643.9214.310.0123.513.857.460.093.813.98170.67<0.0013.964.00456.4<0.001
Seasonal effect1.103.001.890.0218.893.001.910.0630.003.000.000.052.643.0061.22<0.001
Random effects:
Time-S-birth per ID (smooth)3.13113.0082.07<0.00133.3270.00251.9<0.00162.1176.001323.3<0.0015876.01569<0.001
Age per ID (smooth)8.0594.0034.07<0.0010.0061.000.000.010--------
Mother ID (intercept)2.4510.000.00<0.0010.0010.000.000.0016.2812.000.00<0.0010.0112.00.01<0.001
Date (intercept)3.651.000.000.230.001.000.000.280.001.000.00<0.0010.001.000.00<0.001
R2adj (deviance explained)0.39 (62.4%)0.827 (81.7%)0.226 (29.3%)0.319 (39.7%)
N (p-value, full/null comp)545 (<0.001)301 (<0.001)545 (<0.001)545 (<0.001)
  1. Green: classic interaction term derived from a separate model calculation (see ‘Methods’). Statistics for year (categorical control variable) not shown for clarity.

  2. ID: individual; S-birth = sibling birth; ‘:’ = interaction term.

  3. *

    After sibling birth.

  4. Before sibling birth.

Additional files

Supplementary file 1

General additive mixed model results for physiological changes (urinary cortisol and urinary neopterin, levels; all log-transformed) in the older offspring 7 years before and after sibling birth.

Green: classic interaction term derived from a separate model calculation (see ‘Methods’). ID = individual; S-birth = sibling birth; * before = before sibling birth; * early after = 7 and 4.5 months following sibling birth for cortisol and neopterin, respectively; * late after = time following early after. Data points are physiological measures corrected for specific gravity (SG). All smooths are not controlled for age to show cumulative pattern.

https://cdn.elifesciences.org/articles/77227/elife-77227-supp1-v2.docx
Supplementary file 2

Number of individuals and samples/data points (in brackets) for each physiological marker or behavior shown for each sex in relation to sibling birth.

sib birth = sibling birth; total T3 = total triiodothyronine.

https://cdn.elifesciences.org/articles/77227/elife-77227-supp2-v2.docx
Transparent reporting form
https://cdn.elifesciences.org/articles/77227/elife-77227-transrepform1-v2.pdf

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Verena Behringer
  2. Andreas Berghänel
  3. Tobias Deschner
  4. Sean M Lee
  5. Barbara Fruth
  6. Gottfried Hohmann
(2022)
Transition to siblinghood causes a substantial and long-lasting increase in urinary cortisol levels in wild bonobos
eLife 11:e77227.
https://doi.org/10.7554/eLife.77227