Long non-coding RNA profile in banana shrimp, Fenneropenaeus merguiensis and the potential role of lncPV13 in vitellogenesis

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.

bestDEG: a web-based application automatically combines various tools to precisely predict differentially expressed genes (DEGs) from RNA-Seq data

Background

Differential gene expression analysis using RNA sequencing technology (RNA-Seq) has become the most popular technique in transcriptome research. Although many R packages have been developed to analyze differentially expressed genes (DEGs), several evaluations have shown that no single DEG analysis method outperforms all others. The validity of DEG identification could be increased by using multiple methods and producing the consensus results. However, DEG analysis methods are complex and most of them require prior knowledge of a programming language or command-line shell. Users who do not have this knowledge need to invest time and effort to acquire it.

Methods

We developed a novel web application called "bestDEG" to automatically analyze DEGs with different tools and compare the results. A differential expression (DE) analysis pipeline was created combining the edgeR, DESeq2, NOISeq, and EBSeq packages; selected because they use different statistical methods to identify DEGs. bestDEG was evaluated on human datasets from the MicroArray Quality Control (MAQC) project.

Results

The performance of the bestDEG web application with the human datasets showed excellent results, and the consensus method outperformed the other DE analysis methods in terms of precision (94.71%) and specificity (97.01%). bestDEG is a rapid and efficient tool to analyze DEGs. With bestDEG, users can select DE analysis methods and parameters in the user-friendly web interface. bestDEG also provides a Venn diagram and a table of results. Moreover, the consensus method of this tool can maximize the precision or minimize the false discovery rate (FDR), which reduces the cost of gene expression validation by minimizing wet-lab experiments.

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.

Identification, characterization and functional analysis of gonadal long noncoding RNAs in a protogynous hermaphroditic teleost fish, the ricefield eel (Monopterus albus)

Background

An increasing number of long noncoding RNAs (lncRNAs) have been found to play important roles in sex differentiation and gonad development by regulating gene expression at the epigenetic, transcriptional and posttranscriptional levels. The ricefield eel, Monopterus albus, is a protogynous hermaphroditic fish that undergoes a sequential sex change from female to male. However, the roles of lncRNA in the sex change is unclear.

Results

Herein, we performed RNA sequencing to analyse lncRNA expression patterns in five different stages of M. albus development to investigate the roles of lncRNAs in the sex change process. A total of 12,746 lncRNAs (1503 known lncRNAs and 11,243 new lncRNAs) and 2901 differentially expressed lncRNAs (DE-lncRNAs) were identified in the gonads. The target genes of the DE-lncRNAs included foxo1, foxm1, smad3, foxr1, camk4, ar and tgfb3, which were mainly enriched in signalling pathways related to gonadal development, such as the insulin signalling pathway, MAPK signalling pathway, and calcium signalling pathway. We selected 5 highly expressed DE-lncRNAs (LOC109952131, LOC109953466, LOC109954337, LOC109954360 and LOC109958454) for full length amplification and expression pattern verification. They were all expressed at higher levels in ovaries and intersex gonads than in testes, and exhibited specific time-dependent expression in ovarian tissue incubated with follicle-stimulating hormone (FSH) and human chorionic gonadotropin (hCG). The results of quantitative real-time PCR (qRT-PCR) analysis and a dual-luciferase assay showed that znf207, as the gene targeted by LOC109958454, was expressed in multiple tissues and gonadal developmental stages of M. albus, and its expression was also inhibited by the hormones FSH and hCG.

Conclusions

These results provide new insights into the role of lncRNAs in gonad development, especially regarding natural sex changes in fish, which will be useful for enhancing our understanding of sequential hermaphroditism and sex changes in the ricefield eel (M. albus) and other teleosts.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08679-2.