- Research article
- Open Access
- Published:

# Trait heritability in major transitions

*BMC Biology*
**volume 16**, Article number: 145 (2018)

## Abstract

### Background

Increases in biological complexity and the origins of life’s hierarchical organization are described by the “major transitions” framework. A crucial component of this paradigm is that after the transition in complexity or organization, adaptation occurs primarily at the level of the new, higher-level unit. For collective-level adaptations to occur, though, collective-level traits—properties of the group, such as collective size—must be heritable. Since collective-level trait values are functions of lower-level trait values, collective-level heritability is related to particle-level heritability. However, the nature of this relationship has rarely been explored in the context of major transitions.

### Results

We examine relationships between particle-level heritability and collective-level heritability for several functions that express collective-level trait values in terms of particle-level trait values. For clonal populations, when a collective-level trait value is a linear function of particle-level trait values and the number of particles per collective is fixed, the heritability of a collective-level trait is never less than that of the corresponding particle-level trait and is higher under most conditions. For more complicated functions, collective-level heritability is higher under most conditions, but can be lower when the environment experienced by collectives is heterogeneous. Within-genotype variation in collective size reduces collective-level heritability, but it can still exceed particle-level heritability when phenotypic variance among particles within collectives is large. These results hold for a diverse sample of biologically relevant traits.

### Conclusions

Rather than being an impediment to major transitions, we show that, under a wide range of conditions, the heritability of collective-level traits is actually higher than that of the corresponding particle-level traits. High levels of collective-level trait heritability thus arise “for free,” with important implications not only for major transitions but for multilevel selection in general.

## Background

Major transitions, or evolutionary transitions in individuality, are a framework for understanding the origins of life’s hierarchy and of biological complexity [1, 2]. During such a transition, a new unit of evolution emerges from interactions among previously existing units. This new unit, or collective, has traits not present before the transition and distinct from those of the units that comprise it (particles; see [3] for an in-depth discussion of collective-level traits). These collective-level traits are potentially subject to selection. Over the course of the transition, the primary level of selection shifts from the particle (lower-level unit) to the collective (higher-level unit), for example, from cells to multicellular organisms or from individual insects to eusocial societies.

Evolution by natural selection requires heritable variation in phenotypes that affect fitness at the level at which selection occurs [4, 5]. The breeder’s equation of quantitative genetics shows that heritability and strength of selection contribute equally to the adaptive response (see the “Analytical model” section below). When a collective-level trait is exposed to selection, it is collective-level heritability (the heritability of the collective-level trait) that determines the magnitude of the response. Collective-level heritability of traits is thus necessary for collective-level adaptations, but the emergence of collective-level heritability during a major transition has often been assumed to be difficult. For example, Michod considers the emergence of collective-level heritability through conflict mediation a crucial step in major transitions [2, 6, 7]. Simpson says that “From the view of some standard theory, these transitions are impossible,” in part because particle-level heritability greatly exceeds collective-level heritability [8].

Major transitions can be conceptualized as a shift from MLS1 to MLS2, in the sense of Damuth and Heisler [5], as in Okasha [9] (see also Godfrey-Smith [10], Shelton and Michod [11]). In MLS1, properties of the particles are under selection; in MLS2, it is the properties of the collectives. We follow Okasha [9] in referring to the lower-level units in a transition as ‘particles’ and the higher-level units as ‘collectives.’ Although our biological analogies are presented in terms of cells as particles and multicellular organisms as collectives, in principle our model could be extended to any pair of adjacent levels.

According to Michod [6], “…the challenge of ETI [evolutionary transitions in individuality] theory is to explain how fitness at the group level in the sense of MLS2 emerges out of fitness at the group level in the sense of MLS1.” But fitness, or selection, is only half of the breeder’s equation. Predicting the response to selection requires an estimate of heritability.

Whether or not collective-level fitness in MLS2 is a function of particle-level fitness is a matter of some disagreement (for example, Rainey and Kerr say no [12]). However, collective-level *phenotypes* must be functions of particle-level trait phenotypes, unless we accept strong emergence, a philosophical position tantamount to mysticism [13]. The function may be complex and involve cell-cell communication, feedbacks, environmental influences, etc., but it is still a function that is, in principle, predictable from particle-level trait values.

Nevertheless, the relationship between the heritability of particle-level traits and that of collective-level traits has rarely been considered in the context of major transitions, leading Okasha [14] to wonder, “Does variance at the particle level necessarily give rise to variance at the collective level? Does the heritability of a collective character depend somehow on the heritability of particle characters? The literature on multi-level selection has rarely tackled these questions explicitly, but they are crucial.” Similarly, Goodnight [15] says, “...we really do not have a good understanding of what contributes to group heritability, how to measure it, or even how to define it.”

While the role of selection has often been considered in the context of major transitions, the role of trait heritability has been relatively neglected. We examine relationships between particle-level heritability and collective-level heritability for several functions that express collective-level trait values in terms of particle-level trait values. For the simplest (linear) function, we derive an analytical solution for the relationship. For more complex functions, we employ a simulation model to explore the relationship over a range of conditions.

## Results

### Analytical model

