Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2020 Dec 31;16(12):e1009294.
doi: 10.1371/journal.pgen.1009294. eCollection 2020 Dec.

Molecular asymmetry in the cephalochordate embryo revealed by single-blastomere transcriptome profiling

Affiliations

Molecular asymmetry in the cephalochordate embryo revealed by single-blastomere transcriptome profiling

Che-Yi Lin et al. PLoS Genet. .

Abstract

Studies in various animals have shown that asymmetrically localized maternal transcripts play important roles in axial patterning and cell fate specification in early embryos. However, comprehensive analyses of the maternal transcriptomes with spatial information are scarce and limited to a handful of model organisms. In cephalochordates (amphioxus), an early branching chordate group, maternal transcripts of germline determinants form a compact granule that is inherited by a single blastomere during cleavage stages. Further blastomere separation experiments suggest that other transcripts associated with the granule are likely responsible for organizing the posterior structure in amphioxus; however, the identities of these determinants remain unknown. In this study, we used high-throughput RNA sequencing of separated blastomeres to examine asymmetrically localized transcripts in two-cell and eight-cell stage embryos of the amphioxus Branchiostoma floridae. We identified 111 and 391 differentially enriched transcripts at the 2-cell stage and the 8-cell stage, respectively, and used in situ hybridization to validate the spatial distribution patterns for a subset of these transcripts. The identified transcripts could be categorized into two major groups: (1) vegetal tier/germ granule-enriched and (2) animal tier/anterior-enriched transcripts. Using zebrafish as a surrogate model system, we showed that overexpression of one animal tier/anterior-localized amphioxus transcript, zfp665, causes a dorsalization/anteriorization phenotype in zebrafish embryos by downregulating the expression of the ventral gene, eve1, suggesting a potential function of zfp665 in early axial patterning. Our results provide a global transcriptomic blueprint for early-stage amphioxus embryos. This dataset represents a rich platform to guide future characterization of molecular players in early amphioxus development and to elucidate conservation and divergence of developmental programs during chordate evolution.

PubMed Disclaimer

Conflict of interest statement

The authors have declared that no competing interests exist.

Figures

