Figure 4—figure supplement 1.

(A) Heatmap showing log2 fold activation differences for exonic and intronic regions of primary-activated genes for combinations of pou5f3.2, pou5f3.3, and sox3 morpholino-treatments, or Triptolide treatment, compared to controls. Right panel is in the presence of cycloheximide (CHX). (B) Biplot comparing exonic expression levels in cycloheximide-treated control embryos versus embryos also injected with pou5f3.2, pou5f3.3, and sox3 morpholinos. Primary-activated genes with maternal contribution <1 TPM (strictly zygotic) are purple circles, maternal-zygotic genes detected by exonic increases are orange triangles. TPM = transcripts per million. (C) Barplot summarizing the proportion of genes affected by morpholino treatment with cycloheximide on primary-activated genes (left bar), without cycloheximide (middle bar), and all stage 9 activated genes without cycloheximide (right bar). Down = significantly decreased in one of the morpholino treatments, up = significantly increased. (D, F) Biplots showing genes with >2 fold L or S biased activation (upper red and lower blue points, respectively) in control embryos (left panel) versus their activation in pou5f3.2, pou5f3.3, and sox3 morpholino-treated embryos (right panel, maintaining the same color per gene). (E, G) Quantification of the biplots in (D, F) in before-and-after plots. Y-axis is the absolute value of the log2 L vs S activation difference. p Values are from Wilcoxon signed-rank tests (paired). Overlaid boxplots show median, upper and lower quartiles, and 1.5 x interquartile range. (H) Regulatory networks consistent with direct regulation of embryonic gene activation by Pou5f3 and Sox3 (1) versus additional regulation by zygotic factors (2), which likely accounts for genes up-regulated in MO treatments. (I) Stage 8 Pou5f3.3 (left) and Sox3 (right) CUT&RUN coverage near TSSs for genes down-regulated in morpholino-treated embryos with or without cycloheximide (top), genes up-regulated (middle), and genes not significantly affected in any morpholino treatment (bottom). Top enriched motifs for each factor are shown below with p-values from Homer de novo discovery. (J) Aggregate plots of the binding signal in (I), with down-regulated genes further separated into genes down-regulated with morpholino treatment and cycloheximide (1°) or only down-regulated without cycloheximide (2°). p Values are from Kruskal-Wallis tests on summed signal per TSS. (K) Cumulative distributions of distance from a Pou5f3/Sox3-bound regulatory element for genes strongly (≥8 fold) and less strongly (<8 fold) down-regulated in morpholino-treated embryos with or without cycloheximide, compared to up-regulated, unaffected and unactivated genes. p Value is from a Kruskal-Wallis test. (L) Maps showing density of Pou5f3/Sox3-bound regulatory elements around paired homeologous TSSs, divided into elements with differential homeologous L & S binding (left panels) versus both bound (right panels). TSSs are grouped according to L versus S homeolog sensitivity to morpholino treatment. (M) Browser tracks showing CUT&RUN enrichment and ATAC-seq coverage near active homeolog hes3.L and inactive homeolog hes3.S. Seven L-specific high-confidence regulatory regions are highlighted with their homeologous S regions (bold ‘L’), as well as two lower-confidence enhancers, one of which also has weak activity in S, but minimal Pou5f3 or Sox3 binding (labeled ‘LS’).