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

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

Expression Analysis of Moritella viscosa-Challenged Atlantic Salmon Identifies Disease-Responding Genes, MicroRNAs and Their Predicted Target Genes and Pathways.

Moritella viscosa is a bacterial pathogen causing winter-ulcer disease in Atlantic salmon. The lesions on affected fish lead to increased mortality, decreased fish welfare, and inferior meat quality in farmed salmon. MicroRNAs (miRNAs) are small non-coding RNAs involved in post-transcriptional regulation by guiding the miRNA-induced silencing complex to specific mRNA transcripts (target genes). The goal of this study was to identify miRNAs responding to Moritella viscosa in salmon by investigating miRNA expression in the head-kidney and the muscle/skin from lesion sites caused by the pathogen. Protein coding gene expression was investigated by microarray analysis in the same materials. Seventeen differentially expressed guide-miRNAs (gDE-miRNAs) were identified in the head-kidney, and thirty-nine in lesion sites, while the microarray analysis reproduced the differential expression signature of several thousand genes known as infection-responsive. In silico target prediction and enrichment analysis suggested that the gDE-miRNAs were predicted to target genes involved in immune responses, hemostasis, angiogenesis, stress responses, metabolism, cell growth, and apoptosis. The majority of the conserved gDE-miRNAs (e.g., miR-125, miR-132, miR-146, miR-152, miR-155, miR-223 and miR-2188) are known as infection-responsive in other vertebrates. Collectively, the findings indicate that gDE-miRNAs are important post-transcriptional gene regulators of the host response to bacterial infection.

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.

Smoltification (parr-smolt transformation) is a complex developmental process consisting of developmental changes that lead to remodeling of the Atlantic salmon gill. Here, the expression changes of miRNAs and mRNAs were studied by small-RNA sequencing and microarray analysis, respectively, to identify miRNAs and their predicted targets associated with smoltification and subsequent sea water adaptation (SWA). In total, 18 guide miRNAs were identified as differentially expressed (gDE miRNAs). Hierarchical clustering analysis of expression changes divided these into one cluster of 13 gDE miRNAs with decreasing expression during smoltification and SWA that included the miRNA-146, miRNA-30 and miRNA-7132 families. Another smaller cluster that showed increasing expression consisted of miR-101a-3p, miR-193b-5p, miR-499a-5p, miR-727a-3p and miR-8159-5p. The gDE miRNAs were predicted to target 747 of the genes (DE mRNAs), showing expression changes in the microarray analysis. The predicted targets included genes encoding NKA-subunits, aquaporin-subunits, cystic fibrosis transmembrane conductance regulator and the solute carrier family. Furthermore, the predicted target genes were enriched in biological processes associated with smoltification and SWA (e.g., immune system, reactive oxygen species, stress response and extracellular matrix organization). Collectively, the results indicate that remodeling of the gill involves the post-transcriptional regulation of gene expression by the characterized gDE miRNAs.

Expression Analysis in Atlantic Salmon Liver Reveals miRNAs Associated with Smoltification and Seawater Adaptation.

Optimal smoltification is crucial for normal development, growth, and health of farmed Atlantic salmon in seawater. Here, we characterize miRNA expression in liver to reveal whether miRNAs regulate gene expression during this developmental transition. Expression changes of miRNAs and mRNAs was studied by small-RNA sequencing and microarray analysis, respectively. This revealed 62 differentially expressed guide miRNAs (gDE-miRNAs) that could be divided into three groups with characteristic dynamic expression patterns. Three of miRNA families are known as highly expressed in liver. A rare arm shift was observed during smoltification in the Atlantic salmon-specific novel-ssa-miR-16. The gDE-miRNAs were predicted to target 2804 of the genes revealing expression changes in the microarray analysis. Enrichment analysis revealed that targets were significantly enriched in smoltification-associated biological process groups. These included lipid and cholesterol synthesis, carbohydrate metabolism, protein metabolism and protein transport, immune system genes, circadian rhythm and stress response. The results indicate that gDE-miRNAs may regulate many of the changes associated with this developmental transition in liver. The results pave the way for validation of the predicted target genes and further study of gDE-miRNA and their targets by functional assays.

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

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

MicroSalmon: A Comprehensive, Searchable Resource of Predicted MicroRNA Targets and 3'UTR Cis-Regulatory Elements in the Full-Length Sequenced Atlantic Salmon Transcriptome.

