Dual Nature of the Adaptive Immune System in Lampreys
Abstract
Jawless vertebrates use variable lymphocyte receptors (VLR) comprised of leucine-rich-repeat (LRR) segments as counterparts of the immunoglobulin based receptors that jawed vertebrates use for antigen recognition. Highly diverse VLR genes are somatically assembled by the insertion of variable LRR sequences into incomplete germline VLRA and VLRB genes. Here we show that VLRA and VLRB anticipatory receptors are expressed by separate lymphocyte populations through monoallelic VLRA or VLRB assembly in concert with expression of Cytosine deaminase 1 or Cytosine deaminase 2, respectively. Distinctive gene expression profiles for VLRA+ and VLRB+ lymphocytes resemble those of mammalian T and B cells. Although both VLRA and VLRB cells proliferate in response to antigenic stimulation, only the VLRB lymphocytes bind native antigens and differentiate into VLR antibody secreting cells. Conversely, VLRA lymphocytes respond preferentially to a classical T cell mitogen and upregulate their expression of proinflammatory cytokine genes, IL-17 and MIF. The finding of T-like and B-like lymphocytes in lampreys offers new insight into the evolution of adaptive immunity.
The ability to mount specific immune responses to a virtually unlimited variety of antigens is apparently unique to vertebrates. The adaptive immune system in jawed vertebrates is characterized by antigen specific cellular and humoral responses, which together provide durable protective immunity1. This division of labor is conducted by developmentally separate, but functionally intertwined lineages of clonally diverse lymphocytes, named T and B cells because they are generated in the thymus or in the avian bursa of Fabricius and mammalian bone marrow2, 3. For antigen recognition, both populations use the same type of immunoglobulin (Ig)-based receptors. The T cell receptors (TCR) and B cell receptors (BCR) are assembled during lymphocyte differentiation by the recombination of different variable (V), diversity (D) and joining (J) immunoglobulin (Ig) gene segments, imprecise V(D)J splicing, and insertion of non-templated nucleotides at the junctions4–7. The assembly process is usually completed only on one allele, so that each T and B cell expresses a different antigen recognition receptor as the basis for clonal diversity8. Although the TCR and BCR genes are remarkably similar, the antigen binding repertoires of T and B cells are very different, primarily because immature T cells undergo clonal selection on the basis of how well they recognize MHC class I or class II molecules and their self peptide cargos 9–12. This intrathymic selection serves to eliminate T cells that can recognize and damage self tissues. The ensuing collaboration between T and B cells to achieve effective cellular and humoral immunity involves direct cell contact and indirect communication via cytokines and their receptors1.
When and how this complex adaptive immune system arose are subjects of much conjecture. Even the most basal jawed vertebrates, the cartilaginous sharks, skates and rays, have a thymus in which T cell progenitors assemble their TCR genes and hematopoietic tissues wherein B cell progenitors assemble their BCR genes13,14. Although the extant jawless vertebrates, lampreys and hagfish, were noted earlier to have lymphocyte-like cells and to produce antigen-specific agglutinins after immunization, neither a recognizable thymus nor TCR, Ig and MHC genes have been found15–18. More recent studies indicate that these agnathans generate clonally diverse anticipatory receptors via the assembly of genes for variable lymphocyte receptors (VLR) comprised of leucine-rich-repeat (LRR) subunits and an invariant membrane-proximal stalk region 19. The two germline VLR genes, VLRA and VLRB, are incomplete in that they have coding sequences only for portions of the N-terminal LRR and C-terminal LRR subunits and for the stalk region19–23. However, each germline VLR gene is flanked by hundreds of LRR encoding sequences and these are randomly used as templates to add the missing LRR segments needed for completion of a mature VLR gene22. A gene conversion mechanism has been postulated for VLR assembly21, 23, 24, and candidate AID-APOBEC deaminase participants have been identified in the lamprey 21. VLRB assembly has been shown to be confined to one allele19, 23, ensuring that each VLRB lymphocyte expresses a receptor of unique sequence, but it is not yet known whether the same is true for VLRA.
In earlier studies, we found that the VLRB bearing lymphocytes resemble B lymphocytes of jawed vertebrates. They can bind bacterial or erythrocyte antigens and respond by proliferation and differentiation into plasmacytes that secrete multimeric VLRB antibodies specific for either protein or carbohydrate epitopes25, 26. In this study, we used VLRA specific antibodies to identify a second population of lamprey lymphocytes and to determine their VLR gene assembly status, gene expression profile, tissue distribution, morphology and functional responses to antigens and mitogens in comparison with the VLRB lymphocytes. Our findings indicate that VLRA lymphocytes are surprisingly similar to the T cells in jawed vertebrates, and suggest functional interactions between the VLRA and VLRB lymphocyte populations.
Identification of VLRA+ and VLRB+ lymphocyte populations
The invariant stalk regions of the lamprey VLRA and VLRB proteins share less than 11% amino acid sequence identity, a difference that allowed the production of VLRB specific monoclonal antibodies for use in characterizing the VLRB bearing lymphocytes22, 25. In order to identify VLRA bearing lymphocytes, we used a VLRA-Ig Fc fusion protein as an immunogen to produce two rabbit antisera and three mouse monoclonal antibodies. The VLRA specificity of the antibody preparations was established by comparative analysis with a VLRB monoclonal antibody in ELISA and immunoblot assays (Supplementary Fig. 1). Using these VLR specific antibodies for immunofluorescence staining assays, we could identify discrete VLRA+ and VLRB+ lymphocyte populations in the principal lymphoid tissues of lamprey larvae, the blood, kidneys, typhlosole and gill region (Fig. 1a). VLRA bearing lymphocytes were detected in the same relative frequencies by each of the VLRA specific reagents (Supplementary Fig. 2). Notably, the VLRB+ lymphocytes significantly outnumbered VLRA+ lymphocytes in all tissues, except the gill region. The VLRB predominance was most striking in blood and kidneys, wherein the VLRB to VLRA ratio was approximately 8 to 1 (Fig. 1b); examination of adult blood samples indicated persistent VLRB predominance. These findings indicate that VLRA- and VLRB-bearing cells belong to separate lymphocyte populations that are dispersed through the body in a characteristic tissue distribution pattern.
Genetic basis of VLRA versus VLRB expression
To explore the genetic basis for the differential expression of VLR types by the two lymphocyte populations, we used fluorescence activated cell sorting to purify VLRA+, VLRB+ and VLRA/VLRB double-negative (DN) populations of cells having the light scatter characteristics of lymphocytes. The VLRA+ and VLRB+ lymphocytes could be separated cleanly (>95% purity), whereas the DN population included several cell types, among which nucleated thrombocytes constituted the vast majority. When VLR transcription was evaluated for each population by quantitative PCR analysis, the VLRA+ lymphocytes expressed VLRA transcripts exclusively, VLRB+ lymphocytes expressed VLRB transcripts exclusively, and the DN population expressed few or no VLR transcripts of either type. A representative experiment illustrating these results for adult blood cells is shown in Fig. 1c, and the same pattern of selective VLRA and VLRB expression was observed for the sorted populations of blood, kidney and typhlosole cells from larval stage lampreys three to four years of age (Supplementary Fig. 3).
When VLRA and VLRB gene configurations were examined for the VLRA+, VLRB+ and DN cells using paired 5’ and 3’ PCR primers (Fig. 2a and b), we observed mature VLRA assembly only in the VLRA+ lymphocyte population and the same relationship held for VLRB assembly by VLRB+ lymphocytes. Mature and germline VLRA PCR products were observed in equal proportions (Supplementary Fig. 4), indicating that VLRA assembly occurs on one allele, a finding that parallels the monoallelic pattern of VLRB assembly demonstrated previously19, 23. Because the lamprey AID-like deaminases have been postulated to serve a role in the gene conversion mechanism postulated for VLR gene assembly21, we examined the expression of CDA1 and CDA2 in the sorted cell populations. Remarkably, CDA1 expression was observed only in VLRA+ cells, while CDA2 expression was limited to the VLRB+ population (Fig. 2c). These findings suggest that the monoallelic VLRA versus VLRB gene assembly in the two lymphocyte populations is linked to expression of the genes for the CDA1 or CDA2 deaminases, respectively.
Distinctive gene expression profiles
Examination of the VLRA+, VLRB+ and DN populations for expression of a selected panel of genes indicated they have different gene expression profiles, and the distinctions were particularly informative for the two lymphocyte populations. The discriminating genes expressed by VLRA and VLRB lymphocyte populations included ones for transcription factors, Notch1, cytokine/chemokines and their receptors, Toll-like receptors (TLR) and various signaling molecules (Fig. 3, Supplementary Tables 1 – 3 and Supplementary Fig. 5). A surprising number of the genes that were selectively expressed by VLRA lymphocytes are orthologs of genes that are typically, although not exclusively, expressed by T lineage cells in jawed vertebrates. The genes expressed preferentially by VLRA lymphocytes included several transcription factors that may be used for T cell differentiation, GATA2/327, 28, c-Rel29, aryl hydrocarbon receptor (AHR)30 and BCL11b31, a chemokine receptor involved in lymphocyte progenitor homing to the thymus, CCR932, the T cell fate determining molecule, Notch133, a receptor protein tyrosine phosphatase that is essential for T cell development, CD4534, the receptor for interleukin 8 (IL-8), CXCR235, and two pro-inflammatory cytokines made by T cells, IL-1736, 37 and macrophage migration inhibitory factor (MIF)38, 39. The contrasting gene expression profile noted for the VLRB population included genes for the hematopoietic progenitor homing receptor, CXCR440, the herpes virus entry molecule (TNFR 14) which binds to LIGHT on T cells41, two components of the BCR mediated signaling cascade, Syk42and the B cell adaptor protein (BCAP)43, an Ig superfamily member with inhibitory potential, TCRL44, the chemotactic inflammatory cytokine, IL-845, the receptor for IL-1736, and the Toll-like receptor orthologs, TLR2abc, TLR7 and TLR10, the ligation of which is important for B cell activation46–48. This analysis, albeit limited by the incomplete status of the lamprey genomic sequence analysis, reveals gene expression profiles for VLRA and VLRB lymphocytes that are surprisingly similar to those for T and B cells in jawed vertebrates.
VLRA and VLRB lymphocyte responses to antigens
In response to immunization with Bacillus anthracis exosporium, antigen-binding VLRB+ members of the lymphocyte population undergo proliferation and differentiation into plasmacytes that secrete multivalent VLRB antibodies against BclA, a major protein in the outermost layer of B. anthracis spores25. In the present experiments, we found that both the VLRA and VLRB lymphocyte populations responded to this antigen with comparable levels of proliferation (Fig. 4a). In contrast, VLRA lymphocytes that could bind anthrax spores were not demonstrable either before or after immunization (Fig. 4b); the same paradoxical results were obtained when Escherichia coli binding by VLRA and VLRB lymphocyte was examined before and after immunization (data not shown). Equally surprising, the VLRA population failed to secrete VLRA proteins that could bind B. anthracis, E. coli, Salmonella typhimurium, or Streptococcus pneumoniae following immunizations that led to the secretion of specific VLRB antibodies for each of these bacteria (Fig. 4c and Supplementary Fig. 6). In addition, VLRA proteins could not be detected in lamprey plasma by immunoblot analysis (Fig. 4d and Supplementary Fig. 7).
To determine whether an inherent difference in the two types of receptors could explain the failure to secrete VLRA, we transfected cells of the 293T human embryonic kidney cell line with either VLRA or VLRB cDNAs to compare their patterns of intracellular synthesis and secretion. While the VLRB transfected cells secreted their VLRB protein products in the multimeric forms described previously26, the VLRA transfected cells expressed their VLRA protein products on the cell surface, but did not secrete them (Supplementary Fig. 8). These results indicate that VLRA lymphocytes respond to antigenic stimulation in a manner very reminiscent of the T cell response and confirm that VLRB lymphocytes respond to antigens much like B lymphocytes. In addition, the contrast between the abilities of VLRA and VLRB lymphocytes to bind native bacterial surface epitopes strongly suggests they have different receptor repertoires.
Preferential PHA responsiveness of VLRA lymphocytes
Because of the resemblance between the VLRA lymphocytes and T lymphocytes in jawed vertebrates, we compared the responses of VLRA and VLRB lymphocytes to stimulation by the plant lectin, phytohemagglutinin (PHA), a classical T cell mitogen. When larvae were injected with PHA, the VLRA+ lymphocytes responded vigorously by undergoing lymphoblastoid transformation and proliferation to become the predominant lymphocyte population (Fig. 5a and b). The activated VLRA+ cells were found by transmission EM to be large lymphoblasts with much less rough endoplasmic reticulum than the previously described VLRB secreting plasmacytes (Fig. 5c)25. While VLRB lymphocytes also responded to PHA stimulation, their response was much weaker in terms of cell size increase and numbers of proliferating cells in this population. Nevertheless, PHA stimulation enhanced VLRB cell expression of IL-8 transcripts, in contrast to its enhancement of VLRA cell expression of IL-17 and MIF transcripts (Fig. 5d). These findings, coupled with the reciprocal expression of IL-8R by VLRA lymphocytes and IL-17R expression by VLRB lymphocytes (Fig. 3), are indicative of the potential for functional interactions between the two lymphocyte populations.
Discussion
This unforeseen genotypic and functional division of lymphocyte differentiation in lampreys offers an intriguing new piece to the puzzle of how the adaptive immune systems may have evolved. Our initial definition of the VLRA and VLRB lymphocyte populations also offers interesting clues to presently unresolved issues regarding their differentiation and function.
One question that is immediately raised by our findings concerns the anatomical location(s) in which VLRA lymphocytes begin their differentiation. Indeed, one of the reasons why the remarkable resemblance of VLRA lymphocytes to T lymphocytes was so surprising is the thus far fruitless quest to identify a thymus equivalent in lampreys16, 17. However, the relatively high concentration of VLRA lymphocytes in the gill region could favor the branchial region as a potential site for their generation. In addition, our gene profile analysis indicates that VLRA lymphocytes express an array of transcription factors, chemokine receptors and Notch1 that lymphocyte progenitors use for homing to the thymus and T lineage commitment in jawed vertebrates. Analysis of the expression of the corresponding chemokines and Notch1 ligands in the pharyngeal pouch-derived epithelial cells and the neighboring VLRA lymphocytes could therefore facilitate the identification of a thymus equivalent in lampreys.
Solving the issue of where the VLRA and VLRB genes are assembled will also require better understanding of the molecular mechanisms involved in this process. The donor and recipient LRR sequences used in VLR assembly lack recombination signal sequences that characterize the Ig V, D and J gene segments, and the Rag1 and Rag2 recombination activation genes employed for TCR and BCR gene assembly have not been found in lampreys18. The presently known features of VLR gene assembly instead suggest that a gene conversion mechanism is used to copy the donor LRR sequences21, 23, 24. In accord with this hypothesis, two APOBEC-AID family orthologs, CDA1 and CDA2, have been identified in the lamprey. CDA1 has been shown to have an AID-like DNA deaminase mutator function in bacteria and yeast lacking uracil DNA glycosylase activity21.
Although the functional competency of CDA2 has not yet been elaborated, it shares the DNA deaminase catalytic motif with CDA1 and has in addition an extended C-terminus containing a small AT-hook peptide motif that may confer DNA binding specificity. These distinctive features and the restriction of CDA1 expression to the VLRA lymphocyte population and CDA2 expression to VLRB lymphocytes thus imply the potential for selective involvement of these enzymes in VLRA and VLRB assembly during lymphocyte development.
One of the most provocative issues raised by our findings concerns the mode of antigen recognition by VLRA bearing lymphocytes. We expected to find comparable antigen binding and response patterns for the two lymphocyte populations, since the VLRA and VLRB proteins are very similar and have comparable diversity. The proliferative responses of VLRA and VLRB lymphocytes to immunization indeed were found to be equivalent. Surprisingly though, antigen-binding VLRA lymphocytes could not be found before or after immunization, whereas specific antigen-binding VLRB lymphocytes were easily detected in naive animals and immunization increased their frequency. One possible explanation for the inability of VLRA lymphocytes to bind native epitopes on immunogens that induce their proliferation is that VLRA lymphocytes recognize processed antigens and undergo selection in a manner analogous to the T lymphocyte repertoire selection in jawed vertebrates. Lampreys lack the major histocompatibility complex genes for MHC class I and class II molecules that accessory cells use to present peptide fragments to T cells in jawed vertebrates18. However, a precedent for an alternative mechanism for self recognition is provided by the recent identification of a highly polymorphic fusion/histocompatibility gene locus in Botryllus schlosseri49, a colonial tunicate member of a vertebrate sister group in the new chordate phylogeny50. The nature of the antigens recognized by VLRA anticipatory receptors presently remains an unresolved enigma.
The failure of VLRA cells to secrete their receptors as soluble proteins was the most striking functional difference observed for the VLRA and VLRB lymphocytes. This feature could reflect an inability of VLRA cells to develop the necessary secretory machinery. However, when VLRA proteins were produced by a human cell line that can secrete VLRB proteins, the transfected cells produced cell surface VLRA proteins, but did not secrete them. The VLRA proteins, like the TCR proteins, thus appear to be expressed exclusively as transmembrane molecules. The preferential expression of TLR2, TLR7 and TLR10 orthologs by the VLRB lymphocytes suggests that TLR ligands may facilitate activation of this population of lymphocytes in a manner similar to their roles in activating B lymphocytes. The development of appropriate conditions for culturing lamprey lymphocytes and cell lines will facilitate the future analysis of the signaling mechanisms employed by VLRA and VLRB lymphocytes.
Many of the elaborate cytokine networks that T and B cells use to coordinate effective inflammatory responses in jawed vertebrates are not evident in the more basal jawless vertebrates. Agnathan orthologs have not been identified for the interferons and interferon receptors, IL-2 / IL-2R, IL-4 / IL4R, IL-7 / IL-7R and other cytokine/cytokine receptor pairs. Nevertheless, the VLRA lymphocytes preferentially express IL-17 and MIF transcripts and their enhanced expression after activation suggests that VLRA lymphocytes may use these pro-inflammatory cytokines to engage the IL-17R bearing VLRB lymphocytes for defense purposes, in a manner similar to their T cell counterparts. Conversely, activated VLRB lymphocytes upregulate their expression of IL-8 transcripts, suggesting they may use this cytokine to attract and engage IL-8R bearing VLRA lymphocytes and other cell types to amplify the response to potential pathogens.
The evolution around 500 million years ago of two very different anticipatory receptors of comparable diversity in jawless and jawed vertebrates, while conserving similar compartmentalization of lymphocyte differentiation, strongly attests the survival value of adaptive immunity. The evolution of an LRR-based adaptive immune system in lampreys and hagfish, apparently before the second round of whole genomic replication, and an Ig-based adaptive immune system in jawed vertebrates, after a second round of genomic duplication50, suggests that these very different types of anticipatory receptors were convergent evolutionary acquisitions19. The failure so far to find relics of the VLR genes in jawed vertebrates is also consistent with an independent evolution of clonally diverse LRR-based receptors versus Ig-based receptors for otherwise similar adaptive immune systems. An obvious advantage of adaptive immunity over the elaborate invertebrate strategies for innate immunity is the capacity for memory of previous encounters with infectious agents. On the other hand, the development of a randomly generated anticipatory repertoire of great diversity inevitably creates the hazard of self-reactivity and attendant autoimmunity. The orchestration of cellular and humoral immune responses by T cells, whose TCR repertoire is shaped during development to discriminate between self and non-self antigens, provides an important safeguard against autoimmunity in jawed vertebrates. The potential for autoimmunity theoretically would necessitate a similar solution for the adaptive immune system in jawless vertebrates. Our findings thus suggest that dual recognition and response arms with intertwined function were fundamental to the evolution of adaptive immunity in vertebrates.
Methods Summary
Sea lamprey larvae (8–15 cm, 2–4 years of age) were collected from tributaries to Lake Michigan (Lamprey services, Michigan) and housed in sand-lined aquariums at 20°C. Larvae anesthetized with 0.1 g/l MS222 (Sigma) were given intracoelomic injections of antigens or mitogens (25 µg phytohemaglutinin (PHA)-L (Sigma)) prepared in 60 µl of sterile 0.67X PBS. Injections were administered on days 0 and 14, and tissues collected from lampreys euthanized with 1 g/l MS222 on day 28 (unless otherwise indicated). Blood was collected in 30 mM EDTA/0.67X PBS. Buffy coat leukocytes were separated from red blood cells by centrifugation at 50×g. Leukocytes were isolated from kidney, typhlosole, and gills by disrupting tissues between frosted glass slides.
For VLRA antibody production, four VLRA cDNAs isolated from lamprey lymphocytes were cloned in-frame with the constant region of human IgG1 (IgG1-Fc). The VLRA-IgG1-Fc fusion proteins were expressed in HEK-293T cells and purified from tissue culture supernatants by protein A (GE Healthcare) chromatography. Monoclonal anti-VLRA antibodies were produced by immunization of BALB/c mice with VLRA-IgG1-Fc proteins emulsified in TiterMax Gold adjuvant. Lymphocytes from the draining lymph nodes of immunized mice were fused with the Ag8.653 myeloma cell line using PEG-1500 (Roche). Three VLRA-specific hybridoma clones (9A6 (IgG1), 9B3 (IgG1) and 2A5 (IgG1)) were identified by ELISA with VLRA recombinant protein and immunoblotting of VLRA transfectants. Anti-VLRA polyclonal antisera were produced by immunization of rabbits with the VLRA-IgG1-Fc proteins (PickCell Laboratories BV, Netherlands).
Methods
Antigens
Escherichia coli (ATCC 25922 strain), Salmonella enterica Typhimurium (WB99-SR11 strain), and Streptococcus pneumoniae (R36A strain) were provided by David E. Briles and William H. Benjamin, Jr. (University of Alabama at Birmingham). Bacteria were heat-killed by incubation at 60°C for 1 hour. 5×108 heat-killed E. coli, S. typhimurium, and 1×108S. pneumoniae were used for immunization. Bacillus anthracis (spores and exosporium) and Bacillus cereus spores were provided by Charles L. Turnbough Jr. (University of Alabama at Birmingham). Lampreys were immunized with 10 µg of B. anthracis exosporium.
Flow cytometric analysis of VLRA and VLRB lymphocytes
Leukocytes from blood and tissues were stained with anti-VLRA rabbit polyclonal serum (R110 or R111) or mouse mAbs (9A6, 9B3, 2A5, IgG1) and anti-VLRB mouse mAb (4C4, IgG2b) for 10 minutes on ice, washed twice and stained with fluorescent-conjugated secondary Abs. For anti-VLRA mAb (IgG1) and anti-VLRB mAb (4C4, IgG2b) staining combinations, R-phycoerythrin (PE)-conjugated goat anti-mouse IgG1 Abs (Southern Biotech) and FITC-conjugated goat anti-mouse IgG2b Abs (Southern biotech) were used. When staining with anti-VLRA rabbit serum and anti-VLRB mAb (4C4), RPE-conjugated goat anti-rabbit Ig Abs (Southern biotech) and allophycocyanin (APC)-conjugated goat anti-mouse Ig Abs (Southern Biotech) were used. Cells were analyzed on a CyAn ADP high-performance flow cytometer (Dako). Staining and washes were in 0.67X PBS with 1% BSA.
Immunoblotting
Samples were separated on 10% SDS-PAGE gels with or without 5% (vol/vol) 2-mercaptoethanol before transfer onto nitrocellulose membranes. Membranes were blocked overnight with 5% milk and incubated with anti-VLRA or anti-VLRB antibodies for one hour. After 5 washes with PBS-0.5% Tween-20, membranes were incubated with HRP-conjugated goat anti-rabbit or mouse Ig polyclonal Abs (Southern Biotech) and washed. Blots were developed using SuperSignal West Pico chemiluminescent substrate (Pierce, Rockford, IL).
Genomic DNA and PCR
Genomic DNA was extracted from VLRA+, VLRB+ and VLR double negative-sorted cells in blood and tissue lymphocytes using DNeasy kit (QIAGEN). Genomic PCR was carried out using primers VLRA-F + VLRA-R (Expand High Fidelity, Roche and Ex Taq, Takara) or VLRB-F + VLRB-R (Expand Long Template, Roche). Primers are listed in Supplementary Table 1.
Quantitative PCR
Target gene sequences were obtained from the National Center for Biotechnology Information database or the lamprey genome database of the Genome Sequencing Center at Washington University. Cells in the lymphocyte-gate were sorted into VLRA+, VLRB+ and VLR double negative populations on a BD FACS Aria II. RNA was extracted from each population using RNeasy kit with on-column DNA digestion by DNase I (QIAGEN). First-strand cDNA was synthesized with random hexamer primers and Superscript III (Invitrogen). Quantitative real-time PCR (qPCR) was done with SYBR Green on a 7900HT ABI Prism (Applied Biosystems). Three independent qPCR reactions were performed at different times for each sample. The value of the target gene was normalized to β-actin. The normalized value of each gene was compiled into a heat map (z-score x 10). Magenta: z > 0, green: z < 0 and yellow: z = 0 [z = (each value – average)/standard deviation]. Primers are described in Supplementary Table 1.
Proliferation Assay
Lampreys stimulated with antigen or mitogen were injected with 5 µg of 5-ethynyl-2’-deoxyuridine (EdU) (Invitrogen) in 60 µl 0.67X PBS and returned to their tanks for 24 hrs before collection of leukocytes for staining with anti-VLRA or anti-VLRB Abs. Stained cells were fixed in 4% paraformaldehyde for 15 min at room temperature, then resuspended in 100 µL of the 1X saponin-based permeabilization and wash buffer for 10 min on ice. EdU+ cell detection was performed according to the manufacturer’s protocol (Click-iT™ EdU Alexa Fluor 488 Flow Cytometry Assay Kit, Invitrogen). Briefly, cells were incubated with the EdU reaction cocktail, which contains Alexa Fluor 488-azide for 30 min at room temperature, washed twice in 1X saponin-based permeabilization and wash buffer, and analyzed by flow cytometry.
ELISA
for animals immunized with B. anthracis exosporium, 106B. anthracis spores or B. cereus strain T spores per well were applied on poly-L-lysine-coated plates as per Alder MN et al22. For animals immunized with bacteria, 5×106 heat-killed bacteria per well were applied on poly-L-lysine-coated plates. Plates were blocked with 1% BSA in PBS for 1 hr. Serially-diluted plasma was applied for 2 hours at 37°C. VLRs were detected with anti-VLRB (4C4) or anti-VLRA (9A6) mAb and alkaline phosphatase-conjugated goat anti-mouse Ig Abs (Southern Biotech). ELISAs were developed with phosphatase substrate (Sigma) and read at 405 nm (Versamax microplate reader; Molecular Devices, Sunnyvale, CA).
Lymphocyte antigen-binding assay
106 lamprey leukocytes were incubated with 106 Alexa Fluor 488-labeled spores or 107 Alexa Fluor 488-labeled bacteria for 20 min on ice. Spores and bacteria were labeled using the Alexa Fluor 488 Protein Labeling Kit (Invitrogen). Cells were stained with anti-VLRA (R110) and anti-VLRB (4C4) Abs followed by APC-conjugated goat anti-mouse Ig Abs (Southern Biotech) and R-PE-conjugated goat anti-rabbit Ig Abs (Southern Biotech), then analyzed by flow cytometry.
VLR expression in HEK-293T cells
VLRA and VLRB plasmids were transfected into HEK-293T cells cultured in DMEM containing 5% fetal bovine serum (FBS) using linear polyethylenimine (PEI), MW 25,000 (Polysciences, Inc.) at a 3:1 PEI to DNA ratio. Cells were separated from supernatants 48 h post-transfection by centrifugation at 300×g and lysed in 1% NP-40 lysis buffer.
Lymphocyte counts
Lymphocytes were extracted from lamprey kidney and typhlosole by passage through a 70 µm cell strainer (BD Bioscience) to obtain a single cell suspension. Cells were washed and resuspended in 1 ml 0.67X PBS. Total cells were counted in 25 µl of each sample on an Accuri C6 flow cytometer (Accuri Cytometers Inc, Michigan). The total number of cells in the ‘lymphocyte gate’ was calculated by the following formula: total lymphocyte number = number of cells in the lymphocytes gate × 40 (dilution factor: 1000 µl/25 µl). Lymphocytes were then stained with anti-VLRA and VLRB Abs to determine the percentage of VLRA+ and VLRB+ cells in the samples. These percentages were used to calculate the total number of lymphocytes of each type.
Electron microscopy
Blood lymphocytes from naïve or PHA stimulated lamprey were sorted on BD FACS Aria II (BD Bioscience). The VLRA-positive cells were prepared for transmission electron microscopic analysis as described as per Alder MN et al22.
Statistical Analysis
A two-sample Student’s t-test was used for statistical analysis.
VLRA and VLRB expression define distinct lymphocyte populations
a, FACS analysis of lymphocyte-gated cells stained with anti-VLRA (R110) and anti-VLRB (4C4) antibodies. b, Percentages of VLRA and VLRB expressing lymphocytes in lamprey tissues. Larvae (left panel) n = 23, except gill n = 6; adults (right panel) n = 8 (***P < 0.001; **P < 0.01). c, Lymphocyte-gated cells were stained as in (a) and separated based on surface expression of VLRA and VLRB by fluorescence-activated cell sorting (FACS). Sorted cells were analyzed for VLRA (left) and VLRB (right) transcripts by quantitative real-time PCR (qPCR), n = 3. Error bars indicate s.e.m.

