Embryonic exposure to produced water can cause cardiac toxicity and deformations in Atlantic cod (Gadus morhua) and haddock (Melanogrammus aeglefinus) larvae.

Regular discharges of produced water from the oil and gas industry represents the largest direct discharge of effluent into the marine environment worldwide. Organic compound classes typically reported in produced water include saturated hydrocarbons, monoaromatic and polyaromatic hydrocarbons (MAHs, PAHs) as well as oxygenated compounds, such as phenols, acids and ketones. This forms a cocktail of known and suspect toxicants, but limited knowledge is yet available on the sub-lethal toxicity of produced water to cold-water marine fish species. In the present work, we conducted a 4-day exposure of embryos of Atlantic cod (Gadus morhua) and haddock (Melanogrammus aeglefinus) to produced water extracts equivalent to 1:50, 1:500 and 1:5000 times dilutions of raw effluent. No significant reduction in survival or hatching success was observed, however, for cod, hatching was initiated earlier for exposed embryos in a concentration-dependent manner. During recovery, significantly reduced embryonic heart rate was observed for both species. After hatch, larvae subjected to embryonic exposure to produced water extracts were smaller, and displayed signs of cardiotoxicity, jaw and craniofacial deformations. In order to improve risk assessment and regulation of produced water discharges, it is important to identify which produced water components contribute to these effects.

Automatic determination of heart rates from microscopy videos of early life stages of fish.

Toxic effects of organic hydrophobic contaminants include impacts on fish heart rate (HR) and cardiac functioning. Thus, in ecotoxicology as well as aquaculture and even medicine, fish heart functioning plays an important role in application areas. The aim of this study was to assemble a pipeline of image processing and statistical techniques to extract HR information from microscopy videos of the embryo and larval stages of three species of fish (Atlantic cod, haddock, and Atlantic bluefin tuna). The method enables automatic processing for a large number of individuals, saving a significant amount of time compared with manual processing, while simultaneously eliminating the type of errors such a manual process might incur.

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.

Selection in the host structures the microbiota associated with developing cod larvae (Gadus morhua).

Marine fish larvae are immature upon hatching, and share their environment with high numbers of bacteria. The microbial communities associated with developing fish larvae might be structured by other factors than those important in developing terrestrial animals. Here, we analysed the beta (ß)-diversity of the microbiota associated with developing cod larvae and compared it with the bacterial communities in water and live feed by applying pyrosequencing of bar coded v4 16S rDNA amplicons. A total of 15 phyla were observed in the cod larval microbiota. Proteobacteria was the most abundant, followed by Firmicutes, Bacteroidetes and Actinobacteria. The composition and diversity of the cod larval microbiota changed considerably with age. The temporal and spatial patterns of ß-diversity could not be explained by stochastic processes, and did not coincide with changes in the rearing conditions. Furthermore, the larval microbiota was highly distinct from the water and the live feed microbiota, particularly at early developmental stages. However, the similarity between larval and water microbiota increased with age. This study suggests that strong selection in the host structures the cod larval microbiota. The changes in community structure observed with increasing age can be explained by altered selection pressure due to development of the intestinal system.

Variation in embryonic mortality and maternal transcript expression among Atlantic cod (Gadus morhua) broodstock: a functional genomics study.

