Single blastomere expression profiling of Xenopus laevis embryos of 8 to 32-cells reveals developmental asymmetry.

We have measured the expression of 41 maternal mRNAs in individual blastomeres collected from the 8 to 32-cell Xenopus laevis embryos to determine when and how asymmetry in the body plan is introduced. We demonstrate that the asymmetry along the animal-vegetal axis in the oocyte is transferred to the daughter cells during early cell divisions. All studied mRNAs are distributed evenly among the set of animal as well as vegetal blastomeres. We find no asymmetry in mRNA levels that might be ascribed to the dorso-ventral specification or the left-right axis formation. We hypothesize that while the animal-vegetal asymmetry is a consequence of mRNA gradients, the dorso-ventral and left-right axes specifications are induced by asymmetric distribution of other biomolecules, probably proteins.

Understanding the formation of the embryonic body plan is essential in developmental biology. It is generally accepted that mechanisms of cellular differentiation driving development are based on differential transcription and translation of key regulatory elements. Recent development of highly sensitive and reliable molecular diagnostic methods allows for precise measurements of genes' transcriptional activities in very small samples including single cells1. In an excellent study by Guo et al., mouse embryos from 8 to 64 cells were dissociated into individual blastomeres and gene expression profiles were measured using the high-throughput microfluidic quantitative real-time PCR (qPCR) system BioMark2. They found early mouse embryos are composed of three types of cells representing the trophoectoderm, epiblast, and primitive endoderm. The three cell types were readily distinguished by the expression profile of 48 genes. Pluripotency of mouse single blastomeres from preimplantation embryos (1 to 16-cell stages) was also demonstrated with multiplex RT-PCR expression analysis by May et al.3. Expression profiles of trophoectoderm markers, inner cell mass markers and stemness markers have also been measured on human single blastomeres from 5 to 8-cell embryos by Galan et al., using microarray4 and in single putative stem cells from human adult ovarian surface epithelium5.

The Xenopus oocyte has two differentially colored hemispheres known as animal and vegetal. The separation into the hemispheres creates the first developmental axis of the embryo referred to as animal-vegetal. The darker color of the animal hemisphere is due to accumulation of pigmented granules called melanosomes. The vegetal hemisphere stores yolk and has light color6. De novo transcription in the developing Xenopus embryo is silenced until the mid-blastula transition stage (MBT), which takes place after twelve cell divisions post fertilization. All mRNA molecules required for the development into MBT must be present in the oocyte. These were transcribed during oogenesis from solely maternal chromosomes. The maternal mRNAs are asymmetrically distributed along the animal-vegetal axis of the oocyte and direct the specialization of the animal and vegetal parts789. The second embryonic body axis is formed after fertilization and separates the embryo into a dorsal and a ventral part. In Drosophila embryos the mechanism determining the dorsal-ventral axis is triggered by a concentration gradient of gurken maternal mRNA, which is accumulated at the future dorsal site. Absence of gurken mRNA at the ventral site leads to translocation of transcription factor Dorsal from the cytoplasm into the nuclei, where the Dorsal activates genes that specify the ventral part10. The mechanistic details of the dorso-ventral axis formation in zebrafish are not known, but differences between the dorsal and ventral parts appear during the gastrula stage, when the ventral side of the embryo gets thinner than the dorsal side11. The formation of the dorso-ventral axis in Xenopus laevis embryo can be discerned already at the 4-cell stage12. The sperm enters the oocyte through the animal hemisphere and in about 25 minutes the cortical cytoplasm of the oocyte rotates some 30 degrees relative to the inner cytoplasmic mass in opposite direction to the sperm's entrance1314. Dorsal maternal factors, such as the Dishevelled protein (Dvl), move to the future dorsal side in the vegetal hemisphere of the oocyte15. Dvl inactivates gsk3β protein and protects β-catenin from degradation. The accumulation of β-catenin in the dorsal site of the embryo promotes local Wnt signaling, which regulates transcription of zygotic xnr3 and siamois, leading to dorsalization of the embryo1617. The first cell division cuts the oocyte through the point of sperm entrance in the animal pole and the groove elongates to the vegetal pole. The second cell division is perpendicular to the first and separates the embryo into the future dorsal and ventral halves. All four blastomeres contain cytoplasm from both the animal and vegetal hemispheres, but vary in size and pigmentation12. The third and last important body axis of Xenopus laevis embryos creates the right-left asymmetry, which manifests at the gastrula stage when the production of embryonic mRNAs has been initiated. Maternal vg1 protein is more active in the left side of the embryo and activates the expression of xnr1 (Xenopus nodal-related 1). Xnr1 induces pitx2 production, which is required for the left side formation. Injection of vg1 protein to the future right side of the embryo causes left-right deformation manifested by random localization of the heart and of the digestive tube18.