Monoallelic assembly of VLRA and VLRB genes
a,b, Schematic of VLRA (a) and VLRB (b) genes before (top panel) and after (middle panel) gene assembly. Forward and reverse primer locations and predicted sizes of PCR products are indicated. Lymphocytes from the indicated tissues were stained with anti-VLRA (R110) and anti-VLR-B (4C4) antibodies and lymphocyte-gated cells were FACS sorted into three populations: VLRA−/VLRB− (DN), VLRA+ (A+), and VLRB+ (B+). VLRs were amplified from genomic DNA of the sorted lymphocyte populations (bottom panel). Germ-line (GL) and mature (M) products were verified by sequence analysis of representative DNA clones. c, CDA1 and CDA2 expression in sorted lymphocytes was measured by QPCR. Error bars indicate s.e.m., n = 3.

Differential gene expression profiles of VLRA+ and VLRB+ lymphocytes
Lamprey lymphocytes were FACS sorted into three populations on the basis of VLRA and VLRB surface expression. Relative transcript levels of the indicated genes were measured by QPCR and compiled into a heat map as described in Methods.

Antigen-activated VLRA+ lymphocytes do not secrete their receptors
a, Immunization with B. anthracis exosporium induces proliferation of VLRA+ and VLRB+ lymphocyte (n = 5; * P < 0.05). b, VLRA+ lymphocytes do not bind to B. anthracis spores before immunization or 14 days after booster immunization (n = 5; ***P < 0.001). c,d, VLRB, but not VLRA, is secreted. c, B. anthracis-specific VLRA and VLRB reactivity was evaluated by ELISA with anti-VLRA (9A6) and anti-VLRB (4C4) mAbs. d, Immunoblot of VLRA (lane 1) and VLRB (lane 2) HEK-293T cell transfectant lysates, naïve lamprey plasma (lane 3), and immunized plasma (lane 4) under reducing conditions. R110 = anti-VLRA rabbit antiserum; 9B3 = anti-VLRA mAb.