Fig 1
Fig 1. RNA-seq analyses of blastomeres from the 2-cell stage and animal/vegetal tiers from the 8-cell stage.
(A) Images of 2-cell and 8-cell amphioxus embryos, with an illustration of the strategy for RNA-seq. Dashed lines indicate the planes separating blastomeres of each paired sample. The presumptive germ granule-positive (p) and germ granule-negative (n), and animal-tier (a) vs. vegetal-tier (v) blastomeres are indicated. (B) A schematic representation of the experiment using single blastomeres from 2-cell stage embryos (six biological replicates). Blue dot denotes the germ granule in the presumptive germ granule-positive blastomere in each embryo (from 1p to 6p) based on qPCR analysis (S1 Fig); the other cell from the same embryo is designated as the germ granule-negative cell (from 1n to 6n). (C) Hierarchically clustered heatmap based on the Euclidean distances between 2-cell stage samples. Identification of the samples and their batch of origin are shown at the bottom of the heatmap. (D) PCA of the twelve 2-cell stage samples shows strong embryo and batch effects. (E) PCA of the twelve 2-cell stage samples after data normalization to account for individual embryo differences. (F) A schematic representation of the experiment on animal (1a-6a) and vegetal tiers (1v-6v) from 8-cell stage embryos (six biological replicates). Blue dots denote the germ granule that is present in one of the blastomeres of the vegetal tier. (G) Hierarchically clustered heatmap based on the Euclidean distances between 8-cell stage samples. Identification of the samples and batch of origin are shown at the bottom of the heatmap. Batch identities of 5a/v and 6a/v were not recorded. (H) PCA of the twelve 8-cell stage samples shows strong embryo and batch effects. (I) PCA of the twelve 8-cell stage samples after data normalization to account for individual embryo differences.
Fig 2
Fig 2. Identification of differentially enriched transcripts (DETs).
(A-B) MA-plots showing the log2 fold changes between the germ granule-positive and negative cells over the average transcript level of each gene among the 2-cell stage samples analyzed by using DESeq2 (A) or edgeR (B). The DETs (FDR < 0.1) are shown in red. (C) Distribution of FPKM values of DETs based on DESeq2 (orange) and edgeR (blue) analyses at the 2-cell stage. (D-E) MA-plots showing the log2 fold changes between the animal and vegetal tiers over the average transcript level of each gene among the 8-cell stage samples analyzed by using DESeq2 (D) or edgeR (E). (F) Distribution of FPKM values of DETs based on DESeq2 (orange) and edgeR (blue) analyses at the 8-cell stage. (G-J) Venn diagrams show the number of DETs in the germ granule-positive cell (G), germ granule-negative cell (H), animal tier (I) and vegetal tier (J), comparing results between DESeq2 and edgeR analyses. Underlying data are available in S1 Data.
Fig 3
Fig 3. Visualization of the transcriptome data showing asymmetric distributions of transcripts in the early amphioxus embryo.
(A) A scatter plot showing the log2 fold change of each gene between animal (positive y-axis) and vegetal tiers (negative y-axis) at the 8-cell stage, and its corresponding log2 fold change between germ granule-negative (positive x-axis) and positive blastomeres (negative x-axis) at the 2-cell stage (based on DESeq2). Size of the circle depicts FPKM values of the corresponding DET. Color of the circle indicates that the transcript is differentially enriched in a specific embryonic domain. Blue circles (total number = 264) are DETs in the animal tier; red circles (178) are DETs in the vegetal tier; yellow circles (58) are DETs in germ granule-positive blastomere; green circles (4) are DETs in germ granule-negative blastomere; orange circles (57) are DETs in both the germ granule-positive blastomere and the vegetal tier; purple circles (2) are DETs in both the germ granule-negative blastomere and the animal tier. Transcripts that are not enriched are in gray. (B) A scatter plot showing the spatial distributions of nine previously characterized transcripts. (C) Summary information for the nine previously characterized genes. GP, germ granule-positive blastomere. ID, gene model ID from genome assembly v2.0.
Fig 4
Fig 4. Validation of newly identified germ granule-enriched transcripts by WMISH.
(A) Visualization map of 13 germ granule-enriched transcripts based on DESeq2. (B) Annotations of the validated transcripts. (C) WMISH of the corresponding transcripts at 2- and 8-cell stages. The white arrows indicate aggregated signals on the vegetal side. Animal pole is up and vegetal pole is down. (D) Double fluorescence in situ hybridization showed that four newly identified transcripts were co-localized with nanos transcripts in the germ granule. The white arrowheads indicate co-localized signals. The white arrows indicate the compact germ granule signals. Embryo outlines are demarcated by dashed lines. The WMISH images show representative expression patterns (> 80%, n ≥ 5).
Fig 5
Fig 5. Validation of newly identified animal tier-enriched transcripts by WMISH.
(A) The scatter plot shows spatial distributions of ten animal tier-enriched transcripts based on DESeq2. (B) Annotations of the validated transcripts. (C) WMISH of the corresponding transcripts at the 2- and 8-cell stages. The arrows indicate the aggregated signals on the vegetal side. Animal is up and vegetal is down. (D) A schematic representation showing different views of a 4-cell stage amphioxus embryo, (i) animal pole view, (ii) lateral view, and (iii) posterior view. Transcripts enriched in the germ granule and toward the animal pole are illustrated in blue. (E) Double WMISH of zfp665 (198050) and tbx2 (124781) at the 4-cell stage. Embryos are observed from different views as indicated in (D). A surface focal plane of the posterior view (iii’) showing the germ granule-enriched transcripts. The WMISH images show representative expression patterns (> 80%, n ≥ 5).
Fig 6
Fig 6. Gene ontology (GO) enrichment analyses of the DETs.
(A) Results of the GO enrichment network analysis (adjusted p value < 0.1) of animal tier-enriched transcripts. Individual node of the network denotes a specific enriched GO term, including biological process (BP) and molecular function (MF). The GO terms are clustered and divided into different modules that are labeled with different colors. Unclassified GO terms are labeled as singletons. Sizes of the circles indicate numbers of genes in each GO term. (B) Descriptions of the most enriched GO terms within each module for the animal tier-enriched transcripts. The bars indicate -log10 adjusted p values for the corresponding GO terms. (C) Results of the GO enrichment network analysis (adjusted p value < 0.1) of germ granule-positive and/or vegetal tier-enriched transcripts. Individual node of network denotes a specific enriched GO term, including biological process (upper panel) and molecular function (lower panel). Lines connecting different nodes represent their semantic similarity. (D) Descriptions of GO terms within each module with the -log10 adjusted p value. Five GO terms with the smallest adjusted p values from different modules are presented. One manually selected GO term is indicated (*).
Fig 7
Fig 7. Functional validation of four asymmetrically localized transcripts encoding transcription factors.
(A) Phenotypes of the 48 hpf zebrafish embryos injected with amphioxus zfp665 mRNA. (B) Percentages of 48 hpf zebrafish embryos with the indicated phenotypes following injection of gfp (control) or amphioxus zfp665 mRNA. (C) WMISH of pax2.1, krox20, myoD and six3a in gfp (control) or amphioxus zfp665 mRNA-injected zebrafish embryos at the 10 somite stage. (D) Statistical analysis of the lengths of rhombomeres 5 in gfp and zfp665 mRNA injected zebrafish embryos at the 10 somite stage. (E) Phenotype of the 48 hpf zebrafish embryos injected with amphioxus soxB1b mRNA. (F) Percentages of 48 hpf zebrafish embryos with the indicated phenotypes after injection of gfp (control) or amphioxus soxB1b mRNA. (G) WMISH of pax2.1, krox20, myoD and six3a in gfp (control) or amphioxus soxB1b mRNA-injected zebrafish embryos at the 10 somite stage. (H) Statistical analysis of the lengths of rhombomeres 5 in gfp and soxB1b mRNA-injected zebrafish embryos at the 10 somite stage. (I) Embryos injected with amphioxus tbx2 or foxn1/4b mRNA did not show any detectable axial defects compared to control embryos at 72 hpf. DM, dorsal mesoderm; FB, forebrain; MHB, midbrain-hindbrain boundary; r3, rhombomere 3; r5, rhombomere 5. Statistical significance was determined by Student’s t-test. ***p < 0.001. Underlying data are available in S1 Data.
Fig 8
Fig 8. Expression domain of the ventral gene, eve1, was significantly reduced in zfp665 overexpressed zebrafish embryos.
(A-F) WMISH of zebrafish dorsal-ventral genes eve1 (A), vox (B), vent (C), bmp2b (D), gsc (E) and chd (F), in gfp (control) or amphioxus zfp665 mRNA-injected zebrafish embryos. The stages and views of the embryos are indicated below the images. Statistical analyses of eve1, vox, vent, gsc and chd expression domains in embryos injected with gfp and zfp665 mRNA. The expression domains were quantified by measuring the angles of the signals observed from the animal pole. Ectopic expression of vent (arrowhead) was observed in some zfp665-overexpressing embryos (44%, n = 27) (C). The ratio of embryos exhibiting the displayed bmp2b expression patterns is indicated in the bottom right-hand corner (D). V, ventral side; D, dorsal side. Statistical significance was determined by Student’s t-test, ***p < 0.001. Underlying data are available in S1 Data.
Fig 9
Fig 9. Functional validation of zfp665 as a transcriptional repressor in zebrafish embryos.
(A) A schematic representation of the protein domain organization of amphioxus zfp665, showing four consecutive C2H2-type zinc finger domains at the N-terminus. zfp665ΔC denotes the zinc finger region of zfp665 without its C-terminus. (B) A schematic representation of three different zfp665 fusion proteins, including engrailed repressor domain-zfp665ΔC fusion protein (En-zfp665ΔC), VP16 activation domain-zfp665ΔC (VP16-zfp665ΔC) fusion protein and the deletion of C-terminus of zfp665 (zfp665ΔC). (C) Phenotype of the zebrafish embryos injected with one of the three in vitro synthesized mRNA encoding the zfp665 fusion proteins at 24 hpf. (i) and (ii) are high-magnification views of the corresponding boxed areas. The black arrow indicates expanded ventral tissue. (D) Percentages of 24 hpf zebrafish embryos with the indicated phenotypes after injection of gfp (control) or one of the three in vitro synthesized amphioxus mRNAs encoding zfp665 fusion proteins. Only live larvae were scored. (E) WMISH of pax2.1, krox20, myoD and six3a in gfp (control) or indicated amphioxus mRNA-injected zebrafish embryos at the 10 somite stage. (F) Statistical analysis of the lengths of rhombomere 5 at the 10 somite stage in zebrafish embryos injected with gfp or one of the three forms of zfp665 mRNAs. DM, dorsal mesoderm; FB, forebrain; MHB, midbrain-hindbrain boundary; r3, rhombomere 3; r5, rhombomere 5. Statistical significance was determined by Welch’s one-way ANOVA with Games-Howell post hoc test plus Bonferroni correction when the number of groups was greater than 2 or Student’s t-test when the number of groups was 2. ***p < 0.001. Underlying data are available in S1 Data.

