Circulating miRNA measurements are reflective of
cholesterol-based changes in rainbow trout
(Oncorhynchus mykiss)

Tengfei Zhu, Geneviève Corraze, Elisabeth Plagnes-Juan, Sandrine

Skiba-Cassy

► To cite this version:

Tengfei Zhu, Geneviève Corraze, Elisabeth Plagnes-Juan, Sandrine Skiba-Cassy. Circulating miRNA measurements are reflective of cholesterol-based changes in rainbow trout (Oncorhynchus mykiss). PLoS ONE, 2018, 13 (11), pp.e0206727. 10.1371/journal.pone.0206727. hal-01938700

HAL Id: hal-01938700

https://hal.science/hal-01938700

Submitted on 28 Nov 2018

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Circulating miRNA measurements are reflective of cholesterol-based changes in rainbow trout (*Oncorhynchus mykiss*)

Tengfei Zhu, Geneviève Corraze, Elisabeth Plagnes-Juan, Sandrine Skiba-Cassy*

INRA, Univ Pau & Pays Adour, E2S UPPA, UMR 1419, Nutrition Métabolisme Aquaculture, Saint Pée sur Nivelle, France

* sandrine.skiba@inra.fr

<figure><img src="image_2.png" /></figure>

## Abstract

MicroRNAs (miRNAs) are a class of small non-coding RNAs which are known to posttranscriptionally regulate the expression of most genes in both animals and plants. Meanwhile, studies have shown that numbers of miRNAs are present in body fluids including the plasma. Despite the mode of action of these circulating miRNAs still remains unknown, they have been found to be promising biomarkers for disease diagnosis, prognosis and response to treatment. In order to evaluate the potential of miRNAs as non-invasive biomarkers in aquaculture, a time-course experiment was implemented to investigate the postprandial regulation of miRNAs levels in liver and plasma as well as the hepatic expression of genes involved in cholesterol metabolism. We showed that miR-1, miR-33a, miR-122, miR-128 and miR-223 were expressed in the liver of rainbow trout and present at detectable level in the plasma. We also demonstrated that hepatic expression of miR-1, miR-122 and miR-128 were regulated by feed intake and reached their highest levels 12 hours after the meal. Interestingly, we observed that circulating levels of miR-128 and miR-223 are subjected to postprandial regulations similar to that observed in their hepatic counterparts. Statistical correlations were observed between liver and plasma for miR-128 and miR-223 and between hepatic and circulating miR-122, miR-128 and miR-223 and expression of genes related to cholesterol synthesis and efflux or glucose phosphorylation. These results demonstrated that circulating miR-122, miR-128 and miR-223 are potential biomarkers of cholesterol metabolism in rainbow trout.

<figure><img src="image_3.png" /></figure>

## OPEN ACCESS

**Citation:** Zhu T, Corraze G, Plagnes-Juan E, Skiba-Cassy S (2018) Circulating miRNA measurements are reflective of cholesterol-based changes in rainbow trout (*Oncorhynchus mykiss*). PLoS ONE 13(11): e0206727. https://doi.org/10.1371/journal.pone.0206727

**Editor:** José L. Soengas, Universidade de Vigo, SPAIN

**Received:** May 15, 2018

**Accepted:** October 18, 2018

**Published:** November 5, 2018

**Copyright:** © 2018 Zhu et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** T. Zhu gratefully acknowledges the financial assistance provided by the China Scholarship Council (CSC, File No. 201406330071) for his doctoral fellowship.

**Competing interests:** The authors have declared that no competing interests exist.

## Introduction

MicroRNAs (miRNAs) are a class of small non-coding RNAs of 19-24-mer nucleotides, mostly regulating gene expression at posttranscriptional level by mRNA decay or translational repression in animals and plants [1]. The high rate of miRNAs sequence conservation exhibited in many distantly related plant [2] and animal [3] species indicates their important functions in biological processes such as development, differentiation, proliferation, stress responses, apoptosis and oncogenesis [4–9].
Originating from independent transcription units, most miRNA genes are expressed under the control of their own promoters and regulatory sequences, nonetheless, others, located in gene introns are co-regulated with their host genes [4,5]. Initially transcribed by RNA polymerase II, a miRNA gene generates a long transcript called primary miRNA (pri-miRNA) in the nucleus which is then cleaved by a complex formed by RNase III enzyme Drosha and its binding partner DiGeorge syndrome critical region 8 (DGCR8) into a precursor miRNA (pre-miRNA). Following Exportin 5-dependent nuclear export, pre-miRNA is processed by RNase III enzyme Dicer and its interacting partner transactivation-responsive RNA-binding protein (TRBP) to generate a small double-stranded miRNA duplex which contains both the mature miRNA strand (guide strand) and its complementary strand (passenger strand). Finally, the miRNA duplex dissociates and the mature strand is carried by the Argonaute (AGO) protein forming the RNA-induced silencing complex (RISC), which decreases gene expression by mRNA decay or translational repression when perfect or imperfect complementarity occurs between miRNA and the targeted mRNA, respectively [10,11].

Pioneering works from Chim *et al*. found the presence of placental miRNAs in the maternal blood plasma [12] while, for the first time, Lawrie *et al*. suggested the potential utility of miRNAs in serum as non-invasive biomarkers for cancer diagnosis and prognosis [13]. Later, studies reported the presence of miRNAs in other body fluids, including plasma, tears, saliva, urine, breast milk, semen, amniotic fluid, colostrum, bronchial lavage, cerebrospinal fluid, peritoneal fluid, and pleural fluid [14–18], suggesting that circulating miRNAs may be involved in cell-to-cell communication [19]. It is also noteworthy that miRNAs are differently expressed in tissues as well as during developmental stages. For example, miR-122 is enriched in liver [20] whereas miR-133a and miR-133b are preferentially expressed in muscle [21]. miR-427 is reported to be transcribed only during blastula and gastrula stages of *Xenopus* development [22]. Since miRNAs were detected in various body fluids, it has been proposed to use circulating miRNAs as non-invasive biomarkers for disease diagnosis, prognosis and response to treatment [15,23]. Thus, the expression level of miR-141 in serum can thereby serve to distinguish prostate cancer patients from healthy controls [24] while it has also been shown that the expression level of miR-1 in serum is significantly higher in patients with acute myocardial infarction [25,26].

In contrast with mRNA, extracellular miRNAs are quite stable in nuclease- and protease-rich environment. Studies have shown that a majority of extracellular miRNAs can be protected from degradation by interacting with AGO proteins, which have a remarkable stability in extracellular environment [27,28]. Moreover, it has been demonstrated that miRNA stability is increased when miRNAs are engulfed in lipid vesicles, including apoptotic bodies, shedding vesicles, exosomes and high density lipoprotein (HDL) particles [29–33]. Although the selective release mechanism of circulating miRNAs as well as whether they could mediate cell-to-cell signaling remain unclear, they still appear as promising biomarkers thanks to their unique sequences, high stabilities in body fluids and tissue-specificities [34].

MiRNAs also play important roles in the development and physiological processes in fish [35–37]. First and most characterized in zebrafish [38][39][40][35], miRNAs are not much studied in fish compared to mammals. Beyond its importance in aquaculture, rainbow trout (*Oncorhynchus mykiss*) is also a model fish species worldwide studied for research in carcinogenesis, toxicology, comparative immunology, disease ecology, physiology and nutrition [41]. Therefore, studies of miRNAs in rainbow trout are of great interest for these scientific fields. In this context and during the last years, our laboratory and others initiated series of experiments in order to address the physiological importance of miRNAs in rainbow trout. In 2008, 14 miRNAs were identified and first reported in unfertilized eggs and early stage embryos of rainbow trout [42]. Then, tissue specificities, targets predictions and genetic polymorphisms of
miRNAs in rainbow trout were studied [43]. Subsequently, more miRNAs were identified in different rainbow trout tissues [44,45] and miRNAs involvement in development [46] and metabolic processes [47–49] were finally reported in rainbow trout.

Fishmeal and fish oil, the two traditional ingredients in fish feeds, are being replaced by plant ingredients in a much higher proportion to support the sustainable development of aquaculture [50]. Consequently, the nutrient composition of plant-based diet has greatly changed compared with diets mainly based on fishmeal and fish oil [51][52]. For example, cholesterol, which is normally mainly derived from animal ingredients, is almost absent in plant ingredients [53]. Previous studies on fish have shown that the expression of cholesterol related genes is affected when fishmeal and fish oil are replaced by plant ingredients [54–56]. Furthermore, when cholesterol is added to fish diet, plasma cholesterol level and the expression of genes related to cholesterol metabolism are also affected [57,58], indicating a significant effect of diet on cholesterol metabolism in fish.

In mammals, several miRNAs involved in the regulation of cholesterol or lipid metabolism have been identified, including miR-1, -33a, -122, -128 and -223. For instance, miR-1 is reported to repress *lxra*, a master regulator of cholesterol homeostasis known to also target genes related to lipogenesis [59]. miR-33a is shown to act in concert with its host gene sterol regulatory element-binding protein (*srebp*) -2 to target ATP-binding cassette transporter A1 (*abca1*) [60] and to regulate genes involved in cholesterol and lipid metabolism [61]. The inhibition of miR-122 in mice, which represents the most abundant miRNA in mammalian liver [20], revealed that it is a key regulator of cholesterol and fatty acid metabolism [62]. Accordingly with these results, we previously observed that miR-122 inhibition resulted in the decrease of cholesterol in plasma and the increase of gene expression involved in cholesterol degradation and excretion in rainbow trout [49]. Interestingly; miR-128 is shown to alter the expression of multiple genes involved in cellular cholesterol homeostasis through direct targeting of ATP-binding cassette transporters A1 and G1 mRNAs(*abca1* and *abcg1* respectively) and retinoid X receptor (*rxra*) [63]. Finally, the expression of miR-223 has been shown to inhibit cholesterol biosynthesis through the direct repression of sterol enzymes 3-hydroxy-3-methylglutaryl-CoA synthase (*hmgcs*) and methylsterol monooxygenase 1 (*sc4mol*) mRNAs and also to inhibit HDL-cholesterol uptake by scavenger receptor B1 (SR-B1) repression [64,65]. We also observed that miR-223 expression is suppressed, together with the expression of liver X receptor α (*lxra*) in trout fed a plant-based diet while the expression of *srebp-2* is enhanced [56].

Therefore, we decided to evaluate if the abundance of circulating miRNAs could reflect the level of expression of these miRNAs in liver as well as those of related target genes. We undertook a parallel analysis of the hepatic expression and plasma abundance of miRNAs together with the hepatic expression of genes related to cholesterol metabolism in rainbow trout. The further determination of correlations between plasma and hepatic miRNAs as well as between circulating miRNAs and genes then helped us to evaluate the capacity of circulating miRNAs to be used as non-invasive biomarkers of tissue miRNA and gene expression.

## Materials and methods

### Ethics statement

