Genome-Wide Integrated Analysis Revealed Functions of lncRNA-miRNA-mRNA Interaction in Growth of Intermuscular Bones in Megalobrama amblycephala

Effect of runx2b deficiency in intermuscular bones on the regulatory network of lncRNA-miRNA-mRNA

Intermuscular bones (IBs) are mineralized spicules that negatively impact the quality and value of fish products. Runx2b is a crucial modulator in promoting bone formation through regulating osteoblast differentiation. Previous studies suggested that loss of runx2b gene completely inhibited IBs formation in zebrafish. However, how the whole transcriptome, including mRNA and non-coding RNA (ncRNA), affects the IBs development in runx2b-/- zebrafish are not known. The aim of this study was to identify the regulatory networks of differentially expressed (DE) lncRNAs, miRNAs, and mRNAs in zebrafish with and without IBs (runx2b+/+ fish and runx2b-/- fish) utilizing high-throughput sequencing techniques. All together there are 1051 mRNAs, 456 lncRNAs, and 18 miRNAs differentially expressed were found between these two strains. The analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) has highlighted significant pathways linked to the development of IBs, specifically the TGF-beta and Wnt signaling pathways, and a number of genes concentrated on these two signaling pathways related to the formation of IBs. Further, 1989 competing endogenous RNA (ceRNA) networks were created according to the correlation among mRNAs, miRNAs and lncRNAs. The ceRNA networks results revealed 52 ceRNA pairs related to the IBs formation, consisting of 52 mRNAs, 37 lncRNAs, and 6 miRNAs. Of these, we found that dre-miR-2189 was the key element of ceRNA pairs, interacting with 19 mRNAs and 11 lncRNAs, and MSTRG.13175.1 could regulate sp7 expression by interacting with dre-miR-2189 to function in osteogenic differentiation. Subsequent experiments at the cellular level also revealed the interaction mechanism. The outcomes indicated a crucial role of miRNAs and lncRNAs in the development of fish IBs, which offer new views into the functions of ncRNAs involved in IBs formation.

MicroRNA Regulation in Infectious Diseases and Its Potential as a Biosensor in Future Aquaculture Industry: A Review

An infectious disease is the most apprehensive problem in aquaculture as it can lead to high mortality in aquatic organisms and massive economic loss. Even though significant progress has been accomplished in therapeutic, prevention, and diagnostic using several potential technologies, more robust inventions and breakthroughs should be achieved to control the spread of infectious diseases. MicroRNA (miRNA) is an endogenous small non-coding RNA that post-transcriptionally regulates the protein-coding genes. It involves various biological regulatory mechanisms in organisms such as cell differentiation, proliferation, immune responses, development, apoptosis, and others. Furthermore, an miRNA also acts as a mediator to either regulate host responses or enhance the replication of diseases during infection. Therefore, the emergence of miRNAs could be potential candidates for the establishment of diagnostic tools for numerous infectious diseases. Interestingly, studies have revealed that miRNAs can be used as biomarkers and biosensors to detect diseases, and can also be used to design vaccines to attenuate pathogens. This review provides an overview of miRNA biogenesis and specifically focuses on its regulation during infection in aquatic organisms, especially on the host immune responses and how miRNAs enhance the replication of pathogens in the organism. In addition to that, we explored the potential applications, including diagnostic methods and treatments, that can be employed in the aquaculture industry.

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.

Profiling the Spatial Expression Pattern and ceRNA Network of lncRNA, miRNA, and mRNA Associated with the Development of Intermuscular Bones in Zebrafish

Intermuscular bones (IBs) are small spicule-like bones in the muscular septum of fish, which affect their edible and economic value. The molecular mechanism of IB development is still uncertain. Numerous studies have shown that the ceRNA network, which is composed of mRNA, lncRNA, and miRNA, plays an important regulatory role in bone development. In this study, we compared the mRNA, lncRNA, and miRNA expression profiles in different IB development segments of zebrafish. The development of IBs includes two main processes, which are formation and growth. A series of genes implicated in the formation and growth of IBs were identified through gene differential expression analysis and expression pattern analysis. Functional enrichment analysis showed that the functions of genes implicated in the regulation of the formation and growth of IBs were quite different. Ribosome and oxidative phosphorylation signaling pathways were significantly enriched during the formation of IBs, suggesting that many proteins are required to form IBs. Several pathways known to be associated with bone development have been shown to play an important role in the growth of IBs, including calcium, ECM-receptor interaction, Wnt, TGF-β, and hedgehog signaling pathways. According to the targeting relationship and expression correlation of mRNA, lncRNA, and miRNA, the ceRNA networks associated with the growth of IBs were constructed, which comprised 33 mRNAs, 9 lncRNAs, and 7 miRNAs. This study provides new insight into the molecular mechanism of the development of IBs.

