Cis-regulatory networks (CRNs) play a central role in cellular decision making. Like every other biological system, CRNs undergo evolution, which shapes their properties by a combination of adaptive and nonadaptive evolutionary forces. Teasing apart these forces is an important step toward functional analyses of the different components of CRNs, designing regulatory perturbation experiments, and constructing synthetic networks. Although tests of neutrality and selection based on molecular sequence data exist, no such tests are currently available based on CRNs. In this work, we present a unique genotype model of CRNs that is grounded in a genomic context and demonstrate its use in identifying portions of the CRN with properties explainable by neutral evolutionary forces at the system, subsystem, and operon levels. We leverage our model against experimentally derived data from Escherichia coli. The results of this analysis show statistically significant and substantial neutral trends in properties previously identified as adaptive in origin-degree distribution, clustering coefficient, and motifswithin the E. coli CRN. Our model captures the tightly coupled genome-interactome of an organism and enables analyses of how evolutionary events acting at the genome level, such as mutation, and at the population level, such as genetic drift, give rise to neutral patterns that we can quantify in CRNs. Reconstructing a CRN from experimental data, elucidating its dynamic and topological properties, and understanding how these properties emerge during development and evolution are major endeavors in experimental and computational biology (1-5).The complexity of CRNs, coupled with observed "unexpected" trends in their properties, such as scale-freeness (6), high degree of clustering (7), and overrepresented subgraphs (3,(8)(9)(10), has led to several hypotheses of adaptive origins and explanations of CRNs and their properties. Central to most of these studies was the use of simplistic graph-theoretic models, such as randomly rewiring the connectivity of a biological network, to serve as a null model for CRN connectivity maps and their properties (11). However, it has been shown that when subjecting CRNs to the various neutral evolutionary forces and tracing their trajectory in time, many of these topological patterns may simply arise spontaneously due to the forces of mutation, recombination, gene duplication, and genetic drift (10, 12). These studies call into question arguments that were made in favor of adaptive explanations for the emergence and conservation of CRN properties (8,(13)(14)(15) and identify important parameters that may significantly affect the evolution of CRNs from a neutral perspective. Specifically (12), they highlighted the role that promoter length, binding-site size, and population size may play in forming certain topological patterns known as motifs. Nonetheless, a lingering question remains: Which specific parts of a CRN arise due to nonadaptive forces and, moreover, can we quantify these patterns to...