Embryonic temperature has long-term effects on muscle circRNA expression and somatic growth in Nile tilapia

Embryonic temperature has a lasting impact on muscle phenotype in vertebrates, involving complex molecular mechanisms that encompass both protein-coding and non-coding genes. Circular RNAs (circRNAs) are a class of regulatory RNAs that play important roles in various biological processes, but the effect of variable thermal conditions on the circRNA transcriptome and its long-term impact on muscle growth plasticity remains largely unexplored. To fill this knowledge gap, we performed a transcriptomic analysis of circRNAs in fast muscle of Nile tilapia (Oreochromis niloticus) subjected to different embryonic temperatures (24°C, 28°C and 32°C) and then reared at a common temperature (28°C) for 4 months. Nile tilapia embryos exhibited faster development and subsequently higher long-term growth at 32°C compared to those reared at 28°C and 24°C. Next-generation sequencing data revealed a total of 5,141 unique circRNAs across all temperature groups, of which 1,604, 1,531, and 1,169 circRNAs were exclusively found in the 24°C, 28°C and 32°C groups, respectively. Among them, circNexn exhibited a 1.7-fold (log2) upregulation in the 24°C group and a 1.3-fold (log2) upregulation in the 32°C group when compared to the 28°C group. Conversely, circTTN and circTTN_b were downregulated in the 24°C groups compared to their 28°C and 32°C counterparts. Furthermore, these differentially expressed circRNAs were found to have multiple interactions with myomiRs, highlighting their potential as promising candidates for further investigation in the context of muscle growth plasticity. Taken together, our findings provide new insights into the molecular mechanisms that may underlie muscle growth plasticity in response to thermal variation in fish, with important implications in the context of climate change, fisheries and aquaculture.

miR-460b-5p promotes proliferation and differentiation of chicken myoblasts and targets RBM19 gene

The meat production of broilers is crucial to economic benefits of broiler industries, while the slaughter performance of broilers is directly determined by skeletal muscle development. Hence, the broiler breeding for growth traits shows a great importance. As a kind of small noncoding RNA, microRNA (miRNA) can regulate the expression of multiple genes and perform a wide range of regulation in organisms. Currently, more and more studies have confirmed that miRNAs are closely associated with skeletal muscle development of chickens. Based on our previous miR-seq analysis (accession number: PRJNA668199), miR-460b-5p was screened as one of the key miRNAs probably involved in the growth regulation of chickens. However, the regulatory effect of miR-460b-5p on the development of chicken skeletal muscles is still unclear. Therefore, miR-460b-5p was further used for functional validation at the cellular level in this study. The expression pattern of miR-460b-5p was investigated in proliferation and differentiation stages of chicken primary myoblasts. It was showed that the expression level of miR-460b-5p gradually decreased from the proliferation stage (GM 50%) to the lowest at 24 h of differentiation. As differentiation proceeded, miR-460b-5p expression increased significantly, reaching the highest and stabilizing at 72 h and 96 h of differentiation. Through mRNA quantitative analysis of proliferation marker genes, CCK-8 and Edu assays, miR-460b-5p was found to significantly facilitate the transition of myoblasts from G1 to S phase and promote chicken myoblast proliferation. mRNA and protein quantitative analysis of differentiation marker genes, as well as the indirect immunofluorescence results of myotubes, revealed that miR-460b-5p significantly stimulated myotube development and promote chicken myoblast differentiation. In addition, the target relationship was validated for miR-460b-5p according to the dual-luciferase reporter assay and mRNA quantitative analysis, which indicates that miR-460b-5p was able to regulate RBM19 expression by specifically binding to the 3' UTR of RBM19. In summary, miR-460b-5p has positive regulatory effects on the proliferation and differentiation of chicken myoblasts, and RBM19 is a target gene of miR-460b-5p.

The Discovery and Characterization of Conserved and Novel miRNAs in the Different Developmental Stages and Organs of Pikeperch (Sander lucioperca)

Micro RNAs (miRNAs) are short non-coding RNAs that act as post-transcriptional gene expression regulators. Genes regulated in vertebrates include those affecting growth and development or stress and immune response. Pikeperch (Sander lucioperca) is a species that is increasingly being considered for farming in recirculation aquaculture systems. We characterized the pikeperch miRNA repertoire to increase the knowledge of the genomic mechanisms affecting performance and health traits by applying small RNA sequencing to different developmental stages and organs. There were 234 conserved and 8 novel miRNA genes belonging to 104 families. A total of 375 unique mature miRNAs were processed from these genes. Many mature miRNAs showed high relative abundances or were significantly more expressed at early developmental stages, like the miR-10 and miR-430 family, let-7, the miRNA clusters 106-25-93, and 17-19-92. Several miRNAs associated with immune responses (e.g., slu-mir-731-5p, slu-mir-2188-5p, and slu-mir-8159-5p) were enriched in the spleen. The mature miRNAs slu-mir-203a-3p and slu-mir-205-5p were enriched in gills. These miRNAs are similarly abundant in many vertebrates, indicating that they have shared regulatory functions. There was also a significantly increased expression of the disease-associated miR-462/miR-731 cluster in response to hypoxia stress. This first pikeperch miRNAome reference resource paves the way for future functional studies to identify miRNA-associated variations that can be utilized in marker-assisted breeding programs.

Review: Understanding fish muscle biology in the indeterminate growth species pacu (Piaractus mesopotamicus)

The muscle phenotype of fish is regulated by numerous factors that, although widely explored, still need to be fully understood. In this context, several studies aimed to unravel how internal and external stimuli affect the muscle growth of these vertebrates. The pacu (Piaractus mesopotamicus) is a species of indeterminate muscular growth that quickly reaches high body weight. For this reason, it adds great importance to the productive sector, along with other round fish. In this context, we aimed to compile studies on fish biology and skeletal muscle growth, focusing on studies by our research group that used pacu as an experimental model along with other species. Based on these studies, new muscle phenotype regulators were identified and explored in vivo, in vitro, and in silico studies, which strongly contribute to advances in understanding muscle growth mechanisms with future applications in the productive sector.

MicroRNA qPCR normalization in Nile tilapia (Oreochromis niloticus): Effects of acute cold stress on potential reference targets

The Nile tilapia (Oreochromis niloticus) is one of the most important cultured fish worldwide, but tilapia culture is largely affected by low temperatures. Recent studies suggest that microRNAs (miRNAs) regulate cold tolerance traits in fish. In general, qPCR-based methods are the simplest and most accurate forms of miRNA quantification. However, qPCR data heavily depends on appropriate normalization. Therefore, the aim of the present study is to determine whether the expression of previously tested, stably expressed miRNAs are affected by acute cold stress in Nile tilapia. For this purpose, one small nuclear RNA (U6) and six candidate reference miRNAs (miR-23a, miR-25-3, Let-7a, miR-103, miR-99-5, and miR-455) were evaluated in four tissues (blood, brain, liver, and gills) under two experimental conditions (acute cold stress and control) in O. niloticus. The stability of the expression of each candidate reference miRNA was analyzed by four independent methods (the delta Ct method, geNorm, NormFinder, and BestKeeper). Further, consensual comprehensive ranking of stability was built with RefFinder. Overall, miR-103 was the most stable reference miRNA in this study, and miR-103 and Let-7a were the best combination of reference targets. Equally important, Let-7a, miR-23a, and miR-25-3 remained consistently stable across different tissues and experimental groups. Considering all variables, U6, miR-99-5, and miR-455 were the least stable candidates under acute cold stress. Most important, suitable reference miRNAs were validated in O. niloticus, facilitating further accurate miRNA quantification in this species.

The novel circular RNA CircMef2c is positively associated with muscle growth in Nile tilapia

