Next Article in Journal
The Effect of an Innovative Combination of Bergamot Polyphenolic Fraction and Cynara cardunculus L. Extract on Weight Gain Reduction and Fat Browning in Obese Mice
Next Article in Special Issue
Theoretical Analysis and Expression Profiling of 17β-Hydroxysteroid Dehydrogenase Genes in Gonadal Development and Steroidogenesis of Leopard Coral Grouper (Plectropomus leopardus)
Previous Article in Journal
Serum Calretinin and Genetic Variability as a Prognostic and Predictive Factor in Malignant Mesothelioma
Previous Article in Special Issue
Exploring the Effects of Rearing Densities on Epigenetic Modifications in the Zebrafish Gonads
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Institute of Genome Biology, Research Institute for Farm Animal Biology (FBN), 18196 Dummerstorf, Germany
2
Bioinformatics Core Facility, Oslo University Hospital, 0424 Oslo, Norway
3
Department of Infectious Disease Epidemiology, Bernhard Nocht Institute for Tropical Medicine, 20359 Hamburg, Germany
4
Faculty of Agriculture and Environmental Sciences, University of Rostock, 18059 Rostock, Germany
5
Department of Life Sciences and Health, OsloMet—Oslo Metropolitan University, 0167 Oslo, Norway
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2024, 25(1), 189; https://doi.org/10.3390/ijms25010189
Submission received: 10 November 2023 / Revised: 14 December 2023 / Accepted: 19 December 2023 / Published: 22 December 2023
(This article belongs to the Special Issue Fish Genomics and Developmental Biology)

Abstract

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

1. Introduction