We have previously used quantitative real-time PCR (qPCR) expression profiling of cryostat sections (qPCR tomography) of a single oocyte to identify transcripts that form gradients along the animal-vegetal axis89. These transcripts will become differentially distributed among the animal and the vegetal blastomeres formed at the third cell division. The purpose of the present study is to elucidate if there are expression gradients throughout the embryo that also reflect the dorsal-ventral axis specification and the left-right axis formation by measuring expression of selected 41 maternal genes in individual blastomeres collected from Xenopus embryos between 8 to 32-cell stages.

Results

First we extensively optimized our experimental procedure for high extraction and reverse transcription yields, and high qPCR efficiencies. We standardized the protocol for minimum technical variation for the expression profiling of single blastomeres. The vegetal blastomeres store yolk, which is strong inhibitor of biochemical reactions and can induce high variation in the reverse transcription yield and compromise qPCR efficiencies. An RNA spike (in vitro transcribed artificial RNA with 3′ polyA tail and 5′ cap, TATAA Biocenter) was used to validate the reverse transcription reaction and the qPCR. The spike was added to the RNA extracted from the samples together with the random sequence hexamers and oligo-dT primers for reverse transcription. The amount of reverse transcribed spike was then measured along with the endogenous transcripts. For the optimized protocol the standard deviation of the RNA spike Cq values across all samples was below 0.20 cycles, which evidences excellent reproducibility between sample preparations. Embryos from the same mother also showed high concordance. Our protocol for expression profiling is based on collecting the individual blastomeres in random order and we can only keep track from which hemisphere (animal or vegetal) the different blastomeres originated. This was necessary to avoid introducing bias during collection. With this protocol any dorso-ventral (or possibly left-right) asymmetry would be reflected by heterogeneity in expression profiles among the blastomeres collected from the same hemisphere. As markers to probe the formation of the three developmental axes we chose ten maternal transcripts that have been implicated in the dorso-ventral patterning: dvl2, dvl3, lrp6, wnt11, tcf3, gsk3b, ctnnb1( = β-catenin), foxh1, trim36 and axin116192021, one important marker for left-right specification: vg118, and nine genes implicated in animal-vegetal orientation: dazl, cdx1 ( = xcad2), wnt11, vg1, vegt, trim36, ddx25 ( = deadsouth), otx1, and maml19192223. We also included 25 mRNAs that have previously been observed in the mature oocyte: fzd7, bmp2, pias1, foxr1, frat1, mapk8, odc1, 18S rRNA, 5S rRNA, cyc1, acta, tubb, gapdh, eef1a1, RNA polymerase II, U3 snoRNA, par1, oct60, est1, apc, zpc, mos, stat3, and pax689.

