Identification and characterization of 63 MicroRNAs in the Asian seabass Lates calcarifer

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.

Identification of mRNA-miRNA-lncRNA regulatory network associated with the immune response to Aeromonas salmonicides infection in the black rockfish (Sebastes schlegelii)

In aquaculture, Aeromonas salmonicides (A. salmonicida) is a main fish pathogen because of its nearly worldwide distribution, and broad host range. Recently, an increasing number of evidences have uncovered the roles of mRNA-miRNA-lncRNA network in fish diseases. In current study, RNA-seq was conducted in the black rockfish spleen following A. salmonicida infection at 0 h (Sp0 or control) and three different post-infection time-points (2 h: Sp2, 12 h: Sp12 and 24 h: Sp24, respectively) to comprehensively identify differentially expressed (DE) mRNAs, miRNAs and lncRNAs. Enrichment analysis and protein-protein interaction (PPI) analysis of DE mRNAs were performed. Then, expression and correlation analysis for mRNAs and their upstream miRNAs and lncRNAs were conducted. Finally, a total of 1364 mRNAs, 17 miRNAs and 1584 lncRNAs exhibited significantly differential expressions during bacterial infection in the black rockfish spleen. Functional enrichment analysis suggested that they were significantly enriched in several immune-related pathways, including Amino sugar and nucleotide sugar metabolism, Cell adhesion molecules (CAMs), Neuroactive ligand-receptor interaction, Nicotinate and nicotinamide metabolism, Pentose and glucuronate interconversions, Phagosome, Proteasome, etc. Subsequently, 1091 lncRNA-miRNA-mRNA pathways (323 in Sp2, 609 in Sp12 and 207 in Sp24) were constructed including 400 lncRNAs, 69 miRNAs, and 70 mRNAs. Meanwhile, NLRC3/novel-264/LNC_00116154 pathway demonstrated important immune modulating function in the black rockfish against A. salmonicida infection. Finally, the novel mRNA-miRNA-lncRNA sub-networks were established, among which all mRNAs and ncRNAs possessed significant predictive values for further studies for immune responses in the black rockfish.

Characterization of miRNAs in Embryonic, Larval, and Adult Lumpfish Provides a Reference miRNAome for Cyclopterus lumpus

Simple Summary

Lumpfish (Cyclopterus lumpus) is an emergent aquaculture species, and its miRNA repertoire is still unknown. miRNAs are critical post-transcriptional modulators of teleost gene expression. Therefore, a lumpfish reference miRNAome was characterized by small RNA sequencing and miRDeep analysis of samples from different organs and developmental stages. The resulting miRNAome, an essential reference for future expression analyses, consists of 443 unique mature miRNAs from 391 conserved and eight novel miRNA genes. Enrichment of specific miRNAs in particular organs and developmental stages indicates that some conserved lumpfish miRNAs regulate organ and developmental stage-specific functions reported in other teleosts.

Abstract

MicroRNAs (miRNAs) are endogenous small RNA molecules involved in the post-transcriptional regulation of protein expression by binding to the mRNA of target genes. They are key regulators in teleost development, maintenance of tissue-specific functions, and immune responses. Lumpfish (Cyclopterus lumpus) is becoming an emergent aquaculture species as it has been utilized as a cleaner fish to biocontrol sea lice (e.g., Lepeophtheirus salmonis) infestation in the Atlantic Salmon (Salmo salar) aquaculture. The lumpfish miRNAs repertoire is unknown. This study identified and characterized miRNA encoding genes in lumpfish from three developmental stages (adult, embryos, and larvae). A total of 16 samples from six different adult lumpfish organs (spleen, liver, head kidney, brain, muscle, and gill), embryos, and larvae were individually small RNA sequenced. Altogether, 391 conserved miRNA precursor sequences (discovered in the majority of teleost fish species reported in miRbase), eight novel miRNA precursor sequences (so far only discovered in lumpfish), and 443 unique mature miRNAs were identified. Transcriptomics analysis suggested organ-specific and age-specific expression of miRNAs (e.g., miR-122-1-5p specific of the liver). Most of the miRNAs found in lumpfish are conserved in teleost and higher vertebrates, suggesting an essential and common role across teleost and higher vertebrates. This study is the first miRNA characterization of lumpfish that provides the reference miRNAome for future functional studies.

MicroRNAs and Related Cytokine Factors Quickly Respond in the Immune Response of Channel Catfish to Lipopolysaccharides and β-Glucan Stimulation

MicroRNAs (miRNAs) are well-known as powerful regulators of gene expression, with their potential to serve for immunology widely researched in mammals and birds but rarely in fishes. To better understand fish immunology behavior, we herein investigated nine immune-related miRNAs that were reported in other animals, as well as five related cytokine factors and lysozyme (LZM) in the liver, anterior kidney, and spleen of Channel Catfish Ictalurus punctatus after being stimulated by lipopolysaccharides (LPS) and β-glucan. We also predicated the potential targets of these miRNAs via bioinformatics and further investigated nine of them via quantitative real-time PCR. Results showed that expressions of the nine miRNAs were quickly changed in varying extent after stimulation by LPS, especially for miR-122, miR-142a, miR-155, and miR-223, which were significantly changed in spleen, and the same occurred for the LZM and three cytokine factors TNF-α, IFN-γ and TLR2. Compared with LPS, although most of the miRNAs and the cytokine genes were also affected by β-glucan, the extent of the effect was weak. Bioinformatics analysis revealed many immune-related targets of the miRNAs, with some of them reported by previous studies. For the nine investigated target genes, seven targets (77.8%) were significantly upregulated after the stimulation of LPS. It therefore can be inferred that the immune-related miRNAs, LZM, and cytokine factors elicited quick immune responses of Channel Catfish to LPS stimulation as in other animals, but the regulation mechanism of miRNAs might be complex and diverse. This research will contribute to a better understanding will support further immunology research in fishes.

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.

