As important members of the ETS superfamily, the E74-like factor (ELF) transcription factor family regulates gene transcription through a conserved ETS domain and plays critical roles in immune regulation. However, the evolutionary characteristics and functions of this family in lampreys (Lethenteron reissneri) remain unclear. In this study, the ELF gene family of lampreys (Lr-ELF1, Lr-ELF2, Lr-ELF3, and Lr-ELF5) was systematically identified, and their molecular evolutionary features and immune response functions were investigated. Phylogenetic analysis revealed evolutionary characteristics reflecting the transition from jawless to jawed vertebrates. Domain architecture, gene structure, and three-dimensional structural analyses indicated that these genes appear to be conserved among vertebrates, with their three-dimensional structures showing high similarity to the core secondary structural elements of human homologous proteins. Synteny analysis demonstrated significant differences in the genomic neighborhoods of ELF genes between lampreys and jawed vertebrates. Quantitative real-time PCR (qRT-PCR) was performed with three biological and three technical replicates; relative expression levels were calculated using the ΔCt method, and statistical analysis was carried out with GraphPad Prism 9. Quantitative real-time PCR (qRT-PCR) results suggested that the ELF gene family may be involved in immune defense. This study not only enriches our understanding of the evolution of ELF genes but also provides new insights into the roles of lamprey ELFs in immune defense.

Transforming growth factor-β-activated kinase 1 (TAK1) and TAK1-binding proteins (TABs) are critical regulators of innate immune signaling pathways. In this study, TAK1, TAB1, and TAB2 were identified from lampreys (Lethenteron reissneri), their evolutionary characteristics and immune functions were investigated. Multiple sequence alignment analysis, phylogenetic characterization, syntenic analysis, genetic structure analysis, and 3D structure analysis revealed that Lr-TAK1, Lr-TAB1, and Lr-TAB2 possess highly conserved functional domains and evolutionary features shared with vertebrate orthologs. Quantitative real-time PCR demonstrated that all three genes were significantly induced in immune-related tissues following poly(I) and lipopolysaccharide stimulation, indicating their involvement in antiviral and antibacterial immune responses. Furthermore, co-immunoprecipitation assays confirmes that Lr-TAK1 physically interacts with both Lr-TAB1 and Lr-TAB2 in HEK293T cells. These findings provide the first systematic characterization of the TAK1 and TAB gene families in lampreys, highlight the evolutionary conservation of the TAK1-TAB signaling module in early vertebrates, and offer new insights into the origin and evolution of vertebrate innate immunity.

The E2 promoter binding factors (E2Fs) are a group of transcriptional regulators that govern the cell cycle and play crucial roles in various cellular physiological processes, including proliferation and embryonic development. In this study, we identified four homologous genes-Lr-E2F3, Lr-E2F4, Lr-E2F5, and Lr-E2F8-from the lamprey (Lethenteron reissneri) genome database. Phylogenetic tree analysis was conducted to elucidate the evolutionary relationships within the E2F family across different species. Furthermore, analyses of motifs, domains, gene structures, and 3D structures reinforced the conservation of the E2F family. Notably, synteny analysis revealed that the neighboring genes of the Lr-E2Fs exhibited greater diversity compared to those in jawed vertebrates. Activity assays indicated that Lr-E2Fs may be involved in lamprey innate immunity mediated by NF-кB. Additionally, morphological observations of embryos microinjected with Cas9/sgRNA demonstrated that E2F-deficient lamprey embryos displayed embryonic lethality, suggesting that Lr-E2Fs play a significant role in lamprey embryonic development. In summary, our research not only provides new insights into the evolution of Lr-E2Fs but also offers valuable clues regarding their functional roles.

