Comparative analysis of the miRNA-mRNA regulation networks in turbot (Scophthalmus maximus L.) following Vibrio anguillarum infection

Identification and Characterization of Immune-Associated MicroRNAs in Silver Carp (Hypophthalmichthys molitrix) Responding to Aeromonas veronii and LPS Stimulation

The ubiquitous Gram-negative bacterial pathogen Aeromonas veronii (A. veronii) can easily cause inflammatory reactions in aquatic organisms, resulting in high mortality and huge economic losses. MicroRNAs (miRNAs) participate in immune regulation and have certain conserved properties. MiRNAs are involved in the immune responses of a variety of teleost fish infected with bacteria, whereas there is no related report in silver carp (Hypophthalmichthys molitrix). Therefore, we identified the expression profiles of miRNA in silver carp stimulated by A. veronii and LPS. Among them, the quantity of differentially expressed miRNAs (DEmiRNAs) obtained in the silver carp challenge group was 73 (A. veronii) and 90 (LPS). The GO enrichment and analysis of KEGG pathways have shown that the predicted target genes are mainly associated with lipid metabolism and the immune response in silver carp. This indicates the possibility that miRNAs play a role in regulating immune-related pathways. In addition, a total of eight DEmiRNAs validated the accuracy of the sequencing result via quantitative real-time PCR (qRT-PCR). Finally, we selected the silver carp head kidney macrophage cells (HKCs) as model cells and proved that miR-30b-5p can regulate the inflammatory response in silver carp HKCs. This study lays the foundation for exploring miRNA regulation in silver carp during pathogenic bacterial infection. In addition, it provides a reference for the future development of non-coding RNA antibacterial drugs.

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.

Genome-Wide Identification, Molecular Characterization, and Involvement in Response to Abiotic and Biotic Stresses of the HSP70 Gene Family in Turbot (Scophthalmus maximus)

Heat shock proteins 70 (HSP70s) are known to play essential roles in organisms' response mechanisms to various environmental stresses. However, no systematic identification and functional analysis has been conducted for HSP70s in the turbot (Scophthalmus maximus), a commercially important worldwide flatfish. Herein, 16 HSP70 genes unevenly distributed on nine chromosomes were identified in the turbot at the genome-wide level. Analyses of gene structure, motif composition, and phylogenetic relationships provided valuable data on the HSP70s regarding their evolution, classification, and functional diversity. Expression profiles of the HSP70 genes under five different stresses were investigated by examining multiple RNA-seq datasets. Results showed that 10, 6, 8, 10, and 9 HSP70 genes showed significantly up- or downregulated expression after heat-induced, salinity-induced, and Enteromyxum scophthalmi, Vibrio anguillarum, and Megalocytivirus infection-induced stress, respectively. Among them, hsp70 (hspa1a), hspa1b, and hspa5 showed significant responses to each kind of induced stress, and qPCR analyses further validated their involvement in comprehensive anti-stress, indicating their involvement in organisms' anti-stress mechanisms. These findings not only provide new insights into the biological function of HSP70s in turbot adapting to various environmental stresses, but also contribute to the development of molecular-based selective breeding programs for the production of stress-resistant turbot strains in the aquaculture industry.

Derivation and characterization of new cell line from intestine of turbot (Scophthalmus maximus)

A continuous intestine cell line from turbot (Scophthalmus maximus) designated as SMI was established utilizing the tissue explant technique. Primary SMI cell was cultured at 24 °C in a medium with 20% fetal bovine serum (FBS), then subcultured in 10% FBS after 10 passages. Impacts of medium or temperature on the growth of SMI were examined and the results indicated it grew well in DMEM supplemented with 10% FBS at 24 °C. The SMI cell line was subcultured more than 60 times. Karyotyping, chromosome number, and ribosomal RNA genotyping analysis revealed that SMI had a modal diploid chromosome number of 44 and originated from turbot. After being transfected with pEGFP-N1 and FAM-siRNA, a large number of green fluorescence signals were observed in SMI, indicating that SMI could be used as an ideal platform to explore gene function in vitro. In addition, the expression of epithelium-associated genes such as itga6, itgb4, gja1, claudin1, zo-1, and E-cadherin in SMI suggested the SMI had some characteristics of epidermal cells. The upregulation of immune-associated genes such as TNF-β, NF-κB, and IL-1β in SMI after stimulation with pathogen-associated molecular patterns suggested the SMI might exhibit immune functions similar to the intestinal epithelium in vivo.

Identification of stress-related genes by co-expression network analysis based on the improved turbot genome

Turbot (Scophthalmus maximus), commercially important flatfish species, is widely cultivated in Europe and China. With the continuous expansion of the intensive breeding scale, turbot is exposed to various stresses, which greatly impedes the healthy development of turbot industry. Here, we present an improved high-quality chromosome-scale genome assembly of turbot using a combination of PacBio long-read and Illumina short-read sequencing technologies. The genome assembly spans 538.22 Mb comprising 27 contigs with a contig N50 size of 25.76 Mb. Annotation of the genome assembly identified 104.45 Mb repetitive sequences, 22,442 protein-coding genes and 3,345 ncRNAs. Moreover, a total of 345 stress responsive candidate genes were identified by gene co-expression network analysis based on 14 published stress-related RNA-seq datasets consisting of 165 samples. Significantly improved genome assembly and stress-related candidate gene pool will provide valuable resources for further research on turbot functional genome and stress response mechanism, as well as theoretical support for the development of molecular breeding technology for resistant turbot varieties.

Measurement(s)	whole genome sequencing
Technology Type(s)	PacBio long-read and Illumina short-read sequencing technologies