Early life stage mortality is an important issue for Atlantic cod aquaculture, yet the impact of the cod maternal (egg) transcriptome on egg quality and mortality during embryonic development is poorly understood. In the present work, we studied embryonic mortality and maternal transcript expression using eggs from 15 females. Total mortality at 7days post-fertilization (7 dpf, segmentation stage) was used as an indice of egg quality. A 20,000 probe (20K) microarray experiment compared the 7hours post-fertilization (7 hpf, ~2-cell stage) egg transcriptome of the two lowest quality females (>90% mortality at 7 dpf) to that of the highest quality female (~16% mortality at 7 dpf). Forty-three microarray probes were consistently differentially expressed in both low versus high quality egg comparisons (25 higher expressed in low quality eggs, and 18 higher expressed in high quality eggs). The microarray experiment also identified many immune-relevant genes [e.g. interferon (IFN) pathway genes ifngr1 and ifrd1)] that were highly expressed in eggs of all 3 females regardless of quality. Twelve of the 43 candidate egg quality-associated genes, and ifngr1, ifrd1 and irf7, were included in a qPCR study with 7 hpf eggs from all 15 females. Then, the genes that were confirmed by qPCR to be greater than 2-fold differentially expressed between 7 hpf eggs from the lowest and highest quality females (dcbld1, ddc, and acy3 more highly expressed in the 2 lowest quality females; kpna7 and hacd1 more highly expressed in the highest quality female), and the 3 IFN pathway genes, were included in a second qPCR study with unfertilized eggs. While some maternal transcripts included in these qPCR studies were associated with extremes in egg quality, there was little correlation between egg quality and gene expression when all females were considered. Both dcbld1 and ddc showed greater than 100-fold differences in transcript expression between females and were potentially influenced by family. The Atlantic cod ddc (dopa decarboxylase) complete cDNA was characterized, and has a 1461bp open reading frame encoding a 486 amino acid protein that contains all eight residues of the conserved pyridoxal 5'-phosphate binding site including the catalytic lysine. This study provides valuable new information and resources related to the Atlantic cod egg transcriptome. Some of these microarray-identified, qPCR-confirmed, Atlantic cod egg transcripts (e.g. ddc, kpna7) play important roles during embryonic development of other vertebrate species, and may have similar functions in Atlantic cod.

Ontogeny of redox regulation in Atlantic cod (Gadus morhua) larvae.

The reduction potential of a cell is related to its fate. Proliferating cells are more reduced than those that are differentiating, whereas apoptotic cells are generally the most oxidized. Glutathione is considered the most important cellular redox buffer and the average reduction potential (Eh) of a cell or organism can be calculated from the concentrations of glutathione (GSH) and glutathione disulfide (GSSG). In this study, triplicate groups of cod larvae at various stages of development (3 to 63 days post-hatch; dph) were sampled for analyses of GSSG/2GSH concentrations, together with activities of antioxidant enzymes and expression of genes encoding proteins involved in redox metabolism. The concentration of total GSH (GSH+GSSG) increased from 610 ± 100 to 1260 ± 150 µmol/kg between 7 and 14 dph and was then constant until 49 dph, after which it decreased to 810 ± 100 µmol/kg by 63 dph. The 14- to 49-dph period, when total GSH concentrations were stable, coincides with the proposed period of metamorphosis in cod larvae. The concentration of GSSG comprised approximately 1% of the total GSH concentration and was stable throughout the sampling series. This resulted in a decreasing Eh from -239 ± 1 to -262 ± 7 mV between 7 and 14 dph, after which it remained constant until 63 dph. The changes in GSH and Eh were accompanied by changes in the expression of several genes involved in redox balance and signaling, as well as changes in activities of antioxidant enzymes, with the most dynamic responses occurring in the early phase of cod larval development. It is hypothesized that metamorphosis in cod larvae starts with the onset of mosaic hyperplasia in the skeletal muscle at approximately 20 dph (6.8mm standard length (SL)) and ends with differentiation of the stomach and disappearance of the larval finfold at 40 to 50 dph (10-15 mm SL). Thus, metamorphosis in cod larvae seems to coincide with high and stable total concentrations of GSH.

Global transcriptome analysis identifies regulated transcripts and pathways activated during oogenesis and early embryogenesis in Atlantic cod.