Genome assembly and isoform analysis of a highly heterozygous New Zealand fisheries species, the tarakihi (Nemadactylus macropterus)

Although being some of the most valuable and heavily exploited wild organisms, few fisheries species have been studied at the whole-genome level. This is especially the case in New Zealand, where genomics resources are urgently needed to assist fisheries management. Here, we generated 55 Gb of short Illumina reads (92× coverage) and 73 Gb of long Nanopore reads (122×) to produce the first genome assembly of the marine teleost tarakihi [Nemadactylus macropterus (Forster, 1801)], a highly valuable fisheries species in New Zealand. An additional 300 Mb of Iso-Seq reads were obtained to assist in gene annotation. The final genome assembly was 568 Mb long with an N50 of 3.37 Mb. The genome completeness was high, with 97.8% of complete Actinopterygii Benchmarking Universal Single-Copy Orthologs. Heterozygosity values estimated through k-mer counting (1.00%) and bi-allelic SNPs (0.64%) were high compared with the same values reported for other fishes. Iso-Seq analysis recovered 91,313 unique transcripts from 15,515 genes (mean ratio of 5.89 transcripts per gene), and the most common alternative splicing event was intron retention. This highly contiguous genome assembly and the isoform-resolved transcriptome will provide a useful resource to assist the study of population genomics and comparative eco-evolutionary studies in teleosts and related organisms.

Screening and identification of miRNAs regulating Tbx4/5 genes of Pampus argenteus

Background

Silver pomfret (Pampus argenteus) is one of the most widely distributed and economically important pelagic fish species. However, an unique morphological feature of P. argenteus is the loss of pelvic fins, which can increase the energy requirement during food capture to some extent and is therefore not conducive to artificial culture. Tbx4/5 genes are highly conserved regulatory factors that regulate limb development in vertebrates and are in turn regulated by microRNAs (miRNAs). However, the miRNAs that directly regulate the Tbx4/5 genes in P. argenteus remain to be elucidated.

Methods

The Tbx4/5 genes of P. argenteus were first cloned, and the small RNA transcriptomes were sequenced by high-throughput sequencing during the critical period of the fin development at days 1, 7, and 13 of hatching. The miRNAs regulating the Tbx4/5 genes of P. argenteus were subsequently predicted by bioinformatics analysis, and the related miRNAs were verified in vitro using a dual fluorescence reporter system.

Results

A total of 662 miRNAs were identified, of which 257 were known miRNAs and 405 were novel miRNAs were identified. Compared to day 1, 182 miRNAs were differentially expressed (DE) on day 7, of which 77 and 105 miRNAs were downregulated and upregulated, respectively, while 278 miRNAs were DE on day 13, of which 136 and 142 miRNAs were downregulated and upregulated, respectively. Compared to day 13, four miRNAs were DE on day 7, of which three miRNAs were downregulated and one miRNA was upregulated. The results of hierarchical clustering of the miRNAs revealed that the DE genes were inversely expressed between days 1 and 7, and between days 1 and 13 of larval development, indicating that the larvae were in the peak stage of differentiation. However, the number of DE genes between days 7 and 13 of larval development was relatively small, suggesting the initiation of development. The potential target genes of the DE miRNAs were subsequently predicted, and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of target genes were performed. The results suggested that the DE miRNAs were involved in growth, development, and signal transduction pathways, of which the Wnt and Fgfs signaling pathways are known to play important roles in the growth and development of fins. The results of dual fluorescence reporter assays demonstrated that miR-102, miR-301c, and miR-589 had a significant negative regulatory effect on the 3'-UTR of the Tbx4 gene, while miR-187, miR-201, miR-219, and miR-460 had a significant negative regulatory effect on the 3'-UTR of the Tbx5 gene. Altogether, the findings indicated that miRNAs play an important role in regulating the growth and development of pelvic fins in P. argenteus. This study provides a reference for elucidating the interactions between the miRNAs and target genes of P. argenteus in future studies.

Transcriptome-metabolome analysis reveals how sires affect meat quality in hybrid sheep populations

Crossbreeding improves and enhances meat quality and is widely used in sheep production; however, the molecular mechanisms underlying the meat quality of various crossbred sheep remain unknown. In this study, male Southdown, Suffolk and Hu sheep were crossbred with female Hu sheep, and the transcriptomes and metabolomes of the longissimus dorsi muscle of the F1 generation were sequenced to explore how different sire breeds affect meat quality. The results showed that 631 differentially expressed genes and 119 significantly altered metabolites contributed to muscle development characteristics and meat quality-related diversity (P < 0.05). These genes and metabolites were significantly enriched in lipid metabolism pathways, including arachidonic acid metabolism and PPAR signaling. Several candidate genes were associated with muscle growth, such as MYLK3, MYL10, FIGN, MYH8, MYOM3, LMCD1, and FLRT1. Among these, MYH8 and MYL10 participated in regulating muscle growth and development and were correlated with meat quality-related fatty acid levels (|r| > 0.5 and p < 0.05). We selected mRNA from four of these genes to verify the accuracy of the sequencing data via qRT-PCR. Our findings provide further insight into the key genes and metabolites involved in muscle growth and meat quality in hybrid sheep populations.