Expression data were collected for a total of 224 single blastomeres. The data were analyzed with principal component analysis (PCA) and hierarchical clustering, which are the most powerful multivariate methods to classify samples based on the collective expression of multiple genes24. PCA clearly separated cells originating from the animal and from the vegetal hemispheres into clusters (Fig. 1), and indicated that the vegetal blastomeres are more heterogeneous than the animal blastomeres. Neither of the clusters, however, divided into subgroups that would indicate additional asymmetry that could be ascribed to the formation of the other body axes (e.g. dorsal-ventral and left-right asymmetry). The PCA results were confirmed by classifications using SOM (not shown) and hierarchical clustering. The hierarchical clustering is presented in a heatmap (Fig. 2), which is a graphical illustration of all the genes' expressions in all the blastomeres with the level indicated by color, and dendrograms clustering the blastomeres and the genes, respectively, shown on the top and left side. Separate heatmaps were calculated for each embryo from 8-cell (six embryos), 16-cell (seven embryos), and 32-cell (two embryos) stages. Blastomeres of 8 and 16-cell stages separated into well-defined animal and vegetal clusters, but no subclusters that could be ascribed to the formation of dorsal-ventral or left-right body axes can be discerned. Blastomeres from the 32-cell stage also divided into an animal and a vegetal cluster. The 32-cell stage vegetal cluster was diffuse, reflecting heterogeneity among the vegetal blastomeres possibly caused by the blastomeres being arranged into two ventral layers and some transcripts, such as germ plasm determinants, being concentrated in the extreme ventral cortical layer. The genes cluster into four well separated groups at all the developmental stages studied. The first group contains ubiquitously expressed genes represented by 18S rRNA, 5S rRNA, and cyc1. These genes are expressed at high level in all the blastomeres at all stages. A second group is composed of vg1, cdx1, vegt, dazl, wnt11, otx1, trim36, and ddx25, which are most abundant in the vegetal blastomeres. Several of these transcripts have previously been found in the vegetal part of the mature oocyte9. Third and fourth groups of transcripts are more abundant in the animal blastomeres, but subtle difference in expression patterns separate them into two clusters. Fzd7, bmp2, pias1, dvl2, dvl3, lrp6, frat1, mapk8, axin1, est1, and U3 snoRNA are expressed at slightly lower level (indicated by the green shade in the heatmap in Fig. 2). Than the fourth set of genes; ctnnb1, foxh1, odc1, zpc, mos, maml1, eef1a1, RNA polymerase II, foxr1, tcf3, gsk3b, oct60, stat3, par1, acta, tubb, gapdh, apc, and pax6 (indicated by the green shade turning to light red in the heatmap in Fig. 2). The separation of genes into these two groups shows some variability across embryos (data not shown). The separation of the gene transcripts into clusters is even more evident in the principal component analysis (PCA) (Fig. 3). One cluster characterizing genes expressed in animal blastomeres: fzd7, bmp2, pias1, dvl2, dvl3, lrp6, frat1, mapk8, axin1, est1, U3 snoRNA, ctnnb1, foxh1, odc1, zpc, mos, maml1, eef1a1, RNA polymerase II, foxr1, tcf3, gsk3b, oct60, stat3, par1, acta, tubb, gapdh, apc, and pax6; a second cluster of genes expressed in vegatal blastomeres: vg1, cdx1, vegt, dazl, wnt11, otx1, trim36, and ddx25; a third cluster of high expressed genes: 18S rRNA, 5S rRNA, and cyc1.

Figure 1

PCA of profiled blastomeres collected from early developmental stages of Xenopus laevis.

Triangles indicate blastomeres from the 8-cell stage; squares from the 16-cell stage, and circles from the 32-cell stage. Blastomeres originating from the animal hemisphere are shown in red and those from the vegetal hemisphere are shown in blue. Left and right graphs are data from two different females (left: four embryos from 8-cell stage, four embryos from 16-cell stage, and one 32-cell stage embryo, in total 128 blastomeres; right: two embryos from 8-cell stage, three from 16-cell stage, and one from 32-cell stage, in total 96 blastomeres).

Figure 2

Hierarchical clustering of blastomeres and mRNAs of Xenopus laevis early embryos presented as heatmaps.

One embryo from each developmental stage was arbitrarily chosen for the cluster analysis. (A) 8-cell, (B) 16-cell and (C) 32-cell embryo. Green color indicates low expression and red high expression. The dendrograms clustering blastomeres and genes are shown in the top and left side of the heatmap, respectively. In the dendrograms similarity between blastomeres/genes is indicated by the height at which they are joined.

Figure 3

Principal component analysis of maternal genes.

