Full-length transcriptome analysis of Litopenaeus vannamei reveals transcript variants involved in the innate immune system

In silico identification and functional study of long non-coding RNA involved in acute hepatopancreatic necrosis disease caused by Vibrio parahaemolyticus infection in white shrimp, Litopenaeus vannamei

Acute hepatopancreatic necrosis disease (AHPND) caused by toxin-producing Vibrio parahaemolyticus (VpAHPND) has severely affected shrimp production. Long non-coding RNA (lncRNA), a regulatory non-coding RNA, which can play important function in shrimp disease responses. This study aimed to identify and investigate the role of lncRNA involved in VpAHPND infection in Pacific white shrimp, Litopenaeus vannamei. From a total of 368,736 de novo assembled transcripts, 67,559 were identified as putative lncRNAs, and only 72 putative lncRNAs showed differential expression between VpAHPND-infected and normal shrimp. The six candidate lncRNAs were validated for their expression profiles during VpAHPND infection and tissue distribution using RT-qPCR. The role of lnc2088 in response to VpAHPND infection was investigated through RNA interference. The result indicated that the suppression of lnc2088 expression led to an increase in shrimp mortality after VpAHPND infection. To explore the set of genes involved in lnc2088 knockdown, RNA sequencing was performed. A total of 275 differentially expressed transcripts were identified in the hepatopancreas of lnc2088 knockdown shrimp. The expression profiles of five candidate metabolic and immune-related genes were validated in lnc2088 knockdown and VpAHPND-infected shrimp. The result showed that the expression of ChiNAG was significantly increased, while that of NCBP1, WIPF2, and NFKB1 was significantly downregulated in ds2088-injected shrimp. Additionally, the expression of NFKB1, NCBP1 and WIPF2 was significantly increased, whereas that of ChiNAG and CUL5 were significantly decreased after infection with VpAHPND. Our work identified putative lncRNA profiles in L. vannamei in response to VpAHPND infection and investigated the role of lncRNA in shrimp immunity.

Polyploidization of Indotyphlops braminus: evidence from isoform-sequencing

BACKGROUND

Indotyphlops braminus, the only known triploid parthenogenetic snake, is a compelling species for revealing the mechanism of polyploid emergence in vertebrates.

METHODS

In this study, we applied PacBio isoform sequencing technology to generate the first full-length transcriptome of I. braminus, aiming to improve the understanding of the molecular characteristics of this species.

RESULTS

A total of 51,849 nonredundant full-length transcript assemblies (with an N50 length of 2980 bp) from I. braminus were generated and fully annotated using various gene function databases. Our analysis provides preliminary evidence supporting a recent genome duplication event in I. braminus. Phylogenetic analysis indicated that the divergence of I. braminus subgenomes occurred approximately 11.5 ~ 15 million years ago (Mya). The full-length transcript resource generated as part of this research will facilitate transcriptome analysis and genomic evolution studies in the future.

Shrimp laccase degrades polycyclic aromatic hydrocarbons from an oil spill disaster in Brazil: A tool for marine environmental bioremediation

Our work aims to purify, characterize and evaluate a laccase from by-products of the shrimp farming industry (Litopenaeus vannamei) for the degradation of Polycyclic Aromatic Hydrocarbons (PAHs) from 2019 oil spill in Brazilian coast. The enzyme was purified by affinity chromatography and characterized as thermostable, with activity above 90 °C and at alkaline pH. In addition, the laccase was also tolerant to copper, lead, cadmium, zinc, arsenic, hexane and methanol, with significant enzymatic activation in acetone and 10 mM mercury. Concerning PAHs' degradation, the enzyme degraded 42.40 % of the total compounds, degrading >50 % of fluorene, C4-naphthalenes, C3-naphthalenes, C2-naphthalenes, anthracene, acenaphthene, 1-methylnaphthalene and 2-methylnaphthalene. Thus, this laccase demonstrated important characteristics for bioremediation of marine environments contaminated by crude oil spills, representing a viable and ecological alternative for these purposes.

Alternative splicing derived invertebrate variable lymphocyte receptor displays diversity and specificity in immune system of crab Eriocheir sinensis

Variable lymphocyte receptors (VLRs) play vital roles in adaptive immune system of agnathan vertebrate. In the present study, we first discover a novel VLR gene, VLR2, from an invertebrate, the Chinese mitten crab, Eriocheir sinensis. VLR2 has ten different isoforms formed via alternative splicing, which is different from that in agnathan vertebrate with the assembly of LRR modules. The longest isoform, VLR2-L, responds to Gram-positive bacteria Staphylococcus aureus challenge specifically, while shows no response to Gram-negative bacteria Vibrio parahaemolyticus challenge, confirmed by recombinant expression and bacterial binding experiments. Interestingly, VLR2s with short LRRs regions (VLR2-S8 and VLR2-S9) tend to bind to Gram-negative bacteria rather than Gram-positive bacteria. Antibacterial activity assay proves six isoforms of VLR2 have pluralistic antibacterial effects on bacteria which were never reported in invertebrate. These results suggest that the diversity and specificity of VLR2 resulted from alternative splicing and the length of the LRRs region. This pathogen-binding receptor diversity will lay the foundation for the study of immune priming. Furthermore, studying the immune function of VLR2 will provide a new insight into the disease control strategy of crustacean culture.

Genome assembly and isoform analysis of a highly heterozygous New Zealand fisheries species, the tarakihi (Nemadactylus macropterus)