Similar articles

Cited by

References

    1. Stern CD. Evolution of the mechanisms that establish the embryonic axes. Current Opinion in Genetics & Development. 2006;16(4):413–8. 10.1016/j.gde.2006.06.005 WOS:000239614500013. - DOI - PubMed
    1. Lawrence PA, Levine M. Mosaic and regulative development: two faces of one coin. Curr Biol. 2006;16(7):R236–9. Epub 2006/04/04. 10.1016/j.cub.2006.03.016 . - DOI - PubMed
    1. Tao Q, Yokota C, Puck H, Kofron M, Birsoy B, Yan D, et al. Maternal wnt11 activates the canonical wnt signaling pathway required for axis formation in Xenopus embryos. Cell. 2005;120(6):857–71. Epub 2005/03/31. 10.1016/j.cell.2005.01.013 . - DOI - PubMed
    1. Langdon YG, Mullins MC. Maternal and zygotic control of zebrafish dorsoventral axial patterning. Annu Rev Genet. 2011;45:357–77. Epub 2011/09/29. 10.1146/annurev-genet-110410-132517 . - DOI - PubMed
    1. Lu FI, Thisse C, Thisse B. Identification and mechanism of regulation of the zebrafish dorsal determinant. Proc Natl Acad Sci U S A. 2011;108(38):15876–80. Epub 2011/09/14. 10.1073/pnas.1106801108 - DOI - PMC - PubMed

Publication types

Grants and funding

This work was supported by the Academia Sinica (https://www.sinica.edu.tw/en) intramural fund to JKY and YHS, and by the Ministry of Science and Technology, Taiwan, under the grants MOST-102-2311-B-001-011-MY3 and MOST-105-2628-B-001-003-MY3 to JKY, and MOST-107-2321-B-001-017 to YHS. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

LinkOut - more resources

-