Division of Labour: Losing out to improve group fitness

A mathematical model provides clues as to why members of a group divide tasks between them even when specialisation reduces the performance of individuals.
  1. Jos Kramer
  2. Rolf Kümmerli  Is a corresponding author
  1. Department of Quantitative Biomedicine, University of Zurich, Switzerland

In The Wealth of Nations, political economist Adam Smith stated that division of labour is favoured if specialisation leads to an “increase of dexterity in every particular workman” so that the same number of people show a “great increase of the quantity of work". Simply speaking, individuals divide tasks between them when specialisation leads to accelerated returns for everybody involved (Bourke, 2011; Cooper and West, 2018; Michod, 2006). While this concept was initially devised for human societies, it turned out to be applicable to biological systems.

Examples of individuals within social groups specialising in different tasks are plentiful in nature and cover the entire range of biological complexity: from subgroups of bacteria producing different compounds of a biofilm (Dragoš et al., 2018), to ant queens specialising in egg laying and ant workers in brood care (Hölldobler and Wilson, 1990). This scientific concept is intuitive, easy to grasp and solves the challenge of how division of labour can evolve. But, does it cover all scenarios of how division of labour can arise? A recent study showed that in groups where individuals are organized into a particular spatial arrangement (also known as topology) and can only cooperate with direct neighbours, division of labour can evolve even when specialisation makes individuals less rather than more efficient (Yanni et al., 2020). New provoking hypotheses often spur controversy, and these findings are no exception. After all, to work harder but get less out of it seems, at first glance, not compatible with evolutionary theory.

Now, in eLife, Guy Cooper and colleagues from the University of Oxford – Hadleigh Frost, Ming Liu and Stuart West – report the results of a new mathematical model probing this hypothesis (Cooper et al., 2021). The model explored the conditions required for a population of cells, similar to the social groups reported by Rueffler et al., 2012, to divide reproduction and helping between group members. Each cell can differentiate into either a helper, which sacrifices its fecundity to increase the viability of other members, or a reproductive, which has higher fecundity but is less cooperative.

Cooper et al. found that reproductive division of labour can evolve when it improves the overall fitness of the group (Figure 1) and arises through three scenarios: (1) specialisation that results in accelerated returns for individuals; (2) physical characteristics that predispose some cells to become helpers and others to become reproductives; or (3) when specialisation is reciprocal and synergistic, meaning that by losing fitness helpers disproportionally benefit reproductives. The first scenario is the classic case, but how relevant are the other two, and how do they link to group topology?

Simplified model of how division of labour evolves between cells in spatially organized groups.

In the model created by Cooper et al. interactions occur either between helper (yellow) and reproductive (blue) cells, or among generalists (grey). Here, the representative model created for this figure assumes that generalists invest equally into help (h = 1) and personal fecundity (f = 1), whereas specialised cells invest fully into either helping (h = 2; f = 0) or fecundity (f = 2; h = 0). Among generalists (middle), the eventual fecundity of a cell (F) is increased through help from neighbours: in this example, F = f + 0.5 n, where n is the number of neighbours. Reproductives are more efficient than generalists in converting help to fecundity, and this synergistic benefit leads to a higher final fecundity: F = f + n*s, where n represents the number of neighbours and s represents the synergistic benefit. For example, when cells are arranged into a single column, also known as a filament structure (top), and their differentiation is coordinated (left), the fecundity of reproductives (f = 2), which are neighboured by two helpers, increases to five when s = 1.5; the fecundity and help values of individual cells are indicated above and below them, respectively. In filament structures, and also branched structures (bottom), this coordinated differentiation results in the total fecundity (green box) and viability (red box) of the group being higher than in populations that did not specialise and divide their labour (middle). However, when differentiation is not coordinated (right), division of labour decreases total group-level fecundity.

Image credit: Figure created using BioRender.com

Cooper et al. used their mathematical model to investigate three topologies: branching, when cells form tree-like structures; filaments, when cells form a single column; and well-mixed group structures, when each cell is connected to all other cells in the group. They found that division of labour can evolve in all three topologies even when specialisation leads to diminishing returns for individuals (that is, instances where the more support helpers provide the less efficient they are at helping). However, the underlying forces that drive this evolution differ between the three topologies.

In branching topologies, division of labour with diminishing returns arises from physical differences related to where cells are positioned in the group (scenario two): cells with fewer neighbours have reduced costs from refraining to reproduce and specializing into helpers. Conversely, in one-dimensional filaments, division of labour can evolve via reciprocal specialisation (scenario three) provided that the support given to reproductives increases their fecundity more than helpers have given up by specialising. Notably, this scenario also works in well-mixed groups, but only under the stringent assumption that helpers provide extreme levels of support. Altogether, these findings reveal that specific topological constraints are not required for reproductive division of labour to evolve, but allow it to happen under a larger range of conditions – that is, scenarios that differ in terms of the amount of help provided and how much it generates in return.

