Evolutionary origin of peptidoglycan recognition proteins in vertebrate innate immune system
© Montaño et al; licensee BioMed Central Ltd. 2011
Received: 16 June 2010
Accepted: 25 March 2011
Published: 25 March 2011
Innate immunity is the ancient defense system of multicellular organisms against microbial infection. The basis of this first line of defense resides in the recognition of unique motifs conserved in microorganisms, and absent in the host. Peptidoglycans, structural components of bacterial cell walls, are recognized by Peptidoglycan Recognition Proteins (PGRPs). PGRPs are present in both vertebrates and invertebrates. Although some evidence for similarities and differences in function and structure between them has been found, their evolutionary history and phylogenetic relationship have remained unclear. Such studies have been severely hampered by the great extent of sequence divergence among vertebrate and invertebrate PGRPs. Here we investigate the birth and death processes of PGRPs to elucidate their origin and diversity.
We found that (i) four rounds of gene duplication and a single domain duplication have generated the major variety of present vertebrate PGRPs, while in invertebrates more than ten times the number of duplications are required to explain the repertoire of present PGRPs, and (ii) the death of genes in vertebrates appears to be almost null whereas in invertebrates it is frequent.
These results suggest that the emergence of new PGRP genes may have an impact on the availability of the repertoire and its function against pathogens. These striking differences in PGRP evolution of vertebrates and invertebrates should reflect the differences in the role of their innate immunity. Insights on the origin of PGRP genes will pave the way to understand the evolution of the interaction between host and pathogens and to lead to the development of new treatments for immune diseases that involve proteins related to the recognition of self and non-self.
Innate immunity is the ancient defense system of multicellular organisms against microbial infection. The basis of this first line of defense resides in the recognition of unique motifs or components conserved in microorganisms, and absent in the host. The innate immune system uses sets of pattern recognition receptors to recognize such foreign or non-self motifs. Proteins in the immune system can be located intracellularly, on the cell surface, or secreted into the bloodstream, ready to signal the presence of an intruder in every compartment. In systems lacking the adaptive arm of immunity, the pattern recognition concept serves well to explain the general triggering of the system as well as providing receptors for the limited specificity shown by innate immunity [1–4].
Peptidoglycan (PGN) is the major structural component of the cell wall of almost all bacterial species. PGN is a large, repetitive macromolecule that forms the rigid cell wall of bacteria. PGN recognition is mediated by the PGRP (PGN recognition protein) family of receptors [5, 6]. PGRPs are a family of innate immunity pattern recognition molecules that were first discovered in silkworms . There are four loci for each PGRPs in humans [8, 9] (PGRP-S, PGRP-L, PGRP-Iα and PGRP-Iβ), while thirteen loci in Drosophila [10, 11], which encode approximately 17 PGRP proteins through alternative splicing, and seven loci in Anopheles . Several other genomes also show relatively large number of PGRPs in invertebrates, but only up to five in vertebrates (Additional file 1). In invertebrates, the functional divergence of each PGRP molecule is well investigated: Some possess an amidase activity that hydrolyzes the amide bond between the N-acetylmuramic acid and the L-alanine of peptidoglycan, others activate Toll, or Imd pathways to induce an expression of anti-bacteria peptides, induce prophenoloxidase cascade, or directly cause phagocytosis and lysis [13–18]. On the other hand, functions of vertebrate, or mainly mammalian PGRPs, are not fully understood . While PGRP-L has the amidase activity where its role could be to detoxify PGN fragments present in blood and modulate the immune response as insect PGRPs, the PGRP-S, PGRP-Iα and PGRP-Iβ have bacteriostatic and/or bactericidal function [6, 20–22].
Vertebrates have the acquired immunity system in addition to the innate immunity system, while in insects only the latter is a self-defense system. It is of interest whether possessing acquired immunity has effects on the evolution of molecules involved in innate immunity. Here, we investigate the birth and death processes of PGRPs by systematically analyzing PGRP genes from a set of diverse eukaryotes, and discuss the role of selection and diversification of this gene family.
Modes of PGRP evolution
To detect lineage-specific expansions 40 sequences of the PGRP family from 21 vertebrate species and 42 sequences from 6 invertebrate species were studied (Additional file 1 and Additional file 2). Both vertebrate and invertebrate PGRPs have a highly conserved C-terminal region of the PGRP domain with three sub-domains (I, II and III). The sub-domains are determined by sequence conservation and not by their function. The PGRP domain shows a sequence similarity (~35%) with bacteriophage T7 lysozyme, which also has amidase activity, indicating that T7 lysozyme would be the origin of PGRP domains . This ancient origin of PGRP domains is also supported by the similarity of the 3D structure between PGRP-L and T7 lysozyme molecules. However, the orthologous relationship of vertebrate and insect PGRP domains has not been ensured [8, 19] due to the limited number of amino acid sites compared and the great extent of sequence divergence among them. Contrary to the conserved C-terminal PGRP domains, the N-terminal region shows no particular similarities among different PGRPs in invertebrates, and partial similarities among PGRP-S, PGRP-Iα and PGRP-Iβ in vertebrates. Therefore, we used only the PGRP domains for the alignment and tree construction. Due to difficulties in identifying orthologous relationship, we performed phylogenetic analyses of PGRP genes in vertebrates and invertebrates separately by using the neighbor-joining (NJ) and minimum evolution (ME) methods. There was no conflict in the topologies obtained with these methods.
The mentioned observation could be confirmed by using other methods that can predict the number of gene gains and losses. We verified our observations by using the program NOTUNG and EvolMAP [25, 26]. The results with NOTUNG showed a similar tendency of the number of gene gains predicted (Additional file 5 and Additional file 6), however the number of gene losses seems to be overestimated in vertebrates, especially in fish. In NOTUNG, absence of a gene in a particular taxon means gene loss. Thus the number of losses in fish became enormously large. EvolMAP, on the other hand, predicts that gene gains are 8.6 times more frequent than gene losses in invertebrates over all branches. This suggests an expansion of PGRPs in invertebrates. For vertebrates, EvolMAP analysis shows no evidence of expansion or contraction for this gene family (Additional file 7 and Additional file 8). Overall, for the genes and species analyzed here, we find that the number of gains detected in invertebrates is twice the number of gains in vertebrates. Thus we could confirm the large number of gene gains and losses in invertebrates when compared to vertebrates.
Ancestral PGRP genes
Due to the limited number of sites compared and long divergence time of sequences, we could not elucidate the relationship among the ancestors of vertebrate and invertebrate PGRPs from the phylogenetic tree of vertebrate and invertebrate PGRPs, including T7 lysozyme. Thus whether the origin of vertebrate PGRPs is monophyletic or paraphyletic to invertebrate PGRPs remains to be an open problem. However, our analysis clearly shows that for vertebrate PGRPs, the first major divergence took place between PGRP-L and PGRP-S. Therefore, in the following we focus on the vertebrate PGRPs to infer the functions of the ancestral PGRP genes.
To elucidate the function of ancestral PGRP molecules in vertebrates, the amino acid sequences of proto- PGRP-L and PGRP-S molecules were estimated by the maximum-likelihood (ML) method with the JTT substitution matrix . It is known that seven amino acids are responsible for PGRP function [1, 8]. Four amino acid residues (H17, Y46, H122, and C130) are essential for the amidase activity, whereas three (H36, W41, and K128) are important for Zn2+ ligand-binding in the bacteriophage T7 lysozyme. Since all seven amino acids are conserved in both the proto- PGRP-L and PGRP-S sequences, the ancestor of both proto-PGRPs is likely to possess the amidase activity (Additional file 9). While the present-day PGRP-L has reserved its original function of amidase activity , PGRP-S has lost it and instead obtained the bacteriostatic function [21, 22]. On the other hand, the invertebrate PGRPs possessing the amidase activity are paraphyletic to each other. This suggests independent gain or loss of the amidase activity in invertebrate PGRP at an early stage of the evolution.
Selection acting on PGRPs and Evolutionary Rates
Next important question is whether some kind of selection process has acted on each amino acid site of vertebrate PGRP genes that will lead to their functional divergence after gene duplication. We identified positively or negatively selected sites in vertebrate PGRPs using Single Likelihood Ancestor Counting (SLAC) analysis as described in Methods . Only the site 138 in PGRP-L is positively selected among all the vertebrate PGRPs. This site, which is involved in substrate binding, shows a high degree of amino acid variation in PGRP-L of different species.
Average values of the ratios of non-synonymous to synonymous substitutions of PGRP-S, L and I are 0.16, 0.13 and 0.17, respectively, and the overall value for PGRPs is 0.20. This indicates strong functional constraint, suggesting that amino acid sequences of these domains are well conserved among vertebrates.
Although there is a report where a few positively selected sites were observed in the PGRP domain of Drosophila PGRP-LC , we could not apply the above SLAC analysis to invertebrate PGRPs, because the clustering pattern and phylogenetic relationship among sequences are not ensured and the results of SLAC strongly depend on the tree topology (data not shown).
Tests of parallel and convergent evolution of PGRPs.
A. Parallel changes
Pig-S & Cow-S
2, 78, 113
Frog-S2 & Rat-S
Rainbow trout-L1 & stickleback-L
45, 52, 120
Rat-Ia-C & Rat-S
Fugu-L & Pig-L
Chicken-L & Salmon-L
Rainbow trout-L1 & Zebrafish-L1
45, 100, 114
Stickleback-L & Rat-Ia-C
B. Convergent changes
Human-Ib-N & Rat-Ib-C
Cow-S & Frog-L
Pig-S & Frog-S1
Fugu-L & Hagfish-S
Frog-S2 & Rat-S
Mouse-L & Killifish-L
Mouse-L & Rat-Ia-N
Rat-Ia-N & Killifish-L
Carp-L & Hagfish-S
Mouse-Ia-N & Pig-S
Mouse-S & Mouse-Ib-C
Rat-Ib-N & Pig-S
Rainbow trout-L1 & Mouse-Ia-N
Fugu-L & Pig-L
Pig-L & Human-Ia-C
Rat-S & Rainbow trout-L1
Mouse-L & Zebrafish-L1
Rat-S & Cow-S
Human-S & Salmon-L
Pig-S & Salmon-L
Human-S & Chicken-L
Frog-S1 & Mouse-Ia-C
Mouse-L & Mouse-Ia-N
Human-Ia-C & Rat-S
Carp-L & Fugu-L
Zebrafish-L2-N & Medaka-L1
Among the six parallel sites and ten convergent sites, which occurred with a significance level of 1%, sites 45, 78, and 123 may be potential adaptive sites. Site 45 has three independent parallel changes of S (Ser) to A (Ala) only in fish PGRP-L (three spine stickleback, rainbow trout and zebrafish). Serine is the only amino acid residue present in this site in PGRP-L except for these fishes. Moreover, exclusively non-polar character of this site indicates that it may have an important role in the function of these proteins especially because this residue is involved in substrate binding. Site 78 has suffered one parallel and one convergent substitution. The chemical profile of this site in mammal PGRP-S is exclusively hydrophobic. The variability which results from the change of I (Ile) to the polar T (Thr) in pig and cow may suggest an effect on the structure since this amino acid is located in the hydrophobic core. At the site 123 there is a convergent change of A (Ala) or Y (Tyr) to V (Val) in pig and frog PGRP-S, respectively located in the hydrophobic core. Both substitutions create a possible adaptive site in PGRP-S gene.
Evolutionary characteristics of PGRPs in vertebrates and invertebrates
When we compared PGRP evolution in vertebrates and invertebrates, we observed several differences, which are characteristic of each mode of evolution. First of all, despite the similar divergence time of 450 myr for the most recent common ancestor in each vertebrate and invertebrate (insect) species, the phylogenetic trees clearly show different patterns. In invertebrates each PGRP shows a relatively longer branch than those in vertebrates, suggesting a relatively ancient origin of each PGRP in invertebrates. In addition, in invertebrates the clustering pattern was rarely orthologous among PGRP genes, while in vertebrates orthologous relationship was clearly seen. This observation shows that higher rates of birth and death processes are seen in invertebrates than in vertebrates. Although the repertoire of PGRPs in each species may depend on some ecological and biological conditions, a less frequent birth and death process in vertebrates could reflect the presence of acquired immune system.
Consequences of Natural Selection on PGRPs
In vertebrate PGRP proteins we observed changes as consequence of parallel and convergent amino acid substitutions, with significance greater than the random chance expectation. The changes may be either due to conservation of chemical property of amino acids at the site or due to modification of their properties. Convergent or parallel substitutions can provide evidence for the action of natural selection for keeping the function or structure .
PGRP proteins play an important role in innate immunity, which requires updated and immediate responses, because pathogens may change frequently. As a consequence, a high turnover rate is expected to happen. Actually, in invertebrates, the frequent turnover of PGRP repertoire was observed. On the other hand, several motifs indispensable for peptidoglycan recognition should be conserved through evolution of PGRPs in both vertebrates and invertebrates. In addition, we have observed that the amino acid residues that are located on the hydrophobic groove have high degree of conservation and do not show any parallel or convergent amino acid substitution.
Evolution of the PGRP family in vertebrates and invertebrates and functional implications
This study provides for the first time a description of the origin and mode of evolution on vertebrate PGRPs, compared with invertebrates, namely insects, PGRPs.
PGRPs are proposed to be a family of genes that evolved by birth and death process with different rates in vertebrates and invertebrates. In the model of the birth and death process [33, 34], some of the duplicated genes diverge functionally, but others become pseudogenes due to deleterious mutations or are deleted from the genome. The end result of this mode of evolution is a multi-gene family with a mixture of divergent groups of genes and highly homologous genes. We have observed that PGRPs have experienced several rounds of gene duplications and some duplicated genes have been deleted from the genome. This lineage specific birth and death process has been observed both in vertebrates and invertebrates.
The PGRP proteins are involved in innate immunity, which responds to protect the organisms from invading pathogens. Therefore, several motifs in PGRP domains are, of course, indispensable for pathogen recognition and have been conserved through the vertebrate and invertebrate evolution. However, since vertebrates possess acquired immunity, the significance of innate immunity might be more relaxed than in insects whose immune systems depend solely on innate immunity. This difference in the evolutionary patterns could be related to the plasticity of the receptors to detect a broad spectrum of microbial pathogens and it is clearly reflected in the birth and death process of the PGRP molecules in vertebrates and invertebrates.
PGRP gene family reveals an example of genetic and functional variation of which roles in the immune systems are understood through an analysis of comparative genomics. Especially the analysis reveals that the mode of PGRP evolution was characterized by birth and death process. Vertebrates and invertebrates show striking differences in the evolutionary tempo and mode of PGRP genes. Broad repertoire of pathogen recognition proteins is advantageous in invertebrates, due to the absence of adaptive immunity, in contrast to the moderate repertoire in vertebrates. This reveals that the mode of evolution of a system strongly depends on other systems, which interact with the former both directly or indirectly.
The sequences were retrieved from the genomic and EST NCBI database http://www.ncbi.nlm.nih.gov, the Ensembl database http://www.ensembl.org and the TIGR database http://www.tigr.org/tdb/tgi/ using TBLASTN and PSI-BLAST http://www.ncbi.nlm.nih.gov/BLAST/. Search was performed using each exon of human PGRP-S, PGRP-L, PGRP-Iα, and PGRP-Iβ as a probe. Takifugu rubripes PGRP-L full cDNA and genomic sequences were predicted using Genscan and confirmed by sequence analysis. Takifugu rubripes liver tissue was kindly provided by Dr. Shugo Watabe of the University of Tokyo, Japan. The Eptatretus burgeri PGRP cDNA sequence was kindly provided by Dr. Masanori Kasahara of The Graduate University for Advanced Studies (Sokendai), Hayama, Japan.
GenBank accession numbers
Rapid amplification of cDNA ends (RACE)
Based on the partial sequence information of chicken PGRP-L retrieved from databases, we reconstructed the missing 3' region of the transcript using the BD SMART™ RACE cDNA Amplification Kit (Clontech, USA) according to manufacturer's instructions.
Sequences were aligned using the CLUSTALW version 1.83 computer program with its default parameter setting  and manually adjusted using the GeneDoc program version 2.6.002 . Phylogenetical analyses were done using the neighbor-joining (NJ) and minimum evolution (ME) methods in MEGA version 4 . The NJ and ME trees were based on the number of differences, and reliability was assessed by bootstrap values with 1000 replications. The reconciliation between species tree and gene tree along with the confirmation of the gene loss/duplication scenario were determined by using Notung 2.6  and EvolMAP software . To detect positive selection at single amino acid sites the Data Monkey software program was used with its default parameter setting http://www.datamonkey.org/. Poisson and gamma genetic distances were determined by using MEGA version 4 .
Test of convergence
- The following abbreviations used in the manuscript are listed here in alphabetical order:
ME: (minimum evolution)
(Peptidoglycan Recognition Proteins)
(Rapid amplification of cDNA ends)
(Single Likelihood Ancestor Counting).
AMM is grateful to the Japan Society for the Promotion of Science (JSPS) for support. This work was supported in part by a Grant-in-Aid for Scientific Research on Priority Areas (13143202) from The Ministry of Education, Culture, Sports, Science, and Technology of Japan. We thank V. Byrappa for helpful discussions on Fugu PGRPs and Dr O. Sakarya for assistance with EvolMAP software. We thank Dr. A. Noguchi and Dr. Ş. K. Özdemir for critical reviewing of the manuscript. This paper is dedicated to Dr. Fumi Tsujino, who was deceased before completing the manuscript, on 27 October, 2006.
- Steiner H: Peptidoglycan recognition proteins: on and off switches for innate immunity. Immunol Rev. 2004, 198: 83-96. 10.1111/j.0105-2896.2004.0120.x.View ArticlePubMed
- Dziarski R, Gupta D: The peptidoglycan recognition proteins (PGRPs). Genome Biol. 2006, 7: 232-10.1186/gb-2006-7-8-232.View ArticlePubMedPubMed Central
- Girardin SE, Philpott DJ: The role of peptidoglycan recognition in innate immunity. Eur J Immunol. 2004, 34: 1777-1782. 10.1002/eji.200425095.View ArticlePubMed
- Janeway CA, Medzhitov R: Innate immune recognition. Annu Rev Immunol. 2002, 20: 197-216. 10.1146/annurev.immunol.20.083001.084359.View ArticlePubMed
- Aggrawal K, Silverman N: Peptidoglycan recognition in Drosophila. Biochem Soc Trans. 2007, 35: 1496-1500. 10.1042/BST0351496.View ArticlePubMed
- Chaput C, Boneca IG: Peptidoglycan detection by mammals and flies. Microbes Infect. 2007, 9: 637-647. 10.1016/j.micinf.2007.01.022.View ArticlePubMed
- Yoshida H, Kinoshita K, Ashida M: Purification of a peptidoglycan recognition protein from hemolymph of the silkworm, Bombyx mori. J Biol Chem. 1996, 271: 13854-13860. 10.1074/jbc.271.23.13854.View ArticlePubMed
- Kang D, Liu G, Lundström A, Gelius E, Steiner H: A peptidoglycan recognition protein in innate immunity conserved from insects to humans. Proc Natl Acad Sci. 1998, 95: 10078-10082. 10.1073/pnas.95.17.10078.View ArticlePubMedPubMed Central
- Liu C, Xu Z, Gupta D, Dziarski R: Peptidoglycan recognition proteins: a novel family of four human innate immunity pattern recognition molecules. J Biol Chem. 2003, 276: 34686-34694. 10.1074/jbc.M105566200.View Article
- Werner T, Liu G, Kang D, Ekengren S, Steiner H, Hultmark D: A family of peptidoglycan recognition proteins in the fruit fly Drosophila melanogaster. Proc Natl Acad Sci. 2000, 97: 13772-13777. 10.1073/pnas.97.25.13772.View ArticlePubMedPubMed Central
- Werner T, Borge-Renberg K, Mellroth P, Steiner H, Hultmark D: Functional diversity of the Drosophila PGRP-LC gene cluster in the response to lipopolysaccharide and peptidoglycan. J Biol Chem. 2003, 278: 26319-26322. 10.1074/jbc.C300184200.View ArticlePubMed
- Christophides GK, Zdobnov E, Barillas-Mury C, Birney E, Blandin S, Blass C, Brey PT, Collins FH, Danielli A, Dimopoulos G, et al: Immunity-related genes and gene families in Anopheles gambiae. Science. 2002, 298: 159-165. 10.1126/science.1077136.View ArticlePubMed
- Takehana A, Yano T, Mita S, Kotani A, Oshima Y, Kurata S: Peptidoglycan recognition protein (PGRP)-LE and PGRP-LC act synergistically in Drosophila immunity. EMBO J. 2004, 23: 4690-4700. 10.1038/sj.emboj.7600466.View ArticlePubMedPubMed Central
- Bischoff V, Vignal C, Boneca IG, Michel T, Hoffmann JA, Royet J: Function of the Drosophila pattern-recognition receptor PGRP-SD in the detection of Gram-positive bacteria. Nat Immunol. 2004, 5: 1175-1180. 10.1038/ni1123.View ArticlePubMed
- Choe KM, Werner T, Stoven S, Hultmark D, Anderson KV: Requirement for a peptidoglycan recognition protein (PGRP) in Relish activation and antibacterial immune responses in Drosophila. Science. 2002, 296: 359-362. 10.1126/science.1070216.View ArticlePubMed
- Leulier F, Parquet C, Pili-Floury S, Ryu JH, Caroff M, Lee WJ, Mengin-Lecreulx D, Lemaitre B: The Drosophila immune system detects bacteria through specific peptidoglycan recognition. Nat Immunol. 2003, 4: 478-484. 10.1038/ni922.View ArticlePubMed
- Mellroth P, Karlsson J, Steiner H: A scavenger function for a Drosophila peptidoglycan recognition protein. J Biol Chem. 2003, 278: 7059-7064. 10.1074/jbc.M208900200.View ArticlePubMed
- Chang CI, Pili-Floury S, Hervé M, Parquet C, Chelliah Y, Lemaitre B, Mengin-Lecreulx D, Deisenhofer J: A Drosophila pattern recognition receptor contains a peptidoglycan docking groove and unusual L, D-carboxypeptidase activity. PLoS Biol. 2004, 2: E277-10.1371/journal.pbio.0020277.View ArticlePubMedPubMed Central
- Dziarski R, Gupta D: Mammalian PGRPs: novel antibacterial proteins. Cell Microbiol. 2006, 8: 1059-1069. 10.1111/j.1462-5822.2006.00726.x.View ArticlePubMed
- Gelius E, Persson C, Karlsson J, Steiner H: A mammalian peptidoglycan recognition protein with N-acetylmuramoyl-L-alanine amidase activity. Biochem Biophys Res Commun. 2003, 306: 988-994. 10.1016/S0006-291X(03)01096-9.View ArticlePubMed
- Dziarski R, Platt KA, Gelius E, Steiner H, Gupta D: Defect in neutrophil killing and increased susceptibility to infection with non-pathogenic gram-positive bacteria in peptidoglycan recognition protein-S (PGRP-S)-deficient mice. Blood. 2003, 102: 689-697. 10.1182/blood-2002-12-3853.View ArticlePubMed
- Lu X, Wang M, Qi J, Wang H, Li X, Gupta D, Dziarski R: Peptidoglycan recognition proteins are a new class of human bactericidal proteins. J Biol Chem. 2006, 281: 5895-5907. 10.1074/jbc.M511631200.View ArticlePubMed
- Kim MS, Byun M, Oh BH: Crystal structure of peptidoglycan recognition protein LB from Drosophila melanogaster. Nature Immunology. 2003, 4: 787-793. 10.1038/ni952.View ArticlePubMed
- Guan R, Malchiodi EL, Wang Q, Schuck P, Mariuzza RA: Crystal structure of the C-terminal peptidoglycan-binding domain of human peptidoglycan recognition protein Iα. J Biol Chem. 2004, 279: 31873-31882. 10.1074/jbc.M404920200.View ArticlePubMed
- Durand D, Halldorsson BV, Vernot B: A Hybrid Micro-Macroevolutionary Approach to Gene Tree Reconstruction. J Comput Biol. 2006, 13: 320-335. 10.1089/cmb.2006.13.320.View ArticlePubMed
- Sakarya O, Kosik KS, Oakley TH: Reconstructing ancestral genome content based on symmetrical best alignments and Dollo parsimony. Bioinformatics. 2008, 24: 606-612. 10.1093/bioinformatics/btn005.View ArticlePubMed
- Zhang J, Nei M: Accuracies of ancestral amino acid sequences inferred by the parsimony, likelihood, and distance methods. J Mol Evol. 1997, 44 (Suppl 1): S139-S146. 10.1007/PL00000067.View ArticlePubMed
- 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.View Article
- Sackton TB, Lazzaro BP, Schlenke TA, Evans JD, Hultmark D, Clark AG: Dynamic evolution of the innate immune system in Drosophila. Nat Genet. 2007, 39: 1461-1468. 10.1038/ng.2007.60.View ArticlePubMed
- Reiser JB, Teyton L, Wilson IA: Crystal structure of the Drosophila peptidoglycan recognition protein (PGRP)-SA at 1.56 Å resolution. J Mol Biol. 2004, 340: 909-917. 10.1016/j.jmb.2004.04.077.View ArticlePubMed
- Tajima F: Simple methods for testing the molecular clock hypothesis. Genetics. 1993, 135: 599-607.PubMedPubMed Central
- Zhang J, Kumar S: Detection of convergent and parallel evolution at the amino acid sequence level. Mol Biol Evol. 1997, 14: 527-536.View ArticlePubMed
- Nei M, Rooney AP: Concerted and birth-and-death evolution of multigene families. Annu Rev Genet. 2005, 39: 121-152. 10.1146/annurev.genet.39.073003.112240.View ArticlePubMedPubMed Central
- Nei M, Gu X, Sitnikova T: Evolution by the birth-and-death process in multigene families of the vertebrate immune system. Proc Natl Acad Sci USA. 1997, 94: 7799-7806. 10.1073/pnas.94.15.7799.View ArticlePubMedPubMed Central
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, positions-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.View ArticlePubMedPubMed Central
- Nicholas KB, Nicholas HB, Deerfield DWII: GeneDoc: Analysis and visualization of genetic variation. EMBNEW.NEWS. 1997, 4: 14-
- Kumar S, Tamura K, Jakobsen IB, Nei M: MEGA2: Molecular evolutionary genetics analysis software. Bioinformatics. 2001, 17: 1244-1245. 10.1093/bioinformatics/17.12.1244.View ArticlePubMed
- Blair Hedges S, Kumar S: Genomic clocks and evolutionary timescales. Trends Genet. 2003, 19: 200-206. 10.1016/S0168-9525(03)00053-2.View ArticlePubMed
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.