Muscle growth in teleosts is a complex biological process orchestrated by numerous protein-coding genes and non-coding RNAs. A few recent studies suggest that circRNAs are involved in teleost myogenesis, but the molecular networks involved remain poorly understood. In this study, an integrative omics approach was used to determine myogenic circRNAs in Nile tilapia by quantifying and comparing the expression profile of mRNAs, miRNAs, and circRNAs in fast muscle from full-sib fish with distinct growth rates. There were 1947 mRNAs, 9 miRNAs, and 4 circRNAs differentially expressed between fast- and slow-growing individuals. These miRNAs can regulate myogenic genes and have binding sites for the novel circRNA circMef2c. Our data indicate that circMef2c may interact with three miRNAs and 65 differentially expressed mRNAs to form multiple competing endogenous RNA networks that regulate growth, thus providing novel insights into the role of circRNAs in the regulation of muscle growth in teleosts.

Coding and Noncoding Genes Involved in Atrophy and Compensatory Muscle Growth in Nile Tilapia

Improvements in growth-related traits reduce fish time and production costs to reach market size. Feed deprivation and refeeding cycles have been introduced to maximize aquaculture profits through compensatory growth. However, the molecular compensatory growth signature is still uncertain in Nile tilapia. In this study, fish were subjected to two weeks of fasting followed by two weeks of refeeding. The growth curve in refed tilapia was suggestive of a partial compensatory response. Transcriptome profiling of starved and refed fish was conducted to identify genes regulating muscle atrophy and compensatory growth. Pairwise comparisons revealed 5009 and 478 differentially expressed (differential) transcripts during muscle atrophy and recovery, respectively. Muscle atrophy appears to be mediated by the ubiquitin-proteasome and autophagy/lysosome systems. Autophagy-related 2A, F-box and WD repeat domain containing 7, F-box only protein 32, miR-137, and miR-153 showed exceptional high expression suggesting them as master regulators of muscle atrophy. On the other hand, the muscle compensatory growth response appears to be mediated by the continuous stimulation of muscle hypertrophy which exceeded normal levels found in control fish. For instance, genes promoting ribosome biogenesis or enhancing the efficiency of translational machinery were upregulated in compensatory muscle growth. Additionally, myogenic microRNAs (e.g., miR-1 and miR-206), and hypertrophy-associated microRNAs (e.g., miR-27a-3p, miR-29c, and miR-29c) were reciprocally expressed to favor hypertrophy during muscle recovery. Overall, the present study provided insights into the molecular mechanisms regulating muscle mass in fish. The study pinpoints extensive growth-related gene networks that could be used to inform breeding programs and also serve as valuable genomic resources for future mechanistic studies.

Aquaculture omics: An update on the current status of research and data analysis

Aquaculture is the fast-growing agricultural sector and has the ability to meet the growing demand for protein nutritional security for future population. In future aquaculture is going to be the major source of fish proteins as capture fisheries reached at its maximum. However, several challenges need to overcome such as lack of genetically improved strains/varieties, lack of species-specific feed/functional feed, round the year availability of quality fish seed, pollution of ecosystems and increased frequencies of disease occurrence etc. In recent years, the continuous development of high throughput sequencing technology has revolutionized the biological sciences and provided necessary tools. Application of 'omics' in aquaculture research have been successfully used to resolve several productive and reproductive issues and thus ensure its sustainability and profitability. To date, high quality draft genomes of over fifty fish species have been generated and successfully used to develop large number of single nucleotide polymorphism markers (SNPs), marker panels and other genomic resources etc in several aquaculture species. Similarly, transcriptome profiling and miRNAs analysis have been used in aquaculture research to identify key transcripts and expression analysis of candidate genes/miRNAs involved in reproduction, immunity, growth, development, stress toxicology and disease. Metagenome analysis emerged as a promising scientific tool to analyze the complex genomes contained within microbial communities. Metagenomics has been successfully used in the aquaculture sector to identify novel and potential pathogens, antibiotic resistance genes, microbial roles in microcosms, microbial communities forming biofloc, probiotics etc. In the current review, we discussed application of high-throughput technologies (NGS) in the aquaculture sector.

Prediction and Experimental Validation of a New Salinity-Responsive Cis-Regulatory Element (CRE) in a Tilapia Cell Line

Transcriptional regulation is a major mechanism by which organisms integrate gene x environment interactions. It can be achieved by coordinated interplay between cis-regulatory elements (CREs) and transcription factors (TFs). Euryhaline tilapia (Oreochromis mossambicus) tolerate a wide range of salinity and thus are an appropriate model to examine transcriptional regulatory mechanisms during salinity stress in fish. Quantitative proteomics in combination with the transcription inhibitor actinomycin D revealed 19 proteins that are transcriptionally upregulated by hyperosmolality in tilapia brain (OmB) cells. We searched the extended proximal promoter up to intron1 of each corresponding gene for common motifs using motif discovery tools. The top-ranked motif identified (STREME1) represents a binding site for the Forkhead box TF L1 (FoxL1). STREME1 function during hyperosmolality was experimentally validated by choosing two of the 19 genes, chloride intracellular channel 2 (clic2) and uridine phosphorylase 1 (upp1), that are enriched in STREME1 in their extended promoters. Transcriptional induction of these genes during hyperosmolality requires STREME1, as evidenced by motif mutagenesis. We conclude that STREME1 represents a new functional CRE that contributes to gene x environment interactions during salinity stress in tilapia. Moreover, our results indicate that FoxL1 family TFs are contribute to hyperosmotic induction of genes in euryhaline fish.

Genetic regulation and variation of expression of miRNA and mRNA transcripts in fetal muscle tissue in the context of sex, dam and variable fetal weight

BACKGROUND

Impaired skeletal muscle growth in utero can result in reduced birth weight and pathogenesis of intrauterine growth restriction. Fetal and placental growth is influenced by many factors including genetic, epigenetic and environmental factors. In fact, the sex and genotype of the fetus itself, as well as the mother providing it with a suitable environment, influence the growth of the fetus. Hence, our goal was to decipher and elucidate the molecular pathways of developmental processes mediated by miRNAs and mRNAs in fetal muscle tissue in the context of sex, dam, and fetal weight. Therefore, we analyse the variation of miRNA and mRNA expression in relation to these factors. In addition, the coincidence of genetic regulation of these mRNAs and miRNAs, as revealed by expression quantitative trait loci (eQTL) analyses, with sex-, mother- and weight-associated expression was investigated.

METHODS

A three-generation pig F2 population (n = 118) based on reciprocal crossing of German Landrace (DL) and Pietrain (Pi) was used. Genotype information and transcriptomic data (mRNA and miRNA) from longissimus dorsi muscle (LDM) of pig fetuses sampled at 63 days post-conception (dpc) were used for eQTL analyses.

RESULTS

The transcript abundances of 13, 853, and 275 probe-sets were influenced by sex, dam and fetal weight at 63 dpc, respectively (FDR < 5%). Most of significant transcripts affected by sex were located on the sex chromosomes including KDM6A and ANOS1 or autosomes including ANKS1B, LOC100155138 and miR-153. The fetal muscle transcripts associated with fetal weight indicated clearer metabolic directions than maternally influenced fetal muscle transcripts. Moreover, coincidence of genetic regulation (eQTL) and variation in transcript abundance due to sex, dam and fetal weight were identified.

CONCLUSIONS

Integrating information on eQTL, sex-, dam- and weight-associated differential expression and QTL for fetal weight allowed us to identify molecular pathways and shed light on the basic biological processes associated with differential muscle development in males and females, with implications for adaptive fetal programming.

Amino Acids and IGF1 Regulation of Fish Muscle Growth Revealed by Transcriptome and microRNAome Integrative Analyses of Pacu (Piaractus mesopotamicus) Myotubes

