Identification of Glycan-Specific Variable Lymphocyte Receptors Using Yeast Surface Display and Glycan Microarrays.

The jawless vertebrates (lamprey and hagfish) evolved a novel adaptive immune system with many similarities to that found in the jawed vertebrates, including the production of antigen-specific circulating antibodies in response to immunization. However, the jawless vertebrates use leucine-rich repeat (LRR)-based antigen receptors termed variable lymphocyte receptors (VLRs) for immune recognition, instead of immunoglobulin (Ig)-based receptors. VLR genes are assembled in developing lymphocytes through a gene conversion-like process, in which hundreds of LRR gene segments are randomly selected as template donors to generate a large repertoire of distinct antigen receptors, similar to that found within the mammalian adaptive immune system. Here we describe the development of a robust platform using immunized lampreys (Petromyzon marinus) for generating libraries of anti-carbohydrate (anti-glycan) variable lymphocyte receptor B, or VLRBs. The anti-carbohydrate VLRBs are isolated using a yeast surface display (YSD) expression platform and enriched by binding to glycan microarrays through the anti-glycan VLRB. This enables both the initial identification and enrichment of individual yeast clones against hundreds of glycans simultaneously. Through this enrichment strategy a broad array of glycan-specific VLRs can be isolated from the YSD library. Subsequently, the bound yeast cells are directly removed from the microarray, the VLR antibody clone is sequenced, and the end product is expressed as a VLR-IgG-Fc fusion protein that can be used for ELISA, Western blotting, flow cytometry, and immunomicroscopy. Thus, by combining yeast surface display with glycan microarray technology, we have developed a rapid, efficient, and novel method for generating chimeric VLR-IgG-Fc proteins that recognize a broad array of unique glycan structures with exquisite specificity.

Lamprey Variable Lymphocyte Receptor Monoclonal Antibodies for Whole-Cell Surface Antigens.

Lamprey antibodies, the variable lymphocyte receptor B proteins (VLRB), have unique properties that make them promising alternatives to jawed vertebrate immunoglobulin domain antibodies. These leucine-rich repeat proteins exhibit a diversity on par with that of jawed vertebrate antibodies but are structurally completely distinct. VLRB antibodies have been successfully raised to a variety of antigens. A procedure for high-throughput screening of full-length lamprey VLRB libraries using whole cells is described here. Lamprey antibodies against cell surface antigens can be generated and screened quickly using this method.

Identification and evolution of transcription factors RHR gene family (NFAT and RBPJ) involving lamprey (Lethenteron reissneri) innate immunity.

Nuclear factor of activated T cells (NFAT) and recombination signal binding protein (RBP) belong to the family of Rel homology region (RHR) transcription factors which regulate the expression of genes involved in different aspects of the immune response. To gain insights into the evolution and characterisation of RHR genes in lampreys, a jawless vertebrate, four RHR genes, including nuclear factor of activated T cells (NFAT) and recombination signal binding protein for immunoglobulin kappa J region (RBPJ), have been identified and cloned from the lamprey (Lethenteron reissneri) database. Evolutionary relationships of NFAT and RBPJ genes among different species were determined through molecular phylogenetic analysis. Motif, genetic structure, and tertiary structure analyses showed that NFATs and RBPJ are conserved and contain RHD and IPT domains. Moreover, synteny analysis showed that the neighbourhood genes of Lr-NFATs and Lr-RBPJ have undergone significant changes compared to jawed vertebrates. Real-time quantitative results demonstrated that the RHR gene family plays a significant role in immune defence. This study provides a new understanding of the origin and evolution of the RHR gene family in different species.

Identification, molecular evolution, and expression analysis of the transcription factor Smad gene family in lamprey.

