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

MicroRNA Transcriptomes Reveal Prevalence of Rare and Species-Specific Arm Switching Events During Zebrafish Ontogenesis

In metazoans, microRNAs (miRNAs) are essential regulators of gene expression, affecting critical cellular processes from differentiation and proliferation, to homeostasis. During miRNA biogenesis, the miRNA strand that loads onto the RNA-induced Silencing Complex (RISC) can vary, leading to changes in gene targeting and modulation of biological pathways. To investigate the impact of these "arm switching" events on gene regulation, we analyzed a diverse range of tissues and developmental stages in zebrafish by comparing 5p and 3p arms accumulation dynamics between embryonic developmental stages, adult tissues, and sexes. We also compared variable arm usage patterns observed in zebrafish to other vertebrates including arm switching data from fish, birds, and mammals. Our comprehensive analysis revealed that variable arm usage events predominantly take place during embryonic development. It is also noteworthy that isomiR occurrence correlates to changes in arm selection evidencing an important role of microRNA distinct isoforms in reinforcing and modifying gene regulation by promoting dynamics switches on miRNA 5p and 3p arms accumulation. Our results shed new light on the emergence and coordination of gene expression regulation and pave the way for future investigations in this field.

Small RNA sequencing reveals sex-related miRNAs in Collichthys lucidus

Collichthys lucidus (C. lucidus) is an economically important fish species, exhibiting sexual dimorphism in its growth rate. However, there is a lack of research on its underlying sex-related mechanisms. Therefore, small RNA sequencing was performed to better comprehend these sex-related molecular mechanisms. In total, 171 differentially expressed miRNAs (DE-miRNAs) were identified between the ovaries and testes. Functional enrichment analysis revealed that the target genes of DE-miRNAs were considerably enriched in the p53 signaling, PI3K–Akt signaling, and TGF-beta signaling pathways. In addition, sex-related miRNAs were identified, and the expression of miR-430c-3p and miR-430f-3p was specifically observed in the gonads compared with other organs and their expression was markedly upregulated in the testes relative to the ovaries. Bmp15 was a target of miR-430c-3p and was greatly expressed in the ovaries compared with the testes. Importantly, miR-430c-3p and bmp15 co-expressed in the ovaries and testes. This research provides the first detailed miRNA profiles for C. lucidus concerning sex, likely laying the basis for further studies on sex differentiation in C. lucidus.

The Small RNA Universe of Capitella teleta

RNAi is an evolutionarily fluid mechanism with dramatically different activities across animal phyla. One major group where there has been little investigation is annelid worms. Here, the small RNAs of the polychaete developmental model Capitella teleta are profiled across development. As is seen with nearly all animals, nearly 200 microRNAs were found with 58 high-confidence novel species. Greater miRNA diversity was associated with later stages consistent with differentiation of tissues. Outside miRNA, a distinct composition of other small RNA pathways was found. Unlike many invertebrates, an endogenous siRNA pathway was not observed, indicating pathway loss relative to basal planarians. No processively generated siRNA-class RNAs could be found arising from dsRNA precursors. This has a significant impact on RNAi technology development for this group of animals. Unlike the apparent absence of siRNAs, a significant population of piRNAs was observed. For many piRNAs, phasing and ping-pong biogenesis pathways were identified. Interestingly, piRNAs were found to be highly expressed during early development, suggesting a potential role in regulation in metamorphosis. Critically, the configuration of RNAi factors in C. teleta is found in other annelids and mollusks, suggesting that similar biology is likely to be present in the wider clade. This study is the first in providing comprehensive analysis of small RNAs in annelids.

Dynamic Expression and Gene Regulation of MicroRNAs During Bighead Carp (Hypophthalmichthys nobilis) Early Development

