Single molecule real-time sequencing and RNA-seq unravel the role of long non-coding and circular RNA in the regulatory network during Nile tilapia (Oreochromis niloticus) infection with Streptococcus agalactiae

Decoding the fish genome opens a new era in important trait research and molecular breeding in China

Aquaculture represents the fastest-growing global food production sector, as it has become an essential component of the global food supply. China has the world's largest aquaculture industry in terms of production volume. However, the sustainable development of fish culture is hindered by several concerns, including germplasm degradation and disease outbreaks. The practice of genomic breeding, which relies heavily on genome information and genotypephenotype relationships, has significant potential for increasing the efficiency of aquaculture production. In 2014, the completion of the genome sequencing and annotation of the Chinese tongue sole signified the beginning of the fish genomics era in China. Since then, domestic researchers have made dramatic progress in functional genomic studies. To date, the genomes of more than 60 species of fish in China have been assembled and annotated. Based on these reference genomes, evolutionary, comparative, and functional genomic studies have revolutionized our understanding of a wide range of biologically and economically important traits of fishes, including growth and development, sex determination, disease resistance, metamorphosis, and pigmentation. Furthermore, genomic tools and breeding techniques such as SNP arrays, genomic selection, and genome editing have greatly accelerated genetic improvement through the incorporation of functional genomic information into breeding activities. This review aims to summarize the current status, advances, and perspectives of the genome resources, genomic study of important traits, and genomic breeding techniques of fish in China. The review will provide aquaculture researchers, fish breeders, and farmers with updated information concerning fish genomic research and breeding technology. The summary will help to promote the genetic improvement of production traits and thus will support the sustainable development of fish aquaculture.

Genome-Wide Identification and Involvement in Response to Biotic and Abiotic Stresses of lncRNAs in Turbot (Scophthalmus maximus)

Long non-coding RNAs (lncRNAs) play crucial roles in a variety of biological processes, including stress response. However, the number, characteristics and stress-related expression of lncRNAs in turbot are still largely unknown. In this study, a total of 12,999 lncRNAs were identified at the genome-wide level of turbot for the first time using 24 RNA-seq datasets. Sequence characteristic analyses of transcripts showed that lncRNA transcripts were shorter in average length, lower in average GC content and in average expression level as compared to the coding genes. Expression pattern analyses of lncRNAs in 12 distinct tissues showed that lncRNAs, especially lincRNA, exhibited stronger tissue-specific expression than coding genes. Moreover, 612, 1351, 1060, 875, 420 and 1689 differentially expressed (DE) lncRNAs under Vibrio anguillarum, Enteromyxum scophthalmi, and Megalocytivirus infection and heat, oxygen, and salinity stress conditions were identified, respectively. Among them, 151 and 62 lncRNAs showed differential expression under various abiotic and biotic stresses, respectively, and 11 lncRNAs differentially expressed under both abiotic and biotic stresses were selected as comprehensive stress-responsive lncRNA candidates. Furthermore, expression pattern analysis and qPCR validation both verified the comprehensive stress-responsive functions of these 11 lncRNAs. In addition, 497 significantly co-expressed target genes (correlation coefficient (R) > 0.7 and q-value < 0.05) for these 11 comprehensive stress-responsive lncRNA candidates were identified. Finally, GO and KEGG enrichment analyses indicated that these target genes were enriched mainly in molecular function, such as cytokine activity and active transmembrane transporter activity, in biological processes, such as response to stimulus and immune response, and in pathways, such as protein families: signaling and cellular processes, transporters and metabolism. These findings not only provide valuable reference resources for further research on the molecular basis and function of lncRNAs in turbot but also help to accelerate the progress of molecularly selective breeding of stress-resistant turbot strains or varieties.

Transcriptome analysis reveals the neuroactive receptor genes response to Streptococcus agalactiae infection in tilapia, Oreochromis niloticus

