Research Paper Volume 9, Issue 2 pp 573—582

Adaptive sequence convergence of the tumor suppressor ADAMTS9 between small-bodied mammals displaying exceptional longevity

Matthew J. Lambert1, , Christine V. Portfors1, ,

  • 1 School of Biological Sciences, Washington State University, Vancouver, WA 98686, USA

Received: January 6, 2017       Accepted: February 11, 2017       Published: February 26, 2017      

https://doi.org/10.18632/aging.101180
How to Cite

Abstract

Maximum lifespan varies by two orders of magnitude across mammals. How such divergent lifespans have evolved remains an open question, with ramifications that may potentially lead to therapies for age-related diseases in humans. Several species of microbats as well as the naked mole-rat live much longer than expected given their small sizes, show reduced susceptibility to neoplasia, and largely remain healthy and reproductively capable throughout the majority of their extended lifespans. The convergent evolution of extreme longevity in these two groups allows for the opportunity to identify potentially important aging related genes that have undergone adaptive sequence convergence in these long-lived, yet small-bodied species. Here, we have tested 4,628 genes for evidence of convergence between the microbats and naked mole-rat. We find a strong signal of adaptive sequence convergence in the gene A disintegrin-like and metalloprotease with thrombospondin type 1 motifs 9 (ADAMTS9). We also provide evidence that the shared substitutions were driven by selection. Intriguingly, ADAMTS9 is a known inhibitor of the mTor pathway and has been implicated in several aging related processes.

Introduction

For several decades, it has been well recognized that there is strong correlation between lifespan and body mass [1,2], with larger species typically living longer than smaller species. There are, however, several species that violate this general rule, living much longer than expected given their small size and high metabolic rates. Of particular interest are the microbats, several species of which demonstrate longer maximum lifespans than any other mammals when controlling for body size [3]. In addition to their exceptional longevity, microbats appear to be resistant to neoplasia [4,5] and remain healthy and reproductively capable throughout the majority of their lives [6].

Much like the microbats, Heterocephalus glaber (naked mole-rat) lives approximately three times longer than expected given its small size [7], is remarkably resistant to neoplasia [810] and displays no symptoms of aging well into its second decade [11].

Although once thought to be rare, there have been numerous recent studies demonstrating adaptive sequence convergence between a variety of species displaying convergent traits. These studies have highlighted genes that have been repeatedly targeted during the evolution of a given trait. For example, the evolution of echolocation in bats and toothed whales appears to be driven, in part, by common mutations in the genes Prestin [12,13] and Cdh23 [14,15]; cardiac glycoside toxin resistance in numerous disparate vertebrates and invertebrates can be explained by identical amino acid substitutions in the enzyme Na+/K+-ATPase [16]; and the repeated convergent evolution of the sodium ion transporter NaV1.7 in hibernating mammals and mole-rats is believed to impart insensitivity to the accumulation of CO2 [17]. These studies, and many others, demonstrate that common selective pressures can drive common mutations in relevant genes.

The evolution of extreme longevity in microbats and the naked mole-rat is likely attributable to a lack of extrinsic sources of mortality in these species. Bats, being nocturnal and capable of flight, generally contend with few predators. Likewise, the naked mole-rat lives in subterranean burrows where the risk of predation is low. Several theories of aging suggest that a lack of extrinsic sources of mortality will result in selection for longer lifespan [8,18]. For example, according to the antagonistic pleiotropy (AP) theory of aging, a mutation can be beneficial during development, but have late-onset deleterious effects [19]. AP is expected to be more prevalent in species with high levels of extrinsic mortality since most individuals are unlikely to survive long after reaching sexual maturity, therefore there will be little pressure to select against the deleterious effects that manifest later in life. Also, the disposable soma theory of aging suggests that there exists a trade-off between growth/development and repair/maintenance [20]. In species that contend with many predators, it should be beneficial to allocate resources to grow and develop as quickly as possible rather than to invest in repair and maintenance since longevity is already unlikely. According to both theories, for species that contend with numerous extrinsic sources of mortality, the decline in fitness due to aging is minimal, so selection is inefficient at promoting mutations that increase longevity. However, for species that exist in relatively safe niches, like microbats and the naked mole-rat, the strength of selection to delay senescence will be much stronger, as individuals that live longer will have higher lifetime reproductive fitness. We hypothesize that the pressure to delay senescence shared by microbats and naked mole-rat may have led to convergent sequence evolution in key longevity promoting genes. The identification of genes that have undergone convergent evolution in these long-lived species would provide a better understanding of the genetics of longevity and could potentially identify therapeutic targets for cancer and other age-related illnesses.

