## Abstract

Classical T cells, those with αβ T-cell receptors (TCRs), are an important component of the dominant paradigm for self-nonself immune recognition in vertebrates. αβ T cells recognize foreign peptide antigens when they are bound to MHC molecules on the surfaces of antigen-presenting cells. γδ T cells bear a similar receptor, and it is often assumed that these T cells also require specialized antigen-presenting molecules for immune recognition, which we term “indirect antigen recognition.” B-cell receptors, or immunoglobulins, bind directly to antigens without the help of a specialized antigen-presenting molecule. Phylogenetically, it has been assumed that T-cell receptors and the genes that encode them are a monophyletic group, and that “indirect” antigen recognition evolved before the split into two types of TCR. Recently, however, it has been proposed that γδ-TCRs bind directly to antigens, as do immunoglobulins (Ig’s). This calls into question the null hypothesis that indirect antigen recognition is a common characteristic of TCRs and, by extension, the hypothesis that all TCR gene sequences form a monophyletic group. To determine whether alternative explanations for antigen recognition and other historical relationships among TCR genes might be possible, we performed phylogenetic analyses on amino acid sequences of the constant and variable regions which encode the basic subunits of TCR and Ig molecules. We used both maximum-parsimony and genetic distance-based methods and could find no strong support for the hypothesis of TCR monophyly. Analyses of the constant region suggest that TCR γ or δ sequences are the most ancient, implying that the ancestral immune cell was like a modern γδ T cell. From this γδ-like ancestor arose αβ T cells and B cells, implying that indirect antigen recognition is indeed a derived property of αβ-TCRs. Analyses of the variable regions are complicated by strong selection on antigen-binding sequences, but imply that direct antigen binding is the ancestral condition.

## Introduction

In the classical paradigm for T-cell-mediated immune recognition in mammals, an αβ T-cell receptor (TCR) on the surface of a T cell recognizes and binds to a peptide antigen that is bound to a major histocompatibility complex (MHC) molecule on the surface of a specialized antigen-presenting cell (APC). When a T cell encounters a peptide-MHC complex in which either component is foreign, it can respond appropriately, for instance, by lysing the target or by providing help for B-cell antibody responses. However, not all T cells carry the same type of TCR molecule, and only αβ T cells react to the MHC+peptide complex. A second class of T cells, called γδ T cells, is less well understood, perhaps because in humans, these T cells comprise only 1%–5% of peripheral blood lymphocytes (Haas, Pereira, and Tonegawa 1993<$REFLINK> ). However, in other human tissues, especially gut and reproductive epithelia, or in the peripheral blood of other mammals such as ungulates, they may comprise almost half of the T cells (Haas, Pereira, and Tonegawa 1993<$REFLINK> ). The role of γδ T cells in the vertebrate immune system is only just beginning to be understood (Brenner, Strominger, and Krangel 1988<$REFLINK> ; Raulet 1989<$REFLINK> ; Allison and Havran 1991<$REFLINK> ; Kaufmann 1996<$REFLINK> ), but comparative studies in immunology point to an important role in primary defense against infection (Flajnik 1998<$REFLINK> ). Despite this growing body of evidence, it is still sometimes assumed that γδ-TCRs function similarly to αβ-TCRs (Bartl, Baltimore, and Weissman 1994<$REFLINK> ; Flajnik 1998<$REFLINK> ), recognizing antigens only when they are presented by some sort of specialized antigen-presenting molecule analogous to an MHC molecule (we refer to this as “indirect antigen binding”). However, recent studies (Schild et al. 1994<$REFLINK> ; Rock et al. 1994<$REFLINK> ; Chien and Jores 1995<$REFLINK> ) indicate that the process of antigen recognition by γδ T cells may be fundamentally different from that in αβ T cells. In fact, γδ-TCRs appear to recognize antigens in a manner similar to the antigen recognition processes of immunoglobulins, the receptors of B cells. Immunoglobulins bind directly to antigens and do not require specialized antigen processing and presentation as do αβ T cells.