The early development of fish is regulated through dynamic and complex mechanisms involving the regulation of various genes. Many genes are subjected to post-transcriptional regulation by microRNAs (miRNAs). In the Chinese aquaculture industry, the native species bighead carp (Hypophthalmichthys nobilis) is important. However, the genetic regulation related to the early development of bighead carp is unknown. Here, we generated developmental profiles by miRNA sequencing to study the dynamic regulation of miRNAs during bighead carp early development. This study identified 1 046 miRNAs, comprising 312 known miRNAs and 734 uncharacterized miRNAs. Changes in miRNA expression were identified in the six early development stages. An obviously increased expression trend was detected during the development process, with the main burst of activity occurring after the earliest stage (early blastula, DS1). Investigations revealed that several miRNAs were dominantly expressed during the development process, especially in the later stages (e.g., miR-10b-5p, miR-21, miR-92a-3p, miR-206-3p, and miR-430a-3p), suggesting that these miRNAs exerted important functions during embryonic development. The differentially expressed miRNAs (DEMs) and time-serial analysis (profiles) of DEMs were analyzed. A total of 372 miRNAs were identified as DEMs (fold-change >2, and false discovery rate <0.05), and three expression profiles of the DEMs were detected to have co-expression patterns (r > 0.7, and p < 0.05). The broad negative regulation of target genes by miRNAs was speculated, and many development-related biological processes and pathways were enriched for the targets of the DEMs, which might be associated with maternal genome degradation and embryogenesis processes. In conclusion, we revealed the repertoire of miRNAs that are active during early development of bighead carp. These findings will increase our understanding of the regulatory mechanisms of early development of fish.

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

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

MicroRNA roles in regeneration: Multiple lessons from zebrafish

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

The emerging role of microRNAs in fish ovary: A mini review

MicroRNAs (miRNAs) are small non-coding RNAs that regulate gene expression primarily at the post-transcriptional levels. It is now well established that miRNAs are crucial regulators of many developmental and physiological processes, including reproduction. In teleosts, expression profiling studies have shown that miRNAs are expressed in the fish ovary and their levels are regulated during follicle development and by hormones. Using CRISPR/Cas9 mediated gene knockout strategies, several recent studies have provided strong evidence that miR-202 and miR-200s on chromosome 23 play critical roles in regulating ovarian development, oogenesis, and ovulation. In this mini review, we provide a brief overview of canonical miRNA biogenesis and functions; summarize miRNAs that are expressed in fish ovary; and discuss the emerging role of miRNAs in regulating fish ovarian functions.

Single nucleotide polymorphisms affect miRNA target prediction in bovine

Single nucleotide polymorphisms (SNPs) can have significant effects on phenotypic characteristics in cattle. MicroRNAs (miRNAs) are small, non-coding RNAs that act as post-transcriptional regulators by binding them to target mRNAs. In the present study, we scanned ~56 million SNPs against 1,064 bovine miRNA sequences and analyzed, in silico, their possible effects on target binding prediction, primary miRNA formation, association with QTL regions and the evolutionary conservation for each SNP locus. Following target prediction, we show that 71.6% of miRNA predicted targets were altered as a consequence of SNPs located within the seed region of the mature miRNAs. Next, we identified variations in the Minimum Free Energy (MFE), which represents the capacity to alter molecule stability and, consequently, miRNA maturation. A total of 48.6% of the sequences analyzed showed values within those previously reported as sufficient to alter miRNA maturation. We have also found 131 SNPs in 46 miRNAs, with altered target prediction, occurring in QTL regions. Lastly, analysis of evolutionary conservation scores for each SNP locus suggested that they have a conserved biological function through the evolutionary process. Our results suggest that SNPs in microRNAs have the potential to affect bovine phenotypes and could be of great value for genetic improvement studies, as well as production.

Identification of sex differentiation-related microRNA and long non-coding RNA in Takifugu rubripes gonads

Although sex determination and differentiation are key developmental processes in animals, the involvement of non-coding RNA in the regulation of this process is still not clarified. The tiger pufferfish (Takifugu rubripes) is one of the most economically important marine cultured species in Asia, but analyses of miRNA and long non-coding RNA (lncRNA) at early sex differentiation stages have not been conducted yet. In our study, high-throughput sequencing technology was used to sequence transcriptome libraries from undifferentiated gonads of T. rubripes. In total, 231 (107 conserved, and 124 novel) miRNAs were obtained, while 2774 (523 conserved, and 2251 novel) lncRNAs were identified. Of these, several miRNAs and lncRNAs were predicted to be the regulators of the expression of sex-related genes (including fru-miR-15b/foxl2, novel-167, novel-318, and novel-538/dmrt1, novel-548/amh, lnc_000338, lnc_000690, lnc_000370, XLOC_021951, and XR_965485.1/gsdf). Analysis of differentially expressed miRNAs and lncRNAs showed that three mature miRNAs up-regulated and five mature miRNAs were down-regulated in male gonads compared to female gonads, while 79 lncRNAs were up-regulated and 51 were down-regulated. These findings could highlight a group of interesting miRNAs and lncRNAs for future studies and may reveal new insights into the function of miRNAs and lncRNAs in sex determination and differentiation.