There are several ways to estimate heritability, the proportion of phenotypic variation explained by genetic variation. If the strength of selection is known, heritability can be estimated by back-calculating from the breeder’s equation: *R* = *h*^{2}*S*, where *R* is the response to selection, *S* the selection differential, and *h*^{2} the narrow-sense heritability (i.e., the proportion of phenotypic variation explained by additive genetic variation). This can be rearranged as *h*^{2} = *S/R*. Another method is to compare parent and offspring trait values: the slope of the parent-offspring regression is an estimator of heritability [16]. We use the latter method in the simulations described in the next subsection.

Since heritability can be defined as the proportion of phenotypic variance explained by genetic variance, one method of estimation is to partition total variance into its components using an analysis of variance. We employ this approach in an analytical model to derive the relationship between the heritability of a collective-level trait and that of the particle-level trait from which it arises. For the sake of tractability, we begin with the simplest case, assuming that the size (number of particles) of collectives is fixed and that the collective-level trait value is a linear function of the particle-level trait values. We further assume that reproduction is asexual, so the proper measure of heritability is broad-sense heritability, *H*^{2}, which describes the proportion of phenotypic variation explained by all genetic variation [17].

We imagine a population in which collectives are made up of particles and genetically distinct clones are made up of collectives. As a concrete example, we can think of a population of undifferentiated volvocine algae, such as *Gonium*, in which case the particles are cells and the collectives are colonies. Because of asexual reproduction, many genetically-identical collectives may comprise a clone. Genetic variation among clones may arise through mutation or because the population is facultatively sexual, in which case these results will only hold for evolution within the asexual phase (in the *Gonium* example, during the summer bloom that precedes autumn mating and winter dormancy).

Broad-sense heritability is the ratio of genetic variance (*V*_{G}) to total phenotypic variance (*V*_{P}), estimated as the ratio of among-clone variance to total phenotypic variance [17]. Inherent in this concept is that genetically identical individuals are not always phenotypically identical; *V*_{P} includes both genetic and non-genetic variation. Non-genetic variation can arise from maternal effects, environmental (including microenvironmental) effects, and random developmental noise. Phenotypic variation among genetically identical individuals has been extensively documented, including in bacteria [18, 19], unicellular eukaryotes [20], plants [21], animals [17], and volvocine algae [22].

We used an ANOVA framework to estimate heritability as a ratio of sums of squares. Strictly speaking, heritability is a ratio of variances, not of sums of squares. However, the ratios of the relevant sums of squares converge to that of the variances as the number of categories increases, and for all but tiny or genetically uniform biological populations, the difference between the two ratios is negligible (see Additional file 1).

The ANOVA approach yields the following ratio of collective-level heritability, \( {H}_z^2 \), to particle-level heritability, \( {H}_y^2 \), the derivation of which is explained in the “Methods” section:

SSA is the sum of squared deviations from the mean among clones, SS(B/A) that among collectives within clones, and SS(C/B) that among particles within collectives. Given the assumptions of the model (linear function, constant particle number, clonal reproduction), collective-level heritability is therefore never less than particle-level heritability (i.e., the ratio of heritabilities is never less than 1) and is greater unless SS(C/B) = 0, in other words, unless particles within each collective have identical phenotype. Although we have derived this relationship assuming that the collective-level trait value is the average of particle-level trait values, the result holds for any linear function.

The results of the analytical model hold for all but tiny and/or extremely genetically depauperate populations. For example, for a population comprised of ten clones each with ten collectives, the approximation differs from the true value by less than 1%. The number of particles within a collective does not play a role, so our results are relevant even early in a major transition, when the collectives are likely to be small. For most real biological populations, the difference between the true heritability and the sums of squares approximation will be negligible (see Additional file 1 for a simple numerical example).

### Simulation model

The correspondence between particle-level and collective-level trait values is likely to be more complicated than a linear relationship for many interesting and biologically relevant cases. Here, we explore more complicated trait mapping functions using a simulation model. In this model, we consider two sources of non-genetic effects on particle phenotype (Fig. 1), each of which should lower the heritability of both particle- and collective-level traits. The first is intrinsic reproductive stochasticity in particle phenotype, analogous to developmental instability [23]. In the model, we determine the phenotype of daughter cells by sampling from a distribution centered on the parent’s genetic mean, with standard deviation *σ*. As shown in the analytical model above, by averaging out this variation, collectives can gain a heritability advantage over cells.

Our simulation also considers the phenotypic effects of environmental heterogeneity. Here, we model collectives as independently experiencing different environmental conditions that affect the phenotypes of all cells within them in the same manner. To extend the biological analogy offered above, *Gonium* colonies growing near the surface of a pond (where light and CO_{2} are abundant) may form colonies with larger cells than clonemates near the bottom. We implemented this in our model by assigning a size modifier, drawn from a normal distribution centered on 1 with standard deviation *σ*′, to each collective. We then multiplied the phenotype of each particle within the collective by this modifier. This source of phenotypic heterogeneity should reduce the heritability of collectives more than particles, simply because collectives experience a relatively higher frequency of stochastic events than particles do (each collective gets assigned a different size multiplier, but every particle within that collective experiences the same size multiplier).

