44 OCT4 is a fundamental component of the molecular circuitry governing 45 pluripotency in vivo and in vitro. To determine how OCT4 protects the 46 pluripotent lineage from differentiation into trophoblast, we used single cell 47 transcriptomics and quantitative immunofluorescence on blastocysts and 48 established differentially expressed genes and pathways between control and 49 OCT4 null cells. Activation of most pluripotency-associated transcription 50 factors in the early mouse inner cell mass appears independent of OCT4, 51 whereas JAK/STAT signalling requires OCT4, via activation of IL6ST. Single 52 cell deconvolution, diffusion component and trajectory inference dissected the 53 process of differentiation of OCT4 null cells by activating specific gene-54 network and transcription factors. Downregulation of glycolytic and oxidative 55 metabolism was observed. CHIPseq analysis suggests OCT4 directly targets 56 rate-limiting glycolytic enzymes. Concomitant with significant disruption of the 57 STAT3 pathway, oxidative respiration is significantly diminished in OCT4 null 58 cells. Upregulation of the lysosomal pathway detected in OCT4 null embryos 59 is likely attributable to aberrant metabolism. 60 Results 104 Single cell transcriptional profiling reveals divergence of OCT4 null from wild 105 type and heterozygous ICM cells by the mid blastocyst stage 106 Following observation of NANOG protein in OCT4 null blastocysts24,25, we 107 performed whole transcriptome analysis by scRNAseq of ICM cells isolated 108 from Pou5f1 heterozygous inter se mating at embryonic day (E)3.5 to 109 investigate the entire pluripotency network. Quality control, as previously 110 reported27, eliminated inadequate samples, leaving 29 null mutant (MUT), 42 111wild-type (WT) and 16 heterozygous (HET) cells from 4, 5 and 2 mid 112 blastocysts, respectively (Fig.1A, Suppl.Table1). Pou5f1 RNA expression was 113 not detected in MUT embryos, confirming degradation of maternal transcripts 114 ( Fig.1A, Fig.S1A), consistent with lack of OCT4 protein observed at the 115 morula stage13. To characterize global differences and similarities between 116 genotypes t-SNE analysis was performed (Fig.1B, Fig.S1A) using the most 117 variable genes identified in early blastocysts (n=2232, log2FPKM>0.5, 118 logCV2>0.5). MUT cells cluster separately from HET and WT, suggesting 119 major changes in transcriptome, despite no apparent difference in ICM 120 morphology. Interestingly, HET and WT cells clustered together, indicating no 121 more than a negligible effect of reduced Pou5f1 in HET cells in the developing 122 embryo, contrasting with the elevated and more homogeneous expression of 123 Nanog, Klf4 and Estrogen-related-receptor B (Esrrb) previously reported in 124Oct4 HET ESCs28 (Fig.S1B). 125Weighted gene correlation network analysis (WGCNA) allows 126 extraction of modules defined by co-regulated genes, combined with 127 unsupervised clustering (Fig.1C, Fig.S1C). Two main modules emerged, 128 7 clustering cells according to genotype: module 1 co-clust...