Transcription factor small mothers against decapentaplegic (Smad) family SMAD proteins are the essential intracellular signal mediators and transcription factors for transforming growth factor ß (TGF-ß) signal transduction pathway, which usually exert pleiotropic actions on cell physiology, including immune response, cell migration and differentiation. In this study, the Smad family was identified in the most primitive vertebrates through the investigation of the transcriptome data of lampreys. The topology of phylogenetic tree showed that the four Smads (Smad1, Smad3, Smad4 and Smad6) in lampreys were subdivided into four different groups. Meanwhile, homology analysis indicated that most Smads were conserved with typical Mad Homology (MH) 1 and MH2 domains. In addition, Lethenteron reissneri Smads (Lr-Smads) adopted general Smads folding structure and had high tertiary structural similarity with human Smads (H-Smads). Genomic synteny analysis revealed that the large-scale duplication blocks were not found in lamprey genome and neighbor genes of lamprey Smads presented dramatic differences compared with jawed vertebrates. Importantly, quantitative real-time PCR analysis demonstrated that Smads were widely expressed in lamprey, and the expression level of Lr-Smads mRNA was up-regulated with different pathogenic stimulations. Moreover, depending on the weighted gene co-expression network analysis (WGCNA), four Lr-Smads were identified as two meaningful modules (green and gray). The functional analysis of these two modules showed that they might have a correlation with ployI:C. And these genes presented strong positive correlation during the immune response from the results of Pearson's correlation analysis. In conclusion, our results would not only enrich the information of Smad family in jawless vertebrates, but also lay the foundation for immunity in further study.

Novel lamprey antibody recognizes terminal sulfated galactose epitopes on mammalian glycoproteins.

The terminal galactose residues of N- and O-glycans in animal glycoproteins are often sialylated and/or fucosylated, but sulfation, such as 3-O-sulfated galactose (3-O-SGal), represents an additional, but poorly understood modification. To this end, we have developed a novel sea lamprey variable lymphocyte receptor (VLR) termed O6 to explore 3-O-SGal expression. O6 was engineered as a recombinant murine IgG chimera and its specificity and affinity to the 3-O-SGal epitope was defined using a variety of approaches, including glycan and glycoprotein microarray analyses, isothermal calorimetry, ligand-bound crystal structure, FACS, and immunohistochemistry of human tissue macroarrays. 3-O-SGal is expressed on N-glycans of many plasma and tissue glycoproteins, but recognition by O6 is often masked by sialic acid and thus exposed by treatment with neuraminidase. O6 recognizes many human tissues, consistent with expression of the cognate sulfotransferases (GAL3ST-2 and GAL3ST-3). The availability of O6 for exploring 3-O-SGal expression could lead to new biomarkers for disease and aid in understanding the functional roles of terminal modifications of glycans and relationships between terminal sulfation, sialylation and fucosylation.

Analysis of the lamprey genotype provides insights into caspase evolution and functional divergence.

The cysteine-containing aspartate specific proteinase (caspase) family plays important roles in apoptosis and the maintenance of homeostasis in lampreys. We conducted genomic and functional comparisons of six distinct lamprey caspase groups with human counterparts to determine how these expanded molecules evolved to adapt to the changing caspase-mediated signaling pathways. Our results showed that lineage-specific duplication and rearrangement were responsible for expanding lamprey caspases 3 and 7, whereas caspases 1, 6, 8, and 9 maintained a relatively stable genome and protein structure. Lamprey caspase family molecules displayed various expression patterns and were involved in the innate immune response. Caspase 1 and 7 functioned as a pattern recognition receptor with a broad-spectrum of microbial recognition and bactericidal effect. Additionally, caspases 1 and 7 may induce cell apoptosis in a time- and dose-dependent manner; however, apoptosis was inhibited by caspase inhibitors. Thus, these molecules may reflect the original state of the vertebrates caspase family. Our phylogenetic and functional data provide insights into the evolutionary history of caspases and illustrate their functional characteristics in primitive vertebrates.

The evolution and functional characterization of CXC chemokines and receptors in lamprey.

