Global transcriptional response of Saccharomyces cerevisiae to the deletion of SDH3

Background Mitochondrial respiration is an important and widely conserved cellular function in eukaryotic cells. The succinate dehydrogenase complex (Sdhp) plays an important role in respiration as it connects the mitochondrial respiratory chain to the tricarboxylic acid (TCA) cycle where it catalyzes the oxidation of succinate to fumarate. Cellular response to the Sdhp dysfunction (i.e. impaired respiration) thus has important implications not only for biotechnological applications but also for understanding cellular physiology underlying metabolic diseases such as diabetes. We therefore explored the physiological and transcriptional response of Saccharomyces cerevisiae to the deletion of SDH3, that codes for an essential subunit of the Sdhp. Results Although the Sdhp has no direct role in transcriptional regulation and the flux through the corresponding reaction under the studied conditions is very low, deletion of SDH3 resulted in significant changes in the expression of several genes involved in various cellular processes ranging from metabolism to the cell-cycle. By using various bioinformatics tools we explored the organization of these transcriptional changes in the metabolic and other cellular functional interaction networks. Conclusion Our results show that the transcriptional regulatory response resulting from the impaired respiratory function is linked to several different parts of the metabolism, including fatty acid and sterol metabolism.


Background
Mitochondrial respiration plays a central role in energy metabolism in eukaryotic cells. It is involved in generating energy, primarily in the form of ATP, upon oxidation of different carbon sources such as ethanol, pyruvate and diverse organic acids. This functionality is conserved across all eukaryotic cells. Consequently, metabolic and regulatory mechanisms governing the mitochondrial energy generation have many implications for the functioning of the cell as a whole. A primary component of mitochondrial metabolism is the TCA cycle which apart from producing NADH needed in the oxidative phosphorylation also supplies precursors for biomass synthesis, e.g. 2-oxoglutarate and oxaloacetate. The coupling between metabolism and oxidative phosphorylation is also reflected in the tight transcriptional regulation of the TCA cycle, for example, as observed in Saccharomyces cerevisiae during the diauxic shift between fermentative and oxidative metabolism. The TCA cycle is also known to be transcriptionally regulated in response to oxygen and carbon substrate concentrations. The succinate dehydrogenase complex (Sdhp) serves as a link between the TCA cycle and electron transport chain. Specifically, FADH 2 produced during oxidation of succinate to fumarate acts as an electron donor for ubiquinone. Thus, the flux through the succinate dehydrogenase reaction is directly coupled with the respiratory capacity of the cell. The respiratory chain itself is believed to be transcriptionally regulated by oxygen and heme concentrations, and several other not well characterized mechanisms, e.g. in response to osmotic stress [1].
Defects in the respiratory chain and related energy metabolism have been shown to be key factors in inducing several human diseases, including diabetes, obesity and certain types of cancers [2][3][4][5][6]. Consequently, respiratory metabolism has been a major focus for studying varying diseases, and there is much interest in using appropriate eukaryotic model systems that can be used for functional studies. The yeast S. cerevisiae has been widely used as a model eukaryotic organism as many cellular processes are conserved to humans, and comparatively large amounts of genomic, metabolomic, and proteomic data is readily available. Furthermore, S. cerevisiae is industrially used for producing several commodity and high added-value compounds, such as bio-ethanol and insulin. Consequently, there is much interest in developing rational design strategies for metabolically engineering yeast to improve the production of desired compounds. Previous studies have suggested the TCA cycle and energy metabolism as primary targets for certain metabolic engineering problems [7,8].
The Sdhp is a tetramer consisting of two soluble subunits responsible for the dehydrogenase activity, and two hydrophobic subunits that anchor the catalytic subunits to the mitochondrial inner membrane. The SDH3 gene in S. cerevisiae codes for the cytochrome b component of the complex, and as shown previously [9], disruption of the gene leads to a severe growth defect on non-fermentable carbon sources demonstrating its major role in the complex assembly and function. In light of the importance of Sdh3 for improved understanding of human diseases and metabolic engineering applications, we characterized the physiology of sdh3Δ mutant and performed genome wide transcription analysis. The cultivation type was chosen to be batch fermentation with glucose as the sole carbon source so as to have glucose-repressed conditions.