Complete 3'UTRs unambiguously assigned to specific mRNA isoforms from the Atlantic salmon full-length (FL) transcriptome were collected into a 3'UTRome. miRNA response elements (MREs) and other cis-regulatory motifs were subsequently predicted and assigned to 3'UTRs of all FL-transcripts. The MicroSalmon GitHub repository provides all results. RNAHybrid and sRNAtoolbox tools predicted the MREs. UTRscan and the Teiresias algorithm predicted other 3'UTR cis-acting motifs, both known vertebrate motifs and putative novel motifs. MicroSalmon provides search programs to retrieve all FL-transcripts targeted by a miRNA (median number 1487), all miRNAs targeting an FL-transcript (median number 27), and other cis-acting motifs. As thousands of FL-transcripts may be targets of each miRNA, additional experimental strategies are necessary to reduce the likely true and relevant targets to a number that may be functionally validated. Low-complexity motifs known to affect mRNA decay in vertebrates were over-represented. Many of these were enriched in the terminal end, while purine- or pyrimidine-rich motifs with unknown functions were enriched immediately downstream of the stop codon. Furthermore, several novel complex motifs were over-represented, indicating conservation and putative function. In conclusion, MicroSalmon is an extensive and useful, searchable resource for study of Atlantic salmon transcript regulation by miRNAs and cis-acting 3'UTR motifs.

Transcriptome Profiling of Atlantic Salmon Adherent Head Kidney Leukocytes Reveals That Macrophages Are Selectively Enriched During Culture.

The Atlantic salmon (Salmo salar) is an economically important fish, both in aquaculture and in the wild. In vertebrates, macrophages are some of the first cell types to respond to pathogen infection and disease. While macrophage biology has been characterized in mammals, less is known in fish. Our previous work identified changes in the morphology, phagocytic ability, and miRNA profile of Atlantic salmon adherent head kidney leukocytes (HKLs) from predominantly "monocyte-like" at Day 1 of in vitro culture to predominantly "macrophage-like" at Day 5 of culture. Therefore, to further characterize these two cell populations, we examined the mRNA transcriptome profile in Day 1 and Day 5 HKLs using a 44K oligonucleotide microarray. Large changes in the transcriptome were revealed, including changes in the expression of macrophage and immune-related transcripts (e.g. csf1r, arg1, tnfa, mx2), lipid-related transcripts (e.g. fasn, dhcr7, fabp6), and transcription factors involved in macrophage differentiation and function (e.g. klf2, klf9, irf7, irf8, stat1). The in silico target prediction analysis of differentially expressed genes (DEGs) using miRNAs known to change expression in Day 5 HKLs, followed by gene pathway enrichment analysis, supported that these miRNAs may be involved in macrophage maturation by targeting specific DEGs. Elucidating how immune cells, such as macrophages, develop and function is a key step in understanding the Atlantic salmon immune system. Overall, the results indicate that, without the addition of exogenous factors, the adherent HKL cell population differentiates in vitro to become macrophage-like.

A de novo Full-Length mRNA Transcriptome Generated From Hybrid-Corrected PacBio Long-Reads Improves the Transcript Annotation and Identifies Thousands of Novel Splice Variants in Atlantic Salmon.

Atlantic salmon (Salmo salar) is a major species produced in world aquaculture and an important vertebrate model organism for studying the process of rediploidization following whole genome duplication events (Ss4R, 80 mya). The current Salmo salar transcriptome is largely generated from genome sequence based in silico predictions supported by ESTs and short-read sequencing data. However, recent progress in long-read sequencing technologies now allows for full-length transcript sequencing from single RNA-molecules. This study provides a de novo full-length mRNA transcriptome from liver, head-kidney and gill materials. A pipeline was developed based on Iso-seq sequencing of long-reads on the PacBio platform (HQ reads) followed by error-correction of the HQ reads by short-reads from the Illumina platform. The pipeline successfully processed more than 1.5 million long-reads and more than 900 million short-reads into error-corrected HQ reads. A surprisingly high percentage (32%) represented expressed interspersed repeats, while the remaining were processed into 71 461 full-length mRNAs from 23 071 loci. Each transcript was supported by several single-molecule long-read sequences and at least three short-reads, assuring a high sequence accuracy. On average, each gene was represented by three isoforms. Comparisons to the current Atlantic salmon transcripts in the RefSeq database showed that the long-read transcriptome validated 25% of all known transcripts, while the remaining full-length transcripts were novel isoforms, but few were transcripts from novel genes. A comparison to the current genome assembly indicates that the long-read transcriptome may aid in improving transcript annotation as well as provide long-read linkage information useful for improving the genome assembly. More than 80% of transcripts were assigned GO terms and thousands of transcripts were from genes or splice-variants expressed in an organ-specific manner demonstrating that hybrid error-corrected long-read transcriptomes may be applied to study genes and splice-variants expressed in certain organs or conditions (e.g., challenge materials). In conclusion, this is the single largest contribution of full-length mRNAs in Atlantic salmon. The results will be of great value to salmon genomics research, and the pipeline outlined may be applied to generate additional de novo transcriptomes in Atlantic Salmon or applied for similar projects in other species.