TCRs and immunoglobulins (Ig’s) are the most closely related members of a large protein family called the immunoglobulin superfamily (IGSF), many of whose members are involved in immune recognition or cellular adhesion (Hunkapiller and Hood 1989<$REFLINK> ). TCR and Ig molecules are clearly more similar to each other than to other IGSF molecules (Marchalonis, Schluter, and Edmundson 1997<$REFLINK> ), and the genes that encode them display highly similar patterns of organization and rearrangement during transcription and translation (Hunkapiller and Hood 1989<$REFLINK> ) that are not observed in any other genes. In phylogenetic terms, αβ- and γδ-TCR molecules and genes are often thought of as a monophyletic group whose sister group is the Ig’s, and this assumption is one reason that γδ- and αβ-TCRs have been assumed to function similarly. As illustrated in figure 1A, since αβ- TCRs exhibit indirect antigen binding, and since they are closely related to the γδ-TCRs, it has been assumed that this mode of antigen binding evolved before the split between the two. However, it would be equally parsimonious to assume that indirect antigen binding is a derived characteristic of αβ T cells (fig. 1B ). Moreover, it is possible that the phylogenetic relationships in figure 1 are incorrect, specifically, that the TCR genes are not monophyletic. This would lead to different conclusions about the evolution of immune cells and antigen recognition and binding by αβ-TCRs, γδ-TCRs, and also Ig’s. We decided to examine the immunological null hypothesis that the antigen recognition functions of αβ and γδ T cells must be similar by constructing phylogenies of TCR and Ig genes in order to reconstruct the phylogenetic relationships among these molecules. To do so requires phylogenetic studies based on outgroup analysis, because we need to investigate the branching order of the TCR-α, -β, -γ, and -δ sequences, as well as how the TCRs as a whole relate to the immunoglobulin sequences as a whole. Previous studies (e.g., Bernstein et al. 1996<$REFLINK> ; Greenberg et al. 1995, 1996<$REFLINK> ; Schluter, Bernstein, and Marchalonis 1997<$REFLINK> ) have investigated relationships of TCR and Ig genes, but they included few TCR sequences and usually have not included outgroup analysis, precluding inferences about branching order among all six TCR and Ig types.

## Materials and Methods

Amino acid sequences for TCR-α, -β, -γ, and -δ chains, Ig heavy (IgH) and light (IgL) chains (only λ chains were included), and a recently discovered novel antigen receptor in sharks were retrieved from the GenBank, GenPept, and PIR databases (table 1 ). Although TCRs and Ig’s are known to exist in all jawed vertebrates, GenBank’s representation of the various types of sequences among various animal taxa is uneven. The current data set (compiled in January 1998) represents an effort to include as many animal species as possible without undue focus on overrepresented groups such as humans, apes, and rodents. Chondrichthyan sequences are included among all four types of TCR and both types of immunoglobulin (table 1 ).

TCRs and Ig’s are composed of subunits which are themselves composed of homologous regions known as the constant (C) and variable (V) regions. Generally, the C regions comprise the backbone of the receptor molecule, while the V regions include the sites of actual antigen binding. The C and V regions may be under very different selective pressures and could have very different evolutionary histories and so are analyzed separately. IgH’s contain differing numbers of constant regions, depending on the type of Ig molecule they form (e.g., IgM, IgG, etc.). To minimize the effects of such structural variation, the constant regions of IgH sequences are represented by the most membrane-proximal C region, usually the C-4 region of IgM, which is thought to be homologous across taxa (Fellah et al. 1992<$REFLINK> ). The constant and variable regions are joined together by shorter sequences called joining (J) and diversity (D) regions; thus, TCRs and Ig’s are heterodimers composed of one VJC and one VDJC molecule. The J and D regions are highly variable even among TCR lineages of a single individual (Schatz, Oettinger, and Schlissel 1992<$REFLINK> ) and so are unlikely to be phylogenetically useful.