Here we tested for adaptive convergent sequence evolution between microbats and the naked mole-rat in almost 5,000 genes conserved across a wide-range of mammals. We found that A disintegrin-like and metalloprotease with thrombospondin type 1 motifs 9 (ADAMTS9) displays numerous convergent substitutions between the long-lived species that were likely driven by positive selection.

Results

The goal of this study was to determine if any genes display evidence of adaptive sequence convergence between long-lived microbats and the naked mole-rat. We tested for convergence in 4,628 1:1 ortholog groups between 3 microbat species (Myotis lucifugus, Myotis brandtii and Eptesicus fuscus), and the naked mole-rat. Additionally, as outgroup species, we included sequences from 11 small-bodied rodents and laurasiatherian mammals with at least acceptable longevity data curated in the AnAge database [21]. The relationship between the 15 species is shown in the phylogeny in Figure 1a. Importantly, none of these additional species demonstrate exceptional longevity (Fig. 1b). The requirement that sister species to microbats and naked mole-rat be small-bodied and not long-lived allowed us to control for body size and phylogeny while giving us the opportunity to test for convergence in the long-lived species as compared to normal aging species.

Small-bodied mammals that display exceptional longevity. (A) The accepted tree topology for the 15 species included in this study. (B) Maximum longevity residual (left axis, black boxes) and maximum longevity (right axis, grey circles) of the 15 species. Maximum longevity residual (tmax) is the percentage of the expected maximum longevity given adult body size (M), derived from the mammalian allometric equation: tmax = 4.88M0.153 [30].

Figure 1. Small-bodied mammals that display exceptional longevity. (A) The accepted tree topology for the 15 species included in this study. (B) Maximum longevity residual (left axis, black boxes) and maximum longevity (right axis, grey circles) of the 15 species. Maximum longevity residual (tmax) is the percentage of the expected maximum longevity given adult body size (M), derived from the mammalian allometric equation: tmax = 4.88M0.153 [30].

Initially, we constructed maximum likelihood gene trees based on nucleotide alignments for each of the 4,628 ortholog groups. These trees were then parsed programmatically, searching for instances in which the long-lived mammals (naked mole-rat, Myotis davidii, Myotis lucifugus and Eptesicus fuscus) formed a monophyletic group, to the exclusion of their true closest relatives (rodents and carnivores/insectivores for the naked mole-rat and microbats, respectively). We found that only one gene, ADAMTS9, violated the known species tree in such a manner (Fig. 2a), suggesting sequence convergence between the naked mole-rat and microbats occurred during the evolution of this gene. To determine which sites were responsible for the observed grouping of the long-lived species, we used a maximum likelihood approach to reconstruct all ancestral amino acid sequences at interior nodes of the accepted phylogeny. These ancestral sequences were used to map and count all pairwise convergent substitutions (mutations from the ancestral state that lead to identical amino acids) shared between the microbats and naked mole-rat. Although, the microbats and naked mole-rat share 21 convergent substitutions in this gene (Table 1), which is far more than any other instance of sequence convergence described to date, it may be possible that this represents a normal level of background convergence if for example ADAMTS9 is rapidly evolving. However, Castoe et al. [22] developed a method to generate an empirical null distribution for the expected level of convergence for a given gene, provided a sufficient number of taxa are included. Briefly, the number of divergent (mutations that lead to different amino acids at homologous sites in the species being compared) and convergent substitutions are estimated for every pair of branches in a phylogeny. It has previously been demonstrated that the number of divergent substitutions reliably predict the number of convergent substitutions, therefore excess convergence between any two branches would appear as an outlier with a higher convergent to divergent substitution ratio than all other pairwise comparisons between branches [22]. This method allowed us to distinguish between rapid evolution and potential adaptive convergence. To test this, we used the program Grand-conv [23,24], which estimates the posterior numbers of convergent and divergent substitutions shared between all pairs of branches in the given phylogeny. We found that the convergence between the naked mole-rat and microbats exceeded expectations given the phylogeny (Fig. 2b), again suggesting sequence convergence between the long-lived species. In addition, using the methods of Zhang and Kumar [25], we found that the level of convergence between microbats and naked mole-rat exceeds random expectations (p < .001), indicating a significant amount of convergence.

