Full-length transcriptome analysis provides new insights into the innate immune system of Marsupenaeus japonicus

Exploring the innate immune system of Urechis unicinctus: Insights from full-length transcriptome analysis

The Echiura worm Urechis unicinctus refers to a common benthic invertebrate found in the intertidal zone of Huanghai as well as Bohai Bay. U. unicinctus is known to contain various physiologically active substances, making it highly valuable in terms of its edibility, medicinal properties, and economic potential. Nonetheless, the limited study on the immune system of U. unicinctus poses difficulties for its aquaculture and artificial reproduction. Marine invertebrates, including shellfish and U. unicinctus, are thought to primarily depend on their innate immune system for disease protection, owing to the severalinnate immune molecules they possess. Herein, we employed PacBio single-molecule real-time (SMRT) sequencing technology to perform the full-length transcriptome analysis of U. unicinctus individuals under five different conditions (room temperature (RT), low temperature (LT), high temperature (HT), without water (DRY), ultraviolet irradiation (UV)). Concequently, we identified 59,371 unigenes that had a 2,779 bp average length, 2,613 long non-coding RNAs (lncRNAs), 59,190 coding sequences (CDSs), 35,166 simple sequence repeats (SSRs), and 1,733 transcription factors (TFs), successfully annotating 90.58 % (53,778) of the unigenes. Subsequently, key factors associated with immune-related processes, such as non-self-recognition, cellular immune defenses, and humoral immune defenses, were searched. Our study also identified pattern recognition receptors (PRRs) that included 17 peptidoglycan recognition proteins (PGRPs), 13 Gram-negative binding proteins (GNBPs), 18 scavenger receptors (SRs), 74 toll-like receptors (TLRs), and 89 C-type lectins (CLTs). Altogether, the high-quality transcriptome obtained data will offer valuable insights for further investigations into U. unicinctus innate immune response, laying the foundation for subsequent molecular biology studies and aquaculture.

Full-length transcriptome assembly of black amur bream (Megalobrama terminalis) as a reference resource

BACKGROUND

The genus Megalobrama holds significant economic value in China, with M. terminalis (Black Amur bream) ranking second in production within this group. However, lacking comprehensive genomic and transcriptomic data has impeded research progress. This study aims to fill this gap through an extensive transcriptomic analysis of M. terminalis.

METHODS AND RESULTS

We utilized PacBio Isoform Sequencing to generate 558,998 subreads, totaling 45.52 Gb, which yielded 22,141 transcripts after rigorous filtering and clustering. Complementary Illumina short-read sequencing corrected 967,114 errors across these transcripts. Our analysis identified 12,426 non-redundant isoforms, with 11,872 annotated in various databases. Functional annotation indicated 11,841 isoforms matched entries in the NCBI non-redundant protein sequences database. Gene Ontology analysis categorized 10,593 isoforms, revealing strong associations with cellular processes and binding functions. Additionally, 8203 isoforms were mapped to pathways in the Kyoto Encyclopedia of Genes and Genomes, highlighting significant involvement in immune system processes and complement cascades. We notably identified key immune molecules such as alpha-2-macroglobulin and complement component 3, each with multiple isoforms, underscoring their potential roles in the immune response. Our analysis also uncovered 853 alternative splicing events, predominantly involving retained introns, along with 672 transcription factors and 426 long non-coding RNAs.

CONCLUSIONS

The high-quality reference transcriptome generated in this study provides a valuable resource for comparative genomic studies within the Megalobrama genus, supporting future research to enhance aquaculture stocks.

Identification and characterization of olfactory gene families in Macrobrachium rosenbergii based on full-length transcripts and genome sequences