We examine the effect of each of the above sources of phenotypic variation independently for the example of cells (particles) within nascent multicellular organisms (collectives). For a linear relationship, collective size is simply the sum of the sizes of cells within the collective. For both cells and collectives, heritability is assessed by calculating the slope of a linear regression on parent and offspring phenotype [16]. In this simple case, mean collective-level heritability is always greater than or equal to cell-level heritability. Only when *σ* = 0 (i.e., when all cells within a collective have identical phenotype) are cell- and collective-level heritability equal, in agreement with the analytical model. Greater developmental instability for cell size increases the advantage of collective-level heritability over cell-level heritability (Fig. 2a). Larger collectives, which average out cellular stochasticity more effectively, experience a greater increase in heritability than smaller collectives (Fig. 2a). Note that the simulations run in Fig. 2a reflect a very patchy environment in which environmental effects on cell size within collectives are large (*σ*′ = 0.25). While our model is not explicitly spatial, when *σ*′ is high, different collectives experience different environmental effects on their mean cell size, simulating the effects of a patchy environment. Increasing the magnitude of these environmental effects on cell size diminishes the difference in heritability between collectives and cells, but mean collective-level heritability is still greater than cell-level heritability for all parameter combinations (Fig. 2b). The simulations yielding Fig. 2a, b are provided as Additional files 2 and 3, respectively.

The volume of the cellular collective (Figs. 2 and 3a), which is simply the sum of the cell volumes within it, represents the simplest function mapping cellular to multicellular trait values. We now consider more complicated nonlinear functions relating cellular to multicellular trait values, some of which have biological relevance to the evolution of multicellularity.

The first nonlinear collective-level trait we consider is its diameter. Large size is thought to provide a key benefit to nascent multicellular collectives when they become too big to be consumed by gape-limited predators [24, 25]. For a collective that is approximately spherical, the trait that actually determines the likelihood of being eaten is diameter, which is therefore an important component of fitness. For geometric simplicity, we assume that the cells within the collective are pressed tightly together into a sphere. Collective volume (Fig. 3a) and diameter (Fig. 3b) exhibit similar dynamics, with collective-level heritability always exceeding cell-level heritability and being maximized under conditions of strong cell size stochasticity (high *σ*) and no environmental heterogeneity (low *σ*′).

Next, we consider swimming speed as a function of cell radius. We based this simulation on the hydrodynamics model of volvocine green algae derived by Solari et al. [26]. For simplicity, we modeled 32-celled, undifferentiated collectives (GS colonies in [26]), which would be similar to extant algae in the genus *Eudorina*. In this model, the swimming force of cells is independent of cell size, so, as cells get larger, the collective will become heavier (more negatively buoyant) without a corresponding increase in total swimming force, and therefore, its upward swimming speed will decrease. Thus, upward swimming speed is a monotonically declining function of cell radius (Fig. 3c inset), unlike the functions for volume and diameter (Fig. 3a, b insets), both of which are monotonically increasing. Nevertheless, the general behavior of heritability is very similar to the previous ones, and for a wide range of parameter values, the collective-level trait has a higher heritability than the cell-level trait (Fig. 3c).

Next, we consider a function describing a collective’s survival rate in the presence of a predator that can only consume collectives below a certain size. We calculated the survival rate (*c*) as a logistic function of the collective’s radius, effectively assuming that predation efficiency drops off quickly when collectives reach a threshold size (Fig. 3d inset). As with the previous functions (Fig. 3a–c), collective-level heritability is greater than cell-level heritability for much of the trait space and is maximized under conditions of high cellular stochasticity (*σ*) and low environmental heterogeneity (*σ*′; Fig. 3d).

Finally, we consider the case in which the simplifying assumption of constant cell number does not hold. Instead, the number of cells per collective fluctuates around the genetic mean \( \overline{N} \). In this case, each collective reproduces two new collectives, but the number of cells per new collective is a random variable drawn from a normal distribution with mean \( \overline{N} \) and coefficient of variation *CV*_{N} (the coefficient of variation for a normal distribution is the ratio of standard deviation to the mean). We chose to represent variation in the number of cells per collective as *CV*_{N} instead of standard deviation so that the scale of variation would not change with the size of the collective.

Variation in cell number, unlike the developmental and the environmental variation, does not affect the heritability of cells, only that of collectives. Therefore, we expected that increasing *CV*_{N} would decrease the ratio of collective-level to cell-level heritability. The simulation shows that the *CV*_{N} has a strong effect on collective-level heritability (Fig. 4). As *CV*_{N} increases, the ratio of collective- to cell-level heritabilities decreases, falling below one when the magnitude of *σ* is similar to or smaller than that of *CV*_{N} (Fig. 4).

## Discussion

When particles (cells, for example) first form collectives (multicellular colonies, for example), traits inevitably come into being that did not previously exist [3]. Some such traits, such as collective size, will be simple functions of particle-level traits; others will be more complex functions. We are interested in how these new collective-level traits respond to selection, and the breeder’s equation tells us that the response is proportional to the trait’s heritability (for a given strength of selection). We can predict that heritability if we know two things: the heritability of the underlying individual-level trait and relationship between particle-level and collective-level heritability. Estimating the former is straightforward; if we can define the latter, we can in principle predict the efficacy of selection on the new trait.