The Rels, a class of nuclear factor κB (NF-κB) complexes, regulate diverse physiological processes by modulating the transcription of effector genes. IκBs are the critical proteins that inhibit NF-κB nuclear translocation, thereby disrupting NF-κB-mediated signaling pathways. Despite this, the precise role and underlying molecular mechanisms of Rel and IκB transcriptional regulation mediated in lamprey, a member of the oldest surviving vertebrates, remain incompletely understood. In this study, we cloned and identified 4 Rels (designated Lr_Rels) and IκBs (designated Lr_IκBs) from lamprey and explored their sequence structures and evolutionary process, indicating that Lr_Rels and Lr_IκBs represent ancestral lineages in vertebrates, and the dimerization domain (DD) might be crucial for Lr_Rels’ function. Immunoreactivity assays demonstrated a significant induction of Lr_Rel1 expression across various lamprey tissues following LPS and polyinosinic–polycytidylic acid (poly (I:C)) challenge. Functional characterization revealed that Lr_Rel1 mediates the NF-κB signaling through nuclear translocation and sequence-specific recognition, with its activity being inhibited by Lr_IκBs. Furthermore, the Rel homology region (RHR) and transcriptional activation domain (TAD) were identified as key elements for Lr_Rel1 function. Thirteen target genes of Lr_Rel1 were also identified, each containing conserved κB-binding sites within their promoter regions. Our study revealed the cooperation between Lr_Rel and Lr_IκBs, providing insights into the molecular mechanisms of lamprey Rel protein in the immune regulation signaling pathway.

Positive regulatory domain member (PRDM) family proteins play important roles in nervous system development, neural stem cell proliferation and differentiation, and central nervous system inflammation. The unique evolutionary position of the lamprey (Lethenteron reissneri) as one of the oldest jawless vertebrates makes it an ideal animal model for understanding vertebrate evolution. Nevertheless, the evolutionary characteristics of PRDM genes have not yet been demonstrated in lampreys. In this study, we identified PRDM1, 4, 5, 8, 12, 14, and 15 genes in the lamprey genomes and also investigated their evolutionary relationships through phylogenetic analysis. The characterization of PRDM genes appears to be conserved among vertebrates, as indicated by protein structural domain, motif, and 3D structure analysis. Genomic synteny analysis revealed that lamprey PRDM neighbor genes are significantly different from those of jawed vertebrates. Real-time quantitative results demonstrate that the PRDM gene family may be involved in immune defense and spinal cord injury (SCI) repair. This study not only enriches the understanding of PRDM gene evolution but also provides new clues for the Lr-PRDMs' roles in immune defense and SCI.

The inhibitor of kappa B kinase (IKK) complex, a key regulator of the nuclear factor-κB (NF-κB) signaling pathway, plays a pivotal role in immune function. This study identified two homologous genes, Lr-IKK alpha (Lr-IKKα) and Lr-TANK-binding kinase 1 (Lr-TBK1), from the Lethenteron reissneri genome database. Phylogenetic analysis revealed the evolutionary relationships of the IKK family across species. Further examination of motifs, domains, gene structures, and 3D structures confirmed the conservation of Lr-IKKα and Lr-TBK1. In vitro assays suggest that both proteins may be involved in lamprey NF-κB-mediated innate immunity. Luciferase assays demonstrated that Lr-IKKα significantly enhances NF-κB transcriptional activity. Co-immunoprecipitation experiments further indicated that Lr-IKKα may regulate the NF-κB signaling pathway through interaction with lamprey inhibitor of NF-κB molecules (Lr-IκB). Collectively, these findings provide novel insights into the evolutionary mechanisms of Lr-IKKα and Lr-TBK1, while offering critical evidence of their functional roles.

AF4/FMR2 family member (AFF) proteins are a group of transcriptional regulators that can regulate gene transcription and play an important role in cellular physiological processes such as proliferation and differentiation. The transcriptome data of the lamprey spinal cord injury were analyzed in previous research. We then identified a hub gene, Lr-AFF3, from this dataset. Phylogenetic tree analysis determined the evolutionary relationships of the AFF gene family across different species. In addition, analysis of motifs, domains, and 3D structures further confirmed the conservatism of the AFF gene family. In particular, the gene structure of the AFF3 gene was not conserved, possibly because of intron insertion. It was also found that the neighboring genes of the Lr-AFF3 gene had a higher diversity than that in jawed vertebrates through synteny analysis. The results of the MTT and EdU experiments showed that the C-terminal homology domain (CHD) and N-terminal homology domain (NHD) of Lr-AFF3 promoted cell proliferation. In summary, our research will not only provide new insights into the origin and evolution of the AFF gene family in different species, but also provide new clues for the functions of Lr_AFF3.