Chemokines are a large family of soluble peptides guiding cell migration in development and immune defense. They interact with chemokine receptors and are essential for the coordination of cell migration in diverse physiological processes. The CXC subfamily is one of the largest groups in the chemokine family and consists of multiple members. In this study, we identified homologues of three chemokine ligands (CXCL8, CXCL_F5 and CXCL12) and two CXC receptor like molecules (CXCR_L1 and CXCR_L2) in lamprey. Sequence analysis revealed that they share the same genomic organization with their counterparts in jawed vertebrates but synteny was not conserved. Lamprey CXCL8 and CXCL12 have four conserved cysteine residues whilst the CXCL_F5 has two additional cysteine residues. In addition, CXCL_F5 is evolutionarily related to the fish specific CXC chemokine groups previously identified and contains multiple cationic aa residues in the extended C- terminal region. The two CXCRs possess seven transmembrane domains and conserved structural elements for receptor activation and signaling, including the DRYXXI(V)Y motif in TM2, the disulphide bond connecting ECL2 and TM3, the WXP motif in TM6 and NPXXY motif in TM7. The identified CXC chemokines and receptors were constitutively expressed in tissues including the liver, kidney, intestine, heart, gills, supraneural body and primary leukocytes, but exhibited distinct expression patterns. Relatively high expression was detected in the gills for CXCL8, CXCL_F5 and CXCR_L1 and in the supraneural body for CXCL12 and CXCR_L2. All the genes except CXCL12 were upregulated by stimulation with LPS, pokeweed and bacterial infection, and the CXCL8 and CXCL_F5 was induced by poly (I:C). Functional analysis showed that the CXCL8 and CXCL_F5 specifically interacted with CXCR_L1 and CXCR_L2, respectively. Our results demonstrate that the CXC chemokine system had diversified in jawless fish.

Molecular Characterization and Chemotactic Function of CXCL8 in Northeast Chinese Lamprey (Lethenteron morii).

Chemokine-induced chemotaxis of leukocytes is an important part of the innate immunity and has been shown to mediate inflammation in all groups of jawed vertebrates. For jawless vertebrates, hagfish leukocytes are known to show chemotaxis toward mammalian complement anaphylotoxin and Gram-negative bacteria lipopolysaccharide. However, whether chemokines mediate chemotaxis of leukocytes in jawless vertebrates has not been conclusively examined. Here, we show C-X-C motif chemokine ligand 8 (CXCL8, also named interleukin 8) of the Northeast Chinese lamprey (Lethenteron morii) (designated as LmCXCL8) induces chemotaxis in its leukocytes. We identified LmCXCL8 and found it possesses the characteristic N-terminal cysteine residues and GGR (Gly-Gly-Arg) motif. The Lmcxcl8 gene was found to be expressed in all examined tissues, and its expression was inducible in the lamprey challenged by an infectious bacterium, Pseudomonas aeruginosa. A recombinant LmCXCL8 protein elicited concentration-dependent chemotaxis in peripheral blood leukocytes isolated from the Northeast Chinese lamprey. Based on these results, we conclude that LmCXCL8 is a constitutive and inducible acute-phase cytokine that mediates immune defense and trace the chemotactic function of chemokine to basal vertebrates.

Molecular Evolution of Apolipoprotein Multigene Family and the Original Functional Properties of Serum Apolipoprotein (LAL2) in Lampetra japonica.