Using a quantitative genetics framework, we have derived an analytical solution for the relationship between particle-level and collective-level heritability in clonal populations for a limited case. When particle number is constant and the collective-level trait value is a linear function of the particle-level trait values, the organismal heritability turns out to be a simple function of the cell-level heritability. In contrast to claims that particle-level heritability is always higher than collective-level heritability (e.g., [8]), we have shown that collective-level heritability is higher over a wide range of conditions. Because this result depends on the number of clones and the number of colonies within a clone, it may not hold for very small populations or those with little genetic variation. This is not a major limitation, though, since tiny, genetically homogeneous populations are unlikely to be the ones experiencing selectively driven evolutionary transitions in individuality.

This analytical result is a step toward understanding the relationship between heritabilities at two adjacent hierarchical levels, but the assumptions of constant particle number and linear function are restrictive. The simulation model shows that the results are somewhat dependent on the function relating the trait values at the two levels. However, these functions were chosen to be diverse, and the behavior of the relative heritabilities is nevertheless qualitatively similar, increasing with cellular developmental variation (*σ*)*,* decreasing with environmental heterogeneity (*σ*′), and exceeding 1 for most of the parameter space.

Of course, we have not (and cannot) comprehensively explored the universe of possible functions relating collective-level traits to particle-level traits. What we have done is explore a small sample of this space, with functions ranging from extremely simple (volume) to somewhat more complex (swimming speed, survival under predation). We do not claim that the high heritabilities estimated for these collective-level traits would apply to all such traits, and a full accounting of possible functions is beyond the scope of this (or any) study. Rather, we have shown that for at least some such functions, the resulting collective-level traits can have high heritability, and thus be altered by selection, early in an evolutionary transition in individuality.

All four of the collective-level traits in the simulation models are potentially biologically relevant. Volume and diameter are both aspects of size, which can be an important component of fitness both in evolutionary transitions in individuality [27] and in life history evolution [28]. Swimming speed is a measure of motility, which has selective consequences for a wide range of organisms, including many animals and microbes. For planktonic organisms, a positive upward swimming speed provides active control of depth, allowing some control over light intensity (for autotrophs) and prey abundance (for heterotrophs). Survival under predation obviously has important fitness implications for many organisms, and both theoretical and experimental evidence implicate predation as a possible selective pressure driving the evolution of multicellularity. Kirk, for example, suggests that a “predation threshold” above which algae are safe from many filter feeders may have driven the evolution of multicellularity in the volvocine algae [29]. Microbial evolution experiments in the algae *Chlorella* and *Chlamydomonas* have shown that predation can drive the evolution of undifferentiated multicellular clusters [30,31,32].

In our simulations, we examined the effects of three independent sources of phenotypic variation affecting the relative heritability of particle and collective-level traits. Stochastic variation in cell size around the clone’s genetic mean (*σ*) reduces the absolute heritability of cells and collectives by introducing non-heritable phenotypic variation. By averaging across multiple cells, however, collectives reduce the effects of this phenotypic variation, providing them with a relative heritability advantage over cells.

We also considered the effect of environmental heterogeneity in which all of the cells within a collective are affected in the same manner (*σ’*). Collectives are disproportionately affected: each collective is assessed a different size modifier, but all of the cells within these collectives are affected in the same manner. As a result, collectives experience *n-*fold more stochastic events (where *n* is the number of cells per collective), which reduces their heritability relative to cells. The influence of these sources of variation is evident in the contour plots of Fig. 3: the relative heritability of collectives to cells is maximized when the cellular stochastic variation is high and environmental heterogeneity low (lower right corner of the plots).

Finally, we considered variation in the number of particles per collective. Such variation substantially reduces the heritability of a collective-level trait. Even with reasonably large variation in collective size, though, the collective-level trait retains most of the heritability of the particle-level trait on which it is based (for example, ~ 55% at a *CV*_{N} in particle number of 0.25).

A large number of previous studies have addressed heritability in the context of multilevel or group selection. Heritability above the individual level has been called group heritability (e.g., [33]), populational heritability (e.g., [34]), community heritability (e.g., [35]), and heritability of the family mean [16]. These prior treatments differ from ours in one or more of the following respects: they are concerned with the evolution of individual-level traits rather than group level traits (particle- rather than collective-level in our terminology), are based on MLS1 rather than MLS2 models, and are focused on narrow-sense rather than broad-sense heritability. Furthermore, few previous studies have addressed these questions in the context of the major transitions. Without attempting a comprehensive review, we summarize several such studies, and important differences from our own, below.

Queller [33] presents a useful reformulation of the Price equation for selection at two levels:

in which \( \varDelta \overline{G} \) is the change in average trait value, *S*_{b} and *S*_{w} are the selection differentials between collectives and within collectives, respectively, and *h*^{2}_{b} and *h*^{2}_{w} are the heritabilities of the collective-level and individual-level traits, respectively. This formulation partitions the response to selection on a particle-level trait into within- and among-collective change, but the focus is still on particle-level traits. Our focus is on the evolution of collective-level traits. In the terminology of Damuth and Heisler [5], our focus is on MLS2, while Queller’s is on MLS1. In addition, Queller makes no attempt to derive the relationship between collective-level heritability and particle-level heritability.

Michod and Roze [2] have previously modeled the relationship between particle-level and collective-level heritability of fitness during a major transition. However, as Okasha [14] points out, the heritability of fitness only ensures that mean population fitness will increase over time. For selection to result in directional phenotypic change, it is phenotypes that must be heritable. Furthermore, Michod and Roze focused on within-organism genetic change. Our models assume that such change is negligible, as is likely to be true early in a transition, when collectives (e.g., nascent multicellular organisms) presumably include a small number of clonally replicating particles (e.g., cells).

