1. Evolutionary Biology

Late-life fitness gains and reproductive death in Cardiocondyla obscurior ants

  1. Luisa Maria Jaimes-Nino  Is a corresponding author
  2. Jürgen Heinze
  3. Jan Oettler  Is a corresponding author
  1. Zoologie/Evolutionsbiologie, Universität Regensburg, Germany
https://doi.org/10.7554/eLife.74695
Share this article
Cite this article
  1. Luisa Maria Jaimes-Nino
  2. Jürgen Heinze
  3. Jan Oettler
(2022)
Late-life fitness gains and reproductive death in Cardiocondyla obscurior ants
eLife 11:e74695.
https://doi.org/10.7554/eLife.74695
  2. Download BibTeX
  3. Download .RIS

Abstract

A key hypothesis for the occurrence of senescence is the decrease in selection strength due to the decrease in the proportion of newborns from parents attaining an advanced age – the so-called selection shadow. Strikingly, queens of social insects have long lifespans and reproductive senescence seems to be negligible. By lifelong tracking of 99 Cardiocondyla obscurior (Formicidae: Myrmicinae) ant colonies, we find that queens shift to the production of sexuals in late life regardless of their absolute lifespan or the number of workers present. Furthermore, RNAseq analyses of old queens past their peak of reproductive performance showed the development of massive pathology while queens were still fertile, leading to rapid death. We conclude that the evolution of superorganismality is accompanied by ‘continuusparity,’ a life history strategy that is distinct from other iteroparous and semelparous strategies across the tree of life, in that it combines continuous reproduction with a fitness peak late in life.

Editor's evaluation

The article will contribute significantly to our understanding of superorganism development and longevity in ants and will provide testable hypotheses in other ant species and other organisms.

https://doi.org/10.7554/eLife.74695.sa0

Introduction

The phenomenon that social insect queens live exceptionally long compared to solitary insects is widely recognized (Keller and Genoud, 1997; Carey, 2001). Given how prominent this is, however, only patchy information exists about the proximate mechanisms that are involved with the regulation of senescence, that is, a phase marked by an increase in relative mortality and a decrease in relative fecundity with age. Even less is known about the ultimate causes of social insect aging.

The classic trade-off between reproduction and maintenance shapes the life history strategy of species along a continuum between iteroparity (repeated events of reproduction) and semelparity (single event of reproduction) (Hughes, 2017). These strategies shape the way species age, that is, how resources are allocated to maximize fitness. In iteroparous species, fitness decreases after the first reproductive peak(s). Thus, the strength of selection against age-specific mortality decreases with age as the proportion of offspring that come from parents surviving to a specific age becomes smaller with time (Hamilton, 1966; Moorad et al., 2020). This is known as the selection shadow, which begins with maturity (Williams, 1957) and may negatively affect reproductive performance and survival (i.e., senescence). Classic model systems in aging research, such as Drosophila, Caenorhabditis (but see the discussion of quasi-semelparous hermaphrodites, Gems et al., 2021), mice, and humans, are of the iteroparous type, and a plethora of studies have revealed common mechanisms associated with senescence rate (Gems and Partridge, 2013). One prominent evolutionary theory of aging explains senescence by genes with antagonistic pleiotropic effects early and late in life (Williams, 1957). Semelparity instead predicts that organisms optimize their resources to one fitness peak, after which reproductive death occurs, that is, allocation of remaining resources into fecundity and not into maintenance. Thus, selection acts strongly against senescence before the single reproductive event. To understand how investment in reproduction of ant queens changes with chronological age and how social insect queens age, it is vital to investigate where they sit on this parity continuum.

Two dimensions are necessary to understand aging from an evolutionary perspective: the pace and the shape of demographic trajectories (Baudisch, 2011; Baudisch et al., 2019). The pace refers to factors that describe the timescale (e.g., life expectancy), and the shape refers to time-standardized measures of the distribution of mortality and fertility across a life history. Studies that capture the shape of aging (age-specific reproduction and mortality) of ants are scarce (reviewed in Cole, 2009), and based on punctual periods of growth and death of colonies, mostly of long-lived species in which colonies have a single queen (monogyny). Often such field data correspond to less than 20% of the estimated lifespan of the species (Atta cephalotes, Perfecto and Vandermeer, 1993; Atta colombica, Wirth et al., 2003; Pogonomyrmex owyheei, Porter and Jorgensen, 1988; Pogonomyrmex occidentalis, Keeler, 1993; Pogonomyrmex badius, Tschinkel, 2017). These studies have generally failed to capture the end of the queen’s lifespan and thus did not document senescence and lifetime reproductive investment. More complete yearly census data from Pogonomyrmex barbatus showed no relation between reproductive success (number of successfully established offspring colonies) and age (Ingram et al., 2013), but an increase in the production of male and female sexuals with age (Wagner and Gordon, 1999). In contrast, a study on a related species, P. occidentalis, showed no correlation between sexual production and colony size (as a proxy for age) once colonies had initiated sexual reproduction (Cole and Wiernasz, 2000), suggesting that it is difficult to infer the dynamics of age, colony growth, and reproduction from field data. To better understand how aging, senescence, and reproductive investment are related in ants, complete lifetime production data of individual queens are needed.

To study aging patterns and senescence of social insect queens, it is helpful to consider the colony as a superorganism (Wheeler, 1911; Boomsma and Gawne, 2018), analogous to a soma- (i.e., workers) and a germline (i.e., queens), where the investment into both castes is related and affects overall fitness (Bourke, 2007; Kramer and Schaible, 2013). Ant species such as Cardiocondyla obscurior, in which workers are completely sterile and seemingly without any direct reproductive power, exhibit an extreme case of superorganismality. By manipulating colony size, we expected to find trade-offs between lifespan and investment in queen/worker/male offspring. We monitored the lifetime production of individual queens in 99 single-queen colonies maintained with 10, 20, or 30 workers each (Figure 1—figure supplement 1A and B). Worker number corresponds to the colony size variation observed in the field (Schrader et al., 2014, Figure 1—figure supplement 2) and was standardized weekly. Queens whose egg production declined below a rate of ~10 eggs/week exhibited lethargic behavior, were less mobile, left the nest and/or were harassed by workers, and died within a few days to weeks. To assess if senescence was restricted particularly to the end of life, we compared RNAseq data of 18 of such prope mortem (Lat. near death) queens (between 28 and 49 weeks old) and 18 middle-aged queens (between 19 and 21 weeks), which were in their peak of fertility (Figure 1—figure supplement 3A and B). To compare queen and worker mortality, we tracked the survival of 40 workers kept in 40 colonies with either 10 or 20 nestmate workers.

Results

Reproductive strategy

The treatment (varying worker number) did not affect total production of eggs (package ‘generalized linear mixed models using template finder’ v. 1.1.2.3 in R) (Figure 1A, 10 vs. 20 workers: glmmTMB z-value = –0.38, p=0.70 and 10 vs. 30: z-value = –0.96, p=0.34) or worker pupae (Figure 1B, 10 vs. 20 workers: glmmTMB z-value = 0.09, p=0.93 and 10 vs. 30: z-value = –0.39, p=0.70). The treatment also did not affect the lifespan of queens (Figure 1—figure supplement 6A, Cox proportional hazard regression model, likelihood ratio test, X2 = 1.57, p=0.46), which was highly variable across treatments (variation coefficient: 32.2%, Figure 1—figure supplement 6B).

Figure 1 with 8 supplements see all
Download asset Open asset
Productivity of C. obscurior colonies across treatments.

(A) Total number of eggs, (B) worker pupae, and (C) queen pupae (N = 31, 34, and 34 for 10, 20, and 30 worker colonies, respectively). Significant differences are given with **p<0.01 and ***p<0.001. Boxplots depict upper and lower quartile plus 1.5 interquartile range (IQR).

We hypothesized that colonies that experienced a worker shortage would compensate by investing less into the production of new queens as these are larger and therefore more costly to produce. Indeed, queens with 10 workers (n = 31) produced significantly fewer queen pupae than queens with 20 (n = 34) (glmmTMB z-value = 2.81, effect size=1.97, p=0.005) and 30 workers (n = 34) (glmmTMB z-value = 2.58, effect size=1.78, p=0.009, Figure 1C) with no significant differences between 20 and 30 workers (glmmTMB z-value = –0.49, p=0.877). Similar results were obtained when accounting for the differences in biomass between workers and queens (Figure 1—figure supplement 7, Supplementary file 1A). Probably due to difficulties assessing precise egg numbers which are reared in piles, and extremely worker-biased caste ratios (average pupae developed into workers = 0.86), egg counts do not reflect these subtle but significant differences. The median sex ratio (queen/queen + male pupae) across treatments was 0.85 (25 and 75% quantiles = 0.79 and 0.90), and total production of male pupae (two types of males occur in C. obscurior: winged and wingless) was unaffected by the treatment (10 vs. 20: glmmTMB z-value = 1.94, p=0.05 and 10 vs. 30: glmmTMB z-value = 1.52, p=0.13, Figure 1—figure supplement 8). Queens produced very low numbers of winged males during their lifetime (mean = 0.36, median = 0, N = 99).

A first peak in the investment in queen pupae occurred around 15 weeks after the colonies were established (Figure 2A), followed by an increasing queen bias with age (Figure 2B). In general, new queens, which start a new colony, invest first in growing numbers of workers (ergonomic phase) and subsequently in the production of new sexuals, when the colony has reached the threshold required to enter the reproductive phase (Macevicz and Oster, 1976; Oster and Wilson, 1978; Beekman et al., 1998). This shift in caste ratio does not result from a drop in the production of pupae at the end of life. In contrast, pupa production is at its highest just before death (Figure 2—figure supplement 1). Importantly, in C. obscurior this caste ratio shift appeared to be a fixed trait, independent of colony size and queen lifespan. Both queens with short and long lifespans (below and above the mean lifespan of 25 weeks, Figure 2C and D, respectively), equivalent to queens with low and high productivity, exhibited late-life investment into queens.

Figure 2 with 2 supplements see all
Download asset Open asset
Lifetime investment.

(A) Numbers of worker and queen pupae produced over time, (B) queen/(queen + worker) pupae caste ratio produced by queens (n = 31, 34, and 34 for 10, 20, and 30 worker colonies, respectively), (C) pupae caste ratio for queens with lifespan below (n = 44), and (D) above the mean lifespan of 25 weeks, indicated by the dashed line (n = 55). After the queen’s death, eggs and larvae were allowed to develop into pupae for a final count. Therefore, smooth splines extend ca. 4 weeks after queen death.