Experiments were carried out in the INRA experimental facilities (UMR1419 Nutrition, Métabolisme, Aquaculture, Donzacq, France) authorized for animal experimentation by the French veterinary service which is the competent authority (A 64-495-1). Experiments were in strict accordance with EU legal frameworks related to the protection of animals used for scientific research (Directive 2010/63/EU) and according to the National Guidelines for Animal Care of
the French Ministry of Research (decree n°2013-118, February 1st, 2013). Scientists in charge of the experimentation received a training and a personal authorization (N°B64 10 005). In agreement with the ethical committee "Comité d'Ethique Aquitaine Poissons Oiseaux" (C2EA-73), the present study does not need approval by a specific ethical committee since it implies only classical rearing practices with all diets formulated to cover the nutritional requirements of rainbow trout [33]. During this study, fish were daily monitored. If any clinical symptoms (i.e. morphological abnormality, restlessness or uncoordinated movements) were observed, fish were sedated by immersion in 10mg/L benzocaine solution and then euthanized by immersion in a 60mg/L benzocaine solution (anesthetic overdose) during 3 minutes.

## Sampling procedure

Rainbow trout (*Oncorhynchus mykiss*) (50g mean body weight) were reared in our experimental fish farm (INRA, Donzacq, France-permit n°A64-104-1) in open circuit tanks supplied with spring water at 17°C and under natural photoperiod. Fish were fed by hand twice a day with a commercial diet (Skretting, protein 42.5%, lipid 20.5%). At 3h, 8h, 12h, 16h and 24h after the last meal, six fish were anesthetized with benzocaine (30mg/L) and sacrificed by anesthetic overdose for sampling. Blood was removed from the caudal vein into syringes, rinsed with 10% EDTA and centrifuged (3000g, 5min). The recovered plasma was immediately frozen and kept at -80°C for miRNA expression and plasma metabolites analyses. Livers were dissected, weighted and 100 mg of sliced liver were sampled, put on absorbent paper to remove blood before being frozen in liquid nitrogen and kept at -80°C for hepatic gene and miRNA expression analyses.

## Plasma metabolites analysis