The olfactory gene families include odorant binding proteins (OBPs), chemosensory proteins (CSPs), olfactory receptors (ORs), ionotropic receptors (IRs) and gustatory receptors (GRs). To investigate the molecular function of olfactory perception in Macrobrachium rosenbergii, we integrated the full-length transcripts and whole-genome sequences to identify the olfactory gene families. In this study, a total of 38,955 full-length transcripts with an N50 length of 3383 bp were obtained through PacBio SMRT sequencing. Through the annotation of full-length transcripts and whole-genome sequences, several olfactory gene families were identified, including 18 MrORs, 16 MrIRs, 151 MrIGluRs (ionotropic glutamate receptors), 2 MrVIGluRs (variant ionotropic glutamate receptors) and 3 MrCRs (chemosensory receptors). Notably, the CRs were first identified in prawns and shrimps. Additionally, the olfactory gene families in M. nipponense were identified, comprising 4 MnORs, 21 MnIRs, 79 MnIGluRs, 5 MnVIGluRs, 1 MnGR and 1 MnOBP, using the available whole-genome sequences. Meanwhile, the external morphology of the chemical sensory organs of M. rosenbergii was explored, and the presence of plumose setae (PS), hard thorn setae (HTS), bamboo shoot setae (BSS), soft thorn setae (STS) and aesthetascs (AE) on the antennules, HTS and BSS on the second antennae, and PS on the pereiopods were observed by scanning electron microscope. This study provides valuable insights for future functional studies into the olfactory perception of crustaceans and establishes a theoretical basis for molecular design breeding in M. rosenbergii.

Long noncoding RNA profiling in hepatopancreas of Pacific white shrimp and its role in response to white spot syndrome virus infection

Long noncoding RNA (lncRNA) is a potential regulator of biological processes, including immunity, reproduction, and development. Although several transcriptome studies have focused on responses of viral infections in several organisms, the role of lncRNAs in viral responses in shrimp is still unclear. Therefore, this work aimed to identify putative lncRNAs and study their role in white spot syndrome virus (WSSV) infection in white shrimp. The hepatopancreas transcriptome from WSSV infected shrimp was analyzed in silico to identify putative lncRNAs. Among 221,347 unigenes of the de novo assembled transcriptome, 44,539 putative lncRNAs were identified, 32 of which were differentially expressed between WSSV-infected and control shrimp. Five candidate lncRNAs were validated for their expressions in shrimp tissues and in response to WSSV infection. Lnc164 was chosen for further investigation of its role in WSSV infection. Knockdown of lnc164 prolonged survival of shrimp when challenged with WSSV, suggesting a role in shrimp immunity. In addition, lnc164 was not directly involved in the control of total hemocytes and viral loads in hemolymph of WSSV-infected shrimp. A set of lnc164-regulated genes was obtained by RNA sequencing among which 251 transcripts were differentially expressed between lnc164 knockdown and control shrimp. Six immune-related genes were validated for their expression profiles. Our work sheds light on lncRNA profiles in L. vannamei in response to WSSV infection and paves the way to a functional study of lnc164 in host antiviral response.

Comprehensive transcripts analysis based on single-molecule real-time sequencing and Illumina sequencing provides insights into the mining of Toll-like receptor family in Schizothorax lissolabiatus

Schizothorax lissolabiatus is an economically important cold-water fish species in southwestern China. Because of water pollution and habitat destruction, the number of wild populations has dramatically decreased. In this study, we used PacBio single-molecule real-time (SMRT) sequencing and Illumina sequencing to generate the first full-length transcriptome and transcriptome, respectively. A total of 19 310 polished consensus reads (PC) were obtained, with an average length of 1379 bp and an N50 length of 1485 bp. Meanwhile, 12 253 transcripts were successfully annotated as known homologous genes. The pathway annotation indicated that the enrichment and expression of most genes were mainly related to membrane, signal transduction and binding, and immune response. Furthermore, we identified 16 Toll-like receptors (TLRs) by mining the data from the transcripts. Phylogeny analysis showed that S. lissolabiatus TLR genes (slTLRs) supported the classification of TLRs into six families as in other vertebrates. Selection pressure analyses showed that 16 slTLRs revealed purification selection at the overall evolutionary selection. Further, positive selection signals were still detected in eight slTLRs, and most of the positive selection sites were located in the leucine-rich repeat region (LRR domain) associated with the recognition of pathogenic microorganisms, indicating that the function of these slTLR genes may be affected. Tissue specific expression analysis showed all slTLRs are present in kidney, spleen and liver but the relative expression varied among tissues. In conclusion, this study not only provided a valuable resource of transcripts for further research on S. lissolabiatus, but also contributed to improve the current understanding of the evolutionary history of immune-related genes and the TLR gene family in S. lissolabiatus.

