Physiol. Genomics 35: 254-261, 2008.
First published August 26, 2008; doi:10.1152/physiolgenomics.00257.2007
1094-8341/08 $8.00
Received 1 November 2007;
accepted in final form 15 August 2008.
Physiological Genomics 35:254-261 (2008)
1094-8341/08 $8.00 © 2008 American Physiological Society
Response of the adipose tissue transcriptome to dihydrotestosterone in mice
Yonghua Zhang
1,
Ezequiel Calvo
1,
Céline Martel
1,
Van Luu-The
1,
Fernand Labrie
1 and
André Tchernof
1,2
1 Molecular Endocrinology and Oncology Research Center, Laval University Medical Research Center, Québec, Canada
2 Department of Food Science and Nutrition, Laval University, Québec, Canada
 |
ABSTRACT
|
|---|
Androgens have been postulated to be important modulators of adipose tissue metabolism and fat cell function. In the present study, we investigated the response of male and female mice retroperitoneal adipose tissue to the nonaromatizable androgen dihydrotestosterone (DHT). Adipose tissue samples were obtained in gonadectomized animals treated with vehicle (control group), or injected with 0.1 mg DHT 1, 3, 6, 12, 18, and 24 h prior to necropsy. Fourteen animals were pooled at each time point (total 196 animals). Transcripts that were significantly modulated were considered as androgen-responsive genes. Quantitative real-time RT-PCR was used to confirm results from the microarray analysis in a subset of 46 probe sets in male mice and 98 probe sets in female mice. Considering peak time vs. control, we confirmed 74.0 and 63.3% of the modulated genes by PCR in males and females, respectively. Four genes were significantly stimulated in a similar manner by DHT in both sexes, namely metallothionein 1, growth arrest and DNA-damage-inducible 45 gamma, cyclin-dependent kinase inhibitor 1A, and fk506-binding protein 5. All these genes appear to be involved in the regulation of adipocyte differentiation/proliferation and adipogenesis. In conclusion, this study, which evaluated the acute transcriptome response of adipose tissue to DHT in male and female mice, suggests that DHT consistently modulates genes involved in the regulation of adipogenesis in retroperitoneal adipose tissue of both male and female animals.
retroperitoneal; microarray; adipogenesis
 |
INTRODUCTION
|
|---|
SEVERAL STUDIES SUGGEST THAT sex hormones modulate adipose tissue deposition at the level of preadipocyte differentiation and/or proliferation as well as adipocyte lipogenesis and/or lipolysis (2, 36, 40). Cross-sectional (35) and longitudinal (18, 22) studies found that circulating androgen levels are negatively associated with abdominal obesity and visceral fat accumulation in men. A number of studies also reported negative associations between circulating androgens and abdominal adiposity in women (3, 10, 47). On the other hand, female-to-male transsexuals treated with high doses of androgens preferentially accumulate intra-abdominal fat and progressively lose gluteal-femoral fat (15). In women with the polycystic ovary syndrome, a positive association between abdominal obesity and plasma androgen levels is found (13). Triglycerides in retroperitoneal and mesenteric adipose tissue increase after castration in rats but decrease when animals are treated with testosterone (25). Additionally, androgen receptor knockout mice develop late onset obesity through reduced adipose tissue lipolysis in animals of both sexes (16, 42). These findings suggest that androgens may be key steroids in the modulation of fat accumulation and distribution in both humans and rodent models.
The impact of androgens on gene transcription in adipose tissue from male mice has been previously examined by serial analysis of gene expression in one study (7). Although very informative, that approach did not consider temporal expression regulation in response to the hormone. In the present study, using DNA microarrays, we examined the acute response of the adipose tissue transcriptome to dihydrotestosterone (DHT) in gonadectomized male and female mice. We tested the general hypothesis that transcriptional profiling of androgen-responsive genes in both male and female mice would bring new insights on possible signaling pathways linking fat distribution and sex hormones.
 |