In addition to the effect on caste ratio, the treatment had an effect at the colony level. We explored whether the quality of workers was affected by measuring the head width of workers produced over months 3–6 of the queen’s lifetime (approximately five workers per month). Head width of workers was 2% and 3% significantly smaller in colonies with 10 workers than in colonies with 20 (glmmTMB z-value = 2.22, p=0.026) and 30 workers, respectively (glmmTMB z-value = 2.68, p=0.007, Figure 3), but not different between colonies with 20 and 30 workers (glmmTMB z-value = 0.22, p=0.97). This suggests that small colonies lack sufficient numbers of nurse or forager workers, and indeed colonies collected in the field have worker numbers closer to 20 or 30 (Schrader et al., 2014; Figure 1—figure supplement 2).

Figure 3
Download asset Open asset
Worker quality across treatments.

Head width measurements of workers produced by queens of colonies with 10, 20, and 30 workers (n = 160, 112, and 68, respectively). Significant differences are given with * p<0.05 and **p<0.01. Boxplots depict upper and lower quartile plus 1.5 interquartile range (IQR).

After mean-standardizing queen age-specific mortality and fecundity (following Jones et al., 2014), we found that relative fecundity reached its maximum after ~16 weeks, before completion of the median lifespan (~25 weeks), and then decreased (Figure 4). Production of workers tightly followed the curve of egg production. Importantly, relative investment in queen and male pupae reached its maximum late in life (~28 weeks). This pattern is not due to the delay in development from egg to pupa because queen and male development only lasts ~5 and ~ 3 weeks, respectively (Schrempf and Heinze, 2006). Furthermore, C. obscurior ant queens exhibited a below-average level of adult mortality until week 30, after which mortality increased above the average level (Figure 4, Figure 4—figure supplement 1). This indicates maintenance of selection until after the peak of relative investment in sexual offspring. Therefore, queens continue to experience strong selection even at high ages, that is, weeks after they reached the mean lifespan. Monitored workers in colonies with 10 or 20 nestmates did not differ in survival (Cox proportional hazard regression model, likelihood ratio test, X2 = 0.06, p=0.8). Therefore, the mean-standardized age-specific mortality was calculated for the 40 workers. Note that regardless of the differences in timescale, the shape of mean-standardized mortality of workers was similar to that of queens (Figure 4—figure supplement 2). This suggests that aging is a genetically fixed trait expressed by queens and workers alike.

Figure 4 with 2 supplements see all
Download asset Open asset
Relative mortality and fecundity as a function of age.

Mean standardization of age by dividing age-specific mortality and fecundity of queens (n = 99) by their means after maturation (Jones et al., 2014). Survivorship (black dashed line) is depicted on a log scale. The graph uses a Loess smoothing method (span = 0.75) and a confidence interval of 95%. The dashed gray line at y = 1 indicates when relative mortality and fertility are equivalent to mean mortality and fertility.

RNAseq of prope mortem queens

To determine if queens show signs of reproductive senescence and loss of physiological function, we analyzed gene expression data of prope mortem queens exhibiting decreasing egg-laying rates and middle-aged queens that were at their peak reproductive performance. To account for possible effects of fertility, we sampled queens with low, medium, and high egg productivity at 18 weeks of age (Figure 1—figure supplement 4A and B). We subjected the head plus thorax and the gaster (see methods for terminology, Figure 1—figure supplement 5) to RNAseq separately to assess if reproductive tissue shows a different physiological wear and tear than head-thorax tissue. The analyses revealed that head-thorax and gaster tissues showed similar mapping rates to the genome (Figure 5—figure supplement 1A and B), but that gaster samples had a lower GC content on average and more duplicated reads (Figure 5—figure supplement 1C–F) in prope mortem queens compared to middle-aged queens.

Of the 20,006 expressed genes in head-thorax tissue, 3565 (17.8%) genes were differentially expressed between middle-aged and prope mortem queens (after false discovery rate [FDR] adjustment p<0.001, DESeq2, Source data 9). Of these, 1725 genes (48%) were upregulated and 1840 genes (52%) were downregulated in prope mortem queens compared to middle-aged ones. Gene Ontology (GO) term enrichment revealed signs of rapid physiological decay of prope mortem queens, such as reduced translation, proteasomal, ribosomal, and mitochondrial function (Fisher test using the weight01 algorithm, p<0.05, Supplementary file 1B, Figure 5), increased splicing, and transcript processing (Supplementary file 1C, Figure 5—figure supplement 2). Such processes have previously been related to aging in several model organisms (López-Otín et al., 2013); for example, the loss of protein homeostasis (Hipp et al., 2019), the decrease in ribosomal proteins (Walther et al., 2015), alterations in the mitochondrial function (Green et al., 2011), disruption of splicing (Bhadra et al., 2020), and others (Harries et al., 2011). Another characteristic of aging, changes in gene connectivity among gene expression networks found in mice (Southworth et al., 2009), was not affected by age in C. obscurior (calculated using the softConnectivity and the biweight midcorrelation functions on gene networks for middle and prope mortem queens using WGCNA, and modeled using glmmTMB, z-value = −1.7, p=0.09). Principal component analysis (PCA) ordination of the head-thorax tissue separated middle-aged and prope mortem queens by age (PERMANOVA test, F-value = 7.59, p<0.001), but not by fertility (F-value = 1.09, p=0.26) or duplication percentage (Figure 5—figure supplement 3A and B).

Figure 5 with 7 supplements see all
Download asset Open asset
Enriched Gene Ontology (GO) terms downregulated in prope mortem queens compared to middle-aged queens in the head-thorax tissue.

Functional annotation and enrichment analysis using topGO (version 2.46.0) and the weight01 algorithm to calculate significance for (A) biological processes, (B) cellular components, and (C) molecular functions.

In the gaster tissue, 4832 (24.3%) of 19,925 expressed genes were differentially expressed between age groups (after FDR adjustment p<0.001, DESeq2, Source data 10). Of these, 2306 genes were upregulated (48%) and 2526 downregulated (52%) in prope mortem queens compared to middle-aged queens. GO term enrichment likewise showed that many fundamental processes were affected in prope mortem queens, such as DNA damage, telomere maintenance, and enrichment of transcription processes (Supplementary file 1D, Figure 5—figure supplement 4), and among others processes related to protein processing, glycolytic processes, and the Notch signaling pathway were downregulated (Supplementary file 1E, Figure 5—figure supplement 5). In contrast to head-thorax tissue, gene connectivity among gene co-expression networks in the gaster was significantly different (prope mortem queens: median = 69.45; middle-aged queens: median = 52.56) (glmmTMB, z-value = −19.5, p<2e-16), but contrary to what was found for aged mice (Southworth et al., 2009).

The PCA of the 500 most variable expressed genes in gaster tissue shows that the samples group according to age (PERMANOVA test, F-value = 13.91, p<0.001), fertility level (F-value = 1.95, p=0.04), and also the percentage of duplicated reads in the libraries (F-value = 4.83, p=0.002, Figure 5—figure supplement 6A and B). This is not a typical technical artifact (no correlation to sequencing lane, RNA concentration, or quality). Spike-in reads were used as a control for library preparation and showed a positive linear relationship between expected and observed reads independent of age group, tissues, and lanes (Figure 5—figure supplement 7A–C). However, this linear relationship has different slopes among age groups in the gaster samples (Figure 5—figure supplement 7B), indicating biological changes with age pertaining specifically to the gaster.

In spite of these discrepancies between tissue types, 104 GO terms were significantly enriched in both tissues, of these 44 in prope mortem queens (Supplementary file 1F) and 60 in middle-aged queens (Supplementary file 1G). Thus, signs of similar physiological pathologies occur in reproductive and non-reproductive tissue.

Discussion

The near-ubiquitous occurrence of senescence has been explained by two classic prevailing evolutionary theories, mutation accumulation and antagonistic pleiotropy (Medawar, 1941; Williams, 1957). These theories have in common the basic assumption of the existence of a ‘selection shadow’: a decrease in the force of natural selection with age. The selection shadow leads to loss of function and senescence, that is, an increase in relative mortality and a decrease in relative fecundity with age (Maklakov and Chapman, 2019). Originally explained as a consequence of extrinsic mortality, models have shown that the strength of selection is in fact influenced by the proportion of offspring coming from parents that survived to a certain age (Hamilton, 1966; Moorad et al., 2020). Extensive demographic data show a huge diversity of aging patterns across metazoan species, ranging from 20 times the average mortality at terminal age to less than a half in other species (Jones et al., 2014; Cohen, 2018). In some cases, a short phase of senescence is self-evident, for example, in semelparous species such as salmon, where death follows reproduction to provision the next generation with resources.

Here, we show for the first time the shape of aging in a social insect. While fecundity decreases in the ant C. obscurior, reflecting reproductive senescence, investment into sexuals reaches a maximum late in life, regardless of individual fitness (queen lifespan or total egg productivity) and colony size. Males in this species usually transfer an excess amount of sperm (Schrempf and Heinze, 2008), and only one queen showed signs of sperm depletion and produced only males at the end of her life. Therefore, reproductive senescence cannot be explained by sperm depletion. The magnitude of the investment (i.e., number of queen pupae produced) is affected by the number of workers available. In C. obscurior, most new queens were produced by queens older than the mean queen lifespan, indicating that queens continue to experience strong selection at high ages. This is in line with the hypothesis that the strength of selection against age-specific mortality is proportional to the probability for any offspring in the population to be produced by parents of that age and older.