The Emerging Role of Long Non-Coding RNAs in Development and Function of Gilthead Sea Bream (Sparus aurata) Fast Skeletal Muscle

Long non-coding RNAs (lncRNAs) are an emerging group of ncRNAs that can modulate gene expression at the transcriptional or translational levels. In the present work, previously published transcriptomic data were used to identify lncRNAs expressed in gilthead sea bream skeletal muscle, and their transcription levels were studied under different physiological conditions. Two hundred and ninety lncRNAs were identified and, based on transcriptomic differences between juveniles and adults, a total of seven lncRNAs showed potential to be important for muscle development. Our data suggest that the downregulation of most of the studied lncRNAs might be linked to increased myoblast proliferation, while their upregulation might be necessary for differentiation. However, with these data, as it is not possible to propose a formal mechanism to explain their effect, bioinformatic analysis suggests two possible mechanisms. First, the lncRNAs may act as sponges of myoblast proliferation inducers microRNAs (miRNAs) such as miR-206, miR-208, and miR-133 (binding energy MEF < -25.0 kcal). Secondly, lncRNA20194 had a strong predicted interaction towards the myod1 mRNA (ndG = -0.17) that, based on the positive correlation between the two genes, might promote its function. Our study represents the first characterization of lncRNAs in gilthead sea bream fast skeletal muscle and provides evidence regarding their involvement in muscle development.

Prognostic significance of ferroptosis-related genes and their methylation in AML

BACKGROUND

Ferroptosis involves in the development and therapeutic response of various types of tumors. This study aims to explore ferroptosis-related prognostic genes that could further accurately stratify AML patients.

METHODS

We investigated the prognosis significance of ferroptosis-related genes in AML by Univariate and multivariate Cox proportional hazards regression analyses. With the methylation data of TCGA samples, we looked for methylation sites associated with prognostic genes and compared the correlation between methylation and mRNA expression. R software and 'edgeR' packages were used to identify the DEGs between the high-and-low-risk groups divided by the FRPGs prognosis model and then run GO enrichment, KEGG pathway, and PPI network.

RESULTS

We found a prognostic risk model that included AKR1C2 and SOCS1 predicted outcomes in AML patients. Methylation analysis showed that AKR1C2 and SOCS1 are negatively regulated by their methylation, leading to their low expression in AML patients. Besides, both decreased SOCS1 expression and hypermethylation predicted favorable OS and PFS in AML patients. Finally, this prognostic risk model exhibited a close correlation with several clinical features, especially with age (P=0.005), cytogenetic type (P=0.031), risk_cytogenetic (P=0.001), and risk_molecular (P<0.001). Functional enrichment analysis showed that DEGs are most enriched in the regulation of cell death and the PI3K-Akt signaling pathway.

CONCLUSION

AKR1C2 and SOCS1 are promising biomarkers for predicting prognosis in patients with AML.

Genome-Wide Association Study on Reproduction-Related Body-Shape Traits of Chinese Holstein Cows

Reproduction is an important production activity for dairy cows, and their reproductive performance can directly affect the level of farmers' income. To better understand the genomic regions and biological pathways of reproduction-related traits of dairy cows, in the present study, three body shape traits-Loin Strength (LS), Rump Angle (RA), and Pin Width (PW)-were selected as indicators of the reproductive ability of cows, and we conducted genome-wide association analyses on them. The heritability of these three traits was medium, ranging from 0.20 to 0.38. A total of 11 significant single-nucleotide polymorphisms (SNPs) were detected associated with these three traits. Bioinformatics analysis was performed on genes close to the significant SNPs (within 200 Kb) of LS, RA, and PW, and we found that these genes were totally enriched in 20 gene ontology terms and six KEGG signaling pathways. Finally, the five genes CDH12, TARP, PCDH9, DTHD1, and ARAP2 were selected as candidate genes that might affect LS. The six genes LOC781835, FSTL4, ATG4C, SH3BP4, DMP1, and DSPP were selected as candidate genes that might affect RA. The five genes USP6NL, CNTN3, LOC101907665, UPF2, and ECHDC3 were selected as candidate genes that might affect the PW of Chinese Holstein cows. Our results could provide useful biological information for the improvement of body shape traits and contribute to the genomic selection of Chinese Holstein cows.