Variation in thermal physiology can drive the temperature-dependence of microbial community richness

  1. Tom Clegg  Is a corresponding author
  2. Samraat Pawar
  1. Helmholtz Institute for Functional Marine Biodiversity at the University of Oldenburg (HIFMB), Germany
  2. Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research (AWI), Germany
  3. Department of Life Sciences, Silwood Park Campus, Imperial College London, United Kingdom
4 figures and 1 additional file

Figures

How variation in thermal physiology constrains microbial community species richness.

(A) Trait values increase with temperature following the Boltzmann-Arrhenius equation (boltz_maintext), with the shape governed by two parameters: B0 - trait value (r or a) at a reference temperature Tref and E - thermal sensitivity. (B) The joint distribution of E and log(B0) (here with empirically realistic negative covariance) determines how trait distributions vary across temperatures (C). (D) The distribution of trait values in turn determines the probability of feasibility Pfeas (and thus richness; Feas_sp_maintext). Specifically, Pfeas is determined by the proportion of relative growth rates (ri; blue shaded area) that are greater than the bound (solid black line) set by mean interaction strength (a). Populations with relative growth rates below this bound (red shaded area) are unfeasible (cannot persist in the community). All else being equal, the size of the unfeasible region (i.e. richness), decreases with increasing variance in the growth rate distribution (Var(ri)) and increasing interaction strength (which shifts the f(a(T)) bound upwards). (E–H) The effects of varying different aspects of the joint distribution of B0 and E of r and a on the emergent trait distribution across temperatures. Each panel shows the effect of altering the labeled parameter relative to the baseline case (far left), with inset plots showing the effect on the resulting temperature-richness relationship.

The effect of variation in trait thermal performance curves (TPCs) on the temperature-richness relationship in competitive microbial communities.

The analytical predictions (solid lines) are plotted along with the maximum richness reached in the numerical simulations (dots). (A) Mean thermal sensitivity of interactions a determines the direction and steepness of the temperature-richness relationship. (B) Increasing variance of thermal sensitivity increases unimodality. (C) Negative covariance between B0 and E shifts the peak of richness to higher temperatures. Parameter values used were: μr0=0.0,σr02=σa02=0.2,μEr=μEa=0.6,σEr2=σEa2=0.01,σB0,Er=σB0,Ea=0.0.

Sensitivity of theoretical results.

Bar plots show the Pearson correlation coefficient ρ=(x,y)/σxσy of each thermal physiology parameter with the root mean squared error between the theory and numerical simulations across replicate communities. Positive correlations indicate that increasing the parameter value tends to increase the error whilst negative values indicate error decreases the parameter value increases. Each panel shows the effect of a given parameter.

The bacterial temperature-richness relationship predicted by empirically-observed variation in thermal physiology.

(A) The relationship between log(B0) and E for growth rate in the experimental thermal performance curve (TPC) data from Smith et al., 2021. Dots show each pair of B0 and E values estimated for a given species/strain with histograms showing the marginal distributions. Ellipses show the 95% quantiles of the fitted bivariate normal distribution. (B) The actual growth-rate TPCs (solid lines) from the dataset as well as the fitted trait-distributions across temperature (box plots). The dashed line shows the point of minimum variance in growth rates which occurs towards the upper end of the temperature range. (C–D) Analogous plots for the dataset from the literature synthesis (Smith et al., 2019). (E) The analytically- (solid line) and simulation- (points) predicted temperature-richness curves based on the TPC variation seen in both these experimental (blue) and literature-synthesised (red) empirical data. Both are generated using the parameters from their respective fitted distributions and mean normalisation constants of μr0=0.0 and μa0=5.0. We set the normalisation constants such that the magnitude of richness values is not too large to perform the numerical simulations.

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. Tom Clegg
  2. Samraat Pawar
(2024)
Variation in thermal physiology can drive the temperature-dependence of microbial community richness
eLife 13:e84662.
https://doi.org/10.7554/eLife.84662