Early maternal loss leads to short- but not long-term effects on diurnal cortisol slopes in wild chimpanzees

  1. Cédric Girard-Buttoz  Is a corresponding author
  2. Patrick J Tkaczynski
  3. Liran Samuni
  4. Pawel Fedurek
  5. Cristina Gomes
  6. Therese Löhrich
  7. Virgile Manin
  8. Anna Preis
  9. Prince F Valé
  10. Tobias Deschner
  11. Roman M Wittig
  12. Catherine Crockford
  1. Department of Human Behavior, Ecology and Culture, Max Planck Institute for Evolutionary Anthropology, Germany
  2. Taï Chimpanzee Project, Centre Suisse de Recherches Scientifiques, Côte d'Ivoire
  3. Department of Primatology, Max Planck Institute for Evolutionary Anthropology, Germany
  4. Department of Human Evolutionary Biology, Harvard University, United States
  5. Division of Psychology, University of Stirling, United Kingdom
  6. Tropical Conservation Institute, Florida International University, United States
  7. World Wide Fund for Nature, Dzanga Sangha Protected Areas, Central African Republic
  8. Robert Koch Institute, Epidemiology of Highly Pathogenic Microorganisms, Germany
  9. Centre Suisse de Recherches Scientifiques en Côte d'Ivoire, Côte d'Ivoire
  10. Unité de Formation et de Recherche Biosciences, Université Félix Houphouët Boigny, Côte d'Ivoire
  11. Interim Group Primatology, Max Planck Institute for Evolutionary Anthropology, Germany
  12. Institut des Sciences Cognitives, CNRS, France
8 figures, 5 tables and 1 additional file

Figures

Effect of maternal loss on daily urinary cortisol level variations in immature chimpanzees.