Evidence for ADAMTS9 sequence convergence in long-lived mammals. (A) Maximum-likelihood gene tree constructed using the best fitting nucleotide substitution model (GTR+G) demonstrates a monophyletic grouping of the long-lived mammals. (B) Plot of the numbers of divergent versus convergent substitutions occurring between all independent pairs of branches of the species tree. The microbat-naked mole-rat comparison has the highest level of convergence as well as the largest distance from the trendline.

Figure 2. Evidence for ADAMTS9 sequence convergence in long-lived mammals. (A) Maximum-likelihood gene tree constructed using the best fitting nucleotide substitution model (GTR+G) demonstrates a monophyletic grouping of the long-lived mammals. (B) Plot of the numbers of divergent versus convergent substitutions occurring between all independent pairs of branches of the species tree. The microbat-naked mole-rat comparison has the highest level of convergence as well as the largest distance from the trendline.

Table 1. ADAMTS9 convergent sites. The numbers at the top of each column indicate the amino acid position in the alignment.

54236244324362512558565775838884890111013141319148914971597167517651933
Naked mole-ratVRRVPSDQAKVRQRRRERAST
Eptesicus fuscusVRRVPSDQAKVRQRRRERAST
Myotis brandtiiVRRVPSDQAKVRQRR--RAST
Myotis lucifugusVVQVPSDQAKVRQRRRERAST
RatIKKISP--TSIREQTKDQTNL
MouseIKKIYPHRTSIKEQTKDQTNK
Prairie voleIKKISPHRTSIKEQTKDQTNK
HamsterIKKISPHRTNIKEQRKDQTNK
Star-nosed moleIKKISSPRTNVKEQSKDKTKK
Shrew-KKISXLRTNVKEQSXXKTLK
FerretIKKISP-RTNIKEQSKDKTEK
Tree shrewVKKIXPXXTX-KEQSKEKXXK
Hedgehog-SKISYHRXNVKEQDKDKTNK
Opossum-KRVAPHRTNISARRKDKTKK
PikaISKVSPHRTNIKEQMKDKAQK

Notably, the addition of the moderately long-lived rodent, Cavia porcellus (guinea pig), and megabats, Pteropus alecto and Pteropus vampyrus, disrupted the monophyletic grouping of the naked mole-rat and microbats in maximum-likelihood gene trees based on amino acid alignments (Not shown). However, even with the inclusion of these moderately long-lived mammals we still found that the naked mole-rat and microbats share 16 convergent substitutions (Supplemental Table1), which remains significant using the method of Zhang and Kumar [25]. Furthermore, estimates of posterior numbers of divergent versus convergent changes still suggest sequence convergence between the naked mole-rat and microbats even when the megabats and guinea pig were included (Supplemental Figure 1).

Lastly, to test if the observed sequence convergence may be due to positive selection, we employed a unique multistep approach. First, all sites with significant evidence for positive selection along the ancestral microbat branch (indicated by a star in Fig. 1a) were identified with the program TreeSAAP [26], producing a list of 77 sites. These 77 sites were then extracted from the original alignment and used to create a new alignment that consisted only of the sites with evidence of positive selection in microbats. Next, we calculated site-wise log-likelihood values for the convergent gene tree (Fig. 2a) and compared these to site-wise log likelihood values for the accepted species tree (Fig. 1a), using the positively selected sites alignment. If the observed convergence was due to chance, we would expect the positively selected sites to support the accepted species topology. However, we found that the opposite is true, the sites with evidence of positive selection strongly favored the convergent topology (Fig. 3). Indeed, an approximately unbiased test suggests that the convergent tree was much more likely than the accepted species tree (P = .886 and P = .114, respectively), suggesting that the convergence between the long-lived species was driven by selection.