(A) 8-cell stage (four embryos/32 cells analyzed), (B) 16-cell stage (four embryos/64 cells analyzed), and (C) 32-cell stage (32 cells analyzed) embryos. Three clusters are seen in all the developmental stages. First cluster comprises 18S rRNA, 5S rRNA, and cyc1(mitochondrial cytochrome c); second cluster vg1, cdx1 (xcad2), wnt11, dazl, vegt, ddx25 (deadsouth), otx1, and trim36; third cluster (red squares) fzd7, par1, bmp2, pias1, dvl2, dvl3, lrp6, foxr1, fart1, mapk8, odc1, axin1, est1, U3 snoRNA, apc, gapdh, acta, eef1a1, tcf3, zpc, RNA polymerase II, gsk3b, maml1, tubb, ctnnb1 (β-catenin), mos, oct60, foxh1, stat3, and pax6. The animal-vegetal distinction is found along the y-axis (PC2), while differences in expression level of the maternal transcripts is reflected along the x-axis (PC1).

Discussion

Maternal mRNAs are synthesized during oogenesis and translated during the early development of Xenopus laevis. It is generally accepted that the temporal regulation and the spatial distribution of these maternal mRNAs are important for the establishment of the body axes and subsequently the whole body pattern. The animal-vegetal asymmetry is the first body axis formed during oogenesis. The germ layers, endoderm, mesoderm, and ectoderm, are formed along the animal-vegetal axis such that the vegetal part of the oocyte becomes the endoderm and the most distal animal part of the oocyte gives rise to the ectoderm2526. Transcripts of genes coding specific endodermal and mesodermal factors, such as vg1, wnt11, and vegt, are localized in the vegetal hemisphere8922. Genes responsible for primordial germ cell formation, such as dazl and cdx1, are localized towards the extreme vegetal pole9. Previously we described the new method qPCR tomography to quantify maternal mRNA distribution along the animal-vegetal axis formed during oogenesis, and found two dominant profiles: one with animal and one with vegetal mRNA localization8. Subsequent high-resolution qPCR tomography allowed us to discern two vegetal profiles; one encompassing germ plasm determinants (dazl, cdx1, ddx25) and one encompassing the other vegetal genes investigated (vg1, vegt, wnt11, otx1, eg6)9. In this study we show that the animal-vegetal polarizations observed in the oocyte remain in the 8, 16, and 32-cell stage embryos. The vegetal blastomeres had higher levels of dazl, cdx1, ddx25, vg1, vegt, wnt11, otx1, as well as trim36, which was not included in the oocyte study. At the 8 and 16-cell stages the vegetal mRNAs were distributed evenly among the vegetal blastomeres; the determinants of the germ plasm being present in all the vegetal cells. At the 32-cell stage, when the vegetal blastomeres form two layers, heterogeneity among the vegetal blastomers was introduced. 18S rRNA, 5S rRNA, cyc1, acta, tubb, gapdh, eef1a1, RNA polymerase II, U3 snoRNA, oct60, axin1, est1, apc, tcf3, zpc, gsk3b, maml1, ctnnb1, mos, foxh1, stat3, fzd7, par7, bmp2, pias1, dvl2, dvl3, lrp6, foxr1, frat1, mapk8, odc1 and pax6 transcripts were more abundant in the animal blastomeres, and were evenly distributed among the individual cells. Multivariate statistical analysis revealed higher variability among the vegetal blastomeres than among the animal. The difference was particularly pronounced at the 32-cell stage. The heterogeneity seems random, possibly caused by yolk platelets that may interfere with PCR inducing technical noise or introduced by the formation of two ventral layers; but we find no evidence for systematic distribution of transcripts neither along the dorsal-ventral axis nor the left-right axis.

