help button home button Endocrine Society Molecular Endocrinology ENDO 08 Sessions Library
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

Molecular Endocrinology, doi:10.1210/me.2003-0138
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental Data
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow Request Copyright Permission
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Ahluwalia, A.
Right arrow Articles by Waxman, D. J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Ahluwalia, A.
Right arrow Articles by Waxman, D. J.
Molecular Endocrinology 18 (3): 747-760
Copyright © 2004 by The Endocrine Society

Sexual Dimorphism of Rat Liver Gene Expression: Regulatory Role of Growth Hormone Revealed by Deoxyribonucleic Acid Microarray Analysis

Amrita Ahluwalia, Karl H. Clodfelter and David J. Waxman

Division of Cell and Molecular Biology, Department of Biology, Boston University, Boston, Massachusetts 02215

Address all correspondence and requests for reprints to: David J. Waxman, Department of Biology, Boston University, 5 Cummington Street, Boston, Massachusetts 02215. E-mail: djw{at}bu.edu.


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 RESULTS
 DISCUSSION
 MATERIALS AND METHODS
 REFERENCES
 
GH has diverse physiological actions and regulates the tissue-specific expression of numerous genes involved in growth, metabolism, and differentiation. Several of the effects of GH on somatic growth and gene expression are sex dependent and are regulated by pituitary GH secretory patterns, which are sexually differentiated. The resultant sex differences in plasma GH profiles are particularly striking in rodents and are the major determinant of sex differences in pubertal body growth rates and the expression in liver of several cytochrome P450 (CYP) enzymes that metabolize steroids, drugs, and environmental chemicals of importance to endocrinology, pharmacology, and toxicology. DNA microarray analysis was used to identify rat liver-expressed genes that show sexual dimorphism, and to ascertain the role of GH as a regulator of their sexually dimorphic expression. Adult male and female rats were untreated or were treated with GH by 7-d continuous infusion using an Alzet osmotic minipump. Poly(A) RNA was purified from individual livers and Cy3- and Cy5-labeled cDNA probes cohybridized to Pan Rat Liver and 5K Rat Oligonucleotide microarrays representing 5889 unique rat genes. Analysis of differential gene expression profiles identified 37 liver-expressed, female-predominant genes; of these, 27 (73%) were induced by continuous GH treatment of male rats. Moreover, only three of 30 genes up-regulated in male rat liver by continuous GH treatment did not display female-dominant expression. Further analysis revealed that 44 of 49 male-predominant genes (90%) were down-regulated in the livers of continuous GH-treated male rats compared with untreated male rats, whereas only five of 49 genes that were down-regulated in male rats by continuous GH treatment were not male dominant in their expression. Real-time PCR analysis applied to a sampling of 10 of the sexually dimorphic genes identified in the microarray analysis verified their sex- and GH-dependent patterns of regulation. Taken together, these studies establish that GH-regulated gene expression is the major mechanistic determinant of sexually dimorphic gene expression in the rat liver model.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 RESULTS
 DISCUSSION
 MATERIALS AND METHODS
 REFERENCES
 
THE TRADITIONAL GENE-BY-GENE approach to understanding hormone action cannot provide a full understanding of the complex regulatory processes that occur in higher eukoryotes, which may contain upwards of 35,000–40,000 individual genes, of which 10,000 or more may be expressed in any given tissue or cell type. To obtain a global view of biological processes, it is essential to obtain simultaneous readouts from a large number, if not all of its components. Methods for detecting and quantitating gene expression include Northern blots (1), differential display (2), sequencing of cDNA libraries (3, 4), and serial analysis of gene expression (5, 6). Advances in genomics, including the availability of DNA sequences covering large segments of the human and rodent genomes and the development of methods for parallel, high throughput detection and quantitation of gene expression levels, have made it possible to obtain a more global picture of cell-type and tissue-specific responses to hormones. In particular, cDNA and oligonucleotide microarrays (7, 8, 9) can be used to simultaneously monitor the expression patterns of thousands of genes in a single experiment using a systematic global strategy (10).

Using the power of microarray technology, we have undertaken a global analysis of the role of GH in the sexual dimorphism of rat liver gene expression. GH is secreted by the pituitary gland in a sex-dependent manner in both rodents (11) and humans (12). In male rats, GH is released into circulation approximately every 3.5 h, giving a repeated, pulsatile plasma hormone pattern (13) that contrasts with the more frequent, nearly continuous profile of plasma GH exposure that occurs in adult females. The resultant sex differences in plasma GH profiles are particularly striking in rodents, where they serve as a major determinant of the observed sex differences in body growth rates and the expression in liver of several cytochrome P450 (CYP) enzymes that metabolize steroids, drugs, and environmental chemicals of importance to endocrinology, pharmacology, and toxicology (14, 15). Well-studied examples of sex-specific, plasma GH pattern-regulated liver genes include the male-specific androgen 2{alpha}- and 16{alpha}-hydroxylase CYP2C11, which is strongly induced at puberty in male but not female rat liver, and the steroid sulfate 15ß-hydroxylase CYP2C12, which is exclusively expressed in adult female rat liver (16, 17). The sexual dimorphism of liver with respect to expression of these and other liver CYPs is regulated by GH at the level of transcription initiation (18, 19). Sex-dependent expression and GH regulation also characterizes other families of genes involved in steroid or foreign compound metabolism, such as the sulfotransferases (20) and class {alpha} and class µ glutathione-S-transferases (21).

DNA microarrays have been employed to determine the expression profiles of key liver-specific genes in rodents as modulated by insulin signaling (22) and small molecules such as leptin (23). Microarrays have also been used to investigate the impact of hormonal insufficiency (hypophysectomy), GH replacement therapy, and GH overexpression on liver gene expression (24, 25, 26, 27, 28). The impact of aging on gene expression and the normalizing effects of GH treatment have also been investigated (29). In one study, the effects of hypophysectomy and GH treatment on the gene expression patterns in heart, liver, and kidney were evaluated using cDNA microarrays containing 3000 different rat genes (24). cDNA microarrays have also been used to identify genes that are responsive to GH in the GH-deficient dwarf rat model (25). In another study, subtractive hybridization was used to identify 173 potential GH target genes; however, only 41 of those genes could be confirmed in subsequent cDNA array experiments (26). Of note, these studies of GH-regulated gene expression employed cDNA arrays, whose ability to discriminate between closely related DNA sequences may be limited by the cross-hybridization that invariably occurs when a single set of hybridization conditions is applied to the thousands of genes represented on the array. These ambiguities can be avoided by using microarrays constructed with oligonucleotides that are chosen to be gene specific and matched with respect to length, GC content, and melting temperature.

The present study was undertaken to investigate the role of GH in the sexual dimorphism of rat liver gene expression. Using two commercial oligonucleotide microarrays, we have identified 86 genes that show sexual dimorphism in their patterns of expression in rat liver. The large majority of these genes (72 of 86) were also found to be subject to GH regulation, demonstrating that GH is a major determinant of sexually dimorphic liver gene expression.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 RESULTS
 DISCUSSION
 MATERIALS AND METHODS
 REFERENCES
 
Oligonucleotide-based microarray analysis was applied to Poly(A)+ mRNA isolated from livers of untreated adult male, untreated adult female, and continuous GH-treated adult male rats (n = 4 rats/group). Two commercial microarrays representing a total of 5889 unique rat genes were used, one containing 1353 liver-expressed genes (Pan Rat Liver array) and the second containing 5535 genes expressed in various tissues (5K Rat array). Microarray expression data were obtained for four pairs of male and female rat liver cDNA samples, and for four pairs of male and GH-treated male samples, with each of the hybridization pairs (male vs. female and male vs. GH-treated male) analyzed both on the Pan Rat Liver array and the 5K Rat array. Given the uniform length and GC content of each oligonucleotide spotted on the array, the fluorescence hybridization intensity of each microarray spot provides an indication of the relative level of gene expression. Analysis of the frequency distributions of average spot intensities after background subtraction revealed high expression (>5000 fluorescence units) in male rat liver for 7% and approximately 4% of the genes represented on the Pan Rat Liver and 5K Rat arrays, respectively. Intermediate expression (500–5000 fluorescence units) was obtained for approximately 51% and 28% of the genes, and low but still significant expression compared with background (40–500 fluorescence units) for approximately 40% and approximately 63% of the genes on the two respective arrays. Less than 2% of the Pan Rat Liver array genes and less than 6% of the 5K Rat array genes had average spot intensities below twice the average background level.