Strikingly, relative mortality did not increase directly after maturity or after total egg production started decreasing, but after the production of sexual pupae had reached its maximum. IAccordingly, transcriptome profiles of prope mortem queens shortly after this investment peak, which produced decreasing numbers of eggs, revealed signs of a broad range of physiological pathologies. The changes seemed stronger in the gaster (e.g., percentage of duplicated reads), which contains the reproductive organs and most of the digestive system, but to a similar extent occurred in head and thorax, containing most neuronal and muscle tissue. Such a systemic breakdown is expected assuming that the entire physiology is optimized towards a fitness peak. Strikingly, a comparative transcriptomic study of 8-week young queens with fully mature C. obscurior queens at or close to their peak fecundity (18 weeks old) did not find signs indicative of aging, but in comparison to aged Drosophila flies an opposite regulation of processes (e.g., cellular ketone, carbohydrate, and organic acid metabolic processes) and genes (e.g., ref(2)P, emp, P5cr-2, CCHa2, NLaz, Sirt6) involved in aging (von Wyschetzki et al., 2015). Furthermore, a gene co-expression network study using the same data showed higher connectivity in middle-aged queens, indicating increased transcriptional regulation with age (Harrison et al., 2021). Together, this suggests that the physiology of queens is maintained until the fitness peak is reached, at which time they undergo physiological deterioration, while still being reproductively active. This pattern is reminiscent of semelparous species with reproductive death rather than that typical of iteroparous species in which selection against age-specific mortality decreases after a first reproduction event and actuarial senescence unfolds under the selection shadow.

Conclusion

Superorganismality is a major evolutionary transition, and this transition is accompanied by a change in the mode of reproduction. We propose that the evolution of ‘continuusparity’ (Lat.: ‘continuus’ meaning incessant/successive; and ‘parere’ meaning giving birth), that is, the combination of lifelong continuous reproduction and increasing fitness returns late in life, underlies the delay of the selection shadow, the maintenance of selection strength against age-specific mortality, a brief phase of senescence late in life, and finally reproductive death. This is not to be confused with the meaning of the term negligible senescence as actuarial and reproductive senescence clearly occur at the end of life.

‘Continuusparity’ emerges as a combination of iteroparous and semelparous characteristics: reproduction resembling continuous iteroparous species but without the inter-parous nonreproductive breaks, during which nests are built, mating occurs, resources are acquired, etc. The iteroparous solitary ancestor of ants is thought to be related to mud dauber wasps (Sphecidae) and cockroach wasps (Ampulicidae) (Ward, 2014), parasitoid wasps with mass provisioning. This combines with aging/resource allocation patterns of semelparous species (which are in contrast mostly short-lived), optimized towards one reproductive episode at the end of life, followed by reproductive death. Continuous reproduction is possible because no extra time or energy is necessary for further acquisition of resources, brood care, territoriality, etc., because of the extended phenotype, the colony. With time the colony increases in size, and resources increase accordingly, analogous to solitary organisms with lifelong growth (Keller, 1998). Continuous reproduction is further facilitated by the presence of a spermatheca in female insects, which allows for a single mating event and lifelong sperm storage. Thus, the costs of additional matings are zero.

We propose that continuusparity and its effect on the shape of queen aging is a property of superorganismality, and that this life history strategy ultimately underlies the evolution of long lifespans in social insects. For the pace of aging, it is important whether queen and worker interests are aligned, and whether direct and indirect reproductive investments of queens and workers are optimized. With respect to which proximate mechanisms regulate aging in social insects, this framework predicts that there are no genes/pathways with antagonistic pleiotropic effects because there is no ‘later in life’ past the fitness peak. Under this perspective, many questions remain open. Queens appear to have different properties that influence both lifespan and fertility (Kramer et al., 2015), probably determined during larval development (Schultner et al., 2017), and which are key to understanding why some individuals live longer than others. What are these properties, and how are they determined, maintained, and finally terminated? Queens do not senesce until shortly before death, so how is the trade-off between somatic maintenance and reproduction resolved? Using this framework, we can now start to study the proximate regulators that maintain the homeostasis in C. obscurior ant queens, which remain hidden in the excess number of associated processes.

Materials and methods

The species

Request a detailed protocol

C. obscurior is probably the best studied ant species with respect to aging due to the relatively short lifespan of queens (~6 months). Colonies comprise a few queens (body length ~3 mm), a few dozen workers (~2 mm), and nest in small cavities in dead twigs, aborted fruits, rolled leaves, under bark, etc., in trees and shrubs (Schrader et al., 2014). Virgin queens usually mate once with related wingless males inside the natal colony (Heinze and Hölldobler, 1993; Schmidt et al., 2016; Heinze and Hölldobler, 2019), generally stay in the nest, and new colonies are formed by budding of colony fragments. This mode of reproduction from small propagules allows for successful colonization of disturbed habitat in warm climates around the world (Heinze and Delabie, 2005; Heinze, 2017). Various social, environmental, and biotic factors affect the lifespan of queens (Oettler and Schrempf, 2016). Queens that lay more eggs (total output and weekly rate) live longer than less fecund queens, irrespective of body size (Kramer et al., 2015), and thus seem to evade the common trade-off between reproduction and maintenance.

Reproductive strategy

Request a detailed protocol

We set up 138 freshly eclosed queens from stock colonies of a Japanese population (OypB, from the Oonoyama Park in Naha, Okinawa) established in the laboratory since 2011. The experiment took place between January 2019 and January 2020. Queens were allowed to mate with a single wingless male and were placed in nest boxes with either 10, 20, or 30 workers from the maternal colony to establish monogynous colonies, for lifetime production and RNAseq (n = 46 each). These numbers of workers represent the naturally occurring number in the field and correspond to the first, median, and third quantile of number of workers in this population (n = 62, median = 28.5, Figure 1—figure supplement 2). The colony was set up with half of the workers selected from inside of the nest near the brood (younger nurses) of the stock colonies, and the other half from outside the nest (older foragers) in order to minimize a putative effect of worker age on the queen (Giehr et al., 2017). Colonies were kept under a 12 hr dark 22°C/12 hr light 26°C cycle and fed ad libitum three times per week with diluted honey (0.6:1 honey: distilled water), cockroaches, and flies. Once per week workers, eggs, and all pupae (worker, queen, winged, and wingless male) were counted and queen survival was monitored. Additionally, the number of workers was standardized to the assigned treatment, and newly produced sexual pupae produced were removed. C. obscurior workers are sterile, and all produced offspring originated from the focal queen. The number of counted eggs correlates with the production of workers, queens, and the workers and queens together (Figure 2—figure supplement 2A–C, Kendall’s rank correlation test, p<0.001: eggs-worker pupae, τ = 0.70; eggs-queen pupae, τ = 0.59; eggs-worker and queen pupae, τ = 0.73). Pupae might have been counted more precisely than eggs, especially when larger numbers of eggs were produced. Pupae are hardly missed compared to eggs that tend to cluster together. Eggs and worker pupae might have been counted more than once as development lasts a median of 8 and 18 days for eggs and worker pupae, respectively. Colonies were counted ca. 4 weeks after the queen’s death, until the last eggs had developed into pupae. Finally, three colonies (10 worker treatment) were not considered in the analysis as they were accidentally killed, leaving a total of 99 colonies for lifelong tracking and 36 colonies for RNAseq analysis.

Worker aging

Request a detailed protocol

To examine worker aging, 40 focal unmarked worker pupae were set up in individual colonies with 10 or 20 marked workers (n = 20 each). These two treatments were selected because colonies with 20 and 30 workers did not differ in queen productivity. Marking of nonfocal workers was done by clipping the tarsus of the middle right leg. Colonies were set up with brood (5 larvae in the 10 workers colonies, and 10 in the 20 workers colonies), and two wingless adult queens to avoid a queenless period in case one died. The survival of the focal worker was monitored, and the number of marked workers, queens, and larvae was standardized weekly to the assigned treatment. Newly produced pupae were removed. Dead marked workers were replaced with fresh worker pupae, which were marked 1 or 2 days after eclosion to avoid confusion with the non-clipped focal worker. Dead queens were replaced with adult ones.

Offspring investment

Request a detailed protocol

360 freshly eclosed adult workers were sampled monthly for head width measurements (from the third to sixth month of the queen’s life, and up to five workers depending on availability). Workers were dried, pinned, and blindly measured using a Keyence Microscope at 200×. A single worker was chosen randomly and measured 10 times to obtain a proxy for measurement error (mean = 383.61 μm, standard deviation = 5.05 μm).

Statistical tests

Request a detailed protocol

To test for significant differences between treatments, we used generalized linear mixed effects models within the R package glmmTMB (R version 3.5.2, Pinheiro et al., 2011) and a negative binomial distribution for count data. If the count data and caste investment ratios were log transformed, a Gaussian family distribution was used. The dependent variable was analyzed as a function of the fixed effects: treatment (number of workers as a factor); and random effects: stock nest and box of origin, box of set up, setup date. All models were also graphically checked for consistency and model diagnostics were performed using the DHARMA package (R version 0.3.3.0, Hartig, 2020). Caste ratio was calculated as queen over the total caste investment (as Queen*c / [Queen*c + Worker]). The coefficient or correction factor c is used as the dry average weight measurements of queen over workers to the power conversion factor of 0.7 assuming differences in metabolic rates between queens and workers adopting the logic for sex ratio investment (Boomsma, 1989). As this is an assumption, we used different values of c. The results are robust to power conversion values of 0.6–1 (Supplementary file 1A). To test for differences in head width, we used the average of the head width measurements of the workers per time point (each month). Predictions of the data were visualized using the loess method with the geom_smooth function and default span (ggplot2 v.3.3.2). Relative mortality and fecundity as a function of age were mean-standardized by dividing age-specific mortality and fecundity by their means after maturation, following Jones et al., 2014. In contrast to Jones et al., 2014, the whole life range was considered until death, since removing the last 5% of survivorship showed similar results. Age-specific mortality without the mean standardization was also estimated for the 135 queens (99 and 36 ant queens for RNAseq) using a survival Bayesian trajectory analysis (Figure 4—figure supplement 1). Data is available as Source data 1–Source data 7, and the R-script used as Source code 1.

Prope mortem queens selection

Request a detailed protocol

To obtain samples of low, medium, and highly productive queens, 18 queens at age 19–21 weeks were sacrificed for RNAseq based on egg productivity until week 18. Values of weekly egg productivity below the first quantile for the treatment group (colony size) were considered as low, values between the first and the third quantile as medium, and values greater than the third quantile as high. An additional 18 queens were monitored until they showed decreasing fertility (Figure 1—figure supplement 3) and one or more of the following signs of senescence: lethargy, loss of mobility, presence outside the nest, and/or harassment by workers. These senescent queens were also selected based on low, medium, and high fertility, and then sacrificed (28–49 weeks old). Queens were snap-frozen in liquid nitrogen after the head and thorax was separated from the gaster with a blade between the petiolus and post-petiolus in a drop of PBT 0.3% (phosphate buffered saline and Tween 20). During this procedure, queens were manipulated for less than 1 min.