Amino acids (AA) and IGF1 have been demonstrated to play essential roles in protein synthesis and fish muscle growth. The myoblast cell culture is useful for studying muscle regulation, and omics data have contributed enormously to understanding its molecular biology. However, to our knowledge, no study has performed the large-scale sequencing of fish-cultured muscle cells stimulated with pro-growth signals. In this work, we obtained the transcriptome and microRNAome of pacu (Piaractus mesopotamicus)-cultured myotubes treated with AA or IGF1. We identified 1228 and 534 genes differentially expressed by AA and IGF1. An enrichment analysis showed that AA treatment induced chromosomal changes, mitosis, and muscle differentiation, while IGF1 modulated IGF/PI3K signaling, metabolic alteration, and matrix structure. In addition, potential molecular markers were similarly modulated by both treatments. Muscle-miRNAs (miR-1, -133, -206 and -499) were up-regulated, especially in AA samples, and we identified molecular networks with omics integration. Two pairs of genes and miRNAs demonstrated a high-level relationship, and involvement in myogenesis and muscle growth: marcksb and miR-29b in AA, and mmp14b and miR-338-5p in IGF1. Our work helps to elucidate fish muscle physiology and metabolism, highlights potential molecular markers, and creates a perspective for improvements in aquaculture and in in vitro meat production.

Identification of Genes Related to Cold Tolerance and Novel Genetic Markers for Molecular Breeding in Taiwan Tilapia (Oreochromis spp.) via Transcriptome Analysis

Taiwan tilapia is one of the primary species used in aquaculture practices in Taiwan. However, as a tropical fish, it is sensitive to cold temperatures that can lead to high mortality rates during winter months. Genetic and broodstock management strategies using marker-assisted selection and breeding are the best tools currently available to improve seed varieties for tilapia species. The purpose of this study was to develop molecular markers for cold stress-related genes using digital gene expression analysis of next-generation transcriptome sequencing in Taiwan tilapia (Oreochromis spp.). We constructed and sequenced cDNA libraries from the brain, gill, liver, and muscle tissues of cold-tolerance (CT) and cold-sensitivity (CS) strains. Approximately 35,214,833,100 nucleotides of raw sequencing reads were generated, and these were assembled into 128,147 unigenes possessing a total length of 185,382,926 bp and an average length of 1446 bp. A total of 25,844 unigenes were annotated using five protein databases and Venny analysis, and 38,377 simple sequence repeats (SSRs) and 65,527 single nucleotide polymorphisms (SNPs) were identified. Furthermore, from the 38-cold tolerance-related genes that were identified using differential gene expression analysis in the four tissues, 13 microsatellites and 37 single nucleotide polymorphism markers were identified. The results of the genotype analysis revealed that the selected markers could be used for population genetics. In addition to the diversity assessment, one of the SNP markers was determined to be significantly related to cold-tolerance traits and could be used as a molecular marker to assist in the selection and verification of cold-tolerant populations. The specific genetic markers explored in this study can be used for the identification of genetic polymorphisms and cold tolerance traits in Taiwan tilapia, and they can also be used to further explore the physiological and biochemical molecular regulation pathways of fish that are involved in their tolerance to environmental temperature stress.

Variant expression signatures of microRNAs and protein related to growth in a crossbreed between two strains of Nile tilapia (Oreochromis niloticus)

Nile tilapia (Oreochromis niloticus) is a species of worldwide importance for aquaculture. A crossbred lineage was developed through introgressive backcross breeding techniques and combines the high growth performance of the Chitralada (CHIT) lwith attractive reddish color of the Red Stirling (REDS) strains. Since the crossbreed has an unknown genetically improved background, the objective of this work was to characterize expression signatures that portray the advantageous phenotype of the crossbreeds. We characterized the microRNA transcriptome by high throughput sequencing (RNA-seq) and the proteome through mass spectrometry (ESI-Q-TOF-MS) and applied bioinformatics for the comparative analysis of such molecular data on the three strains. Crossbreed expressed a distinct set of miRNAs and proteins compared to the parents. They comprised several microRNAs regulate traits of economic interest. Proteomic profiles revealed differences between parental and crossbreed in expression of proteins associated with glycolisis. Distinctive miRNA and protein signatures contribute to the phenotype of crossbreed.

miR-152 targets pyruvate kinase to regulate the glycolytic activity of pig skeletal muscles and affects pork quality

As a type of non-coding RNA, microRNAs are widely involved in the biological processes of animals. In the present study, the expression of miR-152 in glycolytic muscle fibers (Longissimus thoracis, LT) was lower than that of oxidative muscle fibers (Psoas major, PM). Using dual luciferase assay, miR-152 was shown to target muscle pyruvate kinase (PKM) to perform biological functions. Moreover, overexpression of miR-152 in primary porcine cells inhibited PKM gene expression and reduced lactic acid production in cells, whereas inhibition of miR-152 expression promoted PKM gene expression and increased lactic acid production. Correlation analysis showed that the expression of miR-152 was significantly positively correlated with the ultimate pH of LT after slaughter, while the expression of the PKM gene was significantly negatively correlated with the final pH of LT. In vivo and in vitro experiments discussed herein suggest that miR-152 may affect muscle pH by targeting the expression of the PKM gene. Our findings enrich the understanding of the genetic regulatory network that influences pork quality.

MicroRNA roles in regeneration: Multiple lessons from zebrafish

MicroRNAs (miRNAs) are small noncoding RNAs with pivotal roles in the control of gene expression. By comparing the miRNA profiles of uninjured vs. regenerating tissues and structures, several studies have found that miRNAs are potentially involved in the regenerative process. By inducing miRNA overexpression or inhibition, elegant experiments have directed regenerative responses validating relevant miRNA-to-target interactions. The zebrafish (Danio rerio) has been the epicenter of regenerative research because of its exceptional capability to self-repair damaged tissues and body structures. In this review, we discuss recent discoveries that have improved our understanding of the impact of gene regulation mediated by miRNAs in the context of the regeneration of fins, heart, retina, and nervous tissue in zebrafish. We compiled what is known about the miRNA control of regeneration in these tissues and investigated the links among up-regulated and down-regulated miRNAs, their putative or validated targets, and the regenerative process. Finally, we briefly discuss the forthcoming prospects, highlighting directions and the potential for further development of this field.

In Utero Fetal Weight in Pigs Is Regulated by microRNAs and Their Target Genes

Impaired skeletal muscle growth in utero can result in reduced birth weight and poor carcass quality in pigs. Recently, we showed the role of microRNAs (miRNAs) and their target genes in prenatal skeletal muscle development and pathogenesis of intrauterine growth restriction (IUGR). In this study, we performed an integrative miRNA-mRNA transcriptomic analysis in longissimus dorsi muscle (LDM) of pig fetuses at 63 days post conception (dpc) to identify miRNAs and genes correlated to fetal weight. We found 13 miRNAs in LDM significantly correlated to fetal weight, including miR-140, miR-186, miR-101, miR-15, miR-24, miR-29, miR-449, miR-27, miR-142, miR-99, miR-181, miR-199, and miR-210. The expression of these miRNAs decreased with an increase in fetal weight. We also identified 1315 genes significantly correlated to fetal weight at 63 dpc, of which 135 genes were negatively correlated as well as identified as potential targets of the above-listed 13 miRNAs. These miRNAs and their target genes enriched pathways and biological processes important for fetal growth, development, and metabolism. These results indicate that the transcriptomic profile of skeletal muscle can be used to predict fetal weight, and miRNAs correlated to fetal weight can serve as potential biomarkers of prenatal fetal health and growth.

Expression profile of miRNAs involved in the hepatoprotective effects of curcumin against oxidative stress in Nile tilapia

Curcumin is a polyphenol with antioxidant activity that has been used to protect the health of fish livers. Our previous studies about comparative transcriptome have shown that curcumin can enhance the Nrf2-Keap1 signaling pathway and induce downstream anti-stress genes to maintain cell viability. However, the possible role of miRNAs in the protective mechanism of curcumin is not understood. In this study, the tilapia hepatocyte H₂O₂ stress model was used, and the miRNA expression profile for four groups (control group, curcumin group, H₂O₂ group, and protection group) were established by high-throughput sequencing. In our results, 278-333 types of Oreochromis niloticus miRNAs, 309-543 types of conserved miRNAs, and 535-746 types of novel miRNAs were identified in different samples. Differentially expressed miRNAs were identified by comparing miRNA expression profiles among the four groups. The expression levels were confirmed by q-PCR. The target genes of these differentially expressed miRNAs were predicted, and their functional annotations were enriched by GO and KEGG analysis, which revealed that many target genes were involved in "response to stimulus" and "antioxidant activity" in each pair of groups. Several miRNAs related to oxidative stress showed differential expression. For example, in the H₂O₂ group, the expression of miR-122 was decreased, and the expression of miR-21 and miR-489 increased significantly. In the curcumin group, the expression of miR-153b was decreased, and the expression of miR-200a and miR-29 was increased significantly. miR-153b, miR-200a, and miR-29 may be involved in the regulation of the Nrf2-Keap1 signaling pathway by curcumin. This work might provide insights into the molecular mechanisms of miRNA regulation of curcumin on the prevention and alleviation of liver diseases in fish.