The detailed crosstalk between the neuroendocrine and immune systems in Oreochromis niloticus, an economically important fish, in response to pathogenic infections, remains unclear. This study revealed the head kidney transcriptional profiles of O. niloticus upon infections with Streptococcus agalactiae, a prevalent pathogen known to cause severe meningitis. Twelve cDNA libraries of O. niloticus head kidney, representing four treatment time points (0, 6, 24, and 48 h), were constructed and a total of 2,528 differentially expressed genes were identified based on pairwise comparisons. KEGG pathway analysis revealed a significant enrichment of the 'neuroactive ligand-receptor interaction' pathway (ko04080), with 13 genes exhibiting differential expression during S. agalactiae infection. Among these, six neuroactive receptor genes (lepr, nr3c1, ptger4, thrb, tspo, and β2-ar) were selected, cloned, and characterized. Although these genes are ubiquitously expressed, and in head kidney leukocytes, their expression was mainly observed in T cells, Mo/Mφ, and NCCs, which are characterized by antimicrobial responses. Furthermore, we examined the response patterns of these six neuroactive receptor genes to gram-positive (S. agalactiae) and gram-negative (Aeromonas hydrophila) bacteria in four different tissues. Notably, lepr, ptger4, tspo, and β2-ar were upregulated in all selected tissues in response to S. agalactiae and A. hydrophila infections. However, nr3c1 and thrb were downregulated in response to S. agalactiae infection in the head kidney and spleen, whereas nr3c1 was upregulated, and thrb was unresponsive to A. hydrophila infection. Our findings provide a theoretical foundation for understanding new links between the neuroendocrine and immune systems during bacterial infection in teleost fish.

Applying genetic technologies to combat infectious diseases in aquaculture

Disease and parasitism cause major welfare, environmental and economic concerns for global aquaculture. In this review, we examine the status and potential of technologies that exploit genetic variation in host resistance to tackle this problem. We argue that there is an urgent need to improve understanding of the genetic mechanisms involved, leading to the development of tools that can be applied to boost host resistance and reduce the disease burden. We draw on two pressing global disease problems as case studies-sea lice infestations in salmonids and white spot syndrome in shrimp. We review how the latest genetic technologies can be capitalised upon to determine the mechanisms underlying inter- and intra-species variation in pathogen/parasite resistance, and how the derived knowledge could be applied to boost disease resistance using selective breeding, gene editing and/or with targeted feed treatments and vaccines. Gene editing brings novel opportunities, but also implementation and dissemination challenges, and necessitates new protocols to integrate the technology into aquaculture breeding programmes. There is also an ongoing need to minimise risks of disease agents evolving to overcome genetic improvements to host resistance, and insights from epidemiological and evolutionary models of pathogen infestation in wild and cultured host populations are explored. Ethical issues around the different approaches for achieving genetic resistance are discussed. Application of genetic technologies and approaches has potential to improve fundamental knowledge of mechanisms affecting genetic resistance and provide effective pathways for implementation that could lead to more resistant aquaculture stocks, transforming global aquaculture.

The immune-related circRNA-miRNA-mRNA ceRNA regulatory network in the liver of turbot (Scophthalmus maximus L.) induced by Vibrio anguillarum

Recently, Vibrio anguillarum, a Gram-negative pathogenic bacterium, has been becoming a major constraint on the development of the turbot aquaculture industry because of its characteristics of worldwide distribution, broad host range and potentially devastating impacts. Although the functions of protein-coding mRNAs in the immune response against bacterial infection have been reported, as well as several non-coding RNAs (ncRNAs), such as circular RNAs (circRNAs) and microRNAs (miRNAs), the relationships between mRNAs and ncRNAs in the immune system of turbot liver are still limited during bacterial infection. In present study, the comprehensive analyses of whole-transcriptome sequencing were conducted in turbot liver infected by V. anguillarum. The differential expression was analyzed in the data of circRNAs, miRNAs, and mRNAs. The interactions of miRNA-circRNA pairs and miRNA-mRNA pairs were predicted basing on the negative regulatory relationships between miRNAs and their target circRNAs\mRNAs. The circRNA-related ceRNA regulatory networks were constructed for the analyses of regulated mechanism in turbot immune system. Subsequently, the RT-qPCR was carried out to verify the results of sequencing. Finally, we identified 31 circRNAs, 53 miRNAs and 948 mRNAs with differential expression. Gene set enrichment analyses using Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways showed that innate immunity was principally activated at the early stages of infection, while adaptive immunity was activated after 24 h. Finally, 65 circRNA-miRNA-mRNA pathways were constructed, based on the hypothesis of ceRNA regulatory networks. In conclusion, our findings provide new insights on the underlying immune response to bacterial infection and identify novel target genes for the prevention and control of disease in turbot.