Terminology

Request a detailed protocol

What we refer to as ‘thorax’ actually refers to the thorax plus the fused first abdominal segment, together making up the ‘mesosoma’ in the Hymenoptera. The ‘metasoma’ in Hymenoptera comprises the segments making up the constriction plus the hind end. In the ant subfamily Myrmicinae, this constriction is made of two segments: the petiole corresponds to the second, constricted, abdominal segment, while the post-petiole refers to the third, constricted, abdominal segment. The ‘gaster’ refers to the bulbous posterior part (Figure 1—figure supplement 5).

RNAseq

Request a detailed protocol

Total RNA was extracted using the ReliaPrep kit (Promega) from the 72 samples (36 queens, two samples per queen: head-thorax and gaster). Spike-In RNA Variant Controls (SIRV-Set 3 Lexogen #05101‚ Lot 5746/001492) were spiked to a 2% fraction of the total RNA (measured using Bioanalyzer – Agilent Technologies). 8 of the 72 samples showed RIN values below 7 (gaster samples from older queens that seemed more degraded). For those samples, the concentration of RNA was estimated based on the mean value of the nondegraded gaster samples. Total RNA was amplified using single primer isothermal amplification (SPIA , Ovation RNA-seq System V2, Tecan) prior to cDNA generation. The library preparation and sequencing (100 bp PE) was performed at the Cologne Center for Genomics using Nextera XT sequencing on a NovaSeq 6000 platform. Reads were trimmed with fastp v.0.20.1 to a minimum length 70 and from Nextera adapters. Then, SortMeRNA version 4.2.0 was used to discard undesired rRNA reads using the default database (smr_v4.3_default_db. fasta). Remaining reads were aligned using hisat2 (version 2.1.0) to the newest version of the genome (Cobs.2.1., Errbii et al., 2021). Putative splice sites were obtained using gffread (version 0.12.1), and the extensive transposable elements annotation v.2.1 (Errbii et al., 2021) was considered for the mapping procedure. Samtools (version 1.9) was used to sort and convert.sam into .bam files.

Gene expression analysis

Request a detailed protocol

After filtering genes with 0 values, we used a gene set of 20,006 genes for the head-thorax analysis and 19,925 genes for the gaster. PCA plots were produced to visualize the samples after variance stabilizing transformation. An analysis of the homogeneity of group dispersions (variances) was performed (multivariate analog of Levene’s test for homogeneity of variances), with the function permutest and 999 permutations (vegan package v. 2.5–7) to test for differences in variance among the age groups (middle-aged and prope mortem queens) (betadisper, vegan package).

Subsequently, a nonparametric multivariate ANOVA (PERMANOVA) test was performed (999 permutations) with the design model Head-Thorax_expression ~ AgeGroup + Fertility, with two (middle-aged, prope mortem) and three levels (low, medium, and high fertility) to test for statistical differences in the transcriptomic profiles due to age group (senescent or not) and level of fertility (low, medium, and high) using the adonis function (vegan package), with the default Bray distance method. Age and fertility (average number of laid eggs per week) were scaled and centered. Then, differential expression was analyzed using the design = ~ Eggs per week + Age group, with the R package DESeq2 (v. 1.28.1). Age group was used as categorical variable with two levels (middle-aged and prope mortem), and log2FC were calculated as log2 [prope mortem/middle-aged].

The cutoff threshold of statistical significance (alpha parameter) was set as 0.001 after p-value adjustment with FDR. Functional annotation and enrichment analysis was done using topGO (version 2.46.0) and the weight01 algorithm implemented in the package.

A signed weighted co-expression network was constructed using the WGCNA package (v. 1.70–3), and the count data transformed using variance-stabilizing transformation from the DESeq2 package after excluding genes with 0 read values. For the head-thorax tissue network, the total set of 20,006 genes was used, and the WGCNA was performed with default parameters and a soft threshold power 14, based on the scale-free fit index as recommended in the manual. We compared the connectivity of the two separate networks, one for middle-aged and one for prope mortem queens, with the softConnectivity and the biweight midcorrelation functions.

Data availability

Data generated or analysed during this study are included in the manuscript and supporting file; A source Data file has been provided for Figures 1–5 and accompanying figure supplements. Raw sequencing data have been deposited in SRA under the BioProject accession number PRJNA819887. The Bioproject is part of the NCBI So-Long Umbrella Bioproject: "The So-Long project: Sociality and the Reversal of the Fecundity/Longevity Trade-off.

The following data sets were generated
    1. Jaimes-Nino L
    (2022) NCBI BioProject
    ID PRJNA819887. Cardiocondyla obscurior Raw sequence reads (TaxID: 286306).
    https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA819887

References

  11. Book
    1. Cole BJ
    (2009)
    The Ecological Setting of Social Evolution: The Demography of Ant Populations
    Harvard University Press.
  39. Book
    1. Oster GF
    2. Wilson EO
    (1978)
    Caste and Ecology in the Social Insects
    New Jersey, USA: Princeton University Press.
  41. Book
    1. Pinheiro J
    2. Bates D
    3. Debroy S
    4. Sarkar D
    (2011)
    Linear and Nonlinear Mixed Effects Models
    New Prairie Press.
    1. Ward PS
    (2014) The phylogeny and Evolution of Ants
    Annual Review of Ecology, Evolution, and Systematics 45:23–43.
    https://doi.org/10.1146/annurev-ecolsys-120213-091824

Decision letter

  1. Ehab Abouheif
    Reviewing Editor; McGill University, Canada
  2. Claude Desplan
    Senior Editor; New York University, United States

Our editorial process produces two outputs: i) public reviews designed to be posted alongside the preprint for the benefit of readers; ii) feedback on the manuscript for the authors, including requests for revisions, shown below. We also include an acceptance summary that explains what the editors found interesting or important about the work.

Decision letter after peer review:

Thank you for submitting your article "Late fitness gains explain the delay of the selection shadow in ants" for consideration by eLife. Your article has been reviewed by 2 peer reviewers, and the evaluation has been overseen by a Reviewing Editor and Claude Desplan as the Senior Editor. The reviewers have opted to remain anonymous.

The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.

Essential revisions:

The editors and both reviewers agree that this dataset has potentially important implications for understanding superorganism development and the long lifespans of social insects. Therefore, it will be of interest to a broad audience of biologists. Both reviewers have substantive comments on data interpretation and writing, as well as requests for some additional information and statistical analyses. Specifically:

1) Data interpretation:

a) A particular concern regarding the interpretation of the current data set is that it does not necessarily provide an ultimate explanation for the long life of the queens: an alternative, more conservative, interpretation of the current data set is that it only provides documentation of negligible reproductive senescence. According to Reviewer 2, the long lifespans of queens is strongly correlated too, but is not necessarily caused by, the capacity to reproduce late in life. For example, the long lifespans of queens could be ultimately caused by a shifting of reproductive costs to workers, or by living in protected environments. Therefore, in Reviewer 2's view, the current data set can only address the question of why social insect queens lack reproductive senescence. In other words, why do social insect queens retain the capacity to reproduce late in life? In my view, you (the Authors) should be very careful to distinguish ultimate from proximate causes in this study and be very clear about the precisely question this current data set can or cannot address.

b) Following on point (a) above, Reviewer 3 raises another particular point of concern, which is the lack of consideration for the alternative hypothesis that the drop in fertility later in life is caused by sperm limitation and not reproductive senescence. As Reviewer 3 points out, unless you can rule out this alternative hypothesis, it must be explicitly considered in the main text and would change the interpretation of the data, but not necessarily the relevance of the study, especially if you are willing shift the focus instead on understanding intrinsic properties of superorganism development and less on ageing per se.

c) Another issue of particular concern is the unnecessary focus on queen control. The inflection point, where colonies produce more queen offspring late in life, appears to reflect an intrinsic property of colony (superorganism) development that is independent of the colony-size, lifespan or worker investment. How this inflection is mechanistically controlled is a question you should address in the Discussion and should offer several hypotheses. For example, three possible hypotheses are: (1) An intrinsic colony-level threshold that is be controlled by queens and workers, such that workers change their feeding/foraging behavior in response to a queen pheromone once the queen reaches a certain age; (2) Queen embryos are produced all the time, like in the genus Monomorium, but are killed by workers until the queen reaches a certain age and changes its pheromone status; (3) as you have suggested, the threshold could be genetically determined in the queen regardless of worker input. Therefore, you should remove the focus on the assumption of queen control and provide a more systems-level (superorganismal) explanation for the inflection for producing queen offspring independent of colony size in the discussion.

2) RNAseq data: If you feel that RNAseq data can strengthen your arguments in light of the comments / critiques on data interpretation raised above, then you should add these data to the manuscript. If you decide to do so, please make sure to explain why this is important both in the cover letter and in the response to reviewers.

3) Reviewers 2 and 3 require some missing statistics as well as analyses, which should be straight forward to perform.

4) The manuscript, as presently written, leaves many open questions that need clarification, therefore there is a need for major revision to respond to the queries of Reviewer 2 and 3 (see Reviewer comments below).

5) You should provide additional information about this species' life history (see Reviewer comments below).

Reviewer #1 (Recommendations for the authors):

This is an interesting study looking at the evolution of ageing in social insects using ants as a model. As I haven't seen the initial submission, I have looked at the manuscript and the response to reviewers and I base my suggestions on both documents.

Evolution of ageing remains only partially understood and this field seems to be experiencing a sort of renaissance in recent years with a surge of theoretical advances and new empirical findings. Queens of social insects, and ant queens in particular, have remarkable lifespans and understanding the biology of their long life can help in understanding the biology of ageing in a more general sense.

In this study, the authors focus on following quite a large number of ant (C. obscurior) colonies and provide intriguing data in relation to age-specific mortality and reproduction. The gist of their argument is that the mortality is decreasing with age while reproduction (production of sexuals) is increasing with age, such that there is little evidence of ageing in this species.

