Genome-wide microarray analysis of Atlantic cod (Gadus morhua) oocyte and embryo
BMC Genomics volume 15, Article number: 594 (2014)
Regulation of gene expression plays a central role in embryonic development. Early stages are controlled by gametic transcripts, which are subsequently substituted with transcripts from the genome of the zygote. Transcriptomic analyses provide an efficient approach to explore the temporal gene expression profiles in embryos and to search for the developmental regulators. We report a study of early Atlantic cod development that used a genome-wide oligonucleotide microarray to examine the composition and putative roles of polyadenylated transcripts.
The analyses were carried out in unfertilized oocytes, newly fertilized oocytes and embryos at the stages of mid-blastula transition and segmentation. Numerous genes transcribed in oocytes are involved in multiple aspects of cell maintenance and protection, including metabolism, signal perception and transduction, RNA processing, cell cycle, defense against pathogens and DNA damage. Transcripts found in unfertilized oocytes also encoded a large number of proteins implicated in cell adherence, tight junction and focal adhesion, suggesting high complexity in terms of structure and cellular interactions in embryos prior to midblastula transition (MBT). Prezygotic transcripts included multiple regulators that are most likely involved in developmental processes that take place long after fertilization, such as components of ErbB, hedgehog, notch, retinoid, TGFb, VEGF and Wnt signaling pathways, as well as transcripts involved in the development of nervous system. The major event of MBT was the activation of a large group of histones and other genes that modify chromatin structure preceding massive gene expression changes. A hallmark of events observed during segmentation was the induction of multiple transcription factors, including a large group of homeobox proteins in pace with decay of a large fraction of maternal transcripts. Microarray analyses detected a suite of master developmental regulators that control differentiation and maintenance of diverse cell lineages.
Transcriptome profiling of the early stages in Atlantic cod revealed the presence of transcripts involved in patterning and development of tissues and organs long before activation of the zygotic genome. The switch from maternal to zygotic developmental programs is associated with large-scale modification of chromosomes.
Early ontogeny is associated with dramatic gene expression changes that underlie and determine the developmental processes. Transcription terminates by the end of oogenesis when the maturing oocyte is arrested in the metaphase of its second meiotic division [1, 2]. The oocyte is loaded with maternal mRNAs and proteins that control the cell maintenance and fate and the formation of the body plan prior to the onset of zygotic genome expression [3, 4]. Important transcripts can be also contributed by sperm cell, as was recently shown in Drosophila and mammals [5, 6]. Today, it is generally thought that the combination of determinants deposited by the mother during oogenesis and the inductive signals between different cells trigger the specification of different cell lineages during development of the embryo [7, 8]. Maternal to zygotic transition (MZT) is the key event during embryogenesis marked by the switch of control from the maternal and possibly paternal transcripts to the newly synthesized embryonic gene products [9–11]. Degradation of maternal transcripts and zygotic genome activation is characterized by striking changes in the transcriptome profiles. MZT timing is species-specific according to the extent and form of maternal contributions and generally occurs earlier in mammals [12–15] compared to fish, Drosophila and Xenopus[16–19]. In a number of animal species, MZT roughly coincides with the mid-blastula transition (MBT)  when cells become motile and divide asynchronously. The three germ layers and the body plan of the mature organism are established during gastrulation, and the period is characterized by extensive cell movements and intracellular communications [21, 22]. During the following segmentation stage major events in the formation of tissues and organs take place.
Knowledge of the genetic networks controlling embryogenesis has been obtained principally by mutagenesis screens in model species. Multiple mutations affecting embryonic development have been induced by chemical and insertional mutagenesis resulting in the identification of genes with important roles in development in Drosophila[23–25]. Similarly, large-scale genetic screens in zebrafish have enhanced the overall understanding of critical steps and pathways during embryogenesis, and forward genetics revealed a number of developmentally regulated genes [26–28]. Despite high power, this research strategy encounters limitations because only indispensable genes whose loss cannot be compensated by functionally related genes are found, leaving many important actors undetected. A complementary approach is transcriptome profiling that reveals genes with characteristic temporal expression patterns. The completion of the Atlantic cod whole-genome sequencing project  enabled the development of novel tools for gene expression profiling of this ecologically and commercially important marine species sustaining wild fisheries and aquaculture. DNA microarrays are used for analyses of polyadenylated mRNA and a transcriptome study of Atlantic cod embryogenesis using a cDNA microarray was recently reported . We present herein the use of the Atlantic cod genome-wide oligonucleotide microarray for investigation of transcriptome changes associated with the key events of early development from unfertilized oocytes to late somitogenesis with focus on changes during MZT. Contribution of transcripts with different temporal profiles in diverse processes associated with maintenance and development was assessed and compared.
An overview of oocyte and embryo transcriptome
The microarray analyses of polyadenylated mRNA included four developmental stages: unfertilized oocytes (UFO), oocytes collected at 2 hours post fertilization (2hpf), the midblastula transition (MBT) and segmentation (SGM). The features that showed over 2-fold difference in comparison with reference (adult tissues) in at least one of the analysed stages were selected. These genes were categorized as prezygotic (high expression in UFO) and zygotic (activation after fertilization), and further divided in seven subgroups based on the decreased or increased abundance at the specific stages (Figure 1, Table 1, Additional file 1). Most prezygotic transcripts maintained relatively stable levels either within the whole period examined or until SGM. A small number of transcripts decreased or increased abundance shortly after fertilization. The origin of the transcripts detected at the 2 hpf is uncertain, but the zygote formation in Atlantic cod occurs approximately at 5 hpf . Therefore, most if not all these transcripts identified at 2 hpf most likely should be ascribed to the sperm cell. Genes activated at MBT and SGM comprised respectively 17.2% and 80.4% of all zygotic genes.
Metabolism, cell maintenance, proliferation and protection
Genes with metabolic roles comprised a large part of the developmentally regulated genes. Functional groups involved in degradation of proteins and RNA were represented almost exclusively by prezygotic transcripts that suggested a profound reconstruction of the cellular machinery (Figure 2A). This was in line with zygotic activation of multiple genes encoding components of the key cellular organelles and structures: cell surface, ribosomes and microsomes with the greatest changes observed in the endoplasmic reticulum. Proteins of cytoskeleton, mitochondria and lysosomes were predominantly dependent on maternal transcripts, as well as signalling pathways controlling metabolic processes. The greatest developmental regulation was seen in the pathways of lipid and cholesterol metabolism and PPAR signalling; a sharp induction at SGM was shown by a suite of apolipoproteins (Figure 2B). Proteins involved in cell cycle and apoptosis were encoded predominantly by the maternal transcripts. Cdc25B, which induces mitosis by promoting G2/M phase progression, was abundant during the early stages followed by gradual decreasing levels towards the SGM. Consistently, the elimination of cdc25B transcripts occurs by the end of the MZT in vertebrate and invertebrate species [32, 33]. Other maternal mitotic regulators included cyclins and retinoblastoma-associated protein, one of the key factors that controls the entry into cell cycle. Genes encoding the mitotic checkpoint serine/threonine kinase bub1 ensuring correct chromosomal segregation during the cell division and the wee1-like protein kinase regulating DNA replication prior to mitosis were expressed until MBT and are likely to be involved in securing the integrity of the genome prior to cell division . Several cell cycle related genes were found at 2hpf, including centromere protein J, which participates in centriole duplication , and the positive cell cycle regulator smad nuclear interacting protein. The negative regulator, cyclin-dependent kinase inhibitor 1C was induced at MBT in concordance with reduction of cell proliferation.
Early fish embryos possess a multifaceted defence system. A suite of immune genes was expressed at high levels already in UFO. This group included complement components, cytokines, chemokines and their receptors, IFN and TNF-related genes, together with three negative regulators of immunity from the SOCS family. Several immune genes including myeloperoxidase were activated during SGM. While protection from DNA damage was driven mainly by the maternal transcripts, responses to oxidative stress were switched on later, and seven of eight genes involved in regulation of redox homeostasis were induced during SGM (Figure 2C). Marked developmental regulation was shown by glutathione peroxidases and oxidation resistance protein coding genes.
Cell-to-cell and cell-to-extracellular matrix (ECM) interactions
The maternal transcripts were predominant among cell adhesion molecules and their input was high in other functional groups and pathways implicated in cell contacts and interactions (Figure 3A). Abundance of transcripts for cadherins (11 genes), claudins-12 and 14, and gap junction proteins was high in UFO and decreased sharply at SGM (Figure 3B). This stage was also marked with massive rearrangement of the extracellular space. Expression of two genes encoding the extracellular proteins, lens intrinsic membrane protein 2.3 and collagen XXV, was highly specific for oocytes, and the latter showed 227-fold greater abundance in comparison with adult tissues and remained highly abundant until SGM (Figure 3B). These proteins are scarcely explored and it is unknown whether they are required for the oocyte maturation or embryonic development. Late activation of decorin, laminin beta 1 and collagen type IV witnessed increased complexity and mechanical strength of the extracellular structures (Figure 3B). Up-regulation of collagen degrading matrix metalloproteinase 13 was probably required for remodeling of extracellular matrix. High induction was shown by four genes for keratin 12 that is commonly located in epithelium being involved in interactions between cells and ECM.
Transcriptional repression/activation and chromosomal remodeling
The UFO transcriptome encoded proteins involved in mRNA processing, some of which showed highly specific expression in oocytes. For example U2 small nuclear ribonucleoprotein B displayed 69-fold higher levels compared to adult tissues (Figure 4A). The maternal transcripts of RNA helicase DDX18 encode a protein taking part in spliceosome assembly suggesting that the regulation of mRNA splicing and processing continues even after arrest of transcription in oocytes and early embryos. Transcription repressors and activators seemed to be associated with male and female gametes. UFO included several histone-lysine N-methyltransferases H3, which have been implicated in transcriptional gene silencing, heterochromatin assembly and DNA methylation . UFO transcripts also encoded components of the polycomb repressive complex contributing to the formation of silent chromatin. Sex-comb on midleg-like proteins are required to maintain the transcriptionally repressive state of homeotic genes throughout the development . Further, the transcriptional repressors remodeling and spacing factor 1 (rsf1) and mbt domain containing 1 were present in unfertilized and fertilized oocytes, respectively. A maternally supplied co-activator bromodomain containing 2 (brd2) [40, 41] can be important for the early embryo cell cycle control. Fertilized oocytes contained transcripts for a protein involved in DNA methylation (PWWP domain containing 2 isoform 1). Though a number of genes related to chromosome maintenance and remodeling were present among prezygotic transcripts, major changes took place at MBT and a large fraction of the activated genes encoded histones (66 genes) (Figure 4B). Histone modifications and changes of chromatin architecture enable the formation of transcriptionally active euchromatin in order for zygotic gene expression to take place [20, 42]. Although some of the transcripts coding for histones H1, H2A, H2B, H3 and H4 were present at relatively high levels already in oocytes, most of them markedly increased abundance at MBT and many had biphasic profiles (Figure 4C). Several genes involved in nucleosome and chromatin remodeling, including chromodomain helicase DNA binding protein 8, histone deacetylase 11 and high mobility group protein b2 showed highest expression at SGM (Figure 4A).
Regulation of early cellular differentiation and signaling
Many developmentally regulated genes belong to pathways with crucial roles in embryonic development (Figure 5A). A large part of these genes are transcription factors with different temporal expression profiles (Figure 5B). Transcripts for jun-B that controls differentiation of diverse cell lineages [43, 44], pair box protein 7 and AP-2 involved in formation of eye, limb, neural and cardiac development  were highly abundant in oocytes. UFO also included transcripts for proteins that play important parts in development of heart (myogenic enhancer factor), cartilage and bone (runx2) and erythroid lineage (nuclear factor erythroid derived-2). Four transcription factors detected at 2hpf contained zinc finger domains. Greatest induction in SGM (90.3-fold) was shown by hes-5, a component of Notch pathway  and several homeodomain genes. Many identified transcription factors are unknown but include functional domains present in multiple regulators of development and interestingly, ratios between regulators with domains changed by stages (Figure 5C). While the majority of zinc finger proteins were found among prezygotic transcripts, 73% of transcription factors with homeodomain were activated during SGM, while forkhead proteins occupied an intermediate position.
In addition to transcription factors, the developmentally regulated genes included transcripts for receptors and extracellular proteins assigned to Notch, TGF beta and Wnt signaling pathway (Figure 6A). Major part of them were present in oocytes and only two genes were each activated at MBT and SGM. Greatest differences from adult tissues (51- and 62-fold) were observed in frizzled 8a, a receptor for Wnt proteins, and dorsal-ventral patterning tolloid-like protein. Dorsoventral patterning is also regulated by two genes from TGF pathway: follistatin and noggin. Maternal transcripts encoded three proteins from Smad family, which transmit signals from TGF. Notch are transmembrane proteins that bind jagged ligands controlling differentiation by receiving signals through cell to cell contacts, while deltex is a regulator of notch signaling. Su(H)B encodes the suppressor of hairless, a key transcriptional regulator of Notch pathway. Of note is oocyte expression of genes that regulated differentiation of complex structures, such as connective tissue growth factor (3 genes) and angio-associated migratory cell protein.
A suite of genes known for their roles in neurogenesis was detected already in UFO (Figure 6B). Greatest difference with adult tissues (59-, 34- and 33-fold) was shown by neural cell adhesion molecule L1-like protein (chl1) implicated in cell migration and neuronal positioning, neuropilin and contactin-1a, a neuronal cell adhesion molecule important for the formation of axon connections during the nervous system development [49, 50]. Maternally provided transcripts also included neurogenic differentiation factors (5 genes), which are involved in neuroepithelial stem cell differentiation and neurogenesis, the synaptic protein and receptor of neurotransmitters neuronal pentraxin-1 precursor (nptx1) and two ephrins and ephrin type-B receptors (4 genes), that play a crucial part in migration of axons. Interestingly, switch of ephrin and ephrin receptor isoforms took place at SGM. Maternal transcripts also encoded receptors of the dopamine and acetylcholine neurotransmitters.
Master regulators of embryogenesis
The microarray analyses revealed several genes that control stem cell fates and organ development, which were active at different stages (Figure 7). While transcript for pumilio-2 (proliferation and renewal of stem cells, germ cell development and degradation of maternal mRNAs ) was present in oocytes, snail homolog sna (formation and maintenance of mesoderm during embryogenesis ), RNA-binding protein musashi homolog 2 (proliferation and maintenance of stem cells in the CNS ) were activated at MBT and SGM, respectively. Similar to zebrafish , MBT was marked with an onset of nanor transcription. DiGeorge syndrome critical region gene 8 may be involved in the degradation of maternal transcripts via biogenesis of microRNAs . A role of microRNAs, such as miR430, in the clearance of maternal transcripts during the maternal to zygotic transition has been well documented [56, 57]. Sox-2 (3 genes) is a key regulator of embryonic stem cell pluripotency .
Validation of microarray data
To validate microarray data, real-time qPCR was performed on eight genes: nanor, heat shock protein 70 kDa protein 4 (hsp70), heat shock 90 kDa protein 1 beta isoform a (hsp90ba), stress-induced phosphoprotein 1 (stip1), follistatin (fst), formin-binding protein 4 (fnbp4), keratin-12 (krt12) and ikaros (ikzf) at the embryonic stages of UFO, 2hpf, MBT, 12- and 52-somites. Overall, microarrays and qPCR produced similar results (Figure 8).
Genome sequencing enabled construction of oligonucleotide microarrays that may provide complete coverage of the polyadenylated fraction of transcriptomes; microarray analyses evaluate abundance of mature mRNA, which is capable for translation. While a genome-wide platform was used for evaluating the abundance of mature mRNA during zebrafish development , we report the first study performed with an aquaculture fish species. The main issue was the presentation of pathways and functional groups among the transcripts displayinh different temporal profiles. High complexity of the transcriptome in unfertilized cod oocytes is consistent with similar studies in both invertebrates and vertebrates [10, 59, 60]. Maternally provided mRNA comprised the major part of prezygotic transcripts while the putative paternal contribution was small, but sperm transcripts might play important roles of in the establishing of early embryonic gene expression profiles [5, 61, 62]. Fertilized cod eggs contained a suite of transcripts for proteins involved in chromatin remodeling and regulation of transcription, cell cycle control and cellular transport. However, it is unknown whether these genes have any developmental roles. In general, maternal transcripts support basic requirements of the embryo prior to the onset of zygotic expression. Interestingly, we got evidence that processing of mRNA continues even in absence of transcription that is in line with recent report on large-scale maturation of maternal transcripts in zebrafish embryos [57, 63]. In addition to maintenance of metabolism, cell structure and proliferation, transcripts of oocytes provide immune protection against pathogens and a suite of genes is expressed at higher level in comparison with adult tissues. Maternal transfer of complement factors and their protective roles was reported in wolffish, rainbow trout and zebrafish [64–67]. The female fish also provide offspring with immunoglobulins, lysozymes, protease inhibitors and different types of lectins [68, 69]. The observed prevalence of immune genes involved in signaling suggests that embryos are capable to regulate responses to pathogens. Presentation of multiple signal transduction pathways points to active perception of external cues and complex interactions between early embryos and environment.
Differentiation presumes acquisition of specific properties by cells and increase of their heterogeneity. A large number of transcripts for proteins involved in cell contacts were abundant in UFO being eliminated at SGM. Cadherins are transmembrane cell adhesion proteins that mediate various processes during development including cellular migration and tissue organization . Interestingly, this study identified a large number of cadherin paralogs that are likely involved in cell sorting and tissue morphogenesis . UFO included many transcripts that can be involved in the control of processes taking place long after fertilization, such as components of Wnt, Notch, hedgehog, ErbB, TGF beta and VEGF signaling pathways and markers of specialized cell lines [72, 73]. We also identified multiple transcripts that may regulate neurogenesis or encode proteins known as highly specific for neural tissue, in agreement with a few studies reporting maternal deposition of transcripts later expressed in the CNS (e.g. Drosophila, zebrafish and axolotl [74–76]). This finding can partly be accounted for by the bias in annotation, since a number of genes with pleiotropic functions have been studied mainly in the context of nervous system. Furthermore, some genes could change functions in course of the vertebrate evolution as demonstrated by the identification several genes known as neural specific in mammals were primarily involved in innate antiviral responses in fish [77, 78].
Transcriptome analyses suggested that the onset of zygotic expression is preceded by large scale modification of chromosomes. Histones comprise a major fraction of genes activated during MBT while the number of transcription factors in this group was small. The pre-MBT transcripts encoded several proteins that modify histones and DNA and are known as positive and negative regulators (e.g. myst2, brd2, n6amt1, rsf1, ehmt3, scml-1, 2 and 4). Transcripts coding for histone methyltransferases and members of the polycomb repressors were highly abundant in unfertilized and fertilized oocytes, but showed a decrease in expression after MBT and coincided with the chromatin remodeling prior to the activation of transcription. Preparation of transcriptional machinery to the large-scale activation of gene expression appears a major developmental event that takes place during MBT. In most studied vertebrates this period coincides with the degradation of maternal transcripts and activation of the zygotic genome which takes over the genetic control of embryogenesis [10, 79]. Furthermore, accumulated studies reveal the dynamic nature of chromatin regulation and the importance of its modifications during transitions from maternal to zygotic control of development [20, 42, 80]. Our data are consistent with recent studies reporting the activation of zygotic transcription at MBT in Atlantic cod [30, 81]. As large fraction of genome is transcriptionally inactive, rearrangement of chromatin is essential to provide an access of transcription factors to the cis-regulatory elements . Microarray analyses are insufficient for accurate timing of the onset of transcription. Part of transcripts appear in the polyadenylated fraction due to maturation of maternal RNA [57, 63]. However, given modification of chromosomes during MBT and the size of the the SGM group, it is likely that a large part of mRNA denoted as zygotic was indeed transcribed from the zygotic genome. The SGM group was complex by composition and contained numerous developmental regulators. Massive upregulation of homeobox genes at SGM is consistent with their involvement in the establishment of body plan and formation of anterior-posterior axis of the embryo [83, 84]. Homeobox transcription factors and cell signaling pathways cooperate to pattern tissues and organs and to specify the fate of a variety of cell types. However, none of the functional groups and pathways was restricted to the post-MBT period and all were largely represented among UFO.
Transcriptome profiling of the oocytes and embryos of Atlantic cod with an aid of genome-wide microarray provided an insight in events taking place in early development and the roles of parental and zygotic transcripts. Maternal transcripts are involved in cellular metabolism, signal perception and transduction, defence, communication and contacts between cells. High representation of pathways and genes that control development suggest early cell fate specification and patterning of tissues and organs, especially of the neuronal lineage. The key event of zygotic genome activation at MBT was extensive chromatin rearrangements followed by expression of multiple developmental regulators.
The study was approved by the Norwegian Animal Research Authority and conducted according to the prevailing animal welfare regulations: FOR-1996-01-15-23 (Norway), European Convention for the Protection of Vertebrate Animals used for Experimental and Other Scientific Purposes (Strasbourg, 18.III.1986) and COUNCIL DIRECTIVE of 24 November 1986 on the approximation of laws, regulations and administrative provisions of the Member States regarding the protection of animals used for experimental and other scientific purposes (86/609/EEC).
Atlantic cod eggs and embryos were obtained from farmed fish at the National Cod Breeding Centre (Kraknes, Tromsø, Norway). Eggs were hand stripped, fertilized in vitro and transferred to seawater rearing tanks at an average temperature of 4.5°C and 100% oxygen saturation. The following stages were selected for analyses: 1) unfertilized (UFO) and 2) newly fertilized oocytes, 2 hpf and embryos at 3) mid-blastula (MBT), 4) 12 somites and 5) 52 somites (end of somitogenesis). Embryonic stages were determined based on description of Atlantic cod development with precise timing . Tissues from adult male and female cod were used as a reference in the microarray analyses. Eggs and tissue samples were stored in RNAlater (Ambion, Austin, Texas, USA).
Total RNA was extracted from Atlantic cod eggs and tissues using TRIzol (Life Technologies) and PureLinkTM RNA mini kit (Ambion, Austin, Texas, USA). For each developmental stage, 10 oocytes/embryos were pooled for the analyses. On-column DNase treatment was performed using PureLinkTM DNase (Life Technologies) in order to remove traces of DNA and impurities. The concentration was analyzed by NanoDrop ND-1000 spectrophotometer (Thermo Fisher Scientific, Wilmington, USA). The total RNA quality was assessed with Agilent 2100 Bioanalyzer (Agilent 2100 Bioanalyzer, Agilent Technologies, Waldbronn, Germany) and only the samples of high quality (RIN ≥ 8) were selected for analysis.
The Nofima’s Atlantic cod oligonucleotide microarray (ACIQ-2) produced by Agilent Technologies in the 4 × 44 k format included 60-mer probes to the unique transcripts from Ensembl and Unigene which were annotated by functional categories of GO and pathways of KEGG using bioinformatics package STARS [78, 86]. The genes were assigned to the orthology groups of OrthoDB . Three and two biological replicates of the respectively three first and two last stages were analyzed in a total of 13 microarrays. Reference RNA was prepared by pooling equal amounts of RNA from pyloric caeca, liver, muscle, brain and male and female gonad to identify genes with increased expression in oocytes and embryos or developmentally regulated genes. The common reference design also made possible comparison between stages and finding of stage-specific genes. RNA amplification, labeling and fragmentation were performed using Two-Colour Quick Amp Labeling Kit and Gene Expression Hybridization kit following the manufacturer's instructions (Agilent Technologies). The input of total RNA used in each reaction was 100 ng. Individual samples were compared to the common reference; assignment of fluorescent labels (Cy5 and Cy3) was changed in each hybridization performed at 65°C at the rotation speed of 10 rpm for 17 hours in the oven (Agilent Technologies). The slides were washed with Gene Expression Wash Buffers 1 and 2 as described by the manufacturer and scanning was performed at 5 μm resolution using a GenePix Personal 4100A scanner (Molecular Devices, Sunnyvale, CA, USA). The laser power was manually adjusted and the “auto PMT” was enabled to adjust PMT for each channel such that less than 0.1% of features were saturated and that the mean intensity ratio of the Cy3 and Cy5 signals was close to one. Nofima’s bioinformatic package STARS was used for data processing and mining. After filtration of low quality spots flagged by FE, lowess normalization of log2-expression ratios (ER) was performed. Results for the two last developmental stages were highly similar and these samples were therefore merged and denoted as SGM (segmentation). Features that passed quality control in all samples of at least one stage and showed over 2-fold difference from reference were selected (Additional file 1). Further, the features were assigned to groups with different temporal profiles (Figure 1) according to criteria presented in Table 1 with minor manual editing. Data were submitted to GEO Omnibus (GSE58392).
Quantitative real-time RT-PCR
Eight genes were selected for qPCR analyses based on the results of microarray analyses (Table 2). Gene expression measured by real-time qPCR was performed on the same samples that were used for microarray analysis. Primers were designed with Primer3 software and synthesized by Life Technologies. The amplicon lengths were set to be between 60 and 200 base pairs. The cDNA synthesis was performed on 1.5 μg total RNA using the SuperScript® VILO™ cDNA Synthesis Kit (Life Technologies) in a 20 μl reaction system according to the manufacturer's protocol. The specificity of PCR amplification was confirmed with melting curve analysis. Efficiency was checked from tenfold serial dilutions of cDNA for each primer pair. A 2 × SYBR® Green PCR Mastermix (Roche Diagnostics, Mannheim, Germany), 0.8 mM of each primer, and 4 μl of 1:10 diluted cDNA template were mixed in 12 μl reaction volumes. PCR was performed in duplicates in 96-well optical plates on Light Cycler 480 (Roche Diagnostics, Mannheim, Germany) under the following conditions: 95°C for 5 min (pre-incubation), 95°C for 5 s, 60°C for 15 s, 72°C for 15 s (amplification), followed by 95°C for 5 s and 65°C for 1 min (melting curve). 45 amplification cycles were performed. We tested some of the commonly used reference genes and they all showed substantial variation in expression levels between fertilized oocytes and somitogenesis in agreement with other embryonic studies of Atlantic cod . Hence, we selected six other candidate genes that showed stable expression on microarray throughout the development and tested them as reference genes. Finally, the combination of five most stable genes (Table 2) was used for normalization and the data are given as –ddCt values.
Lubzens E, Young G, Bobe J, Cerdà J: Oogenesis in teleosts: how fish eggs are formed. Gen Comp Endocrinol. 2010, 165: 367-389. 10.1016/j.ygcen.2009.05.022.
DePamphilis ML, Kaneko KJ, Vassilev A: Activation Of Zygotic gene expression in mammals. Gene expression at the beginning of animal development, advances in developmental biology and biochemistry series, Volume 12. Edited by: M.L. D, PM W. 2002, Netherlands: Elsevier Science, 55-84.
Schier AF: The maternal-zygotic transition: death and birth of RNAs. Science. 2007, 316: 406-407. 10.1126/science.1140693.
Mtango NR, Potireddy S, Latham KE: Oocyte quality and maternal control of development. Int Rev Cell Mol Biol. 2008, 268: 223-290.
Fischer BE, Wasbrough E, Meadows LA, Randlet O, Dorus S, Karr TL, Russell S: Conserved properties of Drosophila and human spermatozoal mRNA repertoires. Proc Biol Sci. 2012, 279: 2636-2644. 10.1098/rspb.2012.0153.
Martins RP, Krawetz SA: RNA in human sperm. Asian J Androl. 2005, 7: 115-120. 10.1111/j.1745-7262.2005.00048.x.
Rudel D, Sommer RJ: The evolution of developmental mechanisms. Dev Biol. 2003, 264: 15-37. 10.1016/S0012-1606(03)00353-1.
Lawrence PA, Levine M: Mosaic and regulative development: two faces of one coin. Curr Biol. 2006, 16: R236-R239. 10.1016/j.cub.2006.03.016.
De Renzis S, Elemento O, Tavazoie S, Wieschaus EF: Unmasking activation of the zygotic genome using chromosomal deletions in the Drosophila embryo. PLoS Biol. 2007, 5: e117-10.1371/journal.pbio.0050117.
Mathavan S, Lee SG, Mak A, Miller LD, Murthy KR, Govindarajan KR, Tong Y, Wu YL, Lam SH, Yang H, Ruan Y, Korzh V, Gong Z, Liu ET, Lufkin T: Transcriptome analysis of zebrafish embryogenesis using microarrays. PLoS Genet. 2005, 1: 260-276.
Wieschaus E: Embryonic Transcription and the Control of Developmental Pathways. Genetics. 1996, 142: 5-10.
Flach G, Johnson MH, Braude PR, Taylor RA, Bolton VN: The transition from maternal to embryonic control in the 2-cell mouse embryo. Embo J. 1982, 1: 681-686.
Braude P, Bolton V, Moore S: Human gene expression first occurs between the four- and eight-cell stages of preimplantation development. Nature. 1988, 332: 459-461. 10.1038/332459a0.
Hoffert KA, Anderson GB, Wildt DE, Roth TL: Transition from maternal to embryonic control of development in IVM/IVF domestic cat embryos. Mol Reprod Dev. 1997, 48: 208-215. 10.1002/(SICI)1098-2795(199710)48:2<208::AID-MRD8>3.0.CO;2-W.
Frei RE, Schultz GA, Church RB: Qualitative and quantitative changes in protein synthesis occur at the 8-16-cell stage of embryogenesis in the cow. J Reprod Fertil. 1989, 86: 637-641. 10.1530/jrf.0.0860637.
Kraeussling M, Wagner TU, Schartl M: Highly asynchronous and asymmetric cleavage divisions accompany early transcriptional activity in pre-blastula medaka embryos. PLoS One. 2011, 6: 7-
Kimmel CB, Ballard WW, Kimmel SR, Ullmann B, Schilling TF: Stages of embryonic development of the zebrafish. Dev Dyn. 1995, 203: 253-310. 10.1002/aja.1002030302.
Howe JA, Newport JW: A developmental timer regulates degradation of cyclin E1 at the midblastula transition during Xenopus embryogenesis. Proc Natl Acad Sci U S A. 1996, 93: 2060-2064. 10.1073/pnas.93.5.2060.
Sibon OC, Stevenson VA, Theurkauf WE: DNA-replication checkpoint control at the Drosophila midblastula transition. Nature. 1997, 388: 93-97. 10.1038/40439.
Tadros W, Lipshitz HD: The maternal-to-zygotic transition: a play in two acts. Development. 2009, 136: 3033-3042. 10.1242/dev.033183.
Solnica-Krezel L: Conserved patterns of cell movements during vertebrate gastrulation. Curr Biol. 2005, 15: R213-R228. 10.1016/j.cub.2005.03.016.
Stern CD: Vertebrate gastrulation. Curr Opin Genet Dev. 1992, 2: 556-561. 10.1016/S0959-437X(05)80171-6.
Seeger M, Tear G, Ferres-Marco D, Goodman CS: Mutations affecting growth cone guidance in Drosophila: genes necessary for guidance toward or away from the midline. Neuron. 1993, 10: 409-426. 10.1016/0896-6273(93)90330-T.
Wright AP, Fox AN, Johnson KG, Zinn K: Systematic screening of Drosophila deficiency mutations for embryonic phenotypes and orphan receptor ligands. PLoS One. 2010, 5: 0012288-10.1371/journal.pone.0012288.
Nusslein-Volhard C, Wieschaus E: Mutations affecting segment number and polarity in Drosophila. Nature. 1980, 287: 795-801. 10.1038/287795a0.
Amsterdam A, Hopkins N: Mutagenesis strategies in zebrafish for identifying genes involved in development and disease. Trends Genet. 2006, 22: 473-478. 10.1016/j.tig.2006.06.011.
Haffter P, Granato M, Brand M, Mullins MC, Hammerschmidt M, Kane DA, Odenthal J, van Eeden FJ, Jiang YJ, Heisenberg CP, Kelsh RN, Furutani-Seiki M, Vogelsang E, Beuchle D, Schach U, Fabian C, Nüsslein-Volhard C: The identification of genes with unique and essential functions in the development of the zebrafish, Danio rerio. Development. 1996, 123: 1-36.
Alestrom P, Holter JL, Nourizadeh-Lillabadi R: Zebrafish in functional genomics and aquatic biomedicine. Trends Biotechnol. 2006, 24: 15-21. 10.1016/j.tibtech.2005.11.004.
Star B, Nederbragt AJ, Jentoft S, Grimholt U, Malmstrom M, Gregers TF, Rounge TB, Paulsen J, Solbakken MH, Sharma A, Wetten OF, Lanzén A, Winer R, Knight J, Vogel JH, Aken B, Andersen O, Lagesen K, Tooming-Klunderud A, Edvardsen RB, Tina KG, Espelund M, Nepal C, Previti C, Karlsen BO, Moum T, Skage M, Berg PR, Gjøen T, Kuhl H, et al: The genome sequence of Atlantic cod reveals a unique immune system. Nature. 2011, 477: 207-210. 10.1038/nature10342.
Drivenes O, Taranger GL, Edvardsen RB: Gene expression profiling of Atlantic cod (Gadus morhua) embryogenesis using microarray. Mar Biotechnol. 2012, 14: 167-176. 10.1007/s10126-011-9399-y.
Hall TE, Smith P, Johnston IA: Stages of embryonic development in the Atlantic cod Gadus morhua. J Morphol. 2004, 259: 255-270. 10.1002/jmor.10222.
Wickramasinghe D, Becker S, Ernst MK, Resnick JL, Centanni JM, Tessarollo L, Grabel LB, Donovan PJ: Two CDC25 homologues are differentially expressed during mouse development. Development. 1995, 121: 2047-2056.
Edgar BA, Datar SA: Zygotic degradation of two maternal Cdc25 mRNAs terminates Drosophila's early cell cycle program. Genes Dev. 1996, 10: 1966-1977. 10.1101/gad.10.15.1966.
Tourret J, McKeon F: Tyrosine kinases wee1 and mik1 as effectors of DNA replication checkpoint control. Prog Cell Cycle Res. 1996, 2: 91-97.
Hung LY, Tang CJ, Tang TK: Protein 4.1 R-135 interacts with a novel centrosomal protein (CPAP) which is associated with the gamma-tubulin complex. Mol Cell Biol. 2000, 20: 7813-7825. 10.1128/MCB.20.20.7813-7825.2000.
Roche KC, Wiechens N, Owen-Hughes T, Perkins ND: The FHA domain protein SNIP1 is a regulator of the cell cycle and cyclin D1 expression. Oncogene. 2004, 23: 8185-8195. 10.1038/sj.onc.1208025.
Krishnan S, Horowitz S, Trievel RC: Structure and function of histone H3 lysine 9 methyltransferases and demethylases. Chembiochem. 2011, 12: 254-263. 10.1002/cbic.201000545.
Morris KV: RNA and the Regulation of Gene Expression: A Hidden Layer of Complexity. 2008, Norfolk: Caister Academic Press
Hanai K, Furuhashi H, Yamamoto T, Akasaka K, Hirose S: RSF governs silent chromatin formation via histone H2Av replacement. PLoS Genet. 2008, 4: 1000011-10.1371/journal.pgen.1000011.
Sinha A, Faller DV, Denis GV: Bromodomain analysis of Brd2-dependent transcriptional activation of cyclin A. Biochem J. 2005, 387: 257-269. 10.1042/BJ20041793.
Denis GV, McComb ME, Faller DV, Sinha A, Romesser PB, Costello CE: Identification of transcription complexes that contain the double bromodomain protein Brd2 and chromatin remodeling machines. J Proteome Res. 2006, 5: 502-511. 10.1021/pr050430u.
Vassetzky Y, Hair A, Mechali M: Rearrangement of chromatin domains during development in Xenopus. Genes Dev. 2000, 14: 1541-1552.
Schlingensiepen K-H, Wollnik F, Kunst M, Schlingensiepen R, Herdegen T, Brysch W: The role of Jun transcription factor expression and phosphorylation in neuronal differentiation, neuronal cell death, and plastic adaptationsin vivo. Cell Mol Neurobiol. 1994, 14: 487-505. 10.1007/BF02088833.
Jacobs-Helber SM, Abutin RM, Tian C, Bondurant M, Wickrema A, Sawyer ST: Role of JunB in erythroid differentiation. J Biol Chem. 2002, 277: 4859-4866. 10.1074/jbc.M107243200.
Bamforth SD, Braganca J, Eloranta JJ, Murdoch JN, Marques FI, Kranc KR, Farza H, Henderson DJ, Hurst HC, Bhattacharya S: Cardiac malformations, adrenal agenesis, neural crest defects and exencephaly in mice lacking Cited2, a new Tfap2 co-activator. Nat Genet. 2001, 29: 469-474. 10.1038/ng768.
Torgersen JS, Takle H, Andersen Ø: Differential spatial expression of mef2 paralogs during cardiac development in Atlantic cod (Gadus morhua). Comp Biochem Physiol B Biochem Mol Biol. 2011, 158: 181-187. 10.1016/j.cbpb.2010.11.006.
Kageyama R, Ohtsuka T, Kobayashi T: The Hes gene family: repressors and oscillators that orchestrate embryogenesis. Development. 2007, 134: 1243-1251. 10.1242/dev.000786.
Stafford DA, Monica SD, Harland RM: Follistatin interacts with Noggin in the development of the axial skeleton. Mech Dev. 2014, 134: 78-85.
Schweitzer J, Gimnopoulos D, Lieberoth BC, Pogoda HM, Feldner J, Ebert A, Schachner M, Becker T, Becker CG: Contactin1a expression is associated with oligodendrocyte differentiation and axonal regeneration in the central nervous system of zebrafish. Mol Cell Neurosci. 2007, 35: 194-207. 10.1016/j.mcn.2007.02.018.
Shimoda Y, Watanabe K: Contactins: emerging key roles in the development and function of the nervous system. Cell Adh Migr. 2009, 3: 64-70. 10.4161/cam.3.1.7764.
Racher H, Hansen D: PUF-8, a Pumilio homolog, inhibits the proliferative fate in the Caenorhabditis elegans germline. G3. 2012, 2: 1197-1205. 2012.
Thiery JP, Acloque H, Huang RY, Nieto MA: Epithelial-mesenchymal transitions in development and disease. Cell. 2009, 139: 871-890. 10.1016/j.cell.2009.11.007.
Wuebben EL, Mallanna SK, Cox JL, Rizzino A: Musashi2 is required for the self-renewal and pluripotency of embryonic stem cells. PLoS One. 2012, 7: 4-
Bree RT, McLoughlin S, Jin SW, McMeel OM, Stainier DY, Grealy M, Byrnes L: nanor, a novel zygotic gene, is expressed initially at the midblastula transition in zebrafish. Biochem Biophys Res Commun. 2005, 333: 722-728. 10.1016/j.bbrc.2005.05.168.
Landthaler M, Yalcin A, Tuschl T: The human DiGeorge syndrome critical region gene 8 and Its D. melanogaster homolog are required for miRNA biogenesis. Curr Biol. 2004, 14: 2162-2167. 10.1016/j.cub.2004.11.001.
Giraldez AJ, Mishima Y, Rihel J, Grocock RJ, Van Dongen S, Inoue K, Enright AJ, Schier AF: Zebrafish MiR-430 promotes deadenylation and clearance of maternal mRNAs. Science. 2006, 312: 75-79. 10.1126/science.1122689.
Aanes H, Winata CL, Lin CH, Chen JP, Srinivasan KG, Lee SG, Lim AY, Hajan HS, Collas P, Bourque G, Gong Z, Korzh V, Aleström P, Mathavan S: Zebrafish mRNA sequencing deciphers novelties in transcriptome dynamics during maternal to zygotic transition. Genome Res. 2011, 21: 1328-1338. 10.1101/gr.116012.110.
Avilion AA, Nicolis SK, Pevny LH, Perez L, Vivian N, Lovell-Badge R: Multipotent cell lineages in early mouse development depend on SOX2 function. Genes Dev. 2003, 17: 126-140. 10.1101/gad.224503.
Bonnet A, Bevilacqua C, Benne F, Bodin L, Cotinot C, Liaubet L, Sancristobal M, Sarry J, Terenina E, Martin P, Tosser-Klopp G, Mandon-Pepin B: Transcriptome profiling of sheep granulosa cells and oocytes during early follicular development obtained by laser capture microdissection. BMC Genomics. 2011, 12: 1471-2164.
Shen-Orr SS, Pilpel Y, Hunter CP: Composition and regulation of maternal and zygotic transcriptomes reflects species-specific reproductive mode. Genome Biol. 2010, 11: 2010-2011.
Gilbert I, Bissonnette N, Boissonneault G, Vallee M, Robert C: A molecular analysis of the population of mRNA in bovine spermatozoa. Reproduction. 2007, 133: 1073-1086. 10.1530/REP-06-0292.
Ostermeier GC, Miller D, Huntriss JD, Diamond MP, Krawetz SA: Reproductive biology: delivering spermatozoan RNA to the oocyte. Nature. 2004, 429: 154-
Harvey SA, Sealy I, Kettleborough R, Fenyes F, White R, Stemple D, Smith JC: Identification of the zebrafish maternal and paternal transcriptomes. Development. 2013, 140: 2703-2710. 10.1242/dev.095091.
Ellingsen T, Strand C, Monsen E, Bogwald J, Dalmo RA: The ontogeny of complement component C3 in the spotted wolffish (Anarhichas minor Olafsen). Fish Shellfish Immunol. 2005, 18: 351-358. 10.1016/j.fsi.2004.09.002.
Lovoll M, Kilvik T, Boshra H, Bogwald J, Sunyer JO, Dalmo RA: Maternal transfer of complement components C3-1, C3-3, C3-4, C4, C5, C7, Bf, and Df to offspring in rainbow trout (Oncorhynchus mykiss). Immunogenetics. 2006, 58: 168-179. 10.1007/s00251-006-0096-3.
Wang Z, Zhang S, Wang G, An Y: Complement activity in the egg cytosol of zebrafish Danio rerio: evidence for the defense role of maternal complement components. PLoS One. 2008, 3: 0001463-10.1371/journal.pone.0001463.
Wang Z, Zhang S, Tong Z, Li L, Wang G: Maternal transfer and protective role of the alternative complement components in zebrafish Danio rerio. PLoS One. 2009, 4: 18-
Seppola M, Johnsen H, Mennen S, Myrnes B, Tveiten H: Maternal transfer and transcriptional onset of immune genes during ontogenesis in Atlantic cod. Dev Comp Immunol. 2009, 33: 1205-1211. 10.1016/j.dci.2009.06.013.
Swain P, Nayak SK: Role of maternally derived immunity in fish. Fish Shellfish Immunol. 2009, 27: 89-99. 10.1016/j.fsi.2009.04.008.
Warga RM, Kane DA: A role for N-cadherin in mesodermal morphogenesis during gastrulation. Dev Biol. 2007, 310: 211-225. 10.1016/j.ydbio.2007.06.023.
Gumbiner BM: Regulation of cadherin-mediated adhesion in morphogenesis. Nat Rev Mol Cell Biol. 2005, 6: 622-634. 10.1038/nrm1699.
Kishigami S, Mishina Y: BMP signaling and early embryonic patterning. Cytokine Growth Factor Rev. 2005, 16: 265-278. 10.1016/j.cytogfr.2005.04.002.
Hayward P, Kalmar T, Arias AM: Wnt/Notch signalling and information processing during development. Development. 2008, 135: 411-424. 10.1242/dev.000505.
Hayflick JS, Wolfgang WJ, Forte MA, Thomas G: A unique Kex2-like endoprotease from Drosophila melanogaster is expressed in the central nervous system during early embryogenesis. J Neurosci. 1992, 12: 705-717.
Cadieux B, Chitramuthu BP, Baranowski D, Bennett HP: The zebrafish progranulin gene family and antisense transcripts. BMC Genomics. 2005, 6: 156-10.1186/1471-2164-6-156.
Masi T, Johnson AD: Axbrn-1: a maternal transcript encodes a POU transcription factor that is later expressed in the developing central nervous system of axolotl embryos. Dev Genes Evol. 2001, 211: 449-452. 10.1007/s004270100171.
Krasnov A, Timmerhaus G, Schiøtz BL, Torgersen J, Afanasyev S, Iliev D, Jørgensen J, Takle H, Jørgensen SM: Genomic survey of early responses to viruses in Atlantic salmon, Salmo salar L. Mol Immunol. 2011, 49: 163-174. 10.1016/j.molimm.2011.08.007.
Krasnov A, Kileng O, Skugor S, Jorgensen SM, Afanasyev S, Timmerhaus G, Sommer AI, Jensen I: Genomic analysis of the host response to nervous necrosis virus in Atlantic cod (Gadus morhua) brain. Mol Immunol. 2013, 54: 443-452. 10.1016/j.molimm.2013.01.010.
Alizadeh Z, Kageyama S, Aoki F: Degradation of maternal mRNA in mouse embryos: selective degradation of specific mRNAs after fertilization. Mol Reprod Dev. 2005, 72: 281-290. 10.1002/mrd.20340.
Ho L, Crabtree GR: Chromatin remodelling during development. Nature. 2010, 463: 474-484. 10.1038/nature08911.
Skjaerven KH, Olsvik PA, Finn RN, Holen E, Hamre K: Ontogenetic expression of maternal and zygotic genes in Atlantic cod embryos under ambient and thermally stressed conditions. Comp Biochem Physiol A Mol Integr Physiol. 2011, 159: 196-205. 10.1016/j.cbpa.2011.02.026.
Felsenfeld G: Chromatin unfolds. Cell. 1996, 86: 13-19. 10.1016/S0092-8674(00)80073-2.
Pearson JC, Lemons D, McGinnis W: Modulating Hox gene functions during animal body patterning. Nat Rev Genet. 2005, 6: 893-904. 10.1038/nrg1726.
Wellik DM: Hox genes and vertebrate axial pattern. Curr Top Dev Biol. 2009, 88: 257-278.
Gorodilov YN, Terjesen B, Krasnov A, Takle H: Description of Embryogenesis of Atlantic Cod Gadus morhua. Open Mar Biol J. 2008, 2: 43-53. 10.2174/1874450800802010043.
Krasnov A, Timmerhaus G, Afanasyev S, Jorgensen SM: Development and assessment of oligonucleotide microarrays for Atlantic salmon (Salmo salar L.). Comp Biochem Physiol Part D Genomics Proteomics. 2011, 6: 31-38. 10.1016/j.cbd.2010.04.006.
Waterhouse RM, Tegenfeldt F, Li J, Zdobnov EM, Kriventseva EV: OrthoDB: a hierarchical catalog of animal, fungal and bacterial orthologs. Nucleic Acids Res. 2013, 41: 24-10.1093/nar/gks904.
Johnsen H, Andersen Ø: Sex dimorphic expression of five dmrt genes identified in the Atlantic cod genome. The fish-specific dmrt2b diverged from dmrt2a before the fish whole-genome duplication. Gene. 2012, 505: 221-232. 10.1016/j.gene.2012.06.021.
We would like to thank the Norwegian Research Council (project no: 190371) for providing necessary financial support.
The authors declare that they have no competing interests.
ØA designed the experiment. AŠ performed the gene expression analyses. AK analyzed the data and together with AŠ drafted the manuscript. All authors read, contributed to and approved the final manuscript.
Electronic supplementary material
Additional file 1: Data are log2-Expression Ratio to reference. Column S - groups (prezygotic and zygotic), column T - subgroups defined by the temporal profiles. (XLSX 985 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Škugor, A., Krasnov, A. & Andersen, Ø. Genome-wide microarray analysis of Atlantic cod (Gadus morhua) oocyte and embryo. BMC Genomics 15, 594 (2014). https://doi.org/10.1186/1471-2164-15-594