Growth conditions
All cultures were grown in minimal media prepared according to Verduyn et al. [13] supplemented with 150 mg/liter of uracil to complement auxotrophy for the deletion mutant strain. Before each experiment, cells from 15% glycerol stock preparations were plated on YPD medium (10 g/liter yeast extract, 20 g/liter of peptone, 20 g/liter glucose) and incubated at 30°C for 24-48 h. A single colony was then used to inoculate precultures.

Batch fermentation
Batch experiments for physiological characterization of the strains were carried out in a Braun Biostat-B reactor with a working volume of 1.5 liter and 10 g/liter initial glucose concentration. A constant pH of 5 was maintained via automated addition of 2 M NaOH. Aerobic conditions were established by supplying air at the rate of 1.5 liter/min whereas sparging with pure N 2 established anaerobic conditions. The concentration of CO 2 and O 2 in the exhaust gas was determined by the use of a Brüel & Kjaer acoustic gas analyzer (Brüel & Kjaer, Denmark). For anaerobic cultivation the medium was supplemented with ethanol-dissolved ergosterol and Tween 80 (Sigma-Aldrich) to a final concentration of 0.01 g/liter and 0.42 g/liter, respectively.
Five liter in-house-manufactured bioreactors with 4 liter working volume were used to carry out aerobic batch cultivations for microarray experiments. The initial glucose concentration was increased to 40 g/liter and therefore the concentration of all medium components was accordingly adjusted. The cultivations were carried out in at least triplicate.
Chemostat cultivation C-limited chemostat experiments were also performed in Braun Biostat B with the same process variables as described above for the aerobic batch, at a dilution rate of 0.027 h -1 for the sdh3Δ strain and 0.05 h -1 for the wild type strain (WT). Steady state samples for extracellular metabolite analysis were collected after 3 volume changes.

Inhibitor experiments
Shake flask experiments were also performed to analyse the effects of 3 specific inhibitors of the respiratory chain. The SDH3 mutant and the reference strain were grown in minimal media and treated with antimycin (1 μg/ml), oligomycin (3 μg/ml) or carbonyl cyanide mchlorophenylhydrazone (CCCP) (4.1 μg/ml) (Sigma-Aldrich) during the early exponential phase.

Sampling
During the course of all cultivations 5 ml samples were withdrawn from the bioreactor at regular time intervals for the determination of dry cell weight and extracellular metabolites. Glucose, ethanol, pyruvate, succinate, acetate, fumarate, citrate, oxaloacetate and malate were quantified by HPLC using a Biorad HPX-87H column.

RNA isolation and microarray analysis
Duplicate samples for microarray analysis were collected during the aerobic batch experiments with the sdh3Δ strain when the glucose concentration had reached 20 ± 2 g/liter. Culture aliquots of 20 ml were transferred to falcon tubes already filled with 25 ml of crushed ice and immediately centrifuged at 4000 rpm for 5 minutes. The pellets were rapidly frozen in liquid nitrogen and stored at -80°C until further use. Total RNA was extracted by using the Fast Pro Red Kit (Obiogene, USA) according to the manufacture's instructions with minor modifications. cRNA was synthesized as described in the Affymetrix GeneChip Expression Analysis Manual, after which 15 μg were hybridized to Yeast Genome S98 oligonucleotide arrays (Affymetrix). Microarrays were scanned in Agilent Gene Scanner (Affymetrix), following which raw data processing was performed with the Microarray Suite software, version 5, with a global scaling factor for target intensity equal to 500. Average difference values, representing the absolute hybridization intensities were then calculated for each probe set. The microarray data was first normalized by using dChip software suit (