The molecular mechanisms underlying oogenesis and maternally controlled embryogenesis in fish are not fully understood, especially in marine species. Our aim was to study the egg and embryo transcriptome during oogenesis and early embryogenesis in Atlantic cod. Follicles from oogenesis stages (pre-, early-, and late-vitellogenic), ovulated eggs, and two embryonic stages (blastula, gastrula) were collected from broodstock fish and fertilized eggs. Gene expression profiles were measured in a 44 K oligo microarray consisting of 23,000 cod genes. Hundreds of differentially expressed genes (DEGs) were identified in the follicle stages investigated, implicating a continuous accumulation and degradation of polyadenylated transcripts throughout oogenesis. Very few DEGs were identified from ovulated egg to blastula, showing a more stable maternal RNA pool in early embryonic stages. The highest induction of expression was observed between blastula and gastrula, signifying the onset of zygotic transcription. During early vitellogenesis, several of the most upregulated genes are linked to nervous system signaling, suggesting increasing requirements for ovarian synaptic signaling to stimulate the rapid growth of oocytes. Highly upregulated genes during late vitellogenesis are linked to protein processing, fat metabolism, osmoregulation, and arrested meiosis. One of the genes with the highest upregulation in the ovulated egg is involved in oxidative phosphorylation, reflecting increased energy requirements during fertilization and the first rapid cell divisions of early embryogenesis. In conclusion, this study provides a large-scale presentation of the Atlantic cod's maternally controlled transcriptome in ovarian follicles through oogenesis, ovulated eggs, and early embryos.

Knockdown of the germ cell factor Dead end induces multiple transcriptional changes in Atlantic cod (Gadus morhua) hatchlings.

The RNA binding protein Dead end (DnD) is essential for maintaining viable germ cells in vertebrates and silencing of the gene has been demonstrated to cause sterility in several mammalian and fish species. Here we investigated transcriptome changes in hatched larvae of Atlantic cod induced by DnD knockdown using morpholino oligonucleotides (MO) injected in two-cell embryos. Whereas no fluorescently labeled germ cells were shown in embryos coinjected with dnd MO and nanos3 3'UTR coupled to green fluorescent protein, DnD knockdown had no visible effect on the number and location of Vasa protein positive cells in larvae. However, quantitative real-time RT-PCR (qPCR) revealed decreased vasa, nanos3 and tudor domain containing protein 7 mRNA expression and genome-wide oligonucleotide microarray analyses indicated profound suppression of genes involved in development and regulation of the reproductive system. DnD morphants showed lowered expression of genes encoding proteins involved in lipid, retinoid, cholesterol and steroid metabolism, including those with roles in sex hormone metabolism. Biotransformation of lipophilic compounds appeared suppressed too, as evidenced by down-regulation of several key genes from the phases 1 and 2 detoxification pathways. Effects of DnD silencing were highly pleiotropic and consisted of endocrine and metabolic changes, massive induction of histones and suppression of diverse developmental processes, including erythropoiesis and formation of extracellular matrix. While transient inhibition of dnd mRNA translation did not block development of primordial germ cells until hatch, results suggested that ablation of DnD might have major indirect consequences, including suppression of reproductive functions.

Redox regulation in Atlantic cod (Gadus morhua) embryos developing under normal and heat-stressed conditions.

With regard to predicted oceanic warming, we studied the effects of heat stress on the redox system during embryonic development of Atlantic cod (Gadus morhua), with emphasis on the glutathione balance, activities of key antioxidant enzymes, and their mRNA levels. The embryos were incubated at optimal temperature for development (6 °C) or slightly above the threshold temperature (10 °C). The regulation of all the redox-related parameters measured at optimum development was highly dynamic and complex, indicating the importance of both maternal and zygotic contributions to maintaining redox equilibrium. Development at 10 °C caused a significantly higher mortality at the blastula and early gastrula stages, indicating severe stress. Measures of the glutathione redox couple showed a significantly more reduced state in embryos at 10 °C compared to 6 °C at the post-gastrula stages. Mean normalized expression of nrf2, trxred, g6pd, gclc, nox1, CuZnsod, and mt in embryos kept at 10 °C revealed stage-specific significantly reduced mRNA levels. Activities of antioxidant enzymes changed both during ontogenesis and in response to temperature, but did not correlate with mRNA levels. As the embryos need a tightly regulated redox environment to coordinate between growth and differentiation, these findings suggest that the altered redox balance might participate in inducing phenotypic changes caused by elevated temperature.

Profiling of key apoptotic, stress, and immune-related transcripts during embryonic and postembryonic development of Atlantic cod (Gadus morhua L.).