The emerging regulatory roles of noncoding RNAs in immune function of fish: MicroRNAs versus long noncoding RNAs

The genome could be considered as raw data expressed in proteins and various types of noncoding RNAs (ncRNAs). However, a large portion of the genome is dedicated to ncRNAs, which in turn represent a considerable amount of the transcriptome. ncRNAs are modulated on levels of type and amount whenever any physiological process occurs or as a response to external modulators. ncRNAs, typically forming complexes with other partners, are key molecules that influence diverse cellular processes. Based on the knowledge of mammalian biology, ncRNAs are known to regulate and control diverse trafficking pathways and cellular activities. Long noncoding RNAs (lncRNAs) notably have diverse and more regulatory roles than microRNAs. Expanding these studies on fish has derived the same conclusion with relevance to other species, including invertebrates, explored the potentials to harness such types of RNA to further understand the biology of such organisms, and opened gates for applying recent technologies, such as RNA interference and delivering micromolecules as microRNAs to living cells and possibly to target organs. These technologies should improve aquaculture productivity and fish health, as well as help understand fish biology.

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.

Kisspeptin Influences the Reproductive Axis and Circulating Levels of microRNAs in Senegalese Sole

Kisspeptin regulates puberty and reproduction onset, acting upstream of the brain–pituitary–gonad (HPG) axis. This study aimed to test a kisspeptin-based hormonal therapy on cultured Senegalese sole (G1) breeders, known to have reproductive dysfunctions. A single intramuscular injection of KISS2-10 decapeptide (250 µg/kg) was tested in females and males during the reproductive season, and gonad maturation, sperm motility, plasma levels of gonadotropins (Fsh and Lh) and sex steroids (11-ketotestosterone, testosterone and estradiol), as well as changes in small non-coding RNAs (sncRNAs) in plasma, were investigated. Fsh, Lh, and testosterone levels increased after kisspeptin injection in both sexes, while sperm analysis did not show differences between groups. Let7e, miR-199a-3p and miR-100-5p were differentially expressed in females, while miR-1-3p miRNA was up-regulated in kisspeptin-treated males. In silico prediction of mRNAs targeted by miRNAs revealed that kisspeptin treatment might affect paracellular transporters, regulate structural and functional polarity of cells, neural networks and intracellular trafficking in Senegalese sole females; also, DNA methylation and sphingolipid metabolism might be altered in kisspeptin-treated males. Results demonstrated that kisspeptin stimulated gonadotropin and testosterone secretion in both sexes and induced an unanticipated alteration of plasma miRNAs, opening new research venues to understand how this neuropeptide impacts in fish HPG axis.

The miR-155 regulates cytokines expression by SOSC1 signal pathways of fish in vitro and in vivo

MiR-155 is reported as immune regulated miRNA in mammalian corresponding to immunity, antibacterial and antiviral effects regulation. However, the roles and mechanisms of the miRNA have remained largely undefined. We herein comprehensively investigated the functions of miR-155 in vitro and in vivo by miR-155 mimics, agomir and antagomir in Cyprinus carpio and Ictalurus punctatus, with the target genes in the SOSC1 pathway certified in I. punctatus via luciferase reporter assays. Results showed that the miR-155 regulated the expressions of cytokines, including TNF-α, IFN-γ, IL-1β, IL-6 and IL-10. Further research confirmed SOSC1 as one of the targets of the miRNA, and the JAK1/STAT3/SOSC1 signal pathway involved in the miR-155 effects on the expression of immune cytokines as well. Additionally, the changes of TLR2 in fish may also be related to miR-155 along with its target SOCS1, and the TLR2/MyD88 pathway may partly participate in the effects of the miR-155 on the cytokines. The research here confirmed that the miR-155 can regulate cytokines expression by SOSC1 signal pathways of fish in vitro and in vivo, which would provide resources for understanding and studying about immune regulation in fish.

Dietary Immunostimulant CpG Modulates MicroRNA Biomarkers Associated with Immune Responses in Atlantic Salmon (Salmo salar)

MicroRNAs (miRNAs) are key regulators in fish immune responses. However, no study has previously characterized the impact of polyriboinosinic polyribocytidylic acid (pIC) and formalin-killed typical Aeromonas salmonicida (ASAL) on miRNA expression in Atlantic salmon fed a commercial diet with and without immunostimulant CpG. To this end, first, we performed small RNA deep sequencing and qPCR analyses to identify and confirm pIC- and/or ASAL-responsive miRNAs in the head kidney of salmon fed a control diet. DESeq2 analyses identified 12 and 18 miRNAs differentially expressed in pIC and ASAL groups, respectively, compared to the controls. Fifteen of these miRNAs were studied by qPCR; nine remained significant by qPCR. Five miRNAs (miR-27d-1-2-5p, miR-29b-2-5p, miR-146a-5p, miR-146a-1-2-3p, miR-221-5p) were shown by qPCR to be significantly induced by both pIC and ASAL. Second, the effect of CpG-containing functional feed on miRNA expression was investigated by qPCR. In pre-injection samples, 6 of 15 miRNAs (e.g., miR-181a-5-3p, miR-462a-3p, miR-722-3p) had significantly lower expression in fish fed CpG diet than control diet. In contrast, several miRNAs (e.g., miR-146a-1-2-3p, miR-192a-5p, miR-194a-5p) in the PBS- and ASAL-injected groups had significantly higher expression in CpG-fed fish. Multivariate statistical analyses confirmed that the CpG diet had a greater impact on miRNA expression in ASAL-injected compared with pIC-injected fish. This study identified immune-relevant miRNA biomarkers that will be valuable in the development of diets to combat infectious diseases of salmon.

