Differential gene expression plays a critical role in the development and physiology of multicellular organisms. At a 'systems level' (e.g. at the level of a tissue, organ or whole organism), this process can be studied using gene regulatory network (GRN) models that capture physical and regulatory interactions between genes and their regulators. In the past years, significant progress has been made toward the mapping of GRNs using a variety of experimental and computational approaches. Here, we will discuss gene-centered approaches that we employed to characterize GRNs and describe insights that we have obtained into the global design principles of gene regulation in complex metazoan systems.