# Nme Gene Family Evolutionary History Reveals Pre-Metazoan Origins and High Conservation between Humans and the Sea Anemone, Nematostella vectensis

Thomas Desvignes, Pierre Pontarotti, Julien Bobe

► To cite this version:

Thomas Desvignes, Pierre Pontarotti, Julien Bobe. Nme Gene Family Evolutionary History Reveals Pre-Metazoan Origins and High Conservation between Humans and the Sea Anemone, Nematostella vectensis. PLoS ONE, 2010, 5:11 (e15506), pp.1-12. 10.1371/journal.pone.0015506. hal-01204995

HAL Id: hal-01204995

https://hal.science/hal-01204995v1

Submitted on 4 Oct 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# Nme Gene Family Evolutionary History Reveals Pre-Metazoan Origins and High Conservation between Humans and the Sea Anemone, *Nematostella vectensis*

Thomas Desvignes¹,², Pierre Pontarotti³, Julien Bobe¹*

1 UMR 6632/IFR48, Université de Provence Aix Marseille 1/CNRS, F-13000, Marseille, France, 2 IFREMER, LALR, F-34250, Palavas les flots, France, 3 UMR 6632/IFR48, Université de Provence Aix Marseille 1/CNRS, F-13000, Marseille, France

## Abstract

**Background:** The Nme gene family is involved in multiple physiological and pathological processes such as cellular differentiation, development, metastatic dissemination, and cilia functions. Despite the known importance of Nme genes and their use as clinical markers of tumor aggressiveness, the associated cellular mechanisms remain poorly understood. Over the last 20 years, several non-vertebrate model species have been used to investigate Nme functions. However, the evolutionary history of the family remains poorly understood outside the vertebrate lineage. The aim of the study was thus to elucidate the evolutionary history of the Nme gene family in Metazoans.

**Methodology/Principal Findings:** Using a total of 21 eukaryote species including 14 metazoans, the evolutionary history of Nme genes was reconstructed in the metazoan lineage. We demonstrated that the complexity of the Nme gene family, initially thought to be restricted to chordates, was also shared by the metazoan ancestor. We also provide evidence suggesting that the complexity of the family is mainly a eukaryotic innovation, with the exception of Nme8 that is likely to be a choanoflagellate/metazoan innovation. Highly conserved gene structure, genomic linkage, and protein domains were identified among metazoans, some features being also conserved in eukaryotes. When considering the entire Nme family, the starlet sea anemone is the studied metazoan species exhibiting the most conserved gene and protein sequence features with humans. In addition, we were able to show that most of the proteins known to interact with human NME proteins were also found in starlet sea anemone.

**Conclusion/Significance:** Together, our observations further support the association of Nme genes with key cellular functions that have been conserved throughout metazoan evolution. Future investigations of evolutionarily conserved Nme gene functions using the starlet sea anemone could shed new light on a wide variety of key developmental and cellular processes.

**Citation:** Desvignes T, Pontarotti P, Bobe J (2010) Nme Gene Family Evolutionary History Reveals Pre-Metazoan Origins and High Conservation between Humans and the Sea Anemone, *Nematostella vectensis*. PLoS ONE 5(11): e15506. doi:10.1371/journal.pone.0015506

**Editor:** Sergios-Orestis Kolokotronis, American Museum of Natural History, United States of America

**Received** July 27, 2010; **Accepted** October 5, 2010; **Published** November 11, 2010

**Copyright:** © 2010 Desvignes et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** This work was supported by an INRA - IFREMER PhD fellowship to TD. The research leading to these results has received funding from the European Community's Seventh Framework Programme (FP7/2007-2013) under grant agreement no. 222719 - LIFECYCLE. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

* E-mail: Julien.Bobe@rennes.inra.fr

## Introduction

The Nme family, initially called NDPK or Nm23, was named after the identification of a novel gene associated with low metastatic potential [1]. In humans, NME genes are involved in a wide variety of physiological or pathological cellular processes including development, metastatic potential, ciliary functions, and cell differentiation and proliferation at various tissular and subcellular localization (see [2] for recent review). Despite their critical role in key developmental and pathological processes, the molecular functions of Nme genes remain poorly documented [3]. In vertebrates, Nme genes can be separated in 2 groups – group I and group II – based on their evolutionary history and protein domains [4]. Nme genes of the group I (Nm1-4) originate from a unique gene of the chordate ancestor while Nme genes of the group II (Nm 5-8) are present throughout chordate evolution [4].

Nme-related genes have been sporadically reported in Archaea [5], Eubacteria [6,7], and in several eukaryotic lineages including fungi [8], plants [9], and bilaterians [10,11]. However, the evolutionary history of the family that has led to a repertoire of 5 Nme genes in the chordate ancestor remains poorly understood. Indeed, the complexity of the Nme gene repertoire outside the chordate lineage was previously uncharacterized and existing literature suggested that the complexity of the Nme gene family was much more limited in non-chordate species. In *Dictyostelium discoideum*, two Nme-related proteins, named NdkC-2 and NdkM, and expressed in the cytosol and in the mitochondria, respectively, were used for biochemical and structural studies [12–14]. In *Drosophila melanogaster*, only one Nme-related gene, named *awd*, had been reported and intensively studied for its role in aberrant development [10]. In *Caenorhabditis elegans*, one Nme-related gene had been shown to be associated with severe developmental

<figure><img src="image_2.png" /></figure>
defects [15]. As recently stressed, the use of model species to decipher Nme gene functions is extremely beneficial and needs to be further supported [16]. However, a better understanding of the evolutionary links between the Nme genes that are found in non-vertebrate model species and their mammalian counterparts is required to allow this comparative biology approach. In order to gain insight into putatively conserved key functions of Nme genes that would have been retained throughout evolution, the aim of the present study was thus to characterize gene family complexity and protein features among metazoans. We were able to show that the complexity of the Nme family predates the metazoan radiation. We also provided evidence supporting the association of Nme genes with key cellular functions that have been conserved throughout metazoan evolution.

## Results/Discussion

### Nme gene family complexity predates Metazoan radiation

Using available sequenced genomes of 21 eukaryote species ranging from amoebozoans to humans, we were able to reconstruct, in opisthokonts – the metazoan/choanoflagellate/ fungi phylum – the evolutionary history of the Nme genes that had previously been identified in the chordate ancestor [4]. We were able to show that the metazoan and chordate ancestors share a similar Nme gene repertoire. A similar repertoire was also found in *Monosiga brevicollis*, a choanoflagellate, but not in fungi (Figures 1–5).

All non-vertebrate metazoans species displayed 2 Nme genes of the group I with the exception of *Trichoplax adhaerens*, *Drosophila melanogaster*, *Aedes aegypti*, *Tribolium castaneum*, and *Lottia gigantea*, in which a single gene could be identified (Figure 1A–B). In all studied non-metazoan eukaryotes, two Nme genes of the group I could also be identified with the exception of *Saccharomyces cerevisiae* and *Tetrahymena thermophila*. In starlet sea anemone *Nematostella vectensis*, *Caenorhabditis elegans* and *Branchiostoma floridae*, the two genes are located tandemly and thus originate from a *cis*-duplication of an ancestral gene (Table 1). In addition, a surprisingly well conserved gene synteny of group I Nme sequences was found between sea anemone and humans (Figure 1C). The location of the *Sox8/9* ancestor gene in the vicinity of Nme genes in the sea anemone and the location of *SOX8* and *SOX9* on human chromosomes 17 and 16, respectively, are consistent with the first round of whole genome duplication (1R) that gave rise to *Nme2* and *Nme3/4* in vertebrates [4]. In *Ciona intestinalis*, we have obtained evidence suggesting a duplication of a large portion of genomic DNA resulting in duplicated genes on two different chromosomes (Table 1, Figure S1). Two duplicated Nme genes of the group I originating from a common group I ancestral gene were found on chromosomes 2q and 8q (Figure S1). In all studied non-metazoan eukaryotes, *Capitela teleta*, and *Strongylocentrotus purpuratus*, two Nme genes of the group I can be found on different scaffolds. Because of the limited size of these scaffolds and the level of assembly of corresponding draft genomes, it is not currently possible to speculate on the nature of the duplication event that has led to 2 genes. In *Dictyostelium discoideum*, 2 Nme genes of the group I could be identified on different chromosomes. The topology of the phylogenetic tree displaying – for *Chlamydomonas reinhardtii*, *Nematostella vectensis*, *Capitela teleta*, *Strongylocentrotus purpuratus*, *Branchiostoma floridae*, and *Ciona intestinalis* – the group I Nme sequences more closely related within a species than between species strongly suggests that corresponding duplication events are independent and lineage-specific. It is noteworthy that for the above species, duplicated sequences remained closely related, whereas in other species, *Naegleria gruberi*, *Ustilago maydis*, *M.*