Sequences were aligned using the program PILEUP (Feng and Doolittle 1987<$REFLINK> ) in the GCG package, version 8.1-UNIX (Devereaux, Haeberli, and Smithies 1984<$REFLINK> ). The program CLUSTAL W (Higgins and Sharp 1988<REFLINK> ) was also used and gave similar results, but results presented here are based primarily on the PILEUP alignments. Various gap penalties were tried. A lower gap penalty results in the insertion of more gaps in the alignments and greater similarity among sequences. An extremely low gap penalty tends to eliminate any phylogenetic signal, since aligned sequences are either highly similar or full of missing characters. An excessively high gap penalty is likely to lead to misalignment in cases in which insertions and deletions have occurred. A gap penalty of 2.0 (using PILEUP) guaranteed perfect alignment of two highly conserved motifs common to all TCR and Ig sequences, so this gap penalty was adopted for all alignments. The two motifs are the highly conserved WYRK/Q and YFCA motifs in the second and third framework regions of the V region genes. Overall, then, the alignments were conservative. The V region alignment, including the complementarity determining regions (CDRs), was truncated at a consensus length of 134 amino acids, ending at the highly conserved G-G hinge region between the V and D/J regions in the fourth framework region. Excluding the CDRs resulted in a consensus length of only 64. The C region alignment was truncated at a consensus length of 150 amino acids, the approximate length of the IgL C regions. Two different tree-building methods were used. The parsimony criterion was used to find the set of trees that minimizes the number of changes in character states across the whole tree, ignoring invariant or uninformative characters. Using the program PAUP, version 3.1 (Swofford 1993<REFLINK> ), parsimony trees were found by heuristic search with the branch-and-bound search option and 1,000 iterations. Each iteration starts the search with a different tree (i.e., set of branching patterns) and finds the shortest set of trees based on that starting point. Through dozens of trials and hundreds of hours of tree building with this method, we found that the set of shortest trees was almost always found in the first 10 iterations. Successive approximations character weighting (SACW; Carpenter 1988<$REFLINK> ) was used to further improve tree resolution. Each character was reweighted by its goodness of fit to the most parsimonious trees, and the trees were recalculated. The process was repeated until a constant tree topology was obtained. Thus, each parsimony consensus tree should represent the best fit between tree topology and the current data set. The second tree-building method was based on the use of genetic distance matrices and the neighbor-joining method of Saitou and Nei (1987)<$REFLINK> . These trees were calculated using the PHYLIP (version 3.57c) programs PROTDIST and NEIGHBOR (Felsenstein 1993<$REFLINK> ). One thousand neighbor-joining trees were calculated based on a bootstrap of 1,000 separate genetic distance matrices; the final trees are consensuses of these 1,000 trees. Trees based on outgroup comparison can be used to infer both patterns of relationship among taxa and branching order over time. Several members of the IGSF have been proposed to be representative of a primordial IGSF molecule which probably had only one Ig-like region that would have been ancestral to both the C and the V regions of TCRs and Ig’s (Hunkapiller and Hood 1989<$REFLINK> ). One of these is the TCR-associated cell surface molecule CD3. CD3 exists in three forms, two of which, CD3-γ and CD3-δ, were used as outgroup sequences. Several other molecules were also tested as outgroups, including the cell surface molecules CD3-ϵ, Thy-1, CD4, and CD8 (Williams 1987<REFLINK> ), but this did not result in improved resolution of TCR and Ig topologies. ## Results ### Constant Region Trees The C region alignment was 150 characters long. Of these, 143 characters were informative. Using the parsimony criterion, 73 C region trees of length 2,344 were found. These were reduced to 2 trees of length 2,345 using SACW (fig. 2 ). In these trees, the earliest branches are represented by TCR sequences, specifically TCR-γ. Of the four TCR groups, only TCR-α is monophyletic. The three chondrichthyan TCR-δ sequences appear to be more similar to each other than to their tetrapod homologs, which are associated with IgL. IgH sequences are nested well within the IgL sequences, implying that the light chains are older and gave rise to the heavy chains. The IgL sequences as a whole are polyphyletic, with the chondrichthyan group I (Hfl141, Sktigcve), II (Cpligl, Sktigcvc), and III (Hefigcve) IgL sequences forming the earliest branches. The IgH sequences comprise a monophyletic group that includes shark Novel Antigen Receptor (NAR, Gcu18701), IgNARC (Gcu51450), and IgW (Cpu40560). A second set of trees was calculated based on a bootstrap analysis of genetic distances (fig. 3 ). These trees are probably superior to the parsimony trees, since they include the following monophyletic groups: TCR-α/δ, TCR-β, TCR-γ, and all immunoglobulins. As in the previous trees, the earliest branches are TCR sequences, although in this analysis the oldest type is TCR-α/δ. The shark NAR, IgNARC, and IgW clearly are closely associated with the IgH sequences. The IgL sequences are polyphyletic, with osteichthyan and tetrapod IgL representing branches completely separate from the chondrichthyan group I, II, and III light chains. ### Variable Region The variable regions of T-cell receptors and immunoglobulins are divided into framework (FR) and complementarity-determining (CDR) regions. The latter are highly variable, conferring epitope-specific sensitivity on the receptor, and a priori are expected not to be phylogenetically informative. Therefore, the CDRs were excluded, leaving a V region alignment with a consensus length of 64 characters, of which 55 were informative. Using the parsimony criterion, an initial set of 119 trees of length 894 steps was found. These were reduced to 2 trees of length 896 using SACW (fig. 4 ). The IgH sequences, including IgW (Cpu40560), comprise the basal branch, giving rise to the TCRs as a whole, which in turn eventually give rise to the IgL group. IgL and TCR sequences are more closely related to each other than either is to the IgH sequences. None of the TCR groups is monophyletic. The position of axolotl TCR-β (Amttcrb) is problematic, rendering the IgL sequences paraphyletic. We note that if a second amphibian sequence, Xenopus TCR-β (U60436; Chretien et al. 1997<REFLINK> ), is added to the analysis, the axolotl sequence remains associated with IgL, although the Xenopus sequence groups with all other TCR-β sequences (data not shown). If Amttcrb is removed, TCR-β and IgL become monophyletic (data not shown). However, the tetrapod IgL sequences (those not including the chondrichthyan group I, II, and III sequences; table 1 ) do comprise a monophyletic group. The position of IgNARC (Gcu51450) is unclear: while it is probably associated with the IgH group, it appears as a separate branch intermediate between the IgH and TCR sequences. However, the shark NAR is clearly associated with TCR-γ and embedded in the TCR+IgL group.