Comparative transcriptome profiles of large and small bodied large-scale loaches cultivated in paddy fields

Fish culture in paddy fields is a traditional aquaculture mode, which has a long history in East Asia. Large-scale loach (Paramisgurnus dabryanus) fast growth is suitable for paddy fields aquaculture in China. The objective of this study was to identify differential expression genes (DEGs) in the brain, liver and muscle tissues between large (LG, top 5% of maximum total length) and small (SG, top 5% of minimum total length) groups using RNA-seq. In total, 150 fish were collected each week and 450 fish were collected at twelfth week from three paddy fields for all the experimental. Histological observation found that the muscle fibre diameter of LG loaches was greater than that of SG loaches. Transcriptome results revealed that the high expression genes (HEGs) in LG loaches (fold change ≥ 2, p < 0.05) were mainly concentrated in metabolic pathways, such as "Thyroid hormone signalling pathway", "Citrate cycle (TCA cycle)", "Carbon metabolism", "Fatty acid metabolism", and "Cholesterol metabolism", and the HEGs in SG loaches were enriched in the pathways related to environmental information processing such as "Cell adhesion molecules (CAMs)", "ECM- receptor interaction" and "Rap1 signalling pathway"; cellular processes such as "Tight junction", "Focal adhesion", "Phagosome" and "Adherens junction". Furthermore, IGFs gene family may play an important role in loach growth for their different expression pattern between the two groups. These findings can enhance our understanding about the molecular mechanism of different growth and development levels of loaches in paddy fields.

Identification of differentially expressed miRNAs in serum extracellular vesicles (EVs) of Kazakh sheep at early pregnancy

MiRNAs-containing extracellular vesicles (EVs) possess the unique function of mediating intercellular communication and participating in many biological processes such as post-transcriptional gene regulation of embryo implantation and placental development. In the present study, Illumina small-RNA sequencing was used to identify differentially expressed (DE) miRNAs in serum EVs of pregnant (P) and non-pregnant (NP) Kazakh sheep at Day 17 from mating. The specifically and differentially expressed miRNAs at early pregnancy in sheep were verified by using RT-PCR. The target genes of DE miRNAs were predicted by bioinformatics software, and the functional and pathway enrichment analysis was performed on Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) terms. A total of 562 miRNAs (210 novel miRNAs) were identified by sequencing, of which 57 miRNAs were differentially expressed, 49 were up-regulated, 8 were down-regulated and 22 novel miRNAs were specifically expressed in the pregnant sheep. Eight highly expressed known miRNA (miR-378-3p, miR-320-3p, miR-22-3p, let-7b, miR-423-3p, miR-221, miR-296-3p, miR-147-3p) in pregnant group were down-regulated in the control group. miRNAs-containing pregnancy-related terms and regulatory pathways regulation were enriched using both GO and KEGG analyses. Moreover, we also envisioned a miRNA-mRNA interaction network to understand the function of miRNAs involved in the early pregnancy serum regulatory network. The results of RT-PCR verification confirmed the reliability of small-RNA sequencing. Among them, miR-22-3p and miR-378-3p were significantly differentially expressed (DE) between pregnant sheep and non-pregnant group (p <  0.01). The site at which oar-miR-22-3p binds MAPK3 was determined with a dual-luciferase system. This is the first integrated analysis of the expression profiles of EV-miRNAs and their targets during early pregnancy in ewes. These data identify key miRNAs that influence the implantation of sheep in the early stage of pregnancy, and provide theoretical basis for further molecular regulatory mechanisms research.

Roles of adiponectin and leptin signaling-related microRNAs in the preventive effects of calorie restriction in mammary tumor development

Calorie restriction (CR) is suggested to prevent the development of mammary tumors (MTs); however, the mechanism remains to be clarified. We aimed to determine the microRNA (miRNA) profile in mice applied to 2 different CR protocols; chronic (CCR) and intermittent (ICR) and follow the MT development. In addition, the roles of miRNAs involved in adiponectin and/or leptin signaling pathways were investigated. Mice were divided into 3 groups: ad-libitum (AL), CCR, or ICR, which comprised 3 weeks of AL feeding followed by 1 week of 60% CR in a cyclic manner. Blood and tissue collection were performed at weeks 10, 17/18, 49/50 and 81/82. Long-term CCR provided better protection compared with ICR for MT development with a delay in the MT occurrence. Adiponectin expression in mammary fat pad were significantly higher in CCR group compared with AL. Using GeneChip Array, 250 of 3195 miRNAs were differentially expressed among the dietary groups. Thirteen of 250 miRNAs were related to adiponectin and/or leptin signaling genes. Results were verified by reverse transcription polymerase chain reaction. Specifically, miR-326-3p, miR-500-3p and miR-129-5p, which are related to adiponectin and/or leptin signaling, may play important roles in the preventive effects of CR in MT development and in ageing. Thus, these miRNAs might be putative biomarkers to target for diagnostic and treatment purposes. Novelty: Type of CR and micro RNA interaction is related to ageing. miR-326-3p, miR-500-3p and miR-129-5p expression levels were differentially expressed in MT development and in ageing. The genes associated with adiponectin and/or leptin signaling pathways are regulated by certain miRNAs in the protective effects of CR.

Molecular Mechanisms Regulating Muscle Plasticity in Fish

Growth rates in fish are largely dependent on genetic and environmental factors, of which the latter can be highly variable throughout development. For this reason, muscle growth in fish is particularly dynamic as muscle structure and function can be altered by environmental conditions, a concept referred to as muscle plasticity. Myogenic regulatory factors (MRFs) like Myogenin, MyoD, and Pax7 control the myogenic mechanisms regulating quiescent muscle cell maintenance, proliferation, and differentiation, critical processes central for muscle plasticity. This review focuses on recent advancements in molecular mechanisms involving microRNAs (miRNAs) and DNA methylation that regulate the expression and activity of MRFs in fish. Findings provide overwhelming support that these mechanisms are significant regulators of muscle plasticity, particularly in response to environmental factors like temperature and nutritional challenges. Genetic variation in DNA methylation and miRNA expression also correlate with variation in body weight and growth, suggesting that genetic markers related to these mechanisms may be useful for genomic selection strategies. Collectively, this knowledge improves the understanding of mechanisms regulating muscle plasticity and can contribute to the development of husbandry and breeding strategies that improve growth performance and the ability of the fish to respond to environmental challenges.

MiR-2014-5p and miR-1231-5p regulate muscle growth of Larimichthys crocea by targeting MSTN gene

MicroRNAs (miRNAs) play an important role in regulating gene expression, and myostatin (MSTN) has been widely recognized as a key gene for muscle growth and development. Through high-throughput sequencing to study the effects of starvation on miRNA transcriptomes in Larimichthys crocea muscle tissue, we found that the expression of miR-2014, miR-1231 and miR-1470 were significantly different between fasting and normal feeding Larimichthys crocea. Bioinformatics analysis predicted that miR-2014, miR-1231 and miR-1470 target MSTN mRNA 3'UTR. To verify the accuracy of predictions, we constructed double luciferase plasmids containing MSTN 3'UTR and confirmed that miR-2014-5p and miR-1231-5p can inhibit MSTN expression by targeting MSTN 3'UTR using double luciferase experiments, while miR-1470 is not involved in regulation. Subsequent site-directed mutation experiments reflected the specificity of the target sequence. In addition, quantitative PCR experiments revealed that miR-2014-5p and miR-1231-5p may participate in the regulation of MSTN expression in fasting and refeeding period, respectively. These results implied that miRNA may take part in muscle growth regulation during starvation. It provides some insights into the molecular regulation mechanism of MSTN in response to starvation stress in fish.

