In recent years, the antigen recognition mechanism based on variable lymphocyte receptors (VLRs) was found in agnathan lamprey. To illuminate the genetic basis of immune response of lymphocyte-like cells in the mucosal immune system of lamprey and explore the evolutionary relationship of adaptive immune responses between the jawless and jawed vertebrates, we constructed cDNA libraries of lamprey (Lampetra japonica) gills before and after stimulation, and then performed high-throughput transcriptome sequencing and analysis. Through functional annotation of 88 525 assembled unigenes, 21 704 and 9769 unigenes were annotated in Gene Ontology (GO) and Kyto Encyclopedia of Genes and Genomes (KEGG) databases, respectively. Among 999 unigenes involved in multiple pathways of immune system, 184 unigenes were highly homologous to 51 TCR (T cell receptor) and BCR (B cell receptor) signalling molecules in higher vertebrates, indicating that molecules involved in adaptive immune signalling pathways in higher vertebrates also exist in lampreys. In addition, identification of five VLRA, seven VLRB and four VLRC molecules suggest that at least three types of lymphocyte subsets are distributed in lamprey gill mucosal immune tissues. The results of real-time fluorescence quantitative PCR showed that the expression levels of Lck, Fyn and Zap70 were up-regulated after immune stimulation while those of Syk, Btk and Blnk were not changed significantly, indicating the activation of TCR-like signal transduction pathway after antigen stimulation in lamprey gill tissues. Our studies preliminaryly proved that two parallel adaptive immune systems in jawless and jawed vertebrates have common genetic basis, and also provided valuable clues to the exploration of signalling processes of VLRA+, VLRB+, and VLRC+ lymphocyte-like cells in response to antigens.
[1] Holmgren J, Czerkinsky C. Mucosal immunity and vaccines. Nat Med , 2005, 11(4 Suppl): S45-S53.
[2] Lin SH, Davidson GA, Secombes CJ, Ellis AE. A morphological study of cells isolated from the perfused gill of dab and Atlantic salmon. J Fish Biol , 1998, 53(3): 560-568.
[3] Press CMcL, Evensen Ø. The morphology of the immune system in teleost fishes. Fish Shellfish Immunol , 1999, 9(4): 309-318.
[4] Lin SH, Ellis AE, Davidson GA, Secombes CJ. Migratory, respiratory burst and mitogenic responses of leucocytes isolated from the gills of rainbow trout ( Oncorhynchus mykiss ). Fish Shellfish Immunol , 1999, 9(3): 211-226.
[5] Forey PL, Janvier P. Agnathans and the origin of jawed vertebrates. Nature , 1993, 361(6408): 129-134.
[6] Braun CB. The sensory biology of the living jawless fishes: a phylogenetic assessment. Brain Behav Evol , 1996, 48(5): 262-276.
[7] Baker ME, Funder JW, Kattoula SR. Evolution of hormone selectivity in glucocorticoid and mineralocorticoid receptors. J Steroid Biochem Mol Biol , 2013, 137: 57-70.
[8] Cooper MD, Alder MN. The evolution of adaptive immune systems. Cell , 2006, 124(4): 815-822.
[9] Pancer Z, Amemiya CT, Ehrhardt GRA, Ceitlin J, Gartland GL, Cooper MD. Somatic diversification of variable lymphocyte receptors in the agnathan sea lamprey. Nature , 2004, 430(6996): 174-180.
[10] Pancer Z, Saha NR, Kasamatsu J, Suzuki T, Amemiya CT, Kasahara M, Cooper MD. Variable lymphocyte receptors in hagfish. Proc Natl Acad Sci USA , 2005, 102(26): 9224-9229.
[11] Kasamatsu J, Sutoh Y, Fugo K, Otsuka N, Iwabuchi K, Kasahara M. Identification of a third variable lymphocyte receptor in the lamprey. Proc Natl Acad Sci USA , 2010, 107(32): 14304-14308.
[12] Guo P, Hirano M, Herrin BR, Li JX, Yu CL, Sadlonova A, Cooper MD. Dual nature of the adaptive immune system in lampreys. Nature , 2009, 459(7248): 796-801.
[13] 毛芝娟, 楼丹, 杨波, 吴建军, 吴雄飞, 川合研儿, 陈昌福. 用灭活哈维弧菌菌苗浸泡免疫大黄鱼的研究. 华中农业大学学报, 2004, 23(3): 326-330.
[14] Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng QD, Chen ZH, Mauceli E, Hacohen N, Gnirke A, Rhind N, di Palma F, Birren BW, Nusbaum C, Lindblad-Toh K, Friedman N, Regev A. Full-length transcriptome assembly from RNA-seq data without a reference genome. Nat Biotechnol , 2011, 29(7): 644-652.
[15] Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, Heger A, Hetherington K, Holm L, Mistry J, Sonnhammer ELL, Tate J, Punta M. Pfam: the protein families database. Nucl Acids Res , 2014, 42(D1): D222-D2230.
[16] Götz S, García-Gómez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, Robles M, Talón M, Dopazo J, Conesa A. High-throughput functional annotation and data mining with the Blast2GO suite. Nucl Acids Res , 2008, 36(10): 3420-3435.
[17] Moriya Y, Itoh M, Okuda S, Yoshizawa A, Kanehisa M. KAAS: an automatic genome annotation and pathway reconstruction server. Nucl Acids Res , 2007, 35(Suppl. 2): W182-W185.
[18] Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol , 2009, 10(3): R25.
[19] Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform , 2011, 12: 323-339.
[20] Trapnell C, Williams B A, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotech , 2010, 28(5): 511-515.
[21] Robinson MD, Oshlack A. A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol , 2010, 11(3): R25.
[22] Wang LK, Feng ZX, Wang X, Wang XW, Zhang XG. DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics , 2010, 26(1): 136-1