Overall, I think this is an interesting dataset that provides important information that will advance the field. However, I think the manuscript currently lacks clarity, structure and suffers from poor formulation of ideas in places, and is rather difficult to follow even for an expert in the field. I think that it requires quite a bit of work to sort this out; nevertheless, I think the dataset is, potentially, sufficiently strong to justify publication in a high-impact journal such as eLife following major revision. However, I also have a methodological question (#15) which could be key for the interpretation of the results

I will provide comments in a sequential order, as they appeared while reading the manuscript.

1. You must include the species name either in the title or in the abstract. Currently it reads as if you followed 102 colonies of random species of ants.

2. L 21 advanced not advance.

3. "Furthermore, mortality decreased with age in queens and workers, supporting the hypothesis that aging trajectories are adaptive." I didn't understand this sentence while reading the abstract and reading the reminder of the paper didn't make clearer. What is meant here? Why would decrease in mortality suggest that "aging trajectories" are adaptive? What are "aging trajectories"?

4. "We argue that selection for late life reproduction delays the occurrence of the selection shadow, leads to the apparent absence of senescence in ants and underlies the evolution of long lifespans." This is very confusing – if "selection shadow" is delayed, how come there is an apparent absence of senescence? "Delayed" suggests that senescence will occur but at a later point. Also, this sentence is a bit clumsy and contains grammatical errors.

5. "Together, this highlights the unique aging patterns of social insects in comparison to other iteroparous species across the tree of life." Why? Is there evidence that negligible senescence is unique characteristic of social insects? Also, "aging patterns" is unclear.

6. Introduction: I found it very confusing that the paper on ageing starts with the discussion of reproductive queens and workers of social insects. You need a general introduction paragraph. The most straightforward solution is to take the first paragraph from Conclusions (LL 138-151) and start your main text with this.

7. Introduction: "It is argued that the costs of reproduction (energy intake, brood care) are outsourced to the workers, but this is a necessary and not sufficient condition for the evolution of long lifespans in social insects and does not explain these aging patterns." Why? There is no explanation as to why this is not sufficient, and it is absolutely not clear to me – please explain.

8. LL 45-58 – I like this paragraph and I think it presents good arguments in favour of the novelty provided by this study

9. L 58 – "aging, senescence" – if you treat ageing and senescence as different terms, you have to explain why and provide definitions for each term.

10. L 58 – "investment" – what investment? Reproduction? Somatic maintenance? This is sentence is remarkably unclear.

11. I think that "superorganism" approach, while still debated in the literature is entirely justified in this system.

12. L81 – While not directly relevant as you did not find a correlation, how did you test for causality of correlation? Your experimental treatment arguably had no effect neither on egg production nor on queen lifespan?

13. "Generalizing this finding, we suggest that lateLife investment and its effect on queen aging is an intrinsic property of superorganismality irrespective of the mode of caste determination and social organization." – this is unclear, I thought there is no effect on queen's ageing?

14. I think you should analyse age-specific mortality by selecting the best-fit model in a specialized demography analyses package, such as, for example, BASTA.

15. My understanding is that queens live for 40-50 weeks max (Figure S3). Figure 4 suggests that from week 30 onwards the production of eggs, worker pupae and queen pupae declines. This suggests that while queen mortality declines in late life, so does queen reproduction. So, queens of this species do show reproductive senescence?

Yes, your data suggests that relative investment into reproduction (queen worker ratio) increases with age, but the absolute number of queens declines with age? To me, this suggests an interesting result from the life-history theory perspective – increased investment in reproduction with reduced residual reproductive value, but not necessarily the absence of reproductive senescence. Perhaps I didn't understand the results, please clarify.

16. L155 – this relative mortality is very unclear to me. I think you need proper mortality rate analyses in BASTA and estimate different age-specific mortality parameters.

17. L 156 'random' offspring?

18. "In turn, the short selection shadow might have led to the delay of a distinguishable senescence phase." I think I understand what you mean but this is very unclear, please re-write.

Reviewer #2 (Recommendations for the authors):

Line 42: I'd like clarification on your argument here. Why is this not sufficient in your opinion? I don't see any alternative explanation, and workforce dynamics is central to your experimental design. After all, you are manipulating workforce size to expose trade-offs in queen decisions, so how does it not ultimately come down to the workforce?

Line 72: Why not just say something simple like "compare queen and worker mortality". Throwing in a new term like aging trajectories makes it seem like you are doing something different.

Line 92: I would have liked to see more of your data on the male production, even if they were less common than queens, 15% of your reproductive output is not negligible. Does male production fit the same lateLife pattern as queen production? Do you see the same differences between treatments when you analyze combined alate production? I expect the 10-worker treatment may invest more in males, since they are usually cheaper to produce. I know this isn't the question you are asking but it may explain some of the drop in queen production in that treatment.

Line 154: I don't understand why you seem to downplay the clear influence that the workforce has on this dynamic. From my perspective, your data beautifully shows how these queens transition from worker production to alate production as they age, and that colonies with low worker numbers cannot invest as much into reproduction, which is expected in a species that must establish a robust colony before devoting resources to reproduction. If a queen survives to a relative old age, she will also have a large, established workforce that can raise a lot of alates, so of course we will see longevity selected for at that point (as indicated by your own data, if she reaches that age without a sufficient workforce, her alate production will suffer).

Another important omission in the conclusion is the inevitable problem of sperm limitation. You frame the system as generating a directional push towards negligible senescence but it is more likely stabilizing selection that produces a sufficient workforce and then devotes as much remaining sperm as possible towards queen production. Selection will ensure that an established queen lives long enough to exhausts her sperm storage, which would contribute to the short selection shadow you discuss. The increase in relative mortality at 30 weeks may just represent the optimal lifespan of the colony. In figure 2A we see the peak in queen production at that time point, with both worker and queen production falling dramatically by 40 weeks, potentially signaling sperm depletion. Upon closer inspection of S3, it looks like there were only a few colonies that survived past week 38, but then persisted for another couple months. I have to wonder how reliable the drop in relative mortality is after 40 weeks when it represents only a few colonies. Of course, this is my interpretation based off of the figures, so please include the raw numbers if I am mistaken.

Line 212: Is this species monogamous in the wild? If so, this should be stated because if they are normally polyandrous, having only one mate may introduce sperm limitation.

[Editors' note: further revisions were suggested prior to acceptance, as described below.]

Thank you for resubmitting your work entitled "LateLife fitness gains and reproductive death in Cardiocondyla obscurior ants" for further consideration by eLife. Your revised article has been evaluated by Claude Desplan (Senior Editor) and a Reviewing Editor. The paper is significantly improved but we believe that additional editorial changes must be made to increase the impact of the paper, as described below, before the paper can be accepted for publication.

Although the transcriptomic data were analyzed properly, the conclusions that can be drawn from these data are limited, primarily because the difference between gaster and head/thorax is not easy to interpret. However, we think that these data do add something to the paper by validating some of the conclusions drawn from indirect observations. Therefore, the RNAseq do not necessarily serve the purpose of identifying which mechanisms are failing, but instead provide evidence that the queens are experiencing delayed senescence, something that was not obvious in the previous version. Please keep the interpretation of RNAseq simple and revise the manuscript appropriately to make its purpose clearer (see comment below on Line 82-83 and 178 and 179 below).

Furthermore, the singular focus on the queen is problematic, especially if you are attempting to adopt a superorganism perspective. First, the evidence presented that only the queen is in control is indirect and weak. Therefore, we ask you to maintain a neutral stance on this issue and acknowledge that the workers may play an important role in this process. Second, your arguments do not depend on whether determination is in the maternal, embryonic, or larval stage and we ask you to remove this mention as it is a non-central point in your argument (lines 127-130). Given these two points, we ask you to put back the data on worker longevity (although it can be presented in the supplemental data), because in our view, without direct supporting evidence that it is all about the queen, workers remain an important part of the equation. We therefore ask you to revise lines 127-130 in the Intro and 287-297 in the Conclusions section.

Essential revisions:

Abstract:

Line 23-25: We ask you to add something to describe the inflection point as a threshold, an inflection point that changes in an abrupt, non-linear-like way which distinguishes it from a gradual linear-like increase.

What does "reproductive death" mean?

Introduction:

Line 46: What does "loss of function" refer to specifically?

Line 49: Clarify what "Among others" is referring to here. We assume it means explanations / theories

Lines 55-73: Should also cite: Tschinkel, Insect. Soc. (2017) 64:285-296: DOI 10.1007/s00040-017-0544-0

Line 75: Cite original Wheeler (1911) paper (and if appropriate the Holldobler and Wilson book The Superorganism).

Line 82-83: The RNAseq do not necessarily serve the purpose of identifying which mechanisms are failing, but instead provide evidence that the queens are experiencing delayed senescence. Something that was not obvious in the previous version. Therefore, this should be mentioned in the manuscript.

Short methods:

Integrate into Intro, it's not part of the eLife format to have a "Short Methods" section and shorten it.

Results:

Line 125-126: Please clarify the meaning of "This confirms a previous study, which did not find a causal correlation between queen's lifespan and fecundity (Schrempf et al. 2017)." We are not clear how manipulation of # of workers in the colony leads to this specific conclusion. The result can be interpreted more conservatively as "no relation between queen's lifespan and colony size."

Lines 127-130: As in the previous round of revision, this is not necessary for the impact of the paper and we ask you to remove this, especially in light in the Conclusion section on Lines 281-286, stating it is not important whether it occurs maternally, during embryogenesis, larval development, etc …

Lines 140-141: Here is an opportunity to rule out the sperm depletion hypothesis. It is also redundant with Lines 172-174

Line 82and83 and 177and178: please make the question RNAseq analysis explicit. The question of "how" the queens are dying cannot be answered solely by RNaseq. Instead, focus the question on responding to Reviewer 2/3's critique, which is to distinguish between sperm limitation and senescence and to confirm that there is both reproductive and body senescence. Would we have expected that both head-thorax and gaster to both senesce? We would have expected that head-thorax but not reproductive tissue to senesce. Therefore making the question explicit and some discussion of the predictions and findings would improve this section considerably.

Finally, how confident are we that the expression changes presented really indicate ageing. I would have appreciated more detailed justification. For example, you could include a table with the functions of key genes that are differentially expressed, and then we could see what down or up regulation in old queens would mean. Alternatively, you could choose a couple of examples in the text and expand upon them. For example, we understand that gene co-expression networks are associated with aging, but why? This is a key result that distinguishes thorax-head and gaster transcriptome and is central to your thesis later that thorax-head deteriorates physiologically while gaster maintains some functionality. Therefore, some more justification would allow the reader to better appreciate these changes. See also comment on Lines 247-249.

Discussion:

Abstract and Line 246: Please remove the word "dramatic" before physiological changes, since you only have a 2-point and not 3-point comparison.

Lines 247-249: This statement is not well supported since they were not statistically different: "The changes are stronger in the gaster, which contains the reproductive organs and most of the digestive system, but to a similar extent occur in head and thorax, containing most neuronal and muscle tissue." Please refer explicitly to the evidence supporting this statement.

Line 251-252: Please provide some examples of these genes.

Conclusion:

Lines 287-297: Once again, these predictions are placing singular focus on the queen, without taking workers into account, which in our view, is antithetical to the superorganism perspective. Without direct supporting evidence, you cannot rule out that the interaction between queens and workers is just as important to the superorganism perspective as is the queen herself. In our view, we ask you to maintain a neutral stance on this point and ask that you maintain the data on worker longevity that you removed in the previous version and revise the text.

https://doi.org/10.7554/eLife.74695.sa1

Author response

Essential revisions:

The editors and both reviewers agree that this dataset has potentially important implications for understanding superorganism development and the long lifespans of social insects. Therefore, it will be of interest to a broad audience of biologists. Both reviewers have substantive comments on data interpretation and writing, as well as requests for some additional information and statistical analyses. Specifically:

1) Data interpretation:

a) A particular concern regarding the interpretation of the current data set is that it does not necessarily provide an ultimate explanation for the long life of the queens: an alternative, more conservative, interpretation of the current data set is that it only provides documentation of negligible reproductive senescence. According to Reviewer 2, the long lifespans of queens is strongly correlated too, but is not necessarily caused by, the capacity to reproduce late in life. For example, the long lifespans of queens could be ultimately caused by a shifting of reproductive costs to workers, or by living in protected environments. Therefore, in Reviewer 2's view, the current data set can only address the question of why social insect queens lack reproductive senescence. In other words, why do social insect queens retain the capacity to reproduce late in life? In my view, you (the Authors) should be very careful to distinguish ultimate from proximate causes in this study and be very clear about the precisely question this current data set can or cannot address.