Meta-analysis of differentially-regulated hepatic microRNAs identifies candidate post-transcriptional regulation networks of intermediary metabolism in rainbow trout

MicroRNAs (miRNAs) are small non-coding RNAs which act as post-transcriptional regulators by decreasing targeted mRNA translation and stability. Principally targeting small 3' UTR elements of protein-coding mRNAs through complementary base-pairing, miRNAs are promiscuous regulators of the transcriptome. While potent roles for hepatic miRNAs in the regulation of energy metabolism have emerged in rodent models, comparative roles in other vertebrates remain largely unexplored. Indeed, while several miRNAs are deeply conserved among vertebrates, the acquisition of lineage- and species-specific miRNAs, as well as the rewiring between miRNA-mRNA target relationships beg the question of regulatory and functional conservation and innovation of miRNAs and their targets involved in energy metabolism. Here we provide a meta-analysis of differentially expressed hepatic miRNAs in rainbow trout, a scientifically and economically important teleost species with a 'glucose-intolerant' phenotype. Following exposure to nutritional and social context-dependent metabolic challenges, we analyzed differential miRNA expression from small-RNA-sequencing datasets generated with a consistent bioinformatics pipeline in conjunction with an in silico target prediction of metabolic transcripts and pathways. We provide evidence for evolutionary conserved (let-7, miRNA-27 family) and rewired (miRNA-30 family, miRNA-152, miRNA-722) miRNA-metabolic target gene networks in the context of the salmonid genome. These findings represent important first steps in our understanding of the comparative regulation and function of hepatic miRNAs in rainbow trout energy metabolism. We propose that the identified miRNA families should be prioritized for future comparative functional investigation in the context of hepatic energy- and glucose metabolism in rainbow trout.

Genomics to accelerate genetic improvement in tilapia

Selective breeding of tilapia populations started in the early 1990s and over the past three decades tilapia has become one of the most important farmed freshwater species, being produced in more than 125 countries around the globe. Although genome assemblies have been available since 2011, most of the tilapia industry still depends on classical selection techniques using mass spawning or pedigree information to select for growth traits with reported genetic gains of up to 20% per generation. The involvement of international breeding companies and research institutions has resulted in the rapid development and application of genomic resources in the last few years. GWAS and genomic selection are expected to contribute to uncovering the genetic variants involved in economically relevant traits and increasing the genetic gain in selective breeding programs, respectively. Developments over the next few years will probably focus on achieving a deep understanding of genetic architecture of complex traits, as well as accelerating genetic progress in the selection for growth-, quality- and robustness-related traits. Novel phenotyping technologies (i.e. phenomics), lower-cost whole-genome sequencing approaches, functional genomics and gene editing tools will be crucial in future developments for the improvement of tilapia aquaculture.

miRTil: An Extensive Repository for Nile Tilapia microRNA Next Generation Sequencing Data

Nile tilapia is the third most cultivated fish worldwide and a novel model species for evolutionary studies. Aiming to improve productivity and contribute to the selection of traits of economic impact, biotechnological approaches have been intensively applied to species enhancement. In this sense, recent studies have focused on the multiple roles played by microRNAs (miRNAs) in the post-transcriptional regulation of protein-coding genes involved in the emergence of phenotypes with relevance for aquaculture. However, there is still a growing demand for a reference resource dedicated to integrating Nile Tilapia miRNA information, obtained from both experimental and in silico approaches, and facilitating the analysis and interpretation of RNA sequencing data. Here, we present an open repository dedicated to Nile Tilapia miRNAs: the “miRTil database”. The database stores data on 734 mature miRNAs identified in 11 distinct tissues and five key developmental stages. The database provides detailed information about miRNA structure, genomic context, predicted targets, expression profiles, and relative 5p/3p arm usage. Additionally, miRTil also includes a comprehensive pre-computed miRNA-target interaction network containing 4936 targets and 19,580 interactions.

In-particle stem-loop RT-qPCR for specific and multiplex microRNA profiling

Here we report a novel method of microRNA (miRNA) profiling with particle-based multiplex quantitative reverse transcription polymerase chain reaction (RT-qPCR). To achieve target-specific reaction in a particle, the stem-loop RT primer and forward primer for each target miRNA were chemically immobilized to the particle. Target-specific cDNA synthesis proceeds with the stem-loop RT primer and then qPCR subsequently proceeds with the forward primer to rapidly achieve a quantitative result. High-fidelity multiplex assay was also accomplished in a single PCR process by loading multiple particles for each specific miRNA. The method for primer supply in the particles, involving confinement of the target-specific RT and PCR primers in the matrix of particles, led to the reduction of nonspecific reactions and improved the selectivity of the miRNA assay while minimizing labor in a multiple target assay. Specifically, this particle-based assay enabled the differentiation of mature miRNA from precursor with selectivity of 270:1 in terms of amplification speed. This advanced method also showed good discrimination among highly homologous let-7 family members, with cross-reaction rates of less than 5%. We demonstrated a very simple process of five-plex miRNA profiling in total RNA, and the measured changes in expression level were consistent with those from a conventional singleplex method.

Expression profiles of growth-related genes in two Nile tilapia strains and their crossbred provide insights into introgressive breeding effects

The Nile tilapia (Oreochromis niloticus) is a prominent farmed fish in aquaculture worldwide. Crossbreeding has recently been carried out between the Red-Stirling and the wt Chitralada strains of Nile tilapia, producing a heterotic hybrid (7/8 Chitralada and 1/8 Red-Stirling) that combines the superior growth performance of the Chitralada with the reddish coloration of the Red-Stirling strain. While classical selective breeding and crossbreeding strategies are well known, the molecular mechanisms underlying the phenotypic expression of economically advantageous traits in tilapia remain largely unknown. Molecular investigations have shown that variable expression of growth hormone (gh), insulin-like growth factors (igf1 and 2) and somatolactin (smtla) - components of the growth hormone/insulin-like growth factor (GH/IGF) axis - and myostatin (mstn) genes can affect traits of economic relevance in farmed animals. The aim of this study was to assess and compare the gene expression signature among Chitralada, Red-Stirling and their backcross hybrid in order to gain insights into the effects of introgressive breeding in modulation of the GH/IGF axis. Gene expression analyses in distinct tissues showed that most genes of the GH/IGF axis were up-regulated and mstn was down-regulated in backcross animals in comparison with Red-Stirling and Chitralada animals. These gene expression profiles revealed that backcross animals displayed a distinctive expression signature, which attests to the effectiveness of the introgressive breeding technique. Our findings also suggest that the GH/IGF axis and mstn genes might be candidate markers for fish performance and prove useful within genetic improvement programs aimed at the production of superior-quality tilapia strains using introgressive breeding.

Effects of food deprivation on plasma insulin-like growth factor-1 (Igf1) and Igf binding protein (Igfbp) gene transcription in juvenile cabezon (Scorpaenichthys marmoratus)

The growth hormone (GH)/insulin-like growth factor (Igf) endocrine axis regulates somatic growth in the face of changing environmental conditions. In actinopterygian fishes, food availability is a key modulator of the somatotropic axis, with lower food intake generally depressing liver Igf1 release to diminish growth. Igf1 signaling, however, also involves several distinct IGF binding proteins (Igfbps), and the functional roles of many of these Igfbps in affecting growth during shifting food availability remain uncertain. Here, we tested how complete food deprivation (fasting) affected gene transcription for paralogs of all six types of Igfbps in the liver and fast-twitch skeletal muscle of cabezon (Scorpaenichthys marmoratus), a nearshore marine fish important for recreational fisheries in the eastern North Pacific Ocean. Juvenile cabezon were maintained as either fed (6% mass food⋅g fish wet mass^-1⋅d^-1) or fasted for 14 d. Fasted fish exhibited a lower body condition (K), a depressed mass-specific growth rate (SGR), and reduced plasma concentrations of Igf1. In the liver, fasting reduced the relative abundance of gene transcripts encoding Igfbps igfbp2a and igfbp2b, while significantly elevating mRNA levels for igfbp1a, igfbp1b, igfbp3b, and igfbp4. Fasting also reduced hepatic mRNA levels of GH receptor-1 (ghr1) - but not GH receptor-2 (ghr2) - supporting the idea that changes in liver sensitivity to GH may underlie the decline in plasma Igf1 during food deprivation. In skeletal muscle, fasting downregulated gene transcripts encoding igf1, igfbp2b, igfbp5b, and igfbp6b, while also upregulating mRNAs for igf2 and ghr2. These data demonstrate isoform-specific regulation of Igfbps in liver and skeletal muscle in cabezon experiencing food deprivation and reinforce the idea that the repertoire of duplicated Igfbp genes that evolved in actinopterygian fishes supports a diverse scope of endocrine and paracrine functions.