Although being some of the most valuable and heavily exploited wild organisms, few fisheries species have been studied at the whole-genome level. This is especially the case in New Zealand, where genomics resources are urgently needed to assist fisheries management. Here, we generated 55 Gb of short Illumina reads (92× coverage) and 73 Gb of long Nanopore reads (122×) to produce the first genome assembly of the marine teleost tarakihi [Nemadactylus macropterus (Forster, 1801)], a highly valuable fisheries species in New Zealand. An additional 300 Mb of Iso-Seq reads were obtained to assist in gene annotation. The final genome assembly was 568 Mb long with an N50 of 3.37 Mb. The genome completeness was high, with 97.8% of complete Actinopterygii Benchmarking Universal Single-Copy Orthologs. Heterozygosity values estimated through k-mer counting (1.00%) and bi-allelic SNPs (0.64%) were high compared with the same values reported for other fishes. Iso-Seq analysis recovered 91,313 unique transcripts from 15,515 genes (mean ratio of 5.89 transcripts per gene), and the most common alternative splicing event was intron retention. This highly contiguous genome assembly and the isoform-resolved transcriptome will provide a useful resource to assist the study of population genomics and comparative eco-evolutionary studies in teleosts and related organisms.

Transcriptomics revealed the effect of astaxanthin on apoptosis and immunity of the adult prawn of Exopalaemon carinicauda

Astaxanthin (Axn), a common aquatic feed additive, can enhance immunity, improve the antioxidant capacity of the crustacean and then improve the anti-stress ability of crustaceans. Exopalaemon carinicauda (E. carinicauda) is an economically important fishery species in China that has been found that dietary Axn can significantly increase ACP and AKP compared to a control diet for shrimp hepatopancreas in this study. RNA-sequencing and comparative transcriptomic analyses were utilized to explore changes in E. carinicauda gene expression following Axn feeding. Differential gene expression analyses comparing the control and Axn groups identified 631 transcripts that were differentially expressed following Axn feeding, of which 314 and 317 were respectively upregulated and downregulated. Functional enrichment analyses of these genes revealed their enrichment in 22 Gene Ontology categories and 11 KEGG pathways. In the GO and KEGG enrichment analysis, it was found that dietary astaxanthin can regulate the gene expression level of adult E. carinicauda. Many of the signal pathways enriched by these genes are related to immunity, apoptosis and anti-stress. In addition, through KEGG enrichment analysis, it was found that dietary Axn could also regulate the amino acid metabolism of hepatopancreas of adult E. carinicauda. The comprehensive comparative transcriptomic analysis showed that Axn could improve the hepatopancreatic immunity and anti-apoptosis ability of adult E. carinicauda.

Full-length transcriptome sequencing and comparative transcriptome analysis of Eriocheir sinensis in response to infection by the microsporidian Hepatospora eriocheir

As a new generation of high-throughput sequencing technology, PacBio Iso-Seq technology (Iso-Seq) provides a better alternative sequencing method for the acquisition of full-length unigenes. In this study, a total of 22.27 gigabyte (Gb) subread bases and 128,614 non-redundant unigenes (mean length: 2,324 bp) were obtained from six main tissues of Eriocheir sinensis including the heart, nerve, intestine, muscle, gills and hepatopancreas. In addition, 74,732 unigenes were mapped to at least one of the following databases: Non-Redundant Protein Sequence Database (NR), Gene Ontology (GO), Kyoto Encyclopaedia of Genes and Genomes (KEGG), KEGG Orthology (KO) and Protein family (Pfam). In addition, 6696 transcription factors (TFs), 28,458 long non-coding RNAs (lncRNAs) and 94,230 mRNA-miRNA pairs were identified. Hepatospora eriocheir is the primary pathogen of E. sinensis and can cause hepatopancreatic necrosis disease (HPND); the intestine is the main target tissue. Here, we attempted to identify the key genes related to H. eriocheir infection in the intestines of E. sinensis. By combining Iso-Seq and Illumina RNA-seq analysis, we identified a total of 12,708 differentially expressed unigenes (DEUs; 6,696 upregulated and 6,012 downregulated) in the crab intestine following infection with H. eriocheir. Based on the biological analysis of these DEUs, several key processes were identified, including energy metabolism-related pathways, cell apoptosis and innate immune-related pathways. Twelve selected genes from these DEUs were subsequently verified by quantitative real-time PCR (qRT-PCR) analysis. Our findings enhance our understanding of the E. sinensis transcriptome and the specific association between E. sinensis and H. eriocheir infection.

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 gonad transcriptome analysis of Amur sturgeon Dmrt family genes: identification, characterization, and expression patterns during gonadal differentiation

The regulatory mechanisms that govern sex differentiation in sturgeon are still poorly understood. The doublesex and Mab-3-related transcription factor (Dmrt) gene family is known for its extensive roles in sex determination and differentiation across vertebrates. This study aimed to identify new members of sturgeon Dmrt family genes and core actors in the gonadal differentiation of Amur sturgeon. A full-length gonad transcriptome database was exploited to identify Dmrt gene orthologs. Analyses of phylogenetic relationships and selection pressure were performed, and tissue expression profiles and spatiotemporal expression patterns in gonads were then analyzed using real-time PCR. In total, five Dmrt family genes were identified from the full-length gonad transcriptome, including Dmrt2, DmrtA1, DmrtA2, DmrtB1a, and DmrtB1b. Phylogenetic analysis showed that these genes were clustered into clades corresponding to the doublesex/Mav-3 (DM) genes of vertebrates. Furthermore, the analysis of evolutionary selective pressure indicated that DmrtB1a and DmrtB1b were subject to positive selection, suggesting the existence of adaptive evolution in sturgeon. The extensive tissue expression profiling of each Dmrt family gene revealed typical characteristics. Remarkably, according to a spatiotemporal expression pattern analysis, in later stages, DmrtB1b expression increased in testes and was significantly higher in testes than in ovaries at 24 months after hatching (M) and 36 M. This study provides a genetic resource of full-length Dmrt family genes and increases the understanding of Dmrt functions in sex differentiation in sturgeon.

Full-Length Transcriptome Reconstruction Reveals the Genetic Mechanisms of Eyestalk Displacement and Its Potential Implications on the Interspecific Hybrid Crab (Scylla serrata ♀ × S. paramamosain ♂)

Simple Summary