Based on genetic distances, the neighbor-joining consensus tree (fig. 5 ) also indicates that IgH make up the basal branch and that TCR+IgL compose a monophyletic group. However, this tree differs from the parsimony tree in important respects. As was the case for the C region trees, the distance-based trees are probably superior to the parsimony-based trees, since they include the following monophyletic groups: IgH (including IgW and IgNARC), all IgL, TCR-α/δ, and TCR-γ. Perhaps oddly, the TCR-β sequences are polyphyletic, splitting into two groups, one comprising rabbit (Rabtcbxb), human (Humtcbyy), and axolotl (Amttcrb) sequences, and the other comprising shark (Hfu07624 and Reu75769), mouse (Mustcbxh), and chicken (Chktcrbc) sequences. This may be an artificial consequence of attraction of the Rabtcbxb and Humtcbyy sequences to the problematical sequence Amttcrb (see above), which we found often appeared in odd places on different trees (we did not exclude it because we have no objective grounds for believing this sequence to be odd, and because the addition of the Xenopus TCR-β sequence, as noted above, does not solve the problem [data not shown]). Polyphyly of the TCR-β sequences was also found when the CDR regions were not excluded, using both the maximum-parsimony and the genetic distance-based methods (trees not shown). The positions of IgW and IgNARC are clearly embedded within the IgH groups. However, as in the parsimony tree, NAR is associated with the TCR(+IgL) sequences.

## Discussion

### γδ T Cells as Primordial Immune Cells

Neither the V nor the C region trees support the hypothesis of monophyly of the TCR genes, but they disagree as to the overall branching patterns of the four TCR and two Ig sequence sets. Trees based on C region sequences suggest that TCR sequences evolved earliest, that these gave rise to IgL sequences, and that the last to appear were the IgH sequences. In marked contrast, trees based on V region sequences suggest that IgH sequences are the oldest sequences, followed by the TCRs, and then the IgL. These two hypotheses of relationship among TCR and Ig sequences might seem to be mutually exclusive were it not that the C and V region genes are probably subject to quite different selective regimes, due to the different roles of the C and V regions in the molecule. Topological relationships among the V region sequences could be more reflective of functional similarity than of historical descent, due to convergence of the types of sequences required for direct (as in the Ig’s) versus indirect (as in the αβ-TCRs) antigen binding. By the same token, the C region sequences, which are not involved in antigen recognition and binding, might be less likely to be convergent due to such constraints, and their relationships might be more reflective of phylogenetic history. Furthermore, the evolutionary history of the gene sequences does not necessarily reflect their current roles as members of modern TCR and Ig receptors, since their functions in modern molecules could have arisen after the genes themselves began to diverge.

If, as indicated by the C region trees, TCR-γ and -α/δ sequences represent the earliest immune receptor sequences, and since the TCR-α/δ sequences are highly convergent (Guglielmi at al. 1988<$REFLINK> ), this suggests that γδ-TCRs represent the oldest of the extant immune receptors, whereas αβ-TCRs and Ig’s are more recent. The V region trees, which may reflect functional relationships among antigen-binding sequences, indicate that the IgH’s are the earliest, implying that the earliest immune receptors had Ig-like antigen-binding properties. Despite these important distinctions in V and C region tree topologies, both data sets indicate that direct antigen-recognition is the primitive condition for T and B cells, while indirect antigen recognition is clearly a derived characteristic of αβ T cells. Indeed, it would be startling if this were not the case, since direct ligand recognition and binding is typical of other cell surface adhesion molecules in the IGSF (Hunkapiller and Hood 1989<$REFLINK> ).