Reporter algorithm
In order to identify key proteins involved in the cellular response to SDH3 deletion, we used a systems approach where bio-molecular interaction networks are used as data integration scaffolds [14,15]. Three different interaction networks were used in this study, viz. metabolic network [16], regulators/Transcription Factors (TFs)-'regulated genes' network [15], and 'Gene Ontology' (GO)-'associated genes network' [15]. All three networks can be conceptually seen as a graph where the feature nodes (metabolites, GOs or TFs) are connected with the corresponding functionally linked gene nodes. The p-values for the individual genes were then mapped on each of these networks and the features with significant transcriptional response in their neighbour genes were identified as reporter features, i.e., reporter metabolites, reporter GOs and reporter TFs. The scoring method used for estimating the significance score for the features is summarized in the following three equations [15]: First, p-value for each gene i in the network, p genei , was converted into a z-score by using the inverse normal cumulative distribution function (cdf -1 ). The resulting zscores approximately follow a standard normal distribution. The score for each feature j (z featurej ) was then calculated as the average of the z-scores of its neighbour nodes (equation 2). To assess the significance of each z featurej , it was corrected for the background distribution of z scores, by subtracting the mean (μ N ) and then dividing by the standard deviation (s N ) of 10000 random aggregates of size N selected from the same data; N being the number of neighbours for the feature under consideration. The resulting reporter features can be seen as nodes in the bio-molecular interaction networks around which there are substantial transcriptional changes occurring. These nodes are a metabolite interacting with different enzymes, a transcription factor regulating the transcription of several genes or a common annotation for a group of genes, e.g. GO annotation.
The reporter scores calculated in this study do not require a priori selection of a p-value cut-off at the individual gene level. It may be also of interest to observe the enrichment in the number of significantly changed genes around a particular feature, whereas the cut-off level at the individual gene level is chosen a priori.
The biological meaning of such scoring method is discussed by Oliveira et al [15]. The software for calculating such enrichment scores is freely available at: http://www.cmb.dtu.dk/reporters.

Promoter analysis
A bioinformatics analysis was performed by using Regulatory Sequence Analysis web-tools http://rsat. scmbb.ulb.ac.be/rsat/ [17] in order to identify potential regulatory elements in the promoter region (800 bp region upstream of each start codon) of the genes that showed a significant change in expression following the SDH3 deletion.

Physiological characterization
In order to study the physiological consequences of SDH3 deletion, a comparison between the reference and mutant strain was performed under well-controlled aerobic and anaerobic batch fermentations in minimal medium with glucose as the sole carbon source. The sdh3Δ mutant showed a lower maximum specific growth rate (0.25 h -1 ) in aerobic growth conditions as compared to the reference strain (0.31 h -1 ). Furthermore, the mutant was incapable of utilizing ethanol as a sole carbon source, as also been reported previously [9]. A similar phenotypic difference was observed in anaerobic batch cultivations, where the maximum specific growth rates were 0.27 h -1 and 0.34 h -1 for the mutant and the reference strain, respectively. Since, under both aerobic and anaerobic conditions the cells are mainly having fermentative metabolism, the lower specific growth rate observed in both conditions are likely to be due to the same mechanisms.
The yield coefficients of biomass and secreted metabolites with respect to consumed glucose are listed in Table  1. We observed a slight increase in succinate production in the mutant strain; yield of 0.0033 g/g-glucose was obtained in anaerobic growth conditions whereas no (or undetectable amounts) succinate was produced by the reference strain. The physiology of the mutant was also characterized in aerobic, glucose limited chemostat cultivations at low dilution rates. The results indicate that the respiratory metabolism is non-functional in the mutant, as pronounced fermentative metabolism was observed at dilution rates as low as 0.027 h -1 .

