The phenotypic and functional states of a cell are modulated by a complex interactive molecular hierarchy of multiple omics layers, involving the genome, epigenome, transcriptome, proteome, and metabolome. Spatial omics approaches have enabled the capture of information from different molecular layers directly in the tissue context. However, current technologies are limited to map one to two modalities at the same time, providing an incomplete representation of cellular identity. Such data is inadequate to fully understand complex biological systems and their underlying regulatory mechanisms. Here we present spatial-Mux-seq, a multi-modal spatial technology that allows simultaneous profiling of five different modalities, including genome-wide profiles of two histone modifications and open chromatin, whole transcriptome, and a panel of proteins at tissue scale and cellular level in a spatially resolved manner. We applied this technology to generate multi-modal tissue maps in mouse embryos and mouse brains, which discriminated more cell types and states than unimodal data. We investigated the spatiotemporal relationship between histone modifications, chromatin accessibility, gene and protein expression in neuron differentiation revealing the relationship between tissue organization, function, and gene regulatory networks. We were able to identify a radial glia spatial niche and revealed spatially changing gradient of epigenetic signals in this region. Moreover, we revealed previously unappreciated involvement of repressive histone marks in the mouse hippocampus. Collectively, the spatial multi-omics approach heralds a new era for characterizing tissue and cellular heterogeneity that single modality studies alone could not reveal.