Differences in the fast muscle methylome provide insight into sex-specific epigenetic regulation of growth in Nile tilapia during early stages of domestication

Growth is a complex trait whose variability within a population cannot be explained solely by genetic variation. Epigenetic regulation is often suggested as an important factor shaping the phenotype, but its association with growth can be highly context- and species-dependent. Nevertheless, the mechanisms involved in epigenetic regulation of growth in fish are poorly understood. We have used reduced representation bisulphite sequencing to determine the genome-wide CpG methylation patterns in male and female Nile tilapia of different sizes but at the same early stage of domestication. The average CpG methylation level in the reduced genome representation was 63% across groups but many sites displayed group-specific methylation patterns. The number of differentially methylated (DM) CpGs was much higher when the interaction between sex and weight was included rather than when these factors were considered separately. There were 1128 DM CpGs between large and small females and 970 DM CpGs between large and small males. We have found many growth-related genes associated with DM CpGs, namely map3k5 and akt3 in females and gadd45g and ppargc1a in males. Only 5% of CpG locations associated with growth were common to both sexes. In particular, the autophagy-related gene atg14 displayed a high association of methylation with growth exclusively in males. The sexually dimorphic association between atg14 methylation and growth may uncover novel metabolic mechanisms at play during mouth brooding in Nile tilapia females. Taken together, our data suggest that epigenetic regulation of growth in Nile tilapia involves different gene networks in males and females.

Comparison of microRNAs in adipose and muscle tissue from seven indigenous Chinese breeds and Yorkshire pigs

Elucidation of the pig microRNAome is essential for interpreting functional elements of the genome and understanding the genetic architecture of complex traits. Here, we extracted small RNAs from skeletal muscle and adipose tissue, and we compared their expression levels between one Western breed (Yorkshire) and seven indigenous Chinese breeds. We detected the expression of 172 known porcine microRNAs (miRNAs) and 181 novel miRNAs. Differential expression analysis found 92 and 12 differentially expressed miRNAs in adipose and muscle tissue respectively. We found that different Chinese breeds shared common directional miRNA expression changes compared to Yorkshire pigs. Some miRNAs differentially expressed across multiple Chinese breeds, including ssc-miR-129-5p, ssc-miR-30 and ssc-miR-150, are involved in adipose tissue function. Functional enrichment analysis revealed that the target genes of the differentially expressed miRNAs are associated mainly with signaling pathways rather than metabolic and biosynthetic processes. The miRNA-target gene and miRNA-phenotypic traits networks identified many hub miRNAs that regulate a large number of target genes or phenotypic traits. Specifically, we found that intramuscular fat content is regulated by the greatest number of miRNAs in muscle tissue. This study provides valuable new candidate miRNAs that will aid in the improvement of meat quality and production.

Identification and characterization of microRNAs in the gonads of Crassostrea hongkongensis using high-throughput sequencing

Crassostrea hongkongensis is one of the three most-commonly cultivated oyster species in China. Although microRNAs (miRNAs) expression in the gonads have been widely investigated, few studies of miRNAs in mollusk gonads are available, particularly in oyster. In the present study, we analyzed the miRNAs expressed in the ovaries and testes of C. hongkongensis. We obtained 14,166,409 and 15,133,900 raw reads from the ovaries and testes, respectively, yielding 13,634,997 (ovarian) and 14,494,149 (testicular) 18-35-nt sequences. We mapped these sequences to the C. hongkongensis genome reference sequence, and identified 8,771,717 (ovarian) and 9,926,014 (testicular) sequences corresponding to miRNAs in the Rfam database. After blasting the miRNA sequences against the miRBase database, we identified 50 known mature miRNAs and 53 novel miRNAs. Of these, 27 miRNAs were significantly upregulated in ovaries as compared to the testes, and 43 miRNAs were significantly upregulated in the testes as compared to the ovaries. To validate the differential expression results generated by Illumina sequencing, we used RT-real-time quantitative PCR (RT-qPCR) to characterize the expression patterns of the six most differently expressed miRNAs (lgi-miR-1990, lgi-miR-1986, lgi-miR-263b, lgi-miR-279, lgi-miR-1992, and novel_98) as well as two miRNAs associated with gonad development (lgi-miR-29 and lgi-miR-8). Most of the RT-qPCR miRNA expression patterns were similar to those recovered by high-throughput sequencing with the exceptions of novel_98 and lgi-miR-1992. Gene Ontology (GO) annotations indicated that the multi-organism cellular process GO category was enriched with the target genes of the differentially expressed miRNAs. Additionally, the target genes were enriched in several KEGG pathways, including the ECM-receptor interaction, galactose metabolism, phagosome, and notch signaling pathway. These pathways are involved in gonadal differentiation and the maintenance of gonad function. This identification and characterization of the miRNAs differentially expressed between the ovaries and testes of C. hongkongensis will increase our understanding of the role of miRNAs in gonad differentiation in the oyster.

MiR-34b-5p Mediates the Proliferation and Differentiation of Myoblasts by Targeting IGFBP2

As key post-transcriptional regulators, microRNAs (miRNAs) play an indispensable role in skeletal muscle development. Our previous study suggested that miR-34b-5p and IGFBP2 could have a potential role in skeletal muscle growth. Our goal in this study is to explore the function and regulatory mechanism of miR-34b-5p and IGFBP2 in myogenesis. In this study, the dual-luciferase reporter assay and Western blot analysis showed that IGFBP2 is a direct target of miR-34b-5p. Flow cytometric analysis and EdU assay showed that miR-34b-5p could repress the cell cycle progression of myoblasts, and miR-34b-5p could promote the formation of myotubes by promoting the expression of MyHC. On the contrary, the overexpression of IGFBP2 significantly facilitated the proliferation of myoblasts and hampered the formation of myotubes. Together, our results indicate that miR-34b-5p could mediate the proliferation and differentiation of myoblasts by targeting IGFBP2.

Identification and characteristics of muscle growth-related microRNA in the Pacific abalone, Haliotis discus hannai

BACKGROUND

The Pacific abalone, Haliotis discus hannai, is the most important cultivated abalone in China. Improving abalone muscle growth and increasing the rate of growth are important genetic improvement programs in this industry. MicroRNAs are important small noncoding RNA molecules that regulate post-transcription gene expression. However, no miRNAs have been reported to regulate muscle growth in H. discus hannai.

RESULTS

we profiled six small RNA libraries for three large abalone individuals (L_HD group) and three small individuals (S_HD group) using RNA sequencing technology. A total of 205 miRNAs, including 200 novel and 5 known miRNAs, were identified. In the L_HD group, 3 miRNAs were up-regulated and 7 were down-regulated compared to the S_HD specimens. Bioinformatics analysis of miRNA target genes revealed that miRNAs participated in the regulation of cellular metabolic processes, the regulation of biological processes, the Wnt signaling pathway, ECM-receptor interaction, and the MAPK signaling pathway, which are associated with regulating growth. Bone morphogenetic protein 7 (BMP7) was verified as a target gene of hdh-miR-1984 by a luciferase reporter assay and we examined the expression pattern in different developmental stages.

CONCLUSION

