Regulatory network topology and the genetic architecture of gene expression

In human populations, most of the genetic variance in gene expression can be attributed to trans-acting expression quantitative trait loci (eQTLs) spread across the genome. However, in practice it is difficult to discover these eQTLs, and their cumulative effects on gene expression and complex traits are yet to be fully understood. Here, we assess how properties of the genetic architecture of gene expression constrain the space of plausible gene regulatory networks. We describe a structured causal model of gene expression regulation and consider how it interacts with biologically relevant properties of the gene regulatory network to alter the genomic distribution of expression heritability. Under our model, we find that the genetic architecture of gene expression is shaped in large part by local network motifs and by hub regulators that shorten paths through the network and act as key sources of trans-acting variance. Further, simulated networks with an enrichment of motifs and hub regulators best recapitulate the distribution of cis and trans heritability of gene expression as measured in a recent twin study. Taken together, our results suggest that the architecture of gene expression is sparser and more pleiotropic across genes than would be suggested by naive models of regulatory networks, which has important implications for future studies of complex traits.