The eyestalk is a key organ in crustaceans that produces neurohormones and regulates a range of physiological functions. Eyestalk displacement was discovered in some first-generation (F1) offspring of the novel interspecific hybrid crab (Scylla serrata ♀ × S. paramamosain ♂). To uncover the genetic mechanism underlying eyestalk displacement and its potential implications, high-quality transcriptome was reconstructed using single-molecule real-time (SMRT) sequencing. A total of 37 significantly differential alternative splicing (DAS) events (17 up-regulated and 20 down-regulated) and 1475 significantly differential expressed transcripts (DETs) (492 up-regulated and 983 down-regulated) were detected in hybrid crabs with displaced eyestalks (DH). The most significant DAS events and DETs were annotated as being endoplasmic reticulum chaperone BiP and leucine-rich repeat protein lrrA-like isoform X2. In addition, the top ten significant gene ontology (GO) terms were related to the cuticle or chitin. Overall, this study highlights the underlying genetic mechanisms of eyestalk displacement and provide useful knowledge for mud crab (Scylla spp.) crossbreeding.

Abstract

The lack of high-quality juvenile crabs is the greatest impediment to the growth of the mud crab (Scylla paramamosain) industry. To obtain high-quality hybrid offspring, a novel hybrid mud crab (S. serrata ♀ × S. paramamosain ♂) was successfully produced in our previous study. Meanwhile, an interesting phenomenon was discovered, that some first-generation (F1) hybrid offspring’s eyestalks were displaced during the crablet stage I. To uncover the genetic mechanism underlying eyestalk displacement and its potential implications, both single-molecule real-time (SMRT) and Illumina RNA sequencing were implemented. Using a two-step collapsing strategy, three high-quality reconstructed transcriptomes were obtained from purebred mud crabs (S. paramamosain) with normal eyestalks (SPA), hybrid crabs with normal eyestalks (NH), and hybrid crabs with displaced eyestalks (DH). In total, 37 significantly differential alternative splicing (DAS) events (17 up-regulated and 20 down-regulated) and 1475 significantly differential expressed transcripts (DETs) (492 up-regulated and 983 down-regulated) were detected in DH. The most significant DAS events and DETs were annotated as being endoplasmic reticulum chaperone BiP and leucine-rich repeat protein lrrA-like isoform X2. In addition, the top ten significant GO terms were related to the cuticle or chitin. Overall, high-quality reconstructed transcriptomes were obtained for the novel interspecific hybrid crab and provided valuable insights into the genetic mechanisms of eyestalk displacement in mud crab (Scylla spp.) crossbreeding.

Extent and complexity of RNA processing in honey bee queen and worker caste development

The distinct honeybee (Apis mellifera) worker and queen castes have become a model for the study of genomic mechanisms of phenotypic plasticity. Here we performed a nanopore-based direct RNA sequencing with exceptionally long reads to compare the mRNA transcripts between queen and workers at three points during their larval development. We found thousands of significantly differentially expressed transcript isoforms (DEIs) between queen and worker larvae. These DEIs were formatted by a flexible splicing system. We showed that poly(A) tails participated in this caste differentiation by negatively regulating the expression of DEIs. Hundreds of isoforms uniquely expressed in either queens or workers during their larval development, and isoforms were expressed at different points in queen and worker larval development demonstrating a dynamic relationship between isoform expression and developmental mechanisms. These findings show the full complexity of RNA processing and transcript expression in honey bee phenotypic plasticity.

•

Honeybee caste differentiation has a complexity of RNA processing

•

Isoforms differentially express between queens and workers during larval development

•

Isoforms are formatted by a flexible alternative splicing system

•

Poly(A) tails are negatively correlated with isoform expression

Full-Length Transcriptome Analysis Provides New Insights Into the Diversity of Immune-Related Genes in Portunus trituberculatus

Generally, invertebrates were thought to solely rely on their non-specific innate immune system to fight against invading microorganisms. However, increasing studies have implied that the innate immune response of invertebrates displayed diversity and specificity owing to the hyper-variable immune molecules in organisms. In order to get an insight into the diversity of immune-related genes in Portunus trituberculatus, a full-length transcriptome analysis of several immune-related tissues (hemocytes, hepatopancreas and gills) in P. trituberculatus was performed and the diversity of several immune-related genes was analyzed. The full-length transcriptome analysis of P. trituberculatus was conducted using a combination of SMRT long-read sequencing and Illumina short-read sequencing. A total of 17,433 nonredundant full-length transcripts with average length of 2,271 bp and N50 length of 2,841 bp were obtained, among which 13,978 (80.18%) transcripts were annotated. Moreover, numerous transcript variants of various immune-related genes were identified, including pattern recognition receptors, antimicrobial peptides, heat shock proteins (HSPs), antioxidant enzymes and vital molecules in prophenoloxidase (proPO)-activating system. Based on the full-length transcriptome analysis, open reading frames (ORFs) of four C-type lectins (CTLs) were cloned, and tissue distributions showed that the four CTLs were ubiquitously expressed in all the tested tissues, and mainly expressed in hepatopancreas and gills. The transcription of the four CTLs significantly increased in several immune-related tissues (hemocytes, hepatopancreas and gills) of P. trituberculatus challenged with Vibrio alginolyticus and displayed different profiles. Moreover, the four CTLs displayed distinct bacterial binding and antibacterial activities. The recombinant protein PtCTL-1 (rPtCTL-1) and rPtCTL-3 displayed bacterial binding and antibacterial activities against all tested bacteria. rPtCTL-2 only showed bacterial binding and antibacterial activities against V. alginolyticus. No obvious bacterial binding or antibacterial activities for PtCTL-4 was observed against the tested bacteria. This study enriches the transcriptomic information on P. trituberculatus and provides new insights into the innate immune system of crustaceans. Additionally, our study provided candidates of antibiotic agents for the prevention and treatment of bacteriosis.

The dynamic immune responses of Mandarin fish (Siniperca chuatsi) to ISKNV in early infection based on full-length transcriptome analysis and weighted gene co-expression network analysis

