The composition of the cell nucleus is highly heterogeneous, with different constituents forming complex interactomes. However, the global patterns of these interwoven heterogeneous interactomes remain poorly understood. Here we focus on two different interactomes, chromatin interaction network and gene regulatory network, as a proof-of-principle, to identify heterogeneous interactome modules (HIMs) in the nucleus. Each HIM represents a cluster of gene loci that are in spatial contact more frequently than expected and that are regulated by the same group of transcription factor proteins. We develop a new algorithm MOCHI to facilitate the discovery of HIMs based on network motif clustering in heterogeneous interactomes. By applying MOCHI to five different cell types, we found that HIMs have strong spatial preference within the nucleus and exhibit distinct functional properties. Through integrative analysis, this work demonstrates the utility of MOCHI to identify HIMs, which may provide new perspectives on 3D genome organization and function. co-regulated genes in GRN can be facilitated by long-range chromosomal interactions (Fanucchi et al., 2013) and chromatin interactome has been shown to exhibit strong correlations with GRN (Kosak et al., 2007;Neems et al., 2016). Indeed, network-based representation of both chromatin interactome and GRN has been suggested to consider different subnuclear components holistically (Rajapakse et al., 2010;Chen et al., 2015). The paradigm of viewing the nucleus as a collection of interacting networks among various constituents can potentially be extended to account for other types of related interactomes in the nucleus. However, whether these interactomes, in particular chromatin interactome and GRN, are organized to form functionally relevant, global patterns remains unknown.In this work, as a proof-of-principle, we specifically consider two different types of interactomes in the nucleus: (1) chromatin interactome -a network of chromosomal interactions between different genomic loci -and (2) a GRN where TF proteins bind to the genomic loci to regulate target genes' transcription. Many studies in the past have analyzed the structure and dynamics of chromatin interactomes and GRNs as well as the coordinated binding of transcription factors on folded chromatin (Rao et al., 2014;Tang et al., 2015;Marbach et al., 2016;Ma et al., 2018;Cortini and Filion, 2018). However, the global network level patterns between chromatin interactome and GRN are still unclear, and algorithms that can simultaneously analyze these heterogeneous networks in the nucleus to discover important network structures have not been developed.Here we aim to identify mesoscale network structures where nodes of TFs (from GRN) and gene loci (from both chromatin interactome and GRN) cooperatively form distinct types of modules (i.e., clusters). We develop a new algorithm, MOCHI (MOtif Clustering in Heterogeneous Interactomes), that can effectively uncover such network modules, which we call heterogeneous interactome mod...