Photoreception and transcriptomic response to light during early development of a teleost with a life cycle tightly controlled by seasonal changes in photoperiod.

Light cues vary along the axis of periodicity, intensity and spectrum and perception of light is dependent on the photoreceptive capacity encoded within the genome and the opsins expressed. A global approach was taken to analyze the photoreceptive capacity and the effect of differing light conditions on a developing teleost prior to first feeding. The transcriptomes of embryos and alevins of Atlantic salmon (Salmo salar) exposed to different light conditions were analyzed, including a developmental series and a circadian profile. The results showed that genes mediating nonvisual photoreception are present prior to hatching when the retina is poorly differentiated. The clock genes were expressed early, but the circadian profile showed that only two clock genes were significantly cycling before first feeding. Few genes were differentially expressed between day and night within a light condition; however, many genes were significantly different between light conditions, indicating that light environment has an impact on the transcriptome during early development. Comparing the transcriptome data from constant conditions to periodicity of white light or different colors revealed overrepresentation of genes related to photoreception, eye development, muscle contraction, degradation of metabolites and cell cycle among others, and in constant light, several clock genes were upregulated. In constant white light and periodicity of green light, genes associated with DNA replication, chromatin remodeling, cell division and DNA repair were downregulated. The study implies a direct influence of light conditions on the transcriptome profile at early developmental stages, by a complex photoreceptive system where few clock genes are cycling.

Quantitative proteome profiling reveals molecular hallmarks of egg quality in Atlantic halibut: impairments of transcription and protein folding impede protein and energy homeostasis during early development.

BACKGROUND: Tandem mass tag spectrometry (TMT labeling-LC-MS/MS) was utilized to examine the global proteomes of Atlantic halibut eggs at the 1-cell-stage post fertilization. Comparisons were made between eggs judged to be of good quality (GQ) versus poor quality (BQ) as evidenced by their subsequent rates of survival for 12 days. Altered abundance of selected proteins in BQ eggs was confirmed by parallel reaction monitoring spectrometry (PRM-LC-MS/MS). Correspondence of protein levels to expression of related gene transcripts was examined via qPCR. Potential mitochondrial differences between GQ and BQ eggs were assessed by transmission electron microscopy (TEM) and measurements of mitochondrial DNA (mtDNA) levels. RESULTS: A total of 115 proteins were found to be differentially abundant between GQ and BQ eggs. Frequency distributions of these proteins indicated higher protein folding activity in GQ eggs compared to higher transcription and protein degradation activities in BQ eggs. BQ eggs were also significantly enriched with proteins related to mitochondrial structure and biogenesis. Quantitative differences in abundance of several proteins with parallel differences in their transcript levels were confirmed in egg samples obtained over three consecutive reproductive seasons. The observed disparities in global proteome profiles suggest impairment of protein and energy homeostasis related to unfolded protein response and mitochondrial stress in BQ eggs. TEM revealed BQ eggs to contain significantly higher numbers of mitochondria, but differences in corresponding genomic mtDNA (mt-nd5 and mt-atp6) levels were not significant. Mitochondria from BQ eggs were significantly smaller with a more irregular shape and a higher number of cristae than those from GQ eggs. CONCLUSION: The results of this study indicate that BQ Atlantic halibut eggs are impaired at both transcription and translation levels leading to endoplasmic reticulum and mitochondrial disorders. Observation of these irregularities over three consecutive reproductive seasons in BQ eggs from females of diverse background, age and reproductive experience indicates that they are a hallmark of poor egg quality. Additional research is needed to discover when in oogenesis and under what circumstances these defects may arise. The prevalence of this suite of markers in BQ eggs of diverse vertebrate species also begs investigation.

Heterochiasmy and the establishment of gsdf as a novel sex determining gene in Atlantic halibut.