The potential regulatory role of the lncRNA-miRNA-mRNA axis in teleost fish

Research over the past two decades has confirmed that noncoding RNAs (ncRNAs), which are abundant in cells from yeast to vertebrates, are no longer "junk" transcripts but functional regulators that can mediate various cellular and physiological processes. The dysregulation of ncRNAs is closely related to the imbalance of cellular homeostasis and the occurrence and development of various diseases. In mammals, ncRNAs, such as long noncoding RNAs (lncRNAs) and microRNAs (miRNAs), have been shown to serve as biomarkers and intervention targets in growth, development, immunity, and disease progression. The regulatory functions of lncRNAs on gene expression are usually mediated by crosstalk with miRNAs. The most predominant mode of lncRNA-miRNA crosstalk is the lncRNA-miRNA-mRNA axis, in which lncRNAs act as competing endogenous RNAs (ceRNAs). Compared to mammals, little attention has been given to the role and mechanism of the lncRNA-miRNA-mRNA axis in teleost species. In this review, we provide current knowledge about the teleost lncRNA-miRNA-mRNA axis, focusing on its physiological and pathological regulation in growth and development, reproduction, skeletal muscle, immunity to bacterial and viral infections, and other stress-related immune responses. Herein, we also explored the potential application of the lncRNA-miRNA-mRNA axis in the aquaculture industry. These findings contribute to an enhanced understanding of ncRNA and ncRNA-ncRNA crosstalk in fish biology to improve aquaculture productivity, fish health and quality.

Long non‐coding RNA, a supreme post‐transcriptional immune regulator of bacterial or virus‐driven immune evolution in teleost

The global aquaculture boom, fuelled by a reduction in wild population and detection of novel viruses, has created a demanding market, hence, there is a pressing need to investigate the immune system of fish, further. As the most diverse community of vertebrates and a central contributor to the progressing global aquaculture market, teleost continues to draw vast scientific interest. Recent breakthroughs in multi‐omics technologies have provided a platform to understand the role of long non‐coding RNA (lncRNA) in the host immune system during infection. Emerging evidence shows that teleost lncRNA might have a regulatory role in immune responses, mostly through lncRNA–microRNA (miRNA) sponging. Teleost lncRNA shares a functionally active short sequence complement to target the miRNA which is conserved among the several fish species. Recent report suggests that rhabdovirus exploits a lncRNA in teleost and, to dodge the host immune mechanism and negatively regulate the immune system. This observation reveals the essentiality of lncRNA in pathogen‐driven immunity in teleost. Reports available on the function of teleost lncRNA are still in early stages and experimental verifications are a limiting factor. Unravelling the lncRNA‐mediated immune regulation in fishes could be used against the invading pathogens to strengthen the aquaculture production. This review elaborates on the experimentally identified and functionally characterized lncRNA and its regulatory role in the teleost immune response during infection and pathogen‐driven host immune evolution, which could eventually lead to achieving high standards in aquaculture productivity.

Systematic identification and characterization of lncRNAs and lncRNA-miRNA-mRNA networks in the liver of turbot (Scophthalmus maximus L.) induced with Vibrio anguillarum

