Integrative microRNAome analysis of skeletal muscle of Colossoma macropomum (tambaqui), Piaractus mesopotamicus (pacu), and the hybrid tambacu, based on next-generation sequencing data
BMC Genomics volume 22, Article number: 237 (2021)
Colossoma macropomum (tambaqui) and Piaractus mesopotamicus (pacu) are good fish species for aquaculture. The tambacu, individuals originating from the induced hybridization of the female tambaqui with the male pacu, present rapid growth and robustness, characteristics which have made the tambacu a good choice for Brazilian fish farms. Here, we used small RNA sequencing to examine global miRNA expression in the genotypes pacu (PC), tambaqui (TQ), and hybrid tambacu (TC), (Juveniles, n = 5 per genotype), to better understand the relationship between tambacu and its parental species, and also to clarify the mechanisms involved in tambacu muscle growth and maintenance based on miRNAs expression.
Regarding differentially expressed (DE) miRNAs between the three genotypes, we observed 8 upregulated and 7 downregulated miRNAs considering TC vs. PC; 14 miRNAs were upregulated and 10 were downregulated considering TC vs. TQ, and 15 miRNAs upregulated and 9 were downregulated considering PC vs. TQ. The majority of the miRNAs showed specific regulation for each genotype pair, and no miRNA were shared between the 3 genotype pairs, in both up- and down-regulated miRNAs. Considering only the miRNAs with validated target genes, we observed the miRNAs miR-144-3p, miR-138-5p, miR-206-3p, and miR-499-5p. GO enrichment analysis showed that the main target genes for these miRNAs were grouped in pathways related to oxygen homeostasis, blood vessel modulation, and oxidative metabolism.
Our global miRNA analysis provided interesting DE miRNAs in the skeletal muscle of pacu, tambaqui, and the hybrid tambacu. In addition, in the hybrid tambacu, we identified some miRNAs controlling important molecular muscle markers that could be relevant for the farming maximization.
Colossoma macropomum (tambaqui) and Piaractus mesopotamicus (pacu) are fish species that have widely accepted in the consumer market. In addition, they present desirable characteristics for an intensive breeding environment as rapid growth and optimal adaptation to artificial feeding .
Noncoding RNAs have become a very important tool for carrying out different types of evolutionary and expression profile experiments. In such way, many databases and protocols for annotation have been developed [2,3,4].
It is known that Hybrid individuals possess desirable characteristics for production, such as high growth rate, higher resistance to disease, and higher quality of meat . For this reason, several producers have chosen to cultivate the hybrid between the induced crossing of the tambaqui female and the pacu male, the tambacu. Tambacu, in spite of the few genetic information, has the capacity to be more resistant to parasites and stress . Furthermore, tambacu is considered a fish with great potential for Brazilian aquaculture, since it presents a high and fast growth rate, including the skeletal muscle and higher resistance to low temperatures, which contribute to increasing their rusticity , thus representing an important model to study. The skeletal muscle is directly involved in the growth of the fish, corresponding to about 35–60% of the body weight of the animal . This abundant muscle mass enables the survival of the animals in the aquatic environment and is commercially important for animal production, since it is one of the most important food sources for human diet .
With the focus of massive data analysis, non-coding RNAs have been highly explored in recent years considering their strong impact on controlling several biological processes [10, 11]. Among non-coding RNAs, miRNAs (miRNAs), a well reported class of small non-coding RNAs, are known to perform a fine regulation of gene activity in a post-transcriptional mode upon their association with RNA-induced silencing complex (RISC) and binding via base-complementarity to target mRNAs, resulting in their translational repression and/or in their degradation [12, 13].
Some miRNAs are known to be specifically or highly expressed in the cardiac and/or skeletal muscles and have been dubbed “myo-miRs” . In fish, recent studies have brought information regarding the evolution and genomic organization of myo-miRNAs .
Considering the important role of miRNAs in cell physiology, the aim of our study was to evaluate the global miRNA profile in skeletal muscle of pacu, tambaqui, and the hybrid tambacu, to investigate differentially expressed miRNAs and the ontology analysis considering the comparison between the hybrid with and its parental genotypes, to better understand the molecular signaling pathways involved in tambacu muscle development, growth, and maintenance. In this way, the present work can be an important source of information to supports studies that address the advantages of adopting hybrids for cultivation, which are still scarce .
Early juvenile fish Colossoma macropomum, Piaractus mesopotamicus, and hybrid tambacu were studied. All the extracted RNA samples were analyzed by NanoVue (GE Healthcare), and Bioanalyzer (Agilent Technologies) and only samples with a RIN ≥ 8 were selected for the sequencing procedures.
After preprocessing steps, the number of reads decreased to approximately 50% due to the removal of low-quality reads and the adapters, and the removal of sequences not matching the minimum/maximum size corresponding to miRNAs (Fig. 1, Table 1, and Additional file 1). We obtained a total alignment average of 51,78%, considering all samples analyzed. As a result of featureCounts processing (see Additional files 2 and 3), a count matrix was generated for each miRNA present in the reference (see Additional file 4). As part of the pipeline, Principal Component Analysis (PCA) (Fig. 2a) and Dispersion Analysis (Fig. 2b) were run to cluster the samples based on the expression values and to observe the quality of the data. We observed the samples separated into clusters among samples, indicating that the different groups really correspond to different genotypes. Also, the Dispersion Analysis showed that below a mean read count of 10, the dispersion of the data increased dramatically. Thus, miRNAs presenting a mean read count < 10 were filtered out to keep only miRNAs not presenting high dispersion levels.
Differentially expressed (DE) miRNAs groups were observed in a row-clustered heatmap containing all the comparison genotype pairs (TC vs. PC, TC vs. TQ, and PC vs. TQ) and all the miRNAs detected as DE (Fig. 3). Based on color histogram, it was possible to observe the variations in terms of expression levels between up-regulated and down-regulated miRNAs throughout the three genotypes analyzed. A relative high number of DE miRNAs are observed in the data. It Is possible to detect such a balance of DE miRNAs, where there is no comparison only with up- or down-regulated miRNAs (Figs. 4 and 5).
DE analysis showed 8 upregulated and 7 downregulated miRNAs for TC vs. PC genotype pair (Fig. 4a and d, Table 2 and Additional file 5), 14 upregulated and 10 downregulated miRNAs for TC vs. TQ genotype pair (Fig. 4b and e, Table 3 and Additional file 6), and 15 upregulated and 9 downregulated miRNAs for the PC vs. TQ genotype pair (Fig. 4c and f, Table 4 and Additional file 7). A full list exclusive for DE miRNAs can be observed on Additional file 8.
According to the Venn diagram data (Figs. 6 and 7), we analyzed the miRNA distribution pattern between the parental species and the hybrid, and it was possible to observe some DE miRNA-sharing between the different genotype’s pairs. For example, considering the upregulated set of miRNAs, miR-144-3p and miR-10a-5p appeared to be upregulated in TC compared to both PC and TQ (Figs. 6 and 8). This leads us to think that the upregulation of these miRNAs could be a result of a multifactorial characteristic involving both profiles given that such miRNAs are upregulated only when the contrast involves the hybrid and such miRNAs are not DE between parental species. A similar characteristic was observed in the down-regulated miRNAs, where miR-199-3-3p appears downregulated in TC when compared with both PC and TQ. These miRNAs are considered as TC-exclusive (in terms of expression levels), as they appear to be up- or down-regulated only in the hybrid compared to the parental species, but not between parental species (Fig. 8). Also, we detected some up- or down-regulated miRNAs only in the hybrid, representing a specific characteristic of this genotype.
We then assigned different miRNAs in 6 different categories regarding inheritance characteristics related to expression levels, i.e., (i) TC-exclusive (for those miRNAs differentially expressed only in TC compared with both parental species), (ii) PC-inherited (for those miRNAs with expression patterns similar to PC expression levels), (iii) TQ-inherited (for those miRNAs with expression patterns similar to TQ expression levels), (iv) TC vs. PC-exclusive (for miRNAs that are differentially expressed only when TC and PC are involved in comparisons), (v) PC vs. TQ-exclusive (for miRNAs that are differentially expressed only when PC and TQ are involved in comparisons), and (vi) TC vs. TQ-exclusive (for miRNAs that are differentially expressed only when TC and TQ are involved in comparisons). miRNAs are also discriminated in up- (▲) and downregulated (▼) (Fig. 8).
In the search for target genes involved in DE miRNAs expression, we searched for experimentally validated targets using miRTarBase dataset version 7 . Two networks comprising only validated data were detected (Figs. 9 and 10). A network presenting strong validated interactions between target genes and DE miRNAs expression for Homo sapiens, Mus musculus, Rattus norvegicus, and Danio rerio are presented (Fig. 9). We observed interactions among 31 known miRNAs, and the top five miRNAs with a higher number of interactions were hsa-miR-221-3p, has-miR-27b-3p, hsa-miR-138-5p, mmu-miR-206-3p, and hsa-miR-132-3p, targeting respectively 72, 52, 47, 32, and 32 target genes. Considering the Danio rerio validated interaction data (Fig. 10), we observed that the network presents interactions involving four DE miRNAs and 10 target genes: dre-miR-138-5p (vcana), dre-miR-206-3p (vegfa and jun), dre-miR-499-5p (cyb561d2), and dre-miR-144-3p (lmo2, klfd, gata2a, meis1, klf3, and alas2).
Ontology analysis carried out by using Enrichr , considering the genes validated as targets of the DE miRNAs (considering the Danio rerio validated interactions) showed that the interacting-validated genes are involved in various categories of ontology. For Biological Processes, Cellular Components, and Molecular Functions, the top enriched terms were, respectively, oxygen homeostasis (the most representative gene was 5-aminolevulinate synthase 2 - alas2), nuclear euchromatin (the most representative gene was jun), and transcriptional repressor activity, RNA polymerase II activating transcription factor binding (the most representative genes were jun and lmo2) (Figs. 11, 12 and 13 and Table 5).
Next-generation sequencing has been widely applied for global analysis in several studies involving miRNAs and skeletal muscle. Nachtigall et al. employed next-generation sequencing in Nile Tilapia (Oreochromis niloticus), to obtain information regarding evolutionary pathways with emphasis on muscle miRNAs . The authors identified that there are large syntenic blocks in the genome, possibly being linked to a common function of such miRNAs, and as well specific role of miR-499 expression . Gomes et al. also used next generation sequencing to analyze samples of liver and skin from tambaqui (Colossoma macropomum) to characterize and identify the expression levels of miRNAs and the interaction with target genes having D. rerio as reference . The analysis performed showed that although there are tissue-specific miRNAs expression profile, some miRNAs can be shared between different tissues, as liver and skin. The most expressed miRNAs in both tissues enriched signaling pathways that control several biological processes with a large gene network.
Previous studies in our group have also analyzed the expression of some muscle specific miRNAs in the skeletal muscle of fish involving fasting and re-feeding treatments and muscle development [20, 21]. The authors observed that some miRNAs (miR-1, miR-133, miR-155, miR-206, and miR-499) presented a possible role in the regulation of factors related to muscle cell proliferation and differentiation and with muscle performance and metabolism, modulating the rate of protein synthesis and degradation.
In the present study, global analysis of miRNAs in the three genotypes identified a small number of differentially expressed miRNAs considering the comparison between the genotypes pairs. The comparison between tambacu and P. mesopotamicus (TC vs. PC) showed 15 DE miRNAs, 8 upregulated and 7 downregulated. Considering the comparison between tambacu and C.macropomum (TC vs. TQ), we observed 24 DE miRNAs, 14 upregulated and 10 downregulated, and considering the comparison between P. mesopotamicus and C. macropomum (PC vs. TQ), we observed 24 DE miRNAs, 15 upregulated and 9 downregulated. Interestingly, the majority of the miRNAs were specific for each genotype pair. No miRNA was shared between the 3 genotype pairs, both up and downregulated. Between the upregulated DE miRNAs, 9 overlapped between PC vs. TQ and TC vs. TQ and 2 overlap between TC vs. PC and TC vs. TQ. Considering the downregulated DE miRNAs, 3 overlap between PC vs. TQ and TC vs. TQ, and only one overlaps between TC vs. PC and TC vs. TQ. No miRNA was differentially expressed between all of the three comparison pairs at the same time. Some miRNAs show characteristics that make it possible to track their inheritance pattern by looking into the different groups. For example, some miRNAs appear to be exclusive for the comparison between TC vs. TQ (miR-187, miR-18a, miR-190a, miR-206-3p, miR-221-3p, miR-27b-3p, miR-29a, miR-363-3p, and miR-727-5p), some of which were exclusive in the comparison between PC vs. TQ (miR-146b, miR-212-5p, miR-2188-3p, miR-26a-2-3p, and miR-489). Interestingly, some miRNAs presented a characteristic of being inherited specifically of one genotype, and also, there were some miRNAs that appeared to be up- or down-regulated only in the hybrid (miR-10a-5p, miR-144-3p, and miR-199-3-3p). These miRNAs presented (up or downregulated) behavior, exclusively in the hybrid, not differentially expressed between the parental genotypes.
Ontology analysis by Enrichr, considering the Danio rerio validated interaction data, showed that the experimentally validated target genes of the DE miRNAs were grouped into several Biological Processes, Cellular Components, and Molecular Functions categories. The most enriched processes were oxygen homeostasis, nuclear euchromatin and transcriptional repressor activity, and RNA polymerase II activating transcription factor binding, respectively. The more enriched genes in the signaling pathway involved in these processes were alas2, jun, and lmo2.
miR-206-3p, identified in the present study, was differentially expressed in the comparison pair TC vs. TQ, targets jun and vegfa genes. Jun is a member of the Activator Protein 1 (AP-1) transcription factor family that regulates cell proliferation and differentiation, apoptosis, cellular migration, inflammation, and cell-cell interaction . On the other hand, vegfa is involved in vascular development and new blood vessel formation [23, 24] and stimulates endothelial cell migration by activating AP-1 transcription factor jun . Skeletal muscle is the most abundant source of VEGFA [26, 27], and skeletal-muscle-specific VEGFA knockout changed angiogenesis in muscle fibers .
As jun and vegfa genes are targets of miR-206-3p and, as this miRNA was observed to be differentially expressed in the comparison pair TC vs. TQ, this is an indication that this miRNA can represent a remarkable characteristic for the hybrid, since these genes are involved in the maintenance of blood irrigation that, in turn, controls the oxygen rates in the tissues.
The miRNA involved in the control of the alas2 and lmo2 genes was miR-144-3p. This miRNA appeared upregulated in the comparison pairs TC vs. PC and TC vs. TQ and did not appear in the comparison between the parental genotypes PC vs. TQ.
ALAS2 is one of the isozymes of ALAS (5-aminolevulinate synthase) involved in vertebrate heme biosynthesis. This gene is expressed preferentially in erythroid cell-specific mitochondrial enzymes [29, 30] and catalyzes the biosynthesis of bulk heme for hemoglobin production [31,32,33,34]. ALAS2 is also regulated by hypoxia-inducible factor HIF. Khenchaduri and colaborators observed the overexpression of ALAS2 in cardiac myoblasts submitted to chronic hypoxia, with a corresponding increase in cellular heme levels. They concluded that, similar to erythroid cells, ALAS2 is positively regulated by hypoxia in cardiac myoblasts with an increase in heme levels .
According to Zhang et al., the upregulation of alas2 during hypoxia is directly mediated by a transcription factor, hypoxia-inducible factor 1 (Hif1) , and an increased Hif-1 that occurs under low oxygen concentration conditions also promotes the activation of genes involved in hypoxia, such as vegf .
Moreover, it has been shown that the transcription factor LIM domain only 2 (Lmo2) is required for angiogenesis and is important for angiogenic remodeling of the existing capillary network into mature blood vessels [38,39,40]. Besides, Lmo2 forms common transcription complexes with GATA-binding protein 2 (GATA2) and cooperatively regulate VEGF-induced angiogenesis .
Therefore, it is possible that the interaction between miR-144-3p and alas2 and lmo2 genes can regulate important pathways related to muscle vascularization, oxygen homeostasis, and oxidative metabolism. We believe that this is an important characteristic that probably could provide, for the hybrid tambacu, more resistance to environments with lower oxygen content, once this miRNA was upregulated only in the hybrid, in the comparison pairs TC vs. PC and TC vs. TQ.
In our study, miR-144-3p was identified as DE in TC vs. PC and TC vs. TQ, and the miR-206-3p described as DE only in the comparison pair TC vs. TQ demonstrated that these miRNAs may represent an important species-specific feature. The hybrid tambacu is a very popular fish in Brazil, especially in the south and southeast regions and combines important characteristics for commercial fish farms, such as resistance to low temperatures during the tropical winter, similar to pacu, and faster growth, similar to tambaqui [1, 6]. Also, juvenile tambacu exhibits a total capacity for compensatory growth when submitted to short periods of food deprivation, followed by refeeding, a fact that could contribute to reducing production costs .
Although there are few studies showing the real potential of the hybrid for aquaculture production [42, 43], some of them have been done to better understand the genetics and molecular characteristics of this genotype. Recently, Gomes et al. used Next Generation Sequencing technologies for the analysis of muscle and skin transcriptome of C. macropomum and the hybrid tambacu, observed differences in the specificity and gene expression levels in both the muscle and skin of the hybrid tambacu and the parental tambaqui . The authors demonstrated that in the muscle of tambaqui, up-regulated genes involved in muscle contraction, such as the myosin family, actin and myomesin, and catalytic protein genes. The hybrid tambacu muscle presented upregulated genes controlling oxidative stress, amino acid metabolism, the ubiquitin family, defense, and heat shock protein (HSP) family genes, involved in stress.
In our study, we can emphasize that the interaction of jun and vegfa genes with miR-206-3p (TC vs. TQ) and lmo2 and alas2 genes with miR-144-3p (TC vs. PC and TC vs. TQ) could provide an evolutionary advantage to the hybrid tambacu once these genes are controlling pathways related to blood vessel modulation, oxygen homeostasis, and oxidative metabolism. Taking this into account, the differential regulation of these miRNAs may represent an explanation of why the hybrid has a better adaptation from exposure to hypoxic environments [38, 45]. Besides, it could indicate an important function of miRNAs miR-206-3p and miR-144-3p in the control of the robustness and faster growth of the hybrid tambacu .
In addition, the most significant molecular function enriched by the targets of the up-regulated miRNAs in the tambacu muscle was transcriptional repressor activity. The inhibition of this function reinforces our hypothesis that these miRNAs could be related to the robustness and faster growth of the hybrid, since their activity may promote enhanced transcriptional activity, which combined with translation, is an essential process for muscle hypertrophy and growth .
Using comparative global miRNA analysis in the hybrid tambacu and their parental genotypes Piaractus mesopotamicus and Colossoma macropomum, we demonstrated, for the first time, miRNAs differentially expressed in the skeletal muscle of the genotype comparison pairs (TC vs. PC, TC vs. TQ, and PC vs. TQ). Furthermore, we showed that the hybrid expressed miRNAs that control important molecular muscle markers with potential applications in aquaculture management systems. Our study also provides insight for further investigations involving the validation of these miRNAs and their signaling pathway components in muscle phenotype in these genotypes.
Experimental design, sample collection, Total RNA extraction, and small RNAs sequencing
Pure early juvenile fish Colossoma macropomum, Piaractus mesopotamicus, and hybrid tambacu (Juveniles, n = 5 per genotype) were obtained at the Aquaculture Center of UNESP - CAUNESP, and were maintained in aerated tanks in the Agribusiness Technology Agency (APTA), Presidente Prudente, Sao Paulo, until fish attained ~ 10 g. Fish were kept in a natural photoperiod (12 h light/dark), and water temperature around 28 °C. Dissolved oxygen and pH were monitored daily, and ammonia, nitrite, and nitrate were monitored weekly. Fish were fed with Guabi-Pirá 28 fish feed. Excess benzocaine (concentration over 250 mg/L) was applied to euthanize the animals. Afterward, body weight and total length were determined, and fast muscle samples were collected from the epaxial region near the dorsal fin. Muscle samples were then frozen in liquid nitrogen for molecular analysis.
Total RNA was isolated using TRIZOL (Invitrogen), following the manufacturer’s specifications. The total extracted RNA was treated with DNase (Life Technologies, USA), following the manufacturer’s specifications. All extracted RNA samples were analyzed by NanoVue Spectrophotometer (GE Healthcare) and Bioanalyzer (Agilent Technologies) for checking the total concentration and the RNA Integrity Number (RIN). The RNA extraction resulted in high quality RNA, and the RIN obtained was ≥8. Samples were submitted to the sequencing, comprehending 15 sequenced samples in total. The samples were shipped to LC Sciences (Houston, TX, USA) in dry ice, and the entire process of library preparation and sequencing by Illumina platform was carried out. Libraries preparation were carried out using the TruSeq Small RNA Library Preparation Kit, and samples were run in 50 cycle single-end Illumina HiSeq 2500 Fast Mode. After sequencing, the data were downloaded through a link provided by the company.
The processing pipeline consisted initially of analyzing the data regarding the read quality (Phred Score). For this purpose, the software FastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc) version 0.10.1 was used. After obtaining the sequencing quality information, adapters were removed using the software Fastx-toolkit (http://hannonlab.cshl.edu/fastx_toolkit/index.html). Reads not containing adapters were kept in libraries. Reads that, after removing the adapters, presented less than 17 nucleotides in length were permanently excluded. A filter based on read-quality parameters using the Fastx-toolkit software was also carried out. In such way, 90% of the read content has a phred score of at least 30. Reads not satisfying these quality prerequisites were also permanently excluded. The libraries were then filtered by max length using a custom Perl script (Additional file 9), where reads longer than 26 nucleotides were discarded. The remaining reads were used to obtain an expression profile of known miRNAs, based on known miRNAs lists of Danio rerio, using data from mature and precursor miRNAs obtained from the miRBase database version 22 . For return more alignment than other species, we selected Danio rerio as a reference for this study, which is considered a good biological model.
For the differential expression analysis of miRNAs, mature sequences were aligned against the precursor sequences using GMAP , using standard parameters, to create an annotation file in GFF format, making it possible to locate the mature sequences in the precursor sequences and differentiate alignments between the 5′ and 3′ miRNAs.
After building the annotation file, the filtered reads were aligned against the precursor references using Hisat2  with standard parameters. The software featureCounts  was applied over the alignments to collect the read-count for each miRNA aligned against the precursor reference and matching the GFF annotation file. The counts were then analyzed to detect differentially expressed (DE) miRNAs using the Bioconductor/R package DESeq2 .
The three different genotypes received specific names regarding its popular names: PC - Piaractus mesopotamicus (pacu), TQ - Colossoma macropomum (tambaqui), and TC (hybrid tambacu).
A priori, we performed contrasts between the hybrid and each parental species, always having the hybrid as a base of comparison (TC vs. PC and TC vs. TQ). A posteriori, for further understanding of the traits inherited by each parental, in terms of miRNA expression profiles, we also have run a contrast analysis between the two parental genotypes. Such a second contrast has pacu as a base of comparison (PC vs. TQ). In this way, three contrasts were obtained: TC vs. PC, TC vs. TQ, and PC vs. TQ (see Additional files 5, 6 and 7).
For considering DE characteristics, we established that miRNAs presenting log2FoldChange (log2FC) ≤ − 0.75 and ≥ 0.75, followed by adjusted p-value (padj) ≤ 0.05, were considered differentially expressed as downregulated and upregulated, respectively.
The list of DE miRNAs detected by DESeq2 analysis was applied to the drawing of Venn Diagrams to visualize the DE miRNA sharing characteristics of each genotype by using Venn 2.1.0 (http://bioinfogp.cnb.csic.es/tools/venny/). Venn diagrams were created in two ways, i.e., up-regulated miRNAs and down-regulated miRNAs.
The R package gplots was applied with the DE miRNAs list for obtain heatmaps, which gives us a graphical idea of the expression variations between all the contrasts.
Target genes were detected using MirTarBase  database and the results were applied to ontology analysis by using Enrichr  to retrieve information regarding pathways in which the target genes are involved. The MirTarBase search comprises only experimentally validated interactions classified as “strong evidence”.
Interaction networks were assembled using Cytoscape 3.7.1  based on MirTarBase recovered data for a better observation of the interaction between DE miRNAs and their validated targets. Two sets of networks have been assembled: (i) a network containing validated interactions for Homo sapiens, Mus musculus, Rattus norvegicus, and Danio rerio, and (ii) a network containing only Danio rerio validated interactions.
All bioinformatics steps can be observed in a flowchart presented (Additional file 10).
Availability of data and materials
All raw and processed data have been deposited on the GEO dataset, under access number GSE147532.
- DE miRNAs:
Differentially expressed miRNAs
Principal Component Analysis
- vcana :
- vegfa :
Vascular Endotelial Growth Factor A
- jun :
AP-1 Transcription factor Subunit
- cyb561d2 :
Cytochrome B561 Family Member D2
- lmo2 :
LIM Domain Only 2
- klfd :
Kruppel-like factor d
- gata2a :
GATA-binding protein 2
- meis1 :
Meis Homeobox 1
- klf3 :
Kruppel Like Factor 3
- alas2 :
5-aminolevulinate synthase 2
- hif1 :
Hypoxia-Inducible Factor 1
Brazilian College of Animal Experimentation
Agribusiness Technology Agency
RNA Integrity Number
Araújo-Lima C, Goulding M. Os frutos do tambaqui: ecologia, conservação e cultivo na Amazônia; 1998.
Paschoal AR, Maracaja-Coutinho V, Setubal JC, Simões ZLP, Verjovski-Almeida S, Durham AM. Non-coding transcription characterization and annotation. RNA Biol. 2012;9(3):274–82. https://doi.org/10.4161/rna.19352.
Maracaja-Coutinho V, Paschoal AR, Caris-Maldonado JC, Borges PV, Ferreira AJ, Durham AM. Noncoding RNAs Databases: Current Status and Trends. Methods Mol Biol. 2019;1912:251–85. https://doi.org/10.1007/978-1-4939-8982-9_10.
Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ. miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res. 2006;34(Database issue):140–4.
Bartley DM, Rana K, Immink AJ. The use of inter-specific hybrids in aquaculture and fisheries. Rev Fish Biol Fish. 2000;10(3):325–37. https://doi.org/10.1023/A:1016691725361.
Tavares-dias M, De Moraes FR, Makoto E, Rezende PCB. Changes in blood parameters of hybrid tambacu fish parasitized by Dolops carvalhoi ( Crustacea , Branchiura ), a fish louse. Vet Arh. 2007;77(4):355–63.
A. L. Alves, E. S. Varela, G. V. Moro, and L. N. Kirschnik, “Riscos genéticos da produção de híbridos de peixes nativos,” Embrapa Pesca e Aquicultura, p. 60, 2014.
Listrat A, et al. How muscle structure and composition determine meat quality. Prod Anim. 2015;28(2):125–36.
Pal J, Shukla BN, Maurya AK, Verma HO. A review on role of fish in human nutrition with special emphasis to essential fatty acid. Int J Fish Acquat Stud. 2018;6(2):427–30.
Morris KV, Mattick JS. The rise of regulatory RNA. Nat Rev Genet. 2014;15(6):423–37. https://doi.org/10.1038/nrg3722.
Rinn J, Guttman M. RNA and dynamic nuclear organization. Science (80- ). 2014;345(6202):1240–1. https://doi.org/10.1126/science.1252966.
Bartel DP. MicroRNAs. Cell. 2004;116(2):281–97. https://doi.org/10.1016/S0092-8674(04)00045-5.
Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136(2):215–33. https://doi.org/10.1016/j.cell.2009.01.002.
McCarthy JJ. The MyomiR network in skeletal muscle plasticity. Exerc Sport Sci Rev. 2011;39(3):150–4. https://doi.org/10.1097/JES.0b013e31821c01e1.
Nachtigall PG, Dias MC, Pinhal D. Evolution and genomic organization of muscle microRNAs in fish genomes. BMC Evol Biol. 2014;14(1):1–12.
Chou C-H, Shrestha S, Yang CD, Chang NW, Lin YL, Liao KW, Huang WC, Sun TH, Tu SJ, Lee WH, Chiew MY, Tai CS, Wei TY, Tsai TR, Huang HT, Wang CY, Wu HY, Ho SY, Chen PR, Chuang CH, Hsieh PJ, Wu YS, Chen WL, Li MJ, Wu YC, Huang XY, Ng FL, Buddhakosai W, Huang PC, Lan KC, Huang CY, Weng SL, Cheng YN, Liang C, Hsu WL, Huang HD. miRTarBase update 2018: a resource for experimentally validated microRNA-target interactions. Nucleic Acids Res. 2018;46(D1):D296–302. https://doi.org/10.1093/nar/gkx1067.
Chen EY, et al. Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinformatics. 2013;14(1):128.
Nachtigall PG, Dias MC, Carvalho RF, Martins C, Pinhal D. MicroRNA-499 Expression Distinctively Correlates to Target Genes sox6 and rod1 Profiles to Resolve the Skeletal Muscle Phenotype in Nile Tilapia. PLoS One. 2015;10(3):e0119804.
Gomes F, Watanabe L, Nozawa S, Oliveira L, Cardoso J, Vianez J, Nunes M, Schneider H, Sampaio I. Identification and characterization of the expression profile of the microRNAs in the Amazon species Colossoma macropomum by next generation sequencing. Genomics. 2017;109(2):67–74. https://doi.org/10.1016/j.ygeno.2017.02.001.
Duran BODS, et al. Differential microRNA Expression in Fast- and Slow-Twitch Skeletal Muscle of Piaractus mesopotamicus during Growth. PLoS One. 2015;10(11):e0141967.
de Paula TG, et al. Food restriction increase the expression of mTORC1 complex genes in the skeletal muscle of juvenile pacu (Piaractus mesopotamicus). PLoS One. 2017;12(5):e0177679.
Shaulian E. AP-1 - the Jun proteins: oncogenes or tumor suppressors in disguise? Cell Signal. 2010;22(6):894–9. https://doi.org/10.1016/j.cellsig.2009.12.008.
Dvorak HF, Detmar M, Claffey KP, Nagy JA, Van De Water L, Senger DR. Vascular permeability factor/ vascular endothelial growth factor: an important mediator of angiogenesis in malignancy and inflammation. Int Arch Allergy Immunol. 1995;107(1–3):233–5. https://doi.org/10.1159/000236988.
Choorapoikayil S, Weijts B, Kers R, de Bruin A, den Hertog J. Loss of Pten promotes angiogenesis and enhanced vegfaa expression in zebrafish. Dis Model Mech. 2013;6(5):1159–66. https://doi.org/10.1242/dmm.012377.
Jia J, Ye T, Cui P, Hua Q, Zeng H, Zhao D. AP-1 transcription factor mediates VEGF-induced endothelial cell migration and proliferation. Microvasc Res. 2016;105:103–8. https://doi.org/10.1016/j.mvr.2016.02.004.
Logsdon EA, Finley SD, Popel AS, MacGabhann F. A systems biology view of blood vessel growth and remodelling. J Cell Mol Med. 2014;18(8):1491–508. https://doi.org/10.1111/jcmm.12164.
Yen P, Finley SD, Engel-Stefanini MO, Popel AS. A two-compartment model of VEGF distribution in the mouse. PLoS One. 2011;6(11).
Tang K, Breen EC, Gerber HP, Ferrara NMA, Wagner PD. Capillary regression in vascular endothelial growth factor-deficient skeletal muscle. Physiol Genomics. 2004;18(1):63–9. https://doi.org/10.1152/physiolgenomics.00023.2004.
Yamamoto M, Yew NS, Federspiel M, Dodgson JB, Hayashi N, Engel JD. Isolation of recombinant cDNAs encoding chicken erythroid delta-aminolevulinate synthase. Proc Natl Acad Sci. 1985;82(11):3702–6. https://doi.org/10.1073/pnas.82.11.3702.
Riddle RD, Yamamoto M, Engel JD. Expression of delta-aminolevulinate synthase in avian cells: separate genes encode erythroid-specific and nonspecific isozymes. Proc Natl Acad Sci. 1989;86(3):792–6. https://doi.org/10.1073/pnas.86.3.792.
Nakajima O, Takahashi S, Harigae H, Furuyama K, Hayashi N, Sassa S, Yamamoto M. Heme deficiency in erythroid lineage causes differentiation arrest and cytoplasmic iron overload. EMBO J. 1999;18(22):6282–9. https://doi.org/10.1093/emboj/18.22.6282.
Nakajima O, Okano S, Harada H, Kusaka T, Gao X, Hosoya T, Suzuki N, Takahashi S, Yamamoto M. Transgenic rescue of erythroid 5-aminolevulinate synthase-deficient mice results in the formation of ring sideroblasts and siderocytes. Genes Cells. 2006;11(6):685–700. https://doi.org/10.1111/j.1365-2443.2006.00973.x.
Yamamoto M, Nakajima O. Animal models for X-linked sideroblastic anemia. Int J Hematol. 2000;72(2):157–64.
Hunter GA, Ferreira GC. 5-Aminolevulinate synthase: catalysis of the first step of heme biosynthesis. Cell Mol Biol. 2009;55(1):102–10.
Khechaduri A, Bayeva M, Chang HC, Ardehali H. Heme levels are increased in human failing hearts. J Am Coll Cardiol. 2013;61(18):1884–93. https://doi.org/10.1016/j.jacc.2013.02.012.
Zhang F-L, Shen GM, Liu XL, Wang F, Zhao HL, Yu J, Zhang JW. Hypoxic induction of human Erythroid-specific δ-Aminolevulinate synthase mediated by hypoxia-inducible factor 1. Biochemistry. 2011;50(7):1194–202. https://doi.org/10.1021/bi101585c.
Salnikow K, Kluz T, Costa M, Piquemal D, Demidenko ZN, Xie K, Blagosklonny MV. The regulation of hypoxic genes by calcium involves c-Jun/AP-1, which cooperates with hypoxia-inducible factor 1 in response to hypoxia. Mol Cell Biol. 2002;22(6):1734–41. https://doi.org/10.1128/MCB.22.6.1734-1741.2002.
Yamada Y, Pannell R, Forster A, Rabbitts TH. The oncogenic LIM-only transcription factor Lmo2 regulates angiogenesis but not vasculogenesis in mice. Proc Natl Acad Sci. 2000;97(1):320–4. https://doi.org/10.1073/pnas.97.1.320.
Gratzinger D, Zhao S, West R, Rouse RV, Vogel H, Gil EC, Levy R, Lossos IS, Natkunam Y. The transcription factor LMO2 is a robust marker of vascular endothelium and vascular neoplasms and selected other entities. Am J Clin Pathol. 2009;131(2):264–78. https://doi.org/10.1309/AJCP5FP3NAXAXRJE.
Meng S, Matrone G, Lv J, Chen K, Wong WT, Cooke JP. LIM domain only 2 regulates endothelial proliferation, angiogenesis, and tissue regeneration. J. Am. Heart Assoc. 2016;5(10).
Coma S, Allard-Ratick M, Akino T, van Meeteren LA, Mammoto A, Klagsbrun M. GATA2 and Lmo2 control angiogenesis and lymphangiogenesis via direct transcriptional regulation of neuropilin-2. Angiogenesis. 2013;16(4):939–52. https://doi.org/10.1007/s10456-013-9370-9.
A. de L. Paz, Y. M. Pastrana, and L. V. Brandão, “Food deprivation does not affect growth performance of juvenile tambacu,” Acta Amaz, vol. 48, no. 3, pp. 207–210, 2018.
R. Macedo Silva, M. Tavares-Dias, M. W. Reis Dias, M. K. Reis Dias, and R. das Graças Barbosa Marinho, “Parasitic fauna in hybrid tambacu from fish farms,” Pesqui Agropecu Bras, vol. 48, no. 8, pp. 1049–1057, 2013, DOI: https://doi.org/10.1590/S0100-204X2013000800034.
Gomes F, et al. Comparative analysis of the transcriptome of the Amazonian fish species Colossoma macropomum (tambaqui) and hybrid tambacu by next generation sequencing. PLoS One. 2019;14(2):1–23.
Warren AJ, Colledge WH, Carlton MBL, Evans MJ, Smith AJH, Rabbitts TH. The oncogenic cysteine-rich LIM domain protein Rbtn2 is essential for erythroid development. Cell. 1994;78(1):45–57. https://doi.org/10.1016/0092-8674(94)90571-1.
Melo JSC, Pereira JA. Crescimento do Híbrido Tambacu (Fêmea de Colossoma macropomum X macho de Piaractus mesopotamicus) em criação intensiva. 1994;7:59–75.
Von Walden F. Ribosome biogenesis in skeletal muscle: coordination of transcription and translation. J Appl Physiol. 2019;127(2):591–8. https://doi.org/10.1152/japplphysiol.00963.2018.
Wu TD, Watanabe CK. GMAP: a genomic mapping and alignment program for mRNA and EST sequences. Bioinformatics. 2005;21(9):1859–75. https://doi.org/10.1093/bioinformatics/bti310.
Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12(4):357–60. https://doi.org/10.1038/nmeth.3317.
Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30(7):923–30. https://doi.org/10.1093/bioinformatics/btt656.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.
Shannon P. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504. https://doi.org/10.1101/gr.1239303.
Kilkenny C, Browne WJ, Cuthill IC, Emerson M, Altman DG. Improving bioscience research reporting: the arrive guidelines for reporting animal research. PLoS Biol. 2010;8(6):1–2.
We are thankful to UNESP and CIBIO for providing the structure needed. We are thankful to Rondinelle Simões Salomão for carrying out the morphological classification of the three genotypes analyzed.
BEAF and MDPS were supported by FAPESP Proc. 2015/19985–2, Proc. 2016/05009–4 and Proc. 2017/15708–0. MDPS also was supported by CNPq, Proc. 302656/2015–4 as research scholarship. This work was also supported by the Coordination for the Improvement of Higher Education Personnel (CAPES) [Code 001]
Ethics approval and consent to participate
This study was approved by the Ethics Committee on Animal Use – CEEA/UNESP, Institute of Biosciences, Botucatu, SP (Protocol. 1002 - CEUA) and was performed in accordance with the Ethical Principles in Animal Research adopted by the Brazilian College of Animal Experimentation (COBEA). The experiments were also conducted following the ARRIVE guidelines .
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Read number information. Read number throughout filtering process. Data showing the decrease number of reads according to filtering steps performed.
GFF annotation. GFF file containing annotations of mature sequences over precursor sequences of D. rerio.
FASTA file. FASTA file containing precursor sequences of D. rerio.
Raw and normalized counts. Data showing the raw number of reads obtained from each miRNA detected in libraries, and normalized values.
DESeq2 results. Total DESeq2 results and filtered DESeq2 results for TC vs PC contrasts.
DESeq2 results. Total DESeq2 results and filtered DESeq2 results for TC vs TQ contrasts.
DESeq2 results. Total DESeq2 results and filtered DESeq2 results for PC vs TQ contrasts.
DESeq2 results. DESeq2 results showing only differential expressed miRNAs for TC vs PC, TC vs TQ and PC vs TQ contrasts.
Perl script. Perl script used for filtering reads by size.
Flowchart. Bioinformatics pipeline. (PPTX 50 kb)
About this article
Cite this article
Fantinatti, B.E.A., Perez, E.S., Zanella, B.T.T. et al. Integrative microRNAome analysis of skeletal muscle of Colossoma macropomum (tambaqui), Piaractus mesopotamicus (pacu), and the hybrid tambacu, based on next-generation sequencing data. BMC Genomics 22, 237 (2021). https://doi.org/10.1186/s12864-021-07513-5
- Skeletal muscle