The relationship between sequence convergence and positive selection. Relative support for convergent tree and species tree topologies for all positively selected sites along the ADAMTS9 gene sequence. The values are the difference between site-wise log likelihood scores for the species tree and the convergent tree. Negative values indicate more support for the convergent tree.

Figure 3. The relationship between sequence convergence and positive selection. Relative support for convergent tree and species tree topologies for all positively selected sites along the ADAMTS9 gene sequence. The values are the difference between site-wise log likelihood scores for the species tree and the convergent tree. Negative values indicate more support for the convergent tree.

Discussion

We tested 4,628 genes for evidence of adaptive sequence convergence between long-lived, small-bodied mammals. We found evidence for an enrichment of convergent substitutions between the microbats and the naked mole-rat in the gene ADAMTS9. ADAMTS9 is the most widely conserved member of the ADAMTS family and has recently been reported to be a novel tumor suppressor that is downregulated in several varieties of human cancer [27]. Intriguingly, ADAMTS9 inhibits tumor growth by blocking the mTOR pathway [28], which has long been known to be associated with aging [29]. In addition to its role in tumor suppression, ADAMTS9 has also been implicated in several age-related conditions including arthritis [30], type 2 diabetes [31,32], age-related macular degeneration [33,34] and menopause [35]. Furthermore, in C. elegans the loss of GON-1, the roundworm homolog of ADAMTS9, alters lifespan and promotes dauer formation [36]. These effects are likely due to modified insulin and insulin-like ortholog secretion and altered insulin/IGF-1 signaling, which is also known to contribute to aging [37].

Although, it may be possible that the observed convergent changes shared by microbats and the naked mole-rat may be the product of some non-adaptive force rather than selection for increased longevity, several lines of evidence suggest otherwise. First, the convergent substitutions are distributed along the length of the coding sequence, eliminating gene conversion or alternate exon usage as possible causes. Second, the convergent topology was strongly favored when only sites with evidence of positive selection occurring on the long-lived microbat branch were considered, suggesting that the convergence was indeed driven by selection. Finally, ADAMTS9 has previously been implicated in several aging processes and age-related diseases [27,3037], supporting the hypothesis that modulation of ADAMTS9 function alters lifespan. Together, this evidence suggests that ADAMTS9 has been repeatedly targeted by selection for increased longevity in microbats and the naked mole-rat.

As yet, we can provide no explanation for how the convergent evolution of ADAMTS9 has promoted longevity in microbats and the naked mole-rat. However, it is rather intriguing that ADAMTS9 is a known inhibitor of the mTor pathway. The hyperfunction theory of aging posits that aging is due to the prolonged activity of growth-promoting pathways, such as mTor [38]. According to this theory, aging itself is not a program, but rather aging is a harmful continuation of developmental programs after growth has ended [39]. This hypothesis is supported by several lines of evidence that have identified mTor activity as a driver of senescence [29,4043]. Therefore, it is a distinct possibility that the convergent evolution of ADAMTS9 in long-lived species may have altered this gene’s function such that it now serves to slow the effects of mTor, preventing the declines normally associated with aging. Indeed, the author of the hyperfunction theory has recently speculated that the extraordinary longevity of microbats and the naked mole-rat may be due to selection for a “decelerator of mTor” [44].

Given the role of ADAMTS9 in such a wide range of age-related conditions, its direct effect on lifespan in C. elegans and its ability to modulate mTOR and insulin/IGF-1 signaling, it is likely that the convergent evolution of ADAMTS9 may be, in part, responsible for the exceptional longevity and resistance to neoplasia found in microbats and the naked mole-rat.

Materials and Methods

Ensembl Compara [45] was used to identify all 1:1 orthologs from the species Myotis lucifugus, Pteropus vampyrus, Erinaceus Europaeus, Mustela putorius, Sorex araneus, Mus musculus, Rattus norvegicus, Ochotona princeps, Cavia porcellus and Monodelphis domesticus. This generated a list of 5109 1:1 ortholog groups from these species. Amino acid and coding sequences for the species listed above were downloaded from Ensembl version 84. For the species Myotis brandtii, Eptesicus fuscus, Miniopterus natalensis, Pteropus alecto, Condylura cristata, Mesocricetus auratus, Microtus ochrogaster and Heterocephalus glaber, orthologs were identified with a recriprocal best hit approach. For example, the naked mole-rat orthologs were identified by first using blastp to align the 5,109 Cavia porcellus protein sequences from Ensembl against all naked mole-rat refSeq sequences. The top naked mole-rat hits were then aligned against all Cavia porcellus protein sequences from Ensembl. All instances in which the reciprocal blast alignments identified the same sequences were included for further study, all others were discarded. This reduced the number of 1:1 ortholog groups to 4,978.