Discovery of MicroRNAs from Batrachuperus yenyuanensis Using Deep Sequencing and Prediction of Their Targets

MicroRNAs (miRNAs), a family of ∼22-nucleotide non-coding single-stranded RNA molecules, are considered as key post-transcriptional regulators of gene expression that regulate various biological processes in living organism. Many miRNAs have been identified in animals; however, few have been reported in Hynobiidae species. The present study is aimed to identify a full repertoire of miRNAs in Batrachuperus yenyuanensis (Yenyuan stream salamander), which would significantly increase our knowledge of miRNAs in amphibians. A small RNA library was constructed from B. yenyuanensis and sequenced using deep sequencing. As a result, 1,717,751 clean reads were obtained, representing 356 known and 80 novel miRNAs. Additionally, expression levels of eight randomly selected miRNAs in B. yenyuanensis were confirmed using the stem-loop quantitative real-time reverse transcription PCR. In addition, 13,972 targets were predicted for these identified miRNAs, although the physiological functions of many of these targets remain unknown. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis suggested that the predicted targets are involved in a variety of physiological regulatory functions in B. yenyuanensis. These results provide useful information for further research on the miRNAs involved in the growth and development of B. yenyuanensis, as well as adaptation of this species to its high-altitude habitats.

Low expression levels of plasma miR-141 are associated with susceptibility to gastric cancer

MicroRNAs (miRNAs/miRs) offer great potential as biomarkers for the early detection and prognosis of cancer, and the discovery of miRNAs associated with gastric cancer is required. In the present study, the differences in the plasma expression levels of miR-141 between patients with gastric cancer and healthy controls, and the role of miR-141 in gastric cancer cell oncogenesis were investigated. A follow-up study of 164 patients with gastric cancer who underwent tumor resection was conducted, and comparisons with healthy control subjects were drawn. To investigate the biological functions of miR-141, a series of in vitro and in vivo assays were conducted, including proliferation, wound-healing and Transwell assays, and a xenograft tumor model. The results demonstrated that miR-141 expression was significantly decreased in tumor tissues compared with in healthy tissues (P<0.05). Kaplan-Meier analysis revealed improved survival benefits with increased miR-141 expression (as determined using the log-rank test, P<0.001), and multivariate Cox regression analysis revealed that patients with decreased expression levels of miR-141 carried a greater risk of death (hazard ratio=2.352; 95% CI=1.379-4.012; P=0.002). The downregulation of miR-141 was also associated with WHO staging, particularly for lymph node and distant metastasis. Exogenous overexpression of miR-141 significantly inhibited the proliferative and migratory abilities of the gastric cancer cell line BGC-823. In vivo studies also demonstrated that exogenous overexpression of miR-141 in BGC-823 cells markedly reduced tumor growth in nude mice. The present study revealed that increased miR-141 expression may be a positive prognostic factor, which may be clinically beneficial in patients with gastric cancer.

Lipopolysaccharides-mediated injury to chondrogenic ATDC5 cells can be relieved by Sinomenine via downregulating microRNA-192

Sinomenine (SIN) is an isoquinoline derived from Caulis Sinomenii that has been used to treat rheumatoid arthritis and osteoarthritis for several decades in China. This study aims to reveal the effects of SIN on mouse chondrogenic ATDC5 cells growth and inflammation. SIN was used to treat ATDC5 cells injured by lipopolysaccharides (LPS). The following parameters were determined for evaluating the treatment effects of SIN: cell viability, apoptosis, reactive oxygen species generation, and pro-inflammatory cytokines release. Besides, the expression of LPS-sensitive miRNA (miR-192) and the activation of NF-κB and MAPK signaling were studied to explain SIN's function. SIN with concentration of 30 μM significantly attenuated LPS-induced cell damage via increasing cell viability, inhibiting apoptosis and reactive oxygen species generation, and declining IL-6 and TNF-α release. miR-192 was downregulated by SIN treatment. Restoration of miR-192 expression by miRNA transfection could significantly impede SIN's protective action. Besides, the inhibitory effects of SIN on the activation of NF-κB and MAPK signaling were attenuated by miR-192 overexpression. Furthermore, GDF11 was found to be a target gene of miR-192. LPS-mediated injury to chondrogenic ATDC5 cells can be relieved by SIN via downregulating miR-192 and subsequently repressing the activation of NF-κB and MAPK signaling.

Expanding the miRNA Repertoire in Atlantic Salmon; Discovery of IsomiRs and miRNAs Highly Expressed in Different Tissues and Developmental Stages

MicroRNAs (miRNAs) are important post-transcriptional gene expression regulators. Here, 448 different miRNA genes, including 17 novel miRNAs, encoding for 589 mature Atlantic salmon miRNAs were identified after sequencing 111 samples (fry, pathogen challenged fry, various developmental and adult tissues). This increased the reference miRNAome with almost one hundred genes. Prior to isomiR characterization (mature miRNA variants), the proportion of erroneous sequence variants (ESVs) arising in the analysis pipeline was assessed. The ESVs were biased towards 5’ and 3’ end of reads in unexpectedly high proportions indicating that measurements of ESVs rather than Phred score should be used to avoid misinterpreting ESVs as isomiRs. Forty-three isomiRs were subsequently discovered. The biological effect of the isomiRs measured as increases in target diversity was small (<3%). Five miRNA genes showed allelic variation that had a large impact on target gene diversity if present in the seed. Twenty-one miRNAs were ubiquitously expressed while 31 miRNAs showed predominant expression in one or few tissues, indicating housekeeping or tissue specific functions, respectively. The miR-10 family, known to target Hox genes, were highly expressed in the developmental stages. The proportion of miR-430 family members, participating in maternal RNA clearance, was high at the earliest developmental stage.