There is one important caveat with these results: the models assume that individuals can coordinate who differentiates into a reproductive or a helper cell. Without this coordination, division of labour turns out to be disastrous and reduces the overall reproductive fitness of the group (Figure 1; right). Cooper et al. note that this might limit the scope of spatial topology favouring division of labour in instances where returns for individuals are reduced. After all, coordination is more likely to evolve after division of labour has occurred and cells have already developed specific roles (Liu et al., 2021).

The insights by Cooper et al. are conclusive, novel and expand theoretical knowledge on how division of labour evolved. Nonetheless, empirical scientists might ask whether this is a purely theoretical exercise and whether there are real-world examples for scenarios involving non-accelerating returns. As this is currently unclear, Cooper et al. have provided guidelines for how to test for these scenarios in biological systems which will likely lead to examples in nature surfacing in the near future.


  1. Book
    1. Hölldobler B
    2. Wilson EO
    The Ants
    Harvard University Press.

Article and author information

Author details

  1. Jos Kramer

    Jos Kramer is in the Department of Quantitative Biomedicine, University of Zurich, Zurich, Switzerland

    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-2155-9295
  2. Rolf Kümmerli

    Rolf Kümmerli is in the Department of Quantitative Biomedicine, University of Zurich, Zurich, Switzerland

    For correspondence
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-4084-6679

Publication history

  1. Version of Record published: December 17, 2021 (version 1)
  2. Version of Record updated: January 26, 2022 (version 2)


© 2021, Kramer and Kümmerli

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.


  • 1,089
  • 101
  • 0

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

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. Jos Kramer
  2. Rolf Kümmerli
Division of Labour: Losing out to improve group fitness
eLife 10:e75243.
  1. Further reading

Further reading

    1. Evolutionary Biology
    2. Genetics and Genomics
    Brian PH Metzger, Yeonwoo Park ... Joseph W Thornton
    Research Article

    A protein’s genetic architecture – the set of causal rules by which its sequence produces its functions – also determines its possible evolutionary trajectories. Prior research has proposed that the genetic architecture of proteins is very complex, with pervasive epistatic interactions that constrain evolution and make function difficult to predict from sequence. Most of this work has analyzed only the direct paths between two proteins of interest – excluding the vast majority of possible genotypes and evolutionary trajectories – and has considered only a single protein function, leaving unaddressed the genetic architecture of functional specificity and its impact on the evolution of new functions. Here, we develop a new method based on ordinal logistic regression to directly characterize the global genetic determinants of multiple protein functions from 20-state combinatorial deep mutational scanning (DMS) experiments. We use it to dissect the genetic architecture and evolution of a transcription factor’s specificity for DNA, using data from a combinatorial DMS of an ancient steroid hormone receptor’s capacity to activate transcription from two biologically relevant DNA elements. We show that the genetic architecture of DNA recognition consists of a dense set of main and pairwise effects that involve virtually every possible amino acid state in the protein-DNA interface, but higher-order epistasis plays only a tiny role. Pairwise interactions enlarge the set of functional sequences and are the primary determinants of specificity for different DNA elements. They also massively expand the number of opportunities for single-residue mutations to switch specificity from one DNA target to another. By bringing variants with different functions close together in sequence space, pairwise epistasis therefore facilitates rather than constrains the evolution of new functions.

    1. Evolutionary Biology
    Raphael Aguillon, Mieka Rinsky ... Oren Levy
    Research Article

    The circadian clock enables anticipation of the day/night cycle in animals ranging from cnidarians to mammals. Circadian rhythms are generated through a transcription-translation feedback loop (TTFL or pacemaker) with CLOCK as a conserved positive factor in animals. However, CLOCK’s functional evolutionary origin and mechanism of action in basal animals are unknown. In the cnidarian Nematostella vectensis, pacemaker gene transcript levels, including NvClk (the Clock ortholog), appear arrhythmic under constant darkness, questioning the role of NvCLK. Utilizing CRISPR/Cas9, we generated a NvClk allele mutant (NvClkΔ), revealing circadian behavior loss under constant dark (DD) or light (LL), while maintaining a 24 hr rhythm under light-dark condition (LD). Transcriptomics analysis revealed distinct rhythmic genes in wild-type (WT) polypsunder LD compared to DD conditions. In LD, NvClkΔ/Δ polyps exhibited comparable numbers of rhythmic genes, but were reduced in DD. Furthermore, under LD, the NvClkΔ/Δ polyps showed alterations in temporal pacemaker gene expression, impacting their potential interactions. Additionally, differential expression of non-rhythmic genes associated with cell division and neuronal differentiation was observed. These findings revealed that a light-responsive pathway can partially compensate for circadian clock disruption, and that the Clock gene has evolved in cnidarians to synchronize rhythmic physiology and behavior with the diel rhythm of the earth’s biosphere.