Long noncoding RNAs (lncRNAs), can regulate mRNA by targeting miRNA in a competing endogenous RNA network, have become a hot topic in the research of fish immune mechanism recent years. While in turbot (Scophthalmus maximus L.), an economically important marine fish, there are limited researches about the role of lncRNAs in its immune response to bacterial infection. In this study, a total of 184 differentially expressed lncRNAs (DElncRNAs) were systematically identified and characterized using whole-transcriptome sequencing of the liver of turbot challenged with Vibrioanguillarum at 0 h (control) and three different time points post infection (2 h, 12 h and 24 h, respectively). Subsequently, GO and KEGG signaling pathways of differentially expressed lncRNAs were analyzed to predict their function. We found that lncRNAs in our results were significantly enriched in several immune-related signaling pathways, including the NOD-like receptor signaling pathway, Toll-like receptor signaling pathway, Cytokine-cytokine receptor, MAPK signaling pathway, phagosome, PPAR signaling pathway and the regulation of autophagy. In addition, a total of 492 DE lncRNA - DE miRNA -DE mRNA networks were identified at three different time points post infection, which were consisted of 102 networks at 2 h, 122 networks at 12 h and 81 networks at 24 h post infection, respectively. Noticeably, 92 of these regulated networks were immune-related. These observations suggested that lncRNAs can regulate the expression of immune-related genes in the response to bacterial infection in turbot. Moreover, our findings would provide a new insight into the immune response of turbot to pathogen infection and lay a foundation for future study.

Identification of mRNA-miRNA-lncRNA regulatory network associated with the immune response to Aeromonas salmonicides infection in the black rockfish (Sebastes schlegelii)

In aquaculture, Aeromonas salmonicides (A. salmonicida) is a main fish pathogen because of its nearly worldwide distribution, and broad host range. Recently, an increasing number of evidences have uncovered the roles of mRNA-miRNA-lncRNA network in fish diseases. In current study, RNA-seq was conducted in the black rockfish spleen following A. salmonicida infection at 0 h (Sp0 or control) and three different post-infection time-points (2 h: Sp2, 12 h: Sp12 and 24 h: Sp24, respectively) to comprehensively identify differentially expressed (DE) mRNAs, miRNAs and lncRNAs. Enrichment analysis and protein-protein interaction (PPI) analysis of DE mRNAs were performed. Then, expression and correlation analysis for mRNAs and their upstream miRNAs and lncRNAs were conducted. Finally, a total of 1364 mRNAs, 17 miRNAs and 1584 lncRNAs exhibited significantly differential expressions during bacterial infection in the black rockfish spleen. Functional enrichment analysis suggested that they were significantly enriched in several immune-related pathways, including Amino sugar and nucleotide sugar metabolism, Cell adhesion molecules (CAMs), Neuroactive ligand-receptor interaction, Nicotinate and nicotinamide metabolism, Pentose and glucuronate interconversions, Phagosome, Proteasome, etc. Subsequently, 1091 lncRNA-miRNA-mRNA pathways (323 in Sp2, 609 in Sp12 and 207 in Sp24) were constructed including 400 lncRNAs, 69 miRNAs, and 70 mRNAs. Meanwhile, NLRC3/novel-264/LNC_00116154 pathway demonstrated important immune modulating function in the black rockfish against A. salmonicida infection. Finally, the novel mRNA-miRNA-lncRNA sub-networks were established, among which all mRNAs and ncRNAs possessed significant predictive values for further studies for immune responses in the black rockfish.

Comprehensive CircRNA Profiling and Selection of Key CircRNAs Reveal the Potential Regulatory Roles of CircRNAs throughout Ovarian Development and Maturation in Cynoglossus semilaevis

Simple Summary

CircRNAs: as molecules involved in gene regulation, have become a new research hotspot in the non-coding RNA field. CircRNAs show tissue- or developmental stage-specific patterns of expression and can influence the expression levels of their parental genes. Recent studies have documented the potential biological roles of circRNAs in the growth, development, reproduction and health of humans and animals. Tongue sole (Cynoglossus semilaevis) is a marine flatfish that is an economically important farmed species in China. The commercial aquaculture of tongue sole has developed in the last few years because wild resources have decreased. Reproduction is regulated by brain-pituitary-gonad-liver axis which limits the development of artificial tongue sole culture. However, the roles of circRNAs in the ovarian development and maturation of tongue sole has never been reported. The identification of the potential functions of circRNAs provides a foundation for understanding the genetic mechanisms that regulate oocyte growth and maturation, which will allow the efficiency of tongue sole reproduction to be improved. Moreover, our findings extend the knowledge about a new type of endogenous RNA involved in regulating the ovarian development and maturation of tongue sole.

