- Open Access
Molecular processes during fat cell development revealed by gene expression profiling and functional annotation
Genome Biology volume 6, Article number: R108 (2005)
Large-scale transcription profiling of cell models and model organisms can identify novel molecular components involved in fat cell development. Detailed characterization of the sequences of identified gene products has not been done and global mechanisms have not been investigated. We evaluated the extent to which molecular processes can be revealed by expression profiling and functional annotation of genes that are differentially expressed during fat cell development.
Mouse microarrays with more than 27,000 elements were developed, and transcriptional profiles of 3T3-L1 cells (pre-adipocyte cells) were monitored during differentiation. In total, 780 differentially expressed expressed sequence tags (ESTs) were subjected to in-depth bioinformatics analyses. The analysis of 3'-untranslated region sequences from 395 ESTs showed that 71% of the differentially expressed genes could be regulated by microRNAs. A molecular atlas of fat cell development was then constructed by de novo functional annotation on a sequence segment/domain-wise basis of 659 protein sequences, and subsequent mapping onto known pathways, possible cellular roles, and subcellular localizations. Key enzymes in 27 out of 36 investigated metabolic pathways were regulated at the transcriptional level, typically at the rate-limiting steps in these pathways. Also, coexpressed genes rarely shared consensus transcription-factor binding sites, and were typically not clustered in adjacent chromosomal regions, but were instead widely dispersed throughout the genome.
Large-scale transcription profiling in conjunction with sophisticated bioinformatics analyses can provide not only a list of novel players in a particular setting but also a global view on biological processes and molecular networks.
Obesity, the excess deposition of adipose tissue, is among the most pressing health problems both in the Western world and in developing countries. Growth of adipose tissue is the result of the development of new fat cells from precursor cells. This process of fat cell development, known as adipogenesis, leads to the accumulation of lipids and an increase in the number and size of fat cells. Adipogenesis has been extensively studied in vitro for more than 30 years using the 3T3-L1 preadipocyte cell line as a model. This cell line was derived from disaggregated mouse embryos and selected based on the propensity of these cells to differentiate into adipocytes in culture . When exposed to the appropriate adipogenic cocktail containing dexamethasone, isobutylmethylxanthine, insulin, and fetal bovine serum, 3T3-L1 preadipocytes differentiate into adipocytes .
Experimental studies on adipogenesis have revealed many important molecular mechanisms. For example, two of the CCAAT/enhancer binding proteins (C/EBPs; specifically C/EBPβ and C/EBPδ) are induced in the early phase of differentiation. These factors mediate transcriptional activity of C/EBPα and peroxisome proliferator-activated receptor (PPAR)-gamma (PPARγ) [3, 4]. Another factor, the basic helix-loop-helix (bHLH) transcription factor adipocyte determination and differentiation dependent factor 1/sterol regulatory element binding protein 1 (ADD1/SREBP1c), could potentially be involved in a mechanism that links lipogenesis and adipogenesis. ADD1/SREBP1c can activate a broad program of genes that are involved in fatty acid and triglyceride metabolism in both fat and liver, and can also accelerate adipogenesis . Activation of the adipogenesis process by ADD1/SREBP1c could be effected via direct activation of PPARγ  or through generation of endogenous ligands for PPARγ .
Knowledge of the transcriptional network is far from complete. In order to identify new components involved in fat cell development, several studies using microarrays have been initiated. These studies have used early Affymetrix technology [8–14] or filters , and might have missed many genes that are important to the development of a fat cell. The problem of achieving broad coverage of the developmental transcriptome became evident in a mouse embryo expressed sequence tag (EST) project, which revealed that a significant fraction of the genes are not represented in the collections of genes previously available . Moreover, earlier studies on adipogenesis [8–14] focused on gene discovery for further functional analyses and did not address global mechanisms.
We conducted the present study to evaluate the extent to which molecular processes underlying fat cell development can be revealed by expression profiling. To this end, we used a recently developed cDNA microarray with 27,648 ESTs , of which 15,000 are developmental ESTs representing 78% novel and 22% known genes . We then assayed expression profiles from 3T3-L1 cells during differentiation using biological and technical replicates. Finally, we performed comprehensive bioinformatics analyses, including de novo functional annotation and curation of the generated data within the context of biological pathways. Using these methods we were able to develop a molecular atlas of fat cell development. We demonstrate the power of the atlas by highlighting selected genes and molecular processes. With this comprehensive approach, we show that key loci of transcriptional regulation are often enzymes that control the rate-limiting steps of metabolic pathways, and that coexpressed genes often do not share consensus promoter sequences or adjacent locations on the chromosome.
Expression profiles during adipocyte differentiation
The 3T3-L1 cell line treated with a differentiation cocktail was used as a model to study gene expression profiles during adipogenesis. Three independent time series differentiation experiments were performed. RNA was isolated at the preconfluent stage (reference) and at eight time points after confluence (0, 6, 12, 24, 48 and 72 hours, and 7 and 14 days). Gene expression levels relative to the preconfluent state were determined using custom-designed microarrays with spotted polymerase chain reaction (PCR) products. The microarray developed here contains 27,648 spots with mouse cDNA clones representing 16,016 different genes (UniGene clusters). These include 15,000 developmental clones (the NIA cDNA clone set from the US National Institute of Aging of the National Institutes of Health NIH), 11,000 clones from different brain regions in the mouse (Brain Molecular Anatomy Project [BMAP]), and 627 clones for genes which were selected using the TIGR Mouse Gene Index, Build 5.0 .
All hybridizations were repeated with reversed dye assignment. The data were filtered, normalized, and averaged over biological replicates. Data processing and normalization are described in detail under Materials and methods (see below). Signals at all time points could be detected from 14,368 elements. From these microarray data, we identified 5205 ESTs that exhibited significant differential expression between time points and had a complete profile (P < 0.05, one-way analysis of variance [ANOVA]). Because ANOVA filters out ESTs with flat expression profiles, we used a fold change criteria to select the ESTs for further analysis. We focused on 780 ESTs that had a complete profile over all time points, and that were more than twofold upregulated or downregulated in at least four of those time points. These stringent criteria were necessary to select a subset of the ESTs for in-depth sequence analysis and for examination of the dynamics of the molecular processes. The overlap between the ANOVA and twofold filtered ESTs was 414. All of the data, together with annotations and other files used in the analyses, are available as Additional data files and on our website . The analyses described in the following text were conducted in the set of 780 ESTs.
Validation of expression data
Four lines of evidence support the quality of our data and its consistency with existing knowledge of fat cell biology. First, our array data are consistent with reverse transcriptase (RT)-PCR analysis. We compared the microarray data with quantitative RT-PCR for six different genes (Pparg [number 592, cluster 6], Lpl [number 14, cluster 6], Myc [number 224, cluster 11], Dcn [number 137, cluster 7], Ccna2 [number 26, cluster 5/8], and Klf9 [number 6, cluster 9]) at different time points (Additional data file 9 and on our website ). A high degree of correlation was found (r2 = 0.87), confirming the validity of the microarray data.
Second, statistical analyses of the independent experiments showed that the reproducibility of the generated data is very high. The Pearson correlation coefficient between the replicates was between 0.73 and 0.97 at different time points. The mean coefficient of variation across all genes at each time point was between 0.11 and 0.27. The row data and the details of the statistical analyses can be found in Additional data file 10 and on our website .
Third, comparison between our data and the Gene Atlas V2 mouse data for adipose tissue  shows that the consistency of the two data sets increases with differentiation state (Additional data files 11 and 12, and our on website ). Therefore, this analysis supports the relevance of the chosen cell model to in vivo adipogenesis. Among the 382 transcriptionally modulated genes common in both data sets, 67% are regulated in the same direction at time point zero (confluent pre-adipocyte cell culture). At the final stage of differentiation, the correlation increases up to 72%. If the Gene Atlas expression data are restricted to strongly regulated genes (at least twofold and fourfold change respectively), then the consistency in mature adipocytes rises to 82% (135 genes) and 93% (42 genes), respectively. Out of all 60 tissues in the Gene Atlas V2 mouse, the adipose tissue describes the differentiated state of the 3T3-L1 cells best. Brown fat tissue is the second best match to the differentiated adipocytes (69% of the 382 genes), followed by adrenal gland (66%), kidney (65%), and heart (64%). At each time point in which cell cycle genes were not repressed (12 hours and 24 hours), all tissues had similar correlation to the data set (44-55% for 382 genes).
The fourth line of evidence supporting the quality of our data is that there is clear correspondence between our data and a previously published data set . For a group of 153 genes shared among the two studies, the same upregulation or downregulation was found for 72-89% (depending on time point) of all genes (see Additional data file 13 and our website ). The highest identity (89%) was found for the stage terminally differentiated 3T3-L1 cells, for which the profile is less dependent on the precise extraction time. If the comparison is restricted to expression values that are strongly regulated in both experiments (at least twofold change at day 14, 96 ESTs), then the coincidence at every time point is greater than 90%. Comparisons with this  and two additional data sets [9, 12], and the data pre-processing steps are given in Additional data files 13, 14, 15 and on our website . Note that, because of the differences in the used microarray platforms, availability of the data, normalization methods, and annotations, only 96 genes are shared between all four studies. Of the 780 ESTs monitored in our work, 326 were not detected in the previous studies [8, 9, 12]: 106 RefSeqs (with prefix NM), 43 automatically generated RefSeqs (with prefix XM), and 173 ESTs (Additional data file 16).
Correspondence between transcriptional coexpression and gene function
To examine the relationships between coexpression and gene functions, we first clustered 780 ESTs that were twofold differentially expressed into 12 temporally distinct patterns, containing between 23 and 143 ESTs (Figure 1). ESTs in four of the clusters are mostly upregulated during adipogenesis, whereas genes in the other eight clusters are mostly downregulated.
We then categorized ESTs with available RefSeq annotation and Gene Ontology (GO) term (486 out of 780) for molecular function, cellular component, and biological process (Figure 2). Genes in clusters 5 and 8 are downregulated through the whole differentiation process and upregulated at 12/24 hours. Many of the proteins encoded by these genes are involved in cell cycle processes and were residing in the nucleus (Figure 2). Re-entry into the cell cycle of growth arrested pre-adipocytes is known as the clonal expansion phase and considered to be a prerequisite for terminal differentiation in 3T3-L1 adipocytes . Genes grouped in cluster 2 are highly expressed from 6 hours (onset of clonal expansion) to 3 days (start of the appearance of adipocyte morphology) but are only modestly expressed at the terminal adipocyte differentiation stage. These include a number of genes that encode signaling molecules. Genes increasingly expressed toward the terminal differentiation stage are in clusters 4, 6, and 7, although from different starting values. Some genes in cluster 6 are known players in lipid metabolism and mitochondrial fatty acid metabolism, whereas some genes can be associated with cholesterol biosynthesis and related to extracellular space or matrix in clusters 4 and 7, respectively.
Correspondence between coexpression and targeting by microRNAs
Previous studies suggest that protein production for 10% or more of all human and mouse genes are regulated by microRNAs (miRNAs) [23, 24]. miRNAs are short, noncoding, single-strand RNA species that are found in a wide variety of organisms. miRNAs cause the translational repression or cleavage of target messages . Some miRNAs may behave like small interfering RNAs. It appears that the extent of base pairing between the small RNA and the mRNA determines the balance between cleavage and degradation . Rules for matches between miRNA and target messages have been deduced from a range of experiments  and applied to the prediction and discovery of mammalian miRNA targets [23, 27]. Moreover, it was shown that human miRNA-143 is involved in adipocyte differentiation .
Here we conducted an analysis to determine which of the 780 ESTs differentially expressed during adipocyte differentiation were potential targets of miRNAs and whether there is an over-representation of miRNA targets of coexpressed ESTs clustered in 12 distinct expression patterns. From the 780 ESTs, the 3'-untranslated region (UTR) could be derived for 539. Of these, 518 had at least one exact antisense match for the seven-nucleotide miRNA seed (base 2-8 at the 5' end) from the 234 miRNA sequences (18-24 base pairs [bp]; Additional data file 14). From 395 ESTs with a unique 3'-UTR, 282 (71%) had at least one match over-represented compared with the whole 3'-UTR sequence set (21,396; P < 0.05, by one-sided Fisher's exact test). The distribution of statistically over-represented miRNA motifs in 3'-UTRs across the clusters was variable, with genes grouped in cluster 9 (including many transcriptional regulators) having the most statistically over-represented miRNA motifs and genes in cluster 5 having no detectable motifs (Additional data file 18). The results of the analysis of cluster 9 are given in Figure 3. One of the genes with the most significantly over-represented miRNA motifs in the 3'-UTR is related to the ras family (Figure 3). It was previously shown that human oncogene RAS is regulated by let-7 miRNA . Further potential miRNA target genes from all clusters are given in Additional data files 9, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30.
Molecular atlas of fat cell development derived by de novofunctional annotation of differentially expressed ESTs
In order to functionally characterize the molecular components underlying adipogenesis in detail, comprehensive bioinformatics analyses of 780 differentially expressed ESTs were performed. A total of 659 protein sequences could be derived, and these were subjected to in-depth sequence analytic procedures. The protein sequences have been annotated de novo using 40 academic prediction tools integrated in the ANNOTATOR sequence analysis system. The structure and function was annotated on a sequence segment/domain-wise basis. After extensive literature search and curation using the sequence architecture, 345 gene products were mapped onto known pathways, possible cellular roles, and subcellular localizations (Figure 4) using the PathwayExplorer web service  as well as manual literature and domain-based assignment. The results of the sequence analyses and additional information is available in the supplementary material available on our website  and Additional data files 6, 7, 8.
This molecular atlas of fat cell development provides the first global view of the underlying biomolecular networks and represents a unique resource for deriving testable hypotheses for future studies on individual genes. Below we demonstrate the usefulness of the atlas by highlighting the following: established regulators of fat cell development, recently discovered fat cell gene products, and candidate transcription factors expressed during adipogenesis. The numbering of the genes is given according to the de novo functional annotation (Additional data file 7).
Established regulators of fat cell development
Key transcription factors SREBF1 (Srebf1 [number 119, cluster 9]) and PPARγ (Pparg [number 592, cluster 6]) were highly expressed during the late phase of differentiation. PPARγ  (Pparg [number 592, cluster 6]) is increasing up to about 15-fold. Srebf1 processing is inhibited by insulin-induced gene 1 (Insig1 [number 62, cluster 3/4]) through binding of the SREBP cleavage-activation protein [32, 33]. Insig1 is regulated by Srebf1 and Pparg at the transcriptional level  and the expression of known marker genes of the differentiated adipocyte was increased in parallel with these factors. These include genes from clusters 3, 6, and 9 that are targets of either of these factors: lipoprotein lipase (Lpl [number 14, cluster 6]), c-Cbl-associated protein (Sorbs1 [number 92, cluster 6]), stearoyl-CoA desaturase 1 (Scd1 [number 305, cluster 6]), carnitine palmitoyltransferase II (Cpt2 [number 43, cluster 9]), and acyl-CoA dehydrogenase (Acadm [number 153, clusters 6 and 9]).
Recently discovered fat cell gene products
During the preparation of the manuscript, a number of factors shown to be important to adipocyte function were identified in vivo. All of these factors, which have a possible role in the pathogenesis of obesity and insulin resistance, were highly expressed in the present study. Adipose triglyceride lipase (Pnpla2 [number 157, cluster 6]), a patatin domain-containing triglyceride lipase that catalyzes the initial step in triglyceride hydrolysis , was more than 20-fold upregulated at the terminal differentiation phase. Another example is Visfatin, which is identical to the pre-B cell colony-enhancing factor (Pbef [number 327, cluster 9]). This 52 kDa cytokine has enzymatic function in adipocytes, exerts insulin-mimetic effects in cultured cells, and lowers plasma glucose levels in mice by binding to the insulin receptor [36–38]. The imprinted gene mesoderm-specific transcript (Mest [number 17, cluster 6/9]), which appears to enlarge adipocytes and could be a novel marker of the size of adipocytes , is upregulated during the late stage of 3T3-L1 differentiation.
Members of the Krüppel-like factor (Klf) family, also known as basic transcription element binding proteins, are relevant within the context of adipocyte differentiation. Klf2 was shown to inhibit PPARγ expression and to be a negative regulator of adipocyte differentiation ; Klf5 , Klf6 , and Klf15  have been demonstrated to induce adipocyte differentiation. Whereas Klf9 (Bteb1 [number 6, cluster 9]) was upregulated in the intermediate phase in the present study, Klf4 (number 100, cluster 12), which was shown to exert effects on cell proliferation opposing those of Klf5 , was downregulated. Another twofold upregulated player is Forkhead box O1 (Foxo1 [number 53, cluster 9]), which mediates effects of insulin on the cell. Activation occurs before the onset of terminal differentiation, when Foxo1 becomes dephosphorylated and localizes to the nucleus [44, 45]. The glucocorticoid-induced leucine zipper (Tsc22d3/Gilz [number 173, cluster 2]) functions as a transcriptional repressor of PPARγ and can antagonize glucocorticoid-induced adipogenesis [46, 47]. This is consistent with our observation that Gilz is highly upregulated during the first two days, when dexamethasone is present in the medium, and downregulated at the end of differentiation, when PPARγ is highly induced. C/EBP homologous protein 10 (Ddit3 [number 498, cluster 3]), another type of transcriptional repressor that forms nonfunctional heterodimers with members of the C/EBP family, was early induced and then downregulated. This might be sufficient to restore the transcriptional activity of C/EBPβ and C/EBPδ . The transcription factor insulinoma-associated 1 (Insm1 [number 238, cluster 8]) is associated with differentiation into insulin-positive cells and is expressed during embryo development, where it can bind the PPARγ target Cbl-associated protein (Sorbs1 [number 92, cluster 6]; upregulated after induction) [48, 49].
Candidate transcription factors expressed during adipogenesis
Because knowledge of the transcriptional network during adipogenesis is far from complete, expression profiles have been generated and screened for candidate transcription factors [8, 9, 12]. Here, we identified a number of transcription factors the exhibit distinct kinetic profiles during adipocyte differentiation that were previously not functionally associated with adipogenesis. Two transcription factors were unique to the present study (Zhx3 and Zfp367), and three more were confirmed (Zhx1, Twist1 and Tcf19) and annotated in the pathway context.
We found evidence for a role of the zinc finger and homeobox protein 3 (Zhx3 [number 306, cluster 2]). Zhx3 as well as Zinc finger and homeobox protein 1 (Zhx1 [number 386, cluster 9]) might attach to nuclear factor Y, which in turn binds many CCAAT and Y-box elements . We also provide data regarding the expression of zinc finger protein 367 (Zfp367 [number 320, cluster 8]) during adipogenesis. The molecular function of Zfp367 is as yet uncharacterized.
Additionally, we provide further experimental evidence and pathway context for candidate transcription factors previously identified in microarray screens [9, 12], namely Twist1 and Tcf19. The Twist gene homolog 1 (Twist1 [number 235, cluster 9]) was about two- to threefold upregulated at 0 hours, 72 hours, 7 days, and 14 days. Twist1 is a reversible inhibitor of muscle differentiation . Heterozygous double mutants (Twist1 -/+, Twist2 -/+) exhibit loss of subcutaneous adipose tissue and severe fat deficiency in internal organs . Twist1 is a downstream target of nuclear factor-κB and can repress transcription of tumor necrosis factor-α, which is a potent repressor of adipogenesis [52, 53]. The differential expression during adipogenesis of Tcf19 was also confirmed in the present study. Tcf19 is a transcription regulator that is involved in cell cycle processes at later stages in cell cycle progression . Expression of other regulators that are involved in the same process support this observation. Forkhead box M1 (Foxm1 [number 194, cluster 8]) stimulates the expression of cell cycle genes (for instance the genes encoding cyclin B1 and cyclin B2, and Cdc25B and Cdk1). In addition, TAF10 RNA polymerase II, also known as TATA box binding protein-associated factor (Taf10 [number 518, cluster 8]), is involved in G1/S progression and cyclin E expression .
Correspondence between phenotypic changes and gene expression
In addition to the metabolic networks, the molecular atlas also provides a bird's eye view of other molecular processes, including signaling, the cell cycle, remodeling of the extracellular matrix, and cytoskeletal changes. Changes that occur during adipogenesis (phenotypically seen as rounding of densely packed cells) have aspects in common with other tissue differentiation processes such as endothelial angiogenesis (protease, collagen, and noncollagen molecule secretion)  and specific features. Here we show that phenotypic changes that occur in maturing adipocytes are paralleled by expression of the respective genes.
Extracellular matrix remodeling
Matrix metalloproteinase-2 (MMP-2 [number 342, cluster 2]) was strongly upregulated during the entire process of adipocyte differentiation. Matrix metalloproteinase-2 can cleave various collagen structures and its inhibition can block adipogenesis . Tissue inhibitor of metalloproteinase-2 (Timp2 [number 239, cluster 9]), a known partner of matrix metalloproteinase-2, which balances the activity of the proprotease/protease , was mainly upregulated. Decreased levels of tissue inhibitor of metalloproteinase-3 (number 81, cluster 10; upregulated at 6 hours and repressed after 12 hours) are associated with obese mice . New collagen structures of overexpressed Col6a2 (number 11, cluster 9), Col4a1 (number 58, cluster 2) and Col4a2 (number 303, cluster 2)  are cross-linked by the lysyl oxidase (Lox [number 282, cluster 2]; upregulated during adipogenesis, which is contrary to findings reported by Dimaculangan and coworkers ). Strongly upregulated decorin (Dcn [number 137/623, cluster 7]) and osteoblast specific factor 2 (Postn/Osf-2 [number 183, cluster 7]), as well as proline arginine-rich end leucine-rich repeats (Prelp [number 73/484, cluster 3]; upregulated in the final stages of adipogenesis), attach the matrix to the cell. Matrillin-2 (Matn2 [number 12, cluster 9]; upregulated during adipogenesis) functions as adaptor for noncollagen structures , as does nidogen 2 (Nid2 [number 294, clusters 6 and 9]; increasingly upregulated). Secreted protein acidic and rich in cysteine/osteonectin (SPARC [number 67, cluster 9]; mainly upregulated) and SPARC-like 1 (Sparcl1 [number 154, cluster 9]; upregulated at 0 hours, 72 hours, 7days, and 14 days) can organize extracellular matrix remodeling, inhibit cell cycle progression, and induce cell rounding in cultured cells [63, 64].
Reorganization of the cytoskeleton
Most cytoskeletal proteins are coexpressed in cluster 10 (not repressed from 6 to 12 hours) and might have a common regulatory mechanism. Transcription of actin α (Acta1 [number 445, cluster 10]) and actin γ (Actg1 [number 656, cluster 10]), tubulin α (Tuba4 [number 377, cluster 8]), and tubulin β (Tubb5 [number 110, cluster 8]) were found to diminish during differentiation, which is in agreement with other reports . Myosin light chain 2 (Mylc2b/Mylpf [number 87/88/52/421, cluster 10]), and tropomyosin 1 and 2 (Tpm1/Tpm2 [number 74/68, cluster 10]) are members of the mainly repressed cluster 10. The downregulated transgelin 1 and 2 (Tagln/Tagln2 [number 114/242, cluster 10/8]) as well as fascin homolog 1 (Fscn1 [number 30, cluster 10]) are known actin-bundling proteins [66, 67]. Apparently, their absence decreases the cross-linking of microfilaments in compact parallel bundles. Calponin 2 (Cnn2 [number 7, cluster 10]), a regulator of cytokinesis, is downregulated . The insulin receptor and actin binding proteins filamin α and β (Flna/Flnb [number 506/632, cluster 10]) can selectively inhibit the mitogen-activated protein kinase signaling cascade of the insulin receptor . Finally, the maintenance protein ankycorbin (Rai14 [number 59, cluster 10]) and the cross-linking protein actinin 1 (Actn1 [number 521, cluster 10]) share the mainly repressed expression profile. Tubulin γ 1 (Tubg1 [number 78, cluster 7]; upregulated during adipogenesis, about 42-fold at 6 hours) is not a component of the microtubules like Tuba/Tubb, but it plays a role in organizing their assembly and in establishing cell polarity . Actinin 4 (Actn4 [number 185, cluster 9]; upregulated throughout adipogenesis) differs from Actn1 in its localization. Its expression leads to higher cell motility, and it can be translocated into the nucleus upon phosphatidylinositol 3-kinase inhibition . Adducin 3γ (Add3 [number 50, cluster 9]; permanently upregulated) has different actin-associated cytoskeletal roles.
T-lymphoma invasion and metastasis 1 (Tiam1 [number 159, cluster 2]) is a guanine nucleotide exchange factor of the small GTPase Rac1, which regulates actin cytoskeleton, morphology and adhesion, and antagonizes RhoA signaling [72, 73]. Additionally, the putative constitutive active Rho GTPase ras homolog gene family, member U/Wnt1 responsive Cdc42 homolog (Rhou/Wrch-1 [number 292, clusters 2 and 7]), which has no detectable intrinsic GTPase activity and very high nucleotide exchange capacity, leads to an phenotype of mature adipocyte [74, 75]. Interplay between Rhou and Tiam1, which might reverse Rhou activity through Rac1 signaling , could be a mechanism for regulating cell morphology in adipogenesis.
In summary, the evidence presented above suggests that reduced replenishment of the cytoskeleton with building blocks and the strong transcriptional upregulation of modulating proteins, together with the extracellular remodeling, are responsible for the morphological changes that occur during differentiation of 3T3-L1 cells.
Regulation of metabolic networks at the transcriptional level via key points of pathways
We next used the molecular atlas to derive novel biological insights from the global view of molecular processes. We analyzed transcriptionally regulated genes that are members of 36 different metabolic pathways. Within each pathway, we considered whether these transcriptionally regulated genes occupy key positions, such as a position at the pathway start, which is the typical rate-limiting step where the amount of enzyme is critical , or at some other point of regulation. We found that such key positions are occupied by transcriptionally regulated targets in 27 pathways (an overview is provided in Table 1). Those pathways that are strongly transcriptionally regulated at key points are illustrated in Figure 5 at the time points 0, 24 and 48 hours, and 14 days. For additional time points and images with more detailed information from all investigated pathways, see our website  and Additional data files Additional data file 31 and Additional data file 32.
In the following discussion we present the evidence for transcriptional regulation at key points for five selected metabolic pathways. Further information on other pathways can be found in Additional data files 31, 32, 33, 34, 35, 36, 37, 38 and on our website .
Biosynthesis of the important lipogenic cofactors CoA and NAD(P)+ are transcriptionally regulated at their key enzymes
Coenzyme A (CoA) is the carrier of the fatty acid precursor acetate/malonate [77, 78]. Panthotenate kinase 3 (Pank3 [number 140, cluster 6]; about eightfold upregulated) is responsible for the first and rate-limiting step in converting panthotenate to CoA . Nicotinamide adenine dinucleotide phosphate (reduced form; NADPH) is necessary in reductive reactions for fatty acid synthesis. Pre-B-cell colony-enhancing factor (Visfatin/Pbef1 [number 327, cluster 9]; strongest upregulated in the last three points of the time course in parallel with the emergence of fat droplets) is the rate-limiting enzyme in NAD(P)+ biosynthesis [38, 80]. For reduction of NADP+ to NADPH, two major mechanisms are responsible: the pentose phosphate shunt and the tricarboxylate transport system. Hexose-6P dehydrogenase (H6pd [number 533, cluster 9]; upregulated throughout adipogenesis) is the rate limiting enzyme of the pentose phosphate shunt in the endoplasmic reticulum and provides NADPH to its lumen . In the cytosolic pendant in the pentose phosphate shunt, the transaldolase (Taldo1 [number 160, cluster 3]) is repressed at early stages and is about threefold upregulated at the end of 3T3-L1 differentiation. This expression change appears to switch the shunt between ribose-5-phosphate (for nucleic acid synthesis) and NADPH (for fatty acid production) synthesis at early and late time points, respectively. A similar expression profile is observed for the cytosolic NADP-dependent malic enzyme (Mod1 [number 76, cluster 3]) and the citrate transporter (Slc25a1/Ctp1 [number 209, cluster 3]). Both are part of the tricarboxylate transport system through the mitochondrial membrane. Transcription of the anaplerotic pyruvate carboxylase (Pcx [number 149, cluster 6]; activated by acetyl-CoA) is increasingly upregulated up to 16-fold toward the final two time points.
Fatty acid modification and assimilation is transcriptionally regulated at the rate-limiting steps
The transcriptional expression of stearoyl-CoA desaturase 1 (Scd1 [number 305, cluster 6]), which catalyzes the rate-limiting reaction of monounsaturated fatty acid synthesis and which is an important marker gene of adipogenesis [82, 83], is downregulated at induction but increases up to 60-fold with advancing adipogenesis. In contrast to previous reports , we found that the gene for elongation of long-chain fatty acid (Elovl6 [number 162, cluster 12]) protein, which may be the rate-limiting enzyme of long chain elongation to stearate , is not overexpressed in differentiated 3T3-L1 cells as in adipose tissue. Elovl6 appears repressed during the entire process of adipogenesis in 3T3-L1 cells. Expression of lipoprotein lipase (Lpl [number 14, cluster 6]), the rate-limiting enzyme of extracellular triglyceride-rich lipoprotein hydrolyzation and triglyceride assimilation [85–87], increases with time up to 21-fold in differentiated adipocytes.
Transcriptional regulation of triglyceride and fatty acid degradation is performed at key points
Adipose triglyceride lipase (Pnpla2/Atgl [number 157, cluster 6]) executes the initial step in triglyceride metabolism . Its expression increases strongly with differentiation progression. Acyl-CoA dehydrogenases (Acadm/Acadsb [number 153/220, clusters 6 and 9]) , the rate-limiting enzymes of medium, short and branched chain β-oxidation, are strongly upregulated in the final four time points. In contrast, the acyl-CoA dehydrogenase (Acadvl) of very long chain fatty acids is not in the set of distinctly differentially regulated genes, and exhibits some upregulation at the final two time points. This difference in expression might shift the enrichment from short and medium to long chain fatty acids during adipogensis. Branched chain ketoacid dehydrogenase E1 (Bckdha [number 193, clusters 3 and 9]) is the rate-limiting enzyme of leucine, valine, and isoleucine catabolism and is known to be inhibited by phosphorylation . Its gene shares a similar expression profile with the Acad genes. The elevated degradation of amino acids allows conversion to fatty acids through acetyl-CoA.
Several important nucleotide biosynthetic pathway enzymes follow a cell cycle specific expression profile (strongly repressed except between 12 and 24 hours)
Phosphoribosylpyrophosphate amidotransferse (Ppat [number 287, cluster 11])  is rate-limiting for purin production. Deoxycytidine kinase (Dck [number 363, cluster 8]) is the rate-limiting enzyme of deoxycytidine (dC), deoxyguanosine (dG) and deoxyadenosine (dA) phosphorylation [91–93]. Ribonucleotide reductase M2 (Rrm2 [number 448, cluster 8]) converts ribonucleotides to desoxyribonucleotides [94, 95]. Additionally, thymidine kinase 1 (Tk1 [number 165, cluster 5]) and dihydrofolate reductase (Dhfr [number 161, cluster 5/8]) play important roles in dT and purin biosynthesis during the cell cycle. In contrast, purin degradation is about sixfold upregulated between 6 and 72 hours by the rate-limiting xanthine dehydrogenase (Xdh [number 361, cluster 2]) [96, 97]. These findings are in concordance with those of a previous study , which showed that mitotic clonal expansion is a prerequisite for differentiation of 3T3-L1 preadipocytes into adipocytes. After induction of differentiation, the growth-arrested cells synchronously re-enter the cell cycle and undergo mitotic clonal expansion, as monitored by changes in cellular DNA content . In accord with this experimental evidence, we observed changes in cell cycle genes, most of which were in clusters 5 and 8 (see our website  and Additional data file 37).
Cholesterol biosynthesis is regulated by expression of key steps and whole pathway segments
The synthesis of the early precursor molecule 3-hydroxy-3-methylglutaryl (HMG)-CoA, which might be also used in other metabolic pathways, is transcriptionally controlled at the key enzymes HMG-CoA synthase (Hmgcs1 [number 178, cluster 4]; repressed except in terminal stages) and HMG-CoA reductase (Hmgcr [number 619, cluster 12]; always repressed), which is the rate-limiting enzyme of the cholesterol and mevalonate pathway [98, 99]. After the step of isopentenylpyrophosphate synthesis, cholesterol biosynthesis genes are coexpressed in cluster 4.
Correspondence between coexpression and coregulation
To determine whether coexpressed genes are also coregulated, we analyzed the available promoter sequences of the 780 ESTs. Promoter sequences could be retrieved for 357 genes. Most ESTs are sequenced from the 3' end, and hence it is easier to retrieve the 3'-UTR. Retrieval of promoters is more difficult than retrieval of the 3'-UTR because of experimental problems in extracting full-length cDNAs (and hence transcription start sites) and insufficient computational methods for identifying beginning of the 5'-UTR. We analyzed the occurrences of the binding sites of all transcription factors in vertebrates from the TRANSFAC database. Based on statistical analyses, among transcription factors with binding site motifs described in TRANSFAC  those listed in Table 2 are the most promising candidates for further functional studies on transcriptional regulation.
One example of a functional transcription factor binding site is SREBP-1 in cluster 4. A comparison among clusters showed that cluster 4 has significantly more genes with a SREBP-1 (SRE and E-box motifs ) binding site than all other clusters (P = 0.0484, Fisher's exact test; Table 3). Similarly, a putative SREBP-1 regulatory region is significantly more frequent in the promoters of the genes in cluster 4 compared with all unique sequences in the PromoSer database (P < 0.0289; PromoSer contains 22,549 promoters of 12,493 unique sequences). For a subset of the genes in cluster 4 with predicted SREBP-1 binding site (most genes of the cholesterol biosynthesis pathway), transcriptional regulation with SREBP-1 has been experimentally proven .
Surprisingly, binding sites for the key regulators of adipogenesis, namely PPAR and C/EBP, are not significantly over-represented in any of the promoters of the coexpressed genes. We generated a novel matrix for PPAR using 22 experimentally verified binding sites from the literature and analyzed the promoters of the coexpressed genes and all PromoSer promoters. Again, using this matrix the PPAR binding sites were not significantly over-represented.
Genomic position of coexpressed genes
Finally, we considered whether coexpressed genes also colocalize on the chromosomes. In a broad genomic interval (5 megabases [Mb]) on each mouse chromosome we mapped the ESTs from each cluster. Unexpectedly, our data do not support the observation of the highly significant correlation in the expression and genomic positioning of the genes. A typical example of mapped ESTs to chromosome 10 is illustrated in Figure 6, showing that expression levels of colocalized ESTs are divergent because only two mapped ESTs are members of the same cluster.
Additionally, we analyzed the genomic position of 5,205 ESTs that exhibited significant differentially expression between time points (P < 0.05; one-way ANOVA). These ESTs were grouped in 12 clusters, and we then searched for regions with three or more members in a genomic interval of 500 kilobases (kb). On average, 7 ± 5% of the ESTs from one cluster were colocalized. Comprehensive results of this analysis are accessible within the supplementary website  and Additional data files 42, 43, 44, 45.
In summary, these data do not provide evidence that colocalized genes in the genomic sequence are subject to the same transcriptional regulation (coexpression), as indicated by examples for different processes in other studies .
The data presented here and the functional annotation considerably extend upon previous microarray analyses of gene expression in fat cells [8–14] and demonstrate the extent to which molecular processes can be revealed by global expression profiling in mammalian cells. Our strategy resulted in a molecular atlas of fat cell development and provided the first global view of the underlying biomolecular networks. The molecular atlas and the dissection of molecular processes suggest several important biological conclusions.
First, the data support the notion that there are hundreds of mouse genes involved in adipogenesis that were not previously linked to this process. Out of the 780 selected genes, 326 were not shared with previous studies [8, 9, 12], suggesting that our view of this process is far from complete. Using microarrays enriched with developmental ESTs, we were able not only to identify new components of the transcriptional network but also to map the gene products onto molecular pathways. The molecular atlas we developed is a unique resource for deriving testable hypothesis. For example, we have identified several differentially expressed genes, including recently discovered gene products (Pnpla2, Pbef1, Mest) and transcription factors not previously detected in microarray screens (Zhx3 and Zfp367).
Second, from our global analysis of the potential role of miRNAs in fat cell differentiation, we were able to predict potential target genes for miRNAs in 71% of the 395 genes with a unique 3'-UTR that were differentially expressed during adipocyte differentiation. The distribution of predicted miRNA targets indicated that one miRNA may regulate many genes and that one gene can be regulated by a number of miRNAs. The function of the potential target genes was diverse and included transcription factors, enzymes, transmembrane proteins, and signaling molecules. Genes with the lowest number of over-represented miRNA motifs were cell cycle genes (clusters 5 and 8), whereas genes grouped in cluster 9 exhibited the most over-represented miRNA motifs in relation to the matches in the control set of all available 3'-UTR sequences. Genes in cluster 9 exhibited high expression values at time point 0 and may include genes relevant to the transition from pre-confluent to confluent cells. Genes in cluster 9 also represent molecular components that are involved in other cell processes, including extracellular matrix remodeling, transport, metabolism, and fat cell development (for example, Foxo1 [44, 45]). Genes in other clusters exhibited varying percentages of over-represented miRNA motifs and can be associated with diverse biological processes (Additional data file 6). As an example of functional miRNA targets, we showed that one signaling molecule of the ras family is a potential target of miRNAs, which is consistent with a previous observation in humans, in whom it was shown that the human RAS oncogene is regulated by the let-7 miRNA. This example indicates that the present analysis provides promising candidates ranked according to their significance of over-representation and the number of different miRNAs that might regulate these targets in the specific context of adipocyte differentiation. It should be noted that our analysis included only known miRNAs, suggesting that the number of target sites can be even higher. This striking observation could have implications for post-transcriptional regulation of other developmental processes. Microarrays for the analysis of miRNA expression are becoming available and future studies will shed light on the role of miRNAs in the context of cell differentiation.
Third, we were also able to characterize the mechanisms and gene products involved in the phenotypic changes of pre-adipocytes into mature adipocytes. Although the number of selected genes in this study was limited, we characterized gene products for extracellular matrix remodeling and cytoskeletal changes during adipogenesis. Other molecular components involved in these processes can be identified by mapping the characterized gene products onto curated pathways  and selecting missing candidates for further focused studies. Notably, most of the cytoskeletal proteins are coexpressed in cluster 10 and might have a common regulatory mechanism. Further computational and experimental analyses are needed to verify this hypothesis.
In addition to new information about fat cell development, our comprehensive analysis has provided new general biological insights that could only be derived from such a global analysis. First, we were able to examine at what points metabolic pathways were regulated. The global view of biological processes and networks derived from expression profiles showed that the metabolic networks are transcriptionally regulated at key points, usually the rate-limiting steps. This was the case in 27 out of the 36 metabolic pathways analyzed in this study. During the development of mature adipocytes from pre-adipocytes, distinct metabolic pathways are activated and deactivated by this molecular control mechanism. For example, at the beginning nucleotide metabolism is activated because the cells undergo clonal expansion and one round of the cell cycle (see out website  and Additional data file 37). At the end of development, major metabolic pathways for lipid metabolism are upregulated, including β-cell oxidation and fatty acid synthesis (Figure 5). Cell development is a dramatic process in which the cell undergoes biochemical and morphological changes. In our study signals at every time point could be detected from more than 14,000 ESTs. Thus, regulating metabolic networks at key points represents an energy efficient way to control cellular processes. Metabolic networks might be activated/inactivated in a similar manner in other types of cellular differentiation, such as myogenesis or osteogenesis. It is intriguing to speculate that signaling networks are also transcriptionally regulated at key points. However, as opposed to the metabolic networks, it is difficult to verify this hypothesis because the key points are not clearly identifiable due to both the interwoven nature and the partial incompleteness of the signaling pathways.
A second general biological insight derived from our global analysis is that we found that many genes were upregulated by well known transcription factors that nevertheless lacked anything resembling the established upstream promoter consensus sites. Over-represented binding sites for key regulators such as PPAR and C/EBP were not detectable with the TRANSFAC matrices. Even the use of a matrix based on all currently available experimentally validated sequences, such as PPAR, did not result in a significant hit. Hence, only few sequences contain this motif in their promoters. These results demonstrate either that much more sophisticated methods must be developed or that there are many cases where the current methods do not perform well because other aspects such as chromatin determine the recognition site.
A third general piece of information derived from our global analysis is the finding that coexpressed genes in fat cell development are not clustered in the genome. Previous studies identified a number of such cases in a range of organisms, including yeast , worms , and flies . This observation of significant correlation in the expression and genomic position of genes was recently reported in the mouse . In the present study we could not identify groups of genes with similar expression profiles, for instance within the same cluster within 5 Mb regions on the chromosomes. Our results suggest that such clustering may not be as widespread as may be presumed by extrapolating from previous studies. However, coexpressed genes could have distant locations and still be spatially colocalized due to DNA looping and banding, as was recently shown in a microscopy study . A higher order chromatin structure of the mammalian transcriptome is an emerging concept , and new methods are required to examine the correlation between gene activity and spatial positioning.
The biological insights gained in this study were only possible with in-depth bioinformatics analyses based on segment and domain predictions. Distribution of GO terms permits a first view of the biological processes, molecular functions, and cellular components. However, in our work more than 40% of the ESTs could not be assigned to GO terms. Moreover, detailed information about the specific functions cannot be extracted. For example, the GO term 'DNA binding' could be specified by 'zinc-finger domain binding protein' only by in-depth analyses. Hence, de novo functional annotation of ESTs using integrated prediction tools and subsequent curation of the results based on the available literature is not only necessary to complete the annotation process but also to reveal the actual biological processes and metabolic networks. Although the number of protein sequences to which a GO term can be assigned is steadily increasing, specific and detailed annotation is only possible with de novo functional annotation.
In the present study we demonstrate that, despite the limitations due to mRNA abundance (many thousands of genes are never transcribed above threshold) and insufficient sensitivity, large-scale gene expression profiling in conjunction with sophisticated bioinformatics analyses can provide not only a list of novel players in a particular setting but also a global view on biological processes and molecular networks.
Materials and methods
The microarray developed here contains 27,648 spots with mouse cDNA clones representing 16,016 different genes (UniGene clusters). These include developmental clones (the 15 K NIA cDNA clone set from National Institute of Aging, US National Institutes of Health) and the 11 K clones from different brain regions in the mouse (Brain Molecular Anatomy Project [BMAP]). Moreover, 627 clones for adipose-related genes were selected using the TIGR Mouse Gene Index Build 5.0 . These cDNA clones were obtained from the IMAGE consortium (Research Genetics, Huntsville, AL, USA). The inserts of the NIA and BMAP clones were sequence verified (insert size about 1-1.5 kb). All PCR products were purified using size exclusion vacuum filter plates (Millipore, Billerica, MA, USA) and spotted onto amino-silanated glass slides (UltraGAPS II; Corning, Corning, NY, USA) in a 4 × 12 print tip group pattern. As spotting buffer 50% dimethyl sulfoxide was used. Negative controls (genomic DNA, genes from Arabidopsis thaliana, and dimethyl sulfoxide) and positive controls (Cot1-DNA and salmon sperm DNA) were included in each of the 48 blocks. Samples were bound to the slides by ultraviolet cross-linking at 200 mJ in a Stratalinker (Stratagene, La Jolla, CA, USA).
3T3-L1 cells (American Type Culture Collection number CL-173) were grown in 100 mm diameter dishes in Dulbecco's modified Eagle's medium supplemented with 10% fetal bovine serum, 100 units/ml penicillin, 100 μg/ml streptomycin, and 2 mmol/l L-glutamine in an atmosphere of 5% carbon dioxide at 37°C. Two days after reaching confluence (day 0), cells were induced to differentiate with a two-day incubation of a hormone cocktail [109, 110] (100 μmol/l 3-iso-butyl-1-methylxanthine, 0.25 μmol/l dexamethasone, 1 μg/ml insulin, 0.16 μmol/l pantothenic acid, and 3.2 μmol/l biotin) added to the standard medium described above. After 48 hours (day 2), cells were cultured in the standard medium in the presence of 1 μg/ml insulin, 0.16 μmol/l pantothenic acid, and 3.2 μmol/l biotin until day 14. Nutrition media were changed every second day.
Three independent cell culture experiments were performed. Cells were harvested and total RNA was isolated at the preconfluent stage and at eight time points (0, 6, 12 and 24 hours, and 3, 4, 7 and 14 days) with TRIzol reagent (Invitrogen-Life Technologies; Carlsbad, CA, USA) . For each independent experiment, RNA was pooled from three different culture dishes for each time point and from 24 dishes at the preconfluent stage used as reference. The quality of the RNA was checked using Agilent 2100 Bioanalyzer RNA assays (Agilent Technologies, Palo Alto, CA, USA) by inspection of the 28S and 18S ribosomal RNA intensity peaks.
Labeling and hybridization
The labeling and hybridization procedures used were based on those developed at the Institute for Genomic Research  and detailed protocols can be viewed on the supplementary website . Briefly, 20 μg total RNA from each time point was reverse transcribed in cDNA and indirectly labeled with Cy5 and 20 μg RNA from the preconfluent stage (reference) was indirectly labeled with Cy3, respectively. This procedure was repeated with reversed dye assignment. Slides were prehybridized with 1% bovine serum albumen. Then, 10 μg mouse Cot1 DNA and 10 μg poly(A) DNA was added to the labeled cDNA samples and pair-wise cohybridized onto the slide for 20 hours at 42°C. Following washing, slides were scanned with a GenePix 4000B microarray scanner (Axon Instruments, Sunnyvale, CA, USA) at 10 μm resolution. Identical photo multiplier voltage settings were used in the scanning of the corresponding dye-swapped hybridized slides. The resulting TIFF images were analyzed with GenePix Pro 4.1 software (Axon Instruments).
Data preprocessing and normalization
Data were filtered for low intensity, inhomogeneity, and saturated spots. To obtain expression values for the saturated spots, slides were scanned a second time with lower photomultiplier tube settings and reanalyzed. All spots of both channels were background corrected (by subtraction of the local background). Different sources of systematic (sample, array, dye, and gene effects) and random errors can be associated with microarray experiments . Nonbiological variation must be removed from the measurement values and the random error can be minimized by normalization [114, 115]. In the present study, gene-wise dye swap normalization was applied. Genes exhibiting substantial differences in intensity ratios between technical replicates were excluded from further analysis based on a two standard deviation cutoff. The resulting ratios were log2 transformed and averaged over three independent experiments. The expression profiles were not rescaled in order to identify genes with high expression values. All experimental parameters, images, and raw and transformed data were uploaded to the microarray database MARS  and submitted via MAGE-ML export to a public repository (ArrayExpress , accession numbers A-MARS-1 and E-MARS-2). Differentially expressed genes were first identified using one-way ANOVA (P < 0.05). They were then subjected to a more stringent criterion; specifically, we considered only those genes with a complete temporal profile that were more than twofold upregulated or downregulated at a minimum of four time points. The twofold cutoff for differentially expressed genes was estimated by applying the significance analyses of microarrays method  to the biological replicates and assuming false discovery rate of 5%. In order to capture the dynamics of various processes, only ESTs differentially expressed in at least half of the time points were selected. Data preprocessing was performed with ArrayNorm .
Microarray expression results were confirmed with RT-PCR. cDNA was synthesized from 2.5 μg total RNA in 20 μl using random hexamers and SuperScript III reverse transcriptase (Invitrogen, Carlsbad, CA, USA). The design of LUX™ primers for Pparg, Lpl, Myc, Dec, Ccna2, and Klf9 was done using the Invitrogen web service (for sequences, see Additional data file 9 and our website ). Quantitative RT-PCR analyses for these genes were performed starting with 50 ng reverse transcribed total RNA, with 0.5× Platinum Quantitative PCR SuperMix-UDG (Invitrogen, Carlsbad, CA, USA), with a ROX reference dye, and with a 200 nmol/l concentration of both LUX™ labeled sense and antisense primers (Invitrogen, Carlsbad, CA, USA) in a 25 μl reaction on an ABI PRISM 7000 sequence detection system (Applied Biosystems, Foster City, CA, USA). To measure PCR efficiency, serial dilutions of reverse transcribed RNA (0.24 pg to 23.8 ng) were amplified. Ribosomal 18S RNA amplifications were used to account for variability in the initial quantities of cDNA. The relative quantification for any given gene with respect to the calibrator (preconfluent stage) was determined using the ΔΔCt method and compared with the normalized ratios resulting from microarray experiments.
Clustering and gene ontology classification
Common unsupervised clustering algorithms  were used for clustering expression profiling of 780 selected ESTs, according to the log ratios from all time points. Using hierarchical clustering the boundaries of the clusters were not clearly separable and required arbitrary determination of the branching point of the tree, whereas the results of the clustering using self-organized maps led to clusters with highly divergent number of ESTs (between 3 and 242). We have therefore used the k means algorithm  and Euclidean distance. The number of clusters was varied from k = 1 to k = 20, and predictive power was analyzed with the figure of merit . Subsequently, k = 12 was found to be optimal. To evaluate the results of the k means clustering, principal component analysis  was applied and exhibited low intracluster distances and high intercluster dissimilarities. GO terms and GO numbers for molecular function, biological process, and cellular components were derived from the Gene Ontology database (Gene Ontology Consortium) using the GenPept/RefSeq accession numbers for annotated proteins encoded by selected genes (ESTs). All cluster analyses and visualizations were performed using Genesis .
De novoannotation of ESTs
For each of the 780 selected EST sequences, we attempted to find the corresponding protein sequence. Megablast  searches (word length w = 70, percentage identity p = 95%) against nucleotide databases (in the succession of RefSeq [126, 127], FANTOM , UniGene , nr GenBank, and TIGR Mouse Gene Index  until a gene hit was found) were carried out. For the ESTs still remaining without gene assignment, new Megablast searches were conducted with the largest compilation of RefSeq (including the provisional and automatically generated records [126, 127]). If an EST remained unassigned, then the whole procedure was repeated with blastn . In addition, a blastn search against the ENSEMBL mouse genome  was performed, and ESTs with long stretches (>100 base pairs) of unspecified nucleotides (N) were excluded.
All protein sequences were annotated de novo with academic prediction tools that are integrated into ANNOTATOR, a novel protein sequence analysis system : compositional bias (SAPS , Xnu, Cast , GlobPlot 1.2 ); low complexity regions (SEG ); known sequence domains (Pfam , Smart , Prosite and Prosite pattern  with HMMER, RPS-BLAST , IMPALA , PROSITE-Profile ); transmembrane domains (HMMTOP 2.0 , TOPPRED , DAS-TMfilter , SAPS ); secondary structures (impCOIL , Predator , SSCP [147, 148]); targeting signals (SIGCLEAVE , SignalP-3.0 , PTS1 ); post-translational modifications (big-PI , NMT , Prenylation); a series of small sequence motifs (ELM, Prosite patterns , BioMotif-IMPlibrary); and homology searches with NCBI blast . Further information was retrieved from the databases of Mouse Genome Informatics  and LocusLink .
The promoters were retrieved from PromoSer database  through the gene accession number. PromoSer contains 22,549 promoters for 12,493 unique genes. Nucleotides from 2,000 upstream and 100 downstream of the transcription start site were obtained. With an implementation of the MatInspector algorithm , the Transfac matrices  were checked for binding sites in the promoter regions with a threshold for matrix similarity of 0.85. We counted the number of those gene sequences that were found to carry a predicted transcription factor binding site. As a reference set all unique genes of the PromoSer were reanalyzed. A one-sided χ2 test and a one-sided Fisher's exact test (to improve the statistics for view counts) were performed with the statistical tool R  to determine the clusters with a higher affinity for a transcription factor.
Identification of miRNA target sites in 3'-UTR
All available 3'-UTR sequences (21,396) for mouse genes were derived with EnsMart , using Ensembl gene build for the NCBI m33 mouse assembly. 3'-UTRs for unique genes represented by the 780 selected ESTs were extracted using Ensembl transcript ID. A total of 234 mouse miRNA sequences were derived from the Rfam database . The 3'-UTR sequences were searched for antisense matches to the designated seed region of each miRNA (bases 1-8, 2-8, 1-9, and 2-9 starting from the 5' end). Significantly over-represented miRNA motifs in each cluster in comparison with the remaining motifs in the whole 3'-UTR sequence set were determined using the one-sided Fisher's exact test (significance level: P < 0.05) and miRNA targets of all clusters were analyzed for significantly over-represented miRNAs.
Chromosomal localization analysis
RefSeq sequences for 780 selected ESTs, shown to be more than two times upregulated or downregulated in a minimum of four time points during adipocyte differentiation and clustered according their expression profiles, were mapped onto the chromosomes from the NCBI Mus musculus genome (build 33) using ChromoMapper 2.1.0 software  based on MegaBlast with the following parameters: 99% identity cutoff, word size 32, and E-value (0.001). Colocalized sequences of all selected ESTs and from each of the 12 clusters within a 5 Mb genomic interval were identified. Within the 5Mb genomic intervals of each chromosome with the highest density of mapped ESTs, relative gene expression levels (log2 ratios) of these ESTs at different time points were related to the genomic localization.
Additional data files
The following additional data are provided with the online version of this article: A spot map for the array design (Additional data file 1); a fasta file containing the EST sequences used for the array (Additional data file 2); an Excel file containing expression values for the 780 selected ESTs (Additional data file 3); an Excel file containing expression values for the 5205 ESTs filtered with ANOVA (Additional data file 4); GenePix result files containing raw data (Additional data file 5); images showing the distribution of gene ontology (Additional data file 6); a table listing relevant proteins (Additional data file 7); a fasta file containing sequences of the relevant proteins (Additional data file 8); a pdf file containing real-time RT-PCR data (Additional data file 9); a table including statistical analysis of independent experiments (Additional data file 10); figure showing a comparison with GeneAtlas (Additional data file 11); a table including expression levels from the present study and GeneAtlas (Additional data file 12); figure showing a comparison with the data set reported by Soukas and coworkers Additional data file 13); figure showing a comparison with the data set reported by Ross and coworkers  (Additional data file 14); figure showing a comparison with the data set reported by Burton and coworkers  (Additional data file 15); a table including ESTs unique to the present study (Additional data file 16); a figure showing genes with miRNA motifs in 3'-UTR (Additional data file 17); a figure illustrating the significant over-representation of miRNA motifs in the 3'-UTR of genes in each cluster (Additional data file 18); figures showing the significant over-representation of miRNA motifs in the 3'-UTR from genes in each cluster (Additional data files 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29); a table including over-represented miRNA motifs in the 3'-UTR from genes in the set of 780 selected ESTs (Additional data file 30); text describing regulation of metabolic pathways (Additional data file 31); figure showing regulation of metabolic pathways by key points (Additional data file 32); figure showing the cellular localization of gene products involved in metabolism and their gene expression at different time points (Additional data file 33); figure showing the cellular localization of gene products involved in other biological processes and their gene expression at different time points (Additional data file 30); text describing signaling networks (Additional data file 35); text file describing extracellular matrix remodeling and cytoskeleton reorganization (Additional data file 36); figure showing cell cycle processes (Additional data file 37); figure showing the cholesterol pathway (Additional data file 38); a list of experimental verified binding site for PPAR:RXR and the derived position weight matrix (Additional data file 39); text file containing TRANSFAC matrices for vertebrates (Additional data file 40); a file showing the promoter sequences in fasta format (Additional data file 41); figure showing cluster-wise mapping of 780 ESTs to all chromosomes (Additional data file 42); figure showing expression of colocalized ESTs for each cluster (Additional data file 43); an Excel file showing a statistical analysis of colocalized ESTs for 780 selected ESTs (Additional data file 44); and an Excel file showing a statistical analysis of colocalized ESTs for 5,502 ANOVA selected ESTs (Additional data file 45).
Green H, Meuth M: An established pre-adipose cell line and its differentiation in culture. Cell. 1974, 3: 127-133. 10.1016/0092-8674(74)90116-0.
Macdougald OA, Lane MD: Transcriptional regulation of gene expression during adipocyte differentiation. Annu Rev Biochem. 1995, 64: 345-373. 10.1146/annurev.bi.64.070195.002021.
Yeh WC, Cao Z, Classon M, McKnight SL: Cascade regulation of terminal adipocyte differentiation by three members of the C/EBP family of leucine zipper proteins. Genes Dev. 1995, 9: 168-181.
Tanaka T, Yoshida N, Kishimoto T, Akira S: Defective adipocyte differentiation in mice lacking the C/EBPbeta and/or C/EBPdelta gene. EMBO J. 1997, 16: 7432-7443. 10.1093/emboj/16.24.7432.
Kim JB, Spiegelman BM: ADD1/SREBP1 promotes adipocyte differentiation and gene expression linked to fatty acid metabolism. Genes Dev. 1996, 10: 1096-1107.
Fajas L, Schoonjans K, Gelman L, Kim JB, Najib J, Martin G, Fruchart JC, Briggs M, Spiegelman BM, Auwerx J: Regulation of peroxisome proliferator-activated receptor gamma expression by adipocyte differentiation and determination factor 1/sterol regulatory element binding protein 1: implications for adipocyte differentiation and metabolism. Mol Cell Biol. 1999, 19: 5495-5503.
Kim JB, Wright HM, Wright M, Spiegelman BM: ADD1/SREBP1 activates PPARgamma through the production of endogenous ligand. Proc Natl Acad Sci USA. 1998, 95: 4333-4337. 10.1073/pnas.95.8.4333.
Soukas A, Socci ND, Saatkamp BD, Novelli S, Friedman JM: Distinct transcriptional profiles of adipogenesis in vivo and in vitro. J Biol Chem. 2001, 276: 34167-34174. 10.1074/jbc.M104421200.
Ross SE, Erickson RL, Gerin I, DeRose PM, Bajnok L, Longo KA, Misek DE, Kuick R, Hanash SM, Atkins KB, et al: Microarray analyses during adipogenesis: understanding the effects of Wnt signaling on adipogenesis and the roles of liver X receptor alpha in adipocyte metabolism. Mol Cell Biol. 2002, 22: 5989-5999. 10.1128/MCB.22.16.5989-5999.2002.
Burton GR, Guan Y, Nagarajan R, McGehee RE: Microarray analysis of gene expression during early adipocyte differentiation. Gene. 2002, 293: 21-31. 10.1016/S0378-1119(02)00726-6.
Burton GR, McGehee REJ: Identification of candidate genes involved in the regulation of adipocyte differentiation using microarray-based gene expression profiling. Nutrition. 2004, 20: 109-114. 10.1016/j.nut.2003.09.019.
Burton GR, Nagarajan R, Peterson CA, McGehee REJ: Microarray analysis of differentiation-specific gene expression during 3T3-L1 adipogenesis. Gene. 2004, 329: 167-185. 10.1016/j.gene.2003.12.012.
Jessen BA, Stevens GJ: Expression profiling during adipocyte differentiation of 3T3-L1 fibroblasts. Gene. 2002, 299: 95-100. 10.1016/S0378-1119(02)01017-X.
Gerhold DL, Liu F, Jiang G, Li Z, Xu J, Lu M, Sachs JR, Bagchi A, Fridman A, Holder DJ, et al: Gene expression profile of adipocyte differentiation and its regulation by peroxisome proliferator-activated receptor-gamma agonists. Endocrinology. 2002, 143: 2106-2118. 10.1210/en.143.6.2106.
Guo X, Liao K: Analysis of gene expression profile during 3T3-L1 preadipocyte differentiation. Gene. 2000, 251: 45-53. 10.1016/S0378-1119(00)00192-X.
Ko MS, Kitchen JR, Wang X, Threat TA, Wang X, Hasegawa A, Sun T, Grahovac MJ, Kargul GJ, Lim MK, et al: Large-scale cDNA analysis reveals phased gene expression patterns during preimplantation mouse development. Development. 2000, 127: 1737-1749.
Larkin JE, Frank BC, Gaspard RM, Duka I, Gavras H, Quackenbush J: Cardiac transcriptional response to acute and chronic angiotensin II treatments. Physiol Genomics. 2004, 18: 152-166. 10.1152/physiolgenomics.00057.2004.
Tanaka TS, Jaradat SA, Lim MK, Kargul GJ, Wang X, Grahovac MJ, Pantano S, Sano Y, Piao Y, Nagaraja R, et al: Genome-wide expression profiling of mid-gestation placenta and embryo using a 15,000 mouse developmental cDNA microarray. Proc Natl Acad Sci USA. 2000, 97: 9127-9132. 10.1073/pnas.97.16.9127.
Quackenbush J, Liang F, Holt I, Pertea G, Upton J: The TIGR gene indices: reconstruction and representation of expressed gene sequences. Nucleic Acids Res. 2000, 28: 141-145. 10.1093/nar/28.1.141.
Molecular processes during fat cell development revealed by gene expression profiling and functional annotation. [http://genome.tugraz.at/fatcell]
Su AI, Wiltshire T, Batalov S, Lapp H, Ching KA, Block D, Zhang J, Soden R, Hayakawa M, Kreiman G, et al: A gene atlas of the mouse and human protein-encoding transcriptomes. Proc Natl Acad Sci U S A. 2004, 101: 6062-6067. 10.1073/pnas.0400782101.
Tang QQ, Otto TC, Lane MD: Mitotic clonal expansion: A synchronous process required for adipogenesis. Proc Natl Acad Sci U S A. 2003, 100: 44-49. 10.1073/pnas.0137044100.
Xie X, Lu J, Kulbokas EJ, Golub TR, Mootha V, Lindblad-Toh K, Lander ES, Kellis M: Systematic discovery of regulatory motifs in human promoters and 3' UTRs by comparison of several mammals. Nature. 2005, 434: 338-345. 10.1038/nature03441.
John B, Enright AJ, Aravin A, Tuschl T, Sander C, Marks DS: Human MicroRNA targets. PLoS Biol. 2004, 2: e363-10.1371/journal.pbio.0020363.
Doench JG, Sharp PA: Specificity of microRNA target selection in translational repression. Genes Dev. 2004, 18: 504-511. 10.1101/gad.1184404.
Hutvagner G, Simard MJ, Mello CC, Zamore PD: Sequence-specific inhibition of small RNA function. PLoS Biol. 2004, 2: E98-10.1371/journal.pbio.0020098.
Lewis BP, Shih IH, Jones-Rhoades MW, Bartel DP, Burge CB: Prediction of mammalian microRNA targets. Cell. 2003, 115: 787-798. 10.1016/S0092-8674(03)01018-3.
Esau C, Kang X, Peralta E, Hanson E, Marcusson EG, Ravichandran LV, Sun Y, Koo S, Perera RJ, Jain R, et al: MicroRNA-143 regulates adipocyte differentiation. J Biol Chem. 2004, 279: 52361-52365. 10.1074/jbc.C400438200.
Johnson SM, Grosshans H, Shingara J, Byrom M, Jarvis R, Cheng A, Labourier E, Reinert KL, Brown D, Slack FJ: RAS is regulated by the let-7 microRNA family. Cell. 2005, 120: 635-647. 10.1016/j.cell.2005.01.014.
Mlecnik B, Scheideler M, Hackl H, Hartler J, Sanchez-Cabo F, Trajanoski Z: PathwayExplorer: web service for visualizing high-throughput expression data on biological pathways. Nucleic Acids Res. 2005, 33: W633-W637. 10.1093/nar/gki391.
Tontonoz P, Hu E, Spiegelman BM: Stimulation of adipogenesis in fibroblasts by PPAR gamma 2, a lipid-activated transcription factor. Cell. 1994, 79: 1147-1156. 10.1016/0092-8674(94)90006-X.
Inoue J, Kumagai H, Terada T, Maeda M, Shimizu M, Sato R: Proteolytic activation of SREBPs during adipocyte differentiation. Biochem Biophys Res Commun. 2001, 283: 1157-1161. 10.1006/bbrc.2001.4915.
Yang T, Espenshade PJ, Wright ME, Yabe D, Gong Y, Aebersold R, Goldstein JL, Brown MS: Crucial step in cholesterol homeostasis: sterols promote binding of SCAP to INSIG-1, a membrane protein that facilitates retention of SREBPs in ER. Cell. 2002, 110: 489-500. 10.1016/S0092-8674(02)00872-3.
Kast-Woelbern HR, Dana SL, Cesario RM, Sun L, de Grandpre LY, Brooks ME, Osburn DL, Reifel-Miller A, Klausing K, Leibowitz MD: Rosiglitazone induction of Insig-1 in white adipose tissue reveals a novel interplay of peroxisome proliferator-activated receptor gamma and sterol regulatory element-binding protein in the regulation of adipogenesis. J Biol Chem. 2004, 279: 23908-23915. 10.1074/jbc.M403145200.
Zimmermann R, Strauss JG, Haemmerle G, Schoiswohl G, Birner-Gruenberger R, Riederer M, Lass A, Neuberger G, Eisenhaber F, Hermetter A, Zechner R: Fat mobilization in adipose tissue is promoted by adipose triglyceride lipase. Science. 2004, 306: 1383-1386. 10.1126/science.1100747.
Fukuhara A, Matsuda M, Nishizawa M, Segawa K, Tanaka M, Kishimoto K, Matsuki Y, Murakami M, Ichisaka T, Murakami H, et al: Visfatin: a protein secreted by visceral fat that mimics the effects of insulin. Science. 2005, 307: 426-430. 10.1126/science.1097243.
Kitani T, Okuno S, Fujisawa H: Growth phase-dependent changes in the subcellular localization of pre-B-cell colony-enhancing factor. FEBS Lett. 2003, 544: 74-78. 10.1016/S0014-5793(03)00476-9.
Revollo JR, Grimm AA, Imai S: The NAD biosynthesis pathway mediated by nicotinamide phosphoribosyltransferase regulates Sir2 activity in mammalian cells. J Biol Chem. 2004, 279: 50754-50763. 10.1074/jbc.M408388200.
Banerjee SS, Feinberg MW, Watanabe M, Gray S, Haspel RL, Denkinger DJ, Kawahara R, Hauner H, Jain MK: The Kruppel-like factor KLF2 inhibits peroxisome proliferator-activated receptor-gamma expression and adipogenesis. J Biol Chem. 2003, 278: 2581-2584. 10.1074/jbc.M210859200.
Oishi Y, Manabe I, Tobe K, Tsushima K, Shindo T, Fujiu K, Nishimura G, Maemura K, Yamauchi T, Kubota N, et al: Krüppel-like transcription factor KLF5 is a key regulator of adipocyte differentiation. Cell Metab. 2005, 1: 27-39. 10.1016/j.cmet.2004.11.005.
Li D, Yea S, Li S, Chen Z, Narla G, Banck M, Laborda J, Tan S, Friedman JM, Friedman SL, Walsh MJ: Kruppel-like factor-6 promotes preadipocyte differentiation through histone deacetylase 3-dependent repression of DLK1. J Biol Chem. 2005, 280: 26941-26952. 10.1074/jbc.M500463200.
Mori T, Sakaue H, Iguchi H, Gomi H, Okada Y, Takashima Y, Nakamura K, Nakamura T, Yamauchi T, Kubota N, et al: Role of Kruppel-like factor 15 (KLF15) in transcriptional regulation of adipogenesis. J Biol Chem. 2005, 280: 12867-12875. 10.1074/jbc.M410515200.
Ghaleb AM, Nandan MO, Chanchevalap S, Dalton WB, Hisamuddin IM, Yang VW: Kruppel-like factors 4 and 5: the yin and yang regulators of cellular proliferation. Cell Res. 2005, 15: 92-96.
Nakae J, Kitamura T, Kitamura Y, Biggs WH, Arden KC, Accili D: The forkhead transcription factor foxo1 regulates adipocyte differentiation. Dev Cell. 2003, 4: 119-129. 10.1016/S1534-5807(02)00401-X.
Farmer SR: The forkhead transcription factor Foxo1: a possible link between obesity and insulin resistance. Mol Cell. 2003, 11: 6-8. 10.1016/S1097-2765(03)00003-0.
D'Adamio F, Zollo O, Moraca R, Ayroldi E, Bruscoli S, Bartoli A, Cannarile L, Migliorati G, Riccardi C: A new dexamethasone-induced gene of the leucine zipper family protects T lymphocytes from TCR/CD3-activated cell death. Immunity. 1997, 7: 803-812. 10.1016/S1074-7613(00)80398-2.
Shi X, Shi W, Li Q, Song B, Wan M, Bai S, Cao X: A glucocorticoid-induced leucine-zipper protein, GILZ, inhibits adipogenesis of mesenchymal cells. EMBO Rep. 2003, 4: 374-380. 10.1038/sj.embor.embor805.
Xie J, Cai T, Zhang H, Lan MS, Notkins AL: The zinc-finger transcription factor INSM1 is expressed during embryo development and interacts with the Cbl-associated protein. Genomics. 2002, 80: 54-61. 10.1006/geno.2002.6800.
Zhu M, Breslin MB, Lan MS: Expression of a novel zinc-finger cDNA, IA-1, is associated with rat AR42J cells differentiation into insulin-positive cells. Pancreas. 2002, 24: 139-145. 10.1097/00006676-200203000-00004.
Yamada K, Printz RL, Osawa H, Granner DK: Human ZHX1: cloning, chromosomal location, and interaction with transcription factor NF-Y. Biochem Biophys Res Commun. 1999, 261: 614-621. 10.1006/bbrc.1999.1087.
Hebrok M, Wertz K, Fuchtbauer EM: M-twist is an inhibitor of muscle differentiation. Dev Biol. 1994, 165: 537-544. 10.1006/dbio.1994.1273.
Sosic D, Richardson JA, Yu K, Ornitz DM, Olson EN: Twist regulates cytokine gene expression through a negative feedback loop that represses NF-kappaB activity. Cell. 2003, 112: 169-180. 10.1016/S0092-8674(03)00002-3.
Chae GN, Kwak SJ: NF-kappaB is involved in the TNF-alpha induced inhibition of the differentiation of 3T3-L1 cells by reducing PPARgamma expression. Exp Mol Med. 2003, 35: 431-437.
Ku DH, Chang CD, Koniecki J, Cannizzaro LA, Boghosian-Sell L, Alder H, Baserga R: A new growth-regulated complementary DNA with the sequence of a putative trans-activating factor. Cell Growth Differ. 1991, 2: 179-186.
Metzger D, Scheer E, Soldatov A, Tora L: Mammalian TAF(II)30 is required for cell cycle progression and specific cellular differentiation programmes. EMBO J. 1999, 18: 4823-4834. 10.1093/emboj/18.17.4823.
Novatchkova M, Eisenhaber F: Can molecular mechanisms of biological processes be extracted from expression profiles? Case study: endothelial contribution to tumor-induced angiogenesis. Bioessays. 2001, 23: 1159-1175. 10.1002/bies.10013.
Croissandeau G, Chretien M, Mbikay M: Involvement of matrix metalloproteinases in the adipose conversion of 3T3-L1 preadipocytes. Biochem J. 2002, 364: 739-746. 10.1042/BJ20011158.
Karagiannis ED, Popel AS: A theoretical model of type I collagen proteolysis by matrix metalloproteinase (MMP) 2 and membrane type 1 MMP in the presence of tissue inhibitor of metalloproteinase 2. J Biol Chem. 2004, 279: 39105-39114. 10.1074/jbc.M403627200.
Chavey C, Mari B, Monthouel MN, Bonnafous S, Anglard P, Van Obberghen E, Tartare-Deckert S: Matrix metalloproteinases are differentially expressed in adipose tissue during obesity and modulate adipocyte differentiation. J Biol Chem. 2003, 278: 11888-11896. 10.1074/jbc.M209196200.
Weiner FR, Shah A, Smith PJ, Rubin CS, Zern MA: Regulation of collagen gene expression in 3T3-L1 cells. Effects of adipocyte differentiation and tumor necrosis factor alpha. Biochemistry. 1989, 28: 4094-4099. 10.1021/bi00435a070.
Dimaculangan DD, Chawla A, Boak A, Kagan HM, Lazar MA: Retinoic acid prevents downregulation of ras recision gene/lysyl oxidase early in adipocyte differentiation. Differentiation. 1994, 58: 47-52. 10.1046/j.1432-0436.1994.5810047.x.
Piecha D, Wiberg C, Morgelin M, Reinhardt DP, Deak F, Maurer P, Paulsson M: Matrilin-2 interacts with itself and with other extracellular matrix proteins. Biochem J. 2002, 367: 715-721. 10.1042/BJ20021069.
Brekken RA, Sage EH: SPARC, a matricellular protein: at the crossroads of cell-matrix. Matrix Biol. 2000, 19: 569-580. 10.1016/S0945-053X(00)00105-0.
Bradshaw AD, Sage EH: SPARC, a matricellular protein that functions in cellular differentiation and tissue response to injury. J Clin Invest. 2001, 107: 1049-1054.
Spiegelman BM, Farmer SR: Decreases in tubulin and actin gene expression prior to morphological differentiation of 3T3 adipocytes. Cell. 1982, 29: 53-60. 10.1016/0092-8674(82)90089-7.
Edwards RA, Herrera-Sosa H, Otto J, Bryan J: Cloning and expression of a murine fascin homolog from mouse brain. J Biol Chem. 1995, 270: 10764-10770. 10.1074/jbc.270.18.10764.
Winder SJ, Jess T, Ayscough KR: SCP1 encodes an actin-bundling protein in yeast. Biochem J. 2003, 375: 287-295. 10.1042/BJ20030796.
Hossain MM, Hwang DY, Huang QQ, Sasaki Y, Jin JP: Developmentally regulated expression of calponin isoforms and the effect of h2-calponin on cell proliferation. Am J Physiol Cell Physiol. 2003, 284: C156-C167.
He HJ, Kole S, Kwon YK, Crow MT, Bernier M: Interaction of filamin A with the insulin receptor alters insulin-dependent activation of the mitogen-activated protein kinase pathway. J Biol Chem. 2003, 278: 27096-27104. 10.1074/jbc.M301003200.
Oakley BR: Gamma-tubulin: the microtubule organizer?. Trends Cell Biol. 1992, 2: 1-5. 10.1016/0962-8924(92)90125-7.
Honda K, Yamada T, Endo R, Ino Y, Gotoh M, Tsuda H, Yamada Y, Chiba H, Hirohashi S: Actinin-4, a novel actin-bundling protein associated with cell motility and cancer invasion. J Cell Biol. 1998, 140: 1383-1393. 10.1083/jcb.140.6.1383.
Leeuwen FN, Kain HE, Kammen RA, Michiels F, Kranenburg OW, Collard JG: The guanine nucleotide exchange factor Tiam1 affects neuronal morphology; opposing roles for the small GTPases Rac and Rho. J Cell Biol. 1997, 139: 797-807. 10.1083/jcb.139.3.797.
Sander EE, ten Klooster JP, van Delft S, van der Kammen RA, Collard JG: Rac downregulates Rho activity: reciprocal balance between both GTPases determines cellular morphology and migratory behavior. J Cell Biol. 1999, 147: 1009-1022. 10.1083/jcb.147.5.1009.
Saras J, Wollberg P, Aspenstrom P: Wrch1 is a GTPase-deficient Cdc42-like protein with unusual binding characteristics and cellular effects. Exp Cell Res. 2004, 299: 356-369. 10.1016/j.yexcr.2004.05.029.
Shutes A, Berzat AC, Cox AD, Der CJ: Atypical mechanism of regulation of the Wrch-1 Rho family small GTPase. Curr Biol. 2004, 14: 2052-2056. 10.1016/j.cub.2004.11.011.
Klipp E, Heinrich R, Holzhutter HG: Prediction of temporal gene expression. Metabolic opimization by re-distribution of enzyme activities. Eur J Biochem. 2002, 269: 5406-5413. 10.1046/j.1432-1033.2002.03223.x.
Lynen F: Acetyl coenzyme A and the fatty acid cycle. Harvey Lect. 1952, 48: 210-244.
Ganguly J: Studies on the mechanism of fatty acid synthesis. VII. Biosynthesis of fatty acids from malonyl CoA. Biochim Biophys Acta. 1960, 40: 110-118. 10.1016/0006-3002(60)91320-2.
Song WJ, Jackowski S: Kinetics and regulation of pantothenate kinase from Escherichia coli. J Biol Chem. 1994, 269: 27051-27058.
Rongvaux A, Shea RJ, Mulks MH, Gigot D, Urbain J, Leo O, Andris F: Pre-B-cell colony-enhancing factor, whose expression is up-regulated in activated lymphocytes, is a nicotinamide phosphoribosyltransferase, a cytosolic enzyme involved in NAD biosynthesis. Eur J Immunol. 2002, 32: 3225-3234. 10.1002/1521-4141(200211)32:11<3225::AID-IMMU3225>3.0.CO;2-L.
Clarke JL, Mason PJ: Murine hexose-6-phosphate dehydrogenase: a bifunctional enzyme with broad substrate specificity and 6-phosphogluconolactonase activity. Arch Biochem Biophys. 2003, 415: 229-234. 10.1016/S0003-9861(03)00229-7.
Enoch HG, Catala A, Strittmatter P: Mechanism of rat liver microsomal stearyl-CoA desaturase. Studies of the substrate specificity, enzyme-substrate interactions, and the function of lipid. J Biol Chem. 1976, 251: 5095-5103.
Ntambi JM: Regulation of stearoyl-CoA desaturase by polyunsaturated fatty acids and cholesterol. J Lipid Res. 1999, 40: 1549-1558.
Moon YA, Shah NA, Mohapatra S, Warrington JA, Horton JD: Identification of a mammalian long chain fatty acyl elongase regulated by sterol regulatory element-binding proteins. J Biol Chem. 2001, 276: 45358-45366. 10.1074/jbc.M108413200.
Nilsson-Ehle P: Impaired regulation of adipose tissue lipoprotein lipase in obesity. Int J Obes. 1981, 5: 695-699.
Semenkovich CF, Wims M, Noe L, Etienne J, Chan L: Insulin regulation of lipoprotein lipase activity in 3T3-L1 adipocytes is mediated at posttranscriptional and posttranslational levels. J Biol Chem. 1989, 264: 9030-9038.
Koike T, Liang J, Wang X, Ichikawa T, Shiomi M, Liu G, Sun H, Kitajima S, Morimoto M, Watanabe T, et al: Overexpression of lipoprotein lipase in transgenic Watanabe heritable hyperlipidemic rabbits improves hyperlipidemia and obesity. J Biol Chem. 2004, 279: 7521-7529. 10.1074/jbc.M311514200.
Zhang J, Zhang W, Zou D, Chen G, Wan T, Zhang M, Cao X: Cloning and functional characterization of ACAD-9, a novel member of human acyl-CoA dehydrogenase family. Biochem Biophys Res Commun. 2002, 297: 1033-1042. 10.1016/S0006-291X(02)02336-7.
Harris RA, Hawes JW, Popov KM, Zhao Y, Shimomura Y, Sato J, Jaskiewicz J, Hurley TD: Studies on the regulation of the mitochondrial alpha-ketoacid dehydrogenase complexes and their kinases. Adv Enzyme Regul. 1997, 37: 271-293. 10.1016/S0065-2571(96)00009-X.
Clark DV, MacAfee N: The purine biosynthesis enzyme PRAT detected in proenzyme and mature forms during development of Drosophila melanogaster. Insect Biochem Mol Biol. 2000, 30: 315-323. 10.1016/S0965-1748(00)00005-9.
Bohman C, Eriksson S: Deoxycytidine kinase from human leukemic spleen: preparation and characteristics of homogeneous enzyme. Biochemistry. 1988, 27: 4258-4265. 10.1021/bi00412a009.
Hatzis P, Al Madhoon AS, Jullig M, Petrakis TG, Eriksson S, Talianidis I: The intracellular localization of deoxycytidine kinase. J Biol Chem. 1998, 273: 30239-30243. 10.1074/jbc.273.46.30239.
Sabini E, Ort S, Monnerjahn C, Konrad M, Lavie A: Structure of human dCK suggests strategies to improve anticancer and antiviral therapy. Nat Struct Biol. 2003, 10: 513-519. 10.1038/nsb942.
Wright JA, Chan AK, Choy BK, Hurta RA, McClarty GA, Tagger AY: Regulation and drug resistance mechanisms of mammalian ribonucleotide reductase, and the significance to DNA synthesis. Biochem Cell Biol. 1990, 68: 1364-1371.
Dong Z, Liu LH, Han B, Pincheira R, Zhang JT: Role of eIF3 p170 in controlling synthesis of ribonucleotide reductase M2 and cell growth. Oncogene. 2004, 23: 3790-3801. 10.1038/sj.onc.1207465.
Xu P, Huecksteadt TP, Harrison R, Hoidal JR: Molecular cloning, tissue expression of human xanthine dehydrogenase. Biochem Biophys Res Commun. 1994, 199: 998-1004. 10.1006/bbrc.1994.1328.
Xu P, Huecksteadt TP, Hoidal JR: Molecular cloning and characterization of the human xanthine dehydrogenase gene (XDH). Genomics. 1996, 34: 173-180. 10.1006/geno.1996.0262.
Popplewell PY, Azhar S: Effects of aging on cholesterol content and cholesterol-metabolizing enzymes in the rat adrenal gland. Endocrinology. 1987, 121: 64-73.
Sato R, Takano T: Regulation of intracellular cholesterol metabolism. Cell Struct Funct. 1995, 20: 421-427.
Matys V, Fricke E, Geffers R, Gossling E, Haubrock M, Hehl R, Hornischer K, Karas D, Kel AE, Kel-Margoulis OV, et al: TRANSFAC: transcriptional regulation, from patterns to profiles. Nucleic Acids Res. 2003, 31: 374-378. 10.1093/nar/gkg108.
Kim JB, Spotts GD, Halvorsen YD, Shih HM, Ellenberger T, Towle HC, Spiegelman BM: Dual DNA binding specificity of ADD1/SREBP1 controlled by a single amino acid in the basic helix-loop-helix domain. Mol Cell Biol. 1995, 15: 2582-2588.
Yokoyama C, Wang X, Briggs MR, Admon A, Wu J, Hua X, Goldstein JL, Brown MS: SREBP-1, a basic-helix-loop-helix-leucine zipper protein that controls transcription of the low density lipoprotein receptor gene. Cell. 1993, 75: 187-197. 10.1016/0092-8674(93)90690-R.
Salomonis N, Cotte N, Zambon AC, Pollard KS, Vranizan K, Doniger SW, Dolganov G, Conklin BR: Identifying genetic networks underlying myometrial transition to labor. Genome Biol. 2005, 6: R12-10.1186/gb-2005-6-2-r12.
Cohen BA, Mitra RD, Hughes JD, Church GM: A computational analysis of whole-genome expression data reveals chromosomal domains of gene expression. Nat Genet. 2000, 26: 183-186. 10.1038/79896.
Roy PJ, Stuart JM, Lund J, Kim SK: Chromosomal clustering of muscle-expressed genes in Caenorhabditis elegans. Nature. 2002, 418: 975-979.
Spellman PT, Rubin GM: Evidence for large domains of similarly expressed genes in the Drosophila genome. J Biol. 2002, 1: 5-10.1186/1475-4924-1-5.
Osborne CS, Chakalova L, Brown KE, Carter D, Horton A, Debrand E, Goyenechea B, Mitchell JA, Lopes S, Reik W, Fraser P: Active genes dynamically colocalize to shared sites of ongoing transcription. Nat Genet. 2004, 36: 1065-1071. 10.1038/ng1423.
Oliver B, Misteli T: A non-random walk through the genome. Genome Biol. 2005, 6: 214-10.1186/gb-2005-6-4-214.
Student AK, Hsu RY, Lane MD: Induction of fatty acid synthetase synthesis in differentiating 3T3-L1 preadipocytes. J Biol Chem. 1980, 255: 4745-4750.
Le Lay S, Lefrere I, Trautwein C, Dugail I, Krief S: Insulin and sterol-regulatory element-binding protein-1c (SREBP-1C) regulation of gene expression in 3T3-L1 adipocytes. Identification of CCAAT/enhancer-binding protein beta as an SREBP-1C target. J Biol Chem. 2002, 277: 35625-35634. 10.1074/jbc.M203913200.
Chomczynski P, Sacchi N: Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. Anal Biochem. 1987, 162: 156-159. 10.1016/0003-2697(87)90021-2.
Hegde P, Qi R, Abernathy K, Gay C, Dharap S, Gaspard R, Hughes JE, Snesrud E, Lee N, Quackenbush J: A concise guide to cDNA microarray analysis. Biotechniques. 2000, 29: 548-556.
Kerr MK, Martin M, Churchill GA: Analysis of variance for gene expression microarray data. J Comput Biol. 2000, 7: 819-837. 10.1089/10665270050514954.
Quackenbush J: Microarray data normalization and transformation. Nat Genet. 2002, 32 (Suppl): 496-501. 10.1038/ng1032.
Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J, Speed TP: Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Res. 2002, 30: e15-10.1093/nar/30.4.e15.
Maurer M, Molidor R, Sturn A, Hartler J, Hackl H, Stocker G, Prokesch A, Scheideler M, Trajanoski Z: MARS: Microarray analysis, retrieval and storage system. BMC Bioinformatics. 2005, 6: 101-10.1186/1471-2105-6-101.
Brazma A, Parkinson H, Sarkans U, Shojatalab M, Vilo J, Abeygunawardena N, Holloway E, Kapushesky M, Kemmeren P, Lara GG, et al: ArrayExpress--a public repository for microarray gene expression data at the EBI. Nucleic Acids Res. 2003, 31: 68-71. 10.1093/nar/gkg091.
Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci USA. 2001, 98: 5116-5121. 10.1073/pnas.091062498.
Pieler R, Sanchez-Cabo F, Hackl H, Thallinger GG, Trajanoski Z: ArrayNorm: comprehensive normalization and analysis of microarray data. Bioinformatics. 2004
Quackenbush J: Computational analysis of microarray data. Nat Rev Genet. 2001, 2: 418-427. 10.1038/35076576.
Hartigan JA: Clustering Algorithms. 1975, New York: Wiley & Sons
Yeung KY, Haynor DR, Ruzzo WL: Validating clustering for gene expression data. Bioinformatics. 2001, 17: 309-318. 10.1093/bioinformatics/17.4.309.
Raychaudhuri S, Stuart JM, Altman RB: Principal components analysis to summarize microarray experiments: application to sporulation time series. Pac Symp Biocomput. 2000, 2000: 455-466.
Sturn A, Quackenbush J, Trajanoski Z: Genesis: cluster analysis of microarray data. Bioinformatics. 2002, 18: 207-208. 10.1093/bioinformatics/18.1.207.
Zhang Z, Schwartz S, Wagner L, Miller W: A greedy algorithm for aligning DNA sequences. J Comput Biol. 2000, 7: 203-214. 10.1089/10665270050081478.
Pruitt KD, Katz KS, Sicotte H, Maglott DR: Introducing RefSeq and LocusLink: curated human genome resources at the NCBI. Trends Genet. 2000, 16: 44-47. 10.1016/S0168-9525(99)01882-X.
Pruitt KD, Tatusova T, Maglott DR: NCBI Reference Sequence (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2005, D501-D504.
Okazaki Y, Furuno M, Kasukawa T, Adachi J, Bono H, Kondo S, Nikaido I, Osato N, Saito R, Suzuki H, et al: Analysis of the mouse transcriptome based on functional annotation of 60,770 full-length cDNAs. Nature. 2002, 420: 563-573. 10.1038/nature01266.
Schuler GD: Pieces of the puzzle: expressed sequence tags and the catalog of human genes. J Mol Med. 1997, 75: 694-698. 10.1007/s001090050155.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410. 10.1006/jmbi.1990.9999.
Hubbard T, Barker D, Birney E, Cameron G, Chen Y, Clark L, Cox T, Cuff J, Curwen V, Down T, et al: The Ensembl genome database project. Nucleic Acids Res. 2002, 30: 38-41. 10.1093/nar/30.1.38.
Large Scale Sequence Annotation System. [http://annotator.imp.univie.ac.at/]
Brendel V, Bucher P, Nourbakhsh IR, Blaisdell BE, Karlin S: Methods and algorithms for statistical analysis of protein sequences. Proc Natl Acad Sci USA. 1992, 89: 2002-2006.
Promponas VJ, Enright AJ, Tsoka S, Kreil DP, Leroy C, Hamodrakas S, Sander C, Ouzounis CA: CAST: an iterative algorithm for the complexity analysis of sequence tracts. Complexity analysis of sequence tracts. Bioinformatics. 2000, 16: 915-922. 10.1093/bioinformatics/16.10.915.
Linding R, Russell RB, Neduva V, Gibson TJ: GlobPlot: Exploring protein sequences for globularity and disorder. Nucleic Acids Res. 2003, 31: 3701-3708. 10.1093/nar/gkg519.
Wootton JC, Federhen S: Analysis of compositionally biased regions in sequence databases. Methods Enzymol. 1996, 266: 554-571.
Sonnhammer EL, Eddy SR, Birney E, Bateman A, Durbin R: Pfam: multiple sequence alignments and HMM-profiles of protein domains. Nucleic Acids Res. 1998, 26: 320-322. 10.1093/nar/26.1.320.
Letunic I, Copley RR, Schmidt S, Ciccarelli FD, Doerks T, Schultz J, Ponting CP, Bork P: SMART 4.0: towards genomic data integration. Nucleic Acids Res. 2004, 32: D142-D144. 10.1093/nar/gkh088.
Sigrist CJ, Cerutti L, Hulo N, Gattiker A, Falquet L, Pagni M, Bairoch A, Bucher P: PROSITE: a documented database using patterns and profiles as motif descriptors. Brief Bioinform. 2002, 3: 265-274.
Marchler-Bauer A, Panchenko AR, Shoemaker BA, Thiessen PA, Geer LY, Bryant SH: CDD: a database of conserved domain alignments with links to domain three-dimensional structure. Nucleic Acids Res. 2002, 30: 281-283. 10.1093/nar/30.1.281.
Schaffer AA, Wolf YI, Ponting CP, Koonin EV, Aravind L, Altschul SF: IMPALA: matching a protein sequence against a collection of PSI-BLAST-constructed position-specific score matrices. Bioinformatics. 1999, 15: 1000-1011. 10.1093/bioinformatics/15.12.1000.
Tusnady GE, Simon I: Principles governing amino acid composition of integral membrane proteins: application to topology prediction. J Mol Biol. 1998, 283: 489-506. 10.1006/jmbi.1998.2107.
von Heijne G: Membrane protein structure prediction. Hydrophobicity analysis and the positive-inside rule. J Mol Biol. 1992, 225: 487-494. 10.1016/0022-2836(92)90934-C.
Cserzo M, Eisenhaber F, Eisenhaber B, Simon I: On filtering false positive transmembrane protein predictions. Protein Eng. 2002, 15: 745-752. 10.1093/protein/15.9.745.
Lupas A, Van Dyke M, Stock J: Predicting coiled coils from protein sequences. Science. 1991, 252: 1162-1164.
Frishman D, Argos P: Incorporation of non-local interactions in protein secondary structure prediction from the amino acid sequence. Protein Eng. 1996, 9: 133-142.
Eisenhaber F, Imperiale F, Argos P, Frommel C: Prediction of secondary structural content of proteins from their amino acid composition alone. I. New analytic vector decomposition methods. Proteins. 1996, 25: 157-168. 10.1002/(SICI)1097-0134(199606)25:2<157::AID-PROT2>3.0.CO;2-F.
Eisenhaber F, Frommel C, Argos P: Prediction of secondary structural content of proteins from their amino acid composition alone. II. The paradox with secondary structural class. Proteins. 1996, 25: 169-179. 10.1002/(SICI)1097-0134(199606)25:2<169::AID-PROT3>3.3.CO;2-5.
von Heijne G: A new method for predicting signal sequence cleavage sites. Nucleic Acids Res. 1986, 14: 4683-4690.
Bendtsen JD, Nielsen H, von Heijne G, Brunak S: Improved prediction of signal peptides: SignalP 3.0. J Mol Biol. 2004, 340: 783-795. 10.1016/j.jmb.2004.05.028.
Eisenhaber B, Eisenhaber F, Maurer-Stroh S, Neuberger G: Prediction of sequence signals for lipid post-translational modifications: insights from case studies. Proteomics. 2004, 4: 1614-1625. 10.1002/pmic.200300781.
Eisenhaber B, Bork P, Eisenhaber F: Prediction of potential GPI-modification sites in proprotein sequences. J Mol Biol. 1999, 292: 741-758. 10.1006/jmbi.1999.3069.
Maurer-Stroh S, Eisenhaber B, Eisenhaber F: N-terminal N-myristoylation of proteins: prediction of substrate proteins from amino acid sequence. J Mol Biol. 2002, 317: 541-557. 10.1006/jmbi.2002.5426.
Blake JA, Richardson JE, Bult CJ, Kadin JA, Eppig JT: MGD: the Mouse Genome Database. Nucleic Acids Res. 2003, 31: 193-195. 10.1093/nar/gkg047.
Halees AS, Leyfer D, Weng Z: PromoSer: A large-scale mammalian promoter and transcription start site identification service. Nucleic Acids Res. 2003, 31: 3554-3559. 10.1093/nar/gkg549.
Quandt K, Frech K, Karas H, Wingender E, Werner T: MatInd and MatInspector: new fast and versatile tools for detection of consensus matches in nucleotide sequence data. Nucleic Acids Res. 1995, 23: 4878-4884.
R Project. [http://www.r-project.org]
Kasprzyk A, Keefe D, Smedley D, London D, Spooner W, Melsopp C, Hammond M, Rocca-Serra P, Cox T, Birney E: EnsMart: a generic system for fast and flexible access to biological data. Genome Res. 2004, 14: 160-169. 10.1101/gr.1645104.
Griffiths-Jones S: The microRNA registry. Nucleic Acids Res. 2004, 32: D109-D111. 10.1093/nar/gkh023.
We thank Dr Fatima Sanchez-Cabo for assistance with the statistical analyses, Bernhard Mlecnik for assistance with the miRNA analysis, Dietmar Rieder for chromosomal mapping, Gernot Stocker for support with the computational infrastructure, Roman Fiedler for RT-PCR analysis, and Dr James McNally for discussions and comments on the manuscript. This work was supported by the Austrian Science Fund, Project SFB Biomembranes F718, the bm:bwk GEN-AU projects Bioinformatics Integration Network (BIN), and Genomics of Lipid-Associated Disorders (GOLD).
Hubert Hackl, Thomas Rainer Burkard contributed equally to this work.
Electronic supplementary material
Additional data file 13: A figure showing a comparison with the data set reported by Soukas and coworkers  (PNG 236 KB)
Additional data file 14: A figure showing a comparison with the data set reported by Ross and coworkers  (PNG 336 KB)
Additional data file 15: A figure showing a comparison with the data set reported by Burton and coworkers  (PNG 348 KB)
Additional data file 18: A figure illustrating the significant over-representation of miRNA motifs in the 3'-UTR of genes in each cluster (PNG 99 KB)
Additional data file 19: A figure showing the significant over-representation of miRNA motifs in the 3'-UTR from genes in cluster 1 (PNG 7 KB)
Additional data file 20: A figure showing the significant over-representation of miRNA motifs in the 3'-UTR from genes in cluster 2 (PNG 17 KB)
Additional data file 21: A figure showing the significant over-representation of miRNA motifs in the 3'-UTR from genes in cluster 3 (PNG 9 KB)
Additional data file 22: A figure showing the significant over-representation of miRNA motifs in the 3'-UTR from genes in cluster 4 (PNG 12 KB)
Additional data file 23: A figure showing the significant over-representation of miRNA motifs in the 3'-UTR from genes in cluster 6 (PNG 13 KB)
Additional data file 24: A figure showing the significant over-representation of miRNA motifs in the 3'-UTR from genes in cluster 7 (PNG 7 KB)
Additional data file 25: A figure showing the significant over-representation of miRNA motifs in the 3'-UTR from genes in cluster 8 (PNG 35 KB)
Additional data file 26: A figure showing the significant over-representation of miRNA motifs in the 3'-UTR from genes in cluster 9 (PNG 39 KB)
Additional data file 27: A figure showing the significant over-representation of miRNA motifs in the 3'-UTR from genes in cluster 10 (PNG 24 KB)
Additional data file 28: A figure showing the significant over-representation of miRNA motifs in the 3'-UTR from genes in cluster 11 (PNG 9 KB)
Additional data file 29: A figure showing the significant over-representation of miRNA motifs in the 3'-UTR from genes in cluster 12 (PNG 27 KB)
Additional data file 30: A table including over-represented miRNA motifs in the 3'-UTR from genes in the set of 780 selected ESTs (XLS 28 KB)
Additional data file 33: Figure showing the cellular localization of gene products involved in metabolism and their gene expression at different time points (PDF 2 MB)
Additional data file 34: Figure showing the cellular localization of gene products involved in other biological processes and their gene expression at different time points (PDF 2 MB)
Additional data file 36: A text file describing extracellular matrix remodeling and cytoskeleton reorganization (PDF 48 KB)
Additional data file 39: A list of experimental verified binding site for PPAR:RXR and the derived position weight matrix (PDF 137 KB)
Additional data file 44: An Excel file showing a statistical analysis of colocalized ESTs for 780 selected ESTs (XLS 22 KB)
Additional data file 45: An Excel file showing a statistical analysis of colocalized ESTs for 5,502 ANOVA selected ESTs (XLS 26 KB)
Authors’ original submitted files for images
About this article
Cite this article
Hackl, H., Burkard, T.R., Sturn, A. et al. Molecular processes during fat cell development revealed by gene expression profiling and functional annotation. Genome Biol 6, R108 (2005). https://doi.org/10.1186/gb-2005-6-13-r108
- Metabolic Network
- Additional Data File
- miRNA Target
- Adipocyte Differentiation
- Extracellular Matrix Remodel