Full-length transcriptome analysis provides new insights into the diversity of immune-related genes in the threatened freshwater shellfish Solenaia oleivora

Solenaia oleivora, a valuable and rare bivalve endemic to China, is becoming a threatened freshwater sepcies. However, the lack of research on its genome and immune system will hinder advances in its conservation and artificial breeding. In this study, we obtained the full-length transcriptome of S. oleivora using PacBio sequencing. A total of 21,415 transcripts with an average length of 1,726 bp were generated. Among these transcripts, 12,084 had coding sequences (CDS), of which 8,639 were annotated in 6 databases. The structure analysis identified 625 transcript factors (TFs), 8,005 long non-coding RNAs (lncRNAs), and 5,288 simple sequences repeat (SSRs). Meanwhile, massive immune genes were identified from the transcriptome of S. oleivora. In terms of non-self-identification, 97 transcripts of pattern recognition receptors (PRRs) were discovered, including peptidoglycan recognition proteins (PGRPs), gram-negative bacteria binding proteins (GNBPs), toll-like receptors (TLRs), scavenger receptors (SRs), galectins (GALs), C-type lectins (CLTs), and fibrinogen-related protein (FREPs). For pathogen elimination, 7 transcripts related to antimicrobial peptides, lysozymes, and lysosomal enzymes were identified. Moreover, 33 complement-associated transcripts were found. This study enriched the genome resources of S. oleivora and provided new insights for the study of the immune system of S. oleivora.

LncRNA-mRNA integrated profiling analysis in response to white spot syndrome virus in hepatopancreas in Penaeus japonicus

Penaeus japonicas is an important shrimp species, which is exposed to stressors including a variety of epidemic diseases. To date, little is known about the mechanisms involved in the response to white spot syndrome virus (WSSV) mediated by long non-coding RNAs (lncRNAs). A total of 6544 putative lncRNAs were identified in the hepatopancreas in P. japonicas, which provides a useful lncRNA reference resource for use in future studies. In addition, a total of 444 differentially expressed mRNAs and 457 differentially expressed lncRNAs were identified at 6, 12, and 24 h after WSSV infection in the hepatopancreas of P. japonicas. Functional enrichment analysis showed that the differentially expressed mRNAs were enriched in terms related to immune response and viral infectivity such as defense response, aminopeptidase activity, whereas the differentially expressed lncRNA partner genes were enriched in ubiquitin-dependent protein catabolic process, lipoprotein metabolic process, and antigen processing and presentation. Moreover, several lncRNAs were induced by WSSV infection, indicating these lncRNAs might participate in regulating many immune processes referring to their partner genes. Co-expression analysis of the lncRNAs and their partner genes identified some high lncRNA-mRNA correlations. These results suggest that WSSV stimulates the immune response in the hepatopancreas potentially through an important coding and non-coding gene network, thereby providing valuable information regarding non-coding responses to WSSV in Penaeus species.

Full-length transcriptome sequencing of Heliocidaris crassispina using PacBio single-molecule real-time sequencing

The lack of high-throughput sequencing data makes the research progress of Heliocidaris crassispina slow. Therefore, we used PacBio single-molecule real-time sequencing to generate the first full-length transcriptome. Here, 31,181 isoforms were obtained, with an average length of 2383.20 and a N50 length of 2732 bp. Meanwhile, 764 alternative splicing (AS) events, 5098 long-noncoding RNAs (LncRNAs), 6978 simple sequence repeats (SSRs), and 950 hypothetical transcript factors (TFs) were identified. Moreover, five key innate immune pattern recognition receptors (PRRs), including toll-like receptor (TLR), NACHT domain and leucine-rich repeat (NLR), scavenger receptor cysteine-rich (SRCR), peptidoglycan recognition proteins (PGRP), and gram-negative binding proteins (GNBP), were searched in the transcriptome. In addition, 37 isoforms enriched in KEGG and GO immune systems were also detected. The study provid abundant data support for the current research on H. crassispina.