Abstract

CircRNAs are novel endogenous non-coding small RNAs involved in the regulation of multiple biological processes. However, little is known regarding circRNAs in ovarian development and maturation of fish. Our study, for the first time, provides the genome-wide overview of the types and relative abundances of circRNAs in tongue sole tissues during three ovarian developmental stages. We detected 6790 circRNAs in the brain, 5712 in the pituitary gland, 4937 in the ovary and 4160 in the liver. Some circRNAs exhibit tissue-specific expression, and qRT-PCR largely confirmed 6 differentially expressed (DE) circRNAs. Gene Ontology and KEGG pathway analyses of DE mRNAs were performed. Some DE circRNA parental genes were closely associated with biological processes in key signalling pathways and may play essential roles in ovarian development and maturation. We found that the selected circRNAs were involved in 10 pathways. RNase R digestion experiment and Sanger sequencing verified that the circRNA had a ring structure and was RNase R resistant. qRT-PCR results largely confirmed differential circRNA expression patterns from the RNA-seq data. These findings indicate that circRNAs are widespread in terms of present in production-related tissues of tongue sole with potentially important regulatory roles in ovarian development and maturation.

Revealing New Landscape of Turbot (Scophthalmus maximus) Spleen Infected with Aeromonas salmonicida through Immune Related circRNA-miRNA-mRNA Axis

Simple Summary

In this study, the expression of circRNAs, miRNAs, and mRNA in the immune organs spleen of turbot (Scophthalmus maximus) infected with Aeromonas salmonicida was analyzed by high-throughput sequencing, and circRNA-miRNA-mRNA network was constructed, so as to explore the function of non-coding RNA in the immune system of teleost. A total of 119, 140, and 510 differential expressed circRNAs, miRNAs, and mRNAs were identified in the infected groups compared with the uninfected group. The qRT-PCR verified the reliability and accuracy of the Illumina sequencing data. Fifteen triple networks of circRNA-miRNA-mRNA were presented in the form of “up (circRNA)-down (miRNA)-up (mRNA)” or “down-up-down”. Immune-related genes were also found in these networks. These results indicate that circRNAs and miRNAs may regulate the expression of immune-related genes through the circRNA-miRNA-mRNA regulatory network and thus participate in the immune response of turbot spleen after pathogen infection.

Abstract

Increasing evidence suggests that non-coding RNAs (ncRNA) play an important role in a variety of biological life processes by regulating gene expression at the transcriptional and post-transcriptional levels. Turbot (Scophthalmus maximus) has been threatened by various pathogens. In this study, the expression of circular RNAs (circRNAs), microRNAs (miRNAs), and mRNA in the immune organs spleen of turbot infected with Aeromonas salmonicida was analyzed by high-throughput sequencing, and a circRNA-miRNA-mRNA network was constructed, so as to explore the function of non-coding RNA in the immune system of teleost. Illumina sequencing was performed on the uninfected group and infected group. A total of 119 differential expressed circRNAs (DE-circRNAs), 140 DE-miRNAs, and 510 DE-mRNAs were identified in the four infected groups compared with the uninfected group. Most DE-mRNAs and the target genes of DE-ncRNAs were involved in immune-related pathways. The quantitative real-time PCR (qRT-PCR) results verified the reliability and accuracy of the high-throughput sequencing data. Ninety-six differentially expressed circRNA-miRNA-mRNA regulatory networks were finally constructed. Among them, 15 circRNA-miRNA-mRNA were presented in the form of “up (circRNA)-down (miRNA)-up (mRNA)” or “down-up-down”. Immune-related genes gap junction CX32.2, cell adhesion molecule 3, and CC chemokine were also found in these networks. These results indicate that ncRNA may regulate the expression of immune-related genes through the circRNA-miRNA-mRNA regulatory network and thus participate in the immune response of turbot spleen after pathogen infection.