Micro RNAs (miRNAs) are short non-coding RNA molecules with an average length of 22 nucleotides (nt). Mature miRNAs are processed from hairpin precursor miRNAs (pre-miRNAs), which, in turn, are Drosha-processed products of transcribed primary miRNAs (pri-miRNAs) [1,2]. MicroRNAs play an important role in the post-transcriptional regulation of their target transcripts. By binding to the target mRNA, the mature miRNA mediates their translational repression or degradation [3]. From studies in mammals, it is estimated that more than 60% of the mRNAs have conserved miRNA binding sites [4]. MicroRNAs are thus involved in the regulation of multiple developmental and physiological processes like growth, immune response, and stress response [4,5,6,7]. In teleost fish, there are similarly high numbers of predicted target genes [8], and the high conservation of mature miRNA sequences suggests shared functions across species [9]. Data on the miRNA repertoire in teleosts are still limited, with only 16 out of the over 26,000 teleost fish species being present in miRBase (22 October 2023; http://www.mirbase.org/). Nevertheless, research on teleostean miRNAs has increased significantly (reviewed in [9,10,11]), as has research on fish models including zebrafish (Danio rerio) [12,13] and economically relevant species like Atlantic salmon (Salmo salar) [14,15], rainbow trout (Oncorhynchus mykiss) [16,17], Nile tilapia (Oreochromis niloticus) [18,19], Atlantic cod (Gadus morhua) [20,21], lumpfish (Cyclopterus lumpus) [22], and channel catfish (Ictalurus punctatus) [23,24]. In addition, the recently developed FishmiRNA database also provides miRNAs for black bullhead (Ameiurus melas), Japanese medaka (Oryzias latipes), molly (Poecilia mexicana), European perch (Perca flavescens), stickleback (Gasterosteus aculeatus), and blackfin icefish (Chaenocephalus aceratus) [25]. The studies in fish have, like in higher vertebrates, shown that particular miRNAs are associated with immune response, stress response, and growth and are even supposed to regulate the aforementioned processes. Smaller groups of miRNAs respond, for example, to viral and bacterial infection in Atlantic salmon, rainbow trout, Atlantic cod, and Nile tilapia [20,24,26,27,28]. Many of these miRNAs, including miRNA-146, miRNA-155, and miRNA-223, respond in the same manner to infectious diseases across species, and the predicted targets are key immune genes [10,24,26,27]. Similarly, investigations of the miRNAs associated with smoltification and sea water adaptation in salmonids have revealed groups of miRNAs targeting genes enriched in particular stress and immune response pathways [15,18]. Other studies on teleost fishes like tilapia and zebrafish have identified miRNAs associated with osmoregulatory stress and hypoxia [29,30]. Economically interesting traits like growth are likely regulated in fish by particular miRNAs, as indicated in studies on Nile tilapia [18,19]. Certain miRNAs also respond differently to functional feed designed to elevate robustness to infections. Others are associated with familiar differences in lipid metabolism affecting the fatty acid composition of the fish [31]. Such miRNAs influencing relevant traits may be utilized as biomarkers. Alternatively, differences or variations in expression that affect the miRNA-target gene interaction may be utilized for marker-selected breeding to assure that genetically beneficial variation is included in the breeding stock.
Global aquaculture is experiencing an annual growth of >5%; around 57.5 million tons of finfish are currently produced worldwide, while aquaculture finfish production in Europe is stagnating [32,33]. In order to revive the aquaculture market, production is set to shift towards focusing on high-value fish like pikeperch (Sander lucioperca). This predatory freshwater fish belongs to the family of Percidae and is found in habitats in western Eurasia. Due to its excellent flesh quality, displayed by a low fat content and delicate flavor without small intramuscular bones, pikeperch has a high consumer acceptance [34,35]. Together with its relatively rapid weight gain, these quality features make pikeperch an economically interesting species for intensive aquaculture systems. The global aquaculture production of pikeperch reached 3074 t of live weight in 2020 (FIGIS FAO 2022, https://www.fao.org/fishery/en/figis (accessed on 13 December 2023)), with a European market value up to USD 22.2/kg (whole fish) in countries such as Germany and France [36]. The production of this percid was established in recirculating aquaculture systems (RAS) in Europe at the beginning of the twenty-first century [37]. Its domestication through successful breeding in aquaculture has progressed steadily in recent years, but several bottlenecks still fundamentally limit the production efficiency [38]. The main problems are the species’ fluctuating reproductive performance and relatively slow specific growth, the low survival rates of the larvae due to cannibalism, high deformation rates or the incorrect development of the swim bladder, feed conversion during early development, and the relatively high associated production costs [38,39]. Although RAS in particular, with their defined environment, are predestined for breeding approaches [40,41], to date, no marker-assisted selection or phenotype-based family selection approaches have been applied in pikeperch breeding programs, which could help to overcome some of the problems. However, in the EU in particular, there has been an increasing shift in recent years towards family-specific and genetic breeding approaches for economically relevant fish species such as Atlantic salmon [42,43]. The provision of complete genetic tools, including genome sequences, SNP sets, or ncRNAs, is essential for future breeding approaches [44]. The exploration of the genetic resources of the species is the backbone for potential marker-based selection. In addition to its growing commercial importance, research on the genetic and genomic composition of pikeperch has also gained momentum. With recent publications on its genome [45,46] and transcriptome [47], the systematic exploration of functional genomics is feasible in the context of, for instance, relevant commercial traits. So far, nobody has explored miRNA genes and their expression in pikeperch. The identification and characterization of the miRNAome in pikeperch is the next step and is crucial to understanding the regulatory gene mechanisms and subsequently applying this knowledge to improve the performance of economically important traits in this fish species. The characterization of pikeperch-specific miRNAs would allow for differential expression analyses involving applying the pikeperch miRNAome in comparisons of the conditions important to fish health or other economically interesting traits. Such miRNAs and any variation affecting their miRNA-target gene interactions have, as mentioned above, the potential to be used as biomarkers and in marker-assisted breeding to benefit pikeperch aquaculture production. Such beneficial variants may also, in the future, be efficiently introduced in breeding stock by gene editing methods.
The aim of this study was to identify the pikeperch miRNA repertoire by using small RNA high-throughput sequencing (HTS). Samples representing different developmental stages and organs were investigated to discover both evolutionally conserved teleost miRNAs and any novel pikeperch miRNAs. Applying comparative approaches, we also aimed to disclose miRNAs likely involved in the regulation of development or those serving organ-specific functions.

2. Results and Discussion

2.1. Identification of Conserved and Novel miRNA Genes in Pikeperch

Small RNA libraries from 20 pikeperch samples were sequenced. These comprise six different developmental stages and four organs of adult fish. A total of 172.6 million raw reads were obtained. After processing, 94.1 million high-quality and adapter-trimmed reads were left, ranging from 2.5 to 9.2 million reads per sample (Supplementary Table S1).
An initial miRDeep2 analysis and subsequent manual characterization led to the identification of 242 miRNA genes in the samples investigated. Most genes were located at single unique locations in the pikeperch genome, while seven were identical duplicated genes positioned in two or three different loci. Basic local alignment search tool (BLAST) homology searches against the stem-loop sequences in the miRNA database could classify 234 of the miRNA genes (96.7%) as orthologues to evolutionary conserved miRNAs. The miRNA genes identified from pikeperch were assigned to 104 families which exhibited a wide range of family members. The let-7 (12 family members) and miR-430 (10 family members) families host the most members in the pikeperch miRNAome. Eight of the miRNA genes showed no match with any of the miRNA genes in miRBase but exhibited all the properties expected for a miRNA gene from our identification criteria (see Section 3.3). Additional nucleotide BLAST (BLASTN) searches against other RNA databases ruled out the notions that they were other kinds of non-coding small RNAs or derived from repetitive elements. Thus, they are the only novel miRNAs observed in pikeperch so far, and they were annotated accordingly. An overview of all miRNA precursor sequences, their genomic location, the sequences of the mature 5p and 3p miRNAs originating from each precursor, and the 5p or 3p abundance differences (arm dominance) is given in Supplementary Table S2. The resulting pikeperch miRNAome consists of 375 unique mature miRNA (miRs) sequences (185 mature 5p and 190 mature 3p) (Supplementary File S1). Many of the mature miRNAs were identical, even if processed from different miRNA genes (precursors). The individual annotation in the miRNAome reflects the occurrence of identical mature miRNAs from different loci by adding additional numbers or letters in the suffix of some mature miRNAs. There are, for example, two slu-mir-21 genes (Supplementary Table S2), but the mature 5p is identical in both and is thus annotated as slu-mir-21-1-2-5p in Supplementary File S1. Likewise, the two slu-mir-7 genes have mature 5p’s that can not be distinguished, and the resulting mature miRNA was consequently termed slu-mir-7ab-5p in the unique miRNAome.
Pikeperch’s miRNA gene diversity, with 242 identified miRNA genes from 104 conserved families, is similar to the miRNA gene diversity described for other teleosts such as channel catfish (237 miRNAs and 105 families; [23]), lumpfish (391 miRNAs and 104 families; [22]), and Atlantic cod (219 miRNAs and 110 families; [21]). This finding in pikeperch is therefore in accordance with prior studies pointing out that most miRNAs are evolutionarily highly conserved cross-species non-coding RNAs [48,49]. In addition, we also identified eight novel miRNAs in pikeperch. In many studies, including ours, previously unknown or species-specific miRNAs have been identified in teleosts in addition to the orthologous miRNAs. The emergence of novel miRNAs is an evolutionary ongoing process involving several mechanisms [50]. Lu et al. (2008), in their study on the fruit fly (Drosophila melanogaster), calculated a high “birth and death rate” for new miRNA genes during evolution [51]. Six of the pikeperch-specific miRNAs (from slu-mir-nov1 and -nov3 to -6 and -nov7) were predominantly expressed in the early fingerling stage; one was also expressed in eyed eggs (slu-mir-nov7), and two novel miRNAs (slu-mir-nov2a and -b) were exclusively abundant in eyed eggs and the early larvae stages (Supplementary Table S3). These expression patterns could indicate that their functions are associated with stage-specific development.
Animal miRNAs may be transcribed as individual miRNA genes or, if located in gene clusters, they may be transcribed coordinately [52,53]. These clusters, as defined by miRbase, are composed of at least two physically adjacent miRNA genes (see Materials and Methods) [54]. We identified 113 miRNA genes (47% of all pikeperch miRNA genes) that were grouped into 45 such distinct clusters. Each cluster consists of up to 13 miRNA genes (Supplementary Table S4). The gene clusters are localized throughout the pikeperch genome and can be divided into 31 homo- and 14 hetero-clusters. All the miRNA gene clusters discovered in pikeperch, as well as the miRNA genes in the clusters, are evolutionarily conserved and have orthologues in other teleost species (according to miRbase and MetaMirClust, https://fgfr.ibms.sinica.edu.tw/MetaMirClust (accessed on 28 September 2023)). Like in other vertebrates [13,14,21,55], the miRNA-430 gene cluster stands out as the one with the largest number of clustered genes.

2.2. Mature miRNA Diversity across Sample Groups and Their Relative Abundance within Sample Groups

Assessing diversity in terms of the different unique mature miRNAs in each sample group revealed the presence of 371 to 375 (all) mature miRNAs in eyed eggs, as well as larvae 2 days post hatching (dph), 15 dph, and 26 dph. The spleen from 3 g- and 35 g-fingerlings contains 363 and 350 miRAs, respectively. The adult organs varied from 325 to 368 different mature miRNAs (Table 1). However, all miRNA genes were expressed at least at a basic level in all samples, as there were at least one mature miRNA (either 5p or 3p) detected from each of the 242 miRNA genes. The exception from this were slu-miRNA-1788 and the two novel miRNAs—slu-miRNA-nov2a and b—that were not detected in spleen fingerling 35 g and liver. Liver also lacked the expression of slu-miRNA-726 and the novel miRNA gene slu-miRNA-nov5. These three novel miRNAs were not detected in gills that also lacked slu-miRNA-430a. A complete overview of normalized read counts for all unique mature miRNAs in all samples is given in Supplementary Table S3.
The detection of which mature miRNAs are highly abundant within each sample could provide more information about which miRNAs are biologically important in each group. Such a relative ranking of mature miRNAs by their normalized reads within samples are not hampered by any technical bias as the read numbers used to show the relative abundance are retrieved from the same sample. The read numbers are easily available and can thus directly reflect what the relative abundances of the mature miRNAs in a given sample are. The abundance patterns differed somewhat between (i) the early development from eyed eggs to larvae 26 dph (Figure 1), (ii) the spleen samples from the adult and the two fingerling stages (Figure 2), and (iii) the other organ samples of the adult fish (Figure 3). The top enriched mature miRNAs could be divided into miRNAs that were common to all groups and miRNAs that were highly enriched in certain developmental stages or organs of pikeperch. Figure 1, Figure 2 and Figure 3 display the top 20 most abundant mature miRNAs per sample type (Supplementary Table S5). Forty-nine different mature miRNAs from thirty-seven different families dominated the different samples, representing between 62.3% and 82.8% of the total amount of mature miRNAs per developmental stage or organ. Some mature miRNAs, like slu-mir-21-1-2-5p, slu-mir-22ab-3p, and slu-mir-26a-1-3-5p, were ubiquitously enriched throughout the stages. These mature miRNAs also had similar high abundances in all organs and developmental stages in other teleost fish [14,22]. Together, this indicates housekeeping functions that are conserved across species.
While the composition of the top 20 most abundant mature miRNAs in (i) was quite similar between all four developmental stages (Figure 1), some of these mature miRNAs displayed a relatively high abundance in only these four early developmental stage groups (eyed egg to 26 dph), not the other groups ((ii) and (iii)). These were two slu-mir-10 family members (slu-mir-10b and -c-5p), slu-mir-9-1-4-5p, slu-mir-7ab-5p, slu-mir-1-3p, and slu-mir-206-3p. This could reflect their importance in the regulation of early development. The findings from an analysis of Atlantic salmon support the notion that the mature slu-mir-10b and -c-5p have such important developmental functions, as they constituted more than 50% of the total number of miRNAs in the early developmental stages of Atlantic salmon [14]. In line with this is the finding that they regulate hox genes in Nile tilapia embryonic development [56] and also mediate cell proliferation and differentiation [57]. The mature members slu-mir-9-1-4-5p and slu-mir-7ab-5p have been reported as highly abundant in early developmental stages, as well as in the brain [14,22,58,59,60,61]. It is likely that their high relative abundance in early developmental groups in pikeperch reflect similar functions. Studies have also revealed that orthologs of slu-mir-1-3p and slu-mir-206-3p are important in muscle development and embryonic development [62,63], which would also agree with their high relative abundance observed in the pikeperch developmental groups.
Slu-mir-143-3p was by far the most abundant mature miRNA in the two spleen samples from the fingerling stages, making up 1/6 of the total miRNA read counts in juvenile pikeperch (~17 and ~18%, Figure 2). This indicates that it has a relative abundance in the fingerling spleen that is more than three times higher than what is seen in any other organ. A similarly high abundance was detected in adult spleen (11%). In addition, slu-mir-144-5p displayed a high relative abundance in organ samples from the spleen, while slu-mir-8159-5p, slu-mir-101a-3p, and both of the mature miRNAs processed from slu-miRNA-2188 were predominantly enriched in the spleen of the juvenile fish. The phase of larvae and fingerling stages in fish is correlated with high growth rates and elevated energy demands. Studies on the regulatory role of the miR-143 family in human and fish link it to lipid metabolism [64,65]. In rainbow trout, an increase in miR-143 abundance during early development up to the juvenile stage has been demonstrated [65]. Mir-144-5p shows similarly high abundance in Atlantic cod spleen [21] and is associated with erythropoiesis [66]. Together with miR-101-3p, it is associated with macrophage function [67]. miRNA-2188, which is teleost-specific, is associated with immune responses and monocyte to macrophage maturation in fish [10,27], while the other teleost-specific miRNA that was highly abundant in the fingerling spleen, miR-8159-5p, is associated with immune response in fish [26,68].
The relative abundance of some mature miRNAs in adult pikeperch seems to differ between organs (Figure 3). Remarkably, slu-mir-122-5p accounted for more than one fourth of the total miRNA read count in the liver. Similarly high abundance has been seen in liver samples in all vertebrates, including teleost fish [14,21,22,69]. MicroRNA-122-5p was demonstrated to be a key regulator of cholesterol metabolism in the liver [69]. In addition, slu-mir-192-5p was among the top 20 abundant mature miRNAs in the liver and in the larva developmental stages. Such a high liver abundance has also been reported in other teleost fish [14,21]. In agreement with the notion of it having both developmental and liver-specific functions, it has been reported that miR-192-5p is abundant in the liver of mammals, where it promotes the development of the liver and coordinates energy metabolism [21,70].
Slu-mir-205-5p is highly abundant only in gills. This gill-specific enrichment is also conserved in other fish species [22,29,71]. However, its function in gills is not well known. Similarly, slu-mir-203a-3p is highly abundant in gills and in the developmental stages. Again, this has been reported in other fish species [14,22], but its function is not known. The observed high variation in the relative abundance of mature miRNAs between the different developmental stages and organs suggest that miRNAs are important in organ and tissue differentiation and/or the maintenance of their functions. The conserved miRNAs showed similar abundance patterns in other species, as observed in pikeperch, strongly suggesting that they serve the same conserved functions in pikeperch as well.

2.3. Differential Expression Analysis and qPCR Revealed Additional miRNAs Associated with Particular Life Stages and Organs

The biologically important mature miRNAs at particular life stages or in the different organs are not necessarily the highly abundant ones. Thus, to further investigate whether some mature miRNAs are differentially expressed during pikeperch development beyond what was obvious from the distribution of the most abundant ones, we performed DESeq2 expression analyses including comparisons between (i) the early developmental stages, (ii) the two fingerling stages, and (iii) the adult organs.
The selection of the differentially expressed mature miRNAs detected by the DESeq2 analyses was validated by qPCR, along with some mature miRNAs identified as particularly abundant in a given sample (see Section 2.2) (Supplementary Table S6). To identify stable reference miRNAs, the expression levels of four candidate reference miRNAs were analyzed by qPCR in all 48 samples, and RefFinder suggested that slu-mir-30c-5p is the most stable single gene (stability value of 1.32) and that slu-mir-30c-5p/slu-mir-30-1-2-5p was the best two-gene combination (stability value of 0.58). The two reference miRNAs mir-107-3p and mir-26a-1-3-5p, with proven high stability in Atlantic salmon, were rated as less stable across the samples tested in pikeperch [14,72]. The two genes in the two-gene combination with superior stability were therefore subsequently used as reference genes for our qPCR analysis. The sample-specific high abundances and significant differences from the DESeq2 analysis were mostly confirmed in the qPCR analysis (Figure 2, Supplementary Tables S7 and S8). Furthermore, there was a strong positive correlation between the results from DESeq2 and qPCR (Pearson’s correlation coefficient, r = 0.9) (Figure 4a–c). This strong correlation indicates that the other DESeq2 results, although not validated by qPCR, are also reliable.
A differential expression analysis, facilitated by DESeq2 and involving the comparison of the very early stages and larvae (eyed egg and 2 dph vs. 15 dph and 26 dph), led to the identification of 98 differentially expressed mature miRNAs (Table 2, Supplementary Table S7). This large number was not unexpected, considering the major anatomical changes, the specialization of tissue, and the development of organs occurring in this early phase [28,73,74]. The very early stages revealed 44 mature miRNAs from 25 families, with significantly higher expression compared to the larvae stage. Among these were slu-mir-9-1-4-5p and the miRNA 10 family of mature 5p miRNAs. These displayed a generally higher relative abundances in all the developmental stages compared to the adult organ groups (see Section 2.2). Here, the analysis between early developmental stages also revealed that their expression levels were highest in the very early stages. Another family standing out as highly expressed was the miRNA 430 family (fold change > 50×). This was not unexpected, as they are known to promote the removal of maternal mRNAs [55], and similar expression patterns were also observed in Atlantic salmon in [14]. The clustered miRNAs 25, 93, and 106 also showed an increased expression of their mature miRNAs at this very early stage. This miRNA cluster is known to control the epithelial to mesenchymal transition in vertebrates [75,76]. Another miRNA gene cluster, 17-19-92, was significantly enriched at this very early stage. Similarly high expression was reported in vertebrate embryos in [77], wherein they were associated with skeletal development. The deletion of the two clusters 106-25-93 and 17-19-92 is lethal to embryos in mammals [77,78].
On the other hand, the let-7-5p mature miRNAs was strongly increased in expression (fold changes from 3–20×) in the larvae stage compared to the very early stages. This is in agreement with their function as regulators of the entry into the larvae stage [74,79,80]. We also noted a higher expression in the larvae of slu-mir-133ab-3p, which are clustered with miR-1 and miR-206 (Supplementary Table S4) and act together with these miRNAs to control muscle development (see Section 2.2) [81]. Both the 5p and 3p mature miRNAs of the slu-mir-199 genes were highly increased in the larvae stages, in accordance with findings for other teleost fishes ([14,21,22,82,83]). This miRNA controls sonic hedgehog signaling, which is essential to normal embryonic development [84]. It is clustered with the miRNA-214 genes that are also enriched in larvae. These two miRNAs have also been reported as important in skeletal development in zebrafish [85]. Four other mature miRNAs—the clustered mature miRNAs slu-mir-143-3p and slu-mir-145-5p and the clustered matures slu-mir-192-5p and slu-mir-194b-5p—were also particularly increased in the larvae stage. Again, this seems to be conserved in fish [14,22], and their function in vertebrate development is associated with cardiac progenitor cell differentiation and muscle development, respectively [86,87]. One teleost-specific miRNA, slu-mir-8159-5p, was increased. Despite being associated with immune system and stress response [26], the particular function of slu-mir-8159-5p, at this stage, is not known.
The early development of teleost fish, from fertilized egg, through the larval stages, to juvenile fish, involves a series of major anatomical changes, including, for example, the growth of muscles and the skeleton, the appearance of fins, the development of the nervous system, and the adjustment of the digestive system (as reviewed in [73]). MicroRNAs are instrumental in regulating these anatomical changes, exhibiting high levels of tissue- and cell-specific abundance [28]. Many of the conserved pikeperch mature miRNAs, as well as many of the clustered genes, displayed similar enriched expression profiles, as in other vertebrates at embryo/larvae life stages. This finding suggests they also share the same important regulatory functions in early development, as revealed in other vertebrates. The joint increase in the mature pikeperch miRNAs annotated as clustered genes (Supplementary Table S4) supports the notion that they are transcribed as polycistrons. The detection of the similar expression changes in such polycistrons also lends credence to the expression analysis itself, as each of the mature miRNAs were independently analyzed in the DESeq2 pipeline.
When comparing miRNA expression in the two fingerling stages from spleen, 13 upregulated miRs from 10 families were found in the early stage of 3 g fingerlings. Some of these were also seen in the relative abundance investigation (Section 2.2). Interestingly, one of these, slu-mir-nov6-3p, was among the novel pikeperch miRNAs characterized in this study. Two other miRNAs, slu-mir-2188-5p and slu-mir-8159-5p, are associated with immune system and pathogen responses in several fish species [20,26,68], while a third highly expressed miRNA, slu-mir-216a-1-2-5p, is a conserved mature miRNA that is reportedly associated with viral infection in rainbow trout spleen [88]. The late stage of 35g fingerlings spleen depicted 6 miRNAs from different families as highly expressed. The significantly increased expression of slu-mir-1388-3p and slu-mir-731-5p was also confirmed by qPCR (Figure 2). Orthologs of all these six mature miRNAs (slu-mirs-1388-3p, -730-5p, -7132-5p, -731-5p, -125-1-3-5p, and -99b-5p) have also been associated with pathogen response in fish [10,26,27,88,89]. In summary, most miRNAs with particularly high expression in one of the two fingerling spleen samples were associated with the function of spleen as an important immune organ in fish [90].
The comparisons between adult organs were not only in agreement with the relative abundance distributions (Section 2.2) but also revealed some additional miRNAs that were preferentially expressed in some organs (Supplementary Table S7). There was, for example, an up to 165-fold higher expression of mir-10 family members in the head kidney (slu-mir-10b-5p, -10c-5p) compared to the other organs (Table 2). The enrichment of this family in this central immune organ of bony fish has already been demonstrated in studies on codfish (Gadidae) [21], snakehead (Channidae) [91], and carps (Cyprinidae) [67,92]. The mir-10 family seems to play an important regulatory role in the inflammation and innate immune response of vertebrates and invertebrates through the post-transcriptional regulation of Toll-like receptor signaling and associated transcription factors [93,94].
The abundance of the three most highly expressed mature miRs in the gills (slu-mir-205-5p, slu-mir-203b-3p, and slu-mir-725-3p) was up to 440-fold higher compared to the expression levels in the other organs (Table 2). This organ-specific enrichment is consistent with the findings of many other studies on teleosts and mammals [21,22,29,69,95]. The DE of the dominant 5p-arm of slu-mir-205-5p in the gills, as well as slu-mir-122 in the liver, was also confirmed by qPCR (Figure 3). The liver-enriched miR-122 is, as mentioned in Section 2.2, a key regulator of lipid metabolism in mammals and fish. Menningen et al. (2014) concluded that the involvement of miR-122 in liver-specific metabolic processes could have arisen evolutionarily early in vertebrates [96].

2.4. The miR-462/miR-731 Cluster Is Associated with Hypoxia Stress in Pikeperch

The teleost-specific cluster miR-462/miR-731 is involved in pathogen-associated immune responses and the response to hypoxia stress in bony fish [10,30,97,98]. Previously, we studied the influence of chronic oxygen deficiency on the immune status of pikeperch [99]. To further evaluate the possible conservation of the hypoxia-related functions of the miR-462/miR-731 cluster in pikeperch (Supplementary Table S4), we investigated the expression of the clustered miRNAs under low-oxygen saturation. Maintaining adult pikeperch for one day at low dissolved oxygen (DO) levels (±3.2 mg/L DO) significantly increased the miR-462/miR-731 cluster expression by 0.4-fold compared to that at normoxic levels (±8.3 mg/L DO) (Figure 5). Huang et al. (2015, 2017) found a similar HIF1α-induced expression of the respective cluster members in zebrafish larvae and adult blunt snout bream (Megalobrama amblycephala) 4–24 h after exposure to hypoxic conditions [30,97].

3. Materials and Methods

3.1. Fish Sample Material

The early-stage pikeperch used in this study were bred and reared in RAS at the State Research Centre for Agriculture and Fisheries Mecklenburg-Vorpommern (Hohen Wangelin, Germany); the fingerlings and adult individuals were bred and reared at the Experimental Animal Facility for Aquaculture of the Institute for Farm Animal Biology (FBN, Dummerstorf, Germany). Pikeperch progeny was reared in a RAS system with a total volume of 9 m3. Larvae were kept in fish tanks with a volume of 0.5 m3 at an initial stocking density of 100 larvae/L with a water exchange rate of 20%/h/tank, water temperature of ~15.7 °C, and ~9.2 mg/L dissolved oxygen (DO). The water quality was ensured by continuous purification, UV disinfection, and the daily monitoring of relevant water quality parameters. The photoperiod during hatchery was set at 24 h light (L): 0 h dark (D) until day 45 and subsequently at 17 L:7 D (+1.5 h dusk and dawn). The fingerlings and adult fish were kept in 160 L round tanks at 21.1 °C and ~8.1 mg/L DO and subjected to a 12 L:12 D cycle until sampling.
A total of 52 samples were used for miRNA sequencing and quantitative expression analysis using real-time quantitative PCR (qPCR) (Supplementary Table S1). Samples were selected to cover as many life stages as possible, taking into account the known critical phases of pikeperch rearing mentioned in the Introduction section. This included eyed eggs (six pools of n = 10 each), yolk sac larvae (two days post hatching (dph) and 15 dph: each six pools of n = 3), and larvae (26 dph: six pools of n = 3 each). Furthermore we sampled spleen from two fingerling stages (3g and 35g: n = 6 each), as well as the liver, gills, head kidney, and spleen of four 19-month old adult individuals each (n = 4 each). The sampling and killing methods followed the standards described in the German Animal Welfare Act [§ 4(3) TierSchG]. In detail, the pikeperch were euthanized with an overdose of 2-phenoxyethanol (50 mg/L water), followed by a blow to the head and spine sectioning at the skull level. The collected material was snap-frozen in liquid nitrogen and stored at −80 °C until further investigation. Of the total 52 samples, 2 samples of the developmental stages and adult organs were used for the high-throughput sequencing of the miRNAs (total of 20 samples). The remaining four samples from each of the different developmental stages were used for qPCR expression analysis (total of 24 samples), together with the four organ samples from each of the adult fish (total of 16 samples, including the two samples of each organ also used for HTS (Supplementary Table S1).
In addition, we used eight additional organ samples collected in a prior study by our working group to analyze the expression of the miR-462/731 cluster in response to hypoxia stress [99]. Four pikeperch spleen samples subjected to normal (±8.3 mg/L dissolved oxygen (DO)) or low-oxygen conditions (3.2 mg/L DO) for one day were used for our expression analysis (Supplementary Table S1).
The dissection and sampling protocol was approved by the Committee on the Ethics of Animal Experiments of Mecklenburg-Western Pomerania (Landesamt für Gesundheit und Soziales LAGuS; approval ID: 7221.3-1-009/19).

3.2. Small RNA Isolation, Library Preparation, and Sequencing (Ilumina Technology)

Total RNA was isolated from the 20 samples by using the mirVana miRNA isolation kit (Thermo Fisher Scientific, Bremen, Germany) following the manufacturer’s protocol. The integrity and quality of the total RNA was evaluated using Nanodrop and an Agilent Bioanalyzer 2100 platform (Agilent Technologies, Waldbronn, Germany), employing the Small RNA assay (Agilent Technologies). Only RNA samples with RNA integrity values (RIN value) ≥ 7 were further processed (Supplementary Table S1).
The small RNA libraries were constructed by using the NEBnext® multiplex small RNA Library Prep Set (New England Biolabs, Inc., Ipswich, MA, USA) according to the manufacturer’s protocol. One ug of total RNA was used as the input for the preparation of the libraries. The library preparation included 5′ and 3′ adapter ligation, reverse transcription, and a size selection of 140–150 bp fragments using 6% polyacrylamide gel. An Illumina NextSeq 500 providing 75 bp single end reads was used for sequencing, which was carried out at the Norwegian Sequencing Centre (NSC).
Raw data from the deep sequencing of the samples have been submitted to the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) database under the BioProject accession number PRJNA977610.

3.3. Pre-Processing and miRDeep2 Analysis of Sequencing Data

The quality of the sequenced small RNA raw reads was checked using FastQC (v.0.11.8). The adapter sequence (5′ AGATCGGAAGAGCACACGTCTGAACTCCAGTCAC 3′) was removed using Cutadapt (v.2.3) Python Package (v.3.7.3) [100]. This was followed by the size filtering of the raw reads. Reads that were shorter than 18 nts or longer than 25 nts were discarded. An additional FASTQC analysis after size filtering ensured that the per base sequence quality had phred scores of 32 or more in all samples. The processed high-quality reads were used for downstream analysis applying the miRDeep2 software package (v.2.0.0.7) [101]. The present version of the pikeperch genome assembly (GeneBank accession number: GCA-008315115.2) was used as a reference genome for miRNA identification. All samples were processed independently. During processing, one liver sample was removed due to the uncertainty of its origin (possible contamination). The workflow applied to identify conserved and novel pikeperch miRNA sequences as described in Woldemariam et al. (2019) [14]. In short, it comprises the following: The miRDeep2 tool-specific log-odds score (miRDeep2 score) of ≥1.0 × 102 was used as a cut-off for the detection of the miRNA precursor and the two mature miRNAs (5p and 3p) to prevent the inclusion of false positives. All those with scores above the threshold were thus regarded as putative precursors with their mature miRNAs processed as expected from Dicer cleavage. They were additionally inspected to fulfil the following criteria: (1) detected in at least two independent samples, (2) at least 10 sequence reads of 3p and 5p mature miRNAs mapped without mismatches to the hairpin precursor, (3) the mature miRNAs were paired with the precursor hairpin with 0–4 nt overhang at their 3′ ends, (4) the reads mapped supported a consistent pre-processing of 5′ end (5′ homogeneity), and (5) more than 60% of the bases of the mature sequences paired in the hairpin structure [54]. All miRNA precursor sequences were blasted against the miRNA database (miRNBase, version 22.1) (https://www.mirbase.org/search.shtml (accessed on 8 March 2022)). In the case of returning a significant hit (e-value of ≤1 × 10−5) against any hairpin precursor in the database, they were annotated as a conserved pikeperch miRNA orthologue and further annotated by us according to the miRbase nomenclature guidelines (slu- prefix and the same number as the orthologues in other teleosts) [54,102]. The precursors that fulfilled our selection criteria but did not significantly match any of the miRNbase entries were considered as putative novel miRNAs. They were further analyzed via blastn searches against known RNA databases in GeneBank (http://blast.ncbi.nlm.nih.gov/Blast (accessed on 15 March 2022)), Rfam v.14.9 (https://rfam.xfam.org/search (accessed on 15 March 2022)), RefSeq (Last update: 4 December 2020 https://www.ncbi.nlm.nih.gov/refseq/), and fRNAdb v.4.0 (https://dbarchive.biosciencedbc.jp/en/frnadb/desc.html (accessed on 15 March 2022)). If any of these putative novel precursors showed significant hits against any other kind of smallRNAs, they were removed from the dataset. Finally, the remaining putative novel miRNAs were used as queries in blastn searches against the pikeperch reference genome. Any putative precursor that provided significant BLAST hits (E-value ≤ 1 × 10−5) against multiple loci (>5) in the pikeperch genome were considered repetitive sequence elements and discarded from our dataset. Those that passed all these filtering fit the criteria expected for miRNAs and, consequently, were annotated as novel pikeperch miRNAs.

3.4. Identification of miRNA Gene Clusters

We identified miRNA clusters in pikeperch by comparing the precursor locations within contigs according to the following miRBase definition: (1) two or more miRNAs located no more than 10 kilo bases (kb) apart, (2) the miRNAs are transcribed in the same orientation, and (3) are not separated by a transcription unit in the opposite orientation.

3.5. Prediction of Enriched and Putative Differently Expressed (DE) miRNAs in and among Developmental Stages and Organs

Comparisons to indicate the mature miRNAs highly expressed in certain organs or between developmental stages were carried out using the miRNAome from this study as described in [22]. In short, processed reads from all 20 samples were mapped to the unique mature miRNAs in the miRNAome characterized in this study by applying STAR aligner software (v.2.5.2b). The index for mapping was thus generated from the unique mature pikeperch miRNAs discovered in this study (Supplemental File S1) with the parameter genomeSAindexNbases 6. STAR aligner software (v.2.5.2b) with alignIntronMax1 and default parameters was then used for the mapping of reads. Next, the STAR mapping BAM output files were processed further in R-Studio by using the featureCounts function from the Rsubread package (v.1.34.2) to produce count matrices. The count tables were used as inputs in the DESeq2 R package (v.1.24.0) for differential expression analysis.
For statistical valid expression analysis, we pooled counts from eyed eggs and larvae 2 dph (n = 4) and larvae from 15 dph and 26 dph, (n = 4), respectively. For a comparison of the two fingerling stages (n = 2 vs. n = 2), we combined the DESeq2 prediction, with verification being carried out via subsequent qPCR analysis. To improve the reliability of the DESeq2 comparisons between the organ samples from adult fish, each organ was analyzed against a pool of the other three (n = 2 vs. n = 6), and some results were additionally verified by qPCR. Liver was not included in the DESeq2 analysis as one of the two samples did not meet our quality measures.
Only miRNAs with normalized read calls of ≥30 in the samples compared were included in the analysis. Differential expression between the samples was tested using Wald’s log2 fold change test. MiRNAs with a log2 fold change threshold of at least ≤−1.0 or ≥1.0 and a Benjamini–Hochberg adjusted p-value of ≤0.05 were considered differently expressed. It is well known that, often, only one of the mature miRNAs originating from a precursor is incorporated into Argonaut (RISC complex), while the other is simply degraded. This leads to large disproportions in the read numbers of two mature miRNAs (5p and 3p) originating from the same precursor [103]. The one incorporated into Argonaut is often referred to as the guide-miRNA, while the other is named a passenger-miRNA. As the abundance of such mature pairs from the same precursor was easily accessible in our data, we also filtered out the ones with a clear indication that they were the mature passenger miRNA. This filtering used the following criteria: if a mature miRNA was 10 times less abundant than the dominant, assumed biologically relevant guide miRNA, it was not reported among our differentially expressed mature miRNAs (Supplementary Table S7).

3.6. Stability Testing of Possible Reference miRNAs and RT-qPCR Analysis

The normalized read counts of the mature pikeperch miRNAs were compared across samples to search for reference genes suitable for use as normalizers in the subsequent quantitative real-time PCR analysis (qPCR). The relative standard deviation (RSD, coefficient of variation) of each mature miRNA was calculated and used as a measurement of their stability across the samples used for deep sequencing. The four predicted reference miRNAs with a moderate or high abundance (between 2500 and 130,000 normalized read counts) and a low RSD (≤0.3)—slu-mir-30c-5p, slu-mir-107-3p, slu-mir-26a-1-3-5p, and slu-mir-30-1-2-5p—were chosen, and 48 samples were analyzed by qPCR; the results were used for subsequent stability testing using the analytic tool RefFinder (http://www.leonxie.com/referencegene.php (accessed on 14 September 2022)) [104]. This software combines the output of the programs gNorm [105], Normfinder [106], BestKeeper [107], and Delta-Ct [108]. RefFinder builds the geometric means of the individual rankings of the four tools and thereby classifies the suitability and reliability of the candidate reference miRNAs. The analysis classified slu-mir-30c-5p as the most stable single gene (stability value of 1.32) and slu-mir-30c-5p/slu-mir-30-1-2-5p as the best two-gene combination (stability value of 0.58). This two-gene combination was used as a reference in the qPCR analysis.
The ten mature miRNAs (slu-mir-205-5p, slu-mir-122-5p, slu-mir-1388-3p, slu-mir-731-5p, slu-mir-1388-5p, slu-mir-551-p, slu-mir-202-5p, slu-mir-10b-5p, slu-mir-130a-5p, and slu-mir-130ab-3p) that were predicted to be significantly more expressed or enriched in one or more of the samples were selected for further expression analysis by real-time qPCR to verify the DESeq2 results. In addition, we analyzed the expression profiles of the mature miRNAs slu-mir-462-5p and slu-mir-731-5p to investigate the expression of the miRNA cluster miR-462/731 during hypoxia stress. A list of the primers used in this study is provided in Supplementary Table S6.
Expression analyses were performed on a LightCycler96 system (Roche Diagnostics, Mannheim, Germany) using the SensiFASTTM One-Step qPCR kit (Bioline, Luckenwalde, Germany) in line with the manufacturer’s instructions. Therefore, total RNA (500 ng) was polyadenylated and reverse-transcribed using the miRNA 1st-Strand cDNA Synthesis Kit (Agilent Technologies) according to the manufacturer’s instructions. The resulting cDNA was diluted with 150 µL of RNase-free water and stored at −80 °C until further use. A RT-qPCR reaction was performed using a miRNA-specific forward primer in combination with a universal reverse primer (Agilent Technologies) and 35 cycles of the following reaction conditions: initial denaturation at 95 °C for 5 min, denaturation at 95 °C for 15 s, annealing at 54/60 °C for 10 s, and extension at 72 °C for 20 s. Changes in miRNA expression were calculated based on the ∆Cq quantification model, including the calculation of the primer efficiencies and relative quantity (RQ) normalization using the geometric mean of the two selected reference genes, slu-mir-30c-5p/slu-mir-30-1-2-5p [105,109]. Samples with Cq values higher than 30.0 were excluded from the downstream calculations. For size and quality control, the PCR products were verified via gel electrophoresis. A melting curve analysis validated the amplification of distinct products. Statistical significance (p < 0.05) was evaluated by performing a parametric t-test.

3.7. Statistical Analysis

The qPCR data were analyzed using LightCycler 96 software v. 4.0.1. The graphs were generated using GraphPad Prism (version 9.5.1) software. qPCR data are presented as the means + SEM containing four replicates.

4. Conclusions

In this study, we identified the pikeperch miRNA repertoire, which comprises 234 conserved and 8 novel miRNAs. We profiled their abundance during development (from eyed eggs to the organs of adult fishes). While the expression profiles of many miRNAs were similar in all materials, indicating that they regulate housekeeping cellular functions, others showed developmental and organ-specific expression, indicating the regulation of specialized biological functions. The many similarities in expression and the abundance of certain miRNAs at particular developmental stages or organs with other teleosts evidence that these miRNAs maintain important shared evolutionary conserved regulatory functions. The knowledge gained on the miRNAs in this species paves the way for the further exploration of individual miRNAs by qPCR or HTS analysis to disclose their association with traits important to aquaculture production. The identification of miRNAs associated with traits such as development and growth could contribute to tackling the current problems connected to the critical life stages during pikeperch rearing in aquaculture. These specific miRNAs have the potential to be applied as biomarkers in breeding programs involving the application of marker-assisted selection, especially if family-specific expression differences or genetic variation affect the miRNA–target gene interaction. This could further be used for targeted modifications, such as gene editing (e.g., CRISPR), that may effectively modify traits of interest.

Supplementary Materials

The supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms25010189/s1.

Author Contributions

Conceptualization, T.G. and R.A.; validation, M.V., J.A.N. and R.A.; formal analysis, M.V., T.V., J.A.N., A.R. and R.A.; investigation, M.V. and T.V.; writing—original draft preparation, M.V.; writing—review and editing, M.V., A.R., T.G. and R.A.; visualization, M.V. and A.R.; supervision, T.G. and R.A.; funding acquisition, R.A. and T.G. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the European Maritime and Fisheries Fund (EMFF) and the Ministry of Agriculture and the Environment in Mecklenburg-Western Pomerania (Grant #: MV-II.1-RM-001). This research was also funded by The Research Council of Norway, grant number 280839/E40 (awarded to Rune Andreassen). The publication of this article was funded by the Open Access Fund of the Research Institute for Farm Animal Biology (FBN).

Institutional Review Board Statement

The animal protocol was approved by the Committee on the Ethics of Animal Experiments of Mecklenburg-Western Pomerania (Landesamt für Gesundheit und Soziales LAGuS; approval ID: 7221.3-1-009/19).

Informed Consent Statement

Not applicable.

Data Availability Statement

All sequenced samples have been submitted to the NCBI Sequence Read Archive Centre (SRA) (https://www.ncbi.nlm.nih.gov/sra (accessed on 30 May 2023)) under Bioproject PRJNA977610.

Acknowledgments

We would like to acknowledge J. Krinitskij, A. Theel, and O. Seidel (FBN, Dummerstorf) for their expert technical assistance. We are grateful to M. Stüeken (State Research Centre for Agriculture and Fisheries Mecklenburg-Vorpommern, Hohen Wangelin, Germany) for providing pikeperch.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bartel, D.P. MicroRNAs: Genomics, Biogenesis, Mechanism, and Function. Cell 2004, 116, 281–297. [Google Scholar] [CrossRef] [PubMed]
  2. Finnegan, E.F.; Pasquinelli, A.E. MicroRNA biogenesis: Regulating the regulators. Crit. Rev. Biochem. Mol. Biol. 2013, 48, 51–68. [Google Scholar] [CrossRef] [PubMed]
  3. Guo, H.; Ingolia, N.T.; Weissman, J.S.; Bartel, D.P. Mammalian microRNAs predominantly act to decrease target mRNA levels. Nature 2010, 466, 835–840. [Google Scholar] [CrossRef] [PubMed]
  4. Friedman, R.C.; Farh, K.K.-H.; Burge, C.B.; Bartel, D.P. Most mammalian mRNAs are conserved targets of microRNAs. Genome Res. 2009, 19, 92–105. [Google Scholar] [CrossRef] [PubMed]
  5. Bushati, N.; Cohen, S.M. microRNA functions. Annu. Rev. Cell Dev. Biol. 2007, 23, 175–205. [Google Scholar] [CrossRef] [PubMed]
  6. Stefani, G.; Slack, F.J. Small non-coding RNAs in animal development. Nat. Rev. Mol. Cell Biol. 2008, 9, 219–230. [Google Scholar] [CrossRef] [PubMed]
  7. Sonkoly, E.; Ståhle, M.; Pivarcsi, A. MicroRNAs and immunity: Novel players in the regulation of normal immune function and inflammation. Semin. Cancer Biol. 2008, 18, 131–140. [Google Scholar] [CrossRef]
  8. Ramberg, S.; Andreassen, R. MicroSalmon: A Comprehensive, Searchable Resource of Predicted MicroRNA Targets and 3′UTR Cis-Regulatory Elements in the Full-Length Sequenced Atlantic Salmon Transcriptome. Non-Coding RNA 2021, 7, 61. [Google Scholar] [CrossRef]
  9. Bizuayehu, T.T.; Babiak, I. MicroRNA in Teleost Fish. Genome Biol. Evol. 2014, 6, 1911–1937. [Google Scholar] [CrossRef]
  10. Andreassen, R.; Høyheim, B. miRNAs associated with immune response in teleost fish. Dev. Comp. Immunol. 2017, 75, 77–85. [Google Scholar] [CrossRef]
  11. Herkenhoff, M.E.; Oliveira, A.C.; Nachtigall, P.G.; Costa, J.M.; Campos, V.F.; Hilsdorf, A.W.S.; Pinhal, D. Fishing into the MicroRNA Transcriptome. Front. Genet. 2018, 9, 88. [Google Scholar] [CrossRef] [PubMed]
  12. Presslauer, C.; Tilahun Bizuayehu, T.; Kopp, M.; Fernandes, J.M.O.; Babiak, I. Dynamics of miRNA transcriptome during gonadal development of zebrafish. Sci. Rep. 2017, 7, 43850. [Google Scholar] [CrossRef] [PubMed]
  13. Chen, P.Y.; Manninga, H.; Slanchev, K.; Chien, M.; Russo, J.J.; Ju, J.; Sheridan, R.; John, B.; Marks, D.S.; Gaidatzis, D.; et al. The developmental miRNA profiles of zebrafish as determined by small RNA cloning. Genes Dev. 2005, 19, 1288–1293. [Google Scholar] [CrossRef] [PubMed]
  14. Woldemariam, N.T.; Agafonov, O.; Høyheim, B.; Houston, R.D.; Taggart, J.B.; Andreassen, R. Expanding the miRNA Repertoire in Atlantic Salmon; Discovery of IsomiRs and miRNAs Highly Expressed in Different Tissues and Developmental Stages. Cells 2019, 8, 42. [Google Scholar] [CrossRef] [PubMed]
  15. Shwe, A.; Krasnov, A.; Visnovska, T.; Ramberg, S.; Østbye, T.-K.K.; Andreassen, R. Differential Expression of miRNAs and Their Predicted Target Genes Indicates That Gene Expression in Atlantic Salmon Gill Is Post-Transcriptionally Regulated by miRNAs in the Parr-Smolt Transformation and Adaptation to Sea Water. Int. J. Mol. Sci. 2022, 23, 8831. [Google Scholar] [CrossRef] [PubMed]
  16. Juanchich, A.; Bardou, P.; Rué, O.; Gabillard, J.-C.; Gaspin, C.; Bobe, J.; Guiguen, Y. Characterization of an extensive rainbow trout miRNA transcriptome by next generation sequencing. BMC Genom. 2016, 17, 164. [Google Scholar] [CrossRef]
  17. Cardona, E.; Guyomar, C.; Desvignes, T.; Montfort, J.; Guendouz, S.; Postlethwait, J.H.; Skiba-Cassy, S.; Bobe, J. Circulating miRNA repertoire as a biomarker of metabolic and reproductive states in rainbow trout. BMC Biol. 2021, 19, 235. [Google Scholar] [CrossRef]
  18. Huang, C.W.; Li, Y.H.; Hu, S.Y.; Chi, J.R.; Lin, G.H.; Lin, C.C.; Gong, H.Y.; Chen, J.Y.; Chen, R.H.; Chang, S.J.; et al. Differential expression patterns of growth-related microRNAs in the skeletal muscle of Nile tilapia (Oreochromis niloticus). J. Anim. Sci. 2012, 90, 4266–4279. [Google Scholar] [CrossRef]
  19. Herkenhoff, M.E.; Bovolenta, L.A.; Broedel, O.; Dos Santos, L.D.; de Oliveira, A.C.; Chuffa, L.G.A.; Ribeiro, A.d.O.; Lupi, L.A.; Dias, M.A.D.; Hilsdorf, A.W.S.; et al. Variant expression signatures of microRNAs and protein related to growth in a crossbreed between two strains of Nile tilapia (Oreochromis niloticus). Genomics 2021, 113, 4303–4312. [Google Scholar] [CrossRef]
  20. Eslamloo, K.; Inkpen, S.M.; Rise, M.L.; Andreassen, R. Discovery of microRNAs associated with the antiviral immune response of Atlantic cod macrophages. Mol. Immunol. 2018, 93, 152–161. [Google Scholar] [CrossRef]
  21. Andreassen, R.; Rangnes, F.; Sivertsen, M.; Chiang, M.; Tran, M.; Worren, M.M. Discovery of miRNAs and Their Corresponding miRNA Genes in Atlantic Cod (Gadus morhua): Use of Stable miRNAs as Reference Genes Reveals Subgroups of miRNAs That Are Highly Expressed in Particular Organs. PLoS ONE 2016, 11, e0153324. [Google Scholar] [CrossRef] [PubMed]
  22. Chakraborty, S.; Woldemariam, N.T.; Visnovska, T.; Rise, M.L.; Boyce, D.; Santander, J.; Andreassen, R. Characterization of miRNAs in Embryonic, Larval, and Adult Lumpfish Provides a Reference miRNAome for Cyclopterus lumpus. Biology 2022, 11, 130. [Google Scholar] [CrossRef] [PubMed]
  23. Xu, Z.; Chen, J.; Li, X.; Ge, J.; Pan, J.; Xu, X. Identification and Characterization of MicroRNAs in Channel Catfish (Ictalurus punctatus) by Using Solexa Sequencing Technology. PLoS ONE 2013, 8, e54174. [Google Scholar] [CrossRef] [PubMed]
  24. Tang, X.; Fu, J.; Tan, X.; Shi, Y.; Ye, J.; Guan, W.; Shi, Y.; Xu, M. The miR-155 regulates cytokines expression by SOSC1 signal pathways of fish in vitro and in vivo. Fish Shellfish. Immunol. 2020, 106, 28–35. [Google Scholar] [CrossRef] [PubMed]
  25. Desvignes, T.; Bardou, P.; Montfort, J.; Sydes, J.; Guyomar, C.; George, S.; Postlethwait, J.H.; Bobe, J. FishmiRNA: An Evolutionarily Supported MicroRNA Annotation and Expression Database for Ray-Finned Fishes. Mol. Biol. Evol. 2022, 39, msac004. [Google Scholar] [CrossRef] [PubMed]
  26. Woldemariam, N.T.; Agafonov, O.; Sindre, H.; Høyheim, B.; Houston, R.D.; Robledo, D.; Bron, J.E.; Andreassen, R. miRNAs Predicted to Regulate Host Anti-viral Gene Pathways in IPNV-Challenged Atlantic Salmon Fry Are Affected by Viral Load, and Associated with the Major IPN Resistance QTL Genotypes in Late Infection. Front. Immunol. 2020, 11, 2113. [Google Scholar] [CrossRef] [PubMed]
  27. Ramberg, S.; Krasnov, A.; Colquhoun, D.; Wallace, C.; Andreassen, R. Expression Analysis of Moritella viscosa-Challenged Atlantic Salmon Identifies Disease-Responding Genes, MicroRNAs and Their Predicted Target Genes and Pathways. Int. J. Mol. Sci. 2022, 23, 11200. [Google Scholar] [CrossRef] [PubMed]
  28. Wienholds, E.; Plasterk, R.H.A. MicroRNA function in animal development. FEBS Lett. 2005, 579, 5911–5922. [Google Scholar] [CrossRef]
  29. Su, H.; Fan, J.; Ma, D.; Zhu, H. Identification and Characterization of Osmoregulation Related MicroRNAs in Gills of Hybrid Tilapia under Three Types of Osmotic Stress. Front. Genet. 2021, 12, 526277. [Google Scholar] [CrossRef]
  30. Huang, C.-H.; Chen, N.; Huang, C.-X.; Zhang, B.; Wu, M.; He, L.; Liu, H.; Tang, R.; Wang, W.-M.; Wang, H.-L. Involvement of the miR-462/731 cluster in hypoxia response in Megalobrama amblycephala. Fish Physiol. Biochem. 2017, 43, 863–873. [Google Scholar] [CrossRef]
  31. Østbye, T.-K.K.; Woldemariam, N.T.; Lundberg, C.E.; Berge, G.M.; Ruyter, B.; Andreassen, R. Modulation of hepatic miRNA expression in Atlantic salmon (Salmo salar) by family background and dietary fatty acid composition. J. Fish Biol. 2021, 98, 1172–1185. [Google Scholar] [CrossRef] [PubMed]
  32. FAO. In Brief to the State of World Fisheries and Aquaculture 2022; FAO: Rome, Italy, 2022. [Google Scholar]
  33. FAO. The State of World Fisheries and Aquaculture 2020. Sustainability in Action; FAO: Rome, Italy, 2020. [Google Scholar]
  34. Jankowska, B.; Zakes, Z.; Zmijewski, T.; Szczepkowski, M. A comparison of selected quality features of the tissue and slaughter yield of wild and cultivated pikeperch Sander lucioperca (L.). Eur. Food Res. Technol. 2003, 217, 401–405. [Google Scholar] [CrossRef]
  35. Kestemont, P.; Dabrowski, K.; Summerfelt, R.C. Biology and Culture of Percid Fishes: Principles and Practices; Springer: Dordrecht, The Netherlands, 2015; ISBN 9789401772273. [Google Scholar]
  36. Zakęś, Z. FAO 2023 Cultured Aquatic Species Information Programme. Sander lucioperca; FAO: Rome, Italy, 2012. [Google Scholar]
  37. Fontaine, P.; Kestemont, P.; Teletchea, F.; Wang, N. Percid Fish Culture. From Research to Production; Presses Universitaires de Namur: Namur, Belgium, 2008; ISBN 9782870375822. [Google Scholar]
  38. Policar, T.; Schaefer, F.J.; Panana, E.; Meyer, S.; Teerlinck, S.; Toner, D.; Żarski, D. Recent progress in European percid fish culture production technology—Tackling bottlenecks. Aquac. Int. 2019, 27, 1151–1174. [Google Scholar] [CrossRef]
  39. Steenfeldt, S. Chapter 10: Culture Methods of Pikeperch Early Life Stages. In Biology and Culture of Percid Fishes; Springer: Dordrecht, The Netherlands, 2015; pp. 295–312. ISBN 9789401772273. [Google Scholar]
  40. Liu, Y.; Rosten, T.W.; Henriksen, K.; Hognes, E.S.; Summerfelt, S.; Vinci, B. Comparative economic performance and carbon footprint of two farming models for producing Atlantic salmon (Salmo salar): Land-based closed containment system in freshwater and open net pen in seawater. Aquac. Eng. 2016, 71, 1–12. [Google Scholar] [CrossRef]
  41. Bregnballe, J.; Eurofish; FAO. A Guide to Recirculation Aquaculture: An Introduction to the New Environmentally Friendly and Highly Productive Closed Fish Farming Systems; Food and Agriculture Organization of the United Nations: Rome, Italy, 2022; ISBN 9789251369814. [Google Scholar]
  42. Perry, W.B.; Solberg, M.F.; Besnier, F.; Dyrhovden, L.; Matre, I.H.; Fjelldal, P.G.; Ayllon, F.; Creer, S.; Llewellyn, M.; Taylor, M.I.; et al. Evolutionary drivers of kype size in Atlantic salmon (Salmo salar): Domestication, age and genetics. R. Soc. Open Sci. 2019, 6, 190021. [Google Scholar] [CrossRef] [PubMed]
  43. Lehnert, S.J.; Kess, T.; Bentzen, P.; Kent, M.P.; Lien, S.; Gilbey, J.; Clément, M.; Jeffery, N.W.; Waples, R.S.; Bradbury, I.R. Genomic signatures and correlates of widespread population declines in salmon. Nat. Commun. 2019, 10, 2996. [Google Scholar] [CrossRef] [PubMed]
  44. García-Pérez, I.; Molsosa-Solanas, A.; Perelló-Amorós, M.; Sarropoulou, E.; Blasco, J.; Gutiérrez, J.; Garcia de la Serrana, D. The Emerging Role of Long Non-Coding RNAs in Development and Function of Gilthead Sea Bream (Sparus aurata) Fast Skeletal Muscle. Cells 2022, 11, 428. [Google Scholar] [CrossRef] [PubMed]
  45. Nguinkal, J.A.; Brunner, R.M.; Verleih, M.; Rebl, A.; de Los Ríos-Pérez, L.; Schäfer, N.; Hadlich, F.; Stüeken, M.; Wittenburg, D.; Goldammer, T. The First Highly Contiguous Genome Assembly of Pikeperch (Sander lucioperca), an Emerging Aquaculture Species in Europe. Genes 2019, 10, 708. [Google Scholar] [CrossRef]
  46. de Los Ríos-Pérez, L.; Nguinkal, J.A.; Verleih, M.; Rebl, A.; Brunner, R.M.; Klosa, J.; Schäfer, N.; Stüeken, M.; Goldammer, T.; Wittenburg, D. An ultra-high density SNP-based linkage map for enhancing the pikeperch (Sander lucioperca) genome assembly to chromosome-scale. Sci. Rep. 2020, 10, 22335. [Google Scholar] [CrossRef]
  47. Nguinkal, J.A.; Verleih, M.; De Los Ríos-Pérez, L.; Brunner, R.M.; Sahm, A.; Bej, S.; Rebl, A.; Goldammer, T. Comprehensive Characterization of Multitissue Expression Landscape, Co-Expression Networks and Positive Selection in Pikeperch. Cells 2021, 10, 2289. [Google Scholar] [CrossRef]
  48. Altuvia, Y.; Landgraf, P.; Lithwick, G.; Elefant, N.; Pfeffer, S.; Aravin, A.; Brownstein, M.J.; Tuschl, T.; Margalit, H. Clustering and conservation patterns of human microRNAs. Nucleic Acids Res. 2005, 33, 2697–2706. [Google Scholar] [CrossRef] [PubMed]
  49. Li, S.-C.; Chan, W.-C.; Hu, L.-Y.; Lai, C.-H.; Hsu, C.-N.; Lin, W. Identification of homologous microRNAs in 56 animal genomes. Genom. 2010, 96, 1–9. [Google Scholar] [CrossRef] [PubMed]
  50. Shabalina, S.A.; Koonin, E.V. Origins and evolution of eukaryotic RNA interference. Trends Ecol. Evol. 2008, 23, 578–587. [Google Scholar] [CrossRef] [PubMed]
  51. Lu, J.; Shen, Y.; Wu, Q.; Kumar, S.; He, B.; Shi, S.; Carthew, R.W.; Wang, S.M.; Wu, C.-I. The birth and death of microRNA genes in Drosophila. Nat. Genet. 2008, 40, 351–355. [Google Scholar] [CrossRef] [PubMed]
  52. Baskerville, S.; Bartel, D.P. Microarray profiling of microRNAs reveals frequent coexpression with neighboring miRNAs and host genes. RNA 2005, 11, 241–247. [Google Scholar] [CrossRef] [PubMed]
  53. Lagos-quintana, M.; Rauhut, R.; Lendeckel, W.; Tuschl, T. Identification of novel genes coding for small expressed RNAs. Science 2001, 294, 853–858. [Google Scholar] [CrossRef] [PubMed]
  54. Kozomara, A.; Griffiths-Jones, S. miRBase: Annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 2014, 42, D68–D73. [Google Scholar] [CrossRef]
  55. Giraldez, A.J.; Mishima, Y.; Rihel, J.; Grocock, R.J.; Van Dongen, S.; Inoue, K.; Enright, A.J.; Schier, A.F. Zebrafish MiR-430 promotes deadenylation and clearance of maternal mRNAs. Science 2006, 312, 75–79. [Google Scholar] [CrossRef]
  56. Giusti, J.; Pinhal, D.; Moxon, S.; Campos, C.L.; Münsterberg, A.; Martins, C. MicroRNA-10 modulates Hox genes expression during Nile tilapia embryonic development. Mech. Dev. 2016, 140, 12–18. [Google Scholar] [CrossRef]
  57. Tehler, D.; Høyland-Kroghsbo, N.M.; Lund, A.H. The miR-10 microRNA precursor family. RNA Biol. 2011, 8, 728–734. [Google Scholar] [CrossRef]
  58. Radhakrishnan, B.; Alwin Prem Anand, A. Role of miRNA-9 in Brain Development. J. Exp. Neurosci. 2016, 10, 101–120. [Google Scholar] [CrossRef] [PubMed]
  59. Tessmar-Raible, K.; Raible, F.; Christodoulou, F.; Guy, K.; Rembold, M.; Hausen, H.; Arendt, D. Conserved sensory-neurosecretory cell types in annelid and fish forebrain: Insights into hypothalamus evolution. Cell 2007, 129, 1389–1400. [Google Scholar] [CrossRef] [PubMed]
  60. Giraldez, A.J.; Cinalli, R.M.; Glasner, M.E.; Enright, A.J.; Thomson, J.M.; Baskerville, S.; Hammond, S.M.; Bartel, D.P.; Schier, A.F. MicroRNAs regulate brain morphogenesis in zebrafish. Science 2005, 308, 833–838. [Google Scholar] [CrossRef] [PubMed]
  61. Li, X.; Cassidy, J.J.; Reinke, C.A.; Fischboeck, S.; Carthew, R.W. A microRNA imparts robustness against environmental fluctuation during development. Cell 2009, 137, 273–282. [Google Scholar] [CrossRef] [PubMed]
  62. Mitchelson, K.R.; Qin, W.-Y. Roles of the canonical myomiRs miR-1, -133 and -206 in cell development and disease. World J. Biol. Chem. 2015, 6, 162–208. [Google Scholar] [CrossRef] [PubMed]
  63. Townley-Tilson, W.H.D.; Callis, T.E.; Wang, D. MicroRNAs 1, 133, and 206: Critical factors of skeletal and cardiac muscle development, function, and disease. Int. J. Biochem. Cell Biol. 2010, 42, 1252–1255. [Google Scholar] [CrossRef]
  64. Esau, C.; Kang, X.; Peralta, E.; Hanson, E.; Marcusson, E.G.; Ravichandran, L.V.; Sun, Y.; Koo, S.; Perera, R.J.; Jain, R.; et al. MicroRNA-143 Regulates Adipocyte Differentiation. J. Biol. Chem. 2004, 279, 52361–52365. [Google Scholar] [CrossRef]
  65. Mennigen, J.A.; Skiba-Cassy, S.; Panserat, S. Ontogenesis of expression of metabolic genes and microRNAs in rainbow trout alevins during the transition from the endogenous to the exogenous feeding period. J. Exp. Biol. 2013, 216, 1597–1608. [Google Scholar] [CrossRef]
  66. Rasmussen, K.D.; Simmini, S.; Abreu-Goodger, C.; Bartonicek, N.; Di Giacomo, M.; Bilbao-Cortes, D.; Horos, R.; Von Lindern, M.; Enright, A.J.; O’carroll, D. The miR-144/451 locus is required for erythroid homeostasis. J. Exp. Med. 2010, 207, 1351–1358. [Google Scholar] [CrossRef]
  67. Zhao, Y.; Fan, S.; Yuan, P.; Li, G. Expression characteristics and interaction networks of microRNAs in spleen tissues of grass carp (Ctenopharyngodon idella). PLoS ONE 2022, 17, e0266189. [Google Scholar] [CrossRef]
  68. Cui, J.; Gao, Y.; Chu, Q.; Bi, D.; Xu, T. miRNA-8159 is involved in TLR signaling pathway regulation after pathogen infection by direct targeting TLR13 in miiuy croaker. Fish Shellfish. Immunol. 2017, 66, 531–539. [Google Scholar] [CrossRef] [PubMed]
  69. Jopling, C. Liver-specific microRNA-122: Biogenesis and function. RNA Biol. 2012, 9, 137–142. [Google Scholar] [CrossRef] [PubMed]
  70. Ren, F.-J.; Yao, Y.; Cai, X.-Y.; Fang, G.-Y. Emerging Role of MiR-192-5p in Human Diseases. Front. Pharmacol. 2021, 12, 614068. [Google Scholar] [CrossRef] [PubMed]
  71. Ng, H.M.; Ho, J.C.H.; Nong, W.; Hui, J.H.L.; Lai, K.P.; Wong, C.K.C. Genome-wide analysis of MicroRNA-messenger RNA interactome in ex-vivo gill filaments, Anguilla japonica. BMC Genom. 2020, 21, 208. [Google Scholar] [CrossRef] [PubMed]
  72. Johansen, I.; Andreassen, R. Validation of miRNA genes suitable as reference genes in qPCR analyses of miRNA gene expression in Atlantic salmon (Salmo salar). BMC Res. Notes 2014, 8, 945. [Google Scholar] [CrossRef]
  73. Kunz, Y.W. Developmental Biology of Teleost Fishes; Springer: Dordrecht, The Netherlands, 2004; ISBN 9781402029974. [Google Scholar]
  74. Alberti, C.; Cochella, L. A framework for understanding the roles of miRNAs in animal development. Development 2017, 144, 2548–2559. [Google Scholar] [CrossRef]
  75. Foshay, K.M.; Gallicano, G.I. miR-17 family miRNAs are expressed during early mammalian development and regulate stem cell differentiation. Dev. Biol. 2009, 326, 431–443. [Google Scholar] [CrossRef]
  76. Singh, S.; Saini, H.; Sharma, A.; Gupta, S.; Huddar, V.G.; Tripathi, R. Breast cancer: miRNAs monitoring chemoresistance and systemic therapy. Front. Oncol. 2023, 13, 1155254. [Google Scholar] [CrossRef]
  77. Concepcion, C.P.; Bonetti, C.; Ventura, A. The microRNA-17-92 family of microRNA clusters in development and disease. Cancer J. 2012, 18, 262–267. [Google Scholar] [CrossRef] [PubMed]
  78. Ventura, A.; Young, A.G.; Winslow, M.M.; Lintault, L.; Meissner, A.; Erkeland, S.J.; Newman, J.; Bronson, R.T.; Crowley, D.; Stone, J.R.; et al. Targeted deletion reveals essential and overlapping functions of the miR-17 through 92 family of miRNA clusters. Cell 2008, 132, 875–886. [Google Scholar] [CrossRef] [PubMed]
  79. Reinhart, B.J.; Slack, F.J.; Basson, M.; Pasquinelli, A.E.; Bettinger, J.C.; Rougvie, A.E.; Horvitz, H.R.; Ruvkun, G. The 21-nucleotide let-7 RNA regulates developmental timing in Caenorhabditis elegans. Nature 2000, 403, 901–906. [Google Scholar] [CrossRef] [PubMed]
  80. Abbott, A.L.; Alvarez-Saavedra, E.; Miska, E.A.; Lau, N.C.; Bartel, D.P.; Horvitz, H.R.; Ambros, V. The let-7 MicroRNA family members mir-48, mir-84, and mir-241 function together to regulate developmental timing in Caenorhabditis elegans. Dev. Cell 2005, 9, 403–414. [Google Scholar] [CrossRef] [PubMed]
  81. Horak, M.; Novak, J.; Bienertova-Vasku, J. Muscle-specific microRNAs in skeletal muscle development. Dev. Biol. 2016, 410, 1–13. [Google Scholar] [CrossRef] [PubMed]
  82. Desvignes, T.; Batzel, P.; Sydes, J.; Eames, B.F.; Postlethwait, J.H. miRNA analysis with Prost! reveals evolutionary conservation of organ-enriched expression and post-transcriptional modifications in three-spined stickleback and zebrafish. Sci. Rep. 2019, 9, 3913. [Google Scholar] [CrossRef] [PubMed]
  83. Campos, C.; Sundaram, A.Y.; Valente, L.M.; Conceição, L.E.; Engrola, S.; Fernandes, J.M. Thermal plasticity of the miRNA transcriptome during Senegalese sole development. BMC Genom. 2014, 15, 525. [Google Scholar] [CrossRef] [PubMed]
  84. Richbourg, H.A.; Hu, D.P.; Xu, Y.; Barczak, A.J.; Marcucio, R.S. miR-199 family contributes to regulation of sonic hedgehog expression during craniofacial development. Dev. Dyn. 2020, 249, 1062–1076. [Google Scholar] [CrossRef] [PubMed]
  85. Desvignes, T.; Contreras, A.; Postlethwait, J.H. Evolution of the miR199-214 cluster and vertebrate skeletal development. RNA Biol. 2014, 11, 281–294. [Google Scholar] [CrossRef]
  86. Cordes, K.R.; Sheehy, N.T.; White, M.P.; Berry, E.C.; Morton, S.U.; Muth, A.N.; Lee, T.-H.; Miano, J.M.; Ivey, K.N.; Srivastava, D. miR-145 and miR-143 regulate smooth muscle cell fate and plasticity. Nature 2009, 460, 705–710. [Google Scholar] [CrossRef]
  87. Vacante, F.; Denby, L.; Sluimer, J.C.; Baker, A.H. The function of miR-143, miR-145 and the MiR-143 host gene in cardiovascular development and disease. Vascul. Pharmacol. 2019, 112, 24–30. [Google Scholar] [CrossRef]
  88. Cao, Y.; Wang, D.; Li, S.; Xu, L.; Zhao, J.; Liu, H.; Lu, T.; Zhang, Q. Identification and analysis of differentially expressed microRNAs in rainbow trout (Oncorhynchus mykiss) responding to infectious hematopoietic necrosis virus infection. Dev. Comp. Immunol. 2018, 88, 28–36. [Google Scholar] [CrossRef]
  89. Ordas, A.; Kanwal, Z.; Lindenberg, V.; Rougeot, J.; Mink, M.; Spaink, H.P.; Meijer, A.H. MicroRNA-146 function in the innate immune transcriptome response of zebrafish embryos to Salmonella typhimurium infection. BMC Genom. 2013, 14, 696. [Google Scholar] [CrossRef] [PubMed]
  90. Bjørgen, H.; Koppang, E.O. Anatomy of teleost fish immune structures and organs. Immunogenetics 2021, 73, 53–63. [Google Scholar] [CrossRef] [PubMed]
  91. Tang, X.; Fu, J.; Yao, Y.; Xu, M. Identification and characterization of immune-related microRNAs in hybrid snakehead(Channa maculata♀ × Channa argus♂)after treated by Echinacea purpurea (Linn.) Moench. Fish Shellfish. Immunol. 2023, 135, 108653. [Google Scholar] [CrossRef] [PubMed]
  92. Li, G.; Zhao, Y.; Wen, L.; Liu, Z.; Yan, F.; Gao, C. Identification and characterization of MicroRNAs in the spleen of common carp immune organ. J. Cell. Biochem. 2014, 115, 1768–1778. [Google Scholar] [CrossRef] [PubMed]
  93. Zuo, H.; Liu, X.; Luo, M.; Yang, L.; Zhu, Z.; Weng, S.; He, J.; Xu, X. miR-10c Facilitates White Spot Syndrome Virus Infection by Targeting Toll3 in Litopenaeus vannemei. Front. Immunol. 2021, 12, 733730. [Google Scholar] [CrossRef] [PubMed]
  94. Tahamtan, A.; Teymoori-Rad, M.; Nakstad, B.; Salimi, V. Anti-Inflammatory MicroRNAs and Their Potential for Inflammatory Diseases Treatment. Front. Immunol. 2018, 9, 1377. [Google Scholar] [CrossRef]
  95. Rasal, K.D.; Iquebal, M.A.; Pandey, A.; Behera, P.; Jaiswal, S.; Vasam, M.; Dixit, S.; Raza, M.; Sahoo, L.; Nandi, S.; et al. Revealing liver specific microRNAs linked with carbohydrate metabolism of farmed carp, Labeo rohita (Hamilton, 1822). Genomics 2020, 112, 32–44. [Google Scholar] [CrossRef]
  96. Mennigen, J.A.; Martyniuk, C.J.; Seiliez, I.; Panserat, S.; Skiba-Cassy, S. Metabolic consequences of microRNA-122 inhibition in rainbow trout, Oncorhynchus mykiss. BMC Genom. 2014, 15, 70. [Google Scholar] [CrossRef]
  97. Huang, C.; Chen, N.; Wu, X.; Huang, C.; He, Y.; Tang, R.; Wang, W.; Wang, H. The zebrafish miR-462/miR-731 cluster is induced under hypoxic stress via hypoxia-inducible factor 1α and functions in cellular adaptations. FASEB J. 2015, 29, 4901–4913. [Google Scholar] [CrossRef]
  98. Schyth, B.D.; Bela-Ong, D.B.; Jalali, S.A.H.; Kristensen, L.B.J.; Einer-Jensen, K.; Pedersen, F.S.; Lorenzen, N. Two Virus-Induced MicroRNAs Known Only from Teleost Fishes Are Orthologues of MicroRNAs Involved in Cell Cycle Control in Humans. PLoS ONE 2015, 10, e0132434. [Google Scholar] [CrossRef]
  99. Schäfer, N.; Matoušek, J.; Rebl, A.; Stejskal, V.; Brunner, R.M.; Goldammer, T.; Verleih, M.; Korytář, T. Effects of Chronic Hypoxia on the Immune Status of Pikeperch (Sander lucioperca Linnaeus, 1758). Biology 2021, 10, 649. [Google Scholar] [CrossRef] [PubMed]
  100. Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 2011, 17, 10. [Google Scholar] [CrossRef]
  101. Friedländer, M.R.; Mackowiak, S.D.; Li, N.; Chen, W.; Rajewsky, N. miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 2012, 40, 37–52. [Google Scholar] [CrossRef] [PubMed]
  102. Ambros, V.; Bartel, B.; Bartel, D.P.; Burge, C.B.; Carrington, J.C.; Chen, X.; Dreyfuss, G.; Eddy, S.R.; Griffiths-Jones, S.; Marshall, M.; et al. A uniform system for microRNA annotation. RNA 2003, 9, 277–279. [Google Scholar] [CrossRef] [PubMed]
  103. Medley, J.C.; Panzade, G.; Zinovyeva, A.Y. microRNA strand selection: Unwinding the rules. Wiley Interdiscip. Rev. RNA 2021, 12, e1627. [Google Scholar] [CrossRef] [PubMed]
  104. Xie, F.; Xiao, P.; Chen, D.; Xu, L.; Zhang, B. miRDeepFinder: A miRNA analysis tool for deep sequencing of plant small RNAs. Plant Mol. Biol. 2012, 80, 75–84. [Google Scholar] [CrossRef] [PubMed]
  105. Vandesompele, J.; De Preter, K.; Pattyn, F.; Poppe, B.; Van Roy, N.; De Paepe, A.; Speleman, F. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3, research0034.1. [Google Scholar] [CrossRef] [PubMed]
  106. Andersen, C.L.; Jensen, J.L.; Ørntoft, T.F. Normalization of Real-Time Quantitative Reverse Transcription-PCR Data: A Model-Based Variance Estimation Approach to Identify Genes Suited for Normalization, Applied to Bladder and Colon Cancer Data Sets. Cancer Res. 2004, 64, 5245–5250. [Google Scholar] [CrossRef]
  107. Pfaffl, M.W.; Tichopad, A.; Prgomet, C.; Neuvians, T.P. Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper-Excel-based tool using pair-wise correlations. Biotechnol. Lett. 2004, 26, 509–515. [Google Scholar] [CrossRef]
  108. Silver, N.; Best, S.; Jiang, J.; Thein, S.L. Selection of housekeeping genes for gene expression studies in human reticulocytes using real-time PCR. BMC Mol. Biol. 2006, 7, 33. [Google Scholar] [CrossRef]
  109. Hellemans, J.; Mortier, G.; De Paepe, A.; Speleman, F.; Vandesompele, J. qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biol. 2007, 8, R19. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Top 20 most abundant mature miRNAs in pikeperch. The boxes show the scaled proportion (%) of mature miRNAs in two respective replicates of eyed eggs and three larvae stages (two days post hatching [dph], 15 dph, and 26 dph). Each mature miRNA is labeled in a specific color given in the legend.
Figure 1. Top 20 most abundant mature miRNAs in pikeperch. The boxes show the scaled proportion (%) of mature miRNAs in two respective replicates of eyed eggs and three larvae stages (two days post hatching [dph], 15 dph, and 26 dph). Each mature miRNA is labeled in a specific color given in the legend.
Ijms 25 00189 g001
Figure 2. Top 20 most abundant mature miRNAs in pikeperch. The boxes show the scaled proportion (%) of mature miRNAs in two respective replicates of spleen tissues from two fingerling stages (3 g and 35 g) and adult fish. Each mature miRNA is labeled in a specific color given in the legend.
Figure 2. Top 20 most abundant mature miRNAs in pikeperch. The boxes show the scaled proportion (%) of mature miRNAs in two respective replicates of spleen tissues from two fingerling stages (3 g and 35 g) and adult fish. Each mature miRNA is labeled in a specific color given in the legend.
Ijms 25 00189 g002
Figure 3. Top 20 most abundant mature miRNAs in pikeperch. The boxes show the scaled proportion (%) of mature miRNAs in two respective replicates of head kidney and gill, as well as of one adult fish liver sample. Each mature miRNA is labeled in a specific color given in the legend.
Figure 3. Top 20 most abundant mature miRNAs in pikeperch. The boxes show the scaled proportion (%) of mature miRNAs in two respective replicates of head kidney and gill, as well as of one adult fish liver sample. Each mature miRNA is labeled in a specific color given in the legend.
Ijms 25 00189 g003
Figure 4. The validation of miRNA expression data by qPCR. (a) The table lists ten selected affected mature miRNAs and the compared samples. * developmental stage/organ in which the miRNA is predicted to be significantly enriched (adjusted p-value of ≤0.05, Log2FC ≥ 1). (b) Heatmap displaying differential expression (log10 fold change) based on norm. reads/DESeq2 and qPCR data. (c) Pearson correlation between small RNA-seq (x-axis) and qPCR (y-axis) data for the selected miRNAs. Symbol shape and colors correspond to the respective miRNAs listed beside the plot.
Figure 4. The validation of miRNA expression data by qPCR. (a) The table lists ten selected affected mature miRNAs and the compared samples. * developmental stage/organ in which the miRNA is predicted to be significantly enriched (adjusted p-value of ≤0.05, Log2FC ≥ 1). (b) Heatmap displaying differential expression (log10 fold change) based on norm. reads/DESeq2 and qPCR data. (c) Pearson correlation between small RNA-seq (x-axis) and qPCR (y-axis) data for the selected miRNAs. Symbol shape and colors correspond to the respective miRNAs listed beside the plot.
Ijms 25 00189 g004
Figure 5. Expression of slu-mir-731-5p and slu-mir-462-5p in response to hypoxic stress (evaluated by qPCR). The columns represent the n-fold expression (means + SEM, n = 4) under controlled (light blue) and hypoxic (orange) conditions, normalized to two reference miRNAs. Significant differences in expression levels between the conditions are indicated by asterisks (Student’s t-test; * p < 0.05).
Figure 5. Expression of slu-mir-731-5p and slu-mir-462-5p in response to hypoxic stress (evaluated by qPCR). The columns represent the n-fold expression (means + SEM, n = 4) under controlled (light blue) and hypoxic (orange) conditions, normalized to two reference miRNAs. Significant differences in expression levels between the conditions are indicated by asterisks (Student’s t-test; * p < 0.05).
Ijms 25 00189 g005
Table 1. Distribution of unique mature miRNAs in pikeperch by developmental stage and organ.
Table 1. Distribution of unique mature miRNAs in pikeperch by developmental stage and organ.
Sample GroupmiRs (Total)
Eyed eggs374
Larvae 2 dph375
Larvae 15 dph372
Larvae 26 dph371
Spleen Fingerling 3 g363
Spleen Fingerling 35 g350
Head Kidney368
Spleen366
Liver325
Gills357
Table 2. Top ten differentially expressed mature (assumed guide) miRNAs in pikeperch based on comparative DESeq2 expression analysis.
Table 2. Top ten differentially expressed mature (assumed guide) miRNAs in pikeperch based on comparative DESeq2 expression analysis.
miRNAFold ChangeAdj. p-Value
Eyed eggs and larvae 2 dph vs. larvae 15 dph and 26 dph
  Enriched in eyed eggs and larvae 2 dph:
slu-mir-430b-2-3p268.787.66 × 10−26
slu-mir-430d-1-2-3p235.771.43 × 10−26
slu-mir-430bc-3-3p222.989.14 × 10−27
slu-mir-430bc-2-3p181.614.02 × 10−24
slu-mir-430c-1-3-5p73.144.10 × 10−18
slu-mir-430b-1-4-5p64.959.95 × 10−15
slu-mir-430a-5p59.582.55 × 10−15
slu-mir-nov2ab-5p57.851.68 × 10−15
slu-mir-nov2ab-3p21.871.85 × 10−15
slu-mir-219a-2-3p20.555.57 × 10−18
  Enriched in larvae 15 and 26 dph:
slu-let-7d-1-5p20.049.67 × 10−17
slu-let-7f-1-5p17.991.02 × 10−14
slu-let-7h-1-5p17.384.83 × 10−14
slu-mir-726-3p15.094.52 × 10−23
slu-mir-194a-5p12.491.59 × 10−23
slu-mir-31-5p12.051.40 × 10−16
slu-mir-726-5p11.411.37 × 10−14
slu-mir-146-5p9.401.21 × 10−8
slu-mir-194b-5p9.331.67 × 10−13
slu-let-7b-1-2-5p8.593.52 × 10−10
Spleen fingerlings 3 g versus Spleen fingerlings 35 g
  Enriched in fingerling 3 g
slu-mir-nov6-3p50.401.60 × 10−3
slu-mir-183-1-2-5p46.212.02 × 10−5
slu-mir-8159-5p20.853.15 × 10−3
slu-mir-216c-5p19.591.91 × 10−2
slu-mir-216b-5p19.501.91 × 10−2
slu-mir-216a-1-2-5p15.502.50 × 10−2
slu-mir-375-3p15.358.39 × 10−4
slu-mir-122-5p14.253.36 × 10−5
slu-mir-200b-3p5.113.25 × 10−2
slu-mir-182-5p4.793.29 × 10−2
  Enriched in fingerling 35 g
slu-mir-1388-3p6.682.45 × 10−4
slu-mir-730-5p5.325.11 × 10−4
slu-mir-7132-5p5.155.11 × 10−4
slu-mir-731-5p3.234.14 × 10−2
slu-mir-125-1-3-5p3.221.92 × 10−2
slu-mir-99b-5p3.132.50 × 10−2
Head kidney versus all other adult tissues
  Enriched in head kidney
slu-mir-10b-5p165.581.77 × 10−9
slu-mir-551-3p55.043.44 × 10−3
slu-mir-10c-5p6.903.89 × 10−2
Spleen versus all other adult tissues
  Enriched in spleen
slu-mir-202-5p24.592.03 × 10−2
slu-mir-187-3p20.703.52 × 10−2
slu-mir-2187-3p18.556.82 × 10−4
Gills versus all other adult tissues
  Enriched in gills
slu-mir-205-5p585.853.67 × 10−14
slu-mir-203b-3p579.261.08 × 10−18
slu-mir-725-3p447.411.44 × 10−14
slu-mir-203a-3p398.992.42 × 10−40
slu-mir-206-3p144.422.29 × 10−27
slu-mir-205-3p115.311.45 × 10−15
slu-mir-184ab-3p73.513.44 × 10−11
slu-mir-190b-5p42.672.76 × 10−6
slu-mir-375-3p32.952.86 × 10−2
slu-mir-182-5p25.442.51 × 10−17
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Verleih, M.; Visnovska, T.; Nguinkal, J.A.; Rebl, A.; Goldammer, T.; Andreassen, R. The Discovery and Characterization of Conserved and Novel miRNAs in the Different Developmental Stages and Organs of Pikeperch (Sander lucioperca). Int. J. Mol. Sci. 2024, 25, 189. https://doi.org/10.3390/ijms25010189

AMA Style

Verleih M, Visnovska T, Nguinkal JA, Rebl A, Goldammer T, Andreassen R. The Discovery and Characterization of Conserved and Novel miRNAs in the Different Developmental Stages and Organs of Pikeperch (Sander lucioperca). International Journal of Molecular Sciences. 2024; 25(1):189. https://doi.org/10.3390/ijms25010189

Chicago/Turabian Style

Verleih, Marieke, Tina Visnovska, Julien A. Nguinkal, Alexander Rebl, Tom Goldammer, and Rune Andreassen. 2024. "The Discovery and Characterization of Conserved and Novel miRNAs in the Different Developmental Stages and Organs of Pikeperch (Sander lucioperca)" International Journal of Molecular Sciences 25, no. 1: 189. https://doi.org/10.3390/ijms25010189

APA Style

Verleih, M., Visnovska, T., Nguinkal, J. A., Rebl, A., Goldammer, T., & Andreassen, R. (2024). The Discovery and Characterization of Conserved and Novel miRNAs in the Different Developmental Stages and Organs of Pikeperch (Sander lucioperca). International Journal of Molecular Sciences, 25(1), 189. https://doi.org/10.3390/ijms25010189

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop