Cellular response to genetic and environmental perturbations is often reflected and͞or mediated through changes in the metabolism, because the latter plays a key role in providing Gibbs free energy and precursors for biosynthesis. Such metabolic changes are often exerted through transcriptional changes induced by complex regulatory mechanisms coordinating the activity of different metabolic pathways. It is difficult to map such global transcriptional responses by using traditional methods, because many genes in the metabolic network have relatively small changes at their transcription level. We therefore developed an algorithm that is based on hypothesis-driven data analysis to uncover the transcriptional regulatory architecture of metabolic networks. By using information on the metabolic network topology from genome-scale metabolic reconstruction, we show that it is possible to reveal patterns in the metabolic network that follow a common transcriptional response. Thus, the algorithm enables identification of so-called reporter metabolites (metabolites around which the most significant transcriptional changes occur) and a set of connected genes with significant and coordinated response to genetic or environmental perturbations. We find that cells respond to perturbations by changing the expression pattern of several genes involved in the specific part(s) of the metabolism in which a perturbation is introduced. These changes then are propagated through the metabolic network because of the highly connected nature of metabolism.bioinformatics ͉ reporter metabolites ͉ metabolic subnetworks L inking the genome to its functioning metabolism is of substantial interest not only in studying human diseases (1) but also for identifying metabolic engineering targets in biotechnological applications (2, 3). Transcriptional analysis represents a high-throughput and genome-wide approach for linking the set of expressed genes to functional metabolism of the cell. Indeed, several studies using genome-wide gene-expression analysis have shown that transcriptional regulation plays an important role in regulating metabolism in response to perturbations (4-6). Although many statistical methods and clustering algorithms provide tools to analyze such transcriptomics data (7-9), these methods seldom provide insight into the regulatory architecture of the metabolic networks without intelligent analysis of the results (up͞down-regulation of genes of interest or correlation between genes of interest). This shortcoming is primarily due to the hypothesis that there may be all-to-all interactions among the genes being analyzed, resulting into many biologically nonsignificant results. One of the ways to address this problem is to integrate known biological interactions, e.g., protein-protein interactions, in the analysis of transcription data (10). Such an approach essentially reduces the degrees of freedom in data analysis by using knowledge of molecular interactions occurring in the cell. The organization and functioning of the cell can be ...