Non-orphans are depicted by orange squares, recently orphans (orphaned for less than 2 years) by green triangles, `and non-recently orphaned (orphaned for more than 2 years) by blue circles. Each dot represents the average hourly cortisol level for all individuals of each orphan category. The size of the dot is proportional to the sample size (e.g., the number of data points) contributing to each dot. The orange solid line and green and blue dashed lines depict the model line prediction, and the orange, green, and blue area the 95% credible interval (CI) from the all immature model for non-orphans, recently orphans, and non-recently orphans, respectively. The model lines depict the consistent effect of the interaction between orphan status and the linear term for time of day in the all immature model (estimate: –0.22, 95% CI: [–0.03: –0.48]). The sample size for the all immature model was N = 846 samples and 50 individuals.

Effect of time (in years) since being orphaned on daily urinary cortisol level variations of immature orphan chimpanzees.

Orange squares, blue circles, and green triangles depict cortisol levels for individuals who lost their mother less than 3 years ago, 3–6 years ago, and 6–9 years ago, respectively. Each dot represents the average hourly cortisol level of all individuals for each of the three categorical intervals of years since maternal loss. The size of the dot is proportional to the sample size (number of urine samples collected) for each hour of the day. The orange solid line, blue dashed line, and green dashed lines indicate model line predictions, and the orange, blue, and green light areas indicate the 95% credible interval (CI) for individuals who lost their mother less than 3 years ago, 3–6 years ago, and 6–9 years ago, respectively (immature orphan model). The model line depicts the interaction between years since maternal loss and the cubic term for time of day in the immature orphan model: (estimate and [90% CI]: –0.06 [-0.14: 0.00]), but this effect comprised a relatively large uncertainty since the 95% CI overlapped 0 (Table 2). Note that while the predictor ‘years since maternal loss’ was modeled as a continuous variable in the immature orphan model, for ease of interpretation we depict the results here for three categorical intervals of ‘years since maternal loss’. The sample size for the immature orphan model was N = 393 samples and 17 individuals.

Effect of age at which immature orphans lost their mother on daily urinary cortisol level variations of immature orphan chimpanzees.

Orange square, blue circles, and green triangles depict cortisol levels for individuals who lost their mother when they were less than 5 years, 5–8 years, and 8–12 years of age, respectively. Each dot represents the average hourly cortisol level of all individuals for each of the three categories of age when mother died. The size of the dot is proportional to the sample size (number of urine sample collected) for each hour of the day. The orange solid line and blue and green dashed lines indicate model line predictions, and the orange, blue, and green light areas the 95% credible interval (CI) for individuals who lost their mother when they were less than 5 years, 5–8 years, and 8–12 years of age, respectively (immature orphan model). The model lines depict the consistent effect of the interaction between age when mother died and the quadratic term of time of day in the immature orphan model (estimate and [95% CI]: –0.12 [0.00: –0.23]). Note that while ‘age at maternal loss’ was modeled as a continuous variable in the immature orphan model, for ease of interpretation the model is depicted here for three categorical intervals of ’age when mother died’. The sample size for the immature orphan model was N = 393 samples and 17 individuals.

Appendix 1—figure 1
Diurnal cortisol level variation in immature chimpanzees.

Non-orphans are depicted by orange squares, recently orphans (orphaned for less than 2 years) by green triangles, and non-recently orphaned (orphaned for more than 2 years) in blue circles. Each dot represents one sample. As for Figure 1, the orange solid line and green and blue dashed lines depict the model prediction lines from the all immature model Model 1a for non-orphans, recently orphans, and non-recently orphans, respectively. The model lines depict the consistent effect of the interaction between orphan status and the linear term for time of day in the all immature model Model 1a (estimate: –0.22, 95% credible interval [CI]: [–0.03: –0.48]). The sample size for the all immature model Model 1a was N = 846 samples and 50 individuals.

Appendix 1—figure 2
Posterior predictive check for the all immature model.
Appendix 1—figure 3
Posterior predictive check for the immature orphan model.
Appendix 1—figure 4
Posterior predictive check for the all adult male model.
Appendix 1—figure 5
Posterior predictive check for the adult male orphan model.

Tables

Table 1
Sample size for immature and adult male orphans and non-orphans in each of the four study communities.
CommunityAge classOrphan statusN.IDNo. of samplesMean ± SE no. of sample per individualAge range (years)
Taï EastImmaturesNon-orphans8*11214.0 ± 3.83.8–11.9
Recently orphaned5*11623.2 ± 9.14.1–11.9
Non-recently orphaned9*13615.1 ± 6.46.1–11.8
Adult malesNon-orphans4456114.0 ± 56.213.8–40.7
Orphans3354118.0 ± 18.012.3–19.8
Taï MiddleImmaturesNon-orphans
Recently orphaned
Non-recently orphaned
Adult malesNon-orphans3175.7 ± 0.6717.2–33.7
Orphans
Taï NorthImmaturesNon-orphans1116815.3 ± 4.33.0–12.0
Recently orphaned
Non-recently orphaned273.5 ± 0.510.5–11.3
Adult malesNon-orphans39531.7 ± 9.112.3–20.8
Orphans26130.5 ± 3.512.1–20.4
Taï SouthImmaturesNon-orphans17*17310.2 ± 1.52.8–11.9
Recently orphaned2*189.0 ± 6.04.1–9.2
Non-recently orphaned5*11623.2 ± 7.76.1–11.9
Adult malesNon-orphans7858122.6 ± 18.112.1–45.3
Orphans634357.2 ± 31.512.1–21.9
  1. *

    50 immature individuals were included in this study but one immature in Taï East and two immatures in Taï South were sampled before and after their mother died (i.e., they are counted twice in the table, once as an orphan and once as a non-orphan). Six males were included in the study as both mature and immature individuals.

Table 2
Results of the all immature and the immature orphan models.

These two models tested in immatures for the effect of orphan status (all immature model) and, for immature orphan only, the effect of age when mother died and years since maternal loss (orphan immature model) on cortisol levels and diurnal cortisol slopes. The results presented here are from reduced all immature and immature orphan models after removing the interactions for which the 90% credible interval (CI) overlapped 0. SE indicates the standard error of the estimate for each predictor. The coded level for each categorical predictor is indicated in parentheses. Control predictors are italicized. 95% CI low and 95% CI high indicate the lower and upper limits of the 95% CI. Likewise, 90% CI low and 90% CI high indicate the lower and upper limits of the 90% CI. CIs that do not overlap 0 are indicated in bold. Predictors for which the 95% CI did not overlap 0 are indicated with salmon shading, and predictors for which the 95% CI overlapped 0 but the 90% CI did not overlap 0 are indicated with yellow shading. LCMS: liquid chromatography mass spectrometry. Time of the day2: quadratic term for time of the day. Time of the day3: cubic term for time of the day.

ModelResponsePredictorEstimateSE95% CI low95% CI high90% CI low90% CI high
All immature modelLog urinary cortisol levels (ng/ml SG)Intercept3.690.273.174.233.264.14
Time of the day–0.380.080.530.220.510.25
Time of the day2–0.010.08–0.170.14–0.140.11
Time of the day30.030.04–0.050.12–0.040.10
Orphan category (orphan for less than 2 years)0.010.19–0.360.38–0.300.32
Orphan category (orphan for more than 2 years)–0.140.22–0.580.32–0.500.23
Sex ratio0.070.07–0.060.2–0.040.18
Community size–0.120.06–0.240.010.220.01
Individual sex (male)–0.260.25–0.780.22–0.670.13
Individual age at sample–0.010.17–0.330.37–0.270.29
LCMS method (old)0.060.35–0.620.78–0.500.64
Sin(seasonDate)–0.070.04–0.140–0.13–0.02
Cos(seasonDate)0.010.04–0.060.08–0.050.07
Orphan category (less than 2 years): time of the day–0.220.100.410.030.380.07
Orphan category (more than 2 years): time of the day–0.110.08–0.270.03–0.240.01
Immature orphan modelLog urinary cortisol levels (ng/ml SG)Intercept3.590.263.084.113.174.02
Time of the day–0.440.170.750.050.70.14
Time of the day2–0.010.07–0.160.13–0.130.1
Time of the day30.030.11–0.180.27–0.140.22
Orphan’s age when mother died0.000.24–0.510.46–0.410.38
Years since maternal loss–0.160.16–0.480.14–0.420.09
Sex ratio–0.060.11–0.270.15–0.230.11
Community size–0.050.08–0.220.12–0.190.09
Individual sex (male)0.10.3–0.470.7–0.380.6
LCMS method (old)–0.040.32–0.660.64–0.550.51
Sin(seasonDate)0.020.05–0.090.12–0.070.1
Cos(seasonDate)0.040.05–0.070.14–0.050.12
Years since maternal loss: time of the day0.190.1–0.010.380.020.35
Years since maternal loss: time of the day2–0.050.05–0.150.05–0.140.03
Years since maternal loss: time of the day3–0.060.04–0.150.020.140.00
Age when mother died: time of the day0.080.09–0.090.25–0.060.22
Age when mother died: time of the day2–0.120.060.230.000.210.02
Appendix 1—table 1
Results of the all adult male model testing the effect of maternal loss on cortisol profiles in all mature males.

’All adult male model full’ refers to the full model ran with all the interactions considered. ‘All adult male model reduced’ refers to the reduced model after removing the interactions for which the 90% credible interval (CI) overlapped 0. SE indicates the standard error of the estimate for each predictor. The coded level for each categorical predictor is indicated in parentheses. Control predictors are italicized. 95% CI low and 95% CI high indicate the lower and upper limits of the 95% CI. Likewise, 90% CI low and 90% CI high indicate the lower and upper limits of the 90% CI. CIs that do not overlap 0 are indicated in bold. LCMS: liquid chromatography mass spectrometry. Time of the day2: quadratic term for time of the day. Time of the day3: cubic term for time of the day.

ModelResponsePredictorEstimateSE95% CI low95% CI high90% CI low90% CI high
All adult male model fullLog urinary cortisol levels (ng/ml SG)Intercept3.970.253.444.443.554.36
Time of the day–0.480.09–0.67–0.32–0.63–0.35
Time of the day2–0.070.05–0.170.03–0.150.01
Time of the day30.030.04–0.050.12–0.040.1
Orphan status (yes, orphan)0.020.19–0.360.38–0.290.32
Individual age at sample0.190.14–0.070.49–0.030.42
LCMS method (old)0.170.31–0.410.82–0.310.69
Community size0.060.09–0.130.24–0.10.21
Sex ratio0.040.07–0.110.18–0.090.16
Dominance rank0.040.08–0.120.19–0.090.17
Sin(seasonDate)–0.040.02–0.080–0.070
Cos(seasonDate)–0.010.02–0.050.03–0.040.03
Orphan status (yes): time of the day0.020.04–0.050.1–0.040.09
Orphan status (yes): time of the day20.010.04–0.080.09–0.060.08
Orphan status (yes): time of the day3–0.030.08–0.180.13–0.150.1
All adult male model reducedLog urinary cortisol levels (ng/ml SG)Intercept3.970.243.454.433.554.35
Time of the day–0.490.08–0.67–0.33–0.64–0.36
Time of the day2–0.060.04–0.150.03–0.140.01
Time of the day30.040.04–0.030.12–0.020.11
Orphan status (yes, orphan)0.030.18–0.330.38–0.280.32
Individual age at sample0.190.14–0.070.48–0.020.42
LCMS method (old)0.170.3–0.410.8–0.30.68
Community size0.050.09–0.130.24–0.10.21
Sex ratio0.030.07–0.110.17–0.090.15
Dominance rank0.040.08–0.120.19–0.090.17
Sin(seasonDate)–0.040.02–0.080–0.070
Cos(seasonDate)–0.010.02–0.050.03–0.040.03
Appendix 1—table 2
Results of the adult male orphan model testing, in orphan mature males only, the effect of age at maternal loss on cortisol profiles.

’Adult male orphan model full’ refers to the full model ran with all the interactions considered. ‘Adult male orphan model reduced’ refers to the reduced model after removing the interactions for which the 90% credible interval (CI) overlapped 0.SE indicates the standard error of the estimate for each predictor. The coded level for each categorical predictor is indicated in parentheses. Control predictors are italicized. 95% CI low and 95% CI high indicate the lower and upper limits of the 95% CI. Likewise, 90% CI low and 90% CI high indicate the lower and upper limits of the 90% CI. CIs that do not overlap 0 are indicated in bold. LCMS: liquid chromatography mass spectrometry. Time of the day2: quadratic term for time of the day. Time of the day3: cubic term for time of the day.

ModelResponsePredictorEstimateSE95% CI low95% CI high90% CI low90% CI high
Adult male orphan fullLog urinary cortisol levels (ng/ml SG)Intercept4.070.333.414.713.534.58
Time of the day–0.50.13–0.77–0.24–0.72–0.29
Time of the day2–0.080.07–0.210.05–0.190.03
Time of the day30.040.06–0.090.17–0.060.14
Orphan’s age when mother died–0.120.25–0.620.37–0.530.29
Community size0.040.21–0.340.48–0.270.4
Sex ratio0.020.17–0.310.34–0.250.29
Dominance rank–0.070.18–0.410.3–0.360.24
Orphan’s age at sample0.430.37–0.231.24–0.121.08
LCMS method (old)0.030.31–0.550.66–0.450.54
Sin(seasonDate)–0.010.03–0.080.05–0.070.04
Cos(seasonDate)–0.010.03–0.080.06–0.070.05
Orphan’s age when mother died: time of the day0.030.07–0.110.18–0.090.04
Orphan’s age when mother died: time of the day20.060.04–0.020.14–0.010.12
Orphan’s age when mother died: time of the day3–0.020.04–0.110.06–0.090.15
Adult male orphan reducedLog urinary cortisol levels (ng/ml SG)Time of the day4.070.333.434.723.554.6
Time of the day2–0.510.13–0.77–0.26–0.72–0.3
Time of the day3–0.070.07–0.210.06–0.190.04
Orphan’s age when mother died0.040.06–0.080.17–0.050.14
Community size–0.080.25–0.590.41–0.50.33
Sex ratio0.040.21–0.360.49–0.280.39
Dominance rank0.020.17–0.330.33–0.270.28
Orphan’s age at sample–0.070.18–0.420.31–0.360.24
LCMS method (old)0.430.38–0.241.26–0.121.09
Sin(seasonDate)0.020.31–0.580.66–0.470.54
Cos(seasonDate)–0.010.03–0.080.06–0.070.05
Time of the day–0.010.03–0.080.06–0.070.04
Appendix 1—table 3
List of orphan immatures in the study and information about the adoption by adult individuals in the community.
IdentitySexCommunityAge when mother died (years)Adopted?*Duration of adoption*No. of samples in the studyAge when sampled for the study
BeatriceFEast4.9Yes>3 years385.3–8.6
EmmaFEast4.1Yes>2 years174.1–6.1
EolosMEast3.9Yes>3 years77.0–7.3
ErasmusMEast8.7Unknown318.8–11.8
FatimaFEast6.6Unknown87.8–10.1
GiaFEast2.6Yes17 months810.6–11.9
MaimounaFEast4.7No206.7–8.9
QuarantineFEast5.5Yes>2 years47.6–7.9
RichelieuMEast5.4Unknown6410.7–11.8
WillyMEast10.5Unknown5710.5–11.9
BalooFSouth3.8yes1 year227.3–8.5
CaramelMSouth7.3Unknown48.7 – 9.3
MohanFSouth4.1Yes1.5 year344.1 – 6.5
OscarMSouth4.7Yes>2 years498.1 – 11.9
WalaFSouth4.9Unknown258.1 – 10.9
RoxaneFNorth4.7Unknown410.5 – 11.3
VoltaFNorth3.8Unknown310.7 – 10.7
  1. *

    Data taken from Samuni et al., 2019a.

Additional files

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. Cédric Girard-Buttoz
  2. Patrick J Tkaczynski
  3. Liran Samuni
  4. Pawel Fedurek
  5. Cristina Gomes
  6. Therese Löhrich
  7. Virgile Manin
  8. Anna Preis
  9. Prince F Valé
  10. Tobias Deschner
  11. Roman M Wittig
  12. Catherine Crockford
(2021)
Early maternal loss leads to short- but not long-term effects on diurnal cortisol slopes in wild chimpanzees
eLife 10:e64134.
https://doi.org/10.7554/eLife.64134