Plasma glucose, triglycerides and cholesterol concentrations were measured using commercial kits (BioMerieux, Marcy l'Etoile, France) adapted to microplate format, according to manufacturer's recommendations.

## Hepatic gene expression analysis

For liver samples, total RNA was extracted as previously described [46] using Trizol reagent (Invitrogen, Carlsbad, CA, USA) according to manufacturer's instructions and were quantified by spectrophotometry (absorbance at λ = 260nm). The integrity of ARN samples was assessed by means of agarose gel electrophoresis. 1 µg of total RNA was used for cDNA synthesis. The SuperScript III RNaseH-reverse transcriptase kit (Invitrogen) with oligo dT random primers (Promega, Charbonnières, France) was used to synthesize cDNA (n = 6). The primer sequences used for Reverse Transcription-quantitative Polymerase Chain Reaction RT-qPCR analyses are listed in the Table 1. RT-qPCR gene expression analyses were focused on genes involved in cholesterol metabolism, including ATP-binding cassette transporter A1 (*abca1*), ATP-binding cassette transporter G5 (*abcg5*), ATP-binding cassette transporter G8 (*abcg8*), Acyl-CoA cholesterol acyltransferase 1 (*acat1*), Acyl-CoA cholesterol acyltransferase 2 (*acat2*), Lanosterol 14α-demethylase (*cyp51*), cholesterol 7α-hydroxylase (*cyp7a1*), 7-dehydrocholesterol reductase (*dhcr7*), HMG-CoA reductase (*hmgcr*), HMG-CoA synthase (*hmgcs*), liver X receptor α (*lxra*), sterol regulatory element-binding protein 2 (*srebp-2*) and UDP glycuronosyltransferase 1 A3 (*ugt1a3*); lipogenesis, including fatty acid synthase (*fas*) and sterol regulatory element-binding protein 1c (*srebp-1c*) and glycolysis, glucokinase (*gck*).

Primers were designed to coamplify different paralogs of the same gene.
<table><caption>Table 1. Sequence of the primers in qRT-PCR for gene expression analysis.</caption><thead><tr><th rowspan="2"></th><th rowspan="2">Forward primer</th><th rowspan="2">Reverse primer</th><th colspan="2">Genoscope<fnref n="1" /> or Genbank accession number</th></tr><tr><th>Paralogue 1</th><th>Paralogue 2</th></tr></thead><tbody><tr><td>abca1</td><td>CAGGAAAGACGAGCACCTT</td><td>TCTGCCACCTCACACTTC</td><td>GSONMG00078741001</td><td>GSONMG00074045001</td></tr><tr><td>abcg5</td><td>CACCGACATGGAGACAGAAA</td><td>GACAGATGGAAGGGGATGAA</td><td>GSONMG00075025001</td><td>/</td></tr><tr><td>abcg8</td><td>GATACCAGGGTTCCAGAGCA</td><td>CCAGAAACAGAGGGACCAGA</td><td>GSONMG00075024001</td><td>/</td></tr><tr><td>acat1</td><td>GGCAAGCCTGATGTGGTAGT</td><td>ACCGTGCCATTCTCCTTCTG</td><td>GSONMG00033938001<br>D2OM Loc100702912.1.1</td><td>/</td></tr><tr><td>acat2</td><td>TGCTTGTTGTCCCTGGGTTT</td><td>GTGTGGCTGTGACGTGTTTC</td><td>GSONMG00054446001<br>D2OM THIC.1.1</td><td>/</td></tr><tr><td>cyp51</td><td>CCCGTTGTCAGCTTTACCA</td><td>GCATTGAGATCTTCGTCTTGC</td><td>GSONMG00031182001</td><td>GSONMG00044416001</td></tr><tr><td>cyp7a1</td><td>ACGTCCGAGTGGCTAAAGAG</td><td>GGTCAAAGTGGAGCATCTGG</td><td>GSONMG00066448001<br>AB675933.1</td><td>GSONMG00037174001<br>AB675934.1</td></tr><tr><td>dhcr7</td><td>GTAACCCACCAGACCCAAGA</td><td>CCTCTCCTATGCAGCCAAC</td><td>GSONMG00025402001</td><td>GSONMG00039624001</td></tr><tr><td>fas</td><td>TGATCTGAAGGCCCGTGTCA</td><td>GGGTGACGTTGCCGTGGTAT</td><td>GSONMG00062364001</td><td></td></tr><tr><td>gck</td><td>GCACGGCTGAGATGCTCTTTG</td><td>GCCTTGAACCCTTTGGTCCAG</td><td>GSONMG00033781001</td><td>GSONMG00012878001</td></tr><tr><td>hmgcr</td><td>GAACGGTGAATGTGCTGTGT</td><td>GACCATTTGGGAGCTTGTGT</td><td>GSONMG00016350001</td><td>/</td></tr><tr><td>hmgcs</td><td>AGTGGCAAAGAGAGGGTGTG</td><td>TTCTGGTTGGAGACGAGGAG</td><td>GSONMG00010243001<br>D2OM Loc100712075.1.1</td><td>/</td></tr><tr><td>lxra</td><td>TGCAGCAGCCGTATGTGGA</td><td>GCGGCGGGAGCTTCTTGTTC</td><td>GSONMG00014026001</td><td>GSONMG00064070001</td></tr><tr><td>srebp-1c</td><td>CATGCGCAGGTTGTTTCTT</td><td>GATGTGTTCGTGTGGGACTG</td><td>XM_021624594.1</td><td></td></tr><tr><td>srebp-2</td><td>TAGGCCCCAAAGGGATAAAG</td><td>TCAGACACGACGAGCACAA</td><td>GSONMG00039651001</td><td>GSONMG00061885001</td></tr><tr><td>ugt1a3</td><td>CCACCAGCAAGACAGTCTCA</td><td>CAACAGCACAGTGGCTGACT</td><td>GSONMG00035844001</td><td>/</td></tr></tbody></table>

RT-qPCR assays were performed on the Roche LightCycler 480 II system (Roche Diagnostics, Neuilly sur Seine, France). The assays were carried out using a reaction mix of 6 µL per sample containing 2 µL of 76 times diluted cDNA, 0.24 µL of each primer (10 µM), 3 µL of LightCycler 480 SYBR Green I Master mix and 0.52 µL DNAse/RNAse free water (5 Prime GmbH, Hamburg, Germany). The PCR protocol was initiated at 95°C for 10 min for initial denaturation of the cDNA and hot-start Taq-polymerase activation, followed by 45 cycles of a three-step amplification program (15s at 95°C, 10s at melting temperature Tm (60-65°C), 15s at 72°C), according to the primer set used. Melting curves were systematically monitored (5 s at 95°C, 1 min at 65°C, temperature gradient at 0.11°C/s from 65 to 97°C) at the end of the last amplification cycle to confirm the specificity of the amplification reaction. Each PCR assay included replicate samples (duplicate of reverse transcription and PCR amplification) and negative controls (RT- and cDNA-free samples, respectively). Relative quantification of target gene expression was determined using the E-Method from the LightCycler 480 software (version SW 1.5; Roche Diagnostics). The NORMA-Gene algorithm not using housekeeping genes was used for the normalization by calculating the average deviation from a sample to the mean of specific genes [66] as the expression of the 3 reference genes tested (β-actin, EF1, 18 S) was not consistent between different sampling time points. In all cases, PCR efficiency (E) measured by the slope of a standard curve with serial dilutions of cDNA ranged between 1.8 and 2.0.

## Hepatic miRNA expression analysis

Expression of hepatic microRNAs (miR-1-3p, miR-33a-5p, miR-122-5p, miR-128-3p and miR-223-3p) was assessed by means of the TaqMan Advanced miRNA Assays (A25576,

<footnote n="1">¹ https://www.genoscope.cns.fr/trout</footnote>

<footnote>https://doi.org/10.1371/journal.pone.0206727.t001</footnote>
Applied Biosystems, ThermoFisher) according to manufacturer's instructions. Synthetic cel-miR-39 mimic (miScript miRNA Mimics, QIAGEN) was added to each RNA sample and used as an exogenous control for normalization. It exhibited a constant level in liver samples at the different sampling time points. Briefly, 80 ng of total RNA were used for the poly(A) tailing, ligation and reverse transcription reactions to synthesize the cDNA of all miRNAs followed by a pre-amplification step. qPCR was then performed in a reaction mix of 6 µL containing 2 µL cDNA (200 times diluted), 2.67 µL 2X Fast Advanced Master mix (Applied Biosystems, USA), 0.27 µL TaqMan Advanced miRNA Assay (20X) (Applied Biosystems, USA) and 1.06 µL DNAse/RNAse free water (5 Prime GmbH, Hamburg, Germany). Sequences of the miRNA probes were listed in Table 2.

The PCR protocol was initiated at 95°C for 20s for initial denaturation of the cDNA and the enzyme activation, followed by 50 cycles of a 2 steps amplification program (3s at 95°C for denaturation, 30s at 60°C for annealing). Each PCR assay included replicates for each sample (duplicates of reverse transcription and PCR amplification) and also negative controls (reverse transcriptase free and RNA free samples). Relative quantification of the target miRNAs was determined using the E-Method from the LightCycler 480 software (version SW 1.5; Roche Diagnostics). PCR efficiency measured by the slope of a standard curve with serial dilutions of miRNA cDNA ranged between 1.9 and 2.1.

## Plasmatic miRNA expression analysis

miR-1-3p, miR-33a-5p, miR-122-5p, miR-128-3p and miR-223-3p were detected in plasma samples by means of the TaqMan Advanced miRNA Assays (A25576, Applied Biosystems, ThermoFisher). Synthetic cel-miR-39 mimic (miScript miRNA Mimics, QIAGEN) was added and used as an exogenous control for normalization and showed constant levels in plasma samples.

Total RNA was extracted from plasma samples with the TRIzol LS reagent (Life Technologies, Carlsbad, CA, USA) according to manufacturer's instructions and was quantified by spectrophotometry (absorbance at λ = 260nm). The cDNA synthesis of all miRNAs and the following qPCR steps were performed in the same way as those previously described for hepatic miRNA expression analysis. PCR efficiency measured by the slope of a standard curve with serial dilutions of miRNA cDNA was nearly 2.0.

## Statistical analysis

Results are expressed as means ± SD. Statistical analyses were carried out using one-way ANOVA to detect significant differences, followed by a Turkey test for post hoc analysis. Normality was assessed using the Shapiro-Wilk test, while homogeneity of variance was determined using Levene's test. In cases where data were nonparametric or not homoscedastic, data transformations, such as logarithms, were used to meet ANOVA criteria. For all statistical

<table><caption>Table 2. Sequences of probes used for TaqMan advanced miRNA assays.</caption><thead><tr><td>Assay</td><td>Assay ID</td><td>Probe sequence</td></tr></thead><tbody><tr><td>miR-1-3p</td><td>477820_mir</td><td>TGGAATGTAAAGAAGTATGTAT</td></tr><tr><td>miR-33a-5p</td><td>478347_mir</td><td>GTGCATTGTAGTTGCATTGCA</td></tr><tr><td>miR-39-3p</td><td>478293_mir</td><td>UCACCGGGUGUAAAUCAGCUUG</td></tr><tr><td>miR-122-5p</td><td>477855_mir</td><td>TGGAGTGTGACAATGGTGTTTG</td></tr><tr><td>miR-128-3p</td><td>477892_mir</td><td>TCACAGTGAACCGGTCTCTTT</td></tr><tr><td>miR-223-3p</td><td>477983_mir</td><td>TGTCAGTTTGTCAAATACCCCA</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0206727.t002
analysis, the level of significance was set at *P*<0.05. Pearson correlation coefficients were calculated based on data of normalized genes or miRNA expression calculated by the E-Method from the LightCycler 480 software (version SW 1.5; Roche Diagnostics). All statistical analyses were performed using R software [67].

## Results

The complete dataset of the study is presented within the manuscript and in its Supporting Information files (S1 File).

### Postprandial cholesterol, glucose and triglyceride levels in the plasma of rainbow trout

No significant changes were found for cholesterol, triglycerides and glucose levels between plasma samples collected at different time points after the meal (Fig 1).

### Postprandial hepatic gene expression

The expression of the four main genes involved in cholesterol synthesis (*hmgcr*, *hmgcs*, *cyp51* and *dhcr7*) and the transcription factor *srebp-2* significantly increased reaching their highest level of expression 16h after the meal (*P*<0.05). The expression of these genes returned to basal level 24h after the meal with similar expressions compared with samples collected 3h after the meal (Fig 2). Pearson correlation analysis revealed that the expression of *srebp-2* and its four target genes (*hmgcr*, *hmgcs*, *cyp51*, *dhcr7*) were all positively correlated (Table 3).

As observed for genes involved in cholesterol synthesis, *ugt1a3* responsible for the glucuronidation of bile acid, *abcg5* and *abcg8*, implicated in biliary excretion of cholesterol, exhibited highest level of expression 16h after the meal (*P*<0.05) (Fig 3). However, no significant postprandial changes were found for the expression of *cyp7a1* which is the key enzyme of bile acid synthesis, *abca1* which mediates the efflux of cholesterol and phospholipids to lipid-poor apolipoproteins to form nascent high-density lipoproteins and the transcription factor *lxra* (*P*>0.05) (Fig 3). Meanwhile, the postprandial expression of *ugt1a3*, *abcg5*, and *abcg8* were positively correlated with each other and the expression of *lxra* was found to be positively correlated with *cyp7a1* and *abcg5* as well (Table 3). The expression of *abcg8* positively correlated with all genes involved in cholesterol synthesis, while *abcg5* correlated with all genes involved in cholesterol efflux in a positive way except for *abca1* in a negative way (Table 3). The expression of *abcg5* was positively correlated with plasma cholesterol and triglycerides whereas an

<figure><img src="image_4.png" /><figcaption>**Fig 1.** Postprandial levels of cholesterol (A), triglycerides (B) and glucose (C) in plasma of rainbow trout at 3h, 8h, 12h, 16h and 24h after the meal. Values are means (n = 9), with their standard deviations represented by vertical bars.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0206727.g001
<figure><img src="image_5.png" /><figcaption>**Fig 2.** Postprandial expression of genes involved in cholesterol synthesis including *hmgcs* (A), *hmgcr* (B), *cyp51* (C), *dhcr7* (D), *srebp-2* (E) in the liver of rainbow trout at 3h, 8h, 12h, 16h and 24h after the meal. The NORMA-Gene method was used for the normalization. Values are means (n = 6), with their standard deviations represented by vertical bars. a, b, c Mean values with unlike letters were significantly different (*P*<0.05; One-way analysis of variance, Tukey's test).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0206727.g002

opposite correlation was found between plasma cholesterol and triglycerides and *abca1* (Table 4).

Regarding genes involved in cholesterol esterification, the highest expression of *acat1* and *acat2* was observed 16h after the meal (*P*<0.05) and then decreased to their basal level 24h after feeding (Fig 4). Both *acat1* and *acat2* were positively correlated with each other but also with genes involved in cholesterol synthesis. In addition, *acat1* was found to positively correlate with *ugt1a3*, *abcg5*, and *abcg8*, genes involved in cholesterol efflux (Table 3).

Expression of the lipogenic gene *fas* and its transcription factor *srebp-1c* significantly increased 16h after the meal (*P*<0.05) and then decreased between 16h and 24h. *gck* which is considered as the first enzyme of lipogenesis that induces glucose phosphorylation showed a similar pattern (Fig 5). *fas* was positively correlated with *srebp-1c* and genes involved in cholesterol synthesis, efflux and esterification (Table 3) as well as with the level of cholesterol and triglycerides in the plasma (Table 4).

## Postprandial miRNAs expression

The highest and lowest levels of hepatic miR-1, miR-122 and miR-128 were observed 12h and 16h after the meal, respectively (*P*<0.05). No significant changes were found for miR-33a and miR-223 during the postprandial period (Fig 6).

In the meantime, we analyzed the miRNA abundance in plasma of the same fish. Whereas no postprandial variation was observed for miR-1, miR-33 and miR-122 in plasma samples, our results indicated a significant rise of miR-128 abundance in the plasma at 12h of the postprandial time course before to significantly decrease to its lowest abundance 4h after. Interestingly and despite the absence of difference of miR-233 in postprandial livers, circulating miR-
<table><caption>Table 3. Correlations between gene expressions involved in cholesterol metabolism of rainbow trout.</caption><thead><tr><th rowspan="2">Pearson value<br>n = 29</th><th colspan="5">Cholesterol synthesis</th><th colspan="2">Bile acid synthesis</th><th colspan="3">Cholesterol efflux</th><th colspan="3">Esterification</th><th colspan="3">Lipogenesis</th></tr><tr><th>hmgcs</th><th>hmocr</th><th>cyp51</th><th>dhcr7</th><th>srebp-2</th><th>cyp7a1</th><th>ugt1a3</th><th>abcg5</th><th>abcg8</th><th>abca1</th><th>lxr</th><th>acat1</th><th>acat2</th><th>fas</th><th>srebp-1c</th><th>gck</th></tr></thead><tbody><tr><td>Cholesterol synthesis</td><td>hmgcs</td><td></td><td>0.942**</td><td>0.957**</td><td>0.870**</td><td>0.830**</td><td>-0.1573</td><td>0.1301</td><td>0.1857</td><td>0.565**</td><td>-0.1408</td><td>-0.3116</td><td>0.433*</td><td>0.806**</td><td>0.848**</td><td>0.793**</td><td>0.432*</td></tr><tr><td></td><td>hmocr</td><td></td><td></td><td>0.979**</td><td>0.875**</td><td>0.863**</td><td>-0.1811</td><td>0.361*</td><td>0.2173</td><td>0.667**</td><td>-0.1816</td><td>-0.3186</td><td>0.528**</td><td>0.776**</td><td>0.870**</td><td>0.710**</td><td>0.517**</td></tr><tr><td></td><td>cyp51</td><td></td><td></td><td></td><td>0.898**</td><td>0.895**</td><td>-0.1841</td><td>0.2825</td><td>0.1692</td><td>0.615**</td><td>-0.1133</td><td>-0.3047</td><td>0.487**</td><td>0.775**</td><td>0.841**</td><td>0.729**</td><td>0.471**</td></tr><tr><td></td><td>dhcr7</td><td></td><td></td><td></td><td></td><td>0.742**</td><td>-0.0174</td><td>0.3039</td><td>0.3497</td><td>0.696**</td><td>-0.2184</td><td>-0.1253</td><td>0.501**</td><td>0.566**</td><td>0.840**</td><td>0.641**</td><td>0.518**</td></tr><tr><td></td><td>srebp-2</td><td></td><td></td><td></td><td></td><td></td><td>-0.3235</td><td>0.2252</td><td>-0.0283</td><td>0.533**</td><td>0.0893</td><td>-0.2301</td><td>0.569**</td><td>0.711**</td><td>0.701**</td><td>0.762**</td><td>0.411*</td></tr><tr><td>Bile acid synthesis</td><td>cyp7a1</td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.2908</td><td>0.549**</td><td>0.247</td><td>-0.05</td><td>0.409*</td><td>0.0123</td><td>-0.2376</td><td>-0.1453</td><td>-0.2186</td><td>0.2128</td></tr><tr><td></td><td>ugt1a3</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.493**</td><td>0.583**</td><td>-0.2284</td><td>0.061</td><td>0.491**</td><td>0.0115</td><td>0.377*</td><td>0.0189</td><td>0.630**</td></tr><tr><td>Cholesterol efflux</td><td>abcg5</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.473**</td><td>-0.373*</td><td>0.431*</td><td>0.449*</td><td>0.0093</td><td>0.3446</td><td>0.1457</td><td>0.532**</td></tr><tr><td></td><td>abcg8</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>-0.2559</td><td>-0.1003</td><td>0.450*</td><td>0.289</td><td>0.620**</td><td>0.2953</td><td>0.619**</td></tr><tr><td></td><td>abca1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.0935</td><td>-0.0589</td><td>-0.1452</td><td>-0.373*</td><td>-0.1535</td><td>-0.3369</td></tr><tr><td></td><td>lxr</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.2506</td><td>-0.3376</td><td>-0.3017</td><td>-0.1498</td><td>-0.0056</td></tr><tr><td>Esterification</td><td>acat1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.393*</td><td>0.508**</td><td>0.585**</td><td>0.470**</td><td></td></tr><tr><td></td><td>acat2</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.622**</td><td>0.760**</td><td>0.2371</td><td></td></tr><tr><td>Lipogenesis</td><td>fas</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.757**</td><td>0.603**</td><td></td></tr><tr><td></td><td>srebp-1c</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.3612</td><td></td></tr><tr><td></td><td>gck</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

P-values

* P < 0.05

** P < 0.01

Pearson correlation coefficient larger than zero represents a positive correlation, while Pearson correlation coefficient less than zero represents a negative correlation.

https://doi.org/10.1371/journal.pone.0206727.t003
<figure><img src="image_6.png" /><figcaption>**Fig 3.** Postprandial expression of genes involved in cholesterol efflux including *cyp7a1* (A), *ugt1a3* (B), *abcg5* (C), *abcg8* (D), *abca1* (E) *lxra* (F) in the liver of rainbow trout at 3h, 8h, 12h, 16h and 24h after the meal. The NORMA-Gene method was used for the normalization. Values are means (n = 6), with their standard deviations represented by vertical bars. <sup>a, b</sup> Mean values with unlike letters were significantly different (<i>P</i><0.05; One-way analysis of variance, Tukey's test).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0206727.g003

223 peaked 8h after the meal and remained stable for 4 more hours before to decrease to its lowest level 16h after the meal (*P*<0.05), (Fig 7).

## Correlations between plasmatic parameters and miRNAs in liver and plasma

The expression of miR-33a in liver was positively correlated with the level of cholesterol and triglycerides in plasma (*P*<0.01). On the other hand, circulating miR-223 negatively correlated with the level of cholesterol and triglycerides (*P*<0.05). Meanwhile, the plasma miR-122 and glucose levels showed positive correlation between each other (*P*<0.05) (Table 5).

<table><caption>**Table 4.** Correlations between hepatic gene expressions involved in cholesterol metabolism and plasma parameters of rainbow trout.</caption><thead><tr><th rowspan="2">Pearson value</th><th colspan="5">Cholesterol synthesis</th><th colspan="5">Cholesterol efflux</th><th colspan="2">Esterification</th><th colspan="2">Lipogenesis</th><th>Glycogenesis</th></tr><tr><th><i>hmgcs</i></th><th><i>hmgcr</i></th><th><i>cyp51</i></th><th><i>dhcr7</i></th><th><i>srebp-2</i></th><th><i>cyp7a1</i></th><th><i>ugt1a3</i></th><th><i>abcg5</i></th><th><i>abcg8</i></th><th><i>abca1</i></th><th><i>lxr</i></th><th><i>acat1</i></th><th><i>acat2</i></th><th><i>Fas</i></th><th><i>srebp-1c</i></th><th><i>gck</i></th></tr></thead><tbody><tr><td>n = 30</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>CHOL</td><td>0.160</td><td>0.140</td><td>0.062</td><td>0.141</td><td>-0.182</td><td>0.167</td><td>0.127</td><td>0.571**</td><td>0.183</td><td>-0.477**</td><td>-0.050</td><td>0.076</td><td>0.087</td><td>0.374*</td><td>0.151</td><td>0.051</td></tr><tr><td>TG</td><td>0.249</td><td>0.254</td><td>0.173</td><td>0.108</td><td>-0.003</td><td>-0.067</td><td>0.134</td><td>0.401*</td><td>0.072</td><td>-0.517**</td><td>-0.136</td><td>0.163</td><td>0.285</td><td>0.480**</td><td>0.337</td><td>0.188</td></tr><tr><td>GLU</td><td>0.157</td><td>0.182</td><td>0.178</td><td>0.134</td><td>0.088</td><td>0.176</td><td>0.153</td><td>0.315</td><td>0.197</td><td>0.219</td><td>-0.032</td><td>-0.026</td><td>-0.144</td><td>0.166</td><td>0.023</td><td>0.180</td></tr></tbody></table>

P-values

* P<0.05

** P<0.01.

Pearson correlation coefficient larger than zero represents a positive correlation, while Pearson correlation coefficient less than zero represents a negative correlation.

https://doi.org/10.1371/journal.pone.0206727.t004
<figure><img src="image_7.png" /><figcaption>Fig 4. Postprandial expression of genes involved in cholesterol esterification including *acat1* (A), *acat2* (B) in the liver of rainbow trout at 3h, 8h, 12h, 16h and 24h after the meal. The NORMA-Gene method was used for the normalization. Values are means (n = 6), with their standard deviations represented by vertical bars. a, b, c Mean values with unlike letters were significantly different (*P*<0.05; One-way analysis of variance, Tukey's test).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0206727.g004

## Correlations between hepatic genes and miRNAs in liver and plasma

Regarding the expression of miRNAs in the liver, miR-122, miR-128 and miR-223 presented negative correlations with genes involved in cholesterol synthesis (*P*<0.05), especially the miR-128 whose expression was correlated with four cholesterol synthetic genes, namely *hmgcs*, *hmgcr*, *cyp51* and *dhcr7*. Additionally, miR-128 expression was also negatively correlated with *acat2*, which is involved in cholesterol esterification (*P*<0.05). The expression of miR-223 was negatively correlated with *abcg8* involved in cholesterol efflux (*P*<0.05) as well as *gck* involved in glycolysis (*P*<0.01) (Table 6).

Regarding miRNAs in the plasma, miR-122 was negatively correlated with *hmgcr* and *cyp51* while miR-128 also showed negative correlations with *hmgcs*, all genes involved in cholesterol synthesis (*P*<0.05). Circulating miR-122 was also negatively correlated (*P*<0.05) with the gene expression of *srebp-2* which is a positive regulator of cholesterol synthesis, *abcg5* which is involved in cholesterol efflux and *acat2*. The plasma level of miR-223 was found to be negatively correlated with many genes involved in hepatic cholesterol efflux, including *ugt1a3*, *abcg5* and *abcg8*, but positively correlated with *abca1* (*P*<0.05). Moreover, miR-233 was also detected to negatively correlate with both *fas* and *gck* (*P*<0.05) (Table 6).

<figure><img src="image_8.png" /><figcaption>Fig 5. Postprandial expression of genes involved in lipogenesis including *fas* (A), *srebp-1c* (B) and *gck* (C) in the liver of rainbow trout at 3h, 8h, 12h, 16h and 24h after the meal. The NORMA-Gene method was used for the normalization. Values are means (n = 6), with their standard deviations represented by vertical bars. a, b, c Mean values with unlike letters were significantly different (*P*<0.05; One-way analysis of variance, Tukey's test).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0206727.g005
<figure><img src="image_9.png" /><figcaption>**Fig 6.** Postprandial expression of miRNAs involved in cholesterol metabolism miR-1(A), miR-33a (B), miR-122(C), miR-128(D), and miR-223(E) in the liver of rainbow trout at 3h, 8h, 12h, 16h and 24h after the meal. Indicated miR expression values were normalized to miR-39 levels. Values are means (n = 6), with their standard deviations represented by vertical bars. a, b Mean values with unlike letters were significantly different (*P*<0.05; One-way analysis of variance, Tukey's test).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0206727.g006

<figure><img src="image_10.png" /><figcaption>**Fig 7.** Postprandial expression of miRNAs involved in cholesterol metabolism miR-1(A), miR-33a (B), miR-122(C), miR-128(D), and miR-223(E) in the plasma of rainbow trout at 3h, 8h, 12h, 16h and 24h after the meal. Expression values were normalized to miR-39 levels. Values are means (n = 6), with their standard deviations represented by vertical bars. a, b, c Mean values with unlike letters were significantly different (*P*<0.05; One-way analysis of variance, Tukey's test).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0206727.g007
<table><caption>Table 5. Correlations of plasmatic parameters, miRNAs expression in the liver, and miRNAs level in the plasma of rainbow trout.</caption><thead><tr><th rowspan="2">Pearson value</th><th colspan="2">n = 30</th><th colspan="2">n = 30</th><th colspan="2">n = 30</th><th colspan="2">n = 30</th><th colspan="2">n = 29</th><th colspan="2">n = 30</th></tr><tr><th>L-miR-1</th><th>L-miR-33a</th><th>L-miR-122</th><th>L-miR-128</th><th>L-miR-223</th><th>P-miR-1</th><th>P-miR-33a</th><th>P-miR-122</th><th>P-miR-128</th><th>P-miR-223</th><th>P-miR-122</th><th>P-miR-223</th></tr></thead><tbody><tr><td>CHOL</td><td>0.229</td><td>0.618**</td><td>-0.045</td><td>0.045</td><td>-0.189</td><td>0.026</td><td>0.006</td><td>0.255</td><td>-0.119</td><td>-0.388*</td><td></td><td></td></tr><tr><td>GLU</td><td>0.285</td><td>0.199</td><td>0.157</td><td>0.126</td><td>-0.256</td><td>-0.056</td><td>0.357</td><td>0.369*</td><td>0.321</td><td>0.239</td><td></td><td></td></tr><tr><td>TG</td><td>0.213</td><td>0.671**</td><td>-0.007</td><td>0.084</td><td>-0.163</td><td>0.126</td><td>0.091</td><td>0.127</td><td>-0.093</td><td>-0.401*</td><td></td><td></td></tr></tbody></table>

n = 30 for miRNAs in plasma except for miR-122 (n = 29).

P-values

* P<0.05

** P<0.01.

Pearson correlation coefficient larger than zero represents a positive correlation, while Pearson correlation coefficient less than zero represents a negative correlation.

https://doi.org/10.1371/journal.pone.0206727.t005

## Correlations between hepatic miRNAs and plasmatic miRNAs

Results showed that hepatic and circulating levels of miR-122 and miR-128 positively correlated with each other (P<0.05) whereas, no correlations were found for miR-1 and miR-223 (P>0.05) (Table 7).

<table><caption>Table 6. Correlations of hepatic gene expression involved in cholesterol metabolism, miRNAs expression in the liver, and miRNAs level in the plasma of rainbow trout.</caption><thead><tr><th rowspan="2">Pearson Value</th><th colspan="5">Cholesterol synthesis</th><th colspan="2">Bile acid synthesis</th><th colspan="2">Cholesterol efflux</th><th colspan="2">Esterification</th><th colspan="2">Lipogenesis</th><th>Glycogenesis</th></tr><tr><th>hmgcs</th><th>hmgcr</th><th>cyp51</th><th>dhcr7</th><th>srebp-2</th><th>cyp7a1</th><th>ugt1a3</th><th>abcg5</th><th>abcg8</th><th>abca1</th><th>lxr</th><th>acat1</th><th>acat2</th><th>fas</th><th>srebp-1c</th><th>gck</th></tr></thead><tbody><tr><td>n = 30 L-miR-1</td><td>-0.134</td><td>-0.191</td><td>-0.196</td><td>-0.249</td><td>-0.180</td><td>0.046</td><td>-0.261</td><td>-0.179</td><td>-0.197</td><td>0.223</td><td>-0.128</td><td>-0.191</td><td>-0.268</td><td>-0.084</td><td>-0.018</td><td>-0.331</td></tr><tr><td>n = 30 L-miR-33a</td><td>-0.094</td><td>-0.075</td><td>-0.123</td><td>-0.146</td><td>-0.270</td><td>-0.112</td><td>-0.050</td><td>0.126</td><td>-0.175</td><td>-0.248</td><td>-0.178</td><td>-0.120</td><td>-0.130</td><td>0.222</td><td>0.051</td><td>-0.061</td></tr><tr><td>n = 30 L-miR-122</td><td>-0.378*</td><td>-0.345</td><td>-0.342</td><td>-0.408*</td><td>-0.301</td><td>-0.086</td><td>-0.192</td><td>-0.202</td><td>-0.265</td><td>0.083</td><td>0.050</td><td>-0.257</td><td>-0.308</td><td>-0.275</td><td>-0.298</td><td>-0.240</td></tr><tr><td>n = 30 L-miR-128</td><td>-0.439*</td><td>-0.385*</td><td>-0.376*</td><td>-0.439*</td><td>-0.306</td><td>-0.157</td><td>-0.201</td><td>-0.215</td><td>-0.339</td><td>0.080</td><td>-0.009</td><td>-0.218</td><td>-0.364*</td><td>-0.250</td><td>-0.249</td><td>-0.322</td></tr><tr><td>n = 30 L-miR-223</td><td>-0.415*</td><td>-0.336</td><td>-0.292</td><td>-0.303</td><td>-0.16</td><td>-0.249</td><td>-0.283</td><td>-0.348</td><td>-0.377*</td><td>0.135</td><td>0.162</td><td>-0.035</td><td>-0.268</td><td>-0.351</td><td>-0.204</td><td>-0.501**</td></tr><tr><td>n = 29 P-miR-1</td><td>0.264</td><td>0.331</td><td>0.287</td><td>0.308</td><td>0.247</td><td>-0.188</td><td>0.033</td><td>-0.113</td><td>0.210</td><td>-0.273</td><td>-0.052</td><td>0.099</td><td>0.173</td><td>0.232</td><td>0.147</td><td>0.321</td></tr><tr><td>n = 29 P-miR-33a</td><td>-0.174</td><td>-0.132</td><td>-0.166</td><td>-0.160</td><td>-0.138</td><td>-0.050</td><td>-0.232</td><td>-0.071</td><td>-0.149</td><td>0.043</td><td>0.208</td><td>-0.117</td><td>-0.212</td><td>-0.071</td><td>-0.070</td><td>0.016</td></tr><tr><td>n = 29 P-miR-122</td><td>-0.328</td><td>-0.372*</td><td>-0.424*</td><td>-0.329</td><td>-0.431*</td><td>0.361</td><td>-0.169</td><td>0.378*</td><td>-0.059</td><td>0.054</td><td>0.284</td><td>-0.076</td><td>-0.376*</td><td>-0.261</td><td>-0.190</td><td>-0.100</td></tr><tr><td>n = 29 P-miR-128</td><td>-0.438*</td><td>-0.421*</td><td>-0.408*</td><td>-0.335</td><td>-0.310</td><td>-0.136</td><td>-0.331</td><td>-0.239</td><td>-0.355</td><td>0.310</td><td>0.188</td><td>-0.258</td><td>-0.459*</td><td>-0.326</td><td>-0.289</td><td>-0.327</td></tr><tr><td>n = 29 P-miR-223</td><td>-0.334</td><td>-0.305</td><td>-0.252</td><td>-0.266</td><td>-0.109</td><td>-0.292</td><td>-0.465*</td><td>-0.408*</td><td>-0.414*</td><td>0.456*</td><td>0.120</td><td>-0.173</td><td>-0.287</td><td>-0.383*</td><td>-0.146</td><td>-0.452*</td></tr></tbody></table>

n = 30 for miRNAs in liver, n = 29 for miRNAs in plasma.

p-values

* P<0.05

** P<0.01.

Pearson correlation coefficient larger than zero represents a positive correlation, while Pearson correlation coefficient less than zero represents a negative correlation.

https://doi.org/10.1371/journal.pone.0206727.t006
<table><caption>Table 7. Correlations between miRNAs in the liver and plasma of rainbow trout.</caption><thead><tr><th rowspan="2">Pearson value</th><th rowspan="2"></th><th colspan="5">Liver miR n = 30</th><th colspan="5">Plasma miR n = 29</th></tr><tr><th>L-miR-1</th><th>L-miR-33a</th><th>L-miR-122</th><th>L-miR-128</th><th>L-miR-223</th><th>P-miR-1</th><th>P-miR-33a</th><th>P-miR-122</th><th>P-miR-128</th><th>P-miR-223</th></tr></thead><tbody><tr><td rowspan="4">Liver miR n = 30</td><td>L-miR-1</td><td></td><td>0.596**</td><td>0.455*</td><td>0.597**</td><td>0.304</td><td>-0.074</td><td>0.319</td><td>0.275</td><td>0.397*</td><td>0.191</td></tr><tr><td>L-miR-33a</td><td></td><td></td><td>0.508**</td><td>0.644**</td><td>0.245</td><td>0.106</td><td>0.355</td><td>0.171</td><td>0.336</td><td>0.022</td></tr><tr><td>L-miR-122</td><td></td><td></td><td></td><td>0.890**</td><td>0.588**</td><td>-0.065</td><td>0.526**</td><td>0.182</td><td>0.595**</td><td>0.431*</td></tr><tr><td>L-miR-128</td><td></td><td></td><td></td><td></td><td>0.713**</td><td>-0.102</td><td>0.372*</td><td>0.187</td><td>0.578**</td><td>0.419*</td></tr><tr><td rowspan="5">Plasma miR n = 29</td><td>L-miR-223</td><td></td><td></td><td></td><td></td><td></td><td>0.066</td><td>0.246</td><td>-0.041</td><td>0.390</td><td>0.548**</td></tr><tr><td>P-miR-1</td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.367</td><td>-0.225</td><td>-0.065</td><td>0.047</td></tr><tr><td>P-miR-33a</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.355</td><td>0.708**</td><td>0.622**</td></tr><tr><td>P-miR-122</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.459*</td><td>0.234</td></tr><tr><td>P-miR-128</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.780**</td></tr><tr><td>P-miR-223</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

n = 30 for miRNAs in liver, n = 29 for miRNAs in plasma.

P-values

* P<0.05

** P<0.01.

Pearson correlation coefficient larger than zero represents a positive correlation, while Pearson correlation coefficient less than zero represents a negative correlation.

https://doi.org/10.1371/journal.pone.0206727.t007

Discussion

In the present study, rainbow trout was used as a research model in order to 1) understand the postprandial kinetic of parameters related to cholesterol metabolism, including miRNAs and 2) evaluate if circulating miRNAs could reflect the level of expression of genes in the liver and therefore be used as non-invasive biomarkers. For this purpose, we analyzed parameters at the phenotypic level (cholesterol, triglycerides and glucose in plasma) and transcriptional level (gene expressions in liver and miRNAs expressions in liver and plasma) at several time points after the last meal. Furthermore, the correlations between these parameters were investigated to assess the potential of miRNAs to be used as biomarkers in aquaculture regarding the development of alternative diets devoid of fishmeal and fish oil.

In this study, postprandial levels of cholesterol, triglycerides and glucose in plasma were only slightly affected by the meal. Originated from Acetyl-CoA, cholesterol is synthesized by about 20 enzymes and regulated by the SREBP pathway [68]. Following a meal, the hepatic expression kinetic of genes encoding four important enzymes responsible for cholesterol synthesis (*hmgcs*, *hmgcr*, *cyp51* and *dhcr7*), and their transcription factor *srebp-2* were quite similar. After a relatively stable period, the expression of these genes increased 16h after the meal before to come to basal level 24h after the meal, reflecting synchronous effects on cholesterol synthesis. Consistently with these results, significant positive correlations were observed between these genes, confirming their close cooperation during cholesterol synthesis.

*Abcg5*, *abcg8* and *ugt1a3* are genes encoding proteins involved in hepatic cholesterol efflux and elimination. ABCG5 and ABCG8 form an heterodimer that acts in the liver and intestine as a cholesterol transporter to excrete cholesterol into the bile [69], while UGT1A3 catalyzes its subsequent glucuronidation during bile acid synthesis [70]. The expression of these three genes showed the same postprandial pattern and their expression was also very similar to that of cholesterol synthetic genes. This suggests a close functional relation between ABCG5 and ABCG8 at the transcriptional level in rainbow trout. However, other genes involved in cholesterol efflux, such as, *cyp7a1*, *abca1* and *lxrα* were not significantly affected by feeding. Lxrα has been shown to act as a cholesterol sensor to regulate the expression of *cyp7a1*, *abcg5*, *abcg8* and
*abca1* [71]. Accordingly, in this study, we found a significant correlation between *lxrα* and *cyp7a1* and *abcg5* in rainbow trout.

ACAT1 and ACAT2 are enzymes that catalyze the esterification of cholesterol to cholesterol ester for storage or transportation [72]. Although a peak of expression was recorded for these two genes 16h after the meal, *acat1* and *acat2* expressions did not correlate with each other. Nevertheless, they presented significant correlation with the expression of genes involved in cholesterol biosynthesis and *srebp-2*. This may be attributed to the fact that *acat2* is one of *srebp-2* target genes as previously reported in mouse [73] and plays important roles in cholesterol storage and transportation.

Among lipid metabolism, SREBP-1c is the major regulator of several genes related to fatty acid synthesis amongst which *fas* encodes for the key enzyme of fatty acid synthesis [74,75]. This is supported by the significant positive correlation found in the present study between *fas* and *srebp-1c* gene expression. Meanwhile, the positive correlation between these two genes and the four main genes involved in cholesterol synthesis also suggested a close relationship between fatty acid synthesis and cholesterol synthesis. GCK is known to play a major role in controlling glucose utilization in glycolytic pathways [76] and is also considered as the first enzyme of lipogenesis. The positive correlation between *gck* and cholesterol synthetic and efflux genes may suggest that cholesterol synthesis and efflux are processes requiring ATP which could be supplied by the glycolytic reaction catalyzed by GCK.

At posttranscriptional level, several miRNAs involved in cholesterol metabolism were detected both in liver and plasma and showed postprandial variations. Postprandial changes of miRNAs expression have already been demonstrated in trout liver [47] and skeletal muscle of Chinese perch (*Siniperca chuatsi*) [77]. The present study further indicates that circulating miRNAs levels also vary following trout feeding.

According to previous studies, the functional mechanisms of several miRNAs involved in cholesterol regulation are as follows: miR-1 may promote cellular cholesterol through LXRα suppression [59], miR-33a represses cholesterol efflux by targeting ABCG5, ABCG8, CYP7A1 and ABCA1 [60,78,79], miR-122 may increase some genes in cholesterol synthesis, though this mechanism of activation is still unknown [80], miR-128 increases cellular cholesterol by ABCA1 suppression and SREBP-2 promotion [63] and miR-223 reduces cellular cholesterol and promotes cholesterol efflux by repressing directly *hmgcs* and indirectly *abca1* and *cyp7a1*, respectively [64].

The hepatic expression of miR-1, miR-122 and miR-128 showed very similar evolution after the meal, reaching their highest and lowest expressions 12h and 16h after the meal, respectively. Compared with most genes analyzed in the present study, which reached their highest expression 16h after the meal, the expression of miR-1, miR-122 and miR-128 increased earlier. Although further investigations are needed to highlight a cause to effect relationship, our results suggest that these miRNAs may be involved in the regulation of these aforementioned genes.

The correlation found between miR-122 and cholesterol synthetic genes is in agreement with the effect of miR-122 on cholesterol biosynthesis determined in a previous study using antagonists to inhibit miR-122 *in vivo* [80]. In mouse, it was also shown that miR-223 inhibited the cholesterol synthesis by the direct repression of *hmgcs* [64]. This is in agreement with our results showing that miR-223 negatively correlates with genes involved in cholesterol synthesis. It was also reported that the expression of *srebp-1*, *srebp-2*, *lxrα* and *abca1* was affected by miR-128 in mouse and human cell lines [63]. However, in our study, we did not observe any correlation between miR-128 and these genes in rainbow trout. The negative correlation between miR-122, miR-223 and miR-128 and genes involved in hepatic cholesterol metabolism suggest that these genes may be potential direct targets of miR-122, miR-223 and miR-
128. Further analysis of the complementarity between 3' untranslated region of the genes and seed sequences of the miRNAs are needed to establish this relation. Although none of the genes evaluated in the present study correlate with hepatic and plasma levels of miR-33a, miR-33a is the only miRNA that was strongly and positively correlated with plasma levels of cholesterol and triglycerides, suggesting a role of miR-33a in cholesterol accumulation in fish as previously shown in mammals [60,78].

For the first time, we demonstrate the presence of miRNAs (miR-1, miR-33a, miR-122, miR-128 and miR-223) in the plasma of rainbow trout. Based on mammalian studies, miRNAs present in the plasma may have several origins: the passive release from dead cells, an active secretion *via* microvesicles or an active secretion with RNA-binding proteins [19]. The potential roles of secreted miRNAs have been reported to be involved in the regulation of various physiological and pathological processes in mammals. For example, Mittelbrunn M. *et al*. proposed that miR-335 generated from T cell are transferred to antigen presenting cells and regulate its cellular target genes during immune response [81], suggesting that secreted miRNAs can be used as biomarkers to reflect target gene expression in cells. Despite the biological function and secretory mechanism of extracellular miRNAs are still unclear and remain controversial in mammals, the potential of miRNAs as biomarkers for cancer diagnosis and prognosis has already been verified. For example, miR-29a and miR-92a levels in plasma are significantly higher in patients with advanced-stage colorectal carcinoma than healthy controls [82]. Therefore, we further studied miRNAs in plasma and analyzed the correlations between circulating miRNAs and their hepatic counterparts as well as correlations between circulating miRNAs and genes involved in cholesterol metabolism in order to evaluate the potential of circulating miRNAs to be used as non-invasive biomarkers in rainbow trout.

Interestingly, we found a positive correlation between hepatic expression and plasma level of miR-128 and miR-223 suggesting that circulating miRNAs may be indicative of the expression of their hepatic counterparts in trout. As it is known that 1) circulating miRNAs may be released from either byproducts of cellular activity or selective miRNA export system [34], and 2) liver plays a central role in cholesterol homeostasis [83], the positive correlations that we have found for miR-128 and miR-223 levels between liver and plasma may reflect their active secretion from hepatocyte into plasma for posttranscriptional regulation of cholesterol metabolism in extrahepatic tissues. However, as miRNAs are widely expressed among tissues, their expression should be evaluated in other tissues to validate that liver exhibits the major contribution to circulating levels of miR-128 and miR-223.

As their hepatic counterparts, miR-122, miR-128 and miR-223 in plasma exhibit significant correlations with the expression of genes involved in cholesterol synthesis, efflux or esterification. The correlation between gck and miR-223 was also confirmed at circulating level. Once again, these results highlight the importance of miR-122 and miR-128 in cholesterol synthesis regulation and miR-223 in cholesterol efflux and glycolysis regulation. Circulating miR-223 was the only plasmatic miRNA that negatively correlates with plasma cholesterol level. This correlation was therefore consistent with the negative effect assigned to miR-223 on cholesterol homeostasis in mammals [64]. Altogether, our results strongly suggest that circulating miR-122, miR-128 and miR-223 could represent useful biomarkers of cholesterol and glucose metabolism in rainbow trout.

In conclusion, the present study demonstrates the presence of miRNAs in rainbow trout plasma and showed that hepatic expression of miRNAs and circulating miRNAs level are both subjected to postprandial regulations which could be attributed to nutritional aspects but also to circadian rhythm since analyses were performed over 24 hours. This study also indicates that correlations are found between hepatic and circulating levels of miRNAs as shown for miR-128 and miR-223 and between hepatic gene expression and liver or plasma miRNAs, as
observed for miR-128 and genes involved in cholesterol synthesis. These results support the evidence that circulating miRNAs can reflect the level of expression of their tissue counterparts as well as the expression of genes in rainbow trout and highlight circulating miRNAs as promising non-invasive biomarkers in trout. However, correlations between circulating miRNAs and hepatic gene expression should not presume any direct or indirect regulatory relationship. Therefore, bioinformatics prediction of miRNA/target relationship as well as the implementation of functional studies based on miRNAs inhibitors or mimics or the development of 3'UTR luciferase reporter assays would be necessary to clearly demonstrate the regulatory role of these miRNAs on cholesterol metabolism. Further studies are now also necessary to better understand the mechanisms underlying miRNAs blood secretion and their subsequent molecular and physiological functions in fish. In particular, long term dietary trials with plant-based diets devoid of cholesterol should be conducted in order to improve our knowledges on the role of miRNAs in fish.

**Supporting information**

**S1 File. Plasma metabolites, miRNA expression and gene expression data.** Plasma metabolites excel sheet: postprandial levels of cholesterol (Chol), glucose and triglycerides (TG) in plasma of rainbow trout at 3h, 8h, 12h, 16h and 24h after the meal (n = 6 at each time point). Mir expression excel sheet: data of hepatic expression (Lmir) and plasma abundance (Pmir) of miR-1, miR-33, miR-122, miR-128 and miR-223 measured 3h, 8h, 12h, 16h and 24h after the meal (n = 6 at each time point). miRNA data are normalized with miR-39 used as exogenous spike miRNA. Gene expression excel sheet: data of hepatic expression of genes related to cholesterol synthesis (*hmgcr*, *hmgcs*, *cyp51*, *dhcr7* and *srebp2*), cholesterol efflux (*cyp7a1*, *ugt1a3*, *abcg5*, *abcg8*, *abca1* and *lxr*), cholesterol esterification (*acat1* and *acat2*) and lipogenesis (*fas*, *gck* and *srebp1*) measured 3h, 8h, 12h, 16h and 24h after the meal (n = 6 at each time point). miRNA data are normalized with NORMA-Gene method. (XLSX)

**Acknowledgments**

The authors thank the technical staff (F. Terrier, A. Lanuque and F. Sandres) of the INRA experimental fish farm facilities of Donzacq for fish rearing and thank A. Herman for technical assistance in the laboratory.

**Author Contributions**

**Conceptualization:** Tengfei Zhu, Geneviève Corraze, Sandrine Skiba-Cassy.

**Data curation:** Tengfei Zhu, Elisabeth Plagnes-Juan.

**Formal analysis:** Tengfei Zhu, Elisabeth Plagnes-Juan.

**Investigation:** Tengfei Zhu, Elisabeth Plagnes-Juan.

**Methodology:** Tengfei Zhu, Elisabeth Plagnes-Juan.

**Supervision:** Geneviève Corraze, Sandrine Skiba-Cassy.

**Validation:** Geneviève Corraze, Sandrine Skiba-Cassy.

**Writing – original draft:** Tengfei Zhu.

**Writing – review & editing:** Tengfei Zhu, Geneviève Corraze, Elisabeth Plagnes-Juan, Sandrine Skiba-Cassy.
References

1. Iwakawa H, Tomari Y. The functions of microRNAs: mRNA decay and translational repression. Trends Cell Biol. 2015; 25: 651-665. https://doi.org/10.1016/j.tcb.2015.07.011 PMID: 26437588

2. Axtell MJ, Bartel DP. Antiquity of microRNAs and their targets in land plants. Plant Cell. 2005; 17: 1658-1673. https://doi.org/10.1105/tpc.105.032185 PMID: 15849273

3. Ibáñez-Ventoso Carolina MV, Driscoll M. Sequence relationships among C. elegans, D. melanogaster and human microRNAs highlight the extensive conservation of microRNAs in biology. PLoS One. 2008; 3: e2818. https://doi.org/10.1371/journal.pone.0002818

4. Ambros V. The functions of animal microRNAs. Nature. 2004; 431: 350-5. https://doi.org/10.1038/nature02871 PMID: 15372042

5. Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004; 116: 281-297. https://doi.org/10.1016/S0092-8674(04)00045-5 PMID: 14744438

6. Bushati N, Cohen SM. microRNA Functions. Annu Rev Cell Dev Biol. 2007; 23: 175-205. https://doi.org/10.1146/annurev.cellbio.23.090506.123406 PMID: 17506695

7. Drakaki A, Iliopoulos D. MicroRNA gene networks in oncogenesis. Curr Genomics. 2009; 10: 35-41. https://doi.org/10.2174/138920209787581299 PMID: 19721809

8. Sunkar R, Li Y, Jagadeeswaran G. Functions of microRNAs in plant stress responses. Trends Plant Sci. Elsevier Ltd; 2012; 17: 196-203. https://doi.org/10.1016/j.tplants.2012.01.010 PMID: 22365280

9. Shenoy A, Blelloch RH. Regulation of microRNA function in somatic stem cell proliferation and differentiation. Nat Publ Gr. Nature Publishing Group; 2014; 15: 565-576. https://doi.org/10.1038/nrm3854 PMID: 25118717

10. Winter J, Jung S, Keller S, Gregory RI, Diederichs S. Many roads to maturity: microRNA biogenesis pathways and their regulation. Nat Cell Biol. 2009; 11: 228-234. https://doi.org/10.1038/ncb0309-228 PMID: 19255566

11. Wienholds E, Plasterk RHA. MicroRNA function in animal development. FEBS Lett. 2005; 579: 5911-5922. https://doi.org/10.1016/j.febslet.2005.07.070 PMID: 16111679

12. Chim SSC, Shing TKF, Hung ECW, Leung T, Lau T, Chiu RWK, et al. Detection and characterization of placental microRNAs in maternal plasma. Clin Chem. 2008; 54: 482-490. https://doi.org/10.1373/clinchem.2007.097972 PMID: 18218722

13. Lawrie CH, Gal S, Dunlop M, Pushkaran B, Liggins P, Pulford K, et al. Detection of elevated levels of tumour-associated microRNAs in serum of patients with diffuse large B-cell lymphoma. Br J Haematol. 2008; 141: 672-675. https://doi.org/10.1111/j.1365-2141.2008.07077.x PMID: 18318758

14. Weber JA, Baxter DH, Zhang S, Huang DY, Huang KH, Lee MJ, et al. The microRNA spectrum in 12 body fluids. Clin Chem. 2010; 56: 1733-1741. https://doi.org/10.1373/clinchem.2010.147405 PMID: 20847327

15. Etheridge A, Lee I, Hood L, Galas D, Wang K. Extracellular microRNA: a new source of biomarkers. Mutat Res Mol Mech Mutagen. Elsevier; 2011; 717: 85-90. https://doi.org/10.1016/j.mrfmmm.2011.03.004 PMID: 21402084

16. Hanke M, Hoefig K, Merz H, Feller AC, Kausch I, Jocham D, et al. A robust methodology to study urine microRNA as tumor marker: microRNA-126 and microRNA-182 are related to urinary bladder cancer. Urol Oncol Semin Orig Investig. Elsevier Inc.; 2010; 28: 655-661. https://doi.org/10.1016/j.urolonc.2009.01.027 PMID: 19375957

17. Park NJ, Zhou H, Elashoff D, Henson BS, Kastratovic DA, Abemayor E, et al. Salivary microRNA: discovery, characterization, and clinical utility for oral cancer detection. Clin Cancer Res. 2009; 15: 5473-5477. https://doi.org/10.1158/1078-0432.CCR-09-0736 PMID: 19706812

18. Kosaka N, Izumi H, Sekine K, Ochiya T. microRNA as a new immune-regulatory agent in breast milk. Silence. 2010; 1: 7. https://doi.org/10.1186/1758-907X-1-7 PMID: 20226005

19. Chen X, Liang H, Zhang J, Zen K, Zhang CY. Secreted microRNAs: A new form of intercellular communication. Trends Cell Biol. Elsevier Ltd; 2012; 22: 125-132. https://doi.org/10.1016/j.tcb.2011.12.001 PMID: 22260888

20. Lagos-Quintana M, Rauhut R, Yalcin A, Meyer J, Lendeckel W, Tuschl T. Identification of tissue-specific microRNAs from mouse. Curr Biol. 2002; 12: 735-739. https://doi.org/10.1016/S0960-9822(02)00809-6 PMID: 12007417

21. Sempere LF, Freemantle S, Pitha-rowe I, Moss E, Dmitrovsky E, Ambros V. Expression profiling of mammalian microRNAs uncovers a subset of brain-expressed microRNAs with possible roles in murine and human neuronal differentiation. Genome Biol. 2004; 5: R13. https://doi.org/10.1186/gb-2004-5-3-r13 PMID: 15003116
22. Watanabe T, Takeda A, Mise K, Okuno T, Suzuki T. Stage-specific expression of microRNAs during Xenopus development. FEBS Lett. 2005; 579: 318–324. https://doi.org/10.1016/j.febslet.2004.11.067 PMID: 15642338

23. Cortez MA, Bueso-ramos C, Ferdin J, Lopez-berestein G, Anil K, Calin GA. MicroRNAs in body fluids— the mix of hormones and biomarkers. Nat Rev Clin Oncol. 2011; 8: 467–477. https://doi.org/10.1038/nrclinonc.2011.76 PMID: 21647195

24. Mitchell PS, Parkin RK, Kroh EM, Fritz BR, Wyman SK, Pogosova-agadjanyan EL, et al. Circulating microRNAs as stable blood-based markers for cancer detection. Proc Natl Acad Sci. 2008; 105: 10513–10518. https://doi.org/10.1073/pnas.0804549105 PMID: 18663219

25. Ai J, Zhang R, Li Y, Pu J, Lu Y, Jiao J, et al. Circulating microRNA-1 as a potential novel biomarker for acute myocardial infarction. Biochem Biophys Res Commun. Elsevier; 2010; 391: 73–77. https://doi.org/10.1016/j.bbrc.2009.11.005 PMID: 19896465

26. Cheng Y, Tan N, Yang J, Liu X, Cao X, He P, et al. A translational study of circulating cell-free microRNA-1 in acute myocardial infarction. Clin Sci. 2010; 119: 87–95. https://doi.org/10.1042/CS20090645 PMID: 20218970

27. Turchinovich A, Weiz L, Langheinz A, Burwinkel B. Characterization of extracellular circulating microRNA. Nucleic Acids Res. 2011; 39: 7223–7233. https://doi.org/10.1093/nar/gkr254 PMID: 21609964

28. Arroyo JD, Chevillet JR, Kroh EM, Ruf IK, Pritchard CC, Gibson DF, et al. Argonaute2 complexes carry a population of circulating microRNAs independent of vesicles in human plasma. Proc Natl Acad Sci. 2011; 108: 5003–5008. https://doi.org/10.1073/pnas.1019055108 PMID: 21383194

29. Zernecke A, Bidzhekov K, Noels H, Shagdarsuren E, Denecke B, Hristov M, et al. Delivery of microRNA-126 by apoptotic bodies induces CXCL12-dependent vascular protection. Sci Signal. 2009; 2: ra81–ra81. https://doi.org/10.1126/scisignal.2000610 PMID: 19996457

30. Valadi H, Ekström K, Bossios A, Sjöstrand M, Lee JJ, Lötvall JO. Exosome-mediated transfer of mRNAs and microRNAs is a novel mechanism of genetic exchange between cells. Nat Cell Biol. 2007; 9: 654. https://doi.org/10.1038/ncb1596 PMID: 17486113

31. Muralidharan-Chari V, Clancy JW, Sedgwick A, D'Souza-Schorey C. Microvesicles: mediators of extra-cellular communication during cancer progression. J Cell Sci. 2010; 123: 1603–1611. https://doi.org/10.1242/jcs.064386 PMID: 20445011

32. Vickers KC, Palmisano BT, Shoucri BM, Shamburek RD, Remaley AT. MicroRNAs are transported in plasma and delivered to recipient cells by high-density lipoproteins. Nat Cell Biol. Nature Publishing Group; 2011; 13: 423–433. https://doi.org/10.1038/ncb2210 PMID: 21423178

33. Wagner J, Riwanto M, Besler C, Knau A, Fichtlscherer S, Röxe T, et al. Characterization of levels and cellular transfer of circulating lipoprotein-bound microRNAs. Arterioscler Thromb Vasc Biol. 2013; 33: 1392–1400. https://doi.org/10.1161/ATVBAHA.112.300741 PMID: 23559634

34. Turchinovich A, Weiz L, Burwinkel B. Extracellular miRNAs: The mystery of their origin and function. Trends Biochem Sci. Elsevier Ltd; 2012; 37: 460–465. https://doi.org/10.1016/j.tibs.2012.08.003 PMID: 22944280

35. Wienholds E, Kloosterman WP, Miska E, Alvarez-Saavedra E, Berezikov E, Brujin E de, et al. MicroRNA expression in zebrafish embryonic development. Science. 2005; 309: 310–311. https://doi.org/10.1126/science.1114519 PMID: 15919954

36. Bizuayehu TT, Lanes CFC, Furmanek T, Karlsen BO, Fernandes JMO, Johansen SD, et al. Differential expression patterns of conserved miRNAs and isomiRs during Atlantic halibut development. BMC Genomics. BioMed Central Ltd; 2012; 13: 11. https://doi.org/10.1186/1471-2164-13-11 PMID: 22233483

37. Mishima Y. Widespread roles of microRNAs during zebrafish development and beyond. Dev Growth Differ. 2012; 54: 55–65. https://doi.org/10.1111/j.1440-169X.2011.01306.x PMID: 22150108

38. Lim LP, Glasner ME, Yekta S, Burge CB, Bartel DP. Vertebrate microRNA genes. Science. 2003; 299: 1540. https://doi.org/10.1126/science.1080372 PMID: 12624257

39. Giraldez AJ, Cinalli RM, Glasner ME, Enright AJ, Thomson JM, Baskerville S, et al. MicroRNAs regulate brain morphogenesis in zebrafish. Science. 2005; 308: 833–838. https://doi.org/10.1126/science.1109020 PMID: 15774722

40. Chen PY, Manninga H, Slanchev K, Chien M, Russo JJ, Ju J, et al. The developmental miRNA profiles of zebrafish as determined by small RNA cloning. Genes Dev. 2005; 19: 1288–1293. https://doi.org/10.1101/gad.1310605 PMID: 15937218

41. Thorgaard GH, Bailey GS, Williams D, Buhler DR, Kaattari SL, Ristow SS, et al. Status and opportuni-ties for genomics research with rainbow trout. Comp Biochem Physiol Part B Biochem Mol Biol. 2002; 133: 609–646. https://doi.org/10.1016/S1096-4959(02)00167-7
42. Ramachandra RK, Salem M, Gahr S, Rexroad CE, Yao J. Cloning and characterization of microRNAs from rainbow trout (Oncorhynchus mykiss): their expression during early embryonic development. BMC Dev Biol. 2008; 8: 41. https://doi.org/10.1186/1471-213X-8-41 PMID: 18412968

43. Salem M, Xiao C, Womack J, Rexroad CE, Yao J. A microRNA repertoire for functional genome research in rainbow trout (Oncorhynchus mykiss). Mar Biotechnol. 2010; 12: 410–429. https://doi.org/10.1007/s10126-009-9232-z PMID: 19816740

44. Ma H, Hostuttler M, Wei H, Rexroad CE, Yao J. Characterization of the rainbow trout egg microRNA transcriptome. PLoS One. 2012; 7: e39649. https://doi.org/10.1371/journal.pone.0039649 PMID: 22761856

45. Juanchich A, Bardou P, Rué O, Gabillard J-C, Gaspin C, Bobe J, et al. Characterization of an extensive rainbow trout miRNA transcriptome by next generation sequencing. BMC Genomics. BMC Genomics; 2016; 17: 164. https://doi.org/10.1186/s12864-016-2505-9 PMID: 26931235

46. Mennigen JA, Skiba-Cassy S, Panserat S. Ontogenetic 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. https://doi.org/10.1242/jeb.082248 PMID: 23348939

47. Mennigen JA, Panserat S, Larquier M, Plagnes-Juan E, Medale F, Seiliez I, et al. Postprandial regulation of hepatic microRNAs predicted to target the insulin pathway in rainbow trout. PLoS One. 2012; 7: e38604. https://doi.org/10.1371/journal.pone.0038604 PMID: 22701681

48. Mennigen JA, Plagnes-Juan E, Figueredo-Silva CA, Seiliez I, Panserat S, Skiba-Cassy S. Acute endocrine and nutritional co-regulation of the hepatic omy-miRNA-122b and the lipogenic gene fas in rainbow trout, Oncorhynchus mykiss. Comp Biochem Physiol—B Biochem Mol Biol. Elsevier B.V.; 2014; 169: 16–24. https://doi.org/10.1016/j.cbpb.2013.12.002 PMID: 24333236

49. Mennigen JA, Martyniuk CJ, Seiliez I, Panserat S, Skiba-Cassy S. Metabolic consequences of microRNA-122 inhibition in rainbow trout, Oncorhynchus mykiss. BMC Genomics. 2014; 15: 70. https://doi.org/10.1186/1471-2164-15-70 PMID: 24467738

50. Naylor RL, Hardy RW, Bureau DP, Chiu A, Elliott M, Farrell AP, et al. Feeding aquaculture in an era of finite resources. Proc Natl Acad Sci. National Acad Sciences; 2009; 106: 15103–15110. https://doi.org/10.1073/pnas.0905235106 PMID: 19805247

51. Friedman M. Nutritional value of proteins from different food sources. A review. J Agric Food Chem. ACS Publications; 1996; 44: 6–29. https://doi.org/10.1021/jf9400167

52. Turchini GM, Ng W-K, Tocher DR. Fish oil replacement and alternative lipid sources in aquaculture feeds. Boca Raton, FL, USA: CRC Press, Taylor & Francis group; 2010.

53. Tocher DR, Bendiksen EÅ, Campbell PJ, Bell JG. The role of phospholipids in nutrition and metabolism of teleost fish. Aquaculture. 2008; 280: 21–34. https://doi.org/10.1016/j.aquaculture.2008.04.034

54. Leaver MJ, Villeneuve LA, Obach A, Jensen L, Bron JE, Tocher DR, et al. Functional genomics reveals increases in cholesterol biosynthetic genes and highly unsaturated fatty acid biosynthesis after dietary substitution of fish oil with vegetable oils in Atlantic salmon (Salmo salar). BMC Genomics. 2008; 9: 299. https://doi.org/10.1186/1471-2164-9-299 PMID: 18577222

55. Kortner TM, Gu J, Krogdahl Å, Bakke AM. Transcriptional regulation of cholesterol and bile acid metabolism after dietary soyabean meal treatment in Atlantic salmon (Salmo salar L.). Br J Nutr. 2013; 109: 593–604. https://doi.org/10.1017/S0007114512002024 PMID: 22647297

56. Zhu T, Corraze G, Plagnes-Juan E, Quillet E, Dupont-Nivet M, Sandrine S-C. Regulation of genes related to cholesterol metabolism in rainbow trout (Oncorhynchus mykiss) fed a plant-based diet. Am J Physiol Integr Comp Physiol. 2017; 314: R58–R70. https://doi.org/10.1152/ajpregu.00179.2017 PMID: 28931545

57. Kortner TM, Björkhem I, Krasnov A, Timmerhaus G, Krogdahl Å. Dietary cholesterol supplementation to a plant-based diet suppresses the complete pathway of cholesterol synthesis and induces bile acid production in Atlantic salmon (Salmo salar L.). Br J Nutr. 2014; 111: 2089–2103. https://doi.org/10.1017/S0007114514000373 PMID: 24635969

58. Zhu T, Ai Q, Mai K, Xu W, Zhou H, Liufu Z. Feed intake, growth performance and cholesterol metabolism in juvenile turbot (Scophthalmus maximus L.) fed defatted fish meal diets with graded levels of cholesterol. Aquaculture; 2014; 428: 290–296. https://doi.org/10.1016/j.aquaculture.2014.03.027

59. Zhong D, Huang G, Zhang Y, Zeng Y, Xu Z, Zhao Y, et al. MicroRNA-1 and microRNA-206 suppress LXRα-induced lipogenesis in hepatocytes. Cell Signal. 2013; 25: 1429–1437. https://doi.org/10.1016/j cellsig.2013.03.003 PMID: 23499676

60. Gerin I, Clerbaux LA, Haumont O, Lanthier N, Das AK, Burant CF, et al. Expression of miR-33 from an SREBP2 intron inhibits cholesterol export and fatty acid oxidation. J Biol Chem. 2010; 285: 33652–33661. https://doi.org/10.1074/jbc.M110.152090 PMID: 20732877
61. Najafi-Shoushtari SH, Kristo F, Li Y, Shioda T, Cohen DE, Gerszten RE, et al. MicroRNA-33 and the SREBP host genes cooperate to control cholesterol homeostasis. Science. 2010; 328: 1566–1569.
https://doi.org/10.1126/science.1189123 PMID: 20466882

62. Esau C, Davis S, Murray SF, Yu XX, Pandey SK, Pear M, et al. miR-122 regulation of lipid metabolism revealed by in vivo antisense targeting. Cell Metab. 2006; 3: 87–98. https://doi.org/10.1016/j.cmet.2006.01.005 PMID: 16459310

63. Adlakha YK, Khanna S, Singh R, Singh VP, Agrawal A, Saini N. Pro-apoptotic miRNA-128-2 modulates ABCA1, ABCG1 and RXRα expression and cholesterol homeostasis. Cell Death Dis. 2013; 4: e780.
https://doi.org/10.1038/cddis.2013.301 PMID: 23990020

64. Vickers KC, Landstreet SR, Levin MG, Shoucri BM, Toth CL, Taylor RC, et al. MicroRNA-223 coordinates cholesterol homeostasis. Proc Natl Acad Sci. 2014; 111: 14518–14523. https://doi.org/10.1073/pnas.1215767111 PMID: 25246565

65. Wang L, Jia X-J, Jiang H-J, Du Y, Yang F, Si S-Y, et al. MicroRNAs 185, 96, and 223 repress selective high-density lipoprotein cholesterol uptake through posttranscriptional inhibition. Mol Cell Biol. 2013; 33: 1956–1964. https://doi.org/10.1128/MCB.01580-12 PMID: 23459944

66. Heckmann L, Sørensen PB, Krogh PH, Sørensen JG. NORMA-Gene: A simple and robust method for qPCR normalization based on target gene data. BMC Bioinformatics. 2011; 12: 250. https://doi.org/10.1186/1471-2105-12-250 PMID: 21693017

67. Fox J, Bouchet-Valat M. Rcmdr: R Commander. R package version 2.3–0 [Internet]. 2016. Available: http://www.r-project.org

68. Goldstein JL, DeBose-Boyd RA, Brown MS. Protein sensors for membrane sterols. Cell. 2006; 124: 35–46. https://doi.org/10.1016/j.cell.2005.12.022 PMID: 16413480

69. Wang J, Mitsche MA, Lütjohann D, Cohen JC, Xie X-S, Hobbs HH. Relative roles of ABCG5/ABCG8 in liver and intestine. J Lipid Res. 2015; 56: 319–330. https://doi.org/10.1194/jlr.M054544 PMID: 25378657

70. Barbier O, Trottier J, Kaeding J, Caron P, Verreault M. Lipid-activated transcription factors control bile acid glucuronidation. Mol Cell Biochem. 2009; 326: 3–8. https://doi.org/10.1007/s11010-008-0001-5 PMID: 19130183

71. Lu TT, Repa JJ, Mangelsdorf DJ. Orphan nuclear receptors as eLiXiRs and FiXeRs of sterol metabolism. J Biol Chem. 2001; 276: 37735–37738. https://doi.org/10.1074/jbc.R100035200 PMID: 11459853

72. Chang T, Li B, Chang CCY, Urano Y. Acyl-coenzyme A: cholesterol acyltransferases. Am J Physiol Metab. 2009; 297: E1–E9. https://doi.org/10.1152/ajpendo.90926.2008 PMID: 19141679

73. Horton JD, Shah NA, Warrington JA, Anderson NN, Park SW, Brown MS, et al. Combined analysis of oligonucleotide microarray data from transgenic and knockout mice identifies direct SREBP target genes. Proc Natl Acad Sci. 2003; 100: 12027–12032. https://doi.org/10.1073/pnas.1534923100 PMID: 14512514

74. Horton JD, Goldstein JL, Brown MS. SREBPs: activators of the complete program of cholesterol and fatty acid synthesis in the liver. J Clin Invest. 2002; 109: 1125–1131. https://doi.org/10.1172/JCI15593 PMID: 11994399

75. Latasa MJ, Moon YS, Kim KH, Sul HS. Nutritional regulation of the fatty acid synthase promoter in vivo: sterol regulatory element binding protein functions through an upstream region containing a sterol regulatory element. Proc Natl Acad Sci. 2000; 97: 10619–10624. https://doi.org/10.1073/pnas.180306597 PMID: 10962028

76. Printz R, Magnuson M, Granner D. Mammalian glucokinase. Annu Rev Nutr. 1993; 13: 463–496.
https://doi.org/10.1146/annurev.nu.13.070193.002335 PMID: 8369155

77. Zhu X, Chen D, Hu Y, Wu P, Wang K. The microRNA signature in response to nutrient restriction and refeeding in skeletal muscle of Chinese Perch (Siniperca chuatsi). Mar Biotechnol. 2015; 17: 180–189.
https://doi.org/10.1007/s10126-014-9606-8 PMID: 25403496

78. Li T, Franc JM, Boehme S, Chiang JYL. Regulation of cholesterol and bile acid homeostasis by the cho- lesterol 7α-hydroxylase/steroid response element-binding protein 2/microRNA-33a axis in mice. Hepa- tology. 2013; 58: 1111–1121. https://doi.org/10.1002/hep.26427 PMID: 23536474

79. Bommer GT, MacDougald OA. Regulation of lipid homeostasis by the bifunctional SREBF2-miR33a locus. Cell Metab. 2011; 13: 241–247. https://doi.org/10.1016/j.cmet.2011.02.004 PMID: 21356514

80. Krützfeldt J, Rajewsky N, Braich R, Rajeev KG, Tuschl T, Manoharan M, et al. Silencing of microRNAs in vivo with ‘antagomirs.’ Nature. 2005; 438: 685–689. https://doi.org/10.1038/nature04303 PMID: 16258535

81. Mittelbrunn M, Gutiérrez-Vázquez C, Villarroya-Beltri C, González S, Sánchez-Cabo F, González MÁ, et al. Unidirectional transfer of microRNA-loaded exosomes from T cells to antigen-presenting cells. Nat Commun. 2011; 2: 282. https://doi.org/10.1038/ncomms1285 PMID: 21505438
82. Huang Z, Huang D, Ni S, Peng Z, Sheng W, Du X. Plasma microRNAs are promising novel biomarkers for early detection of colorectal cancer. Int J cancer. 2010; 127: 119-126. https://doi.org/10.1002/ijc.25007 PMID: 19876917

83. Dietschy JM, Turley SD, Spady DK. Role of liver in the maintenance of cholesterol and low density lipoprotein homeostasis in different animal species, including humans. J Lipid Res. 1993; 34: 1637-1659. Available: http://www.jlr.org/content/34/10/1637.long PMID: 8245716