Modulation of hepatic miRNA expression in Atlantic salmon (Salmo salar) by family background and dietary fatty acid composition.

This study finds significant differences in hepatic fatty acid composition between four groups of Atlantic salmon (Salmo salar) consisting of offspring from families selected for high and low capacities to express the delta 6 desaturase isomer b and fed diets with 10% or 75% fish oil. The results demonstrated that hepatic lipid metabolism was affected by experimental conditions (diet/family). The fatty acid composition in the four groups mirrored the differences in dietary composition, but it was also associated with the family groups. Small RNA sequencing followed by RT-qPCR identified 12 differentially expressed microRNAs (DE miRNAs), with expression associated with family groups (miR-146 family members, miR-200b, miR-214, miR-221, miR-125, miR-135, miR-137, miR_nov_1), diets (miR-203, miR-462) or both conditions. All the conserved DE miRNAs have been reported as associated with lipid metabolism in other vertebrates. In silico predictions revealed 37 lipid metabolism pathway genes, including desaturases, transcription factors and key enzymes in the synthesis pathways as putative targets (e.g., srebp-1 and 2, Δ6fad_b and c, hmdh, elovl4 and 5b, cdc42). RT-qPCR analysis of selected target genes showed expression changes that were associated with diet and with family groups (d5fad, d6fad_a, srebp-1). There was a reciprocal difference in the abundance of ssa-miR-203a-3p and srebp-1 in one group comparison, whereas other predicted targets did not reveal any evidence of being negatively regulated by degradation. More experimental studies are needed to validate and fully understand the predicted interactions and how the DE miRNAs may participate in the regulation of hepatic lipid metabolism.

Characterization of miRNAs in Extracellular Vesicles Released From Atlantic Salmon Monocyte-Like and Macrophage-Like Cells.

Cell-derived extracellular vesicles (EVs) participate in cell-cell communication via transfer of molecular cargo including genetic material like miRNAs. In mammals, it has previously been established that EV-mediated transfer of miRNAs can alter the development or function of immune cells, such as macrophages. Our previous research revealed that Atlantic salmon head kidney leukocytes (HKLs) change their morphology, phagocytic ability and miRNA profile from primarily "monocyte-like" at Day 1 to primarily "macrophage-like" at Day 5 of culture. Therefore, we aimed to characterize the miRNA cargo packaged in EVs released from these two cell populations. We successfully isolated EVs from Atlantic salmon HKL culture supernatants using the established Vn96 peptide-based pull-down. Isolation was validated using transmission electron microscopy, nanoparticle tracking analysis, and Western blotting. RNA-sequencing identified 19 differentially enriched (DE) miRNAs packaged in Day 1 versus Day 5 EVs. Several of the highly abundant miRNAs, including those that were DE (e.g. ssa-miR-146a, ssa-miR-155 and ssa-miR-731), were previously identified as DE in HKLs and are associated with macrophage differentiation and immune response in other species. Interestingly, the abundance relative of the miRNAs in EVs, including the most abundant miRNA (ssa-miR-125b), was different than the miRNA abundance in HKLs, indicating selective packaging of miRNAs in EVs. Further study of the miRNA cargo in EVs derived from fish immune cells will be an important next step in identifying EV biomarkers useful for evaluating immune cell function, fish health, or response to disease.

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.