*brevicollis*, and *C. elegans* one sequence is highly divergent in comparison to the other one. Altogether, this strongly suggests independent and lineage specific gene duplications. It was previously shown that a single Nme gene of the group I was present in the vertebrate ancestor [4]. This ancestor gene subsequently duplicated differently in the different vertebrates lineages and resulted in 2 to 5 genes, depending on the species [4]. Together, our data indicate that a single Nme gene of the group I was present in the opisthokont ancestor. Our observations also suggest that a single Nme gene of the group I was present in the eukaryote ancestor. Understanding the selective factors that have led to multiple independent duplications events in the *Nme* family and the role of these proteins in non-vertebrate species would provide major insights into the evolution and functions of the family.

All studied metazoan species displayed a full set of group II Nme genes with the exception of the 4 ecdysozoan species and *T. adhaerens* (Figures 2–5). The 2 other group II Nme genes, *Nme9* and *Nme10*, that have been shown to be eutherian and vertebrate innovations, respectively [4], will not be discussed here. No *Nme7* homolog could be identified in the *C. elegans* genome thus indicating a possible gene loss after the nematode radiation. In *M. brevicollis* the *Nme7* protein is structurally highly divergent as shown but the topology of the phylogenetic tree (Figure 4A) and displays a unique and incomplete domain (Figure 4B). In insects, *Nme7* proteins also displayed specific domain structure (Figure 4B) resulting in a divergent position of the corresponding group in the phylogenetic tree (Figure 4A). As for all phylogenetic analyses reported here, the tree topology remained unchanged whether we used the full length protein sequence or only domains for the phylogenetic reconstruction. Similarly, very divergent *Nme8*-related sequences were identified in insects and *T. adhaerens* (Figure 5). In these 4 species, the *Nme8*-related proteins have lost the 3 *NDPK_TX* domains that are found in all other metazoan species, including the starlet sea anemone (Figure 5). In *M. brevicollis*, the *Nme8* protein does not display typical *Nme8* *NDPK_TX* domains but 2 different domains of the *NDPK* superfamily. It is also noteworthy that the exon/intron structure corresponding to the Thioredoxin *TRX_NDPK* domain is very well conserved among all studied choanoflagellate and metazoan species. In contrast, the exon/intron structure corresponding to the *NDPK* domains is highly divergent in insects, placozoans and choanoflagellates (Figure S2). Together, our data demonstrate that *Nme5*, *Nme6*, *Nme7*, and *Nme8* genes were already present in the genome of the common ancestor of choanoflagellates and metazoans. In non-choanoflagellate/metazoan species, all Nme proteins of the group II were found in low branching eukaryotic lineages such as heteroloboseans, green plants, amoebozoans, alveolates, and fungi (Figures 2–4, Table 2) with the exception of *Nme8*. In contrast, we failed to identify Nme genes of the group II outside the eukaryotic lineage. Together, our results strongly suggest that *Nme5*, *Nme6*, and *Nme7* emerged around Eukaryote radiation. Interestingly, none of the studied non-choanoflagellate/metazoan species display all 3 proteins suggesting lineage specific loss of *Nme5*, *Nme6*, and *Nme7* genes. The *Nme8* protein typically displays 1 Thioredoxin (TRX) domain followed by 2 or 3 complete *NDPK* domains. In non-choanoflagellates/metazoan species investigated, very few sequences displaying one thioredoxin domain could be identified but were never associated with an *NDPK* domain. We thus hypothesize that Thioredoxin domains already existed in the opisthokont ancestor and that *Nme8* emerged in the choanoflagellate/metazoan ancestor by domain shuffling of 2 or 3 *NDPK* domains. This would, however, require further analysis.

<figure><img src="image_3.png" /></figure>
A

<figure><img src="image_4.png" /></figure>

B

<figure><img src="image_5.png" /></figure>

C

Nematostella vectensis
Scaffold 129

<figure><img src="image_6.png" /></figure>

Homo sapiens

Chr 17

<figure><img src="image_7.png" /></figure>

First Round of Whole Genome Duplication (1R) (-744 MY)

Chr 16

<figure><img src="image_8.png" /></figure>
**Figure 1. Phylogenetic tree and exon-intron structure of Group I Nme in Eukaryotes and synteny conservation between human and starlet sea anemone Group I Nme.** A. The phylogenetic tree was constructed from a single multiple sequence alignment. Bootstrap values for neighbor joining, maximum parsimony, and maximum likelihood methods, respectively, are indicated for each node. Asterisks (*) indicate that the node was not recovered by the corresponding phylogenetic method. The consensus midpoint-rooted tree was calculated using the FIGENIX automated phylogenomic annotation pipeline [41]. For each sequence, species and corresponding name are shown. B. Corresponding exon/intron structures of genes coding for the proteins used for the phylogenetic tree estimation. Exon/intron structure was obtained through Ensembl, NCBI or JGI databases. When exon boundaries correspond to identical amino acid positions, the exons are displayed in color. Otherwise, exons are displayed in black. Non-coding exons are shown in grey. Exon size (in nucleotides) is indicated above each box. C. Synteny analysis between *Nematostella vectensis* NmeGp1 genes and human group I NME genes. *Nematostella vectensis* gene names are NCBI Entrez gene symbols. For clarity reasons "NEMVEDRAFT_" was removed from all gene symbols. Duplication time estimations are from TimeTree [46]. doi:10.1371/journal.pone.0015506.g001

Here, we have shown that the Nme gene repertoire of the metazoan ancestor was similar to that of the vertebrate ancestor (Figure 6). In agreement with prior studies reporting major gene losses and genomic rearrangement experienced by ecdysozoans [17–20] we show that Nme genes have highly diverged in this phylogenetic group, resulting in the loss of either functional NDPK domains or entire genes. Furthermore, we demonstrate that the complexity of the family predates metazoan radiation. We also provide evidence suggesting that the complexity of the family is mainly a eukaryotic innovation, with the exception of *Nme8* that is likely to be a choanoflagellate/metazoan innovation. This unexpectedly ancient complexity of the eukaryotic Nme gene family is in striking contrast with the existing hypothesis associating the emergence of Nme family complexity with the differentiation of Bilateralian [21] lineage. It should be stressed that this burst of complexity in the Nme family is much more ancient than what has been reported for many genes in which the expansion of the family is thought to have occurred around the metazoan radiation [18,19,22,23]. This would be in favor of the participation of Nme genes in ancestral functions and would also

be consistent with their known involvement in key biological processes such as cell proliferation and development.

## Highly conserved *Nme* sequence features exist among Metazoans and Eukaryotes

The topology of the phylogenetic tree (Figure 1A) obtained using group I Nme proteins does not match the topology of the tree of life. Indeed, *N. vectensis* sequences appear, along with *B. floridae* and *C. teleta* sequences, as a sister group of the vertebrate and insect sequences, in contrast to *T. adhaerens*, *S. purpuratus*, and *C. intestinalis* sequences that appear much more divergent. The topology of the tree is, however, consistent with the genomic structure of group I Nme genes reported on Figure 1B. A highly conserved exon/intron structure is observed between *T. adhaerens*, *N. vectensis*, *L. gigantea*, and vertebrates, while ecdysozoans, *C. teleta*, *C. intestinalis*, and *S. purpuratus* exhibit a totally different genomic organization that reflects the high gene divergence observed in these species. Interestingly, the *D. discoideum* mitochondrial NdkM shows gene and protein features that are similar (Figure 1B) to its cytosolic paralog NdkC-2, but displays a longer N-terminus