Inhibition of respiration
Shake flask experiments were performed on minimal media to study the effects of three different oxidative phosphorylation inhibitors (antimycin, oligomycin and CCCP) on the SDH3 mutant. Antimycin inhibits the electron transfer and thereby prevents the oxidation of both NADH and succinate. Oligomycin specifically binds to the F 0 subunit of the ATP-synthase complex inhibiting ATP synthesis and ATP-synthesis associated respiration. The last inhibitor used, CCCP, is an uncoupling agent that allows free movement of protons across the mitochondrial inner membrane causing a collapse in the H + gradient that is required for ATP synthesis by the F 0 F 1 -ATPase. The addition of oligomycin and antimycin had no effect on the specific growth rate and biomass production in the sdh3Δ strain while they considerably reduced the final biomass yield in the reference strain and impaired growth on ethanol. The major inhibitory effects in both strains were observed after treatment with CCCP which resulted in reduced specific growth rate and biomass yield in the reference strain, while it completely inhibited growth of sdh3Δ.
Global overview of transcriptional response to SDH3 deletion To understand the cellular processes denominating the observed physiological differences in aerobic batch cultivations we performed a genome-scale transcriptional analysis of the sdh3Δ strain (see methods). We used the previously published data for the reference strain [18] for comparative analysis. Approximately 25% of the yeast genes showed significant changes in their expression levels (p-value ≤ 0.05), indicating a global transcriptional response. After applying a correction for multiple testing by using q-test [19], 307 genes were deemed as significantly changed. We also applied rankproduct test [20] in order to identify the genes that were significantly up/down regulated, yielding additional 127 genes (additional file 1). The distribution of these 434 genes in different functional categories as defined by the MIPS database http://mips.gsf.de/genre/proj/yeast/ is shown in Table 2. The next largest group to the 'unclassified proteins' is the group of genes belonging to RNA synthesis and processing. Moreover, several genes belonging to intracellular protein and mRNA transport (including vacuolar and endoplasmic reticulum related), carbon metabolism, cell cycle, cellular response and sensing, ribosomal biogenesis and DNA metabolism have significantly changed expression levels in the sdh3Δ mutant (Table 2). In general, the majority of the genes are down-regulated in each category, indicating a global repressive effect of the SDH3 deletion. To investigate on the possibility that the slower growth rate of the mutant may lead to decreased expression of several genes, we used the gene expression data of S. cerevisiae grown at several different specific growth rates in chemostat cultivations [21] to compare against the changes observed in the sdh3Δ mutant. It should be noted that due to the different cultivation conditions (i.e. batch vs. chemostat) it is not possible to draw direct conclusions on the growth rate effects in this fashion. Moreover, the effects of growth rate on gene expression are inherently linked with the metabolic state and cell cycle and thus the problem of differentiating between growth rate and perturbation effects is complex [22,23]. Nevertheless, we decided to use a simple extrapolative analysis by using the data from Regenberg et al [21]. Following linear fitting of expression values of each gene with growth rate, we checked whether the direction and magnitude of expression changes in the sdh3Δ mutant concur with those predicted by this linear model. Based on this analysis we postulate that most of the observed changes can not be directly and solely attributed to the difference in growth rate.
Although the classification of differentially expressed genes according to the MIPS functional categories provides a birds-eye view of the transcriptional response, we further investigated whether there are certain biological processes, enzyme functions or cellular locations that responded more significantly than the rest. For this purpose, we integrated the Gene Ontology network with transcription data in order to obtain reporter-Gene Ontology [15] terms. This integration algorithm uses the data from all genes without a priori selection of significant genes, and is capable of identifying weak but biologically significant changes [15]. Thus, reporter-GOs indicate the key biological processes/functions/cellular components that collectively respond significantly to the perturbation (for the complete list of reporter GOs see additional file 2). Reporter-GO terms for the sdh3Δ mutant (Table 3) clearly reflect the large transcriptional changes in nucleolus related genes (148 genes), indicating altered ribosome synthesis and transport due to SDH3 deletion. These changes also appear to be The yield coefficients for biomass, ethanol, glycerol, pyruvate and acetate were determined in the exponential growth phase from linear regression of their concentration as a function of the residual glucose concentration. Succinate yield was calculated at a fermentation end-point. Y sx : grams of biomass per g of glucose consumed. Y ssuc : grams of succinate per g of glucose consumed. Y se : grams of ethanol per g of glucose consumed. Y sac : grams of acetate per g of glucose consumed. Y sg : grams of glycerol per g of glucose consumed. Y spyr : grams of pyruvate per g of glucose consumed. μ max : maximum specific growth rate in h -1 . TA: trace amount. extended to other related processes such as ribosomal assembly and maintenance, nuclear rRNA and mRNA processing and protein localization. Appearance of "nucleus" (1071 genes) as one of the top-scoring reporter-GO terms shows that many of the nuclear genes are responding to the changes induced by defective respiration. This set also includes genes related to meiosis, DNA replication initiation, and DNA repair. These changes might be one of the reasons for reduced growth rate of the deletion mutant, particularly as they are also found to being effected upon changing growth rate in chemostat cultures [21].
Around 20 genes sharing the GO term "transcriptional regulator activity" also scored high. These genes may provide insight into the regulatory pathways that might be responsible for observed transcriptional response in this variety of processes. This set involves not only global transcriptional regulators such as Ume6p, but also specific regulators, especially related to cell cycle and DNA repair.