In humans and other adult mammals, axon regeneration is difficult in axotomized neurons. Therefore, spinal cord injury (SCI) is a devastating event that can lead to permanent loss of locomotor and sensory functions. Moreover, the molecular mechanisms of axon regeneration in vertebrates are not very well understood, and currently, no effective treatment is available for SCI. In striking contrast to adult mammals, many nonmammalian vertebrates such as reptiles, amphibians, bony fishes and lampreys can spontaneously resume locomotion even after complete SCI. In recent years, rapid progress in the development of next-generation sequencing technologies has offered valuable information on SCI. In this review, we aimed to provide a comparison of axon regeneration process across classical model organisms, focusing on crucial genes and signalling pathways that play significant roles in the regeneration of individually identifiable descending neurons after SCI. Considering the special evolutionary location and powerful regenerative ability of lamprey and zebrafish, they will be the key model organisms for ongoing studies on spinal cord regeneration. Detailed study of SCI in these model organisms will help in the elucidation of molecular mechanisms of neuron regeneration across species.

Jawless vertebrates have been studied for the evolutionary origin of adaptive immunity because of their pivotal position in chordate phylogeny. However, the convergent evolution of adaptive immunity in vertebrates is lacking. This study evaluates the similarity of adaptive immune mechanisms between jawless and jawed vertebrates using lamprey cytidine deaminase (CDA). We identified the ancestral gene Lr-CDAs of the AID/APOBEC deaminase family and evaluated its biological function in vivo. Lr-CDA1 deletion affected the assembly of three types of variable lymphocyte receptors (VLRs). We identified a switch-like region in lamprey gVLRs bound to Lr-CDAs, which upon repression, downregulated VLRB expression. Overall, we propose that lampreys have an early form of class switch recombination (CSR) that is mediated by Lr-CDAs and acts on gVLRs, affecting the assembly, maturation, and diversity regulation of VLR genes in Lethenteron reissneri. This CSR process in lampreys is linked to tumorigenesis and chromosomal translocation markers via Lr-CDAs.

The tryptophan-kynurenine (TRP-KYN) pathway is involved in several biological functions, including immunosuppression, inflammatory response, and tumor suppression. Six TRP-KYN pathway-related genes, tryptophan 2,3-dioxygenase (TDO), indoleamine 2,3-dioxygenase 2 (IDO2), aminoadipate aminotransferase (AADAT), glutamate oxaloacetate transaminase 2 (GOT2), kynurenine monooxygenase (KMO), and kynureninase (KYNU) have been identified and cloned from the jawless vertebrate lamprey (Lampetra japonica) to gain insights into their evolution and characterization. Expression distribution showed that the key gene Lj-TDO was highly expressed in the oral gland. Real-time quantitative PCR showed that TRP-KYN pathway-related genes were significantly overexpressed after multi-stimulation. RNA interference showed that Lj-IDO2 knockdown regulated the expression of inflammatory factors. In conclusion, our study successfully clarified the ancestral features and functions of the TRP-KYN pathway, while providing valuable insights into the involvement of this pathway in the immune responses of a jawless vertebrate.

The mammalian testis and ovary possess special immunocompetence, which is central to provide protection against pathogens. However, the innate immune responses to immune challenges in lamprey gonads are poorly understood. In this study, we extracted RNA from testis and ovary tissues of lampreys at 0 hour, 8 hours and 17 days after lipopolysaccharides (LPS) stimulation and performed transcriptome sequencing. While the transcriptome profiles of the two tissues were different for the most part, genes LIP, LECT2, LAL2, GRN, ITLN, and C1q were found to be the most significantly up-regulated genes in both. Quantitative Real-time PCR (qRT-PCR) analysis confirmed that these genes were upregulated after stimulation. Furthermore, immunohistochemical staining showed that these genes in lamprey gonads are expressed in high quantities and have a specific distribution. Taken together, our results suggest that these genes could play an essential role in response of the gonads to LPS induction. This research establishes a basis for investigating the immune mechanism of vertebrate gonads and presents a fresh concept for gaining insight into the evolutionary development of jawless vertebrates.