<figure><img src="image_9.png" /></figure>

**Figure 2. Phylogenetic tree, protein domains, and exon-intron structure of Nme5 in Eukaryotes.** A. The phylogenetic midpoint-rooted tree was constructed as described in Figure 1A. B. Corresponding protein domain structure of Nme proteins. Protein domain information was obtained using Genbank Conserved Domain Database [40]. Parentheses indicate that the domain type is the best hit given by NCBI CDD but is not a specific hit according to CDD default parameters. C. Exon/intron gene structure was obtained as described in Figure 1B. doi:10.1371/journal.pone.0015506.g002

<figure><img src="image_10.png" /></figure>
<figure><img src="image_11.png" /></figure>

**Figure 3. Phylogenetic tree, protein domains, and exon-intron structure of Nme6 in Eukaryotes.** A. The phylogenetic midpoint-rooted tree was constructed as described in Figure 1A. B. Corresponding protein domain structure was obtained as described in Figure 2A. C. Exon/intron gene structure was obtained as described in Figure 1B.
doi:10.1371/journal.pone.0015506.g003

<figure><img src="image_12.png" /></figure>

**Figure 4. Phylogenetic tree and protein domains of Nme7 in Eukaryotes.** A. The phylogenetic midpoint-rooted tree was constructed as described in Figure 1A. B. Corresponding protein domain structure was obtained as described in Figure 2A.
doi:10.1371/journal.pone.0015506.g004

<figure><img src="image_13.png" /></figure>
<figure><img src="image_14.png" /><figcaption>**Figure 5. Phylogenetic tree and protein domains of Nme8 in Metazoans.** A. The phylogenetic midpoint-rooted tree was constructed as described in Figure 1A. B. Corresponding protein domain structure was obtained as described in Figure 2A. C. Synteny analysis between *Nematostella vectensis* Nme8 genes and human *TXNDC3* (*NME8*) genes. *Nematostella vectensis* gene names correspond to the NCBI Entrez Gene identification. For clarity reasons the unchanged part "NEMVEDRAFT_" was removed.
doi:10.1371/journal.pone.0015506.g005</figcaption></figure>

sequence containing a mitochondrial assignation signal. In vertebrates, the Nme4 protein also displays a mitochondrial assignation signal in N-terminus sequence and was shown to be a gnathostome innovation [4]. In contrast, no mitochondrial assignation signal is found in any other Nme protein. This feature is thus likely to be a functional convergence between amoebozoon NdkM and vertebrate Nme4.

As indicated above, orthologs of Human *Nme5*, *Nme6*, *Nme7*, and *Nme8* are found in metazoans (Figure S3). The topology of the phylogenetic tree (Figure 2A) obtained using Nme5 proteins showed that ecdysozoan sequences are highly divergent. Ecdysozoan Nme5 proteins appear to be even more divergent than the early diverging eukaryotes *N. gruberi* and *M. pusilla*. Similarly to group I, gene exon/intron structure and protein domains of *N. vectensis*, *B. floridae*, *S. purpuratus*, *C. teleta*, and *L. gigantea* Nme 5 proteins (Figures 2B and 2C) were highly similar to their human counterpart, while ecdysozoan sequences exhibited different protein length and/or domains. In agreement with these observations was the remarkably conserved exon/intron structure observed between humans and sea anemone *Nme5* genomic sequences (Figure 2C). Similar observations on tree topology, protein domains, and genomic structure were made for Nme6 (Figure 3). For Nme7, no ortholog could be identified in *C. elegans*, while very divergent Nme7 sequences were identified in non-metazoan eukaryote species, *M. brevicollis*, and insects as shown by the topology of the phylogenetic tree (Figure 4A). In contrast to all other studied species, no NDPk7A domain was found (Figure 4B) in *M. brevicollis* and insects sequences whereas they could be identified in several non-opisthokont eukaryote species, thus suggesting a high divergence of ecdysozoan genes within the metazoan lineage. In addition, an incomplete NDPk7B domain was found in *A. aegypti* (Figure 4B). It should be stressed that, in contrast to insects, the Nme7 sequences of *T. adhaerens*, *N. vectensis*, and *B. floridae* were remarkably similar to their human counterpart

in terms of genomic exon/intron structure, protein size, and functional domains (Figures 4B and S4). No synteny analysis could be performed between starlet sea anemone and humans for *Nme5*, *Nme6*, and *Nme7* genes due to the limited number of genes present on corresponding sea anemone scaffolds (Table 2). Similarly to Nme7, no Nme8 ortholog was identified in *C. elegans*. In insects and *T. adhaerens*, the Nme8-related sequences that could be identified were extremely divergent (Figure 5A, Figure S2) and lacked the 3 NDPK_TX domains found in all eumetazoan Nme8 proteins, including the starlet sea anemone (Figure 5B). A conserved synteny was also identified between human and sea anemone (Figure 5C). Together our data suggest that, in contrast to *Nme5* and *Nme6* that have been identified in all investigated metazoan species, *Nme7* and *Nme8* have been lost in *C. elegans*. In insects, very divergent *Nme7* and *Nme8* genes remain. For *Nme8*, the high divergence associated with the loss of the 3 NDPK_TX domains in insects and *T. adhaerens* suggests a fast evolution of the gene possibly associated with a loss of ancestral metazoan Nme8 function.

We have shown that, in addition to the complexity of the Nme family, several highly conserved gene structures and protein domains are also conserved throughout metazoan evolution, some features being also conserved throughout the eukaryotic lineage. When considering the entire Nme family, the starlet sea anemone is the metazoan species exhibiting the most conserved gene and protein sequence features with humans.

## Starlet sea anemone as a model to investigate evolutionarily conserved Nme functions

Several non vertebrate model species, mainly *C. elegans* and *D. melanogaster*, have been used to investigate Nme functions. This fruitful approach has shed light on evolutionary conserved mechanisms involved in Alzheimer's disease [15], ciliary function [24], or epithelial integrity [25]. Nevertheless, the need for studies of