The long lifespans of queens can be explained ultimately as a consequence of higher fitness returns late in life, and could be caused proximately by shifting costs to the workers. This is shown by our study, as queens with 30 workers are capable of rearing more sexuals than queens with 10 workers. Additionally, it is now known that living in a more protected environment does not mean a decrease in mortality (see Moorad Trends Ecol Evol 2020).

We now make clear that aging has two dimensions following Baudisch (2011) and Baudisch and Scott (2019) (now in Line 60). The first dimension “pace” refers to factors describing the time-scale (i.e. the long lifespans of reproductive of social insects). By doing so, we explain that this dataset was not obtained to answer directly questions related to the pace of aging. Instead, this dataset was obtained to answer questions related to the “shape” of aging. So, how is relative mortality and fertility distributed along the ant queens’ life, and how this directly affects the age-specific selection strength. Nevertheless, this is a fundamental start to understand the “pace” of aging. Hopefully with this we now make clear to the readers about the reach and limitations of this study.

b) Following on point (a) above, Reviewer 3 raises another particular point of concern, which is the lack of consideration for the alternative hypothesis that the drop in fertility later in life is caused by sperm limitation and not reproductive senescence. As Reviewer 3 points out, unless you can rule out this alternative hypothesis, it must be explicitly considered in the main text and would change the interpretation of the data, but not necessarily the relevance of the study, especially if you are willing shift the focus instead on understanding intrinsic properties of superorganism development and less on ageing per se.

In line 188 we explain that sperm limitation is not a plausible explanation for the drop of egg production late in life, as we only notice a final increase of solely male pupae production at the end of life in one colony. Queens were mated with wingless males which should transfer an excess amount of sperm (Schrempf and Heinze, 2008), and we show a high production of diploid queen pupae at the end of life.

c) Another issue of particular concern is the unnecessary focus on queen control. The inflection point, where colonies produce more queen offspring late in life, appears to reflect an intrinsic property of colony (superorganism) development that is independent of the colony-size, lifespan or worker investment. How this inflection is mechanistically controlled is a question you should address in the Discussion and should offer several hypotheses. For example, three possible hypotheses are: (1) An intrinsic colony-level threshold that is be controlled by queens and workers, such that workers change their feeding/foraging behavior in response to a queen pheromone once the queen reaches a certain age; (2) Queen embryos are produced all the time, like in the genus Monomorium, but are killed by workers until the queen reaches a certain age and changes its pheromone status; (3) as you have suggested, the threshold could be genetically determined in the queen regardless of worker input. Therefore, you should remove the focus on the assumption of queen control and provide a more systems-level (superorganismal) explanation for the inflection for producing queen offspring independent of colony size in the discussion.

As here suggested, the exact mechanism of determination is not necessary for the discussion of the results. We refrain from discussing the three hypotheses, as we are sure that caste is determined in the embryo. In addition, an older study (Suefuji et al., 2008 Biol Lett.) showed that queens not only produce males faster in multiple queen colonies, but also female sexuals. This was not discussed in that paper, but also points to queen control. Nevertheless, we state clearly that this is not the interesting point, but whether queen and worker interests are aligned (Line 310).

We have a manuscript in preparation that suggests that workers do not discriminate between queen and worker larvae, and that these are also chemically not distinguishable. This supports that interests are aligned.

On a sidenote: Queen larvae in Monomorium are usually culled in the presence of a queen. In our species, and in this experiment, queens are produced all the time.

2) RNAseq data: If you feel that RNAseq data can strengthen your arguments in light of the comments / critiques on data interpretation raised above, then you should add these data to the manuscript. If you decide to do so, please make sure to explain why this is important both in the cover letter and in the response to reviewers.

We have included RNASeq data of 36 queens. These ant queens belonged to the same population survey. The data compare fit queens with queens close to death, and is a substantial extension to the study. Specifically, this is a fundamental point in the discussion of how selection strength is maintained until the sexual investment peak is reached in ant queens.

3) Reviewers 2 and 3 require some missing statistics as well as analyses, which should be straight forward to perform.

We have addressed this point in the revised version.

4) The manuscript, as presently written, leaves many open questions that need clarification, therefore there is a need for major revision to respond to the queries of Reviewer 2 and 3 (see Reviewer comments below).

Responded below.

5) You should provide additional information about this species' life history (see Reviewer comments below).

We have included a “Short methods” section (Line 94-127) which now includes a small section about the species’ life history.

Reviewer #1 (Recommendations for the authors):

[…]

I will provide comments in a sequential order, as they appeared while reading the manuscript.

1. You must include the species name either in the title or in the abstract. Currently it reads as if you followed 102 colonies of random species of ants.

Now the species name is included in the title and abstract.

2. L 21 advanced not advance.

Done.

3. "Furthermore, mortality decreased with age in queens and workers, supporting the hypothesis that aging trajectories are adaptive." I didn't understand this sentence while reading the abstract and reading the reminder of the paper didn't make clearer. What is meant here? Why would decrease in mortality suggest that "aging trajectories" are adaptive? What are "aging trajectories"?

To avoid confusion with the terminology we have decided to use the terminology in Baudisch (2001) and Baudosch and Scott (2019) about aging as the consequence of “pace” and “shape” of demographic trajectories (mortality and fertility), now clarified in Line 60. This part has also changed, as we now calculated relative mortality as a function of age.

4. "We argue that selection for late life reproduction delays the occurrence of the selection shadow, leads to the apparent absence of senescence in ants and underlies the evolution of long lifespans." This is very confusing – if "selection shadow" is delayed, how come there is an apparent absence of senescence? "Delayed" suggests that senescence will occur but at a later point. Also, this sentence is a bit clumsy and contains grammatical errors.

We agree and we have reformulated the abstract. We meant that senescence is apparently absent, as it can only be perceived in a short window of time shortly before death. We have changed this throughout the manuscript.

5. "Together, this highlights the unique aging patterns of social insects in comparison to other iteroparous species across the tree of life." Why? Is there evidence that negligible senescence is unique characteristic of social insects? Also, "aging patterns" is unclear.

We have addressed this now in the introduction, in which we explain that aging theory predicts that senescence has its onset after maturity, and this is seen in iteroparous species. Here we show that the onset of senescence in ant queens is past that point and is delayed.

6. Introduction: I found it very confusing that the paper on ageing starts with the discussion of reproductive queens and workers of social insects. You need a general introduction paragraph. The most straightforward solution is to take the first paragraph from Conclusions (LL 138-151) and start your main text with this.

Aging and reproduction are life-history traits that are interconnected. Now we start with two general paragraphs that bring up both terms.

7. Introduction: "It is argued that the costs of reproduction (energy intake, brood care) are outsourced to the workers, but this is a necessary and not sufficient condition for the evolution of long lifespans in social insects and does not explain these aging patterns." Why? There is no explanation as to why this is not sufficient, and it is absolutely not clear to me – please explain.

It has been argued that, if costs are not overtaken by workers (in terms of protection, energy intake, etc), this will result in increasing levels of extrinsic mortality experienced by ant queens. And this has long been used as an explanation for the long lifespans of social insects. But as already explained by Hamilton and later recapitulated by Moorad and colleagues (now in line 251), age-independent extrinsic mortality cannot explain the pace of aging. Therefore, this cannot be a sufficient explain on why social insects live so long. Additionally, our study shows that an increase in worker force does not affect queens’ lifespan, confirming Moorad et al.

8. LL 45-58 – I like this paragraph and I think it presents good arguments in favour of the novelty provided by this study

Thanks.

9. L 58 – "aging, senescence" – if you treat ageing and senescence as different terms, you have to explain why and provide definitions for each term.

Now we define aging in line 44, and senescence in line 37.

10. L 58 – "investment" – what investment? Reproduction? Somatic maintenance? This is sentence is remarkably unclear.

Now we specify reproductive investment (Line 78).

11. I think that "superorganism" approach, while still debated in the literature is entirely justified in this system.

12. L81 – While not directly relevant as you did not find a correlation, how did you test for causality of correlation? Your experimental treatment arguably had no effect neither on egg production nor on queen lifespan?

We did not test for causality of correlation, we did not find an effect of our treatment in the egg productivity nor queen lifespan.

13. "Generalizing this finding, we suggest that lateLife investment and its effect on queen aging is an intrinsic property of superorganismality irrespective of the mode of caste determination and social organization." – this is unclear, I thought there is no effect on queen's ageing?

We explain in Line 310 that the mode of determination is not relevant. What is important is whether queen and worker interests are aligned, so whether there is conflict over sex and caste of the offspring. This can of course shape the pace of aging.