Infectious pancreatic necrosis virus (IPNV) infection has been a major problem in salmonid aquaculture. Marker-assisted selection of individuals with resistant genotype at the major IPN quantitative trait locus (IPN-QTL) has significantly reduced mortality in recent years. We have identified host miRNAs that respond to IPNV challenge in salmon fry that were either homozygous resistant (RR) or homozygous susceptible (SS) for the IPN-QTL. Small RNA-sequenced control samples were compared to samples collected at 1, 7, and 20 days post challenge (dpc). This revealed 72 differentially expressed miRNAs (DE miRNAs). Viral load (VL) was lower in RR vs. SS individuals at 7 and 20 dpc. However, analysis of miRNA expression changes revealed no differences between RR vs. SS individuals in controls, at 1 or 7 dpc, while 38 "high viral load responding" miRNAs (HVL-DE miRNAs) were identified at 20 dpc. Most of the HVL-DE miRNAs showed changes that were more pronounced in the high VL SS group than in the low VL RR group when compared to the controls. The absence of differences between QTL groups in controls, 1 and 7 dpc indicates that the QTL genotype does not affect miRNA expression in healthy fish or their first response to viral infections. The miRNA differences at 20 dpc were associated with the QTL genotype and could, possibly, contribute to differences in resistance/susceptibility at the later stage of infection. In silico target gene predictions revealed that 180 immune genes were putative targets, and enrichment analysis indicated that the miRNAs may regulate several major immune system pathways. Among the targets of HVL-DE miRNAs were IRF3, STAT4, NFKB2, MYD88, and IKKA. Interestingly, TNF-alpha paralogs were targeted by different DE miRNAs. Most DE miRNAs were from conserved miRNA families that respond to viral infections in teleost (e.g., miR-21, miR-146, miR-181, miR-192, miR-221, miR-462, miR-731, and miR-8159), while eight were species specific. The miRNAs showed dynamic temporal changes implying they would affect their target genes differently throughout disease progression. This shows that miRNAs are sensitive to VL and disease progression, and may act as fine-tuners of both immediate immune response activation and the later inflammatory processes.

Characterization of Differentially Expressed miRNAs and Their Predicted Target Transcripts during Smoltification and Adaptation to Seawater in Head Kidney of Atlantic Salmon.

Smoltification and early seawater phase are critical developmental periods with physiological and biochemical changes in Atlantic salmon that facilitates survival in saltwater. MicroRNAs (miRNAs) are known to have important roles in development, but whether any miRNAs are involved in regulation of gene expression during smoltification and the adaption to seawater is largely unknown. Here, small RNA sequencing of materials from head kidney before, during smoltification and post seawater transfer were used to study expression dynamics of miRNAs, while microarray analysis was applied to study mRNA expression dynamics. Comparing all timepoints, 71 miRNAs and 2709 mRNAs were identified as differentially expressed (DE). Hierarchical clustering analysis of the DE miRNAs showed three major clusters with characteristic expression changes. Eighty-one DE mRNAs revealed negatively correlated expression patterns to DE miRNAs in Cluster I and III. Furthermore, 42 of these mRNAs were predicted as DE miRNA targets. Gene enrichment analysis of negatively correlated target genes showed they were enriched in gene ontology groups hormone biosynthesis, stress management, immune response, and ion transport. The results strongly indicate that post-transcriptional regulation of gene expression by miRNAs is important in smoltification and sea water adaption, and this study identifies several putative miRNA-target pairs for further functional studies.

Characterization of miRNAs in Cultured Atlantic Salmon Head Kidney Monocyte-Like and Macrophage-Like Cells.

Macrophages are among the first cells to respond to infection and disease. While microRNAs (miRNAs) are involved in the process of monocyte-to-macrophage differentiation in mammals, less is known in teleost fish. Here, Atlantic salmon head kidney leukocytes (HKLs) were used to study the expression of miRNAs in response to in vitro culture. The morphological analysis of cultures showed predominantly monocyte-like cells on Day 1 and macrophage-like cells on Day 5, suggesting that the HKLs had differentiated from monocytes to macrophages. Day 5 HKLs also contained a higher percentage of phagocytic cells. Small RNA sequencing and qPCR analysis were applied to examine the miRNA diversity and expression. There were 370 known mature Atlantic salmon miRNAs in HKLs. Twenty-two miRNAs (15 families) were downregulated while 44 miRNAs (25 families) were upregulated on Day 5 vs. Day 1. Mammalian orthologs of many of the differentially expressed (DE) miRNAs are known to regulate macrophage activation and differentiation, while the teleost-specific miR-2188, miR-462 and miR-731 were also DE and are associated with immune responses in fish. In silico predictions identified several putative target genes of qPCR-validated miRNAs associated with vertebrate macrophage differentiation. This study identified Atlantic salmon miRNAs likely to influence macrophage differentiation, providing important knowledge for future functional studies.

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

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

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

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

Discovery of microRNAs associated with the antiviral immune response of Atlantic cod macrophages.