Identification and expression analysis of miRNA in hybrid snakehead by deep sequencing approach and their targets prediction

MicroRNAs (miRNAs) play important regulatory roles in numerous biological processes, but there is no report on miRNAs of hybrid snakehead. In this study, four independent small RNA libraries were constructed from the spleen, liver kidney and muscle of hybrid snakehead. These libraries were sequenced using deep sequencing technology, as result, a total of 1,067,172, 1,152,002, 1,653,941 and 970,866 clean reads from these four libraries were obtained. 252 known miRNAs and 63 putative novel miRNAs in these small RNA dataset were identified. The stem-loop RT-qPCR analysis indicated that eight known miRNAs and two novel miRNAs show different expression in eight different kinds of tissues. For better understanding the functions of miRNAs, 95,947 predicated target genes were analyzed by GO and their pathways, the results indicated that these targets of the identified miRNAs are involved in a broad range of physiological functions.

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.

Differential miRNA expression in the three-spined stickleback, response to environmental changes

miRNAs play important role in the various physiological and evolutionary processes, however, there is no data allowing comparison of evolutionary differences between various ecotypes adapted to different environmental conditions and specimen demonstrating immediate physiological response to the environmental changes. We compared miRNA expression profiles between marine and freshwater stickleback populations of the three-spined stickleback to identify the evolutionary differences. To study the immediate physiological response to foreign environment, we explored the changes induced by transfer of marine sticklebacks into freshwater environment and vice versa. Comparative analysis of changes in miRNA expression suggested that they are driven by three independent factors: (1) non-specific changes in miRNA expression under different environmental conditions; (2) specific response to freshwater conditions in the marine stickleback ecotype; (3) specific response to extreme osmotic conditions for both marine and freshwater ecotypes during the contact with non-native environment. Gene Ontology enrichment analysis of differential expressed miRNA targets supports our current hypothesis.

Identification and characterization of microRNAs expressed in the African malaria vector Anopheles funestus life stages using high throughput sequencing

BACKGROUND

Over the past several years, thousands of microRNAs (miRNAs) have been identified in the genomes of various insects through cloning and sequencing or even by computational prediction. However, the number of miRNAs identified in anopheline species is low and little is known about their role. The mosquito Anopheles funestus is one of the dominant malaria vectors in Africa, which infects and kills millions of people every year. Therefore, small RNA molecules isolated from the four life stages (eggs, larvae, pupae and unfed adult females) of An. funestus were sequenced using next generation sequencing technology.

RESULTS

High throughput sequencing of four replicates in combination with computational analysis identified 107 mature miRNA sequences expressed in the An. funestus mosquito. These include 20 novel miRNAs without sequence identity in any organism and eight miRNAs not previously reported in the Anopheles genus but are known in non-anopheles mosquitoes. Finally, the changes in the expression of miRNAs during the mosquito development were determined and the analysis showed that many miRNAs have stage-specific expression, and are co-transcribed and co-regulated during development.

CONCLUSIONS

This study presents the first direct experimental evidence of miRNAs in An. funestus and the first profiling study of miRNA associated with the maturation in this mosquito. Overall, the results indicate that miRNAs play important roles during the growth and development. Silencing such molecules in a specific life stage could decrease the vector population and therefore interrupt malaria transmission.

Identification of novel microRNA genes in freshwater and marine ecotypes of the three-spined stickleback (Gasterosteus aculeatus)

The three-spined stickleback (Gasterosteus aculeatus L.) is an important model organism for studying the molecular mechanisms of speciation and adaptation to salinity. Despite increased interest to microRNA discovery and recent publication on microRNA prediction in the three-spined stickleback using bioinformatics approaches, there is still a lack of experimental support for these data. In this paper, high-throughput sequencing technology was applied to identify microRNA genes in gills of the three-spined stickleback. In total, 595 miRNA genes were discovered; half of them were predicted in previous computational studies and were confirmed here as microRNAs expressed in gill tissue. Moreover, 298 novel microRNA genes were identified. The presence of miRNA genes in selected 'divergence islands' was analysed and 10 miRNA genes were identified as not randomly located in 'divergence islands'. Regulatory regions of miRNA genes were found enriched with selective SNPs that may play a role in freshwater adaptation.

Identification and profiling of growth-related microRNAs of the swimming crab Portunus trituberculatus by using Solexa deep sequencing

MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs that regulate gene expression by post-transcriptional repression of mRNAs. The swimming crab Portunus trituberculatus is one of the most important crustacean species for aquaculture in China. However, to date no miRNAs have been reported to for modulating growth in P. trituberculatus. To investigate miRNAs involved in the growth of this species, we constructed six small RNA libraries for big individuals (BIs) and small individuals (SIs) from a highly inbred family. Six mixed RNA pools of five tissues (eyestalk, gill, heart, hepatopancreas, and muscle) were obtained. By aligning sequencing data with those for known miRNAs, a total of 404 miRNAs, including 339 known and 65 novel miRNAs, were identified from the six libraries. MiR-100 and miR-276a-3p were among the most prominent miRNA species. We identified seven differentially expressed miRNAs between the BIs and SIs, which were validated using real-time PCR. Preliminary analyzes of their putative target genes and GO and KEGG pathway analyzes showed that these differentially expressed miRNAs could play important roles in global transcriptional depression and cell differentiation of P. trituberculatus. This study reveals the first miRNA profile related to the body growth of P. trituberculatus, which would be particularly useful for crab breeding programs.