<figure><img src="image_15.png" /></figure>
<table><caption>**Table 1.** Group1 Nme proteins: names and symbols by species, accession numbers and corresponding chromosomal location.</caption><thead><tr><th>Species</th><th>Name</th><th>GenBank Acc. no.</th><th>Ensembl or JGI Acc. no.</th><th>Localisation</th><th>Position</th></tr></thead><tbody><tr><td>N. gruberi</td><td>NmeGp1NgA</td><td>XP_002672581</td><td>72250</td><td>Scaffold_52</td><td>70,213–70,711</td></tr><tr><td>N. gruberi</td><td>NmeGp1NgB</td><td>XP_002678151</td><td>60551</td><td>Scaffold_19</td><td>18,232–18,965</td></tr><tr><td>T. thermophila</td><td>NmeGp1</td><td>XP_001018488</td><td></td><td>Scaffold_8254597</td><td>472,318–473,033</td></tr><tr><td>D. discoideum</td><td>NdkC-2</td><td>XP_644519</td><td>DDB0238334</td><td>Chr 2</td><td>3,427,961–3,428,694</td></tr><tr><td>D. discoideum</td><td>NdkM</td><td>XP_641417</td><td>DDB0214817</td><td>Chr 3</td><td>2,740,248–2,741,470</td></tr><tr><td>C. reinhardtii</td><td>NmeGp1CrA</td><td>XP_001698246</td><td>58944</td><td>Chr 16</td><td>365,385–367,897</td></tr><tr><td>C. reinhardtii</td><td>NmeGp1CrB</td><td>XP_001702884</td><td>292075</td><td>Scaffold_26</td><td>140,833–142,142</td></tr><tr><td>U. maydis</td><td>NmeGp1UmA</td><td>XP_759114</td><td>2967</td><td>Contig_1_101</td><td>162,394–163,065</td></tr><tr><td>U. maydis</td><td>NmeGp1UmB</td><td>XP_758923</td><td>2776</td><td>Contig_1_94</td><td>34,001–34,701</td></tr><tr><td>S. cerevisiae</td><td>Ynk1</td><td>NP_012856</td><td>YKL067W</td><td>Chromosome XI</td><td>314,456–314,917</td></tr><tr><td>M. brevicollis</td><td>NmeGp1MbA</td><td>XP_001742597</td><td>19073</td><td>Scaffold_2</td><td>878,808–880,474</td></tr><tr><td>M. brevicollis</td><td>NmeGp1MbB</td><td>XP_001749207</td><td>28652</td><td>Scaffold_28</td><td>66,018–71,920</td></tr><tr><td>T. adhaerens</td><td>NmeGp1</td><td>XP_002115688</td><td>49508</td><td>Scaffold_11</td><td>2,334,568–2,336,179</td></tr><tr><td>N. vectensis</td><td>NmeGp1NvA</td><td>XP_001630018</td><td>169568</td><td>Scaffold_129</td><td>58,094–60,224</td></tr><tr><td>N. vectensis</td><td>NmeGp1NvB</td><td>XP_001630017</td><td>115087</td><td>Scaffold_129</td><td>51,034–52,881</td></tr><tr><td>C. elegans</td><td>NmeGp1CeA</td><td>NP_492761</td><td>F25H2.5.1</td><td>Chr I</td><td>10,553,902–10,553,101</td></tr><tr><td>C. elegans</td><td>NmeGp1CeB</td><td>NP_492819</td><td>F55A3.6</td><td>Chr I</td><td>10,784,780–10,783,989</td></tr><tr><td>T. castaneum</td><td>Awd</td><td>XP_967503</td><td></td><td>Linkage Group 3</td><td>32,141,388–32,142,070</td></tr><tr><td>Group1</td><td>D. melanogaster</td><td>Awd</td><td>NP_476761</td><td>FBpp0085223</td><td>Chr 3R</td><td>27,570,894–27,571,706</td></tr><tr><td>Nme</td><td>A. aegypti</td><td>Awd</td><td>XP_001662512</td><td>AAEL012359-PA</td><td>SuperContig1.684</td><td>87,333–88,186</td></tr><tr><td>L. gigantea</td><td>NmeGp1</td><td></td><td>205662</td><td>Scaffold_3</td><td>3,942,249–3,946,053</td></tr><tr><td>C. teleta</td><td>NmeGp1CtA</td><td></td><td>21725</td><td>Scaffold_758</td><td>80,222–84,888</td></tr><tr><td>C. teleta</td><td>NmeGp1CtB</td><td></td><td>116635</td><td>Scaffold_527</td><td>188,815–189,048</td></tr><tr><td>S. purpuratus</td><td>NmeGp1SpA</td><td>XP_799145</td><td></td><td>Scaffold71291</td><td>17,770–21,754</td></tr><tr><td>S. purpuratus</td><td>NmeGp1SpB</td><td>XP_785384</td><td></td><td>Scaffold42444</td><td>298,840–301,211</td></tr><tr><td>C. intestinalis</td><td>NmeGp1CiA</td><td>XP_002123476</td><td>ENSCINP00000011619</td><td>Chr 8q</td><td>5,908,347–5,908,808</td></tr><tr><td>C. intestinalis</td><td>NmeGp1CiB</td><td>XP_002121438</td><td>ENSCINP00000002194</td><td>Chr 2q</td><td>7,888,026–7,888,562</td></tr><tr><td>B. floridae</td><td>NmeGp1BfA</td><td>XP_002598285</td><td>57540</td><td>Scaffold_24</td><td>2,074,549–2,076,844</td></tr><tr><td>B. floridae</td><td>NmeGp1BfB</td><td>XP_002598283</td><td>69640</td><td>Scaffold_24</td><td>2,064,578–2,066,435</td></tr><tr><td>H. sapiens</td><td>NME1</td><td>NP_937818</td><td>ENSP00000337060</td><td>Chr 17</td><td>49,230,937–49,239,422</td></tr><tr><td>H. sapiens</td><td>NME2</td><td>NP_001018149</td><td>ENSP00000376888</td><td>Chr 17</td><td>49,243,639–49,249,108</td></tr><tr><td>X. tropicalis</td><td>Nme2</td><td>NP_001005140</td><td>ENSXETP00000024764</td><td>Scaffold_673</td><td>77,640–81,451</td></tr><tr><td>Nme2</td><td>D. rerio</td><td>Nme2a</td><td>NP_956264</td><td>ENSDARP00000064338</td><td>Scaffold Zv8_scaffold3117</td><td>528,950–537,135</td></tr><tr><td>D. rerio</td><td>Nme2b1</td><td>NP_571001</td><td>ENSDARP00000099319</td><td>Chr 19</td><td>48,501,402–48,504,320</td></tr><tr><td>D. rerio</td><td>Nme2b2</td><td>NP_571002</td><td>ENSDARP00000098065</td><td>Chr 19</td><td>48,507,488–48,510,647</td></tr><tr><td>H. sapiens</td><td>NME3</td><td>NP_002504</td><td>ENSP00000219302</td><td>Chr 16</td><td>1,820,321–1,821,710</td></tr><tr><td>Nme3</td><td>X. tropicalis</td><td>Nme3</td><td>NP_001005115</td><td>ENSXETP00000022770</td><td>Scaffold_27</td><td>933,363–937,018</td></tr><tr><td>D. rerio</td><td>Nme3</td><td>NP_571003</td><td>ENSDARP00000075112</td><td>Chr 3</td><td>11,995,596–12,045,134</td></tr><tr><td>H. sapiens</td><td>NME4</td><td>NP_005000</td><td>ENSP00000219479</td><td>Chr 16</td><td>447,209–450,759</td></tr><tr><td>Nme4</td><td>X. tropicalis</td><td>Nme4</td><td>NP_001039239</td><td>ENSXETP00000022726</td><td>Scaffold_27</td><td>1,305,654–1,312,714</td></tr><tr><td>D. rerio</td><td>Nme4</td><td>NP_957489</td><td>ENSDARP00000103207</td><td>Chr 3</td><td>10,909,410–10,922,804</td></tr></tbody></table>

Protein names were retrieved from Ensembl and NCBI or proposed according to the evolutionary history of the genes. Chromosomal/genomic location was obtained using Ensembl genome browser, JGI databases, or NCBI Entrez Gene when not available on Ensembl or JGI.
doi:10.1371/journal.pone.0015506.t001

Nme functions in non-vertebrate model species was recently stressed by the scientific community [16]. Here, we show that the starlet sea anemone exhibits a full set of metazoan Nme genes and shares remarkably conserved gene and protein sequence features with humans that were lost in flies and *C. elegans*. The starlet sea anemone is an emerging model [26] offering several biological features such

as separate sexes, inducible spawning, flagellated sperm, and external fertilization. This model species thus offers significant opportunities to investigate Nme gene functions and thus shed new light on Nme functions that remain poorly understood [16].

Nme proteins of the group I are involved in a wide variety of cellular and physiological processes including tumor metastatic

