Spatial clustering and common regulatory elements correlate with coordinated gene expression

Many cellular responses to surrounding cues require temporally concerted transcriptional regulation of multiple genes. In prokaryotic cells, a single-input-module motif with one transcription factor regulating multiple target genes can generate coordinated gene expression. In eukaryotic cells, transcriptional activity of a gene is affected by not only transcription factors but also the epigenetic modifications and three-dimensional chromosome structure of the gene. To examine how local gene environment and transcription factor regulation are coupled, we performed a combined analysis of time-course RNA-seq data of TGF-β treated MCF10A cells and related epigenomic and Hi-C data. Using Dynamic Regulatory Events Miner (DREM), we clustered differentially expressed genes based on gene expression profiles and associated transcription factors. Genes in each class have similar temporal gene expression patterns and share common transcription factors. Next, we defined a set of linear and radial distribution functions, as used in statistical physics, to measure the distributions of genes within a class both spatially and linearly along the genomic sequence. Remarkably, genes within the same class despite sometimes being separated by tens of million bases (Mb) along genomic sequence show a significantly higher tendency to be spatially close despite sometimes being separated by tens of Mb along the genomic sequence than those belonging to different classes do. Analyses extended to the process of mouse nervous system development arrived at similar conclusions. Future studies will be able to test whether this spatial organization of chromosomes contributes to concerted gene expression.Author summary: Cellular responses to environmental stimulation are often accompanied by changes in gene expression patterns. Genes are linearly arranged along chromosomal DNA, which folds into a three-dimensional structure. The chromosome structure affects gene expression activities and is regulated by multiple events such as histone modifications and DNA binding of transcription factors. A basic question is how these mechanisms work together to regulate gene expression. In this study, we analyzed temporal gene expression patterns in the context of chromosome structure both in a human cell line under TGF-β treatment and during mouse nervous system development. In both cases, we observed that genes regulated by common transcription factors have an enhanced tendency to be spatially close. Our analysis suggests that spatial co-localization of genes may facilitate the concerted gene expression.