Reporter TFs analysis
We used a network reconstruction depicting each known yeast transcription factor or regulatory protein connected to all genes known to being affected by these proteins. Using this graph we could identify Reporter Transcription Factors (or more correctly Reporter Regulators) indicating a significant collective change in the expression of the genes regulated by them (Table 4). For the complete list of reporter TFs see additional file 2.
Reporter TFs highlight the regulatory pathways affected following a perturbation, and thus uncover the functional links between the perturbation and the regulatory mechanisms invoked in the cell. Since many TFs/ regulators are regulated at post-transcriptional level, the expression of the reporter TFs/regulators themselves may or may not change per se. Nevertheless, identification of a particular protein as reporter-TF is a strong indication for the change in its activity following SDH3 deletion. The most marked reporter-TF emerging from this analysis is the Hap complex, for which the expression levels of Hap4 and Hap2 are also significantly changed. This confirms the association of Hap4 to the control of respiration [24]. Other reporter-TFs represent potential mediators of the observed global transcriptional changes as their involvement range from the control of metabolism (GCR1, THI2, and SLN1) to cell cycle (FAR1, SWI6, YHP1 and YOX1).

Promoter analysis
We searched for conserved sequence motifs in the upstream region of the genes with significant transcriptional response. With all the 434 genes queried together, no significantly over represented motif was found. However, when up and down regulated genes were queried separately we identified significantly overrepresented sequence motifs. About 40% of up-regulated genes contain a 6-mer "TCGTCA" in their promoter region, while this motif is characteristic for only 20% of the genes in the whole genome. 29% of the significantly down-regulated genes also contain this motif in their  [25]. Binding motifs for certain other transcriptional regulators such as Hsf1p (Heat shock transcription factor) and Ino4p (Inositol requiring) were also identified in the upstream regions of significantly down-regulated genes. Expression levels of ABF1 and INO4 were altered significantly in the SDH3 deletion mutant.

Reporter Metabolite and subnetwork analysis
To gain further insight into the transcriptional response to the deletion of SDH3, we integrated the transcriptome data with the genome-scale yeast metabolic model [16] and thus identified so-called reporter metabolites [14] (methods). This method is based on the hypothesis that expression of genes encoding enzymes must be coordinated in order to satisfy the metabolic demands of the cell and corresponding stoichiometric constraints. The algorithm enables identification of so-called reporter metabolites (metabolites around which the most significant transcriptional changes occur) and enzyme subnetworks that respond significantly to genetic or environmental perturbations [14], in this case deletion of SDH3. Table 5 shows the key reporter metabolites (for the complete list see additional file 2) that were identified and it highlights metabolites involved in the TCA cycle (Isocitrate, Oxalosuccinate, FADH 2 ), in the ergosterol biosynthesis pathway (zymosterol) and in phospholipid biosynthesis (ethanolamine).