This is the first study to demonstrate that miRNAs are related to the muscle growth of H. discus hannai. This information could be used to study the mechanisms of abalone muscle growth. These DE-miRNAs may be useful as molecular markers for functional genomics and breeding research in abalone and closely related species.

Genome-wide identification and differentially expression analysis of lncRNAs in tilapia

Background

Long noncoding RNAs (LncRNAs) play important roles in fundamental biological processes. However, knowledge about the genome-wide distribution and stress-related expression of lncRNAs in tilapia is still limited.

Results

Genome-wide identification of lncRNAs in the tilapia genome was carried out in this study using bioinformatics tools. 103 RNAseq datasets that generated in our laboratory or collected from NCBI database were analyzed. In total, 72,276 high-confidence lncRNAs were identified. The averaged positive correlation coefficient (r_mean = 0.286) between overlapped lncRNA and mRNA pairs showed significant differences with the values for all lncRNA-mRNA pairs (r_mean = 0.176, z statistics = − 2.45, p value = 0.00071) and mRNA-mRNA pairs (r_mean = 0.186, z statistics = − 2.23, p value = 0.0129). Weighted correlation network analysis of the lncRNA and mRNA datasets from 12 tissues identified 21 modules and many interesting mRNA genes that clustered with lncRNAs. Overrepresentation test indicated that these mRNAs enriched in many biological processes, such as meiosis (p = 0.00164), DNA replication (p = 0.00246), metabolic process (p = 0.000838) and in molecular function, e.g., helicase activity (p = 0.000102) and catalytic activity (p = 0.0000612). Differential expression (DE) analysis identified 99 stress-related lncRNA genes and 1955 tissue-specific DE lncRNA genes. MiRNA-lncRNA interaction analysis detected 72,267 lncRNAs containing motifs with sequence complementary to 458 miRNAs.

Conclusions

This study provides an invaluable resource for further studies on molecular bases of lncRNAs in tilapia genomes. Further function analysis of the lncRNAs will help to elucidate their roles in regulating stress-related adaptation in tilapia.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5115-x) contains supplementary material, which is available to authorized users.

Integrated analysis of lncRNA and mRNA expression in rainbow trout families showing variation in muscle growth and fillet quality traits

Muscle yield and quality traits are important for the aquaculture industry and consumers. Genetic selection for these traits is difficult because they are polygenic and result from multifactorial interactions. To study the genetic architecture of these traits, phenotypic characterization of whole body weight (WBW), muscle yield, fat content, shear force and whiteness were measured in ~500 fish representing 98 families from a growth-selected line. RNA-Seq was used to sequence the muscle transcriptome of different families exhibiting divergent phenotypes for each trait. We have identified 240 and 1,280 differentially expressed (DE) protein-coding genes and long noncoding RNAs (lncRNAs), respectively, in fish families exhibiting contrasting phenotypes. Expression of many DE lncRNAs (n = 229) was positively correlated with overlapping, neighboring or distantly located protein-coding genes (n = 1,030), resulting in 3,392 interactions. Three DE antisense lncRNAs were co-expressed with sense genes known to impact muscle quality traits. Forty-four DE lncRNAs had potential sponge functions to miRNAs that affect muscle quality traits. This study (1) defines muscle quality associated protein-coding and noncoding genes and (2) provides insight into non-coding RNAs involvement in regulating growth and fillet quality traits in rainbow trout.

miRNAs regulate acute transcriptional changes in broiler embryos in response to modification of incubation temperature

MicroRNAs are post-transcriptional regulators that play critical roles in diverse biological processes. We hypothesize that miRNAs may be involved in regulating transcriptome responses to changes in embryonic incubation temperature in chickens affecting differentiation and proliferation processes during tissue development. Therefore, we conducted comparative transcriptome profiling of miRNAs to examine altered expression in breast and hind muscle of embryos and day 35 chickens experiencing high (38.8 °C), control (37.8 °C), or low (36.8 °C) embryonic incubation temperature during embryonic day (ED) 7–10 or ED10–13. The results revealed differential expression of miRNAs due to modification of embryonic incubation temperature in a muscle type-specific and a developmental stage-specific manner. The immediate effects of thermal change observed in embryos were substantial compared to the subtle long-term effects in chickens at day 35 post-hatch. Upregulation of miR-133 in breast muscle and downregulation of miR-199a-5p, miR-1915, and miR-638 in hind muscle post ED7–10 high-temperature treatment are functionally associated with myogenesis and body size. ED10–13 low-temperature treatment led to downregulation of let-7, miR-93, and miR-130c that are related to proliferation and differentiation. The results provide insight into the dynamics of miRNA expression at variable embryonic incubation temperatures during developmental processes and indicate a major regulatory role of miRNAs in acute responses to modified environmental conditions that affect remodelling of cells and tissues.

Genome-wide microRNA screening in Nile tilapia reveals pervasive isomiRs' transcription, sex-biased arm switching and increasing complexity of expression throughout development

MicroRNAs (miRNAs) are key regulators of gene expression in multicellular organisms. The elucidation of miRNA function and evolution depends on the identification and characterization of miRNA repertoire of strategic organisms, as the fast-evolving cichlid fishes. Using RNA-seq and comparative genomics we carried out an in-depth report of miRNAs in Nile tilapia (Oreochromis niloticus), an emergent model organism to investigate evo-devo mechanisms. Five hundred known miRNAs and almost one hundred putative novel vertebrate miRNAs have been identified, many of which seem to be teleost-specific, cichlid-specific or tilapia-specific. Abundant miRNA isoforms (isomiRs) were identified with modifications in both 5p and 3p miRNA transcripts. Changes in arm usage (arm switching) of nine miRNAs were detected in early development, adult stage and even between male and female samples. We found an increasing complexity of miRNA expression during ontogenetic development, revealing a remarkable synchronism between the rate of new miRNAs recruitment and morphological changes. Overall, our results enlarge vertebrate miRNA collection and reveal a notable differential ratio of miRNA arms and isoforms influenced by sex and developmental life stage, providing a better picture of the evolutionary and spatiotemporal dynamics of miRNAs.

Fishing Into the MicroRNA Transcriptome

In the last decade, several studies have been focused on revealing the microRNA (miRNA) repertoire and determining their functions in farm animals such as poultry, pigs, cattle, and fish. These small non-protein coding RNA molecules (18-25 nucleotides) are capable of controlling gene expression by binding to messenger RNA (mRNA) targets, thus interfering in the final protein output. MiRNAs have been recognized as the main regulators of biological features of economic interest, including body growth, muscle development, fat deposition, and immunology, among other highly valuable traits, in aquatic livestock. Currently, the miRNA repertoire of some farmed fish species has been identified and characterized, bringing insights about miRNA functions, and novel perspectives for improving health and productivity. In this review, we summarize the current advances in miRNA research by examining available data on Neotropical and other key species exploited by fisheries and in aquaculture worldwide and discuss how future studies on Neotropical fish could benefit from this knowledge. We also make a horizontal comparison of major results and discuss forefront strategies for miRNA manipulation in aquaculture focusing on forward-looking ideas for forthcoming research.

Identifying a Major QTL Associated with Salinity Tolerance in Nile Tilapia Using QTL-Seq

Selection of new lines with high salinity tolerance allows for economically feasible production of tilapias in brackish water areas. Mapping QTLs and identifying the markers linked to salinity-tolerant traits are the first steps in the improvement of the tolerance in tilapia through marker-assisted selection techniques. By using QTL-seq strategy and linkage-based analysis, two significant QTL intervals (chrLG4 and chrLG18) on salinity-tolerant traits were firstly identified in the Nile tilapia. Fine mapping with microsatellite and SNP markers suggested a major QTL region that located at 23.0 Mb of chrLG18 and explained 79% of phenotypic variation with a LOD value of 95. Expression analysis indicated that at least 10 genes (e.g., LACTB2, KINH, NCOA2, DIP2C, LARP4B, PEX5R, and KCNJ9) near or within the QTL interval were significantly differentially expressed in intestines, brains, or gills under 10, 15, or 20 ppt challenges. Our findings suggest that QTL-seq can be effectively utilized in QTL mapping of salinity-tolerant traits in fish. The identified major QTL is a promising locus to improve our knowledge on the genetic mechanism of salinity tolerance in tilapia.

