Evolution of the p53-MDM2 pathway
© The Author(s). 2017
Received: 14 March 2017
Accepted: 26 July 2017
Published: 3 August 2017
The p53 signalling pathway, which controls cell fate, has been extensively studied due to its prominent role in tumor development. The pathway includes the tumor supressor protein p53, its vertebrate paralogs p63 and p73, and their negative regulators MDM2 and MDM4. The p53/p63/p73-MDM system is ancient and can be traced in all extant animal phyla. Despite this, correct phylogenetic trees including both vertebrate and invertebrate species of the p53/p63/p73 and MDM families have not been published.
Here, we have examined the evolution of the p53/p63/p73 protein family with particular focus on the p53/p63/p73 transactivation domain (TAD) and its co-evolution with the p53/p63/p73-binding domain (p53/p63/p73BD) of MDM2. We found that the TAD and p53/p63/p73BD share a strong evolutionary connection. If one of the domains of the protein is lost in a phylum, then it seems very likely to be followed by loss of function by the other domain as well, and due to the loss of function it is likely to eventually disappear. By focusing our phylogenetic analysis to p53/p63/p73 and MDM proteins from phyla that retain the interaction domains TAD and p53/p63/p73BD, we built phylogenetic trees of p53/p63/p73 and MDM based on both vertebrate and invertebrate species. The trees follow species evolution and contain a total number of 183 and 98 species for p53/p63/p73 and MDM, respectively. We also demonstrate that the p53/p63/p73 and MDM families result from whole genome duplications.
The signaling pathway of the TAD and p53/p63/p73BD in p53/p63/p73 and MDM, respectively, dates back to early metazoan time and has since then tightly co-evolved, or disappeared in distinct lineages.
Keywordsp53 MDM Co-evolution Phylogeny
Cancer has been observed in virtually all vertebrates, regardless of body size and lifespan, while cancer-like growths have been reported in protostome invertebrates . In mammals, such as humans and mice, there are protective systems in place. As part of this system, p53, often referred to as the “guardian of the genome”, plays the important role as an anti-cancer protein. p53 is a transcription factor responsible for regulating the fate of the cell, for example during stress and DNA damage . MDM2 is the primary negative regulator of p53, keeping p53 at appropriate levels by ubiquitination in normal functioning cells . Upon stress, p53 is activated and fulfils its role as a tumor suppressor protein, for example by inducing apoptosis. p53, or the p53 pathway, is disabled in roughly half of all human cancers . Consequently, the prominent role of p53 and MDM2 in tumor suppression makes them outstanding targets for drug design , as well as highly interesting for detailed evolutionary studies [6, 7].
p53 shares ancestry with two other transcription factors, p63 and p73, which are paralogs of p53 : p63 is responsible for skin and epithelial development, while p73 plays a role in neuronal development and differentiation . In vertebrates, MDM2 belongs to a family with two members, MDM2 and MDM4. To date, members of the p53/p63/p73 and MDM families have been reported in chordates, but also in non-chordate species, such as Mytilus trossulus (bay mussel) , Ixodes scapularis (deer tick)  and Trichoplax adhaerens, a small (<1 mm) animal that is the only known living representative of the phylum Placozoa . Thus, the ubiquitous presence of both proteins suggests that they were present in the common ancestor of all present-day animals, and we thus refer to these proteins as p53/p63/p73ancestor and MDMancestor, respectively.
Interestingly, the evolutionary history of the p53/p63/p73 family has proven difficult to fully understand, and there is no published phylogenetic tree that agrees with the generally accepted tree of life for animal evolution [13–15]. For the MDM family no comprehensive phylogenetic tree has been published. To investigate the interaction between p53/p63/p73 and MDM, we have re-examined their evolutionary history. We found a strong correlation in the conservation of the interacting domains, p53/p63/p73 TAD and MDM p53/p63/p73BD. Loss of one of the domains is associated with the lack of the other domain, with few exceptions, demonstrating their functional dependence. By utilizing conserved amino acid sequences in domains with retained function, we could infer a phylogenetic relationship of metazoan genes containing p53/p63/p73 TAD and p53/p63/p73BD, respectively. These trees include both vertebrate and invertebrate species, and are consistent with the species evolution for both p53/p63/p73 and MDM. Finally, we have examined the evolution of the p53/p63/p73 TAD domain on a molecular level with regard to protein disorder and regulatory properties. We observed similarities in the phosphorylation pattern of vertebrate p53 and mollusk and annelid p53/p63/p73, which imply that the functional properties of regulation through phosphorylation were present already in the ancestor of deuterostomes (e.g. Chordata) and protostomes (e.g. Mollusca and Arthropoda).
Emergence and loss of domains within the p53/p63/p73 family
Protostomes can be divided into four phyla, where closer ancestry is shared between Annelida (ringed worms) and Mollusca, and between Arthropoda and Nematoda (roundworms), respectively. In species from Annelida and Mollusca, all four p53/p63/p73 domains are conserved, but within the Arthropoda phylum, certain domains have been lost (Fig. 1c). Species in the Arthropoda subphyla Chelicerata (including e.g., scorpions and spiders) and Myriapoda (e.g., millipeds) have a p53/p63/p73 gene that contains all four domains while species from subphyla Hexapoda (e.g., insects), and Crustacea (e.g., crayfish and crabs) contain a truncated p53/p63/p73 gene with the DNA BD and OD. Similarly, in p53/p63/p73 from Nematoda, the TAD and SAM domains have been lost, and only the DNA BD and OD are present (Fig. 1c). All extant phyla of deuterostomes (including Chordata, Hemichordata and Echinodermata) have p53/p63/p73 genes comprising all four domains, which implies that the ancestor of deuterostome species also contained a p53/p63/p73 gene with all domains. Following two whole genome duplications early in the vertebrate lineage , the three paralogs p53, p63, and p73 emerged. p63 and p73 have retained all four domains, while the SAM domain was lost in the p53 lineage and replaced with a C-terminal disordered domain involved in protein-protein interactions .
Duplications within the p53/p63/p73 family
There are several papers that have analyzed the number of p53/p63/p73 genes and which domains these contain in different species in order to understand the p53/p63/p73 evolution [6, 20, 21]. These papers often refer to the genes with the SAM domain in invertebrates as p63/p73 or p63/p73-like and to the ones lacking the SAM domain to p53 or p53-like. To infer such a relationship is however not straightforward since domains are frequently lost during evolution and hence lack of a particular domain in a protein does not confirm close relationship with another protein lacking the same domain. The SAM domain has indeed been lost at multiple occasions during the evolution of the p53/p63/p73 family. A recent study by dos Santos et al. where they published a phylogenetic tree and included duplicates of invertebrates shows that there has been multiple duplications in the evolution of the p53/p63/p73 family . For instance, the choanoflagellate Monosiga brevicollis have one copy of p53/p63/p73 with the SAM domain and one without and these are more similar to each other than to the vertebrate p53, p63 and p73 genes according to the results in dos Santos et al. Furthermore, in several hexapod species in the arthropod lineage the p53/p63/p73 gene has been duplicated at different time points: Aedes aegypti, Anopheles gambiae and Culex quinquefasciatus p53/p63/p73 gene seem to have been duplicated in the ancestor of these species as they cluster together while the p53/p63/p73 gene in Nasonia vitripennis and Tribolium castaneum have been duplicated in two separate events. In another genus of non-vertebrate chordates, Branchiostoma floridae, one of the p53/p63/p73 genes variants has lost the SAM domain while the other has retained it, they do not cluster in the phylogenetic tree, however they are neither located in a way that implies closer relationship to any of the vertebrate p53, p63 or p73 . A more recent duplication of the p53/p63/p73 gene can be found in the chordate but non-vertebrate tunicate Ciona intestinalis. In conclusion, the p53/p63/p73 genes have been duplicated multiple times during the evolution. Furthermore, after the whole genome duplications in the vertebrate lineage leading to fishes, reptiles and mammals the three distinct p53, p63 and p73 genes were retained in the majority of species. However, gene duplications in vertebrates can also be observed, for example, there are 20 copies of p53 in the African elephant Loxodonta africana .
Emergence and loss of domains within the MDM family
Similarly to p53/p63/p73, we performed BLAST searches in metazoan genome databases for MDM, and found 166 unique MDM family genes belonging to 98 species. The MDM protein family consists of four domains, the p53/p63/p73-binding domain (p53/p63/p73BD), the Acidic domain, the zinc binding domain (Zinc BD), and the RING domain (Fig. 1b). An MDM protein comprising all four domains was previously identified in the multicellular Placozoan, Trichoplax adhaerens . The MDM gene is not present in Porifera (sponges), but it can be found within the Cnidaria phylum. However, Cnidaria MDM (i.e., from the species Nematostella vectensis, Hydra vulgaris and Acropora digitifera) lacks the p53/p63/p73BD (Fig. 1c). Since the MDM gene is present in deuterostomes and protostomes, it was consequently present in the common ancestor of extant multicellular animal species. Certain domains of MDM have been lost in the protostome lineage similarly to what we observe for p53/p63/p73 (Fig. 1c). In the Mollusca, Annelida and Arthropoda subphyla Myriapoda and Chelicerata, an MDM gene comprising all four domains was identified. However, in Nematoda, the whole gene has disappeared. In the Arthropoda subphyla Hexapoda and Crustacea, the acidic domain, zinc binding domain and the RING finger domain can be identified in a few species, but not the p53/p63/p73BD. In deuterostome species, all four domains are present in both paralogs, MDM2 and MDM4.
Loss of the TAD domain in p53/p63/p73 correlates with the loss of the p53/p63/p73BD in MDM
The interaction between p53 TAD and MDM2 p53/p63/p73BD is important in mammals, since it is involved in tumor suppression. The origin of the interaction between the domains dates back to the time of early metazoan species . Similar to p53/p63/p73, the MDM gene is not present in Porifera (sponges), but can be found within the Cnidaria phylum. However, the interaction domains in MDM and p53/p63/p73 in Cnidaria are both missing (Fig. 1c). A similar correlation between loss of p53/p63/p73BD in MDM and loss of TAD in p53/p63/p73 was observed in protostomes. For example, species belonging to the Mollusca and Annelida phyla and the Arthropoda subphyla Chelicerata and Myriapoda all contain four p53/p63/p73 domains, as well as the p53/p63/p73BD of MDM. Interestingly, the p53/p63/p73BD in MDM in the Arthropoda subphyla Chelicerata and Myriapoda species Stegodyphus mimosarum (african social velvet spider), Ixodes ricius (castor bean tick), Ixodes scapularis (deer tick), Metaseiulus occidentalis (western predatory mite) and the Strigamia maritima (centipede), is less conserved in length compared to the p53/p63/p73BD in vertebrate, annelid and mollusk species. Likewise, the p53/p63/p73 TAD from these species contains a less conserved MDM binding motif. On the other hand, in the Arthropoda subphyla Hexapoda and Crustacea, we could only find truncated versions of p53/p63/p73 and MDM where the interaction domains is not present. Likewise, all species in the Nematoda phylum lack the whole MDM protein and p53/p63/p73 TAD. By contrast, all deuterostome species contain all MDM domains, as well as the p53/p63/p73 TAD. Thus, we find a clear correlation between presence of p53/p63/p73 TAD and the p53/p63/p73BD in MDM. This suggests a strong and ancient, yet dynamic co-evolution of the interaction domains TAD and p53/p63/p73BD in the p53/p63/p73-MDM regulatory pathway. However, there are a few cases that are not clear, which are detailed below.
Species that might not conform to the co-evolution hypothesis
Phylogeny of proteins containing the interacting domains produces phylogenetic trees that follow the species evolution
Co-localization of genes on paralogons confirms that p53/p63/p73 and MDM2/MDM4 result from whole genome duplications
In local gene duplications, the two duplicated genes are located in the proximity of each other, while after whole genome duplications, the duplicated gene is found on a paralogous block resulting from recombination of chromosomes. The existence of paralogons has been confirmed by comparing the chromosomal location of duplicated human genes with the location of the evolutionary connected genes in invertebrate species as Drosophila melanogaster and Caenorhabditis elegans, which did not undergo whole genome duplications . The duplicated genes were further investigated by phylogenetic and molecular clock analysis to find the time point of the duplication, which was estimated to be around the time of early vertebrate evolution . Present day mammalian p53, p63 and p73, as well as MDM2 and MDM4, have been suggested to result from these two whole genome duplications in the vertebrate lineage, only due to their time point of divergence [13, 28]. That the duplications occur at the time point of the whole genome duplications is supported by our phylogenetic analysis, where the time of duplication happened after the divergence of Vertebrata and Agnatha (Jawless fish). For the p53/p63/p73 family, one copy was subsequently lost, and in case of MDM, two copies were lost. To confirm that the p53/p63/p73 and MDM family genes evolved from whole genome duplication events, we analyzed genes that are co-localized in paralogous chromosomal regions (synteny). p63 is located on chromosome 3, p53 is located on chromosome 17, and p73 is located on chromosome 1. These three regions form a paralogous block , hence supporting that the vertebrate p53/p63/p73 family members arose through whole genome duplications (Fig. 3a). Likewise, the location of MDM2 and MDM4 can be traced to a paralogous block in chromosome 12 and 1, respectively (Fig. 3b). These results strongly suggest that the p53/p63/p73 and MDM genes arose from the whole genome duplications in the vertebrate lineage. In teleost fish, an additional whole genome duplication occurred after the divergence from present day tetrapods , implying that two copies of p53, p63, p73, MDM2 and MDM4, respectively, can be present in some teleost fish species. However, we did not find any instances where the duplicated genes were preserved suggesting they have been lost, which is a common event.
Evolution of phosphorylation sites in the p53-TAD domain
Evolution of residual helicity in p53/p63/p73 TAD domain
The molecular evolution of intrinsically disordered proteins (IDPs) is known to have less constraints and is more prone to insertions and deletions compared to structured domains . However, binding motifs, amino acid composition, and the length of IDPs are generally conserved . Computational analysis of the primary structure in disordered regions in different species can provide some insights with regard to important residues that have persisted in the evolutionary process. Prolines are of interest when considering the residual helicity, since they sterically hinder the continuation of helical structures. The human p53 TAD has two N-terminal prolines and one C-terminal proline present at the respective end of the FxxxWxxL binding motif (Fig. 4). In a recent study  the prolines of human p53 TAD were mutated to alanine to assess the effect of the helical structure on binding affinity to MDM2 p53/p63/73BD and general function of p53. The study revealed that the N-terminal prolines (position 12 and 13) have no effect on binding, while mutation of the C-terminal proline (position 27) results in higher residual helicity and a higher affinity for the p53/p63/73BD. The C-terminal proline is conserved in the vertebrate lineage for p53. Human p63 and p73 also have a proline C-terminal of the TAD binding motif, at position 65 and 24, respectively, while invertebrate p53/p63/p73 TAD lacks a proline in this position (Fig. 4b). Published structures of p53 TAD  (Fig. 4a) and p73 TAD  in complex with MDM2 indicate a helical structure between positions 18-26 and 14-21, respectively. Agadir predictions  of the helical content of TAD from human p53, p63 and p73, as well as for invertebrate p53/p63/p73, indicate a very low helical content, suggesting that the degree of disorder in the free state is preserved in evolution irrespective of the proline (Fig. 4b). However, the conserved C-terminal proline in the vertebrate lineage of p53, p63 and p73 could provide a means for TAD to modulate helicity upon binding and thus the affinity of the interaction .
Explaining the evolution of p53/p63/p73 is challenging since no phylogenetic tree including both vertebrate and invertebrate species, which follows the species evolution, has been published . Phylogenetic trees including deuterostomes [15, 40] present a satisfying evolutionary relationship, while three trees including invertebrates [13, 14, 20] are not consistent with the species evolution. The relationships of species in these three trees are similarly inferred, where deuterostomes, mollusks, and annelids cluster together, while nematodes and arthropods are grouped together in another cluster. This is not in concordance with the evolution of the species, where mollusks, annelids, nematodes and arthropods should cluster together (Fig. 1c), indicating constraints in the gene family. For MDM, a comprehensive phylogenetic study has not been published, however, there are studies involving an evolutionary perspective of the protein family [7, 28]. Here, we present comprehensive phylogenetic trees for both the p53/p63/p73 and MDM family, in which the topology follows the species phylogeny and the whole genome duplications (Fig. 2a, b). We manage to do this by excluding genes that lack the two interacting domains, p53/p63/p73 TAD and MDM p53/p63/p73BD, or essential motifs in these domains.
Our phylogenetic trees of p53/p63/p73 and MDM exclude species belonging to the phylum Nematoda and the Arthropoda subphyla Hexapoda and Crustacea, since these genes lack the complete interaction domains. We believe that this particular limitation of genes is essential for the correct phylogenetic relationship since the species included have an evolutionary conserved p53/p63/p73 TAD: MDM p53/p63/p73BD interaction and hence have more similar constraints. There have been other attempts to create trees of the p53/p63/p73 family with only selected domains in the alignment. For instance dos Santos et al. made an alignment containing only the p53 DNA BD, which is conserved in all p53/p63/p73 family members, but the resulting tree did not follow the species evolution . The TAD domain is intrinsically disordered and has accumulated distinct mutations in different lineages, hence contains valuable evolutionary information. Intrinsically disordered domains can be difficult to align due to the high substitution rates but the conserved FxxxWxxL motif aids in aligning the less conserved regions of the TAD domain. While the TAD domain is only a small part of the whole p53/p63/p73 gene, it is likely that the combination of the TAD and the very conserved folded domains of p53/p63/p73 provides enough information for a correct phylogenetic reconstruction. In the cases where p53/p63/p73 TAD has lost its functional connection to MDM, the substitution rate increased, resulting in sequences that could easily distort a phylogenetic reconstruction.
The human p53/p63/p73BD in MDM2 and MDM4 can both interact with TAD in p53, p63, and p73, respectively . This, together with the interaction between p53/p63/p73 and MDM in bay mussel  implies that the interaction was present in the ancestor of deuterostomes and protostomes. The function of invertebrate p53/p63/p73 (and of p53/p63/p73ancestor) is thought to be protection of the germ line from DNA damage in response to stress , which is similar to the function of vertebrate p53. There is also evidence of leukemic-like disease in mollusks where p53/p63/p73 is up regulated  suggesting that p53/p63/p73 and MDM are involved in cancer in invertebrates as well as in vertebrates. Our data suggests that the TAD domain in mollusk and annelid p53/p63/p73 has a more similar phosphorylation pattern to vertebrate p53 and echinoderm p53/p63/p73 than to the vertebrate p63 and p73 family members. This leads us to propose that at the time of the split of deuterostomes and protostomes, the p53/p63/p73-MDM interaction had p53-like functionality, which has been retained in mollusk and annelid species and in p53 vertebrates. It has been suggested  that the ancestral and invertebrate function of p53/p63/p73 mainly resembles the p63 vertebrate function based on the presence of the conserved SAM domain and a greater sequence similarity between vertebrate p63 and invertebrate p53/p63/p73 . Therefore, we also propose that some functions of p53/p63/p73ancestor are more similar to that of p63 (i.e. the SAM domain functions) and others more similar to p53 (TAD domain functions). It is also possible that other functions not yet analyzed are more similar to p73, since all three family members are equally evolutionarily close to the p53/p63/p73ancestor.
Including all genes that have sequence similarity to MDM in the phylogenetic analysis does not produce a correct relationship according to the species tree. However, similarly to p53/p63/p73, when only species that contain the p53/p63/p73BD are included, the tree is in accordance with the whole genome duplications and species evolution. The MDM family shows highest conservation in the RING domain. The functional role of the RING domain in MDM2, which is conserved in all vertebrate species and jawless fish, is to form heterodimers with MDM4 stimulating MDM2 to ubiquitinate p53 . It has been reported  that MDM4 has no E3-ligase activity, which raises the question whether invertebrate MDM and MDMancestor possess E3-ligase activity.
In conclusion, the signaling pathway of the TAD and p53/p63/p73BD in p53/p63/p73 and MDM, respectively, dates back to the beginning of multicellular life and has since then tightly co-evolved. We have here, by only including genes containing the interaction domains for the first time constructed phylogenetic trees of both p53/p63/p73 and MDM, displaying a relationship in accordance with the whole genome duplications and species evolution.
Identification of p53/p63/p73 genes
p53/p63/p73 was identified in Ensembl using TBLASTN  (www.ensembl.org) and its gene tree (ENSGT00390000015092) was downloaded. In Uniprot (www.uniprot.org) the human p53 sequence was used as query to blast against all metazoan species, all hits were collected. The same search was performed in Ortho DB  (http://orthodb.org/) where all the hits were collected. Additional searches were made in NCBI (www.ncbi.nlm.nih.gov) and at the Reptilian transcriptomes webpage (http://www.reptilian-transcriptomes.org). All retrieved sequences were pooled together and duplicates were removed by using the online programme ElimDupes (www.hiv.lanl.gov/content/sequence/ELIMDUPES/elimdupes.html). The p53 TAD has a well-conserved FxxxWxxL binding motif and previous studies have shown that these are the most critical amino acids for the interaction with MDM2 [37, 45]. In the alignment we kept all sequences containing the TAD and a binding motif resembling the FxxxWxxL in amino acid character. The alignment resulted in 342 sequences from 183 species (Additional file 1: Table S1).
Alignment and phylogenetic tree of p53/p63/p73
The amino acid alignment was done in Guidance  (http://guidance.tau.ac.il) using the MAFFT algorithm with the advanced option max-iterate set to 1000 and pairwise alignment option set to localpair. Gaps where removed with a gap tolerance of 95% with Gap Strip/Squeeze v2.1.0 (http://www.hiv.lanl.gov/content/sequence/GAPSTREEZE/gap.html) and this alignment was used to generate the phylogenetic tree. Alignment of the TAD is presented in Additional file 2: Figure S1 and Sequence Logos of this alignment is presented in Additional file 3: Figure S2. The best-fit model, according to Bayesian information criterion  (BIC) was calculated using MEGA 6  and resulted in the Jones-Taylor-Thornton substitution model (JTT) model with gamma-shaped function (G) (4 categories, fixed alpha to 1.030) together with empirical amino acid equilibrium frequencies (F) and the invariant site model (I). The phylogenetic tree was generated in PhyML 3.0  (http://www.atgc-montpellier.fr/phyml/) using this model with Nearest-Neighbor-Interchange (NNI) improvement and Shimodaira-Hasegawa approximate Likelihood Ratio Test (SH-aLRT) branch support. The tree was rooted against Trichoplax adhaerens (Fig. 2a) (Additional file 4: Figure S3).
Co-localization of p53/p63/p73 genes on the same paralogous block
The human p53 gene (ENSG00000141510) is located at chromosome 17: 7,661,779-7,687,550, the p63 gene (ENSG00000073282) is found at chromosome 3: 189,631,416-189,897,279 and the p73 gene (ENSG00000078900) is located at chromosome 1: 3,652,520-3,736,201 (Fig. 3a). Searching the http://wolfe.ucd.ie/dup/human5.28/ homepage  a paralogous block in the human genome comprises chromosome 17 (5,01-8,10) and chromosome 3 (167,0-187,25), chromosome 1 (0,76-11,92) and chromosome 3(144,91-185,57), which means that all genes belonging to the p53/p63/p73 family are located on or in close proximity of the same paralogous block (Fig. 3a). VAMP2 (ENSG00000220205) is located in the proximity (5 Mb) of p53, and it has a paralog gene, VAMP3, in the proximity of p73 (10 Mb) which further confirms that these genes are a result from the whole genome duplications. The multicellular organism Trichoplax adhaerens contains a single gene of p53/63/73 (TriadG64021) located on scaffold 6 and VAMP2 give a TBLASTN hit on scaffold 6 as well.
Identification of MDM genes
MDM2 was identified using a TBLASTN search in Ensembl  (www.ensembl.org) and its gene tree (ENSGT00530000063539) was downloaded containing 142 MDM2 and MDM4 protein sequences. Additional sequences were collected using TBLASTN human MDM2 (ENST00000258149) as a query. MDM sequences lacking the p53/p63/p73BD were removed. The databases used for browsing and downloading additional sequences were Ensembl Metazoa (www.metazoa.ensembl.org), Pre Ensembl (http://pre.ensembl.org/index.html), NCBI (www.ncbi.nlm.nih.gov), Skatebase  (http://skatebase.org), Elephant Shark Genome project  (http://esharkgenome.imcb.a-star.edu.sg/), Japanese lamprey genome project (http://jlampreygenome.imcb.a-star.edu.sg/), EchinoBase  (www.echinobase.org), MOSAS amphioxus (http://genome.bucm.edu.cn/lancelet/download_data.php), Uniprot (http://www.uniprot.org/), and Botryllus schlosseri genome project  (http://botryllus.stanford.edu/botryllusgenome/). MDM proteins contain a well-conserved RING domain responsible for binding zinc, this RING domain differ from other RING domains in the binding motif. The common motif of zinc binding is Cys3HisCys4, while MDM has a unique motif, Cys2His2Cys4 . Presence of the MDM specific motif in the RING domain was a criterion for keeping the sequence in the alignment. The sequences lacking the p53/p63/p73BD were also removed from the final alignment. The alignment resulted in a total number of 166 MDM sequences from 98 species (Additional file 5: Table S2).
Alignment and phylogenetic tree of MDM
The amino acid alignment was generated in Guidance  (http://guidance.tau.ac.il) using MAFFT algorithm with the advanced option max-iterate set to 1000 and pairwise alignment option set to localpair. The alignment was lightly masked  (0.050) so that 98,9% of the amino acids remained. Gaps where removed with a gap tolerance of 95% with Gap Strip/Squeeze v2.1.0 (http://www.hiv.lanl.gov/content/sequence/GAPSTREEZE/gap.html) and this alignment was used to generate the phylogenetic tree. The alignment of the p53/p63/p73BD is presented in Additional file 6: Figure S4 and Sequence Logos of this alignment is presented in Additional file 7: Figure S5. The best-fit model, according to Bayesian information criterion  (BIC) was calculated using MEGA 6  and resulted in the Jones-Taylor-Thornton substitution (JTT) model with gamma-shaped function (4 categories, fixed alpha to 1.367) (G) together with the invariant site model (I). The phylogenetic tree was calculated using this model in PhyML 3.0  (http://www.atgc-montpellier.fr/phyml/) with Nearest-Neighbor-Interchange (NNI) improvement and Shimodaira-Hasegawa approximate Likelihood Ratio Test (SH-aLRT) branch support. The tree was rooted against Trichoplax adhaerens (Fig. 2b)(Additional file 8: Figure S6).
Co-localization of MDM genes on the same paralogous block
The human MDM2 gene (ENSG00000135679) is located at chromosome 12: 68,808,172-68,850,686 and the MDM4 gene (ENSG00000198625) is located at chromosome 1: 204,516,379-204,558,120 (Fig. 3b). Searching the http://wolfe.ucd.ie/dup/human5.28/ homepage  there is a paralogous block located on chromosome 1 (205,69-211,23) and 12 (70,14-98,25), which is in the proximity where MDM2 and MDM4 genes are located (Fig. 3b). Two other genes called PPP1R12A (ENSG00000058272) and MYF5 (ENSG00000111049) are located in the proximity (12 Mb) of MDM2 and have paralog genes, PPP1R12B (ENSG00000077157) and MYOG (ENSG000001221809) in the proximity of MDM4 (3 Mb). Thus, the genes are all located in the proximity of the same paralogous block, which is a result of whole genome duplications (Fig. 3b). The multicellular organism Trichoplax adhaerens contains an MDM ancestor (TriadG54791) located on scaffold 3:7,103,976-7,107,199. MYOG and PPP1R12A give a TBLAST hit on scaffold 3 as well, TriadG54311 and TriadG54295 respectively.
This work was supported by the Swedish Research Council.
Availability of data and materials
The datasets used and analysed during the current study are available from the corresponding author on reasonable request.
Conceptualizon EÅ, PJ and GH. Methodology EÅ, FS, WO, MG and GH. Analysis EÅ, FS, MG, PJ and GH. Writing EÅ, PJ and GH. All authors read and approved the final manuscript.
Ethics approval and consent to participate
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.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Aktipis CA, Boddy AM, Jansen G, Hibner U, Hochberg ME, Maley CC, et al. Cancer across the tree of life: cooperation and cheating in multicellularity. Philos Trans R Soc Lond B Biol Sci. 2015;370(1673).Google Scholar
- Chen J. The cell-cycle arrest and apoptotic functions of p53 in tumor initiation and progression. Cold Spring Harb Perspect Med. 2016;6:a026104.View ArticlePubMedGoogle Scholar
- Pei D, Zhang Y, Zheng J. Regulation of p53: a collaboration between Mdm2 and MdmX. Oncotarget. 2012;3:228–35.View ArticlePubMedPubMed CentralGoogle Scholar
- Lane D, Levine A. p53 research : the past thirty years and the next thirty years p53 research : the past thirty years and. Cold Spring Harb Perspect Biol. 2010;2:1–11.Google Scholar
- Khoo KH, Hoe KK, Verma CS, Lane DP. Drugging the p53 pathway: understanding the route to clinical efficacy. Nat Rev Drug Discov. 2014;13:217–36.View ArticlePubMedGoogle Scholar
- Belyi VA, Ak P, Markert E, Wang H, Hu W, Puzio-Kuter A, et al. The origins and evolution of the p53 family of genes. Cold Spring Harb Perspect Biol. 2010;2:a001198.View ArticlePubMedPubMed CentralGoogle Scholar
- Momand J, Villegas A, Belyi VA. The evolution of MDM2 family genes. Gene. 2011;486:23–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Harms KL, Chen X. The functional domains in p53 family proteins exhibit both common and distinct properties. Cell Death Differ. 2006;13:890–7.View ArticlePubMedGoogle Scholar
- Dötsch V, Bernassola F, Coutandin D, Candi E, Melino G. P63 and P73, the ancestors of P53. Cold Spring Harb Perspect Biol. 2010;2:1–15.View ArticleGoogle Scholar
- Muttray AF, O’Toole TF, Morrill W, Van Beneden RJ, Baldwin SA. An invertebrate mdm homolog interacts with p53 and is differentially expressed together with p53 and ras in neoplastic Mytilus Trossulus haemocytes. Comp Biochem Physiol B Biochem Mol Biol. 2010;156:298–308.View ArticlePubMedPubMed CentralGoogle Scholar
- Lane DP, Cheok CF, Brown CJ, Madhumalar A, Ghadessy FJ, Verma C. The Mdm2 and p53 genes are conserved in the arachnids. Cell Cycle. 2010;9:748–54.View ArticlePubMedGoogle Scholar
- Lane DP, Cheok CF, Brown C, Madhumalar A, Ghadessy FJ, Verma C. Mdm2 and p53 are highly conserved from placozoans to man. Cell Cycle. 2010;9:540–7.View ArticlePubMedGoogle Scholar
- Gomes Dos Santos H, Nunez-Castilla J, Siltberg-Liberles J. Functional Diversification after Gene Duplication: Paralog Specific Regions of Structural Disorder and Phosphorylation in p53, p63, and p73. PLoS One. 2016;11(3):1–27.Google Scholar
- Rutkowski R, Hofmann K, Gartner A. Phylogeny and function of the invertebrate p53 superfamily. Cold Spring Harb Perspect Biol. 2010;2:a001131.View ArticlePubMedPubMed CentralGoogle Scholar
- Pintus SS, Fomin ES, Oshurkov IS, Ivanisenko VA. Phylogenetic analysis of the p53 and p63/p73 gene families. In Silico Biol. 2007;7:319–32.PubMedGoogle Scholar
- King N, Westbrook MJ, Young SL, Kuo A, Abedin M, Chapman J, et al. The genome of the choanoflagellate Monosiga Brevicollis and the origin of metazoans. Nature. 2008;451:783–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Belyi VA, Levine AJ. One billion years of p53/p63/p73 evolution. Proc Natl Acad Sci U S A. 2009;106:17609–10.View ArticlePubMedPubMed CentralGoogle Scholar
- McLysaght A, Hokamp K, Wolfe KH. Extensive genomic duplication during early chordate evolution. Nat Genet. 2002;31:200–4.View ArticlePubMedGoogle Scholar
- Oldfield CJ, Meng J, Yang JY, Yang MQ, Uversky VN, Dunker a K. Flexible nets: disorder and induced fit in the associations of p53 and 14-3-3 with their partners. BMC Genomics. 2008;9 Suppl 1:S1.View ArticlePubMedGoogle Scholar
- Nedelcu AM, Tan C. Early diversification and complex evolutionary history of the p53 tumor suppressor gene family. Dev Genes Evol. 2007;217:801–6.View ArticlePubMedGoogle Scholar
- Lu W-J, Amatruda JF, Abrams JM. p53 ancestry: gazing through an evolutionary lens. Nat Rev Cancer. 2009;9:758–62.View ArticlePubMedGoogle Scholar
- Sulak M, Fong L, Mika K, Chigurupati S, Yon L, Mongan NP, et al. TP53 copy number expansion is associated with the evolution of increased body size and an enhanced DNA damage response in elephants. elife. 2016;5:e11994.PubMedPubMed CentralGoogle Scholar
- Karakostis K, Ponnuswamy A, Fusée LTS, Bailly X, Laguerre L, Worall E, Vojtesek B, Nylander K, Fåhraeus R. p53 mRNA and p53 protein structures have evolved independently to interact with MDM2; 2015. p. 1–32.Google Scholar
- Brown CJ, Johnson AK, Dunker AK, Daughdrill GW. Evolution and disorder. Curr Opin Struct Biol. 2011;21:441–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Forman-Kay JD, Mittag T. From sequence and forces to structure, function, and evolution of intrinsically disordered proteins. Structure Elsevier Ltd. 2013;21:1492–9.Google Scholar
- Letunic I, Bork P. Interactive tree of life v2: online annotation and display of phylogenetic trees made easy. Nucleic Acids Res. 2011;39:W475–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Dehal P, Boore JL. Two rounds of whole genome duplication in the ancestral vertebrate. PLoS Biol. 2005;3:e314.View ArticlePubMedPubMed CentralGoogle Scholar
- Mendoza M, Mandani G, Momand J. The MDM2 gene family. Biomol Concepts. 2014;5:9–19.View ArticlePubMedPubMed CentralGoogle Scholar
- Meyer A, Van De Peer Y. From 2R to 3R: evidence for a fish-specific genome duplication (FSGD). BioEssays. 2005;27:937–45.View ArticlePubMedGoogle Scholar
- Schon O, Friedler A, Bycroft M, Freund SM, Fersht AR. Molecular mechanism of the interaction between MDM2 and p53. J Mol Biol. 2002;323:491–501.View ArticlePubMedGoogle Scholar
- Iakoucheva LM, Radivojac P, Brown CJ, O’Connor TR, Sikes JG, Obradovic Z, et al. The importance of intrinsic disorder for protein phosphorylation. Nucleic Acids Res. 2004;32:1037–49.View ArticlePubMedPubMed CentralGoogle Scholar
- Teufel DP, Bycroft M, Fersht AR. Regulation by phosphorylation of the relative affinities of the N-terminal transactivation domains of p53 for p300 domains and Mdm2. Oncogene. 2009;28:2112–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee CW, Ferreon JC, Ferreon ACM, Arai M, Wright PE. Graded enhancement of p53 binding to CREB-binding protein (CBP) by multisite phosphorylation. Proc Natl Acad Sci U S A. 2010;107:19290–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Hornbeck PV, Zhang B, Murray B, Kornhauser JM, Latham V, Skrzypek E. PhosphoSitePlus, 2014: mutations, PTMs and recalibrations. Nucleic Acids Res. 2015;43:D512–20.View ArticlePubMedGoogle Scholar
- Light S, Sagit R, Sachenkova O, Ekman D, Elofsson A. Protein expansion is primarily due to indels in intrinsically disordered regions. Mol Biol Evol. 2013;30:2645–53.View ArticlePubMedGoogle Scholar
- Borcherds W, Theillet F, Katzer A, Finzel A, Mishall KM, Powell AT, et al. Disorder and residual helicity alter p53-Mdm2 binding affinity and signaling in cells. Nat Chem Biol. 2014;10(12):1000–2.Google Scholar
- Kussie PH, Gorina S, Marechal V, Elenbaas B, Moreau J, Levine AJ, et al. Structure of the MDM2 oncoprotein bound to the p53 tumor suppressor transactivation domain. Science. 1996;274:948–53.View ArticlePubMedGoogle Scholar
- Shin J-S, Ha J-H, Lee D-H, Ryu K-S, Bae K-H, Park BC, et al. Structural convergence of unstructured p53 family transactivation domains in MDM2 recognition. Cell Cycle. 2015;14(4):533–43.View ArticlePubMedPubMed CentralGoogle Scholar
- Muñoz V, Serrano L. Elucidating the folding problem of helical peptides using empirical parameters. Nat Struct Biol. 1994;1:399–409.View ArticlePubMedGoogle Scholar
- Coffill CR, Lee AP, Siau JW, Chee SM, Joseph TL, Tan YS, et al. The p53 – Mdm2 interaction and the E3 ligase activity of Mdm2 / Mdm4 are conserved from lampreys to humans. Genes Dev. 2016;30(3):281–92.View ArticlePubMedPubMed CentralGoogle Scholar
- Zdzalik M, Pustelny K, Kedracka-Krok S, Huben K, Pecak A, Wladyka B, et al. Interaction of regulators Mdm2 and Mdmx with transcription factors p53, p63 and p73. Cell Cycle. 2010;9:4584–91.View ArticlePubMedGoogle Scholar
- Tanimura S, Ohtsuka S, Mitsui K, Shirouzu K, Yoshimura A, Ohtsubo M. MDM2 interacts with MDMX through their RING finger domains. FEBS Lett. 1999;447:5–9.View ArticlePubMedGoogle Scholar
- Flicek P, Ahmed I, Amode MR, Barrell D, Beal K, Brent S, et al. Ensembl 2013. Nucleic Acids Res. 2013;41:D48–55.View ArticlePubMedGoogle Scholar
- Waterhouse RM, Tegenfeldt F, Li J, Zdobnov EM, Kriventseva EV. OrthoDB: a hierarchical catalog of animal, fungal and bacterial orthologs. Nucleic Acids Res. 2013;41:D358–65.View ArticlePubMedGoogle Scholar
- Li C, Pazgier M, Li C, Yuan W, Liu M, Wei G, et al. Systematic mutational analysis of peptide inhibition of the p53-MDM2/MDMX interactions. J Mol Biol. 2010;398:200–13.View ArticlePubMedPubMed CentralGoogle Scholar
- Penn O, Privman E, Ashkenazy H, Landan G, Graur D, Pupko T. GUIDANCE: a web server for assessing alignment confidence scores. Nucleic Acids Res. 2010;38:W23–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Schwarz G. Estimating the dimension of a model. Ann Stat Institute of Mathematical Statistics. 1978;6:461–4.Google Scholar
- Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Guindon S, Dufayard J-F, 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.View ArticlePubMedGoogle Scholar
- Wang Q, Arighi CN, King BL, Polson SW, Vincent J, Chen C, et al. Community annotation and bioinformatics workforce development in concert--Little Skate Genome Annotation Workshops and Jamborees. Database (Oxford). 2012;0:bar064.Google Scholar
- Venkatesh B, Lee AP, Ravi V, Maurya AK, Lian MM, Swann JB, et al. Elephant shark genome provides unique insights into gnathostome evolution. Nature. 2014;505:174–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Cameron RA, Samanta M, Yuan A, He D, Davidson E. SpBase: the sea urchin genome database and web site. Nucleic Acids Res. 2009;37:D750–4.View ArticlePubMedGoogle Scholar
- Voskoboynik A, Neff NF, Sahoo D, Newman AM, Pushkarev D, Koh W, et al. The genome sequence of the colonial chordate, Botryllus Schlosseri. elife. 2013;2:e00569.View ArticlePubMedPubMed CentralGoogle Scholar
- Privman E, Penn O, Pupko T. Improving the performance of positive selection inference by filtering unreliable alignment regions. Mol Biol Evol. 2012;29(1):1–5.View ArticlePubMedGoogle Scholar