Even before the antigen receptors of T cells were identified and long before it was realized that two classes of T cells existed, it was suggested that the ancestral immune cells were probably T cells (Marchalonis 1977<$REFLINK> ). Since then, several authors (Marchalonis and Schluter 1990<$REFLINK> ; Stewart 1992<$REFLINK> ; Thompson 1995<$REFLINK> ; Flajnik 1998<$REFLINK> ) have suggested that primordial immune cells were γδ-like. Many characteristics of γδ T cells fit the role of a first line of immune defense (Brenner, Strominger, and Krangel 1988<$REFLINK> ; Raulet 1989<$REFLINK> ; Allison and Havran 1991<$REFLINK> ). The primordial immune cell would have to have been capable of self- versus nonself-recognition, as well as of protective effector functions such as cytolysis, and its receptor would likely have been membrane-bound rather than secreted. Both of these capabilities are present in γδ T cells. It might also be expected that the primordial immune cells would have been capable of recognizing a wide array of antigens, many of which could be ingested through the gut. Despite their limited V-gene repertoire in comparison with αβ V genes, γδ TCRs display far greater N region (a region of random nucleotide addition between the V and D regions) variability and are potentially capable of recognizing far more antigens than are αβ T cells. γδ T cells can respond to antigen challenge very quickly, without a requirement for professional APCs, and have been suggested as a first line of immune defense in mammalian immune systems. On the other hand, αβ T cells require professional APC and antigen presentation by MHC molecules and exhibit highly specialized interactions with B cells in the maturation of antibody responses. A characteristic of most receptors in the IGSF is that they are membrane-bound, an ancestral characteristic that is retained in γδ and αβ T cells. However, B cells can secrete immunoglobulins, a trait that is likely derived (Flajnik 1998<$REFLINK> ). In fact, while many characteristics of γδ T cells are suggestive of an ancient and primary role in immune defense, the roles of αβ T cells and B cells imply highly specialized, and probably derived, roles in the immune system, supporting the conclusion that γδ-like T cells gave rise to αβ T cells and B cells. Several lines of structural evidence also support the suggestion that γδ T cells exhibit ancestral antigen-binding characteristics. γδ-TCRs recognize a wider variety of antigens than do αβ-TCRs, including both peptides and phospholigands, and they are especially reactive to conserved molecules such as mycobacterial and heat shock proteins (Constant et al. 1994<$REFLINK> ; Burk, Mori, and DeLibero 1995<$REFLINK> ). An analysis of the length distributions of V region chains, those that actually interact with antigen, in T and B cells indicates greater similarity of γδ-TCRs to Ig’s than to αβ-TCRs (Rock et al. 1994<$REFLINK> ), and this also strongly suggests that γδ-TCRs bind directly to antigens, as do immunoglobulins. αβ T cells display two cell surface markers, CD4 and CD8, that are required for binding of the αβ-TCR to MHC molecules, but most γδ T cells carry neither of these; thus, they are not MHC-restricted, and no alternative presentation molecules have been discovered. The differences in antigen binding are probably reflected in the three-dimensional structures of the two types of TCR. Although the structure of γδ-TCR molecules has not been elucidated, the structure of an αβ-TCR bound to a peptide+MHC complex has been studied (Garcia et al. 1996<$REFLINK> ). The antigen-binding surface of the αβ-TCR is fairly flat and similar to the undulating surfaces where Ig’s bind to proteins, but it is also much smaller than the Ig-binding areas, perhaps bestowing greater complementarity of fit with antigen (Garcia et al. 1996<$REFLINK> ). It may therefore be predicted that the binding surface of the γδ-TCR is also relatively flat but more like an Ig than an αβ-TCR in size, as would be expected for recognition of a wider variety of antigens.

## Conclusions

From a phylogenetic point of view, it is a little surprising that so much effort has been devoted to finding the antigen-presenting molecule for γδ T cells (perusal of immunological literature generated in the last 15 years indicates the huge amount of effort and resources put into this search). It seems that this search has been based on an assumption that current functions of immune receptors, especially in humans (and perhaps mice; Flajnik 1998<$REFLINK> ) must reflect the evolutionary histories of the genes that encode them, as summarized in figure 1 . And yet this study and others, based on comparisons of TCR, IgL, and IgH genes across vertebrate taxa, clearly support an alternative evolutionary history of these sequences, as summarized in figure 6 . The major stumbling block to this scheme is the apparent antiquity of the IgH V genes, since it is difficult to imagine what their early function was if the primordial antigen receptor resembled a γδ T cell and did not use the IgH V genes for its receptors. On the other hand, the relatively “young age” of the IgH C genes supports the hypothesis that secretion of immune receptors (i.e., modern Ig’s) is a derived characteristic. It is now clear that indirect, MHC-restricted antigen recognition must be a derived characteristic of αβ T cells, while γδ T cells exhibit the primitive condition of direct antigen binding. This removes a major stumbling block in several schemes for the evolution of vertebrate immune systems (e.g., Stewart 1992<$REFLINK> ), because it shows that the origins of diverse T-cell characteristics, such as antigen recognition and effector functions, may be separate, independent evolutionary events. Continued phylogenetic analysis of TCR and Ig relationships, especially as more suitable outgroup sequences are obtained from organisms such as cyclostomes or sharks, should continue to improve our understanding of immune system evolution.