Amino acid sequences for each of the 4,978 1:1 ortholog groups were aligned with Muscle [46]. TrimAl [47] was used to filter alignments with mean percentage identity below 60%. This left 4,628 high quality multiple sequence alignments. Pal2Nal [48] was then used to create codon oriented nucleotide alignments for the 4,628 ortholog groups. We then constructed Maximum-likelihood gene trees in PhyML 3.0 [49] using the best fitting nucleotide substitution model for each gene as indicated by jModelTest 2 [50]. A custom perl script was used to analyze the resulting Newick files, searching for monophyletic groupings of the long-lived species.

PAML, specifically the program CodeML [51], was used to reconstruct ancestral amino acid sequences at all interior nodes of the species tree for ADAMTS9. We then used custom Perl scripts to map and count all convergent substitutions occurring between the naked mole-rat and microbat branches. Since detection of convergent and divergent substitutions is critically dependent upon alignment quality we used PRANK [52], a phylogeny aware alignment tool, to generate a second ADAMTS9 multiple sequence alignment. Differences between the Muscle and PRANK alignments were minimal and both recovered the same 21 convergent substitutions between microbats and naked mole-rat.

Both the Muscle and PRANK ADAMTS9 amino acid multiple sequence alignments, the accepted tree topology and the program Grand-conv were used to estimate the posterior numbers of convergent and divergent substitutions. The two alignments produced identical results. To determine which sites were subject to positive selection along the microbat lineage we used the program TreeSAAP. The ADAMTS9 nucleotide multiple sequence alignment and the accepted species tree were used as input. We only considered radical substitutions (categories 6-8) with Z-scores greater than 3.09. PhyML was used to calculate site-wise log likelihood values for both the convergent tree and a constrained species tree using a multiple sequence alignment containing only sites with evidence of positive selection as indicated by TreeSAAP. An approximately unbiased test comparing the likeliness of the convergent and species trees was conducted in Consel [53].

Supplementary Materials

Supplementary File

Author Contributions

MJL and CVP conceived of the study. MJL performed all analyses. MJL and CVP wrote the manuscript.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Funding

Funding provided by Washington State University from the Herbert L. Eastlick Distinguished Professorship to CVP.