Lamprey underwent biliary atresia (BA) at its metamorphosis stage. In contrast to patients with BA who develop progressive disease, lamprey can grow and develop normally, suggesting that lamprey has several adaptations for BA. Here we show that adaptive changes in bile acid and cholesterol metabolism are produced after lamprey BA. Among 1102 differentially expressed genes (DGEs) after BA in lamprey, many are enriched in gene ontology (GO) terms and pathways related to steroid metabolism. We find that among the DGEs related to bile acids and cholesterol metabolism, the expression of cytochrome P450 family 7 subfamily A member 1 (CYP7A1), sodium-dependent taurine cotransport polypeptide (NTCP) are significantly downregulated, whereas nuclear receptor farnesoid X receptor (FXR), multidrug resistance-associated protein 3 (MRP3), 3-hydroxy-3-methylglutaryl-CoA reductase (HMGCR), sterol O-acyltransferase 1 (SOAT1), and ATP binding cassette subfamily A member 1 (ABCA1) are remarkably upregulated. The changes in expression level are also validated by RT-qPCR. Furthermore, the level of high-density lipoprotein-cholesterol (HDL-C) and low-density lipoprotein-cholesterol (LDL-C) in juvenile serum is higher compared to larvae. Taken together, the findings collectively indicate that after BA, lamprey may maintain bile acids and cholesterol homeostasis in liver tissue by inhibiting bile acids synthesis and uptake, promoting its efflux back to circulation, and enhancing cholesterol esterification for storage as lipid droplets and its egress to form nascent HDL (nHDL). Understanding the possible molecular mechanisms of lamprey metabolic adaptation sheds new light on the understanding of the development and treatment of diseases caused by abnormal bile acid and cholesterol metabolism in humans.

Tryptophan is mainly degraded through kynurenine pathway (KP) in vertebrates which is closely related to the nerve and depression, while the studies on immunity is still limited. This study aims to explore the functions of tryptophan in the innate immunity of primitive vertebrate lamprey. MTT (3-(4,5)-dimethylthiahiazo (-z-y1)-3,5-di-phenytetrazoliumromide) assay showed that tryptophan had no obvious effect on cell viability. Tryptophan was transported into leukocytes and degraded via the KP after tryptophan supplement. Tryptophan treatment (T1x and T2x) failed to alter the total antioxidant capacity regardless of stimulation and exposure time. Real-time quantitative PCR and western blotting results revealed that tryptophan was not only able to reduce the expression of pro-inflammatory factors Lj-TNF-α, Lj-IL1β and Lj-NF-κB, but also to upregulate the expression of anti-inflammatory factor Lj-TGF-β independent of stimulation and time. In addition, tryptophan can exert immune tolerance function by inhibiting TLR-MyD88 and promoting (Indoleamine 2, 3-Dioxygenase) IDO-kynurenine-AHR (aryl hydrocarbon receptor) pathways. This study provides a new understanding for tryptophan-kynurenine metabolism and mechanism of immune tolerance function in primitive vertebrate lamprey.

Bone morphogenic protein/retinoic acid inducible neural-specific proteins (BRINPs) and astrotactins (ASTNs) are two members of membrane attack complex/perforin-like (MACPF) superfamily proteins that present high expression in the growing and mature vertebrate neurons. Lamprey has a unique evolutionary status as a representative of the oldest jawless vertebrates, making it an ideal animal model for understanding vertebrate evolution. The evolutionary origins of BRINPs and ASTNs genes in vertebrates, however, have not been shown in lampreys. Here, BRINP and ASTN genes were found in lamprey genomes and the evolutionary relationships of them were investigated by phylogenetic analysis. Protein domains, motifs, genetic structure, and crystal structure analysis revealed that the features of BRINP and ASTN appear to be conserved in vertebrates. Genomic synteny analysis indicated that lamprey BRINP and ASTN neighbor genes differed dramatically from jawed vertebrate. Real-time quantitative results illustrated that the BRINP and ASTN genes family might take part in immune defence and spinal cord injury repair. This study not only enriches a better understanding of the evolution of the BRINP and ASTN genes but also offers a foundation for exploring their roles in the development of the vertebrate central nervous system (CNS).

