Differential transcript isoform usage pre- and post-zygotic genome activation in zebrafish
© Aanes et al.; licensee BioMed Central Ltd. 2013
Received: 26 October 2012
Accepted: 25 April 2013
Published: 15 May 2013
Zebrafish embryos are transcriptionally silent until activation of the zygotic genome during the 10th cell cycle. Onset of transcription is followed by cellular and morphological changes involving cell speciation and gastrulation. Previous genome-wide surveys of transcriptional changes only assessed gene expression levels; however, recent studies have shown the necessity to map isoform-specific transcriptional changes. Here, we perform isoform discovery and quantification on transcriptome sequences from before and after zebrafish zygotic genome activation (ZGA).
We identify novel isoforms and isoform switches during ZGA for genes related to cell adhesion, pluripotency and DNA methylation. Isoform switching events include alternative splicing and changes in transcriptional start sites and in 3’ untranslated regions. New isoforms are identified even for well-characterized genes such as pou5f1, sall4 and dnmt1. Genes involved in cell-cell interactions such as f11r and magi1 display isoform switches with alterations of coding sequences. We also detect over 1000 transcripts that acquire a longer 3’ terminal exon when transcribed by the zygote compared to their maternal transcript counterparts. ChIP-sequencing data mapped onto skipped exon events reveal a correlation between histone H3K36 trimethylation peaks and skipped exons, suggesting epigenetic marks being part of alternative splicing regulation.
The novel isoforms and isoform switches reported here include regulators of transcriptional, cellular and morphological changes taking place around ZGA. Our data display an array of isoform-related functional changes and represent a valuable resource complementary to existing early embryo transcriptomes.
KeywordsZebrafish Mid-blastula Transition Zygotic Genome Activation Alternative Splicing Transcriptional Start Site 3’UTR
During the first ten cell cycles after fertilization, the zebrafish embryo is transcriptionally silent and consists of undifferentiated and rapidly dividing blastomeres . Initiation of transcription at the 10th cell cycle is termed zygotic genome activation (ZGA). At the time of ZGA, development proceeds from a control by mRNAs synthesized during oogenesis and stored in the egg (maternal transcripts) to mRNAs produced by the embryo’s own genome. Following ZGA, blastomeres divide less frequently and more asynchronously, they start to differentiate and migrate to form the three germ layers of the gastrulating embryo [1–4]. Collectively, these transformations characterize the mid-blastula transition (MBT). These fundamental cellular and functional changes occurring during development make zebrafish an attractive model to study transcriptional changes governing between pre-MBT and post-MBT development.
Previous studies have shown essential roles of activation and degradation of maternal transcripts in regulating the MBT and ZGA [5–7]. Signaling pathways involving Bmp, Nodal, Fgf, Wnt and maternal β-catenin are essential for the formation of germ layers and body axis . We and others have shown that the establishment of post-translationally modified histones on specific genomic sites and DNA methylation play a role in transcriptional regulation around the time of ZGA by patterning developmental gene expression [9–12]. However, although several early transcriptomes have recently been published [5, 13, 14], little is known on the isoform-specific dynamics governing developmental transitions around the MBT.
The notion that each gene can give rise to multiple mRNAs has evolved from being reported as a rare phenomenon [15, 16] to include virtually all loci in man [17, 18], and has been shown to be crucial for differentiation , development  and human disease . Isoform switches, defined as a change in the isoform composition of gene products between two conditions (e.g. two developmental stages) are of particular interest since such events are critical for differentiation . The differences between transcript isoforms can affect the coding sequence (CDS) and/or untranslated regions (UTRs) 3' or 5' of the CDS. The former is likely to affect protein function, while the latter may affect translational efficiency, mRNA degradation kinetics and spatial distribution of transcripts [22, 23]. The mechanisms regulating splicing and determining the production of specific transcript isoforms have started to be unveiled and involve cis-elements and trans-acting factors, as well as epigenetic modifications in the proximity of spliced exons [24–26]. A genome-wide landscape of transcript isoforms synthesized during early development, and particularly of the switches occurring between maternal and zygotic isoforms at the time the embryo initiates its own transcriptional program, has been lacking and has hampered a comprehensive appreciation of the transcriptional dynamics occurring at the time of ZGA.
We have used an isoform prediction and quantification program, Cufflinks , to detect novel isoforms and quantify isoform-specific changes from RNA-sequencing (RNA-seq) reads before and after ZGA in zebrafish. We identify numerous novel isoforms related to shifts in transcription start site (TSS), alternative splicing (AS) events and transcription termination sites (TTS), when comparing transcripts of maternal and zygotic origin coming from the same gene. These include transcripts of genes involved in cell-cell interactions, pluripotency control and DNA methylation. Using H3K4 and H3K36 trimethylation (me3) data acquired by chromatin immunoprecipitation and high-throughput sequencing (ChIP-seq), we find that H3K4me3 can form relatively broad domains which cannot distinguish between closely spaced alternative TSSs, unless TSSs are linked to alternative promoters at distant locations. Skipped exons are enriched in H3K36me3 and H3K4me3, extending previous reports on the involvement of these histone modifications in the regulation of isoform usage.
Results and discussion
Mapping of RNA-seq reads achieved from combined use of two aligners
Transcriptome assembly identifies novel genes and isoforms
To identify novel isoforms and genes, the predicted isoforms in the FPKM > 3 dataset were compared to Ensembl annotations (Zv9, version 64). This classifies 8,108 of the transcripts as identical, 12,100 as new isoforms, and 4,174 as novel transcribed regions (NTRs) i.e. possibly novel genes (Figure 2d). This analysis indicates that the assembled zebrafish transcriptome contains numerous novel transcripts, some of which may be encoded by novel genes.
Functional differences between annotated and novel transcripts
The large number of novel isoforms warranted further validation and characterization. We investigated the coding potential of new isoforms and found that nearly all have an open reading frame (ORF) (11,627/12,100; 96%) of >50 amino acids in length (starting with an AUG start codon and ending with a stop codon). Increasing the ORF length requirement results in a slight decrease in the number of ORFs detected (100 amino acids, 89%; 150 amino acids, 81%).
We then compared both the transcript sequence and the inferred coding sequence between novel isoforms and their closest matching annotated transcript from the same locus (as defined by Cufflinks). On average, the novel isoforms have shorter coding sequences than the corresponding annotated coding sequences with a mean difference of −230 bp (median −109 bp). Of 11,529 comparisons, 1,695 of the novel isoforms differ by one amino acid or less, indicating identical coding sequences (this corresponds to 15% of the transcripts with an ORF of ≥50 amino acids). Transcript isoforms with no difference in coding sequences length have on average longer 5’ and/or 3’ UTRs, with a mean length difference of 549 bp (median 256 bp). Therefore, the novel isoforms can be broadly divided into two groups, those with no change in coding sequence length but with longer 5’ and/or 3’ UTR, and those with an alteration in coding sequence length.
We next examined any functional effect these transcript isoforms may have, and used the predicted protein sequences for a domain scan using Pfam, a database of protein families . The detection of a protein domain has been shown to be suggestive of protein functionality . Of the novel isoforms analyzed, 8,949 (74%) contain a sequence with high similarity to an identified protein domain. In comparison, the closest matching annotated transcripts (n = 5,908) have a substantially larger fraction of Pfam domain-containing sequences (5,196; 88%). Further, most of the novel isoforms encode the same functional domains as the annotated isoforms (Figure 2e). However, a substantial fraction loses either all functional domains (n = 1971), or some of them (n = 2017; Figure 2e). Also, a small subgroup (n = 123) of the novel isoforms contains functional domains absent from the annotated isoform (Figure 2e). Taken together, these results suggest that the majority of the novel transcript isoforms are functional, with roles related to their annotated counterparts. However, some lose or gain functional domains, suggesting that they might exhibit different biological functions.
TSS usage and switching
Among the 513 genes which display a change in TSS pattern between the pre- and post-ZGA, dazl represents a subset where the TSSs are defined by different first exons. This type of TSS switch is likely to affect developmental functions as illustrated with dazls changes in the 5’UTR and N-terminus of the protein.
Epigenetic regulation of TSS usage and expression
H3K4me3 enrichment over the TSS is considered as a histone modification associated with transcription permissiveness. While virtually all expressed genes are marked by H3K4me3 at the TSS [39, 40], many genes that are not expressed also harbor H3K4me3 [9, 41]. H3K27me3 promoter occupancy is on the other hand associated with transcriptional repression . To identify a relationship between TSS usage and potentially associated chromatin marks, we used published ChIP-seq peak calling data for H3K4me3 and H3K27me3 in post-MBT stage embryos , and investigated the correlation between these histone modifications and the TSSs mapped in our study.
We find that a subset of 22,534 TSSs from 13,542 genes from our ‘all’ gene annotation dataset are enriched in H3K4me3. Most peaks (18,547/23,778; 78%) were within ±500 bp of a gene. The number of TSSs in the dataset and the number of peaks covering them indicates that multiple TSSs are included within each H3K4me3 peak, as exemplified by dazl and sall4 (Figure S4a,b in Additional file 1). On rare occasions, we observe that H3K4me3 peaks are located over the most frequently used TSS (major TSS) (e.g. actl6a; Figure S4c in Additional file 1). Therefore, at present H3K4me3 peaks cannot be used to predict the major TSS since the mark can overlap multiple TSSs. More accurate determination of TSS usage, e.g. using approaches such as cap analysis of gene expression (CAGE) together with H3K4me3 enrichment data may improve the usefulness of this histone modification in predicting major active start sites among closely located TSSs.
Our RNA-seq data show that 21% of the H3K4me3-marked genes are expressed at very low levels, with 2,822 loci displaying an FPKM value below 1 in both pre- and post-ZGA samples (Additional file 3). For the subgroup of bidirectional promoters identified using the H3K4me3 peaks (n = 1,433 promoter pairs) (see Methods) we find no significant correlation in expression between the bidirectional gene pairs (Pearson’s product–moment correlation <0.01), again supporting H3K4me3 as a mark for predisposition of promoter activation rather than having a direct activating effect, even when presumably sharing the same promoter. However, on average, H3K4me3-marked genes are more strongly expressed than genes without H3K4me3, as described earlier for RefSeq genes (data not shown) .
An important question is whether epigenetic marks can regulate particular subgroups of genes, functionally and spatially. To address the first question, analysis in DAVID of gene ontology (GO) terms for biological functions were obtained using H3K4me3- and H3K27me3-marked genes. This was significant only for one term for the H3K4me3-marked genes; cell cycle (FDR < 0.001; Additional file 6). This is in contrast to H3K27me3-marked genes which were enriched for 21 GO terms with the same cut off (i.e. regulation of transcription, embryonic morphogenesis, pattern specification process and regionalization; Additional file 7). Thus, whereas H3K4me3 mark TSSs of genes linked to many biological functions, H3K27me3 marks mostly the TSSs of genes involved in developmental functions. These results are consistent with a H3K4me3 pre-patterning of promoters of a specific subset of genes pre-MBT (involved in housekeeping and homeostatic functions), and with a view of H3K27me3 repressing developmental genes prior to their activation at or after the MBT [9, 10].
Another intriguing question is whether TSS marking with both H3K4me3 and H3K27me3 exists (so-called ‘bivalency’). We utilized the extensive annotation database of where zebrafish genes are expressed anatomically (anatomical ontology provided by the zebrafish model organism database; ZFIN ) to bring further clarity to this issue. H3K4me3-marked genes are mostly expressed in the entire organism, in contrast to H3K27me3-marked genes which are predominantly annotated with expression in distinct tissue types (in particular central nervous system structures; Additional file 8 and Additional file 9). This is consistent with H3K27me3 analysis of dissected Xenopus embryos , and not supportive of bivalency in zebrafish embryos, as others have suggested .
Taken together, our results show that during zebrafish development, H3K4me3 is a mark of TSSs independent of transcription status. The marked genes show no pronounced tissue-specific preferences and few specialized biological function as a group. In contrast, H3K27me3 marked genes have specific biological functions and are expressed in specialized cells. Our results are in line with the view that H3K4me3 contributes to a transcriptionally permissive chromatin environment for TSS usage; however, the repressive marks seem to contribute to lineage commitment and tissue-specific differentiation . The use of ChIP-seq data and a comprehensive gene annotation based on existing and novel annotations confirm the trends from previous studies based on older annotation and ChIP-chip data [9, 10].
Shifts in alternative splicing during the switch from maternal to zygotic transcripts
Alternative splicing (AS) allows the production of multiple mature mRNAs from the same primary transcript, increasing genetic diversity of a genome by several-fold. We used Cufflinks to identify shifts in AS isoform patterns before and after the ZGA. A total of 2,160 primary transcripts exhibited a significant change in the splice isoform composition (examples include foxh1, hdac5, dnmt3 and gsk3a; Additional file 10). Quality control by visual inspection revealed mostly subtle changes among the significantly changing primary transcripts. Importantly we find, also by visual inspection, that several obvious isoform switches do not pass the significance criteria of Cufflinks. For example, we report f11r, encoding a trans-membrane adhesion protein found in association with tight junctions [46, 47], which displays a switch between the maternal and zygotic stages, with exon 8 skipped only in the zygotic transcripts (Figure 3b). Using reverse transcription (RT)-qPCR, we validated the expression pattern of f11r (Figure S5 in Additional file 1). Exon 8 of f11r is only 13 bp long and is included in the maternal isoform of the f11r transcript (Figure 3c, upper panel). Skipping of exon 8, detected in the zygotic transcriptome (Figure 3c, lower panel), causes a frame-shift with an alteration of the C-terminal 36 amino acids; however two known functional domains in the F11r protein are not directly affected, although indirect structural effects cannot be excluded (Figure 3c). The change in the intracellular C-terminal portion of the protein includes the removal of two terminal valine residues which are involved in binding to PDZ-domain containing proteins . In zebrafish, interference with F11r function causes defects in migration at later stages of development . Interestingly, the shift in f11r isoform usage between pre- and post-MBT stages coincides with the initiation of cell migration taking place after the MBT [1, 2]. It is tempting to speculate that attenuation of binding potential to PDZ domain-containing proteins enables cells to alter their migration patterns [49, 50]. However, functional studies are required to test this hypothesis.
From our dataset, it appears therefore that thousands of splice isoform switches occur between the maternal and zygotic transcriptomes; nonetheless, it is important to emphasize that the majority constitutes subtle changes rather than obvious functional on/off events.
Quantification of AS events
An example of the latter includes dnmt1, encoding DNA methyltransferase 1 responsible for maintenance of DNA methylation patterns, particularly during development (reviewed in ). dnmt1 has approximately 34 exons with three isoforms annotated in zebrafish; however, our analysis detects several additional isoforms (Figure S6 in Additional file 1). Cufflinks predicts 9 dnmt1 isoforms (FPKM > 3), but without full-length sequence data the number of biologically relevant isoforms remains uncertain. One new isoform of dnmt1 is robustly expressed maternally (Additional file 3) and truncated due to skipping of 19 exons (validated by cDNA cloning and sequencing; data shown in Additional file 11). This skipping event removes most of the functional domains involved in DNA methyltransferase activity, but it retains a DMAP1 binding domain (Figure 4b). This domain binds the co-repressor DMAP1 and thereby confers DNMT1 with transcriptional repressor activity [53, 54]. Interestingly, in Xenopus it has been shown that Dnmt1 is needed to repress transcription prior to the MBT  and this function does not rely on the cytosine methyltransferase catalytic activity of Dnmt1 . Zebrafish may therefore be a valuable candidate to shed more light on the functional role of dnmt1 and its isoforms. The novel dnmt1 isoforms discovered in this study may bring clarity to the enigmatic issue of transcriptional repression during cleavage stages in zebrafish [11, 57] and other species .
For the pluripotency-associated gene pou5f1 (also known as oct4) we identify several new splicing events, all of which were confirmed by cDNA cloning and sequencing (sequences are shown in Additional file 12). In zebrafish, there is only one annotated pou5f1 transcript, while in human three arise from AS . We observe novel pou5f1 AS events originating from alternative acceptor sites at exons 2, 3 and 5 (Figure 4c; Figure S7 in Additional file 1). The first event leads to deletion of a nucleotide triplet coding for glutamic acid (validated by one cDNA clone sequence out of 12 random clones analyzed), yet with no effect on the functional domains. The removal of a charged amino acid is however likely to affect the three-dimensional structure of Pou5f1, although functional consequences remain unknown at this stage. The second AS event is located 19 nt prior to the annotated exon 3 (validated by one cDNA clone sequence). This leads to a frame-shift and a truncated Pou5f1 protein, removing both functional domains. Lastly, a 4 nt deletion causes a frame-shift in the exon 5 (validated by 5 of 12 cDNA clones sequences, suggesting higher abundance of this isoform) and leads to a truncation of the homeobox domain. Of note, we also find novel isoforms for other genes important for pluripotency, such as klf4, sall4 and nanog. This indicates that several well-studied genes with central functions in pluripotency are still not fully annotated in zebrafish, and highlights a likelihood of uncovering novel intricacies about pluripotency control mechanisms.
Epigenetic marks associated with exon skipping
AS predominantly occurs co-transcriptionally, which may influence the regulation of splicing by changes of epigenetic signatures such as histone modifications and DNA methylation [24, 25, 60]. We used ChIP-seq to map H3K36me3 occupancy in post-ZGA embryos, and used published peak calling data for H3K4me3 and H3K27me3  to investigate any association between these modified histones and single ES events. We extracted the number of reads spanning each skipped exon (this is direct evidence for skipping of an exon; see Methods); 1,668 single ES events have one or more reads supporting the skipping event, of which most only with a few reads spanning them. Focusing on 446 exons with >5 reads supporting the skipping of these exons at the post-MBT stage, we examined peaks of H3K4, H3K27 and H3K36 trimethylation in their vicinity (see Methods). H3K36me3 peaks and skipped exons co-localize in 194 of the 446 exons (some exons overlap with multiple H3K36me3 peaks), which represent a slight overrepresentation relative to a random set of regions (~145 overlaps on average when 446 exons were drawn at random from all exons a thousand times; Chi-squared test, p < 0.001). H3K4me3 is also slightly overrepresented (113 overlaps versus 71 at random, p < 0.001), whereas H3K27me3 is not overrepresented over skipped exons (~35 peaks overlapped both skipped and non-skipped exons).
Longer 3’UTRs in zygotic than in maternal transcripts
In addition to the extension of last exons, we also observe alternative last exon usage, which may lead to changes in the CDS, the 3’UTR or both. An example of the latter is apparent in magi1, encoding a PDZ domain protein; due to a switch in last exon usage between pre- and post-ZGA, both the CDS and the 3’UTR are changed (Figure 6b; see also RT-qPCR validation in Figure S9 in Additional file 1). The 3’UTR of the post-ZGA-specific isoform contains three additional sequence domains absent in the maternal version, namely an internal ribosome entry site (IRES), a sex-lethal binding site (SXL-BS) and a K-box element. The latter two are known to often co-occur to create spatial and temporal expression patterns through negative regulation of translation [70, 71]. In addition, although the PDZ domains are not directly affected, alteration in protein sequence might affect the protein 3D-structure and in particular the fifth PDZ domain located proximal to the C-terminus . Interestingly, Magi1 interacts with β-catenin , another regulator of cell adhesion crucial for left-right asymmetry of the embryo [8, 74]. Each PDZ domain often displays distinct protein specificity and strikingly, β-catenin binding to MAGI1 depends on the 5th PDZ domain, thus potentially linking the isoform switch reported in our study to the regulation of Wnt signaling through β-catenin .
The transcriptomes assembled in this study contain a high degree of complexity pertaining to TSSs, AS and TTSs. Most of the described novel isoforms encode ORFs with predicted significant functional protein domains, suggesting functional and valid transcripts with a role in development. We report several loci which display a switch-like behavior in isoforms expressed pre-ZGA and post-ZGA stages; these notably include dazl (TSS switch), f11r (exon skipping) and magi1 (alternative last exon usage). Their changes in expression pattern coincide with key developmental processes such as onset of transcription upon ZGA, cell speciation and migration, suggesting functional roles for these isoforms. We also identify validated novel isoforms for dnmt1 and pou5f1, important genes in development and pluripotency control. Moreover, we find increased length of the 3’ terminal exon and 3’UTR after the ZGA in over 1,000 transcripts. Analysis of histone marks suggests that H3K4me3 has little value in marking the major TSS in loci with several start sites. However, we find higher levels of H3K36me3 over skipped exons relative to all exons, suggesting a role of H3K36me3 in zebrafish AS control.
Bioinformatic analyses relied heavily on custom python scripts and the use of publically available tools such as SAMtools  and BEDtools , as well as general data handling in R . Python scripts critical for the analyses are available as supplementary material (Additional file 14). Data were visualized in the Integrative Genomics viewer .
RNA-seq mapping and transcriptome assembly
Embryo collection and RNA sequencing was performed as previously described . Raw data files are available from NCBI’s short read archive (GSE22830). Reads were mapped using Bioscope 1.3 (with default settings) and TopHat 1.3.1 (with default options except: -a = 5, -I = 50, -F = 0). The resulting SAM-files were merged and a python script (best_alignment.py) used to ensure that a read were only present once. We performed transcriptome assembly using Cufflinks (version 1.2.1), with Ensembl annotation as guide (Zv9, annotation version 64), using default settings except; F = 0.0, -j = 0.05, -A = 0.08, --3-overhang-tolerance = 100. Additional programs in the Cufflinks pipeline were used to merge files, compare transcripts to known annotation and detect differences between the stages . We created additional datasets using scripts; FPKM > 3 using extract_gtf.py, TSS’s using get_tss_bed.py and TTS using get_tts.py.
We extracted transcript sequences using gffread in Cufflinks. These were subjected to ORF prediction using getORF (with settings: -minsize 150, -sreverse_sequence No, -find 1) from the EMBOSS suite . The resulting protein sequences were processed by scripts and R handling to retrieve the longest predicted ORF. We used Pfam to identify functional domains (Pfam A)  and domains considered significant by Pfam were used. A script (‘domain_compare.py’) was used to identify if domains were retained, lost or gained relative to the closest matching transcript. Detection of motifs in the 3’UTR was performed using UTRScan .
We used ASTALAVISTA version 2.2  to quantify alternative splicing events. This program takes as input a GTF file (FPKM > 3) and outputs all splicing events in the annotation together with a code representing type of event (exon skipping, intron retention, alternative acceptor site, alternative donor site). We used python scripts to extract and count events from the ASTALAVISTA output (summary.py and events.py).
Detection of length changes in last exon
where p1 and p2 are the proportions of reads in the last part of the exon relative to the first part pre- and post-ZGA, respectively. Exons with a shift towards more reads in the last part post-ZGA will have a positive value. We weighted this value using the log-transformed average of the exons FPKM values pre- and post-ZGA, assuming that higher coverage gives a more robust estimate of change.
Analysis of ChIP-seq data
Embryos for ChIP-seq were collected at 5.3 hpf and chromatin prepared as described . In accordance with institutional, national and international guidelines for early stage (0–5.3 hpf) zebrafish embryos no ethics committee approval was needed for any of the experiments performed in this study. Immunoprecipitation was performed using H3K36me3 (Diagenode 058–050, Denville, NJ) and ChIP enriched DNA and input samples prepared for sequencing according to Illumina standard protocol (#11257047 Rev.A). ChIP-seq reads were mapped using BWA  and peaks detected using CCAT 3.0 (Settings: fragment size = 200, sliding window size = 500, minimum score = 3.0) . Intron and exons from single exon skipping events identified by ASTALAVISTA were extracted with a script (extract_skipped_exon.py) and BEDtools used to count the number of reads supporting each event. Events with >5 reads post-ZGA were considered further. We created databases of all introns and exons using the ‘all’ annotation. The number of overlapping peaks and skipped exons (>1 bp overlap between peak and exon skipping event: from start upstream intron to stop downstream intron) were calculated using BEDtools and the random control using a script (bootstr_rand.py). We utilized Repitools to construct metagens . H3K4me3 and H3K27me3 peak data were from a recently published paper (Pauli et al., 2012). We identified potential bidirectional promoters using ChIPpeakAnno  with the H3K4me3 data  and the ‘all’ gene annotation file. The TSS of either gene in the pair had to be within 1000 bp of the middle of the H3K4me3 peak.
PCR and cloning
Newly assembled and annotated isoforms were validated by RT-(q)PCR. RNA from pools of embryos at different stages was isolated using TRIzol reagent followed by RNA cleanup using Qiagen RNeasy MiniKit (Hilden, Germany) as previously described . Control RNA of kanamycin (Promega #C1381, Madison, Wisconsin) was added prior to the RNA extraction and used as the reference control by qPCR. RNA was reverse transcribed using iScript Select cDNA synthesis Kit (BioRad #170-8896, Hercules, CA) according to the manufactures instructions. Primers were designed to cover specific sequences resulting from alternative splicing (f11r), alternative TSS (dazl), and alternative last exon usage (magi1) (for details see Supplementary Table 1 in Additional file 15). The dynamics of the isoforms was confirmed by RT-qPCR. qPCR was performed with the iCycler MyiQ real time PCR detection system and SYBR Green (BioRad, Hercules, CA). Primers pairs gave no signal in PCRs lacking template (data not shown). Relative expression was determined by the ∆∆-CT method. For pou5f1, both primers were designed in the sequence common for the predicted isoforms (Supplementary Table 1 in Additional file 15). Presence of the isoforms at 3.5 hpf was confirmed by TOPO TA cloning (Invitrogen K45-0001, Carlsberg, CA) and subsequent sequencing of randomly selected clones.
ChIP-seq data for H3K36me3 and input are available under [NCBI GEO GSE40629]. Our RNA-seq data are available under [NCBI GEO GSE22830].
Chromatin immunoprecipitation sequencing
Zygotic genome activation
Transcription start site
Transcription termination site
Fragments per kilobase per million mapped fragments
Open reading frame
False discovery rate
Alternative donor site
Alternative acceptor site
We thank Torbjørn Rognes and Andrew Reiner for valuable discussions on bioinformatic analyses, and Timothy Hughes, Susanne Lorenz and Leonardo Meza-Zapeda for assistance on ChIP-seq and mapping of ChIP-seq reads. This work was supported by the Research Council of Norway, the Norwegian Center for Stem Cell Research and Carlsberg foundation (OØ). HA has a PhD fellowship at the Norwegian School of Veterinary Science.
- Kane DA, Kimmel CB: The zebrafish midblastula transition. Development. 1993, 119: 447-456.PubMedGoogle Scholar
- Keller PJ, Schmidt AD, Wittbrodt J, Stelzer EH: Reconstruction of zebrafish early embryonic development by scanned light sheet microscopy. Science. 2008, 322: 1065-1069. 10.1126/science.1162493.View ArticlePubMedGoogle Scholar
- Ho RK, Kimmel CB: Commitment of cell fate in the early zebrafish embryo. Science. 1993, 261: 109-111. 10.1126/science.8316841.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Aanes H, Winata CL, Lin CH, Chen JP, Srinivasan KG, Lee SG, et al: Zebrafish mRNA sequencing deciphers novelties in transcriptome dynamics during maternal to zygotic transition. Genome Res. 2011, 21: 1328-1338. 10.1101/gr.116012.110.PubMed CentralView ArticlePubMedGoogle Scholar
- Giraldez AJ, Mishima Y, Rihel J, Grocock RJ, Van DS, Inoue K, et al: Zebrafish MiR-430 promotes deadenylation and clearance of maternal mRNAs. Science. 2006, 312: 75-79. 10.1126/science.1122689.View ArticlePubMedGoogle Scholar
- Mathavan S, Lee SG, Mak A, Miller LD, Murthy KR, Govindarajan KR, et al: Transcriptome analysis of zebrafish embryogenesis using microarrays. PLoS Genet. 2005, 1: 260-276.View ArticlePubMedGoogle Scholar
- Chan TM, Longabaugh W, Bolouri H, Chen HL, Tseng WF, Chao CH, et al: Developmental gene regulatory networks in the zebrafish embryo. Biochim Biophys Acta. 2009, 1789: 279-298. 10.1016/j.bbagrm.2008.09.005.View ArticlePubMedGoogle Scholar
- Vastenhouw NL, Zhang Y, Woods IG, Imam F, Regev A, Liu XS, et al: Chromatin signature of embryonic pluripotency is established during genome activation. Nature. 2010, 464: 922-926. 10.1038/nature08866.PubMed CentralView ArticlePubMedGoogle Scholar
- Lindeman LC, Andersen IS, Reiner AH, Li N, Aanes H, Ostrup O, et al: Prepatterning of Developmental Gene Expression by Modified Histones before Zygotic Genome Activation. Dev Cell. 2011, 21: 993-1004. 10.1016/j.devcel.2011.10.008.View ArticlePubMedGoogle Scholar
- Andersen IS, Reiner AH, Aanes H, Alestrom P, Collas P: Developmental features of DNA methylation during activation of the embryonic zebrafish genome. Genome Biol. 2012, 13: R65-10.1186/gb-2012-13-7-r65.PubMed CentralView ArticlePubMedGoogle Scholar
- Andersen IS, Ostrup O, Lindeman LC, Aanes H, Reiner AH, Mathavan S, et al: Epigenetic complexity during the zebrafish mid-blastula transition. Biochem Biophys Res Commun. 2012, 417: 1139-1144. 10.1016/j.bbrc.2011.12.077.View ArticlePubMedGoogle Scholar
- Pauli A, Valen E, Lin MF, Garber M, Vastenhouw NL, Levin JZ, et al: Systematic identification of long noncoding RNAs expressed during zebrafish embryogenesis. Genome Res. 2012, 22: 577-591. 10.1101/gr.133009.111.PubMed CentralView ArticlePubMedGoogle Scholar
- Vesterlund L, Jiao H, Unneberg P, Hovatta O, Kere J: The zebrafish transcriptome during early development. BMC Dev Biol. 2011, 11: 30-10.1186/1471-213X-11-30.PubMed CentralView ArticlePubMedGoogle Scholar
- Rosenfeld MG, Amara SG, Roos BA, Ong ES, Evans RM: Altered expression of the calcitonin gene associated with RNA polymorphism. Nature. 1981, 290: 63-65. 10.1038/290063a0.View ArticlePubMedGoogle Scholar
- Rosenfeld MG, Lin CR, Amara SG, Stolarsky L, Roos BA, Ong ES, et al: Calcitonin mRNA polymorphism: peptide switching associated with alternative RNA splicing events. Proc Natl Acad Sci U S A. 1982, 79: 1717-1721. 10.1073/pnas.79.6.1717.PubMed CentralView ArticlePubMedGoogle Scholar
- Pan Q, Shai O, Lee LJ, Frey BJ, Blencowe BJ: Deep surveying of alternative splicing complexity in the human transcriptome by high-throughput sequencing. Nat Genet. 2008, 40: 1413-1415. 10.1038/ng.259.View ArticlePubMedGoogle Scholar
- Wang ET, Sandberg R, Luo S, Khrebtukova I, Zhang L, Mayr C, et al: Alternative isoform regulation in human tissue transcriptomes. Nature. 2008, 456: 470-476. 10.1038/nature07509.PubMed CentralView ArticlePubMedGoogle Scholar
- Gabut M, Samavarchi-Tehrani P, Wang X, Slobodeniuc V, O'Hanlon D, Sung HK, et al: An alternative splicing switch regulates embryonic stem cell pluripotency and reprogramming. Cell. 2011, 147: 132-146. 10.1016/j.cell.2011.08.023.View ArticlePubMedGoogle Scholar
- Ramani AK, Calarco JA, Pan Q, Mavandadi S, Wang Y, Nelson AC, et al: Genome-wide analysis of alternative splicing in Caenorhabditis elegans. Genome Res. 2011, 21: 342-348. 10.1101/gr.114645.110.PubMed CentralView ArticlePubMedGoogle Scholar
- Cooper TA, Wan L, Dreyfuss G: RNA and disease. Cell. 2009, 136: 777-793. 10.1016/j.cell.2009.02.011.PubMed CentralView ArticlePubMedGoogle Scholar
- Pickering BM, Willis AE: The implications of structured 5' untranslated regions on translation and disease. Semin Cell Dev Biol. 2005, 16: 39-47. 10.1016/j.semcdb.2004.11.006.View ArticlePubMedGoogle Scholar
- Andreassi C, Riccio A: To localize or not to localize: mRNA fate is in 3'UTR ends. Trends Cell Biol. 2009, 19: 465-474. 10.1016/j.tcb.2009.06.001.View ArticlePubMedGoogle Scholar
- Luco RF, Pan Q, Tominaga K, Blencowe BJ, Pereira-Smith OM, Misteli T: Regulation of alternative splicing by histone modifications. Science. 2010, 327: 996-1000. 10.1126/science.1184208.PubMed CentralView ArticlePubMedGoogle Scholar
- Shukla S, Kavak E, Gregory M, Imashimizu M, Shutinoski B, Kashlev M, et al: CTCF-promoted RNA polymerase II pausing links DNA methylation to splicing. Nature. 2011, 479: 74-79. 10.1038/nature10442.View ArticlePubMedGoogle Scholar
- Barash Y, Calarco JA, Gao W, Pan Q, Wang X, Shai O, et al: Deciphering the splicing code. Nature. 2010, 465: 53-59. 10.1038/nature09000.View ArticlePubMedGoogle Scholar
- Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, et al: Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010, 28: 511-515. 10.1038/nbt.1621.PubMed CentralView ArticlePubMedGoogle Scholar
- Trapnell C, Pachter L, Salzberg SL: TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009, 25: 1105-1111. 10.1093/bioinformatics/btp120.PubMed CentralView ArticlePubMedGoogle Scholar
- Roberts A, Trapnell C, Donaghey J, Rinn JL, Pachter L: Improving RNA-Seq expression estimates by correcting for fragment bias. Genome Biol. 2011, 12: R22-10.1186/gb-2011-12-3-r22.PubMed CentralView ArticlePubMedGoogle Scholar
- Roberts A, Pimentel H, Trapnell C, Pachter L: Identification of novel transcripts in annotated genomes using RNA-Seq. Bioinformatics. 2011, 27: 2325-2329. 10.1093/bioinformatics/btr355.View ArticlePubMedGoogle Scholar
- Punta M, Coggill PC, Eberhardt RY, Mistry J, Tate J, Boursnell C, et al: The Pfam protein families database. Nucleic Acids Res. 2012, 40: D290-D301. 10.1093/nar/gkr1065.PubMed CentralView ArticlePubMedGoogle Scholar
- Leoni G, Le PL, Ferre F, Raimondo D, Tramontano A: Coding potential of the products of alternative splicing in human. Genome Biol. 2011, 12: R9-10.1186/gb-2011-12-1-r9.PubMed CentralView ArticlePubMedGoogle Scholar
- Carninci P, Sandelin A, Lenhard B, Katayama S, Shimokawa K, Ponjavic J, et al: Genome-wide analysis of mammalian promoter architecture and evolution. Nat Genet. 2006, 38: 626-635. 10.1038/ng1789.View ArticlePubMedGoogle Scholar
- Xu EY, Moore FL, Pera RA: A gene family required for human germ cell development evolved from an ancient meiotic gene conserved in metazoans. Proc Natl Acad Sci U S A. 2001, 98: 7414-7419. 10.1073/pnas.131090498.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen J, Melton C, Suh N, Oh JS, Horner K, Xie F, et al: Genome-wide analysis of translation reveals a critical role for deleted in azoospermia-like (Dazl) at the oocyte-to-zygote transition. Genes Dev. 2011, 25: 755-766. 10.1101/gad.2028911.PubMed CentralView ArticlePubMedGoogle Scholar
- Haston KM, Tung JY, Reijo Pera RA: Dazl functions in maintenance of pluripotency and genetic and epigenetic programs of differentiation in mouse primordial germ cells in vivo and in vitro. PLoS One. 2009, 4: e5654-10.1371/journal.pone.0005654.PubMed CentralView ArticlePubMedGoogle Scholar
- Collier B, Gorgoni B, Loveridge C, Cooke HJ, Gray NK: The DAZL family proteins are PABP-binding proteins that regulate translation in germ cells. EMBO J. 2005, 24: 2656-2666. 10.1038/sj.emboj.7600738.PubMed CentralView ArticlePubMedGoogle Scholar
- Takeda Y, Mishima Y, Fujiwara T, Sakamoto H, Inoue K: DAZL relieves miRNA-mediated repression of germline mRNAs by controlling poly(A) tail length in zebrafish. PLoS One. 2009, 4: e7513-10.1371/journal.pone.0007513.PubMed CentralView ArticlePubMedGoogle Scholar
- Bernstein BE, Humphrey EL, Erlich RL, Schneider R, Bouman P, Liu JS, et al: Methylation of histone H3 Lys 4 in coding regions of active genes. Proc Natl Acad Sci U S A. 2002, 99: 8695-8700.PubMed CentralView ArticlePubMedGoogle Scholar
- Santos-Rosa H, Schneider R, Bannister AJ, Sherriff J, Bernstein BE, Emre NC, et al: Active genes are tri-methylated at K4 of histone H3. Nature. 2002, 419: 407-411. 10.1038/nature01080.View ArticlePubMedGoogle Scholar
- Mikkelsen TS, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, et al: Genome-wide maps of chromatin state in pluripotent and lineage-committed cells. Nature. 2007, 448: 553-560. 10.1038/nature06008.PubMed CentralView ArticlePubMedGoogle Scholar
- Schuettengruber B, Chourrout D, Vervoort M, Leblanc B, Cavalli G: Genome regulation by polycomb and trithorax proteins. Cell. 2007, 128: 735-745. 10.1016/j.cell.2007.02.009.View ArticlePubMedGoogle Scholar
- Bradford Y, Conlin T, Dunn N, Fashena D, Frazer K, Howe DG, et al: ZFIN: enhancements and updates to the Zebrafish Model Organism Database. Nucleic Acids Res. 2011, 39: D822-D829. 10.1093/nar/gkq1077.PubMed CentralView ArticlePubMedGoogle Scholar
- Akkers RC, van Heeringen SJ, Jacobi UG, Janssen-Megens EM, Francoijs KJ, Stunnenberg HG, et al: A hierarchy of H3K4me3 and H3K27me3 acquisition in spatial gene regulation in Xenopus embryos. Dev Cell. 2009, 17: 425-434. 10.1016/j.devcel.2009.08.005.PubMed CentralView ArticlePubMedGoogle Scholar
- Bogdanovic O, van Heeringen SJ, Veenstra GJ: The epigenome in early vertebrate development. Genesis. 2012, 50: 192-206. 10.1002/dvg.20831.PubMed CentralView ArticlePubMedGoogle Scholar
- Mandell KJ, Parkos CA: The JAM family of proteins. Adv Drug Deliv Rev. 2005, 57: 857-867. 10.1016/j.addr.2005.01.005.View ArticlePubMedGoogle Scholar
- Hamazaki Y, Itoh M, Sasaki H, Furuse M, Tsukita S: Multi-PDZ domain protein 1 (MUPP1) is concentrated at tight junctions through its possible interaction with claudin-1 and junctional adhesion molecule. J Biol Chem. 2002, 277: 455-461.View ArticlePubMedGoogle Scholar
- Gallardo VE, Liang J, Behra M, Elkahloun A, Villablanca EJ, Russo V, et al: Molecular dissection of the migrating posterior lateral line primordium during early development in zebrafish. BMC Dev Biol. 2010, 10: 120-10.1186/1471-213X-10-120.PubMed CentralView ArticlePubMedGoogle Scholar
- Barone V, Heisenberg CP: Cell adhesion in embryo morphogenesis. Curr Opin Cell Biol. 2012, 24: 148-153. 10.1016/j.ceb.2011.11.006.View ArticlePubMedGoogle Scholar
- Siddiqui M, Sheikh H, Tran C, Bruce AE: The tight junction component Claudin E is required for zebrafish epiboly. Dev Dyn. 2010, 239: 715-722. 10.1002/dvdy.22172.View ArticlePubMedGoogle Scholar
- Foissac S, Sammeth M: ASTALAVISTA: dynamic and flexible analysis of alternative splicing events in custom gene datasets. Nucleic Acids Res. 2007, 35: W297-W299. 10.1093/nar/gkm311.PubMed CentralView ArticlePubMedGoogle Scholar
- Jurkowska RZ, Jurkowski TP, Jeltsch A: Structure and function of mammalian DNA methyltransferases. Chembiochem. 2011, 12: 206-222. 10.1002/cbic.201000195.View ArticlePubMedGoogle Scholar
- Rountree MR, Bachman KE, Baylin SB: DNMT1 binds HDAC2 and a new co-repressor, DMAP1, to form a complex at replication foci. Nat Genet. 2000, 25: 269-277. 10.1038/77023.View ArticlePubMedGoogle Scholar
- Lee GE, Kim JH, Taylor M, Muller MT: DNA methyltransferase 1-associated protein (DMAP1) is a co-repressor that stimulates DNA methylation globally and locally at sites of double strand break repair. J Biol Chem. 2010, 285: 37630-37640. 10.1074/jbc.M110.148536.PubMed CentralView ArticlePubMedGoogle Scholar
- Stancheva I, Meehan RR: Transient depletion of xDnmt1 leads to premature gene activation in Xenopus embryos. Genes Dev. 2000, 14: 313-327.PubMed CentralPubMedGoogle Scholar
- Dunican DS, Ruzov A, Hackett JA, Meehan RR: xDnmt1 regulates transcriptional silencing in pre-MBT Xenopus embryos independently of its catalytic function. Development. 2008, 135: 1295-1302. 10.1242/dev.016402.View ArticlePubMedGoogle Scholar
- Schier AF: The maternal-zygotic transition: death and birth of RNAs. Science. 2007, 316: 406-407. 10.1126/science.1140693.View ArticlePubMedGoogle Scholar
- Bogdanovic O, Long SW, van Heeringen SJ, Brinkman AB, Gomez-Skarmeta JL, Stunnenberg HG, et al: Temporal uncoupling of the DNA methylome and transcriptional repression during embryogenesis. Genome Res. 2011, 21: 1313-1327. 10.1101/gr.114843.110.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang X, Dai J: Concise review: isoforms of OCT4 contribute to the confusing diversity in stem cell biology. Stem Cells. 2010, 28: 885-893.PubMed CentralView ArticlePubMedGoogle Scholar
- Vargas DY, Shah K, Batish M, Levandoski M, Sinha S, Marras SA, et al: Single-molecule imaging of transcriptionally coupled and uncoupled splicing. Cell. 2011, 147: 1054-1065. 10.1016/j.cell.2011.10.024.PubMed CentralView ArticlePubMedGoogle Scholar
- Wagner EJ, Carpenter PB: Understanding the language of Lys36 methylation at histone H3. Nat Rev Mol Cell Biol. 2012, 13: 115-126. 10.1038/nrm3274.PubMed CentralView ArticlePubMedGoogle Scholar
- Luco RF, Allo M, Schor IE, Kornblihtt AR, Misteli T: Epigenetics in alternative pre-mRNA splicing. Cell. 2011, 144: 16-26. 10.1016/j.cell.2010.11.056.PubMed CentralView ArticlePubMedGoogle Scholar
- Kuehner JN, Pearson EL, Moore C: Unravelling the means to an end: RNA polymerase II transcription termination. Nat Rev Mol Cell Biol. 2011, 12: 283-294.View ArticlePubMedGoogle Scholar
- Millevoi S, Vagner S: Molecular mechanisms of eukaryotic pre-mRNA 3' end processing regulation. Nucleic Acids Res. 2010, 38: 2757-2774. 10.1093/nar/gkp1176.PubMed CentralView ArticlePubMedGoogle Scholar
- Proudfoot NJ: Ending the message: poly(A) signals then and now. Genes Dev. 2011, 25: 1770-1782. 10.1101/gad.17268411.PubMed CentralView ArticlePubMedGoogle Scholar
- Ji Z, Luo W, Li W, Hoque M, Pan Z, Zhao Y, et al: Transcriptional activity regulates alternative cleavage and polyadenylation. Mol Syst Biol. 2011, 7: 534-PubMed CentralView ArticlePubMedGoogle Scholar
- Hilgers V, Perry MW, Hendrix D, Stark A, Levine M, Haley B: Neural-specific elongation of 3' UTRs during Drosophila development. Proc Natl Acad Sci U S A. 2011, 108: 15864-15869. 10.1073/pnas.1112672108.PubMed CentralView ArticlePubMedGoogle Scholar
- Ji Z, Lee JY, Pan Z, Jiang B, Tian B: Progressive lengthening of 3' untranslated regions of mRNAs by alternative polyadenylation during mouse embryonic development. Proc Natl Acad Sci U S A. 2009, 106: 7028-7033. 10.1073/pnas.0900028106.PubMed CentralView ArticlePubMedGoogle Scholar
- Ulitsky I, Shkumatava A, Jan C, Subtelny AO, Koppstein D, Bell G, et al: Extensive alternative polyadenylation during zebrafish development. Genome Res. 2012Google Scholar
- Turi A, Loglisci C, Salvemini E, Grillo G, Malerba D, D'Elia D: Computational annotation of UTR cis-regulatory modules through Frequent Pattern Mining. BMC Bioinformatics. 2009, 10 (Suppl 6): S25-10.1186/1471-2105-10-S6-S25.PubMed CentralView ArticlePubMedGoogle Scholar
- Lai EC, Burks C, Posakony JW: The K box, a conserved 3' UTR sequence motif, negatively regulates accumulation of enhancer of split complex transcripts. Development. 1998, 125: 4077-4088.PubMedGoogle Scholar
- Stetefeld J, Ruegg MA: Structural and functional diversity generated by alternative mRNA splicing. Trends Biochem Sci. 2005, 30: 515-521. 10.1016/j.tibs.2005.07.001.View ArticlePubMedGoogle Scholar
- Dobrosotskaya IY, James GL: MAGI-1 interacts with beta-catenin and is associated with cell-cell adhesion structures. Biochem Biophys Res Commun. 2000, 270: 903-909. 10.1006/bbrc.2000.2471.View ArticlePubMedGoogle Scholar
- Zhang M, Zhang J, Lin SC, Meng A: beta-Catenin 1 and beta-catenin 2 play similar and distinct roles in left-right asymmetric development of zebrafish embryos. Development. 2012, 139: 2009-2019. 10.1242/dev.074435.View ArticlePubMedGoogle Scholar
- Ranganathan R, Ross EM: PDZ domain proteins: scaffolds for signaling complexes. Curr Biol. 1997, 7: R770-R773. 10.1016/S0960-9822(06)00401-5.View ArticlePubMedGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al: The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009, 25: 2078-2079. 10.1093/bioinformatics/btp352.PubMed CentralView ArticlePubMedGoogle Scholar
- Quinlan AR, Hall IM: BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010, 26: 841-842. 10.1093/bioinformatics/btq033.PubMed CentralView ArticlePubMedGoogle Scholar
- R Core Team. R: A language and environment for statistical computing. 2012, Vienna, Austria: R foundation for statistical computing, ISBN 3-900051-07-0, URLhttp://www.R-project.org/Google Scholar
- Robinson JT, Thorvaldsdottir H, Winckler W, Guttman M, Lander ES, Getz G, et al: Integrative genomics viewer. Nat Biotechnol. 2011, 29: 24-26. 10.1038/nbt.1754.PubMed CentralView ArticlePubMedGoogle Scholar
- Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al: Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012, 7: 562-578.PubMed CentralView ArticlePubMedGoogle Scholar
- Rice P, Longden I, Bleasby A: EMBOSS: the European Molecular Biology Open Software Suite. Trends Genet. 2000, 16: 276-277. 10.1016/S0168-9525(00)02024-2.View ArticlePubMedGoogle Scholar
- Grillo G, Turi A, Licciulli F, Mignone F, Liuni S, Banfi S, et al: UTRdb and UTRsite (RELEASE 2010): a collection of sequences and regulatory motifs of the untranslated regions of eukaryotic mRNAs. Nucleic Acids Res. 2010, 38: D75-D80. 10.1093/nar/gkp902.PubMed CentralView ArticlePubMedGoogle Scholar
- Lindeman LC, Vogt-Kielland LT, Alestrom P, Collas P: Fish'n ChIPs: chromatin immunoprecipitation in the zebrafish embryo. Methods Mol Biol. 2009, 567: 75-86. 10.1007/978-1-60327-414-2_5.View ArticlePubMedGoogle Scholar
- Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25: 1754-1760. 10.1093/bioinformatics/btp324.PubMed CentralView ArticlePubMedGoogle Scholar
- Xu H, Handoko L, Wei X, Ye C, Sheng J, Wei CL, et al: A signal-noise model for significance analysis of ChIP-seq with negative control. Bioinformatics. 2010, 26: 1199-1204. 10.1093/bioinformatics/btq128.View ArticlePubMedGoogle Scholar
- Statham AL, Strbenac D, Coolen MW, Stirzaker C, Clark SJ, Robinson MD: Repitools: an R package for the analysis of enrichment-based epigenomic data. Bioinformatics. 2010, 26: 1662-1663. 10.1093/bioinformatics/btq247.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhu LJ, Gazin C, Lawson ND, Pages H, Lin SM, Lapointe DS, et al: ChIPpeakAnno: a Bioconductor package to annotate ChIP-seq and ChIP-chip data. BMC Bioinformatics. 2010, 11: 237-10.1186/1471-2105-11-237.PubMed CentralView ArticlePubMedGoogle Scholar
- Peterson SM, Freeman JL: RNA isolation from embryonic zebrafish and cDNA synthesis for gene expression analysis. J Vis Exp. 2009,Google Scholar