References

  • 1. Promislow DE. On size and survival: progress and pitfalls in the allometry of life span. J Gerontol. 1993; 48:B115–23. https://doi.org/10.1093/geronj/48.4.B115 [PubMed]
  • 2. de Magalhães JP, Costa J, Church GM. An analysis of the relationship between metabolism, developmental schedules, and longevity using phylogenetic independent contrasts. J Gerontol A Biol Sci Med Sci. 2007; 62:149–60. https://doi.org/10.1093/gerona/62.2.149 [PubMed]
  • 3. Podlutsky AJ, Khritankov AM, Ovodov ND, Austad SN. A new field record for bat longevity. J Gerontol A Biol Sci Med Sci. 2005; 60:1366–68. https://doi.org/10.1093/gerona/60.11.1366 [PubMed]
  • 4. Buckles EL. Chiroptera (bats). In: Fowler ME, Miller RE (eds.). Zoo and wildlife medicine, Volume 8. St. Louis (MO): Elsevier; 2014. p. 285–286.
  • 5. Chu P, Shuo Y, Wang F, Jeng C, Pang V, Chang P, Chin S, Liu C. Spontaneous neoplasms in zoo mammals, birds, and reptiles in Taiwan—a 10 year survey. Anim Biol. 2012; 62:95–110. https://doi.org/10.1163/157075611X616941
  • 6. Brunet-Rossinni A, Wilkinson G. Methods for age estimation and the study of senescence in bats. T. Kunz, S. Parsons (Eds.), Ecological and Behavioral Methods for the Study of Bats, The Johns Hopkins University Press, Baltimore, MD, USA (2009), pp. 315–325.
  • 7. Buffenstein R, Jarvis JU. The naked mole rat--a new record for the oldest living rodent. Sci SAGE KE. 2002; 2002:pe7. https://doi.org/10.1126/sageke.2002.21.pe7 [PubMed]
  • 8. Buffenstein R. The naked mole-rat: a new long-living model for human aging research. J Gerontol A Biol Sci Med Sci. 2005; 60:1369–77. https://doi.org/10.1093/gerona/60.11.1369 [PubMed]
  • 9. Liang S, Mele J, Wu Y, Buffenstein R, Hornsby PJ. Resistance to experimental tumorigenesis in cells of a long-lived mammal, the naked mole-rat (Heterocephalus glaber). Aging Cell. 2010; 9:626–35. https://doi.org/10.1111/j.1474-9726.2010.00588.x [PubMed]
  • 10. Delaney MA, Nagy L, Kinsel MJ, Treuting PM. Spontaneous histologic lesions of the adult naked mole rat (Heterocephalus glaber): a retrospective survey of lesions in a zoo population. Vet Pathol. 2013; 50:607–21. https://doi.org/10.1177/0300985812471543 [PubMed]
  • 11. Edrey YH, Hanes M, Pinto M, Mele J, Buffenstein R. Successful aging and sustained good health in the naked mole rat: a long-lived mammalian model for biogerontology and biomedical research. ILAR J. 2011; 52:41–53. https://doi.org/10.1093/ilar.52.1.41 [PubMed]
  • 12. Liu Y, Cotton JA, Shen B, Han X, Rossiter SJ, Zhang S. Convergent sequence evolution between echolocating bats and dolphins. Curr Biol. 2010; 20:R53–54. https://doi.org/10.1016/j.cub.2009.11.058 [PubMed]
  • 13. Liu Z, Qi FY, Zhou X, Ren HQ, Shi P. Parallel sites implicate functional convergence of the hearing gene prestin among echolocating mammals. Mol Biol Evol. 2014; 31:2415–24. https://doi.org/10.1093/molbev/msu194 [PubMed]
  • 14. Shen YY, Liang L, Li GS, Murphy RW, Zhang YP. Parallel evolution of auditory genes for echolocation in bats and toothed whales. PLoS Genet. 2012; 8:e1002788. https://doi.org/10.1371/journal.pgen.1002788 [PubMed]
  • 15. Lambert M, Nevue A, Portfors C. Contrasting patterns of adaptive sequence convergence among echolocating mammals. Gene. 2017; 605:1–4. https://doi.org/10.1016/j.gene.2016.12.017 [PubMed]
  • 16. Ujvari B, Casewell NR, Sunagar K, Arbuckle K, Wüster W, Lo N, O’Meally D, Beckmann C, King GF, Deplazes E, Madsen T. Widespread convergence in toxin resistance by predictable molecular evolution. Proc Natl Acad Sci USA. 2015; 112:11911–16. https://doi.org/10.1073/pnas.1511706112 [PubMed]
  • 17. Liu Z, Wang W, Zhang TZ, Li GH, He K, Huang JF, Jiang XL, Murphy RW, Shi P. Repeated functional convergent effects of NaV1.7 on acid insensitivity in hibernating mammals. Proc Biol Sci. 2013; 281:20132950. https://doi.org/10.1098/rspb.2013.2950 [PubMed]
  • 18. Ricklefs RE. Evolutionary theories of aging: confirmation of a fundamental prediction, with implications for the genetic basis and evolution of life span. Am Nat. 1998; 152:24–44. https://doi.org/10.1086/286147 [PubMed]
  • 19. Williams GC. Pleiotropy, natural selection and the evolution of senescence. Evolution. 1957; 11:398–411. https://doi.org/10.2307/2406060
  • 20. Kirkwood TB. Evolution of ageing. Nature. 1977; 270:301–04. https://doi.org/10.1038/270301a0 [PubMed]
  • 21. Tacutu R, Craig T, Budovsky A, Wuttke D, Lehmann G, Taranukha D, Costa J, Fraifeld VE, de Magalhães JP. Human Ageing Genomic Resources: integrated databases and tools for the biology and genetics of ageing. Nucleic Acids Res. 2013; 41:D1027–33. https://doi.org/10.1093/nar/gks1155 [PubMed]
  • 22. Castoe TA, de Koning AP, Kim HM, Gu W, Noonan BP, Naylor G, Jiang ZJ, Parkinson CL, Pollock DD. Evidence for an ancient adaptive episode of convergent molecular evolution. Proc Natl Acad Sci USA. 2009; 106:8986–91. https://doi.org/10.1073/pnas.0900233106 [PubMed]
  • 23. Qian C, de Koning AP. Rapid Discovery of Convergent Molecular Evolution Across Entire Phylogenies. 2015. University of Calgary. http://lab.jasondk.io.
  • 24. Qian C, Bryans N, Kruykov I, de Koning AP. Visualization and analysis of statistical signatures of convergent molecular evolution. 2015. University of Calgary. http://lab.jasondk.io.
  • 25. Zhang J, Kumar S. Detection of convergent and parallel evolution at the amino acid sequence level. Mol Biol Evol. 1997; 14:527–36. https://doi.org/10.1093/oxfordjournals.molbev.a025789 [PubMed]
  • 26. Woolley S, Johnson J, Smith MJ, Crandall KA, McClellan DA. TreeSAAP: selection on amino acid properties using phylogenetic trees. Bioinformatics. 2003; 19:671–72. https://doi.org/10.1093/bioinformatics/btg043 [PubMed]
  • 27. Peng L, Yang Z, Tan C, Ren G, Chen J. Epigenetic inactivation of ADAMTS9 via promoter methylation in multiple myeloma. Mol Med Rep. 2013; 7:1055–61. [PubMed]
  • 28. Du W, Wang S, Zhou Q, Li X, Chu J, Chang Z, Tao Q, Ng EK, Fang J, Sung JJ, Yu J. ADAMTS9 is a functional tumor suppressor through inhibiting AKT/mTOR pathway and associated with poor survival in gastric cancer. Oncogene. 2013; 32:3319–28. https://doi.org/10.1038/onc.2012.359 [PubMed]
  • 29. Vellai T, Takacs-Vellai K, Zhang Y, Kovacs AL, Orosz L, Müller F. Genetics: influence of TOR kinase on lifespan in C. elegans. Nature. 2003; 426:620. https://doi.org/10.1038/426620a [PubMed]
  • 30. Yaykasli KO, Hatipoglu OF, Yaykasli E, Yildirim K, Kaya E, Ozsahin M, Uslu M, Gunduz E. Leptin induces ADAMTS-4, ADAMTS-5, and ADAMTS-9 genes expression by mitogen-activated protein kinases and NF-ĸB signaling pathways in human chondrocytes. Cell Biol Int. 2015; 39:104–12. https://doi.org/10.1002/cbin.10336 [PubMed]
  • 31. Keaton JM, Cooke Bailey JN, Palmer ND, Freedman BI, Langefeld CD, Ng MC, Bowden DW. A comparison of type 2 diabetes risk allele load between African Americans and European Americans. Hum Genet. 2014; 133:1487–95. https://doi.org/10.1007/s00439-014-1486-5 [PubMed]
  • 32. Kong X, Zhang X, Xing X, Zhang B, Hong J, Yang W. The Association of Type 2 Diabetes Loci Identified in Genome-Wide Association Studies with Metabolic Syndrome and Its Components in a Chinese Population with Type 2 Diabetes. PLoS One. 2015; 10:e0143607. https://doi.org/10.1371/journal.pone.0143607 [PubMed]
  • 33. Whitmore SS, Braun TA, Skeie JM, Haas CM, Sohn EH, Stone EM, Scheetz TE, Mullins RF. Altered gene expression in dry age-related macular degeneration suggests early loss of choroidal endothelial cells. Mol Vis. 2013; 19:2274–97. [PubMed]
  • 34. Helisalmi S, Immonen I, Losonczy G, Resch MD, Benedek S, Balogh I, Papp A, Berta A, Uusitupa M, Hiltunen M, Kaarniranta K. ADAMTS9 locus associates with increased risk of wet AMD. Acta Ophthalmol. 2014; 92:e410. https://doi.org/10.1111/aos.12341 [PubMed]
  • 35. Pyun JA, Kim S, Cho NH, Koh I, Lee JY, Shin C, Kwack K. Genome-wide association studies and epistasis analyses of candidate genes related to age at menarche and age at natural menopause in a Korean population. Menopause. 2014; 21:522–29. https://doi.org/10.1097/GME.0b013e3182a433f7 [PubMed]
  • 36. Yoshina S, Mitani S. Loss of C. elegans GON-1, an ADAMTS9 Homolog, Decreases Secretion Resulting in Altered Lifespan and Dauer Formation. PLoS One. 2015; 10:e0133966. https://doi.org/10.1371/journal.pone.0133966 [PubMed]
  • 37. Tatar M, Bartke A, Antebi A. The endocrine regulation of aging by insulin-like signals. Science. 2003; 299:1346–51. https://doi.org/10.1126/science.1081447 [PubMed]
  • 38. Blagosklonny MV. Answering the ultimate question “what is the proximal cause of aging?”. Aging (Albany NY). 2012; 4:861–77. https://doi.org/10.18632/aging.100525 [PubMed]
  • 39. Blagosklonny MV. Aging is not programmed: genetic pseudo-program is a shadow of developmental growth. Cell Cycle. 2013; 12:3736–42. https://doi.org/10.4161/cc.27188 [PubMed]
  • 40. Demidenko ZN, Blagosklonny MV. Growth stimulation leads to cellular senescence when the cell cycle is blocked. Cell Cycle. 2008; 7:3355–61. https://doi.org/10.4161/cc.7.21.6919 [PubMed]
  • 41. Demidenko ZN, Zubova SG, Bukreeva EI, Pospelov VA, Pospelova TV, Blagosklonny MV. Rapamycin decelerates cellular senescence. Cell Cycle. 2009; 8:1888–95. https://doi.org/10.4161/cc.8.12.8606 [PubMed]
  • 42. Demidenko ZN, Shtutman M, Blagosklonny MV. Pharmacologic inhibition of MEK and PI-3K converges on the mTOR/S6 pathway to decelerate cellular senescence. Cell Cycle. 2009; 8:1896–900. https://doi.org/10.4161/cc.8.12.8809 [PubMed]
  • 43. Wesierska-Gadek J. mTOR and its link to the picture of Dorian Gray - re-activation of mTOR promotes aging. Aging (Albany NY). 2010; 2:892–93. https://doi.org/10.18632/aging.100240 [PubMed]
  • 44. Blagosklonny MV. Big mice die young but large animals live longer. Aging (Albany NY). 2013; 5:227–33. https://doi.org/10.18632/aging.100551 [PubMed]
  • 45. Vilella AJ, Severin J, Ureta-Vidal A, Heng L, Durbin R, Birney E. EnsemblCompara GeneTrees: Complete, duplication-aware phylogenetic trees in vertebrates. Genome Res. 2009; 19:327–35. https://doi.org/10.1101/gr.073585.107 [PubMed]
  • 46. Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004; 32:1792–97. https://doi.org/10.1093/nar/gkh340 [PubMed]
  • 47. Capella-Gutiérrez S, Silla-Martínez JM, Gabaldón T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics. 2009; 25:1972–73. https://doi.org/10.1093/bioinformatics/btp348 [PubMed]
  • 48. Suyama M, Torrents D, Bork P. PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignments. Nucleic Acids Res. 2006; 34:W609-12. https://doi.org/10.1093/nar/gkl315 [PubMed]
  • 49. Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010; 59:307–21. https://doi.org/10.1093/sysbio/syq010 [PubMed]
  • 50. Darriba D, Taboada GL, Doallo R, Posada D. jModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012; 9:772. https://doi.org/10.1038/nmeth.2109 [PubMed]
  • 51. Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007; 24:1586–91. https://doi.org/10.1093/molbev/msm088 [PubMed]
  • 52. Löytynoja A. Phylogeny-aware alignment with PRANK. Methods Mol Biol. 2014; 1079:155–70. https://doi.org/10.1007/978-1-62703-646-7_10 [PubMed]
  • 53. Shimodaira H, Hasegawa M. CONSEL: for assessing the confidence of phylogenetic tree selection. Bioinformatics. 2001; 17:1246–47. https://doi.org/10.1093/bioinformatics/17.12.1246 [PubMed]