Dan Graur, Reviewing Editor

1

Keywords: T-cell receptor, immunoglobulin, immune system evolution.

2

Address for correspondence and reprints: M. H. Richards, Department of Biological Sciences, Brock University, St. Catharines, Ontario L2S 3A1, Canada. E-mail: mrichard@spartan.ac.brocku.ca.

Fig. 1.—Null hypotheses for phylogenetic relationships among members of the immunoglobulin superfamily (IGSF). It is generally assumed that the T-cell receptors (TCRs) and immunoglobulins (Ig’s) each comprise a monophyletic group and that these two groups are sister groups within the IGSF. Other IGSF molecules bind directly to their ligands, so direct binding must be an ancestral characteristic of the entire group. “Antigen binding” per se arose when the molecules evolved to be part of an immune system. A, Based on this hypothesis of relationship, immunologists often assume that the antigen-binding processes of αβ- and γδ-TCRS are similar: since αβ- TCRs recognize antigen indirectly, so too must γδ-TCRs. B, A second possibility, even with this pattern of relationship, is that indirect antigen binding is a derived characteristic of αβ-TCRs.

Fig. 1.—Null hypotheses for phylogenetic relationships among members of the immunoglobulin superfamily (IGSF). It is generally assumed that the T-cell receptors (TCRs) and immunoglobulins (Ig’s) each comprise a monophyletic group and that these two groups are sister groups within the IGSF. Other IGSF molecules bind directly to their ligands, so direct binding must be an ancestral characteristic of the entire group. “Antigen binding” per se arose when the molecules evolved to be part of an immune system. A, Based on this hypothesis of relationship, immunologists often assume that the antigen-binding processes of αβ- and γδ-TCRS are similar: since αβ- TCRs recognize antigen indirectly, so too must γδ-TCRs. B, A second possibility, even with this pattern of relationship, is that indirect antigen binding is a derived characteristic of αβ-TCRs.

Fig. 2.—Nelson parsimony consensus tree for the constant region genes. The original set of 73 equally parsimonious trees was reduced using successive approximations character weighting to 2 trees of length 2,345 with a consistency index of 0.566. Asterisks indicate the chondrichthyan group I (*), II (**), and III (***) immunoglobulin light-chain genes.

Fig. 2.—Nelson parsimony consensus tree for the constant region genes. The original set of 73 equally parsimonious trees was reduced using successive approximations character weighting to 2 trees of length 2,345 with a consistency index of 0.566. Asterisks indicate the chondrichthyan group I (*), II (**), and III (***) immunoglobulin light-chain genes.

Fig. 3.—Bootstrap neighbor-joining tree for the constant region genes. Numbers indicate bootstrap values over 1,000 trees. Asterisks indicate the chondrichthyan group I (*), II (**), and III (***) immunoglobulin light-chain genes.

Fig. 3.—Bootstrap neighbor-joining tree for the constant region genes. Numbers indicate bootstrap values over 1,000 trees. Asterisks indicate the chondrichthyan group I (*), II (**), and III (***) immunoglobulin light-chain genes.

Fig. 4.—Nelson parsimony consensus tree for the variable region genes. The original set of 119 trees was reduced using successive approximations character weighting to 2 trees of length 896 with a consistency index of 0.535. Asterisks indicate the chondrichthyan group I (*), II (**), and III (***) immunoglobulin light-chain genes.

Fig. 4.—Nelson parsimony consensus tree for the variable region genes. The original set of 119 trees was reduced using successive approximations character weighting to 2 trees of length 896 with a consistency index of 0.535. Asterisks indicate the chondrichthyan group I (*), II (**), and III (***) immunoglobulin light-chain genes.

Fig. 5.—Bootstrap neighbor-joining tree for the variable region genes. Numbers indicate bootstrap values over 1,000 trees. Asterisks indicate the chondrichthyan group I (*), II (**), and III (***) immunoglobulin light-chain genes.

Fig. 5.—Bootstrap neighbor-joining tree for the variable region genes. Numbers indicate bootstrap values over 1,000 trees. Asterisks indicate the chondrichthyan group I (*), II (**), and III (***) immunoglobulin light-chain genes.

Fig. 6.—The evolution of vertebrate immune receptors: evolutionary relationships based on the results of this and other studies (see text for references).