A comparative analysis of heart microRNAs in vertebrates brings novel insights into the evolution of genetic regulatory networks

BACKGROUND

During vertebrate evolution, the heart has undergone remarkable changes that lead to morphophysiological differences in the fully formed heart of these species, such as chamber septation, heart rate frequency, blood pressure, and cardiac output volume. Despite these differences, the heart developmental process is guided by a core gene set conserved across vertebrates. Nonetheless, the regulatory mechanisms controlling the expression of genes involved in heart development and maintenance are largely uncharted. MicroRNAs (miRNAs) have been described as important regulatory elements in several biological processes, including heart biology. These small RNA molecules are broadly conserved in sequence and genomic context in metazoans. Mutations may occur in miRNAs and/or genes that contribute to the establishment of distinct repertoires of miRNA-target interactions, thereby favoring the differential control of gene expression and, consequently, the origin of novel phenotypes. In fact, several studies showed that miRNAs are integrated into genetic regulatory networks (GRNs) governing specific developmental programs and diseases. However, studies integrating miRNAs in vertebrate heart GRNs under an evolutionary perspective are still scarce.

RESULTS

We comprehensively examined and compared the heart miRNome of 20 species representatives of the five major vertebrate groups. We found 54 miRNA families with conserved expression and a variable number of miRNA families with group-specific expression in fishes, amphibians, reptiles, birds, and mammals. We also detected that conserved miRNAs present higher expression levels and a higher number of targets, whereas the group-specific miRNAs present lower expression levels and few targets.

CONCLUSIONS

Both the conserved and group-specific miRNAs can be considered modulators orchestrating the core and peripheral genes of heart GRNs of vertebrates, which can be related to the morphophysiological differences and similarities existing in the heart of distinct vertebrate groups. We propose a hypothesis to explain evolutionary differences in the putative functional roles of miRNAs in the heart GRNs analyzed. Furthermore, we present new insights into the molecular mechanisms that could be helping modulate the diversity of morphophysiology in the heart organ of vertebrate species.

Overview of host miRNA properties and their association with epigenetics, long non-coding RNAs, and Xeno-infectious factors

MicroRNA-derived structures play impressive roles in various biological processes. So dysregulation of miRNAs can lead to different human diseases. Recent studies have extended our comprehension of the control of miRNA function and features. Here, we overview some remarkable miRNA properties that have potential implications for the miRNA functions, including different variants of a miRNA called isomiRs, miRNA arm selection/arm switching, and the effect of these factors on miRNA target selection. Besides, we review some aspects of miRNA interactions such as the interaction between epigenetics and miRNA (different miRNAs and their related processing enzymes are epigenetically regulated by multiple DNA methylation enzymes. moreover, DNA methylation could be controlled by diverse mechanisms related to miRNAs), direct and indirect crosstalk between miRNA and lnc (Long Non-Coding) RNAs as a further approach to conduct intercellular regulation called "competing endogenous RNA" (ceRNA) that is involved in the pathogenesis of different diseases, and the interaction of miRNA activities and some Xeno-infectious (virus/bacteria/parasite) factors, which result in modulation of the pathogenesis of infections. This review provides some related studies to a better understanding of miRNA involvement mechanisms and overcoming the complexity of related diseases that may be applicable and useful to prognostic, diagnostic, therapeutic purposes and personalized medicine in the future.

Molecular insight into regulation of miRNAs in the spleen of zebrafish (Danio rerio) upon pathogenic Streptococcus parauberis infection

