Evolution and divergence of the mammalian SAMD9/SAMD9L gene family
© Lemos de Matos et al.; licensee BioMed Central Ltd. 2013
Received: 6 February 2013
Accepted: 6 June 2013
Published: 12 June 2013
The physiological functions of the human Sterile Alpha Motif Domain-containing 9 (SAMD9) gene and its chromosomally adjacent paralogue, SAMD9-like (SAMD9L), currently remain unknown. However, the direct links between the deleterious mutations or deletions in these two genes and several human disorders, such as inherited inflammatory calcified tumors and acute myeloid leukemia, suggest their biological importance. SAMD9 and SAMD9L have also recently been shown to play key roles in the innate immune responses to stimuli such as viral infection. We were particularly interested in understanding the mammalian evolutionary history of these two genes. The phylogeny of SAMD9 and SAMD9L genes was reconstructed using the Maximum Likelihood method. Furthermore, six different methods were applied to detect SAMD9 and SAMD9L codons under selective pressure: the site-specific model M8 implemented in the codeml program in PAML software and five methods available on the Datamonkey web server, including the Single Likelihood Ancestor Counting method, the Fixed Effect Likelihood method, the Random Effect Likelihood method, the Mixed Effects Model of Evolution method and the Fast Unbiased Bayesian AppRoximation method. Additionally, the house mouse (Mus musculus) genome has lost the SAMD9 gene, while keeping SAMD9L intact, prompting us to investigate whether this loss is a unique event during evolution.
Our evolutionary analyses suggest that SAMD9 and SAMD9L arose through an ancestral gene duplication event after the divergence of Marsupialia from Placentalia. Additionally, selection analyses demonstrated that both genes have been subjected to positive evolutionary selection. The absence of either SAMD9 or SAMD9L genes from some mammalian species supports a partial functional redundancy between the two genes.
To the best of our knowledge, this work is the first study on the evolutionary history of mammalian SAMD9 and SAMD9L genes. We conclude that evolutionary selective pressure has acted on both of these two genes since their divergence, suggesting their importance in multiple cellular processes, such as the immune responses to viral pathogens.
KeywordsSAMD9 SAMD9-like Mammals Evolutionary history Positive selection
The Sterile Alpha Motif Domain-containing 9 (SAMD9) gene is located in chromosome 7q21.2 of the human genome, and is adjacent to its close paralogue, SAMD9-like (SAMD9L), in a head-to-tail position [1, 2] and separated by approximately 12 kb. The physiological functions of both SAMD9 and SAMD9L currently remain poorly understood, but the importance of human SAMD9 has been recently emphasized during the discovery of the genetic cause of a rare life-threatening human disease, normophosphatemic familiar tumoral calcinosis (NFTC) [3, 4]. Patients with NFTC exhibited normal calcium and phosphate metabolism while developing calcified tumorous nodules at their extremities, accompanied by severe gingivitis. Two independent founder genetic events leading to the deleterious mutations in SAMD9 are responsible for the autosomal recessive disease of NFTC [3, 4]. Interestingly, these patients and their kindred are from a culturally isolated ethical group, namely Jewish-Yemenite, suggesting a potential selection pressure associated with this population [3, 4]. In addition to NFTC, misregulated human SAMD9 expression was also shown to be associated with aggressive fibromatosis, breast, and colon cancers .
Human SAMD9 expression can be upregulated by tumor necrosis factor (TNF)  or by type I  and type II interferons (IFNs) , and it is classified as an interferon-stimulated gene (ISG). Recently, an interferon regulatory factor (IRF-1) binding element was identified in the promoter region of the SAMD9 gene in humans , and overexpression of IRF-1 can lead to elevated SAMD9 gene expression . All these observations suggest a key role of SAMD9 as a signalling hub in response to innate immune stimulations. Most importantly, human SAMD9 also has very recently been shown to possess anti-viral properties in cultured cells [8, 9] emphasizing its crucial role in host defence against viral pathogens.
On the other hand, the human SAMD9L gene was shown to exhibit lower expression levels in breast cancer tissue than in normal breast tissue from the same patient . It was also identified to be an inducible gene for type I IFNs (IFNα and β), and in activated human T cells the function of SAMD9L is correlated with its IFN-induced inhibitory effects on cell migration . The murine SAMD9L gene expression was also found to be upregulated by calcitonin , suggesting a potential involvement in calcium homeostasis as well.
Lastly, the human SAMD9 and SAMD9L genes were both classified as myeloid tumor suppressors, as they are localized within a microdeletion cluster associated with myeloid disorders, such as juvenile myelomonocytic leukemia (JMML), acute myeloid leukemia (AML), and myelodysplastic syndrome (MDS) . In another study investigating altered immune responses in patients with metastatic melanoma, both SAMD9 and SAMD9L expression were shown to be significantly reduced in T and B cell populations when compared with those from healthy control individuals . It has been suggested that since these two proteins exhibit considerable sequence similarity, they may function redundantly or in related pathways, but it should be noted that patients with NFTC possess mutations only in SAMD9 and thus it is likely that the two proteins perform non-identical tasks in humans.
Evolutionarily, the orthologous genes for both SAMD9 and SAMD9L are highly conserved in many mammalian genomes, such as rat, primates and rabbit, but not in chicken, frog and fish species, or insects . This suggests that the origin of these two related genes, possibly from an ancestral duplication event, occurred at some point after branching of the mammalian species. In addition, one intriguing fact is that the house mouse genome (Mus musculus, Mumu) has lost the SAMD9 gene while maintaining SAMD9L, after an evolutionary chromosome breakage event .
The absence of SAMD9 from the house mouse (Mumu) genome led us to question if it was a unique event restricted to this taxon and stimulated the study of SAMD9 and SAMD9L evolution and divergence in different mammalian genomes. We have examined the evolutionary history and phylogeny of SAMD9 and SAMD9L, using all the available and complete mammalian genomic sequences of both genes in NCBI and Ensembl databases, in order to obtain a broader understanding of the origin of these two genes. Our deduced phylogenetic tree suggests that SAMD9 and SAMD9L indeed resulted from an ancestral gene duplication event that occurred after the divergence of Marsupialia from Placentalia. At the same time, we applied six different Maximum Likelihood (ML) methods to test for potential positive selective pressures exerted at the gene level, and we also looked for evidence of positive selection at the deduced protein level. The analyses revealed that SAMD9 and SAMD9L, at both the genome and deduced protein sequence levels, were under the effects of what appears to be sustained positive selective pressures. Our results suggest that these two proteins have been selected by long term environmental pressures, such as those exerted by pathogen responses that are under the control of innate immune regulators like the type I interferons.
SAMD9 and SAMD9Lgenes prevalence in mammals
Mammalian SAMD9 and SAMD9L genes accession numbers from species used in phylogenetic and selection analyses
Chromosome 4: 10,302,667-10,307,412a
Chromosome 9: 79,679,836-79,684,587a
Little brown myotis
Scaffold AAPE02063303: 7,766-12,520a
Chromosome 10: 35,728,133-35,732,926a
Chromosome 4: 36,749,161-36,753,927 a
Chromosome 7: 92,731,148-92,735,917 a
Chromosome 7: 92,728,829-92,747,336 a
Northern white-cheeked gibbon
SuperContig GL397261.1: 24,263,901-24,268,665 a
Chromosome 3: 124,130,532-124,147,894 a
Chromosome 7: 83,034,053-83,038,819 a
Chromosome 7: 90,353,240-90,358,009 a
Domestic Guinea pig
Scaffold_257382: 52,686-57,449 a
Canis lupus familiaris
Scaffold GL192585.1: 1,477,672-1,482,429a
Grey short-tailed opossum
West European hedgehog
GeneScaffold_8766: 48,007-52,945 a
Chromosome 10: 35,699,236-35,703,990 a
Chromosome 4: 36,788,011-36,792,765 a
Chromosome 7: 92,759,911-92,778,202 a
Chromosome 8: 54,405,622-54,420,907 a
Chromosome 7: 92,759,368-92,777,682 a
Northern white-cheeked gibbon
SuperContig GL397261.1: 24,263,209-24,320,238 a
Chromosome 3: 124,099,607-124,117,554 a
Chromosome 7: 83,003,315-83,008,287 a
Chromosome 7: 90,382,062-90,397,829 a
African bush elephant
Chromosome 4: 28,180,812-28,185,536 a
Domestic Guinea pig
scaffold_11: 24,689,192-24,742,963 a
Chromosome 6: 3,322,257-3,349,571a
scaffold 194773: 6,206-10,964 a
Special reference has to be made to two particular complete sequences that were included in our evolutionary analyses: the northern white-cheeked gibbon (Nole) SAMD9 and the domestic dog (Calu) SAMD9L. The northern white-cheeked gibbon has no SAMD9 gene currently annotated in Ensembl. However, by comparing SAMD9 sequences of other primates to the gibbon genome in Ensembl using BLAST analysis, we obtained a perfect match with a neighboring designated pseudogene of SAMD9L. Despite this biotype classification, we could not exclude this SAMD9 sequence from being considered as a bona fide gibbon SAMD9 gene. Regarding the domestic dog SAMD9L, this gene is present in NCBI and is annotated in Ensembl, but in this latter database the sequence was missing seventy-four nucleotides when compared to the sequence in NCBI. Thus, for the subsequent analyses we used only the sequence from NCBI. It should also be noted that, despite not being annotated in Ensembl, an incomplete SAMD9 sequence for the domestic dog is available in NCBI. However, when the NCBI sequence (XM_003639470.1) was analyzed by BLAST, it possessed 99 to 100% identity with a non-annotated region of chromosome 14. Since it is a non-complete nucleotide sequence, it was not used further for the study reported here.
When SAMD9 and SAMD9L were mapped in human chromosome 7, orthologous counterparts of both genes were identified in the chimpanzee (Patr), dog (Calu) and rat (Rano), but in the house mouse (Mumu) genome there was only a single genetic correspondence to the SAMD9L open reading frame in chromosome 6 . From what is currently available in Ensembl database, the absence of SAMD9 for the house mouse (Mumu) is confirmed. We checked the other available rodents to confirm the presence or absence of SAMD9 in this specific lineage. In Ensembl there is a single SAMD9 annotation for the thirteen-lined ground squirrel (Ictr). In addition, what appear to be intact SAMD9 genes have been deposited in NCBI database for the brown rat (Rano), the Chinese hamster (Crgr) and the domestic Guinea pig (Capo). On the other hand, like the house mouse (Mumu), the Ord’s kangaroo rat (Dior) does not have SAMD9 gene annotated in Ensembl database.
Complete mammalian SAMD9 and SAMD9Lgene sequences: recombination and phylogenetic analyses
The complete nucleotide coding sequences from SAMD9 and SAMD9L were aligned together (SAMD9 + SAMD9L) and translated into deduced protein sequences (Additional file 1: Figure S1). Before further phylogenetic analyses, we used the software GARD [13, 14] to look for any evidence of recombination in the alignment. Three breakpoints were identified, but only one was strongly supported by the Kishino-Hasegawa (KH) test (Additional file 2: Table S1), which should result in the estimation of a phylogenetic tree for each segment. However, since the breakpoint was located on nucleotide 4755, the genomic segment to the right of the breakpoint was only composed of 150 nucleotides.
In the estimated ML phylogenetic tree (Figure 2), SAMD9 and SAMD9L formed two well defined monophyletic groups, and within each clade we observed a concordant topology with the accepted evolutionary relationships of eutherian mammals  (Additional file 4: Figure S3). Interestingly, the marsupial grey short-tailed opossum (Modo) SAMD9L represented a highly divergent outgroup, even from the remaining SAMD9L species.
A gene duplication event after the split of marsupial and placental mammals originated SAMD9/SAMD9Lgene family
It has been previously suggested that SAMD9 and its paralogous SAMD9L may have originated from a common ancestor by a gene duplication event . In our study, the ML tree (Figure 2) topology supports this view. However, the opossum (Modo) gene annotated as SAMD9L in NCBI database (XM_001378475.1) does not cluster in the placental mammal SAMD9L group. In fact, the opossum sequence can be recognized as being in a basal position. Two highly supported eutherian monophyletic clades in the ML tree, one corresponding to all SAMD9 genes and the other one to all SAMD9L genes, were observed. The most likely evolutionary scenario can be described as following: an ancestral gene is present before the separation of marsupial from placental mammals in the common ancestor that originated the extant SAMD9L gene in the marsupial opossum (Modo) and the ancestral gene of placental SAMD9/SAMD9L gene family. Later, in placental mammals, this ancestral gene suffered an event of gene duplication resulting in the contemporary SAMD9 and SAMD9L genes.
The conservation of similar arrangement of genes in the same relative locations on the chromosomes of different species, denominated as shared synteny, can indicate the existence of a common ancestor. In Ensembl, among the mammalian species where the presence of SAMD9 and/or SAMD9L has been annotated, shared synteny can be readily observed in chromosomes and ‘gene-scaffolds’. The consistent presence of the same common flanking genes (CALCR, CCDC132, CDK6 and HEPACAM2) in different species supports the idea that SAMD9 and SAMD9L are located in highly conserved regions throughout placental mammals’ divergence and diversification (Figure 1).
Inference of positive selection at SAMD9 and SAMD9Lgenes level
SAMD9 and SAMD9L likelihood ratio test (LRT) for four site models from PAML software
M1: nearly neutral
M2: positive selection
M8: beta and ω > 1
M1: nearly neutral
M2: positive selection
M8: beta and ω > 1
Amino acid substitutions can be either conservative or radical, depending on whether they lead to a change in a certain physicochemical property . For the codons identified as being under selection, we investigated the alterations of charge and polarity between mammalian taxa. For SAMD9 all the detected codons (Figure 3) exhibited at least one physicochemical alteration across species and a maximum of five different combinations of properties were identified for codon 331. Primate species SAMD9 amino acid changes were quite conservative, since eleven codons exhibited the same amino acid. Despite the low number of species available for Artiodactyla and Rodentia, we verified in each order a great number of amino acid physicochemical alterations per codon in the SAMD9 genes. In addition, all SAMD9L codons under presumptive selection (Figure 4) exhibited physicochemical alterations across taxa and at least three properties were represented in each codon. A maximum of five different physicochemical properties were identified for codon position 452. In Primates, amino acid substitutions in SAMD9L were once again quite conservative, given that thirteen positions kept the same physicochemical properties even when amino acid substitutions happened. On the contrary, among the four Rodentia species, only three positions in SAMD9L presented the same physicochemical properties, but just one was in fact the same amino acid.
SAMD9 and SAMD9L parameter estimates and likelihood ratio test (LRT) for branch-site model A (PAML)
Branch-site Model A
Foreground branches a
Positively selected sites d
p 0 = 0.693 p 1 = 0.291 p 2a = 0.011 p 2b = 0.005 ω 0 = 0.096 ω 1 = 1.000 ω 2 = 6.192
p 0 = 0.701 p 1 = 0.293 p 2a = 0.004 p 2b = 0.002 ω 0 = 0.096 ω 1 = 1.000 ω 2 = 21.339
p 0 = 0.692 p 1 = 0.291 p 2a = 0.012 p 2b = 0.005 ω 0 = 0.095 ω 1 = 1.000 ω 2 = 2.555
p 0 = 0.696 p 1 = 0.290 p 2a = 0.010 p 2b = 0.004 ω 0 = 0.094 ω 1 = 1.000 ω 2 = 4.818
p 0 = 0.702 p 1 = 0.292 p 2a = 0.005 p 2b = 0.002 ω 0 = 0.096 ω 1 = 1.000 ω 2 = 5.728
p 0 = 0.696 p 1 = 0.286 p 2a = 0.013 p 2b = 0.005 ω 0 = 0.094 ω 1 = 1.000 ω 2 = 8.165
p 0 = 0.729 p 1 = 0.270 p 2a = 0.001 p 2b = 0.000 ω 0 = 0.139 ω 1 = 1.000 ω 2 = 409.279
p 0 = 0.714 p 1 = 0.263 p 2a = 0.016 p 2b = 0.006 ω 0 = 0.138 ω 1 = 1.000 ω 2 = 3.169
p 0 = 0.727 p 1 = 0.268 p 2a = 0.004 p 2b = 0.001 ω 0 = 0.140 ω 1 = 1.000 ω 2 = 11.372
p 0 = 0.717 p 1 = 0.262 p 2a = 0.015 p 2b = 0.006 ω 0 = 0.139 ω 1 = 1.000 ω 2 = 2.244
p 0 = 0.730 p 1 = 0.269 p 2a = 0.013 p 2b = 0.000 ω 0 = 0.140 ω 1 = 1.000 ω 2 = 20.273
p 0 = 0.728 p 1 = 0.269 p 2a = 0.003 p 2b = 0.001 ω 0 = 0.139 ω 1 = 1.000 ω 2 = 16.318
p 0 = 0.725 p 1 = 0.266 p 2a = 0.006 p 2b = 0.002 ω 0 = 0.139 ω 1 = 1.000 ω 2 = 998.998
p 0 = 0.716 p 1 = 0.264 p 2a = 0.015 p 2b = 0.005 ω 0 = 0.138 ω 1 = 1.000 ω 2 = 3.755
p 0 = 0.728 p 1 = 0.268 p 2a = 0.003 p 2b = 0.001 ω 0 = 0.139 ω 1 = 1.000 ω 2 = 38.672
p 0 = 0.727 p 1 = 0.268 p 2a = 0.004 p 2b = 0.001 ω 0 = 0.139 ω 1 = 1.000 ω 2 = 11.843
p 0 = 0.722 p 1 = 0.261 p 2a = 0.012 p 2b = 0.004 ω 0 = 0.139 ω 1 = 1.000 ω 2 = 6.984
p 0 = 0.717 p 1 = 0.264 p 2a = 0.014 p 2b = 0.005 ω 0 = 0.137 ω 1 = 1.000 ω 2 = 7.759
Inference of positive selection at SAMD9 and SAMD9L deduced proteins level
The evaluation of destabilizing radical changes that may occur in specific regions of proteins should complement the information obtained from positive selection analyses at the gene level. Using TreeSAAP software, it is possible to estimate, from a phylogenetic tree, the amino acid properties under selection from the thirty-one available in the software  (see Methods section for full list of the thirty-one properties).
From a previous study, SAMD9 and its paralogue SAMD9L have been identified in a variety of species, namely in human, chimpanzee, dog and rat. However, in the house mouse (Mus musculus, Mumu) genome, SAMD9 was uniquely lost . The same study indicated the absence of both genes in chicken, frog and all currently sequenced fish species, suggesting that the SAMD9/SAMD9L genes originating event had occurred after the mammalian radiation. One of our goals was to intensify the identification of SAMD9 and SAMD9L within different mammalian genomes and also verify whether the loss of mouse SAMD9 was a unique event restricted to this taxon.
Despite the great number of morphological, molecular and phylogenetic studies for the order Rodentia, controversies relating to the divergence times between its major suborders still persist . In a recent study on rodent evolution  some internal rodent branches have been resolved, where three main groups in the phylogenetic tree were supported: the Mouse-related clade, Ctenohystrica clade and the Squirrel-related clade. A scenario has been proposed where the pre-Squirrel-related clade diverged early from the common ancestor followed by a later separation of the pre-Mouse-related and pre-Ctenohystrica clade . We gathered sequences for one or both SAMD9 and SAMD9L genes for species representative of the three clades. The two genes were present in the thirteen-lined ground squirrel (Squirrel-related clade), the domestic Guinea pig (Ctenohystrica clade), the Chinese hamster and the brown rat (Mouse-related clade). Together with the absence of SAMD9 in the house mouse genome, the Ord’s kangaroo rat (Mouse-related clade) also did not have this gene annotated in Ensembl. With the apparent region synteny for the Ord’s kangaroo rat when compared to the other mammals, this absence might just be the case of a genome still to be completely annotated, leaving the house mouse as the only rodent taxon that has lost SAMD9, at least from the currently available genomic sequence database.
A great number of the available mammalian genomes are still not completely annotated. Therefore, we made no assumptions regarding SAMD9 and SAMD9L for those species. Nevertheless, we observed that the fairly well annotated cow and pig genomes (Order Artiodactyla) had no matches or annotations for SAMD9L. This information together with the absence of SAMD9 in the house mouse and the already suggested origin of both genes from a common ancestor by ancient gene duplication  led us to the following hypothesis: in some lineages the presence of both genes might be costly for the genome, resulting in the loss of one of the genes that functionally would be overcome by the remaining paralogue. Although these observations support the potential existence of certain gene redundancy between SAMD9 and SAMD9L, we also note the almost nonexistent recombination between them, despite the proximity in the location of these two genes in the genomes of all the annotated mammalian species. This genetic isolation of the two paralogues does not support the existence of functional redundancy between SAMD9 and SAMD9L. These apparent contradictory hypotheses have to be confirmed with the conduction of functional studies in different species.
With all the available mammalian sequences collected for both SAMD9 and SAMD9L genes, the performed phylogenetic study resulted in a tree with a well-defined monophyletic group per gene gathering solely placental mammals and a single outgroup, the marsupial grey short-tailed opossum. This supported the speculative hypothesis of SAMD9 and SAMD9L resulting from a gene duplication event, more precisely, after the divergence of Marsupialia from Placentalia 147.7 Mya . Despite the common ancestor, when testing for the occurrence of potential positive selection acting at the gene and protein levels, we concluded that SAMD9L is under stronger selection than SAMD9. This is supported by the fact that a higher number of sites at the gene level and of specific lineages were positively selected in SAMD9L than SAMD9. Besides, a greater number of amino acid properties were under selection at the deduced protein level of SAMD9L than SAMD9.
When we examined the amino acid substitutions and changes on physicochemical properties for sites under selection, it was clear, for both proteins, that members of the Rodentia order presented the highest number of divergent alterations for the same codons compared to other mammalian orders. Since it is known that in many proteins the amino acid substitutions caused by positive selection are not random [21, 26], for instances the Primate APOBEC3G residues involved in HIV-1 Vif interaction , we hypothesize that any occurring alteration in rodents or even in other lineages may be the result of consistent arms race between the host and a pathogen stressor. This could be a significant observation, given that anti-viral properties have been already assigned to human SAMD9 in cultured human cells. Specifically, a unique viral gene product, M062 of myxoma virus, was found to antagonize the anti-viral properties of SAMD9 protein in order to permit the replication of this virus in cultured human cells .
Considering the mammalian species included in this study, selection analyses performed on SAMD9 and/or SAMD9L genes for each species individually one may have different results from the obtained in our work, since recombination rates and effective population sizes are expected to differ among species. These species and population specific selection analyses should result in the identification of sites under selection in SAMD9 and/or SAMD9L genes that can be used in genetic population studies by determining parameters like allele and genotype frequencies, and FST and nucleotide diversity values. This contributes to the definition of genotypes that might be favorable or not, for example, to the defence against certain pathogens.
Human SAMD9 and SAMD9L have solely one defined domain, the Sterile Alpha Motif (SAM), a module of about 70 amino acid residues long , specifically 65 amino acids and 66 in SAMD9 and SAMD9L, respectively. SAM domains, one of the most common protein domains found in eukaryotic cells, are protein-protein interaction modules that perform a large number of different functions [29, 30] and are not easily categorized. Indeed, different SAM domains can self-associate, bind to other SAM domains and/or to non-SAM proteins, and even interact with RNA, DNA or lipids . Because of the great variety of known functions, the presence of a SAM domain does not necessarily involve a specific function or pathway, but an array of possible functions. For both human SAMD9 and SAMD9L, no function has yet been assigned to their SAM domains, but for SAMD9 the ability to form SAM polymers has been suggested . From our evolutionary study on both proteins, none of the identified sites or amino acid properties under positive selection overlapped with the deduced SAM domains, demonstrating a high level of conservation among the mammalian species.
Since the origin and evolution of the SAMD9 and SAMD9L genes were first reported, a great number of mammalian genomes have been sequenced, allowing now a more detailed view into the evolutionary history of both genes. Our study supports the previously suggested origin of SAMD9 and SAMD9L from a mammalian ancestral duplication event. Specifically, according to the results from our study, this event occurred after the divergence of Marsupialia from Placentalia. When considering the mostly complete mammalian genomes collected for this study, the apparent loss of SAMD9 or SAMD9L in some species led us to propose that some overlapping functional redundancy exists between the two proteins, despite the almost nonexistent recombination between the two closely located genes from other species. From the positive selection analyses performed, both at gene and protein levels, we demonstrate that SAMD9 and SAMD9L continue to be under long term selective pressure, with even stronger evidence for positive selection in SAMD9L.
Both SAMD9 and SAMD9L genes are upregulated by type I interferon, a classic feature associated with many innate pathogen-response genes called interferon-stimulated genes (ISGs). Indeed, human SAMD9 has already been shown to be a functional inhibitor for at least one viral pathogen, a poxvirus called myxoma virus, that expresses a specific viral inhibitor (M062) that counteracts the anti-viral properties of SAMD9 . Our results suggest that at least the SAMD9 genes may have been under sustained selection pressure exerted by viral pathogens.
Our work is the first complete study to investigate the evolutionary history of mammalian SAMD9 and SAMD9L.
SAMD9 and SAMD9L nucleotide and protein sequences
All the available mammalian SAMD9 and SAMD9L genes coding sequences used in the phylogenetic and positive selection analyses were retrieved from NCBI (http://www.ncbi.nlm.nih.gov) and Ensembl (http://www.ensembl.org/index.html) databases. Next, sequences were aligned with ClustalW  implemented in BioEdit v7.0.9 , followed by visual inspection. Nucleotide sequences translation into protein sequences was performed using also BioEdit.
SAMD9 and SAMD9L genes coding sequences were collected for fifteen and nineteen species, respectively. Based on the Mammal Species of the World database classification (http://www.bucknell.edu/msw3/), representative species of mammalian infraclasses Metatheria (Order Didelphimorphia) and Eutheria (Order Artiodactyla, Carnivora, Chiroptera, Erinaceomorpha, Lagomorpha, Perissodactyla, Primates, Proboscidea, Rodentia and Soricomorpha) were included in this study. Table 1 summarizes the species collected for each gene and their respective accession numbers.
The isoeletric point (pI) of SAMD9 and SAMD9L deduced proteins for different species was estimated using DAMBE (Data Analysis and Molecular Biology and Evolution) .
Recombination and phylogenetic analyses
Recombination can mislead phylogenetic and positive selection analyses , and particularly for SAMD9 and SAMD9L, the genes close location (~12 kb in human genome, for example) might increase the probability of recombination to occur. Therefore, we first performed recombination testing on placental SAMD9 and SAMD9L nucleotide sequences alignments, and also on the alignment of both genes together (SAMD9 + SAMD9L). The software GARD (Genetic Algorithm for Recombination Detection) [13, 14], implemented in the Datamonkey web server , was used to detect possible recombination breakpoints.
For SAMD9 and SAMD9L genes alignments no significant breakpoints were detected while using GARD, thus the complete alignments were used to establish each gene phylogeny. As indicated by the Akaike Information Criterion (AIC) implemented in jModelTest v0.1.1 , the nucleotide substitution model TVM+G was used for SAMD9 tree estimation, while the GTR+G model was the consensus model selected for SAMD9L phylogenetic tree construction. On the other hand, a significant breakpoint was detected when running GARD for the SAMD9+SAMD9L alignment and a phylogenetic tree was estimated for each segment. For the left segment, the AIC in jModelTest indicated GTR+I+G as the best-fit nucleotide substitution model, whereas for the right segment the TPM2uf+G model was indicated as the best for the tree estimation. Also, for the SAMD9+SAMD9L alignment, a phylogenetic tree was estimated without testing recombination. In this case, the jModelTest AIC estimated GTR+I+G model as the best-fit nucleotide substitution model.
To establish mammalian phylogeny for SAMD9, SAMD9L and SAMD9+SAMD9L, based on nucleotide sequences, the Maximum Likelihood (ML) method implemented on GARLI v2.0 (Genetic Algorithm for Rapid Likelihood Inference) was used . The analyses were performed with 1,000,000 generations and 1,000 bootstrap searches. ML trees were displayed using FigTree v1.3.1 (http://tree.bio.ed.ac.uk/).
Codon-based analyses of positive selection
A useful measurement for identifying adaptive protein evolution is the nonsynonymous (d N )/synonymous substitution (d S ) rate (ω = d N /d S ), where values of ω = 1, < 1, and > 1 indicate neutral selection, negative selection, and positive selection, respectively [39, 40]. Naturally, and due to protein structural and functional constraints, ω is expected to be close to 0 and full protein analysis rarely detects positive selection . As a result, several methods, based on models of codon substitution, have been developed to detect adaptive evolution (positive selection) at individual sites in a background of negative selection [42, 43]. We employed six different methods to detect sites under selection, and based on the methodology adopted by several authors [19, 20] only codons identified by at least three of the six used methods were considered to be under positive selection.
To detect selection based on the ratio ω and at the gene-level, for both SAMD9 and SAMD9L, PAML v4.4 (Phylogenetic Analysis by Maximum Likelihood) [16, 17] was used and the codon frequency model F3x4 was fitted to both alignments. In the site-specific models that allow the ratio ω to vary among codons, we performed Likelihood Ratio Tests (LRTs) with 2 degrees of freedom to compare the following models (NS sites): M1 (nearly neutral) with M2 (selection) and M7 (neutral, β distribution of ω < 1) with M8 (selection, β distribution of ω > 1). A significant LRT demonstrates that the selection model fits better than the neutral model [42, 43]. For model M8, a Bayes empirical Bayes (BEB) approach was employed to detect codons with a posterior probability >90% of being under selection . Also the branch-site model A was performed for testing positive selection on individual sites along a specific lineage, called foreground branch, where the other lineages are background branches. In branch-site model A, three ω ratios are assumed for foreground (0 < ω 0 < 1, ω 1 = 1, ω 2 > 1) and two ω ratios for background (0 < ω 0 < 1, ω 1 = 1). The null model is the same as model A, but ω 2 = 1 is fixed. We also used BEB approach to calculate the posterior probability of a specific codon site and to identify those most likely to be under positive selection (posterior probability >90%) .
Both SAMD9 and SAMD9L genes were also analyzed using HyPhy software implemented in the Datamonkey web server . Datamonkey includes three classic ML methods to detect sites under selection: the Single Likelihood Ancestor Counting (SLAC) model, the Fixed Effect Likelihood (FEL) model and the Random Effect Likelihood (REL) model . Besides these three methods, two other recently developed and implemented in the Datamonkey web server were applied to our dataset: the Mixed Effects Model of Evolution (MEME) that allows the distribution of ω to vary from site to site and also from branch to branch at a site, being capable of identifying both episodic and pervasive positive selection , and the Fast Unbiased Bayesian AppRoximation (FUBAR) method that can detect positive selection under a model faster than the existing fixed effects likelihood models through the introduction of an ultra-fast Markov chain Monte Carlo (MCMC) routine and that allows to visualize Bayesian inference for each site . All these methods were run using the best model chosen by AIC on a defined Neighbor-Joining (NJ) phylogenetic tree after running GARD to detect recombination. To avoid a high false-positive rate, due to the reduced number of sequences , sites with p-values <0.1 for SLAC, FEL and MEME models, Bayes Factor >50 for REL model and a posterior probability >0.90 for FUBAR were accepted as candidates for selection.
From the HyPhy software available on the Datamonkey web server, we also run the PARRIS method used to detect if a proportion of sites in the alignment evolve with d N /d S > 1 and that accounts for synonymous rate variation and recombination .
Amino acid-based analyses of positive selection
By using TreeSAAP v3.2 (Selection of Amino Acid Properties based on Phylogenetic Trees)  it was possible to detect selection signatures at the amino acid level, more specifically, positively selected amino acid properties that result in radical structural and functional changes in local regions of the protein (destabilization). Properties that fell into categories 6 through 8 (the most radical values denoting positive destabilizing selection), presented z-score values of 3.09 and higher, and with a probability value of 0.001 were plotted in a sliding window (length = 20).
Thirty-one amino acid properties were evaluated across SAMD9 and SAMD9L phylogenetic trees to identify protein regions that presented evidence of positive destabilization for each property. The thirty-one amino acid properties are the following: alpha-helical tendencies, average number of surrounding residues, beta-structure tendencies, bulkiness, buriedness, chromatographic index, coil tendencies, composition, compressibility, equilibrium constant (ionization of COOH), helical contact area, hydropathy, isoelectric point, long-range non-bonded energy, mean r.m.s. fluctuation displacement, molecular volume, molecular weight, normalized consensus hydrophobicity, partial specific volume, polar requirement, polarity, power to be at the C-terminal, power to be at the middle of alpha-helix, power to be at the N-terminal, refractive index, short and medium range non-bonded energy, solvent accessible reduction ratio, surrounding hydrophobicity, thermodynamic transfer hydrophobicity, total non-bonded energy and turn tendencies.
Giant panda - Ailuropoda melanoleuca
Cow - Bos taurus
Common marmoset - Callithrix jacchus
Domestic dog - Canis lupus familiaris
Domestic Guinea pig - Cavia porcellus
Hoffmann’s two-toed sloth - Choloepus hoffmanni
Chinese hamster - Cricetulus griseus
Nine-banded armadillo - Dasypus novemcinctus
Ord’s kangaroo rat - Dipodomys ordii
Lesser hedgehog tenrec - Echinops telfairi
Horse - Equus caballus
West European hedgehog - Erinaceus europaeus
Domestic cat - Felis catus
Western gorilla - Gorilla gorilla
Human - Homo sapiens
Thirteen-lined ground squirrel - Ictidomys tridecemlineatus
African bush elephant - Loxodonta africana
Rhesus monkey - Macaca mulatta
Grey short-tailed opossum - Monodelphis domestica
House mouse - Mus musculus
Little brown myotis - Myotis lucifugus
Northern white-cheeked gibbon - Nomascus leucogenys
American pika - Ochotona princeps
European rabbit - Oryctolagus cuniculus
Northern greater galago - Otolemur garnettii
Common chimpanzee - Pan troglodytes
Sumatran orangutan - Pongo abelii
Rock hyrax - Procavia capensis
Large flying fox - Pteropus vampyrus
Brown rat - Rattus norvegicus
Common shrew - Sorex araneus
Pig - Sus scrofa
Philippine tarsier - Tarsius syrichta
Alpaca - Vicugna pacos.
The Portuguese Foundation for Science and Technology supported the doctoral fellowship of ALM (SFRH/BD/48566/2008). A research project from the Portuguese Foundation for Science and Technology (PTDC/BIA-BEC/103158/2008) also supported the study. This work was also supported by grant R01 AI080607 from the National Institute of Health to GM. This research has also been assisted by the BIT Core of the University of California, San Diego, Center for AIDS Research (NIH P30 AI036214).
- Li CF, MacDonald JR, Wei RY, Ray J, Lau K, Kandel C, Koffman R, Bell S, Scherer SW, Alman BA: Human sterile alpha motif domain 9, a novel gene identified as down-regulated in aggressive fibromatosis, is absent in the mouse. BMC Genomics. 2007, 8: 92-10.1186/1471-2164-8-92.PubMed CentralPubMedView ArticleGoogle Scholar
- Asou H, Matsui H, Ozaki Y, Nagamachi A, Nakamura M, Aki D, Inaba T: Identification of a common microdeletion cluster in 7q21.3 subband among patients with myeloid leukemia and myelodysplastic syndrome. Biochem Biophys Res Commun. 2009, 383: 245-251. 10.1016/j.bbrc.2009.04.004.PubMedView ArticleGoogle Scholar
- Topaz O, Indelman M, Chefetz I, Geiger D, Metzker A, Altschuler Y, Choder M, Bercovich D, Uitto J, Bergman R, et al: A deleterious mutation in SAMD9 causes normophosphatemic familial tumoral calcinosis. Am J Hum Genet. 2006, 79: 759-764. 10.1086/508069.PubMed CentralPubMedView ArticleGoogle Scholar
- Chefetz I, Ben Amitai D, Browning S, Skorecki K, Adir N, Thomas MG, Kogleck L, Topaz O, Indelman M, Uitto J, et al: Normophosphatemic familial tumoral calcinosis is caused by deleterious mutations in SAMD9, encoding a TNF-alpha responsive protein. J Invest Dermatol. 2008, 128: 1423-1429. 10.1038/sj.jid.5701203.PubMed CentralPubMedView ArticleGoogle Scholar
- Tanaka M, Shimbo T, Kikuchi Y, Matsuda M, Kaneda Y: Sterile alpha motif containing domain 9 is involved in death signaling of malignant glioma treated with inactivated Sendai virus particle (HVJ-E) or type I interferon. Int J Cancer. 2010, 126: 1982-1991.PubMedView ArticleGoogle Scholar
- Hershkovitz D, Gross Y, Nahum S, Yehezkel S, Sarig O, Uitto J, Sprecher E: Functional characterization of SAMD9, a protein deficient in normophosphatemic familial tumoral calcinosis. J Invest Dermatol. 2011, 131: 662-669. 10.1038/jid.2010.387.PubMed CentralPubMedView ArticleGoogle Scholar
- Schoggins JW, Wilson SJ, Panis M, Murphy MY, Jones CT, Bieniasz P, Rice CM: A diverse range of gene products are effectors of the type I interferon antiviral response. Nature. 2011, 472: 481-485. 10.1038/nature09907.PubMed CentralPubMedView ArticleGoogle Scholar
- Liu J, Wennier S, Zhang L, McFadden G: M062 is a host range factor essential for myxoma virus pathogenesis and functions as an antagonist of host SAMD9 in human cells. J Virol. 2011, 85: 3270-3282. 10.1128/JVI.02243-10.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhang LK, Chai F, Li HY, Xiao G, Guo L: Identification of Host Proteins Involved in Japanese Encephalitis Virus Infection by Quantitative Proteomics Analysis. J Proteome Res. 2013, 12 (6): 2666-2678. 10.1021/pr400011k. Epub 2013 May 21PubMedView ArticleGoogle Scholar
- Pappas DJ, Coppola G, Gabatto PA, Gao F, Geschwind DH, Oksenberg JR, Baranzini SE: Longitudinal system-based analysis of transcriptional responses to type I interferons. Physiol Genomics. 2009, 38: 362-371. 10.1152/physiolgenomics.00058.2009.PubMedView ArticleGoogle Scholar
- Jiang Q, Quaynor B, Sun A, Li Q, Matsui H, Honda H, Inaba T, Sprecher E, Uitto J: The Samd9L gene: transcriptional regulation and tissue-specific expression in mouse development. J Invest Dermatol. 2011, 131: 1428-1434. 10.1038/jid.2011.61.PubMed CentralPubMedView ArticleGoogle Scholar
- Critchley-Thorne RJ, Yan N, Nacu S, Weber J, Holmes SP, Lee PP: Down-regulation of the interferon signaling pathway in T lymphocytes from patients with metastatic melanoma. PLoS Med. 2007, 4: e176-10.1371/journal.pmed.0040176.PubMed CentralPubMedView ArticleGoogle Scholar
- Kosakovsky Pond SL, Posada D, Gravenor MB, Woelk CH, Frost SD: GARD: a genetic algorithm for recombination detection. Bioinformatics. 2006, 22: 3096-3098. 10.1093/bioinformatics/btl474.PubMedView ArticleGoogle Scholar
- Kosakovsky Pond SL, Posada D, Gravenor MB, Woelk CH, Frost SD: Automated phylogenetic detection of recombination using a genetic algorithm. Mol Biol Evol. 2006, 23: 1891-1901. 10.1093/molbev/msl051.PubMedView ArticleGoogle Scholar
- Song S, Liu L, Edwards SV, Wu S: Resolving conflict in eutherian mammal phylogeny using phylogenomics and the multispecies coalescent model. Proc Natl Acad Sci USA. 2012, 109: 14942-14947. 10.1073/pnas.1211733109.PubMed CentralPubMedView ArticleGoogle Scholar
- Yang Z: PAML: a program package for phylogenetic analysis by maximum likelihood. Comput Appl Biosci. 1997, 13: 555-556.PubMedGoogle Scholar
- Yang Z: PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24: 1586-1591. 10.1093/molbev/msm088.PubMedView ArticleGoogle Scholar
- Scheffler K, Martin DP, Seoighe C: Robust inference of positive selection from recombining coding sequences. Bioinformatics. 2006, 22: 2493-2499. 10.1093/bioinformatics/btl427.PubMedView ArticleGoogle Scholar
- Wlasiuk G, Nachman MW: Adaptation and constraint at Toll-like receptors in primates. Mol Biol Evol. 2010, 27: 2172-2186. 10.1093/molbev/msq104.PubMed CentralPubMedView ArticleGoogle Scholar
- Areal H, Abrantes J, Esteves PJ: Signatures of positive selection in Toll-like receptor (TLR) genes in mammals. BMC Evol Biol. 2011, 11: 368-10.1186/1471-2148-11-368.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhang J: Rates of conservative and radical nonsynonymous nucleotide substitutions in mammalian nuclear genes. J Mol Evol. 2000, 50: 56-68.PubMedGoogle Scholar
- Woolley S, Johnson J, Smith MJ, Crandall KA, McClellan DA: TreeSAAP: selection on amino acid properties using phylogenetic trees. Bioinformatics. 2003, 19: 671-672. 10.1093/bioinformatics/btg043.PubMedView ArticleGoogle Scholar
- Adkins RM, Gelke EL, Rowe D, Honeycutt RL: Molecular phylogeny and divergence time estimates for major rodent groups: evidence from multiple genes. Mol Biol Evol. 2001, 18: 777-791. 10.1093/oxfordjournals.molbev.a003860.PubMedView ArticleGoogle Scholar
- Churakov G, Sadasivuni MK, Rosenbloom KR, Huchon D, Brosius J, Schmitz J: Rodent evolution: back to the root. Mol Biol Evol. 2010, 27: 1315-1326. 10.1093/molbev/msq019.PubMedView ArticleGoogle Scholar
- Bininda-Emonds OR, Cardillo M, Jones KE, MacPhee RD, Beck RM, Grenyer R, Price SA, Vos RA, Gittleman JL, Purvis A: The delayed rise of present-day mammals. Nature. 2007, 446: 507-512. 10.1038/nature05634.PubMedView ArticleGoogle Scholar
- Hughes AL, Ota T, Nei M: Positive Darwinian selection promotes charge profile diversity in the antigen-binding cleft of class I major-histocompatibility-complex molecules. Mol Biol Evol. 1990, 7: 515-524.PubMedGoogle Scholar
- Zhang J, Webb DM: Rapid evolution of primate antiviral enzyme APOBEC3G. Hum Mol Genet. 2004, 13: 1785-1791. 10.1093/hmg/ddh183.PubMedView ArticleGoogle Scholar
- Ponting CP: SAM: a novel motif in yeast sterile and Drosophila polyhomeotic proteins. Protein Sci. 1995, 4: 1928-1930. 10.1002/pro.5560040927.PubMed CentralPubMedView ArticleGoogle Scholar
- Qiao F, Bowie JU: The many faces of SAM. Sci STKE. 2005, 2005: re7-PubMedGoogle Scholar
- Meruelo AD, Bowie JU: Identifying polymer-forming SAM domains. Proteins. 2009, 74: 1-5. 10.1002/prot.22232.PubMed CentralPubMedView ArticleGoogle Scholar
- Knight MJ, Leettola C, Gingery M, Li H, Bowie JU: A human sterile alpha motif domain polymerizome. Protein Sci. 2011, 20: 1697-1706. 10.1002/pro.703.PubMed CentralPubMedView ArticleGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.PubMed CentralPubMedView ArticleGoogle Scholar
- Hall T: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser. 1999, 41: 95-98.Google Scholar
- Xia X, Xie Z: DAMBE: software package for data analysis in molecular biology and evolution. J Hered. 2001, 92: 371-373. 10.1093/jhered/92.4.371.PubMedView ArticleGoogle Scholar
- Posada D, Crandall KA: The effect of recombination on the accuracy of phylogeny estimation. J Mol Evol. 2002, 54: 396-402.PubMedView ArticleGoogle Scholar
- Pond SL, Frost SD: Datamonkey: rapid detection of selective pressure on individual sites of codon alignments. Bioinformatics. 2005, 21: 2531-2533. 10.1093/bioinformatics/bti320.PubMedView ArticleGoogle Scholar
- Posada D: jModelTest: phylogenetic model averaging. Mol Biol Evol. 2008, 25: 1253-1256. 10.1093/molbev/msn083.PubMedView ArticleGoogle Scholar
- Zwickl DJ: Genetic algorithm approaches for the phylogenetic analysis of large biological sequence datasets under the maximum likelihood criterion. 2006, University of Texas: PhD ThesisGoogle Scholar
- Miyata T, Yasunaga T: Molecular evolution of mRNA: a method for estimating evolutionary rates of synonymous and amino acid substitutions from homologous nucleotide sequences and its application. J Mol Evol. 1980, 16: 23-36. 10.1007/BF01732067.PubMedView ArticleGoogle Scholar
- Nei M, Gojobori T: Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. Mol Biol Evol. 1986, 3: 418-426.PubMedGoogle Scholar
- Crandall KA, Kelsey CR, Imamichi H, Lane HC, Salzman NP: Parallel evolution of drug resistance in HIV: failure of nonsynonymous/synonymous substitution rate ratio to detect selection. Mol Biol Evol. 1999, 16: 372-382. 10.1093/oxfordjournals.molbev.a026118.PubMedView ArticleGoogle Scholar
- Nielsen R, Yang Z: Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene. Genetics. 1998, 148: 929-936.PubMed CentralPubMedGoogle Scholar
- Yang Z, Nielsen R, Goldman N, Pedersen AM: Codon-substitution models for heterogeneous selection pressure at amino acid sites. Genetics. 2000, 155: 431-449.PubMed CentralPubMedGoogle Scholar
- Yang Z, Wong WS, Nielsen R: Bayes empirical bayes inference of amino acid sites under positive selection. Mol Biol Evol. 2005, 22: 1107-1118. 10.1093/molbev/msi097.PubMedView ArticleGoogle Scholar
- Kosakovsky Pond SL, Frost SD: Not so different after all: a comparison of methods for detecting amino acid sites under selection. Mol Biol Evol. 2005, 22: 1208-1222. 10.1093/molbev/msi105.PubMedView ArticleGoogle Scholar
- Murrell B, Wertheim JO, Moola S, Weighill T, Scheffler K, Kosakovsky Pond SL: Detecting individual sites subject to episodic diversifying selection. PLoS Genet. 2012, 8: e1002764-10.1371/journal.pgen.1002764.PubMed CentralPubMedView ArticleGoogle Scholar
- Murrell B, Moola S, Mabona A, Weighill T, Sheward D, Kosakovsky Pond SL, Scheffler K: FUBAR: A Fast, Unconstrained Bayesian AppRoximation for Inferring Selection. Mol Biol Evol. 2013, 30: 1196-1205. 10.1093/molbev/mst030.PubMed CentralPubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.