Okasha [36] considers heritability in MLS1 (which he refers to as group selection 2) and MLS2 (his group selection 1) but does not attempt to derive a relationship between heritabilities at two levels. Knowing the ratio of heritabilities is necessary, though not sufficient, to predict the outcome of opposing selection at two levels and so has important implications for collective-level traits that arise from cooperation among particles. The presumed higher heritability of the particle-level traits has been seen as a problem for the evolution of cooperation that benefits the collective [2, 8, 37,38,39]. Our results show that this problem does not always exist, though we would need to know not only the relative heritabilities but also the relative strengths of selection to predict the outcome of opposing selection at two adjacent levels.

Several previous papers have shown that group-level heritability (collective-level heritability in our terminology) exists and can be substantial. Slatkin [34], for example, showed that one measure of group-level heritability, fraction of total variance between lines, is substantial both in an analytical model and in the *Tribolium* experiments of Wade and McCauley [40]. Under some conditions, the between-line variance of a linear trait such as the one we consider in our analytical model exceeds the within-line variance.

Bijma et al. [41, 42] and Wade et al. [43] showed that variance in the total breeding value of a population can be increased, even to the point of exceeding phenotypic variance, by interactions among individuals. Our model does not consider (or require) interactions among individuals. Further, their model and empirical example are exclusively concerned with individual-level traits (particle-level traits in our terminology), for example, survival days in chickens. They do not estimate group heritability as such and judge that “it is unclear how this parameter should be defined or estimated.”

Goodnight [44] addresses the effect of environmental variance at two scales on the responses to individual and group selection in *Arabidopsis thaliana*. Although Goodnight’s study focused on an individual-level trait (leaf area, thus an MLS1 scenario) of an obligately sexual organism (thus narrow-sense heritability), our results (shown in Fig. 3) showed analogous effects. In both cases, environmental variation at a fine scale (analogous to our *σ*) increased the efficacy of group selection relative to individual selection, while environmental variation among demes (analogous to our *σ*′) had the opposite effect. Wade [45] examines a similar case from a theoretical perspective and finds that increased relative efficacy of group selection results from environmental variation among particles decreasing particle-level heritability.

Goodnight [15] considers the ratio of group-level heritability to individual-level heritability (in the narrow sense) using contextual analysis. Although this paper does not provide a formula to calculate this ratio, its inequality 5 sets a minimum bound (with the assumption that selection at the two levels is in opposition). As in our analyses, Goodnight shows that group-level heritability can exceed individual-level heritability in some circumstances.

Several simplifying assumptions underlie our models, most importantly the genetic identity of particles within collectives. This condition only applies to a subset of the major transitions. Queller recognized two subcategories within Maynard Smith and Szathmáry’s [1] list of transitions, which he called “egalitarian” and “fraternal” transitions [46]. Briefly, egalitarian transitions involve particles that may be genetically distinct, or even from different species, such as the alliance of a bacterium with an Archaean that gave rise to the eukaryotic cell. Fraternal transitions are those in which the particles are genetically similar or identical, such as the origins of eusociality and of most multicellular lineages.

Because of our assumptions of asexual reproduction and genetic identity among particles, we cannot generalize our results to all types of major transitions. Egalitarian transitions will not normally meet this criterion. A possible exception is aggregative multicellularity, as seen in cellular slime molds and myxobacteria, when assortment is so high that fruiting bodies are genetically uniform. This is probably uncommon [47], but it does happen [48, 49]. Transitions in which reproduction of particles is obligately sexual, such as the origins of eusociality, also violate this assumption.

A better fit for our models is clonal multicellularity, which is probably the most common type of major transition. An incomplete list of independent origins of clonal multicellularity includes animals; streptophytes; chytrid, ascomycete, and basidiomycete fungi; florideophyte and bangiophyte red algae; brown algae; peritrich ciliates; ulvophyte green algae; several clades of chlorophyte green algae; and filamentous cyanobacteria [50,51,52,53]. In most cases, the early stages in these transitions probably violated the assumption of uniform particle number per collective, but our simulations show that our main results are robust to reasonable violations of this assumption.

One example that does approximate all of our assumptions is that of the volvocine green algae, an important model system for understanding the evolution of multicellularity. Volvocine algae undergo clonal reproduction only occasionally punctuated by sex, are small enough that within-collective mutation probably has negligible phenotypic effects, and have cell numbers that are under tight genetic control.

## Conclusions

A great deal of work has gone into understanding the selective pressures that may have driven major evolutionary transitions. However, heritability is just as important as the strength of selection in predicting evolutionary outcomes. We have shown that, given some simplifying assumptions, heritability of collective-level traits comes “for free;” that is, it emerges as an inevitable consequence of group formation. Qualitatively, this result holds across a wide range of parameters and for a diverse sample of biologically relevant traits. Collective-level heritability is maximized (relative to particle-level heritability) when phenotypic variation among particles is high and when environmental heterogeneity and variation in collective size are low. Understanding the emergence of trait heritability at higher levels is necessary to model any process involving multilevel selection, so our results are relevant to a variety of other problems.

