Bio-informatics analysis of a gene co-expression module in adipose tissue containing the diet-responsive gene Nnat

Background Obesity causes insulin resistance in target tissues - skeletal muscle, adipose tissue, liver and the brain. Insulin resistance predisposes to type-2 diabetes (T2D) and cardiovascular disease (CVD). Adipose tissue inflammation is an essential characteristic of obesity and insulin resistance. Neuronatin (Nnat) expression has been found to be altered in a number of conditions related to inflammatory or metabolic disturbance, but its physiological roles and regulatory mechanisms in adipose tissue, brain, pancreatic islets and other tissues are not understood. Results We identified transcription factor binding sites (TFBS) conserved in the Nnat promoter, and transcription factors (TF) abundantly expressed in adipose tissue. These include transcription factors concerned with the control of: adipogenesis (Pparγ, Klf15, Irf1, Creb1, Egr2, Gata3); lipogenesis (Mlxipl, Srebp1c); inflammation (Jun, Stat3); insulin signalling and diabetes susceptibility (Foxo1, Tcf7l2). We also identified NeuroD1 the only documented TF that controls Nnat expression. We identified KEGG pathways significantly associated with Nnat expression, including positive correlations with inflammation and negative correlations with metabolic pathways (most prominently oxidative phosphorylation, glycolysis and gluconeogenesis, pyruvate metabolism) and protein turnover. 27 genes, including; Gstt1 and Sod3, concerned with oxidative stress; Sncg and Cxcl9 concerned with inflammation; Ebf1, Lgals12 and Fzd4 involved in adipogenesis; whose expression co-varies with Nnat were identified, and conserved transcription factor binding sites identified on their promoters. Functional networks relating to each of these genes were identified. Conclusions Our analysis shows that Nnat is an acute diet-responsive gene in white adipose tissue and hypothalamus; it may play an important role in metabolism, adipogenesis, and resolution of oxidative stress and inflammation in response to dietary excess.