<figure><img src="image_16.png" /></figure>
<table><caption>**Table 2.** GroupII Nme proteins: names and symbols by species, accession numbers and corresponding chromosomal location.</caption><thead><tr><th>Species</th><th>Name</th><th>GenBank Acc. no.</th><th>Ensembl or JGI Acc. no.</th><th>Localisation</th><th>Position</th></tr></thead><tbody><tr><td>Nme5</td><td>N. gruberi</td><td>Nme5</td><td>XP_002682323</td><td>29950</td><td>Scaffold_3</td><td>679,761-680,378</td></tr><tr><td></td><td>M. pusilla</td><td>Nme5</td><td>XP_003056664</td><td>15249</td><td>Scaffold_3</td><td>169,341-169,943</td></tr><tr><td></td><td>M. brevicollis</td><td>Nme5</td><td>XP_001749876</td><td>38924</td><td>Scaffold_34</td><td>201,106-202,764</td></tr><tr><td></td><td>T. adhaerens</td><td>Nme5</td><td>XP_002112439</td><td>50304</td><td>Scaffold_5</td><td>3,451,699-3,452,497</td></tr><tr><td></td><td>N. vectensis</td><td>Nme5</td><td>XP_001631264</td><td>244029</td><td>Scaffold_105</td><td>447,049-449,367</td></tr><tr><td></td><td>C. elegans</td><td>Nme5</td><td>NP_501212</td><td>R05G6.5</td><td>Chr IV</td><td>7,508,567-7,509,459</td></tr><tr><td></td><td>T. castaneum</td><td>Nme5</td><td>XP_001811234</td><td></td><td>Linkage Group 2</td><td>13,741,320-13,742,593</td></tr><tr><td></td><td>D. melanogaster</td><td>Nme5</td><td>NP_651833</td><td>FBpp0085077</td><td>Chr 3R</td><td>26,707,514-26,709,225</td></tr><tr><td></td><td>A. aegypti</td><td>Nme5</td><td>XP_001662993</td><td>AAEL003030-PA</td><td>SuperContig1.75</td><td>2,400,537-2,515,468</td></tr><tr><td></td><td>L. gigantea</td><td>Nme5</td><td></td><td>226200</td><td>Scaffold_168</td><td>228,601-233,346</td></tr><tr><td></td><td>C. teleta</td><td>Nme5</td><td></td><td>177417</td><td>Scaffold_116</td><td>97,523-98,941</td></tr><tr><td></td><td>S. purpuratus</td><td>Nme5</td><td>XP_790390</td><td></td><td>Scaffold25557</td><td>320,279-323,900</td></tr><tr><td></td><td>C. intestinalis</td><td>Nme5</td><td>NP_001154961</td><td>ENSCINP00000008954</td><td>Chr 7q</td><td>1,031,987-1,036,756</td></tr><tr><td></td><td>B. floridae</td><td>Nme5</td><td>XP_002596479</td><td>61845</td><td>Scaffold_430</td><td>334,501-336,787</td></tr><tr><td></td><td>H. sapiens</td><td>NME5</td><td>NP_003542</td><td>ENSP00000265191</td><td>Chr 5</td><td>137,450,866-137,475,104</td></tr><tr><td></td><td>X. tropicalis</td><td>Nme5</td><td>NP_001072619</td><td>ENSXETP00000008322</td><td>Scaffold_65</td><td>2,613-8,494</td></tr><tr><td></td><td>D. rerio</td><td>Nme5</td><td>NP_001002516</td><td>ENSDARP00000060997</td><td>Chr 14</td><td>6,455,589-6,462,855</td></tr><tr><td>Nme6</td><td>D. discoideum</td><td>Nme6</td><td>XP_629447</td><td>DDB0191701</td><td>Chr 6</td><td>2,282,811-2,283,293</td></tr><tr><td></td><td>C. reinhardtii</td><td>Nme6</td><td>XP_001698136</td><td>139197</td><td>Chr 16</td><td>825,573-828,346</td></tr><tr><td></td><td>U. maydis</td><td>Nme6</td><td>XP_760135</td><td>3988</td><td>Contig_1_139</td><td>15,223-16,271</td></tr><tr><td></td><td>M. brevicollis</td><td>Nme6</td><td>XP_001744091</td><td>15435</td><td>Scaffold_5</td><td>413,751-414,530</td></tr><tr><td></td><td>T. adhaerens</td><td>Nme6</td><td>XP_002114036</td><td>58087</td><td>Scaffold_7</td><td>2,675,108-2,676,755</td></tr><tr><td></td><td>N. vectensis</td><td>Nme6</td><td>XP_001622626</td><td>140474</td><td>Scaffold_525</td><td>60,519-64,497</td></tr><tr><td></td><td>C. elegans</td><td>Nme6</td><td>NP_001021779</td><td>Y48G8AL.15</td><td>Chr I</td><td>1,254,771-1,258,597</td></tr><tr><td></td><td>D. melanogaster</td><td>Nme6</td><td>NP_572965</td><td>FBpp0073750</td><td>Chr X</td><td>14,477,710-14,478,649</td></tr><tr><td></td><td>T. castaneum</td><td>Nme6</td><td>XP_972639</td><td></td><td>Linkage Group 8</td><td>6,083,034-6,083,586</td></tr><tr><td></td><td>A. aegypti</td><td>Nme6</td><td>XP_001648448</td><td>AAEL004107-PA</td><td>SuperContig1.107</td><td>2,016,668-2,033,334</td></tr><tr><td></td><td>L. gigantea</td><td>Nme6</td><td></td><td>112601</td><td>Scaffold_15</td><td>3,040,725-3,042,846</td></tr><tr><td></td><td>C. teleta</td><td>Nme6</td><td></td><td>91319</td><td>Scaffold_32</td><td>665,868-668,578</td></tr><tr><td></td><td>S. purpuratus</td><td>Nme6</td><td>XP_001200902</td><td></td><td>Scaffold35436</td><td>11,834-18,064</td></tr><tr><td></td><td>C. intestinalis</td><td>Nme6</td><td>XP_002129729</td><td>ENSCINP00000027945</td><td>Scaffold_1779</td><td>5,509-6,021</td></tr><tr><td></td><td>B. floridae</td><td>Nme6</td><td>XP_002589414</td><td>279926</td><td>Scaffold_243</td><td>2,320,842-2,322,329</td></tr><tr><td></td><td>H. sapiens</td><td>NME6</td><td>NP_005784</td><td>ENSP00000416658</td><td>Chr 3</td><td>48,334,754-48,342,848</td></tr><tr><td></td><td>X. tropicalis</td><td>Nme6</td><td>NP_001123709</td><td>ENSXETP00000034257</td><td>Scaffold_857</td><td>287,305-293,502</td></tr><tr><td></td><td>D. rerio</td><td>Nme6</td><td>NP_571672</td><td>ENSDARP00000094574</td><td>Chr 20</td><td>18,803,906-18,812,975</td></tr><tr><td>Nme7</td><td>N. gruberi</td><td>Nme7</td><td>XP_002677897</td><td>33146</td><td>Scaffold_20</td><td>272,036-273,256</td></tr><tr><td></td><td>C. reinhardtii</td><td>Nme7</td><td>XP_001702841</td><td>180221</td><td>Chr 12</td><td>8,939,035-8,941,564</td></tr><tr><td></td><td>T. thermophila</td><td>Nme7</td><td>XP_001015884</td><td></td><td>Scaffold_8254649</td><td>122,957-124,735</td></tr><tr><td></td><td>M. brevicollis</td><td>Nme7</td><td>XP_001749106</td><td>11267</td><td>Scaffold_27</td><td>519,462-520,004</td></tr><tr><td></td><td>T. adhaerens</td><td>Nme7</td><td>XP_002108466</td><td>51403</td><td>Scaffold_1</td><td>1,598,731-1,601,855</td></tr><tr><td></td><td>N. vectensis</td><td>Nme7</td><td>XP_001626602</td><td>125694</td><td>Scaffold_222</td><td>17,379-24,092</td></tr><tr><td></td><td>T. castaneum</td><td>Nme7</td><td>XP_974333</td><td></td><td>Linkage Group 7</td><td>18,619,008-18,620,171</td></tr><tr><td></td><td>D. melanogaster</td><td>Nme7</td><td>NP_649926</td><td>FBpp0081561</td><td>Chr 3R</td><td>5,505,663-5,507,224</td></tr><tr><td></td><td>A. aegypti</td><td>Nme7</td><td>XP_001661412</td><td>AAEL011098-PA</td><td>SuperContig1.541</td><td>206,988-220,408</td></tr><tr><td></td><td>L. gigantea</td><td>Nme7</td><td></td><td>187020</td><td>Scaffold_18</td><td>1,179,880-1,185,739</td></tr><tr><td></td><td>C. teleta</td><td>Nme7</td><td></td><td>160391</td><td>Scaffold_422</td><td>28,192-33,650</td></tr><tr><td></td><td>S. purpuratus</td><td>Nme7</td><td>XP_795051</td><td></td><td>Scaffold9766</td><td>108,673-117,014</td></tr><tr><td></td><td>C. intestinalis</td><td>Nme7</td><td>NP_001155162</td><td>ENSCINP00000025129</td><td>Chr 1p</td><td>3,423,461-3,424,231</td></tr><tr><td></td><td>B. floridae</td><td>Nme7</td><td>XP_002588622</td><td>287848</td><td>Scaffold_254</td><td>342,459-356,708</td></tr><tr><td></td><td>H. sapiens</td><td>NME7</td><td>NP_037462</td><td>ENSP00000356785</td><td>Chr 1</td><td>169,101,769-169,337,205</td></tr></tbody></table>