Identification and Characterization of Long Non-coding RNAs in the Intestine of Olive Flounder (Paralichthys olivaceus) During Edwardsiella tarda Infection

Long non-coding RNAs (lncRNAs) play widespread roles in fundamental biological processes, including immune responses. The olive flounder (Paralichthys olivaceus), an important economical flatfish widely cultured in Japan, Korea, and China, is threatened by infectious pathogens, including bacteria, viruses, and parasites. However, the role of lncRNAs in the immune responses of this species against pathogen infections is not well-understood. Therefore, in this study, we aimed to identify lncRNAs in the intestine of olive flounder and evaluate their differential expression profiles during Edwardsiella tarda infection, which is an important zoonotic and intestinal pathogen. A total of 4,445 putative lncRNAs were identified, including 3,975 novel lncRNAs and 470 annotated lncRNAs. These lncRNAs had shorter lengths and fewer exons compared with mRNAs. In total, 115 differentially expressed lncRNAs (DE-lncRNAs) were identified during E. tarda infection. To validate the expression pattern of lncRNAs, six DE-lncRNAs were randomly selected for quantitative real-time PCR. The co-located and co-expressed mRNAs of DE-lncRNAs were predicted, which were used to conduct the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. The target genes of DE-lncRNAs enriched numerous immune-related processes and exhibited a strong correlation with immune-related signaling pathways. To better understand the extensive regulatory functions of lncRNAs, the lncRNA–miRNA–mRNA regulatory networks were constructed, and two potential competing endogenous RNA (ceRNA) networks, LNC_001979-novel_171-Potusc2 and LNC_001979-novel_171-Podad1, were preliminarily identified from the intestine of olive flounders for the first time. In conclusion, this study provides an invaluable annotation and expression profile of lncRNAs in the intestine of olive flounder infected with E. tarda; this forms a basis for further studies on the regulatory function of lncRNAs in the intestinal mucosal immune responses of olive flounder.

Knockdown of 91 H Suppresses the Tumorigenesis of Osteosarcoma via Inducing Methylation of CDK4 Promoter

Background:

Osteosarcoma is the most leading primary malignancy of the bone in adolescents all over the world. Long non-coding RNA (lncRNA) 91 H has been reported to participated in multiple cancers. Meanwhile, lncRNA 91 H has been proved to be upregulated in osteosarcoma. However, the function of 91 H in osteosarcoma remains unclear.

Methods:

Gene and protein expressions in osteosarcoma cells were detected by qRT-PCR and western blot, respectively. Cell viability was tested by CCK-8 assay. Ki67 staining was used to measure cell proliferation. Cell apoptosis and cycle were assessed by flow cytometry. In addition, transwell assay was used to detect cell migration and invasion. Furthermore, Methylation-specific PCR (MSP) was performed to test the methylation of CDK4 promoter. Finally, xenograft mice model was established to explore the role of 91 H in osteosarcoma in vivo.

Results:

Knockdown of 91 H significantly inhibited the growth of osteosarcoma cells via inducing the cell apoptosis. In addition, 91 H siRNA notably suppressed the migration and invasion of osteosarcoma cells. Meanwhile, knockdown of 91 H inhibited the progression of osteosarcoma via inducing methylation of CDK4 promoter. Furthermore, 91 H knockdown obviously induced G1 arrest in osteosarcoma cells via inhibition of PCNA and Cyclin D1. Finally, knockdown of 91 H notably inhibited the tumor growth of osteosarcoma in vivo.

Conclusion:

knockdown of 91 H suppressed the tumorigenesis of osteosarcoma via inducing methylation of CDK4 promoter in vitro and in vivo. Thus, 91 H may serve as a new target for the treatment of osteosarcoma.