14. I think you should analyse age-specific mortality by selecting the best-fit model in a specialized demography analyses package, such as, for example, BASTA.

We prefer to show the relative mortality as a function of age as in Jones et al. (2014), as this gives valuable information about the specific age at which mortality increases over average life-time mortality (when mortality passes the threshold of 1 in Figure 4.) As suggested by the reviewer we also now include age-specific mortality of the best-model fitted using BaSTA and the estimated parameters in the supplement (Figure 4 —figure supplement 1, Supplementary File 1H and 1I), which does not add additional insights.

15. My understanding is that queens live for 40-50 weeks max (Figure S3). Figure 4 suggests that from week 30 onwards the production of eggs, worker pupae and queen pupae declines. This suggests that while queen mortality declines in late life, so does queen reproduction. So, queens of this species do show reproductive senescence?

Yes, now better shown in Figure 2 —figure supplement 1.

Yes, your data suggests that relative investment into reproduction (queen worker ratio) increases with age, but the absolute number of queens declines with age? To me, this suggests an interesting result from the life-history theory perspective – increased investment in reproduction with reduced residual reproductive value, but not necessarily the absence of reproductive senescence. Perhaps I didn't understand the results, please clarify.

Yes, we agree with you. We hope now it becomes clear from the text that queens do experience reproductive senescence.

16. L155 – this relative mortality is very unclear to me. I think you need proper mortality rate analyses in BASTA and estimate different age-specific mortality parameters.

See comment 14.

17. L 156 'random' offspring?

The word ‘random’ has been removed as it did not add new information.

18. "In turn, the short selection shadow might have led to the delay of a distinguishable senescence phase." I think I understand what you mean but this is very unclear, please re-write.

We have rewritten the respective parts. The RNAseq data reveal a clear senescence phase, evident at the end of life and not in middle aged queens. Clearly, only a time course analysis can reveal whether senescence occurs earlier, the fecundity data however speak against a prolonged phase of senescence.

Reviewer #2 (Recommendations for the authors):

Line 42: I'd like clarification on your argument here. Why is this not sufficient in your opinion? I don't see any alternative explanation, and workforce dynamics is central to your experimental design. After all, you are manipulating workforce size to expose trade-offs in queen decisions, so how does it not ultimately come down to the workforce?

As answered in the point 7 for reviewer 1: It has been argued that, if costs are not overtaken by workers (in terms of protection, energy intake, etc), this will result in increasing levels of extrinsic mortality experienced by ant queens. And this has been long used as an explanation for the long lifespans of social insects, but as explained by Hamilton and later recapitulated by Moorad and colleagues (now in line 251), age-independent extrinsic mortality cannot explain the pace of aging. Therefore, this cannot be a sufficient explanation why social insects live so long. Additionally, our study shows that an increase in worker force does not affect queens’ lifespan.

Line 72: Why not just say something simple like "compare queen and worker mortality". Throwing in a new term like aging trajectories makes it seem like you are doing something different.

We have removed the worker data in order to focus on the queen’s data.

Line 92: I would have liked to see more of your data on the male production, even if they were less common than queens, 15% of your reproductive output is not negligible. Does male production fit the same lateLife pattern as queen production? Do you see the same differences between treatments when you analyze combined alate production? I expect the 10-worker treatment may invest more in males, since they are usually cheaper to produce. I know this isn't the question you are asking but it may explain some of the drop in queen production in that treatment.

We have now included relative male production in Figure 4, and it follow the same late pattern as queen production. The lifetime production of wingless males is very low (mean: 0.36, median: 0, N = 99 queens). We did not test for winged males (there is a male polyphenism in Cardiocondyla) specifically as there were almost no produced. Further, it is not clear whether wingless males are low-cost investments, as males are under strong sexual selection, adapted for fighting with rivals, have longish lifespans and lifelong spermatogenesis.

Line 154: I don't understand why you seem to downplay the clear influence that the workforce has on this dynamic. From my perspective, your data beautifully shows how these queens transition from worker production to alate production as they age, and that colonies with low worker numbers cannot invest as much into reproduction, which is expected in a species that must establish a robust colony before devoting resources to reproduction. If a queen survives to a relative old age, she will also have a large, established workforce that can raise a lot of alates, so of course we will see longevity selected for at that point (as indicated by your own data, if she reaches that age without a sufficient workforce, her alate production will suffer).

We have now included in Line 262 “The magnitude of such investment (i.e. number of queen pupae produced) is affected by the number of workers available”, to acknowledge the role of worker size, so this message is not lost from the discussion.

Another important omission in the conclusion is the inevitable problem of sperm limitation. You frame the system as generating a directional push towards negligible senescence but it is more likely stabilizing selection that produces a sufficient workforce and then devotes as much remaining sperm as possible towards queen production. Selection will ensure that an established queen lives long enough to exhausts her sperm storage, which would contribute to the short selection shadow you discuss. The increase in relative mortality at 30 weeks may just represent the optimal lifespan of the colony. In figure 2A we see the peak in queen production at that time point, with both worker and queen production falling dramatically by 40 weeks, potentially signaling sperm depletion. Upon closer inspection of S3, it looks like there were only a few colonies that survived past week 38, but then persisted for another couple months. I have to wonder how reliable the drop in relative mortality is after 40 weeks when it represents only a few colonies. Of course, this is my interpretation based off of the figures, so please include the raw numbers if I am mistaken.

The decrease of queen and worker production due to sperm limitation is highly unlikely as we only found one case of sperm depletion in the whole data set (now reported in line 154). In the new version of the manuscript, we have included that in this species queens generally mate once, and that ergatoid males transfer an excess amount of sperm making very uncommon the cases of sperm depletion.

After revision of the relative mortality as a function of age according to Jones et al. 2014, we find an increase of mortality above average after 30 weeks in ant queens, corroborated by the suggested BaSTA analysis. Additionally, raw numbers of exact number of surviving queens at each age are reported in the Source Data 7.

Line 212: Is this species monogamous in the wild? If so, this should be stated because if they are normally polyandrous, having only one mate may introduce sperm limitation.

We have now included this information in Line 98 as it follows: “Virgin queens usually mate once with related wingless males inside the natal colony (Heinze and Hölldobler 1993; Schmidt et al. 2016; Heinze and Hölldobler 2019), generally stay in the nest, and new colonies are formed by budding of colony fragments.”