Progranulin (PGRN) is an autocrine growth factor that regulates cell proliferation, migration, wound healing, and tissue repair in mammals. Lamprey is the most primitive of the extant vertebrates and is regarded as the survivor of a once flourishing group of paleozoic vertebrates, with a history of more than 500 million years. To date, the evolutionary dynamics and the underlying function of the PGRNs remain largely unclear in lamprey. Here, we screened four genes encoding PGRNs from the genomes of Lethenteron reissneri and Petromyzon marinus , including one long form (named Lr-PGRN-L) and three short forms (named Lr-PGRN-S1, Lr-PGRN-S2, and Lr-PGRN-S3), and performed phylogenetic tree, functional domain, and synteny analyses to identify the evolutionary history of the four Lr-PGRNs. In addition, the expressions of the four Lr-pgrn family genes and the immune response against various pathogenic challenges were also investigated. We found that these genes were widely distributed in various tissues of lamprey and performed a variety of functions. Moreover, our results suggest that Lr-PGRN-S1 induces cell migration and proliferation, and is involved in repair after skin and spinal cord injury under appropriate conditions. Our findings are valuable because they improve the understanding of the evolutionary relationship of vertebrate pgrn genes, as well as providing new insights into the diverse and important roles of Lr-PGRNs.

Transcription factors (TFs), key regulators for gene expression, play varied yet crucial roles throughout cellular polarization, migration, proliferation, differentiation and apoptosis. An important function of TFs is acting in embryogenesis and organogenesis. To identify the candidate TFs in the progression of lamprey early embryogenesis, datasetGSE76037 was downloaded from the Gene Expression Omnibus (GEO) database and an integrated bioinformatic analysis was performed. Our findings revealed a total of 152 TFs in the dataset. The function enrichment analysis showed that these genes were mainly enriched in transcription from RNA polymerase II promoter, cell differentiation, embryonic digestive tract morphogenesis and so on. Hierarchical clustering analysis suggested the expression of TFs was significantly different during early embryogenesis. Moreover, volcano plots and Venn diagrams analysis identified 36 key TFs, which were considered to play an important role during embryogenesis. The weighted correlation network analysis (WGCNA) was constructed and the Pearson correlation coefficient was performed, indicating these TFs might involve in early development of lamprey germline by synergistically regulating each other. The result was confirmed by real-time polymerase chain reaction and Western blotting analysis. In conclusion, differentially expressed genes identified in the present study help us understand the molecular mechanisms underlying the lamprey embryogenesis and provide candidate TFs for further study of vertebrate embryonic development.

Toll/interleukin-1 receptor domain-containing adaptor molecule (TICAM) genes respond to infections. We identified TICAM-a and TICAM-b in Lampetra japonica and investigated their evolutionary history and potential function via comparative genomics and molecular evolution analyses. They are arranged in tandem and evolved from a multi-exon to a single-exon structure. Lj-TICAM-a and Lj-TICAM-b might be the ancestral gene of the vertebrate TICAM genes. Lj-TICAM-b arose via a lamprey-specific tandem duplication event. Both genes are expressed in many tissues during an immune response, and exhibit different responses to peptidoglycan, indicating their functional divergence. Simultaneous overexpression of both proteins activated nuclear factor κB expression and co-immunoprecipitation assays indicated that they might form a complex for signal transduction. However, unlike in mammals, the TICAM-dependent signaling pathway in lamprey might rely on TRAF3 rather than on TRAF6. These results suggest that both Lj-TICAM-a and Lj-TICAM-b play a role in host defenses.

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.

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.

Transcription factors are a class of proteins that regulate gene transcription and expression by binding to gene-specific sequences and play an essential role in regulating the biological activities of cells. The RHR (Rel-homology region) transcription factor family is the primary member of the IF (immunoglobulin fold) transcription factor superfamily, whose members contain the conserved Rel domain and IPT (immunoglobulin-like fold) domain. As an ancient transcription factor family, the RHR family continues differentiation on gene gain and loss through gene duplication, mutation, and silencing, accompanied with the evolution of diverse species. Natural selection has led to different rates of evolution among members of the family, and some domains of the protein family have shown unique mechanisms of evolution. However, the current reviews about the origin and differentiation of RHR family are rare. In this review, we summarize the research results on the distribution, classification, function, and evolution of the members of the RHR family in order to provide a reference and new idea for studying the evolution mechanism of the whole transcription factor family and the evolutionary relationship among species.