Atlantic Halibut (Hippoglossus hippoglossus) has a X/Y genetic sex determination system, but the sex determining factor is not known. We produced a high-quality genome assembly from a male and identified parts of chromosome 13 as the Y chromosome due to sequence divergence between sexes and segregation of sex genotypes in pedigrees. Linkage analysis revealed that all chromosomes exhibit heterochiasmy, i.e. male-only and female-only meiotic recombination regions (MRR/FRR). We show that FRR/MRR intervals differ in nucleotide diversity and repeat class content and that this is true also for other Pleuronectidae species. We further show that remnants of a Gypsy-like transposable element insertion on chr13 promotes early male specific expression of gonadal somatic cell derived factor (gsdf). Less than 4.5 MYA, this male-determining element evolved on an autosomal FRR segment featuring pre-existing male meiotic recombination barriers, thereby creating a Y chromosome. Our findings indicate that heterochiasmy may facilitate the evolution of genetic sex determination systems relying on linkage of sexually antagonistic loci to a sex-determining factor.

Cardiac dysfunction affects eye development and vision by reducing supply of lipids in fish.

Developing organisms are especially vulnerable to environmental stressors. Crude oil exposure in early life stages of fish result in multiple functional and developmental defects, including cardiac dysfunction and abnormal and smaller eyes. Phenanthrene (Phe) has a reversible impact on cardiac function, and under exposure Phe reduces cardiac contractility. Exposure to a known L-type channel blocker, nicardipine hydrochloride (Nic) also disrupts cardiac function and creates eye deformities. We aimed to investigate whether cardiac dysfunction was the major underlying mechanism of crude oil-, Phe- and Nic-induced eye malformations. We exposed Atlantic haddock (Melanogrammus aeglefinus) early embryos to Nic and crude oil (Oil) and late embryos/early larvae to Phe exposure. All three exposures resulted in cardiac abnormalities and lead to severe, eye, jaw and spinal deformities at early larval stages. At 3 days post hatching, larvae from the exposures and corresponding controls were dissected. Eyes, trunk, head and yolk sac were subjected to lipid profiling, and eyes were also subjected to transcriptomic profiling. Among most enriched pathways in the eye transcriptomes were fatty acid metabolism, calcium signaling and phototransduction. Changes in lipid profiles and the transcriptome suggested that the dysfunctional and abnormal eyes in our exposures were due to both disruption of signaling pathways and insufficient supply of essential fatty acids and other nutrients form the yolk.

The salmon louse genome: Copepod features and parasitic adaptations.

Copepods encompass numerous ecological roles including parasites, detrivores and phytoplankton grazers. Nonetheless, copepod genome assemblies remain scarce. Lepeophtheirus salmonis is an economically and ecologically important ectoparasitic copepod found on salmonid fish. We present the 695.4 Mbp L. salmonis genome assembly containing ≈60% repetitive regions and 13,081 annotated protein-coding genes. The genome comprises 14 autosomes and a ZZ-ZW sex chromosome system. Assembly assessment identified 92.4% of the expected arthropod genes. Transcriptomics supported annotation and indicated a marked shift in gene expression after host attachment, including apparent downregulation of genes related to circadian rhythm coinciding with abandoning diurnal migration. The genome shows evolutionary signatures including loss of genes needed for peroxisome biogenesis, presence of numerous FNII domains, and an incomplete heme homeostasis pathway suggesting heme proteins to be obtained from the host. Despite repeated development of resistance against chemical treatments L. salmonis exhibits low numbers of many genes involved in detoxification.

Comparison of anadromous and landlocked Atlantic salmon genomes reveals signatures of parallel and relaxed selection across the Northern Hemisphere.