[Editors' note: further revisions were suggested prior to acceptance, as described below.]

Although the transcriptomic data were analyzed properly, the conclusions that can be drawn from these data are limited, primarily because the difference between gaster and head/thorax is not easy to interpret. However, we think that these data do add something to the paper by validating some of the conclusions drawn from indirect observations. Therefore, the RNAseq do not necessarily serve the purpose of identifying which mechanisms are failing, but instead provide evidence that the queens are experiencing delayed senescence, something that was not obvious in the previous version. Please keep the interpretation of RNAseq simple and revise the manuscript appropriately to make its purpose clearer (see comment below on Line 82-83 and 178 and 179 below).

We realize that there are problems with the data. For one, as you point out, it is futile to speculate which processes fail in old queens, as there are many. Second, we only scratch the surface in a non-model insect and rely on assumptions. Third, we only compare two time points, thus ‘changes’ could be transient. Lastly, we deduce a pattern from two different studies. Harrison et al. (2020) compare young with middle-aged queens just before their reproductive peak and find no signs of senescence, rather the contrary. Here, we compare middle-aged at their reproductive peak with prope mortem queens and show strong signs senescence. We hope to have revised the manuscript accordingly to make the purpose of the RNAseq data clear.

Furthermore, the singular focus on the queen is problematic, especially if you are attempting to adopt a superorganism perspective. First, the evidence presented that only the queen is in control is indirect and weak. Therefore, we ask you to maintain a neutral stance on this issue and acknowledge that the workers may play an important role in this process. Second, your arguments do not depend on whether determination is in the maternal, embryonic, or larval stage and we ask you to remove this mention as it is a non-central point in your argument (lines 127-130). Given these two points, we ask you to put back the data on worker longevity (although it can be presented in the supplemental data), because in our view, without direct supporting evidence that it is all about the queen, workers remain an important part of the equation. We therefore ask you to revise lines 127-130 in the Intro and 287-297 in the Conclusions section.

We have changed the respective section. Please understand, that coming from a background nested in the vast diversity of myrmecological life, the last author struggles with generalizations. Indeed, the focus is unnecessary and does not add to the argument.

Essential revisions:

Abstract:

Line 23-25: We ask you to add something to describe the inflection point as a threshold, an inflection point that changes in an abrupt, non-linear-like way which distinguishes it from a gradual linear-like increase.

We now describe the non-linear increase of sexual production: “By life-long tracking of 99 Cardiocondyla obscurior (Formicidae: Myrmicinae) ant colonies, we find that queens shift to the production of sexuals in late life regardless of their absolute lifespan or worker investment”

What does "reproductive death" mean?

We now state that it is the development of massive pathology following reproductive effort as:

“Furthermore, RNAseq analyses of old queens past that fitness peak showed the development of massive pathology while still being fertile, leading to rapid post-reproductive death.”

Introduction:

Line 46: What does "loss of function" refer to specifically?

Now should be clear as “…and may lead to loss of reproductive performance and survival (i.e. senescence).”

Line 49: Clarify what "Among others" is referring to here. We assume it means explanations / theories.

Done.

Lines 55-73: Should also cite: Tschinkel, Insect. Soc. (2017) 64:285-296: DOI 10.1007/s00040-017-0544-0.

Done.

Line 75: Cite original Wheeler (1911) paper (and if appropriate the Holldobler and Wilson book The Superorganism).

Done.

Line 82-83: The RNAseq do not necessarily serve the purpose of identifying which mechanisms are failing, but instead provide evidence that the queens are experiencing delayed senescence. Something that was not obvious in the previous version. Therefore, this should be mentioned in the manuscript.

We included now the purpose of the RNAseq analysis to provide evidence of senescence in this prope mortem queens as:

“To assess if senescence was restricted particularly to the end of life we compared RNAseq data of 18 of such prope mortem (Lat. near death) queens (between 28-29 weeks old) and 18 middle-aged queens (between 19-21 weeks) which were in their peak of fertility (Figure 1 —figure supplement 3A-B.)”

Short methods:

Integrate into Intro, it's not part of the eLife format to have a "Short Methods" section and shorten it.

Done.

Results:

Line 125-126: Please clarify the meaning of "This confirms a previous study, which did not find a causal correlation between queen's lifespan and fecundity (Schrempf et al. 2017)." We are not clear how manipulation of # of workers in the colony leads to this specific conclusion. The result can be interpreted more conservatively as "no relation between queen's lifespan and colony size."

We decided to remove this sentence as it does not add to the discussion.

Lines 127-130: As in the previous round of revision, this is not necessary for the impact of the paper and we ask you to remove this, especially in light in the Conclusion section on Lines 281-286, stating it is not important whether it occurs maternally, during embryogenesis, larval development, etc …

We deleted those lines.

Lines 140-141: Here is an opportunity to rule out the sperm depletion hypothesis. It is also redundant with Lines 172-174

We deleted lines 172-174 and left the ones in the result section.

Line 82and83 and 177and178: please make the question RNAseq analysis explicit. The question of "how" the queens are dying cannot be answered solely by RNaseq. Instead, focus the question on responding to Reviewer 2/3's critique, which is to distinguish between sperm limitation and senescence and to confirm that there is both reproductive and body senescence. Would we have expected that both head-thorax and gaster to both senesce? We would have expected that head-thorax but not reproductive tissue to senesce. Therefore making the question explicit and some discussion of the predictions and findings would improve this section considerably.

We changed the Introduction, and explain now in the beginning of the Results section the purpose of the RNAseq analysis:

“To determine if queens show signs of reproductive senescence and loss of physiological function, we analyzed gene expression data of prope mortem queens exhibiting decreasing egg laying rates, and middle-aged queens that were at their peak reproductive performance.“

About the comparison between tissues, we included:

“We subjected the head plus thorax and the gaster (see methods for terminology, Figure 1 —figure supplement 5) to RNAseq separately, to asses if reproductive tissue shows different physiological wear and tear than head-thorax tissue”

Finally, how confident are we that the expression changes presented really indicate ageing. I would have appreciated more detailed justification. For example, you could include a table with the functions of key genes that are differentially expressed, and then we could see what down or up regulation in old queens would mean. Alternatively, you could choose a couple of examples in the text and expand upon them. For example, we understand that gene co-expression networks are associated with aging, but why? This is a key result that distinguishes thorax-head and gaster transcriptome and is central to your thesis later that thorax-head deteriorates physiologically while gaster maintains some functionality. Therefore, some more justification would allow the reader to better appreciate these changes. See also comment on Lines 247-249.

We considered adding such table of most important genes (besides the complete table now in Source data file 9 and 10), but realize that it is impossible to avoid cherry picking. This is the reason why we present a functional enrichment analysis, to give a clear perspective of the functions involved. Even so, the GOTerm results is quite extensive (see Supplement file 1B-E) and picking just a few functions to describe the changes with age does not seem very useful.

Discussion:

Abstract and Line 246: Please remove the word "dramatic" before physiological changes, since you only have a 2-point and not 3-point comparison.

Done.

Lines 247-249: This statement is not well supported since they were not statistically different: "The changes are stronger in the gaster, which contains the reproductive organs and most of the digestive system, but to a similar extent occur in head and thorax, containing most neuronal and muscle tissue." Please refer explicitly to the evidence supporting this statement.

We have added that the tissues show differences in the number of duplicated reads (which does not seem to be a technical artifact). Nevertheless, we state that processes occurring in both tissues are similar in the Discussion as:

“Given these discrepancies between tissue types, 104 GO terms were significantly enriched in both tissues, of these 44 in prope mortem queens (Supplementary File 1F) and 60 in middle-aged queens (Supplementary File 1G). Thus, signs of similar physiological pathologies occur in reproductive and non-reproductive tissue.”

Line 251-252: Please provide some examples of these genes.

Done.

Conclusion:

Lines 287-297: Once again, these predictions are placing singular focus on the queen, without taking workers into account, which in our view, is antithetical to the superorganism perspective. Without direct supporting evidence, you cannot rule out that the interaction between queens and workers is just as important to the superorganism perspective as is the queen herself. In our view, we ask you to maintain a neutral stance on this point and ask that you maintain the data on worker longevity that you removed in the previous version and revise the text.

We reintegrated the data on worker longevity to maintain the perspective on the superorganism, and agree there could be indirect interactions between queen and worker longevity.

https://doi.org/10.7554/eLife.74695.sa2

Article and author information

Author details

  1. Luisa Maria Jaimes-Nino

    Zoologie/Evolutionsbiologie, Universität Regensburg, Regensburg, Germany
    Contribution
    Data curation, Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review and editing
    For correspondence
    jaimes.luisa@outlook.com
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-1186-1837

  2. Jürgen Heinze

    Zoologie/Evolutionsbiologie, Universität Regensburg, Regensburg, Germany
    Contribution
    Conceptualization, Funding acquisition, Project administration, Supervision, Validation, Writing – review and editing
    Competing interests
    No competing interests declared

  3. Jan Oettler

    Zoologie/Evolutionsbiologie, Universität Regensburg, Regensburg, Germany
    Contribution
    Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing – original draft, Writing – review and editing
    For correspondence
    joettler@gmail.com
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-8539-6029

Funding

Deutsche Forschungsgemeinschaft (OE549/2-2)

  • Jan Oettler
  • Jürgen Heinze

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank Vera Ermer, Benjamin Dofka, Julia Haschlar, Lena-Marie Süß, and Judith Weber for help with the experiment, Eva Schultner, Tomer Czaczkes, Boris Kramer, Andrew Bourke, two anonymous reviewers for comments. This study was funded by the Deutsche Forschungsgemeinschaft (OE549/2-2). No funding sources were involved in study design, data collection and interpretation, or the decision to submit the work for publication.

Senior Editor

  1. Claude Desplan, New York University, United States

Reviewing Editor

  1. Ehab Abouheif, McGill University, Canada

Publication history

  1. Preprint posted: February 25, 2021 (view preprint)
  2. Received: October 13, 2021
  3. Accepted: March 21, 2022
  4. Version of Record published: April 6, 2022 (version 1)

Copyright

© 2022, Jaimes-Nino et al.

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.

Metrics

  • 800
    Page views
  • 128
    Downloads
  • 3
    Citations

Article citation count generated by polling the highest count across the following sources: PubMed Central, Crossref, Scopus.

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. Luisa Maria Jaimes-Nino
  2. Jürgen Heinze
  3. Jan Oettler
(2022)
Late-life fitness gains and reproductive death in Cardiocondyla obscurior ants
eLife 11:e74695.
https://doi.org/10.7554/eLife.74695

Categories and tags

Further reading

    1. Developmental Biology
    2. Evolutionary Biology

    Mutation of vsx genes in zebrafish highlights the robustness of the retinal specification network

    Joaquín Letelier, Lorena Buono ... Juan R Martínez-Morales
    Research Article

    Genetic studies in human and mice have established a dual role for Vsx genes in retina development: an early function in progenitors’ specification, and a later requirement for bipolar-cells fate determination. Despite their conserved expression patterns, it is currently unclear to which extent Vsx functions are also conserved across vertebrates, as mutant models are available only in mammals. To gain insight into vsx function in teleosts, we have generated vsx1 and vsx2 CRISPR/Cas9 double knockouts (vsxKO) in zebrafish. Our electrophysiological and histological analyses indicate severe visual impairment and bipolar cells depletion in vsxKO larvae, with retinal precursors being rerouted toward photoreceptor or Müller glia fates. Surprisingly, neural retina is properly specified and maintained in mutant embryos, which do not display microphthalmia. We show that although important cis-regulatory remodelling occurs in vsxKO retinas during early specification, this has little impact at a transcriptomic level. Our observations point to genetic redundancy as an important mechanism sustaining the integrity of the retinal specification network, and to Vsx genes regulatory weight varying substantially among vertebrate species.

    1. Evolutionary Biology
    2. Genetics and Genomics

    The impact of local genomic properties on the evolutionary fate of genes

    Yuichiro Hara, Shigehiro Kuraku
    Research Article

    Functionally indispensable genes are likely to be retained and otherwise to be lost during evolution. This evolutionary fate of a gene can also be affected by factors independent of gene dispensability, including the mutability of genomic positions, but such features have not been examined well. To uncover the genomic features associated with gene loss, we investigated the characteristics of genomic regions where genes have been independently lost in multiple lineages. With a comprehensive scan of gene phylogenies of vertebrates with a careful inspection of evolutionary gene losses, we identified 813 human genes whose orthologs were lost in multiple mammalian lineages: designated ‘elusive genes.’ These elusive genes were located in genomic regions with rapid nucleotide substitution, high GC content, and high gene density. A comparison of the orthologous regions of such elusive genes across vertebrates revealed that these features had been established before the radiation of the extant vertebrates approximately 500 million years ago. The association of human elusive genes with transcriptomic and epigenomic characteristics illuminated that the genomic regions containing such genes were subject to repressive transcriptional regulation. Thus, the heterogeneous genomic features driving gene fates toward loss have been in place and may sometimes have relaxed the functional indispensability of such genes. This study sheds light on the complex interplay between gene function and local genomic properties in shaping gene evolution that has persisted since the vertebrate ancestor.

    1. Evolutionary Biology
    2. Plant Biology

    Environment as a limiting factor of the historical global spread of mungbean

    Pei-Wen Ong, Ya-Ping Lin ... Cheng-Ruei Lee
    Research Article

    While the domestication process has been investigated in many crops, the detailed route of cultivation range expansion and factors governing this process received relatively little attention. Here using mungbean (Vigna radiata var. radiata) as a test case, we investigated the genomes of more than one thousand accessions to illustrate climatic adaptation’s role in dictating the unique routes of cultivation range expansion. Despite the geographical proximity between South and Central Asia, genetic evidence suggests mungbean cultivation first spread from South Asia to Southeast, East, and finally reached Central Asia. Combining evidence from demographic inference, climatic niche modeling, plant morphology, and records from ancient Chinese sources, we showed that the specific route was shaped by the unique combinations of climatic constraints and farmer practices across Asia, which imposed divergent selection favoring higher yield in the south but short-season and more drought-tolerant accessions in the north. Our results suggest that mungbean did not radiate from the domestication center as expected purely under human activity, but instead the spread of mungbean cultivation is highly constrained by climatic adaptation, echoing the idea that human commensals are more difficult to spread through the south-north axis of continents.