Identification and characterization of immune-related microRNAs in blunt snout bream, Megalobrama amblycephala

MicroRNAs (miRNAs) play vital roles in diverse biological processes, including in immune response. Blunt snout bream (Megalobrama amblycephala) is a prevalent and important commercial endemic freshwater fish species in China's intensive polyculture systems. To identify immune-related miRNAs of M. amblycephala, two small RNA (sRNA) libraries from immune tissues with or without lipopolysaccharide (LPS) stimulation were constructed and sequenced using the high-throughput sequencing technology. Totally, 16,425,543 and 15,076,813 raw reads, corresponding to 14,156,755 and 13,445,869 clean reads, were obtained in the normal and infected libraries, respectively. A total of 324 miRNAs, including 218 known miRNAs and 106 putative novel miRNAs were identified by bioinformatic analysis. We analyzed differentially expressed miRNAs between two libraries using pairwise comparison. 113 (34.88%) miRNAs were found to be significantly differentially expressed between two libraries, with 63 (55.75%) exhibiting elevated expression in LPS stimulation sample. Thereinto, a number of known miRNAs were identified immune-related. Real-time quantitative PCR (RT-qPCR) were implemented for 12 miRNAs of two samples, and agreement was confirmed between the sequencing and RT-qPCR data. Target genes likely regulated by these differentially expressed miRNAs were predicted using computational prediction. The functional annotation of target genes by Gene Ontology enrichment (GO) and Kyoto Encyclopedia of Genes and Genomes pathway analysis (KEGG) indicated that a majority of differential miRNAs might involved in immune response. To our knowledge, this is the first comprehensive study of miRNAs in response to LPS stimulation in M. amblycephala, even in fish. These results deepened our understanding of the role of miRNAs in the intricate host's immune system, and should be useful to develop new control strategies for host immune defense against various bacterial invasions in M. amblycephala.

Micromanaging metabolism-a role for miRNAs in teleost energy metabolism

MicroRNAs (miRNAs) are small, non-protein coding RNA sequences, which are found in most eukaryotes. Since their initial discovery, miRNAs have emerged as important regulators of many biological processes. One of the most important processes profoundly regulated by miRNAs is energy metabolism. Traditionally, metabolic functions of miRNAs have been studied in genome-sequenced mammalian organisms, especially the mouse model. However, partially driven by commercial interest in aquaculture, increasingly feasible large-scale molecular techniques have resulted in the characterization of miRNA repertoires, and importantly, several genome sequences of several (commercially important) teleost species, which also hold important roles as research models in the comparative physiology of energy metabolism. This review aims to introduce the recent advances in miRNA research in teleost fish and to describe the current knowledge of miRNA function in teleost energy metabolism. The most pressing research needs and questions to determine metabolic roles of miRNAs in teleost models are presented, as well as applicable technical approaches and current bottlenecks. Rainbow trout, which possess the advantages of newly available molecular tools and a long history as comparative research model in teleost energy metabolism, are discussed as a promising research model to address these questions.

Dynamics of gene expression patterns during early development of the European seabass (Dicentrarchus labrax)

Larval and embryonic stages are the most critical period in the life cycle of marine fish. Key developmental events occur early in development and are influenced by external parameters like stress, temperature, salinity, and photoperiodism. Any failure may cause malformations, developmental delays, poor growth, and massive mortalities. Advanced understanding of molecular processes underlying marine larval development may lead to superior larval rearing conditions. Today, the new sequencing and bioinformatic methods allow transcriptome screens comprising messenger (mRNA) and microRNA (miRNA) with the scope of detecting differential expression for any species of interest. In the present study, we applied Illumina technology to investigate the transcriptome of early developmental stages of the European seabass ( Dicentrarchus labrax). The European seabass, in its natural environment, is a euryhaline species and has shown high adaptation processes in early life phases. During its embryonic and larval phases the European seabass lives in a marine environment and as a juvenile it migrates to coastal zones, estuaries, and lagoons. Investigating the dynamics of gene expression in its early development may shed light on factors promoting phenotypic plasticity and may also contribute to the improvement and advancement of rearing methods of the European seabass, a species of high economic importance in European and Mediterranean aquaculture. We present the identification, characterization, and expression of mRNA and miRNA, comprising paralogous genes and differentially spliced transcripts from early developmental stages of the European seabass. We further investigated the detection of possible interactions of miRNA with mRNA.

Temperature during early development has long-term effects on microRNA expression in Atlantic cod

Background

Environmental temperature has serious implications in life cycle of aquatic ectotherms. Understanding the molecular mechanisms of temperature acclimation and adaptation of marine organisms is of the uttermost importance for ecology, fisheries, and aquaculture, as it allows modeling the effects of global warming on population dynamics. Regulatory molecules are major modulators of acclimation and adaptation; among them, microRNAs (miRNAs) are versatile and substantial contributors to regulatory networks of development and adaptive plasticity. However, their role in thermal plasticity is poorly known. We have asked whether the temperature and its shift during the early ontogeny (embryonic and larval development) affect the miRNA repertoire of Atlantic cod (Gadus morhua), and if thermal experience has long-term consequences in the miRNA profile.

Results