Background
The WHO estimates there will be 2.3 billion overweight or obese adults world-wide by 2015; including 60% of adults and alarmingly 25% of children in developing countries http://www.who.int. This overweight and obesity epidemic is predicted to affect half the global population by 2025 [1,2], largely because of increased urbanization, adaptation of the western diet, and sedentary lifestyle, particularly in China and India.
A serious and nearly ubiquitous complication of obesity is defective insulin signalling in target tissuesskeletal muscle, liver, adipose tissue and brain. Insulin resistance unites obesity with a major disease clusterdiabetes, fatty liver and heart disease -through compensatory hyperinsulinaemia leading to pancreatic b-cell failure, and complicating dyslipidaemia and hypertension, major risk factors for CVD. T2D affects 150 million people world-wide, and is estimated to double by 2025 [1,2]. CVD is predicted to be the leading cause of death globally by 2020 [3,4]. Dementia and cancer are also more common in the obese [4,5]. Therefore, understanding the mechanisms underlying obesity and insulin resistance is of great significance, as the need for new and more effective treatments grows.
Insulin resistance in obesity is causally associated with macrophage influx and low-grade, metabolic inflammation in adipose tissue. Macrophage-driven chronic inflammation in adipose tissue from the obese is considered to be triggered by endoplasmic reticulum (ER) and oxidative stress, together with hypoxia consequent on adiposity expansion and inadequate blood supply. Cellular stress leads to activation of stress kinesis and inflammatory effectors such as Jnk and Nfkb, with systemic release of cytokine and fatty acid mediators, which inhibit insulin signalling in target tissues [1][2][3][4][5][6][7].
While appetite, energy expenditure and metabolism are all critically regulated by hypothalamic brain circuits, white adipose tissue (WAT) is the major storage site for fat in the form of triglycerides, and therefore forms a repository balancing energy intake and expenditure. The overall adipose tissue burden of an individual is a combination of adiposity number (i.e. adipogenesis) and lipid load per cell (i.e. lipogenesis), and both processes may contribute to the adipose tissue mass in obesity.
Nnat is a paternally expressed imprinted gene first discovered in the developing brain and it is implicated in antenatal brain development and foetal growth [6][7][8][9]. In adult rodents and humans Nnat is highly expressed in neurons in many brain regions, including hypothalamic nuclei concerned with appetite regulation, and is concerned in nutrient sensing, and responds to feeding/ fasting and leptin [10,11]. Nnat is also highly expressed in adult white adipose tissue and aortic endothelial cells, where expression is increased in obese and diabetic rodents [12,13]. Nnat is moreover expressed in pancreatic b-cell lines, where expression is glucose-dependent, and reduced in diabetic rodents [14][15][16]. By contrast Nnat expression in WAT is substantially reduced in lipodystrophic, in lean S6kko mice and in mice fed conjugated linoleic acid, which induces weight loss [16][17][18]. Recently a consistent association between single nucleotide polymorphisms in the human Nnat gene and severe obesity in both children and adults has been demonstrated [10]. This observation suggests that altered Nnat function in adipose tissue or the hypothalamic regions regulating energy homeostasis may be involved in the pathogenesis of human obesity.
Nnat has three exons, and resides within a large intron of Blcap [7,8], tissue specific enhancers lie well upstream of the CpG-rich Nnat promoter. Knowledge about the upstream regulation of Nnat is confined to the β-cell and control by NeuroD1, a transcription factor important for neuronal and endocrine cell differentiation and survival [14].
Knowledge concerning downstream effecter mechanisms through which Nnat operates, and its physiological functions are limited. Nnat is an amphipathic proteolipid. It has two isoforms, α of 81 and b of 54 amino acids, generated by alternative splicing of the middle exon, encoding a membrane-binding domain [13][14][15]19]. Both isoforms reside in the endoplasmic ER membranes, the b-isoform is implicated in ER stress and b-cell apoptosis [15], it is considered to regulate ER calcium ATPase through homology to other proteolipid. Nnat expression in vitro is associated with calcium-induced 3T3-L1 cell adipogenesis; glucose-stimulated, calcium-induced, insulin secretion from b-cells; activation of Pi3k, Erk, mTor and calcium signalling in medulloblastoma; Nfkb-regulated inflammation in aortic endothelial cells and protection against mitochondrial toxins and ionophores in PC12 cells and resolution of injury after ischemic/reperfusion of the kidney [12][13][14]19,20].
The physiological roles of Nnat in the brain, pancreas, adipose and other tissues are not understood. Through genome-wide gene expression analysis we sought to improve our understanding of the mechanisms of Nnat regulation and its biological functions. First we carried out an integrated bioinformatics analysis of the conserved Nnat promoter in seven species and derived a candidate list of potential Nnat transcription factors. Then we interrogated multi-tissue gene expression data to look for primary evidence on the applicability of potential Nnat transcriptional regulators. We then used Gene Set Enrichment Analysis [21] to understand which biological functions are most highly correlated with Nnat expression. Finally, we generated a mouse proteinprotein interaction network in adipose tissue to probe the biological relationships among genes forming a co-expression network with Nnat.