We investigated the profiles of 25 genes involved in apoptosis (bcl-x2, casp3, casp8, ccar1, mcl1, and tpt1), immunity (bty, cathl, ifng, il1b, il6, il8, il10, lyzg, and tfa), oxidative stress (cat, gpx4, gsh-px, hsp70, hsp90a, and sod1), and stress axis (crh, pomc, grl1, and mlr) during Atlantic cod development and compared the mRNA transcript levels between samples from farmed (FB) and wild broodstock (WB) using real-time polymerase chain reaction. The suitability of nine endogenous housekeeping genes and an external standard (luciferase) as reference genes was also evaluated. The cycle threshold values of all housekeeping genes differed significantly throughout Atlantic cod development. Fertilization and hatching rates were significantly higher in WB group (95 ± 1.8% and 89 ± 2.8%, respectively) compared with FB (75 ± 3.4% and 66 ± 3.2%, respectively). Eleven target genes, namely, ccar1, casp3, bcl-x2, mcl-1, cat, gsh-px, hsp70, sod1, lyzg, il8, and grl were expressed in both groups at fertilization stage, indicating their maternal transfer. Among them, transcripts of gsh-px were more abundant in WB eggs, while the expression of hsp70 was significantly higher in FB eggs. In FB larvae, expression of cat, hsp70, hsp90a, pomc, mlr, grl1, bclx2, and il6 was significantly higher at hatching and the expression of cat, gpx4, casp3 and ccar1 was significantly higher at first feeding stages, than in WB group. These findings give an insight into the expressional changes in certain category of genes involved in the embryonic development of Atlantic cod, which may eventually determine the ultimate quality of the larvae.

Expression of vasa and nanos3 during primordial germ cell formation and migration in Atlantic cod (Gadus morhua L.).

Primordial germ cells (PGCs), progenitors of gametes, are specified very early in embryonic development and undergo an active migration to the site where the future gonads will form. While the developmental pattern of PGCs during embryogenesis has been documented in few model teleost fishes, there is currently no information available for any representative of Superorder Paracanthopterygii. This includes Atlantic cod (Gadus morhua), which is a historically important food fish in both fisheries and aquaculture industries. In the present study, we cloned and characterized vasa and nanos3 and used them as germ cell markers in Atlantic cod. Sequencing results showed prospective vasa and nanos3 mRNA contained the domains used to describe their respective protein family. Furthermore, phylogenetic analysis using the amino acid sequence placed Atlantic cod Vasa distinct from representatives of three other taxonomic Superorders. Atlantic cod Nanos3 was placed with other homologues from the Nanos3 subfamily. Expression of both genes was detected from the first cleavage division; both were specifically expressed in Atlantic cod PGCs from the 32-cell stage. While nanos3 expression ceased during early somitogenesis, vasa was strongly expressed throughout embryonic development. Using vasa as a marker, we described the Atlantic cod PGC migration pattern. We demonstrated that Atlantic cod PGCs migrate ventral to the trunk mesoderm. With the exception of Pacific herring (Clupea pallasii), PGCs in other described teleost fishes migrate lateral to the trunk. The results from this study are the first step toward understanding germ line formation in Atlantic cod.

Low impact of exposure to environmentally relevant doses of 226Ra in Atlantic cod (Gadus morhua) embryonic cells.

The aim of this study was to investigate whether (226)Ra, a radionuclide present in produced water from oil platforms in the North Sea and other offshore drilling areas, could affect vulnerable early life stages of Atlantic cod (Gadus morhua). Blastula-stage embryonic cells (EC) from fertilized eggs of Atlantic cod were isolated and exposed to environmental relevant concentrations of (226)Ra and transcription of selected genes quantified. The results showed a weak, but significant up-regulation of GPx3 and HSP70 transcripts after 48 h of exposure to 2.11 Bq/L. In EC exposed to three (226)Ra concentrations (2.11, 23 and 117 Bq/L) for 12 h, metallothionein, HSP90AA, thioredoxin and caspase 8 were significantly up-regulated in cells exposed to 117 Bq/L, whereas thioredoxin was also significantly up-regulated in EC exposed to 23 Bq/L. When EC were exposed to the same (226)Ra concentrations for 48 h, only heme oxygenase was significantly up-regulated in the 23 Bq/L exposure group. The results suggest that environmentally relevant activities of (226)Ra may induce oxidative stress and apoptosis in fish ECs. Exposure of Atlantic cod EC to Cd, selected as a model toxicant, supported the ability of EC around blastula stage to respond to toxicants by altered transcription. Due to dilution, environmentally relevant concentrations of radionuclides present in produced water would be expected to pose a minor threat to early life stages of fish.