We characterized miRNA during different developmental stages and in juvenile tissues using next generation sequencing. We identified 389 putative miRNA precursor loci, 120 novel precursor miRNAs, and 281 mature miRNAs. Some miRNAs showed stage- or tissue-enriched expression and miRNAs, such as the miR-17 ~ 92 cluster, myomiRs (miR-206), neuromiRs (miR-9, miR-124), miR-130b, and miR-430 showed differential expression in different temperature regimes. Long-term effect of embryonic incubation temperature was revealed on expression of some miRNAs in juvenile pituitary (miR-449), gonad (miR-27c, miR-30c, and miR-200a), and liver (let-7 h, miR-7a, miR-22, miR-34c, miR-132a, miR-192, miR-221, miR-451, miR-2188, and miR-7550), but not in brain. Some of differentially expressed miRNAs in the liver were confirmed using LNA-based rt-qPCR. The effect of temperature on methylation status of selected miRNA promoter regions was mostly inconclusive.

Conclusions

Temperature elevation by several degrees during embryonic and larval developmental stages significantly alters the miRNA profile, both short-term and long-term. Our results suggest that a further rise in seas temperature might affect life history of Atlantic cod.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1503-7) contains supplementary material, which is available to authorized users.

Polyphenols and the modulation of gene expression pathways: can we eat our way out of the danger of chronic disease?

Plant-derived dietary polyphenols may improve some disease states and promote health. Experimental evidence suggests that this is partially attributable to changes in gene expression. The rational use of bioactive food components may therefore present an opportunity to activate or repress selected gene expression pathways and, consequently, to manage or prevent disease. It remains to be determined whether this use of bioactive food components can be done safely. This article reviews the associated controversies and limitations of polyphenol therapy. There is a paucity of clinical data on the rational use of polyphenols, including a lack of knowledge on effective dosage, actual chemical formulations, bioavailability, distribution in tissues, the effect of genetic variations, differences in gut microflora, the synergistic (or antagonistic) effects observed in extracts, and the possible interaction between polyphenols and lipid domains of cell membranes that may alter the function of relevant receptors. The seminal question of why plants make substances that benefit humans remains unanswered, and there is still much to learn in terms of correlative versus causal effects of human exposure to various nutrients. The available data strongly suggest significant effects at the molecular level that represent interactions with the epigenome. The advent of relatively simple technologies is helping the field of epigenetics progress and facilitating the acquisition of multiple types of data that were previously difficult to obtain. In this review, we summarize the molecular basis of the epigenetic regulation of gene expression and the epigenetic changes associated with the consumption of polyphenols that illustrate how modifications in human nutrition may become relevant to health and disease.

MicroRNA in teleost fish

MicroRNAs (miRNAs) are transcriptional and posttranscriptional regulators involved in nearly all known biological processes in distant eukaryotic clades. Their discovery and functional characterization have broadened our understanding of biological regulatory mechanisms in animals and plants. They show both evolutionary conserved and unique features across Metazoa. Here, we present the current status of the knowledge about the role of miRNA in development, growth, and physiology of teleost fishes, in comparison to other vertebrates. Infraclass Teleostei is the most abundant group among vertebrate lineage. Fish are an important component of aquatic ecosystems and human life, being the prolific source of animal proteins worldwide and a vertebrate model for biomedical research. We review miRNA biogenesis, regulation, modifications, and mechanisms of action. Specific sections are devoted to the role of miRNA in teleost development, organogenesis, tissue differentiation, growth, regeneration, reproduction, endocrine system, and responses to environmental stimuli. Each section discusses gaps in the current knowledge and pinpoints the future directions of research on miRNA in teleosts.

Molecular responses to toxicological stressors: profiling microRNAs in wild Atlantic salmon (Salmo salar) exposed to acidic aluminum-rich water

Atlantic salmon (Salmo salar) is among the most sensitive organisms toward acidic, aluminum exposure. Main documented responses to this type of stress are a combination of hypoxia and loss of blood plasma ions. Physiological responses to stress in fish are often grouped into primary, secondary and tertiary responses, where the above mentioned effects are secondary responses, while primary responses include endocrine changes as measurable levels of catecholamines and corticosteroids. In this study we have exposed young (14 months) Atlantic salmon to acid/Al water (pH ≈ 5.6, Al(i) ≈ 80 μg L⁻¹) for 3 days, and obtained clear and consistent decrease of Na⁺ and Cl⁻ ions, and increases of glucose in blood plasma, hematocrit and P(CO₂) in blood. We did not measure plasma cortisol (primary response compound), but analyzed effects on microRNA level (miRNA) in muscle tissue, as this may represent initial markers of primary stress responses. miRNAs regulate diverse biological processes, many are evolutionarily conserved, and hundreds have been identified in various animals, although only in a few fish species. We used a novel high-throughput sequencing (RNA-Seq) method to identify miRNAs in Atlantic salmon and specific miRNAs as potential early markers for stress. A total of 18 miRNAs were significantly differentially expressed (FDR<0.1) in exposed compared to control fish, four down-regulated and 14 up-regulated. An unsupervised hierarchical clustering of significant miRNAs revealed two clusters representing exposed and non-exposed individuals. Utilizing the genome of the zebrafish and bioinformatic tools, we identified 224 unique miRNAs in the Atlantic salmon samples sequenced. Additional laboratory studies focusing on function, stress dose-responses and temporal expression of the identified miRNAs will facilitate their use as initial markers for stress responses.

Characterization of novel precursor miRNAs using next generation sequencing and prediction of miRNA targets in Atlantic halibut

BACKGROUND

microRNAs (miRNAs) are implicated in regulation of many cellular processes. miRNAs are processed to their mature functional form in a step-wise manner by multiple proteins and cofactors in the nucleus and cytoplasm. Many miRNAs are conserved across vertebrates. Mature miRNAs have recently been characterized in Atlantic halibut (Hippoglossus hippoglossus L.). The aim of this study was to identify and characterize precursor miRNA (pre-miRNAs) and miRNA targets in this non-model flatfish. Discovery of miRNA precursor forms and targets in non-model organisms is difficult because of limited source information available. Therefore, we have developed a methodology to overcome this limitation.