An essential question in developmental biology is how amphibians establish the dorsal-ventral axis. The site of sperm entry determines the future ventral half of the embryo. The dorsal part forms on the opposite side and is determined by the accumulation of maternal dorsal determinants during the cortical rotation272829. Consequently, blastomeres of the early Xenopus embryos can be separated into a dorsal and a ventral group with equal number of animal and vegetal cells in each. In Drosophila the dorsalizing determinants have been shown to be maternal mRNA molecules that are asymmetrically distributed within the oocyte and regulate translation that leads to the polarization of the embryo3031. In our study we tested if the dorsalizing factors in Xenopus embryogenesis are also mRNAs coding for dorsal specific proteins within the Wnt pathway. We found that among the maternal transcripts that have been implicated in the induction of dorso-ventral asymmetry ctnnb1, dvl2, dvl3, lrp6, gsk3b, tcf3, foxh1, and axin1, are present predominantly in the animal blastomeres, while trim36 and xwnt11 are more abundant in the vegetal blastomeres. This supports the idea that the majority of the dorsal-ventral inducing molecules is localized in the dorsal animal blastomeres of the 16-cell stage embryo32, although we found no evidence of asymmetric distribution of maternal mRNAs among the dorsal or ventral blastomeres. In contrast to Tao et al.33, we found wnt11 distributed evenly between the dorsal and ventral cells. Our results rather support the finding of Schroeder et al.34, who found maternal wnt11 mRNA symmetrically distributed along the dorsal-ventral axis, and proposed that spatially regulated translation leads to asymmetric distribution of the Wnt11 protein. None of the 41 maternal mRNAs included in this study shows dorso-ventral polarization, although dorso-ventral asymmetry is observed starting from the 4-cell stage12. Of course, the dorso-ventral asymmetry could be induced by entirely different genes, but since our panel includes ten genes (dvl2, dvl3, lrp6, wnt11, tcf3, gsk3b, ctnnb1, foxh1, trim36, and axin1) previously implicated in dorsal-ventral patterning16192021, this would be surprising. Rather, we suggest the dorsal-ventral asymmetry in Xenopus laevis is due to early accumulation of dorsal proteins that regulate signaling pathways, which after the midblastula transition activate the expression of downstream dorsalizing factors at the dorsal site of the embryo. We speculate these proteins are present already in the oocyte and accumulate in the future dorsal site during the cortical rotation. Our set includes only one, gene (vg1) previously shown to be involved in the left-right specification18. We find no evidence of asymmetric distribution of vg1 or of any other of the studied transcripts that could be implicated in left-right specification.

Methods

Ethics statement

This study was carried out in accordance with the Act No 246/1992 Coll., on the protection of animals against cruelty. Official permission was issued to Faculty of Science, Charles University in Prague by the Central Commission for Animal Welfare under the Ministry of Agriculture of Czech Republic (accreditation No.24773/2008-10001, date of expiry 10.12.2013). The experiment was also approved by the institutional ethics committee.

In vitro fertilization and collection of single blastomeres

Two Xenopus laevis females were injected with 450 U of human chorionic gonadotrophin hormone (hCG) and kept overnight at room temperature. Ovulated oocytes were obtained 12 hours after stimulation by manual squeezing. Testes of Xenopus laevis males were removed and homogenized in L-15 Leibowitźs medium with 15% of fetal bovine serum. Oocytes were in vitro fertilized by sperm suspension and 1× MMR medium (Marćs Modified Ringers medium: 0.1 M NaCl, 2.0 mM KCl, 1 mM MgCl2, 2 mM CaCl2, 5 mM HEPES, pH 7.7) was added to cover the fertilized oocytes. Jelly coats were removed by cysteine treatment (2.2% in MMR, for 8 minutes) at 20 minutes after fertilization. During cell division the embryos were transferred to Ca2+ and Mg2+ free medium (88 mM NaCl, 1 mM KCl, 2.4 mM NaHCO3, 7.5 mM Tris, pH 7.7) and kept there for 60 minutes in order to separate single blastomeres. Vitelline membranes of the embryos containing 8, 16, and 32 cells were manually torn and individual blastomeres were collected. The blastomeres were collected separately from animal and vegetal hemispheres, but in random order to avoid introducing any technical bias and subjective dorso-ventral specification before measurement.

Isolation of total RNA and reverse transcription