Microarray datasets and meta-analysis
In this study, we used four different microarray gene expression datasets from WAT of C57BL/6J mice, as described below, to complete a meta-analysis. The first dataset designated BAIR was obtained as part of the Biological Atlas of Insulin Resistance (BAIR, http:// www.bair.org.uk) project, a multi-lab collaborative effort to define gene expression changes involved in insulin signalling and resistance. Male mice were weaned at 3 weeks of age and fed a control chow diet containing 5% w/w fat, 19% protein and 3.5% fibre (S&K Universal Ltd, Hull, UK). At 5 weeks of age, a group of mice was transferred to 40% w/w high-fat diet, of which saturated fat is 40%, containing 32% pig fat, 8% casein fat, 19% protein, 21% glucose (Special Diet Services, Witham, UK) ad libitum, whereas the control group remained on chow. Mice fat-fed for 2 days, 8 days, 3 weeks and 15 weeks, and with age-matched controls, were killed and epididymal WAT was obtained for microarray gene expression studies, employing Affymetrix mouse genome 4302 microarray. From these four time points a total of 40 Affymetrix gene expression raw CEL files were produced and RMA normalized by using the Affy package http://www.bioconductor.org; present/absent call information was also calculated based on Affymetrix MAS5 algorithm. This dataset has been deposited in EBI ArrayExpress Database with accession number E-BAIR-12 http://www.ebi.ac.uk/microarray-as/ae/.
The second dataset was obtained from the Gene Expression Omnibus (GEO, http://www.ncbi.nlm.nih. gov/geo/, accession ID: GSE4671), which contains 28 Affymetrix mouse genome 4302 chips. Nine-week old male mice in this experiment were fed with either a control chow diet or a chow diet that also contained 0.5% conjugated linoleic acid (CLA). At each of seven time points (day 1, 2, 3, 4, 7, 10, 17), retroperitoneal WAT was obtained from two control mice and two test mice. Microarray data were processed using GCRMA algorithm.
The third dataset was also acquired from GEO with accession number GSE8831, in which 20 female and 15 male C57BL/6J mice, fed by ad libitum, varied in body weight and insulin sensitivity were studied. Fasting blood glucose and serum insulin concentrations were measured 2-4 days prior to collection of adipose tissue. RNA from perigonadal adipose tissue was extracted and processed for analysis on Affymetrix Mu74Av2 microarrays. RMA normalization was applied to the 35 raw CEL files, and MAS5 present/absent call information was calculated as well.
The fourth dataset consisted of 13 Affymetrix MG-U74av2 chips for epididymal WAT extracted from six-week old C57BL/6J male mice from the Diabetes Genome Anatomy Project DGAP http://www.diabetesgenome.org. There were four high-fat feeding mice (derives 55% calories from fat, 21% calories from protein, 24% calories from carbohydrates), nine low-fat feeding control mice (derives 14% calories from fat, 25% calories from protein and 61% calories from carbohydrates). Gene expression data was RMA normalized with present/absent call information.
In order to identify the gene co-vary with Nnat in WAT, we calculated Pearson correlation coefficients between the expression of Nnat and all other genes in each of the above four datasets separately, and chose genes showing the highest correlation with Nnat expression. Subsequently, rather than merely combining the p-values of the Pearson correlation coefficients, we used weighted unbiased Pearson statistical values, i.e. we normalized the statistical values to approximate a normal distribution with unit variance and zero mean. Some necessary filters were designed for quality control: (a) probesets with expression value lower than 100 or (b) absolute Pearson value less than 0.2 or (c) present/ absent call less than 30% were excluded. Each gene was represented by the probesets which had the largest absolute Pearson correlation in each dataset in order to remove multiple annotated probesets for the same gene. In this way we obtained 2841, 3061, 2148 and 1201 unique genes within the four datasets respectively, giving 6120 unique genes in total.
The following formula was applied to calculate the combined Z-score A two-tailed P-value can be derived by the normal distribution for each gene from this combined Z score . Bonfferroni corrected FDR < 0.05 was applied, yielding the top 34 genes which are most strongly co-expressed with Nnat. For comparison we also applied Fisher's Z-transform to combine the transformed Pearson correlation coefficients in the meta-analysis, we got similar results. For the multi-species analysis, seven genes, Acp5, Tcta, Znrf2, Serpina3c, S3-12, Sgce and Lrg1, which are not conserved in rat, were excluded.

Promoter analysis
To delineate the Nnat promoter, we identified the conserved non-coding region of the Nnat locus in seven mammalian genera (Pan, Bos, Canis, Mus, Sus, Rattus and human). This corresponds to 1200 bp upstream and 800 bp downstream of the most utilized transcription start points (TSP) of Nnat, identified through the presence of CAGED transcripts in mice and human. For example, mouse Nnat has six transcripts, with two promoters having respectively 36 and 28 CAGE tags, while human NNAT has four transcripts with one promoter having 29 CAGE tags. For these analyses we used the Genome VISTA genome alignment tool http://www-gsd. lbl.gov/vista/ in conjunction with Genomatix Gene2promoter tool (Genomatix Software Inc., Ann Arbor, MI, USA). The Genomatix MatInspector programme with default parameter settings was used to detect common Transcription Factor Binding Sites (TFBS) between the seven promoter sequences of the seven genera. The occurrence of the TFBSs within all vertebrate genes was used as a control in MatInspector. Significant common TFBSs among the seven species were defined by having a p-value < 0.05 and presence in six or more species, compared with the control sets in the Genomatix database. We further explored the TF genes in each TFBSs family, examining expression data to determine whether they were expressed in human or mouse adipose, hypothalamus and pancreas islet (using BioGPS database http://biogps.gnf.org and our own BAIR data as reference sources). Promoter analysis was also applied to these 27 fully conserved Nnat co-expression genes (plus Nnat) all together, focusing on human, mouse and rat only. Again, their expression status in adipose, hypothalamus and pancreatic islets were included as well.

Gene set enrichment analysis
In order to inspect which biological pathways are significantly correlated with Nnat expression in adipose tissue, we employed Gene Set Enrichment Analysis, using the GSEA package http://www.broadinstitute.org/gsea/ for each of the above four datasets separately. We chose Pearson correlation and permuted 1000 times to detect statistical significance. Expression data were pre-filtered by using the present/absent call and background correction, similarly to the pre-process in meta-analysis above. In the GSEA database there were 200 KEGG pathway sets and 639 canonical pathway sets. We first focused on those approximately 100 KEGG gene sets which have more than 15 genes and fewer than 500 genes; the results returned were significant in the above four datasets as well as an additional dataset which included 45 sleep mice model with hypothalamus tissue (GEO: GSE6514). We then worked on those canonical pathway sets which passed the same filter only on BAIR fat-fed mice and human fat dataset about 20 lean and 19 obese person [22].

Protein-protein interaction network analysis
By using BioNetBuilder, a cytoscape plugin package http://err.bio.nyu.edu/cytoscape/bionetbuilder/ with supported databases DIP, BIND, Prolinks, KEGG, HPRD, The BioGrid and GO, we thus generated a mouse protein-protein interaction network which had 11777 unique nodes (genes) and 107813 unique edges (links). For the 27 fully conserved Nnat co-expression genes, we searched each of their direct correlation genes (first neighbours) in the PPI network, again filtered by expression value > 200 in adipose tissue by BioGPS.

Nnat promoter analysis
We sought to investigate mechanisms of transcriptional regulation of Nnat and of co-varying genes in adipose tissue. The neuronatin gene is located on human chromosome 20 and mouse chromosome 2. To define the Nnat promoter we described the conserved non-coding region in seven mammalian species, including mouse and human, corresponding to 1200 bp upstream and 800 bp downstream of the most utilized transcription start points (TSP) of Nnat identified through the presence of CAGED transcripts in mice and humans. We used Genomatix software to determine the presence of binding sites for TF families in the conserved Nnat promoter, across seven mammalian species. We then filtered these potential regulators based on their expression in mouse WAT, and compared their expression also in hypothalamus and pancreatic b-cells in which Nnat is known to be expressed. Table 1 shows the gene expression behaviour for members of TFBS families scoring as significant for presence in the Nnat promoter (p-value < 0.05), and present in at least six of the Nnat promoter sequences; expression values were filtered on the basis of BioGPS expression bigger than 200, which therefore focuses attention on the most abundantly expressed TFs.
The TFBS for the V$NEUR family (including NeuroD1 and NeuroD2) was found to be highly significant in six out of seven promoter sequences, but it was not found in the rat promoter sequences in our defined region. However, V$NEUR does exist farther upstream from the rat Nnat transcription start site (within 2.5kb). NeuroD1, the only confirmed transcriptional regulator of Nnat expression [14] is not one of the most highly expressed TFs in adipose, though its expression value is similar across all three tissues shown in Table 1, it is highly expressed in cerebellum.
The binding sites for the key adipose tissue regulator Pparg are the most significant in the seven Nnat promoter sequences. Of the TFs with conserved binding sites it is also by far the most highly expressed in WAT, and has a crucial role in adipocyte metabolism. Other important TFs with conserved binding sites in the Nnat promoter and which are highly expressed in WAT include several genes which have crucial roles in the control of: adipogenesis (Klf15, Irf1, Creb1, Egr2, Gata3); lipogenesis (Mlxipl, Srebp1c); inflammation (Jun, Stat3); insulin signalling and diabetes susceptibility (Foxo1, Tcf7l2). Interestingly, the Rreb1 has a binding site in all seven Nnat promoters, and is expressed at high levels in adipose and hypothalamus. Rreb1 can enhance the transcriptionactivating function of NeuroD1, a property that requires both the physical interaction of Rreb1 with NeuroD1, and binding of Rreb1 to DNA [23]. Notable by their absence are binding sites for the C/EBP family of TFs, clearly suggesting that this family (several members of which are required for adipocyte differentiation) are not concerned with Nnat regulation, and may be down-stream of Nnat.

Identification of biological pathways associated with Nnat expression
In order to determine which biological pathways are significantly correlated with Nnat expression in WAT and hypothalamus, we employed Gene Set Enrichment Analysis. Initially we focused on KEGG pathways to investigate those processes that are most highly associated with Nnat expression in WAT and hypothalamus. In WAT, there are five KEGG pathways significantly correlated with Nnat expression (FDR < 0.05) in at least two of the four mouse WAT datasets (Table 2), these five pathways are all significant in BAIR fat-fed data, and include oxidative phosphorylation, ribosome, proteasome, cholera infection and neuroactive ligand receptor interaction, these last two indicating perturbation of ion transport and neuronal function. Furthermore, oxidative phosphorylation, ribosome and cholera infection are significant in the mouse hypothalamus dataset as well. Additional file 1, Figure S1 illustrates the enrichment plots and positively and negatively correlated genes in these three common pathways. Additional file 1, Table S1 lists the core enriched genes that correlate with Nnat in these three common pathways in both WAT and hypothalamus. GSEA on the canonical pathways (Additional file 1, Table S2) reveals that Nnat expression in mice on a high-fat diet and obese models displayed highly positive correlations with inflammatory pathways (complement and coagulation, cytokine and integrin signalling), and strongly negative correlations with metabolic pathways (most prominently oxidative phosphorylation, pyruvate metabolism, glycolysis and gluconeogenesis) and protein turnover (including proteasome, ribosome and mRNA processing). In human adipose tissue NNAT mRNA expression is decreased (Additional file 1, Table S3). In consequence NNAT expression displayed strong positive correlation with genes involved in protein turnover and fatty acid metabolism/degradation, and negatively correlated with genes involved in adhesion, complement and coagulation and cytokine. Additional file 1, Table S4 lists the common pathways correlated with Nnat expression between BAIR fat-fed mice and human obese WAT, and these are concerned with cytokine signalling and cell-cell interaction. Additional file 2, Table S5 contains the complete list of genes associated with Nnat in each KEGG pathway in mouse models.

Identification of genes that co-vary with Nnat
The genes whose expression co-varies with Nnat are likely to share regulatory mechanisms and, perhaps, biological functions. We therefore identified genes whose expression is highly correlated with Nnat in a meta-analysis of WAT gene expression data. Table 3 shows the correlation of gene expression behaviour with that of Nnat, in four microarray datasets; entries with a null value indicate that this gene is missing in the specific dataset; their fold change information in BAIR fat-fed mice and human with obesity are shown in Table 4. Additional file 1, Figure S2 gives the hierarchical clustering for these 28 genes in BAIR fat-fed dataset, here genes with expression most closely associated with Nnat include Aqp1, Sncg, Sulf2 and Cxcl9. Additional file 1, Figure S3 illustrates the expressions of the top four most correlated genes (Gstt1, Ccdc80, Hfe Sod3) with Nnat in BAIR and GSE6571 datasets. Among these genes is only one transcription factor, Ebf1 which has previously been suggested to have a regulatory role in adipogenesis [24]. Another of the Nnat-correlated genes, Lgals12, is a member of the galectin family of betagalactoside-binding proteins, and is a major regulator of adipose tissue development [25]. In order to investigate whether these 27 genes share conserved TFBS with Nnat as identified in Table 1, we used Genomatix tools to analyze the promoters of all 28 genes in human, mouse and rat separately. We found that nine TFBS families in Table 1, V$CREB, V$EVI1, V $GATA, V$KLFS, V$NR2F, V$RXRF, V$SORY, V$STAT, V$HAND, which are conserved in the promoter region of Nnat from seven species, are significant TFBS of all these 28 genes in human, mouse and rat, indicating these TFBSs are likely to be important in WAT. We also found that V$RREB and V$TEAF are highly significant TFBSs (p-value < 3.31e-10, p-value < 6.22e-8, p-value < 5.58e-7) in most of the 28 genes, while V$PPAR is highly significant (p-value < 5.7e-4, p-value < 4.14e-6, p-value < 5.7e-4) half in of the 28 genes (Table 5).
We used the Genomatix Promoter Module Library and were able to identify highly conserved models of five transcription factor elements highly conserved across the Nnat promoter. Among these promoter sequences, we got 20 significant models which consist of elements from five transcription factor family. Figure  1 shows TFBPs for V$EVI1 -V$EGRF -V$E2FF -V $NRF1 -V$LEFF, on the Nnat promoter close together in a similar region for six species separately just upstream of the major TSP. V$EVI1 -V$EGRF -V $E2FF and -V$LEFF are in Table 1 and TF members of these families are expressed in WAT, hypothalamus and pancreas. V$EVI1 -V$EGRF and -V$LEFF in Table 5 indicate conservation in the promoters of genes that co-vary with Nnat. There are important transcription factors in these families, for example Egr2 and Tcf7l2, which are implicated in stress responses and T2D; E2f which is concerned with adipogenesis and Evi1, which activates Stat3, a transcription factor concerned with leptin signalling in the hypothalamus.

Discussion
This study explored the behaviour of the diet-responsive gene Nnat in multiple gene expression datasets from mouse adipose tissue, and subsequently analysed in silico the promoter of Nnat in seven mammalian species. Rosen and MacDougald [24] have pointed out that at least 100 TFs are expressed in adipocyte, and our analyses here have sought to focus attention on those which could play a part in the regulation of expression of Nnat and its co-varying genes. Based on our study we find 50 TFBSs conserved across seven species on the Nnat promoter, and expression of the TFs corresponding to these in adipose tissue, many of which have known roles in adipocyte processes including, adipogenesis (Klf15, Irf1, Creb1, Egr2, Gata3); lipogenesis (Mlxipl, Srebp1c); inflammation (Jun, Stat3); insulin signalling and diabetes susceptibility (Foxo1, Tcf7l2). These TFs are also well expressed in the hypothalamus and pancreas. NeuroD1 is the only confirmed transcription factor for Nnat, whose expression regulates in the developing pancreas [14]. Binding sites for the TF family V$NEUR (which includes NeuroD1 and NeuroD2) are predicted in six of the seven mammalian Nnat proximal promoter regions examined here (all except rat), and are also present further upstream in the rat Nnat promoter region. Members of this family are expressed at reasonable levels in all three mouse tissues analysed here, so it remains possible that the V$NEUR family participates in Nnat regulation in several tissues. NeuroD1 is considered to be a weak activator of transcription, but is known to cooperate with other transcription factors, including Rreb1 and Sp1 [23,26], in the stimulation of gene expression. Binding sites for both of these families are conserved in the Nnat promoter across mammalian species, as well as in the promoters for the Nnat co-varying genes analysed above. Indeed this property requires both the physical interaction of Rreb1 with NeuroD1, and binding of Rreb1 to DNA. Direct evidence for such a mechanism of regulation requires experimental validation. The expression of Nnat remains strong in the developing pituitary of NeuroD1 knockout mice [14], so it is clear that different TFs (or combinations thereof) are likely to be of varying importance in regulating Nnat expression in different tissues. The common TFBS among the 28 (Nnat and its co-vary) genes in three species respectively. The last three columns give the number of promoter sequences among the 28 which share the same TFBS in their promoter sequences. Pparg is both necessary and sufficient for adipocyte differentiation of fibroblast/pre-adipocyte cell lines [24] and the presence, in all the Nnat promoters analysed here, of Pparg binding sites is consistent with it playing a role in regulation of Nnat expression in WAT. Nnat expression is strongly reduced in WAT of a mouse model in which a dominant negative version of Pparg is expressed [17]. Activation of Pparg with Pparg agonists improves insulin signalling and decreases inflammation. This is consistent with the view that Nnat protects against toxins and is concerned with resolution of tissue injury. Indeed our unpublished results in 3T3L1 cells, where we have over expressed Nnat indicate an anti-inflammatory function. However there is currently no evidence on whether Pparg plays a direct role in the regulation of Nnat. Likewise, in mice expressing an inactive mutant version of Srebp1c, Nnat expression in WAT is severely reduced as well [17]. In both of these mouse models complex changes in WAT gene expression make it difficult to draw conclusions as to the regulation of particular genes by specific TFs. Pparg is known to operate in a cellular network with other TFs of the EBF, C/EBP, KLF and IRF families in regulating adipogenesis [24,[27][28][29]. Our bioinformatical analysis suggests that Pparg could be a candidate for Nnat regulation, meanwhile the absence of V $SREB binding sites in the Nnat promoter does not support a direct role for Srebp1c.
The pattern of gene expression observed in the mouse WAT upon high fat feeding in experiments analysed here suggests a critical in oxidative stress; indeed sentinel indicators of this Gstt1, Ccdc80, Hfe and Sod3 all concerned with cellular management of free radicals are the genes most highly correlated with Nnat in the meta-analysis (Additional file 1, Figure S3). Genes concerned with adipogenesis (Ebf1, Fzd4 and Lgals12) and inflammation (Sncg, Cxcl9 and Aoc3) are also highly correlated in mouse WAT upon high fat feeding (Additional file 1, Figure S2). This is consistent with data indicating that Nnat expression favours adipogenesis and inflammation. Of the aforementioned adipogenic TFs, only Ebf1 expression correlates well with Nnat, but the conserved Nnat promoter does not possess Ebf1 binding sites (TF family V$NOLF). Thus it is unlikely that Ebf1 regulates Nnat expression; instead they may both be targets of the same TFs. By treating the Nnat co-varying genes as a set, and carrying out promoter analysis on this set, we discovered that many wide-acting TFs have binding sites in all of these genes. This does not itself indicate which, if any, of these TFs may regulate the Nnat co-varying genes, but does draw attention to the families V$HIFF (example: HIF1a) V$RREB (example: Rreb1) and V$TEAF (example: Tead4) as being the most highly over-represented TFBS in this group of genes.
Our gene set enrichment analysis showed that expression of Nnat is strongly negatively correlated with energy metabolism (most notably oxidative phosphorylation) and positively correlated with inflammation. Repression of oxidative phosphorylation gene expression is one of the major changes seen in human diabetes [21], and a metabolic/inflammatory network is thought to underlie metabolic syndrome in susceptible mice [30]. Although in the human data NNAT expression is not correlated with oxidative phosphorylation in the GSEA in a more general analysis of gene expression in this human data the oxidative phosphorylation gene set is enriched. We would expect that major changes in gene set regulation require the participation of wide-acting TFs and/or extensive networks of TFs and their cofactors.
We found 27 genes that show strong co-variation with Nnat. 21 of these can be mapped to the PPI network, and eight of these (Anxa6, Cav1, Cd151, Gstt1, Map1lc3a, Npdc1, Npr3, and Sncg) are expressed in adipose tissue, pancreas and hypothalamus. App and Tcf7l2, genes implicated in metabolic inflammation, Alzheimer's disease and diabetes are shown linked to their respective first neighbours (Figure 2). Sncg along with Nnat is a leptin responsive gene [11].

Conclusions
In conclusion we identify numerous TFBSs, conserved among mammals, on the Nnat promoter, and expression of TFs corresponding to these in adipose tissue, hypothalamus and pancreas. Many of these have known roles in adipocyte, as well as pancreatic and hypothalamic processes. We observe clustering of TFBSs for Egr1 and Tcf7l2, which are implicated in stress responses and T2D; E2f which is concerned with adipogenesis and Evi1, a developmental gene highly expressed in Figure 2 Mouse Nnat co-expression gene tissue-specific PPI network. All Nnat co-expression genes in Table 3 and their first neighbours in mouse PPI network showing expressed in WAT, hypothalamus and pancreas, generated a tissue-specific network with 151 nodes and 447 links. The eight genes in yellow are Nnat co-varying genes in Table 3. Two obesity relative genes App and Tcf7l2 and their direct associated genes linked by red edges. pancreatic islets. GSEA analysis suggests the roles for Nnat in oxidative stress, resolution of inflammation, energy metabolism, mainly oxidative phosphorylation, and protein turnover. Genes that co-vary with Nnat, such as Sncg, Gstt1, Hfe and Ebf1, share some of these TFBSs and form functional networks. The recent demonstration that SNPs in Nnat are found in obese humans suggest that altered Nnat function in adipose tissue or the hypothalamus may be involved in the pathogenesis of obesity via such mechanisms.