Mandarin fish (Siniperca chuatsi) been seriously harmed by infectious spleen and kidney necrosis virus (ISKNV) in recent years, but the early immune response mechanism of infection is still unknown. Here, we performed RNA sequencing on the spleens of mandarin fish infected with ISKNV at 0, 12, 24, 48, and 72 h post-infection (hpi) using short-read Illumina RNA sequencing and long-read Pacific Biosciences isoform sequencing to generate a full-length transcriptome. The immune responses of mandarin fish infected with ISKNV at the molecular level were characterized by RNA-seq analysis and weighted gene co-expression network analysis (WGCNA). A total of 26,528 full-length transcript sequences were obtained. There were 2,729 (1,680 up-regulated and 1,112 down-regulated), 1,874 (1,136 up-regulated and 738 down-regulated), 2,032 (1,158 up-regulated and 847 down-regulated), and 4,176 (2,233 up-regulated and 1,943 down-regulated) differentially expressed genes (DEGs) in mandarin fish at 12, 24, 48, and 72 hpi, compared with uninfected fish, respectively. A total of four modules of co-expressed DEGs identified by WGCNA were significantly positively correlated to the four time points after infection, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that the immune-related DEGs in all these modules were mainly enriched in Phagosome, Endocytosis, Herpes simplex infection, and Cytokine-cytokine receptor interaction pathways. Further analysis showed that oher signaling pathways, including CAMs, NOD-like receptor and ER protein processing, Intestinal immune network for IgA production, TLR pathway, and Apoptosis significantly enriched in four modules corresponding to 12, 24, 48, and 72 hpi respectively, had specifically participated in the immune response. Hub genes identified based on the high-degree nodes in the WGCN, including CAM3, IL-8, CCL21, STING, SNX1, PFR and TBK1, and some DEGs such as MHCI, MHCII, TfR, STING, TNF α, TBK1, IRF1, and NF-kB, BCR, IgA and Bcl-XL had involved in dynamic molecular response of mandarin fish to ISKNV infection. In sum, this study provides a set of full-length transcriptome of the spleen tissue of mandarin fish for the first time and revealed a group of immune genes and pathways involved in different temporal responses to ISKNV infection, which has implications for resource conservation and aiding the development of strategies to prevent virus early infection for mandarin fish.

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.

Comprehensive transcriptome characterization of Grus japonensis using PacBio SMRT and Illumina sequencing

The red-crowned crane (Grus japonensis) is an endangered species distributed across southeast Russia, northeast China, Korea, and Japan. Here, we sequenced for the first time the full-length unreferenced transcriptome of red-crowned crane mixed samples using a PacBio Sequel platform. A total of 359,136 circular consensus sequences (CCS) were obtained via clustering to remove redundancy. A total of 303,544 full-length non-chimeric sequences were identified by judging whether CCS contained 5' and 3' adapters, and the poly(A) tail. Eight samples were sequenced using Illumina, and PacBio sequencing data were corrected according to the collected Illumina data to obtain more accurate full-length transcripts. A total of 4,100 long non-coding RNAs, 13,115 simple sequences repeat loci and 29 transcription factor families were identified. The expression of lncRNAs and TFs in pancreas was lowest comparing with other tissues. Many enriched immune-related transmission pathways (MHC and IL receptors) were identified in the spleen. This study will contribute to a better understanding of the gene structure and post-transcriptional regulatory network, and provide references for future studies on red-crowned cranes.

A proteomics investigation of 'immune priming' in Penaeus vannamei as shown by isobaric tags for relative and absolute quantification

Invertebrates are considered completely dependent on their innate immunity to defend themselves against pathogens as they lack an adaptive immunity. However, a growing body of evidence has indicated a specific acquired immunity called 'immune priming' may exist. The Pacific white shrimp, Penaeus vannamei is one of the most economically important shrimp species in the world. In the previous research, we investigated the hepatopancreas immune response of shrimp immunized with trans -vp28 gene Synechocystis sp. PCC6803 at the protein level. In this study, on the basis of the previous research, the shrimp were then challenged with WSSV, and hepatopancreas analyzed using isobaric tags for relative and absolute quantification (i TRAQ) labeling. In total, 308 differentially expressed proteins (DEPs) were identified including 84 upregulated and 224 downregulated. Upregulated proteins such as calmodulin B and calreticulin, and downregulated proteins such as calnexin, and signaling pathways like Ras, mTOR were differentially expressed in both studies. Data from this study are more significant than previous work and indicate increased sensitivity to WSSV after immunization with trans-vp28 gene Synechocystis sp. PCC6803. In addition, selected DEPs (upregulated: A0A3R7QHH6 and downregulated: A0A3R7PEF6, A0A3R7MGX8, A0A423TPJ4, and A0A3R7QCC2) were randomly analyzed using parallel reaction monitoring (PRM). These data preliminarily confirm immune priming in P. vannamei, and show that the initial stimulation with trans -vp28 gene Synechocystis sp. PCC6803 regulate P. vannamei immune responses and they provide shrimp with enhanced immune protection against secondary stimulation.

Full-length transcriptome analysis reveals the mechanism of acupuncture at PC6 improves cardiac function in myocardial ischemia model

Background

The pathological process of myocardial ischemia (MI) is very complicated. Acupuncture at PC6 has been proved to be effective against MI injury, but the mechanism remains unclear. This study investigated the mechanism that underlies the effect of acupuncture on MI through full-length transcriptome.

Methods

Adult male C57/BL6 mice were randomly divided into control, MI, and PC6 groups. Mice in MI and PC6 group generated MI model by ligating the left anterior descending (LAD) coronary artery. The samples were collected 5 days after acupuncture treatment.

Results