Discussion
Results from the physiological characterization indicate that the SDH3 deletion mutant shows a growth defect under both aerobic and anaerobic conditions. This observation was partly expected for aerobic conditions and could be attributed to a major reduction in the respiratory capacity of the mutant strain, which is also supported by the carbon limited chemostat experiments. However, as the flux through the TCA cycle of the reference strain at the studied growth conditions is very low [26] due to glucose repression of respiration, it clearly points out that Sdh3 plays an additional role than being involved in respiration. This hypothesis is confirmed by our finding that there is a similar reduction in specific growth rate under anaerobic conditions. This additional role of Sdh3 may be different at the two different growth conditions. Thus, at anaerobic growth conditions it is likely that Sdh3 simply plays a role in redox balancing, whereas at aerobic growth conditions where there is predominantly fermentative metabolism, it may still be important to have functional respiration in order to prevent oxidative damage in the cells. Treatment with respiratory inhibitors antimycin and oligomycin had no effect on growth and final biomass yield of the sdh3Δ mutant while it severely affected biomass production of the reference strain. These additional experiments with mitochondrial inhibitors demonstrated that SDH3 deletion is sufficient to impair the electron transport system preventing the oxidation of NADH and succinate. However, the general down regulation of respiratory genes clearly implies that the cells are not attempting to restore their respiratory capacity. This is in line with the results obtained for respiratory deficient petite cells grown on non-repressive raffinose [27].
Off gas profiles during the batch experiments showed that the oxygen uptake rate in sdh3Δ is~80% lower compared to the reference strain. Since oxygen is not used as final electron acceptor in the respiratory chain, as inferred from the experiments with mitochondrial inhibitors, the oxygen consumed is used in other cellular pathways. In yeast species non respiratory oxygen consumption is favoured when respiratory functions are impaired or absent, however, little is known about the non-respiratory oxygen uptake of aerobically growing cells.
During anaerobiosis, S. cerevisiae strongly increases glycerol production to provide for non-respiratory oxidation of NADH to NAD + , and consistent with the hypothesis of an altered redox metabolism originating from the deletion of SDH3 and consequent impairment of oxidative phosphorylation, a higher glycerol production was observed in aerobic cultivations of the mutant strain (similar yield as obtained in anaerobic cultivation). Thus, the deletion strain seems to readjust its metabolic fluxes mimicking growth under anaerobic conditions.
In an attempt to uncover the molecular mechanisms behind the observed global transcriptional response, we searched for significantly over-represented sequence motifs in the promoter regions of genes with significantly changed expression in the SDH3 mutant. This analysis revealed a binding site for Abf1p in 40% of the significantly up-regulated genes. Several different chromatin-related events such as DNA replication [28], gene silencing [29], chromatin remodeling [30], nucleotide excision repair [31], and gene activation and repression [25] involve direct binding of Abf1p. The genes that are reported to be transcriptionally regulated by Abf1p are involved in a multitude of cellular processes including carbon source regulation [32,33] nitrogen utilization [34,35] sporulation [36], meiosis [37], and ribosomal biogenesis [30,38]. Since these processes are also significantly affected in the SDH3 mutant, Abf1p appears to be one of the major mediators of the transcriptional response following the SDH3 disruption and consequent respiratory impairment. This is further supported by the fact that ABF1 transcription is significantly down regulated by more than four-fold in the SDH3 mutant, whereas under usual conditions ABF1 is found to be present in abundance [28,39]. Down regulation of ABF1 under anaerobic conditions also supports its role as a mediator of transcriptional changes in response to the absence of respiration. The binding motif recognized by Ino4p was also identified in the promoter analysis. Ino4p is a known regulator of genes involved in phospholipid synthesis [40,41], some of which also showed altered expression in the SDH3 mutant (ACC1, CHO1, EKI1 and DPP1).
The reporter metabolites indicate specific parts of metabolism where significant transcriptional regulation is exerted, either to maintain homeostasis and/or to adjust the metabolite levels to altered demands. Two of the top scoring reporter metabolites suggested an interesting connection between mitochondrial impairment and lipid metabolism. Zymosterol is the  Only the top-20 reporter metabolites are shown. Number of neighbours indicates the number of enzyme-coding genes whose products catalyze reactions involving that metabolite either as a substrate or as a product. immediate precursor of ergosterol, the major yeast sterol, which has a functional similarity to cholesterol in higher eukaryotic membranes, and is only synthesised in yeast under aerobic conditions. Our data show up-regulation of AUS1, an important sterol uptake gene normally expressed under anaerobic conditions, and down regulation of several key genes in the ergosterol biosynthesis pathway (ERG11, ERG25, ERG27, ERG5, ERG9) due to the deletion of SDH3. This information indicates the presence of a link that connects mitochondrial respiration and/or TCA cycle activity to the sterol biosynthetic pathways. We postulate that one of the mediators of this functional link is Hap1p which is known to positively regulate the AUS1 expression that is up-regulated in the SDH3 deletion mutant. Since the effect of mitochondrial impairment on lipid metabolism has been a focus of several studies related to metabolic diseases such as diabetes, our results may help in identifying similar regulatory mechanisms in human cells.
Ethanolamine is a precursor for the synthesis of phosphatydilethanolamine (PE), one of the major phospholipids components in cellular membranes of S. cerevisiae, and its role is crucial for growth when mitochondrial function is required [42]. In humans impairment or alterations in the phospholipids metabolism are suggested to represent important aspects of certain diseases [43]. Expression of several genes involved in PE biosynthesis such as EKI1 (ethanolamine kinase/choline kinase), EPT1 (dyacil-glycerol-ethanolamine-phosphotrasferase), and CHO1 (phosphatidylserine synthase) is known to be transcriptionally regulated by Ino4p, and they are all down-regulated in the sdh3Δ strain. As found in the promoter analysis, the binding motif of Ino4p was significantly over-represented in the promoter region of some of the significantly downregulated genes, making it another potential mediator of the signal between respiratory chain and lipid metabolism.
Certain other metabolites including isocitrate, oxalosuccinate, and FADH 2 involved in the TCA cycle also emerge as reporter metabolites, and together with the altered expression of many genes related to the TCA and glyoxylate cycles (Figure 1) they possibly suggest alterations of Acetyl-CoA metabolism (see additional file 4 for supporting information).
In accordance with our results, transcription profiling of a collection of mutants defective in each of the TCA cycle genes also showed transcriptional changes in a large number of genes, albeit under non-repressive conditions [44]. Taken together, our results suggest that nuclear gene signalling is responsive to TCA cycle dysfunction. Many of the transcriptional changes following the  deletion of SDH3 are located in fatty acid and sterol metabolism. What is the regulatory link that connects the respiratory dysfunction to changes in fatty acid metabolism and other pathways? Towards investigating this question, inferences and postulates based on our analysis are summarized in the qualitative regulatory model represented in Figure 2. In S. cerevisiae the heme-activated complex Hap2/3/4/5 is known to regulate the expression of many of the TCA cycle and respiratory chain genes, although neither the heme nor the oxygen regulation of the Hap-complex is clearly understood. The Hap-complex appears as one of the transcription factors with altered expression level in the SDH3 mutant. Moreover several of the genes regulated by the Hap-complex were also responsive to this change as indicated by Reporter-TF analysis. Thus, it may be hypothesized that the deletion of the SDH3 gene (and perhaps respiratory deficiency in general) might be in the first place responsible for causing the down-regulation of genes coding for the Hap-complex, which in turn would explain the down-regulation of respiratory genes and some of the TCA cycle genes. The Hap-complex has also been found to downregulate the expression of heme biosynthesis pathway [24]. Consequently, the up-regulation of HAP1, another hemeresponsive transcription factor, can be attributed to reduced heme biosynthesis. Hap1p is involved in the regulation of several pathways including sterol biosynthesis and uptake, and is essential during anaerobic or heme-deficient growth.
AUS1 (coding for sterol uptake protein) and several genes involved in the ergosterol, fatty acid and phospholipid biosynthesis show differential gene expression in the SDH3 mutant that is consistent with the changes observed in the reference strain under anaerobic conditions [45]. Our analysis suggests that Ino4p may play a key role in the regulation of these genes following SDH3 deletion, which is consistent with the known regulatory targets of Ino4p. However, it is unclear how the Ino4p regulation is triggered.

Conclusion
Overall, we present a genome-wide transcription analysis of SDH3 deleted S. cerevisiae mutant grown under glucose-repressed conditions. By using an integrative systems biology approach we reconstructed a qualitative regulatory map that links mitochondrial respiration with fatty acid and sterol metabolism. Several bioinformatics and systems biology algorithms were employed to elucidate the transcriptional hot-spots in different cellular networks, including metabolic networks and transcriptional regulatory networks that significantly respond to the deletion of SDH3. Several regulatory proteins including the Hap-complex, Hap1, Ino4, Abf1 and Sln1 appear to play an important role in mediating the cellular response following respiratory defect at the succinate dehydrogenase node.

Additional File 3
Promoter analysis. List of all the significantly changed genes that have Abf1 binding motif in their promoter region.