RNA-seq and microRNA-seq analysis of Japanese flounder (Paralichthys olivaceus) larvae treated by thyroid hormones

Transcriptome analysis of the Japanese eel (Anguilla japonica) during larval metamorphosis

BACKGROUND

Anguillid eels spend their larval period as leptocephalus larvae that have a unique and specialized body form with leaf-like and transparent features, and they undergo drastic metamorphosis to juvenile glass eels. Less is known about the transition of leptocephali to the glass eel stage, because it is difficult to catch the metamorphosing larvae in the open ocean. However, recent advances in rearing techniques for the Japanese eel have made it possible to study the larval metamorphosis of anguillid eels. In the present study, we investigated the dynamics of gene expression during the metamorphosis of Japanese eel leptocephali using RNA sequencing.

RESULTS

During metamorphosis, Japanese eels were classified into 7 developmental stages according to their morphological characteristics, and RNA sequencing was used to collect gene expression data from each stage. A total of 354.8 million clean reads were generated from the body and 365.5 million from the head, after the processing of raw reads. For filtering of genes that characterize developmental stages, a classification model created by a Random Forest algorithm was built. Using the importance of explanatory variables feature obtained from the created model, we identified 46 genes selected in the body and 169 genes selected in the head that were defined as the "most characteristic genes" during eel metamorphosis. Next, network analysis and subsequently gene clustering were conducted using the most characteristic genes and their correlated genes, and then 6 clusters in the body and 5 clusters in the head were constructed. Then, the characteristics of the clusters were revealed by Gene Ontology (GO) enrichment analysis. The expression patterns and GO terms of each stage were consistent with previous observations and experiments during the larval metamorphosis of the Japanese eel.

CONCLUSION

Genome and transcriptome resources have been generated for metamorphosing Japanese eels. Genes that characterized metamorphosis of the Japanese eel were identified through statistical modeling by a Random Forest algorithm. The functions of these genes were consistent with previous observations and experiments during the metamorphosis of anguillid eels.

mRNA-miRNA sequencing reveals mechanisms of 2,2'-dipyridyl disulfide-induced thyroid disruption in Japanese flounder (Paralichthys olivaceus)

This study was conducted to evaluate the thyroid-disrupting effects of 2,2'-dipyridyl disulfide using Japanese flounder (Paralichthys olivaceus) as an animal model and to reveal the underlying mechanisms from the perspective of miRNA-mRNA interactions. The results indicated that 2,2'-dipyridyl disulfide exposure decelerated the metamorphic progress of P. olivaceus, suggesting its thyroid-disrupting property as an antagonist. Furthermore, radioimmunoassays, thyroid histological observation, real-time polymerase chain reaction, and mRNA sequencing showed that 2,2'-dipyridyl disulfide exposure exerted its thyroid-disrupting effects on larval and juvenile P. olivaceus by targeting multiple processes and pathways involved in the thyroid system, including peripheral metabolism of thyroid hormones, the thyroid hormone synthesis pathway, and the thyroid hormone/thyroid hormone receptor signaling pathway. In particular, global upregulation of the gene expression of three deiodinases caused decreases in thyroid hormone levels after 2,2'-dipyridyl disulfide exposure that are believed to be responsible for the inhibition of metamorphosis in P. olivaceus. Finally, miRNA sequencing suggested that several evolutionarily conserved miRNAs play important roles in the mechanism of 2,2'-dipyridyl disulfide-induced thyroid disruption. Specifically, overexpression of pny-miR-723a and pny-miR-216a resulted in upregulation of deiodinase 1 mRNA levels in the 2,2'-dipyridyl disulfide exposure group. This study provides the first evidence that 2,2'-dipyridyl disulfide has thyroid-disrupting properties and is also the first study remarking on the roles of miRNA-mRNA interactions in the action mechanisms of thyroid disruptors.

Identification of Potential Blind-Side Hypermelanosis-Related lncRNA–miRNA–mRNA Regulatory Network in a Flatfish Species, Chinese Tongue Sole (Cynoglossus semilaevis)

Blind-side hypermelanosis has emerged as a major concern in commercial rearing environments of the flatfish aquaculture industry. To date, the underlying molecular mechanisms are not well understood. To fill this gap, in this study, whole transcriptomic sequencing and analyses were performed using normal skins and hypermelanic skins of the blind side of Chinese tongue sole (Cynoglossus semilaevis). Differentially expressed long non-coding RNAs (DElncRNAs), miRNAs (DEmiRNAs), and differentially expressed genes as well as their competing endogenous RNA (ceRNA) networks were identified. A total of 34 DElncRNAs, 226 DEmiRNAs, and 610 DEGs were identified. Finally, lncRNA–miRNA–mRNA regulatory networks (involving 29 DElncRNAs, 106 DEmiRNAs, and 162 DEGs) associated with blind-side hypermelanosis were constructed. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of 162 DEGs in ceRNA networks identified DEGs (e.g., oca2, mc1r, and ihhb) in pigmentation-related biological processes and DEGs (e.g., ca4, glul, and fut9) in nitrogen metabolism, glycosphingolipid biosynthesis, and folate biosynthesis pathways, as well as their corresponding DElncRNAs and DEmiRNAs to potentially play key regulatory roles in blind-side hypermelanosis. In conclusion, this is the first study on the ceRNA regulatory network associated with blind-side hypermelanosis in flatfish. These new findings expand the spectrum of non-coding regulatory mechanisms underpinning blind-side hypermelanosis, which facilitates the further exploration of molecular regulatory mechanisms of malpigmentation in flatfish.

Integrated microRNA-mRNA analysis reveals the roles of microRNAs in the muscle fat metabolism of Yanbian cattle

Fat deposition is an important economic trait in farm animals. However, it is difficult to genetically improve intramuscular fat deposition via trait-based cattle breeding. The main objectives of this study were to analyze the factors about beef flavor, and to detect functional microRNA (miRNA, miR) associated with intramuscular fat deposition in Yanbian cattle. Longissimus dorsi samples from six steers were separated into high- and low-fat groups (n = 3 each) based on the marbling score, and transcriptomic analysis was performed using miRNA sequencing. A total of 33 miRNAs and 38 genes were found to be differentially expressed in the high- and low-fat groups. Quantitative real-time polymerase chain reaction was performed to validate the sequencing results. Integrated miRNA-mRNA analysis revealed that miRNA-associated target genes were primarily associated with skeletal muscle development. However, some of the miRNAs (miR-424 etc.) and genes (ATF3 etc.) were also associated with fat metabolism. A targeted relationship between miR-22-3p and the WFIKKN2 gene and its involvement in adipocyte differentiation were confirmed experimentally. The study findings may provide potential candidate molecular targets for the selection of cattle with improved meat quality.