The results showed that treatment by acupuncture improved cardiac function, decreased myocardial infraction area, and reduced the levels of cTnT and cTnI. Based on full-length transcriptome sequencing, 5083 differential expression genes (DEGs) and 324 DEGs were identified in the MI group and PC6 group, respectively. These genes regulated by acupuncture were mainly enriched in the inflammatory response pathway. Alternative splicing (AS) is a post-transcriptional action that contributes to the diversity of protein. In all samples, 8237 AS events associated with 1994 genes were found. Some differential AS-involved genes were enriched in the pathway related to heart disease. We also identified 602 new genes, 4 of which may the novel targets of acupuncture in MI.

Conclusions

Our findings suggest that the effect of acupuncture on MI may be based on the multi-level regulation of the transcriptome.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13020-021-00465-8.

Integrative Transcriptomic Analysis Reveals the Immune Mechanism for a CyHV-3-Resistant Common Carp Strain

Anti-disease breeding is becoming the most promising solution to cyprinid herpesvirus-3 (CyHV-3) infection, the major threat to common carp aquaculture. Virus challenging studies suggested that a breeding strain of common carp developed resistance to CyHV-3 infection. This study illustrates the immune mechanisms involved in both sensitivity and anti-virus ability for CyHV3 infection in fish. An integrative analysis of the protein-coding genes and long non-coding RNAs (lncRNAs) using transcriptomic data was performed. Tissues from the head kidney of common carp were extracted at days 0 (the healthy control) and 7 after CyHV-3 infection (the survivors) and used to analyze the transcriptome through both Illumina and PacBio sequencing. Following analysis of the GO terms and KEGG pathways involved, the immune-related terms and pathways were merged. To dig out details on the immune aspect, the DEGs were filtered using the current common carp immune gene library. Immune gene categories and their corresponding genes in different comparison groups were revealed. Also, the immunological Gene Ontology terms for lncRNA modulation were retained. The weighted gene co-expression network analysis was used to reveal the regulation of immune genes by lncRNA. The results demonstrated that the breeding carp strain develops a marked resistance to CyHV-3 infection through a specific innate immune mechanism. The featured biological processes were autophagy, phagocytosis, cytotoxicity, and virus blockage by lectins and MUC3. Moreover, the immune-suppressive signals, such as suppression of IL21R on STAT3, PI3K mediated inhibition of inflammation by dopamine upon infection, as well as the inhibition of NLRC3 on STING during a steady state. Possible susceptible factors for CyHV-3, such as ITGB1, TLR18, and CCL4, were also revealed from the non-breeding strain. The results of this study also suggested that Nramp and PAI regulated by LncRNA could facilitate virus infection and proliferation for infected cells respectively, while T cell leukemia homeobox 3 (TLX3), as well as galectin 3 function by lncRNA, may play a role in the resistance mechanism. Therefore, immune factors that are immunogenetically insensitive or susceptible to CyHV-3 infection have been revealed.

Full-length transcriptome sequencing and identification of immune-related genes in the critically endangered Hucho bleekeri

Hucho bleekeri is a glacial relict and critically endangered fish restricted to the Yangtze River drainage in China. The lack of basic genomic information and immune characteristics will hinder the way toward protecting this species. In the present study, we conducted the first transcriptome analysis of H. bleekeri using the combination of SMRT and Illumina sequencing technology. Transcriptome sequencing generated a total of 93,330 non-redundant full-length unigenes with a mean length of 3072 bp. A total of 92,472 (99.08%) unigenes were annotated in at least one of the Nr protein, Swiss-Prot, KEGG, KOG, GO, Nt and Pfam databases. KEGG analysis showed that a total of 7240 unigenes belonging to 28 immune pathways were annotated to the immune system category. Meanwhile, differentially expressed genes between mucosa-associated tissues (skin, gill and hindgut) and systemic-immune tissues (spleen, head kidney and liver) were obtained. Importantly, genes participating in diverse immune signalling pathways and their expression profiles in H. bleekeri were discussed. In addition, a large number of long non-coding RNAs (lncRNAs) and simple sequence repeats (SSRs) were obtained in the H. bleekeri transcriptome. The present study will provide basic genomic information for H. bleekeri and for further research on analysing the characteristics of both the innate and adaptive immune systems of this critically endangered species.

Global tissue transcriptomic analysis to improve genome annotation and unravel skin pigmentation in goldfish

Goldfish is an ornamental fish with diverse phenotypes. However, the limited genomic resources of goldfish hamper our understanding of the genetic basis for its phenotypic diversity. To provide enriched genomic resources and infer possible mechanisms underlying skin pigmentation, we performed a large-scale transcriptomic sequencing on 13 adult goldfish tissues, larvae at one- and three-days post hatch, and skin tissues with four different color pigmentation. A total of 25.52 Gb and 149.80 Gb clean data were obtained using the PacBio and Illumina platforms, respectively. Onto the goldfish reference genome, we mapped 137,674 non-redundant transcripts, of which 5.54% was known isoforms and 78.53% was novel isoforms of the reference genes, and the remaining 21,926 isoforms are novel isoforms of additional new genes. Both skin-specific and color-specific transcriptomic analyses showed that several significantly enriched genes were known to be involved in melanogenesis, tyrosine metabolism, PPAR signaling pathway, folate biosynthesis metabolism and so on. Thirteen differentially expressed genes across different color skins were associated with melanogenesis and pteridine synthesis including mitf, ednrb, mc1r, tyr, mlph and gch1, and xanthophore differentiation such as pax7, slc2a11 and slc2a15. These transcriptomic data revealed pathways involved in goldfish pigmentation and improved the gene annotation of the reference genome.

Full-length transcriptome and long non-coding RNA profiling of whiteleg shrimp Penaeus vannamei hemocytes in response to Spiroplasma eriocheiris infection