Gene expression profiling of Atlantic cod (Gadus morhua) embryogenesis using microarray.

Atlantic cod (Gadus morhua) is a fish species of high importance, as a key species in a range of Northern ecosystems, in fisheries, and as an emerging species in aquaculture. So far, little is known about the transcriptional activity during early developmental stages of Atlantic cod. Hence, we decided to use a cDNA microarray covering 7,000 genes to analyze the temporal activity of the transcriptome during cod embryogenesis. Twelve different embryonic time points were selected, covering major developmental stages and processes such as maternally derived mRNA, blastula, gastrula, segmentation, hatching, and first-feeding larval stage. The microarray analysis revealed a highly dynamic transcriptional profile, showing for the first time the differential expression of thousands of known and unknown genes during Atlantic cod embryogenesis. These initial findings will serve as an important baseline for future in-depth studies of candidate genes involved in development, reproductive control, disease resistance, growth, nutrient digestion, and metabolism.

Ontogenetic expression of maternal and zygotic genes in Atlantic cod embryos under ambient and thermally stressed conditions.

The embryonic stages of Atlantic cod (Gadus morhua) are especially sensitive to incubation temperature. The purpose of the present study was to follow the ontogenetic expression of selected genes of maternal (pou2 and nanog) and zygotic origin (hsp70, hsp90α and stip1), in Atlantic cod embryos under ambient and thermally stressed conditions. The study also investigated how reference genes can be applied to studies on embryonic development, when maternal genes are degraded and the zygotic transcription stabilizes. Three batches of eggs were reared and gene expression profiles from the reference and target genes were determined. The embryos were reared at ambient 6 °C, and 10 °C for continuous long-term and acute short-term heat exposure. Both pou2 and nanog showed reduced expression whereas the zygotic and reference genes showed increased expression until stabilizing at gastrulation, when a normalized ontogenetic expression profile of target genes could be generated. pou2 and nanog were not affected by thermal stress. In contrast, hsp70 and hsp90α were upregulated after short-term heat exposure at the early blastula (hsp70 only), late blastula, 50% epiboly and 90% epiboly stages (hsp90α only). Long-term heat exposure of Atlantic cod embryos upregulated both hsp70 (90% epiboly) and hsp90α (90% epiboly and 20-somites). The results suggest that a cellular defense mechanism is activated even in the earliest stages of embryonic development, a period critical to developmental temperature.

Differential spatial expression of mef2 paralogs during cardiac development in Atlantic cod (Gadus morhua).

The myogenic enhancer factor 2 (Mef2) transcription factors are known for their role in the control of cardiac development. Here we describe the spatial and temporal expression patterns of five Atlantic cod mef2 genes designated as mef2a, mef2cI, mef2cII, mef2dI and mef2dII during cardiogenesis. Whole mount in situ hybridization showed that mef2a and mef2dI were expressed in both cardiac ring and cone prior to looping morphogenesis, while mef2dII expression was only detectable in the cardiac ring. The mef2cI and mef2cII paralogs displayed different spatial expression patterns in the heart tube with a venous and arterial pole preference, respectively. After the cardiac loop formation mef2cI was expressed in cells of the ventricle and lateral arteries, while mef2cII appeared more abundant and was also present in the atrium. Larvae raised at constant 8 °C showed malformed morphology of the lateral arteries, and the transcription of both mef2c variants was highly elevated compared to those kept at 4 °C. Acute temperature stress also resulted in deviations in the expression of the mef2c paralogs, and the treated embryos displayed defects in the developing heart, including impaired fusion of the bilateral primordia and truncated heart tubes.