METHODS
|
|---|
Mice and RNA sample preparation.
Male and female, 14- to 15-wk-old, C57BL/6 mice were obtained from Charles River Canada (St-Constant, Québec, Canada). Fourteen animals were pooled for one replica at each time point (1, 3, 6, 12, 18, and 24 h), and fourteen mice were pooled to generate two replica of the control group in each sex. A total of 196 animals were used in our study. Mice were housed (1 per cage) in plastic cages under 12 h light-dark cycles; the animals had access to Lab Rodent Diet no. 5002 and were given tap water ad libitum. Bilateral gonadectomy (castration for male mice and ovariectomy for female mice) were performed under isoflurane-induced anesthesia 7 days prior to death and organ collection for the seven gonadectomized groups (n = 14 mice per group). DHT (0.1 mg) was injected subcutaneously 1, 3, 6, 12, 18, and 24 h prior to death in DHT-treated groups. One group was retained as control and received a vehicle solution [0.4% (wt/vol) Methocel A15LV Premium; 5% ethanol]. Retroperitoneal adipose tissue samples from all mice of the same group were collected, rapidly trimmed, snap-frozen in liquid nitrogen and stored at –80°C prior to mRNA extraction. These fat samples were pooled to eliminate interindividual variation and to extract a sufficient amount of mRNA. For blood DHT measurements, gonadectomized C57BL/6 male and female mice were injected with 0.1 mg DHT 1, 3, 6, 12, 18, or 24 h prior to blood sampling (n = 12 mice/time point/sex). Control mice (n = 24 mice/sex) were injected with the vehicle alone (5% ethanol-0.4% methylcellulose) 24 h prior to necropsy. Mice under isoflurane anesthesia were exsanguinated by cardiac venipuncture. Blood samples were processed for serum preparation and kept frozen at –80°C until use for the determination of serum DHT concentrations by gas chromatography mass spectrometry as previously described (23) with the exception of the serum volume used (0.15–0.3 ml in mice instead of 0.75 ml in humans). Serum samples from two mice of the same group were pooled for each determination of DHT levels (n = 6 for DHT groups; n = 12 for the control group). The experiment was conducted in an animal facility approved by the Canadian Council on Animal Care (CCAC) and the Association for Assessment and Accreditation of Laboratory Animal Care. The study was performed in accordance with the CCAC Guide for Care and Use of Experimental Animals. Institutional approval was obtained.
Target preparation.
Total mRNA was isolated with TRIzol (Invitrogen, Burlington, Ontario) following the manufacturer's protocol. Total mRNA (20 µg) was converted to cDNA by incubation with 400 U SuperScript II reverse transcriptase (Invitrogen) using a T7- oligo-(dt)24 primer [5'-GGCCAGTGAATTGTAATACGACTCACTATAGGGAGGCGG-(dt)24-3'] 1x first-strand buffer (50 mM Tris·HCl pH 8.3, 75 mM KCl, 3 mM MgCl2, 10 mM DTT) and 0.5 mM dNTPs at 42°C for 1 h. Second-strand synthesis was performed using 40 U DNA polymerase I (Invitrogen), 10 U Escherichia coli DNA ligase (Invitrogen), 2 U RNase H (Invitrogen), 1x reaction buffer [18.8 mM Tris·HCl pH 8.3, 90.6 mM KCl, 4.6 mM MgCl2, 3.8 mM DTT, 0.15 mM NAD, 10 mM (NH4)2SO4] and 0.2 mM dNTPs at 16°C for 2 h. cDNAs were blunt-ended with 10 U of T4 DNA polymerase (Invitrogen) incubated 5 min at 16°C. cDNA was extracted with phenol-chloroform using phase lock gels (Brinkman, Mississauga, Ontario, Canada), ethanol precipitated, and resuspended in 10 µl of DEPC-treated H2O. cDNA was in vitro transcribed using a T7 BioArray High Yield RNA Transcript Labeling Kit (Enzo Diagnostics, Farmingdale, NY) to produce biotinylated cRNA. The mixture (20 µl final volume) was incubated at 37°C for 5 h, with gentle mixing every 30 min. Labeled cRNA was purified with an RNeasy Mini Kit (Qiagen, Valencia, CA) according to the manufacturer's protocol. Purified cRNA was fragmented to 30–200 mer using a fragmentation buffer (100 mM potassium acetate, 30 mM magnesium acetate, 40 mM Tris-acetate pH 8.1), for 20 min at 94°C. The quality of total RNA, cDNA synthesis, cRNA amplification and cRNA fragmentation was monitored by microcapillary electrophoresis (Bioanalyzer 2100; Agilent Technologies, Palo Alto, CA).
Microarray hybridization and scanning.
Microarray analysis was performed using Murine GeneChip U74Av2 GeneChips (Affymetrix, Santa Clara, CA), which allow the analysis of the expression level of >12,400 mouse genes and expressed sequence tags (ESTs) on each array in the present study. Fragmented cRNA (15 µg) was incubated with 1x hybridization buffer (0.1 mg/ml herring sperm DNA, 0.5 mg/ml acetylated BSA, 5 nM control oligonucleotide B2) and 1x eukaryotic hybridization control solution (1.5 pM BioB, 5 pM BioC, 25 pm BioD, and 100 pM cre), for 16 h at 45°C with constant rotation (60 rpm). Microarrays were processed using an Affymetrix GeneChip Fluidic Station 400 (protocol EukGE-WS2Av4). Staining was performed with streptavidin-conjugated phycoerythrin (SAPE) followed by amplification with a biotinylated anti-streptavidin antibody and by a second round of SAPE solution. GeneChips were scanned using an Agilent GeneArray Scanner (Agilent Technologies). Signal intensities for β-actin and GAPDH genes were used as internal quality controls. The ratio of fluorescent intensities for the 5' and 3' ends of these housekeeping genes was <2. Scanned images were analyzed with Microarray Suite 5.0 (MAS, Affymetrix). Tissues from the 14 animals at each time point were pooled, and the cRNA probe corresponding to each time point was hybridized once on separate microarrays. Tissues from the 14 animals of the control group were pooled, and the cRNA probe from the control group was hybridized on two microarrays (duplicate). A total of 16 microarrays were used in our experiments. The microarray data have been deposited in the Gene Expression Omnibus repository (series number GSE9631).
Real-time RT-PCR.
First-strand cDNA synthesis was performed using 5 µg of RNA extracted from each time pool in a reaction containing 200 units of Superscript II RNase H-reverse transcriptase (Invitrogen), 300 ng of oligo dT18, 500 µM dNTP, 10 mM DTT, and 34 units of human RNase inhibitor (Amersham Pharmacia) in a final volume of 50 µl. The resulting products were then treated with 1 µg of RNase A for 30 min at 37°C and purified thereafter with Qiaquick PCR purification kits (Qiagen). For quantitative PCR analyses, a Light-Cycler PCR (Roche Diagnostics) was used to confirm the expression level of cDNAs measured in the microarray method. The FastStart DNA Master SYBR green kit (Roche Diagnostics) was used in a final reaction volume of 20 µl containing 3 mM MgCl2 and cDNA corresponding to 20 ng total RNA. The PCR was carried out according to the following conditions: 95°C/10 min, 50 cycles of (95°C/10 s, 55–65°C/5 s, 72°C/11 s) with a temperature transition of 3°C/s. PCR results were calculated and normalized using the second-derivative and double-correction method as described previously (29), and with use of the housekeeping gene ATP synthase O subunit (Atp5o). This gene was found to display remarkably stable expression levels from embryonic life through adulthood in various tissues (48). Messenger RNA levels were expressed as the number of copies per microgram of total RNA using a standard curve of crossing points vs. logarithm of quantity. The standard curve was established using known cDNA amounts of 10, 102, 103, 104, 105, and 106 copies of Atp5o. The Light Cycler 3.5 program was provided by the manufacturer (Roche). Expression of the housekeeping gene was constant at all time points and not significantly influenced by DHT treatment.
Data processing and statistical analysis.
Raw data of microarrays were collected and analyzed using the Affymetrix MAS 5.0 software. To control for variability of the measurement, which is greater at low expression level (50), we established a relationship between fold changes and expression levels as described previously (31). To select differentially expressed genes, we used the following criteria: 1) the fold change (FC) had to be higher than the limit fold change (LFC) over the time course; and 2) at least two "present calls" (P value <0.04) had to be obtained across all arrays. The LFC was defined as LFC = a + b/X, where X is the minimum average difference intensity (ADI, i.e., absolute expression level) for a given gene. The parameters a and b were estimated based on the distribution of ratios calculated from replicated chips (1.9 and 60, respectively, in the present experiment). The FC was calculated with the following equation: FC = Max ADI/Min ADI over the time course and control. Average ADI of two technique replica of control were used in data analysis. To identify genes to be validated by quantitative real-time RT-PCR, genes were further filtered following additional criteria: 1) ADI of a given gene >400 on >60% of arrays; 2) absolute FC on peak time (peak time ADI/control ADI)
2.0. Annotations of interesting probe sets were collected from different available databases using the NetAffx online tool (http://www.affymetrix.com/analysis/index.affx), and their identity was verified with BLAST (http://www.ncbi.nlm.nih.gov). The Ingenuity Pathway Analysis System (Ingenuity Systems, www.ingenuity.com) was also used to visualize gene expression data in the context of biological pathways.
Validation analyses were performed using peak time vs. control comparisons to examine the direction of fold changes with the microarray and RT-PCR methods. Since data for all time points were also compared with those of real-time RT-PCR for the genes selected, Spearman rank nonparametric correlation analysis was also used to compare time courses obtained with microarrays and real time RT-PCR. Time points from the microarray analysis were correlated to time points from the real-time RT-PCR analysis and DHT responses that generated correlation coefficients with a P
0.1 were considered similar. This permissive cutoff was chosen to avoid underestimation of time course reproducibility.
 |
RESULTS
|
|---|
In both male and female mice, serum DHT reached the maximal level 1 h after treatment and decreased thereafter to reach baseline at 24 h (Fig. 1). The overall objective of the data mining approach was the identification of candidate genes within adipose tissue that were sensitive to DHT treatment. Globally, based on the criteria mentioned in METHODS, we found that 94 probe sets were significantly modulated by DHT on chip arrays from male mice representing 93 genes (Supplemental Table S1).1
A total of 116 probe sets were significantly modulated by DHT on chip arrays from female mice representing 112 genes (Supplemental Table S2). These genes exhibited gene expression changes that were
2.0-fold comparing peak time point to control.

View larger version (14K):
[in this window]
[in a new window]
|
Fig. 1. Serum dihydrotestosterone (DHT) was measured at each time point after subcutaneous injection with 0.1 mg DHT in each male and female mouse as well as each of control gonadectomized (GDX) animals. Data are presented as means ± SE.
|
|
The distributions of peaking times for significantly modulated probe sets were examined (Fig. 2). In male mice, 38 (40%) probe sets presented peak times at 1 h. A total of 49 (42.2%) and 20 (17.2%) probe sets in female mice presented peak times at 18 and 24 h, respectively. In male mice, 33 (35.1%) upregulated probe sets and 61 (64.9%) downregulated probe sets were identified. In female mice, 51 (44%) probe sets were upregulated and 65 (56%) probe sets were downregulated.

View larger version (19K):
[in this window]
[in a new window]
|
Fig. 2. Distribution of peaking times for probe sets that were significantly modulated in male and female retroperitoneal adipose tissue (AT) samples.
|
|
Among the microarray data, a total of 21 genes were found to be significantly modulated in adipose tissue samples from both male and female mice. Sixteen of them changed in a similar direction in response to DHT. Five genes showed opposite regulation trends in males and females, including LOC677213, Tnnc2, Ddx6, Mgea5, and Myh1. Four of the latter genes were downregulated in male mice and upregulated in female mice. Only the Mgea5 gene was upregulated in male mice and downregulated in female mice (data not shown).
To confirm gene expression responses by an independent method, quantitative real-time RT-PCR validation was conducted. We validated subsets of 46 and 98 transcripts in male and female mice, respectively (Supplemental Tables S1 and S2). Using peak time vs. control comparisons, we confirmed 74.0 and 63.3% of modulated genes by PCR in male and female mice, respectively. We also used a more stringent approach to precisely compare time courses between the two methods. Using P
0.1 as a cut-off for significance, we performed Spearman nonparametric correlation analysis, and 17 genes (37.0%) showed reproducible time courses in PCR validation in male mice (Table 1). In female mice, 26 genes (26.5%) showed a Spearman rank correlation coefficient with a P value <0.10, indicating identical time courses in PCR validation (Table 2). Among upregulated genes with positive PCR validations in male mice, five were related to cell cycle; four to DNA replication, recombination, and repair; 10 to cell death; seven to cellular development; and eight to cellular growth and proliferation. Among downregulated genes, five were related to carbohydrate metabolism, three to cellular compromise, two to protein degradation, four to cell death, and three to cell-to-cell signaling and interaction. In female mice, among upregulated genes with positive PCR validations, nine were related to cellular development, eight to cellular growth and proliferation, eight to cell death, seven to cellular assembly and organization, and four to cellular compromise. Among downregulated genes, seven were related to lipid metabolism, nine to molecular transport, three to nucleic acid metabolism, 12 to small molecule biochemistry, and six to cell cycle.
Table 3 shows that four genes generated similar responses to the DHT treatment in both sexes, namely FK506 binding protein 5 (Fkbp5), cyclin-dependent kinase inhibitor 1A (p21 or Cdkn1a), growth arrest and DNA-damage-inducible 45 gamma (Gadd45g), and metallothionein 1 (Mt1). The DHT response of all these genes was validated by quantitative real-time RT-PCR. They all showed highly significant Spearman rank correlation coefficients with P values of 0.0001, 0.014, 0.003, and 0.023 in fat samples from male mice and 0.0001, 0.01, 0.07, and 0.0001 in fat samples from female mice. These genes were all upregulated by DHT in both sexes (Table 3 and Fig. 3).

View larger version (20K):
[in this window]
[in a new window]
|
Fig. 3. Effect of DHT after 1, 3, 6, 12, 18, or 24 h of treatment or vehicle (GDX) on mRNA expression levels of 4 androgen-responsive genes in retroperitoneal AT from both male and female mice. Microarray results (solid line) were validated by quantitative real-time RT-PCR (dotted line). Data are shown in relative expression units for the results from the microarray analysis and in relative mRNA abundance for the results of RT-PCR validation.
|
|
 |
DISCUSSION
|
|---|
This work represents the first microarray investigation on the effects of DHT treatment in retroperitoneal adipose tissue from a rodent model. Serum DHT measures obtained over the course of the experiment suggest that our study should be considered as an acute response experiment. A significant modulation of the adipose tissue transcriptome was observed in response to DHT in both male and female mice. A number of genes involved in distinct cellular functions were modulated in a sex-specific manner. A total of 21 genes were found to be simultaneously modulated in fat samples from both male and female mice. Sixteen of them showed similar regulatory patterns, and five genes showed opposite regulatory patterns.
Bolduc et al. (7) previously examined the effects of DHT on adipose tissue in male mice using serial analysis of gene expression at two time points (3 and 24 h). They found that various pathways were regulated by DHT in retroperitoneal adipose tissue. Of note, DHT treatment indicated a stimulation of genes involved in fatty acid and triacylglycerol production as well as lipolysis in fat samples. They concluded that almost all aspects of cell function were affected by DHT through the modulation of gene expression. Much like the latter study, several cell functions were affected by DHT in the present analysis. For example, in male mice, we found that some genes related to cell cycle, DNA replication and recombination, as well as repair, molecular transport, and cell death were upregulated. Other genes involved in cellular function and maintenance were downregulated. In female mice, some genes related to connective tissue development and function or tissue morphology were upregulated, and some genes involved in cell cycle and cell death were downregulated.
Only four genes showed identical regulatory patterns in both male and female mice. These genes were Fkbp5, Cdkn1a, Gadd45g, and Mt1. All these genes, when upregulated, appear to be related to a negative regulation of adipocyte differentiation and adipogenesis. Fkbp5 is considered as a glucocorticoid receptor heat shock protein 90-associated cochaperone and plays a central role in steroid hormone signaling (33). Fkbp51 overexpression induces glucocorticoid resistance in squirrel monkeys by reducing hormone-binding affinity (38). Our finding indirectly suggests that upregulated expression of Fkbp5 in response to DHT treatment may mediate inhibitory effects of DHT on glucocorticoid-induced fat cell differentiation and lipogenesis. Androgens stimulate lipolysis in adipose tissue and induce an antiadipogenic effect, at least in primary cultured preadipocytes, in adipocyte precursor cells and in 3T3-L1 preadipocytes (11, 43, 44).
The Gadd45 family of genes is involved in cell cycle and programmed cell death. Specifically, proteins of the Gadd45 family play a role in the inhibition of cellular growth (17, 49) and promote apoptotic cell death (19, 21, 28, 45). They are powerful suppressors of cell growth (20). Testosterone and DHT are thought to inhibit adipocyte differentiation through an AR-mediated nuclear translocation of β-catenin and activation of downstream Wnt signaling (43). Studies in a different model have shown that Gadd45a plays a role in the negative regulation of cell development through the suppression of β-catenin (20). Elevated expression of the Gadd45g gene in the present study may relate to reduced adipogenesis in response to DHT through this pathway.
Regarding MT1, obesity is induced in older MT-null transgenic mice (5), which are characterized by increases in body weight and white adipose tissue that are similar to that of obese diabetic (db/db) mice. Under cold exposure and induction by catecholamines, the MT1 gene is abundantly detected in brown adipose tissue of rats (4, 6). Secretion or expression of MT protein is also observed in white adipose tissue from mice (46), dogs (14), and humans (12). All these observations indicate that expression of MT genes is linked to the regulation of energy balance. Our data on increased MT1 in response to DHT relate androgens to adipose tissue metabolism and energy homeostasis.
As an adipogenic transcription factor, PPAR
induces growth arrest and is involved in adipocyte differentiation through cyclin-dependent kinase inhibitors p18 and p21 (30). p27/p21 double knockout mice become obese and show adipocyte hyperplasia and related metabolic consequences (32). Significantly upregulated Cdkn1a expression was detected in the present study, which suggests that DHT treatment may decrease adipocyte number and lead to reductions in retroperitoneal adiposity in male and female mice. We hereby provide evidence indicating that DHT signaling involves cross talks between adipocyte proliferation, apoptosis, and differentiation pathways as well as energy homeostasis in retroperitoneal fat. These phenomena appear to be independent of the sex of the animal.
We identified five genes that were simultaneously and significantly modulated in fat samples from both male and female mice but showed opposite regulatory patterns. LOC677213, Tnnc2, Ddx6, and Myh1 were downregulated in male mice and upregulated in female mice. These genes are involved in the differentiation of skeletal muscle cells (9, 34) or immunity (1). Mgea5 was upregulated in male mice and downregulated in female mice. The meningioma-expressed antigen 5 (hyaluronidase) gene (Mgea5) encodes a β-O-linked N-acetylglucosaminidase (O-GlcNAcase) and is involved in Type 2 diabetes as a biological and positional candidate gene (8).
Several other genes that were significantly modulated in the microarray analysis and positively validated by quantitative real-time RT-PCR in both male and female retroperitoneal fat samples were myogenic in nature. These genes included myosin light polypeptide kinase (Mylk), myosin heavy polypeptide 4 (Myh4), enolase 3 (Eno3), phosphofructokinase (Pfkm), and myosin heavy polypeptide 1 (Myh1). The procedure for fat sampling was carefully reviewed to ascertain that the sampling included only retroperitoneal adipose tissue and that the tissues were carefully trimmed. Similar findings of myogenic genes in fat were also reported in other studies (24, 41, 51). Adipose tissue apparently shows mesenchymal plasticity for myogenic differentiation (24), and adipocyte precursors and skeletal myoblasts are capable of differentiating along similar developmental lineages (39). The use of whole tissue samples prevents us from reaching conclusions on which cell type is responsible for the androgenic response observed and should be acknowledged. Nevertheless, Lin et al. (26) found that fat pad weight and total lipid content of myostatin-null mice were significantly lowered compared with wild-type mice and suggested that increased muscle development in myostatin knockout mice resulted in reduced adipogenesis and leptin secretion. Rebbapragada et al. (37) also found that BMP7-induced adipogenesis is restrained by myostatin in both mesenchymal precursor cells and preadipocytes. These previous reports and the present findings support the notion that a strong relationship between myogenic and adipose genes may exist in adipose tissue of both male and female mice.
In this study, we used six time points to examine trends of gene modulation in mice from both sexes. Significant modulations were observed with microarray data using criteria described in METHODS. With respect to validation by RT-PCR, most previous studies used only two time points on a limited number of genes. In the present study, a very large number of genes were validated by RT-PCR on all six time points and in the vehicle group. Using a rather stringent validation procedure, our analysis indicates that microarray results are moderately reproducible compared with quantitative real-time RT-PCR over several time points. Indeed, only 37.0 and 26.5% of transcripts in male and female mice respectively showed identical time courses. Using only peak time vs. control comparisons to validate the direction of fold changes in gene expression, valid findings increased to 74.0 and 63.3% of modulated genes in male and female mice, respectively. Because we used the pooled samples, our analysis may have contributed to yield lower than expected reproducibility rates. The use of several arrays would have allowed us to use repeated-measures ANOVA to identify significantly modulated genes. This may have led to higher reproducibility rates when comparing with RT-PCR results. We suggest, however, that the use of multiple time points and a rather stringent filtering procedure decreased the likelihood of identified false positives. Our results nevertheless suggest that careful validation on a different platform remains an absolute necessity (27).
In conclusion, our study indicates that microarray technology can be used to identify androgen-responsive genes in retroperitoneal adipose tissue from male and female mice. We found that DHT consistently stimulates genes that are associated with the regulation of adipogenesis in adipose tissue. This study also provides evidence to support the hypothesis that a myogenic response may be triggered by DHT in fat.
 |
GRANTS
|
|---|
André Tchernof is funded by the Canadian Institutes of Health Research. This work was supported by Genome Quebec and Genome Canada.
 |
ACKNOWLEDGMENTS
|
|---|
We thank all the research assistants involved in the ATLAS project for skilful technical assistance. We thank Karine Petit and Sophie Blanchette for contributions to the initial analysis of microarray data.
 |
FOOTNOTES
|
|---|
Address for reprint requests and other correspondence: A. Tchernof, Molecular Endocrinology and Oncology Research Center, Dept. of Food Science and Nutrition, Laval Univ. Medical Research Center, 2705 Laurier Blvd. (T3-67), Québec, (Québec), Canada G1V 4G2 (e-mail: andre.tchernof{at}crchul.ulaval.ca).
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
1 The online version of this article contains supplemental material. 
 |
REFERENCES
|
|---|
- Akao Y, Matsuda Y. Identification and chromosome mapping of the mouse homologue of the human gene (DDX6) that encodes a putative RNA helicase of the DEAD box protein family. Cytogenet Cell Genet 75: 38–44, 1996.[CrossRef][Web of Science][Medline]
- Anderson LA, McTernan PG, Barnett AH, Kumar S. The effects of androgens and estrogens on preadipocyte proliferation in human adipose tissue: influence of gender and site. J Clin Endocrinol Metab 86: 5045–5051, 2001.[Abstract/Free Full Text]
- Armellini F, Zamboni M, Castelli S, Robbi R, Mino A, Todesco T, Bergamo-Andreis IA, Bosello O. Interrelationships between intraabdominal fat and total serum testosterone levels in obese women. Metabolism 43: 390–395, 1994.[CrossRef][Web of Science][Medline]
- Beattie JH, Black DJ, Wood AM, Trayhurn P. Cold-induced expression of the metallothionein-1 gene in brown adipose tissue of rats. Am J Physiol Regul Integr Comp Physiol 270: R971–R977, 1996.[Abstract/Free Full Text]
- Beattie JH, Wood AM, Newman AM, Bremner I, Choo KH, Michalska AE, Duncan JS, Trayhurn P. Obesity and hyperleptinemia in metallothionein (-I and -II) null mice. Proc Natl Acad Sci USA 95: 358–363, 1998.[Abstract/Free Full Text]
- Beattie JH, Wood AM, Trayhurn P, Jasani B, Vincent A, McCormack G, West AK. Metallothionein is expressed in adipocytes of brown fat and is induced by catecholamines and zinc. Am J Physiol Regul Integr Comp Physiol 278: R1082–R1089, 2000.[Abstract/Free Full Text]
- Bolduc C, Larose M, Yoshioka M, Ye P, Belleau P, Labrie C, Morissette J, Raymond V, Labrie F, St-Amand J. Effects of dihydrotestosterone on adipose tissue measured by serial analysis of gene expression. J Mol Endocrinol 33: 429–444, 2004.[Abstract/Free Full Text]
- Cameron EA, Martinez-Marignac VL, Chan A, Valladares A, Simmonds LV, Wacher N, Kumate J, McKeigue P, Shriver MD, Kittles R, Cruz M, Parra EJ. MGEA5-14 polymorphism and type 2 diabetes in Mexico City. Am J Hum Biol 19: 593–596, 2007.[CrossRef][Web of Science][Medline]
- Davis RL, Weintraub H, Lassar AB. Expression of a single transfected cDNA converts fibroblasts to myoblasts. Cell 51: 987–1000, 1987.[CrossRef][Web of Science][Medline]
- De Pergola G, Triggiani V, Giorgino F, Cospite MR, Garruti G, Cignarelli M, Guastamacchia E, Giorgino R. The free testosterone to dehydroepiandrosterone sulphate molar ratio as a marker of visceral fat accumulation in premenopausal obese women. Int J Obes Relat Metab Disord 18: 659–664, 1994.[Web of Science][Medline]
- De Pergola G, Xu XF, Xang SM, Giorgino R, Bjorntorp P. Up-regulation of androgen receptor binding in male rat fat pad adipose precursor cells exposed to testosterone: study in a whole cell assay system. J Steroid Biochem Mol Biol 37: 553–558, 1990.[CrossRef][Web of Science][Medline]
- Do MS, Nam SY, Hong SE, Kim KW, Duncan JS, Beattie JH, Trayhurn P. Metallothionein gene expression in human adipose tissue from lean and obese subjects. Horm Metab Res 34: 348–351, 2002.[CrossRef][Web of Science][Medline]
- Dunaif A. Insulin resistance and the polycystic ovary syndrome: mechanism and implications for pathogenesis. Endocr Rev 18: 774–800, 1997.[Abstract/Free Full Text]
- Eisele I, Wood IS, German AJ, Hunter L, Trayhum P. Adipokine gene expression in dog adipose tissues and dog white adipocytes differentiated in primary culture. Horm Metab Res 37: 474–481, 2005.[CrossRef][Web of Science][Medline]
- Elbers JM, Asscheman H, Seidell JC, Megens JA, Gooren LJ. Long-term testosterone administration increases visceral fat in female to male transsexuals. J Clin Endocrinol Metab 82: 2044–2047, 1997.[Abstract/Free Full Text]
- Fan W, Yanase T, Nomura M, Okabe T, Goto K, Sato T, Kawano H, Kato S, Nawata H. Androgen receptor null male mice develop late-onset obesity caused by decreased energy expenditure and lipolytic activity but show normal insulin sensitivity with high adiponectin secretion. Diabetes 54: 1000–1008, 2005.[Abstract/Free Full Text]
- Fornace AJ Jr, Alamo I Jr, Hollander MC. DNA damage-inducible transcripts in mammalian cells. Proc Natl Acad Sci USA 85: 8800–8804, 1988.[Abstract/Free Full Text]
- Gapstur SM, Gann PH, Kopp P, Colangelo L, Longcope C, Liu K. Serum androgen concentrations in young men: a longitudinal analysis of associations with age, obesity, and race. The CARDIA male hormone study. Cancer Epidemiol Biomarkers Prev 11: 1041–1047, 2002.[Abstract/Free Full Text]
- Harkin DP, Bean JM, Miklos D, Song YH, Truong VB, Englert C, Christians FC, Ellisen LW, Msheswaran S, Oliner JD, Haber DA. Induction of GADD45 and JNK/SAPK-dependent apoptosis following inducible expression of BRCA1. Cell 97: 575–585, 1999.[CrossRef][Web of Science][Medline]
- Hildesheim J, Belova GI, Tyner SD, Zhou X, Vardanian L, Fornace AJ Jr. Gadd45a regulates matrix metalloproteinases by suppressing DeltaNp63alpha and beta-catenin via p38 MAP kinase and APC complex activation. Oncogene 23: 1829–1837, 2004.[CrossRef][Web of Science][Medline]
- Hildesheim J, Bulavin DV, Anver MR, Alvord WG, Hollander MC, Vardanian L, Fornace AJ Jr. Gadd45a protects against UV irradiation-induced skin tumors, and promotes apoptosis and stress signaling via MAPK and p53. Cancer Res 62: 7305–7315, 2002.[Abstract/Free Full Text]
- Khaw KT, Barrett-Connor E. Lower endogenous androgens predict central adiposity in men. Ann Epidemiol 2: 675–682, 1992.[Medline]
- Labrie F, Belanger A, Belanger P, Berube R, Martel C, Cusan L, Gomez J, Candas B, Castiel I, Chaussade V, Deloche C, Leclaire J. Androgen glucuronides, instead of testosterone, as the new markers of androgenic activity in women. J Steroid Biochem Mol Biol 99: 182–188, 2006.[CrossRef][Web of Science][Medline]
- Lee GH, Kemp DM. Human adipose-derived stem cells display myogenic potential and perturbed function in hypoxic conditions. Biochem Biophys Res Commun 341: 882–888, 2006.[CrossRef][Web of Science][Medline]
- Li M, Bjorntorp P. Effects of testosterone on triglyceride uptake and mobilization in different adipose tissues in male rats in vivo. Obes Res 3: 113–119, 1995.[Web of Science][Medline]
- Lin J, Armold HB, Della-Fera MA, Azain MJ, Hartzell DL, Baile CA. Myostatin knockout in mice increases myogenesis and decreases adipogenesis. Biochem Biophys Res Commun 291: 701–706, 2002.[CrossRef][Web of Science][Medline]
- Liu YZ, Dvornyk V, Lu Y, Shen H, Lappe JM, Recker RR, Deng HW. A novel pathophysiological mechanism for osteoporosis suggested by an in vivo gene expression study of circulating monocytes. J Biol Chem 280: 29011–29016, 2005.[Abstract/Free Full Text]
- Lu B, Yu H, Chow C, Li B, Zheng W, Davis RJ, Flavell RA. GADD45gamma mediates the activation of the p38 and JNK MAP kinase pathways and cytokine production in effector TH1 cells. Immunity 14: 583–590, 2001.[CrossRef][Web of Science][Medline]
- Luu-The V, Paquet N, Calvo E, Cumps J. Improved real-time RT-PCR method for high-throughput measurements using second derivative calculation and double correction. Biotechniques 38: 287–293, 2005.[Web of Science][Medline]
- Morrison RF, Farmer SR. Role of PPARgamma in regulating a cascade expression of cyclin-dependent kinase inhibitors, p18(INK4c) and p21(Waf1/Cip1), during adipogenesis. J Biol Chem 274: 17088–17097, 1999.[Abstract/Free Full Text]
- Mutch DM, Berger A, Mansourian R, Tytz A, Roberts MA. The limit fold change model: a practical approach for selecting differentially expressed genes from microarray data. BMC Bioinformatics 3: 17, 2002.[CrossRef][Medline]
- Naaz A, Holsberger DR, Iwamoto GA, Nelson A, Kiyokawa H, Cooke PS. Loss of cyclin-dependent kinase inhibitors produces adipocyte hyperplasia and obesity. FASEB J 18: 1925–1927, 2004.[Abstract/Free Full Text]
- Nuber UA, Kriaucionis S, Roloff TC, Guy J, Selfridge J, Steinhoff C, Schulz R, Lipkowitz B, Ropers HH, Holmes MC, Bird A. Up-regulation of glucocorticoid-regulated genes in a mouse model of Rett syndrome. Hum Mol Genet 14: 2247–2256, 2005.[Abstract/Free Full Text]
- Parmacek MS, Leiden JM. Structure, function, and regulation of troponin C. Circulation 84: 991–1003, 1991.[Free Full Text]
- Pasquali R, Casimirri F, Cantobelli S, Melchionda N, Morselli Labate AM, Fabbri R, Capelli M, Bortoluzzi L. Effect of obesity and body fat distribution on sex hormones and insulin in men. Metabolism 40: 101–104, 1991.[CrossRef][Web of Science][Medline]
- Price TM, O'Brien SN, Welter BH, George R, Anandjiwala J, Kilgore M. Estrogen regulation of adipose tissue lipoprotein lipase–possible mechanism of body fat distribution. Am J Obstet Gynecol 178: 101–107, 1998.[CrossRef][Web of Science][Medline]
- Rebbapragada A, Benchabane H, Wrana JL, Celeste AJ, Attisano L. Myostatin signals through a transforming growth factor beta-like signaling pathway to block adipogenesis. Mol Biol Cell 23: 7230–7242, 2003.[CrossRef]
- Reynolds PD, Ruan Y, Smith DF, Scammell JG. Glucocorticoid resistance in the squirrel monkey is associated with overexpression of the immunophilin FKBP51. J Clin Endocrinol Metab 84: 663–669, 1999.[Abstract/Free Full Text]
- Rink C, Roy S, Khanna S, Rink T, Bagchi D, Sen CK. Transcriptome of the subcutaneous adipose tissue in response to oral supplementation of type 2 Leprdb obese diabetic mice with niacin-bound chromium. Physiol Genomics 27: 370–379, 2006.[Abstract/Free Full Text]
- Roncari DA, Van RL. Promotion of human adipocyte precursor replication by 17beta-estradiol in culture. J Clin Invest 62: 503–508, 1978.[Web of Science][Medline]
- Sasao N, Hirayama E, Kim J. Characterization of heterokaryons between skeletal myoblasts and preadipocytes: myogenic potential of 3T3-L1 preadipocytes. Eur J Cell Biol 82: 97–103, 2003.[CrossRef][Web of Science][Medline]
- Sato T, Matsumoto T, Yamada T, Watanabe T, Kawano H, Kato S. Late onset of obesity in male androgen receptor-deficient (AR KO) mice. Biophys Res Commun 300: 167–171, 2003.[CrossRef][Web of Science][Medline]
- Singh R, Artaza JN, Taylor WE, Braga M, Yuan X, Gonzalez-Cadavid NF, Bhasin S. Testosterone inhibits adipogenic differentiation in 3T3-L1 cells: nuclear translocation of androgen receptor complex with beta-catenin and T-cell factor 4 may bypass canonical Wnt signaling to down-regulate adipogenic transcription factors. Endocrinology 147: 141–154, 2006.[Abstract/Free Full Text]
- Singh R, Artaza JN, Taylor WE, Gonzalez-Cadavid NF, Bhasin S. Androgens stimulate myogenic differentiation and inhibit adipogenesis in C3H 10T1/2 pluripotent cells through an androgen receptor-mediated pathway. Endocrinology 144: 5081–5088, 2003.[Abstract/Free Full Text]
- Takekawa M, Saito H. A family of stress-inducible GADD45-like proteins mediate activation of the stress-responsive MTK1/MEKK4 MAPKKK. Cell 95: 521–530, 1998.[CrossRef][Web of Science][Medline]
- Trayhurn P, Duncan JS, Wood AM, Beattie JH. Metallothionein gene expression and secretion in white adipose tissue. Am J Physiol Regul Integr Comp Physiol 279: R2329–R2335, 2000.[Abstract/Free Full Text]
- Turcato E, Zamboni M, De Pergola G, Armellini F, Zivelonghi A, Bergamo-Andreis IA, Giorgino R, Bosello O. Interrelationships between weight loss, body fat distribution and sex hormones in pre- and postmenopausal obese women. J Intern Med 241: 363–372, 1997.[CrossRef][Web of Science][Medline]
- Warrington JA, Nair A, Mahadevappa M, Tsyganskaya M. Comparison of human adult and fetal expression and identification of 535 housekeeping/maintenance genes. Physiol Genomics 2: 143–147, 2000.[Abstract/Free Full Text]
- Zhan Q, Antinore MJ, Wang XW, Carrier F, Smith ML, Harris CC, Fornace AJ Jr. Association with Cdc2 and inhibition of Cdc2/Cyclin B1 kinase activity by the p53-regulated protein Gadd45. Oncogene 18: 2892–2900, 1999.[CrossRef][Web of Science][Medline]
- Zhang Y, Bosse Y, Marceau P, Biron S, Lebel S, Richard D, Vohl MC, Tchernof A. Gene expression variability in subcutaneous and omental adipose tissue of obese men. Gene Expr 14: 35–46, 2007.[CrossRef][Web of Science][Medline]
- Zuk PA, Zhu M, Ashjian P, De Ugarte DA, Huang JI, Mizuno H, Alfonso ZC, Fraser JK, Benhaim P, Hedrick MH. Human adipose tissue is a source of multipotent stem cells. Mol Biol Cell 13: 4279–4295, 2002.[Abstract/Free Full Text]
Copyright © 2008 by the American Physiological Society.