Characterization and expression analysis of microRNAs in the tube foot of sea cucumber Apostichopus japonicus

Identification and functional analysis of Dmrt1 gene and the SoxE gene in the sexual development of sea cucumber, Apostichopus japonicus

Members of the Doublesex and Mab-3-related transcription factor (Dmrt) gene family handle various vital functions in several biological processes, including sex determination/differentiation and gonad development. Dmrt1 and Sox9 (SoxE in invertebrates) exhibit a very conserved interaction function during testis formation in vertebrates. However, the dynamic expression pattern and functional roles of the Dmrt gene family and SoxE have not yet been identified in any echinoderm species. Herein, five members of the Dmrt gene family (Dmrt1, 2, 3a, 3b and 5) and the ancestor SoxE gene were identified from the genome of Apostichopus japonicus. Expression studies of Dmrt family genes and SoxE in different tissues of adult males and females revealed different expression patterns of each gene. Transcription of Dmrt2, Dmrt3a and Dmrt3b was higher expressed in the tube feet and coelomocytes instead of in gonadal tissues. The expression of Dmrt1 was found to be sustained throughout spermatogenesis. Knocking-down of Dmrt1 by means of RNA interference (RNAi) led to the downregulation of SoxE and upregulation of the ovarian regulator foxl2 in the testes. This indicates that Dmrt1 may be a positive regulator of SoxE and may play a role in the development of the testes in the sea cucumber. The expression level of SoxE was higher in the ovaries than in the testes, and knocking down of SoxE by RNAi reduced SoxE and Dmrt1 expression but conversely increased the expression of foxl2 in the testes. In summary, this study indicates that Dmrt1 and SoxE are indispensable for testicular differentiation, and SoxE might play a functional role during ovary differentiation in the sea cucumber.

Identification of molecular markers for superior quantitative traits in a novel sea cucumber strain by comparative microRNA-mRNA expression profiling

To investigate the adaptability of Apostichopus japonicus (A. japonicus) strain "Anyuan No. 1" in the South China Sea, field monitoring and microRNA-mRNA integrated analyses were conducted between "Anyuan No. 1" and a regular A. japonicus population from Wendeng (Shandong Province, as a control) in the Xiapu farming area in Fujian Province, China. The results showed that "Anyuan No. 1" exhibited greater body weight increase and a higher number of papillae compared to the control during two and a half months of field monitoring. Comparative microRNA (miRNA) and mRNA transcriptome analyses identified 12 differentially expressed miRNAs (DEMs) and 165 differentially expressed genes (DEGs) in "Anyuan No. 1" compared to the control. Long-chain specific acyl-CoA dehydrogenase (ACADL), transmembrane protein 251 (TMEM251), dehydrogenase/reductase SDR family protein 7-like (Dhrs7), insulin-like growth factor-binding protein 7 (IGFBP-7), CDK5 regulatory subunit-associated protein 1 (CDK5RAP1), visual pigment-like receptor peropsin, 39S ribosomal protein, miR-10, miR-153, miR-7, and miR-3529 were identified as gene and miRNA candidates correlated with superior economic traits in "Anyuan No. 1". Collectively, "Anyuan No. 1" is suitable for large-scale cultivation extension due to its better adaptability to the South China Sea area. Furthermore, we identified "miR10-ACADL" as a potential module for further molecular marker-assisted selective breeding of A. japonicus.

Identification and characterization of microRNAs from the tube foot in the sea urchin Strongylocentrotus intermedius

MicroRNAs (miRNAs) play critical roles in regulating many bio-processes of eukaryotes. The sea urchin Strongylocentrotus intermedius (an important fishery resource) is of great economic importance in Japan, North Korea, Russia, and China. In the current study, miRNAs of tube foot in S. intermedius were firstly identified and characterized. Data in this study can provide more genomic information for the further understanding of the complex regulation network in sea urchins and present a new way for monitoring the health status of cultured sea urchins.

Long Non-Coding RNAs (lncRNAs) of Sea Cucumber: Large-Scale Prediction, Expression Profiling, Non-Coding Network Construction, and lncRNA-microRNA-Gene Interaction Analysis of lncRNAs in Apostichopus japonicus and Holothuria glaberrima During LPS Challenge and Radial Organ Complex Regeneration

Long non-coding RNA (lncRNA) structurally resembles mRNA but cannot be translated into protein. Although the systematic identification and characterization of lncRNAs have been increasingly reported in model species, information concerning non-model species is still lacking. Here, we report the first systematic identification and characterization of lncRNAs in two sea cucumber species: (1) Apostichopus japonicus during lipopolysaccharide (LPS) challenge and in heathy tissues and (2) Holothuria glaberrima during radial organ complex regeneration, using RNA-seq datasets and bioinformatics analysis. We identified A. japonicus and H. glaberrima lncRNAs that were differentially expressed during LPS challenge and radial organ complex regeneration, respectively. Notably, the predicted lncRNA-microRNA-gene trinities revealed that, in addition to targeting protein-coding transcripts, miRNAs might also target lncRNAs, thereby participating in a potential novel layer of regulatory interactions among non-coding RNA classes in echinoderms. Furthermore, the constructed coding-non-coding network implied the potential involvement of lncRNA-gene interactions during the regulation of several important genes (e.g., Toll-like receptor 1 [TLR1] and transglutaminase-1 [TGM1]) in response to LPS challenge and radial organ complex regeneration in sea cucumbers. Overall, this pioneer systematic identification, annotation, and characterization of lncRNAs in echinoderm pave the way for similar studies and future genetic, genomic, and evolutionary research in non-model species.

High throughput sequencing of small RNAs transcriptomes in two Crassostrea oysters identifies microRNAs involved in osmotic stress response