Fold increases or decreases in gene expression compared with the untreated male group were scored as male:female and male:male + GH ratios. Genes found to be reproducibly expressed in either a sex-dependent manner (male:female >=1.5 or <=0.66) or a GH-regulated manner (male:male + GH >=1.5 or <=0.66) are listed in Tables 1–4GoGo GoGoGoGoGoGo. These regulated genes varied approximately 900-fold in their relative levels of expression after background correction, with the urinary pheromone-binding protein {alpha}-2u globulin being the most highly expressed gene in male liver (averaged normalized fluorescence intensity = 45,540), and heme-oxygenase 3 being the least expressed (average normalized fluorescence intensity = 51) (c.f. Table 1Go, right column). The log distributions of the male:female and male:male + GH ratios were compared with 0 on a per-gene basis to examine their statistical significance. Genes identified as sex-dependent and/or GH-regulated and whose distributions surpassed the significance level of P < 0.05 are shown in bold in the tables. A summary of the numbers of genes identified in each category of regulation is presented in Table 5Go.


View this table:
[in this window]
[in a new window]
 
Table 1. Male-Dominant Genes Down-Regulated in Male Rat Liver by Continuous GH Treatment

 

View this table:
[in this window]
[in a new window]
 
Table 1A. Continued

 

View this table:
[in this window]
[in a new window]
 
Table 2. Female-Dominant Genes Up-Regulated in Male Rat Liver by Continuous GH Treatment

 

View this table:
[in this window]
[in a new window]
 
Table 3A. Male-Dominant Genes Not Responsive to Continuous GH Treatment in Male Rat Liver (mean ± SE, n = 4 or 8 Individual Livers/Group)1

 

View this table:
[in this window]
[in a new window]
 
Table 3B. Female-Dominant Genes Not Responsive to Continuous GH Treatment in Male Rat Liver (mean ± SE, n = 4 or 8 Individual Livers/Group)1

 

View this table:
[in this window]
[in a new window]
 