Most Atlantic salmon (Salmo salar L.) populations follow an anadromous life cycle, spending early life in freshwater, migrating to the sea for feeding, and returning to rivers to spawn. At the end of the last ice age ~10,000 years ago, several populations of Atlantic salmon became landlocked. Comparing their genomes to their anadromous counterparts can help identify genetic variation related to either freshwater residency or anadromy. The objective of this study was to identify consistently divergent loci between anadromous and landlocked Atlantic salmon strains throughout their geographical distribution, with the long-term aim of identifying traits relevant for salmon aquaculture, including fresh and seawater growth, omega-3 metabolism, smoltification, and disease resistance. We used a Pool-seq approach (n = 10-40 individuals per population) to sequence the genomes of twelve anadromous and six landlocked Atlantic salmon populations covering a large part of the Northern Hemisphere and conducted a genomewide association study to identify genomic regions having been under different selection pressure in landlocked and anadromous strains. A total of 28 genomic regions were identified and included cadm1 on Chr 13 and ppargc1a on Chr 18. Seven of the regions additionally displayed consistently reduced heterozygosity in fish obtained from landlocked populations, including the genes gpr132, cdca4, and sertad2 on Chr 15. We also found 16 regions, including igf1 on Chr 17, which consistently display reduced heterozygosity in the anadromous populations compared to the freshwater populations, indicating relaxed selection on traits associated with anadromy in landlocked salmon. In conclusion, we have identified 37 regions which may harbor genetic variation relevant for improving fish welfare and quality in the salmon farming industry and for understanding life-history traits in fish.

Illuminating the planktonic stages of salmon lice: A unique fluorescence signal for rapid identification of a rare copepod in zooplankton assemblages.

Monitoring of planktonic salmon louse (Lepeophtheirus salmonis salmonis) abundance and parameterization of key life-history traits has been hindered by labour-intensive and error-prone quantification using traditional light microscopy. Fluorescence illumination has been proposed as a means of improving visualization, but prior to this study adequate investigation of the relevant fluorescence profiles and measurement conditions has not been undertaken. We investigated the fluorescence profiles of L. salmonis and non-target copepod spp. with excitation and emission matrices (200-600 nm) and identified unique fluorescence signals. Fluorescence microscopy using excitation wavelengths of 470 ± 40 nm, and emission wavelengths of 525 ± 50 nm, showed that after 90 days of formalin storage salmon lice have a mean fluorescence intensity that is 2.4 times greater than non-target copepods (copepodid and adult stages). A 7-day heat treatment of 42°C in formalin increased the difference between salmon louse copepodids and non-target copepods to a factor of 3.6, eliminating the need for prolonged storage. Differences in the fluorescence signal and endogenous fluorophores were investigated with respect to variation in sea lice species, age, stage and host fish origin. Under the conditions outlined in this paper, the fluorescence signal was found to be a reliable means of visualizing and differentiating salmon lice from non-target zooplankters. Adaptation of the fluorescence signal would greatly expedite traditional methods of enumerating salmon louse larvae in plankton samples and could provide a means of automated detection.

Insulin-like 3 affects zebrafish spermatogenic cells directly and via Sertoli cells.

Pituitary hormones can use local signaling molecules to regulate target tissue functions. In adult zebrafish testes, follicle-stimulating hormone (Fsh) strongly increases the production of insulin-like 3 (Insl3), a Leydig cell-derived growth factor found in all vertebrates. Little information is available regarding Insl3 function in adult spermatogenesis. The Insl3 receptors Rxfp2a and 2b were expressed by type A spermatogonia and Sertoli and myoid cells, respectively, in zebrafish testis tissue. Loss of insl3 increased germ cell apoptosis in males starting at 9 months of age, but spermatogenesis appeared normal in fully fertile, younger adults. Insl3 changed the expression of 409 testicular genes. Among others, retinoic acid (RA) signaling was up- and peroxisome proliferator-activated receptor gamma (Pparg) signaling was down-regulated. Follow-up studies showed that RA and Pparg signaling mediated Insl3 effects, resulting in the increased production of differentiating spermatogonia. This suggests that Insl3 recruits two locally active nuclear receptor pathways to implement pituitary (Fsh) stimulation of spermatogenesis.