Apolipoprotein (APO) genes represent a large family of genes encoding various binding proteins associated with plasma lipid transport. Due to the long divergence history, it remains to be confirmed whether these genes evolved from a common ancestor through gene duplication and original function, and how this evolution occurred. In this study, based on the phylogenetic tree, sequence alignment, motifs, and evolutionary analysis of gene synteny and collinearity, APOA, APOC, and APOE in higher vertebrates may have a common ancestor, lamprey serum apolipoprotein LAL1 or LAL2, which traces back to 360 million years ago. Moreover, the results of immunofluorescence, immunohistochemistry, and flow cytometry show that LAL2 is primarily distributed in the liver, kidney, and blood leukocytes of lampreys, and specifically localized in the cytoplasm of liver cells and leukocytes, as well as secreted into sera. Surface plasmon resonance technology demonstrates that LAL2 colocalizes to breast adenocarcinoma cells (MCF-7) or chronic myeloid leukemia cells (K562) associated with lamprey immune protein (LIP) and further enhances the killing effect of LIP on tumor cells. In addition, using quantitative real-time PCR (Q-PCR) and western blot methods, we found that the relative mRNA and protein expression of lal2 in lamprey leukocytes and sera increased significantly at different times after stimulating with Staphylococcus aureus, Vibrio anguillarum, and Polyinosinic-polycytidylic acid (Poly I:C). Moreover, LAL2 was found to recognize and bind to gram-positive bacteria (Staphylococcus aureus and Bacillus cereus) and gram-negative bacteria (Escherichia coli) and play an important role in the antibacterial process. All in all, our data reveals a long, complex evolutionary history for apolipoprotein genes under different selection pressures, confirms the immune effect of LAL2 in lamprey sera against pathogens, and lays the foundation for further research regarding biological functions of lamprey immune systems.

A novel protein upstream stimulatory factor 2 identified in lamprey, Lethenteron reissneri.

Upstream stimulatory factors are kinds of multi-functional transcription factors, which are expressed in eukaryotes widely, including Upstream stimulatory factor 1 (USFl) and upstream stimulatory factor 2 (USF2). USF protein has a typical basic helix-loop-helix leucine zipper (b-HLH-LZ) structure, which is involved in cell cycle, cell proliferations, glucose and lipid metabolism, and other biochemical processes. Although the USF family is an important regulator of cellular processes, little is known about the USF genes of lampreys, especially their evolutionary relationships, expression profiles, and biological functions. Here, an upstream stimulatory factor 2 (USF2) homolog from lamprey (Lethenteron reissneri) was identified and characterized (designated as L-USF2) because it is closer to USF2 subfamily than to USF1 subfamily. The cDNA fragment of L-USF2 has an open reading frame (ORF) of 765-bp length, encodes 254 amino acids, and contains an HLH domain at the c-terminal of amino acids. Meanwhile, motifs and genetic structure analysis reveal that USF2 gene exons are conserved. Moreover, the 3D structure analysis indicates that L-USF2 adopts the general USF2 folding and has a high structural similarity with H-USF2. The synteny results showed that the L-USF2 adjacent gene changed greatly compared with the jaw vertebrates. By real-time quantitative experiment and Western blot analysis, we found that L-USF2 gene played a significant role in the immune responses. This study not only provides us with a further understanding of the evolution and function of the USF gene family but also provides a basis for exploring its immune responses and immune defenses in lampreys.

VLRs expression were significantly affected by complement C3 knockdown morphants in Lampetra morii.

The complement component 3 of the lamprey, a jawless vertebrate, functions as an opsonin during the phagocytosis of rabbit red cells. Furthermore, lamprey C3 may be activated and cleaved into C3b, which is attached to the surface of target cells in the cytolytic process. However, the mechanism mediating the biological function of C3 in the lamprey is unknown. To our knowledge, this study is the first to show that variable lymphocyte receptors (VLRs) expression were significantly affected by complement C3 knockdown morphants in Lampetra morii. We identified the C3 gene in the lamprey genome based on its orthologs, conserved synteny, functional domains, phylogenetic tree, and conserved motifs. Additionally, we determined the optimal infection concentration of Aeromonas hydrophila to perform immune stimulation experiments in the lamprey larvae. The quantitative real-time polymerase chain reaction and immunofluorescence analyses revealed that the expression of Lampetra morii C3 (lmC3) was significantly upregulated in the larvae infected with 107 CFU/mL of A. hydrophila. The lmC3 morphants (lmC3 MO) of lamprey larvae were generated by morpholino-mediated knockdown. The lmC3 MO larvae were highly susceptible to A. hydrophila infection, which indicated that lmC3 is critical in lamprey immune response. The expression of a selected panel of orthologous genes was comparatively analyzed in the infected wild type, infected lmC3 MO, infected control MO, uninfected wild type and uninfected lmC3 MO one-month-old ammocoete larvae. The knockdown of lmC3 strongly affected the expression of VLRA+/VLRB+/VLRC+-associated genes, which was also confirmed by immunohistochemical analysis. Thus, VLR expression were significantly affected by complement C3 knockdown morphants in Lampetra morii.