METHODS

Genomic DNA and small transcriptome of Atlantic halibut were sequenced using Roche 454 pyrosequencing and SOLiD next generation sequencing (NGS), respectively. Identified pre- miRNAs were further validated with reverse-transcription PCR. miRNA targets were identified using miRanda and RNAhybrid target prediction tools using sequences from public databases. Some of miRNA targets were also identified using RACE-PCR. miRNA binding sites were validated with luciferase assay using the RTS34st cell line.

RESULTS

We obtained more than 1.3 M and 92 M sequence reads from 454 genomic DNA sequencing and SOLiD small RNA sequencing, respectively. We identified 34 known and 9 novel pre-miRNAs. We predicted a number of miRNA target genes involved in various biological pathways. miR-24 binding to kisspeptin 1 receptor-2 (kiss1-r2) was confirmed using luciferase assay.

CONCLUSION

This study demonstrates that identification of conserved and novel pre-miRNAs in a non-model vertebrate lacking substantial genomic resources can be performed by combining different next generation sequencing technologies. Our results indicate a wide conservation of miRNA precursors and involvement of miRNA in multiple regulatory pathways, and provide resources for further research on miRNA in non-model animals.

In vivo screening of zebrafish microRNA responses to bacterial infection and their possible roles in regulating immune response genes after lipopolysaccharide stimulation

Micro (mi)RNAs are abundant small noncoding RNAs found in plants and animals, the regulatory functions of which are not fully understood in fish. To identify potential miRNAs, we screened an miRNA microarray with total RNA from zebrafish infected with Vibrio harveyi and another from uninfected zebrafish. Six miRNAs were obtained from the microarray screening. We studied miRNA expression patterns of 2 miRNAs (miR-122 and miR-194) after bacterial infection of transgenic zebrafish (containing tilapia hepcidin (TH)2-3) and non-transgenic zebrafish from which the 2 miRNAs were obtained from the microarray experiment. The results indicated that miR-122 and miR-194 were higher in PBS-injected zebrafish compared with TH2-3 zebrafish or wild-type (WT) zebrafish after V. harveyi infection. Overexpression of miRNAs (miR-122, miR-192, and miR-194a) was seen in zebrafish liver (ZFL) cells after lipopolysaccharide (LPS) treatment and in untreated fish. Our results showed that after 24 h of doxycycline treatment without LPS stimulation, interleukin (IL)-22, lysozyme, toll-like receptor (TLR)1, TLR3, TLR4a, and tumor necrosis factor (TNF)-α gene expressions were, respectively, upregulated by ~14-, 22-, 2.2-, 13-, 200-, and 38-fold in miR-122-transfected compared with non-transfected (WT) ZFL cells. In cells transfected with miR-192 and treated with LPS after 8-12 h, IL-22, lysozyme, TLR1, TLR3, TLR4a, and TNF-α expressions significantly differed between WT and miR-192-overexpressing ZFL cells. However, we observed significantly higher IL-22 expression levels after 12 h of LPS treatment in miR-192-transfected ZFL cells compared with non-transfected cells. In contrast, IL-22, lysozyme, and TNF-α were markedly upregulated (>100-fold) after miR-194a transfection and overexpression in ZFL cells and treatment with LPS. Our cloning and expression analyses indicated that miR-122, miR-192, and miR-194a play important roles in zebrafish immunology.

Identification of common carp (Cyprinus carpio) microRNAs and microRNA-related SNPs

BACKGROUND

MicroRNAs (miRNAs) exist pervasively across viruses, plants and animals and play important roles in the post-transcriptional regulation of genes. In the common carp, miRNA targets have not been investigated. In model species, single-nucleotide polymorphisms (SNPs) have been reported to impair or enhance miRNA regulation as well as to alter miRNA biogenesis. SNPs are often associated with diseases or traits. To date, no studies into the effects of SNPs on miRNA biogenesis and regulation in the common carp have been reported.

RESULTS

Using homology-based prediction combined with small RNA sequencing, we have identified 113 common carp mature miRNAs, including 92 conserved miRNAs and 21 common carp specific miRNAs. The conserved miRNAs had significantly higher expression levels than the specific miRNAs. The miRNAs were clustered into three phylogenetic groups. Totally 394 potential miRNA binding sites in 206 target mRNAs were predicted for 83 miRNAs. We identified 13 SNPs in the miRNA precursors. Among them, nine SNPs had the potential to either increase or decrease the energy of the predicted secondary structures of the precursors. Further, two SNPs in the 3' untranslated regions of target genes were predicted to either disturb or create miRNA-target interactions.

CONCLUSIONS

The common carp miRNAs and their target genes reported here will help further our understanding of the role of miRNAs in gene regulation. The analysis of the miRNA-related SNPs and their effects provided insights into the effects of SNPs on miRNA biogenesis and function. The resource data generated in this study will help advance the study of miRNA function and phenotype-associated miRNA identification.

Postprandial regulation of hepatic microRNAs predicted to target the insulin pathway in rainbow trout