ReCodLiver0.9: Overcoming Challenges in Genome-Scale Metabolic Reconstruction of a Non-model Species.

The availability of genome sequences, annotations, and knowledge of the biochemistry underlying metabolic transformations has led to the generation of metabolic network reconstructions for a wide range of organisms in bacteria, archaea, and eukaryotes. When modeled using mathematical representations, a reconstruction can simulate underlying genotype-phenotype relationships. Accordingly, genome-scale metabolic models (GEMs) can be used to predict the response of organisms to genetic and environmental variations. A bottom-up reconstruction procedure typically starts by generating a draft model from existing annotation data on a target organism. For model species, this part of the process can be straightforward, due to the abundant organism-specific biochemical data. However, the process becomes complicated for non-model less-annotated species. In this paper, we present a draft liver reconstruction, ReCodLiver0.9, of Atlantic cod (Gadus morhua), a non-model teleost fish, as a practicable guide for cases with comparably few resources. Although the reconstruction is considered a draft version, we show that it already has utility in elucidating metabolic response mechanisms to environmental toxicants by mapping gene expression data of exposure experiments to the resulting model.

Transcriptomic analysis of dead end knockout testis reveals germ cell and gonadal somatic factors in Atlantic salmon.

BACKGROUND: Sustainability challenges are currently hampering an increase in salmon production. Using sterile salmon can solve problems with precocious puberty and genetic introgression from farmed escapees to wild populations. Recently sterile salmon was produced by knocking out the germ cell-specific dead end (dnd). Several approaches may be applied to inhibit Dnd function, including gene knockout, knockdown or immunization. Since it is challenging to develop a successful treatment against a gene product already existing in the body, alternative targets are being explored. Germ cells are surrounded by, and dependent on, gonadal somatic cells. Targeting genes essential for the survival of gonadal somatic cells may be good alternative targets for sterility treatments. Our aim was to identify and characterize novel germ cell and gonadal somatic factors in Atlantic salmon. RESULTS: We have for the first time analysed RNA-sequencing data from germ cell-free (GCF)/dnd knockout and wild type (WT) salmon testis and searched for genes preferentially expressed in either germ cells or gonadal somatic cells. To exclude genes with extra-gonadal expression, our dataset was merged with available multi-tissue transcriptome data. We identified 389 gonad specific genes, of which 194 were preferentially expressed within germ cells, and 11 were confined to gonadal somatic cells. Interestingly, 5 of the 11 gonadal somatic transcripts represented genes encoding secreted TGF-ß factors; gsdf, inha, nodal and two bmp6-like genes, all representative vaccine targets. Of these, gsdf and inha had the highest transcript levels. Expression of gsdf and inha was further confirmed to be gonad specific, and their spatial expression was restricted to granulosa and Sertoli cells of the ovary and testis, respectively. Finally, we show that inha expression increases with puberty in both ovary and testis tissue, while gsdf expression does not change or decreases during puberty in ovary and testis tissue, respectively. CONCLUSIONS: This study contributes with transcriptome data on salmon testis tissue with and without germ cells. We provide a list of novel and known germ cell- and gonad somatic specific transcripts, and show that the expression of two highly active gonadal somatic secreted TGF-ß factors, gsdf and inha, are located within granulosa and Sertoli cells.

Comparative transcriptome analysis of trout skin pigment cells.

