Variation in thermal physiology can drive the temperature-dependence of microbial community richness
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: - trait value ( or ) at a reference temperature and - thermal sensitivity. (B) The joint distribution of and (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 (and thus richness; Feas_sp_maintext). Specifically, is determined by the proportion of relative growth rates (; blue shaded area) that are greater than the bound (solid black line) set by mean interaction strength (). 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 () and increasing interaction strength (which shifts the bound upwards). (E–H) The effects of varying different aspects of the joint distribution of and of and 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 determines the direction and steepness of the temperature-richness relationship. (B) Increasing variance of thermal sensitivity increases unimodality. (C) Negative covariance between and shifts the peak of richness to higher temperatures. Parameter values used were: .

Sensitivity of theoretical results.
Bar plots show the Pearson correlation coefficient 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 and for growth rate in the experimental thermal performance curve (TPC) data from Smith et al., 2021. Dots show each pair of and 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 and . We set the normalisation constants such that the magnitude of richness values is not too large to perform the numerical simulations.