MicroRNAs (miRNAs) constitute a group of small non-coding RNAs (~22 nucleotides) and one of their main functions is to regulate the immune responses. Gram-positive bacterium, Streptococcus parauberis is the main causative agent of "Streptococcosis" in wide range of fish species. In this study, we performed high throughput sequencing analysis to identify the miRNA profile against S. parauberis infection in the spleen of zebrafish (Danio rerio). Overall, 349 known and 151 novel miRNAs were discovered. Among them, 12 known miRNAs (dre-miR-34b, dre-miR-135a, dre-miR-200b-5p, dre-miR-146b, dre-miR-31, dre-miR-17a-3p, dre-miR-222a-3p, dre-miR-731, dre-miR-301b-3p and dre-miR-30a-3p) and 9 novel miRNAs were differentially expressed (DE) in the spleen of S. parauberis challenged zebrafish. The identified 12 DE miRNAs were predicted to regulate 721 target genes. We confirmed the miRNA expression results by validating selected known and novel DE miRNAs using qRT-PCR. Gene Ontology (GO), Kyoto Encyclopedia of Genes (KEGG) pathway analysis and miRNA-mRNA interactions implies that specific target genes of DE miRNAs are associated with immune responses. The enriched pathways included Toll-like receptor (TLR), C-type lectin, NOD-like receptor, and RIG-I-like receptor signaling pathways, etc. Especially, dre-miR-200b-5p, dre-miR-146b, dre-miR-731, dre-miR-222a-3p, and dre-miR-34b were able to target potential immune-related genes such as il10, irak1, traf6, hspa8 and ikbke upon S. parauberis challenge. Thus, overall results could lay a foundation to understand the underlying immune regulatory role of miRNAs in response to pathogenic S. parauberis infection in teleosts.

Analysis of conserved miRNAs in cynomolgus macaque genome using small RNA sequencing and homology searching

MicroRNAs (miRNAs) are important regulators that fine-tune diverse cellular activities. Cynomolgus macaques (Macaca fascicularis) are used extensively in biomedical and pharmaceutical research; however, substantially fewer miRNAs have been identified in this species than in humans. Consequently, we investigated conserved miRNA profiles in cynomolgus macaques by homology searching and small RNA sequencing. In total, 1,455 high-confidence miRNA gene loci were identified, 408 of which were also confirmed by RNA sequencing, including 73 new miRNA loci reported in cynomolgus macaques for the first time. Comparing miRNA expression with age, we found a positive correlation between sequence conservation and expression levels during miRNA evolution. Additionally, we found that the miRNA gene locations in cynomolgus macaque genome were very flexible. Most were embedded in intergenic spaces or introns and clustered together. Several miRNAs were found in certain gene locations, including 64 exon-resident miRNAs, six splice-site-overlapping miRNAs (SO-miRNAs), and two pairs of distinct mirror miRNAs. We also identified 78 miRNA clusters, 68 of which were conserved in the human genome, including 10 large miRNA clusters predicted to regulate diverse developmental and cellular processes in cynomolgus macaque. Thus, this study not only expands the number of identified miRNAs in cynomolgus macaques but also provides clues for future research on the differences in miRNA repertoire between macaques and humans.

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

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

What's the target: understanding two decades of in silico microRNA-target prediction

MOTIVATION

Since the initial discovery of microRNAs as post-transcriptional, regulatory key players in the 1990s, a total number of $2656$ mature microRNAs have been publicly described for Homo sapiens. As discovery of new miRNAs is still on-going, target identification remains to be an essential and challenging step preceding functional annotation analysis. One key challenge for researchers seems to be the selection of the most appropriate tool out of the larger multiverse of published solutions for a given research study set-up.

RESULTS

In this review we collectively describe the field of in silico target prediction in the course of time and point out long withstanding principles as well as recent developments. By compiling a catalog of characteristics about the 98 prediction methods and identifying common and exclusive traits, we signpost a simplified mechanism to address the problem of application selection. Going further we devised interpretation strategies for common types of output as generated by frequently used computational methods. To this end, our work specifically aims to make prospective users aware of common mistakes and practical questions that arise during the application of target prediction tools.

AVAILABILITY

An interactive implementation of our recommendations including materials shown in the manuscript is freely available at https://www.ccb.uni-saarland.de/mtguide.

Conservation and novelty in the microRNA genomic landscape of hyperdiverse cichlid fishes