BACKGROUND: Enormous variability in skin colour and patterning is a characteristic of teleost fish, including Salmonidae fishes, which present themselves as a suitable model for studying mechanisms of pigment patterning. In order to screen for candidate genes potentially involved in the specific skin pigment pattern in marble trout (labyrinthine skin pattern) and brown trout (spotted skin pattern), we conducted comparative transcriptome analysis between differently pigmented dermis sections of the adult skin of the two species. RESULTS: Differentially expressed genes (DEGs) possibly associated with skin pigment pattern were identified. The expression profile of 27 DEGs was further tested with quantitative real-time PCR on a larger number of samples. Expression of a subset of ten of these genes was analysed in hybrid (marble x brown) trout individuals and compared with the complexity of their skin pigment pattern. A correlation between the phenotype and the expression profile assessed for hybrid individuals was detected for four (gja5, clcn2, cdkn1a and tjp1) of the ten candidate genes tested. The potential role of these genes in skin pigment pattern maintenance is discussed. CONCLUSIONS: Our results indicate that the maintenance of different pigment patterns in trout is dependent upon specific communication-involving gap junctions, tight junctions and ion channels-between chromatophores present in differentially pigmented skin regions.

Transcriptome-wide mapping of signaling pathways and early immune responses in lumpfish leukocytes upon in vitro bacterial exposure.

We performed RNA sequencing, identified components of the immune system and mapped early immune responses of lumpfish (Cyclopterus lumpus) leukocytes following in vitro exposure to the pathogenic bacterium Vibrio anguillarum O1. This is the first characterization of immune molecules in lumpfish at the gene level. In silico analyses revealed that genes encoding proteins involved in pathogen recognition, cell signaling and cytokines in mammals and teleosts are conserved in lumpfish. Unique molecules were also identified. Pathogen recognition components include 13 TLRs, several NLRs and complement factors. Transcriptome-wide analyses of immune responses 6 and 24 hours post bacterial exposure revealed differential expression of 9033 and 15225 genes, respectively. These included TLR5S, IL-1ß, IL-8, IL-6, TNFα, IL-17A/F3, IL-17C and several components of the complement system. The data generated will be valuable for comparative studies and make an important basis for further functional analyses of immune and pathogenicity mechanisms. Such knowledge is also important for design of immunoprophylactic measures in lumpfish, a species of fish now farmed intensively for use as cleaner-fish in Atlantic salmon (Salmo salar) aquaculture.

Vgll3 and the Hippo pathway are regulated in Sertoli cells upon entry and during puberty in Atlantic salmon testis.

Vgll3 is linked to age at maturity in Atlantic salmon (Salmo salar). However, the molecular mechanisms involving Vgll3 in controlling timing of puberty as well as relevant tissue and cell types are currently unknown. Vgll3 and the associated Hippo pathway has been linked to reduced proliferation activity in different tissues. Analysis of gene expression reveals for the first time that vgll3 and several members of the Hippo pathway were down-regulated in salmon testis during onset of puberty and remained repressed in maturing testis. In the gonads, we found expression in Sertoli and granulosa cells in males and females, respectively. We hypothesize that vgll3 negatively regulates Sertoli cell proliferation in testis and therefore acts as an inhibitor of pubertal testis growth. Gonadal expression of vgll3 is located to somatic cells that are in direct contact with germ cells in both sexes, however our results indicate sex-biased regulation of vgll3 during puberty.

Novel adverse outcome pathways revealed by chemical genetics in a developing marine fish.

Crude oil spills are a worldwide ocean conservation threat. Fish are particularly vulnerable to the oiling of spawning habitats, and crude oil causes severe abnormalities in embryos and larvae. However, the underlying mechanisms for these developmental defects are not well understood. Here, we explore the transcriptional basis for four discrete crude oil injury phenotypes in the early life stages of the commercially important Atlantic haddock (Melanogrammus aeglefinus). These include defects in (1) cardiac form and function, (2) craniofacial development, (3) ionoregulation and fluid balance, and (4) cholesterol synthesis and homeostasis. Our findings suggest a key role for intracellular calcium cycling and excitation-transcription coupling in the dysregulation of heart and jaw morphogenesis. Moreover, the disruption of ionoregulatory pathways sheds new light on buoyancy control in marine fish embryos. Overall, our chemical-genetic approach identifies initiating events for distinct adverse outcome pathways and novel roles for individual genes in fundamental developmental processes.