## Methods

### Analytical model

Treating particles and collectives separately, the phenotype of particle *k* in collective *j* within clone *i* can be expressed as

where *m* is the mean genotypic value of all clones, *A*_{i} is the deviation of clone *i* from *m*, *B*_{j(i)} is the deviation of collective *j* from the mean of clone *i*, and *C*_{k(ij)} is the deviation of particle *k* from the mean of collective *j* within clone *i*. The model in (2) describes a nested ANOVA framework, in which the sums of squared deviations from the population mean are partitioned into among-clone, among collectives within clone, and within-collective components. The among-clone component, the sum of squared deviations of *A* from *m*, is

where *a, b,* and *c* are the number of clones, collectives within a clone, and particles within a collective, respectively. The sum of squared deviations of collectives within clones is

that among particles within collectives is

and the total sum of squares is

Broad-sense heritability of a particle-level trait, \( {H}_y^2 \), is the ratio of genetic variance to total phenotypic variance:

We now turn our attention to collective-level traits. The phenotype of collective *j* within clone *i* can be expressed as

where *μ* is the mean genetic value of all clones, *α*_{i} is the deviation of clone *i* from *μ*, and *β*_{j(i)} is the deviation of collective *j* from the mean of clone *i*. The sum of squared deviations of *α* from *μ* is

The sum of squares among colonies within clones is

and the total sum of squares is

Broad-sense heritability of a collective-level trait, \( {H}_z^2 \), is the ratio of genetic variance to total phenotypic variance,

If collective-level trait value is the average of cell-level trait values, *z*_{ij} = *y*_{ij∙}, \( {\overline{z}}_{i\bullet }={\overline{y}}_{i\bullet \bullet } \), and \( {\overline{z}}_{\bullet \bullet }={\overline{y}}_{\bullet \bullet \bullet } \). Thus SSα = *c*SSA, and SS(β/α) = *c*SS(B/A). Substituting into (12),

we get

The ratio of collective-level heritability to particle-level heritability is thus

This result holds for any linear function. The substitution that gets us from (12) to (13) introduces the constant *c*, which scales both numerator and denominator and therefore cancels out. Different linear functions would change the magnitude of the constant relating SSα to *c*SSA and SS(β/α) to *c*SS(B/A) but not the fact that numerator and denominator are scaled by the same constant.

### Simulation model

All simulation models were carried out using custom Python scripts. As in the analytical model, particles grow in clonal collectives, which reproduce by forming two new collectives, each with as many particles as its parent. The initial population is founded by ten genetically distinct clones, each of which has a different genetically determined mean particle phenotype (spaced evenly between 1 and 2). These are grown for at least seven generations, resulting in at least 127 collective-level reproductive events per genotype and 127*n* (where *n* is particle number per collective) particle-level reproductive events per genotype. For each function, we calculated the relative heritability of collective- to cell-level traits for 32-celled collectives across 1024 combinations of *σ* and *σ*′ ranging from 0 to 0.25.

For the simulation of collective volume, we assume a linear relationship between cell volume and collective volume, specifically that collective volume is the sum of cell volumes. The simulation of collective volume is provided as Additional file 4.

For the simulation of collective diameter, we assume that the particles within the collective are pressed tightly together into a sphere, allowing us to calculate collective diameter as \( d=2{\left(\frac{3V}{4\pi}\right)}^{\frac{1}{3}} \), where *V* is the sum of the cell volumes within the collective. The simulation of collective diameter is provided as Additional file 5.

For the simulation of swimming speed, the function relating cell radius to upward swimming speed in Solari et al. [26] (Eq. 4 from [26]) can be simplified to

where *f* is average effective upward swimming force per cell, *N* is the number of cells per collective, *η*_{w} is water viscosity, *r* is the average radius of cells in the collective, and *Δρ*_{c} is the density difference between cells and water. Additional file 6 provides a more detailed description of the derivation of Eq. 15.

Using the numerical values in Solari et al [26], *η*_{w} = 0.01 g/cm·s, *Δρ*_{c} = 0.047 g/cm^{3}, and *f* = 2.4 × 10^{−7} g cm/s^{2}, so

The simulation of swimming speed is provided as Additional file 7.

For the simulation of survival rate in the presence of a predator, we calculated the survival rate (*c*) as a logistic function of the collective’s radius:

The simulation of survival under predation is provided as Additional file 8.

For the simulation of variable cell number, we calculated the relative heritability of size (volume) for collectives and cells across 1024 combinations of *σ* and *CV*_{N} ranging from 0 to 0.25 with a mean cell number (\( \overline{N}\Big) \) of 32. The simulation of variable cell number is provided as Additional file 9.

## References

- 1.
Maynard Smith J, Szathmáry E. The major transitions in evolution. Oxford: Oxford University Press; 1995.

- 2.
Michod RE, Roze D. Transitions in individuality. Proc Biol Sci. 1997;264:853–7.

- 3.
Ratcliff WC, Herron M, Conlin PL, Libby E. Nascent life cycles and the emergence of higher-level individuality. Philos Trans R Soc B Biol Sci. 2017;372:20160420.

- 4.
Lewontin RC. The units of selection. Annu Rev Ecol Syst. 1970;1:1–18.