Proteomic and microRNA Transcriptome Analysis revealed the microRNA-SmyD1 network regulation in Skeletal Muscle Fibers performance of Chinese perch

Fish myotomes are comprised of anatomically segregated fast and slow muscle fibers that possess different metabolic and contractile properties. Although the expression profile properties in fast and slow muscle fibers had been investigated at the mRNA levels, a comprehensive analysis at proteomic and microRNA transcriptomic levels is limited. In the present study, we first systematically compared the proteomic and microRNA transcriptome of the slow and fast muscles of Chinese perch (Siniperca chuatsi). Total of 2102 proteins were identified in muscle tissues. Among them, 99 proteins were differentially up-regulated and 400 were down-regulated in the fast muscle compared with slow muscle. MiRNA microarrays revealed that 199 miRNAs identified in the two types of muscle fibers. Compared with the fast muscle, the 32 miRNAs was up-regulated and 27 down-regulated in the slow muscle. Specifically, expression of miR-103 and miR-144 was negatively correlated with SmyD1a and SmyD1b expression in fast and slow muscles, respectively. The luciferase reporter assay further verified that the miR-103 and miR-144 directly regulated the SmyD1a and SmyD1b expression by targeting their 3′-UTR. The constructed miRNA-SmyD1 interaction network might play an important role in controlling the development and performance of different muscle fiber types in Chinese perch.

Genome-Wide QTL Analysis Identified Significant Associations Between Hypoxia Tolerance and Mutations in the GPR132 and ABCG4 Genes in Nile Tilapia

Exposure to hypoxia induces both acute and chronic stress responses, which plays an important role in health of cultured organisms including growth, reproduction, immunity, and other energy demanding activities. Application of advanced genomic technologies allows rapid identification of hypoxia trait-associated genes and precise selection of superior brood stocks with high tolerance in tilapia. By applying QTL-seq and double-digest restriction-site associated DNA sequencing (ddRAD-seq) techniques, we identified four genome-wide significant quantitative trait loci (QTLs) for hypoxia tolerance and many suggestive QTLs in Nile tilapia. These QTLs explained 6.6-14.7% of the phenotypic variance. Further analysis revealed that single nucleotide polymorphisms (SNPs) in exons of both GPR132 and ABCG4 genes located in genome-wide QTL intervals were significantly associated with hypoxia-tolerant traits. Expression analysis of both genes suggested that they were strong candidate genes involved into hypoxia tolerance in tilapia. Our findings suggest that both QTL-seq and ddRAD-seq techniques can be effectively utilized in QTL mapping of hypoxia traits in fish. Our data supply a basis for further marker-assisted selection of super lines with a high level of tolerance against low oxygen stress in the tilapia.

Transcriptome analysis of the sea cucumber (Apostichopus japonicus) with variation in individual growth

The sea cucumber (Apostichopus japonicus) is an economically important aquaculture species in China. However, the serious individual growth variation often caused financial losses to farmers and the genetic mechanisms are poorly understood. In the present study, the extensively analysis at the transcriptome level for individual growth variation in sea cucumber was carried out. A total of 118946 unigenes were assembled from 255861 transcripts, with N50 of 1700. Of all unigenes, about 23% were identified with at least one significant match to known databases. In all four pair of comparison, 1840 genes were found to be expressed differently. Global hypometabolism was found to be occurred in the slow growing population, based on which the hypothesis was raised that growth retardation in individual growth variation of sea cucumber is one type of dormancy which is used to be against to adverse circumstances. Besides, the pathways such as ECM-receptor interaction and focal adhesion were enriched in the maintenance of cell and tissue structure and communication. Further, 76645 SSRs, 765242 SNPs and 146886 ins-dels were detected in the current study providing an extensive set of data for future studies of genetic mapping and selective breeding. In summary, these results will provides deep insight into the molecular basis of individual growth variation in marine invertebrates, and be valuable for understanding the physiological differences of growth process.

Pou1f1, the key transcription factor related to somatic growth in tilapia (Orechromis niloticus), is regulated by two independent post-transcriptional regulation mechanisms

This study aims to determine the post-transcriptional regulation mechanism of the transcription factor pou1f1 (pou class 1 homeobox 1), which is the key gene for pituitary development, somatic growth in vertebrates, and transcription of several hormone genes in teleost fish. MicroRNA miR-223-3p was identified as a bona fide target of pou1f; overexpression of miR-223-3p in primary pituitary cells led to the down-regulation of pou1f1 and downstream genes, and inhibition of miR-223-3p led to the up-regulation of pou1f1 in Nile tilapia dispersed primary pituitary cells. An adenylate-uridylate-rich element (AU-Rich element) was found in the 3'UTR of pou1f1 mRNA, and deletion of the AU-Rich element led to slower mRNA decay and therefore more protein output. A potential mutual relationship between miR-223-3p and the AU-rich element was also investigated, and the results demonstrated that with or without the AU-Rich element, miR-223-3p induced the up-regulation of a reporter system under serum starvation conditions, indicating that miR-223-3p and the AU-Rich element function independent of each other. This study is the first to investigate the post-transcriptional mechanism of pou1f1, which revealed that miR-223-3p down-regulated pou1f1 and downstream gene expressions, and the AU-Rich element led to rapid decay of pou1f1 mRNA. MicroRNA miR-223-3p and the AU-Rich element co-regulated the post-transcriptional expression of pou1f1 independently in Nile tilapia, demonstrating that pou1f1 is under the control of a dual post-transcription regulation mechanism.

IGF-1/IGF-1R/hsa-let-7c axis regulates the committed differentiation of stem cells from apical papilla

Insulin-like growth factor-1 (IGF-1) and its receptor IGF-1R play a paramount role in tooth/bone formation while hsa-let-7c actively participates in the osteogenic differentiation of mesenchymal stem cells. However, the interaction between IGF-1/IGF-1R and hsa-let-7c on the committed differentiation of stem cells from apical papilla (SCAPs) remains unclear. In this study, human SCAPs were isolated and treated with IGF-1 and hsa-let-7c over/low-expression viruses. The odonto/osteogenic differentiation of these stem cells and the involvement of mitogen-activated protein kinase (MAPK) pathway were subsequently investigated. Alizarin red staining showed that hsa-let-7c low-expression can significantly promote the mineralization of IGF-1 treated SCAPs, while hsa-let-7c over-expression can decrease the calcium deposition of IGF-1 treated SCAPs. Western blot assay and real-time reverse transcription polymerase chain reaction further demonstrated that the expression of odonto/osteogenic markers (ALP, RUNX2/RUNX2, OSX/OSX, OCN/OCN, COL-I/COL-I, DSPP/DSP, and DMP-1/DMP-1) in IGF-1 treated SCAPs were significantly upregulated in Let-7c-low group. On the contrary, hsa-let-7c over-expression could downregulate the expression of these odonto/osteogenic markers. Moreover, western blot assay showed that the JNK and p38 MAPK signaling pathways were activated in Let-7c-low SCAPs but inhibited in Let-7c-over SCAPs. Together, the IGF-1/IGF-1R/hsa-let-7c axis can control the odonto/osteogenic differentiation of IGF-1-treated SCAPs via the regulation of JNK and p38 MAPK signaling pathways.

Estrogen regulation of gene expression in the teleost fish immune system

Elucidating the mechanisms of estrogens-induced immunomodulation in teleost fish is of great importance due to the observed worldwide continuing decrease in pristine environments. However, little is know about the immunotoxicological consequences of exposure to these chemicals in fish, or of the mechanisms through which these effects are mediated. In this review, we summarize the results showing estrogens (natural or synthetic) acting through estrogen receptors and regulating specific target genes, also through microRNAs (miRNAs), leading to modulation of the immune functioning. The identification and characterization of miRNAs will provide new opportunities for functional genome research on teleost immune system and can also be useful when screening for novel molecule biomarkers for environmental pollution.