Spiroplasma eriocheiris (S. eriocheiris) infection causes a significant economic loss in Penaeus vannamei (P. vannamei) culture industry. However, the response of P. vannamei hemocytes to S. eriocheiris infection has not been extensively studied. In this study, we conducted full-length transcriptome and long non-coding RNA (lncRNA) analyses of P. vannamei hemocytes by a challenge test with S. eriocheiris. Following assembly and annotation, there were 8077 high-quality unigenes. A total of 1168 differentially expressed genes (DEGs) were obtained, including 792 up-regulated and 376 down-regulated genes by differential expression analysis. Gene ontology (GO) enrichment analysis showed that the up-regulated DEGs were mainly clustered into immune system process, defense response, cell cycle and organelle organization. On the other hand, the down-regulated DEGs included that genes that were mainly clustered into metabolic processes related to organic compounds, metabolic process and cellular metabolic process. Protein-protein interaction (PPI) network analysis of DEGs indicated that the pivotal gene interactions were connected to stress response, immune system process and cell cycle. The lncRNA analysis identified multiple lncRNAs, which were highly co-expressed with the immune-related genes, such as lncRNA transcript-12631 and transcript-12631, suggesting that lncRNAs may be involved in the regulation of immune defense in shrimp hemocytes. Additionally, 20 hub unigenes and putative lncRNAs related to immune system were validated by quantitative real-time PCR (qRT-PCR), validating the reliability of RNA-Seq. This study revealed a close connection between the immune and metabolic systems of S. eriocheiris infected P. vannamei.

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

As invertebrates, shrimp are generally thought to solely rely on their innate immune system to combat invading pathogens. Recently, an increasing number of studies have revealed that the innate immune response of invertebrates exhibits diversity and specificity based on their diverse immune molecules. Herein, a full-length transcriptome analysis of several immune-related tissues (hepatopancreas, gill, hemocytes, stomach and intestine) in the kuruma shrimp (Marsupenaeus japonicus) was conducted to identify immune-related molecules with a focus on transcript variations. In total, 11,222 nonredundant full-length transcripts with an N50 length of 5174 were obtained, and most of these transcripts (94.84%) were successfully annotated. In addition, a total of 147 long noncoding RNAs (lncRNAs) were also predicted. Importantly, transcript variants of several vital immune-related genes were observed, including twenty-five alpha-2-macroglobulins (α2-Ms), ten Toll-like receptors (TLRs), six C-type lectins (CTLs), five M-type lectins (MTLs) and three Down syndrome cell adhesion molecules (Dscams). Furthermore, 509 nonredundant full-length transcripts were predicted to be generated from alternative splicing (AS) events, which contribute to the diversity of immune molecules. Overall, our study provides valuable data on the full-length transcripts of M. japonicus, which will facilitate the exploration of immune molecules in this species. Moreover, numerous transcript variants of immune molecules detected in this study provide clues for further investigating the diversity and specificity of the innate immune response in shrimp.

Full Transcriptome Analysis of Callus Suspension Culture System of Bletilla striata

Background

Bletilla striata has been widely used in the pharmacology industry. To effectively produce the secondary metabolites through suspension cultured cells of B. striata, it is important to exploring the full-length transcriptome data and the genes related to cell growth and chemical producing of all culture stages. We applied a combination of Real-Time Sequencing of Single Molecule (SMRT) and second-generation sequencing (SGS) to generate the complete and full-length transcriptome of B. striata suspension cultured cells.

Methods

The B. striata transcriptome was formed in de novo way by using PacBio isoform sequencing (Iso-Seq) on a pooled RNA sample derived from 23 samples of 10 culture stages, to explore the potential for capturing full-length transcript isoforms. All unigenes were obtained after splicing, assembling, and clustering, and corrected by the SGS results. The obtained unigenes were compared with the databases, and the functions were annotated and classified.

Results and conclusions

A total of 100,276 high-quality full-length transcripts were obtained, with an average length of 2530 bp and an N50 of 3302 bp. About 52% of total sequences were annotated against the Gene Ontology, 53,316 unigenes were hit by KOG annotations and divided into 26 functional categories, 80,020 unigenes were mapped by KEGG annotations and clustered into 363 pathways. Furthermore, 15,133 long-chain non-coding RNAs (lncRNAs) were detected. And 68,996 coding sequences were identified based on SSR analysis, among which 31 pairs of primers selected at random were amplified and obtained stable bands. In conclusion, our results provide new full-length transcriptome data and genetic resources for identifying growth and metabolism-related genes, which provide a solid foundation for further research on its growth regulation mechanisms and genetic engineering breeding mechanisms of B. striata.

Reconstruction and functional annotation of Ascosphaera apis full-length transcriptome utilizing PacBio long reads combined with Illumina short reads

Ascosphaera apis is a widespread fungal pathogen of honeybee larvae that results in chalkbrood disease, leading to heavy losses for the beekeeping industry in China and many other countries. This work was aimed at generating a full-length transcriptome of A. apis using PacBio single-molecule real-time (SMRT) sequencing. Here, more than 23.97 Gb of clean reads was generated from long-read sequencing of A. apis mycelia, including 464,043 circular consensus sequences (CCS) and 394,142 full-length non-chimeric (FLNC) reads. In total, we identified 174,095 high-confidence transcripts covering 5141 known genes with an average length of 2728 bp. We also discovered 2405 genic loci and 11,623 isoforms that have not been annotated yet within the current reference genome. Additionally, 16,049, 10,682, 4520 and 7253 of the discovered transcripts have annotations in the Non-redundant protein (Nr), Clusters of Eukaryotic Orthologous Groups (KOG), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Moreover, 1205 long non-coding RNAs (lncRNAs) were identified, which have less exons, shorter exon and intron lengths, shorter transcript lengths, lower GC percent, lower expression levels, and fewer alternative splicing (AS) evens, compared with protein-coding transcripts. A total of 253 members from 17 transcription factor (TF) families were identified from our transcript datasets. Finally, the expression of A. apis isoforms was validated using a molecular approach. Overall, this is the first report of a full-length transcriptome of entomogenous fungi including A. apis. Our data offer a comprehensive set of reference transcripts and hence contributes to improving the genome annotation and transcriptomic study of A. apis.

Full-length transcriptome sequencing combined with RNA-seq analysis revealed the immune response of fat greenling (Hexagrammos otakii) to Vibrio harveyi in early infection

