Fig. 3

a Experimental setup: RNA-seq time series included eight time points, four replicates for the wild type, and two replicates for each MZspg, MZsox19b, and double mutant MZsox19bspg (see Fig. ^S4a for the details). Time series, where material was collected in the same experiment, have the same color. b Heatmap of all 4643 zygotic transcripts in the indicated genotypes. 2766 transcripts were downregulated in the mutants (groups A–G),1877 transcripts were not (group H). Zygotic transcripts were sorted by ascending switch time and switch p-value in the wild type. n number of transcripts. c Enrichment in Gene Ontology terms (DAVID). Top four categories per group and example genes are shown. GO: enrichment for “structural molecule activity” in group A was due to the battery of eight keratins, activated by Klf17 in the epithelial layer⁴⁶. d Non-additive (group A) and compensatory (group B) effects of Pou5f3 and Sox19b on the earliest zygotic transcription. Mean zygotic transcription profiles for the groups A and B, relative to 2.5 hpf. e Chordin/BMPs ratio is increased in the double mutant, but not in the single mutants throughout the time curve. Note that the transcription start of bmp2b and bmp7a, but not Chordin is delayed in the double mutant compared to the WT. In d, e: n(wt) = 4, n(MZspg) = 2, n(MZsox19b) = 2, n(double) = 2, where n is a number of biologically independent experiments. Data are presented as mean values ± SEM. Source data for b, d, and e are provided as a Source Data file. Source data for b–d are provided as Dataset ^S1.