The individual blastomeres collected from the 8 and 16-cell stages were homogenized in 200 μl and blastomeres from the 32-cell stage in 100 μl of TRIzol reagent (Invitrogen), and immediately frozen at −80°C. After thawing, the samples were carefully vortexed for 2 minutes to pulp yolk. Manufacturer's instructions were followed during isolation of total RNA, and glycogen was added to enhance RNA precipitation yield. RNA from individual cells from the 8 and 16-cell stages was dissolved in 20 μl of RNase-DNase free water. Samples from the 32-cell stage were dissolved in 12 μl of water. RNA concentration was measured using the Nanodrop ND1000 quantification system. cDNA was produced with SuperScript™ III Reverse transcriptase kit (Invitrogen). 10 ng of total RNA, 0.5 μl of a mixture of oligo-dT and random hexamers (mixture 1:1, 50 μM each), 0.5 μl of dNTPs (10 mM each), and 0.5 μl of RNA spike (in vitro transcribed artificial RNA with 3′ polyA tail and 5′ cap, TATAA Biocenter) were mixed with sterile water to a total volume of 6.5 μl. Samples were incubated for 5 min at 70°C, followed by 20 sec at 25°C, and then cooled at 4°C for 1 min. 100 U of SuperScript III reverse transcriptase, 20 U of RNaseOUT™ (recombinant ribonuclease inhibitor, Invitrogen), 0.5 μl of 0.1 M DTT and 2 μl of 5× first strand synthesis buffer was added to a final volume of 10 μl. The following temperature gradient profile was used to synthesize cDNA: 5 min at 25°C, 60 min at 50°C, 15 min at 55°C, and 15 min at 75°C. cDNA samples were diluted 8 times to a final volume of 80 μl and were stored at −20°C.

Primer design and quantitative PCR

Primer assays for the selected 41 maternal genes (fzd7, otx1, bmp2, pias1, dvl2, dvl3, lrp6, trim36, foxr1, frat1, mapk8, odc1, 18S rRNA, 5S rRNA, cyc1, acta, tubb, gapdh, maml1, eef1a1, RNA polymerase II, U3 snoRNA, par1, oct60, ddx25, dazl, cdx1, wnt11, vg1, vegt, axin1, est1, apc, tcf3, zpc, gsk3b, ctnnb1, mos, foxh1, stat3, pax6) were designed with Primer3Plus (http://www.bioinformatics.nl/cgi-bin/primer3plus/primer3plus.cgi). Annealing temperature was set to 60°C and the length of qPCR products was around 120bp. Specificity of the assays was tested in silico using blast (http://blast.ncbi.nlm.nih.gov/Blast.cgi). Primer sequences are available upon request. qPCR and melting curve validations were measured on the real-time CFX384 cycler system (BioRad). qPCR mix contained 2 μl of cDNA, 0.5 μl of forward and reverse primers (mixture 1:1, 10 μM each), 5 μl of SYBR Green JumpStart™ Taq Ready Mix™ (Sigma), and water in a total volume of 10 μl. The cycling program was: 2 min at 95°C for activation of the polymerase, followed by 40 cycles of denaturation at 95°C for 15 sec, annealing at 60°C for 20 sec and elongation at 72°C for 30 sec. Post PCR melting curves were measured from 65 to 95°C in 0.5°C intervals to validate the formation of expected PCR products.

Data analysis

Data were analyzed with GenEx (MultiD, version 5.3.). Off scale data were removed during pre-processing using a cut off at 36 cycles and outliers were identified with Grubb's test. Missing data were imputed using the average Cq of the other blastomeres from the same hemisphere in the same embryo (roughly 6% of the data was imputed). All data were normalized to the spike and converted to relative quantities (relative to the highest Cq for each gene, corresponding to arbitrary assigning an expression of 1 to the least expressed sample). Last step of the pre-processing was to transform the data to log2 scale. 8-cell, 16-cell and 32-cell embryos were analyzed separately. The combined expression of the individual blastomeres adds up to the expression of the original oocyte, which is expected, since de novo expression does not start until MBT. Three multivariate statistical methods were used for classification: principal component analysis (PCA), Kohonen self-organizing map (SOM), and hierarchical clustering presented as heatmaps. PCA of individual blastomeres (Fig. 1) was performed on data autoscaled along genes, while heatmaps (Fig. 2) and PCA of genes (Fig. 3) were performed on data autoscaled along samples35.

Author Contributions

M.F. performed the experiments and together with M.K. analyzed the obtained data. R.S. and M.K. conceived and designed the experiments. All authors co-wrote the manuscript and approved the final paper.