Table 3C. Female-Dominant Genes Suppressed in Male Rat Liver by Continuous GH Treatment (mean ± SE, n = 4 or 8 Individual Livers/Group1

 

View this table:
[in this window]
[in a new window]
 
Table 4A. Non-Sex-Specific Genes Induced by Continuous GH Treatment of Male Rats (mean ± SE, n = 4 or 8 Individual Liver Comparisons)

 

View this table:
[in this window]
[in a new window]
 
Table 4B. Non-Sex-Specific Genes Suppressed by Continuous GH Treatment of Male Rats (mean ± SE, n = 4 or 8 Individual Liver Comparisons)

 

View this table:
[in this window]
[in a new window]
 
Table 5. Summary of the Gene Lists

 
Forty-nine of the approximately 6000 genes examined were found to be up-regulated in adult male as compared with adult female liver. The expression of 44 of these genes (90%) was down-regulated in male liver by continuous GH treatment (Table 1Go), whereas five of the 49 male-predominant genes were apparently not responsive to GH treatment (Tables 3AGo and 5Go). Four other genes were down-regulated by continuous GH treatment but did not display a detectable sex-difference in expression (Tables 4BGo and 5Go).

Thirty-seven genes were found to be expressed in a female-predominant fashion. Of these, 27 genes (73%) were both female predominant and GH regulated, as demonstrated by their up-regulation in continuous GH-treated male liver as compared with untreated male liver (Table 2Go). Three other genes were induced in male rat liver by GH treatment, but in contrast to the other continuous GH-inducible genes, these three genes were not expressed at a significantly higher level in female compared with male liver (Table 4AGo). Nine other genes were expressed in a female-predominant fashion but were not induced in male liver by continuous GH treatment (Table 2BGo), whereas one female-predominant gene was suppressed after GH treatment (Tables 3CGo and 5Go).

The data presented in Table 5Go, which plots sex specificity vs. GH regulation, lie heavily on the diagonal, indicating that GH is an important determinant of the sex dependence of liver gene expression. Analysis of these data showed the distribution to be highly significant (P << 0.001 by Pearson’s {chi}2 analysis), supporting the conclusion that GH is a major regulatory determinant for both the male-predominant and the female-predominant genes. Furthermore, analysis of the quantitative relationship between sex specificity and GH regulation revealed a linear correlation (r = 0.814) when male:female and the corresponding male:male + GH expression ratios were plotted for a group of 278 genes selected on the basis of their highly reproducible ratios between the different liver pairs (P < 0.1 for n = 4 or n = 8) (Fig. 1Go). A similar correlation was obtained by analysis of male:female and male:male + GH ratio distributions for approximately all 6000 genes represented on the arrays (data not shown).



View larger version (17K):
[in this window]
[in a new window]
 
Fig. 1. Scatterplot of the Male:Female vs. Male:Male + GH Expression Ratios for Subset of Microarray Genes

Shown are male:female (x-axis) and the corresponding male:male + GH expression ratios (y-axis) plotted on a logarithmic scale with the origin at (1 1 ) for 278 genes that were selected based on their highly reproducible expression ratios (P < 0.1) in either n = 4 or n = 8 arrays (genes represented on one, or both, microarray platforms, respectively). Genes mapping to the upper right quadrant (e.g. Cyp2c11, Cyp2c13, st1e2) are male specific and GH suppressed, and genes mapping to the lower left quadrant (e.g. A1bg, Cyp2c12, Sult2a1) are female specific and GH induced. The best-fit line (y = 1.017x - 0.074; r = 0.814) shows a slope close to 1 and an intercept close to 0, with a high correlation between the degree of sex specificity and GH responsiveness.

 
QPCR analysis was carried out to quantitate the relative expression levels of 17 genes selected from genes represented on the microarrays (see QPCR primer sequences published as supplemental data on The Endocrine Society’s Journals Online web site at http://mend.endojournals.org). These analyses used RNA samples prepared from a set of male, female, and continuous GH-treated male rat livers (n = 4–8 livers/group) separate from the set used in the microarray studies. Expression patterns determined by QPCR for 10 of the genes shown in Tables 1Go and 2Go are presented in Table 6Go, A and B. In all 10 cases, the overall pattern of expression revealed by QPCR was in good agreement with the results of the microarray experiments. Moreover, for seven of the 10 genes, male:female and/or male:male + GH expression ratios determined by QPCR were substantially higher (Table 1Go genes) or substantially lower (Table 2Go genes) than the corresponding ratios determined by microarray analysis (Table 6Go, A and B). In the case of one of the genes, CYP4A2, the higher expression ratio obtained by QPCR reflects the specificity of the PCR primers, in so far as the oligonucleotide used to detect CYP4A2 on the microarray cross-hybridizes with CYP4A3, whose expression in liver is not sex dependent or hormonally regulated (30). However, in the case of the other CYP genes, the microarray primers are gene specific (Fig. 2Go and data not shown). Expression ratios determined by microarray analysis may therefore understate the true extent of differential regulation of those genes that are highly regulated, as quantitated by QPCR, in agreement with another report (31).


View this table:
[in this window]
[in a new window]
 
Table 6. QPCR Analysis of Select Regulated Genes Identified by Microarray Analysis (mean ± SE, n = 4 or 8 Individual Liver Comparisons)

 


View larger version (34K):
[in this window]
[in a new window]
 
Fig. 2. Sequence Specificity of the Arrayed 50-Nucleotide-Long Oligomers Used to Detect Rat CYP2C Family Members

Shown are portions of each rat CYP2C mRNA sequence (cDNA nucleotide positions as shown on the left, and GenBank accession numbers shown on the right, as indicated). Also shown are portions of the 50-nucleotide sequences spotted on the array for CYP2C11, 2C13, and 2C23, with nucleotides that are identical shown as dots and mismatches to the CYP2C oligonucleotide sequence as indicated. A, Comparisons to the 3'-most 26 nucleotide segment of the CYP2C13 oligonucleotide spotted on the array (boxed); B, comparisons to the 3'-most 28 nucleotide segment of the CYP2C23 oligonucleotide. Sequences of each full-length CYP2C cDNA and of the indicated CYP2C oligonucleotides were aligned using DNAStar Megalign software (Madison, WI). Each of the oligonucleotides is seen to display multiple sequence differences with each of the six other CYP2C family members. Bars at top represent relative degree of sequences identity between the cDNAs for each of the positions shown. Percent sequence identities to the nontarget CYP2C genes over the entire 50-nucleotide length of each CYP2C oligonucleotide were as follows: 46–68% (CYP2C7 oligo), 38–46% (CYP2C11 oligo), 40–62% (CYP2C12 oligo), 40–76% (CYP2C13 oligo), and 36–56% (CYP2C23 oligo).

 
QPCR analysis was applied to seven other genes represented on the microarray. Five of these genes (Table 6CGo) exhibited higher expression in female and in GH-treated male liver than in untreated male liver by microarray analysis, but did not make the male:female and male:male +GH threshold ratio of 0.66 for inclusion in Table 2Go. In all five cases, QPCR verified the overall pattern of expression seen in the microarray analysis. Indeed, QPCR indicated a substantially higher degree of regulation in several cases: most notably, activated leukocyte cell adhesion molecule. QPCR also verified the unique pattern of expression seen for two other genes (Table 6DGo): hydroxysteroid-17ß dehydrogenase type 2, which was female predominant but suppressed by continuous GH treatment (Table 3CGo), and growth response protein CL-6, which was sex independent and induced by GH but did not meet the criteria for inclusion in Table 4AGo.


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 RESULTS
 DISCUSSION
 MATERIALS AND METHODS
 REFERENCES
 
Oligonucleotide-based microarrays were used to characterize the sex dependence of liver gene expression and to identify the role that GH plays in establishing and maintaining this sexual dimorphism. Traditional methods of analyzing GH-regulated genes have employed GH-deficient rodent models, e.g. hypophysectomized rats, and have identified genes such as CYPs, sulfotransferases, and glutathione-S-transferases as being dependent on GH for expression. Some liver-expressed genes are direct targets of GH action (e.g. IGF-I), whereas others, including the GH-dependent liver CYP enzymes, are regulated by GH in a complex fashion and in a manner that depends on the temporal pattern of plasma GH stimulation, which is sexually dimorphic. The major goal of this study was to obtain a global view of the sexual dimorphism of liver gene expression, and in particular, to ascertain whether GH is a major factor regulating this expression. Our findings demonstrate that the large majority of male-predominant rat liver genes identified by microarray analysis (44 of 49 genes identified) are under the control of GH, as shown by their down-regulation in male rat liver after treatment with GH for 7 d in a continuous manner (i.e. female-like plasma GH pattern) (Fig. 3Go and Table 5Go). The majority of female-predominant genes identified were also shown to be responsive to the effects of GH (27 of 37 genes), as demonstrated by their up-regulation in male rat liver after continuous GH treatment. However, nine other female-dominant genes did not apparently exhibit this GH response, suggesting that their female-predominant pattern of expression reflects other endocrine regulatory factors. In addition, three genes were expressed at similar levels in male and female liver, yet were induced by continuous GH treatment, and four genes that were not sexually dimorphic were suppressed by GH treatment. These observations indicate that the actions of GH are not restricted to sexually dimorphic liver gene expression and may include regulation of genes that are not gender specific.



View larger version (36K):
[in this window]
[in a new window]
 
Fig. 3. Venn Diagram Representations of Liver Gene Expression Profiles of Adult Male and Female Rats, and of Male Rats Treated with GH by Continuous Infusion

A, The intersection highlights 27 genes that are up-regulated in common in female and continuous GH-treated male liver as compared with untreated male rat liver (Table 2Go). Ten genes were female dominant but not induced by GH treatment (Table 3GoGo, B and C), and three genes were induced by GH treatment but were not sex specific (Table 4AGo). B, The intersection highlights 44 genes that are down-regulated in common in female and continuous GH-treated male liver as compared with untreated male rat liver (Table 1Go). Five genes (Table 3AGo) were male dominant in their expression but were not responsive to GH treatment of males, whereas four other genes were suppressed by GH treatment but were not sex specific, and one gene was suppressed by GH treatment but was female dominant (Tables 4BGo and 3CGo).

 
Using traditional methods to monitor gene expression, sexual dimorphism and/or GH regulation have been demonstrated for several of the female-predominant genes identified in the present microarray study. These include CYP2C12 (16, 17), 17ß-hydroxysteroid dehydrogenase type 2 (32), alcohol dehydrogenase (33), glutathione-S-transferase Yc2 (21), {alpha}1B-glycoprotein (34), prolactin receptor (35), and several hydroxysteroid sulfotransferases (20). Male-predominant genes listed in Table 1Go and shown by traditional methods to be subject to GH regulation include CYP2C11 (16, 17), CYP2C13 (36), CYP3A2 (37), CYP3A18 (38), CYP4A2 (30), corticosteroid 11ß-dehydrogenase (39), carbonic anhydrase 3 (40), and sulfotransferases ST1E2 and 1C1 (20, 41). One of the male-predominant genes, CYP2C22 (42), was confirmed to be male predominant, but unexpectedly, was not found to be suppressed by continuous GH treatment (Table 3AGo). Thus, there may be multiple mechanisms for regulation of male-predominant liver gene expression. Earlier studies identified two distinct classes of GH-regulated, male-specific genes (15). Members of one class, represented by CYP2C11, are obligatorily dependent on pulsatile GH for expression, whereas members of the second class, which includes CYP2A2, CYP3A2, and CYP4A2, do not exhibit this GH dependence, as indicated by their full expression in hypophysectomized male rat liver and by their up-regulation to near male levels after hypophysectomy in female rat liver (16). However, both classes of male-specific genes are down-regulated by continuous GH treatment, a response that is distinct from that exhibited by CYP2C22 in the present study. Additional studies are required to further characterize the male-predominant, GH-regulated genes presently identified to determine which of these classes they belong to, and to elucidate the mechanisms that underlie their distinct responses to GH ablation by hypophysectomy.

cDNA-based microarrays were used in several earlier studies to identify target genes for GH action in the paradigm of liver gene expression (24, 25, 26, 27). Although cDNA and oligonucleotide-based microarrays exhibit similar sensitivities for detecting changes in gene expression (31, 43), cDNA microarrays can be limited by their ability to discriminate between closely related DNA sequences because of the cross-hybridization that invariably occurs when a single arbitrary set of hybridization conditions is applied uniformly across the thousands of genes represented on the cDNA array (c.f. Refs. 44 and 45). The present study was carried out using microarrays constructed from 50-nucleotide-long oligonucleotides chosen to maximize gene specificity and matched with respect to length, GC content, melting temperature and designed to be free from secondary structure and self-annealing tendencies. These oligonucleotides were generally effective in distinguishing between closely related cDNAs within complex gene families and subfamilies, as exemplified by CYPs belonging to subfamily CYP2C. This gene subfamily includes seven closely related rat genes (75–86% nucleotide identity), all of which are expressed in liver, and six of which (CYPs 2C7, 2C11, 2C12, 2C13, 2C22, and 2C23) were represented on the arrays used in the present study. The percent identity between any given 50-nucleotide-long CYP2C oligomer on the array and nontarget CYP2C transcripts was generally between 40% and 70%, and in all cases was devoid of complementary sequences greater than 15 contiguous bases (Fig. 2Go and data not shown). Cross-hybridization is unlikely under these conditions (43). As noted above, several rat CYP2C genes are GH regulated and sexually dimorphic in their expression, with CYPs 2C11, 2C13, and 2C22 being male dominant and CYPs 2C7 and 2C12 female dominant. These patterns of expression were verified in the present study, supporting the specificity of the oligonucleotide probes spotted on the arrays.

The two commercial microarrays used in this study encompassed a total of 5889 rat genes, enabling us to identify both known genes and novel genes that display sexual dimorphism and/or regulation by GH, including six sequences not represented in GenBank. Unlike a previous study of GH-regulated liver genes (24), the present study specifically sought to identify liver-expressed genes that are both sex dependent and GH regulated. Many of the sex-dependent, GH-regulated genes that we identified are involved either in steroid synthesis or in the metabolism of steroids and foreign compounds. Others encode serum proteins, structural proteins, receptors, membrane trafficking proteins, and transporters. Several of the genes identified code for signaling molecules such as nuclear factors and phosphatases, some of which may potentially help mediate the effects of GH on sex-dependent liver gene expression. GH-regulated signaling molecules that are more highly expressed in female than male liver include the H-ras-like tumor suppressor H-rev 107, fibroblast growth factor 5 and several hormone, cytokine, and G protein-coupled receptors (Table 2Go). Potential signaling molecules found to be GH-regulated and more highly expressed in males include phosphodiesterase I/nucleotide pyrophosphatase 2, the protein kinase phosphorylase kinase {gamma}-2, muscarinic acetylcholine receptor M2, and the cell growth regulator and gap junction component connexin 30 (Table 1Go). Further investigation will be required to ascertain whether these changes in mRNA expression are accompanied by corresponding changes in protein levels, and to elucidate the precise role that these proteins play in liver physiology and in the sexual dimorphism of liver function.

Many of the genes presently identified displayed moderate regulation by GH (i.e. 1.5- to 3-fold differences in expression in response to changes in GH status), as compared with the more than 100-fold differences in gene expression seen with the classic sexually dimorphic, GH pattern-regulated genes CYP2C11 and CYP2C12. These striking quantitative differences in the extent of sexual dimorphism and GH regulation were verified by QPCR (Table 6Go). This finding suggests that intrinsic signals distinct from GH, both inductive and repressive, contribute to the observed sexually dimorphic pattern of liver gene expression. QPCR also verified the expression patterns for several genes identified in the microarray studies that were not previously known to be sexually dimorphic and GH-regulated (transcobalamin I precursor, phosphodiesterase I/nucleotide pyrophosphatase 2, acyl-coenzyme A (CoA) synthase 5, and H-ras-like tumor suppressor H-rev 107). Several other genes that were found by microarray analysis to be more highly expressed in female than in male liver, but did not meet our threshold ratio criteria for inclusion in Table 2Go, were also analyzed by QPCR and shown to be female predominant and induced by continuous GH treatment (Table 6CGo). The absence of these five genes from Table 2Go serves as an example of false negatives, which occur frequently in microarray analysis. False positives are also an unavoidable feature of microarray analysis. However, this possibility may be minimized by the requirement that the genes listed in Tables 1–4GoGoGoGoGoGoGo show a consistent pattern of regulation in at least three of four independent pairs of male:female livers and a corresponding number of male:GH-treated male livers, and by the fact that a majority of the genes listed in Tables 1–4GoGoGoGoGoGoGo gave expression ratios that were tightly clustered in the four pairs of livers analyzed, and thus highly significant (table entries shown in bold). Nevertheless, it will be important to validate the patterns of gene expression profiling obtained in the present microarray studies through the use of alternative, more traditional methods, such as the QPCR analyses presented in Table 6Go.

Whereas the microarray analysis reported here is by no means comprehensive, it has led to the identification of several distinct groups (classes) of liver-expressed genes, each of which is unique with respect to its responsiveness to sex-dependent hormonal factors, of which GH is seen to be the major, but apparently not the sole determinant (Table 5Go). Conceivably, genes within each of the classes represented by Tables 1Go and 2Go are likely to share common regulatory mechanisms that enable them to respond to continuous GH treatment by gene up-regulation (Table 2Go) or down-regulation (Table 1Go). As GH-regulated, sexual dimorphic liver gene expression is primarily, if not exclusively, controlled at the level of transcription initiation (18, 19), future studies designed to elucidate common promoter region sequences and upstream regulatory elements may provide further insight into the mechanisms that underlie the sexual dimorphism of hepatic gene expression.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 RESULTS
 DISCUSSION
 MATERIALS AND METHODS
 REFERENCES
 
Animals
Adult male and female Fischer 344 rats (8–10 wk of age) (n = 4/group) (Taconic, Inc., Germantown, NY) were untreated or were treated with rat GH (rGH-B-14-SIAFP, National Institute of Diabetes and Digestive and Kidney Diseases), given as a continuous infusion using an Alzet osmotic minipump delivering 2 µg GH/100 g body weight·h for 7 d. Intact, untreated male and female rats were killed and livers were collected, frozen in liquid N2 and stored at -80 C.

RNA Isolation
Total liver RNA was isolated from frozen liver samples using Trizol reagent (Invitrogen Life Technologies, Carlsbad, CA). Poly(A) RNA was then isolated using a commercial kit (Ambion, Austin, TX). RNA samples having an A260/A280 ratio of at least 1.8 were deemed suitable for use in subsequent analyses. Total liver RNA and Poly(A) RNA samples were treated with deoxyribonuclease I (Ambion, Austin, TX) as described previously (46) and their integrity was checked on Northern blots probed for CYP2C6 mRNA (47), which is expressed at high levels in both male and female rat liver.

Oligonucleotide Microarray Experiments
The expression of sex-dependent and GH-regulated liver mRNAs was investigated using two commercial DNA microarrays, Pan Rat Liver Oligonucleotide array and 5K Rat Oligonucleotide array (MWG Biotech, Inc., High Point, NC). These arrays are respectively comprised of oligonucleotide sequences 50 nucleotides in length selected to represent unique (gene specific) sequences, generally derived from the coding sequence or 3'-untranslated region of each mRNA. The Pan Rat Liver and 5K Rat arrays used respectively represent 1353 unique liver-expressed genes (1408 features, including replicates; Pan Rat Liver array) and 5535 unique rat genes (5760 features, including replicates; 5K Rat array). Together, these microarrays represent 5889 unique rat genes, of which 999 genes are represented on both arrays. Poly(A) RNA purified from individual livers was labeled with Cy3-deoxyuridine triphosphate (dUTP) or Cy5-dUTP (PerkinElmer Life Sciences Inc., Boston, MA) in a reverse transcription reaction using 1 µg Poly(A) RNA and Superscript II (ribonuclease H mutant) enzyme (Invitrogen Life Technologies, Carlsbad, CA) at 39 C for 2 h. Unincorporated deoxynucleotide triphosphates were removed using the Qiaquick PCR-purification kit (QIAGEN, Valencia, CA). Cy5-labeled male liver cDNA was mixed with Cy3-labeled female or Cy3-labeled GH-treated male liver cDNA followed by cohybridization of the mixed cDNA population to the microarray for 18 h at 42 C. Arrays were washed according to the manufacturer’s protocol. RNA from each of four individual sets of male, female and GH-treated male liver was converted to cDNA and labeled in duplicate. One series of eight hybridizations (four hybridizations with male + female liver cDNA, and four hybridizations with male + GH-treated male liver cDNA; n = 4 livers/treatment group) was carried out with the Pan Rat Liver arrays, and a second series containing the same eight hybridizations (n = 4 livers/group) was carried out with the 5K Rat arrays. Dye swapping experiments using a single pair of male and female RNA samples, respectively labeled with Cy3- and Cy5-dUTP, were performed to verify the independence of the microarray results from the nature of the cyanine dye label.

Data Acquisition and Analysis
Fluorescent cDNA bound to each microarray was detected with a GenePix 4000B array scanner (Axon Instruments, Foster City, CA) using a fixed laser power intensity and the minimum photomultiplier tube setting necessary to avoid bleaching while allowing for detection of features with faint fluorescent intensities. The manufacturer’s GenePix 3.0 software package was used to locate individual spots, quantitate the Cy3- and Cy5-fluorescence intensity at each spot, and determine background signal intensities. Data from spots that were determined to result from hybridization anomalies or microarray errors (e.g. dust particles, fibers, or spotting anomalies including bleeding around the actual spotting area) were excluded from the analysis. The data were imported into GeneSpring 6.0 (Silicon Genetics, Redwood City, CA). Fluorescence intensity values were calculated by subtracting the local median background value from the median foreground value for each spot. The use of median values instead of mean values reduced the effect of outliers generated by imperfections in GenePix’s automatic spot detection protocol. Average background values were determined to be 40 ± 8 (Cy5 channel) and 51 ± 13 (Cy3 channel), mean ± SD. Lowess normalization, using 20% of the data for smoothing, was applied to the Cy3 channel for each microarray spot on each array using GeneSpring software. Lowess-normalized Cy3 values below 10 were set to a baseline value of 10. The ratio of the normalized channels (Cy5/Cy3) as reported by GeneSpring was used for expression analysis. All genes included in our analysis and listed in Tables 1–4GoGoGoGoGoGoGo exhibited expression differences of 1.5-fold or higher (i.e. at least 50% greater expression in the Cy5 sample than in the Cy3 sample) or 0.66-fold or lower that were reproducible in microarray slides representing at least three of four male:female or male:male+ GH liver pairs, or at least two of three liver pairs for those genes where only three of the four microarrays gave data of suitable technical quality. The threshold values of >=1.5 and <=0.66 were chosen so as to not exclude genes whose expression is known to be sex dependent and/or GH regulated by other, independent studies, but whose microarray expression ratios did not satisfy a 2.0-threshold value. Genes with reproducible expression ratios between 0.8 and 1.2 were not considered to be sex specific (male:female comparisons) or GH regulated (male:male + GH comparisons). For those genes that are represented on both microarrays, the mean ratio and SE values were calculated based on n = 8 replicates. Expression ratios for all other genes were calculated based on n = 4 unless one of the replicates returned no data for that spot. Average fluorescence intensity values reported in Tables 1–4GoGoGoGoGoGoGo are normalized values and were corrected for background fluorescent intensity. Expression ratios for the approximately 6000 genes examined on the 16 microarray slides used in this study are available as published as supplemental data on The Endocrine Society’s Journals Online web site at http://mend.endojournals.org. The data are also available for query or download from the Gene Expression Omnibus (GEO) web site at NCBI (http://www.ncbi.nlm.nih.gov/geo).

Statistical Analysis
A one-sample t test (GeneSpring 6.0 software) was applied to the distribution of natural logs of the ratios for each gene. The t test implemented by the GeneSpring software package calculates the P value for the distribution of natural logs as compared with 0. A filter (P < 0.05) was applied to the P values to determine the statistical significance of each gene’s differential expression for each DNA array type (Pan Rat Liver and 5K Rat array). Expression ratios that met the P < 0.05 criterion are shown in bold in Tables 1–4GoGoGoGoGoGoGo. Multiple testing correction methods, such as Bonferroni or Holm step-down, were not applied to the P values because these options, available in the GeneSpring 6.0 software package, depend heavily on the independence of each gene’s expression and thus filter out many bona fide regulated genes to avoid all type I errors; they are thus too restrictive in their effort to avoid false positives (48).

A Pearson’s {chi}2 test was applied using MATLAB 6.1.0.450 (The MathWorks, Inc., Natick, MA) on a table of the gene counts for each observed relation between sex dominance and GH regulation (Table 5Go) to validate the correlation between the two experimental paradigms (sex specificity and GH regulation). The numbers for these relationships were set into a 3 x 3 table and the test was applied to determine whether the observed distribution differs significantly from the expected normal distribution. Any element in the table less than five was increased by pseudo-counts to five to fit the assumptions of the {chi}2 test (i.e. no observation should be less than 5). The value of the non-GH-responsive, non-sex-specific genes was set at an underestimated value of 2475, to take into account potential experimental error. A {chi}2 value of 2993, corresponding to P << 0.001, was calculated. The same high level of significance was found when a second calculation for the {chi}2 test was made underestimating the value of the non-GH-responsive, non-sex-specific genes as 20, to reduce the bias on the expected distribution due to the large difference between 2475 and the other table entries. An equally high level of significance was found when a {chi}2 test for a 2 x 2 table of only up-regulated and down-regulated genes for each condition was calculated.

Real-Time PCR Studies
The expression profiles of select genes represented on the microarrays were examined by a reverse transcription real-time QPCR method using an ABI Prism 7700HT instrument (Applied Biosystems, Foster City, CA) and SYBR green detection. cDNA was prepared from individual livers, and each cDNA was assayed in triplicate PCRs. Briefly, 1 µg of total liver RNA was treated with deoxyribonuclease I and reverse transcribed with oligo deoxythymidine20-VN primers (V = any nucleotide but thymidine) using the GeneAmp RNA-PCR kit (Applied Biosystems) in a reaction volume of 20 µl as recommended by the manufacturer and described earlier (46). cDNA samples were used at a dilution of 1:100 in QPCR assays using the SYBR Green mix kit (Applied Biosystems, Foster City, CA). PCRs (4–5 µl total volume per well) were carried out in triplicate in individual wells of a 384-well plate. With each cycle of PCR, more SYBR Green dye molecules bind to the newly synthesized double-stranded DNA amplicon, resulting in an increase in fluorescence in real-time. Gene-specific primers were designed using Primer Express software (Applied Biosystems). Primer sequences used for the 17 genes shown in Table 6Go are published as supplemental data on The Endocrine Society’s Journals Online web site at http://mend.endojournals.org. Primer specificity and the extent of possible cross-hybridization with other genes in GenBank was evaluated by basic local alignment and search tool analysis (BLAST). Primer concentrations of 300–600 nM were generally found to be optimal. Assays were performed using liver cDNA prepared from n = 4–8 individual male, female and continuous GH-treated male livers, all separate from the livers used in the microarray analysis. QPCR amplification of 18S rRNA was carried out in parallel to normalize the input cDNA for each PCR. Products of mock reverse transcription reactions for each of the liver RNA samples were included as negative controls to verify the absence of amplification signal attributable to contamination by genomic DNA.


    ACKNOWLEDGMENTS
 
The authors gratefully acknowledge Binyamin Berkovits and Christopher Wiwi for assistance with primer design and QPCR analysis of GH target genes.


    FOOTNOTES
 
This work was supported in part by NIH Grant DK33765 (to D.J.W.). K.H.C. received support from the Superfund Basic Research Center at Boston University, NIH Grant 5-P42-ES07381.

Abbreviations: CoA, Coenzyme A; CYP, cytochrome P450; dUTP, deoxyuridine triphosphate; QPCR, quantitative PCR.

Received for publication April 15, 2003. Accepted for publication December 10, 2003.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 RESULTS
 DISCUSSION
 MATERIALS AND METHODS
 REFERENCES
 

  1. Alwine JC, Kemp DJ, Stark GR 1977 Method for detection of specific RNAs in agarose gels by transfer to diazobenzyloxymethyl-paper and hybridization with DNA probes. Proc Natl Acad Sci USA 74:5350–5354[Abstract/Free Full Text]
  2. Liang P, Pardee AB 1992 Differential display of eukaryotic messenger RNA by means of the polymerase chain reaction. Science 257:967–971[Abstract/Free Full Text]
  3. Adams MD, Kelley JM, Gocayne JD, Dubnick M, Polymeropoulos MH, Xiao H, Merril CR, Wu A, Olde B, Moreno RF, et al 1991 Complementary DNA sequencing: expressed sequence tags and human genome project. Science 252:1651–1656[Abstract/Free Full Text]
  4. Okubo K, Hori N, Matoba R, Niiyama T, Fukushima A, Kojima Y, Matsubara K 1992 Large scale cDNA sequencing for analysis of quantitative and qualitative aspects of gene expression. Nat Genet 2:173–179[CrossRef][Medline]
  5. Velculescu VE, Zhang L, Vogelstein B, Kinzler KW 1995 Serial analysis of gene expression. Science 270:484–487[Abstract/Free Full Text]
  6. Hibi K, Liu Q, Beaudry GA, Madden SL, Westra WH, Wehage SL, Yang SC, Heitmiller RF, Bertelsen AH, Sidransky D, Jen J 1998 Serial analysis of gene expression in non-small cell lung cancer. Cancer Res 58:5690–5694[Abstract/Free Full Text]
  7. Schena M, Shalon D, Heller R, Chai A, Brown PO, Davis RW 1996 Parallel human genome analysis: microarray-based expression monitoring of 1000 genes. Proc Natl Acad Sci USA 93:10614–10619[Abstract/Free Full Text]
  8. Lockhart DJ, Dong H, Byrne MC, Follettie MT, Gallo MV, Chee MS, Mittmann M, Wang C, Kobayashi M, Horton H, Brown EL 1996 Expression monitoring by hybridization to high-density oligonucleotide arrays. Nat Biotechnol 14:1675–1680[CrossRef][Medline]
  9. Lipshutz RJ, Fodor SP, Gingeras TR, Lockhart DJ 1999 High density synthetic oligonucleotide arrays. Nat Genet 21:20–24[CrossRef][Medline]
  10. Brown PO, Botstein D 1999 Exploring the new world of the genome with DNA microarrays. Nat Genet 21:33–37[CrossRef][Medline]
  11. Jansson JO, Eden S, Isaksson O 1985 Sexual dimorphism in the control of growth hormone secretion. Endocr Rev 6:128–150[Abstract/Free Full Text]
  12. Veldhuis JD 1996 Gender differences in secretory activity of the human somatotropic (growth hormone) axis. Eur J Endocrinol 134:287–295[Abstract/Free Full Text]
  13. Tannenbaum GS, Martin JB 1976 Evidence for an endogenous ultradian rhythm governing growth hormone secretion in the rat. Endocrinology 98:562–570[Abstract/Free Full Text]
  14. Shapiro BH, Agrawal AK, Pampori NA 1995 Gender differences in drug metabolism regulated by growth hormone. Int J Biochem Cell Biol 27:9–20[CrossRef][Medline]
  15. Waxman DJ, Chang TKH 1995 Cytochrome P-450: structure, mechanism, and biochemistry. In: Ortiz de Montellano PR, ed. Cytochrome P-450: structure, mechanism, and biochemistry. New York: Plenum Press; 391–417
  16. Waxman DJ 1992 Regulation of liver-specific steroid metabolizing cytochromes P450: cholesterol 7{alpha}-hydroxylase, bile acid 6ß-hydroxylase, and growth hormone-responsive steroid hormone hydroxylases. J Ster Biochem Mol Biol 43:1055–1072[CrossRef]
  17. Mode A 1993 Sexually differentiated expression of genes encoding the P4502C cytochromes in rat liver—a model system for studying the action of growth hormone. J Reprod Fertil Suppl 46:77–86[Medline]
  18. Sundseth SS, Alberta JA, Waxman DJ 1992 Sex-specific, growth hormone-regulated transcription of the cytochrome P450 2C11 and 2C12 genes. J Biol Chem 267:3907–3914[Abstract/Free Full Text]
  19. Legraverend C, Mode A, Westin S, Strom A, Eguchi H, Zaphiropoulos PG, and Gustafsson J-A 1992 Transcriptional regulation of rat P-450 2C gene subfamily members by the sexually dimorphic pattern of growth hormone secretion. Mol Endocrinol 6:259–266[Abstract/Free Full Text]
  20. Klaassen CD, Liu L, Dunn RT 1998 Regulation of sulfotransferase mRNA expression in male and female rats of various ages. Chem Biol Interact 109:299–313[CrossRef][Medline]
  21. Srivastava PK, Waxman DJ 1993 Sex-dependent expression and growth hormone regulation of class {alpha} and class µ glutathione S-transferase mRNAs in adult rat liver. Biochem J 294:159–165
  22. Tobe K, Suzuki R, Aoyama M, Yamauchi T, Kamon J, Kubota N, Terauchi Y, Matsui J, Akanuma Y, Kimura S, Tanaka J, Abe M, Ohsumi J, Nagai R, Kadowaki T 2001 Increased expression of the sterol regulatory element-binding protein-1 gene in insulin receptor substrate-2(-/-) mouse liver. J Biol Chem 276:38337–38340[Abstract/Free Full Text]
  23. Ferrante Jr AW, Thearle M, Liao T, Leibel RL 2001 Effects of leptin deficiency and short-term repletion on hepatic gene expression in genetically obese mice. Diabetes 50:2268–2278[Abstract/Free Full Text]
  24. Flores-Morales A, Stahlberg N, Tollet-Egnell P, Lundeberg J, Malek RL, Quackenbush J, Lee NH, Norstedt G 2001 Microarray analysis of the in vivo effects of hypophysectomy and growth hormone treatment on gene expression in the rat. Endocrinology 142:3163–3176[Abstract/Free Full Text]
  25. Thompson BJ, Shang CA, Waters MJ 2000 Identification of genes induced by growth hormone in rat liver using cDNA arrays. Endocrinology 141:4321–4324[Abstract/Free Full Text]
  26. Gardmo C, Swerdlow H, Mode A 2002 Growth hormone regulation of rat liver gene expression assessed by SSH and microarray. Mol Cell Endocrinol 190:125–133[CrossRef][Medline]
  27. Lu C, Kasik J, Stephan DA, Yang S, Sperling MA, Menon RK 2001 Grtp1, a novel gene regulated by growth hormone. Endocrinology 142:4568–4571[Abstract/Free Full Text]
  28. Olsson B, Bohlooly YM, Brusehed O, Isaksson OG, Ahren B, Olofsson SO, Oscarsson J, Tornell J 2003 Bovine growth hormone-transgenic mice have major alterations in hepatic expression of metabolic genes. Am J Physiol Endocrinol Metab 285:E504–E511
  29. Tollet-Egnell P, Flores-Morales A, Stahlberg N, Malek RL, Lee N, Norstedt G 2001 Gene expression profile of the aging process in rat liver: normalizing effects of growth hormone replacement. Mol Endocrinol 15:308–318[Abstract/Free Full Text]
  30. Sundseth SS, Waxman DJ 1992 Sex-dependent expression and clofibrate inducibility of cytochrome P450 4A fatty acid {omega}-hydroxylases. Male specificity of liver and kidney CYP4A2 mRNA and tissue-specific regulation by growth hormone and testosterone. J Biol Chem 267:3915–3921[Abstract/Free Full Text]
  31. Yuen T, Wurmbach E, Pfeffer RL, Ebersole BJ, Sealfon SC 2002 Accuracy and calibration of commercial oligonucleotide and custom cDNA microarrays. Nucleic Acids Res 30:e48
  32. Akinola LA, Poutanen M, Vihko R 1996 Cloning of rat 17ß-hydroxysteroid dehydrogenase type 2 and characterization of tissue distribution and catalytic activity of rat type 1 and type 2 enzymes. Endocrinology 137:1572–1579[Abstract]
  33. Crabb DW, Roepke J 1987 Loss of growth hormone-dependent characteristics of rat hepatocytes in culture. In Vitro Cell Dev Biol 23:303–307[Medline]
  34. Gardmo C, Persson B, Mode A 2001 Cloning of a novel growth hormone-regulated rat complementary deoxyribonucleic acid with homology to the human {alpha}1B-glycoprotein, characterizing a new protein family. Endocrinology 142:2695–2701[Abstract/Free Full Text]
  35. Robertson JA, Haldosen LA, Wood TJ, Steed MK, Gustafsson JA 1990 Growth hormone pretranslationally regulates the sexually dimorphic expression of the prolactin receptor gene in rat liver. Mol Endocrinol 4:1235–1239[Abstract/Free Full Text]
  36. McClellan-Green PD, Negishi M, Goldstein JA 1989 Characterization of a cDNA for rat P-450g, a highly polymorphic, male-specific cytochrome in the P-450IIC subfamily. Biochemistry 28:5832–839[CrossRef][Medline]
  37. Waxman DJ, Dannan GA, Guengerich FP 1985 Regulation of rat hepatic cytochrome P-450: age-dependent expression, hormonal imprinting, and xenobiotic inducibility of sex-specific isoenzymes. Biochemistry 24:4409–4417[CrossRef][Medline]
  38. Robertson GR, Farrell GC, Liddle C 1998 Sexually dimorphic expression of rat CYP3A9 and CYP3A18 genes is regulated by growth hormone. Biochem Biophys Res Commun 242:57–60[CrossRef][Medline]
  39. Agarwal AK, Monder C, Eckstein B, White PC 1989 Cloning and expression of rat cDNA encoding corticosteroid 11 ß-dehydrogenase. J Biol Chem 264:18939–8943[Abstract/Free Full Text]
  40. Jeffery S, Carter ND, Clark RG, Robinson IC 1990 The episodic secretory pattern of growth hormone regulates liver carbonic anhydrase III. Studies in normal and mutant growth-hormone-deficient dwarf rats. Biochem J 266:69–74[Medline]
  41. Nagata K, Ozawa S, Miyata M, Shimada M, Gong DW, Yamazoe Y, Kato R 1993 Isolation and expression of a cDNA encoding a male-specific rat sulfotransferase that catalyzes activation of N-hydroxy-2-acetylaminofluorene. J Biol Chem 268:24720–4725[Abstract/Free Full Text]
  42. Nagata K, Shimada M, Yamazoe Y, Kato R 1991 Regulations of male-dominant P-450Md mRNA in rat liver by hormonal factors and xenobiotics. Jpn J Pharmacol 57:123–126[Medline]
  43. Kane MD, Jatkoe TA, Stumpf CR, Lu J, Thomas JD, Madore SJ 2000 Assessment of the sensitivity and specificity of oligonucleotide (50mer) microarrays. Nucleic Acids Res 28:4552–4557[Abstract/Free Full Text]
  44. Li J, Pankratz M, Johnson JA 2002 Differential gene expression patterns revealed by oligonucleotide versus long cDNA arrays. Toxicol Sci 69:383–390[Abstract/Free Full Text]
  45. Rae JM, Johnson MD, Lippman ME, Flockhart DA 2001 Rifampin is a selective, pleiotropic inducer of drug metabolism genes in human hepatocytes: studies with cDNA and oligonucleotide expression arrays. J Pharmacol Exp Ther 299:849–857[Abstract/Free Full Text]
  46. Ahluwalia A, Hurteau JA, Bigsby RM, Nephew KP 2001 DNA methylation in ovarian cancer. II. Expression of DNA methyltransferases in ovarian cancer cell lines and normal ovarian epithelial cells. Gynecol Oncol 82:299–304[CrossRef][Medline]
  47. Waxman DJ 1991 Rat hepatic P450IIA and P450IIC subfamily expression using catalytic, immunochemical, and molecular probes. Methods Enzymol 206:249–267[Medline]
  48. Draghici S 2002 Statistical intelligence: effective analysis of high-density microarray data. Drug Discov Today 7:S55–S63



This article has been cited by other articles:


Home page
Mol. Endocrinol.Home page
R. D. Meyer, E. V. Laz, T. Su, and D. J. Waxman
Male-Specific Hepatic Bcl6: Growth Hormone-Induced Block of Transcription Elongation in Females and Binding to Target Genes Inversely Coordinated with STAT5
Mol. Endocrinol., November 1, 2009; 23(11): 1914 - 1926.
[Abstract] [Full Text] [PDF]


Home page
Mol. Endocrinol.Home page
E. V. Laz, A. Sugathan, and D. J. Waxman
Dynamic in Vivo Binding of STAT5 to Growth Hormone-Regulated Genes in Intact Rat Liver. Sex-Specific Binding at Low- But Not High-Affinity STAT5 Sites
Mol. Endocrinol., August 1, 2009; 23(8): 1242 - 1254.
[Abstract] [Full Text] [PDF]


Home page
J. Biol. Chem.Home page
I. M. Bur, A. M. Cohen-Solal, D. Carmignac, P.-Y. Abecassis, N. Chauvet, A. O. Martin, G. T. J. van der Horst, I. C. A. F. Robinson, P. Maurel, P. Mollard, et al.
The Circadian Clock Components CRY1 and CRY2 Are Necessary to Sustain Sex Dimorphism in Mouse Liver Metabolism
J. Biol. Chem., April 3, 2009; 284(14): 9066 - 9073.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Regul. Integr. Comp. Physiol.Home page
J. S. Gilbert and M. J. Nijland
Sex differences in the developmental origins of hypertension and cardiorenal disease
Am J Physiol Regulatory Integrative Comp Physiol, December 1, 2008; 295(6): R1941 - R1952.
[Abstract] [Full Text] [PDF]


Home page
Mol. Endocrinol.Home page
V. Wauthier and D. J. Waxman
Sex-Specific Early Growth Hormone Response Genes in Rat Liver
Mol. Endocrinol., August 1, 2008; 22(8): 1962 - 1974.
[Abstract] [Full Text] [PDF]


Home page
Toxicol SciHome page
J. B. Silkworth, E. A. Carlson, C. McCulloch, K. Illouz, S. Goodwin, and T. R. Sutter
Toxicogenomic Analysis of Gender, Chemical, and Dose Effects in Livers of TCDD- or Aroclor 1254-Exposed Rats Using a Multifactor Linear Model
Toxicol. Sci., April 1, 2008; 102(2): 291 - 309.
[Abstract] [Full Text] [PDF]


Home page
Cancer Res.Home page
A. B. Rogers, E. J. Theve, Y. Feng, R. C. Fry, K. Taghizadeh, K. M. Clapp, C. Boussahmain, K. S. Cormier, and J. G. Fox
Hepatocellular Carcinoma Associated with Liver-Gender Disruption in Male Mice
Cancer Res., December 15, 2007; 67(24): 11536 - 11546.
[Abstract] [Full Text] [PDF]


Home page
Drug Metab. Dispos.Home page
R. N. Dhir, C. Thangavel, and B. H. Shapiro
Attenuated Expression of Episodic Growth Hormone-Induced CYP2C11 in Female Rats Associated with Suboptimal Activation of the Jak2/Stat5B and Other Modulating Signaling Pathways
Drug Metab. Dispos., November 1, 2007; 35(11): 2102 - 2110.
[Abstract] [Full Text] [PDF]


Home page
Physiol. GenomicsHome page
K. H. Clodfelter, G. D. Miles, V. Wauthier, M. G. Holloway, X. Zhang, P. Hodor, W. J. Ray, and D. J. Waxman
Role of STAT5a in regulation of sex-specific gene expression in female but not male mouse liver revealed by microarray analysis
Physiol Genomics, September 11, 2007; 31(1): 63 - 74.
[Abstract] [Full Text] [PDF]


Home page
J. Nutr.Home page
I. Ueki and M. H. Stipanuk
Enzymes of the Taurine Biosynthetic Pathway Are Expressed in Rat Mammary Gland
J. Nutr., August 1, 2007; 137(8): 1887 - 1894.
[Abstract] [Full Text] [PDF]


Home page
Physiol. GenomicsHome page
E. M. Santos, V. L. Workman, G. C. Paull, A. L. Filby, K. J. W. Van Look, P. Kille, and C. R. Tyler
Molecular basis of sex and reproductive status in breeding zebrafish
Physiol Genomics, July 18, 2007; 30(2): 111 - 122.
[Abstract] [Full Text] [PDF]


Home page
EndocrinologyHome page
E. V. Laz, M. G. Holloway, C.-S. Chen, and D. J. Waxman
Characterization of Three Growth Hormone-Responsive Transcription Factors Preferentially Expressed in Adult Female Liver
Endocrinology, July 1, 2007; 148(7): 3327 - 3337.
[Abstract] [Full Text] [PDF]


Home page
EndocrinologyHome page
H. Jiang, Y. Wang, M. Wu, Z. Gu, S. J. Frank, and R. Torres-Diaz
Growth Hormone Stimulates Hepatic Expression of Bovine Growth Hormone Receptor Messenger Ribonucleic Acid through Signal Transducer and Activator of Transcription 5 Activation of a Major Growth Hormone Receptor Gene Promoter
Endocrinology, July 1, 2007; 148(7): 3307 - 3315.
[Abstract] [Full Text] [PDF]


Home page
Mol. Endocrinol.Home page
D. J. Waxman and C. O'Connor
Growth Hormone Regulation of Sex-Dependent Liver Gene Expression
Mol. Endocrinol., November 1, 2006; 20(11): 2613 - 2629.
[Abstract] [Full Text] [PDF]


Home page
J Mol EndocrinolHome page
M. L Rise, S. E Douglas, D. Sakhrani, J. Williams, K V. Ewart, M. Rise, W. S Davidson, B. F Koop, and R. H Devlin
Multiple microarray platforms utilized for hepatic gene expression profiling of GH transgenic coho salmon with and without ration restriction.
J. Mol. Endocrinol., October 1, 2006; 37(2): 259 - 282.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Endocrinol. Metab.Home page
K. E. Govoni, S. K. Lee, R. B. Chadwick, H. Yu, Y. Kasukawa, D. J. Baylink, and S. Mohan
Whole genome microarray analysis of growth hormone-induced gene expression in bone: T-box3, a novel transcription factor, regulates osteoblast proliferation
Am J Physiol Endocrinol Metab, July 1, 2006; 291(1): E128 - E136.
[Abstract] [Full Text] [PDF]


Home page
Mol. Endocrinol.Home page
K. H. Clodfelter, M. G. Holloway, P. Hodor, S.-H. Park, W. J. Ray, and D. J. Waxman
Sex-Dependent Liver Gene Expression Is Extensive and Largely Dependent upon Signal Transducer and Activator of Transcription 5b (STAT5b): STAT5b-Dependent Activation of Male Genes and Repression of Female Genes Revealed by Microarray Analysis
Mol. Endocrinol., June 1, 2006; 20(6): 1333 - 1351.
[Abstract] [Full Text] [PDF]


Home page
Mol. Endocrinol.Home page
M. G. Holloway, E. V. Laz, and D. J. Waxman
Codependence of Growth Hormone-Responsive, Sexually Dimorphic Hepatic Gene Expression on Signal Transducer and Activator of Transcription 5b and Hepatic Nuclear Factor 4{alpha}
Mol. Endocrinol., March 1, 2006; 20(3): 647 - 660.
[Abstract] [Full Text] [PDF]


Home page
J. Pharmacol. Exp. Ther.Home page
C. Cheung, A.-M. Yu, C.-S. Chen, K. W. Krausz, L. G. Byrd, L. Feigenbaum, R. J. Edwards, D. J. Waxman, and F. J. Gonzalez
Growth Hormone Determines Sexual Dimorphism of Hepatic Cytochrome P450 3A4 Expression in Transgenic Mice
J. Pharmacol. Exp. Ther., March 1, 2006; 316(3): 1328 - 1334.
[Abstract] [Full Text] [PDF]


Home page
J. Biol. Chem.Home page
T. Landsman and D. J. Waxman
Role of the Cytokine-induced SH2 Domain-containing Protein CIS in Growth Hormone Receptor Internalization
J. Biol. Chem., November 11, 2005; 280(45): 37471 - 37480.
[Abstract] [Full Text] [PDF]


Home page
Drug Metab. Dispos.Home page
M. F. Paine, S. S. Ludington, M.-L. Chen, P. W. Stewart, S.-M. Huang, and P. B. Watkins
DO MEN AND WOMEN DIFFER IN PROXIMAL SMALL INTESTINAL CYP3A OR P-GLYCOPROTEIN EXPRESSION?
Drug Metab. Dispos., March 1, 2005; 33(3): 426 - 433.
[Abstract] [Full Text] [PDF]


Home page
Mol. Endocrinol.Home page
E. Rico-Bautista, C. J. Greenhalgh, P. Tollet-Egnell, D. J. Hilton, W. S. Alexander, G. Norstedt, and A. Flores-Morales
Suppressor of Cytokine Signaling-2 Deficiency Induces Molecular and Metabolic Changes that Partially Overlap with Growth Hormone-Dependent Effects
Mol. Endocrinol., March 1, 2005; 19(3): 781 - 793.
[Abstract] [Full Text] [PDF]


Home page
J. Biol. Chem.Home page
C. A. Wiwi and D. J. Waxman
Role of Hepatocyte Nuclear Factors in Transcriptional Regulation of Male-specific CYP2A2
J. Biol. Chem., February 4, 2005; 280(5): 3259 - 3268.
[Abstract] [Full Text] [PDF]


Home page
J EndocrinolHome page
S Eleswarapu and H Jiang
Growth hormone regulates the expression of hepatocyte nuclear factor-3 gamma and other liver-enriched transcription factors in the bovine liver
J. Endocrinol., January 1, 2005; 184(1): 95 - 105.
[Abstract] [Full Text] [PDF]


Home page
EndocrinologyHome page
P. Kotokorpi, C. Gardmo, C. S. Nystrom, and A. Mode
Activation of the Glucocorticoid Receptor or Liver X Receptors Interferes with Growth Hormone-Induced akr1b7 Gene Expression in Rat Hepatocytes
Endocrinology, December 1, 2004; 145(12): 5704 - 5713.
[Abstract] [Full Text] [PDF]


Home page
Mol. Cell. ProteomicsHome page
E. V. Laz, C. A. Wiwi, and D. J. Waxman
Sexual Dimorphism of Rat Liver Nuclear Proteins: Regulatory Role Of Growth Hormone
Mol. Cell. Proteomics, December 1, 2004; 3(12): 1170 - 1180.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Gastrointest. Liver Physiol.Home page
F. R. Simon, J. Fortune, M. Iwahashi, I. Qadri, and E. Sutherland
Multihormonal regulation of hepatic sinusoidal Ntcp gene expression
Am J Physiol Gastrointest Liver Physiol, October 1, 2004; 287(4): G782 - G794.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental Data
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow Request Copyright Permission
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Ahluwalia, A.
Right arrow Articles by Waxman, D. J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Ahluwalia, A.
Right arrow Articles by Waxman, D. J.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Endocrinology Endocrine Reviews J. Clin. End. & Metab.
Molecular Endocrinology Recent Prog. Horm. Res. All Endocrine Journals