bmp15l, figla, smc1bl, and larp6l are preferentially expressed in germ cells in Atlantic salmon (Salmo salar L.).

Atlantic salmon is a valuable commercial aquaculture species that would benefit economically and environmentally by controlling precocious puberty and preventing escapees from reproducing with wild populations. One solution to both these challenges is the production of sterile individuals by inhibiting the formation of germ cells, but achieving this requires more information on the specific factors that control germ cell formation. Here, we identified and characterized novel factors that are preferentially expressed in Atlantic salmon germ cells by screening for gonad-specific genes using available adult multi-tissue transcriptomes. We excluded genes with expression in tissues other than gonads based on quantity of reads, and then a subset of genes was selected for verification in a multi-tissue PCR screen. Four gonad-specific genes (bmp15l, figla, smc1bl, and larp6l) were chosen for further characterization, namely: germ cell specificity, investigated by comparing mRNA abundance in wild-type and germ cell-free gonads by quantitative real-time PCR, and cellular location, visualized by in situ hybridization. All four genes were expressed in both testis and ovary, and preferentially within the germ cells of both sexes. These genes may be essential players in salmon germ cell development, and could be important for future studies aiming to understand and control reproduction. Mol. Reprod. Dev. 84: 76-87, 2017. © 2016 Wiley Periodicals, Inc.

Atlantic salmon populations reveal adaptive divergence of immune related genes - a duplicated genome under selection.

BACKGROUND: Populations of Atlantic salmon display highly significant genetic differences with unresolved molecular basis. These differences may result from separate postglacial colonization patterns, diversifying natural selection and adaptation, or a combination. Adaptation could be influenced or even facilitated by the recent whole genome duplication in the salmonid lineage which resulted in a partly tetraploid species with duplicated genes and regions. RESULTS: In order to elucidate the genes and genomic regions underlying the genetic differences, we conducted a genome wide association study using whole genome resequencing data from eight populations from Northern and Southern Norway. From a total of ~4.5 million sequencing-derived SNPs, more than 10 % showed significant differentiation between populations from these two regions and ten selective sweeps on chromosomes 5, 10, 11, 13-15, 21, 24 and 25 were identified. These comprised 59 genes, of which 15 had one or more differentiated missense mutation. Our analysis showed that most sweeps have paralogous regions in the partially tetraploid genome, each lacking the high number of significant SNPs found in the sweeps. The most significant sweep was found on Chr 25 and carried several missense mutations in the antiviral mx genes, suggesting that these populations have experienced differing viral pressures. Interestingly the second most significant sweep, found on Chr 5, contains two genes involved in the NF-KB pathway (nkap and nkrf), which is also a known pathogen target that controls a large number of processes in animals. CONCLUSION: Our results show that natural selection acting on immune related genes has contributed to genetic divergence between salmon populations in Norway. The differences between populations may have been facilitated by the plasticity of the salmon genome. The observed signatures of selection in duplicated genomic regions suggest that the recently duplicated genome has provided raw material for evolutionary adaptation.

Expression profiling identifies Sertoli and Leydig cell genes as Fsh targets in adult zebrafish testis.

