Zebrafish mRNA sequencing deciphers novelties in transcriptome dynamics during maternal to zygotic transition - PubMed (original) (raw)
. 2011 Aug;21(8):1328-38.
doi: 10.1101/gr.116012.110. Epub 2011 May 9.
Cecilia L Winata, Chi Ho Lin, Jieqi P Chen, Kandhadayar G Srinivasan, Serene G P Lee, Adrian Y M Lim, Hajira Shreen Hajan, Philippe Collas, Guillaume Bourque, Zhiyuan Gong, Vladimir Korzh, Peter Aleström, Sinnakaruppan Mathavan
Affiliations
- PMID: 21555364
- PMCID: PMC3149499
- DOI: 10.1101/gr.116012.110
Zebrafish mRNA sequencing deciphers novelties in transcriptome dynamics during maternal to zygotic transition
Håvard Aanes et al. Genome Res. 2011 Aug.
Abstract
Maternally deposited mRNAs direct early development before the initiation of zygotic transcription during mid-blastula transition (MBT). To study mechanisms regulating this developmental event in zebrafish, we applied mRNA deep sequencing technology and generated comprehensive information and valuable resources on transcriptome dynamics during early embryonic (egg to early gastrulation) stages. Genome-wide transcriptome analysis documented at least 8000 maternal genes and identified the earliest cohort of zygotic transcripts. We determined expression levels of maternal and zygotic transcripts with the highest resolution possible using mRNA-seq and clustered them based on their expression pattern. We unravel delayed polyadenylation in a large cohort of maternal transcripts prior to the MBT for the first time in zebrafish. Blocking polyadenylation of these transcripts confirms their role in regulating development from the MBT onward. Our study also identified a large number of novel transcribed regions in annotated and unannotated regions of the genome, which will facilitate reannotation of the zebrafish genome. We also identified splice variants with an estimated frequency of 50%-60%. Taken together, our data constitute a useful genomic information and valuable transcriptome resource for gene discovery and for understanding the mechanisms of early embryogenesis in zebrafish.
Figures
Figure 1.
Overview of mRNA-seq data and mapping to the zebrafish genome. (A) Zebrafish embryonic developmental stages analyzed in our study. Libraries were generated from eggs, 1-cell, 16-cell, 128-cell, 3.5-hpf, and 5.3-hpf embryos. The graph depicts general expression patterns during these developmental stages. (B) Total number of mRNA-seq reads mapped in each developmental-stage library. Differences in the number of generated reads were observed, but the proportions of mapped reads were similar between all stages. (Blue bars) The number of reads generated from sequencing; (red bars) the number of reads mapped to the genome; (green bars) the number of reads mapped to Ensembl-annotated genes. (C–F) Mapping of mRNA-seq reads as viewed in the UCSC Genome Browser. Horizontal red bars at the bottom of each panel represent annotated exons (tallest boxes), UTRs (half-sized boxes), and introns (lines with arrowheads). (Arrowheads on bars) Orientation of transcript from 5′ to 3′. Mapped reads are represented by red histograms. (C) Strand-specific mapping in a region of chromosome 5, showing reads mapping to mobkl1aa in the negative, and grsf1 in the positive strands. (D) Reads mapping to the htt gene, containing 67 exons. Note the sharp exon–intron boundary of the mapping. (E) Mapping of exon–exon junction spanning reads (red horizontal lines) to klf4 locus, allowing the derivation of splice junctions. (F) Mapping of reads to regions beyond annotated 3′ (zic3) or 5′ UTRs (foxa3), indicated by arrows. Di-tags (blue horizontal bars at top of panel) provide support for this mapping.
Figure 2.
Expression clusters of early developmental genes. (A) Heatmap showing distinct expression profiles of different clusters at each developmental stage. Values were scaled for each cluster; color intensity represents expression level relative to its own cluster's average. (B–D) Graphs showing mean expression of all genes in each cluster from the egg to the 5.3-hpf stage. (B) Two clusters of the maternal supercluster, showing two different degradation patterns. (C) Two clusters of the pre-MBT supercluster, both characterized by a significant increase at the 16-cell stage, with either decrease (pre-MBT1) or a continued increase after 3.5 hpf (pre-MBT2). (D) Zygotic supercluster consisting of genes activated at 3.5 hpf (MBT) and 5.3 hpf (post-MBT), or maternal–zygotic cluster genes that showed stable accumulation of transcripts during pre-MBT and with a sudden increase at MBT and post-MBT.
Figure 3.
Validation of mRNA-seq expression clusters. (A–C) RT-qPCR of the maternal gene cldng and the zygotic genes id1 and hspb1 show similar expression patterns to those detected by mRNA-seq. WISH confirmed the expression pattern of these genes as observed in mRNA-seq. (D) RT-qPCR of the pre-MBT gene xbp1 showed a different expression profile from that of mRNA-seq when random primers (r.p.) were used, while it showed good correlation when oligo d(T) primers were used. (E) WISH of xbp1 demonstrates the presence of maternal transcript in the egg and the 1-cell stage and confirmed r.p. RT-qPCR data. (F) Measurement of xbp1 poly(A) tail length by RT-PCR in the egg and 1-cell, 16-cell, and 5.3-hpf stages. Total RNA was G/I-tailed and amplified with the _xbp1_-specific forward primer (F) and a universal oligo d(C) reverse primer (U). A control RT-PCR with _xbp1_-specific F and R primers gave a 180-bp fragment (S). The increasing poly(A) tail appears as a smear that increases gradually from the 1-cell stage to the 5.3-hpf stage, indicating the presence of an extended poly(A) tail.
Figure 4.
Inhibition of cytoplasmic polyadenylation by 3′-dA. (A) Schematic of treatments at various developmental time frames. Horizontal arrows represent exposure to 3′-dA. (B–D) Control embryos showing normal development. (Red arrows) The level of marginal cells, signifying the progress of epiboly. Note the presence of embryonic shield (S) at the 6-hpf stage. (E–G) Embryos treated with 3′-dA from the 1-cell to 3.5-hpf stage, showing a severe delay in development. (H) U-rich motif found in genes of the pre-MBT supercluster using the MEME algorithm. (I) Distance between Hex and CPE in different groups of genes. The boxplot shows the distribution of the shortest distance between any Hex and CPE pair in the 3′ UTR of genes in maternal–zygotic and pre-MBT clusters, as well as in a randomly chosen group of genes (random). Notice the longer distance between the two elements in the pre-MBT group compared to the other two groups.
Figure 5.
Discovery of NTRs. (A) Genome-wide distribution of NTRs plotted at their mapped chromosomal positions. The color code depicts at which stage an NTR was detected. (B) Novel exons of ENSDARG00000035540, some of which correspond to exons of its human homolog, BRWD1. Splice junctions were also identified among these novel exons. RT-PCR performed using reverse primer at annotated exon 2 (R) and forward primers at each of the 12 NTRs (inset) showed that at least 11 of these NTRs are expressed and are part of the transcription unit. (C) Novel exon in foxa found 5′ of annotated transcription start site. Note the presence of the 5′ annotated exon at 3.5 hpf and 5.3 hpf (white arrowheads). Di-tag (blue horizontal lines) provides support for an isoform that includes the novel exon. Junction mapping at 5.3 hpf (green horizontal lines) indicated the presence of two distinct transcript isoforms. RT-PCR using a reverse primer (R) designed on the 3′ distal exon and forward primers on either 5′ annotated (*) or NTR (+) confirmed the existence of a new 5′ distal exon (inset). The 5′ annotated exon (*) is skipped in the maternal mRNAs (egg to 128-cell stage).
Figure 6.
Analysis of splice variants. (A) The number of alternatively spliced transcripts (blue bars) and alternatively spliced exons (red bars) at each stage. Both were stable through all stages. (B) Location of alternatively spliced exons in all transcripts. Most spliced exons were located between exons, but with a substantial fraction corresponding to distal 5′ exons. (C) Splice variants in msna at 5.3 hpf. Junction mapping shows skipping of exons 1 to 8, 2 to 8, and 3 to 11 (blue arrows). Di-tags (red bars with arrowheads) support the presence of several splice isoforms.
Similar articles
- Cytoplasmic polyadenylation-mediated translational control of maternal mRNAs directs maternal-to-zygotic transition.
Winata CL, Łapiński M, Pryszcz L, Vaz C, Bin Ismail MH, Nama S, Hajan HS, Lee SGP, Korzh V, Sampath P, Tanavde V, Mathavan S. Winata CL, et al. Development. 2018 Jan 8;145(1):dev159566. doi: 10.1242/dev.159566. Development. 2018. PMID: 29229769 - m6A-dependent maternal mRNA clearance facilitates zebrafish maternal-to-zygotic transition.
Zhao BS, Wang X, Beadell AV, Lu Z, Shi H, Kuuspalu A, Ho RK, He C. Zhao BS, et al. Nature. 2017 Feb 23;542(7642):475-478. doi: 10.1038/nature21355. Epub 2017 Feb 13. Nature. 2017. PMID: 28192787 Free PMC article. - Evolution of maternal and zygotic mRNA complements in the early Drosophila embryo.
Atallah J, Lott SE. Atallah J, et al. PLoS Genet. 2018 Dec 17;14(12):e1007838. doi: 10.1371/journal.pgen.1007838. eCollection 2018 Dec. PLoS Genet. 2018. PMID: 30557299 Free PMC article. - Zygotic genome activation during the maternal-to-zygotic transition.
Lee MT, Bonneau AR, Giraldez AJ. Lee MT, et al. Annu Rev Cell Dev Biol. 2014;30:581-613. doi: 10.1146/annurev-cellbio-100913-013027. Epub 2014 Aug 11. Annu Rev Cell Dev Biol. 2014. PMID: 25150012 Free PMC article. Review. - Transcriptional Regulation During Zygotic Genome Activation in Zebrafish and Other Anamniote Embryos.
Wragg J, Müller F. Wragg J, et al. Adv Genet. 2016;95:161-94. doi: 10.1016/bs.adgen.2016.05.001. Epub 2016 Jul 4. Adv Genet. 2016. PMID: 27503357 Review.
Cited by
- Quantitative proteomics reveals the dynamic proteome landscape of zebrafish embryos during the maternal-to-zygotic transition.
Fang F, Chen D, Basharat AR, Poulos W, Wang Q, Cibelli JB, Liu X, Sun L. Fang F, et al. iScience. 2024 May 8;27(6):109944. doi: 10.1016/j.isci.2024.109944. eCollection 2024 Jun 21. iScience. 2024. PMID: 38784018 Free PMC article. - Unveiling Gene Expression Dynamics during Early Embryogenesis in Cynoglossus semilaevis: A Transcriptomic Perspective.
Cheng X, Jiang W, Wang Q, Liu K, Dai W, Liu Y, Shao C, Li Q. Cheng X, et al. Life (Basel). 2024 Apr 15;14(4):505. doi: 10.3390/life14040505. Life (Basel). 2024. PMID: 38672775 Free PMC article. - Cell-type-specific mRNA transcription and degradation kinetics in zebrafish embryogenesis from metabolically labeled single-cell RNA-seq.
Fishman L, Modak A, Nechooshtan G, Razin T, Erhard F, Regev A, Farrell JA, Rabani M. Fishman L, et al. Nat Commun. 2024 Apr 10;15(1):3104. doi: 10.1038/s41467-024-47290-9. Nat Commun. 2024. PMID: 38600066 Free PMC article. - Interwoven processes in fish development: microbial community succession and immune maturation.
Auclert LZ, Chhanda MS, Derome N. Auclert LZ, et al. PeerJ. 2024 Mar 27;12:e17051. doi: 10.7717/peerj.17051. eCollection 2024. PeerJ. 2024. PMID: 38560465 Free PMC article. Review. - Mature mRNA processing that deletes 3' end sequences directs translational activation and embryonic development.
Takada Y, Fierro L, Sato K, Sanada T, Ishii A, Yamamoto T, Kotani T. Takada Y, et al. Sci Adv. 2023 Nov 24;9(47):eadg6532. doi: 10.1126/sciadv.adg6532. Epub 2023 Nov 24. Sci Adv. 2023. PMID: 38000026 Free PMC article.
References
- Alestrom P, Holter JL, Nourizadeh-Lillabadi R 2006. Zebrafish in functional genomics and aquatic biomedicine. Trends Biotechnol 24: 15–21 - PubMed
- Amores A, Force A, Yan YL, Joly L, Amemiya C, Fritz A, Ho RK, Langeland J, Prince V, Wang YL, et al. 1998. Zebrafish hox clusters and vertebrate genome evolution. Science 282: 1711–1714 - PubMed
- Amsterdam A, Hopkins N 2006. Mutagenesis strategies in zebrafish for identifying genes involved in development and disease. Trends Genet 22: 473–478 - PubMed
Publication types
MeSH terms
Substances
LinkOut - more resources
Full Text Sources
Molecular Biology Databases