- 5.
Damuth J, Heisler IL. Alternative formulations of multilevel selection. Biol Philos. 1988;3:407–30.

- 6.
Michod RE. On the transfer of fitness from the cell to the multicellular organism. Biol Philos. 2005;20:967–87.

- 7.
Michod RE. Evolutionary transitions in individuality: multicellularity and sex. In: Calcott B, Sterelny K, editors. The major transitions evolution revisited. Cambridge: MIT Press; 2011. p. 167–97.

- 8.
Simpson C. How many levels are there? How insights from evolutionary transitions in individuality help measure the hierarchical complexity of life. In: Calcott B, Sterelny K, editors. The major transitions evolution revisited. Cambridge: MIT Press; 2011. p. 199–225.

- 9.
Okasha S. Multilevel selection and the major transitions in evolution. Philos Sci. 2005;72:1013–25.

- 10.
Godfrey-Smith P. Darwinian populations and natural selection. Oxford: Oxford University Press; 2009.

- 11.
Shelton DE, Michod RE. Philosophical foundations for the hierarchy of life. Biol Philos. 2010;25:391–403.

- 12.
Rainey PB, Kerr B. Conflicts among levels of selection as fuel for the evolution of individuality. In: Calcott B, Sterelny K, editors. The major transitions evolution revisited. Cambridge: MIT Press; 2011. p. 141–62.

- 13.
Bedau MA. Weak emergence. Noûs. 1997;31:375–9.

- 14.
Okasha S. Evolution and the levels of selection. Oxford: Oxford University Press; 2006.

- 15.
Goodnight CJ. Multilevel selection: the evolution of cooperation in non-kin groups. Popul Ecol. 2005;47:3–12.

- 16.
Falconer DS. Introduction to quantitative genetics. New York: Ronald Press Co.; 1960.

- 17.
Lynch M, Walsh B. Genetics and analysis of quantitative traits. Sunderland: Sinauer Associates, Inc.; 1998.

- 18.
Elowitz MB, Levine AJ, Siggia ED, Swain PS. Stochastic gene expression in a single cell. Science. 2002;297:1183–6.

- 19.
McAdams HH, Arkin A. Stochastic mechanisms in gene expression. Proc Natl Acad Sci U S A. 1997;94:814–9.

- 20.
Blake WJ, Kærn M, Cantor CR, Collins JJ. Noise in eukaryotic gene expression. Nature. 2003;422:633–7.

- 21.
Burton GW, DeVane EH. Estimating heritability in tall fescue (

*Festuca Arundinaceae*) from replicated clonal material. Agron J. 1953;45:478–81. - 22.
Herron MD, Ghimire S, Vinikoor CR, Michod RE. Fitness trade-offs and developmental constraints in the evolution of soma: an experimental study in a volvocine alga. Evol Ecol Res. 2014;16:203–21.

- 23.
Palmer AR. Fluctuating asymmetry analyses: a primer. In: Markow TA, editor. Developmental Instability: Its Origins and Evolutionary Implications. Dordrecht: Kluwer; 1994. p. 335–64.

- 24.
Stanley SM. An ecological theory for the sudden origin of multicellular life in the late Precambrian. Proc Natl Acad Sci U S A. 1973;70:1486–9.

- 25.
Bell G. The origin and early evolution of germ cells as illustrated by the Volvocales. In: Halvorson HO, Monroy A, editors. The origin and evolution of sex. New York: Alan R. Liss; 1985.

- 26.
Solari CA, Kessler JO, Michod RE. A hydrodynamics approach to the evolution of multicellularity: flagellar motility and germ-soma differentiation in volvocalean green algae. Am Nat. 2006;167:537–54.

- 27.
Bonner JT. Perspective: the size-complexity rule. Evolution. 2004;58:1883–90.

- 28.
Peters RH. The ecological implications of body size. Cambridge: Cambridge University Press; 1983.

- 29.
Kirk DL. Volvox: molecular-genetic origins of multicellularity. Cambridge: Cambridge University Press; 1998.

- 30.
Boraas ME, Seale DB, Boxhorn JE. Phagotrophy by a flagellate selects for colonial prey: a possible origin of multicellularity. Evol Ecol. 1998;12:153–64.

- 31.
Herron MD, Borin JM, Boswell JC, Walker J, Chen I-CK, Knox CA, et al. De novo origins of multicellularity in response to predation. bioRxiv. 2018. https://doi.org/10.1101/247361.

- 32.
Becks L, Ellner SP, Jones LE, Hairston NG. Reduction of adaptive genetic diversity radically alters eco-evolutionary community dynamics. Ecol Lett. 2010;13:989–97.

- 33.
Queller DC. Quantitative genetics, inclusive fitness, and group selection. Am Nat. 1992;139:540–58.

- 34.
Slatkin M. Populational heritability. Evolution. 1981;35:859–71.

- 35.
Shuster SM, Lonsdorf EV, Wimp GM, Bailey JK, Whitham TG. Community heritability measures the evolutionary consequences of indirect genetic effects on community structure. Evolution. 2006;60:991–1003.

- 36.
Okasha S. The concept of group heritability. Biol Philos. 2003;18:445–61.

- 37.
Michod RE. Evolution of the individual. Am Nat. 1997;150(Suppl):S5–21.