Antibacterial activity of four recombinant carbohydrate recognition domain proteins identified from the kuruma shrimp Marsupenaeus japonicus

The carbohydrate recognition domain (CRD) is the key component of C-type lectins (CTLs) with the capacity to recognize and eliminate invading pathogens. Herein, the recombinant proteins of four CRDs identified from the kuruma shrimp, Marsupenaeus japonicus, were produced and purified by an Escherichia coli expression system and affinity chromatography. Bacterial binding and antibacterial assays showed that the four CRDs displayed various bacterial binding and antibacterial activities against different bacteria. Among the four recombinant CRDs, His-CRD2-3 exhibited the broadest spectrum of bacterial binding and antibacterial activities against gram-negative bacteria (Vibrio parahaemolyticus, V. alginolyticus and V. harveyi) and gram-positive bacteria (Staphylococcus aureus and Micrococcus lysodeikticus). Moreover, the four recombinant CRDs showed different capacities to regulate the expression of several immune effector genes (MjCTL3, MjCTL4, MjCTL, Mjily and Mjsty), among which His-CRD2-3 displayed broader and stronger inductive effects on these immune effector genes. This study indicated that the four CRDs participated in immune defense by binding and killing bacteria and regulating the transcription of other immune effector genes. In addition, our results suggested that His-CRD2-3 might be a promising agent for the prevention and treatment of bacteriosis.

SMRT Sequencing of the Full-Length Transcriptome of the Coelomactra antiquata

Coelomactra antiquata is an important aquatic economic shellfish with high medicinal value. However, because C. antiquata has no reference genome, a lot of molecular biology research cannot be carried out, so the analysis of its transcripts is an important step to study the regulatory genes of various substances in C. antiquata. In the present study, we conducted the first full-length transcriptome analysis of C. antiquata by using PacBio single-molecule real-time (SMRT) sequencing technology. The results identified a total of 39,209 unigenes with an average length of 2,732 bp, 23,338 CDSs, 251 AS events, 9,881 lncRNAs, 20,106 SSRs, and 2,316 TFs. Subsequently, 59.22% (23,220) of the unigenes were successfully annotated, of which 23,164, 18,711, 15,840, 13,534, and 13,474 unigenes could be annotated using NR, Swiss-prot, KOG, GO, and KEGG databases, respectively. This study lays the foundation for the follow-up research of molecular biology and provides a reference for studying the more medicinal value of C. antiquata.

Characterization of the foxl2 gene involved in the vtg expression in mud crab (Scylla paramamosain)

Forkhead box protein L2 (Foxl2) is involved in multiple physiological processes, such as ovarian development, granulosa cell differentiation, ovarian follicle development, and oocyte growth. In this study, a Spfoxl2 gene encoded 530 amino acid protein with characteristic forkhead (FH) domain was identified from transcriptome data of mud crab Scylla paramamosain and validated the accuracy by PCR technology. Meanwhile, the orthologues of the Spfoxl2 gene in other 14 crustacean species were identified with the same method. Further multiple sequence alignment analysis revealed the Foxl2 was highly conserved, especially in the FH domain, even completely identical in several species. Besides, the semi-quantitative PCR (Sq-PCR) result showed Spfoxl2 gene was mainly expressed in the gonad (testis and ovary). Further quantitative real-time PCR (qRT-PCR) result demonstrated its expression level in the testis was significantly higher than that in the ovary (p < 0.01). In addition, the qRT-PCR result showed that in zoea V, megalopa, and larval I, the expression level of Spfoxl2 in megalopa is the highest. In addition, a putative Foxl2 binding site was identified on the promoter region of Spvtg, and knockdown of Spfoxl2 mediated by RNAi technology increased the expression of Spvtg in the ovary, suggesting Spfoxl2 might be the upstream negative regulator of Spvtg. Overall, this study provided new insights into the role of Spfoxl2 in ovary development through regulating Spvtg expression in S. paramamosain.