<figure><img src="image_17.png" /></figure>
<table><caption>**Table 2.** Cont.</caption><thead><tr><th>Species</th><th>Name</th><th>GenBank Acc. no.</th><th>Ensembl or JGI Acc. no.</th><th>Localisation</th><th>Position</th></tr></thead><tbody><tr><td><i>X. tropicalis</i></td><td>Nme7</td><td>NP_988903</td><td>ENSXETP00000005150</td><td>Scaffold_169</td><td>1,646,165-1,680,298</td></tr><tr><td><i>D. rerio</i></td><td>Nme7</td><td>NP_571004</td><td>ENSDARP00000073091</td><td>Chr 6</td><td>31,135,867-31,196,533</td></tr><tr><td><i>Nme8</i></td><td>M. brevicollis</td><td>Nme8</td><td>XP_001746342</td><td>Scaffold_12</td><td>606,463-609,773</td></tr><tr><td></td><td><i>T. adhaerens</i></td><td>Nme8</td><td>XP_002110931</td><td>Scaffold_3</td><td>3,392,460-3,393,060</td></tr><tr><td></td><td><i>N. vectensis</i></td><td>Nme8</td><td>XP_001634297</td><td>Scaffold_61</td><td>602,301-615,863</td></tr><tr><td></td><td><i>T. castaneum</i></td><td>TRX-Nme8</td><td>XP_972627</td><td>Linkage Group 5</td><td>12,145,420-12,148,468</td></tr><tr><td></td><td><i>D. melanogaster</i></td><td>TRX-Nme8</td><td>NP_572772</td><td>Chr X</td><td>11,884,151-11,886,801</td></tr><tr><td></td><td><i>A. aegypti</i></td><td>TRX-Nme8</td><td>XP_001652618</td><td>SuperContig1.245</td><td>180,440-198,050</td></tr><tr><td></td><td><i>L. gigantea</i></td><td>Nme8</td><td></td><td>Scaffold_7</td><td>645,416-655,143</td></tr><tr><td></td><td><i>C. teleta</i></td><td>Nme8</td><td></td><td>Scaffold_940</td><td>8,079-11,527</td></tr><tr><td></td><td><i>S. purpuratus</i></td><td>Nme8</td><td>XP_001181827</td><td>Scaffold66657</td><td>21,472-41,768</td></tr><tr><td></td><td><i>C. intestinalis</i></td><td>Nme8</td><td>NP_001027618</td><td>Chr 9q</td><td>3,694,058-3,703,935</td></tr><tr><td></td><td><i>B. floridae</i></td><td>Nme8</td><td>XP_002597926</td><td>Scaffold_152</td><td>66,164-77,826</td></tr><tr><td></td><td><i>H. sapiens</i></td><td>NME8</td><td>NP_057700</td><td>ENSP00000199447</td><td>37,888,199-37,940,003</td></tr><tr><td></td><td><i>X. tropicalis</i></td><td>Nme8</td><td>NP_001121456</td><td>Scaffold_664</td><td>444,670-467,036</td></tr><tr><td></td><td><i>D. rerio</i></td><td>Nme8</td><td>NP_001082944</td><td>Chr9</td><td>51,493,472-51,569,898</td></tr></tbody></table>

Protein names were retrieved from Ensembl and NCBI or proposed according to the evolutionary history of the genes. Chromosomal/genomic location was obtained using Ensembl genome browser, JGI databases, or NCBI Entrez Gene when not available on Ensembl or JGI.
doi:10.1371/journal.pone.0015506.t002

potential. Using a reciprocal BLASTP strategy, the starlet sea anemone genome was searched for homologs of human proteins known to interact with group I NME proteins. Among the 48 human proteins known to interact with NME1, NME2, NME3, or NME4, 44 had a homolog in sea anemone (Table S1). For instance, homologs of proteins involved in cancer and cell cycle control, such as MIF [27] and Rac1 [28], were clearly identified in sea anemone (Table S1). In addition, NME1 and NME2 have been demonstrated to regulate the expression of specific genes such as *c-MYC* [29,30] and *p53* [31]. Interestingly, *c-MYC* and *p53* homologs could also be identified in the starlet sea anemone genome thus suggesting that at least some down-stream targets of Nme proteins are also present in this species (Table S1). Interestingly, the importance of p53 in sea anemone development was recently stressed and found to be similar to its known function in vertebrate development [32].

As previously documented, most Nme proteins of the group II, with the exception of Nme6, have been associated with ciliary functions. They have been shown to play critical roles in spermatogenesis [33–35], sperm motility [2], development [36], and human conditions associated with primary ciliary dyskinesia [37]. The existence of orthologs in non-metazoan species was however unsuspected with the exception of the report of Nme7 in *C. reinhardtii* and *T. thermophila* [24]. In the present study, we have demonstrated that *Nme5*, *Nme6*, *Nme7*, and *Nme8* genes were present in the metazoan ancestor, while *Nme5*, *Nme6*, and *Nme7* were most likely present in the eukaryote ancestor. To date, only 7 proteins are known to interact with NME proteins of the group II. We were however able to identify homologs for 6 of these interacting partners in starlet sea anemone (Table S1).

As indicated above, functional evidence exist demonstrating the importance of Nme gene for key biological processes, including development, cell proliferation, ciliary function, and cancer. We have shown here that the complexity of the Nme family predates metazoan radiation and that all Nme proteins display functional domains that have been conserved throughout evolution. Using

the starlet sea anemone we were able to show that most proteins known to interact with human NME were also found in the eumetazoan ancestor. Together, these observations suggest a participation of Nme genes in key cellular functions that have been conserved throughout evolution. In this context, the starlet sea anemone that exhibits a full set of highly conserved Metazoan group II Nme genes and appropriate biological features - such as separate sex, flagellated sperm, and asymmetrical expression patterns during development -offers major opportunities to investigate Nme functions.

## Conclusion

In summary, we demonstrated that the complexity of the Nme gene family initially thought to be restricted to chordates was also shared by the Metazoan ancestor. We also provide evidence suggesting that the complexity of the family is mainly a eukaryotic innovation, with the exception of *Nme8* that is likely to be a choanoflagellate/metazoan innovation. Remarkably conserved gene structures, genomic linkage, and protein domains were identified among metazoans, some features being also conserved in eukaryotes. When considering the entire Nme family, the starlet sea anemone is the studied metazoan species exhibiting the most conserved gene and protein sequence features with humans. In addition, we were able to show that most of the proteins known to interact with human NME proteins were also found in the starlet sea anemone. Together, our observations further support the association of Nme genes with key cellular functions that have been conserved throughout metazoan evolution.

## Materials and Methods

### Sequence analysis

All Nme sequences were identified using the following genome assemblies: human (*Homo sapiens*, Assembly GRCh37), xenopus (*Xenopus tropicalis*, Assembly V.4.1), zebrafish (*Danio rerio*, Assembly ZV8), tunicate (*Ciona intestinalis*, Assembly V.2.0), florida lancelet

<figure><img src="image_18.png" /></figure>
<figure><img src="image_19.png" /><figcaption>**Figure 6.** Schematic depiction of Nme genes evolution in Opisthokonts. The gene repertoire is shown at the root of the tree for the opisthokont ancestor. For all lineages, duplication events are shown and cis-duplications (CD) indicated. For the vertebrate lineage, 1R and 3R whole genome duplications events are shown and duplication events of group I Nme genes redrawn from [4]. Nme 9 and Nme10 that have been shown to be vertebrate innovations [4] are not displayed here for clarity reasons. The gene repertoire is given for all studied species. For clarity reasons, the complete name of the gene is not given and only the gene number is given and the color kept consistent with the corresponding gene in the opisthokont ancestor.</figcaption></figure>