Rainbow trout are carnivorous fish and poor metabolizers of carbohydrates, which established this species as a model organism to study the comparative physiology of insulin. Following the recent characterisation of key roles of several miRNAs in the insulin action on hepatic intermediary metabolism in mammalian models, we investigated the hypothesis that hepatic miRNA expression is postprandially regulated in the rainbow trout and temporally coordinated in the context of insulin-mediated regulation of metabolic gene expression in the liver. To address this hypothesis, we used a time-course experiment in which rainbow trout were fed a commercial diet after short-term fasting. We investigated hepatic miRNA expression, activation of the insulin pathway, and insulin regulated metabolic target genes at several time points. Several miRNAs which negatively regulate hepatic insulin signaling in mammalian model organisms were transiently increased 4 h after the meal, consistent with a potential role in acute postprandial negative feed-back regulation of the insulin pathway and attenuation of gluconeogenic gene expression. We equally observed a transient increase in omy- miRNA-33 and omy-miRNA-122b 4 h after feeding, whose homologues have potent lipogenic roles in the liver of mammalian model systems. A concurrent increase in the activity of the hepatic insulin signaling pathway and the expression of lipogenic genes (srebp1c, fas, acly) was equally observed, while lipolytic gene expression (cpt1a and cpt1b) decreased significantly 4 h after the meal. This suggests lipogenic roles of omy-miRNA-33 and omy-miRNA-122b may be conserved between rainbow trout and mammals and that these miRNAs may furthermore contribute to acute postprandial regulation of de novo hepatic lipid synthesis in rainbow trout. These findings provide a framework for future research of miRNA regulation of hepatic metabolism in trout and will help to further elucidate the metabolic phenotype of rainbow trout.

Differential expression patterns of conserved miRNAs and isomiRs during Atlantic halibut development

BACKGROUND

MicroRNAs (miRNAs) play a major role in animal ontogenesis. Size variants of miRNAs, isomiRs, are observed along with the main miRNA types, but their origin and possible biological role are uncovered yet. Developmental profiles of miRNAs have been reported in few fish species only and, to our knowledge, differential expressions of isomiRs have not yet been shown during fish development. Atlantic halibut, Hippoglossus hippoglossus L., undergoes dramatic metamorphosis during early development from symmetrical pelagic larval stage to unsymmetrical flatfish. No data exist on role of miRNAs in halibut metamorphosis.

RESULTS

miRNA profiling using SOLiD deep sequencing technology revealed a total of 199 conserved, one novel antisense, and one miRNA* mature form. Digital expression profiles of selected miRNAs were validated using reverse transcription quantitative PCR. We found developmental transition-specific miRNA expression. Expression of some miRNA* exceeded the guide strand miRNA. We revealed that nucleotide truncations and/or additions at the 3' end of mature miRNAs resulted in size variants showing differential expression patterns during the development in a number of miRNA families. We confirmed the presence of isomiRs by cloning and Sanger sequencing. Also, we found inverse relationship between expression levels of sense/antisense miRNAs during halibut development.

CONCLUSION

Developmental transitions during early development of Atlantic halibut are associated with expression of certain miRNA types. IsomiRs are abundant and often show differential expression during the development.

Analysis of the Asian seabass transcriptome based on expressed sequence tags

Analysis of transcriptomes is of great importance in genomic studies. Asian seabass is an important fish species. A number of genomic tools in it were developed, while large expressed sequence tag (EST) data are lacking. We sequenced ESTs from nine normalized cDNA libraries and obtained 11 431 high-quality ESTs. We retrieved 8524 ESTs from dbEST database and analyzed all 19 975 ESTs using bioinformatics tools. After clustering, we obtained 8837 unique sequences (2838 contigs and 5999 singletons). The average contig length was 574 bp. Annotation of these unique sequences revealed that 48.9% of them showed significant homology to RNA sequences in GenBank. Functional classification of the unique ESTs identified a broad range of genes involved in different functions. We identified 6114 putative single-nucleotide polymorphisms and 634 microsatellites in ESTs. We discovered different temporal and spatial expression patterns of some immune-related genes in the Asian seabass after challenging with a pathogen Vibrio harveyi. The unique EST sequences are being used in developing a cDNA microarray to examine global gene expression and will also facilitate future whole-genome sequence assembly and annotation of Asian seabass and comparative genomics.

A high-resolution linkage map for comparative genome analysis and QTL fine mapping in Asian seabass, Lates calcarifer

Background

High density linkage maps are essential for comparative analysis of synteny, fine mapping of quantitative trait loci (QTL), searching for candidate genes and facilitating genome sequence assembly. However, in most foodfish species, marker density is still low. We previously reported a first generation linkage map with 240 DNA markers and its application to preliminarily map QTL for growth traits in Asian seabass (Lates calcarifer). Here, we report a high-resolution linkage map with 790 microsatellites and SNPs, comparative analysis of synteny, fine-mapping of QTL and the identification of potential candidate genes for growth traits.

Results

A second generation linkage map of Asian seabass was developed with 790 microsatellite and SNP markers. The map spanned a genetic length of 2411.5 cM, with an average intermarker distance of 3.4 cM or 1.1 Mb. This high density map allowed for comparison of the map with Tetraodon nigroviridis genome, which revealed 16 synteny regions between the two species. Moreover, by employing this map we refined QTL to regions of 1.4 and 0.2 cM (or 400 and 50 kb) in linkage groups 2 and 3 in a population containing 380 progeny; potential candidate genes for growth traits in QTL regions were further identified using comparative genome analysis, whose effects on growth traits were investigated. Interestingly, a QTL cluster at Lca371 underlying growth traits of Asian seabass showed similarity to the cathepsin D gene of human, which is related to cancer and Alzheimer's disease.

Conclusions

We constructed a high resolution linkage map, carried out comparative mapping, refined the positions of QTL, identified candidate genes for growth traits and analyzed their effects on growth. Our study developed a framework that will be indispensable for further identification of genes and analysis of molecular variation within the refined QTL to enhance understanding of the molecular basis of growth and speed up genetic improvement of growth performance, and it also provides critical resource for future genome sequence assembly and comparative genomics studies on the evolution of fish genomes.