Lamprey VLRB participates in pathogen detection, VLRB/L-BLNK/L-NF-κB (B-like cells) signal transduction, and development.

In jawed vertebrates, B cell receptors (BCR) are primary pathogen detectors that activate downstream signaling pathways to express adaptive immune effectors. In jawless vertebrates, the variable lymphocyte receptors (VLR) B positive lymphocytes can express and secrete specific VLRB molecules in an analogous manner to that of immunoglobulins by B cells in jawed vertebrates. Our study is the first to demonstrate the possibility of incubation of fertilized eggs and artificial breeding of Lampetra morii larvae throughout their life cycle under laboratory condition. We also found that VLRB, lamprey B-cell linker (L-BLNK), and lamprey nuclear factor-kappa B (L-NF-κB) play key roles in early larval development. Aeromonas hydrophila was found to be a lethal pathogen of L. morii larvae causing rapid infection at a concentration of 107 cfu/mL qRT-PCR results revealed that gene expression levels of VLRB, L-BLNK, and L-NF-κB were up-regulated significantly. Ten-day infection trials showed that VLRB, L-BLNK, and L-NF-κB are crucial for lamprey immune response. Furthermore, the expression levels of L-BLNK and L-NF-κB were down-regulated drastically both at mRNA and protein levels after bacterial infection than in the naive group of VLRB morphants. A similar expression pattern of VLRB and L-BLNK was found in L-NF-κB morphants post bacterial infection. The results were strikingly different in the other two morphants. The VLRB and L-NF-κB expression levels were found to be down-regulated at mRNA and protein levels by less than 30% and 45%, respectively, in L-BLNK morphants compared to those in the naive group. These results indicate that L-BLNK and L-NF-κB might participate in VLRB-mediated immune response. Additionally, in VLRB morphants, the mRNA expression levels of some genes, especially the ones expressed in VLRB+ lymphocytes but not in VLRA+ lymphocytes, were found to be affected. Therefore, these findings of B-like lymphocytes in lamprey offer key evidence with regard to the evolution of adaptive immunity.

Characterization of lamprey (Lampetra japonica) tnfr10-like gene: A potential granulocyte marker molecule and its immune functions.

Tumor necrosis factor receptor superfamily (TNFRSF) is an ancient protein superfamily. By binding to tumor necrosis factor (TNF), it can participate in inflammatory response, apoptosis, lymphocyte homeostasis and tissue development. Seven TNFR members have previously been identified in lampreys but detailed functions of TNFR members are not yet to be resolved. Here, we demonstrate some of the distinguishing features of TNFR10-like member which belongs to TNFRSF. The immunohistochemical results indicate that the TNFR10-like protein is abundant in vascular epithelial cells of the lamprey typhlosole and gills. The expression of tnfr10-like gene has a significantly increased at transcription level after Vibrio anguillarum, Staphylococcus aureus and Poly (I:C) stimulation. Notably, TNFR10-like is specifically expressed in the granulocytes of lamprey peripheral blood and supraneural body. Besides, overexpression tnfr10-like gene in HEK-293 T cells cause a decrease in cell activity and able to activate nuclear transcription factor-κB (NF-κB). Together, these results imply that L-TNFR10-like may play a vital role as a potential marker in lamprey granulocytes and may also be involved in regulation of immune response mediated by itself.

Complement component C1q plays a critical role in VLRA/VLRC-mediated immune response.