Fat greenling (Hexagrammos otakii) is an important commercial marine fish species cultured in northeast Asia, but its available gene sequences are limited. Vibrio harveyi is a causative agent of vibriosis in fat greenling and also causes severe losses to the aquaculture industry in China. In order to obtain more high-quality transcript information and investigate the early immune response of fat greenling against V. harveyi, the fish were artificially infected with V. harveyi, and five sampling points were set within 48 h. Iso-Seq combined with RNA-Seq were applied in the comprehensive transcriptome analysis of V. harveyi-infected fat greenling. Total 42,225 consensus isoforms were successfully extracted from the result of Iso-Seq, and more than 19,000 ORFs were predicted. In addition, total three modules were identified by WGCNA which significantly positive correlated to the infection time, and the KEGG analysis showed that the immune-related genes in these modules mainly enriched in TLR signaling pathway, NF-κB signaling pathway and Endocytosis. The activation of inflammation and endocytosis was the most significant characteristics of fat greenling immune response during the early infection. Based on the WGCNA, a series of high-degree nodes in the networks were identified as hub genes. The protein structures of cold-inducible RNA-binding protein (CIRBP), poly [ADP-ribose] polymerase 1 (PARP1) and protein arginine N-methyl transferase 1 (PRMT1) were subsequently found to be highly conserved in vertebrate, and the gene expression pattern of CIRBP, PARP1, PRMT1 and a part of TLR/NF-κB pathway-related genes indicated that these proteins might have similar biological functions in regulation of inflammatory response in teleost fish. The results of this study provided the first systematical full-length transcriptome profile of fat greenling and characterized its immune responses in early infection of V. harvey, which will serve as the foundation for further exploring the molecular mechanism of immune defense against bacterial infection in fat greenling.

Transcriptome analysis and immune-related genes expression reveals the immune responses of Macrobrachium rosenbergii infected by Enterobacter cloacae

Macrobrachium rosenbergii is an important cultural species in China and other Southeast Asian countries. However, Enterobacter cloacae infection has caused a great economic loss in M. rosenbergii culture industry. The immune responses of M. rosenbergii to the E. cloacae infection is not fully characterized. To investigate the immune response of M. rosenbergii against E. cloacae, we performed transcriptome analysis of the M. rosenbergii hepatopancreas with and without E. cloacae infection using RNA-seq. After assembly and annotation, 29,731 high quality unigenes were obtained from RNA-seq data. Differential expression analysis revealed the existence of 2498 significantly differently expressed genes (DEGs) at 12 h post infection, with 1365 up-regulated and 1133 down-regulated genes. Among these DEGs, some well-known immune-related genes were up-regulated significantly, including C-type lectin 1, lectin 3, anti-lipopolysaccharide factor 2, Cu/Zn superoxide dismutase and heat shock protein 70. GO analysis demonstrated 24 biological process subcategories, 14 cellular component subcategories, and 12 molecular function subcategories that were enriched among these DEGs, and some DEGs were clustered into immune related subcategories such as immune system process, response to stimulus, biological adhesion, and antioxidant activity. These DEGs were enriched into 216 KEGG pathways including a core set of immune correlated pathways notably in phagosome and lysosome. In addition, 5 up-regulated and 5 down-regulated immune-related DEGs were selected for further validation by quantitative real-time PCR and the results showed consistence with the RNA-seq data. Additionally, the expression level of six selected immune-related genes (ALF2, CLEC1, LEC3, hemocyanin1, HSP70 and SOD) based on the transcriptomic data were monitored at different point of time in hepatopancreas, gill, hemolymph and intestine. Results revealed these immune-related genes were significantly up-regulated in different tissues from 6 to 24 h after E. cloacae infection. Overall, these results provided valuable information for further studying the immune response of M. rosenbergii against E. cloacae infection.

Full-length transcriptome sequencing and comparative transcriptomic analysis to uncover genes involved in early gametogenesis in the gonads of Amur sturgeon (Acipenser schrenckii)

BACKGROUND

Sturgeons (Acipenseriformes) are polyploid chondrostean fish that constitute an important model species for studying development and evolution in vertebrates. To better understand the mechanisms of reproduction regulation in sturgeon, this study combined PacBio isoform sequencing (Iso-Seq) with Illumina short-read RNA-seq methods to discover full-length genes involved in early gametogenesis of the Amur sturgeon, Acipenser schrenckii.

RESULTS

A total of 50.04 G subread bases were generated from two SMRT cells, and herein 164,618 nonredundant full-length transcripts (unigenes) were produced with an average length of 2782 bp from gonad tissues (three testes and four ovaries) from seven 3-year-old A. schrenckii individuals. The number of ovary-specific expressed unigenes was greater than those of testis (19,716 vs. 3028), and completely different KEGG pathways were significantly enriched between the ovary-biased and testis-biased DEUs. Importantly, 60 early gametogenesis-related genes (involving 755 unigenes) were successfully identified, and exactly 50% (30/60) genes of those showed significantly differential expression in testes and ovaries. Among these, the Amh and Gsdf with testis-biased expression, and the Foxl2 and Cyp19a with ovary-biased expression strongly suggested the important regulatory roles in spermatogenesis and oogenesis of A. schrenckii, respectively. We also found the four novel Sox9 transcript variants, which increase the numbers of regulatory genes and imply function complexity in early gametogenesis. Finally, a total of 236,672 AS events (involving 36,522 unigenes) were detected, and 10,556 putative long noncoding RNAs (lncRNAs) and 4339 predicted transcript factors (TFs) were also respectively identified, which were all significantly associated with the early gametogenesis of A. schrenckii.

CONCLUSIONS

Overall, our results provide new genetic resources of full-length transcription data and information as a genomic-level reference for sturgeon. Crucially, we explored the comprehensive genetic characteristics that differ between the testes and ovaries of A. schrenckii in the early gametogenesis stage, which could provide candidate genes and theoretical basis for further the mechanisms of reproduction regulation of sturgeon.