MicroRNAs (miRNAs) play crucial roles in the post-transcriptional control of messenger RNA (mRNA). These miRNA-mRNA regulatory networks are present in nearly all organisms and contribute to development, phenotypic divergence, and speciation. To examine the miRNA landscape of cichlid fishes, one of the most species-rich families of vertebrates, we profiled the expression of both miRNA and mRNA in a diverse set of cichlid lineages. Among these, we found that conserved miRNAs differ from recently arisen miRNAs (i.e. lineage specific) in average expression levels, number of target sites, sequence variability, and physical clustering patterns in the genome. Furthermore, conserved miRNA target sites tend to be enriched at the 5' end of protein-coding gene 3' UTRs. Consistent with the presumed regulatory role of miRNAs, we detected more negative correlations between the expression of miRNA-mRNA functional pairs than in random pairings. Finally, we provide evidence that novel miRNA targets sites are enriched in genes involved in protein synthesis pathways. Our results show how conserved and evolutionarily novel miRNAs differ in their contribution to the genomic landscape and highlight their particular evolutionary roles in the adaptive diversification of cichlids.

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

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

Selection of suitable candidate genes for miRNA expression normalization in Yellow River Carp (Cyprinus carpio. var)

Yellow River carp is widely cultivated in the world due to its economic value in aquaculture, and the faster growth of females compared to males. It is believed that microRNAs (miRNA) are involved in gonadal differentiation and development. qPCR is the most preferred method for miRNA functional analysis. Reliable reference genes for normalization in qRT-PCR are the key to ensuring the accuracy of this method. The aim of present research was to evaluate as well as identify the efficacy of reference genes for miRNA expression using qRT-PCR in Yellow River carp. Nine ncRNAs (miR-101, miR-23a, let7a, miR-26a, miR-146a, miR-451, U6, 5S, and 18S) were chosen and tested in four sample sets: (1) different tissues in adult carp, (2) different tissues in juvenile carp, (3) different early developmental stages of carp, and (4) different developmental stages of carp gonads. The stability and suitability values were calculated using NormFinder, geNorm, and BestKeeper software. The results showed that 5S was a suitable reference gene in different tissues of adult and juvenile carp. The genes 5S, 18S, and U6 were the most stable reference genes in the early developmental stages of carp. Let-7a and miR-23a were considered as the suitable reference genes in the development of gonads. All these reference genes were subsequently validated using miR-430. The results showed that genes 5S and 18S were the most suitable reference genes to normalize miRNA expression under normal growth conditions in early different developmental stages. The genes Let-7a, and miR-23a were the most suitable in different developmental stages. The present study is the first comprehensive study of the stability of miRNA reference genes in Yellow River carp, providing valuable as well as basic data for investigating more accurate miRNA expression during gonadal differentiation and development of carp.

Identification and expression profiling analysis of microRNAs in Nile tilapia (Oreochromis niloticus) in response to Streptococcus agalactiae infection

MicroRNAs (miRNAs) play vital regulatory roles in various biological processes, including in immune responses. Nile tilapia (Oreochromis niloticus) is an important commercial fish species in China. To identify immune-related miRNAs of O. niloticus, 4 libraries from liver during S. agalactiae infection (0 h, 5 h, 50 h, and 7 d) were sequenced by high-throughput sequencing technology in tilapia. We obtained 10,703,531, 11,507,163, 11,180,179 and 13,408,414 clean reads per library, respectively. In our results, a total of 482 miRNAs were identified through bioinformatic analysis, including 220 conserved miRNAs and 262 putative novel miRNAs. Moreover, 21 (4.36%), 50 (10.37%), and 46 (9.54%) miRNAs were significantly differentially expressed at 5 h, 50 h and 7 d, respectively. In addition, 6939 target genes regulated by these differentially expressed miRNAs were predicted, and their functional annotations were predicted by Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, which revealed that a majority of differentially expressed miRNAs were involved in apoptotic process, metabolic process, and immune responses. Finally, Real-time quantitative PCR experiments were performed for 7 miRNAs by stem-loop RT-PCR, and a general agreement was confirmed between the sequencing and RT-qPCR data. To our understanding, this is the first report of comprehensive identification of O. niloticus miRNAs being differentially regulated in liver related to S. agalactiae infection. This work provides an opportunity for further understanding of the molecular mechanisms of miRNA regulation in O. niloticus host-pathogen interactions, and genetic resources for molecular assistant selection for disease resistant breeding program.