Development of Atlantic cod (Gadus morhua) exposed to produced water during early life stages: Effects on embryos, larvae, and juvenile fish.

Produced water (PW) contains numerous toxic compounds of natural origin, such as dispersed oil, metals, alkylphenols (APs), and polycyclic aromatic hydrocarbons (PAHs). In addition, PW also contains many different chemicals which have been added during the oil production process. In the study described here, cod were exposed to real PW collected from an oil production platform in the North Sea. This was done in order to best recreate the most realistic field-exposure regime in which fish will be affected by a wide range of chemicals. The biological effects found in this study therefore cannot be assigned to one group of chemicals alone, but are the result of exposure to the complex chemical mixture found in real PW. Since APs are well known to cause endocrine disruption in marine organisms, we focused our chemical analysis on APs in an attempt to better understand the long-term effects of APs from PW on the biology of fish. In this study, cod were exposed to several concentrations of real PW and 17ß-oestradiol (E(2)), a natural oestrogen, at different developmental stages. Cod were exposed to PW either during the embryo and early larvae stage (up to 3 months of age) or during the early juvenile stage (from 3 to 6 months of age). Results showed that, in general, APs bioconcentrate in fish tissue in a dose and developmental stage dependent manner during PW exposure. However, juveniles appeared able to effectively metabolise the short chain APs. Importantly, PW exposure had no effect on embryo survival or hatching success. However, 1% PW clearly interfered with the development of normal larval pigmentation. After hatching most of the larvae exposed to 1% PW failed to begin feeding and died of starvation. This inability to feed may be linked to the increased incidence of jaw deformities seen in these larvae. In addition, cod exposed to 1% PW, had significantly higher levels of the biomarkers vitellogenin and CYP1A in plasma and liver, respectively. No similar effects were seen in cod exposed to either 0.1% or 0.01% PW.

Ribosomal RNA gene sequences confirm that protistan endoparasite of larval cod Gadus morhua is Ichthyodinium sp.

An enigmatic protistan endoparasite found in eggs and larvae of cod Gadus morhua and turbot Psetta maxima was isolated from Baltic cod larvae, and DNA was extracted for sequencing of the parasite's small subunit ribosomal RNA (SSU rRNA) gene. The endoparasite has previously been suggested to be related to Ichthyodinium chabelardi, a dinoflagellate-like protist that parasitizes yolk sacs of embryos and larvae of a variety of fish species. Comparison of a 1535 bp long fragment of the SSU rRNA gene of the cod endoparasite showed absolute identity with I. chabelardi, demonstrating that the 2 parasites are very closely related, if not identical. This finding is discussed in relation to some morphological differences that appear to exist between I. chabelardi and the cod endoparasite.

Embryonic stem cells isolated from Atlantic cod (Gadus morhua) and the developmental expression of a stage-specific transcription factor ac-Pou2.

The establishment of embryonic stem cell cultures and the identification of molecular markers for undifferentiated embryonic stem cells (ESC) as well as differentiated cells types will open new opportunities in the study of developmental biology and for developing embryonic in vitro models of the ecologically and economically important fish specie Atlantic cod (Gadus morhua). We report here that cod blastula cells express a Class V POU gene known to be highly expressed in embryonic cell populations of vertebrates. The cod transcript, designated Atlantic cod-Pou2 (ac-Pou2), can be used as a genetic marker for cod blastula cells in vivo and in vitro. Using a quantitative real-time PCR approach, we found that the ac-Pou2 transcript was downregulated before the egg reached the stage of gastrulation, the starting point of extensive cell differentiation. We also demonstrate the culturing of ESC isolated from cod blastula stage eggs. The cod ESC exhibited in vitro characteristics of pluripotency described for both mammalian ESC and fish ES-like cells (medaka, zebrafish, seabream, sea perch and rainbow trout). Cod ESC in culture expressed ac-pou2, differentiated spontaneously and had the ability to form embryoid bodies following retinoic acid treatment. The ESC could also be directed to differentiate.