Spermatogonial stem cells are quiescent, undergo self-renewal or differentiating divisions, thereby forming the cellular basis of spermatogenesis. This cellular development is orchestrated by follicle-stimulating hormone (FSH), through the production of Sertoli cell-derived factors, and by Leydig cell-released androgens. Here, we investigate the transcriptional events induced by Fsh in a steroid-independent manner on the restart of zebrafish (Danio rerio) spermatogenesis ex vivo, using testis from adult males where type A spermatogonia were enriched by estrogen treatment in vivo. Under these conditions, RNA sequencing preferentially detected differentially expressed genes in somatic/Sertoli cells. Fsh-stimulated spermatogonial proliferation was accompanied by modulating several signaling systems (i.e. Tgf-ß, Hedgehog, Wnt and Notch pathways). In silico protein-protein interaction analysis indicated a role for Hedgehog family members potentially integrating signals from different pathways during fish spermatogenesis. Moreover, Fsh had a marked impact on metabolic genes, such as lactate and fatty acid metabolism, or on Sertoli cell barrier components. Fish Leydig cells express the Fsh receptor and one of the most robust Fsh-responsive genes was insulin-like 3 (insl3), a Leydig cell-derived growth factor. Follow-up work showed that recombinant zebrafish Insl3 mediated pro-differentiation effects of Fsh on spermatogonia in an androgen-independent manner. Our experimental approach allowed focusing on testicular somatic genes in zebrafish and showed that the activity of signaling systems known to be relevant in stem cells was modulated by Fsh, providing promising leads for future work, as exemplified by the studies on Insl3.

The two-step development of a duplex retina involves distinct events of cone and rod neurogenesis and differentiation.

Unlike in mammals, persistent postembryonic retinal growth is a characteristic feature of fish, which includes major remodeling events that affect all cell types including photoreceptors. Consequently, visual capabilities change during development, where retinal sensitivity to different wavelengths of light (photopic vision), -and to limited photons (scotopic vision) are central capabilities for survival. Differently from well-established model fish, Atlantic cod has a prolonged larval stage where only cone photoreceptors are present. Rods do not appear until juvenile transition (metamorphosis), a hallmark of indirect developing species. Previously we showed that whole gene families of lws (red-sensitive) and sws1 (UV-sensitive) opsins have been lost in cod, while rh2a (green-sensitive) and sws2 (blue-sensitive) genes have tandem duplicated. Here, we provide a comprehensive characterization of a two-step developing duplex retina in Atlantic cod. The study focuses on cone subtype dynamics and delayed rod neurogenesis and differentiation in all cod life stages. Using transcriptomic and histological approaches we show that different opsins disappear in a topographic manner during development where central to peripheral retina is a key axis of expressional change. Early cone differentiation was initiated in dorso-temporal retina different from previously described in fish. Rods first appeared during initiation of metamorphosis and expression of the nuclear receptor transcription factor nr2e3-1, suggest involvement in rod specification. The indirect developmental strategy thus allows for separate studies of cones and rods development, which in nature correlates with visual changes linked to habitat shifts. The clustering of key retinal genes according to life stage, suggests that Atlantic cod with its sequenced genome may be an important resource for identification of underlying factors required for development and function of photopic and scotopic vision.

First feed affects the expressions of microRNA and their targets in Atlantic cod.

To our knowledge, there is no report on microRNA (miRNA) expression and their target analysis in relation to the type of the first feed and its effect on the further growth of fish. Atlantic cod (Gadus morhua) larvae have better growth and development performance when fed natural zooplankton as a start-feed, as compared with those fed typical aquaculture start-feeds. In our experiment, two groups of Atlantic cod larvae were fed reference feed (zooplankton, mostly copepods, filtered from a seawater pond) v. aquaculture feeds: enriched rotifers (Brachionus sp.) and later brine shrimp (Artemia salina). We examined the miRNA expressions of six defined developmental stages as determined and standardised by body length from first feeding for both diet groups. We found eight miRNA (miR-9, miR-19a, miR-130b, miR-146, miR-181a, miR-192, miR-206 and miR-11240) differentially expressed between the two feeding groups in at least one developmental stage. We verified the next-generation sequencing data using real-time RT-PCR. We found 397 putative targets (mRNA) to the differentially expressed miRNA; eighteen of these mRNA showed differential expression in at least one stage. The patterns of differentially expressed miRNA and their putative target mRNA were mostly inverse, but sometimes also concurrent. The predicted miRNA targets were involved in different pathways, including metabolic, phototransduction and signalling pathways. The results of this study provide new nutrigenomic information on the potential role of miRNA in mediating nutritional effects on growth during the start-feeding period in fish larvae.