- 38.
Buss LW. The uniqueness of the individual revisited. In: Jackson JBC, Buss LW, Cook RE, editors. Population biology and evolution of clonal organisms; 1986. p. 467–505.

- 39.
Buss LW. The evolution of individuality. Princeton: Princeton University Press; 1987.

- 40.
Wade MJ, McCauley DE. Group selection: the phenotypic and genotypic differentiation of small populations. Evolution. 1980;34:799–812.

- 41.
Bijma P, Muir WM, Ellen ED, Wolf JB, Van Arendonk JAM. Multilevel selection 2: estimating the genetic parameters determining inheritance and response to selection. Genetics. 2007;175:289–99.

- 42.
Bijma P, Muir WM, Van Arendonk JAM. Multilevel selection 1: quantitative genetics of inheritance and response to selection. Genetics. 2007;175:277–88.

- 43.
Wade MJ, Bijma P, Ellen ED, Muir W. Group selection and social evolution in domesticated animals. Evol Appl. 2010;3:453–65.

- 44.
Goodnight CJ. The influence of environmental variation on group and individual selection in a cress. Evolution. 1985;39:545–58.

- 45.
Wade M. Adaptation in metapopulations: how interaction changes evolution. Chicago: University of Chicago Press; 2016.

- 46.
Queller DC. Cooperators since life began. Q Rev Biol. 1997;72:184–8.

- 47.
Sathe S, Kaushik S, Lalremruata A, Aggarwal RK, Cavender JC, Nanjundiah V. Genetic heterogeneity in wild isolates of cellular slime mold social groups. Microb Ecol. 2010;60:137–48.

- 48.
Gilbert OM, Foster KR, Mehdiabadi NJ, Strassmann JE, Queller DC. High relatedness maintains multicellular cooperation in a social amoeba by controlling cheater mutants. Proc Natl Acad Sci U S A. 2007;104:8913–7.

- 49.
Gilbert OM, Queller DC, Strassmann JE. Discovery of a large clonal patch of a social amoeba: implications for social evolution. Mol Ecol. 2009;18:1273–81.

- 50.
Parfrey LW, Lahr DJG. Multicellularity arose several times in the evolution of eukaryotes. BioEssays. 2013;35:339–47.

- 51.
Niklas KJ, Newman SA. The origins of multicellular organisms. Evol Dev. 2013;15:41–52.

- 52.
Herron MD, Rashidi A, Shelton DE, Driscoll WW. Cellular differentiation and individuality in the “minor” multicellular taxa. Biol Rev. 2013;88:844–61.

- 53.
Melkonian M. Phylum Chlorophyta. In: Margulis L, Corliss JO, Melkonian M, Chapman DJ, editors. Handbook of Protoctista. Boston: Jones and Bartlett Publishers; 1990. p. 597–660.

## Acknowledgements

We would like to thank Sam Brown, Peter Conlin, Michael Doebeli, David Garton, Rick Michod, and Deborah Shelton for helpful discussions.

### Funding

This work was supported by grants from NASA (NNA17BB05A, NNX15AR33G), NSF (DEB-1457701, DEB-1723293, DEB-1456652), and the John Templeton Foundation (43285).

### Availability of data and materials

All data generated or analyzed during this study are included in this published article [and its Additional files 1, 2, 3, 4, 5, 6, 7, 8, and 9].

## Author information

### Affiliations

### Contributions

MDH conceived the project, developed the analytical model, contributed to the simulation models, and contributed to writing the manuscript. SAZ-D and WCR developed the simulation models and contributed to writing the manuscript. All authors read and approved the final manuscript.

### Corresponding author

Correspondence to Matthew D. Herron.

## Ethics declarations

### Ethics approval and consent to participate

Not applicable

### Consent for publication

Not applicable

### Competing interests

The authors declare that they have no competing interests.

### Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Additional files

### Additional file 1:

Ratios of sums of squares as approximations for ratios of variances. Justifies the use of sums of squares as approximations for variances and gives a numerical example showing that inaccuracies introduced by these approximations will normally be negligible. (DOCX 20 kb)

### Additional file 2:

Python script for fixed effect of the environment (standard deviation *σ*′ = 0.25), and varying degree of developmental instability for particle size *σ*. This script generates Fig. 2a. (PY 9 kb)

### Additional file 3:

Python script for fixed developmental instability (*σ* = 0.25) and varying between-collective environmental effects on cell size. This script generates Fig. 2b. (PY 9 kb)

### Additional file 4:

Python script for simulation of collective volume. This script generates Fig. 3a. (PY 13 kb)

### Additional file 5:

Python script for simulation of collective diameter. This script generates Fig. 3b. (PY 13 kb)

### Additional file 6:

Detailed description of the derivation of Eq. 15. (DOCX 25 kb)

### Additional file 7:

Python script for simulation of swimming speed. This script generates Fig. 3c. (PY 13 kb)

### Additional file 8:

Python script for simulation of survival under predation. This script generates Fig. 3d. (PY 13 kb)

### Additional file 9:

Simulation of variable cell number. This script generates Fig. 4. (PY 12 kb)

## Rights and permissions

**Open Access** This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

## About this article

#### Received

#### Accepted

#### Published

#### DOI

### Keywords

- Evolution
- Heritability
- Major transitions
- Multicellularity
- Quantitative genetics
- Simulations