Fig. 6.—The evolution of vertebrate immune receptors: evolutionary relationships based on the results of this and other studies (see text for references).

Table 1 Sequences Used in Phylogenetic Analyses

We thank Laurence Packer and several anonymous reviewers for their very helpful comments and Dan Graur for his patience. This work was supported by NIH operating grants AR-39282 and AI-38583 to J.L.N. We also thank Michael Parker of the FHCRC Biocomputing Shared Resource Center and acknowledge its support from National Cancer Institute award NCI P30 CA15704.

## References

Allison, J. P., and W. L. Havran.
1991
. The immunobiology of T cells with invariant γδ antigen receptors.
Annu. Rev. Immunol.

9
:
679
–705.
Bartl, S., D. Baltimore, and I. L. Weissman.
1994
. Molecular evolution of the vertebrate immune system. Proc. Natl. Acad. Sci. USA 91:10769–10770.
Bernstein, R. M., S. F. Schluter, S. Shen, and J. J. Marchalonis.
1996
. A new high molecular weight immunoglobulin class from the carcharhine shark: implications for the properties of the primordial immunoglobulin. Proc. Natl. Acad. Sci. USA 93:3289–3293.
Brenner, M. B., J. L. Strominger, and M. S. Krangel.
1988
. The γδ T cell receptor.

43
:
133
–192.
Burk, M. R., L. Mori, and G. DeLibero.
1995
. Human Vγ9-Vδ2 cells are stimulated in a cross-reactive fashion by a variety of phosphorylated metabolites.
Eur. J. Immunol.

25
:
2052
–2058.
Carpenter, J. M.
1988
Chien, Y.-H., and R. Jores.
1995
. T cells with B-cell-like properties.
Curr. Biol.

5
:
1116
–1118.
Chretien, I., A. Marcuz, J. Fellah, J. Charlemagne, and L. Du Pasquier.
1997
. The T cell receptor beta genes of Xenopus. Eur.
J. Immunol.

27
:
763
–771.
Constant, P., F. Davodeau, M.-A. Peyrat, Y. Poquet, G. Puzo, M. Bonneville, and J.-J. Fournie.
1994
. Stimulation of human γδ T cells by nonpeptidic mycobacterial ligands. Science 264:267–270.
Devereaux, J., P. Haeberli, and O. Smithies.
1984
. A comprehensive set of sequence analysis programs for the VAX.
Nucleic Acids Res.

12
:
387
–395.
Fellah, J. S., M. V. Wiles, J. Charlemagne, and J. Schwager.
1992
. Evolution of vertebrate IgM: complete amino acid sequences of the constant region of Ambystoma mexicanum μ chain deduced from cDNA sequences.
Eur. J. Immunol.

22
:
2595
–2601.
Felsenstein, J.
1993
. PHYLIP (phylogeny inference package). Version 3.57c. Distributed by the author, Department of Genetics, University of Washington, Seattle.
Feng, D.-F., and R. F. Doolittle.
1987
. Progressive sequence alignment as a prerequisite to correct phylogenetic trees.
J. Mol. Evol.

25
:
351
–360.
Flajnik, M. F.
1998
. Churchill and the immune system of ectothermic vertebrates.
Immunol. Rev.

166
:
5
–14.
Garcia, K. C., M. Degan, R. L. Stanfield, A. Brunmark, M. R. Jackson, P. A. Peterson, L. Teyton, and I. A. Wilson.
1996
. An αβ T cell receptor structure at 2.5 Å and its orientation in the TCR-MHC complex. Science 274:209–219.
Greenberg, A. S., D. Avila, M. Hughes, A. L. Hughes, E. McKinney, and M. F. Flajnik.
1995
. A novel antigen receptor gene family that undergoes rearrangement and extensive somatic diversification in sharks. Nature 374:168–173.
Greenberg, A. S., A. L. Hughes, J. Guo, D. Avila, E. C. McKinney, and M. F. Flajnik.
1996
. A novel “chimeric” antibody class in cartilaginous fish: IgM may not be the primordial immunoglobulin.
Eur. J. Immunol.

26
:
1123
–1129.
Guglielmi, P. E., F. Davi, L. D’Auriol, J. Bories, J. Dausset, and A. Bensussan.
1988
. Use of a variable α region to create a functional T-cell receptor δ chain. Proc. Natl. Acad. Sci. USA 85:5634–5638.
Haas, W., P. Pereira, and S. Tonegawa.
1993
. Gamma/delta cells.
Annu. Rev. Immunol.