Increasing evidence suggests that microRNAs post-transcriptionally regulate gene expression and are involved in responses to biotic and abiotic stress. However, the role of miRNAs involved in osmotic plasticity remains largely unknown in marine bivalves. In the present study, we performed low salinity challenge with two Crassostrea species (C. gigas and C. hongkongensis), and conducted high-throughput sequencing of four small RNA libraries constructed from the gill tissues. A total of 202 and 87 miRNAs were identified from C. gigas and C. hongkongensis, respectively. Six miRNAs in C. gigas and two in C. hongkongensis were differentially expressed in response to osmotic stress. The expression profiles of these eight miRNAs were validated by qRT-PCR. Based on GO enrichment and KEGG pathway analysis, genes associated with microtubule-based process and cellular component movement were enriched in both species. In addition, five miRNA-mRNA interaction pairs that showed opposite expression patterns were identified in the C. hongkongensis, Differential expression analysis identified the miRNAs that play important regulatory roles in response to low salinity stress, providing insights into molecular mechanisms that are essential for salinity tolerance in marine bivalves.

Expression analysis of microRNAs related to the skin ulceration syndrome of sea cucumber Apostichopus japonicus

MicroRNAs (miRNAs) are small non-coding RNAs that are involved in many biological processes. To investigate the miRNAs related to skin ulceration syndrome (SUS) of Apostichopus japonicus, small RNA libraries of body wall, intestine, respiratory tree and coelomocytes from healthy and diseased A. japonicus were sequenced on Illumina Hiseq 2000 platform. A total of 247 conserved and 10 novel miRNAs were identified across all libraries. After pair-wise comparisons, 215 miRNAs in body wall, 36 in intestine, 2 in respiratory tree and 38 in coelomocytes showed significant expression differences. Further analyses were conducted on some tissue-specific differentially expressed miRNAs: miR-8 and miR-486-5p in body wall, miR-200-3p, let-7-5p and miR-125 in intestine, miR-278a-3p and bantam in respiratory, miR-10a and miR-184 in coelomocytes. Notably, these miRNAs in some species were reported to function in various physiological or pathological processes associated with immune regulations. Using stem-loop quantitative real time PCR, six representative miRNAs in four tissues were selected to validate the sequencing results. The Pearson's correlation coefficient (R) of the six miRNAs ranged from 0.777 to 0.948, which confirmed the consistency and accuracy between these two approaches. This study provides comprehensive expression and regulation patterns of functional miRNAs in different tissues and gives insights into the tissue-specific immune response mechanisms in SUS-infected A. japonicus.

Identification and Characterization of MicroRNAs from Longitudinal Muscle and Respiratory Tree in Sea Cucumber (Apostichopus japonicus) Using High-Throughput Sequencing

MicroRNAs (miRNAs), as a family of non-coding small RNAs, play important roles in the post-transcriptional regulation of gene expression. Sea cucumber (Apostichopus japonicus) is an important economic species which is widely cultured in East Asia. The longitudinal muscle (LTM) and respiratory tree (RPT) are two important tissues in sea cucumber, playing important roles such as respiration and movement. In this study, we identified and characterized miRNAs in the LTM and RPT of sea cucumber (Apostichopus japonicus) using Illumina HiSeq 2000 platform. A total of 314 and 221 conserved miRNAs were identified in LTM and RPT, respectively. In addition, 27 and 34 novel miRNAs were identified in the LTM and RPT, respectively. A set of 58 miRNAs were identified to be differentially expressed between LTM and RPT. Among them, 9 miRNAs (miR-31a-3p, miR-738, miR-1692, let-7a, miR-72a, miR-100b-5p, miR-31b-5p, miR-429-3p, and miR-2008) in RPT and 7 miRNAs (miR-127, miR-340, miR-381, miR-3543, miR-434-5p, miR-136-3p, and miR-300-3p) in LTM were differentially expressed with foldchange value being greater than 10. A total of 14,207 and 12,174 target genes of these miRNAs were predicted, respectively. Functional analysis of these target genes of miRNAs were performed by GO analysis and pathway analysis. This result provided in this work will be useful for understanding biological characteristics of the LTM and RPT of sea cucumber and assisting molecular breeding of sea cucumber for aquaculture.

Expression profiling analysis of the microRNA response of Cynoglossus semilaevis to Vibrio anguillarum and other stimuli

To investigate the roles of microRNAs (miRNA) of Cynoglossus semilaevis in response to Vibrio anguillarum that were previously identified using high-throughput sequencing, microarray analyses was performed on three small RNA libraries (CG, NOSG, and HOSG) prepared from C. semilaevis immune tissues. In total, of 1279 designed probes, 739 (57.78 %) were detectable. The expression levels of these miRNAs were analyzed using pairwise comparisons among the three libraries, and a total of 99 miRNAs were observed to be significantly differentially expressed. The expression patterns of 10 differentially expressed miRNAs were validated by real-time quantitative PCR (RT-qPCR). In addition, expression of miR-142-5p, miR-223, and miR-181a in response to V. anguillarum at numerous time-points in four tissues, as well as the responses to lipopolysaccharide (LPS), polyinosinic:polycytidylic acid (poly I:C), peptidoglycan (PGN), and red-spotted grouper nervous necrosis virus (RGNNV) in head kidney cells, were studied by qRT-PCR. Taken together, all of the expression profiles showed significant differences compared to the control group; both similarities and differences in the expression responses to the same pathogen were observed. Collectively, these findings highlighted the putative roles for miRNAs in the context of the innate immune response of C. semilaevis exposing to pathogens and that further studies are needed to understand the molecular mechanisms of miRNA regulation in C. semilaevis host-pathogen interactions.