MicroRNAs (miRNAs) are known to play important immunoregulatory roles in teleosts, although miRNAs involved in the antiviral immune response of Atlantic cod (Gadus morhua) were previously uncharacterised. Using deep sequencing and qPCR, the present study was conducted to identify miRNAs responsive to the viral mimic, polyriboinosinic polyribocytidylic acid (pIC) in Atlantic cod macrophages. Macrophage samples isolated from Atlantic cod (n=3) and treated with pIC or phosphate buffered saline (PBS control) for 24 and 72h were used for miRNA profiling. Following deep sequencing, DESeq2 analyses identified four (miR-731-3p, miR-125b-3-3p, miR-150-3p and miR-462-3p) and two (miR-2188-3p and miR-462-3p) significantly differentially expressed miRNAs at 24 and 72h post-stimulation (HPS), respectively. Sequencing-identified miRNAs were subjected to qPCR validation using a larger number of biological replicates (n=6) exposed to pIC or PBS over time (i.e. 12, 24, 48 and 72 HPS). As in sequencing, miR-731-3p, miR-462-3p and miR-2188-3p showed significant up-regulation by pIC. The sequencing results were not qPCR-validated for miR-125b-3-3p and miR-150-3p as up- and down-regulated miRNAs at 24 HPS, respectively; however, qPCR results showed significant up-regulation in response to pIC stimulation at later time points (i.e. 48 and/or 72 HPS). We also used qPCR to assess the expression of other miRNAs that were previously shown as immune responsive in other vertebrates. qPCR results at 48 and/or 72 HPS revealed that miR-128-3-5p, miR-214-1-5p and miR-451-3p were induced by pIC, whereas miR-30b-3p and miR-199-1-3p expression were repressed in response to pIC. The present study identified ten pIC-stimulated miRNAs, suggesting them as important in antiviral immune responses of Atlantic cod macrophages. Some pIC-responsive miRNAs identified in this study were predicted to target putative immune-related genes of Atlantic cod (e.g. miR-30b-3p targeting herc4), although the regulatory functions of these miRNAs need to be validated by future studies.

Identification of differentially expressed Atlantic salmon miRNAs responding to salmonid alphavirus (SAV) infection.

BACKGROUND: MicroRNAs (miRNAs) control multiple biological processes including the innate immune responses by negative post-transcriptional regulation of gene expression. As there were no studies on the role(s) of miRNAs in viral diseases in Atlantic salmon, we aimed to identify miRNAs responding to salmonid alphavirus (SAV) infection. Their expression were studied at different time points post infection with SAV isolates associated with different mortalities. Furthermore, the genome sequences of the identified miRNAs were analysed to reveal putative cis-regulatory elements, and, finally, their putative target genes were predicted. RESULTS: Twenty differentially expressed miRNAs (DE miRNAs) were identified. The expression of the majority of these increased post infection with maximum levels reached after the viral load were stabilized or decreasing. On the other hand, some miRNAs (e.g. the miRNA-21 family) showed decreased expression at the early time points post infection. There were significant differences in the temporal expression of individual miRNA associated with different SAV isolates. Target gene prediction in SAV responsive immune network genes showed that seventeen of the DE miRNAs could target 24 genes (e.g. IRF3, IRF7). Applying the Atlantic salmon transcriptome as input 28 more immune network genes were revealed as putative targets (e.g. IRF5, IRF4). The majority of the predicted target genes promote inflammatory response. The upstream sequences of the miRNA genes revealed a high density of cis-regulatory sequences known as binding sites for immune network transcription factors (TFs). A high expression in the late phase could therefore be due to increased transcription promoted by immune response activated TFs. Based on the in silico target predictions, we discuss their putative roles as early promotors or late inhibitors of inflammation. We propose that the differences in expressions associated with different SAV isolates could contribute to their differences in mortality rates. CONCLUSIONS: This study represents the first steps in exploring miRNAs important in viral-host interaction in Atlantic salmon. We identified several miRNAs responding to SAV infection. Some likely to prohibit harmful inflammation while other may promote an early immune response. Their predicted functions need to be validated and further studied in functional assays to fully understand their roles in immune homeostasis.

miRNAs associated with immune response in teleost fish.

