Accurate specification of female and male germ cells during embryonic development is critical for sexual reproduction. Primordial germ cells (PGCs) are the bipotential precursors of mature gametes that commit to an oogenic or spermatogenic fate in response to sex-determining cues from the fetal gonad. The critical processes required for PGCs to integrate and respond to signals from the somatic environment in gonads are not well understood. In this study, we developed the first single-nucleus multiomics map of chromatin accessibility and gene expression during murine PGC development in both XX and XY embryos. Profiling of cell-type-specific transcriptomes and regions of open chromatin from the same cell captured the molecular signatures and gene networks underlying PGC sex determination. Joint RNA and ATAC data for single PGCs resolved previously unreported PGC subpopulations and cataloged a multimodal reference atlas of differentiating PGC clusters. We discovered that regulatory element accessibility precedes gene expression during PGC development, suggesting that changes in chromatin accessibility may prime PGC lineage commitment prior to differentiation. Similarly, we found that sexual dimorphism in chromatin accessibility and gene expression increased temporally in PGCs. Combining single-nucleus sequencing data, we computationally mapped the cohort of transcription factors that regulate the expression of sexually dimorphic genes in PGCs. For example, the gene regulatory networks of XX PGCs are enriched for the transcription factors, TFAP2c, TCFL5, GATA2, MGA, NR6A1, TBX4, and ZFX. Sex-specific enrichment of the forkhead-box and POU6 families of transcription factors was also observed in XY PGCs. Finally, we determined the temporal expression patterns of WNT, BMP, and RA signaling during PGC sex determination, and our discovery analyses identified potentially new cell communication pathways between supporting cells and PGCs. Our results illustrate the diversity of factors involved in programming PGCs toward a sex-specific fate.

The first cell-fate decision is the process by which cells of an embryo take on distinct lineage identities for the first time, thus representing the beginning of developmental patterning. Here, we demonstrate that the molecular chaperone heat shock protein A2 (HSPA2), a member of the 70 kDa heat shock protein (HSP70) family, is asymmetrically expressed in the late 2-cell stage of mouse embryos. The knockdown of Hspa2 in one of the 2-cell blastomeres prevented its progeny predominantly towards the inner cell mass (ICM) fate. In contrast, the overexpression of Hspa2 in one of the 2-cell blastomeres did not induce the blastomere to differentiate towards the ICM fate. Furthermore, we demonstrated that HSPA2 interacted with CARM1 and its levels correlated with ICM-associated genes. Collectively, our results identify HSPA2 as a critical early regulator of the first cell-fate decision in mammalian 2-cell embryos.