doi:10.1371/journal.pone.0015506.g006

(Branchiostoma floridae, Assembly V.2.0), purple sea urchin (Strongylocentrotus purpuratus, Assembly NCBI V.2.1), fruit fly (Drosophila melanogaster, Assembly BDGP5), yellow fever mosquito (Aedes aegypti, Assembly AaegL1), red flour beetle (Tribolium castaneum, Assembly Tcas 3.0), nematode (Caenorhabditis elegans, Assembly WS214), polychaete worm (Capitella teleta, Assembly V1.0), owl limpet (Lottia gigantea, Assembly V.1.0), starlet sea anemone (Nematostella vectensis, Assembly V.1.0), placozoa (Trichoplax adhaerens, Assembly GrellBS-1999 V.1.0), marine choanoflagellate (Monosiga brevicollis, Assembly V.1.0), fungi (Saccharomyces cerevisiae, Assembly EF 2; and Ustilago maydis, Assembly 1), amoebozoa (Dictyostelium discoideum, Assembly V.2.1), alveolate (Tetrahymena thermophila, Assembly 1.1), green plants (Micromonas pusilla, Assembly V.2.0; and Chlamydomonas reinhardtii, Assembly V.4.0) and heterolobosean (Naegleria gruberi, Assembly V.1.0). A large number of sequences were obtained from the NCBI NR database using human or zebrafish protein sequences as a query. When more than one sequence was obtained, the RefSeq one was preferentially selected.

When no RefSeq sequence was available, the longest sequence was used. When no sequences were available in the NR database, BLASTP was used on the Ensembl [38] and DoE Joint Genome Institute databases. The chromosomal localization of *Nme* genes was established using the Ensembl genome browser or JGI gene information, or when not available, using the UCSC Genome Bioinformatics BLAT [39] and the NCBI Sequence Viewer. Exon/intron structure was obtained from the Ensembl, NCBI, or JGI databases. The protein domain structure of Nme proteins was obtained from the GenBank Conserved Domain Database [40].

## Phylogenetic analysis of *Nme* proteins

Phylogenetic reconstructions were performed using the automated genomic annotation platform FIGENIX [41]. For each phylogenetic tree reconstruction, all selected protein sequences were added to a single multiple sequence alignment. Sequence alignment was performed automatically by the FIGENIX pipeline using MUSCLE v3.6 [42,43]. The pipeline used is based on three

<figure><img src="image_20.png" /></figure>
different methods of phylogenetic tree reconstruction (i.e. neighbour-joining (NJ), maximum parsimony (MP), and maximum likelihood (ML)). The substitution model was calculated from data for ML while BLOSUM was used for NJ. Bootstrapping was carried out to assess node support with 1000 pseudoreplicates [44]. Support values were mapped onto a midpoint-rooted 50% majority rule consensus tree for each optimality criterion. Bootstrap values are reported for the nodes that are present in all three phylogenetic reconstruction methods. Asterisks denote the absence of a node for a given phylogenetic method.

## Synteny analysis

The synteny relationships of starlet sea anemone and human *Nme* genes were analyzed by reciprocal BLASTP on the NCBI NR database using surrounding genes of *Nme* genes in *Nematostella vectensis*. Homologous genes were considered in the analysis only when reciprocal BLASTP returned the couple as best hit. For the *Ciona intestinalis* paralogy analysis, synteny relationships were obtained using the Synteny Database [45] and putative paralogs were validated by reciprocal BLASTP on the NCBI NR database.

## Identification of *Nme* partners homologs