Full-length transcriptome analysis of Coptis deltoidea and identification of putative genes involved in benzylisoquinoline alkaloids biosynthesis based on combined sequencing platforms

The study carry out comprehensive transcriptome analysis of C. deltoidea and exploration of BIAs biosynthesis and accumulation based on UHPLC-MS/MS and combined sequencing platforms. Coptis deltoidea is an important medicinal plant with a long history of medicinal use, which is rich in benzylisoquinoline alkaloids (BIAs). In this study, Ultra performance liquid chromatography-electrospray ionization tandem mass spectrometry (UHPLC-ESI-MS/MS) and combined sequencing platforms were performed for exploration of BIAs biosynthesis, accumulation and comprehensive transcriptome analysis of C. deltoidea. By metabolism profiling, the accumulation of ten BIAs was analyzed using UHPLC-MS/MS and different contents were observed in different organs. From transcriptome sequencing result, we applied single-molecule real-time (SMRT) sequencing to C. deltoidea and generated a total of 75,438 full-length transcripts. We proposed the candidate biosynthetic pathway of tyrosine, precursor of BIAs, and identified 64 full length-transcripts encoding enzymes putatively involved in BIAs biosynthesis. RNA-Seq data indicated that the majority of genes exhibited relatively high expression level in roots. Transport of BIAs was also important for their accumulation. Here, 9 ABC transporters and 2 MATE transporters highly homologous to known alkaloid transporters related with BIAs transport in roots and rhizomes were identified. These findings based on the combined sequencing platforms provide valuable genetic information for C. deltoidea and the results of transcriptome combined with metabolome analysis can help us better understand BIAs biosynthesis and transport in this medicinal plant. The information will be critical for further characterization of C. deltoidea transcriptome and molecular-assisted breeding for this medicinal plant with scarce resources.

Characterization of Full-Length Transcriptome Sequences and Splice Variants of Lateolabrax maculatus by Single-Molecule Long-Read Sequencing and Their Involvement in Salinity Regulation

Transcriptome complexity plays crucial roles in regulating the biological functions of eukaryotes. Except for functional genes, alternative splicing and fusion transcripts produce a vast expansion of transcriptome diversity. In this study, we applied PacBio single-molecule long-read sequencing technology to unveil the whole transcriptome landscape of Lateolabrax maculatus. We obtained 28,809 high-quality non-redundant transcripts, including 18,280 novel isoforms covering 8,961 annotated gene loci within the current reference genome and 3,172 novel isoforms. A total of 10,249 AS events were detected, and intron retention was the predominant AS event. In addition, 1,359 alternative polyadenylation events, 3,112 lncRNAs, 29,609 SSRs, 365 fusion transcripts, and 1,194 transcription factors were identified in this study. Furthermore, we performed RNA-Seq analysis combined with Iso-Seq results to investigate salinity regulation mechanism at the transcripts level. A total of 518 transcripts were differentially expressed, which were further divided into 8 functional groups. Notably, transcripts from the same genes exhibited similar or opposite expression patterns. Our study provides a comprehensive view of the transcriptome complexity in L. maculatus, which significantly improves current gene models. Moreover, the diversity of the expression patterns of transcripts may enhance the understanding of salinity regulatory mechanism in L. maculatus and other euryhaline teleosts.

Conjoint Analysis of SMRT- and Illumina-Based RNA-Sequencing Data of Fenneropenaeus chinensis Provides Insight Into Sex-Biased Expression Genes Involved in Sexual Dimorphism

Fenneropenaeus chinensis (F. chinensis) is one of the most commercially important cultured shrimps in China. The adult F. chinensis exhibit sexual dimorphism in growth and body color. In this research, we profiled the whole transcriptome of F. chinensis by using single molecule real-time-based full-length transcriptome sequencing. We further performed Illumina-based short reads RNA-seq on muscle and gonad of two sexes to detect the sex-biased expression genes. In muscle, we observed significantly more female-biased transcripts. With the differentially expressed transcripts (DETs) in muscle, some pathways related to the energy metabolism were enriched, which may be responsible for the difference of growth. We also digged out a pathway named porphyrin and chlorophyll metabolism. It was speculated to relevant to the difference of body color between the two sexes of shrimp. Interestingly, almost all DETs in these pathways were female-biased expression in muscle, which could explain the phenomenon of better growth performance and darker body color in female. In gonad, several pathways involved in reproduction were enriched. For instance, some female-biased DETs participated in the arachidonic acid metabolism, which was reported crucial in female reproduction. In conclusion, our studies identified abundant sex-biased expression transcripts and important pathways involved in sexual dimorphism by using the RNA-seq method. It provided a basis for future researches on the sexual dimorphism of F. chinensis.

Comparative proteomic investigation of Marsupenaeus japonicus hepatopancreas challenged with Vibrio parahaemolyticus and white spot syndrome virus

This study aimed to use isobaric tags (IBTs) to investigate the immune response of the hepatopancreas of Marsupenaeus japonicas infected with Vibrio parahaemolyticus or white spot syndrome virus (WSSV). Liquid chromatography-tandem mass spectrometry and protein sequencing identified 1005 proteins. Among them, 109 proteins were upregulated and 94 were downregulated after V. parahaemolyticus infection. After WSSV infection, 130 proteins were identified as differentially abundant, including 88 that were upregulated and 42 were downregulated. Fifty-four proteins were identified as differentially abundant after both V. parahaemolyticus and WSSV infection. A number of proteins related to cytoskeletal processes, including actin and myosin, and apoptosis-related proteins were upregulated in shrimp after V. parahaemolyticus and WSSV infection, indicating that phagocytosis and apoptosis may be involved in the response to in V. parahaemolyticus or WSSV infection. Quantitative real-time PCR was carried out to verify the reliability of the proteomic data. These data provide a basis to characterize the immunity-related processes of shrimp in response to infection with WSSV or V. parahaemolyticus.