In jawless vertebrates, the lamprey complement component C1q (LC1q) acts as a lectin and activates lamprey complement component C3 (LC3) in association with mannose-binding lectin (MBL)-associated serine protease (MASP) via the lectin pathway. Furthermore, LC1q may interact with variable lymphocyte receptor B (VLRB) in a complex with antigens and mediate the activation of LC3, leading to cytolysis. In the present study, we found, for the first time, that LC1q plays a critical role in VLRA/VLRC-mediated immune response. Escherichia coli, Shigella flexneri, Aeromonas hydrophila, Pseudomonas plecoglossicida, Aeromonas allosaccharophila, P. luteola, Brevundimonas diminuta, and Bacillus cereus were isolated from infected Lampetra morii in our laboratory and identified using the 16s rRNA method. A. hydrophila was confirmed as a rapidly spreading lethal pathogen in the larvae of L. morii and was used in subsequent immune stimulation experiments. The results of real-time quantitative polymerase chain reaction (Q-PCR) and immunofluorescence analyses indicated that the RNA and protein expression levels of LC1q were upregulated following exposure to 107 cfu/mL of A. hydrophila, compared to the levels of the naïve group. We obtained LC1q morphants (LC1q MO) of lamprey larvae by morpholino-mediated knockdowns. We found that LC1q played key roles in the embryonic development of lamprey. The median lethal time (LT50) of LC1q MO larvae was 2 d after being exposed to the pathogens, whereas the LT50 of control MO was 5 d. The drastic decrease in LT50 values after LC1q knockdown implies that LC1q plays a critical role in lamprey immune response. Gene expression profiles of LC1q-deficient A. hydrophila, control MO A. hydrophila, wild type A. hydrophila, and naive 1-month-old ammocoetes larvae were compared by examining the expression levels of a selected panel of orthologous genes. It is worth mentioning that LC1q MO affected the VLRA+/VLRC + population genes but did not affect the VLRB + populations. Immunohistochemical data indicated that LC1q deficiency also affected VLRA and VLRC but not VLRB. Thus, LC1q plays a critical role in VLRA/VLRC-mediated immune response in lamprey.

Cytidine deaminase 2 is required for VLRB antibody gene assembly in lampreys.

The antibodies of jawless vertebrates consist of leucine-rich repeat arrays encoded by somatically assembled VLRB genes. It is unknown how the incomplete germline VLRB loci are converted into functional antibody genes during B lymphocyte development in lampreys. In Lampetra planeri larvae lacking the cytidine deaminase CDA2 gene, VLRB assembly fails, whereas the T lineage-associated VLRA and VLRC antigen receptor gene assemblies occur normally. Thus, CDA2 acts in a B cell lineage-specific fashion to support the somatic diversification of VLRB antibody genes. CDA2 is closely related to activation-induced cytidine deaminase (AID), which is essential for the elaboration of immunoglobulin gene repertoires in jawed vertebrates. Our results thus identify a convergent mechanism of antigen receptor gene assembly and diversification that independently evolved in the two sister branches of vertebrates.

Development of smart anti-glycan reagents using immunized lampreys.

Studies on the expression of cellular glycans are limited by a lack of sensitive tools that can discriminate specific structural features. Here we describe the development of a robust platform using immunized lampreys (Petromyzon marinus), which secrete variable lymphocyte receptors called VLRBs as antibodies, for generating libraries of anti-glycan reagents. We identified a wide variety of glycan-specific VLRBs detectable in lamprey plasma after immunization with whole fixed cells, tissue homogenates, and human milk. The cDNAs from lamprey lymphocytes were cloned into yeast surface display (YSD) libraries for enrichment by multiple methods. We generated VLRB-Ig chimeras, termed smart anti-glycan reagents (SAGRs), whose specificities were defined by microarray analysis and immunohistochemistry. 15 VLRB antibodies were discovered that discriminated between linkages, functional groups and unique presentations of the terminal glycan motif. The development of SAGRs will enhance future studies on glycan expression by providing sequenced, defined antibodies for a variety of research applications.