PHA preferentially stimulates VLRA+ lymphocytes
a,b,In vivo PHA stimulation. Lampreys were given 25 µg of PHA by intracoelomic injection on day 0. Proliferation (a) and absolute numbers (b) of VLRA+ and VLRB+ cells were determined on days 9 and 15 after injection; (a)n = 4; (b)n = 7; *P < 0.05; **P < 0.01; ***P < 0.001). c, Transmission electron microscopy imaging of unstimulated (left panel) and PHA stimulated (day 9) (right panel) VLRA+ cells. PHA stimulated cells have a larger diameter [naïve mean = 3.8±0.1µm (n = 7) versus PHA stimulated mean = 7.4±0.3µm (n = 7, P <0.01)] and increased cytoplasm. d, Cytokine expression by VLRA+ and VLRB+ lymphocytes. FACS sorted blood lymphocytes were isolated from PHA stimulated (day 7, n = 3) or unstimulated lampreys (n = 3) and analyzed for expression of IL-17 (left panel), MIF (middle panel), and IL-8 (right panel) by QPCR. Error bars indicate s.e.m.

Footnotes
Author Information
The authors declare no competing financial interests.
Acknowledgements
We thank Matthew N. Alder and Götz R. A. Ehrhardt for helpful suggestions and discussion; Charles L. Turnbough, Jr. for providing B. anthracis spores and exosporium; David E. Briles and William H. Benjamin, Jr. for E. coli, S. pneumoniae, and S. typhimurium; Hong Yi for help with electron microscopy; Sommer A. Durham and Robert E. Karaffa, II for help with cell sorting; Marsha Flurry for help with preparation of figures. This work is supported by National Institutes of Health and Georgia Research Alliance.
References
- 1. Fundamental Immunology2008Sixth EditionPhiladelphiaLippincott Williams & Wilkins
- 2. Delineation of the Thymic and Bursal Lymphoid Systems in the ChickenNature1965205143146[PubMed][Google Scholar]
- 3. The cellular basis of immunological responses. A synthesis of some current viewsLancet19692367371[PubMed][Google Scholar]
- 4. Somatic generation of antibody diversityNature1983302575581[PubMed][Google Scholar]
- 5. A human T cell-specific cDNA clone encodes a protein having extensive homology to immunoglobulin chainsNature1984308145149[PubMed][Google Scholar]
- 6. Isolation of cDNA clones encoding T cell-specific membrane-associated proteinsNature1984308149153[PubMed][Google Scholar]
- 7. Mechanism and control of V(D)J recombination at the immunoglobulin heavy chain locusAnnu. Rev. Immunol200624541570[PubMed][Google Scholar]
- 8. A modification of Jerne's theory of antibody production using the concept of clonal selectionAust. J. Sci1957206769[Google Scholar]
- 9. Cytotoxic thymus-derived lymphocytes in cerebrospinal fluid of mice with lymphocytic choriomeningitisJ. Exp. Med197313812661269[PubMed][Google Scholar]
- 10. Role of self-peptides in positively selecting the T-cell repertoireNature19903446567[PubMed][Google Scholar]
- 11. Selection of the T-cell repertoire: receptor-controlled checkpoints in T-cell developmentAdv. Immunol200484201238[PubMed][Google Scholar]
- 12. Perspective on antigen processing and presentationImmunol. Rev200218586102[PubMed][Google Scholar]
- 13. Comparative analyses of immunoglobulin genes: surprises and portentsNat. Rev. Immunol20022688698[PubMed][Google Scholar]
- 14. Reconstructing immune phylogeny: new perspectivesNat. Rev. Immunol20055866879[PubMed][Google Scholar]
- 15. The Evolution of the Immune Response. 3. Immunologic Responses in the LampreyJ. Exp. Med196412011511168[PubMed][Google Scholar]
- 16. The pharyngeal lymphoid tissue of lampreys. A morpho-functional equivalent of the vertebrate thymus?Thymus1988115965[PubMed][Google Scholar]
- 17. The evolutionary history of lymphoid organsNat. Immunol20078131135[PubMed][Google Scholar]
- 18. Lamprey lymphocyte-like cells express homologs of genes involved in immunologically relevant activities of mammalian lymphocytesProc. Natl. Acad. Sci. U S A2002991435614361[PubMed][Google Scholar]
- 19. Somatic diversification of variable lymphocyte receptors in the agnathan sea lampreyNature2004430174180[PubMed][Google Scholar]
- 20. Variable lymphocyte receptors in hagfishProc. Natl. Acad. Sci. U S A200510292249229[PubMed][Google Scholar]
- 21. Evolution and diversification of lamprey antigen receptors: evidence for involvement of an AID-APOBEC family cytosine deaminaseNat. Immunol20078647656[PubMed][Google Scholar]
- 22. Diversity and function of adaptive immune receptors in a jawless vertebrateScience200531019701973[PubMed][Google Scholar]
- 23. Antigen-receptor genes of the agnathan lamprey are assembled by a process involving copy choiceNat. Immunol20078206213[PubMed][Google Scholar]
- 24. The evolution of adaptive immune systemsCell2006124815822[PubMed][Google Scholar]
- 25. Antibody responses of variable lymphocyte receptors in the lampreyNat. Immunol20089319327[PubMed][Google Scholar]
- 26. Structure and specificity of lamprey monoclonal antibodiesProc. Natl. Acad. Sci. U S A200810520402045[PubMed][Google Scholar]
- 27. GATA3 and the T-cell lineage: essential functions before and after T-helper-2-cell differentiationNat. Rev. Immunol20099125135[PubMed][Google Scholar]
- 28. Transcriptome analysis of hagfish leukocytes: a framework for understanding the immune system of jawless fishesDev. Comp. Immunol2004289931003[PubMed][Google Scholar]
- 29. c-Rel-dependent priming of naive T cells by inflammatory cytokinesImmunity200523445458[PubMed][Google Scholar]
- 30. Control of T(reg) and T(H)17 cell differentiation by the aryl hydrocarbon receptorNature20084536571[PubMed][Google Scholar]
- 31. Bcl11b is required for differentiation and survival of alphabeta T lymphocytesNat. Immunol20034533539[PubMed][Google Scholar]
- 32. Coordination between CCR7-and CCR9-mediated chemokine signals in prevascular fetal thymus colonizationBlood200610825312539[PubMed][Google Scholar]
- 33. Regulation of lymphoid development, differentiation, and function by the Notch pathwayAnnu. Rev. Immunol200523945974[PubMed][Google Scholar]
- 34. Normal B lymphocyte development but impaired T cell maturation in CD45-exon6 protein tyrosine phosphatase-deficient miceCell199374143156[PubMed][Google Scholar]
- 35. Human T lymphocytes and mast cells differentially express and regulate extra-and intracellular CXCR1 and CXCR2Exp. Dermatol200413520525[PubMed][Google Scholar]
- 36. IL-17 family cytokines and the expanding diversity of effector T cell lineagesAnnu. Rev. Immunol200725821852[PubMed][Google Scholar]
- 37. Lamprey (Lethenteron japonicum) IL-17 upregulated by LPS-stimulation in the skin cellsImmunogenetics200759873882[PubMed][Google Scholar]
- 38. Requirement for T cells in the production of migration inhibitory factorJ. Exp. Med197514213061311[PubMed][Google Scholar]
- 39. Molecular cloning of a cDNA encoding a human macrophage migration inhibitory factorProc. Natl. Acad. Sci. U S A19898675227526[PubMed][Google Scholar]
- 40. The chemokine receptor CXCR4 is required for the retention of B lineage and granulocytic precursors within the bone marrow microenvironmentImmunity199910463471[PubMed][Google Scholar]
- 41. Balancing co-stimulation and inhibition with BTLA and HVEMNat. Rev. Immunol20066671681[PubMed][Google Scholar]
- 42. B cell signaling and tumorigenesisAnnu. Rev. Immunol200523415445[PubMed][Google Scholar]
- 43. BCAP: the tyrosine kinase substrate that connects B cell receptor to phosphoinositide 3-kinase activationImmunity200013817827[PubMed][Google Scholar]
- 44. Inhibitory signaling potential of a TCR-like molecule in lampreyEur. J. Immunol200939571579[PubMed][Google Scholar]
- 45. In vivo expression of interleukin-8, and regulated on activation, normal, T-cell expressed, and secreted, by human germinal centre B lymphocytesImmunology2003110296303[PubMed][Google Scholar]
- 46. Control of B-cell responses by Toll-like receptorsNature2005438364368[PubMed][Google Scholar]
- 47. Pathogen recognition and innate immunityCell2006124783801[PubMed][Google Scholar]
- 48. Lamprey TLRs with properties distinct from those of the variable lymphocyte receptorsJ. Immunol2007178397406[PubMed][Google Scholar]
- 49. Isolation and characterization of a protochordate histocompatibility locusNature2005438454459[PubMed][Google Scholar]
- 50. The amphioxus genome and the evolution of the chordate karyotypeNature200845310641071[PubMed][Google Scholar]