Validated human NME partners were obtained though NCBI Entrez Gene Interactions (http://www.ncbi.nlm.nih.gov/gene) information. Homologous genes in *Nematostella vectensis* were identified by reciprocal BLASTP on the NCBI NR database. Only BLASTP hits with an *E*-value lower than $10^{-10}$ were considered significant.

## Supporting Information

**Figure S1** *Ciona intestinalis* genomic region parology relationships between chromosomes 2q and 8q. For *Ciona intestinalis* parology analysis, synteny relationships were inferred using the Synteny Database [45] and putative paralogs were validated by reciprocal BLASTP on NCBI NR databases. (TIF)

**Figure S2** Exon/intron structure of *Nme8* genes. Exon/intron structure was obtained through Ensembl, NCBI, or JGI databases. When exon boundaries correspond to similar amino

acid positions, the exons are displayed in color. Otherwise, exons are displayed in black. Non-coding exons are shown in grey. Numbers indicate exon size in nucleotides. (TIF)

**Figure S3** Phylogenetic reconstruction of the Nme protein family in eumetazoans. Phylogenetic tree was constructed from a single multiple alignment. Bootstrap values for neighbor joining, maximum parsimony, and maximum likelihood methods, respectively, are indicated for each node. * indicates that the node does not exist in the corresponding tree. The consensus tree was calculated using the FIGENIX [41] automated phylogenomic annotation pipeline. Only *Homo sapiens*, *Danio rerio*, *Ciona intestinalis* and *Nematostella vectensis* sequences were used in this phylogenetic tree reconstruction because of the highly divergent ecdysozoans sequences greatly modifying the tree topology. (TIF)

**Figure S4** Exon/intron structure of *Nme7* genes. Exon/intron structure was obtained through Ensembl, NCBI, or JGI databases. When exon boundaries correspond to similar amino acid positions, the exons are displayed in color. Otherwise, exons are displayed in black. Non-coding exons are shown in grey. Numbers indicate exon size in nucleotides. (TIF)

**Table S1** BLASTP hits of human NME partners against starlet sea anemone sequences. Partners of human NME proteins were listed from NCBI Entrez Gene interaction information. BLASTP hits were considered significant for *e*-values lower than $10^{-10}$. Accession numbers and BLASTP *e*-values are also given for c-Myc and p53, two downstream targets of human NME1 and NME2 proteins. (XLS)

## Acknowledgments

Authors thank Alexis Fostier for helpful discussions.

## Author Contributions

Conceived and designed the experiments: TD JB PP. Performed the experiments: TD. Analyzed the data: TD JB PP. Contributed reagents/ materials/analysis tools: JB PP. Wrote the paper: TD JB. Performed in silico analyses and prepared the figures: TD.

## References

1. Steeg PS, Bevilacqua G, Kopper L, Thorgeirsson UP, Talmadge JE, et al. (1988) Evidence for a novel gene associated with low tumor metastatic potential. *J Natl Cancer Inst* 80: 200-204.

2. Boissan M, Dabernat S, Peuchant E, Schlattner U, Lascu I, et al. (2009) The mammalian Nm23/NDPK family: from metastasis control to cilia movement. *Mol Cell Biochem* 329: 51-62.

3. Postel EH (2003) Multiple biochemical activities of NM23/NDP kinase in gene regulation. *J Bioenerg Biomembr* 35: 31-40.

4. Desvignes T, Pontarotti P, Fauvel C, Bobe J (2009) Nme protein family evolutionary history, a vertebrate perspective. *BMC Evol Biol* 9: 256.

5. Polosina YY, Jarrell KF, Fedorov OV, Kostyukova AS (1998) Nucleoside diphosphate kinase from haloalkaliphilic archaeon *Natronobacterium magadii*: purification and characterization. *Extremophiles* 2: 333-338.

6. Hama H, Almaula N, Lerner CG, Inouye S, Inouye M (1991) Nucleoside diphosphate kinase from *Escherichia coli*; its overproduction and sequence comparison with eukaryotic enzymes. *Gene* 105: 31-36.

7. Lu Q, Zhang X, Almaula N, Mathews CK, Inouye M (1995) The gene for nucleoside diphosphate kinase functions as a mutator gene in *Escherichia coli*. *J Mol Biol* 254: 337-341.

8. Yang M, Jarrett SG, Craven R, Kaetzel DM (2009) YNK1, the yeast homolog of human metastasis suppressor NM23, is required for repair of UV radiation- and etoposide-induced DNA damage. *Mutat Res* 660: 74-78.

9. Hammargren J, Sundstrom J, Johansson M, Bergman P, Knorpp C (2007) On the phylogeny, expression and targeting of plant nucleoside diphosphate kinases. *Physiologia Plantarum* 129: 79-89.

10. Rosengard AM, Krutzsch HC, Shearn A, Biggs JR, Barker E, et al. (1989) Reduced Nm23/Awd protein in tumour metastasis and aberrant *Drosophila* development. *Nature* 342: 177-180.

11. Ogawa K, Takai H, Ogiwara A, Yokota E, Shimizu T, et al. (1996) Is outer arm dynein intermediate chain 1 multifunctional? *Mol Biol Cell* 7: 1895-1907.

12. Troll H, Winckler T, Lascu I, Muller N, Saurin W, et al. (1993) Separate nuclear genes encode cytosolic and mitochondrial nucleoside diphosphate kinase in *Dictyostelium discoideum*. *J Biol Chem* 268: 25469-25475.

13. Cherfils J, Morera S, Lascu I, Veron M, Janin J (1994) X-ray structure of nucleoside diphosphate kinase complexed with thymidine diphosphate and Mg2+ at 2-A resolution. *Biochemistry (Mosc)* 33: 9062-9069.

14. Lascu L, Giartosio A, Ransac S, Erent M (2000) Quaternary structure of nucleoside diphosphate kinases. *J Bioenerg Biomembr* 32: 227-236.

15. Napolitano F, D'Angelo F, Bimonte M, Perrina V, D'Ambrosio C, et al. (2008) A differential proteomic approach reveals an evolutionary conserved regulation of Nme proteins by Fe65 in *C. elegans* and mouse. *Neurochem Res* 33: 2547-2553.

16. Mehta A, Orchard S (2009) Nucleoside diphosphate kinase (NDPK, NM23, AWD): recent regulatory advances in endocytosis, metastasis, psoriasis, insulin release, fetal erythroid lineage and heart failure; translational medicine exemplified. *Mol Cell Biochem* 329: 3-15.

17. Takahashi T, McDougall C, Troscianko J, Chen WC, Jayaraman-Nagarajan A, et al. (2009) An EST screen from the annelid *Pomatoeceras lamarckii* reveals patterns of gene loss and gain in animals. *BMC Evol Biol* 9: 240.

18. Nikolaidis N, Chalkia D, Watkins DN, Barrow RK, Snyder SH, et al. (2007) Ancient origin of the new developmental superfamily DANGER. *PLoS ONE* 2: e204.

19. Ryan JF, Mazza ME, Pang K, Matus DQ, Baxevanis AD, et al. (2007) Pre-bilaterian origins of the Hox cluster and the Hox code: evidence from the sea anemone, *Nematostella vectensis*. *PLoS ONE* 2: e153.

20. Yamada A, Pang K, Martindale MQ, Tochinai S (2007) Surprisingly complex T-box gene complement in diploblastic metazoans. *Evol Dev* 9: 220-230.

21. Ishikawa N, Shimada N, Takagi Y, Ishijima Y, Fukuda M, et al. (2003) Molecular evolution of nucleoside diphosphate kinase genes: conserved core structures and multiple-layered regulatory regions. *J Bioenerg Biomembr* 35: 7-18.

<figure><img src="image_21.png" /></figure>
22. Degnan BM, Vervoort M, Larroux C, Richards GS (2009) Early evolution of metazoan transcription factors. Curr Opin Genet Dev 19: 591-599.

23. Srivastava M, Simakov O, Chapman J, Fahey B, Gauthier ME, et al. (2010) The *Amphimedon queenslandica* genome and the evolution of animal complexity. Nature 466: 720-726.

24. Ikeda T (2010) NDP kinase 7 is a conserved microtubule-binding protein preferentially expressed in ciliated cells. Cell Struct Funct 35: 23-30.

25. Woolworth JA, Nallamothu G, Hsu T (2009) The *Drosophila* metastasis suppressor gene *Nm23* homolog, *awd*, regulates epithelial integrity during oogenesis. Mol Cell Biol 29: 4679-4690.

26. Putnam NH, Srivastava M, Hellsten U, Dirks B, Chapman J, et al. (2007) Sea anemone genome reveals ancestral cumetazoan gene repertoire and genomic organization. Science 317: 86-94.

27. Jung H, Seong HA, Ha H (2008) Direct interaction between NM23-H1 and macrophage migration inhibitory factor (MIF) is critical for alleviation of MIF-mediated suppression of p53 activity. J Biol Chem 283: 32669-32679.

28. Otsuki Y, Tanaka M, Yoshii S, Kawazoe N, Nakaya K, et al. (2001) Tumor metastasis suppressor nm23H1 regulates Rac1 GTPase by interaction with Tiam1. Proc Natl Acad Sci U S A 98: 4385-4390.

29. Postel EH, Berberich SJ, Flint SJ, Ferrone CA (1993) Human *c-myc* transcription factor PuF identified as nm23-H2 nucleoside diphosphate kinase, a candidate suppressor of tumor metastasis. Science 261: 478-480.

30. Postel EH, Berberich SJ, Rooney JW, Kaetzel DM (2000) Human NM23/ nucleoside diphosphate kinase regulates gene expression through DNA binding to nuclease-hypersensitive transcriptional elements. J Bioener Biomembr 32: 277-284.

31. Jung H, Seong HA, Ha H (2007) NM23-H1 tumor suppressor and its interacting partner STRAP activate p53 function. J Biol Chem 282: 35293-35307.

32. Pankow S, Bamberger C (2007) The p53 tumor suppressor-like protein nvp63 mediates selective germ cell death in the sea anemone *Nematostella vectensis*. PLoS ONE 2: e782.

33. Munier A, Feral C, Milon L, Pinon VP, Gyapay G, et al. (1998) A new human *nm23* homologue (*nm23-H5*) specifically expressed in testis germinal cells. FEBS Lett 434: 289-294.

34. Hwang KC, Ok DW, Hong JC, Kim MO, Kim JH (2003) Cloning, sequencing, and characterization of the murine *nm23-M5* gene during mouse spermatogenesis and spermiogenesis. Biochem Biophys Res Commun 306: 198-207.

35. Choi YJ, Cho SK, Hwang KC, Park C, Kim JH, et al. (2009) *Nm23-M5* mediates round and elongated spermatid survival by regulating GPX-5 levels. FEBS Lett 583: 1292-1298.

36. Vogel P, Read R, Hansen GM, Freay LC, Zambrowicz BP, et al. (2010) Situs inversus in Dpcd/Poll-/-, Nme7-/-, and Pkd111-/- mice. Vet Pathol 47: 120-131.

37. Duriez B, Duquesnoy P, Escudier E, Bridoux AM, Escalier D, et al. (2007) A common variant in combination with a nonsense mutation in a member of the thioredoxin family causes primary ciliary dyskinesia. Proc Natl Acad Sci U S A 104: 3336-3341.

38. Flicek P, Aken BL, Ballester B, Beal K, Bragin E, et al. (2010) Ensembl's 10th year. Nucleic Acids Res 38: D557-D562.

39. Kent WJ (2002) BLAT—the BLAST-like alignment tool. Genome Res 12: 656-664.

40. Marchler-Bauer A, Anderson JB, Chitsaz F, Derbyshire MK, DeWeese-Scott C, et al. (2009) CDD: specific functional annotation with the Conserved Domain Database. Nucleic Acids Res 37: D205-D210.

41. Gouret P, Vitiello V, Balandraud N, Gilles A, Pontarotti P, et al. (2005) FIGENIX: intelligent automation of genomic annotation: expertise integration in a new software platform. BMC Bioinformatics 6: 198.

42. Edgar RC (2004) MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics 5: 113.

43. Edgar RC (2004) MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 32: 1792-1797.

44. Felsenstein J (1985) Confidence-limits on phylogenies - An approach using the bootstrap. Evolution 39: 783-791.

45. Catchen JM, Conery JS, Postlethwait JH (2009) Automated identification of conserved synteny after whole-genome duplication. Genome Res 19: 1497-1505.

46. Hedges SB, Dudley J, Kumar S (2006) TimeTree: a public knowledge-base of divergence times among organisms. Bioinformatics 22: 2971-2972.

<figure><img src="image_22.png" /></figure>