MicroRNAs (miRNAs) have been identified as important post transcriptional regulators of gene expression. In higher vertebrates, a subset of miRNAs has been identified as important regulators of a number of key genes in immune system gene networks, and this paper review recent studies on miRNAs associated with immune response in teleost fish. Challenge studies conducted in several species have identified differently expressed miRNAs associated with viral or bacterial infection. The results from these studies point out several miRNAs that are likely to have evolutionary conserved functions that are related to immune response in teleost fish. Changed expression levels of mature miRNAs from the five miRNA genes miRNA-462, miRNA-731, miRNA-146, miRNA-181 and miRNA-223 are observed following viral as well as bacterial infection in several teleost fish. Furthermore, significant changes in expression of mature miRNAs from the five genes miRNA-21, miRNA-155, miRNA-1388, miRNA-99 and miRNA-100 are observed in multiple studies of virus infected fish while changes in expression of mature miRNA from the three genes miRNA-122, miRNA-192 and miRNA-451 are observed in several studies of fish with bacterial infections. Interestingly, some of these genes are not present in higher vertebrates. The function of the evolutionary conserved miRNAs responding to infection depends on the target gene(s) they regulate. A few target genes have been identified while a large number of target genes have been predicted by in silico analysis. The results suggest that many of the targets are genes from the host's immune response gene networks. We propose a model with expected temporal changes in miRNA expression if they target immune response activators/effector genes or immune response inhibitors, respectively. The best way to understand the function of a miRNA is to identify its target gene(s), but as the amount of genome resources for teleost fish is limited, with less well characterized genomes and transcriptomes, identifying the true target genes of the miRNAs associated with the immune response is a challenge. Identifying such target genes by applying new methods and approaches will likely be the next important step to understand the function of the miRNAs associated with immune response in teleost fish.

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.

BACKGROUND: Atlantic cod (Gadus morhua) is among the economically most important species in the northern Atlantic Ocean and a model species for studying development of the immune system in vertebrates. MicroRNAs (miRNAs) are an abundant class of small RNA molecules that regulate fundamental biological processes at the post-transcriptional level. Detailed knowledge about a species miRNA repertoire is necessary to study how the miRNA transcriptome modulate gene expression. We have therefore discovered and characterized mature miRNAs and their corresponding miRNA genes in Atlantic cod. We have also performed a validation study to identify suitable reference genes for RT-qPCR analysis of miRNA expression in Atlantic cod. Finally, we utilized the newly characterized miRNA repertoire and the dedicated RT-qPCR method to reveal miRNAs that are highly expressed in certain organs. RESULTS: The discovery analysis revealed 490 mature miRNAs (401 unique sequences) along with precursor sequences and genomic location of the miRNA genes. Twenty six of these were novel miRNA genes. Validation studies ranked gmo-miR-17-1-5p or the two-gene combination gmo-miR25-3p and gmo-miR210-5p as most suitable qPCR reference genes. Analysis by RT-qPCR revealed 45 miRNAs with significantly higher expression in tissues from one or a few organs. Comparisons to other vertebrates indicate that some of these miRNAs may regulate processes like growth, lipid metabolism, immune response to microbial infections and scar damage repair. Three teleost-specific and three novel Atlantic cod miRNAs were among the differentially expressed miRNAs. CONCLUSIONS: The number of known mature miRNAs was considerably increased by our identification of miRNAs and miRNA genes in Atlantic cod. This will benefit further functional studies of miRNA expression using deep sequencing methods. The validation study showed that stable miRNAs are suitable reference genes for RT-qPCR analysis of miRNA expression. Applying RT-qPCR we have identified several miRNAs likely to have important regulatory functions in particular organs.

Identification of Endogenous Controls for Use in miRNA Quantification in Human Cancer Cell Lines.

BACKGROUND: miRNAs play important roles in multiple biological processes, and deregulation has been linked to several human diseases, including cancer. Studying changes in miRNA expression in cancer development is commonly performed in vitro in human cancer cell lines using quantitative polymerase chain reaction (qPCR), a method requiring the use of a robust reference gene that displays stable expression across all samples. MATERIALS AND METHODS: Using the NormFinder software, a selection of commonly used endogeneous controls and miRNAs were tested in six human cancer cell lines to identify for the most suitable gene for use as a reference. RESULTS: The frequently used endogenous control U6B small nuclear RNA (RNU6B) was found to be an unsuitable reference for normalization. The most suitable single endogeneous control identified was miR-25-3p, whereas the best combination of two endogeneous controls was miR-25-3p and miR-93-5p. CONCLUSION: We identified a single and a pair of miRNAs suitable for use as endogenous controls when performing qPCR-based miRNA expression analyses in human cancer cell lines.