Perforin evolved from a gene duplication of MPEG1, followed by a complex pattern of gene gain and loss within Euteleostomi
© D'Angelo et al.; licensee BioMed Central Ltd. 2012
Received: 24 November 2011
Accepted: 2 May 2012
Published: 2 May 2012
The pore-forming protein perforin is central to the granule-exocytosis pathway used by cytotoxic lymphocytes to kill abnormal cells. Although this mechanism of killing is conserved in bony vertebrates, cytotoxic cells are present in other chordates and invertebrates, and their cytotoxic mechanism has not been elucidated. In order to understand the evolution of this pathway, here we characterize the origins and evolution of perforin.
We identified orthologs and homologs of human perforin in all but one species analysed from Euteleostomi, and present evidence for an earlier ortholog in Gnathostomata but not in more primitive chordates. In placental mammals perforin is a single copy gene, but there are multiple perforin genes in all lineages predating marsupials, except birds. Our comparisons of these many-to-one homologs of human perforin show that they mainly arose from lineage-specific gene duplications in multiple taxa, suggesting acquisition of new roles or different modes of regulation. We also present evidence that perforin arose from duplication of the ancient MPEG1 gene, and that it shares a common ancestor with the functionally related complement proteins.
The evolution of perforin in vertebrates involved a complex pattern of gene, as well as intron, gain and loss. The primordial perforin gene arose at least 500 million years ago, at around the time that the major histocompatibility complex-T cell receptor antigen recognition system was established. As it is absent from primitive chordates and invertebrates, cytotoxic cells from these lineages must possess a different effector molecule or cytotoxic mechanism.
Cytotoxic lymphocytes (CLs) is a collective term for natural killer (NK) and cytotoxic T lymphocytes (CTL). As the name suggests, these cells are cytotoxic towards virally infected, neoplastic or foreign cells. The two major mechanisms they use to elicit apoptosis in target cells involve (1) cell surface death receptors and their ligands (e.g., Fas/Fas-ligand) and (2) the granule-exocytosis pathway . The latter involves the targeted secretion of specialised secretory lysosomes (granules) from CLs into the immunological synapse, a cleft formed at the site of CL-target cell contact . The granules contain the granzyme family of serine proteases, effectors that cleave cytoplasmic proteins to induce apoptosis, and perforin, a membrane pore forming protein that is required for entry of the granzymes into target cells [3, 4].
CTL form part of the adaptive immune system in all jawed vertebrates but not in earlier chordates [5, 6]. The lamprey, a jawless vertebrate, has an unconventional adaptive immune system which does not use the major histocompatibility complex (MHC) or T cell receptor (TCR) recognition system, but a more primitive leucine rich repeat-containing antigen receptor . This species appears to have CTL-like leukocytes but whether they are armed with granule mediated cytotoxic machinery is unknown. NK cells, by contrast, are more difficult to define than CTL, but appear to have evolved earlier. There is evidence for cells with NK properties in the tunicate Botryllus schlosseri, and the Ciona intestinalis genome includes homologs of some NK cell receptors [5, 8]. More basic cytotoxic NK-like cells have been described in earlier divergent invertebrates such as earthworms . How similar these cells are to conventional mammalian NKs, including their mechanisms of killing, remains to be seen.
Perforin (gene symbol PRF1) is essential and central to the granule-exocytosis pathway in mammals. Effective CL induction of apoptosis requires both granzymes and perforin, although at high concentrations perforin alone can kill cells by causing necrosis, whereas granzymes are ineffective without perforin to translocate them into the target cell cytoplasm. This is highlighted by the human autosomal recessive disease familial hemophagocytic lymphohistiocytosis type 2 (FHL2), caused by mutations in the perforin gene . CTL from these patients cannot kill Fas-deficient target cells and so do not have an active granule-exocytosis pathway .
Perforin forms circular pores in the plasma membrane of target cells by a mechanism involving at least three steps: (1) perforin monomers bind to the membrane via their C2 domains in a calcium dependent manner; (2) monomers polymerise into a ring, mediated in part by salt bridging between residues in adjacent N-terminal membrane attack complex/perforin (MACPF) domains; (3) two clusters of α-helices within each MACPF domain rearrange into anti-parallel β-strands that puncture and span the membrane, creating an aqueous pore [11–16]. The mechanism of (3) and the order of (2) and (3) are inferred from structural similarity to the well-studied cholesterol-dependent cytolysin family of proteins as well as experimental observations of the perforin pore [16, 17].
The MACPF domain has been identified in 12 human proteins and is named after the six best characterised members found in the immune system: five of the terminal complement components (C6, C7, C8α, C8β and C9) that form the membrane attack complex (MAC), and perforin [18, 19]. The MAC is formed when C5b, C6, C7, C8 (a complex of C8α, C8β and C8γ) assemble on foreign cell membranes, which then recruits multiple C9 monomers to polymerise and insert into the membrane [20, 21]. Perforin has long been compared to C9 as they are both able to polymerise and insert into membranes and the pores formed look similar by transmission electron microscopy [22–24].
The only other MACPF domain-containing protein known to be involved in the human immune system is macrophage expressed gene 1 protein (also referred to as mps1 and mpg-1, here-in the gene and protein is abbreviated as MPEG1), produced by macrophages . Besides the MACPF domain MPEG1 contains one or more additional domains with no identified relationship to known protein folds, and a C-terminal transmembrane anchor . MPEG1 is an ancient gene with homologs in species from one of the earliest metazoan lineages, the phylum Porifera (sponges), Amphimedon queenslandica and Suberites domuncula[26, 27]. The homolog from S. domuncula is the best studied MPEG1 gene and is part of an ancient toll-like receptor pathway that is upregulated by lipopolysaccharide . This role in innate immunity, along with its expression in macrophages, has led to the hypothesis that MPEG1 clears phagocytosed Gram-negative bacteria . Indeed, recent evidence shows that the isolated MACPF domain from MPEG1 of the Pacific oyster Crassostrea gigas has anti-microbial activity against both Gram-positive and Gram-negative bacteria .
Here we trace the origins and evolution of the perforin gene to gain insight into the evolution of the granule-exocytosis pathway. Using a variety of approaches including linked gene comparisons, BLAST searches and protein phylogenetic trees we have catalogued all of the available perforin homologs. These data suggest that the perforin-dependent granule-exocytosis pathway originated in jawed vertebrates (Gnathostomata), at around the same time as true CTLs. In addition, we present evidence that MPEG1 is the precursor of perforin.
Identification of human perforin homologs
To search for perforin homologs, the human perforin protein sequence [Refseq:NP_005032.2] was used to query genome and protein databases on NCBI and ENSEMBL using both tBLASTn and BLASTp . Hits representing true perforin homologs were distinguished from other MACPF domain-containing proteins/genes in two ways. Firstly, full-length sequences were subjected to domain prediction in PFAM, any sequences not possessing both a MACPF domain and C2 domain were discarded. Secondly, partial sequences without sufficient data to reasonably contain both domains were used to query the non-redundant protein database by BLAST, any sequences with top hits to proteins other than perforin were discarded. In some cases, annotated protein sequences in ENSEMBL required amendment to fix non-canonical splice sites and obtain valid start and/or stop codons. A complete table of accessions used in this study is in Additional file 1: Table S1, and a complete file of amended protein sequences used in this study is in Additional file 2: Figure S1.
Genome scaffolds were viewed in ENSEMBL and NCBI and redrawn for ease of comparison. Orientation of genes were made relative to perforin, which was always represented in the positive orientation, except where it was necessary to be drawn in the negative orientation for comparison to other scaffolds in the figure.
Assembly of Ornithorhynchus anatinuscontigs
We exported the contig sequences containing perforin genes and used these as queries with megaBLAST on the trace archive (Ornithorhynchus anatinus – other) to search for end sequences of BAC and fosmid clones. We used the paired ends of matching clones to query the reference genome sequence to find matching contigs and by comparing the orientation of the clone hits we could infer the orientation of the various contigs. We assembled ten previously disconnected contigs in one instance, and two contigs in another instance, shown in Additional file 3: Figure S2. To search for genes on these contigs we viewed them in ENSEMBL. We found annotated partial genes that matched ADAMTS14 on some contigs, and where there were no annotated genes we searched for homologous protein sequences using tBLASTx and compiled these into a single protein sequence in the order suggested by our assembly.
Protein sequence alignment and phylogenetic trees
Fish perforin protein sequences were aligned using ClustalW and manually edited to minimise gaps and align conserved structural elements in bioedit (version 184.108.40.206) [31, 32]. The alignment is available in Additional file 4: Figure S3. The Bayesian inference tree was constructed using MrBayes (version 3.2.1) with the WAG amino acid substitution model, invariant sites and a gamma distribution [33, 34]. 100,000 generations were run with trees sampled every 100 generations and the final 50% majority rule tree was calculated after discarding the first 25% of trees as burnin. Trees were displayed using FigTree (version 1.3.1) .
To align the MACPF domains of perforin, C6 and MPEG1, domain boundaries were chosen based on three rounds of PSI-BLAST using full-length murine MPEG1 as the probe, together with information of the perforin and C8 structures [16, 36]. Sequences truncated to the MACPF boundaries were then aligned with ClustalW, positions with gaps were removed and phylogenetic trees were constructed with MrBayes as above. The alignment (including gaps) is available in Additional file 5: Figure S4. Phylogenetic trees have been deposited in TreeBASE and can be accessed here: http://purl.org/phylo/treebase/phylows/study/TB2:S12411.
Search criteria for perforin orthologs
Conservation of the human perforin locus extends to the last common ancestor (LCA) of Mammaliabut not to earlier vertebrate lineages
Perforin is present in Gnathostomata but not earlier Chordata, and many species have multiple perforin genes
Tracing perforin linked genes through bony vertebrate lineages we created maps of perforin loci and inferred the locus configurations in various major lineages (Figure 3). This analysis shows that the perforin gene has undergone a striking degree of duplication and repositioning within vertebrate genomes, being found at no fewer than 5 distinct loci across various lineages, and more if individual fish species are considered. The most conserved locus is the STX1b locus, which emerged before the LCA of Euteleostomi, as it is found in both teleost fish and tetrapods (Figure 3). This locus is present in amphibians and reptiles but absent in birds and therian mammals. It is also present in platypus, where we again linked two contigs using the approach described above (Figure 3, Additional file 3: Figure S2). The orientation of the O. anatinus perforin 1.3 gene relative to STX4 is opposite to other examples of this locus but this may be a result of incorrect assembly rather than a true gene inversion as there are many gaps in this short contig. The second contig contains the FUS gene, present at this locus in A. carolinensis which further confirms the arrangement of this locus in the platypus.
In Figure 4 we have summarised the numbers of full length genes, apparent pseudogenes and partial genes present in the available bony vertebrate genomes. All therian mammals have a single perforin gene and no pseudogenes, with the exception of the hedgehog (Erinaceus europaeus), which also has a full length pseudogene, and the hyrax (Procavia capensis) which also has a partial pseudogene. By contrast, most other lineages within Euteleostomi have species with multiple perforin genes. Birds appear to be an exception with evidence for only one perforin gene in each of chicken (partial), turkey (possible pseudogene, see discussion) and mallard duck (partial) genomes, and the aforementioned exception of zebra finch (no perforin gene). All sequenced teleost fish have at least three full length perforin genes and the frog Xenopus tropicalis has the largest number of perforin genes at 11. We investigated these cases of perforin gene duplication in more detail.
The LCA of Teleosteipossessed multiple perforin genes
We compiled all the full length perforin protein sequences from the five available fish genome sequences (Danio rerio, Takifugu rubripes, Tetraodon nigroviridis, Gasterosteus aculeatus and Oryzias latipes) as well as nine sequences from cloning projects from five other fish species (Paralichthys olivaceus, Oncorhynchus mykiss Salmo salar, Ctenopharyngodon idella and Carassius auratus langsdorfii) giving a total of 38 full length fish perforin proteins [38–41].
We also noted that each clade grouped genes found at syntenic loci. Two clades, with genes at PAQR4 and MPEG1 loci (termed here ‘PAQR4 clade’, ‘MPEG1 clade’ etc.), have representatives in all fish genome data available, demonstrating that they were present in the LCA of Teleostei (Figure 3). Others have representatives in only a few species, suggesting that they either arose via a later duplication or that they were present in the LCA of Teleostei but were subsequently lost in some lineages. We see evidence for both of these processes. For example (1) the MYADML2 clade is present in T. rubripes and G. aculeatus but no other species, suggesting it was acquired in the Percomorpha lineage and subsequently lost from T. nigroviridis and O. latipes; and (2) the STX1b clade is found in D. rerio, T. rubripes and T. nigroviridis as well as some tetrapods (showing that it must have been present in the LCA of Teleostei) but is lost from Smegmamorpha.
There are also nine perforin sequences from fish species without genome assemblies on this tree (black branches, Figure 5). This allowed us to make predictions about the loci of these genes based on their clustering. Of these, six fall in the PAQR4 clade and one in the STX1b clade. C. auratus langsdorfii perforin 1.3 clusters with a group of linked D. rerio perforins (1.2, 1.7 & 1.8) at the S100v2 locus, a locus which we originally believed to be D. rerio specific but this additional evidence suggests that it may have originated earlier, possibly in the LCA of the Cyprinidae family. The other sequence, S. salar perforin 1.2 does not make a strong cluster with any group and may represent a new locus, perhaps specific to the Atlantic salmon.
We were interested to know about the expression of these multiple fish perforin genes. We cross-referenced our collection of perforin genes with available EST data. In this way we found evidence that seventeen fish perforin genes are expressed in a wide range of tissues, and that more than one perforin gene is expressed in most species (Figure 5). Looking at how these expressed genes clustered we noticed that while all major clades contained genes that were expressed in one or more species, only for members of the MPEG1 clade was expression evident in all species examined. It is tempting to speculate that this locus contains a highly - or widely - expressed perforin gene (which would lead to a high probability of representation in EST datasets). Conflicting with this is the fact that most of the manually cloned perforin genes fall within the PAQR4 clade and none fall in the MPEG1 clade. In any case these findings raise the interesting question of why fish possess multiple perforin paralogs.
The additional perforin genes in X. tropicalisare recent, lineage specific duplications
The massive expansion of perforin genes in X. tropicalis intrigued us. As shown in Figure 3, nine of these are found on two scaffolds. Both of these are likely to be STX1b loci, supported by upstream gene information in comparison to A. carolinensis, and both scaffolds finish shortly after perforin genes and are missing downstream sequence that would confirm this hypothesis. The remaining two are all found alone on small scaffolds with no linked genes to provide context.
To further assess the origin of these multiple genes we looked at their exon/intron patterns. All of the nine full length genes in X. tropicalis contain a second CDS intron (phase 0) in addition to the completely conserved phase 2 intron (Figure 6). This intron is specific to this species and demonstrates conclusively that these multiple perforin genes are a result of lineage specific duplications, occurring subsequent to intron gain. Genome information from additional amphibians would allow us to narrow the timeframe of this event but at this stage we can only conclude that it occurred sometime after the divergence of Amphibia and Amniota. EST evidence exists for two full length and one partial gene, again raising the question of why some species need many perforin genes when one is sufficient in mammals.
Perforin most likely evolved from an MPEG1-like ancestor
Zebrafish complement, perforin and MPEG1 intron numbers
To further examine the evolutionary relationship between perforin, MPEG1 and complement proteins we curated and aligned these protein sequences from human, mouse, lizard and fugu, and produced a phylogenetic tree. The terminal complement components and perforin all contain an EGF-like domain after their MACPF domain. While MPEG1 has a cysteine rich region in a similar position, this is not believed to be an EGF-like domain . We therefore restricted our alignment to the MACPF domain, the boundaries of which were chosen based on PSI-BLAST results and structural data and then aligned with ClustalW. As shown in Figure 7b, the branch lengths of the tree show that the MACPF domain of perforin is more closely related to C6 than to MPEG1. This is also true of trees where C6 is substituted for C7 (data not shown).
Perforin is a critical protein in the granule-exocytosis pathway of CLs, therefore tracing the evolution of this gene yields insights into the evolution of the pathway itself (discussed below). We have characterised all of the perforin genes from available databases and find many species with multiple perforin genes. Classifying these duplicated genes led to the discovery of multiple gene duplication events in different lineages, and insights into the origins of perforin itself.
MPEG1: the precursor of perforin?
As no other proteins share the entire domain structure of perforin, the origin of the perforin gene is obscure. In terms of structure and function it has always been compared to C9 because this terminal complement component also polymerises and inserts into membranes, although unlike perforin it requires a primer (the C5b-C8 oligomer) to initiate membrane binding. Although the complement proteins appear to have descended from a common ancestral gene, a precursor of the perforin gene has not yet been identified .
We suggest here that perforin originated from the ancient gene MPEG1 on the basis of similar gene structure and chromosomal co-location of these genes in ancient bony vertebrate lineages. It is more likely that perforin originated from MPEG1 than from any other extant MACPF family member because no other genes encoding MACPF domain-containing proteins share these genetic characteristics. However, given that the MACPF domains of perforin and complement component C6 are more closely related than those of perforin and MPEG1, it is likely that perforin and MPEG1 are separated by an intermediate with features of both perforin and a terminal complement component. At present no available genome data contains remnants of such a gene.
The origins of MPEG1 itself are presently obscure, but its presence as an intronless gene in one of the most ancient metazoan lineages (e.g. Porifera), and absence in other eukaryotic kingdoms, coupled with the existence of MACPF proteins in prokaryotes, suggests a mechanism of horizontal gene transfer from prokaryotes to early metazoans.
Do birds possess functional perforin genes?
Perforin genes appear in all available bony vertebrate genomes with the surprising exception of zebra finch (Taeniopygia guttata). The presence of other genes of the granule-exocytosis pathway in zebra finch (granzyme A), as well as perforin’s absolute conservation in Euteleostomi, predicts that any tetrapod without a perforin gene would be severely immunocompromised, yet there are no reports of zebra finches having such a defect. We therefore suggest that this species does indeed have a perforin gene but it has not yet been uncovered by sequencing projects. The turkey (Meleagris gallopavo) perforin gene also appears to be affected by a sequencing error; the gene contains a frameshift caused by a single nucleotide deletion. ENSEMBL corrects this by modelling a frameshift intron (2 bp) to maintain the coding sequence, but the possibility remains that this perforin gene in turkey is a pseudogene, in which case the species would be immunocompromised unless a functional paralog has been missed. The best studied bird species, the chicken (G. gallus) has only a partial perforin gene sequence annotated, with a gap in the assembly where the remainder of the gene should be found. Chicken perforin expression has been assessed in the context of Marek’s disease, where mRNA was shown by real time PCR to be upregulated in the spleens of infected chickens, although another study involving infectious bursal disease saw no upregulation in bursal mononuclear cells [44, 45]. This suggests that at least one bird species has a functional perforin gene but ultimately more data, particularly EST or biological data, is required to make firm conclusions about the perforin gene in birds.
Multiple perforin genes for multiple functions?
It is well established that a whole-genome duplication occurred in the Teleostei lineage, and that many genomic rearrangements are possible when reverting from a tetraploid to a diploid state [46, 47]. Duplications of genes functioning in the immune system have been noted before in teleost fish, for instance an expansion of the CC chemokine genes has been described . However perforin is an extreme example with a single gene in the Gnathostomata LCA giving rise to nine in Danio rerio. Paralogs are likely to have arisen from the whole-genome duplication and from local and distal segmental duplications. Multiple gene loss as well as intron gain and loss events have also occurred, highlighting the instability of the perforin gene in the teleost fish.
One of the most interesting questions that the duplicated fish perforin genes poses is why fish require so many, when in therian mammals a single gene is sufficient for the function of cytotoxic lymphocytes. Simple explanations of this phenomenon are that perforin has evolved additional molecular functions in fish (unlikely as the protein domains are highly conserved), and/or it has a broader tissue distribution and role or a more complex regulatory pattern than in mammals. For example, perforin is essentially restricted to two hematopoetic cell types in mammals (CTL and NK cells) and is controlled by similar circuitry. In fish it may be present in three or more immune cell types, with each paralog being restricted to a particular type. Alternatively, all perforin paralogs may be present in the same cell type, with each responding to distinct developmental or environmental signals. Indeed one study characterizing perforin from C. auratus langsdorfii notes that only one of the three perforin genes cloned was up-regulated by allo-reactive stimulation as measured by real time PCR, suggesting that the other two are not important for CTL cytotoxicity . However it has been shown in mice that the perforin protein can be translated on demand from a stored pool of mRNA so increased mRNA expression may not be required to increase levels of perforin protein . Nevertheless, the study by Nakanishi and colleagues, as well a wealth of EST data which we have mapped to our collection of perforin genes, demonstrates that multiple perforin genes per species are expressed at least at the mRNA level. Being present at distant loci, it is likely that these paralogs would be controlled by different regulatory elements, and therefore may have different tissue distribution. The role of these extra perforins will remain obscure until molecular and cellular studies examine their functions.
Emergence of the perforin-mediated cytotoxic pathway
Cell-mediated cytotoxicity has been described in many species but the mechanisms of killing are not well understood except in mammals. Invertebrates can determine self-non-self and NK-like cells have been implicated in rejecting non-self cells . In addition, a primitive form of adaptive immunity is present in the lamprey, a jawless vertebrate. This species possesses T cell-like lymphocytes that would be expected to kill target cells . Our data suggests that ancient cytotoxic cells from species earlier than Gnathostomata are unlikely to have an active granule-exocytosis/perforin pathway, so invertebrate NK-like cells and lamprey CTL therefore must kill their targets in another fashion. This pathway seems to have arisen concurrently with the MHC-TCR antigen presentation system, suggesting primordial CTLs possessed perforin, whereas NKs acquired this cytotoxic machinery after their inception.
The pore-forming protein perforin is the only component that is absolutely required for the granule-exocytosis pathway that cytotoxic lymphocytes deploy to eliminate deleterious cells. The perforin gene is present in Gnathostomata but not earlier species. Perforin evolved from a duplication of the related gene MPEG1, and shares a common ancestor with the terminal complement components. Surprisingly, we find that most bony vertebrate species predating placental mammals have multiple perforin genes, of unknown function. These findings also indicate that cytotoxic cells from invertebrates and jawless vertebrates must use alternative proteins or pathways to kill their target cells.
= Major histocompatibility complex
= T-cell receptor
= Macrophage expressed gene 1
= Cytotoxic lymphocyte
= Natural kill cell
= Cytotoxic T lymphocyte
= Familial hemophagocytic lymphohistiocytosis type 2
= Membrane attack complex
= Membrane attack complex/perforin
= Complement component 6
= Epidermal growth factor-like
= Untranslated region
= Expressed sequence tag
= Coding sequence
= Last common ancestor.
We thank B.T. Porebski and the Orchard supercomputer at Monash University for computational resources and assistance. This work was supported by the National Health and Medical Research Council of Australia and the Australian Research Council. M.E.D. is a recipient of the Australian Postgraduate Award, M.A.D. is an NHMRC Career Development Fellow, J.C.W. is an ARC Federation Fellow and honorary NHMRC Principal Research Fellow.
- Kagi D, Vignaux F, Ledermann B, Burki K, Depraetere V, Nagata S, Hengartner H, Golstein P: Fas and perforin pathways as major mechanisms of T cell-mediated cytotoxicity. Science. 1994, 265 (5171): 528-530. 10.1126/science.7518614.PubMedView ArticleGoogle Scholar
- Yannelli J, Sullivan J, Mandell G, Engelhard V: Reorientation and fusion of cytotoxic T lymphocyte granules after interaction with target cells as determined by high resolution cinemicrography. J Immunol. 1986, 136 (2): 377-382.PubMedGoogle Scholar
- Podack ER, Konigsberg PJ: Cytolytic T cell granules. Isolation, structural, biochemical, and functional characterization. J Exp Med. 1984, 160 (3): 695-710. 10.1084/jem.160.3.695.PubMedView ArticleGoogle Scholar
- Masson D, Nabholz M, Estrade C, Tschopp J: Granules of cytolytic T-lymphocytes contain two serine esterases. EMBO J. 1986, 5 (7): 1595-1600.PubMedPubMed CentralGoogle Scholar
- Azumi K, De Santis R, De Tomaso A, Rigoutsos I, Yoshizaki F, Pinto MR, Marino R, Shida K, Ikeda M, Arai M, et al: Genomic analysis of immunity in a Urochordate and the emergence of the vertebrate immune system: “waiting for Godot”. Immunogenetics. 2003, 55 (8): 570-581. 10.1007/s00251-003-0606-5.PubMedView ArticleGoogle Scholar
- Kasahara M, Suzuki T, Pasquier LD: On the origins of the adaptive immune system: novel insights from invertebrates and cold-blooded vertebrates. Trends Immunol. 2004, 25 (2): 105-111. 10.1016/j.it.2003.11.005.PubMedView ArticleGoogle Scholar
- Pancer Z, Amemiya CT, Ehrhardt GRA, Ceitlin J, Larry Gartland G, Cooper MD: Somatic diversification of variable lymphocyte receptors in the agnathan sea lamprey. Nature. 2004, 430 (6996): 174-180. 10.1038/nature02740.PubMedView ArticleGoogle Scholar
- Khalturin K, Becker M, Rinkevich B, Bosch TC: Urochordates and the origin of natural killer cells: identification of a CD94/NKR-P1-related receptor in blood cells of Botryllus. Proc Natl Acad Sci U S A. 2003, 100 (2): 622-627. 10.1073/pnas.0234104100.PubMedPubMed CentralView ArticleGoogle Scholar
- Engelmann P, Kiss J, Csongei V, Cooper EL, Nemeth P: Earthworm leukocytes kill HeLa, HEp-2, PC-12 and PA317 cells in vitro. J Biochem Biophys Methods. 2004, 61 (1–2): 215-227.PubMedView ArticleGoogle Scholar
- Stepp SE, Dufourcq-Lagelouse R, Le Deist F, Bhawan S, Certain S, Mathew PA, Henter J-I, Bennett M, Fischer A, de Saint Basile G, et al: Perforin Gene Defects in Familial Hemophagocytic Lymphohistiocytosis. Science. 1999, 286 (5446): 1957-1959. 10.1126/science.286.5446.1957.PubMedView ArticleGoogle Scholar
- Voskoboinik I, Thia MC, Fletcher J, Ciccone A, Browne K, Smyth MJ, Trapani JA: Calcium-dependent plasma membrane binding and cell lysis by perforin are mediated through its C2 domain: A critical role for aspartate residues 429, 435, 483, and 485 but not 491. J Biol Chem. 2005, 280 (9): 8426-8434.PubMedView ArticleGoogle Scholar
- Hadders MA, Beringer DX, Gros P: Structure of C8alpha-MACPF reveals mechanism of membrane attack in complement immune defense. Science. 2007, 317 (5844): 1552-1554. 10.1126/science.1147103.PubMedView ArticleGoogle Scholar
- Rosado CJ, Buckle AM, Law RH, Butcher RE, Kan WT, Bird CH, Ung K, Browne KA, Baran K, Bashtannyk-Puhalovich TA, et al: A common fold mediates vertebrate defense and bacterial attack. Science. 2007, 317 (5844): 1548-1551. 10.1126/science.1144706.PubMedView ArticleGoogle Scholar
- Slade DJ, Lovelace LL, Chruszcz M, Minor W, Lebioda L, Sodetz JM: Crystal structure of the MACPF domain of human complement protein C8 alpha in complex with the C8 gamma subunit. J Mol Biol. 2008, 379 (2): 331-342. 10.1016/j.jmb.2008.03.061.PubMedPubMed CentralView ArticleGoogle Scholar
- Baran K, Dunstone M, Chia J, Ciccone A, Browne KA, Clarke CJP, Lukoyanova N, Saibil H, Whisstock JC, Voskoboinik I, et al: The Molecular Basis for Perforin Oligomerization and Transmembrane Pore Assembly. Immunity. 2009, 30 (5): 684-695. 10.1016/j.immuni.2009.03.016.PubMedView ArticleGoogle Scholar
- Law RH, Lukoyanova N, Voskoboinik I, Caradoc-Davies TT, Baran K, Dunstone MA, D’Angelo ME, Orlova EV, Coulibaly F, Verschoor S, et al: The structural basis for membrane binding and pore formation by lymphocyte perforin. Nature. 2010, 468 (7322): 447-451. 10.1038/nature09518.PubMedView ArticleGoogle Scholar
- Rossjohn J, Feil SC, McKinstry WJ, Tweten RK, Parker MW: Structure of a cholesterol-binding, thiol-activated cytolysin and a model of its membrane form. Cell. 1997, 89 (5): 685-692. 10.1016/S0092-8674(00)80251-2.PubMedView ArticleGoogle Scholar
- Rosado CJ, Kondos S, Bull TE, Kuiper MJ, Law RH, Buckle AM, Voskoboinik I, Bird PI, Trapani JA, Whisstock JC, et al: The MACPF/CDC family of pore-forming toxins. Cell Microbiol. 2008, 10 (9): 1765-1774. 10.1111/j.1462-5822.2008.01191.x.PubMedPubMed CentralView ArticleGoogle Scholar
- Kondos SC, Hatfaludi T, Voskoboinik I, Trapani JA, Law RH, Whisstock JC, Dunstone MA: The structure and function of mammalian membrane-attack complex/perforin-like proteins. Tissue Antigens. 2010, 76 (5): 341-351. 10.1111/j.1399-0039.2010.01566.x.PubMedView ArticleGoogle Scholar
- Bhakdi S, Tranum-Jensen J: Molecular nature of the complement lesion. Proc Natl Acad Sci U S A. 1978, 75 (11): 5655-5659. 10.1073/pnas.75.11.5655.PubMedPubMed CentralView ArticleGoogle Scholar
- Podack ER, Tschoop J, Muller-Eberhard HJ: Molecular organization of C9 within the membrane attack complex of complement. Induction of circular C9 polymerization by the C5b-8 assembly. J Exp Med. 1982, 156 (1): 268-282. 10.1084/jem.156.1.268.PubMedPubMed CentralView ArticleGoogle Scholar
- Dourmashkin RR, Deteix P, Simone CB, Henkart P: Electron microscopic demonstration of lesions in target cell membranes associated with antibody-dependent cellular cytotoxicity. Clin Exp Immunol. 1980, 42 (3): 554-560.PubMedPubMed CentralGoogle Scholar
- Podack ER, Dennert G: Assembly of two types of tubules with putative cytolytic function by cloned natural killer cells. Nature. 1983, 302 (5907): 442-445. 10.1038/302442a0.PubMedView ArticleGoogle Scholar
- Young JD, Cohn ZA, Podack ER: The ninth component of complement and the pore-forming protein (perforin 1) from cytotoxic T cells: structural, immunological, and functional similarities. Science. 1986, 233 (4760): 184-190. 10.1126/science.2425429.PubMedView ArticleGoogle Scholar
- Spilsbury K, O’Mara MA, Wu WM, Rowe PB, Symonds G, Takayama Y: Isolation of a novel macrophage-specific gene by differential cDNA analysis. Blood. 1995, 85 (6): 1620-1629.PubMedGoogle Scholar
- Wiens M, Korzhev M, Krasko A, Thakur NL, Perovic-Ottstadt S, Breter HJ, Ushijima H, Diehl-Seifert B, Muller IM, Muller WE: Innate immune defense of the sponge Suberites domuncula against bacteria involves a MyD88-dependent signaling pathway. Induction of a perforin-like molecule. J Biol Chem. 2005, 280 (30): 27949-27959. 10.1074/jbc.M504049200.PubMedView ArticleGoogle Scholar
- Hemmrich G, Miller DJ, Bosch TC: The evolution of immunity: a low-life perspective. Trends Immunol. 2007, 28 (10): 449-454. 10.1016/j.it.2007.08.003.PubMedView ArticleGoogle Scholar
- Podack ER, Deyev V, Shiratsuchi M: Pore formers of the immune system. Adv Exp Med Biol. 2007, 598: 325-341. 10.1007/978-0-387-71767-8_23.PubMedView ArticleGoogle Scholar
- He X, Zhang Y, Yu Z: An Mpeg (macrophage expressed gene) from the Pacific oyster Crassostrea gigas: molecular characterization and gene expression. Fish Shellfish Immunol. 2011, 30 (3): 870-876. 10.1016/j.fsi.2011.01.009.PubMedView ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.PubMedView 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 (22): 4673-4680. 10.1093/nar/22.22.4673.PubMedPubMed CentralView ArticleGoogle Scholar
- Hall TA: 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
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19 (12): 1572-1574. 10.1093/bioinformatics/btg180.PubMedView ArticleGoogle Scholar
- Whelan S, Goldman N: A general empirical model of protein evolution derived from multiple protein families using a maximum-likelihood approach. Mol Biol Evol. 2001, 18 (5): 691-699. 10.1093/oxfordjournals.molbev.a003851.PubMedView ArticleGoogle Scholar
- Rambaut A: FigTree: Tree Figure Drawing Tool Version 1.3.1. 2009, Software available at, [http://tree.bio.ed.ac.uk/software/figtree/]Google Scholar
- Lovelace LL, Cooper CL, Sodetz JM, Lebioda L: Structure of human C8 protein provides mechanistic insight into membrane pore formation by complement. J Biol Chem. 2011, 286 (20): 17585-17592. 10.1074/jbc.M111.219766.PubMedPubMed CentralView ArticleGoogle Scholar
- Kumar S, Hedges SB: A molecular timescale for vertebrate evolution. Nature. 1998, 392 (6679): 917-920. 10.1038/31927.PubMedView ArticleGoogle Scholar
- Hwang JY, Ohira T, Hirono I, Aoki T: A pore-forming protein, perforin, from a non-mammalian organism, Japanese flounder, Paralichthys olivaceus. Immunogenetics. 2004, 56 (5): 360-367.PubMedView ArticleGoogle Scholar
- Athanasopoulou S, Marioli D, Mikrou A, Papanastasiou AD, Zarkadis IK: Cloning and characterization of the trout perforin. Fish Shellfish Immunol. 2009, 26 (6): 908-912. 10.1016/j.fsi.2009.03.021.PubMedView ArticleGoogle Scholar
- Leong JS, Jantzen SG, von Schalburg KR, Cooper GA, Messmer AM, Liao NY, Munro S, Moore R, Holt RA, Jones SJ, et al: Salmo salar and Esox lucius full-length cDNA sequences reveal changes in evolutionary pressures on a post-tetraploidization genome. BMC Genomics. 2010, 11: 279-10.1186/1471-2164-11-279.PubMedPubMed CentralView ArticleGoogle Scholar
- Toda H, Araki K, Moritomo T, Nakanishi T: Perforin-dependent cytotoxic mechanism in killing by CD8 positive T cells in ginbuna crucian carp, Carassius auratus langsdorfii. Dev Comp Immunol. 2011, 35 (1): 88-93. 10.1016/j.dci.2010.08.010.PubMedView ArticleGoogle Scholar
- Mondragon-Palomino M, Pinero D, Nicholson-Weller A, Laclette JP: Phylogenetic analysis of the homologous proteins of the terminal complement complex supports the emergence of C6 and C7 followed by C8 and C9. J Mol Evol. 1999, 49 (2): 282-289. 10.1007/PL00006550.PubMedView ArticleGoogle Scholar
- Suzuki MM, Satoh N, Nonaka M: C6-like and C3-like molecules from the cephalochordate, amphioxus, suggest a cytolytic complement system in invertebrates. J Mol Evol. 2002, 54 (5): 671-679. 10.1007/s00239-001-0068-z.PubMedView ArticleGoogle Scholar
- Sarson AJ, Abdul-Careem MF, Read LR, Brisbin JT, Sharif S: Expression of cytotoxicity-associated genes in Marek’s disease virus-infected chickens. Viral Immunol. 2008, 21 (2): 267-272. 10.1089/vim.2007.0094.PubMedView ArticleGoogle Scholar
- Rauf A, Khatri M, Murgia MV, Saif YM: Expression of perforin-granzyme pathway genes in the bursa of infectious bursal disease virus-infected chickens. Dev Comp Immunol. 2011, 35 (5): 620-627. 10.1016/j.dci.2011.01.007.PubMedView ArticleGoogle Scholar
- Postlethwait JH: The zebrafish genome in context: ohnologs gone missing. J Exp Zool B Mol Dev Evol. 2007, 308 (5): 563-577.PubMedView ArticleGoogle Scholar
- Jaillon O, Aury JM, Brunet F, Petit JL, Stange-Thomann N, Mauceli E, Bouneau L, Fischer C, Ozouf-Costaz C, Bernot A, et al: Genome duplication in the teleost fish Tetraodon nigroviridis reveals the early vertebrate proto-karyotype. Nature. 2004, 431 (7011): 946-957. 10.1038/nature03025.PubMedView ArticleGoogle Scholar
- Peatman E, Liu Z: Evolution of CC chemokines in teleost fish: a case study in gene duplication and implications for immune diversity. Immunogenetics. 2007, 59 (8): 613-623. 10.1007/s00251-007-0228-4.PubMedView ArticleGoogle Scholar
- Fehniger TA, Cai SF, Cao X, Bredemeyer AJ, Presti RM, French AR, Ley TJ: Acquisition of murine NK cell cytotoxicity requires the translation of a pre-existing pool of granzyme B and perforin mRNAs. Immunity. 2007, 26 (6): 798-811. 10.1016/j.immuni.2007.04.010.PubMedView ArticleGoogle Scholar
- Muller WE, Muller IM: Origin of the metazoan immune system: identification of the molecules and their functions in sponges. Integr Comp Biol. 2003, 43 (2): 281-292. 10.1093/icb/43.2.281.PubMedView ArticleGoogle Scholar
- Guo P, Hirano M, Herrin BR, Li J, Yu C, Sadlonova A, Cooper MD: Dual nature of the adaptive immune system in lampreys. Nature. 2009, 459 (7248): 796-801. 10.1038/nature08068.PubMedPubMed CentralView 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.