Identification, recombinant expression and immunological study of Lja-SHP2 in Lampetra japonica.

The protein tyrosine phosphatase SHP2 of higher vertebrates, encoded by ptpn11 gene, catalyzes the dephosphorylation of tyrosine phosphorylation site, and plays regulatory roles in various signaling pathways by cooperating with other protein tyrosine kinase. Previous studies have shown that SHP2 plays an important role in the activation and signal transduction of T and B cells in higher vertebrates. To study the role of a SHP2 homologous molecule of lampreys (Lja-SHP2) in immune response, we cloned and expressed the open reading frame sequence of Lja-SHP2 gene in prokaryotic expression vector pET-32a. The recombinant protein was successfully expressed in E. coli and the rabbit-derived polyclonal antibody was prepared. Lampetra japonica were immunized with mixed bacteria, and the mRNA and protein of Lja-SHP2 in immune-related cells and tissues were detected by real-time quantitative PCR and Western blotting after immunization. The Lja-SHP2 mRNA and protein were not significantly affected in leukocytes and supraneural myeloid bodies, but up-regulated significantly in gill tissues (P<0.05) after challenged by mixed bacteria, which indicated that Lja-SHP2 mainly participates in the immune response of gill tissues after mixed bacteria stimulation. To further investigate whether Lja-SHP2 level was affected in three lymphocyte subsets, the B-cell mitogen lipopolysaccharide (LPS) and T-cell mitogen phytohaemagglutinin (PHA) were employed to boost the immune response in L. japonica. LPS immune stimulation increased Lja-SHP2 in leucocytes significantly compared with the control group, and but had a marginal effect on Lja-SHP2 expression in gills and supraneural myeloid bodies. PHA immune stimulation could up-regulate Lja-SHP2 level in leukocytes, gill tissues and supraneural myeloid bodies. The change of Lja-SHP2 was especially dramatical in leukocytes, which was about 2.5 times higher than that in the control group, suggesting that Lja-SHP2 is involved in the lamprey immune response mediated by PHA. Consistent with the previous finding that PHA could induce the activation of VLRA+ lymphocytes, our results showed that Lja-SHP2 might be included in the immune response of VLRA+ lymphocytes mediated by PHA in gills. This research will benefit exploring the functions of Lja-SHP2 in the immune response of lamprey and will provide clues for understanding the phylogenesis of SHP2 molecular family, and its roles in the early occurrence and evolution of adaptive immune system in higher vertebrates.

Comparative transcriptomic analysis provides insights into immune responses of lamprey larvae under three pathogens infections.

Based on the large number of deaths of lamprey larvae infected by pathogens in the culture process, the infection response of lamprey larvae to different pathogens was studied. Thus, high-throughput sequencing is firstly used to analyze the lamprey larvae infected by Aeromonas hydrophila, Vibrosplendidus, and Saprolegnia parasitica. A high quality of 338,628,426 Reads were obtained and 203,829 transcripts were assembled. A total of 158,001 Unigenes were annotated by Blast comparison, and a relatively comprehensive transcriptome database of lamprey larvae was established. Based on the transcriptome of lamprey larvae infected by three aquatic pathogens, classification results of COG and GO show that the genes classified into signal transduction mechanism accounted for the largest proportion, and 2316 differentially expressed genes were identified and screened out. After infection by the three pathogens, 16 genes and 19 genes expression were up-regulated and down-regulated, respectively. KEGG pathway analysis of differentially expressed genes showed that lamprey larvae were mainly affected their metabolism and oxidative phosphorylation pathways after infection by two gram-negative bacteria, Vibrosplendidus and Aeromonas hydrophila. As a fungus, Saprolegnia parasitica mainly affected ribosome synthesis in lamprey larvae. The results were validated by detecting the expression levels of 13 DEGs at 24 h after pathogens treatment using Q-PCR. The results provide valuable information for further research into the development of immunity and the innate immune response against pathogen invasion during the artificial propagation of lamprey larvae.