11
:
637
–685.
Hamers-Casterman, C., T. Atarhouch, S. Muyldermans, G. Robinson, C. Hamers, E. Bajyana Songa, N. Bendahman, and R. Hamers.
1993
. Naturally occurring antibodies devoid of light chains. Nature 363:446–448.
Higgins, D. G., and P. M. Sharp.
1988
. CLUSTAL: a package for performing multiple sequence alignments on a microcomputer. Gene 73:237–244.
Hunkapiller, T., and L. Hood.
1989
. Diversity of the immunoglobulin gene superfamily.

44
:
1
–63.
Kaufmann, S. H. E.
1996
. γ/δ and other unconventional T lymphocytes: what do they see and what do they do? Proc. Natl. Acad. Sci. USA 93:2272–2279.
Klein, J.
1998
. In an immunological twilight zone. Proc. Natl. Acad. Sci. USA 95:11504–11505.
Marchalonis, J. J.
1977
. Immunity in evolution. Harvard University Press, Cambridge, Mass.
Marchalonis, J. J., and S. F. Schluter.
1990
. Origins of immunoglobulins and immune recognition molecules. Bioscience 40:758–768.
Marchalonis, J. J., S. F. Schluter, R. M. Bernstein, and V. S. Hohman.
1998
. Antibodies of sharks: revolution and evolution.
Immunol. Rev.

166
:
103
–122.
Marchalonis, J. J., S. F. Schluter, and A. B. Edmundson.
1997
. The T-cell receptor as immunoglobulin: paradigm regained.
Proc. Soc. Exp. Biol. Med.

216
:
303
–318.
Rast, J. P., M. K. Anderson, T. Ota, R. T. Litman, M. Margittai, J. J. Shamblott, and G. W. Litman.
1994
. Immunoglobulin light chain class multiplicity and alternative organizational forms in early vertebrate phylogeny. Immunogenetics 40:83–99.
Rast, J. P., M. K. Anderson, S. J. Strong, C. Luer, R. T. Litman, and G. W. Litman.
1997
. α, β, γ, and δ T cell antigen receptor genes arose early in vertebrate phylogeny. Immunity 6:1–11.
Rast, J. P., and G. W. Litman.
1998
. Towards understanding the evolutionary origins and early diversification of rearranging antigen receptors.
Immunol. Rev.

166
:
79
–86.
Raulet, D. H.
1989
. The structure, function, and molecular genetics of the γ/δ T cell receptor.
Annu. Rev. Immunol.

7
:
175
–207.
Rock, E. P., P. R. Sibbald, M. M. Davis, and Y.-H. Chien.
1994
. CDR3 length in antigen-specific immune receptors.
J. Exp. Med.

179
:
323
–328.
Saitou, N., and M. Nei.
1987
. The neighbor-joining method: a new method for reconstructing phylogenetic trees.
Mol. Biol. Evol.

4
:
406
–425.
Schatz, D. G., M. A. Oettinger, and M. S. Schlissel.
1992
. V(D)J recombination: molecular biology and regulation.
Annu. Rev. Immunol.

10
:
359
–383.
Schild, H., N. Macaddat, C. Litzenberger, E. W. Erich, M. M. Davis, J. A. Bluestone, L. Matis, R. K. Draper, and Y. Chien.
1994
. The nature of major histocompatibility complex recognition by γδ T-cells. Cell 76:29–37.
Schluter, S. F., R. M. Bernstein, and J. J. Marchalonis.
1997
. Molecular origins and evolution of immunoglobulin heavy-chain genes of jawed vertebrates. Immunol. Today 18:543–549.
Shen, S. X., R. M. Bernstein, S. F. Schluter, and J. J. Marchalonis.
1996
. Heavy-chain variable regions in carcharhine sharks: development of a comprehensive model for the evolution of VH domains among the gnathanstomes.
Immunol. Cell Biol.

74
:
357
–364.
Stewart, J.
1992
. Immunoglobulins did not arise in evolution to fight infection. Immunol. Today 13:396–399.
Swofford, D. L.
1993
. PAUP: phylogenetic analysis using parsimony. Version 3.1. Illinois Natural History Survey, Champaign.
Takihara, Y., J. Reimann, E. Michalopoulos, E. Ciccone, M. Lorenzo, and T. W. Mak.
1989
. Diversity and structure of human T cell receptor δ chain genes in peripheral blood γ/δ-bearing T lymphocytes.
J. Exp. Med.

169
:
393
–405.
Thompson, C. B.
1995
. New insights into V(D)J recombination and its role in the evolution of the immune system. Immunity 3:531–539.
Williams, A. F.
1987
. A year in the life of the immunoglobulin superfamily. Immunol. Today 8:298–303.