Evolution and functional divergence of NLRPgenes in mammalian reproductive systems
© Tian et al; licensee BioMed Central Ltd. 2009
Received: 8 April 2009
Accepted: 14 August 2009
Published: 14 August 2009
NLRPs (Nucleotide-binding oligomerization domain, Leucine rich Repeat and Pyrin domain containing Proteins) are members of NLR (Nod-like receptors) protein family. Recent researches have shown that NLRP genes play important roles in both mammalian innate immune system and reproductive system. Several of NLRP genes were shown to be specifically expressed in the oocyte in mammals. The aim of the present work was to study how these genes evolved and diverged after their duplication, as well as whether natural selection played a role during their evolution.
By using in silico methods, we have evaluated the evolution and functional divergence of NLRP genes, in particular of mouse reproduction-related Nlrp genes. We found that (1) major NLRP genes have been duplicated before the divergence of mammals, with certain lineage-specific duplications in primates (NLRP7 and 11) and in rodents (Nlrp1, 4 and 9 duplicates); (2) tandem duplication events gave rise to a mammalian reproduction-related NLRP cluster including NLRP2, 4, 5, 7, 8, 9, 11, 13 and 14 genes; (3) the function of mammalian oocyte-specific NLRP genes (NLRP4, 5, 9 and 14) might have diverged during gene evolution; (4) recent segmental duplications concerning Nlrp4 copies and vomeronasal 1 receptor encoding genes (V1r) have been undertaken in the mouse; and (5) duplicates of Nlrp4 and 9 in the mouse might have been subjected to adaptive evolution.
In conclusion, this study brings us novel information on the evolution of mammalian reproduction-related NLRPs. On the one hand, NLRP genes duplicated and functionally diversified in mammalian reproductive systems (such as NLRP4, 5, 9 and 14). On the other hand, during evolution, different lineages adapted to develop their own NLRP genes, particularly in reproductive function (such as the specific expansion of Nlrp4 and Nlrp9 in the mouse).
The innate immune system is an ancestral and ubiquitous system of defense against microbial infection and other potential threats to the host. The first mammalian molecules shown to be involved in innate immune recognition of, and defense against, microbial pathogens are the Toll-like receptors (TLRs), which constitute the main sensors for detection of extracellular microbes . Recent research has shown that Nod-like receptors (NLRs) act as intracellular regulators of bacterial-induced inflammation [2–4].
The NLR protein family is separated into several subfamilies based on the different N-terminal effecter domains . The NLRP subfamily (known also as NALP family) is a new identified NLR group characterized by a PYRIN domain. Initially, research on the function of NLRP proteins focused on their roles in apoptotic and inflammatory signaling pathways via the formation of a large signaling platform (named inflammasome) and the activation of caspases in innate immunity [4, 6, 7]. Interestingly, more and more researches have recently revealed that some NLRP genes play roles in mammalian reproduction. For instance, the mouse Nlrp5 (known also as MATER) was one of the first identified mammalian maternal effect genes, i.e. it encodes mRNA required for successful development of a fertilized oocyte . In particular, the targeted invalidation of mouse oocyte-specific Nlrp5 leads to the sterility of females due to a blockage at the two-cell stage in the development of embryos . In vitro knockdown experiments in mouse fertilized eggs also revealed that a decrease of mouse germ-cell-specific Nlrp14 transcripts results in an arrest of development between the one-cell and eight-cell stages of the embryos [9, 10]. In the mouse, exclusive duplications of Nlrp4 and Nlrp9 have been detected with specific expression profiles restricted to the oocyte , whereas the mouse Nlrp4a and the cattle NLRP9 have also been reported to be expressed in the testis [11, 12]. In humans, the mutations of NLRP7 are found to cause recurrent hydatidiform moles, spontaneous abortions, stillbirths and intrauterine growth retardation [13, 14]. Additionally, the recent expression analyses of NLRP genes in the human and the rhesus macaque (Macaca mulatta) have shown that most if not all NLRP genes are expressed in primate gametes and early embryos, suggesting a role of NLRPs in primate pre-implantation development [15, 16].
Overall, it seems that the NLRP genes play roles not only in innate immunity but also in the reproductive system of mammals. Moreover, Nalp5 and Nalp14 invalidation data in the mouse suggests that there is no functional compensation of each of these oocyte-specific genes. In this context, one question is to study how these reproduction-related NLRP genes evolved and diverged in function after their duplication, as well as whether natural selection played a role during their evolution.
The availability of several completely sequenced vertebrate genomes allows us to use a phylogenomic approach to identify the NLRP orthologues from different mammals and to evaluate the evolutionary features of NLRPs in different lineages. In this paper, we focus on the evolution and functional divergence of reproduction-related NLRPs, as well as lineage-specific expansion of Nlrp4 and Nlrp9 in rodents. The aim of the present work is to provide new information in understanding how NLRP genes evolve in reproductive system and how different mammals adapted to develop their own NLRP gene copies particularly in reproductive biology.
Phylogeny and syntenic comparison of NLRPs
A total of 83 mammalian NLRP amino acid sequences with chicken NLRP3-like as outgroup were used for the reconstruction of phylogenetic trees. A total of 122 positions are included in the final dataset.
In rodents, three Nlrp genes show lineage-specific duplications including Nlrp1, Nlrp4 and Nlrp9 (Figure 1). More specifically, Nlrp1 exclusively duplicated in the mouse, Nlrp9 duplicated independently in the rodents after the separation of mouse and rat, whereas the duplication of Nlrp4 is likely to have occurred both before and after this separation.
Unexpectedly, NLRP8 and NLRP13 seem to be universal in mammals except for rodents. We speculate that these two genes might have been lost during evolution of the rodents. Unfortunately, our tBLASTn search (using the human NLRP8 and 13 respectively as queries against the mouse genome) failed to find any trace of pseudogene of neither NLRP8 nor NLRP13 within the mouse conserved syntenic region. We found that, on the human chromosome 19, these two genes are located side by side between NLRP4 and NLRP5, whereas the mouse Nlrp4 gene has been expansively duplicated and the duplicates are located around Nlrp5 on the mouse chromosome 7. In rodents, the biological significance of the probable loss of NLRP8 and 13 on one hand, and the parallel expansion of Nlrp4 on the other hand, remains unknown.
Remarkably, our phylogenetic analyses identified a well supported (bootstrap values >90% by ME, NJ and ML methods) reproduction-related clade including nine NLRP proteins: NLRP2, 4, 5, 7, 8, 9, 11, 13 and 14 (Figure 1). All these NLRP genes have been shown to be expressed in the human oocyte and pre-implantation embryos . NLRP4, 5, 8, 9 and 14 are exclusively expressed in germ-cells, especially in the oocyte of mammals [8–12, 17, 18]. The mutations of NLRP7 are associated with abnormal embryo development in the human [13, 14]. The expression of mouse Nlrp2 has been found to be highly enriched in fully grown oocytes but diminished in the 2-cell embryos upon embryo genome activation . In both the human  and the rhesus macaque , NLRP11 and 13 have been found specifically or preferentially expressed in the oocyte, with a similar expression pattern to other oocyte-specific NLRP genes, i.e., enriched in maturing oocytes and then progressively diminished in embryos, indicating an exclusive maternal origin of these transcripts. Unfortunately, the biological study of these two genes is not sufficient to address their putative roles in reproduction.
We identified the other NLRP members (NLRP1, 3, 6, 10 and 12) as non-reproduction-related NLRP genes, because they are not in the reproduction-related clade in the phylogenetic trees and ubiquitously expressed. These genes have been shown to participate in inflammatory and immune responses by regulating the activation of other signaling factors, such as NF-κB, caspase-1 and IL-1β [20–24]. In contrast, little is known about the function of the reproduction-related NLRP proteins in inflammation and immunity, except that the human NLRP2, 4 and 7 may be able to inhibit IL-1β and/or NF-κB [25–28], and thus may contribute to modulate the inflammatory response. It is known that inflammation and bacterial infection can cause infertility, ectopic pregnancy and abortion . Thus, NLRP proteins might play roles both in innate immunity and in reproductive biology.
In the reproduction-related NLRP cluster referred by the phylogenetic analysis above, there are several oocyte-specific NLRP genes including NLRP4, 5, 8, 9 and 14. To investigate whether functional divergence have occurred following the duplications of these oocyte-specific NLRPs, four oocyte-expressed NLRP genes (NLRP4, 5, 9 and 14) were estimated for functional divergence by using the DIVERGE program . NLRP8, which was not found in rodents, was therefore not included in this analysis in order to compare the level of divergence across all candidate species.
Pairwise comparison of functional divergence of mouse oocyte-expressed Nlrp genes
0.65 ± 0.07
0.61 ± 0.05
0.96 ± 0.08
0.89 ± 0.06
0.68 ± 0.08
0.67 ± 0.06
0.76 ± 0.08
0.65 ± 0.06
0.64 ± 0.09
0.58 ± 0.07
0.49 ± 0.10
0.45 ± 0.07
Furthermore, the important amino acid residues, responsible for functional divergence, were predicted by calculating the site-specific profile based on posterior analysis for all pairs of clusters with functional divergence. In order to extensively reduce false positive, cutoff values, i.e., the minimal posterior probabilities for RFD (Residues with predicted Functional Divergence) were established empirically by progressively removing the highest scoring residues from the alignment until θ dropped to zero. As shown in Table 1, the least RFD (33 residues, covering 5.1% of a total of 650 aligned sites) are observed between NLRP5 and NLRP14, and the most RFD (118 residues, covering 18.2% of a total of 650 aligned sites) are observed between NLRP4 and NLRP14. In general, the RFD are detected in all three functional domains (NACHT, PYRIN and LRRs) of NLRP proteins, implying that shifted functional constrains might have acted on each protein domain.
Segmental duplications concerning Nlrp4 and V1rgenes in the mouse genome
In the mouse genome, Nlrp4 duplicates (except for Nlrp4f and 4g) are located between 6.00 Mb and 27.32 Mb on the chromosome 7. Interestingly, all these Nlrp4 duplicates are found next to the vomeronasal 1 receptor encoding genes (V1r), such as V1rd and V1re genes (Figure 3a). The phylogeny and genomic location of these V1r genes reveal that V1r genes might have a similar duplication model to their neighbors, Nlrp4 genes (Figure 3a). Thus, one may hypothesize that Nlrp4 may have not duplicated one by one, but were rather involved in segmental duplications encompassing several different genes, such as V1r genes.
Vomeronasal organ, which is well developed in rodents in comparison with human, detects pheromones and other chemical signals implicated in innate reproductive and social behaviours between the members of the same species [31, 32]. The V1r gene family is supposed to have emerged during the divergence of placental mammals (80-110 MYA), but many subfamilies including V1rd and V1re originated only in rodents after their separation from the primates , which is consistent with the rodent-specific duplication of Nlrp4. Thus, the phylogenetic relationship and genomic location (Figure 3a) of these V1r genes, that are located in the segmental duplication region, could help us to better understand the segmental duplications. Overall, as shown in Figure 3b, a presumed itinerary of segmental duplication was deduced by combining phylogeny with genomic location of both Nlrp4 duplicates and V1r genes. Although we can not confirm the order of duplication events by the present data, we can state that segmental duplication events characterize the evolution of this genomic region.
Positive selection on Nlrp4 and Nlrp9paralogues in the mouse
Several examples have shown that genes in duplication blocks could be maintained by positive selection [34–37]. Here, we implemented PAML4  to investigate the model of selective force acting on mouse Nlrp4 and Nlrp9 duplicates during their evolution.
Parameter estimates and likelihood scores for site models in PAML
Estimates of parameters
Positively selected sites (BEB)
ρ0 = 0.33
ρ0 = 0.28,(ρ1 = 0.67), ρs = 0.05
ω s = 3.93
(M2a vs M1a)
2 sites > 95%
p = 0.04, q = 0.02
ρ0 = 0.34,(ρ1= 0.66),
p = 0.81, q = 6.69
(M8 vs M7)
ρ0 = 0.94,ρs = 0.06,
p = 0.20, q = 0.09,
ω s = 3.64
(M8 vs M8a)
4 sites > 95%
133C, 153S, 657N, 677V
ρ0 = 0.36043
ρ0 = 0.28, ρ1 = 0.66, ρs = 0.06,
ω 2 = 10.17
(M2a vs M1a)
10 sites > 95%, 2 sites > 99%
106H, 123 E, 356F, 360L, 378R, 513E, 536V, 571V, 661L, 948A
p = 0.02, q = 0.01
ρ0 = 0.36,(ρ1 = 0.64),
p = 0.01, q = 98.97
(M8 vs M7)
ρ0 = 0.94, ρs = 0.06,
p = 0.01, q = 0.01,
ω s = 10.20
(M8 vs M8a)
11 sites > 95%, 6 sites > 99%
106H, 123 E, 278L, 356F, 360L, 378R, 513E, 536V, 571V, 661L, 948A
Subsequently, we wanted to know if the target of positive selection had changed during evolution in these two gene clades. We compared the potential sites subjected to positive selection between the two datasets (Table 2). In Nlrp4_mouse dataset, 2 sites (133C and 677V) are identified as positively selected sites at p > 95% level by both M2a and M8 models, and 2 more sites (153S and 657N) are identified only by the M8 model. In Nlrp9_mouse dataset, 10 sites (106H, 123 E, 356F, 360L, 378R, 513E, 536V, 571V, 661L and 948A) are identified as potential targets of positive selection by both M2a and M8 models, and 1 more sites (278L) is identified only by the M8 model. We note that no common sites are shared by the two datasets as positively selected sites at p > 95% level, suggesting a dramatic shift in the target of positive selection between the Nlrp4 and Nlrp9 genes in the mouse. When the 3D structures or models of these proteins are available in the future, we could investigate whether these positively selected sites are located in functionally pivotal regions.
Independent evolution of reproduction-related NLRPgenes in different mammalian lineages
Phylogenetic analyses show that NLRP genes have been originated and duplicated before the divergence of mammals. Certain genes are well conserved during their evolution, such as NLRP3 (82% identity between the human and the mouse), whereas many genes involved in reproduction have rapidly evolved resulting in higher sequence divergence among lineages, such as NLRP5 (50% identity between the human and the mouse) and NLRP14 (62% identity between the human and the mouse). The pairwise estimates of evolutionary divergence between the reproduction-related paralogues are ranged from 61% to 75% (except that the pairwise distance between NLRP2 and its recent primate-specific duplicate, NLRP7, is 45%), which is obviously higher than the one between the non-reproduction-related paralogues (ranged from 43% to 69%). Moreover, gene duplication and/or gene loss is found to have occurred independently in different mammalian lineages. On the one hand, Nlrp4 (originated before the divergence of primates and rodents) and Nlrp9 have been extensively duplicated in rodents, and other lineage-specific gene duplication events concern the specific duplication of Nlrp1 (three paralogues) in the mouse as well as the unique origin of NLRP7 and NLRP11 in primates. On the other hand, Nlrp8 and Nlrp13 are seemed to be lost from the genomes of rodents. Interestingly, the major gene duplication and gene loss events are found to be associated with the reproduction-related genes (such as NLRP4, 7, 8, 9, 11 and 13), implying that reproduction-related genes have undertaken a fast and diverged evolution among different mammalian lineages. Given note that major reproduction-related NLRP genes are germ-cell specific, whereas the other non-reproduction-related genes are usually expressed in multiple tissues from a single organism. Recent large-scale gene expression studies have shown that the tissue specificity of genes is correlated positively with gene evolution rates [39–41]. Furthermore, our result reinforces the hypothesis that reproduction-related genes are most highly divergent and evolve more rapidly than genes expressed in other organs [42–45].
Functional divergence of the oocyte-specific genes
Gene duplication is thought to be the essential source of gene novelty, with new or altered functions, as shown by widespread existence of gene families. Among NLRP genes, there are at least four oocyte specifically expressed in mammals, including NLRP4, 5, 9 and 14. These genes have been found to be restricted expressed in oocytes and early embryos [9, 15], suggesting their important roles in oogenesis and preimplantation embryo development. Additionally, their temporal expression in oocytes and degradation during preimplantation development coincide with the timing of gene expression transition from maternal to zygotic control, indicating their roles as maternal effect genes.
In the present work, by using an in silico method (with the help of the program DIVERGE), we demonstrate that the functional divergence could have occurred between each pair of oocyte-specific NLRP genes. This result is supported by certain published experimental data. For example, the knock-out of mouse Nlrp5 induced female infertility due to a blockage of early embryonic development ; and the injection of siRNA against Nlrp14 into fertilized eggs resulted in arrested development of embryos between 1-cell and 8-cell stages . These results suggest that the expression of other oocyte-specific NLRP genes is not able to compensate the absence of neither NLRP5 nor NLRP14 genes. Additionally, recent research has shown that NLRP14 is also expressed in the testis and its mutation might cause spermatogenic failure in the human . Thus, we presume that the expressional divergence of NLRP5 and NLRP14 during the development of mouse germ-cells and embryos may lead to functional specialization. On the other hand, target invalidation/inhibition of NLRP5 and 14 did not reveal any deleterious effect on ovarian folliculogenesis but rather in early embryonic development, suggesting other genes might play the similar roles during folliculogenesis and meiotic maturation. For these latter functions, plausible candidates could be other oocyte-specific NLRP genes, such as NLRP4 and NLRP9. To investigate this possibility, further functional studies such as knock-out models or other targeted inhibition experiments on these genes should be carried out.
Gene duplication is thought to be a major driving force in enabling the evolution of tissue specialization. By using microarray gene expression data from mammals , fruit flies and yeast , it has recently shown that as gene family size enlarges, there is a general trend for paralogous genes with decreased breadth and increased specificity of expression. In particular, by studying the relationship between gene family size and expression breadth of 1249 protein families in the mouse, Freilich and his colleagues  have recently demonstrated that duplicates that arose through post-multicellularity duplication events show a tendency to become more specifically expressed, supporting the view suggested by the subfunctionalization model . In our case, especially with the evidence that NLRP5 and NLRP14 genes are significantly involved in different developmental stages of embryos, we hypothesize that the functional divergence of oocyte-specific NLRP duplicates might be derived by the expression specialization of duplicate genes, which could provide the mammals reproductive advantages such as in adaptation.
Segmental duplication in the mouse genome
The phylogenetic research indicates the rodent-specific expansion of Nlrp4 and NLRP9, and the furthermore genomic analyses find that these Nlrp duplicates, especially Nlrp4, have been duplicated together with other genes, such as V1r copies. This segmental duplication is not rare during the mouse and human genomic evolution. It has recently revealed that the segmental duplication constitutes about 5% of mouse genome, and its distribution is in a highly non-random fashion . Interestingly, the duplication blocks account for 32% of the first 50 Mb of the chromosome 7, where the segmental duplications concerning Nlrp4 and V1r genes are located. Nlrp9 has been also duplicated in this region, but its duplication scenario is not clear and might be due to subsequent recombination.
As found in the human , the segmental duplication regions in the mouse might also be "hot spots" for the occurrence of non-allelic homologous recombination, leading to genomic mutations such as deletion, duplication, inversion or translocation [53–55]. Thus, these instable regions might have a biological significance in form of genome evolution [56, 57].
Adaptive diversification after gene duplications in Nlrp4 and Nlrp9
In the mouse, Nlrp4 and Nlrp9 are specifically extensive in gene copies. By using Site-Models implemented in PAML4, we evaluated the variation of selective pressure acting on Nlrp4 and Nlrp9 duplicates, respectively. Our result shows that positive selection is significant in both datasets, which consistent with other reports referring that the rapid divergence of the reproductive genes may be driven by positive selection [43, 58]. In particular, in the mouse, several gene cluster concerning lineage-specific expansion of reproduction-related genes have been found to be under adaptive evolution (positive selection), such as the Sva (seminal vesicle autoantigen) gene cluster , the Prame (or oogenesin) gene cluster , the Psg (pregnancy-specific glycoprotein) gene cluster , the Rhox4 gene cluster  and the Svs (seminal vesicle secretion proteins) gene cluster . The lineage-independent expansion and subsequently rapid evolution of such genes might contribute to speciation or adaptation, such as genetic barriers between species and hybrid incompatibilities, or provide the species specification in reproductive processes such as sperm competition, host immunity to pathogens, and manipulation of female/male reproductive physiology and behavior. In the case of NLRP4 and 9 genes, their biological functions are still unclear, so we can not evaluate the significance of their specific duplications in the reproduction of the mouse. However, the similar evolutionary mechanism driven by positive selection suggests that after duplication, these NLRP gene copies diverged and acquired abilities (probably by subfunctionalization) to adapt the new environment. In the subsequent research, the adaptive evolution of the other reproduction-related NLRP genes could be further evaluated. Hopefully, when the biological functions of these NLRP genes are elucidated, the nature of the selective pressure acting on NLRPs will be better understood.
Interestingly, the evolution of V1r genes, the other genes involved in the segmental duplication with Nlrp4 and 9, is also driven by positive selection in rodents . V1R proteins are thought to be responsible for the detection of pheromones that induce innate reproductive behaviors between members of the same species [31, 32]. Thus, it is presumed that the adaptive evolution of the V1r genes might play an important role in reinforcing pre-zygotic barriers among species of rodents .
NLRP genes have originated and duplicated mainly before the divergence of mammals. During evolution, NLRPs, especially the reproduction-related NLRP genes, have undergone a fast and independent diversification in different mammalian lineages. The expansion of reproduction-related NLRP genes has been associated with functional divergence after duplication, suggesting that each NLRP oocyte-specific gene might play an essential role in oogenesis and early preimplantation embryo development. The mouse-specific expansion of Nlrp4 and Nlrp9, concerning segmental duplication events, has been driven by positive evolution. Our founding suggests that the duplication and functional divergence of NLRP might provide mammals advantages in reproductive biology.
Molecular phylogenetic and syntenic analyses
The protein sequences of all the 14 known human NLRPs were retrieved from the GenBank http://www.ncbi.nlm.nih.gov/. The NLRP proteins from other species were searched by PSI-BLAST  with human NLRP protein sequences as queries against the protein databases (NCBI: RefSeq protein databases) of chimpanzee (Pan troglodytes), dog (Canis familiaris), cattle (Bos taurus), mouse (Mus musculus), rat (Rattus norvegicus) and chicken (Gallus gallus). The predicted coding sequences of the best hit proteins were retrieved when the hits presented more than 80% in length to be aligned with the query sequence (with E values < 10-100). These settings could effectively detect the potential NLRP members from different species but avoid involving the relative NLR proteins from other protein subfamilies. After removal of redundant sequences, the initial data set (Additional file 2) for NLRP phylogenetic studies includes 83 protein sequences from 6 mammals and 1 sequence from chicken.
Analyses of the orthologous and paralogous relationships among different species were carried out by combining the phylogenetic reconstruction of the gene family with the syntenic comparison. Multiple sequence alignments were performed using the Clustal W algorithm  and then manually edited. The alignment was submitted as Online Additional file 3. All alignment gap sites were eliminated before phylogenetic analyses. The phylogenetic trees were reconstructed with Neighbor Joining (NJ), Minimum-Evolution (ME) and Maximum Parsimony (MP) methods implemented in MEGA4 , as well as with Maximum Likelihood (ML) method implemented in PhyML . The consensus phylogenetic tree  was generated by the fusion of the independent trees reconstructed by the four methods. In all analyses, chicken NLRP protein (Q5F3J4) was treated as the root of all the mammalian NLRPs. The bootstrap values  were estimated with 1000 replications. The syntenic comparison is based on Ensembl http://www.ensembl.org/index.html utilities such as "orthologue prediction" and "view syntenic regions".
Test for functional divergence of oocyte-specific NLRPgenes
The program DIVERGE  was used to estimate type I functional divergence [70, 71] between oocyte-expressed NLRP paralogues (NLRP4, 5, 9 and 14). Type I sites represent amino acid residues conserved in one clade (gene cluster) but highly variable in another, suggesting that these residues have been subjected to different functional constraints. Statistically, this functional divergence between two clades is measured by the coefficient of functional divergence, θ, ranging from 0 to 1. A null hypothesis of θ = 0 indicates that the evolutionary rate is virtually the same between two duplicate genes at each site [70, 71]. If the null hypothesis was rejected, a site-specific profile was then used to predict the critical amino acid residues most likely to be responsible for the detected functional divergence.
A set of 33 protein sequences from mammalian NLRP4, 5, 9 and 14 is included in the analysis. Sequences of rat NLRP4d, 9a and 14 were excluded from the alignment because of their shorter lengths. Before the use of DIVERGE, our sequence dataset was examined to satisfy the conditions recommended by Gu et Vander Velden (2002), which permit us to obtain higher efficiency of detecting functional divergence-related residues: 1) each cluster has more than four amino acid sequences; (2) all pairwise sequence identities are <90%; and (3) multiple alignment is reliable. The phylogenetic tree used for DIVERGE was reconstructed by MEGA4  with ME method. The coefficients of functional divergence (θ) between gene clusters are calculated by Model-Free Estimation (MFE) and Maximum-Likelihood Estimation under two-state model (MLE). To detect amino acid residues reflecting functional divergence, all four reproduction-related NLRP gene clades were compared to each other.
Investigation of mouse specific segmental duplications
In order to investigate the predicted segmental duplications concerning Nlrp4 duplicates on the mouse chromosome 7, a tBLASTn search  was implemented by using the mouse Nlrp4a (the longest Nlrp4 sequence) as a query against the mouse genome. This method permitted us to identify both expressed Nlrp4 genes and traces of the pseudogenes. Then we combined phylogenetic analyses and genomic mapping to identify the segmental duplication . The phylogenetic tree of mouse Nlrp4 genes and pseudogenes was reconstructed by using nucleotide sequences instead of protein sequences. The phylogenetic tree was also reconstructed for vomeronasal 1 receptor proteins (V1R), because their encoding genes are located paralleling to Nlrp4 duplicates in the mouse genome and they are presumed to have duplicated together with Nlrp4. The accession numbers of the V1R protein sequences used here are listed in Additional file 4. The genomic organization of all the concerned genes and pseudogenes in the mouse is inferred from Ensembl (release 49).
To examine whether the duplicates of Nlrp4 and Nlrp9 in the mouse have been submitted to adaptive evolution, an analysis of variation in selective pressure following gene duplication events was carried out with the CODEML program implemented in PAML4 . The alignments were resulted from Clustal W and PAL2NAL . The shorter mouse Nlrp4g was excluded in order to obtain more informative sites.
Two datasets named Nlrp4_mouse (including 7 sequences) and Nlrp9_mouse (including 3 sequences) were investigated for different Site Models  implemented in PAML4. In this study, we employed three pairs of models including M1a (NearlyNeutral: 0< ω0 <1 and ω1 = 1) versus M2a (PositiveSelection: 0< ω0 <1, ω1 = 1 and ωs >1) , M7 (beta: 0< ω <1) versus M8 (beta and ω: 0< ω <1 and ωs >1) , and M8a (beta &ωs = 1: 0< ω <1 and ωs = 1) versus M8 . LRTs (Likelihood Ratio Tests) were used to test for significant differences in the fit of the models incorporating selection relative to their (nested) counterparts that did not allow positive selection . These tests provide a useful series of metrices for interpreting the significance of the results and a degree of protection against false positives . Bayes empirical Bayes (BEB) method  implemented in PAML4 was used to estimate posterior probabilities of selection on each codon. Furthermore, a newly described model MEC  was also employed on the Selecton Server  to compare with the results from other models for positive selection. The advantage of the MEC model over the other models used here is that by treating specifically each amino-acid replacement, Ka is computed differently. So under the MEC model, a position with radical replacements will obtain a higher Ka value than a position with more moderate replacements .
Xin Tian is funded by an INRA postdoctoral fellowship. We are grateful to John Williams for attention to the English language version, and to Pierre Pontarotti and Anthony Levasseur for helpful discussion.
- Barton GM, Medzhitov R: Toll-like receptors and their ligands. Current topics in microbiology and immunology. 2002, 270: 81-92.PubMedGoogle Scholar
- Chamaillard M, Girardin SE, Viala J, Philpott DJ: Nods, Nalps and Naip: intracellular regulators of bacterial-induced inflammation. Cellular microbiology. 2003, 5 (9): 581-592. 10.1046/j.1462-5822.2003.00304.x.View ArticlePubMedGoogle Scholar
- Fritz JH, Ferrero RL, Philpott DJ, Girardin SE: Nod-like proteins in immunity, inflammation and disease. Nature immunology. 2006, 7 (12): 1250-1257. 10.1038/ni1412.View ArticlePubMedGoogle Scholar
- Meylan E, Tschopp J, Karin M: Intracellular pattern recognition receptors in the host response. Nature. 2006, 442 (7098): 39-44. 10.1038/nature04946.View ArticlePubMedGoogle Scholar
- Ting JP, Lovering RC, Alnemri ES, Bertin J, Boss JM, Davis BK, Flavell RA, Girardin SE, Godzik A, Harton JA, et al: The NLR gene family: a standard nomenclature. Immunity. 2008, 28 (3): 285-287. 10.1016/j.immuni.2008.02.005.PubMed CentralView ArticlePubMedGoogle Scholar
- Martinon F, Gaide O, Petrilli V, Mayor A, Tschopp J: NALP inflammasomes: a central role in innate immunity. Seminars in immunopathology. 2007, 29 (3): 213-229. 10.1007/s00281-007-0079-y.View ArticlePubMedGoogle Scholar
- Tschopp J, Martinon F, Burns K: NALPs: a novel protein family involved in inflammation. Nat Rev Mol Cell Biol. 2003, 4 (2): 95-104. 10.1038/nrm1019.View ArticlePubMedGoogle Scholar
- Tong ZB, Gold L, Pfeifer KE, Dorward H, Lee E, Bondy CA, Dean J, Nelson LM: Mater, a maternal effect gene required for early embryonic development in mice. Nature genetics. 2000, 26 (3): 267-268. 10.1038/81547.View ArticlePubMedGoogle Scholar
- Hamatani T, Falco G, Carter MG, Akutsu H, Stagg CA, Sharov AA, Dudekula DB, VanBuren V, Ko MS: Age-associated alteration of gene expression patterns in mouse oocytes. Human molecular genetics. 2004, 13 (19): 2263-2278. 10.1093/hmg/ddh241.View ArticlePubMedGoogle Scholar
- Horikawa M, Kirkman NJ, Mayo KE, Mulders SM, Zhou J, Bondy CA, Hsu SY, King GJ, Adashi EY: The mouse germ-cell-specific leucine-rich repeat protein NALP14: a member of the NACHT nucleoside triphosphatase family. Biology of reproduction. 2005, 72 (4): 879-889. 10.1095/biolreprod.104.033753.View ArticlePubMedGoogle Scholar
- Dade S, Callebaut I, Paillisson A, Bontoux M, Dalbies-Tran R, Monget P: In silico identification and structural features of six new genes similar to MATER specifically expressed in the oocyte. Biochemical and biophysical research communications. 2004, 324 (2): 547-553. 10.1016/j.bbrc.2004.09.086.View ArticlePubMedGoogle Scholar
- Ponsuksili S, Brunner RM, Goldammer T, Kuhn C, Walz C, Chomdej S, Tesfaye D, Schellander K, Wimmers K, Schwerin M: Bovine NALP5, NALP8, and NALP9 genes: assignment to a QTL region and the expression in adult tissues, oocytes, and preimplantation embryos. Biology of reproduction. 2006, 74 (3): 577-584. 10.1095/biolreprod.105.045096.View ArticlePubMedGoogle Scholar
- Murdoch S, Djuric U, Mazhar B, Seoud M, Khan R, Kuick R, Bagga R, Kircheisen R, Ao A, Ratti B, et al: Mutations in NALP7 cause recurrent hydatidiform moles and reproductive wastage in humans. Nature genetics. 2006, 38 (3): 300-302. 10.1038/ng1740.View ArticlePubMedGoogle Scholar
- Qian J, Deveault C, Bagga R, Xie X, Slim R: Women heterozygous for NALP7/NLRP7 mutations are at risk for reproductive wastage: report of two novel mutations. Human mutation. 2007, 28 (7): 741-10.1002/humu.9498.View ArticlePubMedGoogle Scholar
- McDaniel P, Wu X: Identification of oocyte-selective NLRP genes in rhesus macaque monkeys (Macaca mulatta). Molecular reproduction and development. 2009, 76 (2): 151-159. 10.1002/mrd.20937.View ArticlePubMedGoogle Scholar
- Zhang P, Dixon M, Zucchelli M, Hambiliki F, Levkov L, Hovatta O, Kere J: Expression analysis of the NLRP gene family suggests a role in human preimplantation development. PLoS ONE. 2008, 3 (7): e2755-10.1371/journal.pone.0002755.PubMed CentralView ArticlePubMedGoogle Scholar
- Dalbies-Tran R, Papillier P, Pennetier S, Uzbekova S, Monget P: Bovine mater-like NALP9 is an oocyte marker gene. Molecular reproduction and development. 2005, 71 (4): 414-421. 10.1002/mrd.20298.View ArticlePubMedGoogle Scholar
- Tong ZB, Bondy CA, Zhou J, Nelson LM: A human homologue of mouse Mater, a maternal effect gene essential for early embryonic development. Human reproduction (Oxford, England). 2002, 17 (4): 903-911. 10.1093/humrep/17.4.903.View ArticleGoogle Scholar
- Evsikov AV, Graber JH, Brockman JM, Hampl A, Holbrook AE, Singh P, Eppig JJ, Solter D, Knowles BB: Cracking the egg: molecular dynamics and evolutionary aspects of the transition from the fully grown oocyte to embryo. Genes & development. 2006, 20 (19): 2713-2727. 10.1101/gad.1471006.View ArticleGoogle Scholar
- Grenier JM, Wang L, Manji GA, Huang WJ, Al-Garawi A, Kelly R, Carlson A, Merriam S, Lora JM, Briskin M, et al: Functional screening of five PYPAF family members identifies PYPAF5 as a novel regulator of NF-kappaB and caspase-1. FEBS letters. 2002, 530 (1-3): 73-78. 10.1016/S0014-5793(02)03416-6.View ArticlePubMedGoogle Scholar
- Manji GA, Wang L, Geddes BJ, Brown M, Merriam S, Al-Garawi A, Mak S, Lora JM, Briskin M, Jurman M, et al: PYPAF1, a PYRIN-containing Apaf1-like protein that assembles with ASC and regulates activation of NF-kappa B. The Journal of biological chemistry. 2002, 277 (13): 11570-11575. 10.1074/jbc.M112208200.View ArticlePubMedGoogle Scholar
- Martinon F, Burns K, Tschopp J: The inflammasome: a molecular platform triggering activation of inflammatory caspases and processing of proIL-beta. Molecular cell. 2002, 10 (2): 417-426. 10.1016/S1097-2765(02)00599-3.View ArticlePubMedGoogle Scholar
- Wang L, Manji GA, Grenier JM, Al-Garawi A, Merriam S, Lora JM, Geddes BJ, Briskin M, DiStefano PS, Bertin J: PYPAF7, a novel PYRIN-containing Apaf1-like protein that regulates activation of NF-kappa B and caspase-1-dependent cytokine processing. The Journal of biological chemistry. 2002, 277 (33): 29874-29880. 10.1074/jbc.M203915200.View ArticlePubMedGoogle Scholar
- Wang Y, Hasegawa M, Imamura R, Kinoshita T, Kondo C, Konaka K, Suda T: PYNOD, a novel Apaf-1/CED4-like protein is an inhibitor of ASC and caspase-1. International immunology. 2004, 16 (6): 777-786. 10.1093/intimm/dxh081.View ArticlePubMedGoogle Scholar
- Bruey JM, Bruey-Sedano N, Newman R, Chandler S, Stehlik C, Reed JC: PAN1/NALP2/PYPAF2, an inducible inflammatory mediator that regulates NF-kappaB and caspase-1 activation in macrophages. The Journal of biological chemistry. 2004, 279 (50): 51897-51907. 10.1074/jbc.M406741200.View ArticlePubMedGoogle Scholar
- Fiorentino L, Stehlik C, Oliveira V, Ariza ME, Godzik A, Reed JC: A novel PAAD-containing protein that modulates NF-kappa B induction by cytokines tumor necrosis factor-alpha and interleukin-1beta. The Journal of biological chemistry. 2002, 277 (38): 35333-35340. 10.1074/jbc.M200446200.View ArticlePubMedGoogle Scholar
- Fontalba A, Gutierrez O, Fernandez-Luna JL: NLRP2, an inhibitor of the NF-kappaB pathway, is transcriptionally activated by NF-kappaB and exhibits a nonfunctional allelic variant. J Immunol. 2007, 179 (12): 8519-8524.View ArticlePubMedGoogle Scholar
- Kinoshita T, Wang Y, Hasegawa M, Imamura R, Suda T: PYPAF3, a PYRIN-containing APAF-1-like protein, is a feedback regulator of caspase-1-dependent interleukin-1beta secretion. The Journal of biological chemistry. 2005, 280 (23): 21720-21725. 10.1074/jbc.M410057200.View ArticlePubMedGoogle Scholar
- Slim R, Mehio A: The genetics of hydatidiform moles: new lights on an ancient disease. Clinical genetics. 2007, 71 (1): 25-34. 10.1111/j.1399-0004.2006.00697.x.View ArticlePubMedGoogle Scholar
- Gu X, Velden Vander K: DIVERGE: phylogeny-based analysis for functional-structural divergence of a protein family. Bioinformatics (Oxford, England). 2002, 18 (3): 500-501. 10.1093/bioinformatics/18.3.500.View ArticleGoogle Scholar
- Halpern M: The organization and function of the vomeronasal system. Annual review of neuroscience. 1987, 10: 325-362. 10.1146/annurev.ne.10.030187.001545.View ArticlePubMedGoogle Scholar
- Keverne EB: The vomeronasal organ. Science (New York, NY). 1999, 286 (5440): 716-720.View ArticleGoogle Scholar
- Shi P, Bielawski JP, Yang H, Zhang YP: Adaptive diversification of vomeronasal receptor 1 genes in rodents. Journal of molecular evolution. 2005, 60 (5): 566-576. 10.1007/s00239-004-0172-y.View ArticlePubMedGoogle Scholar
- Samonte RV, Eichler EE: Segmental duplications and the evolution of the primate genome. Nature reviews. 2002, 3 (1): 65-72.View ArticlePubMedGoogle Scholar
- Darbo E, Danchin EG, Mc Dermott MF, Pontarotti P: Evolution of major histocompatibility complex by "en bloc" duplication before mammalian radiation. Immunogenetics. 2008, 60 (8): 423-438. 10.1007/s00251-008-0301-7.View ArticlePubMedGoogle Scholar
- Bailey JA, Eichler EE: Primate segmental duplications: crucibles of evolution, diversity and disease. Nature reviews. 2006, 7 (7): 552-564. 10.1038/nrg1895.View ArticlePubMedGoogle Scholar
- Johnson ME, Viggiano L, Bailey JA, Abdul-Rauf M, Goodwin G, Rocchi M, Eichler EE: Positive selection of a gene family during the emergence of humans and African apes. Nature. 2001, 413 (6855): 514-519. 10.1038/35097067.View ArticlePubMedGoogle Scholar
- Yang Z: PAML 4: phylogenetic analysis by maximum likelihood. Molecular biology and evolution. 2007, 24 (8): 1586-1591. 10.1093/molbev/msm088.View ArticlePubMedGoogle Scholar
- Duret L, Mouchiroud D: Determinants of substitution rates in mammalian genes: expression pattern affects selection intensity but not mutation rate. Molecular biology and evolution. 2000, 17 (1): 68-74.View ArticlePubMedGoogle Scholar
- Winter EE, Goodstadt L, Ponting CP: Elevated rates of protein secretion, evolution, and disease among tissue-specific genes. Genome research. 2004, 14 (1): 54-61. 10.1101/gr.1924004.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang L, Li WH: Mammalian housekeeping genes evolve more slowly than tissue-specific genes. Molecular biology and evolution. 2004, 21 (2): 236-239. 10.1093/molbev/msh010.View ArticlePubMedGoogle Scholar
- Singh RS: Patterns of species divergence and genetic theories of speciation. Topics in Population Biology and Evolution. Edited by: Wohrmann K, Jain S. 1990, New York: Springer-Verlag, 231-265.View ArticleGoogle Scholar
- Swanson WJ, Vacquier VD: The rapid evolution of reproductive proteins. Nature reviews. 2002, 3 (2): 137-144.View ArticlePubMedGoogle Scholar
- Torgerson DG, Kulathinal RJ, Singh RS: Mammalian sperm proteins are rapidly evolving: evidence of positive selection in functionally diverse genes. Molecular biology and evolution. 2002, 19 (11): 1973-1980.View ArticlePubMedGoogle Scholar
- Torgerson DG, Singh RS: Rapid evolution through gene duplication and subfunctionalization of the testes-specific alpha4 proteasome subunits in Drosophila. Genetics. 2004, 168 (3): 1421-1432. 10.1534/genetics.104.027631.PubMed CentralView ArticlePubMedGoogle Scholar
- Westerveld GH, Korver CM, van Pelt AM, Leschot NJ, Veen van der F, Repping S, Lombardi MP: Mutations in the testis-specific NALP14 gene in men suffering from spermatogenic failure. Human reproduction (Oxford, England). 2006, 21 (12): 3178-3184. 10.1093/humrep/del293.View ArticleGoogle Scholar
- Huminiecki L, Wolfe KH: Divergence of spatial gene expression profiles following species-specific gene duplications in human and mouse. Genome research. 2004, 14 (10A): 1870-1879. 10.1101/gr.2705204.PubMed CentralView ArticlePubMedGoogle Scholar
- Gu Z, Rifkin SA, White KP, Li WH: Duplicate genes increase gene expression diversity within and between species. Nature genetics. 2004, 36 (6): 577-579. 10.1038/ng1355.View ArticlePubMedGoogle Scholar
- Freilich S, Massingham T, Blanc E, Goldovsky L, Thornton JM: Relating tissue specialization to the differentiation of expression of singleton and duplicate mouse proteins. Genome biology. 2006, 7 (10): R89-10.1186/gb-2006-7-10-r89.PubMed CentralView ArticlePubMedGoogle Scholar
- Force A, Lynch M, Pickett FB, Amores A, Yan YL, Postlethwait J: Preservation of duplicate genes by complementary, degenerative mutations. Genetics. 1999, 151 (4): 1531-1545.PubMed CentralPubMedGoogle Scholar
- She X, Cheng Z, Zollner S, Church DM, Eichler EE: Mouse segmental duplication and copy number variation. Nature genetics. 2008, 40 (7): 909-914. 10.1038/ng.172.PubMed CentralView ArticlePubMedGoogle Scholar
- Emanuel BS, Shaikh TH: Segmental duplications: an 'expanding' role in genomic instability and disease. Nature reviews. 2001, 2 (10): 791-800. 10.1038/35093500.View ArticlePubMedGoogle Scholar
- Giglio S, Broman KW, Matsumoto N, Calvari V, Gimelli G, Neumann T, Ohashi H, Voullaire L, Larizza D, Giorda R, et al: Olfactory receptor-gene clusters, genomic-inversion polymorphisms, and common chromosome rearrangements. American journal of human genetics. 2001, 68 (4): 874-883. 10.1086/319506.PubMed CentralView ArticlePubMedGoogle Scholar
- Saglio G, Storlazzi CT, Giugliano E, Surace C, Anelli L, Rege-Cambrin G, Zagaria A, Jimenez Velasco A, Heiniger A, Scaravaglio P, et al: A 76-kb duplicon maps close to the BCR gene on chromosome 22 and the ABL gene on chromosome 9: possible involvement in the genesis of the Philadelphia chromosome translocation. Proceedings of the National Academy of Sciences of the United States of America. 2002, 99 (15): 9882-9887. 10.1073/pnas.152171299.PubMed CentralView ArticlePubMedGoogle Scholar
- Shaw CJ, Bi W, Lupski JR: Genetic proof of unequal meiotic crossovers in reciprocal deletion and duplication of 17p11.2. American journal of human genetics. 2002, 71 (5): 1072-1081. 10.1086/344346.PubMed CentralView ArticlePubMedGoogle Scholar
- Armengol L, Marques-Bonet T, Cheung J, Khaja R, Gonzalez JR, Scherer SW, Navarro A, Estivill X: Murine segmental duplications are hot spots for chromosome and gene evolution. Genomics. 2005, 86 (6): 692-700. 10.1016/j.ygeno.2005.08.008.View ArticlePubMedGoogle Scholar
- Jiang Z, Tang H, Ventura M, Cardone MF, Marques-Bonet T, She X, Pevzner PA, Eichler EE: Ancestral reconstruction of segmental duplications reveals punctuated cores of human genome evolution. Nature genetics. 2007, 39 (11): 1361-1368. 10.1038/ng.2007.9.View ArticlePubMedGoogle Scholar
- Clark NL, Aagaard JE, Swanson WJ: Evolution of reproductive proteins from animals and plants. Reproduction (Cambridge, England). 2006, 131 (1): 11-22.View ArticleGoogle Scholar
- Osawa M, Horiuchi H, Tian W, Kaneko M: Divergent evolution of the prolactin-inducible protein gene and related genes in the mouse genome. Gene. 2004, 325: 179-186. 10.1016/j.gene.2003.10.013.View ArticlePubMedGoogle Scholar
- Birtle Z, Goodstadt L, Ponting C: Duplication and positive selection among hominin-specific PRAME genes. BMC genomics. 2005, 6: 120-10.1186/1471-2164-6-120.PubMed CentralView ArticlePubMedGoogle Scholar
- McLellan AS, Zimmermann W, Moore T: Conservation of pregnancy-specific glycoprotein (PSG) N domains following independent expansions of the gene families in rodents and primates. BMC evolutionary biology. 2005, 5: 39-10.1186/1471-2148-5-39.PubMed CentralView ArticlePubMedGoogle Scholar
- Jackson M, Watt AJ, Gautier P, Gilchrist D, Driehaus J, Graham GJ, Keebler J, Prugnolle F, Awadalla P, Forrester LM: A murine specific expansion of the Rhox cluster involved in embryonic stem cell biology is under natural selection. BMC genomics. 2006, 7: 212-10.1186/1471-2164-7-212.PubMed CentralView ArticlePubMedGoogle Scholar
- Karn RC, Clark NL, Nguyen ED, Swanson WJ: Adaptive evolution in rodent seminal vesicle secretion proteins. Molecular biology and evolution. 2008, 25 (11): 2301-2310. 10.1093/molbev/msn182.PubMed CentralView ArticlePubMedGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic acids research. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralView ArticlePubMedGoogle 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 research. 1994, 22 (22): 4673-4680. 10.1093/nar/22.22.4673.PubMed CentralView ArticlePubMedGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Molecular biology and evolution. 2007, 24 (8): 1596-1599. 10.1093/molbev/msm092.View ArticlePubMedGoogle Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Systematic biology. 2003, 52 (5): 696-704. 10.1080/10635150390235520.View ArticlePubMedGoogle Scholar
- Felsenstein J: PHYLIP (Phylogeny Inference Package) version 3.6. 2004, Distributed by the author Department of Genome Sciences, University of Washington, SeattleGoogle Scholar
- Felsenstein J: Confidence limits on phylogenies: An approach using the bootstrap. Evolution. 1985, 39: 783-791. 10.2307/2408678.View ArticleGoogle Scholar
- Gu X: Maximum-likelihood approach for gene family evolution under functional divergence. Molecular biology and evolution. 2001, 18 (4): 453-464.View ArticlePubMedGoogle Scholar
- Gu X: Statistical methods for testing functional divergence after gene duplication. Molecular biology and evolution. 1999, 16 (12): 1664-1674.View ArticlePubMedGoogle Scholar
- Suyama M, Torrents D, Bork P: PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignments. Nucleic acids research. 2006, W609-612. 10.1093/nar/gkl315. 34 Web Server
- Nielsen R, Yang Z: Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene. Genetics. 1998, 148 (3): 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 (1): 431-449.PubMed CentralPubMedGoogle Scholar
- Swanson WJ, Nielsen R, Yang Q: Pervasive adaptive evolution in mammalian fertilization proteins. Molecular biology and evolution. 2003, 20 (1): 18-20.View ArticlePubMedGoogle Scholar
- Anisimova M, Nielsen R, Yang Z: Effect of recombination on the accuracy of the likelihood method for detecting positive selection at amino acid sites. Genetics. 2003, 164 (3): 1229-1236.PubMed CentralPubMedGoogle Scholar
- Yang Z, Wong WS, Nielsen R: Bayes empirical bayes inference of amino acid sites under positive selection. Molecular biology and evolution. 2005, 22 (4): 1107-1118. 10.1093/molbev/msi097.View ArticlePubMedGoogle Scholar
- Doron-Faigenboim A, Pupko T: A combined empirical and mechanistic codon model. Molecular biology and evolution. 2007, 24 (2): 388-397. 10.1093/molbev/msl175.View ArticlePubMedGoogle Scholar
- Stern A, Doron-Faigenboim A, Erez E, Martz E, Bacharach E, Pupko T: Selecton 2007: advanced models for detecting positive and purifying selection using a Bayesian inference approach. Nucleic acids research. 2007, W506-511. 10.1093/nar/gkm382. 35 Web Server
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.