- Research article
- Open Access
Analysis of the coding sequences of clownfish reveals molecular convergence in the evolution of lifespan
BMC Evolutionary Biologyvolume 19, Article number: 89 (2019)
Standard evolutionary theories of aging postulate that reduced extrinsic mortality leads to evolution of longevity. Clownfishes of the genus Amphiprion live in a symbiotic relationship with sea anemones that provide protection from predators. We performed a survey and identified at least two species with a lifespan of over 20 years. Given their small size and ease of captive reproduction, clownfish lend themselves as experimental models of exceptional longevity. To identify genetic correlates of exceptional longevity, we sequenced the transcriptomes of Amphiprion percula and A. clarkii and performed a scan for positively-selected genes (PSGs).
The PSGs that we identified in the last common clownfish ancestor were compared with PSGs detected in long-lived mole rats and short-lived killifishes revealing convergent evolution in processes such as mitochondrial biogenesis. Among individual genes, the Mitochondrial Transcription Termination Factor 1 (MTERF1), was positively-selected in all three clades, whereas the Glutathione S-Transferase Kappa 1 (GSTK1) was under positive selection in two independent clades. For the latter, homology modelling strongly suggested that positive selection targeted enzymatically important residues.
These results indicate that specific pathways were recruited in independent lineages evolving an exceptionally extended or shortened lifespan and point to mito-nuclear balance as a key factor.
The lifespan of vertebrate species spans more than two orders of magnitude, from a few months for annual killifishes  to several centuries for the Greenland shark . Understanding the genetic architecture underlying these differences is a major challenge that may deliver new insights into the mechanisms controlling evolution of lifespan and human longevity.
Next-generation sequencing technology can provide genome-scale sequence information for a large number of species and has revolutionized evolutionary genomics. A particularly useful approach to identify the genetic architecture of evolutionary novelties is the analysis of positive selection. A frequently used method for detecting positive selection relies on the comparison of the sequence of protein-coding genes in related clades where one of the clades evolved the trait of interest, in this specific case exceptional lifespan. To date, several different mammalian taxa/clades were analysed with this approach with the purpose of identifying sequence changes associated to evolution of longevity: the elephant, the bowhead whale, bats and mole-rats [3,4,5,6,7,8]. In addition to the analysis of positive selection, other phylogeny-based methods have been used to investigate genetic changes related to the evolution of longevity, e.g. [9,10,11,12]. The different analyses delivered interesting candidate genes and pathways that underwent accelerated molecular evolution in coincidence with evolution of exceptional lifespan. A major drawback of all these approaches – including the analysis of positive selection – is that long-lived mammals are difficult or impossible to be kept in captivity and manipulated experimentally. This creates the need for a long-lived vertebrate that is small in size, easily adaptable to captive life, can be bred in large numbers and therefore could become a convenient experimental model organism.
The standard evolutionary theories of aging – namely the Mutation Accumulation , Antagonistic Pleiotropy  and Disposable Soma theory  – predict that low extrinsic mortality conditions lead to the evolution of slow senescence and increased lifespan. Some examples that confirm these theories are the exceptional longevity of vertebrate species under low predation risk since they are chemically protected [16, 17], adapted to an arboreal life  or found in protected environments such as caves , respectively. On the other hand, annual fishes of the genus Nothobranchius provide an example of how increased extrinsic mortality conditions lead to the evolution of accelerated senescence and a short lifespan [20,21,22]. Analysis of positive selection in annual killifishes revealed a potential link between the evolution of genes governing mitochondrial biogenesis and the evolution of lifespan .
All clownfish species (genus Amphiprion) evolved a specific adaptation that allows them to live in symbiosis with sea anemones. Symbiosis evolved in the last common ancestor of the clownfish and clownfish represent a monophyletic group in the Pomacentridae family (damselfishes) . In the Indo-Pacific Ocean, clownfishes are found in association with one or more sea anemone species and a large variation in host usage exists [25,26,27]. Fish that feel threatened by predators immediately seek protection by the anemone’s tentacles; without that symbiosis, fish are readily attacked by predators [28,29,30]. Therefore, clownfishes are protected from predators through the presence of anemones . Hence, the overall mortality rate of the clownfish is low as compared to other coral reef fishes or other tropical species of Pomacentridae of the same size [28, 31,32,33,34].
All clownfishes are born as males and develop, through protandrous hermaphroditism, into females: in a colony, only the dominant pair contributes to the reproduction of the colony . Other individuals of the colony are non-breeding males. Studies in the wild have shown that natural mortality of adult clownfishes can be very low: during the period 2011–2013, the average biannual mortality rate per capita varied, depending on the study site, between 0.18 and 0.49 for juveniles, 0.09 and 0.44 for males, and 0.19 and 0.55 for females . Predatory pressure differs in different stages of adulthood and is increased for non-breeding males .
These fishes are small in size (less than 10 cm for the smallest species) and the closely-related species A. percula and A. ocellaris are popular and hardy aquarium fishes, are bred in large numbers for the aquarium trade  and are subject to selective breeding to fix specific pigmentation patterns so that a number of different captive strains are available. For these reasons, the clownfish could become the first experimental model for long-lived vertebrates.
In order to identify the genetic basis of adaptations linked to the clownfishes’ exceptional lifespan, we performed a positive selection analysis. This analysis requires the identification of the closest related taxon that does not possess the trait of interest in order to exclude events of positive selection that predate the evolution of this trait .
Other species of damselfish evolved an inter-specific mutualistic relationship with branching corals [39, 40]. In this case, corals are used by fishes as shelter that can provide protection from predators and a safe area for egg laying [41, 42]. Among the family Pomacentridae, Chromis viridis shows an interesting relationship with a wide range of scleractinians [43, 44]. Despite the presence of a favourable microhabitat, C. viridis are predated by a wide range of generalist predator species. Hixon and Carr  suggested there is a clear relationship among transient and benthic predators and damselfish mortality: damselfishes that search for protection in the shelter from transient predators are susceptible to attack by resident benthic predators and vice versa. In the presence of both groups of predators, mortality increases dramatically due to the lack of available refuge that expose Chromis to intense predation . Therefore, Chromis viridis represents a well-suited outgroup for our analysis because it shares with clownfish several general traits linked to benthic life and symbiosis with corals but it is subject to much higher predation rates (Fig. 1).
There are no reliable estimates of C. viridis lifespan in captivity. This is mainly due to the fact that these animals build large schools and it is not possible to identify them individually, as opposed to clownfishes that are normally kept as pairs and, if multiple pairs are maintained in the aquarium, each pair associates to a different symbiont making their distinction easy. For the same reason, the presence of a specific clownfish on a specific anemone can be assessed in multiple years to obtain the rates of per annum mortality , which is not possible for C. viridis. However, it is certain that tropical Chromis fishes undergo severe predation in the post-settlement phase . In one experiment where the total number of the population was assessed, a decrease in the population size of almost 80–90% in one year was documented. Notably, this not only affects juvenile fish, but also adults, indicating high adult mortality. This, combined with a very rapid growth (80% of maximum size reached within the first year) clearly indicates that these animals are short-lived in the wild . Indeed, they are considered a priori a model for short-lived reef inhabitants . Mortality data and age structure are available for the temperate species Chromis chromis (that, due to lower temperatures, may be longer-lived than tropical species). In this case, the natural adult annual mortality was estimated to be 0.72 and only < 7% of the individuals reached ages higher than 7 years .
Lifespan-data confirms that clownfish evolved exceptional long lifespans
Several anecdotic reports are present in the hobbyist literature on the exceptional lifespan of the clownfish. In order to obtain independent data on the captive lifespan of clownfish species, in 2016 we distributed a questionnaire to researchers working with clownfishes and to public aquaria across Europe (Table 1 and Additional file 1: Table S1) resulting in lifespan data for, in total, 114 individuals. Additionally, we surveyed existing literature. For six different species, at least one individual was reported to have lived more than 10 years and for two different species, A. melanopus and A. ocellaris, we obtained records of animals alive and actively spawning at an age of over 20 years, confirming hobbysts’ reports. This indicates that even the longest-lived individuals observed in captivity did not show reproductive senescence and were not approaching the limit of their lifespan.
More systematic data could be obtained for the species A. ocellaris (the most common species in the aquarium trade). The oldest cohort for which a record was available comprised 27 fish born in 2008 of which 25 were still alive in 2016.
We conclude that there is solid evidence that at least the species A. ocellaris and A. melanopus can live in captivity for more than two decades, making them the first teleost model of exceptional longevity.
Genes that were positively selected on the ancestral clownfish branch are enriched for aging-relevant functions
In order to perform genome-wide scans for positive selection, we obtained the transcriptomes of the species A. clarkii and A. percula based on our sequencing using methods previously described for the killifishes . Furthermore, we assembled clownfish transcriptomes from public read data of A. bicinctus, A. ocellaris and A. sebae. As the closest-related non-symbiotic species, we additionally sequenced the transcriptome of Chromis viridis, a closely-related Pomacentrid of 10 cm maximum size, i.e. comparable to that of the smallest clownfish species (Table 1), which is very abundant in coral reefs and lacks adaptations for symbiosis with sea anemones. Although data on C. viridis mortality in the wild are not available, the higher predation pressure on this species (see Background) makes it a suitable outgroup for our analysis. More distant outgroups were a selection of species from the series Ovalentaria, whose genomes are available in GenBank (see also ). We analysed positive selection on the branch leading to the last common ancestor (LCA) of all clownfish species (Fig. 1).
A total of 157 positively selected genes (PSGs) of 14,214 analyzed genes were identified in the LCA of the clownfishes (Additional file 1: Table S2). We tested for overrepresentation of gene ontology (GO, FDR < 0.1) and observed 19 biological processes enriched for PSGs (Table 2, Additional file 1: Table S3). A majority of these processes are of particular interest for aging research: altogether nine enriched processes are linked to the metabolism of xenobiotics, detoxification or glutathione metabolism, respectively. Interestingly, these processes were shown to be strongly up-regulated in experimental conditions favoring longevity such as dietary restriction and inhibition of the somatotropic axis making the animals more resistant to toxins [48,49,50,51]. Furthermore, experimental manipulation of mitochondrial translation, another enriched process, is known to increase lifespan in C. elegans  and variations in the expression of these genes is associated to lifespan variation in killifish  and mouse . We also tested formally the relationship between the results of the positive selection analysis on the LCA of the clownfishes and aging. For this, we checked whether processes whose gene expression was shown to be associated with maximum longevity in an analysis across 33 diverse mammal species , were also enriched during the positive selection analysis, which was the case (p = 0.06).
The positive selection analysis provides not only candidate genes but also candidate amino acids for follow-up studies. To exemplify this, we performed protein homology modeling for GSTK1 starting from the publicly available structures of the human dimeric apoenzyme (PDB 3RPP, ) and the rat dimeric enzyme with the bound GSH substrate (PDB 1R4W; ). The latter was used to assess on a structural basis the relationship of the six positively selected sites in the clownfish with those that are known to be involved in the enzyme’s function . Interestingly, also the LCA of Nothobranchius shows positive selection in GSTK1 contains, in addition, one site with high probability of positive selection in the LCA of Nothobranchius (Glu167, blue in Fig. 2). The selected site in Nothobranchius, however, is structurally remote to the functionally relevant sites. In contrast, we found that in clownfish two of three sites that were predicted with high probability to be positively selected (≥ 95%, Phe60, Met63, red in Fig. 2) and one of three sites with lower probability (41%, His64, orange in Fig. 2) belong to the same α-helical stretch of amino acids that lines the substrate access channel, contribute to the dimer interface (Asn61, Tyr65, Asp69, green in Fig. 2) as well as to the substrate binding sites (Lys62, turquoise in Fig. 2), respectively . The third site with a high probability to be positively selected is Glu88 (brown in Fig. 2). Glu88 is one of four amino acids at the entrance of the substrate access channel and situated in close proximity to Pro55, Pro56 and Pro87 (black in Fig. 2). The latter three are also part of the substrate access channel . We found another site positively selected with a lower probability in close proximity to the dimer interface (Lys177, orange in Fig. 2). This positive selection at particular positions related to enzymatic function invites the speculation that it might have a bearing on the enzymatic activity of the clownfish GSTK1.
Meta-analysis suggests that adaptation of mitochondrial biogenesis is a key player in evolution of lifespans
Recent observations of similar genes and pathways found to be affected by positive selection, both, in very long- and short-lived species led to hypotheses of antiparallel evolution acting on these entities [56, 57]. This means that functionally opposite selection pressures with regard to the tradeoff between fast growth and a long lifespan can result in adaptations of the same genes and pathways – in opposite functional directions.
This applies particularly for genes involved in mitochondrial biogenesis: Functions like “Mitochondrial large/small ribosomal subunit” (GO:0005762/ GO:0005763) and “Mitochondrial respiratory chain complex I” (GO:0005747) were found to be enriched for PSGs on ant branches that are associated with a 10- to 100-fold increase in lifespan [58, 59]. Similar genes, associated with mitochondrial functions, were found to be under positive selection in the bats P. poliocephalus and M. lucifugus that can reach lifespans of more than 20 and 30 years, respectively [60, 61]. On the other hand, also extreme reduction of lifespan on three killifish branches (Nothobranchius family) was associated with an enrichment of PSGs for a gene set that stands explicitly for mitochondrial biogenesis (, p < 10− 6).
We examined this hypothesis of lifespan-associated, antiparallel evolution of mitochondrial biogenesis genes by testing whether the PSGs identified on the clownfish ancestor are enriched in the same gene set that was used in the killifish study. Furthermore, we reanalyzed in the same regard PSGs that were identified across mole-rat branches on which lifespans were remarkably prolonged. In both cases – clownfish ancestor as well as long-lived mole-rats – we found again an enrichment of PSGs for genes involved in mitochondrial biogenesis (p = 0.007 and p = 0.097, respectively, Additional file 1: Table S7).
We further tested the hypothesis by using Fisher’s method to combine enrichment p-values across the results of the mentioned positive selection analyses in short-lived killifishes, long-lived mole-rats, and long-lived clownfishes. In this meta-analysis, 34 genes exhibited a signature of positive selection (FDR < 0.1) across species (Additional file 1: Tables S4, S5 and S6). Among the genes involved in mitochondrial biogenesis were TFB2M and MTERF, that are necessary for mitochondrial transcription, FASTKD5 and FASTKD2 whose gene products are required for the biogenesis of mitochondrial ribosomes, , as well as RARS2 coding for a mitochondrial tRNA-synthetase. Again, the 34 genes that were significant in the meta-analysis across three different analysis associated with lifespan-changes were enriched mitochondrial biogenesis gene set (p = 1.05*10− 5, Additional file 1: Table S7).
Among the other 15 PSGs genes showing evidence for positive selection, both, in the clownfish LCA and in meta-analysis were, e.g., LAMP2 and CD63 (also called LAMP3) which code for major protein components of the lysosomal membrane [64, 65]. In addition, CD63 gene expression was shown to predict the malignancy grade of many different tumor types [66,67,68,69,70] and the artificial prevention of the decrease of LAMP2 gene expression during aging in mice results in considerably reduced cell damage, as well as in liver functions in old mice that are indistinguishable from those in young mice . Finally, another interesting example that was identified as significant, both, in the clownfish LCA and in the meta-analysis, is GSTK1 encoding a glutathione-S-transferase that localizes to the peroxisome. GSTK1 was shown to be associated with diabetes type 2 which is another major aging related disease [72, 73].
As a negative control, we searched for PSGs on evolutionary branches that were not associated with drastic changes of lifespan. For this, we used the closely-related sister-taxons of the three above mentioned species: C. viridis as sister-taxon of the clownfish-ancestor, C. porcellus (guinea pig) for the mole-rats and Aphysemion striatum for killifishes. 188, 124 and 27 PSGs were identified, respectively (FDR < 0.05, Additional file 1: Tables S8, S9 and S10). On none of these branches, however, an enrichment of PSGs against the mitochondrial biogenesis gene set or any related gene ontology category was found. Neither were the PSGs that were identified as significant (FDR < 0.1) in the meta-analysis across the three branches enriched for the mitochondrial biogenesis gene set (Additional file 1: Table S11) (Table 3).
We have corroborated the evidence for exceptional longevity of clownfishes in captivity. The species A. ocellaris is bred in captivity and commercially available in large numbers and it could become the first broadly-available laboratory model for extended lifespan and exceptional longevity.
Analysis of positive selection has shown evolutionary convergence, both, with the exceptionally short-lived genus Nothobranchius and with exceptionally-long lived mole rats.
In particular, clownfishes and mole rats both show positive selection in two key proteins of the lysosome: LAMP2 and CD63. These results are consistent with the conserved up-regulation across tissues and species of genes coding for lysosomal proteins [74, 75], the widespread accumulation of lysosomal aggregates observed during aging  and earlier findings that associated selection on lysosomal genes with evolution of mammalian longevity . This suggests that lysosomal function is of key importance for evolution of exceptional longevity. Another interesting example of convergent evolution is GSTK1, which is positively selected in both the exceptionally-long and exceptionally-short lived fish clades. GSTK1 is involved in glutathione metabolism. Since detailed structures of this protein are available [54, 55], homology modelling was possible and it strongly suggests that positive selection targeted positions that are involved in the enzymatic function of the encoded protein. This speculation would have to be tested experimentally.
Finally, prominent signs of convergence were observed for genes and pathways related to biogenesis of mitochondrially-encoded proteins with the remarkable observation that MTERF is under positive selection in all three taxa. These findings parallel experimental evidence obtained in laboratory species where reduced expression of mitochondrially-encoded proteins prolongs lifespan [52, 53] and point to the key importance of mito-nuclear balance in the evolution of animal longevity. Further, the observation that the same pathway is under positive selection both in exceptionally short- and long-lived species indicate that the same genetic architecture underlies both evolution of longevity and evolution of compressed lifespan.
The clownfish represents the first vertebrate model of exceptional longevity that can be easily cultivated in captivity.
Analysis of positive selection on protein-coding genes reveals convergence on lysosomal protein with long-lived mole rats and convergence on mito-nuclear balance genes with both mole-rats and the exceptionally short-lived annual killifishes.
These data indicate that the same genetic pathways were repeatedly recruited in response to evolution of lifespan in either direction.
Clownfish lifespan estimation
The determination of the clownfish lifespan was performed through the distribution of an internet-based questionnaire to zoos and aquariums worldwide, requesting information on clownfish demographic details: (1) the various clownfish species maintained in captivity, (2) the number of individuals for each species, (3) if each individual is captive bred or not, (4) the year of acquisition and, if not still alive, death, and (5) the sex of each individual, if determined. The questionnaire was circulated in 2016 to international associations and organizations of zoos and public aquariums such as the European Association of Zoos and Aquaria (EAZA), the Association of Zoos and Aquariums (AZA), the European Union of Aquarium Curators (EUAC) and the World Association of Zoos and Aquariums (WAZA). Responses to our questionnaire were received from 5 zoos and aquariums as well as two private entities resulting in lifespan data for 114 individuals (see Additional file 1: Table S1 and Acknowledgments).
Experimental fish and sampling
Sub adult Amphiprion percula (total length, 45.2 ± 1.2 mm; Wt, 1.6 ± 0.1 g, n = 12), Amphiprion clarkii (total length, 46.4 ± 5.1 mm; Wt, 2.3 ± 0.9 g, n = 12) and Chromis viridis (total length, 43.0 ± 1.6 mm; Wt, 1.3 ± 0.1 g, n = 12), were used. Animals were acquired from local dealers and subjected to acclimation during one month in the facilities of the Marine Aquarium at the University of Murcia (Spain). Fishes were kept in groups under exactly the same conditions (temperature, 27 ± 1 °C; salinity, 24 ± 1, pH, 8 ± 0.2; dissolved oxygen, 6.5 ± 0.2 mg/L) and fed ad libitum four times a day a standard low-fat diet to match their requirements (composed by Mysis shrimp, enriched Artemia nauplii and red plankton).
Fishes were euthanized by exposure to the anesthetic benzocaine hydrochloride (400 mg l− 1) for 10 min following the cessation of gill movement. Brains, livers and samples of skeletal muscle were collected for analyses. For each species, three whole brains were frozen in liiquid nitrogen and stored at − 80 °C prior to molecular determinations.
The animal procedures were approved by responsible authorities (A13160603, from the Consejeria de Agua, Agricultura, Ganaderia y Pesca, Comunidad Autonoma de la Region de Murcia, Spain).
Coding sequence data
Our analysis comprised five clownfish species (A. ocellaris, A. clarkii, A. bicinctus, A. percula, A. sebae), C. viridis representing the non-symbiotic sister-taxon of the Amphiprion genus and nine more distantly related outgroup species (Stegastes partitus, Pundamilia nyererei, Maylandia zebra, Oryzias latipes, Xiphophorus maculatus, Poecilia formosa, Fundulus heteroclitus, Nothobranchius furzeri, Aphyosemion striatum). mRNA sequences of the outgroups were obtained from RefSeq along with their coding sequence annotation (Additional file 1: Table S12). For A. ocellaris, A. bicinctus and A. sebae, we downloaded read data from the short read archive (Bio projects PRJNA374650, PRJNA261388 and PRJNA285007, respectively). For A. clarkii, A. percula and C. viridis, we performed novel RNA-seq using samples from three individuals per species, the Qiagen RNeasy Mini Kit for purification, an Illumina HiSeq 2500 sequencing device (rapid run mode) and a paired end sequencing strategy with a read length of 150 base pairs. The number of resulting read pairs per sample was between 26 and 47 million reads (see Additional file 1: Table S13 for more details). The reads of the clownfishes and C. viridis were preprocessed using SeqPrep with minimum adapter length of five as well as a demanded minimum read length of 50. De novo transcriptome assemblies for these species were performed using FRAMA with Stegastes partitus as reference species . For the clownfishes and C. viridis the longest isoform was chosen to represent the gene. For the outgroups, in cases in which multiple isoforms per gene were annotated based on the reference, all of them were used in subsequent analyses. The assembly completeness of all examined species was estimated using BUSCO , was 90–100% (Additional file 1: Table S12).
Identification of positively selected genes
To scan on a genome-wide scale for genes under positive selection, we fed the coding sequences of the described species set into the PosiGene pipeline . Stegastes partitus was used as PosiGene’s anchor species. Orthology was determined by PosiGene via best bidirectional BLAST searches [80, 81] against Stegastes partitus. Within the analysis of positive selection PosiGene determined among others alignments of the orthologous genes using PRANK , reconstructed a phylogenetic tree of the examined species with DNAML from the PHYLIP package  and applied the branch-site test of positive selection  by using PAML . The branches of the last common ancestor of the clownfishes and C. viridis were tested separately for genes under positive selection (Additional file 1: Table S2). FDR < 0.05 was used as threshold for significance. Assembled sequence data, visualizations of alignments and positively selected sites were made available for download (see Additional file 2).
Furthermore, PSGs for A. striatum and C. porcellus were determined using the same sequences, species and PosiGene settings as described in the original killifish and mole-rat study, respectively [8, 23].
We determined enrichments for GO categories using Fisher’s exact test based on the R package GOstats (Table 2 and Additional file 1: Table S3). The resulting p-values were corrected using the Benjamini-Hochberg method . We used throughout the manuscript 0.1 as significance threshold. Applying implementations of Benjamini-Hochberg correction to raw p-values often result in same FDR values for multiple elements despite that they have different p-values, e.g. in row 3 to row 16 of Table 2 and Additional file 1: Table S3. The reason for this behavior is that most implementations of this method want to ensure that the sortation of elements by FDR is the same as by p-value. Therefore, they calculate the respective FDR values for the p-values p1 ≤ p2 ≤ … ≤ pn starting with FDRn and then set FDRi:= FDRi + 1, for 1 ≤ i < n in case that otherwise (according to the formula of the FDR) FDRi would be greater than FDRi + 1. By chance this happens, e.g., in the case of Table 2/ Additional file 1: Table S3 from FDR3 to FDR16 several times in a row.
Enrichment for mitochondrial biogenesis genes was tested using Fisher’s exact test (for the clownfish LCA, C. viridis, A. striatum and C. porcellus) and the union set of the genes in the following five mitochondrial related GO terms: GO:0000959, 0032543, 0045333, 0033108, 0070584 (Additional file 1: Table S5). The same GO terms were used in our previous study  to test for enrichment. In the case of the re-analyzation of the mole-rat study  with regard to mitochondrial biogenesis, Fisher’s exact test was not applicable since enrichment for multiple branches (in total 11) with different background gene sets had to be tested at the same time. Therefore, a resampling approach was used for p-value calculation analogous to the one that was used to test the three killifish branches at the same time for enrichment of the mitochondrial biogenesis gene set .
The relationship between aging and the results from the positive selection analysis on the clownfish LCA was tested based on a list of 41 gene ontology functions whose gene expressions were found to be associated with maximum lifespan residual in a study across 33 mammals . This means, that this study found the expression of the genes in those gene ontologies to be correlated or anti-correlated with the difference between observed maximum lifespan and expected maximum lifespan based on the weight of the respective species. Out of these 41 functions, 11 were tested during our enrichment analysis. We used Fisher’s method to compute the p-values of these 11 gene ontology processes from our enrichment analysis into a test statistic. The latter was tested against a null distribution that was empirically estimated from 100,000 random drawings of the p-values from enrichment analysis.
To identify genes that show signs of positive selection across multiple evolutionary branches on which the lifespan was altered considerably, we combined p-values from this study with those of two previous studies using Fisher’s method  (Additional file 1: Table S4, S5 and S6). In all three studies, PosiGene was used to determine p-values. The first study searched for genes under positive selection on 11 rodent branches in which the lifespan was presumably extended – most of them in the clade of the African mole-rat family that covers the longest-lived known rodents . The second study examined three branches of the Nothobranchius genus on which lifespan was presumably reduced  – the genus covers the shortest-lived vertebrate species that can be held in captivity . As significance threshold for meta-analysis p-values, 0.1 was chosen throughout the manuscript.
Protein homology modeling
Homology modelling of the clownfish GSTK1 was carried out with SWISS–MODEL (http://swissmodel.expasy.org; [89, 90] using the crystal structures of the dimeric apoform of the human mitochondrial GSTK1 (PDB 3rpp; ) and the substrate bound dimer of the rat enzyme PDB 1r4w; ). No further optimization was applied to the resulting models. Visualization, superimposition of the respective crystal structures and the models as well as rendering was carried out using CHIMERA . Figure 2c shows one subunit of the modelled dimer of the Clownfish GSTK1 model. Additional file 3: Figure S1 provides an overview over both subunits.
False discovery rate
Glutathione S-transferase kappa 1
Last common ancestor
Mitochondrial transcription termination factor 1
Positively selected gene
Vrtilek M, Zak J, Polacik M, Blazek R, Reichard M. Longitudinal demographic study of wild populations of African annual killifish. Sci Rep. 2018;8(1):4774.
Nielsen J, Hedeholm RB, Heinemeier J, Bushnell PG, Christiansen JS, Olsen J, Ramsey CB, Brill RW, Simon M, Steffensen KF, et al. Eye lens radiocarbon reveals centuries of longevity in the Greenland shark (Somniosus microcephalus). Science. 2016;353(6300):702–4.
Finch TM, Zhao N, Korkin D, Frederick KH, Eggert LS. Evidence of positive selection in mitochondrial complexes I and V of the African elephant. PLoS One. 2014;9(4):e92587.
Seim I, Fang X, Xiong Z, Lobanov AV, Huang Z, Ma S, Feng Y, Turanov AA, Zhu Y, Lenz TL, et al. Genome analysis reveals insights into physiology and longevity of the Brandt's bat Myotis brandtii. Nat Commun. 2013;4:2212.
Keane M, Semeiks J, Webb AE, Li YI, Quesada V, Craig T, Madsen LB, van Dam S, Brawand D, Marques PI, et al. Insights into the evolution of longevity from the bowhead whale genome. Cell Rep. 2015;10(1):112–22.
Kim EB, Fang X, Fushan AA, Huang Z, Lobanov AV, Han L, Marino SM, Sun X, Turanov AA, Yang P, et al. Genome sequencing reveals insights into physiology and longevity of the naked mole rat. Nature. 2011;479(7372):223–7.
Fang X, Nevo E, Han L, Levanon EY, Zhao J, Avivi A, Larkin D, Jiang X, Feranchuk S, Zhu Y, et al. Genome-wide adaptive complexes to underground stresses in blind mole rats Spalax. Nat Commun. 2014;5:3966.
Sahm A, Bens M, Szafranski K, Holtze S, Groth M, Gorlach M, Calkhoven C, Muller C, Schwab M, Kraus J, et al. Long-lived rodents reveal signatures of positive selection in genes associated with lifespan. PLoS Genet. 2018;14(3):e1007272.
Li Y, de Magalhaes JP. Accelerated protein evolution analysis reveals genes and pathways associated with the evolution of mammalian longevity. Age (Dordr). 2013;35(2):301–14.
Muntane G, Farre X, Rodriguez JA, Pegueroles C, Hughes DA, de Magalhaes JP, Gabaldon T, Navarro A. Biological processes modulating longevity across Primates: a phylogenetic genome-phenome analysis. Mol Biol Evol. 2018;35(8):1990–2004.
Aledo JC, Li Y, de Magalhaes JP, Ruiz-Camacho M, Perez-Claros JA. Mitochondrially encoded methionine is inversely related to longevity in mammals. Aging Cell. 2011;10(2):198–207.
Fushan AA, Turanov AA, Lee SG, Kim EB, Lobanov AV, Yim SH, Buffenstein R, Lee SR, Chang KT, Rhee H, et al. Gene expression defines natural changes in mammalian lifespan. Aging Cell. 2015;14(3):352–65.
Medawar PB. An unsolved problem of biology: Printed lecture: University College London; 1952.
Williams GC. Pleiotropy, Natural Selection, and the Evolution of Senescence. Evolution. 1957;11:398–411.
Kirkwood TB. Evolution of ageing. Nature. 1977;270(5635):301–4.
Hossie TJ, Hassall C, Knee W, Sherratt TN. Species with a chemical defence, but not chemical offence, live longer. J Evol Biol. 2013;26(7):1598–602.
Blanco MA, Sherman PW. Maximum longevities of chemically protected and non-protected fishes, reptiles, and amphibians support evolutionary hypotheses of aging. Mech Ageing Dev. 2005;126(6–7):794–803.
Shattuck MR, Williams SA. Arboreality has allowed for the evolution of increased longevity in mammals. Proc Natl Acad Sci U S A. 2010;107(10):4635–9.
Voituron Y, de Fraipont M, Issartel J, Guillaume O, Clobert J. Extreme lifespan of the human fish (Proteus anguinus): a challenge for ageing mechanisms. Biol Lett. 2011;7(1):105–7.
Tozzini ET, Dorn A, Ng'oma E, Polacik M, Blazek R, Reichwald K, Petzold A, Watters B, Reichard M, Cellerino A. Parallel evolution of senescence in annual fishes in response to extrinsic mortality. BMC Evol Biol. 2013;13:77.
Blazek R, Polacik M, Kacer P, Cellerino A, Rezucha R, Methling C, Tomasek O, Syslova K, Terzibasi Tozzini E, Albrecht T, et al. Repeated intraspecific divergence in life span and aging of African annual fishes along an aridity gradient. Evolution. 2017;71(2):386–402.
Cellerino A, Valenzano DR, Reichard M. From the bush to the bench: the annual Nothobranchius fishes as a new model system in biology. Biol Rev Camb Philos Soc. 2016;91(2):511–33.
Sahm A, Bens M, Platzer M, Cellerino A. Parallel evolution of genes controlling mitonuclear balance in short-lived annual fishes. Aging Cell. 2017;16(3):488–96.
Litsios G, Sims CA, Wuest RO, Pearman PB, Zimmermann NE, Salamin N. Mutualism with sea anemones triggered the adaptive radiation of clownfishes. BMC Evol Biol. 2012;12:212.
Fautin DG, Allen GR. Field guide to anemonefishes and their Host Sea anemones. Perth: Western Australia Museum; 1992.
Elliott JK, Mariscal RN. Ontogenetic and interspecific variation in the protection of anemonefishes from sea anemones. J Exp Mar Biol Ecol. 1997;208(1):57–72.
Ollerton J, McCollin D, Fautin DG, Allen GR. Finding NEMO: nestedness engendered by mutualistic organization in anemonefish and their hosts. Proc Biol Sci. 2007;274(1609):591–8.
Buston PM. Mortality is associated with social rank in the clown anemonefish (Amphiprion percula). Mar Biol. 2003;143:811–5.
Elliott JK, Elliott JM, Mariscal RN. Host selection, location, and association behaviors of anemonefishes in field settlement experiments. Mar Biol. 1995;122(3):377–89.
Mariscal RN. The nature of the symbiosis between indo-Pacific anemone fishes and sea anemones. Mar Biol. 1970;6(1):58–65.
Buston PM, García MB. An extraordinary life span estimate for the clown anemonefish Amphiprion percula. J Fish Biol. 2007;70(6):1710–9.
Aldenhoven JM. Local variation in mortality rates and life-expectancy estimates of the coral reef fish Centropyge bicolor (Pisces: Pomacanthidae). Mar Biol. 1986;92:237–44.
Eckert GJ. Estimates of adult and juvenile mortality for labrid fishes at one tree reef, great barrier reef. Mar Biol. 1987;95:161–71.
Munro JL, Williams DM. Assessment and management of coral reef fisheries: biological, environmental, and socioeconomic aspects. In: Proceedings of the 5th International Coral Reef Congress, vol. 4; 1985. p. 545–81.
Maison KA, Graham KS. Status review report: orange clownfish (Amphiprion percula). Washington, D.C.: Pacific Islands Fisheries Science Center, National Marine Fisheries Service, National Oceanic and Atmospheric Administration, US Department of Commerce; 2016.
Salles OC, Maynard JA, Joannides M, Barbu CM, Saenz-Agudelo P, Almany GR, Berumen ML, Thorrold SR, Jones GP, Planes S. Coral reef fish populations can persist without immigration. Proc Biol Sci. 2015;282(1819):20151311.
Olivotto I, Planas M, Simões N, Holt GJ, Avella MA, Calado R. Advances in breeding and rearing marine ornamentals. J World Aquacult Soc. 2011;42(2):135–66.
Sahm A, Platzer M, Cellerino A. Outgroups and positive selection: the Nothobranchius furzeri case. Trends Genet. 2016;32(9):523–5.
Garcia-Herrera N, Ferse SCA, Kunzmann A, Genin A. Mutualistic damselfish induce higher photosynthetic rates in their host coral. J Exp Biol. 2017;220(10):1803–11.
Holbrook SJ, Brooks AJ, Schmitt RJ, Stewart HL. Effects of sheltering fish on growth of their host corals. Mar Biol. 2008;155(5):521–30.
Sweatman H. The timing of settlement by larval Dascyllus aruanus: some consequences for larval habitat selection. In: Proc 5th Int Coral Reef Conf 5; 1985. p. 367–72.
Liberman T, Genin A, Loya Y. Effects on growth and reproduction of the coral Stylophora pistillata by the mutualistic damselfish Dascyllus marginatus. Mar Biol. 1995;121:741–6.
Ben-Tzv IO, Abelson A, Polak O, Kiflawi M. Habitat selection and the colonization of new territories by Chromis viridis. Journal of Fish Biology. 2008;73(4):1005–18.
Lecchini D, Nakamura Y, Grignon J, Tsuchiya M. Evidence of density-independent mortality in a settling coral reef damselfish, Chromis viridis. Ichthyological Research. 2006;53(3):298–300.
Hixon MA, Carr MH. Synergistic predation, density dependence, and population regulation in marine fish. Science. 1997;277(5328):946–9.
Wantiez L, Thollot P. Settlement, post-settlement mortality and growth of the damselfish Chromis fumea (Pisces: Pomacentridae) on two artificial reefs in New Caledonia (south-West Pacific Ocean). J Mar Biol Assoc U K. 2001;80(6):1111–8.
Dulčić J, Kraljević M. Age, growth and mortality of damselfish (Chromis chromis L.) in the eastern middle Adriatic. Fish Res. 1995;22(3):255–64.
McElwee JJ, Schuster E, Blanc E, Piper MD, Thomas JH, Patel DS, Selman C, Withers DJ, Thornton JM, Partridge L, et al. Evolutionary conservation of regulated longevity assurance mechanisms. Genome Biol. 2007;8(7):R132.
Plank M, Wuttke D, van Dam S, Clarke SA, de Magalhaes JP. A meta-analysis of caloric restriction gene expression profiles to infer common signatures and regulatory mechanisms. Mol BioSyst. 2012;8(4):1339–49.
Amador-Noguez D, Dean A, Huang W, Setchell K, Moore D, Darlington G. Alterations in xenobiotic metabolism in the long-lived little mice. Aging Cell. 2007;6(4):453–70.
Steinbaugh MJ, Sun LY, Bartke A, Miller RA. Activation of genes involved in xenobiotic metabolism is a shared signature of mouse models with extended lifespan. Am J Physiol Endocrinol Metab. 2012;303(4):E488–95.
Houtkooper RH, Mouchiroud L, Ryu D, Moullan N, Katsyuba E, Knott G, Williams RW, Auwerx J. Mitonuclear protein imbalance as a conserved longevity mechanism. Nature. 2013;497(7450):451–7.
Baumgart M, Priebe S, Groth M, Hartmann N, Menzel U, Pandolfini L, Koch P, Felder M, Ristow M, Englert C, et al. Longitudinal RNA-Seq analysis of vertebrate aging identifies mitochondrial complex I as a small-molecule-sensitive modifier of lifespan. Cell Syst. 2016;2(2):122–32.
Wang B, Peng Y, Zhang T, Ding J. Crystal structures and kinetic studies of human kappa class glutathione transferase provide insights into the catalytic mechanism. Biochem J. 2011;439(2):215–25.
Ladner JE, Parsons JF, Rife CL, Gilliland GL, Armstrong RN. Parallel evolutionary pathways for glutathione transferases: structure and mechanism of the mitochondrial class kappa enzyme rGSTK1-1. Biochemistry. 2004;43(2):352–61.
Valenzano DR, Benayoun BA, Singh PP, Zhang E, Etter PD, Hu CK, Clement-Ziza M, Willemsen D, Cui R, Harel I, et al. The African turquoise killifish genome provides insights into evolution and genetic architecture of lifespan. Cell. 2015;163(6):1539–54.
Sahm A, Cellerino A. (Anti-)parallel evolution of lifespan. Aging (Albany NY). 2017;9(10):2018–9.
Jemielity S, Chapuisat M, Parker JD, Keller L. Long live the queen: studying aging in social insects. Age (Dordr). 2005;27(3):241–8.
Roux J, Privman E, Moretti S, Daub JT, Robinson-Rechavi M, Keller L. Patterns of positive selection in seven ant genomes. Mol Biol Evol. 2014;31(7):1661–85.
Shen YY, Liang L, Zhu ZH, Zhou WP, Irwin DM, Zhang YP. Adaptive evolution of energy metabolism genes and the origin of flight in bats. Proc Natl Acad Sci U S A. 2010;107(19):8666–71.
Tacutu R, Craig T, Budovsky A, Wuttke D, Lehmann G, Taranukha D, Costa J, Fraifeld VE, de Magalhaes JP. Human ageing genomic resources: integrated databases and tools for the biology and genetics of ageing. Nucleic Acids Res. 2013;41(Database issue):D1027–33.
Sahm A, Bens M, Platzer M, Cellerino A. Parallel evolution of genes controlling mitonuclear balance in short-lived annual fishes. Aging Cell. 2017;16(3):488-496.
Antonicka H, Shoubridge EA. Mitochondrial RNA granules are centers for posttranscriptional RNA processing and ribosome biogenesis. Cell Rep. 2015;10(6):920-932.
Eskelinen EL. Roles of LAMP-1 and LAMP-2 in lysosome biogenesis and autophagy. Mol Asp Med. 2006;27(5–6):495–502.
Berditchevski F, Odintsova E. Tetraspanins as regulators of protein trafficking. Traffic. 2007;8(2):89–96.
Sordat I, Decraene C, Silvestre T, Petermann O, Auffray C, Pietu G, Sordat B. Complementary DNA arrays identify CD63 tetraspanin and alpha3 integrin chain as differentially expressed in low and high metastatic human colon carcinoma cells. Lab Investig. 2002;82(12):1715–24.
Sauer G, Kurzeder C, Grundmann R, Kreienberg R, Zeillinger R, Deissler H. Expression of tetraspanin adaptor proteins below defined threshold values is associated with in vitro invasiveness of mammary carcinoma cells. Oncol Rep. 2003;10(2):405–10.
Zhijun X, Shulan Z, Zhuo Z. Expression and significance of the protein and mRNA of metastasis suppressor gene ME491/CD63 and integrin alpha5 in ovarian cancer tissues. Eur J Gynaecol Oncol. 2007;28(3):179–83.
Kwon MS, Shin SH, Yim SH, Lee KY, Kang HM, Kim TM, Chung YJ. CD63 as a biomarker for predicting the clinical outcomes in adenocarcinoma of lung. Lung Cancer. 2007;57(1):46–53.
Lai X, Gu Q, Zhou X, Feng W, Lin X, He Y, Cao J, Liu P, Zhang H, Zheng X. Decreased expression of CD63 tetraspanin protein predicts elevated malignant potential in human esophageal cancer. Oncol Lett. 2017;13(6):4245–51.
Zhang C, Cuervo AM. Restoration of chaperone-mediated autophagy in aging liver improves cellular maintenance and hepatic function. Nat Med. 2008;14(9):959–65.
Gao F, Fang Q, Zhang R, Lu J, Lu H, Wang C, Ma X, Xu J, Jia W, Xiang K. Polymorphism of DsbA-L gene associates with insulin secretion and body fat distribution in Chinese population. Endocr J. 2009;56(3):487–94.
Sharma M, Gupta S, Singh K, Mehndiratta M, Gautam A, Kalra OP, Shukla R, Gambhir JK. Association of glutathione-S-transferase with patients of type 2 diabetes mellitus with and without nephropathy. Diabetes Metab Syndr. 2016;10(4):194–7.
de Magalhaes JP, Curado J, Church GM. Meta-analysis of age-related gene expression profiles identifies common signatures of aging. Bioinformatics. 2009;25(7):875–81.
Aramillo Irizar P, Schauble S, Esser D, Groth M, Frahm C, Priebe S, Baumgart M, Hartmann N, Marthandan S, Menzel U, et al. Transcriptomic alterations during ageing reflect the shift from cancer to degenerative diseases in the elderly. Nat Commun. 2018;9(1):327.
Kurz T, Terman A, Gustafsson B, Brunk UT. Lysosomes and oxidative stress in aging and apoptosis. Biochim Biophys Acta. 2008;1780(11):1291–303.
Bens M, Sahm A, Groth M, Jahn N, Morhart M, Holtze S, Hildebrandt TB, Platzer M, Szafranski K. FRAMA: from RNA-seq data to annotated mRNA assemblies. BMC Genomics. 2016;17:54.
Simao FA, Waterhouse RM, Ioannidis P, Kriventseva EV, Zdobnov EM. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics. 2015;31(19):3210–2.
Sahm A, Bens M, Platzer M, Szafranski K. PosiGene: automated and easy-to-use pipeline for genome-wide detection of positively selected genes. Nucleic Acids Res. 2017;45(11):e100.
Overbeek R, Fonstein M, D'Souza M, Pusch GD, Maltsev N. The use of gene clusters to infer functional coupling. Proc Natl Acad Sci U S A. 1999;96(6):2896–901.
Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, Madden TL. BLAST+: architecture and applications. BMC Bioinformatics. 2009;10:421.
Loytynoja A, Goldman N. Phylogeny-aware gap placement prevents errors in sequence alignment and evolutionary analysis. Science. 2008;320(5883):1632–5.
Felsenstein J. PHYLIP (Phylogeny Inference Package) version 3.6. Distributed by the author. Seattle: Department of Genome Sciences, University of Washington; 2005.
Zhang J, Nielsen R, Yang Z. Evaluation of an improved branch-site likelihood method for detecting positive selection at the molecular level. Mol Biol Evol. 2005;22(12):2472–9.
Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007;24(8):1586–91.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol. 1995;57(1):289–300.
Fisher RA. Statistical methods for research workers; 1932.
Valdesalici S, Cellerino A. Extremely short lifespan in the annual fish Nothobranchius furzeri. Proc Biol Sci. 2003;270(Suppl 2):S189–91.
Arnold K, Bordoli L, Kopp J, Schwede T. The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling. Bioinformatics. 2006;22(2):195–201.
Biasini M, Bienert S, Waterhouse A, Arnold K, Studer G, Schmidt T, Kiefer F, Gallo Cassarino T, Bertoni M, Bordoli L, et al. SWISS-MODEL: modelling protein tertiary and quaternary structure using evolutionary information. Nucleic Acids Res. 2014;42(Web Server issue):W252–8.
Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, Ferrin TE. UCSF Chimera--a visualization system for exploratory research and analysis. J Comput Chem. 2004;25(13):1605–12.
Moyer JT. Longevity of the anemonefish Amphiprion clarkii at Miyake-Jima, Japan with Notes on Four Other Species. Copeia. 1986;1986:135–9.
Froese R, Pauly D: FishBase. World Wide Web electronic publication 2018, www.fishbase.org, (06/2018).
We thank Emilio Cortéz Melendrenas of the Acuario del la Universidad de Murcia for help in fish husbandry and the following individuals at zoos or aquariums for responding to our questionnaire: Vicky Béduneau (Océarium du Croisic, Le Croisic, France), Nicolas Hirel (Aquarium Mare Nostrum, Montpellier, France), Thomas Ziegler (Koelner Zoo, Köln, Germany), Nikolaj Meyer (Skansen-Akvariet, Stockholm, Sweden), Markus Dernjatin (Sea Life Helsinki, Helsinki, Finland). Personal communications on clownfish longevity were received also from Prof. Ike Olivotto (University of Ancona) and Prof. Hellen Thaler (Innsbruck). We thank Cornelia Luge, Ivonne Görlich and Marco Groth (CF DNA sequencing at Leibniz Institute on Aging - Fritz Lipmann Institute) for conducting Illumina sequencing. We thank Debra Weih for proofreading the manuscript. We thank Matthias Platzer and Steve Hoffmann for helpful discussions and support.
The project was partially supported by internal grant of Scuola Normale Superiore SNS17_B_CELLERINO and SNS_RB_CELLERINO to AC that provided basic support to the research activity that were used to support collection and analysis of the data.
Availability of data and materials
The datasets generated and/or analysed during the current study are available in the European Nucleotide Archive as study PRJEB27750 (https://www.ebi.ac.uk/ena/data/view/PRJEB27750).
Ethics approval and consent to participate
The animal procedures were approved by responsible authorities (A13160603, from the Consejeria de Agua, Agricultura, Ganaderia y Pesca, Comunidad Autonoma de la Region de Murcia, Spain).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Clownfish lifespan questionnaire results. Table S2. PosiGene results for positively selected genes on the phylogenetic branch representing the last common ancestor of the clownfishes (genus Amphiprion). Table S3. Enrichment test results of the biological gene ontology processes enriched for positively selected genes. Table S4. Meta analysis using Fisher’s method of positive selection across three analyses of phylogenetic branches on which lifespan changed considerably. Table S5. Meta analysis using Fisher’s method of positive selection across three phylogenetic branches of the Nothobranchius genus on which the lifespan was reduced considerably. Table S6. Meta analysis using Fisher’s method of positive selection across eleven phylogenetic rodent branches on which the lifespan was reduced considerably. Table S7. Genes that were regarded as mitochondrial biogenesis related from five gene ontology terms. Table S8. PosiGene results for positively selected genes on the phylogenetic branch representing Chromis viridis. Table S9. PosiGene results for positively selected genes on the phylogenetic branch representing Cavia porcellus. Table S10. PosiGene results for positively selected genes on the phylogenetic branch representing Aphysemion striatum. Table S11. Meta analysis using Fisher’s method of positive selection across three analyses of phylogenetic branches on which the lifespan did not change considerably. Table S12. Assembly and sequence statistics. Table S13. Samples that were sequenced to create genome/transcriptome assemblies. (XLSX 6532 kb)
Supplement data. (DOCX 12 kb)
Figure S1. Homology modelling of Clownfish GSTK1. Ribbon representation of the model dimer for the clownfish enzyme as derived from SWISS-MODEL in grey, superimposed onto the dimeric structure of the substrate bound rat GSTK1 (PDB 1r4w; ) used as template in light green. The pairwise root mean square deviation for the Cα positions between the model and 1r4w amounts to 0.52 Å as determined with the CHIMERA Matchmaker tool. The GSH substrate in the rat enzyme structure is depicted in light purple. (PNG 1250 kb)